BJ1111-Chapter6二維單元及二維問題分析

第六章二維單元及二維問題分析孫會(huì)本章內(nèi)容一、二維單元（重點(diǎn)、掌握）二、等參單元（重點(diǎn)、掌握）三、平面應(yīng)力問題（重點(diǎn)，熟悉）四、軸對(duì)稱問題（了解）五、ANSYS中的二維單元（重點(diǎn)，掌握）六、ANSYS應(yīng)用（重點(diǎn)，掌握）七、結(jié)果驗(yàn)證（重點(diǎn)，熟悉）一、二維單元1、矩形單元2、二次四邊形單元3、線性三角形單元4、二次三角形單元5、軸對(duì)稱單元1、矩形單元以直散熱片為例當(dāng)溫度沿X方向和Y方向均發(fā)生變化時(shí)，此時(shí)需采用二維單元描述散熱片的溫度分布情況。圖

使用矩形單元描述二維溫度分布的例子1、矩形單元（局部坐標(biāo)）對(duì)任意矩形單元：圖

典型的矩形單元

得到各個(gè)形函數(shù)1、矩形單元（局部坐標(biāo)）進(jìn)一步可推廣到任一變量：1、矩形單元（自然坐標(biāo)）自然坐標(biāo)：無(wú)量綱，且局部坐標(biāo)系x，y的原點(diǎn)與自然坐標(biāo)點(diǎn)=－1，＝－1一致。圖

自然坐標(biāo)下的四邊形單元

1、矩形單元二維矩形單元形函數(shù)的性質(zhì)：形函數(shù)在相應(yīng)節(jié)點(diǎn)處的值為1，而在其它節(jié)點(diǎn)處的值為0；形函數(shù)的和為1。2、二次四邊形單元8節(jié)點(diǎn)二次四邊形單元：4節(jié)點(diǎn)四邊形單元的高階單元，適合對(duì)曲線形邊界問題建模。與線性單元相比，對(duì)于同樣數(shù)目的單元，二次單元要提供的結(jié)果更精確。圖8節(jié)點(diǎn)的二次四邊形單元

2、二次四邊形單元8節(jié)點(diǎn)二次單元的一般形式：形函數(shù)的推導(dǎo)：利用節(jié)點(diǎn)值，求解方程組2、二次四邊形單元各角節(jié)點(diǎn)的形函數(shù)：中間節(jié)點(diǎn)的形函數(shù)：3、線性三角形單元雙線性矩形單元最主要的缺點(diǎn)：不能很好地滿足彎曲邊界的要求。相比之下，描述二維溫度分布的三角形單元就能較好地描述彎曲邊界。圖

使用三角形單元描述二維溫度分布

3、線性三角形單元三角形區(qū)域內(nèi)的獨(dú)立變量，如溫度的變化：圖

三角形單元

為求解未知系數(shù)，需利用節(jié)點(diǎn)值：3、線性三角形單元A是三角形單元的面積：3、線性三角形單元（整體坐標(biāo)）3、線性三角形單元三角形形函數(shù)的基本性質(zhì)：形函數(shù)在相應(yīng)節(jié)點(diǎn)處的值為1，而在其它節(jié)點(diǎn)處的值為0；形函數(shù)的和為1。3、線性三角形單元（自然/面積坐標(biāo)）設(shè)三角形區(qū)域內(nèi)有一坐標(biāo)為（X，Y）的點(diǎn)P，將點(diǎn)P和節(jié)點(diǎn)i、j、k相連，則這個(gè)三角形的面積就分為3個(gè)更小的面積A1、A2、A3。圖

三角形單元的自然（面積）坐標(biāo)

3、線性三角形單元（自然/面積坐標(biāo)）一般而言，三角形單元的自然(面積)坐標(biāo)、、定義：注意只有兩個(gè)自然坐標(biāo)是獨(dú)立的。因?yàn)槿切巫匀唬娣e）坐標(biāo)與形函數(shù)Si、Sj、Sk完全相同：4、二次三角形單元相對(duì)于線性三角形單元，二次三角形單元可以提供更為準(zhǔn)確的結(jié)果。由自然坐標(biāo)表示的二次三角形單元的形函數(shù)：圖

二次三角形單元

5、軸對(duì)稱單元實(shí)際工程應(yīng)用中存在一類特殊的三維問題，其幾何形狀和載荷都關(guān)于軸對(duì)稱，這類問題可應(yīng)用二維軸對(duì)稱單元來(lái)分析。圖

軸對(duì)稱單元模型

5、軸對(duì)稱單元（三角形單元）對(duì)于三角形線性單元，其上任一未知量：以r和z坐標(biāo)的形式重新表達(dá)上面的形函數(shù)，這種坐標(biāo)常用于軸對(duì)稱三角形單元中。5、軸對(duì)稱單元（三角形單元）典型軸對(duì)稱三角形單元及其坐標(biāo)如圖所示。用r和z替換X和Y得到如下形函數(shù)：圖

軸對(duì)稱三角形單元

5、軸對(duì)稱單元（矩形單元）矩形單元形函數(shù)：用r代替x，用z代替y，并利用如下關(guān)系：圖

軸對(duì)稱矩形單元

5、軸對(duì)稱單元（矩形單元）軸對(duì)稱矩形單元的形函數(shù)：二、等參單元一維問題中的等參公式和等參單元： 使用單一一組參數(shù)（如形函數(shù)）定義u，v，T等未知變量，并使用同樣的參數(shù)（同一形函數(shù)）表示幾何關(guān)系對(duì)于二維單元也存在相似情況。二、等參單元以固體力學(xué)問題為例：?jiǎn)卧獌?nèi)任意一點(diǎn)的位移:單元內(nèi)任意一點(diǎn)的位置:圖平面應(yīng)力問題中的四邊形單元

三、平面應(yīng)力問題1、基本概念與相關(guān)公式2、單元?jiǎng)偠染仃嚕ㄡ槍?duì)三角形等參單元）3、單元載荷矩陣（針對(duì)三角形等參單元）4、單元?jiǎng)偠染仃嚕ㄡ槍?duì)四邊形等參單元）1、基本概念與相關(guān)公式微元法材料體內(nèi)任一點(diǎn)圖

任一點(diǎn)上的應(yīng)力分量外力作用產(chǎn)生內(nèi)力產(chǎn)生應(yīng)力、應(yīng)變?nèi)绾蚊枋霾牧象w內(nèi)任一點(diǎn)的應(yīng)力、應(yīng)變狀態(tài)？1、基本概念與相關(guān)公式(1)任一點(diǎn)的應(yīng)力狀態(tài)

采用6個(gè)獨(dú)立的應(yīng)力分量表示，其中XX、YY、ZZ正應(yīng)力，XY、YZ、XZ剪應(yīng)力。

三維問題變?yōu)槠矫鎽?yīng)力問題：Z方向沒有施加力備注：在應(yīng)力分量的表示中，第一個(gè)下角標(biāo)對(duì)應(yīng)應(yīng)力所在面的法向方向，第二個(gè)下角標(biāo)對(duì)應(yīng)應(yīng)力的方向。1、基本概念與相關(guān)公式圖

平面應(yīng)力狀態(tài)1、基本概念與相關(guān)公式(2)任一點(diǎn)的應(yīng)變狀態(tài)

采用6個(gè)獨(dú)立的應(yīng)力分量表示，其中

XX、

YY、

ZZ正應(yīng)變，

XY、

YZ、

XZ剪應(yīng)變。

Z方向沒有位移w=0三維問題變?yōu)槠矫鎽?yīng)變問題1、基本概念與相關(guān)公式(3)應(yīng)力和應(yīng)變的關(guān)系---虎克定律

E彈性模量；泊松比；G彈性剪切模量。1、基本概念與相關(guān)公式對(duì)于平面應(yīng)力問題，虎克定律變?yōu)椋簯?yīng)變和位移的關(guān)系：2、單元?jiǎng)偠染仃嚥捎米钚】倓?shì)能法： 對(duì)于由n個(gè)單元m個(gè)節(jié)點(diǎn)構(gòu)成的穩(wěn)定系統(tǒng)，平衡位置上產(chǎn)生的位移總會(huì)使系統(tǒng)的總勢(shì)能最?。簩?duì)應(yīng)單元?jiǎng)偠染仃噷?duì)應(yīng)單元載荷矩陣2、單元?jiǎng)偠染仃噷?duì)于平面應(yīng)力問題，應(yīng)變能：關(guān)鍵：構(gòu)造應(yīng)變矩陣的表達(dá)形式2、單元?jiǎng)偠染仃嚕ㄈ切螁卧┮跃€性三角形單元為例：等參公式2、單元?jiǎng)偠染仃嚕ㄈ切螁卧┣箨P(guān)于節(jié)點(diǎn)位移的微分：

得到單元?jiǎng)偠染仃嘖(e)：V是單元的體積3、單元載荷矩陣（三角形單元）對(duì)于最小總勢(shì)能法： 對(duì)于由n個(gè)單元m個(gè)節(jié)點(diǎn)構(gòu)成的穩(wěn)定系統(tǒng)，平衡位置上產(chǎn)生的位移總會(huì)使系統(tǒng)的總勢(shì)能最?。簩?duì)應(yīng)單元?jiǎng)偠染仃噷?duì)應(yīng)單元載荷矩陣關(guān)鍵：獲得功的表達(dá)形式3、單元載荷矩陣（三角形單元）集中載荷Q：

3、單元載荷矩陣（三角形單元）分量為px和py的分布載荷所做的功：如使用三角形單元 表示位移，則分布載荷所做的功：u和v：x和y方向的位移；A：分布載荷作用范圍的面積，其大小為單元厚度t和分布載荷作用邊長(zhǎng)的乘積。3、單元載荷矩陣----推導(dǎo)過(guò)程推導(dǎo)過(guò)程：3、單元載荷矩陣----推導(dǎo)過(guò)程3、單元載荷矩陣----推導(dǎo)過(guò)程3、單元載荷矩陣沿k－i邊：沿i－j邊：沿j－k邊：4、四邊形等參單元上面針對(duì)平面應(yīng)力問題，采用最小總勢(shì)能法，以線性三角形單元為例，推導(dǎo)獲得了單元?jiǎng)偠染仃嚭蛦卧d荷矩陣。實(shí)際上，對(duì)于四邊形等參單元，其方法類似，只不過(guò)過(guò)程更加繁瑣，這里給出最終結(jié)果。4、單元?jiǎng)偠染仃嚕ㄋ倪呅螁卧┧倪呅蔚葏卧膯卧獎(jiǎng)偠染仃嚕簍e

單元厚度；四、軸對(duì)稱問題軸對(duì)稱問題：幾何形狀和載荷關(guān)于某個(gè)軸對(duì)稱，可用二維軸對(duì)稱單元進(jìn)行分析。軸對(duì)稱單元?jiǎng)偠染仃嚨耐茖?dǎo)步驟與平面應(yīng)力問題類似。差別在于，平面應(yīng)力問題采用直角坐標(biāo)系；而軸對(duì)稱問題使用柱坐標(biāo)系。軸對(duì)稱三角形單元的剛度矩陣：七、ANSYS中的二維單元ANSYS提供了許多二維單元，這些單元大多數(shù)是基于線性、二次四邊形和三角形形函數(shù)的。下面舉例說(shuō)明一些二維結(jié)構(gòu)力學(xué)單元。五、ANSYS中的二維單元PLANE2：二維6節(jié)點(diǎn)三角形單元每個(gè)節(jié)點(diǎn)有2個(gè)自由度，即在節(jié)點(diǎn)的x和y方向都能平移。圖ANSYS中的二維單元

五、ANSYS中的二維單元PLANE42：二維4節(jié)點(diǎn)四邊形單元，常用于固體力學(xué)問題建模。每個(gè)節(jié)點(diǎn)有2個(gè)自由度，即每個(gè)節(jié)點(diǎn)都能在x和y方向平移。圖ANSYS中的二維單元

五、ANSYS中的二維單元PLANE82：二維8節(jié)點(diǎn)四邊形單元，常用于二維結(jié)構(gòu)問題的建模。PLANE42的高階形式，可用于對(duì)含曲線邊界的問題建模，且計(jì)算精度更高。每個(gè)節(jié)點(diǎn)有2個(gè)自由度，即能在x和y方向平移。圖ANSYS中的二維單元

五、ANSYS中的二維單元使用高階單元：可獲得更好的結(jié)果和更高的精度；計(jì)算時(shí)間通常會(huì)更長(zhǎng)，這是因?yàn)樗婕暗降膯卧仃嚁?shù)值積分也更多。六、ANSYS應(yīng)用設(shè)有一支撐書架的鋼托架，E＝29106lb/in2，＝0.3。該托架上表面承受均布載荷，左端固定。試在給定的載荷和約束下，繪制托架變形后的形狀，并確定托架的主應(yīng)力和vonMises應(yīng)力。六、ANSYS應(yīng)用圖

鋼托架示意圖六、ANSYS應(yīng)用具體過(guò)程軟件演示。七、結(jié)果驗(yàn)證基本原理基于靜力平衡條件。具體操作可取