mitcil+bu 0 =
- mit 2mzwhi +wn" m.u 0.
=
solucion: M = ext
i dext
=
i x ext
=
Next 2 [Link]+
+
wr2. 0.
=
x+2zwn.x whe +
0.
=
x -
=
2swn=-ton" awnie
2
wn2-1
-
xi -
=
bwn +
12 = -
gwn-wn
①
sobreamortiguado(((1 =
car)
Xe,2 = -
wh wh
= --
5 1 -
e = Aeit -> Belzt
ult) =
=
-
e(- wn + wnr- 1 t
-> Be
(wn-wnt
u(0) A B
= +
Mo
=
B
=
No
= -
A
wntwart
i A.(zwn wn) el-
= +
t
13 ( -swn-wn) el wn-warst
i(0) f(xyn
= +
wnt) A(a(wn wns) -
-
+
mol-swn-wn ) no =
=2awnrg mo(-zwn-wn Fan) +
=io
A i -uo(-zwn-rn
=
zwnt Wo-A
li-ew(-zwn-wns. e(-wn
+wn) ↳ (-wn-wn)
I I I
-
7
Ii-rot-wn-wneI
I
uit) =
I
coen
-
C
->
zwnt zwnet
& critico
amortiguamiento (3 1)
=
X1,2 wh
=
- z
=>u(t) Ae
=
-
4wnt te- swnt
+
u(0) A
=
No
=
i(t) = -
mozwnewht_wnte-hwht-Behwnt
ind Bl-zwnte-
=
- wntteswnt)-mogwon
S in
=
B i
=
wwn
+
wnt
=
Mit) Noe wnt+ Lin
=
+mozwn) te
③
subamortiguado(3)
xeiz =
-zwn (wi
=
WD
--
zwn I WDi
uct) xe(wn wiltbtzwn -wbilt
=
+
eit =
Los(t) isenct) +
uit) wnt) A(0s(wit) +isen(wit)) B(cos(wit) isen(-wit))
-
+
e
+
=
e(t) 4wnt) Acos(wit) Bren(wit))
-
e
= +
u(0) A
=
no
=
eit) 4wnt(Acos(wit) Ben(wit))
-
wn e
= +
- +
e
-
zwht) -
Awasen (wist) +
BWDcos (wit))
-
lico) =
-
wn
- e
-
1
-
4wnt(uncos(wit) Ben(wit)) +
+ o
/
e-zwnt)- - (wist)
mowissen (wit)) +
BWDC)
--
wn Mo+ is
BaD=
B i + 3 wnwo
=
--
WD
=rect)
eqwnt((rocos(wost) (3 cen(wot
>
+
Graficando los 3 cosos en Mathab:
() =
en
-
05
en
(5) -30
2.
=
Ve-
-
· =3
=0,004674
=>
TD 1 0,132 el periodo depende
=
->
=
30 de Los ciclos
se
*
asume io=in (V 1) =
wn 2
==
47,124/r)
0,133
k =
wh2 m
47,1272.15
=
=33309,9/N/m2)
1
2mwn2
=
2.15-47,124.0,004674
=>
/No
-
3.5
Problema 3
Encuentre la ecuación de desplazamiento y velocidad de un sistema de 1GDL mostrado en Figura 1. Considere
masa igual a 30 [kg], rigideces k1 = 10000 [N/m], k2 = 8000 [N/m], k3 = 5000 [N/m] y fuerza armónica P(t) =
50 sin(50t). Suponga que para el instante t=5[s], la posición es de -0.2 [m] y la velocidad es 5 [m/s]. Desprecie
amortiguamiento.
Indique:
• Frecuencia natural del sistema.
• Posición y velocidad inicial.
• Posición y velocidad máxima.
• ¿Qué pasaría si P(t) = 50 cos(50t)? ¿Cómo cambiarían los resultados?
Sol.
%Datos iniciales
m=30; %kg
k1=10000; %N/m
k2=8000; %N/m
k3=5000; %N/m
p0=50; %N
w=50; %rad/s
En primer lugar, se calcula la rigidez equivalente del sistema:
keq=k1+k2*k3/(k2+k3) %N/m
keq = 1.3077e+04
Luego, hacemos sumatoria de fuerzas, se obtiene la ecuación del sistema:
1
Sabiendo que
wn=sqrt(keq/m) %rad/s
wn = 20.8782
Despejamos el sistema para las condiciones de borde:
t0=5;
syms A B
eq1= A*cos(wn*t0)+B*sin(wn*t0)+p0/keq/(1-(w/wn)^2)*sin(t0*w)==-0.2;
eq2= -A*wn*sin(wn*t0)+B*wn*cos(wn*t0)+p0/keq*w/(1-(w/wn)^2)*cos(w*t0)==5;
sol = vpasolve([eq1, eq2], [A, B])
sol = struct with fields:
A: 0.30908455157202815007148654543092
B: -0.048550771785708385660618767769552
Reemplazamos y graficamos
t=0:0.01:5;
%desplazamiento
u=0.30908455157202815007148654543092*cos(wn*t)-0.0485507717857083856606187677
69552*sin(wn*t)+p0/keq/(1-(w/wn)^2)*sin(t*w);
umax=max(u)
umax = 0.3136
umin=min(u)
umin = -0.3136
plot(t,u,0,0.30908455157202815007148654543092,'o')
yline([umax umin],'--',{'Max','Min'})
legend('Desplazamiento en el tiempo','Desplazamiento inicial = 0.309
[m]','Desplazamiento máximo = 0.314 [m]')
xlabel('Tiempo [s]')
ylabel('Desplazamiento [m]')
grid on
2
%velocidad
v=-0.30908455157202815007148654543092*wn*sin(wn*t)-0.048550771785708385660618
767769552*wn*cos(wn*t)+p0/keq*w/(1-(w/wn)^2)*cos(w*t);
v0=-0.048550771785708385660618767769552*wn+p0/keq*w/(1-(w/wn)^2)
v0 = -1.0540
vmax=max(v)
vmax = 6.5721
vmin=min(v)
vmin = -6.5684
plot(t,v,0,v0,'o')
yline([vmax vmin],'--',{'Max','Min'})
xlabel('Tiempo [s]')
ylabel('Velocidad [m/s]')
legend('Velocidad en el tiempo','Velocidad inicial = -1.054 [m/
s]','Velocidad máxima = 6.572 [m/s]')
grid on
3
Luego, si se aplicara una fuerza cosenoidal:
syms C D
eq3= C*cos(wn*t0)+D*sin(wn*t0)+p0/keq/(1-(w/wn)^2)*cos(t0*w)==-0.2;
eq4= -C*wn*sin(wn*t0)+D*wn*cos(wn*t0)+p0/keq*w/(1-(w/wn)^2)*cos(w*t0)==5;
sol2 = vpasolve([eq3, eq4], [C, D])
sol2 = struct with fields:
C: 0.30834796037054792232233867644716
D: -0.049194505239525763576913951181234
%Desplazamiento
u2=0.30834796037054792232233867644716*cos(wn*t)-0.049194505239525763576913951
181234*sin(wn*t)+p0/keq/(1-(w/wn)^2)*cos(t*w);
4
u2_0=0.30834796037054792232233867644716+p0/keq/(1-(w/wn)^2)
u2_0 = 0.3075
u2max=max(u2)
u2max = 0.3130
u2min=min(u2)
u2min = -0.3129
%Velocidad
v2=-0.30834796037054792232233867644716*wn*sin(wn*t)-0.04919450523952576357691
3951181234*wn*cos(wn*t)-p0/keq*w/(1-(w/wn)^2)*sin(w*t);
v2_0=-0.049194505239525763576913951181234*wn
v2_0 = -1.0271
v2max=max(v2)
v2max = 6.5569
v2min=min(v2)
v2min = -6.5584
En base a estos resultados, podemos observar que tanto el desplazamiento como velocidad inicial disminuyen
un 0.5% y 2.6%, respectivamente. Por otro lado, se observa que sus valores máximos disminuyen en un 0.2%
aproximadamente, lo que podría deberse a solo un error asociado a los decimales.
abs(0.30908455157202815007148654543092-u2_0)/
0.30908455157202815007148654543092*100 %diferencia despl inicial porcentual
ans = 0.4996
abs(umax-u2max)/umax*100 %diferencia despl max porcentual
ans = 0.1966
abs((v0-v2_0)/v0)*100 %diferencia vel inicial porcentual
ans = 2.5552
abs(vmax-v2max)/vmax*100 %diferencia vel max porcentual
ans = 0.2315
