第三章-桿梁結(jié)構(gòu)的有限元分析原理-課件

第3章桿系結(jié)構(gòu)的有限元法

2020/10/281第3章桿系結(jié)構(gòu)的有限元法

2020/10/281桿、梁?jiǎn)卧攀鲇懻摋U梁?jiǎn)卧陀伤鼈兘M成的平面和空間桿梁結(jié)構(gòu)系統(tǒng).從構(gòu)造上來(lái)說(shuō)其長(zhǎng)度遠(yuǎn)大于其截面尺寸的一維構(gòu)件承受軸力或扭矩的桿件成為桿桿梁?jiǎn)栴}都有精確解（且是唯一的）承受橫向力和彎矩的桿件稱為梁平面桁架平面剛架連續(xù)梁空間剛架空間桁架等承受軸力或扭矩的桿件稱為桿將承受橫向力和彎矩的桿件稱為梁變截面桿和彎曲桿件2020/10/282桿、梁?jiǎn)卧攀鲇懻摋U梁?jiǎn)卧陀伤鼈兘M成的平面和空間桿梁結(jié)構(gòu)系精品資料2020/10/283精品資料2020/10/283本章主要內(nèi)容2020/10/284本章主要內(nèi)容2020/10/2843.1有限元分析的完整過(guò)程E1=E2=2E7PaA1=A2=2cm2l1=l2=10cmP3為10N作用下二桿結(jié)構(gòu)的變形。E1、A1E2、A2說(shuō)明：u1、u2、u2分別表示節(jié)點(diǎn)1、2、3的水平位移2020/10/2853.1有限元分析的完整過(guò)程E1=E2=2E7PaA1=A2=1）用標(biāo)準(zhǔn)化的分段小單元來(lái)逼近原結(jié)構(gòu)2）尋找能夠滿足位移邊界條件的許可位移場(chǎng)3）基于位移場(chǎng)的最小勢(shì)能原理來(lái)求解基本變量為：節(jié)點(diǎn)位移內(nèi)部各點(diǎn)位移應(yīng)變應(yīng)力(1)(3)(2)問(wèn)題的解題思路2020/10/2861）用標(biāo)準(zhǔn)化的分段小單元來(lái)逼近原結(jié)構(gòu)節(jié)點(diǎn)位移內(nèi)部各點(diǎn)位移應(yīng)變完整的求解過(guò)程1）結(jié)構(gòu)離散化該構(gòu)件由兩根桿件做成，因此可以自然離散成2個(gè)桿單元。假定以這類單元位移的特征為兩個(gè)端點(diǎn)位移，就這兩個(gè)離散單元給出節(jié)點(diǎn)編號(hào)和單元編號(hào)。單元1：i=1，j=2單元2：i=2，j=32020/10/287完整的求解過(guò)程1）結(jié)構(gòu)離散化2020/10/287單元位移模式：u(x)=a0+a1x單元節(jié)點(diǎn)條件：u(0)=u1，u(l)=u2

將式（b）代入式（a），從而得2）單元分析（a）（b）2020/10/2882）單元分析（a）（b）2020/10/288回代得寫(xiě)成矩陣形式為其中Ni，Nj是形函數(shù)。形函數(shù)矩陣說(shuō)明：u表示位移列陣ue表示單元位移2020/10/289回代得根據(jù)幾何方程可得單元應(yīng)變的表達(dá)單元應(yīng)變寫(xiě)成矩陣形式為簡(jiǎn)記為幾何函數(shù)矩陣或者是應(yīng)變轉(zhuǎn)換矩陣2020/10/2810根據(jù)幾何方程可得單元應(yīng)變的表達(dá)幾何函數(shù)矩陣或者是應(yīng)變轉(zhuǎn)換矩陣根據(jù)物理方程可得單元應(yīng)力的表達(dá)單元應(yīng)力寫(xiě)成矩陣形式為簡(jiǎn)記為單元應(yīng)力矩陣或者是應(yīng)力轉(zhuǎn)換矩陣節(jié)點(diǎn)位移列陣2020/10/2811根據(jù)物理方程可得單元應(yīng)力的表達(dá)單元應(yīng)力矩陣或者是應(yīng)力轉(zhuǎn)換矩陣單元e勢(shì)能的表達(dá)說(shuō)明積分域，P1、P2、分別表示作用單元e上的節(jié)點(diǎn)在u1、u2的力2020/10/2812單元e勢(shì)能的表達(dá)說(shuō)明積分域，P1、P2、分別表示寫(xiě)成矩陣形式為單元e剛度矩陣單元e節(jié)點(diǎn)力列陣2020/10/2813寫(xiě)成矩陣形式為單元e剛度矩陣單元e節(jié)點(diǎn)力列陣2020/10/在得到各個(gè)單元的勢(shì)能表達(dá)式后，需要進(jìn)行離散單元的裝配，以求出整個(gè)系統(tǒng)的總勢(shì)能，對(duì)于該系統(tǒng)，總勢(shì)能包括兩個(gè)單元部分3）離散單元的裝配2020/10/2814在得到各個(gè)單元的勢(shì)能表達(dá)式后，需要進(jìn)行離散單元的裝配，以求出處理邊界條件是獲取可能位移場(chǎng)，將左端的約束條件，即u1=0代入上式可以得到簡(jiǎn)化的勢(shì)能表達(dá)式4）邊界條件的處理2020/10/2815處理邊界條件是獲取可能位移場(chǎng)，將左端的約束條件，即u1=0代由于上式是基于許可位移場(chǎng)的表達(dá)的系統(tǒng)勢(shì)能，這是由全部節(jié)點(diǎn)位移分段所插值出的位移場(chǎng)為全場(chǎng)許位移場(chǎng)，且基本未知量為節(jié)點(diǎn)位移，根據(jù)最小勢(shì)能原理（即針對(duì)未知位移求一階變分）有5）建立剛度方程2020/10/2816由于上式是基于許可位移場(chǎng)的表達(dá)的系統(tǒng)勢(shì)能，這是由全部節(jié)點(diǎn)位5將結(jié)構(gòu)參數(shù)和外載荷代入上式有求解得（單位m）6）求解節(jié)點(diǎn)位移2020/10/2817將結(jié)構(gòu)參數(shù)和外載荷代入上式有6）求解節(jié)點(diǎn)位移2020/10/7）計(jì)算單元應(yīng)變2020/10/28187）計(jì)算單元應(yīng)變2020/10/28188）計(jì)算單元應(yīng)力2020/10/28198）計(jì)算單元應(yīng)力2020/10/2819對(duì)于單元?jiǎng)菽艿谋磉_(dá)，對(duì)其取極值有具體地對(duì)于單元1，有其中R1是節(jié)點(diǎn)1的支反力，P2是單元1的節(jié)點(diǎn)2所受的力，即單元2對(duì)該節(jié)點(diǎn)的作用力，將前面求得的節(jié)點(diǎn)位移代入上式可得支反力大小。9）計(jì)算支反力2020/10/2820對(duì)于單元?jiǎng)菽艿谋磉_(dá)，對(duì)其取極值有9）計(jì)算支反力2020/10以上是一個(gè)簡(jiǎn)單結(jié)構(gòu)有限元方法求解得完整過(guò)程，對(duì)于復(fù)雜結(jié)構(gòu)，其求解過(guò)程完全相同，由于每一個(gè)步驟都具備標(biāo)準(zhǔn)化和規(guī)范性的特征，所以可以在計(jì)算機(jī)上編程而自動(dòng)實(shí)現(xiàn)。討論1：對(duì)于一個(gè)單元的勢(shì)能取極值，所得到的方程為節(jié)點(diǎn)的位移和節(jié)點(diǎn)力之間的關(guān)系，也稱為單元的平衡關(guān)系，由此可以求出每一個(gè)單元所受的節(jié)點(diǎn)力。2020/10/2821以上是一個(gè)簡(jiǎn)單結(jié)構(gòu)有限元方法求解得完整過(guò)程，對(duì)于復(fù)雜結(jié)構(gòu)，其討論2：由前面的步驟，我們也可以直接將各個(gè)單元的剛度矩陣按照節(jié)點(diǎn)編號(hào)的對(duì)應(yīng)位置來(lái)進(jìn)行裝配，即在未處理邊界條件之前，先形成整體剛度矩陣。其物理意義是，表示在未處理邊界條件前的基于節(jié)點(diǎn)描述的總體平衡關(guān)系。在對(duì)該方程進(jìn)行位移邊界條件的處理后就可以求解，這樣與先處理邊界條件再求系統(tǒng)勢(shì)能的最小值所獲得的方程完全相同。2020/10/2822討論2：由前面的步驟，我們也可以直接將各個(gè)單元的剛度矩陣按照

小結(jié)有限元分析的基本步驟及表達(dá)式1、物體幾何區(qū)域的離散化2、單元的研究（所有力學(xué)信息都用節(jié)點(diǎn)位移）來(lái)表達(dá)3、裝配集成4、邊界條件的處理并求解節(jié)點(diǎn)位移5、支反力的求取以及其它力學(xué)量（應(yīng)力、應(yīng)變及位移三大物理量）的計(jì)算2020/10/2823小結(jié)有限元分析的基本步驟及表達(dá)式1、物體幾有限元分析的基本步驟及表達(dá)式2020/10/2824有限元分析的基本步驟及表達(dá)式2020/10/2824一拉壓桿單元圖2.1拉壓桿單元示意圖

設(shè)桿單元長(zhǎng)度為，橫截面面積為，單元材料的彈性模量為，在局部坐標(biāo)系中桿端荷載分別為和，桿端位移分別為和，單元上的軸向分布荷載為。3.2局部坐標(biāo)下的桿單元分析2020/10/2825一拉壓桿單元圖2.1拉壓桿單元示意圖設(shè)桿單元長(zhǎng)度為用結(jié)點(diǎn)位移表示單元上任意截面的位移。對(duì)拉壓桿單元，可以取其位移為一次多項(xiàng)式，即

由位移的邊界條件：可得系數(shù)、為：這樣，截面任意一點(diǎn)的位移為：

用矩陣表示為：

其中

（3-1）

（3-2）

①單元位移模式。2020/10/2826用結(jié)點(diǎn)位移表示單元上任意截面的位移。對(duì)拉壓桿單元，可以取根據(jù)材料力學(xué)中應(yīng)變的定義，有

這里為應(yīng)變矩陣。由虎克定律，其應(yīng)力為：

（3-4）

（3-3）

②

進(jìn)行應(yīng)力、應(yīng)變分析其中2020/10/2827根據(jù)材料力學(xué)中應(yīng)變的定義，有（3-4）（利用虛位移原理求單元?jiǎng)偠染仃?，設(shè)桿端i、j分別產(chǎn)生虛位移、，則由此引起的桿軸截面任意位置的虛位移為：對(duì)應(yīng)的虛應(yīng)變?yōu)椋?/p>

根據(jù)虛位移原理虛功方程（力乘以虛位移得虛功、外力虛功等于變形虛功），有：

將上式整理得：

（3-5）

（3-6）

③

求單元?jiǎng)偠染仃?020/10/2828（3-5）（3-6）③求單元?jiǎng)偠染仃?020/式中：為局部坐標(biāo)系下單元結(jié)點(diǎn)荷載矩陣。

記

則可以得到拉壓桿單元的單元?jiǎng)偠确匠虨椋?/p>

這里為局部坐標(biāo)系下的單元?jiǎng)偠染仃嚕瑸榫植孔鴺?biāo)系下等效結(jié)點(diǎn)荷載矩陣，但值得指出的是：分布荷載中可以包含集中荷載。根據(jù)定義，可以進(jìn)一步求得單元?jiǎng)偠染仃嚍椋?/p>

（3-10）

（3-7）

（3-8）

（3-9）

等效結(jié)點(diǎn)荷載2020/10/2829式中二扭轉(zhuǎn)桿單元圖2

扭轉(zhuǎn)桿單元示意圖設(shè)扭轉(zhuǎn)桿單元的長(zhǎng)度為，截面慣性矩為，剪切模量為，桿端扭矩分別為、，桿端扭轉(zhuǎn)角分別為、，單元上的分布荷載集度為，則任意截面的扭轉(zhuǎn)角為位移函數(shù)求得如（一）式中為局部坐標(biāo)系下扭轉(zhuǎn)桿單元的結(jié)點(diǎn)位移矩陣

（3-11）

2020/10/2830二扭轉(zhuǎn)桿單元圖2扭轉(zhuǎn)桿單元示意圖設(shè)扭轉(zhuǎn)桿單元由材料力學(xué)可知，截面扭矩為：

式中：

我們利用極小勢(shì)能原理來(lái)進(jìn)行單元分析，桿單元的勢(shì)能用泛函表示為：

內(nèi)力勢(shì)能外力勢(shì)能其中為局部坐標(biāo)系下扭轉(zhuǎn)桿單元的結(jié)點(diǎn)集中荷載矩陣

（3-12）

2020/10/2831內(nèi)力勢(shì)能外力勢(shì)能其中由極小勢(shì)能原理，取上述泛函的變分，可得：

或者寫(xiě)為：

設(shè)：

可得扭轉(zhuǎn)桿單元的單元?jiǎng)偠确匠虨椋?/p>

可以看到，其形式與拉壓桿單元的單元?jiǎng)偠确匠掏耆恢?。同樣，由上式可以進(jìn)一步求得其局部坐標(biāo)系下得單元?jiǎng)偠染仃嚍椋?/p>

（3-13a）

（3-13b)

（3-14、3-15）

（3-16）

（3-17）

等效結(jié)點(diǎn)荷載2020/10/2832（3-13a）（3-13b)（3-14、3-15）三只計(jì)彎曲的桿單元

設(shè)桿單元的長(zhǎng)度為，截面慣性矩為，彈性模量為，桿端集中剪力為、，桿端集中彎矩分別為、，桿端橫向位移為、，桿端扭轉(zhuǎn)角分別為、，在單元上分布有荷載集度為的豎向分布荷載和集度為的分布力偶。2020/10/2833三只計(jì)彎曲的桿單元2020/10/2833根據(jù)梁的平截面假定可知平面純彎梁?jiǎn)卧妮S向應(yīng)變?yōu)椋哼@里利用平截面假設(shè)（變形后橫截面仍保持平面，與縱線正交）如圖：

材料力學(xué)基礎(chǔ)知識(shí)2020/10/2834根據(jù)梁的平截面假定可知平面純彎梁?jiǎn)卧妮S向應(yīng)變?yōu)椋哼@里利用平

取撓曲線方程為的三次多項(xiàng)式，即單元上任意一點(diǎn)的撓度為：

根據(jù)單元的位移邊界條件：時(shí)：時(shí)：可以得到式中的待定系數(shù)結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為

（3-18）

2020/10/2835可以得到式中的待定系數(shù)結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分

將系數(shù)a、b、c、d代入式，并將撓曲線方程用矩陣形式表示為：

式中

為形函數(shù)矩陣，其中：

上式為平面彎曲單元的形函數(shù)。

（3-19）

（3-20）

2020/10/2836（3-19）（3-20）2020/10/2836

根據(jù)式（2-19）確定的單元位移場(chǎng)，可得單元上某一點(diǎn)得曲率為：

截面的彎矩為：

這里：

為平面彎曲桿單元的應(yīng)變矩陣。根據(jù)虛位移原理。有：

（3-21）

2020/10/2837（3-21）2020/10/2837

則平面彎曲桿單元的單元?jiǎng)偠确匠虨椋?/p>

其中的單元?jiǎng)偠染仃嚳捎墒剑?-23）求得為：

記：（3-23）（3-22）等效結(jié)點(diǎn)荷載（3-24）（3-25）2020/10/2838記：（3-23）（3-22）等效結(jié)四平面一般桿單元（考慮拉伸、彎曲、不考慮扭轉(zhuǎn)）

桿單元的長(zhǎng)度為，截面面積為，截面慣性矩為，彈性模量為，單元的、端各有三個(gè)力為、、和、、，其對(duì)應(yīng)的位移為、、和、、，建立如圖所示的局部坐標(biāo)系，各物理量的正向如圖中所標(biāo)。2020/10/2839四平面一般桿單元（考慮拉伸、彎曲、不考慮扭轉(zhuǎn)）桿單元的

設(shè)單元上沒(méi)有荷載作用，首先考慮軸向力的作用，由于桿端軸力、只引起桿端軸向位移、，根據(jù)拉壓桿單元的單元?jiǎng)偠确匠?，有?/p>

則結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為：（3-26）（3-27）2020/10/2840則結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為：（3-2

其次，桿端彎矩、和桿端剪力、只與桿端的轉(zhuǎn)角位移、和桿端的橫向位移、有關(guān)系，根據(jù)只計(jì)彎曲桿單元的單元?jiǎng)偠确匠蹋ㄗ⒁?，由于不考慮單元上的荷載作用，故方程式中的等效結(jié)點(diǎn)荷載等于零）可得：結(jié)構(gòu)力學(xué)相關(guān)知識(shí)2020/10/2841其次，桿端彎矩、和桿端剪力這樣，上述表達(dá)式合并在一起，寫(xiě)成矩陣形式如下：

可以將上式簡(jiǎn)寫(xiě)為：（3-28）（3-29）2020/10/2842這樣，上述表達(dá)式合并在一起，寫(xiě)成矩陣形式如下：可以其中單元?jiǎng)偠染仃嚕?-30）2020/10/2843其中單元?jiǎng)偠染仃嚕?-30）2020/10/2843五空間受力桿件單元（考慮扭轉(zhuǎn)、拉伸、彎曲）對(duì)空間桿件單元，除了桿端力和結(jié)點(diǎn)位移數(shù)目較平面單元多外，其分析方法與平面桿單元類似（包含拉伸、扭轉(zhuǎn)、兩個(gè)方向彎曲）

設(shè)局部坐標(biāo)系的軸為單元的形心主軸，橫截面的兩個(gè)主軸分別為軸和軸（如圖所示）。設(shè)桿橫截面面積為，桿單元長(zhǎng)度為，在平面內(nèi)抗彎剛度為，在平面內(nèi)的抗彎剛度為，桿件的抗扭剛度為。2020/10/2844五空間受力桿件單元（考慮扭轉(zhuǎn)、拉伸、彎曲）對(duì)空間桿件

空間剛架有6個(gè)位移分量和6個(gè)結(jié)點(diǎn)力分量，設(shè)局部坐標(biāo)系下它們分別為純軸向拉壓純扭轉(zhuǎn)2020/10/2845空間剛架有6個(gè)位移分量和6個(gè)結(jié)點(diǎn)力分量，設(shè)局部坐標(biāo)系下它xoy面內(nèi)彎曲xoz面內(nèi)彎曲2020/10/2846xoy面內(nèi)彎曲xoz面內(nèi)彎曲2020/10/2846（3-31）2020/10/2847（3-31）2020/10/2847其中的單元?jiǎng)偠染仃嚳蓪?xiě)為將式（2-31）寫(xiě)成矩陣的形式有（3-32）（3-33）2020/10/2848其中的單元?jiǎng)偠染仃嚳蓪?xiě)為將式（2-31）寫(xiě)成矩陣的形式有（3六單元?jiǎng)偠染仃嚨男再|(zhì)

①

單元?jiǎng)偠染仃嚍閷?duì)稱矩陣，其元素

②

單元?jiǎng)偠染仃囍械拿總€(gè)元素代表單位桿端位移引起的桿端力。其中的任意元素的物理意義是第個(gè)桿端位移分量等于1（其余位移分量等于0）時(shí)，所引起的第個(gè)桿端力的分量值。2020/10/2849六單元?jiǎng)偠染仃嚨男再|(zhì)①單元?jiǎng)偠染仃?/p>

③一般單元的單元?jiǎng)偠染仃囀瞧娈惥仃?，它的元素組成的行列式等于零，即。根據(jù)奇異矩陣的性質(zhì)，沒(méi)有逆矩陣。也就是說(shuō)，如果給定桿端位移，根據(jù)（2-29）或（2-31）式可以求出桿端力的惟一解，但反過(guò)來(lái)，如果已知桿端力，則不能根據(jù)來(lái)確定桿端位移的惟一解。因?yàn)榧词乖跅U端力已知的情況下，由于單元兩端無(wú)任何約束，因此除出桿端自身變形外，還可以發(fā)生任意的剛體位移。舉例來(lái)說(shuō)，如果物體處于靜止?fàn)顟B(tài)，我們可以說(shuō)其處于平衡狀態(tài)，但反過(guò)來(lái)，如果物體處于平衡狀態(tài)，則我們不能說(shuō)其一定處于靜止。2020/10/2850③一般單元的單元?jiǎng)偠染仃囀瞧娈惥仃?，它的元素組成

④

單元?jiǎng)偠染仃嚲哂蟹謮K的性質(zhì)，即可以用子矩陣表示。用虛線把分為四個(gè)子矩陣，把和各分為兩個(gè)子矩陣，因此，又可以寫(xiě)為：

這里：

或

或

或

或

用子矩陣形式表示單元?jiǎng)偠染仃嚭蛦卧獎(jiǎng)偠确匠?，可以使其表達(dá)的物理意義更加明顯。在單元?jiǎng)偠染仃囍?，其任意子矩陣表示桿端力和桿端位移之間的關(guān)系。（3-34）2020/10/2851④單元?jiǎng)偠染仃嚲哂蟹謮K的性質(zhì)，即可以用子矩陣表3.3桿系結(jié)構(gòu)的整體分析（整體坐標(biāo)系）一平面問(wèn)題坐標(biāo)轉(zhuǎn)換矩陣圖

平面問(wèn)題桿端力轉(zhuǎn)換示意圖一般情況下，在進(jìn)行單元分析時(shí)是在局部坐標(biāo)下完成的。對(duì)于某一單元而，如果局部坐標(biāo)系與整體坐標(biāo)系不一致，則有單元分析的物理量必須通過(guò)坐標(biāo)轉(zhuǎn)換到整體坐標(biāo)系中，然后再進(jìn)行整體坐標(biāo)系下的分析2020/10/28523.3桿系結(jié)構(gòu)的整體分析（整體坐標(biāo)系）一平面問(wèn)題坐標(biāo)轉(zhuǎn)換

這里表示由軸到軸的角，角度轉(zhuǎn)動(dòng)的正負(fù)由右手定則確定，本書(shū)中以順時(shí)針?lè)较蜣D(zhuǎn)動(dòng)為正。在兩個(gè)坐標(biāo)系中，力偶分量保持不變，即有：

同理，對(duì)于端的桿端力，有：

根據(jù)力的投影定理，將整體坐標(biāo)下的桿端力分別投影到局部坐標(biāo)下，有如下關(guān)系（3-35）（3-36）（3-36）2020/10/2853根據(jù)力的投影定理，將整體坐標(biāo)下的桿端力分別投影到將這些式子用矩陣形式可表示為：上式可以簡(jiǎn)寫(xiě)成：這即為兩種坐標(biāo)系下單元桿端力的坐標(biāo)變換式。其坐標(biāo)轉(zhuǎn)換矩陣為：（3-38）（3-39）（3-40）2020/10/2854將這些式子用矩陣形式可表示為：上式可以簡(jiǎn)寫(xiě)成：這即為兩種坐標(biāo)從坐標(biāo)轉(zhuǎn)換矩陣的表達(dá)式可以看出，為正交矩陣，其逆矩陣等于其轉(zhuǎn)置矩陣，即有：

并且有：

式中為單位矩陣。

同樣的推導(dǎo)，可以得到兩種坐標(biāo)系下的桿端位移之間的轉(zhuǎn)換關(guān)系為：

這里和分別為局部坐標(biāo)系和整體坐標(biāo)系下的桿端位移矩陣，為前面介紹的轉(zhuǎn)換矩陣。整體坐標(biāo)下的桿端力為局部坐標(biāo)下的桿端力為（3-42）（3-41）（3-43）（3-44）（3-45）2020/10/2855從坐標(biāo)轉(zhuǎn)換矩陣的表達(dá)式可以看出，為正交因此可得：

上式兩邊同乘以，可以得到：

設(shè)

可得：

上式即為整體坐標(biāo)系下的單元?jiǎng)偠确匠?。?-47）（3-46）剛度矩陣轉(zhuǎn)換）2020/10/2856因此可得：（3-47）（3-46）剛度矩陣轉(zhuǎn)換）2020/1二空間問(wèn)題的坐標(biāo)變換（空間問(wèn)題）

考慮結(jié)點(diǎn)i在局部坐標(biāo)下的桿端力與在整體坐標(biāo)系的桿端力的關(guān)系

2020/10/2857二空間問(wèn)題的坐標(biāo)變換（空間問(wèn)題）考慮結(jié)點(diǎn)i在局部坐標(biāo)下設(shè)軸與x、y、z軸的方向余弦分別為：，，則將桿端力、、向軸投影，可以求得桿端力，即：

同理可以求得：

2020/10/2858設(shè)軸與x、y、z軸的方向余弦分別為：2020/10/2用矩陣形式可以表示為：

上式即為結(jié)點(diǎn)i的桿端力在局部坐標(biāo)系和整體坐標(biāo)系下的轉(zhuǎn)換關(guān)系，其中的矩陣

稱為關(guān)系矩陣。與上面的推導(dǎo)類似，同樣可以推出以、、表示、、，以及對(duì)于結(jié)點(diǎn)j的相對(duì)應(yīng)的轉(zhuǎn)換關(guān)系，其中轉(zhuǎn)換關(guān)系矩陣都是。綜上所述，整體坐標(biāo)系下單元桿端力矩陣與局部坐標(biāo)系下單元桿端力矩陣具有如下的關(guān)系表達(dá)式：（3-48）2020/10/2859用矩陣形式可以表示為：綜上所述，整體坐標(biāo)系下單元桿其中的為：

稱為空間坐標(biāo)系的單元轉(zhuǎn)換矩陣，它是一個(gè)正交矩陣，即：

對(duì)于桿端位移，同樣可推導(dǎo)出在兩種坐標(biāo)系中的轉(zhuǎn)換關(guān)系：

這樣，可得空間桿件單元在整體坐標(biāo)系中單元?jiǎng)偠确匠虨椋?/p>

其中表示空間單元在整體坐標(biāo)系中的單元?jiǎng)偠染仃嚒＃?-49）（3-50）（3-51）（3-52）（3-53）2020/10/2860其中的為：（3-49）（3-50）（3-51）（3三桿系結(jié)構(gòu)的整體分析對(duì)桿系結(jié)構(gòu)進(jìn)行單元分析，僅僅是有限元分析中的第一步。我們的目的是要對(duì)整個(gè)結(jié)構(gòu)進(jìn)行分析，研究結(jié)構(gòu)的整體性能。因此，在對(duì)結(jié)構(gòu)的各單元分析完成后，必須將單元分析的結(jié)果進(jìn)行整合，對(duì)結(jié)構(gòu)進(jìn)行整體分析。整體分析的過(guò)程實(shí)際上是如何將單元分析的結(jié)果進(jìn)行有效組合，建立整體剛度方程并求解結(jié)點(diǎn)位移的過(guò)程。根據(jù)對(duì)結(jié)點(diǎn)位移的編碼方式，可以采用“先處理法”和“后處理法”來(lái)建立整體剛度方程。2020/10/2861三桿系結(jié)構(gòu)的整體分析對(duì)桿系結(jié)構(gòu)進(jìn)行單元分析，僅僅是有限1后處理法

所謂后處理法，就是由單元?jiǎng)偠染仃囆纬烧w剛度矩陣，建立剛度方程后再引入支承條件，進(jìn)而求解結(jié)點(diǎn)位移的方法。運(yùn)用這種方法時(shí)，假設(shè)所有結(jié)點(diǎn)位移均為未知量，按照順序統(tǒng)一進(jìn)行編碼，如圖所示的平面桿件單元。2020/10/28621后處理法所謂后處理法，就是由單元?jiǎng)偠染仃囆纬烧w剛度矩結(jié)點(diǎn)位移矩陣為：

結(jié)點(diǎn)荷載矩陣為：

求出各單元?jiǎng)偠确匠毯?，根?jù)平衡條件和位移連續(xù)條件，可以建立整個(gè)結(jié)構(gòu)的位移法方程：2020/10/2863結(jié)點(diǎn)位移矩陣為：2020/10/2863簡(jiǎn)寫(xiě)成：

這里為結(jié)構(gòu)的整體剛度矩陣，有：

注意，在建立方程的過(guò)程中，我們假設(shè)所有結(jié)點(diǎn)都有位移。因此整個(gè)結(jié)構(gòu)在外力作用下，除了發(fā)生彈性變形外，還可能發(fā)生剛體平動(dòng)位移，這樣各結(jié)點(diǎn)位移不能唯一確定。這說(shuō)明整體剛度矩陣為一奇異矩陣，不能求逆矩陣，即根據(jù)整體剛度方程可得到無(wú)窮多個(gè)解。2020/10/2864簡(jiǎn)寫(xiě)成：注意，在建立方程的過(guò)程中，我們假設(shè)所有結(jié)實(shí)際上，在圖所示剛架中，結(jié)點(diǎn)1和結(jié)點(diǎn)4均為固定端，其三個(gè)位移分量均為0，即有：

這樣，將上述支承條件引入到方程中，對(duì)整體剛度方程進(jìn)行修改，可得：

2020/10/2865實(shí)際上，在圖所示剛架中，結(jié)點(diǎn)1和結(jié)點(diǎn)4均為固定端，其三個(gè)位移對(duì)上述方程進(jìn)行化簡(jiǎn)，可以得到兩組方程：

和

這樣，利用第1式可以求得結(jié)點(diǎn)位移和，再根據(jù)第2式可以求得支座反力和。2020/10/2866對(duì)上述方程進(jìn)行化簡(jiǎn)，可以得到兩組方程：這樣，利用第1式可以求2先（前）處理法所謂先處理法，就是先引入支承條件，根據(jù)支承條件僅對(duì)未知的自由結(jié)點(diǎn)位移分量編號(hào)，得到的位移矩陣中不包含已知的約束位移分量，即可以直接得到方程求解自由結(jié)點(diǎn)位移分量。2020/10/28672先（前）處理法所謂先處理法，就是先引入支承條件，根據(jù)如圖所示的平面剛架結(jié)構(gòu)ABCD，由于在A處和D處均為固定端，其位移為0，故位移編碼均為0，在C處為鉸接，故BC桿在C端的角位移與DC桿在C端的角位移不相同，因此在C處編兩個(gè)結(jié)點(diǎn)3和4，但結(jié)點(diǎn)3和4的橫向位移和豎向位移相同，故采用相同的編號(hào)，各結(jié)點(diǎn)位移編碼如圖所示。圖先處理法位移編碼示意圖2020/10/2868如圖所示的平面剛架結(jié)構(gòu)ABCD，由于在A處和D處均為固定3.4等效結(jié)點(diǎn)荷載與邊界條件的處理非結(jié)點(diǎn)等效荷載和邊界條件的處理是有限元分析中必須考慮的的兩個(gè)重要方面。由于只考慮結(jié)點(diǎn)荷載，因此必須將作用于單元上的非結(jié)點(diǎn)荷載轉(zhuǎn)換到節(jié)點(diǎn)上。有限元的結(jié)點(diǎn)荷載來(lái)自兩部分（1）作用于結(jié)點(diǎn)處的集中力、集中力偶，前面多次提到的直接疊加到結(jié)點(diǎn)上即可（整體坐標(biāo)下）（2）非結(jié)點(diǎn)荷載首先須在局部坐標(biāo)下等效到結(jié)點(diǎn)荷載，然后再轉(zhuǎn)到整體坐標(biāo)系下的2020/10/28693.4等效結(jié)點(diǎn)荷載與邊界條件的處理非結(jié)點(diǎn)等效荷載和邊界條件1非結(jié)點(diǎn)荷載的處理在前面的分析中，我們已經(jīng)介紹了求等效結(jié)點(diǎn)荷載的方法，如3-7式、3-15式、3-22式分別可用來(lái)求不同情況下的等效結(jié)點(diǎn)荷載。

此外，可以這樣來(lái)考慮：第一步，在局部坐標(biāo)系下求單元的固端力。對(duì)于某個(gè)單元，我們假定單元的兩端均固定，然后根據(jù)靜力平衡求得固定端的反力。第二步，根據(jù)單元固端力求單元的等效結(jié)點(diǎn)荷載。根據(jù)局部坐標(biāo)系與整體坐標(biāo)系單元桿端力的變換式，固端內(nèi)力在兩種坐標(biāo)系下的變換形式可以寫(xiě)成：

因此，整體坐標(biāo)系下的等效結(jié)點(diǎn)荷載矩陣可以由下式計(jì)算：

（3-7）（3-15）（3-22）（拉壓桿）純彎桿扭轉(zhuǎn)桿2020/10/28701非結(jié)點(diǎn)荷載的處理在前面的分析中，我們已經(jīng)介紹了求等效結(jié)平面剛架單元固端力（遵行結(jié)構(gòu)力學(xué)）2020/10/2871平面剛架單元固端力（遵行結(jié)構(gòu)力學(xué)）2020/10/28712020/10/28722020/10/28722邊界條件的處理

（1）鉸結(jié)點(diǎn)在桿系結(jié)構(gòu)中，除了剛性結(jié)點(diǎn)外，通常會(huì)遇到一些桿件通過(guò)鉸結(jié)點(diǎn)與其它桿件聯(lián)結(jié)，如下圖所示桿件系統(tǒng)，有4根桿件匯交于D點(diǎn)，其中BD桿在D端通過(guò)鉸支座與其它桿件鉸接，其余3根桿為剛性接觸。對(duì)于這樣的鉸結(jié)點(diǎn)，具有如下的性質(zhì)：①鉸結(jié)點(diǎn)上各桿具有相同的線位移，但截面的轉(zhuǎn)角位移不相同；②結(jié)點(diǎn)上具有鉸接桿端不承受彎矩作用。如下圖所示結(jié)構(gòu)中，BD桿在D端的桿端彎矩為0，只有CD、ED、GD桿在結(jié)點(diǎn)D上與外彎矩保持平衡。2020/10/28732邊界條件的處理（1）鉸結(jié)點(diǎn)2020/10/2873對(duì)于這樣的結(jié)點(diǎn)，我們?cè)趯?duì)其進(jìn)行單元?jiǎng)澐謺r(shí)，通?？紤]在D處設(shè)置2個(gè)結(jié)點(diǎn)。按照先處理法，對(duì)圖示結(jié)構(gòu)進(jìn)行位移編碼，如圖2b所示。圖鉸結(jié)點(diǎn)的處理示意圖ANSYS通過(guò)結(jié)點(diǎn)耦合實(shí)現(xiàn)2020/10/2874對(duì)于這樣的結(jié)點(diǎn)，我們?cè)趯?duì)其進(jìn)行單元?jiǎng)澐謺r(shí)，通?？紤]在D處2彈性支承點(diǎn)在實(shí)際工程中，有時(shí)會(huì)遇到彈性支承的情況（如圖），這時(shí)一般將彈性支座看作是在結(jié)構(gòu)約束點(diǎn)沿約束方向的一個(gè)彈簧，彈簧的剛度系數(shù)為，在數(shù)值上等于使彈簧支座沿約束方向產(chǎn)生單位位移時(shí)所需施加的力。ANSYS引入彈簧單元即可2020/10/28752彈性支承點(diǎn)ANSYS引入彈簧單元即可2020/10/28具體做法可以歸結(jié)為：先解除彈性支承點(diǎn)約束，在i處給一個(gè)結(jié)點(diǎn)號(hào)，形成總剛度矩陣，然后在總剛度矩陣中將第i行的主元素加上彈性支承的剛度系數(shù)，此時(shí)第行變?yōu)椋阂陨系姆治鲆策m用與角位移為彈性約束的情況。若有多個(gè)彈性支座，可同時(shí)引入，即只需將相應(yīng)的主對(duì)角線元素加上相應(yīng)的彈性剛度系數(shù)即可。2020/10/2876具體做法可以歸結(jié)為：先解除彈性支承點(diǎn)約束，在i處給一個(gè)結(jié)點(diǎn)號(hào)例13.5桿系結(jié)構(gòu)有限元計(jì)算實(shí)例2020/10/2877例13.5桿系結(jié)構(gòu)有限元計(jì)算實(shí)例2020/10/281結(jié)構(gòu)的離散化與編號(hào)2020/10/28781結(jié)構(gòu)的離散化與編號(hào)2020/10/28782各個(gè)單元的矩陣描述結(jié)構(gòu)包括有斜桿，所以必須在總體坐標(biāo)下對(duì)節(jié)點(diǎn)位移進(jìn)行表達(dá)，所推導(dǎo)的單元?jiǎng)偠染仃囈惨M(jìn)行變換2020/10/28792各個(gè)單元的矩陣描述結(jié)構(gòu)包括有斜桿，所以3建立整體剛度方程1.將所得到的各個(gè)單元?jiǎng)偠染仃嚢垂?jié)點(diǎn)編號(hào)進(jìn)行組裝，可以形成整體剛度矩陣;2.同時(shí)將所有節(jié)點(diǎn)載荷也進(jìn)行組裝。

2020/10/28803建立整體剛度方程1.將所得到的各個(gè)單元?jiǎng)偠染仃嚢垂?jié)點(diǎn)編4邊界條件的處理及剛度方程求解2020/10/28814邊界條件的處理及剛度方程求解2020/10/28812020/10/28822020/10/28825各單元應(yīng)力的計(jì)算2020/10/28835各單元應(yīng)力的計(jì)算2020/10/28836支反力的計(jì)算將節(jié)點(diǎn)位移的結(jié)果代入整體剛度方程中2020/10/28846支反力的計(jì)算將節(jié)點(diǎn)位移的結(jié)果代入整體剛度方程中2020對(duì)于單元2：取i=1，j=2，則，故

對(duì)于單元1：取i=3，j=1，則c=1，s=0，故

對(duì)于單元3：取i=2，j=3，則c=0，s=1，故

例22020/10/2885對(duì)于單元2：取i=1，j=2，則，整體編號(hào)，對(duì)號(hào)入座得總剛2020/10/2886整體編號(hào)，對(duì)號(hào)入座得總剛2020/10/2886消除邊界條件后有求解左式即可得到相關(guān)位移、應(yīng)力、應(yīng)變等等2020/10/2887消除邊界條件后有求解左式即可得到相關(guān)位移、應(yīng)力、應(yīng)變例3其中E1=3e4MPa,E2=3e4MPa解(1)單元?jiǎng)澐郑⒕植孔鴺?biāo)系和整體坐標(biāo)系，并對(duì)數(shù)據(jù)進(jìn)行整理，對(duì)單元和結(jié)點(diǎn)編號(hào)2020/10/2888例3其中E1=3e4MPa,E2=3e4MPa解(1)單元?jiǎng)?2)求局部坐標(biāo)系下的單元?jiǎng)偠染仃?，由于單?和2的尺寸完全一樣，因此其單元?jiǎng)偠染仃囈粯樱瑸榍笳w坐標(biāo)系下的單元?jiǎng)偠染仃?，由于單?和整體坐標(biāo)為-90角，單元2的坐標(biāo)和整體坐標(biāo)完全一樣，單元的轉(zhuǎn)換矩陣為2020/10/2889(2)求局部坐標(biāo)系下的單元?jiǎng)偠染仃?，由于單?和整體坐標(biāo)系下的單元?jiǎng)偠染仃?020/10/2890整體坐標(biāo)系下的單元?jiǎng)偠染仃?020/10/2890整體坐標(biāo)系下的總剛度矩陣2020/10/2891整體坐標(biāo)系下的總剛度矩陣2020/10/2891（4）求總結(jié)點(diǎn)荷載

（a)等效結(jié)點(diǎn)荷載：

首先求各單元在局部坐標(biāo)下的固端內(nèi)力，對(duì)單元1（局部y坐標(biāo)向左），其q=-12kN/m，l=5m，對(duì)于單元2有F=8kN，得到局部坐標(biāo)下的結(jié)點(diǎn)力向量利用轉(zhuǎn)換坐標(biāo)矩陣，得到整體坐標(biāo)下的單元等效結(jié)點(diǎn)荷載為（b)直接作用在結(jié)點(diǎn)上的荷載（C)整體坐標(biāo)下的結(jié)點(diǎn)荷載向量2020/10/2892（4）求總結(jié)點(diǎn)荷載利用轉(zhuǎn)換坐標(biāo)矩陣，得到整體坐標(biāo)下的單元等（5）建立整體平衡方程2020/10/2893（5）建立整體平衡方程2020/10/2893（6）求解2020/10/2894（6）求解2020/10/2894（7）計(jì)算各單元的桿端內(nèi)力，首先從求出的結(jié)點(diǎn)位移中取出各單元在整體坐標(biāo)下的桿端位移，有然后計(jì)算桿端力（桿端力一定要恢復(fù)到局部坐標(biāo)下完成）2020/10/2895（7）計(jì)算各單元的桿端內(nèi)力，首先從求出的結(jié)點(diǎn)位移中取出各單2020/10/28962020/10/2896（c）剪力圖

（b）彎矩圖（d）軸力圖（a）位移曲線2020/10/2897（c）剪力圖（b）彎例42020/10/2898例42020/10/2898(1)單元?jiǎng)澐?，建立局部坐?biāo)系和整體坐標(biāo)系，其中局部坐標(biāo)下的x方向如箭頭所示，y方向向下，z軸依據(jù)右手法則（2)形成局部坐標(biāo)下的單剛2020/10/2899(1)單元?jiǎng)澐郑⒕植孔鴺?biāo)系和整體坐標(biāo)系，其中局部坐標(biāo)下的將以上數(shù)據(jù)代入單剛矩陣且有2020/10/28100將以上數(shù)據(jù)代入單剛矩陣且有2020/10/281002020/10/281012020/10/281012020/10/281022020/10/28102（3）單元位移分量：依據(jù)先處理法：即先引入支承條件，根據(jù)支承條件僅對(duì)未知的自由結(jié)點(diǎn)位移分量進(jìn)行編號(hào)，得到位移分量中不包含已知約束的位移分量（4）整體剛度矩陣：依據(jù)位移分量的對(duì)應(yīng)情況，可以得到整體剛度矩陣，疊加同位置的單剛（對(duì)應(yīng)零位移分量去掉）2020/10/28103（3）單元位移分量：依據(jù)先處理法：即先引入支承條件，根據(jù)依據(jù)位移分量的對(duì)應(yīng)情況，可以得到坐標(biāo)下的結(jié)點(diǎn)力（對(duì)應(yīng)零位移分量去掉）2020/10/28104依據(jù)位移分量的對(duì)應(yīng)情況，可以得到坐標(biāo)下的結(jié)點(diǎn)力（對(duì)應(yīng)零位移分2020/10/281052020/10/28105第3章桿系結(jié)構(gòu)的有限元法

2020/10/28106第3章桿系結(jié)構(gòu)的有限元法

2020/10/281桿、梁?jiǎn)卧攀鲇懻摋U梁?jiǎn)卧陀伤鼈兘M成的平面和空間桿梁結(jié)構(gòu)系統(tǒng).從構(gòu)造上來(lái)說(shuō)其長(zhǎng)度遠(yuǎn)大于其截面尺寸的一維構(gòu)件承受軸力或扭矩的桿件成為桿桿梁?jiǎn)栴}都有精確解（且是唯一的）承受橫向力和彎矩的桿件稱為梁平面桁架平面剛架連續(xù)梁空間剛架空間桁架等承受軸力或扭矩的桿件稱為桿將承受橫向力和彎矩的桿件稱為梁變截面桿和彎曲桿件2020/10/28107桿、梁?jiǎn)卧攀鲇懻摋U梁?jiǎn)卧陀伤鼈兘M成的平面和空間桿梁結(jié)構(gòu)系精品資料2020/10/28108精品資料2020/10/283本章主要內(nèi)容2020/10/28109本章主要內(nèi)容2020/10/2843.1有限元分析的完整過(guò)程E1=E2=2E7PaA1=A2=2cm2l1=l2=10cmP3為10N作用下二桿結(jié)構(gòu)的變形。E1、A1E2、A2說(shuō)明：u1、u2、u2分別表示節(jié)點(diǎn)1、2、3的水平位移2020/10/281103.1有限元分析的完整過(guò)程E1=E2=2E7PaA1=A2=1）用標(biāo)準(zhǔn)化的分段小單元來(lái)逼近原結(jié)構(gòu)2）尋找能夠滿足位移邊界條件的許可位移場(chǎng)3）基于位移場(chǎng)的最小勢(shì)能原理來(lái)求解基本變量為：節(jié)點(diǎn)位移內(nèi)部各點(diǎn)位移應(yīng)變應(yīng)力(1)(3)(2)問(wèn)題的解題思路2020/10/281111）用標(biāo)準(zhǔn)化的分段小單元來(lái)逼近原結(jié)構(gòu)節(jié)點(diǎn)位移內(nèi)部各點(diǎn)位移應(yīng)變完整的求解過(guò)程1）結(jié)構(gòu)離散化該構(gòu)件由兩根桿件做成，因此可以自然離散成2個(gè)桿單元。假定以這類單元位移的特征為兩個(gè)端點(diǎn)位移，就這兩個(gè)離散單元給出節(jié)點(diǎn)編號(hào)和單元編號(hào)。單元1：i=1，j=2單元2：i=2，j=32020/10/28112完整的求解過(guò)程1）結(jié)構(gòu)離散化2020/10/287單元位移模式：u(x)=a0+a1x單元節(jié)點(diǎn)條件：u(0)=u1，u(l)=u2

將式（b）代入式（a），從而得2）單元分析（a）（b）2020/10/281132）單元分析（a）（b）2020/10/288回代得寫(xiě)成矩陣形式為其中Ni，Nj是形函數(shù)。形函數(shù)矩陣說(shuō)明：u表示位移列陣ue表示單元位移2020/10/28114回代得根據(jù)幾何方程可得單元應(yīng)變的表達(dá)單元應(yīng)變寫(xiě)成矩陣形式為簡(jiǎn)記為幾何函數(shù)矩陣或者是應(yīng)變轉(zhuǎn)換矩陣2020/10/28115根據(jù)幾何方程可得單元應(yīng)變的表達(dá)幾何函數(shù)矩陣或者是應(yīng)變轉(zhuǎn)換矩陣根據(jù)物理方程可得單元應(yīng)力的表達(dá)單元應(yīng)力寫(xiě)成矩陣形式為簡(jiǎn)記為單元應(yīng)力矩陣或者是應(yīng)力轉(zhuǎn)換矩陣節(jié)點(diǎn)位移列陣2020/10/28116根據(jù)物理方程可得單元應(yīng)力的表達(dá)單元應(yīng)力矩陣或者是應(yīng)力轉(zhuǎn)換矩陣單元e勢(shì)能的表達(dá)說(shuō)明積分域，P1、P2、分別表示作用單元e上的節(jié)點(diǎn)在u1、u2的力2020/10/28117單元e勢(shì)能的表達(dá)說(shuō)明積分域，P1、P2、分別表示寫(xiě)成矩陣形式為單元e剛度矩陣單元e節(jié)點(diǎn)力列陣2020/10/28118寫(xiě)成矩陣形式為單元e剛度矩陣單元e節(jié)點(diǎn)力列陣2020/10/在得到各個(gè)單元的勢(shì)能表達(dá)式后，需要進(jìn)行離散單元的裝配，以求出整個(gè)系統(tǒng)的總勢(shì)能，對(duì)于該系統(tǒng)，總勢(shì)能包括兩個(gè)單元部分3）離散單元的裝配2020/10/28119在得到各個(gè)單元的勢(shì)能表達(dá)式后，需要進(jìn)行離散單元的裝配，以求出處理邊界條件是獲取可能位移場(chǎng)，將左端的約束條件，即u1=0代入上式可以得到簡(jiǎn)化的勢(shì)能表達(dá)式4）邊界條件的處理2020/10/28120處理邊界條件是獲取可能位移場(chǎng)，將左端的約束條件，即u1=0代由于上式是基于許可位移場(chǎng)的表達(dá)的系統(tǒng)勢(shì)能，這是由全部節(jié)點(diǎn)位移分段所插值出的位移場(chǎng)為全場(chǎng)許位移場(chǎng)，且基本未知量為節(jié)點(diǎn)位移，根據(jù)最小勢(shì)能原理（即針對(duì)未知位移求一階變分）有5）建立剛度方程2020/10/28121由于上式是基于許可位移場(chǎng)的表達(dá)的系統(tǒng)勢(shì)能，這是由全部節(jié)點(diǎn)位5將結(jié)構(gòu)參數(shù)和外載荷代入上式有求解得（單位m）6）求解節(jié)點(diǎn)位移2020/10/28122將結(jié)構(gòu)參數(shù)和外載荷代入上式有6）求解節(jié)點(diǎn)位移2020/10/7）計(jì)算單元應(yīng)變2020/10/281237）計(jì)算單元應(yīng)變2020/10/28188）計(jì)算單元應(yīng)力2020/10/281248）計(jì)算單元應(yīng)力2020/10/2819對(duì)于單元?jiǎng)菽艿谋磉_(dá)，對(duì)其取極值有具體地對(duì)于單元1，有其中R1是節(jié)點(diǎn)1的支反力，P2是單元1的節(jié)點(diǎn)2所受的力，即單元2對(duì)該節(jié)點(diǎn)的作用力，將前面求得的節(jié)點(diǎn)位移代入上式可得支反力大小。9）計(jì)算支反力2020/10/28125對(duì)于單元?jiǎng)菽艿谋磉_(dá)，對(duì)其取極值有9）計(jì)算支反力2020/10以上是一個(gè)簡(jiǎn)單結(jié)構(gòu)有限元方法求解得完整過(guò)程，對(duì)于復(fù)雜結(jié)構(gòu)，其求解過(guò)程完全相同，由于每一個(gè)步驟都具備標(biāo)準(zhǔn)化和規(guī)范性的特征，所以可以在計(jì)算機(jī)上編程而自動(dòng)實(shí)現(xiàn)。討論1：對(duì)于一個(gè)單元的勢(shì)能取極值，所得到的方程為節(jié)點(diǎn)的位移和節(jié)點(diǎn)力之間的關(guān)系，也稱為單元的平衡關(guān)系，由此可以求出每一個(gè)單元所受的節(jié)點(diǎn)力。2020/10/28126以上是一個(gè)簡(jiǎn)單結(jié)構(gòu)有限元方法求解得完整過(guò)程，對(duì)于復(fù)雜結(jié)構(gòu)，其討論2：由前面的步驟，我們也可以直接將各個(gè)單元的剛度矩陣按照節(jié)點(diǎn)編號(hào)的對(duì)應(yīng)位置來(lái)進(jìn)行裝配，即在未處理邊界條件之前，先形成整體剛度矩陣。其物理意義是，表示在未處理邊界條件前的基于節(jié)點(diǎn)描述的總體平衡關(guān)系。在對(duì)該方程進(jìn)行位移邊界條件的處理后就可以求解，這樣與先處理邊界條件再求系統(tǒng)勢(shì)能的最小值所獲得的方程完全相同。2020/10/28127討論2：由前面的步驟，我們也可以直接將各個(gè)單元的剛度矩陣按照

小結(jié)有限元分析的基本步驟及表達(dá)式1、物體幾何區(qū)域的離散化2、單元的研究（所有力學(xué)信息都用節(jié)點(diǎn)位移）來(lái)表達(dá)3、裝配集成4、邊界條件的處理并求解節(jié)點(diǎn)位移5、支反力的求取以及其它力學(xué)量（應(yīng)力、應(yīng)變及位移三大物理量）的計(jì)算2020/10/28128小結(jié)有限元分析的基本步驟及表達(dá)式1、物體幾有限元分析的基本步驟及表達(dá)式2020/10/28129有限元分析的基本步驟及表達(dá)式2020/10/2824一拉壓桿單元圖2.1拉壓桿單元示意圖

設(shè)桿單元長(zhǎng)度為，橫截面面積為，單元材料的彈性模量為，在局部坐標(biāo)系中桿端荷載分別為和，桿端位移分別為和，單元上的軸向分布荷載為。3.2局部坐標(biāo)下的桿單元分析2020/10/28130一拉壓桿單元圖2.1拉壓桿單元示意圖設(shè)桿單元長(zhǎng)度為用結(jié)點(diǎn)位移表示單元上任意截面的位移。對(duì)拉壓桿單元，可以取其位移為一次多項(xiàng)式，即

由位移的邊界條件：可得系數(shù)、為：這樣，截面任意一點(diǎn)的位移為：

用矩陣表示為：

其中

（3-1）

（3-2）

①單元位移模式。2020/10/28131用結(jié)點(diǎn)位移表示單元上任意截面的位移。對(duì)拉壓桿單元，可以取根據(jù)材料力學(xué)中應(yīng)變的定義，有

這里為應(yīng)變矩陣。由虎克定律，其應(yīng)力為：

（3-4）

（3-3）

②

進(jìn)行應(yīng)力、應(yīng)變分析其中2020/10/28132根據(jù)材料力學(xué)中應(yīng)變的定義，有（3-4）（利用虛位移原理求單元?jiǎng)偠染仃?，設(shè)桿端i、j分別產(chǎn)生虛位移、，則由此引起的桿軸截面任意位置的虛位移為：對(duì)應(yīng)的虛應(yīng)變?yōu)椋?/p>

根據(jù)虛位移原理虛功方程（力乘以虛位移得虛功、外力虛功等于變形虛功），有：

將上式整理得：

（3-5）

（3-6）

③

求單元?jiǎng)偠染仃?020/10/28133（3-5）（3-6）③求單元?jiǎng)偠染仃?020/式中：為局部坐標(biāo)系下單元結(jié)點(diǎn)荷載矩陣。

記

則可以得到拉壓桿單元的單元?jiǎng)偠确匠虨椋?/p>

這里為局部坐標(biāo)系下的單元?jiǎng)偠染仃?，為局部坐?biāo)系下等效結(jié)點(diǎn)荷載矩陣，但值得指出的是：分布荷載中可以包含集中荷載。根據(jù)定義，可以進(jìn)一步求得單元?jiǎng)偠染仃嚍椋?/p>

（3-10）

（3-7）

（3-8）

（3-9）

等效結(jié)點(diǎn)荷載2020/10/28134式中二扭轉(zhuǎn)桿單元圖2

扭轉(zhuǎn)桿單元示意圖設(shè)扭轉(zhuǎn)桿單元的長(zhǎng)度為，截面慣性矩為，剪切模量為，桿端扭矩分別為、，桿端扭轉(zhuǎn)角分別為、，單元上的分布荷載集度為，則任意截面的扭轉(zhuǎn)角為位移函數(shù)求得如（一）式中為局部坐標(biāo)系下扭轉(zhuǎn)桿單元的結(jié)點(diǎn)位移矩陣

（3-11）

2020/10/28135二扭轉(zhuǎn)桿單元圖2扭轉(zhuǎn)桿單元示意圖設(shè)扭轉(zhuǎn)桿單元由材料力學(xué)可知，截面扭矩為：

式中：

我們利用極小勢(shì)能原理來(lái)進(jìn)行單元分析，桿單元的勢(shì)能用泛函表示為：

內(nèi)力勢(shì)能外力勢(shì)能其中為局部坐標(biāo)系下扭轉(zhuǎn)桿單元的結(jié)點(diǎn)集中荷載矩陣

（3-12）

2020/10/28136內(nèi)力勢(shì)能外力勢(shì)能其中由極小勢(shì)能原理，取上述泛函的變分，可得：

或者寫(xiě)為：

設(shè)：

可得扭轉(zhuǎn)桿單元的單元?jiǎng)偠确匠虨椋?/p>

可以看到，其形式與拉壓桿單元的單元?jiǎng)偠确匠掏耆恢?。同樣，由上式可以進(jìn)一步求得其局部坐標(biāo)系下得單元?jiǎng)偠染仃嚍椋?/p>

（3-13a）

（3-13b)

（3-14、3-15）

（3-16）

（3-17）

等效結(jié)點(diǎn)荷載2020/10/28137（3-13a）（3-13b)（3-14、3-15）三只計(jì)彎曲的桿單元

設(shè)桿單元的長(zhǎng)度為，截面慣性矩為，彈性模量為，桿端集中剪力為、，桿端集中彎矩分別為、，桿端橫向位移為、，桿端扭轉(zhuǎn)角分別為、，在單元上分布有荷載集度為的豎向分布荷載和集度為的分布力偶。2020/10/28138三只計(jì)彎曲的桿單元2020/10/2833根據(jù)梁的平截面假定可知平面純彎梁?jiǎn)卧妮S向應(yīng)變?yōu)椋哼@里利用平截面假設(shè)（變形后橫截面仍保持平面，與縱線正交）如圖：

材料力學(xué)基礎(chǔ)知識(shí)2020/10/28139根據(jù)梁的平截面假定可知平面純彎梁?jiǎn)卧妮S向應(yīng)變?yōu)椋哼@里利用平

取撓曲線方程為的三次多項(xiàng)式，即單元上任意一點(diǎn)的撓度為：

根據(jù)單元的位移邊界條件：時(shí)：時(shí)：可以得到式中的待定系數(shù)結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為

（3-18）

2020/10/28140可以得到式中的待定系數(shù)結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分

將系數(shù)a、b、c、d代入式，并將撓曲線方程用矩陣形式表示為：

式中

為形函數(shù)矩陣，其中：

上式為平面彎曲單元的形函數(shù)。

（3-19）

（3-20）

2020/10/28141（3-19）（3-20）2020/10/2836

根據(jù)式（2-19）確定的單元位移場(chǎng)，可得單元上某一點(diǎn)得曲率為：

截面的彎矩為：

這里：

為平面彎曲桿單元的應(yīng)變矩陣。根據(jù)虛位移原理。有：

（3-21）

2020/10/28142（3-21）2020/10/2837

則平面彎曲桿單元的單元?jiǎng)偠确匠虨椋?/p>

其中的單元?jiǎng)偠染仃嚳捎墒剑?-23）求得為：

記：（3-23）（3-22）等效結(jié)點(diǎn)荷載（3-24）（3-25）2020/10/28143記：（3-23）（3-22）等效結(jié)四平面一般桿單元（考慮拉伸、彎曲、不考慮扭轉(zhuǎn)）

桿單元的長(zhǎng)度為，截面面積為，截面慣性矩為，彈性模量為，單元的、端各有三個(gè)力為、、和、、，其對(duì)應(yīng)的位移為、、和、、，建立如圖所示的局部坐標(biāo)系，各物理量的正向如圖中所標(biāo)。2020/10/28144四平面一般桿單元（考慮拉伸、彎曲、不考慮扭轉(zhuǎn)）桿單元的

設(shè)單元上沒(méi)有荷載作用，首先考慮軸向力的作用，由于桿端軸力、只引起桿端軸向位移、，根據(jù)拉壓桿單元的單元?jiǎng)偠确匠?，有?/p>

則結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為：（3-26）（3-27）2020/10/28145則結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為：（3-2

其次，桿端彎矩、和桿端剪力、只與桿端的轉(zhuǎn)角位移、和桿端的橫向位移、有關(guān)系，根據(jù)只計(jì)彎曲桿單元的單元?jiǎng)偠确匠蹋ㄗ⒁?，由于不考慮單元上的荷載作用，故方程式中的等效結(jié)點(diǎn)荷載等于零）可得：結(jié)構(gòu)力學(xué)相關(guān)知識(shí)2020/10/28146其次，桿端彎矩、和桿端剪力這樣，上述表達(dá)式合并在一起，寫(xiě)成矩陣形式如下：

可以將上式簡(jiǎn)寫(xiě)為：（3-28）（3-29）2020/10/28147這樣，上述表達(dá)式合并在一起，寫(xiě)成矩陣形式如下：可以其中單元?jiǎng)偠染仃嚕?-30）2020/10/28148其中單元?jiǎng)偠染仃嚕?-30）2020/10/2843五空間受力桿件單元（考慮扭轉(zhuǎn)、拉伸、彎曲）對(duì)空間桿件單元，除了桿端力和結(jié)點(diǎn)位移數(shù)目較平面單元多外，其分析方法與平面桿單元類似（包含拉伸、扭轉(zhuǎn)、兩個(gè)方向彎曲）

設(shè)局部坐標(biāo)系的軸為單元的形心主軸，橫截面的兩個(gè)主軸分別為軸和軸（如圖所示）。設(shè)桿橫截面面積為，桿單元長(zhǎng)度為，在平面內(nèi)抗彎剛度為，在平面內(nèi)的抗彎剛度為，桿件的抗扭剛度為。2020/10/28149五空間受力桿件單元（考慮扭轉(zhuǎn)、拉伸、彎曲）對(duì)空間桿件

空間剛架有6個(gè)位移分量和6個(gè)結(jié)點(diǎn)力分量，設(shè)局部坐標(biāo)系下它們分別為純軸向拉壓純扭轉(zhuǎn)2020/10/28150空間剛架有6個(gè)位移分量和6個(gè)結(jié)點(diǎn)力分量，設(shè)局部坐標(biāo)系下它xoy面內(nèi)彎曲xoz面內(nèi)彎曲2020/10/28151xoy面內(nèi)彎曲xoz面內(nèi)彎曲2020/10/2846（3-31）2020/10/28152（3-31）2020/10/2847其中的單元?jiǎng)偠染仃嚳蓪?xiě)為將式（2-31）寫(xiě)成矩陣的形式有（3-32）（3-33）2020/10/28153其中的單元?jiǎng)偠染仃嚳蓪?xiě)為將式（2-31）寫(xiě)成矩陣的形式有（3六單元?jiǎng)偠染仃嚨男再|(zhì)

①

單元?jiǎng)偠染仃嚍閷?duì)稱矩陣，其元素

②

單元?jiǎng)偠染仃囍械拿總€(gè)元素代表單位桿端位移引起的桿端力。其中的任意元素的物理意義是第個(gè)桿端位移分量等于1（其余位移分量等于0）時(shí)，所引起的第個(gè)桿端力的分量值。2020/10/28154六單元?jiǎng)偠染仃嚨男再|(zhì)①單元?jiǎng)偠染仃?/p>

③一般單元的單元?jiǎng)偠染仃囀瞧娈惥仃?，它的元素組成的行列式等于零，即。根據(jù)奇異矩陣的性質(zhì)，沒(méi)有逆矩陣。也就是說(shuō)，如果給定桿端位移，根據(jù)（2-29）或（2-31）式可以求出桿端力的惟一解，但反過(guò)來(lái)，如果已知桿端力，則不能根據(jù)來(lái)確定桿端位移的惟一解。因?yàn)榧词乖跅U端力已知的情況下，由于單元兩端無(wú)任何約束，因此除出桿端自身變形外，還可以發(fā)生任意的剛體位移。舉例來(lái)說(shuō)，如果物體處于靜止?fàn)顟B(tài)，我們可以說(shuō)其處于平衡狀態(tài)，但反過(guò)來(lái)，如果物體處于平衡狀態(tài)，則我們不能說(shuō)其一定處于靜止。2020/10/28155③一般單元的單元?jiǎng)偠染仃囀瞧娈惥仃嚕脑亟M成

④

單元?jiǎng)偠染仃嚲哂蟹謮K的性質(zhì)，即可以用子矩陣表示。用虛線把分為四個(gè)子矩陣，把和各分為兩個(gè)子矩陣，因此，又可以寫(xiě)為：

這里：

或

或

或

或

用子矩陣形式表示單元?jiǎng)偠染仃嚭蛦卧獎(jiǎng)偠确匠?，可以使其表達(dá)的物理意義更加明顯。在單元?jiǎng)偠染仃囍?，其任意子矩陣表示桿端力和桿端位移之間的關(guān)系。（3-34）2020/10/28156④單元?jiǎng)偠染仃嚲哂蟹謮K的性質(zhì)，即可以用子矩陣表3.3桿系結(jié)構(gòu)的整體分析（整體坐標(biāo)系）一平面問(wèn)題坐標(biāo)轉(zhuǎn)換矩陣圖

平面問(wèn)題桿端力轉(zhuǎn)換示意圖一般情況下，在進(jìn)行單元分析時(shí)是在局部坐標(biāo)下完成的。對(duì)于某一單元而，如果局部坐標(biāo)系與整體坐標(biāo)系不一致，則有單元分析的物理量必須通過(guò)坐標(biāo)轉(zhuǎn)換到整體坐標(biāo)系中，然后再進(jìn)行整體坐標(biāo)系下的分析2020/10/281573.3桿系結(jié)構(gòu)的整體分析（整體坐標(biāo)系）一平面問(wèn)題坐標(biāo)轉(zhuǎn)換

這里表示由軸到軸的角，角度轉(zhuǎn)動(dòng)的正負(fù)由右手定則確定，本書(shū)中以順時(shí)針?lè)较蜣D(zhuǎn)動(dòng)為正。在兩個(gè)坐標(biāo)系中，力偶分量保持不變，即有：

同理，對(duì)于端的桿端力，有：

根據(jù)力的投影定理，將整體坐標(biāo)下的桿端力分別投影到局部坐標(biāo)下，有如下關(guān)系（3-35）（3-36）（3-36）2020/10/28158根據(jù)力的投影定理，將整體坐標(biāo)下的桿端力分別投影到將這些式子用矩陣形式可表示為：上式可以簡(jiǎn)寫(xiě)成：這即為兩種坐標(biāo)系下單元桿端力的坐標(biāo)變換式。其坐標(biāo)轉(zhuǎn)換矩陣為：（3-38）（3-39）（3-40）2020/10/28159將這些式子用矩陣形式可表示為：上式可以簡(jiǎn)寫(xiě)成：這即為兩種坐標(biāo)從坐標(biāo)轉(zhuǎn)換矩陣的表達(dá)式可以看出，為正交矩陣，其逆矩陣等于其轉(zhuǎn)置矩陣，即有：

并且有：

式中為單位矩陣。

同樣的推導(dǎo)，可以得到兩種坐標(biāo)系下的桿端位移之間的轉(zhuǎn)換關(guān)系為：

這里和分別為局部坐標(biāo)系和整體坐標(biāo)系下的桿端位移矩陣，為前面介紹的轉(zhuǎn)換矩陣。整體坐標(biāo)下的桿端力為局部坐標(biāo)下的桿端力為（3-42）（3-41）（3-43）（3-44）（3-45）2020/10/28160從坐標(biāo)轉(zhuǎn)換矩陣的表達(dá)式可以看出，為正交因此可得：

上式兩邊同乘以，可以得到：

設(shè)

可得：

上式即為整體坐標(biāo)系下的單元?jiǎng)偠确匠?。?-47）（3-46）剛度矩陣轉(zhuǎn)換）2020/10/28161因此可得：（3-47）（3-46）剛度矩陣轉(zhuǎn)換）2020/1二空間問(wèn)題的坐標(biāo)變換（空間問(wèn)題）

考慮結(jié)點(diǎn)i在局部坐標(biāo)下的桿端力與在整體坐標(biāo)系的桿端力的關(guān)系

2020/10/28162二空間問(wèn)題的坐標(biāo)變換（空間問(wèn)題）考慮結(jié)點(diǎn)i在局部坐標(biāo)下設(shè)軸與x、y、z軸的方向余弦分別為：，，則將桿端力、、向軸投影，可以求得桿端力，即：

同理可以求得：

2020/10/28163設(shè)軸與x、y、z軸的方向余弦分別為：2020/10/2用矩陣形式可以表示為：

上式即為結(jié)點(diǎn)i的桿端力在局部坐標(biāo)系和整體坐標(biāo)系下的轉(zhuǎn)換關(guān)系，其中的矩陣

稱為關(guān)系矩陣。與上面的推導(dǎo)類似，同樣可以推出以、、表示、、，以及對(duì)于結(jié)點(diǎn)j的相對(duì)應(yīng)的轉(zhuǎn)換關(guān)系，其中轉(zhuǎn)換關(guān)系矩陣都是。綜上所述，整體坐標(biāo)系下單元桿端力矩陣與局部坐標(biāo)系下單元桿端力矩陣具有如下的關(guān)系表達(dá)式：（3-48）2020/10/28164用矩陣形式可以表示為：綜上所述，整體坐標(biāo)系下單元桿其中的為：

稱為空間坐標(biāo)系的單元轉(zhuǎn)換矩陣，它是一個(gè)正交矩陣，即：

對(duì)于桿端位移，同樣可推導(dǎo)出在兩種坐標(biāo)系中的轉(zhuǎn)換關(guān)系：

這樣，可得空間桿件單元在整體坐標(biāo)系中單元?jiǎng)偠确匠虨椋?/p>

其中表示空間單元在整體坐標(biāo)系中的單元?jiǎng)偠染仃?。?-49）（3-50）（3-51）（3-52）（3-53）2020/10/28165其中的為：（3-49）（3-50）（3-51）（3三桿系結(jié)構(gòu)的整體分析對(duì)桿系結(jié)構(gòu)進(jìn)行單元分析，僅僅是有限元分析中的第一步。我們的目的是要對(duì)整個(gè)結(jié)構(gòu)進(jìn)行分析，研究結(jié)構(gòu)的整體性能。因此，在對(duì)結(jié)構(gòu)的各單元分析完成后，必須將單元分析的結(jié)果進(jìn)行整合，對(duì)結(jié)構(gòu)進(jìn)行整體分析。整體分析的過(guò)程實(shí)際上是如何將單元分析的結(jié)果進(jìn)行有效組合，建立整體剛度方程并求解結(jié)點(diǎn)位移的過(guò)程。根據(jù)對(duì)結(jié)點(diǎn)位移的編碼方式，可以采用“先處理法”和“后處理法”來(lái)建立整體剛度方程。2020/10/28166三桿系結(jié)構(gòu)的整體分析對(duì)桿系結(jié)構(gòu)進(jìn)行單元分析，僅僅是有限1后處理法

所謂后處理法，就是由單元?jiǎng)偠染仃囆纬烧w剛度矩陣，建立剛度方程后再引入支承條件，進(jìn)而求解結(jié)點(diǎn)位移的方法。運(yùn)用這種方法時(shí)，假設(shè)所有結(jié)點(diǎn)位移均為未知量，按照順序統(tǒng)一進(jìn)行編碼，如圖所示的平面桿件單元。2020/10/281671后處理法所謂后處理法，就是由單元?jiǎng)偠染仃囆纬烧w剛度矩結(jié)點(diǎn)位移矩陣為：

結(jié)點(diǎn)荷載矩陣為：

求出各單元?jiǎng)偠确匠毯?，根?jù)平衡條件和位移連續(xù)條件，可以建立整個(gè)結(jié)構(gòu)的位移法方程：2020/10/28168結(jié)點(diǎn)位移矩陣為：2020/10/2863簡(jiǎn)寫(xiě)成：

這里為結(jié)構(gòu)的整體剛度矩陣，有：

注意，在建立方程的過(guò)程中，我們假設(shè)所有結(jié)點(diǎn)都有位移。因此整個(gè)結(jié)構(gòu)在外力作用下，除了發(fā)生彈性變形外，還可能發(fā)生剛體平動(dòng)位移，這樣各結(jié)點(diǎn)位移不能唯一確定。這說(shuō)明整體剛度矩陣為一奇異矩陣，不能求逆矩陣，即根據(jù)整體剛度方程可得到無(wú)窮多個(gè)解。2020/10/28169簡(jiǎn)寫(xiě)成：注意，在建立方程的過(guò)程中，我們假設(shè)所有結(jié)實(shí)際上，在圖所示剛架中，結(jié)點(diǎn)1和結(jié)點(diǎn)4均為固定端，其三個(gè)位移分量均為0，即有：

這樣，將上述支承條件引入到方程中，對(duì)整體剛度方程進(jìn)行修改，可得：

2020/10/28170實(shí)際上，在圖所示剛架中，結(jié)點(diǎn)1和結(jié)點(diǎn)4均為固定端，其三個(gè)位移對(duì)上述方程進(jìn)行化簡(jiǎn)，可以得到兩組方程：

和

這樣，利用第1式可以求得結(jié)點(diǎn)位移和，再根據(jù)第2式可以求得支座反力和。2020/10/28171對(duì)上述方程進(jìn)行化簡(jiǎn)，可以得到兩組方程：這樣，利用第1式可以求2先（前）處理法所謂先處理法，就是先引入支承條件，根據(jù)支承條件僅對(duì)未知的自由結(jié)點(diǎn)位移分量編號(hào)，得到的位移矩陣中不包含已知的約束位移分量，即可以直接得到方程求解自由結(jié)點(diǎn)位移分量。2020/10/281722先（前）處理法所謂先處理法，就是先引入支承條件，根據(jù)如圖所示的平面剛架結(jié)構(gòu)ABCD，由于在A處和D處均為固定端，其位移為0，故位移編碼均為0，在C處為鉸接，故BC桿在C端的角位移與DC桿在C端的角位移不相同，因此在C處編兩個(gè)結(jié)點(diǎn)3和4，但結(jié)點(diǎn)3和4的橫向位移和豎向位移相同，故采用相同的編號(hào)，各結(jié)點(diǎn)位移編碼如圖所示。圖先處理法位移編碼示意圖2020/10/28173如圖所示的平面剛架結(jié)構(gòu)ABCD，由于在A處和D處均為固定3.4等效結(jié)點(diǎn)荷載與邊界條件的處理非結(jié)點(diǎn)等效荷載和邊界條件的處理是有限元分析中必須考慮的的兩個(gè)重要方面。由于只考慮結(jié)點(diǎn)荷載，因此必須將作用于單元上的非結(jié)點(diǎn)荷載轉(zhuǎn)換到節(jié)點(diǎn)上。