下载文档原格式
第N章水杯热分析案例这个案例是使用ANSYS WORKBENCH 热分析模块功能进行演示。
通过对杯子模型加载温度荷载来分析出其温度分布状况。
1.5.1案例介绍此案了使用一个铜合金材料的水杯模型，在内表面施加100℃的温度载荷在外表面施加对流传热系数来模拟一下当水杯装满热水时候的温度以及热流分布状况以演示AMSYS WORKBENCH 热分析模块的基本操作过程。
1.5.2启动Workbench并建立分析项目(1) 首先打开ANSYS WORKBENCH 14.0。
双击Toolbox（工具箱）→Analysis System （分析系统）→Steady-State Therma（热分析模块）打开，如图-1所示。
图-1 打开热分析模块图-2 选择材料1.5.3定义材料参数（1）双击打开A2项目（Engineering Data）。
（2）打开后在菜单栏中点击图标。
然后在Engineering Data Sources项目中单击A3项General Materials（一般材料）→单击Outline of General Materials中的A7项Copper Alloy（铜合金）中的B7项的图标来选择上此材料。
（3）材料选择完成后单击（回到项目管理区）图标。
1.5.4导入模型（1）双击项目A3项Geometry（模型）进入DM模块来导入我们需要分析的模型。
如图-3所示。
图-3 打开DM模块图-5 导入模型（2）单击File（文件）→Import External Geometry File（导入模型文件）。
如图-5所示。
（3）找到模型文件单击选择后单击“打开”按钮。
如图-6所示。
图-6 打开模型图-7 划分网格1.5.5划分网格（1）双击项目文件A4项Model（网格）。
如图-7所示。
打开的模型如图-8所示。
图-8 杯子的模型图-9 刷新网格（2）我们先使用程序自动划分网格，通过查看网格质量再决定是否进行网格控制。
(完整)ANSYS热分析详解编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整)ANSYS热分析详解）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为(完整)ANSYS热分析详解的全部内容。
第一章简介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量〕等。
热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。
二、ANSYS的热分析•在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。
•ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。
•ANSYS热分析包括热传导、热对流及热辐射三种热传递方式.此外，还可以分析相变、有内热源、接触热阻等问题。
三、ANSYS 热分析分类•稳态传热:系统的温度场不随时间变化•瞬态传热：系统的温度场随时间明显变化四、耦合分析•热－结构耦合•热－流体耦合•热－电耦合•热－磁耦合•热－电－磁－结构耦合等第二章 基础知识一、符号与单位 W/m 2—℃ 二、传热学经典理论回顾热分析遵循热力学第一定律,即能量守恒定律:● 对于一个封闭的系统（没有质量的流入或流出〕PE KE U W Q ∆+∆+∆=-式中： Q —— 热量；W -- 作功；∆U ——系统内能；∆KE ——系统动能;∆PE —-系统势能；●对于大多数工程传热问题：0＝＝PE KE ∆∆； ●通常考虑没有做功:0=W ， 则：U Q ∆=； ● 对于稳态热分析：0=∆=U Q ，即流入系统的热量等于流出的热量;●对于瞬态热分析：dt dU q =,即流入或流出的热传递速率q 等于系统内能的变化。
ANSYS热分析指南(第三章)第三章稳态热分析3.1稳态传热的定义ANSYS/Multiphysics，ANSYS/Mechanical，ANSYS/FLOTRAN和ANSYS/Professional这些产品支持稳态热分析。
稳态传热用于分析稳定的热载荷对系统或部件的影响。
通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。
也可以在所有瞬态效应消失后，将稳态热分析作为瞬态热分析的最后一步进行分析。
稳态热分析可以计算确定由于不随时间变化的热载荷引起的温度、热梯度、热流率、热流密度等参数。
这些热载荷包括：对流辐射热流率热流密度（单位面积热流）热生成率（单位体积热流）固定温度的边界条件稳态热分析可用于材料属性固定不变的线性问题和材料性质随温度变化的非线性问题。
事实上，大多数材料的热性能都随温度变化，因此在通常情况下，热分析都是非线性的。
当然，如果在分析中考虑辐射，则分析也是非线性的。
3.2热分析的单元ANSYS和ANSYS/Professional中大约有40种单元有助于进行稳态分析。
有关单元的详细描述请参考《ANSYS Element Reference》，该手册以单元编号来讲述单元，第一个单元是LINK1。
单元名采用大写，所有的单元都可用于稳态和瞬态热分析。
其中SOLID70单元还具有补偿在恒定速度场下由于传质导致的热流的功能。
这些热分析单元如下：表3-1二维实体单元表3-2三维实体单元表3-3辐射连接单元表3-4传导杆单元表3-5对流连接单元表3-6壳单元表3-7耦合场单元表3-8特殊单元3.3热分析的基本过程ANSYS热分析包含如下三个主要步骤：前处理：建模求解：施加荷载并求解后处理：查看结果以下的内容将讲述如何执行上面的步骤。
首先，对每一步的任务进行总体的介绍，然后通过一个管接处的稳态热分析的实例来引导读者如何按照GUI路径逐步完成一个稳态热分析。
最后，本章提供了该实例等效的命令流文件。
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。
/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 element MP,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 tank WPROTA,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 tank ASEL,A,LOC,Y,0 ! Select area at remote Y edge of tank CM,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 thickness LESIZE,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 meshing VMESH,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-chosen Newton-Raphson optionTUNIF,450 ! 给结点统一的温度；Uniform starting temperature at all nodesCSYS,1 ! 1 —Cylindrical with Z as the axis of rotationNSEL,S,LOC,X,RI1 ! Nodes on inner tank surfaceSF,ALL,CONV,250/144,450 ! 为结点指定表面载荷；对流；Convection（对流）；load at all nodesCMSEL,,AREMOTE ! 选择子集组合；Select AREMOTE componentNSLA,,1 ! Nodes belonging to AREMOTED,ALL,TEMP,450 ! 设定边界温度条件Temperature constraints at those nodesWPROTA,0,-90 ! Rotate working plane to pipe axis CSWPLA,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 pipe ESLN ! 选择单元NSLE ! 选择结点/SHOW,,,1 ! 向量显示；Vector mode/TITLE,Thermal flux vectors at pipe/tank junctionPLVECT,TF ! Plot thermal flux（热通量）vectors FINISH。
第10 章热分析典型工程实例本章要点拉伸特征旋转特征扫掠特征混合特征孔特征壳特征本章案例某型号手机电池的散热分析冷库复合隔热板热量流动分析电子元器件散热装置温度分析10.1 工程实例1——某型号手机电池的散热分析该算例为某型手机电池的散热分析，如图10-1为某型号手机背面的照片，图中可见手机的电池的位置。
在手机工作时，电池可向外传递热量。
使用手机的读者应该都体会过手机电池发热的现象，特别是在长时间接打电话时，这种现象尤为明显。
本实例对某型号手机进行分析，电池的标准电压为3.7V，电池容量为750mAh。
试求手机开机状态下外壳的温度分布。
手机的各部分材料性能参数如表10.1所示。
图10-1 手机背面照片在计算分析过程中我们将手机看做三个组成部分：塑料外壳、手机内部材料和手机电池。
忽略手机内部线路和芯片，可以将手机电池看做唯一热源。
简化后的手机模型如图10-2所示，图中单位均为cm。
本实例拟采用Solid Tet 10node 87单元进行分析。
由于电池功率和环境温度均可视为恒定不变，因此分析类型为稳态。
图10-2 简化后的手机模型由电池的电压和电流可以算得电池的功率：==⨯=P UI 3.70.75 2.775W电池的体积为：3=⨯⨯=V0.040.010.050.00002m电池的发热量：3==Q P/V138750W/m——附带光盘“Ch10\实例10－1_start”——附带光盘“Ch10\实例10－1_end”——附带光盘“A VI\Ch10\10－1.avi”1、定义分析文件名1、选择Utility Menu>File>Change Jobname，在弹出的单元增添对话框中输入Example10-1，然后点击OK按钮。
2、选择Main Menu>Preferences，弹出Preferences for GUI Filtering对话框，点选Thermal复选框，单击OK按钮关闭该对话框。
实例1：某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。
几何参数：筒外径30 feet总壁厚 2 inch不锈钢层壁厚0.75 i nch玻纤层壁厚 1 inch铝层壁厚0.25 i nch筒长200 feet导热系数不锈钢8.27 B TU/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度进行分析,如图示。
以下分别列出log文件和菜单文件。
/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,,,2eshape,2 !设定为映射网格划分mat,1amesh,1mat,2amesh,2mat,3amesh,3/SOLUSFL,11,CONV,HAIR,,TAIR !施加空气对流边界SFL,1,CONV,HSEA,,TSEA !施加海水对流边界SOLVE/POST1PLNSOL ！输出温度彩色云图finish菜单操作:1.U tility Menu>File>change jobename, 输入Steady1；2.U tility Menu>File>change title,输入Steady-state thermal analysis of submarine；3.在命令行输入:/units, BFT；4.M ain Menu: Preprocessor；5.M ain Menu: Preprocessor>Element Type>Add/Edit/Delete,选择PLANE55；6.M ain Menu: Preprocessor>Material Prop>-Constant-Isotropic,默认材料编号为1,在KXX框中输入8.27,选择APPL Y,输入材料编号为2,在KXX框中输入0.028,选择APPL Y,输入材料编号为3,在KXX框中输入117.4；7.M ain Menu: Preprocessor>-Modeling->Create>-Areas-Circle>By Dimensions ,在RAD1中输入15,在RAD2中输入15-(.75/12),在THERA1中输入-0.5,在THERA2中输入0.5,选择APPL Y,在RAD1中输入15-(.75/12),在RAD2中输入15-(1.75/12),选择APPL Y,在RAD1中输入15-(1.75/12),在RAD2中输入15-2/12,选择OK；8.M ain Menu: Preprocessor>-Modeling->Operate>-Booleane->Glue>Area,选择PICK ALL；9.M ain Menu: Preprocessor>-Meshing-Size Contrls>-Lines-Picked Lines,选择不锈钢层短边,在NDIV框中输入4,选择APPL Y,选择玻璃纤维层的短边,在NDIV框中输入5,选择APPL Y,选择铝层的短边,在NDIV框中输入2,选择APPL Y,选择四个长边,在NDIV中输入16；10.Main Menu: Preprocessor>-Attributes-Define>Picked Area,选择不锈钢层,在MAT框中输入1,选择APPL Y,选择玻璃纤维层,在MA T框中输入2,选择APPL Y,选择铝层,在MAT框中输入3,选择OK；11.Main Menu: Preprocessor>-Meshing-Mesh>-Areas-Mapped>3 or 4 sided,选择PICK ALL；12.Main Menu: Solution>-Loads-Apply>-Thermal-Convection>On lines,选择不锈钢外壁,在V ALI框中输入80,在V AL2I框中输入44.5,选择APPL Y,选择铝层内壁,在V ALI框中输入2.5,在V AL2I框中输入70,选择OK；13.Main Menu: Solution>-Solve-Current LS；14.Main Menu: General Postproc>Plot Results>-Contour Plot-Nodal Solu,选择Temperature。
对流换热系数定义：流体与固体表面之间的换热能力，比如说，物体表面与附近空气温差1℃，单位时间单位面积上通过对流与附近空气交换的热量。
单位为W/(m^2·℃)。
表面对流换热系数的数值与换热过程中流体的物理性质、换热表面的形状、部位、表面与流体之间的温差以及流体的流速等都有密切关系。
物体表面附近的流体的流速愈大，其表面对流换热系数也愈大。
如人处在风速较大的环境中，由于皮肤表面的对流换热系数较大，其散热（或吸热）量也较大。
对流换热系数可用经验公式计算，通常用巴兹公式计算。
的大致量级：空气自然对流 5 ～25 ，气体强制对流 20 ～100。
实例2:一钢铸件及其砂模的横截面尺寸如图所示:铸钢的热物理性能如下表所示:初始条件:铸钢的温度为2875o F，砂模的温度为80o F；砂模外边界的对流边界条件：对流系数0.014Btu/hr.in2.o F，空气温度80o F；求3个小时后铸钢及砂模的温度分布。
/Title, Casting Solidification!进入前处理/prep7et,1,plane55 ！定义单元mp,dens,1,0.254 ！定义砂模热性能mp,kxx,1,0.025mp,c,1,0.28mptemp,1,0,2643,2750,2875 ！定义铸钢的热性能mpdata,kxx,2,1.44,1.54,1.22,1.22mpdata,enth,2,0,128.1,163.8,174.2mpplot,kxx,2mpplot,enth,2save！创建几何模型k,1,0,0,0k,2,22,0,0k,3,10,12,0k,4,0,12,0/pnum,kp,1/pnum,line,1/pnum,area,1/Triad,ltopkplota,1,2,3,4saverectng,4,22,4,8aplotaovlap,alladele,3aplotsave！划分网格esize,1amesh,5mat,2aplotamesh,4eplot/pnum,elem/number,1save！进入加载求解/SOLUantype,trans ！设定为瞬态分析esel,s,mat,,2 ！设定铸钢的初始温度nsle,s/replotic,all,temp,2875esel,inve ！设定砂模的初始温度nsle,s/replotic,all,temp,80allselsavelplotsfl,1,CONV,0.014,,80 ！设定砂模外边界对流sfl,3,CONV,0.014,,80sfl,4,CONV,0.014,,80/psf,conv,2time,3 ！设定瞬态分析时间kbc,1 ！设定为阶越的载荷autots,on ！打开自动时间步长deltim,0.01,0.001,0.25 ！设定时间步长timint,on ！打开时间积分tintp,,,,1 ！将THETA设定为1 outres,all,all ！输入每个子步的结果solve！进入后处理/post26/pnum,node,1/number,0eplotnsol,2,204,temp,center ！设定铸钢中心点温度随时间的变量plvar,2 ！绘制温度～时间曲线savefinish菜单操作:1.Utility Menu>File>Change Title, 输入Casting Solidification;2.定义单元类型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete, Add, Quad 4node55；3.定义砂模热性能:Main Menu>Preprocessor>Material Props>Isotropic, 默认材料编号1, 在Density(DENS)框中输入0.054, 在Thermal conductivity (KXX)框中输入0.025, 在Specific heat(C)框中输入0.28;4.定义铸钢热性能温度表:Main Menu>Preprocessor>Material Props>-Temp Dependent->TempTable, 输入T1=0,T2=2643, T3=2750, T4=2875;5.定义铸钢热性能:Main Menu>Preprocessor>Material Props>-Temp Dependent ->Prop Table,选择Th Conductivity,选择KXX, 输入材料编号2,输入C1=1.44, C2=1.54, C3=1.22, C4=1.22,选择Apply, 选择Enthalpy,输入C1=0, C2=128.1, C3=163.8, C4=174.2;6.创建关键点:Main Menu>Preprocessor>-Modeling->Create>Keypoints>In Active CS,输入关键点编号1,输入坐标0,0,0, 输入关键点编号2, 输入坐标22,0,0, 输入关键点编号3, 输入坐标10,12,0, 输入关键点编号4, 输入坐标0,12,0;7.创建几何模型:Main Menu>Preprocessor>-Modeling->Create>-Areas->Arbitrary>ThroughKPs,顺序选取关键点1,2,3,4;8.Main Menu>Preprocessor>-Modeling->Create>-Areas->Rectangle>By Dimension,输入X1=4,X2=22,Y1=4,Y2=8;9.进行布尔操作:MainMenu>Preprocessor>-Modeling->Operate>-Booleans->Overlap>Area,Pick all；10.删除多余面:Main Menu>Preprocessor>-Modeling->Delete>Area and Below,311.保存数据库:在Ansys Toolbar中选取SA VE_DB;12.定义单元大小:Main Menu>Preprocessor>-Meshing->Size Cntrls>-Global->Size, 在Element edge length框中输入1;13.对砂模划分网格:Main Menu>Preprocessor>-Meshing->Mesh>-Areas->Free,选择砂模;14.对铸钢划分网格:Main Menu>Preprocessor>-Attributes->Define>Default Attribs, 在Material number菜单中选择2;15.Main Menu>Preprocessor>-Meshing->Mesh>-Areas->Free,选择铸钢;16.定义分析类型:Main Menu>Solution>-Analysis Type->New Analysis, 选择Transient;17.选择铸钢上的节点:Utility Menu>Select>Entities, 选择element,mat,输入2,选择Apply,选择node, attached to element,选择OK;18.定义铸钢的初始温度:Main Menu>Solution>-Loads->Apply>Initial Condit’n>Define, 选择Pick all,选择temp, 输入2875, OK;19.选择砂模上的节点:Utility Menu>Select>Entities,Nodes, inverse20.定义砂模的初始温度:Main Menu>Solution>-Loads->Apply>Initial Condit’n>Define, 选择Pick all, 选择temp, 输入80, OK;21.Utility Menu>Select>Everything；22.Utility Menu>Plot>Lines；23.定义对流边界条件:Main Menu>Solution>-Loads->Apply>-Thermal->Converction>OnLines,选择砂模的三个边界1,3,4, 在file coefficent框中输入80, 在Bulk temperature框中输入, 80;24.设定瞬态分析时间选项:Main Menu>Solution>Load Step Opts>Time/Frequenc>Time-Time Step,Time at end of load step 3Time Step size 0.01Stepped or ramped b.c. SteppedAutomatic time stepping onMinimun time Step size 0.001Maximum time step size 0.2525.设置输出:Main Menu>Solution>Load Step Opts>Output Ctrls>DB/Results File, 在Filewrite frequency框中选择Every substep;26.求解:Main Menu>Solution>-Solve->Current LS；27.进入后处理: Main Menu>Timehist Postproc；28.定义铸钢中心节点的温度变量:Main Menu>Timehist Postproc>Define Variables, Add,Nodal DOF result,2,204；29.绘制节点温度随时间变化曲线：Main Menu>Timehist Postproc>Graph Variable,2。
