1、专题研究基于飞参的尾翼载荷谱编制的均值法汤庆辉1,黎泽金,任智勇1,张海东1(1.航空工业洪都,江西南昌,330 0 2 4;2.航空工业成飞,四川 成都,6 1 0 0 9 2)摘要:介绍了基于飞参的尾翼载荷谱编制的均值法,通过对飞参数据的处理获取尾翼载荷谱数据,用均值法编制飞机尾翼载荷谱。此方法编制的尾翼载荷谱可以为飞机尾翼的寿命分析和试验提供载荷谱依据。关键词:飞参;均值法;尾翼载荷谱Mean Method for Compiling Tail Load Spectrum Based on Flight ParametersTang Qinghui,Li Zejin,Ren Zhiyon
2、g,Zhang Haidong(1.AVIC HONGDU,Nanchang,Jiangxi,330024;2.AVIC CHENGFEI,Chengdu,Sichuan,610092)Abstract:The paper introduces the mean method for compiling tail load spectrum based on flight parameters.The tail load spectrum is obtained by processing flight parameter data and the mean method is used to
3、compile the aircraft tail load spectrum.The tail load spectrum compiled using this method can provide a loadspectrum basis for life analysis and testing of aircraft tail.Key words:Flight parameters;Mean method;Tail load spectrum翼载荷谱提供了数据依据。0引言1飞参数据在飞机载荷谱实测结束前,重心过载谱及起落架地谱可以根据规范给定的标准谱进行编制,而规范未提及考核后机身一
4、尾翼结构的尾翼载荷谱的编制。随着飞机气动载荷计算技术的不断发展及研制阶段气动风洞试验数据库的丰富,为根据飞参数据编制尾序号参数名称1相对时间2左垂尾位置3右垂尾位置4左副翼位置5右副翼位置6方向舵位置7左前缘襟翼内侧位置8右前缘襟翼内侧位置9俯仰速率10滚转速率11偏航速率在飞行部队调取了若干个起落的飞机飞参数据,分属于飞行大纲的典型飞行课目。各架次的飞参个数必须满足用于尾翼载荷谱编制。主要的参数可见表1。表1 用于尾翼载荷谱编制的飞行参数单位序号h:ms1213141516。171819o/s20/s21/s22参数名称全机总油量大气机气压高度大气机M数大气机指示空速大气机真侧滑角大气机真攻
5、角左发动机低压转速右发动机低压转速左发动机高压转速右发动机高压转速飞机俯仰角度单位序号kg23m24马赫25km/h2627。28%29%30俯仰角加速度%31%3233参数名称飞机横滚角度飞机磁航向角度飞机轴向加速度飞机法向加速度飞机侧向加速度法向过载表决值侧向过载表决值滚转角加速度偏航角加速度后机身法向过载单位m/sm/sm/s/s2/s2/s2g作者简介:汤庆辉,男,1 9 7 0 年9 月出生,研究员级高级工程师。教练机2 0 2 3.NO.29教练机TRA/NERAVIC数据库或软件,经计算获得飞参采集时刻点的平尾、2尾翼载荷谱数据垂尾和机翼的气动载荷及压心。某型飞机根据飞参计根据表
6、1 的飞参数据,利用气动载荷风洞试验的时间平尾载荷h:m:skg0:11:15-43.140:11:155.740:11:15-17.550:11:15-4.090:11:16-6.04算得出的平尾、垂尾载荷情况见表2。表2 经计算获得的各翼面载荷及压心弦向压心展向压心mm11.2262.14811.8551.63311.1172.21010.5742.86310.7242.606垂尾载荷kg152.62111.1294.5490.3710.29弦向压心m10.46310.29710.27910.3738.642展向压心m2.0162.0372.0382.0312.2952.1平尾、垂尾根部弯
7、矩平尾根部切面弯矩是通过平尾总载和展向压心经计算获得:Mipo=Ppo(Z,-Zo)式中:Mip一第i采集点平尾弯矩;Pp一第i采集点平尾根部切面载荷;Z一第i采集点平尾载荷的展向压心;Zo一平尾根部切面展向坐标。垂尾根部切面弯矩是通过垂尾总载和展向压心经计算获得:Micu=Picw(Y,-Yo)式中:Micw一第i采集点垂尾根部切面弯矩;Picu一第i采集点垂尾载荷;Y一第i采集点垂尾载荷的展向压心;Y。一垂尾根部切面展向坐标。2.2后机身过载当量化飞机飞行过程中,飞参记录到的过载(三向加速度)是对应于当时飞机重量下的过载(三向加速度),为了便于分析及试验,需将随时间和重量变化的过载n当量为
8、某一重量(通常用基本设计重量Go)下的过载nyicq。利用飞参中全机总油量及其他外挂情况可以知道随时间变化的飞机重量Gi,通过各时刻的飞机重量与基本设计重量的比值经计算获得各时刻的当量过载:njieg=G;nyGo3飞参数据预处理飞参数据的预处理主要包括对平尾、垂尾弯矩历程进行去伪码、峰谷值检测和滤波处理等。通过各参数的连续性及各参数变化的协调性来判断参数是否是伪码,然后将伪码进行删除的数据处理为去伪码;峰谷值检测主要是将尾翼弯矩随时间历程处理成满足相邻数据符合谷/峰或峰/谷形式。滤波主要是将尾翼弯矩随时间历程处理成满足相邻数据符合峰/谷差值大于滤波门槛值,滤波门槛值产生的结构应力水平应小于材
9、料的疲劳极限。某型飞机平尾根部切面在承受1 0 0 0 0 Nm的弯矩情况下,平尾根部应力水平为1 5 8 MPa,可以将7 0 0 Nm的弯矩(平尾根部应力为1 1 MPa,小于根部结构材料的疲劳极限)作为滤波的门槛值。4平尾、垂尾弯矩谱4.1平尾、垂尾弯矩历程雨流计数处理对每个架次中的平尾、垂尾弯矩历程,用雨流计数法进行计数处理。某型飞机某课目的平尾弯矩平均计数结果见表3。同理也可获得相应课目垂尾弯矩的平均计数结果。4.2平尾、垂尾弯矩总谱将各典型课目的平尾、垂尾平均弯矩谱与一个训练周期中典型课目的架次数相乘后叠加,便可得到平尾、垂尾弯矩总谱,计算公式:P=BQ式中:Q为课目矩阵;B为各课
10、目平尾、垂尾弯矩平均谱;P为平尾、垂尾弯矩总谱。4.3尾翼弯矩谱超越频次曲线根据飞参得到的平尾的均值非负、均值为负及垂尾左、右面谱可通过材料等寿图(见图1),按等损伤10教练机2 0 2 3.NO.2专题研究表3某型飞机某课目平尾弯矩平均雨流计数结果MaMm-7000-6300-6300-5600-5600-4900-4900-4200-4200-3500-3500-2800-2800-2100-2100-1400-1400-700-700-00-700700-14001400-21002100-28002800-35003500-42004200-49004900-56005600-6300
11、6300-70007000-77007700-84008400-9100原理折算成谷值为零的常谷值谱。折算方法为:PA=2Pa=P,-P.+P.uNPPP图1 材料等寿图式中:P一破坏载荷,根据静力试验结果反推得到;Pm一均值;P一幅值;Pa一折算后的幅值;PA一折算后的峰值。经折算后的平尾常谷值谱见表4(均值为非负)和表5(均值为负),1 0 0 0 飞行小时平尾谱的超越频数701400-0-70014000.20.20.20.21.40.621.6628.63.820.411.829.88205153.810.84.614.67.87.8687.45.644.633.21.81.60.80
12、.40.60.20.22100-28003500-4200-4900-210028000.210.60.81.62.21.40.61.41.62.832.61.20.20.60.60.42xPixPa5600-6300-7000-7700-8400-350042000.80.60.60.40.41.20.20.42.81.21.61PNm49000.20.20.40.61.20.40.20.60.410.20.4560000.20.40.20.40.210.40.40.2曲线见图2;同理可得1 0 0 0 飞行小时垂尾左、右面谱的超越频数曲线,见图3。弯矩/Nm)均值非负01均值为负图2 某型
13、飞机1 0 0 0 飞行小时平尾弯矩谱超越频次曲线弯矩/(Nm)左面谱01右面谱图3某型飞机1 0 0 0 飞行小时垂尾弯矩谱超越频次曲线63000.22270000.23377000.20.40.64484000.60.40.2LOG(NC)56LOG(NC)5691000.4教练机2 0 2 3.NO.211教练机TRA/NERAVIC表4 1 0 0 0 飞行小时平尾弯矩常谷值谱(均值为非负)序号峰值17002140032100谷值常谷峰值0700.0700704.71400709.5频次累积频次序号峰值12598.378240.714968.565642.44739.250673.9谷
14、值常谷峰值467000474200484900频次累积频次7006342.5-35007450.3-28007498.9109.1157.2187.83838.63729.53572.4434445序号峰值102-7003-1400490056006300-1400-7000谷值常谷峰值-700-700.0-1400-704.7-2100-709.56216.76258.16300.0表5 1 0 0 0 飞行小时平尾弯矩常谷值谱(均值为负)频次累积频次序号峰值15232.8315634.531320.0300401.764315.3269081.7287.0461.059.34645.943
15、58.93897.98889904546479100980010500-2800-35002800-8400-7700-7000谷值常谷峰值-7700-5035.0-8400-5069.9-3500-6135.516197.216298.316400.747.0150.920.6频次累积频次31.46155.86.36124.4709.36118.1218.4171.420.6.4243444.4高载截取和低载截除叫由于高载迟滞效应,在载荷时间历程中,最高载荷的选择对于得到一个可靠的裂纹扩展预测将是一个关键。根据以往的经验,一般选取1 0 0 0 飞行小时内超越数为1 0 次的载荷作为最高级载
16、荷。为缩短试验周期,应去除几乎不造成疲劳损伤的低载荷循环。5尾翼载荷前、后压心统计5.1尾翼刚心线根据尾翼有限元模型,通过有限元分析的位移结果进行分析,可以获得尾翼若干个切面的刚心位置。对这些切面压心可以进行拟合并获取拟合刚心线方程。某型飞机垂尾的5 个切面的刚心位置见表6。表中展切面号切面展向位置Y(mm)1切面11352切面15353切面19355.2尾翼载荷前、后压心统计建立局部坐标系:坐标原点为刚心拟合线与飞机对称面(垂尾:水平基准面)的交点;z轴(垂尾y轴)12教练机2 0 2 3.NO.2:-700-5600-1400-6300-2100-7000-4933.1-4966.6-50
17、00.6846.7353.0314.2表6 某型飞机垂尾刚心位置刚心位置X(mm)切面号切面展向位置Ymm)101304切面104405切面10710为刚心拟合线,沿展向为正;y轴向上为正(垂尾z轴向右为正),x轴符合右手定则。将每个采集点的载荷压心坐标转换到局部坐标系下的坐标,并判断其符号,7669.66822.96469.9868788向位置是指到飞机水平基准面的距离,刚心位置是到机头的距离。由5 个切面组成的刚心线如图4 所示,拟合的刚心线方程为:X=0.7131Y+9331.5(R 2=0.9 9 9 6),同理可以获得平尾的刚心线方程。X11400112001100010800106
18、001040010200100000630077007000-8400-8400-91005001000图4 某型飞机垂尾刚心线刚心位置X(mm)233511000278211310-13863.6-14994.4-15088.61500200049.010.314.5R2=0.99962500300073.724.714.5Y专题研究正为后压心,负为前压心。某型飞机垂尾坐标转换公式为:Xjt=(Xo-9331.5)cos0-Yosino后压心课目名称百分比/%课目142.84课目251.61课目345.00平均44.77式中,0 为刚心拟合线与飞机水平基准线的夹角。对所有起落垂尾载荷压心进行
19、统计后,得到各课目的前、后压心比例及平均压心,见表7。表7 各课目垂尾前、后压心比例及平均压心统计结果坐标/mmy坐标/mm11392.92200.311459.12248.911394.92197.511557.12283.9前压心百分比/%x坐标/mm57.1611272.648.3911227.255.0011245.955.23y坐标/mm2081.62099.22051.710902.61990.1尾翼不同载荷与后机身过载的一种组合称之为6尾翼一后机身试验载荷谱确定一种试验载荷状态。后机身过载峰谷值在平尾弯矩谱6.1尾翼载荷谱确定在进行尾翼一后机身试验时,尾翼载荷通过分布载荷施加,载
20、荷压心有很多种不同情况,按飞参中各对应参数获得的载荷压心施加难以实现,有必要进行工程简化。简化原则:1)以弯矩历程为主编制尾翼载荷谱,保证施加的载荷历程产生的根部弯矩历程与飞参数据等效;2)认为扭矩和剪力对构件产生的损伤较小,总体上应与飞参数据相当,有利于试验实施。为便于按弯矩等效确定平尾载荷谱,将前、后压心分别设为(xgo,zo)和(xio,zo),展向压心zo为统计出的前、后压心展向坐标的平均值;同理可分别设定垂尾前、后压心(xgo,yo)和(xio,yo),展向压心yo为统计出的前、后压心展向坐标的平均值。通过弯矩等效转换可以获得一个周期的平尾载荷谱,同理也可以获取垂尾左、右面载荷谱。6
21、.2尾翼、后机身试验谱载荷状态确定由于尾翼载荷谱分别有前、后压心两种情况,可以从飞参数据中挑选出接近于尾翼前、后压心附近的一系列参数用于分析尾翼谱载荷状态。通过这些参数计算出的平尾、垂尾的载荷分布,最后确定平尾和垂尾前、后压心情况下的载荷状态。6.3尾翼、后机身试验谱编排6.3.1尾翼一后机身试验谱课目顺序尾翼、后机身试验谱是按飞一续一飞谱型来编制的,课目顺序可参考重心过载谱的课目顺序。6.3.2各课目内尾翼载荷谱和后机身过载编排统计计数时跟随取出。为编制试验谱,需进行平尾(后机身过载)和垂尾载荷的配对。分别对各课目载荷进行匹配。一个课目按平尾对称谱前、后压心,垂尾左面谱前、后压心,垂尾右面谱
22、前、后压心6 种情况的各级载荷进行组合。垂尾有载荷视为非对称载荷,并有左、右平尾载荷按对后机身严重扭转配对。平尾对称情况无垂尾载荷,并且左、右平尾载荷相同。平尾和垂尾前、后压心比例按统计结果确定。垂尾载荷向左为正,即正载荷产生机身的航后扭矩。左面谱和右面谱的垂尾载荷仅方向相反,大小和频数相同。按以上原则进行编排后,某型飞机某课目平尾对称谱前、后压心载荷和垂尾左、右面谱的前、后压心的载荷组合结果分别见表8 表1 0。表8 某型飞机某课目平尾对称载荷组合号载荷峰载荷谷后机身过后机身过前压心后压心值(N)值(N)144822-7470329884-672313-149414-74715-224116
23、06.4尾翼一后机身试验谱的实施载峰值载谷值频次频次-97124.63-89641.04-82174.06-82171.28-67231.81-67233.38-59762.88-59763.981.270.661.460.831.290.970.510.790.30.20.31.90.20.81.20.70.40.20.42.40.21.01.60.9教练机2 0 2 3.NO.213教练机TRAINERAVIC表9 某型飞机某课目垂尾左面谱载荷组合左平尾序号峰值(N)13735237353373543735左平尾谷值(N)峰值(N)747-2988747-2988747-1494747-1
24、494右平尾谷值(N)-5229-5229-4482-4482右平尾垂尾载荷垂尾载荷后机身峰值(N)谷值(N)0-51361027-61641027-61641027-8218后机身过载峰值过载谷值心频次3.460.683.460.683.460.683.460.68前压2.10.70.010.3后压心频次1.00.30.0030.135363738序号左平尾峰值(N)谷值(N)13735237353373543735.353637381)编排后每种载荷级的频次为小数,同载荷级累计满1 施加1 次;2)非对称垂尾载荷为前压心时,平尾载荷按后压心施加,垂尾载荷为后压心时,平尾载荷按前压心施加;3
25、)载荷级为平尾对称情况时,左、右平尾载荷对称施加;4)平尾、垂尾、后机身载荷同时施加;5)平尾外推载荷按平尾对称情况后压心处理、垂尾外推载荷按后压心处理;6)载荷循环由谷到峰,执行到第i级峰值载荷Si时,则判断可执行的下一级(i+1)的谷值,当Si+1谷S;谷时,载荷退到S,所对应的谷值;否则退到Si1的谷值。14941494149414941494149414941494-747-747-747-747左平尾峰值(N)747-2988747-2988747-2988747-2988-7474482-7474482-7474482-74744824482448244824482表1 0 某型飞
26、机某课目垂尾右面谱载荷组合右平尾右平尾垂尾载荷垂尾载荷后机身谷值(N)峰值(N)-522915409-522910273-52299246-52298218-5229-4109-5229-5136-5229-6164-5229-92467结语依据上述方法所编制的尾翼载荷谱,可以为后机身一尾翼结构关键部位的疲劳分析及模拟件耐久性试验乃至后机身一尾翼组合体耐久性试验提供载荷谱依据,进而为完成飞机后机身一尾翼结构抗疲劳设计和寿命评定奠定基础。参考文献1蒋祖国,田丁栓,周占廷.飞机结构载荷/环境谱M.北京:电子工业出版社,2 0 1 2:1 6 0.-5229-5229-5229-5229102731
27、02731027312327-3082-1027-154099246谷值(N)-10273-5136-4109-1027-7191-8218-9246-123277)试验循环周期之间的衔接均从载荷零开始,最后回到载荷零值。1.411.411.411.41过载峰值过载谷值心频次3.131.893.131.893.131.893.131.894.001.494.001.494.001.491.410.27(收稿日期:2 0 2 3-0 4-2 3)0.270.270.270.27后机身0.40.50.10.1前压0.10.20.30.60.61.10.30.10.20.30.10.1后压心频次0.10.10.10.30.30.50.10.114教练机2 0 2 3.NO.2