有限元法概念意義與應(yīng)用(編輯修改稿)
2025-07-21 02:09 本頁面
【文章內(nèi)容簡(jiǎn)介】 算域上總體的基函數(shù)可以看為由每個(gè)單元基函數(shù)組成的，則整個(gè)計(jì)算域內(nèi)的解可以看作是由所有單元上的近似解構(gòu)成。在河道數(shù)值模擬中，常見的有限元計(jì)算方法是由變分法和加權(quán)余量法發(fā)展而來的里茲法和伽遼金法、最小二乘法等。根據(jù)所采用的權(quán)函數(shù)和插值函數(shù)的不同，有限元方法也分為多種計(jì)算格式。從權(quán)函數(shù)的選擇來說，有配置法、矩量法、最小二乘法和伽遼金法，從計(jì)算單元網(wǎng)格的形狀來劃分，有三角形網(wǎng)格、四邊形網(wǎng)格和多邊形 網(wǎng)格，從插值函數(shù)的精度來劃分，又分為線性插值函數(shù)和高次插值函數(shù)等。不同的組合 同樣構(gòu)成不同的有限元計(jì)算格式。對(duì)于權(quán)函數(shù)，伽遼金(Galerkin)法是將權(quán)函數(shù)取為逼近函數(shù)中的基函數(shù)；最小二乘法是令權(quán)函數(shù)等于余量本身，而內(nèi)積的極小值則為對(duì)代求系數(shù)的平方誤差最?。辉谂渲梅ㄖ?，先在計(jì)算域 內(nèi)選取N個(gè)配置點(diǎn)。令近似解在選定的N個(gè)配置點(diǎn)上嚴(yán)格滿足微分方程，即在配置點(diǎn)上令方程余量為0。插值函數(shù)一般由不同次冪的多項(xiàng)式組成，但也有采用三角函數(shù)或指數(shù)函數(shù)組成的乘積表示，但最常用的多項(xiàng)式插值函數(shù)。有限元插值函數(shù)分為兩大類，一類只要求插值多項(xiàng)式本身在插值點(diǎn)取已知值，稱為拉格朗日(Lagrange)多項(xiàng)式插值；另一種不僅要求插值多項(xiàng)式本身，還要求它的導(dǎo)數(shù)值在插值點(diǎn)取已知值，稱為哈密特(Hermite)多項(xiàng)式插值。單元坐標(biāo)有笛卡爾直角坐標(biāo)系和無因次自然坐標(biāo)，有對(duì)稱和不對(duì)稱等。常采用的無因次坐標(biāo)是一種局部坐標(biāo)系，它的定義取決于單元的幾何形狀，一維看作長(zhǎng)度比，二維看作面積比，三維看作體積比。在二維有限元中，三角形單元應(yīng)用的最早，近來四邊形等參元的應(yīng)用也越來越廣。對(duì)于二維三角形和四邊形電源單元，常采用的插值函數(shù)為有Lagrange插值直角坐標(biāo)系中的線性插值函數(shù)及二階或更高階插值函數(shù)、面積坐標(biāo)系中的線性插值函數(shù)、二階或更高階插值函數(shù)等。三 有限元單元法的基本思路彈性力學(xué)解法的問題彈性力學(xué)解法的問題在于：不論是應(yīng)力函數(shù)解法數(shù)解法、扭轉(zhuǎn)函數(shù)解法、撓曲函數(shù)解法、還是基于最小勢(shì)能原還是基于最小勢(shì)能原理的瑞利－李茲等方法，其困難在于如何給出一個(gè)在全求解區(qū)給出一個(gè)在全求解區(qū)域上均成立的試探函數(shù)。在有限單元法里在有限單元法里，這個(gè)問題通過定義分片插值的位移或應(yīng)力函數(shù)得到了巧妙的解決。對(duì)于任意單元對(duì)于任意單元(i,j,m)，以結(jié)點(diǎn)位移以結(jié)點(diǎn)位移(u,u ,u )為待定系數(shù)，可以給出該單元的插值函數(shù)： 線性代數(shù)方程組的求解在數(shù)學(xué)上是極其容易的。 也就是說有限元法通過單元離散和最小勢(shì)能原理小勢(shì)能原理，避開了微分方程直接求避微分方程直接求解在數(shù)學(xué)上的困難，把定解條件下的微分方程組的求解巧妙地轉(zhuǎn)化為線性方程組的運(yùn)算，實(shí)現(xiàn)了任何復(fù)雜彈性力學(xué)問題輕易分析計(jì)算。四 有限元單元法求解問題的的基本步驟1. (1)建立積分方程，根據(jù)變分原理或方程余量與權(quán)函數(shù)正交化原理，建立與微分方程初邊值問題等價(jià)的積分表達(dá)式，這是有限元法的出發(fā)點(diǎn)。 (2) 區(qū)域單元剖分，根據(jù)求解區(qū)域的形狀及實(shí)際問題的物理特點(diǎn)，將區(qū)域剖分為若干相互連接、不重疊的單元。區(qū)域單元?jiǎng)澐质遣捎糜邢拊椒ǖ那捌跍?zhǔn)備工作，這部分工作量比較大，除了給計(jì)算單元和節(jié)點(diǎn)進(jìn)行編號(hào)和確定相互之間的關(guān)系之外，還要表示節(jié)點(diǎn)的位置坐標(biāo)，同時(shí)還需要列出自然邊界和本質(zhì)邊界的節(jié)點(diǎn)序號(hào)和相應(yīng)的邊界值。3) 確定單元基函數(shù)，根據(jù)單元中節(jié)點(diǎn)數(shù)目及對(duì)近似解精度的要求，選擇滿足一定插值條 件的插值函數(shù)作為單元基函數(shù)。有限元方法中的基函數(shù)是在單元中選取的，由于各單元具有規(guī)則的幾何形狀，在選取基函數(shù)時(shí)可遵循一定的法則。 (4) 單元分析：將各個(gè)單元中的求解函數(shù)用單元基函數(shù)的線性組合表達(dá)式進(jìn)行逼近；再將近似函數(shù)代入積分方程，并對(duì)單元區(qū)域進(jìn)行積分，可獲得含有待定系數(shù)(即單元中各節(jié)點(diǎn)的參數(shù)值)的代數(shù)方程組，稱為單元有限元方程。 (5) 總體合成：在得出單元有限元方程之后，將區(qū)域中所有單元有限元方程按一定法則進(jìn)行累加，形成總體有限元方程。 (6) 邊界條件的處理：一般邊界條件有三種形式，分為本質(zhì)邊界條件(狄里克雷邊界條件)、自然邊界條件(黎曼邊界條件)、混合邊界條件(柯西邊界條件)。對(duì)于自然邊界條件， 一般在積分表達(dá)式中可自動(dòng)得到滿足。對(duì)于本質(zhì)邊界條件和混合邊界條件，需按一定法 則對(duì)總體有限元方程進(jìn)行修正滿足。 (7) 解有限元方程：根據(jù)邊界條件修正的總體有限元方程組，是含所有待定未知量的封閉方程組，采用適當(dāng)?shù)臄?shù)值計(jì)算方
電大資料相關(guān)推薦
鄂ICP備17016276號(hào)-1