ANSYS热分析.doc

No BoundariesANSYS热分析指南第一章简 介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量等。热分析在许多工程应用中扮演重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等。二、ANSYS的热分析 在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。 ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。 ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题。三、ANSYS 热分析分类 稳态传热：系统的温度场不随时间变化 瞬态传热：系统的温度场随时间明显变化四、耦合分析 热结构耦合 热流体耦合 热电耦合 热磁耦合 热电磁结构耦合等第二章基础知识一、符号与单位项目国际单位英制单位ANSYS代号长度mft时间ss质量Kglbm温度oF力Nlbf能量（热量）JBTU功率（热流率）WBTU/sec热流密度W/m2BTU/sec-ft2生热速率W/m3BTU/sec-ft3导热系数W/m-BTU/sec-ft-oFKXX对流系数W/m2-BTU/sec-ft2-oFHF密度Kg/m3lbm/ft3DENS比热J/Kg-BTU/lbm-oFC焓J/m3BTU/ft3ENTH二、传热学经典理论回顾热分析遵循热力学第一定律，即能量守恒定律:l 对于一个封闭的系统（没有质量的流入或流出 式中:Q 热量;W 作功;系统内能;系统动能；系统势能；l 对于大多数工程传热问题：；l 通常考虑没有做功：, 则：；l 对于稳态热分析：，即流入系统的热量等于流出的热量；l 对于瞬态热分析：，即流入或流出的热传递速率q等于系统内能的变化。三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：，式中为热流密度（W/m2），为导热系数（W/m-），“-”表示热量流向温度降低的方向。2、热对流热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述：，式中h为对流换热系数（或称膜传热系数、给热系数、膜系数等），为固体表面的温度，为周围流体的温度。3、热辐射热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温度越高，单位时间辐射的热量越多。热传导和热对流都需要有传热介质，而热辐射无须任何介质。实质上，在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬波尔兹曼方程来计算：，式中为热流率，为辐射率（黑度），为斯蒂芬波尔兹曼常数，约为5.6710-8W/m2.K4，A1为辐射面1的面积，为由辐射面1到辐射面2的形状系数，为辐射面1的绝对温度，为辐射面2的绝对温度。由上式可以看出，包含热辐射的热分析是高度非线性的。四、稳态传热如果系统的净热流率为，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为（以矩阵形式表示）式中：为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；为节点温度向量；为节点热流率向量，包含热生成；ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件，生成、以及。五、瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理，瞬态热平衡可以表达为（以矩阵形式表示）：式中:为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；为比热矩阵,考虑系统内能的增加;为节点温度向量；为温度对时间的导数;为节点热流率向量，包含热生成。六、线性与非线性如果有下列情况产生，则为非线性热分析：、材料热性能随温度变化，如K(T),C(T)等；、边界条件随温度变化，如h(T)等；、含有非线性单元；4、 考虑辐射传热非线性热分析的热平衡矩阵方程为：七、边界条件、初始条件ANSYS热分析的边界条件或初始条件可分为七种：温度、热流率、热流密度、对流、辐射、绝热、生热。八、热分析误差估计 仅用于评估由于网格密度不够带来的误差； 仅适用于SOLID或SHELL的热单元（只有温度一个自由度）； 基于单元边界的热流密度的不连续； 仅对一种材料、线性、稳态热分析有效； 使用自适应网格划分可以对误差进行控制。第三章稳态传热分析一、稳态传热的定义稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数二、热分析的单元热分析涉及到的单元有大约40种，其中纯粹用于热分析的有14种：线性：LINK32两维二节点热传导单元LINK33三维二节点热传导单元LINK34二节点热对流单元LINK31二节点热辐射单元二维实体：PLANE55四节点四边形单元PLANE77八节点四边形单元PLANE35三节点三角形单元PLANE75四节点轴对称单元PLANE78八节点轴对称单元三维实体SOLID87六节点四面体单元SOLID70八节点六面体单元SOLID90二十节点六面体单元壳SHELL57四节点点MASS71有关单元的详细解释，请参阅ANSYS Element Reference Guide三、ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤： 前处理：建模 求解：施加载荷计算 后处理：查看结果1、建模、确定jobname、title、unit；、进入PREP7前处理，定义单元类型，设定单元选项；、定义单元实常数；、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化；、创建几何模型并划分网格，请参阅ANSYS Modeling and Meshing Guide。2、施加载荷计算1、 定义分析类型l 如果进行新的热分析：Command: ANTYPE, STATIC, NEWGUI: Main menuSolution-Analysis Type-New AnalysisSteady-statel 如果继续上一次分析，比如增加边界条件等：Command: ANTYPE, STATIC, RESTGUI: Main menuSolutionAnalysis Type-Restart2、 施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件) ：a、恒定的温度通常作为自由度约束施加于温度已知的边界上。Command Family:DGUI：Main MenuSolution-Loads-Apply-Thermal-Temperatureb、热流率热流率作为节点集中载荷，主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷)，如果输入的值为正，代表热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时，尤其要注意。此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。Command Family: FGUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flowc、对流对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可施加于实体和壳模型上，对于线模型，可以通过对流线单元LINK34考虑对流。Command Family: SFGUI：Main MenuSolution-Loads-Apply-Thermal-Convectiond、热流密度热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRAN CFD计算得到时，可以在模型相应的外表面施加热流密度。如果输入的值为正，代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面，但ANSYS仅读取最后施加的面载进行计算。Command Family: FGUI：Main MenuSolution-Loads-Apply-Thermal-Heat Fluxe、生热率生热率作为体载施加于单元上，可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。Command Family: BFGUI：Main MenuSolution-Loads-Apply-Thermal-Heat Generat、确定载荷步选项对于一个热分析，可以确定普通选项、非线性选项以及输出控制。a. 普通选项 时间选项：虽然对于稳态热分析，时间选项并没有实际的物理意义，但它提供了一个方便的设置载荷步和载荷子步的方法。Command: TIMEGUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps 每载荷步中子步的数量或时间步大小：对于非线性分析，每一载荷步需要多个子步。Command: NSUBSTGUI: Main MenuSolution-Load Step Opts-Time/FrequencTime and SubstpsCommand: DELTIMGUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step 递进或阶越选项：如果定义阶越(stepped)选项，载荷值在这个载荷步内保持不变；如果为递进(ramped)选项，则载荷值由上一载荷步值到本载荷步值随每一子步线性变化。Command: KBCGUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substpsb. 非线性选项 迭代次数：本选项设置每一子步允许的最多的迭代次数。默认值为25，对大数热分析问题足够。Command: NEQITGUI: Main MenuSolution-Load Step Opts-NolinearEquilibrium Iter 自动时间步长: 对于非线性问题，可以自动设定子步间载荷的增长，保证求解的稳定性和准确性。Command: AUTOTSGUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps 收敛误差：可根据温度、热流率等检验热分析的收敛性。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: Main MenuSolution-Load Step Opts-NolinearLine Search 预测矫正：本选项可激活每一子步第一次迭代对自由度求解的预测矫正。Command: PREDGUI: Main MenuSolution-Load Step Opts-NolinearPredictorc. 输出控制 控制打印输出：本选项可将任何结果数据输出到*.out 文件中。Command: OUTPRGUI: Main MenuSolution-Load Step Opts-Output CtrlsSolu Printout 控制结果文件：控制*.rth的内容。Command: OUTRES GUI: Main MenuSolution-Load Step Opts-Output CtrlsDB/Results File、确定分析选项a. Newton-Raphson选项（仅对非线性分析有用）Command: NROPT GUI: Main MenuSolutionAnalysis Optionsb. 选择求解器：可选择如下求解器中一个进行求解： Frontal solver(默认） Jacobi Conjugate Gradient(JCG) solver JCG out-of-memory solver Incomplete Cholesky Conjugate Gradient(ICCG) solver Pre-Conditioned Conjugate Gradient Solver(PCG) Iterative(automatic solver selection option)Command: EQSLV GUI: Main MenuSolutionAnalysis Options注意：热分析可选用Iterative选项进行快速求解，但如下情况除外： 热分析包含SURF19或SURF22或超单元； 热辐射分析； 相变分析 需要restart an analysisc. 确定绝对零度：在进行热辐射分析时，要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度，此值应设定为273；如果使用的是华氏度，则为460。Command: TOFFSTGUI: Main MenuSolutionAnalysis Options5、 保存模型: 点击ANSYS工具条SAVE_DB。、求解Command: SOLVEGUI: Main MenuSolutionCurrent LS3、后处理ANSYS将热分析的结果写入*.rth文件中，它包含如下数据：基本数据： 节点温度导出数据： 节点及单元的热流密度 节点及单元的热梯度 单元热流率 节点的反作用热流率 其它对于稳态热分析，可以使用POST1进行后处理，关于后处理的完整描述，可参阅ANSYS Basic Analysis Procedures Guide。进入POST1后，读入载荷步和子步：Command: SET GUI: Main MenuGeneral Postproc-Read Results-By Load Step可以通过如下三种方式查看结果： 彩色云图显示Command: PLNSOL, PLESOL, PLETAB等GUI: Main MenuGeneral PostprocPlot ResultsNodal Solu, Element Solu, Elem Table 矢量图显示Command: PLVECTGUI: Main MenuGeneral PostprocPlot ResultsPre-defined or Userdefined 列表显示Command: PRNSOL, PRESOL, PRRSOL等GUI: Main MenuGeneral PostprocList ResultsNodal Solu, Element Solu, Reaction Solu详细过程请参阅ANSYS Basic Analysis Procedures Guide。实例1：某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。几何参数：筒外径30feet总壁厚2inch不锈钢层壁厚0.75inch玻纤层壁厚1inch铝层壁厚0.25inch筒长200feet导热系数不锈钢8.27BTU/hr.ft.oF玻纤0.028BTU/hr.ft.oF铝117.4BTU/hr.ft.oF边界条件空气温度70oF海水温度44.5oF空气对流系数2.5BTU/hr.ft2.oF海水对流系数80BTU/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=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. Utility MenuFilechange jobename, 输入Steady1；2. Utility MenuFilechange title,输入Steady-state thermal analysis of submarine；3. 在命令行输入:/units, BFT；4. Main Menu: Preprocessor；5. Main Menu: PreprocessorElement TypeAdd/Edit/Delete,选择PLANE55；6. Main Menu: PreprocessorMaterial 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/12),选择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,选择四个长边,在NDIV中输入16；10. Main Menu: Preprocessor-Attributes-DefinePicked 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-ConvectionOn 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。实例2一圆筒形的罐有一接管,罐外径为3英尺，壁厚为0.2英尺，接管外径为0.5英尺，壁厚为0.1英尺，罐与接管的轴线垂直且接管远离罐的端部。如图所示：罐内流体温度为华氏450度,与罐壁的对流换热系数年为250BUT/hr-ft2-oF,接管内流体的温度为华氏100度,与管壁的对流换热系数随管壁温度而变。接管与罐为同一种材料，它的热物理性能如下表所示：温度70200300400500oF密度0.2850.2850.2850.2850.285lbm/in3导热系数8.358.909.359.810.23Btu/hr-ft-oF比热0.1130.1170.1190.1220.125Btu/lbm-oF对流系数*426405352275221Btu/hr-ft2-oF*接管内壁对流系数求罐与接管的温度分布。以下分别列出LOG文件及菜单操作/prep7/title,Steady-state thermal analysis of pipe junction/units,bin！使用英制单位et,1,90！定义热单元mp,dens,1,.285！密度mptemp,70,200,300,400,500！建立温度表mpdata,kxx,1,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12！导热系数mpdata,c,1,0.133,0.177,0.119,0.122,0.125！比热mpdata,hf,2,426/144,405/144,352/144,275/144,221/144！接管对流系数！定义几何模型参数ri1=1.3！罐内半径ro1=1.5！罐外半径z1=2！罐长ri2=0.4！接管内半径ro2=0.5！接管外半径z2=2！接管长！建立几何模型cylind,ri1,ro1,z1,90！1/4罐体wprota,0,-90！将工作平面旋转到垂直于接管轴线cylind,ri2,ro2,z2,-90！1/4接管wpstyl,defa！将工作平面恢复到默认状态vovlap,1,2！进行OVERLAP布尔操作/pnum,volu,1！打开实体编号/view,-3,-1,1！定义显示角度/type,4/title, Volumes used in building pipe/tank junctionvplot！显示实体vdele,3,4,1！删除多余实体！划分网格asel,loc,z,z1！选择罐上Z=Z1的面 asel,a,loc,y,0！添加选择罐上Y=0的面cm,aremote,area！创建名为AREMOTE的面组/pnum,area,1/pnum,line,1/title,lines showing the portion being modeledaplot/noeraselplot/eraseaccat,all！组合罐远端的面及线，为映射划分网！格作准备lccat,12,7lccat,10,5lesize,20,4！在接管壁厚方向分4等分lesize,40,6！在接管长度方向分6等分lesize,6,4！在罐壁厚方向分4等分allsel！选择EVERYTHINGesize,0.4！设定默认的单元大小mshape,0,3d！选择3D映射网格mshkey,1save！保存数据文件vmesh,all！划分网格，产生节点与单元/pnum,defa/title, elements in portion being modeledeplot！显示单元finish！加载求解/soluantype,static！定义为稳态分析nropt,auto！设置求解选项为Program-chosen ！Newton-Raphsontunif,450！设定初始所有节点温度csys,1！变为柱坐标nsel,s,loc,x,ri1！选择罐内表面的节点sf,all,conv,250/144,450！定义对流边界条件cmsel,aremote！选择AREMOTE面组nsla,1！选择属于AREMOTE面组的节点d,all,temp,450！定义节点温度wprota,0,-90！将工作平面旋转到垂直于接管轴线cswpla,11,1！创建局部柱坐标nsel,s,loc,x,ri2！选择接管内壁的节点sf,all,conv,-2,100！定义对流边界条件allsel！选择EVERYTHING/pbc,temp,1！显示所有温度约束/psf,conv,2！显示所有对流边界/title,Boundary conditionsnplot！显示节点wpstyle,defa！工作平面恢复默认状态csys,0！变为直角坐标autots,on！打开自动步厂长nsubst,50！设定子步数量kbc,0！设定为阶越outpr,nsol,last！设置输出solve！进行求解finish！进入后处理/post1/title,Temperature contrours at pipe/tank junctionplnsol,temp！显示温度彩色云图finish/exit,all菜单操作1、 设定标题:Utility MenuFileChange Title,输入Steady-State analysis of pipe junction,选择OK;2、 设定单位制:在命令提示行输入/UNITS,BIN;3、 定义单元类型:Main MenuPreprocesorElement TypeAdd/Edit/Delete,选择Thermal Solid, Bricck 20 node 90号单元;4、 定义材料属性(1) Main MenuPreprocessorMaterial Props-Constant-Isotropic,默认材料编号1,在DENSITY框中输入0.285;(2) Main MenuPreprocessorMaterial Props-Temp Dependent-Temp Table,输入温度70,200,300,400,500;(3) Main MenuPreprocessorMaterial Props-Temp Dependent-Prop Table,选择导热系数KXX,材料编号为1,输入与温度表对应的导热系数8.35/12,8.9/12,9.35/12,9.8/12,10.23/12,选择APPLY;(4) 选择比热C,材料编号为1,输入0.113,0.117,0.119,0.122,0.125,选择APPLY;(5) 选择对流系数HF,材料编号为2,输入426/144,405/144,352/144,275/144, 221/144，选择OK。5、 定义几何模型参数：Utility MenuParametersScalar Parameters,输入ri1=1.3,ro1=1.5,z1=2,ri2=0.4,ro2=0.5,z2=2;6、 建立几何模型(1) Main MenuPreprocessor-Modeling-Create-Volumes-CylinderBy Dimensions, Outer radius框中输入ro1,Optional inner radium框中输入ri1,Z coordinates框中输入0和Z1,Ending angle框中输入90;(2) Utility MenuWorkPlaneOffset WP by Increments,在XY,YZ,ZX框中输入0,-90;(3) Main MenuPreprocessor-Modeling-Create-Volumes-CylinderBy Dimensions; Outer radius框中输入ro2, Optional inner radium框中输入ri2, Z coordinates框中输入0和Z2,Starting angle框中输入-90,Ending angle框中输入0;(4) Utility MenuWorkPlaneAlign WP withGlobal Cartesian；7、 进行布尔操作:Main MenuPreprocessor-Modeling-Operate-Booleans- Overlap Volumes,选择Pick All；8、 观察几何模型(1) Utility MenuPlotCtrlsNumbering,打开volumes;(2) Utility MenuPlotCtrlsView Direction, 在Coords of view point框中输入-3,-1,1；9、 删除多余实体Main MenuPreprocessor-Modeling-DeleteVolume and Below,在命令输入行输入3,4回车；10、 创建组AREMOTE(1) Utility MenuSelectEntities,选择Area, By location, Z Coordinates, 在Min, Max框中输入Z1,选择APPLY,Y Coordinates, 在Min, Max框中输入0,OK;(2) Utility MenuSelectComp/AssemblyCreate Component,在Component name框中输入AREMOTE, 在Components is made of菜单中选择AREA；11、 组合面及线(1) Main MenuPreprocessor-Meshing-Mesh-Volumes-Mapped-Concatenate-Area,选择Pick all；(2) Main MenuPreprocessor-Meshing-Mesh-Volumes-Mapped-Concatenate-Lines,在命令行中输入12,7回车,选择APPLY,在命令行中输入10,5回车,OK；12、 设定网格密度(1) Main MenuPreprocessor-Meshing-Size CntrlsPicked Lines,选择线6和20,OK,在No. of element divisions框中输入4,OK；(2) Main MenuPreprocessor-Meshing-Size CntrlsPicked Lines,选择线40,OK,在No. of element divisions框中输入6,OK；(3) Utility MenuSelectEverything；(4) Main MenuPreprocessor-Meshing-Size Cntrls-Global-Size,在element edge length框中输入0.4,OK；13、 划分网格:Main MenuPreprocessor-Meshing-Mesh-Volumes-Mapped4 to 6 sides,选择Pick All；14、 定义求解类型及选项(1) Main MenuSolution-Analysis Type-New Analysis,选择Steady-State；(2) Main MenuSolution-Analysis Options,选择Program-chosen；15、 施加对流载荷(1) Utility MenuWorkPlaneChange Active CS toGlobal Cylindrical；(2) Utility MenuSelectEntities,选择Nodes, By location, X,在Min, Max框中输入ri1,OK；(3) Main MenuSolution-Loads-Apply-Thermal-ConvectionOn Nodes,选择Pick All, 输入250/144及450,OK；16、 在AREMOTE组上施加温度约束(1) Utility MenuSelectComp/AssemblySelect Comp/Assembly,选aremote;(2) Utility MenuSelectEntities,选择Nodes, Attached to, On the Area all, OK；(3) Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes,选择Pick all,输入45,OK；17、 施加与温度有关的对流边界条件(1) Utility MenuWorkPlaneOffset WP by Increments,在XY,YZ,ZX Angles框中输入0,-90,OK;(2) Utility MenuWorkPlaneLocal Coordinate SystemsCreate Local CSAt WP Origin,在Type of coordinate system菜单中,选择Cylindrical 1,OK；(3) Utility MenuSelect Entities,选择Nodes, By location, X, 在Min, Max框中输入ri2,OK；(4) Main MenuSolution-Loads-Apply-Thermal-ConvectionOn Nodes,选择Pick All,在Film coefficient框中输入-2,在Bulk temperature框中输入100,OK;(5) Utility MenuSelectEverything；(6) Utility MenuPlotCtrlsSymbols,在Show pres and convect as菜单中选择Arrow, OK;(7) Utility MenuPlotNodes；18、 恢复工作平面及坐标系统(1) Utility MenuWorkPlaneChange Active CS toGlobal Cartesian；(2) Utility MenuWorkPlaneAlign WP withGlobal Cartesian；19、 设定载荷步选项: Main MenuSolution-Load Step Options-Time/FrequencTime and Substeps,在Number of substeps框中输入50,设置Automatic time stepping为On；20、 求解:Main MenuSolution-Solve-Current LS21、 显示温度分布彩色云图: Main MenuGeneral PostprocPlot Results-Contour Plot-Nodal Solu,选择Temperature TEMP。ANSYS Verification Manual中关于稳态热分析的实例：VM58Centerline temperature of a heat 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 finVM99Temperature distribution in a trapezoidal finVM100Heat conductivity across a chimney sectionVM101Temperature distribution in a short solid cylinderVM102Cylinder with temperature dependent conductivityVM103Thin plate with a central heat sourceVM105Heat generation 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 profileVM193Adaptive analysis of two-dimensional heat transfer with convection第四章瞬态传热分析一、瞬态传热分析的定义瞬态热分析用于计算一个系统的随时间变化的温度场及其它热参数。在工程上一般用瞬态热分析计算温度场，并将之作为热载荷进行应力分析。瞬态热分析的基本步骤与稳态热分析类似。主要的区别是瞬态热分析中的载荷是随时间变化的。为了表达随时间变化的载荷，首先必须将载荷时间曲线分为载荷步。载荷时间曲线中的每一个拐点为一个载荷步，如下图所示。对于每一个载荷步，必须定义载荷值及时间值，同时必须选择载荷步为渐变或阶越。二、瞬态热分析中的单元及命令瞬态热分析中使用的单元与稳态热分析相同。要了解每个单元的详细说明，请参阅ANSYS Element Reference Guide。要了解每个命令的详细功能，请参阅ANSYS Command Reference Guide。三、ANSYS 瞬态热分析的主要步骤 建模 加载求解 后处理四、建模 确定jobname、title、units, 进入PREP7； 定义单元类型并设置选项； 如果需要，定义单元实常数； 定义材料热性能：一般瞬态热分析要定义导热系数、密度及比热； 建立几何模型； 对几何模型划分网格。关于建模及划分网格，请参阅ANSYS Modeling and Meshing Guide。五、加载求解1、定义分析类型 如果第一次进行分析，或重新进行分析GUI: Main MenuSolutionAnalysis TypeNew AnalysisTransientCommand: ANTYPE,TRANSIENT,NEW 如果接着上次的分析继续进行（例如增加其它载荷）GUI: Main MenuSolutionAnalysis TypeRestartCommand: ANTYPE,TRANSIENT,REST2、获得瞬态热分析的初始条件、定义均匀温度场如果已知模型的起始温度是均匀的，可设定所有节点初始温度Command:TUNIFGUI: Main Menu Solution-Loads-SettingsUniform Temp如果不在对话框中输入数据，则默认为参考温度，参考温度的值默认为零，但可通过如下方法设定参考温度：Command:TREFGUI: Main Menu Solution-Loads-SettingsReference Temp注意：设定均匀的初始温度，与如下的设定节点的温度（自由度）不同Command:DGUI: Main MenuSolution-Loads-Apply-Thermal-TemperatureOn Nodes初始均匀温度仅对分析的第一个子步有效；而设定节点温度将保持贯穿整个瞬态分析过程，除非通过下列方法删除此约束：Command:DDELEGUI: Main Menu Solution-Loads-Delete-Thermal-TemperatureOn Nodes、设定非均匀的初始温度在瞬态热分析中，节点温度可以设定为不同的值：Command:ICGUI: Main Menu SolutionLoadsApply-Initial ConditnDefine如果初始温度场是不均匀的且又是未知的，就必须首先作稳态热分析确定初始条件： 设定载荷（如已知的温度、热对流等） 将时间积分设置为OFF：Command: TIMINT, OFFGUI: Main Menu PreprocessorLoads-Load Step Opts-Time/FrequencTime Integration 设定一个只有一个子步的，时间很小的载荷步（例如0.001）：Command: TIMEGUI: Main Menu PreprocessorLoads-Load S