Scribd
Upload
21 views6 pages

Lab5 Matlab

The document contains code for calculating various properties of mechanical vibrations, including: 1) Finding the roots of a polynomial coefficient matrix to determine natural frequencies of vibration. 2) Deriving equations for angular velocity, linear and angular acceleration as functions of time using derivatives. 3) Evaluating the equations numerically by plugging in values for time. 4) Repeating the process for a second example involving sinusoidal vibration.

Uploaded by

Viorel Mirzenco
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
21 views6 pages

Lab5 Matlab

The document contains code for calculating various properties of mechanical vibrations, including: 1) Finding the roots of a polynomial coefficient matrix to determine natural frequencies of vibration. 2) Deriving equations for angular velocity, linear and angular acceleration as functions of time using derivatives. 3) Evaluating the equations numerically by plugging in values for time. 4) Repeating the process for a second example involving sinusoidal vibration.

Uploaded by

Viorel Mirzenco
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd

1a

>> Coef=[1,1.5,0.75,-125*pi/180];
>> t=roots(Coef)

t=

-1.1606 + 1.1443i
-1.1606 - 1.1443i
0.8213 + 0.0000i
1b
>> syms t;
>> fe=t^3+1.5*t^2+0.75*t;
>> omega=diff(fe)

omega =

3*t^2 + 3*t + 3/4

>> R=25;
>> OM=(pi*R)/8;
>> OO1=(2*pi*R)/4;
>> MO1=OO1-OM

MO1 =

29.4524

>> Vm=omega*MO1

Vm =

(75*pi*(3*t^2 + 3*t + 3/4))/8

>> An=omega^2*MO1

An =

(75*pi*(3*t^2 + 3*t + 3/4)^2)/8

>> epsilon=diff(omega)

epsilon =

6*t + 3

>> At=epsilon*MO1

At =
(75*pi*(6*t + 3))/8

>> Am=sqrt(An^2+At^2)

Am =

((5625*pi^2*(6*t + 3)^2)/64 + (5625*pi^2*(3*t^2 + 3*t + 3/4)^4)/64)^(1/2)

>> t=0.8213;
>> omega =3*t^2 + 3*t + 3/4

omega =

5.2375

>> Vm =(75*pi*(3*t^2 + 3*t + 3/4))/8

Vm =

154.2571

>> An =(75*pi*(3*t^2 + 3*t + 3/4)^2)/8

An =

807.9219

>> epsilon =6*t + 3

epsilon =

7.9278

>> At =(75*pi*(6*t + 3))/8

At =

233.4930

>> Am =((5625*pi^2*(6*t + 3)^2)/64 + (5625*pi^2*(3*t^2 + 3*t + 3/4)^4)/64)^(1/2)

Am =

840.9856
1c
at
am
Vm

an
ɛ
ω

2a
Fail functia
function[fe]=functie(t);
fe=16*(sin(pi*t))^2-125*pi/180;

Programul
t=fzero('functie',pi/2)

t=

1.8796
2b
>> syms t;
>> fe=16*(sin(pi*t))^2;
>> omega=diff(fe)

omega =

32*pi*cos(pi*t)*sin(pi*t)

>> R=25;
>> OO1=(2*pi*R)/4

OO1 =

39.2699

>> Vo=omega*OO1

Vo =

400*pi^2*cos(pi*t)*sin(pi*t)

>> An=omega^2*OO1

An =

12800*pi^3*cos(pi*t)^2*sin(pi*t)^2

>> epsilon=diff(omega)

epsilon =

32*pi^2*cos(pi*t)^2 - 32*pi^2*sin(pi*t)^2

>> At=epsilon*OO1

At =

-(25*pi*(32*pi^2*sin(pi*t)^2 - 32*pi^2*cos(pi*t)^2))/2

>> Ao=sqrt(An^2+At^2)

Ao =

((625*pi^2*(32*pi^2*sin(pi*t)^2 - 32*pi^2*cos(pi*t)^2)^2)/4 +
163840000*pi^6*cos(pi*t)^4*sin(pi*t)^4)^(1/2)

>> t=1.8796;
>> omega =32*pi*cos(pi*t)*sin(pi*t)
omega =

-34.5011

>> Vo =400*pi^2*cos(pi*t)*sin(pi*t)

Vo =

-1.3549e+03

>> An =12800*pi^3*cos(pi*t)^2*sin(pi*t)^2

An =

4.6744e+04

>> epsilon =32*pi^2*cos(pi*t)^2 - 32*pi^2*sin(pi*t)^2

epsilon =

229.6841

>> At =-(25*pi*(32*pi^2*sin(pi*t)^2 - 32*pi^2*cos(pi*t)^2))/2

At =

9.0197e+03

>> Ao =((625*pi^2*(32*pi^2*sin(pi*t)^2 - 32*pi^2*cos(pi*t)^2)^2)/4


+ 163840000*pi^6*cos(pi*t)^4*sin(pi*t)^4)^(1/2)

Ao =

4.7606e+04
2c
at

ao
Vo
an
ɛ
ω

You might also like