第八講-有限元法

等效積分形式：與原有微分方程和定解條件完全等價(jià)。加權(quán)余量法：對(duì)場(chǎng)函數(shù)進(jìn)行近似，令加權(quán)余量等于零。伽遼金法：加權(quán)函數(shù)與場(chǎng)函數(shù)的試探函數(shù)（基函數(shù)、形函數(shù)）相同。小結(jié)：伽遼金法是有限元法中使用最為普遍的。第八講_有限元法(8)（5）伽遼金法簡(jiǎn)單地說(shuō)，將近似解的試探函數(shù)作為權(quán)函數(shù)。更簡(jiǎn)潔的形式：伽遼金法的一般表達(dá)式引入變分等效積分形式第八講_有限元法(8)靜態(tài)線彈性有限元定解問(wèn)題真實(shí)位移的變分，連續(xù)可導(dǎo)。在給定位移的邊界上，虛應(yīng)變高斯定律張量形式矩陣形式8/9/20133第八講_有限元法(8)等效積分形式：與原有微分方程和定解條件完全等價(jià)。加權(quán)余量法：對(duì)場(chǎng)函數(shù)進(jìn)行近似，令加權(quán)余量等于零。伽遼金法：加權(quán)函數(shù)與場(chǎng)函數(shù)的試探函數(shù)（基函數(shù)、形函數(shù)）相同。小結(jié)：伽遼金法是有限元法中使用最為普遍的。第八講_有限元法(8)

基本概念偏微分方程和偏微分方程組：一個(gè)未知函數(shù)及其偏導(dǎo)數(shù)組成的方程叫偏微分方程，兩個(gè)以上未知函數(shù)及其偏導(dǎo)數(shù)組成的方程組叫偏微分方程組。方程組中未知函數(shù)和方程個(gè)數(shù)相等，叫封閉的偏微分方程組（或完全的）。偏微分方程的階和偏微分方程組的階：方程中偏導(dǎo)數(shù)的最高階次叫偏微分方程的階;偏微分方程組的階是方程組中各偏微分方程的階數(shù)之和。第八講_有限元法(8)線性、非線性和擬線性偏微分方程:a)方程中所有出現(xiàn)未知函數(shù)或其偏導(dǎo)數(shù)的項(xiàng)都是未知函數(shù)的一次式的方程叫線性方程b)未知函數(shù)項(xiàng)或未知函數(shù)偏導(dǎo)數(shù)項(xiàng)不是一次式的方程叫非線性方程；c)非線性方程中所有未知函數(shù)的最高階偏導(dǎo)數(shù)是一次式的方程叫擬線性方程。第八講_有限元法(8)齊次和非齊次偏微分方程偏微分方程分類(lèi)：Auxx+2Buxy+Cuyy+Dux+Euy+f=0Δ（=B2–AC）>0：雙曲型，波動(dòng)方程.Δ（=B2–AC）=0：拋物型，熱傳導(dǎo)方程.Δ（=B2–AC）<0：橢圓型，

位勢(shì)方程.第八講_有限元法(8)定解問(wèn)題：偏微分方程+定解條件（邊界條件+初始條件）a)初值問(wèn)題n階方程有n個(gè)初始條件，初始條件偏導(dǎo)數(shù)的最高階次是n–1.b)邊值問(wèn)題i:第一類(lèi)邊界條件（Dirichlet）,ii:第二類(lèi)邊界條件（Neuman）,iii:第三類(lèi)邊界條件（Robin）.第八講_有限元法(8)

偏微分方程的解a)解：使偏微分方程兩端恒等的有定義的函數(shù)叫偏微分方程的解.b)通解:對(duì)于n階方程，未知函數(shù)有m個(gè)自變量，其通解由n個(gè)獨(dú)立的滿足一定可微要求的函數(shù)組成，且每個(gè)函數(shù)有m–1個(gè)自變量.

uxyz=0,u=f1(x,y)+f2(x,z)+f3(y,z)

utt–c2uxx=0,u=f1(x+ct)+f2(x–ct).

c)定解問(wèn)題的解：滿足邊界條件和初始條件的通解.第八講_有限元法(8)

線性疊加原理解的存在、唯一性和穩(wěn)定性解的性質(zhì)a）橢圓形方程的極值只能在邊界達(dá)到。解在內(nèi)部沒(méi)有弱間斷，解在邊界上間斷，在內(nèi)部也是充分光滑的，邊界條件是封閉的.b）雙曲型方程沒(méi)有像橢圓形方程那樣的極值原理，解在內(nèi)部可以有弱間斷，邊界條件不是封閉的。第八講_有限元法(8)

泛函與變分（Functionalandvariation）1．泛函函數(shù)的函數(shù)a)兩端固定的曲線長(zhǎng)度：b)彈性桿的總勢(shì)能：c)溫度場(chǎng)泛函：式中f,u,T叫做泛函的容許函數(shù)：滿足一定邊界條件和連續(xù)性的所有函數(shù)有限元法的基本原理曲線長(zhǎng)度總勢(shì)能溫度場(chǎng)泛函第八講_有限元法(8)

變分定義

a)容許函數(shù)的變分有限元法的基本原理第八講_有限元法(8)第八講_有限元法(8)i）

泛函的值由1個(gè)自變量的函數(shù)確定ii）泛函的值由有3個(gè)自變量的函數(shù)確定iii）泛函的值由有3個(gè)自變量的2個(gè)函數(shù)確定第八講_有限元法(8)d）變分運(yùn)算第八講_有限元法(8)3．變分問(wèn)題a)函數(shù)的極值問(wèn)題（無(wú)約束和約束）b)變分問(wèn)題：求泛函的極值函數(shù)c)泛函極值函數(shù)的必要條件d)歐拉方程第八講_有限元法(8)有限元法的基本原理第八講_有限元法(8)e)定解問(wèn)題與變分問(wèn)題此式即桿的平衡方程，它就是變分的歐拉方程。有限元法的基本原理i）

固定邊界變分問(wèn)題與基本邊界條件兩端約束的彈性桿問(wèn)題：第八講_有限元法(8)ii）非固定邊界變分問(wèn)題與自然邊界條件邊界條件：泛函：此式即桿的平衡方程第八講_有限元法(8)一端約束（指定位移）的彈性桿解法1：Lagrange乘子法構(gòu)造新泛函iii）含有約束條件的變分問(wèn)題Lagrange乘子法與原定解問(wèn)題完全等效，代價(jià)就是引入了新的求解變量第八講_有限元法(8)解法2：罰函數(shù)法構(gòu)造新泛函有限元法的基本原理當(dāng)k無(wú)窮大時(shí)，則滿足第一類(lèi)邊界條件。不引入新的求解變量。第八講_有限元法(8)小結(jié)：1.定解問(wèn)題（微分方程加定解條件）等價(jià)于相應(yīng)泛函取極值。2.泛函取極值就是有限元方法的理論基礎(chǔ)，將微分形式變成了積分形式。3.不是所有的定解問(wèn)題都存在相應(yīng)的泛函。4.不存在泛函的定解問(wèn)題，可以直接用更廣義的加權(quán)余量法。第八講_有限元法(8)泛函的變分彈性長(zhǎng)桿的定解問(wèn)題微分方程定解條件對(duì)應(yīng)泛函第八講_有限元法(8)近似函數(shù)表示的微分方程的殘差邊界條件的殘差其思想是使由近似函數(shù)表示的微分方程的殘差和邊界條件的殘差的加權(quán)積分為零

2．加權(quán)余量法直接從微分方程出發(fā)的一種積分方法。假設(shè)未知函數(shù)采用近似函數(shù)表達(dá)：有限元法的基本原理第八講_有限元法(8)伽遼金加權(quán)余量法弱形式與變分結(jié)果一致。分部積分伽遼金加權(quán)余量法

3．加權(quán)余量法弱形式：泛函取極值變分等于零可以看出，第八講_有限元法(8)取加權(quán)函數(shù)的試探函數(shù)與近似解的試探函數(shù)相同

加遼金有限元法解題過(guò)程1）構(gòu)造加權(quán)殘差積分方程2）離散化3）單元分析4）總體分析5）建立和求解有限元方程組4．加遼金加權(quán)余量法第八講_有限元法(8)弱形式1）構(gòu)造加權(quán)余量積分方程第八講_有限元法(8)b）單元的總體節(jié)點(diǎn)編號(hào)和局部節(jié)點(diǎn)編號(hào)，單元e=I，總體節(jié)點(diǎn)編號(hào)：1，2，局部節(jié)點(diǎn)編號(hào)：1，2；單元e=III，總體節(jié)點(diǎn)編號(hào)：3，4，局部節(jié)點(diǎn)編號(hào)：1，2。2）離散化第八講_有限元法(8)a)在自然坐標(biāo)系中構(gòu)造單元近似解：b）構(gòu)造加權(quán)函數(shù)c）單元坐標(biāo)：3）單元分析加權(quán)函數(shù)、近似解試探函數(shù)、坐標(biāo)插值函數(shù)的類(lèi)型一致第八講_有限元法(8)d）單元平衡方程第八講_有限元法(8)a）建立選擇矩陣：4）總體分析第八講_有限元法(8)b)組集單元?jiǎng)偠染仃嘽)組集等效節(jié)點(diǎn)載荷d）解以節(jié)點(diǎn)為未知量的方程組第八講_有限元法(8)熱傳導(dǎo)問(wèn)題的有限元方法第八講_有限元法(8)1）傅里葉定律熱傳導(dǎo)方程1.一維問(wèn)題熱流密度與溫度梯度成正比。q:單位時(shí)間、單位面積流過(guò)的熱量單位：W/(m·K)第八講_有限元法(8)2）平衡方程Q=cpmΔT比熱容：cp1千克的物質(zhì)的溫度上升（或下降）1攝氏度所需的能量。單位：W·s/(kg·K)擴(kuò)散率的單位：m^2/s第八講_有限元法(8)x=13m,t=1y,theta=0.1x=13m,t=6667y,theta=0.98第八講_有限元法(8)本構(gòu)關(guān)系第八講_有限元法(8)第八講_有限元法