基于LS-dyna模擬鈑金件撞擊壁面

教程目標：

熟悉接觸定義

了解怎樣設置時間步長

問題描述：

以速度為10m/s撞擊一個100X100X200的盒子模型，通過撞擊吸收的能量可以觀察在dyna中如何定義接觸算法，將盒子做成四分之一模型，邊界條件和尺寸如下圖所示：

材料屬性：

密度：7850kg/m3

楊氏模量：210GPa

泊松比：0.3

屈服極限：230Mpa

切線模量：500Mpa

材料設置：

設置鈑金的材料，點擊Model > Keywrd. MAT > 024-PIECEWISE_LINEAR_PLASTICTY，

設置運動墻的材料，將運動墻看做剛體，點擊MAT > 020-RIGID，CMO中選擇1，表示約束都在全局坐標系下進行的；CON1選擇5，表示約束Y和Z的平移方向；CON2選擇7表示約束所有旋轉方向。

設置單元屬性：

因為本例是殼單元，因此截面應是shell，點擊SECTION>SHELL，Title中填入crashbox，SHRF=0.8333，NIP=5，T1-T4=1.5表示殼單元的厚度是1.5mm，點擊accept，然后點擊NewID，改變Title為plate，改變厚度為2最后點擊accept，done。

賦予零件材料和截面：

點擊PART>PART，點擊對話框右上角，選擇零件并賦予對應的截面和材料。

設置運動參數：

點擊Model > Keywrd>DEFINE > CURVE，命名為velocity，在A1和O1分別輸入0，1，insert；1，1，insert，點擊Accept，done。

點擊BOUNDARY > PRESCRIED_MOTION_RIGID，在DOF中選擇1表示可以在X方向運動；VAD選擇0表示速度，如果選擇1表示加速度，如果選擇2表示位移。在LCID選擇剛剛定義曲線，在SF中填入因子系數，表示與曲線相乘，正負號表示方向。

設置接觸：

點擊Model > Keywrd>CONTACT > AUTOMATIC_ SINGLE_ SURFACE ，這是推薦設置，這里主要作用自身接觸。設置SSTYP=5，表示所有零件都被包括。填寫FS和FD=0.1，表示靜、動摩擦系數都是0.1。設置VDC=20，表示粘度系數；soft=1，表示接觸比較圓滑；IGNORE=2，表示模型中允許初始穿透。

在盒子尾部需要連接一個剛性面，用焊接的方式連接。點擊CONTACT > TIED_SHELL_EDGE_TO_SURFACE_OFFSET，先在SSTYP和MSTYP都選擇3，然后在SSID和MSID分別選擇1和2表分別指兩個零件。

設置終止時間：

點擊Model > Keywrd>control > TERMINATION，設置ENDTIM=0.015，隨后點擊Accept和Done。

設置輸出文件：

點擊Model > Keywrd>DATABASE > BINARY_D3PLOT，設置DT=2e-4，隨后點擊Accept和Done。

點擊Model > Keywrd>DATABASE > > ASCII_option，激活GLSTAT，設置DT=1e-6，隨后點擊Accept和Done。

設置控制卡片：

點擊Model > Keywrd>control > ENERGY，所有參數都選擇2，表示所有能量都要計算。

點擊Model > Keywrd> control > ACCURACY，設置OSU=1，INN=4。

點擊Model > Keywrd> control > SHELL，設置ESORT=1。

后處理：

打開d3plot文件，點擊Post>ASCII，選擇glstat*，并點擊load，選擇8-Sliding Energy，plot，即可查看滑移能，1-Kinetic Energy 可以查看動能。

關注公眾號：CAE備忘錄，回復 crash 可獲得結果k文件。