applied sciences
Article Parametric Analysis on Landing Gear Strut Friction of Light Aircraft for Touchdown Performance
Shengyong Gan 1, Xingbo Fang 1 and Xiaohui Wei 2,*
1 Key Laboratory of Fundamental Science for National Defense-Advanced Design Technology of Flight Vehicle, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China; [email protected] (S.G.); [email protected] (X.F.) 2 State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China * Correspondence: [email protected]
Abstract: The aim of this paper is to obtain the strut friction–touchdown performance relation for designing the parameters involving the strut friction of the landing gear in a light aircraft. The numerical model of the landing gear is validated by drop test of single half-axle landing gear, which is used to obtain the energy absorption properties of strut friction in the landing process. Parametric studies are conducted using the response surface method. Based on the design of the experiment results and response surface functions, the sensitivity analysis of the design variables is implemented. Furthermore, a multi-objective optimization is carried out for good touchdown performance. The results show that the proportion of energy absorption of friction load accounts for more than 35% of the total landing impact energy. The response surface model characterizes well for the landing response, with a minimum fitting accuracy of 99.52%. The most sensitive variables for the four
landing responses are the lower bearing width and the wheel moment of inertia. Moreover, the max overloading of sprung mass in LC-1 decreases by 4.84% after design optimization, which illustrates
Citation: Gan, S.; Fang, X.; Wei, X. that the method of analysis and optimization on the strut friction of landing gear is efficient for Parametric Analysis on Landing Gear improving the aircraft touchdown performance. Strut Friction of Light Aircraft for Touchdown Performance. Appl. Sci. Keywords: light aircraft; landing gear; response surface method (RSM); sensitivity analysis (SA) 2021, 11, 5445. https://doi.org/ 10.3390/app11125445
Academic Editor: Rosario Pecora 1. Introduction The aircraft landing gear plays a crucial role in the takeoff, landing impact energy- Received: 6 May 2021 absorbing, and ground operations [1]. In the shock absorber strut, the friction force is an Accepted: 9 June 2021 elemental force of the total strut force and converts kinetic energy into internal energy in Published: 11 June 2021 the landing process [2]. Therefore, friction can affect landing performance, especially in the light aircraft landing gear [3], which has a smaller structure stiffness and total strut force. Publisher’s Note: MDPI stays neutral To improve the landing performance by designing friction in the strut, a dynamic model with regard to jurisdictional claims in of the landing gear needs to be established to obtain the effect of friction load on the soft published maps and institutional affil- landing buffering; it is also required for the designer to determine the influence of landing iations. gear structure parameters on strut friction. The dynamic modeling and analysis of a landing gear with an oleo-pneumatic shock absorber has received relatively large attention in the literature [4–6]. The analysis model with sprung and unsprung mass vertical degrees of freedom is established in Milwitzky [7]. Copyright: © 2021 by the authors. Various studies have focused on predicting landing gear dynamics behavior during land- Licensee MDPI, Basel, Switzerland. ing [8–10] based on the numerical model. Elemental forces, including gas spring stiffness, This article is an open access article oleo damping, and friction, are expressed in the literature. In particular, the friction force is distributed under the terms and calculated as the absorber strut is regarded as a rigid body [9,11,12] or flexible body [13]. conditions of the Creative Commons Attribution (CC BY) license (https:// Response surface methodology (RSM) is an assemblage of statistical and mathematical creativecommons.org/licenses/by/ techniques utilized to develop, improve, and optimize processes of industrial produc- 4.0/). tion [14,15]. It is commonly used in the design process of some complex objects such as
Appl. Sci. 2021, 11, 5445. https://doi.org/10.3390/app11125445 https://www.mdpi.com/journal/applsci Appl. Sci. 2021, 11, 5445 2 of 18
planet lander [16,17] and aircraft [18,19]. RSM can reduce the computational expense for analysis by constructing a substitute model to approximately conform to the numerical model response based on the design of experiment results [20]. To acquire the relation be- tween the design parameters of landing gear with the strut friction force and the preferable touchdown performance, parametric studies on the half-axle aircraft main landing gear with oleo-pneumatic shock absorber are carried out using RSM. The studies are organized into two parts. Section2 presents the dynamics model of the main landing gear, including the dynamics equation and the elemental force expression. The method for analyzing normal forces in bearing contact area considering strut flexibility is expressed in Section 2.2. Moreover, a drop test is conducted to validate the numerical model in two landing con- ditions. Section3 first illustrates the effect of strut friction on landing gear touchdown performance. Second, Sections 3.2 and 3.3 present the response surface functions for the landing responses and the sensitivity analysis of the parameters, respectively. Finally, Section 3.4 presents a multi-objective optimization design of the parameters of the landing gear to improve the landing performance based on response surface functions.
2. Landing Gear Modeling and Validating 2.1. Landing Gear Dynamics Model In this subsection, the aircraft half-axle main landing gear (MLG) is selected as the research object, and a two-degrees-of-freedom spring damping model is established, as shown in Figure1. The aircraft landing buffer response is simplified as the dynamic response of two masses under spring damping constraint in this model [7]. The two masses are the sprung mass and unsprung mass. The sprung mass is the airframe structure and the structural mass above the outer cylinder of the shock absorber, and the unsprung mass is the other structure of the landing gear below the shock absorber piston rod. The oleo-pneumatic shock absorber consists of upper and lower chambers separated by orifices and the metering pin. The upper part of the upper chamber is filled with pressurized nitrogen to provide a gas spring, and the other spaces of the two chambers are filled with oil to provide oil damping. The tire force acting point of the half-axle landing gear is not Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 18 on the strut axis, which will augment the normal force in the bearing area. In this way, the half-axle landing gear can fully reflect the impact of the friction load on the landing buffer.
Figure 1. Schematic of the landing gear with two degrees of freedom.
The gas spring stiffness force is related to initial gas volume closely, the equation can be expressed as 𝑉 𝐹 = 𝐴 𝑃 −𝑃 , (2) 𝑉 − 𝐴 𝑆
where 𝐴 is the gas compressed area, 𝑃 is the gas initial pressure, 𝑉 is the gas initial volume, 𝑆 is the shock absorber stroke, 𝛾 is the gas polytropic exponent, and 𝑃 is the atmospheric pressure. The oil damping force is proportional to stroke slide velocity squared and parameters of the orifice, which is presented by the equation 𝐴 𝐹 =𝜌 𝑠 |𝑠 |, (3) 2𝐶 𝐴
where 𝜌 is the mass density of oil, 𝐴 is the oil compressed area of the absorber, 𝐴 is the orifice bypass area, 𝐶 is the oil damping discharge coefficient, and 𝑆 is the stroke slide velocity. The friction forces in the strut contain the journal friction force and seal friction force. The journal friction force is induced by the normal force acting on the upper and lower bearing area. The seal friction force result from the friction of internal seals in the shock absorber depends on the internal gas pressure. The friction forces in the shock strut are described by equation
𝐹 =𝐹 +𝐹 , (4)
where 𝐹 is the journal friction force and 𝐹 is the seal friction force. The seal friction force depends on the internal gas pressure [11] and is expressed as
𝐹 =−𝜇 𝐹 sgn( 𝑠 ), (5)
where 𝜇 is the seal friction coefficient and sgn is the signum function. The journal friction force is the product of the friction coefficient and the normal force. Due to the poor lubricating properties of hydraulic oil and the shape of the bearing surfaces in the shock absorber, it can be assumed that the shock–strut is under a dry fric- tion condition [7]. Coulomb’s law is the simplest and most utilized friction force model since it only requires one input parameter, i.e., the kinetic coefficient of friction [21]. The Appl. Sci. 2021, 11, 5445 3 of 18
According to the mathematical model shown in Figure1, the dynamic equilibrium governing equation of motion for the main landing gear is written as .. m1z1 = m1g − Fa − Fh − Ff , .. (1) m2z2 = −FV + m2g + Fa + Fh + Ff ,
where m1 and m2 indicate the sprung mass and unsprung mass, respectively; z1 and z2 denote the vertical displacement of the sprung mass and unsprung mass, respectively; Fa, Fh, and Ff represent the gas spring stiffness force, oil damping force, and friction force in the absorber strut, respectively. FV is the ground vertical force acting on the tire. The gas spring stiffness force is related to initial gas volume closely, the equation can be expressed as γ V0 Fa = Aa P0 − Patm , (2) V0 − AaS
where Aa is the gas compressed area, P0 is the gas initial pressure, V0 is the gas initial volume, S is the shock absorber stroke, γ is the gas polytropic exponent, and Patm is the atmospheric pressure. The oil damping force is proportional to stroke slide velocity squared and parameters of the orifice, which is presented by the equation
3 A . . = h Fh ρ 2 2 s s , (3) 2Cd Ad where ρ is the mass density of oil, A is the oil compressed area of the absorber, A is h . d the orifice bypass area, Cd is the oil damping discharge coefficient, and S is the stroke slide velocity. The friction forces in the strut contain the journal friction force and seal friction force. The journal friction force is induced by the normal force acting on the upper and lower bearing area. The seal friction force result from the friction of internal seals in the shock absorber depends on the internal gas pressure. The friction forces in the shock strut are described by equation Ff = Fn f + Fs f , (4)
where Fn f is the journal friction force and Fs f is the seal friction force. The seal friction force depends on the internal gas pressure [11] and is expressed as
. Fs f = −µs f Fasgn(s), (5)
where µs f is the seal friction coefficient and sgn is the signum function. The journal friction force is the product of the friction coefficient and the normal force. Due to the poor lubricating properties of hydraulic oil and the shape of the bearing surfaces in the shock absorber, it can be assumed that the shock–strut is under a dry friction condition [7]. Coulomb’s law is the simplest and most utilized friction force model since it only requires one input parameter, i.e., the kinetic coefficient of friction [21]. The model with the Stribeck effect reveals that the friction force decreases continuously with the increase of relative velocity from zero velocity, which can be described as
| | − v δ F = F + (F − F )e vs sgn(v) + σv, n f C S C (6) FS = µs N, FC = µk N,
where FS and FC represent the magnitude of static friction and Coulomb friction, respec- tively. µs denotes the static coefficient of friction which is higher than the kinetic coefficient µk. v is the relative velocity, vs is the Stribeck velocity, and δ is an exponent which depends on the geometry of the contacting surfaces, often considered to be equal to 2. σ indicates the viscous friction coefficient, which is neglected in this work as the dry friction is considered. Appl. Sci. 2021, 11, x FOR PEER REVIEW 4 of 18
model with the Stribeck effect reveals that the friction force decreases continuously with the increase of relative velocity from zero velocity, which can be described as
| | 𝐹 =(𝐹 +(𝐹 −𝐹 )𝑒 )sgn(𝑣)+𝜎𝑣, (6) 𝐹 =𝜇 𝑁, 𝐹 =𝜇 𝑁,
where 𝐹 and 𝐹 represent the magnitude of static friction and Coulomb friction, respec- tively. 𝜇 denotes the static coefficient of friction which is higher than the kinetic coeffi- cient 𝜇 . 𝑣 is the relative velocity, 𝑣 is the Stribeck velocity, and 𝛿 is an exponent Appl. Sci. 2021, 11, 5445 4 of 18 which depends on the geometry of the contacting surfaces, often considered to be equal to 2. 𝜎 indicates the viscous friction coefficient, which is neglected in this work as the dry friction is considered. FigureFigure 2a2 ashows shows the the classic classic shape shape of of the the Stri Stribeckbeck curve, curve, which which contains contains discontinuity discontinuity atat zero zero velocity. velocity. This This discontinuity bringsbrings aboutabout several several numerical numerical issues issues during during a dynamic a dy- namicsimulation. simulation. To eliminateTo eliminate the the numerical numerical issues issues at at zero zero velocity, velocity, a a finite finite slope slope model model is is establishedestablishedto to replace replace the the discontinuity discontinuity at zeroat zero velocity, velocity, as shown as shown in Figure in Figure2b. The 2b. model The modelutilized utilized in this in work this work can be can expressed be expressed by the by following the following equations: equations: 𝑣 𝐹 v sgn(𝑣) if |𝑣| <𝑣 , F S sgn(v) if |v| < v 0, 𝑣v 0 Fn𝐹 f == |v|−v (7)(7) −| | 0 v δ (FC + (FS − FC)e s )sgn(v) if |v| ≥ v0, (𝐹 +(𝐹 −𝐹 )𝑒 ) sgn( 𝑣) fi |𝑣| ≥𝑣 ,
wherewhere 𝑣 v 0 isis the the tolerance tolerance velocity. velocity.
(a) (b)
FigureFigure 2. Representation 2. Representation of Stribeck of Stribeck curves: curves: (a) The (a classic) The classicStribeck Stribeck curve; ( curve;b) Stribeck (b) Stribeckcurve with- curve outwithout discontinuity. discontinuity.
TheThe ground ground vertical vertical force force results results from from the the compression compression of of the the tires tires of of the the landing landing gear gear afterafter touching touching the the ground. ground. A A semi-empirical semi-empirical computational modelmodel [[22]22] cancan be be described described by bythe the equation equation . FV = 1 + CTz2 f (z2), (8) 𝐹 = (1+𝐶 z )𝑓(𝑧 ), (8) where C is the tire vertical damping deformation coefficient, and f (z ) is the tire vertical where 𝐶 T is the tire vertical damping deformation coefficient, and 𝑓(𝑧2) is the tire verti- static force corresponding to the tire compression amount. cal static force corresponding to the tire compression amount. The ground longitudinal force is the friction load caused by the relative rotation of the The ground longitudinal force is the friction load caused by the relative rotation of landing gear tire and the ground. Its magnitude is related to the ground vertical force and the landing gear tire and the ground. Its magnitude is related to the ground vertical force the ground friction coefficient, which are given by and the ground friction coefficient, which are given by
FD = µwFV, (9) 𝐹 =𝜇 𝐹 , (9)
where µw is the ground friction coefficient with typical values ranges from 0.4 to 0.9, which depends on tire angular velocity, tire-ground contact pressure, and runway condition.
2.2. Bearing Normal Force Analysis The effects of ground forces on the absorber strut contain the normal force on the bearing area between the outer cylinder and the piston rod, and the strut elastic deforma- tion. Additional bending moments will be generated at the bearing area to ensure that the deformation of the outer cylinder and the piston rod remains consistent within the supporting area. It is embodied by the reaction force on the upper and lower surfaces of the bearing, which increases the total normal force on the bearing area. Figure3 is a schematic diagram of the internal force analysis for the half-axle main landing gear in the xoz plane and the yoz plane under the ground vertical force and longitudinal force. Appl. Sci. 2021,, 11,, xx FORFOR PEERPEER REVIEWREVIEW 5 of 18
where 𝜇 isis thethe groundground frictionfriction coefficientcoefficient withwith typicaltypical valuesvalues rangesranges fromfrom 0.40.4 toto 0.9,0.9, which depends on tire angular velocity, tire-ground contact pressure, and runway condi- tion.tion.
2.2. Bearing Normal Force Analysis The effects of ground forces on the absorber strut contain the normal force on the bearing area between the outer cylinder and the piston rod, and the strut elastic defor- mation. Additional bending moments will be generated at the bearing area to ensure that thethe deformationdeformation ofof thethe outerouter cylindercylinder andand thethe pistonpiston rodrod remainsremains consistentconsistent withinwithin thethe supporting area. It is embodied by the reaction force on the upper and lower surfaces of
Appl. Sci. 2021, 11, 5445 thethe bearing,bearing, whichwhich increasesincreases thethe totaltotal normalnormal forceforce onon thethe bearingbearing area.area. FigureFigure 33 isis aa5 sche-sche- of 18 matic diagram of the internal force analysis for the half-axle main landing gear in the xoz plane and the yoz planeplane underunder thethe groundground verticalvertical forceforce andand longitudinallongitudinal force.force.
FigureFigure 3.3. SchematicSchematic ofof mainmain landinglanding geargear (MLG)(MLG) loadload analysis.analysis.
N𝑁 andandand M 𝑀in inin Figure FigureFigure3 denote33 denotedenote the thethe normal normalnormal forces forcforc andes and additional additional bending bending moments, moments, the subscriptsthethe subscriptssubscripts “A” “A”“A” and andand “B” “B” mean“B” meanmean the actingthethe actingacting points popoints pointints pointpoint A and AA andand point pointpoint B, the B,B, superscripts thethe superscriptssuperscripts “x” and“x” “andy” indicate“y” indicate the xoz theplane xoz planeplane andthe andandyoz thetheplane, yoz plane, respectively. respectively.FT is the 𝐹 internal is is thethe internalinternal force acting forceforce onacting the intersectionon the intersection point of point the upper of the andupper lower and torque lower link.torque link. ToTo simplifysimplify the the calculation calculation process, process, the the followingllowing following assumptionsassumptions assumptions areare made aremade made forfor thethe for land-land- the landinginging geargear gear strutstrut strut inin thisthis in this work.work. work. TheThe The outerouter outer cylindercylinder cylinder isis is simplifiedsimplified simplified toto to aa a cantilevercantilever cantilever beam beam model,model,model, andand thethe piston piston rod rod is is simplified simplified to to an an overhanging overhanging beam beam model. model. The The flexible flexible deformation deformation of theof the outer outer cylinder cylinder and theand piston the piston rod is consideredrod is considered under the under small the deformation small deformation assumption, as- whichsumption, uses which the equivalent uses the equivalent area moment area of mo inertiament ofof theinertia simplified of the simplified uniform beam uniform for calculation,beam for calculation, as shown inas Figureshown4 in. Figure 4.
FigureFigure 4.4. SchematicSchematic ofof uniformuniform beambeam withwith equivalentequivalent areaarea momentmoment ofof inertia.inertia.
I𝐼1 andandand I𝐼2 in inin Figure FigureFigure4 44represent representrepresent the thethe equivalent equivalentequivalent area areaarea moment momentmoment of ofof inertia inertiainertia of ofof the thethe outer outerouter cylindercylinder andand thethe pistonpiston rod,rod, respectively.respectively. InIn thethe xozxoz planeplaneplane andandand thethethe yozyoz plane,plane,plane, thethethe outerouterouter cylinder and the piston rod are subjected to the external forces and the internal forces at the bearing area. The deformation and rotation angle of the outer cylinder and the piston rod at the upper and lower bearing area can be calculated using the external force and the contact internal force. The deformation and rotation angle on the piston rod are recorded as wA, 0 wB, θA, and θB. The deformation and rotation angle on the piston rod are written as wA , 0 0 0 wB , θA , and θB , respectively. According to the coordinated relationship of deformation, it is concluded that 0 0 0 0 wA = wA , wB = wB , θA = θA , θB = θB . (10) Appl. Sci. 2021, 11, x FOR PEER REVIEW 6 of 18
cylinder and the piston rod are subjected to the external forces and the internal forces at the bearing area. The deformation and rotation angle of the outer cylinder and the piston rod at the upper and lower bearing area can be calculated using the external force and the contact internal force. The deformation and rotation angle on the piston rod are recorded as 𝑤 , 𝑤 , 𝜃 , and 𝜃 . The deformation and rotation angle on the piston rod are written as 𝑤 , 𝑤 , 𝜃 , and 𝜃 , respectively. According to the coordinated relationship of de- formation, it is concluded that