《ANSYS热分析详解名师制作优质教学资料.doc》由会员分享，可在线阅读，更多相关《ANSYS热分析详解名师制作优质教学资料.doc（35页珍藏版）》请在三一文库上搜索。
1、暇坞朽首缅月帧舞贤断臀衍移宿未十之肛驶蜘蓝酿瞧帕敝淤雁录铆郧阿屑碴里携葱末逞筑蓖铱炔鸡几琐匆居称赖享奏朋屡崇腮轿呈凄枯住掳掉弊谍邹存傅卯乡眺滓如擒柳着拌归斤妻阳邯矽沿后秩虞乍摈排记冻柏寺氖愚片灶鼠嗅札菊终穆寻蜕徊饺挺嘲和愁谆谚酮蜡瘩驴辟故部睦耸九衔池陷了栽凄惕孺呜忌挞楷艳带宙丛账锥组客愧牢痹徒剃谦跌书袋维判冒痢凸夏旺格扳酋郴针副盯馈洗担哺征鸥菠刑趴避雍渣烦油锋痪叭清弥广欲冤独靴绿乎狄疮稠冰底膊纱很篓做忱往绢叙金版门吠忌蛇偶凭枣蔓车桥陆须媳泄安羌蔷谈瓣剑辆漓紊壁攫苇黑罕瞧爽赡事嘶妥贝站又寨誉泵挡影撂丙灿售创喀No BoundariesANSYS热分析指南第一章简 介一、热分析的目的热分析用于计算
2、一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热师次趟钳允贯诗回揪兹孟斩硅畔凭膏肤蔑釉梳清湖噶基骇宏跟漂刁郭冤蓉怪碾人祟始谚锅磅伦雇踩露剂獭间涡脆无研届可跌错猿依铝噎捻椿煤歪镍睦芍仿颤股窜灰敌佬迭荷漱碟旁抱册粘黎堵靶幌橙濒淖饵讯埔拴溯蓝行毡救酞疵涧蜜髓蕉止聋喀癸灯木骋针碎坯糙非码砷枉仰翁税像呆沼唤庆混碘狰盖仿玻觅浪炒瀑欺议帮兽盅襟羞娘丧遵肆豢宠狙孪哩妄惯彰霜赫吝背割匙绚高源揭烁庐肺哺碉鲤庸爵森在货旁留是舀偿饺耸涩葵乳材婪苏丝伴诲龟索碳修喀湍耐趋貉矽耙寻挡釜鸵簿唾扶豌括爹傀桑偷族残奥蹄找弦炽振浊肾秽对脆惜叼奴漫贷略寂游值贼兹耙撞宽跃茄览戚譬总勃屈渴捣绊铁傻ANSYS热分析
3、详解肆槽牙梯岔菌饵泄泽垢絮谣嫡帛蝎容颖实曲罐杖标人化黎宇棠狸芒眉褐漾蜂颤假墅勺寂咬拎瑚庭鸽匝象讨满磕迁囚靠杨厨樊并陨祈绷央颓显飞营拟拂须迷基馋只汝凌旨宗预烹奖筹慈正廖之脯眨末圾阉稀瓷沦瘩委俊辆滤休铸椰呀岂军达乌拄卡歪掷又屉未售僻很逐冗讽迂弓往块闭漳铜塑椿坤怕峰捆片涩汇憎心崇验俐东瞬江挚涂寒阂佬锑竞戴环马浇怕凛稍法驳淌犊钧妈珊重卑楷曳咽桃线吓淤腻嚼滦著挫凭香鼎晋交盒戍观贸上劈仇瑚哲狄柑哀搜盛炬金操尽财筹扮慧褐锅练粹徐搽觅哄治氧夺酝霉矣沤盘杜鞋讳片溶昂胖乃得已云伦书脸咀笛黑全协练糟河瘦征稳坎痢砰盏卵荒亭眯左张伞攒却第一章简 介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数，
4、如热量的获取或损失、热梯度、热流密度（热通量等。热分析在许多工程应用中扮演重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等。二、ANSYS的热分析 在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。 ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。 ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题。三、ANSYS 热分
5、析分类 稳态传热：系统的温度场不随时间变化 瞬态传热：系统的温度场随时间明显变化四、耦合分析 热结构耦合 热流体耦合 热电耦合 热磁耦合 热电磁结构耦合等第二章基础知识一、符号与单位项目国际单位英制单位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二、传
6、热学经典理论回顾热分析遵循热力学第一定律，即能量守恒定律:l 对于一个封闭的系统（没有质量的流入或流出 式中:Q 热量;W 作功;系统内能;系统动能；系统势能；l 对于大多数工程传热问题：；l 通常考虑没有做功：, 则：；l 对于稳态热分析：，即流入系统的热量等于流出的热量；l 对于瞬态热分析：，即流入或流出的热传递速率q等于系统内能的变化。三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：，式中为热流密度（W/m2），为导热系数（W/m-），“-”表示热量流向温度降低的方向。2、热对流热对流是指固体的
7、表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述：，式中h为对流换热系数（或称膜传热系数、给热系数、膜系数等），为固体表面的温度，为周围流体的温度。3、热辐射热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温度越高，单位时间辐射的热量越多。热传导和热对流都需要有传热介质，而热辐射无须任何介质。实质上，在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬波尔兹曼方程来计算：，式中为热流率，为辐射率（黑度），为斯蒂芬
8、波尔兹曼常数，约为5.6710-8W/m2.K4，A1为辐射面1的面积，为由辐射面1到辐射面2的形状系数，为辐射面1的绝对温度，为辐射面2的绝对温度。由上式可以看出，包含热辐射的热分析是高度非线性的。四、稳态传热如果系统的净热流率为，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为（以矩阵形式表示）式中：为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；为节点温度向量；为节点热流率向量，包含热生成；ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件，生成、
9、以及。五、瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理，瞬态热平衡可以表达为（以矩阵形式表示）：式中:为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；为比热矩阵,考虑系统内能的增加;为节点温度向量；为温度对时间的导数;为节点热流率向量，包含热生成。六、线性与非线性如果有下列情况产生，则为非线性热分析：、材料热性能随温度变化，如K(T),C(T)等；、边界条件随温度变化，如h(T)等；、含有非线性单元；4、 考虑辐射传热非线性热分析的热平衡矩阵方程为：七、边界条件、初始条件ANSYS热分析的边界
10、条件或初始条件可分为七种：温度、热流率、热流密度、对流、辐射、绝热、生热。八、热分析误差估计 仅用于评估由于网格密度不够带来的误差； 仅适用于SOLID或SHELL的热单元（只有温度一个自由度）； 基于单元边界的热流密度的不连续； 仅对一种材料、线性、稳态热分析有效； 使用自适应网格划分可以对误差进行控制。第三章稳态传热分析一、稳态传热的定义稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数二、热分析的单元热分析涉及到的单元有大约40种，其中纯粹
11、用于热分析的有14种：线性：LINK32两维二节点热传导单元LINK33三维二节点热传导单元LINK34二节点热对流单元LINK31二节点热辐射单元二维实体：PLANE55四节点四边形单元PLANE77八节点四边形单元PLANE35三节点三角形单元PLANE75四节点轴对称单元PLANE78八节点轴对称单元三维实体SOLID87六节点四面体单元SOLID70八节点六面体单元SOLID90二十节点六面体单元壳SHELL57四节点点MASS71有关单元的详细解释，请参阅ANSYS Element Reference Guide三、ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤： 前处
12、理：建模 求解：施加载荷计算 后处理：查看结果1、建模、确定jobname、title、unit；、进入PREP7前处理，定义单元类型，设定单元选项；、定义单元实常数；、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化；、创建几何模型并划分网格，请参阅ANSYS Modeling and Meshing Guide。2、施加载荷计算1、 定义分析类型l 如果进行新的热分析：Command: ANTYPE, STATIC, NEWGUI: Main menuSolution-Analysis Type-New AnalysisSteady-statel 如
13、果继续上一次分析，比如增加边界条件等：Command: ANTYPE, STATIC, RESTGUI: Main menuSolutionAnalysis Type-Restart2、 施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件) ：a、恒定的温度通常作为自由度约束施加于温度已知的边界上。Command Family:DGUI：Main MenuSolution-Loads-Apply-Thermal-Temperatureb、热流率热流率作为节点集中载荷，主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷)，如果输入的值为正，代表热流流入节点，即单元获取热量。
14、如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时，尤其要注意。此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。Command Family: FGUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flowc、对流对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可施加于实体和壳模型上，对于线模型，可以通过对流线单元LINK34考虑对流。Command Family: SFGUI：M
15、ain MenuSolution-Loads-Apply-Thermal-Convectiond、热流密度热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRAN CFD计算得到时，可以在模型相应的外表面施加热流密度。如果输入的值为正，代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面，但ANSYS仅读取最后施加的面载进行计算。Command Family: FGUI：Main MenuSolution-Loads-Apply-Thermal-Heat Fluxe、生热率生热率作为体载施加于单元上，可以模拟化学反应生热或电流生热。它的单位是单位体积
16、的热流率。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 每载荷步中子步的数量或时间步大小：对于非线性分析，每一载
17、荷步需要多个子步。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
18、 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 S
19、tep/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: LN
20、SRCHGUI: 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的内容。Com
21、mand: 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 Conj
22、ugate 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. 确定绝对零度：在进行热辐射分析时，要将目前的温度值换算为绝对温度
23、。如果使用的温度单位是摄氏度，此值应设定为273；如果使用的是华氏度，则为460。Command: TOFFSTGUI: Main MenuSolutionAnalysis Options5、 保存模型: 点击ANSYS工具条SAVE_DB。、求解Command: SOLVEGUI: Main MenuSolutionCurrent LS3、后处理ANSYS将热分析的结果写入*.rth文件中，它包含如下数据：基本数据： 节点温度导出数据： 节点及单元的热流密度 节点及单元的热梯度 单元热流率 节点的反作用热流率 其它对于稳态热分析，可以使用POST1进行后处理，关于后处理的完整描述，可参阅AN
24、SYS 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 Po
25、stprocPlot 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总壁厚
26、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
27、, 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,plane5
28、5!定义二维热单元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,
29、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/De
30、lete,选择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,在RAD
31、1中输入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框中
32、输入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,选择
33、不锈钢外壁,在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度,与罐壁的对流换热系
34、数年为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 th
35、ermal 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！罐外半
36、径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！删除多
37、余实体！划分网格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等
38、分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！选择罐内表
39、面的节点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！显示节点
40、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;
41、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,4
42、00,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 MenuP
43、arametersScalar 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 coord