搜索一個說明有限元法數(shù)學原理的簡單實例
2017-03-23 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
在數(shù)學層面說明有限元法的數(shù)學原理并不是一件容易的事情,因為當涉及到二維或者三維維度時,一個簡單的問題都需要結(jié)合很多原理和公式來演算。為了讓大家對有限元法的數(shù)學原理有個簡單的認識和了解,今天給大家介紹一個通俗易懂的例子來說明有限元法的(FEM)的數(shù)學思想,希望對大家有所啟發(fā)。(在我共享的資料中,其中科羅拉多州立大學的有限元學習資料在開頭就有一個講得非常詳細的桁架結(jié)構(gòu)的有限元法的數(shù)學講解,大家可以仔細閱讀做為補充)。文末鏈接又給大家提供了一個相對簡單但又通俗全面的有限元法的說明資料,大家可以下載閱讀。
如下圖所示,一個圓截面桿件一端固定,另一端施加拉力P,計算桿件的變形以及受到的應力。
模型是一個簡單的桿件,可以簡化為1D單元,1D單元有不同類型,本例設(shè)定為桁架單元(truss element),只能承受軸向力。如下圖所示,為了方便計算,我們把桿件劃分為2個單元和3個節(jié)點(節(jié)點a,b,c)。
其中每個節(jié)點都受到相應的外力(由外部載荷P引起)和內(nèi)力l(由內(nèi)部應力引起)。當模型處于靜力學平衡時(static equilibrium),節(jié)點力(作用于節(jié)點的內(nèi)力和外力的合力)必須等于0,每個節(jié)點的受力平衡如下圖所示:
假設(shè)桿件變形過程中伸長量很小,對于單元1(element 1)的應變有:
其中
和
分別為節(jié)點a和b的位移(displacement),L為單元1(element 1)的長度。假設(shè)材料為彈性材料,楊氏模量為E(Young's Modulus),則由材料力學可得單元1的應力:
作用于節(jié)點a的軸向力等于應力乘以橫截面積(cross-sectional area),所以可以得到內(nèi)力、材料屬性以及位移的關(guān)系如下式:
由力學平衡可得:
即:
同樣的方法,運用力學平衡和材料力學公式可以得到節(jié)點b的平衡式如下:
節(jié)點c的平衡式如下:
將節(jié)點a,b,c的平衡式寫成矩陣形式,得:
這樣,我們就得到了桿件的平衡方程,由于節(jié)點a固定,所以位移定于0,又聯(lián)立Pb=Pc=P,這樣我們就可以得到節(jié)點b和節(jié)點c的位移以及Pa,算出位移之后,我們就可以返回去算出應力值。
其中EA/L就是我們常說的剛度,當然在很多時候,每個單元的剛度并不一樣,如本例,如果兩個單元的長度不一樣,則剛度也不一樣。當剛度不一樣時,設(shè)單元1和單元2的剛度分別為K1和K2,此時平衡方程有:
上面這個例子簡單的說明了有限元法的數(shù)學原理,有限元法的核心思想就是離散化分割為有限個單元,將無窮轉(zhuǎn)化為有限。實際問題的計算并沒有這么簡單,但思想是類似的,對力學問題來說,總可以列出形如下所示的外力和內(nèi)力平衡的方程:
更深入詳細的資料可以參看我前面共享的資料。
附件:http://pan.baidu.com/s/1eS5dCI2
相關(guān)標簽搜索：一個說明有限元法數(shù)學原理的簡單實例 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓
