① ls-dyna演算法的詳細介紹
以前有這方面很詳細的資料，我簡單整理一下，用來回答你的問題
LS-DYNA中的接觸類型大體上可以分為五大類：
One-Way Contact （單向接觸）
Two-Way Contact（ 雙向接觸）
Single Contact（單面接觸）
Entity
Tied Contac（固－連接觸）
在以上接觸類型中，前四種接觸類型的接觸演算法均採用罰函數法。
固－連接觸有的採用的罰函數法，有的採用動約束法，少部分採用分布參數法。
Tied Contac
(1)Translational DOF only, No Failure, No Offset
這種類型接觸採用動態約束演算法。
如下兩個命令是常用的固連接觸
*Contact_Tied_Nodes_To_Surface(6)
*Contact_Tied_Surface_To_Surface(2)
(2)Translational DOF only, No Failure, With Offset
這種接觸採用罰函數演算法，
與上述接觸類型2、6對應的為
*Contact_Tied_Nodes_To_Surface_OFFSET(O6)
*Contact_Tied_Surface_To_Surface_OFFSET(O2)
(3)Translational DOF & Rotational DOF, With Failure, No Offset
採用動態約束演算法。
(4)Translational DOF & Rotational DOF, With Failure, With Offset
罰函數法。
(5)Translational DOF Only, With Failure, With Offset
動態約束演算法。
作爆炸分析採用的是僅滑動接觸演算法，當炸葯和金屬的單元節點對應時不會出現節點穿透，可是這樣就限制了單元的劃分。當炸葯和金屬的單元節點不對應會出現穿透，而且增加接觸剛度的值也控制不了，這個時候還有其它控制穿透的辦法嗎？
下面僅滑動接觸關鍵字：
*CONTACT_SLIDING_ONLY
$ SSID MSID SSTYP MSTYP ***OXID MBOXID SPR MPR
2 1 0 0 0 0
$ FS FD DC VC VDC PENCHK BT DT
0.0 0.0 0.0 0.0 0.0 0 0.01.0000E+20
$ SFS SFM SST MST SFST SFMT FSF VSF
1.0 1.0 0.0 0.0 1.0 1.0 1.0 1.0
*CONTROL_CONTACT
$ SLSFAC RWPNAL ISLCHK SHLTHK PENOPT THKCHG ORIEN ENMASS
1.00 2 0 1 0 3 0
$ USRSTR USRFRC N***CS INTERM XPENE SSTHK ECDT TIEDPRJ
0 0 10 0 4.0 0 0 0
$ SFRIC DFRIC EDC VFC TH TH_SF PEN_SF
0.0 0.0 0.0 0.0 0.0 0.0 0.0
$ IGNORE FRCENG SKIPRWG OUTSEG SPOTSTP SPOTDEL
1 0 0 0 0 0