Ejemplo S Integral de Duhamel
Ejemplo S Integral de Duhamel
Ejemplo S Integral de Duhamel
EJEMPLOS
Integral de Duhamel para un sistema sin amortiguamiento
Ejemplo 7.11
Calcular la respuesta dinmica del tanque de agua de la Figura 7.3, el cual est sujeto a una carga explosiva cuyo histograma de fuerza se muestra en la misma figura.
histograma de carga p(t) p(t) 96.6 k k=2700 k/ft
w=96.6 k
0.025 s fs
Figura 7.1
0.025 s
Solucin Para la resolucin de este problema se utiliza a continuacin Mathcad 2000, el cual es un programa de anlisis matemtico que hace ms fcil la resolucin de integrales de este tipo. Clculos adicionales Gravedad [ft/s2]: Frecuencia natural: Periodo natural:
g : 32.3
n
k g w
n
Tn : 2
Primera fase, para 0<t<0.025 El histograma de carga para esta fase est definido por la recta, de coordenadas [a,b] y [c,d], para la evaluacin de la carga en funcin del tiempo. a: 0 b: 0 c : 0.025 d : 96.6 la ecuacin de la recta resultante es: d b p ( x) : ( x a) b 3864.0 x c a Evaluacin de A(t) y B(t) para el clculo de la integral de Duhamel:
72
t
A(t ) :
p ( x ) cos(
x)dx
B(t ) :
p ( x ) sen (
x)dx
t ) B(t ) cos(
t)
f (t ) : k u (t )
La respuesta de velocidad es:
2
v (t ) :
d dt
u (t )
Segunda fase, para 0.025<t<0.05 El histograma de carga para esta fase est definido por la recta, de coordenadas [e,g] y [h,i], para la evaluacin de la carga en funcin del tiempo. e: 0 g : 96.6 h : 0.025 i: 0 la ecuacin de la recta resultante es: i g p( ) : ( e) g 3864.0 96.6 h e las condiciones iniciales para esta fase son: tiempo inicial [s]: desplazamiento inicial[ft]: velocidad inicial [ft/s]: Evaluacin de 3C(t) y D(t) para el clculo de la integral de Duhamel:
j j
tr : 0.025
u (tr )
v (tr )
3.271 10
0.385
C( j) :
p(
cos(
)d
D( j ) :
p(
sen (
)d
j)
32.3 w n
C ( j ) sen (
j ) D( j ) cos(
j)
fuerza ( j ) : k res ( j )
La respuesta de velocidad es:
vel( j ) :
d dj
res ( j )
Tercera fase, para vibracin libre t>0.05 las condiciones iniciales para esta fase son: tiempo inicial [s]: desplazamiento inicial[ft]:
2 3
0.017
(t)=v(t) C(t)=C(j)
73
velocidad inicial [ft/s]: La respuesta de desplazamiento para vibracin libre es: vel(0.05 tr ) reslib ( s ) : sen (
n
vel(0.05
tr )
0.563
s ) res (0.05
tr )
cos(
s)
flib( s ) : k reslib ( s )
Respuesta mxima:
Fuerza[k]=69.214
en un tiempo [s]=0.0772
Respuesta de Desplazamiento
0.04 0.03 0.02
desplazamiento [ft]
0.01 0 0.01 0.02 0.03 0.04 0 0.05 0.1 tiempo [s] 0.15 0.2 0.25
74
Calcular la respuesta dinmica del tanque de agua de la Figura 7.4 que tiene una razn de amortiguamiento =5%, el cual est sujeto a una carga explosiva cuyo histograma de fuerza se muestra en la misma figura.
histograma de carga p(t) p(t) 96.6 k k=2700 k/ft
w=96.6 k
0.025 s fs
Figura 7.2
0.025 s
Solucin Clculos adicionales Gravedad [ft/s2]: Frecuencia natural: Periodo natural: Razn de amortiguamiento: Frecuencia de amortiguamiento:
g : 32.3
n
k g w
n
Tn : 2
: 0.05
D
Primera fase, para 0<t<0.025 El histograma de carga para esta fase est definido por la recta, de coordenadas [a,b] y [c,d], para la evaluacin de la carga en funcin del tiempo. a: 0 b: 0 c : 0.025 d : 96.6 la ecuacin de la recta resultante es: d b p ( x) : ( x a) b 3864.0 x c a Evaluacin de A(t) y B(t) para el clculo de la integral de Duhamel:
75
x t t
A(t ) :
p ( x)
e e
n n
cos(
x)dx
B( t ) :
p ( x)
e e
n n
x t
sen (
x)dx
t ) B(t ) cos(
t)
f (t ) : k u (t )
La respuesta de velocidad es:
5
v (t ) :
d dt
u (t )
Segunda fase, para 0.025<t<0.05 El histograma de carga para esta fase est definido por la recta, de coordenadas [e,g] y [h,i], para la evaluacin de la carga en funcin del tiempo. e: 0 g : 96.6 h : 0.025 i: 0 la ecuacin de la recta resultante es: i g p( ) : ( e) g 3864.0 96.6 h e las condiciones iniciales para esta fase son: tiempo inicial [s]: desplazamiento inicial[ft]: velocidad inicial [ft/s]: Evaluacin de 6C(t) y D(t) para el clculo de la integral de Duhamel:
j
tr : 0.025
u (tr )
v (tr )
3.211 10
0.376
C( j ) :
p(
e
)
n n
j j
cos(
)d
D( j ) :
p(
e
)
n n
sen (
)d
u (tr ) cos(
j)
v (tr )
D
u (tr )
sen (
j)
32.3 C ( j ) sen ( w D
j ) D( j ) cos(
j)
fuerza ( j ) : k res ( j )
La respuesta de velocidad es:
vel( j ) :
d dj
res ( j )
Tercera fase, para vibracin libre t>0.05 las condiciones iniciales para esta fase son:
5 6
(t)=v(t) C(t)=C(j)
76
tiempo inicial [s]: desplazamiento inicial[ft]: velocidad inicial [ft/s]: La respuesta de desplazamiento para vibracin libre es:
reslib ( s ) : e
n
0.017 0.52
tr )
res (0.05
tr )
cos(
s)
vel(0.05
res (0.05
tr )
sen (
s)
flib( s ) : k reslib ( s )
las graficas de respuesta en las tres fases son:
Respuesta de Fuerza Elstica
100 75 50
Respuesta mxima:
Fuerza[k]=64.1402
en un tiempo [s]=0.0758
Respuesta de Desplazamiento
desplazamiento [ft]
0.0075 0 0.0075 0.015 0.0225 0.03 0 0.1 0.2 0.3 0.4 tiempo [s] 0.5 0.6
77
Respuesta mxima:
Desplazamiento[ft]=0.023756
en un tiempo [s]=0.075