ANSYS热辐射第5节.pptx

QUST15 热辐射分析热辐射分析n5.1 热辐射分析简介n 辐射是一种通过电磁波传递能量的方式。电磁波以光速传播且无需任何介质。热辐射仅为电磁波谱中的一小段。因为由于热辐射引起的热流与物体表面绝对温度的四次方成正比，因此热辐射分析是高度非线性的。n 辐射热传递是通过电磁波传递热能的方法。热辐射的电磁波波长为0.1 到 100 mm。QUST2k不象其他热传递方式需要介质，辐射在真空中(如外层空间)效率最高。k对于半透明体(如玻璃),辐射是三维实体现象,因为辐射从体中发散出。k对于不透明体，辐射主要是平面现象因为几乎所有内部辐射都被实体吸收了。kANSYS 可以模拟不透明体间的辐射，所以我们将讨论范围限制在平面辐射现象上。QUST3n平面总发射率,是平面在所有方向使用所有波长发射热的能力。这是一个无量纲数值。n平面在所有方向用所有波长发射的总能量(热流单位)由施蒂芬-波斯曼定律确定:Diffuse SurfaceSpectral SurfaceqqQUST4n平面可以理想化为散射或反射装置。散射装置会将辐射均匀反射到所有方向，而不管辐射源的方位:n通常情况下，平面可以被理想化为散射或反射面。反射平面会将辐射以近乎镜象的方式反射:n没有实际的平面是真正的散射或反射面。比较灰暗的平面接近散射面，高度抛光的平面接近反射面。QUST5n为了简化计算，平面的辐射特性可以在所有的波长和方向平均因此，在散射和反射平面之间没有差别。n两平面间的辐射热传递与它们平面绝对温度差的四次方成正比:QUST6n1.黑体n黑体是理想化的平面，用来与实际平面进行比较。n黑体的特性:k黑体吸收所有的偶然辐射(没有反射),不管波长和方向。k黑体为纯粹的发射器。对于给定的波长和温度，没有平面比黑体发射更多的能量。k黑体是纯粹的散射发射器;辐射在所有方向均一致。k因此，对于黑体:B=B=1QUST7n.灰体k实际平面叫做“灰体”因为他们不象黑体。n3.ANSYS和辐射nANSYS中关于辐射的重要假设和方法:kANSYS 认为辐射是平面现象，因此适合用不透明平面建模。kANSYS 不直接计入平面反射率。考虑到效率，假设平面吸收率和发射率相等(a=)。因此，只有发射率特性需要在ANSYS辐射分析中定义。kANSYS 不自动计入发射率的方向特性，也不允许发射率定义随波长变化。发射率可以在某些单元中定义为温度的函数。QUST8到此为止我们只讨论了单独的辐射平面。但是，在研究实际问题时，我们通常要考虑多个辐射平面的相互作用。要考虑的平面越多，问题越复杂:QUST9n4.形状因子n“形状因子”由相互辐射的两个平面(i和j)定义。它的定义是由于从一个平面(i)发射的辐射能偶然施加到另一个平面(j)上而得到。n两个平面的形状因子是面积，方向和距离的函数。nAi、Aj表示表面I与表面J的面积，表示面单元dAi与面单元dAj之间的距离，i表示面单元dAj的法线Nj与两面但与连线的夹角，Ni表示面单元dAi的法线；Nj表示面单元dAj的法线QUST10QUST11k角系数Fij具有以下特征：k与面距离的平方成反比；k与cosi成正比，在i=0时表面I辐射出最大的热量。k与cosj成正比，在j=0时表面I辐射出最大的热量。QUST12k若两表面彼此“看不见”（即cosi=0且cosj），则彼此的角系数等于0；k根据相互作何用原理，从任何平面发射的能量必须守恒，对于任意两表面均有AiFij=Aj Fjik热辐射计算方程：k两个表面之间的热辐射计算公式为：QUST13n上式中各参量的物理意义如下：kQ：表面I的传热率k ：Stefan-Bolzman常数k ：有效热辐射率；k ：表面I的面积；、k ：表面I与表面J的绝对温度值。、QUST14n5.2 热辐射问题分析nANSYS提供了三种方法分析热辐射问题：n（1）用LINK31，辐射线单元，分析两个点或多对点之间的热辐射n（2）用表面效应单元SURF19 或SURF22，分析点对面的热辐射n（3）用AUX12，热辐射矩阵生成器，分析面与面之间的热辐射n以上三种方法既可用于稳态热分析，也可用于瞬态传热分析。QUST15n1.点-点问题n模拟两节点或多对节点的热辐射问题时，应使用热辐射线单元LINK31。LINK31作为两节点非线性单元可以计算两点之间因辐射引起的热交换。使用此单元需指明以下常数：n有效的辐射表面面积n角系数n辐射率nStefan-bolzmann常数。QUST16n2.点-面问题n模拟点面之间的热辐射问题，通常采用表面效应单元-SURF151和SURF152。其中SURF151用于2D单元和SURF152用于3D单元。n表面效应单元利用实体表面的节点形成单元，并且直接覆盖在实体单元的表面。使用表面效应单元可以更QUST17n灵活的在实体表面施加热载荷。例如，热流密度与对流可以施加在同一外表面，但ANSYS在计算过程中仅读取最后施加的面载荷进行计算。为避免ANSYS只读取一种载荷，可以利用实体单元承受热流密度，而表面效应单元承受对流载荷。QUST18nSURF151n实常数:FORMF(角系数)、SBCONST（Stefan-Bolzmann常数）n材料属性：DENS（密度）、EMIS（辐射率）n表面载荷：对流、热流密度。n体载荷：生热率QUST19nSURF152n实常数:FORMF(角系数)、SBCONST（Stefan-Bolzmann常数）n材料属性：DENS（密度）、EMIS（辐射率）n表面载荷：对流、热流密度。n体载荷：生热率QUST20n注意：n热辐射分析要注意温度的单位制，因为计算热辐射使用的温度单位是绝对温度。如果在加载时使用的是华氏温度，就要设置460 的差值；如果为摄氏温度，差值为273。QUST21n3.使用LINK31 单元的注意事项nLINK31 是一个两节点非线性线单元，用于计算由辐射引起的两点之间的热传递。此单元要求键入如下的实常数：n有效的热辐射面积n形状系数n辐射率nStefan-Boltzmann常数 QUST22n4.使用表面效应单元的注意事项n表面效应单元可以方便地分析点与面之间的辐射传热。SURF19 用于两维模型，SURF22 用于三维模型。单元应设置为包含辐射KEYOPT（9）。QUST23n5.使用AUX12辐射矩阵生成器n此方法用于计算多个辐射面之间的辐射传热。这种方法生成辐射面之间形状系数矩阵，并将此矩阵作为超单元用于热分析。nAUX12 方法由三个步骤组成，即首先定义辐射面，然后生成辐射矩阵，最后在热分析中使用辐射矩阵。下面介绍使用AUX12 求解热辐射问题的步骤：QUST24n(1)定义辐射面n首先建立有限元模型n在辐射面上覆盖一层SHELL57辐射面往往是3D模型中的面或2D模型中的边，因此在辐射表面用SHELL57（3D）或LINK32（2D）划分网格。最好的方法是先选择辐射表面的节点，然后用如下方法创建SHELL57 或LINK32 单元：n命令：ESURFnGUI：Main Menu Preprocessor Modeling Create Elements Surf/Contact Surf-suf General Surface QUST25n注意：n辐射面上的SHELL57 或LINK32 单元与节点必须与实体单元相吻合，否则计算的结果是不正确的。n生成的SHELL57 或LINK32 单元的取向也很重要。AUX12 假定辐射的方向是SHELL57 的+Z 向或LINK32 的+Y 向。因此在生成SHELL57 或LJNK32 单元时要注意节点的排列顺序。n如果所分析的系统是开放的，即一个面所辐射的热能未被模型中其它的面吸收，则必须定义一个空间节点，用于吸收损失的辐射热量。这个节点的位置是任意的。对于封闭的系统，不应定义空间节点。QUST26n(2)生成辐射矩阵n生成辐射矩阵首先需要进入AUX12，然后选择所需节点和单元，确定模型的维数。n进入AUX12n命令：/AUX12nGUI：Main Menu Radiation Matrixn选择组成辐射面的节点和单元n比较方便的方法是根据单元类型选择单元，并选择单元上的节点。n确定模型是3D 还是2Dn命令：GEOMQUST27nGUI：Main Menu Radiation Matrix Other SettingnAUX12 用不同的算法计算2D 或3D 模型的形状系数。AUX12 默认为3D。2D 分为纯平面或轴对称，默认为纯平面。n定义每个辐射面的辐射率n命令：EMISnGUI：Main Menu Radiation Matrix Emissivitiesn辐射面的辐射率ANSYS 默认为1。n定义Stefan-Boltzmann 常数n命令：STEFnGUI：Main Menu Radiation Matrix Other SettingsnStefan-Boltzmann 常数QUST28n(3)确定状系数n命令：VTYPEnGUI：Main Menu Radiation Matrix Write Matrixn用什么方式计算形形状系数，ANSYS 提供两种选择，即选择是隐藏还是非隐藏方法。非隐藏方法计算每两个单元之间的形状系数，无论它们之间有无障碍；隐藏方法默认）用一种隐藏线算法判断两辐射面之间是否“可见”，如果可见则计算形状系 QUST29n(4)定义空间节点：n命令：SPACEnGUI：Main Menu Radiation Matrix Other Settingsn如为开放系统.n(5)计算辐射矩阵并写入jobename.sub 文件n命令：WRITEnGUI：Main Menu Radiation Matrix Write Matrixn(6)选择所有的节点和单元 QUST305.3 热辐射问题分析1n1.问题描述n如图所示一黑体，表面积1m2，形状系数和辐射率为1，温度为2000，周围环境温度为0，求黑体的辐射热流率。nn黑体结构示意图QUST31n2.问题分析k选择LINK31热辐射单元进行求解。n3.建立模型(1)过虑菜单k简化菜单（过虑菜单）操作如下：kGUI：Main Menu Referenck在弹出的Referenc for GUI Filtering对话框中，选择 Thermal。单击OK。k命令：/COM,ThermalQUST32n(2)选择单元k选择热分析单元，操作如下：kGUI：Main Menu Preprocessor Element Type Add/Edit/Deletek选择Thermal Link3D radiation 31选项，在Element type reference number 文本框中输入如1。QUST33kGUI：Main Menu PreprocessorReal contantsAdd/edit/delete 命令，点击Add按钮，在Element type for real contants对话框，单击OK按钮，进行如图设置QUST34n(3)建立模型n首先进入Define Material Model Behavior对话框，操作如下：nGUI：Main Menu Preprocessor Modeling Createnodeson active cs 命令，在弹出的对话框中，在Node Numbers 文本框中输入1，在X，Y，Z文本框中依次输入0、0、0。n点击Apply按钮，在Node Numbers 文本框中输入2，在X，Y，Z文本框中依次输入0、0、0。QUST35nGUI：Main Menu Preprocessor Modeling CreateElementAuto NumberedThru Nodes，出现Element from nodes菜单，在文本框中输入1，2，单击OK。nUnitityFileSave as命令 进行存盘。QUST364.加载求解n选择Transient分析，操作如下：nGUI：Main Menu Preprocessor Loads Analysis Type New Analysisn选择Steady-state 分析，单击OK。n进入Time and time Step Optios 对话框，操作如下：nGUI：Main Menu Solution Load Step Opts Time/Frequenc Time Time Step命令，进行如下图设置 QUST37QUST38n择GUI：Main Menu Solution Radiation OptsSolution Opt命令，在弹出的对话框里在STEF文本框中输入5.6E-8，在Toffset 文本框中输入273，其余选项采用默认设置，如下图所示，单击OK关闭对话框。nGUI：Main Menu SolutionDefine loads applythermal temperatureon nodes命令，在弹出的对话框中，在文本框中输入1，在Lab2 DOFs to be constrained列表框中选择TEMP选项，在Value Load Temp Value文本框中输入2000，单击OK。QUST39nGUI：Main Menu SolutionDefine loads applythermal temperatureon nodes命令，在弹出的对话框中，在文本框中输入2，在Lab2 DOFs to be constrained列表框中选择TEMP选项，在Value Load Temp Value文本框中输入0，单击OK。nGUI：Main Menu Solutionsolve Current LS命令，单击OK，进行求解。QUST40QUST41n查看求解结果nGUI：Main Menu General PostprocRead resultlast setnGUI：Main Menu General Postproclist resultsReaction Solu。在弹出的对话框中，在Lab Item to be listed 对话框中选择 Heat flow Heat选项。单击OK。QUST42n6.命令流文件n/Filename,exercise1n/title,radiant energy emissionnKeyw PR_THERM,1n/PREP7nET,1,LINK32nR,1,1,1,1nN,1nN,2nE,1,2nALLSELnFINISHn/SOLUQUST43nANTYPE,STATICnTIME,1nAUTOTS,ONnKBC,1nSTEF,5.67e-8nTOFFSET,273nD,1,TEMP,2000nD,2,TEMP,0nSOLVEnFINISHn/POST1nSET,LASTnPRRSOLnFINISHn/EXITQUST445.4 热辐射问题分析2n1.问题描述n两个等长度的同轴长圆柱体，截面如图所示，内外圆柱体的初始温度分别为1000、100，两者之间的辐射率为1，两圆柱体的材料性能参数如表所示，求两圆柱体间的辐射热流率。QUST45204080100QUST46导热系数K1导热系数K2密1密度2比热C1比热C2703078002700480350QUST47n2.问题分析n该问题为平面热辐射问题，整个热辐射问题为封闭系统，不需设置空间节点，在计算过程中使用AUX12辐射矩阵生成器，选择PLANE55 热分析单元，进行求解。QUST48n3.建立模型n（1）过滤菜单k简化菜单（过滤菜单）操作如下：kGUI：Main Menu Referenck在弹出的Referenc for GUI Filtering对话框中，选择 Thermal。单击OK。k命令：/COM,ThermalQUST49(2)选择单元k选择热分析单元，操作如下：kGUI：Main Menu Preprocessor Element Type Add/Edit/Deletek选择Thermal SolidQuad 4node 55选项，在Element type reference number 文本框中输入如1。k点击OK QUST50n(3)定义材料性能参数kGUI：Main Menu Preprocessor Material Props Material Models命令，出现Define Material Model Behavior 对话框。k在弹出的定义材料热传导率对话框中的KXX 栏键入“70k命令：MPDATA,KXX,1,70 QUST51n(4)定义密度nGUI：Main Menu Preprocessor Material Props Thermal Densityn在弹出密度定义对话框中的DENS栏键入“7800”。n命令：MPDATA,DENS,1,7800QUST52n(5)定义比热容nGUI：Main Menu Preprocessor Material Props Thermal Specific Heatn在弹出的定义比热容对话框中的C栏键入“480”。n命令：MPDATA,C,1,480 QUST53n(6)定义材料2的材料特性k在Define Material Model Behavior 对话框上单击Material New Model 命令Define Material ID 对话框，在文本框中输入2，单击OK。k定义材料热传导率对话框中的KXX 栏键入“70”。k命令：MPDATA,KXX,2,70QUST54n(7)定义密度nGUI：Main Menu Preprocessor Material Props Thermal Densityn在弹出密度定义对话框中的DENS栏键入“1000”。n命令：MPDATA,DENS,2,1000n(8)定义比热容nGUI：Main Menu Preprocessor Material Props Thermal Specific Heatn在弹出的定义比热容对话框中的C栏键入“350”。n命令：MPDATA,C,2,350QUST554.建立实体模型n(1)几何模型nGUI：Preprocessor Modeling Create Areas Circle Partial Annulusn在弹出的对话框中进行如下设置：WPX=0,WPY=0,Rad-1=0.1,Theta-1=0,Rad-2=0.2,Theta-2=90。n单击Apply，进行如下设置：WPX=0,WPY=0,Rad-1=0.4,Theta-1=0,Rad-2=0.5,Theta-2=90。QUST56n(2)网格划分n在设定网格尺寸参数，操作如下：kGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked Linesk在 Element Size On 菜单中选择线2、4k在Element Size On Picked Lines 对话框中，在NDIV 文本框中输入“5”，单击OKQUST57kGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked Linesk在 Element Size On 菜单中选择线1、3，在Element Size On Picked Lines 对话框中，在NDIV 文本框中输入“10”，单击OKkGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked LinesQUST58k在 Element Size On 菜单中选择线6、8，在Element Size On Picked Lines 对话框中，在NDIV 文本框中输入“6”，单击OKkGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked Linesk在 Element Size On 菜单中选择线5、7，在Element Size On Picked Lines 对话框中，在NDIV 文本框中输入“30”，单击OKkGUI：Main Menu Preprocessor Meshing Size CntrlsManualSize Lines Picked LinesQUST59n(3)网格划分nGUI：Main Menu Preprocessor MeshingMesh Tool，出现Mesh Tool菜单，在shape选项组中选中Quad和Mapped选项，单击Mesh，出现Mesh Areas，在文本框中输入1，单击OK。nGUI：Utility Menu Select EverythingnGUI：Main Menu Preprocessor MeshingMesh attributes Default attribs，出现Meshing attributes 菜单，在MAT下拉列表框中选择2，单击OK.QUST60kGUI：Main Menu Preprocessor MeshingMesh Tool，出现Mesh Tool菜单，单击Mesh，出现Mesh Areas，在文本框中输入2，单击OK。kGUI：Main Menu Preprocessor ModelingReflectAreas 命令，出现Reflect Areas菜单，单击Pick All 按钮，在弹出的对话框中，选择Y-Z Plane X 单选按钮，在NOELEM下拉列表框中选择Areas and mesh选项，在IMOVE下拉列表框中选择Copied选项，单击OK。kGUI：Main Menu Preprocessor ModelingReflectAreas 命令，出现Reflect Areas菜单，单击Pick All 按钮，在弹出的对话框中，选择X-Z Plane Y 单选按钮，其余采用默认值，单击OK。n注意:以上两步间不能采用Apply 按钮过渡QUST61nGUI：Main Menu Numbering Ctrls Merge Items 命令，出现Merge Coincident or Equivalently Defined Items 对话框，在Label Type of Item to be merge 在下拉列表框中选择Keypoints选项，单击OK。n选择UTILITY MenuSelect everythying 命令。QUST62n生成AUX12辐射矩阵n 选择Main Menu Radiation Opt Radiosity MethSolution Opt命令，出现Radiation Solution Options对话框，在STEF 文本框中输入5.67E-8，在TOFFSET文本框中输入0.001，在SPCTEMP/SPCNOD下拉列表框中选择temperature选择，在Value文本框中输入20，其余采用默认值，单击OK。如图所示QUST63QUST64n加载求解n（1）选择Transient分析，操作如下：nGUI：Main Menu Preprocessor Loads Analysis Type New Analysisn选择Steady-state 分析，单击OK。nGUI：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。QUST65nGUI：Main Menu SolutionDefine loads applythermaltemperatureon lines，出现Apply TEMP on lines 菜单，在文本框中输入5、13、21、29，单击Ok按钮，出现Apply TEMP on lines对话框，在LAB2列表框中选择TEMP，在VALUE文本框中输入“100”，选择KEXPND选项，使其状态为“Yes”nGUI：Main Menu SolutionDefine loads applythermaltemperatureon lines，出现Apply TEMP on lines 菜单，在文本框中输入3、11、19、27，单击Ok按钮，出现Apply TEMP on lines对话框，在LAB2列表框中选择TEMP，在VALUE文本框中输入“1000”，选择KEXPND选项，使其状态为“Yes”。QUST66QUST67n选择UTILITY MenuSelect everythying 命令。nGUI：Main Menu Solution Load Step Opts Time/Frequenc Time Time Step命令，进行如下图设置nGUI：Main Menu Solutionsolve Current LS命令，单击OK，进行求解。QUST68QUST696.查看求解结果(1)nGUI：Main Menu General PostprocRead resultlast setnUtility Menu SelectEntities 命令，在对话框中选择Element选择，在第2个下拉列表框中选择By Attributes 选项，在第3个下拉列表框中选择Material num单选按钮，在Min,Max,Inc 文本框中输入1，在第5个下拉列表框中选择From Full单选按钮，如下图所示nUtility Menu SelectEntities 命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attached to 选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择From Full，如下图所示QUST70QUST71nGUI：Main Menu General PostprocPlot ResultsContour Plot SelectEntities 命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attached to 选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择From Full。nUI：Main Menu General PostprocPlot ResultsContour Plot SelectEntities 命令，在对话框中选择Element选择，在第2个下拉列表框中选择By Attributes 选项，在第3个下拉列表框中选择Material num单选按钮，在Min,Max,Inc 文本框中输入1，在第5个下拉列表框中选择From Full单选按钮。nUtility Menu SelectEntities 命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attached to 选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择From Full。QUST76nGUI：Main Menu General PostprocList ResultsReaction Solu命令，在弹出的对话框中保持默认值，单击OK，显示内圆柱体上所有节点的辐射换热量，如图所示。QUST77n7.命令流n/PREP7nET,1,PLANE55nMP,DENS,1,7800nMP,KXX,1,70nMP,C,1,480nMP,DENS,2,2700nMP,KXX,2,30nMP,C,2,350nCYL4,0,0,0.1,90,0.2nCYL4,0,0,0.4,90,0.5nQUST78n APLOTn/PNUM,LINE,1n/PNUM,AREA,1nLSEL,S,2,4,2nLESIZE,ALL,5nLSEL,S,1,3,2nLESIZE,ALL,10nLSEL,S,6,8,2nLESIZE,ALL,6nLSEL,S,5,7,2nLESIZE,ALL,30nMSHAPE,0,2DQUST79 MSHKEY,1kAMESH,1kAPLOTkMAT,2kAMESH,2kARSYM,X,ALLkARSYM,Y,ALLkNUMMRG,NODEkNUMMRG,KPkEPLOTkALLSELQUST80n/AUX12nSTEF,5.67E-8nTOFFSET,273nRADOPT,0.001nSPCTEMP,1,20nFINISHQUST81k /SOLUkANTYPE,STATICkLPLOTkLSEL,S,1,7,6kLSEL,A,9,15,6kLSEL,A,17,23,6kLSEL,A,25,31,6kSFL,ALL,RDSF,1,1kLSEL,S,5,29,8kDL,ALL,TEMP,100,1kLSEL,S,3,27,8kDAL,ALL,TMEP,1000,1kALLSELkDETIM,0.5QUST82kSOLVEkFINISHk/POST1kSET,LASTkESEL,S,MAT,1kNSLE,SkPLNSOL,TEMPkESEL,S,MAT,2kNSLE,SkPLNSOL,TEMPkESEL,S,MAT,1n QUST83kNSLE,SkPRRSOLkESEL,S,MAT,2kNSLE,SkPRRSOLkALLSELkFINISH