INGENIERA ANTISSMICA

1) RESPUESTA A UNA FUERZA EXCITADORA

Para el sistema de un grado de libertad mostrado en la Tabla 1 calcular la respuesta
dinmica en desplazamiento, velocidad y aceleracin para la fuerza excitadora
mostrada en la Figura 1.

La respuesta ser calculada para un valor de fuerza Po y una duracin de la excitacin (tr)
relacionada con el periodo (T) del sistema de un grado de libertad, mostrado en la Tabla 2.
El mtodo a emplear es mostrado en la Tabla 3.

SOLUCIN
MTODO NUMRICO (aceleracin lineal)
La constante de rigidez del sistema durante el comportamiento elstico es:
= 423 /
423 .
El coeficiente de amortiguacin = = 0.05 2 31 = 3.66
9.81

El desplazamiento inicial y la velocidad inicial son 0 = 0 = 0. La aceleracin inicial es:

(0)
= =0

El periodo natural es T = 2**/ = 0.543 seg.
Para un tr = 10.86 seg
Seleccionamos un intervalo = 0.05

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

La rigidez efectiva para la ecuacin es

= + 62 + 3

= 423 + 63.16
2 +
33.66
= 8226.6
0.05 0.05

El incremento efectivo de la excitacin dado por la ecuacin es

6
= + ( + 3 ) + (3 + )
t 2

= + 390.18 + 9.57
El incremento de velocidad dado por la ecuacin es
3
= 3
2
= 60 3 0.025
Este mtodo supone aproximaciones significativas: la aceleracin vara linealmente durante
el intervalo de tiempo y las propiedades de amortiguacin y rigidez del sistema se
evalan al inicio de cada incremento de tiempo y permanecen constantes. Estas no
suposiciones no introducen errores, salvo que sea pequeo. Generalmente el intervalo
de tiempo debe ser no mayor a un dcimo del periodo natural de la estructura.
RESPUESTA DINMICA DE DESPLAZAMIENTO, VELOCIDAD Y ACELERACIN
USANDO MTODOS NUMRICOS
F Y ( m) K

(m/s) ( m/s2) ( tnf/m) (Tnf) (Tnf)
0 0 0 0 423.00 8226.6 0.04 0.04144 0.00001
0.04 5.04E-06 3.02E-04 1.21E-02 423.00 8226.6 0.04 0.27504 0.00003
0.08 3.85E-05 1.10E-03 1.98E-02 423.00 8226.6 0.04 0.65989 0.00008
0.12 1.19E-04 2.12E-03 2.10E-02 423.00 8226.6 0.04 1.06921 0.00013
0.17 2.49E-04 3.04E-03 1.57E-02 423.00 8226.6 0.04 1.37543 0.00017
0.21 4.16E-04 3.57E-03 5.76E-03 423.00 8226.6 0.04 1.48962 0.00018
0.25 5.97E-04 3.58E-03 -5.37E-03 423.00 8226.6 0.04 1.38681 0.00017
0.29 7.65E-04 3.09E-03 -1.43E-02 423.00 8226.6 0.04 1.11031 0.00013
0.33 9.00E-04 2.28E-03 -1.83E-02 423.00 8226.6 0.04 0.75468 0.00009
0.37 9.92E-04 1.41E-03 -1.64E-02 423.00 8226.6 0.04 0.43383 0.00005
0.41 1.04E-03 7.57E-04 -9.63E-03 423.00 8226.6 0.04 0.24488 0.00003
0.46 1.07E-03 5.12E-04 -2.12E-04 423.00 8226.6 0.04 0.23905 0.00003
0.50 1.10E-03 7.25E-04 8.76E-03 423.00 8226.6 0.04 0.40836 0.00005

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

0.54 1.15E-03 1.31E-03 1.46E-02 423.00 8226.6 0.04 0.69123 0.00008

0.58 1.24E-03 2.06E-03 1.56E-02 423.00 8226.6 0.04 0.99435 0.00012
0.62 1.36E-03 2.74E-03 1.17E-02 423.00 8226.6 0.04 1.22313 0.00015
0.66 1.51E-03 3.14E-03 4.44E-03 423.00 8226.6 0.04 1.31100 0.00016
0.70 1.67E-03 3.16E-03 -3.80E-03 423.00 8226.6 0.04 1.23837 0.00015
0.75 1.82E-03 2.81E-03 -1.04E-02 423.00 8226.6 0.04 1.03622 0.00013

Respuesta de un sistema de un grado de libertad

0.025

0.02

0.015

0.01
y (m), v(m/s) , a( m/s2)

0.005

0
0 2 4 6 8 10 12
-0.005

-0.01

-0.015
Y ( cm) y ( m/s) ( m/s2)
-0.02

-0.025
t (s)

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

Para un tr = 0.543
Seleccionamos un intervalo = 0.05

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

Para un tr = 0.181
Seleccionamos un intervalo = 0.01

MTODO ANALITICO
Se emple la integral de Duhamel
Para un tr = 10.86 seg
Seleccionamos un intervalo = 0.05

Tr 10.86 C1 1.4629E-06 C12 0.00195625

B 0.05 C2 -0.00016843 C13 9.7935E-05
Wn 11.5712 C3 0.00195625 C14 -0.00193669
C1 1.4629E-06 C4 -1.6906E-05 C15 0.00029283
Wd 11.5567269 C5 -1.6906E-05 C16 0.00205741
( ) 0.33473167 C6 0.00016843 C17 0.02245804
133.557938 C7 -1.6906E-05 C18 0
C8 -0.00016843 C19 -0.02266456
C9 0.02124492 C20 0.00450461
C10 -0.00194732 C21 0.0222124
C11 8.0539E-05

RESPUESTA DINMICA DE DESPLAZAMIENTO, VELOCIDAD Y ACELERACIN

USANDO MTODOS ANALTICOS

t Y (m) Y (m/s) Y (m/s2)

0 -1.54433E-05 8.9349E-06 0.00205741
0.05 8.06774E-05 6.53452E-05 0.00221868
0.1 0.000176894 0.000120309 0.00236885
0.15 0.000273201 0.000173857 0.00250843
0.2 0.000369596 0.000226019 0.00263787
0.25 0.000466075 0.000276828 0.00275763

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

0.3 0.000562635 0.000326311 0.00286814

0.35 0.000659272 0.000374499 0.00296983
0.4 0.000755982 0.000421421 0.00306309
0.45 0.000852764 0.000467104 0.00314831
0.5 0.000949613 0.000511578 0.00322588
0.55 0.001046527 0.000554868 0.00329614
0.6 0.001143502 0.000597003 0.00335944
0.65 0.001240537 0.000638009 0.00341613
0.7 0.001337629 0.000677911 0.00346651
0.75 0.001434775 0.000716736 0.0035109
0.8 0.001531972 0.000754509 0.0035496

Para un tr = 0.543
Seleccionamos un intervalo = 0.05
Tr 0.53 C1 2.9258E-05 C12 0.03912508
B 0.05 C2 -0.00336855 C13 0.0019587
Wn 11.5712 C3 0.03912508 C14 -0.03873383
C1 2.9258E-05 C4 -0.00033812 C15 0.00585652
Wd 11.5567269 C5 -0.00033812 C16 0.04114823
( ) 0.33473167 C6 0.00336855 C17 0.44916078
133.557938 C7 -0.00033812 C18 0
C8 -0.00336855 C19 -0.45329111
C9 0.02124492 C20 0.0900921
C10 -0.03894638 C21 0.44424795
C11 0.00161079

RESPUESTA DINMICA DE DESPLAZAMIENTO, VELOCIDAD Y ACELERACIN

USANDO MTODOS ANALTICOS
t Y (m) Y (m/s) Y (m/s2)
0 -0.000308867 0.000178698 0.04114823
0.05 0.001613548 0.001306905 0.04437353
0.1 0.00353787 0.002406175 0.04737706
0.15 0.00546402 0.003477131 0.05016855

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

0.2 0.007391923 0.00452039 0.05275737

0.25 0.009321506 0.005536556 0.05515256
0.3 0.0112527 0.006526227 0.05736281
0.35 0.013185436 0.007489989 0.05939651
0.4 0.015119649 0.00842842 0.06126173
0.45 0.017055277 0.009342088 0.06296624
0.5 0.018992258 0.010231552 0.06451752
0.55 0.020397152 0.011437549 0.1002751
0.6 0.020430648 0.011201478 0.09606349
0.65 0.020463003 0.010968464 0.09200103
0.7 0.020494252 0.010738547 0.08808305
0.75 0.020524429 0.010511763 0.084305
0.8 0.020553566 0.010288143 0.08066246
0.85 0.020581696 0.010067712 0.07715113
0.9 0.02060885 0.009850492 0.07376683
0.95 0.020635057 0.009636501 0.07050552
1 0.020660346 0.009425753 0.06736326
1.05 0.020684748 0.009218256 0.0643362
1.1 0.020708288 0.009014019 0.06142065

Para un tr = 0.181
Seleccionamos un intervalo = 0.01

Tr 0.181 C1 8.7773E-05 C12 0.11737524

B 0.05 C2 -0.01010566 C13 0.00587611
Wn 11.5712 C3 0.11737524 C14 -0.11620149
C1 8.77735E-05 C4 -0.00101437 C15 0.01756957
Wd 11.55672695 C5 -0.00101437 C16 0.12344469
( ) 0.334731674 C6 0.01010566 C17 1.34748233
133.5579378 C7 -0.00101437 C18 0
C8 -0.01010566 C19 -1.35987332
C9 0.02124492 C20 0.27027631
C10 -0.11683915 C21 1.33274384
C11 0.00483236

RESPUESTA DINMICA DE DESPLAZAMIENTO, VELOCIDAD Y ACELERACIN

USANDO MTODOS ANALTICOS
t Y (m) Y (m/s) Y (m/s2)
0 -0.000926601 0.000536094 0.12344469
0.005 -0.00035014 0.000878517 0.12444307

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

0.01 0.00022638 0.001220054 0.12543453

0.015 0.000802958 0.001560708 0.1264191
0.02 0.001379596 0.001900481 0.1273968
0.025 0.001956292 0.002239374 0.12836767
0.03 0.002533047 0.00257739 0.12933174
0.035 0.00310986 0.00291453 0.13028903
0.04 0.003686731 0.003250796 0.13123959
0.045 0.004263659 0.00358619 0.13218343
0.05 0.004840645 0.003920715 0.13312059
0.055 0.005417688 0.004254371 0.1340511
0.06 0.005994788 0.004587161 0.13497498
0.065 0.006571945 0.004919087 0.13589227
0.07 0.007149158 0.00525015 0.136803
0.075 0.007726428 0.005580353 0.1377072
0.08 0.008303753 0.005909697 0.13860489
0.085 0.008881134 0.006238184 0.1394961
0.09 0.009458571 0.006565816 0.14038087
0.095 0.010036063 0.006892596 0.14125922
0.1 0.01061361 0.007218524 0.14213118

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

2) RESPUESTA A UN SISMO
Los mtodos para hallar la respuesta dinmica de un sistema de un grado de libertad
resuelven la ecuacin:

Donde es la respuesta en desplazamiento del sistema para cada instante de tiempo.

Para el caso de un sismo, la fuerza puede ser considerada como una fuerza inercial igual
a siendo la aceleracin del suelo o la base. Entonces

Se solicita hallar la respuesta dinmica en desplazamiento, velocidad y aceleracin de un

sistema de un grado de libertad de periodo T = 2.0 s, masa 1.0 kg y amortiguamiento =5%,
frente al registro del sismo del 03 de octubre de 1974 obtenido en la estacin Parque de la
Reserva (PQR). Esta estacin, ahora deshabilitada, contaba con un acelergrafo analgico
que usaba una frecuencia de muestreo de 50 muestras por segundo. El mtodo y la
direccin del registro a utilizar se muestra en la Tabla 3.

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

Utilizando el programa DEGTRA (Ordaz y Montoya, 2013) calcular el espectro de respuesta

de aceleraciones relativas (AR) para la direccin del registro utilizado, luego, variando el
periodo T (0.05, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, 0.90, 1.00, 1.10, 1.20, 1.30,
1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 2.50, 3.00, 3.50, 4.00, 4.50, 5.00) del sistema de un
grado de libertad calcular la respuesta mxima de aceleracin y esbozar el espectro de
respuesta de aceleraciones relativas y compararlo con el obtenido con el programa
DEGTRA. La Tabla 4 muestra los parmetros a utilizar en el programa DEGTRA.

MTODO NUMRICO

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

Respuesta de un sistema de un grado de libertad a una fuerza

ssmica
0.0150
Desplazamiento ( cm)

0.0100
0.0050
0.0000
0.0000 10.0000 20.0000 30.0000 40.0000 50.0000 60.0000 70.0000
-0.0050
-0.0100
-0.0150
t(s)

Respuesta de la velocidad a una fuerza ssmica

1.0000
Velocidad cm/s

0.5000

0.0000
0.0000 20.0000 40.0000 60.0000 80.0000 100.0000 120.0000
-0.5000

-1.0000
t(s)

Aceleracion
150.0000

100.0000
Aceleracion (cm/s2)

50.0000

0.0000
0.0000 10.0000 20.0000 30.0000 40.0000 50.0000 60.0000 70.0000
-50.0000

-100.0000

-150.0000
t(s)

PROGRAMA DEGTRA

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

Acelermetro del sismo del 03 de octubre de 1974 obtenido en la estacin Parque

de la Reserva (PQR).

ESPECTRO DE RESPUESTA (Degtra)

600.00
Aceleracion Relativa (cm/s2)

500.00

400.00

300.00

200.00

100.00

0.00
0.00 0.40 0.80 1.20 1.60 2.00 2.40 2.80 3.20 3.60 4.00 4.40 4.80 5.20 5.60
T (s)

Del espectro de aceleraciones relativas obtenidas se puede observar que el pico ms alto
es cuando el periodo es 0.39 y la aceleracin 518.81 cm/s2.

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

ESPECTRO DE RESPUESTA USANDO MTODO NUMRICOS (aceleracin lineal)

Se vari el periodo (0.05, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, 0.90, 1.00, 1.10,
1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 2.50, 3.00, 3.50, 4.00, 4.50, 5.00); Se
obtuvo para cada periodo una aceleracin mxima como se describe a continuacin:

Masa (Kg) 1
T (s) 0.05 Wn= 125.663706
Amortiguamiento (%) 0.05

Por el mtodo numrico de Newmark (aceleracin lineal), se obtuvo el espectro de

respuesta en donde el pico esta cuando el periodo es 0.4 s para una aceleracin de 508.85
cm/s2.
Comparando el espectro de respuesta obtenido de del programa DEGTRA y por mtodos
numricos, se puede notar que Ambos presentan el pico de aceleraciones para un mismo
periodo (0.39 s -0.40s). Teniendo el espectro del mtodo numrico tiene una aceleracin

INGENIERA ANTISSMICA | UNMSM

 INGENIERA ANTISSMICA

de 508.85 cm/s2 aproximadamente, mientras que con el programa DEGTRA el pico de

aceleracin del espectro tiene un valor 518.81 cm/s2.

INGENIERA ANTISSMICA | UNMSM