Ansys第36例熱應(yīng)力分析間接法實(shí)例—液體管路

1、第36例熱應(yīng)力分析(間接法)實(shí)例一液體管路本例介紹了利用間接法進(jìn)行熱應(yīng)力計(jì)算的方法和步驟：首先進(jìn)行熱分析得到結(jié)構(gòu)節(jié)點(diǎn)溫度分布，然后把溫度作為載荷施加到結(jié)構(gòu)上并進(jìn)行結(jié)構(gòu)分析。36.1 概述利用間接法計(jì)算熱應(yīng)力，首先進(jìn)行熱分析，然后進(jìn)行結(jié)構(gòu)分析。熱分析可以是瞬態(tài)的，也可以是穩(wěn)態(tài)的，需要將熱分析求得的節(jié)點(diǎn)溫度作為體載荷施加到結(jié)構(gòu)上。當(dāng)熱分析是瞬態(tài)的時(shí)，需要找到溫度梯度最大的時(shí)間點(diǎn)，并將該時(shí)間點(diǎn)的結(jié)構(gòu)溫度場(chǎng)作為體載荷施加到結(jié)構(gòu)IC由于間接法可以使用所有熱分析和結(jié)構(gòu)分析的功能，所以對(duì)于大多數(shù)情況都推薦使用該萬法。間接法進(jìn)行熱應(yīng)力計(jì)算的主要步驟如下。36.1.1 熱分析瞬態(tài)熱分析的過程在前例已經(jīng)介紹過，下

2、面介紹穩(wěn)態(tài)熱分析。A穩(wěn)態(tài)熱分析用于研究穩(wěn)定的熱載荷對(duì)結(jié)構(gòu)的影響，有時(shí)還用于瞬態(tài)熱分析時(shí)計(jì)算初始溫度場(chǎng)。穩(wěn)態(tài)熱分析王要步驟如下。土0J1 .建模穩(wěn)態(tài)熱分析的建模過程與其他分析相似，包括定義單元類型、定義單元實(shí)常數(shù)、定義材C料特性、建立幾何模型和劃分網(wǎng)格等。但需注意的是：穩(wěn)態(tài)熱分析必須定義材料的導(dǎo)熱系數(shù)。2 .施加載荷和求解(1)指定分析類型。MainMenuSolutionAnalysisTypNewAnalysis,選擇Static.(2)施加載荷?？梢允┘拥妮d荷有恒定的溫度、熱流率、對(duì)流、熱流密度、生熱率，MainMenufSolutionfDefineLoadsfApplyfThermal

3、.(3)設(shè)置載荷步選項(xiàng)。普通選項(xiàng)包括時(shí)間(用于定義載荷步和子步)、每一載荷步的子步數(shù)，以及階躍選項(xiàng)等，MainMenuSolutionLoadStepOpts-Time/FrequencTime一TimeStep.非線性選項(xiàng)包括：一迭代次數(shù)(默認(rèn)25),MainMenuSolutionLoadStepOpts-NonlinearEquilibriumIter;打開自動(dòng)時(shí)間步長，MainMenuSolutionLoadStepOpts-Time/FrequencfTimefTimeStep等.輸出選項(xiàng)包括：控制打印的輸出，MainMenuSolutionLoadStepOptsOutputCtr

4、lsfSoluPrintout;控制結(jié)果文彳的輸出，MainMenuSolutionfLoadStepOptsfOutputCtrls一DB/ResultsFileo(4)設(shè)置分析選項(xiàng)。牛頓拉普森選項(xiàng)(僅對(duì)非線性分析有用)。選擇求解器。確定絕對(duì)零度(在熱輻射分析時(shí)使用)。(5)求解。3 .查看結(jié)果穩(wěn)態(tài)熱分析使用POST1后處理器查看結(jié)果，對(duì)于熱應(yīng)力計(jì)算而言，查看結(jié)果不是必要的。36.1.2 結(jié)構(gòu)分析與一般的結(jié)構(gòu)分析過程大致相同，不同的有以下幾點(diǎn)。4一I)(1)轉(zhuǎn)換單兀。/將熱單元轉(zhuǎn)換為相應(yīng)的結(jié)構(gòu)單元，并設(shè)置需要的單元選項(xiàng)，/MainMenufPreprocessorfElementTypefS

5、witchElemTypeoO日，2)設(shè)置結(jié)構(gòu)分析的材料屬性，例如，熱膨脹系數(shù)、彈性模量、泊松比等。設(shè)置前處理細(xì)節(jié)，例如，節(jié)點(diǎn)耦合等，MainMenufPreprocessorfMaterialPropsMaterialModels。7j(3)讀入熱分析得到的節(jié)點(diǎn)溫度，工*MainMenufSolutionfDefineLoadsfApply-TOStructural一TemperaturefFromThermAnaly。IkJ(4)指定參考溫度，當(dāng)計(jì)算熱膨脹大小時(shí)，溫度差等于節(jié)點(diǎn)溫度減去參考溫度，MainMenufSolutionfDefineLoadsfSettingsfReference

6、Tempo36.2 問題描述某液體管路內(nèi)部通有液體，外部包裹有保溫層，保溫層與空氣接觸如圖36-1所示。已知管路由鑄鐵制造，其導(dǎo)熱系數(shù)為70W/(m.C),彈性模量為2xl011N/m2,泊松比為0.3,熱膨脹系數(shù)為1.2x10-6/C；保溫層的導(dǎo)熱系數(shù)為0.02W/(m.C),彈性模量為2x1010N/m2,泊松比為0.4,熱膨脹系數(shù)為1.2x10-6/C；管路內(nèi)液體壓力為0.3MPa,溫度為70C,對(duì)流換熱系數(shù)為1W/(m2.C)；空氣溫度為-40C,對(duì)流換熱系數(shù)為0.5W/(m2.C)。試分析管路內(nèi)熱應(yīng)力的情況。V)36.3 分析步驟C)jJ*36.3.1 改變?nèi)蝿?wù)名Wz拾取菜單Util

7、ityMenufFile-ChangeJobname,彈出如圖36-2所示的對(duì)話框，在“/(、FILNAM”文本框中輸入EXAMPLE36單擊“OK”按鈕。圖36-2改變?nèi)蝿?wù)名對(duì)話框36.3.2 選擇單元類型拾取菜單MainMenufPreprocessorfElementTypAdd/Edit/Delete,彈出ElementTypes”對(duì)話框，單擊其“Add,”按鈕，彈出如圖36-3所示的對(duì)話框，在左側(cè)列表中選“ThermalSolid”,在右側(cè)列表中選“8node77,單擊OK按鈕，返回到ElementTypes”對(duì)話框，選中“TypelPLANE7Z,單擊“Options”按鈕，彈出如

8、圖36-4所示的對(duì)話框，選擇“K3”“Close”下拉列表框?yàn)锳xisymmetric,單擊OK按鈕，然后單擊如圖36-3所示對(duì)話框的按鈕。圖36-3單元類型庫對(duì)話框圖36-4單元選項(xiàng)對(duì)話框廣36.3.3 定義材料模型拾取菜單MainMenuPreprocessorMaterialPropsMaterialModels,彈出如圖36-5Jfjhr所示的對(duì)話框，在右側(cè)列表中依次拾取Thermal、Conductivity、Isotropic,彈出如圖36-6所示的對(duì)話框，在KXX文本框中輸入70（導(dǎo)熱系數(shù)），單擊OK”按鈕。單擊如X圖36-5所示對(duì)話框的菜單項(xiàng)MaterialNewModel,單

9、擊所彈出DefineMaterialID”對(duì)話框中的“OK”按鈕，然后重復(fù)定義材料模型1時(shí)的各步驟，定義材料模型2的導(dǎo)熱系數(shù)為0.02,最后關(guān)閉如圖36-5所示的對(duì)話框。圖36-5材料模型對(duì)話框餐圖36-6材料特性對(duì)話框36.3.4 創(chuàng)建矩形面，模擬金屬管路和保溫層拾取采單MainMenu-Preprocessor-Modeling-Create-Areas-RectanglefByDimension,彈出如圖36-7所示的對(duì)話框，在“XI,X2文本框中分別輸入0.14,0.15,在“Y1,Y2文本框中分別輸入0,1,單擊“Apply”按鈕。再次彈出如圖36-7所示的對(duì)話框,在“刈X2”文本框

10、中分別輸入0.15,0.2,在“Y1,Y2”文本框中分別輸入0,1,單擊“按鈕。OK”按鈕。圖36-7創(chuàng)建矩形對(duì)話框JX36.3.5 黏結(jié)面iry拾取菜單MainMenuPreprocessorModelingOperateBooleansfGlue-JAreas,彈出拾取窗口，單擊“PickAll”按鈕。36.3.6 顯小面號(hào)拾取菜單UtilityMenufPlotCtrlsfNumbering,在所彈出:PlotNumberingControls”對(duì)話框中，將AreaNumbers（面號(hào)）打開，單擊OK按鈕。36.3.7 劃分單元JrA拾取菜單MainMenuPreprocessorMes

11、hingMeshTool,彈出如圖36-8所示的對(duì)話框。選擇ElementAttributesAreas”，單擊下拉列表框后面的“Set”鈕，彈出拾取窗口，選擇面1,單擊拾取窗口中的OK按鈕，彈出AreasAttributes話框，選擇“MAT”下拉列表框?yàn)?,單擊“Apply”按鈕，再次彈出拾取窗口，選擇面()單擊拾取窗口中的“OK按鈕，選擇“AreasAttributes對(duì)話框的“MAT”下拉列表框?yàn)榘磳?duì)3,2,單擊“OK按鈕。圖36-8單元工具對(duì)話框2f（SlcbalzIESIZEGobal.IlflH斗士h小心3%城口悵（jpplifsodytounsized5lines|SIZEEe

12、mcmcdg,elerlfi0皿ND!VNdcfelementdmons-Q（usedonlyitele-nented9Eng*hS1ZLisblsrkcrzero）OKCaredHdp圖36-9單元尺寸對(duì)話框單擊SizeControls區(qū)域中Global”后面的Set”按鈕，彈出如圖36-9所示的對(duì)話L1框，在“SIZE文本框中輸入0.01,單擊“OK”按鈕。在如圖36-8所示對(duì)話框中的“Mesh”區(qū)域，選擇單元形狀為“Quad”（四邊形），選擇劃分單元的方法為“Mapped”（映射），單擊“Mesh”按鈕，彈出拾取窗口，拾取面1和面3,單擊“OK按鈕。36.3.8 施加對(duì)流邊界條件r、拾取

13、菜單MainMenuSolutionDefineLoadsfApplyfThermalfConvectionfOnLines,彈出拾取窗口，拾取最左面一條直線（即x坐標(biāo)最小的直線），單擊“OK”按鈕，彈Jrf出如圖36-10所示的對(duì)話框，在“VAL1”文本框中輸入1,在“VAL2I”文本框中輸入70,單擊“OK按鈕。再次執(zhí)行命令，彈出拾取窗口，拾取最右面一條直線（即x坐標(biāo)最大的+小斗一,，O）:，、人、直線），單擊“OK”按鈕，彈出如圖、36-10所示的對(duì)話框，在“VAL1”文本框中輸入0.5,在“VAL2r文本框中輸入-40,單擊”O(jiān)K“按鈕。圖36-10施加對(duì)流邊界條件對(duì)話框436.3.9

14、求解r拾取菜單框中的“OKMainMenufSolutionfSoke-CurrentLS,單擊“SolveCurrentLoadStep”對(duì)話按鈕，當(dāng)出現(xiàn)Solutionisdone!”提示時(shí)，求解結(jié)束。以下進(jìn)行結(jié)構(gòu)分析。36.3.10 轉(zhuǎn)換熱單元為結(jié)構(gòu)單元拾取菜單MainMenufPreprocessorfElementTypefSwitchElemType,彈出如圖36-11所示的對(duì)話框，選擇Changeelementtype”下拉歹U表框?yàn)門hermaltoStruc,單擊OK按鈕。圖36-11轉(zhuǎn)換單元類型對(duì)話框36.3.11 設(shè)定單元軸對(duì)稱選項(xiàng)拾取菜單MainMenufPreproc

15、essorfElementTypAdd/Edit/Delete彈出ElementTypes”對(duì)話框，單擊其“Options”按鈕，彈出如圖36-12所示的對(duì)話框，選擇“K3”下拉列表框?yàn)锳xisymmetric,單擊OK按鈕，然后單擊“ElementTypes”對(duì)話中1的Close按鈕。圖36-12單元選項(xiàng)對(duì)話框36.3.12 定義村料模型拾取菜單MainMenuPreprocessorMaterialPropsMaterialModels,彈出如圖36-5所示的對(duì)話框，在左側(cè)列疊沖選中MaterialModelNumber1”,在右側(cè)列表中依農(nóng)拾取“Structural、Linear、Ela

16、stic、Isotropic,彈出如圖36-13所示的對(duì)話框,在EX文本框中輸入2e11、（彈性模量），在“PRXY文本框中輸入0.3（泊松比），單擊“OK”按鈕再在如圖36_5所示對(duì)話框的右側(cè)列表中依次拾取Structural、ThermalExpansion、SecantCoefficient、“Isotropic，彈出如圖36-14所示的對(duì)話框，在“ALPX文本框中輸入1.2e-5（線膨脹系數(shù)），然后單擊“OK”按鈕。FALinearh-otrcpicPropertieslorMateriNumber1圖36-13材料特性對(duì)話框（一）0圖36-14材料特性對(duì)話框（二）選擇如圖36-5所示

17、對(duì)話框左側(cè)列表中的MaterialModelNumber2,用與MaterialModelNumber1同樣的方法，定義彈性模量EX為2e10、泊松比PRXY為0.4、線膨脹系數(shù)ALPX為1.2e-6。最后關(guān)閉如圖36-5所示的對(duì)話框。36.3.13 施加約束拾取菜單MainMenufSolutionfDefineLoadsfApplyfStructuralfDisplacementfOn4條直線），單擊“OK”按Lines,彈出拾取窗口，拾取水平方向的所有直線（即長度較短的鈕，彈出如圖36-15所示的對(duì)話框，在列表中選擇“UY”，單擊“OK”按鈕。圖36-15施加約束對(duì)話框36.3.14 施

18、加壓力載荷/拾取菜單MainMenufSolution-DefinefLoadsApplyStructuralfPressurefOn拓人仃，人仃日I，一日A/4田+，拓Lines,彈出拾取窗口，拾取最左面一條直線（即x坐標(biāo)最小的直線），單擊OK按鈕，彈出如圖36-16所示的對(duì)話框，在“VALUE文本框中輸入3e5,單擊“OK按鈕。圖36-16施加壓力載荷對(duì)話框36.3.15 讀入熱分析所得到的節(jié)點(diǎn)溫度拾取菜單MainMenu-Solution-DefineLoads-Apply-Structural-TemperaturefFrom或者單擊“Browse”按鈕在ANSYSX作文件夾中選擇同一

19、文件，單擊“OK按鈕。圖36-17從熱分析讀入節(jié)點(diǎn)溫度對(duì)話框ThermAnaly,彈出如圖36-17所示的對(duì)話框，在“Fame”文本框中輸入EXAMPLE36Roth36.3.16 指定參考溫度拾取菜單MainMenufSolutionfDefineLoadsfSettingsfReferenceTemp,彈出如圖在計(jì)算熱膨脹大小時(shí)，溫度差等于節(jié)點(diǎn)溫度減去參考溫度，參考溫度的默認(rèn)值為0。a*36-18所示的對(duì)話框，在“TREF文本框中輸入20,單擊“OK按鈕。圖36-18定義參考溫度對(duì)話框36.3.17 求解拾取菜單MainMenufSolutionfSoke-CurrentLS,單擊“Sol

20、veCurrentLoadStep”對(duì)話框中的OK按鈕，當(dāng)出現(xiàn)Solutionisdone!”提示時(shí)，求解結(jié)束，即可查看即可查看結(jié)果。36.3.18 查看結(jié)果，用等圖線顯小VonMises應(yīng)力拾取菜單MainMenuGeneralPostprocPlotResultsContourPlotNodalSolu,彈出“ContourNodalSolutionData”對(duì)話框，在列表中依次選擇NodalSolutionfStressfvonMisesSEQV（第四強(qiáng)度理論的等效應(yīng)力），單擊“OK按鈕。結(jié)果如圖36-19所示，可以看出，VonMises應(yīng)力的最大值為0.588x108Pa,即58.8M

21、Pa。ANSYSJAN11201622:10*27圖36-19液體管路的VonMises應(yīng)力第37例熱應(yīng)力分析（直接法）實(shí)例-液體管路3MLSOLIRIONSTEP-1SUB-1TTME-1SRQV(AVG)DMX=.450E-045MN=.d59忙于tn3MX:.弓gE芭*UH.459E+0.166E4(J3.2H7E+D3.407E+0B.52BE+0B.lOiEfD6.227E+Q6.3A7E+0B.469E+0B.568ErOa與間接法相比，直接法步驟更簡本例介紹了利用直接法進(jìn)行熱應(yīng)力計(jì)算的方法和步驟,I捷，但要使用耦合單元進(jìn)行非線性分析，計(jì)算量較大。37.1概述與一般的分析過程相似，

22、只是要使用耦合單元，并同時(shí)施加結(jié)構(gòu)分析和熱分析的約束和載荷。JlX.37.2問題描述為了與間接法進(jìn)行熱應(yīng)力計(jì)算相對(duì)比，本例采用了與第36例相同的模型一液體管路。7.3分析步驟37.3.1 改變?nèi)蝿?wù)名拾取菜單UtilityMenuFile一ChangeJobname,彈出如圖37-1所示的對(duì)話框，在“/FILNAM”文本框中輸入EXAMPLE37單擊”O(jiān)K“按鈕。圖37-1改變?nèi)蝿?wù)名對(duì)話框37.3.2 選擇單元類型拾取菜單MainMenufPreprocessorfElementTypefAdd/Edit/Delete,彈出ElementTypes”對(duì)話框，單擊其“Add.。，”按鈕，彈出如圖3

23、7-2所示的對(duì)話框，在左側(cè)列表中選CoupledField，在右側(cè)列表中選VectorQuad13，單擊OK按鈕，返回至1NElementTypes對(duì)話框，選中“TypelPLANE13,單擊“Options”按鈕，彈出如圖37-3所示的對(duì)話框，選VJ擇“K1”下拉列表框?yàn)閁XUYTEMPAZ,選擇K3”下拉列表框?yàn)锳xisymmetric”，單擊“OK”按鈕，單擊如圖37-2所示對(duì)話框中的“Close”按鈕。圖37-2:單元類型庫對(duì)話框APLANEl3elementtype叩tinn零圖37-3單元選項(xiàng)對(duì)話框37.3.3 定義材料模型1拾取采單MainMenufPreprocessorfMa

24、terialPropsfMaterialModels,彈出如圖37-4所示的對(duì)話框，在右側(cè)列表中依次拾取Thermal、Conductivity、Isotropic,彈出如圖37-5所示的對(duì)話框，在“KXX文本框中輸入70（導(dǎo)熱系數(shù)），單擊“OK”按鈕；再在右側(cè)列表中依次拾取Structural、Linear、Elastic、Isotropic,彈出如圖37-6所示的對(duì)話框，在“EX”文本框中輸入2ell（彈性模量），在“PRXY文本框中輸入0.3（泊松比），單擊“OK”按鈕；再在右側(cè)列表中依次拾取LStructural、ThermalExpansion、SecantCoefficient”、

25、“Isotropic，彈出如圖37-7所示的對(duì)話框，在“ALPX文本框中輸入1.2e-5（線膨脹系數(shù)），單擊“OK按鈕。單擊如圖37-4。所示對(duì)話框的菜單項(xiàng)MaterialNewModel,單擊所彈出“DefineMaterialID”對(duì)話框中的“OK”按鈕，然后重復(fù)定義材料模型1時(shí)的各步驟，定義材料模型2的導(dǎo)熱系數(shù)KXX為0.02、彈性模量EX為2e10、泊松比PRXY為0.4、線膨脹系數(shù)ALPX為1.2e-6。最后關(guān)閉如圖37-4所示的對(duì)話框。圖37-4材料模型對(duì)話框0圖37-5材料特性對(duì)話框（一）圖37-6材料特性對(duì)話框（二）37-7材料特性對(duì)話框（三）37.3.4 創(chuàng)建矩形面，模擬金屬

26、管路和保溫層拾取菜單MainMenu-Preprocessor-Modeling-Create-Areas-RectanglefByDimension,彈出如圖37-8所示的對(duì)話框，在“X1,x2”文本框中分別輸入0.14,0.15,在“Y1,Y2文本框中分別輸入0,1,單擊“Apply”按鈕。再次彈出如圖37-8所示的對(duì)話框，在“Xl,X2”文本框中分別輸入0.15,0.2,在“Y1,Y2”文本框中分別輸入0,1,單擊”O(jiān)K“按鈕。圖37-8.創(chuàng)建矩形對(duì)話框37.3.5黏結(jié)面拾取菜單MainMenufPreprocessorfModelingfOperatefBooleansfGluefAr

27、eas。彈kT出拾取窗口，單擊“PickAll”按鈕。37.3.6 顯示面號(hào)拾取菜單UtilityMenuPlotCtrlsNumbering,將所彈出PlotNumberingControls”對(duì)話框中的AreaNumbers（面3）打開，單擊OK按鈕。37.3.7 劃分單元拾取菜單MainMenufMeshingfMeshTool,彈出如圖37-9所示的對(duì)話框。選擇ElementAttributes”下拉列表框?yàn)椤癆reas”，單擊下拉列表框后面的“Set”按鈕，彈出拾取窗口，選擇面1,單擊拾取窗口中的OK按鈕，彈出AreasAttributes”對(duì)話框，選擇MAT”下拉列表框?yàn)?,單擊“

28、Apply”按鈕，再次彈出拾取窗口，選擇面3,單擊拾取窗口中的“OK”按鈕，選擇“AreasAttributes對(duì)話框的“MAT”下拉列表框?yàn)?,單擊OK按鈕。單擊“SizeControls區(qū)域中Global”后面的Set”按鈕，弓t出如圖37-10所示的對(duì)話框，在“SIZE文本框中輸入0.01,單擊“OK按鈕。在如圖37-9所示對(duì)話框的“Mesh”區(qū)域，選擇單元形狀為“Quad（四邊形），選擇劃分單元的方法為“Mapped”（映射），單擊“Mesh”按鈕，彈出拾取窗口，拾取面1和面3,單擊“OK按鈕。圖37-9單元工具對(duì)話框陌在Globalx坐標(biāo)最小的直線），單擊“OK按鈕,A引口MlElA

29、mentsi?i=eferiertanddivijion&(appliesonlySI7FFlfEnnigpFiragrhND1VNo,of-kinent：iivisions-unlyifclernicredijeergth,SIZE.、Hankorzrrc)Cancel圖37-10單元尺寸對(duì)話框37.3.8 流邊界條件施加對(duì)流拾取菜單MainMenufSolutionfDefineLoadsApplyThermalConvection一OnLines,彈出拾取窗口，拾取最左面一條直線（即彈出如圖37-11所示的對(duì)話框，在“VAL文本框中輸入1（對(duì)流換熱系數(shù)），在“VAL2I”文本框中輸入70（溫度），單擊“OK”按鈕，再次執(zhí)行命令，彈出拾取窗口，拾取最右面一條直線（即x坐標(biāo)最大的直線），單擊“OK”按鈕，彈出如圖37-11所示的對(duì)話框，在“VAL1”文本框中輸入0.5（對(duì)流換熱系數(shù)），在“VAL2I”文本框中輸入-40（溫度），單擊“OK按圖37-11施加對(duì)流邊界條件對(duì)話框G3.3.9施加約束、九u拾取菜單MainMenufSolution