實用有限單元法計算及其程序的自動生成

1、61530Y 實用有限單元法計算及其程序的自動生成第五講 FEPG系統(tǒng)的前后處理數(shù)據(jù)結(jié)構(gòu)及程序調(diào)試目標(biāo) l 掌握FEPG的前處理數(shù)據(jù)文件-PRE文件的基本功能和語法l 掌握FEPG的后處理數(shù)據(jù)文件-POS文件的基本功能和語法l 掌握FEPG的輸入數(shù)據(jù)結(jié)構(gòu)和離散數(shù)據(jù)文件的存儲格式提綱1、 PRE文件的基本內(nèi)容和功能l PRE文件的三段內(nèi)容。第一段給出坐標(biāo)系以及各個物理場的總體信息。第二段給出各個物理場的單元類型信息及單元的材料參數(shù)名。關(guān)鍵字 element。第三段信息給出各個物理場單元的材料參數(shù)值。關(guān)鍵字 matedata。l 段落之間用#分隔2、 PRE文件的自動生成和修改l PRE文件是前處

2、理的數(shù)據(jù)文件，由他給出有限元的數(shù)據(jù)結(jié)構(gòu)l 通過PRE文件，F(xiàn)EPG系統(tǒng)生成與前后處理子系統(tǒng)接口的接口文件l Pre文件用來給出有限元問題的所有輸入數(shù)據(jù)信息l Pre文件由系統(tǒng)根據(jù)gcn, gio, pde文件自動生成，但內(nèi)容常需要修改3、 POS文件的結(jié)構(gòu)和語法4、 ID0, ELEM0, DISP0, COOR0等文件的結(jié)構(gòu)及命名原則輸入數(shù)據(jù)節(jié)點坐標(biāo)數(shù)據(jù)節(jié)點規(guī)格數(shù)據(jù)節(jié)點指定位移和節(jié)點載荷數(shù)據(jù)節(jié)點初始位移數(shù)據(jù)節(jié)點初始速度數(shù)據(jù)節(jié)點初始加速度數(shù)據(jù)單元信息表格coor0文件disp1文件id0文件disp0文件disp2文件disp3文件第1種材料單元節(jié)點信息數(shù)據(jù)第1種材料單元材料性質(zhì)數(shù)據(jù)第n種材料

3、單元材料性質(zhì)數(shù)據(jù)第n種材料單元節(jié)點信息數(shù)據(jù)elem0文件5、 GID造型的一些簡單技巧例題六考慮孔周受力問題。長方區(qū)域的左下角點和右上角點坐標(biāo)分別是(0,0), (2,1). 中間挖去兩個半徑0.2的圓，圓心分別是(0.5,0.5), (1.5,0.5)。其他邊值條件和材料參數(shù)與例題二相同。計算孔周的受力情況。例題七非線性瞬態(tài)熱傳導(dǎo)問題¡ 見SAMPLES¡ Heat->2dxy->nonlinear-dynamic-heat¡ 要點：pos文件及動畫顯示練習(xí)五如圖所示， 一個內(nèi)徑0.5m, 外徑5m的彈性圓環(huán)， 內(nèi)圓邊界作用脈沖壓強p后，產(chǎn)生彈性波。圓環(huán)外為某種均勻介質(zhì)。忽略材料自重影響，求彈性波在介質(zhì)中的傳播過程。彈性球環(huán)的彈性模量為1.0E10N/m2, 泊松比 0.3, 密度 3000kg/m3, 脈沖壓強p = 100N/m2. 時間步長取0.01SEC. 介質(zhì)的彈性模量為1.0E