基于全隐式欧拉算法的粘塑性损伤本构模型数值实现方法

本发明提供了一种基于全隐式欧拉算法的粘塑性损伤本构模型数值实现方法，涉及火箭发动机的技术领域，该方法包括：基于全隐式欧拉算法，对预先构建的推力室用铜合金粘塑性损伤本构模型进行数值化计算，确定Jacobian矩阵的解析式和一致性切线算子；基于所述Jacobian矩阵的解析式和所述一致性切线算子，确定粘塑性损伤本构方程的材料子程序；利用预先确定的粘塑性损伤本构模型输入参数，调用所述粘塑性损伤本构方程的材料子程序，实现推力室结构分析，该方法解决了现有推力室粘塑性损伤本构模型的数值实现方法缺乏的问题，并且达到了提高计算稳定性和准确性的技术效果。