有限法的基礎(chǔ)知識(shí)

有限法的基礎(chǔ)知識(shí)第1頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月有限元的基礎(chǔ)理論3.1加權(quán)余量法3.2里茲方法3.3虛功原理3.4最小位能原理和最小余能原理第2頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月基礎(chǔ)知識(shí)——1.變分問(wèn)題在科學(xué)研究和工程實(shí)踐中，常常需要確定某一函數(shù)的極大值和極小值，這類計(jì)算已為大家所熟悉。但是，實(shí)際中常常需要確定另一類特殊的量（即所謂泛函）的極大值和極小值問(wèn)題，這就是變分法要處理的問(wèn)題。舉例：最速降線問(wèn)題。通俗的講泛函就是函數(shù)的函數(shù)，是變量與函數(shù)的關(guān)系。變分問(wèn)題就是泛函求極值的問(wèn)題。第3頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月2.微分方程的等效積分形式工程或物理學(xué)中的許多問(wèn)題，通常是以未知場(chǎng)函數(shù)應(yīng)滿足的微分方程和邊界條件的形式提出的，一般地表示為未知函數(shù)u應(yīng)滿足微分方程組域可以是體積域、面積域等。同時(shí)未知函數(shù)u還應(yīng)滿足邊界條件，是域的邊界。第4頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月由于微分方程組在域內(nèi)中的每一點(diǎn)都必須為零，因此就有是函數(shù)向量，它是一組和微分方程個(gè)數(shù)相等的任意函數(shù)。其中*式是與微分方程組完全等效的積分形式。可以證明，若積分方程*對(duì)于任意的都成立，則微分方程必然在域內(nèi)任意一點(diǎn)都滿足。第5頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月同理，假如邊界條件也同時(shí)在邊界上每一點(diǎn)都得到滿足，則對(duì)于任一組函數(shù)，下式應(yīng)當(dāng)成立。上式即為原微分方程的等效積分形式。問(wèn)題是有限元法不是在整個(gè)求解域上假設(shè)近似函數(shù)，而是在各個(gè)單元上分片假設(shè)近似函數(shù)。第6頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月等效積分的“弱”形式通過(guò)適當(dāng)提高對(duì)任意函數(shù)的連續(xù)性要求，以降低對(duì)微分方程場(chǎng)函數(shù)的u的連續(xù)性要求所建立的等效積分形式稱為微分方程的等效積分“弱”形式。從形式上看“弱”形式對(duì)函數(shù)u的連續(xù)性降低了，但對(duì)實(shí)際物理問(wèn)題卻常常較原始的微分方程更逼近真正解，原因在于微分方程往往對(duì)解提出了過(guò)分平滑的要求。第7頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月3.1加權(quán)余量法

（WRM--Weightedresidualmethod）加權(quán)余量法：指采用使余量的加權(quán)函數(shù)為零求得微分方程近似解的方法。加權(quán)余量法是求解微分方程近似解的一種有效的方法。顯然，任何獨(dú)立的完全函數(shù)集都可以作為權(quán)函數(shù)。按照對(duì)權(quán)函數(shù)的不同選擇得到不同的加權(quán)余量計(jì)算方法，主要有：配點(diǎn)法、子域法、最小二乘法、力矩法和伽遼金法。其中伽遼金法的精度最高。第8頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月3.2里茲方法里茲方法：如果微分方程具有線性和自伴隨的性質(zhì)，那么它不僅可以建立它的等效積分形式，并利用加權(quán)余量法求其近似解，而且還可以建立與之相等效的變分原理，從而得到的另一種近似求解方法。第9頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月3.3虛功原理（平衡方程和幾何方程的等效積分“弱”形式）變形體的虛功原理：變形體中任意滿足平衡的力系在任意滿足協(xié)調(diào)條件的變形狀態(tài)上作的虛功等于零，即體系外力的虛功與內(nèi)力的虛功之和等于零。虛功原理是虛位移原理和虛應(yīng)力原理的總稱。他們都可以認(rèn)為是與某些控制方程相等效的積分“弱”形式。第10頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月虛位移原理虛位移原理是平衡方程和力的邊界條件的等效積分的“弱”形式；虛位移原理的力學(xué)意義：如果力系是平衡的，則它們?cè)谔撐灰坪吞搼?yīng)變上所作的功的總和為零。反之，如果力系在虛位移（及虛應(yīng)變）上所作的功的和等于零，則它們一定滿足平衡方程。所以，虛位移原理表述了力系平衡的必要而充分條件。一般而言，虛位移原理不僅可以適用于線彈性問(wèn)題，而且可以用于非線性彈性及彈塑性等非線性問(wèn)題。但是否適用所有的問(wèn)題呢？第11頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月虛應(yīng)力原理虛應(yīng)力原理是幾何方程和位移邊界條件的等效積分“弱”形式。虛應(yīng)力原理的力學(xué)意義：如果位移是協(xié)調(diào)的（即在內(nèi)部連續(xù)可導(dǎo)），則虛應(yīng)力和虛邊界約束反力在他們上面所作的功的總和為零。反之，如果上述虛力系在他們上面所作的功的和為零，則它們一定是滿足協(xié)調(diào)的。所以，虛應(yīng)力原理表述了位移協(xié)調(diào)的必要而充分條件。虛應(yīng)力原理可以應(yīng)用于線彈性以及非線性彈性等不同的力學(xué)問(wèn)題。第12頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月但是必須指出，無(wú)論是虛位移原理還是虛應(yīng)力原理，他們所依賴的幾何方程和平衡方程都是基于小變形理論的，所以他們不能直接應(yīng)用于基于大變形理論的力學(xué)問(wèn)題。第13頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月3.4最小位能原理和最小余能原理明確：最小位能原理建立在虛位移原理基礎(chǔ)上，而最小余能原理建立在虛應(yīng)力原理基礎(chǔ)上。最小位能原理是指在所有可能位移中，真實(shí)位移使系統(tǒng)總位能取最小值。總位能是指彈性體變形位能和外力位能之和。最小余能原理是指在所有的應(yīng)力中，真實(shí)應(yīng)力使系統(tǒng)的總余能取最小值?？傆嗄苁侵笍椥泽w余能和外力余能總和。第14頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月意義一般而言，利用最小位能原理求得位移近似解的彈性變形能是精確解變形能的下界，即近似的位移場(chǎng)在總體上偏小，也就是說(shuō)結(jié)構(gòu)的計(jì)算模型顯得偏于剛硬；而利用最小余能原理求得的應(yīng)力近似解的彈性余能是精確解余能的上界，即近似的應(yīng)力解在總體上偏大，結(jié)構(gòu)的計(jì)算模型偏于柔軟。當(dāng)分別利用這兩個(gè)極值原理求解同一問(wèn)題時(shí)，我們將獲得這個(gè)問(wèn)題的上界和下界，可以較準(zhǔn)確地估計(jì)所得近似解的誤差，這對(duì)工程計(jì)算具有實(shí)際意義。第15頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月第四節(jié)有限元法的求解方法與步驟

4.1有限元的基本思想有限元是一種以計(jì)算機(jī)為手段，通過(guò)離散化將研究對(duì)象變換成一個(gè)與原始結(jié)構(gòu)近似的數(shù)學(xué)模型，再經(jīng)過(guò)一系列規(guī)范化的步驟以求解應(yīng)力、應(yīng)變、位移等參數(shù)的數(shù)值計(jì)算方法。

第16頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月什么是離散化呢？所謂離散化是將一個(gè)連續(xù)體分割成若干個(gè)通過(guò)節(jié)點(diǎn)相連的單元，這樣一個(gè)有無(wú)限個(gè)自由度的結(jié)構(gòu)就變換成一個(gè)具有有限個(gè)自由度的近似結(jié)構(gòu)。該過(guò)程還包括對(duì)單元和節(jié)點(diǎn)進(jìn)行編碼以及局部坐標(biāo)系和整體坐標(biāo)系的確定。第17頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月什么是單元和節(jié)點(diǎn)呢？在有限元法中，將求解的工程結(jié)構(gòu)看成是由許多小的、彼此用點(diǎn)聯(lián)接的基本構(gòu)件，如桿和梁、板和殼組成，這些基本構(gòu)件稱為單元，單元與單元之間的聯(lián)接點(diǎn)稱為節(jié)點(diǎn)。第18頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月4.2有限元法在結(jié)構(gòu)工程中的應(yīng)用

根據(jù)研究對(duì)象的不同，有限元法中采用的單元形式也不相同，常見(jiàn)的有以下幾種：1.桁架桿單元：主要應(yīng)用于受軸向力作用的桿和桿系，如桁架結(jié)構(gòu)；2.

剛架桿單元：用于梁及剛架結(jié)構(gòu)分析；3.三角形平面單元：主要用于彈性力學(xué)中平面應(yīng)力問(wèn)題和平面應(yīng)變問(wèn)題的有限元分析；第19頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月4.三棱圓環(huán)單元：用于軸對(duì)稱問(wèn)題的有限元分析；5.等參數(shù)單元：用于一些具有曲線輪廓的復(fù)雜結(jié)構(gòu)。其特點(diǎn)是能簡(jiǎn)化復(fù)雜結(jié)構(gòu)的單元?jiǎng)澐止ぷ?，又能滿足同樣精度的要求時(shí)大大減少使用的單元數(shù)，成功地解決許多二維和三維彈性力學(xué)問(wèn)題。第20頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月4.3有限元法求解問(wèn)題的基本步驟1.結(jié)構(gòu)離散化對(duì)整個(gè)結(jié)構(gòu)進(jìn)行離散化，將其分割成若干個(gè)單元，單元間彼此通過(guò)節(jié)點(diǎn)相連；

2.求出各單元的剛度矩陣

是由單元節(jié)點(diǎn)位移量求單元節(jié)點(diǎn)力向量的轉(zhuǎn)移矩陣，其關(guān)系式為第21頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月3.集成總體剛度矩陣[K]并寫出總體平衡方程：總體剛度矩陣[K]是由整體節(jié)點(diǎn)位移向量求整體節(jié)點(diǎn)力向量的轉(zhuǎn)移矩陣，其關(guān)系式為，此即為總體平衡方程。

第22頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月確定總體剛度矩陣的方法有三種：

1）直接利用總體剛度系數(shù)的定義在求出整體結(jié)構(gòu)中各節(jié)點(diǎn)力與節(jié)點(diǎn)位移關(guān)系的基礎(chǔ)上獲得總體剛度矩陣。此方法旨在簡(jiǎn)單情況下才能采用。2）集成法將整體坐標(biāo)下的單元?jiǎng)偠染仃囘M(jìn)行迭加而得。這里所說(shuō)的迭加不是簡(jiǎn)單的相加，而是將下角標(biāo)相同的總體剛度系數(shù)相加，然后按總碼的順序?qū)μ?hào)入座。

第23頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月3）利用節(jié)點(diǎn)間的剛度系數(shù)直接寫出總體剛度矩陣總體剛度矩陣對(duì)角線上的剛度系數(shù)等于在節(jié)點(diǎn)i匯交的幾個(gè)單元的剛度系數(shù)之和；非對(duì)角線上的剛度系數(shù)等于聯(lián)結(jié)節(jié)點(diǎn)i與節(jié)點(diǎn)j間幾個(gè)單元的剛度系數(shù)之和。

第24頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月4.引入支撐條件，求出各節(jié)點(diǎn)的位移節(jié)點(diǎn)的支撐條件有兩種：一種是節(jié)點(diǎn)n沿某個(gè)方向的位移為零，另一種是節(jié)點(diǎn)n沿某個(gè)方向的位移為一給定值。5.求出各單元內(nèi)的應(yīng)力和應(yīng)變。

第25頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月第五節(jié)有限元法求解實(shí)例

下面通過(guò)一維單元的簡(jiǎn)單實(shí)例說(shuō)明具體過(guò)程。

例1一根由兩段組成的階梯軸，一端固定，另一端承受一個(gè)軸向載荷F3。這兩段的橫截面積分別為和，長(zhǎng)度分別為和，彈性模量分別為和，如圖1-1所示。求出這兩段的應(yīng)力和應(yīng)變。已知數(shù)據(jù)分別為F3=100N

第26頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月圖1-11

①2②3A(1)E(2)A(2)E(2)L(1)

L(2)

①②2Φ1F2F3Φ2Φ3F1F3第27頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月解：1）離散化把這根階梯軸看成是由兩個(gè)單元組成的，節(jié)點(diǎn)選在截面積突變處，兩個(gè)單元的連接處是一個(gè)節(jié)點(diǎn)，該階梯軸的兩端視為另外兩個(gè)節(jié)點(diǎn)，所以整個(gè)結(jié)構(gòu)共有三個(gè)節(jié)點(diǎn)。這根軸是一維結(jié)構(gòu)，并只受軸向載荷，因此各單元內(nèi)只有軸向位移。三個(gè)節(jié)點(diǎn)位置的位移量分別記為、、。在整個(gè)結(jié)構(gòu)中節(jié)點(diǎn)載荷及節(jié)點(diǎn)位移均用大寫字母標(biāo)記，其角標(biāo)為節(jié)點(diǎn)在總體結(jié)構(gòu)中的編碼，簡(jiǎn)稱總碼。

第28頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月2）求單元?jiǎng)偠染仃囅旅娣治瞿车冉孛鎲卧╡）。當(dāng)兩端分別承受兩個(gè)軸向力和作用時(shí)的位移情況。根據(jù)材料力學(xué)的知識(shí)可知，在兩端節(jié)點(diǎn)i、j處的位移量和與軸向力和的關(guān)系式為第29頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月注意在分析單元?jiǎng)偠染仃嚂r(shí)，載荷F和位移等參數(shù)的上角標(biāo)為該單元的編碼，下角標(biāo)為該單元內(nèi)節(jié)點(diǎn)的局部編碼。上兩式可寫成：或簡(jiǎn)寫為：第30頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月式中—為單元?jiǎng)偠染仃嚮騿卧匦跃仃?，其階數(shù)等于單元中所包含的節(jié)點(diǎn)數(shù)；

—為單元節(jié)點(diǎn)力向量（列陣）；

第31頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月

——為單元節(jié)點(diǎn)位移向量（列陣），也為單元自由度列陣；將單元?jiǎng)偠染仃嚫膶懗删仃嚨臉?biāo)準(zhǔn)形式，則

第32頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月（1-5）矩陣中任意一個(gè)元素都稱為單元?jiǎng)偠认禂?shù)，它表示該單元內(nèi)除節(jié)點(diǎn)j產(chǎn)生單位位移外，其余各節(jié)點(diǎn)的位移均為零時(shí)在節(jié)點(diǎn)i處所引起的載荷Fi。

第33頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月3)總體剛度矩陣的集成和總體平衡方程的寫出

該階梯軸上三個(gè)節(jié)點(diǎn)位移和三個(gè)節(jié)點(diǎn)軸向力分別組成該整體結(jié)構(gòu)節(jié)點(diǎn)位移向量和節(jié)點(diǎn)軸向力

同理，這兩向量間的轉(zhuǎn)換關(guān)系可表示為或

第34頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月式中的轉(zhuǎn)移矩陣稱為總體剛度矩陣或總體特性矩陣，其階數(shù)等于總體結(jié)構(gòu)中的節(jié)點(diǎn)總數(shù)。[K]中的元素稱為總體剛度系數(shù)，它表示在整體結(jié)構(gòu)中除了節(jié)點(diǎn)j產(chǎn)生單位位移外，其余各節(jié)點(diǎn)的位移均為零時(shí)在節(jié)點(diǎn)i處所引起的載荷Fi。

第35頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月求出總體剛度矩陣時(shí)進(jìn)行總體分析的主要任務(wù)是一旦獲得總體剛度矩陣，可以很容易地寫出總體平衡方程。求總體剛度矩陣[K]的方法主要有兩種：一是直接法，即根據(jù)總體剛度系數(shù)的定義求解；另一種方法是集成法，即由各單元?jiǎng)偠染仃嚽罂傮w剛度矩陣。下面分別說(shuō)明。第36頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月根據(jù)剛度系數(shù)的定義，當(dāng)本結(jié)構(gòu)中的節(jié)點(diǎn)2和節(jié)點(diǎn)3位移量均為零時(shí)，要使節(jié)點(diǎn)1產(chǎn)生單位位移，在節(jié)點(diǎn)1處所需施加的載荷為，此即為K11；當(dāng)節(jié)點(diǎn)1、3固定，節(jié)點(diǎn)2產(chǎn)生單位位移時(shí)，在節(jié)點(diǎn)1處所引起的載荷為，此即為K12；

（1）直接法求總體剛度矩陣[K]

第37頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月當(dāng)節(jié)點(diǎn)1、2固定，節(jié)點(diǎn)3產(chǎn)生單位位移時(shí)，在節(jié)點(diǎn)1處不會(huì)引起載荷，因此K13=0；當(dāng)節(jié)點(diǎn)1、3固定，節(jié)點(diǎn)2產(chǎn)生單位位移時(shí)，要在節(jié)點(diǎn)2施加的載荷為，此即為K22；還可以按此方法依次寫出其余各總體剛度系數(shù)。

第38頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月因此，總體剛度矩陣為

這種方法具有概念清晰的特點(diǎn)，但是在分析復(fù)雜結(jié)構(gòu)時(shí)運(yùn)算極其復(fù)雜，因而限制了它的應(yīng)用。

第39頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月（2）用集成法求總體剛度矩陣[K]

這種方法從單元?jiǎng)偠染仃嚦霭l(fā)，根據(jù)迭加原理，利用剛度系數(shù)集成的方法獲得總體剛度矩陣。這樣，首先要寫出各單元的剛度矩陣。由式（1-4），當(dāng)單元(e)分別為(1)和(2)時(shí)，兩個(gè)單元的剛度矩陣分別為

第40頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月注意雖然結(jié)構(gòu)的總變形是由兩個(gè)單元變形的迭加，但是總體剛度矩陣[K]并不是兩個(gè)單元?jiǎng)偠染仃嚭偷暮?jiǎn)單迭加。

在進(jìn)行這項(xiàng)工作之前，一定要分清節(jié)點(diǎn)的兩種編碼方式：一種為節(jié)點(diǎn)的局部編碼，單元（e）的兩個(gè)局部碼分別為和；另一種為節(jié)點(diǎn)總碼，即對(duì)結(jié)構(gòu)中全部節(jié)點(diǎn)進(jìn)行統(tǒng)一的編碼。在本例中，階梯軸的三個(gè)節(jié)點(diǎn)總碼分別記為1、2、3。

第41頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月該階梯軸兩種編碼的對(duì)應(yīng)關(guān)系為

局部碼

總碼

1223一定要注意，在單元?jiǎng)偠染仃囍懈髟匕淳植看a排列，而總體剛度矩陣[K]中，各元素按總碼排列。

第42頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月集成[K]的步驟為：

（1）將原單元?jiǎng)偠染仃囍械母飨禂?shù)總碼進(jìn)行標(biāo)記，則

第43頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月（2）將角標(biāo)相同的系數(shù)相加，并按總碼的順序排列，則總體剛度矩陣為總體剛度方程為：

從上式可以看出，用兩種方法獲得的總體剛度矩陣相同。

第44頁(yè)，課件共52頁(yè)，創(chuàng)作于2023年2月（4）引入支撐條件，計(jì)算節(jié)點(diǎn)位移

上式中的未知量仍不能求出，因?yàn)閇K]是一個(gè)奇異矩陣，必須引入支撐條件。在本