第三章桿梁結構的有限元原理

1、1桿、梁單元概述討論桿梁單元和由它們組成的平面和空間桿梁結構系統(tǒng).o從構造上來說其長度遠大于其截面尺寸的一維構件o承受軸力或扭矩的桿件成為桿o桿梁問題都有精確解（且是唯一的）o承受橫向力和彎矩的桿件稱為梁o平面桁架 平面剛架 連續(xù)梁 空間剛架 空間桁架等o承受軸力或扭矩的桿件稱為桿o將承受橫向力和彎矩的桿件稱為梁o變截面桿和彎曲桿件2本章主要內容33.1有限元分析的完整過程E1=E2=2E7PaA1=A2=2cm2l1=l2=10cmP3為10N作用下二桿結構的變形。4E1、A1E2、A2說明：說明：u1、u2、u2分別表示分別表示節(jié)點節(jié)點1、2、3的水平位移的水平位移1）用標準化的分段小單元

2、來逼近原結構2）尋找能夠滿足位移邊界條件的許可位移場3）基于位移場的最小勢能原理來求解 基本變量為：節(jié)點位移內部各點位移應變應力(1)(3)(2)5問題的解題思路完整的求解過程1）結構離散化 該構件由兩根桿件做成，因此可以自然離散成2個桿單元。假定以這類單元位移的特征為兩個端點位移，就這兩個離散單元給出節(jié)點編號和單元編號。 單元1：i=1，j=2 單元2：i=2，j=36 單元位移模式：u(x)=a0+a1x 單元節(jié)點條件：u(0)=u1， u(l)=u2 將式（b）代入式（a）， 從而得01,1,2jiieuuauaijl72）單元分析（a）（b）回代得 寫成矩陣形式為01( )1jiiei

3、jeeiijju xaa xuuuxlxxuullN uN uiiujujeuNNuNuu其中Ni，Nj是形函數(shù)。形函數(shù)矩陣形函數(shù)矩陣8說明：說明：u表示位移列陣表示位移列陣 ue表示單元位移表示單元位移根據(jù)幾何方程可得單元應變的表達單元應變寫成矩陣形式為簡記為211d1dxjieuuuauuxll11 1iiiujuejjuuNNuuleBu幾何函數(shù)矩陣或者是應變轉換矩陣幾何函數(shù)矩陣或者是應變轉換矩陣9根據(jù)物理方程可得單元應力的表達單元應力寫成矩陣形式為簡記為ddxxjieuEEEuuxleSu單元應力矩陣或者是應力轉換矩陣單元應力矩陣或者是應力轉換矩陣1 1iiiujuejjuuEE NN

4、uul節(jié)點位移列陣10單元e勢能的表達 1 122T1 1220TT1 1220TT1 1221d21d21d212eeeeeeTeeleeeleeeeee eUWPuPuBuSuAxPuPuuB EBu AxPuPuuB EBu A lPuPu 11說明說明 積分域，積分域， P1、P2、分別表示作用單元分別表示作用單元e上的節(jié)點在上的節(jié)點在u1、u2的力的力eeSueeBueeUW分別為形變勢能和外力勢能寫成矩陣形式為TT1 1221112122211121222TT1211111 1121212eeee ee eeeeeeeeeeeeeeeeuB EBu A lPuPuuuuuEA lP

5、PuullEAEAuulluuPPuuEAEAllKF 單元單元e剛度矩陣剛度矩陣單元單元e節(jié)點力列陣節(jié)點力列陣12在得到各個單元的勢能表達式后，需要進行離散單元的裝配，以求出整個系統(tǒng)的總勢能，對于該系統(tǒng)，總勢能包括兩個單元部分 121T112T221T12T2112211222211121233112233221122111111123112110022012equ K uuK uP uP uEAEAEAEAuuuulllluuRuuFuuuuEAEAEAEAllllEAEAllEAEAuuull 1122213212222332200uuEAEAuRFulluuEAEAll133）離散單元

6、的裝配處理邊界條件是獲取可能位移場，將左端的約束條件，即u1=0代入上式可以得到簡化的勢能表達式 121T112T221T12T21221222223322332212102eu K uuK uP uP uEAEAEAuullluuFuuEAEAll 144）邊界條件的處理由于上式是基于許可位移場的表達的系統(tǒng)勢能，這是由全部節(jié)點位移分段所插值出的位移場為全場許位移場，且基本未知量為節(jié)點位移，根據(jù)最小勢能原理（即針對未知位移求一階變分）有12212222233220EAEAEAullluFEAEAll155）建立剛度方程120e 將結構參數(shù)和外載荷代入上式有求解得（單位m）22222223322

7、30EAEAulluFEAEAll233102E41110uu232.5E47.5E4uu166）求解節(jié)點位移7）計算單元應變1111211 12.5E3iiujujuNNuuul22222311 15E3iiujujuNNuuul178）計算單元應力111121 10.05MpaiiujujuE NNuuEul2222231 10.1MpaiiujujuE NNuuEul18對于單元勢能的表達，對其取極值有具體地對于單元1，有其中R1是節(jié)點1的支反力，P2是單元1的節(jié)點2所受的力，即單元2對該節(jié)點的作用力，將前面求得的節(jié)點位移代入上式可得支反力大小。eeeK uF1111221111uREA

8、uPl199）計算支反力以上是一個簡單結構有限元方法求解得完整過程，對于復雜結構，其求解過程完全相同，由于每一個步驟都具備標準化和規(guī)范性的特征，所以可以在計算機上編程而自動實現(xiàn)。討論1：對于一個單元的勢能取極值，所得到的方程為節(jié)點的位移和節(jié)點力之間的關系，也稱為單元的平衡關系，由此可以求出每一個單元所受的節(jié)點力。20討論2：由前面的步驟，我們也可以直接將各個單元的剛度矩陣按照節(jié)點編號的對應位置來進行裝配，即在未處理邊界條件之前，先形成整體剛度矩陣。其物理意義是，表示在未處理邊界條件前的基于節(jié)點描述的總體平衡關系。在對該方程進行位移邊界條件的處理后就可以求解，這樣與先處理邊界條件再求系統(tǒng)勢能的最

9、小值所獲得的方程完全相同。KUF21 小結 有限元分析的基本步驟及表達式1、物體幾何區(qū)域的離散化2、單元的研究（所有力學信息都用節(jié)點位移）來表達3、裝配集成4、邊界條件的處理并求解節(jié)點位移5、支反力的求取以及其它力學量（應力、應變及位移三大物理量）的計算eBueFNFeSueeeK uF22有限元分析的基本步驟及表達式2324一拉壓桿單元q (x)FiFjuiujyxij圖2.1 拉壓桿單元示意圖設桿單元長度為 ，橫截面面積為 ，單元材料的彈性模量為 ，在局部坐標系中桿端荷載分別為 和 ，桿端位移分別為 和 ，單元上的軸向分布荷載為 。 lAEiFjFiuju( )q x3.2 局部坐標下的桿

10、單元分析25用結點位移表示單元上任意截面的位移。對拉壓桿單元，可以取其位移為一次多項式，即 由位移的邊界條件：可得系數(shù) 、 為： 這樣，截面任意一點的位移 為：用矩陣表示為： 其中 u( )u xabx(0)iuu( )ju luabiaujiuubl( )(1)ijxxu xuulliiijjijjuuN uN uNNuN1ixNl =, TjijxNuul （3-1） （3-2） 單元位移模式。26 根據(jù)材料力學中應變的定義，有 這里 為應變矩陣。由虎克定律，其應力為： 11ijdudBBdxdxll NB11ll BEEBS （3-4） （3-3） 進行應力、應變分析EEEll SB其中

11、27利用虛位移原理求單元剛度矩陣，設桿端i、j分別產生虛位移 、 ，則由此引起的桿軸截面任意位置的虛位移為：對應的虛應變?yōu)椋焊鶕?jù)虛位移原理虛功方程（力乘以虛位移得虛功、外力虛功等于變形虛功），有： 將上式整理得： iujuTiiuuuNN B 000( )lTdllTTWq xdxWAdxEAdxFN BB 外變00( )TllTTTdq xdxEA dxFNBB （3-5） （3-6） 求單元剛度矩陣28式中 ：為局部坐標系下單元結點荷載矩陣。記 則可以得到拉壓桿單元的單元剛度方程為： 這里 為局部坐標系下的單元剛度矩陣 ， 為局部坐標系下等效結點荷載矩陣，但值得指出的是：分布荷載 中可以包

12、含集中荷載。根據(jù)定義，可以進一步求得單元剛度矩陣為： TdijFFF0( )lTEq xdxFN0lTEA dxkBBdEFFFk kEF( )q x1111EAlk （3-10） （3-7） （3-8） （3-9） 等效結點荷載29二扭轉桿單元m (x)MiMjyxijji圖2 扭轉桿單元示意圖設扭轉桿單元的長度為 ，截面慣性矩為 ，剪切模量為 ，桿端扭矩分別為 、 ，桿端扭轉角分別為 、 ，單元上的分布荷載集度為 ，則任意截面的扭轉角為 位移函數(shù)求得如（一） lIGiMjMij( )m x(1)ijxxll N式中 為局部坐標系下扭轉桿單元的結點位移矩陣Tij （3-11） 30由材料力學

13、可知，截面扭矩為：式中： 我們利用極小勢能原理來進行單元分析，桿單元的勢能用泛函表示為： dMGIGIdxB11ddxll NB00001()( )21( )2llTTpdllTTTddMdxm xdxdxGI dxm xdx FBBNF內力勢能外力勢能其中 為局部坐標系下扭轉桿單元的結點集中荷載矩陣TdijMMF （3-12） 31由極小勢能原理，取上述泛函的變分 ，可得： 或者寫為：設： 可得扭轉桿單元的單元剛度方程為： 可以看到，其形式與拉壓桿單元的單元剛度方程完全一致。同樣，由上式可以進一步求得其局部坐標系下得單元剛度矩陣為： 0p 00( )llTTTdGI dxm xdxBBNF0

14、0()( )llTTdGI dxm xdxBBNF0lTGI dxkBB0( )lTEm xdxFNdEFFk 1111GIlk （3-13a） （3-13b) （3-14、3-15） （3-16） （3-17） 等效結點荷載32三只計彎曲的桿單元q (x)MiMjFyiFyjvivjyxm (x)ijji設桿單元的長度為 ，截面慣性矩為 ，彈性模量為 ，桿端集中剪力為 、 ，桿端集中彎矩分別為 、 ，桿端橫向位移為 、 ，桿端扭轉角分別為 、 ，在單元上分布有荷載集度為 的豎向分布荷載和集度為 的分布力偶。lIEyiFyjFiMjMivjvij( )q x( )m x33根據(jù)梁的平截面假定可

15、知平面純彎梁單元的軸向應變?yōu)椋哼@里利用平截面假設（這里利用平截面假設（變形后橫截面仍保持平面，與縱線正交）如圖： 22223322dxvdEyEdxvdydxvdEIQdxvdEIMdxdv應變和應力公式：彎曲公式： 材料力學基礎知識34取撓曲線方程為 的三次多項式，即單元上任意一點的撓度為： 根據(jù)單元的位移邊界條件： 時： 時： TiijjvvTdyiiyjjFMFMFx23vabxcxdx0 x iidvvvdxxl jjdvvvdx222323212312121ijiijjiijjavbvcvvlllldvvllll 可以得到式中的待定系數(shù)結點位移矩陣和結點荷載矩陣分別為 （3-18）

16、35將系數(shù)a、b、c、d代入式，并將撓曲線方程用矩陣形式表示為： 式中 為形函數(shù)矩陣，其中： 上式為平面彎曲單元的形函數(shù)。 2322323210000100323112121iijjvvxxxvllllllllN1234NNNNN232122322323342323221 (1)32 xxxxNNxllllxxxxNNllll （3-19） （3-20） 36 根據(jù)式（2-19）確定的單元位移場，可得單元上某一點得曲率為： 截面的彎矩為： 這里： 為平面彎曲桿單元的應變矩陣。 根據(jù)虛位移原理。有： 2222d vddxdxNBTTMEIEIEIBB222322326124661226dxxxx

17、dxllllllll NB000( )( )llTdlTTdWq xdxm xdxWdxEI dxNNFBB外變 （3-21） 37 則平面彎曲桿單元的單元剛度方程為： 其中的單元剛度矩陣可由式（2-23）求得為： 00( )( )llEdq xdxm xdxdxNFN0lTEI dxkBBdEFFk 2232212612664621261266264llllllEIlllllllk 記：（3-23）（3-22）等效結點荷載（3-24）（3-25）38四平面一般桿單元（考慮拉伸、彎曲、不考慮扭轉）MiMjFxiFyiFxjFyjuiviujvjyxji桿單元的長度為 ，截面面積為 ，截面慣性矩

18、為 ，彈性模量為 ，單元的 、 端各有三個力為 、 、和 、 、 ，其對應的位移為 、 、 和 、 、 ，建立如圖所示的局部坐標系，各物理量的正向如圖中所標。lAIEijxiFyiFiMxjFyjFjMiuivijujvj39 TiiijjjuvuvTxiyiixjyjjFFMFFMF設單元上沒有荷載作用，首先考慮軸向力的作用，由于桿端軸力 、 只引起桿端軸向位移 、 ，根據(jù)拉壓桿單元的單元剛度方程，有： xiFxjFiujuxiijEAEAFuullxjijEAEAFuull 則結點位移矩陣和結點荷載矩陣分別為： （3-26）（3-27）40 其次，桿端彎矩 、 和桿端剪力 、 只與桿端的轉

19、角位移 、 和桿端的橫向位移 、 有關系，根據(jù)只計彎曲桿單元的單元剛度方程（注意，由于不考慮單元上的荷載作用，故方程式中的等效結點荷載 等于零）可得：iMjMyiFyjFijivjvEF3232126126yiiijjEIEIEIEIFvvllll3232126126yjiijjEIEIEIEIFvvllll 226462iiijjEIEIEIEIMvvllll226264jiijjEIEIEIEIMvvllll結構力學相關知識41這樣，上述表達式合并在一起，寫成矩陣形式如下： 323222323222000012612600646200000012612600626400ixiiyiiijx

20、jjyjjjEAEAllEIEIEIEIuFllllvFEIEIEIEIMlllluFEAEAllvFEIEIEIEIMllllEIEIEIEIllllFk 可以將上式簡寫為：（3-28）（3-29）42其中單元剛度矩陣323222323222000012612600646200000012612600626400eEAEAllEIEIEIEIllllEIEIEIEIllllkEAEAllEIEIEIEIllllEIEIEIEIllll（3-30）43五 空間受力桿件單元（考慮扭轉、拉伸、彎曲）MxiMyjFxiFyiFxjFyjyxziMziFziMyiMxjMzjFzjj對空間桿件單元，除

21、了桿端力和結點位移數(shù)目較平面單元多外，其分析方法與平面桿單元類似（包含拉伸、扭轉、兩個方向彎曲）x設局部坐標系的 軸為單元的形心主軸，橫截面的兩個主軸分別為 軸和 軸（如圖所示）。設桿橫截面面積為 ，桿單元長度為 ，在 平面內抗彎剛度為 ，在 平面內的抗彎剛度為 ，桿件的抗扭剛度為 。 xyzAlx zyEIx yzEIEI44空間剛架有6個位移分量和6個結點力分量，設局部坐標系下它們分別為TiiixiyizijjjxjyjzjuvwuvwTxiyizixiyizixjyjzjxjyjzjFFFMMMFFFMMMFo純軸向拉壓,1111iju uEAklo純扭轉,1111ijGlkI 45ox

22、oy面內彎曲22z,32212612664621261266264vizivjzjllllllEIkllllllloxoz面內彎曲22,32212612664621261266264yiwjyjwiyllEIllllklllllll46323232320000000000000000000000000000000000000126126126126zzzzxiyiyyyyzixiyizixjyjzjxjyjzjEAEAllEIEIEIEIFllllFEIEIEIEIllllFGIGIllMMMFFFMMM222232323232220000000000000000000000000000000

23、000000000000000000000000000006462646212612612612662646yyyyzzzzzzzzyyyyyyyyEIEIEIEIllllEIEIEIEIllllEAEAllEIEIEIEIllllEIEIEIEIllllGIGIllEIEIEIEIllll220000000264iiixiyizijjjxjyjzjzzzzuvwuvwEIEIEIEIllll （3-31）47其中的單元剛度矩陣可寫為Fk將式（2-31）寫成矩陣的形式有323232322222300000000000000000000000000000000000000000000000000

24、000000000000000012612612612664626462126zzzzyyyyyyyyzzzzzEAEAllEIEIEIEIllllEIEIEIEIllllGIGIllEIEIEIEIllllEIEIEIEIllllEAEAllEIEl232323222220000000000000000000000000000000000000012612612662646264zzzyyyyyyyyzzzzIEIEIlllEIEIEIEIllllGIGIllEIEIEIEIllllEIEIEIEIllll（3-32）（3-33）48六單元剛度矩陣的性質 單元剛度矩陣 為對稱矩陣，其元素 單

25、元剛度矩陣 中的每個元素代表單位桿端位移引起的桿端力。其中的任意元素 的物理意義是第 個桿端位移分量等于1（其余位移分量等于0）時，所引起的第 個桿端力的分量值。k()ijjikkijkijkji49 一般單元的單元剛度矩陣 是奇異矩陣，它的元素組成的行列式等于零，即 。根據(jù)奇異矩陣的性質， 沒有逆矩陣。也就是說，如果給定桿端位移 ，根據(jù)（2-29）或（2-31）式可以求出桿端力 的惟一解，但反過來，如果已知桿端力 ，則不能根據(jù) 來確定桿端位移 的惟一解。因為即使在桿端力已知的情況下，由于單元兩端無任何約束，因此除出桿端自身變形外，還可以發(fā)生任意的剛體位移。舉例來說，如果物體處于靜止狀態(tài)，我們

26、可以說其處于平衡狀態(tài)，但反過來，如果物體處于平衡狀態(tài)，則我們不能說其一定處于靜止。 k0kkF1()dkFF50 單元剛度矩陣 具有分塊的性質，即可以用子矩陣表示 。用虛線把 分為四個子矩陣，把 和 各分為兩個子矩陣，因此，又可以寫為： 這里： 或 或 或 或 用子矩陣形式表示單元剛度矩陣和單元剛度方程，可以使其表達的物理意義更加明顯。在單元剛度矩陣 中，其任意子矩陣 表示桿端力 和桿端位移 之間的關系。 kkkFiiiijijjijjjFkkFkkTixiyiiFFMFTixiyizixiyiziFFFMMMFTjxjyjjFFMFTjxjyjzjxjyjzjFFFMMMFTiiiiuvTi

27、iiixiyiziuvwTjjjjuvTjjjjxjyjzjuvwkrskrFs（3-34）513.3桿系結構的整體分析（整體坐標系）一平面問題坐標轉換矩陣圖 平面問題桿端力轉換示意圖ijoyxMiMjFxiFyiFxjFyjyxFxiFxjFyiFyjMjMiijoyxyx一般情況下，在進行單元分析時是在局部坐標下完成的。對于某一單元而，如果局部坐標系與整體坐標系不一致，則有單元分析的物理量必須通過坐標轉換到整體坐標系中，然后再進行整體坐標系下的分析52 這里 表示由 軸到 軸的角，角度轉動的正負由右手定則確定，本書中以順時針方向轉動為正。在兩個坐標系中，力偶分量保持不變，即有： 同理，對于

28、 端的桿端力，有： c o ss ins inc o sx ix iy iy ix iy iFFFFFF xxiiMMjcossinsincosxjxjyjyjxjyjjjFFFFFFMM 根據(jù)力的投影定理，將整體坐標下的桿端力分別投影到局部坐標下，有如下關系（3-35）（3-36）（3-36）53將這些式子用矩陣形式可表示為：cossin0000sincos0000001000000cossin0000sincos0000001xixiyiyiiixjxjyjyjjjFFFFMMFFFFMM FTFcossin0000sincos0000001000000cossin0000sincos00

29、00001T上式可以簡寫成：這即為兩種坐標系下單元桿端力的坐標變換式。其坐標轉換矩陣為：（3-38）（3-39）（3-40）54從坐標轉換矩陣 的表達式可以看出， 為正交矩陣，其逆矩陣等于其轉置矩陣，即有： 并且有： 式中 為單位矩陣。 同樣的推導，可以得到兩種坐標系下的桿端位移之間的轉換關系為： 這里 和 分別為局部坐標系和整體坐標系下的桿端位移矩陣， 為前面介紹的轉換矩陣。 TT1TTT1TTIIT T整體坐標下的桿端力為( )TexiyiixjyjjFFFMFFM局部坐標下的桿端力為( )TexiyiixjyjjFFFMFFM（3-42）（3-41）（3-43）（3-44）（3-45）5

30、5因此可得：上式兩邊同乘以 ，可以得到： 設 可得： 上式即為整體坐標系下的單元剛度方程。 T Fk T 1T1TFTk T Tk T TkTk TFk （3-47）（3-46）剛度矩陣轉換）56二 空間問題的坐標變換（空間問題） oFxiFyiFxjFyjyxxyzijFzjFziz考慮結點i在局部坐標下 的桿端力 與在整體坐標系 的桿端力 的關系 oxyz,xiyiziFFFoxyz,xiyiziFFF57設 軸與x、y、z軸的方向余弦分別為： ， ， 則將桿端力 、 、 向 軸投影，可以求得桿端力 ，即：同理可以求得： xcos( , )xxlx xcos( , )xylx ycos(

31、, )xzlx zxiFyiFziFxxiFxixi xxyi xyzi xzFF lF lF lyixi yxyi yyzi yzFF lF lF lzixi zxyi zyzi zzFF lF lF l58用矩陣形式可以表示為：上式即為結點i的桿端力在局部坐標系和整體坐標系下的轉換關系，其中的矩陣稱為關系矩陣。與上面的推導類似，同樣可以推出以 、 、 表示 、 、 ，以及對于結點j的相對應的轉換關系，其中轉換關系矩陣都是 。 xixxxyxzxiyiyxyyyzyizizxzyzzziFlllFFlllFFlllFxxxyxzyxyyyzzxzyzzlllllllllxiMyiMziMxi

32、MyiMziM 綜上所述，整體坐標系下單元桿端力矩陣與局部坐標系下單元桿端力矩陣具有如下的關系表達式：FT F（3-48）59其中的 為：稱為空間坐標系的單元轉換矩陣，它是一個正交矩陣，即：對于桿端位移，同樣可推導出在兩種坐標系中的轉換關系： 這樣，可得空間桿件單元在整體坐標系中單元剛度方程為： 其中表示空間單元在整體坐標系中的單元剛度矩陣。TT0000000000001TTTT Fk TkTk T（3-49）（3-50）（3-51）（3-52） （3-53）60三桿系結構的整體分析對桿系結構進行單元分析，僅僅是有限元分析中的第一步。我們的目的是要對整個結構進行分析，研究結構的整體性能。因此，

33、在對結構的各單元分析完成后，必須將單元分析的結果進行整合，對結構進行整體分析。整體分析的過程實際上是如何將單元分析的結果進行有效組合，建立整體剛度方程并求解結點位移的過程。根據(jù)對結點位移的編碼方式，可以采用“先處理法”和“后處理法”來建立整體剛度方程。 611后處理法所謂后處理法，就是由單元剛度矩陣形成整體剛度矩陣，建立剛度方程后再引入支承條件，進而求解結點位移的方法。運用這種方法時，假設所有結點位移均為未知量，按照順序統(tǒng)一進行編碼，如圖所示的平面桿件單元。1 (1 2 3)2 (4 5 6)3 (7 8 9)4 (10 11 12)213xyO62結點位移矩陣為：結點荷載矩陣為：求出各單元剛

34、度方程后，根據(jù)平衡條件和位移連續(xù)條件，可以建立整個結構的位移法方程：1234111222333444TTuvuvuvuv1234111222333444TTxyxyxyxyFFMFFMFFMFFMFFFFF11223344iiijjijjiiijjijjiiijjijjkkFkkkkFkkkkFkkF00000063 簡 寫成： 這里 為結構的整體剛度矩陣，有： FK K11121314212223243132333441424344KKKKKKKKKKKKKKKKK注意，在建立方程的過程中，我們假設所有結點都有位移。因此整個結構在外力作用下，除了發(fā)生彈性變形外，還可能發(fā)生剛體平動位移，這樣各

35、結點位移不能唯一確定。這說明整體剛度矩陣為一奇異矩陣，不能求逆矩陣，即根據(jù)整體剛度方程可得到無窮多個解。 64實際上，在圖所示剛架中，結點1和結點4均為固定端，其三個位移分量均為0，即有：這樣，將上述支承條件引入到方程中，對整體剛度方程進行修改，可得：1114440uvuv11121314121222324223132333433414243444KKKKKKKKKKKKKKKKFFFF0065對上述方程進行化簡，可以得到兩組方程： 和 222232332333KKKKFF12214334KKFF00這樣，利用第1式可以求得結點位移 和 ，再根據(jù)第2式可以求得支座反力 和 。 231F4F66

36、2先（前）處理法所謂先處理法，就是先引入支承條件，根據(jù)支承條件僅對未知的自由結點位移分量編號，得到的位移矩陣中不包含已知的約束位移分量，即可以直接得到方程求解自由結點位移分量。67如圖所示的平面剛架結構ABCD，由于在A處和D處均為固定端，其位移為0，故位移編碼均為0，在C處為鉸接，故BC桿在C端的角位移與DC桿在C端的角位移不相同，因此在C處編兩個結點3和4，但結點3和4的橫向位移和豎向位移相同，故采用相同的編號，各結點位移編碼如圖所示。圖先處理法位移編碼示意圖1 (0 0 0)2 (1 2 3)3 (4 5 6)4 (4 5 7)ABCD5 (0 0 0)213xyO683.4等效結點荷載

37、與邊界條件的處理非結點等效荷載和邊界條件的處理是有限元分析中必須考慮的的兩個重要方面。由于只考慮結點荷載，因此必須將作用于單元上的非結點荷載轉換到節(jié)點上。有限元的結點荷載來自兩部分（1）作用于結點處的集中力、集中力偶，前面多次提到的 直接疊加到結點上即可（整體坐標下）（2）非結點荷載首先須在局部坐標下等效到結點荷載 ， 然后再轉到整體坐標系下的 ( )edF( )eEF( )eEF691非結點荷載的處理在前面的分析中，我們已經(jīng)介紹了求等效結點荷載的方法，如3-7式、3-15式、3-22式分別可用來求不同情況下的等效結點荷載。此外，可以這樣來考慮：第一步，在局部坐標系下求單元的固端力 。對于某個

38、單元，我們假定單元的兩端均固定，然后根據(jù)靜力平衡求得固定端的反力。第二步，根據(jù)單元固端力求單元的等效結點荷載 。根據(jù)局部坐標系與整體坐標系單元桿端力的變換式，固端內力在兩種坐標系下的變換形式可以寫成：因此，整體坐標系下的等效結點荷載矩陣可以由下式計算： fFTffFTFEFTEf FTF0( )lTEq xdxFN0( )lTEm xdxFN00( )( )llEdq xdxm xdxdxNFN（3-7）（3-15）（3-22）（拉壓桿）純彎桿扭轉桿70平面剛架單元固端力（遵行結構力學）71722邊界條件的處理 （1）鉸結點）鉸結點在桿系結構中，除了剛性結點外，通常會遇到一些桿件通過鉸結點與其

39、它桿件聯(lián)結，如下圖所示桿件系統(tǒng)，有4根桿件匯交于D點，其中BD桿在D端通過鉸支座與其它桿件鉸接，其余3根桿為剛性接觸。對于這樣的鉸結點，具有如下的性質：鉸結點上各桿具有相同的線位移，但截面的轉角位移不相同；結點上具有鉸接桿端不承受彎矩作用。如下圖所示結構中，BD桿在D端的桿端彎矩為0，只有CD、ED、GD桿在結點D上與外彎矩保持平衡。73對于這樣的結點，我們在對其進行單元劃分時，通?？紤]在D處設置2個結點。按照先處理法，對圖示結構進行位移編碼，如圖2b所示。ABCDEHGF1 (0 0 1)2 (2 3 4)3 (5 6 7)4 (5 6 8)5 (0 0 0)6 (9 10 11)7 (0

40、0 0)8 (12 13 14)9 (15 16 17)圖鉸結點的處理示意圖 ANSYS通過結點耦合實現(xiàn)742彈性支承點彈性支承點在實際工程中，有時會遇到彈性支承的情況（如圖），這時一般將彈性支座看作是在結構約束點沿約束方向的一個彈簧，彈簧的剛度系數(shù)為 ， 在數(shù)值上等于使彈簧支座沿約束方向產生單位位移時所需施加的力。 kk12i-1ii+1nkANSYS引入彈簧單元即可75具體做法可以歸結為：先解除彈性支承點約束，在i處給一個結點號，形成總剛度矩陣，然后在總剛度矩陣中將第i行的主元素加上彈性支承的剛度系數(shù)，此時第行變?yōu)椋?23iiiiiinKKKKkK以上的分析也適用與角位移為彈性約束的情況。

41、若有多個彈性支座，可同時引入，即只需將相應的主對角線元素加上相應的彈性剛度系數(shù)即可。76例1 3.5桿系結構有限元計算實例1 結構的離散化與編號 772 各個單元的矩陣描述 結構包括有斜桿，所以必須在總體坐標下對節(jié)點位移進行表達，所推導的單元剛度矩陣也要進行變換22222222-eeeeTeeeccsccscsscssE AkT k Tlccsccscsscss783 建立整體剛度方程 1.將所得到的各個單元剛度矩陣按節(jié)點編號進行組裝，可以形成整體剛度矩陣;2.同時將所有節(jié)點載荷也進行組裝。 794 邊界條件的處理及剛度方程求解 80815 各單元應力的計算 826 支反力的計算 將節(jié)點位移的

42、結果代入整體剛度方程中83對于單元2：取i=1，j=2，則 ， 故 x 1 2 3 1 1 p 45o y 對于單元1：取i=3，j=1，則c=1，s=0，故 111333111111131110-1 00000-1 0100000kkE Aklkk22,22sc22211122222212221-1 -11-111-1-111-12 21-1 -11kkE Akkkl對于單元3：取i=2，j=3，則c=0，s=1，故 333333233333232200-0-001-0-10-000-0-101kkE Aklkk84例 2整體編號，對號入座得總剛1 22111121322 332122231

43、31 331323311111102 22 22 22 21111002 22 22 22 21111002 22 22 22 211111012 22 22 22 2100010000101kkkKkkkkkk8511112222333311111102 22 22 22 21111002 22 22 22 21111002 22 22 22 211111012 22 22 22 2100010000101xyxyxyFuFvFuFvFuFv11112233111102 22 22 211102 22 22 2111112 22 22 20001xyyyFuFvFvFv 86消除邊界條件后有 求解左式即可得到相關位移、應力、應變等等87例3其中E1=3e4MPa,E2=3e4MPa解(1)單元劃分，建立局部坐標系和整體坐標系，并對數(shù)據(jù)進行整理，對單元和結點編號4444234300 10 ,100 10 ,61230 10 ,12 10EAEIllEIEIll88(2) 求局部坐標系下的單元剛度矩陣，由于單元1和2的尺寸完全一樣，因此其單元剛度矩陣一樣，為(1)(2)4300003000001230012300