Ansys第例瞬态热分析实例一水箱

1、第 33 例 瞬态热分析实例水箱本例介绍了利用ANSYS进行瞬态热分析的方法和步骤、瞬态热分析时材料模型所包含的 内容，以及模型边界条件和初始温度的施加方法。概述热分析是计算热应力的基础， 热分析分为稳态热分析和瞬态热分析， 稳态热分析将在后 面两个例子中介绍，本例介绍瞬态热分析。瞬态热分析的定义瞬态热分析用于计算系统随时间变化的温度场和其他热参数。 一般用瞬态热分析计算温 度场，并找到温度梯度最大的时间点，将此时间点的温度场作为热载荷来进行应力计算。嚼态热分析的步骤瞬态热分析包括建模、施加载荷和求解、查看结果等几个步骤。1建模瞬态热分析的建模过程与其他分析相似， 包括定义单元类型、 定义单元

2、实常数、 定义材 料特性、建立几何模型和划分网格等。注意：瞬态热分析必须定义材料的导热系数、密度和比热。2. 施加载荷和求解(1) 指定分析类型，Main Menu 宀 SolutionAnalysis TypeNew Analysis, 选择 Transient 。(2) 获得瞬态热分析的初始条件。定义均匀的初始温度场：Main Menu 宀 Solution 宀 Define LoadsSett in gsUniformTemp初始温度仅对第一个子步有效，而用Main Menu宀Solution 宀Define Loads ApplyTThermal Temperature命令施加的温度在

3、整个瞬态热分析过程中均不变，应注意二者的 区别。定义非均匀的初始温度场：如果非均匀的初始温度场是已知的，可以用Main MenutSolution tDefine Loads tApplytInitial ConditntDefine 即 IC 命令施加。非均匀的初始温度场一般是未知的， 此时必须先进行行稳态分析确定该温度场。 该稳态分析与一般的 稳态分析相同。注意：要设定载荷(如已知的温度、热对流等) ，将时间积分关闭，选择 Main MenutSolution tLoad Step Opts tTime/Frequenc tTime Integration tAmplitude Decay

4、; 设定只有一个子步，时间很短(如的载荷步，Main Menu f Solution f Load Step Opts fTime/Frequenc f Timef Time Step。(3) 设置载荷步选项。普通选项包括每一载荷步结束的时间、每一载荷步的子步数、阶跃选项等，选择MainMenif Solution f Load Step Opts f Time/Frequenc f Time-Time Step.非线性选项包括：迭代次数(默认 25),选择 Main Menuf Solution f Load Step Opts f Non li near f Equilibrium Ite

5、r;打开自动时间步长，选择Ma in Menu f Solution f LoadStep Optsf Time/Frequenc f Timef Time Step :将时间积分打开，选择 Main Menuf Solution f Load Step Opts f Time/Frequenc f Time Integrationf Amplitude Decay.输出选项包括： 控制打印的输出， 选择 Main Menuf Solution f Load Step Optsf Output Ctrls f Solu Prin tout;结果文件的输出，选择 Ma in Menif Solu

6、tion f Load Step Optsf Output Ctrls f DB/Results File.(4) 如果均匀的初始温度场是通过稳态分析施加的，则须删除稳态分析时施加的温度载荷， 选择 Main Menu f Solution f Define Loads f Delete f Thermal f Temperature 。(5) 求解。3. 查看结果可以用POST26或 POST1查看结果。问题描述图33-1所示为一个温度为 500C的铁块和一个温度为 400 C的铜块，突然放入温度为 20C的完全绝热的水箱中。忽略水的流动，试分析 1h后铜块和铁块的最高温度，以及铜块 和铁块

7、的温度变化情况。材料热物理性能参数如表33-1所示。图33-1 水箱示意图表33-1 材料热物理性能参数热性能参数铜铁水导热系数/W(m. C)383702密度/kg/m388897833966比热 /J/(kg.C)3904484185分析步骤拾取菜单 Utility Menu 宀File宀Cha nge Job name,在所弹出对话框的“ /FILNAM ”文 本框中输入 EXAMPLE3,单击“ 0K”按钮。图 33-2 改变任务名对话框拾取菜单 Ma in Menu PreprocessorEleme nt Type Add/Edit/Delete ，弹出如图 33-3所示的对话框，

8、单击“Add”按钮，弹出如图33-4所示的对话框，在左侧列表中选“Thermal Solid ”，在右侧列表中选“ 8node 77 ”，单击“ 0K按钮，返回到如图33-3所示的对话框，单击“ Close ”按钮关闭对话框。图 33-3单元类型对话框图 33-4单元类型库对话框拾取菜单 Main Menu宀Preprocessor 宀Material Props Material Models，弹出如图33-5 所示的对话框,在右侧列表中依次拾取“ Thermal ”、“Conductivity ”、“Isotropic ”, 弹出如图33-6所示的对话框，在“ KXX文本框中输入383 （

9、导热系数），单击“ OK”按钮； 再次拾取如图 33-5 所示对话框中右侧列表的“ Specific Heat ”项,弹出如图 33-7 所示的 对话框，在“ C”文本框中输入 390 （比热），单击“ 0K按钮；再次拾取如图33-5所示对话框中右侧列表的“ Density ”项，弹出如图33-8所示的对话框，在“ DENS文本框中输入 8889 （密度），单击“ 0K按钮，于是定义好了材料模型1 （铜）。单击如图 33-5 所示对话框的菜单项Material-New Model,单击弹出的“ DefineMaterial TD ”对话框中的“ 0K按钮，然后重复定义材料模型1时的各步骤，定义

10、材料模型 2（铁）的导热系数为 70,比热为 448,密度为 7833。重复定义材料模型 2 时的各步骤, 定义材料模型 3（水）的导热系数为 2,比热为 4185,密度为 996。最后关闭如图 33-5 所示 的对话框。略图 33-5 材料模型对话框图 33-6 定义导热系数对话框图 33-7 定义比热对话框拾取菜单 Main Menu 宀 Preprocessor 宀 Modeling Create Areas Rectangle TBy Dimension ,弹出如图33-9所示的对话框，在“X1, X2文本框中分别输入 0,，在“Y1, Y2文本框中分别输入 0,单击“Apply ”按

11、钮，再次弹出如图33-9所示的对话框，在“Xl, X2文本框中分别输入， 在“ Y1, Y2”文本框中分别输入，单击“ Apply ”按钮，再 次弹出如图33-9所示的对话框，在“ Xl , X2”文本框中分别输入，在“ Y1, Y2”文本框中 分别输入，单击“ 0K按钮。拾取菜单 UtilityMenu PlotCtrls 宀 Numbering，在弹出的Plot Numbering Controls对话框中，将 Area numbers （面号）打开，单击“ OK按钮。图 33-8 定义密度对话框图 33-9 创建矩形对话框拾取菜单 Main Menu 宀 Preprocessor 宀 M

12、odeling 宀 Operate 宀 BooleansOverlap 宀Areas, 弹出拾取窗口，单击“ Pick All ”按钮。拾取菜单 Main MenuPreprocessorMeshingMeshTool，弹出如图 33-10 所示的对话框。选择“ Element Attributes 的下拉列表框为“Areas ”，单击下拉列表框后面的“Set”按钮，弹出拾取窗口，选择面2，单击拾取窗口中的 “OK按钮，弹出.” Areas Attributes 对话框，选择“ MAT下拉列表框为1,单击“Apply ”按钮，再次弹出拾取窗口，选择面3,单击拾取窗口中的“ OK按钮，选择“ A

13、reas Attributes对话框的“ MAT下拉列表框为2,单击Apply 按钮，再次窗口，择面4，单击拾取窗口中的“OK按钮，选择“Areas Attributes 对话框的“ MAT下拉列表框为3，单击OK按钮。单击“ Size Controls 区域中“ Global 后面的“ Set按钮，弹出如图 33-11所示的对话框，在“SIZE”文本中输入，单击“OK按钮；在如图33-10所示对话框的Mesh区域， 选择单元形状为“Quad”（四边形），选择划分单元的方法为 “ Map ped（映射），单击“ Mesh” 按钮，弹出拾取窗口，拾取面 2和面3,单击“ OK按钮。图 33-10

14、划分单元工具对话框图 33-11单元尺寸对话框单击“ Size Controls ”区域中“ Global ”后面的“ Set”按钮，弹出如图 33-11所示的 对话框，在” SIZE文本框中输入，单击“ OK按钮；在如图33-10所示对话框的“ Mesh” 区域，选择单元形状为Quad （四边形）， 选择划分单元的方法为“ Free ”（自由），单击“Mesh”按钮，弹出拾取窗口，拾取面4,单击“OK按钮，关闭如图33-10所示的对话框。拾取菜单 Main Menu宀Solution 宀Analysis Type New Analysis，在弹出的对话框中 选择“Type of Analys

15、is 为Transient ”，单击OK按钮，在随后弹出的 “Transient Analysis 对话框中，单击“ OK按钮。以下步骤进行稳态分析,以得到初始温度场。拾取菜单 Main Menu SolutionLoad Step Opts Time/FrequencTime IntegrationTAmplitude Decay,弹出如图33-12所示的对话框，将“ TIMINT”关闭，单击“ OK按钮。图 33-12 时间积分控制对话框拾取菜单 Main Menu Solution Load Step Opts Time/Frequenc TimeTime Step, 弹出如图33-13

16、所示的对话框，在“ TIME”文本框中输入，在“ DELTIM Time step size ” 文本框中输入，单击“ OK按钮。图 33-13 确定载荷步时间和时间步长对话框拾取菜单 Main Menu 宀 Solution 宀 Define Loads Apply Thermal Temperature OnAreas，弹出拾取窗口，拾取面2,单击“OK按钮，弹出如图33-14所示的对话框，在“Lab2” 列表中选择“ TEMP，在“ VALUE文本框中输入 400, 单击“ Apply ”按钮。再次弹出拾 取窗口，拾取面 3,单击“ OK按钮，弹出如图 33-14所示的对话框，在“ VA

17、LUE文本框 中输入500,单击“ Apply ”按钮。再次弹出拾取窗口，拾取面4,单击“ OK按钮，弹出如图33-14所示的对话框，在“ VALUE文本框中输入20，单击OK “按钮。拾取菜单 Main Menu Solution Solve Current LS,单击 “ Solve Current Load Step 对话框中的OK”按钮。图 33-14 在面上施加温度载荷对话框以下步骤进行瞬态分析。拾取菜单 Main Menu Solution Load Step Opts Time/Frequenc TimeTime Step , 弹出如图33-13所示的对话框，在“TIME文本框中

18、输入 3600,在“DELTIMTime step size ” 文本框中输入 50,选择“ AUTOTS 为“ OM,在“ DELTIM Minimum time step size ”文本 框中输入10,在“ DELTIMMaximumtime step size ”文本框中输入 200，单击“ OK”按钮。拾取菜单 Main Menu 宀 Solution 宀 Load Step Opts Time/Frequenc Time Integration-Amplitude Decay，弹出如图 33-12 所示的对话框，将“ TIMJNT打开，单击“ OK按钮。拾取菜单 Main Menu

19、 宀 Solution 宀 Define Loads Delete Thermal Temperature TOn Areas，弹出拾取窗口，单击“ Pick All ”按钮。拾取菜单 Main Menir Solution T Load Step OptsT Output Ctrls T DB/Results File , 弹出“ Controls for Database and Results File Writing ”对话框, 选择下拉列表框 “ Item ” 为“ All Items ”，选中“ Every substep ”，单击“ OK” 按钮。菜单 Main MenutSol

20、ution tSolve tCurrent LS ，单击” Solve Current Load Step ”对话框中的“ OK按钮。以下步骤查看结果。拾取菜单 Main Menut General Postproc t Plot Results tContour Plot tNodal Solu ， 在列表中依次选择Nodal Solution 宀 DOFSolution 宀 Nodal Temperature ”，单击OK按 钮，结果如图 33-15 所示。图 3600 s 后的温度分布情况拾取菜 Main MenuTimeHist PostproDefine Variables,弹出 “

21、Defined TimeHistory Variables 对话框，单击Add”按钮，弹出Add Time-History Variables ” 对话框，选择Type of Variable ”为Nodal DOFresult ”，单击OK按钮，弹出拾取窗口，在拾取窗口的文本窗口中输入106 （位于铜块中心），单击“ OK按钮，弹出如图33-16所示的对话框，在“ Name文本框中输入 TEMP_106单击“ OK按钮。关闭“ Defined Time-History Variables对话框。图 33-16 定义节点数据对话框拾取菜单 Main Menu TimeHist Postpro

22、Graph Variables，弹出Graph Time-History Variables ”对话框，在“NVAR1文本框中输入2,单击“OK按钮，结果如图33-17所示。图 33-17 节点 106 处的温度变化情况第 34 例 在结构上直接施加温度载荷进行热应力分析实例 双金属簧片本例介绍了在结构单元上直接施加温度载荷进行热应力分析的方法和步骤，并使用解析解对有限元分析结果进行了验证。概述当一个结构发生温度变化时， 它会发生膨胀或收缩。 如果其膨胀或收缩受到限制， 或者 由于温度分布不均而导致膨胀和收缩的程度不同，就会产生热应力。ANSYS提供了三种计算热应力的方法，即在结构上直接施加温

23、度载荷、间接法和直接法。 本例介绍第一种方法， 该方法主要用于结构温度场为已知且均匀分布的场合。 分析类型为静 力学分析，温度场的施加采用 Main Menu宀Solution 宀 Define Loads Apply Structural TTemperature命令。此时，节点温度作为静力学分析的体载荷而不是节点自由度。 问题描述及解析解温度调节器所使用的双金属簧片的计算简图如图 34-1 所示。两簧片等厚度、等宽度， 材料不同，二者紧密连成一体，其左端可视为固定端。第1 个金属片为钢制，其弹性模量E1=2x1011Pa,线膨胀系数，a 1=10 x10-6/ C；第2个金属片为铜制，其弹

24、性模量E2=x1011Pa,线膨胀系数a2=16xl0-7 C。簧片各部分尺寸为长度L=40mm高度h=。试分析当温度升高 t=IOO c时，簧片自由端的挠度。图34-1双金属簧片根据材料力学的知识，容易得出簧片自由端的挠度为?= 6 （a2- a12?tElE2L2 =X 10-4m h （ E12 + E 12+14E 1 E2）分析步骤拾取菜单 Utility Menu File Cha nge Job name,在所弹出对话框的“ /FILNAM ” 文本框中输入 EXAMPLE34单击“ 0K按钮。拾取菜单 Ma in Menu PreprocessorEleme nt Type A

25、dd/Edit/Delete ，弹出如图 34-2所示的对话框，单击“Add”按钮， 弹出如图34-3所示的对话框，在左侧列表中选“Structural Solid ”，在右侧列表中选“ 20node 186”，单击“ OK按钮，返回到如图 34-2 所示的对话框，单击“ Close ”按钮。图34-2 单元类型对话框图34-3 单元类型库对话框拾取菜单 Main Menu宀Preprocessor 宀Material Props Material Models ，弹出如图34-4 所示的对话框，在右侧列表中依次拾取“ Structural ”、“ Lin ear ”、“ Elastic ”、

26、Isotropic ”，弹出如图34-5所示的对话框，在“ EX文本框中输入 2e11（弹性模量），在“ PRXY文本框中输入（泊松比），单击“ OK按钮；再在图34-4所示对 话框的右侧列表 中依次拾取，Structural 下面的“ Thenual Expansion ”、“ Secant Coefficient ”、“Isotropic ”，弹出如图34-6所示的对话框，在“ALPX文本框中输入10e-6（线膨胀系数），单击“OK按钮。单击如图34-4所示对话框的菜单项 Material宀NewModel, 单击所弹出“ Define Material ID”对话框中的“ OK按钮，然后

27、重复定义材料模型 1时的 各步骤，定义材料模型2的弹性模量为，泊松比为，线膨胀系数为16e-6，最后关闭如图34-4 所示的对话框。图34-4 材料模型对话框图34-5 材料特性对话框（一）图34-6 材料特性对话框（二）拾取菜单Main Menu 宀 Preprocessor 宀 Modeling Create Volumes Block ByDimension，弹出如图34-7所示的对话框，在“ X1, X2”文本框中分别输入0,在“ Y1,Y2文本框中分别输入 0, 5，在“ 21 , 22”文本框中分别输入 0,单击“ Apply ”按钮。再次弹出如图34-7所示的对话框，在“ X1,

28、 X2文本框中分别输入0,在“ Y1, Y2”文本框中分别输入，在“ Z1，22”文本框中分别输入 0，单击“ OK按钮。图 34-7 创建块对话框5 黏结体拾取菜单 Main MenuPreprocessor Modeling OperateBooleansGlueVolumes. 弹出拾取窗口,单击“ Pick All ”按钮。拾取菜单 Main Menu Pteprocessor Meshing MeshTool，弹出如图 34-8 所示的对话 框。选择“ Element Attributes ”的下拉列表框为“ Volumes”，单击列表框后面的“ Set ” 按钮，弹出拾取窗口.选择

29、上面的六面体，单击拾取窗口中的“OK按钮，弹出“ VolumeAttributes ”对话框，选择“ MAT下拉列表框为1,单击“ Apply ”按钮，再次弹出拾取窗 口，选择下面的六面体，单击拾取窗口中的“OK按钮，选择“ Volume Attributes ”对话框的“ MAT下拉列表框为 2,单击“ OK按钮。单击“ Size Controls ”区域中“ Global 后面的“ Set”按钮，弹出如图 34-9所示的 对话框，在“ SIZE”文本框中输入，单击“ OK按钮。图 34-8 划分单元工具对话框图 34-9 单元尺寸对话框在图34-8所示对话框的“ Mesh”区域，选择单元形

30、状为“ Hex”（六面体），选择划分单 元的方法为“Mapped（映射），单击“Mesh”按钮，弹出拾取窗口，单击“Pick All ”按钮， 最后单击如图34-8所示对话框中的“ Close ”按钮。拾取菜单 Utility Menu宀PlotCtrls 宀Pan Zoom Rotate，在弹出的对话框中，依次单击“ Iso ”、“ Fit ”按钮，或者单击图形窗口右侧显示控制工具条上的剧按钮。拾取菜单 Main Menu 宀 Solution 宀 Define Loads Apply StructuralDisplacementTOn Areas，弹出拾取窗口，拾取 x=0平面，单击“ O

31、K按钮，弹出如图 34-10所示的对话 框，在列表中选择“ All DOF ”，单击“ OK按钮。图 34-10 在面上施加约束对话框拾取菜单 Main MenutSolution tDefine Loads tApplytStructural tTemperature t On Volumes,弹出拾取窗口，单击“ Pick All ”按钮，弹出如图34-11所示的对话框，在“ VAL1” 文本框中输入100,单击“ OK按钮。Y2文本框中分别输入 0, 5，在“ 21 , 22”文本框中分别输入 0,单击“ Apply ”按钮。图 34-11 在体上施加温度载荷对话框拾取菜单 Ma in

32、Menu Solution Define LoadsSetti ngs-Refere nee Temp，弹出如图34-12所示的对话框，在“ TREF文本框中输入0，单击“ OK按钮。图34-12 定义参考温度对话框拾取菜单 Main Menu SolutionSolve Current LS,单击Solve Current Load Step对话框中的“OK按钮，当出现“Solution is done! ”提示时，求解结束，即可查看结果。拾取菜单 Main Menu 宀 General Postproc Plot Results Deformed Shape,弹出如图34-13所示的对话框

33、，选中“ Def+undef edge ”(变形+未变形的模型边界)，单击“ OK按 钮，结果如图34-14所示。图34-13 显示变形对话框图34-14 簧片的变形可见，簧片自由端的位移即簧片的最大位移为x 10-4mm，与理论解基本一致，从而，证明有限元解具有较高的精度第35例基于实测温度场进行热应力分析实例-转炉托圈本例通过实例介绍了利用 ANSYS参数化设计语言(APDL)编制程序，将实测温度场施加 到有限元模型上，并对结构的热应力进行有限元分析，得到了结构各部分热应力的精确分布。 概述托圈是转炉的重要传动和承载部件，工作过程中的热物理过程十分复杂，其温度场分布通过解析法或数值法进行计

34、算几乎是不可能的，通常托圈的温度场只能是通过现场测试得到，但由于托圈结构比较复杂，其热应力计算一般都是用有限元方法进行的，因此存在着如何将实测得到的温度场施加到有限元模型上的问题。某转炉托圈形状可近似为如图35-l (a)所示的回转体，图35-1(b)为其横截面形状，已知r0=、r1=、h=,经测量托圈在 0、90。、180。、270。四个横截面处 A、B C D四点的温 度如表35-1所示，沿托圈圆周和横截面高度方向温度近似线性分布，现计算其热应力。图35-1转炉托圈的横截面表35-1 托圈温度分布横截面位置A点BCD01418812280MSHKEY,1映射网格90171115135891

35、801389012379270212142163135由于托圈的温度按线性规律分布，所以有限元分析时可以由0、90、180、270四个横截面的温度值通过线性插值得到其他截面上A、B C、D四点的温度值，再由该截面上A B、C D四点的温度值线性插值得到截面上其他点的温度值。于是，由若干个测量点的 温度值插值得到整个托圈上所有点的温度值。使用ANSYS寸托圈热应力进行分析时，可以利用其参数化设计语言（APDL）编制程序来模拟上述线性插值过程，然后把所计算得到的各个节点的温度值作为载荷施加到有限元模型上。程序框图如图35-2所示。命令流/CLEAR!清除数据库，新建分析/FILNAME,EXAMP

36、LE35!定义任务名为” EXAMPLE35 “R0=!定义参数，尺寸R1 =H=TH=/PREP7!ET, 1, PLANE182!ET, 2, SOLID185MP,EX,1,2E11MP,PRXY,1,!MP,ALPX,1,!RECTING,R0,R0+TH,-H/2,H/2RECTING,R1-TH,R1,-H/2,H/2RECTING,R0,R1,H/2-TH,H/2RECTING,R0,R1,-H/2+TH,-H/2AOVLAP,ALL!MSHAPE,0!进入预处理器选择单元类型!定义材料模型，弹性模量泊松比线膨胀系数!创建矩形面MSHKEY,1映射网格交叠所有面六面体单元ESIZ

37、E,单元边长度挤出段数旋转挤出FINI/SOLUNSEL,S,LOC,X,0在选择的节点上施加位移约束K,101,0,1EXTOPT, ESIZE,10VROTAT,ALL,100,101,360退出预处理器进入求解器D,ALL,UXNSEL,S,LOC,Y,0D,ALL,UYNSEL,S,LOC,Z,0D,ALL,UZALLS ! 选择所有TA_0=141 ! 定义参数，实测温度 TB_0=88TC_0=122TD_0=80TA_90=171TB_90=115TC_90=135TD_90=89TA_180=138TB_180=90TC_180=123TD_180=79TA_270=212TB_270=142AMESH,ALL对面划分单元K,100创建关键点EXTOP