CN109783953A - A new method for calculating the landing dynamic load of aircraft - Google Patents

A new method for calculating the landing dynamic load of aircraft Download PDF

Info

Publication number
CN109783953A
CN109783953A CN201910059439.4A CN201910059439A CN109783953A CN 109783953 A CN109783953 A CN 109783953A CN 201910059439 A CN201910059439 A CN 201910059439A CN 109783953 A CN109783953 A CN 109783953A
Authority
CN
China
Prior art keywords
aircraft
landing
dynamic load
load
gear
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910059439.4A
Other languages
Chinese (zh)
Other versions
CN109783953B (en
Inventor
梁磊
崔云化
李婉
种小雷
王观虎
顾强康
王乐凡
郭天雄
石鑫刚
刘国栋
苏力海
岑国平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Air Force Engineering University of PLA
Original Assignee
Air Force Engineering University of PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Air Force Engineering University of PLA filed Critical Air Force Engineering University of PLA
Priority to CN201910059439.4A priority Critical patent/CN109783953B/en
Publication of CN109783953A publication Critical patent/CN109783953A/en
Application granted granted Critical
Publication of CN109783953B publication Critical patent/CN109783953B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Testing Of Balance (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of calculation methods of type aircraft landing dynamic load, in order to study load response rule of pavement slab under the action of high tire pressure, big quality, short take-off and landing, according to the mechanical characteristic and the characteristics of motion of type aircraft, analyze undercarriage-wheel-road face coupling, and consider that airplane aerodynamic changes, it is established that the aircraft floor kinetic model of 5DOF.Differential equation group is solved by four step Runge-Kutta, studies the road area load variation issue with space at any time.Result of study shows: the road face dynamic load factor that type aircraft main landing gear generates can reach 1.4, and aircraft subsidence velocity is the principal element for influencing landing load.

Description

A kind of calculation method of type aircraft landing dynamic load
Technical field
The present invention relates to aeronautical technologies and airport engineering field, specifically, being related to a kind of type aircraft landing dynamic load Calculation method.
Background technique
With the rapid development of China's aviation and Defence business, large quantities of type aircrafts such as J-20, Y-20, C919, An- 225, J-31 etc., into the people visual field.These type aircrafts often have the characteristics that big weight, high tire pressure, short take-off and landing.By In the change of aircraft parameter drop characteristic with, great variety occurs more in the past for ground load characteristic.And the machine that China is existing In the rigid pavement design specification of field, only with constant dynamic load factor KdTo reflect aircraft to the characteristic of road area load: at road face end Portion when tire pressure q >=1.08MPa, takes Kd=1.25;When q < 1.08MPa, K is takend=1.20, the novel aircraft tire pressure in China is remote Much larger than 1.25MPa.It is this to gain knowledge by elastic system theory and static(al), dynamic load is equivalent to dead load, then analyze The method of face structural stress, aircraft weight is little, tire pressure is small, in the lower situation of subsidence velocity it is substantially reasonable but right Type aircraft is still indicated obviously to lack scientific and applicability with static behavior.Currently, domestic and foreign scholars are dynamic to aircraft floor Mechanics and road (road) face Structural Dynamics all achieve certain research achievement in respective field, but still there are many queries: Such as the ground load measured data of type aircraft is extremely lacked, meets the aircraft and road plane system coupling machine of actual state It manages still unclear.Therefore, research type aircraft load action rule and road face mechanical response, it is fast to be in response to national defence aviation industry Fast development situation and airfield pavement design requirement Basic Problems urgently to be resolved.
Summary of the invention
For type aircraft landing to the load characteristic problem of airfield runway, in analysis type aircraft landing characteristics and aircraft On the basis of each subsystem, aircraft floor dynamics mathematical model is established, considers that different landing conditions are got off the plane and road plane system Coupling mechanism, the present invention propose a kind of calculation method of type aircraft landing dynamic load.
Its technical solution is as follows:
A kind of calculation method of type aircraft landing dynamic load, comprising the following steps:
Step 1, the sliding race stage force analysis of aircraft landing.
Step 2, the aircraft floor kinetic model for establishing 5DOF.
Step 3 establishes mechanical balance equation in each freedom degree direction.
Step 4 solves differential equation group by four step Runge-Kutta.
Further, the results showed that the road face dynamic load factor that type aircraft main landing gear generates can reach 1.4.
Further, aircraft main landing gear active force accounting can reach 98% or more, and nose-gear landing load is smaller.
Further, aircraft subsidence velocity is the principal element for influencing landing load, and landing quality is secondly, tire pressure and pitch angle Road face dynamic load is influenced little.
The present invention provides a kind of calculation methods of aircraft dynamic load simplified, touched briefly on the essentials, and are convenient for engineering staff Researching and designing meets engineering calculation requirement.
Detailed description of the invention
Fig. 1 is aircraft force analysis figure;
Fig. 2 is aircraft floor kinetic model;
Fig. 3 is body nodal point vertical displacement change curve;
Fig. 4 is aircraft pitch angle change curve;
Fig. 5 is coasting distance change curve;
Fig. 6 is nose-gear landing loads change curve;
Fig. 7 is main landing gear landing loads change curve;
Fig. 8 is main landing gear buffer function accounting curve;
Fig. 9 is the changing rule of main landing gear landing load under the influence of landing quality;
Figure 10 is the changing rule of main landing gear dynamic load factor under the influence of landing quality;
Figure 11 is the changing rule of main landing gear dynamic load factor under the influence of tire pressure;
Figure 12 is the changing rule of main landing gear dynamic load factor under the influence of subsidence velocity;
Figure 13 is the changing rule of main landing gear dynamic load factor under the influence of pitch angle;
Figure 14 is the changing rule of nose-gear dynamic load factor under the influence of pitch angle.
Specific embodiment
Technical solution of the present invention is described in more detail with reference to the accompanying drawings and detailed description.
During aircraft landing, aircraft gliding speed can be divided into both vertically as well as horizontally, the subsidence velocity of vertical direction, Shock effect is generated to road face, forms dynamic load;Horizontal velocity makes wheel and road face generate relative motion, forms frictional force, rubs Wipe the size that power size depends on coefficient of friction and vertical load.
Aircraft in alighting run, self gravity W, air drag D, lift L, face main landing gear vertical force Fz1, frictional force Fx1, face nose-gear vertical force Fz2, frictional force Fx2Collective effect under, be decelerated to and normally slide speed Degree, system of then flying off the runway, aircraft force analysis are as shown in Figure 1.
In sliding race on runway, ground faces to be transmitted on body when the active force of aircraft by undercarriage aircraft, therefore Undercarriage has an important influence aircraft floor kinetic characteristics.Undercarriage is mainly by pillar, buffer, wheel system, support Or the part such as extension and retraction system forms, major function is support and buffer function, to improve the stress feelings on aircraft vertical direction Condition.When calculating, Landing Gear System is reduced to undercarriage buffer and wheel system.
The effect of buffer can be reduced to oil liquid damping action power and air spring force in undercarriage.Aero tyre It can be divided into two parts effect, i.e. tire flexibility power and damping force, as shown in Figure 2.Aircraft by fuselage system, main landing gear system and Nose-gear system composition.m0, m1, m2The respectively quality of airframe system, the quality of main landing gear system, nose-gear The quality of system.J0For the rotary inertia of airframe system.s1, s2The respectively center of gravity of airplane is to nose-gear, main gear wheels The distance of axis center.k1, k2, k3, k4Respectively main landing gear suspension coefficient of elasticity, main landing gear tire coefficient of elasticity are preceding Undercarriage suspension coefficient of elasticity, nose-gear tire spring rate.c1, c2, c3, c4Respectively main landing gear suspension Damped coefficient, main landing gear tire damped coefficient, nose-gear suspension damped coefficient, nose-gear tire damped coefficient.θ For airframe pitch angle.z0, z1, z2Respectively airframe center of gravity, main landing gear suspension center of gravity and nose-gear hang center of gravity The vertical displacement at place.
To airframe system m0, formula (1) is obtained by the Rigid Body in Rotation With differential equation:
In formula,For the angular speed of airframe system.For the angular acceleration of airframe system.For The vertical velocity of airframe center of gravity, main landing gear suspension center of gravity and nose-gear suspension center of gravity.Fh1, Fa1It is risen based on respectively Fall frame buffer damping action power and air spring force.Fh2, Fa2Respectively nose landing gear bumper damping action power and air bullet Spring force.
There are the differential equation (2) on fuselage vertical direction:
In formula, x0For the horizontal displacement on vector,For the course speed of airframe center of gravity.Center of gravity is hung for airframe center of gravity, main landing gear and nose-gear hangs the vertical velocity of center of gravity And the course acceleration of airframe center of gravity, g are acceleration of gravity, kLFor comprehensive lift coefficient, k is takenL=0.5CLρ S, wherein Lift factor CL=2 π (θ-θ0), θ0For aerofoil profile zero-lift angle, ρ is atmospheric density, and S is wing area.
To main landing gear suspension m1Analysis, establishes differential equation of motion (3) in the vertical direction:
In formula, Fk1, Ff1The elastic force and damping action power of main landing gear tire respectively.
To nose-gear suspension m2Analysis, establishes differential equation of motion (4) in the vertical direction:
In formula, Fk2, Ff2The elastic force and damping action power of nose-gear tire respectively.
On the sliding race direction of aircraft landing, differential equation of motion (5) are integrally established to aircraft:
In conjunction with the empirical equation (6) of tire spring rate, (7)
k2=34.514+0.387p0+21.816p0z1 (6)
k4=34.514+0.387p0+21.816p0z2 (7)
In formula, p0For the initial tire pressure of aircraft.μ is between longitudinal sliding motion coefficient of friction, with straight skidding rate Relationship Comparison is complicated, can be provided by empirical equation (8):
In formula, SgFor slip rate, is defined as:
VxFor aircraft taxi speed, VωFor tire linear velocity, r is tire radius, and δ is squeegee action amount, and ω is tyre revolution Moving-wire speed.
Simultaneous (1)~(9) differential equation obtains five yuan of differential equation groups (10) of second order:
In formula, sign is sign function:
Then dynamic load of the tire to road face are as follows:
q1, q2For road face incentive action.
Initial Value Problems of One Order System of Ordinary Differential can be turned to by substitution of variable to the initial-value problem of differential equation of higher order It is calculated.
If m=2,3,4 ..., there are m rank Initial Value Problems For Ordinary Differential Equations:
Enable y1=z, y2=z ' ..., ym=z(m-1), then above formula just turns to first-order ordinary differential equation system:
To the Solve problems of first order differential equation system, since quadravalence Runge-Kutta (imperial lattice-library tower) method being capable of basis The required precision in each stage changes step-length, and programs and use and be all easier to, extensive in practical engineering application, therefore the present invention adopts It is solved with the method.
The dynamic differential equation group of landing period is subjected to variable conversion, by five yuan of differential equations of second order in step 2 Group depression of order turns to following formula.It enables:
Take q1, q2It is 0, five yuan of differential equation group conversions of second order in step 1 are as follows:
By numerical value computing differential equation group, numerical solution is obtained, aircraft landing characteristic is as shown in Fig. 3-Fig. 8, each factor shadow Sound gets off the plane dynamic load properties change curve as shown in Fig. 9-Figure 14.
The foregoing is only a preferred embodiment of the present invention, the scope of protection of the present invention is not limited to this, it is any ripe Know those skilled in the art within the technical scope of the present disclosure, the letter for the technical solution that can be become apparent to Altered or equivalence replacement are fallen within the protection scope of the present invention.

Claims (2)

1.一种新型飞机着陆动荷载的计算方法,其特征在于,将四阶Runge-Kutta法应用于求解飞机地面动力学问题,包括以下步骤:1. the calculation method of a novel aircraft landing dynamic load is characterized in that, the fourth-order Runge-Kutta method is applied to solve the aircraft ground dynamics problem, comprises the following steps: 步骤1、飞机着陆滑跑阶段受力分析;Step 1. The force analysis of the aircraft during the landing and rolling phase; 步骤2、建立5自由度飞机地面动力学模型;Step 2. Establish a 5-DOF aircraft ground dynamics model; 步骤3、建立各自由度方向上的力学平衡方程;Step 3. Establish the mechanical balance equation in the direction of each degree of freedom; 步骤4、通过四阶Runge-Kutta法求解微分方程组。Step 4. Solve the system of differential equations by the fourth-order Runge-Kutta method. 2.根据权利要求1所述的新型飞机着陆动荷载的计算方法,其特征在于,得到以下研究结论:2. the calculation method of novel aircraft landing dynamic load according to claim 1, is characterized in that, obtains following research conclusion: ①飞机主起落架产生的道面动载系数能达到1.4以上;①The dynamic load coefficient of the road surface generated by the main landing gear of the aircraft can reach more than 1.4; ②飞机主起落架作用力占比达到98%以上,前起落架着陆荷载较小;②The main landing gear force of the aircraft accounts for more than 98%, and the landing load of the front landing gear is small; ③飞机下沉速度是影响着陆荷载的主要因素,着陆质量其次,胎压与俯仰角对道面动载影响不大。③ The sinking speed of the aircraft is the main factor affecting the landing load, the landing quality is second, and the tire pressure and pitch angle have little effect on the dynamic load of the pavement.
CN201910059439.4A 2019-01-22 2019-01-22 Novel aircraft landing dynamic load calculation method Expired - Fee Related CN109783953B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910059439.4A CN109783953B (en) 2019-01-22 2019-01-22 Novel aircraft landing dynamic load calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910059439.4A CN109783953B (en) 2019-01-22 2019-01-22 Novel aircraft landing dynamic load calculation method

Publications (2)

Publication Number Publication Date
CN109783953A true CN109783953A (en) 2019-05-21
CN109783953B CN109783953B (en) 2023-04-18

Family

ID=66501836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910059439.4A Expired - Fee Related CN109783953B (en) 2019-01-22 2019-01-22 Novel aircraft landing dynamic load calculation method

Country Status (1)

Country Link
CN (1) CN109783953B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110341980A (en) * 2019-07-11 2019-10-18 中国人民解放军空军工程大学 Analysis method of aircraft taking off and landing at plateau
CN112173158A (en) * 2020-09-25 2021-01-05 中国直升机设计研究所 Landing/ship-borne load calculation method for wheeled landing gear helicopter
CN112597641A (en) * 2020-12-10 2021-04-02 上海宇航系统工程研究所 Carrier landing stability optimization method
CN114386157A (en) * 2020-10-16 2022-04-22 中航西飞民用飞机有限责任公司 Damping design method for main landing gear buffer of airplane
CN116167249A (en) * 2023-04-23 2023-05-26 民航机场规划设计研究总院有限公司 Dynamic load calculation method, device and storage medium for aircraft asymmetrical landing

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201104515D0 (en) * 2011-03-17 2011-05-04 Messier Dowty Ltd Method and system for determining friction coefficient µ for an aircraft landing event
US8180504B1 (en) * 2009-05-21 2012-05-15 Nance C Kirk Aircraft landing gear compression rate monitor and method to increase aircraft landing weight limitation
CN103995917A (en) * 2014-04-17 2014-08-20 中国航空工业集团公司沈阳飞机设计研究所 Undercarriage load simulation method for landing impact of full scale aircraft
CN105138805A (en) * 2015-09-29 2015-12-09 中国航空工业集团公司沈阳飞机设计研究所 Load simulation method for cataplane landing gear

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8180504B1 (en) * 2009-05-21 2012-05-15 Nance C Kirk Aircraft landing gear compression rate monitor and method to increase aircraft landing weight limitation
GB201104515D0 (en) * 2011-03-17 2011-05-04 Messier Dowty Ltd Method and system for determining friction coefficient µ for an aircraft landing event
CN103995917A (en) * 2014-04-17 2014-08-20 中国航空工业集团公司沈阳飞机设计研究所 Undercarriage load simulation method for landing impact of full scale aircraft
CN105138805A (en) * 2015-09-29 2015-12-09 中国航空工业集团公司沈阳飞机设计研究所 Load simulation method for cataplane landing gear

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
梁磊等: "基于ADAMS仿真确定飞机着陆道面动荷载", 《西南交通大学学报》 *
牟丹等: "起落架四点布局无人机着陆动力学分析", 《机械设计与制造工程》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110341980A (en) * 2019-07-11 2019-10-18 中国人民解放军空军工程大学 Analysis method of aircraft taking off and landing at plateau
CN110341980B (en) * 2019-07-11 2022-11-01 中国人民解放军空军工程大学 Airplane plateau take-off and landing load reduction use analysis method
CN112173158A (en) * 2020-09-25 2021-01-05 中国直升机设计研究所 Landing/ship-borne load calculation method for wheeled landing gear helicopter
CN114386157A (en) * 2020-10-16 2022-04-22 中航西飞民用飞机有限责任公司 Damping design method for main landing gear buffer of airplane
CN114386157B (en) * 2020-10-16 2024-12-17 中航西飞民用飞机有限责任公司 Damping design method for main landing gear buffer of airplane
CN112597641A (en) * 2020-12-10 2021-04-02 上海宇航系统工程研究所 Carrier landing stability optimization method
CN112597641B (en) * 2020-12-10 2022-07-22 上海宇航系统工程研究所 Carrier landing stability optimization method
CN116167249A (en) * 2023-04-23 2023-05-26 民航机场规划设计研究总院有限公司 Dynamic load calculation method, device and storage medium for aircraft asymmetrical landing

Also Published As

Publication number Publication date
CN109783953B (en) 2023-04-18

Similar Documents

Publication Publication Date Title
CN109783953A (en) A new method for calculating the landing dynamic load of aircraft
Khapane Simulation of asymmetric landing and typical ground maneuvers for large transport aircraft
CN112560191B (en) Turboprop power airplane thrust resistance system with slipstream influence correction and performance analysis method
CN103324778A (en) Ground load determination method of multi-fulcrum airplane
CN109711008A (en) A method for calculating the envelope of the center of gravity of an aircraft
CN102117362B (en) Light airplane horizontal tail design load determination method under slipstream influence
Khapane Gear walk instability studies using flexible multibody dynamics simulation methods in SIMPACK
Xie et al. Linearization method of nonlinear aeroelastic stability for complete aircraft with high-aspect-ratio wings
CN106599466B (en) Method for evaluating road runway flatness based on vibration acceleration
CN113051662A (en) CFD and DATCOM-based pneumatic modeling and performance evaluation method for folding wingtip variant aircraft
CN103577649B (en) The defining method of cargo hold floor load when transport class aircraft cargo drops
Guo et al. Research on aircraft take-off and landing performance based on flight simulation
Deulgaonkar et al. Design and drag analysis of fixed wing unmanned aerial vehicle for high lift
CN106650077B (en) A dynamic response analysis method for an elastic aircraft wake vortex encounter
Dewey et al. III. NACA TN 3466
CN214138956U (en) Aeromagnetic geophysical prospecting unmanned aerial vehicle wing tip mounting device with adjustable posture
CN108108527B (en) Theoretical calculation method for vertical stiffness ratio of aircraft landing gear
Musaj et al. Numerical and experimental investigation of the aerodynamics of an unconventional w-leading edge reversed delta wing in ground effect
Portapas et al. Simulated pilot-in-the-loop testing of handling qualities of the flexible wing aircraft
CN108001706B (en) Large-span aircraft wing elastic deformation calculation method
CN107798153B (en) Method for determining acting force of landing gear on airplane in landing and landing processes of ski-jump deck
CN107330476B (en) A classification method for highway airstrips based on applicable aircraft types
Khapane Simulation of landing gear dynamics using flexible multi-body methods
CN113642093B (en) Separated ejection type aircraft landing gear modeling method
Nikodem et al. Design of a retractable landing gear for the sagitta demonstrator UAV

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230418