Konstrukcijos tikrinių ir priverstinių virpesių tyrimas

3. Tirkinių kampinių dažnių ir tikrinių formų apskaičiavimas MATLAB programa.
%Mechanines sistemos neslopinamu tikriniu virpesiu
%apskaiciavimas MATLAB programine iranga
clc
close all
clear all
K=[];
M=[];
omega=[];
y=[];
%Konstrukcijos laisves laipsniuskaicius
N=6;
%Spyruokliu standumai, Nm
k1=10000;
k2=6000;
k3=7000;
k4=6000;
k5=7000;
k6=7000;
%Mases, kg
m1=35;
m2=35;
m3=11;
%Atstumai, m
x1=0.48;
x2=0.16;
x3=-0.48;
x4=0.32;
x5=-0.32;
x6=-0.32;
x7=0.32;
x8=-0.32;
x9=0.32;
x10=-0.32;
%Inercijos momentai, kg*m^2
I1=50;
I2=46;
I3=13;
%Standumo matrica
K=[k1+k2+k3 k1*x1+k2*x2+k3*x3 -k2-k3 -k2*x7-k3*x8 0 0;
k1*x1+k2*x2+k3*x3 k1*x1^2+k2*x2^2+k3*x3^2 -k2*x2-k3*x3 -k2*x2*x7-k3*x3*x8 0 0;
-k2-k3 -k2*x2-k3*x3 k2+k3+k4+k5 k2*x7+k3 +k4*x4+k5*x5 -k4-k5 -k4*x9-k5*x10;
-k2*x7-k3*x8 -k2*x2*x7-k3*x3*x8 k2*x7+k3*x8+k4*x4+k5*x5 k2*x7^2+k3*x8^2+k4*x4^2+ 2 -k4*x4-k5*x5 -k4*x4*x9-k5*x5*x10;
0 0 -k4-k5 -k4*x4-k5*x5 k4+k5+k6 k4*x9+k5*x10+k6*x6;
0 0 -k4*x9-k5*x10 -k4*x4*x9-k5*x5*x10 k4*x9+k5*x10+k6*x6 k4*x9^2+k5* x6^2];
%Masiu matrica
M=[m1 0 0 0 0 0;
0 I1 0 0 0 0;
0 0 m2 0 0 0;
0 0 0 I2 0 0;
0 0 0 0 m3 0;
0 0 0 0 0 I3];
disp('Laisves laipsniu skaicius:');
N
disp('Standumo matrica:');
K
disp('Masiu matrica:');
M
[Y,omega]=eig(K,M);
omega=diag(sqrt(omega));
disp('Savieji kampiniai dazniai:');
omega
%Tikriniu vektoriu normavimas
F=Y'*M*Y;
f=diag(F);
f=sqrt(f);
for i=1:N
Y(1:N,i)=Y(1:N,i)/f(i);
end;
disp('Normuoti savieji vektoriai:');
Y
Matlab komandiniame lange gauti tokie rezultatai:
Laisves laipsniu skaicius:
N =
6
Standumo matrica:
K =
1.0e+004 *
2.3000 0.2400 -1.3000 0.0320 0 0
0.2400 0.4070 0.2400 -0.1382 0 0
-1.3000 0.2400 2.6000 -0.0640 -1.3000 0.0320
0.0320 -0.1382 -0.0640 0.2662 0.0320 -0.1331
0 0 -1.3000 0.0320 2.0000 -0.2560
0 0 0.0320 -0.1331 -0.2560 0.2048
Masiu matrica:
M =
35 0 0 0 0 0
0 50 0 0 0 0
0 0 35 0 0 0
0 0 0 46 0 0
0 0 0 0 11 0
0 0 0 0 0 13
Savieji kampiniai dazniai:
omega =
4.1759
6.1083
12.3786
15.6865
29.8256
46.5982
Normuoti savieji vektoriai:
Y =
0.0181 0.0481 -0.0132 -0.0853 0.1347 -0.0182
-0.0743 -0.0914 -0.0563 -0.0542 0.0030 0.0013
0.0147 0.0646 -0.0265 -0.1045 -0.0837 0.0744
-0.1150 0.0499 0.0773 -0.0068 0.0007 -0.0018
0.0004 0.0584 -0.0478 -0.0614 -0.1005 -0.2671
-0.0861 0.1249 -0.1982 0.1153 0.0241 0.0271
Žinant tikrinius normuotus vektorius galime pavaizduoti grafiškai tikrąsias sistemos formas.
4. Konstrukcijos amplitudės ir fazės dažninės charakteristikos apskaičiavimas.
%Programa daugelio masiu sistemos priverstiniu virpesiu
%ADCH ir FDCH apskaiciuoti
clc
close all
clear all
K=[];
M=[];
C=[];
FS=[];
FC=[];
w=[];
A=[];
B=[];
U=[];
UC=[];
US=[];
AMPL=[];
PHASE=[];
%Konstrukcijos laisves laipsniu skaicius
N=6;
%Spyruokliu standumai, Nm
k1=10000;
k2=6000;
k3=7000;
k4=6000;
k5=7000;
k6=7000;
%Mases, kg
m1=35;
m2=35;
m3=11;
%Atstumai, m
x1=0.48;
x2=0.16;
x3=-0.48;
x4=0.32;
x5=-0.32;
x6=-0.32;
x7=0.32;
x8=-0.32;
x9=0.32;
x10=-0.32;
%Inercijos momentai, kg*m^2
I1=50;
I2=46;
I3=13;
%Maziausias savasis kampinis daznis
omega1=4.1759;
%Dinamiskumo koeficientas
miu=12;
%Standumo matrica
K=[k1+k2+k3 k1*x1+k2*x2+k3*x3 -k2-k3 -k2*x7-k3*x8 0 0;
k1*x1+k2*x2+k3*x3 k1*x1^2+k2*x2^2+k3*x3^2 -k2*x2-k3*x3 -k2*x2*x7-k3*x3*x8 0 0;
-k2-k3 -k2*x2-k3*x3 k2+k3+k4+k5 k2*x7+k3*x8+k4*x4+k5*x5 -k4-k5 -k4*x9-k5*x10;
-k2*x7-k3*x8 -k2*x2*x7-k3*x3*x8 k2*x7+k3*x8+k4*x4+k5*x5 k2*x7^2+k3*x8^2+k4*x4^2+k5*x5^2 -k4*x4-k5*x5 -k4*x4*x9-k5*x5*x10;
0 0 -k4-k5...

Daugiau informacijos...

★ Klientai rekomenduoja

Šį rašto darbą rekomenduoja mūsų klientai. Ką tai reiškia?

Mūsų svetainėje pateikiama dešimtys tūkstančių skirtingų rašto darbų, kuriuos įkėlė daugybė moksleivių ir studentų su skirtingais gabumais. Būtent šis rašto darbas yra patikrintas specialistų ir rekomenduojamas kitų klientų, kurie po atsisiuntimo įvertino šį mokslo darbą teigiamai. Todėl galite būti tikri, kad šis pasirinkimas geriausias!