ANSYSWORKBENCH疲劳分析指南_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
ANSYSWORKBENCH疲劳分析指南
疲​劳​分​析​指​南
阅读已结束，如果下载本文需要使用
想免费下载本文？
你可能喜欢&& 查看话题
ansys 简单的瞬态热分析问题
本人ANSYS 新手，最近在学习瞬态热分析，一个简单的长方形一边设置50度，环境温度25度， 模拟简单的热传导过程。
但是无论把时间设置多长，长方形另一边的温度曲线总是斜线上升，而不是趋于稳态的曲线~~ 请大神指点是哪点出错了？
ET,1,PLANE13
KEYOPT,1,1,4
MP,DENS,1,7800
MP,C,1,460
MP,KXX,1,66.6
MP,ALPX,1,1.06E-5
MP,EX,1,206E9
MP,PRXY,1,0.3
RECTNG,0,,0,100/10000
ESIZE,10/10000
DL,4,,UX,0,1
DL,4,,UY,0,1
DL,4,,TEMP,50
LSEL,S,EXT
SFL,ALL,CONV,10,,25,
ALLSEL,ALL
ANTYPE, TRANS
TIME,36000
DELTIM,000
OUTRES,,ALL
PLNSOL,TEMP,,0,
呵呵，14.0热分析第一道例题，清华弄的教程 问题解决了 原因是 加载热载荷的时候没有使用KBC,1 :sweat: 将KEYOPT,1,1,4改为keyopt,1,1,2后就可以。KP2仅是temp自由度，kp4中也含有temp自由度，但是计算结果却是不一样。 : Originally posted by 不死蜗牛 at
将KEYOPT,1,1,4改为keyopt,1,1,2后就可以。KP2仅是temp自由度，kp4中也含有temp自由度，但是计算结果却是不一样。 您好， 请教您一个问题
&&在瞬态热分析中，直接法，采用 平面13单元， 单元option 设置为 UX, UY, AZ, TEMP,
& &但是 采用间接法时候， 平面问题 thermal&&solid 单元约束自由度只有 UX, UY ，& &这个会导致 间接法和 直接法 做出来的结果不同吗？
& &因为我直接法和间接 法 做出来的 结果不同，&&我不知道 原因在哪里？谢谢啊， : Originally posted by 仰望星空90 at
您好， 请教您一个问题
&&在瞬态热分析中，直接法，采用 平面13单元， 单元option 设置为 UX, UY, AZ, TEMP,
& &但是 采用间接法时候， 平面问题 thermal&&solid 单元约束自由度只有 UX, UY ，& &这个会导致 间接 ... 有差距是正常的，但不知道你的差多少？ 间接法 做的 热应力&&峰值在 0.6 GPa，&&直接法 的 热应力峰值在1.2GPa， 还有，当 网格划分比较密时，&&结果会完全不同，差距也很大，&&请问您， 就是 两种单元 约束自由度时候， 间接法是&&UX,&&UY&&直接法是 UX, UY, AZ, TEMP, 这个会导致 结果不同吗？&&网格太密，结果会完全不同？&&谢谢啊 怎么用ansys模拟工字钢在火灾条件下的升温曲线？问题有二，第一是怎么在热力学条件下建立工字钢模型？（现在我是用的关键点连线的方法）第二十怎么设置温度荷载随时间变化（按照火灾升温曲线，函数是T0+345log10（8*t+1),T0是室温，t是分钟），求大神们指教ANSYS热分析指南(第四章)
发表时间： && 作者: 安世亚太&&来源: e-works
关键字: &&&
& 第四章瞬态热分析 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 Menu&Preferences中对菜单进行过滤。然后进入前处理器（PREP7）完成以下工作： 定义单元类型 定义需要的单元实常数 定义材料属性 建立几何实体 划分网格 《ANSYS Modeling and Meshing Guide》中对本部分有详细说明。 4.5施加荷载和求解 在瞬态分析中，施加荷载的第一步是定义分析类型，然后为分析建立初始条件。 4.5.1指定分析类型 在这一步中，可以如下指定分析类型： GUI: Main Menu&Solution&Analysis Type&New Analysis&Transient 如果是一个新的分析，执行命令：ANTYPE,TRANSIENT,NEW 如果是重新启动以前的分析，比如，附加一个荷载。执行命令：ANTYPE,TRANSIENT,REST。（条件是先前分析的jobname.ESAV、jobname.DB等文件是可以利用的） 4.5.2为分析建立初始条件 瞬态热分析的初始条件来自于对应的一个稳态计算结果，或者直接为所有节点设定初始温度。 4.5.2.1设置均匀的初始温度 如果已知模型起始时的环境温度，可用下面的方法来设定所有节点的初始温度： 命令：TUNIF GUI:Main Menu& Solution&-Loads-&Settings&Uniform Temp 如果不在对话框中输入数据，则默认为参考温度，参考温度的值默认为零，可以如下设定参考温度： 命令：TREF GUI:Main Menu&Solution&-Loads-&Settings&Reference Temp 注意：设定均匀的初始温度，与下面的设定节点温度（自由度）不同。 命令：D GUI:Main Menu&Solution&-Loads-&Apply&-Thermal-&On Nodes 初始均匀温度仅对分析的第一个子步有效；而设定节点温度将使节点温度在整个瞬态分析过程等于指定值，除非通过下列方法删除此约束： 命令：DDELE GUI: Main Menu& Solution&-Loads-&Delete&-Thermal-Temperature&On Nodes 4.5.2.1设置非均匀的初始温度 在瞬态热分析（不是稳态热分析）中，可以指定一个和一组初始温度不均匀的节点，方法如下： 命令：IC GUI: Main Menu& Solution&Loads&Apply&-Initial Condit"n&Define 还可以对某些节点设定非均匀的初始温度，同时再设定其它节点的初始温度为均匀初始温度。要做到这点，只需要在为选择的节点定义不均匀温度之前，先定义均匀的温度就行了。 用以下命令可显示具有非均匀初始温度的节点： 命令：ICLIST GUI: Main Menu& Preprocessor&Loads&Apply&Initial Condit"n&List Picked 如果初始温度场是不均匀的且又是未知的，就必须首先作稳态热分析确定初始条件，步骤如下：指定相应的稳态分析荷载，如：温度约束，对流换热等。 关闭瞬态效应 命令：TIMINT,OFF,THERRM GUI： Main Menu&Preprocessor&Load&Time/Frequenc&Time-Time Integration 定义通常较小的一个时间值（如：1E-6秒） 命令：TIME GUI： Main Menu&Preprocessor&Load&Time/Frequenc&Time-Time Step 定义斜坡或阶越荷载，如果使用斜坡荷载，则就必须考虑相应的时间内产生的温度梯度效应。 命令：KBC GUI： Main Menu&Preprocessor&Load&Time/Frequenc&Time-Time Step 写荷载步文件 命令：LSWRITE GUI： Main Menu&Preprocessor&Load&Write LS File 对于第二个载荷步，要记住删除所有固定温度边界条件，除非能够判断那些节点上的温度确实在整个瞬态分析过程中都保持不变。同时，记住执行TIMINT,ON, THERM命令以打开瞬态效应。更多的细节，请见《ANSYS Commands Reference》中对D、DDELE、LSWRITE、SF、TIME和TIMINT等命令的详细描述。 分页 4.5.3设置荷载步选项 对热分析可以设置通用选项，非线性选项和输入控制。 4.5.3.1设置时间步的策略 对于瞬态热分析，既可以用多个载荷步完成（对于阶跃或渐变边界条件），也可以只用一个载荷步、采用表格边界条件（对于随时间任意变化在边界条件）并由一个数组参数定义时间点。表格边界条件方式仅适用于仅传热单元、热电单元、热表面效应单元、热流体单元以及这些类型单元的部分组合。 如果采用载荷步的方法，则按下述步骤进行： 1．设定每一载荷步结束时的时间： 命令：TIME GUI: Main　Menu& Solution&-Load Step Opts-Time/Frequenc&Time and Time Step 2. 设定载荷变化方式。如果载荷在这个载荷步是恒定的，需要设为阶越选项；如果载荷值随时间线性变化，则要设定为渐变选项： 命令：KBC GUI： Main Menu&Preprocessor&Load&Time/Frequenc&Time-Time Step 3．定义在本载荷步结束时的载荷数值（相关的命令及菜单路径参见表3-9）。 4．将载荷步信息写入载荷步文件： 命令：LSWRITE GUI： Main Menu&Preprocessor&Load&Write LS File 5．对于其它载荷步，重复步骤1～4即可，直到所有的载荷都已经写入到荷载步文件中。如要删除部分载荷（非温度约束），最好将其设置为在一个微小的时间段中值变为零，而不是直接删除。 如果采用表格参数定义载荷，按如下步骤进行： 1. 如《ANSYS Basic Analysis Procedures Guide》中的“采用表格数组参数施加载荷”所述，用TABLE类型的数组参数定义载荷特性（例如，载荷与时间的关系）。 2. 打开自动时间步长功能（AUTOTS,ON），定义时间步长（DELTIM）或子步数。 3. 定义时间步重置选项。可以选择在求解中不重置时间步,或基于一个已定义好的时间（关键时间）数组重置时间步,或基于一个新的关键时间数组重置时间步。 命令：TSRES GUI： Main Menu&Preprocessor&Load&Time/Frequenc&Time-Time Step GUI： Main Menu&Preprocessor&Load&Time/Frequenc&Time and Substeps 如果选择用新数组并交互时运行，此时程序要求填写一个n×1的关键时间数组。如果以批处理方式运行，则必须在执行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被忽略。 定义关键时间数组的方式如下： 命令：*DIM GUI： Utility Menu&Parameters&Array Parameters&Define/Edit 在关键时间数组中，时间值必须是升序排列的，并且不能超过由TIME命令定义的载荷步结束时间。在求解过程中,时间步可能会在数组定义的关键时刻点被重置.重置的大小基于命令DELTIM,DTIME或NSUBST,NSBSTP设置的初始时间步尺寸或子步数。 4．用一个与关键时间数组类似的n×1数组参数来指定将哪些时刻的计算结果写入结果文件。可以就利用关键时间数组，或用一个不同的数组。如果是交互式运行程序，可在此时创建一个数组或采用已有数组，如果是批处理方式运行程序，则必须在OUTRES命令之前定义该数组。 命令：OUTRES GUI：Main Menu&Solution&-Load Step Opts-Ouput Ctrls&DB/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命令的描述。按下述方式打开或关闭求解控制： 命令：SOLCONTROL GUI：Main Menu&Solution&Solution Ctrl 时间选项 该选项定义载荷步的结束时间，缺省情况下，第一个荷载步结束的时间是1.0，此后的荷载步对应的时间强逐次加1.0。 命令：TIME GUI:　Main Menu& Solution&-Load Step Opts-Time/Frequenc&Time and Substps 每载荷步中子步的数量或时间步大小 对于非线性分析，每一载荷步需要多个子步。缺省情况下每个荷载步有一个子步。 对于瞬态分析，在热梯度较大的区域（如淬火体的表面），热流方向的最大单元尺寸和能够得到好结果的最小时间步长有一个关系。在时间步保持不变的时候，更多的单元通常会得到更好的结果；但是，在网格尺寸不变的时候，子步越多，结果反而会变得更差。当采用自动时间步和代中间节点的二次单元时，ANSYS建议使用者根据输入的荷载来控制最大的时间步长，根据下面的关系来定义最小的时间步长：
其中为在热梯度最大处沿热流方向的单元长度，为扩散率，它等于导热系数除以密度与比热的乘积（）。当采用有中间节点的单元时，如果违反上述关系式，ANSYS的计算会出现不希望的振荡，计算出的温度会在物理上超出可能的范围。如果不采用带中间节点的单元，则一般不会计算出振荡的温度分布，那么上述建议的最小时间步长就有些保守。 注意：不要采用特别小的时间步长，特别是当建立初始条件时。在ANSYS中，很小的数可能导致计算错误，比如：当一个问题的时间量级为的时候，时间步长为1×10－10时就可能产生数值错误。 命令：NSUBST或DELTIM GUI: Main Menu&Solution&-Load Step Opts-&Time/Frequenc&Time and Substps 分页 4.5.4非线性选项 对于单场非线性热分析，ANSYS允许三种求解选项：FULL选项对应于缺省的全N-R算法；Quasi选项对应于在非线性热问题求解过程中有选择性地重构热矩阵－只有当非线性材料的性质改变量较大（用户控制）时，才重构热矩阵，该选项在时间步间不执行平衡迭代，材料性质根据载荷步开始时的温度来确定；Linear选项只在每个载荷步的第一个时间步构建一个热矩阵，它只适用于进行快速求解以得到一个近似的结果。 在ANSYS中，这些选项可通过THOPT命令来选择，Quasi和Linear选项直接组集热矩阵，只有ICCG和JCG求解器支持这种求解，可用EQSLV命令选择这些求解器。 对于Quasi求解选项，必须定义用于矩阵重构的材料参数改变容差，缺省的容差为0.05，对应于材料参数变化5％。Quasi选项设置一个单一的固定材料表以及在最高和最低温度之间等分的温度指针，用以计算随温度变化的材料性质。因此，采用该选项时，必须为固定材料表定义温度指针数（缺省为64）以及最高和最低温度（缺省为MPTEMP命令定义的最高和最低温度）。还可用THOPT命令定义其它非线性载荷选项。 命令：THOPT GUI:Main Menu&Solution&Analysis Options 只有存在非线性时，才需要定义非线性荷载步选项，包括： 平衡迭代次数 本选项设置每一子步允许的最大迭代次数，默认值为25，对大多数非线性热分析问题已经足够。如果打开求解控制（SOLCONTROL,ON），则缺省的迭代数介于15到26之间，根据具体的物理问题而变化。 命令：NEQIT GUI: Main Menu& Solution&-Load step opts&Nonlinear&Equilibrium Iter 自动时间步长 在瞬态分析中也称为的时间步优化，它使程序自动确定子步间的载荷增量。同时，它根据分析模型的响应情况，自动增、减时间步大小。在瞬态分析中，响应检测基于热特征值。对于THOPT,Quasi选项，时间步的修正也基于求解过程中的材料参数变化情况。如果特征值小，就采用大的时间步，反之亦然。在确定下一时间步长时，上一时间步中所进行的平衡迭代数量也是要考虑的依据之一，同时也要考虑非线性单元的状态变化。对于大多数问题，都应该打开自动时间步长功能并设置积分时间步长的上下限，上下限的设置可用NSUBST或DELTIM命令或下面相应的菜单路径，这种设置有助于控制时间步长的变化量。 GUI:　Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time and Substps 设置自动时间步选项： 命令：AUTOTS GUI:　Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time and Substps 调整自动时间步长中的默认参数值 命令：TINTP GUI: Main Menu&　Preprocessor& Loads &-Load Step Opts-Time/Frequenc& Time Integration 时间积分效应 该选项决定了是否包括结构惯性力，热容之类的瞬态效应。 注意：在瞬态分析时，时间积分效应缺省是打开的，如果将其设为OFF，ANSYS将进行一个稳态分析。指定时间积分效应，使用： 命令：TIMINT GUI: Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time Integration 瞬态积分参数 瞬态积分参数：此参数控制时间积分方案的性质并设定自动时间步长控制标准，详情请参考《ANSYS Theory Reference》。为尽量减少计算结果中的误差，可将此参数（THETA值）设为1。 命令：TIMINT GUI: Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time Integration 对收敛容差，求解结束，线性搜索，预测－矫正等选项的设置参见本书第三章。 4.5.5输出控制 本节内容可参见本书第三章。 4.5.6求解 本节内容可参见本书第三章。 4.6后处理 ANSYS提供两种后处理方式。 4.6.1通用后处理（POST1） 本节内容可参见本书第三章。 4.6.2时间历程后处理（POST26） 时间历程后处理器POST26对随时间变化的变量进行操作，ANSYS为每一个变量安排一个编号，第一号固定为时间。在时间历程后处理中首先要定义变量。 命令：NSOL,ESOL或RFORCE GUI: Main Menu&TimeHist Postproc&Define Variables 然后就可以绘制这些变量随时间变化的曲线： 命令：PLVAR GUI: Main Menu&TimeHist Postproc&Graph Variables 或列表输出： 命令：PRVAR GUI: Main Menu&TimeHist Postproc&List Variables 或仅列出极值： 命令：EXTREM GUI: Main Menu&TimeHist Postproc&List Extremes 通过在时间历程后处理器中观察模型的一些临界时刻点，就能够进一步在通用后处理器中这些时刻的结果进行后处理。此外，POST26还提供许多其它功能，更多的细节请参阅《ANSYS 基本过程指南》。 分页 4.7相变问题 ANSYS热分析最强大的功能之一就是可以分析相变问题，例如凝固或熔化等。含有相变问题的热分析是一个非线性的瞬态问题，典型的相变应用有： 金属浇铸：确定相变过程中不同点处的温度分布、相变发生的时间长度、浇铸的热效率等； 合金生产：由化学差异而不是物理差异导致相变； 热处理问题； 相变问题是一个非线性的瞬态热分析，线性的瞬态热分析与非线性的瞬态热分析之间唯一的不同是：需要考虑潜热，即在相变过程吸收或释放的热量。ANSYS通过定义材料的焓随温度变化来考虑潜热（如下图所示）。 图4-2材料焓随温度的变化
焓的单位是J/m3，是密度与比热的乘积对温度的积分：
求解相变问题，应当设定足够小的时间步长，并将自动时间步长设置为ON，以使程序在相变前、相变中、相变后自动调整时间步长； 选用低阶的热单元，例如PLANE55或SOLID70。如果必须选用高阶单元，请利用相应的单元选项设置“Diagonalized Specific heat matrix”（对于低阶单元，这是缺省值）； 在设定瞬态积分参数时，请将THETA值设置为1(默认为0.5)，以便瞬态时间积分采用欧拉向后差分算法； 命令：TINTP GUI:　Main Menu& Solution&-Load and Step Opts-&Time/Frequence&Time Intergation&THETA 线性搜索将有助于加速相变问题的求解。 命令：LNSRCH GUI: Main Menu&Solution&-Load and Step Opts-&Nonlinear&Line Search 4.8瞬态热分析的实例1 4.8.1问题描述 一个30公斤重、温度为70℃的铜块,以及一个20公斤重、温度为80℃的铁块，突然放入温度为20℃、盛满了300升水的、完全绝热的水箱中,如图所示。过了一个小时，求铜块与铁块的最高温度(假设忽略水的流动)。材料热物理性能如下：
4.8.2　菜单操作过程 4.8.2.1设置分析标题 1、选择“Utility Menu&File&Change Jobname”，输入文件名Transient1。 2、选择“Utility Menu&File&Change Title” 输入Thermal Transient Exercise 1。 4.8.2.2定义单元类型 1、选择“Main Menu&Preprocessor”，进入前处理。 2、选择“Main Menu&Preprocesor&Element Type&Add/Edit/Delete”。选择热平面单元plane77。 4.8.2.3定义材料属性 1、选择“Main Menu&Preprocessor&Material Props&Material Models”，在弹出的材料定义窗口中顺序双击Thermal选项。 2、点击Conductivity，Isotropic，在KXX框中输入383；点击Density，在DENS框中输入8898；点击Specific Heat，在C框中输入390。 3、在材料定义窗口中选择Material&New Model，定义第二种材料。 4、点击Conductivity，Isotropic，在KXX框中输入70；点击Density，在DENS框中输入7833；点击Specific Heat，在C框中输入448。 5、在材料定义窗口中选择Material&New Model，定义第三种材料。 6、点击Conductivity，Isotropic，在KXX框中输入.61；点击Density，在DENS框中输入996；点击Specific Heat，在C框中输入4185。 4.8.2.4创建几何模型 1、选择“Main Menu&Preprocessor&-Modeling-&Create&-Areas-&Retangle&By 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 Menu&Preprocessor&-Modeling-&Operate&Booleans&Overlap”，选择Pick All。 3、选择“Utility Menu&Plotctrls&Numbering&Areas, on”。 4、选择“Utility Menu&Plot&Areas”。 4.8.2.5划分网格 1、选择“Main Menu&Preprocessor&-Attributes-&Define-&All Areas”，选择材料1。 2、选择“Main Menu&Preprocessor&Meshing-&Size Cntrls-&-Manualsize-&-Global-&Size”，输入单元大小0.02。 3、选择“Main Menu&Preprocessor&Meshing-&Mesh-&-Areas-&Mapped&3 or 4 sided”，选择铜块。 4、选择“Main Menu&Preprocessor&-Attributes-&Define-&All Areas”，选择材料2。 5、选择“Main Menu&Preprocessor&Meshing-&Mesh-&-Areas-&Mapped&3 or 4 sided”，选择铁块。 6、选择“Main Menu&Preprocessor&-Attributes-&Define-&All Areas”，选择材料3。 7、选择“Main Menu&Preprocessor&Meshing-&Size Cntrls-&-Manualsize-&-Global-&Size”，输入单元大小0.05。 8、选择“Main Menu&Preprocessor&Meshing-&Mesh-&-Areas-&Free”，选择水箱。 9、选择“Utility Menu&Plot&Area”。 4.8.2.6进行稳态分析设置初始条件 1、选择“Main Menu&Solution&-Analysis Type-&New Analysis”，选择Transient，定义为瞬态分析。 2、选择“Main Menu&Solution&-Load Step Opts&Time/Frenquenc&Time Integration”，将TIMINT设定为 off，首先进行稳态分析。 3、选择“Main Menu&Solution&-Load Step Opts&Time/Frenquenc&Time-Time Step”，设定TIME为0.01、DELTIM也为0.01 4、选择“Utility Menu: Select&Entities”，在对话框中自上而下依次选择：Elements，By Attributes，Material num，在“Min， Max”框中输入3，选择From Full，点击APPLY；选择选择Nodes，Attached to, Element，点击OK。 5、选择“Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes”，选择Pick All, 输入20。 6、选择“Utility Menu: Select&Entities”，在对话框中自上而下依次选择：Elements，By Attributes，Material num，在“Min， Max”框中输入2，选择From Full，点击APPLY；选择选择Nodes，Attached to, Element，点击OK。 7、选择“Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes”，选择Pick All, 输入80。 8、选择“Utility Menu&Select&Entities”，在对话框中自上而下依次选择：Elements，By Attributes，Material num，在“Min， Max”框中输入1，选择From Full，点击APPLY；选择选择Nodes，Attached to, Element，点击OK。 9、选择“Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes”，选择Pick All, 输入70。 10、选择“Utility Menu&Select Everything”。 11、Main Menu&Solution&-Solve-&Current LS”。 分页 4.8.2.7进行瞬态分析 1、选择“Main Menu&Solution&-Load Step Opts&Time/Frenquenc&Time-Time Step”，设定TIME=3600,DELTIM=26, 最小、最大时间步长分别为2, 200, 将Autots设置为ON。 2、选择“Main Menu&Solution&-Load Step Opts&Time/Frenquenc&Time Integration”，将TIMINT设置为ON。 3、选择“Main Menu&Solution&-Loads-&Delete&-Thermal-&Temperature&On Nodes”，选择Pick All，删除稳态分析定义的节点温度。 4、选择“Main Menu&Solution&-Load Step Opts&Output Ctrls-&DB/Results”，选择Every Substeps。 5、选择“Main Menu&Solution&-Solve-&Current LS”。 4.8.2.8后处理 1、选择“Main Menu&TimeHist PostPro”，进入POST26。 2、选择“Main Menu&TimeHist PostPro&Define Variables”，点击Add，选择Solution summary，点击OK，在User specified label框中输入dtime，选择Solution Items和Step Time，点击OK定义子步时间为2号变量。 3、选择“Main Menu&TimeHist PostPro&Define Variables”，点击Add，选择Nodal result，点击OK，在User specified label框中输入T_Copper，在Node number框中输入node(0.5,0)，点击OK定义3号变量。同理可以定义其他节点解。 4、选择“Main Menu&TimeHist PostPro&Graph Virables”，输入变量代号，显示各变量随时间变化的曲线。 5、选择“Main Menu&General Postproc”，进入POST1。 6、选择“Main Menu&General Postproc&-Read Results-&Last set”。 7、选择“Main Menu&General Postproc&Plot result&Nodal Solution”，选择temperature。 4.8.3　等效的命令流方法 /filename,transient1 /title, Thermal Transient Exercise 1 !进入前处理 /prep7 et,1,plane77! 定义单元类型 mp,kxx,1,383! 定义材料热性能参数 mp,dens,1,8889!1～铜，2～铁，3～水 mp,c,1,390 mp,kxx,2,70 mp,dens,2,7837 mp,c,2,448 mp,kxx,3,0.61 mp,dens,3,996 mp,c,3,4185 !创建几何实体 rectnag,0,0.6,0,0.5 rectang,0.15,0.225,0.225,0.27 rectang,0.6-0.2-0.058,0.6-0.2,0.225,0.225+0.044 aovlap,all!布尔操作 /pnum,area,1 aplot !划分网格 aatt,1,1,1 eshape,2 esize,0.02 amesh,2 aatt,2,1,1 amesh,3 aatt,3,1,1 eshape,3 esize,0.05 amesh,4 /pnum,mat,1 eplot finish !加载求解 /solu antype,trans timint,off!先作稳态分析,确定初始条件 time,0.01!设定只有一个子步的时间很小的载荷步 deltim,0.01 esel,s,mat,,3 nsle,s d,all,temp,20 esel,s,mat,,2 nsle,s d,all,temp,80 esel,s,mat,,1 nsle,s d,all,temp,70 allsel solve!得到初始温度分布 !进行瞬态分析 time,3600 timint,on!打开时间积分 deltim,26,2,200!设置时间步长，最大及最小时间步长 autots,on!打开自动时间步长 ddelet,all,temp!删除稳态分析中定义的节点温度 outres,all,1!将每个子步的值写入数据库文件 solve finish save !进入POST26后处理 /post26 solu,2,dtime,,dtime!2～每一子步采用的时间步长 nsol,3,node(0.5,0),temp,,T_Copper!3～铜块的中心点 nsol,4,node(0.371,0.247,0),temp,,T_Iron!4～铁块的中心点 nsol,5,node(30,0,0),temp,,T_H2O_Bot!5～水箱的底部 nsol,6,node(30,50,0),temp,,T_H2O_Top!6～水箱的顶部 nsol,7,node(0,25,0),temp,,T_H2O_Left!7～水箱的左部 nsol,8,node(60,25,0),temp,,T_H2O_Right!8～水箱的右部 Plvar,2 plvar,3,4,5,6,7,8 finish !进入POST1后处理 /post1!设置为最后一个载荷子步 set,last esel,s,mat,,1 nsle,s plnsol,temp esel,s,mat,,2 nsle,s plnsol,temp finish 分页 4.9瞬态热分析的实例2 4.9.1问题描述 矩形的一个边温度恒定100℃。对边施加对流边界载荷，对流换热系数随固体表面温度而变化，如下表。初始温度为100℃,求解此矩形60秒温度分布的变化。
表面温度（℃）
对流换热系数
30 4.9.2　命令流方法 finish /clear *dim,cnvtab,table,3,,,TEMP ! 定义表格矩阵参数cnvtab，变量为TEMP cnvtab(1,0) = 20.0,50.0,100.0
cnvtab(1,1) = 10.0,20.0,30.0 /prep7 esize,0.5 et,1,55 rect,0,2,0,1 amesh,1 MP,KXX,,1.0 MP,DENS,,10.0 MP,C,,100.0 lsel,s,loc,x,0 dl,all,,temp,100 lsel,s,loc,x,2 sfl,all,conv,%cnvtab%,,20 alls /psf,conv,hcoef,2! 显示对流边界 /pnum,tabn,on! show table names lplot fini /solu antyp,trans kbc,1 DELTIM,.1,.05,10 time,60 tunif,100 outres,all,all solve finish /post1 set,last sflist,all!对流系数值 /pnum,tabn,off /psf,conv,hcoef,2! 显示对流边界. /pnum,sval,1! 显示表格参数的数值. eplot /pnum,sval,0 plns,temp fini /post26 nsol,2,2,temp prva,2 plva,2 finish 4.10《ANSYS Verification Manual》中与瞬态热分析相关的实例
VM28- Transient　Heat Transfer in an Infinite Slab
VM94- Heat Generating Plate
VM104 -Liquid-solid phase change
VM109 -Temperature gradient across a solid cylinder
VM110 -Transient temperature distribution in a slab
VM110 - Transient Temperature Distribution in a Slab
VM111 - Cooling of a Spherical Body
VM112 - Cooling of a Spherical Body
VM113 - Transient　Temperature Distribution in an Orthotropic Metal Bar
VM114 - Temperature　Response to a Linearly Rising Surface Temperature
VM115 - Thermal Response of a Heat-generating Slab
VM116 - Heat-conducting Plate with Sudden Cooling
VM159 - Temperature Controlled Heater
VM192 - Cooling of a Billet by Radiation
责任编辑：蒋汉桥}