有限元法基礎(chǔ)-4單元和插值函數(shù)的構(gòu)造PPT課件

1 第四章單元與插值函數(shù) 4 1面積坐標4 2Lagrange單元4 3Serendipity單元4 4體積坐標4 5Hermite插值 2 4 單元與插值函數(shù) 通過變分法或加權(quán)余量法建立有限元方程時 首先是在確定單元形狀后 在單元域內(nèi)假設(shè)場函數(shù)的試解 本章重點介紹構(gòu)造單元插值函數(shù)規(guī)范化形式的兩類自然坐標的建立方法和特點構(gòu)造單元插值函數(shù)的兩類方法的步驟和特點 有限元法基礎(chǔ) 3 4 單元與插值函數(shù) 關(guān)鍵概念自然坐標面積坐標體積坐標Lagrange單元Serendipity單元 有限元法基礎(chǔ) 4 4 單元與插值函數(shù) 廣義坐標有限元法的存在的問題 1 建立單元插值函數(shù)方法繁瑣2 形成單元矩陣過于復雜 有限元法基礎(chǔ) 5 4 單元與插值函數(shù) 單元插值函數(shù)的構(gòu)造與求解問題的微分方程無關(guān)插值函數(shù)的構(gòu)造方法與單元形狀有關(guān)與單元節(jié)點數(shù)量與位置有關(guān)與單元節(jié)點DOF的類型和數(shù)量有關(guān) 有限元法基礎(chǔ) 6 4 單元與插值函數(shù) 有限元法基礎(chǔ) 7 4 1面積坐標 定義在三角形內(nèi)任意一點P的位置由其三角形子域的面積與三角形面積的比值確定 即其中A為三角形面積 為的面積 為的面積 為的面積 有限元法基礎(chǔ) 8 4 1面積坐標 記則三角形內(nèi)的點P表示為稱為面積坐標 有限元法基礎(chǔ) 9 面積坐標的性質(zhì)1 與j m平行的線上具有相同的Li 4 1面積坐標 有限元法基礎(chǔ) 10 4 1面積坐標 2 角點坐標為i 1 0 0 j 01 0 m 0 0 1 3 形心坐標為4 三角形三條邊的坐標為j m邊 Li 0 m i邊 Lj 0 i j邊 Lm 05 三個坐標只有2個是獨立的 有限元法基礎(chǔ) 11 4 1面積坐標 面積坐標與直角坐標的關(guān)系三角形單元的面積三角形內(nèi)任意點P x y 有限元法基礎(chǔ) 12 4 1面積坐標 有限元法基礎(chǔ) 13 4 1面積坐標 面積坐標的微積分運算1 導數(shù) 有限元法基礎(chǔ) 14 4 1面積坐標 2 面積分3 i j邊長為l的線積分 有限元法基礎(chǔ) 15 4 1面積坐標 例 有限元法基礎(chǔ) 16 4 1面積坐標 例 均質(zhì)等厚單元的自重 有限元法基礎(chǔ) 17 4 1面積坐標 用面積坐標給出的單元的插值函數(shù)以面積坐標作為三角形單元的自然坐標 表示的插值函數(shù) 對每一個節(jié)點來講 插值函數(shù)是對稱的 有限元法基礎(chǔ) 18 4 1面積坐標 1 線性單元 3節(jié)點三角形單元根據(jù)形函數(shù)的特點這樣可用過其他兩節(jié)點的直線方程來構(gòu)成 例如節(jié)點1 可用2 3邊的直線方程來構(gòu)成插值函數(shù) 即 有限元法基礎(chǔ) 19 4 1面積坐標 2 二次單元 6節(jié)點三角形單元節(jié)點1 節(jié)點4 通用表達式 角節(jié)點中節(jié)點注 有限元法基礎(chǔ) 20 4 1面積坐標 2 三次單元 10節(jié)點三角形單元節(jié)點1節(jié)點4節(jié)點10 有限元法基礎(chǔ) 21 4 2Lagrange單元 單元場函數(shù)的插值表示為插值函數(shù)滿足下列性質(zhì) 有限元法基礎(chǔ) 22 4 2Lagrange單元 一維Lagrange插值1 總體坐標下的位移插值函數(shù)對于n個節(jié)點的一維單元 節(jié)點坐標為多項式插值可達n 1階 即 有限元法基礎(chǔ) 23 4 2Lagrange單元 當 2時令 則引進無量綱坐標 有限元法基礎(chǔ) 24 4 2Lagrange單元 2 自然坐標下的位移插值函數(shù)對于n個節(jié)點的一維單元 節(jié)點坐標為多項式插值可達n 1階 即或 有限元法基礎(chǔ) 25 4 2Lagrange單元 當n 2時 通式 有限元法基礎(chǔ) 26 4 2Lagrange單元 當n 3時 有限元法基礎(chǔ) 27 4 2Lagrange單元 二維Lagrange單元二維Lagrange單元的場插值函數(shù)由一維Lagrange插值分別在兩個方向插值 即場插值函數(shù)為 有限元法基礎(chǔ) 28 4 2Lagrange單元 可以證明 有限元法基礎(chǔ) 29 4 2Lagrange單元 一次單元 4節(jié)點單元雙線性插值 有限元法基礎(chǔ) 30 4 2Lagrange單元 二次單元 9節(jié)點單元角節(jié)點邊中節(jié)點內(nèi)部節(jié)點 有限元法基礎(chǔ) 31 4 2Lagrange單元 三維Lagrange單元單元節(jié)點場插值函數(shù) 有限元法基礎(chǔ) 32 4 2Lagrange單元 可以證明 有限元法基礎(chǔ) 33 4 2Lagrange單元 Lagrange單元族的特點1 插值函數(shù)構(gòu)造方便2 高次單元內(nèi)部節(jié)點過多 影響計算效率 有限元法基礎(chǔ) 平面單元內(nèi)部節(jié)點數(shù) n 1 m 1 34 4 3Serendipity單元 Serendipity單元族單元的節(jié)點僅配置在角點和邊界上 不改變精度的情況下 減少內(nèi)節(jié)點 Irons等首先提出 按字面意思是意外發(fā)現(xiàn)的 但有規(guī)律可循 有限元法基礎(chǔ) 35 4 3Serendipity單元 Serendipity單元族 有限元法基礎(chǔ) 36 4 3Serendipity單元 Serendipity插值函數(shù)的構(gòu)造4節(jié)點單元的插值函數(shù)與Lagrange單元相同 有限元法基礎(chǔ) 37 4 3Serendipity單元 8節(jié)點單元在邊中點的插值函數(shù) 有限元法基礎(chǔ) 38 4 3Serendipity單元 顯然 有限元法基礎(chǔ) 39 4 3Serendipity單元 Serendipity插值函數(shù)的一般構(gòu)造方法 有限元法基礎(chǔ) 40 4 3Serendipity單元 Serendipity單元族 有限元法基礎(chǔ) 41 4 3Serendipity單元 二次平面單元 8節(jié)點單元 有限元法基礎(chǔ) 42 4 3Serendipity單元 劃線法構(gòu)造插值函數(shù)二次平面單元 8節(jié)點單元由 得 有限元法基礎(chǔ) 43 4 3Serendipity單元 平面四邊形8節(jié)點單元插值函數(shù) 有限元法基礎(chǔ) 44 4 3Serendipity單元 三維Serendipity單元族1 線性單元 8節(jié)點單元 有限元法基礎(chǔ) 45 4 3Serendipity單元 2 二次單元 20節(jié)點單元插值函數(shù)完備到二次 有限元法基礎(chǔ) 46 4 3Serendipity單元 Serendipity插值與Lagrange插值的差異Serendipity插值函數(shù)的多項式表示與Lagrange插值相比少都是二次完備 沒達到三次完備 有限元法基礎(chǔ) 47 4 3Serendipity單元 有限元法基礎(chǔ) 48 4 3Serendipity單元 Lagrange插值項 有限元法基礎(chǔ) 49 4 3Serendipity單元 Serendipity插值項 有限元法基礎(chǔ) 50 4 3Serendipity單元 Serendipity插值與Lagrange插值的差異 有限元法基礎(chǔ) 51 4 3Serendipity單元 在邊 二次變化在邊 二次變化單元的每條邊上有三個節(jié)點 因此插值是協(xié)調(diào)的 有限元法基礎(chǔ) 52 4 4體積坐標 定義 四面體中任一點P的位置由下列參數(shù)確定 有限元法基礎(chǔ) 53 4 4體積坐標 四面體單元族 有限元法基礎(chǔ) 54 4 4體積坐標 1 線性單元 4節(jié)點單元2 二次單元 10節(jié)點單元角節(jié)點邊中節(jié)點 有限元法基礎(chǔ) 55 4 4體積坐標 3 三次單元 20節(jié)點單元角節(jié)點邊中節(jié)點面內(nèi)節(jié)點 有限元法基礎(chǔ) 56 4 4體積坐標 體積坐標的微分復合函數(shù)求導公式 有限元法基礎(chǔ) 57 4 4體積坐標 體積坐標的積分 有限元法基礎(chǔ) 58 4 5Hermite插值 特點1 多項式插值2 節(jié)點參數(shù)包含導數(shù) 例如3 在插值點上導數(shù)也連續(xù)4 0階連續(xù)的Hermite插值就是Lagrange插值 有限元法基礎(chǔ) 59 4 5Hermite插值 例 一維問題有n個節(jié)點 包含有2n個節(jié)點未知參數(shù) 可唯一確定2n 1次多項式的插值函數(shù) 有限元法基礎(chǔ) 60 4 5Hermite插值 容易得到其中為Lagrange插值函數(shù) 有限元法基礎(chǔ) 61 4 5Hermite插值 采用局部無量綱坐標時 端點為 有限元法基礎(chǔ) 62 4 5Hermite插值 在無量綱坐標下 插值函數(shù)曲線 有限元法基礎(chǔ) 63 4 5Hermite插值 二階Hermite插值函數(shù) 有限元法基礎(chǔ) 64 4 5Hermite插值 例 2節(jié)點梁單元 有限元法基礎(chǔ) 65 4 6五面體單元 有限元法基礎(chǔ) 線性五面體單元單元內(nèi)任意一點的坐標用面積坐