民用航空发动机在翼时间控制方法
**** Hidden Message ***** 第8卷第3期<BR>2008年6月<BR>交通运输工程学报V01.8 No.3<BR>Journal of Traffic and Transportation Engineering Jun·2008<BR>文章编号:1671-1637(2008)03—0028—05<BR>0<BR>民用航空发动机在翼时间控制方法<BR>戎翔,左洪福,张海军<BR>(南京航空航天大学民航学院,江苏南京210016)<BR>摘要:为了通过视情维修策略控制民航发动机在翼时间,分析了发动机维修与在翼时间的特点,<BR>建立了反映在翼时间与状态监测变量和故障强度之间关系的比例强度模型,确立了性能衰退下的<BR>在翼时间控制限与单位时间总成本期望值最小时的最优预防性维修间隔,并进行了CF6—80C2A5<BR>型发动机运营数据的实例验证。结果表明:故障强度小,则在翼时间控制限低;预防性维修与事后<BR>维修成本率会影响最优预防性维修间隔;完备的发动机拆换与在翼历史记录、维修成本、适航与制<BR>造商要求、寿命件限制等数据是提高在翼时间控制准确性的关键。<BR>关键词:航空发动机;视情维修;比例强度模型;在翼时间<BR>中图分类号:V267;V37 文献标识码:A<BR>Control method of time on wing for civil aeroengine<BR>Rong Xiang,Zuo Hong-fu,Zhang Hai—jun<BR>(School of Civil Aviation,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,Jiangsu,China)<BR>Abstract:In order to apply condition based maintenance(CBM)to control civil aeroengine time<BR>on wing(TOW),the features of aeroengine CBM and TOW were discussed,and the proportional<BR>intensity model was established to describe the relation among TOW,condition monitoring<BR>variables and failure intensity.TOW contr01 1imit due to performance deterioration and optimal<BR>preventive maintenance(PM)interval under minimizing expected total maintenance cost per unit<BR>time were implemented.The operation data of a CF6—80C2A5 engine were studied.It is pointed<BR>that the】ess failure intensity iS,the lower TOW controllimit iS;moreover,the cost rate between<BR>PM and corrective maintenance can affect optimal PM interval;sufficient engine operation data<BR>are the keys of accurate TOW control,such as engine removal and on wing history records,<BR>maintenance cost,requirements of airworthiness and manufacturer,and time limits of life limited<BR>parts etc.3 tabs,4 figs,14 refs.<BR>Key words:aeroengine;condition based maintenance;proportional intensity model;time on wing<BR>Author resumes:Rong Xiang(1979一),male,doctoral student,+86—25-84895772—803,rongxiang@<BR>nuaa.edu.cn;Zuo Hong-fu(1959-),male,PhD,professor,+86-25—84891098,FInS@nuaa.edu.cn,<BR>引<BR>目前航空公司对航空发动机广泛采用视情维修<BR>(CBM)策略,根据发动机的气路性能、振动、孔探、滑<BR>油等状态监控信息,做出维修决策,而维修决策对发<BR>动机在翼时间也有很大影响‘1‘引。现在航空公司的发<BR>动机工程师主要是依靠工程经验来控制发动机在翼<BR>时间,通常是由排气温度(exhaust gas temperature,<BR>收稿日期:2007—12—26<BR>基金项目:国家自然科学基金委员会与中国民用航空总局联合项目(60672164);国家863计划项目(2006AA042427)<BR>作者简介:戎翔(1979一),男,江苏丹阳人,南京航空航天火学博士研究生,从事民航维修工程研究。<BR>导师简介:左洪福(1959一),男,湖南茶陵人,南京航空航天大学教授,博士。<BR>吉口<BR>万方数据<BR>第3期戎翔,等:民用航空发动机在翼时间控制方法29<BR>EGT)衰退率来预测性能衰退引起的发动机剩余在翼<BR>时间,该方法主要局限于单一的EGT参数,很少考<BR>虑其他性能参数,更没有考虑可靠性和维修数据,预<BR>测结果不够准确。另外,国内外还有关于发动机寿命<BR>管理的研究,该模型可以处理由不同种类数<BR>据组成及量纲参差不一的变量,来确定系统寿命与<BR>相关协变量间的关系,已在航海、电子、军事等领域<BR>的可靠性或维修分析方面得到了应用。民航发<BR>动机CBM所得的监测、维修与可靠性数据也具有<BR>类似特点,可以选择该模型,因此,本文在分析历史<BR>失效数据基础上,针对CBM策略,以发动机性能参<BR>数为状态协变量来建立Weibull型PI模型,进而研<BR>究性能衰退下的发动机拆换决策控制限和维修成本<BR>最小时的最优预防性维修间隔,再结合适航指令<BR>(airworthiness directive,AD)、服务通告(service<BR>bulletin,SB)和时控件(1ife limited parts,LI.P)等<BR>约束,实现适合航空公司工程实际的民航发动机维<BR>修决策和在翼时间控制。<BR>1发动机视情维修理论<BR>针对发动机故障ElO-lZ3,可采取两类维修活动:预<BR>防性维修(PM)与故障后的修复性维修(corrective<BR>maintenance,CM)。民航发动机视情维修CBM就<BR>是根据发动机状态,进行维修决策,合理安排维修活<BR>动,尽量避免“维修过剩”或者“维修不足”。其中<BR>PM一般可在航线完成,而CM是故障后的拆换<BR>(removal)和送修(shop visit),这些拆换可以是计划<BR>性或者非计划性,因此,视情维修决策包含维修活动<BR>和维修间隔等内容,对在翼时间有很大的影响。<BR>发动机在翼时间(TOW或time since installa—<BR>tion,TSI)就是从某次初始装机在翼到本次在翼的<BR>某时刻之间经历的连续工作时间,单位是飞行小时。<BR>累积在翼时间(time since new,TSN)是发动机出厂<BR>后首次全新装机到目前时刻之间的工作时间,单位<BR>同上。图1中,£Tow。和t,。舢是某台发动机前两次拆<BR>换时的在翼时间,tT。№是第3次装机后的某个在翼<BR>时间,此时仍未拆换;tT啪的含义也类似理解。<BR>从以上分析可看出,有大量的因素影响到发动<BR>机拆换,因此,在翼时间也就与众多外部运行环境和<BR>内部损伤等因素有关。<BR>2发动机在翼时间控制<BR>2.1比例强度模型<BR>2.1.1 可修系统的可靠性建模<BR>可修系统的可靠物<BR>性和维修建模理论主萋<BR>要有更新过程(renewal墓<BR>process)和非齐次Poisson<BR>过程(nonhomogeneous<BR>Poisson process,<BR>NHPP) ,。<BR>df<BR>\上,<BR>式中:N(£)为(o,≠]区间内观察到的故障次数;Z(£)<BR>为失效过程在t时刻的历史信息。<BR>在此基础上,描述失效过程行为的PI模型为<BR>vet x(£)]=矾(£)expEyTx(£)] (2)<BR>式中:X(t)为协变量向量;y为协变量系数向量;<BR>Vo(£)为基本强度函数,也可以看作风险。<BR>‰(f)已知时式(2)为全参数PI模型,否则为半<BR>参数PI模型,前者的基本强度函数的参数和协变量<BR>系数可同时估计得到,参数估计的准确率相对较高,<BR>故本文选择该类模型。<BR>协变量分为外部协变量和内部协变量,外部协<BR>变量是指其变化会引起系统故障征兆变化的变量,<BR>如工作负荷、环境温度等;内部协变量是指可反映系<BR>统故障征兆变化的变量。由此进一步求解系统的可<BR>靠性函数和失效概率密度函数为<BR>, rf 、<BR>REt x(£)]一exp{一I rE.1 X(u)]du} (3)<BR>、J 0 ,<BR>几l x(£)]=出f X(t)]exp{一l出I X(u)3du}(4) 、J 0 ,<BR>式中:vet}x(£)]为考虑协变量后的强度函数。<BR>2.1.2 Weibull过程型PI模型的求解方法<BR>(1)样本符合Weibull分布的检验方法<BR>对于航空发动机来说,基本强度函数可以采用<BR>幂率强度函数(Weibull类型)[“]。Weibull分布的<BR>可靠度函数为<BR>R(£)一1一F(f)一exp『-一f三门(5)<BR>万方数据<BR>30 交通运输工程学报2008生<BR>式中:F(£)为失效概率分布函数;p为形状参数;叩为<BR>尺度参数。对式(5)两边取两次对数,可得到<BR>In{一ln)一卢ln(t)一卢ln(7/) (6)<BR>由式(6)可知,对于时间样本t,绘制In(t)与<BR>In{一ln)之间的函数曲线,若样本服从<BR>Weibull分布,则该曲线应近似地拟合为一条斜率<BR>大于0的直线,直线的斜率即为形状参数的估计值,<BR>其反映了Weibull分布的基本特征,是用于检验样<BR>本数据是否符合Weibull分布规律的最直观方法。<BR>(2)Weibull过程型PI模型的一般形式<BR>Weibull过程型全参数PI模型为<BR>让£I x(£)]一卫f三11 exp (7)<BR>71、吖/<BR>式(7)中未知参数/9、叩、y的估计需要采用极大<BR>似然估计来实现ET]。<BR>(3)似然函数的构造<BR>似然函数的一般形式为<BR>L(/9,'7,y)一IIf(tl,/9,叩,y)屯R(屯,卢,T/,y)1一屯=<BR>Ⅱ口(£;,卢,'7,y)ⅡR(tJ,卢,叩,y) (8)<BR>式中:咒为样本总数;&为截尾的状态量,值为0表示截<BR>尾,为1表示失效;q为观测到失效时间的样品个数。对<BR>于发动机,性能衰退引起的发动机拆换是真正的失效,寿<BR>命件到寿、服务通告SB和适航指令AD、外来物损伤等<BR>拆换影响因素造成的发动机拆换是右截尾数据。<BR>如果协变量个数为P,且相互无关,并都与时间<BR>相关,则对应的PI模型的对数似然函数为<BR>tn。<BR>单位时间内总成本的期望可以表示为<BR>(11)<BR>式中:肚(z)为t内的平均失效数。对于Weibull类型<BR>础)一卜0(洲“=(爿/ (12)<BR>J 、,,<BR>将式(12)代人式(11),有<BR>E(手)=CpM+CcMt即矿(13)<BR>当13>1时,经函数求极值可知,存在如下nM,<BR>使得上式左边最小。TrM一虻丽CpM]9 (14)<BR>若已知维修成本,再将前述的PI模型估计出的<BR>参数代入上式,就可以得到成本最小时的最优预防<BR>性维修间隔,以此为维修决策的推荐值。<BR>2.3其他在翼时间控制因素<BR>适航当局或制造商发布的AD/SB中可能规定<BR>发动机使用了多少时间后要将其拆下做某项维修工<BR>作,则此种拆换因素决定的发动机维修时间为<BR>TAD=min{TAl,TA2,⋯,ni} (15)<BR>式中:TAD为该台发动机受AD/SB影响的最短维修<BR>时间;TA。为第i个AD/SB规定的维修时间。<BR>发动机中时控件(LLP)都有一个设计使用时<BR>间,该发动机受LLP影响的最短维修时间为<BR>万方数据<BR>第3期戎翔,等:民用航空发动机在翼时间控制方法31<BR>TLLP=min{TLl,TI.2,⋯,T“) (16)<BR>式中:T。i为第i个LLP限制的使用时间。<BR>3 实例计算<BR>本文选择某航空公司1999~2004年内CF6-<BR>80c2A5型发动机的45条拆换历史记录(表1),记录<BR>含有拆换时刻的参数:发动机序号、截尾指示量艿、发<BR>动机排气温度偏差值X。、燃油消耗量偏差值五、高<BR>压转子转速偏差值x3、滑油压力偏差值X4、低压转<BR>子振动偏差值x5、高压转子振动偏差值x6、起飞排<BR>气温度裕度值为、起飞外界温度限制值x8。<BR>从Weibull分布。针对,<BR>tTow样本,由上述的检1<BR>验方法,计算结果见薯。<BR>图2,由图2可知,失效早:<BR>时间t与F(t)之间的一一3<BR>对数曲线可近似拟合_4<BR>为一条直线,表明在翼<BR>时间服从WeibuU分布图2 WeibuII检验<BR>规律。Fig.2 Weibull test<BR>然后比较各协变量的相关系数,计算结果见<BR>表2,在建模时须剔除相关性明显的协变量而保留<BR>首先检验在翼时间tTOW变量的基本分布是否服相关系数较小的协变量,使各参数间独立性较好。<BR>裹1发动机拆换历史记录<BR>Tab.1 Engine removal history records<BR>序号tTOW/h 占X1/℃ X2/% X3/% X4/Pa x5/% X6/% X7∥c X8/℃<BR>149 7 095 1 9.122 3.372 1.783 47 492.760 0.314 O.458 9.278 32.746<BR>243 7 801 1 24.369 4.747 2.092 60 193.350 0.028 O.518 —0.162 29.956<BR>594 7 820 0 4.870 2.238 1.479 —29 641.605 O.163 O.871 1.509 30.447<BR>27l 3 282 O 22.688 7.579 2.123 ——7 384.545 0.645 0.992 —10.376 26.997<BR>裹2协变量之间的相关性<BR>Tab.2 Correlation among covariates<BR>协变量X1 X2 Xs X4 X5 X6 X7<BR>X2 0.628<BR>X3 0.164 0.272<BR>X4 —0.069 —0.190 0.354<BR>X5 0.248 0.235 0.223 一O.145<BR>X6 0.082 —0.058 —0.014 0.153 —0.026<BR>X7 —0.733 —0.046 —0.040 0.232 —0.146 0.034<BR>X8 —0.732 一O.465 —0.040 0.232 一O.145 0.037 1.000<BR>本文选择排气温度偏差和高压转子振动偏差这<BR>2个参数为协变量,根据前面介绍的最大似然估计<BR>方法,估计出相关参数,并得到相应回归方程,其中<BR>分,、y。分别为排气温度偏差协变量系数估计值及高<BR>压转子振动偏差协变量系数估计值。<BR>卢22.880,叩一4 767.081<BR>yl一一0.059 3,y。=一1.913 可ItIx(£)]一揣(赢)1‘%。·<BR>exp(--0.059 3X1—1.913X6)(17)<BR>‘3.1考虑故障强度的在翼时间控制<BR>式(17)是整个发动机机队故障强度的体现,将<BR>表1中数据代入式(17),得到对应的最小、最大故障<BR>强度为0.000 011 7、0.000 582 0,这样可以建立发<BR>动机拆换控制限为<BR>fYl=≥:y^xl(£)=15.888—1.8801n(t)<BR>{ 至. ㈣,<BR>【Y2=芝:y^X^(£)一11.979—1.8801n(t)<BR>h‘。。。=——l<BR>图3中两条对数曲线即为发动机拆换控制限,<BR>由图3可知,故障强度小,则在翼时间控制限就越<BR>低,也就意味着控制越保守。将某台发动机的性能<BR>参数X。和X。的数值序列,经过协变量系数的加权<BR>后在图3中画点,点位于下限下方表示在翼安全,位<BR>于上限上方意味性能衰退,需要拆换,位于中间属于<BR>临近失效的过渡状态,需要重点监控。拆换时发动<BR>机所用时间就是在翼时间£TOW。<BR>3.2考虑成本的在翼时间控制<BR>将估计出的模型参数代入式(14),令成本率口<BR>为C,M/ccM,可得到图4的曲线,该曲线反应了预防<BR>性维修成本和事后维修成本对最优预防性维修间隔<BR>图3发动机拆换控制限<BR>Fig.3 Engine removal<BR>control limits<BR>图4成本率影响下的TPM<BR>Fig.4 TeM effeeted by<BR>cost rate<BR>万方数据<BR>32 交通运输工程学报2008年<BR>T,M的影响。预防性维修成本所占比例高,则<BR>大;事后维修成本所占比例高,则TrM小。<BR>3.3在翼时间综合控制与分析<BR>TpM Rong Xiang,Zuo Hong—fu,Zhang Hai-jun,et a1.Research<BR>由表3可知,对于LLP到寿、AD/SB指令等拆<BR>换原因,具有强制性要求,其所规定的时间限就是发<BR>动机的在翼时间;对于性能衰退引起的发动机拆换,<BR>考虑的性能参数较多,在翼时间控制难度较大,利用<BR>PI模型,考虑故障强度或成本的在翼时间具有一定<BR>的柔性;要提高发动机在翼时间控制的精度,一方面<BR>要做好性能参数的收集整理,建立更准确的状态与<BR>机队可靠性模型,另一方面也要更准确把握LLP到<BR>寿等限制。<BR>表3发动机拆换后在翼时间统计<BR>Tab.3 trow<BR>statistics after engine removals<BR>序号拆换时tTOW/h 发动机拆换原因<BR>271 7 721 性能衰退<BR>272 4 656 LLP到寿<BR>275 7 947 LLP到寿<BR>281 6 198 AD指令<BR>582 7 335 性能衰退<BR>601 5 988 性能衰退<BR>608 9 882 性能衰退<BR>4 结语<BR>本文探索了一种新的民航发动机在翼时间控制<BR>方法,该方法考虑了发动机拆换的历史记录、状态监<BR>控参数、发动机可靠性、关键部件寿命限制、制造商<BR>和适航要求等内容,并通过CF6-80C2A5型发动机<BR>为例,验证了该思路的有效性。本文方法控制过程<BR>简单,结果明了,对航空公司发动机维修工作有重要<BR>的实际指导意义。但是,状态变化趋势与在翼时间<BR>预测等问题还有待于进一步考虑到模型中,以期为<BR>发动机寿命管理提供更好的借鉴。<BR>参考文献:<BR>References:<BR>梁剑.基于成本优化的民用航空发动机视情维修决策<BR>研究.南京:南京航空航天大学,2005.<BR>左洪福,张海军,戎翔.基于比例风险模型的航空发动机视<BR>情维修决策.航空动力学报,2006,21(4):716—721.<BR>Zuo Hong-fu,Zhang Hai-jun,Rong Xiang.Condition based<BR>aero-engine maintenance decision method using proportional<BR>hazards model.Journal of Aerospace Power,2006,<BR>21(4):716—721.(in Chinese)<BR>on an integrated model—based CBM decision support system<BR>for aeroengineEqI}Wang Wen-bin,Jia Xi-sheng,Huang Hongzhong.<BR>Proceeding of the First International Conference on MRin—<BR>tenance Engineering.Beijing:Science Press。2006:216—221.<BR>富涛,许春生.在翼航空发动机剩余寿命预测.中国民航<BR>飞行学院学报,2006,17(3):18-21.<BR>Fu Tao,Xu Chun-sheng.Residual life prediction for aero—<BR>engine on wing[-J].Journal of Civil Aviation Flight University<BR>of China,2006,17(3):18-21.(in Chinese)<BR>李建国,叶新农,徐可君.涡扇发动机寿命管理研究rJ].海军<BR>航空工程学院学报,2005,20(4):468-472.<BR>Li Jian—guo,Ye Xin—nong,Xu Ke-jun.Life management re—<BR>search of military turbofan engine.Journal of Naval Aero—<BR>nautical Engineering Institute,2005,20(4):468—472.(in<BR>Chinese)<BR>Vittal S。Hajela P,Joshi A.Review of approaches to gas tur—<BR>bine life management∥AIAA.Proceeding of 10th<BR>AIAA/ISSMO Multidisciplinary Analysis and Optimization<BR>Conference.New York:AIAA,2004:I-11.<BR>Jiang S T,Landers T L,Rhoads T R.Assessment of repair—<BR>able-system reliability using proportional intensity models:a<BR>review.IEEE Transactions on Reliability,2006,55(2):<BR>328—336.<BR>Jiang S T。Landers T L,Rhoads T R.Semi—parametric pro—<BR>portional intensity models robustness for right-censored re—<BR>current failure dataFJ].Reliability Engineering and System<BR>Safety.2005。90(1):91—98.<BR>Percy D F,Kobbacy K A H,Ascher HE.Using proportionalintensities<BR>models to schedule preventive-maintenance inter—<BR>vals.IMA Journal of Mathematics Applied in Business<BR>and Industry,1998,9(3):289—302.<BR>Naeem M,Singh R,Probert D.Consequences of aero-engine<BR>deteriorations for military aircraft.Applied Energy,<BR>2001,70(2):103-133.<BR>.<BR>Naval Research Logistics,2005,52(3):276—284.<BR>∥ASME.Proceeding of Turbo Expo 2005:<BR>50th ASME International Gas Turbine&Aeroengine Teehni—<BR>cal Congress.New York:ASME,2005:1-13.<BR>r13]Weckman G R,Shell R L,Marvel J H.Modeling the relia—<BR>bility of repairable systems in the aviation industry.Com—<BR>puters and Industrial Engineering,2001,40(1/2):51—63.<BR> GUO Huai—rui,Liao Hai—tao,Zhao Wen-biao,et a1.A new<BR>stochastic model for systems under general repairs.IEEE<BR>Transactions on Reliability。2007,56(1):182—188.谢谢
············ 好好学习天天上上回复 1# 航空 的帖子
南京航空航天大学民航学院,江苏南京210016) 先下下来再说 好东西要学习 值得下载,谢谢
页:
[1]