[发明专利]材料超弹性本构在LS-DYNA中的嵌入方法

权利要求书

1.一种材料超弹性本构在LS-DYNA中的嵌入方法，其特征在于，包括以下步骤：

步骤1、从主程序中提取当前时间步长内的储存于历史变量中的变形梯度张量F；

步骤2、根据材料的本构模型形式计算材料变形后与变形前的体积比J与右柯西格林应变张量C，确定应变张量不变量表示的弹性应变能函数W；

步骤3、若材料没有应变率效应，则计算第二皮拉霍夫应力张量S，若材料具有应变率效应，则根据相邻时间步长的变形梯度张量与时间步长dt确定右柯西-格林应变率张量不变量确定应变率张量不变量表示的黏性应变能函数W_e，然后计算第二皮拉霍夫应力张量S；

步骤4、确定真实柯西应力张量σ，将柯西应力张量σ传入主程序输出，进入下一时间步；

重复上述步骤1～步骤4，直至计算终止，最后输出终止时刻的结果。

2.根据权利要求1所述的材料超弹性本构在LS-DYNA中的嵌入方法，其特征在于，步骤1变形张量F矩阵表达形式为：

其中，F_ij为变形梯度张量的分量，i＝1,2,3，j＝1,2,3。

3.根据权利要求1所述的材料超弹性本构在LS-DYNA中的嵌入方法，其特征在于，材料没有应变率效应时，步骤3计算第二皮拉霍夫应力张量S为：

4.根据权利要求1所述的材料超弹性本构在LS-DYNA中的嵌入方法，其特征在于，材料具有应变率效应时，步骤3计算第二皮拉霍夫应力张量S为：

5.根据权利要求1所述的材料超弹性本构在LS-DYNA中的嵌入方法，其特征在于，步骤4确定真实柯西应力张量σ为：

σ＝J^-1·F·S·F^T。