刘 伟 党西军
(航空工业西飞,西安 710089)
国内大中型民用飞机机翼挂装副油箱构形较少,挂装副油箱时需满足适航条款的相关规定,在副油箱载荷计算过程中,考虑CCAR-25部第25.301(c)条:如果载荷作用下的变形会显著改变外部载荷或内部载荷的分布,必须考虑载荷分布变化的影响[1],是一个相对复杂的问题。
大中型民用类飞机一般将机翼简化为悬臂梁,采用升力线理论确定机翼结构弹性变形的影响[2-8]。《运输类飞机载荷计算程序包》是符合CCAR-25部适航条款,适用于亚音速、大展弦比飞机的通用载荷计算软件,该程序采用上述弹性修正的方法,解决了结构弹性变形对机翼载荷分布的影响,具备计算尾翼、机翼、机身等部件载荷的能力,但不具备机翼挂装副油箱结构载荷的计算能力。
本文在《运输类飞机载荷计算程序包》基础上,主要借鉴机翼结构弹性修正方法,开展副油箱静载荷计算研究,增加副油箱气动力弹性修正计算、惯性力计算、副油箱质量实时更新等功能模块,形成了满足适航条款的机翼外挂物静载荷计算能力,拓展了《运输类飞机载荷计算程序包》计算功能。
求解飞机在不同高度、不同速度、不同重量和不同受载情况下带外挂物飞行时的运动参数,考虑飞机机动引起外挂物自身攻角和侧滑角的变化以及机翼结构弹性修正引起的攻角增量,分别得到副油箱气动载荷与惯性载荷,两者迭加筛选出总载荷。
外挂物风洞测力试验得出外挂物气动力,一般未考虑机翼的真实刚度,为更精确的得出外挂物气动力,应考虑真实机翼的弹性变形对外挂物自身气动迎角的影响,考虑外挂物自身运动及机翼结构弹性变形的气动力计算如下:
[xryrzr]=[xwgywgzwg]-[xcyczc]
(1)
(2)
[Vx·wgVy·wgVz·wg]=[VxVyVz]+[ΔVxΔVyΔVz]
(3)
(4)
(5)
(6)
式中:
xwg、ywg、zwg为外挂质心处坐标;xcyczc为飞机质心处坐标;ωx、ωy、ωz为飞机运动角速度;VxVyVz为飞机运动速度;VxwgVywgVzwg为外挂运动速度;αwg为外挂攻角;αs为机翼弹性变形引起的外挂物当地攻角增量;βwg为外挂侧滑角;δ为飞机相关舵面偏角;dxj为自定义插值函数;T为插值类型;ci_ID、mi_ID为外挂物气动力数据集;Fxq、Fyq、FzqFxq、Fyq、Fzq为外挂物气动力;Mxq、Myq、Mzq为外挂物气动力矩。
在实际的工程计算中,由于弹性变形对惯性载荷的影响较小,可以忽略弹性变形对惯性载荷的影响,将飞机视为刚体计算[9-10],计算公式如下:
Fxw=-(nx+Δnx)Gw
(7)
Fyw=-(ny+Δny)Gw
(8)
Fzw=-(nz+Δnz)Gw
(9)
Mxw=-Ixwεx+(Iyw-Izw)ωyωz+Iyzw(ω2y-ω2z)
+Ixzw(εz+ωxωy)+Ixyw(εy-ωxωz)
(10)
Myw=-Iywεy+(Izw-Ixw)ωxωz+Ixzw(ω2z-ω2x)
+Ixyw(εx+ωyωz)+Iyzw(εz-ωxωy)
(11)
Mzw=-Izwεz+(Ixw-Iyw)ωxωy+Ixyw(ω2x-ω2y)
+Izyw(εy+ωxωz)+Ixzw(εx-ωzωy)
(12)
Δnx=[-(ω2y+ω2z)(xwg-xc)+(ωxωy-εz)
(ywg-yc)+(ωxωz+εy)(zwg-zc)]/g
(13)
Δny=[-(ω2x+ω2z)(ywg-yc)+(ωyωz-εx)
(zwg-zc)+(ωxωy+εz)(xwg-xc)]/g
(14)
Δnz=[-(ω2x+ω2y)(zwg-zc)+(ωxωz-εy)
(xwg-xc)+(ωzωy+εx)(ywg-yc)]/g
(15)
式中:
Gw为外挂物的重量;FxwFywFzw为外挂物惯性力;Mxw、Myw、Mzw为外挂物惯性力矩;Δnx、Δny、Δnz为外挂物引起的附加过载;Ixw、Iyw、Izw为外挂物质量块对自身质心的惯性矩;Ixyw、Iyzw、Ixzw为外挂物质量块对自身质心的惯性积;εx、εy、εz为飞机角加速度;ωx、ωy、ωz为飞机角速度;nx、ny、nz为飞机质心处过载系数。
根据全机质量、燃油顺序、耗油率,计算出当前工况下副油箱的燃油总质量,插值副油箱燃油百分比组成的数据矩阵,计算出燃油质量特性,与副油箱本体合并为副油箱总质量,副油箱典型质量数据矩阵见表1。
表1 副油箱质量数据矩阵
CCAR-25部中没有针对类似机翼挂装副油箱结构的适航条款,计算载荷时,参考机体结构的适航符合性,按照最大包容原则,选取下列条款进行分析和计算,详见表2,其中突风载荷必须通过动态分析方法确定,本文未考虑。
表2 适航条款
按照上述的适航条款确定出过载-速度包线,高度-速度包线,重量构形,与襟翼收上/放下对称机动、襟翼收上滚转机动、偏航机动、发动机停车、地面情况等工况组合成多种载荷计算情况。
副油箱载荷设计严重情况按照载荷需求的不同设置三个载荷点分别选取包线,每一点载荷按照单力素极值和组合力素极值等28个力素选取严重载荷包线。典型副油箱结构示意图见图1。
注:B点为单独副油箱的载荷,用于校核副油箱本体及其连接;C点为副油箱和挂架的载荷,用于校核挂架本体及其连接;D点为副油箱、挂架、过渡梁三者的总载荷,用于校核过渡梁本体及其与机翼连接。图1 外挂物示意图
在充分消化吸收《运输类飞机载荷计算程序包》源代码基础上,借鉴该程序机翼结构弹性修正,数据采用载数据库管理、计算情况采用代码控制等经验,增加副油箱气动弹性修正计算、惯性力计算和副油箱质量计算模块,实现了副油箱静载荷计算,详细计算流程见图2。
以某民机机翼挂装副油箱静载荷计算为例,计算了4 740种载荷工况,考虑了弹性修正和未弹性修正的副油箱静载荷,根据包线筛选的结果,仅列出部分副油箱本体严重工况计算结果,载荷数据见表3和表4。
表3 B点修正载荷
图2 外挂物载荷计算流程
表4 B点未修正载荷
由上计算结果可以看出:考虑弹性修正和未弹性修正副油箱载荷相差在1%以内,载荷变化不大。
经分析弹性修正和未修正载荷变化不大,与该飞机机动性小,机翼刚度大、副油箱气动面小,惯性力占比大等因素有关,后续载荷计算时,可以进行设计说明,不必考虑弹性修正。
将外挂物载荷计算方法程序化,借鉴机翼结构弹性修正、数据采用数据库管理、计算情况采用代码控制等经验,在工程应用成熟的《运输类飞机载荷计算程序包》的基础上,增加了外挂物载荷计算模块,使外挂物载荷计算满足CCAR-25部的适航条款要求,扩展了《运输类飞机载荷计算程序包》的能力。
目前,只能给出副油箱的集中力,还需单独计算副油箱分布载荷,校核副油箱本体,后续需要在该外挂物计算模块上,增加副油箱分布力计算方法,使分布载荷能直接用于CAE软件,校核副油箱本体。