EPDM包覆层热解状态下的力学行为及本构模型

为了准确描述固体火箭发动机三元乙丙(EPDM)包覆层在不同热解状态下的力学特性,通过热重分析的方法得到EPDM的热解温度范围,并采用准静态单轴拉伸实验的方法获得EPDM在不同热解温度下的应力-应变曲线。以Ogden(n=2)模型为基础,分别对不同热解状态下的力学曲线进行建模,并通过对实验数据的拟合获得模型参数。结果表明,在初始热解状态下,材料只发生了少量的失水和气体的逃逸,其力学行为可通过粘超弹本构模型进行描述;随着材料热解程度增加,基体材料和填充纤维均遭到破坏,其力学行为表现为超弹脆性材料特性。     

HTPB复合固体推进剂粘弹性应变能及非线性本构模型

为了准确表征HTPB复合固体推进剂在有限变形条件下的力学性能,针对推进剂粘弹性应变能及本构模型进行研究。提出了推进剂粘弹性应变能函数和非线性本构方程的一般形式,并通过一元非线性回归方法拟合不同应变率下的拉伸试验数据,得到了材料参数关于应变率的函数,并由此建立了推进剂单轴拉伸变形下的应变能函数和本构方程,预测了不同应变率下的应力曲线,与试验结果和已有模型的预测结果进行了对比。结果表明,材料参数与应变率之间呈现幂函数关系;推进剂应变能密度随变形量的增大呈非线性单调增长,同一变形条件下,应变率越高,推进剂的应变能密度越大;本构方程可准确描述推进剂拉伸变形的应力应变关系,且尤其适用于表征低应变率下,材料在有限变形内的粘弹特性。     

EPDM绝热层的粘超弹本构模型

为了准确描述固体火箭发动机内三元乙丙(EPDM)绝热层在有限变形下的力学特性,主要通过多步松弛和单轴拉伸2种实验方法,获得平衡响应曲线和拉伸曲线。分别采用Ogden模型、Mooney-Rivlin模型,对平衡响应曲线进行拟合,并以此为基础,引入应变率参数μ来描述EPDM在单轴拉伸实验中的率相关特性,从而建立了EPDM在准静态单轴拉伸条件下的本构模型。研究结果表明,本文所建立的模型能很好地描述EPDM在准静态条件下的力学特性,通过该模型,也可在有限应变率范围内预测EPDM的应力响应特性。     

双脉冲发动机EDPM软隔层含率相关能量限制器的黏超弹本构模型

在软隔层式双脉冲发动机中的设计过程中,软隔层的变形及破坏是一个很重要的问题。传统的表征橡胶类材料的超弹本构模型不能准确描述软隔层在工作过程中的变形及失效行为。为此,在Volokh等提出的含能量限制器的黏超弹本构模型的基础上,引入含率相关项的能量限制器,使其能够描述EPDM软隔层在单轴拉伸时的断裂伸长率的率相关性。为获取模型参数,进行了单轴拉伸试验,并采用最小二乘法拟合出模型参数。最后,利用所建的本构模型与试验结果进行了比较,比较结果表明模型能较好地预测EPDM在单轴拉伸时的失效行为。     

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

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

Kevlar层合板超高速撞击数值建模及参数识别

针对Kevlar层合板在超高速撞击条件下的力学特性的复杂性及参数不准确对数值仿真的不利影响,提出一种使用正交各向异性本构模型,Chang-Chang复合材料失效模型和有限元方法对Kevlar层合板的超高速撞击力学特性进行建模并运用连续响应面技术对材料失效模型参数进行识别的方法.然后使用该方法建立了与现有文献中的试验工况相对应的数值模型,并把Kevlar层合板材料失效模型的4个参数作为优化变量,以仿真结果与试验结果的二乘残差最小作为优化目标,使用连续响应面技术建立参数优化模型,由此识别这些参数.结果表明,该建模方法能正确描述Kevlar层合板的超高速撞击力学特性,采用识别后的参数进行仿真计算可以显著提高数值模型的仿真精度.     

固体推进剂动态力学行为研究进展

从动态力学表征、损伤机理、动态断裂和动态本构4个方面,对固体推进剂动态力学行为进行了综述。分析表明,动态力学实验重点关注固体推进剂的初始弹性模量、屈服应力、强度等参数,填充颗粒破裂是动态损伤的主要形态,动态断裂主要集中动态起裂研究,并对3类动态本构模型性进行了评述。在此基础上,梳理了下一步研究的重点,认为深化实验研究、力学数值模拟和含损伤热粘-超弹性动态本构模型将是下一步研究的重点。     

HTPB复合固体推进剂老化损伤的CT研究

针对HTPB复合固体推进剂老化存在的材料性能劣化问题，采用了计算机层析识别技术研究HTPB复合固体推进剂材料老化损伤．通过HTPB复合固体推进剂损伤的数学模型，定量描述了材料的老化问题，并将其代入推进剂粘弹本构方程．结果表明．老化损伤的粘弹本构方程能够一定程度上描述损伤老化引起的材料劣化问题。     

三维五向编织复合材料宏细观力学性能分析

基于三维五向编织复合材料的细观结构胞元,采用细观有限元方法建立了材料宏观等效弹性性能的力学分析模型。模型数值结果与拉伸试验数据基本吻合,说明模型可用于该材料的宏观力学性能预测。经分析其细观应力分布及变形,表明基于胞元分析的力学模型能有效预测含周期性结构材料的细观力学响应。由于该模型能给出较为合理的细观应力应变场,为进一步开展三维五向编织复合材料的失效分析奠定了基础。
     

柔性接头弹性件压缩剪切联合加载试验研究及数值分析

通过单轴拉伸与平面拉伸试验,得到了柔性接头弹性件橡胶材料“三阶五项式”超弹本构模型参数,通过有限计算证明了该模型的优越性.为了获得柔性接头弹性件材料剪切模量随容压的变化规律,设计了双轴加载试验装置和征柔性接头弹性件的三重片试件,进行了压缩剪切联合加载试验及有限元数值模拟,发现橡胶材料的剪切模量随压强的高略有增加.数值模...     

基于CRKSPH的炸药巴西压裂仿真研究

使用守恒型再生核粒子(CRKSPH)方法,结合弹脆性损伤本构模型,对PBX9501炸药标准巴西试验断裂过程进行了仿真,并对炸药内裂纹的产生和分布进行了预测。计算中,使用CRKSPH方法,可有效地提高计算精度和避免人工粘性误差;本构模型基于PBX9501炸药材料细观结构特征和破坏模式,并使用损伤参数D(单位面积上微缺陷的密度)来表征材料的破坏程度;计算结果中,试件断裂形式和中心点应力结果与文献中试验结果符合较好,而且加载速率对试件的断裂形式和断裂时间均有一定的影响。研究表明,CRKSPH方法在计算PBX9501炸药断裂损伤及裂纹扩展问题上有一定的工程实用价值。     

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

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

考虑应变率和温度响应的少烟NEPE推进剂粘弹性本构模型

为研究少烟NEPE推进剂力学性能的应变率及温度相关性,使用万能材料试验机分别在不同应变率(4.17×10-4～4.17×10-1s-1)和温度(-40～50℃)下测试了推进剂的力学性能,建立了考虑应变率和温度响应的本构模型。结果表明,少烟NEPE推进剂具有应变率硬化特性,其拉伸应力与对数应变率呈线性关系;降低温度使推进剂的定伸应力和模量增大,升温则相反。结合少烟NEPE推进剂的线性对数应变率效应和温度变化特性,建立了粘弹性本构模型,该模型用多蠕变模式与非线性弹簧的组合来反映力学性能的率相关性,用率相关模型与温度函数的乘积形式来描述力学性能的温度相关性。模型预测与实验曲线对比表明,所建模型在实验温度、应变率、0.1～1.0应变范围内预测的准确性较好,其百分误差小于24%。     

固体推进剂杆状药柱的动态响应

本文用三参数分数阶导数模型描述了了固体推进剂的粘弹特性，用Laplace变换和反奕换求解一维杆件的频响数物脉冲应函数，分析了极点部分和截部分在响应中所占比例，计算了在简谐激励了杆件的瞬态和稳态响应，并就加速度稳态响应部分与试验做了对比。     

Experimental research on rate dependent constitutive relation of double-basepropellant under impact load

为研究双基推进剂在准静态和高应变率下的压缩力学行为,利用万能材料试验机和分离式Hopkinson压杆(SHPB),对双基推进剂材料进行了单轴压缩实验,得到材料在10-4～103 s-1应变率下的应力-应变曲线.实验结果表明,双基推进剂是应变率敏感材料,屈服应力和初始弹性模量的对教与应变率的对数近似呈线性关系,且表现出韧脆转化现象.利用朱-王-唐非线性粘弹性本构关系,采用最小二乘法拟合了本构材料参数.研究表明,朱-王-唐本构模型能较好地描述双基推进剂在不同应变率条件下的力学行为.     

双基推进剂高应变率型本构模型的实验研究

为研究双基推进剂在准静态和高应变率下的压缩力学行为,利用万能材料试验机和分离式Hopkinson压杆(SHPB),对双基推进剂材料进行了单轴压缩实验,得到材料在10-4～103s-1应变率下的应力-应变曲线。实验结果表明,双基推进剂是应变率敏感材料,屈服应力和初始弹性模量的对数与应变率的对数近似呈线性关系,且表现出韧脆转化现象。利用朱-王-唐非线性粘弹性本构关系,采用最小二乘法拟合了本构材料参数。研究表明,朱-王-唐本构模型能较好地描述双基推进剂在不同应变率条件下的力学行为。     

基于GRNN的三元乙丙橡胶薄膜 粘接界面力学性能参数反演

基于广义回归神经网络(GRNN),提出了一种获取粘接界面的内聚强度和断裂韧性的数值方法。在ABAQUS中建立了双悬臂夹层梁(DCSB)的三维有限元模型(FEM),引入双线性内聚力模型描述界面单元的本构关系,采用不同的内聚参数计算加载点的载荷-位移关系。将有限元计算结果整理后作为广义回归神经网络的训练样本,将实验结果作为输入,预测了粘接界面的力学性能参数。结果表明,采用GRNN和Hooke-Jeeves算法获取粘接界面的力学性能参数分别需要进行25轮次和62轮次有限元分析,与GRNN和修正梁模型对应的加载点载荷-位移曲线和实验曲线之间的最大相对误差分别为1.12%和6.67%,基于GRNN的反演方法能够快速、准确地获取粘三元乙丙橡胶薄膜接界面的力学性能参数。     

固体火箭发动机绝热层脱粘的脉冲热像检测分析

为了给红外热像法在固体火箭发动机绝热层检测中的应用提供科学依据,用数值模拟法分析了绝热层脱粘脉冲热像检测的基本规律。得到了决策参数与结构参数之间的多变量关系,并对部分检测规律进行了实验验证。最大温差和最大对比度与脱粘尺寸之间的关系,可用分段线性函数近似描述,随着脱粘尺寸的增大,最大温差和最大对比度都增加;最大温差和最大对比度与绝热层厚度的关系是非线性函数关系,随着绝热层厚度增加,最大温差和最大对比度迅速下降;钢壳厚度的增加,对从绝热层一侧的单面法检测有利。所得数据和结论为固体火箭发动机绝热层脱粘的脉冲热像检测提供了指导。     

用三维非线性粘弹模型分析药柱结构

一个三维非线性粘弹模型结构分析有限元程序业已完成。该模型可运用于分析准静态条件下承受有限应变的各向同性、近乎不可压缩的粘弹材料,包括材料非线性影响.以热流变的简单性假设来建立与温度有关的模型。该模型已用于分析固化后冷却过程中的固体火箭发动机。当变形有限时,发现非线性模型预估的药柱应力水平高于线性粘弹模型预估的药柱应力水平.     

含能材料变泊松比粘弹性本构及其数值实现方法

为研究含能材料变泊松比粘弹力学行为,基于热流变简单材料假设,建立了一种采用时间相关泊松比的三维热粘弹性本构方程;利用积分算法对该本构进行数值离散,给出了相应的应力更新方法;基于Abaqus软件UMAT接口技术编制了相应材料子程序,完成了算例分析。结果表明,数值解与对应解析解吻合良好,并通过了合理性算例验证,分析方法及算法程序有效,可为后续三维装药精细结构完整性分析提供支持;粘弹材料泊松比具有时间相关性,其较小变化对计算结果影响显著,尤其是在短时加载或状态转变过程。