Shimmy analysis of light weight aircraft nose wheel landing gear

Giridharan, V.; Sivakumar, S.

doi:10.21595/vp.2022.22988

Vibroengineering Procedia

Browse Procedia

Published: 12 December 2022

Check for updates

Shimmy analysis of light weight aircraft nose wheel landing gear

V. Giridharan¹

S. Sivakumar²

^{1, 2}Department of Aerospace Engineering, SRMIST, Kattankulathur, India

Corresponding Author:

S. Sivakumar

Cite the article Download PDF

Downloads 786

CrossRef Citations 1

Abstract

This paper presents mathematical modeling and analysis of shimmy oscillations for a light weight airplane single wheel nose landing gear. Shimmy is a self-excited oscillation which occurs usually on the nose wheel landing gear assembly during ground maneuvers which is governed by the dynamic characteristics of the landing gear and tires. Shimmy oscillation may lead to reduce the fatigue life of the landing gear and fuselage structure. So, the study of dynamic response and stability boundaries of landing gear plays crucial role while designing of airplanes. In earlier studies of vehicle shimmy only 3 degrees of freedom (DOF) considered such as torsional mode, lateral bending mode and tire lateral deformation. In this work along with above mentioned DOF, two more additional DOF introduced such as axial vibration of strut and tire in order to include the effect of vertical dynamics on shimmy model. Gyroscopic coupling effect also included in the model to study its influence on shimmy. Analysis carried out to determine critical velocity region for occurrence of shimmy and to investigate the effectiveness of ground unevenness on the landing gear system for two different runway conditions such as flat runway and random roughness runway. The results are more helpful to study significant interaction between the different parameters of landing gear and to represent stability boundaries.

1. Introduction

Lightweight airplanes are extensively used for surveillance, reconnaissance, defense and military. Development of airplanes integrating with a critical component like landing gear is a challenging task. Airplanes may be operated in both prepared and unprepared runways for takeoff, landing and ground maneuvers. So, efficient, good, robust and easy to maintain landing gear is required to ensure safety and comfort. A literature review is made on the previously published articles for wheel shimmy. In [1] the maximum total force with the random runway roughness of landing gears are obtained. Stability analysis of main landing gears are illustrated in [2]. In [3] stability of the system determined based on eigenvalues. Analysis of nonlinear differential equations using linearization techniques is presented in [4]. Linear and non-linear stability analysis of the nose landing gear model discussed in [5]. Percentages of the stable regions are computed by stability analysis conducted in different planes with different parameters given in [6]. Modeling of shimmy analysis is presented in [7]. Stability analysis for airplane with landing gear and shimmy damper illustrated in [8]. From the stability analysis limit of the shimmy amplitude and the magnitude of initial disturbance to initiate the self-excited oscillation are determined in [9].Various analytical and graphical methods to determine the stability margins for nonlinear systems are presented in [10].Vehicle shimmy with the effect of nonlinear tire forces, and the steering mechanism are investigated in [11] with the tire lateral forces are modeled by using cubic equation rather than popular magic formula. In [12] the potential effect of the two structural modifications to the parameters design ranges for the shimmy damper are investigated. Landing gear dynamics with torsional DOF and freeplay is investigated in [13]. Tire is modeled using the elastic string theory. Forces acting on the wheel are transmitted from the ground through the tire, and these forces also deflect the tire. In [14] active landing gears random vibration analysis are presented. Active landing gear responses are compared to the passive system. In [15] lateral response of the landing gear investigated with the effect of ground unevenness using mathematical model. Summary of the literature survey on solving landing gear vibration problems are presented in [16]. In this paper, mathematical model is developed to analyze the shimmy oscillation of nose wheel landing gear with consideration of total 5 DOF such as torsional angle $ψ$ , lateral bending angle $δ$ , axial displacement of strut $x_{s}$ , axial displacement of tire $x_{t}$ and slip angle $α$ . Dynamics equations derived using second Lagrange equation. The stability boundaries of the model are determined based on the critical velocity regions for occurrence of shimmy. Dynamic responses of the model on two different runways in both stable and unstable regions of velocity are compared.

2. Dynamic model of nose landing gear

Simple trailing wheel system shown in Fig. 1 represents side view, top view, and front view along with direction of vertical load, forward velocity, torsion angle and lateral bending angle. Mathematical model of nose wheel landing gear shimmy due to axial, torsional and lateral vibration is shown in Fig. 2. The model of consists of wheel with tire of radius $R$ fitted with axle of the shock strut. The single point contact tire model is considered for analysis. The vertical load $F_{z}$ acts on landing gear comes from airframe and aircraft forward velocity $V$ is given as inputs to model. Dynamic model considers five DOF such as the shimmy angle of the nose wheel landing gear $ψ$ , lateral bending angle $δ$ , the axial displacement of the tire $x_{t}$ , the axial displacement of the strut $x_{s}$ and slip angle $α$ of the leading contact point. $K_{ψ}$ , $K_{δ}$ , $K_{s}$ , $C_{ψ}$ , $C_{δ}$ , $C_{s}$ represents the torsional, bending, vertical stiffness and damping coefficients for nose landing gear strut respectively. The damping and stiffness coefficient for tire are mentioned as $C_{t}$ and $K_{t}$ . $M_{s}$ and $M_{u}$ represents sprung mass corresponding to nose landing gear and unsprung mass of wheel assembly.

Fig. 1Simple trailing wheel system

a) Top view

b) Side view

c) Front view

Fig. 2Mathematical model of landing gear with axial, torsional and lateral vibration

Based on Lagrange’s principle, equations of motion of damped multi degree-of-freedom system are written as shown in Eq. (1-9):

1

\frac{d}{d t} (\frac{\partial T}{{\dot{\partial x}}_{i}}) - (\frac{\partial T}{{\partial x}_{i}}) + \frac{\partial U}{{\partial x}_{i}} + \frac{\partial D}{{\dot{\partial x}}_{i}} = Q_{i}, i = 1,2, 3, \dots, n,

where $T$ is kinetic energy of the model, $U$ is potential energy of the model, $D$ is dissipative potential function of the model, $Q_{i}$ is the generalized force to which each DOF of the system is subjected, and $x_{i}$ represents the generalized coordinate of the system, which is expressed $[x_{1}, x_{2}, x_{3}, x_{4}]$ as ${[ψ, y, x}_{t}, x_{s}$ ]:

2

T = \frac{1}{2} I_{z} {\dot{ψ}}^{2} + \frac{1}{2} I_{x} {\dot{δ}}^{2} + \frac{1}{2} M_{u} {\dot{x}}_{t}^{2} + \frac{1}{2} M_{s} {\dot{x}}_{s}^{2},

3

U = \frac{1}{2} K_{ψ} ψ^{2} + \frac{1}{2} K_{δ} δ^{2} + \frac{1}{2} K_{t} ({x_{t} - q)}^{2} + \frac{1}{2} K_{s} {(x_{s} - x_{t})}^{2},

4

D = \frac{1}{2} C_{ψ} {\dot{ψ}}^{2} + \frac{1}{2} C_{δ} {\dot{δ}}^{2} + \frac{1}{2} C_{s} ({{\dot{x}}_{s} - {\dot{x}}_{t})}^{2} + \frac{1}{2} C_{t} ({{\dot{x}}_{t} - \dot{q})}^{2} .

The differential equations of motion of the landing gear model is derived based on the second Lagrange equation:

5

I_{z} \ddot{ψ} + C_{ψ} \dot{ψ} + K_{ψ} ψ + \frac{I_{p}}{B} * \frac{v}{R} * ({\dot{x}}_{t}) + F_{R} (R τ + n) - F_{v} A τ = 0,

where $F_{R} = [c_{1} α + c_{2} (F_{z_{o}} - F_{z})] α + c_{3} α^{3} F_{R} = K_{1} α + K_{2} (A τ ψ + x_{t} - q) α + c_{3} α^{3}$ , $c_{1}$ - $c_{3}$ are coefficients for tire formula, and $F_{z o}$ is the nominal load of the tire, then $K_{1} = c_{1} + c_{2} F_{z 0}$ and $K_{2} = c_{2} k_{t}$ , $K_{3} = c_{3}$ , and:

6

I_{x} \ddot{δ} + C_{δ} \dot{δ} + K_{δ} δ + \frac{I_{p}}{B} * \frac{v}{R} * ({\dot{x}}_{t}) - F_{z} e_{e f f} s i n θ + l_{g} F_{R} c o s θ c o s τ = 0,

where $e_{e f f} = e c o s τ + (R + e s i n τ) t a n τ$ , $θ = ψ c o s τ$ , and:

7

M_{u} {\ddot{x}}_{t} - C_{s} ({\dot{x}}_{s} - {\dot{x}}_{t}) {+ C}_{t} ({\dot{x}}_{t} - \dot{q}) - K_{s} (x_{s} - x_{t}) + K_{t} (x_{t} - q) - \frac{I_{p}}{B} \frac{v}{R} \dot{ψ} + {(K}_{t} A τ ψ) = 0,

8

M_{s} {\ddot{x}}_{s} + C_{s} ({\dot{x}}_{s} - {\dot{x}}_{t}) - K_{s} (x_{s} - x_{t}) = 0 .

The relation governing stretched string tire model can be expressed as:

9

\dot{α} + \frac{v}{σ} α - v ψ + a \dot{ψ} = 0 .

3. Numerical analysis

Based on the mathematical model the dynamic response obtained using MATLAB/Simulink. The developed Simulink model is shown in Fig. 3. Using solver ode 45 equations of motion are resolved. Analysis of nose wheel landing gear carried out by considering initial value of the torsion angle as 0.001 radians on two different runway conditions such as flat runway and random roughness runway. The random runway roughness can be modeled by PSD function. By measuring the surface profile PSD function can be determined, Quality of runway is graded based on roughness variance and roughness index in the domain of power spectral density as illustrated in [17]. By increasing the roughness variance the runway characteristics changes from very good to very poor. Runway grade considered in this work is very poor grade with roughness variance = 0.032 m; roughness inde $x =$ 0.127 rad/m The parameters used for analysis are taken from [7] and the tire parameters are used from [9] which are listed in Table 1.

Fig. 3Simulink model for Nose wheel landing gear

Table 1Landing gear and tire parameters

Parameter	Value	Unit
Forward velocity, $V$	0-80	m/s
Vertical force, $F_{z}$	8000-12000	N
Sprung mass $M_{s}$	750	Kg
Unsprung mass $M_{u}$	38	kg
Caster length, $e$	0.1	m
Caster angle, $τ$	0.08	rad
Distance from the wheel center to strut, $A$	0.068	m
Distance from the wheel center to the hinge center, $l g$ / $B$	0.48	m
Strut torsional stiffness, $K_{ψ}$	100000	Nm./rad
Strut torsional damping, $C_{ψ}$	45	Nm.s/rad
Strut vertical stiffness, $K_{s}$	31600	N/m
Strut vertical damping, $C_{s}$	2300	Ns/m
Strut lateral bending stiffness, $K_{δ}$	3240000	Nm./r
Strut lateral bending damping, $C_{δ}$	1	Nm.s/rad
Strut moment of inertia about $z$ -axis, $I_{z}$	1	kg.m2
Strut moment of inertia about $y$ -axis, $I_{y}$	1	kg.m2
Strut moment of inertia about $x$ -axis, $I_{x}$	1	kg.m2
Tire parameters
Radius of nose wheel, $R$	0.36	m
Contact patch length, $a$	0.1	m
Relaxation length, $σ = 3 * a$	0.3	m
Vertical stiffness of tyre $K_{t}$	260	KN/m
Vertical damping of tyre $C_{t}$	4.066	KN.s/m
Tire coefficient $C_{1}$	–30.30	KN/rad
Tire coefficient $C_{2}$	–7.96	1/rad
Tire coefficient $C_{3}$	3250	KN/rad^3
Nominal load of tire $F_{z o}$	4700	N

4. Results and discussion

Parameters given in Table 1 are used for analysis of nose landing gear. The shimmy model solved in Matlab/Simulink. Dynamic-response to study the effect of runway roughness and its interaction with shimmy amplitude are presented. Fig. 4 shows that velocity range of shimmy will be enlarged when sprung mass $M_{s}$ increases for the system with sprung mass of 750 kg, it can be observed from figure 4a that the shimmy angle amplitude is 0.045 radian for flat runway whereas for random roughness runway shimmy angle amplitude is 0.15 radian. In both the runways velocity less than 54 km/h the system is stable and from 54 km/h to 152 km/h, the system is unstable, and system becomes stable above 152 km/h. It is found that the critical velocity of instability is from 54 km/h to 152 km/h for the case of NLG configuration of sprung mass 750 kg. Similarly for the system with sprung mass of 1000 kg, it can be observed from figure 4a that the shimmy angle amplitude is 0.18 radian for flat runway whereas for random roughness runway shimmy angle amplitude is 0.24 radian. In both the runways velocity less than 33 km/h the system is stable and from 33 km/h to 210 km/h, the system is unstable, and system becomes stable above 210 km/h. The critical velocity of instability is from 33 km/h to 210 km/h for the case of NLG configuration of sprung mass 1000 kg. These results show that the amplitude of shimmy angle increases dramatically for random roughness runway when compared to the flat runway even though critical velocity region have less effect on runway roughness.

Fig. 4Aircraft velocity vs shimmy angle for two different mass

a) Flat runway

b) Random roughness runway

Fig. 5Time response for torsional angle ψ

Fig. 6Phase portrait plot for torsional angle ψ

The study is made to compare the effect of two different runways at stable and unstable velocity regions. For system with sprung mass of 1000 kg, time response of shimmy angle on flat runway and random roughness runway at velocity 30 km/h and 100 km/h are shown in Fig. 5(a-d). From the Fig. 5 it is noted that in the critical velocity region time response diverges and the system exhibits unstable nature. Phase portrait plot of shimmy angle $ψ$ shown in Fig. 6(a-d). At the forward velocity of 30 km/h, the limit cycle is approached towards inside for both runways and at 100 km/h the limit cycle moves from inside to outwards on both runways. From the results shown in Figs. 7-14 it is observed that landing gear systems other shimmy influencing parameters such as lateral bending angle, tire displacement, strut displacement and slip angle makes nose landing gear system stable in the velocity lower than critical velocity region and the system becomes unstable upto higher critical velocity irrespective of the grade of runways.

Fig. 7Time response for bending angle δ

Fig. 8Phase portrait plot for bending angle δ

Fig. 9Time response for tire displacement xt

Fig. 10Phase portrait plot for tire displacement xt

Fig. 11Time response for strut displacement, xs

Fig. 12Phase portrait plot for strut displacement, xs

5. Conclusions

In this paper mathematical model of an aircraft nose landing gear developed with consideration of the vertical load acting on the wheel and the gyroscopic effects due to a rotation of wheel. Interaction of torsion, lateral, axial vibration of sprung mass, axial vibration of un sprung mass are investigated.

Fig. 13Time response for slip angle α

Fig. 14Phase portrait plot for slip angle α

The dynamic behaviour of shimmy model is studied for $ψ$ , $δ$ , $x_{t}$ , $x_{s}$ , $α$ and are compared on two different runways in both stable and unstable regions of velocity. Based on the simulation results of mathematical model the following conclusions are made:

1. Critical velocity region for occurrence of shimmy enlarged when sprung mass $M_{s}$ , increases.

2. Amplitude of self-excited shimmy angle also increases by increasing sprung mass $M_{s}$ .

3. The amplitude of shimmy angle increases dramatically for random roughness runway when compared to the flat runway.

4. Dynamic responses decays when the velocity is in stable region and diverges when the velocity is in unstable region.

5. Phase portrait shows that the trajectory spiral up inwards for stable velocity regions and outward for unstable velocity regions.

6. Fatigue life of landing gear and fuselage structure may be improved if airplane is operated at velocity lower than critical velocity during ground roll.

7. Future work to study the influence of various landing gear parameters on shimmy vibration and to analyze stability nature.

References

S. Sivakumar, T. Selvakumaran, and B. Sanjay, “Investigation of random runway effect on landing of an aircraft with active landing gears using nonlinear mathematical model,” Journal of Vibroengineering, Vol. 23, No. 8, pp. 1785–1799, Dec. 2021, https://doi.org/10.21595/jve.2021.21915

Publisher
Ijm Igo Besselink, “Shimmy of aircraft main landing gears,” Ph.D. Thesis, Technische Universiteit Delft, 2000.

Search CrossRef
C. Arreaza, K. Behdinan, and J. W. Zu, “Linear stability analysis and dynamic response of shimmy dampers for main landing gears,” Journal of Applied Mechanics, Vol. 83, No. 8, pp. 1–10, Aug. 2016, https://doi.org/10.1115/1.4033482

Publisher
R. Morgan, “Linearization and stability analysis of nonlinear problems,” Rose-Hulman Undergraduate Mathematics Journal, Vol. 16, No. 2, pp. 67–91, 2015.

Search CrossRef
P. Thota, B. Krauskopf, and M. H. Lowenberg, “Shimmy in a nonlinear model of an aircraft nose landing gear with non-zero rake angle,” in 6th EUROMECH (European Mechanics Society) Nonlinear Dynamics Conference, 2008.

Search CrossRef
G. Somieski, “Shimmy analysis of a simple aircraft nose landing gear model using different mathematical methods,” Aerospace Science and Technology, Vol. 1, No. 8, pp. 545–555, Dec. 1997, https://doi.org/10.1016/s1270-9638(97)90003-1

Publisher
H. Georgieva, “Modeling of shimmy oscillations in aircraft landing gear,” MATEC Web of Conferences, Vol. 133, p. 01007, 2017, https://doi.org/10.1051/matecconf/201713301007

Publisher
N. K. Sura and S. Suryanarayan, “Lateral stability of aircraft nose-wheel landing gear with closed-loop shimmy damper,” Journal of Aircraft, Vol. 46, No. 2, pp. 505–509, Mar. 2009, https://doi.org/10.2514/1.37626

Publisher
H. Wei, J.-W. Lu, S.-Y. Ye, and H.-Y. Lu, “Bifurcation analysis of vehicle shimmy system exposed to road roughness excitation,” Journal of Vibration and Control, Vol. 28, No. 9-10, pp. 1045–1056, May 2022, https://doi.org/10.1177/1077546320987791

Publisher
G. Somieski, “An Eigenvalue method for calculation of stability and limit cycles in nonlinear systems,” Nonlinear Dynamics, Vol. 26, No. 1, pp. 3–22, 2001, https://doi.org/10.1023/a:1017384211491

Publisher
L. Dai and Q. Han, “Stability and Hopf bifurcation of a nonlinear model for a four-wheel-steering vehicle system,” Communications in Nonlinear Science and Numerical Simulation, Vol. 9, No. 3, pp. 331–341, Jun. 2004, https://doi.org/10.1016/s1007-5704(02)00084-9

Publisher
Yang Chen, “Time response simulation for nose landing gear shimmy damper design,” Thesis of Master of Applied Science, University of Toronto, 2018.

Search CrossRef
Freeplay, Elmas Atabay, and I. Ozkol, “On dynamics of a landing gear mechanism with torsional freeplay,” Global Journal of Researches in Engineering Aerospace Engineering, Vol. 12, No. 1, p. 16, 2012.

Search CrossRef
S. Sivakumar and A. P. Haran, “Aircraft random vibration analysis using active landing gears,” Journal of Low Frequency Noise, Vibration and Active Control, Vol. 34, No. 3, pp. 307–322, Jun. 2015, https://doi.org/10.1260/0263-0923.34.3.307

Publisher
N. K. Sura and S. Suryanarayan, “Lateral response of nose-wheel landing gear system to ground-induced excitation,” Journal of Aircraft, Vol. 44, No. 6, pp. 1998–2005, Nov. 2007, https://doi.org/10.2514/1.28854

Publisher
J. Pritchard, “Overview of landing gear dynamics,” Journal of Aircraft, Vol. 38, No. 1, pp. 130–137, Jan. 2001, https://doi.org/10.2514/2.2744

Publisher
F. Tyan, Yu Hong, Shun-Hsu Tu, and W. Jeng, “Generation of random road profiles,” Journal of Advanced Engineering, Vol. 4, No. 2, pp. 151–156, 2009.

Search CrossRef

Cited by

Shimmy vibration analysis of unmanned aircraft coupled with landing gears

S. Sivakumar | L. R. Ganapathy Subramanian | V. Giridharan

(2024)

About this article

Received

15 October 2022

Accepted

03 November 2022

Published

12 December 2022

SUBJECTS

Vibration in transportation engineering

DOI

https://doi.org/10.21595/vp.2022.22988

Keywords

torsion

bending

axial displacement

gyroscopic effect and critical velocity

Acknowledgements

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous article in issue Previous Next article in issue Next

Research article

2022 06 30

Vibration characteristics analysis of face gear transmission system considering gyroscopic effect

Xiaozhen Li, Kuo Su, Tongjie Li, Liyong Zhang

Research article

2021 09 09

Investigation of random runway effect on landing of an aircraft with active landing gears using nonlinear mathematical model

S. Sivakumar, T. Selvakumaran, B. Sanjay

Research article

2018 06 30

Non-linear vibration analysis of oleo pneumatic landing gear at touchdown impact

S. Sivakumar, M. Syedhaleem

Research article

2017 09 30

Analytical coupling characterization of multi-stage planetary gear free vibration considering flexible structure

Aiqiang Zhang, Jing Wei, Datong Qin, Shaoshuai Hou

V. Giridharan and S. Sivakumar, “Shimmy analysis of light weight aircraft nose wheel landing gear,” Vibroengineering PROCEDIA, Vol. 47, pp. 8–15, Dec. 2022, https://doi.org/10.21595/vp.2022.22988

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/vp.2022.22988
UR  - https://doi.org/10.21595/vp.2022.22988
TI  - Shimmy analysis of light weight aircraft nose wheel landing gear
T2  - Vibroengineering PROCEDIA
AU  - Giridharan, V.
AU  - Sivakumar, S.
PY  - 2022
DA  - 2022/12/12
PB  - JVE International Ltd.
SP  - 8-15
VL  - 47
SN  - 2345-0533
SN  - 2538-8479
ER  - 

Copy Ris

Copied to clipboard!

@article{Giridharan_2022,
	doi = {10.21595/vp.2022.22988},
	url = {https://doi.org/10.21595/vp.2022.22988},
	year = 2022,
	month = {dec},
	publisher = {{JVE} International Ltd.},
	volume = {47},
	pages = {8--15},
	author = {V. Giridharan and S. Sivakumar},
	title = {Shimmy analysis of light weight aircraft nose wheel landing gear},
	journal = {Vibroengineering {PROCEDIA}}
}

Copy Bibtex

Copied to clipboard!

[1]V. Giridharan and S. Sivakumar, “Shimmy analysis of light weight aircraft nose wheel landing gear,” Vibroengineering PROCEDIA, vol. 47, pp. 8–15, Dec. 2022, doi: 10.21595/vp.2022.22988.

Copy IEEE

Copied to clipboard!

Giridharan, V., and S. Sivakumar. “Shimmy Analysis of Light Weight Aircraft Nose Wheel Landing Gear.” Vibroengineering PROCEDIA 47 (December 12, 2022): 8–15. https://doi.org/10.21595/vp.2022.22988.

Copy Chicago

Copied to clipboard!

Shimmy analysis of light weight aircraft nose wheel landing gear

Abstract

1. Introduction

2. Dynamic model of nose landing gear

3. Numerical analysis

4. Results and discussion

5. Conclusions

References

Cited by

About this article

Related Articles