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

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

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

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

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

通過(guò)泛函取駐值得有限元方程

單元?jiǎng)偠染仃嚞F(xiàn)在是12頁(yè)\一共有47頁(yè)\編輯于星期一10.1Kirchhoff板單元有限元法基礎(chǔ)13二.非協(xié)調(diào)矩形板單元

每節(jié)點(diǎn)有3DOF，4節(jié)點(diǎn)單元共12個(gè)節(jié)點(diǎn)DOF。

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

按廣義坐標(biāo)有限元法，在Pascal三角形中選取12項(xiàng)多項(xiàng)式

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

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

表示為矩陣形式

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

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

位移模式代表剛體位移

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

三種類(lèi)型：1）2）3）

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

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

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