這是一種彈性常數突變的損傷模擬方法，我編寫子程序，模擬層合板沖擊損傷，用C3D8R，發現ALLAE過大，是因為沙漏？我把中間一部分（沖擊頭附近）改為C3D8，發現穩定時間增量步變小，算得較慢。