kaolvwendu new.k

      目前我在做一個模擬火炸藥撞擊感度測試的仿真，其中主要內容是把火炸藥薄片夾在上下擊柱間并置于套筒中，對上擊柱賦予3kg的質量，隨后讓它以不同的速度與炸藥薄片進行碰撞。大致如圖所示。

      建模后上下擊柱留出了一點間隙，如圖：

      目前的第一階段，我不考慮炸藥的化學反應，只考慮材料的壓縮生熱，因此我對中間的綠色藥片部分使用的是JOHOSON_COOK關鍵字，搭配EOS_GRUNEISEN狀態方程；上下擊柱和套筒都使用一種鋼材，也用JOHOSON_COOK關鍵字，搭配EOS_GRUNEISEN狀態方程。

       但是在實際計算中，計算過程會出現復雜聲速的報錯，

 *** Error 40510 (SOL+510)

     complex sound speed in solid element # 97194 cycle 14118

 *** Error 40510 (SOL+510)

     complex sound speed in solid element # 97214 cycle 14118

   14118 t 6.6122E+01 dt 3.28E-06 write d3plot file            08/03/23 15:03:52

   14118 t 6.6122E+01 dt 3.28E-06 write d3plot file            08/03/23 15:03:52

結果文件導入prepost中后是這樣的情況（其他部分都比較正常，我只貼出藥片和下擊柱部分）：

上面是藥片，這是下擊柱：

我也嘗試過更換材料模型，將中間的綠色藥片也換成跟摩擦板一樣的鋼材質，結果并沒有什么變化。

而當我把中間的藥片換成MAT_PIECEWISE_LINEAR_PLASTICITY材料模型后，則會出現負體積，但是看起來好了很多：

 *** Error 40509 (SOL+509)

     negative volume in solid element # 93174 cycle 13982

   13982 t 6.4856E+01 dt-1.01E-04 write d3plot file            08/03/23 15:39:30

   13982 t 6.4856E+01 dt-1.01E-04 write d3plot file            08/03/23 15:39:30

我想解決這個問題，但又不知從何著手，下面是我用的材料模型具體參數，我用的單位制是cm-g-us

上面這個是鋼質摩擦板的參數。

這個是hmx炸藥的參數。

這個是摩擦板的狀態方程。

這個是藥片的狀態方程。

我現在完全不知道問題出在哪里，細化過網格后并沒有改善，之前使用其他材料模型我記得有成過，但是沒有保存（畢竟藥片用其他材料意義不大），希望可以得到大佬幫助，指出問題所在。k文件附上，希望能得到回應，想做生熱方面內容的朋友也可以下載查看。