有限元仿真分析實驗一、實驗?zāi)康耐ㄟ^剛性球與薄板的碰撞仿真實驗，學(xué)習(xí)有限元方法的基本思想與建模仿真的實現(xiàn)過程，并以此實踐相關(guān)有限元軟件的使用方法。本實驗使用HyperMesh軟件進(jìn)行建模、網(wǎng)格劃分和建立約束及載荷條件，然后使用LS-DYNA軟件進(jìn)行求解計算和結(jié)果后處理，計算出鋼球與金屬板相撞時的運(yùn)動和受力情況，并對結(jié)果進(jìn)行可視化。二、實驗軟件HyperMesh、LS-DYNA三、實驗基本原理 本實驗?zāi)M剛性球撞擊薄板的運(yùn)動和受力情況。仿真分析主要可分為數(shù)據(jù)前處理、求解計算和結(jié)果后處理三個過程。前處理階段任務(wù)包括：建立分析結(jié)構(gòu)的幾何模型，劃分網(wǎng)格、建立計算模型，確定并施加邊界條件。四、實驗步驟1、按照點線面的順序創(chuàng)建球和板的幾何模型（1）建立球的模型：在坐標(biāo) (0,0,0)建立臨時節(jié)點，以臨時節(jié)點為圓心，畫半徑為5mm的球體。（2）建立板的模型：在tool-translate面板下node選擇臨時節(jié)點，選擇Y-axis,magnitude輸入5.5，然后點擊translate+，return;再在2Dplanes-square 面板上選擇Y-axis,B選擇上一步移下來的那個節(jié)點，surface only ,size=30。2、畫網(wǎng)格（1）畫球的網(wǎng)格：以球模型為當(dāng)前part，在2Datuomesh面板下，surfs選擇前面建好的球面，element size設(shè)為0.5mm，mesh type選擇quads，選擇elems to current comp,first order， interactive。（2）畫板的網(wǎng)格：做法和設(shè)置同上。3、對球和板賦材料和截面屬性（1）給球賦材料屬性：在materials面板內(nèi)選擇20號剛體，設(shè)置Rho為2.000e-08,E為200000,NU為0.30。（2）給球賦截面屬性：屬性選擇SectShll,thickness設(shè)置為0.1，QR設(shè)為0。（3）給板賦材料屬性：材料選擇MATL1，其他參數(shù)：Rho為2.000e-08,E為100000,Nu為0.30，選擇Do Not Export。（4）給板賦截面屬性：截面選擇SectShll，thickness設(shè)為0.2。其他參數(shù):SHRE為8.333-01，QR為0，T1為0.2。（5）給板設(shè)置沙漏控制：在Properties-Create面板下Card image選擇HourGlass，IHQ為4，QM為0.100。更新平板。4、加載邊界條件（1）將板上最外面的四行節(jié)點分別建成4個set。（2）建立一個load collector。（3）Analysisconstraints面板中，設(shè)置SIZE為1，nodes通過by sets選擇set_1、set_2、set_3、set_4，然后點擊creat即可，邊界條件加載完畢。5、建立載荷條件（給球一個3mm的位移）（1）建立一個plot: postxy plots-plots-creat plot,然后點擊return;(2)在post-xy plots-edit curves面板中輸入X0，0，0.0001,Y0,3。(3)給剛性球一個3mm的沿y正方向的位移：card image設(shè)為PrcrRgd,DOF為2,VAD為2,LCID為1,SF為1,option選擇Rigid。6、接觸處理（1）做兩個用于接觸的segment：在Analysis-set_segment面板中，Card image選擇setSegment，elems選擇球這個part。重復(fù)操作對平板創(chuàng)建segment。要同時保證球的setsegment的方向朝外，plane的setsegment方向朝上。（2）建立接觸：對稱接觸：在Analysis-interfaces面板中,type選擇SurfaceToSuface,Card image選擇SingleSurface。master的contactsurfs選擇球的setgment;slave的contactsurfs選擇平板的setgment。點擊edit后，設(shè)置FS為10，F(xiàn)D為10，SFS為100，SFM為100，Automatic一項勾選OneWay。非對稱接觸：在Analysis-interfaces面板中,type選擇SingleSurface,Card image選擇SingleSurface。master的contactsurfs選擇球的setgment;slave的contactsurfs選擇平板的setgment。點擊edit后，設(shè)置FS為10，F(xiàn)D為10，SFS為100，SFM為100，Automatic一項勾選smooth。7、定義控制卡片在Analysis-control cards面板中，(1)選擇Control_Enegy,將hgen設(shè)置為2；(2) 選擇Control_Termination,將ENDTIM設(shè)為0.0001s;(3) 選擇Control_Time_step,將DTINIT設(shè)為1*10-6s，將TSSFAC設(shè)置為0.6;(4) 選擇DATABASE_BINARY_D3PLOT,將DT設(shè)置為5*10-6;(5) 選擇DATABASE_OPTION，將MATSUM設(shè)置為1*10-6，將RCFORC設(shè)置為1*10-6。8、刪除臨時節(jié)點在geom中選擇temp nodes，點擊node，選擇all，然后點擊clear。9、節(jié)點重新排號。tool中選擇renumber，然后選擇all，最后點擊renumber。10、將文件導(dǎo)出成KEY文件。11、生成的KEY文件導(dǎo)入LS-DYNA中，并運(yùn)行。12、打開后處理程序Ls-Prepost，點擊file-open-ls-dyna binary plot，選擇計算得到的d3plot結(jié)果。13、點擊下方運(yùn)行按鈕，觀看碰撞動畫。選擇右邊的ASCII按鈕，選擇rcforc后，點擊左邊的Load按鈕，再選擇Sl或者M(jìn)a，然后選擇最下方的Resultant force,最后點擊Plot，觀看沖擊力隨時間的變化曲線。五、實驗結(jié)果1、對稱接觸情況對稱接觸定義下，沖擊過程中薄板的彈塑性變形過程如圖1至圖5所示。圖1圖2圖3圖4圖5沖擊過程中剛性球與薄板層接觸力時間歷程如圖6所示。圖62、非對稱接觸情況非對稱接觸定義下，沖擊過程中薄板的彈塑性變形過程如圖7至圖11所示。圖7圖8圖9圖10圖11沖擊過程中剛性球與薄板層接觸力時間歷程如圖12所示。圖12六、實驗總結(jié)根據(jù)對稱接觸定義與非對稱接觸定義的對比實驗計算過程和實驗結(jié)果可以看出，對稱接觸定義下剛性球撞擊薄板的沖擊現(xiàn)象更加明顯，接觸力更大，計算所花費的時間也更長，這與理論上雙接觸面與單接觸面的結(jié)論是相同的。