第五章偏微分方程的有限元法演示

（優(yōu)選）第五章偏微分方程的有限元法目前一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@2/104第五章偏微分方程的有限元法有限元法（FEA，F(xiàn)initeElementAnalysis，F(xiàn)EM）有限元法的基本思想是用較簡單的問題代替復雜問題，然后再對簡單問題進行求解的數(shù)值計算方法。有限元法將求解域看成是由許多被稱為有限元的小的互連子域組成，對每一單元假定一個較簡單的近似解，然后推導求解這個域總的滿足條件，從而得到問題的解。這個解不是準確解，而是近似解。有限元不僅計算精度高，而且能適應(yīng)各種復雜形狀，因而成為行之有效的數(shù)值計算方法。有限元法于上世紀50年代首先在力學領(lǐng)域-----飛機結(jié)構(gòu)的靜、動態(tài)特性分析中得到應(yīng)用，隨后很快廣泛的應(yīng)用于求解熱傳導、電磁場、流體力學等連續(xù)性問題。有限元法主要用于求解拉普拉斯方程和泊松方程所描述的各類物理場中。目前二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@3/104第五章偏微分方程的有限元法有限元法---變分原理基于變分原理的有限元法是逼近論、偏微分方程、變分與泛函分析的巧妙結(jié)合?；谧兎衷淼挠邢拊ㄒ宰兎衷頌榛A(chǔ)，把所要求解的微分方程定解問題，首先轉(zhuǎn)化為相應(yīng)的變分問題，即泛函求極值問題；它將求解域看成是由許多稱為有限元的小的互連子域組成，然后利用剖分插值，對每一單元假定一個合適的(較簡單的）近似解，把離散化的變分問題轉(zhuǎn)化為普通多元函數(shù)的極值問題，然后推導求解這個域總的滿足條件(邊界條件），即最終歸結(jié)為一組多元的代數(shù)方程組，求解代數(shù)方程組，就得到待求邊值問題的數(shù)值解。目前三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@4/104第五章偏微分方程的有限元法有限元法--加權(quán)余數(shù)法自從1969年以來，某些學者在流體力學中應(yīng)用加權(quán)余數(shù)法中的迦遼金法或最小二乘法等同樣獲得了有限元方程，因而有限元法可應(yīng)用于以任何微分方程所描述的各類物理場中，而不再要求這類物理場和泛函的極值問題有所聯(lián)系。加權(quán)余數(shù)法的核心思想是：近似解與解析解相比會存在誤差R，但是可以通過一個準則使R盡量小，求解這個等式，就可以得到待定常數(shù)的值，也就得到了近似解。目前四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@5/104第五章偏微分方程的有限元法有限元法特點有限元法的物理意義直觀明確，理論完整可靠。因為變分原理描述了支配物理現(xiàn)象的物理學中的最小作用原理（如力學中的最小勢能原理）。優(yōu)異的解題能力。有限元法對邊界幾何形狀復雜以及媒質(zhì)物理性質(zhì)變異等復雜物理問題求解上，有突出優(yōu)點：①不受幾何形狀和媒質(zhì)分布的復雜程度限制。②不必單獨處理第二、三類邊界條件。③離散點配置比較隨意，通過控制有限單元剖分密度和單元插值函數(shù)的選取，可以充分保證所需的數(shù)值計算精度。目前五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@6/1045.1泛函與變分原理數(shù)學上，通常自變量與因變量間的關(guān)系稱為函數(shù)，而泛函則是函數(shù)集合的函數(shù)，也就是函數(shù)的函數(shù)，即自變量為函數(shù)，而不是變量。5.1.1泛函的定義泛函通常是指一種定義域為函數(shù)，而值域為實數(shù)的“函數(shù)”。設(shè)C是函數(shù)的集合，B是實數(shù)集合。如果對C中的任一元素y(x)，在B中都有一個元素J與之對應(yīng)，則稱J為y(x)的泛函，記為J[y(x)]。目前六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@7/1045.1泛函與變分原理例5.1.1質(zhì)點在重力作用下，沿一條光滑的從A點到B點的曲線運動，如圖所示。求下落時間最短的曲線。曲線上任一小段線元長度為：ABxyOx0x1捷線問題目前七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@8/1045.1泛函與變分原理線元處的質(zhì)點速度為ABxyOx0x1ds線元下落時間為從A點到B點的下落時間為目前八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@9/1045.1泛函與變分原理5.1.2函數(shù)的變分設(shè)y(x)是泛函J定義域內(nèi)任一函數(shù)，如果y(x)變化為新函數(shù)Y(x)，且Y(x)屬于泛函J的定義域，則Y(x)與y(x)之差為函數(shù)y(x)的變分。變分δy是x的函數(shù)，它不同于函數(shù)的增量Δy。性質(zhì)：函數(shù)求導與求變分可以交換次序目前九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@10/1045.1泛函與變分原理5.1.3泛函的變分定義最簡泛函F(x,y,y’)稱為泛函的“核函數(shù)”泛函的變分最簡泛函:核函數(shù)只包含自變量x、未知函數(shù)y(x)以及導數(shù)y’(x)目前十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@11/1045.1泛函與變分原理利用二元函數(shù)的泰勒展開目前十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@12/1045.1泛函與變分原理其中分別稱為泛函的一階變分和二階變分。目前十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@13/1045.1泛函與變分原理泛函取極值的必要條件：一階變分為零性質(zhì)：對于最簡泛函，變分運算可以與積分、微分運算交換次序目前十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@14/1045.1泛函與變分原理5.1.4泛函的極值問題泛函的一階變分利用1泛函的極值問題的間接解法轉(zhuǎn)化為微分方程：歐拉方程目前十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@15/1045.1泛函與變分原理對于駐定問題，兩邊界固定這就是最簡泛函的歐拉方程，等價于泛函取極值的必要條件。把變分問題轉(zhuǎn)化微分方程的定解問題（邊值問題）來求解。目前十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@16/1045.1泛函與變分原理對于例求下落時間最短的軌跡利用最簡泛函的歐拉方程。目前十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@17/1045.1泛函與變分原理代入歐拉方程目前十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@18/1045.1泛函與變分原理變換得到進一步化簡得到積分目前十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@19/1045.1泛函與變分原理做變量替換得而目前十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@20/1045.1泛函與變分原理對上式積分得到這樣就得到了下落時間最短曲線的參數(shù)方程式中常數(shù)c1和c2由始末兩點位置確定練習：畫出經(jīng)過(0,0)和(1,1)的下落時間最短曲線。連接兩個點上凹的唯一一段旋輪線目前二十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@21/1045.1泛函與變分原理2泛函的極值問題的直接解法基本做法：瑞利--里茲(Rayleigh-Ritz)法(1)選定一組具有相對完備性的基函數(shù)，構(gòu)造一個線性組合的近似函數(shù)(2)將含有n個待定系數(shù)的構(gòu)造函數(shù)作為近似的極值函數(shù)，代入泛函目前二十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@22/1045.1泛函與變分原理(3)為了求泛函的極值，按照多元函數(shù)取極值的必要條件(4)求解以上方程組，求出就可以得到極值函數(shù)的近似解(5)再將含有n+1個待定系數(shù)的函數(shù)作為近似極值函數(shù)，重復(2)~(4)，就可以得到極值函數(shù)新的近似解。如果連續(xù)兩次所得到的結(jié)果接近，就認為最后得到的函數(shù)就是極值函數(shù)的近似解。目前二十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@23/1045.1泛函與變分原理例5.1.2求下列泛函的極值函數(shù)。解：為了滿足邊界條件，取基函數(shù)為近似函數(shù)為目前二十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@24/1045.1泛函與變分原理當n=1時代入泛函取極值目前二十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@25/1045.1泛函與變分原理計算得到近似函數(shù)同理n=2時利用歐拉方程，得到的精確解目前二十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@26/1045.1泛函與變分原理目前二十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@27/1045.1泛函與變分原理泛函的極值問題可以通過變分運算產(chǎn)生一個微分方程和相應(yīng)的邊界條件，即歐拉方程，其解對應(yīng)于最簡泛函的極值函數(shù)。也就是泛函的極值問題可以等價為在一定邊界條件下求解微分方程問題。

變分原理

通過求解一個相應(yīng)的泛函的極小函數(shù)而得到偏微分方程邊值問題的解。

有限元法正是里茲法與有限差分法相結(jié)合的成果，它取長補短地在理論上以變分為基礎(chǔ)，在具體方法構(gòu)造上又利用了有限差分法網(wǎng)格離散化處理的思想。目前二十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@28/1045.1泛函與變分原理

20世紀60年代初首次提出結(jié)構(gòu)力學計算有限元概念的克拉夫（Clough）教授形象地將其描繪為：“有限元法=RayleighRitz法＋分片函數(shù)”。有限元法是RayleighRitz法的一種局部化情況。不同于求解滿足整個定義域邊界條件的允許函數(shù)的RayleighRitz法（往往是困難的），有限元法將函數(shù)定義在簡單幾何形狀（如二維問題中的三角形或任意四邊形）的單元域上（分片函數(shù)），且不考慮整個定義域的復雜邊界條件，這是有限元法優(yōu)于其它近似方法的原因之一。目前二十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@29/1045.2基于變分原理的有限元法對于具有不同物理性質(zhì)和數(shù)學模型的問題，有限元法的基本做法是相同的，只是具體公式推導和運算求解不同。有限元法基本做法首先把待求的偏微分方程邊值問題轉(zhuǎn)化為等價的變分問題。然后通過有限單元剖分的離散處理，構(gòu)造一個分片解析的有限元子空間。通過構(gòu)造近似函數(shù)，把變分問題近似地轉(zhuǎn)化為有限元子空間中的多元函數(shù)極值問題，由此直接利用RayleighRitz法探求變分問題的近似解（極值函數(shù)解），以此作為所求邊值問題的近似解。

目前二十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@30/1045.2基于變分原理的有限元法有限元法具體求解步驟

建立積分方程根據(jù)變分原理或方程余量與權(quán)函數(shù)正交化原理，建立與微分方程初邊值問題等價的積分表達式，這是有限元法的出發(fā)點。區(qū)域單元剖分根據(jù)求解區(qū)域的形狀及實際問題的物理特點，將區(qū)域剖分為若干相互連接、不重疊的單元。區(qū)域單元劃分是采用有限元方法的前期準備工作，這部分工作量比較大，除了給計算單元和節(jié)點進行編號和確定相互之間的關(guān)系之外，還要表示節(jié)點的位置坐標，同時還需要列出自然邊界和本質(zhì)邊界的節(jié)點序號和相應(yīng)的邊界值。目前三十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@31/1045.2基于變分原理的有限元法確定單元基函數(shù)根據(jù)單元中節(jié)點數(shù)目及對近似解精度的要求，選擇滿足一定插值條件的插值函數(shù)作為單元基函數(shù)。有限元方法中的基函數(shù)是在單元中選取的，由于各單元具有規(guī)則的幾何形狀，在選取基函數(shù)時可遵循一定的法則。

單元分析將各個單元中的求解函數(shù)用單元基函數(shù)的線性組合表達式進行逼近；再將近似函數(shù)代入積分方程，并對單元區(qū)域進行積分，可獲得含有待定系數(shù)(即單元中各節(jié)點的參數(shù)值)的代數(shù)方程組，稱為單元有限元方程?？傮w合成在得出單元有限元方程之后，將區(qū)域中所有單元有限元方程按一定法則進行累加，形成總體有限元方程。目前三十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@32/1045.2基于變分原理的有限元法

邊界條件的處理一般邊界條件有三種形式，對于第二類邊界條件，一般在積分表達式中可自動得到滿足。對于第一類邊界條件和第三類邊界條件，需按一定法則對總體有限元方程進行修正滿足。解有限元方程根據(jù)邊界條件修正的總體有限元方程組，是含所有待定未知量的封閉方程組，采用適當?shù)臄?shù)值計算方法求解，可求得各節(jié)點的函數(shù)值。目前三十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@33/1045.2基于變分原理的有限元法

有限元分析可分成三個階段：

前置處理、計算求解和后置處理。前置處理是建立有限元模型，完成單元網(wǎng)格劃分；后置處理則是采集處理分析結(jié)果，使用戶能簡便提取信息，了解計算結(jié)果。目前三十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@34/1045.2基于變分原理的有限元法1.求解區(qū)域離散

離散單元基本要求：各單元只能在頂點處相交。不同單元在邊界處相連，既不能相互分離又不能相互重疊。各單元節(jié)點編號循序應(yīng)一致，一律按逆時針方向，從最小節(jié)點號開始。同一單元節(jié)點編號相差不能太懸殊，對多區(qū)域的編號，按區(qū)域連續(xù)編號。把求解區(qū)域分割成有限個單元體的集合。單元體形狀原則上是任意的，一般取有規(guī)則形體。有限元法計算步驟目前三十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@35/1045.2基于變分原理的有限元法三角單元是經(jīng)常使用的單元剖分方法，剖分時應(yīng)注意幾下幾點：

三角形不能重疊。不能把一個三角形的頂點取為相鄰三角形的邊上。剖分的三角形應(yīng)該避免鈍角。三角形不可過于狹長，最長邊一般不大于最短邊的3倍。三角形三邊之比盡量接近1。不能把一個三角形跨越不同的介質(zhì)。每個三角形最多只有一個邊在邊界上。三角形單元面積越小，計算精度越高目前三十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@36/1045.2基于變分原理的有限元法把求解區(qū)域劃分m個三角形有限單元，共有n個節(jié)點在有限單元e(j,k,l)上進行分片線性插值，插值函數(shù)為

2.選擇近似函數(shù)目前三十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@37/1045.2基于變分原理的有限元法在單元節(jié)點上求解以上方程組可以得到3.求解單元形函數(shù)目前三十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@38/1045.2基于變分原理的有限元法同理可以求出目前三十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@39/1045.2基于變分原理的有限元法目前三十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@40/1045.2基于變分原理的有限元法則插值函數(shù)可以寫為單元形函數(shù)（基函數(shù)）目前四十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@41/1045.2基于變分原理的有限元法三角元e插值函數(shù)可以改寫為矩陣形式目前四十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@42/1045.2基于變分原理的有限元法下面以泊松方程為例討論有限元解法所對應(yīng)的泛函為4.建立單元特征式難點：尋找與微分方程對應(yīng)的泛函目前四十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@43/1045.2基于變分原理的有限元法在第e個三角元的泛函由于目前四十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@44/1045.2基于變分原理的有限元法改寫為矩陣形式目前四十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@45/1045.2基于變分原理的有限元法其中同理目前四十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@46/1045.2基于變分原理的有限元法三角元e的泛函其中目前四十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@47/1045.2基于變分原理的有限元法改寫Ke到所有n個節(jié)點，即把擴充部分添零，以方便總體矩陣的處理其中目前四十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@48/1045.2基于變分原理的有限元法求解區(qū)域上的總體泛函其中變分問題被離散化的多元二次函數(shù)的極值問題5.建立系統(tǒng)有限元方程目前四十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@49/1045.2基于變分原理的有限元法根據(jù)多元函數(shù)極值理論得到第i點有限元方程即求解上述有限元方程（線性代數(shù)方程組），就可以得到節(jié)點上的函數(shù)值。

目前四十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@50/1045.2基于變分原理的有限元法獲得有限元方程之后，就可以選擇各種方法求解相應(yīng)的代數(shù)方程組，常用方法有高斯消去法、列元素消去法、迭代法等等。在變分問題中第二類、第三類邊界條件已經(jīng)自然包含在泛函達到極值的要求中，不必單獨處理，稱為自然滿足的邊界條件，只需考慮第一類強加邊界條件，強加邊界條件的處理方法因代數(shù)方程組的解法而異。6.有限元方程求解與邊界條件處理目前五十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@51/1045.2基于變分原理的有限元法迭代法求解：凡是遇到邊界節(jié)點所對應(yīng)的方程均不迭代，節(jié)點值始終保持給定值，不必單獨處理邊界。直接法求解：節(jié)點m為邊界，函數(shù)值um=u0，處理方法為，把對角元素的特征元素設(shè)置為1，即kmm=1,然后把m行與m列的其它元素全部設(shè)置為0，方程的等式右邊改為給定的函數(shù)值u0，其它元素則要減去該節(jié)點處理前對應(yīng)的m列的特征系數(shù)kim與u0的乘積。目前五十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@52/1045.2基于變分原理的有限元法

例5.2.1一個邊長為1的二維正方形靜電場域，電位函數(shù)為φ(x,y)，邊界條件如圖所示，試用有限元法確定二維靜電場域的電位分布。解：該二維靜電場域的電位函數(shù)φ(x,y),可以用下列第一類邊界條件的偏微分方程描述：目前五十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@53/1045.2基于變分原理的有限元法按照右圖進行三角形單元剖分，單元編號按照從左到右，從下到上的順序編號。節(jié)點編號：1(0,0)2(0,1)3(0.5,0.5)4(1,0)5(1,1)三角形單元編號：e(j,k,l)單元內(nèi)頂點按逆時針編號1(1,3,2)2(1,4,3)3(2,3,5)4(3,4,5)第一類邊界條件目前五十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@54/1045.2基于變分原理的有限元法對于三角元1(1,3,2)目前五十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@55/1045.2基于變分原理的有限元法擴展到全部節(jié)點目前五十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@56/1045.2基于變分原理的有限元法同樣，對于三角元2(1,4,3)目前五十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@57/1045.2基于變分原理的有限元法擴展到全部節(jié)點目前五十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@58/1045.2基于變分原理的有限元法同樣，對于三角元3(2,3,5)目前五十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@59/1045.2基于變分原理的有限元法擴展到全部節(jié)點目前五十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@60/1045.2基于變分原理的有限元法同樣，對于三角元4(3,4,5)目前六十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@61/1045.2基于變分原理的有限元法擴展到全部節(jié)點目前六十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@62/1045.2基于變分原理的有限元法全部節(jié)點K目前六十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@63/1045.2基于變分原理的有限元法有限元方程可以采用迭代法和直接解法，求解此線性代數(shù)方程組。目前六十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@64/1045.2基于變分原理的有限元法迭代法：迭代公式為代入初值目前六十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@65/1045.2基于變分原理的有限元法直接解法：需要處理第一類邊界條件目前六十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@66/1045.2基于變分原理的有限元法差分法得到差分遞推公式目前六十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@67/1045.2基于變分原理的有限元法差分網(wǎng)格利用邊界條件目前六十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@68/1045.3matlab有限元法工具箱大型通用有限元商業(yè)軟件國外軟件ANSYS、ADINA、ABAQUS、MSC等國內(nèi)軟件FEPG、JFEX、KMAS等

Matlab偏微分方程工具箱（PDEToolbox）提供了利用有限元法、圖形界面求解偏微分方程的計算環(huán)境。PDEtool有較大的局限性，可以求解特殊PDE問題，比如只能求解二階PDE問題，并且不能解決偏微分方程組，但是它提供了GUI界面，從繁雜的編程中解脫出來了，同時還可以通過File->SaveAs直接生成M代碼。目前六十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@69/1045.3matlab有限元法工具箱PDEToolbox求解偏微分方程類型

1橢圓型方程（Elliptic）2拋物線型方程（Parabolic）3雙曲型方程（Hyperbolic）目前六十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@70/1045.3matlab有限元法工具箱

4特征值方程（Eigenmodes）上述微分方程中c、a、d、f在橢圓型方程中可以為函數(shù)，但在其它方程中必須為常數(shù)。目前七十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@71/1045.3matlab有限元法工具箱PDEToolbox邊界條件

1狄里赫利條件（Didchlet）2諾依曼條件（Neumann）n為邊界上的單位外法線矢量，h、r、q、g可以為函數(shù)目前七十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@72/1045.3matlab有限元法工具箱PDEToolbox啟動

1啟動2界面目前七十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@73/1045.3matlab有限元法工具箱PDEToolbox菜單

Options打開或關(guān)閉柵格調(diào)整柵格大小打開或關(guān)閉捕捉柵格功能繪圖軸的坐標范圍打開或關(guān)閉繪圖方軸關(guān)閉幫助信息圖形縮放選擇應(yīng)用模式重新顯示圖形目前七十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@74/1045.3matlab有限元法工具箱PDEToolbox菜單

Draw進入繪圖模式對角點繪矩形固定中心繪矩形矩形對角點繪橢圓固定中心繪橢圓繪多邊形旋轉(zhuǎn)已選圖形將幾何描述矩陣輸出到主工作空間目前七十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@75/1045.3matlab有限元法工具箱PDEToolbox菜單

Boundary進入邊界模式對已選邊界輸入條件顯示邊界區(qū)域標識開關(guān)顯示子區(qū)域標識開關(guān)刪除已選的子域邊界刪除所有的子域邊界將分解幾何矩陣、邊界條件矩陣輸出到主工作空間目前七十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@76/1045.3matlab有限元法工具箱PDEToolbox菜單

PDE進入偏微分方程模式顯示子區(qū)域標識開關(guān)調(diào)整PDE參數(shù)和類型將PDE參數(shù)輸出到主工作空間目前七十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@77/1045.3matlab有限元法工具箱PDEToolbox菜單

Mesh輸入網(wǎng)格模式初始化三角形網(wǎng)格加密當前三角形網(wǎng)格優(yōu)化網(wǎng)格退回上一步用數(shù)字化的顏色顯示網(wǎng)格質(zhì)量，大于0.6可接受顯示網(wǎng)格節(jié)點標識顯示三角形網(wǎng)格標識修改網(wǎng)格生成參數(shù)輸出網(wǎng)格矩陣到主工作空間目前七十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@78/1045.3matlab有限元法工具箱PDEToolbox菜單

Solve對已經(jīng)定義的偏微分方程求解調(diào)整解PDE的參數(shù)輸出解到主工作空間Plot顯示圖形解繪圖參數(shù)設(shè)置輸出動畫目前七十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@79/1045.3matlab有限元法工具箱PDEToolbox求解步驟

求解區(qū)域設(shè)置應(yīng)用模式設(shè)置輸入邊界條件微分方程參數(shù)設(shè)定網(wǎng)格剖分初值和誤差設(shè)置解方程圖形解顯示參數(shù)設(shè)置File->SaveAs直接生成M代碼目前七十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@80/1045.3matlab有限元法工具箱例5.3.1如圖帶有矩形孔（0.1*0.8）的金屬板（1*1.6），金屬板左側(cè)保持在100℃，右側(cè)熱量可以向環(huán)境定常流動，上下側(cè)及內(nèi)孔保持絕熱，初始溫度為0℃。求t=0.1、0.3、0.5、1.5s時金屬板溫度分布

解：此問題可以表示為如下定解問題目前八十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@81/1045.3matlab有限元法工具箱

求解區(qū)域設(shè)置

提示符輸入pdetool1選擇畫矩形2畫矩形3雙擊矩形，彈出對話框，輸入準確矩形參數(shù)目前八十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@82/1045.3matlab有限元法工具箱

同樣畫出矩形孔，利用兩個矩形運算得到求解區(qū)域圖形運算目前八十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@83/1045.3matlab有限元法工具箱

2.應(yīng)用模式設(shè)置目前八十三頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@84/1045.3matlab有限元法工具箱

3輸入邊界條件1點擊，顯示邊界2雙擊邊界，彈出邊界條件窗口h=1,r=100目前八十四頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@85/1045.3matlab有限元法工具箱

輸入每個邊的邊界條件紅色：Dirichlet藍色：Neumann目前八十五頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@86/1045.3matlab有限元法工具箱

4.微分方程參數(shù)設(shè)定1點擊，設(shè)置方程，彈出窗口2拋物線型d=1c=1a=0f=0目前八十六頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@87/1045.3matlab有限元法工具箱

5.網(wǎng)格剖分點擊，網(wǎng)格剖分點擊，加密網(wǎng)格目前八十七頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@88/1045.3matlab有限元法工具箱

6.初值和誤差設(shè)置單擊Solve菜單中的Paramenters目前八十八頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@89/1045.3matlab有限元法工具箱

7.解方程點擊，解方程目前八十九頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@90/1045.3matlab有限元法工具箱

8.圖形解顯示參數(shù)設(shè)置目前九十頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@91/1045.3matlab有限元法工具箱

不同時刻溫度分布目前九十一頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@92/1045.3matlab有限元法工具箱

不同時刻溫度分布目前九十二頁\總數(shù)一百零五頁\編于七點HarbinInstituteofTechnologyYangkunkyang@93/1045.3matlab有限元法工具箱

MATLAB除了提供有限元工具箱求解二階偏微分方