有限元理論與方法-第7講

1、講授內(nèi)容第7講（屬7周）2.應變血車確定了單元位移后，可以很方便地利用幾何方程和物理方程求得單元的應變和應力。作為平面問 題，單元內(nèi)具' 3個應變分量 X、y yy （各符號的意義見附錄 1）,用矩陣表示為XXv£ y-y xy u v將（2.1.4）式代入上式中，得到ebi 0 bj 0 bm 01 一一一-e0 G 0 cj 0 cms2Aci bi c j b j cm bm£ B Se式中B稱為應變矩陣，寫為分塊形式，即B=Bi Bj Bm而其子陣為bi 01Bi 0 ci（i,j,m）2A ci bi(2-1-7)(2-1-8)(2-1-9)（也即應力梯度

2、較大）的0又稱為初應變3節(jié)點三角形單元的 B是常量陣，所以稱為常應變單元。在應變梯度較大 部位，單元劃分應適當密集，否則將不能反映應變的真實變化而導致較大的誤差。上述應變中包括與應力有關(guān)的應變和與應力無關(guān)的應變兩部分，無關(guān)的應變X0龜y0xy0G由溫度變化、收縮、晶體生長等因素引起，對工程結(jié)構(gòu)一般只考慮溫度應變，無論線性和非線性溫 度，計算時可近似地采用平均溫度T Ti Tj Tm 3Tref3式中，Ti、Tj、Tm分別為節(jié)點i、j、m的溫度，Tref為參考溫度。對于平面應力問題，溫度 T引起的初始應變?yōu)槠渲校盀榫€膨脹系數(shù)。由于溫度變化在各向同性介質(zhì)中不引起剪切變形, 都指各向同性介質(zhì)。I對

3、平面應力問題，溫度 T引起的初始應變?yōu)殡?1當不考慮溫度的影響時，當前溫度即為參考溫度。響。3.單元應力根據(jù)物理方程，對平面應力問題，取應變分量xy由上式解出xyE2(1式中，D為彈性矩陣,所以對。=0。以后所述問題，除非特別說明,以后所述問題，除非特別說明，不考慮溫度影yE2(1Exyxy(Txyxy(2-1-10)(2-1-11)取決于彈性常數(shù)E和因?qū)⑹剑?-1- 7）代入式（2-1-10）得(T(2-1-12)SiSj Sm(2-1-13)bibiCiCi(i, j,m)Ci12bi式中，S為應力矩陣，反映了單元應力與單元節(jié)點位移之間的關(guān)系。由于單元應力和應變分量為常量, 所以單元邊界上

4、有應力階越，隨單元劃分變密，突變將減小。對平面應變問題，有四個應力分量:cy、衣y和(Zo取應變分量由應變分量解出以、(y> xy,彈性矩陣為1E1Ey2(1xy1Exy(2-1-14)E(1(1)(1 2 )2(1根據(jù)物理方程可以求解各應力分量。4.單元剛度矩陣單元節(jié)點力為Fe,節(jié)點虛位移為* e TS Fe(e*)e,節(jié)點虛應變?yōu)?0e,平面單元的厚度為to應用虛位移原理* Teddxdy* TA £ D 杷XdyB £代入上式整理得到FeaBtDB tdxdy 空可見單元剛度矩陣為對于三節(jié)點三角形單元,KeBTDBtdxdyA(2-1-15)面積為 A,所取為線

5、性位移模式，單元剛度矩陣為KeBTStA進一步表示為KiiKe KjiKmiKij Kjj KmjKim Kjm K mm對平面應力問題有K 日K rs 5412 AbrbsCrbs1F1"2CrCsbrCsbrCsCrCs1Crbs1 2r i, j,m; s i, j,m (2-1-16)-brbs2單元剛度矩陣表達單元抵抗變形的能力，其元素值為單位位移所引起的節(jié)點力，與普通彈簧的剛度系數(shù)具有同樣的物理本質(zhì)。例如子塊KijKrs12221121其中：上標1表示x方向自由度，2表示y方向自由度，后一上標代表單位位移的方向，前一上標代表單位位移引起的節(jié)點力方向。如k11表示j節(jié)點產(chǎn)生

6、單位水平位移時在i節(jié)點引起的水平節(jié)點力分量，k21表示j節(jié)點產(chǎn)生單位水平位移時在i節(jié)點引起的豎直節(jié)點力分量，其余類推。單元剛度矩陣為對稱矩陣。由于單元可有任意的剛體位移，給定的節(jié)點力不能惟一地確定節(jié)點位 移，可知單元剛度矩陣不可求逆，具有奇異性。5.等效節(jié)點載荷有限單元法分析只采用節(jié)點載荷，作用于單元上的非節(jié)點載荷都必須移置為等效節(jié)點載荷。可依 照靜力等效原則，即原載荷與等效節(jié)點載荷在虛位移上所作的虛功相等，求等效節(jié)點載荷。(1)集中力的移置。設單元ijm內(nèi)坐標為(x, y)的任意一點M受有集中載荷f=fx fyT,移置為等效 節(jié)點載荷Pe=Xi YiXjYjXmYmT。假想單元發(fā)生了虛位移，

7、其中，M點虛位移為u*=NT(%)e,其中(町為單元節(jié)點虛位移。按照靜力等效原則有(s*)e)TPe(2-1-17)Pe NT f(2)體力的移置。設單元承受有分布體力,單位體積的體力記為TPeaN qtdxdyq = qx qyT,其等效節(jié)點荷載為(2-1-18)(3)面力的移置。設在單元的某一個邊界上作用有分布的面力，單位面積上的面力為 在此邊界上取微面積 tds,對整個邊界面積分，得到P=Px pyT，PelN T Ptds(2-1-19)例2-1求單元在以下受力情況下的等效節(jié)點荷載: 力p、圖2-3所示jm邊受x方向線性分布力。y方向的重力為G、圖2-2所示ij邊受x方向均布圖2-2均

8、布力圖2-3線性分布力1求解說明利用上述公式求等效節(jié)點載荷，當原載荷是分布體力或面力時，進行積分運算是比較繁瑣的。但 在線性位移模式下，可以按照靜力學中力的分解原理直接求出等效節(jié)點載荷，上述三種情況等效節(jié)點 荷載分別為e GPe 0101013e 11Pe ptl - 0 - 0 0 022Pe-ptl- 0 0 2 0 1 02336 .整體分析結(jié)構(gòu)的整體分析就是將離散后的所有單元通過節(jié)點連接成原結(jié)構(gòu)物進行分析，分析過程是將所有 單元的單元剛度方程組集成總體剛度方程，引進邊界條件后求解整體節(jié)點位移向量?？傮w剛度方程實際上就是所有節(jié)點的平衡方程，由單元剛度方程組集總體剛度方程應滿足以下兩 個原

9、則：(1)各單元在公共節(jié)點上協(xié)調(diào)地彼此連接，即在公共節(jié)點處具有相同的位移。由于基本未知量為整體節(jié)點位移向量，這一點已經(jīng)得到滿足。(2)結(jié)構(gòu)的各節(jié)點離散出來后應滿足平衡條件， 也就是說，環(huán)繞某一節(jié)點的所有單元作用于該節(jié)點 的節(jié)點力之和應與該節(jié)點的節(jié)點載荷平衡。每一節(jié)點統(tǒng)一使用整體節(jié)點編號(如圖2-4所示)，第4單元節(jié)點編號i、j、m統(tǒng)一依次改為8、7、5。確定各單元的大域變換矩陣，如第G44單元為0 0 00 0 00 0 0其中，I為2X2階單位矩陣。圖2-4節(jié)點與單元編號求出各單元剛度矩陣，利用大域變換法求出Z構(gòu)整體剛度矩陣K,引入邊界條件，得到結(jié)構(gòu)的節(jié)點平衡方程K 6 P(2-1-19)進

10、而，求解節(jié)點位移、單元應力和應變。例2-2如圖2-5所示，一懸臂梁，自由端受合力為P的均布力作用，梁厚t=1,月1/3,求節(jié)點位移。P/2XP/2圖2-5結(jié)構(gòu)與離散求解說明結(jié)構(gòu)為平面應力問題，劃分為2個三角形單元、，有四個節(jié)點0）、（2, 0）、| （2, 1）、（0, 1）。單元、節(jié)點順序分別取3、1、2和對單兀：bi=0 , bj= 1 , bm=1 , Ci=2 , Cj=0 , Cm=一 2 o對單兀：bi=0 , bj=1 , bm= 1 , Ci=-2 , Cj=0 , Cm=2 o 單元的剛度矩陣為400240122021、2、3、4,坐標分別為（0, 1、3、4,剛度矩陣完全一

11、樣。21212 3E 023K K42322122174413利用大域變換法求出整體剛度矩陣為700133242212003E 232 014 1324427402120120003213214200212 00327421413節(jié)點荷載向量為位移向量為P=0 0 0 P/2 0 P/20 0 0T32200曲U1 V1 U2 V2 U3 V3 U4 V4T由結(jié)構(gòu)平衡方程求得節(jié)點位移為1.888.991.508.42 00T圖2-6六節(jié)點三角形單元7 .平面問題高次單元如前所述，三節(jié)點三角形單元因其位移模式是線性函數(shù)，應變與應力在單元內(nèi)都是常量，而彈性 體實際的應力場是隨坐標而變化的。因此，這種單元在各單元間邊界上應力有突變，存在一定誤差。為了更好地逼近實際的應變與應力狀態(tài)，提高單元本身的計算精度，可以增加單元節(jié)點而采用更高階 次的位移模式，稱為平面問題高次單元。如六節(jié)點三角形單元、矩形單元等。這里只介紹六節(jié)點三角 形單元與矩形單元的位移模式，其它單元的位移模式和具體求解步驟與三節(jié)點三角形單元類似，且應 用較少，不在贅述。圖2-7四節(jié)點矩形單元六節(jié)點三角形單元如圖2-6所示，在三角形單元各邊中點處增加一個節(jié)點，則每個單元有六個節(jié) 點，共有12個自由度。