-
1 引 言
固体火箭发动机是火箭弹、导弹等武器的常用动力装置。随着现代武器远程化战术要求的不断提高，固体推进剂的能量要求和长径比也越来越高。然而，由于固体火箭发动机中装有高能推进剂，其安全问题受到高度关注。近些年来国内外在固体火箭武器试飞过程以及战争中，多次出现由于推进剂药柱结构完整性破坏导致的发动机爆炸事故，造成重大人员财产损失。分析因药柱结构完整性破坏而引起事故的原因，可能是在装药设计中缺乏准确有效的装药结构完整性分析方法，没有精确掌握推进剂材料的损伤失效特性。
国外在研究推进剂等粘弹性材料损伤失效特性方面进展较快，基于累积损伤理论，建立了不同的损伤模型或失效准则来描述固体推进剂等粘弹性材料的损伤失效过程。基于Miner[1]在研究金属材料循环载荷下提出的线性累积损伤理论，Laheru[2]和Bills[3]在研究粘弹性材料受载于恒定载荷下时，认为损伤与时间呈线性关系，提出了失效时间关于应力的幂率形式的损伤失效模型。后来Desmorat等[4]研究认为损伤不再是线性变化，而是以频率关于损伤的函数来表征损伤失效过程。E. J. S. Duncan[5]在研究端羟基聚丁二烯（HTPB）推进剂时，建立了一种包含累积损伤的单轴非线性粘弹性本构方程，并成功应用该方程表征常温下整个应力过程中材料的损伤失效过程。近年，Sullivan[6]利用Kachanov[7]提出的损伤演化定律研究粘弹性材料在循环载荷下的损伤失效情况，并给出损伤参数S的拟合方程。Kunz[8,9]阐述了线性累积损伤模型能很好地预测推进剂在某些载荷条件下的失效情况，提出了适用于Laheru线性累积模型参数确定的方法，以及利用Kachanov[7]损伤模型，得到了固体推进剂在恒定应变率、蠕变等载荷下的损伤失效过程。孟红磊等[10]将E. J. S. Duncan建立的累积损伤模型在常温下应用于双基推进剂，结果显示该累积损伤模型能够很好判断双基推进剂材料的失效情况。史佩[11]利用连续损伤力学理论并耦合线性累积损伤，建立了复合固体推进剂蠕变累积损伤模型，得到了推进剂在蠕变下的损伤失效曲线。梁蔚等[12]在研究循环载荷下HTPB推进剂的疲劳特性时，基于损伤力学理论和粘弹性理论，建立了HTPB推进剂含温度效应的疲劳损伤三阶段演化模型。
综上所述，目前国内外众多学者提出了各种损伤失效模型，认为材料的加载过程，是一个损伤不断累积直至材料出现宏观失效的过程，其失效过程不仅与最大应力、最大应变有关，而且与材料的整个加载应力作用的历史密切相关。然而，这些损伤失效模型中考虑温度变量的却很少，都停留在常温下对推进剂材料损伤失效模型的研究。我国疆域广阔，地域之间温差较大，固体推进剂在服役期间长期承受温度载荷的作用，经研究发现[13,14]外部环境温度严重影响药柱的疲劳失效特性。因此,本研究拟针对HTPB推进剂在不同温度和不同应变率下的失效特性进行探索，基于E. J. S. Duncan[5]提出的累积损伤模型，建立一种含时温等效因子的失效准则，为HTPB推进剂装药结构完整性分析计算和装药结构设计以及推进剂的储存提供必要的依据。
-
2 实验方法与实验结果
-
2.1 实验试件
实验对象为HTPB推进剂，主要成分为氧化剂高氯酸铵（AP）、粘合剂HTPB、铝粉、黑索今（RDX）等，试件形状为哑铃型板条状，其尺寸如图1所示。
-
2.2 应力松弛实验
首先将材料试样密封包装后进行保温3 h左右，待保温箱中温度达到实验所需温度时，利用QJ211B型电子万能实验机进行材料的应力松弛实验。按GJB770B-2005方法413.4规定，实验温度选择为233.15，253.15，273.15，293.15，323.15，343.15 K；松弛应力测试时间：2,4,8,20,40,80,200,600,1000 s；初始恒定应变：5%。每组温度下应力松弛实验重复3次，结果取平均值。其中HTPB推进剂在293.15 K下的应力松弛实验曲线如图2所示。
-
2.3 单轴等速拉伸实验
与应力松弛实验同样的保温方法，在低温233.15，253.15，273.15 K下进行2，20，100，500 mm·min-1的等速拉伸破坏试验和在高温293.15，323.15，343.15 K进行0.5，2，20，100 mm·min-1的等速拉伸破坏试验，同样利用QJ211B型电子万能实验机记录整个实验过程中试件的应力和应变随时间的变化情况，同组实验重复3次，取平均值为当组最终结果。其中HTPB推进剂在293.15 K温度水平下以不同拉伸速度的应力‑应变曲线如图3所示，HTPB推进剂在不同温度水平下以20 mm·min-1拉伸速度的应力‑应变曲线图4所示。
图3 293.15 K，不同拉伸速度下HTPB推进剂的应力‑应变曲线
Fig.3 The stress‑strain curves of HTPB propellant under different tensile speeds at 293.15 K
图4 HTPB推进剂在不同温度于20 mm·min-1拉伸速度下的应力‑应变曲线
Fig.4 The stress‑strain curves of HTPB propellant at tensile speed of 20 mm·min-1 under different temperatures
由图3可见，HTPB推进剂的力学性能明显受应变率影响，应变率越高，屈服应力越高。由图4可知，HTPB推进剂力学特性具有温度相关性，低温下（233.15 K）推进剂最大应力可以达到2.16 MPa，然而相对高温（343.15 K）时，最大应力只有0.68 MPa，可见外部环境温度对HTPB推进剂力学性能影响很大，尤其是在低温下，HTPB推进剂的力学特性在温度之间变化很大。
-
3 时温等效原理及时温等效因子获取
大量研究发现，复合固体推进剂是一种典型的非线性粘弹性聚合物。粘弹性聚合物材料的力学性能的时间相关性和温度相关性之间存在某种转换或等价关系，即时间-温度等效原理[15]，因此时温等效原理也适用于复合固体推进剂。
松弛模量‑温度曲线方法[16]是研究时温等效原理常用方法，即采用同一对数坐标描述不同温度下松弛模量与时间的关系，选取Ts作为参考温度，将对应温度T下的模量曲线随对数时间轴平移至与参考温度Ts重合位置，需要移动的量记为logαT，其数学关系式如下：
即
式中,
由于温度变化会引起材料松弛模量和密度的变化，而模量又随单位体积内所含物质的多少而变化，考虑到这些因素，因此需要对式（2）作出ρT/ρsTs的修正(ρ表示实验温度（T）下材料密度，ρs表示参考温度（Ts）下的材料密度)。但是对于HTPB推进剂材料而言，温度引起材料密度变化量非常小，因此可以忽略不计，式(2)修正后为：
将修正量移至等式左边，取对数，
将不同温度下的松弛应力-时间曲线用对数坐标表示；再根据式（4）取Ts=293.15 K为参考温度，在同一对数坐标系中对不同温度下的模量曲线进行温度修正，得到不同温度水平下的模量曲线，如图5所示。将其他温度所对应的模量曲线平移至与参考温度所对应的模量曲线重叠，即可得到松弛模量主曲线，分别记录不同温度下的模量曲线平移到参考温度下松弛模量曲线的距离值，即为时温等效因子的对数值logαT，重复三次取平均值，得到时温等效因子对数值，见表1。
-
4 基于累积损伤的失效准则研究
-
4.1 失效准则模型
E.J.S.Duncan[5]在研究HTPB推进剂时，提出了一种累积损伤模型，孟红磊[10]也用此模型研究了双基推进剂，该模型表征如下：
式中,σ(t)是应力的时间函数，λ和β是材料的累积损伤系数。此模型认为刚开始材料还没有受力时，D(t)=0,此时材料没有损伤；当力作用材料一段时间后，D(t)=1，此时材料已失效。但是他们均未考虑温度对推进剂失效过程的影响，因此，为了准确表述HTPB推进剂在不同温度下损伤失效过程，需要建立考虑温度及应变率的基于累积损伤模型的失效准则。
HTPB推进剂线粘弹性本构模型的卷积分形式如下[14]：
式中，E(t)为松弛模量，MPa；ξ为折算时间,s；τ为积分变量。其中将松弛模量表示为n阶Prony级数形式，对于等速拉伸过程，应变率保持不变，
（7）
式中，E∞为稳态松弛模量，MPa;n为松弛模量prony级数;Ei则为第i项系数;τi为第i项松弛特征时间，s。根据式（7），可以发现应力与折算时间在HTPB推进剂线粘弹性阶段存在一种函数关系，同时根据许进升[17]关于时温等效原理在HTPB推进剂非线性阶段的适用研究，可以知道HTPB推进剂应力‑物理时间关系同时也可以通过应力‑折算时间关系表示，故基于时温等效原理，引入折算时间ξ，ξ=t/αT，将原模型中积分实际时间更换为折算时间，使模型考虑温度效应，则新失效准则如下：
-
4.2 模型参数获取
-
4.2.1 参数率相关性
孟红磊[18]在研究改性双基推进剂在不同拉伸速率下的屈服和破坏情况时，提出修正型累积损伤模型，即式（9），认为损伤模型参数a与应变率呈幂率关系。
基于孟红磊这种思想，假设失效准则中累积损伤系数λ与应变率的关系也呈幂率关系，即
由此可知在同一温度下，失效准则模型（10）有三个未知数β，a，b。故选用温度293.15 K下三组实验数据（0.5，2，100 mm·min-1）代入式（10）求取模型参数。由于从实验得到的三组实验数据是应力‑时间曲线，通过折算时间与实际物理时间的关系，即ξ=t/αT（αT由第3节求出），可以得到应力‑折算时间曲线，但是仍不便带入式（10）中进行数值计算。为了方便计算，在这里引入一种常用的应力-折算时间函数[19]：
式中，M0为应力时间曲线的初始斜率,MPa·s-1；ξf为材料折算失效时间，s；σf为临界失效应力，MPa。将温度为293.15 K下拉伸速度为0.5， 2, 100 mm·min-1的应力随折算时间变化值用公式（11）拟合，拟合结果如图6所示,分布点表示实验获得的应力随折算时间变化过程，红色实线表示实验数据根据式(11)的形式拟合得到的函数曲线。由图6可见，拟合情况较好，该方程能较好地表示HTPB推进剂在等速拉伸情况下的应力折算时间变化过程。其中各拉伸速度下函数方程（11）的参数M0值、ξf值、σf值见表2。
a. 0.5 mm·min-1
b. 2 mm·min-1
c. 100 mm·min-1
图6 利用函数方程（11）拟合的拉伸试验结果
Fig.6 Fitting results of tensile tests by function equation（11）
表2 不同拉伸速度下函数方程（9）的参数值
Table 2 The parameters of function equation（9）under different tensile speeds
v / mm·min-1 M0 / MPa·s-1 ξf / s σf / MPa 0.5 0.000328 3282 0.66 2 0.00228 868 0.85 100 0.2247 15.4 1.13
v is the tensile speed. M0 is initial slope of the stress‑reduced time curve. ξf is the reduced time. σf is the failure stress.
将表2中各项参数代入式(11)中，可以得到各拉伸速度下应力随折算时间曲线变化曲线的数值方程，将数值方程代入到式（10）中，以失效准则模型数值为1,运用MATLAB数值积分和插值运算方法，反算得到参数β=1.82，从而得到不同拉伸速度下的λ值，将λ的值拟合成关于应变率的幂率函数，拟合曲线如图7所示，其中a=0.69682,b=-0.824,即：
-
4.2.2 参数温度相关性
假设β的值不随温度的变化而变化，认为λ同时也是温度的函数，其它温度下λ的获取方法与293.15 K下相同，再通过同样的拟合方法，以幂率形式表示λ（将温度323.15 K下实验数据用做失效准则验证，在此不拟合），可以得到不同温度下的a，b值，见表3。
表3 不同温度下的a,b值
Table 3 The value of a and b under different temperatures
T / K a b 233.15 0.00077472 -0.77145 253.15 0.0121563 -0.79871 273.15 0.21749 -0.80818 293.15 5.02728 -0.92094 343.15 8.80415 -1.016
利用origin将表3中不同温度下a，b值分别用含温度的函数关系式拟合表示。
当考虑到应变率和温度对HTPB推进剂失效的影响时，失效准则的表达式为：
-
4.3 失效准则验证
根据已经确定的失效准则模型，即式（12），结合323.15 K温度下的2，100 mm·min-1拉伸实验数据，可以得到相应拉伸速度下HTPB推进剂材料损伤演化曲线，如图8所示。
a. 2 mm·min-1
b. 100 mm·min-1
图8 HTPB推进剂在323.15 K下不同拉伸速度的损伤演化曲线
Fig.8 Damage evolution curves of HTPB propellant under different speeds at 323.15 K （The black point is the experimental failure time and the red point is the failure time predicted by failure criterion）
由图8可见，损伤在载荷作用初期，损伤演化很慢，一段时间的损伤模型数值都接近于0，分析其原因为是HTPB推进剂处于弹性阶段，损伤较小，在损伤演化曲线中基本不显示；随后伴随着加载时间的增长，应力不断增大，应变在10%左右时，损伤开始明显出现，这与常武军[20]和职世君等[21]对HTPB推进剂的细观损伤研究，HTPB推进剂的损伤大约发生在其应变10%左右的结果基本吻合；随着加载时间延续，损伤累积增加，直至材料发生失效，此时失效准则模型数值为1。
对于在323.15 K下，以2 mm·min-1的拉伸速度得到的推进剂实验失效折算时间为
拉伸速度为2 mm·min-1时:
拉伸速度为100 mm·min-1时:
通过上述相对误差计算方法，由失效准则预测的HTPB推进剂失效时间的相对误差见表4。
表4 失效准则预测的失效时间相对误差
Table 4 The relative error of failure time predicted by failure criterion
T / K relative error / % T / K relative error / % 2 mm·min-1 20 mm·min-1 100 mm·min-1 500 mm·min-1 0.5 mm·min-1 2 mm·min-1 20 mm·min-1 100 mm·min-1 233.15 1.89 8.10 16.1 3.97 293.15 >20 1.69 2.42 13.1 253.15 8.75 14.1 13.0 18.1 323.15 7.50 15.8 10.1 15.3 273.15 14.3 17.5 14.2 >20 343.15 4.51 9.85 12.7 10.8
由表4可见，由失效准则得到的失效时间相对误差小于20%（除293.15 K下500 mm·min-1和20 ℃下0.5 mm·min-1，分析这两组实验误差较大的原因可能是由于两次拟合造成较大的误差）。说明该基于累积损伤模型的失效准则可以很好地预测HTPB推进剂在不同应变率及温度下的损伤失效过程和失效时间，具备描述HTPB推进剂损伤演化的能力，可作为普遍适用的失效准则用于推进剂的失效预测。
-
5 结 论
结合其他学者关于其他固体推进剂的损伤失效研究，基于累积损伤模型，针对HTPB推进剂建立了考虑温度和应变率因素的失效准则，得到如下结论：
（1） HTPB推进剂的力学性能具有明显的率相关和温度相关性，因此在研究推进剂的损伤失效特性时，必须考虑外界环境温度和应变率对其的影响。
（2） 基于累积损伤模型，并结合不同环境温度和不同拉伸速率的单轴等速拉伸破坏实验以及应力松弛实验，建立了预测HTPB推进剂损伤演化情况和失效时间的失效准则。在323.15 K下，结合2 mm·min-1和100 mm·min-1单轴拉伸实验数据对所建立的失效准则进行验证，结果发现误差低于20%，说明该失效准则能在低温233.15~273.15 K，拉伸速度2~500 mm·min-1和高温293.15~343.15 K，拉伸速度0.5~100 mm·min-1条件下预测HTPB推进剂的失效情况，能够为固体推进剂装药结构完整性分析方法、推进剂失效研究、及推进剂贮存提供一定依据。
-
参考文献
-
1
Miner M A. Cumulative damage in fatigue[J]. J Applied Mechanics Trans Asme,1945, 67：159-164.
-
2
Laheru K L.Development of a generalized failure criterion for viscoelastic materials[J]. Journal of Propulsion & Power, 1992, 8(4):756-759.
-
3
Bills K W. Observations on the linear cumulative damage concept[C]//Proceedings of the JANNAF Structures and Mechanical Behavior Subcommittee Meeting, Chemical Propulsion Information Agency Publication. 1980, 331: 57-59.
-
4
Desmorat R.Damage and fatigue: Continuum damage mechanics modeling for fatigue of materials and structures[J]. Revue Européenne de Génie Civil, 2006, 10(6-7):849-877.
-
5
Duncan E J S, Margetson J.A nonlinear viscoelastic theory for solid rocket propellants based on a cumulative damage approach[J].Propellants,Explosives, Pyrotechnics, 1998, 23(2):94–104.
-
6
Sullivan R W. Development of a viscoelastic continuum damage model for cyclic loading[J]. Mechanics of Time‑Dependent Materials, 2008, 12(4):329-342.
-
7
Kachanov L M. Introduction to continuum damage mechanics[M].Dordrecht: Martinus Nijhoff, 1986.
-
8
Kunz R. Characterization of solid propellant for linear cumulative damage modeling[C]//Aiaa/asme/sae/asee Joint Propulsion Conference & Exhibit. 2013.
-
9
Kunz R. Continuum damage mechanics modeling of solid propellant[C]// Aiaa/asme/sae/asee Joint Propulsion Conference & Exhibit. 2013.
-
10
孟红磊, 赵秀超, 鞠玉涛,等.基于累积损伤的双基推进剂强度准则及实验[J]. 推进技术, 2011, 32(1):109-112
MENG Hong‑lei, ZHAO Xiu‑chao, JU Yu‑tao, et al. Strength criteria and experiments of double‑base propellant based on cumulative damage [J].Propulsion Technology, 2011, 32(1):109-112.
-
11
史佩, 曲凯, 张旭东. 基于连续损伤模型的复合固体推进剂力学性能研究[J].海军航空工程学院学报, 2010, 25(6):662-666.
SHI Pei, QU Kai, ZHANG Xu‑dong.Mechanical properties of composite solid propellants based on continuous damage model[J]. Journal of Naval Aeronautical Engineering Institute, 2010, 25(6):662-666.
-
12
梁蔚, 吕庆山, 陈雄,等.温度对HTPB推进剂疲劳特性的影响[J]. 含能材料, 2017, 25(3):184-190.
LIANG Wei, LÜ Qing‑shan, CHEN Xiong, et al. Influence of temperature on fatigue characteristics of HTPB propellants[J]. Journal of Energetic Materials(Hanneng Cailiao), 2017, 25(3): 184-190.
-
13
Yıldırım H C, Ozupek Š. Structural assessment of a solid propellant rocket motor: Effects of aging and damage[J]. Aerospace Science & Technology, 2011, 15(8):635-641.
-
14
Chyuan S W.Dynamic analysis of solid propellant grains subjected to ignition pressurization loading[J]. Journal of Sound & Vibration, 2003, 268(3):465-483.
-
15
杨挺青.粘弹性力学[M]. 武汉：华中理工大学出版社, 1990.
YANG Ting‑qing. Viscoelastic mechanics[M].Wuhan:Huazhong University of Science and Technology Press, 1990.
-
16
杨挺青.黏弹性理论与应用[M].北京：科学出版社, 2004.
YANG Ting‑qing. Viscoelastic theory and applications[M]. Beijing:Science Press, 2004.
-
17
许进升.复合推进剂热粘弹性本构模型实验及数值仿真研究[D].南京：南京理工大学, 2013.
XU Jin‑sheng. Experimental study and numerical simulation of the thermal viscoelastic constitutive model of composite propellant[D].Nanjing: Nanjing University of Science and Technology, 2013.
-
18
孟红磊.改性双基推进剂装药结构完整性数值仿真方法研究[D]. 南京：南京理工大学, 2011.
MENG Hong‑lei. Research on numerical simulation method for structural integrity of modified double‑base propellant charge[D].Nanjing: Nanjing University of Science and Technology, 2011.
-
19
Saenz L P. Discussion of "equation for the stress‑strain curve of concrete" by Desayi and Krishnan[J]. Aci Journal, 1964, 61:1229-1235.
-
20
常武军, 鞠玉涛, 王蓬勃.HTPB推进剂脱湿与力学性能的相关性研究[J].兵工学报, 2012, 33(3):261-266.
CHANG Wu‑jun, JU Yu‑tao, WANG Peng‑bo. Correlation of dehydration and mechanical properties of HTPB propellants[J]. Acta Armamentari, 2012, 33(3):261-266.
-
21
职世君, 曹付齐, 申志彬,等.复合固体推进剂颗粒脱湿损伤参数反演[J]. 推进技术, 2016, 37(10):1977-1983.
ZHI Shi‑jun, CAO Fu‑qi, SHEN Zhi‑bin, et al. Inversion of dewettingdamage parameters of composite solid propellant particles [J].Propulsion Technology, 2016, 37(10): 1977-1983.
-
1
摘要
为建立考虑温度及应变率效应的端羟基聚丁二烯（HTPB）推进剂失效准则，通过不同温度下的HTPB推进剂应力松弛实验（233.15，253.15，273.15，293.15，323.15，343.15 K）得到了HTPB推进剂的时温等效因子，基于累积损伤理论和线性粘弹性理论，建立了含时温等效因子αT的推进剂失效准则，结合不同温度、不同速率下的单轴拉伸实验数据获取了失效准则损伤参数。利用该失效准则预测了不同温度和拉伸速度下推进剂材料的损伤演化特性和临界失效时间，与实验结果对比分析发现，失效准则预测相对误差低于20%，表明该失效准则能在低温233.15~273.15 K，拉伸速度2~500 mm·min-1和高温293.15~343.15 K，拉伸速度0.5~100 mm·min-1的条件下预测HTPB推进剂的失效情况。
Abstract
To establish a failure criterion for hydroxyl‑terminated polybutadiene (HTPB) propellant which considers the effect of temperature and strain rate. The time‑temperature shift factor of HTPB propellant was obtained by stress relaxation tests at different temperatures (233.15，253.15，273.15，293.15，323.15，343.15 K）of HTPB propellant. Based on cumulative damage theory and linear viscoelastic theory, a failure criterion of propellant with time‑temperature shift factor αT was established. The damage parameters of failure criterion were obtained by combining the uniaxial tensile test data at different temperatures and different strain rates. The failure criterion was used to predict the damage evolution characteristics and critical failure time of propellant material at different temperatures and strain rates. Compared with the experimental results, it is found that the relative error predicted by the failure criterion is less than 20%, indicating that the failure criterion can predict the failure condition of HTPB propellant in the range of low temperature 233.15-273.15 K at tensile speed of 2-500 mm·min-1 and in the range of high temperature 293.15-343.15 K at tensile speed of 0.5-100 mm·min-1.
LI Hui,XU Jin‑sheng,ZHOU Chang‑sheng,et al. Failure Criterion Related to Temperature for HTPB Propellant[J]. Chinese Journal of Energetic Materials（Hanneng Cailiao）,2018,26(9):732-738.
