Diseño Estructural 2012

UNIVERSIDAD AUTONOMA DE QUERETARO
APUNTES DE DISEO ESTRUCTURAL
UNIVERSIDAD AUTNMA DE QUERTARO
Facultad de Ingeniera Ingeniera Civil
APUNTES DE
DISEO ESTRUCTURAL
SEXTO SEMESTRE
DR. GUADALUPE MOISS ARROYO CONTRERAS
Divisin de Estudios de Posgrado Centro Universitario, Cerro de las Campanas, C.P. 76010, Stgo. de Quertaro, Qro. Tel.: (42) 192-12-00, EXT. 6071, Email: marroyoc@uaq.mx
Facultad de Ingeniera
Arroyo M.
DISEO ESTRUCTURAL
(COMPORTAMIENTO, ANLISIS Y DISEO DE ESTRUCTURAS SOMETIDAS A EXCITACIONES DINMICAS, TIPO SISMO, VIENTO, OLEAJE, TRFICO)
GUADALUPE MOISS ARROYO CONTRERAS Profesor de la Facultad de Ingeniera, UAQ Tel.: 2-16-35-99 ext. 113 y 143, Fax.: 2-16-35-99 ext. 103 E-mail: marroyoc@sunserver.uaq.mx El objetivo de este curso es aprender los mtodos y algoritmos en que est basada la filosofa del anlisis y diseo de estructuras sometidas a excitaciones dinmicas, como el sismo, viento y oleaje. Adems, aprender una metodologa para el clculo de la probabilidad de falla o de seguridad implcita en este tipo de sistemas estructurales.
CONTENIDO TEMTICO
I. DINMICA ESTRUCTURAL I.1 Sistemas lineales de un grado de libertad I.2 Sistemas lineales discretos de varios grados de libertad I.3 Nociones sobres sistemas continuos I.4 Nociones sobre sistemas no-lineales I.5 Vibraciones aleatorias (procesos gaussianos estacionarios) I.6 Respuesta aleatoria de sistemas lineales de un grado de libertad I.7 Respuesta aleatoria de sistemas lineales de varios grados de libertad II. EXCITACIONES DINMICAS (sismo, viento y oleaje) II.1 Sismo II.1.1 Caractersticas de los temblores II.1.2 Sistemas estructurales empleados en edificios II.1.3 Anlisis y diseo ssmico de edificios (esttico y dinmico) II.2 Viento II.2.1 Caractersticas del viento II.2.2 Acciones de viento en estructuras II.2.3 Accin dinmica del viento en estructuras II.3 Oleaje II.3.1 Teora de la ola lineal II.3.2 Tratamiento estocstico de la ola II.3.3 Modelo de fuerza de oleaje II.3.4 Anlisis estocstico de la respuesta de un sistema mecnico, tipo estructura marina, sujeto a oleaje estacionario II.3.5 Aproximaciones casi-esttica y por el mtodo modal III. SOLUCIONES NUMRICAS DE LA ECUACIN DIFERENCIAL DE EQUILIBRIO DINMICO.
Arroyo M.
IV. NORMATIVIDAD
Arroyo M.
BIBLIOGRAFA 1. Cloug R.W. and Penzien J., Dynamics of Structures, Mc. Graw-Hill, 1975. 2. Thomson W.T., Theory of vibration with applications, fourth edition, Prentice Hall, Englewood cliffs, New Jersey, 1993. 3. Bathe K.J. and Wilson E.L., Numerical Methods in Finite Element Analysis, Prentice-Hall, 1976. 4. Newmark y Rosenblueth E., Fundamentals of earthquake engineering Prentice Hall Inc., traduccin editorial Diana, Mxico, 1976. 5. Bazan E. y Meli R., Manual de diseo ssmico de edificios, Limusa-Noriega, 1990. 6. Sarpkaya T. and Isaacson M., Mechanics of Wave forces on offshore structures, Van Nostrand Reinhold Company, USA, 1981. 7. THOPT-CHRISTENSEN, P. & MUROTSU, Y. "Application of Structural Systems Reliability Theory", Springer-Verlag Berlin, Heidelberg, 1986. 8. THOFT-CHRISTENSEN, P. & BAKER MICHAEL J., "Structural Reliability Theory and its Applications", Springer-Verlag Berlin, Heidelberg, 1982. 9. ANG Alfred H-S. & TANG Wilson H., "Probability Concepts in Engineering Planningand Design", Volume II: Decision, Risk, and Reliability, Chapters 6 y 7, John Wiley & Sons, Inc., 1984. 10. BJERAGER, P., Reliability Analysis of Structural Systems", Phi. D Thesis, The Technical University of Denmark. 11. AUGUSTI G., BARATTA A. and CASCIATI F., Probabilistic Methods in Structural Engineering", Chapman and Hall 1984.
Arroyo M.
CRITERIOS DE DISEO CARGAS MUERTAS CARGAS VIVAS ACCIONES (FUERZAS) Clasificacin en base a la duracin con que obran sobre la estructura: - Acciones permanentes (gravitacionales) * Carga muerta * Empuje esttico (presiones hidrostticas) * Deformaciones y desplazamientos impuestos (presfuerzo, movimientos diferenciales) - Acciones variables * Carga viva * Cambios de temperatura (fuerzas trmicas) * Carga de transito en puentes * Cargas de impacto (LRFD-A4.2) * Cargas longitudinales (Frenaje, Choque de un barco contra un muelle) * Fuerzas centrfugas (puentes curvos) - Acciones accidentales * Sismo * Viento * Oleaje * Explosiones * Cargas de piso: hielo, nieve, lluvia (encharcamiento)
Arroyo M.
UNIVERSID DAD AUTON NOMA DE QU UERETARO
AP PUNTES DE DISEO D ESTRUCTURAL
CAR RGAS MU UERTAS PE ESO PRO OPIO Comenta arios al Ttulo T Sex xto, Seguridad Estructural de las Construcc C ciones, del d Reglame ento de C Construc ccin para a el Distr rito Fede eral, Serie es del II-UNAM, No ES-1, 19 991.
Facultad d de Ingeniera
Arr royo M.
Pesos Volumtric V cos de Ma ateriales para Con nstruccin n, Valores de dise eo segn n el RDF-76. .
Arr royo M.
CO ONCEPTO OS DE TEORA T D PROB DE BABILIDA ADES Variable aleatoria a (v.a.), es e aquella a cuyo va alor no pu uede ser fijado con n precisi n omento de tomar una u desic cin. en el mo Las probabilidad des de que una v.a. v cont tinua tom me difere entes valo ores, est tn descritas s por la fu uncin de e distribucin de p probabilid dades (o funcin de d densid dad de proba abilidades s)
f X (x ) .
Func cin de de ensidad d de probab bilidades Si la v. a. es continua, la a funcin n es tal que q la pr robabilida ad de la variable se encuentr re en un intervalo i diferencial alrededor de cie erto valor r
x0 , por r tanto, va ale:
d dx d dx P x0 < X < x0 + = f x ( x0 )dx 2 2

anto, la probabilida ad de que la varia able alea atroira x to ome un valor v men nor Por lo ta que ciert ta cantida ad X m , re esulta:
P[ X < X m ] =
xm
f x ( x)dx = FX ( X m )
donde FX ( X m ) se s le llam ma funcin n de distr ribucin acumulad a da de probabilidad des de la var riable ale eatoria. Facultad d de Ingeniera Arr royo M.
Por otro lado, la probabilidad de que X resulte mayor que cierto valor
X m , es:
P[X > X m ] =
xm
f X ( x)dx = 1
xm
f X ( x)dx = 1 FX ( X m )
CARACTERISTICAS MS IMPORTANTES DE LA FUNCIN DE DENSIDAD DE PROBABILIDADES - La posicin del centroide del rea bajo la curva, primer momento o la media o la esperanza matemtica de la variable aleatroria (v.a.) es:
E [x ] = m X =
x f X ( x ) dx
- El momento de inercia del rea bajo la curva, o la variancia, o segundo momento:

2 Var ( X ) = x = 2 ( ) x m x
f X ( x ) dx
Se llama desviacin estandar a la raiz cuadrada de la variancia:
X = Var ( X )
Se llama coeficiente de variacin a la relacin:
CVX =
X
mX
La Var ( X ) , X y el CVX miden de distinta manera el grado de incertidumbre que se tiene a cerca de la variable aleatoria (v.a.).
Arroyo M.
Ley de distribucin Normal La variable aleatoria x sigue un ley N (mX , X ) , si la funcin de densidad de probabilidades se escribe como:
f X ( x) =
x 2
1 x mx 2 x
Funcin de distribucin acumulada
FX ( x) =
f x ( x) dx =
e 2
1 x mx 2 x
dx
Ley de distribucin Normal estandar si x sigue una ley normal, y la v.a. u es igual a:
u=
x mX
x = u X + mX
x = m X + u X xmx = mX + 2 X
xmn = mX 2 X
entonces la v. a. u sigue un ley normal estandar, N(0,1), es decir, una distribucin normal con media cero y desviacin estandar 1.
u (u ) =
1 2u2 e 2
Arroyo M.
U (u ) =
u (u ) du u=
1 2
u2 2
dx
Si requie ere obten ner la func cin de distribuci d n de pro obabilidad des acum mulada pa ara valores de d x, Fx ( x0 ) , se proce p com mo sigue: :
x mx x0 mx FX ( x0 ) = P[ X < X 0 ] = P < x x x0 mx FX ( x0 ) = P[u < u0 ] = = (u0 ) x
Si requie ere obten ner la probabilidad d acumula ada para valores n negativos s de u se proce co omo sigue e: Facultad d de Ingeniera Arr royo M.
U ( u ) = 1 U (u )
FRMULA PARA A OBTEN NER LOS S VALORES DE LA A FUNCI IN DE BUCIN NORMAL L ESTAN NDAR AC CUMULAD DA, DISTRIB
U (u ) .
Arr royo M.
APLICACIN DE LA PROBABILIDAD EN EL DISEO PARA TOMAR EN CUENTA LA INCERTIDUMBRE EN EL CONOCIMIENTO DE LAS ACCIONES, LOS REGLAMENTOS ESPECIFICAN VALORES DE DISEO QUE SON CONSERVADORES (O SEA, VALORES QUE TIENEN UNA PROBABILIDAD PEQUEA DE SER EXCEDIDOS DEL LADO DESFAVORABLE, Valores nominales). SEA LA ACCIN DEBIDA AL PESO PROPIO DE UNA ESTRUCTURA DE CONCRETO. ESTA DEPENDE DE LAS DIMENSIONES DE LOS ELEMENTOS Y DEL PESO VOLUMTRICO DEL CONCRETO. EL PESO VOLUMTRICO DEL CONCRETO SIMPLE ES DE ACUERDO A LOS MANUALES, 2,200 Kg / m3
EJEMPLO, CONSTRUCCIN DE LA CORTINA DE UNA PRESA. DE ACUERDO A LOS DATOS PROPORCIONADOS POR VARIAS PLANTAS DEL VALLE DE MXICO, PARA LA ELABORACIN DE CONCRETOS ENTRE 200 Y 300 Kg / cm 2 de resistencia, Y PARA DIFERENTES BANCOS DE GRAVA Y ARENA, SE TIENE PARA LA V.A. X DEL PESO VOLUMTRICO, QUE:
mx = 2137 kg / m3
x = 39.4 kg / m3
CVx = 0.01846
SE BUSCA OBTENER EL VALOR NOMINAL DE LA V.A. X, CON PORCENTIL 98. O SEA: EL VALOR DE LA V.A. QUE TIENE 98% DE PROBABILIDADES DE NO SER EXCEDIDA O EL 2% DE PROBABILIDADES DE SER REBASADO.
X m f ( x ) dx = 0.98 P [ X X m ] = X
Facultad de Ingeniera Arroyo M.
USANDO TABLAS SE PROCEDE COMO SIGUE:
U=
X m X
X 98 = X mx =m X +U X X 98 = X mx = m X +2 X
EN LAS TABLAS SE TIENE PARA: UN REA DE 0.97725 UN VALOR DE LA V.A. U= 2, Y UN REA DE 0.98214 UN VALOR DE LA V.A. U= 2.1 INTERPOLANDO LINEALMENTE SE TIENE: PARA UN REA DE 0.98 UN VALOR DE LA V.A. U = 2.0562 Con la expresin aproximada se tiene: U = 2.054
POR LO TANTO, EL VALOR CARACTERSTICO O NOMINAL BUSCADO, O VALOR MXIMO PROBABLE DE LA V.A. X, ES:
X 98 =m X +U X = 2,137+2.054(39.4)=2,218kg / m3
ESTE VALOR DE 2,218 Kg / m3 DEBE USARSE EN EL CALCULO DEL PESO PROPIO, PARA LA REVISIN DE ESFUERZOS EN EL CONCRETO Y EN EL SUELO, POR EFECTO DE LAS CARGAS GRAVITACIONALES.
X m f ( x ) dx = 0.02 P [ X X m ] = X
DE LA MISMA MANERA, USANDO TABLAS SE PROCEDE COMO:
Arroyo M.
U=
X mX
X 98 = mX + U X
POR SIMETRA DE LA FUNCIN DE DENSIDAD DE PROBABILIDADES DE LA V.A. U, SE TIENE QUE: PARA UN REA DE 0.02 UN VALOR DE LA V.A. U = -2.054
POR LO TANTO, EL VALOR CARACTERSTICO O NOMINAL BUSCADO, O VALOR MNIMO PROBABLE DE LA V.A. X, ES:
X 02 =mX +U X = 2,1372.054(39.4)=2,056kg / m3
ESTE VALOR DE 2,056 Kg / m3 DEBE USARSE PARA LA REVISIN DE LA POSIBILIDAD DE VOLTEO DE LA CORTINA ANTE EL EMPUJE HIDROSTTICO. EN ESTE CASO ES CONVENIENTE CALCULAR EL PESO PROPIO DE LA CORTINA CON EL VALOR MNIMO QUE PUDIERA ALCANZAR.
Arroyo M.
LEY DE DISTRIBUCIN LOG-NORMAL La v.a. Y sigue una ley Log-normal, si X sigue una ley normal
U=
log Y0 mX
Arroyo M.
Estadisticamente, se tiene que: Media de una muesta:
x=
x
n
Variancia de una muestra:

2 ( ) x x i
S2 =
Arroyo M.
FUNCIN DE DENSIDAD DE PROBABILIDADES CONJUNTA DE DOS VARIABLES ALEATORIAS X e Y

f X ,Y ( x , y ) = 1 2 x y ( x m x )2 2 ( x m x )( y m y ) ( y m y )2 1 exp + 2 2 2 2 x y y 2(1 ) 1 x
DE n VARIABLES X
f X ( X ) = f X 1 , X 2 ,..., X n ( x1 , x 2 ,..., x n ) =
donde:
[C x ]1 (2 )n / 2
1/ 2
T 1 [x m x ] [C x ] [x m x ] exp 2
x1 mx1 x 2 m x2 . ; [x mx ] = . . xn mxn
[C x ] = Ci j = E [(xi mx )(x j mx
i
)] = C
2 x , i= j
i
xi x j
,i j
C xi x j = xi x j xi x j
Arroyo M.
FUNCIN DE VARIAS VARIABLES ALEATORIAS NORMALES Para variables aleatorias independientes (v. a.i.) no correlacionadas de tiene:
Si g (x ) = ai X i
i =1
E (g (x )) = ai E [ X i ] = ai X i
i =1 i =1
Var ( g ( x )) = ai2 Var [X i ]

i =1
Si g ( x ) = X i
i =1
E ( g ( x )) = E [X i ] = X i
i =1 i =1
n 2 2 Var ( g ( x )) m x j xi i =1 j =1 j i
n
para n = 3 se tiene:
2 2 2 2 2 2 x21 + mx x22 + mx x23 Var ( g ( x )) mx mx mx mx 2 3 1 3 1 2 2 2 2 2 2 2 2 2 2 2 2 2 m m CV m m m CV m m m CV = mx + + x x x x x x x x x x x 2 3 1 1 1 3 2 2 1 2 3 3 2 2 2 2 2 2 m m CV CV CV = mx + + x x x x x 2 3 1 1 2 3
CV ( g (x )) =
Var ( g ( x )) CVx2 + CVx2 + CVx2 1 2 3 E (g ( x ))

1/ 2
n CV ( g (x )) = CVx2 i i =1
Arroyo M.
CARGA A VIVA EN N EDIFIC CIOS
Donde
es la carga viv va, CV, m media, wa CV insta antnea y wm CV V mxima a. Arr royo M.
AP PUNTES DE D DISEO ESTRUCTURAL
Arr royo M.
MODELO PROB BABILIST TA PARA A LA DE ETERMIN NACIN DE LAS S CARGA AS VIVAS, ESPECIF E FICADAS S POR EL L RC-DF F-76.
d la car rga viva en el re ea de infl luencia d de una columna. Distribucin de VO. OBT TENER UNA U CA ARGA UN NIFORME EMENTE E DISTRI IBUIDA DE D OBJETIV DISEO O. VARIAB BLES: a es e un re ea eleme ental,
ai
reas ele ementale es, s fue erza inter rna provo ocada por r la carga a total, en una secc cin crtic ca. El efec cto de la carga a total,
w
w*
un conju unto de ca argas act tuando so obre rea as elementales, la a carga
carga un niformemente distr ribuida de e diseo, , tal que la fuerza interna que q tenga un p produce, na probab bilidad pr re-estable ecida de ser exced dida.
a) OB BTENER UNA FU UERZA IN NTERNA NOMINA AL, MXI IMA PRO OBABLE, O CA ARACTRISTIA S .
*
- MODELO O ANAL LTICO O Y PR ROBAB BILIST TA

La fuerza interna que prod duce cad da carga e es:
Arr royo M.
S i = a z i wi
La fuerza total es:
S = a z i wi
i =1
donde
wi
es la intensidad de carga por unidad de rea, se considera una variable aleatoria (v.a.), a un rea unitaria, se considera una variable determinista (v. d.),
un factor que depende de la posicin del rea considerada, del tipo de estructura y de la fuerza interna en estudio, s se considera una v. a. Se supone que la carga otras cargas, es decir, no correlacionadas. Dado que las v. a. de s son:
zi
wi , sobre cada rea elemental, es independiente de las
wi son v. a. normales e Independientes (v. a.n.i.), es decir
wi son normles independientes (no correlacionadas) y s una funcin lineal de estas v. a. wi , se tiene que la media (o esperanza) y la variancia
S = a z i wi
i =1
Si g (s ) = a zi wi
i =1
E (g (s )) = S = ms = a zi E [wi ] = a zi wi
i =1 n i =1
S = a zi mwi = a mw zi
i =1 i =1
Arroyo M.

n
Var (g (s )) = = a
2 s
Var[z w ] = a z
2 i =1 i i i =1 2 i n i =1
2 wi
2 = a 2 w zi2
dado que a es una superficie unitaria a = a =1 .

2
Var ( g (s )) = a
2
2 w
z
i =1
2 i
= a
2 w
z
i =1
2 i
= s2
2 n 2 s = a w z i i =1
1/ 2
De acuerdo a la expresin de exedencia obtenida anteriormente, el valor caracterstico de la fuerza interna
s* es:
u=
s ms
s * = ms + u s = ms + k s
2 n 2 s = amw zi + k a w zi i =1 i =1
n * 1/ 2
La carga uniformemente distribuida nominal
w* es:
Arroyo M.

n * = n w = n a mw z i + i =1 a zi a zi
s*
2 n 2 k a w z i i =1
1/ 2
i =1
i =1
= mw +
n 2 zi k i =1 n a zi
i =1
1/ 2
Si n es el nmero de reas elementales que componen el rea de influencia AT o aquella rea en que z tiene un valor significativo, se tiene que el trmino constante de la desviacin estandar de w se puede expresar como:
n 2 zi n k i=1 n n n a zi
i =1
1/ 2
n 2 zi i =1 n k n na zi
i =1
1/ 2
k (mz 2 ) AT mz
1/ 2
z
donde:
i =1
= mz ,
z
i =1
2 i
= mz 2 y na = AT .
k mz 2 AT mz
Por tanto
w* = mw +
[ ]
1/ 2
w.
haciendo
[m ] C=
z
2
1/ 2
mz
se tiene que:
Arroyo M.
w* = mw +
1 / 2 w* = mw + (kC w )AT
donde k es un factor de excedencia o factor de la probabilidad de excedencia
k C w AT
C vara entre 1.1 y 1.2, se puede considerar constante unitaria w. La expresin para la CV (carga viva) tiene la forma:
mw y w son parmetros de la funcin de distribucin de la carga
w* = a + b AT1/ 2
Los valores de a y b se pueden deducir estadsticamente y probabilsticamente. Estos valores estn especificado por los reglamentos. El reglamento de construcciones para el DF propone para la CV de diseo las expresiones siguientes
w = 180 + 420 A
*
1 2 1 2
= 100 + 420 A
Arroyo M.
Los valo ores de a y b se pueden de educir es stadstica y probab bilsticam mente. Est tos valores estn e esp pecificado os por los s reglame entos.
Arr royo M.
NOCIONES DE CONFIABILIDAD
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
I.
CONCEPTOS BSICOS DE DINMICA
VIBRACIN puede ser definida como un movimiento oscilatorio DINMICA implica el estudio del cuerpo en movimiento y las fuerzas que causan el movimiento La dinmica se divide en dos partes: a) CINEMTICA es el estudio de la geometra de movimiento sin referencia a la fuerza o masa b) CINTICA es el estudio de la relacin entre el movimiento de los cuerpos y las fuerzas que actan sobre ellos. RESEA HISTRICA GALILEO (1564 - 1642), hizo la primera contribucin a la dinmica SIR ISAAC NEWTON (1642 - 1727). Con las experiencias de Galileo, sobre cuerpos uniformemente acelerados, condujeron a Newton a enunciar las leyes fundamentales del movimiento. HUYGENS (1629 - 1695). En 1584 Galileo concibi el principio del pndulo y en 1657 Huygens aplic el principio al reloj. BARN JOHN WILLIAM STRUTT, LORD RAYLEIG (1842 1919). A causa de las mquinas de alta velocidad, es hasta fines del siglo XIX, que Rayleigh organiz y desarroll la teora de vibracin mecnica. La teora actual est basada en sus contribuciones.
Arroyo M.
I.1 ECUACIONES DIFERENCIALES DE MOVIMIENTO DE LA CINEMTICA DE UNA PARTCULA

DESPLAZAMIENTO LINEAL es la distancia que un punto ha recorrido sobre una trayectoria, medida a partir de un origen seleccionado convenientemente (vector). - S el origen es fijo el desplazamiento es absoluto, - S el origen est en movimiento desplazamiento es relativo. el
VELOCIDAD INSTANTNEA DE UN PUNTO es la razn del tiempo de cambio del desplazamiento s :
ds v = = dt s
(1)
vx
dx = = dt x
ACELERACIN INSTANTNEA es el cambio de la razn de tiempo de la velocidad:
d v a = = = dt v s
d vx = vx = x ax = dt
(2)
Despejando dt de la ec. 1 y remplazndola en la ec. 2 se tiene:
a =v
dv a d s = vd v ds
(3)
ax = vx
d vx a d x = vx d vx dx
Arroyo M.
I.2 ECUACIONES DIFERENCIALES ROTACIN DE UN CUERPO RGIDO
DE
MOVIMIENTO
DE
DESPLAZAMIENTO ANGULAR es el cambio en posicin angular de una lnea respecto al eje z (dada como medida de lnea de referencia del movimiento de un cuerpo rgido girando en el plano x-y), seleccionada convenientemente. El ngulo entre las lneas AB Y AB es un desplazamiento angular de la lnea AB: . La longitud del vector desplazamiento angular. VELOCIDAD ANGULAR
AB
AB
es proporcional a la magnitud del
d z = = dt z
(4) ACELERACIN ANGULAR
d = = dt z
(5) Despejando dt de la ecuacin 4 y remplazndola en la ec. 5, se tiene: Facultad de Ingeniera Arroyo M.
d d
d = d
(6)
Arroyo M.
I.3 MOV VIMIENT TO CUR RVILNE EO
V = l
cin del v vector La direcc
V es tang gente a la a trayecto oria del p ndulo.

V es pro oporciona al a la tasa de tie empo de cambio del
La magn nitud del vector
desplaza amiento angular, a r El vector
l es la a distancia a de A a P.
Arr royo M.
I.4 COO ORDENA ADAS (GENER ( RALIZAD DAS)

ccin de las l coord denadas es e un aspecto mu uy import tante en el e estudio oy La selec solucin de probl lemas en n dinmica. mplo: las s coorden nadas car rtesianas s Se tiene por ejem Las s coordenadas po olares Las s coordenadas cilndricas Las s coordenadas es sfricas DENADAS S GENE ERALIZA ADAS so on un ju uego de coordenadas que q COORD describe e el movim miento ge eneral e identifica la sujeci n.
El pn ndulo tien ne dos g grados de e libertad d S el pn ndulo es desplaza ado por un u ngulo o permane ecer en dicho pla ano z y .
en el plano z y , el movimien nto
ndulo es s desplaz zado por un ngu ulo en el plano S el p z x permane ecer en dicho pla ano .
z x , el movimien nto
dente que el mo ovimiento o en la direcci n de una coordenada es Es evid independ diente de e la otra. MODO PRINCIP P PAL DE MOVIMIE M ENTO es s el movim miento en el cual l solamen nte una coor rdenada vara. Facultad d de Ingeniera Arr royo M.
COORDENDA PRINCIPAL. Se llama coordenada principal, si solamente esta coordenada es necesaria para definir el movimiento. En sistemas lineales, todo movimiento puede ser descrito por la superposicin de modos principales.
Arroyo M.
II.
MOVIMIENTO OSCILATORIO
Vibraciones. Existen dos clases de vibracin: La vibracin libre y la vibracin forzada. Vibracin libre. La vibracin libre es cuando la oscilacin del sistema se debe nicamente a la accin de fuerzas inherentes del sistema. Un sistema bajo vibracin libre, vibrar en una o ms de sus frecuencias naturales, las cuales son propiedades del sistema dinmico establecido por su masa y distribucin de rigidez. Vibracin forzada. La vibracin forzada se da bajo la excitacin de fuerzas externas. Cuando la excitacin es oscilatoria, el sistema es forzado a vibrar con la frecuencia de excitacin. De acuerdo a lo anterior, se puede decir que el clculo de las frecuencias naturales es de suma importancia en el estudio de vibraciones. Los sistemas vibrando estn sujetos al amortiguamiento en mayor o menor grado, ya que la energa es disipada por friccin y algunas otras resistencias. El nmero de coordenadas independientes requerido para describir el movimiento de un sistema es llamado GRADOS DE LIBERTAD DEL SISTEMA. II.1 MOVIMIENTO ARMNICO Movimiento peridico. El movimiento peridico se tiene cuando el movimiento se repite en iguales intervalos de tiempo t. Si el movimiento est descrito por la funcin tiempo x(t), entonces para todo movimiento peridico se tiene:
x(t ) = x(t + T )
donde T es el periodo de oscilacin. La frecuencia es:
f = 1/ T
El movimiento armnico es la forma ms simple del movimiento peridico.
Arroyo M.
Sea el movimiento armnico generado por el sistema masa-resorte, mostrado en la figura 1, siguiente.
Figura 1. Movimiento armnico. El desplazamiento vertical X de la masa se puede obtener como sigue:
2 x(t ) = A sen t T
donde A es la amplitud de la oscilacin. Movimiento armnico circular.
(1)
El movimiento armnico obtenido como una proyeccin de un punto P movindose circularmente, se muestra en la figura siguiente:
Figura 2. Movimiento armnico producido por un movimiento angular. El desplazamiento x puede ser expresado como:
Arroyo M.
x(t ) = A sen(t )
(2)
donde A es la amplitud de oscilacin y la velocidad angular de la lnea OP. La velocidad angular , se mide generalmente en radianes por segundo, y est referida a la frecuencia circular, ya que el movimiento se repite cada 2 radianes:
2 = 2 f T
(3)
La velocidad del movimiento armnico se obtiene derivando x respecto al tiempo, ec 2, es decir:
dx(t ) d ( A sen(t )) = dt dt . x = Acost = Asen(t + / 2)

La aceleracin del movimiento armnico es:
(4) (5)
.. x = 2 Asent = 2 Asen(t + )
(6)
La aceleracin del movimiento, expresada en la ecuacin anterior 6, se puede representar en funcin del desplazamiento x (ec. 2):
.. x = 2 x
(7)
La ecuacin diferencial homognea lineal de segundo orden del movimiento armnico resulta ser:
.. 2 x + n x=0
(8)
La solucin a esta ecuacin diferencial se obtiene como a continuacin se describe.
Arroyo M.
Solucin general:
x(t ) = A sen(nt ) + B cos(nt )

si t = 0, x(0) = B
(9) (10)
Derivando x respeto al tiempo, en la ecuacin de la solucin general (ec. 9) se tiene:
. x = n A cosnt n Bsen(nt )
(11)
. si t = 0, x (0) = n A . x (0)
(12)
A=
Substituyendo los valores de A y B en la ecuacin de la solucin general, se tiene:
x(t ) =
. x (0)
sen(nt ) + x(0) cos(nt )
(13)
Arroyo M.
II.2 ECUACIN DE MOVIMIENTO DE UN SISTEMA MASA - RESORTE Sea el sistema Masa Resorte, de un grado de libertad, mostrado en la figura 3 siguiente.
Figura 3. Sistema Masa Resorte de un grado de libertad. Por equilibrio esttico de las fuerzas que estn actuando en la masa del diagrama de cuerpo libre mostrado en la figura 3a, se tiene:
= 0,
W k = 0,
W = k
(14)
donde es el desplazamiento o alargamiento del resorte en la direccin x, provocado por la carga W. De acuerdo a la segunda ley de Newton se tiene que:
W = k = mg
Por equilibrio dinmico de la masa de la figura 3b se tiene que:
(15)
mx =
..
..
= W k ( + x ) = (W k ) kx
= W W kx
m x = kx
Facultad de Ingeniera (16) Arroyo M.
La aceleracin del sistema masa resorte es:
x=
..
k x m
(17)
FRECUENCIA NATURAL
fn .
La frecuencia natural es una propiedad del sistema, definida como:
fn =
1 Tn
(18)
Expresada en funcin de la frecuencia angular se tiene:
n =
fn =
2 T
o T =
n
(19)
1 Tn
fn =
n 2
Por otro lado, la aceleracin del movimiento armnico y la aceleracin del sistema masa resorte se expresan respectivamente como:
.. 2 x = n x
x=
..
k x m
(20)
Por lo tanto, se tiene que la frecuencia angular natural del sistema es:
n =
k m
(21)
La frecuencia natural de vibrar y el periodo natural de vibrar son respectivamente: Facultad de Ingeniera Arroyo M.
fn = Tn =
n 1 = 2 2
1 2 = fn n
m k
k m
(22)
Tn = 2
(23)
Dado que k = mg (ec. 15), la frecuencia natural, expresada en la ec. 22, se puede representar como:
fn = fn =
1 2
g
1 k = m 2
mg m
1 2
(24)
Arroyo M.
EJEMPLO 1. Obtener la frecuencia natural del sistema estructural tipo viga en voladizo, mostrado en la figura 4.
Figura 4. Viga en voladizo.
= f P=
1 P k
l3 1 l3 3EI = 3EI P = 3EI (mg ) = k mg , donde k = l 3 y P = mg = W

La frecuencia natural de vibrar de la viga es por tanto:
fn =
1 2
1 2
g 1 = 2
3EI ml 3
3EI g 1 = l 3 mg 2
3EI ml 3
fn =
Viga simplemente apoyada:
k1 Pl Pl 3 = + 4GA 48EI z
Arroyo M.
II.3 VIBRACIONES LIBRES CON AMORTIGUAMIENTO VISCOSO Ecuacin de equilibrio amortiguamiento: dinmico de un sistema masa resorte con
.. m x + Fd + kx = F (t )
donde F (t ) es la funcin excitadora y
(25)
Fd la fuerza de amortiguamiento.
Si la fuerza de amortiguamiento viscoso se expresa como:
. Fd = c x
(26)
donde c es una constante de proporcionalidad. La ecuacin de equilibrio dinmico se puede entonces representar como:
.. . m x + c x + kx = F (t )
.. . m x + c x + kx = 0
Si se supone una solucin de la forma
st
(27)
Si F (t ) = 0 , corresponde fsicamente a una vibracin libre amortiguada: (28)
x = e st
se tiene: (29)
. .. st x = e , x = se y x = s 2 e st
Remplazando estos valores en la ecuacin de vibracin libre amortiguada, ec. 28, se tiene:
(ms
.. . m x + c x + kx = ms 2 e st + cse st + ke st = 0
2
+ cs + k e st = 0
(30)
Arroyo M.
Dado que la solucin de la ecuacin diferencial, tiene entonces que:
e st , no puede valer cero, se
ms 2 + cs + c s2 + s+ m c s2 + s+ m
k=0 k =0 m
(31)
2 = 0
La ec. 3, se le conoce como la ecuacin caracterstica, donde los valores de s se pueden obtener como sigue. - la solucin a una expresin de segundo orden, ec 31, est dada por:
A1 x 2 + B2 x + C 3 = 0 B B 2 4 AC B = 2A 2A B2
x =
(2 A)2
4 AC 4 A2
2 B C B 2A A 2A
- por lo tanto para s se tiene:
s2 +
c k s+ =0 m m
donde A1 = 1, B2 =
c k y C 3 = , remplazando los valores en la ecuacin m m anterior se obtienen los valores de s:
Arroyo M.
c c k s1, 2 = 2m 2m m
(32)
La solucin general de la ecuacin de vibracin libre amortiguada est dada por:
x = Ae s1t + Be s 2 t
c 2m + k c m 2m
2
x = Ae
+ Be
c 2m
k c m 2m
x = Ae
x=e
2 c k c t t 2 m m 2m
+ Be
2 c k c t t 2 m m 2m
c t 2m
Ae
k c m 2m
+ Be
2 c k t m 2m
(33)
Si
c 2m
es
2
menor
que
k m
el
exponente
se
vuelve
imaginario
k c i m 2m y como:
Arroyo M.

2
i c t 2m x=e Ae
i e
k c m 2m
+ Be
k c m 2m
c 2 2m m k
t 2 2 k c c = cos k t i sen t m m m m 2 2
(34)
La solucin obtenida contiene trminos oscilatorios.
Arroyo M.
Se hace notar que si:
c k > m se trata de un sistema sobre amortiguado, 2m c k < m se tiene un sistema sub-amortiguado, movimiento oscilatorio, 2m c k = m se tiene un amortiguamiento crtico, anulndose la raz o el 2m radical.
2 2
Amortiguamiento crtico

2 cc k = = 2 n 2m m
cc = n 2m
cc 2m
(35)
cc =
4km = 2 km = 2m n n =
Si se expresa un amortiguamiento cualquiera en funcin del amortiguamiento crtico se tiene: adimensional.
= c / cc , donde
se conoce como razn de amortiguamiento
Si se expresan ahora las races en trminos de se tiene:
c c = c = n 2m 2m
Por consiguiente, las races expresadas en trminos de , son:
(36)
Arroyo M.
c k c s1, 2 = 2m m 2m = n
( n )2 n2
s1, 2 = 2 1 n
La condicin depende ahora de , por consiguiente si:
(37)
> 1
0 < < 1
caso sobre-amortiguado caso sub-amortiguado, movimiento oscilatorio caso crtico, se anula la raz o el radical.
= 1
S1, 2 = S1, 2
( )2
1 n
= i 1 (
)2
=1 =0
s1, 2
n
s1, 2
=1 =i
Si Figura 5 Representacin grfica
Arroyo M.
II.4 MOVIMIENTO OSCILATORIO (movimiento libre amortiguado). Caso sub-amortiguado
< 1.
La Solucin general est dada por:
x = Ae s1t + Be s2t
s1, 2 = n
( n )2 n2
(38a)
1
2
s1, 2 = n i n 1 2
x=e
n t
(i Ae
) t
n
+ Be ( i
) t
n
(38b)
x = e nt C1sen 1 2 n t + C 2 cos 1 2 n t
(39)
x = X 1e nt sen 1 2 n t + 1
De acuerdo al teorema de Euler - Moivre de los nmeros complejos y del clculo diferencial e integral, se tiene:
eiy + e iy eint + e int cos y = = 2 2 eiy e iy e int e int sin y = = 2 2

e iy = cos y i sin y
O
e iy = cos y i sin y
Arroyo M.
e i n t + e i n t e i n t e i n t u (t ) = A +B 2 2
e iy + e iy cos y = 2 cos 1 2 n t x = e nt C1sen 1 2 n t + C2
_____________________________________________________________
donde las constantes C1 y C 2 depende de las condiciones inciales de . x (0) y x (0) , es decir, para t = 0 se obtiene que:
x (0) = e 0 (C1 sen 0 + C 2 cos 0 )
x ( 0) = C 2
Para la derivada de x respecto a t, se tiene para t = 0 lo que sigue:
(40)
. t n C x (t ) = e 1 n e
( ( 1
nt
n cos 1 2 n t C 2 1 2 n sen 1 2 n t
( C1 sen 1 2 n t + C 2 cos 1 2 n t )
2
. x (0) = (1)C1
C1 =
( 1
n n (C 2 ) = C1 1 2 n x (0) n
. x (0) + x (0) n
n 1 2
(41)
Finalmente el valor de x es: Facultad de Ingeniera Arroyo M.
. nt x(0) + n x (0) 2 2 + x=e sen 1 t x ( 0 ) cos 1 t n n n 1 2
(42)
= donde d
= n 1 2
es la frecuencia de las oscilaciones.
Figura 6 Oscilacin amortiguada < 1.0 .
Arroyo M.
II.5 MOVIMIENTO NO OSCILATORIO Caso sobre-amortiguado
> 1.
La Solucin general est dada por:
x(t ) = Ae s1t + Be s 2 t
s1, 2 = n
(43)
2 n
( n )2
s1, 2 = 2 1 n
Como
(44)
> 1
las dos races son reales:
x(t ) = Ae (
2 1 nt
+ Be (
2 1 nt
(45)
donde las constantes depende de las condiciones inciales, es decir, para t = 0 se obtiene que:
x(0) = A + B
Para la derivada de x respecto a t, se tiene para t = 0 lo que sigue:
(46)
. x = A e (
2 1 nt
+ Be (
+ 2 1 n 2 1 ) n t 2 1 n
(47)
Arroyo M.
. 2 2 x ( 0) = A + 1 n + B 1 n
2 Si se multiplica la ec. 46 por + 1 n se obtiene:
(48)
2 2 x ( 0) + 1 n = A + 1 n 2 + B + 1 n
(49)
Restando a la ec.48 esta ltima ecuacin 49, se obtiene finalmente el valor de B:
. x ( 0 ) x ( 0) + 2 1 = B 2 1 + 2 1
n
[(
(
)(
n
)] n
= B 2 2 1
. 2 x (0) 1 n x ( 0) B= 2 2 1 n
(50)
De la ecuacin 46 se obtiene el valor de A:
x(0) = A + B
. 2 x ( 0) + 1 n x ( 0) A = x(0) B = x(0) + 2 2 1 n
(51)
Simplificando trminos se obtiene:
Arroyo M.
. 2 x ( 0) + 1 n x (0) A = x ( 0) + 2 1 2 n . 2 1 n + x ( 0) + 2 2 1 n
A =
x ( 0) 2
2 1 n x ( 0)
A =
. x ( 0) + + 2 n
2 1 n x ( 0) 21
(52)
Finalmente, se tiene la solucin remplazando los valores de A y B en la ec. 45.
x = Ae ( +
. x (0) +
2 1 nt
+ Be (
1 nt
x =
2 2
1 x (0) n 1
e ( +
2 1 nt
. 2 x ( 0) 1 x ( 0) n
2 2 1 n
e (
2 1 nt
(53)
Arroyo M.
Figura 7 Movimiento aperidico > 1.0 .
Arroyo M.
II.6 VIBRACIN EXCITADA ARMNICAMENTE La ecuacin diferencial de la Vibracin armnica forzada se puede escribir como:
.. . m x + c x + kx = F0 Sen t
(54)
Si se supone una solucin particular de la forma
x = Xsen( t )
se tiene:
x = X sen( t ) . x = X cos( t ) .. x = 2 X sen( t )

Remplazando estos valores en la ecuacin diferencial, ec. 54, se tiene:
(55)
m 2 X sen( t ) + c X cos( t ) + k Xsen( t ) = F0 Sen t

(56)
Si ( t ) = 0 , entonces t = . Substituyendo este valor en la ecuacin anterior, ec. 56, se tiene:
c X = F0 Sen
Si ( t ) = se tiene:
(57)
, entonces t =
+ . Substituyendo este valor en la ec. 56
m 2 X + kX = F0 Sen + 2
(k m )X
2
= F0 Cos
(58)
Arroyo M.
Con las ecuaciones 57 y 58, se plantea el sistema de ecuaciones siguiente.
c X = F0 Sen
(k m )X = F Cos
2 0
(59)
Elevando al cuadrado ambas ecuaciones y sumndolas se obtiene:
(k m ) X + (c )
2 2 2 2 2 2
X 2 = F02 Sen2 + Cos 2

2
)
(60)
((k m ) + (c ) )X
X = de donde: F0
2 2
= F02
(k m ) + (c )
(61)
Dividiendo ambos trminos de las ecs. 59 se tiene:
c X F0 Sen = = tan k m 2 X F0 Cos
tan =
c k m 2
(62)
Dividiendo entre k los trminos de los cocientes de las ec. 61 y 62, se tiene el valor de X y la fase, respectivamente:
X =
(k m ) + (c )
2 2
F0
1 F X = 0 2 2 k m 2 c 1 k + k
(63)
Arroyo M.
Y la fase como una expresin adimensional:
tan =
c k m 2
c cc 2mn c 2 cc k n k k tan = = = 2 = 2 m 2 2 1 1 2 1 1 k n n n
donde:
2 n =
(64)
cc = 2mn
k m es la frecuencia natural de oscilacin no amortiguada, c = el amortiguamiento crtico, ec.35, y cc un factor de
amortiguamiento. La amplitud expresada en forma adimensional es:
Xk = F0
1 m 2 c 1 k + k
2 2
1
2 2 (65) 2 1 2 + n n
La resonancia se calcula para
= 1: n
X=
F0 F = 0 2 k cn
Donde
c c 2 = 2 = k cc n 2 k = cn
c 2 = 2mn n
c = m 2 = c k n n
Arroyo M.
c c = c = n cc = 2mn 2m 2m
Figura 8 Representacin grfica de las ecuaciones 64 y 65 Finalmente, se tiene la ecuacin diferencial de la Vibracin armnica forzada:
.. . m x + c x + kx = F0 Sen t
.. c . k F x+ x+ x = 0 Sen t m m m
Haciendo
c c c = c = 2 n 2m 2m m
(36)
Arroyo M.
.. . F 2 x + 2n x + n x = 0 Sen t m
con la solucin completa (homognea y particular):
sen( t ) F x(t ) = 0 + X 1e n t sen 1 2 nt + 1 2 2 2 k 1 + 2 n n

donde la solucin particular es:
x = Xsen( t )
y la solucin homognea:
x = X 1e n t sen 1 2 nt + 1
. nt x (0) + n x(0) 2 2 x=e sen 1 t x ( 0 ) cos 1 t + n n n 1 2

sen( t ) F x(t ) = 0 2 2 2 k 1 2 + n n . nt x (0) + n x (0) 2 2 +e sen 1 n t + x(0) cos 1 n t n 1 2
Arroyo M.
II DINMICA ESTRUCTURAL Respuesta de sistemas estructurales de uno y de varios de libertad, ante excitacin ssmica. II.1 Sistema Estructural tipo marco de un grado de libertad, Caso 1.
Modelo de un sistema estructural tipo marco, de un grado de libertad.
Acelerograma19sep85EW
200
Aceleracin(cm/s2)
19 sep85 150 100 50 0 0 50 100 150 25 50 75 100 125 150
200
Tiempo (s)
Figura 5. Registros acerlerogrficos del sismo SCT-85
Arroyo M.
Magnitud 5.6, Guerrero
Arroyo M.
Historia de aceleraciones en la base
Cuando el sistema est sujeto a un movimiento en su base la masa estar en oscilacin y se generarn sobre ella tres tipos de fuerzas: a) La fuerza de inercia FI b) La fuerza elstica FE c) La fuerza de amortiguamiento FA Fuerza de inercia. De acuerdo con el principio de DAlambert, la fuerza de inercia es proporcional a la masa y a la aceleracin total que sta sufre, sea:
FI = m u T
(1)
donde: m es la masa y
uT
la aceleracin total igual a:
uT = u0 + u
(2)
donde u 0 , la aceleracin del terreno y respecto al terreno.
u,
la aceleracin de la masa con
Fuerza elstica. Es la fuerza que se genera en la columna por su rigidez lateral al tratar de ser desplazada respecto al terreno, sea:
FE = k u
donde:
k , es la rigidez lateral de las columna.
(3)
u, es el desplazamiento relativo de la masa respecto al desplante de la estructura (terreno). Fuerza de amortiguamiento. Es la fuerza que trata de restablecer el equilibrio de la estructura en vibracin, sea:
FA = c u
(4) Arroyo M.
donde: c, es el coeficiente de amortiguamiento, con relacin al terreno.
u , es la velocidad de la masa
La ecuacin de equilibrio dinmico se escribe como:
FI + FA + FE = 0
Sustituyendo las ecuaciones (1), (2) y (3) en la ecuacin (4), tenemos:
(4)
m u T +c u + ku = 0
como:
(5)
uT = u &&0 + u &&
, remplazando en la ecuacin (5) :
m u 0 + u + c u + ku = 0
m u + c u + ku = mu 0
dividiendo entre siguiente:
(6)
m,
resulta la ecuacin diferencial de equilibrio dinmico
c k u + u + u = u 0 m m
(7)
A partir de la ecuacin homognea de la ecuacin diferencial anterior, se tiene que:
u0= 0,
Arroyo M.
c k u + u+ u = 0 m m
(8) )
stantes las cons
c k y , represen ntan conc ceptos re elacionad dos con la l vibracin m m
libre del sistema, sea:
n =
k k 2 n = m m
(9) (
donde: n es la frecuenc cia circula ar (angula ar) del sis stema no amortigu uado, o se ea, aquella con c la que oscila el e sistema a cuando o se le imp pone un d desplaza amiento y se suelta. iguamien nto es nulo el siste ema describe un m movimient to armnico Cuando el amorti simple, con c la fre ecuencia
y con periodo o igual a:
Tn =
n
de la ec cuacin (9), resu ulta que el perio odo es:
Remplazando
Tn = 2
m k
(10)
Figura a 6. Vibra acin lib bre, Movimiento armnico a o simple. El amort tiguamien nto repres senta la disipacin d n de energa. Amortig guamient to crtico o. Facultad d de Ingeniera Arr royo M.
En vibracin libre amortiguada, es aquel para el cual el sistema despus de desplazado volvera a su posicin de reposo sin oscilar, sea:
ccr = 2 km
Constante de amortiguamiento.
(11)
La constante de amortiguamiento puede expresarse como una fraccin del crtico:
c ccr
c 2 km
k = mn , por tanto: m
(12)
Por otro lado:
km = m
c 2n m
c de la ecuacin (8), resulta: m
(13) (12.a)
La constante
c = 2n m
La ecuacin diferencial (7) puede, entonces, escribirse como:
c k u + u + u = u 0 m m
u + 2n u + u = u 0
2 n
(14)
Arroyo M.
La solucin de la ecuacin diferencial, cuando la estructura parte del reposo, se obtiene como la superposicin de la respuesta a una serie de impulsos diferenciales, en la forma llamada integral de Duhamel:
u(t ) =
vlida para:
n t
t
0
u0 ( )e(t )sen(t )d
(15)
1 2 1 , en estructuras de ingeniera civil < 0.1
Un caso particular de funcin excitadora continua, es:
Uo = a sen o t
(16)
donde: o , es la frecuencia circular de vibracin de movimiento del terreno y a es una constante que define la amplitud de la aceleracin. La solucin se presenta en la siguiente figura:
Respuesta ssmica de las estructuras
Arroyo M.
II.2 Siste emas Es structurales de va arios gra ados de libertad l
Caso cl sico que e puede idealizars se en es sta forma es un edificio, si imtrico, en que las masas s se consid deran con ncentrada as en ca ada piso y las col lumnas (los resortes) represe entan la ri igidez lat teral de c cada entre episo. En el equilibrio de e cada en ntrepiso interviene i en los mismos tipo os de fue erzas:
{FI }+ {FA }+ {FR } = 0

donde:
(17)
{FI } = [M ]{u &&T } {FA } = [C ]{u &} {FR } = [K L ]{u}

{FI }
[M ]
v vector de fuerzas de d inercia a, m matriz dia agonal de e masas,
siendo:
{u &&T } vector v de acelerac ciones tot tales, {FA } vector v de fuerzas de d amortiguamien nto,
[C ]
{u} vec ctor de de esplazam mientos.

{FE } [K s ]
{u &}
matriz de m e amortiguamiento o, usualm mente se conside era igual en to odos los entrepiso os y por tanto t se vuelve v un na consta ante, v vector de velocidad, v vector de fuerzas en e cada entrepiso e o, m matriz de rigidez del d sistem ma (no dia agonal),
Arr royo M.
Sustituyendo estas fuerzas en la ecuacin de equilibrio dinmico, se tiene:
[M ]{u &&}+ [C ]{u &}+ [K L ]{u} = [M ]{u &&0 }
(18)
II.2.1 Vibracin Libre. El primer paso, para la solucin de esta ecuacin diferencial es el caso de vibracin libre con amortiguamiento nulo; para determinar con buena aproximacin los periodos de vibracin y formas modales:
[M ]{u &&} + [K L ]{u} = {0}
(19)
Dado un impulso, la estructura vibrar libremente, adoptando una configuracin de desplazamientos que se denomina forma modal, con una frecuencia de vibracin peculiar de cada modo. El nmero de modos es igual al nmero de grados de libertad del sistema. En vibracin libre, el movimiento es armnico simple:
{u} = {a}sent
donde:
(20)
{u} es un vector de desplazamientos, {a} es un vector de constantes de amplitud.
Derivando 2 veces a la ecuacin (20), se tiene:
{u &} = {a}cos t
2 2 {u &&} = {a}sent = {u}
(21)
Sustituyendo en la ecuacin diferencial (19):
2 [M ]{a}sent + [K L ]{a}sent = {0}

Reordenando y eliminando Facultad de Ingeniera
(22)
sent , se tiene:
Arroyo M.
[[K ] [M ]]{a} = {0}

2 L
23)
La solucin no trivial es posible cuando:
[[K ] [M ]] = {0}
2 L
o sea, para que las vibraciones libres de amplitud finita (o diferentes de cero) sean posibles.
(24)
Resolviendo el determinante de la ec. 24 se obtiene un polinomio de orden n; dando solucin al polinomio, con mtodos numricos, se pueden conocer las n frecuencias correspondientes a los n modos naturales de vibracin del sistema.
Dados los n valores de se pueden determinar los n valores de a que corresponden a cada frecuencia natural y que determinan la forma modal correspondiente. La ecuacin (23) se puede representar como:
{}
[[K ] [M ]]{a } = {0}

L 2 i i
[Ei ]{ai } = {0}

donde la matriz
i
(25)
[Ei ] = [K L ] 2 [M ] , depende de la frecuencia i.
La forma del sistema vibrando puede ser resuelta para todos los desplazamientos en trminos de una cualquier coordenada:
a1i = 1 a i 2 a3 i {a} = . . . a ni
(26)
Arroyo M.
en forma explcita la ecuacin (25), queda expresada:

i i e11 e12 i i e21 e22 i i e31 e32 . . . . . . e i e i n1 n 2 i e13 i e23 i e33 .
. .
i en 3
. . . e1in 1 0 i a . . . e2 2 i n 0 i a3i 0 . . . e3 n . . . . . = . . . . . . . . . . . . . i . . . enn a2i 0
(26.a)
i e11 i E21 NGL 1,1
{ }
[ ]
0 1 = {a } i 2i NGL1,1 {0} E22 NGL 1, NGL 1

i 12 1, NGL 1
{E }
(26.b)
de donde se obtiene (o se deduce) que:
{E }
i 21 n 1,1
i + E22
[ ]
n 1, n 1
{a2i }n 1,1 = {0}n 1,1
(27.a) (27.b)
i i {a2i } = {0} e11 + E12
{ }
La ecuacin (27.a) se puede resolver simultneamente para las amplitudes de desplazamiento:
{E }
i 22
i 21 n 1,1
i + E22
[ ]
n 1, n 1
{a2i }n 1,1 = {0}n 1,1

(28.a)
[E ]{a
2i
i } = {E21 }
i 1 i ]n1,n1{E21 {a2i }n1,1 = [E22 }n1,1

La ecuacin (27.b) es redundante, por el hecho que el vector de frecuencias se satisface idnticamente.
(28.b)
Arroyo M.
Por conveniencia el vector de
{a} es expresado sin dimensin, al dividir todos

i :
los componentes por un componente de referencia (el ms grande aref = ai ,mx ). El vector resultante es llamado la forma modal i-sima,
1i 1 a i 2 2i 3i a3i 1 {i } = . = . a ref . . . . a ni ni
(29)
Arroyo M.
II.2.2 An nlisis de e la respuesta dinmica d de un sis stema es structura al de vari ios grados de libert tad, Mto odo de su uperposicin mo odal.
Para cua alquier co omponente modal vector de e la forma modal
{i } multtiplicado por la am mplitud modal
{ui }, los desplaz zamiento os estn dados d por r el

yi .
(30)
{ui } = {i }yi
el desplazamie ento total est dado d como la su uma de los componentes moda ales:
{u} = {1}y1 + {2 }y2 + {3}y3 + ... + {n }yn = {i }yi

i =1
Donde n son los grados de d libertad d. esin ante erior se puede p expresar en n forma matricial: m La expre
{u} = [ ]{y}
coordenadas gen neralizada as
(31) (
en este caso la a matriz de form mas moda ales [ ] sirve pa ara trans sformar las
yi , en n coordenadas ge eomtrica as u .
Arr royo M.
Gracias a la ortogonalidad de la matriz normal arbitraria
[ ] , para evaluar cualquier coordenada
{i }T , sea:
yi ; la ecuacin (31) puede ser reducida, premultiplicando por
{i }T [M ]{u} = {n }T [M ][]{y}
(32)
{i }T [M ]{u} = {i }T [M ]{1}y1 + {i }T [M ]{2 }y2 + ... + {i }T [M ]{n }yn (33)

Debido a la propiedad de ortogonalidad con respecto a las masas, todos los trminos desaparecen excepto el que corresponde al modo
i , {i } :
{i }T [M ]{u} = {i }T [M ]{i }yi

de donde:
{ i }T [M ]{u} yi = {i }T [M ]{i }
Sea la ecuacin de movimiento, sin amortiguamiento:
(34)
Desacoplamiento de la ecuacin de movimiento, sin amortiguamiento.
[M ]{u &&}+ [K L ]{u} = {p (t )}

- Sea:
(35)
La ecuacin (35) puede ser simplificada o desacoplada, aprovechando la propiedad de ortogonalidad de las coordenadas normales, como sigue:
{u} = [ ]{y} y {u y &&} = [ ]{& &}

Arroyo M.
- Remplazando estas equivalencias en la ecuacin de movimiento anterior, se tiene:
[M ][ ]{& y &} + [K L ][ ]{y} = {p (t )}

- multiplicando la ecuacin (36) por
(36)
{i }T , se tiene:
(37)
T T {i }T [M ][ ]{& y &} + {i } [K L ][ ]{y} = {i } {p (t )}
- Gracias a la ortogonalidad de la forma modal, se obtiene que:

T T {i }T [M ]{i }& y &i + { i } [K L ]{ i }y i = { i } {p (t )}
(38)
Ahora bien, sean las variables generalizadas definidas a travs de las equivalencias siguientes:
M i = {i } [M ]{i } , Masa generalizada de coordenadas normales,

T
K i = {i } [K L ]{i } , Rigidez generalizada, para el modo i , y

T
(39)
pi (t ) = {i } {p (t )} , Fuerza generalizada para el modo i .

T
Remplazando las variables generalizadas, ecs. 39, en la ecuacin 38, sta se puede representar como:
Mi& y &i + K i yi = pi (t )
(40)
que es una ecuacin de movimiento, de un grado de libertad, para el modo i . Por otro lado se tiene que, para vibracin libre:
[[K ] [M ]]{a} = {0}

2 L
Arroyo M.
[K L ]{i } = i2 [M ]{i }
T
(41)
i } , se tiene: Pre-multiplicando ambos lados de la ecuacin (41) por {
{i }T [K ]{i } = i2 {i }T [M ]{i }
K i = i2 M i
(42)
La expresin (42) relaciona la rigidez generalizada para el modo i con la masa generalizada a travs de la frecuencia de vibracin, por lo anterior se tiene:
i2 =
Ki Mi
De acuerdo a lo anteriormente expuesto, la ecuacin (40) se pude simplificar, expresndose como:
Mi& y &i + K i yi = pi (t ) y & &i + Ki 1 yi = pi (t ) Mi Mi 1 2 y + y = pi (t ) &&i i i Mi
(43)
Este procedimiento es utilizado para obtener una ecuacin independiente de un grado de libertad, para cada modo de vibrar de la estructura. O sea, el uso de la coordenada normal (generalizada) sirve para transformar las ecuaciones de movimiento de un conjunto de n ecuaciones diferenciales simultneas, que estn acopladas, dada la matriz de rigideces llena, es decir, con trminos fuera de la diagonal, a un sistema de n ecuaciones independientes de coordenadas normales y de un grado de libertad. La respuesta dinmica puede ser, entonces, obtenida resolviendo separadamente para la respuesta de cada coordenada (modal) normal y despus Facultad de Ingeniera Arroyo M.
superponindolas de acuerdo a la ec. 31, para obtener la respuesta en las coordenadas originales. Este procedimiento es llamado el Mtodo de superposicin modal II.2.3 Vibracin de un sistema estructural excitado, con amortiguamiento. Sean las ecuaciones desacopladas de movimiento, con amortiguamiento:
[M ]{u &&}+ [C ]{u &} + [K L ]{u} = {p(t )}

[ ]
n } , esta ecuacin (44) se transforma en: modal {
T
(44)
coordenadas generalizadas y desplazamientos (ec. 31): {u} = {y} , as como sus derivadas, adems, pre-multiplicando por la traspuesta de la n-sima forma
Nuevamente, introduciendo en la ecuacin 44, la equivalencia entre
T T T {i }T [M ][ ]{& y &} + {i } [C ][ ]{y & } + {i } [K L ][ ]{y} = {i } {p (t )} (45)
Por otro lado, se tienen las condiciones de ortogonalidad:
{ } [M ]{ } = 0 { } [K ]{ } = 0
T j i T j L i
para
j i,
(46)
Asumiendo que esta condicin de ortogonalidad tambin es vlida para el amortiguamiento:
{ } [C ]{ } = 0 para
T j i
ji
(47)
Por consiguiente, la ecuacin (45) puede re-escribirse como:

T T T {i }T [M ]{i }& y &i + {i } [C ]{i }y & i + {i } [K L ]{i }yi = {i } {p (t )}
Mi& y &i + Ci y & i + K i yi = pi (t )

y & &i +
Ci K 1 y pi (t ) & i + i yi = Mi Mi Mi
(48)
donde:
1 Ki {i }T [K L ]{i } = i2 = Mi Mi Ci 1 {i }T [C ]{i } = 2i = Mi Mi
pi (t ) = {i } {p (t )}
T
(49)
Dividiendo entre M n e introduciendo estas nuevas igualdades en la ecuacin (48), se tiene:
y & &i +
Ci K 1 y pi (t ) & i + i yi = Mi Mi Mi 1 pi (t ) Mi
(51)
2 y & &i + 2 i y & i + i yi =
Arroyo M.
II.2.4 Resumen del mtodo o procedimiento de superposicin modal. La transformacin de coordenada normal, es la base del mtodo de superposicin modal del anlisis dinmico. Este mtodo puede ser usado para evaluar la respuesta dinmica de cualquier estructura lineal, para el cual los desplazamientos han sido expresados en trminos de un conjunto de n coordenadas discretas y donde el amortiguamiento puede ser expresado por relaciones de amortiguamiento modal. Paso1. Ecuacin de movimiento. Para esta clase de sistema, las ecuaciones de movimiento pueden ser expresadas como:
[M ]{u &&}+ [C ]{u &} + [K L ]{u} = {p(t )}

Paso 2. Forma modal y anlisis de frecuencia.
(21)
Para vibraciones libres sin amortiguamiento, esta ecuacin matricial puede ser reducida por la ecuacin (de valores y vectores caractersticos) caracterstica:
[[K ] [M ]]{a} = {0}

2 L
(22)
A partir de la cual la matriz de formas modales de vibracin frecuencia puede ser determinado.
, y el vector de
Paso 3. Masa, rigidez, amortiguamiento y carga generalizadas. Con cada vector de forma modal n , la masa, la rigidez, el amortiguamiento y la carga generalizados, pueden transformarse en constantes generalizadas:
M i = {i } [M ]{i }
T
K n = {i } [K L ]{i } = i2 M i
T
Cn = {m } [C ]{n } = 2 n n M n
T
Arroyo M.
pi (t ) = {i } {p (t )}
T
n = 1, 2, ..., N
(23)
Paso 4. Ecuaciones desacopladas de movimiento. La ecuacin de movimiento para cada modo puede, entonces, ser reescrita usando la masa, rigidez, amortiguamiento y fuerza generalizados, junto con la frecuencia modal modal
n , y un valor especificado de la relacin de amortiguamiento
, como sigue:
2 y & &n + 2 n y & n + n yn =
1 pn (t ) Mn
(24)
donde:
n =
c c c = = ccr 2n m 2 km
Paso 5. Respuesta modal para la carga. El resultado de la etapa 4 es un conjunto de n ecuaciones independientes de movimiento, una para cada modo de vibrar. Estas ecuaciones de un grado de libertad pueden ser resueltas por cualquier mtodo que se quiera, dependiendo del tipo de carga. La expresin de respuesta general dada por la integral de Duhamel para cada modo es:
yn (t ) =
donde: tanto:
1 M n Dn
t0
pn ( )e n n (t ) sen Dn (t )d
(25)
Dn = n 1 2
, si < 0.1 , para fines prcticos:
Dn = n ,
por
Arroyo M.
yn (t ) =
t 1 n n (t ) ( ) p e senn (t )d n t M nn 0
Paso 6. Vibraciones libres modal. La ecuacin anterior es aplicable para un sistema que est en reposo en el tiempo t = 0 . Si la velocidad inicial y el desplazamiento inicial no son cero, una respuesta de vibracin libre debe ser aadida a la expresin de la integral de Duhamel para cada modo. La respuesta de vibracin libre amortiguada esta dada, para cada modo, por:
yn (t ) =
y & ennt
(0) + n n yn (0)
Dn
sen D n t + y n (0) cos D n t (26)
Donde
yn (0 ) y y & n (0 ) representan el desplazamiento y velocidad modal inicial (en el tiempo cero). Estos pueden ser obtenidos del desplazamiento vn (0 ) y &n (0) inciales, especificados y expresados en las coordenadas velocidad v
geomtricas originales como sigue, para cada componente modal:
{ n }T [M ]{v(0)} {n }T [M ]{v(0)} = yn (0) = T Mn {n } [M ]{n }

T { n }T [M ]{v &(0)} {n } [M ]{v &(0)} = y & n (0) = Mn {n }T [M ]{n }
(27)
Paso 7. Respuesta. Desplazamiento en coordenadas geomtricas.
Arroyo M.
Cuando la respuesta para cada modo yn (t ) ha sido determinada de la ecuacin (25) y/o (26); el desplazamiento expresado en coordenadas geomtricas esta dado por la transformacin de coordenadas normales:
{v(t )} = []{y(t )}
(28)
{u(t )} = {1}y1 (t ) + {2 }y2 (t ) + {3 }y3 (t ) + ... + {n }yn (t )

Es importante hacer notar que para la mayora de tipos de carga la contribucin de varios modos es generalmente ms grande para frecuencias bajas y tiende a decrecer para las frecuencias ms altas. Por consiguiente usualmente no es necesario intuir todos los modos ms altos de vibracin en el proceso de superposicin. Paso 8. Respuesta de fuerza elstica. La historia de desplazamiento de la estructura puede ser considerada como medida bsica de su respuesta ante carga dinmica. En general otros parmetros de respuesta tal como los esfuerzos o fuerzas desarrolladas en varios componentes estructurales pueden ser evaluados directamente de los desplazamientos. Por ejemplo las fuerzas elsticas f s (t ) , las cuales resisten la deformacin de la estructura, estn dadas directamente por:
f s (t ) = [K L ]{u (t )} = [K L ][ ]{y (t )}
(29)
f s (t ) = [K L ]{1}y1 (t ) + [K L ]{2 }y2 (t ) + ... + [K L ]{n }yn (t )

y sustituyendo como sigue:
[K L ]{n } = n2 [M ]{n } , la ecuacin (29) se puede reescribe

(30)
2 f s (t ) = 12 [M ]{1}y1 (t ) + 2 [M ]{2 }y2 (t ) + ... + n2 [M ]{n }yn (t )
expresando la ecuacin (30) en forma matricial, se tiene: Facultad de Ingeniera Arroyo M.
2 f s (t ) = [M ][ ] n yn (t )
(31)
donde n yn t representa un vector de amplitudes modales multiplicadas respectivamente por su frecuencia modal al cuadrado.
2
( )}
En la ecuacin 31, la fuerza elstica asociada con cada componente modal (ec. 29), ha sido reemplazada por una expresin equivalente de fuerza-inercial-modal.
Arroyo M.
RESPUESTA SSMICA DE SISTEMAS LINEALES ESTRUCTRUALES Ecuacin de movimiento:

2 u && + 2 n u & + n u = u &&g
Una vez que se ha evaluado la respuesta de deformacin realizando un anlisis dinmico de la estructura, las fuerzas internas pueden determinarse mediante un anlisis esttico del sistema estructural, para cada instante de tiempo t , utilizando las fuerzas estticas equivalentes:
f s (t ) = ku (t )
Expresada esta ecuacin en trminos de la masa se tiene:
f s (t ) = m 2u (t ) = mA(t )
Donde A t es llamada seudo-aceleracin o aceleracin espectral del sistema, cuya respuesta puede ser calculada a partir de la respuesta del desplazamiento
()
u (t ) .
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.
I. ANLISIS TRIDIMENSIONAL SIMPLIFICADO PARA EDIFICIOS SOMETIDOS A CARGAS LATERALES Teniendo como base el mtodo de rigideces se puede realizar un anlisis simplificado para estructuras en el espacio del tipo edificio. Los edificios se pueden idealizar como un conjunto de subestructuras (marcos y muros) planas y verticales unidas entre s mediante los sistemas de piso (figura 1). Los sistemas de piso, paralelas entre s y generalmente conformados con losas, tienen como caracterstica principal que pueden servir como un diafragma en su entrepiso, es decir, que la rigidez de la losa es muy grande comparada a la de las vigas que la rodean por lo puede ser considerada como un cuerpo rgido (indeformable) en su plano, y por lo tanto las deformaciones de los elementos dependern de los desplazamientos de las losas.
Trabes C Columnas
1 2 3
Losa PLANTA FIGURA 1 Sistema Estructural tipo edificio. X
Las subestructuras que conforman el edificio pueden ser compuestas por vigas, columnas, muros de cortante, contraventeos, etc., para los cuales el mtodo se puede utilizar, pero algunos elementos tienen que tomarse en una forma especial, que queda fuera del alcance de este trabajo, por lo cual se limitar a marcos formados por vigas y columnas. Este mtodo permite realizar anlisis de edificios sometidos a carga lateral nicamente, obteniendo como resultado la distribucin de las fuerzas laterales en cada uno de los marcos planos que conforman el edificio, as como el desplazamiento lateral de stos en la direccin de su plano.
Arroyo M.
Las cargas laterales a las que se someten las estructuras pueden ser de diferente origen, sismo, viento, empuje de tierra, o de agua, oleaje, etc. En edificios por lo general son las dos primeras, siendo las cargas producidas por sismo las que tomaremos en cuenta en este trabajo por ser un tipo de carga sumamente significativa en varias regiones de nuestro pas y del mundo; para la aplicacin del mtodo se deben utilizar las cargas ssmicas obtenidas de acuerdo al mtodo esttico propuesto por el Reglamento de Construcciones del D.F. Las hiptesis fundamentales en las que se basa el mtodo son:
- El anlisis es del tipo Elstico Lineal - Los desplazamientos se consideran pequeos. - Las losas de los niveles se consideran infinitamente rgidas ante cargas horizontales y se desplazan como un cuerpo rgido, considerando tres grados de libertad por entrepiso (figura 2).
Z y z x
FIGURA 2 Grados de libertad en cada nivel de entrepiso o losa.
- Los marcos se modelan y analizan como marcos planos en dos dimensiones, y estn formados exclusivamente por trabes (vigas) y columnas, y cada uno de stos trabaja en forma independiente en su eje (figura 3).
Arroyo M.
FIGURA 3 Marco plano.
- Se desprecia la rigidez torsional de trabes normales a los marcos. - Las masas se concentran en los niveles y en el centro de masa de cada uno se aplican las fuerzas ssmicas (figura 4).
Fx
C.M. Fy
2 C.M.: Centro de Masa 3
FIGURA 4 Centro de masa en planta.
Arroyo M.
I.1 Matriz de rigidez lateral de un marco en 2d.

El primer paso para la elaboracin del anlisis es obtener la matriz de rigidez lateral de cada uno de los marcos planos que forman la estructura, lo cual se realiza partiendo de la matriz de rigidez de cada marco.
F n Fi F
2
n i 2 1
F 1
FIGURA 5 Marco K j . El procedimiento implica ensamblar convenientemente la matriz de rigidez de cada uno de los marcos K j , que forman el edificio, de tal manera que las rigideces ante los desplazamientos horizontales ocupen las ltimas posiciones de dicha matriz. As se puede obtener la matriz simplificada, condensndola de la manera como a continuacin se deduce. Sea la ecuacin de equilibrio esttico del marco j:
{F } j
{F } {0} [K 11 ] = [K ] j {d } j = 1 = = {F2 } {F2 } [K 21 ]
[K 12 ] { 1 } , [K 22 ] { 2 } lat
(1)
donde las nicas fuerzas que actan en la estructura son las fuerzas horizontales {F2 } . Desarrollando la ecuacin 1, se tiene:
[K11 ]{1} + [K12 ]{ 2 }

= {0}
(2)
Arroyo M.
[K 21 ]{1} + [K 22 ]{ 2 }
{1}
{F2 }
(3)
Despejando {1} de la ecuacin 2, se tiene:
= [K11 ] [K12 ]{ 2 }
1
(4)
Sustituyendo la ec. 4 en la ec. 3, resulta:

[K 21 ][K11 ] [K12 ]{ 2 } + [K 22 ]{ 2 } = {F2 }
1
[ [K
dado que:
21
][K11 ]1 [K12 ] + [K 22 ]]{ 2 } = {F2 }
[F ] j = [K L ] j [ ] j
KL j
[ ]
nn , nn
= [K 22 ] [K 21 ][K11 ] [K12 ]
1
(5) en donde [K L j ] es la matriz de rigidez lateral del marco j.
I.2 Relacin entre desplazamientos horizontales del marco j y desplazamientos de entrepiso.

Los desplazamientos horizontales (en cada nivel i) de cualquier marco j, en direccin a su sistema de ejes locales, se pueden obtener en funcin de los desplazamientos de las losas referenciados a un sistema global. De la figura 6 se pueden obtener los valores del desplazamiento en cada nivel i del marco j dando como resultado:
ij = cos j xi + sen j y i + r j i
{ i } j
(6)
= cos j { xi } + sen j { y i } + r j { i }
Arroyo M.
y i sen
Y
yi
Pj (x j, y )
j
+ Marco r j j
xi x i cos
j
i
j
yi
rj
C.M.
xi
i
X
FIGURA 6 Desplazamientos del entrepiso y marco j. Ahora bien, dado que el valor de j no cambia para cada entrepiso y los desplazamientos en cada entrepiso dependen nicamente del valor de los desplazamientos de su losa se puede generalizar la expresin anterior para cada marco j como:
{}
(7) donde:
[[cos] [sen]
{ x } [r ]] j y } = [T ] j { } { { }
ne, es el nmero de entrepisos del marco.

[ cos] , es una matriz diagonal con elementos iguales a cos j , de orden ne x ne. [ sen] , es una matriz diagonal con elementos iguales a sen j , de orden ne x ne.
[r ] , es una matriz diagonal con elementos rj (distancia euclidiana del marco al

origen), de orden ne x ne.
{ x }, { y } y {} , son vectores de desplazamientos de la losa de orden ne.

{T}j , es la matriz de transformacin del marco j, unin de las anteriores, y por lo tanto de orden ne x 3ne.
{ } j , es el desplazamiento horizontal del marco j en su direccin, de orden ne x

1. I.2.1 Obtencin de r j . En forma grfica, de la figura 7 se tiene:
x jsen j
Pj (xj ,y j)
MARCO j
R A
yj cos j
d
Y
G C
rj
O
yj
xj
FIGURA 7 Obtencin de la distancia r.
y j cos j = x j sen j + d
De la figura 7 se tiene: Del triagulo ROG: Del triagulo ABC:
d = x j sen j y j cos j
r j = b cos j
d = b cos j rj = d
r j = x j sen j y j cos j
(8)
Arroyo M.
De forma anlitica se puede obtener tambin rj . Como el punto Pj x j , y j pertenece a la recta del eje del marco j, se tiene la relacin siguiente:
yj =m xj + b
sen j m = tan = j Donde cos j
Figura 8 Eje del marco j. Del triangulo rectngulo ROG se obtiene el valor de b:
(a)
(b)
cos j =
(c)
rj b
b=
rj cos j
Substituyendo las ecuaciones (b) y (c) en (a) se tiene:
yj =
sen j cos j
xj +
rj cos j
Despejando rj de la ecuacin anterior:
rj = x j sen j y j cos j
Representando el valor de en valor absoluto, se obtiene:
rj = x j sen j y j cos j
Arroyo M.
I.3 Obtencin de la matriz de rigidez tridimensional simplificada.

A partir de la matriz de rigidez lateral de cada marco es posible ensamblar la matriz de rigidez tridimensional simplificada, la cual permite resolver el sistema de ecuaciones que relaciona las fuerzas laterales que actan en la estructura con los desplazamientos de cada entrepiso (en ambas direcciones y el giro), lo que nos permite obtener tanto los desplazamientos laterales en la direccin de cada marco, as como el valor de la fuerza lateral que debe soportar. Por Ley de Hooke, las fuerzas horizontales aplicadas en cada uno de los marcos j, son:
{Fi }j
= [K L ] j { i }j = [K L ] j [T ] j { }
x1 1 y1 2 y2 x2 . + sen j . + rj . . . . n yn xn
(9)
F1 F 2 . = [K L ] j cos j . .Fn j
Por Equilibrio Esttico las fuerzas totales aplicadas en las losas deben ser la sumatoria de las fuerzas de cada marco, tomando en cuenta las componentes en planta de stas:
{Fi } = [T ]Tj {Fi }j

j =1
nm
(10)
donde: nm es el nmero de marcos. Sustituyendo la ecuacin (4.29) en la ecuacin (4.30), se tiene:
{Fi }
donde: [K T ] = simplificada. Facultad de Ingeniera
nm T j
[T ] [K ] [T ] { } = [K ] { }
j =1 T j L j j T j
j j
nm
(11)
[T ] [K ] [T ]
j =1 L
es
la
matriz
de
rigideces
tridimensional
Arroyo M.
Con la finalidad de observar cada uno de los trminos e ir armando la matriz de rigidez tridimensional simplificada del sistema se efectuarn las operaciones correspondientes: La porcin de {Fx} resulta:
{Fxi } = {Fi } j cos j

j =1
nm
Tomando en cuenta la ecuacin (9), se tiene:
{Fx } = [K L ]j {cos2 j { x } + cos j sen j { y }+ cos j rj { }}

mn j =1
[K xx ]
haciendo:
[K ] {cos a }
nm 2 j =1 L j j nm j =1 nm L j j j
[K xy ] = [K ] {cos a sena }
[K xr ]
=
[K ] {cos a r }
j =1 L j j j
entonces:
{Fx }
{Fx }
= [K xx ]{ x }+ K xy y + [K xr ]{ }
=
[ ]{ }
[ [K xx ] + [K xy ]+ [K xr ]] { }
{ x } y } {
La porcin de {Fx} resulta:
{F } = {F } sen
nm y j =1 i j
Tomando en cuenta la ecuacin (9), se tiene:
Arroyo M.
{F } = [K ] {cos
nm y j =1 L j
sen j { x } + sen 2 j { y }+ sen j r j { }
Haciendo:
[K ] = [K ] [K ] = [K ] sen
yx xy nm 2 yy j =1 L j nm yr j =1 L j
[K ] = [K ] sen r
por lo tanto:
y yx x
j j
{F } = [K ]{ } + [K ]{ }+ [K ]{ }
yy y yr
{F } = [[K ] + [K ] + [K ]]
y yx yy yr
{ x } y } { } {
La parte del momento {Mo} resulta:
{M o } = {Fi }j rj
j =1
nm
Con la ecuacin (9), se tiene:
{M o }
[K ] {cos r { } +
nm j =1 L j j j x
sen j rj { y } + rj { }
2
[K rx ] = [K xr ]
Haciendo:
[K ] = [K ]
ry yr nm j =1
[K rr ] = [K L ]rj 2
por lo tanto:
{M o } = [K rx ]{ x } +
[K ]{ } + {K }{ }
ry y rr
El sistema de ecuaciones, en forma matricial, resulta:
Arroyo M.
{Fx } [K xx ] {Fy } = K yx [K rx ] {M o }
{ } [K ] [K ] [ ] [K ] [K ] { } [K ] [K ] } {
xy xr x yy yr y ry rr
{FS } = [K T ]{ S }
donde [K T ] es la matriz de rigideces tridimensional simplificada del sistema estructural.
[K xx ]
[K ] {cos a }
nm 2 j =1 L j j nm j =1 nm L j j j
[K xy ] = [K ] {cos a sena }
[K xr ]
yx
[K ] {cos a r }
j =1 L j j j
[K ] = [K ] [K ] = [K ] sen
xy nm 2 yy j =1 L j nm yr
[K ] = [K ] sen r
[K rx ] = [K xr ]
j =1 L j
j j
[K ] = [K ]
ry yr nm j =1
[K rr ] = [K L ]rj 2
Arroyo M.
I.4 Desplazamientos y fuerzas ssmicas en cada marco.

Con la matriz de rigidez tridimensional se puede resolver el sistema con las fuerzas aplicadas sobre las losas, resultantes del mtodo esttico. Cabe sealar que el trmino correspondiente a los momentos sobre las losas Moi se obtiene con respecto al origen del sistema global de referencia, tomando en cuenta que las fuerzas ssmicas se aplican en el centro de cortante, la evaluacin tanto del centro de cortante como del trmino de momento se explicar con mayor detalle ms adelante. Dando solucin al sistema de ecuaciones, se obtienen los desplazamientos de las losas de cada entrepiso con sus tres componentes:
x1 x2 = . , . xn y1 y 2 = . . y n 1 2 { } = . . n
{ x }
{ }
y
Con estos desplazamientos se obtienen los desplazamientos laterales de cada marco en su direccin mediante:
{}j
= cos j { x } + sen j { y } + rj { }
Finalmente, se obtienen las fuerzas laterales que deben resistir cada uno de los marcos con:
{F }
= [K L ] j {}j
I.5 Centros de torsin, de cortante, de excentricidades y de las fuerzas de diseo.
las
Para la obtencin de las fuerzas de diseo de acuerdo al reglamento del Distrito Federal de 1987, es necesario determinar la excentricidad entre el centro de torsin y el centro de cortante, para poder obtener el nuevo centro de cortante y obtener los momentos de diseo para el anlisis tridimensional.
Arroyo M.
I.5.1 Centro de Torsin Considerando que las fuerzas estn aplicadas en el centro de torsin para no provocar giros en las losas, se puede plantear lo siguiente: A) Si el sismo acta en la direccin X, se tiene que:
{Fx } {0},
{F } = {0}
y
y { } =
{0}
El sistema de ecuaciones de equilibrio, resulta:
[K xx ] {Fx } = K {0} yx
[K ] { } [ ] [K ] } {
xy x yy y
Dando solucin al sistema se obtienen los desplazamientos {x} y {y}, con los cuales se determinan los momentos que hay que aplicar a la losa para que no gire:
{M o } = [K rx ]{ x } + [K ry ]{ y }
La distancia YCTi del centro de torsin para cada nivel y, se obtiene como:
YCTi =
M oi Fxi
b) Si el sismo acta en la direccin Y, se tiene:
{Fx }
= {0},
{F } {0}
y
{ }
= {0}
El sistema de ecuaciones de equilibrio ahora resulta:
[K xx ] {0} = K {F } yx y
[K ] { } } [ ] [K ] {
xy x yy y
Similarmente al inciso (a) se determina el valor del momento con los valores de {x} y {y}, para que la losa no gire:
Arroyo M.
{M o }
= [K rx ]{ x } + K ry { y}
[ ]
La distancia XCTi del centro de torsin para cada nivel y, se obtiene como:
X CTi =
M oi Fyi
I.5.2 Centro de masa (CM). El centro de masa para cada nivel i, se obtiene por primeros momentos, y normalmente coincide con el centro de rea, pero existen casos en que hay masas concentradas en algunos puntos de algn nivel. I.5.3 centro de cortante (CC). Conocido el centro de masa de cada nivel i, se puede conocer el centro de cortante como sigue:
X CC i =
F X
j =i j
nn
CM j
Fj Vi
Vi
Vi 8. Donde nn, es el nmero de Figura entrepiso. niveles

I.5.4 Excentricidades
YCC J =
F Y
j =i
nn
j CM j
Cortantes
de
Calculadas entre el centro de torsin y el centro de cortante:
e x i = X CCi - X CTi ,
e y i = YCCi YCTi
Para diseo, stas deben incrementarse para tomar en cuenta un incremento de excentricidad tanto por una distribucin accidental de las masas como por razones dinmicas. De acuerdo al reglamento del Distrito Federal de 1987 los valores de diseo son: Facultad de Ingeniera Arroyo M.
Para X:
Para Y:
exd1i = 1.5exi + 0.1bx
e yd 1i = 1.5e yi + 0.1b e yd 2 i = e y i 0.1b
exd 2i = exi 0.1b

donde:
- las excentricidades estn medidas con respecto al centro de torsin y de acuerdo a la convencin del signo del momento. - b, es la mayor dimensin, en planta, del nivel y medida perpendicularmente a la direccin del mismo.
I.5.5 Coordenadas del centro de cortante referidas al origen, para cada excentricidad de diseo Dado que tanto para la obtencin de las excentricidades se obtiene la diferencia de las coordenadas del centro de torsin de las de cortante, y stas se incrementan para diseo de acuerdo al R.D.F., para la realizacin del anlisis es necesario obtener las nuevas coordenadas del centro de cortante, tomando en cuenta las nuevas excentricidades de diseo:
X 'CC 1i = X CTi + exd 1i , X 'CC 2 i = X CTi + exd 2 i , Y 'CC1i = YCTi + e yd 1i
Y 'CC 2 i = YCTi + e yd 2 i
I.5.6 Momentos de diseo, para cada excentricidad de diseo En la direccin X: En la direccin Y:

1i Vi X CC 1i +1 Vi +1 M D1i = X CC 2 i Vi X CC 2 i +1 Vi +1 M D 2 i = X CC
1i V i YCC 1i +1 V i +1 M D1i = YCC 2 i V i YCC 2 i +1 V i +1 M D 2 i = YCC

I.5.7 fuerzas resultantes de diseo En la direccin X:
En la direccin Y:
Arroyo M.
{Fx } {0} {Fy } = {0} M {0} o1
{ }
xd1
{Fx } {0} , {Fy } = {0} M {0} o2
{ }
xd 2
{Fx } = {0} {F y } {0} M {0} o1
{ }
yd1
{Fx } = {0} , {F y } {0} M {0} o2
{ }
yd 2
Hasta aqu es la contribucin del sismo en direcciones independientes.
Arroyo M.
FUERZAS SSMICAS CASO ESTTICO De acuerdo al modelo esttico propuesto por el reglamento, las fuerzas y cortntes que actan en una estructura tipo marco se pueden obtener como sigue.
Se supone una distribucin lineal de aceleraciones de respuesta, en los diferentes niveles de la estructura, partiendo de una aceleracin cero en la base hasta un mximo
am
en la punta.
1) La fuerza en cualquier entrepiso o nivel de piso es:
Fi = mi ai =
Donde
wi ai g
mi
(a) es la masa total del nivel i,
wi
es el peso total del nivel i,
g la
gravedad y
ai
la aceleracin del nivel i.
2) A partir de la variacin de la aceleracin se obtiene por triangulos semejantes la siguiente expresin:
Arroyo M.
ai a h = m ai = i am ghi gH H
Donde
H es la altura total del edificio y hi la altura del nivel i.
Remplazando esta aceleracin del nivel i en la ecuacin (a) se tiene que:
Fi = mi ai =
wi hi am g H
(b)
3) El cortante en la base se obtiene sumando los efectos de las fuerzas de cada uno de los niveles, como sigue:
V = Fi =
wi hi a am = m g H gH
wh
i i
(c)
4) La aceleracin mxima a m se puede obtener en funcin del cortante, despejndola de la ecuacin (c):
am =
VgH wi hi
(d)
5) Substituyendo la ecuacin (d) en la ecuacin (b), se tiene la fuerza ssmica del nivel i:
Fi =
wi hi w h VgH am = i i = g H g H wi hi
wi h V wi hi
(e)
6) El cortante en la base V se obtiene de acuerdo al modelo esttico propuesto por el Reglamento de Construcciones como:
V=
C W Q'
Arroyo M.
C es un coeficiente ssmico, Q' un factor de reduccin que depende del factor de comportamiento Q y W = wi el peso total del edificio.
Donde, 7) Finalmente, las fuerzas ssmicas que se aplican en los centros de masa de cada nivel se pueden obtener como:
wi hi wi hi C CW Fi = V= Q ' W = Q ' w h wi hi w h w h i i i i i i = C wi wi h ; Q ' wi hi c a0 Q'
Arroyo M.
Arroyo M.
Arroyo M.
Arroyo M.

Diseño Estructural 2012

Cargado por

Copyright:

Formatos disponibles

Diseño Estructural 2012

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Diseño Estructural 2012

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSIDAD AUTNMA DE QUERTARO

Facultad de Ingeniera Ingeniera Civil

DR. GUADALUPE MOISS ARROYO CONTRERAS

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSID DAD AUTON NOMA DE QU UERETARO

AP PUNTES DE DISEO D ESTRUCTURAL

UNIVERSID DAD AUTON NOMA DE QU UERETARO

AP PUNTES DE DISEO D ESTRUCTURAL

UNIVERSID DAD AUTON NOMA DE QU UERETARO

AP PUNTES DE DISEO D ESTRUCTURAL

x0 , por r tanto, va ale:

d dx d dx P x0 < X < x0 + = f x ( x0 )dx 2 2

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

- El momento de inercia del rea bajo la curva, o la variancia, o segundo momento:

Se llama desviacin estandar a la raiz cuadrada de la variancia:

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

Funcin de distribucin acumulada

UNIVERSID DAD AUTON NOMA DE QU UERETARO

AP PUNTES DE DISEO D ESTRUCTURAL

x mx x0 mx FX ( x0 ) = P[ X < X 0 ] = P < x x x0 mx FX ( x0 ) = P[u < u0 ] = = (u0 ) x

UNIVERSID DAD AUTON NOMA DE QU UERETARO

AP PUNTES DE DISEO D ESTRUCTURAL

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

USANDO TABLAS SE PROCEDE COMO SIGUE:

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

Estadisticamente, se tiene que: Media de una muesta:

Variancia de una muestra:

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

FUNCIN DE DENSIDAD DE PROBABILIDADES CONJUNTA DE DOS VARIABLES ALEATORIAS X e Y

UNIVERSIDAD AUTONOMA DE QUERETARO

APUNTES DE DISEO ESTRUCTURAL

Var ( g ( x )) = ai2 Var [X i ]

Var ( g ( x )) CVx2 + CVx2 + CVx2 1 2 3 E (g ( x ))

UNIVERSID DAD AUTON NOMA DE QU UERETARO

AP PUNTES DE DISEO D ESTRUCTURAL

CARGA A VIVA EN N EDIFIC CIOS

UNIVERSID DAD AUTON NOMA DE QU UERETARO

AP PUNTES DE D DISEO ESTRUCTURAL

UNIVERSID DAD AUTON NOMA DE QU UERETARO