1、 1 / 36第三章 稳态热分析3.1 稳态传热的定义ANSYS/Multiphysics，ANSYS/Mechanical，ANSYS/FLOTRAN 和ANSYS/Professional 这些产品支持稳态热分析。稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。也可以在所有瞬态效应消失后，将稳态热分析作为瞬态热分析的最后一步进行分析。稳态热分析可以计算确定由于不随时间变化的热载荷引起的温度、热梯度、热流率、热流密度等参数。这些热载荷包括：对流辐射热流率热流密度（单位面积热流）热生成率（单位体积热流）固定温度的边界条件稳态热分析可
2、用于材料属性固定不变的线性问题和材料性质随温度变化的非线性问题。事实上，大多数材料的热性能都随温度变化，因此在通常情况下，热分析都是非线性的。当然，如果在分析中考虑辐射，则分析也是非线性的。3.2 热分析的单元ANSYS 和 ANSYS/Professional 中大约有 40 种单元有助于进行稳态分析。有关单元的详细描述请参考ANSYS Element Reference，该手册以单元编号来讲述单元，第一个单元是 LINK1。单元名采用大写，所有的单元都可用于稳态和瞬态热分析。其中 SOLID70 单元还具有补偿在恒定速度场下由于传质导致的热流的功能。这些热分析单元如下： 表 3-1 二维实
3、体单元单元 维数 形状及特点 自由度PLANE35 二维 六节点三角形单元 温度（每个节点）PLANE55 二维 四节点四边形单元 温度（每个节点）PLANE75 二维 四节点谐单元 温度（每个节点）PLANE77 二维 八节点四边形单元 温度（每个节点）PLANE38 二维 八节点谐单元 温度（每个节点）2 / 36表 3-2 三维实体单元单元 维数 形状及特点 自由度SOLID70 三维 八节点六面体单元 温度（每个节点）SOLID87 三维 十节点四面体单元 温度（每个节点）SOLID90 三维 二十节点六单元 温度（每个节点）表 3-3 辐射连接单元单元 维数 形状及特点 自由度LIN
4、K31 二维或三维 二节点线单元 温度（每个节点）表 3-4 传导杆单元单元 维数 形状及特点 自由度LINK32 二维 二节点线单元 温度（每个节点）LINK33 三维 二节点线单元 温度（每个节点）表 3-5 对流连接单元单元 维数 形状及特点 自由度LINK34 三维 二节点线单元 温度（每个节点）表 3-6 壳单元单元 维数 形状及特点 自由度SHELL57 三维 四节点四边形单元 温度（每个节点）表 3-7 耦合场单元单元 维数 形状及特点 自由度PLANE13 二维 四节点热应力耦合单元温度、结构位移、电位、磁矢量位CONTACT48 二维 三节点热应力接触单元 温度、结构位移CO
5、NTACT49 三维 热应力接触单元 温度、结构位移FLUID116 三维 二或四节点热流单元 温度、压力SOLID5 三维八节点热应力和热电单元温度、结构位移、电位、磁标量位SOLID98 三维十节点热应力和热电单元温度、结构位移、电位、磁矢量位PLANE67 二维 四节点热电单元 温度、电位LINK68 三维 两节点热电单元 温度、电位SOLID69 三维 八节点热电单元 温度、电位3 / 36SHELL157 三维 四节点热电单元 温度、电位表 3-8 特殊单元单元 维数 形状及特点 自由度MASS71 一维到三维 一个节点的质量单元 温度COMBINE37 一维 四节点控制单元温度、结
6、构位移、转动、压力SURF151 二维 二到四节点面效应单元 温度SURF152 三维 四到九节点面效应单元 温度MATRIX50由包括在超单元中的单元类型决定没有固定形状的矩阵或辐射矩阵超单元由包括在超单元中的单元类型决定INFIN9 二维 二节点无限边界单元 温度、磁矢量位INFIN47 三维 四节点无限边界单元 温度、磁矢量位COMBINE14 一维到三维 两节点弹簧阻尼单元温度、结构位移、转动、压力COMBINE39 一维 两节点非线性弹簧单元温度、结构位移、转动、压力COMBINE40 一维 两节点组合单元温度、结构位移、转动、压力.3 热分析的基本过程ANSYS 热分析包含如下三个
7、主要步骤：前处理：建模求解：施加荷载并求解后处理：查看结果以下的内容将讲述如何执行上面的步骤。首先，对每一步的任务进行总体的介绍，然后通过一个管接处的稳态热分析的实例来引导读者如何按照 GUI 路径逐步完成一个稳态热分析。最后，本章提供了该实例等效的命令流文件。3.4 建模建立一个模型的内容包括：首先为分析指定 jobname 和 title；然后在前处理器（PREP7）中定义单元类型，单元实常数，材料属性以及建立几何实体。ANSYS Modeling and Meshing Guide中对本部分有详细说明。对于热分析有：定义单元类型4 / 36命令：ETGUI：Main MenuPrepro
8、cessorElement TypeAdd/Edit/Delete定义固定材料属性命令：MPGUI：Main MenuPreprocessorMaterial PropsMaterial ModelsThermal定义温度相关的材料属性，首先要定义温度表，然后定义对应的材料属性值。通过下面的方法定义温度表命令：MPTEMP 或 MPTEGN，然后定义对应的材料属性，使用 MPDATAGUI：Main MenuPreprocessorMaterial Props Material ModelsThermal对于温度相关的对流换热系数也是通过上述的 GUI 路径和命令来定义的。注意-如果以多项式的
9、形式定义了与温度相关的膜系数，则在定义其它具有固定属性的材料之前，必须定义一个温度表。创建几何模型及划分划分网格的过程，请参阅ANSYS Modeling and Meshing Guide3.5 施加荷载和求解在这一步骤中，必须指定所要进行的分析类型及其选项，对模型施加荷载，定义荷载选项，最后执行求解。3.5.1 指定分析类型在这一步中，可以如下指定分析类型：GUI: Main MenuSolutionNew AnalysisSteady-state(static)命令：ANTYPE,STATIC,NEW如果是重新启动以前的分析，比如，附加一个荷载。命令：ANTYPE,STATIC,rest
10、。（条件是先前分析的 jobname.ESAV、jobname.DB 等文件是可以利用的）3.5.2 施加荷载5 / 36可以直接在实体模型（点、线、面、体）或有限元模型（节点和单元）上施加载荷和边界条件，这些载荷和边界条件可以是单值的，也可以是用表格或函数的方式来定义复杂的边界条件，详见ANSYS 基本分析过程指南。可以定义以下五种热载荷：3.5.2.1 恒定的温度(TEMP)通常作为自由度约束施加于温度已知的边界上。3.5.2.2 热流率（HEAT）-可以的话，就避免了吧，它是用于提高精度的补充。热流率作为节点集中载荷，主要用于线单元（如传导杆、辐射连接单元等）模型中，而这些线单元模型通常
11、不能直接施加对流和热流密度载荷。如果输入的值为正，表示热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上，则温度约束条件优先。注意-如果在实体单元的某一节点上施加热流率，则此节点周围的单元应该密一些；特别是与该节点相连的单元的导热系数差别很大时，尤其要注意，不然可能会得到异常的温度值。因此，只要有可能，都应该使用热生成或热流密度边界条件，这些热荷载即使是在网格较为粗糙的时候都能得到较好的结果。3.5.2.3 对流（CONV）对流边界条件作为面载施加于分析模型的外表面上，用于计算与模型周围流体介质的热交换，它仅可施加于实体和壳模型上。对于线单元模型，可以通过对流杆单元 LINK34
12、 来定义对流。3.5.2.4 热流密度（HEAT）热流密度也是一种面载荷。当通过单位面积的热流率已知或通过 FLOTRAN CFD 的计算可得到时，可以在模型相应的外表面或表面效应单元上施加热流密度。如果输入的值为正，表示热流流入单元。热流密度也仅适用于实体和壳单元。单元的表面可以施加热流密度也可以施加对流，但 ANSYS 仅读取最后施加的面载进行计算。3.5.2.5 热生成率（HGEN）热生成率作为体载施加于单元上，可以模拟单元内的热生成，比如化学反应生热或电流生热。它的单位是单位体积的热流率。下表总结了在热分析中的载荷类型：表 3-9 热荷载类型6 / 36载荷类型 类别 命令 族 GUI
13、 路径温度 (TEMP) 约束 D Main MenuSolution-Loads-Apply -Thermal-Temperature热流率 (HEAT) 力 F Main MenuSolution-Loads-Apply -Thermal-Heat Flow对流 (CONV), 热流密度 (HFLUX)面载荷 SFMain MenuSolution-Loads-Apply -Thermal-ConvectionMain MenuSolution-Loads-Apply -Thermal-Heat Flux热生成率 (HGEN) 体载荷 BF Main MenuSolution-Loads-
14、Apply -Thermal-Heat Generat下表详细列出了热分析中用于施加载荷，删除载荷，对载荷进行操作、列表的所以命令：表 3-10 热荷载相关的命令载荷类型实体或有限元模型实体施加 删除列表显示 运算 设置温度 实体模型关键点DK DKDELE DKLIST DTRAN -“有限元模型节点 D DDELE DLIST DSCALE DCUMTUNIF热流率实体模型关键点FK FKDELE FKLIST FTRAN -“有限元模型节点 F FDELE FLIST FSCALE FCUM对流, 热流密度实体模型 线 SFL SFLDELESFLLIST SFTRAN SFGRAD“实
15、体模型 面 SFA SFADELESFALIST SFTRAN SFGRAD“有限元模型节点 SF SFDELE SFLIST SFSCALESFGRADSFCUM“有限元模型单元 SFE SFEDELESFELIST SFSCALESFBEAMSFCUMSFFUNSFGRAD生热 实体模 关 BFK BFKDELEBFKLIST BFTRAN -7 / 36率 型 键点“实体模型 线 BFL BFLDELEBFLLIST BFTRAN -“实体模型 面 BFA BFADELEBFALIST BFTRAN -“实体模型 体 BFV BFVDELEBFVLIST BFTRAN -“有限元模型节点
16、 BF BFDELE BFLIST BFSCALEBFCUM“ “单元 BFE BFEDELEBFELIST BFSCALEBFCUM3.5.3 采用表格和函数边界条件除了一般的使用表格来定义边界条件的方法，本节讨论热分析中特有的一些问题。关于定义表参数的详细叙述，请参考ANSYS APDL Programmers Guide。本节内容对单元类型没有特别的限制。下表列出了热分析中能够用于每一种边界条件的自变量：表 3-11 荷载边界条件及其自变量热边界条件 命令族 自变量固定温度 D TIME, X, Y, Z热流 F TIME, X, Y, Z, TEMP对流换热系数 (对流) SF TIM
17、E, X, Y, Z, TEMP, VELOCITY环境温度 (对流) SF TIME, X, Y, Z热流密度 SF TIME, X, Y, Z, TEMP热生成 BF TIME, X, Y, Z, TEMP流率 SFE TIME压力 D TIME, X, Y, Z后面有一个例题详细介绍在一个稳态热分析中如何采用表格边界条件。为了使用更加灵活的热传导系数，可以使用函数的方式来定义边界条件。有关这种用法的详细说明，可以参考ANSYS Basic Analysis Procedures Guide。除了上述自变量外，函数边界条件还可用下面的参数作为函数的自变量：表面温度（TS）（SURF151、
18、SURF152 单元的表面温度）密度（ ）（材料属性 DENS）8 / 36比热（材料属性 C）导热率（材料属性 kxx）导热率（材料属性 kyy）导热率（材料属性 kzz）粘度（材料属性 ）辐射率（材料属性 ）3.5.4 定义载荷步选项对于一个热分析，可以确定通用选项、非线性选项以及输出控制。下表列出了热分析中可能用到的载荷步选项：表 3-12 分析中的载荷步选项选项 命令 GUI 路径通用选项时间 TIMEMain MenuSolution-Load Step Opts-Time/FrequencTime-Time Step时间步数 NSUBSTMain MenuSolution -Loa
19、d Step Opts-Time/FrequencTime and Substps时间步长 DELTIMMain MenuSolution -Load Step Opts-Time /FrequencTime-Time Step阶跃或斜坡加载 KBCMain MenuSolution -Load Step Opts-Time /FrequencTime -Time Step非线性选项最大平衡迭代数 NEQITMain MenuSolution -Load Step Opts-NonlinearEquilibrium Iter自动时间步长 AUTOTSMain MenuSolution -Loa
20、d Step Opts-Time/FrequencTime-Time Step收敛容差 CNVTOLMain MenuSolution -Load Step Opts-NonlinearConvergence Crit求解中断选项 NCNV Main MenuSolution -Load Step Opts 9 / 36-NonlinearCriteria to Stop线性搜索选项 LNSRCHMain MenuSolution -Load Step Opts -NonlinearLine Search预测矫正因子 PREDMain MenuSolution -Load Step Opts
21、-NonlinearPredictor输出控制选项打印输出 OUTPR Main MenuSolution -Load Step Opts-Output CtrlsSolu Printout数据库和结果文件输出 OUTRESMain MenuSolution -Load Step Opts-OutputCtrlsDB/Results File结果外插 ERESXMain MenuSolution -Load Step Opts-OutputCtrlsIntegration Pt3.5.5 通用选项时间选项该选项定义载荷步的结束时间，虽然对于稳态热分析来说，时间选项并没有实际的物理意义，但它提供
22、了一个方便的设置载荷步和载荷子步的方法。缺省情况下，第一个荷载步结束的时间是 1.0，此后的荷载步对应的时间强逐次加 1.0。每载荷步中子步的数量或时间步大小对于非线性分析，每一载荷步需要多个子步。缺省情况下每个荷载步有一个子步。阶跃或斜坡加载如果定义阶跃载荷，则载荷值在这个载荷步内保持不变；如果为斜坡加载，则载荷值在当前载荷步的每一子步内线性变化。3.5.6 非线性选项如果存在非线性则需要定义非线性荷载步选项，包括平衡迭代次数本选项设置每一子步允许的最大迭代次数，默认值为 25，对大多数非线性热分析问题已经足够。10 / 36自动时间步长对于非线性问题，可以自动设定子步间载荷的增量，保证求解
23、的稳定性和准确性。收敛容差只要运算满足所说明的收敛判据，程序就认为它收敛，收敛判据可以基于温度、也可以是热流率，或二者都有。在实际定义时，需要说明一个典型值(CNVTOL 命令的 VALUE 域)和收敛容差(TOLER 域)，程序将 VALUE*TOLER 的值视为收敛判据。例如，如说明温度的典型值为 500，容差为 0.001，那么收敛判据则为 0.5 度。对于温度，ANSYS 将连续两次平衡迭代之间节点上温度的变化量（ ）与收敛准则进行比较来判断是否收敛。就上面的例子来说，如果在某两次平衡迭代间，每个节点的温度变化都小于 0.5 度，则认为求解收敛。对于热流率，ANSYS 比较不平衡载荷矢
24、量与收敛标准。不平衡载荷矢量表示所施加的热流与内部（计算）热流之间的差值。ANSYS 公司推荐 VALUE 值由缺省确定，TOLER 的值缺省为 1.0e-3。求解结束选项假如在规定平衡迭代数内，其解并不收敛，那么 ANSYS 程序会根据用户设置的终止选项，来决定程序停止计算或是继续进行下一个载荷步。线性搜索设置本选项可使 ANSYS 用 Newton-Raphson 方法进行线性搜索预测矫正本选项在每一子步的第一次迭代时，对自由度求解进行预测矫正。3.5.6.1 用图形跟踪收敛进行非线性热分析时，ANSYS 在每次平衡迭代完成后，都计算收敛范数，并与相应的收敛标准比较。不管是使用在批处理还是
25、交互式方式的方法，都可以在计算过程中，使用图形求解跟踪（GST）来显示计算的收敛范数和收敛标准。在交互式时，缺省为图形求解跟踪（GST）打开，批处理运行时，缺省为 GST 关闭。使用下面的方法可以，可打开或关闭 GST：命令：/GST11 / 36GUI：Main MenuSolutionLoad Step Opts-Output CtrlsGrph Solu Track下图是一个典型的 GST 图形。图 3-1 使用 GST 追踪收敛范数3.5.7 输出控制可以控制下列三种输出：控制打印输出本选项控制将何种结果数据输出到打印输出文件（jobname.out）中。控制数据库和结果文件输出该选项
26、控制将何种结果数据输出到结果文件（jobname.rth）中。外推结果该选项可将单元积分点结果拷贝到节点上，而不是按常规的方式外推到节点上（缺省采用外推方式）。3.5.8 定义分析选项12 / 36可考虑的分析选项有：Newton-Raphson 选项。该选项仅对非线性分析有用，用以定义在求解过程中切线矩阵的更新频率，有四种选择：1Program-chosen （程序选择，此为默认值，在热分析中建议采用）2Full（完全法）3Modified（修正法）4Initial Stiffness（初始刚度法）注意-对于单物理场非线性热分析，ANSYS 通常采用全 N-R 算法。要定义该选项，或打开/关
27、闭 N-R 自适应下降功能（只对全 N-R 法有效），方发如下：命令：NROPTGUI：Main MenuSolutionUnabridged MenuAnalysis Options选择求解器ANSYS 中可以选择下列的求解器：1Sparse 求解器（静态和全瞬态分析的默认求解器）2Frontal 求解器3Jacobi Conjugate Gradient(JCG) 求解器4JCG out-of-memory 求解器5Incomplete Cholesky Conjugate Gradient(ICCG) 求解器6Pre-Conditioned Conjugate Gradient (PCG
28、) 求解器7PCG out-of-memory 求解器8Algebraic Multigrid (AMG) 求解器9Distributed Domain Solver (DDS) 求解器10Iterative（程序自动选择求解器）13 / 36注意-AMG 和 DDS 求解器属并行求解器，需要单独的 ANSYS 产品支持。在ANSYS Advanced Analysis Techniques Guide中对并行求解有更详细描述。选择求解器的方法如下：命令：EQSLV+-GUI: Main MenuSolutionAnalysis Options注意：对于不含超单元（辐射分析中用 AUX12 可
29、产生超单元）的热分析模型，可选用 Iterative（快速求解）求解器，但对于含相变的传热问题，则不建议采用（可用 sparse 或 frontal 求解器）。该求解器在解算过程中不生成Jobname.EMAT 和 Jobname.EROT 文件。定义温度偏移温度偏移为当前所采用温度系统的零度与绝对零度之间的差值。温度偏移包含在相关单元（诸如有辐射效应或蠕变特性的的单元）计算中。偏移温度输入可以是摄氏度，也可以是华氏度，在进行热辐射分析时，要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度，此值应设定为 273；如果使用的是华氏度，则为 460。在后处理中，不同的温度可以用同样的方法进
30、行处理。设置温度偏移的方式如下：命令：TOFFSTGUI：Main MenuSolutionAnalysis Options3.5.9 保存模型在完成了加载和指定分析类型后，通常建议保存数据库文件，以备在求解过程中由于计算机系统故障而导致数据丢失后能够恢复数据。命令: SAVEGUI: 点击 ANSYS 工具条 SAVE_DB3.5.10 求解命令: SOLVEGUI: Main MenuSolutionCurrent LS14 / 363.5.11 后处理ANSYS 将热分析的结果写入热结果文件 jobname.rth 中，该文件包含如下数据:基本数据：节点温度导出数据：节点及单元的热流密度
31、(TFX, TFY, TFZ, TFSUM)节点及单元的热梯度(TGX, TGY, TGZ, TGSUM)单元热流率节点的反作用热流率其它可以用通用后处理器 POST1 进行后处理，下面将讲述在热分析中典型的后处理功能。关于后处理的完整描述，可参阅ANSYS Basic Analysis Procedures Guide。注意：在通用后处理器中查看结果时，数据库必须与结果有相同的模型（可以使用命令 RESUME 恢复模型）。此外，结果文件 jobname.rth 必须可用。3.5.12 读入结果进入 POST1 后，首先应读入想要看的载荷步和子步的计算结果：命令: SETGUI: Main M
32、enuGeneral Postproc-Read Results-By Load Step可通过编号选择要读入的载荷步，可以直接读入第一载荷步、或最后载荷步、或下一载荷步等。如果是使用 GUI，将会出现一个对话框提示选择要读入的荷载步。用 SET命令的 TIME 域可读入指定时刻的计算结果，如在指定时刻无计算结果，则程序根据附近时间点的值线性插值计算得到此时刻的结果。3.5.13 查看结果图 3-215 / 36结果显示云图彩色云图显示命令：PLESOL，PLETAB 或 PLNSOLGUI：Main MenuGeneral PostprocPlot ResultsElement SoluMa
33、in MenuGeneral PostprocPlot ResultsElem TableMain MenuGeneral PostprocPlot ResultsNodal Solu矢量图显示命令：PLVECTGUI：Main MenuGeneral PostprocPlot ResultsPre-defined or Userdefined图 3-2 矢量结果显示16 / 36列表显示命令：PRESOL，PRNSOL，PRRSOLGUI：Main MenuGeneral PostprocList ResultsElement SolutionMain MenuGeneral Postpro
34、cList ResultsNodal SolutionMain MenuGeneral PostprocList ResultsReaction Solu3.6 稳态热分析的实例 1带接管的圆筒罐本例讲述了如何逐步对一个带接管的圆筒罐进行稳态热分析，包括批处理的方式和 GUI 的方式。3.6.1 问题描述（恒定的温度，热流率，对流，热流密度，热生成率）本例题的主要部分为一个圆筒形罐，其上沿径向有一材料一样的接管（如图 4 所所示），罐内流动着 450F（232C）的高温流体，接管内流动着100F（38 C）的低温流体，两个流体区域由薄壁管隔离。罐的对流换热系数为 250Btu/hr-ft2-o
35、F（1420watts/m 2-K），接管的对流换热系数随管壁温度而变，它的热物理性能如表 3-13 所示。要求计算罐与接管的温度分布。注意：本例只是很多可能的热分析中的一个，并不是所有的热分析都完全按照与本例相同的步骤。材料属性及其周围的环境条件决定了一个分析应该包括哪些步骤。17 / 36表 3-13 实例的材料属性温度（TMP） 70 200 300 400 500 oF密度(DENS) 0.285 0.285 0.285 0.285 0.285 lb/in3导热系数() 8.35 8.90 9.35 9.8 10.23 Btu/hr-ft-oF比热 0.113 0.117 0.119
36、0.122 0.125 Btu/lb-oF对流换热系数(CON) 426 405 352 275 221 Btu/hr-ft2-oF图 3-3 圆柱罐与接管的相接模型（所有单位均为英制）3.6.2 分析方法取四分之一对称模型进行分析。假定罐体足够长，使其端部温度能保持常数华氏 450 度。同样的假设也用于 Y=0 的平面上。建模时，先定义两个圆柱体，再进行“overlap”布尔运算。采用映射网格划分（全六面体网格），分网时可能会出现警告信息说有扭曲单元存在，但可以不理会该警告，因为所产生的扭曲单元远离所关心的区域（两个柱体相交处）。由于材料性质与温度相关，该分析需要多个子载荷步（本例用了 50
37、 个子载荷步），同时，采用了自动时间步长功能。求解完毕后，温度云图和热流密度向量图详细显示了计算结果。3.6.3 菜单操作过程3.6.3.1 设置分析标题18 / 361、选择“Utility MenuFileChange Title”。2、输入“Steady-State analysis of pipe junction”，点击“OK”。3.6.3.2 设置单位制在命令提示行输入/UNITS,BIN（该命令无法通过菜单完成）。3.6.3.3 定义单元类型1、选择“Main MenuPreprocesorElement TypeAdd/Edit/Delete”。2、点击 Add,打开单元类型库
38、对话框。3、选择 Thermal Solid，Brick 20 node 90 号单元，点击 OK 和 Close 关闭单元选择菜单。3.6.3.4 定义材料属性1、选择“Main MenuPreprocessorMaterial PropsMaterial Models”在弹出的材料定义窗口中顺序双击 Thermal 选项。2、点击 Density，弹出一个对话框，在 DENS 框中输入 0.285，材料编号 1出现在材料定义窗口的左边。3、在材料定义窗口中顺序双击 Conductivity，Isotropic，弹出一个对话框。4、点击 Add Temperature 按钮四次，增加四列。5
39、、在 T1 到 T5 域，分别输入温度值： 70，200，300，400，500，选择温度行，用 Ctrl-C 拷贝温度值。6、在 KXX 框中，按温度的顺序，序列输入下列值（注意，各 KXX 值都要除以 12，以保证单位制一致）： 8.35/12,8.9/12,9.35/12,9.8/12,10.23/12。(为什么除以 12)导热系数7、在材料定义窗口中双击 Specific Heat，弹出一个对话框8、点击 Add Temperature 按钮四次，增加四列9、将鼠标置于 T1 域，用 Ctrl-v 粘贴 5 个温度值10、在 C 框中，按温度的顺序，序列输入下列值0.113，0.117
40、，0.119，0.122，0.125（比热）19 / 3611、在材料定义窗口中选择 MaterialNew Model，建立新材料号 212、在材料定义窗口，双击 Convection 或 Film Coef13、点击 Add Temperature 按钮四次，增加四列14、将鼠标置于 T1 域，用 Ctrl-v 粘贴 5 个温度值15、在 HF（对流换热系数）框中，按温度的顺序，序列输入下列值（注意，各 HF 值都要除以 144，以保证单位制一致）426/144，405/144，352/144，275/144， 221/144；16、点击 Graph 按钮，查看对流换热系数与温度的关系曲线
41、，然后点击 OK17、在材料定义窗口中选择 MaterialExit 退出材料定义窗口18、在工具栏点击 SAVE_DB 保存数据库。3.6.3.5 定义几何模型参数选择“Utility MenuParametersScalar Parameters”，输入ri1=1.3，ro1=1.5，z1=2，ri2=0.4，ro2=0.5，z2=2；然后点击 Close。3.6.3.6 创建几何模型1、选择“Main MenuPreprocessor-Modeling-Create-Volumes-CylinderBy Dimensions”,在弹出菜单的 Outer radius 框中输入ro1，Op
42、tional inner radium 框中输入 ri1，Z coordinates 框中输入 0 和Z1，Ending angle 框中输入 90。2、选择“Utility MenuWorkPlaneOffset WP by Increments”，在“XY,YZ,ZXAngles”框中输入 0,-90。3、选择“Main MenuPreprocessor-Modeling-Create-Volumes-CylinderBy Dimensions”，Outer radius 框中输入 ro2 ， Optional inner radium 框中输入 ri2，Z coordinates 框中输
43、入 0 和 Z2，Starting angle 框中输入-90，Ending angle 框中输入 0。4、选择“Utility MenuWorkPlaneAlign WP withGlobal Cartesian”。3.6.3.7 进行布尔运算选择“Main MenuPreprocessor-Modeling-Operate-Booleans-OverlapVolumes”，选择 Pick All。20 / 363.6.3.8 观察几何模型1、选择“Utility MenuPlotCtrlsNumbering”，将 volumes 设置为 ON。2、选择“Utility MenuPlotCt
44、rlsView Direction”，在“Coords of view point”框中输入-3,-1,1。3、在工具栏点击 SAVE_DB 保存数据库。3.6.3.9 删除多余实体选择“Main MenuPreprocessor-Modeling-DeleteVolume and Below”，拾取第三，四号体，或在命令输入行输入 3,4 回车。3.6.3.10 创建组件 AREMOTE本步将选择圆罐的 Y,Z 端面，并将它们定义为一个组件 AREMOTE。1、选择“Utility MenuSelectEntities”，打开选择实体对话框。2、在对话框中自上而下依次选择：Area，By l
45、ocation，Z Coordinates，在“Min， Max”框中输入 Z1，选择 From Full，点击 APPLY。3、接下来选择 Y Coordinates，在“Min， Max”框中输入 0，选择 Also Select，点击 OK。4、选择“Utility MenuSelectComp/AssemblyCreate Component”，在“Component name”框中输入 AREMOTE，在“Components is made of”菜单中选择 AREA3.6.3.11 将线叠加在面上显示1、选择“Utility MenuPlotCtrlsNumbering”，打开
46、Area 和 Line 的编号2、选择“Utility MenuPlotAreas”3、选择“Utility MenuPlotCtrlsErase Options”，将 Erase between Plots 按钮设置成 Off4、选择“Utility MenuPlotLines”5、选择“Utility MenuPlotCtrlsErase Options”，将 Erase between Plots 按钮设置成 On3.6.3.12 连接面及线21 / 36为划分映射式网格，连接端部的面和线。1、选择“Main MenuPreprocessor-Meshing-Concatenate-Ar
47、ea”，选择 Pick all。2、选择“Main MenuPreprocessor-Meshing-Mesh-Concatenate-Lines”，拾取 12 和 7 号线（或在命令行中输入 12,7 并回车），点击APPLY；拾取 10 和 5 号线（或在命令行中输入 10,5 并回车），点击 OK。3.6.3.13 设定网格密度1、选择“Main MenuPreprocessor-Meshing-Size CntrlsPicked Lines”，选择线 6 和 20，点击 OK，在 No. of element divisions 框中输入4，点击 OK。2、选择“Main MenuPr
48、eprocessor-Meshing-Size CntrlsPicked Lines”，选择线 40，点击 OK，在 No. of element divisions 框中输入 6，点击 OK。3、选择“Utility MenuSelectEverything”。4、选择“Main MenuPreprocessor-Meshing-Size Cntrls-Global-Size”，在 element edge length 框中输入 0.4，点击 OK。3.6.3.14 划分网格1、选择“Main MenuPreprocessor-MeshTool”，弹出“MeshTool”对话框，用选择 H
49、ex，Mapped。2、点击 Mesh，选择 Pick All。ANSYS 对实体划分网格后，将会弹出一个对话框显示对单元形状的警告，点击 Close 将其关闭。3、选择“Utility MenuPlotCtrlsNumbering”，将 Line, Area, and Volume 的按钮设置为 Off，点击 OK。4、在工具栏点击 SAVE_DB 保存数据库。3.6.3.15 定义求解类型及选项1、选择“Main MenuSolution-Analysis Type-New Analysis”，点击OK 以选择缺省设置（Steady-State）。2、选择“Main MenuSolution-Analysis Options”，点击 OK 以接受Newton- Raphson option 的缺省设置（Program-chosen）。3.6.3.16 设定均一的起始温度22 / 36选择“Main MenuSolution- Loads-Apply -Thermal- Temperature Uniform Temp”，在弹出窗口的 Uniform temperature 框中输入 450。3.6.3.17 施加对流载荷
