文章快速检索  
  高级检索
基于广义线性模型的固体推进剂贮存寿命评估
洪东跑1, 王英华1, 管飞1, 马小兵2    
1. 中国运载火箭技术研究院, 北京 100076;
2. 北京航空航天大学 可靠性与系统工程学院, 北京 100191
摘要:研究了固体推进剂在贮存过程中性能退化对贮存寿命与可靠性的影响,提出了一种基于广义线性模型的固体推进剂贮存寿命与可靠性评估方法.首先结合固体推进剂退化数据的特点,利用位置刻度模型来描述固体推进剂抗拉强度与贮存时间的关系,建立了固体推进剂的贮存寿命模型.然后利用广义线性模型,给出了贮存寿命模型系数的极大似然估计.根据极大似然估计的渐近正态性,利用Fisher信息矩阵,给出了贮存寿命与可靠度近似置信下限.结合算例表明该方法综合利用固体推进剂在不同贮存时间下的性能数据,改善了固体推进剂贮存寿命与可靠性评估精度.
关键词固体推进剂     抗拉强度     贮存寿命     广义线性模型     位置刻度模型    
Storage life assessment for solid propellant based on generalized linear model
HONG Dongpao1 , WANG Yinghua1, GUAN Fei1, MA Xiaobing2     
1. China Academy of Launch Vehicle Technology, Beijing 100076, China;
2. School of Reliability and Systems Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100191, China
Abstract:To describe the influence of performance degradation on storage life and reliability for solid propellant, a method of storage life and reliability assessment based on generalized linear model was proposed. According to the characteristic of solid propellant degradation data, the local scale model was introduced to relate the solid propellant tensile strength to storage time, and the storage life and reliability model for solid propellant was established. The maximum likelihood estimations of the model coefficients were obtained with the generalized linear model for Poisson distribution. Furthermore, the approximated confidence interval for the reliability function was obtained with Fisher information matrix according to the normal law. The instance analysis shows that the performance data under varied storage time was used synthetically to enlarge the sample size and increase the information for storage reliability assessment. Thus, the storage life and reliability assessment precision was improved
Key words: solid propellant     tensile strength     storage life     generalized linear model     location scale model    

固体推进剂是影响发动机贮存寿命与可靠性的主要因素,对固体推进剂的贮存寿命与可靠性进行评估具有重要的意义.传统的寿命与可靠性评估通常是基于产品的失效或寿命终态特征进行的,而随着产品设计、制造水平的提高以及新技术、新材料与新工艺的不断应用,固体推进剂正逐步朝着高可靠、长寿命的方向发展.固体推进剂在允许的时间和成本内难以通过发动机解剖、试车和实弹试射等试验来获得有效的寿命与可靠性数据[1].因此,传统的寿命与可靠性评估方法难以满足固体推进剂的贮存评估要求.然而,在贮存过程中,固体推进剂的外观性能、燃烧性能、力学性能、密度等会随着贮存时间发生变化[2].针对固体推进剂的这一特性,国内外学者致力于研究固体推进剂的非破坏性评估方法[1, 2, 3, 4],这类方法通过固体推进剂老化机理研究其性能随贮存时间的变化规律,从而对其贮存寿命与可靠性进行评估.

固体发动机及固体推进剂平贮件贮存试验均表明固体推进剂的力学性能会随贮存时间不断退化.退化是能够引起产品性能发生变化的一种物理或化学过程,这一变化随着贮存时间逐渐发展[5].当退化量超出规定的阈值时,会导致推进剂失效.通过定期对固体推进剂平贮件的力学性能进行检测,可获得其在贮存过程中的力学性能退化数据.通过对退化数据进行建模分析,描述产品的性能随贮存时间的变化规律,并外推确定产品的失效时间,可在产品失效数据少甚至无失效的情况下,实现对高可靠、长寿命产品的贮存寿命与可靠性评估[6, 7, 8, 9].

本文结合固体推进剂在贮存过程中的失效模式和失效机理,通过研究固体推进剂的力学性能在贮存过程中随贮存时间的变化规律,给出了一种利用性能退化数据的固体推进剂贮存寿命与可靠性评估方法.

1 固体推进剂性能退化

固体推进剂在贮存过程中失效主要是由其贮存环境的影响和内在因素的变化造成的.在贮存过程中,固体推进剂的性能包括力学性能、燃烧性能和密度等会发生退化,并最终导致推进剂失效,其中力学性能的退化主要表现为推进剂抗拉强度下降.在推进剂浇铸成型后,在固化过程中抗拉强度会有一个明显上升过程,稳定后,随着贮存时间缓慢开始下降.抗拉强度下降属于正常老化现象,通常是由于推进剂黏合剂分子结构网络链条断裂,使聚合物分子离解,变为更小的单元,增加了推进剂的流动性,降低了强度[2].

为了获得固体推进剂力学性能退化数据,需要在贮存过程中对推进剂的力学性能进行测试.为了便于对固体推进剂进行性能测试且不破坏固体发动机结构,在实际应用中,通常将推进剂方坯作为平贮件与固体发动机一起贮存.通过对方坯进行性能测试以反映固体推进剂的性能随贮存时间的变化情况.有研究表明受固化条件、贮存条件、应力状态条件等因素影响,方坯预测寿命与发动机推进剂实际寿命之间存在一定差异,但其在贮存过程中失效模式和失效机理是一致的,同时大多数性能参数变化趋势也是相同的[10].鉴于推进剂的长寿命要求,利用方坯预测寿命的偏差在工程应用可接受范围.因此,在国内外工程应用中,通常利用方坯来分析推进剂的贮存性能,同时会采取相应的措施,减少方坯和发动机推进剂贮存相关性能的差异性.

本文选取抗拉强度作为固体推进剂性能研究对象,通过研究其随贮存时间的变化规律,对固体推进剂的贮存寿命进行分析.取n个固体推进剂方坯作为平贮件与固体发动机一起贮存,在贮存过程中,随机选取固体推进剂方坯进行力学性能试验,以测试固体推进剂贮存后的抗拉强度.由于该试验为破坏性试验,故每个方坯只能进行一次试验.根据试验的顺序对n个固体推进剂方坯的退化数据进行排序:

其中tixi分别为第i个方坯的贮存时间和抗拉强度.

2 利用退化数据的可靠性建模

在寿命与可靠性分析中,位置-刻度模型是一种常用的可靠性模型.它已普遍应用于电子元器件、机械产品的疲劳、金属材料断裂及材料老化等寿命试验中[11].利用位置-刻度模型来描述固体推进剂抗拉强度随贮存时间的变化规律.记抗拉强度变量为X,则位置-刻度模型为

其中,μ(t)σ(t)为位置参数和刻度参数;ε为分布函数为G(x)的随机变量,其中G(x)与位置参数及刻度参数无关.对式(2)进行变换可得

在贮存过程中,固体推进剂的抗拉强度会不断退化,即位置参数μ(t)是关于贮存时间t的函数.在工程应用中,通常利用变换线性模型来描述固体推进剂的性能与贮存时间的关系:

其中a,b为待估参数.

由于固体推进剂的生产工艺较为稳定,不同推进剂方坯抗拉强度的一致性较好,而且刻度参数受贮存时间影响较小,可假设σ(t)为常数,记为σ.则式(2)的位置-刻度模型可变换为

3 基于广义线性模型的可靠性评估 3.1 参数估计

利用式(5)位置-刻度模型来描述固体推进剂抗拉强度与贮存时间的关系,假设ε的分布函数G(x)已知,记其密度函数为g(x).对于式(1)的退化数据,令Y=ln X,yi=ln/xi,则Y的密 度函数为则样本y1,y2,…,yn对应的似然函数为

利用数值求解方法,由式(6)可得模型未知参数的极大似然估计.为了增强对退化数据的适应性并改善参数估计效率.利用广义线性模型求解未知参数的极大似然估计.

假设固体推进剂抗拉强度X服从双参数Weibull分布,则Y服从标准极值分布,其分布函数为),密度函数为,代入式(6)可得对数似然函数

σ已知时,可以利用Poisson分布广义线性模型来获得参数ab的极大似然估计[11, 12, 13].

由式(7)的对数似然函数对σ求导有

已知,,可得wi的估计i和ηi的估计i,由式(8)可得

如此循环迭代直到参数估计没有显著变化[14],可得参数估计,.

3.2 贮存寿命与可靠性评估

给定贮存时间t,把和(,)代入式(3)可得固体推进剂的贮存可靠度估计:

给定贮存可靠度R,由(10)可得固体推进剂的贮存寿命估计:

在工程应用中,通常还需要获得寿命与可靠度置信下限.令θ=(a,b,σ)T,由Wald统计量可知,极大似然估计具有渐近正态性:

其中Ι()为Fisher信息矩阵的估计.记h(θ)为θ的函数,由式(12)可得

其中Jh(θ)=,利用式(13)可得h(θ)的近似区间估计[15].记固体推进剂的贮存可靠度为R(x,θ|t),取h(θ)=,给定置信水平1-α,可得h(θ)的置信下限:

其中,σh=;z1-α为标准正态分布的1-α分位点.通过对(14)进行变换,可得固体推进剂在贮存时间t时的贮存可靠度的置信下限:

同理,结合贮存寿命与可靠度之间的单调关系,由式(15)通过迭代可得贮存寿命置信下限tL.

4 数值算例

为研究固体推进剂在贮存过程中的性能变化规律,制作了3个固体推进剂方坯进行贮存,并对方坯进行力学性能试验,如表 1所示.该型固体推进剂的极限抗拉强度为0.55 MPa.

表 1 固体推进剂方坯抗拉强度Table 1 Tensile strength of solid propellant carton
序号贮存年限/a抗拉强度/MPa
110.710
250.680
3150.650

为了验证本文给出方法的合理性与可行性,假设固体推进剂的抗拉强度服从双参数Weibull分布,在表 1的数据基础上,利用数值模拟方法,随机产生一组数据如表 2所示.

表 2 固体推进剂数值模拟数据Table 2 Simulation data of solid propellant
序号贮存年限/a抗拉强度/MPa
110.710,0.703
230.687,0.686
350.681,0.676,0.668
470.670,0.693,0.665
590.663,0.661
6110.659,0.658,0.645
7130.655,0.653,0.649
8150.650,0.644

利用表 2的固体推进剂数值模拟数据,可得参数的极大似然估计=-0.338 2,=-0.032 0,=0.023 7.给定置信水平γ=0.9,由式(15)分别可得固体推进剂在15 a和20 a的贮存可靠性下限0.995 7和0.994 1,同时可得固体推进剂贮存可靠度为0.995的贮存寿命置信下限17.2 a.

5 结 论

1) 利用位置-刻度模型来描述固体推进剂抗拉强度随贮存时间的变化规律,为度量推进剂性能随贮存时间的变化提供了一种可行的方法.

2) 利用广义线性模型,对贮存寿命模型系数估计进行整体统计推断,有效地增强了对性能数据的适应性,提高了贮存寿命与可靠性评估精度.

参考文献
[1] 张昊,庞爱民,彭松.固体推进剂贮存寿命非破坏性评估方法-老化特征参数监测法[J].固体火箭技术,2005,28(4):271-275.Zhang H,Pang A M,Peng S.Nondestructive assessment approaches to storage life of solid propellants(I)-aging characteristic parameter surveillance method[J].Journal of Solid Rocket Technology,2005,28(4):271-275(in Chinese).
Cited By in Cnki (15)
[2] 张仕念,易当详,宋亚男,等.固体推进剂多失效模式相关的贮存可靠性评估[J].固体火箭技术,2007,30(6):525-528.Zhang S N,Yi D X,Song Y N,et al.Evaluation on storage reliability of solid propellant based on correlative failure modes[J].Journal of Solid Rocket Technology,2007,30(6):525-528(in Chinese).
Cited By in Cnki (2) | Click to display the text
[3] Fillerup J,Pritchard R.Service life prediction technology program,ADA397950[R].Edwards:Air Force Research Laboratory,1998.
[4] 杨根,赵永俊,张炜,等.HTPB推进剂贮存期预估模型研究[J].固体火箭技术,2006,29(4):283-285.Yang G,Zhao Y J,Zhang W,et al.Investigation on prediction model for storage life of HTPB propellants[J].Journal of Solid Rocket Technology,2006,29(4):283-285(in Chinese).
Cited By in Cnki (8) | Click to display the text
[5] 尤琦,赵宇,胡广平,等.基于时序模型的加速退化数据可靠性评估[J].系统工程理论与实践,2011,31(2):328-332.You Q,Zhao Y,Hu G P,et al.Reliability assessment using accelerated degradation data based on time series model[J].Systems Engineering Society of China,2011,31(2):328-332(in Chinese).
Cited By in Cnki (12)
[6] 马小兵,王晋忠,赵宇.基于伪寿命分布的退化数据可靠性评估方法[J].系统工程与电子技术,2011,33(1):228-232.Ma X B,Wang J Z,Zhao Y.Reliability assessment using constant-stress accelerated degradation data based on pseudo life distribution[J].Systems Engineering and Electronics,2011,33(1):228-232(in Chinese).
Cited By in Cnki (11) | Click to display the text
[7] 尤琦,赵宇,马小兵.产品性能可靠性评估的时序分析方法[J].北京航空航天大学学报,2009,35(5):644-648.You Q,Zhao Y,Ma X B.Performance reliability assessment for products based on time series analysis[J].Journal of Beijing University of Aeronautics and Astronautics,2009,35(5):644-648(in Chinese).
Cited By in Cnki (13) | Click to display the text
[8] Huang W,Duane L.An alternative degradation reliability modeling approach using maximum likelihood estimation[J].IEEE Transactions on Reliability,2005,54(2):310-317.
Click to display the text
[9] Sun Q,Zhou J L,Zhong Z,et al.Gauss-Poisson joint distribution model for degradation failure[J].IEEE Transactions on Plasma Sciences,2004,32(5):1864-1868.
Click to display the text
[10] 张昊,庞爱民,彭松.方坯药预测寿命与发动机推进剂药柱实际寿命差异研究[J].固体火箭技术,2005,28(1):53-56.Zhang H,Pang A M,Peng S.Investigation on the differences between predicted lifetime of propellant carton and actual lifetime of motor propellant grain[J].Journal of Solid Rocket Technology,2005,28(1):53-56(in Chinese).
Cited By in Cnki (10)
[11] 洪东跑,马小兵,赵宇.利用变环境数据的可靠性灵敏度分析[J].兵工学报,2011,32(3):341-347.Hong D P,Ma X B,Zhao Y.Reliability sensitivity analysis using variable environment test data[J].Acta Armamentarii,2011,32(3):341-347(in Chinese).
Cited By in Cnki (2) | Click to display the text
[12] 洪东跑,马小兵,赵宇.利用变环境数据的Weibull分布可靠性综合评估[J].北京航空航天大学学报,2012,38(11):1487-1491.Hong D P,Ma X B,Zhao Y.Integrated reliability assessment for Weibull distribution using varied environment data[J].Journal of Beijing University of Aeronautics and Astronautics,2012,38(11):1485-1491(in Chinese).
Cited By in Cnki (1) | Click to display the text
[13] Uusipaikka E.Confidence intervals in generalized regression models[M].London:Chapman & Hall/CRC,2009.
[14] 洪东跑,马小兵,赵宇.基于比例风险模型的可靠性综合评估[J].系统工程与电子技术,2010,32(10):2132-2135.Hong D P,Ma X B,Zhao Y.Integrated reliability assessment based on proportional risk model[J].Systems Engineering and Electronics,2010,32(10):2132-2135(in Chinese).
Cited By in Cnki (6) | Click to display the text
[15] Kalbfleisch J D,Prentice R L.The statistical analysis of failure time data[M].Hoboken:John Wiley & Sons,2002.
http://dx.doi.org/10.13700/j.bh.1001-5965.2014.0030
北京航空航天大学主办。
0

文章信息

洪东跑, 王英华, 管飞, 马小兵
HONG Dongpao, WANG Yinghua, GUAN Fei, MA Xiaobing
基于广义线性模型的固体推进剂贮存寿命评估
Storage life assessment for solid propellant based on generalized linear model
北京航空航天大学学报, 2015, 41(1): 29-32
Journal of Beijing University of Aeronautics and Astronsutics, 2015, 41(1): 29-32.
http://dx.doi.org/10.13700/j.bh.1001-5965.2014.0030

文章历史

收稿日期:2014-01-15
收稿日期:2014-04-30
网络出版日期:2014-06-24

相关文章

工作空间