有限元法緒論

1、1第一章 有限元法緒論(xln)共三十八頁(yè)2第一節(jié)有限元法的基本思路第二節(jié)有限元法的計(jì)算(j sun)步驟第四節(jié) 有限元法的用途(yngt)第三節(jié) 有限元法的成就第一章 有限元法緒論共三十八頁(yè)3 有限元法的基本思路可以歸結(jié)為：將連續(xù)系統(tǒng)分割成有限個(gè)分區(qū)或單元(dnyun)，對(duì)每個(gè)單元(dnyun)提出一個(gè)近似解，再將所有單元(dnyun)按標(biāo)準(zhǔn)方法組合成一個(gè)與原有系統(tǒng)近似的系統(tǒng)。 等截面直桿在自重作用下的材料力學(xué)解答有限元法緒論(xln)第一節(jié)有限元法的基本思路 受自重作用的等截面直桿 共三十八頁(yè) 受自重作用的等截面直桿如圖所示，桿的長(zhǎng)度為L(zhǎng)，截面積為A，彈性模量為E，單位長(zhǎng)度的重量為q，桿的

2、內(nèi)力為N。試求：桿的位移(wiy)分布，桿的應(yīng)變和應(yīng)力。共三十八頁(yè)等截面(jimin)直桿在自重作用下的有限元法解答離散(lsn)后的直桿 1、離散化 如右圖所示，將直桿劃分成n個(gè)有限段，有限段之間通過(guò)一個(gè)鉸接點(diǎn)連接。稱(chēng)兩段之間的連接點(diǎn)為結(jié)點(diǎn)，稱(chēng)每個(gè)有限段為單元。 第i個(gè)單元的長(zhǎng)度為L(zhǎng)i，包含第i，i+1個(gè)結(jié)點(diǎn)。共三十八頁(yè)2、用單元(dnyun)節(jié)點(diǎn)位移表示單元(dnyun)內(nèi)部位移 第i個(gè)單元中的位移用所包含的結(jié)點(diǎn)位移來(lái)表示，其中(qzhng)為第i結(jié)點(diǎn)的位移，第i個(gè)單元的應(yīng)變?yōu)閼?yīng)力為，內(nèi)力為為第i結(jié)點(diǎn)的坐標(biāo)。(*)共三十八頁(yè)3、把外載荷(zi h)集中到節(jié)點(diǎn)上 把第i單元和第i+1單元重量

3、(zhngling)的一半集中到第i+1結(jié)點(diǎn)上。 集中單元重量 4、建立結(jié)點(diǎn)的力平衡方程 對(duì)于第i+1結(jié)點(diǎn)，由力的平衡方程可得：令，并將(*)代入得：令 共三十八頁(yè)根據(jù)(gnj)約束條件，對(duì)于(duy)第n+1個(gè)結(jié)點(diǎn)， 建立所有結(jié)點(diǎn)的力平衡方程，可以得到由n+1個(gè)方程構(gòu)成的方程組，可解出n+1個(gè)未知的結(jié)點(diǎn)位移。共三十八頁(yè)9有限元法緒論(xln)第二節(jié)有限元法的計(jì)算(j sun)步驟 有限元法的計(jì)算步驟歸納為以下三個(gè)基本步驟：網(wǎng)格劃分，單元分析，整體分析。1.網(wǎng)格劃分 有限元法的基礎(chǔ)是用有限個(gè)單元體的集合來(lái)代替原有的連續(xù)體。因此首先要對(duì)彈性體進(jìn)行必要的簡(jiǎn)化，再將彈性體劃分為有限個(gè)單元組成的離散體

4、。單元之間通過(guò)單元節(jié)點(diǎn)相連接。由單元、結(jié)點(diǎn)、結(jié)點(diǎn)連線(xiàn)構(gòu)成的集合稱(chēng)為網(wǎng)格。共三十八頁(yè)10有限元法緒論(xln) 通常把三維實(shí)體劃分成4面體或6面體單元的網(wǎng)格，平面(pngmin)問(wèn)題劃分成三角形或四邊形單元的網(wǎng)格。四面體四節(jié)點(diǎn)單元 三維實(shí)體的四面體單元?jiǎng)澐?共三十八頁(yè)11有限元法緒論(xln)六面體8節(jié)點(diǎn)(ji din)單元 三維實(shí)體的六面體單元?jiǎng)澐?共三十八頁(yè)三角形3節(jié)點(diǎn)(ji din)單元 平面問(wèn)題的三角形單元(dnyun)劃分 共三十八頁(yè)四邊形4節(jié)點(diǎn)(ji din)單元 平面問(wèn)題的四邊形單元(dnyun)劃分 共三十八頁(yè)2.單元分析 對(duì)于彈性力學(xué)問(wèn)題，單元分析，就是建立各個(gè)單元的節(jié)點(diǎn)位移和

5、節(jié)點(diǎn)力之間的關(guān)系式。 由于將單元的節(jié)點(diǎn)位移作為基本變量，進(jìn)行單元分析首先要為單元內(nèi)部(nib)的位移確定一個(gè)近似表達(dá)式，然后計(jì)算單元的應(yīng)變、應(yīng)力，再建立單元中節(jié)點(diǎn)力與節(jié)點(diǎn)位移的關(guān)系式。 以平面問(wèn)題的三角形3結(jié)點(diǎn)單元為例。如下(rxi)圖所示，單元有三個(gè)結(jié)點(diǎn)I、J、M，每個(gè)結(jié)點(diǎn)有兩個(gè)位移u、v和兩個(gè)結(jié)點(diǎn)力U、V。共三十八頁(yè)三角形3結(jié)點(diǎn)(ji din)單元 單元的所有(suyu)結(jié)點(diǎn)位移、結(jié)點(diǎn)力，可以表示為結(jié)點(diǎn)位移向量 結(jié)點(diǎn)位移 結(jié)點(diǎn)力 單元的結(jié)點(diǎn)位移和結(jié)點(diǎn)力之間的關(guān)系：共三十八頁(yè)3.整體分析 對(duì)由各個(gè)單元組成的整體進(jìn)行分析，建立節(jié)點(diǎn)外載荷與結(jié)點(diǎn)位移的關(guān)系，以解出結(jié)點(diǎn)位移，這個(gè)過(guò)程為整體分析。再以

6、彈性力學(xué)(l xu)的平面問(wèn)題為例，如下圖所示，在邊界結(jié)點(diǎn)i上受到集中力 作用(zuyng)。 結(jié)點(diǎn)i是三個(gè)單元的結(jié)合點(diǎn)，因此要把這三個(gè)單元在同一結(jié)點(diǎn)上的結(jié)點(diǎn)力匯集在一起建立平衡方程。整體分析 共三十八頁(yè)i結(jié)點(diǎn)(ji din)的結(jié)點(diǎn)(ji din)力：i結(jié)點(diǎn)(ji din)的平衡方程：共三十八頁(yè)18 隨著有限元法的不斷應(yīng)用和計(jì)算機(jī)的高速發(fā)展,很多工程軟件已成為(chngwi)有效的科學(xué)計(jì)算手段并大量用于世紀(jì)工程中. 如：ANSYS ADINA ABAQUS SAP5 SUPERSAP NASTRAN 這些軟件雖僅是計(jì)算工具(gngj)，但他們的發(fā)展也為有限元法的應(yīng)用拓寬了領(lǐng)域，趨近于多元化，多

7、角度多功能的計(jì)算。特別是前后處理的不斷更新，優(yōu)化，更使軟件運(yùn)用的得心應(yīng)手。有限元法緒論第三節(jié)有限元法的成就共三十八頁(yè)19 Ansys軟件是融結(jié)構(gòu)、熱、流體、電磁、聲學(xué)于一體的大型(dxng)通用有限元分析軟件，可廣泛用于核工業(yè)、鐵道、石油化工、航空航天、機(jī)械制造、能源、汽車(chē)交通、國(guó)防軍工、電子、土木工程、造船、生物醫(yī)學(xué)、輕工、地礦、水利、日用家電等一般工業(yè)及科學(xué)研究。Ansys 簡(jiǎn)介(jin ji)有限元法緒論共三十八頁(yè)20 該軟件可在大多數(shù)計(jì)算機(jī)及操作系統(tǒng)中運(yùn)行，從PC機(jī)到工作站直至巨型計(jì)算機(jī)，Ansys文件在其所有的產(chǎn)品系列和工作平臺(tái)上均兼容。Ansys軟件多物理(wl)場(chǎng)耦合的功能,允許

8、在統(tǒng)一模型上進(jìn)行各式各樣的耦合計(jì)算,如：熱流體耦合,磁電耦合,以及電磁流體熱耦合，確保了所有的Ansys用戶(hù)的多領(lǐng)域多變工程問(wèn)題的求解。Ansys基于Motif的菜單系統(tǒng)是用戶(hù)能夠通過(guò)對(duì)話(huà)框、下拉菜單和子菜單進(jìn)行數(shù)據(jù)輸入和功能選擇，為用戶(hù)使用Ansys提供“導(dǎo)航”。有限元法緒論(xln)共三十八頁(yè)21該程序最為顯著特點(diǎn)是加入交互式操作方式,即能在分析前用交互式圖形為驗(yàn)證模型(mxng)的幾何形狀,材料及邊界條件。分析之后又能立即用圖形交互的檢查分析結(jié)果。2 該程序提供一個(gè)循序漸進(jìn)的能力表,包括高級(jí)結(jié)構(gòu)非線(xiàn)性,電磁場(chǎng),計(jì)算流體力學(xué),設(shè)計(jì)優(yōu)化一般接觸面,自適應(yīng)網(wǎng)格(wn )劃分，大應(yīng)變,有限轉(zhuǎn)動(dòng)動(dòng)

9、能和參數(shù)化建模.一、Ansys軟件作為一個(gè)大型通用的有限元程序,其功能以為全 世界公認(rèn), Ansys軟件特點(diǎn)如下:有限元法緒論共三十八頁(yè)22軟件(run jin)的設(shè)計(jì)分析和優(yōu)化功能也可很方便地應(yīng)用到CAD系統(tǒng)生成的模型上。新增強(qiáng)的功能,使程序更靈活更好用,計(jì) 算速度更快。 3 為用戶(hù)使用本程序提供直觀(zhgun)的途徑界面,并配有完整的聯(lián) 機(jī)聲明及幫助系統(tǒng).供用戶(hù)完成深入應(yīng)用課題。有限元法緒論共三十八頁(yè)23二、 Ansys軟件提供了一個(gè)(y )不斷改進(jìn)的功能清單， 具體包括： 結(jié)構(gòu)高度非線(xiàn)性仿真 Ansys采用(ciyng)了牛頓-拉普森迭代求解，并為了增強(qiáng) 問(wèn)題的收斂性，提供了自適應(yīng)下降

10、、線(xiàn)性搜索、自動(dòng) 載荷步、二分法及弧長(zhǎng)法等一系列命令?？梢杂?jì)算由 大的位移、應(yīng)變及有限轉(zhuǎn)動(dòng)引起的結(jié)構(gòu)幾何非線(xiàn)性問(wèn) 題、與時(shí)間有關(guān)的材料非線(xiàn)性問(wèn)題以及接觸引起的狀 態(tài)非線(xiàn)性問(wèn)題。有限元法緒論共三十八頁(yè)24 熱分析 Ansys熱分析基于能量守恒原理的熱平衡方程， 用有限元法計(jì)算各節(jié)點(diǎn)的溫度，并導(dǎo)出其它熱物理參 數(shù)。包括熱傳導(dǎo)、熱對(duì)流及熱輻射三種傳導(dǎo)方式。此 外，還可以分析相變、有內(nèi)熱源、接觸(jich)熱阻等問(wèn)題。 熱分析用于計(jì)算一個(gè)系統(tǒng)或部件的溫度分布，如熱量 獲取或損失、熱梯度、熱流密度等。有限元法緒論(xln)共三十八頁(yè)25 電磁分析 Ansys可分析電磁場(chǎng)的多方面問(wèn)題，如電感、電 容、磁通

11、量密度、渦流、電場(chǎng)分布、磁力線(xiàn)、力、運(yùn) 動(dòng)效應(yīng)、電路和能量損失等?？捎糜谟行У胤治鱿旅?所列的各類(lèi)設(shè)備：電力發(fā)電機(jī)、變壓器、螺線(xiàn)管傳動(dòng)(chundng) 器、電動(dòng)機(jī)、磁成像系統(tǒng)、圖象顯示設(shè)備傳感器、磁 懸浮裝置、波導(dǎo)、開(kāi)關(guān)等。有限元法緒論(xln)共三十八頁(yè)26有限元法緒論(xln) 設(shè)計(jì)優(yōu)化 Ansys提供了兩種優(yōu)化方法，它們可以處理大多 數(shù)的優(yōu)化問(wèn)題。零階方法是一個(gè)很完美的處理方法， 可以很有效的處理大多數(shù)的工程問(wèn)題。一階方法基于 目標(biāo)函數(shù)對(duì)設(shè)計(jì)變量的敏感程度，因此更加適合于精 確的優(yōu)化分析。優(yōu)化中Ansys采用一系列的分析-評(píng)估 -修正的循環(huán)過(guò)程，這個(gè)過(guò)程重復(fù)(chngf)進(jìn)行直到所有設(shè)

12、計(jì) 滿(mǎn)足要求為止。共三十八頁(yè)27 計(jì)算流體(lit)動(dòng)力學(xué)分析 Ansys程序中的FLOTRAN CFD分析功能是一個(gè)用于 分析二維及三維流體流動(dòng)場(chǎng)的先進(jìn)的工具，可解決如 下的問(wèn)題： 作用于氣動(dòng)翼型上的升力和阻力； 超音速?lài)姽苤械牧鲌?chǎng)； 彎管中流體的復(fù)雜的三維流動(dòng)； 計(jì)算發(fā)動(dòng)機(jī)排氣系統(tǒng)中氣體(qt)壓力及溫度分布； 研究管路系統(tǒng)中熱的層化及分離； 使用混合流研究來(lái)估計(jì)熱沖擊的可能性； 用自然對(duì)流分析估計(jì)電子封裝芯片的熱性能； 對(duì)含有多種流體的熱交換器進(jìn)行研究。有限元法緒論共三十八頁(yè)28 利用Ansys參數(shù)設(shè)計(jì)語(yǔ)言（APDL）的擴(kuò)展宏命令功能 APDL有參數(shù)、數(shù)組參數(shù)、表達(dá)式和函數(shù)(hnsh)、

13、分支和 循環(huán)、重復(fù)和縮寫(xiě)、宏以及用戶(hù)程序等功能。有限元法緒論(xln)共三十八頁(yè)29 本世紀(jì)中葉，隨著生產(chǎn)力的發(fā)展，簡(jiǎn)單的實(shí)驗(yàn)?zāi)M(mn)和解析求解方法，已經(jīng)遠(yuǎn)遠(yuǎn)滿(mǎn)足不了工業(yè)發(fā)展的需要，人們?cè)趯ふ乙环N新的方法來(lái)更好的指導(dǎo)生產(chǎn)。有限元法緒論(xln)第四節(jié)有限元法的用途共三十八頁(yè)30有限元法緒論(xln) 隨著計(jì)算機(jī)的發(fā)展一種現(xiàn)代計(jì)算方法迅速發(fā)展起來(lái)，這就是有限元法。它是50年代首先在連續(xù)體力學(xué)領(lǐng)域(ln y)飛機(jī)結(jié)構(gòu)靜、動(dòng)態(tài)特性分析中應(yīng)用的一種有效的數(shù)值分析方法，隨后很快就廣泛地應(yīng)用于求解熱傳導(dǎo)、電磁場(chǎng)、流體力學(xué)等連續(xù)性問(wèn)題。目前被廣泛地應(yīng)用在航空、造船、機(jī)械、建筑、水利、鐵道、橋梁、石油、化

14、工、冶金、采礦、汽車(chē)等很多工程領(lǐng)域。共三十八頁(yè)31 下面我們將結(jié)合實(shí)例對(duì)有限元方法(fngf)的應(yīng)用進(jìn)行講解。有限元法緒論(xln)超彈性密封壓縮靜力分析用于靜態(tài)荷載.可以考慮結(jié)構(gòu)的非線(xiàn)性行為，例如：大變形、大應(yīng)變、應(yīng)力剛化、接觸、塑性、超彈及蠕變等.共三十八頁(yè)32有限元法緒論(xln)動(dòng)力分析包括質(zhì)量和阻尼效應(yīng)。模態(tài)分析，用于計(jì)算固有頻率和振型。諧響應(yīng)分析，用于確定結(jié)構(gòu)對(duì)輻值已知、頻率按正弦曲線(xiàn)變化的荷載的響應(yīng)。瞬態(tài)動(dòng)力分析，用于確定結(jié)構(gòu)對(duì)隨時(shí)間任意變化的載荷的響應(yīng)，而且可以包含非線(xiàn)性特性。其它結(jié)構(gòu)分析功能。譜分析。隨機(jī)振動(dòng)。特征值屈曲(q q)。子結(jié)構(gòu)，子模型。模態(tài)動(dòng)畫(huà)共三十八頁(yè)33有限元

15、法緒論(xln)進(jìn)行顯式動(dòng)力分析模擬以慣性力為主的大變形分析。用于模擬碰撞、擠壓和快速(kui s)成形等。汽車(chē)碰撞試驗(yàn)的沖擊分析共三十八頁(yè)34有限元法緒論(xln)熱分析用于確定物體(wt)中的溫度分布。熱分析考慮的物理量是：熱量的獲取和損失、溫度梯度、熱流量。可模擬三種熱傳遞方式：熱傳導(dǎo)、熱對(duì)流、熱輻射。 熨斗的瞬態(tài)熱分析共三十八頁(yè)35有限元法緒論(xln)電磁分析用于計(jì)算電磁設(shè)備中的磁場(chǎng)。穩(wěn)態(tài)和低頻電磁場(chǎng)分析模擬(mn)由直流電源，低頻交流電或低頻瞬時(shí)信號(hào)引起的磁場(chǎng)。共三十八頁(yè)36有限元法緒論(xln)高頻電磁場(chǎng)分析模擬電磁波傳播裝置。例如(lr)：微波及RF無(wú)源組件、波導(dǎo)、同軸連接器。電磁場(chǎng)分析中考慮的物理量是：S-參數(shù)、Q-因子、反射波損耗、電解質(zhì)及傳導(dǎo)損耗。同軸電纜中的電場(chǎng)共三十八頁(yè)37有限元法緒論(xln)流體動(dòng)力分析(CFD)用于確定流體中的流動(dòng)(lidng)狀態(tài)和溫度。模擬層流和湍流，可壓縮和不可壓縮流體，以及多組份流。應(yīng)用：航空航天、電子元件封裝、汽車(chē)設(shè)計(jì)。典型的物理量是：速度、壓力、溫度和對(duì)流換熱系數(shù)。導(dǎo)管中的流速足球上的壓力分布共三十八頁(yè)內(nèi)容摘要1。因此首先要對(duì)彈性體進(jìn)行必要的簡(jiǎn)化，再將彈性體劃分為有限個(gè)單元組成的離散體。Ansys基于Motif的菜單系統(tǒng)是用戶(hù)能夠通過(guò)對(duì)話(huà)框、下拉菜單和子