Tarea: Relación de carga dinámica

Función cosenoidal

Santiago Salcedo Díaz

1000464950

Juan Sebastián Suarez Méndez

1003565112

Universidad Nacional de Colombia – Sede Bogotá

Facultad de Ingeniería
Análisis dinámico de estructuras
2019331
Ing. Maritzabel Molina Herrera
 26 de febrero de 2024

Rótula y manejo matricial para elementos tipo viga

1. Planteamiento del problema

Para el caso de estudio se tiene la siguiente geometría

Ilustración 1: Geometría y distribución de la edificación a analizar.

Partiendo de una relación de carga dinámica cosenoidal, realice

a) la demostración de su solución considerando que para un t=0 tiene un
Uo y Úo
b) Compare la relación de Umax, Fsmax, Fd max ty Fi max de los
siguientes casos
 Ome =0.1 wn
 Ome =0.5 wn
 Ome =0.9 wn
 Ome =1 wn
 Ome =1.1 wn
 Ome =2 wn
c) Compare el comportamiento de U,Fs,Fd,Fi entre la aceleración
cosenoidal con respecto a la sinusoidal para
 Ome =0.5 wn
 Ome =1 wn
 Ome =2 wn
d) Compare el comportamiento de U,Fs, Fi entre la aceleración cosenoidal
de la solución exacta determinada por MATLAB y la obtenida por SAP
para los siguientes casos
 Ome =0.5 wn
  Ome =1 wn
 Ome =2 wn

2. Demostración de solución a carga dinámica cosenoidal, para

t=0

Para realizar la demostración de la solución a carga dinámica cosenoidal con

las condiciones iniciales dadas es necesario comprender que la ecuación a
resolver es la del movimiento amortiguado forzado, en la cual la fuerza externa
que actúa sobre el sistema es regida por la función coseno, asi la ecuación a
resolver es la que se presenta a continuación:

m Ü +c U̇ + kU =F0 cos ( Ωt )
La cual al ser divida por m para dejar la segunda derivada sin factor se
reacomoda de la siguiente manera:

2 F0
Ü +2 ξ wn U̇ +w n U = cos(Ωt)
m
Teniendo en cuenta la forma de la ecuación se hace necesario encontrar una
solución homogénea y una solución particular, en donde la solución total de la
ecuación será igual a la suma de estas dos soluciones.
 Solución particular:
Para encontrar la solución particular lo primero es considerar que la forma de
esta y sus derivadas son las siguientes:

U p =A p cos ( Ωt )+ B p sen ( Ωt )

U̇ p =Ω ( − A p sen ( Ωt ) + B p cos ( Ωt ) )

Ü p =−Ω ( A p cos ( Ωt )+ B p sen ( Ωt ) )

2

Asi al reemplazar esta solución particular en la ecuación del movimiento

forzado amortiguado, se tiene:

F0
−Ω ( A p cos ( Ωt ) + B p sen ( Ωt ) ) +2 ξ wn Ω ( −A p sen ( Ωt ) + B p cos ( Ωt )) + w n ( A p cos ( Ωt ) + B p sen ( Ωt ) )=
2 2
cos (Ω t)
m
Para mayor facilidad primero se resuelven los valores asociados a seno y luego
los asociados a coseno:
o sen ( Ωt )

−Ω ( B p sen ( Ωt ) ) +2 ξ wn Ω ( − A p sen ( Ωt ) ) + wn ( B p sen ( Ωt ) )=0

2 2

2 2
−Ω ( B p ) +2 ξ wn Ω ( −A p ) + wn ( B p )=0
 2 2
−Ω ( B p ) 2 ξ w n Ω (− A p ) wn ( B p )
+ + =0
w2n w2n w2n
Teniendo en cuenta la siguiente igualdad se realiza las respectivas
simplificaciones:

Ω
=r
wn
2
−r ( B p ) +2 ξ r (− A p ) + ( B p )=0

B p ( 1−r 2) −2 ξ r A p =0
o cos ( Ωt )
F0
−Ω ( A p cos ( Ωt ) ) +2 ξ wn Ω ( B p cos ( Ωt ) ) + wn ( A p cos ( Ωt ) ) =
2 2
cos (Ω t)
m

2 2 F0
−Ω ( A p ) +2 ξ wn Ω ( B p ) +w n ( A p )=
m
2 2
−Ω ( A p ) 2 ξ wn Ω ( B p ) w n ( A p ) F
2
+ 2
+ 2
= 02
wn wn wn m wn
2
−r ( A p ) +2 ξ r ( B p )+ ( A p ) =U est

A p ( 1−r 2 ) +2 ξ r B p=U est

Ahora al tener las ecuaciones deducidas de los factores de seno y coseno, se
realiza una resta de ambas teniendo en cuenta un factor por el cual se
multiplicará cada ecuación, asi se tiene:

( A ¿ ¿ p ( 1−r 2 ) +2 ξ r B p=U est )∗(2 ξ r )¿

(B ¿ ¿ p ( 1−r 2 )−2 ξ r A p =0)∗(1−r 2 )¿

Al realizar esta resta entre las dos ecuaciones, se tiene la siguiente ecuación,
de la cual se despejará el factor B de la ecuación particular y posteriormente se
hallará el factor A:
2 2
B p ( 1−r ) −( 2 ξ r ) B p=2 ξ rU est
2

( 2 ξr ) U est
B p=
[ ( 2 ξr ) +( 1−r ) ]
2 2 2

( 1−r 2 ) U est
Ap=
[ ( 2ξr ) +( 1−r ) ]
2 2 2
Al colocar estos factores en la solución particular se obtiene la siguiente
ecuación:

( 1−r 2 ) cos ( Ω ) + ( 2 ξr ) sin ⁡(Ω)

U p= U est
[ ( 2ξr ) + 1−r ]
2
( 2 2
)
Esta ecuación será simplificada a través de los numero complejos, para esto es
necesario realizar una representación en el siguiente plano:

Al realizar un breve análisis tenemos que el ángulo phi es igual a:

2 ξr
tan(φ)= 2
1−r
Ahora para encontrar la simplificación de esta solución particular es necesario
encontrar la proyección del vector Up en la parta real de plano, de esta manera
por medio de sumas de vectores y teniendo en cuenta que esta proyección se
halla por medio del coseno del ángulo formado entre el vector Up y la línea de
proyección de los reales, se obtiene la siguiente simplificación de la solución
particular:

U est cos ( φ−Ωt )

U p=
√ ( 2 ξr ) + ( 1−r )
2 2 2

Al derivar esta solución particular se obtienen las siguientes derivadas:

U est sin ( φ−Ωt ) Ω

U̇ p =
√ ( 2 ξr ) + ( 1−r )
2 2 2

2
−U est cos ( φ−Ωt ) Ω
Ü p =
√ ( 2ξr ) +( 1−r )
2 2 2

 Solución homogénea
Para encontrar la solución homogénea lo primero a considerar es que la
solución tiene la siguiente forma, la cual proviene de la teoría de las
ecuaciones diferenciales, además sus derivas son las mostradas:
 U h =[ A d cos ( w d t ) + Bd sin ( wd t ) ] e
−ξ w nt

U̇ h =w d [ −A d sin ( wd t ) +B d cos ( wd t ) ] e
−ξ w nt
−ξ w n U h

Ü h =−w d U h−ξ w n w d [− A d sin ( w d t ) + Bd cos ( w d t ) ] e

2 −ξ w n t
−ξ w n U̇ h

Ahora teniendo en cuenta las condiciones iniciales se obtiene los factores Ad y

Bd ; asi de esta manera para t=0 se sabe que U o =0, por tanto:

U est cos ( φ−Ω∗0 )

U o =[ Ad cos ( w d∗0 ) + Bd sin ( wd∗0 ) ] e−ξ w ∗0 + n

√ (2 ξr ) +( 1−r )
2 2 2

U est cos ( φ )
U o =[ A d ] +
√ ( 2 ξr ) + ( 1−r )
2 2 2

U est cos ( φ )
Ad =U o−
√( 2 ξr ) +( 1−r )
2 2 2

Para obtener en facto B de esta ecuación homogénea se reemplazará en la

primera derivada de U que, para un tiempo inicial t=0, U̇ =U̇ o , asi se tiene:

U est sen ( φ−Ω∗0 )

U̇ o =w d [− Ad sin ( wd∗0 ) + B d cos ( w d∗0 ) ] e−ξ w ∗0−ξ w n U h+
n

√ ( 2 ξr ) + ( 1−r )
2 2 2

U est sen ( φ )
U̇ o =w d [ B d ]−ξ w n ( [ A d ]) + Ω
√ ( 2 ξr ) + ( 1−r )
2 2 2

( )√
U est cos ( φ ) U est sen ( φ )
U̇ o =w d [ B d ]−ξ w n U o − + Ω
√ ( 2 ξr ) + ( 1−r2 )
2 2
( 2 ξr ) + ( 1−r 2 )
2 2

( )√
U est cos ( φ ) U est sen ( φ )
U̇ o =w d [ B d ]−ξ w n U o − + Ω
√ ( 2 ξr ) + ( 1−r )
2 2 2
( 2 ξr )2+ ( 1−r 2 )
2

( )
U est cos ( φ ) U est sen ( φ )
w d B d=U̇ o +ξ wn U o− − Ω
√ ( 2ξr )2 + ( 1−r ) 2 2
√ ( 2 ξr )2+ ( 1−r 2 )
2

U̇ o+ ξ wn U o U est ξ wn cos ( φ ) +sen ( φ ) Ω

Bd = −
wd wd
√
( 2 ξr )2 + ( 1−r 2 )
2

( )
U̇ o+ ξ wn U o U est ξ wn cos ( φ )−Ω sin ⁡( φ)
Bd = −
wd wd
√
( 2 ξr )2+ ( 1−r 2 )
2
Al tener estos dos factores se reemplazan en la ecuación homogénea,
obteniendo asi lo siguiente solución:

U
h=¿
[( U o−
U est cos( φ )
√ (2 ξr ) +(1−r )
2 2 2 ) (
cos (w d t ) +
U̇ o+ ξw n U o
wd
−U est
(ξ wn cos ( φ ) +Ω sin ⁡(φ)
√ ( 2 ξr ) + ( 1−r )
2 22 )) ]
sin ⁡(w d t ) e
−ξ wn t
¿

Por último, cabe resaltar que al ser un movimiento forzado amortiguado, la

solución es la suma de solución particular y homogénea:

U ¿ U p+U h

3. Modificación del código en Mathlab

El código de MATLAB se desarrolla partir del ejercicio realizado en la clase para

fuerza externa sinusoidal por etapas, teniendo en cuenta que para este caso
deben modificarse la solución particular de desplazamiento, velocidad y
aceleración, además de las constantes homogéneas A y B.
Las propiedades de la estructura, características de carga, condiciones
iniciales, características dinámicas de la estructura y constantes de
movimiento se conservan iguales:

Por su parte, de acuerdo con la solución mostrada en el inciso 2 del presente

documento, las constante A y B para la primera etapa son:
Mientras que las constantes A y B para la segunda etapa se calculan como

En el proceso de análisis se establece la solución particular para fuerza externa

cosenoidal:

Del este ciclo se solicitan los vectores de fuerza de rigidez (Fs),

amortiguamiento (Fd) e inercial (Fi), el factor de carga dinámica y valores
máximos de fuerzas y desplazamientos.
Estos resultados se presentan en gráficas de fuerzas y desplazamientos.

4. Comparación de Umax, Fsmax, Fd max ty Fi max para

diferentes frecuencias de fuerza externa.

Del código de MATLAB presentado en el punto anterior, se calculan los valores

máximos de desplazamiento, fuerza de rigidez (Fs), fuerza de amortiguamiento
(Fd), fuerza inercial (Fi) y fuerza externa (Fext). Adicionalmente, teniendo la
rigidez de la estructura, la magnitud de la fuerza y el desplazamiento máximo,
se obtiene el factor dinámico de carga FDC, calculado como el Umax/Uest.
Tabla 1: Comparación de Umax, Fsmax, Fd max ty Fi max para diferentes frecuencias de fuerza
externa cosenoidal.

Umax Fsmax Fdmax Fimax Fextmax FDC Fs/Fext

0.1 0.008938 12864.996 648.621 7082.677 7082.465 1.816 1.816
0.5 0.006721 9673.888 958.732 9141.377 7082.649 1.366 1.366
0.9 0.028798 41450.080 3803.898 35518.372 7082.482 5.852 5.852
1 0.049122 70703.181 7075.684 70791.580 7082.649 9.983 9.983
1.1 0.025647 36915.050 3998.228 42916.014 7082.664 5.212 5.212
2 0.002915 4196.044 602.390 9671.571 7082.649 0.592 0.592

De acuerdo con la tabla 2, las condiciones más críticas para la estructura se

alcanzan cuando la frecuencia de fuerza externa es igual a la frecuencia
natural de la estructura, como era de esperarse, ya que se presenta una
condición cercana a la resonancia.
En cuanto a los desplazamientos, se observa que las mayores magnitudes se
presentan para frecuencias de fuerza externas cercanas a la frecuencia natural
de la estructura, en donde, para un r de 0.9 el FDC es de 5.852, para r=1 el
FDC es de 9.983 y para r=1.1 es de 5.212. Si embargo, la tendencia no es
simétrica para valores de Ome mayores y menores que Wn, hay mayores
desplazamientos cuando la frecuencia externa es menor a la natural. Se
observa que a pesar de que los valores de r de 0.1 y 2 tienen un delta similar
con respecto a r=1, el factor dinámico de carga para r=0.1 es 3 veces mayor
que el FDC para r=2.

Ilustración 2: Valores de fuerzas máximas externa y en estructura para diferentes valores de r.

En cuanto a las fuerzas actuantes de la estructura, se observa que, cuando la

frecuencia de fuerza externa Ome es menor a Wn, predomina la fuerza de
rigidez sobre la fuerza de amortiguamiento y la fuerza inercial, mientras que
Cuando Ome es mayor a Wn, hay una influencia mucho más marcada de la
fuerza inercial sobre la estructura. Por otro lado, cuando Ome es igual a Wn,
tanto la fuerza de rigidez, como la inercial, predominan con valores muy
cercanos, del orden de 10 veces la fuerza externa.
La fuerza de amortiguamiento presenta una participación muy baja para todas
las frecuencias de fuerza externa. Alcanza un valor máximo en la condición
crítica, muy cercano al de la fuerza externa.
Se identifica además que, el factor dinámico de carga es equivalente a la
relación entre la fuerza de rigidez de la estructura y la fuerza externa. Este
factor es importante porque indica la participación de la rigidez de la
estructura en la moderación de un fenómeno externo, con el cual se puede
buscar controlar los efectos desde el diseño de la estructura.

5. Comparación de U,Fs,Fd,Fi entre la aceleración cosenoidal

con respecto a la sinusoidal para diferentes frecuencias de
fuerza externa.

Para los dos casos de fuerza externa, el de aceleración cosenoidal y

aceleración senoidal, se analizan los desplazamientos y fuerzas de la
estructura para valores de r de 0.5, 1 y 2, a fin de comparar el comportamiento
de la estructura en los diferentes casos.
 Ome=0.5*wn
Fuerza externa

Ilustración 3: Fuerza función Sinusoidal y cosenoidal. Ome=0.5Wn.

Cabe resaltar que ambas fuerzas conservan la misma amplitud y

frecuencia, dado que se aplican para una misma región y estructura. Se
diferencian únicamente en el desplazamiento de un cuarto del periodo T Ω
en la escala del tiempo.
Desplazamiento

Ilustración 4: Desplazamiento función Sinusoidal y cosenoidal. Ome=0.5Wn.

En ambos casos, el desplazamiento ocurre en a la dirección opuesta a la

aplicación de la carga, es decir, en la dirección de una de las fuerzas de la
estructura. Posterior a los 4 segundos de aplicación, comienza a reducirse
hasta llegar al equilibrio estático.
Las mayores magnitudes se presentan en un rango de tiempo entre los 0.07
y 0.15 segundos y ambas funciones tienen un comportamiento no uniforme
durante los primero 1.3 segundos de aplicación de la carga, sin embargo, la
tendencia de la función sinusoidal es en general más uniforme que la
función cosenoidal.
Ilustración 5: Graficas de desplazamiento función Sinusoidal y cosenoidal. Ome=0.5Wn.

De acuerdo con las soluciones de desplazamiento, se observa que la

tendencia poco uniforme de la función con aceleración cosenoidal en los
primero 1.5 segundos se debe a que la solución particular tiene mayor
magnitud inicial, además de que no coincide en algunos tramos con la
dirección de desplazamiento de la solución homogénea.
Fuerzas

Ilustración 6: Graficas de fuerzas función Sinusoidal y cosenoidal. Ome=0.5Wn.

En las fuerzas de la estructura se observa en ambos casos una participación

predominante de la fuerza de rigidez que se opone a la fuerza externa
aplicada y que presenta una magnitud considerablemente mayor. La fuerza
inercial, por su parte, actúa en la misma dirección de la fuerza externa; esta
tiene mayor influencia en los primeros segundos y en la etapa de
movimiento libre, mientras que durante los 1.5 y 4 segundos tiene una
magnitud baja respecto a la fuerza externa y fuerza de rigidez. En cuanto a
la fuerza de amortiguamiento, esta actúa con un patrón independiente,
diferente al de la fuerza de rigidez y al de la fuerza externa; es constante
durante la aplicación de la carga y en la etapa de movimiento libre se
reduce paulatinamente.
Dado que el desplazamiento total en la función del coseno es mayor al
inicio de la carga, resulta un mayor aporte de la fuerza de rigidez en la
estructura que en el caso de la aceleración sinusoidal.
 Ome=1*wn

Ilustración 7: Fuerza función Sinusoidal y cosenoidal. Ome=1*Wn.

Para este caso se aplica el doble de ciclos de carga respecto al caso anterior
por el incremento en la frecuencia Ome, en donde este valor es igual a la
frecuencia natural de la estructura. La grafica de aceleración sinusoidal se
desplaza un cuarto del periodo T Ω en la escala del tiempo con respecto a la
función de aceleración cosenoidal.

Ilustración 8: Desplazamiento función Sinusoidal y cosenoidal. Ome=1*Wn.

A diferencia del caso Ome=0.5, en ambas funciones el desplazamiento es el

mínimo en cuanto inicia la aplicación de carga, hasta los 1.5 segundos, en
donde el desplazamiento máximo se mantiene constante, alcanzado valores
del orden de 0.05 y -0.05 m, para posteriormente entrar en etapa de
movimiento libre a los 4 segundos y reducir el desplazamiento.
Se observa que el desplazamiento ocurre en dirección contraria a la fuerza
externa, con lo cual se genera la fuerza de rigidez de la estructura.

Ilustración 9: Graficas de desplazamiento función Sinusoidal y cosenoidal. Ome=1*Wn.

En las gráficas de desplazamiento para aceleración sinusoidal y cosenoidal,

el desplazamiento de la ecuación particular es contrario al de la ecuación
homogénea, lo que ocasiona la reducción de desplazamientos para los
primeros 1.5 segundos de aplicación de la carga.
 Ilustración 10: Graficas de fuerzas función Sinusoidal y cosenoidal. Ome=1*Wn.

En las fuerzas de la estructura se observa en ambos casos una participación

predominante tanto de la fuerza de rigidez, como de la fuerza inercial, que
actúan opuestas entre sí y presentan además magnitudes mucho mayores
respecto a la fuerza externa, del orden de 10 veces. Para este caso, ni la
fuerza de rigidez, ni la inercial son opuestas a la fuerza actuante, la opuesta
es la de amortiguamiento, con magnitudes similares a las de la fuerza
externa.
En los primeros 1.5 segundos, el aporte de las fuerzas de la estructura no es
tan representativo, compartiendo la misma tendencia del desplazamiento
previamente ilustrado. Sin embargo, estas fuerzas adquieren un papel
importante con el transcurso de los segundos de la aplicación de carga y
una vez que el edificio deja de estar sometido a la fuerza externa.

 Ome=2*wn

Ilustración 11: Fuerza función Sinusoidal y cosenoidal. Ome=2*Wn.

Con el aumento de la frecuencia de la fuerza externa a 2 veces la

frecuencia natural Wn, a la estructura se la aplican un total de 59 ciclos de
carga La grafica en los casos de aceleración sinusoidal y cosenoidal. Cabe
resaltar que una gráfica esta desplazada respecto de la otra un cuarto
periodo T Ω en la escala del tiempo.
 Ilustración 12: Desplazamiento función Sinusoidal y cosenoidal. Ome=2*Wn.

De acuerdo con la gráfica de desplazamientos, las mayores diferencias para

los casos de fuerza con aceleración sinusoidal y cosenoidal están en los
intervalos de tiempo iniciales (primer segundo de aplicación de carga) y una
vez que el cuerpo deja de estar sometido a la fuerza externa. Se observa
que, para el caso de aceleración sinusoidal, el edificio está sometido a un
desplazamiento alterno, en donde luego de tener un amplio
desplazamiento, pasa uno muy corto, vuelve a un amplio desplazamiento y
sigue ese patrón hasta después de transcurrido 1 segundo de carga, en
donde se logra estabilizar. Además, una vez que se deja de aplicar la carga,
se presentan desplazamientos más grandes que los que se presentan con
aplicación de carga entre los segundos 3 y 4, posteriormente, ese
desplazamiento se regula por las fuerzas de la estructura.
En el caso de la aplicación de fuerza con aceleración cosenoidal, también se
presenta ese desplazamiento alterno, con la diferencia de que los
desplazamientos presentan una tendencia a tener mayores magnitudes
hacia el costado negativo. Una vez que deja de actuar la carga, comienza el
proceso de amortiguamiento, sin desplazamientos fuera del rango
alcanzado entre los segundos 3 y 4.
Ilustración 13: Graficas de desplazamiento función Sinusoidal y cosenoidal. Ome=2*Wn.

En las gráficas de desplazamiento para aceleración sinusoidal y cosenoidal,

el desplazamiento de la ecuación particular es contrario al de la ecuación
homogénea algunos instantes, pero también es igual en otros, esto pasa en
intervalos de tiempo muy pequeños en el primer segundo de la aplicación
de la carga, lo que ocasiona que la estructura tenga el extraño
comportamiento de desplazamientos alternos (grandes-pequeños). Pasado
ese intervalo de tiempo, el aporte de la solución particular empieza a ser
más determinante, mientras que el aporte de la solución homogénea pasa a
ser casi nulo.

Ilustración 14: Graficas de fuerzas función Sinusoidal y cosenoidal. Ome=2*Wn.

 En las fuerzas de la estructura se observa en ambos casos una participación
predominante de la fuerza inercial, que actúan en la dirección de la fuerza
externa y con una magnitud mayor, mientras que la fuerza de rigidez actúa
en la dirección contraria y con una magnitud cercana al 70% de la fuerza
externa. En cuanto a la fuerza de amortiguamiento, esta actúa con un
patrón independiente, diferente al de la fuerza de rigidez y al de la fuerza
externa; es constante durante la aplicación de la carga y en la etapa de
movimiento libre se reduce paulatinamente.

6. Comparación del comportamiento de U,Fs, Fi entre la

aceleración cosenoidal de la solución exacta determinada por
MATLAB y la obtenida por SAP.

Teniendo en cuenta que se desea realizar la comparación entre los resultados

obtenidos por el programa de análisis SAP y la obtenida por medio del código
de Matlab, se analizaran los desplazamientos, fuerzas de amortiguamiento y
fuerza inercial de la estructura para los valores de r de 0.5, 1 y 2. Cabe resaltar
que, para hallar las fuerzas de SAP se extrae el vector de velocidades para el
calculo de las fuerzas de amortiguamiento y se multiplican por el valor c,
mientras que para la fuerza inercial se extraen del programa las aceleraciones
y son multiplicadas por la masa total de la estructura.
A continuación, se presentarán de manera grafica los resultados a analizar,
esto para cada uno de valores de r.
 Ome=0.5*wn

Ilustración 15: Grafica de desplazamiento de SAP y Matlab. Ome=0.5*Wn.

 Ilustración 16: Grafica de fuerza de amortiguamiento de SAP y Matlab. Ome=0.5*Wn.

Ilustración 17: Grafica de fuerza inercial de SAP y Matlab. Ome=0.5*Wn.

 Ome=1*wn

Ilustración 18: Grafica de desplazamiento de SAP y Matlab. Ome=1*Wn.

 Ilustración 19: Grafica de fuerza de amortiguamiento de SAP y Matlab. Ome=1*Wn.

Ilustración 20: Grafica de fuerza inercial de SAP y Matlab. Ome=1*Wn.

 Ome=2*wn

Ilustración 21: Grafica de desplazamiento de SAP y Matlab. Ome=2*Wn.

 Ilustración 22: Grafica de fuerza de amortiguamiento de SAP y Matlab. Ome=2*Wn.

Ilustración 23: Grafica de fuerza inercial de SAP y Matlab. Ome=2*Wn.

Al realizar una comparación por medio de las graficas se evidencias grandes

similitudes entre los desplazamientos y las dos fuerzas para los tres valores de
r analizados, por medio de ambos programas. Al analizar un poco más a fondo
cada una de estas graficas se logra captar como las del valor de r=2, son las
que mayores discrepancias tienen entre los resultados obtenidos por SAP y
Matlab.
Asi para este valor de r=2, se denota como tanto en las gráficas de
desplazamientos como en las dos de fuerzas después del segundo 4 los valores
de obtenidos por medio de SAP son un poco mayores, asi mismo en la grafica
de fuerza inercial (Fi) se tienen como distorsiones, ya que a pesar de llegar a la
misma tendencia en cuanto a fuerzas estas fuerzas presentan un cambio de
tiempo liguero.
Después de este análisis de las tendencias en cuanto a las graficas se
revisaran los valores máximos de desplazamiento y el facto dinámico de carga
para cada uno de los tres valores de r, esto teniendo en cuenta la comparación
entre SAP y Matlab, de esta manera se tiene:
Tabla 3: Comparación de Umax y FDC para diferentes frecuencias de fuerza externa cosenoidal,
comparando SAP y Matlab.

Umax FDC
Sap Matlab Dif (SAP-Mat) Sap Matlab Dif (SAP-Mat)
0.5 0.006699 0.00672108 -0.0000221 1.361 1.366 -0.00448701
1 0.048827 0.0491221 -0.0002951 9.923 9.983 -0.059969912
2 0.002888 0.002915265 -0.0000273 0.587 0.592 -0.005540709

A partir de esta tabla se puede determinar que el mayor cambio en cuanto a

estos valores máximo entre los dos programas de análisis se da para el valor
de r=1, esto en cuanto a Umax, por tanto, también el FDC. Además de este
análisis de desplazamientos máximos se realiza un análisis de fuerzas
máximas, las cuales son:
Tabla 4: Comparación de Fd max y Fi max para diferentes frecuencias de fuerza externa
cosenoidal, comparando SAP y Matlab.

Fdmax Fimax
Sap Matlab Dif (SAP-Mat) Sap Matlab Dif (SAP-Mat)
0.5 957.0790043 958.7322641 -1.653259798 9100.732137 9141.376745 -40.64460876
1 7030.600647 7075.683829 -45.08318212 70314.11475 70791.57991 -477.4651572
2 597.3652529 602.3896658 -5.024412917 9540.451762 9671.570697 -131.1189353

Con esta información se puede determinar que al igual que en los

desplazamientos los mayores cambios entre los resultados obtenidos por SAP y
Matlab se presentan cuando el valor de r=1, esto teniendo cambios de mas de
10 veces con respecto al que menos presentan cambios, el cual es cuando
r=0.5; además es importante resaltar que el segundo que mas presenta
cambios es cuando r=2, siendo en este mas del triple de variación con
respecto a cuando r=0.5.
De este análisis también es importante resaltar que los valores mayores mas
altos tanto para desplazamientos como fuerzas se presentan cuando se realiza
el análisis por Matlab, esto debido a que la diferencia siempre da negativa.

7. Anexo

 Código de Matlab Implementado.

 Hoja de cálculo de Excel.