利用分析結(jié)構(gòu)靜力

1、利用分析結(jié)構(gòu)靜力第1頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2本章的目標(biāo)是，假如已將幾何模型劃分網(wǎng)格, 應(yīng)如何加載、求解.LessonObjectives目 標(biāo) 第1課. 載荷6-1.列表和分類載荷。6-2.在實(shí)體模型上完成下列操作:a.加載.b.校驗(yàn)載荷.c.刪除載荷.第2課. 求解6-3.描述求解過程.第3課. 結(jié)果后處理6-4.描述ANSYS后處理中觀看結(jié)果的各種功能.6-5.描述靜力分析結(jié)果后處理的五個(gè)步驟. 6-6 實(shí)例第2頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一3第 1 課 加 載 第3頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一4

2、載荷分類Objective6-1.列表和分類載荷ANSYS中的載荷可分為:自由度DOF - 定義節(jié)點(diǎn)的自由度（ DOF ） 值 (結(jié)構(gòu)分析_位移、熱分析_ 溫度、電磁分析_磁勢等)集中載荷 - 點(diǎn)載荷 (結(jié)構(gòu)分析_力、熱分析_ 熱導(dǎo)率、電磁分析_ magnetic current segments)面載荷 - 作用在表面的分布載荷 (結(jié)構(gòu)分析_壓力、熱分析_熱對流、電磁分析_magnetic Maxwell surfaces等)體積載荷 - 作用在體積或場域內(nèi) (熱分析_ 體積膨脹、內(nèi)生成熱、電磁分析_ magnetic current density等)慣性載荷 - 結(jié)構(gòu)質(zhì)量或慣性引起的載荷

3、 (重力、角速度等)第4頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一5加載Objective6-2a.加載.可在實(shí)體模型或 FEA 模型 (節(jié)點(diǎn)和單元) 上加載.在關(guān)鍵點(diǎn)處約束實(shí)體模型沿線均布的壓力在關(guān)鍵點(diǎn)加集中力在節(jié)點(diǎn)處約束FEA 模型沿單元邊界均布的壓力在節(jié)點(diǎn)加集中力第5頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一6加載 (續(xù))幾何模型加載獨(dú)立于有限元網(wǎng)格. 重新劃分網(wǎng)格或局部網(wǎng)格修改不影響載荷.加載的操作更加容易 ,尤其是在圖形中直接拾取時(shí).直接在實(shí)體模型加載的優(yōu)點(diǎn):Guidelines第6頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一7加載 (續(xù)

4、)無論采取何種加載方式,ANSYS求解前都將載荷轉(zhuǎn)化到有限元模型.因此, 加載到實(shí)體的載荷將自動(dòng)轉(zhuǎn)化到 其所屬的節(jié)點(diǎn)或單元上。實(shí)體模型加載到實(shí)體的載荷自動(dòng)轉(zhuǎn)化到其所屬的節(jié)點(diǎn)或單元上FEA 模型沿線均布的壓力 均布?jí)毫D(zhuǎn)化到以線為邊界的各單元上第7頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一8加載 (續(xù))注意到這是 很長的菜單, 對于結(jié)構(gòu)分析，部分菜單呈暗淡灰色，表示不屬于結(jié)構(gòu)分析的范疇。 (ANSYS 可由模型中的單元類型識(shí)別分析類型) 實(shí)體模型加載:Main Menu: Solution -Loads- Apply 步驟1. .2. .3. . 說明: 可通過在preferen

5、ces 中選擇適當(dāng)?shù)姆治鲱愋瓦^濾菜單中的選項(xiàng)。第8頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一9加載 (續(xù))輸入一個(gè)壓力值即為 均布載荷, 兩個(gè)數(shù)值 定義坡度壓力說明：壓力數(shù)值為正表示其方向指向表面Main Menu: Solution -Loads- Apply Pressure On Lines加載面力載荷拾取Line第9頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一10加載 (續(xù))VALI = 500VALI = 500VALJ = 1000VALI = 1000VALJ = 500500L3500L31000500L31000500坡度壓力載荷沿起始關(guān)鍵點(diǎn)(I

6、) 線性變化到第二個(gè)關(guān)鍵點(diǎn) (J)。 如果加載后坡度的方向相反, 將兩個(gè)壓力數(shù)值顛倒即可。加載面力載荷（續(xù)）第10頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松11加載 (續(xù))軸對稱載荷可加載到具有對稱軸的3-D 結(jié)構(gòu)上。3-D 軸對稱結(jié)構(gòu)可用一2-D 軸對稱模型描述。加載軸對稱載荷10” 直徑5” 半徑軸對稱模型3-D 結(jié)構(gòu)對稱軸第11頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一12加載 (續(xù))加載 軸對稱載荷, 注意以下方面:載荷數(shù)值 (包括輸出的反力) 基于360度轉(zhuǎn)角的3-D結(jié)構(gòu)。在右圖中，軸對稱模型中的載荷是3-D結(jié)構(gòu)均布面力載荷的總量。

7、Total Force = 2pr = 47,124 lb.準(zhǔn)則3-D 結(jié)構(gòu)2-D 有限元模型Axis of symmetry第12頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一13加載 (續(xù))在關(guān)鍵點(diǎn)加載位移約束: 加載約束載荷Main Menu: Solution -Loads- Apply -Structural- Displacement On Keypoints +procedure1. .2. .3. . Expansion option 可使相同的載荷加在位于兩關(guān)鍵點(diǎn)連線的所有節(jié)點(diǎn)上拾取keypoints例要固定一邊，只要拾取關(guān)鍵點(diǎn)6、7,并設(shè)置 all DOFs =

8、 0 和 KEXPND = yes.K6K7第13頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一14加載 (續(xù))加載約束載荷（續(xù)）在線和面上加載位移約束: Main Menu: Solution -Loads- Apply -Structural- Displacement On Lines + OR On Areas+步驟1. .2. .3. . 拾取 lines 拾取areas第14頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一15校驗(yàn)載荷實(shí)體模型載荷顯示在幾何模型上 (體、面、線或關(guān)鍵點(diǎn))有限元模型載荷在畫節(jié)點(diǎn)或單元時(shí)顯示通過 plotting畫出載荷: Util

9、ity Menu: PlotCtrls Symbols .或通過 listing列表載荷:Utility Menu: List Loads步驟1. .2. .3. . Objective6-2b. 校驗(yàn)載荷第15頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一16將載荷轉(zhuǎn)化到有限元模型上說明: 只有到求解初始化時(shí)，才將模型中的載荷自動(dòng)轉(zhuǎn)化到有限元模型中的節(jié)點(diǎn)和單元上。Procedure1. .2. .3. . 下面將載荷轉(zhuǎn)化到節(jié)點(diǎn)和單元上，不進(jìn)行求解:Main Menu: Solution -Loads-Operate這些選項(xiàng)出現(xiàn)的信息大致相同第16頁，共58頁，2022年，5月20

10、日，13點(diǎn)51分，星期一17刪除載荷Main Menu: Solution -Loads- Delete All Load Data 選項(xiàng)可同時(shí)刪除模型中的任一類載荷。individual entities by picking 選項(xiàng)只刪除模型選定的載荷。而.Procedure1. .2. .3. . 6-2c. 刪除載荷Objective第17頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一18刪除載荷（續(xù)）當(dāng)刪除實(shí)體模型時(shí), ANSYS 將自動(dòng)刪除其上所有的載荷實(shí)體模型FEA 模型l刪除線上的均布?jí)毫?自動(dòng)刪除以線為邊界的各單元均布?jí)毫?記住這 一關(guān)系?第18頁，共58頁，202

11、2年，5月20日，13點(diǎn)51分，星期一19刪除載荷（續(xù)）兩關(guān)鍵點(diǎn)的擴(kuò)展位移約束載荷例外：刪除兩點(diǎn)的約束 只刪除了兩角點(diǎn)（ CORNER ）約束， 而加載時(shí)擴(kuò)展的 （ inside ） 節(jié)點(diǎn)約束必須手工刪除.實(shí)體 模型FEA 模型l第19頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一20第 2 課 求 解 第20頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一21求解求解結(jié)果保存在數(shù)據(jù)庫中并輸出到結(jié)果文件 (Jobname.RST, Jobname.RTH, Jobname.RMG, or Jobname.RFL)結(jié)果文件結(jié)果數(shù)據(jù)數(shù)據(jù)庫求解器結(jié)果輸入數(shù)據(jù)第21頁，共58頁，

12、2022年，5月20日，13點(diǎn)51分，星期一2022/9/322模態(tài)分析諧波分析瞬態(tài)分析頻譜分析子結(jié)構(gòu)分析 求解前使用該對話框?qū)⒏鱾€(gè)控制項(xiàng)設(shè)置好，就可以求解了。命令：SOLCONTROL，Key1，Key2，Key3，VtolGUI：MainMenuSolutionAnalysis TypeSolution Ctrl第22頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一23求解時(shí)模型是否準(zhǔn)備就緒?在求解初始化前，應(yīng)進(jìn)行分析數(shù)據(jù)檢查，包括下面內(nèi)容:統(tǒng)一的單位單元類型和選項(xiàng)材料性質(zhì)參數(shù)考慮慣性時(shí)應(yīng)輸入材料密度熱應(yīng)力分析時(shí)應(yīng)輸入材料的熱膨脹系數(shù)實(shí)常數(shù) (單元特性)單元實(shí)常數(shù)和材料類型的設(shè)置

13、實(shí)體模型的質(zhì)量特性 (Preprocessor Operate Calc Geom Items)模型中不應(yīng)存在的縫隙殼單元的法向節(jié)點(diǎn)坐標(biāo)系集中、體積載荷面力方向溫度場的分布和范圍熱膨脹分析的參考溫度 (與 ALPX 材料特性協(xié)調(diào)?)第23頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一24求解過程:1. 求解前保存數(shù)據(jù)庫2. 將Output 窗口提到最前面觀看求解信息3. Main Menu: Solution -Solve-Current LS.進(jìn)行求解Objective6-3.描述求解過程Procedure1. .2. .3. . 第24頁，共58頁，2022年，5月20日，13

14、點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松25進(jìn)行求解(續(xù))在求解過程中，應(yīng)將OUTPUT窗口提到最前面。 ANSYS 求解過程中的一系列信息都將顯示在此窗口中，主要信息包括:模型的質(zhì)量特性- 模型質(zhì)量是精確的 - 質(zhì)心和 質(zhì)量矩的值有一定誤差。單元矩陣系數(shù) - 當(dāng)單元矩陣系數(shù)最大/最小值的比率 1.0E8 時(shí)將預(yù)示模型中的材料性質(zhì)、實(shí)常數(shù)或幾何模型可能存在問題。當(dāng)比值過高時(shí)，求解可能中途退出。模型尺寸和求解統(tǒng)計(jì)信息。匯總文件和大小。第25頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松26進(jìn)行求解(續(xù))沒有獲得結(jié)果的原因是什么? 往往是求解輸入的模型不完整或存在錯(cuò)誤，

15、典型原因有: 約束不夠! (通常出現(xiàn)的問題)。當(dāng)模型中有非線性單元 (如縫隙 gaps、滑塊sliders、鉸hinges、索cables等），整體或部分結(jié)構(gòu)出現(xiàn)崩潰或“松脫”。材料性質(zhì)參數(shù)有負(fù)值, 如密度或瞬態(tài)熱分析時(shí)的比熱值。未約束鉸接結(jié)構(gòu)，如兩個(gè)水平運(yùn)動(dòng)的梁單元在豎直方向沒有約束。屈曲 - 當(dāng)應(yīng)力剛化效應(yīng)為負(fù)（壓）時(shí)，在載荷作用下整個(gè)結(jié)構(gòu)剛度弱化。如果剛度減小到零或更小時(shí)，求解存在奇異性，因?yàn)檎麄€(gè)結(jié)構(gòu)已發(fā)生屈曲。第26頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一27第 3 課 結(jié)果后處理 第27頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一28結(jié)果的繪圖和列表A

16、NSYS 有兩個(gè)后處理器:通用后處理器 (即 “POST1”) 只能觀看整個(gè)模型在某一時(shí)刻的結(jié)果 (如：結(jié)果的照相 “snapshot”). 時(shí)間歷程后處理器 (即 “POST26”) 可觀看模型在不同時(shí)間的結(jié)果。 但此后處理器只能用于處理瞬態(tài)和/或動(dòng)力分析結(jié)果。在這一課只討論通用后處理器Objective6-4.介紹ANSYS后處理功能第28頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松29結(jié)果的繪圖和列表(續(xù))靜力分析結(jié)果后處理的步驟主要包括:1.繪變形圖2.變形動(dòng)畫3.支反力列表4.應(yīng)力等值線圖5.網(wǎng)格密度檢查GuidelinesObjective6-5.

17、介紹靜力分析結(jié)果后處理的五個(gè)步驟第29頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松30繪變形圖繪出結(jié)構(gòu)在靜力作用下的變形結(jié)果: Main Menu: General Postprocessor Plot Results Deformed Shape.第30頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松31變形動(dòng)畫以動(dòng)畫方式模擬結(jié)構(gòu)在靜力作用下的變形過程: Utility Menu: PlotCtrls Animate Deformed Shape.第31頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開

18、松32支反力在任一方向，支反力總和必等于在此方向的載荷總和。節(jié)點(diǎn)反力列表:Main Menu: General Postprocessor List Results Reaction Solution.第32頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松33應(yīng)力等值線應(yīng)力等值線方法可清晰描述一種結(jié)果在整個(gè)模型中的變化，可以快速確定模型中的 “危險(xiǎn)區(qū)域”。顯示應(yīng)力等值線 :Main Menu: General Postprocessor Plot Results -Contour Plot- Nodal Solution.第33頁，共58頁，2022年，5月20日

19、，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松34應(yīng)力等值線動(dòng)畫結(jié)果動(dòng)畫 :Utility Menu: PlotCtrls Animate Deformed Results第34頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松35關(guān)于PowerGraphics的說明PowerGraphics 特點(diǎn): 快速重畫、圖形輪廓分明。模型顯示光滑、具有相片的真實(shí)感。支持單元類型（ lines、pipes、elbows、 contact 等單元）和幾何實(shí)體 （ lines、areas、 volumes 等）。PowerGraphics 打開 (缺省)PowerGraphics 關(guān)

20、閉第35頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松36檢查網(wǎng)格精度由于網(wǎng)格密度影響分析結(jié)果的精度，因此有必要驗(yàn)證網(wǎng)格的精度是否足夠。有三種方法進(jìn)行網(wǎng)格精度檢查：1.觀察（ Visual inspection ）2.誤差估計(jì)3.將網(wǎng)格加密一倍，重新求解并比較兩者結(jié)果。注意: 有些情況下這種做法不適用。第36頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松37觀察畫出非平均（ unaveraged ）應(yīng)力等值線，例如，畫出單元應(yīng)力而不是節(jié)點(diǎn)應(yīng)力。顯示每個(gè)單元的應(yīng)力尋找單元應(yīng)力變化大的區(qū)域，這些區(qū)域應(yīng)進(jìn)行網(wǎng)格加密。（在MeshTool

21、中對網(wǎng)格加密非常方便， MeshTool 將在后面的內(nèi)容介紹。）Averaged stress contoursUnaveraged stress contours第37頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松38誤差估計(jì)ANSYS 對平均應(yīng)力和非平均應(yīng)力采用幾種不同的誤差計(jì)算方法，誤差估計(jì)只在進(jìn)入后處理前PowerGraphics 被關(guān)閉的情況下進(jìn)行。 (如果進(jìn)入后處理后關(guān)閉 PowerGraphics則ANSYS將重新計(jì)算誤差因子。) 關(guān)閉 PowerGraphics，應(yīng)力等值線圖可顯示應(yīng)力分布和最大最小值范圍，這可表明誤差的大小。通過畫出結(jié)構(gòu)能誤差的

22、等值線圖，可顯示誤差較大的區(qū)域 - 這些區(qū)域需要網(wǎng)格加密。畫出所有單元的應(yīng)力偏差圖，可給出每個(gè)單元的應(yīng)力誤差值。 (平均應(yīng)力和非平均應(yīng)力不同)第38頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一機(jī)械工程學(xué)院 王開松39應(yīng)力平均FEA的計(jì)算結(jié)果包括通過計(jì)算直接得到的初始量和導(dǎo)出量。任一節(jié)點(diǎn)處的DOF結(jié)果 (UX、UY、TEMP等) 是初始量。 它們只是在每個(gè)節(jié)點(diǎn)計(jì)算出來的初始值。其它量，如應(yīng)力應(yīng)變，是由DOF 結(jié)果通過單元計(jì)算導(dǎo)出而得到的。因此，在給定節(jié)點(diǎn)處，可能存在不同的應(yīng)力值。這是由以與此節(jié)點(diǎn)相連的不同單元計(jì)算而產(chǎn)生的?！肮?jié)點(diǎn)結(jié)果”（nodal solution）畫出 的是在節(jié)點(diǎn)

23、處導(dǎo)出量的平均值，而 “單元結(jié)果”（ element solution ） 畫出非平均量。 在同一個(gè)節(jié)點(diǎn)處，相臨單元計(jì)算出的結(jié)果不同。第39頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一40應(yīng)力平均(續(xù))在彈性模量不同的材料交界處，應(yīng)力分量會(huì)不連續(xù)。 (PowerGraphics 自動(dòng)考慮到這一點(diǎn)并對此界面不進(jìn)行平均處理。)在多數(shù)情況下，畫出平均應(yīng)力圖，但有時(shí)要畫出：非平均單元應(yīng)力顯示不連續(xù)的應(yīng)力平均的節(jié)點(diǎn)應(yīng)力顯示連續(xù)的應(yīng)力第40頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一41應(yīng)力平均(續(xù))在不同厚度的殼單元的交界處，大多數(shù)應(yīng)力會(huì)不連續(xù) 。(PowerGraphics

24、 自動(dòng)考慮到這一點(diǎn)并對此界面不進(jìn)行平均處理。)非平均單元應(yīng)力顯示不連續(xù)的應(yīng)力平均的節(jié)點(diǎn)應(yīng)力顯示掩蓋了應(yīng)力的不連續(xù)性第41頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一42應(yīng)力平均(續(xù))在殼單元構(gòu)成的尖角或連接處，某些應(yīng)力分量不連續(xù)。非平均單元應(yīng)力顯示不連續(xù)的應(yīng)力平均的節(jié)點(diǎn)應(yīng)力顯示掩蓋了應(yīng)力的不連續(xù)性第42頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2022/9/3436-6 簡單實(shí)例6.6.1 平面問題：板中圓孔的應(yīng)力集中實(shí)例1：如圖所示為承受雙向拉伸的板件，其中心位置有一個(gè)小圓孔，尺寸(mm)如圖所示。彈性模量E2105N/mm2，泊松比v=0.3拉伸載荷:q=20

25、N/mm平板的厚度：t=20mm解題思路分析：1.屬于平面應(yīng)力問題2.中心帶孔，應(yīng)使用8節(jié)點(diǎn)四邊形單元或三角形單元3.注意單位：尺寸mm，力N，故應(yīng)力N/mm24.最大變形約為0.001mm（忽略孔的影響）,最大應(yīng)力在孔的頂部和底部，大小約為3.9N/mm2，即3.9MPa。依次檢驗(yàn)有限元的分析結(jié)果。第43頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2022/9/344一. 相關(guān)設(shè)置設(shè)置jobname為bracket、Title為”a Example for Bracket”，且過濾參數(shù)為Structural;二.建立網(wǎng)格模型1.創(chuàng)建幾何模型：在XY平面內(nèi)建立一個(gè)矩形和圓，并用布

26、爾subtract得到幾何模型；2.定義單元類型：選擇Structural Solid選項(xiàng)下的Quad 8node 82單元，確定返回最后自動(dòng)得到單元類型為PLANE82。它是8節(jié)點(diǎn)的四邊形單元，是平面4節(jié)點(diǎn)單元PLANE42的高階形式，更適合有曲線邊界的模型。對于本問題，我們需要有厚度的平面應(yīng)力單元，只需單擊單元類型表中的options按鈕，彈出PLANE82選項(xiàng)設(shè)置窗口，如圖所示，在K3對應(yīng)的方框中選擇”Plane strs w/thk“，使得可以設(shè)置板的厚度。第44頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2022/9/3453.定義實(shí)常數(shù)：PreprocessorRea

27、l ConstantsAdd/Edit/Delete。彈出如下對話框,輸入板的厚度204.定義材料特性：EX200000，PRXY0.35.定義網(wǎng)格尺寸：網(wǎng)格邊長25，劃自由網(wǎng)格。如下頁所示 PreprocessorMeshingSize CntrlsManual SizeAll Areas PreprocessorMeshingMeshAreasFree或用MeshTool6. 保存工作第45頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2022/9/346三. 加載和求解1.定義分析類型：SolutionAnalysis TypeNew Analysis ，設(shè)為Static，即

28、結(jié)構(gòu)靜態(tài)分析；2.施加約束：SolutionDefine LoadsApplyStructuralDisplacementon Lines，用鼠標(biāo)點(diǎn)選模型最左側(cè)邊，并全部約束(All DOF)3.施加載荷：板右側(cè)邊緣上有一個(gè)背離平板的20N/mm的均布線載荷，則均布?jí)毫€載荷除以板厚20mm1N/mm2。對對模型右側(cè)邊施加1的均布表面壓力。SolutionDefine LoadsApplyStructuralPressureOn Lines4. 求解：SolutionSolveCurrent LS網(wǎng)格模型圖約束、載荷模型圖第46頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2022

29、/9/347三.檢查計(jì)算結(jié)果，觀察收斂情況,決定是否修改網(wǎng)格模型1.節(jié)點(diǎn)最大應(yīng)力檢查 由于最大應(yīng)力點(diǎn)的應(yīng)力值有解析解，因此可以檢查該點(diǎn)的應(yīng)力值。首先顯示節(jié)點(diǎn)編號(hào)，找出孔部對應(yīng)的節(jié)點(diǎn)：Utility MenuPlotNodes（如果沒有顯示節(jié)點(diǎn)編號(hào)，則Utility MenuPlotCtrlsNumbering，打開節(jié)點(diǎn)編號(hào)），記下與圓頂部對應(yīng)的節(jié)點(diǎn)編號(hào)。如圖所示 列出應(yīng)力值：General PostprocList ResultsNodal SolutionStress，檢查所要考察節(jié)點(diǎn)的SEQV值（等效應(yīng)力值），查結(jié)果可知其大小為2.67N/mm2，與手工計(jì)算的結(jié)果3.9Mpa有較大的差別，

30、因此需要在孔的周圍采用更小的網(wǎng)格尺寸才能獲得更為精確的解。第47頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2022/9/3482. 修改網(wǎng)格模型：網(wǎng)格模型的修改方式很多，這里選擇將孔周圍的單元進(jìn)行網(wǎng)格細(xì)化。運(yùn)用MeshTool的Refine完成，細(xì)化級(jí)別1（稍作細(xì)化）。結(jié)果如下頁圖所示。第48頁，共58頁，2022年，5月20日，13點(diǎn)51分，星期一2022/9/349孔周圍要細(xì)化的單元細(xì)化后的網(wǎng)格3.重新計(jì)算：這時(shí)檢查孔頂部最大等效應(yīng)力（編號(hào)可能會(huì)改變）3.39N/mm2四. 后處理1.繪制變形圖：General PostprocPlot ResultsDeformed ShapeDef+undeformed從圖中可以看出孔的變形情況、整體變形情況，并且從圖中左上角說明得知，最大位移0.00124mm2.繪制等效應(yīng)力云圖： General PostprocPlot ResultsContour plotNodal Solution，在彈出的窗口中選von Mises SEQV，得到等效應(yīng)力云圖如下頁所示。第49頁，共58頁，20