固体火箭发动机药柱大变形数值分析

综合考虑固体发动机药柱的近似不可压缩性和几何非线性,基于不可压缩材料的粘弹性本构关系,应用完全拉格朗日(T.L)法的虚功方程,推导了二维粘弹性几何非线性有限元的增量列式,编写了有限元程序并与算例进行了对比,并对星型药柱在受压力载荷作用下的结构进行了分析。结果表明,该方法和程序可用于所有泊松比粘弹性问题的计算,尤其适用于推进剂药柱的应力分析。     

固体火箭发动机药柱三维结构非线性分析

非线性分析是药柱结构分析中的难点，针对其对结构分析的影响，基于不可压缩材料的粘弹性本构关系，应用完全拉格朗日（T．L）法的虚功方程，综合考虑药柱的近似不可压缩性和几何非线性，推导了三维粘弹性几何非线性有限元增量方程，编写了有限元程序对星型药柱在受压力载荷以及固化降温载荷作用下的结构进行了分析，并与线性计算的结果进行了对比，结果表明，在大载荷作用下，非线性对药柱结构分析的影响比较显著，计算时考虑非线性的影响是必要的。     

常应变率下某固体推进剂非线性粘弹性本构关系研究

以工程可应用性和简单性为目的,开展某固体推进剂非线性粘弹性本构关系研究。首先从积分型线性粘弹性本构关系出发,通过松弛实验获得了某固体推进剂材料的松弛模量,建立了线性粘弹性本构关系。再引入随应变率变化的非线性松弛模量(含非线性弹簧和非线性粘壶)概念,建立了考虑应变率效应的积分型非线性粘弹性本构关系。最后,借助于常应变率拉伸实验曲线,提出参数识别策略,得到了三参数的非线性松弛模量。多种常应变率拉伸实验结果表明,所建立的非线性粘弹性本构关系不但与线性粘弹性本构关系具有很好的衔接性,而且因只含有3个物理意义明确的材料参数,极易于工程应用。     

含损伤非线性粘弹性本构模型及数值仿真应用

提出了一种含累积损伤的非线性粘弹性本构方程,通过数值展开,获得该应力更新方程和切线模量矩阵,利用FORTRAN语言编制了自定义材料本构模型子程序,将子程序融合到ABAQUS软件中进行有限元数值仿真。通过验证试验证明了所建立的非线性粘弹性本构方程及数值方法的准确性和实用性,为装药结构完整性数值仿真提供了较为精确的本构模型及材料子程序。     

考虑粘弹性泊松比的固体推进剂蠕变型本构模型

在粘弹性材料中,泊松比是时间和温度相关的函数。提出了一种新的考虑时间和温度相关泊松比的推进剂蠕变型本构模型。为了建立合适的本构模型,采用拉普拉斯变换方法对经典蠕变型本构模型进行了改进,将时间无关的弹性泊松比变换成时间相关的粘弹性泊松比。基于时温等效原理,将泊松比的时间和温度相关特性利用WLF方程进行建立联系。基于蠕变型本构模型的增量形式将该模型开发到有限元程序MSC.MARC中。针对星孔发动机点火增压工况进行了结构完整性分析。仿真结果表明,粘弹性泊松比对结构的应力应变影响介于初始和平衡泊松比对结构分析的影响之间。并且该蠕变型本构模型与松弛型本构模型的分析结果一致。     

基于分数阶三参数模型的HTPB推进剂粘弹性力学响应分析

为拓展分数阶粘弹性模型在端羟基聚丁二烯(HTPB)推进剂中的应用,研究分数阶粘弹性模型预示推进剂频域下的动态模量变化。推导了三维下的分数阶三参数模型,并结合Grünwald-Letnikov定义编写UMAT子程序。为标定本构参数,推导了带有预加载的蠕变、松弛响应的解析解,并采用遗传算法分别标定了0.4 MPa和0.6 MPa应力下蠕变实验以及10%和30%应变下松弛实验的本构参数。计算得到的有限元数值与解析解、解析解与实验值的时程相对误差均小于5%。采用分数阶三参数模型动态模量对-5、25、60℃下的动态力学实验(DMA)结果进行拟合得到本构参数,动态模量与实验结果误差小于2%,同时随着温度升高,本构参数均减小,体现出HTPB推进剂升温软化现象。结果表明,分数阶三参数模型能精确地描述不同温度下多种频率的HTPB推进剂模量变化,结合二次开发本构,可为后续HTPB推进剂复杂边界下的力学响应的研究提供材料本构模型。     

复合固体推进剂非线性粘弹本构方程的微观力学分析

在填充颗粒为弹性硬颗粒的条件下，本文成功地将Ｅｓｈｅｌｂｙ的等效夹杂理论及相应的平均化方法推广应用于研究复合固体推进剂的非线性粘弹性本构方程，结果表明，复合固体推进剂的非线性本构关系可以用基体的本构关系与一个增强系数的 积β来表示，β只与颗粒的形状，体积分数及泊松比有关。本文的结论具有重要的工程应用价值。     

一种粒子增强体复合材料的线粘弹性模型

定量分析了粒子增强体复合材料的粘弹性力学性能。首先,依据复合材料的细观力学理论与Laplace变换理论,在Tan提出的粘弹性模型基础上,发展了一种改进的粒子增强体复合材料的三维线粘弹性模型,克服了Tan的模型只能刻画界面脱湿的缺点;其次,采用非线性有限元方法,实现了粒子增强体复合材料在常应变率下粘弹性本构关系的数值模拟。结果表明,该模型与程序可用于定量化分析多种情况下粘弹性问题的应力-应变关系,如不同粒径、温度、常应变率及粒子体积分数等参数,并得出粒径大的固体颗粒与粒子体积分数高的复合材料都容易发生界面脱湿。     

HTPBⅠ型裂纹粘弹性有限元分析

在简化模型的基础上，推导了复合固体推进剂中Ⅰ型拉伸裂纹的裂尖区域应力、应变与位移的关系式，并进行了粘弹性有限元方法分析。结果表明，有限元分析能获得具一定精度的结果，并与理论分析基本吻合。     

适用于固体推进剂药柱应力分析的一种粘弹性有限元法

本文根据Ｈｅｒｒｍａｎｎ变分原理导出一种适用于不可压缩和近似不可压粘弹性的新的本构关系和相应的有限元公式，开发了用于平面问题分析的计算机程序ＶＦＡＰＩＮＰ。计算结果表明，本方法和程序可用于所有泊松比的粘弹性问题的计算，尤其适用于推进剂药柱的应力分析。