第2講-CAE技術(shù)基本求解過程

11.ANSYS2.ADINA3.MARC4.ABAQUS5.ADMAS6.PRO/E-CAD7.PhoenicsCFD8.Fluent6.3CFD我校購置的正版工程軟件2第2章CAE技術(shù)基本求解過程一、有限單元法基本知識二、線性分析有限元法的基本計算步驟三、非線性分類及有限元法基本流程第2講3一、有限單元法的基本概念

有限元分析（FEA）

利用數(shù)學近似的方法對真實物理系統(tǒng)（幾何和載荷工況）進行模擬。還利用簡單而又相互作用的元素，即單元，就可以用有限數(shù)量的未知量去逼近無限未知量的真實系統(tǒng)。4F一、有限單元法的基本概念輪齒有限元模型網(wǎng)格節(jié)點邊界單元有限單元節(jié)點數(shù)是有限的單元數(shù)目也是有限的5按有限元規(guī)則，下面單元編號正確與否?kij3節(jié)點單元6節(jié)點單元kijlmn單元編號那些是正確的？1.k,j,i2.i,j,k3.i,k,j4.j,k,i單元編號那些是正確的？1.k,j,i,n,m,l2.i,j,k,l,m,n3.i,n,j,m,l,k4.j,k,l,m,n,i5.i,k,m,j,l,n6下面單元編號正確與否?53121021675438單元編號那些是正確的？1.5,3,12,102.3,12,10,53.10,12,3,54.10,5,3,12單元編號那些是正確的？1.1,2,3,4,6,7,8,52.1,2,3,4,8,5,6,73.3,4,1,2,7,8,5,64.2,3,4,1,6,7,8,54節(jié)點單元8節(jié)點單元7注意：單元編號逆時針方向!!i53kj2112106754388網(wǎng)格劃分912345678910!*K,1,,,,K,2,25,,,

K,3,25,20,,

K,4,14,22.5,,

K,5,11,25,,

K,6,9,30,,

K,7,8,39,,

K,8,6,43,,

K,9,3.5,47,,K,10,0,47,,

!關(guān)鍵點1．物體離散化10位移等值線圖11剪應(yīng)力等值線圖12等效應(yīng)力圖132．單元特性分析1)選擇位移模式選擇節(jié)點位移作為基本未知量時稱為位移法；選擇節(jié)點力作為基本未知量時稱為力法；取一部分節(jié)點力和一部分節(jié)點位移作為基本未知量時稱為混合法，結(jié)構(gòu)物離散化后位移、應(yīng)變和應(yīng)力等由節(jié)點位移來表示2)分析單元的力學性質(zhì)位移法力法混合法幾何方程物理方程力和位移的方程式單元剛度矩陣易于實現(xiàn)計算自動化143)計算等效節(jié)點力3．單元組集結(jié)構(gòu)力的平衡條件邊界條件Kq=f表面力體積力集中力等效地移到節(jié)點上各個單元按原來的結(jié)構(gòu)重新聯(lián)接起來整體的有限元方程整體結(jié)構(gòu)的剛度矩陣節(jié)點位移列陣載荷列陣4求解未知節(jié)點位移15二、線性分析有限元法的基本計算步驟1單元特性的導(dǎo)出方法(1)直接方法(2)虛功原理法(3)能量變分原理方法(4)迦遼金法。16(1)直接方法17Fe=Keqe1819(2)虛功原理法1）設(shè)定位移函數(shù)u(x,y)=a1+a2x+a3yv(x,y)=a4+a5x+a5yd=Nqe位移分量形狀函數(shù)單元內(nèi)位移平面問題202)由位移函數(shù)求應(yīng)變彈性力學理論應(yīng)變矩陣213)由虎克定律，通過應(yīng)變求應(yīng)力彈性矩陣224)由虛功原理求單元剛度矩陣外力做的虛功內(nèi)力做的虛功單元剛度矩陣233.能量變分原理方法l）最小勢能原理彈性體受外力作用產(chǎn)生變形時用變分法求能量泛函的極值方法就是能量變分原理系統(tǒng)總勢能變形勢能勢能降低虛功原理泛函函數(shù)的函數(shù)242）能量變分原理的應(yīng)用彎曲產(chǎn)生的應(yīng)變能撓度分布載荷直梁的基本微分方程外載荷作功252有限單元法的解題步驟1)單元剖分和插值函數(shù)的確定2)單元特性分析應(yīng)變：=Bqe應(yīng)力：=D=DBqe虛功原理節(jié)點力與位移的關(guān)系fe=Keqe節(jié)點位移：d=Nqe263)單元組集4)解有限元方程5)計算應(yīng)力結(jié)構(gòu)平衡方程組進行邊界條件處理整體結(jié)構(gòu)平衡方程組F=Kq總的載荷列陣整體結(jié)構(gòu)的剛度矩陣整體結(jié)構(gòu)所有節(jié)點的位移列陣解出節(jié)點位移273.用三角形單元進行靜力分析的實例12cm6cm100N/cm

厚度t=0.1cm,u=0,E=2.1e7N/cm2281單元劃分12對稱性取一半300N300N12341xyijky2xijk293組成整體剛度矩陣4邊界條件處理5線性方程組的建立與求解6單元應(yīng)力分量的計算2計算單元剛度矩陣K1K2K30ANSYS軟件驗證：平面應(yīng)力問題1324300N300N單元1單元2單元1PLANE2SXSYSXY-285.71-1285.7-285.71單元2

SXSYSXY.00000-714.29-714.2931三、非線性分類及有限元法基本流程狀態(tài)變化；幾何非線性；材料非線性。非線性結(jié)構(gòu)非線性結(jié)構(gòu)行為的普通例子32三、非線性分類及有限元法基本流程狀態(tài)變化；幾何非線性；材料非線性。非線性結(jié)構(gòu)幾何非線性33固體力學問題，從本質(zhì)上講是非線性的，線性假設(shè)僅是實際問題中的一種簡化

假設(shè)節(jié)點位移無限小；材料的應(yīng)力與應(yīng)變關(guān)系滿足虎克定律：加載時邊界條件的性質(zhì)保持不變。線性彈性體系非線性問題否則三、非線性分類及有限元法基本流程34非線性問題材料非線性幾何非線性由于材料的應(yīng)力與應(yīng)變關(guān)系的非線性引起的結(jié)構(gòu)的位移使體系的受力狀態(tài)發(fā)生了顯著變化，以至不能采用線性體系的分析方法①大位移小應(yīng)變問題;②大位移大應(yīng)變問題;③結(jié)構(gòu)的變位引起外載荷大小、方向或邊界支承條件的變化等問題。35有限單元法分析非線性問題三個基本步驟組1．單元分析2．整體組集3．非線性方程組的求解材料非線性問題，使用材料的非線性本構(gòu)關(guān)系幾何非線性問題，計算位移高階導(dǎo)數(shù)的應(yīng)變－位移矩陣整體剛度方程寫成增量形式增量法迭代法混合法36增量逐步解法的基本思想假定t時刻的解已知，t為選擇的時間增量，在t+t時刻，有：t時刻的解為已知外載荷節(jié)點力向量單元應(yīng)力引起的節(jié)點力向量37節(jié)點力增量向量材料和幾何條件的切向剛度矩陣節(jié)點位移增量向量38修正Newton法的迭代公式迭代步數(shù)初始值t時刻的解迭代所用的初始值正是t時刻的解滿足精度指標時迭代終止39Area=b*h=12.5mm2Izz=b*h3/12=26.0417