Ansys 例瞬态热分析实例一水箱
ansys稳态及瞬态热分析.ppt
数； {T}为节点温度向量； {Q}为节点热流率向量，包含热生成； ANSYS利用模型几何参数、材料热性能参数以及所施加的边界 条件，生成[K] 、 {T}以及{Q} 。
2001年10月1日 2023/11/13
*ANSYS培训教程 – 版本 5.5 – XJTU MSSV By: Haich Gao (011001)
Guidelines Them-16
第五讲、瞬态传热
瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统 的温度、热流率、热边界条件以及系统内能随时间都有明显变化。 根据能量守恒原理，瞬态热平衡可以表达为（以矩阵形式表示）：
[C]{T}+[K]{T}={Q}
式中: [K]为传导矩阵，包含导热系数、对流系数及辐射率和形状 系数； [C]为比热矩阵,考虑系统内能的增加; {T}为节点温度向量；
2001年10月1日 2023/11/13
*ANSYS培训教程 – 版本 5.5 – XJTU MSSV By: Haich Gao (011001)
Them-15
第四讲、稳态传热
如果系统的净热流率为０，即流入系统的热量加上系统自身产生的 热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳 态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析 的能量平衡方程为（以矩阵形式表示）
2001年10月1日 2023/11/13
*ANSYS培训教程 – 版本 5.5 – XJTU MSSV By: Haich Gao (011001)
Guidelines Them-19
第八讲、热分析误差估计
瞬态传热分析实例
设定一个只有一个子步的，时间很小的载荷步（例如0.001） 命令：TIME GUI：Main Menu > Solution > Analysis Type > Sol'n Controls > basic 4）写入载荷步文件 命令：LSWRITE GUI：Main Menu > Solution > Load Step Opts > Write LS File 或先求解： GUI：Main Menu > Solution > Solve > Current LS
POST1 后处理
读入结果数据
命令：SET
GUI：Main Menu > General Postproc > Read Results > By Time/Freq
进入POST1 后，可以读出某一时间点的结果。如果设定的时间点不在任何一个子步的时间点上，ANSYS 会进行线性插值。
此外还可以读出某一载荷步的结果：
3
在弹出对话框中，单击OK得到圆面。
4
命令：CYL4, 0.075, 0.025, 0.01
*
相减
根据ANSYS建模过程中面序号赋值原理，直接可以肯定圆面序号为2，矩形序号为1，因此采用直接键入命令建实体模型：
命令：asba,1,2
*
下面介绍网格尺寸的设定（SmartSize 方式）：
设定网格尺寸并划分网格
建立实体模型
根据本例所用模型，首先需要创建矩形，然后是圆，最后在矩形板中央减去（Substract）圆。下面介绍建立实体模型的操作：
创建矩形
命令：RECTNG,0,0.15,0,0.05,
*
创建圆面
ansys瞬态热分析教程及实例
大家好
QUST
18
自动时间步长选项 命令：AUTOTS GUI：Main Menu > Solution > Analysis Type
> Sol'n Controls 打开后求解过程中将自动调整时间步长。
大家好
QUST
19
时间积分选项 命令：TIMINT GUI：Main Menu > Solution > Load Step
QUST
6
大家好
5. 瞬态传热分析
均匀初始温度：如果整个模
型的初始温度为均匀且非0
2
，使用下列菜单指定:
3 4 1
QUST
7
大家好
1
2 3
5. 瞬态传热分析
非均匀的初始温度
如果模型的初始温度分布
已知但不均匀，使用这些
菜单将初始条件施加在特
4
定节点上
5
QUST
8
大家好
5. 瞬态传热分析
ANSYS 瞬态传热分析的主要步骤 1.建立有限元模型 2.施加载荷并求解 3.求解 4.查看分析结果
QUST
5
大家好
5. 瞬态传热分析
在瞬态分析中，载荷步和子步的定义与非线性稳态分 析十分类似。载荷定义的每个载荷步的终点，并可以 随时间阶跃或渐进的施加。
每个载荷步的求解是在子步上得到。子步长根据时间 积分步长得到。
自动时间步 (ATS) 同样适用于瞬态分析，可以简化 ITS选择。 ITS选择将影响到瞬态分析的精度和非线性收敛性 (如果存在)。
Restart
大家好
QUST
10条件分为两种情况：其一， 初始温度场已知；其二，初始温度场未知。
瞬态热分析实例(二维)
模型[1]热传导问题：如图，110R cm =，220R cm =，密度为36000/kg m ，比热容为220/()J kg K ⋅，热传导率为6/()W m K ⋅，初始温度为300℃，突然放入30℃的液体中冷却，这种液体对流换热系数2120/()h W m K =⋅。
计算：（1）第1秒和第60秒这两个时刻温度分布情况；（2）内外边在60秒内温度变化。
1．设置环境① 设置分析模块。
本例是温度分布分析，所以只需要选择热分析模块，这样就可以把结构分析模块、电磁场分析模块和流体分析模块的菜单都过滤掉。
设置如图② 设置单位在命令行输入“/units,SI ”，SI 为设定为国际单位制。
必须注意：[1] 秦宇.ANSYS 11.0基础与实例教程[M]. 北京,化学工业出版社,2009:318-330ANSYS程序不会为你的分析假定一个单位制，除了磁场分析以外，你可以使用任何单位制，只要你能保证你输入的所有数据都是按照这个单位制进行的。
也就是说，单位制在所有输入数据中应该保持一致。
使用/UNITS命令，你可以在ANSYS数据库中进行标记来表示你使用的单位制。
但是请注意，这个命令并不将一个单位转化为另一个。
它仅仅只作为对分析的一个评论记录。
什么意思呢？就是/UNITS只是个标记，告诉别人程序的单位制，即使程序中没有使用这种单位制，它也不能将这种单位制转化为自己标记的那个单位制。
所以，如果你要让ANSYS的单位为国际单位制，你在输入物理量之前，先将所有的物理量转换为国际单位制，如：原先你的图纸上均为毫米，比如一个矩形截面尺寸是400mm*500mm，那么，你在建模之前先转化为0.4m*0.5m然后输入的长度为0.4和0.5，ANSYS只知道你输入的是0.4和0.5，它不知道你的单位是什么。
2．定义单元类型和材料属性①选择单元类型。
如图：我们选择【Quad 4node 55】即选择了PLANE55单元。
下面介绍一下PLANE55单元，我们直接从ANSYS帮助文档中摘录。
ansys技巧总结_如何进行盛水水箱的模态分析
d,all,pres,0
alls
asel,u,loc,y,width,
sfa,all,,fsi !定义流固耦合界面
alls
solv
fini
/post1
set,first
plnsol.u,sum,2,1
fini
fini
/solu antype,2
modopt,unsym,10 !非对称模态提取方法处理流固耦合问题
eqslv,front
mxpand,10,,,1
nsel,s,loc,x,
nsel,a,loc,x,length
nsel,r,loc,y
d,all,,,,,,ux,uy,uz,
length=1
width=0.6
height=0.8
/prep7
et,1,63
et,2,30 !选用FLUID30单元，用于流固耦合问题
r,1,0.01
mp,ex,1,2e11
mp,nuxy,1,0.3
mp,dens,1,7800
mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水
mp,sonc,2,1400
mp,mu,2,
!
block,,length,,width,,height
esize,0.1 mFra bibliotekhkey,1
!
type,1
mat,1
real,1
asel,u,loc,y,width
amesh,all
alls
!
type,2
Ansys-第33例瞬态热分析实例一水箱
第33例瞬态热分析实例——水箱本例介绍了利用ANSYS进行瞬态热分析的方法和步骤、瞬态热分析时材料模型所包含的内容，以及模型边界条件和初始温度的施加方法。
33.1概述热分析是计算热应力的基础，热分析分为稳态热分析和瞬态热分析，稳态热分析将在后面两个例子中介绍，本例介绍瞬态热分析。
33.1.1 瞬态热分析的定义瞬态热分析用于计算系统随时间变化的温度场和其他热参数。
一般用瞬态热分析计算温度场，并找到温度梯度最大的时间点，将此时间点的温度场作为热载荷来进行应力计算。
33.1.2 嚼态热分析的步骤瞬态热分析包括建模、施加载荷和求解、查看结果等几个步骤。
1．建模瞬态热分析的建模过程与其他分析相似，包括定义单元类型、定义单元实常数、定义材料特性、建立几何模型和划分网格等。
注意：瞬态热分析必须定义材料的导热系数、密度和比热。
2.施加载荷和求解(1)指定分析类型，Main Menu→Solution→Analysis Type→New Analysis,选择Transient。
(2)获得瞬态热分析的初始条件。
定义均匀的初始温度场：Main Menu→Solution→Define Loads→Settings→Uniform Temp，初始温度仅对第一个子步有效，而用Main Menu→Solution→Define Loads→Apply→Thermal→Temperature命令施加的温度在整个瞬态热分析过程中均不变，应注意二者的区别。
定义非均匀的初始温度场：如果非均匀的初始温度场是已知的，可以用Main Menu→Solution→Define Loads →Apply→Initial Condit'n→Define即IC命令施加。
非均匀的初始温度场一般是未知的，此时必须先进行行稳态分析确定该温度场。
该稳态分析与一般的稳态分析相同。
注意：要设定载荷（如已知的温度、热对流等），将时间积分关闭，选择Main Menu→Solution→Load Step Opts →Time/Frequenc→Time Integration→Amplitude Decay;设定只有一个子步，时间很短（如(0.01s)的载荷步，Main Menu→Solution→Load Step Opts→Time/Frequenc→Time→Time Step。
ansys热分析实例教程
Temperature distribution in a CylinderWe wish to compute the temperature distribution in a long steel cylinder with inner radius 5 inches and outer radius 10 inches. The interior of the cylinder is kept at 75 deg F, and heatis lost on the exterior by convection to a fluid whose temperature is 40 deg F. The convection coefficient is 0.56 BTU/hr-sq.in-F and the thermal conductivity for steel is 0.69 BTU/hr-in-F.1. Start ANSYS and assign a job name to the project. Run Interactive -> set working directory and jobname.2. Preferences -> Thermal will show -> OK3. Recognize symmetry of the problem, and a quadrant of a section through the cylinder is created using ANSYS area creation tools. Preprocessor -> Modeling -> Create -> Areas -> Circle -> Partial annulusThe following geometry is created.4. Preprocessor -> Element Type -> Add/Edit/Delete -> Add -> Thermal Solid -> Solid 8 node 77 -> OK -> Close5. Preprocessor -> Material Props -> Isotropic -> Material Number 1 -> OKEX = 3.E7 (psi)DENS = 7.36E-4 (lb sec^2/in^4)ALPHAX = 6.5E-6PRXY = 0.3KXX = 0.69 (BTU/hr-in-F)6. Mesh the area and refine using methods discussed in previous examples.7. Preprocessor -> Loads -> Apply -> Temperatures -> NodesSelect the nodes on the interior and set the temperature to 75.8. Preprocessor -> Loads -> Apply -> Convection -> LinesSelect the lines defining the outer surface and set the convection coefficient to 0.56 and the fluid temp to 40.9. Preprocessor -> Loads -> Apply -> Heat Flux -> LinesTo account for symmetry, select the vertical and horizontal lines of symmetry and set the heat flux to zero.10. Solution -> Solve current LS11. General Postprocessor -> Plot Results -> Nodal Solution -> TemperaturesThe temperature on the interior is 75 F and on the outside wall it is found to be 45. These results can be checked using results from heat transfer theory.BackThermal Stress of a Cylinder using Axisymmetric ElementsA steel cylinder with inner radius 5 inches and outer radius 10 inches is 40 inches long and has spherical end caps. The interior of the cylinder is kept at 75 deg F, and heat is lost on the exterior by convection to a fluid whose temperature is 40 deg F. The convection coefficient is 0.56 BTU/hr-sq.in-F. Calculate the stresses in the cylinder caused by the temperature distribution.The problem is solved in two steps. First, the geometry is created, the preference set to'thermal', and the heat transfer problem is modeled and solved. The results of the heat transfer analysis are saved in a file 'jobname.RTH' (Results THermal analysis) when you issue a save jobname.db command.Next the heat transfer boundary conditions and loads are removed from the mesh, the preference is changed to 'structural', the element type is changed from 'thermal' to 'structural', and the temperatures saved in 'jobname.RTH' are recalled and applied as loads.1. Start ANSYS and assign a job name to the project. Run Interactive -> set working directory and jobname.2. Preferences -> Thermal will show -> OK3. A quadrant of a section through the cylinder is created using ANSYS area creation tools.4. Preprocessor -> Element Type -> Add/Edit/Delete -> Add -> Solid 8 node 77 -> OK ->Options -> K3 Axisymmetric -> OK5. Preprocessor -> Material Props -> Isotropic -> Material Number 1 -> OKEX = 3.E7 (psi)DENS = 7.36E-4 (lb sec^2/in^4)ALPHAX = 6.5E-6PRXY = 0.3KXX = 0.69 (BTU/hr-in-F)6. Mesh the area using methods discussed in previous examples.7. Preprocessor -> Loads -> Apply -> Temperatures -> NodesSelect the nodes on the interior and set the temperature to 75.8. Preprocessor -> Loads -> Apply -> Convection -> LinesSelect the lines defining the outer surface and set the coefficient to 0.56 and the fluid temp to 40.9. Preprocessor -> Loads -> Apply -> Heat Flux -> LinesSelect the vertical and horizontal lines of symmetry and set the heat flux to zero.10. Solution -> Solve current LS11. General Postprocessor -> Plot Results -> Nodal Solution -> TemperatureThe temperature on the interior is 75 F and on the outside wall it is found to be 43.12. File -> Save Jobname.db13. Preprocessor -> Loads -> Delete -> Delete All -> Delete All Opts.14. Preferences -> Structural will show, Thermal will NOT show.15. Preprocessor -> Element Type -> Switch Element Type -> OK (This changes the element to structural)16. Preprocessor -> Loads -> Apply -> Displacements -> Nodes(Fix nodes on vertical and horizontal lines of symmetry from crossing the lines of symmetry.)17. Preprocessor -> Loads -> Apply -> Temperature -> From Thermal AnalysisSelect Jobname.RTH (If it isn't present, look for the default 'file.RTH' in the root directory)18. Solution -> Solve Current LS19. General Postprocessor -> Plot Results -> Element Solution - von Mises StressThe von Mises stress is seen to be a maximum in the end cap on the interior of the cylinder and would govern a yield-based design decision.Back。
ANSYS瞬态热分析--样件加热
ANSYS瞬态热分析--样件加热⽬录1. 要求⼀块0.8m*0.4m*0.2m厚的钢板，在300℃的炉⼦中进⾏加热，其材料参数如下表所⽰:材料参数序号名称参数1导热系数452弹性模量210E93⽐热5403泊松⽐0.34密度78502. ANSYS有限元分析ANSYS单位W/m℃对流换热系数：取202.1 APDL建模!程序头FINISH/CLEAR/FILENAME, TSR/PREP7!设置单元ET,1,SOLID278!设置材料MP,KXX,1,45!定义材料的导热系数MP,C,1,540!定义⽐热容MP,DENS,1,7850!密度!建⽴模型BLC4,0,0,0.8,0.4,0.2!设置属性TYPE,1MAT,1!设置⽹格划分LESIZE,1,,,25LESIZE,2,,,25LESIZE,9,,,20VMESH,ALL!设置初始温度TUNIF,22,!对流换热系数和环境温度SFA,ALL,1,CONV,100,300!设置温度偏移TOFFST,273,!设置求解FINISH/SOLANTYPE,4!瞬态分析KBC,1TRNOPT,FULL!完全法LUMPM,0DELTIM,60,30,300!默认时间间隔60s，最⼩30s，最⼤300s OUTRES,ALL,ALL!结果完全输出TIME,6000!100分钟SOLVE!求解FINISH/POST1PLNSOL,TEMP,,0!显⽰温度云图!显⽰中间⼀点的温度随时间的曲线FINISH/POST26SOLU,191,NCMITSTORE,MERGEFILLDATA,191,,,,1,1REALVAR,191,191NSOL,2,NODE(0.4,0.2,0.2),TEMP,,TEMP_3STORE,MERGEXVAR,1PLVAR,2,SAVE最后时刻温度云图温度随时间变化曲线。
Ansys作业——瞬态热分析研究
Ansys作业—瞬态热分析问题描述瞬态热分析实例1长方形地板，几何参数及其边界条件如图3-6 所示.板地宽度为5cm，其中间有一个半径为1cm 地圆孔.板地初始温度为20℃，将其右侧突然置于温度为20℃且对流换热系数为100W/M2℃地流体中，左端置于温度为500℃地温度场，试计算：⏹（1）第1s 和第50s板内地温度分布情况.⏹（2）整个板在前50s内地温度变化过程.⏹（3）圆孔边缘A点处温度随时间变化曲线.1.建立有限元模型首先建立瞬态传热分析所需地有限元模型，选择单元.(1) 选择热分析单元，操作如下：GUI：Main Menu > Preprocessor > Element Type > Add/Edit/Delete在弹出地对话框中，单击Add.在单元类型库对话框中选择Plane55单元.单击OK.命令：ET,1,PLANE55(2) 定义材料属性首先进入Define Material Model Behavior对话框，操作如下：GUI：Main Menu > Preprocessor > Material Props下面定义瞬态热分析所需地材料参数，如热传导率、比热容及材料密度：定义热传导GUI：Main Menu > Preprocessor > Material Props > Thermal > Conductivity > Isotropic在弹出地定义材料热传导率对话框中地KXX 栏键入“5”.命令：MPDATA,KXX,1,,5定义比热容GUI：Main Menu > Preprocessor > Material Props > Thermal > Specific Heat在弹出地定义比热容对话框中地C栏键入“200”.命令：MPDATA,C,1,,200定义密度GUI：Main Menu > Preprocessor > Material Props > Thermal > Density在弹出密度定义对话框中地DENS栏键入“5000”.命令：MPDATA,DENS,1,,5000材料属性定义完毕.(3)建立实体模型根据本例所用模型，首先需要创建矩形，然后是圆，最后在矩形板中央减去（Substract）圆.下面介绍建立实体模型地操作：创建矩形命令：RECTNG,0,0.15,0,0.05创建圆面其操作如下：GUI：Main Menu > Preprocessor > Modeling > Create > Circle > By dimensions 在弹出对话框中，单击OK得到圆面.命令：CYL4, 0.075, 0.025, 0.01相减根据ANSYS建模过程中面序号赋值原理，直接可以肯定圆面序号为2，矩形序号为1，因此采用直接键入命令建实体模型：命令：asba,1,2图1 实体模型(4) 设定网格尺寸并划分网格图2 有限元模型2.施加载荷并求解(1)选择分析类型GUI：Main Menu > Preprocessor > Loads > Analysis Type > New Analysis 选择Transient 分析，单击OK.在弹出地子对话框中单击OK.(2)定义初始条件板地初始温度为20℃，设置初始温度操作如下：GUI：Main Menu > Solution > Define Loads > Apply > Initial Condit'n > Define在弹出地拾取对话框中，单击Pick All.弹出Define Initial Conditions对话框，命令：IC,All,TEMP,20,(3)定义热约束定义对流边界GUI：Utility Menu > Plot > lines命令：LPLOT定义对流载荷/边界首先进入Apply Conv on lines 对话框，操作如下：GUI：Main Menu > Solution > Define Loads > Apply > Thermal > Convection > On Lines在弹出对话框中，键入Film coefficient和Bulk Temperature值.命令：SFL,L2,CONV,100, ,20,定义稳态热边界在边线上定义稳态热边界，操作如下：GUI：Main Menu > Solution > Define Loads > Apply > Thermal > Temperature > On Lines在弹出对话框中，键入边界温度为“500”.单击OK.命令：DL,L4, ,TEMP,500,1图3施加边界条件和载荷(4)设置时间及时间步进参数GUI：Main Menu > Solution > Load Step Opts > Time/Frequenc > Time –time steps在Time and time Step Optios对话框中按照下图设置，其余选项不变.TIME,50AUTOTS,1DELTIM,1,0.1,2.5,1KBC,0下面确认时间积分选项打开.操作如下：GUI：Main Menu > Solution > Load Step Opts > Time/Frequenc > Time Integration保持默认设置，单击OK.命令：TIMINT,1TINTP,0.005, , ,-1,0.5,-1(5)设置输出控制GUI：Main Menu > Solution > Load Step Opts > Output Ctrls > DB/Results File 在File Write Frequency栏选项 Every Substep.命令：OUTRES,ALL,ALL,(6)求解GUI：Main Menu > Solution > Solve > Current LS.命令：Solve3.查看分析结果（1）POST1后处理，显示第1s和和第50s时温度场.图4显示第1s时地温度场重复上述操作，读入第50s时结果并显示温度场.下图所示为第50s温度场.图5显示第50s时地温度场(2)POST26 后处理需要查看A点随时间变化地温度场，首先必须定义A点温度变量，然后才能显示A点处.定义A点温度变量显示A 点温度随时间变化曲线，单击“GRAPH VARIABLE 按钮，得到A 点温度随时间变化曲线.图6A点温度-时间变化曲线4.命令流文件!进入前处理/prep7 /COM,Thermal Et,1,plane55Mp,dens,1,5000Mp,kxx,1,5Mp,c,1,200Save!创建几何模型Rectng,0,0.15,0,0.05CYL4,0.075,0.025,0.01Asba,1,2Save!划分网格SMRT,3Amesh,ALLSave!进入加载求解/soluAntype,transTRNOT,FULLLUMPM,0Ic,all,temp,20SaveLplotSFL,2,conv,100,,20 DL,4,,TEMP,500,0Time,50Kbc,0Autots,onDeltim,1,0.1,2.5,1Timint,onTINTP,0.005,,,-1,0.5,-1Outres,all,allSolveFINISH/POST1!*/EFACET,1PLNSOL, TEMP,, 0瞬态传热分析实例21.问题描述：一个直径为0.12m，温度为900摄氏度地钢球突然放入盛满了水地，完全绝热地横截面直径和孤傲度均为0.6m地圆柱体水箱中，水地温度为20，材料参数如图所示，求解10分钟后钢球与水地温度场分布（钢球置于水箱正中央）.2. 问题分析属于瞬间热传导问题，研究对象为钢球和水，根据轴对称性，在求解过程中取钢球和水中心纵截面地1/4建立几何模型，选择PLANE55轴对称单元进行分析.3.建立模型(1)添加标题添加标题，操作如下GUI：Utility Menu > File > Change Title在弹出地Referenc for GUI Filtering对话框中，选择 Thermal.单击OK.(2) 选择单元在单元类型库对话框中选择Plane55单元.单击OK.在Element behavior K3下拉框中选择Axisymmetric选项，其余默认.(3)定义材料属性1⏹定义热传导率在弹出地定义材料热传导率对话框中地KXX 栏键入“70”.⏹定义比热容在弹出地定义比热容对话框中地C栏键入“448”⏹定义密度在弹出密度定义对话框中地DENS栏键入“7800”.(4)定义材料属性2⏹定义热传导率在弹出地定义材料热传导率对话框中地KXX 栏键入“0.61”.⏹定义比热容在弹出地定义比热容对话框中地C栏键入“4185”⏹定义密度在弹出密度定义对话框中地DENS栏键入“1000”.(5)建立实体模型图7平面实体模型4.设定网格尺寸并划分网格图8有限元模型5．加载及求解求解之前首先要选择分析类型，然后定义边界条件及其载荷步选项，最后计算.首先选择分析类型.图9求解后显示图6.查看求解结果在POST1后处理器中查看结果图10总体温度场等值线图图11水地温度场等值线图图12球体温度场等值线图图13球心温度随时间地变化关系图7.命令流文件/Title,Transient slab ProblemKEYW, PR_THERM, !指定分析类型为热分析/prep7ET,1,Plane55KEYOPT,1,3,1MP,KXX,1,70MP,DENS,1,7800MP,C,1,448MP,KXX,2,0.61MP,DENS,2,1000MP,C,2,4185RECTANG,0,0.3,O,0.3PCIRC,0.06,0,0,90AOVERLAP,ALLNUMCMP,AREANUMCMP,LINE/TITLE,GEOMETRIC MODELAPLOTMSHAPE,0,2DMSHKEY,1LPLOT/PNUM,LINE,1/PNUM,AREA,1/REPLOTLSEL,S,LINE,,4,5LESIZE,ALL,,,30,0.1LSEL,S,LINE,,6,7 !选择线段6，7LESIZE,ALL,,,32,0.1 !设置单元数量LSEL,S,LINE,,3 !选择线段3LESIZE,ALL,,,30 !设置单元数量ALLSELLCCAT,2,1 !将线段12联接在一起MAT,1 !选择材料参考号1AMESH,1 !对面1进行网格划分MAT,2 !选择材料参考号1AMESH,2 !对面1进行网格划分/TITLE, ELEMENTS IN MODELEPLOTFINISH/SOLUANTYPE,TRANSTIMINT,OFFTIME,0.01DELTIM,0.01ESEL,S,MAT,,1NSLE,SD,ALL,TEMP,900ESEL,S,MAT,,2NSLE,SD,ALL,TEMP,20ALLSELSOLVETIMINT,ONTIME,600DELTIM,26,2,200AUTOS,ONDDELE,ALL,TEMPOUTRES,,ALLSOLVEFINISH版权申明本文部分内容，包括文字、图片、以及设计等在网上搜集整理.版权为个人所有This article includes some parts, including text, pictures, and design. Copyright is personal ownership.GMsIa。
Ansys-第36例--热应力分析(间接法)实例—液体管路
第36例热应力分析（间接法）实例一液体管路本例介绍了利用间接法进行热应力计算的方法和步骤：首先进行热分析得到结构节点温度分布，然后把温度作为载荷施加到结构上并进行结构分析。
36.1概述利用间接法计算热应力，首先进行热分析，然后进行结构分析。
热分析可以是瞬态的，也可以是稳态的，需要将热分析求得的节点温度作为体载荷施加到结构上。
当热分析是瞬态的时，需要找到温度梯度最大的时间点，并将该时间点的结构温度场作为体载荷施加到结构上。
由于•间接法可以使用所有热分析和结构分析的功能，所以对「•大多数情况都推荐使用该方法。
间接法进行热应力计算的主要步骤如下。
热分析瞬态热分析的过程在前例已经介绍过，下面介绍稳态热分析。
稳态热分析用于研究稳定的热载荷对结构的影响，有时还用于瞬态热分析时计算初始温度场。
稳态热分析主要步骤如bo1.建模稳态热分析的建模过程与其他分析相似，包括定义单元类型、定义单元实常数、定义材料特性、建立几何模型和划分网格等。
但需注意的是：稳态热分析必须定义材料的导热系数。
2.施加载荷和求解⑴指定分析类型。
Main Menu-*Solution-*Analysis Type~*New Analysis,选择Static.⑵施加载荷。
nJ■以施加的载荷有恒定的温度、热流率、对流、热流密度、生热率，Main Menu-*Solution-*Define Loadsf Apply—Thermal.（3）设置载荷步选项。
普通选项包括时间〔用于定义载荷步和子步）、每一载荷步的子步数，以及阶跃选项等，Main Menu-*Solution-*Load Step Opts—k Time/Frequenc->Time->Time Step.非线性选项包括:一迭代次数（默认25）, Main Menu-*Solution-*Load Step Opts-* Nonlinear-* Equilibrium Iter;翻开自动时间步长，Main Menu-*Solution-* Load Step Opts —Time/Frequenc—Time—Time Step 等.图36-11转换单元类型对话框设定单元轴对称选项拾取菜单Main Menu —Preprocessor—Element Type —Add/Edit/Delete 弹出“ElementTypes”对话框，单击其“Options”按钮，弹出如图36-12所示的对话框，选择“K3”下拉列表框为“Axisymmetric",单击“0K"按钮，然后单击"Element Types M对话框的“Close"按钮。
用 ANSYS 模拟瞬态热传导问题
用ANSYS 模拟瞬态热传导问题问题描述：一个30 公斤重、温度为70℃的铜块,以及一个20 公斤重、温度为80℃的铁块，突然放入温度为20℃、盛满了300 升水的、完全绝热的水箱中。
过了一个小时，求铜块与铁块的最高温度(假设忽略水的流动)。
材料热物理性能如下:菜单操作:1、Main Menu→Preprocessor，进入前处理;择单元plane77;3、定义材料Main Menu→Preprocessor→Material Props→Thermal→Conductivity Isotropic定义材料1(铜)的KXX 等于383、DENS 等于8889、C 等于1390;同理，Material→New Model，按上述方法添加铁和水定义材料2(铁)的KXX 等于70、DENS 等于7837、C 等于448;定义材料3(水)的KXX 等于0.61、DENS 等于996、C 等于41854、建立几何模型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;进行布尔操作Main Menu→Preprocessor→Modeling→Operate→Booleans→Overlap, areas，选择Pick All;Glue, areas，选择Pick All;5、Utility Menu: Plotctrls→Numbering→Areas, on;显示不同面积区域。
6、网格划分Main Menu→Preprocessor→Mesh tool→Set，Areas→Picked Areas→选择材料1(铜) 选Mat 1，Apply 。
ANSYS瞬态热分析---零件在水中冷却
ANSYS瞬态热分析---零件在⽔中冷却⽬录1. 案例⼀个温度为300℃的铜环和⼀个温度为200℃的铁环，放置到22℃的⽔中进⾏淬⽕。
⽔桶为铁质的圆形。
分析中忽略⽔的流动。
材料参数热性能铜铁⽔导热系数（W/(m℃）383700.61密度（kg/m388897833996⽐热容（J/(kg.℃）3904484185因为忽略了⽔的流动，故可将⽔看为⼀种材料，直接进⾏分析2. APDL分析finish/clear!==============第⼀步：进⾏稳态⼆维轴对称分析===============================/filname,quench/prep7et,1,plane55 ! 选择单元类型keyopt,1,3,1 ! 设置为轴对称分析! 定义材料热性能参数mp,kxx,1,383mp,dens,1,8889mp,c,1,390 ! 定义铜的热物理性能参数mp,kxx,2,70mp,dens,2,7837mp,c,2,448 ! 定义铁的热物理性能参数mp,kxx,3,0.61mp,dens,3,996mp,c,3,4185 ! 定义⽔的热物理性能参数rectng,0,0.08,0,0.01rectng,0.08,0.1,0,0.01rectng,0.1,0.12,0,0.01rectng,0.12,0.14,0,0.01rectng,0.14,0.15,0,0.01rectng,0.14,0.15,0.01,0.055rectng,0.14,0.15,0.055,0.1rectng,0.14,0.15,0.1,0.15 ! 建⽴铁箱矩形rectng,0.08,0.1,0.01,0.055rectng,0.1,0.12,0.01,0.055 ! 建⽴铜环rectng,0.08,0.1,0.055,0.1 ! 建⽴铁环rectng,0,0.08,0.01,0.055rectng,0.12,0.14,0.01,0.055rectng,0,0.08,0.055,0.1rectng,0.10,0.12,0.055,0.1rectng,0.12,0.14,0.055,0.1rectng,0,0.08,0.1,0.15rectng,0.08,0.1,0.1,0.15rectng,0.10,0.12,0.1,0.15rectng,0.12,0.14,0.1,0.15 ! 建⽴⽔⼏何模型aglue,all ! 粘接各矩形esize,0.003,0, ! 定义单元划分尺⼨asel,s,,,1,21,20asel,a,,,23,29,1aatt,2,1,1 ! 附于铁箱和铁环的材料属性asel,s,,,30,33,3aatt,1,1,1 ! 附于铜环的材料属性asel,s,,,31,32,1asel,a,,,34,39,1asel,a,,,22aatt,3,1,1 ! 附于⽔的材料属性allsel,allmshkey,2 ! 定义映射单元划分类型amesh,all ! 划分单元asel,s,,,1,21,20asel,a,,,23,28,1asel,a,,,31,32,1asel,a,,,34,39,1asel,a,,,22nsla,,1d,all,temp,20 ! 施加铁箱和⽔的温度约束条件asel,s,,,30,33,3nsla,,1d,all,temp,300 ! 施加铜环的温度约束条件asel,s,,,29nsla,,1d,all,temp,200 ! 施加铁环的温度约束条件allsel,all/soluantype,trans ! 设定为瞬态分析timint,off ! 关掉时间积分, 即设定为稳态分析time,0.01 ! 设定只有⼀个⼦步的时间很⼩的载荷步deltim,0.01 ! 设定求解时间⼦步solve ! 求解!==============第⼆步：进⾏瞬态⼆维轴对称分析=============================== time,3600 ! 设定求解时间timint,on ! 打开时间积分deltim,26,2,200 ! 设置时间步长，最⼤及最⼩时间步长autots,on ! 打开⾃动时间步长ddelet,all,temp ! 删除稳态分析中定义的节点温度outres,all,1 ! 将每个⼦步的值写⼊数据库⽂件solve ! 求解finish/post1plnsol,temp,,0样件上的某点随时间变化的曲线。
Ansys作业——瞬态热分析
Ansys作业—瞬态热分析问题描述瞬态热分析实例1⏹长方形的板，几何参数及其边界条件如图3-6 所示。
板的宽度为5cm，其中间有一个半径为1cm 的圆孔。
板的初始温度为20℃，将其右侧突然置于温度为20℃且对流换热系数为100W/M2℃的流体中，左端置于温度为500℃的温度场，试计算：⏹（1）第1s 和第50s板内的温度分布情况。
⏹（2）整个板在前50s内的温度变化过程。
⏹（3）圆孔边缘A点处温度随时间变化曲线。
1.建立有限元模型首先建立瞬态传热分析所需的有限元模型，选择单元。
(1) 选择热分析单元，操作如下：GUI：Main Menu > Preprocessor > Element Type > Add/Edit/Delete在弹出的对话框中，单击Add。
在单元类型库对话框中选择Plane55单元。
单击OK。
命令：ET,1,PLANE55(2) 定义材料属性首先进入Define Material Model Behavior对话框，操作如下：GUI：Main Menu > Preprocessor > Material Props下面定义瞬态热分析所需的材料参数，如热传导率、比热容及材料密度：定义热传导GUI：Main Menu > Preprocessor > Material Props > Thermal > Conductivity > Isotropic 在弹出的定义材料热传导率对话框中的KXX 栏键入“5”。
命令：MPDATA,KXX,1,,5定义比热容GUI：Main Menu > Preprocessor > Material Props > Thermal > Specific Heat在弹出的定义比热容对话框中的C栏键入“200”。
命令：MPDATA,C,1,,200定义密度GUI：Main Menu > Preprocessor > Material Props > Thermal > Density在弹出密度定义对话框中的DENS栏键入“5000”。
ANSYS热分析指南——ANSYS瞬态热分析24页word
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 Menu>Preferences中对菜单进行过滤。
瞬态传热分析实例
瞬态传热分析实例 ansys
h
1
4. 瞬态传热分析
瞬态传热分析的基本步骤与稳态热分析类似，主要 的区别是瞬态传热分析中的载荷是随时间变化的。 为了表达随时间变化的载荷，首先必须将载荷－时 间曲线分为载荷步。载荷－时间曲线中的每一个拐 点为一个载荷步，如图3-9 所示。对于每一个载荷 步，必须定义载荷值及时间值，同时必须选择载荷 步为Ramped 方式变化或Stepped 方式变化。
1. 问题描述
实例类型：ANSYS结构分析。 分析类型：瞬态传热分析。
单元类型：PLANE55 ANSYS功能示例：实体建模包括基本的建模操
作；定义比热容；施加瞬态热载荷；设置瞬态热 载荷分析选项；显示模型温度等值线图；显示节 点温度随时间变化曲线。
h
23
h
24
长方形的板，几何参数及其边界条件如图3-6 所示。板的宽度为5cm，其中间有一个半径 为1cm 的圆孔。板的初始温度为20℃，将其 右侧突然置于温度为20℃且对流换热系数为 100W/M2℃的流体中，左端置于温度为 500℃的温度场，试计算：
Apply > Structural > Temperature > Uniform Temp
如果已知模型的起始温度是均匀的，可设定所有节 点初始温度。
h
12
2. 设定参考温度 命令：TREF GUI：Main Menu > Solution > Define Loads
> Settings > Reference Temp
Time/Frequenc > Time Integration
(瞬态)ansys热分析例题2.doc
问题描述：一个30公斤重、温度为70℃的铜块,以及一个20公斤重、温度为80℃的铁块，突然放入温度为20℃、盛满了300升水的、完全绝热的水箱中,如图所示。
过了一个小时，求铜块与铁块的最高温度(假设忽略水的流动)。
材料热物理性能如下：热性能单位制铜铁水导热系数W/m℃383 37 0.61密度Kg/m 8889 7833 996比热J/kg℃390 448 4185菜单操作过程：一、设置分析标题1、选择“Utility Menu>File>Change Jobname”，输入文件名Transient1。
2、选择“Utility Menu>File>Change Title”输入Thermal Transient Exercise 1。
二、定义单元类型1、选择“Main Menu>Preprocessor”，进入前处理。
2、选择“Main Menu>Preprocesor>Element Type>Add/Edit/Delete”。
选择热平面单元plane77。
三、定义材料属性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，定义第三种材料。
四个ANSYS热分析经典例子
实例1：某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。
几何参数：筒外径30 feet总壁厚2 inch不锈钢层壁厚0.75inch玻纤层壁厚1 inch铝层壁厚0.25inch筒长200 feet导热系数不锈钢8.27BTU/hr.ft.o F玻纤0.028 BTU/hr.ft.o F铝117.4 BTU/hr.ft.o F边界条件空气温度70 o F海水温度44.5 o F空气对流系数2.5 BTU/hr.ft2.o F海水对流系数80 BTU/hr.ft2.o F沿垂直于圆筒轴线作横截面，得到一圆环，取其中1 度进行分析，如图示。
/filename，Steady1/title，Steady-state thermal analysis of submarine/units，BFTRo=15 !外径(ft)Rss=15-(0.75/12) !不锈钢层内径ft)Rins=15-(1.75/12) !玻璃纤维层内径(ft)Ral=15-(2/12) !铝层内径(ft)Tair=70 !潜水艇内空气温度Tsea=44.5 !海水温度Kss=8.27 !不锈钢的导热系数(BTU/hr.ft.oF)Kins=0.028 !玻璃纤维的导热系数(BTU/hr.ft.oF)Kal=117.4 !铝的导热系数(BTU/hr.ft.oF)Hair=2.5 !空气的对流系数(BTU/hr.ft2.oF)Hsea=80 !海水的对流系数(BTU/hr.ft2.oF)prep7et，1，plane55 !定义二维热单元mp，kxx，1，Kss !设定不锈钢的导热系数mp，kxx，2，Kins !设定玻璃纤维的导热系数mp，kxx，3，Kal !设定铝的导热系数pcirc，Ro，Rss，-0.5，0.5 !创建几何模型pcirc，Rss，Rins，-0.5，0.5pcirc，Rins，Ral，-0.5，0.5aglue，allnumcmp，arealesize，1，，，16 !设定划分网格密度lesize，4，，，4lesize，14，，，5lesize，16，，，2Mshape，2 !设定为映射网格划分mat，1amesh，1mat，2amesh，2mat，3amesh，3/SOLUSFL，11，CONV，HAIR，，TAIR !施加空气对流边界SFL，1，CONV，HSEA，，TSEA !施加海水对流边界SOLVE/POST1PLNSOL ！输出温度彩色云图finish实例2一圆筒形的罐有一接管，罐外径为3 英尺，壁厚为0.2 英尺，接管外径为0.5 英尺，壁厚为0.1英尺，罐与接管的轴线垂直且接管远离罐的端部。
一个经典的ansys热分析实例(流程序)
/PREP7/TITLE,Steady-state thermal analysis of pipe junction/UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches)! /SHOW, ! Specify graphics driver for interactive runET,1,90 ! Define 20-node, 3-D thermal solid elementMP,DENS,1,.285 ! Density = .285 lbf/in^3MPTEMP,,70,200,300,400,500 ! Create temperature tableMPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12! 指定与温度相对应的数据材料属性；导热系数；Define conductivity valuesMPDATA,C,1,,.113,.117,.119,.122,.125! Define specific heat values（比热）MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144! Define film coefficient;除144是单位问题，上面的除12也是单元问题! Define parameters for model generationRI1=1.3 ! Inside radius of cylindrical tankRO1=1.5 ! Outside radiusZ1=2 ! LengthRI2=.4 ! Inside radius of pipeRO2=.5 ! Outside pipe radiusZ2=2 ! Pipe lengthCYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tankWPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axisCYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipeWPSTYL,DEFA ! 重新安排工作平面的设置；另外WPSTYL,STAT to list the status of the working plane；；Return working plane to default settingBOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warningVOVLAP,1,2 ! 交迭体；Overlap the two cylinders/PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on/VIEW,,-3,-1,1/TYPE,,4 ! 精确面的显示；Precise hidden display/TITLE,Volumes used in building pipe/tank junctionVPLOTVDELE,3,4,,1 ! 修剪一些体与体相关的体的因素都删掉；Trim off excess volumes! Meshing 网格划分ASEL,,LOC,Z,Z1 ! Select area at remote Z edge of tankASEL,A,LOC,Y,0 ! Select area at remote Y edge of tankCM,AREMOTE,AREA ! 为面建立数组；Create area component called AREMOTE/PNUM,AREA,1/PNUM,LINE,1/TITLE,Lines showing the portion being modeledAPLOT/NOERASE ! 预防抹去LPLOT ! Overlay line plot on area plot/ERASEACCAT,ALL ! 连接面和线的准备映射；Concatenate areas and lines at remote tank edgesLCCAT,12,7LCCAT,10,5LESIZE,20,,,4 ! 4 divisions through pipe thicknessLESIZE,40,,,6 ! 6 divisions along pipe lengthLESIZE,6,,,4 ! 4 divisions through tank thicknessALLSEL ! Restore full set of entitiesESIZE,.4 ! Set default element size线的默认划分数MSHAPE,0,3D ! Choose mapped brick meshMSHKEY,1 ! 映射网格SAVE ! Save database before meshingVMESH,ALL ! Generate nodes and elements within volumes/PNUM,DEFA ! 重新安排数字规格/TITLE,Elements in portion being modeledEPLOTFINISH/COM, *** Obtain solution ***/SOLUANTYPE,STATIC ! Steady-state analysis typeNROPT,AUTO ! 自动选择牛顿-拉普森Program-chosenNewton-Raphson optionTUNIF,450 ! 给结点统一的温度；Uniform starting temperature at all nodesCSYS,1 ! 1 —Cylindrical with Z as the axis of rotation NSEL,S,LOC,X,RI1 ! Nodes on inner tank surfaceSF,ALL,CONV,250/144,450 ! 为结点指定表面载荷；对流；Convection（对流）；load at all nodesCMSEL,,AREMOTE ! 选择子集组合；Select AREMOTE component NSLA,,1 ! Nodes belonging to AREMOTED,ALL,TEMP,450 ! 设定边界温度条件Temperature constraints at those nodesWPROTA,0,-90 ! Rotate working plane to pipe axisCSWPLA,11,1 ! 在工作区声明本地的圆柱体系；Define local cylindrical c.s at working planeNSEL,S,LOC,X,RI2 ! Nodes on inner pipe surfaceSF,ALL,CONV,-2,100 ! 这里的-2表示材料2；；Temperature-dep. convection load at those nodesALLSEL/PBC,TEMP,,1 ! 边界符号的显示Temperature b.c. symbols on/PSF,CONV,,2 ! Convection symbols on 箭头显示/TITLE,Boundary conditionsNPLOTWPSTYL,DEFACSYS,0AUTOTS,ON ! Automatic time steppingNSUBST,50 ! Number of substepsKBC,0 ! Ramped loading (default)OUTPR,NSOL,LAST ! 显示最后一次的结点约束；Optional command for solution printoutSOLVEFINISH/COM, *** Review results ***/POST1/EDGE,,1 ! Displays only the "edges（刀口, 利刃, 锋, 优势, 边缘, 优势, 尖锐）" of an object；Edge display/PLOPTS,INFO,ON ! Legend column on/PLOPTS,LEG1,OFF ! Legend header off 圆柱数列的头部/WINDOW,1,SQUARE ! SQUA, form largest square window within the current graphics area；Redefine window size/TITLE,Temperature contours at pipe/tank junctionPLNSOL,TEMP ! Plot temperature contoursCSYS,11NSEL,,LOC,X,RO2 ! Nodes and elements at outer radius of pipeESLN ! 选择单元NSLE ! 选择结点/SHOW,,,1 ! 向量显示；Vector mode/TITLE,Thermal flux vectors at pipe/tank junctionPLVECT,TF ! Plot thermal flux（热通量）vectorsFINISH。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第33例瞬态热分析实例——水箱本例介绍了利用ANSYS进行瞬态热分析的方法和步骤、瞬态热分析时材料模型所包含的内容，以及模型边界条件和初始温度的施加方法。
33.1概述热分析是计算热应力的基础，热分析分为稳态热分析和瞬态热分析，稳态热分析将在后面两个例子中介绍，本例介绍瞬态热分析。
33.1.1 瞬态热分析的定义瞬态热分析用于计算系统随时间变化的温度场和其他热参数。
一般用瞬态热分析计算温度场，并找到温度梯度最大的时间点，将此时间点的温度场作为热载荷来进行应力计算。
33.1.2 嚼态热分析的步骤瞬态热分析包括建模、施加载荷和求解、查看结果等几个步骤。
1．建模瞬态热分析的建模过程与其他分析相似，包括定义单元类型、定义单元实常数、定义材料特性、建立几何模型和划分网格等。
注意：瞬态热分析必须定义材料的导热系数、密度和比热。
2.施加载荷和求解(1)指定分析类型，Main Menu→Solution→Analysis Type→New Analysis,选择 Transient。
(2)获得瞬态热分析的初始条件。
定义均匀的初始温度场：Main Menu→Solution→Define Loads→Settings→Uniform Temp，初始温度仅对第一个子步有效，而用Main Menu→Solution→Define Loads→Apply→Thermal→Temperature命令施加的温度在整个瞬态热分析过程中均不变，应注意二者的区别。
定义非均匀的初始温度场：如果非均匀的初始温度场是已知的，可以用Main Menu →Solution→Define Loads→Apply→Initial Condit'n→Define即IC命令施加。
非均匀的初始温度场一般是未知的，此时必须先进行行稳态分析确定该温度场。
该稳态分析与一般的稳态分析相同。
注意：要设定载荷（如已知的温度、热对流等），将时间积分关闭，选择Main Menu →Solution→Load Step Opts→Time/Frequenc→Time Integration→Amplitude Decay;设定只有一个子步，时间很短（如(0.01s)的载荷步， Main Menu→Solution →Load Step Opts→Time/Frequenc→Time→Time Step。
(3)设置载荷步选项。
普通选项包括每一载荷步结束的时间、每一载荷步的子步数、阶跃选项等，选择Main Menu→Solution→Load Step Opts→Time/Frequenc→Time-Time Step.非线性选项包括：迭代次数（默认25），选择Main Menu→Solution→Load StepOpts→Nonlinear→Equilibrium Iter;打开自动时间步长，选择Main Menu→Solution→Load Step Opts→Time/Frequenc→Time→Time Step：将时间积分打开，选择Main Menu→Solution→Load Step Opts→Time/Frequenc→Time Integration →Amplitude Decay.输出选项包括：控制打印的输出，选择Main Menu→Solution→Load Step Opts →Output Ctrls→Solu Printout; 结果文件的输出，选择Main Menu→Solution →Load Step Opts→Output Ctrls→DB/Results File.(4)如果均匀的初始温度场是通过稳态分析施加的，则须删除稳态分析时施加的温度载荷，选择 Main Menu→Solution→Define Loads→Delete→Thermal→Temperature。
(5)求解。
3.查看结果可以用POST26或POST1查看结果。
33.2问题描述图33-1所示为一个温度为500℃的铁块和一个温度为400℃的铜块，突然放入温度为20℃的完全绝热的水箱中。
忽略水的流动，试分析1h后铜块和铁块的最高温度，以及铜块和铁块的温度变化情况。
材料热物理性能参数如表33-1所示。
图33-1 水箱示意图表33-1 材料热物理性能参数33.3分析步骤拾取菜单Utility Menu→File→Change Jobname，在所弹出对话框的“[/FILNAM]”文本框中输入EXAMPLE33，单击“OK”按钮。
图33-2 改变任务名对话框拾取菜单Main Menu→Preprocessor→Element Type→Add/Edit/Delete，弹出如图33-3所示的对话框，单击“Add…”按钮，弹出如图33-4所示的对话框，在左侧列表中选“Thermal Solid”，在右侧列表中选“8node 77”，单击“OK”按钮，返回到如图33-3所示的对话框，单击“Close”按钮关闭对话框。
图33-3 单元类型对话框图33-4 单元类型库对话框拾取菜单Main Menu→Preprocessor→Material Props→Material Models，弹出如图33-5所示的对话框，在右侧列表中依次拾取“Thermal”、“Conductivity”、“Isotropic”，弹出如图33-6所示的对话框，在“KXX”文本框中输入383（导热系数），单击“OK”按钮；再次拾取如图33-5所示对话框中右侧列表的“Specific Heat”项，弹出如图33-7所示的对话框，在“C”文本框中输入390（比热），单击“OK”按钮；再次拾取如图33-5所示对话框中右侧列表的“Density”项，弹出如图33-8所示的对话框，在“DENS”文本框中输入8889（密度），单击“OK”按钮，于是定义好了材料模型1（铜）。
单击如图33-5所示对话框的菜单项Material--New Model，单击弹出的“Define Material TD”对话框中的“OK”按钮，然后重复定义材料模型1时的各步骤，定义材料模型2（铁）的导热系数为70，比热为448，密度为7833。
重复定义材料模型2时的各步骤，定义材料模型3（水）的导热系数为2，比热为4185，密度为996。
最后关闭如图33-5所示的对话框。
略图33-5 材料模型对话框图33-6 定义导热系数对话框图33-7 定义比热对话框拾取菜单 Main Menu→Preprocessor→ Modeling→Create→ Areas→Rectangle→By Dimension，弹出如图33-9所示的对话框，在“X1，X2”文本框中分别输入0，0.6，在“Y1，Y2”文本框中分别输入0，0.5，单击“Apply”按钮，再次弹出如图33-9所示的对话框，在“Xl, X2”文本框中分别输入0.15, 0.225，在“Y1, Y2”文本框中分别输入 (0.225，0.27)，单击“Apply”按钮，再次弹出如图33-9所示的对话框，在“Xl，X2”文本框中分别输入0.342，0.42)，在“Y1，Y2”文本框中分别输入0.225，0.27，单击“OK”按钮。
拾取菜单Utility Menu→PlotCtrls→Numbering，在弹出的"Plot Numbering Controls "对话框中，将Area numbers（面号）打开，单击“OK”按钮。
图33-8 定义密度对话框图33-9 创建矩形对话框拾取菜单Main Menu→Preprocessor→Modeling→Operate→Booleans→Overlap→Areas,弹出拾取窗口，单击“Pick All”按钮。
拾取菜单Main Menu→Preprocessor→Meshing→MeshTool，弹出如图33-10所示的对话框。
选择“Element Attributes”的下拉列表框为“Areas”，单击下拉列表框后面的“Set”按钮，弹出拾取窗口，选择面2，单击拾取窗口中的“OK”按钮，弹出.”Areas Attributes"对话框，选择“MAT”下拉列表框为1，单击“Apply”按钮，再次弹出拾取窗口，选择面3，单击拾取窗口中的“OK”按钮，选择“Areas Attributes”对话框的“MAT”下拉列表框为2，单击“Apply”按钮，再次窗口，择面4，单击拾取窗口中的“OK”按钮，选择“Areas Attributes”对话框的“MAT”下拉列表框为3，单击"OK"按钮。
单击“Size Controls”区域中“Global”后面的“Set”按钮，弹出如图33-11所示的对话框，在“SIZE”文本中输入0.01，单击“OK”按钮；在如图33-10所示对话框的"Mesh"区域，选择单元形状为“Quad”（四边形），选择划分单元的方法为“Mapped”（映射），单击“Mesh”按钮，弹出拾取窗口，拾取面2和面3，单击“OK”按钮。
图33-10 划分单元工具对话框图33-11 单元尺寸对话框单击“Size Controls”区域中“Global”后面的“Set”按钮，弹出如图33-11所示的对话框，在" SIZE"文本框中输入0.03，单击“OK”按钮；在如图33-10所示对话框的“Mesh”区域，选择单元形状为"Quad "（四边形），选择划分单元的方法为“Free”（自由），单击“Mesh”按钮，弹出拾取窗口，拾取面4，单击“OK”按钮，关闭如图33-10所示的对话框。
拾取菜单Main Menu→Solution→Analysis Type→New Analysis，在弹出的对话框中选择“Type of Analysis”为“Transient”，单击“OK”按钮，在随后弹出的“Transient Analysis”对话框中，单击“OK”按钮。
以下步骤进行稳态分析，以得到初始温度场。
拾取菜单Main Menu→Solution→Load Step Opts→Time/Frequenc→Time Integration→Amplitude Decay，弹出如图33-12所示的对话框，将“TIMINT”关闭，单击“OK”按钮。
图33-12 时间积分控制对话框．拾取菜单Main Menu→Solution→Load Step Opts→Time/Frequenc→Time→Time Step,弹出如图33-13所示的对话框，在“TIME”文本框中输入0.01，在“DELTIM Time step size”文本框中输入0.01，单击“OK”按钮。