有限元法介紹

關于有限元法介紹第1頁，講稿共40頁，2023年5月2日，星期三一、有限元法的概念二、基本計算步驟三、發(fā)展與應用第2頁，講稿共40頁，2023年5月2日，星期三基本思想第3頁，講稿共40頁，2023年5月2日，星期三一、幾個基本概念有限元法是把要分析的連續(xù)體假想地分割成有限個單元所組成的組合體，簡稱離散化。這些單元僅在頂角處相互聯(lián)接，稱這些聯(lián)接點為節(jié)點。離散化的組合體與真實彈性體的區(qū)別在于：組合體中單元與單元之間的聯(lián)接除了結點之外再無任何關聯(lián)。但是這種聯(lián)接要滿足變形協(xié)調條件，即不能出現(xiàn)裂縫，也不允許發(fā)生重疊。顯然，單元之間只能通過結點來傳遞內力。通過結點來傳遞的內力稱為節(jié)點力，作用在結點上的荷載稱為節(jié)點荷載。當連續(xù)體受到外力作用發(fā)生變形時，組成它的各個單元也將發(fā)生變形，因而各個結點要產生不同程度的位移，這種位移稱為節(jié)點位移。第4頁，講稿共40頁，2023年5月2日，星期三單元節(jié)點節(jié)點力節(jié)點荷載節(jié)點位移第5頁，講稿共40頁，2023年5月2日，星期三在有限元中，常以結點位移作為基本未知量。并對每個單元根據(jù)分塊近似的思想，假設一個簡單的函數(shù)近似地表示單元內位移的分布規(guī)律，再利用力學理論中的變分原理或其他方法，建立結點力與位移之間的力學特性關系，得到一組以結點位移為未知量的代數(shù)方程，從而求解結點的位移分量。然后利用插值函數(shù)確定單元集合體上的場函數(shù)。顯然，如果單元滿足問題的收斂性要求，那么隨著縮小單元的尺寸，增加求解區(qū)域內單元的數(shù)目，解的近似程度將不斷改進，近似解最終將收斂于精確解。第6頁，講稿共40頁，2023年5月2日，星期三有限元法的基本計算步驟物體離散化單元特性分析分析單元的力學性質選擇位移模式計算等效節(jié)點力單元組集求解未知節(jié)點位移第7頁，講稿共40頁，2023年5月2日，星期三物體離散化單元選擇：應根據(jù)連續(xù)體的形狀選擇最能完滿地描述連續(xù)體形狀的單元。常見的單元有：桿單元，梁單元，三角形單元，矩形單元，四邊形單元，曲邊四邊形單元，四面體單元，六面體單元以及曲面六面體單元等等。單元劃分：進行單元劃分，單元劃分完畢后，要將全部單元和結點按一定順序編號，每個單元所受的荷載均按靜力等效原理移植到結點上，并在位移受約束的結點上根據(jù)實際情況設置約束條件。第8頁，講稿共40頁，2023年5月2日，星期三3維實體的4面體單元劃分平面的三角形單元劃分3維實體的6面體單元劃分第9頁，講稿共40頁，2023年5月2日，星期三將某個工程結構離散為由各種單元組成的計算模型，這一步又稱作單元剖分或網(wǎng)格劃分。離散后單元于單元之間利用單元的節(jié)點相互連接起來；單元節(jié)點的設置、性質、數(shù)目等應視問題的性質，描述變形形態(tài)的需要和計算進度而定。用有限元分析計算所獲得的結果只是近似的。如果劃分單元數(shù)目非常多而又合理，則所獲得的結果就與實際情況相符合。網(wǎng)格劃分的好壞將直接影響到計算結果的準確性和計算進度，甚至會因為網(wǎng)格劃分不合理而導致計算不收斂。網(wǎng)格的劃分主要取決于專業(yè)知識和經驗積累。一個水平高的FEA工程師，80%的時間是用在網(wǎng)格劃分上。對于一般的問題，各種FEA均能自動的進行合理的網(wǎng)格劃分。hypermesh最目前好的劃分網(wǎng)格工具。第10頁，講稿共40頁，2023年5月2日，星期三單元特性分析分析單元的力學性質根據(jù)單元的材料性質、形狀、尺寸、節(jié)點數(shù)目、位置及其含義等，找出單元節(jié)點力和節(jié)點位移的關系式，這是單元分析中的關鍵一步。此時需要應用彈性力學中的幾何方程和物理方程來建立力和位移的方程式，從而導出單元剛度矩陣，這是有限元法的基本步驟之一。第11頁，講稿共40頁，2023年5月2日，星期三平面問題的三角形單元的例子：單元有三個結點I、J、M，每個結點有兩個位移u、v和兩個結點力U、V。節(jié)點位移節(jié)點力取決于材料性質、形狀、尺寸第12頁，講稿共40頁，2023年5月2日，星期三選擇位移模式：在反映力和位移的關系式中，依據(jù)那一個量是未知量，可建立不同的模型。位移法：選擇節(jié)點位移作為基本未知量稱為位移法；力法：選擇節(jié)點力作為基本未知量時稱為力法；混合法：取一部分節(jié)點力和一部分節(jié)點位移作為基本未知量時稱為混合法。位移法易于實現(xiàn)計算自動化，所以，在有限單元法中位移法應用范圍最廣。第13頁，講稿共40頁，2023年5月2日，星期三計算等效節(jié)點力：將外在的負載力等效到各個節(jié)點上。物體離散化后，假定力是通過節(jié)點從一個單元傳遞到另一個單元。但是，對于實際的連續(xù)體，力是從單元的公共邊傳遞到另一個單元中去的。因而，這種作用在單元邊界上的表面力、體積力和集中力都需要等效的移到節(jié)點上去，也就是用等效的節(jié)點力來代替所有作用在單元上得力。彈性體有限元模型第14頁，講稿共40頁，2023年5月2日，星期三有限元法的基本計算步驟單元組集利用結構力的平衡條件和邊界條件把各個單元按原來的結構重新連接起來，形成整體的有限元方程。對由各個單元組成的整體進行分析，建立節(jié)點外載荷與結點位移的關系，以解出節(jié)點位移，這個過程為整體分析。i節(jié)點的節(jié)點力：

i節(jié)點的平衡方程：集中力單元節(jié)點力第15頁，講稿共40頁，2023年5月2日，星期三單元組集

最終，將所有單元組合起來得到整體的方程：[K]{δ}={R}

[K]——整體剛度矩陣；{δ}——全部結點位移組成的列陣；{R}——全部結點荷載組成的列陣。

在位移法中，只有{δ}是未知的，求解該線性方程組就可得到各結點的位移。將結點位移代入相應方程中可求出單元的應力分量。有限元法不僅可以求結構體的位移和應力，還可以對結構體進行穩(wěn)定性分析和動力分析。例如，結構體的整體動力方程：

[M]{δ}+[C]{δ}+[K]{δ}={F}

[M]——整體質量矩陣；[C]——整體阻尼矩陣；[K]——整體剛度矩陣；{δ}——整體結點位移向量；{F}——整體結點荷載向量。

求出結構的自激振動頻率、振型等動力響應，以及動變形和動應力等。第16頁，講稿共40頁，2023年5月2日，星期三求解未知節(jié)點位移

可以根據(jù)方程組的具體特點來選擇合適的計算方法。傳統(tǒng)有限元分析的數(shù)值計算方法之中，有直接計算法（DirectSolver）與迭代法（Iterative）兩種。由于在過去的經驗中，迭代法一直無法直接而有效的保證數(shù)值計算的收斂性，因此，直接計算法在多數(shù)有限元素分析軟件中，仍然是一種主流的計算方法。第17頁，講稿共40頁，2023年5月2日，星期三線性靜力分析的基本矩陣方程單元剛度矩陣

[K]=剛度矩陣{F}=力向量(已知){u}=由{F}引起的未知位移向量

第18頁，講稿共40頁，2023年5月2日，星期三總體剛度矩陣第19頁，講稿共40頁，2023年5月2日，星期三F=10000磅L=10英寸YINL=10英寸面積A=1.5英寸2彈性模量E=10×106磅/英寸2面積A=2.0英寸2彈性模量E=10×106磅/英寸2例子：1、建立結構有限元模型

L=10英寸A=1.5英寸2L=10英寸A=2.0英寸2結點3#結點1#結點2#單元②單元①第20頁，講稿共40頁，2023年5月2日，星期三2、形成單元剛度矩陣3、總裝剛度矩陣第21頁，講稿共40頁，2023年5月2日，星期三4、施加邊界條件5、施加作用載荷第22頁，講稿共40頁，2023年5月2日，星期三6、求解矩陣方程7、計算單元應力第23頁，講稿共40頁，2023年5月2日，星期三三、有限元法的發(fā)展與與應用在大力推廣CAD技術的今天，從自行車到航天飛機，所有的設計制造都離不開有限元分析計算，F(xiàn)EA在工程設計和分析中將得到越來越廣泛的重視。國際上早在20世紀50年代末、60年代初就投入大量的人力和物力開發(fā)具有強大功能的有限元分析程序。其中最為著名的是由美國國家宇航局（NASA）在1965年委托美國計算科學公司和貝爾航空系統(tǒng)公司開發(fā)的NASTRAN有限元分析系統(tǒng)。該系統(tǒng)發(fā)展至今已有幾十個版本，是目前世界上規(guī)模最大、功能最強的有限元分析系統(tǒng)。第24頁，講稿共40頁，2023年5月2日，星期三目前應用較多的通用有限元軟件如下表所列：另外還有許多針對某類問題的專用有限元軟件，例如金屬成形分析軟件Deform、Autoform。軟件名稱簡介MSC/Nastran著名結構分析程序，最初由NASA研制MSC/Dytran動力學分析程序MSC/Marc非線性分析軟件ANSYS通用結構分析軟件ADINA非線性分析軟件ABAQUS非線性分析軟件第25頁，講稿共40頁，2023年5月2日，星期三第26頁，講稿共40頁，2023年5月2日，星期三轉向機構支架的強度分析（劉道勇，東風汽車工程研究院動，用MSC/Nastran完成）第27頁，講稿共40頁，2023年5月2日，星期三金屬成形過程分析（用Deform軟件完成）:分析金屬成形過程中的各種缺陷

型材擠壓成形的分析。型材在擠壓成形的初期，容易產生形狀扭曲螺旋齒輪成形過程的分析T形鍛件的成形分析第28頁，講稿共40頁，2023年5月2日，星期三焊接殘余應力分析（用Sysweld完成）結構與焊縫布置焊接過程的溫度分布與軸向殘余應力第29頁，講稿共40頁，2023年5月2日，星期三有限元方法的應用BMW曲軸的感應淬火（用SysWeld軟件完成）第30頁，講稿共40頁，2023年5月2日，星期三有限元方法的應用復雜形狀工件的組織轉變預測（石偉，用NSHT3D完成）：預測工件的組織分布和機械性能二分之一工件的有限元模型淬火3.06min時的溫度分布淬火3.06min時的馬氏體分布第31頁，講稿共40頁，2023年5月2日，星期三其他的應用還包括：電磁學、流體力學、電磁場等等等等等等等等等…第32頁，講稿共40頁，2023年5月2日，星期三1、從單純結構力學計算發(fā)展到求解許多物理場問題有限元分析方法最早是從結構化矩陣分析發(fā)展而來，逐步推廣到板、殼和實體等連續(xù)體固體力學分析，實踐證明這是一種非常有效的數(shù)值分析方法。有限元方法已發(fā)展到流體力學、溫度場、電傳導、磁場、滲流和聲場等問題的求解計算，最近又發(fā)展到求解幾個交叉學科的問題。

例如當氣流流過一個很高的鐵塔產生變形，而塔的變形又反過來影響到氣流的流動……這就需要用固體力學和流體動力學的有限元分析結果交叉迭代求解，即所謂"流固耦合"的問題。發(fā)展趨勢：第33頁，講稿共40頁，2023年5月2日，星期三2、由求解線性工程問題進展到分析非線性問題線性理論已經遠遠不能滿足設計的要求。例如：航天和動力工程的高溫部件存在熱變形和熱應力，要考慮材料的非線性問題；諸如塑料、橡膠和復合材料等各種新材料的出現(xiàn)，只有采用非線性有限元算法才能解決。非線性的數(shù)值計算是很復雜的，很難為一般工程技術人員所掌握。為此近年來國外一些公司花費了大量的人力和投資開發(fā)諸如MARC、ABQUS和ADINA等專長于求解非線性問題的有限元分析軟件，并廣泛應用于工程實踐。第34頁，講稿共40頁，2023年5月2日，星期三3、增強可視化的前置建模和后置數(shù)據(jù)處理功能隨著數(shù)值分析方法的逐步完善，尤其是計算機運算速度的飛速發(fā)展，整個計算系統(tǒng)用于求解運算的時間越來越少，而數(shù)據(jù)準備和運算結果的表現(xiàn)問題卻日益突出。在現(xiàn)在的工程工作站上，求解一個包含10萬個方程的有限元模型只需要用幾十分鐘。工程師在分析計算一個工程問題時有80%以上的精力都花在數(shù)據(jù)準備和結果分析上。第35頁，講稿共40頁，2023年5月2日，星期三與CAD軟件的無縫集成當今有限元分析系統(tǒng)的另一個特點是與通用CAD軟件的集成使用，即：在用CAD軟件完成部件和零件的造型設計后，自動生成有限元網(wǎng)格并進行計算，如果分析的結果不符合設計要求則重新進行造型和計算，直到滿意為止，從而極大地提高了設計水平和效率。當今所有的商業(yè)化有限元系統(tǒng)商都開發(fā)了和著名的CAD軟件（例如Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、Bentley和AutoCAD等）的接口。第36頁，講稿共40頁，2023年5月2日，星期三有限元軟件使用的基本步驟三個步驟：前處理（PREPROCESSION），求解（SOLUTION），后處理POSTPROCESSION）前處理單元選擇材料定義幾何建模網(wǎng)格劃分模型局部調整施加荷載第37頁，講稿