Universidad Técnica de Oruro

Facultad Nacional de Ingeniería

Ingeniería Civil

DINAMICA ESTRUCTURAL CIV 3308

PROBLEMA 1

Determine el comportamiento dinámico del tanque de agua modelado en la figura, bajo las características idealizadas dadas, sujeto
a una carga dinámica P(t) cuya variación se muestra en la figura. Asumiendo que el coeficiente de amortiguación es el 20% de la
amortiguación crítica y asumiendo la variación del tiempo igual a intervalos de 𝜋/[Link] un cuadro de los valores de tiempo,
desplazamiento, velocidad y aceleración desde t=0 hasta t=𝜋 procesado por computador. Así mismo presente este cuadro en
forma gráfica. Determine el máximo desplazamiento en lo alto del tanque y determine la máxima fuerza cortante. Altura=14 m

SOLUCION

𝑊 22[𝑡𝑛] 𝑠2
𝑚= = = 0.022 [𝑡𝑛 ∗ ]
𝑔 980 [𝑐𝑚] 𝑐𝑚
𝑠2
𝑐 20 ∗ 𝐶𝑐𝑟𝑖𝑡
𝜉= = = 0.2
𝐶𝑐𝑟𝑖𝑡 100 ∗ 𝐶𝑐𝑟𝑖𝑡

16 𝑠
𝑐 = 2𝑚𝑤𝜉 = 2 ∗ 0.022 ∗ √ ∗ 0.2 = 0.24 [𝑡𝑛 ∗ ]
0.022 𝑐𝑚

𝑡𝑛
𝑘 = 16[ ]
𝑐𝑚

0.022𝑥̈ (𝑡) + 0.24𝑥̇ (𝑡) + 16𝑥(𝑡) = 𝑃(𝑡)

 Función complementaria
0.022𝑠 2 + 0.24𝑠 + 16 = 0
𝑠 = −5,339 ± 26,158𝑖

𝑥𝑐 (𝑡) = 𝑒 −5,339𝑡 (𝐴 cos 26,158 𝑡 + 𝐵 sin 26,158𝑡)

 Integral Particular
Para P(t)
2𝜋 2𝜋
𝑇= 𝑤= =2
𝑤 𝜋
Del grafico
𝑃(𝑡) ∗ 3 = 9 − 9 cos 2𝑡

𝑃(𝑡) = 3 − 3 cos 2𝑡
La integral particular podría ser:

𝑥𝑝 (𝑡) = 𝐶 + 𝐷 cos 2𝑡 + 𝐸 sin 2𝑡

𝑥𝑝̇ (𝑡) = −2𝐷 sin 2𝑡 + 2𝐸 cos 2𝑡
𝑥𝑝̈ (𝑡) = −4𝐷 cos 2𝑡 − 4𝐸 sin 2𝑡

Sustituyendo estos valores en la ecuación de movimiento tenemos:

0.022(−4𝐷 cos 2𝑡 − 4𝐸 sin 2𝑡) + 0.24(−2𝐷 sin 2𝑡 + 2𝐸 cos 2𝑡) + 16(𝐶 + 𝐷 cos 2𝑡 + 𝐸 sin 2𝑡) = 3 − 3 cos 2𝑡

(15,912𝐷 + 0.479𝐸) cos 2𝑡 + (−0.479𝐷 + 15,912𝐸) sin 2𝑡 + 16𝐶 = 3 − 3 cos 2𝑡

16𝐶 = 3
15,912𝐷 + 0.479𝐸 = −3
−0.479𝐷 + 15,912𝐸 = 0

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

Luego 𝐶 = 0.1875 𝐷 = −0.1884 𝐸 = −0.0057

La solución completa es

𝑥(𝑡) = 𝑒 −5,339𝑡 (𝐴 cos 26,158 𝑡 + 𝐵 sin 26,158𝑡) + 0.1875 − 0.1884 cos 2𝑡 − 0.0057 sin 2𝑡

Entonces

𝑥̇ (𝑡) = −5.339𝑒 −5,339𝑡 (𝐴 cos 26,158 𝑡 + 𝐵 sin 26,158𝑡)

+ 𝑒 −5,339𝑡 (−26.158𝐴 sin 26,158 𝑡 + 26.158𝐵 cos 26,158𝑡) + 0.3768 sin 2𝑡 − 0.0114 cos 2𝑡

El sistema esta inicialmente en reposo

𝑥(0) = 0 = 𝐴 + 0.1875 − 0.1884 𝑦 𝐴 = 0.0009

𝑥̇ (0) = 0 = −5.339𝐴 + 26,158𝐵 − 0.0114 cos 2𝑡 𝑦 𝐵 = 0.0006

 Ecuación del movimiento

Por lo tanto para 𝟎 < 𝒕 < 𝝅

𝑥(𝑡) = 𝑒 −5,339𝑡 (0.0009 cos 26,158 𝑡 + 0.0006 sin 26,158𝑡) + 0.1875 − 0.1884 cos 2𝑡 − 0.0057 sin 2𝑡

𝑥̇ (𝑡) = 𝑒 −5,339𝑡 (0.0109 cos 26,158 𝑡 − 0.0267 sin 26,158𝑡) − 0.0114 cos 2𝑡 + 0.3768 sin 2𝑡

𝑥̈ (𝑡) = 𝑒 −5,339𝑡 (−0.7578 cos 26,158 𝑡 − 0.1421 sin 26,158𝑡) + 0.7536 cos 2𝑡 + 0.0228 sin 2𝑡

 Desplazamiento Máximo
Para que x sea máximo 𝑥̇ = 0

𝑥̇ (𝑡) = 𝑒 −5,339𝑡 (0.0109 cos 26,158 𝑡 − 0.0267 sin 26,158𝑡) − 0.0114 cos 2𝑡 + 0.3768 sin 2𝑡 = 0
𝑡 = 1,5859 𝑠

Reemplazando ese valor en la ecuación del desplazamiento tenemos

𝑥𝑚𝑎𝑥 = 𝑒 −5,339∗1.5859 (0.0009 cos 26,158 ∗ 1.5859 + 0.0006 sin 26,158 ∗ 1.5859) + 0.1875
− 0.1884 cos 2 ∗ 1.5859 − 0.0057 sin 2 ∗ 1.5859

𝑥𝑚𝑎𝑥 = 0.376[𝑐𝑚]
 Fuerza Cortante Máxima
𝑡𝑛
𝑄𝑚𝑎𝑥 = 𝑘 ∗ 𝑥𝑚𝑎𝑥 = 16 [ ] ∗ 0.3760[𝑐𝑚]
𝑐𝑚

𝑄𝑚𝑎𝑥 = 6.0160[𝑡𝑛]

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

t[s] x(t)[cm] x*(t) [cm/s] x**(t) [cm/s2]

0 0.0000 0.0000 -0.0042
0.1963 0.0112 0.1438 0.6412
0.3927 0.0501 0.2599 0.6237
0.5890 0.1101 0.3430 0.3388
0.7854 0.1818 0.3764 0.0221
0.9817 0.2543 0.3525 -0.2711
1.1781 0.3167 0.2745 -0.5178
1.3744 0.3594 0.1547 -0.6873
1.5708 0.3759 0.0114 -0.7534
1.7671 0.3637 -0.1337 -0.7049
1.9635 0.3247 -0.2584 -0.5490
2.1598 0.2649 -0.3438 -0.3095
2.3562 0.1932 -0.3768 -0.0228
2.5525 0.1207 -0.3525 0.2673
2.7489 0.0583 -0.2745 0.5168
2.9452 0.0156 -0.1547 0.6875
3.1416 -0.0009 -0.0114 0.7536

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

PROBLEMA 2

Determine el comportamiento de la torre mostrada en la figura sujeta a una carga dincamica impulsiva de duración 0.80
seg. Asuma una amortiguación igual al 10% de la amortiguación critica, asumiendo la variación del tiempo igual a
intervalos de 0.08 seg. Encuentre el máximo: desplazamiento y velocidad y determine la máxima fuerza cortante.
Altura=17 m procesado por computador

SOLUCION

𝑊 16[𝑡𝑛] 𝑠2
𝑚= = = 0.0163 [𝑡𝑛 ∗ ]
𝑔 980 [𝑐𝑚] 𝑐𝑚
𝑠 2

𝑐 10 ∗ 𝐶𝑐𝑟𝑖𝑡
𝜉= = = 0.1
𝐶𝑐𝑟𝑖𝑡 100 ∗ 𝐶𝑐𝑟𝑖𝑡

12 𝑠
𝑐 = 2𝑚𝑤𝜉 = 2 ∗ 0.0163 ∗ √ ∗ 0.1 = 0.0884 [𝑡𝑛 ∗ ]
0.0163 𝑐𝑚

𝑡𝑛
= 12[ ]
𝑐𝑚

0.0163𝑥̈ (𝑡) + 0.0884𝑥̇ (𝑡) + 12𝑥(𝑡) = 𝑃(𝑡)

 Función complementaria
0.0163𝑠 2 + 0.0884𝑠 + 12 = 0
𝑠 = −2.7117 ± 26,9971𝑖

𝑥𝑐 (𝑡) = 𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡)

 Integral Particular
Para 0<t<0,16
Del grafico
𝑃(𝑡) 7
= 𝑡
2 0.16

𝑃(𝑡) = 87.5𝑡
La integral particular podría ser:

𝑥𝑝 (𝑡) = 𝐶 + 𝐷𝑡
𝑥𝑝̇ (𝑡) = 𝐷
𝑥𝑝̈ (𝑡) = 0

Sustituyendo estos valores en la ecuación de movimiento tenemos:

0.0163(0) + 0.0884(𝐷) + 12(𝐶 + 𝐷𝑡) = 87.5𝑡

0.0884𝐷 + 12𝐶 + 12𝐷𝑡 = 87.5𝑡

0.0884𝐷 + 12𝐶 = 0
12𝐷 = 87.5

Luego 𝐶 = −0.0537 𝐷 = 7.2917

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL
 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

La solución completa es

𝑥(𝑡) = 𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡) − 0.0537 + 7.2917𝑡

Entonces

𝑥̇ (𝑡) = −2.7117𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡)

+ 𝑒 −2.7117𝑡 (−26.9971𝐴 sin 26,9971 𝑡 + 26.9971𝐵 cos 26,9971𝑡) + 7.2917

Para 0.16<t<0,32
Del grafico
𝑃(𝑡)
=7
2

𝑃(𝑡) = 14
La integral particular podría ser:

𝑥𝑝 (𝑡) = 𝐶
𝑥𝑝̇ (𝑡) = 0
𝑥𝑝̈ (𝑡) = 0

Sustituyendo estos valores en la ecuación de movimiento tenemos:

0.0163(0) + 0.0884(0) + 12(𝐶) = 14

12𝐶 = 14

Luego 𝐶 = 1.1667

La solución completa es

𝑥(𝑡) = 𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡) + 1.1667

Entonces

𝑥̇ (𝑡) = −2.7117𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡)

+ 𝑒 −2.7117𝑡 (−26.9971𝐴 sin 26,9971 𝑡 + 26.9971𝐵 cos 26,9971𝑡)

Para 0.32<t<0,64
Del grafico
𝑃(𝑡) 14
=7− (𝑡 − 0.32)
2 0.32

𝑃(𝑡) = 42 − 87.5𝑡
La integral particular podría ser:

𝑥𝑝 (𝑡) = 𝐶 + 𝐷𝑡
𝑥𝑝̇ (𝑡) = 𝐷
𝑥𝑝̈ (𝑡) = 0

Sustituyendo estos valores en la ecuación de movimiento tenemos:

0.0163(0) + 0.0884(𝐷) + 12(𝐶 + 𝐷𝑡) = 42 − 87.5𝑡

0.0884𝐷 + 12𝐶 + 12𝐷𝑡 = 42 − 87.5𝑡

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

0.0884𝐷 + 12𝐶 = 42
12𝐷 = −87.5

Luego 𝐶 = 3.5537 𝐷 = −7.2917

La solución completa es

𝑥(𝑡) = 𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡) + 3.5537 − 7.2917𝑡

Entonces

𝑥̇ (𝑡) = −2.7117𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡)

+ 𝑒 −2.7117𝑡 (−26.9971𝐴 sin 26,9971 𝑡 + 26.9971𝐵 cos 26,9971𝑡) − 7.2917

El sistema esta inicialmente en reposo

𝑥(0) = 0 = 𝐴 + 0.1875 − 0.1884 𝑦 𝐴 = 0.0009

𝑥̇ (0) = 0 = −5.339𝐴 + 26,158𝐵 − 0.0114 cos 2𝑡 𝑦 𝐵 = 0.0006

Para 0.64<t<0,8
Del grafico
𝑃(𝑡) 7
= −7 + (𝑡 − 0.64)
2 0.16

𝑃(𝑡) = −70 + 87.5𝑡

La integral particular podría ser:

𝑥𝑝 (𝑡) = 𝐶 + 𝐷𝑡
𝑥𝑝̇ (𝑡) = 𝐷
𝑥𝑝̈ (𝑡) = 0

Sustituyendo estos valores en la ecuación de movimiento tenemos:

0.0163(0) + 0.0884(𝐷) + 12(𝐶 + 𝐷𝑡) = 42 − 87.5𝑡

0.0884𝐷 + 12𝐶 + 12𝐷𝑡 = −70 + 87.5𝑡

0.0884𝐷 + 12𝐶 = −70

12𝐷 = 87.5

Luego 𝐶 = −5.8870 𝐷 = 7.2917

La solución completa es

𝑥(𝑡) = 𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡) − 5.8870 + 7.2917𝑡

Entonces

𝑥̇ (𝑡) = −2.7117𝑒 −2.7117𝑡 (𝐴 cos 26,9971 𝑡 + 𝐵 sin 26,9971𝑡)

+ 𝑒 −2.7117𝑡 (−26.9971𝐴 sin 26,9971 𝑡 + 26.9971𝐵 cos 26,9971𝑡) + 7.2917

 Condiciones iniciales para cada tramo

Para 𝟎 < 𝒕 < 𝟎. 𝟏𝟔
El sistema esta inicialmente en reposo

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

𝑥(0) = 0 = 𝐴 − 0.0537 𝑦 𝐴 = 0.0537

𝑥̇ (0) = 0 = −2.7117𝐴 + 26,9971𝐵 + 7.2917 𝑦 𝐵 = −0.2647

Para las condiciones iniciales del otro tramo:

𝑥(0.16) = 𝑒 −2.7117∗0.16 (0.0537 cos 26,9971 ∗ 0.16 − 0.2647 sin 26,9971 ∗ 0.16) − 0.0537 + 7.2917 ∗ 0.16
𝑥(0.16) = 1.2581[𝑐𝑚]

𝑥̇ (0.16) = 𝑒 −2.7117∗0.16 (−7.2917 cos 26,9971 ∗ 0.16 + 0.7320 sin 26,9971 ∗ 0.16) + 7.2917
𝑐𝑚
𝑥̇ (0.16) = 9.5387[ ]
𝑠
Para 𝟎. 𝟏𝟔 < 𝒕 < 𝟎. 𝟑𝟐
𝑥(0.16) = 𝑒 −2.7117∗0.16 (𝐴 cos 26,9971 ∗ 0.16 + 𝐵 sin 26,9971 ∗ 0.16) + 1.1667 = 1.2581

𝑥̇ (0.16) = −2.7117𝑒 −2.7117∗0.16 (𝐴 cos 26,9971 ∗ 0.16 + 𝐵 sin 26,9971 ∗ 0.16)

+ 𝑒 −2.7117𝑡 (−26.9971𝐴 sin 26,9971 ∗ 0.16 + 26.9971𝐵 cos 26,9971 ∗ 0.16) = 9.5387
Resolviendo el sistema
𝐴 = 0.4628 𝑦 𝐵 = −0.3445

Para las condiciones iniciales del otro tramo

𝑥(0.32) = 𝑒 −2.7117∗0.32 (0.4628 cos 26,9971 ∗ 0.32 − 0.3445 sin 26,9971 ∗ 0.32) + 1.1667
𝑥(0.32) = 0.9270[𝑐𝑚]

𝑥̇ (0.32) = 𝑒 −2.7117∗0.32 (−10,5555 cos 26,9971 ∗ 0.32 − 11.5601 sin 26,9971 ∗ 0.32)
𝑐𝑚
𝑥̇ (0.32) = −0.3003[ ]
𝑠
Para 𝟎. 𝟑𝟐 < 𝒕 < 𝟎. 𝟔𝟒
𝑥(0.32) = 𝑒 −2.7117∗0.32 (𝐴 cos 26,9971 ∗ 0.32 + 𝐵 sin 26,9971 ∗ 0.32) + 3.5537 − 7.2917 ∗ 0.32 = 0.9270

𝑥̇ (0.32) = −2.7117𝑒 −2.7117∗0.32 (𝐴 cos 26,9971 ∗ 0.32 + 𝐵 sin 26,9971 ∗ 0.32)

+ 𝑒 −2.7117𝑡 (−26.9971𝐴 sin 26,9971 ∗ 0.32 + 26.9971𝐵 cos 26,9971 ∗ 0.32) − 7.2917 = −0.3003
Resolviendo el sistema
𝐴 = 0.1073 𝑦 𝐵 = −0.8805

Para las condiciones iniciales del otro tramo

𝑥(0.64) = 𝑒 −2.7117∗0.64 (0.1073 cos 26,9971 ∗ 0.64 − 0.8805 sin 26,9971 ∗ 0.64) + 3.5537 − 7.2917 ∗ 0.64
𝑥(0.64) = −0.9578[𝑐𝑚]

𝑥̇ (0.64) = 𝑒 −2.7117∗0.64 (−24.0619 cos 26,9971 ∗ 0.64 − 0.5091 sin 26,9971 ∗ 0.64) − 7.2917
𝑐𝑚
𝑥̇ (0.64) = −7.1993[ ]
𝑠
Para 𝟎. 𝟔𝟒 < 𝒕 < 𝟎. 𝟖
𝑥(0.64) = 𝑒 −2.7117∗0.64 (𝐴 cos 26,9971 ∗ 0.64 + 𝐵 sin 26,9971 ∗ 0.64) − 5.8870 + 7.2917 ∗ 0.64 = −0.9578

𝑥̇ (0.64) = −2.7117𝑒 −2.7117∗0.64 (𝐴 cos 26,9971 ∗ 0.64 + 𝐵 sin 26,9971 ∗ 0.64)

+ 𝑒 −2.7117𝑡 (−26.9971𝐴 sin 26,9971 ∗ 0.64 + 26.9971𝐵 cos 26,9971 ∗ 0.64) + 7.2917 = −7.1993
Resolviendo el sistema
𝐴 = −2.8956 𝑦 𝐵 = −1.4871

 Ecuación del movimiento

Para 𝟎 < 𝒕 < 𝟎. 𝟏𝟔
UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL
 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

𝑥(𝑡) = 𝑒 −2.7117∗𝑡 (0.0537 cos 26,9971 𝑡 − 0.2647 sin 26,9971𝑡) − 0.0537 + 7.2917𝑡

𝑥̇ (𝑡) = 𝑒 −2.7117∗𝑡 (−7.2918 cos 26,9971 𝑡 − 0.7320 sin 26,9971𝑡) + 7.2917

𝑥̈ (𝑡) = 𝑒 −2.7117∗𝑡 (0.0123 cos 26,9971 𝑡 + 198.8410 sin 26,9971𝑡)

Para 𝟎. 𝟏𝟔 < 𝒕 < 𝟎. 𝟑𝟐

𝑥(𝑡) = 𝑒 −2.7117𝑡 (0.4628 cos 26,9971 𝑡 − 0.3445 sin 26,9971𝑡) + 1.1667

𝑥̇ (𝑡) = 𝑒 −2.7117𝑡 (−10.5555 cos 26,9971 𝑡 − 11.5601 sin 26,9971𝑡)

𝑥̈ (𝑡) = 𝑒 −2.7117𝑡 (−283.4653 cos 26,9971 𝑡 + 316.3147 sin 26,9971𝑡)

Para 𝟎. 𝟑𝟐 < 𝒕 < 𝟎. 𝟔𝟒

𝑥(𝑡) = 𝑒 −2.7117𝑡 (0.1073 cos 26,9971 𝑡 − 0.8805 sin 26,9971𝑡) + 3.5537 − 7.2917𝑡

𝑥̇ (𝑡) = 𝑒 −2.7117𝑡 (−24.0619 cos 26,9971 𝑡 − 0.5091 sin 26,9971𝑡) − 7.2917

𝑥̈ (𝑡) = 𝑒 −2.7117𝑡 (51.5035 cos 26,9971 𝑡 + 650.9825 sin 26,9971𝑡)

Para 𝟎. 𝟔𝟒 < 𝒕 < 𝟎. 𝟖

𝑥(𝑡) = 𝑒 −2.7117𝑡 (−2.8956 cos 26,9971 𝑡 − 1.4871 sin 26,9971𝑡) − 5.8870 + 7.2917𝑡

𝑥̇ (𝑡) = 𝑒 −2.7117𝑡 (−32.2954 cos 26,9971 𝑡 + 82.20541 sin 26,9971𝑡) + 7.2917

𝑥̈ (𝑡) = 𝑒 −2.7117𝑡 (2306.8820 cos 26,9971 𝑡 + 648.9655 sin 26,9971𝑡)

 Desplazamiento Máximo
Para que 𝑥̇ sea máximo 𝑥̈ = 0
𝟎. 𝟏𝟔 < 𝒕 < 𝟎. 𝟑𝟐

𝑥̇ (𝑡) = 𝑒 −2.7117𝑡 (−10.5555 cos 26,9971 𝑡 − 11.5601 sin 26,9971𝑡) = 0

𝑡 = 0.2053 𝑠
Reemplazando ese valor en la ecuación de la velocidad tenemos

𝑥𝑚𝑎𝑥 = 𝑒 −2.7117𝑡 (0.4628 cos 26,9971 ∗ 0.2053 − 0.3445 ∗ 0.2053) + 1.1667

𝑥𝑚𝑎𝑥 = 1.4957[𝑐𝑚]
 Velocidad Máxima
Viendo la gráfica vemos que la velocidad máxima es

𝑥̇ 𝑚𝑎𝑥 = 14.304[𝑐𝑚/𝑠]

 Fuerza Cortante Máxima

𝑡𝑛
𝑄𝑚𝑎𝑥 = 𝑘 ∗ 𝑥𝑚𝑎𝑥 = 12 [ ] ∗ 1.4957[𝑐𝑚]
𝑐𝑚

𝑄𝑚𝑎𝑥 = 17.9480[𝑡𝑛]

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

t[s] x(t)[cm] x*(t) [cm/s] x**(t) [cm/s2]

0 0.0000 -0.0001 0.0123
0.0800 0.3284 10.0624 133.0891
0.1600 1.2581 9.5388 -119.0355
0.2400 1.3685 -6.5755 -112.8758
0.3200 0.9270 -0.3003 178.0852
0.4000 0.9218 -5.5332 -219.1969
0.4800 -0.0109 -13.3945 80.6574
0.5600 -0.6564 -2.9732 69.9030
0.6400 -0.9578 -7.1993 -114.7830
0.7200 -1.0960 9.9541 323.4805
0.8000 0.1868 14.3004 -214.9960

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

x**(t) [cm/s2]
00,400

00,300

00,200
x(t) [cm/s2]

00,100

00,000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
-00,100

-00,200

-00,300
t[s]

PROBLEMA 3

Determine el comportamiento dinámico del tanque de agua modelado en la figura bajo las características idealizadas
dadas, sujeto a una carga dinámica P(t) cuya variación se muestra en la figura. Asumiendo que el coeficiente de
amortiguación es el 20% de la amortiguación crítica y asumiendo la variación del tiempo igual a intervalos de 𝜋/64 .
Presente un cuadro de los valores de tiempo, desplazamiento, velocidad y aceleración desde t=0 hasta t= 𝜋/4
procesado por computador. Así mismo presente este cuadro en forma gráfica. Determine el máximo desplazamiento
en lo alto del tanque y determine la máxima fuerza cortante. Altura=18 m

𝑊 20[𝑡𝑛] 𝑠2
𝑚= = = 0.0204 [𝑡𝑛 ∗ ]
𝑔 980 [𝑐𝑚] 𝑐𝑚
𝑠2

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

𝑐 20 ∗ 𝐶𝑐𝑟𝑖𝑡
𝜉= = = 0.2
𝐶𝑐𝑟𝑖𝑡 100 ∗ 𝐶𝑐𝑟𝑖𝑡

15 𝑠
𝑐 = 2𝑚𝑤𝜉 = 2 ∗ 0.0204 ∗ √ ∗ 0.2 = 0.2213 [𝑡𝑛 ∗ ]
0.0204 𝑐𝑚

𝑡𝑛
𝑘 = 15[ ]
𝑐𝑚

0.0204(𝑡) + 0.2213𝑥̇ (𝑡) + 15𝑥(𝑡) = 𝑃(𝑡)

 Función complementaria
0.0204𝑠 2 + 0.2213𝑠 + 15 = 0
𝑠 = −5,4240 ± 26,5683𝑖

𝑥𝑐 (𝑡) = 𝑒 −5,4240𝑡 (𝐴 cos 26,5683 𝑡 + 𝐵 sin 26,5683𝑡)

 Integral Particular
Para P(t)
Del grafico
𝜋𝑡 𝜋𝑡
𝑃(𝑡) = 𝑃𝑜 cos = 5 cos
2𝑡𝑑 2𝜋/4

𝑃(𝑡) = 5 cos 2𝑡
La integral particular podría ser:

𝑥𝑝 (𝑡) = 𝐶 cos 2𝑡 + 𝐷 sin 2𝑡

𝑥𝑝̇ (𝑡) = −2𝐶 sin 2𝑡 + 2𝐷 cos 2𝑡
𝑥𝑝̈ (𝑡) = −4𝐶 cos 2𝑡 − 4𝐷 sin 2𝑡

Sustituyendo estos valores en la ecuación de movimiento tenemos:

0.0204(−4𝐶 cos 2𝑡 − 4𝐷 sin 2𝑡) + 0.2213(−2𝐶 sin 2𝑡 + 2𝐷 cos 2𝑡) + 15(𝐶 cos 2𝑡 + 𝐷 sin 2𝑡) = 5 cos 2𝑡

(14.9184𝐶 + 0.4426𝐷) cos 2𝑡 + (−0.4426𝐶 + 14.9184𝐷) sin 2𝑡 = 5 cos 2𝑡

14.9184𝐶 + 0.4426𝐷 = 5
−0.4426𝐶 + 14.9184𝐷 = 0

Luego 𝐶 = 0.3349 𝐷 = 0.0099

La solución completa es

𝑥𝑐 (𝑡) = 𝑒 −5,4240𝑡 (𝐴 cos 26,5683 𝑡 + 𝐵 sin 26,5683𝑡) + 0.3349 cos 2𝑡 + 0.0099 sin 2𝑡

Entonces

𝑥̇ (𝑡) = −5.4240𝑒 −5,4240𝑡 (𝐴 cos 26,5683 𝑡 + 𝐵 sin 26,5683𝑡)

+ 𝑒 −5,4240𝑡 (−26.5683𝐴 sin 26,5683 𝑡 + 26.5683𝐵 cos 26,5683𝑡) − 0.6697 sin 2𝑡 + 0.0199 cos 2𝑡

El sistema esta inicialmente en reposo

𝑥(0) = 0 = 𝐴 + 0.3349 𝑦 𝐴 = −0.3349

𝑥̇ (0) = 0 = −5.4240𝐴 + 26,5683𝐵 + 0.0199 𝑦 𝐵 = −0.0691

 Ecuación del movimiento

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

Por lo tanto para 𝟎 < 𝒕 < 𝝅/𝟒

𝑥(𝑡) = 𝑒 −5,4240𝑡 (−0.3349 cos 26,5683 𝑡 − 0.0691 sin 26,5683𝑡) + 0.3349 cos 2𝑡 + 0.0099 sin 2𝑡

𝑥̇ (𝑡) = 𝑒 −5,4240𝑡 (−0.0194 cos 26,5683 𝑡 + 9.2725 sin 26,5683𝑡) + 0.0199 cos 2𝑡 − 0.6698 sin 2𝑡

𝑥̈ (𝑡) = 𝑒 −5,4240𝑡 (246.4602 cos 26,5683 𝑡 − 49.7795 sin 26,5683𝑡) − 1.3396 cos 2𝑡 − 0.0396 sin 2𝑡

 Desplazamiento Máximo
Para que x sea máximo 𝑥̇ = 0

𝑥̇ (𝑡) = 𝑒 −5,4240𝑡 (−0.0194 cos 26,5683 𝑡 + 9.2725 sin 26,5683𝑡) + 0.0199 cos 2𝑡 − 0.6698 sin 2𝑡 = 0
𝑡 = 0.1173 𝑠

Reemplazando ese valor en la ecuación del desplazamiento tenemos

𝑥𝑚𝑎𝑥 = 𝑒 −5,4240∗0.1173 (−0.3349 cos 26,5683 ∗ 0.1173 − 0.0691 sin 26,5683 ∗ 0.1173) + 0.3349 cos 2 ∗ 0.1173
+ 0.0099 sin 2 ∗ 0.1173

𝑥𝑚𝑎𝑥 = 0.5043[𝑐𝑚]

 Fuerza Cortante Máxima

𝑡𝑛
𝑄𝑚𝑎𝑥 = 𝑘 ∗ 𝑥𝑚𝑎𝑥 = 15 [ ] ∗ 0.5043[𝑐𝑚]
𝑐𝑚

𝑄𝑚𝑎𝑥 = 7.5645[𝑡𝑛]

t[s] x(t)[cm] x*(t) [cm/s] x**(t) [cm/s2]

0 0.0000 0.0004 245.1206
0.0491 0.2156 6.8041 11.6256
0.0982 0.4791 2.6661 -140.7930
0.1473 0.4531 -3.0760 -65.2205
0.1963 0.2783 -3.0399 54.8231
0.2454 0.2097 0.2734 58.9763
0.2945 0.2680 1.5202 -9.7617
0.3436 0.3117 0.0122 -39.8754
0.3927 0.2720 -1.3907 -11.5834
0.4418 0.2041 -1.1288 17.6257
0.4909 0.1713 -0.2511 12.9237
0.5400 0.1667 -0.0962 -5.9998
0.5890 0.1508 -0.5884 -10.7470
0.6381 0.1121 -0.9109 -1.4130
0.6872 0.0693 -0.7778 5.2891
0.7363 0.0372 -0.5531 2.6652
0.7854 0.0111 -0.5516 -2.1761

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL

 Universidad Técnica de Oruro
Facultad Nacional de Ingeniería
Ingeniería Civil

UNIV: KEVIN PADILLA LOPEZ DINAMICA ESTRUCTURAL