塑性加工過(guò)程CAE02FEM基礎(chǔ)課件

第三章板料成形的有限元基礎(chǔ)§3.1序言§3.2桿和梁?jiǎn)卧?.2二維問(wèn)題§3.3有限元求解技術(shù)§3.4板殼單元第三章板料成形的有限元基礎(chǔ)§3.1序言§3.1序言§3.1.1

基本概念§3.1.2矩陣代數(shù)的回顧§3.1.3彈簧元§3.1序言§3.1.1基本概念§3.1序言有限元方法或有限元分析的基本思想就是將一個(gè)復(fù)雜的物體分解成許多易于處理的小片來(lái)處理。這種思路在日常生活和工程實(shí)際中經(jīng)常被使用?！?.1.1

基本概念搭積木游戲建筑§3.1序言有限元方法或有限元分析的基本思想就是§3.1.1基本概念圓面積的近似求解：注意：一個(gè)復(fù)雜的連續(xù)的物體可以被許多小片（單元）來(lái)近似代替?！?.1.1基本概念圓面積的近似求解：注意：一個(gè)復(fù)雜的連續(xù)工程中有限元的應(yīng)用機(jī)械、航空、土木、汽車(chē)工程結(jié)構(gòu)分析（靜態(tài)/動(dòng)態(tài)、線(xiàn)性/非線(xiàn)性）熱/液體流動(dòng)電磁地質(zhì)力學(xué)生物力學(xué)…齒輪式彈性軸接的模型§3.1.1基本概念工程中有限元的應(yīng)用齒輪式彈性軸接的模型§3.1.1基本概念有限元的簡(jiǎn)要發(fā)展歷史1943——Courant1956——Tuner,Clough,Martin,Topp(Stiffness)1960——Clough(“FiniteElement”平面問(wèn)題)1970s——在大型機(jī)上得到應(yīng)用1980s——在微機(jī)上得到應(yīng)用，前后處理軟件1990s——大型機(jī)構(gòu)系統(tǒng)分析§3.1.1基本概念有限元的簡(jiǎn)要發(fā)展歷史§3.1.1基本概念§3.1.1基本概念易拉罐的跌落測(cè)試§3.1.1基本概念易拉罐的跌落測(cè)試結(jié)構(gòu)分析中的有限元（過(guò)程）將結(jié)構(gòu)體劃分成小片（單元，節(jié)點(diǎn)）形成描繪物理特性的單元?jiǎng)偠葘卧M裝成一個(gè)整體結(jié)構(gòu)的近似方程組求解已經(jīng)引入位置物理量（位移）的方程組計(jì)算單元中用戶(hù)所關(guān)心的物理量（應(yīng)變，應(yīng)力）§3.1.1基本概念結(jié)構(gòu)分析中的有限元（過(guò)程）§3.1.1基本概念計(jì)算機(jī)前處理（建立有限元模型，載荷和約束）有限元求解器（組裝和求解系統(tǒng)方程）后處理（顯示計(jì)算結(jié)果）商業(yè)化有限元軟件包ANSYS，NASTRAN，ALGOR(通用目的)ABAQUS(非線(xiàn)性動(dòng)力分析)PATRAN，HyperMesh(前后處理)LS-DYNA(碰撞動(dòng)力分析)DynaForm(前后處理)…§3.1.1基本概念計(jì)算機(jī)§3.1.1基本概念本章有限元基礎(chǔ)課的目的理解有限元分析的基本原理和思路掌握本章中所涉及到的單元模型的推導(dǎo)和適用范圍對(duì)一個(gè)給定的問(wèn)題能夠建立適當(dāng)?shù)挠邢拊Ｐ湍軌蚪忉尣⒃u(píng)價(jià)有限元分析結(jié)果的優(yōu)劣（知道問(wèn)題的物理意思）明白有限元的局限性（不要錯(cuò)誤地應(yīng)用有限元—數(shù)值工具）§3.1.1基本概念本章有限元基礎(chǔ)課的目的§3.1.1基本概念線(xiàn)性代數(shù)方程組系統(tǒng)：矩陣形式：A稱(chēng)為n×n矩陣，x和b分別是n維的列向量§3.1.1矩陣代數(shù)的回顧線(xiàn)性代數(shù)方程組系統(tǒng)：§3.1.1矩陣代數(shù)的回顧矩陣的加法減法：矩陣的乘法：矩陣的轉(zhuǎn)置：對(duì)稱(chēng)矩陣：?jiǎn)挝痪仃嚕骸?.1.1矩陣代數(shù)的回顧矩陣的加法減法：§3.1.1矩陣代數(shù)的回顧矩陣行列式的值：奇異矩陣：如果detA=0，那么系統(tǒng)存在問(wèn)題（非唯一解，發(fā)散等）矩陣的逆：§3.1.1矩陣代數(shù)的回顧矩陣行列式的值：§3.1.1矩陣代數(shù)的回顧線(xiàn)性方程組系統(tǒng)的求解技術(shù)：高斯消去法迭代法正定矩陣：對(duì)于所有非零向量X，有XTAX>0，A為正定矩陣正定矩陣為非奇異矩陣矩陣的導(dǎo)數(shù)和積分：§3.1.1矩陣代數(shù)的回顧線(xiàn)性方程組系統(tǒng)的求解技術(shù)：§3.1.1矩陣代數(shù)的回顧彈簧元：§3.1.3彈簧單元§3.1.3彈簧單元考慮彈簧元的力的平衡條件：節(jié)點(diǎn)i：節(jié)點(diǎn)j：矩陣形式：注意：K為對(duì)稱(chēng)矩陣，K是非奇異的還是奇異的？§3.1.3彈簧單元考慮彈簧元的力的平衡條件：§3.1.3彈簧單元彈簧元系統(tǒng)：?jiǎn)卧?：?jiǎn)卧?：§3.1.3彈簧單元彈簧元系統(tǒng)：§3.1.3彈簧單元對(duì)整個(gè)系統(tǒng)進(jìn)行單元?jiǎng)偠染仃嚨慕M裝：考慮節(jié)點(diǎn)力的平衡條件：節(jié)點(diǎn)1：節(jié)點(diǎn)2：節(jié)點(diǎn)3：矩陣形式：K為該彈簧系統(tǒng)的剛度矩陣（結(jié)構(gòu)矩陣）§3.1.3彈簧單元對(duì)整個(gè)系統(tǒng)進(jìn)行單元?jiǎng)偠染仃嚨慕M裝：§3.1.3彈簧單元單元?jiǎng)偠染仃嚵硪环N組裝方法：分別擴(kuò)大單元1和2的剛度矩陣這與根據(jù)節(jié)點(diǎn)力平衡得出的矩陣是一樣的。§3.1.3彈簧單元單元?jiǎng)偠染仃嚵硪环N組裝方法：§3.1.3彈簧單元引入邊界條件和力的條件：假設(shè)u1=0，F(xiàn)2=F3=P，那么未知量為：u2，u3和F1求解方程可得：§3.1.3彈簧單元引入邊界條件和力的條件：§3.1.3彈簧單元檢查計(jì)算結(jié)果：結(jié)構(gòu)變形后的形狀外力平衡有關(guān)彈簧元的注意事項(xiàng)：適于剛度分析的計(jì)算不適合用于彈簧本身的應(yīng)力分析計(jì)算在彈簧元的橫向是否具有剛度，彈簧元是否具有扭轉(zhuǎn)剛度§3.1.3彈簧單元檢查計(jì)算結(jié)果：§3.1.3彈簧單元例子1：已知：求：(a)整體剛度矩陣(b)節(jié)點(diǎn)2和3的位移(c)節(jié)點(diǎn)1和4的支反力(d)彈簧2的力§3.1.3彈簧單元例子1：§3.1.3彈簧單元（a）問(wèn)：分別求出單元?jiǎng)偠染仃嚕簡(jiǎn)卧?：?jiǎn)卧?：?jiǎn)卧?：組裝后：§3.1.3彈簧單元（a）問(wèn)：分別求出單元?jiǎng)偠染仃嚕骸?.1.3彈簧單元組裝后：注意到整體剛度矩陣是對(duì)稱(chēng)并帶狀分布的。該系統(tǒng)的平衡方程為：§3.1.3彈簧單元組裝后：§3.1.3彈簧單元（b）問(wèn)：將邊界條件（u1=u4=0）應(yīng)用到平衡方程中，去掉1行1列，4行4列后：求解得：（c）問(wèn)：從平衡方程組中的1和4可得：§3.1.3彈簧單元（b）問(wèn)：§3.1.3彈簧單元（d）問(wèn)：彈簧元2的平衡方程為：式中i=2，j=3，可以計(jì)算出彈簧力為：§3.1.3彈簧單元（d）問(wèn)：§3.1.3彈簧單元例子2：?jiǎn)栴}描述：對(duì)上述一個(gè)具有任意彈簧元節(jié)點(diǎn)和單元的系統(tǒng)，求其整體剛度矩陣?！?.1.3彈簧單元例子2：§3.1.3彈簧單元單元拓?fù)潢P(guān)系：上表中為每個(gè)彈簧元的局部節(jié)點(diǎn)號(hào)和整體節(jié)點(diǎn)號(hào)的對(duì)應(yīng)關(guān)系。然后依次求出每個(gè)彈簧元的剛度矩陣：§3.1.3彈簧單元單元拓?fù)潢P(guān)系：§3.1.3彈簧單元組裝后：整體剛度矩陣是對(duì)稱(chēng)并帶狀分布的。§3.1.3彈簧單元§3.1.3彈簧單元§3.2桿和梁?jiǎn)卧?.2.1

線(xiàn)性靜力分析§3.2.2桿單元§3.2.3梁?jiǎn)卧?.2桿和梁?jiǎn)卧?.2.1線(xiàn)性靜力分析大部分結(jié)構(gòu)分析問(wèn)題都可以看作是線(xiàn)性靜力分析問(wèn)題，它們都基于以下假設(shè)：1.小變形（加載方式不會(huì)因?yàn)樽冃味淖儯?.線(xiàn)彈性材料（不存在塑性和破裂）3.靜力載荷（在結(jié)構(gòu)上的載荷是慢速平穩(wěn)地施加上去的）線(xiàn)性靜力分析可以解決結(jié)構(gòu)分析問(wèn)題中大部分的問(wèn)題，對(duì)于大多數(shù)的結(jié)構(gòu)分析問(wèn)題，靜力分析可以得到一個(gè)近似的結(jié)果。線(xiàn)性靜力分析是非線(xiàn)性分析的基礎(chǔ)?！?.2.1線(xiàn)性靜力分析大部分結(jié)構(gòu)分析問(wèn)題都可以看作是線(xiàn)性靜力分析問(wèn)題，它們都基于以桿單元：§3.2.2桿單元桿單元：§3.2.2桿單元?jiǎng)偠染仃嚒苯臃ǎ杭僭O(shè)位移u沿著桿的軸向線(xiàn)性分布：K即為桿單元的剛度系數(shù)，桿單元和彈簧元類(lèi)似?！?.2.2桿單元?jiǎng)偠染仃嚒苯臃ǎ骸?.2.2桿單元單元?jiǎng)偠染仃嚍椋簡(jiǎn)卧钠胶夥匠探M為：節(jié)點(diǎn)自由度：對(duì)一維的桿單元，每個(gè)節(jié)點(diǎn)就只有一個(gè)自由度剛度矩陣K中系數(shù)的物理意義：K中第j列的系數(shù)表示在節(jié)點(diǎn)j施加單位位移而其他節(jié)點(diǎn)固定不動(dòng)的時(shí)候，桿上所承受的力?！?.2.2桿單元單元?jiǎng)偠染仃嚍椋骸?.2.2桿單元?jiǎng)偠染仃嚒?guī)推導(dǎo)方法：定義兩個(gè)形函數(shù)：位移u沿著桿的軸向線(xiàn)性分布：可得應(yīng)變?yōu)椋築為單元的應(yīng)變-位移矩陣：§3.2.2桿單元?jiǎng)偠染仃嚒?guī)推導(dǎo)方法：§3.2.2桿單元單元應(yīng)力為：桿單元上的應(yīng)變能為：桿單元2個(gè)節(jié)點(diǎn)所做的功為：對(duì)于保守系統(tǒng)，U=W

，故有：§3.2.2桿單元單元應(yīng)力為：§3.2.2桿單元上式等價(jià)于：其中k即為單元?jiǎng)偠染仃嚿鲜龇椒ň褪钦?guī)的推導(dǎo)過(guò)程，該方法也可以用來(lái)推導(dǎo)其他類(lèi)型的單元?jiǎng)偠染仃?。單元?jiǎng)偠染仃囈部梢酝ㄟ^(guò)其他嚴(yán)格的方法獲得，比如最小勢(shì)能原理，伽遼金方法等。于是我們可以得到桿單元的剛度矩陣：上式結(jié)果和直接法的結(jié)果一樣?！?.2.2桿單元上式等價(jià)于：§3.2.2桿單元例子：?jiǎn)栴}描述：在節(jié)點(diǎn)2處施加F，求桿1和桿2上的應(yīng)力。求解方法：用1-D桿單元。單元1和單元2的剛度矩陣分別為：假設(shè)在節(jié)點(diǎn)2處是用一個(gè)無(wú)摩擦的鉸鏈將桿1和桿2所連接?！?.2.2桿單元例子：§3.2.2桿單元組裝后：載荷和邊界條件為：刪除第1行第1列和第3行第3列后得：§3.2.2桿單元組裝后：§3.2.2桿單元最后可以得到單元1內(nèi)的應(yīng)力：同理，可以得到單元2內(nèi)的應(yīng)力：負(fù)號(hào)表示單元2所受的是壓應(yīng)力。§3.2.2桿單元最后可以得到單元1內(nèi)的應(yīng)力：§3.2.2桿單元注意：1)在這個(gè)例子中，根據(jù)一維線(xiàn)性理論所計(jì)算出來(lái)的單元1和2中的應(yīng)力是精確解，所以如果我們將單元細(xì)化不會(huì)提高精度。2)如果對(duì)于階梯桿結(jié)構(gòu)，其中的A要采用橫截面的平均面積。3)為了得到單元1和2中的應(yīng)力，我們首先要得到節(jié)點(diǎn)的位移，因?yàn)椴捎玫幕谖灰茍?chǎng)的有限元法。§3.2.2桿單元注意：§3.2.2桿單元二維空間上的桿單元：注意：在線(xiàn)彈性理論的前提下，橫向位移對(duì)桿單元的拉伸沒(méi)有貢獻(xiàn)?！?.2.2桿單元二維空間上的桿單元：§3.2.2桿單元坐標(biāo)變換：矩陣形式：§3.2.2桿單元坐標(biāo)變換：§3.2.2桿單元對(duì)于2個(gè)節(jié)點(diǎn)的桿單元有：節(jié)點(diǎn)力也采用相同的變換方法：二維空間下的局部坐標(biāo)系下的剛度矩陣：§3.2.2桿單元對(duì)于2個(gè)節(jié)點(diǎn)的桿單元有：§3.2.2桿單元擴(kuò)大后：矩陣形式：K為整體坐標(biāo)系下的剛度矩陣§3.2.2桿單元擴(kuò)大后：§3.2.2桿單元K的顯式表達(dá)式為：?jiǎn)卧獞?yīng)力：§3.2.2桿單元K的顯式表達(dá)式為：§3.2.2桿單元三維空間上的桿單元：先在局部坐標(biāo)系下求出單元?jiǎng)偠染仃嚕缓笤俎D(zhuǎn)換并組裝到整體坐標(biāo)系下進(jìn)行計(jì)算?！?.2.2桿單元三維空間上的桿單元：§3.2.2桿單元簡(jiǎn)單的平面梁?jiǎn)卧夯镜牧豪碚摚骸?.2.3梁?jiǎn)卧?jiǎn)單的平面梁?jiǎn)卧骸?.2.3梁?jiǎn)卧?.3二維問(wèn)題§3.3.1

基本理論回顧§3.3.2二維問(wèn)題的有限元§3.3二維問(wèn)題§3.3.1基本理論回顧應(yīng)力和應(yīng)變：在特定的條件下，應(yīng)力應(yīng)變狀態(tài)可以簡(jiǎn)化，因此，一般的三維問(wèn)題可以簡(jiǎn)化成二維問(wèn)題來(lái)分析§3.3.1基本理論回顧應(yīng)力和應(yīng)變：§3.3.1基本理論回顧平面應(yīng)力：平面應(yīng)變：§3.3.1基本理論回顧平面應(yīng)力：§3.3.1基本理論回顧本構(gòu)關(guān)系（平面應(yīng)力）：本構(gòu)關(guān)系（平面應(yīng)變）：§3.3.1基本理論回顧本構(gòu)關(guān)系（平面應(yīng)力）：§3.3.1基本理論回顧應(yīng)變與位移的關(guān)系：邊界條件：§3.3.1基本理論回顧應(yīng)變與位移的關(guān)系：§3.3.1基本理論回顧常應(yīng)變單元（CSTorT3）位移函數(shù)：§3.3.2二維問(wèn)題的有限元常應(yīng)變單元（CSTorT3）§3.3.2二維問(wèn)題的有限位移插值函數(shù)：§3.3.2二維問(wèn)題的有限元位移插值函數(shù)：§3.3.2二維問(wèn)題的有限元應(yīng)變：CST單元?jiǎng)偠染仃嚍椋簍為單元厚度，k為一個(gè)6×6的對(duì)稱(chēng)矩陣?！?.3.2二維問(wèn)題的有限元應(yīng)變：§3.3.2二維問(wèn)題的有限元雙線(xiàn)性四邊形單元（Q4）位移函數(shù)：§3.3.2二維問(wèn)題的有限元雙線(xiàn)性四邊形單元（Q4）§3.3.2二維問(wèn)題的有限元例子：§3.3.2二維問(wèn)題的有限元例子：§3.3.2二維問(wèn)題的有限元網(wǎng)格劃分：§3.3.2二維問(wèn)題的有限元網(wǎng)格劃分：§3.3.2二維問(wèn)題的有限元計(jì)算結(jié)果：§3.3.2二維問(wèn)題的有限元計(jì)算結(jié)果：§3.3.2二維問(wèn)題的有限元載荷的變換：集中力(點(diǎn)載荷)，面力(壓力載荷)，體力(重力)為三種主要的外載荷。面力和體力都需要經(jīng)過(guò)變換以后才能加到有限元中。變換的基本思想就是：等效功?！?.3.2二維問(wèn)題的有限元載荷的變換：§3.3.2二維問(wèn)題的有限元其中，t為單元厚度，L為單元邊長(zhǎng)，un為垂直于邊界AB的位移分量。對(duì)于Q4單元來(lái)說(shuō)，有：§3.3.2二維問(wèn)題的有限元§3.3.2二維問(wèn)題的有限元應(yīng)力計(jì)算：?jiǎn)卧獌?nèi)的應(yīng)力計(jì)算公式為：B為應(yīng)變與節(jié)點(diǎn)位移之間的關(guān)系矩陣，d為計(jì)算后的節(jié)點(diǎn)位移向量。我們可以計(jì)算出單元內(nèi)部任意一點(diǎn)的應(yīng)力，包括單元中心和節(jié)點(diǎn)位置。應(yīng)力分布的等值線(xiàn)可以通過(guò)后處理軟件顯示出來(lái)?！?.3.2二維問(wèn)題的有限元應(yīng)力計(jì)算：§3.3.2二維問(wèn)題的有限元討論：1）需要了解每個(gè)單元類(lèi)型的特性：T3和Q4：線(xiàn)性位移，常量應(yīng)變和應(yīng)力T6和Q8：二次位移，線(xiàn)性的應(yīng)變和應(yīng)力2）對(duì)于一個(gè)特定的問(wèn)題選擇一種合適的單元毫無(wú)疑問(wèn)，應(yīng)該盡可能地采用高階單元和細(xì)網(wǎng)格。3）要避免單元具有大形狀比和尖角：§3.3.2二維問(wèn)題的有限元討論：§3.3.2二維問(wèn)題的有限元討論：4）單元之間要互相連通：在有限元模型中不能存在間隙或自由單元?！?.3.2二維問(wèn)題的有限元討論：§3.3.2二維問(wèn)題的有限元§3.3有限元求解技術(shù)§3.3.1

方程組求解§3.3.2有限元方法的實(shí)質(zhì)§3.3.3數(shù)值誤差§3.3.4有限元求解的收斂性§3.3有限元求解技術(shù)§3.3.1方程組求解直接法（高斯消去法）：求解時(shí)間是和NB2成正比（N為矩陣的維數(shù)，B為矩陣的帶寬）適合小到中等規(guī)模的問(wèn)題，小帶寬的問(wèn)題。對(duì)于多重載荷的情況，容易處理。迭代法：求解時(shí)間事先不可預(yù)知降低對(duì)存儲(chǔ)空間的要求適合大型規(guī)模問(wèn)題的求解（帶寬可以很大，收斂快）對(duì)于不同載荷的工況需要重新求解。§3.3.1方程組求解直接法（高斯消去法）：§3.3.1方程組求解高斯消去法例子：回代求解：§3.3.1方程組求解高斯消去法例子：§3.3.1方程組求解迭代法（高斯-賽德?tīng)枺├樱簭囊粋€(gè)初始解X(0)開(kāi)始按照以下公式計(jì)算：一直到X滿(mǎn)足以下條件，迭代終止。其中ε為收斂控制的容忍誤差?！?.3.1方程組求解迭代法（高斯-賽德?tīng)枺├樱骸?.3.1方程組求解有限元模型—是基于很多近似以后對(duì)實(shí)際結(jié)構(gòu)的一個(gè)數(shù)學(xué)模型。實(shí)際結(jié)構(gòu)—具有無(wú)限個(gè)節(jié)點(diǎn)，所以具有無(wú)限個(gè)自由度。有限元模型—具有有限個(gè)節(jié)點(diǎn)，所以具有有限個(gè)自由度。位移場(chǎng)是被有限個(gè)節(jié)點(diǎn)位移的值所控制的：§3.3.2有限元方法的實(shí)質(zhì)有限元模型—是基于很多近似以后對(duì)實(shí)際結(jié)構(gòu)的一個(gè)數(shù)學(xué)模型?！?剛度效果：有限元模型比實(shí)際結(jié)構(gòu)要更剛一些。一般來(lái)說(shuō)，位移解要比實(shí)際情況偏小一些。所以，有限元計(jì)算的位移是精確解的一個(gè)下限：§3.3.2有限元方法的實(shí)質(zhì)剛度效果：§3.3.2有限元方法的實(shí)質(zhì)誤差≠錯(cuò)誤（模型和計(jì)算）。誤差的類(lèi)型：模型誤差（梁，板，…，理論）離散誤差（有限，分段，…

）數(shù)值誤差（在求解有限元方程組的時(shí)候）例子（數(shù)值誤差）：§3.3.3數(shù)值誤差誤差≠錯(cuò)誤（模型和計(jì)算）。§3.3.3數(shù)值誤差有限元方程組為：如果K2<<K1的話(huà)，方程組是奇異的，病態(tài)的：§3.3.3數(shù)值誤差有限元方程組為：§3.3.3數(shù)值誤差如果K2>>K1的話(huà)，方程組是非病態(tài)的：在有限元模型中不同部分的剛度相差很大的話(huà)有可能導(dǎo)致有限元方程組的病態(tài)，產(chǎn)生很大的誤差。病態(tài)方程組的情況下，很小的輸入變化量（右端項(xiàng)向量）求解時(shí)會(huì)引起很大的變化量?！?.3.3數(shù)值誤差如果K2>>K1的話(huà)，方程組是非病態(tài)的：§3.3.3數(shù)隨著有限元模型中的網(wǎng)格不斷地被加密和細(xì)化，有限元的計(jì)算結(jié)果將收斂于被求解的問(wèn)題的數(shù)學(xué)模型的精確解。有限元網(wǎng)格自適應(yīng)類(lèi)型：h-refinement:減小單元尺寸（h指單元的尺寸）。p-refinement:增加單元多項(xiàng)式插值函數(shù)的階次（p指更高階次的多項(xiàng)式插值函數(shù)）。r-refinement:重新劃分有限元網(wǎng)格。hp-refinement:h方法和p方法互相結(jié)合的方法?！?.3.4有限元求解的收斂性隨著有限元模型中的網(wǎng)格不斷地被加密和細(xì)化，有限元的計(jì)算結(jié)果將網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性§3.4板殼單元§3.4.1

板理論§3.4.2板單元§3.4.3殼和殼單元§3.4板殼單元§3.4.1板理論板理論：平板橫向加載彎曲效應(yīng)為主注意：應(yīng)用范圍：剪力墻地板架子等§3.4.1板理論板理論：§3.4.1板理論板單元的力和彎矩：應(yīng)力：§3.4.1板理論板單元的力和彎矩：§3.4.1板理論板單元的力和應(yīng)力之間的關(guān)系：?jiǎn)挝婚L(zhǎng)度的彎矩：?jiǎn)挝婚L(zhǎng)度的扭矩：?jiǎn)挝婚L(zhǎng)度的剪力：最大彎曲應(yīng)力：最大彎曲應(yīng)力總是在在中性層沒(méi)有彎曲應(yīng)力（與梁?jiǎn)卧?lèi)似）§3.4.1板理論板單元的力和應(yīng)力之間的關(guān)系：§3.4.1板理論薄板理論（Kichhoff板理論）：直法線(xiàn)假設(shè)（與梁?jiǎn)卧?lèi)似）：與中性層垂直的直法線(xiàn)在變形后仍然與中性層保持垂直，也就是說(shuō)沒(méi)有剪切變形：位移：§3.4.1板理論薄板理論（Kichhoff板理論）：§3.4.1板理論應(yīng)變：注意：在中性層沒(méi)有拉伸變形。應(yīng)力（平面應(yīng)力）：剪力和彎矩：§3.4.1板理論應(yīng)變：厚板理論（Mindlin板理論）：如果板厚不夠薄，即t/L>1/10，應(yīng)該采用厚板理論。直線(xiàn)假設(shè)：該理論認(rèn)為在變形過(guò)程中截面角度發(fā)生改變，即：與中性層垂直的直線(xiàn)在變形后與中性層不再保持垂直，但仍然是直線(xiàn)。§3.4.1板理論厚板理論（Mindlin板理論）：§3.4.1板理論新的獨(dú)立變量：θx和θy：分別為直線(xiàn)沿x軸和y軸的旋轉(zhuǎn)角度，變形前該直線(xiàn)與中性層垂直。新的關(guān)系：注意：如果我們引入以下條件，厚板理論又可以退化成薄板理論。§3.4.1板理論新的獨(dú)立變量：§3.4.1板理論Kichhoff板單元：四節(jié)點(diǎn)四邊形單元每個(gè)節(jié)點(diǎn)的自由度：每個(gè)單元的z方向的位移為：§3.4.2板單元Kichhoff板單元：§3.4.2板單元Kichhoff板單元是一種不協(xié)調(diào)單元，但是其剛度矩陣仍舊可以寫(xiě)成：B為應(yīng)變與位移的關(guān)系矩陣，E為應(yīng)力應(yīng)變關(guān)系矩陣?！?.4.2板單元Kichhoff板單元是一種不協(xié)調(diào)單元，但是其剛度矩陣仍舊Mindlin板單元：4-nodes四邊形單元8-nodes四邊形單元每個(gè)節(jié)點(diǎn)的自由度：每個(gè)單元有：三個(gè)獨(dú)立的位移場(chǎng)對(duì)Q4，w(x,y)是線(xiàn)性的，對(duì)Q8，w(x,y)是二次。§3.4.2板單元Mindlin板單元：§3.4.2板單元離散的Kichhoff板單元（DKT）：三角形板單元，推導(dǎo)開(kāi)始于6-nodes三角形單元：每個(gè)角節(jié)點(diǎn)的自由度：每個(gè)中節(jié)點(diǎn)的自由度：總自由度為21。然后引入邊界條件：減少自由度。§3.4.2板單元離散的Kichhoff板單元（DKT）：§3.4.2板單可以得到：每個(gè)節(jié)點(diǎn)的自由度：總自由度為9（DKT單元）。DKT板單元的w(x,y)也是不協(xié)調(diào)位移函數(shù)，DKT板單元收斂速度快，有效?！?.4.2板單元可以得到：§3.4.2板單元?dú)ぃ壕哂袕澢娴谋〗Y(jié)構(gòu)例子：雞蛋殼（自然界的奇跡）容器，管道，罐汽車(chē)車(chē)身屋頂，穹頂結(jié)構(gòu)的建筑，等等§3.4.1殼與殼單元?dú)ぃ壕哂袕澢娴谋〗Y(jié)構(gòu)§3.4.1殼與殼單元?dú)ぶ械牧Γ耗ち?彎矩（板只能承受彎矩）例子：圓柱形容器§3.4.1殼與殼單元?dú)ぶ械牧Γ骸?.4.1殼與殼單元?dú)だ碚摚罕だ碚摵駳だ碚撛诹W(xué)中，殼理論是一種在列式和分析方面最復(fù)雜的理論之一。（前蘇聯(lián)的科學(xué)家貢獻(xiàn)很大）：工程≠技能需要很強(qiáng)的分析技巧§3.4.1殼與殼單元§3.4.1殼與殼單元平板殼單元（應(yīng)用最廣泛）：比較：桿單元+簡(jiǎn)單梁?jiǎn)卧?一般的梁?jiǎn)卧卧?平板單元=平板殼單元每個(gè)節(jié)點(diǎn)的自由度：Q4或Q8殼單元：§3.4.1殼與殼單元平板殼單元（應(yīng)用最廣泛）：§3.4.1殼與殼單元彎曲殼單元：基于殼理論最一般的殼單元平板殼單元和板單元是其子集列式和推導(dǎo)極其復(fù)雜§3.4.1殼與殼單元彎曲殼單元：§3.4.1殼與殼單元第三章板料成形的有限元基礎(chǔ)§3.1序言§3.2桿和梁?jiǎn)卧?.2二維問(wèn)題§3.3有限元求解技術(shù)§3.4板殼單元第三章板料成形的有限元基礎(chǔ)§3.1序言§3.1序言§3.1.1

基本概念§3.1.2矩陣代數(shù)的回顧§3.1.3彈簧元§3.1序言§3.1.1基本概念§3.1序言有限元方法或有限元分析的基本思想就是將一個(gè)復(fù)雜的物體分解成許多易于處理的小片來(lái)處理。這種思路在日常生活和工程實(shí)際中經(jīng)常被使用?！?.1.1

基本概念搭積木游戲建筑§3.1序言有限元方法或有限元分析的基本思想就是§3.1.1基本概念圓面積的近似求解：注意：一個(gè)復(fù)雜的連續(xù)的物體可以被許多小片（單元）來(lái)近似代替。§3.1.1基本概念圓面積的近似求解：注意：一個(gè)復(fù)雜的連續(xù)工程中有限元的應(yīng)用機(jī)械、航空、土木、汽車(chē)工程結(jié)構(gòu)分析（靜態(tài)/動(dòng)態(tài)、線(xiàn)性/非線(xiàn)性）熱/液體流動(dòng)電磁地質(zhì)力學(xué)生物力學(xué)…齒輪式彈性軸接的模型§3.1.1基本概念工程中有限元的應(yīng)用齒輪式彈性軸接的模型§3.1.1基本概念有限元的簡(jiǎn)要發(fā)展歷史1943——Courant1956——Tuner,Clough,Martin,Topp(Stiffness)1960——Clough(“FiniteElement”平面問(wèn)題)1970s——在大型機(jī)上得到應(yīng)用1980s——在微機(jī)上得到應(yīng)用，前后處理軟件1990s——大型機(jī)構(gòu)系統(tǒng)分析§3.1.1基本概念有限元的簡(jiǎn)要發(fā)展歷史§3.1.1基本概念§3.1.1基本概念易拉罐的跌落測(cè)試§3.1.1基本概念易拉罐的跌落測(cè)試結(jié)構(gòu)分析中的有限元（過(guò)程）將結(jié)構(gòu)體劃分成小片（單元，節(jié)點(diǎn)）形成描繪物理特性的單元?jiǎng)偠葘卧M裝成一個(gè)整體結(jié)構(gòu)的近似方程組求解已經(jīng)引入位置物理量（位移）的方程組計(jì)算單元中用戶(hù)所關(guān)心的物理量（應(yīng)變，應(yīng)力）§3.1.1基本概念結(jié)構(gòu)分析中的有限元（過(guò)程）§3.1.1基本概念計(jì)算機(jī)前處理（建立有限元模型，載荷和約束）有限元求解器（組裝和求解系統(tǒng)方程）后處理（顯示計(jì)算結(jié)果）商業(yè)化有限元軟件包ANSYS，NASTRAN，ALGOR(通用目的)ABAQUS(非線(xiàn)性動(dòng)力分析)PATRAN，HyperMesh(前后處理)LS-DYNA(碰撞動(dòng)力分析)DynaForm(前后處理)…§3.1.1基本概念計(jì)算機(jī)§3.1.1基本概念本章有限元基礎(chǔ)課的目的理解有限元分析的基本原理和思路掌握本章中所涉及到的單元模型的推導(dǎo)和適用范圍對(duì)一個(gè)給定的問(wèn)題能夠建立適當(dāng)?shù)挠邢拊Ｐ湍軌蚪忉尣⒃u(píng)價(jià)有限元分析結(jié)果的優(yōu)劣（知道問(wèn)題的物理意思）明白有限元的局限性（不要錯(cuò)誤地應(yīng)用有限元—數(shù)值工具）§3.1.1基本概念本章有限元基礎(chǔ)課的目的§3.1.1基本概念線(xiàn)性代數(shù)方程組系統(tǒng)：矩陣形式：A稱(chēng)為n×n矩陣，x和b分別是n維的列向量§3.1.1矩陣代數(shù)的回顧線(xiàn)性代數(shù)方程組系統(tǒng)：§3.1.1矩陣代數(shù)的回顧矩陣的加法減法：矩陣的乘法：矩陣的轉(zhuǎn)置：對(duì)稱(chēng)矩陣：?jiǎn)挝痪仃嚕骸?.1.1矩陣代數(shù)的回顧矩陣的加法減法：§3.1.1矩陣代數(shù)的回顧矩陣行列式的值：奇異矩陣：如果detA=0，那么系統(tǒng)存在問(wèn)題（非唯一解，發(fā)散等）矩陣的逆：§3.1.1矩陣代數(shù)的回顧矩陣行列式的值：§3.1.1矩陣代數(shù)的回顧線(xiàn)性方程組系統(tǒng)的求解技術(shù)：高斯消去法迭代法正定矩陣：對(duì)于所有非零向量X，有XTAX>0，A為正定矩陣正定矩陣為非奇異矩陣矩陣的導(dǎo)數(shù)和積分：§3.1.1矩陣代數(shù)的回顧線(xiàn)性方程組系統(tǒng)的求解技術(shù)：§3.1.1矩陣代數(shù)的回顧彈簧元：§3.1.3彈簧單元§3.1.3彈簧單元考慮彈簧元的力的平衡條件：節(jié)點(diǎn)i：節(jié)點(diǎn)j：矩陣形式：注意：K為對(duì)稱(chēng)矩陣，K是非奇異的還是奇異的？§3.1.3彈簧單元考慮彈簧元的力的平衡條件：§3.1.3彈簧單元彈簧元系統(tǒng)：?jiǎn)卧?：?jiǎn)卧?：§3.1.3彈簧單元彈簧元系統(tǒng)：§3.1.3彈簧單元對(duì)整個(gè)系統(tǒng)進(jìn)行單元?jiǎng)偠染仃嚨慕M裝：考慮節(jié)點(diǎn)力的平衡條件：節(jié)點(diǎn)1：節(jié)點(diǎn)2：節(jié)點(diǎn)3：矩陣形式：K為該彈簧系統(tǒng)的剛度矩陣（結(jié)構(gòu)矩陣）§3.1.3彈簧單元對(duì)整個(gè)系統(tǒng)進(jìn)行單元?jiǎng)偠染仃嚨慕M裝：§3.1.3彈簧單元單元?jiǎng)偠染仃嚵硪环N組裝方法：分別擴(kuò)大單元1和2的剛度矩陣這與根據(jù)節(jié)點(diǎn)力平衡得出的矩陣是一樣的?！?.1.3彈簧單元單元?jiǎng)偠染仃嚵硪环N組裝方法：§3.1.3彈簧單元引入邊界條件和力的條件：假設(shè)u1=0，F(xiàn)2=F3=P，那么未知量為：u2，u3和F1求解方程可得：§3.1.3彈簧單元引入邊界條件和力的條件：§3.1.3彈簧單元檢查計(jì)算結(jié)果：結(jié)構(gòu)變形后的形狀外力平衡有關(guān)彈簧元的注意事項(xiàng)：適于剛度分析的計(jì)算不適合用于彈簧本身的應(yīng)力分析計(jì)算在彈簧元的橫向是否具有剛度，彈簧元是否具有扭轉(zhuǎn)剛度§3.1.3彈簧單元檢查計(jì)算結(jié)果：§3.1.3彈簧單元例子1：已知：求：(a)整體剛度矩陣(b)節(jié)點(diǎn)2和3的位移(c)節(jié)點(diǎn)1和4的支反力(d)彈簧2的力§3.1.3彈簧單元例子1：§3.1.3彈簧單元（a）問(wèn)：分別求出單元?jiǎng)偠染仃嚕簡(jiǎn)卧?：?jiǎn)卧?：?jiǎn)卧?：組裝后：§3.1.3彈簧單元（a）問(wèn)：分別求出單元?jiǎng)偠染仃嚕骸?.1.3彈簧單元組裝后：注意到整體剛度矩陣是對(duì)稱(chēng)并帶狀分布的。該系統(tǒng)的平衡方程為：§3.1.3彈簧單元組裝后：§3.1.3彈簧單元（b）問(wèn)：將邊界條件（u1=u4=0）應(yīng)用到平衡方程中，去掉1行1列，4行4列后：求解得：（c）問(wèn)：從平衡方程組中的1和4可得：§3.1.3彈簧單元（b）問(wèn)：§3.1.3彈簧單元（d）問(wèn)：彈簧元2的平衡方程為：式中i=2，j=3，可以計(jì)算出彈簧力為：§3.1.3彈簧單元（d）問(wèn)：§3.1.3彈簧單元例子2：?jiǎn)栴}描述：對(duì)上述一個(gè)具有任意彈簧元節(jié)點(diǎn)和單元的系統(tǒng)，求其整體剛度矩陣?！?.1.3彈簧單元例子2：§3.1.3彈簧單元單元拓?fù)潢P(guān)系：上表中為每個(gè)彈簧元的局部節(jié)點(diǎn)號(hào)和整體節(jié)點(diǎn)號(hào)的對(duì)應(yīng)關(guān)系。然后依次求出每個(gè)彈簧元的剛度矩陣：§3.1.3彈簧單元單元拓?fù)潢P(guān)系：§3.1.3彈簧單元組裝后：整體剛度矩陣是對(duì)稱(chēng)并帶狀分布的?！?.1.3彈簧單元§3.1.3彈簧單元§3.2桿和梁?jiǎn)卧?.2.1

線(xiàn)性靜力分析§3.2.2桿單元§3.2.3梁?jiǎn)卧?.2桿和梁?jiǎn)卧?.2.1線(xiàn)性靜力分析大部分結(jié)構(gòu)分析問(wèn)題都可以看作是線(xiàn)性靜力分析問(wèn)題，它們都基于以下假設(shè)：1.小變形（加載方式不會(huì)因?yàn)樽冃味淖儯?.線(xiàn)彈性材料（不存在塑性和破裂）3.靜力載荷（在結(jié)構(gòu)上的載荷是慢速平穩(wěn)地施加上去的）線(xiàn)性靜力分析可以解決結(jié)構(gòu)分析問(wèn)題中大部分的問(wèn)題，對(duì)于大多數(shù)的結(jié)構(gòu)分析問(wèn)題，靜力分析可以得到一個(gè)近似的結(jié)果。線(xiàn)性靜力分析是非線(xiàn)性分析的基礎(chǔ)。§3.2.1線(xiàn)性靜力分析大部分結(jié)構(gòu)分析問(wèn)題都可以看作是線(xiàn)性靜力分析問(wèn)題，它們都基于以桿單元：§3.2.2桿單元桿單元：§3.2.2桿單元?jiǎng)偠染仃嚒苯臃ǎ杭僭O(shè)位移u沿著桿的軸向線(xiàn)性分布：K即為桿單元的剛度系數(shù)，桿單元和彈簧元類(lèi)似?！?.2.2桿單元?jiǎng)偠染仃嚒苯臃ǎ骸?.2.2桿單元單元?jiǎng)偠染仃嚍椋簡(jiǎn)卧钠胶夥匠探M為：節(jié)點(diǎn)自由度：對(duì)一維的桿單元，每個(gè)節(jié)點(diǎn)就只有一個(gè)自由度剛度矩陣K中系數(shù)的物理意義：K中第j列的系數(shù)表示在節(jié)點(diǎn)j施加單位位移而其他節(jié)點(diǎn)固定不動(dòng)的時(shí)候，桿上所承受的力?！?.2.2桿單元單元?jiǎng)偠染仃嚍椋骸?.2.2桿單元?jiǎng)偠染仃嚒?guī)推導(dǎo)方法：定義兩個(gè)形函數(shù)：位移u沿著桿的軸向線(xiàn)性分布：可得應(yīng)變?yōu)椋築為單元的應(yīng)變-位移矩陣：§3.2.2桿單元?jiǎng)偠染仃嚒?guī)推導(dǎo)方法：§3.2.2桿單元單元應(yīng)力為：桿單元上的應(yīng)變能為：桿單元2個(gè)節(jié)點(diǎn)所做的功為：對(duì)于保守系統(tǒng)，U=W

，故有：§3.2.2桿單元單元應(yīng)力為：§3.2.2桿單元上式等價(jià)于：其中k即為單元?jiǎng)偠染仃嚿鲜龇椒ň褪钦?guī)的推導(dǎo)過(guò)程，該方法也可以用來(lái)推導(dǎo)其他類(lèi)型的單元?jiǎng)偠染仃?。單元?jiǎng)偠染仃囈部梢酝ㄟ^(guò)其他嚴(yán)格的方法獲得，比如最小勢(shì)能原理，伽遼金方法等。于是我們可以得到桿單元的剛度矩陣：上式結(jié)果和直接法的結(jié)果一樣。§3.2.2桿單元上式等價(jià)于：§3.2.2桿單元例子：?jiǎn)栴}描述：在節(jié)點(diǎn)2處施加F，求桿1和桿2上的應(yīng)力。求解方法：用1-D桿單元。單元1和單元2的剛度矩陣分別為：假設(shè)在節(jié)點(diǎn)2處是用一個(gè)無(wú)摩擦的鉸鏈將桿1和桿2所連接?！?.2.2桿單元例子：§3.2.2桿單元組裝后：載荷和邊界條件為：刪除第1行第1列和第3行第3列后得：§3.2.2桿單元組裝后：§3.2.2桿單元最后可以得到單元1內(nèi)的應(yīng)力：同理，可以得到單元2內(nèi)的應(yīng)力：負(fù)號(hào)表示單元2所受的是壓應(yīng)力?！?.2.2桿單元最后可以得到單元1內(nèi)的應(yīng)力：§3.2.2桿單元注意：1)在這個(gè)例子中，根據(jù)一維線(xiàn)性理論所計(jì)算出來(lái)的單元1和2中的應(yīng)力是精確解，所以如果我們將單元細(xì)化不會(huì)提高精度。2)如果對(duì)于階梯桿結(jié)構(gòu)，其中的A要采用橫截面的平均面積。3)為了得到單元1和2中的應(yīng)力，我們首先要得到節(jié)點(diǎn)的位移，因?yàn)椴捎玫幕谖灰茍?chǎng)的有限元法。§3.2.2桿單元注意：§3.2.2桿單元二維空間上的桿單元：注意：在線(xiàn)彈性理論的前提下，橫向位移對(duì)桿單元的拉伸沒(méi)有貢獻(xiàn)。§3.2.2桿單元二維空間上的桿單元：§3.2.2桿單元坐標(biāo)變換：矩陣形式：§3.2.2桿單元坐標(biāo)變換：§3.2.2桿單元對(duì)于2個(gè)節(jié)點(diǎn)的桿單元有：節(jié)點(diǎn)力也采用相同的變換方法：二維空間下的局部坐標(biāo)系下的剛度矩陣：§3.2.2桿單元對(duì)于2個(gè)節(jié)點(diǎn)的桿單元有：§3.2.2桿單元擴(kuò)大后：矩陣形式：K為整體坐標(biāo)系下的剛度矩陣§3.2.2桿單元擴(kuò)大后：§3.2.2桿單元K的顯式表達(dá)式為：?jiǎn)卧獞?yīng)力：§3.2.2桿單元K的顯式表達(dá)式為：§3.2.2桿單元三維空間上的桿單元：先在局部坐標(biāo)系下求出單元?jiǎng)偠染仃?，然后再轉(zhuǎn)換并組裝到整體坐標(biāo)系下進(jìn)行計(jì)算。§3.2.2桿單元三維空間上的桿單元：§3.2.2桿單元簡(jiǎn)單的平面梁?jiǎn)卧夯镜牧豪碚摚骸?.2.3梁?jiǎn)卧?jiǎn)單的平面梁?jiǎn)卧骸?.2.3梁?jiǎn)卧?.3二維問(wèn)題§3.3.1

基本理論回顧§3.3.2二維問(wèn)題的有限元§3.3二維問(wèn)題§3.3.1基本理論回顧應(yīng)力和應(yīng)變：在特定的條件下，應(yīng)力應(yīng)變狀態(tài)可以簡(jiǎn)化，因此，一般的三維問(wèn)題可以簡(jiǎn)化成二維問(wèn)題來(lái)分析§3.3.1基本理論回顧應(yīng)力和應(yīng)變：§3.3.1基本理論回顧平面應(yīng)力：平面應(yīng)變：§3.3.1基本理論回顧平面應(yīng)力：§3.3.1基本理論回顧本構(gòu)關(guān)系（平面應(yīng)力）：本構(gòu)關(guān)系（平面應(yīng)變）：§3.3.1基本理論回顧本構(gòu)關(guān)系（平面應(yīng)力）：§3.3.1基本理論回顧應(yīng)變與位移的關(guān)系：邊界條件：§3.3.1基本理論回顧應(yīng)變與位移的關(guān)系：§3.3.1基本理論回顧常應(yīng)變單元（CSTorT3）位移函數(shù)：§3.3.2二維問(wèn)題的有限元常應(yīng)變單元（CSTorT3）§3.3.2二維問(wèn)題的有限位移插值函數(shù)：§3.3.2二維問(wèn)題的有限元位移插值函數(shù)：§3.3.2二維問(wèn)題的有限元應(yīng)變：CST單元?jiǎng)偠染仃嚍椋簍為單元厚度，k為一個(gè)6×6的對(duì)稱(chēng)矩陣?！?.3.2二維問(wèn)題的有限元應(yīng)變：§3.3.2二維問(wèn)題的有限元雙線(xiàn)性四邊形單元（Q4）位移函數(shù)：§3.3.2二維問(wèn)題的有限元雙線(xiàn)性四邊形單元（Q4）§3.3.2二維問(wèn)題的有限元例子：§3.3.2二維問(wèn)題的有限元例子：§3.3.2二維問(wèn)題的有限元網(wǎng)格劃分：§3.3.2二維問(wèn)題的有限元網(wǎng)格劃分：§3.3.2二維問(wèn)題的有限元計(jì)算結(jié)果：§3.3.2二維問(wèn)題的有限元計(jì)算結(jié)果：§3.3.2二維問(wèn)題的有限元載荷的變換：集中力(點(diǎn)載荷)，面力(壓力載荷)，體力(重力)為三種主要的外載荷。面力和體力都需要經(jīng)過(guò)變換以后才能加到有限元中。變換的基本思想就是：等效功?！?.3.2二維問(wèn)題的有限元載荷的變換：§3.3.2二維問(wèn)題的有限元其中，t為單元厚度，L為單元邊長(zhǎng)，un為垂直于邊界AB的位移分量。對(duì)于Q4單元來(lái)說(shuō)，有：§3.3.2二維問(wèn)題的有限元§3.3.2二維問(wèn)題的有限元應(yīng)力計(jì)算：?jiǎn)卧獌?nèi)的應(yīng)力計(jì)算公式為：B為應(yīng)變與節(jié)點(diǎn)位移之間的關(guān)系矩陣，d為計(jì)算后的節(jié)點(diǎn)位移向量。我們可以計(jì)算出單元內(nèi)部任意一點(diǎn)的應(yīng)力，包括單元中心和節(jié)點(diǎn)位置。應(yīng)力分布的等值線(xiàn)可以通過(guò)后處理軟件顯示出來(lái)?！?.3.2二維問(wèn)題的有限元應(yīng)力計(jì)算：§3.3.2二維問(wèn)題的有限元討論：1）需要了解每個(gè)單元類(lèi)型的特性：T3和Q4：線(xiàn)性位移，常量應(yīng)變和應(yīng)力T6和Q8：二次位移，線(xiàn)性的應(yīng)變和應(yīng)力2）對(duì)于一個(gè)特定的問(wèn)題選擇一種合適的單元毫無(wú)疑問(wèn)，應(yīng)該盡可能地采用高階單元和細(xì)網(wǎng)格。3）要避免單元具有大形狀比和尖角：§3.3.2二維問(wèn)題的有限元討論：§3.3.2二維問(wèn)題的有限元討論：4）單元之間要互相連通：在有限元模型中不能存在間隙或自由單元。§3.3.2二維問(wèn)題的有限元討論：§3.3.2二維問(wèn)題的有限元§3.3有限元求解技術(shù)§3.3.1

方程組求解§3.3.2有限元方法的實(shí)質(zhì)§3.3.3數(shù)值誤差§3.3.4有限元求解的收斂性§3.3有限元求解技術(shù)§3.3.1方程組求解直接法（高斯消去法）：求解時(shí)間是和NB2成正比（N為矩陣的維數(shù)，B為矩陣的帶寬）適合小到中等規(guī)模的問(wèn)題，小帶寬的問(wèn)題。對(duì)于多重載荷的情況，容易處理。迭代法：求解時(shí)間事先不可預(yù)知降低對(duì)存儲(chǔ)空間的要求適合大型規(guī)模問(wèn)題的求解（帶寬可以很大，收斂快）對(duì)于不同載荷的工況需要重新求解?！?.3.1方程組求解直接法（高斯消去法）：§3.3.1方程組求解高斯消去法例子：回代求解：§3.3.1方程組求解高斯消去法例子：§3.3.1方程組求解迭代法（高斯-賽德?tīng)枺├樱簭囊粋€(gè)初始解X(0)開(kāi)始按照以下公式計(jì)算：一直到X滿(mǎn)足以下條件，迭代終止。其中ε為收斂控制的容忍誤差?！?.3.1方程組求解迭代法（高斯-賽德?tīng)枺├樱骸?.3.1方程組求解有限元模型—是基于很多近似以后對(duì)實(shí)際結(jié)構(gòu)的一個(gè)數(shù)學(xué)模型。實(shí)際結(jié)構(gòu)—具有無(wú)限個(gè)節(jié)點(diǎn)，所以具有無(wú)限個(gè)自由度。有限元模型—具有有限個(gè)節(jié)點(diǎn)，所以具有有限個(gè)自由度。位移場(chǎng)是被有限個(gè)節(jié)點(diǎn)位移的值所控制的：§3.3.2有限元方法的實(shí)質(zhì)有限元模型—是基于很多近似以后對(duì)實(shí)際結(jié)構(gòu)的一個(gè)數(shù)學(xué)模型?！?剛度效果：有限元模型比實(shí)際結(jié)構(gòu)要更剛一些。一般來(lái)說(shuō)，位移解要比實(shí)際情況偏小一些。所以，有限元計(jì)算的位移是精確解的一個(gè)下限：§3.3.2有限元方法的實(shí)質(zhì)剛度效果：§3.3.2有限元方法的實(shí)質(zhì)誤差≠錯(cuò)誤（模型和計(jì)算）。誤差的類(lèi)型：模型誤差（梁，板，…，理論）離散誤差（有限，分段，…

）數(shù)值誤差（在求解有限元方程組的時(shí)候）例子（數(shù)值誤差）：§3.3.3數(shù)值誤差誤差≠錯(cuò)誤（模型和計(jì)算）。§3.3.3數(shù)值誤差有限元方程組為：如果K2<<K1的話(huà)，方程組是奇異的，病態(tài)的：§3.3.3數(shù)值誤差有限元方程組為：§3.3.3數(shù)值誤差如果K2>>K1的話(huà)，方程組是非病態(tài)的：在有限元模型中不同部分的剛度相差很大的話(huà)有可能導(dǎo)致有限元方程組的病態(tài)，產(chǎn)生很大的誤差。病態(tài)方程組的情況下，很小的輸入變化量（右端項(xiàng)向量）求解時(shí)會(huì)引起很大的變化量。§3.3.3數(shù)值誤差如果K2>>K1的話(huà)，方程組是非病態(tài)的：§3.3.3數(shù)隨著有限元模型中的網(wǎng)格不斷地被加密和細(xì)化，有限元的計(jì)算結(jié)果將收斂于被求解的問(wèn)題的數(shù)學(xué)模型的精確解。有限元網(wǎng)格自適應(yīng)類(lèi)型：h-refinement:減小單元尺寸（h指單元的尺寸）。p-refinement:增加單元多項(xiàng)式插值函數(shù)的階次（p指更高階次的多項(xiàng)式插值函數(shù)）。r-refinement:重新劃分有限元網(wǎng)格。hp-refinement:h方法和p方法互相結(jié)合的方法。§3.3.4有限元求解的收斂性隨著有限元模型中的網(wǎng)格不斷地被加密和細(xì)化，有限元的計(jì)算結(jié)果將網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解的收斂性網(wǎng)格自適應(yīng)類(lèi)型例子：§3.3.4有限元求解