Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
PROBLEMA 3
Determine el comportamiento dinmico del tanque de agua modelado en la figura, bajo las caractersticas idealizadas dadas, sujeto a una carga
dinmica P(t) cuya variacin se muestra en la figura. Asumiendo que el coeficiente de amortiguacin es del 20% de la amortiguacin crtica y
asumiendo la variacin del tiempo igual a intervalos de /32. Presente un cuadro de los valores de tiempo, desplazamiento, velocidad y
aceleracin desde t=0 hasta t=/2 procesado por computadora. As mismo presente este cuadro en forma grfica. Determine el mximo
desplazamiento en lo alto del tanque y determine la mxima fuerza cortante. Altura= 13m.
DATOS
w 18000 ( kg)
kg
k 12000
cm
cm
g 980 2
s
"Problema [Link]"
SOLUCION
1) DIGRAMA DE CUERPO LIBRE Y ECUACION DE EQUILIBRIO
Ecuacin de equilibrio:
fi fc fs p ( t)
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
2) CALCULO DE LA MASA "m"
w 18000
m 18.367
g 980
3) CALCULO DEL COEFICIENTE DE AMORTIGUACION
k 12000
25.5604
m 18.3673
Ccr 2 m 2 18.367 25.56 938.953
La amortiguacin es el 20% de la amortiguacin crtica:
c 0.2 Ccr 187.791
4) ECUACION DE MOVIMIENTO
d d2
X ( t) X ( t) X ( t) X ( t)
dt 2
dt
Ecc_mov ( t) m X ( t) c X ( t) k X ( t) P ( t)
d2 d
Ecc_mov ( t) 12000 X ( t) 18.367 X ( t) 187.791 X ( t) P ( t)
2 dt
dt
Consideramos el sistema como vibracin libre:
2 2
mS cS k 0 187.791S
18.367S
12000 0
Soluciones de la ecuacin
2
c c 2 5.112 25.044i
S1
2m 2m w 5.112
Entonces sabemos que:
wd 25.044
2
c c 2 5.112 25.044i
S2
2m 2m
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
4) FUNCION COMPLEMENTARIA
Para vibracin libre con amortiguacin tenemos:
w t
Xc ( t) e ( A cos( wd t) B sin ( wd t) )
5.112 t
Xc ( t) e ( A cos( 25.044t
) B sin ( 25.044t
))
5) INTEGRAL PARTICULAR
T 2
2 wt 4
T
2
T
wt
Xp ( t) C sin ( 2 t)
Xp ( t) C sin ( 2 t)
Derivando tenemos la velocidad
d
Xp ( t) Xp ( t)
dt
Xp ( t) 2 C cos ( 2 t)
Derivando una vez ms tenemos la aceleracin
d2
Xp ( t) Xp ( t)
2
dt
Xp ( t) 4 C sin ( 2 t)
Reemplazamos las derivadas y la integral particular en la ecuacin de movimiento
Ec_mp ( t) m Xp ( t) c Xp ( t) k Xp ( t) Pp ( t)
Ec_mp ( t) 12000 C sin ( 2 t) 73.469 C sin ( 2 t) 375.581 C cos ( 2 t) Pp ( t)
Ec_mp ( t) simplify 375.581 C cos ( 2 t) 11926.531 C sin ( 2 t) Pp ( t)
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
6) COMPORTAMIENTO DE LA CARGA Pp(t)
td
Para Pp(t): 2
t
Pp ( t) Po sin
td
Pp ( t) Po sin ( 2 t)
Para t= /4
Entonces:
Pp 3 Po 3 Po 3
4
Entonces el comportamiento de la carga es:
t
Pp ( t) Po sin
td
Pp ( t) 3 sin ( 2 t)
Reemplazamos Pp(t) en la ecuacin para el valor de C
375.581 C cos ( 2 t) 11926.531 C sin ( 2 t) Pp ( t)
375.581 C cos ( 2 t) 11926.531 C sin ( 2 t) 3 sin ( 2 t)
Agrupando terminos:
4
C 11926.531C
sin ( 2 t) 3 sin ( 2 t) solveC 2.515 10
Xp ( t) C sin ( 2 t)
Xp ( t) C sin ( 2 t)
4
Xp ( t) 2.515 10 sin ( 2 t)
7) COMPORTAMIENTO DE LA ESTRUCTURA:
X ( t) Xc ( t) Xp ( t)
El desplazamiento ser:
4 5.112 t
X ( t) 2.515 10 sin ( 2 t) e ( A cos( 25.044t
) B sin( 25.044t
))
4 5.112 t
X ( t) 2.515 10 sin ( 2 t) e ( A cos( 25.044t
) B sin ( 25.044t
))
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
La velocidad ser:
d
X ( t) X ( t)
dt
4 5.112 t 5.112 t
X ( t) 5.03 10 cos( 2 t) e ( 25.044 A sin ( 25.044t
) 25.044B
cos( 25.044t
) ) 5.112 e ( A cos( 25.044t
) B sin ( 25.044t
))
Usamos las condiciones iniciales para calcular A Y B:
En t=0 el sistema est en reposo:
X ( 0) 0A 0
A X ( 0) 0 solve A 0 0
4
X ( 0) 0 5.112 A 25.044B
5.03 10 0
5 5
B X ( 0) 0 solveB 0.204 A 2.008 10 2.008 10
Reemplazamos estos valores en la educacin del comportamiento
Desplazamiento
4 5.112 t
X ( t) 2.515 10 sin ( 2 t) e ( A cos( 25.044t
) B sin ( 25.044t
))
4 5 5.112 t
X ( t) 2.515 10 sin ( 2 t) 2.008 10 sin ( 25.044t
) e
4 5 5.112 t
X ( t) simplify 2.515 10 sin ( 2 t) 2.008 10 sin ( 25.044t
) e
Velocidad
d
X ( t) X ( t)
dt
4 4 5.112 t 4 5.112 t
X ( t) 5.03 10 cos( 2 t) 1.027 10 sin ( 25.044t
) e 5.03 10 cos( 25.044t
) e
4 4 5.112 t 4 5.112 t
X ( t) simplify 5.03 10 cos( 2 t) 5.03 10 cos( 25.044t
) e 1.027 10 sin ( 25.044t
) e
Aceleracin
d2
X ( t) X ( t)
2
dt
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
3 5.112 t 4 5.112 t 3 5.112 t
X( t) 1.006 10 sin ( 2 t) 0.013 sin ( 25.044t
) e 5.249 10 sin ( 25.044t
) e 5.143 10 cos( 25.044t
) e
5.112 t 5.112 t
X( t) simplify 0.005 cos( 25.044t
) e 0.001 sin ( 2 t) 0.012 sin ( 25.044t
) e
Finalmente traficamos el desplazamiento, velocidad y aceleracin desde t= 0 hasta t= /2
t 0
32 2
t
t X ( t) X ( t) X ( t)
0 0 0 0 5.14310-3
1/32 0.098 4.13910-5 7.68810-4 2.00110-3
1/16 0.196 1.03510-4 3.90310-4 -4.33310-3
3/32 0.295 1.35810-4 3.87110-4 2.34410-3
1/8 0.393 1.78910-4 4.12110-4 -1.99110-3
5/32 0.491 2.09610-4 2.37810-4 -6.98410-4
3/16 0.589 2.31510-4 2.10910-4 -5.910-4
7/32 0.687 2.47310-4 9.6110-5 -1.35610-3
1/4 0.785 2.51210-4 -4.83810-6 -7.83410-4
9/32 0.884 2.46710-4 -9.28410-5 -1.0610-3
5/16 0.982 2.32410-4 -1.95710-4 -9.41810-4
11/32 1.08 2.0910-4 -2.78410-4 -7.97810-4
3/8 1.178 1.77910-4 -3.55510-4 -7.43110-4
13/32 1.276 1.39710-4 -4.18810-4 -5.43210-4
7/16 1.374 9.62410-5 -4.64310-4 -3.88110-4
15/32 1.473 4.90710-5 -4.93610-4 -1.99110-4
1/2 1.571 -6.52310-9 -5.0310-4 3.80510-6
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
t
TIEMPO - DESPLAZAMIENTO
4
3 10
t X ( t)
0 0 0
1/32 0.098 4.13910-5 4
2.499988 10
1/16 0.196 1.03510-4
3/32 0.295 1.35810-4
1/8 0.393 1.78910-4 4
1.999977 10
5/32 0.491 2.09610-4
3/16 0.589 2.31510-4
X(t)
7/32 0.687 2.47310-4 1.499965 10
4
1/4 0.785 2.51210-4
9/32 0.884 2.46710-4
5/16 0.982 2.32410-4 5
9.999533 10
11/32 1.08 2.0910-4
3/8 1.178 1.77910-4
5
13/32 1.276 1.39710-4 4.999417 10
7/16 1.374 9.62410-5
15/32 1.473 4.90710-5
9
1/2 1.571 -6.52310-9 7 10
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
t [RAD]
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
t
t X ( t) TIEMPO - VELOCIDAD
0 0 0 3
1 10
1/32 0.098 7.68810-4
1/16 0.196 3.90310-4
3/32 0.295 3.87110-4 4
6.66667 10
1/8 0.393 4.12110-4
5/32 0.491 2.37810-4
4
3/16 0.589 2.10910-4 3.33333 10
7/32 0.687 9.6110-5
1/4 0.785 -4.83810-6
X'(t)
9/32 0.884 -9.28410-5
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
5/16 0.982 -1.95710-4
11/32 1.08 -2.78410-4
4
3/8 1.178 -3.55510-4 3.33333 10
13/32 1.276 -4.18810-4
7/16 1.374 -4.64310-4
4
15/32 1.473 -4.93610-4 6.66667 10
1/2 1.571 -5.0310-4
3
1 10
t [RAD]
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
t
TIEMPO - ACELERACION
t X ( t)
3
0 0 5.14310-3 6 10
1/32 0.098 2.00110-3
1/16 0.196 -4.33310-3
3
3/32 0.295 2.34410-3 4 10
1/8 0.393 -1.99110-3
5/32 0.491 -6.98410-4
3
3/16 0.589 -5.910-4 2 10
7/32 0.687 -1.35610-3
1/4 0.785 -7.83410-4
X''(t)
9/32 0.884 -1.0610-3
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
5/16 0.982 -9.41810-4
11/32 1.08 -7.97810-4
3/8 1.178 -7.43110-4 2 10
3
13/32 1.276 -5.43210-4
7/16 1.374 -3.88110-4
15/32 1.473 ... 3
4 10
1/2 1.571
3
6 10
t [RAD]
Univ. Lazzo Titichoca Toms A.
Universidad Tcnica de Oruro
Facultad Nacional de Ingeniera
Ingeniera Civil
El mximo desplazamiento ser donde la velocidad se hace 0 que por la grfica podemos observar que est cerca de t=0.8
Given t 0.8
X ( t) 0
tmax Find ( t) 0.779252032
4
max X ( tmax) 2.5125 10
El cortante mximo ser:
1 1 4 4
Qmax k max 1.2 10 2.512 10
2 2
Qmax 1.507 ( ton )
Univ. Lazzo Titichoca Toms A.