热分析.DOC_Ansys中文手册_doc_大学课件预览_高等教育资讯网
Ansys中文手册:热分析.DOC
分类: 格式: 日期:日
第一章 简 介一、热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量〕等。热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。二、ANSYS的热分析在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。ANSYS热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外,还可以分析相变、有内热源、接触热阻等问题。三、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密度Kg/m3lbm/ft3DENS比热J/Kg-℃BTU/lbm-oFC焓J/m3BTU/ft3ENTH二、传热学经典理论回顾热分析遵循热力学第一定律,即能量守恒定律:对于一个封闭的系统(没有质量的流入或流出〕式中: Q ―― 热量;W ―― 作功;――系统内能;――系统动能;――系统势能;对于大多数工程传热问题:;通常考虑没有做功:,则:;对于稳态热分析:,即流入系统的热量等于流出的热量;对于瞬态热分析:,即流入或流出的热传递速率q等于系统内能的变化。三、热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律:,式中为热流密度(W/m2),为导热系数(W/m-℃),“-”表示热量流向温度降低的方向。2、热对流热对流是指固体的表面与它周围接触的流体之间,由于温差的存在引起的热量的交换。热对流可以分为两类:自然对流和强制对流。热对流用牛顿冷却方程来描述:,式中h为对流换热系数(或称膜传热系数、给热系数、膜系数等),为固体表面的温度,为周围流体的温度。3、热辐射热辐射指物体发射电磁能,并被其它物体吸收转变为热的热量交换过程。物体温度越高,单位时间辐射的热量越多。热传导和热对流都需要有传热介质,而热辐射无须任何介质。实质上,在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射,系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬―波尔兹曼方程来计算:,式中为热流率,为辐射率(黑度),为斯蒂芬-波尔兹曼常数,约为5.67×10-8W/m2.K4,A1为辐射面1的面积,为由辐射面1到辐射面2的形状系数,为辐射面1的绝对温度,为辐射面2的绝对温度。由上式可以看出,包含热辐射的热分析是高度非线性的。四、稳态传热如果系统的净热流率为0,即流入系统的热量加上系统自身产生的热量等于流出系统的热量:q流入+q生成-q流出=0,则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳态热分析的能量平衡方程为(以矩阵形式表示)式中:为传导矩阵,包含导热系数、对流系数及辐射率和形状系数; 为节点温度向量; 为节点热流率向量,包含热生成;ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件,生成、以及。五、瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理,瞬态热平衡可以表达为(以矩阵形式表示):式中: 为传导矩阵,包含导热系数、对流系数及辐射率和形状系数; 为比热矩阵,考虑系统内能的增加; 为节点温度向量; 为温度对时间的导数; 为节点热流率向量,包含热生成。六、线性与非线性如果有下列情况产生,则为非线性热分析:①、材料热性能随温度变化,如K(T),C(T)等;②、边界条件随温度变化,如h(T)等;③、含有非线性单元;考虑辐射传热非线性热分析的热平衡矩阵方程为:七、边界条件、初始条件ANSYS热分析的边界条件或初始条件可分为七种:温度、热流率、热流密度、对流、辐射、绝热、生热。八、热分析误差估计仅用于评估由于网格密度不够带来的误差;仅适用于SOLID或SHELL的热单元(只有温度一个自由度);基于单元边界的热流密度的不连续;仅对一种材料、线性、稳态热分析有效;使用自适应网格划分可以对误差进行控制。第三章 稳态传热分析一、稳态传热的定义稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前,进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数二、热分析的单元热分析涉及到的单元有大约40种,其中纯粹用于热分析的有14种:线性:
两维二节点热传导单元 LINK33
三维二节点热传导单元 LINK34
二节点热对流单元 LINK31
二节点热辐射单元二维实体: PLANE55
四节点四边形单元
八节点四边形单元
三节点三角形单元
四节点轴对称单元
八节点轴对称单元三维实体 SOLID87
六节点四面体单元
八节点六面体单元
二十节点六面体单元壳
MASS71 有关单元的详细解释,请参阅《ANSYS Element Reference Guide》三、ANSYS稳态热分析的基本过程ANSYS热分析可分为三个步骤:前处理: 建模求解: 施加载荷计算后处理: 查看结果1、建模①、确定jobname、title、unit;②、进入PREP7前处理,定义单元类型,设定单元选项;③、定义单元实常数;④、定义材料热性能参数,对于稳态传热,一般只需定义导热系数,它可以是恒定的,也可以随温度变化;⑤、创建几何模型并划分网格,请参阅《ANSYS Modeling and Meshing Guide》。2、施加载荷计算定义分析类型如果进行新的热分析:Command,ANTYPE,STATIC,NEWGUI,Main menu&Solution&-Analysis Type-&New Analysis&Steady-state如果继续上一次分析,比如增加边界条件等:Command,ANTYPE,STATIC,RESTGUI,Main menu&Solution&Analysis Type-&Restart施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件),a、恒定的温度通常作为自由度约束施加于温度已知的边界上。Command Family: DGUI:Main Menu&Solution&-Loads-Apply&-Thermal-Temperatureb、热流率热流率作为节点集中载荷,主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷),如果输入的值为正,代表热流流入节点,即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意:如果在实体单元的某一节点上施加热流率,则此节点周围的单元要密一些,在两种导热系数差别很大的两个单元的公共节点上施加热流率时,尤其要注意。此外,尽可能使用热生成或热流密度边界条件,这样结果会更精确些。Command Family,FGUI:Main Menu&Solution&-Loads-Apply&-Thermal-Heat Flowc、对流对流边界条件作为面载施加于实体的外表面,计算与流体的热交换,它仅可施加于实体和壳模型上,对于线模型,可以通过对流线单元LINK34考虑对流。Command Family,SFGUI:Main Menu&Solution&-Loads-Apply&-Thermal-Convectiond、热流密度热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRAN CFD计算得到时,可以在模型相应的外表面施加热流密度。如果输入的值为正,代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面,但ANSYS仅读取最后施加的面载进行计算。Command Family,FGUI:Main Menu&Solution&-Loads-Apply&-Thermal-Heat Fluxe、生热率生热率作为体载施加于单元上,可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。Command Family,BFGUI:Main Menu&Solution&-Loads-Apply&-Thermal-Heat Generat③、确定载荷步选项对于一个热分析,可以确定普通选项、非线性选项以及输出控制。a,普通选项时间选项:虽然对于稳态热分析,时间选项并没有实际的物理意义,但它提供了一个方便的设置载荷步和载荷子步的方法。Command,TIMEGUI,Main Menu&Solution&-Load Step Opts-Time/Frequenc&Time-Time Step/Time and Substps每载荷步中子步的数量或时间步大小:对于非线性分析,每一载荷步需要多个子步。Command,NSUBSTGUI,Main Menu&Solution&-Load Step Opts-&Time/Frequenc&Time and SubstpsCommand,DELTIMGUI,Main Menu&Solution&-Load Step Opts-&Time/Frequenc&Time-Time Step递进或阶越选项:如果定义阶越(stepped)选项,载荷值在这个载荷步内保持不变;如果为递进(ramped)选项,则载荷值由上一载荷步值到本载荷步值随每一子步线性变化。Command,KBCGUI,Main Menu&Solution&-Load Step Opts-Time/Frequenc&Time-Time Step/Time and Substpsb,非线性选项迭代次数:本选项设置每一子步允许的最多的迭代次数。默认值为25,对大数热分析问题足够。Command,NEQITGUI,Main Menu&Solution&-Load Step Opts-Nolinear&Equilibrium Iter自动时间步长,对于非线性问题,可以自动设定子步间载荷的增长,保证求解的稳定性和准确性。Command,AUTOTSGUI,Main Menu&Solution&-Load Step Opts-Time/Frequenc&Time-Time Step/Time and Substps收敛误差:可根据温度、热流率等检验热分析的收敛性。Command,CNVTOLGUI,Main Menu&Solution&-Load Step Opts-Nolinear&Convergence Crit求解结束选项:如果在规定的迭代次数内,达不到收敛,ANSYS可以停止求解或到下一载荷步继续求解。Command,NCNVGUI,Main Menu&Solution&-Load Step Opts-Nolinear&Criteria to Stop线性搜索:设置本选项可使ANSYS用Newton-Raphson方法进行线性搜索。Command,LNSRCHGUI,Main Menu&Solution&-Load Step Opts-Nolinear&Line Search预测矫正:本选项可激活每一子步第一次迭代对自由度求解的预测矫正。Command,PREDGUI,Main Menu&Solution&-Load Step Opts-Nolinear&Predictor输出控制控制打印输出:本选项可将任何结果数据输出到*.out 文件中。Command,OUTPRGUI,Main Menu&Solution&-Load Step Opts-Output Ctrls&Solu Printout控制结果文件:控制*.rth的内容。Command,OUTRESGUI,Main Menu&Solution&-Load Step Opts-Output Ctrls&DB/Results File④、确定分析选项Newton-Raphson选项(仅对非线性分析有用)Command,NROPTGUI,Main Menu&Solution&Analysis Options选择求解器:可选择如下求解器中一个进行求解:Frontal solver(默认)Jacobi Conjugate Gradient(JCG) solverJCG out-of-memory solverIncomplete Cholesky Conjugate Gradient(ICCG) solverPre-Conditioned Conjugate Gradient Solver(PCG)Iterative(automatic solver selection option)Command,EQSLVGUI,Main Menu&Solution&Analysis Options注意:热分析可选用Iterative选项进行快速求解,但如下情况除外:热分析包含SURF19或SURF22或超单元;热辐射分析;相变分析需要restart an analysis确定绝对零度:在进行热辐射分析时,要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度,此值应设定为273;如果使用的是华氏度,则为460。Command,TOFFSTGUI,Main Menu&Solution&Analysis Options保存模型,点击ANSYS工具条SAVE_DB。⑥、求解Command,SOLVEGUI,Main Menu&Solution&Current LS3、后处理ANSYS将热分析的结果写入*.rth文件中,它包含如下数据:基本数据:节点温度导出数据:节点及单元的热流密度节点及单元的热梯度单元热流率节点的反作用热流率其它对于稳态热分析,可以使用POST1进行后处理,关于后处理的完整描述,可参阅《ANSYS Basic Analysis Procedures Guide》。进入POST1后,读入载荷步和子步:Command,SETGUI,Main Menu&General Postproc&-Read Results-By Load Step可以通过如下三种方式查看结果:彩色云图显示Command,PLNSOL,PLESOL,PLETAB等GUI,Main Menu&General Postproc&Plot Results&Nodal Solu,Element Solu,Elem Table矢量图显示Command,PLVECTGUI,Main Menu&General Postproc&Plot Results&Pre-defined or Userdefined列表显示Command,PRNSOL,PRESOL,PRRSOL等GUI,Main Menu&General Postproc&List Results&Nodal Solu,Element Solu,Reaction Solu详细过程请参阅《ANSYS Basic Analysis Procedures Guide》。实例1:某一潜水艇可以简化为一圆筒,它由三层组成,最外面一层为不锈钢,中间为玻纤隔热层,最里面为铝层,筒内为空气,筒外为海水,求内外壁面温度及温度分布。几何参数: 筒外径
不锈钢层壁厚 0.75 inch
玻纤层壁厚
200 feet导热系数 不锈钢
8.27 BTU/hr.ft.oF
0.028 BTU/hr.ft.oF
117.4 BTU/hr.ft.oF边界条件 空气温度
空气对流系数 2.5 BTU/hr.ft2.oF
海水对流系数 80 BTU/hr.ft2.oF沿垂直于圆筒轴线作横截面,得到一圆环,取其中1度进行分析,如图示。以下分别列出log文件和菜单文件。/filename,Steady1/title,Steady-state thermal analysis of submarine/units,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,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.5aglue,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菜单操作:Utility Menu&File&change jobename,输入Steady1;Utility Menu&File&change title,输入Steady-state thermal analysis of submarine;在命令行输入:/units,BFT;Main Menu,Preprocessor;Main Menu,Preprocessor&Element Type&Add/Edit/Delete,选择PLANE55;Main Menu,Preprocessor&Material Prop&-Constant-Isotropic,默认材料编号为1,在KXX框中输入8.27,选择APPLY,输入材料编号为2,在KXX框中输入0.028,选择APPLY,输入材料编号为3,在KXX框中输入117.4;Main Menu,Preprocessor&-Modeling-&Create&-Areas-Circle&By Dimensions,在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;Main Menu,Preprocessor&-Modeling-&Operate&-Booleane-&Glue&Area,选择PICK ALL;Main Menu,Preprocessor&-Meshing-Size Contrls&-Lines-Picked Lines,选择不锈钢层短边,在NDIV框中输入4,选择APPLY,选择玻璃纤维层的短边,在NDIV框中输入5,选择APPLY,选择铝层的短边,在NDIV框中输入2,选择APPLY,选择四个长边,在NDIV中输入16;Main Menu,Preprocessor&-Attributes-Define&Picked Area,选择不锈钢层,在MAT框中输入1,选择APPLY,选择玻璃纤维层,在MAT框中输入2,选择APPLY,选择铝层,在MAT框中输入3,选择OK;Main Menu,Preprocessor&-Meshing-Mesh&-Areas-Mapped&3 or 4 sided,选择PICK ALL;Main Menu,Solution&-Loads-Apply&-Thermal-Convection&On lines,选择不锈钢外壁,在VALI框中输入80,在VAL2I框中输入44.5,选择APPLY,选择铝层内壁,在VALI框中输入2.5,在VAL2I框中输入70,选择OK;Main Menu,Solution&-Solve-Current LS;Main Menu,General Postproc&Plot Results&-Contour Plot-Nodal Solu,选择Temperature。实例2一圆筒形的罐有一接管,罐外径为3英尺,壁厚为0.2英尺,接管外径为0.5英尺,壁厚为0.1英尺,罐与接管的轴线垂直且接管远离罐的端部。如图所示:罐内流体温度为华氏450度,与罐壁的对流换热系数年为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 thermal 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
!罐外半径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
!删除多余实体!划分网格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等分allsel
!选择EVERYTHINGesize,0.4
!设定默认的单元大小mshape,0,3d
!选择3D映射网格mshkey,1
!保存数据文件vmesh,all
!划分网格,产生节点与单元/pnum,defa/title,elements in portion being modeledeplot
!显示单元finish!加载求解/solu
antype,static
!定义为稳态分析nropt,auto
!设置求解选项为Program-chosen!Newton-Raphsontunif,450
!设定初始所有节点温度csys,1
!变为柱坐标nsel,s,loc,x,ri1
!选择罐内表面的节点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
!显示节点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菜单操作设定标题:Utility Menu&File&Change Title,输入Steady-State analysis of pipe junction,选择OK;设定单位制:在命令提示行输入/UNITS,BIN;定义单元类型:Main Menu&Preprocesor&Element Type&Add/Edit/Delete,选择Thermal Solid,Bricck 20 node 90号单元;定义材料属性Main Menu&Preprocessor&Material Props&-Constant-&Isotropic,默认材料编号1,在DENSITY框中输入0.285;Main Menu&Preprocessor&Material Props&-Temp Dependent-&Temp Table,输入温度70,200,300,400,500;Main Menu&Preprocessor&Material Props&-Temp Dependent-&Prop Table,选择导热系数KXX,材料编号为1,输入与温度表对应的导热系数8.35/12,8.9/12,9.35/12,9.8/12,10.23/12,选择APPLY;选择比热C,材料编号为1,输入0.113,0.117,0.119,0.122,0.125,选择APPLY;选择对流系数HF,材料编号为2,输入426/144,405/144,352/144,275/144,221/144,选择OK。定义几何模型参数:Utility Menu&Parameters&Scalar Parameters,输入ri1=1.3,ro1=1.5,z1=2,ri2=0.4,ro2=0.5,z2=2;建立几何模型Main Menu&Preprocessor&-Modeling-&Create&-Volumes-&Cylinder&ByDimensions,Outer radius框中输入ro1,Optional inner radium框中输入ri1,Z coordinates框中输入0和Z1,Ending angle框中输入90;Utility Menu&WorkPlane&Offset WP by Increments,在XY,YZ,ZX框中输入0,-90;Main Menu&Preprocessor&-Modeling-&Create&-Volumes-&Cylinder&ByD Outer radius框中输入ro2,Optional inner radium框中输入ri2,Z coordinates框中输入0和Z2,Starting angle框中输入-90,Ending angle框中输入0;Utility Menu&WorkPlane&Align WP with&Global Cartesian;进行布尔操作:Main Menu&Preprocessor&-Modeling-&Operate&-Booleans-&Overlap &Volumes,选择Pick All;观察几何模型Utility Menu&PlotCtrls&Numbering,打开Utility Menu&PlotCtrls&View Direction,在Coords of view point框中输入-3,-1,1;删除多余实体Main Menu&Preprocessor&-Modeling-&Delete&Volume and Below,在命令输入行输入3,4回车;创建组AREMOTEUtility Menu&Select&Entities,选择Area,By location,Z Coordinates,在Min,Max框中输入Z1,选择APPLY,Y Coordinates,在Min,Max框中输入0,OK;Utility Menu&Select&Comp/Assembly&Create Component,在Component name框中输入AREMOTE,在Components is made of菜单中选择AREA;组合面及线Main Menu&Preprocessor&-Meshing-&Mesh&-Volumes-&Mapped&-Concatenate-&Area,选择Pick all;Main Menu&Preprocessor&-Meshing-&Mesh&-Volumes-&Mapped&-Concatenate-&Lines,在命令行中输入12,7回车,选择APPLY,在命令行中输入10,5回车,OK;设定网格密度Main Menu&Preprocessor&-Meshing-&Size Cntrls&Picked Lines,选择线6和20,OK,在No,of element divisions框中输入4,OK;Main Menu&Preprocessor&-Meshing-&Size Cntrls&Picked Lines,选择线40,OK,在No,of element divisions框中输入6,OK;Utility Menu&Select&Everything;Main Menu&Preprocessor&-Meshing-&Size Cntrls&-Global-&Size,在element edge length框中输入0.4,OK;划分网格:Main Menu&Preprocessor&-Meshing-&Mesh&-Volumes-&Mapped&4 to 6 sides,选择Pick All;定义求解类型及选项Main Menu&Solution&-Analysis Type-&New Analysis,选择Steady-State;Main Menu&Solution&-Analysis Options,选择Program-chosen;施加对流载荷Utility Menu&WorkPlane&Change Active CS to&Global Cylindrical;Utility Menu&Select&Entities,选择Nodes,By location,X,在Min,Max框中输入ri1,OK;Main Menu&Solution&-Loads-&Apply&-Thermal-&Convection&On Nodes,选择Pick All,输入250/144及450,OK;在AREMOTE组上施加温度约束Utility Menu&Select&Comp/Assembly&Select Comp/Assembly,选Utility Menu&Select&Entities,选择Nodes,Attached to,On the Area all,OK;Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes,选择Pick all,输入45,OK;施加与温度有关的对流边界条件Utility Menu&WorkPlane&Offset WP by Increments,在XY,YZ,ZX Angles框中输入0,-90,OK;Utility Menu&WorkPlane&Local Coordinate Systems&Create Local CS&At WP Origin,在Type of coordinate system菜单中,选择Cylindrical 1,OK;Utility Menu&Select Entities,选择Nodes,By location,X,在Min,Max框中输入ri2,OK;Main Menu&Solution&-Loads-&Apply&-Thermal-&Convection&On Nodes,选择Pick All,在Film coefficient框中输入-2,在Bulk temperature框中输入100,OK;Utility Menu&Select&Everything;Utility Menu&PlotCtrls&Symbols,在Show pres and convect as菜单中选择Arrow,OK;Utility Menu&Plot&Nodes;恢复工作平面及坐标系统Utility Menu&WorkPlane&Change Active CS to&Global Cartesian;Utility Menu&WorkPlane&Align WP with&Global Cartesian;设定载荷步选项,Main Menu&Solution&-Load Step Options-&Time/Frequenc&Time and Substeps,在Number of substeps框中输入50,设置Automatic time stepping为On;求解:Main Menu&Solution&-Solve-&Current LS显示温度分布彩色云图,Main Menu&General Postproc&Plot Results&-Contour Plot-&Nodal Solu,选择Temperature TEMP。《ANSYS Verification Manual》中关于稳态热分析的实例:VM58
Centerline temperature of a heat generating wireVM92
Insulted wall temperatureVM93
Temperature dependent conductivityVM94
Heat generating plateVM95
Heat transfer from a cooling spineVM96
Temperature distribution in a short solid cylinderVM97
Temperature distribution along a straight finVM98
Temperature distribution along a tapered finVM99
Temperature distribution in a trapezoidal finVM100
Heat conductivity across a chimney sectionVM101
Temperature distribution in a short solid cylinderVM102
Cylinder with temperature dependent conductivityVM103
Thin plate with a central heat sourceVM105
Heat generation coil with temperature dependent conductivityVM108
Temperature gradient across a solid cylinderVM118
Centerline temperature of a heat generating wireVM160
Solid cylinder with harmonic temperature loadVM161
Heat flow from a insulated pipeVM162
Cooling of a circular fin of rectangular profileVM193
Adaptive analysis of two-dimensional heat transfer withconvection第四章 瞬态传热分析一、瞬态传热分析的定义瞬态热分析用于计算一个系统的随时间变化的温度场及其它热参数。在工程上一般用瞬态热分析计算温度场,并将之作为热载荷进行应力分析。瞬态热分析的基本步骤与稳态热分析类似。主要的区别是瞬态热分析中的载荷是随时间变化的。为了表达随时间变化的载荷,首先必须将载荷~时间曲线分为载荷步。载荷~时间曲线中的每一个拐点为一个载荷步,如下图所示。对于每一个载荷步,必须定义载荷值及时间值,同时必须选择载荷步为渐变或阶越。二、瞬态热分析中的单元及命令瞬态热分析中使用的单元与稳态热分析相同。要了解每个单元的详细说明,请参阅《ANSYS Element Reference Guide》。要了解每个命令的详细功能,请参阅《ANSYS Command Reference Guide》。三、ANSYS 瞬态热分析的主要步骤建模加载求解后处理四、建模确定jobname、title、units,进入PREP7;定义单元类型并设置选项;如果需要,定义单元实常数;定义材料热性能:一般瞬态热分析要定义导热系数、密度及比热;建立几何模型;对几何模型划分网格。关于建模及划分网格,请参阅《ANSYS Modeling and Meshing Guide》。五、加载求解1、定义分析类型如果第一次进行分析,或重新进行分析GUI,Main Menu&Solution&Analysis Type&New Analysis&TransientCommand,ANTYPE,TRANSIENT,NEW如果接着上次的分析继续进行(例如增加其它载荷)GUI,Main Menu&Solution&Analysis Type&RestartCommand,ANTYPE,TRANSIENT,REST2、获得瞬态热分析的初始条件①、定义均匀温度场如果已知模型的起始温度是均匀的,可设定所有节点初始温度Command: TUNIFGUI,Main Menu& Solution&-Loads-&Settings&Uniform Temp如果不在对话框中输入数据,则默认为参考温度,参考温度的值默认为零,但可通过如下方法设定参考温度:Command: TREFGUI,Main Menu& Solution&-Loads-&Settings&Reference Temp注意:设定均匀的初始温度,与如下的设定节点的温度(自由度)不同Command: DGUI,Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes初始均匀温度仅对分析的第一个子步有效;而设定节点温度将保持贯穿整个瞬态分析过程,除非通过下列方法删除此约束:Command: DDELEGUI,Main Menu& Solution&-Loads-&Delete&-Thermal-Temperature&On Nodes②、设定非均匀的初始温度在瞬态热分析中,节点温度可以设定为不同的值:Command: ICGUI,Main Menu& Solution&Loads&Apply&-Initial Condit'n&Define如果初始温度场是不均匀的且又是未知的,就必须首先作稳态热分析确定初始条件:设定载荷(如已知的温度、热对流等)将时间积分设置为OFF:Command,TIMINT,OFFGUI,Main Menu& Preprocessor&Loads&-Load Step Opts-Time/Frequenc&Time Integration设定一个只有一个子步的,时间很小的载荷步(例如0.001):Command,TIMEGUI,Main Menu& Preprocessor&Loads&-Load Step Opts-Time/Frequenc&Time and Substps写入载荷步文件:Command: LSWRITEGUI,Main Menu& Preprocessor&Loads&Write LS File 或先求解:Command: SOLVEGUI,Main Menu& Solution&Solve&Current LS注意:在第二载荷步中,要删去所有设定的温度,除非这些节点的温度在瞬态分析与稳态分析相同。3、设定载荷步选项①、普通选项时间:本选项设定每一载荷步结束时的时间:Command,TIMEGUI,Main Menu& Solution&-Load Step Opts-Time/Frequenc&Time and Substps每个载荷步的载荷子步数,或时间增量对于非线性分析,每个载荷步需要多个载荷子步。时间步长的大小关系到计算的精度。步长越小,计算精度越高,同时计算的时间越长。根据线性传导热传递,可以按如下公式估计初始时间步长:其中为沿热流方向热梯度最大处的单元的长度,为导温系数,它等于导热系数除以密度与比热的乘积()。Command,NSUBST or DELTIMGUI,Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time and Substps如果载荷在这个载荷步是恒定的,需要设为阶越选项;如果载荷值随时间线性变化,则要设定为渐变选项:Command,KBCGUI,Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time and Substps②、非线性选项迭代次数:每个子步默认的次数为25,这对大多数非线性热分析已经足够。Command,NEQITGUI,Main Menu& Solution&-Load step opts&Nonlinear&Equilibrium Iter自动时间步长:本选项为ON时,在求解过程中将自动调整时间步长。Command,AUTOTSGUI,Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time and Substps时间积分效果:如果将此选项设定为OFF,将进行稳态热分析。Command,TIMINTGUI,Main Menu& Solution&-Load Step Opts-&Time/Frequenc&Time Integration③、输出选项控制打印输出:本选项可将任何结果数据输出到*.out 文件中Command,OUTPRGUI,Main Menu&Solution&-Load Step Opts-&Output Ctrls&Solu Printout控制结果文件:控制*.rth的内容Command,OUTRESGUI,Main Menu&Solution&-Load Step Opts-&Output Ctrls&DB/Results File(、存盘求解六、后处理ANSYS提供两种后处理方式:POST1,可以对整个模型在某一载荷步(时间点)的结果进行后处理;Command,POST1GUI,Main Menu&General Postproc.POST26,可以对模型中特定点在所有载荷步(整个瞬态过程)的结果进行后处理。Command,POST26GUI,Main Menu&TimeHist Postproc1、用POST1进行后处理进入POST1后,可以读出某一时间点的结果:Command,SETGUI,Main Menu&General Postproc&Read Results&By Time/Freq如果设定的时间点不在任何一个子步的时间点上,ANSYS会进行线性插值。此外还可以读出某一载荷步的结果:GUI,Main Menu&General Postproc&Read Results&By Load Step然后就可以采用与稳态热分析类似的方法,对结果进行彩色云图显示、矢量图显示、打印列表等后处理。2、用POST26进行后处理首先要定义变量:Command,NSOL or ESOL or RFORCEGUI,Main Menu&TimeHist Postproc&Define Variables然后就可以绘制这些变量随时间变化的曲线:Command,PLVARGUI,Main Menu&TimeHist Postproc&Graph Variables或列表输出:Command,PRVARGUI,Main Menu&TimeHist Postproc&List Variables此外,POST26还提供许多其它功能,如对变量进行数学操作等,请参阅《ANSYS Basic Analysis Procedures Guide》七、相变问题ANSYS热分析最强大的功能之一就是可以分析相变问题,例如凝固或熔化等。含有相变问题的热分析是一个非线性的瞬态的问题:相变问题需要考虑熔融潜热,即在相变过程吸收或释放的热量。ANSYS通过定义材料的焓随温度变化来考虑熔融潜热(如图所示)。焓的单位是J/m3,是密度与比热的乘积对温度的积分:求解相变问题,应当设定足够小的时间步长,并将自动时间步长设置为ON;选用低阶的热单元,例如PLANE55或SOLID70。如果必须选用高阶单元,请将单元选项KEYOPT(1)设置为1:Command,keyopt(1)=1GUI,Main Menu&Prepocessor&Element Type&Add/Edit/Delete&Options&-Specific heat matrix-&Diagonalized在设定瞬态积分参数时,请将THETA值设置为1(默认为0.5):Command,TINTPGUI,Main Menu& Solution&-Load and Step Opts-&Time/Frequence&Time intergration&THETA线性搜索将有助于加速相变问题的求解。Command,LNSRCHGUI,Main Menu&Solution&-Load and Step Opts-&Nonlinear&Line Search实例1:一个30公斤重、温度为70℃的铜块,以及一个20公斤重、温度为80℃的铁块,突然放入温度为20℃、盛满了300升水的、完全绝热的水箱中,如图所示。过了一个小时,求铜块与铁块的最高温度(假设忽略水的流动)。材料热物理性能如下:热性能单位制铜铁水导热系数W/m℃38370.61密度Kg/m388897833996比热J/kg℃3904484185以下列出log文件及菜单操作说明/filename,transient1
/title,Thermal Transient Exercise 1!进入前处理/prep7
et,1,plane77
! 定义单元类型mp,kxx,1,383
! 定义材料热性能参数mp,dens,1,8889
!1~铜,2~铁,3~水mp,c,1,390
mp,kxx,2,70mp,dens,2,7837mp,c,2,448mp,kxx,3,0.61mp,dens,3,996mp,c,3,4185rectnag,0,0.6,0,0.5
!创建几何实体rectang,0.15,0.225,0.225,0.27rectang,0.6-0.2-0.058,0.6-0.2,0.225,0.225+0.044aovlap,all
!布尔操作/pnum,area,1aplotaatt,1,1,1
!划分网格eshape,2esize,0.02amesh,2aatt,2,1,1amesh,3aatt,3,1,1eshape,3esize,0.05amesh,4/pnum,mat,1eplotfinish!加载求解/solu
antype,transtimint,off
!先作稳态分析,确定初始条件time,0.01
!设定只有一个子步的时间很小的载荷步deltim,0.01esel,s,mat,,3nsle,sd,all,temp,20esel,s,mat,,2nsle,sd,all,temp,80esel,s,mat,,1nsle,sd,all,temp,70allselsolve
!得到初始温度分布time,3600
!进行瞬态分析timint,on
!打开时间积分deltim,26,2,200
!设置时间步长,最大及最小时间步长autots,on
!打开自动时间步长ddelet,all,temp
!删除稳态分析中定义的节点温度
outres,all,1
!将每个子步的值写入数据库文件solvefinishsave!进入POST26后处理
solu,2,dtime,,dtime
!2~每一子步采用的时间步长nsol,3,node(0.5,0),temp,,T_Copper !3~铜块的中心点nsol,4,node(0.371,0.247,0),temp,,T_Iron
!4~铁块的中心点nsol,5,node(30,0,0),temp,,T_H2O_Bot
!5~水箱的底部 nsol,6,node(30,50,0),temp,,T_H2O_Top
!6~水箱的顶部 nsol,7,node(0,25,0),temp,,T_H2O_Left
!7~水箱的左部 nsol,8,node(60,25,0),temp,,T_H2O_Right
!8~水箱的右部
Plvar,2plvar,3,4,5,6,7,8finish!进入POST1后处理/post1
!设置为最后一个载荷子步set,lastesel,s,mat,,1nsle,splnsol,tempesel,s,mat,,2nsle,splnsol,tempfinish菜单操作:Utility Menu,File&Change Jobname,输入文件名Transient1;Utility Menu,File&Change Title,输入Thermal Transient Exercise 1;Main Menu&Preprocessor,进入前处理;Main Menu&Preprocessor&Element Type&Add/Edit/Delete,选择单元plane77;Main Menu&Preprocessor&Material Props&Isotropic,定义材料1(铜)的KXX等于 383、DENS等于8889、C等于1390;Main Menu&Preprocessor&Material Props&Isotropic,定义材料2(铁)的KXX等于 70、DENS等于7837、C等于448;Main Menu&Preprocessor&Material Props&Isotropic,定义材料3(水)的KXX等于0.61、DENS等于996、C等于4185;Main Menu&Preprocessor&-Modeling-&Create&-Areas-&Retangle&By Dimensions,输入X1=0,Y1=0,X2=0.6,Y2=0.5,选择Apply,输入X1=0.15,Y1=0.225,X2= 0.225,Y2=0.27,选择Apply,输入X1=0.6-0.2-0.058,Y1=0.225,X2=0.6-0.2,Y2=0.225+0.044,选择OK;Main Menu&Preprocessor&-Modeling-&Operate&Booleans&Overlap,选择Pick All;Utility Menu,Plotctrls&Numbering&Areas,on;Utility Menu:Plot&Areas;Main Menu&Preprocessor&-Attributes-&Define-&All Areas,选择材料1;Main Menu&Preprocessor&Meshing-&Size Cntrls-&-Manualsize-&-Global-&Size,输入单元大小0.02;Main Menu&Preprocessor&Meshing-&Mesh-&-Areas-&Mapped&3 or 4 sided,选择铜块;Main Menu&Preprocessor&-Attributes-&Define-&All Areas,选择材料2;Main Menu&Preprocessor&Meshing-&Mesh-&-Areas-&Mapped&3 or 4 sided,选择铁块;Main Menu&Preprocessor&-Attributes-&Define-&All Areas,选择材料3;Main Menu&Preprocessor&Meshing-&Size Cntrls-&-Manualsize-&-Global-&Size,输入单元大小0.05;Main Menu&Preprocessor&Meshing-&Mesh-&-Areas-&Free,选择水箱;Utility Menu&Plot&Area;Main Menu&Solution,进入加载求解;Main Menu&Solution&-Analysis Type-&New Analysis,选择Transient,定义为瞬态分析;Main Menu&Solution&-Load Step Opts&Time/Frenquenc&Time Integration,TIMINT,off,首先进行稳态分析;Main Menu&Solution&-Load Step Opts&Time/Frenquenc&Time-Time Step,设定TIME为0.01、DELTIM也为0.01;Utility Menu,Select&Element&mat,输入3,选择Apply,选择Nodes&Attached to,选择Element;Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes,选择Pick All,输入20;Utility Menu,Select&Element&mat,输入2,选择Apply,选择Nodes&Attached to,选择Element;Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes,选择Pick All,输入80;Utility Menu,Select&Element&mat,输入1,选择Apply,选择Nodes&Attached to,选择Element;Main Menu&Solution&-Loads-&Apply&-Thermal-&Temperature&On Nodes,选择Pick All,输入70;Utility Menu,Select EverythingMain Menu&Solution&-Solve-&Current LSMain Menu&Solution&-Load Step Opts&Time/Frenquenc&Time-Time Step,设定TIME=3600,DELTIM=26,最小、最大时间步长分别为2,200,将Autots设置为ON;Main Menu&Solution&-Load Step Opts&Time/Frenquenc&Time Integration,将TIMINT设置为ON;Main Menu&Solution&-Loads-&Delete&-Thermal-&Temperature&On Nodes,选择Pick All,删除稳态分析定义的节点温度;Main Menu&Solution&-Load Step Opts&Output Ctrls-&DB/Results,选择Every Substeps;Main Menu&Solution&-Solve-&Current LS;Main Menu&TimeHist PostPro,进入POST26;Main Menu&TimeHist PostPro&Define Variables&Add&Solution summary,在User specified label框中输入dtime,选择Solution Items&Step Time,选择OK,Add&Nodal result&,在Node number框中输入node(0.5,0),在User specified label框中输入T_Copper,同样输入其它节点;Main Menu&TimeHist PostPro&Graph Virables,输入变量代号,显示各变量随时间变化的曲线;Main Menu&General Postproc,进入POST1;Main Menu&General Postproc&-Read Results-&Last setUtility Menu&Select&Element&mat,输入1,选择Apply,Nodes&Attach to,选择Element,OK;Main Menu&General Postproc&Plot result&Nodal Solution,选择Utility Menu&Select&Element&mat,输入2,选择Apply,Nodes&Attach to,选择Element,OK;Main Menu&General Postproc&Plot result&Nodal Solution,选择实例2:一钢铸件及其砂模的横截面尺寸如图所示:砂模的热物理性能如下表所示:单位制导热系数(KXX)Btu/hr.in.oF0.025密度(DENS)lbm/in30.254比热(C)Btu/lbm.oF0.28铸钢的热物理性能如下表所示:单位制0oF2643oF2750oF2875oF导热系数Btu/hr.in.oF1.441.541.221.22焓Btu/in30128.1163.8174.2初始条件:铸钢的温度为2875oF,砂模的温度为80oF;砂模外边界的对流边界条件:对流系数0.014Btu/hr.in2.oF,空气温度80oF;求3个小时后铸钢及砂模的温度分布。/Title,Casting Solidification!进入前处理/prep7et,1,plane55
!定义单元mp,dens,1,0.054
!定义砂模热性能mp,kxx,1,0.025mp,c,1,0.28mptemp,1,0,75
!定义铸钢的热性能mpdata,kxx,2,1.44,1.54,1.22,1.22mpdata,enth,2,0,128.1,163.8,174.2mpplot,kxx,2mpplot,enth,2save!创建几何模型k,1,0,0,0
k,2,22,0,0k,3,10,12,0k,4,0,12,0/pnum,kp,1/pnum,line,1/pnum,area,1/Triad,ltopkplota,1,2,3,4
saverectng,4,22,4,8aplotaovlap,alladele,3aplotsave!划分网格esize,1amesh,5mat,2aplotamesh,4eplot/pnum,elem/number,1save!进入加载求解/SOLUantype,trans
!设定为瞬态分析esel,s,mat,,2
!设定铸钢的初始温度nsle,s/replotic,all,temp,2875esel,inve
!设定砂模的初始温度nsle,s/replotic,all,temp,80allselsavelplotsfl,1,CONV,0.014,,80
!设定砂模外边界对流sfl,3,CONV,0.014,,80sfl,4,CONV,0.014,,80/psf,conv,2time,3
!设定瞬态分析时间kbc,1
!设定为阶越的载荷autots,on
!打开自动时间步长deltim,0.01,0.001,0.25
!设定时间步长timint,on
!打开时间积分tintp,,,,1
!将THETA设定为1outres,all,all
!输入每个子步的结果solve!进入后处理/post26/pnum,node,1/number,0eplotnsol,2,204,temp,center
!设定铸钢中心点温度随时间的变量plvar,2
!绘制温度~时间曲线savefinish菜单操作:Utility Menu&File&Change Title,输入Casting S定义单元类型:Main Menu&Preprocessor&Element Type&Add/Edit/Delete,Add,Quad 4node 55;定义砂模热性能:Main Menu&Preprocessor&Material Props&Isotropic,默认材料编号1,在Density(DENS)框中输入0.054,在Thermal conductivity (KXX)框中输入0.025,在Specific heat(C)框中输入0.28;定义铸钢热性能温度表:Main Menu&Preprocessor&Material Props&-Temp Dependent-&Temp Table,输入T1=0,T2=0,T4=2875;定义铸钢热性能:Main Menu&Preprocessor&Material Props&-Temp Dependent -&Prop Table,选择Th Conductivity,选择KXX,输入材料编号2,输入C1=1.44,C2=1.54,C3=1.22,C4=1.22,选择Apply,选择Enthalpy,输入C1=0,C2=128.1,C3=163.8,C4=174.2;创建关键点:Main Menu&Preprocessor&-Modeling-&Create&Keypoints&In Active CS,输入关键点编号1,输入坐标0,0,0,输入关键点编号2,输入坐标22,0,0,输入关键点编号3,输入坐标10,12,0,输入关键点编号4,输入坐标0,12,0;创建几何模型:Main Menu&Preprocessor&-Modeling-&Create&-Areas-&Arbitrary&Through KPs,顺序选取关键点1,2,3,4;Main Menu&Preprocessor&-Modeling-&Create&-Areas-&Rectangle&By Dimension,输入X1=4,X2=22,Y1=4,Y2=8;进行布尔操作:Main Menu&Preprocessor&-Modeling-&Operate&-Booleans-&Overlap&Area,Pick all;删除多余面:Main Menu&Preprocessor&-Modeling-&Delete&Area and Below,3保存数据库:在Ansys Toolbar中选取SAVE_DB;定义单元大小:Main Menu&Preprocessor&-Meshing-&Size Cntrls&-Global-&Size,在Element edge length框中输入1;对砂模划分网格:Main Menu&Preprocessor&-Meshing-&Mesh&-Areas-&Free,选择砂模;对铸钢划分网格:Main Menu&Preprocessor&-Attributes-&Define&Default Attribs,在Material number菜单中选择2;Main Menu&Preprocessor&-Meshing-&Mesh&-Areas-&Free,选择铸钢;定义分析类型:Main Menu&Solution&-Analysis Type-&New Analysis,选择T选择铸钢上的节点:Utility Menu&Select&Entities,选择element,mat,输入2,选择Apply,选择node,attached to element,选择OK;定义铸钢的初始温度:Main Menu&Solution&-Loads-&Apply&Initial Condit’n&Define,选择Pick all,选择temp,输入2875,OK;选择砂模上的节点:Utility Menu&Select&Entities,Nodes,inverse定义砂模的初始温度:Main Menu&Solution&-Loads-&Apply&Initial Condit’n&Define,选择Pick all,选择temp,输入80,OK;Utility Menu&Select&Everything;Utility Menu&Plot&Lines;定义对流边界条件: Main Menu&Solution&-Loads-&Apply&-Thermal-&Converction&On Lines,选择砂模的三个边界1,3,4,在file coefficent框中输入80,在Bulk temperature框中输入,80;设定瞬态分析时间选项:Main Menu&Solution&Load Step Opts&Time/Frequenc&Time-Time Step,Time at end of load step 3Time Step size
0.01Stepped or ramped b.c. SteppedAutomatic time stepping onMinimun time Step size 0.001Maximum time step size 0.25设置输出:Main Menu&Solution&Load Step Opts&Output Ctrls&DB/Results File,在File write frequency框中选择E求解:Main Menu&Solution&-Solve-&Current LS;进入后处理,Main Menu&Timehist Postproc;定义铸钢中心节点的温度变量: Main Menu&Timehist Postproc&Define Variables,Add,Nodal DOF result,2,204;绘制节点温度随时间变化曲线:Main Menu&Timehist Postproc&Graph Variable,2。《ANSYS Verification Manual》中关于瞬态热分析的实例:VM104
Liquid-solid phase changeVM109
Temperature gradient across a solid cylinderVM110
Transient temperature distribution in a slabVM111
Cooling of a spherical bodyVM112
Cooling of a spherical bodyVM113
Transient temperature distribution in an orthotropic metal barVM114
Temperature response to a linearly rising surface temperatureVM115
Thermal response of heat generating slabVM116
Heat conducting plate with sudden cooling第五章 热辐射一、什么是热辐射辐射是一种通过电磁波传递能量的方式。电磁波以光速传播且无需任何介质。热辐射仅为电磁波谱中的一小段。因为由于热辐射引起的热流与物体表面绝对温度的四次方成正比,因此热辐射分析是高度非线性的。二、分析热辐射问题ANSYS提供了三 种方法分析热辐射问题:用LINK31,辐射线单元,分析两个点或多对点之间的热辐射;用表面效应单元SURF19或SURF22,分析点对面的热辐射;用AUX12,热辐射矩阵生成器,分析面与面之间的热辐射以上三种方法既可用于稳态热分析,也可用于瞬态热分析。热辐射分析要注意温度的单位制,因为计算热辐射使用的温度单位是绝对温度。如果在加载时使用的是华氏温度,就要设置460的差值;如果为摄氏温度,差值为273。 Command:TOFFST GUI,Main Menu&Preprocessor&Loads&Analysis Options GUI,Main Menu&Solution&Analysis Options三、使用LINK31―辐射线单元LINK31是一个两节点非线性线单元,用于计算由辐射引起的两点之间的热传递。此单元要求输入如下的实常数:有效的热辐射面积;形状系数辐射率Stefan-Boltzmann 常数四、使用表面效应单元表面效应单元可以方便地分析点与面之间的辐射传热。SURF19用于两维模型,SURF22用于三维模型。单元应设置为包含辐射KEYOPT(9)。五、使用AUX12―辐射矩阵生成器此方法用于计算多个辐射面之间的辐射传热。这种方法生成辐射面之间形状系数矩阵,并将此矩阵作为超单元用于热分析。AUX12方法由三个步骤组成:定义辐射面生成辐射矩阵在热分析中使用辐射矩阵定义辐射面在PREP7中建模、划分网格。辐射面往往是3D模型中的面或2D模型中的边,如下图所示:在辐射表面用SHELL57(3D)或LINK32(2D)划分网格。最好的方法是先选择辐射表面的节点,然后用如下方法创建SHELL57或LINK32单元:Command,ESURFGUI,Main Menu&Preprocessor&Create&Elements&On Free Surf注意:辐射面上的SHELL57或LINK32单元与节点必须与实体单元相吻合,否则计算的结果是不正确的。生成的SHELL57或LINK32单元的取向也很重要。AUX12假定辐射的方向是SHELL57的+Z向或LINK32的+Y向。因此在生成SHELL57或LJNK32单元时要注意节点的排列顺序。如图所示:(3)、如果所分析的系统是开放的,即一个面所辐射的热能未被模型中其它的面吸收,则必须定义一个空间节点,用于吸收损失的辐射热量。这个节点的位置是任意的。对于封闭的系统,不应定义空间节点。2、生成辐射矩阵(1)、进入AUX12:Command,/AUX12GUI,Main Menu&Radiation Matrix(2)、选择组成辐射面的节点和单元。比较方便的方法是根据单元类型选择单元,并选择单元上的节点。有关选择目标的技巧,请参阅《ANSYS Basic Analysis Guide》(3)、确定模型是3D还是2D:Command,GEOMGUI,Main Menu&Radiation Matrix&Other SettingAUX12用不同的算法计算2D或3D模型的形状系数。AUX12默认为3D。2D分为纯平面或轴对称,默认为纯平面。(4)定义每个辐射面的辐射率(默认为1):Command,EMISGUI,Main Menu&Radiation Matrix&Emissivities(5)、定义Stefan-Boltzmann常数(默认为英制单位0.199E-10Btu/hr-in2-R4)Command,STEFGUI,Main Menu&Radiation Matrix&Other Settings(6)、确定用什么方式计算形状系数:Command,VTYPEGUI,Main Menu&Radiation Matrix&Write Matrix选择是隐藏还是非隐藏方法:非隐藏方法计算每两个单元之间的形状系数,无论它们之间有无障碍;隐藏方法(默认)用一种隐藏线算法判断两辐射面之间是否“可见”,如果可见则计算形状系数。(7)如为开放系统,定义空间节点:Command,SPACEGUI,Main Menu&Radiation Matrix&Other Settings(8)计算辐射矩阵并写入jobename.sub文件:Command,WRITEGUI,Main Menu&Radiation Matrix&Write Matrix如果要打印此矩阵,在执行上述命令之前输入MPRINT,1。(9)选择所有的节点和单元3、在热分析中使用辐射矩阵(1)、重新进入PREP7,定义一个新的单元类型MATRIX50(超单元);(2)、将单元类型指向超单元:Command,TYPEGUI,Main Menu&Preprocessor&-Modeling-Create&Elements&Elem Attributes(3)、读入超单元矩阵:Command,SEGUI,Main Menu&Preprocessor&-Modeling-Create&Elements&-Superelements-From,SUB File(4)、不选择或删除用于生成辐射矩阵的SHELL57或LINK32单元:Command,EDELEGUI,Main Menu&Preprocessor&-Modeling-Delete&Elements(5)、进入/SOLUTION施加其它热载荷并求解。目 录第一章 简介……………………………………………………………………….1 一、热分析的目的…………………………………………………………1 二、ANSYS的热分析………………………………………………………1 三、ANSYS热分析分类……………………………………………………1耦合分析…………………………………………………………….1第二章 基础知识…………………………………………………………………2符号与单位………………………………………………………….2传热学经典理论回顾………………………………………………2热传递的方式………………………………………………………3稳态传热……………………………………………………………3瞬态传热……………………………………………………………4线性与非线性………………………………………………………4边界条件、初始条件…………………………………………………4热分析误差估计……………………………………………………4第三章 稳态传热分析……………………………………………………………5稳态传热的定义……………………………………………………5热分析的单元………………………………………………………5ANSYS稳态热分析的基本过程……………………………………5实例1………………………………………………………………………9实例2………………………………………………………………………12第四章 瞬态传热分析…………………………………………………………20瞬态传热分析的定义………………………………………………20瞬态热分析的单元及命令…………………………………………20ANSYS瞬态热分析的主要步骤……………………………………20建模…………………………………………………………………20加载求解……………………………………………………………21后处理………………………………………………………………23相变问题…………………………………………………………..23实例1……………………………………………………………………24实例2…………………………………………………………………….25
课件名称:课件分类:机械课件类型:教学课件文件大小:2.04MB下载次数:39评论次数:11用户评分:6.8
4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22.}