《ANSYS热辐射学习教程.pptx》由会员分享，可在线阅读，更多相关《ANSYS热辐射学习教程.pptx（83页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、1不象其他热传递方式需要介质，辐射在真空中(如外层空间)效率最高。对于半透明体(如玻璃),辐射是三维实体现象,因为辐射从体中发散出。对于不透明体，辐射主要是平面现象因为几乎所有内部辐射都被实体吸收了。ANSYS 可以模拟不透明体间的辐射，所以我们将讨论范围限制在平面辐射现象上。第1页/共83页2平面总发射率,是平面在所有方向使用所有波长发射热的能力。这是一个无量纲数值。平面在所有方向用所有波长发射的总能量(热流单位)由施蒂芬-波斯曼定律确定:Diffuse SurfaceSpectral Surfaceqq第2页/共83页3平面可以理想化为散射或反射装置。散射装置会将辐射均匀反射到所有方向，而
2、不管辐射源的方位:通常情况下，平面可以被理想化为散射或反射面。反射平面会将辐射以近乎镜象的方式反射:没有实际的平面是真正的散射或反射面。比较灰暗的平面接近散射面，高度抛光的平面接近反射面。第3页/共83页4为了简化计算，平面的辐射特性可以在所有的波长和方向平均。因此，在散射和反射平面之间没有差别。两平面间的辐射热传递与它们平面绝对温度差的四次方成正比:第4页/共83页51.黑体黑体是理想化的平面，用来与实际平面进行比较。黑体的特性:黑体吸收所有的偶然辐射(没有反射),不管波长和方向。黑体为纯粹的发射器。对于给定的波长和温度，没有平面比黑体发射更多的能量。黑体是纯粹的散射发射器;辐射在所有方向均
3、一致。因此，对于黑体:B=B=1第5页/共83页6.灰体实际平面叫做“灰体”因为他们不象黑体。3.ANSYS和辐射ANSYS中关于辐射的重要假设和方法:ANSYS 认为辐射是平面现象，因此适合用不透明平面建模。ANSYS 不直接计入平面反射率。考虑到效率，假设平面吸收率和发射率相等(a=)。因此，只有发射率特性需要在ANSYS辐射分析中定义。ANSYS 不自动计入发射率的方向特性，也不允许发射率定义随波长变化。发射率可以在某些单元中定义为温度的函数。第6页/共83页7到此为止我们只讨论了单独的辐射平面。但是，在研究实际问题时，我们通常要考虑多个辐射平面的相互作用。要考虑的平面越多，问题越复杂:
4、第7页/共83页84.形状因子“形状因子”由相互辐射的两个平面(i和j)定义。它的定义是由于从一个平面(i)发射的辐射能偶然施加到另一个平面(j)上而得到。两个平面的形状因子是面积，方向和距离的函数。Ai、Aj表示表面I与表面J的面积，表示面单元dAi与面单元dAj之间的距离，i表示面单元dAj的法线Nj与两面但与连线的夹角，Ni表示面单元dAi的法线；Nj表示面单元dAj的法线第8页/共83页9第9页/共83页10角系数Fij具有以下特征：与面距离的平方成反比；与cosi成正比，在i=0时表面I辐射出最大的热量。与cosj成正比，在j=0时表面I辐射出最大的热量。第10页/共83页11若两表
5、面彼此“看不见”（即cosi=0且cosj），则彼此的角系数等于0；根据相互作何用原理，从任何平面发射的能量必须守恒，对于任意两表面均有AiFij=Aj Fji热辐射计算方程：两个表面之间的热辐射计算公式为：第11页/共83页12上式中各参量的物理意义如下：Q：表面I的传热率 ：Stefan-Bolzman常数 ：有效热辐射率；：表面I的面积；、：表面I与表面J的绝对温度值。、第12页/共83页135.2 热辐射问题分析ANSYS提供了三种方法分析热辐射问题：（1）用LINK31，辐射线单元，分析两个点或多对点之间的热辐射（2）用表面效应单元SURF19 或SURF22，分析点对面的热辐射（3
6、）用AUX12，热辐射矩阵生成器，分析面与面之间的热辐射以上三种方法既可用于稳态热分析，也可用于瞬态传热分析。第13页/共83页141.点-点问题模拟两节点或多对节点的热辐射问题时，应使用热辐射线单元LINK31。LINK31作为两节点非线性单元可以计算两点之间因辐射引起的热交换。使用此单元需指明以下常数：有效的辐射表面面积角系数辐射率Stefan-bolzmann常数。第14页/共83页152.点-面问题模拟点面之间的热辐射问题，通常采用表面效应单元-SURF151和SURF152。其中SURF151用于2D单元和SURF152用于3D单元。表面效应单元利用实体表面的节点形成单元，并且直接覆
7、盖在实体单元的表面。使用表面效应单元可以更第15页/共83页16灵活的在实体表面施加热载荷。例如，热流密度与对流可以施加在同一外表面，但ANSYS在计算过程中仅读取最后施加的面载荷进行计算。为避免ANSYS只读取一种载荷，可以利用实体单元承受热流密度，而表面效应单元承受对流载荷。第16页/共83页17SURF151实常数:FORMF(角系数)、SBCONST（Stefan-Bolzmann常数）材料属性：DENS（密度）、EMIS（辐射率）表面载荷：对流、热流密度。体载荷：生热率第17页/共83页18SURF152实常数:FORMF(角系数)、SBCONST（Stefan-Bolzmann常数
8、）材料属性：DENS（密度）、EMIS（辐射率）表面载荷：对流、热流密度。体载荷：生热率第18页/共83页19注意：热辐射分析要注意温度的单位制，因为计算热辐射使用的温度单位是绝对温度。如果在加载时使用的是华氏温度，就要设置460 的差值；如果为摄氏温度，差值为273。第19页/共83页203.使用LINK31 单元的注意事项LINK31 是一个两节点非线性线单元，用于计算由辐射引起的两点之间的热传递。此单元要求键入如下的实常数：有效的热辐射面积形状系数辐射率Stefan-Boltzmann常数 第20页/共83页214.使用表面效应单元的注意事项表面效应单元可以方便地分析点与面之间的辐射传热
9、。SURF19 用于两维模型，SURF22 用于三维模型。单元应设置为包含辐射KEYOPT（9）。第21页/共83页225.使用AUX12辐射矩阵生成器此方法用于计算多个辐射面之间的辐射传热。这种方法生成辐射面之间形状系数矩阵，并将此矩阵作为超单元用于热分析。AUX12 方法由三个步骤组成，即首先定义辐射面，然后生成辐射矩阵，最后在热分析中使用辐射矩阵。下面介绍使用AUX12 求解热辐射问题的步骤：第22页/共83页23(1)定义辐射面首先建立有限元模型在辐射面上覆盖一层SHELL57辐射面往往是3D模型中的面或2D模型中的边，因此在辐射表面用SHELL57（3D）或LINK32（2D）划分网
10、格。最好的方法是先选择辐射表面的节点，然后用如下方法创建SHELL57 或LINK32 单元：命令：ESURFGUI：Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf-suf General Surface 第23页/共83页24注意：辐射面上的SHELL57 或LINK32 单元与节点必须与实体单元相吻合，否则计算的结果是不正确的。生成的SHELL57 或LINK32 单元的取向也很重要。AUX12 假定辐射的方向是SHELL57 的+Z 向或LINK32 的+Y 向。因此在生成SHELL57 或LJNK32
11、单元时要注意节点的排列顺序。如果所分析的系统是开放的，即一个面所辐射的热能未被模型中其它的面吸收，则必须定义一个空间节点，用于吸收损失的辐射热量。这个节点的位置是任意的。对于封闭的系统，不应定义空间节点。第24页/共83页25(2)生成辐射矩阵生成辐射矩阵首先需要进入AUX12，然后选择所需节点和单元，确定模型的维数。进入AUX12命令：/AUX12GUI：Main Menu Radiation Matrix选择组成辐射面的节点和单元比较方便的方法是根据单元类型选择单元，并选择单元上的节点。确定模型是3D 还是2D命令：GEOM第25页/共83页26GUI：Main Menu Radiatio
12、n Matrix Other SettingAUX12 用不同的算法计算2D 或3D 模型的形状系数。AUX12 默认为3D。2D 分为纯平面或轴对称，默认为纯平面。定义每个辐射面的辐射率命令：EMISGUI：Main Menu Radiation Matrix Emissivities辐射面的辐射率ANSYS 默认为1。定义Stefan-Boltzmann 常数命令：STEFGUI：Main Menu Radiation Matrix Other SettingsStefan-Boltzmann 常数第26页/共83页27(3)确定状系数命令：VTYPEGUI：Main Menu Radia
13、tion Matrix Write Matrix用什么方式计算形形状系数，ANSYS 提供两种选择，即选择是隐藏还是非隐藏方法。非隐藏方法计算每两个单元之间的形状系数，无论它们之间有无障碍；隐藏方法默认）用一种隐藏线算法判断两辐射面之间是否“可见”，如果可见则计算形状系 第27页/共83页28(4)定义空间节点：命令：SPACEGUI：Main Menu Radiation Matrix Other Settings如为开放系统.(5)计算辐射矩阵并写入jobename.sub 文件命令：WRITEGUI：Main Menu Radiation Matrix Write Matrix(6)选择
14、所有的节点和单元 第28页/共83页295.3 热辐射问题分析11.问题描述如图所示一黑体，表面积1m2，形状系数和辐射率为1，温度为2000，周围环境温度为0，求黑体的辐射热流率。黑体结构示意图第29页/共83页302.问题分析选择LINK31热辐射单元进行求解。3.建立模型(1)过虑菜单简化菜单（过虑菜单）操作如下：GUI：Main Menu Referenc在弹出的Referenc for GUI Filtering对话框中，选择 Thermal。单击OK。命令：/COM,Thermal第30页/共83页31(2)选择单元选择热分析单元，操作如下：GUI：Main Menu Prepro
15、cessor Element Type Add/Edit/Delete选择Thermal Link3D radiation 31选项，在Element type reference number 文本框中输入如1。第31页/共83页32GUI：Main Menu PreprocessorReal contantsAdd/edit/delete 命令，点击Add按钮，在Element type for real contants对话框，单击OK按钮，进行如图设置第32页/共83页33(3)建立模型首先进入Define Material Model Behavior对话框，操作如下：GUI：Mai
16、n Menu Preprocessor Modeling Createnodeson active cs 命令，在弹出的对话框中，在Node Numbers 文本框中输入1，在X，Y，Z文本框中依次输入0、0、0。点击Apply按钮，在Node Numbers 文本框中输入2，在X，Y，Z文本框中依次输入0、0、0。第33页/共83页34GUI：Main Menu Preprocessor Modeling CreateElementAuto NumberedThru Nodes，出现Element from nodes菜单，在文本框中输入1，2，单击OK。UnitityFileSave as
17、命令 进行存盘。第34页/共83页354.加载求解选择Transient分析，操作如下：GUI：Main Menu Preprocessor Loads Analysis Type New Analysis选择Steady-state 分析，单击OK。进入Time and time Step Optios 对话框，操作如下：GUI：Main Menu Solution Load Step Opts Time/Frequenc Time Time Step命令，进行如下图设置 第35页/共83页36第36页/共83页37择GUI：Main Menu Solution Radiation Opts
18、Solution Opt命令，在弹出的对话框里在STEF文本框中输入5.6E-8，在Toffset 文本框中输入273，其余选项采用默认设置，如下图所示，单击OK关闭对话框。GUI：Main Menu SolutionDefine loads applythermal temperatureon nodes命令，在弹出的对话框中，在文本框中输入1，在Lab2 DOFs to be constrained列表框中选择TEMP选项，在Value Load Temp Value文本框中输入2000，单击OK。第37页/共83页38GUI：Main Menu SolutionDefine loads
19、applythermal temperatureon nodes命令，在弹出的对话框中，在文本框中输入2，在Lab2 DOFs to be constrained列表框中选择TEMP选项，在Value Load Temp Value 文本框中输入0，单击OK。GUI：Main Menu Solutionsolve Current LS命令，单击OK，进行求解。第38页/共83页39第39页/共83页40查看求解结果GUI：Main Menu General PostprocRead resultlast setGUI：Main Menu General Postproclist results
20、Reaction Solu。在弹出的对话框中，在Lab Item to be listed 对话框中选择 Heat flow Heat选项。单击OK。第40页/共83页416.命令流文件/Filename,exercise1/title,radiant energy emissionKeyw PR_THERM,1/PREP7ET,1,LINK32R,1,1,1,1N,1N,2E,1,2ALLSELFINISH/SOLU第41页/共83页42ANTYPE,STATICTIME,1AUTOTS,ONKBC,1STEF,5.67e-8TOFFSET,273D,1,TEMP,2000D,2,TEMP,
21、0SOLVEFINISH/POST1SET,LASTPRRSOLFINISH/EXIT第42页/共83页435.4 热辐射问题分析21.问题描述两个等长度的同轴长圆柱体，截面如图所示，内外圆柱体的初始温度分别为1000、100，两者之间的辐射率为1，两圆柱体的材料性能参数如表所示，求两圆柱体间的辐射热流率。第43页/共83页44204080100第44页/共83页45导热系数K1导热系数K2密1密度2比热C1比热C2703078002700480350第45页/共83页462.问题分析该问题为平面热辐射问题，整个热辐射问题为封闭系统，不需设置空间节点，在计算过程中使用AUX12辐射矩阵生成器，
22、选择PLANE55 热分析单元，进行求解。第46页/共83页473.建立模型（1）过滤菜单简化菜单（过滤菜单）操作如下：GUI：Main Menu Referenc在弹出的Referenc for GUI Filtering对话框中，选择 Thermal。单击OK。命令：/COM,Thermal第47页/共83页48(2)选择单元选择热分析单元，操作如下：GUI：Main Menu Preprocessor Element Type Add/Edit/Delete选择Thermal SolidQuad 4node 55选项，在Element type reference number 文本框中
23、输入如1。点击OK 第48页/共83页49(3)定义材料性能参数GUI：Main Menu Preprocessor Material Props Material Models命令，出现Define Material Model Behavior 对话框。在弹出的定义材料热传导率对话框中的KXX 栏键入“70”。命令：MPDATA,KXX,1,70 第49页/共83页50(4)定义密度GUI：Main Menu Preprocessor Material Props Thermal Density在弹出密度定义对话框中的DENS栏键入“7800”。命令：MPDATA,DENS,1,7800第
24、50页/共83页51(5)定义比热容GUI：Main Menu Preprocessor Material Props Thermal Specific Heat在弹出的定义比热容对话框中的C栏键入“480”。命令：MPDATA,C,1,480 第51页/共83页52(6)定义材料2的材料特性在Define Material Model Behavior 对话框上单击Material New Model 命令Define Material ID 对话框，在文本框中输入2，单击OK。定义材料热传导率对话框中的KXX 栏键入“70”。命令：MPDATA,KXX,2,70第52页/共83页53(7)
25、定义密度GUI：Main Menu Preprocessor Material Props Thermal Density在弹出密度定义对话框中的DENS栏键入“1000”。命令：MPDATA,DENS,2,1000(8)定义比热容GUI：Main Menu Preprocessor Material Props Thermal Specific Heat在弹出的定义比热容对话框中的C栏键入“350”。命令：MPDATA,C,2,350第53页/共83页544.建立实体模型(1)几何模型GUI：Preprocessor Modeling Create Areas Circle Partial
26、Annulus在弹出的对话框中进行如下设置：WPX=0,WPY=0,Rad-1=0.1,Theta-1=0,Rad-2=0.2,Theta-2=90。单击Apply，进行如下设置：WPX=0,WPY=0,Rad-1=0.4,Theta-1=0,Rad-2=0.5,Theta-2=90。第54页/共83页55(2)网格划分在设定网格尺寸参数，操作如下：GUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked Lines在 Element Size On 菜单中选择线2、4在Element Size On Picke
27、d Lines 对话框中，在NDIV 文本框中输入“5”，单击OK第55页/共83页56GUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked Lines在 Element Size On 菜单中选择线1、3，在Element Size On Picked Lines 对话框中，在NDIV 文本框中输入“10”，单击OKGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked Lines第56页/共83页57在 Element S
28、ize On 菜单中选择线6、8，在Element Size On Picked Lines 对话框中，在NDIV 文本框中输入“6”，单击OKGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked Lines在 Element Size On 菜单中选择线5、7，在Element Size On Picked Lines 对话框中，在NDIV 文本框中输入“30”，单击OKGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked
29、 Lines第57页/共83页58(3)网格划分GUI：Main Menu Preprocessor MeshingMesh Tool，出现Mesh Tool菜单，在shape选项组中选中Quad和Mapped选项，单击Mesh，出现Mesh Areas，在文本框中输入1，单击OK。GUI：Utility Menu Select EverythingGUI：Main Menu Preprocessor MeshingMesh attributes Default attribs，出现Meshing attributes 菜单，在MAT下拉列表框中选择2，单击OK.第58页/共83页59GUI：
30、Main Menu Preprocessor MeshingMesh Tool，出现Mesh Tool菜单，单击Mesh，出现Mesh Areas，在文本框中输入2，单击OK。GUI：Main Menu Preprocessor ModelingReflectAreas 命令，出现Reflect Areas菜单，单击Pick All 按钮，在弹出的对话框中，选择Y-Z Plane X 单选按钮，在NOELEM下拉列表框中选择Areas and mesh选项，在IMOVE下拉列表框中选择Copied选项，单击OK。GUI：Main Menu Preprocessor ModelingReflec
31、tAreas 命令，出现Reflect Areas菜单，单击Pick All 按钮，在弹出的对话框中，选择X-Z Plane Y 单选按钮，其余采用默认值，单击OK。注意:以上两步间不能采用Apply 按钮过渡第59页/共83页60GUI：Main Menu Numbering Ctrls Merge Items 命令，出现Merge Coincident or Equivalently Defined Items 对话框，在Label Type of Item to be merge 在下拉列表框中选择Keypoints选项，单击OK。选择UTILITY MenuSelect everyth
32、ying 命令。第60页/共83页61生成AUX12辐射矩阵 选择Main Menu Radiation Opt Radiosity MethSolution Opt命令，出现Radiation Solution Options对话框，在STEF 文本框中输入5.67E-8，在TOFFSET文本框中输入0.001，在SPCTEMP/SPCNOD下拉列表框中选择temperature选择，在Value文本框中输入20，其余采用默认值，单击OK。如图所示第61页/共83页62第62页/共83页63加载求解（1）选择Transient分析，操作如下：GUI：Main Menu Preprocesso
33、r Loads Analysis Type New Analysis选择Steady-state 分析，单击OK。GUI：Main Menu SolutionDefine loads applythermalradiationon lines命令，出现Apply RDSF on lines 菜单，在文本框中输入1、7、9、15、17、23、25、31，单击Ok按钮，出现Apply RDSF on lines 对话框，在VALUE Emissivity文本框中输入1，在VALUE2的文本框中输入1，单击OK。第63页/共83页64GUI：Main Menu SolutionDefine load
34、s applythermaltemperatureon lines，出现Apply TEMP on lines 菜单，在文本框中输入5、13、21、29，单击Ok按钮，出现Apply TEMP on lines对话框，在LAB2列表框中选择TEMP，在VALUE文本框中输入“100”，选择KEXPND选项，使其状态为“Yes”GUI：Main Menu SolutionDefine loads applythermaltemperatureon lines，出现Apply TEMP on lines 菜单，在文本框中输入3、11、19、27，单击Ok按钮，出现Apply TEMP on lin
35、es对话框，在LAB2列表框中选择TEMP，在VALUE文本框中输入“1000”，选择KEXPND选项，使其状态为“Yes”。第64页/共83页65第65页/共83页66选择UTILITY MenuSelect everythying 命令。GUI：Main Menu Solution Load Step Opts Time/Frequenc Time Time Step命令，进行如下图设置GUI：Main Menu Solutionsolve Current LS命令，单击OK，进行求解。第66页/共83页67第67页/共83页686.查看求解结果(1)GUI：Main Menu Gener
36、al PostprocRead resultlast setUtility Menu SelectEntities 命令，在对话框中选择Element选择，在第2个下拉列表框中选择By Attributes 选项，在第3个下拉列表框中选择Material num单选按钮，在Min,Max,Inc 文本框中输入1，在第5个下拉列表框中选择From Full单选按钮，如下图所示Utility Menu SelectEntities 命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attached to 选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择From
37、 Full，如下图所示第68页/共83页69第69页/共83页70GUI：Main Menu General PostprocPlot ResultsContour Plot SelectEntities 命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attached to 选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择From Full。UI：Main Menu General PostprocPlot ResultsContour Plot SelectEntities 命令，在对话框中选择Element选择，在第2个下拉列表框中选择By Att
38、ributes 选项，在第3个下拉列表框中选择Material num单选按钮，在Min,Max,Inc 文本框中输入1，在第5个下拉列表框中选择From Full单选按钮。Utility Menu SelectEntities 命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attached to 选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择From Full。第74页/共83页75GUI：Main Menu General PostprocList ResultsReaction Solu 命令，在弹出的对话框中保持默认值，单击OK，显示内圆柱体
39、上所有节点的辐射换热量，如图所示。第75页/共83页767.命令流/PREP7ET,1,PLANE55MP,DENS,1,7800MP,KXX,1,70MP,C,1,480MP,DENS,2,2700MP,KXX,2,30MP,C,2,350CYL4,0,0,0.1,90,0.2CYL4,0,0,0.4,90,0.5第76页/共83页77 APLOT/PNUM,LINE,1/PNUM,AREA,1LSEL,S,2,4,2LESIZE,ALL,5LSEL,S,1,3,2LESIZE,ALL,10LSEL,S,6,8,2LESIZE,ALL,6LSEL,S,5,7,2LESIZE,ALL,30MS
40、HAPE,0,2D第77页/共83页78 MSHKEY,1AMESH,1APLOTMAT,2AMESH,2ARSYM,X,ALLARSYM,Y,ALLNUMMRG,NODENUMMRG,KPEPLOTALLSEL第78页/共83页79/AUX12STEF,5.67E-8TOFFSET,273RADOPT,0.001SPCTEMP,1,20FINISH第79页/共83页80 /SOLUANTYPE,STATICLPLOTLSEL,S,1,7,6LSEL,A,9,15,6LSEL,A,17,23,6LSEL,A,25,31,6SFL,ALL,RDSF,1,1LSEL,S,5,29,8DL,ALL,TEMP,100,1LSEL,S,3,27,8DAL,ALL,TMEP,1000,1ALLSELDETIM,0.5第80页/共83页81SOLVEFINISH/POST1SET,LASTESEL,S,MAT,1NSLE,SPLNSOL,TEMPESEL,S,MAT,2NSLE,SPLNSOL,TEMPESEL,S,MAT,1 第81页/共83页82NSLE,SPRRSOLESEL,S,MAT,2NSLE,SPRRSOLALLSELFINISH第82页/共83页83感谢您的观赏！第83页/共83页