彈性力學(xué)平面問題的有限單元法

第6章平面問題的有限單元法目錄1目錄主要內(nèi)容概述§6-1基本量及基本方程的矩陣表示§6-2有限單元法的概念§6-3單元的位移模式與解答的收斂性§6-4單元的應(yīng)變列陣和應(yīng)力列陣§6-5單元的結(jié)點力列陣與勁度矩陣§6-6荷載向結(jié)點移置單元的結(jié)點荷載列陣2目錄主要內(nèi)容§6-7結(jié)構(gòu)的整體分析結(jié)點平衡方程組§6-8解題的具體步驟單元的劃分§6-9計算成果的整理§6-10計算實例31.有限元法（FiniteElementMethod）

2.FEM的特點

概述（1）具有通用性和靈活性。首先將連續(xù)體變換為離散化結(jié)構(gòu)，然后再利用分片插值技術(shù)與虛功原理或變分方法進行求解。簡稱FEM，是彈性力學(xué)的一種近似解法。43.FEM簡史

（2）對同一類問題，可以編制出通用程序，應(yīng)用計算機進行計算。（3）只要適當(dāng)加密網(wǎng)格，就可以達到工程要求的精度。1943年柯朗（Courant）第一次提出了FEM的概念。FEM是上世紀(jì)中期才出現(xiàn)，并得到迅速發(fā)展和廣泛應(yīng)用的一種數(shù)值解法。51970年后，F(xiàn)EM被引入我國，并很快地得到應(yīng)用和發(fā)展。1956年，特納（Turner）等人提出了FEM。

20世紀(jì)50年代，平面問題的FEM建立，并應(yīng)用于工程問題。1960年克勞夫（Clough）提出了FEM的名稱。20世紀(jì)60年代后，F(xiàn)EM應(yīng)用于各種力學(xué)問題和非線性問題，并得到迅速發(fā)展。65.本章介紹平面問題的FEM4.FEM的主要導(dǎo)出方法

應(yīng)用靜力方法或變分方法導(dǎo)出。僅敘述按位移求解的方法。且一般都以平面應(yīng)力問題來表示。7§6.1基本量和基本方程的矩陣表示

本章無特別指明，均表示為平面應(yīng)力問題的公式。

采用矩陣表示,可使公式統(tǒng)一、簡潔，且便于編制程序。8§6.1基本量和基本方程的矩陣表示基本物理量：體力:位移函數(shù):應(yīng)變:應(yīng)力:結(jié)點位移列陣:結(jié)點力列陣:面力:9§6.1基本量和基本方程的矩陣表示物理方程:FEM中應(yīng)用的方程：幾何方程:其中D為彈性矩陣，對于平面應(yīng)力問題是:10§6.1基本量和基本方程的矩陣表示

--結(jié)點虛位移;--對應(yīng)的虛應(yīng)變。ij虛功方程:其中:在FEM中，用結(jié)點的平衡方程代替平衡微分方程，后者不再列出。11§6.2有限單元法的概念3.整體分析。

FEM的概念，可以簡述為：采用有限自由度的離散單元組合體模型去描述實際具有無限自由度的考察體，是一種在力學(xué)模型上進行近似的數(shù)值計算方法。

其理論基礎(chǔ)是分片插值技術(shù)與變分原理。1.將連續(xù)體變換為離散化結(jié)構(gòu)；

2.單元分析；FEM的分析過程：12§6.2有限單元法的概念結(jié)構(gòu)力學(xué)研究的對象是離散化結(jié)構(gòu)。如桁架，各單元（桿件）之間除結(jié)點鉸結(jié)外，沒有其他聯(lián)系（圖（a））。彈力研究的對象，是連續(xù)體（圖（b）)。1.

結(jié)構(gòu)離散化－－將連續(xù)體變換為離散化結(jié)構(gòu)13§6.2有限單元法的概念將連續(xù)體變換為離散化結(jié)構(gòu)（圖（c））：即將連續(xù)體劃分為有限多個、有限大小的單元，并使這些單元僅在一些結(jié)點處用絞連結(jié)起來，構(gòu)成所謂‘離散化結(jié)構(gòu)’。14§6.2有限單元法的概念

圖（c）與圖(a）相比，兩者都是離散化結(jié)構(gòu)；區(qū)別是，桁架的單元是桿件，而圖（c）的單元是三角形塊體（注意：三角形單元內(nèi)部仍是連續(xù)體）。例如：將深梁劃分為許多三角形單元，這些單元僅在角點用鉸連接起來。15§6.2有限單元法的概念2.單元分析每個三角形單元仍然假定為連續(xù)的、均勻的、各向同性的完全彈性體。因單元內(nèi)部仍是連續(xù)體，應(yīng)按彈性力學(xué)方法進行分析。取各結(jié)點位移

為基本未知量。然后對每個單元,分別求出各物理量,并均用來表示。16§6.2有限單元法的概念（1）應(yīng)用插值公式,由單元結(jié)點位移，求單元的位移函數(shù)這個插值公式稱為單元的位移模式，為：單元分析的主要內(nèi)容：17§6.2有限單元法的概念（4）應(yīng)用虛功方程，由單元的應(yīng)力，求出單元的結(jié)點力，表示為（3）應(yīng)用物理方程，由單元的應(yīng)變，求出單元的應(yīng)力，表示為（2）應(yīng)用幾何方程，由單元的位移函數(shù)d，

求出單元的應(yīng)變，表示為18§6.2有限單元法的概念－－單元對結(jié)點的作用力，與數(shù)值相同,方向相反，作用于結(jié)點。

--結(jié)點對單元的作用力，作用于單元，稱為結(jié)點力，以正標(biāo)向為正。19§6.2有限單元法的概念（5）將每一單元中的各種外荷載，按虛功等效原則移置到結(jié)點上，化為結(jié)點荷載，表示為

20§6.2有限單元法的概念

為已知值,是用結(jié)點位移表示的值。通過求解聯(lián)立方程，得出各結(jié)點位移值，從而求出各單元的應(yīng)變和應(yīng)力。

各單位移置到i結(jié)點上的結(jié)點荷載其中表示對圍繞i結(jié)點的單元求和；3.整體分析各單元對i結(jié)點的結(jié)點力作用于結(jié)點i上的力有：21§6.2有限單元法的概念

3.整體分析

2.對單元進行分析

1.將連續(xù)體變換為離散化結(jié)構(gòu)（1）單元的位移模式（2）單元的應(yīng)變列陣（4）單元的結(jié)點力列陣（5）單元的等效結(jié)點荷載列陣建立結(jié)點平衡方程組，求解各結(jié)點的位移。（3）單元的應(yīng)力列陣歸納起來，F(xiàn)EM分析的主要步驟：22§6.2有限單元法的概念思考題

1.桁架的單元為桿件，而平面體的單元為三角形塊體，在三角形內(nèi)仍是作為連續(xù)體來分析的。前者可用結(jié)構(gòu)力學(xué)方法求解，后者只能用彈性力學(xué)方法求解，為什么？2.在平面問題中，是否也可以考慮其它的單元形狀，如四邊形單元？23§6.3單元的位移模式與解答的收斂性應(yīng)用插值公式，可由求出位移。

首先必須解決：由單元的結(jié)點位移來求出單元的位移函數(shù)

FEM是取結(jié)點位移為基本未知數(shù)的。問題是如何求應(yīng)變、應(yīng)力。

這個插值公式表示了單元中位移的分布形式，因此稱為位移模式。24§6.3單元的位移模式與解答的收斂性

插值公式（a）在結(jié)點應(yīng)等于結(jié)點位移值。由此可求出

泰勒級數(shù)展開式中，低次冪項是最重要的。所以三角形單元的位移模式，可取為：（a）25§6.3單元的位移模式與解答的收斂性將式（a）按未知數(shù)歸納為:或用矩陣表示為:（b）26§6.3單元的位移模式與解答的收斂性N－－稱為形（態(tài)）函數(shù)矩陣。（c）27§6.3單元的位移模式與解答的收斂性

A為三角形的面積（圖示坐標(biāo)系中，按逆時針編號），有：其中:28§6.3單元的位移模式與解答的收斂性三結(jié)點三角形單元的位移模式，略去了2次以上的項，因而其誤差量級是且其中只包含了的1次項，所以在單元中的分布如圖（a）所示，的分布如圖（b）、（c）所示。(a)(b)(c)129§6.3單元的位移模式與解答的收斂性

所以當(dāng)單元趨于很小時，即時，為了使FEM之解逼近于真解。則為了保證FEM收斂性,位移模式應(yīng)滿足下列條件：

FEM中以后的一系列工作，都是以位移模式為基礎(chǔ)的。

30§6.3單元的位移模式與解答的收斂性

因為當(dāng)單元時，單元中的位移和應(yīng)變都趨近于基本量－－剛體位移和常量位移。（1）位移模式必須能反映單元的剛體位移。

（2）位移模式必須能反映單元的常量應(yīng)變。31§6.3單元的位移模式與解答的收斂性可見剛體位移項在式（a）中均已反映。與剛體位移相比，將式（a）寫成32§6.3單元的位移模式與解答的收斂性（3）位移模式應(yīng)盡可能反映位移的連續(xù)性。即應(yīng)盡可能反映原連續(xù)體的位移連續(xù)性。在三角形單元內(nèi)部，位移為連續(xù)；在兩單元邊界ij

上，之間均為線性變化，也為連續(xù)。對式（a）求應(yīng)變，得：可見常量應(yīng)變也已反映。33§6.3單元的位移模式與解答的收斂性（1）和（2）是必要條件，而加上（3）就為充分條件。為了保證FEM的收斂性：34§6.3單元的位移模式與解答的收斂性思考題1.應(yīng)用泰勒級數(shù)公式來選取位移模式，為什么必須從低次項開始選取？2.試考慮：將結(jié)構(gòu)力學(xué)解法引入到求解連續(xù)體的問題時，位移模式的建立是一個關(guān)鍵性工作，它使得單元(連續(xù)體)內(nèi)部的分析工作都有可能進行了。

35§6.4單元的應(yīng)變列陣和應(yīng)力列陣其中，單元中的位移函數(shù)用位移模式表示為36§6.4單元的應(yīng)變列陣和應(yīng)力列陣應(yīng)用幾何方程，求出單元的應(yīng)變列陣：37§6.4單元的應(yīng)變列陣和應(yīng)力列陣再應(yīng)用物理方程，求出單元的應(yīng)力列陣：B——稱為應(yīng)變矩陣，用分塊矩陣表示，S——稱為應(yīng)力轉(zhuǎn)換矩陣，寫成分塊形式為38§6.4單元的應(yīng)變列陣和應(yīng)力列陣

對于線性位移模式，求導(dǎo)后得到的應(yīng)變和應(yīng)力，均成為常量，因此，稱為常應(yīng)變（應(yīng)力）單元。應(yīng)變和應(yīng)力的誤差量級是其精度比位移低一階，且相鄰單元的應(yīng)力是跳躍式的。39§6.4單元的應(yīng)變列陣和應(yīng)力列陣思考題1.如果在位移模式中取到泰勒級數(shù)中的二次冪項，略去高階小量，試考慮位移、應(yīng)變和應(yīng)力的誤差量級。40§6.5單元的結(jié)點力列陣和勁度矩陣

現(xiàn)在來考慮其中一個單元：

在FEM中，首先將連續(xù)體變換為離散化結(jié)構(gòu)的模型。41§6.5單元的結(jié)點力列陣和勁度矩陣（2）單元與周圍的單元在邊界上已沒有聯(lián)系，只在結(jié)點互相聯(lián)系。（1）將作用于單元上的各種外荷載，按靜力等效原則移置到結(jié)點上去，化為等效結(jié)點荷載。故單元內(nèi)已沒有外荷載。42§6.5單元的結(jié)點力列陣和勁度矩陣假想將單元與結(jié)點i

切開，則：

其數(shù)值與相同，而方向相反。以沿正坐標(biāo)向為正。對單元而言，這是作用于單元上的“外力”。

結(jié)點作用于單元上的力，稱為結(jié)點力，單元作用于結(jié)點的力，為：43§6.5單元的結(jié)點力列陣和勁度矩陣按虛功方程，在虛位移上，外力的虛功等于應(yīng)力的虛功。而其內(nèi)部有應(yīng)力作用，考察已與結(jié)點切開后的單元，則此單元上作用有外力－－結(jié)點力，應(yīng)用虛功方程，求單元的結(jié)點力：44§6.5單元的結(jié)點力列陣和勁度矩陣假設(shè)發(fā)生一組結(jié)點虛位移則單元內(nèi)任一點（x,y）的虛位移為單元內(nèi)任一點（x,y）的虛應(yīng)變?yōu)榇胩摴Ψ匠蹋涸趩卧校饬ΓńY(jié)點力）在虛位移（結(jié)點虛位移）上的虛功，等于應(yīng)力在虛應(yīng)變上的虛功，即：45§6.5單元的結(jié)點力列陣和勁度矩陣其中與無關(guān)，故式(a)成為式(b)是由應(yīng)力求結(jié)點力的一般公式。因為是獨立的任意的虛位移，虛功方程對任意的均應(yīng)滿足，可得出代入

(b)46§6.5單元的結(jié)點力列陣和勁度矩陣式（c）是由結(jié)點位移求結(jié)點力的一般公式，－－稱為單元的勁度矩陣K其中：再將應(yīng)力公式代入上式，得（c）（d）47§6.5單元的結(jié)點力列陣和勁度矩陣對于三角形單元，B矩陣內(nèi)均為常數(shù)，有

代入B，D，得出k

如書中（6-37）及（6-38）所示。48§6.5單元的結(jié)點力列陣和勁度矩陣（1）是6×6的方陣，中每一個元素都表示單元各結(jié)點沿坐標(biāo)方向發(fā)生單位位移時所引起的結(jié)點力。（2）由反力互等定理，所以是對稱矩陣，以對角線為對稱軸。單元勁度矩陣k的性質(zhì)：（3）當(dāng)單元作剛體平移時，如三角形內(nèi)不產(chǎn)生應(yīng)力和應(yīng)變，結(jié)點力也為0。49§6.5單元的結(jié)點力列陣和勁度矩陣（4）由（3）可導(dǎo)出行列式。（5）

的元素與單元的形狀和方位等有關(guān)，但與單元的大小和剛體的平動以及作度轉(zhuǎn)動無關(guān)。即有：中每一行（或列）的元素之和為零（其中第1、3、5元素之和或2、4、6元素之和也為0）。50§6.5單元的結(jié)點力列陣和勁度矩陣（書中P.117頁），以直角三角形單元為例，計算了應(yīng)力轉(zhuǎn)換矩陣S和單元勁度矩陣。

從例題中可以看出，將單元邊界上的應(yīng)力向結(jié)點移置，化為作用于結(jié)點上的力，正好就是結(jié)點力。在FEM中，單元邊界之間的聯(lián)系和相互作用力，都向結(jié)點簡化，歸結(jié)成為結(jié)點的鉸結(jié)和結(jié)點力。

思考題例題試求出書中例題的位移模式。51§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣在FEM中，須將作用于單元中的外荷載向結(jié)點移置，化為等效結(jié)點荷載，52§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣（2）變形體靜力等效原則－－在任意的虛位移上，使原荷載與移置荷載的虛功相等。

1、等效原則（1）剛體靜力等效原則－－使原荷載與移置荷載的主矢量以及對同一點的主矩也相同。剛體靜力等效原則只從運動效應(yīng)來考慮，得出移置荷載不是唯一的解；變形體的靜力等效原則考慮了變形效應(yīng)，在一定的位移模式下，其結(jié)果是唯一的，且也滿足了前者條件的。所以在FEM中，采用變形體的靜力等效原則。53§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣假設(shè)發(fā)生一組結(jié)點虛位移，

2、集中力的移置公式

原荷載作用于單元中任一點M為單位厚度上的作用力；移置荷載作用于結(jié)點則M點的虛位移為54§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣使移置荷載的虛功等于原荷載的虛功：對于任意的虛位移，虛功方程都必須滿足，得：55§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣3、單元邊界上面力的移置公式

應(yīng)用式，將代之為并在邊界上積分，得:56§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣

應(yīng)用式，將代之為并對單元域A積分，得

4、單元內(nèi)體力的移置公式

當(dāng)位移模式為線性函數(shù)時，由虛功方程得出的移置荷載，與按剛體靜力等效原則得出的結(jié)點荷載相同。57§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣若干荷載的移置如下圖所示：58§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣若干荷載的移置如下圖所示：59§6.6荷載向結(jié)點移置單元的結(jié)點載荷列陣思考題1.試導(dǎo)出書中例題的荷載移置公式。

60§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組在單元分析中，從單元的結(jié)點位移→求位移分布→求應(yīng)變→求應(yīng)力→求結(jié)點力，為單元的內(nèi)力分析；外荷載移置到結(jié)點荷載，為單元的外力分析。

假設(shè)將結(jié)點i與周圍的單元切開，則圍繞i結(jié)點的每個單元對i結(jié)點有結(jié)點力(）的作用,也有外荷載移置的結(jié)點荷載（）的作用。下面考慮整體分析。61§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組對某一個單元，其中是對圍繞i結(jié)點的單元求和。

i

結(jié)點的平衡條件為

62§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組是單元結(jié)點的局部編號；

是整體結(jié)點的整體編號。

代入式，可表示為將式按整體結(jié)點編號排列，得整個結(jié)構(gòu)的平衡方程組。63§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組

考慮結(jié)構(gòu)的約束條件后，從式求出，就可以求出各單元的位移和應(yīng)力。－－整體結(jié)點位移列陣，－－整體結(jié)點荷載列陣，－－整體勁度矩陣。64§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組例1列出圖示結(jié)構(gòu)i

結(jié)點的平衡條件。65§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組結(jié)點i聯(lián)系四個單元①～④，對i點列方程：66§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組整體勁度矩陣是由單元勁度矩陣的元素集合而成，稱總剛的組裝。組裝原則：對號入座整體勁度矩陣的特點：（1）對稱性；（2）奇異性；（3）稀疏性；（4）非零元素成帶狀分布67§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組例2（見書中P.121）設(shè)有對角受壓的正方形薄板，荷載沿厚度均勻分布，為2N/m，試用有限單元法求解薄板內(nèi)的應(yīng)力。68§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組解：（1）建立有限單元法的計算模型：考慮對稱性：4個單元，6個結(jié)點平面應(yīng)力問題69§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組（2）列出整體結(jié)構(gòu)的結(jié)點平衡方程組K=FL70§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組（3）建立結(jié)構(gòu)的整體結(jié)點荷載列陣FL整體結(jié)點位移列陣：71§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組（4）求出各單元的勁度矩陣k72§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組應(yīng)用式（6-38），可得兩種單元的勁度矩陣都是：73§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組（5）求出整體勁度矩陣K建立結(jié)點的整體編碼與單元內(nèi)的局部編碼之間的關(guān)系：74§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組75§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組整體勁度矩陣的組裝：76§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組77§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組（6）結(jié)點平衡方程組（a）的簡化與求解由于有位移邊界條件u1=u2=u4=v4=v5=v6=0，未知的整體結(jié)點位移列陣簡化為：整體荷載列陣簡化為：78§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組整體平衡方程：79§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組求解后，得結(jié)點位移：80§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組（7）由位移求應(yīng)力：對于單元Ⅰ，Ⅱ，Ⅳ：對于單元Ⅲ：81§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組各單元的應(yīng)力如下：82§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組83§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組

例題4圖（a）所示的深梁，在跨中受集中力F的作用，若取試用有限單元法求解跨中的位移。84§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組解1.將圖（a）劃分網(wǎng)格，化為離散化結(jié)構(gòu)，如圖（b）所示。由于結(jié)構(gòu)具有對稱性，可取1/2部分進行分析，如圖（c）所示。F/285§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組2.

中，只有兩個未知結(jié)點位移其余的結(jié)點位移均為零。

未知的結(jié)點位移列陣是對應(yīng)的結(jié)點荷載列陣是3.下面我們直接來建立對應(yīng)于未知結(jié)點位移的平衡方程式，圖（c）86§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組4.對于三角形單元，按照結(jié)點的局部編號結(jié)點力一般公式是87§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組當(dāng)且結(jié)點的局部編號如圖時，單元的單元勁度矩陣均如書中所示。對于單元，結(jié)點的局部編號與整體編號的關(guān)系是將書中的k和結(jié)點編號代入式，有88§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組

其中由上式，得出

I單元中不存在，而89§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組對于Ⅱ單元，結(jié)點的局部編號與整體編號的關(guān)系是。再將書中的k代入式（c），得90§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組其中由上式，可得Ⅱ單元的結(jié)點力5.將各單元的結(jié)點力代入式（a）（b）得從上兩式解出結(jié)點位移值，91§6.7結(jié)構(gòu)的整體分析結(jié)點平衡方程組92§6.8解題的具體步驟單元的劃分有限單元法的具體計算步驟：

1、劃分單元網(wǎng)格，對單元和結(jié)點編號。2、選定直角坐標(biāo)系，按程序要求填寫和輸入有關(guān)信息。單元內(nèi)的ijm的局部編號應(yīng)按書中規(guī)定的右手規(guī)則編號。否則會使三角形的面積出現(xiàn)負號等問題。93§6.8解題的具體步驟單元的劃分

3、使用已編好的程序進行上機計算。事先須將有限單元法的公式，計算方法和步驟都編入程序。4、對成果進行整理、分析。

對第1和第4步的工作，也盡可能讓計算機執(zhí)行，以減少人工的工作量。如自動劃分網(wǎng)格，整理成果等。94§6.8解題的具體步驟單元的劃分關(guān)于單元的劃分，注意幾點：（8）結(jié)構(gòu)具有凹槽或孔洞等應(yīng)力集中處等。（1）單元大小問題；（2）單元在不同部位的合理布置問題；（3）三角形三個內(nèi)角最好較接近；（4）利用對稱性和反對稱性；（5）厚度突變之處和材料不同之處；（6）載荷作用（集中力或突變分布載荷）處；（7）水利閘壩工程問題；95§6.9計算成果的整理

在有限單元法中，位移的精度較高，其誤差量級是，即與單元尺度的二次冪成正比。應(yīng)力的誤差量級是，即與單元的大小成正比。

96§6.9計算成果的整理三結(jié)點三角形單元的應(yīng)力的成果，不但應(yīng)力的精度較低，而且還產(chǎn)生了所謂應(yīng)力的波動性。對于結(jié)點位移的成果，可以直接采用。97§6.9計算成果的整理

應(yīng)力的波動性在三結(jié)點三角形單元中較為顯著。由于計算出的應(yīng)力的精度較低。假設(shè)Ⅰ單元的應(yīng)力成果為，其中為真解，為誤差。則由于在結(jié)點都列出了平衡方程并令其滿足，從而使相鄰的Ⅱ單元的應(yīng)力趨近于。這就產(chǎn)生了應(yīng)力的波動性。

原因是，98§6.9計算成果的整理為了提高應(yīng)力的精度，解決應(yīng)力波動性問題，可以采用兩種應(yīng)力成果的整理方法：一般地講，兩相鄰單元平均法的精度較好，因為它涉及的區(qū)域范圍較小。

（1）兩相鄰單元平均法。

（2）繞結(jié)點平均法。99§6.9計算成果的整理在受面力邊界線附近，求得的應(yīng)力誤差較大?？刹捎孟蛲獠逯档姆椒ǎɡ龗佄锞€插值）來解決。

100§6.9計算成果的整理為了提高應(yīng)力的精度，可以采用兩種方法。是加密網(wǎng)格，減少單元的尺寸，以提高應(yīng)力的精度。是可以采用較多結(jié)點的單元，并使位移模式中包含一些高冪次的項，從而提高位移和應(yīng)力的精度。二一101§6.10計算實例書中應(yīng)用三結(jié)點三角形單元，計算了下列例題：1.楔形體受自重及齊頂水壓力。2.簡支梁受均布荷載。3.圓孔附近的應(yīng)力集中。102§6.11等參單元誤差的產(chǎn)生：1.單元計算誤差；2.總剛度方程求解誤差。

克服誤差的途徑：1.加密網(wǎng)格；2.提高單元的計算精度。

高精度的單元模式需滿足的要求：1.單元形狀能較好地適應(yīng)復(fù)雜的幾何形狀（曲邊或曲面單元）；2.位移采用高階插值多項式，更加逼近真實位移；

但要構(gòu)造同時滿足以上要求的單元模式比較困難，怎么辦？103§6.11等參單元公認的做法是：

1.先在局部坐標(biāo)系中對規(guī)則形狀的單元（母單元）按高階插值多項式來構(gòu)造形狀函數(shù)，