有限元方法空間問題

有限元方法空間問題第一頁，共四十八頁，2022年，8月28日1空間問題——三維應(yīng)力狀態(tài)實際問題本質(zhì)上都是立體的、空間的。對承受載荷的彈性體，應(yīng)有三維應(yīng)力狀態(tài)。對彈性體內(nèi)每點的位移，有u、v、w分別代表對應(yīng)空間坐標系x、y、z方向的位移。u、v、w本身也代表彈性體內(nèi)的位移場，即它們都是物體內(nèi)有效的空間坐標的函數(shù)，一般可以表示為：u=u(x,y,z)，v=v(x,y,z)，w=w(x,y,z)。第二頁，共四十八頁，2022年，8月28日回顧空間問題的幾何方程為：第三頁，共四十八頁，2022年，8月28日按有限元的習(xí)慣寫法——算子形式，為：第四頁，共四十八頁，2022年，8月28日應(yīng)力矢量定義為：物理方程為：彈性矩陣[D]的一般形式為教材中(4-4)式。第五頁，共四十八頁，2022年，8月28日2簡單四面體單元2．1形狀函數(shù)一般的三維結(jié)構(gòu)，都可以劃分成很多小的四面體，為四面體單元。大量的小四面體單元拼合起來，可以逼近任意形狀的實際三維結(jié)構(gòu)體。簡單四面體單元如下圖，其中4個節(jié)點編號設(shè)為k、l、m、n。單元變形時，各節(jié)點都有沿x、y、z的3項位移，單元有4個節(jié)點，共有12項節(jié)點位移，合起來以列陣表示為：第六頁，共四十八頁，2022年，8月28日第七頁，共四十八頁，2022年，8月28日對于這種簡單的四面體單元，其內(nèi)部位移可假設(shè)為坐標的線性函數(shù)．為滿足完備性條件，應(yīng)取為第八頁，共四十八頁，2022年，8月28日上式含12個a參數(shù)，可以由單元的12項節(jié)點位移確定。將4個節(jié)點的坐標值代入(4．5a)中的u式，在k、l、m、n4個節(jié)點上，分別有第九頁，共四十八頁，2022年，8月28日由式(4．5b)求出al、a2、a3與a4，再代回式(4．5a)，整理后得：同理，用v式可求得a5到a8

，用w求得a9到a12,為：第十頁，共四十八頁，2022年，8月28日用矩陣記法統(tǒng)一表達為：[N]為形狀函數(shù)矩陣，可表示為：[I]為三階單位矩陣，而各節(jié)點的形狀函數(shù)可按下式計算得到，即第十一頁，共四十八頁，2022年，8月28日如記矩陣為四面體單元的體積，其他系數(shù)皆可由[L]確定，如第十二頁，共四十八頁，2022年，8月28日為矩陣第一行各元素的代數(shù)余子式。同樣可以確定al、bl、cl、dl…an、bn、cn、dn等，它們是矩陣[L]第二、三、四行元素的代數(shù)余子式。第十三頁，共四十八頁，2022年，8月28日在通常的右手坐標系xyz中，按上式計算時，四面體單元的4個節(jié)點排列的順序應(yīng)按右手規(guī)則，以使體積V為正。即由n點看klm平面，應(yīng)使k、l、n為逆時針排列。第十四頁，共四十八頁，2022年，8月28日簡單四面體單元內(nèi)，位移是坐標的線性函數(shù)，單元體的任一三角形界面，變形后仍保持為一平面，且由該面上3個節(jié)點的位移決定。因而相鄰兩單元的三角形交界面上，在變形過程中，其位移是一致的，即兩相鄰單元的位移在交界面上是連續(xù)的，單元滿足相容性條件。簡單四面體單元的形狀函數(shù)滿足完備性又滿足相容性要求，因而用此單元分析三維變形問題時，能收斂于精確解。第十五頁，共四十八頁，2022年，8月28日2．2單元剛陣將表達式(4．6)代入幾何關(guān)系式(4．2)，經(jīng)過微分運算，可以得到單元內(nèi)應(yīng)變?yōu)槠渲袘?yīng)變矩陣[B]是形狀函數(shù)矩陣經(jīng)微分算子矩陣作用所得的結(jié)果。[B]中任一個子矩陣[Bi]的顯式應(yīng)為：第十六頁，共四十八頁，2022年，8月28日由V及bi、ci、di等式可見前式，這里[Bi]的每項元素都是由節(jié)點坐標決定的常數(shù)。因而簡單四面體單元內(nèi)，各點的應(yīng)變都是一樣的，這是一種常應(yīng)變單元。第十七頁，共四十八頁，2022年，8月28日單元內(nèi)應(yīng)變?yōu)槌Ｖ?，按物理方程，單元?nèi)的應(yīng)力也是常值。當然，一般受力情況下，三維體內(nèi)有限大小的四面體內(nèi)的應(yīng)力并不是常值，用常應(yīng)力單元來代替它，只是近似的。對此單元，單元間的應(yīng)力是不連續(xù)的。只有當單元劃分得較小時，單元內(nèi)的應(yīng)力才會接近于常值，此時計算的應(yīng)力在單元間的不連續(xù)才會比較小，因而可以作為真實應(yīng)力分布的近似。一般，把這種單元應(yīng)力的計算值作為單元中心一點的應(yīng)力近似值是比較適當?shù)?。第十八頁，共四十八頁?022年，8月28日計算單元剛度矩陣的公式如前仍為：這里Ve為單元體積．由于簡單四面體單元為常應(yīng)變單元，故積分結(jié)果為：第十九頁，共四十八頁，2022年，8月28日按節(jié)點分塊，此單元剛陣可以表示為：其中任一個子矩陣為：第二十頁，共四十八頁，2022年，8月28日2．3載荷分配三維彈性體內(nèi)如受有均布的體積力(如重力)作用，對于這種簡單的四面體單元，可以逐個單元計算出整個單元的全部體積力，再平均分配到4個節(jié)點上，即每個節(jié)點分配1／4的單元體積力。如果單元的某個表面作用有均布的面積力(如氣體壓力)，也可將此面上的全部面積力平均分配到相應(yīng)的3個節(jié)點上，即每個節(jié)點分配到三角面上面積力總和的1／3。如果體積力、面積力不是均布的，則不應(yīng)平均分配，而應(yīng)按做功相等的原則等效分配。第二十一頁，共四十八頁，2022年，8月28日注意前泛函P計算公式中關(guān)于外力功的表達，有：為e單元內(nèi)分布體積力和分布面積力分配到單元節(jié)點的載荷，{q}和{p}分別為單位體積力和單位面積力。第二十二頁，共四十八頁，2022年，8月28日在討論了單元剛度矩陣及外載荷向節(jié)點的移置的計算公式后，由單元剛度陣去形成整體剛度矩陣及整體節(jié)點載荷向量的形成過程與前平面問題相比無邏輯上的差別。只要在“對號入座”的過程中注意在此每單元有4個節(jié)點，每節(jié)點有3個自由度即可。關(guān)于邊界條件、約束情況的處理也與前完全類似。因此，下面再討論某種有限元的時候，我們就不再討論“單元”到“整體”的形成過程了。我們知道：只要在單元水平講清該單元的節(jié)點數(shù)、節(jié)點自由度、形狀函數(shù)及單元剛度陣等單元的特征就已經(jīng)完全表達清楚了此單元的特點了。第二十三頁，共四十八頁，2022年，8月28日簡單四面體單元公式簡單，但是精度比較低，單元內(nèi)應(yīng)力為常值且單元間應(yīng)力不連續(xù)。為得到一定準確度的結(jié)果，往往要求將單元劃分比較小，增加了整個問題求解的自由度，總的計算效益是不理想的。但是，能很好的逼近任意幾何形狀是它的突出優(yōu)點，是目前大型CAD軟件所附的結(jié)構(gòu)分析模塊的網(wǎng)格自動劃分功能的常用單元。對三維問題的有限元分析，一般多采用復(fù)雜一些的、精度高一些的單元(如后的三維等參單元)，其綜合效益會更好。第二十四頁，共四十八頁，2022年，8月28日3軸對稱問題前已講過軸對稱問題。其結(jié)構(gòu)幾何特征是旋轉(zhuǎn)體，即幾何形狀對稱于中心軸。如果旋轉(zhuǎn)體所受的載荷也對稱于中心軸，則其變形也是對稱于此軸的。工程中常見的旋轉(zhuǎn)軸、輪盤、受均勻壓力的旋轉(zhuǎn)體容器等，都屬于軸對稱問題。第二十五頁，共四十八頁，2022年，8月28日如圖4·2，取柱坐標作為參考系。結(jié)構(gòu)受載荷而產(chǎn)生軸對稱變形時，其位移、應(yīng)變、應(yīng)力都與角坐標q無關(guān)，而只是徑向坐標r與軸向坐標z的函數(shù)。陰影部分為通過中心軸的平截面(子午面)。軸對稱變形的每個子午面的變形在柱坐標系內(nèi)是完全一樣的。因而，結(jié)構(gòu)雖處于三維應(yīng)力狀態(tài)，但可以只研究其任一個子午面內(nèi)的情況。第二十六頁，共四十八頁，2022年，8月28日用位移法，就是只研究這個代表截面的位移求得一個截面的位移分布，也就有了整個三維結(jié)構(gòu)內(nèi)的位移分布，從而可以求得體內(nèi)任一點的應(yīng)變及應(yīng)力。這樣，一個三維問題，就可以轉(zhuǎn)化為一個二維問題。

由于結(jié)構(gòu)的變形是對稱于中心軸的，因而子午面內(nèi)各點都只有沿徑向r的位移u和沿軸向z的位移w，一般應(yīng)為截面坐標r，z的函數(shù)，即第二十七頁，共四十八頁，2022年，8月28日軸對稱問題中，上述截面內(nèi)任一點p，實際上代表一個半徑為r的圓周(圖4-2)，當此圓周上各點都有徑向位移u時，圓周被拉伸，多出一個環(huán)向應(yīng)變eq。有：第二十八頁，共四十八頁，2022年，8月28日全部應(yīng)變的4項分量與兩項位移分量之間的幾何關(guān)系（幾何方程），以矩陣表示為:第二十九頁，共四十八頁，2022年，8月28日軸對稱問題的4項應(yīng)力分量，以列陣表示為：軸對稱問題的應(yīng)力與應(yīng)變間的物理關(guān)系仍寫為：其物理關(guān)系矩陣[D]仍見教材p69.。第三十頁，共四十八頁，2022年，8月28日4軸對稱問題的簡單三角形單元4．1形狀函數(shù)軸對稱問題的分析，轉(zhuǎn)化為對其任一個子午面的分析，可將此截面剖分為許多三角形單元，可構(gòu)造與前平面問題類似的簡單三角形單元。單元有3個節(jié)點，每節(jié)點有沿r及z的兩項位移u及w。單元有6個自由度。單元節(jié)點位移可以列陣表示為：第三十一頁，共四十八頁，2022年，8月28日單元內(nèi)位移場由節(jié)點位移插值表示為：如假定位移為坐標的線性函數(shù)，形函數(shù)矩陣[N]與平面三角形單元的完全相同，只不過需將其中的坐標x改為r，y改為z，即：第三十二頁，共四十八頁，2022年，8月28日其中D為單元的三角形面積，其他系數(shù)為：由平面問題的分析可知，這種形狀函數(shù)是滿足單元收斂的充分必要條件，故有限元分析結(jié)果能收斂于真解。第三十三頁，共四十八頁，2022年，8月28日4．2應(yīng)變與應(yīng)力將假定的位移代入式(4．12)，得到單元內(nèi)應(yīng)變?yōu)椋簩?yīng)變矩陣[B]按節(jié)點分塊表示為：由(4．12)，得到應(yīng)變矩陣[B]中任一子矩陣[Bi]為：第三十四頁，共四十八頁，2022年，8月28日其中bi、ci及D如前，而按物理關(guān)系式，有應(yīng)力

注意軸對稱問題三角形單元的形函數(shù)雖與平面問題三角形單元相同，但其應(yīng)變、應(yīng)力則不相同的。這里不僅有環(huán)向應(yīng)變eq及環(huán)向應(yīng)力sq，而且單元內(nèi)應(yīng)、應(yīng)力并非常值，是r、z的函數(shù)。第三十五頁，共四十八頁，2022年，8月28日4．3單元剛陣由于一個三角形單元實際上代表一個環(huán)狀的單元體，計算單元剛陣時，Ve為圖4－3所示的環(huán)狀單元體積域。由于單元為一旋轉(zhuǎn)體，其體積元素應(yīng)為：故單元剛度陣：計算式為：第三十六頁，共四十八頁，2022年，8月28日4．4節(jié)點載荷軸對稱結(jié)構(gòu)常受有內(nèi)外表面的分布壓力以及結(jié)構(gòu)旋轉(zhuǎn)時由離心力而形成的體積分布力，這些分布載荷均應(yīng)在每個單元內(nèi)，按做功等效的原則，分配到相應(yīng)的節(jié)點上，形成單元節(jié)點載荷，再疊加成整個結(jié)構(gòu)的載荷{Q}。這與前方法是一致的。但這里表面壓力與離心力的分配計算都有其特殊性，現(xiàn)說明如下：第三十七頁，共四十八頁，2022年，8月28日1．表面壓力如圖，此面分布力分配到e單元各節(jié)點的單元節(jié)點載荷為：注意到受載荷作用面為旋轉(zhuǎn)面，其面積元素可寫為：第三十八頁，共四十八頁，2022年，8月28日則有：其中為截面中三角形單元上受面積力作用的一個邊。對圖4－4問題，由于形狀函數(shù)Nn在lm邊上的值為零，故作用在lm邊上的分布載荷分配在n節(jié)點上的載荷亦為零，故此單元節(jié)點載荷為：第三十九頁，共四十八頁，2022年，8月28日工程結(jié)構(gòu)常受有面分布的壓力．其方向與受力面垂直，當單元受力邊與中心軸z有任意夾角。(圖4—4)時，如巳知法向壓力值a，則徑向及軸向力可計算求得，有:對于圖4—5所示的單元，單元節(jié)點載荷為第四十頁，共四十八頁，2022年，8月28日其形狀函數(shù)在lm邊的表達式為：完成積分運算后得：可見在軸對稱情況下，單元lm面上的均布載荷p向l、m節(jié)點分配的載荷并不相等的。第四十一頁，共四十八頁，2022年，8月28日還應(yīng)注意，由于單元上的一個節(jié)點實際上代表一個節(jié)圓，因此，這里分配到單元節(jié)點上的各項載荷，都是沿節(jié)圓均勻分布力的總和，而不一定是合力。這兩個單元節(jié)點載荷的和應(yīng)等于作用在圓環(huán)面均布載荷的總和，即第四十二頁，共四十八頁，2022年，8月28日2．離心力如旋轉(zhuǎn)體繞中心軸z轉(zhuǎn)動的角速度為w，則體內(nèi)任一點處單位體積的離心力大小為，沿半徑方向而指向外，是一種體積分布力。按體積分布力移置公式，e單元的離心力分配到e單元各節(jié)點的單元節(jié)點載荷為：第四十三頁，共四十八頁，2022年，8月28日為子午截面上單元三角形面積域。將形狀函數(shù)矩陣[N]的分塊表達式代入上式，可得：可見離心力只分配為徑向的節(jié)點載荷，而軸向節(jié)點載荷項皆為零。第四十四頁，共四十八頁，2022年，8月28日4．5求積問題計算單元剛度矩陣以及單元節(jié)點載荷時，常要計算三角形域內(nèi)的面積積分，其形式為一般情況下，上式應(yīng)采