循环加载下复合推进剂的非线性热粘弹性本构模型

复合推进剂在动态加载下伴随着自热效应，即存在明显的温度升高现象。为了定量描述这种热力耦合特性，基于不可逆热力学理论建立了非线性热粘弹性本构模型和相应的数值算法，编写了基于Abaqus的用户子程序，并通过实验验证了数值算法的正确性。研究结果发现，非线性本构模型成功地描述了复合推进剂在不同加载模式下的非线性力学行为和自热效应，自热效应对复合推进剂的力学行为有着直接影响，与不考虑自热效应相比，循环加载下应力-应变滞回圈的斜率降低显著，表明力学性能已受到自身温度升高的影响。结果表明，考虑自热效应的热力耦合本构模型能更准确地描述复合推进剂的循环力学行为，滞回圈的形状与实验结果更为吻合，斜率误差在10%以内。     

基于热力学理论的固体推进剂粘弹-粘损伤本构模型

为了描述固体推进剂在不同应变率和围压环境下的非线性力学特性，首先通过假设推进剂非线性力学特性由损伤导致，基于不可逆热力学框架，推导出粘弹-粘损伤本构模型。在构建粘损伤模型时，以线性粘弹性应变能密度为损伤驱动力，并且引入了损伤历史、应变率和围压效应对于损伤增长的影响。然后利用文献中HTPB推进剂的围压实验数据对一维形式下的本构模型进行了参数获取、验证和预测误差分析。在获取损伤萌发参数S0时，基于时间-压强等效原理，构建了损伤萌发参数S0主曲线。最后采用NEPE推进剂单轴拉伸实验验证了本构模型对于当前固体推进剂大变形非线性力学性能的适用性。结果表明，损伤萌发参数S0随着围压和应变率的增加而增加。在应变率和围压的双重作用下，在相对压强5.516MPa，0.24s-1条件下的S0是相对压强0MPa，6×10-4s-1条件下数值的10.7倍。另外，模型对于HTPB推进剂抗拉强度的最大预测误差为6.15%，模型预测结果与两种实验数据重合较好，表明建立的粘弹-粘损伤本构模型可以很好地预测HTPB推进剂在不同应变率和不同围压环境下的力学响应和当前NEPE推进剂的大变形非线性力学行为，可为点火增压载荷下固体推进剂药柱结构完整性数值分析提供理论基础。