有限元簡介,應用

1、有限元方法綜述0 0引言隨著現(xiàn)代科學技術的發(fā)展，人們正在不斷建造更為快速的交通工具、更大規(guī) 模的建筑物、更大跨度的橋梁、更大功率的發(fā)電機組和更為精密的機械設備。這一切都要求工程師在設計階段就能精確地預測出產(chǎn)品和工程的技術性能，需要對結構的靜、動力強度以及溫度場、流場、電磁場和滲流等技術參數(shù)進行分析計算。 這些都可歸結為求解物理問題的控制偏微分方程式往往是不可能的。近年來在計算機技術和數(shù)值分析方法支持下發(fā)展起來的有限元分析(FEAFEA，F(xiàn)initeFinite ElementElementAnalysisAnalysis)方法則為解決這些復雜的工程分析計算問題提供了有效的途徑。有限元法的基本思

2、想是將連續(xù)的求解區(qū)域離散為一組有限個單元且按一定 方式相互聯(lián)結在一起的單元組合體。由于單元能按不同的聯(lián)結方式進行組合，且 單元本身可以有不同形狀，因此可以模型化幾何形狀復雜的求解區(qū)域。顯然，隨著單元數(shù)目的增加，即單元尺寸的縮小，或者隨著單元自由度的增加及插值函數(shù) 精度的提高，解的近似程度將不斷改進。如果單元是滿足收斂要求的，則近似解 最后將收斂于精確解。1 1有限元方法發(fā)展 1 1.1 1有限元方法起源有限元方法是受內(nèi)外動力的綜合作用而產(chǎn)生的。19431943年，柯朗在美國數(shù)學 學會公報上發(fā)表了平衡和振動問題的變分解法一文，這篇文章實際上是他 19411941年在美國數(shù)學學會演講的書面稿， 在

3、其中柯朗提出了有限元法的核心思想。 大約與柯朗同時，工程師阿格瑞斯在另一個領域獨立地提出了有限元法。柯朗和阿格瑞斯各自在數(shù)學和工程學領域獨立提出了有限元法，他們分別開創(chuàng)了有限元法的數(shù)學傳統(tǒng)和工程學傳統(tǒng)。有限元法被提出來以后，經(jīng)過一段時間的沉寂期，終于在二十世紀五十年代 和六十年代初有了很大的發(fā)展。主要表現(xiàn)為在代數(shù)表達形式、單元劃分、單元類 型選擇和解的收斂性研究上取得的突破。19601960年，克勞夫在平面應力分析中 的有限元的論文中，第一次從數(shù)學上說明了將定義域劃分成有限的單元能夠成 功的原因：他表明對一些特定類型的單元來說， 隨著單元尺寸的減小，近似解將 收斂到精確解，這就在某些情況下證明

4、了有限元法的收斂性。 并在此文中第一次 提出了“有限元法”這個名稱，這個名稱一直沿用至今，標志著有限元法早期發(fā) 展階段的結束。1 1. 2 2國內(nèi)外發(fā)展現(xiàn)狀目前，有限元方法在國外的發(fā)展主要表現(xiàn)為：第一，建立了嚴格的數(shù)學和工 程學基礎；第二，應用范圍擴展到了結構力學以外的領域；第三，收斂性得到了 進一步研究，形成了系統(tǒng)的誤差估計理論；第四，發(fā)展起了相應的、較為完備的 商業(yè)軟件包。我國的力學工作者為有限元方法的初期發(fā)展做出了許多貢獻，其中比較著名的有：陳伯屏（結構矩陣方法），錢令希（余能原理），錢偉長（廣義變分原理）， 胡海昌（廣義變分原理），馮康（有限單元法理論）。遺憾的是由于當時環(huán)境所致， 我

5、國有限元方法的研究工作受到阻礙，有限元理論的發(fā)展也逐漸與國外拉開了距 離。近年來，隨著國家對發(fā)展自主 CAECAE （ComputerComputer aidedaided EngineeringEngineering CACA日平臺 已經(jīng)愈發(fā)重視，國內(nèi)CAECAE的研究已經(jīng)逐漸走出低迷狀態(tài)，獲得了一定的發(fā)展， 而且有限元技術不再僅僅停留在高校中，而是更多的走向了企業(yè)。1 1. 3 3未來發(fā)展趨勢在大力推廣CADCAD技術的今天，從自行車到航天飛機，所有的設計制造都離 不開有限元分析計算，有限元法在工程設計和分析中將得到越來越廣泛的重視。 目前以分析、優(yōu)化和仿真為特征的 CAECAE技術在世界

6、范圍內(nèi)蓬勃發(fā)展。它通過先 進的CAECAE技術快速有效地分析產(chǎn)品的各種特性、揭示結構各類參數(shù)變化對產(chǎn)品 性能的響，進行設計方案的修改和調(diào)整，使產(chǎn)品達到性能和質(zhì)量上的最優(yōu)， 原材 料消耗最低。因此，基于計算機的分析、優(yōu)化和仿真的CAECAE技術的研究和應用， 是高質(zhì)量、高水平、低成本產(chǎn)品設計與開發(fā)的保證。當今國際上FEAFEA方法和軟件發(fā)展呈現(xiàn)出以下一些趨勢特征：1 1、從單純的結 構力學計算發(fā)展到求解許多物理場問題；2 2、由求解線性工程問題進展到分析非 線性問題；3 3、增強可視化的前置建模和后置數(shù)據(jù)處理功能；4 4、與CADCAD軟件的無縫集成；5 5、在WinWin teltel平臺上的

7、發(fā)展；6 6、擴大仿真分析的范圍；7 7、智能化、 人性化和網(wǎng)絡化。2 2有限元方法基本思想有限元方法（FEMFEM）的基礎是變分原理和加權余量法，其基本求解思想是把 計算域劃分為有限個互不重疊的單元，在每個單元內(nèi)，選擇一些合適的節(jié)點作為 求解函數(shù)的插值點，將微分方程中的變量改寫成由各變量或其導數(shù)的節(jié)點值與所 選用的插值函數(shù)組成的線性表達式，借助于變分原理或加權余量法，將微分方程離散求解。采用不同的權函數(shù)和插值函數(shù)形式，便構成不同的有限元方法。有限元方法 最早應用于結構力學，后來隨著計算機的發(fā)展慢慢用于流體力學的數(shù)值模擬。在 有限元方法中，把計算域離散剖分為有限個互不重疊且相互連接的單元，在每

8、個單元內(nèi)選擇基函數(shù)，用單元基函數(shù)的線形組合來逼近單元中的真解，整個計算域上總體的基函數(shù)可以看為由每個單元基函數(shù)組成的， 則整個計算域內(nèi)的解可以看 作是由所有單元上的近似解構成。 在河道數(shù)值模擬中，常見的有限元計算方法是 由變分法和加權余量法發(fā)展而來的里茲法和伽遼金法、最小二乘法等。根據(jù)所采 用的權函數(shù)和插值函數(shù)的不同，有限元方法也分為多種計算格式。從權函數(shù)的選 擇來說，有配置法、矩量法、最小二乘法和伽遼金法，從計算單元網(wǎng)格的形狀來 劃分，有三角形網(wǎng)格、四邊形網(wǎng)格和多邊形網(wǎng)格，從插值函數(shù)的精度來劃分，又 分為線性插值函數(shù)和高次插值函數(shù)等。不同的組合同樣構成不同的有限元計算格式。對于權函數(shù)，伽遼金

9、（GalerkinGalerkin）法是將權函數(shù)取為逼近函數(shù)中的基函數(shù)； 最小二乘法是令權函數(shù)等于余量本身，而內(nèi)積的極小值則為對代求系數(shù)的平方誤 差最?。辉谂渲梅ㄖ?，先在計算域 內(nèi)選取N N個配置點。令近似解在選定的N N個 配置點上嚴格滿足微分方程，即在配置點上令方程余量為0 0。插值函數(shù)一般由不同次幕的多項式組成，但也有采用三角函數(shù)或指數(shù)函數(shù)組成的乘積表示，但最常用的多項式插值函數(shù)。有限元插值函數(shù)分為兩大類，一類只要求插值多項式本身在插值點取已知 值，稱為拉格朗日（LagrangLagrang多項式插值；另一種不僅要求插值多項式本身， 還要求它的導數(shù)值在插值點取已知值，稱為哈密特（Herm

10、iteHermite）多項式插值。單元坐標有笛卡爾直角坐標系和無因次自然坐標，有對稱和不對稱等。常采用的無因次坐標是一種局部坐標系，它的定義取決于單元的幾何形狀，一維看作長度比， 二維看作面積比，三維看作體積比。在二維有限元中，三角形單元應用的最早， 近來四邊形等參元的應用也越來越廣。 對于二維三角形和四邊形電源單元， 常采 用的插值函數(shù)為有LagrangeLagrange插值直角坐標系中的線性插值函數(shù)及二階或更高階 插值函數(shù)、面積坐標系中的線性插值函數(shù)、二階或更高階插值函數(shù)等。3 3有限元方法求解基本步驟對于不同物理性質(zhì)和數(shù)學模型的問題，有限元求解法的基本步驟是相同的， 只是具體公式推導和運

11、算求解不同。有限元請求解問題的基本步驟通常為：1 1）建立積分方程。根據(jù)實際問題近似確定求解域的物理性質(zhì)和幾何區(qū)域。2 2）求解域離散化。將求解域近似為具有不同大小和形狀且彼此相連的有限 個單元組成的離散域，習慣上稱為有限元網(wǎng)絡劃分。顯然單元越小（網(wǎng)絡越細） 則離散域的近似程度越好、計算結果也越精確，但是計算量及誤差都將增大，因 此求解域的離散化是有限元法的核心技術之一。3 3）確定狀態(tài)變量及控制方法。一個具體的物理問題通?？梢杂靡唤M包含問 題狀態(tài)變量邊界條件的微分方程式表示， 為適合有限元求解，通常將微分方程化 為等價的泛函數(shù)形式。4 4）單元分析。對單元構造一個適合的近似解，即推導有限單元

12、的列式，其中包括選擇合理的單元坐標系，建立單元式函數(shù)，以某種方法給出單元各狀態(tài)變 量的離散關系，從而形成單元矩陣（結構力學中稱剛度陣或柔度陣）。為保證問題求解的收斂性，單元推導有許多原則要遵循。對工程應用而言，重要的是應注 意每一種單元的解題性能與約束。 例如，單元形狀應以規(guī)則為好，畸形時不僅精 度低，而且有缺值的危險，將導致無法求解。5 5）總體合成。將單元總裝形成離散域的總矩陣方程，反映對近似求解域的離散域的要求，即單元函數(shù)的連續(xù)性要滿足一定的連續(xù)條件?？傃b是在相鄰單元 結點進行，狀態(tài)變量及其導數(shù)連續(xù)性建立在結點處。6 6）求解有限元方程及結果解釋。有限元法最終導致聯(lián)立方程組。聯(lián)立方程 組

13、的求解可用直接法、選代法和隨機法。求解結果是單元結點處狀態(tài)變量的近似 值。對于計算結果的質(zhì)量，將通過與設計準則提供的允許值比較來評價并確定是 否需要重復計算。4 4總結有限元技術誕生至今已半個多世紀，隨著計算機技術與軟件技術的不斷發(fā)展 壯大，有限元的技術手段與應用范圍已經(jīng)不可同日而語，隨之而來也誕生了數(shù)以百計的有限元軟件商。如今的有限元市場是一個群雄割據(jù)的年代， 據(jù)不完全統(tǒng)計 全球有超過200200種仿真分析的軟件在被企業(yè)所使用著，昔日的巨頭命運也各不相 同，更有著無數(shù)的新貴崛起進入人們的視野。 種類繁多的有限元軟件為人類探索 未知提供了工具，而有限元軟件已經(jīng)逐漸的被更多的工程人員所接受，在有

14、限元軟件中進行著產(chǎn)品的仿真、分析與優(yōu)化，同時工程和產(chǎn)品的自主創(chuàng)新也為有限元 發(fā)展提供了強大的動力。參考文獻1林群,微分方程數(shù)值解法基礎教程，科學出版社,2003.6.2徐長發(fā)，李紅，偏微分方程數(shù)值解法，華中理工大學出版社，2005.53魯建霞，茍惠芳,有限元法的基本思想與發(fā)展過程J,機械管理開發(fā)，2009(2):7475.4趙月坤，高常，一維非齊次邊值問題的有限元法，棗莊學院學報，2011 (2):4346出師表兩漢：諸葛亮先帝創(chuàng)業(yè)未半而中道崩殂， 今天下三分，益州疲弊，此誠危急存亡之秋也。然侍衛(wèi)之臣 不懈于內(nèi)，忠志之士忘身于外者，蓋追先帝之殊遇，欲報之于陛下也。誠宜開張圣聽，以光 先帝遺德，

15、恢弘志士之氣，不宜妄自菲薄，引喻失義，以塞忠諫之路也。宮中府中，俱為一體；陟罰臧否，不宜異同。若有作奸犯科及為忠善者，宜付有司論其 刑賞，以昭陛下平明之理；不宜偏私，使內(nèi)外異法也。侍中、侍郎郭攸之、費祎、董允等，此皆良實，志慮忠純，是以先帝簡拔以遺陛下：愚 以為宮中之事，事無大小，悉以咨之，然后施行，必能裨補闕漏，有所廣益。將軍向寵，性行淑均，曉暢軍事，試用于昔日，先帝稱之曰能”是以眾議舉寵為督:愚以為營中之事，悉以咨之，必能使行陣和睦，優(yōu)劣得所。親賢臣，遠小人，此先漢所以興隆也； 親小人，遠賢臣，此后漢所以傾頹也。 先帝在時, 每與臣論此事，未嘗不嘆息痛恨于桓、 靈也。侍中、尚書、長史、參軍，此悉貞良死節(jié)之臣， 愿陛下親之、信之，則漢室之隆，可計日而待也 .FT臣本布衣，躬耕于南陽，茍全性命于亂世，不求聞達于諸侯。先帝不以臣卑鄙，猥自枉 屈，三顧臣于草廬之中，咨臣以當世之事，由是感激，遂許先帝以驅馳。后值傾覆，受任于 敗軍之際，奉命于危難之間，爾來二十