1,第3章桿系結(jié)構(gòu)的有限元法,2,桿、梁單元概述,討論桿梁單元和由它們組成的平面和空間桿梁結(jié)構(gòu)系統(tǒng).從構(gòu)造上來說其長度遠(yuǎn)大于其截面尺寸的一維構(gòu)件承受軸力或扭矩的桿件成為桿桿梁問題都有精確解（且是唯一的）承受橫向力和彎矩的桿件稱為梁平面桁架平面剛架連續(xù)梁空間剛架空間桁架等承受軸力或扭矩的桿件稱為桿將承受橫向力和彎矩的桿件稱為梁變截面桿和彎曲桿件,3,本章主要內(nèi)容,4,3.1有限元分析的完整過程,E1=E2=2E7Pa,A1=A2=2cm2,l1=l2=10cm,P3為10N作用下二桿結(jié)構(gòu)的變形。,E1、A1,E2、A2,說明：u1、u2、u2分別表示節(jié)點(diǎn)1、2、3的水平位移,5,1）用標(biāo)準(zhǔn)化的分段小單元來逼近原結(jié)構(gòu)2）尋找能夠滿足位移邊界條件的許可位移場3）基于位移場的最小勢能原理來求解基本變量為：,問題的解題思路,6,完整的求解過程,1）結(jié)構(gòu)離散化該構(gòu)件由兩根桿件做成，因此可以自然離散成2個桿單元。假定以這類單元位移的特征為兩個端點(diǎn)位移，就這兩個離散單元給出節(jié)點(diǎn)編號和單元編號。單元1：i=1，j=2單元2：i=2，j=3,7,單元位移模式：u(x)=a0+a1x單元節(jié)點(diǎn)條件：u(0)=u1，u(l)=u2將式（b）代入式（a），從而得,2）單元分析,（a）,（b）,8,回代得寫成矩陣形式為,其中Ni，Nj是形函數(shù)。,形函數(shù)矩陣,說明：u表示位移列陣ue表示單元位移,9,根據(jù)幾何方程可得單元應(yīng)變的表達(dá)單元應(yīng)變寫成矩陣形式為簡記為,幾何函數(shù)矩陣或者是應(yīng)變轉(zhuǎn)換矩陣,10,根據(jù)物理方程可得單元應(yīng)力的表達(dá)單元應(yīng)力寫成矩陣形式為簡記為,單元應(yīng)力矩陣或者是應(yīng)力轉(zhuǎn)換矩陣,節(jié)點(diǎn)位移列陣,11,單元e勢能的表達(dá),說明積分域，P1、P2、分別表示作用單元e上的節(jié)點(diǎn)在u1、u2的力,12,寫成矩陣形式為,單元e剛度矩陣,單元e節(jié)點(diǎn)力列陣,13,在得到各個單元的勢能表達(dá)式后，需要進(jìn)行離散單元的裝配，以求出整個系統(tǒng)的總勢能，對于該系統(tǒng)，總勢能包括兩個單元部分,3）離散單元的裝配,14,處理邊界條件是獲取可能位移場，將左端的約束條件，即u1=0代入上式可以得到簡化的勢能表達(dá)式,4）邊界條件的處理,15,由于上式是基于許可位移場的表達(dá)的系統(tǒng)勢能，這是由全部節(jié)點(diǎn)位移分段所插值出的位移場為全場許位移場，且基本未知量為節(jié)點(diǎn)位移，根據(jù)最小勢能原理（即針對未知位移求一階變分）有,5）建立剛度方程,16,將結(jié)構(gòu)參數(shù)和外載荷代入上式有求解得（單位m）,6）求解節(jié)點(diǎn)位移,17,7）計(jì)算單元應(yīng)變,18,8）計(jì)算單元應(yīng)力,19,對于單元勢能的表達(dá)，對其取極值有具體地對于單元1，有其中R1是節(jié)點(diǎn)1的支反力，P2是單元1的節(jié)點(diǎn)2所受的力，即單元2對該節(jié)點(diǎn)的作用力，將前面求得的節(jié)點(diǎn)位移代入上式可得支反力大小。,9）計(jì)算支反力,20,以上是一個簡單結(jié)構(gòu)有限元方法求解得完整過程，對于復(fù)雜結(jié)構(gòu)，其求解過程完全相同，由于每一個步驟都具備標(biāo)準(zhǔn)化和規(guī)范性的特征，所以可以在計(jì)算機(jī)上編程而自動實(shí)現(xiàn)。討論1：對于一個單元的勢能取極值，所得到的方程為節(jié)點(diǎn)的位移和節(jié)點(diǎn)力之間的關(guān)系，也稱為單元的平衡關(guān)系，由此可以求出每一個單元所受的節(jié)點(diǎn)力。,21,討論2：由前面的步驟，我們也可以直接將各個單元的剛度矩陣按照節(jié)點(diǎn)編號的對應(yīng)位置來進(jìn)行裝配，即在未處理邊界條件之前，先形成整體剛度矩陣。其物理意義是，表示在未處理邊界條件前的基于節(jié)點(diǎn)描述的總體平衡關(guān)系。在對該方程進(jìn)行位移邊界條件的處理后就可以求解，這樣與先處理邊界條件再求系統(tǒng)勢能的最小值所獲得的方程完全相同。,22,小結(jié)有限元分析的基本步驟及表達(dá)式,1、物體幾何區(qū)域的離散化2、單元的研究（所有力學(xué)信息都用節(jié)點(diǎn)位移）來表達(dá)3、裝配集成4、邊界條件的處理并求解節(jié)點(diǎn)位移5、支反力的求取以及其它力學(xué)量（應(yīng)力、應(yīng)變及位移三大物理量）的計(jì)算,23,有限元分析的基本步驟及表達(dá)式,24,一拉壓桿單元,圖2.1拉壓桿單元示意圖,設(shè)桿單元長度為，橫截面面積為，單元材料的彈性模量為，在局部坐標(biāo)系中桿端荷載分別為和，桿端位移分別為和，單元上的軸向分布荷載為。,3.2局部坐標(biāo)下的桿單元分析,25,用結(jié)點(diǎn)位移表示單元上任意截面的位移。對拉壓桿單元，可以取其位移為一次多項(xiàng)式，即由位移的邊界條件：可得系數(shù)、為：這樣，截面任意一點(diǎn)的位移為：用矩陣表示為：其中,（3-1）,（3-2）,單元位移模式。,26,根據(jù)材料力學(xué)中應(yīng)變的定義，有這里為應(yīng)變矩陣。由虎克定律，其應(yīng)力為：,（3-4）,（3-3）,進(jìn)行應(yīng)力、應(yīng)變分析,其中,27,利用虛位移原理求單元剛度矩陣，設(shè)桿端i、j分別產(chǎn)生虛位移、，則由此引起的桿軸截面任意位置的虛位移為：對應(yīng)的虛應(yīng)變?yōu)椋焊鶕?jù)虛位移原理虛功方程（力乘以虛位移得虛功、外力虛功等于變形虛功），有：將上式整理得：,（3-5）,（3-6）,求單元剛度矩陣,28,式中：為局部坐標(biāo)系下單元結(jié)點(diǎn)荷載矩陣。記則可以得到拉壓桿單元的單元剛度方程為：這里為局部坐標(biāo)系下的單元剛度矩陣，為局部坐標(biāo)系下等效結(jié)點(diǎn)荷載矩陣，但值得指出的是：分布荷載中可以包含集中荷載。根據(jù)定義，可以進(jìn)一步求得單元剛度矩陣為：,（3-10）,（3-7）,（3-8）,（3-9）,等效結(jié)點(diǎn)荷載,29,二扭轉(zhuǎn)桿單元,圖2扭轉(zhuǎn)桿單元示意圖,設(shè)扭轉(zhuǎn)桿單元的長度為，截面慣性矩為，剪切模量為，桿端扭矩分別為、，桿端扭轉(zhuǎn)角分別為、，單元上的分布荷載集度為，則任意截面的扭轉(zhuǎn)角為位移函數(shù)求得如（一）,式中為局部坐標(biāo)系下扭轉(zhuǎn)桿單元的結(jié)點(diǎn)位移矩陣,（3-11）,30,由材料力學(xué)可知，截面扭矩為：式中：我們利用極小勢能原理來進(jìn)行單元分析，桿單元的勢能用泛函表示為：,內(nèi)力勢能,外力勢能,其中為局部坐標(biāo)系下扭轉(zhuǎn)桿單元的結(jié)點(diǎn)集中荷載矩陣,（3-12）,31,由極小勢能原理，取上述泛函的變分，可得：或者寫為：設(shè)：可得扭轉(zhuǎn)桿單元的單元剛度方程為：可以看到，其形式與拉壓桿單元的單元剛度方程完全一致。同樣，由上式可以進(jìn)一步求得其局部坐標(biāo)系下得單元剛度矩陣為：,（3-13a）,（3-13b),（3-14、3-15）,（3-16）,（3-17）,等效結(jié)點(diǎn)荷載,32,三只計(jì)彎曲的桿單元,設(shè)桿單元的長度為，截面慣性矩為，彈性模量為，桿端集中剪力為、，桿端集中彎矩分別為、，桿端橫向位移為、，桿端扭轉(zhuǎn)角分別為、，在單元上分布有荷載集度為的豎向分布荷載和集度為的分布力偶。,33,根據(jù)梁的平截面假定可知平面純彎梁單元的軸向應(yīng)變?yōu)椋?這里利用平截面假設(shè)（變形后橫截面仍保持平面，與縱線正交）如圖：,材料力學(xué)基礎(chǔ)知識,34,取撓曲線方程為的三次多項(xiàng)式，即單元上任意一點(diǎn)的撓度為：根據(jù)單元的位移邊界條件：時：時：,可以得到式中的待定系數(shù),結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為,（3-18）,35,將系數(shù)a、b、c、d代入式，并將撓曲線方程用矩陣形式表示為：式中為形函數(shù)矩陣，其中：上式為平面彎曲單元的形函數(shù)。,（3-19）,（3-20）,36,根據(jù)式（2-19）確定的單元位移場，可得單元上某一點(diǎn)得曲率為：截面的彎矩為：這里：為平面彎曲桿單元的應(yīng)變矩陣。根據(jù)虛位移原理。有：,（3-21）,37,則平面彎曲桿單元的單元剛度方程為：其中的單元剛度矩陣可由式（2-23）求得為：,記：,（3-23）,（3-22）,等效結(jié)點(diǎn)荷載,（3-24）,（3-25）,38,四平面一般桿單元（考慮拉伸、彎曲、不考慮扭轉(zhuǎn)）,桿單元的長度為，截面面積為，截面慣性矩為，彈性模量為，單元的、端各有三個力為、和、，其對應(yīng)的位移為、和、，建立如圖所示的局部坐標(biāo)系，各物理量的正向如圖中所標(biāo)。,39,設(shè)單元上沒有荷載作用，首先考慮軸向力的作用，由于桿端軸力、只引起桿端軸向位移、，根據(jù)拉壓桿單元的單元剛度方程，有：,則結(jié)點(diǎn)位移矩陣和結(jié)點(diǎn)荷載矩陣分別為：,（3-26）,（3-27）,40,其次，桿端彎矩、和桿端剪力、只與桿端的轉(zhuǎn)角位移、和桿端的橫向位移、有關(guān)系，根據(jù)只計(jì)彎曲桿單元的單元剛度方程（注意，由于不考慮單元上的荷載作用，故方程式中的等效結(jié)點(diǎn)荷載等于零）可得：,結(jié)構(gòu)力學(xué)相關(guān)知識,41,這樣，上述表達(dá)式合并在一起，寫成矩陣形式如下：,可以將上式簡寫為：,（3-28）,（3-29）,42,其中單元剛度矩陣,（3-30）,43,五空間受力桿件單元（考慮扭轉(zhuǎn)、拉伸、彎曲）,對空間桿件單元，除了桿端力和結(jié)點(diǎn)位移數(shù)目較平面單元多外，其分析方法與平面桿單元類似（包含拉伸、扭轉(zhuǎn)、兩個方向彎曲）,設(shè)局部坐標(biāo)系的軸為單元的形心主軸，橫截面的兩個主軸分別為軸和軸（如圖所示）。設(shè)桿橫截面面積為，桿單元長度為，在平面內(nèi)抗彎剛度為，在平面內(nèi)的抗彎剛度為，桿件的抗扭剛度為。,44,空間剛架有6個位移分量和6個結(jié)點(diǎn)力分量，設(shè)局部坐標(biāo)系下它們分別為,純軸向拉壓,純扭轉(zhuǎn),45,xoy面內(nèi)彎曲,xoz面內(nèi)彎曲,46,（3-31）,47,其中的單元剛度矩陣可寫為,將式（2-31）寫成矩陣的形式有,（3-32）,（3-33）,48,六單元剛度矩陣的性質(zhì),單元剛度矩陣為對稱矩陣，其元素單元剛度矩陣中的每個元素代表單位桿端位移引起的桿端力。其中的任意元素的物理意義是第個桿端位移分量等于1（其余位移分量等于0）時，所引起的第個桿端力的分量值。,49,一般單元的單元剛度矩陣是奇異矩陣，它的元素組成的行列式等于零，即。根據(jù)奇異矩陣的性質(zhì)，沒有逆矩陣。也就是說，如果給定桿端位移，根據(jù)（2-29）或（2-31）式可以求出桿端力的惟一解，但反過來，如果已知桿端力，則不能根據(jù)來確定桿端位移的惟一解。因?yàn)榧词乖跅U端力已知的情況下，由于單元兩端無任何約束，因此除出桿端自身變形外，還可以發(fā)生任意的剛體位移。舉例來說，如果物體處于靜止?fàn)顟B(tài)，我們可以說其處于平衡狀態(tài)，但反過來，如果物體處于平衡狀態(tài)，則我們不能說其一定處于靜止。,50,單元剛度矩陣具有分塊的性質(zhì)，即可以用子矩陣表示。用虛線把分為四個子矩陣，把和各分為兩個子矩陣，因此，又可以寫為：這里：或或或或用子矩陣形式表示單元剛度矩陣和單元剛度方程，可以使其表達(dá)的物理意義更加明顯。在單元剛度矩陣中，其任意子矩陣表示桿端力和桿端位移之間的關(guān)系。,（3-34）,51,3.3桿系結(jié)構(gòu)的整體分析（整體坐標(biāo)系）,一平面問題坐標(biāo)轉(zhuǎn)換矩陣,圖平面問題桿端力轉(zhuǎn)換示意圖,一般情況下，在進(jìn)行單元分析時是在局部坐標(biāo)下完成的。對于某一單元而，如果局部坐標(biāo)系與整體坐標(biāo)系不一致，則有單元分析的物理量必須通過坐標(biāo)轉(zhuǎn)換到整體坐標(biāo)系中，然后再進(jìn)行整體坐標(biāo)系下的分析,52,這里表示由軸到軸的角，角度轉(zhuǎn)動的正負(fù)由右手定則確定，本書中以順時針方向轉(zhuǎn)動為正。在兩個坐標(biāo)系中，力偶分量保持不變，即有：同理，對于端的桿端力，有：,根據(jù)力的投影定理，將整體坐標(biāo)下的桿端力分別投影到局部坐標(biāo)下，有如下關(guān)系,（3-35）,（3-36）,（3-36）,53,將這些式子用矩陣形式可表示為：,上式可以簡寫成：,這即為兩種坐標(biāo)系下單元桿端力的坐標(biāo)變換式。其坐標(biāo)轉(zhuǎn)換矩陣為：,（3-38）,（3-39）,（3-40）,54,從坐標(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）,55,因此可得：上式兩邊同乘以，可以得到：設(shè)可得：上式即為整體坐標(biāo)系下的單元剛度方程。,（3-47）,（3-46）,剛度矩陣轉(zhuǎn)換）,56,二空間問題的坐標(biāo)變換（空間問題）,考慮結(jié)點(diǎn)i在局部坐標(biāo)下的桿端力與在整體坐標(biāo)系的桿端力的關(guān)系,57,設(shè)軸與x、y、z軸的方向余弦分別為：，則將桿端力、向軸投影，可以求得桿端力，即：同理可以求得：,58,用矩陣形式可以表示為：上式即為結(jié)點(diǎn)i的桿端力在局部坐標(biāo)系和整體坐標(biāo)系下的轉(zhuǎn)換關(guān)系，其中的矩陣稱為關(guān)系矩陣。與上面的推導(dǎo)類似，同樣可以推出以、表示、，以及對于結(jié)點(diǎn)j的相對應(yīng)的轉(zhuǎn)換關(guān)系，其中轉(zhuǎn)換關(guān)系矩陣都是。,綜上所述，整體坐標(biāo)系下單元桿端力矩陣與局部坐標(biāo)系下單元桿端力矩陣具有如下的關(guān)系表達(dá)式：,（3-48）,59,其中的為：稱為空間坐標(biāo)系的單元轉(zhuǎn)換矩陣，它是一個正交矩陣，即：對于桿端位移，同樣可推導(dǎo)出在兩種坐標(biāo)系中的轉(zhuǎn)換關(guān)系：這樣，可得空間桿件單元在整體坐標(biāo)系中單元剛度方程為：其中表示空間單元在整體坐標(biāo)系中的單元剛度矩陣。,（3-49）,（3-50）,（3-51）,（3-52）,（3-53）,60,三桿系結(jié)構(gòu)的整體分析,對桿系結(jié)構(gòu)進(jìn)行單元分析，僅僅是有限元分析中的第一步。我們的目的是要對整個結(jié)構(gòu)進(jìn)行分析，研究結(jié)構(gòu)的整體性能。因此，在對結(jié)構(gòu)的各單元分析完成后，必須將單元分析的結(jié)果進(jìn)行整合，對結(jié)構(gòu)進(jìn)行整體分析。整體分析的過程實(shí)際上是如何將單元分析的結(jié)果進(jìn)行有效組合，建立整體剛度方程并求解結(jié)點(diǎn)位移的過程。根據(jù)對結(jié)點(diǎn)位移的編碼方式，可以采用“先處理法”和“后處理法”來建立整體剛度方程。,61,1后處理法,所謂后處理法，就是由單元剛度矩陣形成整體剛度矩陣，建立剛度方程后再引入支承條件，進(jìn)而求解結(jié)點(diǎn)位移的方法。運(yùn)用這種方法時，假設(shè)所有結(jié)點(diǎn)位移均為未知量，按照順序統(tǒng)一進(jìn)行編碼，如圖所示的平面桿件單元。,62,結(jié)點(diǎn)位移矩陣為：結(jié)點(diǎn)荷載矩陣為：求出各單元剛度方程后，根據(jù)平衡條件和位移連續(xù)條件，可以建立整個結(jié)構(gòu)的位移法方程：,63,簡寫成：這里為結(jié)構(gòu)的整體剛度矩陣，有：,注意，在建立方程的過程中，我們假設(shè)所有結(jié)點(diǎn)都有位移。因此整個結(jié)構(gòu)在外力作用下，除了發(fā)生彈性變形外，還可能發(fā)生剛體平動位移，這樣各結(jié)點(diǎn)位移不能唯一確定。這說明整體剛度矩陣為一奇異矩陣，不能求逆矩陣，即根據(jù)整體剛度方程可得到無窮多個解。,64,實(shí)際上，在圖所示剛架中，結(jié)點(diǎn)1和結(jié)點(diǎn)4均為固定端，其三個位移分量均為0，即有：這樣，將上述支承條件引入到方程中，對整體剛度方程進(jìn)行修改，可得：,65,對上述方程進(jìn)行化簡，可以得到兩組方程：和,這樣，利用第1式可以求得結(jié)點(diǎn)位移和，再根據(jù)第2式可以求得支座反力和。,66,2先（前）處理法,所謂先處理法，就是先引入支承條件，根據(jù)支承條件僅對未知的自由結(jié)點(diǎn)位移分量編號，得到的位移矩陣中不包含已知的約束位移分量，即可以直接得到方程求解自由結(jié)點(diǎn)位移分量。,67,如圖所示的平面剛架結(jié)構(gòu)ABCD，由于在A處和D處均為固定端，其位移為0，故位移編碼均為0，在C處為鉸接，故BC桿在C端的角位移與DC桿在C端的角位移不相同，因此在C處編兩個結(jié)點(diǎn)3和4，但結(jié)點(diǎn)3和4的橫向位移和豎向位移相同，故采用相同的編號，各結(jié)點(diǎn)位移編碼如圖所示。,圖先處理法位移編碼示意圖,68,3.4等效結(jié)點(diǎn)荷載與邊界條件的處理,非結(jié)點(diǎn)等效荷載和邊界條件的處理是有限元分析中必須考慮的的兩個重要方面。由于只考慮結(jié)點(diǎn)荷載，因此必須將作用于單元上的非結(jié)點(diǎn)荷載轉(zhuǎn)換到節(jié)點(diǎn)上。,有限元的結(jié)點(diǎn)荷載來自兩部分（1）作用于結(jié)點(diǎn)處的集中力、集中力偶，前面多次提到的直接疊加到結(jié)點(diǎn)上即可（整體坐標(biāo)下）（2）非結(jié)點(diǎn)荷載首先須在局部坐標(biāo)下等效到結(jié)點(diǎn)荷載，然后再轉(zhuǎn)到整體坐標(biāo)系下的,69,1非結(jié)點(diǎn)荷載的處理,在前面的分析中，我們已經(jīng)介紹了求等效結(jié)點(diǎn)荷載的方法，如3-7式、3-15式、3-22式分別可用來求不同情況下的等效結(jié)點(diǎn)荷載。此外，可以這樣來考慮：第一步，在局部坐標(biāo)系下求單元的固端力。對于某個單元，我們假定單元的兩端均固定，然后根據(jù)靜力平衡求得固定端的反力。第二步，根據(jù)單元固端力求單元的等效結(jié)點(diǎn)荷載。根據(jù)局部坐標(biāo)系與整體坐標(biāo)系單元桿端力的變換式，固端內(nèi)力在兩種坐標(biāo)系下的變換形式可以寫成：因此，整體坐標(biāo)系下的等效結(jié)點(diǎn)荷載矩陣可以由下式計(jì)算：,（3-7）,（3-15）,（3-22）,（拉壓桿）,純彎桿,扭轉(zhuǎn)桿,70,平面剛架單元固端力（遵行結(jié)構(gòu)力學(xué)）,71,72,2邊界條件的處理,（1）鉸結(jié)點(diǎn)在桿系結(jié)構(gòu)中，除了剛性結(jié)點(diǎn)外，通常會遇到一些桿件通過鉸結(jié)點(diǎn)與其它桿件聯(lián)結(jié)，如下圖所示桿件系統(tǒng)，有4根桿件匯交于D點(diǎn)，其中BD桿在D端通過鉸支座與其它桿件鉸接，其余3根桿為剛性接觸。對于這樣的鉸結(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上與外彎矩保持平衡。,73,對于這樣的結(jié)點(diǎn)，我們在對其進(jìn)行單元劃分時，通?？紤]在D處設(shè)置2個結(jié)點(diǎn)。按照先處理法，對圖示結(jié)構(gòu)進(jìn)行位移編碼，如圖2b所示。,圖鉸結(jié)點(diǎn)的處理示意圖,ANSYS通過結(jié)點(diǎn)耦合實(shí)現(xiàn),74,2彈性支承點(diǎn)在實(shí)際工程中，有時會遇到彈性支承的情況（如圖），這時一般將彈性支座看作是在結(jié)構(gòu)約束點(diǎn)沿約束方向的一個彈簧，彈簧的剛度系數(shù)為，在數(shù)值上等于使彈簧支座沿約束方向產(chǎn)生單位位移時所需施加的力。,ANSYS引入彈簧單元即可,75,具體做法可以歸結(jié)為：先解除彈性支承點(diǎn)約束，在i處給一個結(jié)點(diǎn)號，形成總剛度矩陣，然后在總剛度矩陣中將第i行的主元素加上彈性支承的剛度系數(shù)，此時第行變?yōu)椋?以上的分析也適用與角位移為彈性約束的情況。若有多個彈性支座，可同時引入，即只需將相應(yīng)的主對角線元素加上相應(yīng)的彈性剛度系數(shù)即可。,76,例1,3.5桿系結(jié)構(gòu)有限元計(jì)算實(shí)例,77,1結(jié)構(gòu)的離散化與編號,78,2各個單元的矩陣描述,結(jié)構(gòu)包括有斜桿，所以必須在總體坐標(biāo)下對節(jié)點(diǎn)位移進(jìn)行表達(dá)，所推導(dǎo)的單元剛度矩陣也要進(jìn)行變換,79,3建立整體剛度方程,1.將所得到的各個單元剛度矩陣按節(jié)點(diǎn)編號進(jìn)行組裝，可以形成整體剛度矩陣;2.同時將所有節(jié)點(diǎn)載荷也進(jìn)行組裝。,80,4邊界條件的處理及剛度方程求解,81,82,5各單元應(yīng)力的計(jì)算,83,6支反力的計(jì)算,將節(jié)點(diǎn)位移的結(jié)果代入整體剛度方程中,84,對于單元2：取i=1，j=2，則，故,對于單元1：取i=3，j=1，則c=1，s=0，故,對于單元3：取i=2，j=3，則c=0，s=1，故,例2,85