-
0 引言
-
从国内外飞机机型的发展过程可知,为了增加运载能力和降低运输成本[1],研发和使用超大型飞机(包括客机、货机)是必然趋势,美国、欧盟和俄罗斯等技术发达国家相继推出过波音747、A380、C5A等机型。该类超大型飞机起飞重量为300t~400t,而我国目前研制的运-20飞机最大起飞重量为220t,立项中的宽体客机也仅有230t。随着国力增强和航空事业发展,超大型客机的研制必定会提上日程,针对一些关键技术进行预先研究非常必要。
-
就超大型民用飞机而言,除了飞机各部件都采用更大的结构尺寸以外,其主要特点是普遍采用五点式(多轮多支柱)起落架布置。除了机头的前起落架,还在机身和机翼下布置了主起落架,如图1所示。这可以显著减少起落架和机场跑道的载荷,保证机场跑道的通过性[2]。但也对有严格适航规章要求的起落架地面操作载荷设计,提出了新的挑战。
-
民用飞机起落架地面操作载荷,是飞机在跑道上进行滑行、刹车、转弯等地面操作情况时,起落架承受的地面载荷[3]。根据适航条款第25.471条和第25.489条的规定,地面操作情况下飞机在地面载荷、重力和惯性力的作用下处于平衡状态[4],是一个静力学求解问题。三点式飞机的起落架地面载荷计算方法目前已比较成熟,在《飞机设计手册》上有详细描述[5]。中国商飞基于要求编制的算法和软件,已经得到了CAAC(中国民用航空局)的认可,并在中国商飞相关机型的载荷设计和取证过程中使用。五点式飞机增加了两个起落架,引入了两套新求解变量和几何协调要求,原模型和方法不再适用,需要根据其特性研究出新的模型和方法。
-
图1 五点式起落架布局
-
通过文献查阅可知,该技术作为五点式超大型飞机设计的重点和难点,国内目前几乎没有相关的详细设计经验或资料,理论研究也不多,国外由于技术保密等原因公开文献更少。以前受计算软件和分析条件的限制,国内外普遍将多支柱起落架等效成一个起落架或轮胎考虑[6]。比如1970年,美国航空航天局(NASA)在综合考虑模型的复杂性和计算速度的要求后,将机翼和机身上的主起落架等效成一个起落架,建立了波音747的地面模型[7]。类似方法在国内某专利上也有出现。在20世纪80年代以后,随着计算机技术的发展,开始考虑飞机六自由度的地面运动,用数学建模方法开发了飞机的地面运行仿真程序[8-9],而在国内,通过建立动态响应模型,或者利用多体动力学软件,研究着落载荷的文献较多,地面操作载荷的相关研究几乎没有[10] [11] [12]。
-
本文将对基于三点式飞机起落架地面操作载荷计算模型和算法进行扩展,建立五点式飞机的计算模型和算法,并利用参考机型的设计数据进行试算,验证它的有效性和准确性。
-
1 概述
-
起落架地面操作载荷,是飞机在跑道上进行滑行、刹车、转弯等地面操作工况时,地面反馈的载荷,包括航向(Di)、侧向(Si)和垂向(Vi)三个分力。假设单个起落架上各轮胎的受力情况相同,载荷作用于轮胎与地面的虚拟接触点,如图2所示。
-
图2 起落架地面载荷
-
地面操作情况时,假设飞机和起落架结构为刚体,不考虑气动力。每个起落架上的轮胎变形和缓冲支柱压缩行程,由起落架垂直载荷决定。物理模型如图3所示,起落架压缩包括“缓冲器压缩”和“轮胎压缩”,是两个类似弹簧的串联形式。
-
图3 物理模型
-
根据适航条款规定,飞机在重力、惯性力和地面载荷的作用下,处于平衡状态。惯性力由条款规定的过载系数(侧向SLF、航向DLF、垂向VLF)乘以重力获得。航向载荷和侧向载荷,依据条款适用垂向载荷乘以摩擦系数活动,定义轮胎与地面的阻力载荷摩擦系数uj、侧向载荷摩擦系数usj也依据条款要求。
-
参考《飞机设计:第九册》[5],简化后的全机载荷平衡方程如式(1)~式(3)所示。三个方程先后表示:垂直载荷平衡、俯仰力矩平衡和滚转力矩平衡。
-
式中,W表示飞机重力;C i表示轮胎接地点至重心水平航向距离;T i表示轮胎接地点至重心侧向航向距离;ZA表示重心离地高度;V i表示起落架编号为i的垂直载荷;i表示起落架编号,如图2所示;n表示起落架个数,可以等于3或5。
-
任意一组起落架地面载荷对应一种飞机姿态,可用俯仰角、滚转角、重心高度ZA表示,该对应关系用几何协调方程来表示。
-
根据已知的起落架压缩变形公式和全机几何尺寸参数,可求得如式4所表示的几何协调方程。
-
三点式飞机的几何协调方程已知,可见参考文献5。五点式飞机协调方程则根据三点式模型扩展而来,两者的相互关系详见后文介绍。
-
总结可知,起落架地面操作载荷计算,就是已知飞机设计参数(几何协调方程)和工况参数(重量重心、过载系数、摩擦系数),求解一组满足全机载荷平衡方程的起落架垂直载荷。三点式飞机可以通过迭代法求解该问题,五点式飞机的计算方法暂时没有。
-
2 三点式飞机的迭代算法
-
三点式飞机起落架地面操作载荷的迭代算法的如流程图4所示。
-
图4 三点式迭代算法
-
具体的迭代步骤如下:
-
第一步,根据飞机重量W,给定姿态参数和载荷的大概初值:俯仰角α0、滚转角δ0、起落架垂直载荷V i0(i=1、2、3分别表示前、左、右起落架)。
-
第二步,根据V i0和三点式飞机几何协调方程fT计算姿态参数,并计算飞机重心处受到的合力矩。
-
第三步,根据全机平衡条件,更新起落架垂直载荷。为满足滚转力矩平衡,更新左、右主起落架载荷;为满足俯仰力矩平衡,更新前起落架载荷,并按比例更新主起落架载荷;为了满足垂直载荷平衡,再次按比例更新起落架垂直载荷。
-
第四步,迭代收敛判断,根据垂直载荷V i0迭代变化的差值是否满足精度EPS,决定继续返回迭代,还是结束。
-
3 五点式飞机的迭代算法
-
3.1 全机平衡方程
-
五点式飞机的全机平衡方程与三点式基本相同,差别在于增加了两个机身起落架的垂直载荷,其全机平衡方程可参考式(1)~式(3)获得。
-
3.2 几何协调方程
-
建立五点式飞机几何协调方程是解决问题的难点,为此将五点式飞机起落架布局虚拟分解成两个三点式布局,前起落架与机翼起落架构成一个三点式构型,与机身起落架构成另一个三点式构型,分别定义为1T构型和2T构型。分解示意图如图5所示。
-
图5 五点式构型分解示意图
-
五点式飞机的姿态,即可用两个三点式飞机的姿态参数表示,如式(15)、式(16)所示。
-
将式(15)、式(16)合并后,得到五点式构型的几何协调关系公式定义如下。
-
式中:α1,δ1,α2,δ2分别表示两个三点式构型飞机的俯仰角和滚转角;C i、T i为起落架轮胎接地点至重心的水平航向和侧向距离。
-
当两个三点式飞机构型的姿态角相同,并满足全机的平衡方程时,对应的地面载荷即为所取载荷。两个三点式构型姿态相同的这个条件,我们称之为姿态协调条件。
-
3.3 “载荷平衡-姿态协调”迭代算法
-
相比三点式飞机的迭代求解过程,五点式飞机多了两个机身起落架垂直载荷为未知量,但也增加了姿态协调条件。对应三点式的“载荷平衡”条件,建立“载荷平衡-姿态协调”收敛条件,通过建立内嵌的循环迭代即可。如图6所示,五点式的算法外部迭代判断为载荷平衡条件,内嵌迭代判断为姿态协调条件。
-
具体步骤描述如下:
-
第一步,根据飞机重量设定姿态和起落架载荷初值。
-
第二步,使用式(17)计算1T和2T两个构型的姿态参数和几何参数。
-
第三步,判断两个构型的姿态角是否相同(姿态协调)。如果满足精度要求,则进行第四步。否则,更新起落架载荷后返回第二步,直到姿态角相同。
-
第四步,根据更新的几何参数和起落架载荷,使用式(1)~式(3)计算重心处力矩,根据力矩更新分配起落架载荷。
-
第五步,判断载荷变化是否满足精度,即重心力矩是否满足精度(载荷平衡),决定是否退出迭代。
-
3.4 试算和验证
-
该算法对迭代次数要求较高,必须使用计算机实现,为此本文使用Fortran语言根据迭代算法编制了相应的计算软件,参考国外某型飞机的数据,计算得到的结果差异情况如图6所示。其中差异在2%左右,最大五点滑行刹车情况也仅为3.5%,可以基本认为该算法和结果准确,达到了与国外的一致计算水平。
-
图6 五点式飞机迭代算法
-
图7 试算结果差异
-
4 结论
-
本文研究了五点式飞机起落架地面操作载荷的计算方法,在三点式飞机的“载荷平衡”条件迭代算法上,首次提出了把五点式飞机虚拟分解成两个三点式构型的,建立五点式飞机几何协调方程的方法,确定了载荷还应满足的“姿态协调”条件。根据全机平衡方程和几何协调方程,建立了“载荷平衡-姿态协调”的循环迭代算法,并使用Fortran语言编制了计算软件,通过参考机型的设计参数进行了试算验证,表明了该方法准确可靠,具有较大工程价值和意义。
-
参考文献
-
[1] 范玉青,张丽华.超大型复合材料机体部件应用技术的新进展——飞机制造技术的新跨越[J].航空学报,2009,30(3):534-543.
-
[2] 杨肖敏.飞机地面漂浮性评估与分析[D/OL].南京:南京航空航天大学,2002[2020-02-25].https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CMFD&dbname=CMFD9904&filename=2002000413.nh&v=MDcxNDRITE80SHRYTnJKRWJQSVI4ZVgxTHV4 WVM3RGgxVDNxVHJXTTFGckNVUjdxZlplUnFGeXJo VmI3T1YxMjc=.
-
[3] 聂宏,魏小辉.大型民用飞机起落架关键技术[J].南京航空航天大学学报,2008(4):427-432.
-
[4] 中国民用航空规章编委会.中国民用航空规章第25部:运输类飞机适航标准[S].[S.l.:s.n.],2009:45.
-
[5] 飞机设计手册总编委会.飞机设计手册第9册:载荷、强度和刚度[M].北京:航空工业出版社,2002:130-135.
-
[6] 张明.飞机地面动力学若干关键技术研究[D/OL].南京:南京航空航天大学,2009[2020-02-25].https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CDFD&dbname=CDFD0911&filename=1011254021.nh&v=Mjk4MzRFYlBJUjhlWDFMdXhZUzdEaDFUM3 FUcldNMUZyQ1VSN3FmWmVScUZ5M2xVcjdNVkYyNkg3RzlHdEhPcnA=.
-
[7] BOLTON B,ROGERS J.Landing and Taxi Dynamic Loads Analysis of the Boeing Model 747 Airplane[C].AIAA/ASME 9th Structures,Structural Dynamics and Materials Conference.[S.l.:s.n],1968.
-
[8] PI W S,YAMANE J R,SMITH M J C.Generic Aircraft Ground Operation Simulation[C].27th AIAA Structures,Structural Dynamics and Materials Conference.[S.l.:s.n],1986.
-
[9] PRITCHARD J.Overview of landing gear dynamics[J].Journal of Aircraft,2001,38(1):130-137.
-
[10] ZHANG M,NIE H,WEI X H,et al.Modeling and Simulation of Aircraft Anti-skid Braking and Steering Using Co-simulation Method[J].Compel International Journal for Computation and Mathematics in Electrical and Electronic Engineering.2009,28(6):1471-1488.
-
[11] 顾宏斌.飞机地面运行的动力学模型[J].航空学报,2001,22(2):163-167.
-
[12] WEI X H,LIU C L,LIU X Y,et al.Improved model of landing-gear drop dynamics[J].Journal of Aircraft,2014,51(2):695-700.
-
摘要
为了研究五点式起落架布局的民用飞机起落架地面载荷计算方法,根据中国民用航空规章相关的适航条款要求,以三点式布局的起落架地面操作载荷计算模型和迭代算法为基础,建立了五点式的计算模型和迭代算法,使用国外某型飞机的输入数据进行试算和对比,验证了算法的有效性和准确性,完成了国内该领域的前瞻性探索,具有工程价值和意义。
Abstract
To study the ground handling loads calculation method of landing gear for five-struct civil aircraft, based on the model and iterative algorithm of tri-struct configuration according to CCAR25, a calculation model and iterative algorothm for five-struct were created. After calculation and comparison by using the design data for some foreign airplane, the validity and accuracy of the algorithm were proved,and it has fairly great practical value in engineering application.
Keywords
Five-struct ; mutil-struct ; landing gear ; ground load ; ultra-big airplane