ANSYS热分析指南——ANSYS瞬态热分析.docVIP

ANSYS热分析指南(第四章)第四章瞬态热分析4.1瞬态传热的定义ANSYS/Multiphysics，ANSYS/Mechanical，ANSYS/FLOTRAN和ANSYS/Professional这些产品支持瞬态热分析。瞬态热分析用于计算一个系统的随时间变化的温度场及其它热参数。在工程上一般用瞬态热分析计算温度场，并将之作为热载荷进行应力分析。许多传热应用热处理问题，喷管，引擎堵塞，管路系统，压力容器等，都包含瞬态热分析。瞬态热分析的基本步骤与稳态热分析类似。主要的区别是瞬态热分析中的载荷是随时间变化的。为了表达随时间变化的载荷，可使用提供的函数工具描述载荷时间曲线并将该函数作为载荷施加（请参考ANSYS Basic Porcedures Guide中的“施加函数边界条件载荷”），或将载荷时间曲线分为载荷步。载荷时间曲线中的每一个拐点为一个载荷步，如下图所示:图4-1 用荷载步定义时变荷载对于每一个载荷步，必须定义载荷值及时间值，同时还需定义其它载荷步选项，如：载荷步为渐变或阶跃、自动时间步长等，定义完一个载荷步的所有信息后，将其写为载荷步文件，最后利用载荷步文件统一求解。本章对一个铸件的分析的实例对此有进一步说明。4.2瞬态热分析中使用的单元和命令瞬态热分析中使用的单元与稳态热分析相同，第三章对单元有简单的描述。要了解每个单元的详细说明，请参阅ANSYS Element Reference。要了解每个命令的详细功能，请参阅ANSYS Commands Reference。4.3瞬态热分析的过程瞬态热分析的过程为：建模施加荷载并求解在后处理中查看结果以下的内容将讲述瞬态分析的基本步骤，由于并不是每个瞬态分析的过程都一致，因此本书先对整个过程进行了一般的讲解，再进行实例的分析。4.4建模建立一个模型首先应为分析指定jobname和title。如果是运行的是GUI，可以在Main MenuPreferences中对菜单进行过滤。然后进入前处理器（PREP7）完成以下工作：定义单元类型定义需要的单元实常数定义材料属性建立几何实体划分网格ANSYS Modeling and Meshing Guide中对本部分有详细说明。4.5施加荷载和求解在瞬态分析中，施加荷载的第一步是定义分析类型，然后为分析建立初始条件。4.5.1指定分析类型在这一步中，可以如下指定分析类型：GUI: Main MenuSolutionAnalysis TypeNew AnalysisTransient如果是一个新的分析，执行命令：ANTYPE,TRANSIENT,NEW如果是重新启动以前的分析，比如，附加一个荷载。执行命令：ANTYPE,TRANSIENT,REST。（条件是先前分析的jobname.ESAV、jobname.DB等文件是可以利用的）4.5.2为分析建立初始条件瞬态热分析的初始条件来自于对应的一个稳态计算结果，或者直接为所有节点设定初始温度。4.5.2.1设置均匀的初始温度如果已知模型起始时的环境温度，可用下面的方法来设定所有节点的初始温度：命令：TUNIFGUI:Main Menu Solution-Loads-SettingsUniform Temp如果不在对话框中输入数据，则默认为参考温度，参考温度的值默认为零，可以如下设定参考温度：命令：TREFGUI:Main MenuSolution-Loads-SettingsReference Temp注意：设定均匀的初始温度，与下面的设定节点温度（自由度）不同。命令：DGUI:Main MenuSolution-Loads-Apply-Thermal-On Nodes初始均匀温度仅对分析的第一个子步有效；而设定节点温度将使节点温度在整个瞬态分析过程等于指定值，除非通过下列方法删除此约束：命令：DDELEGUI: Main Menu Solution-Loads-Delete-Thermal-TemperatureOn Nodes4.5.2.1设置非均匀的初始温度在瞬态热分析（不是稳态热分析）中，可以指定一个和一组初始温度不均匀的节点，方法如下：命令：ICGUI: Main Menu SolutionLoadsApply-Initial ConditnDefine还可以对某些节点设定非均匀的初始温度，同时再设定其它节点的初始温度为均匀初始温度。要做到这点，只需要在为选择的节点定义不均匀温度之前，先定义均匀的温度就行了。用以下命令可显示具有非均匀初始温度的节点：命令：ICLISTGUI: Main Menu PreprocessorLoadsApplyInitial ConditnList Picked如果初始温度场是不均匀的且又是未知的，就必须首先作稳态热分析确定初始条件，步骤如下：指定相应的稳态分析荷载，如：温度约束，对流换热等。关闭瞬态效应命令：TIMINT,OFF,THERRMGUI： Main MenuPreprocessorLoadTime/FrequencTime-Time Integration定义通常较小的一个时间值（如：1E-6秒）命令：TIMEGUI： Main MenuPreprocessorLoadTime/FrequencTime-Time Step定义斜坡或阶越荷载，如果使用斜坡荷载，则就必须考虑相应的时间内产生的温度梯度效应。命令：KBCGUI： Main MenuPreprocessorLoadTime/FrequencTime-Time Step写荷载步文件命令：LSWRITEGUI： Main MenuPreprocessorLoadWrite LS File对于第二个载荷步，要记住删除所有固定温度边界条件，除非能够判断那些节点上的温度确实在整个瞬态分析过程中都保持不变。同时，记住执行TIMINT,ON, THERM命令以打开瞬态效应。更多的细节，请见ANSYS Commands Reference中对D、DDELE、LSWRITE、SF、TIME和TIMINT等命令的详细描述。4.5.3设置荷载步选项对热分析可以设置通用选项，非线性选项和输入控制。4.5.3.1设置时间步的策略对于瞬态热分析，既可以用多个载荷步完成（对于阶跃或渐变边界条件），也可以只用一个载荷步、采用表格边界条件（对于随时间任意变化在边界条件）并由一个数组参数定义时间点。表格边界条件方式仅适用于仅传热单元、热电单元、热表面效应单元、热流体单元以及这些类型单元的部分组合。如果采用载荷步的方法，则按下述步骤进行：1设定每一载荷步结束时的时间：命令：TIMEGUI: MainMenu Solution-Load Step Opts-Time/FrequencTime and Time Step2. 设定载荷变化方式。如果载荷在这个载荷步是恒定的，需要设为阶越选项；如果载荷值随时间线性变化，则要设定为渐变选项：命令：KBCGUI： Main MenuPreprocessorLoadTime/FrequencTime-Time Step3定义在本载荷步结束时的载荷数值（相关的命令及菜单路径参见表3-9）。4将载荷步信息写入载荷步文件：命令：LSWRITEGUI： Main MenuPreprocessorLoadWrite LS File5对于其它载荷步，重复步骤14即可，直到所有的载荷都已经写入到荷载步文件中。如要删除部分载荷（非温度约束），最好将其设置为在一个微小的时间段中值变为零，而不是直接删除。如果采用表格参数定义载荷，按如下步骤进行：1. 如ANSYS Basic Analysis Procedures Guide中的“采用表格数组参数施加载荷”所述，用TABLE类型的数组参数定义载荷特性（例如，载荷与时间的关系）。2. 打开自动时间步长功能（AUTOTS,ON），定义时间步长（DELTIM）或子步数。3. 定义时间步重置选项。可以选择在求解中不重置时间步,或基于一个已定义好的时间（关键时间）数组重置时间步,或基于一个新的关键时间数组重置时间步。命令：TSRESGUI： Main MenuPreprocessorLoadTime/FrequencTime-Time StepGUI： Main MenuPreprocessorLoadTime/FrequencTime and Substeps如果选择用新数组并交互时运行，此时程序要求填写一个n1的关键时间数组。如果以批处理方式运行，则必须在执行TSRES命令之前定义一个数组，其将时间步重置为由DELTIM或NSUBST命令定义的初始值。如果在应用时间步重置数组（TSRES命令）的同时又采用了另外的时间值数组（OUTRES命令的FREQ = %array%），则需确认：如果FREQ数组的时间值比在TSRES数组中所对应的最接近的时间值大，则所大的数值至少应为由DELTIM或NSUBST命令定义的初始时间步增量。例如，如果FREQ数组的时间值为1.5, 2, 10, 14.1,和15，TSRES数组的时间值为1, 2, 10, 14,和16（在这些值处时间步将重新开始），初始时间步增量DTIME =0.2，则程序将停止运算，因为在这种设置下，在时间为14时时间步将重置，那么下一个时间步至少为14.2，所要求的FREQ数组时间值14.1并不存在。注意：TSRES命令只有在设置了AUTOTS,ON的情况下才有效，如果采用固定时间步长（AUTOTS,OFF）,则TRES被忽略。定义关键时间数组的方式如下：命令：*DIMGUI： Utility MenuParametersArray ParametersDefine/Edit在关键时间数组中，时间值必须是升序排列的，并且不能超过由TIME命令定义的载荷步结束时间。在求解过程中,时间步可能会在数组定义的关键时刻点被重置.重置的大小基于命令DELTIM,DTIME或NSUBST,NSBSTP设置的初始时间步尺寸或子步数。4用一个与关键时间数组类似的n1数组参数来指定将哪些时刻的计算结果写入结果文件。可以就利用关键时间数组，或用一个不同的数组。如果是交互式运行程序，可在此时创建一个数组或采用已有数组，如果是批处理方式运行程序，则必须在OUTRES命令之前定义该数组。命令：OUTRESGUI：Main MenuSolution-Load Step Opts-Ouput CtrlsDB/Results File注意：只有在采用下列仅传热单元、热电单元、热表面效应单元、流体单元（FLUID166）或这些单元的组合的情况下，才能使用TSRES命令和相应的时间步策略：LINK31、LINK32、LINK33、PLANE35、MATEIX50、PLANE55、SHELL57、PLANE67（只有热自由度）、LINK68（只有热自由度）、SOLID69（只有热自由度）、SOLID70、MASS71、PLANE75、PLANE77、SOLID87、SOLID90、FLUID116、SURF151、SURF152、SHELL157（只有热自由度）、TARGE169、TARGE170、CONTA171、CONTA172、CONTA173、CONTA174。4.5.3.2通用选项求解控制选项该选项打开或关闭ANSYS内部的求解控制功能，如果打开，则用户通常只需定义子步数（NSUBST）或时间步长（DELTIM），以及载荷步结束时间（TIME），其它的求解控制命令将由程序自动设置为其最佳值，详见ANSYS命令手册中对SOLCONTROL命令的描述。按下述方式打开或关闭求解控制：命令：SOLCONTROLGUI：Main MenuSolutionSolution Ctrl时间选项该选项定义载荷步的结束时间，缺省情况下，第一个荷载步结束的时间是1.0，此后的荷载步对应的时间强逐次加1.0。命令：TIMEGUI:Main Menu Solution-Load Step Opts-Time/FrequencTime and Substps每载荷步中子步的数量或时间步大小对于非线性分析，每一载荷步需要多个子步。缺省情况下每个荷载步有一个子步。对于瞬态分析，在热梯度较大的区域（如淬火体的表面），热流方向的最大单元尺寸和能够得到好结果的最小时间步长有一个关系。在时间步保持不变的时候，更多的单元通常会得到更好的结果；但是，在网格尺寸不变的时候，子步越多，结果反而会变得更差。当采用自动时间步和代中间节点的二次单元时，ANSYS建议使用者根据输入的荷载来控制最大的时间步长，根据下面的关系来定义最小的时间步长：其中为在热梯度最大处沿热流方向的单元长度，为扩散率，它等于导热系数除以密度与比热的乘积（）。当采用有中间节点的单元时，如果违反上述关系式，ANSYS的计算会出现不希望的振荡，计算出的温度会在物理上超出可能的范围。如果不采用带中间节点的单元，则一般不会计算出振荡的温度分布，那么上述建议的最小时间步长就有些保守。注意：不要采用特别小的时间步长，特别是当建立初始条件时。在ANSYS中，很小的数可能导致计算错误，比如：当一个问题的时间量级为的时候，时间步长为11010时就可能产生数值错误。命令：NSUBST或DELTIMGUI: Main MenuSolution-Load Step Opts-Time/FrequencTime and Substps4.5.4非线性选项对于单场非线性热分析，ANSYS允许三种求解选项：FULL选项对应于缺省的全N-R算法；Quasi选项对应于在非线性热问题求解过程中有选择性地重构热矩阵只有当非线性材料的性质改变量较大（用户控制）时，才重构热矩阵，该选项在时间步间不执行平衡迭代，材料性质根据载荷步开始时的温度来确定；Linear选项只在每个载荷步的第一个时间步构建一个热矩阵，它只适用于进行快速求解以得到一个近似的结果。在ANSYS中，这些选项可通过THOPT命令来选择，Quasi和Linear选项直接组集热矩阵，只有ICCG和JCG求解器支持这种求解，可用EQSLV命令选择这些求解器。对于Quasi求解选项，必须定义用于矩阵重构的材料参数改变容差，缺省的容差为0.05，对应于材料参数变化5。Quasi选项设置一个单一的固定材料表以及在最高和最低温度之间等分的温度指针，用以计算随温度变化的材料性质。因此，采用该选项时，必须为固定材料表定义温度指针数（缺省为64）以及最高和最低温度（缺省为MPTEMP命令定义的最高和最低温度）。还可用THOPT命令定义其它非线性载荷选项。命令：THOPTGUI:Main MenuSolutionAnalysis Options只有存在非线性时，才需要定义非线性荷载步选项，包括：平衡迭代次数本选项设置每一子步允许的最大迭代次数，默认值为25，对大多数非线性热分析问题已经足够。如果打开求解控制（SOLCONTROL,ON），则缺省的迭代数介于15到26之间，根据具体的物理问题而变化。命令：NEQITGUI: Main Menu Solution-Load step optsNonlinearEquilibrium Iter自动时间步长在瞬态分析中也称为的时间步优化，它使程序自动确定子步间的载荷增量。同时，它根据分析模型的响应情况，自动增、减时间步大小。在瞬态分析中，响应检测基于热特征值。对于THOPT,Quasi选项，时间步的修正也基于求解过程中的材料参数变化情况。如果特征值小，就采用大的时间步，反之亦然。在确定下一时间步长时，上一时间步中所进行的平衡迭代数量也是要考虑的依据之一，同时也要考虑非线性单元的状态变化。对于大多数问题，都应该打开自动时间步长功能并设置积分时间步长的上下限，上下限的设置可用NSUBST或DELTIM命令或下面相应的菜单路径，这种设置有助于控制时间步长的变化量。GUI:Main Menu Solution-Load Step Opts-Time/FrequencTime and Substps设置自动时间步选项：命令：AUTOTSGUI:Main Menu Solution-Load Step Opts-Time/FrequencTime and Substps调整自动时间步长中的默认参数值命令：TINTPGUI: Main MenuPreprocessor Loads -Load Step Opts-Time/Frequenc Time Integration时间积分效应该选项决定了是否包括结构惯性力，热容之类的瞬态效应。注意：在瞬态分析时，时间积分效应缺省是打开的，如果将其设为OFF，ANSYS将进行一个稳态分析。指定时间积分效应，使用：命令：TIMINTGUI: Main Menu Solution-Load Step Opts-Time/FrequencTime Integration瞬态积分参数瞬态积分参数：此参数控制时间积分方案的性质并设定自动时间步长控制标准，详情请参考ANSYS Theory Reference。为尽量减少计算结果中的误差，可将此参数（THETA值）设为1。命令：TIMINTGUI: Main Menu Solution-Load Step Opts-Time/FrequencTime Integration对收敛容差，求解结束，线性搜索，预测矫正等选项的设置参见本书第三章。4.5.5输出控制本节内容可参见本书第三章。4.5.6求解本节内容可参见本书第三章。4.6后处理ANSYS提供两种后处理方式。4.6.1通用后处理（POST1）本节内容可参见本书第三章。4.6.2时间历程后处理（POST26）时间历程后处理器POST26对随时间变化的变量进行操作，ANSYS为每一个变量安排一个编号，第一号固定为时间。在时间历程后处理中首先要定义变量。命令：NSOL,ESOL或RFORCEGUI: Main MenuTimeHist PostprocDefine Variables然后就可以绘制这些变量随时间变化的曲线：命令：PLVARGUI: Main MenuTimeHist PostprocGraph Variables或列表输出：命令：PRVARGUI: Main MenuTimeHist PostprocList Variables或仅列出极值：命令：EXTREMGUI: Main MenuTimeHist PostprocList Extremes通过在时间历程后处理器中观察模型的一些临界时刻点，就能够进一步在通用后处理器中这些时刻的结果进行后处理。此外，POST26还提供许多其它功能，更多的细节请参阅ANSYS 基本过程指南。4.5.4非线性选项对于单场非线性热分析，ANSYS允许三种求解选项：FULL选项对应于缺省的全N-R算法；Quasi选项对应于在非线性热问题求解过程中有选择性地重构热矩阵只有当非线性材料的性质改变量较大（用户控制）时，才重构热矩阵，该选项在时间步间不执行平衡迭代，材料性质根据载荷步开始时的温度来确定；Linear选项只在每个载荷步的第一个时间步构建一个热矩阵，它只适用于进行快速求解以得到一个近似的结果。在ANSYS中，这些选项可通过THOPT命令来选择，Quasi和Linear选项直接组集热矩阵，只有ICCG和JCG求解器支持这种求解，可用EQSLV命令选择这些求解器。对于Quasi求解选项，必须定义用于矩阵重构的材料参数改变容差，缺省的容差为0.05，对应于材料参数变化5。Quasi选项设置一个单一的固定材料表以及在最高和最低温度之间等分的温度指针，用以计算随温度变化的材料性质。因此，采用该选项时，必须为固定材料表定义温度指针数（缺省为64）以及最高和最低温度（缺省为MPTEMP命令定义的最高和最低温度）。还可用THOPT命令定义其它非线性载荷选项。命令：THOPTGUI:Main MenuSolutionAnalysis Options只有存在非线性时，才需要定义非线性荷载步选项，包括：平衡迭代次数本选项设置每一子步允许的最大迭代次数，默认值为25，对大多数非线性热分析问题已经足够。如果打开求解控制（SOLCONTROL,ON），则缺省的迭代数介于15到26之间，根据具体的物理问题而变化。命令：NEQITGUI: Main Menu Solution-Load step optsNonlinearEquilibrium Iter自动时间步长在瞬态分析中也称为的时间步优化，它使程序自动确定子步间的载荷增量。同时，它根据分析模型的响应情况，自动增、减时间步大小。在瞬态分析中，响应检测基于热特征值。对于THOPT,Quasi选项，时间步的修正也基于求解过程中的材料参数变化情况。如果特征值小，就采用大的时间步，反之亦然。在确定下一时间步长时，上一时间步中所进行的平衡迭代数量也是要考虑的依据之一，同时也要考虑非线性单元的状态变化。对于大多数问题，都应该打开自动时间步长功能并设置积分时间步长的上下限，上下限的设置可用NSUBST或DELTIM命令或下面相应的菜单路径，这种设置有助于控制时间步长的变化量。GUI:Main Menu Solution-Load Step Opts-Time/FrequencTime and Substps设置自动时间步选项：命令：AUTOTSGUI:Main Menu Solution-Load Step Opts-Time/FrequencTime and Substps调整自动时间步长中的默认参数值命令：TINTPGUI: Main MenuPreprocessor Loads -Load Step Opts-Time/Frequenc Time Integration时间积分效应该选项决定了是否包括结构惯性力，热容之类的瞬态效应。注意：在瞬态分析时，时间积分效应缺省是打开的，如果将其设为OFF，ANSYS将进行一个稳态分析。指定时间积分效应，使用：命令：TIMINTGUI: Main Menu Solution-Load Step Opts-Time/FrequencTime Integration瞬态积分参数瞬态积分参数：此参数控制时间积分方案的性质并设定自动时间步长控制标准，详情请参考ANSYS Theory Reference。为尽量减少计算结果中的误差，可将此参数（THETA值）设为1。命令：TIMINTGUI: Main Menu Solution-Load Step Opts-Time/FrequencTime Integration对收敛容差，求解结束，线性搜索，预测矫正等选项的设置参见本书第三章。4.5.5输出控制本节内容可参见本书第三章。4.5.6求解本节内容可参见本书第三章。4.6后处理ANSYS提供两种后处理方式。4.6.1通用后处理（POST1）本节内容可参见本书第三章。4.6.2时间历程后处理（POST26）时间历程后处理器POST26对随时间变化的变量进行操作，ANSYS为每一个变量安排一个编号，第一号固定为时间。在时间历程后处理中首先要定义变量。命令：NSOL,ESOL或RFORCEGUI: Main MenuTimeHist PostprocDefine Variables然后就可以绘制这些变量随时间变化的曲线：命令：PLVARGUI: Main MenuTimeHist PostprocGraph Variables或列表输出：命令：PRVARGUI: Main MenuTimeHist PostprocList Variables或仅列出极值：命令：EXTREMGUI: Main MenuTimeHist PostprocList Extremes通过在时间历程后处理器中观察模型的一些临界时刻点，就能够进一步在通用后处理器中这些时刻的结果进行后处理。此外，POST26还提供许多其它功能，更多的细节请参阅ANSYS 基本过程指南。4.7相变问题ANSYS热分析最强大的功能之一就是可以分析相变问题，例如凝固或熔化等。含有相变问题的热分析是一个非线性的瞬态问题，典型的相变应用有：金属浇铸：确定相变过程中不同点处的温度分布、相变发生的时间长度、浇铸的热效率等；合金生产：由化学差异而不是物理差异导致相变；热处理问题；相变问题是一个非线性的瞬态热分析，线性的瞬态热分析与非线性的瞬态热分析之间唯一的不同是：需要考虑潜热，即在相变过程吸收或释放的热量。ANSYS通过定义材料的焓随温度变化来考虑潜热（如下图所示）。图4-2材料焓随温度的变化焓的单位是J/m3，是密度与比热的乘积对温度的积分：求解相变问题，应当设定足够小的时间步长，并将自动时间步长设置为ON，以使程序在相变前、相变中、相变后自动调整时间步长；选用低阶的热单元，例如PLANE55或SOLID70。如果必须选用高阶单元，请利用相应的单元选项设置“Diagonalized Specific heat matrix”（对于低阶单元，这是缺省值）；在设定瞬态积分参数时，请将THETA值设置为1(默认为0.5)，以便瞬态时间积分采用欧拉向后差分算法；命令：TINTPGUI:Main Menu Solution-Load and Step Opts-Time/FrequenceTime IntergationTHETA线性搜索将有助于加速相变问题的求解。命令：LNSRCHGUI: Main MenuSolution-Load and Step Opts-NonlinearLine Search4.8瞬态热分析的实例14.8.1问题描述一个30公斤重、温度为70的铜块,以及一个20公斤重、温度为80的铁块，突然放入温度为20、盛满了300升水的、完全绝热的水箱中,如图所示。过了一个小时，求铜块与铁块的最高温度(假设忽略水的流动)。材料热物理性能如下：热性能单位制铜铁水导热系数W/m38370.61密度Kg/m388897833996比热J/kg39044841854.8.2菜单操作过程4.8.2.1设置分析标题1、选择“Utility MenuFileChange Jobname”，输入文件名Transient1。2、选择“Utility MenuFileChange Title” 输入Thermal Transient Exercise 1。4.8.2.2定义单元类型1、选择“Main MenuPreprocessor”，进入前处理。2、选择“Main MenuPreprocesorElement TypeAdd/Edit/Delete”。选择热平面单元plane77。4.8.2.3定义材料属性1、选择“Main MenuPreprocessorMaterial PropsMaterial Models”，在弹出的材料定义窗口中顺序双击Thermal选项。2、点击Conductivity，Isotropic，在KXX框中输入383；点击Density，在DENS框中输入8898；点击Specific Heat，在C框中输入390。3、在材料定义窗口中选择MaterialNew Model，定义第二种材料。4、点击Conductivity，Isotropic，在KXX框中输入70；点击Density，在DENS框中输入7833；点击Specific Heat，在C框中输入448。5、在材料定义窗口中选择MaterialNew Model，定义第三种材料。6、点击Conductivity，Isotropic，在KXX框中输入.61；点击Density，在DENS框中输入996；点击Specific Heat，在C框中输入4185。4.8.2.4创建几何模型1、选择“Main MenuPreprocessor-Modeling-Create-Areas-RetangleBy Dimensions”，输入X1=0, Y1=0, X2=0.6, Y2=0.5, 点击Apply；输入X1=0.15, Y1=0.225, X2= 0.225,Y2=0.27, 点击Apply；输入X1=0.6-0.2-0.058, Y1=0.225, X2=0.6-0.2, Y2=0.225+0.044, 选择OK。2、选择“Main MenuPreprocessor-Modeling-OperateBooleansOverlap”，选择Pick All。3、选择“Utility MenuPlotctrlsNumberingAreas, on”。4、选择“Utility MenuPlotAreas”。4.8.2.5划分网格1、选择“Main MenuPreprocessor-Attributes-Define-All Areas”，选择材料1。2、选择“Main MenuPreprocessorMeshing-Size Cntrls-Manualsize-Global-Size”，输入单元大小0.02。3、选择“Main MenuPreprocessorMeshing-Mesh-Areas-Mapped3 or 4 sided”，选择铜块。4、选择“Main MenuPreprocessor-Attributes-Define-All Areas”，选择材料2。5、选择“Main MenuPreprocessorMeshing-Mesh-Areas-Mapped3 or 4 sided”，选择铁块。6、选择“Main MenuPreprocessor-Attributes-Define-All Areas”，选择材料3。7、选择“Main MenuPreprocessorMeshing-Size Cntrls-Manualsize-Global-Size”，输入单元大小0.05。8、选择“Main MenuPreprocessorMeshing-Mesh-Areas-Free”，选择水箱。9、选择“Utility MenuPlotArea”。4.8.2.6进行稳态分析设置初始条件1、选择“Main MenuSolution-Analysis Type-New Analysis”，选择Transient，定义为瞬态分析。2、选择“Main MenuSolution-Load Step OptsTime/FrenquencTime Integration”，将TIMINT设定为 off，首先进行稳态分析。3、选择“Main MenuSolution-Load Step OptsTime/FrenquencTime-Time Step”，设定TIME为0.01、DELTIM也为0.014、选择“Utility Menu: SelectEntities”，在对话框中自上而下依次选择：Elements，By Attributes，Material num，在“Min， Max”框中输入3，选择From Full，点击APPLY；选择选择Nodes，Attached to, Element，点击OK。5、选择“Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes”，选择Pick All, 输入20。6、选择“Utility Menu: SelectEntities”，在对话框中自上而下依次选择：Elements，By Attributes，Material num，在“Min， Max”框中输入2，选择From Full，点击APPLY；选择选择Nodes，Attached to, Element，点击OK。7、选择“Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes”，选择Pick All, 输入80。8、选择“Utility MenuSelectEntities”，在对话框中自上而下依次选择：Elements，By Attributes，Material num，在“Min， Max”框中输入1，选择From Full，点击APPLY；选择选择Nodes，Attached to, Element，点击OK。9、选择“Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes”，选择Pick All, 输入70。10、选择“Utility MenuSelect Everything”。11、Main MenuSolution-Solve-Current LS”。4.8.2.7进行瞬态分析1、选择“Main MenuSolution-Load Step OptsTime/FrenquencTime-Time Step”，设定TIME=3600,DELTIM=26, 最小、最大时间步长分别为2, 200, 将Autots设置为ON。2、选择“Main MenuSolution-Load Step OptsTime/FrenquencTime Integration”，将TIMINT设置为ON。3、选择“Main MenuSolution-Loads-Delete-Thermal-TemperatureOn Nodes”，选择Pick All，删除稳态分析定义的节点温度。4、选择“Main MenuSolution-Load Step OptsOutput Ctrls-DB/Results”，选择Every Substeps。5、选择“Main MenuSolution-Solve-Current LS”。4.8.2.8后处理1、选择“Main MenuTimeHist PostPro”，进入POST26。2、选择“Main MenuTimeHist PostProDefine Variables”，点击Add，选择Solution summary，点击OK，在User specified label框中输入dtime，选择Solution Items和Step Time，点击OK定义子步时间为2号变量。3、选择“Main MenuTimeHist PostProDefine Variables”，点击Add，选择Nodal result，点击OK，在User specified label框中输入T_Copper，在Node number框中输入node(0.1875,0.2475,0)，点击OK定义3号变量。同理可以定义其他节点解。4、选择“Main MenuTimeHist PostProGraph Virables”，输入变量代号，显示各变量随时间变化的曲线。5、选择“Main MenuGeneral Postproc”，进入POST1。6、选择“Main MenuGeneral Postproc-Read Results-Last set”。7、选择“Main MenuGeneral PostprocPlot resultNodal Solution”，选择temperature。4.8.3等效的命令流方法/filename,transient1/title, Thermal Transient Exercise 1!进入前处理/prep7et,1,plane77! 定义单元类型mp,kxx,1,383! 定义材料热性能参数mp,dens,1,8889!1铜，2铁，3水mp,c,1,390mp,kxx,2,70mp,dens,2,7837mp,c,2,448mp,kxx,3,0.61mp,dens,3,996mp,c,3,4185!创建几何实体rectnag,0,0.6,0,0.5rectang,0.15,0.225,0.225,0.27rectang,0.6-0.2-0.058,0.6-0.2,0.225,0.225+0.044aovlap,all!布尔操作/pnum,area,1aplot!划分网格aatt,1,1,1eshape,2esize,0.02amesh,2aatt,2,1,1amesh,3aatt,3,1,1eshape,3esize,0.05amesh,4/pnum,mat,1eplotfinish!加载求解/soluantype,transtimint,off!先作稳态分析,确定初始条件time,0.01!设定只有一个子步的时间很小的载荷步deltim,0.01esel,s,mat,3nsle,sd,all,temp,20esel,s,mat,2nsle,sd,all,temp,80esel,s,mat,1nsle,sd,all,temp,70allselsolve!得到初始温度分布!进行瞬态分析time,3600timint,on!打开时间积分deltim,26,2,200!设置时间步长，最大及最小时间步长autots,on!打开自动时间步长ddelet,all,temp!删除稳态分析中定义的节点温度outres,all,1!将每个子步的值写入数据库文件solvefinishsave!进入POST26后处理/post26solu,2,dtime,dtime!2每一子步采用的时间步长nsol,3,node(0.1875,0.2475,0),temp,T_Copper!3铜块的中心点nsol,4