ANSYS有限元基礎(chǔ)教程課件-第2章

第2章平面問(wèn)題的有限元法有限元分析的3個(gè)步驟：

離散化單元分析整體分析第2章平面問(wèn)題的有限元法有限元分析的3個(gè)步驟：離散化12.1結(jié)構(gòu)的離散化實(shí)例：將一個(gè)受力的連續(xù)體離散化離散化:把連續(xù)的結(jié)構(gòu)看成由有限個(gè)單元組成的集合體。目的：建立有限元計(jì)算模型以三角形單元為例2.1結(jié)構(gòu)的離散化實(shí)例：將一個(gè)受力的連續(xù)體2注意事項(xiàng)對(duì)稱性的利用

如果結(jié)構(gòu)與載荷都有對(duì)稱性可資利用，可取其中的一半或1/4等作為分析對(duì)象，能減少很多工作量。

圖為平面薄板的離散化模型注意事項(xiàng)對(duì)稱性的利用圖為平面薄板的離散化模型32）節(jié)點(diǎn)的布置：①集中載荷的作用點(diǎn)，②分布載荷強(qiáng)度的突變點(diǎn)，③分布載荷與自由邊界的分界點(diǎn)，④支承點(diǎn)，⑤厚度不同或材料不同的區(qū)域等都應(yīng)取為節(jié)點(diǎn)。3）對(duì)于重要的或應(yīng)力變化急劇的部位，單元應(yīng)劃得小些，對(duì)于次要的和應(yīng)力變化緩慢的部位，單元可劃得大些，“中間地帶”以大小逐漸變化的單元來(lái)過(guò)渡。2.節(jié)點(diǎn)的選擇和單元的劃分1）單元形狀和尺寸可自由調(diào)整。2）節(jié)點(diǎn)的布置：①集中載荷的作用點(diǎn)，②分布載荷強(qiáng)度的突變點(diǎn)43.節(jié)點(diǎn)的編號(hào)在節(jié)點(diǎn)編號(hào)時(shí)，應(yīng)注意盡量使同一單元的相鄰節(jié)點(diǎn)的號(hào)碼差值盡可能地小些，以便縮小剛度矩陣的帶寬，節(jié)約計(jì)算機(jī)存儲(chǔ)。

（a）（b）如圖，（a）與（b）單元?jiǎng)澐窒嗤?，（b）的編號(hào)要比（a）的編號(hào)為好，即節(jié)點(diǎn)應(yīng)順短邊編號(hào)為好。3.節(jié)點(diǎn)的編號(hào)（a）52.2單元分析單元分析的主要任務(wù)是推導(dǎo)單元節(jié)點(diǎn)位移與單元節(jié)點(diǎn)力之間的轉(zhuǎn)換關(guān)系，實(shí)質(zhì)上就是求出單元?jiǎng)偠染仃?。單元分析的步驟（實(shí)施過(guò)程）：2.2單元分析單元分析的主要任務(wù)是推導(dǎo)6三角形單元節(jié)點(diǎn)位移向量：

三角形單元節(jié)點(diǎn)力向量：

從離散化的結(jié)構(gòu)中任取一個(gè)三角形單元e—單元首先對(duì)節(jié)點(diǎn)編碼：稱為局部碼1.位移函數(shù)的概念三角形單元節(jié)點(diǎn)位移向量：三角形單元節(jié)點(diǎn)力向量：從離散化的7廣泛使用多項(xiàng)式來(lái)構(gòu)造位移函數(shù)。將單元中的位移分布假定是坐標(biāo)的簡(jiǎn)單函數(shù)，稱為位移函數(shù)。設(shè)單元內(nèi)任意一點(diǎn)的位移含有6個(gè)待定參數(shù)，稱為廣義坐標(biāo)。

（2-6）廣泛使用多項(xiàng)式來(lái)構(gòu)造位移函數(shù)。將單元中的位移8求形函數(shù)三角形單元的面積：求形函數(shù)三角形單元的面積：9令式中順序輪換簡(jiǎn)寫為：

矩陣形式：Ni，Nj，Nm為形函數(shù)，是關(guān)于坐標(biāo)x、y的線性函數(shù)[N]為形函數(shù)矩陣

常數(shù)令式中順序輪換簡(jiǎn)寫為：矩陣形式：Ni，Nj，Nm為形函數(shù)102.位移函數(shù)收斂準(zhǔn)則在有限元法中，把能夠滿足條件（1）、（2）的單元，稱為完備單元；滿足條件（3）的單元，稱為協(xié)調(diào)單元。

（3）位移函數(shù)應(yīng)盡可能反映位移的連續(xù)性。要求所選擇的位移函數(shù)既能使單元內(nèi)部的位移保持連續(xù)，又能使相鄰單元之間的位移保持連續(xù)，后者是指單元之間不出現(xiàn)互相脫離和互相嵌入的現(xiàn)象。（2）位移函數(shù)必須能反映單元的常量應(yīng)變。在位移函數(shù)中的一次項(xiàng)就是提供單元中的常量應(yīng)變的。（1）位移函數(shù)必須能反映單元的剛體位移。常數(shù)項(xiàng)就是用于提供剛體位移的。2.位移函數(shù)收斂準(zhǔn)則在有限元法中，把能夠滿足條件（1）、（11圖2-4相鄰三角形單元的位移協(xié)調(diào)性分析圖2-4相鄰三角形單元的位移協(xié)調(diào)性分析121）要考慮到解的收斂性，即要考慮到完備性和協(xié)調(diào)性的要求。巴斯卡三角形3.選擇單元位移函數(shù)的一般原則3）多項(xiàng)式中的項(xiàng)數(shù)必須等于或稍大于單元邊界上的外節(jié)點(diǎn)的自由度數(shù)。通常是取項(xiàng)數(shù)與單元的外節(jié)點(diǎn)的自由度數(shù)相等。2）模式應(yīng)該與局部坐標(biāo)系的方位無(wú)關(guān)，這一性質(zhì)稱為幾何各向同性。1）要考慮到解的收斂性，即要考慮到完備性和協(xié)調(diào)性的要求。巴斯132.2.2單元應(yīng)變（2-10）（2-11）幾何方程：

2.2.2單元應(yīng)變（2-10）幾何方程：14[B]稱作幾何矩陣，是常數(shù)矩陣。因此，三角形單元是常應(yīng)變單元。[B]稱作幾何矩陣，是常數(shù)矩陣。因此，三角形單元是常應(yīng)變單元152.2.3單元應(yīng)力物理方程：?jiǎn)卧獞?yīng)力：令：[S]稱為單元應(yīng)力矩陣

由于三角形單元中的[D]，[B]矩陣都是常數(shù)矩陣，所以[S]矩陣也是常數(shù)矩陣。也就是說(shuō)，三角形單元內(nèi)的應(yīng)力分量也是常量。

2.2.3單元應(yīng)力物理方程：?jiǎn)卧獞?yīng)力：令：[S]稱162.2.4單元?jiǎng)偠染仃囂摴Ψ匠蹋喝切螁卧?.2.4單元?jiǎng)偠染仃囂摴Ψ匠蹋喝切螁卧?7令：?jiǎn)卧胶夥匠蹋▌偠确匠蹋﹩卧獎(jiǎng)偠染仃嚸總€(gè)分塊矩陣均為2×2階方陣。三角形單元的剛度矩陣為6×6階方陣。令：?jiǎn)卧胶夥匠蹋▌偠确匠蹋﹩卧獎(jiǎng)偠染仃嚸總€(gè)分塊矩陣均為2×18單元?jiǎng)偠染仃囉腥缦滦再|(zhì)：①每一個(gè)元素物理意義：是單位節(jié)點(diǎn)位移分量所引起的節(jié)點(diǎn)力分量。

②是對(duì)稱矩陣。③每一行（或列）元素之和為零。是奇異矩陣，④的元素決定于單元的形狀、大小、方位和彈性常數(shù)，而與單元的位置無(wú)關(guān)，即不隨單元（或坐標(biāo)軸）的平行移動(dòng)或作（n為整數(shù)）角度的轉(zhuǎn)動(dòng)而改變。

單元?jiǎng)偠染仃嚢垂?jié)點(diǎn)寫成分塊形式：?jiǎn)卧獎(jiǎng)偠染仃囉腥缦滦再|(zhì)：?jiǎn)卧獎(jiǎng)偠染仃嚢垂?jié)點(diǎn)寫成19

20例2-1如圖2-6所示平面應(yīng)力情形的直角三角形單元，直角邊長(zhǎng)均為，厚度為t，彈性模量為E，泊松比為，求單元?jiǎng)偠染仃?。圖2-6直角三角形單元解：（1）求例2-1如圖2-6所示平面應(yīng)力情形的直角三角形單元21（2）求（3）求（4）求（2）求（3）求（4）求222.3整體分析結(jié)構(gòu)的整體分析就是將離散后的所有單元通過(guò)節(jié)點(diǎn)連接成原結(jié)構(gòu)物進(jìn)行分析。——總體剛度矩陣?！麄€(gè)結(jié)構(gòu)上節(jié)點(diǎn)位移列陣——整個(gè)結(jié)構(gòu)上節(jié)點(diǎn)力列陣——整體平衡方程（整體剛度方程）分析過(guò)程：是將所有單元平衡方程組集在一起，形成總體平衡方程，引進(jìn)邊界條件后，求解整體節(jié)點(diǎn)位移向量。2.3整體分析結(jié)構(gòu)的整體分析就是將離散23平面問(wèn)題的整體剛度矩陣是由相關(guān)單元的單元?jiǎng)偠染仃囍械姆謮K矩陣集合而成，按節(jié)點(diǎn)編號(hào)對(duì)號(hào)入座，即剛度集成法。它的集成規(guī)律有下列幾點(diǎn)：1）先對(duì)每個(gè)單元求出其單元?jiǎng)偠染仃?，以分塊形式按節(jié)點(diǎn)編號(hào)順序排列。2）將單元?jiǎng)偠染仃嚁U(kuò)大階數(shù)為2n×2n，并將單元?jiǎng)偠染仃囍械姆謮K矩陣按局部碼與總碼的對(duì)應(yīng)關(guān)系，搬到擴(kuò)大后的矩陣中，形成單元貢獻(xiàn)矩陣。3）將所有單元貢獻(xiàn)矩陣同一位置上的分塊矩陣簡(jiǎn)單疊加成總體剛度矩陣中的一個(gè)子矩陣，各行各列都按以上步驟即形成總體剛度矩陣。總體剛度矩陣為2n×2n階，亦即n×n階分塊矩陣，n為節(jié)點(diǎn)總數(shù)。平面問(wèn)題的整體剛度矩陣是由相關(guān)單元的單元?jiǎng)偠染?4例2-2用剛度集成法求圖2-7所示結(jié)構(gòu)的整體剛度矩陣1）找出各單元局部碼與總碼的對(duì)應(yīng)關(guān)系解：2）分別寫出各個(gè)單元的分塊矩陣3）形成各單元貢獻(xiàn)陣?yán)?-2用剛度集成法求圖2-7所示結(jié)構(gòu)的整體剛度矩陣1）254）將擴(kuò)大后的4個(gè)單元貢獻(xiàn)陣相疊加，就得到總體剛度矩陣

整體剛度矩陣有以下一些性質(zhì)：1）整體剛度矩陣是對(duì)稱矩陣。2）整體剛度矩陣的主對(duì)角線上的元素總是正的。3）整體剛度矩陣是一個(gè)稀疏陣。4）整體剛度矩陣是一個(gè)奇異陣。4）將擴(kuò)大后的4個(gè)單元貢獻(xiàn)陣相疊加，就得到總體剛度矩陣整體26例2-3已知如圖a）所示的懸臂深梁，在右端面作用著均布拉力，其合力為P。采用如圖b）所示簡(jiǎn)單網(wǎng)格，設(shè)，厚度為t。試求節(jié)點(diǎn)位移。a）b）解：１）求出各單元?jiǎng)偠染仃嚴(yán)?-3已知如圖a）所示的懸臂深梁，在右端面作用著均布27２）形成各單元貢獻(xiàn)矩陣３）形成總體剛度矩陣

２）形成各單元貢獻(xiàn)矩陣３）形成總體剛度矩陣28ANSYS有限元基礎(chǔ)教程課件-第2章29４）形成整體載荷列陣

５）形成整體節(jié)點(diǎn)位移列陣

６）形成整體平衡方程

４）形成整體載荷列陣５）形成整體節(jié)點(diǎn)位移列陣６）形成整30７）引入邊界條件，求節(jié)點(diǎn)位移“化1置0法”處理：若已知節(jié)點(diǎn)在方向位移為零，則令

中的元素為1，而第行和列的其余元素都為零，中的第個(gè)元素變?yōu)榱?。若已知?jié)點(diǎn)在方向位移為零，則令

中的元素為1，而第行和列的其余元素都為零，中的第個(gè)元素變?yōu)榱?。７）引入邊界條件，求節(jié)點(diǎn)位移“化1置0法”處理：若已知節(jié)點(diǎn)312.4

有限元法解題過(guò)程與算例

有限元法的具體解題過(guò)程為：1）將結(jié)構(gòu)進(jìn)行離散化，包括單元?jiǎng)澐?、?jié)點(diǎn)編號(hào)、單元編號(hào)、節(jié)點(diǎn)坐標(biāo)計(jì)算、位移約束條件的確定。2）等效節(jié)點(diǎn)力的計(jì)算。3）剛度矩陣的計(jì)算。4）建立整體平衡方程，引入約束條件，求解節(jié)點(diǎn)位移。5）應(yīng)力計(jì)算。2.4有限元法解題過(guò)程與算例有限元法的具體解題過(guò)程32例2-4如圖2-20a所示兩端固支的矩形深梁，跨度為2a，梁高為a，厚度為t，已知E，，承受均布?jí)毫，試用有限元法求解此平面應(yīng)力問(wèn)題。a）b）圖2-20矩形深梁例2-4如圖2-20a所示兩端固支的矩形深梁，跨度為2a33解：利用對(duì)稱性，可取梁的一半分析，例如右半。1.劃分單元并準(zhǔn)備原始數(shù)據(jù)2.計(jì)算單元?jiǎng)偠染仃?31324解：利用對(duì)稱性，可取梁的一半分析，例如右半。1.劃分單元并343.集成整體剛度矩陣依照各單元局部編號(hào)與整體編號(hào)的對(duì)應(yīng)關(guān)系，兩個(gè)單元的貢獻(xiàn)矩陣分別為再集成整體剛度矩陣3.集成整體剛度矩陣依照各單元局部編號(hào)與整體354.處理載荷，形成整體平衡方程整體節(jié)點(diǎn)載荷列陣為組成結(jié)構(gòu)整體平衡方程4.處理載荷，形成整體平衡方程整體節(jié)點(diǎn)載荷列陣為組成結(jié)構(gòu)365.引入位移邊界條件，求解節(jié)點(diǎn)位移由于5.引入位移邊界條件，求解節(jié)點(diǎn)位移由于376.應(yīng)力計(jì)算在整體分析中求得節(jié)點(diǎn)位移之后，為了計(jì)算結(jié)構(gòu)上任意一點(diǎn)的應(yīng)變或應(yīng)力，應(yīng)該又返回到單元分析中去。計(jì)算單元①的應(yīng)力矩陣由整體節(jié)點(diǎn)位移向量獲取單元節(jié)點(diǎn)位移向量6.應(yīng)力計(jì)算計(jì)算單元①的應(yīng)力矩陣由整體節(jié)點(diǎn)位移向量獲取單元38計(jì)算應(yīng)力計(jì)算應(yīng)力392.5單元等效節(jié)點(diǎn)力一、單元自重三角形單元的厚度為t，重度為，面積為，自重沿y軸負(fù)方向。受自重載荷情形的等效節(jié)點(diǎn)力為單元重量的1/3。

圖2-14三角形單元

介紹幾種常用載荷作用下的等效節(jié)點(diǎn)力2.5單元等效節(jié)點(diǎn)力一、單元自重三角形單元40二、均布面力三角形單元的厚度為t，ij邊的長(zhǎng)度為，集度為圖2-15受均布面力三角形單元相當(dāng)于把作用于ij邊上的表面力按靜力等效平均分配到該邊兩端的節(jié)點(diǎn)上。二、均布面力三角形單元的厚度為t，ij邊的41三、線性分布面力三角形單元的厚度為t，ij邊的長(zhǎng)度為，表面力在點(diǎn)集度為圖2-16受線性分布面力三角形單元

相當(dāng)于將總載荷的2/3分配給點(diǎn)，1/3分配給點(diǎn)。

三、線性分布面力三角形單元的厚度為t，ij422.6邊界條件的處理常用的且比較方便的做法是以某種方法引入已知的節(jié)點(diǎn)位移（包括零位移約束），而保持方程原有的數(shù)目不變，只是修

和中某些元素，以避免計(jì)算機(jī)存儲(chǔ)作大的變動(dòng)。2.6邊界條件的處理常用的且比較方便的做43

1、“化1置0法”----為節(jié)點(diǎn)總碼編號(hào)1、“化1置0法”----為節(jié)點(diǎn)總碼編號(hào)44設(shè)已知節(jié)點(diǎn)位移為，當(dāng)引進(jìn)上述已知節(jié)點(diǎn)位移后，方程變成例如：為說(shuō)明這一過(guò)程，現(xiàn)考察一個(gè)只有四個(gè)方程的簡(jiǎn)單例子設(shè)已知節(jié)點(diǎn)位移為，當(dāng)452、“乘大數(shù)法”設(shè)已知節(jié)點(diǎn)位移為，當(dāng)引進(jìn)上述已知節(jié)點(diǎn)位移后，方程變成2、“乘大數(shù)法”設(shè)已知節(jié)點(diǎn)位移為462.7計(jì)算結(jié)果的整理繞節(jié)點(diǎn)平均法：

ABCDEF把環(huán)繞該節(jié)點(diǎn)的各單元應(yīng)力加以平均，視為該節(jié)點(diǎn)的應(yīng)力。2.7計(jì)算結(jié)果的整理繞節(jié)點(diǎn)平均法：AB47兩單元平均法：把相鄰兩單元應(yīng)力的平均值作為公共邊中點(diǎn)的應(yīng)力。

ABCDEFGHI78910兩單元平均法：把相鄰兩單元應(yīng)力的平均值作為公共邊中點(diǎn)的應(yīng)力。48采用上述兩種應(yīng)力平均法時(shí)應(yīng)注意幾點(diǎn)：（1）只有當(dāng)相連單元具有相同厚度和材料時(shí)平均法才有意義。（3）位于結(jié)構(gòu)邊界點(diǎn)的應(yīng)力不應(yīng)該用平均法求得，若用繞節(jié)點(diǎn)平均法則因其相連單元太少而不能得到較佳的近似值。這種情況往往改用內(nèi)部應(yīng)力點(diǎn)外推的辦法，去求它的近似值。（2）為了使繞節(jié)點(diǎn)平均法得來(lái)的應(yīng)力能夠較好地表示節(jié)點(diǎn)處的實(shí)際應(yīng)力，環(huán)繞該節(jié)點(diǎn)的各個(gè)單元的面積不應(yīng)相差太大。采用上述兩種應(yīng)力平均法時(shí)應(yīng)注意幾點(diǎn)：（3）位于結(jié)構(gòu)邊界點(diǎn)的應(yīng)492.8矩形單元

設(shè)有矩形單元，其邊長(zhǎng)分別為和，矩形的兩邊分別與軸平行。取矩形的四個(gè)角點(diǎn)作為節(jié)點(diǎn)。單元節(jié)點(diǎn)位移向量為：

2.8矩形單元設(shè)有矩形單元50在單元分析中為了計(jì)算上的方便和簡(jiǎn)化，我們引用一個(gè)無(wú)量綱的局部坐標(biāo)系，局部坐標(biāo)系的原點(diǎn)取在矩形的形心上，軸分別與整體坐標(biāo)軸平行，它們之間的坐標(biāo)變換為在局部坐標(biāo)系中，四個(gè)節(jié)點(diǎn)坐標(biāo)分別是即為（-1，-1），（1，-1），（1，1），（-1，1）。在單元分析中為了計(jì)算上的方便和簡(jiǎn)化，我們引用一個(gè)無(wú)量綱的局部511、位移模式

形函數(shù)雙線性模式從中求出合并寫成:（2-74）1、位移模式形函數(shù)雙線性模式從中求出合并寫成:（2-7452形函數(shù)寫成矩陣形式形函數(shù)矩陣（2-74）合并寫成:形函數(shù)寫成矩陣形式形函數(shù)矩陣（2-74）合并寫成:532、單元應(yīng)變寫成矩陣形式四節(jié)點(diǎn)矩形單元不再是常應(yīng)變單元。幾何方程將位移函數(shù)帶入上式復(fù)合函數(shù)求導(dǎo)法則2、單元應(yīng)變寫成矩陣形式四節(jié)點(diǎn)矩形單元不再是常應(yīng)變單元。幾何543、單元應(yīng)力四節(jié)點(diǎn)矩形單元不再是常應(yīng)力單元。根據(jù)物理方程：3、單元應(yīng)力四節(jié)點(diǎn)矩形單元不再是常應(yīng)力單元。根據(jù)物理方程：554、單元?jiǎng)偠染仃噷懗煞謮K形式：

4、單元?jiǎng)偠染仃噷懗煞謮K形式：565、單元等效節(jié)點(diǎn)力（1）對(duì)于單元的自重W，載荷列陣為

即移置于每一節(jié)點(diǎn)的載荷都是四分之一的自重。（2）如果單元在一個(gè)邊界上受有三角形分布的表面力，在該邊界上一個(gè)節(jié)點(diǎn)處為零，而在另一個(gè)節(jié)點(diǎn)處為最大，則將總表面力的三分之一移置到前一個(gè)節(jié)點(diǎn)，三分之二移置到后一個(gè)節(jié)點(diǎn)。5、單元等效節(jié)點(diǎn)力（1）對(duì)于單元的自重W，載荷列陣為（2）576、整體平衡方程將各單元的、和都擴(kuò)大到整個(gè)彈性體自由度的維數(shù)，再進(jìn)行疊加，便可得到整個(gè)彈性體的平衡方程，它仍具有如下的形式引入位移約束條件，解上述線性方程組可得節(jié)點(diǎn)位移，進(jìn)而可求各單元應(yīng)力。6、整體平衡方程將各單元的、587、矩形單元與三角形單元的比較：

3、但矩形單元也存在明顯的缺點(diǎn)：從單元的幾何形狀看，矩形單元比三角形單元的適應(yīng)性要差。

2、在彈性體中,若用相同數(shù)目的節(jié)點(diǎn)時(shí)，矩形單元比三角形單元能更好地反映應(yīng)力急劇變化的情況，所以計(jì)算精度高。1、矩形單元為雙線性位移模式，所以單元的應(yīng)力、應(yīng)變分量都不是常量。7、矩形單元與三角形單元的比較：3、但矩形59第2章平面問(wèn)題的有限元法有限元分析的3個(gè)步驟：

離散化單元分析整體分析第2章平面問(wèn)題的有限元法有限元分析的3個(gè)步驟：離散化602.1結(jié)構(gòu)的離散化實(shí)例：將一個(gè)受力的連續(xù)體離散化離散化:把連續(xù)的結(jié)構(gòu)看成由有限個(gè)單元組成的集合體。目的：建立有限元計(jì)算模型以三角形單元為例2.1結(jié)構(gòu)的離散化實(shí)例：將一個(gè)受力的連續(xù)體61注意事項(xiàng)對(duì)稱性的利用

如果結(jié)構(gòu)與載荷都有對(duì)稱性可資利用，可取其中的一半或1/4等作為分析對(duì)象，能減少很多工作量。

圖為平面薄板的離散化模型注意事項(xiàng)對(duì)稱性的利用圖為平面薄板的離散化模型622）節(jié)點(diǎn)的布置：①集中載荷的作用點(diǎn)，②分布載荷強(qiáng)度的突變點(diǎn)，③分布載荷與自由邊界的分界點(diǎn)，④支承點(diǎn)，⑤厚度不同或材料不同的區(qū)域等都應(yīng)取為節(jié)點(diǎn)。3）對(duì)于重要的或應(yīng)力變化急劇的部位，單元應(yīng)劃得小些，對(duì)于次要的和應(yīng)力變化緩慢的部位，單元可劃得大些，“中間地帶”以大小逐漸變化的單元來(lái)過(guò)渡。2.節(jié)點(diǎn)的選擇和單元的劃分1）單元形狀和尺寸可自由調(diào)整。2）節(jié)點(diǎn)的布置：①集中載荷的作用點(diǎn)，②分布載荷強(qiáng)度的突變點(diǎn)633.節(jié)點(diǎn)的編號(hào)在節(jié)點(diǎn)編號(hào)時(shí)，應(yīng)注意盡量使同一單元的相鄰節(jié)點(diǎn)的號(hào)碼差值盡可能地小些，以便縮小剛度矩陣的帶寬，節(jié)約計(jì)算機(jī)存儲(chǔ)。

（a）（b）如圖，（a）與（b）單元?jiǎng)澐窒嗤?，（b）的編號(hào)要比（a）的編號(hào)為好，即節(jié)點(diǎn)應(yīng)順短邊編號(hào)為好。3.節(jié)點(diǎn)的編號(hào)（a）642.2單元分析單元分析的主要任務(wù)是推導(dǎo)單元節(jié)點(diǎn)位移與單元節(jié)點(diǎn)力之間的轉(zhuǎn)換關(guān)系，實(shí)質(zhì)上就是求出單元?jiǎng)偠染仃?。單元分析的步驟（實(shí)施過(guò)程）：2.2單元分析單元分析的主要任務(wù)是推導(dǎo)65三角形單元節(jié)點(diǎn)位移向量：

三角形單元節(jié)點(diǎn)力向量：

從離散化的結(jié)構(gòu)中任取一個(gè)三角形單元e—單元首先對(duì)節(jié)點(diǎn)編碼：稱為局部碼1.位移函數(shù)的概念三角形單元節(jié)點(diǎn)位移向量：三角形單元節(jié)點(diǎn)力向量：從離散化的66廣泛使用多項(xiàng)式來(lái)構(gòu)造位移函數(shù)。將單元中的位移分布假定是坐標(biāo)的簡(jiǎn)單函數(shù)，稱為位移函數(shù)。設(shè)單元內(nèi)任意一點(diǎn)的位移含有6個(gè)待定參數(shù)，稱為廣義坐標(biāo)。

（2-6）廣泛使用多項(xiàng)式來(lái)構(gòu)造位移函數(shù)。將單元中的位移67求形函數(shù)三角形單元的面積：求形函數(shù)三角形單元的面積：68令式中順序輪換簡(jiǎn)寫為：

矩陣形式：Ni，Nj，Nm為形函數(shù)，是關(guān)于坐標(biāo)x、y的線性函數(shù)[N]為形函數(shù)矩陣

常數(shù)令式中順序輪換簡(jiǎn)寫為：矩陣形式：Ni，Nj，Nm為形函數(shù)692.位移函數(shù)收斂準(zhǔn)則在有限元法中，把能夠滿足條件（1）、（2）的單元，稱為完備單元；滿足條件（3）的單元，稱為協(xié)調(diào)單元。

（3）位移函數(shù)應(yīng)盡可能反映位移的連續(xù)性。要求所選擇的位移函數(shù)既能使單元內(nèi)部的位移保持連續(xù)，又能使相鄰單元之間的位移保持連續(xù)，后者是指單元之間不出現(xiàn)互相脫離和互相嵌入的現(xiàn)象。（2）位移函數(shù)必須能反映單元的常量應(yīng)變。在位移函數(shù)中的一次項(xiàng)就是提供單元中的常量應(yīng)變的。（1）位移函數(shù)必須能反映單元的剛體位移。常數(shù)項(xiàng)就是用于提供剛體位移的。2.位移函數(shù)收斂準(zhǔn)則在有限元法中，把能夠滿足條件（1）、（70圖2-4相鄰三角形單元的位移協(xié)調(diào)性分析圖2-4相鄰三角形單元的位移協(xié)調(diào)性分析711）要考慮到解的收斂性，即要考慮到完備性和協(xié)調(diào)性的要求。巴斯卡三角形3.選擇單元位移函數(shù)的一般原則3）多項(xiàng)式中的項(xiàng)數(shù)必須等于或稍大于單元邊界上的外節(jié)點(diǎn)的自由度數(shù)。通常是取項(xiàng)數(shù)與單元的外節(jié)點(diǎn)的自由度數(shù)相等。2）模式應(yīng)該與局部坐標(biāo)系的方位無(wú)關(guān)，這一性質(zhì)稱為幾何各向同性。1）要考慮到解的收斂性，即要考慮到完備性和協(xié)調(diào)性的要求。巴斯722.2.2單元應(yīng)變（2-10）（2-11）幾何方程：

2.2.2單元應(yīng)變（2-10）幾何方程：73[B]稱作幾何矩陣，是常數(shù)矩陣。因此，三角形單元是常應(yīng)變單元。[B]稱作幾何矩陣，是常數(shù)矩陣。因此，三角形單元是常應(yīng)變單元742.2.3單元應(yīng)力物理方程：?jiǎn)卧獞?yīng)力：令：[S]稱為單元應(yīng)力矩陣

由于三角形單元中的[D]，[B]矩陣都是常數(shù)矩陣，所以[S]矩陣也是常數(shù)矩陣。也就是說(shuō)，三角形單元內(nèi)的應(yīng)力分量也是常量。

2.2.3單元應(yīng)力物理方程：?jiǎn)卧獞?yīng)力：令：[S]稱752.2.4單元?jiǎng)偠染仃囂摴Ψ匠蹋喝切螁卧?.2.4單元?jiǎng)偠染仃囂摴Ψ匠蹋喝切螁卧?6令：?jiǎn)卧胶夥匠蹋▌偠确匠蹋﹩卧獎(jiǎng)偠染仃嚸總€(gè)分塊矩陣均為2×2階方陣。三角形單元的剛度矩陣為6×6階方陣。令：?jiǎn)卧胶夥匠蹋▌偠确匠蹋﹩卧獎(jiǎng)偠染仃嚸總€(gè)分塊矩陣均為2×77單元?jiǎng)偠染仃囉腥缦滦再|(zhì)：①每一個(gè)元素物理意義：是單位節(jié)點(diǎn)位移分量所引起的節(jié)點(diǎn)力分量。

②是對(duì)稱矩陣。③每一行（或列）元素之和為零。是奇異矩陣，④的元素決定于單元的形狀、大小、方位和彈性常數(shù)，而與單元的位置無(wú)關(guān)，即不隨單元（或坐標(biāo)軸）的平行移動(dòng)或作（n為整數(shù)）角度的轉(zhuǎn)動(dòng)而改變。

單元?jiǎng)偠染仃嚢垂?jié)點(diǎn)寫成分塊形式：?jiǎn)卧獎(jiǎng)偠染仃囉腥缦滦再|(zhì)：?jiǎn)卧獎(jiǎng)偠染仃嚢垂?jié)點(diǎn)寫成78

79例2-1如圖2-6所示平面應(yīng)力情形的直角三角形單元，直角邊長(zhǎng)均為，厚度為t，彈性模量為E，泊松比為，求單元?jiǎng)偠染仃?。圖2-6直角三角形單元解：（1）求例2-1如圖2-6所示平面應(yīng)力情形的直角三角形單元80（2）求（3）求（4）求（2）求（3）求（4）求812.3整體分析結(jié)構(gòu)的整體分析就是將離散后的所有單元通過(guò)節(jié)點(diǎn)連接成原結(jié)構(gòu)物進(jìn)行分析?！傮w剛度矩陣。——整個(gè)結(jié)構(gòu)上節(jié)點(diǎn)位移列陣——整個(gè)結(jié)構(gòu)上節(jié)點(diǎn)力列陣——整體平衡方程（整體剛度方程）分析過(guò)程：是將所有單元平衡方程組集在一起，形成總體平衡方程，引進(jìn)邊界條件后，求解整體節(jié)點(diǎn)位移向量。2.3整體分析結(jié)構(gòu)的整體分析就是將離散82平面問(wèn)題的整體剛度矩陣是由相關(guān)單元的單元?jiǎng)偠染仃囍械姆謮K矩陣集合而成，按節(jié)點(diǎn)編號(hào)對(duì)號(hào)入座，即剛度集成法。它的集成規(guī)律有下列幾點(diǎn)：1）先對(duì)每個(gè)單元求出其單元?jiǎng)偠染仃?，以分塊形式按節(jié)點(diǎn)編號(hào)順序排列。2）將單元?jiǎng)偠染仃嚁U(kuò)大階數(shù)為2n×2n，并將單元?jiǎng)偠染仃囍械姆謮K矩陣按局部碼與總碼的對(duì)應(yīng)關(guān)系，搬到擴(kuò)大后的矩陣中，形成單元貢獻(xiàn)矩陣。3）將所有單元貢獻(xiàn)矩陣同一位置上的分塊矩陣簡(jiǎn)單疊加成總體剛度矩陣中的一個(gè)子矩陣，各行各列都按以上步驟即形成總體剛度矩陣?？傮w剛度矩陣為2n×2n階，亦即n×n階分塊矩陣，n為節(jié)點(diǎn)總數(shù)。平面問(wèn)題的整體剛度矩陣是由相關(guān)單元的單元?jiǎng)偠染?3例2-2用剛度集成法求圖2-7所示結(jié)構(gòu)的整體剛度矩陣1）找出各單元局部碼與總碼的對(duì)應(yīng)關(guān)系解：2）分別寫出各個(gè)單元的分塊矩陣3）形成各單元貢獻(xiàn)陣?yán)?-2用剛度集成法求圖2-7所示結(jié)構(gòu)的整體剛度矩陣1）844）將擴(kuò)大后的4個(gè)單元貢獻(xiàn)陣相疊加，就得到總體剛度矩陣

整體剛度矩陣有以下一些性質(zhì)：1）整體剛度矩陣是對(duì)稱矩陣。2）整體剛度矩陣的主對(duì)角線上的元素總是正的。3）整體剛度矩陣是一個(gè)稀疏陣。4）整體剛度矩陣是一個(gè)奇異陣。4）將擴(kuò)大后的4個(gè)單元貢獻(xiàn)陣相疊加，就得到總體剛度矩陣整體85例2-3已知如圖a）所示的懸臂深梁，在右端面作用著均布拉力，其合力為P。采用如圖b）所示簡(jiǎn)單網(wǎng)格，設(shè)，厚度為t。試求節(jié)點(diǎn)位移。a）b）解：１）求出各單元?jiǎng)偠染仃嚴(yán)?-3已知如圖a）所示的懸臂深梁，在右端面作用著均布86２）形成各單元貢獻(xiàn)矩陣３）形成總體剛度矩陣

２）形成各單元貢獻(xiàn)矩陣３）形成總體剛度矩陣87ANSYS有限元基礎(chǔ)教程課件-第2章88４）形成整體載荷列陣

５）形成整體節(jié)點(diǎn)位移列陣

６）形成整體平衡方程

４）形成整體載荷列陣５）形成整體節(jié)點(diǎn)位移列陣６）形成整89７）引入邊界條件，求節(jié)點(diǎn)位移“化1置0法”處理：若已知節(jié)點(diǎn)在方向位移為零，則令

中的元素為1，而第行和列的其余元素都為零，中的第個(gè)元素變?yōu)榱恪Ｈ粢阎?jié)點(diǎn)在方向位移為零，則令

中的元素為1，而第行和列的其余元素都為零，中的第個(gè)元素變?yōu)榱?。７）引入邊界條件，求節(jié)點(diǎn)位移“化1置0法”處理：若已知節(jié)點(diǎn)902.4

有限元法解題過(guò)程與算例

有限元法的具體解題過(guò)程為：1）將結(jié)構(gòu)進(jìn)行離散化，包括單元?jiǎng)澐?、?jié)點(diǎn)編號(hào)、單元編號(hào)、節(jié)點(diǎn)坐標(biāo)計(jì)算、位移約束條件的確定。2）等效節(jié)點(diǎn)力的計(jì)算。3）剛度矩陣的計(jì)算。4）建立整體平衡方程，引入約束條件，求解節(jié)點(diǎn)位移。5）應(yīng)力計(jì)算。2.4有限元法解題過(guò)程與算例有限元法的具體解題過(guò)程91例2-4如圖2-20a所示兩端固支的矩形深梁，跨度為2a，梁高為a，厚度為t，已知E，，承受均布?jí)毫，試用有限元法求解此平面應(yīng)力問(wèn)題。a）b）圖2-20矩形深梁例2-4如圖2-20a所示兩端固支的矩形深梁，跨度為2a92解：利用對(duì)稱性，可取梁的一半分析，例如右半。1.劃分單元并準(zhǔn)備原始數(shù)據(jù)2.計(jì)算單元?jiǎng)偠染仃?31324解：利用對(duì)稱性，可取梁的一半分析，例如右半。1.劃分單元并933.集成整體剛度矩陣依照各單元局部編號(hào)與整體編號(hào)的對(duì)應(yīng)關(guān)系，兩個(gè)單元的貢獻(xiàn)矩陣分別為再集成整體剛度矩陣3.集成整體剛度矩陣依照各單元局部編號(hào)與整體944.處理載荷，形成整體平衡方程整體節(jié)點(diǎn)載荷列陣為組成結(jié)構(gòu)整體平衡方程4.處理載荷，形成整體平衡方程整體節(jié)點(diǎn)載荷列陣為組成結(jié)構(gòu)955.引入位移邊界條件，求解節(jié)點(diǎn)位移由于5.引入位移邊界條件，求解節(jié)點(diǎn)位移由于966.應(yīng)力計(jì)算在整體分析中求得節(jié)點(diǎn)位移之后，為了計(jì)算結(jié)構(gòu)上任意一點(diǎn)的應(yīng)變或應(yīng)力，應(yīng)該又返回到單元分析中去。計(jì)算單元①的應(yīng)力矩陣由整體節(jié)點(diǎn)位移向量獲取單元節(jié)點(diǎn)位移向量6.應(yīng)力計(jì)算計(jì)算單元①的應(yīng)力矩陣由整體節(jié)點(diǎn)位移向量獲取單元97計(jì)算應(yīng)力計(jì)算應(yīng)力982.5單元等效節(jié)點(diǎn)力一、單元自重三角形單元的厚度為t，重度為，面積為，自重沿y軸負(fù)方向。受自重載荷情形的等效節(jié)點(diǎn)力為單元重量的1/3。

圖2-14三角形單元

介紹幾種常用載荷作用下的等效節(jié)點(diǎn)力2.5單元等效節(jié)點(diǎn)力一、單元自重三角形單元99二、均布面力三角形單元的厚度為t，ij邊的長(zhǎng)度為，集度為圖2-15受均布面力三角形單元相當(dāng)于把作用于ij邊上的表面力按靜力等效平均分配到該邊兩端的節(jié)點(diǎn)上。二、均布面力三角形單元的厚度為t，ij邊的100三、線性分布面力三角形單元的厚度為t，ij邊的長(zhǎng)度為，表面力在點(diǎn)集度為圖2-16受線性分布面力三角形單元

相當(dāng)于將總載荷的2/3分配給點(diǎn)，1/3分配給點(diǎn)。

三、線性分布面力三角形單元的厚度為t，ij1012.6邊界條件的處理常用的且比較方便的做法是以某種方法引入已知的節(jié)點(diǎn)位移（包括零位移約束），而保持方程原有的數(shù)目不變，只是修

和中某些元素，以避免計(jì)算機(jī)存儲(chǔ)作大的變動(dòng)。2.6邊界條件的處理常用的且比較方便的做102

1、“化1置0法”----為節(jié)點(diǎn)總碼編號(hào)1、“化1置0法”----為節(jié)點(diǎn)總碼編號(hào)103設(shè)已知節(jié)點(diǎn)位移為，當(dāng)引進(jìn)上述已知節(jié)點(diǎn)位移后，方程變成例如：為說(shuō)明這一過(guò)程，現(xiàn)考察一個(gè)只有四個(gè)方程的簡(jiǎn)單例子設(shè)已知節(jié)點(diǎn)位移為，當(dāng)1042、“乘大數(shù)法”設(shè)已知節(jié)點(diǎn)位移為，當(dāng)引進(jìn)上述已知節(jié)點(diǎn)位移后，方程變成2、“乘大數(shù)法”設(shè)已知節(jié)點(diǎn)位移為1052.7計(jì)算結(jié)果的整理繞節(jié)點(diǎn)平均法：

ABCDEF把環(huán)繞該節(jié)點(diǎn)的各單元應(yīng)力加以平均，視為該節(jié)點(diǎn)的應(yīng)力。2.7計(jì)算結(jié)果的整理繞節(jié)點(diǎn)平均法：AB106兩單元平均法：把相鄰兩單元應(yīng)力的平均值作為公共邊中