Ejemplo Sobreamortiguado y Amortiguacion Critica

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 15

UNIVERSIDAD NACIONAL DE SAN ANTONIO ABAD DEL

CUSCO

INGENIERÍA MECÁNICA
Trabajo: ejercicios sobre amortiguados y críticamente
amortiguados.
Curso: Análisis Vibracional
Docente: Arturo Macedo Silva
Alumno: 141537 Huisa Fernandez, Toly Kuper

Cusco – 2021-I

Sistema sobre amortiguado


b). 3 ẍ+ 12 ẋ +9 x=0
solución:
De la siguiente ecuación general se tiene: m ẍ + c ẋ+ k x=0
c 12 2 3
ζ= = = √ =1.1547>1 por lo tanto es un caso sobreamortiguado
2 √ mk 2 √ 3∗9 3
2
2
√ ζ −1=
√( 2 √3
3
−1= √
) 3
3

k 9
ω n=
√ √ m
= =√ 3
3

x0 ω n ( ζ + √ ζ 2−1 ) + ẋ0 −x 0 ωn ( ζ −√ ζ 2 −1 )− ẋ 0
x ( t )=
( 2 ωn √ ζ −1 2 ) 2

e (−ζ +√ ζ −1)ω t + n

( 2 ωn √ ζ −12 ) e ¿¿

( 2 √33 + √33 )+ ẋ ( 2 √3 3 − √33 )− ẋ


( ) ( )
( x 0)( √3) 0 −2 √3 √3 −( x 0 ) ( √3 ) 0 −2 √ 3 √ 3
( + ) √3 t ( − )√3 t
3 3 3 3
x ( t )= e + e
3 3
2( √ 3)( √ ) 2( √ 3)( √ )
3 3

3 x 0 + ẋ 0 −t −x 0− ẋ 0 −3 t
x ( t )= [ 2 ] [e +
2
e ]

b.1) para: x 0=0 ẋ 0=1


3(0)+(1) −t −(0)−(1) −3 t
x ( t )= [ 2 ] [
e +
2
e ]
Respuesta en el tiempo: x ( t )=0.5 e−t −0.5 e−3 t
Procedemos a comprobar si nuestra ecuación en el tiempo es correcta tanto para la posición
inicial como para la velocidad inicial.

x ( t=0 )=0.5 e−(0)−0.5 e−3 (0 )


x ( t=0 )=x 0=0.5−0.5=0

x 0=0

−1 −t 3 −3 t
ẋ ( t )= e + e
2 2
−1 −(0 ) 3 −3(0)
ẋ ( t=0 )= e + e
2 2
−1 3
ẋ ( t=0 )= + =1
2 2
ẋ 0=1
clear
clc
m = 3;%masa (kg)
c = 12;%constante de amortiguamiento (N.s/m)
k = 9;%constante de resorte (N/m)
wn = sqrt(k/m);%frecuebcia natural (1/s)
cc = 2*m*wn;
pfi = c/cc;%factor de amortiguamiento
M =(sqrt((pfi^2)-1));
xo = 0;%posicion inicial
vo = 1;%velocidad inicial
t = 0:0.01:10;%intervalo de tiempo
C1 = (xo.*wn.*(pfi+M)+vo)/(2.*wn.*M);
C2 = (-xo.*wn*(pfi-M)-vo)/(2.*wn.*M);
xt = C1.*exp((-pfi+M).*wn.*t)+C2.*exp((-pfi-M).*wn.*t);
figure(1)
plot(t,xt)
grid

>> wn

wn =

1.7321
>> pfi

pfi =

1.1547

b.2) para: x 0=0 ẋ 0=0


3(0)+(0) −t −(0)−(0) −3t
x ( t )= [ 2 ] [
e +
2 ]e

Respuesta en el tiempo: x ( t )=0


clear
clc
m = 3;%masa (kg)
c = 12;%constante de amortiguamiento (N.s/m)
k = 9;%constante de resorte (N/m)
wn = sqrt(k/m);%frecuebcia natural (1/s)
cc = 2*m*wn;
pfi = c/cc;%factor de amortiguamiento
M =(sqrt((pfi^2)-1));
xo = 0;%posicion inicial
vo = 0;%velocidad inicial
t = 0:0.01:10;%intervalo de tiempo
C1 = (xo.*wn.*(pfi+M)+vo)/(2.*wn.*M);
C2 = (-xo.*wn*(pfi-M)-vo)/(2.*wn.*M);
xt = C1.*exp((-pfi+M).*wn.*t)+C2.*exp((-pfi-M).*wn.*t);
figure(1)
plot(t,xt)
grid
b.3) para: x 0=1 ẋ 0=¿−1

3( 1)+(−1) −t −( 1)−(−1) −3 t
x ( t )= [ 2 ] [
e +
2 ]
e

Respuesta en el tiempo: x ( t )=e−t


Procedemos a comprobar si nuestra ecuación en el tiempo es correcta tanto para la posición
inicial como para la velocidad inicial.

x ( t=0 )=e−(0)
x ( t=0 )=x 0=1

x 0=1

ẋ ( t )=−e−t

ẋ ( t=0 )=−e−(0)
ẋ ( t=0 )=−1
ẋ 0=−1

La ecuación cumple con las condiciones iniciales


clear
clc
m = 3;%masa (kg)
c = 12;%constante de amortiguamiento (N.s/m)
k = 9;%constante de resorte (N/m)
wn = sqrt(k/m);%frecuebcia natural (1/s)
cc = 2*m*wn;
pfi = c/cc;%factor de amortiguamiento
M =(sqrt((pfi^2)-1));
xo = 1;%posicion inicial
vo = -1;%velocidad inicial
t = 0:0.01:10;%intervalo de tiempo
C1 = (xo.*wn.*(pfi+M)+vo)/(2.*wn.*M);
C2 = (-xo.*wn*(pfi-M)-vo)/(2.*wn.*M);
xt = C1.*exp((-pfi+M).*wn.*t)+C2.*exp((-pfi-M).*wn.*t);
figure(1)
plot(t,xt)
grid
Sistema críticamente amortiguado:
c.) 2 ẍ+ 8 ẋ +8 x=0

solución:
De la siguiente ecuación general se tiene: m ẍ + c ẋ+ k x=0
c 8
ζ= = =1 por lo tanto es un caso criticamente amortiguado
2 √ mk 2 √ 2∗8

x ( t )=[ x 0 + ( ẋ 0 +ω n x 0 ) t ] e−ω t
n

8
ω n=
√ 2
=2

c.1) para: x 0=0 ẋ 0=1

−2 t
x ( t )=[ (0)+ ((1)+(2)(0) ) t ] e

Respuesta en el tiempo: x ( t )=t e−2 t


Procedemos a comprobar si nuestra ecuación en el tiempo es correcta tanto para la posición
inicial como para la velocidad inicial.

x ( t=0 )=( 0 ) e−2 (0 )=0


x 0=0

ẋ ( t )=e−2 t −2 t e−2 t

ẋ ( t=0 )=e−2 (0 )−2 ( 0 ) e−2( 0)


ẋ 0=1

La ecuación cumple con las condiciones iniciales


clear
clc
m = 2;%masa (kg)
c = 8;%constante de amortiguamiento (N.s/m)
k = 8;%constante de resorte (N/m)
wn = sqrt(k/m);%frecuebcia natural (1/s)
cc = 2*m*wn;
pfi = c/cc;%factor de amortiguamiento
xo = 0;%posicion inicial
vo = 1;%velocidad inicial
t = 0:0.01:10;%intervalo de tiempo
M = (xo+(vo+wn*xo)*t);
xt =M.*(exp(-wn*t));
figure(1)
plot(t,xt)
grid
c.2) para: x 0=0 ẋ 0=0

x ( t )=[ (0)+ (( 0)+(2)(0) ) t ] e−2t

Respuesta en el tiempo: x ( t )=0


clear
clc
m = 2;%masa (kg)
c = 8;%constante de amortiguamiento (N.s/m)
k = 8;%constante de resorte (N/m)
wn = sqrt(k/m);%frecuebcia natural (1/s)
cc = 2*m*wn;
pfi = c/cc;%factor de amortiguamiento
xo = 0;%posicion inicial
vo = 0;%velocidad inicial
t = 0:0.01:10;%intervalo de tiempo
M = (xo+(vo+wn*xo)*t);
xt =M.*(exp(-wn*t));
figure(1)
plot(t,xt)
grid
c.3) para: x 0=1 ẋ 0=−1
−2t
x ( t )=[ (1)+ ( (−1)+(2)(1)) t ] e

Respuesta en el tiempo: x ( t )=(1+t)e−2 t


Procedemos a comprobar si nuestra ecuación en el tiempo es correcta tanto para la posición
inicial como para la velocidad inicial.

x ( t=0 )=( 1+0 ) e−2 (0 )=1


x 0=1

ẋ ( t )=−e−2t −2 t e−2 t

ẋ ( t=0 )=−e−2 (0)−2(0)e−2(0)


ẋ 0=−1
clear
clc
m = 2;%masa (kg)
c = 8;%constante de amortiguamiento (N.s/m)
k = 8;%constante de resorte (N/m)
wn = sqrt(k/m);%frecuebcia natural (1/s)
cc = 2*m*wn;
pfi = c/cc;%factor de amortiguamiento
xo = 1;%posicion inicial
vo = -1;%velocidad inicial
t = 0:0.01:10;%intervalo de tiempo
M = (xo+(vo+wn*xo)*t);
xt =M.*(exp(-wn*t));
figure(1)
plot(t,xt)
grid

Comentario:
Se logró graficar las gráficas para las ecuaciones, cabe mencionar que se comprobó que las
condiciones iniciales fueran las correctas a la hora de reemplazar t=0 para las ecuaciones
de posicionamiento y velocidad, también se comprobó mediante las gráficas de la misma
manera. También mencionar que para un xo = 0 y vo = 0 las gráficas arrojaron una línea
horizontal y esto tiene mucho sentido ya que la masa o cuerpo está en reposo.

Resolución de los ejercicios de forma analítica en manuscrito.

También podría gustarte