含损伤粘弹性本构及其在有限元分析中的实现

为研究固体发动机药柱的损伤粘弹性行为,针对一种耦合损伤的三维非线性粘弹性本构模型,采用积分算法对其进行了数值离散,导出了本构方程的增量有限元形式。针对该增量型本构方程对应的非线性方程组,建立了采用newton-raphson迭代法求解的应力更新方法,并基于abaqus软件和fortran语言开发了相应材料本构子程序,最后应用算例对所设计数值方法及程序进行了验证与分析。结果表明,本文有限元解与对应解析解吻合良好,保持了很高的精度,所设计方法及程序有效,可应用于实际型号发动机药柱损伤结构分析。