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