Scribd
Upload
20 views

Aero Elasticity

This MATLAB code defines variables and calculates matrices related to modeling the dynamics of a vehicle. It defines parameters like mass and coefficients, calculates inertia matrices, and determines the eigenvalues of the system dynamics matrix Q for increasing vehicle velocities. The eigenvalues become more negative and complex as the velocity increases, indicating greater stability.

Uploaded by

AkinoJohnkennedy
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
20 views

Aero Elasticity

This MATLAB code defines variables and calculates matrices related to modeling the dynamics of a vehicle. It defines parameters like mass and coefficients, calculates inertia matrices, and determines the eigenvalues of the system dynamics matrix Q for increasing vehicle velocities. The eigenvalues become more negative and complex as the velocity increases, indicating greater stability.

Uploaded by

AkinoJohnkennedy
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 7

1.

MATLAB CODE:
clear all; clc;
m = 100;
c = 1.7;
s = 7;
xcm = 0.5*c;
xf = 0.5*c;
e = xf/c-0.25;
K_k = 1.2277e7;
K_t = 7.2758e5;
Mthetadot = -1.2;
roll = 1.225;
aw = 2*pi;
M = (m*c^2/(2*xcm))-m*c;
b=c/2;
a11=(m*s^3*c)/3 + M*s^3/3;
a12 = ((m*s^2)/2)*(c^2/2 - c*xf) - M*s^2*xf/2;
a21 = a12;
a22= m*s*(c^3/3 - c^2*xf + xf^2*c) + M*(s*xf^2);
A=[a11,a12;a21,a22];
b11=c*s^3*aw/6;
b12=0;
b21=-e*c^2*s^2*aw/4;
b22=-c^3*s/8*Mthetadot;
B=[b11 b12; b21 b22]
c11=0;
c12=c*s^2*aw/4;
c21=c11;
c22=-e*c^2*s*aw/2;
C=[c11 c12; c21 c22]
D=zeros(2)
e11=K_k;
e12=0;
e21=e12;
e22=K_t;
E=[e11 e12; e21 e22]
V=32;
Q=[zeros(2) eye(2); -A\(roll*V^2*C+E) -A\(roll*V*B+D)];
lambda=eig(Q);
output:
AT v=32
B = 610.6209
0
-55.6101 5.1586
C = 0 130.8473
0 -15.8886
D=0
0
0
0
E =12277000
0
0
727580
lambda = -0.6332 +25.1244i
-0.6332 -25.1244i
-0.3353 +49.6901i
-0.3353 -49.6901i

Q=
1.0e+03 *
0
0 0.0010
0
0
0
0 0.0010
-0.6316 -0.0084 -0.0012 0.0000
0.0000 -2.4692 0.0076 -0.0007
>> [v,d]=eig(Q)
v=
-0.0010 - 0.0395i -0.0010 + 0.0395i 0.0000 - 0.0001i 0.0000 + 0.0001i
0.0041 + 0.0000i 0.0041 - 0.0000i -0.0001 - 0.0201i -0.0001 + 0.0201i
0.9938 + 0.0000i 0.9938 + 0.0000i 0.0046 + 0.0001i 0.0046 - 0.0001i
-0.0034 + 0.1033i -0.0034 - 0.1033i 0.9998 + 0.0000i 0.9998 + 0.0000i
d=
-0.6332 +25.1244i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -0.6332 -25.1244i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -0.3353 +49.6901i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.3353 -49.6901i
AT V=64,
lambda =
-1.3888
-1.3888
-0.5483
-0.5483
Q=

+25.0908i
-25.0908i
+47.5484i
-47.5484i

1.0e+03 *
0
0 0.0010
0
0
0
0 0.0010
-0.6316 -0.0338 -0.0025 0.0000
0.0000 -2.2606 0.0152 -0.0014
>> [v,d]=eig(Q)
v=
-0.0021 - 0.0387i -0.0021 + 0.0387i 0.0000 - 0.0004i 0.0000 + 0.0004i
0.0091 + 0.0002i 0.0091 - 0.0002i -0.0002 - 0.0210i -0.0002 + 0.0210i
0.9729 + 0.0000i 0.9729 + 0.0000i 0.0207 + 0.0008i 0.0207 - 0.0008i
-0.0174 + 0.2273i -0.0174 - 0.2273i 0.9996 + 0.0000i 0.9996 + 0.0000i
d=
-1.3888 +25.0908i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -1.3888 -25.0908i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -0.5483 +47.5484i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.5483 -47.5484i
AT V=96,
lambda =
-2.5133 +24.9684i
-2.5133 -24.9684i
-0.3923 +43.8003i

-0.3923 -43.8003i
Q=
1.0e+03 *
0
0 0.0010
0
0
0
0 0.0010
-0.6316 -0.0760 -0.0037 0.0000
0.0000 -1.9128 0.0228 -0.0021
>> [v,d]=eig(Q)
v=
-0.0036 - 0.0362i -0.0036 + 0.0362i -0.0001 + 0.0013i
0.0161 + 0.0009i 0.0161 - 0.0009i 0.0002 + 0.0228i
0.9136 + 0.0000i 0.9136 + 0.0000i -0.0583 - 0.0058i
-0.0631 + 0.3998i -0.0631 - 0.3998i -0.9980 + 0.0000i

-0.0001 - 0.0013i
0.0002 - 0.0228i
-0.0583 + 0.0058i
-0.9980 + 0.0000i

d=
-2.5133 +24.9684i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -2.5133 -24.9684i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -0.3923 +43.8003i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.3923 -43.8003i
AT V = 128,
lambda =
-4.7176 +24.2309i
-4.7176 -24.2309i
0.8434 +38.4369i
0.8434 -38.4369i

Q=
1.0e+03 *
0
0 0.0010
0
0
0
0 0.0010
-0.6316 -0.1351 -0.0049 0.0000
0.0000 -1.4260 0.0304 -0.0028
>> [v,d]=eig(Q)
v=
0.0058 + 0.0300i 0.0058 - 0.0300i
-0.0261 - 0.0049i -0.0261 + 0.0049i
-0.7537 + 0.0000i -0.7537 + 0.0000i
0.2428 - 0.6095i 0.2428 + 0.6095i

0.0012 - 0.0038i 0.0012 + 0.0038i


0.0006 - 0.0257i 0.0006 + 0.0257i
0.1455 + 0.0440i 0.1455 - 0.0440i
0.9880 + 0.0000i 0.9880 + 0.0000i

d=
-4.7176 +24.2309i 0.0000 + 0.0000i
0.0000 + 0.0000i -4.7176 -24.2309i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i

0.0000 + 0.0000i
0.0000 + 0.0000i
0.8434 +38.4369i
0.0000 + 0.0000i

0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.8434 -38.4369i

AT v=160,
lambda =
-8.2329 +19.6889i
-8.2329 -19.6889i
3.3902 +33.1394i
3.3902 -33.1394i
Q=
0
0 1.0000
0
0
0
0 1.0000
-631.6412 -211.1150 -6.1575 0.0000
0.0000 -800.1396 38.0318 -3.5280
>> [v,d]=eig(Q)
v=
-0.0066 + 0.0243i -0.0066 - 0.0243i
-0.0152 - 0.0364i -0.0152 + 0.0364i
-0.4244 - 0.3299i -0.4244 + 0.3299i
0.8420 + 0.0000i 0.8420 + 0.0000i

0.0076 - 0.0062i 0.0076 + 0.0062i


0.0029 - 0.0282i 0.0029 + 0.0282i
0.2325 + 0.2296i 0.2325 - 0.2296i
0.9446 + 0.0000i 0.9446 + 0.0000i

d=
-8.2329 +19.6889i 0.0000 + 0.0000i
0.0000 + 0.0000i -8.2329 -19.6889i
0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i

0.0000 + 0.0000i
0.0000 + 0.0000i
3.3902 +33.1394i
0.0000 + 0.0000i

0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
3.3902 -33.1394i

At v = 192,
lambda =
4.2758 +28.7405i
4.2758 -28.7405i
-18.7731 + 0.0000i
-1.4011 + 0.0000i
Q=
0
0 1.0000
0
0
0
0 1.0000
-631.6412 -304.0056 -7.3890 0.0000
0.0000 -35.1581 45.6381 -4.2336
>> [v,d]=eig(Q)
v=
-0.0182
-0.0043
-0.1608
-0.8445

+ 0.0029i -0.0182 - 0.0029i 0.0180 + 0.0000i -0.2547 + 0.0000i


+ 0.0287i -0.0043 - 0.0287i -0.0501 + 0.0000i 0.5221 + 0.0000i
- 0.5097i -0.1608 + 0.5097i -0.3379 + 0.0000i 0.3568 + 0.0000i
+ 0.0000i -0.8445 + 0.0000i 0.9397 + 0.0000i -0.7316 + 0.0000i

d=
4.2758 +28.7405i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i

0.0000 + 0.0000i 4.2758 -28.7405i 0.0000 + 0.0000i 0.0000 + 0.0000i


0.0000 + 0.0000i 0.0000 + 0.0000i -18.7731 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -1.4011 + 0.0000i
At v =224,
lambda =
-39.1714 + 0.0000i
2.7629 +26.2666i
2.7629 -26.2666i
20.0859 + 0.0000i
Q=
0
0 1.0000
0
0
0
0 1.0000
-631.6412 -413.7855 -8.6205 0.0000
0.0000 868.9109 53.2445 -4.9392
>> [v,d]=eig(Q)
v=
0.0056 + 0.0000i
-0.0249 + 0.0000i
-0.2207 + 0.0000i
0.9750 + 0.0000i

0.0029 - 0.0280i 0.0029 + 0.0280i -0.0161 + 0.0000i


-0.0249 - 0.0045i -0.0249 + 0.0045i 0.0470 + 0.0000i
0.7426 + 0.0000i 0.7426 + 0.0000i -0.3236 + 0.0000i
0.0482 - 0.6669i 0.0482 + 0.6669i 0.9449 + 0.0000i

d=
-39.1714 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i
2.7629 +26.2666i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i 0.0000 + 0.0000i


0.0000 + 0.0000i 0.0000 + 0.0000i
2.7629 -26.2666i 0.0000 + 0.0000i
0.0000 + 0.0000i 20.0859 + 0.0000i

At v =256,
lambda =
-52.9952 + 0.0000i
1.4834 +25.6470i
1.4834 -25.6470i
34.5315 + 0.0000i
Q=
1.0e+03 *
0
0 0.0010
0
0
0
0 0.0010
-0.6316 -0.5405 -0.0099 0.0000
0.0000 1.9121 0.0609 -0.0056
>> [v,d]=eig(Q)
v=
0.0034 + 0.0000i
-0.0186 + 0.0000i
-0.1821 + 0.0000i
0.9831 + 0.0000i
d=

0.0019 - 0.0332i 0.0019 + 0.0332i -0.0070 + 0.0000i


-0.0201 - 0.0017i -0.0201 + 0.0017i 0.0281 + 0.0000i
0.8536 + 0.0000i 0.8536 + 0.0000i -0.2422 + 0.0000i
0.0147 - 0.5193i 0.0147 + 0.5193i 0.9698 + 0.0000i

-52.9952 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i
1.4834 +25.6470i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i 0.0000 + 0.0000i


0.0000 + 0.0000i 0.0000 + 0.0000i
1.4834 -25.6470i 0.0000 + 0.0000i
0.0000 + 0.0000i 34.5315 + 0.0000i

At v =288,
lambda =
-65.1343 + 0.0000i
0.7147 +25.4558i
0.7147 -25.4558i
46.2710 + 0.0000i
>> Q
Q=
1.0e+03 *
0
0 0.0010
0
0
0
0 0.0010
-0.6316 -0.6840 -0.0111 0.0000
0.0000 3.0943 0.0685 -0.0064
>> [v,d]=eig(Q)
v=
0.0025 + 0.0000i
-0.0151 + 0.0000i
-0.1625 + 0.0000i
0.9866 + 0.0000i

0.0010 - 0.0356i 0.0010 + 0.0356i -0.0044 + 0.0000i


-0.0165 - 0.0009i -0.0165 + 0.0009i 0.0212 + 0.0000i
0.9058 + 0.0000i 0.9058 + 0.0000i -0.2038 + 0.0000i
0.0105 - 0.4218i 0.0105 + 0.4218i 0.9788 + 0.0000i

d=
-65.1343 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i
0.7147 +25.4558i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i 0.0000 + 0.0000i


0.0000 + 0.0000i 0.0000 + 0.0000i
0.7147 -25.4558i 0.0000 + 0.0000i
0.0000 + 0.0000i 46.2710 + 0.0000i

At v =320,
lambda =
-76.4353 + 0.0000i
0.1851 +25.3691i
0.1851 -25.3691i
56.6940 + 0.0000i
Q=
1.0e+03 *
0
0 0.0010
0
0
0
0 0.0010
-0.6316 -0.8445 -0.0123 0.0000
0.0000 4.4156 0.0761 -0.0071
>> [v,d]=eig(Q)
v=
0.0020 + 0.0000i 0.0003 - 0.0368i 0.0003 + 0.0368i -0.0032 + 0.0000i

-0.0129 + 0.0000i -0.0140 - 0.0005i -0.0140 + 0.0005i 0.0173 + 0.0000i


-0.1509 + 0.0000i 0.9337 + 0.0000i 0.9337 + 0.0000i -0.1827 + 0.0000i
0.9885 + 0.0000i 0.0107 - 0.3558i 0.0107 + 0.3558i 0.9830 + 0.0000i
d=
-76.4353 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i
0.1851 +25.3691i
0.0000 + 0.0000i
0.0000 + 0.0000i

0.0000 + 0.0000i 0.0000 + 0.0000i


0.0000 + 0.0000i 0.0000 + 0.0000i
0.1851 -25.3691i 0.0000 + 0.0000i
0.0000 + 0.0000i 56.6940 + 0.0000i

You might also like