ANSYS热分析详解解析.pdf
《ANSYS热分析详解解析.pdf》由会员分享,可在线阅读,更多相关《ANSYS热分析详解解析.pdf(36页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、 第一章 简 介 一、热分析的目的 热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量等。热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。二、ANSYS的热分析 在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED 五种产品中包含热分析功能,其中 ANSYS/FLOTRAN不含相变热分析。ANSYS 热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。ANSYS 热分析包括热传导、热
2、对流及热辐射三种热传递方式。此外,还可以分析相变、有内热源、接触热阻等问题。三、ANSYS 热分析分类 稳态传热:系统的温度场不随时间变化 瞬态传热:系统的温度场随时间明显变化)四、耦合分析 热结构耦合 热流体耦合 热电耦合 热磁耦合 热电磁结构耦合等 第二章 基础知识 一、符号与单位 项目 国际单位 英制单位 ANSYS 代号 长度 m ft 时间 s s 质量 Kg lbm 温度 oF 力 N lbf 能量(热量)J BTU 功率(热流率)W BTU/sec 热流密度 W/m2 BTU/sec-ft2 生热速率 W/m3 BTU/sec-ft3 导热系数 W/m-BTU/sec-ft-oF
3、 KXX 对流系数 W/m2-BTU/sec-ft2-oF HF 密度 Kg/m3 lbm/ft3 DENS 比热 J/Kg-BTU/lbm-oF C 焓 J/m3 BTU/ft3 ENTH 二、传热学经典理论回顾 热分析遵循热力学第一定律,即能量守恒定律:对于一个封闭的系统(没有质量的流入或流出 PEKEUWQ 式中:Q 热量;W 作功;U系统内能;KE系统动能;【PE系统势能;对于大多数工程传热问题:0PEKE;通常考虑没有做功:0W,则:UQ;对于稳态热分析:0UQ,即流入系统的热量等于流出的热量;对于瞬态热分析:dtdUq,即流入或流出的热传递速率 q 等于系统内能的变化。三、热传递的
4、方式 1、热传导 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温 度梯度而引起的内能的交换。热传导遵循付里叶定律:dxdTkq,式中q为热流密度(W/m2),k为导热系数(W/m-),“-”表示热量流向温度降低的方向。2、热对流 热对流是指固体的表面与它周围接触的流体之间,由于温差的存在引起的热量的交换。热对流可以分为两类:自然对流和强制对流。热对流用牛顿冷却方程来描述:)(BSTThq,式中 h 为对流换热系数(或称膜传热系数、给热系数、膜系数等),TS为固体表面的温度,TB为周围流体的温度。3、热辐射-热辐射指物体发射电磁能,并被其它物体吸收转变为热的热量交换过程。物
5、体温度越高,单位时间辐射的热量越多。热传导和热对流都需要有传热介质,而热辐射无须任何介质。实质上,在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射,系统中每个物体同时辐射并吸收热量。它们之 间的净热量 传递可以 用斯蒂芬 波尔兹曼 方程来计算:qA FTT1121424(),式中q为热流率,为辐射率(黑度),为斯蒂芬波尔兹曼常数,约为10-8W/,A1为辐射面 1 的面积,F12为由辐射面 1 到辐射面 2 的形状系数,T1为辐射面 1 的绝对温度,T2为辐射面 2 的绝对温度。由上式可以看出,包含热辐射的热分析是高度非线性的。四、稳态传热 如果系统的净热流率为,即流入系
6、统的热量加上系统自身产生的热量等于流出系统的热量:q流入+q生成-q流出=0,则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为(以矩阵形式表示)K TQ 式中:K为传导矩阵,包含导热系数、对流系数及辐射率和形状系数;T为节点温度向量;Q为节点热流率向量,包含热生成;ANSYS 利用模型几何参数、材料热性能参数以及所施加的边界条件,生成 K、T以及 Q。五、瞬态传热 瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理,瞬态热平衡可以表达为(以矩阵形式表示):C TKTQ 式中:
7、K为传导矩阵,包含导热系数、对流系数及辐射率和形状系数;C为比热矩阵,考虑系统内能的增加;T为节点温度向量;T为温度对时间的导数;Q为节点热流率向量,包含热生成。六、线性与非线性 如果有下列情况产生,则为非线性热分析:、材料热性能随温度变化,如 K(T),C(T)等;、边界条件随温度变化,如 h(T)等;、含有非线性单元;、考虑辐射传热 非线性热分析的热平衡矩阵方程为:C TTK TTQ T 七、边界条件、初始条件 ANSYS 热分析的边界条件或初始条件可分为七种:温度、热流率、热流密度、对流、辐射、绝热、生热。八、热分析误差估计 仅用于评估由于网格密度不够带来的误差;仅适用于 SOLID 或
8、 SHELL 的热单元(只有温度一个自由度);基于单元边界的热流密度的不连续;仅对一种材料、线性、稳态热分析有效;使用自适应网格划分可以对误差进行控制。!第三章 稳态传热分析 一、稳态传热的定义 稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前,进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数 二、热分析的单元 热分析涉及到的单元有大约 40 种,其中纯粹用于热分析的有 14 种:线性:LINK32 两维二节点热传导单元 LINK33 三维二节点热传导单元 LINK34 二节点热对流单元 LI
9、NK31 二节点热辐射单元 二维实体:PLANE55 四节点四边形单元|PLANE77 八节点四边形单元 PLANE35 三节点三角形单元 PLANE75 四节点轴对称单元 PLANE78 八节点轴对称单元 三维实体 SOLID87 六节点四面体单元 SOLID70 八节点六面体单元 SOLID90 二十节点六面体单元 壳 SHELL57 四节点 点 MASS71 有关单元的详细解释,请参阅ANSYS Element Reference Guide 三、ANSYS稳态热分析的基本过程、ANSYS 热分析可分为三个步骤:前处理:建模 求解:施加载荷计算 后处理:查看结果 1、建模、确定 jobn
10、ame、title、unit;、进入 PREP7 前处理,定义单元类型,设定单元选项;、定义单元实常数;、定义材料热性能参数,对于稳态传热,一般只需定义导热系数,它可以是恒定的,也可以随温度变化;、创建几何模型并划分网格,请参阅ANSYS Modeling and Meshing Guide。2、施加载荷计算、。、定义分析类型 如果进行新的热分析:Command:ANTYPE,STATIC,NEW GUI:Main menuSolution-Analysis Type-New AnalysisSteady-state 如果继续上一次分析,比如增加边界条件等:Command:ANTYPE,STA
11、TIC,REST GUI:Main menuSolutionAnalysis Type-Restart、施加载荷 可以直接在实体模型或单元模型上施加五种载荷(边界条件):a、恒定的温度 通常作为自由度约束施加于温度已知的边界上。,Command Family:D GUI:Main MenuSolution-Loads-Apply-Thermal-Temperature b、热流率 热流率作为节点集中载荷,主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷),如果输入的值为正,代表热流流入节点,即单元获取热量。如果温度与热流率同时施加在一节点上则 ANSYS 读取温度值进行计算。注意
12、:如果在实体单元的某一节点上施加热流率,则此节点周围的单元要密一些,在两种导热系数差别很大的两个单元的公共节点上施加热流率时,尤其要注意。此外,尽可能使用热生成或热流密度边界条件,这样结果会更精确些。Command Family:F GUI:Main MenuSolution-Loads-Apply-Thermal-Heat Flow c、对流 对流边界条件作为面载施加于实体的外表面,计算与流体的热交换,它仅可施加于实体和壳模型上,对于线模型,可以通过对流线单元 LINK34 考虑对流。Command Family:SF GUI:Main MenuSolution-Loads-Apply-Th
13、ermal-Convection)d、热流密度 热流密度也是一种面载。当通过单位面积的热流率已知或通过 FLOTRAN CFD 计算得到时,可以在模型相应的外表面施加热流密度。如果输入的值为正,代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面,但 ANSYS 仅读取最后施加的面载进行计算。Command Family:F GUI:Main MenuSolution-Loads-Apply-Thermal-Heat Flux e、生热率 生热率作为体载施加于单元上,可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。Command Family:BF
14、GUI:Main MenuSolution-Loads-Apply-Thermal-Heat Generat GUI:Main MenuSolution-Load Step Opts-NolinearPredictor c.输出控制 控制打印输出:本选项可将任何结果数据输出到*.out 文件中。Command:OUTPR*GUI:Main MenuSolution-Load Step Opts-Output CtrlsSolu Printout 控制结果文件:控制*.rth 的内容。Command:OUTRES GUI:Main MenuSolution-Load Step Opts-Outp
15、ut CtrlsDB/Results File、确定分析选项 a.Newton-Raphson 选项(仅对非线性分析有用)Command:NROPT GUI:Main MenuSolutionAnalysis Options b.选择求解器:可选择如下求解器中一个进行求解:Frontal solver(默认)Jacobi Conjugate Gradient(JCG)solver *JCG out-of-memory solver Incomplete Cholesky Conjugate Gradient(ICCG)solver Pre-Conditioned Conjugate Gradi
16、ent Solver(PCG)Iterative(automatic solver selection option)Command:EQSLV GUI:Main MenuSolutionAnalysis Options 注意:热分析可选用 Iterative 选项进行快速求解,但如下情况除外:热分析包含 SURF19 或 SURF22 或超单元;热辐射分析;相变分析 需要 restart an analysis c.:d.确定绝对零度:在进行热辐射分析时,要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度,此值应设定为 273;如果使用的是华氏度,则为 460。Command:TOF
17、FST GUI:Main MenuSolutionAnalysis Options、保存模型:点击 ANSYS 工具条 SAVE_DB。、求解 Command:SOLVE GUI:Main MenuSolutionCurrent LS 3、后处理 ANSYS 将热分析的结果写入*.rth 文件中,它包含如下数据:基本数据:节点温度%导出数据:节点及单元的热流密度 节点及单元的热梯度 单元热流率 节点的反作用热流率 其它 对于稳态热分析,可以使用 POST1 进行后处理,关于后处理的完整描述,可参阅ANSYS Basic Analysis Procedures Guide。进入 POST1 后,
18、读入载荷步和子步: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:PLVECT GUI:Main MenuGeneral PostprocPlot ResultsPre-defined or Userdefined 列表显示 C
19、ommand:PRNSOL,PRESOL,PRRSOL等 GUI:Main MenuGeneral PostprocList ResultsNodal Solu,Element Solu,Reaction Solu 详细过程请参阅ANSYS Basic Analysis Procedures Guide。实例1:某一潜水艇可以简化为一圆筒,它由三层组成,最外面一层为不锈钢,中间为玻纤隔热层,最里面为铝层,筒内为空气,筒外为海水,求内外壁面温度及温度分布。几何参数:筒外径 30 feet 总壁厚 2 inch 不锈钢层壁厚 inch 玻纤层壁厚 1 inch 铝层壁厚 inch 筒长 200 f
20、eet 导热系数 不锈钢 BTU/玻纤 BTU/铝 BTU/边界条件 空气温度 70 oF 海水温度 oF 空气对流系数 BTU/海水对流系数 80 BTU/沿垂直于圆筒轴线作横截面,得到一圆环,取其中 1 度进行分析,如图示。不锈钢玻璃纤维铝R15 feet海水空气(以下分别列出 log 文件和菜单文件。/filename,Steady1/title,Steady-state thermal analysis of submarine /units,BFT Ro=15 !外径(ft)Rss=15-12)!不锈钢层内径 ft)Rins=15-12)!玻璃纤维层内径(ft)Ral=15-(2/1
21、2)!铝层内径(ft)Tair=70 !潜水艇内空气温度 Tsea=!海水温度 Kss=!不锈钢的导热系数(BTU/!玻璃纤维的导热系数 (BTU/!铝的导热系数(BTU/!空气的对流系数(BTU/!海 水 的 对 流 系 数(BTU/!定义二维热单元 mp,kxx,1,Kss !设定不锈钢的导热系数 mp,kxx,2,Kins !设定玻璃纤维的导热系数 mp,kxx,3,Kal !设定铝的导热系数 pcirc,Ro,Rss,!创建几何模型 pcirc,Rss,Rins,pcirc,Rins,Ral,aglue,all numcmp,area lesize,1,16 !设定划分网格密度 les
22、ize,4,4】lesize,14,5 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,TSEA!施加海水对流边界 SOLVE /POST1 PLNSOL !输出温度彩色云图 finish 菜单操作:1.Utility MenuFilechange jobename,输入 Steady1;2.Utility MenuFilechange title,输 入Steady-state therm
23、al analysis of submarine;3.4.在命令行输入:/units,BFT;5.Main Menu:Preprocessor;6.Main Menu:PreprocessorElement TypeAdd/Edit/Delete,选择 PLANE55;7.Main Menu:PreprocessorMaterial Prop-Constant-Isotropic,默认材料编号为 1,在KXX 框中输入,选择 APPLY,输入材料编号为 2,在 KXX 框中输入,选择 APPLY,输入材料编号为 3,在 KXX 框中输入;8.Main Menu:Preprocessor-Mod
24、eling-Create-Areas-CircleBy Dimensions,在 RAD1中输入 15,在 RAD2 中输入 15-(.75/12),在 THERA1 中输入,在 THERA2 中输入,选择 APPLY,在 RAD1 中输入 15-(.75/12),在 RAD2 中输入 15-12),选择 APPLY,在 RAD1 中输入 15-12),在 RAD2 中输入 15-2/12,选择 OK;9.Main Menu:Preprocessor-Modeling-Operate-Booleane-GlueArea,选择 PICK ALL;10.Main Menu:Preprocessor
25、-Meshing-Size Contrls-Lines-Picked Lines,选择不锈钢层短边,在 NDIV 框中输入 4,选择 APPLY,选择玻璃纤维层的短边,在 NDIV 框中输入 5,选择APPLY,选择铝层的短边,在NDIV框中输入2,选择APPLY,选择四个长边,在NDIV中输入16;11.Main Menu:Preprocessor-Attributes-DefinePicked Area,选择不锈钢层,在 MAT框中输入 1,选择 APPLY,选择玻璃纤维层,在 MAT 框中输入 2,选择 APPLY,选择铝层,在 MAT框中输入 3,选择 OK;12.Main Menu:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 分析 详解 解析
限制150内