角形單元的有限元法.ppt

《角形單元的有限元法.ppt》由會員分享，可在線閱讀，更多相關(guān)《角形單元的有限元法.ppt（98頁珍藏版）》請在裝配圖網(wǎng)上搜索。

1、第五章 三角形單元的有限元法,5.1 基本思想,把整體結(jié)構(gòu)離散為有限個單元，研究單元的平衡和變形協(xié)調(diào)；再把這有限個離散單元集合還原成結(jié)構(gòu)，研究離散結(jié)構(gòu)的平衡和變形協(xié)調(diào)。劃分的單元大小和數(shù)目根據(jù)計算精度和計算機能力來確定。,彈性懸臂板剖分與集合,單元、節(jié)點需編號,（1-1）,2、單元內(nèi)任意點的體積力列陣qV,（1-2）,1、單元表面或邊界上任意點的表面力列陣qs,5.2 基本力學量矩陣表示,3、單元內(nèi)任意點的位移列陣f,（1-3）,4、單元內(nèi)任意點的應變列陣 ,（1-4）,5、單元內(nèi)任意點的應力列陣,（1-5）,6、幾何方程,,（1-6）,將上式代入式(1-4），,7、物理方程矩陣式,（1-7）

2、,式中 E、彈性模量、泊松比。,上式可簡寫為,,（1-8）,其中,,對于彈性力學的平面應力問題，物理方程的矩陣形式可表示為：,(1-9）,,矩陣D稱為彈性矩陣。對于平面應變問題，將式(1-9）中的E換為 ，換為 。,5.3 位移函數(shù)和形函數(shù),1、位移函數(shù)概念 由于有限元法采用能量原理進行單元分析，因而必須事先設定位移函數(shù)。 “位移函數(shù)”也稱 “位移模式”，是單元內(nèi)部位移變化的數(shù)學表達式，設為坐標的函數(shù)。 一般而論，位移函數(shù)選取會影響甚至嚴重影響計算結(jié)果的精度。在彈性力學中，恰當選取位移函數(shù)不是一件容易的事情；但在有限元中，當單元劃分得足夠小時，把位移函數(shù)設定為簡單的多項式就可以獲得相

3、當好的精確度。這正是有限單元法具有的重要優(yōu)勢之一。,不同類型結(jié)構(gòu)會有不同的位移函數(shù)。這里，仍以平面問題三角形單元（圖1-2）為例，說明設定位移函數(shù)的有關(guān)問題。,圖1-2是一個三節(jié)點三角形單元，其節(jié)點i、j、m按逆時針方向排列。每個節(jié)點位移在單元平面內(nèi)有兩個分量：,,（1-10）,一個三角形單元有3個節(jié)點（以 i、j、m為 序），共有6個節(jié)點位移分量。其單元位移或單元節(jié)點位移列陣為：,2、位移函數(shù)設定,本問題選位移函數(shù)（單元中任意一點的位移與節(jié)點位移的關(guān)系）為簡單多項式：,(1-12）,式中：a1、a2、、a6待定常數(shù)，由單元位移的 6個分量確定。a1、a4代表剛體位移，a2、 a3 、 a5

4、、 a6 代表單元中的常應變，而且，位移函數(shù)是連續(xù)函數(shù)。,(1-11）,,,選取位移函數(shù)應考慮的問題,,（1）位移函數(shù)的個數(shù) 等于單元中任意一點的位移分量個數(shù)。本單元中有u和v，與此相應，有2個位移函數(shù)；,（3）位移函數(shù)中待定常數(shù)個數(shù) 待定常數(shù)個數(shù)應等于單元節(jié)點自由度總數(shù)，以便用單元節(jié)點位移確定位移函數(shù)中的待定常數(shù)。本單元有6個節(jié)點自由度，兩個位移函數(shù)中共包含6個待定常數(shù)。,（2）位移函數(shù)是坐標的函數(shù) 本單元的坐標系為：x、y；,（4）位移函數(shù)中必須包含單元的剛體位移。,（5）位移函數(shù)中必須包含單元的常應變。,（6）位移函數(shù)在單元內(nèi)要連續(xù)。相鄰單元間要盡 量協(xié)調(diào)。,條件（4）、（

5、5）構(gòu)成單元的完備性準則。 條件（6）是單元的位移協(xié)調(diào)性條件。 理論和實踐都已證明，完備性準則是有限元解收斂于真實解的必要條件。單元的位移協(xié)調(diào)條件構(gòu)成有限元解收斂于真實解的充分條件。 容易證明，三角形三節(jié)點常應變單元滿足以上必要與充分條件。,,（7）位移函數(shù)的形式 一般選為完全多項式。為實現(xiàn)（4）（6）的要求，根據(jù)Pascal三角形由低階到高階按順序、對稱地選?。欢囗検降捻棓?shù)等于（或稍大于）單元節(jié)點自由度數(shù)。,例：平面應力矩形板被劃分為若干三角形單元。,對任一單元，如單元，取位移函數(shù)：,、、、單元的位移函數(shù)都是,可以看出： 位移函數(shù)在單元內(nèi)是連續(xù)的；,以、的邊界26為例,兩條

6、直線上有兩個點重合，此兩條直線必全重合。,位移函數(shù)在單元之間的邊界上也連續(xù)嗎？是。,3、形函數(shù),形函數(shù)是用單元節(jié)點位移分量來描述位移函數(shù)的插值函數(shù)。,(1-13）,（1）形函數(shù)確定,現(xiàn)在，通過單元節(jié)點位移確定位移函數(shù)中的待定常數(shù)a1、a2、、a6 。設節(jié)點i、j、m的坐標分別為（xi、yi）、（ xj、yj ）、（ xm、ym ），節(jié)點位移分別為（ui、vi）、 （uj、vj） 、 （um、vm）。將它們代入式（1-12），有,從式(1-13）左邊3個方程中解出待定系數(shù)a1、a2、a3為,,,,(1-14）,式中， A為三角形單元的面積，有,,(1-15）,特別指出：為使求得面積的值為正值

7、，本單元節(jié)點號的次序必須是逆時針轉(zhuǎn)向，如圖所示。至于將哪個節(jié)點作為起始節(jié)點i，則沒有關(guān)系。,將式(1-14）代入式(1-12）的第一式，整理后得,,同理,,(1-16）,式中,,,,,(1-16）,令,,,(1-18）,位移模式(1-16）可以簡寫為,,(1-19）,式(1-19）中的Ni、Nj、Nm是坐標的函數(shù)，反應了單元的位移形態(tài)，稱為單元位移函數(shù)的形函數(shù)。數(shù)學上它反應了節(jié)點位移對單元內(nèi)任一點位移的插值，又稱插值函數(shù)。,,(1-16）,用形函數(shù)把式(1-16）寫成矩陣，有,縮寫為,(1-20）,,形函數(shù)是有限單元法中的一個重要函數(shù)，它具有以下性質(zhì)：,N為形函數(shù)矩陣，寫成分塊形式：,,(1-

8、21）,其中子矩陣,(1-22）,I是22的單位矩陣。,（2）形函數(shù)性質(zhì),性質(zhì)1 形函數(shù)Ni在節(jié)點i上的值等于1，在其它節(jié)點 上的值等于0。對于本單元，有,（i、j、m）,性質(zhì)2 在單元中任一點，所有形函數(shù)之和等于1。對 于本單元，有,,,圖1-3,？？？,,,,,,圖1-4,也可利用行列式代數(shù)余子式與某行或列元素乘積的性質(zhì)（等于行列式值或0）證明。,性質(zhì)3 在三角形單元的邊界ij上任一點（x，y），有,證,圖1-5,（1）,性質(zhì)4 形函數(shù)在單元上的面積分和在邊界上的線積分公式為,,(1-23）,式中 為 邊的長度。,5.4 單元應變和應力,根據(jù)幾何方程(1-6）和位移函數(shù)(1

9、-16）可以求得單元應變。,1、單元應變,對位移函數(shù)（式(1-16））,,(1-24）,(1-16）,求導后代入式(1-6），得到應變和節(jié)點位移的關(guān)系式。,,,上式簡寫一般式：,,(1-25）,式中， B單元應變矩陣。,對本問題，維數(shù)為36。它的分塊形式為：,,子矩陣,,(1-26）,,由于 與x、y無關(guān)，都是常量，因此B矩陣也是常量。單元中任一點的應變分量是B矩陣與單元位移的乘積，因而也都是常量。因此，這種單元被稱為常應變單元。,,2、單元應力,將式(1-25）代入物理方程式(1-8），得 單元應力,,(1-27）,也可寫為,,(1-28）,其中：S稱為單元應力矩陣，并有,

10、(1-29）,這里，D是33矩陣，B是36矩陣，因此S也是36矩陣。它可寫為分塊形式,,,,(1-30）,將彈性矩陣（式(1-9）） 和應變矩陣（式(1-26））代入，得子矩陣Si,由式(1-29）,,,(1-31）,,,式(1-31）是平面應力的結(jié)果。對于平面應變問題，只要將上式中的E換成 ，換成 即得。,,(1-32）,由于同一單元中的D、B矩陣都是常數(shù)矩陣，所以S矩陣也是常數(shù)矩陣。也就是說，三角形三節(jié)點單元內(nèi)的應力分量也是常量。 當然，相鄰單元的bi、ci(i，j，m)一般不完全相同，因而具有不同的應力，這就造成在相鄰單元的公共邊上存在著應力突變現(xiàn)象。但是隨著網(wǎng)格的細分，這種突

11、變將會迅速減小，收斂于平衡被滿足。,5.5 單元平衡方程,1、 單元應變能,對于平面應力問題中的三角形單元，設單元厚度為h 。,將式(1-25）和(1-8）代入上式進行矩陣運算，并注意到彈性矩陣D的對稱性，有,,,應變能 U為,,,由于和T是常量，提到積分號外，上式可寫成,引入矩陣符號k，且有,(1-33a）,式(1-33a）是針對平面問題三角形單元推出的。注意到其中hdxdy的實質(zhì)是任意的微體積dv，于是得計算k的一般式。,(1-33）,式(1-33）不僅適合于平面問題三角形單元，也是計算各種類型單元k的一般式。,5.6節(jié)中將明確k的力學意義是單元剛度矩陣。式(1-33）便是計算單元剛度矩陣

12、的基本矩陣式。它適合于各種類型的單元。,單元應變能寫成,(1-34）,2、 單元外力勢能,單元受到的外力一般包括體積力、表面力和集中力。自重屬于體積力范疇。表面力指作用在單元表面的分布載荷，如風力、壓力，以及相鄰單元互相作用的內(nèi)力等。,(1-33）,（1） 體積力勢能,單位體積中的體積力如式(1-35）所示。,單元上體積力具有的勢能Vv為,注意到式(1-20）,有,（2） 表面力勢能,面積力雖然包括單元之間公共邊上互相作用的分布力，但它們屬于結(jié)構(gòu)內(nèi)力，成對出現(xiàn)，集合時互相抵消，在結(jié)構(gòu)整體分析時可以不加考慮，因此單元分析時也就不予考慮。,,現(xiàn)在，只考慮彈性體邊界上的表面力，它只在部分單元上形成表

13、面力（右下圖）。設邊界面上單位面積受到的表面力如下式：,l單元邊界長度 h單元厚度 A表面力作用面積,qs,qs沿厚度均勻分布，則單元表面力的勢能Vs為,（3） 集中力勢能,當結(jié)構(gòu)受到集中力時，通常在劃分單元網(wǎng)格時就把集中力的作用點設置為節(jié)點。于是單元集中力Pc的勢能Vc為,,（4）總勢能,,把(1-35）式中原括符內(nèi)的部分用列陣Fd代替，,綜合以上諸式，單元外力的總勢能V為,(1-35）,Fd具有和相同的行、列數(shù)。則,(1-36）,由單元的應變能U(1-34）和外力勢能V(1-36），可得單元的總勢能,(1-37）,將式(1-37）代入，,根據(jù)彈性力學最小勢能原理：結(jié)構(gòu)處于穩(wěn)定平衡的必要和充

14、分條件是總勢能有極小值。,3、單元平衡方程,于是有，,式(1-38）是從能量原理導出的單元平衡方程。這個方程表達了單元力與單元位移之間的關(guān)系。其中，F(xiàn)d和單元節(jié)點力F具有相同的意義。,(1-38）,即得單元平衡方程,5.6 單元剛度矩陣,平衡方程(1-38）中的矩陣k是單元力和單元位移關(guān)系間的系數(shù)矩陣，代表了單元的剛度特性，稱為單元剛度矩陣。單元剛度矩陣的體積為nj nj， nj 是單元位移總數(shù)。其一般計算公式為：,1、一般計算公式,它與單元應變矩陣B和彈性矩陣D有關(guān)。,,對于平面應力三角形單元，應變矩陣B是常數(shù)矩陣，同時彈性矩陣D也是常數(shù)矩陣，于是式(1-33）可以化簡為,式中A表示三角形單

15、元的面積。h是單元厚度。,2、平面問題三角形單元剛度矩陣,（1）平面應力三角形單元,將式(1-9）和(1-26）代入上式，,即得平面應力三角形單元剛度矩陣。寫成分塊形式，有,,(1-40）,式(1-40)中子矩陣krs為22矩陣，有,,,(1-41）,（2）平面應變?nèi)切螁卧?對于平面應變問題，須將上式中的E換為 ， 換為 ，于是有,,,其中，bi(j,m)、ci(j,m)是形函數(shù)式(1-16）中的系數(shù)(式2-17)。,(1-42）,平面問題的單元剛度矩陣k不隨單元（或坐標軸）的平行移動或作n角度（n為整數(shù)）的轉(zhuǎn)動而改變。 由公式(1-41）、(1-42）知，krs矩陣和其中的br

16、、cr 、 bs、cs （r、s=i、j、m）有關(guān)。 單元平移時， bi、ci不變。,(3) 三角形單元剛度矩陣與坐標系無關(guān), 單元轉(zhuǎn)動時， bi、ci不變。 當單元旋轉(zhuǎn)時，各節(jié)點的編號保持不變。如圖1-7所示，圖a所示的單元旋轉(zhuǎn)時，到達圖b所示位置。,可以證明，這兩種情形的k是相同的。 其實，推演公式(1-40）、(1-41）、(1-42）時并沒有規(guī)定坐標系的方位，當坐標系旋轉(zhuǎn)任意角度時，也不影響剛度矩陣的結(jié)果。因此，平面問題的單元剛度矩陣可以認為是結(jié)構(gòu)坐標系中的單元剛度矩陣，沒有坐標變換問題。,(1-38）,（1）單元剛度矩陣中每個元素有明確的物理意義例如，kij表示單元第j個

17、自由度產(chǎn)生單位位移(j=1)，其他自由度固定（=0）時，在第i個自由度產(chǎn)生的節(jié)點力Fi。,,,,,主對角線上元素kii(i=1,nj)恒為正值。,,3、單元剛度矩陣性質(zhì),（2）k的每一行或每一列元素之和為零,以上式中第i行為例，,當所有節(jié)點沿x向或y向 都產(chǎn)生單位位移時，,單元作平動運動，無應變，也無應力。則有：,即：k的每一行元素之和為零。根據(jù)對稱性，每一列元素之和也為零。,,（3）k是對稱矩陣 由k各元素的表達式，可知k具有對稱性。,njnj,,對于主對角線元素對稱。對稱表達式：,kij = kji,證明, kij表示當單元位移中第j個元素為1(j=1)其余元素為零時，引起的單元力中的

18、第i個節(jié)點力Fi, kji表示當單元位移中第i個元素為1(i=1)其余元素為零時，引起的單元力中的第j個節(jié)點力Fj,由虛功原理，得,kij = kji,（4）單元剛度矩陣是奇異矩陣 即k的行列式為零（由行列式性質(zhì)） 。 單元剛度矩陣是在單元處于平衡狀態(tài)的前提下得出的。單元作為分離體看待，作用在它上面的外力（單元力）必定是平衡力系。然而，研究單元平衡時沒有引入約束。承受平衡力系作用的無約束單元，其變形是確定的，但位移不是確定的。所以出現(xiàn)性質(zhì)（3）中的“平動問題”，即單元可以發(fā)生任意的剛體運動。從數(shù)學上講，方程(1-28）的解不是唯一的或不能確定的。由此，單元剛度矩陣一定是奇異的。,（5）

19、單元剛度矩陣是常量矩陣,單元力和單元位移成線性關(guān)系是基于彈性理論的結(jié)果。,4、例：平面應力直角三角形單元剛度矩陣,圖1-8示出一平面應力直角三角形單元，直角邊長分別為a、b，厚度為h，彈性模量為E，泊松比為，計算單元剛度矩陣。,圖1-8,第一步：計算bi、ci和單元 面積A。,圖1-8,（1-17）,表2-1 單元節(jié)點坐標和bi、ci值（i、j、m）,參數(shù),節(jié)點,單元面積: A=ab/2, 計算步驟,第二步：求子矩陣 由式(1-41），算得,,,,其他從略。,第三步：形成k 將kii等按式(1-40）組集成k 。,,(1-43a）,2i-1 2i 2j-1 2j

20、2m-1 2m,紅色號碼是單元位移（1、 2、）在結(jié)構(gòu)中對應的節(jié)點位移的序號。,,,,,i,j,m,i,j,m,i、j、m表示單元中3個節(jié)點在結(jié)構(gòu)系統(tǒng)中的編號。,當a=b時，即等腰直角三角形單元，有,,(1-43b）,i j m,,,,,i,j,m,5.7 等價節(jié)點力,從前面單元分析可以看出：單元平衡所用到的的量均要屬于節(jié)點的量，如單元位移、單元力。載荷亦應如此，必須將體積力、表面力轉(zhuǎn)化到節(jié)點上去，成為等價節(jié)點力（載荷）。在第2.5節(jié)中已經(jīng)得到了公式(1-35）和(1-36） 。,這里，F(xiàn)d就是體積力、表面力和集中力之和的總等價節(jié)點力。,,,(1-44）,把總等價節(jié)

21、點力 Fd 分解成體積力、表面力和集中力的等價節(jié)點力之和，有,FV單元上體積力的等價節(jié)點力 FS單元上表面力的等價節(jié)點力 pC單元上節(jié)點上的集中力,注意到式(1-35），得體積力等價節(jié)點力計算公式：,表面力的等價節(jié)點力計算公式：,(1-45）,,(1-46）,1、體積力的等價節(jié)點力,2、表面力的等價節(jié)點力,,,3、等價節(jié)點力計算舉例,（1）單元自重,圖1-9所示平面應力三角形單元，單元厚度為h。單元單位體積自重為，自重指向y軸的負方向。,(1-45）,, 計算式,注意到形函數(shù)的性質(zhì)4：,(1-23）,得自重荷載的等價節(jié)點力,,根據(jù)體積力和式(1-45）、(1-21）、(1-22），得,(1

22、-47）,上式表明：自重載荷的等價節(jié)點力為單元重量的1/3。,（2）均布面力,單元邊界上作用了均勻的分布力，如圖1-10所示，其集度為qs。,(1-46）,(1-21）,根據(jù)式(1-46）、(1-21）和(1-22）, 計算式,,注意到形函數(shù)性質(zhì)4 ：,(1-23）,得,,(1-48）,,(1-22）,均勻分布力的等價節(jié)點力為,式(1-48）表明：在ij邊上受均布面力的平面問題三角形單元，其等價節(jié)點力等于將均布面力合力之半簡單地簡化到i、j節(jié)點上，方向與分布力方向相同。m節(jié)點上為零。,(1-48）,,（3）線性分布面力,,表面力集度在i點為qsx qsyT，而在j點為0。設坐標軸s的原點取在j

23、點，沿ji為正向， 。,ij邊上任一點的面力集度qs,,,,在ij邊上有：,將qs和上式代入式(1-46），有,由形函數(shù)的性質(zhì)3：,,,,(1-49）,式(1-49）表明：ij邊受線性分布面力： i點為qsx, qsyT，j點為0 時，其等價節(jié)點力可將總載荷的2/3分配給i點，1/3分配給j點，m點為零得出。,體積力和表面力向節(jié)點的移置符合靜力等效原理的前提條件是：線性位移模式。,5.7 系統(tǒng)分析,5.7.1 坐標系,研究各離散單元集合成整體結(jié)構(gòu)，集合整體結(jié)構(gòu)的平衡和變形協(xié)調(diào)，建立整體結(jié)構(gòu)平衡方程。,單元分析時采用的坐標系成為局部坐標或單元坐標（單元剛度矩陣的通用性）。而結(jié)構(gòu)系

24、統(tǒng)分析時，必須在統(tǒng)一的坐標系內(nèi)進行（各力學量才能疊加），稱為“結(jié)構(gòu)坐標”或“整體坐標”，如圖1-13所示。,,單元坐標系下，單元位移、單元力、單元剛度矩陣表示為：,整體坐標系下，單元位移、單元力、單元剛度矩陣表示為：,如何從單元坐標轉(zhuǎn)化為結(jié)構(gòu)坐標將在第4章中討論。,5.7.2 整體剛度矩陣,假設整體結(jié)構(gòu)被劃分為ne個單元和n個節(jié)點，在整體坐標系下，對于每個單元均有：,將上述這些方程集合起來（整體坐標下疊加），便可得到整個結(jié)構(gòu)的平衡方程。為此，需要將k、、F體積膨脹，分別擴大為n1n1、n11和n11的矩陣才能相加。膨脹后，原有節(jié)點號對應位置的元素不變，而其它元素均為零。,組裝方法：建立一個體積

25、為n1n1的方陣，按單元序號依次把結(jié)構(gòu)坐標單元剛度矩陣的元素放入該方陣中。 放入方法：（1）按單元節(jié)點編碼對號入座； （2）同位置元素累加。,式中：K為整體剛度矩陣，為整體節(jié)點位移列陣；P為整體等價節(jié)點荷載列陣。如下：,（1-50）,i,j,m,i,j,m,,,,,,,,,,例：平面三角單元,雙行雙列,,5.7.3 結(jié)構(gòu)剛度矩陣特性,1、結(jié)構(gòu)剛度矩陣元素的力學意義,把方程(1-50）寫開，,,=1,,(1-51）,2、結(jié)構(gòu)剛度矩陣是對稱矩陣 已知單元剛度矩陣是對稱矩陣（5.7節(jié)），用單元剛度矩陣組集 結(jié)構(gòu)剛度矩陣的 過程又沒有破壞 其對稱性，結(jié)構(gòu) 剛度矩陣必然也 是對稱的。當然，

26、 對稱性也可以通 過虛功原理得到證明。,結(jié)構(gòu)剛度矩陣中的任一元素kij是j為單位位移（ j =1），其它位移為零時的Pi。,,3、結(jié)構(gòu)剛度矩陣主對角線上的元素恒為正值 由性質(zhì)(1)可知，任一主對角線上元素kii是使節(jié)點位移i為一單位位移，其它節(jié)點位移為零時必須在第i號位移方向施加的力Pi。它的方向自然應與位移方向相同，因而是正值。,,4、結(jié)構(gòu)剛度矩陣是一個稀疏矩陣,,稀疏矩陣指：存在大量零元素。非零元素稀疏排列。,矩陣的每一列都有很多零元素?？疾炀仃囍械趈列。,再分析圖（1-14）。設節(jié)點b發(fā)生單位位移j=1，其它位移為零時， j只能在與點節(jié)b有直接聯(lián)系的 q 、 r節(jié)點引起節(jié)點力，不能

27、在其它節(jié)點引起節(jié)點 力。所以式（1-52） 中，只有和q、p、r、 b節(jié)點位移的相關(guān)元素 才不為零，其余的元素 都是零元素。,任一元素kij是j=1（其它=0）引起的Pi（i=1、2）,（1-52）,b,其它各列的情況也是類似的。 結(jié)構(gòu)的節(jié)點總數(shù)通常都比直接環(huán)繞于任何一個節(jié)點的節(jié)點數(shù)大得多，因而，結(jié)構(gòu)剛度矩陣中很大一部分元素是零，即所謂的稀疏矩陣。,5、結(jié)構(gòu)剛度矩陣是一個奇異矩陣,從單元剛度矩陣的奇異性討論中知，處于靜力平衡狀態(tài)的無約束單元可以發(fā)生任意的剛體位移。與單元剛度矩陣是奇異矩陣的理由一樣，無約束結(jié)構(gòu)的結(jié)構(gòu)剛度矩陣K也是奇異矩陣，即K的行列式為零。,5.7.6 引入支承約束的結(jié)構(gòu)節(jié)

28、點平衡方程,6、結(jié)構(gòu)剛度矩陣是常量矩陣,結(jié)構(gòu)剛度矩陣是常量矩陣。結(jié)構(gòu)的節(jié)點力和節(jié)點位移成線性關(guān)系都是基于彈性理論的結(jié)果。,（1-53）,用平衡方程（1-53）是解不出結(jié)構(gòu)的節(jié)點位移的，因為結(jié)構(gòu)剛度矩陣是奇異矩陣。因此，必須引入約束，排除任何剛體位移，使結(jié)構(gòu)為幾何不變體系。,方程（1-53）中的剛度矩陣K和節(jié)點荷載向量列陣P可分割為約束和自由兩部分：,式中，Pr是支承反力，約束位移,自由,約束,(1-55）,(1-56）,,展開（154），有：,Kff引入約束后的結(jié)構(gòu)剛度矩陣。它通對K引入約束后獲得，具體方法: 從無約束的結(jié)構(gòu)剛度矩陣K中刪去與受約束位移號對應的行和列，再將矩陣壓縮排列成nn階方

29、陣，即為約化后的結(jié)構(gòu)剛度矩陣Kff 。 Kff這是一個非奇異矩陣，它存在逆矩陣。,方程(1-55）是引入約束后的結(jié)構(gòu)節(jié)點平衡方程，用于計算結(jié)構(gòu)所有非剛性約束節(jié)點的節(jié)點位移。而方程(1-60）可以用來計算結(jié)構(gòu)所有受剛性約束節(jié)點的反力。,,,(1-61）,由式(1-55）即可解出全部未知的節(jié)點位移：,5.7.7 節(jié)點位移和單元力的解答,1、結(jié)構(gòu)節(jié)點位移,2、支座反力,把解出的f代入(1-56），即得支座反力Pr：,關(guān)于方程 (1-61）的解算方法，當Kff采用本章中上述方法組集時，可直接采用結(jié)構(gòu)力學中的高斯（Gauss）法求解。,(1-56）,至此，我們可以看出：系統(tǒng)分析的主要任務是： （

30、1）組集引入約束后的結(jié)構(gòu)剛度矩陣Kff； （2）求解式(1-55）給出的線性代數(shù)方程組。算出全部未知的節(jié)點位移。,至于支座反力的計算，實際計算時，根本不去組集式(1-56）中的矩陣Krf，不用式(1-62）而是直接對支承點使用節(jié)點平衡方程計算。,3、單元力,(1-62）,（1）全解公式,全解是指結(jié)構(gòu)在未經(jīng)簡化的實際荷載作用下（全解系統(tǒng)）的解答。,余解是結(jié)構(gòu)在節(jié)點荷載作用下（余解系統(tǒng)）的解答。節(jié)點荷載包括原來給定的節(jié)點荷載和由分布荷載簡化而來的等價節(jié)點荷載。,特解是單元從實際分布荷載到等價節(jié)點力過程中（特解系統(tǒng)）導致的單元節(jié)點力。,通過式(1-55）解出獲得全部未知節(jié)點位移，據(jù)此找出任何一個

31、單元的單元位移（單元節(jié)點位移），然后可以求得單元力。,圖(1-15）,=,+,（b）特解系統(tǒng),（c）余解系統(tǒng),（2）余解,,（3）特解,把單元等價節(jié)點力反向作用在單元的相應節(jié)點，即為單元力的特解。,4、結(jié)論,（1）全解的節(jié)點位移等于余解的節(jié)點位移 因為，特解給出的節(jié)點位移總是為零。,（2）對不承受分布荷載的單元，內(nèi)力由余解直接給出。,（3）對承受分布荷載作用的單元，其內(nèi)力由余解和特解的迭加給出。,F全e = F余e Fqe (1-64）,等價節(jié)點力,,3.9 舉例,如圖1-17a所示兩端固支的矩形深梁，跨度為2a，梁高為a，截面寬度為h，已知， =0、E，承受均布壓力q。試用有限元

32、法解此平面應力問題。,,,,圖1-17,(a),利用對稱性，可取梁的一半分析，共分2個單元，4個結(jié)點，支座約束如(a)圖。,解 1、劃分單元。,,（C)自由度,,,等腰直角三角形單元dde k,,2、單元剛度矩陣,由書公式（1-46b）,3 4 5 6 1 2,,,,,,1,3,,,,,5 6 3 4 7 8,,3,3,,,3、結(jié)構(gòu)剛度矩陣,引入約束前，結(jié)構(gòu)剛度矩陣是88階的。受約束的自由度號為1，3，4，5，7，8。刪除相應行和列后，結(jié)構(gòu)剛度矩陣剩下2行2列，相應的自由度號是：2，6。,3/2,-1,1,-1,1/2,+,4、結(jié)構(gòu)節(jié)點荷載向量,一般情況

33、節(jié)點荷載向量有8個元素。引入約束后剩下2個，相應序號為：2，6。,5、結(jié)構(gòu)無約束節(jié)點平衡方程,解得,根據(jù)結(jié)點位移找單元位移 ，代入（1-38）式可求各單元的單元力F 。,（1-38）,利用單元力F計算支點反力。,將單元位移代入式（1-28），其中S可由式（1-30）、（1-31）求得，從而可求出單元應力。,（1-28）,習題：算出式（1-28）的結(jié)果.,有限元分析之大腕版 一定要選最變態(tài)的題目，什么材料非線性啊，幾何非線性啊，接觸非線性啊，多物理場耦合啊，都給他弄進去。是個模型就幾百萬個單元，上千萬個節(jié)點，畫個剖面圖就要十幾個小時。再整一并行機群，TOP500的，張口就是 High Performance Computation，一口地道的倫敦腔，倍兒有面子。題目一扔進去就跑個把月，你要是一個星期以內(nèi)出結(jié)果，你都不好意思和別人打招呼。你說這樣一趟算下來要發(fā)多少Paper? 10篇？10篇？就1篇！你還別嫌少，說不定人家還發(fā)在會議上。你得琢磨牛人的心理啊，有能耐算這樣題目的人，根本就不在乎多發(fā)一篇兩篇文章。什么叫大牛知道么？大牛就是不求灌水，但求經(jīng)典。,