《汽車結(jié)構(gòu)有限元分析第二講有限元基礎(chǔ)理論》由會員分享,可在線閱讀,更多相關(guān)《汽車結(jié)構(gòu)有限元分析第二講有限元基礎(chǔ)理論(33頁珍藏版)》請在裝配圖網(wǎng)上搜索。
1、課件僅作為學習交流之用,不能用 于商業(yè)用途 第二講 有限元基礎(chǔ)理論 及平面問題有限元方法 講述以下問題 - 1.有限元與力學關(guān)系 2.回顧 -材料力學研究對象與研究方法 3.強度 問題 、 剛度 問題、 穩(wěn)定性 問題 4.點的應力狀態(tài) -空間問題 5.廣義 Hooke定律 6.彈性力學的基本方程 7.彈性力學問題分類 8.三大方程、三類問題、三種解法 9.平面問題 10.平面問題的有限元方法 1.有限元與力學關(guān)系 彈性力學與理論力學區(qū)別:理論力學研究對象是質(zhì)點、質(zhì) 點系與 剛體(質(zhì)點系力學與剛體力學) 。 材料力學與彈性力學研究 變形體 。 力學分支眾多: 材料力學、結(jié)構(gòu)力學、彈性力學、板殼力
2、 學、塑性力學、 斷裂力學、損傷力學、復合材料力學 、 結(jié) 構(gòu)穩(wěn)定性理論、振動理論、流體力學 、結(jié)構(gòu)動力學等 ; 有限元方法是以力學理論為基礎(chǔ),是一種現(xiàn)代數(shù)值計算方 法,是一種解決工程實際問題的數(shù)值計算工具,是現(xiàn)代設(shè) 計與分析方法的支柱! 2.回顧 -材料力學研究對象與研究方法 研究各種工程結(jié)構(gòu):常見的如下結(jié)構(gòu)元件(構(gòu)件) : ( 1) 桿、桿系、梁、柱, (長 寬和高) -材料力學 ( 2) 板 (中厚板 )、殼, (厚 長與寬) -扳殼力學 ( 3) 三維體 , -彈性力學 截面法是處理固體力學問題的最基本的方法: 通過外力(作用力和約束力)與內(nèi)力(應力)平衡求構(gòu)件的響應, 通過本構(gòu)(物理
3、)關(guān)系求變形(位移與應變), 最重要的是材料力學中的 平截面法 ,其中尤以梁的平截面假設(shè)最 為重要。 -簡化計算! 平截面假設(shè) 初始與梁的中性軸垂直的平面 ,在變形后仍垂直于 軸線 , 并且在垂直軸線方向上無變形; 梁的基本方程: 2 2 dx wd EI M 1 2 2 dx wd 2max 6 bh M )4(2 22 ayhIQ bh Q 2 3 max max max I yM 3.研究工程結(jié)構(gòu)在使用狀態(tài)下的 安 全性 、 可靠性、使用性等 ,實現(xiàn) 結(jié)構(gòu)的功能與性能。 強度 問題 (應力值不超過許用值 ) ; 剛度 問題 (變形不太大 ); 穩(wěn)定性 問題(不失穩(wěn)); 振動 問題(量值在
4、限制范圍); 碰撞問題(安全生存空間); 4 .點的應力狀態(tài) -空間問題 彈性問題 應力只取決于應變狀態(tài),與達到該狀態(tài)的過程無關(guān) 。 九個應力分量,九個應變分量(獨立變量各六個)。 單元體研究方法。 zzyzx yzyyx xzxyx zzyzx yzyyx xzxyx 2 1 2 1 2 1 2 1 2 1 2 1 6.彈性力學的基本方程 -三大方程 物理方程 x=2Gx + xy = Gxy y=2Gy + yz = Gyz z=2Gz + zx = Gzx 0 Xzyx zxyxx 0 Yzyx zyyxy 0 Zzyx zyzxz 平衡方程 x u x x v y u xy y v y
5、 y w z v yz z w z z u x w zx 幾何方程 5.各向同性彈性體 廣義 Hooke定律 EE zyx x xyxy E 12 EE xzy y yzyz E 12 EE yxz z zxzx E 12 彈性力學有 15個基本方程 : 3個平衡方程; 6個幾何方程; 6個本構(gòu)方程; 15個基本未知量 : 3個位移分量; 6個應力分量; 6個應變分量; * 加適當邊界條件。 彈性力學問題解法 -三種解法(位移法、應力法、混合法) 物理方程 應力 平衡微分方程 靜力邊界條件 變形 (位移與應變 ) 變形協(xié)調(diào)方程 (或位移單值連續(xù) ) 位移邊界條件 以位移作為未知數(shù) 幾何方程求應
6、變 物理方程求應力 位移解法 聯(lián)立求解 彈性力學問題分類 -三類邊界問題 靜力邊界問題 位移邊界問題 混合邊界問題 S u S ( X , Y , Z ) ( X ,Y ,Z ) 由位移表示的平衡微分方程 其中 是 Lplace算子 靜力邊界條件使用位移表示 位移邊界條件 0)( XxGuG 2 0)( YyGvG 2 0)( ZGwG 2 z 2 2 2 2 2 22 zyx 9. 平面問題 平面應變 物體是一柱體,軸向方向很長 所有外力(體積力和面力)都平行 于橫截面作用,且沿軸線大小不變 平面應力 沿 z方向的厚度 t均勻且很小 所有外力均作用在板的周邊和板內(nèi), 平行于板面作用,且沿厚度
7、不變 x y z y t/ 2t/ 2 z x y 平面應變特點 ( 1)位移 u=u(x,y) v=v(x,y) w = 0 ( 2)應變 平面內(nèi), x、 y、 xy 0,均為 x、 y的函數(shù); 平面外, z=xz=yz =0; ( 3)應力 z=(x+y) 平面問題的協(xié)調(diào)方程 01 )( yxzz E yx r xy xyyx 2 2 2 2 2 2 平面應力特點 ( 1)應力 在 z = 的面上各點沒有任何應力 z=zx =zy =0 在面內(nèi): x、 y、 xy 0 ( 2)應變 2 t xyxyyxyyxx E EE EE 12 yxz E xz=yz=0 ( 3) 位移 u=u(x,
8、y) v=v(x,y) w 0 平面問題平衡微分方程 0 Xyx yxx 0 Yyx yxy 平面問題幾何方程 y v y x v y u xy x u x 10. 有 限 元 方 法 概 念 平面問題的有限元法 用彈性力學經(jīng)典解法解決實際問題的主要困難在于求解偏微分方程 的復雜性 , 而有限元方法則將原來連續(xù)的彈性體離散化 , 其中最簡 單的就是采用三角形單元對彈性體進行劃分 。 把整個求解區(qū)域分成許多個有限小區(qū)域 , 這些小區(qū)域稱之為單元 。 在每個單元上構(gòu)造近似位移函數(shù) , 即進行所謂的分片插值 。 在每一個單元上求勢能 。 將所有單元上的勢能加起來得彈性體的總勢能 。 最后應用最小勢能
9、原理求解單元節(jié)點位移 。 對每個三角形單元選擇最簡單的線性函數(shù)為位移模式, 單元中任一點的位移可以通過 3個結(jié)點的位移進行插值運 算,這樣整個區(qū)域中無限多個未知位移量就可以用有限 個節(jié)點來表示,從而避免了求解覆蓋整個區(qū)域的位移函 數(shù)的困難。平面問題的有限元法,不僅可用來解決實際 問題,而且通過其相對簡單的概念,可以詳細了解用有 限元法對一般彈性體進行應力分析的基本原理和方法步 驟,了解有限元法的性能特點,使用中應注意的問題, 從而為學習后續(xù)各章節(jié)打下基礎(chǔ)。 i j m x y ( x , y ) u v 下面就以平面三角形單元闡明有限元的基本概念 單元位移模式 每個節(jié)點在單元平面內(nèi)有兩個位移分
10、量 , 相應有兩個自由度: 一個三角形單元有三個節(jié)點,共 6個節(jié)點位移分量,其單元節(jié)點位移 列陣可表示為: 位移模式可取為最簡單的線性函數(shù),包含 6個待定常數(shù) 、 。 Tiii vu ),( mji TmmjjiiTTmTjTie vuvuvu 1 6 3321 321 321 yxu yxu yxu mm jjj iii 3654 654 654 yxv yxv yxv mm jjj iii 一種簡單的線性位移函數(shù)為: 式中 、 、 為 6個待定常數(shù) , 可以由單元的節(jié)點位移確定 。 設(shè)節(jié)點 的坐標分別為 ( , )、 ( , ) 、 ( , ) , 其節(jié)點位移為 , , 將它們代入上式得:
11、 聯(lián)立求解上述公式左邊的 6個方程,可以求出待定常數(shù) : 整理后得 : yxv yxu 654 321 1 6 ix iy jx jy mx my ),(),(),( mmjjii vuvuvu 、 3321 321 321 yxu yxu yxu mm jjj iii 3654 654 654 yxv yxv yxv mm jjj iii )()()(2 1 mmmmjjjjiiii uycxbauycxbauycxbaAu )()()(2 1 mmmmjjjjiiii vycxbavycxbavycxbaAv 單元形函數(shù) 函數(shù) 表示單元內(nèi)部的位移分布形態(tài),故 可稱為單元的形態(tài)函數(shù),簡稱為形
12、函數(shù)。 得到由節(jié)點位移表達單元內(nèi)任一點位移的插值公 式,即位移模式的另一形式。 )()()(2 1 mmmmjjjjiiii uycxbauycxbauycxbaAu )()()(2 1 mmmmjjjjiiii vycxbavycxbavycxbaAv ),()(2 1 mjiycxbaAN iiii ),( mjivNvNvNv uNuNuNu mmjjii mmjjii iN iN jN mN 單元應變和應力 m m j j i i mmjjii mji mji xy y x v u v u v u bcbcbc ccc bbb A 000 000 2 1 eB mji BBBB ee
13、SBD mji SSSBDS ),( 2/)1(2/)1()1(2 2 mji bc cb cb A EBDS ii ii ii ii 單元平衡方程 整個結(jié)構(gòu)處于平衡狀態(tài),所劃分出的 一個小單元體同樣處于平衡狀態(tài),而結(jié)構(gòu) 的平衡條件可通過節(jié)點的平衡條件表示。 有限元的任務就是要建立和求解整個彈性 體的節(jié)點位移和節(jié)點力之間關(guān)系的平衡方 程。為此首先要建立每一個單元的節(jié)點位 移和節(jié)點力之間關(guān)系的平衡方程。單元平 衡方程可以利用最小勢能原理建立,也可 以利用虛功原理求解。 單元節(jié)點力列陣 : 單元節(jié)點虛位移列陣: 單元內(nèi)部引起的虛應變 : 根據(jù)虛功原理:外力虛功等于內(nèi)力虛功 。 所以節(jié)點力 在節(jié)點的
14、虛位移上所作的虛功應等于單元內(nèi)部應力在虛應 變上所作的虛功 。 這就是單元保持平衡狀態(tài)所必須滿足的 條件 , 即單元的平衡條件 。 Tmmjjiie YXYXYXF Tmmjjii vuvuvu * Tzyx * t d xd yF TeeT * eTe t d x d yBDBF eee kF t dx dyBDBk Te e mmmjmi jmjjji imijii mji T m T j T i e kkk kkk kkk AtBBBD B B B k 單元剛度矩陣 利用虛功方程來建立剛度方程,其實質(zhì)就是 單元的平衡方程 。 單元剛度矩陣具有以下性質(zhì): (1) 單元剛度矩陣中每個元素有明
15、確的物理意義。其物理意 義是單位節(jié)點位移分量所引起的節(jié)點力。例如, 是表示 當單元第 n個自由度產(chǎn)生單位位移而其它自由度固定時, 在第 m個自由度產(chǎn)生的節(jié)點力。 (2) 是對稱矩陣。其元素之間有如下關(guān)系: ,這個特 性是由彈性力學中功的互等定理所決定的。 ( 3) 是奇異矩陣。其每一行每一列元素之和均為零,物 理意義就是:在無約束的條件下,單元可作剛體運動。 根據(jù)行列式性質(zhì),可知值也為零。 mnk srrs kk ek ek 單元等效節(jié)點載荷 外載荷必須作用在節(jié)點上,而實際的 外載荷又往住并不是通過節(jié)點作用的。 因此,必須將這些非節(jié)點載荷按一定原 則移置到節(jié)點上,即所謂等效節(jié)點載荷 處理。這種
16、移置必須滿足靜力等效原則 。 處理單元內(nèi)的集中力、體力和單元邊界上 的分布力 ,慣性力則作用在整個結(jié)構(gòu)上。 總剛度矩陣 當以有限個單元通過有限個節(jié)點連接而成的組 合體來代替實際的連續(xù)體結(jié)構(gòu)而受力變形時,顯然 它們必須滿足整個結(jié)構(gòu)的變形連續(xù)條件和平衡條件。 在整體分析中,利用節(jié)點為分析對象,根據(jù)各 節(jié)點的靜力平衡條件,即可建立起組合體所有節(jié)點 的靜力平衡方程式。把它們匯集在一起,得到的平 衡方程組就代表了整個結(jié)構(gòu)的平衡條件。進行整體 分析,即是將各個單元的平衡方程集合在一起,得 到結(jié)構(gòu)的整體平衡方程。 K為結(jié)構(gòu)的整體剛度矩陣,一般稱為總剛度矩陣, 其維數(shù)為 2n 2n??蓪懗煞謮K形式。 RK T
17、TnTTT 321 TTnTTT RRRRR 321 解題步驟與算例 (1)首先繪出結(jié)構(gòu)幾何簡圖,在此基礎(chǔ)上將結(jié)構(gòu)離散化。平面問題 采用三角形單元 (其他形狀單元以后講述 ),所以其離散就是將 計算對象劃分成許多三角形單元。包括:進行節(jié)點編號、單元 編號,任選一直角坐標系,定出所有節(jié)點的坐標值等等。確定 載荷和邊界約束條件,將各單元所受的非節(jié)點載荷,包括體力、 面力以及可能有的集中力按虛功等效原則移置到節(jié)點上,并將 各節(jié)點上的這些載荷(包括直接作用在節(jié)點上的集中載荷)分 別按相同方向 疊加等。 (2)其次進行單元分析、組集總剛度矩陣、求單元應力和節(jié)點應力。 前處理 計算 后處理 nnnnn n
18、 n R R R nKKK KKK KKK 2 1 2 1 21 22221 11211 平面問題的離散化 單元類型的選擇 單元的大小 單元有密有疏 不同厚度或不同材料處,應取作為單元的邊界線 平面問題的有限元法,不僅有實際意義, 而且通過其相對簡單的概念,可以詳細了解用 有限元法對一般彈性體進行應力分析的基本原 理和方法步驟,了解有限元法的性能特點,使 用中應注意的問題,從而為學習以后各章打下 基礎(chǔ)。 有限元解法的三個主要步驟就是: 離散化、單元分析、整體分析 。 平面高階單元 四節(jié)點矩形單元: 為了提高有限單元法計算結(jié)果的精 度,除了增加單元數(shù)目外,還常采用具有較高次位移函數(shù)的單元。 等參數(shù)單元: 三角形單元和矩形單元的位移模式和坐換 變換式都采用了相同的形函數(shù)。例平面四節(jié)點任意四邊形等參單 元。 xyyxyxv xyyxyxu 8765 4321 )( )( , , i i i uNu 4 1 )()( , i i i vNv 4 1 )()( , 合肥工業(yè)大學 車輛工程系 第二講結(jié)束語 溫故而知新