热分析(ansys教程).ppt

《热分析(ansys教程).ppt》由会员分享，可在线阅读，更多相关《热分析(ansys教程).ppt（99页珍藏版）》请在三一办公上搜索。

1、热分析,热分析的目的 热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度(热通量)等热分析在许多工程应用中扮演重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等,ANSYS的热分析,在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题,

2、ANSYS的热分析分类,ANSYS的热分析分类 稳态传热：系统的温度场不随时间变化瞬态传热：系统的温度场随时间明显变化与热有关的耦合分析 热结构耦合热流体耦合热电耦合热磁耦合热电磁结构耦合等,热分析的符号与单位,表征物体吸收的热量，为一个体系的内能与体系的体积和外界施加于体系的压强的乘积之和,传热学经典理论回顾,热传递的方式,1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：qn=-k*(dT/dx)，式中qn为热流密度（W/m2），k为导热系数(W/m-)，“-”表示热量流向温度降低的方向。2、热对流热对流是指固体的

3、表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述：qn=h*(TS-TB)，式中h为对流换热系数(或称膜传热系数、给热系数、膜系数等)，TS为固体表面的温度，TB为周围流体的温度。,热传递的方式(续),3、热辐射热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温度越高，单位时间辐射的热量越多。热传导和热对流都需要有传热介质，而热辐射无须任何介质。实质上，在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬波尔兹曼

4、方程来计算：q=A1F12(T14-T24)，式中q为热流率，为辐射率（黑度），为斯蒂芬波尔兹曼常数，约为5.6710-8W/m2.K4，A1为辐射面1的面积，F12为由辐射面1到辐射面2的形状系数，T1为辐射面1的绝对温度，T2为辐射面2的绝对温度。由上式可以看出，包含热辐射的热分析是高度非线性的。,稳态传热,如果系统的净热流率为，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为(以矩阵形式表示):KT=Q式中：K为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；T

5、为节点温度向量；Q为节点热流率向量，包含热生成；ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件，生成K、T以及Q。,瞬态传热,瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理，瞬态热平衡可以表达为(以矩阵形式表示)：C+KT=Q式中:K为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；C为比热矩阵,考虑系统内能的增加;T为节点温度向量；为温度对时间的导数;Q为节点热流率向量，包含热生成。,线性与非线性,如果有下列情况产生，则为非线性热分析：材料热性能随温度变化，如K(T),C(T)等；边界条

6、件随温度变化，如h(T)等；含有非线性单元；考虑辐射传热非线性热分析的热平衡矩阵方程为：C(T)+K(T)T=Q(T),边界条件、初始条件,ANSYS热分析的边界条件或初始条件可分为七种：温度：模型区温度已知热流率：热流率已知的点对流：表面的热传递给周围的流体通过对流。输入对流换热系数h和环境流体的 平均温度Tb热辐射：通过辐射产生热传递的面.输入辐射系数，Stefan-Boltzmann常数，“空间节点”的温度作为可选项输入绝热面：“完全绝热”面，该面上不发生热传递热通量：单位面积上的热流率已知的面热生成率：体的生热率已知的区域,热分析误差估计,仅用于评估由于网格密度不够带来的误差；仅适用于

7、SOLID或SHELL的热单元(只有温度一个自由度)；基于单元边界的热流密度的不连续；仅对一种材料、线性、稳态热分析有效；使用自适应网格划分可以对误差进行控制。,稳态传热分析,稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数,热分析的单元,热分析涉及到的单元有大约40种，其中纯粹用于热分析的有14种：线性：LINK32两维二节点热传导单元 LINK33三维二节点热传导单元 LINK34二节点热对流单元 LINK31二节点热辐射单元二维实体：PL

8、ANE55四节点四边形单元 PLANE77八节点四边形单元 PLANE35三节点三角形单元 PLANE75四节点轴对称单元 PLANE78八节点轴对称单元,热分析的单元(续),三维实体：SOLID87六节点四面体单元 SOLID70八节点六面体单元 SOLID90二十节点六面体单元壳：SHELL57四节点点：MASS71,稳态热分析的基本过程,ANSYS热分析可分为三个步骤：前处理：建模求解：施加载荷计算后处理：查看结果,建模,确定jobname、title、unit；进入PREP7前处理，定义单元类型，设定单元选项；定义单元实常数；定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可

9、以是恒定的，也可以随温度变化；创建几何模型并划分网格。,几何尺寸(模型),既可用ANSYS建立模型，也可用其它方法建好模型后导入模型建好后，以上两种建模方法的具体过程将不再显示,划分网格,首先定义单元属性:单元类型,实常数,材料属性.单元类型下表给出了常用的热单元类型每个结点只有一个自由度:温度,划分网格(续),材料属性必须输入导热系数,KXX如果施加了内部热生成率，则需指定比热(C)ANSYS提供的材料库(/ansys57/matlib)包括几种常用材料的结构属性 和热属性,但是建议用户创建、使用自己的材料库把优先设置为“热分析”，使材料模型图形用户界面只显示材料的热属性实常数主要应用于壳单

10、元和线单元,划分网格(续),划分网格存储数据文件使用 MeshTool 划分网格，使用缺省的智能网格划分级别6可以生成很好的初始网格至此完成前处理，下面开始求解,施加载荷计算,定义分析类型如果进行新的热分析：Command:ANTYPE,STATIC,NEWGUI:Main menuSolution-Analysis Type-New AnalysisSteady-state如果继续上一次分析，比如增加边界条件等：Command:ANTYPE,STATIC,RESTGUI:Main menuSolutionAnalysis Type-Restart,施加载荷计算(续),施加载荷可以直接在实体模

11、型或单元模型上施加五种载荷(边界条件)：a、恒定的温度通常作为自由度约束施加于温度已知的边界上。Command Family:DGUI：Main MenuSolution-Loads-Apply-Thermal-Temperature,施加载荷计算(续),b、热流率热流率作为节点集中载荷，主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷)，如果输入的值为正，代表热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率

12、时，尤其要注意。此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。Command Family:FGUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flow,施加载荷计算(续),c、对流对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可施加于实体和壳模型上，对于线模型，可以通过对流线单元LINK34考虑对流。Command Family:SFGUI：Main MenuSolution-Loads-Apply-Thermal-Convection,施加载荷计算(续),d、热流密度热流密度也是一种面载。当通过单位面积的热流率已

13、知或通过FLOTRAN CFD计算得到时，可以在模型相应的外表面施加热流密度。如果输入的值为正，代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面，但ANSYS仅读取最后施加的面载进行计算。Command Family:FGUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flux,施加载荷计算(续),e、生热率生热率作为体载施加于单元上，可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。Command Family:BFGUI：Main MenuSolution-Loads-Apply-Thermal-He

14、at Generat,施加载荷计算(续),定载荷步选项对于一个热分析，可以确定普通选项、非线性选项以及输出控制。a.普通选项时间选项：虽然对于稳态热分析，时间选项并没有实际的物理意义，但它提供了一个方便的设置载荷步和载荷子步的方法。Command:TIMEGUI:Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps,施加载荷计算(续),每载荷步中子步的数量或时间步大小：对于非线性分析，每一载荷步需要多个子步。Command:NSUBSTGUI:Main MenuSolution-Load St

15、ep Opts-Time/Frequenc Time and SubstpsCommand:DELTIMGUI:Main MenuSolution-Load Step Opts-Time/Frequenc Time-Time Step,施加载荷计算(续),递进或阶越选项：如果定义阶越(stepped)选项，载荷值在这个载荷步内保持不变；如果为递进(ramped)选项，则载荷值由上一载荷步值到本载荷步值随每一子步线性变化。Command:KBCGUI:Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Sub

16、stps,施加载荷计算(续),b.非线性选项迭代次数：本选项设置每一子步允许的最多的迭代次数。默认值为25，对大数热分析问题足够。Command:NEQITGUI:Main MenuSolution-Load Step Opts-NolinearEquilibrium Iter,施加载荷计算(续),自动时间步长:对于非线性问题，可以自动设定子步间载荷的增长，保证求解的稳定性和准确性。Command:AUTOTSGUI:Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps,施加载荷计算(续),收

17、敛误差：可根据温度、热流率等检验热分析的收敛性。Command:CNVTOLGUI:Main MenuSolution-Load Step Opts-NolinearConvergence Crit,施加载荷计算(续),求解结束选项：如果在规定的迭代次数内，达不到收敛，ANSYS可以停止求解或到下一载荷步继续求解。Command:NCNVGUI:Main MenuSolution-Load Step Opts-NolinearCriteria to Stop,施加载荷计算(续),线性搜索：设置本选项可使ANSYS用Newton-Raphson方法进行线性搜索。Command:LNSRCHGUI

18、:Main MenuSolution-Load Step Opts-NolinearLine Search,施加载荷计算(续),预测矫正：本选项可激活每一子步第一次迭代对自由度求解的预测矫正。Command:PREDGUI:Main MenuSolution-Load Step Opts-NolinearPredictor,施加载荷计算(续),c.输出控制控制打印输出：本选项可将任何结果数据输出到*.out 文件中。Command:OUTPRGUI:Main MenuSolution-Load Step Opts-Output CtrlsSolu Printout,施加载荷计算(续),控制结果

19、文件：控制*.rth的内容。Command:OUTRES GUI:Main MenuSolution-Load Step Opts-Output CtrlsDB/Results File,施加载荷计算(续),确定分析选项a.Newton-Raphson选项(仅对非线性分析有用)Command:NROPT GUI:Main MenuSolutionAnalysis Options,施加载荷计算(续),b.选择求解器：可选择如下求解器中一个进行求解：Frontal solver(默认）Jacobi Conjugate Gradient(JCG)solverJCG out-of-memory sol

20、ver Incomplete Cholesky Conjugate Gradient(ICCG)solverPre-Conditioned Conjugate Gradient Solver(PCG)Iterative(automatic solver selection option),施加载荷计算(续),Command:EQSLV GUI:Main MenuSolutionAnalysis Options注意：热分析可选用Iterative选项进行快速求解，但如下情况除外：热分析包含SURF19或SURF22或超单元；热辐射分析；相变分析 需要restart an analysis,施加载

21、荷计算(续),c.确定绝对零度：在进行热辐射分析时，要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度，此值应设定为273；如果使用的是华氏度，则为460。Command:TOFFSTGUI:Main MenuSolutionAnalysis Options,施加载荷计算(续),保存模型:点击ANSYS工具条SAVE_DB。求解Command:SOLVEGUI:Main MenuSolutionCurrent LS,后处理,ANSYS将热分析的结果写入*.rth文件中，它包含如下数据：基本数据：节点温度导出数据：节点及单元的热流密度节点及单元的热梯度单元热流率节点的反作用热流率其它,后

22、处理(续),对于稳态热分析，可以使用POST1进行后处理，进入POST1后，读入载荷步和子步：Command:SET GUI:Main MenuGeneral Postproc-Read Results-By Load Step,后处理(续),可以通过如下三种方式查看结果：彩色云图显示Command:PLNSOL,PLESOL,PLETAB等GUI:Main Menu General Postproc Plot ResultsNodal Solu,Element Solu,Elem Table,后处理(续),对3-D 实体模型绘制云图时，选项isosurfaces(等值面)是非常有用的.用/C

23、TYPE 命令或 Utility Menu PlotCtrls Style Contours Contour Style.,后处理(续),矢量图显示Command:PLVECTGUI:Main Menu General Postproc Plot ResultsPre-defined or Userdefined,后处理(续),列表显示Command:PRNSOL,PRESOL,PRRSOL等GUI:Main Menu General Postproc List ResultsNodal Solu,Element Solu,Reaction Solu,后处理(续),检查结果是否正确温度是否在预

24、期的范围内?在指定温度和热流边界的基础上，估计预期的范围网格大小是否满足精度?和受力分析一样，可以画出非均匀分布的温度梯度(单元解)并找出高梯度的单元.这些区域可作为重新定义网格时的参考若节点温度梯度（平均的）和单元温度梯度（非平均的）之间的差别很大，则可能是网格划分太粗糙,练习-稳态传热分析,某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。几何参数：筒外径30 feet 总壁厚2 inch 不锈钢层壁厚0.75 inch 玻纤层壁厚1 inch 铝层壁厚0.25 inch 筒长200 feet导热

25、系数：不锈钢8.27BTU/hr.ft.oF 玻纤0.028BTU/hr.ft.oF 铝117.4BTU/hr.ft.oF,边界条件,边界条件：空气温度 70oF 海水温度 44.5oF空气对流系数 2.5 BTU/hr.ft2.oF海水对流系数 80 BTU/hr.ft2.oF沿垂直于圆筒轴线作横截面，得到一圆环，取其中1度进行分析，如图示。,log文件,/filename,Steady1/title,Steady-state thermal analysis of submarine/units,BFTRo=15!外径(ft)Rss=15-(0.75/12)!不锈钢层内径ft)Rins=1

26、5-(1.75/12)!玻璃纤维层内径(ft)Ral=15-(2/12)!铝层内径(ft)Tair=70!潜水艇内空气温度Tsea=44.5!海水温度Kss=8.27!不锈钢的导热系数(BTU/)Kins=0.028!玻璃纤维的导热系数(BTU/),log文件(续),Kal=117.4!铝的导热系数(BTU/)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,

27、Kal!设定铝的导热系数pcirc,Ro,Rss,-0.5,0.5!创建几何模型pcirc,Rss,Rins,-0.5,0.5pcirc,Rins,Ral,-0.5,0.5aglue,all,log文件(续),numcmp,arealesize,1,16!设定划分网格密度lesize,4,4lesize,14,5lesize,16,2eshape,2!设定为映射网格划分mat,1amesh,1mat,2amesh,2mat,3amesh,3,log文件(续),/SOLUSFL,11,CONV,HAIR,TAIR!施加空气对流边界SFL,1,CONV,HSEA,TSEA!施加海水对流边界SOLV

28、E/POST1PLNSOL！输出温度彩色云图finish,菜单操作,菜单操作:1.Utility MenuFilechange jobname,输入Steady1；2.Utility MenuFilechange title,输入Steady-state thermal analysis of submarine；3.在命令行输入:/units,BFT；4.Main Menu:Preprocessor；,菜单操作(续),5.Main Menu:Preprocessor Element Type Add/Edit/Delete,选择PLANE55；6.Main Menu:Preprocessor

29、 Material Prop-Constant-Isotropic,默认材料编号为1,在KXX框中输入8.27,选择APPLY,输入材料编号为2,在KXX框中输入0.028,选择APPLY,输入材料编号为3,在KXX框中输入117.4；,菜单操作(续),7.Main Menu:Preprocessor-Modeling-Create-Areas-CircleBy Dimensions 在RAD1中输入15,在RAD2中输入15-(.75/12),在THERA1中输入-0.5,在THERA2中输入0.5,选择APPLY;在RAD1中输入15-(.75/12),在RAD2中输入15-(1.75/1

30、2),选择APPLY;在RAD1中输入15-(1.75/12),在RAD2中输入15-2/12,选择OK；,菜单操作(续),8.Main Menu:Preprocessor-Modeling-Operate-Booleane-GlueArea,选择PICK ALL；,菜单操作(续),9.Main Menu:Preprocessor-Meshing-Size Contrls-Lines-Picked Lines选择不锈钢层短边,在NDIV框中输入4,选择APPLY;,菜单操作(续),选择玻璃纤维层的短边,在NDIV框中输入5,选择APPLY;选择铝层的短边,在NDIV框中输入2,选择APPLY;

31、选择四个长边,在NDIV中输入16；,菜单操作(续),10.Main Menu:Preprocessor-Attributes-Define Picked Area选择不锈钢层,在MAT框中输入1,选择APPLY;选择玻璃纤维层,在MAT框中输入2,选择APPLY;选择铝层,在MAT框中输入3,选择OK；,菜单操作(续),11.Main Menu:Preprocessor-Meshing-Mesh-Areas-Mapped3 or 4 sided,选择PICK ALL；,菜单操作(续),12.Main Menu:Solution-Loads-Apply-Thermal-Convection O

32、n lines选择不锈钢外壁,在VALI框中输入80,在VAL2I框中输入44.5,选择APPLY;,菜单操作(续),选择铝层内壁,在VALI框中输入2.5,在VAL2I框中输入70,选择OK；,菜单操作(续),13.Main Menu:Solution-Solve-Current LS；14.Main Menu:General PostprocPlot Results-Contour Plot-Nodal Solu,选择Temperature。,ANSYS Verification Manual中关于稳态热分析的实例,VM58Centerline temperature of a heat

33、generating wireVM92Insulted wall temperatureVM93Temperature dependent conductivityVM94Heat generating plateVM95Heat transfer from a cooling spineVM96Temperature distribution in a short solid cylinderVM97Temperature distribution along a straight finVM98Temperature distribution along a tapered finVM99

34、Temperature distribution in a trapezoidal finVM100Heat conductivity across a chimney section,ANSYS Verification Manual中关于稳态热分析的实例(续),VM101Temperature distribution in a short solid cylinderVM102Cylinder with temperature dependent conductivityVM103Thin plate with a central heat sourceVM105Heat generat

35、ion coil with temperature dependent conductivityVM108Temperature gradient across a solid cylinderVM118Centerline temperature of a heat generating wireVM160Solid cylinder with harmonic temperature loadVM161Heat flow from a insulated pipeVM162Cooling of a circular fin of rectangular profileVM193Adapti

36、ve analysis of two-dimensional heat transfer with convection,练习-瞬态传热分析,一个30公斤重、温度为70的铜块,以及一个20公斤重、温度为80的铁块，突然放入温度为20、盛满了300升水的、完全绝热的水箱中,如图所示。过了一个小时，求铜块与铁块的最高温度(假设忽略水的流动)。材料热物理性能如下:,练习-瞬态传热分析(续),log文件,/filename,transient1/title,Thermal Transient Exercise 1!进入前处理/prep7et,1,plane77!定义单元类型mp,kxx,1,383!定

37、义材料热性能参数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,log文件(续),rectnag,0,0.6,0,0.5!创建几何实体rectang,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,1aplotaatt,1,1,1!划分网格eshape,2esize,0.02amesh,2aat

38、t,2,1,1amesh,3aatt,3,1,1eshape,3esize,0.05amesh,4,log文件(续),/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,80,log文件(续),esel,s,mat,1nsle,sd,all,temp,70allselsolve!得到初始温度分布time,3600!进行

39、瞬态分析timint,on!打开时间积分deltim,26,2,200!设置时间步长，最大及最小时间步长autots,on!打开自动时间步长ddelet,all,temp!删除稳态分析中定义的节点温度outres,all,1!将每个子步的值写入数据库文件solvefinishsave,log文件(续),!进入POST26后处理/post26solu,2,dtime,dtime!2每一子步采用的时间步长nsol,3,node(0.1875,0.2475,0),temp,T_Copper!3铜块的中心点nsol,4,node(0.371,0.247,0),temp,T_Iron!4铁块的中心点ns

40、ol,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,2plvar,3,4,5,6,7,8finish,log文件(续),!进入POST1后处理/post1!设置为最后一个载荷子步set,lastesel,s,mat,1nsle,splnsol,tempesel,s,mat,2nsle,splnsol,te

41、mpfinish,菜单操作,1、Utility Menu:FileChange Jobname,输入文件名Transient1；2、Utility Menu:FileChange Title,输入Thermal Transient Exercise 1；,19.16.2 菜单操作(续),3、Main MenuPreprocessor，进入前处理；4、Main MenuPreprocessorElement Type Add/Edit/Delete,选择单元plane77；,菜单操作(续),5、Main MenuPreprocessorMaterial Props Isotropic,定义材料1

42、（铜）的KXX等于 383、DENS等于8889、C等于1390；,菜单操作(续),6、Main MenuPreprocessorMaterial Props Isotropic,定义材料2（铁）的KXX等于 70、DENS等于7837、C等于448；7、Main MenuPreprocessorMaterial Props Isotropic,定义材料3（水）的KXX等于0.61、DENS等于996、C等于4185；,菜单操作(续),8、Main MenuPreprocessor-Modeling-Create-Areas-Retangle By Dimensions输入X1=0,Y1=0,

43、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.,菜单操作(续),9、Main MenuPreprocessor-Modeling-Operate Booleans Overlap,选择Pick All；10、Utility Menu:PlotctrlsNumberingAreas,on；11、Utility Menu:PlotAreas；,菜单操作(续),12、Main Me

44、nuPreprocessor-Attributes-Define-All Areas,选择材料1；13、Main MenuPreprocessorMeshing-Size Cntrls-Manualsize-Global-Size,输入单元大小0.02；14、Main MenuPreprocessorMeshing-Mesh-Areas-Mapped 3 or 4 sided,选择铜块；,菜单操作(续),15、Main MenuPreprocessor-Attributes-Define-All Areas,选择材料2；16、Main MenuPreprocessorMeshing-Mesh-

45、Areas-Mapped 3 or 4 sided,选择铁块；17、Main MenuPreprocessor-Attributes-Define-All Areas,选择材料3；18、Main MenuPreprocessorMeshing-Size Cntrls-Manualsize-Global-Size,输入单元大小0.05；19、Main MenuPreprocessorMeshing-Mesh-Areas-Free,选择水箱；,菜单操作(续),菜单操作(续),20、Utility MenuPlotArea；21、Main MenuSolution，进入加载求解；22、Main Me

46、nuSolution-Analysis Type-New Analysis,选择Transient，定义为瞬态分析；,菜单操作(续),23、Main MenuSolution-Load Step Opts Time/FrenquencTime Integration,TIMINT,off，首先进行稳态分析；,菜单操作(续),24、Main MenuSolution-Load Step Opts Time/FrenquencTime-Time Step，设定TIME为0.01、DELTIM也为0.01；,菜单操作(续),25、Utility Menu:SelectElementmat,输入3,选

47、择Apply,选择NodesAttached to,选择Element；,菜单操作(续),26、Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes,选择Pick All,输入20；,菜单操作(续),27、Utility Menu:SelectElementmat,输入2,选择Apply,选择NodesAttached to,选择Element；28、Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes,选择Pick All,输入80；29、Utility Menu:Se

48、lectElementmat,输入1,选择Apply,选择NodesAttached to,选择Element；30、Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes,选择Pick All,输入70；31、Utility Menu:Select Everything,菜单操作(续),菜单操作(续),32、Main MenuSolution-Solve-Current LS33、Main MenuSolution-Load Step Opts Time/FrenquencTime-Time Step，设定TIME=3600,DEL

49、TIM=26,最小、最大时间步长分别为2,200,将Autots设置为ON；,菜单操作(续),34、Main MenuSolution-Load Step Opts Time/FrenquencTime Integration,将TIMINT设置为ON；,菜单操作(续),35、Main MenuSolution-Loads-Delete-Thermal-TemperatureOn Nodes,选择Pick All，删除稳态分析定义的节点温度；36、Main MenuSolution-Load Step OptsOutput Ctrls-DB/Results,选择Every Substeps；,

50、菜单操作(续),37、Main MenuSolution-Solve-Current LS；38、Main MenuTimeHist PostPro，进入POST26；39、Main MenuTimeHist PostProDefine Variables AddSolution summary，在User specified label框中输入dtime，选择Solution ItemsStep Time，选择OK，AddNodal DOF result，在Node number框中输入node(0.1875,0.2475,0),在User specified label框中输入T_Copp