DA-42飞机减震器逆向重建与外载荷计算.pdf

第42 卷第1期2024年1月DO1:10.20096/j.xhxb.1008-9233.2024.01.002西安航空学院学报Journal of Xian Aeronautical InstituteVol.42 No.1Jan.2024DA-42飞机减震器逆向重建与外载荷计算牟正，贺强（中国民用航空飞行学院航空工程学院，四川广汉6 18 30 7）摘要：为了逆向重建DA-42飞机起落架减震器，首先通过拆解测量实物，进行减震器的结构模型设计；然后根据飞机总体参数和减震器结构模型，结合经验公式和Adams软件开展减震器载荷的理论计算和仿真分析，得到飞机在三点着陆工况、俯仰角为2 6 10 及与各俯仰角分别对应的右滚转角为0 1、2 3的工况下的载荷结果及减震器最大轴向力变化规律；最后通过与整机实验的对比证明了计算结果的可靠性。研究结果为DA-42飞机起落架减震器国产化提供一定的支持。关键词：DA-42飞机；减震器；起落架；Adams中图分类号：V226Reverse Reconstruction and External Load Calculation of DA-42文献标识码：AAircraft Shock Absorbers文章编号：10 0 8-9 2 33(2 0 2 4)0 1-0 0 0 7-0 7MU Zheng,HE Qiang(College of Aviation Engineering,Civil Aviation Flight University of China,Guanghan 618307,China)Abstract:In order to reconstruct the landing gear shock absorber of DA-42 aircraft in reverse,thestructural model of the shock absorber is designed by disassembling the measuring object.Then,according to the overall parameters of the aircraft and the structure model of the shock absorber,combined with the empirical formula and Adams software,the theoretical calculation andsimulation analysis of the shock absorber load are carried out.The load results and the maximumaxial force of the shock absorber are obtained under the three-point landing scenarios,with pitchangles of 2,6,10 and right roll angles of 0,1,2,3 corresponding to each pitch anglerespectively.Finally,the reliability of the calculated results is proved by comparison with the wholemachine experiment.The results provide some support for the localization of DA-42 landing gearshock absorbers.Keywords:DA-42 aircraft;shock absorbers;landing gear;Adams收稿日期：2 0 2 3-0 9-0 5基金项目：四川省通用航空器维修工程技术研究中心资助课题（GAMRC2021YB14)作者简介：牟正（1995一），男，四川成都人，硕士研究生，主要从事适航与安全研究。80引言DA-42飞机是国内中级教练机培训中的主流机型。在DA-42飞机服役过程中,不可避免的会对飞机起落架减震器造成损伤，因此需要对减震器进行换件来保证减震持续适航1，而减震器这类航材被国外原始设备制造商（OriginalEquipmentManufacturer,OEM)垄断，购买时会出现采购周期长、金额高甚至难以购买的情况，开展DA-42飞机起落架减震器国产化替代对于保证我国通用航空高质量发展具有重要意义。在国产化替代过程中，计算DA-42飞机减震器外载荷是设计的先决条件2 。减震器是一个包含油、气、橡胶和不锈钢等多种材料的系统，在机身和地面之间承接着多种工况下的受力，这使其载荷的计算复杂多变。目前，国内外对此研究很多采用有限元分析、动力学仿真和落震试验等方法。Sonowal等人3 分别对不同金属材料的减震器进行受力分析和有限元仿真，得知高强度不锈钢模型具有更小的应力和变形。Ark等人4根据减震器冲程速度和飞机下沉速度，设计减震器模糊控制器系统，旨在减少飞机因着陆冲击而产生的振动。Kang等人5 提出了一种压力损失的磁流变减震器设计模型，并以某机型减震器参数设计了磁流变减震器原型，最后结合落震试验，研究了磁流变减震器工作时的运动学特性。Ahmad等人6 根据适航标准，对一个全质量为16 0 0 kg飞机的前起落架减震器组件进行设计，并通过有限元分析对减震器的初步设计进行结构完整性验证。Pecora7在Matlab环境中设计并实现了一套关联程序，整合减震器结构参数与初始条件到动态冲击方程，并模拟落震试验，得到与真实试验相符的结果。祝世兴等人8 1通过建立动力学数学模型与流固耦合有限元模型，对比分析了不同激励下减震器的动力学特性。齐浩等人9建立数字化样机，对无人机减震器在着陆过程中进行动态特性研究。吴卫国等人10 1使用随机模拟方法，结合下沉速度、下沉高度和路面不平度等因素，分析了飞机着陆起落架减震器的随机响应特征。王永庆等人11 对舰载机滑跃起飞过程中减震器的动力学特性进行研究，并结合甲板风场、发动机性能、起落架重量等若干影响因素进行讨论，得出设计所需重点考虑的关键因素。赵知辛等人12 1基于变密度法和优化准则法对某机型前起落架减震器进行拓扑优化设计，在满足西安航空学院学报刚度、强度的条件下实现质量的优化。周进13 利用LMS Virtual.Lab构建包含拦阻系统的舰载机虚拟样机，并结合Automation技术，实现了某舰载机着舰段自动化仿真。综上所述，对DA-42飞机减震器外载荷的计算需要通过理论计算、仿真模拟和整机实验来确定。本文通过对DA-42飞机减震器逆向重建，结合经验公式和Adams软件对DA-42飞机减震器在各着陆工况下进行理论计算、仿真计算，并讨论各工况下减震器最大轴向力变化规律，最后与整机试验结果对比验证。1DA-42飞机参数1.1 DA-42总体参数DA-42飞机起落架采用前三点式布置，飞机总体参数见表1。表1飞机总体参数参数参数值最大着陆重量W/kg1999质心高度h/m1.2质心主轮距6/m0.18前主轮距a/m1.735主轮距c/m2.95重力加速度g/(ms-2)9.8则前起落架停机载荷F，和主起落架停机载荷Fm可由下式计算：F,=MgXbaFu=1 Meg(a-b)2a1.2DA-42减震器逆向重建对DA-42飞机减震器进行拆解，测量结构具体尺寸，并根据机械设计手册设置公差，DA-42减震器均为油-气混合型,DA-42飞机减震器结构示意图如图 1 所示。(1)前起落架减震器结构DA-42前起落架减震器主要由外筒、限油筒、内筒、端盖、底座、防扭臂、密封件（包含密封圈、涨圈、刮油圈、防尘圈)等零部件构成,如图1(a)所示。（2）主起落架减震器结构DA-42主起落架减震器主要由内筒、内筒端盖、外筒、外筒端盖、中心螺栓、节流阀、密封件（方形密封圈、备用密封圈、内筒密封圈、外筒密封圈)第42 卷(1)(2)第1期等零部件组成，如图1(b)所示。限油简密封圈外筒内筒涨圈限油筒内筒端盖端盖密封圈限油筒涨圈端盖涨圈刮油圈防尘圈底座密封圈底座(a)图1DA-42飞机减震器结构示意图2DA-42飞机减震器外载荷计算2.1理论计算依据飞机设计手册1414，着陆载荷的计算应依减震器压缩时间、机轮起转时间、起落架自振周期等数据进行当量重量、使用功量、垂直过载的经验公式计算确定。主起落架着陆当量重量Wm由经验公式算得：Wm=0.5W(3)前起落架着陆当量重量W，由下式计算：W.牛Wa+6其中，为轮胎与地面摩擦系数。前起落架减震器系统吸收的使用功量A，和主起落架减震器系统吸收的使用功量Am分别由公式式（5）、公式（6）确定。A.=W.V,2其中,V,为下沉速度,取3.0 5ms-1。起落架垂直过载ny由下式确定。KiK,:An=n.(Su.gu+ou).W.g其中：K1为考虑减震系统填充参数容差影响的放大系数，取1.0 51.1；K为考虑起转与回弹引起支柱变形影响的放大系数，取1.1；为减震系统使用功功量曲线效率系数，取0.6 50.7；Su为减震器的使用行程；qu为使用行程的垂直方向位移传递系数；ou为轮胎压缩量。对于两点水平着陆主起落架垂直载荷Py.m和水平载荷Px.m由下式计算。Py.m=ny.m W.mg牟正，等：DA-42飞机减震器逆向重建与外载荷计算充气阀顶座密封圈方形密封圈备用密封圈防扭臂底座密封圈底座9Px.m=-0.4Py,y.m在飞机俯仰角为姿态下着陆的主起落架垂直顶座载荷Pyl.m和水平载荷Pxl.m由下式计算。内筒Pyl.m=Py.m coso-Px.m sino中心螺栓Pxl.m=Px.mcos0+Py.msino锁止螺母节流阀外简(h)(4)(5)（6)(7)(8)(9)(10)(11)对于三点水平着陆前起落架垂直载荷Py，n 和水平载荷Px.n由下式计算。Py.n=ny.n W,gPx.,=-0.4 Py.n在飞机着陆过程中，由于机轮会吸收一部分能量，减震器所承受的轴向载荷会略小于起落架垂直载荷，考虑减震器设计时应保留一定的安全余量，本文将减震器轴向受力等同于起落架垂直受力计算。2.2仿真计算对于油-气混合的减震器，其工作中的轴向载荷F，主要由空气弹簧力Fa、油液阻尼力F。、结构间摩擦力F和结构限制力F。组成。F,=Fa+F。+Fr+Fs其空气弹簧力F。可以由下式计算15。V。F.=A.(P。(V。-A a Sm其中：A。为减震器压气横截面面积；P。为减震器充气压强；V。为减震器全伸长时空气体积；Sm为减震器离下止点的轴向位移；为多变指数；Patm为大气压强。油液阻尼力F。可由下式计算16 PoAF。=2(C.A.3/s/其中:po为油液密度;A。为压油面积；Ca为油液流量系数；Ah为油孔面积；为内筒相对于外筒速度。结构间摩擦力F，可表示为Ff=Fn+Fr2Fa=mf.TT其中：Fn为密封件摩擦力；m为密封件摩擦系数；f。为密封件接触面正压力；Fz2为库伦摩擦力；b为库伦摩擦系数；对于前减震器,N和Nt分别为内筒对外筒的作用力和内筒对限油筒的作用力，对于主减震器,N.为内筒对外筒作用力,N.为零。减震器参数如表2 所示。(12)(13)(14)-Patm)(15)(16)(17)(18)(19)10前减震器空载氮气压力Po1/MPa主减震器最大伸长量s2/cm主减震器空载氮气压力Po2/MPa油液密度 po/(kg m-3)大气压力Patm/MPa多变指数油液流量系数Ca库伦摩擦系数Mb减震器在设计时会适当保留一些行程余量，在正常情况下，其结构限制力F。0。2.2.1动力学仿真场景建立DA-42动力学仿真场景如图2 所示。将主、前起落架减震器模型在SolidWorks软件中完成DA-42整机装配，利用“Motion插件进行运动算例计算,将结果以“adm”格式导人Adams 软件中，在Adams/View中添加合适的轮胎和大地模型，调整大地坐标系和飞机坐标系一致，如图2(a)所示。选择减震器金属材料为40 CrMnSiMoVa，密封圈材料为聚四氟乙烯，通过减震器中空气弹簧力Fa、油液阻尼力F。、结构间摩擦力Fr和结构限制力F。的计算，设置合理的减震器刚度系数和阻尼系数，为简化计算量，本次仿真采用刚性机体假设，考虑刚体响应中非弹性质量分布影响（如内筒、防扭臂、机轮等），对零件进行真实质量的添加并设置好约束、接触力、重力等参数，如图2(b)。调整飞机姿态以2.5下滑角接近地面，进场速度为130 kmh-l，下沉速度为3.0 5ms-1,模型设定的仿真时间到10 s为止。2.2.2模型仿真与结果分析动力学仿真场景搭建完成后，对飞机着陆工况下可行的最大范围内俯仰角和滚转角进行组合，选择三点着陆工况（Pitch angle d e g r e e s a n d r o llangle0 degrees,Po_Ro）、俯仰角为2、6 10 及与各俯仰角分别对应的右滚转角为0 1、2 3的工况进行动力学仿真。本文仅展示PO_RO、P2 _R 2、P6_R2、P10 _R 2 工况下右主减震器（MainRight）、左主减震器（MainLeft)和前减震器（Nose)的仿真结果，包含减震器轴向载荷（AxialForce）随时间（Time）变化关系（图3（a）、图4（a）、图5（a）、图6（a）、减震器位移（Displacement）随时间（Time）变化关系（图3（b）、图4（b）、图5（b)图6（b)）、减震器功量图(图 7)。西安航空学院学报表2 减震器参数参数前减震器最大伸长量si/cm第42 卷参数值201.6151.98600.101 31.10.890.2(a)(h)图2 DA-42动力学仿真场景在减震器轴向载荷/位移随时间变化的关系中，工况PO_RO（图3）飞机前轮、主轮在0.3s时刻同时接地，因此前、主减震器轴向载荷/位移同时发生变化，在减震器首次压缩达到最大值后开始反方向伸展，经多次反复的压缩伸展，吸收消耗着陆动能，最终在停机载荷/位移附近趋于稳定，从图中可以看出前减震器的停机载荷和位移分别约为2 0 34N和2 2 mm，主减震器的停机载荷和位移分别约为8 787 N和 42 mm。工况P2_R2、P6 _R 2、P10 _R 2 均是以飞机右主轮先接地，紧接着左主轮和前轮接地，因此，右主减震器、左主减震器和前减震器轴向载荷/位移根据飞机着陆姿态不同先后发生变化，同样也在停机载荷/位移附近趋于稳定。在减震器轴向载荷变化中，前减震器（除了工况PO_R0)会出现从一150 N至一32 0 N的载荷，这是由于飞机主轮先着陆，此时前减震器内筒已伸展至下止点，主轮的着陆会使在前起落架的弹簧下质量（轮胎、内筒、端盖、底座、防扭臂、减摆器等)对前减震器产生一定的冲击载荷。在减震器位移变化中,前减震器位移为零的点表示飞机着陆过程中，前轮未触地，前起落架减震器内筒处于下止点位置。减震器功量图（图7)是根据减震器轴向载荷与位移关系绘制而成，图7（a）、图7（b）、图7（c）、图7(d)分别对应工况PO_RO、P2 _R 2、P6 _R 2、P1O _R2的载荷位移关系，由图可知，功量图相对饱满，表明该工况下减震器具有较好的缓冲性能。第1期3000025000A/Q.10200001500010000三5000XT0-5000L0.0 1.0 2.03.04.05.06.0 7.0 8.09.010.0t/s(a)P2_R23500030.0002500020000150001000050000-5.000L0.01.02.03.04.05.06.07.08.09.010.0t/s(a)P6_R235000300002500020000150001000050000-50000.0 1.02.0 3.04.05.06.07.08.0 9.010.0t/s(a)P10_R235000300002500020000150001000050000-5 0000.01.02.03.04.05.06.07.08.09.010.0t/s(a)PO_RO32.500300002750025000N/22500200001750015 00012.50010000750050002.5000-2.500-5000-20牟正，等：DA-42飞机减震器逆向重建与外载荷计算PO-ROPO-RO150-.Nose-Muin图3PO_RO150-NoseMainl.eft.-.MainRight图4P2_R2150-Nosce-Mainlef.MainRight图5P6_R2150NoseMainleftMainRight图6P10_R2MainNose020406080100120140 160减震器位移/mm(a)11.Nose100-Main500-500.01.0 2.03.04.05.06.07.08.09.010.0t/s(b)P2_R2100500-500.01.0 2.03.04.05.06.07.08.09.010.0t/s(b)P6_R2100500-500.0 1.02.0 3.04.05.06.07.08.09.010.0t/s(b)P10_R2100500-500.01.02.03.04.0 5.06.07.08.09.0 10.0t/s(b)P2_R235000.-MainLeft32.50030000275002500022.50020000175001500012.500100007500500025000-2.500-5000Nose.MainkRightMaihlef-NosC-MainLeft.MainRight-Nose-Mainl.eft.MainRight.-MainRightNose02040减震器位移/mm(b)608010012012N/22.50020 00017 50015 00012.50010000750050002.5000-2.500-10010203040506070减震器位移/mm(c)2.3各着陆工况减震器最大轴向力分析基于仿真数据，对DA-42飞机在各着陆工况下减震器最大轴向力进行整理分析，各工况着陆减震器最大轴向力汇总如图8 所示，可以得到以下结论：(1)DA-42飞机在较小的俯仰角下着陆时,前减震器所受最大轴向载荷高于主减震器。(2)对于同一俯仰角下不同的滚转角着陆工况，前减震器轴向受力基本保持一致。(3)随着俯仰角度的增加，前减震器最大的轴向受力呈下降趋势，而主减震器最大轴向受力则呈上升趋势。(4)在单轮着陆的情况下，首先触地的主减震器轴向受力会高于另一侧的主减震器，但其受力与对应俯仰角的两点水平着陆主减震器受力总体相近。这符合CCAR23规章2 3.48 3条的要求，原因是单轮着陆在未达到极限载荷之前，会产生一个很大的恢复力矩，使飞机倾向于另一侧，变成两点水平着陆情况。3500032.500号30000275002500022.50020000175001500012500100007.5005000R0R1R2R3R0R1R2R3R0R1R2R3R0P2PO着陆工况图8各工况着陆减震器最大轴向力汇总西安航空学院学报P6_R235000.MainLeft32.500300002750025000第42 卷P10_R2.MainRightNose8090100 110图7减震器功量图2.4试验结果对比通过对DA-42飞机在P2_RO工况下进行整机试验，试验所测数据与理论计算、仿真计算的结果误差小于10%，三种方式的结果对比情况如图9所示。30 00025 00020 0001500010000-仿真计算理论计算试验结果5000PO_RO3100030.00029000280002700026 00025 000-仿真计算理论计算试验结果PO_ROP2_ROMainRightMainTeftVNoseP6P1035000F32.50030000275002500022.50020000175001500012.50010000750050002.5000-2.500-50003结论本文面对DA-42飞机减震器国产化需求，通过理论、仿真的方法分别对飞机在三点着陆工况、俯仰角为2、6 10 及与各俯仰角分别对应的右滚转角为0、1、2 3的工况下减震器轴向载荷进行计算，并将计算结果与整机试验数据对比，验证了结果的可靠性，为DA-42飞机起落架减震器国产化提供一定的支持。.-MainLeft020图9对比情况-MainRight4060减震器位移/mm(d)P2_RO卷陆上况Nose80100P6_ROP6_RO120PIO_ROPIO_RO第1期1蒋启登.基于结构分区与应变电桥解耦的支柱式主起落架载荷测量J.北京航空航天大学学报，2 0 2 1，47(5):919-927.2何绪飞，艾剑良，宋智桃.民机起落架摆振仿真与虚拟适航验证J.机械工程学报，2 0 18，54（14）：17 9-18 4.3SONOWAL P,PANDEY K M,SHARMA K K.Design and static analysis of landing gear shockabsorber of commercial aircraftJ.Materials Today:Proceedings,2021,45:6712-6717.4 ARK A E,BILGI B.Fuzzy control of a landing gearsystem with oleo-pneumatic shock absorber to reduceaircraft vibrations by landing impacts J.AircraftEngineering and Aerospace Technology,2023,95(2):284-291.5KANG B H,HWANG J H,CHOI S B.A new designmodel of an MR shock absorber for aircraft landinggear systems considering major and minor pressurelosses:experimental validation J.Applied Sciences,2021,11(17):7895.6AHMAD M A,SHAH S I A,SHAMS T A,et al.Comprehensive design of an oleo-pneumatic noselanding gear strutJ.Proceedings of the Institution ofMechanical Engineers,Part G:Journal of AerospaceEngineering,2021,235(12):1605-1622.7PECORA R.Arational numerical method forsimulation of drop-impact dynamics of oleo-pneumatic牟正，等：DA-42飞机减震器逆向重建与外载荷计算参考文献13landing gearJ.Applied Sciences,2021,11(9):4136.8祝世兴，李佩窈，祝恒佳，等，飞机起落架双腔式油气缓冲器非线性动力学特性研究J.机床与液压，2 0 2 1,49(24):31-36.9齐浩，王泽河，朱华娟，等。飞机起落架落震动力学建模及仿真分析J.机床与液压，2 0 2 1，49（8）：141-146.10吴卫国，孙建桥，冷永刚，等.飞机起落架动力学建模及着陆随机响应分析J.航空学报，2 0 16，37（4）：1228-1239.11王永庆，于浩，施岩.舰载机滑跃起飞动力学与运动学特性J.航空学报，2 0 2 1，42（8）：158-16 9.12赵知辛，王琨，汪杰，等.飞机起落架的动力学分析与拓扑优化研究.机械设计与制造，2 0 2 110）：8 1-8 5.13周进.舰载机起降段设计使用载荷谱编制及载荷谱加重研究D.南京：南京航空航天大学，2 0 19.14高泽迥.飞机设计手册第14册起飞着陆系统设计M.北京：航空工业出版社，2 0 0 2：6 0 5-6 15.15 DING YW,WEI XH,NIEH,et al.Dischargecoefficient calculation method of landing gear shockabsorber and its influence on drop dynamics J.Journal of Vibroengineering,2018,20(7):2550-2562.16 JIAO F J.Oil damping energy loss analysis of landinggear shock absorberJJ.Proceedings of the Institutionof Mechanical Engineers,Part G:Journal of AerospaceEngineering,2019,233(8):3096-3106.责任编辑：吴振松