El documento presenta el uso de la integral de Duhamel para determinar la respuesta dinámica de una estructura sometida a cargas variables en el tiempo. Explica cómo la integral de Duhamel permite calcular la respuesta como la suma de los efectos de impulsos infinitesimales aplicados en diferentes momentos del tiempo. Además, describe cómo aplicar métodos numéricos como Simpson para aproximar la integral y así obtener la respuesta dinámica del sistema para cargas arbitrarias tanto en sistemas amortiguados como no amortiguados. Finalmente, presenta cál
El documento presenta el uso de la integral de Duhamel para determinar la respuesta dinámica de una estructura sometida a cargas variables en el tiempo. Explica cómo la integral de Duhamel permite calcular la respuesta como la suma de los efectos de impulsos infinitesimales aplicados en diferentes momentos del tiempo. Además, describe cómo aplicar métodos numéricos como Simpson para aproximar la integral y así obtener la respuesta dinámica del sistema para cargas arbitrarias tanto en sistemas amortiguados como no amortiguados. Finalmente, presenta cál
El documento presenta el uso de la integral de Duhamel para determinar la respuesta dinámica de una estructura sometida a cargas variables en el tiempo. Explica cómo la integral de Duhamel permite calcular la respuesta como la suma de los efectos de impulsos infinitesimales aplicados en diferentes momentos del tiempo. Además, describe cómo aplicar métodos numéricos como Simpson para aproximar la integral y así obtener la respuesta dinámica del sistema para cargas arbitrarias tanto en sistemas amortiguados como no amortiguados. Finalmente, presenta cál
El documento presenta el uso de la integral de Duhamel para determinar la respuesta dinámica de una estructura sometida a cargas variables en el tiempo. Explica cómo la integral de Duhamel permite calcular la respuesta como la suma de los efectos de impulsos infinitesimales aplicados en diferentes momentos del tiempo. Además, describe cómo aplicar métodos numéricos como Simpson para aproximar la integral y así obtener la respuesta dinámica del sistema para cargas arbitrarias tanto en sistemas amortiguados como no amortiguados. Finalmente, presenta cál
1) Determinar la respuesta dinámica de una estructura conociendo sus
características físicas y mecánicas. 2) Establecer una diferenciación entre la respuesta dinámica de sistema amortiguado al 5% y del mismo sistema asumiendo una amortiguación de 0 %. 3) Realizar un análisis de los resultados obtenidos de las respuestas dinámicas obtenidas para la carga dada. 3. Marco teórico. Dado que en muchas ocasiones las estructuras están sometidas a cargas que no son periódicas o armónicas. La integral de Duhamel o integral de convolución tiene como propósito calcular la respuesta de una estructura sometida a una excitación arbitraria. El procedimiento consiste en tratar el efecto de la fuerza como la superposición de impulsos infinitesimales como se indica en la figura 1.
Figura no 1
La respuesta al cabo de un instante genérico será igual a la suma (integral) de
los efectos producidos por los impulsos elementales P(τ) aplicados hasta ese instante. Cuando un sistema como el mostrado en la Figura 2 se somete a una excitación arbitraria expresada en términos de fuerza, como la indicada en la Figura 2, es posible dividirla en una serie de impulsos que se aplican en el tiempo т y que tienen una duración dт.
Figura No 2. A) Fig izq: Sistema lineal amortiguado.
B) Fig der: Excitación arbitraria. Al integrar el efecto de cada uno de estos impulsos diferenciales variando r; se obtiene para el caso sin amortiguamiento:
Ecuación 2. Ecuación de Duhamel para sistemas amortiguados.
Estas integrales se conocen como integrales de convolución o de Duhamel, y
corresponden a la solución particular del sistema. Si hay condiciones iniciales hay necesidad de adicionarles la solución homogénea, las siguientes ecuaciones, respectivamente:
Cabe destacar que la integral de Duhamel normalmente se la utiliza para tratar
impulsos o efectos transitorios, ya que para condiciones de una carga armónica en régimen ya se cuenta con la solución general tratada en los temas de movimientos oscilatorios armónicos amortiguados y no amortiguados.
La integral de Duhamel es un caso particular de la Integral de Convolución entre
dos funciones, la de carga y la de la respuesta a un impulso unitario. Aquellos casos en que la variación de la carga no es una función sencilla como para aproximarla por alguno de los casos cuya solución se conoce, la solución puede obtenerse evaluando la integral de Duhamel por algún procedimiento numérico (método de los trapecios, Simpson, etc.). Estrictamente, la Integral de Duhamel sólo resulta conveniente para calcular la respuesta en un instante dado perfectamente definido, es decir para un instante “t” dado. 3.1 Integral de Duhamel bajo el método de Metodos Numéricos.
3.1.1 Para sistemas NO amortiguados.
Es posible convertir la ecuación 1 en:
Donde:
Ecuación 3.
Ecuación 4.
Considerando primero la integración de la ecuación 3, calculamos el producto de
la siguiente ecuación:
𝑍(т) = F(т) ∗ cos(𝑤𝑡)
A intervalos de tiempo constantes Δ т como se muestra en la figura 3:
Figura 3. Discretizacion de la función para integración. Por lo tanto, el valor de Zn y Cn serán:
𝑍𝑛 = F(t = n ∗ Δ ∗ т) ∗ cos(w ∗ n ∗ Δ ∗ т) 𝐶𝑛 = C(t = n ∗ Δ ∗ т)
Cn se calcula aproximadamente como el área bajo la función Z(т), lo cual
puede lograrse utilizando la regla de Simpson, así:
La cual es válida para n=2,4,6…En forma consecutiva se convierte en:
El valor de Sn se calcula en forma analógica, solo que el cálculo de Zn se utiliza
sen (w*n* Δ*t) en vez de cos (w*n* Δ*t). Una vez calculado los valores de Cn y Sn, la respuesta para el punto n, se obtiene por medio de:
Es recomendable tomar como intervalo Δт un valor del orden de un ventiavo del
periodo de vibración del sistema (Δт = 𝑇/20). En la medida que se utilice un intervalo de tiempo menor, mayor será la precisión.
3.1.2 Para sistemas amortiguados.
Para el caso con amortiguamiento, la ecuación:
En: Donde 𝑤𝑑 = 𝑤 ∗ √1 − Ɛ2 . Es posible, entonces, convertir la ecuación 2 en:
Donde:
Por lo tanto, el valor de 𝑍𝑛 = F(t = nΔт) ∗ cos(wnΔт) y 𝐶𝑛 = C(t = nΔт). Cn se
calcula aproximadamente como el área bajo la función Z(т), lo cual puede lograrse utilizando la regla de Simpson, tomando en cuenta la disminución exponencial de la amplitud. La cual en forma recursiva es:
La cual es válida para n = 2, 4, 6… .El valor de Sn se calcula en forma análoga,
solo que en el cálculo de Zn se utiliza sen(𝑤𝑑 nΔт) en vez de cos(𝑤𝑑 nΔт). Una vez se han calculado los valores de Cn y Sn, la respuesta para el punto n se obtiene por medio de:
De igual manera que para el caso no amortiguado, es recomendable tomar corno
intervalo Δт un valor del orden de un veinteavo del período de vibración del 𝑇 sistema (Δт = 20). En la medida que se utilice un intervalo de tiempo menor, mayor será la precisión. 4. Cálculos y resultados.
Dadas las siguientes condiciones para el ejercicio, se determinó la respuesta
dinámica del sistema considerando métodos numéricos aplicados a la integral de Duhamel, así mismo, se calculó la respuesta dinámica con la aplicación directa de la integral de Duhamel para los 2 primeros tramos del diagrama de Fuerza VS tiempo bajo condiciones no amortiguadas. A continuación, los resultados:
4.1. Sistema No amortiguado.
Datos:
200 F t (kN) (seg) 150 150 0 0 100 150 0.5 -150 1.5 50 0 2 Fuerza (kN)
0 0 2.5 0 0.5 1 1.5 2 2.5 -50 DATOS INICIALES m (Kg) 10 -100 Rigidez (N/m) 9.87 -150 0% -200 Tiempo (s)
- De los resultados obtenidos de la aplicación directa de la integral de
Duhamel en los tramos 1 y 2, y comparándolos con los resultados de la integral de Duhamel por métodos numéricos (los tramos comprendidos entre 0 y 1 segundos), se logra concluir que la proximidad entre ambos resultados es muy alta, esto debido a los cortos intervalos de tiempo ∆τ que se usaron para el desarrollo de la integral por el método numérico. Como una medida porcentual del error relativo se puede estimar dicho parámetro en relación a 2 datos de la integral directa Vs integral por métodos numéricos, valores que corresponden a las amplitudes en el tiempo t=0.15 s entre ambos métodos de Duhamel:
𝑽. 𝒓𝒆𝒂𝒍 − 𝑽. 𝑪𝒂𝒍𝒄 ⌊𝟓𝟎. 𝟓𝟑𝒎 − 𝟓𝟑. 𝟑𝟎⌋
𝑬𝒓𝒓𝒐𝒓 = ∗ 𝟏𝟎𝟎% = ∗ 𝟏𝟎𝟎% = 𝟓. 𝟒𝟖𝟐% 𝑽. 𝒓𝒆𝒂𝒍 𝟓𝟎. 𝟓𝟑𝒎
- Bajo condiciones de amortiguamiento la amplitud del sistema tiende
a disminuir de manera oscilante – sub amortiguada-, cumplimento parámetros como el periodo calculado (T=2 s) y los demás parámetros físicos, como la frecuencia natural y la frecuencia amortiguada, después de la nulidad de la fuerza a partir del segundo 2. - Bajo condiciones de no amortiguamiento la amplitud del sistema tiende a llevar una oscilación constante luego de la nulidad de la fuerza a partir del segundo 2, cumpliendo además con todos los demás parámetros oscilatorios calculados para dicho sistema.
- Dado que las fuerza son muy grandes en relación a la pequeña
masa que tiene el sistema, las amplitudes del mismo son muy grandes, concluyendo que existe una relación intrínseca entre las variaciones en las amplitudes y las fuerzas que actúan sobre el sistema y las masas en el mismo.
