2022年热分析 .pdf

《2022年热分析 .pdf》由会员分享，可在线阅读，更多相关《2022年热分析 .pdf（36页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、No BoundariesANSYS热 分 析 指 南第一章简 介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量等。热分析在许多工程应用中扮演重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等。二、ANSYS的热分析在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、 ANSYS/ED五 种 产 品 中 包 含 热 分 析 功 能 , 其 中ANSYS/FLOTRAN不含相变热分析。ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算

2、各节点的温度，并导出其它热物理参数。ANSYS热分析包括热传导 、 热对流 及 热辐射 三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题。三、ANSYS 热分析分类稳态传热：系统的温度场不随时间变化瞬态传热：系统的温度场随时间明显变化四、耦合分析热结构耦合热流体耦合热电耦合热磁耦合热电磁结构耦合等名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南第二章基础知识一

3、、符号与单位项目国际单位英制单位ANSYS代号长度m ft 时间s s 质量Kg lbm 温度oF 力N lbf 能量（热量）J BTU 功率（热流率）W BTU/sec 热流密度W/m2BTU/sec-ft2生热速率W/m3BTU/sec-ft3导热系数W/m- BTU/sec-ft-oF KXX 对流系数W/m2-BTU/sec-ft2-oF HF 密度Kg/m3lbm/ft3DENS 比热J/Kg- BTU/lbm-oF C 焓J/m3BTU/ft3ENTH 二、传热学经典理论回顾热分析遵循热力学第一定律，即能量守恒定律: 对于一个封闭的系统（没有质量的流入或流出PEKEUWQ式中 :

4、Q 热量 ; W 作功 ; U系统内能; KE系统动能；PE系统势能；对于大多数工程传热问题：0PEKE；通常考虑没有做功：0W, 则：UQ；对于稳态热分析：0UQ，即流入系统的热量等于流出的热量；对于瞬态热分析：dtdUq， 即流入或流出的热传递速率q 等于系统内能的变化。三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：dxdTkq，式中q为热流名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - -

5、第 2 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南密度（ W/m2） ，k为导热系数（W/m- ）， “ -”表示热量流向温度降低的方向。2、热对流热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述：)(BSTThq，式中h 为对流换热系数（或称膜传热系数、给热系数、膜系数等），TS为固体表面的温度，TB为周围流体的温度。3、热辐射热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温度越高，单位时间辐射的热量越多。热传导

6、和热对流都需要有传热介质，而热辐射无须任何介质。实质上，在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并 吸 收 热 量 。 它 们 之 间 的 净 热 量 传 递 可 以 用 斯 蒂 芬 波 尔 兹 曼 方 程 来 计 算 ：qA FTT1121424()，式中q为热流率，为辐射率（黑度），为斯蒂芬波尔兹曼常数，约为5.67 10-8W/m2.K4， A1为辐射面1 的面积，F12为由辐射面1 到辐射面2 的形状系数，T1为辐射面1 的绝对温度，T2为辐射面2 的绝对温度。由上式可以看出，包含热辐射的热分析是高度非线性的。四、稳态传热如果系统的

7、净热流率为，即流入系统的热量加上系统自身产生的热量等于流出系统的热量： q流入+q生成-q流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为（以矩阵形式表示）KTQ式中：K为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；T为节点温度向量；Q为节点热流率向量，包含热生成；ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件，生成K、T以及Q。五、瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、 热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理， 瞬态热平衡可以表达为 （以矩阵形式表示） ：C

8、 TKTQ式中 : K为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；C为比热矩阵,考虑系统内能的增加; T为节点温度向量；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南T为温度对时间的导数; Q为节点热流率向量，包含热生成。六、线性与非线性如果有下列情况产生，则为非线性热分析：、材料热性能随温度变化，如K(T),C(T)等；、边界条件随温度变化，如h(T) 等；、

9、含有非线性单元；、考虑辐射传热非线性热分析的热平衡矩阵方程为：C TTK TTQ T七、边界条件、初始条件ANSYS热分析的边界条件或初始条件可分为七种：温度、热流率、热流密度、对流、辐射、绝热、生热。八、热分析误差估计仅用于评估由于网格密度不够带来的误差；仅适用于SOLID或 SHELL的热单元（只有温度一个自由度）；基于单元边界的热流密度的不连续；仅对一种材料、线性、稳态热分析有效；使用自适应网格划分可以对误差进行控制。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共

10、36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南第三章稳态传热分析一、稳态传热的定义稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数二、热分析的单元热分析涉及到的单元有大约40 种，其中纯粹用于热分析的有14 种：线性：LINK32 两维二节点热传导单元LINK33 三维二节点热传导单元LINK34 二节点热对流单元LINK31 二节点热辐射单元二维实体： PLANE55 四节点四边形单元PL

11、ANE77 八节点四边形单元PLANE35 三节点三角形单元PLANE75 四节点轴对称单元PLANE78 八节点轴对称单元三维实体SOLID87 六节点四面体单元SOLID70 八节点六面体单元SOLID90 二十节点六面体单元壳SHELL57 四节点点MASS71 有关单元的详细解释，请参阅ANSYS Element Reference Guide三、ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤：前处理：建模求解：施加载荷计算后处理：查看结果1、建模、确定jobname 、 title 、 unit ；、进入PREP7 前处理，定义单元类型，设定单元选项；、定义单元实常数；

12、、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化；、创建几何模型并划分网格，请参阅ANSYS Modeling and Meshing Guide 。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南2、施加载荷计算、定义分析类型如果进行新的热分析：Command: ANTYPE, STATIC, NEW GUI: Main menu

13、Solution-Analysis Type-New AnalysisSteady-state 如果继续上一次分析，比如增加边界条件等：Command: ANTYPE, STATIC, REST GUI: Main menuSolutionAnalysis Type-Restart 、施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件) ：a、恒定的温度通常作为自由度约束施加于温度已知的边界上。Command Family:D GUI ： Main MenuSolution-Loads-Apply-Thermal-Temperature b、热流率热流率作为节点集中载荷，主要用于线单

14、元 模型中 (通常线单元模型不能施加对流或热流密度载荷)，如果输入的值为正，代表热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意 ：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时，尤其要注意。此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。Command Family: F GUI ： Main MenuSolution-Loads-Apply-Thermal-Heat Flow c、对流对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可

15、施加于实体和壳模型上，对于线模型，可以通过对流线单元LINK34考虑对流。Command Family: SF GUI ： Main MenuSolution-Loads-Apply-Thermal-Convection d、热流密度热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRAN CFD计算得到时，可以在模型相应的外表面施加热流密度。如果输入的值为正，代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面，但ANSYS仅读取最后施加的面载进行计算。Command Family: F GUI ： Main MenuSolution-Loads-

16、Apply-Thermal-Heat Flux e、生热率生热率作为体载施加于单元上，可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。Command Family: BF GUI ： Main MenuSolution-Loads-Apply-Thermal-Heat Generat 、确定载荷步选项对于一个热分析，可以确定普通选项、非线性选项以及输出控制。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 36 页 - - - - - - - - - No Bo

17、undariesANSYS热 分 析 指 南a. 普通选项时间选项：虽然对于稳态热分析，时间选项并没有实际的物理意义，但它提供了一个方便的设置载荷步和载荷子步的方法。Command: TIME GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps 每载荷步中子步的数量或时间步大小：对于非线性分析，每一载荷步需要多个子步。Command: NSUBST GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime and Substps

18、 Command: DELTIM GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step 递进或阶越选项：如果定义阶越(stepped) 选项， 载荷值在这个载荷步内保持不变；如果为递进(ramped) 选项，则载荷值由上一载荷步值到本载荷步值随每一子步线性变化。Command: KBC GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps b. 非线性选项迭代次数：本选项设置每一子步允许的最多的迭代次数。默

19、认值为25，对大数热分析问题足够。Command: NEQIT GUI: Main MenuSolution-Load Step Opts-NolinearEquilibrium Iter 自动时间步长: 对于非线性问题，可以自动设定子步间载荷的增长，保证求解的稳定性和准确性。Command: AUTOTS GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps 收敛误差：可根据温度、热流率等检验热分析的收敛性。Command: CNVTOL GUI: Main MenuSolutio

20、n-Load Step Opts-NolinearConvergence Crit 求解结束选项：如果在规定的迭代次数内，达不到收敛，ANSYS可以停止求解或到下一载荷步继续求解。Command: NCNV GUI: Main MenuSolution-Load Step Opts-NolinearCriteria to Stop 线性搜索：设置本选项可使ANSYS用 Newton-Raphson方法进行线性搜索。Command: LNSRCH GUI: Main MenuSolution-Load Step Opts-NolinearLine Search 预测矫正：本选项可激活每一子步第一

21、次迭代对自由度求解的预测矫正。Command: PRED GUI: Main MenuSolution-Load Step Opts-NolinearPredictor c. 输出控制控制打印输出：本选项可将任何结果数据输出到*.out 文件中。Command: OUTPR GUI: Main MenuSolution-Load Step Opts-Output CtrlsSolu Printout 控制结果文件：控制*.rth 的内容。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第

22、 7 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南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(J

23、CG) 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 或超单元；热辐射分析

24、；相变分析需要restart an analysis c. 确定绝对零度：在进行热辐射分析时，要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度，此值应设定为273；如果使用的是华氏度，则为460。Command: TOFFST GUI: Main MenuSolutionAnalysis Options 、保存模型: 点击 ANSYS工具条SAVE_DB 。、求解Command: SOLVE GUI: Main MenuSolutionCurrent LS 3、后处理ANSYS将热分析的结果写入*.rth文件中，它包含如下数据：基本数据：节点温度导出数据：节点及单元的热流密度节点及单

25、元的热梯度单元热流率节点的反作用热流率其它对于稳态热分析，可以使用POST1 进行后处理，关于后处理的完整描述，可参阅 ANSYS Basic Analysis Procedures Guide 。进入 POST1 后，读入载荷步和子步：Command: SET GUI: Main MenuGeneral Postproc-Read Results-By Load Step 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 36 页 - - - - - - - - - No

26、 BoundariesANSYS热 分 析 指 南可以通过如下三种方式查看结果：彩色云图显示Command: PLNSOL, PLESOL, PLETAB等GUI: Main MenuGeneral PostprocPlot ResultsNodal Solu, Element Solu, Elem Table矢量图显示Command: PLVECT GUI: Main MenuGeneral PostprocPlot ResultsPre-defined or Userdefined 列表显示Command: PRNSOL, PRESOL, PRRSOL等GUI: Main MenuGene

27、ral PostprocList ResultsNodal Solu, Element Solu, Reaction Solu详细过程请参阅ANSYS Basic Analysis Procedures Guide 。实例1：某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。几何参数：筒外径30 feet 总壁厚2 inch 不锈钢层壁厚0.75 inch 玻纤层壁厚1 inch 铝层壁厚0.25 inch 筒长200 feet 导热系数不锈钢8.27 B TU/hr.ft.oF 玻纤0.028 BT

28、U/hr.ft.oF 铝117.4 BTU/hr.ft.oF 边界条件空气温度70 oF 海水温度44.5 oF 空气对流系数2.5 BTU/hr.ft2.oF 海水对流系数80 BTU/hr.ft2.oF 沿垂直于圆筒轴线作横截面，得到一圆环，取其中1 度进行分析,如图示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南不锈钢玻璃纤维铝R15 feet海水空气以下分别列出

29、log 文件和菜单文件。/filename, Steady1 /title, Steady-state thermal analysis of submarine /units, BFT Ro=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

30、.4 !铝的导热系数 (BTU/hr.ft.oF) Hair=2.5 !空气的对流系数(BTU/hr.ft2.oF) Hsea=80 !海水的对流系数(BTU/hr.ft2.oF) /prep7 et,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.5 pcirc,Rins,Ral,-0.5,0.5 aglue,all numcmp,area lesize,1,1

31、6 !设定划分网格密度lesize,4,4 lesize,14,5 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南lesize,16,2 eshape,2 !设定为映射网格划分mat,1 amesh,1 mat,2 amesh,2 mat,3 amesh,3 /SOLU SFL,11,CONV ,HAIR,TAIR !施加空气对流边界SFL,1,CONV,HSEA,TS

32、EA !施加海水对流边界SOLVE /POST1 PLNSOL ！输出温度彩色云图finish 菜单操作: 1.Utility MenuFilechange jobename, 输入 Steady1 ；2.Utility MenuFilechange title,输入 Steady-state thermal analysis of submarine；3.在 命令行输入:/units, BFT ；4.M ain Menu: Preprocessor；5.M ain Menu: PreprocessorElement TypeAdd/Edit/Delete, 选择 PLANE55 ；6.M a

33、in Menu: PreprocessorMaterial Prop-Constant-Isotropic, 默认材料编号为1,在 KXX 框中输入 8.27,选择 APPLY,输入材料编号为2,在 KXX 框中输入0.028,选择 APPLY,输入材料编号为 3,在 KXX 框中输入 117.4；7.M ain Menu: Preprocessor-Modeling-Create-Areas-CircleBy Dimensions ,在 RAD1 中输入 15,在 RAD2 中输入15-(.75/12),在 THERA1中输入 -0.5,在 THERA2中输入0.5,选择APPLY,在 RA

34、D1 中输入 15-(.75/12),在 RAD2 中输入 15-(1.75/12), 选择 APPLY,在 RAD1 中输入 15-(1.75/12), 在 RAD2 中输入 15-2/12,选择 OK；8.M ain Menu: Preprocessor-Modeling-Operate-Booleane-GlueArea,选择 PICK ALL ；9.M ain Menu: Preprocessor-Meshing-Size Contrls-Lines-Picked Lines,选择不锈钢层短边,在NDIV框中输入4,选择 APPLY,选择玻璃纤维层的短边,在 NDIV框中输入5,选择

35、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-Loa

36、ds-Apply-Thermal-ConvectionOn lines,选择不锈钢外壁,在VALI 框中输入80,在 VAL2I 框中输入44.5,选择 APPLY,选择铝层内壁,在 VALI 框中输入2.5,在 VAL2I 框中输入 70,选择 OK ；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南13.Main Menu: Solution-Solve-Curren

37、t 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 度,与管壁的对流换热系数随管壁温度而变。接管与罐为同一种材料，它的热物理性能如下表所示：温度70 200 300 400 500 oF 密度

38、0.285 0.285 0.285 0.285 0.285 lbm/in3导热系数8.35 8.90 9.35 9.8 10.23 Btu/hr-ft-oF 比热0.113 0.117 0.119 0.122 0.125 Btu/lbm-oF 对流系数 * 426 405 352 275 221 Btu/hr-ft2-oF *接管内壁对流系数求罐与接管的温度分布。以下分别列出LOG 文件及菜单操作/prep7 /title,Steady-state thermal analysis of pipe junction /units,bin ！使用英制单位et,1,90 ！定义热单元mp,dens

39、,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 ！比热名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南mpdata,hf,2,426/1

40、44,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,

41、1 ！定义显示角度/type,4 /title, V olumes used in building pipe/tank junction vplot ！显示实体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 modeled aplot /noerase lplot /erase accat,al

42、l ！组合罐远端的面及线，为映射划分网！格作准备lccat,12,7 lccat,10,5 lesize,20,4 ！在接管壁厚方向分4 等分lesize,40,6 ！在接管长度方向分6 等分lesize,6,4 ！在罐壁厚方向分4 等分allsel ！选择 EVERYTHING esize,0.4 ！设定默认的单元大小mshape,0,3d ！选择 3D 映射网格mshkey,1 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 36 页 - - - - - - - -

43、 - No BoundariesANSYS热 分 析 指 南save ！保存数据文件vmesh,all ！划分网格，产生节点与单元/pnum,defa /title, elements in portion being modeled eplot ！显示单元finish ！加载求解/solu antype,static ！定义为稳态分析nropt,auto ！设置求解选项为Program-chosen ！Newton-Raphson tunif,450 ！设定初始所有节点温度csys,1 ！变为柱坐标nsel,s,loc,x,ri1 ！选择罐内表面的节点sf,all,conv,250/144,

44、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 conditions nplot ！显示节点wpstyl

45、e,defa ！工作平面恢复默认状态csys,0 ！变为直角坐标autots,on ！打开自动步厂长nsubst,50 ！设定子步数量kbc,0 ！设定为阶越outpr,nsol,last ！设置输出solve ！进行求解finish ！进入后处理/post1 /title,Temperature contrours at pipe/tank junction plnsol,temp ！显示温度彩色云图finish /exit,all 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第

46、14 页，共 36 页 - - - - - - - - - No BoundariesANSYS热 分 析 指 南菜单操作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 MenuPre

47、processorMaterial 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/

48、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-

49、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 inne

50、r 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 MenuPlo