1、ANSYS的热分析第一章 简介一、热分析之目的热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量等。热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。二、ANSYS之热分析 在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。 ANSYS热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。 ANSYS热分析包括热传
2、导、热对流及热辐射三种热传递方式。此外,还可以分析相变、有内热源、接触热阻等问题。三、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密
3、度Kg/m3lbm/ft3DENS比热J/Kg-BTU/lbm-oFC焓J/m3BTU/ft3ENTH二、传热学经典理论回顾热分析遵循热力学第一定律,即能量守恒定律: 对于一个封闭的系统(没有质量的流入或流出式中: Q热量;W作功;系统内能;系统动能;系统势能; 对于大多数工程传热问题:; 通常考虑没有做功:,则:; 对于稳态热分析:,即流入系统的热量等于流出的热量; 对于瞬态热分析:,即流入或流出的热传递速率q等于系统内能的变化。三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律:,式中为热流密度(W/
4、m2),为导热系数(W/m-),“-”表示热量流向温度降低的方向。2、热对流热对流是指固体的表面与它周围接触的流体之间,由于温差的存在引起的热量的交换。热对流可以分为两类:自然对流和强制对流。热对流用牛顿冷却方程来描述:,式中h为对流换热系数(或称膜传热系数、给热系数、膜系数等),为固体表面的温度,为周围流体的温度。3、热辐射热辐射指物体发射电磁能,并被其它物体吸收转变为热的热量交换过程。物体温度越高,单位时间辐射的热量越多。热传导和热对流都需要有传热介质,而热辐射无须任何介质。实质上,在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射,系统中每个物体同时辐射并吸收热量。它
5、们之间的净热量传递可以用斯蒂芬波尔兹曼方程来计算:,式中为热流率,为辐射率(黑度),为斯蒂芬波尔兹曼常数,约为5.6710-8W/m2.K4,A1为辐射面1的面积,为由辐射面1到辐射面2的形状系数,为辐射面1的绝对温度,为辐射面2的绝对温度。由上式可以看出,包含热辐射的热分析是高度非线性的。四、稳态传热如果系统的净热流率为,即流入系统的热量加上系统自身产生的热量等于流出系统的热量:q流入+q生成-q流出=0,则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为(以矩阵形式表示)式中:为传导矩阵,包含导热系数、对流系数及辐射率和形状系数; 为节点温度向量; 为节
6、点热流率向量,包含热生成;ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件,生成、以及。五、瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理,瞬态热平衡可以表达为(以矩阵形式表示):式中: 为传导矩阵,包含导热系数、对流系数及辐射率和形状系数; 为比热矩阵,考虑系统内能的增加; 为节点温度向量; 为温度对时间的导数; 为节点热流率向量,包含热生成。六、线性与非线性如果有下列情况产生,则为非线性热分析:、材料热性能随温度变化,如K(T),C(T)等;、边界条件随温度变化,如h(T)等;、含有
7、非线性单元;4、 考虑辐射传热非线性热分析的热平衡矩阵方程为:七、边界条件、初始条件ANSYS热分析的边界条件或初始条件可分为七种:温度、热流率、热流密度、对流、辐射、绝热、生热。八、热分析误差估计 仅用于评估由于网格密度不够带来的误差; 仅适用于SOLID或SHELL的热单元(只有温度一个自由度); 基于单元边界的热流密度的不连续; 仅对一种材料、线性、稳态热分析有效; 使用自适应网格划分可以对误差进行控制。第三章稳态传热分析一、稳态传热的定义稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前,进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定
8、的热载荷引起的温度、热梯度、热流率、热流密度等参数二、热分析的单元热分析涉及到的单元有大约40种,其中纯粹用于热分析的有14种:线性: LINK32 两维二节点热传导单元 LINK33 三维二节点热传导单元 LINK34 二节点热对流单元 LINK31 二节点热辐射单元二维实体: PLANE55 四节点四边形单元 PLANE77 八节点四边形单元 PLANE35 三节点三角形单元 PLANE75 四节点轴对称单元 PLANE78 八节点轴对称单元三维实体 SOLID87 六节点四面体单元 SOLID70 八节点六面体单元 SOLID90 二十节点六面体单元壳 SHELL57 四节点点 MASS
9、71 有关单元的详细解释,请参阅ANSYSElementReferenceGuide三、ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤: 前处理: 建模 求解: 施加载荷计算 后处理: 查看结果1、建模、确定jobname、title、unit;、进入PREP7前处理,定义单元类型,设定单元选项;、定义单元实常数;、定义材料热性能参数,对于稳态传热,一般只需定义导热系数,它可以是恒定的,也可以随温度变化;、创建几何模型并划分网格,请参阅ANSYSModelingandMeshingGuide。2、施加载荷计算1、 定义分析类型 如果进行新的热分析:Command:ANTYPE,S
10、TATIC,NEWGUI:MainmenuSolution-AnalysisType-NewAnalysisSteady-state 如果继续上一次分析,比如增加边界条件等:Command:ANTYPE,STATIC,RESTGUI:MainmenuSolutionAnalysisType-Restart2、 施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件):a、恒定的温度通常作为自由度约束施加于温度已知的边界上。CommandFamily:DGUI:MainMenuSolution-Loads-Apply-Thermal-Temperatureb、热流率热流率作为节点集中载荷,
11、主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷),如果输入的值为正,代表热流流入节点,即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意:如果在实体单元的某一节点上施加热流率,则此节点周围的单元要密一些,在两种导热系数差别很大的两个单元的公共节点上施加热流率时,尤其要注意。此外,尽可能使用热生成或热流密度边界条件,这样结果会更精确些。CommandFamily:FGUI:MainMenuSolution-Loads-Apply-Thermal-HeatFlowc、对流对流边界条件作为面载施加于实体的外表面,计算与流体的热交换,它仅可施加于实体
12、和壳模型上,对于线模型,可以通过对流线单元LINK34考虑对流。CommandFamily:SFGUI:MainMenuSolution-Loads-Apply-Thermal-Convectiond、热流密度热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRANCFD计算得到时,可以在模型相应的外表面施加热流密度。如果输入的值为正,代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面,但ANSYS仅读取最后施加的面载进行计算。CommandFamily:FGUI:MainMenuSolution-Loads-Apply-Thermal-HeatF
13、luxe、生热率生热率作为体载施加于单元上,可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。CommandFamily:BFGUI:MainMenuSolution-Loads-Apply-Thermal-HeatGenerat、确定载荷步选项对于一个热分析,可以确定普通选项、非线性选项以及输出控制。a.普通选项 时间选项:虽然对于稳态热分析,时间选项并没有实际的物理意义,但它提供了一个方便的设置载荷步和载荷子步的方法。Command:TIMEGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTime-TimeStep/TimeandSu
14、bstps 每载荷步中子步的数量或时间步大小:对于非线性分析,每一载荷步需要多个子步。Command:NSUBSTGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTimeandSubstpsCommand:DELTIMGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTime-TimeStep 递进或阶越选项:如果定义阶越(stepped)选项,载荷值在这个载荷步内保持不变;如果为递进(ramped)选项,则载荷值由上一载荷步值到本载荷步值随每一子步线性变化。Command:KBCGUI:MainMenu
15、Solution-LoadStepOpts-Time/FrequencTime-TimeStep/TimeandSubstpsb.非线性选项 迭代次数:本选项设置每一子步允许的最多的迭代次数。默认值为25,对大数热分析问题足够。Command:NEQITGUI:MainMenuSolution-LoadStepOpts-NolinearEquilibriumIter 自动时间步长:对于非线性问题,可以自动设定子步间载荷的增长,保证求解的稳定性和准确性。Command:AUTOTSGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTime-TimeSt
16、ep/TimeandSubstps 收敛误差:可根据温度、热流率等检验热分析的收敛性。Command:CNVTOLGUI:MainMenuSolution-LoadStepOpts-NolinearConvergenceCrit 求解结束选项:如果在规定的迭代次数内,达不到收敛,ANSYS可以停止求解或到下一载荷步继续求解。Command:NCNVGUI:MainMenuSolution-LoadStepOpts-NolinearCriteriatoStop 线性搜索:设置本选项可使ANSYS用Newton-Raphson方法进行线性搜索。Command:LNSRCHGUI:MainMenuS
17、olution-LoadStepOpts-NolinearLineSearch 预测矫正:本选项可激活每一子步第一次迭代对自由度求解的预测矫正。Command:PREDGUI:MainMenuSolution-LoadStepOpts-NolinearPredictorc. 输出控制 控制打印输出:本选项可将任何结果数据输出到*.out文件中。Command:OUTPRGUI:MainMenuSolution-LoadStepOpts-OutputCtrlsSoluPrintout 控制结果文件:控制*.rth的内容。Command:OUTRESGUI:MainMenuSolution-Loa
18、dStepOpts-OutputCtrlsDB/ResultsFile、确定分析选项a. Newton-Raphson选项(仅对非线性分析有用)Command:NROPTGUI:MainMenuSolutionAnalysisOptionsb. 选择求解器:可选择如下求解器中一个进行求解: Frontalsolver(默认) JacobiConjugateGradient(JCG)solver JCGout-of-memorysolver IncompleteCholeskyConjugateGradient(ICCG)solver Pre-ConditionedConjugateGradie
19、ntSolver(PCG) Iterative(automaticsolverselectionoption)Command:EQSLVGUI:MainMenuSolutionAnalysisOptions注意:热分析可选用Iterative选项进行快速求解,但如下情况除外: 热分析包含SURF19或SURF22或超单元; 热辐射分析; 相变分析 需要restartananalysisc. 确定绝对零度:在进行热辐射分析时,要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度,此值应设定为273;如果使用的是华氏度,则为460。Command:TOFFSTGUI:MainMenuSolu
20、tionAnalysisOptions5、 保存模型:点击ANSYS工具条SAVE_DB。、求解Command:SOLVEGUI:MainMenuSolutionCurrentLS3、后处理ANSYS将热分析的结果写入*.rth文件中,它包含如下数据:基本数据: 节点温度导出数据: 节点及单元的热流密度 节点及单元的热梯度 单元热流率 节点的反作用热流率 其它对于稳态热分析,可以使用POST1进行后处理,关于后处理的完整描述,可参阅ANSYSBasicAnalysisProceduresGuide。进入POST1后,读入载荷步和子步:Command:SETGUI:MainMenuGeneral
21、Postproc-ReadResults-ByLoadStep可以通过如下三种方式查看结果: 彩色云图显示Command:PLNSOL,PLESOL,PLETAB等GUI:MainMenuGeneralPostprocPlotResultsNodalSolu,ElementSolu,ElemTable 矢量图显示Command:PLVECTGUI:MainMenuGeneralPostprocPlotResultsPre-definedorUserdefined 列表显示Command:PRNSOL,PRESOL,PRRSOL等GUI:MainMenuGeneralPostprocListRe
22、sultsNodalSolu,ElementSolu,ReactionSolu详细过程请参阅ANSYSBasicAnalysisProceduresGuide。实例1:某一潜水艇可以简化为一圆筒,它由三层组成,最外面一层为不锈钢,中间为玻纤隔热层,最里面为铝层,筒内为空气,筒外为海水,求内外壁面温度及温度分布。几何参数: 筒外径 30 feet 总壁厚 2 inch 不锈钢层壁厚 0.75 inch 玻纤层壁厚 1 inch 铝层壁厚 0.25 inch 筒长 200 feet导热系数 不锈钢 8.27 BTU/hr.ft.oF 玻纤 0.028 BTU/hr.ft.oF 铝 117.4 BT
23、U/hr.ft.oF边界条件 空气温度 70 oF 海水温度 44.5 oF 空气对流系数 2.5 BTU/hr.ft2.oF 海水对流系数 80 BTU/hr.ft2.oF沿垂直于圆筒轴线作横截面,得到一圆环,取其中1度进行分析,如图示。以下分别列出log文件和菜单文件。/filename,Steady1/title,Steady-statethermalanalysisofsubmarine/units,BFTRo=15 !外径(ft)Rss=15-(0.75/12) !不锈钢层内径ft)Rins=15-(1.75/12) !玻璃纤维层内径(ft)Ral=15-(2/12) !铝层内径(f
24、t)Tair=70 !潜水艇内空气温度Tsea=44.5 !海水温度Kss=8.27 !不锈钢的导热系数Kins=0.028 !玻璃纤维的导热系数Kal=117.4 !铝的导热系数Hair=2.5 !空气的对流系数Hsea=80 !海水的对流系数/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.5ag
25、lue,all numcmp,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. UtilityMenuFilechangejobename,输入Steady1;2. UtilityMenuFilechangeti
26、tle,输入Steady-statethermalanalysisofsubmarine;3. 在命令行输入:/units,BFT;4. MainMenu:Preprocessor;5. MainMenu:PreprocessorElementTypeAdd/Edit/Delete,选择PLANE55;6. MainMenu:PreprocessorMaterialProp-Constant-Isotropic,默认材料编号为1,在KXX框中输入8.27,选择APPLY,输入材料编号为2,在KXX框中输入0.028,选择APPLY,输入材料编号为3,在KXX框中输入117.4;7. MainM
27、enu:Preprocessor-Modeling-Create-Areas-CircleByDimensions,在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. MainMenu:Preprocessor-Modeling-Operate-Booleane-GlueArea,选择PICKALL;9. MainMenu:P
28、reprocessor-Meshing-SizeContrls-Lines-PickedLines,选择不锈钢层短边,在NDIV框中输入4,选择APPLY,选择玻璃纤维层的短边,在NDIV框中输入5,选择APPLY,选择铝层的短边,在NDIV框中输入2,选择APPLY,选择四个长边,在NDIV中输入16;10. MainMenu:Preprocessor-Attributes-DefinePickedArea,选择不锈钢层,在MAT框中输入1,选择APPLY,选择玻璃纤维层,在MAT框中输入2,选择APPLY,选择铝层,在MAT框中输入3,选择OK;11. MainMenu:Preproces
29、sor-Meshing-Mesh-Areas-Mapped3or4sided,选择PICKALL;12. MainMenu:Solution-Loads-Apply-Thermal-ConvectionOnlines,选择不锈钢外壁,在VALI框中输入80,在VAL2I框中输入44.5,选择APPLY,选择铝层内壁,在VALI框中输入2.5,在VAL2I框中输入70,选择OK;13. MainMenu:Solution-Solve-CurrentLS;14. MainMenu:GeneralPostprocPlotResults-ContourPlot-NodalSolu,选择Temperature。实例2一圆筒形的罐有一接管,罐外径为3英尺,壁厚为0.2英尺,接管外径为0.5英尺,壁厚为0.1英尺,罐与接管的轴线垂直且接管远离罐的端部。如图所示:罐内流体温度为华氏450度,与罐壁的对流换热系数年为250BUT/hr-ft2-oF,接管内流体的温度为华氏100度
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1