1、有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義分析的對象的一些行為計(jì)算出的幾何項(xiàng)求解的自由度及應(yīng)力反作用力或節(jié)點(diǎn)力有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義21.分析的對象的一些基本的行為分析的對象的一些基本的行為:重力方向總是豎直向下的重力方向總是豎直向下的離心力總是沿徑向向外的離心力總是沿徑向向外的沒有一種材料能抵抗沒有一種材料能抵抗 1,000,000 psi 的應(yīng)力的應(yīng)力軸對稱的物體幾乎沒有為零的軸對稱的物體幾乎沒有為零的 環(huán)向應(yīng)力環(huán)向應(yīng)力彎曲載荷造成的應(yīng)力使一側(cè)受壓，另一側(cè)受拉彎曲載荷造成的應(yīng)力使一側(cè)受壓，另一側(cè)受拉如果只有一個(gè)載荷施加在結(jié)構(gòu)上，檢驗(yàn)結(jié)果比較容易如果只有一個(gè)載荷施加在結(jié)構(gòu)上，

2、檢驗(yàn)結(jié)果比較容易. 如果有多個(gè)載荷，可單獨(dú)如果有多個(gè)載荷，可單獨(dú)施加一個(gè)或幾個(gè)載荷分別檢驗(yàn)，然后施加所有載荷檢驗(yàn)分析結(jié)果施加一個(gè)或幾個(gè)載荷分別檢驗(yàn)，然后施加所有載荷檢驗(yàn)分析結(jié)果.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義33.檢驗(yàn)求解的自由度及應(yīng)力：確認(rèn)施加在模型上的載荷環(huán)境是合理的.確認(rèn)模型的運(yùn)動(dòng)行為與預(yù)期的相符 - 無剛體平動(dòng)、無剛體轉(zhuǎn)動(dòng)、無裂縫等.確認(rèn)位移和應(yīng)力的分布與期望的相符，或者利用物理學(xué)或數(shù)學(xué)可以解釋.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義4模型所有的反作用力應(yīng)該與施加的點(diǎn)力、壓力和慣性力平衡.在所有約束節(jié)點(diǎn)的豎直方向的反作用力.在所有約束節(jié)點(diǎn)水平方向的反作用力必須與水平方向的載

3、荷平衡.所有約束節(jié)點(diǎn)的反作用力矩必須與施加的載荷平衡.注意包含在約束方程中自由度的反力，不包括由這個(gè)約束方程傳遞的力.必須與施加的豎直方向的載荷平衡有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義5在任意選取的單元字集中的節(jié)點(diǎn)力，應(yīng)與作用在結(jié)構(gòu)此部分的已知載荷向平衡，除非節(jié)點(diǎn)的符號(hào)約定與自由體圖上所示的相反.未選擇的單元上的豎直方向的節(jié)點(diǎn)總力.必須與被選擇的單元上施加的豎直方向的載荷平衡注意包含在約束方程中自由度的反力，不包括由這個(gè)約束方程傳遞的力.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義7ANSYS通用后處理包含網(wǎng)格離散誤差估計(jì)通用后處理包含網(wǎng)格離散誤差估計(jì).

4、誤差估計(jì)是依據(jù)沿單元內(nèi)邊界的應(yīng)力或熱流的不連續(xù)性，是平均與未平均節(jié)點(diǎn)應(yīng)力間的差值誤差估計(jì)是依據(jù)沿單元內(nèi)邊界的應(yīng)力或熱流的不連續(xù)性，是平均與未平均節(jié)點(diǎn)應(yīng)力間的差值.Elem 1Elem 2savg = 1100savg = 1100s = 1200s = 1200s = 1000s = 1000savg = 1200savg = 1200s = 1300s = 1300s = 1100s = 1100(節(jié)點(diǎn)的節(jié)點(diǎn)的 ss 是積分點(diǎn)的外插）是積分點(diǎn)的外插）有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義8誤差估計(jì)作用條件：誤差估計(jì)作用條件：線性靜力結(jié)構(gòu)分析及線性穩(wěn)態(tài)熱分析線性靜力結(jié)構(gòu)分析及線性穩(wěn)態(tài)熱分析大

5、多數(shù)大多數(shù) 2-D 或或 3-D 實(shí)體或殼單元實(shí)體或殼單元PowerGraphics off誤差信息誤差信息:能量百分比誤差能量百分比誤差 sepc單元應(yīng)力偏差單元應(yīng)力偏差 sdsg單元能量偏差單元能量偏差 serr應(yīng)力上、下限應(yīng)力上、下限 smnb smxb有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義9能量百分比誤差是對所選擇的單元的位移、應(yīng)力能量百分比誤差是對所選擇的單元的位移、應(yīng)力、溫度或熱流密度的粗略估計(jì)、溫度或熱流密度的粗略估計(jì). 它可以用于比較承它可以用于比較承受相似載荷的相似結(jié)構(gòu)的相似模型受相似載荷的相似結(jié)構(gòu)的相似模型.這個(gè)值的通常應(yīng)該在這個(gè)值的通常應(yīng)該在10%以下以下. 如果不選擇

6、其他單如果不選擇其他單元，而只選擇在節(jié)點(diǎn)上施加點(diǎn)載荷或應(yīng)力集中處元，而只選擇在節(jié)點(diǎn)上施加點(diǎn)載荷或應(yīng)力集中處的單元，誤差值有時(shí)會(huì)達(dá)到的單元，誤差值有時(shí)會(huì)達(dá)到50%或以上或以上.SEPC 2 %PowerGraphic offMain menu general postproc plot results deformed shape 選 ：Def+undefedge有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義10要檢驗(yàn)?zāi)硞€(gè)位置的網(wǎng)格離散應(yīng)力誤差，可要檢驗(yàn)?zāi)硞€(gè)位置的網(wǎng)格離散應(yīng)力誤差，可以列出或繪制應(yīng)力偏差以列出或繪制應(yīng)力偏差.某一個(gè)單元的應(yīng)力偏差是此單元上全部節(jié)某一個(gè)單元的應(yīng)力偏差是此單元上全部節(jié)點(diǎn)的六

7、個(gè)應(yīng)力分量值與此節(jié)點(diǎn)的平均應(yīng)力點(diǎn)的六個(gè)應(yīng)力分量值與此節(jié)點(diǎn)的平均應(yīng)力值之差的最大值值之差的最大值.應(yīng)力偏差：應(yīng)力偏差：節(jié)點(diǎn)節(jié)點(diǎn)n的應(yīng)力矢量：的應(yīng)力矢量：所關(guān)心位置上的應(yīng)力偏差值450 psi (30,000 psi 應(yīng)力的1.5%)察看應(yīng)力偏差：Plot Results Element Solu Error Estimation Stress deviation (SDSG) inaninsss= neNiinanNne=1ss有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l平均應(yīng)力為4421 （nodal solution）l應(yīng)力偏差為689.598l誤差=689.598/4421=15.53%（局

8、部細(xì)化）有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義12每個(gè)單元的另一種誤差值是能量誤差. 它與單元上節(jié)點(diǎn)應(yīng)力差值有關(guān)的, 用于計(jì)算選擇的單元的能量百分比誤差.察看能量誤差：Plot Results Element Solu Error Estimation Energy error (ENER). 應(yīng)力誤差矢量應(yīng)變矩陣單元的應(yīng)力單元體積的能量誤差單元其中：sss=DvolievoldDeivolTi)(211單元數(shù)整個(gè)模型的能量誤差：Nr 1=nriiee有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義13應(yīng)力上下限可以確定由于網(wǎng)格離散誤差對模型的應(yīng)力最大值的影響.顯示或列出的應(yīng)力上下限包括:估計(jì)的上限

9、 - SMXB估計(jì)的下限 - SMNB應(yīng)力上下限限并不是估計(jì)實(shí)際的最高或最小應(yīng)力。它定義了一個(gè)確信范圍。 如果沒有其他的確鑿的驗(yàn)證，就不能認(rèn)為實(shí)際的最大應(yīng)力低于 SMXB.例如:SMX=32750是節(jié)點(diǎn)解的實(shí)際值 SMXB=33200是估計(jì)的上限X stress SMAX 32,750 psiSMXB 33,200 psi(difference 450 psi 1.5 %)max()min(najmmxbjnajmmnbjssssss=有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義采用plane42單元網(wǎng)格局部細(xì)化與未細(xì)化能量百分比誤差局部細(xì)化未細(xì)化DisplacementDMX=0.88E-03S

10、EPC=14.442DMX=0.803E-03應(yīng)力偏差Element Solution(SDSG)SMN=63.453SMX=426.86SDSGSMN=64.528SMX=689.589能量誤差估計(jì)Element Solution(SERR)SMN=0.365E-03SMX=0.600595SERRSMN=0.005173SMX=0.38503應(yīng)力上下限Nodal solution(SEQV)SMN=725.21SMNB=720.133SMX=4579SMXB=4623SEQVSMN=773.769SMNB=708.94SMX=4421SMXB=4999有限元分析及應(yīng)用講義有限元分析及應(yīng)用講

11、義15P方法應(yīng)用控制：P方法用于線彈性結(jié)構(gòu)分析實(shí)體和殼體。P單元由以下5種單元：2-D Quadrilateral (Plane145)2-D Triangle (Plane146)3-D Brick (Solid 147)3-D Tetrehedron (Solid 148)3-D Shell (Solid 150)規(guī)定規(guī)定 0.1% 局部應(yīng)力差，使用局部應(yīng)力差，使用p方法計(jì)算的最方法計(jì)算的最大大X方向應(yīng)力約為方向應(yīng)力約為 34,700 psi(比普通比普通h方法高出大約方法高出大約 5% )P 單元的位移形函數(shù)u=a1+a2x+a3y+a4x2+a5xy+a6y2v=a7+a8x+a9y+

12、 a10 x2+a11xy+a12y2P方法的優(yōu)點(diǎn)：如果使用 p-方法 進(jìn)行結(jié)構(gòu)分析，可以依靠p單元自動(dòng)調(diào)整單元多項(xiàng)式階數(shù)2-8），達(dá)到收斂到設(shè)定的精度. 對這種方法的相信程度，與使用經(jīng)驗(yàn)有關(guān).有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義lE=30e6 lb/in2lV=0.29lThick=0.25inl在節(jié)點(diǎn)0，5，0處的收斂標(biāo)準(zhǔn)設(shè)為1%2010pres=-100lb/in2ininR=5 in有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義延伸網(wǎng)格劃分映射網(wǎng)格劃分層狀網(wǎng)格劃分有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l將一個(gè)二維網(wǎng)格延伸生成一個(gè)三維網(wǎng)格；三維網(wǎng)格

13、生成后去掉二維網(wǎng)格l步驟：l1.先生成橫截面l2.指定網(wǎng)格密度并對面進(jìn)行網(wǎng)格劃分l3.拖拉面網(wǎng)格生成體網(wǎng)格l 指定單元屬性l 拖拉，完成體網(wǎng)格劃分。l4.釋放已選的平面單元有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義機(jī)翼沿著長度方向輪廓一致，且它的橫截面由直線和樣條曲線定義。機(jī)翼的一端固定在機(jī)體上，另一端為懸空的自由端。采樣點(diǎn):A(0,0,0) B(2,0,0) C(2.3,0.2,0) D(1.9,0.45,0) E(1,0.25,0)xyz210斜度=0.25彈性模量Ex=38E03 psi泊松比：0.3密度：D=1.033e-3 slugs/in3有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義截面

14、寬度：10mm 截面形狀：正六變形手柄長度： 20cm 桿長 ： 7.5cm導(dǎo)角半徑： 1cm 彈性模量：2.07E11pa 延伸網(wǎng)格劃分：作業(yè)有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義有兩種主要的網(wǎng)格劃分方法: 自由劃分和映射劃分.自由劃分無單元形狀限制.網(wǎng)格無固定的模式.適用于復(fù)雜形狀的面和體.映射劃分面的單元形狀限制為四邊形,體的單元限制為六面體 (方塊).通常有規(guī)則的形式,單元明顯成行.僅適用于 “規(guī)則的” 面和體, 如 矩形和方塊.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義自由網(wǎng)格易于生成; 不須將復(fù)雜形狀的體分解為規(guī)則形狀的體.體單元僅包含四面體網(wǎng)格, 致使單元數(shù)量較多.僅高階 (10

15、-節(jié)點(diǎn)) 四面體單元較滿意, 因此DOF(自由度)數(shù)目可能很多. 映射網(wǎng)格通常包含較少的單元數(shù)量.低階單元也可能得到滿意的結(jié)果,因此DOF(自由度)數(shù)目較少.面和體必須形狀 “規(guī)則”, 劃分的網(wǎng)格必須滿足一定的準(zhǔn)則.難于實(shí)現(xiàn), 尤其是對形狀復(fù)雜的體.網(wǎng)格劃分的優(yōu)缺點(diǎn)：有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義自由網(wǎng)格自由網(wǎng)格自由網(wǎng)格是面和體網(wǎng)格劃分時(shí)的缺省設(shè)自由網(wǎng)格是面和體網(wǎng)格劃分時(shí)的缺省設(shè)置置.生成自由網(wǎng)格比較容易生成自由網(wǎng)格比較容易:導(dǎo)出導(dǎo)出 MeshTool 工具工具, 劃分方式設(shè)為自由劃分方式設(shè)為自由劃分劃分.推薦使用智能網(wǎng)格劃分推薦使用智能網(wǎng)格劃分 進(jìn)行自由網(wǎng)格劃進(jìn)行自由網(wǎng)格劃分分,

16、激活它并指定一個(gè)尺寸級(jí)別激活它并指定一個(gè)尺寸級(jí)別. 存儲(chǔ)存儲(chǔ)數(shù)據(jù)庫數(shù)據(jù)庫.按按 Mesh 按鈕開始劃分網(wǎng)格按鈕開始劃分網(wǎng)格.按拾取器中按拾取器中 Pick All 選擇所有實(shí)體選擇所有實(shí)體 (推推薦薦).或使用命令或使用命令 VMESH,ALL 或或 AMESH,ALL.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義映射網(wǎng)格劃分映射網(wǎng)格劃分由于面和體必須滿足一定的要求由于面和體必須滿足一定的要求,生成映射網(wǎng)格不如生成自生成映射網(wǎng)格不如生成自由網(wǎng)格容易由網(wǎng)格容易 :面必須包含面必須包含 3 或或 4 條線條線 (三角形或四邊形三角形或四邊形).體必須包含體必須包含4, 5, 或或 6 個(gè)面?zhèn)€面 (四面

17、體四面體, 三棱柱三棱柱, 或六面體或六面體).對邊的單元分割必須匹配對邊的單元分割必須匹配.對三角形面或四面體對三角形面或四面體, 單元分割數(shù)必須為偶數(shù)單元分割數(shù)必須為偶數(shù).有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l因而 ,映射網(wǎng)格劃分包含以下三個(gè)步驟:l保證 “規(guī)則的形狀, 即, 面有 3 或4 條邊, 或 體有 4, 5, 或 6 個(gè)面.l指定尺寸和形狀控制l生成網(wǎng)格0有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義1.保證規(guī)則的形狀保證規(guī)則的形狀在許多情況下在許多情況下, 模型的幾何形狀上有多于模型的幾何形狀上有多于4條邊的面條邊的面,有多有多于于6個(gè)面的體個(gè)面的體. 為了將它們轉(zhuǎn)換成規(guī)則的形

18、狀為了將它們轉(zhuǎn)換成規(guī)則的形狀, 您可能您可能進(jìn)行如下的一項(xiàng)或兩項(xiàng)操作進(jìn)行如下的一項(xiàng)或兩項(xiàng)操作:把面把面 (或體或體) 切割成小的切割成小的, 簡單的形狀簡單的形狀.連接兩條或多條線連接兩條或多條線 (或面或面) 以減少總的邊數(shù)以減少總的邊數(shù).有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l切割 （divide可以通過布爾減運(yùn)算實(shí)現(xiàn).l您可以使用工作平面, 一個(gè)面, 或一條線 作為切割工具.l有時(shí), 生成一條新的線或面會(huì)比移動(dòng)或定向工作平面到正確的方向容易得多.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l連接 操作是生成一條新線 (為網(wǎng)格劃分) , 它通過連接兩條或多條線以減少構(gòu)成面的線數(shù).l使用 LC

19、CAT 命令或 Preprocessor -Meshing- Concatenate Lines, 然后拾取須連接的線.l對面進(jìn)行連接, 使用 ACCAT 命令或Preprocessor -Meshing- Concatenate Areasl若兩條線或兩個(gè)面 相切交匯可考慮用加 (布爾) 運(yùn)算連接這兩條線使其成為一個(gè)由4條邊構(gòu)成的面Concatenate 有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義L1 and L2 are added.New Line L#L4 and L5 are concatenated.產(chǎn)生四條線圍成的面，適于網(wǎng)格劃分產(chǎn)生四條線圍成的面，適于網(wǎng)格劃分New concat

20、enated line.原始邊的單元數(shù)為原始邊的單元數(shù)為4條條連接邊的單元數(shù)為連接邊的單元數(shù)為8條條有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l您也可以簡單地通過一個(gè)面上的3個(gè)或4個(gè)角點(diǎn) 暗示 一個(gè)連接. 此時(shí), ANSYS 內(nèi)在地 生成一個(gè)連接.l在MeshTool中選擇Quad shape 和 Map 網(wǎng)格.l將 3/4 sided 變?yōu)?Pick corners.l按 Mesh 鍵, 拾取面, 然后拾取 3 或 4 角點(diǎn)形成一規(guī)則的形狀.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l使用連接時(shí)注意:l它僅僅是一個(gè)網(wǎng)格劃分操作,因而應(yīng)為網(wǎng)格劃分前的最后一步, 在所有的實(shí)體建模之后. 這是因?yàn)?，?jīng)

21、連接操作得到的實(shí)體不能在后續(xù)的實(shí)體建模操作中使用.l可以通過刪除產(chǎn)生的線或面 “undo(取消)” 一個(gè)連接.l連接面 (為在體上映射網(wǎng)格) 通常比較復(fù)雜,因?yàn)槟矐?yīng)該連接一些線. 只有在對相鄰的兩個(gè)4邊形面作連接時(shí)其中的線會(huì)自動(dòng)連接.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義指定尺寸和形狀控制這是映射網(wǎng)格劃分3個(gè)步驟中的第2步.選擇單元形狀非常簡單. 在 MeshTool中,對面的網(wǎng)格劃分選擇 Quad,對體的網(wǎng)格劃分選擇 Hex, 點(diǎn)擊 Map.其中通常采用的尺寸控制和級(jí)別如下:線尺寸 LESIZE 級(jí)別較高.若指定了總體單元尺寸, 它將用于 “未給定尺寸的” 線.缺省的單元尺寸 DESIZ

22、E僅在未指定ESIZE時(shí)用于 “未給定尺寸的” 線上.(智能網(wǎng)格劃分 無效.)有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義l若您指定線的分割數(shù), 切記:l對邊的分割數(shù)必須匹配, 但您只須指定一邊的分割數(shù). 映射網(wǎng)格劃分器 將把分割數(shù)自動(dòng)傳送到它的對邊.l如果模型中有連接(Concatenate)線, 只能在原始(輸入)線上指定分割數(shù),而不能在合成線上指定分割數(shù).每條初始線上指定6份分割.此線上將自動(dòng)使用12 份分割 (合成線的對邊).其它兩條線上會(huì)采用幾 份分割 呢? (后面的演示將會(huì)回答這一問題.)有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義生成映射網(wǎng)格只要保證了規(guī)則的形狀 并指定了合適的份數(shù), 生

23、成網(wǎng)格將非常簡單. 只須按MeshTool中的 Mesh 鍵, 然后按拾取器中的 Pick All 或選擇需要的實(shí)體即可.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義說明用自由及映射網(wǎng)格對輪模型進(jìn)行混合的網(wǎng)格劃分.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義1. 按指定的工作目錄，以 “wheelb-3d為作業(yè)名，進(jìn)入ANSYS.或清除 ANSYS 數(shù)據(jù)庫，改換作業(yè)名為 “wheelb-3d”:Utility Menu File Clear & Start New .Utility Menu File Change Jobname .2. 恢復(fù) “wheelb.db1” 數(shù)據(jù)庫文件:Utili

24、ty Menu File Resume from 選擇 “wheelb.db” 數(shù)據(jù)庫文件, 然后選擇 OK或用命令:RESUME,wheelb,db13. 進(jìn)入前處理器，用工作平面切分體 :Main Menu Preprocessor -Modeling- Operate -Booleans- Divide Volu by WrkPlane +拾取Pick AllUtility Menu Plot Volumes或用命令:/PREP7VSBW,1VPLOT有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義4. 平移工作平面到19號(hào)關(guān)鍵點(diǎn):Utility Menu WorkPlane Offset WP

25、 to Keypoints +選擇如圖所示的19號(hào)關(guān)鍵點(diǎn), 然后選擇 OK或用命令:KWPAVE,19有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義5. 以工作平面切分體:Main Menu Preprocessor -Modeling- Operate -Booleans- Divide Volu by WrkPlane +拾取Pick AllUtility Menu Plot Volumes或用命令:VSBW,4VPLOT有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義6. 關(guān)閉工作平面，設(shè)置總體單元尺寸為 0.25:Utility Menu WorkPlane Display Working Plan

26、eMain Menu Preprocessor MeshTool 設(shè)置大小控制為Global，按 Set設(shè)置SIZE = 0.25按OK或用命令:WPSTYLEESIZE,0.257. 用SOLID45單元，對四個(gè)外部的體進(jìn)行映射網(wǎng)格劃分 (TYPE 1):Main Menu Preprocessor MeshTool 在Shape下選擇 “Hex” 和 “Mapped” :按Mesh拾取四個(gè)外部的體 (體號(hào) 1, 2, 3, 和 5)按OK或用命令:MSHAPE,0,3DMSHKEY,1VMESH,1,3,1VMESH,5有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義有限元分析

27、及應(yīng)用講義8. 用SOLID95單元，對內(nèi)部的體進(jìn)行自由網(wǎng)格劃分 (TYPE 2):Main Menu Preprocessor MeshTool 單擊單元屬性Element Attributes下的 Set :TYPE = “2 SOLID95”, 然后選擇 OK設(shè)置大小控制為Global，按 Set設(shè)置SIZE = 0.2按OK在Shape 下選擇 “Tet和“Free” :按Mesh拾取內(nèi)部的體 (體號(hào) 6)按OK或用命令:TYPE,2ESIZE,0.2MSHAPE,1,3DMSHKEY,0VMESH,6有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義9. 將 SOLID95單元轉(zhuǎn)變?yōu)?SOLI

28、D92單元:Main Menu Preprocessor -Meshing - Modify Mesh Change Tets .按OK或用命令:TCHG,95,9210.選擇并畫出 SOLID95 四面體單元:Utility Menu Select Entities .選擇 “Elements”, “By Attributes”, “Elem type num”設(shè)置 Min,Max,Inc = 2按OKUtility Menu Plot Elements或用命令:ESEL,S,TYPE,2EPLOT有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義11.選擇 “全部實(shí)體并保存數(shù)據(jù)庫 :Utility

29、Menu Select EverythingUtility Menu Plot ElementsUtility Menu File Save as 輸入數(shù)據(jù)庫文件名“wheelb-3d-mesh.db”, 然后選擇 OK或用命令:ALLSEL,ALLEPLOTSAVE,wheelb-3d-mesh,db11.退出ANSYS:在工具條中選擇 “QUIT”選擇 “Quit - No Save!”按OK或用命令:FINISH/EXIT,NOSAVE有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義適用于2D情況，生成線性過渡的自由網(wǎng)格平行于邊線方向的單元尺寸相當(dāng)垂直于邊線方向的單元尺寸和數(shù)目急劇變化當(dāng)分析要求

30、邊界單元高精度時(shí)，層狀網(wǎng)格很有用有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義GUI：Main Menu: Preprocessor MeshTool Layer Set button指定：線上的單元尺寸，線上兩端單元的比率和內(nèi)部網(wǎng)格層的厚度。線間距比率space）,對層狀劃分一般取1.0內(nèi)部網(wǎng)格層厚度layer1線上單元尺寸系數(shù)： size factor=2 沿線生成兩行尺寸均勻的單元外部網(wǎng)格層厚度layer2這層的單元尺寸會(huì)從layer1緩慢增加到總體單元尺寸，layer2的厚度可以用一個(gè)網(wǎng)格過渡系數(shù) 如：Transition factor=2 生成大約等于前面垂直于線網(wǎng)格2倍尺寸的單元有限元分

31、析及應(yīng)用講義有限元分析及應(yīng)用講義47在完成此章的學(xué)習(xí)之后，給出一個(gè)已經(jīng)劃分好網(wǎng)格的模型的數(shù)據(jù)庫文件，我們應(yīng)該能夠使用耦合或約束方程來建立節(jié)點(diǎn)自由度之間的聯(lián)系第一講第一講 耦合耦合定義耦合設(shè)置定義耦合設(shè)置說明耦合的三種普遍應(yīng)用說明耦合的三種普遍應(yīng)用.采用采用3種不同的方法建立耦合關(guān)系種不同的方法建立耦合關(guān)系.第二講第二講 約束方程約束方程定義定義“約束方程約束方程”說明約束方程的四種普遍應(yīng)用說明約束方程的四種普遍應(yīng)用采用四種不同的方法生成約束方程采用四種不同的方法生成約束方程.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義48l耦合是使一組節(jié)點(diǎn)具有相同的自由度值.l除了自由度值是由求解器計(jì)算而非用戶指

32、定外，與約束相類似。l例如:如果節(jié)點(diǎn)1和節(jié)點(diǎn)2在UX方向上耦合， 求解器將計(jì)算節(jié)點(diǎn)1的UX值并簡單地把該值賦值給節(jié)點(diǎn)2的UX。l一個(gè)耦合設(shè)置是一組被約束在一起，有著同一方向的節(jié)點(diǎn) (即一個(gè)自由度)。l一個(gè)模型中可以定義多個(gè)耦合，但一個(gè)耦合中只能包含一個(gè)方向的自由度。有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義491.施加對稱性條件：耦合自由度常被用來實(shí)施移動(dòng)或循環(huán)對稱條件.考慮在均勻軸向壓力下的空心長圓柱體，此3D結(jié)構(gòu)可用下面右圖所示的2D軸對稱模型表示.xyxy123451112131415由于結(jié)構(gòu)的對稱性，上面的一排結(jié)點(diǎn)在軸向由于結(jié)構(gòu)的對稱性，上面的一排結(jié)點(diǎn)在軸向上的位移應(yīng)該相同上的位移應(yīng)該相

33、同有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義502. 無摩擦的界面 如果滿足下列條件，則可用耦合自由度來模擬接觸面:表面保持接觸, 此分析是幾何線性的小變形）忽略摩擦在兩個(gè)界面上，節(jié)點(diǎn)是一一對應(yīng)的.通過僅耦合垂直于接觸面的移動(dòng)來模擬接觸. 優(yōu)點(diǎn):分析仍然是線性的無間隙收斂性問題有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義513.鉸接耦合可用來模擬力耦松弛，例如鉸鏈、無摩擦滑動(dòng)器、萬向節(jié)考慮一個(gè)考慮一個(gè)2D的梁模型，每個(gè)節(jié)點(diǎn)上有三個(gè)的梁模型，每個(gè)節(jié)點(diǎn)上有三個(gè)自由度自由度ux、uy和和rotz，A點(diǎn)為一鉸鏈連接點(diǎn)為一鉸鏈連接。將同一位置節(jié)點(diǎn)的自由度。將同一位置節(jié)點(diǎn)的自由度ux、uy耦合起耦合起來。來。1

34、2A節(jié)點(diǎn)節(jié)點(diǎn)1和節(jié)點(diǎn)和節(jié)點(diǎn)2處于同一位置處于同一位置，但為于清楚，但為于清楚起見，在圖上起見，在圖上分開顯示。分開顯示。.為了模擬鉸接，將同一位置兩個(gè)節(jié)點(diǎn)為了模擬鉸接，將同一位置兩個(gè)節(jié)點(diǎn)的移動(dòng)自由度耦合起來，而不耦合轉(zhuǎn)的移動(dòng)自由度耦合起來，而不耦合轉(zhuǎn)動(dòng)自由度動(dòng)自由度有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義52在循環(huán)對稱切面上的對應(yīng)位置實(shí)在循環(huán)對稱切面上的對應(yīng)位置實(shí)施自由度耦合。施自由度耦合。用耦合施加循環(huán)對稱性用耦合施加循環(huán)對稱性有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義53進(jìn)入創(chuàng)建耦合關(guān)系的菜單路徑進(jìn)入創(chuàng)建耦合關(guān)系的菜單路徑:Main Menu: Preprocessor Coupling /

35、 Ceqn Couple DOFs2. 單擊單擊OK1. 拾取將要耦合的結(jié)點(diǎn)拾取將要耦合的結(jié)點(diǎn)3. 輸入耦合設(shè)置參考輸入耦合設(shè)置參考號(hào)，選擇自由度卷號(hào)，選擇自由度卷標(biāo)標(biāo).4. 單擊單擊OK.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義54在零偏移量的一組節(jié)點(diǎn)之間生成附加耦合關(guān)系在零偏移量的一組節(jié)點(diǎn)之間生成附加耦合關(guān)系:Main Menu: Preprocessor Coupling / Ceqn Gen w/Same Nodes3. 單擊單擊OK1. 輸入現(xiàn)存耦合輸入現(xiàn)存耦合設(shè)置的參考號(hào)設(shè)置的參考號(hào).2. 對每個(gè)設(shè)置指定對每個(gè)設(shè)置指定新的自由度卷標(biāo)新的自由度卷標(biāo).有限元分析及應(yīng)用講義有限元分析及應(yīng)

36、用講義55在同一位置的節(jié)點(diǎn)之間自動(dòng)生成耦合關(guān)系在同一位置的節(jié)點(diǎn)之間自動(dòng)生成耦合關(guān)系:Main Menu: Preprocessor Coupling / Ceqn Coincident Nodes1. 指定自由度卷標(biāo)指定自由度卷標(biāo).2. 指定節(jié)點(diǎn)位置的指定節(jié)點(diǎn)位置的容差容差3. 單擊單擊OK有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義56在此練習(xí)中，將用耦合在此練習(xí)中，將用耦合/約束選項(xiàng)在兩部分間產(chǎn)約束選項(xiàng)在兩部分間產(chǎn)生耦合生耦合DOF設(shè)置來模擬設(shè)置來模擬接觸問題接觸問題1. 恢復(fù)數(shù)據(jù)庫恢復(fù)數(shù)據(jù)庫cpnorm.db1，并在圖形窗，并在圖形窗口中畫單元口中畫單元.2. 在重合節(jié)點(diǎn)的所有節(jié)點(diǎn)對上建立在

37、重合節(jié)點(diǎn)的所有節(jié)點(diǎn)對上建立UY耦合耦合關(guān)系關(guān)系a. 選擇耦合重合的結(jié)點(diǎn)選擇耦合重合的結(jié)點(diǎn).b. 拾取拾取UY3. 求解并進(jìn)行后處理求解并進(jìn)行后處理有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義57約束方程定義節(jié)點(diǎn)自由度之間的線性關(guān)系約束方程的特點(diǎn)自由度卷標(biāo)的任意組合.任意節(jié)點(diǎn)號(hào).任意實(shí)際的自由度方向在不同的節(jié)點(diǎn)上ux可能不同.例Constant = Coef1 * DOF1 + Coef2 * DOF2 + .有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義581.連接不同的網(wǎng)格:實(shí)體與實(shí)體的界面2D或3D相同或相似的單元類型單元面在同一表面上，但結(jié)點(diǎn) 位置不重合有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義59l建立轉(zhuǎn)動(dòng)自由度和移動(dòng)自由度之間的關(guān)系2.連接不相似的單元類型連接不相似的單元類型:殼與實(shí)體殼與實(shí)體垂直于殼或?qū)嶓w的梁垂直于殼或?qū)嶓w的梁.有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義603. 建立剛性區(qū)在某些特殊情況下，全剛性區(qū)給出了約束方程的另一種應(yīng)用全剛性區(qū)和部分剛性區(qū)的約束方程都可由程序自動(dòng)生成4. 過盈裝配與接觸耦合相似，但在兩個(gè)界面之間允許有過盈量或穿透典型方程:0.01 = UX (node 51) - UX (node 251)有限元分析及應(yīng)用講義有限元分析及應(yīng)用講義61人工建立約束方程的菜單路徑：人工建立約束方程的菜單路徑：Main Menu: Preproc