有限元分析理論基本

1、有限元分析概念有限元法：把求解區(qū)域看作由許多小的在節(jié)點處相互連接的單元（子域）所構(gòu)成，其模型給出基本方程的分片（子域）近似解，由于單 元（子域）可以被分割成各種形狀和大小不同的尺寸，所以它能很好地 適應(yīng)復雜的幾何形狀、復雜的材料特性和復雜的邊界條件有限元模型：它是真實系統(tǒng)理想化的數(shù)學抽象。由一些簡單形狀的 單元組成，單元之間通過節(jié)點連接，并承受一定載荷。有限元分析：是利用數(shù)學近似的方法對真實物理系統(tǒng)（幾何和載荷 工況）進行模擬。并利用簡單而又相互作用的元素，即單元，就可以用 有限數(shù)量的未知量去逼近無限未知量的真實系統(tǒng)。線彈性有限元是以理想彈性體為研究對象的，所考慮的變形建立在 小變形假設(shè)的基礎(chǔ)

2、上。在這類問題中，材料的應(yīng)力與應(yīng)變呈線性關(guān)系， 滿足廣義胡克定律；應(yīng)力與應(yīng)變也是線性關(guān)系，線彈性問題可歸結(jié)為求 解線性方程問題，所以只需要較少的計算時間。如果采用高效的代數(shù)方 程組求解方法，也有助于降低有限元分析的時間。線彈性有限元一般包括線彈性靜力學分析與線彈性動力學分析兩 方面。非線性問題與線彈性問題的區(qū)別：1）非線性問題的方程是非線性的，一般需要迭代求解；2）非線性問題不能采用疊加原理；3）非線性問題不總有一致解，有時甚至沒有解。有限元求解非線性問題可分為以下三類:1）材料非線性問題材料的應(yīng)力和應(yīng)變是非線性的，但應(yīng)力與應(yīng)變卻很微小，此時 應(yīng)變與位移呈線性關(guān)系，這類問題屬于材料的非線性問題

3、。由于從 理論上還不能提供能普遍接受的本構(gòu)關(guān)系，所以，一般材料的應(yīng)力 與應(yīng)變之間的非線性關(guān)系要基于試驗數(shù)據(jù)，有時非線性材料特性可 用數(shù)學模型進行模擬，盡管這些模型總有他們的局限性。在工程實 際中較為重要的材料非線性問題有：非線性彈性（包括分段線彈性） 彈塑性、粘塑性及蠕變等。2）幾何非線性問題幾何非線性問題是由于位移之間存在非線性關(guān)系引起的。當物體的位移較大時，應(yīng)變與位移的關(guān)系是非線性關(guān)系。研究 這類問題一般都是假定材料的應(yīng)力和應(yīng)變呈線性關(guān)系。它包括大位 移大應(yīng)變及大位移小應(yīng)變問題。如結(jié)構(gòu)的彈性屈曲問題屬于大位移 小應(yīng)變問題，橡膠部件形成過程為大應(yīng)變問題。3）非線性邊界問題在加工、密封、撞擊等

4、問題中，接觸和摩擦的作用不可忽視， 接觸邊界屬于高度非線性邊界。平時遇到的一些接觸問題，如齒輪傳動、沖壓成型、軋制成型、 橡膠減振器、緊配合裝配等，當一個結(jié)構(gòu)與另一個結(jié)構(gòu)或外部邊界 相接觸時通常要考慮非線性邊界條件。實際的非線性可能同時出現(xiàn)上述兩種或三種非線性問題。有限元理論基礎(chǔ)有限元方法的基礎(chǔ)是變分原理和加權(quán)余量法，其基本求解思想是把計算域劃分為有限個互不重疊的單元，在每個單元內(nèi)，選擇一些合適的節(jié)點作為求解函數(shù)的插值點，將微分方程中的變量改寫成由各變量或其 導數(shù)的節(jié)點值與所選用的插值函數(shù)組成的線性表達式，借助于變分原理或加權(quán)余量法，將微分方程離散求解。采用不同的權(quán)函數(shù)和插值函數(shù)形式，便構(gòu)成不

5、同的有限元方法。1.加權(quán)余量法：是指采用使余量的加權(quán)函數(shù)為零求得微分方程近似解的方法稱為 加權(quán)余量法。(Weighted residual method WRM )是一種直接從 所需求解的微分方程及邊界條件出發(fā)，尋求邊值問題近似解的數(shù)學方 法。加權(quán)余量法是求解微分方程近似解的一種有效的方法。設(shè)問題的控制微分方程為：在 V 域內(nèi)L(u)f 0(5.1.1)在 S 邊界上B(u) g 0(5.1.2)式中：L、B分別為微分方程和邊界條件中的微分算子；f、g為與未知函數(shù)u無關(guān)的已知函數(shù)域值；為問題待求的未知函數(shù)當小J川加權(quán)余“法以解日寸. 廿先在水種或JL-M；立一個這屆I敬 衣. 一般共方如下相式

6、：力RgkCN = NC(5 1)一 上 1 /=1式中： C,付定冢牧. 也可稱為廣義生標；A;雙門 完衣泌才女買的右把也兀關(guān)的電函 故、由于力 一般只W將求自1效11的玨:似筋. 因此將式心I >)代入式(5 1 1)珀式(，L2)后將得不到滿尺. 音也；RZ -L(fi-f在V域內(nèi)4 -Bi fll-g在S線界上""顯然用 %反映亍送函數(shù)與H次解之間的諒妾.它T門分捌將做內(nèi)部和邊界余/ _若在歧'內(nèi)引入內(nèi)部權(quán)函數(shù)叫，在邊界S上312之界權(quán)國救 % 則可延立11個消除合同一的條件，一般可農(nóng)示為：| . WdV + | HRpdS = 0 (i = L2.L

7、 e)，* s不同的權(quán)函數(shù)”;和反映了不同的酒除余量的近則“ 從上 式可以得到求解付定系數(shù)炬蹤C的代 執(zhí)方程組,一經(jīng)解科付定系效.由式( I 3建P可得所需求解邊伍向虺的近似解,由于試函數(shù)H的不同. 余仲4和 q 可有如下三種餅況.依此加權(quán)余量。去可分為：1,內(nèi)部法試函數(shù)浙兄垃界條件. 也叩 火伊叱時消除余貨的條件成為 WkRiV- 0 (frlj.L E)(5.1.6)2.邊界法試函數(shù)滿尺控制方程，也即 尺=喇-六。此時消除余位的條件為： 開，/屋¥5 = 00 =L2,L .11)為里混合法試函數(shù)不調(diào)昆拄制方程和邊界條件.此時用式(5.5)未消除余旨混合法對于試函數(shù)的選取最方便，

8、但在相同精度條件下，工作量最大。對內(nèi)部法和邊界法必須使基函數(shù)事先滿足一定條件，這對復雜結(jié)構(gòu) 分析往往有一定困難，但試函數(shù)一經(jīng)建立，其工作量較小。無論采用何種方法，在建立試函數(shù)時均應(yīng)注意以下幾點：(1)試函數(shù)應(yīng)由完備函數(shù)集的子集構(gòu)成。已被采用過的試函數(shù)有暮級 數(shù)、三角級數(shù)、樣條函數(shù)、貝賽爾函數(shù)、切比雪夫和勒讓德多項式等等。(2)試函數(shù)應(yīng)具有直到比消除余量的加權(quán)積分表達式中最高階導數(shù)低一階的導數(shù)連續(xù)性。(3)試函數(shù)應(yīng)與問題的解析解或問題的特解相關(guān)聯(lián)。若計算問題具有 對稱性，應(yīng)充分利用它。顯然，任何獨立的完全函數(shù)集都可以作為權(quán)函數(shù)。按照對權(quán)函數(shù)的不同選擇得到不同的加權(quán)余量計算方法，主要有：配點法、子

9、域法、最 小二乘法、力矩法和伽遼金法。其中伽遼金法的精度最高。下而以內(nèi)涉法為例. 介紹按權(quán)的數(shù)分類時加權(quán)余呈的五種基本 方法口對內(nèi)部法來說，消除余量的筑一格式足二 0。= L2,L Ji)L 子域法(Subdonwiii Method)此法甘先將求解城丫劃分成ii個子域 匕, 在缶個吁 域內(nèi)令權(quán)函數(shù) 等于 L nr在子城之外;風權(quán)畫數(shù)為 卷. 也即.11俏內(nèi)) 代外) 如果在各個子歧里分別造型試留效，那名它的水番在形式上將 再似于有限元法2 配點法(Collocation Method)子域法是令余量在一個子城上的總和為安口而配點法是防 余量在指 定的n個 點上尋于4, 這些點制;為配比。此法

10、的權(quán)函數(shù)為：WI:=P-POQO:誑D送娛(優(yōu)拉克)網(wǎng)數(shù)，它的定義為：I 0 X 7= A.沆工止I CO A J /工-0=4r JQV A, <7 , bPf P分別代表求角咫內(nèi)任一點和配點3由于此法只在配點上保證余量為令.因此不冊要作積分計芥. 所以是筑肅單的方門權(quán)余過法3 .枇小二乘法(Least Square Method)本 法通過使在整個水廨墳上余苗的 平方和 取極 小 來 戲立消 除余量咱勺條件、若記余里平方和為I(C),艮p 7(C)= ：a；m= i寓R© .產(chǎn)，«爐則極1EL條件為：粵二J噲)，?”=o&R由此可見.本法權(quán)函數(shù)為：=-f

11、"L2工 g-i4 .伽遼金法(QalB依uMdhod)本法是使余胃與每一個 基 函數(shù)正 交. 也即以基函數(shù)作為 權(quán) 函 數(shù)%二%；? = I2L JI)當試函數(shù)期包會挺個元備國數(shù)算時，用本法必可求浮精確解i5,矩法(Method of Moment)本法與伽遼金法粕似* 也是用完備函數(shù)聯(lián)作權(quán)的數(shù)俱本法的 杈國 效與儂 遼金5去又力區(qū)利，它與試 函 數(shù)無關(guān)口 消除余登的條件是從筌開始的各階坦為零，因此對一維問題 %=婷 U =12 L網(wǎng)對二組問題 % =產(chǎn) C = i2L 其余類推這五種丞本方法在彳寸定系數(shù)足繆多”象做高階近似)時, 其精 成彼此相近口 但對低階近1以(n領(lǐng)小下次于U

12、下，后三種的餅皮 要高于前兩種L1基本方法羋例為說明上述基本視念. 以圖所示甘裁而號卉羿. 爻滿跨均布而 放作用,求懸劈蠟E的絳向位移;為例，說明基本方法的百1 口 囪示梁的控制方程為：y - - - Q(X = O )V a .若抽試國數(shù)為：9h«/+A J14產(chǎn)？+ 26/。(*)dy d3y n ,7 - -7-7 - 0 (x = 0L© 阪不昭驗證具滿足邊界條件，也即 為=稟 而拄制方程的內(nèi)部余量 號為:% 二 E/r(120v+ 24/)-q因此本問汜工隔內(nèi)郃法下面汾地用基本方法必仃求解.子域法解由于試色數(shù)僅一個待定用數(shù).因此只鏘用一個子球（孑于全城）即 可，捎

13、除余量的條件為：(120a + 24/)-7cA = 0由此可解詞:亡_ qS Ell代百出式可'碑川 7/& 42EI品點法懈同上所還.只需璉一"V'配點來咫工消除余圣的條件若令:技1L,c="可科：c二q於二型 114 巨衛(wèi) *57,B2若令則評C 44砒可見不JH的配點結(jié)果是不一樣的 最小二年法解此時消除余量的烝件為：I：號駕比=。貿(mào)匚|120工十即L q_|，|_&1120匯十四1_|心=0可得Ci.OlOj 7g ?=g W*應(yīng)如遼金法解此時，苑=d - d _+ 26戶F消除余量的條件為：N.Rrdx 0Jo11由此可得：6警蹲二

14、處咨L1矩法解由于只個傳定常數(shù),因此消除余管條件只需等次矩即可.此時顯然與子域法完全相同0本 例各方法的 精度比較本問儂的精確解由梁1 立核計算可得力：、_ / = 0.125" SEI EI由此可得，上述各方法對本例計鼻的誤差依次為：-333% ; 1.75% (222%) ; 13.9% ： 0.96%; -333%上面22.2 %為式竹村轉(zhuǎn)果2、虛功原理平衡方程和幾何方程的等效積分“弱”形式虛功原理包含虛位移原理和虛應(yīng)力原理，是虛位移原理和虛應(yīng)力原理 的總稱。他們都可以認為是與某些控制方程相等效的積分“弱”形式。虛 功原理：變形體中任意滿足平衡的力系在任意滿足協(xié)調(diào)條件的變形狀態(tài)

15、上 作的虛功等于零，即體系外力的虛功與內(nèi)力的虛功之和等于零。虛位移原理是平衡方程和力的邊界條件的等效積分的“弱”形式；虛應(yīng)力原理是幾何方程和位移邊界條件的等效積分“弱”形式。虛位移原理的力學意義：如果力系是平衡的，則它們在虛位移和虛應(yīng) 變上所作的功的總和為零。反之，如果力系在虛位移(及虛應(yīng)變)上所作 的功的和等于零，則它們一定滿足平衡方程。所以，虛位移原理表述了力系平衡的必要而充分條件。一般而言，虛位移原理不僅可以適用于線彈性 問題，而且可以用于非線性彈性及彈塑性等非線性問題。虛應(yīng)力原理的力學意義：如果位移是協(xié)調(diào)的，則虛應(yīng)力和虛邊界約束 反力在他們上面所作的功的總和為零。反之，如果上述虛力系在

16、他們上面 所作的功的和為零，則它們一定是滿足協(xié)調(diào)的。所以，虛應(yīng)力原理表述了 位移協(xié)調(diào)的必要而充分條件。虛應(yīng)力原理可以應(yīng)用于線彈性以及非線性彈性等不同的力學問題。但是必須指出，無論是虛位移原理還是虛應(yīng)力原理，他們所依賴的幾何方程 和平衡方程都是基于小變形理論的，他們不能直接應(yīng)用于基于大變形理論 的力學問題。3、最小總勢能法應(yīng)變能：作用在物體上的外載荷會引起物體變形，變形期間外力所做的功以彈性能的形式儲存在物體中，即為應(yīng)變能。由n個單元和m個節(jié)點組成的物體的總勢能為總應(yīng)變能和外力所 做功的差：mFUi1n(e)e 1i最小勢能原理：對于一個穩(wěn)定的系統(tǒng),相對于平衡位置發(fā)生的位移總會使系統(tǒng)的總勢能最小

17、，即:nm(e)F|Ji 0, i=1,2,3, ,nUi Ui e 1Ui i i有限元法的收斂性有限元法是一種數(shù)值分析方法，因此應(yīng)考慮收斂性問題。有限元法的收斂性是指：當網(wǎng)格逐漸加密時，有限元解答的序列收 斂到精確解；或者當單元尺寸固定時，每個單元的自由度數(shù)越多，有限 元的解答就越趨近于精確解。有限元的收斂條件 包括如下四個方面：1）單元內(nèi)，位移函數(shù)必須連續(xù)。多項式是單值連續(xù)函數(shù)，因此選 擇多項式作為位移函數(shù)，在單元內(nèi)的連續(xù)性能夠保證。2）在單元內(nèi)，位移函數(shù)必須包括常應(yīng)變項。每個單元的應(yīng)變狀態(tài) 總可以分解為不依賴于單元內(nèi)各點位置的常應(yīng)變和由各點位置決定的 變量應(yīng)變。當單元的尺寸足夠小時，單

18、元中各點的應(yīng)變趨于相等，單元 的變形比較均勻，因而常應(yīng)變就成為應(yīng)變的主要部分。 為反映單元的應(yīng) 變狀態(tài)，單元位移函數(shù)必須包括常應(yīng)變項。3）在單元內(nèi)，位移函數(shù)必須包括剛體位移項。一般情況下，單元 內(nèi)任一點的位移包括形變位移和剛體位移兩部分。形變位移與物體形狀及體積的改變相聯(lián)系，因而產(chǎn)生應(yīng)變；剛體位移只改變物體位置，不改 變物體的形狀和體積，即剛體位移是不產(chǎn)生變形的位移?？臻g一個物體 包括三個平動位移和三個轉(zhuǎn)動位移，共有六個剛體位移分量。由于一個單元牽連在另一些單元上，其他單元發(fā)生變形時必將帶動 單元做剛體位移，由此可見，為模擬一個單元的真實位移，假定的單元位移函數(shù)必須包括剛體位移項。4）位移函數(shù)

19、在相鄰單元的公共邊界上必須協(xié)調(diào)。對一般單元而言，協(xié)調(diào)性是指相鄰單元在公共節(jié)點處有相同的位移，而且沿單元邊界也有相同的位移，也就是說，要保證不發(fā)生單元的相互脫離開裂和相互侵入 重疊。要做到這一點，就要求函數(shù)在公共邊界上能由公共節(jié)點的函數(shù)值 唯一確定。對一般單元，協(xié)調(diào)性保證了相鄰單元邊界位移的連續(xù)性。但是，在板殼的相鄰單元之間，還要求位移的一階導數(shù)連續(xù)，只有 這樣，才能保證結(jié)構(gòu)的應(yīng)變能是有界量?？偟恼f來，協(xié)調(diào)性是指在相鄰單元的公共邊界上滿足連續(xù)性條件。前三條又叫完備性條件，滿足完備條件的單元叫完備單元； 第四條 是協(xié)調(diào)性要求，滿足協(xié)調(diào)性的單元叫協(xié)調(diào)單元；否則稱為非協(xié)調(diào)單元。完備性要求是收斂的必要條

20、件，四條全部滿足，構(gòu)成收斂的充分必要條 件。在實際應(yīng)用中，要使選擇的位移函數(shù)全部滿足完備性和協(xié)調(diào)性要求 是比較困難的，在某些情況下可以放松對協(xié)調(diào)性的要求。需要指出的是，有時非協(xié)調(diào)單元比與它對應(yīng)的協(xié)調(diào)單元還要好，其 原因在于近似解的性質(zhì)。假定位移函數(shù)就相當于給單元施加了約束條 件，使單元變形服從所加約束，這樣的替代結(jié)構(gòu)比真實結(jié)構(gòu)更剛一些。 但是，這種近似結(jié)構(gòu)由于允許單元分離、重疊，使單元的剛度變軟了， 或者形成了（例如板單元在單元之間的繞度連續(xù)，而轉(zhuǎn)角不連續(xù)時，剛 節(jié)點變?yōu)檩^接點）對于非協(xié)調(diào)單元，上述兩種影響有誤差相消的可能， 因此利用非協(xié)調(diào)單元有時也會得到很好的結(jié)果。在工程實踐中，非協(xié)調(diào)元必須

21、通過“小片試驗后”才能使用應(yīng)力的單元平均或節(jié)點平均處理方法最簡單的處理應(yīng)力結(jié)果的方法是取相鄰單元或圍繞節(jié)點各單元應(yīng) 力的平均值。? 1.取相鄰單元應(yīng)力的平均值這種方法最常用于3節(jié)點三角形單元中。這種最簡單而又相當實用的單元得到的應(yīng)力解在單元內(nèi)是常數(shù)?？梢詫⑵淇醋魇菃卧獌?nèi)應(yīng) 力的平均值，或是單元形心處的應(yīng)力。由于應(yīng)力近似解總是在精確 解上下振蕩，可以取相鄰單元應(yīng)力的平均值作為此兩個單元合成的 較大四邊形單元形心處的應(yīng)力。如2單元的情況下，取平均應(yīng)力可以采用算術(shù)平均，即平均應(yīng)力=（單元1的應(yīng)力+單元2的應(yīng)力）/2。也可以采用精確一些的面積加權(quán)平均，即平均應(yīng)力=單元1應(yīng)力X單元1的面積+單元2應(yīng)力X

22、單元2面積/ （單元1面積+單元2面積）當相鄰兩單元面積相差不大時，兩者的結(jié)果基本相同。在單元劃分時應(yīng)避免相鄰兩單元的面積相差太多，從而使求解的誤差相近°一般而言，3節(jié)點三角形單元的最佳應(yīng)力點是單元的中心點，此 點的應(yīng)力具有1階的精度。? 2.取圍繞節(jié)點各單元應(yīng)力的平均值首先計算圍繞該節(jié)點（i）周圍的相關(guān)單元在該節(jié)點出的應(yīng)力值 然后以他們的平均值作為該節(jié)點的最后應(yīng)力值 其中，1m是圍繞在i節(jié)點周圍的全部單元。取平均值時也可進行面積加權(quán)。有限元法求解問題的基本步驟1 .結(jié)構(gòu)離散化對整個結(jié)構(gòu)進行離散化，將其分割成若干個單元，單元間彼此通過節(jié)點相連；2 .求生各單元的剛度矩陣KK是由單元節(jié)

23、點位移量求單元節(jié)點力向量F產(chǎn)）的轉(zhuǎn)移矩陣, 其關(guān)系式為:F產(chǎn)K3 .集成總體剛度矩陣K并寫由總體平衡方程：總體剛度矩陣K是由整體節(jié)點位移向量求整體節(jié)點力向量的轉(zhuǎn)移矩陣，其關(guān)系式為F= K ,此即為總體平衡方程。4 .引入支撐條件，求由各節(jié)點的位移節(jié)點的支撐條件有兩種：一種是節(jié)點n沿某個方向的位移為零，另一種是節(jié)點n沿某個方向的位移為一給定值。5 .求生各單元內(nèi)的應(yīng)力和應(yīng)變。對于有限元方法，其基本思路和解題步驟可歸納為(1)建立積分方程，根據(jù)變分原理或方程余量與權(quán)函數(shù)正交化原理, 建立與微分方程初邊值問題等價的積分表達式，這是有限元法的出發(fā) 點°(2)區(qū)域單元剖分，根據(jù)求解區(qū)域的形狀及

24、實際問題的物理特點，將區(qū)域剖分為若干相互連接、不重疊的單元。區(qū)域單元劃分是采用有限元 方法的前期準備工作，這部分工作量比較大，除了給計算單元和節(jié)點進行編號和確定相互之間的關(guān)系之外，還要表示節(jié)點的位置坐標，同時還 需要列出自然邊界和本質(zhì)邊界的節(jié)點序號和相應(yīng)的邊界值。(3)確定單元基函數(shù)，根據(jù)單元中節(jié)點數(shù)目及對近似解精度的要求， 選擇滿足一定插值條件的插值函數(shù)作為單元基函數(shù)。有限元方法中的基函數(shù)是在單元中選取的，由于各單元 具有規(guī)則的幾何形狀，在選取基 函數(shù)時可遵循一定的法則。(4)單元分析：將各個單元中的求解函數(shù)用單元基函數(shù)的線性組合表 達式進行逼近；再將 近似函數(shù)代入積分方程，并對單元區(qū)域進行

25、積分， 可獲得含有待定系數(shù)(即單元中各節(jié)點 的參數(shù)值)的代數(shù)方程組，稱為單 元有限元方程。(5)總體合成：在得出單元有限元方程之后，將區(qū)域中所有單元有限 元方程按一定法則進行累加，形成總體有限元方程。(6)邊界條件的處理：一般邊界條件有三種形式，分為本質(zhì)邊界條件 (狄里克雷邊界條件 卜自然邊界條件(黎曼邊界條件)、混合邊界條件(柯 西邊界條件)。對于自然邊界條件，一般在積分表達式中可自動得到滿足。對于本質(zhì)邊界條件和混合邊界條件，需按一定法則對總體有限元方程進行修正滿足。(7)解有限元方程：根據(jù)邊界條件修正的總體有限元方程組，是含所有待定未知量的封閉 方程組，采用適當?shù)臄?shù)值計算方法求解，可求得

26、各節(jié)點的函數(shù)值。單元剛度矩陣的特性單元剛度矩陣無論在局部坐標系中還是在整體坐標系中都具有相同的三個特性：1)對稱性由材料力學中的位移互等定理可知，對一個構(gòu)件，作用在點j的力引起點i的繞度等于有同樣大小而作用于點i的力引起的點j的繞度，即kij=kji，表明單元剛度矩陣是一個對稱矩陣。2) 奇異性無逆陣的矩陣就叫做奇異矩陣，其行列式的值為 0,即|k|=0 ,這 一點可以從例題直接得到驗證。其物理意義是引入支撐條件之前，單元可平移。3) 分塊性有前面所講的內(nèi)容可以看出，矩陣k可以用虛線分成四塊，因此 可寫成如下的分塊形式,(e)1(e)1(e)kk 仆1112kk2122式中kmn（e）局部坐標

27、系中單元（e）按局部碼標記的節(jié)點 m、n之 間的剛度子矩陣剛架結(jié)構(gòu)中非節(jié)點載荷的處理的方法在剛架結(jié)構(gòu)以及其他較復雜的結(jié)構(gòu)上，他們所受的載荷可以直接作 用在節(jié)點上，又可以不直接作用在節(jié)點上而作用于單元節(jié)點間的其他位 置上。后一種情況下的載荷稱為非節(jié)點載荷。有限元分析時，總體剛度 方程中所用到的力向量是節(jié)點力向量。因此在進行整體分析前應(yīng)當進行載荷的移植，將作用于單元上的力移植到節(jié)點上。 移植時按靜力等 效的原則進行。處理非節(jié)點載荷一般可直接在整體坐標系內(nèi)進行，其過程為：1）將各桿單元看成一根兩端固定的梁， 分別求出兩個固定端的約束 反力。其結(jié)果可直接利用材料力學的公式求得；2）將各固定端的約束反力

28、變號，按節(jié)點進行集成，獲得各節(jié)點的等 效載荷總體剛度矩陣的集成法使用剛度矩陣獲得的方法獲得總體剛度矩陣。在此將其擴展到由整體坐標系中的單元剛度矩陣的子矩陣集成總體剛度矩陣。步驟如下：1）對一個有n個節(jié)點的結(jié)構(gòu)，將總體剛度矩陣K劃分為n xn各子 區(qū)間，然后按節(jié)點總碼的順序進行編號；2）將整體坐標系中單元剛度矩陣的各子矩陣根據(jù)其下標的兩個總碼 對號入座，寫在總體剛度矩陣相應(yīng)的子區(qū)間；3）同一子區(qū)間內(nèi)的子矩陣相加， 成為總體剛度矩陣中的相應(yīng)的子矩 陣。總體剛度矩陣的特性1）對稱性：因為由此特性，在計算機中只需存儲其上三角部分；2）奇異性：物理意義仍為在無約束的情況下，整個結(jié)構(gòu)可做剛體運 動；3）稀

29、疏性：K中有許多零子矩陣，而且在非零子矩陣中還有大量 的零元素，這種矩陣稱為稀疏矩陣。大型結(jié)構(gòu)的總體剛度矩陣一般都是 稀疏矩陣；4）分塊性：平面問題離散化時的規(guī)定1）單元之間只在節(jié)點處相連；2）所有的節(jié)點都為較接點；3）單元之間的力通過節(jié)點傳遞；4）外載荷都要移植到節(jié)點上；5）在節(jié)點位移或某一分量可以不計之處，就必須在該節(jié)點安置一個 較支座或相應(yīng)的連桿支座。通過以上的規(guī)定來建立平面有限元分析模型。結(jié)構(gòu)對稱性的利用規(guī)律一般來說，作用在對稱結(jié)構(gòu)上的載荷系統(tǒng)分為對稱的、反對稱的和一般的三種情況。1.結(jié)構(gòu)對稱，載荷對稱或反對稱這種情況下，對稱面上的邊界條件可按以下規(guī)則確定：A.在不同的對稱面上，將位移

30、分量區(qū)分為對稱分量和反對稱分量；B.將載荷也按不同的對稱面分別區(qū)分為對稱分量和反對稱分量；C.對于同一個對稱面，如載荷是對稱的，則對稱面上位移的反對稱分量為零，如載荷是反對稱的，則對稱面上位移的對稱分量為零。如果所分析的結(jié)構(gòu)對稱，但載荷是不對稱的，也不是反對稱的，這時可以將這種結(jié)構(gòu)系統(tǒng)簡化成載荷為對稱和 /或反對稱情況的組合，仍 可以簡化分析過程，提高分析的綜合效率。如圖a所示，結(jié)構(gòu)對稱，載荷一般，可將其載荷分解為圖b和圖c的組合。圖b為對稱結(jié)構(gòu)，載荷對x、y軸均為對稱，圖c為結(jié)構(gòu)對稱， 載荷對x軸反對稱、對y軸對稱，此時可取相同的四分之一進行研究， 分別施加對稱面上節(jié)點的邊界條件，進行兩次分

31、析計算，并將計算結(jié)果 迭加起來，即可得到原結(jié)構(gòu)四分之一的解答， 進而得出整個結(jié)構(gòu)的解答。利用結(jié)構(gòu)的對稱性取某一部分建立有限元模型時，往往會產(chǎn)生約束 不足現(xiàn)象。例如，若取上例中圖c的四分之一建立有限元時，根據(jù)上述分析， 在兩對稱面上應(yīng)加水平放置的滾動較支座，因此模型在垂直方向存在剛體位移。對這種約束不足問題，利用有限元分析時，必須增加附加約束, 以消除模型的剛體位移。在本例中，垂直方向可以用剛度很小的桿單元 或邊界彈簧單元連接到模型某節(jié)點上，使得既消除了模型的剛體位移， 又不致于因附加的桿單元或邊界彈簧單元剛度太大而影響結(jié)構(gòu)原有的 變形狀態(tài)。單元形態(tài)的選擇原則單元形態(tài)包括單元形狀、邊中節(jié)點的位置

32、、細長比等，在結(jié)構(gòu)離散 化過程中必須合理選擇。一般來說，為了保證有限元分析的精度，必須 是單元的形態(tài)盡可能的規(guī)則。對于三角形單元，三條邊長盡量接近，不應(yīng)出現(xiàn)大的鈍角、大的邊 長。這是因為根據(jù)誤差分析，應(yīng)力和位移的誤差都和單元的最小內(nèi)角的 正弦成反比。因而，等邊三角形單元的形態(tài)最好，它與等腰直角三角形 單元的誤差之比為sin45 :sin60 =1:1.23。但是為了適應(yīng)彈性體邊界， 以及單元由小到大逐漸過渡，不可能是所有的三角形單元都接近等邊三 角形。實際上，常常使用等腰直角三角形。對于矩形單元來說，細長比不宜過大。細長比是指單元最大尺寸和 最小尺寸之比。最優(yōu)細長比在很大程度上取決于不同方向上

33、位移梯度的 差別。梯度較大的方向，單元尺寸要小些，梯度小的方向，單元尺寸可 以大一些；如果各方向上位移梯度大致相同，則細長比越接近1,精度越高。有文獻推薦，一般情況下，為了得到較好的位移結(jié)果，細長比不 應(yīng)超過7;為了獲得較好的應(yīng)力結(jié)果，細長比不應(yīng)超過3。一般情況下， 正方形單元的形態(tài)最好。對于一般的四邊形單元應(yīng)避免過大的邊長比，過大的邊長比會導致 病態(tài)的方程組。邊界條件的確定確定邊界條件是建立有限元模型的重要一環(huán)，合理確定有限元模型 的邊界條件是成功地進行結(jié)構(gòu)有限元分析的基本要求。一般情況下，建模對象的邊界條件是明確的。根據(jù)力學模型的邊界 條件可以很容易確定其有限元模型的邊界條件。例如電線桿插入地基的 一端為固定端，橋梁一端為固定較支座，另一端為滾動較支座。但是，在機械工程中，建模對象往往是整個結(jié)構(gòu)中的一部分，在建 立有限元模型，確定其邊界條件時，必須考慮其余部分的影響。這方面 主要考慮如下兩類問題。1 .邊界位置的確定在建立連續(xù)彈性體局部區(qū)域的有限元模型時，往往取該局部區(qū)域為隔離體，取其隔離邊界條件為零位移約束，并通過試探校正確定零位移 邊界條件的位置。例如，進行齒輪齒有限元分析時，取一個輪齒的局部區(qū)域為隔離體，如圖所示，設(shè)定PQRS的邊界條件為零位移約束，