HTPB推进剂定应变老化损伤本构特性研究

为研究复合固体推进剂在定应变老化条件下的力学特性及本构行为,开展了HTPB推进剂方坯定应变加速老化试验。考虑到推进剂内部缺陷分布的随机性,建立了定应变老化条件下推进剂的统计损伤本构模型。在构建损伤变量方程时,通过引入初始损伤系数,将推进剂的老化过程等效看作是一种损伤过程。基于定应变加速老化试样的单轴拉伸试验数据,对模型进行了验证,结果显示理论模型与试验结果吻合较好,能准确地描述推进剂在定应变老化过程中的力学行为。通过对模型参数进一步分析发现:初始损伤系数方程能区分化学老化与定应变对推进剂的损伤作用。化学老化对推进剂的损伤效应随老化时间呈指数规律增大;同一老化时间段内,定应变对推进剂的损伤作用呈指数规律增大。化学老化不但影响本构模型曲线在损伤段强度的大小,还影响其形状特性,在55℃条件下,老化时间小于284d时,试样拉伸曲线具有明显的屈服区;老化时间大于284d时,没有明显的屈服区。定应变仅降低损伤段强度的大小,不改变其形状特性;在定应变老化过程中,损伤应变阈值随老化时间及定应变水平基本呈现线性关系下降。     

HTPB固化胶片的超弹性本构模型

为描述HTPB固化胶片的力学特性,采用单轴拉伸试验方法,研究了载荷作用下的应力响应特征,建立了超弹性本构模型,并拟合试验数据得到超弹性本构参数;对比分析了Mooney-Rivlin和不同阶次的Ogden型应变能函数。结果表明,Mooney-Rivlin和二阶以上的Ogden型拟合的本构参数都能很好地表征HTPB固化胶片拉伸载荷下的力学行为,且Ogden型超弹性本构可用于预测HTPB固化胶片轻微压缩载荷下的力学特性。     

金属丝网橡胶的本构关系

通过对圆环形金属丝网橡胶成型制造过程的分析,在圆环嵌套模型基础上,结合干摩擦非线性理论和小曲梁模型建立了金属丝网橡胶材料的非线性本构关系,对不同相对密度的金属丝网橡胶试件进行了静态压缩实验,通过实验数据确定了本构关系方程的各项系数并实验结果进行比较.研究结果表明:建立的本构关系模型可以很好地描述金属丝网橡胶的非线性力学特性,实验结果能够较好地符合理论结果,并且该模型反映了相对密度、网格大小和网格宽度等参数的影响,从而实现了对材料结构刚度的预估.      

针刺陶瓷基复合材料损伤本构模型及构件应力分析

基于连续介质损伤力学方法提出了适用于针刺陶瓷基复合材料的各向异性损伤本构模型.该模型考虑了拉伸与剪切损伤的影响.通过试验得到了材料的拉伸和剪切的应力-应变曲线,采用曲线拟合获得了本构模型的参数.应用用户子程序技术将本构模型写入商用有限元软件,计算所得拉伸和剪切的应力-应变曲线与试验曲线比较吻合,最大误差分别为5.62%和1.47%.使用该损伤本构模型及弹性模型计算了航空发动机尾喷管调节片在气动和温度载荷下的力学响应,两者计算所得应力分布相似,但损伤本构模型计算值明显小于弹性模型计算值.      

复合推进剂线粘弹本构方程的细观力学分析(Ⅰ)球形颗粒增强效应分析

从基体和微粒本构关系出发, 应用Ｅｓｈｅｌｂｙ 等效夹杂理论及宏观平均理论的分析方法, 确定复合体的本构关系, 从细观层次揭示弹性微粒对粘弹基体的增强效应。     

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

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

HTPB推进剂拉伸力学行为的应变速率相关超弹本构模型

为研究HTPB推进剂拉伸力学行为的应变速率相关性,采用万能试验机和液压试验机,在室温下开展了不同应变速率(1.2×10~(-4)~80s~(-1))的拉伸实验。结果表明:给定应变对应的应力随应变率对数双线性增加,1s~(-1)为双线性关系的转折点;随应变的减小,HTPB推进剂的应变率敏感性线性增强。在Mohotti建立的模型基础上,结合拉伸力学行为的双线性率相关特性及应变率敏感性的应变依赖性,提出改进的应变速率相关超弹本构模型。该模型以超弹性元件作为基础描述参考应变率下的力学行为,率相关元件乘入超弹性元件描述率相关特性。模型预测与实验曲线对比表明:所提出的率相关超弹本构模型能够描述HTPB推进剂在1.2×10~(-4)~80s~(-1)应变率、30%应变范围内的拉伸力学行为。     

TiC 颗粒增强钛基复合材料的动态损伤本构

分析了TiC颗粒增强钛基复合材料的微损伤演化规律,建立了含损伤演化的动态本构模型。TiC颗粒增强钛基复合材料在拉伸载荷作用下,微裂纹以翼型裂纹形式扩展,基于平面翼型裂纹扩展模型,建立了二维动态损伤本构关系,并退化到一维拉伸状态,假设微裂纹成核规律满足Weibull分布,得到了一维拉伸应力状态下能够反映TiC颗粒增强钛基复合材料的损伤演化规律的宏细观相结合的动态本构关系。模型计算结果与试验结果吻合较好。     

耦合蠕变损伤的Chaboche粘塑性本构方程的应用

为了改善 Chaboche粘塑性统一本构方程描述第三阶段蠕变的能力 ,按照 Lemaitre有效应力的概念 ,采用 Kachanov损伤演化方程 ,推导了耦合各向同性损伤因子的 Chaboche粘塑性流动方程和硬化方程。当耦合损伤时 ,内变量演化方程形式不变 ,并不显含损伤因子 D。将之用于镍基高温合金 Udimet72 0 Li的蠕变描述 ,得到了该合金的蠕变损伤参数。通过对双轴试样的有限元计算 ,发现了 Kachanov损伤演化方程对单轴和多轴蠕变的预测并不协调。在损伤演化方程中引入多轴损伤因子得到了与实验数据基本一致的蠕变寿命结果。     

基于修正N-S方程本构关系的气体动理论耦合方法

随着稀薄程度的增加,Navier-Stokes方程的线性本构关系难以正确描述稀薄气体输运特性,高阶非线性本构关系往往数学形式极为复杂,对数值求解造成稳定性差等问题。为了发展适宜于近空间飞行器气动特性分析的高超声速稀薄流动模拟方法,本文利用求解Boltzmann模型方程的气体动理论统一算法(Gas Kineitc Unified Algorithm, GKUA)对应力张量、热流等宏观量数值积分求解的优势,提出了一种基于数值修正N-S方程本构关系的气体动理论耦合方法。通过将GKUA获得的应力张量及热流用于修正N-S方程的本构关系,实现了存在局部稀薄效应的流动模拟,并且通过可压缩平板边界层、圆柱绕流问题的数值模拟,验证了方法的有效性。