有限元法講稿

《有限元法講稿》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。

1、平面問題的有限單元法有限單元法是隨著計算機(jī)的出現(xiàn)而發(fā)展起來的一種有效數(shù)值計算方法，有限單元法出現(xiàn)于40年代，被應(yīng)用于飛機(jī)結(jié)構(gòu)分析，有限元這個術(shù)語是1956年首先使用的。目前已廣泛地用于工程結(jié)構(gòu)的力學(xué)分析中。第一節(jié)基本概念一、實質(zhì)理想化連續(xù)體―――――――單元集合體（解析模擬、逼近求解區(qū)域）無限自由度有限個自由度有限單元法首先把結(jié)構(gòu)劃分成許多單元,在一定的簡化假設(shè)前提下,研究單元的力學(xué)特性，即單元分析；然后把各單元綜合起來,把局部的力學(xué)特性擴(kuò)展到整體,即整體分析37；最后導(dǎo)出一組以結(jié)構(gòu)結(jié)點位移為未知量的代數(shù)方程組。通過求解方程組而得到單元的結(jié)點位移值，就可近似計算出結(jié)構(gòu)任意一點的受力

2、狀態(tài)。這種以結(jié)點位移為基本未知量的計算方法稱為有限單元位移法。二、理論基礎(chǔ)彈性力學(xué)：變分原理能量原理基本方程：幾何方程、物理方程1.平面問題的幾何方程37這就是彈性力學(xué)平面問題的幾何方程，它給出了某一點的位移與該點應(yīng)變之間的關(guān)系。反映了變形協(xié)調(diào)關(guān)系。2.平面問題的物理方程即：給出了力與變形之間的關(guān)系，稱為平面問題物理方程，是針對平面應(yīng)力問題推導(dǎo)出的。對于平面應(yīng)變問題，只需將公式中的換成，把換成即可。這樣，彈性矩陣就變?yōu)椋?7（2.8）這就是適用于平面應(yīng)變問題的彈性矩陣。3.彈性體的能量原理（1）應(yīng)變能在彈性范圍內(nèi)，對于平面問題，某個面積A的應(yīng)變能可以用下式表示寫成矩陣的形式為37（

3、2）外力勢能外力勢能用矩陣的形式可表示為式中??作用在彈性體i點的外力分量；??i點的位移分量（3）彈性體的總勢能彈性體在外力作用的總勢能定義為應(yīng)變能和荷載勢能之和，即（4）最小勢能原理37單元的眾多的結(jié)點位移中,須滿足條件：即的一組位移才是真正的位移。最小勢能原理反映的是平面體的平衡方程。4.變分原理(1)基本概念設(shè)一泛函Π,它由積分形式確定,一般形式為:Π是隨函數(shù)u變化而變化的函數(shù)，稱為u的泛函。37（2）歐拉公式設(shè)一維泛函使此泛函取得極值的函數(shù)y(x)必定滿足下列方程式(推導(dǎo)從略)其中:;;;;這就是歐拉方程。滿足歐拉方程的函數(shù)y(x)能使泛函取極值。換言之,欲求能使泛函取極

4、值的函數(shù)y(x)就必須求解歐拉方程這一偏微分方程。(3)(位移)變分的近似解法＿＿瑞雷-里茲法37最小勢能原理給彈性力學(xué)問題提供了這樣一種近似解法：設(shè)定一個位移分量的表達(dá)式,使其滿足位移邊界條件和連續(xù)可微條件(這是容許位移所要滿足的條件),并包含若干待定系數(shù),用位移分量來表示變形能。這樣,總勢能即由位移和變形的泛函變?yōu)檫@些待定系數(shù)的函數(shù)。然后根據(jù)極值條件方程來確定這些待定系數(shù),從而求得近似解。這就是瑞雷法。里茲把位移函數(shù)假定為:其中:＿＿＿待定系數(shù)；＿＿＿假定的某種函數(shù)。這些函數(shù)可以是冪函數(shù)()也可以是其他合適的函數(shù)。把這些位移函數(shù)代入總勢能泛函,總勢能泛函即由位移的泛函變?yōu)檫@些待

5、定系數(shù)的函數(shù)。然后根據(jù)極值條件來確定這些待定系數(shù)以求得近似解。這就是瑞雷-里茲法。37這2n個方程正好解出這2n個待定系數(shù),從而求得位移函數(shù)的近似解。(4)用有限單元法解微分方程37三、分析的一般過程圖1.1（a）（b）1、劃分單元―――連續(xù)體的離散化構(gòu)成了以單元的集合體來近似代替原來的連續(xù)體的有限元分析計算模型2、建立位移模型373、單元分析單元的剛度方程可表示為，（1.1）式中：??單元結(jié)點力向量；??就是單元的剛度矩陣；它反映了單元結(jié)點力與單元結(jié)點位移之間的關(guān)系。單元分析的結(jié)果就建立單元的單元剛度方程和求出單元剛度矩陣。4、整體分析整個結(jié)構(gòu)的剛度方程可以由最小勢能原理來建立。

6、總剛度方程為（1.2）式中：??結(jié)構(gòu)的等效結(jié)點荷載向量；??結(jié)構(gòu)的總剛度矩陣；??結(jié)構(gòu)的結(jié)點位移向量，由單元的結(jié)點位移向量集成；整體分析的目的就是建立整個結(jié)構(gòu)的剛度方程，以求解未知的結(jié)點位移及計算單元應(yīng)力。375、求解節(jié)點力位移6、通過位移求解應(yīng)力、應(yīng)變四、單元類型把結(jié)構(gòu)離散成有限個單元時，可以選擇不同的單元形狀，圖1.2列出了工程中常用的幾種單元的形狀,自左至右分別是：三結(jié)點三角形單元,六結(jié)點三角形單元、四結(jié)點矩形單元、四結(jié)點等參數(shù)單元和八結(jié)點曲邊等參數(shù)單元。這些單元的特性和分析方法將在下面討論。圖1.2常用平面單元的形狀37第一節(jié)常應(yīng)變?nèi)切螁卧?、位移模式和形函?shù)圖5.3.1

7、37設(shè)、為單元中任一點沿軸和軸方向的位移，假定單元中各點的位移成線性變化，即（3.1）這就是我們所假定的三角形單元的線性位移模式。把上式寫成矩陣形式，即：（3.1a）或?qū)懗桑?.2）當(dāng)時,37當(dāng)時當(dāng)時同理（3.3a）（3.3b）求解以上兩式的方程組可得待定系數(shù)，即得到。這里略去推導(dǎo)過程中的繁瑣運算，最后可得到矩陣N為：（3.9）其中：A為單元面積；37（3.11）（3.12）37（3.13）令（3.14）、、稱為單元的形函數(shù)。矩陣Ｎ稱為單元的形函數(shù)矩陣。由式（3.9）