【LSDYNA操作小技巧七】LSDYNA中具体材料历史变量的定义与输出问题

       众所周知，针对硬脆材料的模拟通常采用JHC或JH-2本构模型模拟，玻璃类材质、混凝土等是这类硬脆材料的典型代表。相同的是，无论JHC还是JH2本构模型都内置有材料的累积损伤模型用以准确描述硬脆材料在收到破坏时的内部单元失效变化，反映在宏观上即产生裂纹、断裂、损伤等。研究上述诸如裂纹等材料大变形问题对于深入认知本构模型的作用机理就显得尤为必要。

       LSDYNA作为专门针对材料非线性、大变形问题而开发的大型显示动力学仿真软件，其拥有相当齐全的材料库，针对其中材料的本构关系都有非常明确的软件设置与输出步骤，方便用户对本构的快速设置甚至二次开发。本帖以典型的材料大变形问题磨粒切削加工引起的加工损伤为案例讲解损伤历史变量的定义与输出步骤。

       在我们没有更好的文献参考时，官方仿真软件的帮助文档是一个最佳的选择。以下提供在LSDYNA中具体查阅历史变量定义的链接如下，

History Variables for Certain Material Models — Welcome to the LS-DYNA support site.

       基于此，必须首先定义目标历史变量如图1所示。具体通过，在关键字 *DATABASE_EXTENT_BINARY中定义NEIPH 或 NEIPS 实现。这里NEIPH 或 NEIPS 只能填入具体的数字加以定义材料具体哪些的历史变量。

图1

       具体以此为例：为输出磨粒切削加工引起的工件损伤，首先找到对应的工件材料编号110（工件为K9玻璃，采用JH-2本构），可知损伤因定义为2如图1所示，故在NEIPH输入数字2即可，若为了查看其他材料历史变量，对应输出1-4即可实现。

图2

       完成定义后，完成模型建模输出为k文件，在LSPP中进行损伤云图的输出。图3为定义的损伤在后处理中的查看方式，图4为输出的损伤云图，红色粒子表示损伤。

图3

图4

       以上，基于此案例，可在LSDYNA软件中定义并输出所给材料的历史变量。