眾所周知，針對硬脆材料的模擬通常采用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軟件中定義并輸出所給材料的歷史變量。