空間實(shí)體單元演示文稿

空間實(shí)體單元演示文稿現(xiàn)在是1頁\一共有41頁\編輯于星期四（優(yōu)選）空間實(shí)體單元現(xiàn)在是2頁\一共有41頁\編輯于星期四(b)圖w為長方體單元，可以類似平面四節(jié)點(diǎn)矩形單元進(jìn)行分析。

(c)圖為任意八節(jié)點(diǎn)六面體單元，可以類似平面四節(jié)點(diǎn)任意四邊形等參元進(jìn)行分析。

(d)圖為20節(jié)點(diǎn)曲邊六面體單元，可以類似平面八節(jié)點(diǎn)曲邊四邊形等參元進(jìn)行分析。8.24節(jié)點(diǎn)四面體常應(yīng)變單元1、位移模式

如圖8-2所示，取四面體的4個(gè)頂點(diǎn)i,j,m,n

為節(jié)點(diǎn)。每一個(gè)節(jié)點(diǎn)有3個(gè)位移分量，即（8-1）現(xiàn)在是3頁\一共有41頁\編輯于星期四單元節(jié)點(diǎn)位移向量為（8-2）

與平面問題式（2-12）類似，假定單元內(nèi)一點(diǎn)的位移分量為坐標(biāo)的線性函數(shù)（8-3）將式（8-3）的第1式應(yīng)用于4個(gè)結(jié)點(diǎn)，則ijmnxyz圖8-2現(xiàn)在是4頁\一共有41頁\編輯于星期四（8-4）由此可解出a1～a4，再代回到式（8-3）的第1式，與式（2-19）的第1式類似，有（8-5）式中形函數(shù)具有與式（2-18）類似的形式：（8-6）其中現(xiàn)在是5頁\一共有41頁\編輯于星期四（8-7）（8-8）現(xiàn)在是6頁\一共有41頁\編輯于星期四在式（8-8）中，V為四面體的體積。為使其計(jì)算值不為負(fù)，單元的節(jié)點(diǎn)（i,j,m,n）編號(hào)次序應(yīng)遵循右手法則。(p4)

采用同樣的方法，可得（8-9）（8-10）將式（8-5）、（8-9）（8-10）統(tǒng)一用矩陣式表示，可得與平面問題式（2-20）類似的公式（8-11）現(xiàn)在是7頁\一共有41頁\編輯于星期四式中[N]為單元形函數(shù)矩陣，其維數(shù)為3×12。進(jìn)一步可寫為與平面問題式（2-21）、（2-22）類似的子塊形式（8-12）其中，子矩陣（8-13）式中，I為3階單位矩陣。2、應(yīng)變矩陣

現(xiàn)在是8頁\一共有41頁\編輯于星期四在空間問題中，每點(diǎn)有6個(gè)應(yīng)變分量。幾何方程為：（8-14）將式（8-11）～（8-13）和（8-6）代入上式，得（8-15）式中（8-16）現(xiàn)在是9頁\一共有41頁\編輯于星期四上式（8-15）、（8-16）與平面問題式（2-24）～（2-26）類似。與平面問題三節(jié)點(diǎn)三角形單元相同，在四節(jié)點(diǎn)四面體單元中，[B]的元素都是常量，因此是常應(yīng)變單元。2、應(yīng)力矩陣

三維問題的應(yīng)力應(yīng)變關(guān)系也可寫為式（2-8）的矩陣形式（8-17）

與平面問題不同，這里和分別由6個(gè)分量組成，彈性矩陣[D]是一個(gè)6×6的矩陣：（8-18）現(xiàn)在是10頁\一共有41頁\編輯于星期四（8-19）（8-20）現(xiàn)在是11頁\一共有41頁\編輯于星期四

將式（8-20）所表示的[D]和式（8-15）、（8-16）所表示的[B]代入式（8-22），并將[S]定成分塊矩陣的形式，有將式（8-15）代入式（8-17），得（8-21）應(yīng)力矩陣[S]為（8-22）

由于[D]、[B]都是常數(shù)矩陣，因此應(yīng)力矩陣[S]也是常數(shù)矩陣。也就是說，單元中的應(yīng)力分量也是常數(shù)。（8-23）現(xiàn)在是12頁\一共有41頁\編輯于星期四式中（8-24）其中（8-25）3、單元?jiǎng)偠染仃?/p>

仿照平面問題中的推導(dǎo)，可得單元平衡方程現(xiàn)在是13頁\一共有41頁\編輯于星期四（8-26）單元?jiǎng)偠染仃嚲哂信c式（2-33）類似的形式（8-27）式中，[k]是一個(gè)12×12的矩陣。由于[B]、[D]都是常數(shù)矩陣，所以[k]也是一個(gè)常量矩陣。并且（8-28）寫成分塊矩陣的形式，有（8-29）現(xiàn)在是14頁\一共有41頁\編輯于星期四式中子矩陣[krs]為3×3的矩陣（8-30）4、等價(jià)節(jié)點(diǎn)力向量

式（8-26）中的單元等價(jià)節(jié)點(diǎn)力也包括體積力、表面力、集中力幾部分。體積力與表面力的計(jì)算公式與平面三角形單元公式（2-36）、（2-37）類似：（8-31）現(xiàn)在是15頁\一共有41頁\編輯于星期四（8-32）對于簡單情形，也可采用靜力等效原則簡化計(jì)算。

進(jìn)一步的整體平衡方程的建立（即結(jié)構(gòu)剛度矩陣、結(jié)構(gòu)等價(jià)節(jié)點(diǎn)力列陣的組集）、位移約束條件的引入、線性方程組的求解等，和平面問題有限元法一樣，不再贅述。8.3二十結(jié)點(diǎn)六面體等參數(shù)單元

由于精度高，容易適應(yīng)不同邊界，在平面問題中常選用了八節(jié)點(diǎn)四邊形等參數(shù)單元。與此類似，在三維問題中,常選用二十節(jié)點(diǎn)六面體等參數(shù)單元。現(xiàn)在是16頁\一共有41頁\編輯于星期四如圖8-3所示，在整體坐標(biāo)系的二十節(jié)點(diǎn)六面體的實(shí)際單元與中心在局部坐標(biāo)的原點(diǎn)、邊長為2的立方體基本單元相對應(yīng)。1234567891011121314151617181920????????????????????xyz6121812345789101113141516171920????????????????????=-1=1=1=-1=1圖8-3實(shí)際單元基本單元1、位移模式、形函數(shù)和坐標(biāo)變換式=-1現(xiàn)在是17頁\一共有41頁\編輯于星期四（8-33）（8-34）形函數(shù)的表達(dá)式如下：

（8-35）式中位移模式和坐標(biāo)變換式可寫為如下形式：現(xiàn)在是18頁\一共有41頁\編輯于星期四根據(jù)幾何方程，單元中的應(yīng)變?yōu)?、應(yīng)變矩陣（8-36）其中（8-37）應(yīng)變矩陣的子矩陣[Bi]為：現(xiàn)在是19頁\一共有41頁\編輯于星期四（8-38）現(xiàn)在是20頁\一共有41頁\編輯于星期四根據(jù)復(fù)合函數(shù)求導(dǎo)規(guī)則，有（8-39）式中[J]-1為雅可比矩陣[J]的逆矩陣。[J]的表達(dá)式為（8-40）現(xiàn)在是21頁\一共有41頁\編輯于星期四3、應(yīng)力矩陣單元中的應(yīng)力為單元應(yīng)力矩陣[S]為（8-41）4、單元?jiǎng)偠染仃噯卧獎(jiǎng)偠染仃嚳蓪懗桑?-42）[k]是一個(gè)60×60的矩陣，式（8-42）通常采用高斯法進(jìn)行積分?，F(xiàn)在是22頁\一共有41頁\編輯于星期四5、等價(jià)節(jié)點(diǎn)力矩陣等價(jià)節(jié)點(diǎn)力計(jì)算公式如下：（1）體積力

設(shè)單位體積力是，則等價(jià)節(jié)點(diǎn)力為（8-43）（2）表面力現(xiàn)在是23頁\一共有41頁\編輯于星期四設(shè)某邊界面上作用表面力，則等價(jià)節(jié)點(diǎn)力為（8-44）設(shè)該邊界面對應(yīng)基本單元=±1的面。由數(shù)學(xué)公式，結(jié)構(gòu)坐標(biāo)下曲面微元dS對應(yīng)于單元坐標(biāo)下的微元面積式為（8-45）現(xiàn)在是24頁\一共有41頁\編輯于星期四則式（8-44）可寫為單元坐標(biāo)系下的積分公式（8-46）以上為=±1的表面力計(jì)算公式。對于其他表面力，可類似處理。式（8-46）通常也采用高斯法進(jìn)行積分計(jì)算。現(xiàn)在是25頁\一共有41頁\編輯于星期四8.4空間軸對稱單元

許多工程構(gòu)件，其幾何形狀、約束條件及所受的荷載都對稱于某一軸，因而所有的位移、應(yīng)變和應(yīng)力分量也都對稱于該軸。這類問題稱空間軸對稱問題。

對空間軸對稱問題，采用圓柱坐標(biāo)系，r表示徑向坐標(biāo)，z表示軸向坐標(biāo)，任一對稱面為rz面。在有限元分析時(shí)，采用繞對稱軸旋轉(zhuǎn)一周的軸對稱環(huán)形單元。把軸對稱的工程構(gòu)件與rz坐標(biāo)平面正交的截面劃分為一些三角形，例如其中一個(gè)是ijm（圖8-4）。由這些三角形旋轉(zhuǎn)一周就構(gòu)成上述環(huán)形單元?，F(xiàn)在是26頁\一共有41頁\編輯于星期四圖8-4軸對稱構(gòu)件及三角形環(huán)形單元r(u)z(w)ij現(xiàn)在是27頁\一共有41頁\編輯于星期四當(dāng)然，也可以劃分為矩形，構(gòu)成矩形截面的環(huán)形單元?？傊瑔卧孛嫘螤?，可采用三角形、四邊形等平面有限元法所用單元形狀。本節(jié)只討論三角形截面的環(huán)形單元，其它截面形狀環(huán)形單元可參照本節(jié)方法進(jìn)行分析。

各個(gè)環(huán)形單元（后簡稱單元）間用位于三角形頂點(diǎn)的環(huán)形鉸聯(lián)系起來，就成了離散結(jié)構(gòu)模型。我們把環(huán)形鉸在rz平面的交點(diǎn)（即三角形頂點(diǎn)）稱為單元的節(jié)點(diǎn)（如i,j,m

）。因此對軸對稱問題進(jìn)行有限元分析只需在rz平面劃分網(wǎng)格，就像平面問題在xy平面上劃分網(wǎng)格一樣。照此思路，軸對稱空間問題被大大簡化。現(xiàn)在是28頁\一共有41頁\編輯于星期四1、位移模式

對于如圖8-4所示軸對稱三角形環(huán)形單元，是由如圖8-5所示平面上的三角形繞對稱軸軸回旋一周得到的。考慮到問題的軸對稱特點(diǎn)，可以借助于平面問題三節(jié)點(diǎn)三角形單元位移模式進(jìn)行單元分析。圖8-5rzijmrirjrm

由于軸對稱，環(huán)向位移恒等于零。只有徑（r）向位移和軸（z）向位移，它們恰在rz坐標(biāo)平面上。設(shè)徑向位移為u，軸向位移為w。對于圖8-5所示情形，借助平面問題的三角形單元，取位移模式為？現(xiàn)在是29頁\一共有41頁\編輯于星期四（8-47）代入節(jié)點(diǎn)位移后，可解出a1～a6，再代入上式，得（8-48）其中，形函數(shù)（8-49）式中A,ai,bi,ci，與平面問題三角形單元的對應(yīng)公式（2-15）、（2-17）一致，區(qū)別僅僅是將那里的x,y換成這里的r,z。

式（8-48）也可寫為式（2-20）同樣的形式現(xiàn)在是30頁\一共有41頁\編輯于星期四（8-50）式中，[Ni]=NiI(i,j,m)，其中I為2階單位矩陣。2、應(yīng)變矩陣根據(jù)彈性力學(xué)理論，空間軸對稱問題的幾何方程為（8-51）現(xiàn)在是31頁\一共有41頁\編輯于星期四與平面問題中只含有3個(gè)應(yīng)變分量不同，這里含有4個(gè)應(yīng)變分量。將u,w的表達(dá)代入式（8-51），得（8-52）式中（8-53）其中（8-54）現(xiàn)在是32頁\一共有41頁\編輯于星期四由式（8-52）～（8-54）可見，矩陣中含有變量r,z，因此它不是常數(shù)矩陣。即軸對稱問題的三角形環(huán)形單元不是常應(yīng)變單元。3、應(yīng)力矩陣根據(jù)彈性力學(xué)理論，空間軸對稱問題的應(yīng)力-應(yīng)變關(guān)系為（8-55）式中[D]是軸對稱問題的彈性矩陣現(xiàn)在是33頁\一共有41頁\編輯于星期四（8-56）將式（8-52）代入式（8-55），得（8-57）應(yīng)力矩陣

[S]=[D][B]顯然，除rz外，單元中其他應(yīng)力分量不是常數(shù)?，F(xiàn)在是34頁\一共有41頁\編輯于星期四4、單元?jiǎng)偠染仃囕S對稱問題的單元?jiǎng)偠染仃嚳捎墒剑?-27）計(jì)算由于被積函數(shù)與無關(guān)，故在三角形截面的環(huán)單元的積分可簡化為在三角形截面上的積分：（8-58）

式（8-58）中，[B]不是常數(shù)，而由式（8-53）、（8-54）確定。式（8-53）、（8-54）中，ai,bi,ci是（P32）（1）一般公式（2）近似計(jì)算現(xiàn)在是35頁\一共有41頁\編輯于星期四可見，單元?jiǎng)偠染仃囉?jì)算比平面三角形單元麻煩。須對其進(jìn)行數(shù)值計(jì)算或近似計(jì)算?，F(xiàn)在，討論面積分?Ag(r,z)drdz的計(jì)算問題。雖為常數(shù)，由式（2-17）計(jì)算，但fi是r,z的函數(shù)。式（8-58）被積函數(shù)[B]T[D][B]r是r,z的函數(shù)，暫簡寫為g(r,z)。式（8-58）可寫為圖8-6rzijm對于一個(gè)典型三角形截面ijm（圖8-6），過j點(diǎn)作垂直于r軸的直線，將三角形ijm，分成兩個(gè)三角形：A1，A2。于是，有A2A1現(xiàn)在是36頁\一共有41頁\編輯于星期四（8-59）

完成上述積分可采用數(shù)值方法，有專門的二重?cái)?shù)值積分程序可供引用。

也可采用近似方法完成：當(dāng)單元較小時(shí)，常把各個(gè)單元中的r