ANSYS瞬态热分析教程及实例

1、QUST 精品课件 1 n4. 瞬态传热分析 QUST 精品课件 2 4. 瞬态传热分析 n瞬态传热分析的基本步骤与稳态热分析类似，主要的 区别是瞬态传热分析中的载荷是随时间变化的。为了 表达随时间变化的载荷，首先必须将载荷时间曲线 分为载荷步。载荷时间曲线中的每一个拐点为一个 载荷步，如图3-9 所示。对于每一个载荷步，必须定 义载荷值及时间值，同时必须选择载荷步为Ramped 方式变化或Stepped 方式变化。 n如果需要知道系统受随时间变化(或不变)的载荷和边 界条件时的响应，就需要进行“瞬态分析” 。 QUST 精品课件 3 4. 瞬态传热分析 QUST 精品课件 4 5. 瞬态传热

2、分析 nANSYS 缺省是渐进加载的。渐进加载可以提 高瞬态求解的适应性，如果有非线性时可以提 高收敛性。 QUST 精品课件 5 5. 瞬态传热分析 QUST 精品课件 6 5. 瞬态传热分析 n在瞬态分析中，载荷步和子步的定义与非线性稳态分 析十分类似。载荷定义的每个载荷步的终点，并可以 随时间阶跃或渐进的施加。 n每个载荷步的求解是在子步上得到。子步长根据时间 积分步长得到。 n 自动时间步 (ATS) 同样适用于瞬态分析，可以简化 ITS选择。 nITS选择将影响到瞬态分析的精度和非线性收敛性 (如果存在)。 QUST 精品课件 7 均匀初始温度：如果整个模如果整个模 型的初始温度为均

3、匀且非型的初始温度为均匀且非0 ，使用下列菜单指定，使用下列菜单指定: 2 1 3 4 5. 瞬态传热分析 QUST 精品课件 8 非均匀的初始温度 如果模型的初始温度分布如果模型的初始温度分布 已知已知但但不均匀不均匀，使用这些使用这些 菜单将初始条件施加在特菜单将初始条件施加在特 定节点上定节点上 2 1 3 5 4 5. 瞬态传热分析 QUST 精品课件 9 5. 瞬态传热分析 nANSYS 瞬态传热分析的主要步骤 n1.建立有限元模型 n2.施加载荷并求解 n3.求解 n4.查看分析结果 QUST 精品课件 10 n择分析类型 n进行瞬态传热分析需要首先需要定义分析类型及其相 关选项。

4、下面介绍分析类型及其选项的设定： n进行第一次分析或者重新进行分析 n命令：ANTYPE,TRANSIENT,NEW nGUI：Main Menu Solution Analysis Type New Analysis Transient n延续上一次分析 n命令：ANTYPE,TRANSIENT,REST nGUI：Main Menu Solution Analysis Type Restart QUST 精品课件 11 n定义瞬态传热分析的初始条件 n瞬态传热分析的初始条件分为两种情况：其一， 初始温度场已知；其二，初始温度场未知。 QUST 精品课件 12 n已知初始温度场已知初始温度场

5、 n如果初始温度场是已知的，则定义过程比较简单，定 义过程如下： n1. 定义均匀温度场 n命令：TUNIF nGUI：Main Menu Solution Define Loads Apply Structural Temperature Uniform Temp n如果已知模型的起始温度是均匀的，可设定所有节点 初始温度。 QUST 精品课件 13 n2. 设定参考温度 n命令：TREF nGUI：Main Menu Solution Define Loads Settings Reference Temp n如果不在对话框中键入数据，则默认为参考温 度，参考温度的值默认为零。 QUST

6、精品课件 14 n3. 设置节点温度 n命令：D nGUI：Main Menu Solution Define Loads Apply Thermal Temperature On Nodes QUST 精品课件 15 n3. 设定非均匀的初始温度 n命令：IC nGUI：Main Menu Solution Define Loads Apply Initial Conditn Define n在瞬态传热分析中，节点温度可以通过此项设 定为不同的值。 QUST 精品课件 16 n初始温度场未知初始温度场未知 n如果初始温度场是不均匀的且又是未知的，就 必须首先作稳态热分析确定初始条件，下面介

7、绍分析选项的设定： n1）施加载荷（如已知的温度、热对流等） n2）关闭时间积分 n命令：TIMINT, OFF nGUI：Main Menu Solution Load Step Opts Time/Frequenc Time Integration QUST 精品课件 17 n设定一个只有一个子步的，时间很小的载荷步（例如 0.001） n命令：TIME nGUI：Main Menu Solution Analysis Type Soln Controls basic n4）写入载荷步文件 n命令：LSWRITE nGUI：Main Menu Solution Load Step Opts

8、 Write LS File n或先求解： nGUI：Main Menu Solution Solve Current LS QUST 精品课件 18 n非线性选项常用选项如下： n）迭代次数选项选项 n命令：NEQIT nGUI：Main Menu Preprocessor Loads Load Step Opts Nonlinear Equilibrium Iter n每个子步默认的次数为25，这对大多数非线性 热分析已经足够。 QUST 精品课件 19 n自动时间步长选项 n命令：AUTOTS nGUI：Main Menu Solution Analysis Type Soln Cont

9、rols n打开后求解过程中将自动调整时间步长。 QUST 精品课件 20 n时间积分选项 n命令：TIMINT nGUI：Main Menu Solution Load Step Opts Time/Frequenc Time Integration n如果将此选项设定为OFF，将进行稳态热分析。 QUST 精品课件 21 n求解 n命令：SOLVE nGUI：Main Menu Solution Current LS QUST 精品课件 22 POST1 后处理 n读入结果数据 n命令：SET nGUI：Main Menu General Postproc Read Results By

10、Time/Freq n进入POST1 后，可以读出某一时间点的结果。如果 设定的时间点不在任何一个子步的时间点上， ANSYS 会进行线性插值。 n此外还可以读出某一载荷步的结果： nGUI：Main Menu General Postproc Read Results By Load Step QUST 精品课件 23 n3.3 瞬态传热分析实例 1. 问题描述 k实例类型：ANSYS结构分析。 k分析类型：瞬态传热分析。 k单元类型：PLANE55 kANSYS功能示例：实体建模包括基本的建模操作； 定义比热容；施加瞬态热载荷；设置瞬态热载荷 分析选项；显示模型温度等值线图；显示节点温 度

11、随时间变化曲线。 QUST 精品课件 24 QUST 精品课件 25 n长方形的板，几何参数及其边界条件如图3-6 所示。板的宽度为5cm，其中间有一个半径为 1cm 的圆孔。板的初始温度为20，将其右 侧突然置于温度为20且对流换热系数为 100W/M2的流体中，左端置于温度为500 的温度场，试计算： n（1）第1s 和第50s板内的温度分布情况。 n（2）整个板在前50s内的温度变化过程。 n（3）圆孔边缘A点处温度随时间变化曲线。 QUST 精品课件 26 n2.建立有限元模型 k首先建立瞬态传热分析所需的有限元模型 k 选择单元. (1) 选择热分析单元，操作如下：选择热分析单元，操

12、作如下： kGUI：Main Menu Preprocessor Element Type Add/Edit/Delete k 在弹出的对话框中，单击Add。在单元类型 库对话框中选择Plane55单元。单击OK。 k命令：ET,1,PLANE55 QUST 精品课件 27 n(2) 定义材料属性 k首先进入Define Material Model Behavior 对话框，操作如下： kGUI：Main Menu Preprocessor Material Props 下面定义瞬态热分析所需的材料参数， 如热传导率、比热容及材料密度： QUST 精品课件 28 n定义热传导率 k GUI：

13、Main Menu Preprocessor Material Props Thermal Conductivity Isotropic k 在弹出的定义材料热传导率对话框中的KXX 栏键 入“5”。 k 命令：MPDATA,KXX,1,5 QUST 精品课件 29 n定义比热容 k GUI：Main Menu Preprocessor Material Props Thermal Specific Heat k 在弹出的定义比热容对话框中的C栏键入 “200”。 k 命令：MPDATA,C,1,200 QUST 精品课件 30 n定义密度 k GUI：Main Menu Preprocess

14、or Material Props Thermal Density k 在弹出密度定义对话框中的DENS栏键入 “5000”。 k 命令：MPDATA,DENS,1,5000 材料属性定义完毕. QUST 精品课件 31 n(3)建立实体模型 k 根据本例所用模型，首先需要创建矩形，然后 是圆，最后在矩形板中央减去（Substract）圆。 下面介绍建立实体模型的操作： k 创建矩形 k 命令：RECTNG,0,0.15,0,0.05, QUST 精品课件 32 n创建圆面 k 其操作如下： k GUI：Main Menu Preprocessor Modeling Create Circle

15、 By dimensions k 在弹出对话框中，单击OK得到圆面。 k 命令：CYL4, 0.075, 0.025, 0.01 QUST 精品课件 33 n相减 k根据ANSYS建模过程中面序号赋值原理，直 接可以肯定圆面序号为2，矩形序号为1，因 此采用直接键入命令建实体模型： k命令：asba,1,2 QUST 精品课件 34 n(4) 设定网格尺寸并划分网格 下面介绍网格尺寸的设定（SmartSize 方 式）： k设定网格尺寸参数并划分网格,通过 SmartSize控制网格密度，操作如下： kGUI：Main Menu Preprocessor Meshing MeshTool k选

16、择SmartSize=3。单击Mesh。单击拾取 对话框中Pick All按钮。得到网格图。 QUST 精品课件 35 k命令：SMRT, 3 kAMESH,All k保存数据库,其操作如下： kGUI：Toolbar SAVE-DB k命令：SAVE QUST 精品课件 36 n4.施加载荷并求解 k 求解之前首先要选择分析类型，然后定义 边界条件及其载荷步选项，最后计算。首先 选择分析类型。 (1) 选择分析类型 k 选择Transient分析，操作如下： QUST 精品课件 37 kGUI：Main Menu Preprocessor Loads Analysis Type New An

17、alysis k选择Transient 分析，单击OK。采用ANSYS 默认设置，在弹出的子对话框中单击OK。 k命令：ANTYPE,4 kTRNOPT,FULL kLUMPM,0 QUST 精品课件 38 n(2)定义初始条件 k板的初始温度为20，设置初始温度操作 如下： kGUI：Main Menu Solution Define Loads Apply Initial Conditn Define k在弹出的拾取对话框中，单击Pick All。弹 出Define Initial Conditions对话框， k命令：IC,All,TEMP,20, QUST 精品课件 39 n(3)定义

18、热约束 k 瞬态传热分析中的载荷是随着时间发生 变化的。对于每一个载荷步都需要指明载 荷值及时间值， k还需要指定载荷步选项，如加载方式是 Ramped方式还是Stepped方式。 QUST 精品课件 40 n需要施加流体载荷和板的传热载荷。首 先定义对流载荷。 定义对流边界 k为了便于设定对流边界，首先显示边框图， 操作如下： kGUI：Utility Menu Plot lines k命令：LPLOT QUST 精品课件 41 k定义对流载荷/边界首先进入Apply Conv on lines 对话框，操作如下： kGUI：Main Menu Solution Define Loads A

19、pply Thermal Convection On Lines k在弹出对话框中，键入Film coefficient和Bulk Temperature值。 k命令：SFL,L2,CONV,100, ,20, QUST 精品课件 42 n定义稳态热边界 k在边线上定义稳态热边界，操作如下： kGUI：Main Menu Solution Define Loads Apply Thermal Temperature On Lines k在弹出对话框中，键入边界温度为“500”。 单击OK。 k命令：DL,L4, ,TEMP,500,1 QUST 精品课件 43 n(4) 设置时间及时间步进参数

20、 n设定步进参数使ANSYS 在较短的时间内达到可以接受的收 敛精度。 k首先进入Time and time Step Optios 对话框，操作如下： kGUI：Main Menu Solution Load Step Opts Time/Frequenc Time time steps k在Time and time Step Optios对话框中按照下图 设置，其余选项不 变。 kTIME,50 kAUTOTS,1 kDELTIM,1,0.1,2.5,1 kKBC,0 QUST 精品课件 44 n下面确认时间积分选项打开。操作如下： kGUI：Main Menu Solution Loa

21、d Step Opts Time/Frequenc Time Integration k保持默认设置，单击OK。 k命令：TIMINT,1 kTINTP,0.005, , ,-1,0.5,-1 QUST 精品课件 45 n(5) 设置输出控制 n输出控制用于决定ANSYS 要等多长时间向数 据库文件写入一次结果。这里设时间间隔为 1s。输出控制选项设置，操作如下： kGUI：Main Menu Solution Load Step Opts Output Ctrls DB/Results File k在File Write Frequency栏选项 Every Substep。 k命令：OUT

22、RES,ALL,ALL, QUST 精品课件 46 n(6) 求解 kGUI：Main Menu Solution Solve Current LS. k命令：Solve QUST 精品课件 47 n5.查看分析结果 k根据题目要求，首先在POST1 后处理中查 看第1s 和50s 时温度分布，然后动画显示温 度场变化情况， k最后是通过POST26 查看A点处温度随时间 变化曲线。 n下面解释查看分析结果的详细操作： QUST 精品课件 48 n(1) POST1后处理 k显示第1s和和第50s时温度场 k首先读入分析结果文件，然后分别显示第 1s和第50s时模型温度场。 首先读入第1s时结

23、果文件，操作如下： QUST 精品课件 49 k GUI：Main Menu General Postproc By Time/Freq k 在弹出的Read Results by time or Frequency对话框中进行如下图设置。单击 OK。显示第1s 时温度场，操作如下： k GUI：Main Menu General Postproc Plot Results Nodal Solu k 取ANSYS 默认设置，单击OK。显示第1s 时温度场。 k 命令：PLNSOL,TEMP, ,0, QUST 精品课件 50 QUST 精品课件 51 n重复上述操作，读入第50s时结果并显示温

24、度 场。下图所示为第50s温度场。 QUST 精品课件 52 n显示温度场变化动画 n显示150 秒内温度变化情况，操作如下： nGUI：Utility Menu PlotCtrls Animate Animate Over Time n设置显示150s内温度变化动画。 n命令：ANTIME,20,0.5, ,1,2,1,50 n动画如下transient1.avi QUST 精品课件 53 n(2) POST26 后处理 k 需要查看A点随时间变化的温度场，首 先必须定义A点温度变量，然后才能显示A 点处 k 温度随时间变化曲线。下面介绍查看A 点分析结果的详细操作： k定义A点温度变量 k

25、GUI：Main Menu TimeHist Postpro Define Variables QUST 精品课件 54 k单击Add按钮，弹出Add Time-History Variable对话框，选择温度作为变量）。选 择一点A k命令：NSOL,2,7,TEMP, ,TEMP-2 k显示A点温度-时间变化曲线 k显示A 点温度随时间变化曲线，单击 “GRAPH VARIABLE 按钮，得到A 点温度 随时间变化曲线。 k命令：PLVAR,TEMP-2 QUST 精品课件 55 n3.3.5 3.3.5 瞬态传热分析实例瞬态传热分析实例2 2 n(1)(1)问题描述：问题描述： k 一个

26、直径为一个直径为0.12m，温度为，温度为900摄氏度的钢球突摄氏度的钢球突 然放入盛满了水的，完全绝热的横截面直径和孤然放入盛满了水的，完全绝热的横截面直径和孤 傲度均为傲度均为0.6m的圆柱体水箱中，水的温度为的圆柱体水箱中，水的温度为20， 材料参数如图所示，求解材料参数如图所示，求解10分钟后钢球与水的温分钟后钢球与水的温 度场分布（钢球置于水箱正中央）。度场分布（钢球置于水箱正中央）。 QUST 精品课件 56 材料密度 kg/m3 导热系数K W/m 比热C J/kg. 水10000.614185 钢780070448 材料性能参数 QUST 精品课件 57 n2. 问题分析 k

27、属于瞬间热传导问题，研究对象为钢 球和水，根据轴对称性，在求解过程中 取钢球和水中心纵截面的1/4建立几何模 型，选择PLANE55轴对称单元进行分析。 几何模型 QUST 精品课件 58 n3.建立模型 (1)添加标题 k添加标题，操作如下GUI：Utility Menu File Change Title k标题名：Transient thermal analysis to a sphere. k命令：/Title,Transient slab Problem QUST 精品课件 59 简化菜单（过虑菜单）操作如下： kGUI：Main Menu Referenc k在弹出的Referen

28、c for GUI Filtering对话框中，选 择 Thermal。单击OK。 k命令：/COM,Thermal QUST 精品课件 60 (2) 选择单元 k选择热分析单元，操作如下： kGUI：Main Menu Preprocessor Element Type Add/Edit/Delete k在弹出的对话框中，单击Add。在单元类型库对 话框中选择Plane55单元。单击OK。 k命令：ET,1,PLANE55 QUST 精品课件 61 n单击Element Type 对话框中的Options按钮， 出现PLANE55 Element Type options 对话框， 在Ele

29、ment behavior K3下拉框中选择 Axisymmetric 选项，其余默认。单击“OK”。 QUST 精品课件 62 n(3)定义材料属性1 k首先进入Define Material Model Behavior 对话框，操作如下： kGUI：Main Menu Preprocessor Material Props Material Models k下面定义瞬态热分析所需的材料参数，如 热传导率、比热容及材料密度： QUST 精品课件 63 n定义热传导率 kGUI：Main Menu Preprocessor Material Props Thermal Conductivit

30、y Isotropic k在弹出的定义材料热传导率对话框中的KXX 栏键入“70”。 k命令：MPDATA,KXX,1,70 QUST 精品课件 64 n定义比热容 kGUI：Main Menu Preprocessor Material Props Thermal Specific Heat k在弹出的定义比热容对话框中的C栏键入 “448”。 k命令：MPDATA,C,1,448 QUST 精品课件 65 n定义密度 kGUI：Main Menu Preprocessor Material Props Thermal Density k在弹出密度定义对话框中的DENS栏键入 “7800”。

31、 k命令：MPDATA,DENS,1,7800 k材料属性定义完毕。 QUST 精品课件 66 n(4)定义材料属性2 k 在首先进入Define Material Model Behavior 对话框，选择MaterialNew Model，出现Define Material ID对话框。在文本框中输入材料参考号 2，单击“OK”。 QUST 精品课件 67 n定义热传导率 kGUI：Main Menu Preprocessor Material Props Thermal Conductivity Isotropic k在弹出的定义材料热传导率对话框中的KXX 栏键 入“0.61”。 k命

32、令：MPDATA,KXX,1,0.61 QUST 精品课件 68 n定义比热容 kGUI：Main Menu Preprocessor Material Props Thermal Specific Heat k在弹出的定义比热容对话框中的C栏键入“4185”。 k命令：MPDATA,C,1,4185 QUST 精品课件 69 n定义密度 kGUI：Main Menu Preprocessor Material Props Thermal Density k在弹出密度定义对话框中的DENS栏键入“1000”。 k命令：MPDATA,DENS,1, 1000 k材料属性定义完毕。 QUST 精品

33、课件 70 n(5)建立实体模型 kGUI：Preprocessor Modeling Create Areas Rectangle By Dimensions k创建矩形 k命令：RECTNG,0, 0.3, 0, 0.3, QUST 精品课件 71 n创建圆面 n其操作如下： kGUI：Main Menu Preprocessor Modeling Create Arcs By Dimensions k在弹出对话框中，单击OK得到圆面。 k命令：PCIRC, 0.06, ,0,90, QUST 精品课件 72 n合并 kGUI：Main Menu Preprocessor Modeling

34、 Operate Booleans overlapAreas k出现overlap Areas 对话框，单击Pick All 关闭对话 框。 QUST 精品课件 73 n压缩编号 kGUI：Main Menu Preprocessor Numbering Ctrls Compress Numbers，在 Compress Numbers对话框，在Label Item to be compressed 下拉列表中选择All 选择， 单击“OK”关闭对话框。 QUST 精品课件 74 n4.设定网格尺寸并划分网格 n设定网格尺寸参数，操作如下： kGUI：Main Menu Preprocesso

35、r Meshing Size Cntrls ManualSize LinesPicked Lines k在 Element Size On 菜单中选择线4、5 k在Element Size On Picked Lines 对话框中， 在NDIV 文本框中输入“30”，在SPACE Spacing ratio文本框中输入“0.1”，单击OK QUST 精品课件 75 k在 Element Size On 菜单中选择线6、7 k在Element Size On Picked Lines 对话框中， 在NDIV 文本框中输入“32”，在SPACE Spacing ratio文本框中输入“0.1”，单

36、击OK k在 Element Size On 菜单中选择线3 k在Element Size On Picked Lines 对话框中， 在NDIV 文本框中输入“30”， QUST 精品课件 76 n保存数据库,其操作如下： kGUI：Toolbar SAVE-DB k命令：SAVE QUST 精品课件 77 n划分网格，操作如下： kGUI：Main Menu Preprocessor Meshing Mesh Volumes Mapped ConcatenateLines， 出现Concatenate Lines菜单，在文本框中输入2， 1，单击OK kGUI：Main Menu Prep

37、rocessor Meshing Mesh attributes Default attribs，出现Meshing attributes 菜单，在MAT下拉列表框中选择1，单 击OK kGUI：Main Menu Preprocessor Meshing Mesh Tool，出现Mesh Tool菜单，在shape选项 组中选中Quad和Mapped选项，单击Mesh，出现 Mesh Areas，在文本框中输入1，单击OK。 QUST 精品课件 78 kGUI：Utility Menu Select Everything kGUI：Main Menu Preprocessor Meshing

38、Mesh attributes Default attribs， 出现Meshing attributes 菜单，在MAT下拉 列表框中选择2，单击OK. kGUI：Main Menu Preprocessor MeshingMesh Tool，出现Mesh Tool菜单， 单击Mesh，出现Mesh Areas，在文本框中 输入2，单击OK。 QUST 精品课件 79 k添加标题，操作如下GUI：Utility Menu File Change Title k标题名：Element in model kGUI：Utility Menu FileSave as，输 入文件名为“exercise

39、11.db” QUST 精品课件 80 n(5)加载 n求解之前首先要选择分析类型，然后定义边界 条件及其载荷步选项，最后计算。首先选择分 析类型。 n选择Transient分析，操作如下： nGUI：Main Menu Preprocessor Loads Analysis Type New Analysis QUST 精品课件 81 n选择Transient 分析，单击OK。采用ANSYS默 认设置，在弹出的子对话框中单击OK。 n命令：ANTYPE,4 nTRNOPT,FULL nLUMPM,0 QUST 精品课件 82 k进入Time and time Step Optios 对话框，

40、操作如 下： kGUI：Main Menu Solution Load Step Opts Time/Frequenc Time Integrationamplitude Decay k在Time Integration Controls对话框中：GAMMA： 0.005 kTheta：1；OSLM：0.5；TOL：0.2 QUST 精品课件 83 k选择体 k在UI：Utility Menu Select Entities 选择： ElementBy AttributesMaterial Num，在文本框 中输入1，单击OK。 k在UI：Utility Menu Select Everyth

41、ing 选择： Nodesattached to Element，单击OK。 kGUI：Main Menu Solution Define Loads Apply Thermal Temperature On Lines QUST 精品课件 84 k定义热约束 kGUI：Main Menu Solution Define Loads Apply Thermal Temperature On Nodes， k在弹出对话框中，单击Pick all，在Apply temp on nodes对话框中，在LAB2中，选择TEMP选项，在 文本框中输入如900，单击OK。 QUST 精品课件 85 n选择

42、体 k在GUI：Utility Menu Select Everything k在GUI：Utility Menu Select Entities 选择： ElementBy AttributesMaterial Num，在文本框 中输入2，单击OK。 k在GUI：Utility Menu Select Everything 选择： Nodesattached to Element，单击OK。 QUST 精品课件 86 kGUI：Main Menu Solution Define Loads Apply Thermal Temperature On Nodes k在弹出对话框中，单击Pick

43、all，在Apply temp on nodes对话框中，在LAB2中，选择 TEMP选项，在文本框中输入如20，单击OK。 k在GUI：Utility Menu Select Everything QUST 精品课件 87 n6.求解 kGUI：Main Menu Solution Solve Current LS ， 进行求解，求解结束，点击OK，结束 kGUI：Main Menu Solution Load Step Opts Time/Frequenc Time IntegrationAmplitude decay，出现Time Integration Controls对话框，将 TI

44、MINT选项设置为“ON”，单击OK按钮 kGUI：Main Menu Solution Analysis Type Soln Controls，选择Basic选项， n如图所示 QUST 精品课件 88 QUST 精品课件 89 kGUI：Main Menu Solution Solve Current LS ，进行求解，求解结束，点击 OK，结束。 kGUI：Utility Menu FileSAVE as 名称： exercise2 QUST 精品课件 90 n7. 查看求解结果 k等值线查看 kGUI：Main Menu General PostprocRead resultlast

45、set kGUI：Utility Menu Select Entities 选择： ElementBy AttributesMaterial Num，在 文本框中输入1，单击OK。 QUST 精品课件 91 k GUI：Utility Menu Select Everything 选择：Nodesattached to Element，单击 OK。 k GUI：Main Menu General PostprocPlot ResultsContour Plot Select Everything k 在GUI：Utility Menu Select Entities 选择：ElementBy

46、AttributesMaterial Num，在文本框中输入2，单击OK。 QUST 精品课件 93 k在GUI：Utility Menu Select Everything 选择：Nodesattached to Element，单击OK。 kGUI：Main Menu General PostprocPlot ResultsContour Plot PlotCtrls Graphs Modify Axis 如图进行设置。 nGUI：Utility Menu PlotCtrlsGraphs Modify Curves命令，在/GTHK Thickness of curves 下拉列表框中选择

47、Triple选项，单击OK。 QUST 精品课件 95 QUST 精品课件 96 nGUI：Utility Menu PlotCtrlsStylesColorsGraph Colors命令， 在Curve Graph number 1 下拉列表框中选择 黄色，单击OK。 nGUI：Utility MenuselectEverything 命令 QUST 精品课件 97 kGUI：Main MenuTimehist PostproDefine Variable命令，出现defined Time-history variable， 单击Add按钮，出现Add Time-history variab

48、le对 话框，选中Nodal DOF result，单击OK。 k出现Define Nodal Data 菜单，在文本框中输入1， 单击OK，出现Define Nodal Data对话框，进行如 下图的设置 QUST 精品课件 98 QUST 精品课件 99 kGUI：Main MenuTimehist PostproGraph Variable 命令，出现Graph Time-history Variable 对话框，在NAVR1 1st variable to graph 文本框中 输入2，单击OK，显示球心温度随时间的变化关系 图 QUST 精品课件 100 n8.命令流文件 k /FI

49、LNAME, EXERCISE1 k /TITLE,TRANSIENT THERMAL ANALYSIS TO A SPHERE kKEYW, PR_THERM, !指定分析类型为热分析 QUST 精品课件 101 n/PREP7 !进入前处理器 nET,1,PLANE55!定义单元类型 nKEYOPT,1,3,1!定义单元关键项 nMP, KXX,1,70 !定义钢的导热系数 nMP, DENS,1,7800!定义钢的密度 nMP, C,1,448!定义钢的比热容 nMP, KXX, 2, 0.61!定义水的导热系数 nMP, DENS, 2,1000!定义水的密度 nMP, C, 2, 4185!定义水的比热容 QUST 精品课件 102 kRECTNG, 0, 0.3, 0, 0.3!生成矩形面 kPCIRC, 0.06,0, 0, 90!生成1/4圆面 kAOVERLAP, ALL!面叠加操作 kNUMCMP, AREA!压缩面编号 kNUMCMP, LINE!压