第一章_常用數(shù)值分析方法§4_有限差分法與有限單元法

1、 3:43 2022-3-51/45X.Z.LinW Y 3:43 2022-3-52/38X.Z.Lin 有限差分方法是數(shù)值計算中應用非常廣泛的一種有限差分方法是數(shù)值計算中應用非常廣泛的一種方法，是求解微分方程的主要方法之一。其實質(zhì)就是方法，是求解微分方程的主要方法之一。其實質(zhì)就是以有限差分代替無限微分、以差分代數(shù)方程代替微分以有限差分代替無限微分、以差分代數(shù)方程代替微分方程、以數(shù)值計算代替數(shù)學推導的過程，從而將連續(xù)方程、以數(shù)值計算代替數(shù)學推導的過程，從而將連續(xù)函數(shù)離散化。以有限的、離散的數(shù)值代替連續(xù)的函數(shù)函數(shù)離散化。以有限的、離散的數(shù)值代替連續(xù)的函數(shù)分布。分布。 4.1.1 概述概述W Y

2、 3:43 2022-3-53/38X.Z.Lin1、構成差分格式、構成差分格式 差分方程通常是一組數(shù)量較多的線性代數(shù)方程（即：線性方差分方程通常是一組數(shù)量較多的線性代數(shù)方程（即：線性方程組）。其求解方法有下列兩種程組）。其求解方法有下列兩種：（：（1）精確法，又稱直接法，精確法，又稱直接法，即消元法；（即消元法；（2）近似法，又稱間接法，即迭代法。）近似法，又稱間接法，即迭代法。2、求解差分方程、求解差分方程3、對所得到的數(shù)值解進行對所得到的數(shù)值解進行精度與收斂性分析和檢驗精度與收斂性分析和檢驗。 首先選擇網(wǎng)格布局、差分形式和步長；其次，以有限差分首先選擇網(wǎng)格布局、差分形式和步長；其次，以有

3、限差分代替無限微分，即以代替無限微分，即以 代替代替dx以差商以差商 代代替微商（導數(shù)）替微商（導數(shù)） ，以差分方程代替微分方程及邊界條件。，以差分方程代替微分方程及邊界條件。xxx122121yyyxxx21xxx dydxW Y 3:43 2022-3-54/38X.Z.Lin 建立差分方程是有限差分法的關鍵環(huán)節(jié)。導出差分建立差分方程是有限差分法的關鍵環(huán)節(jié)。導出差分方程的途徑可有兩種：方程的途徑可有兩種：（1）從微分方程出發(fā)，以泰勒級數(shù)截斷，從有限差分）從微分方程出發(fā)，以泰勒級數(shù)截斷，從有限差分的數(shù)學含義去建立有限差分和差分方程。的數(shù)學含義去建立有限差分和差分方程。（2）從由網(wǎng)格所劃分的單

4、元體的能量平衡分析出發(fā)、）從由網(wǎng)格所劃分的單元體的能量平衡分析出發(fā)、由積分方由積分方 程去建立差分方程，該方法又稱程去建立差分方程，該方法又稱單元體平衡法單元體平衡法。 兩種方法各具特色，但無論采取何種差分方程的推兩種方法各具特色，但無論采取何種差分方程的推導方法，在建立差分方程前，均需對所論區(qū)域進行離散導方法，在建立差分方程前，均需對所論區(qū)域進行離散化化 。W Y 3:43 2022-3-55/38X.Z.Lin合理選擇網(wǎng)格布局及步長合理選擇網(wǎng)格布局及步長 在實施有限差分法中首先在如圖在實施有限差分法中首先在如圖4.1所示的求解區(qū)域內(nèi)，將自所示的求解區(qū)域內(nèi)，將自變量變量x，y 分別沿分別沿x

5、，y軸方向的連續(xù)變化，離散為軸方向的連續(xù)變化，離散為x0，x1 ，x2，xn及及y0，y1 ，y2，yn個不連續(xù)點形成離散化網(wǎng)格；網(wǎng)格交點稱個不連續(xù)點形成離散化網(wǎng)格；網(wǎng)格交點稱為為結點結點(或節(jié)點或節(jié)點)，依次將結點編號，與區(qū)域自變量離散化相對應，依次將結點編號，與區(qū)域自變量離散化相對應，區(qū)域內(nèi)函數(shù)也將同時被離散化區(qū)域內(nèi)函數(shù)也將同時被離散化 。 離散化后各相鄰離散點之間的距離，或離散化單離散化后各相鄰離散點之間的距離，或離散化單元的長度稱為元的長度稱為步長步長，步長的大小可以是常量，也可以，步長的大小可以是常量，也可以是變量。是變量。網(wǎng)格的粗細與是否均勻網(wǎng)格的粗細與是否均勻，要根據(jù)求解區(qū)域物，

6、要根據(jù)求解區(qū)域物理場的實際分布和對結果所要求的精確度而定。理場的實際分布和對結果所要求的精確度而定。 一般說來，對均質(zhì)、形狀簡單且規(guī)則、物理量變一般說來，對均質(zhì)、形狀簡單且規(guī)則、物理量變化不劇烈的物體或求解精度要求不高時，可采用化不劇烈的物體或求解精度要求不高時，可采用等等步長、大步長步長、大步長，即采用均勻網(wǎng)格；而對形狀復雜、組，即采用均勻網(wǎng)格；而對形狀復雜、組分不同、物理量變化劇烈的物體，或求解精度要求較分不同、物理量變化劇烈的物體，或求解精度要求較高時，則采用高時，則采用小步長、變步長小步長、變步長。圖圖4.1 求解區(qū)域離散化求解區(qū)域離散化 W Y 3:43 2022-3-56/38X.

7、Z.Lin 另外，對一些較復雜的問題，在選擇網(wǎng)格與步長前，往往要對所論區(qū)另外，對一些較復雜的問題，在選擇網(wǎng)格與步長前，往往要對所論區(qū)域的物理場作出粗略估計，然后以較粗的網(wǎng)格、較大的步長計算出參考性域的物理場作出粗略估計，然后以較粗的網(wǎng)格、較大的步長計算出參考性物理場，根據(jù)這一參考性物理場再選擇合理的離散化網(wǎng)格物理場，根據(jù)這一參考性物理場再選擇合理的離散化網(wǎng)格。 離散化網(wǎng)格的布局，要根據(jù)所要求解的問題的性質(zhì)及求解離散化網(wǎng)格的布局，要根據(jù)所要求解的問題的性質(zhì)及求解要求確定。一般說來，有兩種方法：要求確定。一般說來，有兩種方法：（1）物理劃分法物理劃分法：這種方法是根據(jù)問題的物理特性劃分，如建筑：這

8、種方法是根據(jù)問題的物理特性劃分，如建筑物墻壁內(nèi)外層面磚、普通磚和內(nèi)灰泥層組成；若擬求各層界面壁物墻壁內(nèi)外層面磚、普通磚和內(nèi)灰泥層組成；若擬求各層界面壁溫，則離散化時應按不同材料組分劃分區(qū)域。溫，則離散化時應按不同材料組分劃分區(qū)域。圖圖4.2 扇形網(wǎng)格和三角形網(wǎng)格扇形網(wǎng)格和三角形網(wǎng)格 （2）幾何劃分法：幾何劃分法：以幾何區(qū)域以幾何區(qū)域形狀為依據(jù)來劃分，如對矩形區(qū)形狀為依據(jù)來劃分，如對矩形區(qū)域可采用矩形離散化網(wǎng)格，非矩域可采用矩形離散化網(wǎng)格，非矩形區(qū)域可采用三角形、四角形或形區(qū)域可采用三角形、四角形或其他形狀的網(wǎng)格，以適應溫度場其他形狀的網(wǎng)格，以適應溫度場分布的要求。分布的要求。W Y 3:43

9、2022-3-57/38X.Z.Lin將微分方程轉化為差分方程將微分方程轉化為差分方程 微分方程轉化為差分方程實際上就是以差分代替微微分方程轉化為差分方程實際上就是以差分代替微分、以差商代替微商的過程，是以有限小量去代替無限分、以差商代替微商的過程，是以有限小量去代替無限微量的近似化過程。微量的近似化過程。 方法：寫出微分方程中各微分與微商所對應方法：寫出微分方程中各微分與微商所對應的差分與差商形式，代入原微分方程即可。的差分與差商形式，代入原微分方程即可。W Y 3:43 2022-3-58/38X.Z.Lin某物理量的有限增量某物理量的有限增量 分分類類按按組組成成分分按按階數(shù)階數(shù)分分f2

10、一階差分一階差分二階差分二階差分n 階差分階差分fnf一階差分一階差分二階差分二階差分iiiffff1,2,1,21121() ()()2f if if if iiiiiiiifffffffffff 一階差分一階差分二階差分二階差分,1b iiifff2,111212() ()()2b ib ib ib iiiiiiiifffffffffff 一階差分一階差分二階差分二階差分1111,1122222iiiiiic iiifffffffff2,11,221 11 11 11 1112 22 22 22 2() () ()2c ic ic ic iiiiiiiifffffffffff 向前差分向前

11、差分向后差分向后差分中心差分中心差分W Y 3:43 2022-3-59/38X.Z.Lin函數(shù)的差分與自變量差分之比函數(shù)的差分與自變量差分之比 用差分代替微分方程中的微分，用差商代替微分方程中的用差分代替微分方程中的微分，用差商代替微分方程中的微商，即可將微分方程轉化為差分方程。微商，即可將微分方程轉化為差分方程。 差分方程通常是一個線性方程組，利用以前介紹的直接法差分方程通常是一個線性方程組，利用以前介紹的直接法（消元法）或間接法（迭代法）即可解之，從而得到原微分方（消元法）或間接法（迭代法）即可解之，從而得到原微分方程的解。程的解。W Y 3:43 2022-3-510/38X.Z.Li

12、n 有限單元法（又稱為有限元素法，簡稱有限元法），是有限單元法（又稱為有限元素法，簡稱有限元法），是20世紀世紀50年代初年代初才出現(xiàn)的才出現(xiàn)的 一種新的數(shù)值分析方法，最早應用于一種新的數(shù)值分析方法，最早應用于航空航天領域航空航天領域，主要用于，主要用于力學力學與結構分析與結構分析中，中， 20 世紀世紀 70 年以來被應用到年以來被應用到傳熱學傳熱學計算中。與有限差分法相計算中。與有限差分法相比較，有限元法的準確性和穩(wěn)定性都比較好，且由于其單元的靈活性，使它比較，有限元法的準確性和穩(wěn)定性都比較好，且由于其單元的靈活性，使它更適應于數(shù)值求解非線性熱傳導問題以及具有不規(guī)則幾何形狀與邊界，特別更適

13、應于數(shù)值求解非線性熱傳導問題以及具有不規(guī)則幾何形狀與邊界，特別是要求同時得到熱應力場的各種復雜導熱問題；有限元法在傳熱學中的應用是要求同時得到熱應力場的各種復雜導熱問題；有限元法在傳熱學中的應用正處于開拓與發(fā)展階段，迄今為止，其應用已波及熱傳導、正處于開拓與發(fā)展階段，迄今為止，其應用已波及熱傳導、 對流傳熱及換熱對流傳熱及換熱器設計與計算。器設計與計算。4.1.1 概述概述 從從“有限元有限元”的名字出現(xiàn)到今天，經(jīng)歷了幾十年的發(fā)展，其基本理論已的名字出現(xiàn)到今天，經(jīng)歷了幾十年的發(fā)展，其基本理論已經(jīng)日趨完善，復雜非線性問題的各種算法得到很大的發(fā)展，并且在工程領域經(jīng)日趨完善，復雜非線性問題的各種算法

14、得到很大的發(fā)展，并且在工程領域（如：結構力學、熱傳導、電磁場、流體力學等連續(xù)域問題）得到廣范的應（如：結構力學、熱傳導、電磁場、流體力學等連續(xù)域問題）得到廣范的應用。用。 W Y 3:43 2022-3-511/38X.Z.Lin （1）假想把連續(xù)系統(tǒng)（包括桿系，連續(xù)體，連續(xù)介質(zhì)）假想把連續(xù)系統(tǒng)（包括桿系，連續(xù)體，連續(xù)介質(zhì)）分割分割成數(shù)目成數(shù)目有限的單元，單元之間只在數(shù)目有限的指定點（稱為結點）處相互連接，有限的單元，單元之間只在數(shù)目有限的指定點（稱為結點）處相互連接，構成一個構成一個單元集合體單元集合體來代替原來的連續(xù)系統(tǒng)。在結點上引進等效載荷來代替原來的連續(xù)系統(tǒng)。在結點上引進等效載荷（或邊

15、界條件），代替實際作用于系統(tǒng)上的外載荷（或（或邊界條件），代替實際作用于系統(tǒng)上的外載荷（或 邊界條件）。邊界條件）。 （2）對每個單元由）對每個單元由分塊近似分塊近似的思想，按一定的規(guī)則（由力學關系或的思想，按一定的規(guī)則（由力學關系或選擇一個簡單函數(shù)）建立待求選擇一個簡單函數(shù)）建立待求未知量與結點相互作用（力）之間的關系未知量與結點相互作用（力）之間的關系（力位移、熱量溫度、電壓電流等）。（力位移、熱量溫度、電壓電流等）。 （3）把所有單元的這種特性關系按一定的條件（變形協(xié)調(diào)條件、連）把所有單元的這種特性關系按一定的條件（變形協(xié)調(diào)條件、連續(xù)條件或變分原理及能量原理）續(xù)條件或變分原理及能量原理）

16、集合集合起來，引入邊界條件，構成一組以起來，引入邊界條件，構成一組以結點變量（位移、溫度、電壓等）為未知量的結點變量（位移、溫度、電壓等）為未知量的代數(shù)方程組代數(shù)方程組， 解之就可得解之就可得到有限個節(jié)點處的待求變量到有限個節(jié)點處的待求變量 。 可見：有限元法實質(zhì)上是把具有無限個自由度的連續(xù)系可見：有限元法實質(zhì)上是把具有無限個自由度的連續(xù)系統(tǒng)，理想化為只有有限個自由度的單元集合體，使問題轉化統(tǒng)，理想化為只有有限個自由度的單元集合體，使問題轉化為適合于數(shù)值求解的結構型問題。為適合于數(shù)值求解的結構型問題。 W Y 3:43 2022-3-512/38X.Z.Lin圖圖4.3v為沖壓位移為沖壓位移W

17、 Y 3:43 2022-3-513/38X.Z.Lin圖圖4.4W Y 3:43 2022-3-514/38X.Z.Lin直接剛度法直接剛度法 在剛提出有限元法的時候采用的是直接剛度法它源在剛提出有限元法的時候采用的是直接剛度法它源于結構分析的剛度法。因剛度法只能處理于結構分析的剛度法。因剛度法只能處理些比較簡單的些比較簡單的實際問題，現(xiàn)在已很少使用了，但它對我們理解和明確有實際問題，現(xiàn)在已很少使用了，但它對我們理解和明確有限元法的一些物理概念是很有幫助的。所以我們首先通過限元法的一些物理概念是很有幫助的。所以我們首先通過一個例子來介紹直接剛度法，同時說明有限元法求解的一一個例子來介紹直接剛

18、度法，同時說明有限元法求解的一般步驟。般步驟。 W Y 3:43 2022-3-515/38X.Z.Lin 考慮一個變截面桿，如圖考慮一個變截面桿，如圖4.5所所示，桿的一端固定，另一端承受示，桿的一端固定，另一端承受 P1000N的載荷桿的頂部寬的載荷桿的頂部寬w1=2cm， 桿的底部寬桿的底部寬w2=1cm ，桿，桿的厚度的厚度t=0.125cm，長度，長度L 10cm、桿的彈性模量桿的彈性模量E10.4106MPa。試分析該桿沿長度方向不同位置的試分析該桿沿長度方向不同位置的變形情況，假設桿的質(zhì)量可以忽略變形情況，假設桿的質(zhì)量可以忽略不計。不計。圖圖4.5 受軸向載荷的變截面桿受軸向載荷

19、的變截面桿有限元基本概念與一般步驟有限元基本概念與一般步驟W Y 3:43 2022-3-516/38X.Z.Lin圖圖4.6 將桿劃分為單元和結點將桿劃分為單元和結點 先將求解的問題分解為結點和單元。為簡單起見，將桿劃分成五個結先將求解的問題分解為結點和單元。為簡單起見，將桿劃分成五個結點和四個單元（如圖點和四個單元（如圖4.6所示）。所示）。 給定的變截而桿簡化為四個獨立的部分，每部分的截面面積恒定（為給定的變截而桿簡化為四個獨立的部分，每部分的截面面積恒定（為組成該單元的兩個結點處的面積的平均值）。組成該單元的兩個結點處的面積的平均值）。(1) 求解域離散化求解域離散化 W Y 3:43

20、 2022-3-517/38X.Z.Lin圖圖4.7 具有均勻截面的固體單元在力具有均勻截面的固體單元在力F作用下的變形作用下的變形 單元中的平均應力為：單元中的平均應力為： 平均正應變?yōu)椋浩骄龖優(yōu)椋?在彈性范圍內(nèi)，由在彈性范圍內(nèi)，由Hooke定律定律E可得：可得： lAEklkFeqeq注意：上式與線性彈簧等式注意：上式與線性彈簧等式F=kx相似。相似。 因此可因此可用彈簧的變形來模擬固態(tài)單元的變形，從而得：用彈簧的變形來模擬固態(tài)單元的變形，從而得：（keq稱為等價剛度）稱為等價剛度）W Y 3:43 2022-3-518/38X.Z.Lin圖圖4.8 結點受力分析結點受力分析 將上述結

21、果擴展到整個變截面桿上，則：桿可將上述結果擴展到整個變截面桿上，則：桿可以用一個由四個彈簧（五個結點）串聯(lián)組成的模型以用一個由四個彈簧（五個結點）串聯(lián)組成的模型來表示，每個單元模型的彈性行為可以用等價線性來表示，每個單元模型的彈性行為可以用等價線性彈簧來表示：彈簧來表示： 根據(jù)靜態(tài)平衡的要求：作用在每個結點上的力的根據(jù)靜態(tài)平衡的要求：作用在每個結點上的力的總和為零，因此得到如下方程：總和為零，因此得到如下方程： )(2)()()(1111iiiiiiavgiieqiuulEAAuulEAuukFlEAAkiieq2)(1等價單元剛度：W Y 3:43 2022-3-519/38X.Z.Lin

22、將這些方程將這些方程進行變化可得進行變化可得到：到： 矩陣形式為：矩陣形式為： 在載荷矩陣在載荷矩陣中區(qū)分施加力中區(qū)分施加力與反作用力也與反作用力也是必要的。因是必要的。因此矩陣式又可此矩陣式又可寫為寫為 ：W Y 3:43 2022-3-520/38X.Z.Lin 反作用力矩陣反作用力矩陣 剛度矩陣剛度矩陣 位移矩陣位移矩陣 載荷矩陣載荷矩陣 對我們所討論的問對我們所討論的問題，因為桿的頂端固定，題，因為桿的頂端固定，結點結點1的位移應力的位移應力0，即，即 u10，將此邊界條件用于，將此邊界條件用于上式，得到如下矩陣：上式，得到如下矩陣： 解此方程組即可得到各結點的位移解此方程組即可得到各

23、結點的位移 W Y 3:43 2022-3-521/38X.Z.Lin討論討論 一般單元的單元剛度矩陣，并討論總體剛度矩陣的集成。一般單元的單元剛度矩陣，并討論總體剛度矩陣的集成。 先以一個單元（含兩個結點）作為研究對象，其傳輸力先以一個單元（含兩個結點）作為研究對象，其傳輸力可用圖可用圖4.9表示。表示。圖圖4.9 單元傳輸力單元傳輸力FiFiFi+1Fi+1=)()(111iieqiiieqiuukFuukF確定每個單元的方程：確定每個單元的方程： 單元剛單元剛度矩陣度矩陣W Y 3:43 2022-3-522/38X.Z.Lin集成單元：集成單元： 首先寫出各單元剛度矩陣及其在總體剛度矩

24、陣中的位置如下：首先寫出各單元剛度矩陣及其在總體剛度矩陣中的位置如下：單元單元 1單元單元 2單元單元 3單元單元 4組裝（相加），即得總體剛度矩陣：組裝（相加），即得總體剛度矩陣：W Y 3:43 2022-3-523/38X.Z.Lin施加邊界條件和載荷施加邊界條件和載荷 桿的頂端固定，應滿足邊界條件桿的頂端固定，應滿足邊界條件 u1=0載荷載荷P 施加施加在結點在結點5上。應用這些條件得到下式上。應用這些條件得到下式 ： 剛度矩陣剛度矩陣 位移矩陣位移矩陣 載荷矩陣載荷矩陣 有限元分析的一般格式有限元分析的一般格式（固體力學問題）（固體力學問題）解此方程組即可得到各結點的位移解此方程組即

25、可得到各結點的位移 W Y 3:43 2022-3-524/38X.Z.Lin桿的截面沿桿的截面沿 y 軸的變化可以描述為：軸的變化可以描述為： 計算得：計算得：1()2iieqAA Ekl由等價單元剛度（剛度系數(shù)）公式：由等價單元剛度（剛度系數(shù)）公式：計算得：計算得：各單元剛度矩陣為：各單元剛度矩陣為：W Y 3:43 2022-3-525/38X.Z.Lin解此矩陣方程（線性方程組），即可得到各結點的位移值：解此矩陣方程（線性方程組），即可得到各結點的位移值：將單元剛度矩陣組合，得到總體剛度矩陣：將單元剛度矩陣組合，得到總體剛度矩陣： 應用邊界條件并施加載荷，得到：應用邊界條件并施加載荷，

26、得到： W Y 3:43 2022-3-526/38X.Z.Lin因為結點的位移已知、上式也可直接從應力因為結點的位移已知、上式也可直接從應力應變關系得到：應變關系得到： F 對分析結果進行后處理，我們可以得到其他一些信息。對分析結果進行后處理，我們可以得到其他一些信息。如：每個單元中的平均正應力。這些值可以從下式得到：如：每個單元中的平均正應力。這些值可以從下式得到： 對本例，可計算出各單元的平均正應力值如下：對本例，可計算出各單元的平均正應力值如下： W Y 3:43 2022-3-527/38X.Z.Lin 有限元法程序總體可分為三個組成部分：有限元法程序總體可分為三個組成部分：前處理部

27、分、有前處理部分、有限元分析本體程限元分析本體程 序、后處理部分序、后處理部分。 有限元法的實現(xiàn)必須通過計算機，全部有限元法的計算原有限元法的實現(xiàn)必須通過計算機，全部有限元法的計算原理和數(shù)值方法集中反映在有限元法的程序中，因此有限元法的理和數(shù)值方法集中反映在有限元法的程序中，因此有限元法的程序極為重要。它應具有程序極為重要。它應具有分析準確可靠、計算效率高、使用方分析準確可靠、計算效率高、使用方便、易于擴充和修改便、易于擴充和修改等特點：等特點：（1）特點）特點（2）結構）結構W Y 3:43 2022-3-528/38X.Z.Lin 對于一個實際的工程問題，離散模型的數(shù)據(jù)文件十分龐大，靠人對

28、于一個實際的工程問題，離散模型的數(shù)據(jù)文件十分龐大，靠人工處理和生成一般是不可能的。為了解決這一問題，有限元分析程序必工處理和生成一般是不可能的。為了解決這一問題，有限元分析程序必須有前處理程序。須有前處理程序。前處理程序前處理程序是根據(jù)使用者提供的對計算模型外形及網(wǎng)是根據(jù)使用者提供的對計算模型外形及網(wǎng)格要求的簡單數(shù)據(jù)描述，自動或半自動地生成離散模型的數(shù)據(jù)文件，并格要求的簡單數(shù)據(jù)描述，自動或半自動地生成離散模型的數(shù)據(jù)文件，并要生成網(wǎng)格圖供使用者檢查和修改。這部分程序的功能很大程度上決定要生成網(wǎng)格圖供使用者檢查和修改。這部分程序的功能很大程度上決定了程序使用的方便性。了程序使用的方便性。 有限元分

29、析本體程序有限元分析本體程序是有限元分析程序的核心，它根據(jù)離散模型是有限元分析程序的核心，它根據(jù)離散模型的數(shù)據(jù)文件進行有限元分析；有限元分析的原理和采用的數(shù)值方法集的數(shù)據(jù)文件進行有限元分析；有限元分析的原理和采用的數(shù)值方法集中于此，因此它是有限元分忻準確可靠的關鍵。選用計算方法的合理中于此，因此它是有限元分忻準確可靠的關鍵。選用計算方法的合理與否決定了有限元分析程序的計算效率和結果的精度及可靠性與否決定了有限元分析程序的計算效率和結果的精度及可靠性 。 W Y 3:43 2022-3-529/38X.Z.Lin 同樣，有限元分析程序的計算結果也是針對離散模型得到的。例同樣，有限元分析程序的計算

30、結果也是針對離散模型得到的。例如靜力平衡問題可以得到離散模型各結點的位移、各單元的應力等，如靜力平衡問題可以得到離散模型各結點的位移、各單元的應力等，輸出的文本文件量很輸出的文本文件量很 大，但卻不易得到所分析對象的全貌，例如位移大，但卻不易得到所分析對象的全貌，例如位移哪里最大、應力集中發(fā)生在什么部位以及變化趨勢如何等。因此一個哪里最大、應力集中發(fā)生在什么部位以及變化趨勢如何等。因此一個使用方便的有限元分析程序木不僅要有可供選擇輸出內(nèi)容的文本文件，使用方便的有限元分析程序木不僅要有可供選擇輸出內(nèi)容的文本文件，還需有結果的圖形顯示，如位移圖、等應力線圖或截面應力分布圖等；還需有結果的圖形顯示，

31、如位移圖、等應力線圖或截面應力分布圖等；這部分程序稱這部分程序稱后處理程序后處理程序。與前處理程序相似，后處理程序?qū)Τ绦蚴?。與前處理程序相似，后處理程序?qū)Τ绦蚴褂玫姆奖阈杂信e足輕重的作用。用的方便性有舉足輕重的作用。 有限元分析程序的三個組成部分對于一個較好的用于實際問題分有限元分析程序的三個組成部分對于一個較好的用于實際問題分析的有限元程序來說，前后處理的析的有限元程序來說，前后處理的程序量程序量常常超出有限元分析的本體常常超出有限元分析的本體程序，前后處理功能越強，程序的使用就越方便。有限元分析程序中程序，前后處理功能越強，程序的使用就越方便。有限元分析程序中前后處理程序一般可占到全部程序

32、條數(shù)的前后處理程序一般可占到全部程序條數(shù)的 2/3 4/5。有的近期發(fā)展的。有的近期發(fā)展的通用程序更注重程序的通用程序更注重程序的“包裝包裝”和使用功能。有限元分析本體程序以和使用功能。有限元分析本體程序以外部分的比例更高。外部分的比例更高。 W Y 3:43 2022-3-530/38X.Z.Lin ANSYS（Analysis System）世界著名力學分析專家、匹茲堡大學J.Swansan教授創(chuàng)立的SASI（Swansan Analysis System Inc.）的大型通用有限元分析軟件世界最權威的有限元產(chǎn)品。 SAP （Structural Analysis Program）美國加州

33、大學伯克利分校M.J.Wilson教授的線性靜、動力學結構分析程序。 NASTRAN（NASA Structural Analysis ）美國國家航空宇航局（NASA）的結構分析程序。 ADINA（A Finite Element Program for Automatic Dynamic Incremental Nonlinear Analysis ）美國麻省理工學院機械工程系的自動動力增量非線性分析有限元程序。 IDEAS（Integrate Design Engineering Analysis System）美同SDRC公司的機械通用軟件，集成化設計工程分析系統(tǒng)。集設計、分析、數(shù)控加工

34、、塑料模具設計和測試數(shù)據(jù)為一體的工作站用軟件。 AGOBJ美國AGOBJ 公司在SAP5和ADINA 有限元分析程序的基礎上針對微機平臺開發(fā)的通用有限元分忻系統(tǒng)。W Y 3:43 2022-3-531/38X.Z.Lin ANSYS軟件是軟件是20世紀世紀70年代由美國年代由美國ANSYS公司開發(fā)的一套功能公司開發(fā)的一套功能強大的有限元通用分析程序，具有強大的前處理、求解和后處理能強大的有限元通用分析程序，具有強大的前處理、求解和后處理能力，目前廣泛應用于航空、航天、汽車、力，目前廣泛應用于航空、航天、汽車、 船舶、土木、電子、機械船舶、土木、電子、機械等科學應用領域。等科學應用領域。ANSY

35、S把把CAD、CAE、CAM技術集成于一身，技術集成于一身，可以滿足用戶從設計、計算、制造全過程的使用要求。程序的開放可以滿足用戶從設計、計算、制造全過程的使用要求。程序的開放性是衡量其能力的一個重要標準，性是衡量其能力的一個重要標準， ANSYS自帶的自帶的APDL語言可以直語言可以直接以接以ANSYS程序為平臺，為用戶提供二次開發(fā)，其語法簡單，易學程序為平臺，為用戶提供二次開發(fā)，其語法簡單，易學易用，是用戶進行深入研究的強大輔助工具。易用，是用戶進行深入研究的強大輔助工具。 ANSYS程序可以采用程序可以采用APDL命令流輸入建模分析，也可以采用命令流輸入建模分析，也可以采用GUI交互式輸

36、入方式交互式輸入方式 建模。前一種方式最適合科研人員對同一個產(chǎn)建模。前一種方式最適合科研人員對同一個產(chǎn)品進行參數(shù)分析，修改方便，耗用計算機資源少，而且高級使用者品進行參數(shù)分析，修改方便，耗用計算機資源少，而且高級使用者還可以結合行業(yè)的特殊問題進行擴展開發(fā)，滿足用戶的特殊需求。還可以結合行業(yè)的特殊問題進行擴展開發(fā)，滿足用戶的特殊需求。 后一種方式操作簡單、直觀，適合于初學者；同時，對已確定的模后一種方式操作簡單、直觀，適合于初學者；同時，對已確定的模型進行分析快捷高效。型進行分析快捷高效。W Y 3:43 2022-3-532/38X.Z.Lin1結構靜力分析結構靜力分析 用來求解穩(wěn)態(tài)外載引起的

37、系統(tǒng)或部件的位移、應變、應力和用來求解穩(wěn)態(tài)外載引起的系統(tǒng)或部件的位移、應變、應力和力。靜力分析很適合求解慣性和阻尼對結構的影響并不顯著的問力。靜力分析很適合求解慣性和阻尼對結構的影響并不顯著的問題，如確定結構中應力集中現(xiàn)象。題，如確定結構中應力集中現(xiàn)象。ANSYS程序中的靜力分析不僅程序中的靜力分析不僅可以進行線性分析，而且可以進行非線性分析，如塑性、蠕變、可以進行線性分析，而且可以進行非線性分析，如塑性、蠕變、膨脹、大膨脹、大 變形、大應變及接觸分析。變形、大應變及接觸分析。2結構動力學分析結構動力學分析 結構動力學分析用來求解隨時間變化的載荷對結構或部件的結構動力學分析用來求解隨時間變化的

38、載荷對結構或部件的影響。與靜力分析不同，影響。與靜力分析不同， 動力學分析要考慮隨時間變化的力載荷動力學分析要考慮隨時間變化的力載荷以及對阻尼和慣性的影響。以及對阻尼和慣性的影響。ANSYS可進行的結構動力學分析類型可進行的結構動力學分析類型包括：瞬態(tài)動力學分析、模態(tài)分析、諧波響應分析及隨機振動響包括：瞬態(tài)動力學分析、模態(tài)分析、諧波響應分析及隨機振動響應分析。應分析。 W Y 3:43 2022-3-533/38X.Z.Lin3結構非線性分析結構非線性分析 結構非線性導致結構或部件的響應隨外載荷不成比例的變化。結構非線性導致結構或部件的響應隨外載荷不成比例的變化。ANSYS可求解靜態(tài)和瞬態(tài)非線

39、性間題，包括材料非線性、幾何非可求解靜態(tài)和瞬態(tài)非線性間題，包括材料非線性、幾何非線性、單元非線性線性、單元非線性3種。種。 （1）幾何非線性）幾何非線性 ：主要包括大變形、大應變、應力強化、旋轉軟主要包括大變形、大應變、應力強化、旋轉軟化、非線性屈曲等問題的分析化、非線性屈曲等問題的分析（2）材料非線性：）材料非線性：主要包括彈塑性（雙線性隨動硬化、雙線性各主要包括彈塑性（雙線性隨動硬化、雙線性各向同性硬化、多線性隨動硬化、多線性各向同性硬化、非線性隨動向同性硬化、多線性隨動硬化、多線性各向同性硬化、非線性隨動硬化、非線性各向同性硬化、非均勻各向異性、速率相關塑性、硬化、非線性各向同性硬化、非

40、均勻各向異性、速率相關塑性、 復合彈塑性等）、非線性彈性（分段線性彈性）、超彈性（各種橡復合彈塑性等）、非線性彈性（分段線性彈性）、超彈性（各種橡膠、膠、Ni-Ti合金等）、合金等）、 粘彈性各種玻璃、塑料等）、粘塑性（高粘彈性各種玻璃、塑料等）、粘塑性（高溫金屬等）、蠕變、膨脹（核材料、巖土和混凝土材料等）。溫金屬等）、蠕變、膨脹（核材料、巖土和混凝土材料等）。（3）單元非線性）單元非線性 ：W Y 3:43 2022-3-534/38X.Z.Lin4運動學分析運動學分析 ANSYS 可以分析大型三維柔體運動。當運動的積累起主要影響可以分析大型三維柔體運動。當運動的積累起主要影響作用時，可使

41、用這些功能分析復雜結構在空間的運動特性，并確定結作用時，可使用這些功能分析復雜結構在空間的運動特性，并確定結構中由此產(chǎn)生的應力、應變和變形。構中由此產(chǎn)生的應力、應變和變形。 6電磁場分析電磁場分析 主要用于電磁場問題的分析，如電感、電容、磁通量密度、渦流、主要用于電磁場問題的分析，如電感、電容、磁通量密度、渦流、電場分布、磁力線分布、力、運動效應、電路和能量損失等，還可以電場分布、磁力線分布、力、運動效應、電路和能量損失等，還可以應用于應用于 調(diào)節(jié)器、發(fā)電機、磁體、電解槽及無損檢測裝置等的設計和分調(diào)節(jié)器、發(fā)電機、磁體、電解槽及無損檢測裝置等的設計和分析領域析領域 。 5熱分析熱分析 ANSYS

42、軟件可處理熱傳遞的軟件可處理熱傳遞的3種基本類型：傳質(zhì)、對流和輻射，對種基本類型：傳質(zhì)、對流和輻射，對熱傳遞的熱傳遞的3種類型均可進行穩(wěn)態(tài)和瞬態(tài)、線性和非線性分析。種類型均可進行穩(wěn)態(tài)和瞬態(tài)、線性和非線性分析。ANSYS熱分析還具有可以模擬材料固化和熔解過程中的相變分析能力，以及熱分析還具有可以模擬材料固化和熔解過程中的相變分析能力，以及模擬熱與結構應力之間的熱結構耦合分析能力。模擬熱與結構應力之間的熱結構耦合分析能力。 W Y 3:43 2022-3-535/38X.Z.Lin7流體動力學分析流體動力學分析 ANSYS流體單元能進行流體動力學分析，分析類型可以為瞬態(tài)或穩(wěn)流體單元能進行流體動力學

43、分析，分析類型可以為瞬態(tài)或穩(wěn)態(tài)，分析結果可以是每個節(jié)點的壓力和通過每個單元的流率。并且可以態(tài)，分析結果可以是每個節(jié)點的壓力和通過每個單元的流率。并且可以利用后處理功能產(chǎn)生壓力、流率和溫度分布的圖形顯示。另外，還可以利用后處理功能產(chǎn)生壓力、流率和溫度分布的圖形顯示。另外，還可以使用三維表面效應單元和熱流管單元模擬結構的流體繞流，并包括對使用三維表面效應單元和熱流管單元模擬結構的流體繞流，并包括對流換熱效應。流換熱效應。9壓電分析壓電分析 用于分析二維或二維結構對用于分析二維或二維結構對AC(交流交流)、DC(直流直流)或任意隨時間變化或任意隨時間變化的電流或機械裁荷的響應。這種分析類型可用于換熱

44、器、振蕩器、諧振的電流或機械裁荷的響應。這種分析類型可用于換熱器、振蕩器、諧振器、麥克風等部件及其他電子設備的結構動態(tài)性能分析，可進行器、麥克風等部件及其他電子設備的結構動態(tài)性能分析，可進行 4 種類種類型的分析：靜態(tài)分析、模態(tài)分析、諧波響應分析、瞬態(tài)響應分析。型的分析：靜態(tài)分析、模態(tài)分析、諧波響應分析、瞬態(tài)響應分析。8、聲場分析、聲場分析 軟件的聲學功能用來研究含流體的介質(zhì)中的聲波傳播，或分析浸在軟件的聲學功能用來研究含流體的介質(zhì)中的聲波傳播，或分析浸在流體中的固體結構的動態(tài)特性。這些功能可用來確定音響話筒的頻率響流體中的固體結構的動態(tài)特性。這些功能可用來確定音響話筒的頻率響應，研究音樂大廳的聲場強度分布，應，研究音樂大廳的聲場強度分布， 或預測水對振動船體的阻尼效應?；蝾A測水對振動船體的阻尼效應。W Y 3:43 2022-3-536/38X.Z.Lin1多物理場耦合分析多物理場耦合分析 考慮兩個或多個物理場之間的相互作用對于正確求解很多實際工程考慮兩個或多個物理場