有限元法基礎(chǔ)平板彎曲問題詳解演示文稿

有限元法基礎(chǔ)平板彎曲問題詳解演示文稿現(xiàn)在是1頁\一共有47頁\編輯于星期一（優(yōu)選）有限元法基礎(chǔ)平板彎曲問題現(xiàn)在是2頁\一共有47頁\編輯于星期一10.平板彎曲問題關(guān)鍵概念C1類板單元C0類板單元非協(xié)調(diào)板單元協(xié)調(diào)板單元Ks奇異性條件Ke非奇異性條件DKT板單元有限元法基礎(chǔ)3現(xiàn)在是3頁\一共有47頁\編輯于星期一10.平板彎曲問題有限元法基礎(chǔ)4Z

XY中面板的特點：在一個方向的尺度遠遠小于其他兩個方向，中面是平面，只承受橫向載荷?，F(xiàn)在是4頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)5一.基本方程Kirchhoff假設(shè)1）變形前垂直于中面的直線段，變形后依然垂

直于中面，并且忽略它的伸縮變形2）忽略厚度方向的應(yīng)力，即現(xiàn)在是5頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)6板中任意點的位移表示為三維問題二維問題現(xiàn)在是6頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)7定義廣義應(yīng)變和廣義內(nèi)力廣義應(yīng)力應(yīng)變關(guān)系抗彎剛度現(xiàn)在是7頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)8應(yīng)力與廣義內(nèi)力的關(guān)系平衡方程以中面撓度w表示的微分方程現(xiàn)在是8頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)9邊界條件1）固支類邊界2）簡支類邊界3）給定力邊界現(xiàn)在是9頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)10現(xiàn)在是10頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)11最小勢能原理以上廣義應(yīng)變是撓度w的二階導(dǎo)數(shù)關(guān)系，基于此理論的板單元是C1類連續(xù)問題?，F(xiàn)在是11頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)12有限元列式

設(shè)插值函數(shù)為

通過泛函取駐值得有限元方程

單元剛度矩陣現(xiàn)在是12頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)13二.非協(xié)調(diào)矩形板單元

每節(jié)點有3DOF，4節(jié)點單元共12個節(jié)點DOF。

現(xiàn)在是13頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)14插值函數(shù)

按廣義坐標有限元法，在Pascal三角形中選取12項多項式

現(xiàn)在是14頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)15

現(xiàn)在是15頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)16以節(jié)點DOF表示插值函數(shù)

表示為矩陣形式

現(xiàn)在是16頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)17以自然坐標表示

現(xiàn)在是17頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)18收斂性檢查1）

位移模式代表剛體位移

沿Z向的平移和繞y軸和X軸的轉(zhuǎn)動2）位移模式代表常曲率

滿足完備性要求現(xiàn)在是18頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)193）單元間連續(xù)性檢查

單元邊界為x=常數(shù)或y=常數(shù)，w是三次變化曲線。以2－3邊為例，可以由4個參數(shù)完全確定。在2－3邊的法向?qū)?shù)為

為三次x變化,而在邊界上只有2個參數(shù)。

法向?qū)?shù)不連續(xù)現(xiàn)在是19頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)204）由于在單元間邊界上法向?qū)?shù)不連續(xù)，所以插值函數(shù)是非協(xié)調(diào)的；5）單元不滿足收斂準則，但是可以驗證該單元通過補片試驗（PatchTest），故當單元剖分不斷縮小時，計算結(jié)果還是能收斂于精確解。

通過補片試驗實際驗算現(xiàn)在是20頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)21例:均布載荷下四邊固支方形薄板,利用對稱性取四分之一板計算現(xiàn)在是21頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)22例:載荷作用下方形薄板,利用對稱性取四分之一板計算注：由于是非協(xié)調(diào)元，位移解并補滿足下界條件現(xiàn)在是22頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)23三.3節(jié)點三角形非協(xié)調(diào)板單元共有3×3＝9個DOF三次完備多項式

ijm10項現(xiàn)在是23頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)24插值函數(shù)

面積坐標剛體位移常應(yīng)變現(xiàn)在是24頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)25坐標變換代入節(jié)點坐標求出系數(shù)，得到形函數(shù)

現(xiàn)在是25頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)26位移插值函數(shù)的特點

插值函數(shù)包含有完備的線性項和二次項，能正確反映剛體位移和常應(yīng)變；在單元邊界上，w是三次變化，可由兩端節(jié)點的w和w,s唯一確定，w是協(xié)調(diào)的；在單元邊界上，w,n是二次變化的,不能由兩端節(jié)點的w,n確定，w,n是非協(xié)調(diào)的。

現(xiàn)在是26頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)27

Irons等已證明如果單元網(wǎng)格是由3組等間距直線產(chǎn)生的，單元能夠通過補片試驗，并收斂于解析解?，F(xiàn)在是27頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)283節(jié)點三角板元四.協(xié)調(diào)單元思路：在邊界(如i-j)上尋找校正函數(shù)，具有性質(zhì)1）在全部邊界上2）在j-m,i-m邊上3)在i-j上，按二次變化，且在中點上取1

單元邊界上w,n二次變化非協(xié)調(diào)元現(xiàn)在是28頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)29插直函數(shù)w是非協(xié)調(diào)元的產(chǎn)值函數(shù),為待定常數(shù)。目的：調(diào)整使在單元邊界中點處的w,n等于兩端節(jié)點的w,n的平均值，也即使得邊界上法向?qū)?shù)線性化，可由兩端點的值唯一確定。

現(xiàn)在是29頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)30

的確定線性化要求，在邊界中點處原插值函數(shù)計算出的各邊界中點值原插值函數(shù)計算的邊界中點平均值現(xiàn)在是30頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)31校正函數(shù)可以驗證以上函數(shù)滿足校正函數(shù)的要求，即在全部邊界上等于零，在i-m和j-m邊法向?qū)?shù)為零，在i-j邊上二次變化。令現(xiàn)在是31頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)32單元特點單元協(xié)調(diào)性完全滿足隨著單元尺寸不斷減小，解能單調(diào)收斂于精確解有高階校正函數(shù)，要提高數(shù)值積分階次實際計算時，單元往往過于剛硬現(xiàn)在是32頁\一共有47頁\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)33例：簡支方板受中心集中力現(xiàn)在是33頁\一共有47頁\編輯于星期一協(xié)調(diào)薄板元列式的其他方法1）組合單元法

將四個三角形單元組合為一個四邊形單元，選用特殊插值函數(shù)，使之滿足連續(xù)性要求，并凝聚內(nèi)部節(jié)點2）多節(jié)點參數(shù)法

引入高階導(dǎo)數(shù)項作為節(jié)點DOF，以提高邊界的協(xié)調(diào)性，例如10.1Kirchhoff板單元有限元法基礎(chǔ)34現(xiàn)在是34頁\一共有47頁\編輯于星期一Reissner-Mindlin變形假設(shè)

變形前垂直于中面的直線段，變形后仍然保持為直線段，但不在垂直于中面。10.2Mindlin板單元有限元法基礎(chǔ)35現(xiàn)在是35頁\一共有47頁\編輯于星期一廣義應(yīng)變變分原理10.2Mindlin板單元有限元法基礎(chǔ)36一般取k=5/6現(xiàn)在是36頁\一共有47頁\編輯于星期一位移插值10.2Mindlin板單元有限元法基礎(chǔ)37現(xiàn)在是37頁\一共有47頁\編輯于星期一應(yīng)變-節(jié)點DOF矩陣10.2Mindlin板單元有限元法基礎(chǔ)38現(xiàn)在是38頁\一共有47頁\編輯于星期一有限元方程由泛函取極值條件得單元剛度矩陣10.2Mindlin板單元有限元法基礎(chǔ)39現(xiàn)在是39頁\一共有47頁\編輯于星期一邊界條件

三種類型：1）2）3）

給定位移屬于強制邊界條件，給定內(nèi)力屬于自然邊界條件10.2Mindlin板單元有限元法基礎(chǔ)40現(xiàn)在是40頁\一共有47頁\編輯于星期一剪切自鎖

與Timoshenko梁單元一樣Mindlin板元中剪切能量引入后，存在罰因子現(xiàn)象解決辦法有減縮積分、假設(shè)應(yīng)變等方法多變量有限元也是常見的處理方法10.2Mindlin板單元有限元法基礎(chǔ)41現(xiàn)在是41頁\一共有47頁\編輯于星期一積分方案目標：保證K非奇異性和Ks奇異性保證K非奇異性的必要條件M單元數(shù)；ng高斯積分點數(shù)；d應(yīng)變分量數(shù)；N系統(tǒng)的獨立DOF數(shù)。N＝節(jié)點總數(shù)×每節(jié)點DOF數(shù)－給定約束數(shù)10.2Mindlin板單元有限元法基礎(chǔ)42現(xiàn)在是42頁\一共有47頁\編輯于星期一對Mindlin板單元，保證K非奇異性的必要條件nb

和ns分別為Kb和Ks的高斯積分點數(shù)；db和ds分別為Kb和Ks的應(yīng)變分量數(shù)，db＝3，ds＝2。保證Ks奇異性的必要條件10.2Mindlin板單元有限元法基礎(chǔ)43現(xiàn)在是43頁\一共有47頁\編輯于星期一積分方案10.2Mindlin板單元有限元法