UUM 523E Advanced Flight Dynamics Assignment 4

(Numerical Linearization)

Çağrı Ege Altunkaya 110160104

08 January 2021
 Problem:Write a generic m-file code for the numerical linearization around the trim conditions.
1.Define your linearized model as follows,

Use body-fixed axes system for translational dynamics. Find the linearized models around the trim condition at the
reference c.g. location (0.35c) for,

a. M = 0.6, h = 10000f t, γ = 0°, φ = 0° ((steady-state wings level flight)

 
−0.0103 −0.0743 0 0 0 0 0 0
 0.0554
 −1.0099 0.0176 0 0 0 0 0 
 −4.672 183.269 −1.0215 1 0 0 0.0025 0 
 
−9.8034 −0.2446 0 0 0 0 0 0 
A: 

 0 −0.002 0 0 −0.2861 −0.1773 0.0508 0 

 0 0 0 0 4.7984 −3.015 0.0197 1 
 0 0 −0.0029 0 −195.8697 0.3113 −0.9757 0.0249
0 0 0 0 9.7871 0 0 0

 
0.0631 0.0641 0 0
 0
 −0.455 0 0 

 0
 −0.2276 0 0 

 0 0 0 0 
B:
 0

 0 0.0536 0.1455 

 0
 0 −0.8561 0.1741 

 0 0 −0.0112 −0.0813
0 0 0 0

Eigenvalues:
p1,2 = −0.6470 ± 3.2331i
p3 = −2.8123
p4 = −2.9357
p5 = 0.7877
p6,7 = −0.0086 ± 0.0875i
p8 = −0.0471

A21 and A12 matrices are nearly zero matrices, and decoupled characteristics is observed which is expected in
straight flight trim; therefore, any coupling effect in modes is not expected. The poles of the system are located in
left-half plane except of one of them which is spiral mode pole. Any perturbation in bank angle causes spiral divergence
motion of aircraft.

1
 b. M = 0.6, h = 10000f t, γ = 0°, β = 0°, φ = 60° ((steady-state turn flight)

 
−0.0085 −0.0451 −0.0008 0 −0.0432 −0.0004 0 0
 −0.0066 −1.0307 0.0168 0 −0.0022 0.0019 0.0001 0 
 
−13.5791 182.2 −1.024 0.5 0 −0.0334 0.0054 0.062 
 
 −9.7814 −0.3503 0 0 −0.6067 0 0 0.0865
A: 
 0.0431
 −0.0083 −0.0005 0 −0.2841 −0.2071 0.0502 0  

 0 0 0.0413 0 14.1585 −3.3777 −0.0436 1  
 0 0 −0.0065 −0.866 −195.3721 0.6773 −0.8872 0.0358
0 −8.4555 0 −0.0859 4.8818 0 0 0

 
0.0631 0.0579 0 0
 0
 −0.4298 0 0 

 0
 −0.2241 0 0 

 0 0 0 0 
B:
 0

 0 0.0518 0.1481 

 0
 0 −0.8915 0.1698 

 0 0 −0.0087 −0.0812
0 0 0 0
Eigenvalues:
p1,2 = −0.7182 ± 3.4568i
p3 = −0.7182
p4 = 0.7326
p5,6 = −0.0101 ± 1079i
p7 = −2.7848
p8 = −3.0933

A21 and A12 matrices are different than zero and they are in non-negligible condition. Coupled characteristics is
observed which is expected in steady turn trim; therefore, coupling effect in longitudinal and lateral/directional modes
is expected. The poles of the system are located in left-half plane except of one of them which is spiral mode pole.
Any perturbation in bank angle causes spiral divergence motion of aircraft.

2
 2. Define your linearized model as follows,

Use wind axes system for translational dynamics. Find the linearized models around the trim condition at the refer-
ence c.g. location (0.35c) for,

a. M = 0.6, h = 10000f t, γ = 0°, φ = 0° ((steady-state wings level flight)

 
−0.0157 −0.0005 0 0 0 0 0 0
 6.02
 −1.0085 3.4712 0 0 0 0 0 
−0.0999 0.9304 −1.0215 1 0 −0.0003 0.0025 0 
 
−9.8065 0 0 0 0 0 0 0 
A:
−0.3376

 −0.0005 0.003 0 −0.2875 −39.0700 11.8773 0 

 0 0 0 0 0.0244 −3.0150 0.0197 1 
 0 0 −0.0029 0 −0.9941 0.3113 −0.9757 0.0249
0 0 0 0 0.0497 0 0 0

 
0.0631 0.0527 0 0
 0
 −0.0023 0 0 

 0
 −0.2276 0 0 

 0 0 0 0 
B:
 0

 0 0.0003 0.0007 

 0
 0 −0.8561 0.1741 

 0 0 −0.0112 −0.0813
0 0 0 0
Eigenvalues:
p1,2 = −0.6491 ± 3.5245i
p3 = 0.7879
p4,5 = −0.0107 ± 0.0875i
p6 = −0.0442
p7 = −2.8122
p8 = −2.9359

A21 and A12 matrices are nearly zero matrices, and decoupled characteristics is observed which is expected in
straight flight trim; therefore, any coupling effect in modes is not expected. The poles of the system are located in
left-half plane except of one of them which is spiral mode pole. Any perturbation in bank angle causes spiral divergence
motion of aircraft.

3
 b. M = 0.6, h = 10000f t, γ = 0°, β = 0°, φ = 60° ((steady-state turn flight)

 
−0.0210 −0.0006 −0.0004 0 −0.0002 −0.0002 0 0
−15.5947 −1.0218 3.3207 0 0.0009 0.3715 0.0130 0 
 
 −0.5283 0.9273 −1.0240 0.5 0 −0.0333 0.0054 0.0227
 
 −9.7873 0.0029 0 0 −0.0011 0 0 0.0862
A: 
 8.2578
 −0.0066 −0.0562 0 −0.2867 −43.7680 11.9487 0 

 0 0 0.0413 0 0.0719 −3.3776 −0.0436 1 
 0 0 −0.0065 −0.866 −0.9916 0.6773 −0.8872 0.0131
−0.6054 −0.0429 0 −0.0859 0.0248 0 0 0

 
0.0629 0.0271 0 0
 0
 −0.0022 0 0 

 0
 −0.2241 0 0 

 0 0 0 0 
B:
 0

 0 0.0003 0.0008 

 0
 0 −0.8915 0.1698 

 0 0 −0.0087 −0.0812
0 0 0 0

Eigenvalues:
p1,2 = −0.7149 ± 3.7708i
p3 = 0.7354
p4,5 = −0.0111 ± 0.1073i
p6 = −0.0151
p7 = −2.7837
p8 = −3.1032

A21 and A12 matrices are different than zero and they are in non-negligible condition. Coupled characteristics is
observed which is expected in steady turn trim; therefore, coupling effect in longitudinal and lateral/directional modes
is expected. The poles of the system are located in left-half plane except of one of them which is spiral mode pole.
Any perturbation in bank angle causes spiral divergence motion of aircraft.