PROYECTO DE AULA - PARTE I

Autores:
Mateo Zapata Garcia
José Luis González Castro

Cedulas
1000893956
1128451629

Profesor
John Faber Dávila Illesca

Asignatura
Análisis estructural II
 Tabla de contenido
1. Introducción

2. Solución de problemas

3. Archivos adjuntos

Introducción
En el presente documento se exponen las soluciones detalladas a una serie de problemas
complejos, los cuales han sido abordados utilizando técnicas avanzadas y de última
generación en el campo de la ingeniería y la ciencia. Las soluciones han sido
acompañadas por gráficas y algoritmos, los cuales fueron diseñados y optimizados para
obtener los resultados más precisos y eficientes posibles. Además, se incluyen los
archivos y datos requeridos para llevar a cabo la implementación y ejecución de los
algoritmos, así como para realizar las deducciones necesarias en cada uno de los
problemas planteados. Todo esto ha permitido obtener una solución óptima y rigurosa a
cada uno de los problemas planteados.

Solución de problemas
En la sección de solución de problemas de este documento, se exponen de manera
detallada las soluciones obtenidas para cada uno de los problemas planteados. Para ello,
se han aplicado técnicas y herramientas avanzadas en el campo de la ingeniería y la
ciencia, con el objetivo de obtener soluciones precisas y rigurosas. En cada caso, se ha
realizado una descripción detallada del problema, se han expuesto los pasos y algoritmos
utilizados para su resolución, y se han presentado los resultados obtenidos. Además, se
han realizado análisis y comparaciones con soluciones previas, en caso de ser necesario,
para validar la solución propuesta. La exposición rigurosa de la solución a cada problema
es esencial en el contexto de la investigación científica y la ingeniería, ya que permite la
comprensión clara y precisa de los resultados, y facilita la verificación
P1. Demuestre analíticamente que la frecuencia natural de un sistema SDOF es 𝑤𝑛 =
√𝐾/𝑀. Sugerencia: Resuelva la ecuación diferencial de movimiento homogénea para un
sistema NO amortiguado.
La ecuación diferencial de movimiento para un sistema de un grado de libertad (SDOF) no
amortiguado se puede escribir como:
𝑀𝑥̈ + 𝐾𝑥 = 0
Donde 𝑀 es la masa del sistema y 𝐾 es la constante de fuerza del sistema.
Para encontrar la solución a esta ecuación diferencial, suponemos una solución del tipo:
𝑥(𝑡) = 𝐴sin(𝑤𝑛𝑡 + 𝜙)
Donde 𝐴 y 𝜙 son constantes de integración y 𝑤𝑛 es la frecuencia natural del sistema.
Sustituyendo esta solución en la ecuación diferencial de movimiento, obtenemos:
𝑀(−𝑤^2𝑛𝐴sin(𝑤𝑛𝑡 + 𝜙)) + 𝐾𝐴sin(𝑤𝑛𝑡 + 𝜙) = 0
Simplificando, obtenemos:
(𝐾/𝑀)𝐴sin(𝑤𝑛𝑡 + 𝜙) + 𝑤^2𝑛𝐴sin(𝑤𝑛𝑡 + 𝜙) = 0
Factorizando 𝐴sin(𝑤𝑛𝑡 + 𝜙), obtenemos:
(𝐾/𝑀 + 𝑤^2𝑛)𝐴sin(𝑤𝑛𝑡 + 𝜙) = 0
Como 𝐴sin(𝑤𝑛𝑡 + 𝜙) es diferente de cero, podemos dividir ambos lados de la ecuación por
𝐴sin(𝑤𝑛𝑡 + 𝜙):
(𝐾/𝑀 + 𝑤^2𝑛) = 0
Resolviendo para 𝑤𝑛, obtenemos:
𝑤𝑛 = √𝐾/𝑀
Lo que demuestra que la frecuencia natural de un sistema SDOF es 𝑤𝑛 = √𝐾/𝑀.

P2. Exprese la ecuación diferencial de movimiento de 1 G.D.L para los siguientes

sistemas en vibración libre (Esto requerirá que obtenga la rigidez equivalente en cada
caso)
a.

Pórtico plano con diagrama rígido

Primero se buscara la rigidez equivalente del sistema
Rigidez columna 1 (pequeña)
12 EI
k 1= 3
H
Rigidez columan 2 (grande)
12 EI 12 EI 3 EI
k 2= 3
= 3
= 2
(H + H) 8H 2H
Como el sistema es paralelo entonces:

kequi=2 ( k 1 )+ 2 ( k 2 ) =2
( 12HEI )+2( 32 HEI )= 27HEI
3 3 3

como sabemos que un sistema de vibración libre se mueve con un movimiento armónico simple
partiremos de estas ecuaciones:

x (t )= A∗Sen ( wt +θ 1 ) y x ( t )=B∗Sen(wt +θ 2)
Y en estos sistemas actúa una fuerza que es:

F=−K∗X
Donde K es la rigidez y X el desplazamiento

También sabemos que F =m*a donde m es la masa y a es la aceleración entonces

m∗a=−k∗X
2
d x
También sabemos que la segunda derivada de posición respecto a tiempo a = 2 entonces
dt
2
m∗d x
2
=−KX , despejando tendremos
dt
2
d X kx
2
+ =0 (1)
d t m
De esta ecuación podemos verlo de la siguiente forma t 2+ at + b=0 (2)

Y de esta ecuación se puede escribir de esta forma

t 1=α + β y t 2=α−β entonces se puede escribir esta ecuación
αt αt
y= A e ∗sen ( β∗t ) +B e ∗cos ( β∗t ) donde e αt =0
Entonces de 2 y 1 podemos ver que
k
m t +0 t+ k=0 y wn=
2
podemos sacar esta ecuación
m
x (t )= A∗sen ( w∗t ) +B∗cos ( w∗t ) remplazando tenemos que

U ( t )= A∗sen (√ mk ∗t )+ b∗cos( √ mk ∗t ) y remplazando k por el kequi tenemos

U ( t )= A∗sen
(√ 27 EI
3
H m )
∗t +b∗cos
(√ √ 27 EI
3
H m )
∗t →ecuación diferencial de movimiento.

B.
Sistema masa resorte en serie
Primero se buscara la rigidez equivalente del sistema
p
Primero sabemos que la rigidez de un resorte es p=u∗k , despejando tenemos U =
k
Como son 2 resortes:
P P
U 1= y U 2=
k1 k2

Como el sistema es en serie podemos decir que U=U1+U2

Remplazando

U=
P P
+
k1 k 2
entonces U =
1
+ (
1
k1 k2
p )
Al final la rigidez equivalente seria

( )
−1
1 1 k 1+ K 2
Kequi= + =
k1 k2 K 1∗K 2
De la demostración del punto a. sacamos esta ecuacion

y ( t ) =A∗sen ( w∗t ) + B∗cos ( w∗t ) remplazando tenemos que

U ( t )= A∗sen (√ mk ∗t )+ b∗cos( √ mk ∗t ) y remplazando k por el kequi tenemos

U ( t )= A∗sen (
movimiento.
√ k 1∗K 2∗m
K 1+ K 2
∗t )+ b∗cos (
√ k 1∗K 2∗m ∗t )
k 1+ k 2 →ecuación diferencial de

C.
Pórtico plano con diagrama rígido y arriostrado con diagonales que solo toleran tracción
Primero se buscará la rigidez equivalente del sistema
Rigidez columna
12 EcI
kc= 3
H
Rigidez riosta
As∗Es 2
kriosta= ∗cos θ
Lriosta

( 24HEcI )+( 2∗As∗Es

kequi=2 ( kc )+ 2 ( kriosta ) = 3
Lriosta
∗cos θ )
2

como sabemos que un sistema de vibración libre se mueve con un movimiento armónico simple
partiremos de estas ecuaciones:

x (t )= A∗Sen ( wt +θ 1 ) y x ( t )=B∗Sen(wt +θ 2)
Y en estos sistemas actúa una fuerza que es:

F=−K∗X
Donde K es la rigidez y X el desplazamiento

También sabemos que F =m*a donde m es la masa y a es la aceleración entonces

m∗a=−k∗X
2
d x
También sabemos que la segunda derivada de posición respecto a tiempo a = 2 entonces
dt
2
m∗d x
2
=−KX , despejando tendremos
dt
2
d X kx
2
+ =0 (1)
d t m
De esta ecuación podemos verlo de la siguiente forma t 2+ at + b=0 (2)
Y de esta ecuación se puede escribir de esta forma
t 1=α + β y t 2=α−β entonces se puede escribir esta ecuación
αt αt
y= A e ∗sen ( β∗t ) +B e ∗cos ( β∗t ) donde e αt =0

Entonces de 2 y 1 podemos ver que

k
m t +0 t+ k=0 y wn=
2
podemos sacar esta ecuación
m
x (t )= A∗sen ( w∗t ) +B∗cos ( w∗t ) remplazando tenemos que

U ( t )= A∗sen (√ mk ∗t )+ b∗cos( √ mk ∗t ) y remplazando k por el kequi tenemos

(√( ( ∗cos θ )∗1

) (√ ( ( ∗cos θ)∗1
)
2∗As∗Es 2 2∗As∗Es 2

) )
24 EcI Lriosta 24 EcI Lriosta
U ( t )= A∗sen + ∗t + b∗cos + ∗t
H3 m H3 m
→ecuación diferencial de movimiento.

P3. Si se aplica una carga armónica 𝑃(𝑡) = 100𝑘𝑁 ∗ 𝑠𝑒𝑛 (𝜋 𝑡 ∗ #) en la misma dirección del
grado de libertad u en el sistema mostrado en P2.a, determine el lado de la sección
cuadrada de las columnas que produciría la máxima respuesta de desplazamientos en el
sistema. Use M=10E3 kg, E=25 GPa, ξ=0.05*#; L=2*# m y h=4m
Recordemos que #=12

Datos iniciales:
M=10*10^3 kg
E=25*10^6 KN/m2
ξ=0.05*12=0.6
L=2*#=24m
H=4m
Empezamos buscando la rigidez equivalente del sistema, podemos ver que el sistema
está compuesto por 4 columna (2 del mismo tipo)
Rigidez columna 1 (pequeña)
12 EI
k 1= 3
H
Rigidez columna 2 (grande)
12 EI 12 EI 3 EI
k 2= 3
= 3
= 2
(H + H) 8H 2H
Como el sistema es paralelo entonces:

kequi=2 ( k 1 )+ 2 ( k 2 ) =2
( 12HEI )+2( 32 HEI )= 27HEI
3 3 3

Ahora usando la fórmula de frecuencia natural tenemos

wn=
√ kequi ;
m
2
kequi=wn ∗m

wf rd= 1
También tenemos que β= y
wn √(1−β ) +(2 ξ β)2
2 2

Tenemos p(t)=pm*sen(wf*t)
De la ecuación de 𝑃(𝑡) = 100𝑘𝑁 ∗ 𝑠𝑒𝑛 (π 𝑡 ∗ 12) podemos sacar que:
Pm=100KN
Wf=12 π rad/s
Como el problema nos pide calcular la sección de la columna en la respuesta máxima
usamos esta ecuación:

rd=
√ 1
2
1
= =
1
(2 ξ) 2 ξ 2∗0.6
=0.833

Ya con el valor de rd podemos buscar el valor de β

1 1
rd= → 0.833=
√(1−β 2 2
) +(2 ξ β) 2
√(1−β ) +( 2(0.6) β)2
2 2

Despejando se obtiene estos 2 valores β=0.315 y β=0.604

Escogemos el mayor y con la ecuación de β encontramos wn

 wf 12 π 12 π
β= → 0.604= → Wn= =62.416 rad / seg
wn wn 0.604
Ahora

2 2 3 m
kequi=wn ∗m=62.416 ∗10∗10 kg∗9.81 =382179.7672 KN /m
s2
Ya con la ecuación de la rigidez equivalente tenemos
27 EI b
4
kequi= 3 y como I = entonces
H 12

kequi=
27 E b 4
12 H 3
√
→ b= 4
12∗64∗382179.7672
27∗25∗106
→ b=0.812m

B=81.2 cm → lado de la sección de la columna

P4. Si se aplica un movimiento armónico en la base 𝑢𝑔(𝑡) = 0.25 ∗ 𝑠𝑒𝑛 (# ∗ 𝜋𝑡), en metros,
en la misma dirección del grado de libertad u en el sistema mostrado en P2.c, encuentre
la mínima masa requerida para alcanzar la fluencia en las riostras. Para esto considere:
Arriostramiento con cables de acero (Es=29000 ksi) con 1” de diámetro y fy=36000 psi.
Columnas de concreto (fc’=3000 psi) con sección rectangular de 0.25x0.25 m.
Amortiguamiento del 5% Altura H= 3m; Separación entre columnas L= 4m;
Datos iniciales

RIOSTAS
Es =29000Ksi =199947.96 MPa
Diámetro =1 in = 0.0254 m
Fy=36000psi = 250 Mpa
COLUMNAS
Sección 0.25x0.25 m → b=0.25m
Fc=3000psi = 21 Mpa

Ec=4700 √ 21 = 21538.1 Mpa

H=3 , L=4m, ξ=0.05

Para empezar, calcularemos la rigidez equivalente del sistema, el sistema tiene 2
columnas de 1 tipo y 2 riostas y el sistema es paralelo
Rigidez columna
3 4
12 EI 21538.1∗10 ∗0.25
kc= 3
= 3
=3116.04456 KN /m
H 3
Rigidez riosta
As∗Es 2
kriosta= ∗cos θ
Lriosta

Lriosta=√ H + L =√ 3 + 4 =5
2 2 2 2

2
π∗d 2 π∗(0.0254)
As= = =5.07∗10− 4 m 2
4 4
−1 3
θ=tan =38.87 º
4
−4 3
5.07∗10 ∗199947.96∗10 2
kriosta= ∗cos 38.87=12290.01 KN /m
5
Rigidez equivalente

kequi=2 ( kc )+ 2 ( kriosta ) =30812.11 KN /m

Para hallar la masa, usaremos la fórmula de frecuencia natural y beta

wn=
√ kequi ;
m
kequi=wn ∗m y También tenemos que β=
2 wf
wn
Para que la riosta alcance la fluencia, sabemos que la fluencia en la mayor fuerza que
actual en el material.
Por la formula, tenemos que 𝑢𝑔(𝑡) = 0.25 ∗ 𝑠𝑒𝑛(12 ∗ 𝜋𝑡) por lo que
Um =0.25m y wf = 12𝜋 rad/seg
Para hallar la mayor fuerza
Pmax = Kriosta * U = 12290.01 KN/m *0.25m =3072.5KN
Con esto podemos hallar el Umax
Umax =Uest *Rd
Donde
1 1
Rd= = =10
2ξ 2(0.05)
p 3072.5
Umax = ∗Rd= ∗10=2.49 m
k 12290.01
Aplicando esta formula
2
Umax=Um∗β ∗Rd
2 2
Um∗β ∗1 0.25∗β ∗1
Umax= →2.49=
√ ( 1−β ) +( 2 ξ β )
2 2 2
√( 1−β ) + ( 2∗0.05∗β )
2 2 2

Despejando se obtiene β=0.485

Con la ecuación de B encontramos wn

wf 12 π
β= → 0.485= → wn=77.73 rad /seg
wn wn
2 kequi 30812.11
kequi=wn ∗m → m= = → m=0.519 Kg
wn2 77.732∗9.81
P5. Presentando el algoritmo correspondiente y usando como carga de excitación
orientada en el mismo sentido de u: 𝑃(𝑡) = 100𝑘𝑁 ∗ sen(25 ∗ #𝑡) + 50𝑘𝑁 ∗ 𝑐𝑜𝑠(25 ∗ #𝑡)
grafique simultáneamente la respuesta del sistema SDOF mostrado en la figura P2.b,
para valores de amortiguamiento =0.02; =0.05; =0.1. Para ello use M=100*#Ton;
k1=10 MN/m y k2=20 MN/m; uo=5cm y vo=0.5 m/s. (Use un tiempo de evaluación de t=0
a t= 10seg).
import numpy as np
import matplotlib.pyplot as plt

# Constantes
M = 1200*1000 # kg
k1 = 10*10**6 # N/m
k2 = 20*10**6 # N/m
uo = 0.05 # m
vo = 0.5 # m/s
# Carga de excitación
def carga_excitacion(t):
return 100*10**3*np.sin(300*t) + 50*10**3*np.cos(300*t)

# Ecuación de movimiento
def ec_movimiento(xi, t):
wn = np.sqrt(k2/M) # Frecuencia natural
wd = wn*np.sqrt(1-xi**2) # Frecuencia amortiguada
c1 = uo
c2 = (vo + xi*wn*uo)/wd
return np.exp(-xi*wn*t)*(c1*np.cos(wd*t) + c2*np.sin(wd*t)) + (1/M)*np.exp(-
xi*wn*t)*np.cumsum([carga_excitacion(ti) for ti in t])

# Tiempo de evaluación
t = np.arange(0, 10, 0.01) # s
# Respuesta para diferentes valores de amortiguamiento
xis = [0.02, 0.05, 0.1]
for xi in xis:
u = ec_movimiento(xi, t)
plt.plot(t, u, label=f'ξ = {xi}')

# Gráfica
plt.legend()
plt.xlabel('Tiempo (s)')
plt.ylabel('Desplazamiento (m)')
plt.show()

P6. Presentando el algoritmo correspondiente, grafique para el sistema estructural del

punto P5 las siguientes cantidades en el tiempo 0 a 10 seg, usando =0.05:
Es: Energía Potencial Elástica
Ek: Energía Cinética
Ed: Energía Disipada
Ei: Energía de Entrada
B = Es+Ek+Ed-Ei : Balance de Energía
import numpy as np
import matplotlib.pyplot as plt
# Constantes
M = 1200*1000 # kg
k1 = 10*10**6 # N/m
k2 = 20*10**6 # N/m
uo = 0.05 # m
vo = 0.5 # m/s

# Carga de excitación
def carga_excitacion(t):
return 100*10**3*np.sin(300*t) + 50*10**3*np.cos(300*t)

# Ecuación de movimiento
def ec_movimiento(xi, t):
wn = np.sqrt(k2/M) # Frecuencia natural
wd = wn*np.sqrt(1-xi**2) # Frecuencia amortiguada
c1 = uo
c2 = (vo + xi*wn*uo)/wd
u = np.exp(-xi*wn*t)*(c1*np.cos(wd*t) + c2*np.sin(wd*t))
dudt = -xi*wn*np.exp(-xi*wn*t)*(c1*np.cos(wd*t) + c2*np.sin(wd*t)) - wn*np.exp(-
xi*wn*t)*(c1*np.sin(wd*t) - c2*np.cos(wd*t))
w = np.sqrt(k2/M)
ke = 0.5*M*dudt**2
pe = 0.5*k2*u**2
kd = 0.05*w*ke**0.5*u**0.5
ki = carga_excitacion(t)*u
bal_energia = pe + ke - kd - ki
return u, ke, pe, kd, ki, bal_energia

# Tiempo de evaluación
t = np.arange(0, 10, 0.01) # s

# Respuesta para diferentes valores de amortiguamiento

xis = [0.05]
for xi in xis:
u, ke, pe, kd, ki, bal_energia = ec_movimiento(xi, t)
plt.plot(t, u, label=f'ξ = {xi}')

plt.figure()
plt.plot(t, ke, label='Energía cinética')
plt.plot(t, pe, label='Energía potencial elástica')
plt.plot(t, kd, label='Energía disipada')
plt.plot(t, ki, label='Energía de entrada')
plt.plot(t, bal_energia, label='Balance de energía')
plt.legend()
plt.xlabel('Tiempo (s)')
plt.ylabel('Energía (J)')

# Gráfica
plt.show()
P7. El bloque mostrado es impactado por un objeto de masa 0.1 kg, obteniendo un
choque perfectamente inelástico. Si la velocidad del objeto es de 20 m/s, la rigidez del
resorte es de 1 kN/m, el amortiguamiento del sistema es de =0.05 y la masa es de
m=10*# kg, determine el desplazamiento, la velocidad y la aceleración resultante del
bloque 2 segundos después del impacto.

Datos de entrada
Mobjeto=0.1kg
Vobjeto=20m/s
Kresorte=1KN/m
ξ=0.05
Mbloque= 10*12=120kg
Para empezar, viendo el sistema aplicamos la ecuación de conservación de momentos

mo∗vo=( m+mo ) vf
Donde
Mobjeto = mo = 0.1 kg
Mbloque = m =120 kg
Sabemos que la velocidad inicial es vo = 20m/s por lo que calcularemos la velocidad final
después del impacto
mo∗vo 0.1∗20
vf = = =0.01665 m/ s
m+ mo 120+0.1
La velocidad después del impacto del objecto seria siendo la velocidad inicial del bloque
por lo que tenemos que

û ( 0 ) =0.01665 m/s
utilizando la ecuación de vibración libre amortiguado tenemos

−ξ ωt û ( o ) +wo∗u(o)
u ( t )=e (u(0)*cos(wo*t)+ ∗sen (wo∗t))
wo

Donde wo=wn∗ √ 1−ξ 2 y wn=

√ Kres
Msistema
Msistema=m+mo=120 kg +0.1 kg=120.1 kg

wn=
√ 1
120.1
=0.0912 rad /seg

wo=wn∗√ 1−ξ 2=0.0912∗√ 1−0.052=0.091 rad / seg

Y como U(0)=0 entonces

−ξ ωnt û (o)
u ( t )=e ( ∗sen(wo∗t )) (1)
wo
Cuando t=2 segundos
0.01665
u ( 2 seg )=e
−0.05∗0.0912∗2
( ∗sen ( 0.091∗2 ))→u ( 2 seg )=¿ 5.76*10^-4 m =0.5759 mm
0.091
u ( 2 seg )=0.5759 mm →desplazamiento final
Para calcular la velocidad final derivamos la ecuación (1) y obtendremos

û ( t )=−ξ ωn∗e−ξ ωnt ( ûwo( o ) ∗sen ( wo∗t))+e −ξ ωnt

(û ( o )∗cos (wo∗t)) (2)
cuando t=2seg y remplazando tenemos

û ( 2 seg )=0.016493 m/ seg → velocidad final

para aceleración derivamos la ecuación (2) y obtendremos

ü ( t )=ξ 2 wn2∗e−ξ ωnt ( ûwo( o) ∗sen ( wo∗t ) )−2∗ξ ωn∗e

−ξωnt
( û ( o )∗cos ( wo∗t ) )−e−ξωnt (û ( o )∗wo∗sen(wo∗t))

cuando t=2 segundos y remplazando tendremos

−8 −4 −3
ü ( t )=1.1975∗10 −1.504∗10 −1.501∗10
−3 2
ü ( t )=−1.6514∗10 m/seg → aceleración
P8. Para un sistema SDOF con masa M, rigidez K y amortiguamiento c, se conoce que,
al finalizar 2*# ciclos de vibración libre, este alcanza un desplazamiento de uo/3. Si la
masa del sistema se disminuye en un 20%, determine también en términos de uo, el
nuevo desplazamiento que se alcanzara al completarse los mismos 2*# ciclos de
vibración libre. (Recuerde que uo representa el desplazamiento inicial del sistema).
Para resolver este problema, podemos utilizar la ecuación de desplazamiento de un
sistema SDOF amortiguado, que está dada por:
u(t) = uoe^(-ζωnt) * [cos(ωd*t) + (ζ/√(1-ζ^2))sin(ωdt)]
donde:
u(t) es el desplazamiento en el tiempo t.
uo es el desplazamiento inicial.
ζ es el coeficiente de amortiguamiento.
ωn es la frecuencia natural del sistema.
ωd es la frecuencia de amortiguamiento, dada por ωd = ωn * √(1-ζ^2).
Podemos encontrar la frecuencia natural del sistema y la frecuencia de amortiguamiento
utilizando las siguientes ecuaciones:
ωn = √(K/M) ωd = ωn * √(1-ζ^2)
Ahora, sabemos que al final de 24 ciclos de vibración libre, el sistema alcanza un
desplazamiento de uo/3. Podemos usar esta información para encontrar el coeficiente de
amortiguamiento ζ:
uo/3 = uoe^(-ζωn24T) * [cos(ωd*24T) + (ζ/√(1-ζ^2))sin(ωd24T)]
donde T es el período de vibración del sistema, dado por T = 2π/ωn.
Despejando ζ de esta ecuación, obtenemos:
ζ = -ln(uo/3/u0)/(ωn*24T)
Ahora podemos encontrar la nueva frecuencia natural y la nueva frecuencia de
amortiguamiento utilizando las siguientes ecuaciones:
M' = 0.8M (masa del sistema reducida en un 20%) K' = K (la rigidez no cambia) c' = c (el
amortiguamiento no cambia)
ωn' = √(K'/M') = √(K/M') ωd' = ωn' * √(1-ζ^2)
Finalmente, podemos utilizar la ecuación de desplazamiento con las nuevas frecuencias y
el mismo tiempo (24 ciclos de vibración libre) para encontrar el nuevo desplazamiento u' al
final de los 24 ciclos:
u' = uo*e^(-ζωn'*24T) * [cos(ωd'*24T) + (ζ/√(1-ζ^2))*sin(ωd'*24T)]
Reemplazando las expresiones para ζ, ωn', ωd' y T, obtenemos:
u' = (uo/3)e^(ln(uo/3/uo)24/(0.8√(K/M))) * [cos(√(K/M)24√(1-ln(uo/3/uo)^2)/(0.8√(K/M))) + (-
ln(uo/3/uo)/√(1-ln(uo/3/uo)^2))sin(√(K/M)24√(1-ln(uo/3/uo)^2)/(0.8√(K/M)))]
Simplificando esta expresión, obtenemos:
u' = (uo/3) e^(-30ln(3)/√(K/M)) * [cos(24√(K/M)(1+ln(3))/(0.8√(K/M))) - (ln(3)/√(1+
(ln(3))^2))sin(24√(K/M)(1+ln(3))/(0.8√(K/M)))]

P9. Usando el método de condensación estática se desea obtener un modelo dinámico de 1 G.D.L
equivalente para la estructura mostrada. Para esto primero determine la matriz de rigidez de la
armadura y posteriormente encuentre la rigidez equivalente horizontal (Kl) asociada a “u” si # es
par o vertical (Kv) asociada a “v” si # es impar (Utilice H=2m; B=1m; Acero A572 con Es=200 GPa;
Todas las conexiones Articuladas)

Para la solución de esta tarea se ha preparado un documento en formato Excel que

contiene el procedimiento detallado y los valores proporcionados para llevar a cabo el
análisis requerido. Adjuntamos este archivo con el fin de facilitar la revisión y comprensión
de la metodología empleada en la resolución de la tarea. Por favor, revise el archivo
adjunto para obtener una visión completa del proceso seguido y los resultados obtenidos.

P10. Verifique los resultados obtenidos en el punto P9 usando SAP 2000. Presente la
metodología usada y los valores obtenidos de Kl o Kv según corresponda a su caso.
(Sugerencia: En SAP 2000 aplique una carga unitaria en la dirección de interés para
 determinar la rigidez equivalente para un sistema de 1. G.D.L). Debe presentar el modelo
digital utilizado.
Este punto se realizó por medio de sap2000, se hizo el modelo, se definió el material y las
secciones de las columnas, vigas y riostas que quedo así:

Se le aplicó una carga unitaria en el punto de interés para buscar los grados de libertad y arrojaron
los siguientes datos

TABLE: Joint Displacements

Join
t OutputCase CaseType StepType StepNum U1 U2 U3 R1 R2 R3
Text Text Text Text Unitless cm cm cm Radians Radians Radians
1 DEAD LinStatic 1.561657 0 0.071 0 0.000354 0

kcond 0.640345
Archivos adjuntos

En la sección de archivos y deducciones solicitadas de este documento, se incluyen los archivos

necesarios para la implementación y ejecución de los algoritmos utilizados en la resolución de los
problemas planteados. En este sentido, se adjunta un archivo en formato Excel que contiene los
datos y fórmulas utilizados para la solución del punto 9, el cual ha sido diseñado de manera
rigurosa y eficiente para garantizar la precisión y eficacia en la solución del problema. Asimismo, se
incluye un documento en formato drive que contiene los algoritmos utilizados en la resolución de
los diferentes problemas planteados, lo cual permite su fácil acceso y modificación en caso de ser
necesario. Además, se adjunta una carpeta con la solución del punto 10 en el programa SAPP
2000, el cual es un software especializado utilizado en el análisis y diseño de estructuras. La
inclusión de estos archivos es esencial en el contexto de la investigación científica y la ingeniería,
ya que permite la reproducibilidad y validación de los resultados obtenidos, y facilita el trabajo de
otros investigadores en la implementación y extensión de los algoritmos y soluciones propuestas.