現(xiàn)代設(shè)計理論與方法ANSYS部分

1、2021-12-15現(xiàn)代設(shè)計理論與方法現(xiàn)代設(shè)計理論與方法II. ANSYS軟件及其應(yīng)用軟件及其應(yīng)用2021-12-152 第第1章章 有限元基本理論有限元基本理論 第第2章章 ANSYS功能簡介功能簡介 第第3章章 ANSYS基本過程基本過程 第第4章章 ANSYS入門與準(zhǔn)備入門與準(zhǔn)備 第第5章章 模型輸入及修復(fù)模型輸入及修復(fù) 第第6章章 坐標(biāo)系坐標(biāo)系 第第7章章 選擇、組件與部件選擇、組件與部件 第第8章章 實體建模技術(shù)實體建模技術(shù) 第第9章章 布爾操作布爾操作 第第10章章 單元屬性單元屬性 第第11章章 網(wǎng)格劃分網(wǎng)格劃分 第第12章章 加載求解技術(shù)加載求解技術(shù) 第第13章章 后處理技術(shù)后

2、處理技術(shù) 第第14章章 結(jié)構(gòu)非線性分析結(jié)構(gòu)非線性分析 第第15章章 模態(tài)分析模態(tài)分析 第第16章章 耦合和約束方程耦合和約束方程 第第17章章 APDL基礎(chǔ)基礎(chǔ) 第第18章章 子模型子模型 第第19章章 熱分析熱分析 第第20章章 熱熱-應(yīng)力耦合分析應(yīng)力耦合分析第一章第一章 有限元基本理論有限元基本理論平衡方程幾何方程物理方程邊界條件物理系統(tǒng)有限元離散單元的位移場(假定單元內(nèi)位移函數(shù))單元節(jié)點關(guān)系求解區(qū)域的位移場、應(yīng)力場簡單化1.1 有限元分析有限元分析 (FEA)有限元分析 是利用數(shù)學(xué)近似的方法對真實物理系統(tǒng)（幾何和載荷工況）進(jìn)行模擬。它利用簡單而又相互作用的元素，即單元，用有限數(shù)量的未知量

3、去逼近無限未知量的真實系統(tǒng)。1.2 有限單元法的基本思想有限單元法的基本思想 將連續(xù)的結(jié)構(gòu)離散成有限個單元，并在每一單元中將連續(xù)的結(jié)構(gòu)離散成有限個單元，并在每一單元中設(shè)定有限個節(jié)點，將連續(xù)體看作只在節(jié)點處相連接設(shè)定有限個節(jié)點，將連續(xù)體看作只在節(jié)點處相連接的一組單元的集合體。的一組單元的集合體。 選定場函數(shù)的節(jié)點值作為基本未知量，并在每一單選定場函數(shù)的節(jié)點值作為基本未知量，并在每一單元中假設(shè)一近似插值函數(shù)，以表示單元中場函數(shù)的元中假設(shè)一近似插值函數(shù)，以表示單元中場函數(shù)的分布規(guī)律。分布規(guī)律。 利用力學(xué)中的某種變分原理去建立用以求節(jié)點未知利用力學(xué)中的某種變分原理去建立用以求節(jié)點未知量的有限單元法方程

4、，將一個連續(xù)域中有限自由度量的有限單元法方程，將一個連續(xù)域中有限自由度問題化為離散域中有限自由度問題。問題化為離散域中有限自由度問題。1.3 物理系統(tǒng)舉例物理系統(tǒng)舉例幾何體 載荷 物理系統(tǒng)結(jié)構(gòu)熱電磁000ZzyxYzyxXzyxzyzxzzyyxyzxyxx1.3.1 平衡方程平衡方程1.3.2 幾何方程幾何方程xwzuzwzvywyvyuxvxuzxzyzyxyx1.3.3 物理方程物理方程(本構(gòu)方程本構(gòu)方程)zxzxzzyzyzyyxyxyxxGGeGGeGGe222211E12EGzyxe拉梅系數(shù)體積應(yīng)變剪切模量1.3.4 邊界條件邊界條件ZnmlYnmlXnmlzyzxzzyyxyzx

5、yxxwwvvuu應(yīng)力邊界條件位移邊界條件1.4 有限元模型有限元模型真實系統(tǒng)有限元模型 有限元模型 是真實系統(tǒng)理想化的數(shù)學(xué)抽象。1.5 自由度自由度(DOFs)自由度(DOFs) 用于描述一個物理場的響應(yīng)特性。結(jié)構(gòu) DOFs 結(jié)構(gòu) 位移 熱 溫度 電 電位 流體 壓力 磁 磁位 問題 自由度ROTZUYROTYUXROTXUZ1.6 節(jié)點和單元節(jié)點和單元節(jié)點:空間中的坐標(biāo)位置，具有一定自由度和存在相互物理作用。單元: 一組節(jié)點自由度間相互作用的數(shù)值、矩陣描述（稱為剛度或系數(shù)矩陣)。單元有線、面或?qū)嶓w以及二維或三維的單元等種類。有限元模型由一些簡單形狀的單元組成，單元之間通過節(jié)點連接，并承受一

6、定載荷。載荷載荷1.6 節(jié)點和單元節(jié)點和單元 (續(xù)續(xù))信息是通過單元之間的公共節(jié)點傳遞的。分離但節(jié)點重疊的單元A和B之間沒有信息傳遞（需進(jìn)行節(jié)點合并處理）具有公共節(jié)點的單元之間存在信息傳遞 .AB.AB.1 node2 nodes1.6 節(jié)點和單元節(jié)點和單元 (續(xù)續(xù))節(jié)點自由度是隨連接該節(jié)點 單元類型 變化的。JIIJJKLILKIPOMNKJIL三維桿單元 (鉸接) UX, UY, UZ三維梁單元UX, UY, UZ, ROTX, ROTY, ROTZ二維或軸對稱實體單元UX, UY三維四邊形殼單元UX, UY, UZ, ROTX, ROTY, ROTZ三維實體熱單元TEMPJPOMNKJI

7、L三維實體結(jié)構(gòu)單元UX, UY, UZ1.7 單元形函數(shù)單元形函數(shù)FEA僅僅求解節(jié)點處的DOF值。單元形函數(shù)是一種數(shù)學(xué)函數(shù)，規(guī)定了從節(jié)點DOF值到單元內(nèi)所有點處DOF值的計算方法。因此，單元形函數(shù)提供出一種描述單元內(nèi)部結(jié)果的“形狀”。單元形函數(shù)描述的是給定單元的一種假定的特性。單元形函數(shù)與真實工作特性吻合好壞程度直接影響求解精度。真實的二次曲線.節(jié)點單元 二次曲線的線性近似 (不理想結(jié)果).21.7 單元形函數(shù)單元形函數(shù)(續(xù)續(xù))節(jié)點單元 DOF值二次分布.1節(jié)點 單元 線性近似(更理想的結(jié)果)真實的二次曲線.3節(jié)點單元二次近似 (接近于真實的二次近似擬合) (最理想結(jié)果).41.7 單元形函數(shù)

8、單元形函數(shù)(續(xù)續(xù)) DOF值可以精確或不太精確地等于在節(jié)點處的真實值可以精確或不太精確地等于在節(jié)點處的真實解，但單元內(nèi)的平均值與實際情況吻合得很好。解，但單元內(nèi)的平均值與實際情況吻合得很好。 這些平均意義上的典型解是從單元這些平均意義上的典型解是從單元DOFs推導(dǎo)出來推導(dǎo)出來的（如：結(jié)構(gòu)應(yīng)力、熱梯度）。的（如：結(jié)構(gòu)應(yīng)力、熱梯度）。 nodal solutionelement solutionUX, UY, UZ, ROTX, ROTY, ROTZ、E1.7 單元形函數(shù)單元形函數(shù)(續(xù)續(xù)) 如果單元形函數(shù)不能精確描述單元內(nèi)部的如果單元形函數(shù)不能精確描述單元內(nèi)部的DOFs，就不能很好地得到導(dǎo)出數(shù)據(jù)，

9、因為這些導(dǎo)出數(shù)據(jù)是就不能很好地得到導(dǎo)出數(shù)據(jù)，因為這些導(dǎo)出數(shù)據(jù)是通過單元形函數(shù)推導(dǎo)出來的。通過單元形函數(shù)推導(dǎo)出來的。 當(dāng)選擇了某種單元類型時，也就十分確定地選擇并當(dāng)選擇了某種單元類型時，也就十分確定地選擇并接受接受該種單元類型所假定的單元形函數(shù)。該種單元類型所假定的單元形函數(shù)。 在選定單元類型并隨之確定了形函數(shù)的情況下，必在選定單元類型并隨之確定了形函數(shù)的情況下，必須確保分析時有須確保分析時有足夠足夠數(shù)量的單元和節(jié)點來精確描述數(shù)量的單元和節(jié)點來精確描述所要求解的問題。所要求解的問題。1.8 直桿受自重作用的拉伸問題直桿受自重作用的拉伸問題)(11111iiexxxxxxxxeiluuiuuNNu

10、xxuuiiiiiiiii1.8 直桿受自重作用的拉伸問題直桿受自重作用的拉伸問題(續(xù)續(xù)) 就整個直桿來說，位移函數(shù)就整個直桿來說，位移函數(shù)U(x)是未知的，但對每是未知的，但對每一單元可以近似地假設(shè)一位移函數(shù)，它在結(jié)點上等一單元可以近似地假設(shè)一位移函數(shù)，它在結(jié)點上等于結(jié)點位移。此處，假設(shè)單元中的位移按線性分布于結(jié)點位移。此處，假設(shè)單元中的位移按線性分布 ，即：，即：1.8 直桿受自重作用的拉伸問題直桿受自重作用的拉伸問題(續(xù)續(xù)) 有了位移插值函數(shù)，就可以按材料力學(xué)公式求出應(yīng)有了位移插值函數(shù)，就可以按材料力學(xué)公式求出應(yīng)變和應(yīng)力用節(jié)點位移表示的公式：變和應(yīng)力用節(jié)點位移表示的公式：iiiiiilu

11、uEixixluudxduixE)(112)()()(1111iiiiiiiillqluuEAluuEA 1.8 直桿受自重作用的拉伸問題直桿受自重作用的拉伸問題(續(xù)續(xù)) 外載荷與結(jié)點的平衡方程外載荷與結(jié)點的平衡方程為第i個結(jié)點上承受的外載荷2)(1iillqTqaaEAaEAaEAuuupkuuqauuuqauu)()2()2(432243432321.8 直桿受自重作用的拉伸問題直桿受自重作用的拉伸問題(續(xù)續(xù)) 假定將直桿分割成假定將直桿分割成3個單元，每個單元長為個單元，每個單元長為a=L/3，則對結(jié)點則對結(jié)點2，3，4列出的平衡方程為：列出的平衡方程為：1101210122/aEATk

12、qaqaqap1.8 直桿受自重作用的拉伸問題直桿受自重作用的拉伸問題(續(xù)續(xù))EAqaqEAqaEAqauuu222242832521.8 直桿受自重作用的拉伸問題直桿受自重作用的拉伸問題(續(xù)續(xù)) 聯(lián)立求解線性代數(shù)方程組得：聯(lián)立求解線性代數(shù)方程組得：1.9 有限單元法解題的一般步驟有限單元法解題的一般步驟 結(jié)構(gòu)的離散化結(jié)構(gòu)的離散化 選擇位移模式選擇位移模式 建立平衡方程建立平衡方程 求解節(jié)點位移求解節(jié)點位移 計算單元中的應(yīng)力和應(yīng)變計算單元中的應(yīng)力和應(yīng)變1.9.1 結(jié)構(gòu)的離散化結(jié)構(gòu)的離散化 將分析的結(jié)構(gòu)物分割成有限個單元體，使相鄰的將分析的結(jié)構(gòu)物分割成有限個單元體，使相鄰的單元體單元體僅僅在節(jié)點

13、處相連接，而以如此單元的結(jié)合在節(jié)點處相連接，而以如此單元的結(jié)合體去代替原來的結(jié)構(gòu)。體去代替原來的結(jié)構(gòu)。1.9.2 選擇位移模式選擇位移模式(形函數(shù)形函數(shù)) 首先對單元假設(shè)一個位移差值函數(shù)，或稱之為位移首先對單元假設(shè)一個位移差值函數(shù)，或稱之為位移模式，得到用節(jié)點位移表示單元體內(nèi)任一點的唯一模式，得到用節(jié)點位移表示單元體內(nèi)任一點的唯一的關(guān)系式的關(guān)系式 有了位移模式，就可利用幾何關(guān)系和應(yīng)力有了位移模式，就可利用幾何關(guān)系和應(yīng)力-應(yīng)變關(guān)系應(yīng)變關(guān)系表出用單元節(jié)點位移表示單元中應(yīng)變和應(yīng)力的表達(dá)表出用單元節(jié)點位移表示單元中應(yīng)變和應(yīng)力的表達(dá)式式eNueeBDB1.9.3 三角形單元的形函數(shù)三角形單元的形函數(shù) 基

14、本假定：基本假定：假定單元內(nèi)假定單元內(nèi)的位移可以用一個比較的位移可以用一個比較簡單的函數(shù)來表示，如簡單的函數(shù)來表示，如線性插值函數(shù)。這在單線性插值函數(shù)。這在單元劃分比較密的情況下元劃分比較密的情況下是合理可行的。是合理可行的。yaxaau321yaxaav6541.9.3 三角形單元的形函數(shù)三角形單元的形函數(shù)(續(xù)續(xù)) 將三角形單元的將三角形單元的3個頂點的個頂點的2個方向位移代個方向位移代入位移函數(shù)可求出入位移函數(shù)可求出6個待定系數(shù)。即可用節(jié)個待定系數(shù)。即可用節(jié)點的位移表示內(nèi)部任意一點的位移：點的位移表示內(nèi)部任意一點的位移： emmjjiimjimjiBvuvuvuNNNNNNvuu00000

15、01.9.4 建立平衡方程建立平衡方程 可利用最小勢能原理建立結(jié)構(gòu)的節(jié)點載荷和節(jié)點位可利用最小勢能原理建立結(jié)構(gòu)的節(jié)點載荷和節(jié)點位移之間的關(guān)系式，即結(jié)構(gòu)的平衡方程移之間的關(guān)系式，即結(jié)構(gòu)的平衡方程pk1.9.5 求解結(jié)點位移v將邊界條件代入線性代數(shù)方程組 后，經(jīng)解算可求得所有未知的結(jié)點位移。pk1.9.6 計算單元中的應(yīng)變和應(yīng)力計算單元中的應(yīng)變和應(yīng)力 依據(jù)求得的結(jié)點位移，由依據(jù)求得的結(jié)點位移，由可求得單元中任一點的應(yīng)變和應(yīng)力?？汕蟮脝卧腥我稽c的應(yīng)變和應(yīng)力。eeBDB 結(jié)構(gòu)的離散化結(jié)構(gòu)的離散化 用有限元法對結(jié)構(gòu)進(jìn)行應(yīng)力分析時，首先要將結(jié)構(gòu)進(jìn)用有限元法對結(jié)構(gòu)進(jìn)行應(yīng)力分析時，首先要將結(jié)構(gòu)進(jìn)行離散化。即

16、將一個連續(xù)體看成由有限個單元組成的行離散化。即將一個連續(xù)體看成由有限個單元組成的體系。彈性力學(xué)平面問題中最常見的單元是三角形單體系。彈性力學(xué)平面問題中最常見的單元是三角形單元。元。 所有作用在單元上的載荷都按靜力等效的原則移置到所有作用在單元上的載荷都按靜力等效的原則移置到結(jié)點上，并在受幾何約束的結(jié)點處設(shè)置相應(yīng)的鉸支座結(jié)點上，并在受幾何約束的結(jié)點處設(shè)置相應(yīng)的鉸支座。這樣就得到了用以代替原來彈性體的有限單元計算。這樣就得到了用以代替原來彈性體的有限單元計算模型。模型。 位移模式位移模式 取一個典型的三角形單元進(jìn)行力學(xué)分析。在有限單元位移法取一個典型的三角形單元進(jìn)行力學(xué)分析。在有限單元位移法中，假

17、設(shè)結(jié)點上的位移是基本未知量。為了能用單元的結(jié)點中，假設(shè)結(jié)點上的位移是基本未知量。為了能用單元的結(jié)點位移表示單元中的應(yīng)變和應(yīng)力分量，必須假定一個位移模式位移表示單元中的應(yīng)變和應(yīng)力分量，必須假定一個位移模式，也就是說根據(jù)單元的結(jié)點位移去構(gòu)造單元上的位移插值函，也就是說根據(jù)單元的結(jié)點位移去構(gòu)造單元上的位移插值函數(shù)。數(shù)。位移模式（續(xù)）位移模式（續(xù)） 位移插值函數(shù)位移插值函數(shù) 采用線性插值，即假定單元上的位移分量是坐標(biāo)的線采用線性插值，即假定單元上的位移分量是坐標(biāo)的線性函數(shù)：性函數(shù)： 它們可以由結(jié)點位移確定如下：它們可以由結(jié)點位移確定如下：位移模式(續(xù))聯(lián)立求解上述方程，可得：聯(lián)立求解上述方程，可得：位

18、移模式(續(xù))其中：其中：而：而：是三角形是三角形ijm的面積。的面積。位移模式(續(xù))于是可以得到：于是可以得到：其中：其中：同理得：同理得：位移模式(續(xù))可以將位移模式改寫為矩陣模式：可以將位移模式改寫為矩陣模式： 單元中的應(yīng)變和應(yīng)力 有了單元的位移模式，就可以借助平面問題的幾何和有了單元的位移模式，就可以借助平面問題的幾何和物理方程，導(dǎo)出用單于的結(jié)點位移表示單元中的應(yīng)變物理方程，導(dǎo)出用單于的結(jié)點位移表示單元中的應(yīng)變和應(yīng)力分量的公式。和應(yīng)力分量的公式。 由：由：單元中的應(yīng)變和應(yīng)力(續(xù))得到：得到：或簡寫為：或簡寫為：單元中的應(yīng)變和應(yīng)力(續(xù))將應(yīng)變代入物理方程：將應(yīng)變代入物理方程：可得：可得：

19、即為用單元中的結(jié)點位移表示單元中應(yīng)力的關(guān)系式。即為用單元中的結(jié)點位移表示單元中應(yīng)力的關(guān)系式。單元中的應(yīng)變和應(yīng)力(續(xù)) 式中式中D為彈性矩陣，對于平面應(yīng)力問題，矩陣為：為彈性矩陣，對于平面應(yīng)力問題，矩陣為： 單元的總勢能單元的總勢能 我們已經(jīng)知道由各個單元的位移模式就形成了整個結(jié)我們已經(jīng)知道由各個單元的位移模式就形成了整個結(jié)構(gòu)的位移模式。按彈性力學(xué)最小勢能原理，結(jié)構(gòu)中最構(gòu)的位移模式。按彈性力學(xué)最小勢能原理，結(jié)構(gòu)中最接近于真實解的位移應(yīng)該是使結(jié)構(gòu)總勢能取得最小值接近于真實解的位移應(yīng)該是使結(jié)構(gòu)總勢能取得最小值的那組位移函數(shù)。的那組位移函數(shù)。 由于在位移函數(shù)公式中，結(jié)點位移為自變量，這樣就由于在位移函

20、數(shù)公式中，結(jié)點位移為自變量，這樣就使一個泛函的極值問題變?yōu)橐粋€多元函數(shù)的極值問題使一個泛函的極值問題變?yōu)橐粋€多元函數(shù)的極值問題。為此我們來討論單元的總勢能關(guān)于結(jié)點位移的表達(dá)。為此我們來討論單元的總勢能關(guān)于結(jié)點位移的表達(dá)式。式。 每一個單元的總勢能由該單元的應(yīng)變能以及此單元上每一個單元的總勢能由該單元的應(yīng)變能以及此單元上所有外力的勢能組成。所有外力的勢能組成。 單元的應(yīng)變能單元的應(yīng)變能 平面應(yīng)力狀態(tài)下，設(shè)物體厚度為平面應(yīng)力狀態(tài)下，設(shè)物體厚度為h，則單元中的應(yīng)變能，則單元中的應(yīng)變能為：為：單元的應(yīng)變能（續(xù)）將將和和Bi代入上式，應(yīng)用矩陣相乘的轉(zhuǎn)置的逆序法則代入上式，應(yīng)用矩陣相乘的轉(zhuǎn)置的逆序法則,注

21、意到彈性矩陣注意到彈性矩陣D的對稱性的對稱性,有：有：單元的應(yīng)變能（續(xù)）因為矩陣因為矩陣B及及D的元素都是常量，所以可記：的元素都是常量，所以可記：單元的應(yīng)變能（續(xù)）從而單元的應(yīng)變能可寫為：從而單元的應(yīng)變能可寫為：利用利用=Be，有：，有：單元的應(yīng)變能（續(xù)）注意到注意到B=Bi Bj Bm，記子矩陣，記子矩陣 單元上體積力的勢能單元上體積力的勢能物體中常見的體力為旋轉(zhuǎn)離心體力和重力。在平面問題中物體中常見的體力為旋轉(zhuǎn)離心體力和重力。在平面問題中，體積力在，體積力在z軸方向的分力為零，設(shè)單元體積中的體積軸方向的分力為零，設(shè)單元體積中的體積力為：力為： 單元上體積力具有的勢能為：單元上體積力具有的

22、勢能為： 單元上表面力的勢能單元上表面力的勢能設(shè)物體邊界上一單元某邊上受到表面力的作用，單位長度上設(shè)物體邊界上一單元某邊上受到表面力的作用，單位長度上所受到的表面力為：所受到的表面力為：則單元上表面力的勢能為：則單元上表面力的勢能為： 單元節(jié)點上集中力的勢能單元節(jié)點上集中力的勢能如果彈性物體受到集中力如果彈性物體受到集中力Re 的作用，通常劃分單元網(wǎng)格的作用，通常劃分單元網(wǎng)格時都在集中力的作用點設(shè)置結(jié)點。設(shè)某單元時都在集中力的作用點設(shè)置結(jié)點。設(shè)某單元3個結(jié)點上所個結(jié)點上所受到的集中力為：受到的集中力為：于是該單元上集中力的勢能是：于是該單元上集中力的勢能是： 單元中的總勢能單元中的總勢能 綜合

23、前面的幾種情況，可以得到單元中的總勢能為：綜合前面的幾種情況，可以得到單元中的總勢能為： 單元中的總勢能單元中的總勢能 分別引進(jìn)單元體積力，表面力，集中力向量如下：分別引進(jìn)單元體積力，表面力，集中力向量如下： 單元中的總勢能單元中的總勢能 則單元中的總勢能可以表示為：則單元中的總勢能可以表示為： 物體中的總勢能物體中的總勢能 把各單元的總勢能疊加起來，就可得到整個彈性體的總勢把各單元的總勢能疊加起來，就可得到整個彈性體的總勢能。為了便于疊加和歸并，需將單元剛度矩陣表達(dá)式能。為了便于疊加和歸并，需將單元剛度矩陣表達(dá)式(2-18)作適當(dāng)?shù)母膶?。作適當(dāng)?shù)母膶憽?假設(shè)結(jié)構(gòu)離散化后共有假設(shè)結(jié)構(gòu)離散化后共

24、有n個結(jié)點，將編號為個結(jié)點，將編號為 l的結(jié)點位移記的結(jié)點位移記為：為： 則結(jié)構(gòu)的結(jié)點位移向量：則結(jié)構(gòu)的結(jié)點位移向量： 是一個是一個2n維的列向量。維的列向量。物體中的總勢能（續(xù)）物體中的總勢能（續(xù)） 可將單元剛度矩陣式用補(bǔ)零的辦法由可將單元剛度矩陣式用補(bǔ)零的辦法由6X6的矩陣擴(kuò)大到的矩陣擴(kuò)大到2nX2n的矩陣的矩陣物體中的總勢能（續(xù)）物體中的總勢能（續(xù)）如果在物體上劃分的單元總數(shù)是如果在物體上劃分的單元總數(shù)是e0,再引進(jìn)再引進(jìn)結(jié)構(gòu)的總剛度陣：結(jié)構(gòu)的總剛度陣：物體總勢能就可寫為：物體總勢能就可寫為：物體中的總勢能（續(xù)）物體中的總勢能（續(xù)） 代入約束條件后的彈性體總勢能可以寫為：代入約束條件后的

25、彈性體總勢能可以寫為： 空間問題的有限單元法空間問題的有限單元法 用有限單元法求解彈性力學(xué)空間問題，首先也要將連續(xù)的空用有限單元法求解彈性力學(xué)空間問題，首先也要將連續(xù)的空間物體用一系列的單元離散化。間物體用一系列的單元離散化。 空間問題中，最簡單的是四面體單元。離散的空間結(jié)構(gòu)是這空間問題中，最簡單的是四面體單元。離散的空間結(jié)構(gòu)是這些單元只在節(jié)點處以空間鉸相互連接的集合體。些單元只在節(jié)點處以空間鉸相互連接的集合體?？臻g問題的有限單元法（續(xù)）空間問題的有限單元法（續(xù)）位移模式位移模式 空間問題中，每一個結(jié)點有空間問題中，每一個結(jié)點有3個位移分量，單元結(jié)點個位移分量，單元結(jié)點位移向量由位移向量由12

26、個分量組成，分別表示為：個分量組成，分別表示為：位移模式（續(xù)）位移模式（續(xù)）假定單元內(nèi)的位移分量為坐標(biāo)的線性函數(shù)假定單元內(nèi)的位移分量為坐標(biāo)的線性函數(shù):位移模式（續(xù)）位移模式（續(xù)）將上式中的第一式應(yīng)用于將上式中的第一式應(yīng)用于4個結(jié)點，則有：個結(jié)點，則有：位移模式（續(xù)）位移模式（續(xù)） 由上式可解出由上式可解出a1,a2,a3和和a4再代回位移分量的表達(dá)式再代回位移分量的表達(dá)式，可得：，可得： 式中：式中： 為形函數(shù)，其中：為形函數(shù)，其中：位移模式（續(xù)）位移模式（續(xù)）位移模式（續(xù)）位移模式（續(xù)）用同樣的方法，可以得到：用同樣的方法，可以得到：合并合并,的表達(dá)式，可以將單元內(nèi)任一點的位移寫為：的表達(dá)式

27、，可以將單元內(nèi)任一點的位移寫為：單元中的應(yīng)變和應(yīng)力單元中的應(yīng)變和應(yīng)力 在空間問題中，每點有在空間問題中，每點有6個應(yīng)變分量，由幾何關(guān)系：個應(yīng)變分量，由幾何關(guān)系：將,的表達(dá)式代入上式，得到:式中：單元中的應(yīng)變和應(yīng)力單元中的應(yīng)變和應(yīng)力(續(xù)續(xù))單元中的應(yīng)變和應(yīng)力單元中的應(yīng)變和應(yīng)力(續(xù)續(xù)) 可以看出，應(yīng)變矩陣可以看出，應(yīng)變矩陣B中的元素都是常量，從而單中的元素都是常量，從而單元中的應(yīng)變都是常量，故線性位移模式的四面體單元元中的應(yīng)變都是常量，故線性位移模式的四面體單元是常應(yīng)變單元。是常應(yīng)變單元。 由應(yīng)力由應(yīng)力-應(yīng)變關(guān)系，得到單元中的應(yīng)力為：應(yīng)變關(guān)系，得到單元中的應(yīng)力為： 式中式中D為一般空間問題的彈性矩

28、陣為一般空間問題的彈性矩陣 從下面從下面D的表達(dá)式可以看出，單元中的應(yīng)力都是常的表達(dá)式可以看出，單元中的應(yīng)力都是常數(shù)。數(shù)。單元中的應(yīng)變和應(yīng)力單元中的應(yīng)變和應(yīng)力(續(xù)續(xù)) 單元剛度矩陣和結(jié)點載單元剛度矩陣和結(jié)點載荷向量荷向量 仿照平面問題的推導(dǎo)，可以得到四面體單元的剛度矩陣：仿照平面問題的推導(dǎo)，可以得到四面體單元的剛度矩陣： 分塊形式：分塊形式：單元剛度矩陣和結(jié)點載荷單元剛度矩陣和結(jié)點載荷向量（續(xù)）向量（續(xù)）式中子矩陣可以表達(dá)為：式中子矩陣可以表達(dá)為：其中：其中：單元剛度矩陣和結(jié)點載荷單元剛度矩陣和結(jié)點載荷向量（續(xù)）向量（續(xù)）經(jīng)過與平面問題中同樣的推導(dǎo)，單元的體積力向量和表經(jīng)過與平面問題中同樣的推

29、導(dǎo)，單元的體積力向量和表面力向量可以用下列公式計算：面力向量可以用下列公式計算：經(jīng)疊加，組合，得有限元支配方程：經(jīng)疊加，組合，得有限元支配方程：代入約束條件，可解出結(jié)點位移向量，從而就可以求出代入約束條件，可解出結(jié)點位移向量，從而就可以求出各單元的應(yīng)變和應(yīng)力。各單元的應(yīng)變和應(yīng)力。第二章第二章 ANSYS軟件的功能簡介軟件的功能簡介 前處理模塊分析計算模塊后處理模塊三大模塊結(jié)構(gòu)分析：靜動力、非線性熱分析（滲流分析）流體動力學(xué)分析電磁場分析聲場分析壓電分析多物理場耦合分析通用后處理模塊時間歷程后處理模塊實體建模網(wǎng)格劃分 ANSYS在部分工業(yè)領(lǐng)域中的應(yīng)用如下：在部分工業(yè)領(lǐng)域中的應(yīng)用如下： 航空航天航

30、空航天 汽車工業(yè)汽車工業(yè) 生物醫(yī)學(xué)生物醫(yī)學(xué) 橋梁、建筑橋梁、建筑 電子產(chǎn)品電子產(chǎn)品 重型機(jī)械重型機(jī)械 微機(jī)電系統(tǒng)微機(jī)電系統(tǒng) 運動器械運動器械 ANSYS/Multiphysics 包括所有工程學(xué)科的所有性能包括所有工程學(xué)科的所有性能 ANSYS/Multiphysics有三個主要的組成產(chǎn)品有三個主要的組成產(chǎn)品 ANSYS/Mechanical - ANSYS/機(jī)械機(jī)械-結(jié)構(gòu)及熱結(jié)構(gòu)及熱 ANSYS/Emag -ANSYS電磁學(xué)電磁學(xué) ANSYS/FLOTRAN - ANSYS計算流體動力學(xué)計算流體動力學(xué) 其它產(chǎn)品其它產(chǎn)品: ANSYS/LS-DYNA -高度非線性結(jié)構(gòu)問題高度非線性結(jié)構(gòu)問題 D

31、esignSpace CAD環(huán)境下，適合快速分析容易使環(huán)境下，適合快速分析容易使用的設(shè)計和分析工具用的設(shè)計和分析工具 ANSYS/ProFEA Pro/ENGINEER的的ANSYS 分分析接口。析接口。2.1 前處理模塊前處理模塊PREP7 實體建模方式之一實體建模方式之一 ：自頂向下：自頂向下 先建高級圖元，如圓柱、圓錐等自動生成相關(guān)的面、線及關(guān)鍵點自頂向下問題：如何保證各實體的連接、交叉等關(guān)系？布爾操作2.1 前處理模塊前處理模塊PREP7(續(xù)續(xù)) 實體建模方式之二實體建模方式之二 ：自底向上：自底向上 關(guān)鍵點自底向上線面體2.1 前處理模塊前處理模塊PREP7(續(xù)續(xù)) 實體建模方式的選

32、擇：實體建模方式的選擇：自頂向下建?？梢蕴岣呓５男?，但在需要利用自頂向下建模可以提高建模的效率，但在需要利用布爾操作時比較難以掌握；布爾操作時比較難以掌握；自底向上建?？梢詼p少出錯的機(jī)會，但效率較低。自底向上建模可以減少出錯的機(jī)會，但效率較低。建議：先學(xué)習(xí)自底向上建議：先學(xué)習(xí)自底向上后學(xué)習(xí)自頂向下后學(xué)習(xí)自頂向下 網(wǎng)格劃分方法：延伸劃分與映射劃分網(wǎng)格劃分方法：延伸劃分與映射劃分延伸劃分將一個二維網(wǎng)格延伸成一個三維網(wǎng)格映射劃分將幾何模型分解成簡單的幾部分，然后選擇合適的單元屬性和網(wǎng)格控制，生成映射網(wǎng)格。2.1 前處理模塊前處理模塊PREP7(續(xù)續(xù)) 網(wǎng)格劃分方法：自由劃分和自適應(yīng)劃分網(wǎng)格劃分方

33、法：自由劃分和自適應(yīng)劃分自由劃分可對復(fù)雜模型直接劃分，避免了用戶對各個部分分別劃分然后進(jìn)行組裝時各部分網(wǎng)格不匹配帶來的麻煩。自適應(yīng)劃分先生成具有邊界條件的實體模型，用戶指示程序自動地生成有限元網(wǎng)格，分析、估計網(wǎng)格的離散誤差，然后重新定義網(wǎng)格大小，再次分析計算、估計網(wǎng)格的離散誤差，直至誤差低于用戶定義的值或達(dá)到用戶定義的求解次數(shù)。2.1 前處理模塊前處理模塊PREP7(續(xù)續(xù))2.2 求解模塊求解模塊SOLUTION結(jié)構(gòu)靜力分析用來求解外載荷引起的位移、應(yīng)力和反作用力。靜力分析很適合求解慣性和阻尼對結(jié)構(gòu)的影響并不顯著的問題。ANSYS程序中的靜力分析不僅可以進(jìn)行線性分析，而且可以進(jìn)行非線性分析，如

34、塑性、蠕變、膨脹、大變形、大應(yīng)變、應(yīng)力剛化及接觸分析。用來求解隨時間變化的載荷對結(jié)構(gòu)或部件的影響。與靜力分析不同，動力分析要考慮隨時間變化的力載荷以及它對阻尼和慣性的影響。結(jié)構(gòu)動力分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù)) 結(jié)構(gòu)動力分析的類型：結(jié)構(gòu)動力分析的類型：模態(tài)分析模態(tài)分析：計算線性結(jié)構(gòu)的自振頻率及振型。：計算線性結(jié)構(gòu)的自振頻率及振型。譜分析譜分析：是模態(tài)分析的擴(kuò)展，用于計算由于隨機(jī)：是模態(tài)分析的擴(kuò)展，用于計算由于隨機(jī)振動引起的結(jié)構(gòu)應(yīng)力和應(yīng)變振動引起的結(jié)構(gòu)應(yīng)力和應(yīng)變(也叫作響應(yīng)譜或也叫作響應(yīng)譜或PSD)。諧響應(yīng)分析諧響應(yīng)分析：確定

35、線性結(jié)構(gòu)對隨時間按正弦曲線：確定線性結(jié)構(gòu)對隨時間按正弦曲線變化的載荷的響應(yīng)。變化的載荷的響應(yīng)。瞬態(tài)動力學(xué)分析瞬態(tài)動力學(xué)分析：確定結(jié)構(gòu)：確定結(jié)構(gòu) 對隨時間任意變化的載荷的響應(yīng)，對隨時間任意變化的載荷的響應(yīng)， 可以考慮與靜力分析相同的結(jié)構(gòu)可以考慮與靜力分析相同的結(jié)構(gòu) 非線性行為。非線性行為。結(jié)構(gòu)非線性導(dǎo)致結(jié)構(gòu)或部件的響應(yīng)隨外載荷不成比例變化。ANSYS程序可以求解靜態(tài)和瞬態(tài)非線性問題，包括材料非線性、幾何非線性和單元非線性三種。結(jié)構(gòu)非線性分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))ANSYS程序可以分析大型柔體運動。當(dāng)運動的積累影響起主要作用時，可使用這些功能分析復(fù)雜結(jié)構(gòu)在空間中的運動特性

36、，并確定結(jié)構(gòu)中由此產(chǎn)生的應(yīng)力、應(yīng)變和變形。動力學(xué)分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))用于計算線性屈曲載荷并確定屈曲模態(tài)形狀 (結(jié)合瞬態(tài)動力學(xué)分析可以實現(xiàn)非線性屈曲分析) 。特征屈曲分析斷裂分析復(fù)合材料分析疲勞分析專項分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))ANSYS除了提供標(biāo)準(zhǔn)的隱式動力學(xué)分析以外，還提供了顯式動力學(xué)分析模塊ANSYS/LS-DYNA。顯式動力學(xué)分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù)) 顯式動力學(xué)分析的特點：顯式動力學(xué)分析的特點：用于模擬非常大的變形，用于模擬非常大的變形，慣性力占支配地位，

37、并考慣性力占支配地位，并考慮所有的非線性行為。慮所有的非線性行為。它的顯式方程求解沖擊、它的顯式方程求解沖擊、碰撞、快速成型等問題，碰撞、快速成型等問題，是目前求解這類問題最有是目前求解這類問題最有效的方法。效的方法。 熱分析之后往往進(jìn)行結(jié)構(gòu)分熱分析之后往往進(jìn)行結(jié)構(gòu)分析，計算由于熱膨脹或收縮析，計算由于熱膨脹或收縮不均勻引起的應(yīng)力。不均勻引起的應(yīng)力。 ANSYS熱分析功能熱分析功能:相變相變(熔化及凝固熔化及凝固)內(nèi)熱源內(nèi)熱源(如電阻發(fā)熱等如電阻發(fā)熱等)三種熱傳遞方式三種熱傳遞方式(熱傳導(dǎo)、熱熱傳導(dǎo)、熱對流、熱輻射對流、熱輻射)ANSYS熱分析計算物體的穩(wěn)態(tài)或瞬態(tài)溫度分布，以及熱量的獲取或損失

38、、熱梯度、熱通量等。熱分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))磁場分析中考慮的物理量是磁通量密度、磁場密度、磁力、磁力矩、阻抗、電感、渦流、能耗及磁通量泄漏等。磁場可由電流、永磁體、外加磁場等產(chǎn)生。磁場分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))磁場分析的類型磁場分析的類型:靜磁場分析靜磁場分析：計算直流電：計算直流電(DC)或永磁體產(chǎn)生的磁場。或永磁體產(chǎn)生的磁場。交變磁場分析交變磁場分析：計算由于交流電：計算由于交流電(AC)產(chǎn)生的磁場。產(chǎn)生的磁場。瞬態(tài)磁場分析瞬態(tài)磁場分析：計算隨時間隨機(jī)：計算隨時間隨機(jī)變化的電流或外界引起

39、的磁場。變化的電流或外界引起的磁場。用于計算電阻或電容系統(tǒng)的電場。 典型的物理量有電流密度、電荷密度、電場及電阻熱等。電場分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))用于微波及RF無源組件，波導(dǎo)、雷達(dá)系統(tǒng)、同軸連接器等分析。高頻電磁場分析同軸電纜中的電場 (EFSUM) 2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))用于確定流體的流動及熱行為流體分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù)) 流體分析的分類：流體分析的分類：CFD -ANSYS/FLOTRAN 提供強(qiáng)大的計算流體提供強(qiáng)大的計算流體動力學(xué)分析功能，包括不可壓縮或可壓縮流體、動力學(xué)分析功能，包括不可壓縮或可壓縮流體

40、、層流及湍流，以及多組份流等。應(yīng)用于：航空航層流及湍流，以及多組份流等。應(yīng)用于：航空航天，電子元件封裝，汽車設(shè)計。天，電子元件封裝，汽車設(shè)計。聲學(xué)分析聲學(xué)分析 - 考慮流體介質(zhì)與周圍固體的相互作用考慮流體介質(zhì)與周圍固體的相互作用, 進(jìn)行聲波傳遞或水下結(jié)構(gòu)的動力學(xué)分析等。例如進(jìn)行聲波傳遞或水下結(jié)構(gòu)的動力學(xué)分析等。例如：揚聲器、汽車內(nèi)部、聲納。典型的物理量是：壓：揚聲器、汽車內(nèi)部、聲納。典型的物理量是：壓力分布、位移和自振頻率。力分布、位移和自振頻率。容器內(nèi)流體容器內(nèi)流體 分析分析 - 考慮容器內(nèi)的非流動流體的考慮容器內(nèi)的非流動流體的影響?？梢源_定由于晃動引起的靜水壓力。例如：影響。可以確定由于晃

41、動引起的靜水壓力。例如：油罐，其它液體容器。油罐，其它液體容器。流體動力學(xué)耦合流體動力學(xué)耦合分析分析 - 在考慮流體約束質(zhì)量的動在考慮流體約束質(zhì)量的動力響應(yīng)基礎(chǔ)上，力響應(yīng)基礎(chǔ)上，在結(jié)構(gòu)動力學(xué)分析中使用流體耦合在結(jié)構(gòu)動力學(xué)分析中使用流體耦合單元。單元。2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))考慮兩個或多個物理場之間的相互作用。如果兩個物理場之間相互影響，單獨求解一個物理場是不可能得到正確結(jié)果的，因此你需要一個能夠?qū)蓚€物理場組合到一起求解的分析軟件。例如：在壓電力分析中，需要同時求解電壓分布（電場分析）和應(yīng)變（結(jié)構(gòu)分析）。耦合場分析2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))v典型耦

42、合場分析：熱應(yīng)力分析流體結(jié)構(gòu)相互作用感應(yīng)加熱（電磁熱）壓電分析（電場和結(jié)構(gòu)）聲學(xué)分析（流體和結(jié)構(gòu)）熱-電分析靜電-結(jié)構(gòu)分析兩根熱膨脹系數(shù)不同的棒焊接在一起，圖示為加熱后的變形。2.2 求解模塊求解模塊SOLUTION(續(xù)續(xù))2.3 后處理模塊后處理模塊POST1和和POST26這個模塊對前面的分析結(jié)果能以圖形形式顯示和輸出。例如，計算結(jié)果（如應(yīng)力）在模型上的變化情況可用等值線圖表示，不同的等值線顏色代表了不同的值（如應(yīng)力值） 。POST1這個模塊用于檢查在一個時間段或子步歷程中的結(jié)果，如節(jié)點位移、應(yīng)力或支反力。這些結(jié)果能通過繪制曲線或列表查看。繪制一個或多個變量隨頻率或其他量變化的曲線，有助于

43、形象化地表示分析結(jié)果。POST26第三章第三章 ANSYS的基本過程的基本過程 一個典型的一個典型的ANSYS分析過程可分為以下分析過程可分為以下3個步驟：個步驟：前處理求解后處理A1XYZ3.1 前處理前處理前處理指定工程名稱和分析標(biāo)題定義單位定義單元類型定義單元常數(shù)創(chuàng)建橫截面定義材料特性創(chuàng)建有限元模型定義分析類型求解控制加載3.1.1 指定工程名稱和分析標(biāo)題指定工程名稱和分析標(biāo)題更改工程名定義分析標(biāo)題3.1.2 定義單位定義單位 除電磁分析以外，你不必為除電磁分析以外，你不必為ANSYS設(shè)置單位系設(shè)置單位系統(tǒng)。簡單地確定你將采用的單位制，然后保證所統(tǒng)。簡單地確定你將采用的單位制，然后保證所

44、有輸入數(shù)據(jù)均采用該種單位制就可以。有輸入數(shù)據(jù)均采用該種單位制就可以。 你確定的單位制將影響尺寸、實常數(shù)、材料特性你確定的單位制將影響尺寸、實常數(shù)、材料特性和載荷等的輸入值。和載荷等的輸入值。 ANSYS 不進(jìn)行單位換算不進(jìn)行單位換算! 它只是簡單地接受所輸它只是簡單地接受所輸入的數(shù)據(jù)而不會懷疑它們的合法性。入的數(shù)據(jù)而不會懷疑它們的合法性。 命令命令/UNITS允許你指定單位制允許你指定單位制, 但它只是一個紀(jì)但它只是一個紀(jì)錄設(shè)計，從而讓使用你模型的用戶知道你所用的錄設(shè)計，從而讓使用你模型的用戶知道你所用的單位。單位。3.1.2 定義單位定義單位(續(xù)續(xù)) 使用使用/UNITS命令可以設(shè)置系統(tǒng)單位

45、，沒有相應(yīng)的命令可以設(shè)置系統(tǒng)單位，沒有相應(yīng)的GUI。 USER：用戶自定義單位，是缺省設(shè)置：用戶自定義單位，是缺省設(shè)置 SI：國際單位制，：國際單位制，m, kg, s, BFT：以英尺為基礎(chǔ)的單位制，：以英尺為基礎(chǔ)的單位制，ft, slug, s, F CGS：cm, g, s, MPA：mm, mg, s, BIN：以英寸為基礎(chǔ)的單位制：以英寸為基礎(chǔ)的單位制in, lbm, s, F3.1.3 定義單元類型定義單元類型BEAMCIRCUitCOMBINationCONTACtFLUIDHF(High Frequency)HYPERelasticINFINiteLINKMASSMATRIXM

46、ESHPIPEPLANEPRETS(pretension)SHELLSOLIDSOURCeSURFaceTARGEtTRANSducerUSERVISCOelastic3.1.4 定義單元常數(shù)定義單元常數(shù) 單元實常數(shù)是由單元類型單元實常數(shù)是由單元類型的特性決定的，如梁單元的特性決定的，如梁單元的橫截面特性。并不是所的橫截面特性。并不是所有的單元類型都需要實常有的單元類型都需要實常數(shù)，同類型的不同單元也數(shù)，同類型的不同單元也可以有不同的實常數(shù)。可以有不同的實常數(shù)。指定單元的實常數(shù)號3.1.5 創(chuàng)建橫截面創(chuàng)建橫截面創(chuàng)建梁的橫截面3.1.6 定義材料特性定義材料特性定義材料特性指定單元材料號3.1.

47、7 定義分析類型求解控制定義分析類型求解控制定義分析類型求解控制基本設(shè)置瞬態(tài)設(shè)置求解選項非線性設(shè)置求解終止的高級控制 包括：自由度約束、力包括：自由度約束、力、表面分布載荷、體積、表面分布載荷、體積載荷、慣性載荷、耦合載荷、慣性載荷、耦合場載荷場載荷載荷步：僅指可求得載荷步：僅指可求得解的載荷設(shè)置。解的載荷設(shè)置。子步：是指在一個載子步：是指在一個載荷步中每次增加的步長荷步中每次增加的步長，主要是為了在瞬態(tài)分，主要是為了在瞬態(tài)分析和非線性分析中提高析和非線性分析中提高分析精度和收斂性。子分析精度和收斂性。子步也稱作時間步，代表步也稱作時間步，代表一段時間。一段時間。3.1.8 加載加載3.2 求

48、解求解求解當(dāng)前載荷步求解某載荷步3.3 通用后處理器通用后處理器畫出分析的結(jié)果用列表的形式列出分析的結(jié)果查詢某些結(jié)點或者單元處的應(yīng)力值以及其它分析選項Deformed Shape表示畫出變形后的形狀。有如下選項：3.3.1 畫出分析的結(jié)果畫出分析的結(jié)果3.3.2 畫出節(jié)點的結(jié)果畫出節(jié)點的結(jié)果位移轉(zhuǎn)角3.3.3 求解自由度結(jié)果求解自由度結(jié)果正應(yīng)力和剪應(yīng)力主應(yīng)力應(yīng)力強(qiáng)度平均等效應(yīng)力3.3.4 求解應(yīng)力結(jié)果求解應(yīng)力結(jié)果正應(yīng)變和剪應(yīng)變主應(yīng)變應(yīng)變強(qiáng)度平均等效應(yīng)變3.3.5 求解總應(yīng)變結(jié)果求解總應(yīng)變結(jié)果求解能量彈性應(yīng)變?nèi)渥兤渌鼞?yīng)變正應(yīng)變和剪應(yīng)變主應(yīng)變應(yīng)變強(qiáng)度平均等效應(yīng)變3.3.6 其它求解結(jié)果其它求解結(jié)果

49、3.3.7 圖形圖形輸出選項輸出選項只畫出變形后的圖形畫出變形前后的圖形畫出變形后的圖形和變形前的邊界圖3.4 時間歷程后處理器時間歷程后處理器適用于：瞬態(tài)分析的后處理。3.5 懸壁梁算例懸壁梁算例使用ANSYS分析一個工字懸壁梁，如圖所示。PPoint ALH求解在力P作用下點A處的變形，已知條件如下：P = 4000 lbL = 72 inI = 833 in4E = 29 E6 psi橫截面積 (A) = 28.2 in2H = 12.71 in見beam.txt1.啟動啟動 ANSYS： 以交互模式進(jìn)入以交互模式進(jìn)入ANSYS，工作文件，工作文件名為名為beam。2.創(chuàng)建基本模型創(chuàng)建基

50、本模型：(使用帶有兩個關(guān)鍵點的線模擬梁，使用帶有兩個關(guān)鍵點的線模擬梁，梁的高度及橫截面積將在單元的實常數(shù)中設(shè)置梁的高度及橫截面積將在單元的實常數(shù)中設(shè)置)a. Main Menu: Preprocessor -Modeling- Create Keypoints In Active CS.b.輸入關(guān)鍵點編號輸入關(guān)鍵點編號 1，輸入，輸入x,y,z坐標(biāo)坐標(biāo) 0,0,0，選擇，選擇 Applyc. 輸入關(guān)鍵點編號輸入關(guān)鍵點編號 2，輸入，輸入x,y,z坐標(biāo)坐標(biāo)72,0,0，選擇，選擇 OKd.Main Menu: Preprocessor -Modeling- Create -Lines- Line

51、s Straight Linee. 選取兩個關(guān)鍵點，在拾取菜單中選擇選取兩個關(guān)鍵點，在拾取菜單中選擇OK3.存儲存儲ANSYS數(shù)據(jù)庫數(shù)據(jù)庫：Toolbar: SAVE_DB4.設(shè)定分析模塊設(shè)定分析模塊：a.Main Menu: Preferencesb.選擇選擇 Structural.c. 選擇選擇 OK.5. 設(shè)定單元類型相應(yīng)選項設(shè)定單元類型相應(yīng)選項：a.Main Menu:PreprocessorElement TypeAdd/Edit/Deleteb.選擇選擇 Add . . .c. 左邊單元庫列表中選擇左邊單元庫列表中選擇 Beam.d.在右邊單元列表中選擇在右邊單元列表中選擇 2D

52、elastic (BEAM3).e. 選擇選擇 OK 接受單元類型并關(guān)閉對話框接受單元類型并關(guān)閉對話框.f. 選擇選擇 Close 關(guān)閉單元類型對話框關(guān)閉單元類型對話框.6. 定義實常數(shù)定義實常數(shù)：a.Main Menu: Preprocessor Real Constantsb.選擇選擇 Add . . .c. 選擇選擇 OK 定義定義BEAM3的實常數(shù)的實常數(shù).d.選擇選擇 Help 得到有關(guān)單元得到有關(guān)單元 BEAM3的幫助的幫助.e. 查閱單元描述查閱單元描述.f. File Exit 退出幫助系統(tǒng)退出幫助系統(tǒng).g.在在AREA框中輸入框中輸入 28.2 （橫截面積）橫截面積）.h.在

53、在IZZ框中輸入框中輸入 833 （慣性矩）慣性矩）.i. 在在HEIGHT框中輸入框中輸入 12.71 （梁的高度）梁的高度）.j. 選擇選擇 OK 定義實常數(shù)并關(guān)閉對話框定義實常數(shù)并關(guān)閉對話框.k.選擇選擇 Close 關(guān)閉實常數(shù)對話框關(guān)閉實常數(shù)對話框.7.定義材料屬性定義材料屬性：a.Preprocessor Material Props -Constant- Isotropicb.選擇選擇 OK to 定義材料定義材料 1.c. 在在EX框中輸入框中輸入29e6（彈性模量）（彈性模量）.d.選擇選擇OK 定義材料屬性并關(guān)閉對話框定義材料屬性并關(guān)閉對話框.8. 保存保存ANSYS數(shù)據(jù)庫文

54、件數(shù)據(jù)庫文件 beamgeom.db：a.Utility Menu: File Save asb.輸入文件名輸入文件名 beamgeom.db.c. 選擇選擇 OK 保存文件并退出對話框保存文件并退出對話框.9.對幾何模型劃分網(wǎng)格對幾何模型劃分網(wǎng)格：a.Main Menu: Preprocessor MeshToolb.選擇選擇 Mesh.c. 拾取拾取 line.d.在拾取對話框中選擇在拾取對話框中選擇 OK.e. (可選可選) 在在MeshTool對話框中選擇對話框中選擇 Close.10.保存保存ANSYS數(shù)據(jù)庫到文件數(shù)據(jù)庫到文件 beammesh.db：a.Utility Menu:

55、File Save asb.輸入文件名：輸入文件名： beammesh.db.c. 選擇選擇 OK 保存文件并退出對話框保存文件并退出對話框.11.施加載荷及約束施加載荷及約束：a. Main Menu: Solution -Loads- Apply -Structural- Displacement On Nodesb.拾取最左邊的節(jié)點，在拾取菜單中選擇拾取最左邊的節(jié)點，在拾取菜單中選擇 OK.c. 選擇選擇All DOF.d. 選擇選擇 OK. （如果不輸入任何值，位移約束默認(rèn)為（如果不輸入任何值，位移約束默認(rèn)為0）e. Main Menu: Solution -Loads- Apply

56、-Structural- Force/Moment On Nodesf. 拾取最右邊的節(jié)點，在選取對話框中選擇拾取最右邊的節(jié)點，在選取對話框中選擇OK.g.選擇選擇 FY.h.在在 VALUE框中輸入框中輸入 -4000，選擇，選擇 OK. 12.保存數(shù)據(jù)庫文件到保存數(shù)據(jù)庫文件到 beamload.db：a.Utility Menu: File Save asb.輸入文件名輸入文件名 beamload.db.c. 選擇選擇OK保存文件并關(guān)閉對話框保存文件并關(guān)閉對話框.13.進(jìn)行求解進(jìn)行求解：a.Main Menu: Solution -Solve- Current LSb.查看狀態(tài)窗口中的信息

57、查看狀態(tài)窗口中的信息, 然后選擇然后選擇 File Closec. 選擇選擇 OK開始計算開始計算.d.當(dāng)出現(xiàn)當(dāng)出現(xiàn) “Solution is done!” 提示后，選擇提示后，選擇OK關(guān)閉此窗口關(guān)閉此窗口.14.進(jìn)入通用后處理讀取分析結(jié)果進(jìn)入通用后處理讀取分析結(jié)果:Main Menu: General Postproc -Read Results- First Set15.圖形顯示變形圖形顯示變形：a. Main Menu: General Postproc Plot Results Deformed Shapeb.在對話中選擇在對話中選擇 deformed and undeformed.c

58、. 選擇選擇 OK.16.(可選可選) 列出反作用力列出反作用力：a. Main Menu: General Postproc List Results Reaction Solub.選擇選擇 OK 列出所有項目，并關(guān)閉對話框列出所有項目，并關(guān)閉對話框.c. 看完結(jié)果后，選擇看完結(jié)果后，選擇File Close 關(guān)閉窗口關(guān)閉窗口.17.退出退出ANSYS：a.工具條工具條: Quitb.選擇選擇Quit - No Save!c. 選擇選擇 OK.18.與彈性梁理論計算對比與彈性梁理論計算對比： ya = (PL3)/(3EI) = 0.0206 in.兩個結(jié)果一致兩個結(jié)果一致.3.6 車床刀具

59、算例車床刀具算例說明說明 下面求解一個車床刀具模型的下面求解一個車床刀具模型的 3-D應(yīng)力分析。通過下列步驟在通應(yīng)力分析。通過下列步驟在通用后處理器中查看結(jié)果：用后處理器中查看結(jié)果：1) 畫位侈畫位侈2) 列反力列反力3) 畫畫von Mises 應(yīng)力應(yīng)力4) 動態(tài)顯示動態(tài)顯示von Mises應(yīng)力應(yīng)力施加不同的約束重新求解，并與施加不同的約束重新求解，并與第一次的結(jié)果進(jìn)行比較。第一次的結(jié)果進(jìn)行比較。第一次求解載荷及材料特性第一次求解載荷及材料特性:第二次求解時的荷載及材料特性：第二次求解時的荷載及材料特性：1. 用用“cutter”作為作業(yè)名，進(jìn)入作為作業(yè)名，進(jìn)入ANSYS。2. 設(shè)置設(shè)置

60、GUI 優(yōu)先選擇優(yōu)先選擇 為為“結(jié)構(gòu)結(jié)構(gòu)”: Main Menu Preferences 選擇選擇 “Structural”, 然后按然后按OK3. 選擇選擇 “No defeaturing”，輸入，輸入 “cutter.igs” IGES 文件文件: Utility Menu File Import IGES 選擇選擇 “No defeaturing”,然后按然后按 OK 選擇選擇“cutter.igs,然后按然后按OK 或用命令或用命令:/AUX15 IOPTN,IGES,NODEFEAT IGESIN,cutter,igs, 4. 按一定比例轉(zhuǎn)換模型，把厘米單位轉(zhuǎn)換英寸：按一定比例轉(zhuǎn)換模