ANSYS热分析详解解析.doc
《ANSYS热分析详解解析.doc》由会员分享,可在线阅读,更多相关《ANSYS热分析详解解析.doc(202页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-dateANSYS热分析详解解析第一章第一章简 介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量等。热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。二、ANSYS的热分析 在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal
2、、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。 ANSYS热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。 ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外,还可以分析相变、有内热源、接触热阻等问题。三、ANSYS 热分析分类 稳态传热:系统的温度场不随时间变化 瞬态传热:系统的温度场随时间明显变化四、耦合分析 热结构耦合 热流体耦合 热电耦合 热磁耦合 热电磁结构耦合等第二章基础知识一、符号与单位项目国际单位英制单位ANSYS代号长度mft时间ss质量Kglbm温度oF
3、力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 对于稳态热分析:,即流入系统的
4、热量等于流出的热量;l 对于瞬态热分析:,即流入或流出的热传递速率q等于系统内能的变化。三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律:,式中为热流密度(W/m2),为导热系数(W/m-),“-”表示热量流向温度降低的方向。2、热对流热对流是指固体的表面与它周围接触的流体之间,由于温差的存在引起的热量的交换。热对流可以分为两类:自然对流和强制对流。热对流用牛顿冷却方程来描述:,式中h为对流换热系数(或称膜传热系数、给热系数、膜系数等),为固体表面的温度,为周围流体的温度。3、热辐射热辐射指物体发射电磁
5、能,并被其它物体吸收转变为热的热量交换过程。物体温度越高,单位时间辐射的热量越多。热传导和热对流都需要有传热介质,而热辐射无须任何介质。实质上,在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射,系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬波尔兹曼方程来计算:,式中为热流率,为辐射率(黑度),为斯蒂芬波尔兹曼常数,约为5.6710-8W/m2.K4,A1为辐射面1的面积,为由辐射面1到辐射面2的形状系数,为辐射面1的绝对温度,为辐射面2的绝对温度。由上式可以看出,包含热辐射的热分析是高度非线性的。四、稳态传热如果系统的净热流率为,即流入系统的热量加上系
6、统自身产生的热量等于流出系统的热量:q流入+q生成-q流出=0,则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为(以矩阵形式表示)式中:为传导矩阵,包含导热系数、对流系数及辐射率和形状系数;为节点温度向量;为节点热流率向量,包含热生成;ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件,生成、以及。五、瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理,瞬态热平衡可以表达为(以矩阵形式表示):式中:为传导矩阵,包含导热系数、对流系数及辐射率和形状系数;为比
7、热矩阵,考虑系统内能的增加;为节点温度向量;为温度对时间的导数;为节点热流率向量,包含热生成。六、线性与非线性如果有下列情况产生,则为非线性热分析:、材料热性能随温度变化,如K(T),C(T)等;、边界条件随温度变化,如h(T)等;、含有非线性单元;、 考虑辐射传热非线性热分析的热平衡矩阵方程为:七、边界条件、初始条件ANSYS热分析的边界条件或初始条件可分为七种:温度、热流率、热流密度、对流、辐射、绝热、生热。八、热分析误差估计 仅用于评估由于网格密度不够带来的误差; 仅适用于SOLID或SHELL的热单元(只有温度一个自由度); 基于单元边界的热流密度的不连续; 仅对一种材料、线性、稳态热
8、分析有效; 使用自适应网格划分可以对误差进行控制。第三章稳态传热分析一、稳态传热的定义稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前,进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数二、热分析的单元热分析涉及到的单元有大约40种,其中纯粹用于热分析的有14种:线性:LINK32两维二节点热传导单元LINK33三维二节点热传导单元LINK34二节点热对流单元LINK31二节点热辐射单元二维实体:PLANE55四节点四边形单元PLANE77八节点四边形单元PLANE35三节点三角形单元PLANE
9、75四节点轴对称单元PLANE78八节点轴对称单元三维实体SOLID87六节点四面体单元SOLID70八节点六面体单元SOLID90二十节点六面体单元壳SHELL57四节点点MASS71有关单元的详细解释,请参阅ANSYS Element Reference Guide三、ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤: 前处理:建模 求解:施加载荷计算 后处理:查看结果1、建模、确定jobname、title、unit;、进入PREP7前处理,定义单元类型,设定单元选项;、定义单元实常数;、定义材料热性能参数,对于稳态传热,一般只需定义导热系数,它可以是恒定的,也可以随温度变化
10、;、创建几何模型并划分网格,请参阅ANSYS Modeling and Meshing Guide。2、施加载荷计算、 定义分析类型l 如果进行新的热分析:Command: ANTYPE, STATIC, NEWGUI: Main menuSolution-Analysis Type-New AnalysisSteady-statel 如果继续上一次分析,比如增加边界条件等:Command: ANTYPE, STATIC, RESTGUI: Main menuSolutionAnalysis Type-Restart、 施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件) :a、恒定
11、的温度通常作为自由度约束施加于温度已知的边界上。Command Family:DGUI:Main MenuSolution-Loads-Apply-Thermal-Temperatureb、热流率热流率作为节点集中载荷,主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷),如果输入的值为正,代表热流流入节点,即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意:如果在实体单元的某一节点上施加热流率,则此节点周围的单元要密一些,在两种导热系数差别很大的两个单元的公共节点上施加热流率时,尤其要注意。此外,尽可能使用热生成或热流密度边界条件,这样结果会更
12、精确些。Command Family: FGUI:Main MenuSolution-Loads-Apply-Thermal-Heat Flowc、对流对流边界条件作为面载施加于实体的外表面,计算与流体的热交换,它仅可施加于实体和壳模型上,对于线模型,可以通过对流线单元LINK34考虑对流。Command Family: SFGUI:Main MenuSolution-Loads-Apply-Thermal-Convectiond、热流密度热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRAN CFD计算得到时,可以在模型相应的外表面施加热流密度。如果输入的值为正,代表热流流入单
13、元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面,但ANSYS仅读取最后施加的面载进行计算。Command Family: FGUI:Main MenuSolution-Loads-Apply-Thermal-Heat Fluxe、生热率生热率作为体载施加于单元上,可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。Command Family: BFGUI:Main MenuSolution-Loads-Apply-Thermal-Heat Generat、确定载荷步选项对于一个热分析,可以确定普通选项、非线性选项以及输出控制。a. 普通选项 时间选项:虽然对于
14、稳态热分析,时间选项并没有实际的物理意义,但它提供了一个方便的设置载荷步和载荷子步的方法。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-Lo
15、ad 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 MenuS
16、olution-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 求解结束选项:
17、如果在规定的迭代次数内,达不到收敛,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-L
18、oad 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 G
19、UI: 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: EQS
20、LV GUI: Main MenuSolutionAnalysis Options注意:热分析可选用Iterative选项进行快速求解,但如下情况除外: 热分析包含SURF19或SURF22或超单元; 热辐射分析; 相变分析 需要restart an analysisc. 确定绝对零度:在进行热辐射分析时,要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度,此值应设定为273;如果使用的是华氏度,则为460。Command: TOFFSTGUI: Main MenuSolutionAnalysis Options、 保存模型: 点击ANSYS工具条SAVE_DB。、求解Command:
21、 SOLVEGUI: Main MenuSolutionCurrent LS3、后处理ANSYS将热分析的结果写入*.rth文件中,它包含如下数据:基本数据: 节点温度导出数据: 节点及单元的热流密度 节点及单元的热梯度 单元热流率 节点的反作用热流率 其它对于稳态热分析,可以使用POST1进行后处理,关于后处理的完整描述,可参阅ANSYS Basic Analysis Procedures Guide。进入POST1后,读入载荷步和子步:Command: SET GUI: Main MenuGeneral Postproc-Read Results-By Load Step可以通过如下三种方
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 分析 详解 解析
限制150内