— UNIVERSIDAD CATOLICA facultad de ingenieria civil . 0,98 ‘feos -0.98 ‘ @= 0.40 “ANALIS| D. EDIFI ao EL CODIGO (DECRETO 1400) Qure (ricer Caen oF . CARLOS RAMIRO VALLECILLA B septiembre i6 DE 1993ANALISIS DINAMICO DE UN EDIFICIO DE TRES PI PROFESOR: CARLOS VALLECILLA B. OBER IVO. Analizar en la direccién X el edificio mostrado 25m [30/20 DATO! Ubicacién: Cali CARGAS: ‘25m Cubierta: 0.5 ton/m* t Pisos: 0.7 ton/m? sored ra conereto: {i £7c=210 kg/cm B=13000/210=188387ke/em* =2.000.000 ton/m* PROCEDIMIENTO 1) Obtencién de la masa por piso. cubierta: W=0.5*10#11=55 ton ; M3 = 55/9.81 =5.60 ton - seg*/m piso: W=0.7*10*11=77 ton 2)0btencion de los coeficientes de rigidez. — i . weeV 3 Se supone que Ja inercia de la e placa de entrepiso es mucho més grande que la inercia de las co- lumnas. Se supone ademas que los pérticos paralelos al eje X son iguales. : [Link] piso 130.940.4712 =.0016 mt Gols. 2 y 3 piso 1=0,84/12' =0.000875 m TINA 4 =T1/9.81 =7.85 ton-ses* /m rio seve =~ eeePRIMER PISO ieEt/P =1422 tom/m -— KO a 3 ne 1get/t =1037 tom/m 2 YO VP . Ja Ka 1037 1037 1037, — — at ek, £7377 —— oe, — —_— — 7 1442 tae, eae VIGA- PRIMER PISO SEGUNDO PISO PAGINA 2La matriz de rigidez de edificio es: 22131-9333 ° Ke -9333 sess -9333 0 -9333 ~ 9333 3)Ecuacién del movimiento La ecuscién que describe e1 movimiento es isin amortigua- miento): (4G) + HOI= es] ecuacién (a) es decir: 7.88 0 o][%1) 22181 -9333 oo x1 F(tyt © ~7.95 0 |/¥a| +|-9980 sees -9399 |fx2} = | ecey2 ° o “s.60}{ Xs o = -9383 9383 } [xa] = [| Fre )3| Introduciendo un cambio de variable: @)]= (len ecuacion (b) aes} + Fels) = Feces] muitiplicando pe la traspuesta de ($l: (6) (6) (1) ©) fem) offeTfo Tin] = @)fere) — ecuacion ce» y normalizando la matriz (@Jde manera tal que: Tee) =(t) (matriz identidad) se obtiene: (Wl) = (ut) Entonces: (on) + (} Ga) =0 ( vibraciones libres) ecuacién (d) Nétese que esta es la ecuacién para sistemas con un grado de libertad.8s la ecuacién gue describe e1 movimiento de las magas de cada unc de los modos de vibracién del sistema PAGINA &4) Obtencién de las frecuencias naturales y sus correspondientes periodos de vibracién Es oportuno mostrar como se obtiene el determinante que permite calcular las frecuencias y periodos fundamentales de los tres modos de vibracion Para un sismtema de tres grados de libertad vibrando en vibra- cién libre se tienen tres ecuaciones diferenciales de la for- ma: Mi M1 + Kl. X1 + K12 X2 + K13 X3 =0 M2 H2 + K21 Xl + K2Z XZ + K23 xB = ecuacién (e) mg Xo + K81 X1 + K82 X2 + KS3 X38 =0 Como cada ecuacién describe un movimiento armnica ( cada modo de vibracién lo es ) se tiene una solucién del tipo: = Am cos (wt +8) Eouacién (£) Am= amplitud del movimiento $= Angulo de fase. Por lo tanto: ‘Am = -w* Am cos (wt +) Sustituyendo en la ecuacién (e) y cancelando el cos(wt+$), se tiene: Am1(K11-M1 w* ) + K12 Am2 * K13 Am3 =O Ami K21 + (K22-M2w* )AmzZ + K23 Am3 9 =0 Ecuacién (g) Ami K31 + =K92 Am2 + (K33-M3w* )Am3 =0 La ecuacién (g) puede escribirse en términos matriciales asi (K11 - M1 wy K12 K13 Am3 0 K21 (K2z - M2 we) 23. am2} =| 0 KS1 K82_— (K33-M3 w*)J | AmS, ° PAGINA 4Cuya solucién implica un problema de los valores propios de una ecuacién. Como se recordaré dichos valores se obtienen haciendo el determinante de la matriz igual a cero. 0 sea: ( Kii~ Mi we) kz K13 21 (K22 = Mawe) — K23 =o 31 Kaz (K33-Maw* ) Por lo que sustituyendo los valores numéricos se tiene: (22181-7.85 we) -9333 ° -9333 (18666-7.85 w*) -9333 =0 ° -9333 (9333-5.6 we) resolviendo el determinante se obtiene: w~ 6946 w* +12.11362174 E 6 w*-3254,35424 E 6 =0 Cuyas raices son: w, #1 =928.91 W1=18,14 seg-1 TL = 0,346 seg Wp" 2 =2842.67 w 2 48.04 seg-1 TZ = 0.181 see Wy $4192.07 w 3 =64.75 seg-1 TS = 0.096 seg Recuérdese: T=(2 n)/w PERIOD DE VIBRACION. ‘T1S0.347 SEG T2=0.191 SEG T3=0.096 SEG PRIMER MORO SEGUNDO _MODO TERCER MODO PAGINA 55) Obtencién de la matriz modal. La matriz modal contiene los valores numéricos que representan le relacion que existe entre los desplazamientos de las diferen- tes masas de la estructura vibrando segtin cada uno de los modos. Nétese que no es posible, de la solucién de la ecuacién (g), en ~ contrar los verdaderos valores de las smplitudes. PRIMER MODO. Para encontrar los valores % del primer modo se tiene: (w)*=927.10. Sustituyende en (g) los valores numéricos se encuentra: (22131-7.85*327.10) Am1 9333 Am2 -0 -9333 Am1 (18666-7.85*327.1)Am2 -9333Am3 =0 0 9333 AmZ + = (9333-5. 6*327.1)Am3=0 Efectuando: 19563 Am1 -9333 Am2 -0 Am3 — =0 -93339Am1 +16096Am2 -9333AmS =O OAm1 © -9333 Amz +7501 Am3 =0 Si se supone que Ami=1, entonces: Am2=2.10 y Am3=2.61. Graficamente: aye eel / 2.81 p10 a1 = [2.10 00 1,00 a SEGUNDO MODO Igualmente, al sustituir en ls ecuacién (g) el valor de la fre- cuencia natural del segundo modo w*=2345 y efectuar las operacio— nes, se tiene: PAGINA 6$723 Am1 - 9333 Am2- © Am3 =0 ~8333 Am2 + 266 Am2 -939§" Ans =0 0 Ami ~8333Am2 -3789 Am3 =0 81 se supone que Ami = 1 , se obtier Am2=0.40 y Am3=-0.98 Gréficamente: -0.98 @= | 0.40 1.00 TERCER MODO Sustituyendo en la ecuacién (g) el valor de la frecuencia na- tural del tecer modo (w3)"=4192 se tiene: -10776 Ami -9333 Am2 = +0 Am3_— =O -9333 Ami + 258 Am2 -9333Am3 =0 0 Am1 -9333 Am2 -14142Am3 =0 Suponiendo que Ami = 1 se tiene: Am2=-1,15 y Am3=0.76. Gréficamente: 0.76 @ | -1.15 1.00 PAGINA 7Por consiguiente la matriz modal es: 2.61 -0.98 0.76 a | 2.10 0.40 -2.16] , 1.00 1.00 1.00 6) Eouaciones del cédigo colombiano de construcciones sismo-re - sistentes. i) CORFICIENTE SISMICO MODAL En el articulo A.5.2.1 se define el coeficiente sismico modal como: Cem= Sam / R (A.5-1) Csm=Valor del espectro eldstico de aceleraciones para el pe- riodo de vibracién Tm. Para Cali se obtiene de la tabla A.9.1 (zona de riesgo aismico alto,pérticos en conereto reforzado) el valor R=6 Para Sam se tienen dos expresiones: sam-1.2 avs 1 / 1° (4.2.1) Sam = 2.5 Aa I (4.2.2) Pare Cali de la tabla A.2-1: Av= 0.25 Aaz 0.25 I=coeficiente de importancia. Se supone igual a 1,1 S=Coeficiente de sitio. Se supone igual a 1.2. Nétese: Sam=1.2%0.25¥1.21.1/T"= 0.396/1/* Sam=2.540.2541.1=0.6875 a -396/1 bor lo que T=0.437 seg. igualande: 0.687: PAGINA 8Representacién grdfica del coeficiente siemico de disefio. Ubicacién del proyecto : Cali. (25AaT)/R (12 AyST/72/3)/R T (39) De acuerdo con lo anterior se tiene para los modos de vibra- T1=0.346 SEG ( 0.346 SEG < 0.437 SEG ) Entonces: Sam=2.5%0.25*1.i1=0.6875 Cem=0.6875/6 =0.114 SEGUNDO MODO Tr 131 SEG ( 0.181 < 0.437 SEG ) Entonces: Sam=2.5%0.25%1.1 = 0.6875 Cam =0.8875/6 =0.114 TERCER MODO TS=0.096 SEG ( 0.096 < 0.487 SEG ). Entonces: Sam=2.5%*0.25*1.1=0. 6875 Cem=0.6875/6 = 0.114 ii)CORTANTE MODAL EN LA BASE De acuerdo con el cédigo (seccién A.5.3): "La parte del cortante en la base contribuida por el modo m, Vm, debe determinarse de acuerdo con la siguiente formula: ‘am Win PAGINA 9Donde: Win = Peso efectivo modal determinado de acuerdo con la f6rmu- la siguiente: (2 Wi &im > tin = ——_______ E Wi Sim? Sim= Amplitud de desplazamiento del nivel i del edificio cuando esta vibrando en el modo m CORTANTE MODAL EN LA BASE: PRIMER MODO DE VIBRACION. los cdélculos para obtener el cortante modal en la base debido al primer modo de vibracién se pueden organizar como sigue 791.12 For consiguiente el cortante modal en la base para el pri- mer modo de vibracién es: Vm = 0.114% 184.69 = 21.05 ton Notese que la participacién del primer modo sustancialmen- te més alta que la de los otros dos modos. (88%). PAGINA 10CORTANTE MODAL EN LA BASE: SEGUNDO MODO : pep] [arf | ote [oe | ee | 17 1.00 77.00 77.00 209 53.90 142.12 ( 53.90 )# fim = =20.44 ton 142.12 Por consiguiente el cortante modal en la base para el segundo modo de vibracién es: Vm= 0.114*20.44 = 2.33 ton CORTANTE MODAL EN LA BASE: TERCER MODO. 3 2 ; : ( 30.25 )* im = 4.34 ton 210.68 El cortante modal en la base es, para el tercer modo: Vm=0.114%4.34=0.50 ton Nétese que la suma de los cortantes modales en la base debe ser igual al peso del edifici W = Wm=1 + Wm=2 +Wm=3 = 184.69 + 20.44 + 4.35=209.48 ton . PAGINA 17ii) DISTRIBUCION EN LA ALTURA DEL CORTANTE MODAL. segin el cédigo “ la fuerza cortante modal Fxm , en cualquier nivel debe determinarse de acuerdo con la siguiente formula: Fxm = Cvx Vm (A.5-6) Donde: Wx ® xm Cvx = ( A.5-7) Wi ® im PRIMER MODO: FUERZAS MODALBS POR PISO (NIVELES): El cdlculo de las fuerzae Fam por piso se puede llevar a cabo como se indica en el cuadro siguiente: PISO We fxm Wi @im 3 55 2.61 143.55 0.376 2 77 2.10 161.70 0.423 7 1.00 TT.00 0.201 4.23 z 209 382.25 1.00 21.04) “Recuérdese que Vm=21.05 ton para el primer modo. Recuérdese que el cortante en la base para el segundo mo- do es Vm=2.33 ton. PAGINA 12El cortante en la base para el tercer modo es 0.50 ton 7) Representacién gréfica de los resultados PRIMER MODO = j= 7,91 TON == 8,90 TON’ 4,23 7T0N! ——m 21,05T CORTANTE EN LA BASE SEGUNDO_MODQ 0,627 3,50 TON CORTANTE 233 TON EN LA BASE PAGINA I3 FUERZAS DE PT FUERZAS DE PISOTERCER MODO —s 050 7) Valores de disetio del cédigo Las fuerzas cortantes de piso méximas se obtienen sumando los cortantes de piso de cada uno de los modos de vibracién Ee evidente que las fuerzas asi calculadas son muy desfavora— bles pues es poco probable que los tree modos de vibracién o- curran simultaneamente. Por esta razon el cédigo define los va- lores de disefio como sigue: “los valores de disefio del cortante en la base, los cortantes en cada nivel y los desplazamientos deben determinarse combinando ios valores [Link] combinacién debe hacerse tomando la raiz cusdrada de la suma de los cuadradoe de cada valor modal". De acuerdo con esto y para los valores indicados, ee obtiene: RUERZAS MODALES CORTANTES MODALES PRIMER MODO: F3= 7.91 TON V3= 7.91 TON F2= 8.90 TON VZ=7.91 +8.90 =16.81 TON Fl= 4.23 TON V1=7.9148.9+4.23=21.04 TON PAGINA 14UNDO MODO F3= -2.45 TON V3=-2.45 TON F2= 1.28 ‘TON V2e-2.4541.28 = -1.17 TON Fi= 3.50 TON Vis ~2,454+1.20+3.5=2.33 TON TERCER MODO F3= 0.69 TON Va= 0.89 TON F2= -1.46 TON V2=0.69-1.46= -0.77 TON Fi= 1.27 TON Viz 0,69-1.46+1,27=0.50 TON Resumiendo: | e.3: | 31 ras ee en te De acuerdo con esto, las fuerzas por piso son: F (tercer piso) = 8.31 ton FUERZAS DE DISENO F (segundo piso)= 16.87-8.31 =8.56 ton POR PISO F (primer piso j= 21.04-16.87=4.17 ton Graficamente: 8,31T 8,56T Z| 477 wae PAGINA 15Como se supuso que los pérticos en el sentido X son iguales basta dividir por tree cada fuerza de piso pare obtener la fuerza de portico (existen tres porticos en ese sentido) Si loe porticos fueran diferentes, seria necesario dietribuir las fuerzas por portico de acuerdo con la rigidez al desplaza miento de cada uno de ellos. Es decir PORTICO EJE UNO PORTICO EJE DOS | PORTICO EJE TRES 8.91/3=2.77T 8.31/3 =2.77.T 8.56/38 =2.85 T 4.17/3 =1.39 T Graficamente para el pértico del eje 1 se tiene: —~— 2,77T —m 2,85 T —~ 139T 8) Desplazamientos ‘fodales eldsticos Segin el cédigo, e] desplazamiento modal elastico del nivel x , Sxem, se obtiene de acuerdo con la ecuacién: g Tm? Fxm bxem AIS. 4m? We PAGINA 16PRIMER MODO, DESPLAZAMIENTOS MODALES ELASTICOS Tm = 0.346 seg PISO | Fxm | Wx | g Tm Fxm | 4 1 Wx éxem 7.91] 55 9.29 2171 0.0042 m 10.45 3040 0.0034 m SEGUNDO MODO: DESPLAZAMIENTOS MODALES ELASTICOS. ‘Tm=0.131 seg ~0.00019m 0.00007m 0.00019 P1so | Fxm Wx g Tm Fxm | 4 n* Wx 3 0.69 | 55 0.06 2171 0.000027m 2 1.46 | 77 <0. 3040 | -0.000043m| 1 1.27 | 77 0.41 3040 0.000362m| VALORES DE DISENO: DESPLAZAMIENTOS MODALES. TERCER PISO bxem = \[( 0.0042 )* + (-0.00018)* + 0.000027)" =0.0042 m SEGUNDO PISO xem = Y'(0.0084)* +(0.00007)* + (-0.000043)*=. 0.0034 m PAGINA 17PRIMER PISO Sxem = (0,0016)* +(0.00019)* + (0.000362)* = 0.00165 m DESPLAZAMIENTOS ELASTICOS: ——~ 042em a / —-/ 034M / / / 0,165 cm LIMITES DE LA DERIVA: De acuerdo con la definicién del cédigo: 8x= Cd 6xe — (A.8.1) Sx= desplazamiento horizontal real de la estructura tomando en cuenta los efectos ineldsticos. Cd= 5 ( tabla A.9-1 )- Zona de riesgo sismico alto. De acuerdo con la definicién del codigo: La deriva de disefio del piso Ax , es la diferencia en- tre los desplazamientos horizontales reales (en el ran- go ineldstico ) 5x, de dos niveles consecutivos y no de- be exceder e1 0.015 % de la altura del piso.” PAGINA 18LIMITES DE LA DERIVA: x COMENTARIOS: A) QUE SUCEDE SI SE TOMA LA INERCIA REAL DE LAS VIGAS. Si se quiere efectuar un cdélculo més detallado de las rigideces al desplazamiento de los pérticos del edificio es necesario en- tonces considerar la inercia de las vigas y aplicar una carga su- puesta en cada uno de los pisos. Al dividir la carga por el des - plazamiento que produce, se obtiene la rigidez al desplazamiento. En este ejemplo se supuso que todos los pérticos en el sentido X son iguales por lo que la rigidez al desplazamiente de un pér~ tico se multiplica por tres ( hay trea pérticos en la direccién X ) para obtener la matriz de rigidez del edificio en e] sentido X. Si los pérticos fueran diferentes, se sumaria la rigidez de cada uno de ellos. PORTICO DEL EJE 1 Para efectos de célculo se supone una carga de 6000 kg en el primer piso y se calculan los desplazamientos mediante un mé~ todo eldstico conocido (Cross, Kani, matricial o cualquier o- tro). PAGINA 19A continuacién se muestran los resultados obtenidos para los desplazamientos. VIGAS:30/35 PRIMER PISO ft <— a= 20647 | ——m 78814 17M / {Babbar daa 3 O TIO AOD Jock CRI SA OEDD DISEAO FSTRUCTURAL. ING. CARLOS RANTRD VALLECILLA E BEA TODA II TOOT ANAL ISES MODAL nas DE 2993 ANINY DISPLACEMENTS aOINT x Disp. ¥ pis) Zz Rate 1 9.000907 0, 000:54 +0 00005 2 0.009587 0, 000042 9.000004 3 9.090000 ~0,00019% -0,000058 4 [Link] 6.000255 9.006200 5 0.004195 9.000025 9.090107 6 0.000000 -0. 000% 9. 0002204 ? CRaGEy 9.0000% 9.000118 8 0.0) -9,900008 =, 000089 y 0.081 0.900099 9.900198 20 9-000000 9.000000 9.000000 uu 9-06 9.000009 -@.000000 a 0.00 9.000000 9.000000 PAGINA 20Nétese que al aplicar una fuerza de 6000 kg en e} primer piso se produce un desplazamiento de 0.086 om por lo que la rigidez al desplazamiento del piso es: z K1li= §000/0.086 =58139 kg/ cm =5814 ton/m ( contra 7377 ton/m que se obtuvieron al suponer las vigas con inercia infinita )}. Para el cAlculo de K21 y K31 se aplica al pértico en el pri- mer piso une carga de 5814 ton (que debe producir un despla ~ miento de 1 metro).Efectuando los célculos correspondientes se obtiene: : K21=-2984 TON/M K31= 464 TON/M SEGUNDO PISO Siguiendo con el procedimiento indicado en el pértico del'eje 1 se tiene: K82=-2192 ton/m K22=4545 ton/m KL 2964 ton/m TERCER PISO Iguaimente: Ka: 563 ton/m K23=-2132 ton/m K13=464 ton/m Multiplicande por tres las rigideces del pértico se obtiene ja matriz de rigidez del edificio (son tres pérticos parale~ al eje X) PAGINA 21174420 = «BBS 1392 K = | -8892 13635-6396 1392 -6396 4689 Al suponer las vigas de entrepiso con inercia infinita se obtuvo: 22131 © -9333 ° K =] -9339 18666 = -9333 ° -9333 9333 Valores mayores como era de esperarse ( con esa suposicién la es- tructura es més rigida). Sustituyendo en la ecuacién de valores propios se obtiene: (17442-7.85 w*) 8892 1392 -BB92, (13835-7.85 wi) -6398 =O 1392 -8396 —-( 4689-5.6 w*) Resolviendo el determinante se obtiene: w°-4796.12 w*+4916343.14 w*-471.6977516 E 6 = 0 Cuyas raices son: (wi)*=106.83 win10.33 Ti=0.608 seg (w2)?=1304.46 w2=96.11 T2=0.174 seg (w3 3384.83 wa=58.17 T3=0.108 seg Comparnacién: T1=0.347 seg 220.131 seg Vigas con inercia infinita T3=0.096 seg PAGINA 22B) CUALES SON LAS FUERZAS QUE SE OBTIENEN CON EL METODO DE LA FUERZA HORIZONTAL EQUIVALENTE . Como se sabe, el métotio de la fuerza horizontal equivalente se basa en el primer modo de vibracié[Link] parametros de cal culo son: i) Perfodo de vibracién de la estructura De acuerdo con el cédigo: = 0.08” * 3/4 0.08 (8/2 0.38 seg ( contra 0.546 ceg para Es decir: el primer modo de vibracién ). Ai) Coeficiente K K = 1 ( Periodo menor que 0.5 seg ) 411) Coeficiente sismico de disefio Ce=0.114 iv) Peso del edificio W= 55 + 2x77 =209 ton v) Cortante en la base. V=Cst W = 0.114209 = 28.83 ton vi) Fuerza de piso Cvx | Fpiso | Frértico 0.402] 9.58 | 3.19 0.387 22 | 3.07 0.211] 5.03 | 1.68 z | 209 1094.5] 1.00] 23 L 4 Graficamente: & PAGINA 23PORTICO EJE 1 \34om -277QK,_Ese —_ Tee o4zem lopoem 2850K J o34em 3 FUERZA HORIZONTAL EQUIVALENTE ANALISIS MODAL BIBLIOGRAFIA ~CODIGO COLOMBIANO DE CONSTRUCCIONES [Link]- CRETO 1400 -DYNAMICS OF STRUCTURES. RAY W. CLOUGHJ. PENZIEN. McGRAW- HILL. -CONFERENCIAS DE ANALISIS DINAMICO. LUIS ENRIQUE GARCIA. U- NIANDES. ~CURSO DE ANALISIS DINAMICO. LUIS GUILLERMO AYCARDY. UNIVER- SIDAD NACIONAL. MAGISTER DE ESTRUCTURAS. -NOTAS DE DINAMICO. DIVISION DE EDUCACION CONTINUA. U.N.A.M. MEXICO. -VIBRATION PROBLEMS IN ENGINEERING. [Link] WILEY SONS. PAGINA 24