有限元方法空間問題

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

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

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

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