ANSYS(热辐射_第5节)

5热辐射分析5热辐射分析5.1热辐射分析简介辐射是一种通过电磁波传递能量的方式。电磁波以光速传播且无需任何介质。热辐射仅为电磁波谱中的一小段。因为由于热辐射引起的热流与物体 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 面绝对温度的四次方成正比，因此热辐射分析是高度非线性的。辐射热传递是通过电磁波传递热能的方法。热辐射的电磁波波长为0.1到100mm。不象其他热传递方式需要介质，辐射在真空中(如外层空间)效率最高。对于半透明体(如玻璃),辐射是三维实体现象,因为辐射从体中发散出。对于不透明体，辐射主要是平面现象因为几乎所有内部辐射都被实体吸收了。ANSYS可以模拟不透明体间的辐射，所以我们将讨论范围限制在平面辐射现象上。平面可以理想化为散射或反射装置。散射装置会将辐射均匀反射到所有方向，而不管辐射源的方位:通常情况下，平面可以被理想化为散射或反射面。反射平面会将辐射以近乎镜象的方式反射:没有实际的平面是真正的散射或反射面。比较灰暗的平面接近散射面，高度抛光的平面接近反射面。为了简化计算，平面的辐射特性可以在所有的波长和方向平均。因此，在散射和反射平面之间没有差别。两平面间的辐射热传递与它们平面绝对温度差的四次方成正比:1.黑体黑体是理想化的平面，用来与实际平面进行比较。黑体的特性:黑体吸收所有的偶然辐射(没有反射),不管波长和方向。黑体为纯粹的发射器。对于给定的波长和温度，没有平面比黑体发射更多的能量。黑体是纯粹的散射发射器;辐射在所有方向均一致。因此，对于黑体:αB=εB=1.灰体实际平面叫做“灰体”因为他们不象黑体。3.ANSYS和辐射ANSYS中关于辐射的重要假设和方法:ANSYS认为辐射是平面现象，因此适合用不透明平面建模。ANSYS不直接计入平面反射率。考虑到效率，假设平面吸收率和发射率相等(a=ε)。因此，只有发射率特性需要在ANSYS辐射分析中定义。ANSYS不自动计入发射率的方向特性，也不允许发射率定义随波长变化。发射率可以在某些单元中定义为温度的函数。到此为止我们只讨论了单独的辐射平面。但是，在研究实际问题时，我们通常要考虑多个辐射平面的相互作用。要考虑的平面越多，问题越复杂:4.形状因子“形状因子”由相互辐射的两个平面(i和j)定义。它的定义是由于从一个平面(i)发射的辐射能偶然施加到另一个平面(j)上而得到。两个平面的形状因子是面积，方向和距离的函数。Ai、Aj表示表面I与表面J的面积，γ表示面单元dAi与面单元dAj之间的距离，θi表示面单元dAj的法线Nj与两面但与连线的夹角，Ni表示面单元dAi的法线；Nj表示面单元dAj的法线角系数Fij具有以下特征：与面距离的平方成反比；与cosθi成正比，在θi=0时表面I辐射出最大的热量。与cosθj成正比，在θj=0时表面I辐射出最大的热量。若两表面彼此“看不见”（即cosθi=0且cosθj），则彼此的角系数等于0；根据相互作何用原理，从任何平面发射的能量必须守恒，对于任意两表面均有AiFij=AjFji热辐射计算方程：两个表面之间的热辐射计算公式为：上式中各参量的物理意义如下：Q：表面I的传热率：Stefan-Bolzman常数：有效热辐射率；：表面I的面积；、：表面I与表面J的绝对温度值。、5.2热辐射问题分析ANSYS提供了三种方法分析热辐射问题：（1）用LINK31，辐射线单元，分析两个点或多对点之间的热辐射（2）用表面效应单元SURF19或SURF22，分析点对面的热辐射（3）用AUX12，热辐射矩阵生成器，分析面与面之间的热辐射以上三种方法既可用于稳态热分析，也可用于瞬态传热分析。1.点-点问题模拟两节点或多对节点的热辐射问题时，应使用热辐射线单元LINK31。LINK31作为两节点非线性单元可以计算两点之间因辐射引起的热交换。使用此单元需指明以下常数：有效的辐射表面面积角系数辐射率Stefan-bolzmann常数。2.点-面问题模拟点面之间的热辐射问题，通常采用表面效应单元-SURF151和SURF152。其中SURF151用于2D单元和SURF152用于3D单元。表面效应单元利用实体表面的节点形成单元，并且直接覆盖在实体单元的表面。使用表面效应单元可以更灵活的在实体表面施加热载荷。例如，热流密度与对流可以施加在同一外表面，但ANSYS在计算过程中仅读取最后施加的面载荷进行计算。为避免ANSYS只读取一种载荷，可以利用实体单元承受热流密度，而表面效应单元承受对流载荷。SURF151实常数:FORMF(角系数)、SBCONST（Stefan-Bolzmann常数）材料属性：DENS（密度）、EMIS（辐射率）表面载荷：对流、热流密度。体载荷：生热率SURF152实常数:FORMF(角系数)、SBCONST（Stefan-Bolzmann常数）材料属性：DENS（密度）、EMIS（辐射率）表面载荷：对流、热流密度。体载荷：生热率注意：热辐射分析要注意温度的单位制，因为计算热辐射使用的温度单位是绝对温度。如果在加载时使用的是华氏温度，就要设置460的差值；如果为摄氏温度，差值为273。3.使用LINK31单元的 注意事项 软件开发合同注意事项软件销售合同注意事项电梯维保合同注意事项软件销售合同注意事项员工离职注意事项 LINK31是一个两节点非线性线单元，用于计算由辐射引起的两点之间的热传递。此单元要求键入如下的实常数：有效的热辐射面积形状系数辐射率Stefan-Boltzmann常数4.使用表面效应单元的注意事项表面效应单元可以方便地分析点与面之间的辐射传热。SURF19用于两维模型，SURF22用于三维模型。单元应设置为包含辐射KEYOPT（9）。5.使用AUX12—辐射矩阵生成器此方法用于计算多个辐射面之间的辐射传热。这种方法生成辐射面之间形状系数矩阵，并将此矩阵作为超单元用于热分析。AUX12方法由三个步骤组成，即首先定义辐射面，然后生成辐射矩阵，最后在热分析中使用辐射矩阵。下面介绍使用AUX12求解热辐射问题的步骤：(1)定义辐射面首先建立有限元模型在辐射面上覆盖一层SHELL57辐射面往往是3D模型中的面或2D模型中的边，因此在辐射表面用SHELL57（3D）或LINK32（2D）划分网格。最好的方法是先选择辐射表面的节点，然后用如下方法创建SHELL57或LINK32单元：命令：ESURFGUI：MainMenu>Preprocessor>Modeling>Create>Elements>Surf/Contact>Surf-suf>GeneralSurface注意：辐射面上的SHELL57或LINK32单元与节点必须与实体单元相吻合，否则计算的结果是不正确的。生成的SHELL57或LINK32单元的取向也很重要。AUX12假定辐射的方向是SHELL57的+Z向或LINK32的+Y向。因此在生成SHELL57或LJNK32单元时要注意节点的排列顺序。如果所分析的系统是开放的，即一个面所辐射的热能未被模型中其它的面吸收，则必须定义一个空间节点，用于吸收损失的辐射热量。这个节点的位置是任意的。对于封闭的系统，不应定义空间节点。(2)生成辐射矩阵生成辐射矩阵首先需要进入AUX12，然后选择所需节点和单元，确定模型的维数。进入AUX12命令：/AUX12GUI：MainMenu>RadiationMatrix选择组成辐射面的节点和单元比较方便的方法是根据单元类型选择单元，并选择单元上的节点。确定模型是3D还是2D命令：GEOMGUI：MainMenu>RadiationMatrix>OtherSettingAUX12用不同的算法计算2D或3D模型的形状系数。AUX12默认为3D。2D分为纯平面或轴对称，默认为纯平面。定义每个辐射面的辐射率命令：EMISGUI：MainMenu>RadiationMatrix>Emissivities辐射面的辐射率ANSYS默认为1。定义Stefan-Boltzmann常数命令：STEFGUI：MainMenu>RadiationMatrix>OtherSettingsStefan-Boltzmann常数(3)确定状系数命令：VTYPEGUI：MainMenu>RadiationMatrix>WriteMatrix用什么方式计算形形状系数，ANSYS提供两种选择，即选择是隐藏还是非隐藏方法。非隐藏方法计算每两个单元之间的形状系数，无论它们之间有无障碍；隐藏方法默认）用一种隐藏线算法判断两辐射面之间是否“可见”，如果可见则计算形状系(4)定义空间节点：命令：SPACEGUI：MainMenu>RadiationMatrix>OtherSettings如为开放系统.(5)计算辐射矩阵并写入jobename.sub文件命令：WRITEGUI：MainMenu>RadiationMatrix>WriteMatrix(6)选择所有的节点和单元5.3热辐射问题分析11.问题描述如图所示一黑体，表面积1m2，形状系数和辐射率为1，温度为2000℃，周围环境温度为0℃，求黑体的辐射热流率。黑体结构示意图2.问题分析选择LINK31热辐射单元进行求解。3.建立模型(1)过虑菜单简化菜单（过虑菜单）操作如下：GUI：MainMenu>Referenc在弹出的ReferencforGUIFiltering对话框中，选择Thermal。单击OK。命令：/COM,Thermal(2)选择单元选择热分析单元，操作如下：GUI：MainMenu>Preprocessor>ElementType>Add/Edit/Delete选择ThermalLink>3Dradiation31选项，在Elementtypereferencenumber文本框中输入如1。GUI：MainMenu>Preprocessor>Realcontants>Add/edit/delete命令，点击Add按钮，在Elementtypeforrealcontants对话框，单击OK按钮，进行如图设置(3)建立模型首先进入DefineMaterialModelBehavior对话框，操作如下：GUI：MainMenu>Preprocessor>Modeling>Create>nodes>onactivecs命令，在弹出的对话框中，在NodeNumbers文本框中输入1，在X，Y，Z文本框中依次输入0、0、0。点击Apply按钮，在NodeNumbers文本框中输入2，在X，Y，Z文本框中依次输入0、0、0。GUI：MainMenu>Preprocessor>Modeling>Create>Element>AutoNumbered>ThruNodes，出现Elementfromnodes菜单，在文本框中输入1，2，单击OK。Unitity>File>Saveas命令进行存盘。4.加载求解选择Transient分析，操作如下：GUI：MainMenu>Preprocessor>Loads>AnalysisType>NewAnalysis选择Steady-state分析，单击OK。进入TimeandtimeStepOptios对话框，操作如下：GUI：MainMenu>Solution>LoadStepOpts>Time/Frequenc>Time–TimeStep命令，进行如下图设置择GUI：MainMenu>Solution>RadiationOpts>SolutionOpt命令，在弹出的对话框里在STEF文本框中输入5.6E-8，在Toffset文本框中输入273，其余选项采用默认设置，如下图所示，单击OK关闭对话框。GUI：MainMenu>Solution>Defineloads>apply>thermal>temperature>onnodes命令，在弹出的对话框中，在文本框中输入1，在Lab2DOFstobeconstrained列表框中选择TEMP选项，在ValueLoadTempValue文本框中输入2000，单击OK。GUI：MainMenu>Solution>Defineloads>apply>thermal>temperature>onnodes命令，在弹出的对话框中，在文本框中输入2，在Lab2DOFstobeconstrained列表框中选择TEMP选项，在ValueLoadTempValue文本框中输入0，单击OK。GUI：MainMenu>Solution>solve>CurrentLS命令，单击OK，进行求解。查看求解结果GUI：MainMenu>GeneralPostproc>Readresult>lastsetGUI：MainMenu>GeneralPostproc>listresults>ReactionSolu。在弹出的对话框中，在LabItemtobelisted对话框中选择HeatflowHeat选项。单击OK。6.命令流文件/Filename,exercise1/title,radiantenergyemissionKeywPR_THERM,1/PREP7ET,1,LINK32R,1,1,1,1N,1N,2E,1,2ALLSELFINISH/SOLUANTYPE,STATICTIME,1AUTOTS,ONKBC,1STEF,5.67e-8TOFFSET,273D,1,TEMP,2000D,2,TEMP,0SOLVEFINISH/POST1SET,LASTPRRSOLFINISH/EXIT5.4热辐射问题分析21.问题描述两个等长度的同轴长圆柱体，截面如图所示，内外圆柱体的初始温度分别为1000、100，两者之间的辐射率为1，两圆柱体的材料性能参数如表所示，求两圆柱体间的辐射热流率。204080100导热系数K1导热系数K2密ρ1密度ρ2比热C1比热C27030780027004803502.问题分析该问题为平面热辐射问题，整个热辐射问题为封闭系统，不需设置空间节点，在计算过程中使用AUX12辐射矩阵生成器，选择PLANE55热分析单元，进行求解。3.建立模型（1）过滤菜单简化菜单（过滤菜单）操作如下：GUI：MainMenu>Referenc在弹出的ReferencforGUIFiltering对话框中，选择Thermal。单击OK。命令：/COM,Thermal(2)选择单元选择热分析单元，操作如下：GUI：MainMenu>Preprocessor>ElementType>Add/Edit/Delete选择ThermalSolid>Quad4node55选项，在Elementtypereferencenumber文本框中输入如1。点击OK(3)定义材料性能参数GUI：MainMenu>Preprocessor>MaterialProps>MaterialModels命令，出现DefineMaterialModelBehavior对话框。在弹出的定义材料热传导率对话框中的KXX栏键入“70”。命令：MPDATA,KXX,1,,70(4)定义密度GUI：MainMenu>Preprocessor>MaterialProps>Thermal>Density在弹出密度定义对话框中的DENS栏键入“7800”。命令：MPDATA,DENS,1,,7800(5)定义比热容GUI：MainMenu>Preprocessor>MaterialProps>Thermal>SpecificHeat在弹出的定义比热容对话框中的C栏键入“480”。命令：MPDATA,C,1,,480(6)定义材料2的材料特性在DefineMaterialModelBehavior对话框上单击Material>NewModel命令DefineMaterialID对话框，在文本框中输入2，单击OK。定义材料热传导率对话框中的KXX栏键入“70”。命令：MPDATA,KXX,2,,70(7)定义密度GUI：MainMenu>Preprocessor>MaterialProps>Thermal>Density在弹出密度定义对话框中的DENS栏键入“1000”。命令：MPDATA,DENS,2,,1000(8)定义比热容GUI：MainMenu>Preprocessor>MaterialProps>Thermal>SpecificHeat在弹出的定义比热容对话框中的C栏键入“350”。命令：MPDATA,C,2,,3504.建立实体模型(1)几何模型GUI：Preprocessor>Modeling>Create>Areas>Circle>PartialAnnulus在弹出的对话框中进行如下设置：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。(2)网格划分在设定网格尺寸参数，操作如下：GUI：MainMenu>Preprocessor>Meshing>SizeCntrls>ManualSize>Lines>PickedLines在ElementSizeOn菜单中选择线2、4在ElementSizeOnPickedLines对话框中，在NDIV文本框中输入“5”，单击OKGUI：MainMenu>Preprocessor>Meshing>SizeCntrls>ManualSize>Lines>PickedLines在ElementSizeOn菜单中选择线1、3，在ElementSizeOnPickedLines对话框中，在NDIV文本框中输入“10”，单击OKGUI：MainMenu>Preprocessor>Meshing>SizeCntrls>ManualSize>Lines>PickedLines在ElementSizeOn菜单中选择线6、8，在ElementSizeOnPickedLines对话框中，在NDIV文本框中输入“6”，单击OKGUI：MainMenu>Preprocessor>Meshing>SizeCntrls>ManualSize>Lines>PickedLines在ElementSizeOn菜单中选择线5、7，在ElementSizeOnPickedLines对话框中，在NDIV文本框中输入“30”，单击OKGUI：MainMenu>Preprocessor>Meshing>SizeCntrls>ManualSize>Lines>PickedLines(3)网格划分GUI：MainMenu>Preprocessor>Meshing>MeshTool，出现MeshTool菜单，在shape选项组中选中Quad和Mapped选项，单击Mesh，出现MeshAreas，在文本框中输入1，单击OK。GUI：UtilityMenu>Select>EverythingGUI：MainMenu>Preprocessor>Meshing>Meshattributes>Defaultattribs，出现Meshingattributes菜单，在MAT下拉列表框中选择2，单击OK.GUI：MainMenu>Preprocessor>Meshing>MeshTool，出现MeshTool菜单，单击Mesh，出现MeshAreas，在文本框中输入2，单击OK。GUI：MainMenu>Preprocessor>Modeling>Reflect>Areas命令，出现ReflectAreas菜单，单击PickAll按钮，在弹出的对话框中，选择Y-ZPlaneX单选按钮，在NOELEM下拉列表框中选择Areasandmesh选项，在IMOVE下拉列表框中选择Copied选项，单击OK。GUI：MainMenu>Preprocessor>Modeling>Reflect>Areas命令，出现ReflectAreas菜单，单击PickAll按钮，在弹出的对话框中，选择X-ZPlaneY单选按钮，其余采用默认值，单击OK。注意:以上两步间不能采用Apply按钮过渡GUI：MainMenu>NumberingCtrls>MergeItems命令，出现MergeCoincidentorEquivalentlyDefinedItems对话框，在LabelTypeofItemtobemerge在下拉列表框中选择Keypoints选项，单击OK。选择UTILITY>Menu>Select>everythying命令。生成AUX12辐射矩阵选择MainMenu>RadiationOpt>RadiosityMeth>SolutionOpt命令，出现RadiationSolutionOptions对话框，在STEF文本框中输入5.67E-8，在TOFFSET文本框中输入0.001，在SPCTEMP/SPCNOD下拉列表框中选择temperature选择，在Value文本框中输入20，其余采用默认值，单击OK。如图所示加载求解（1）选择Transient分析，操作如下：GUI：MainMenu>Preprocessor>Loads>AnalysisType>NewAnalysis选择Steady-state分析，单击OK。GUI：MainMenu>Solution>Defineloads>apply>thermal>radiation>onlines命令，出现ApplyRDSFonlines菜单，在文本框中输入1、7、9、15、17、23、25、31，单击Ok按钮，出现ApplyRDSFonlines对话框，在VALUEEmissivity文本框中输入1，在VALUE2的文本框中输入1，单击OK。GUI：MainMenu>Solution>Defineloads>apply>thermal>temperature>onlines，出现ApplyTEMPonlines菜单，在文本框中输入5、13、21、29，单击Ok按钮，出现ApplyTEMPonlines对话框，在LAB2列表框中选择TEMP，在VALUE文本框中输入“100”，选择KEXPND选项，使其状态为“Yes”GUI：MainMenu>Solution>Defineloads>apply>thermal>temperature>onlines，出现ApplyTEMPonlines菜单，在文本框中输入3、11、19、27，单击Ok按钮，出现ApplyTEMPonlines对话框，在LAB2列表框中选择TEMP，在VALUE文本框中输入“1000”，选择KEXPND选项，使其状态为“Yes”。选择UTILITYMenu>Select>everythying命令。GUI：MainMenu>Solution>LoadStepOpts>Time/Frequenc>Time–TimeStep命令，进行如下图设置GUI：MainMenu>Solution>solve>CurrentLS命令，单击OK，进行求解。6.查看求解结果(1)GUI：MainMenu>GeneralPostproc>Readresult>lastsetUtilityMenu>Select>Entities命令，在对话框中选择Element选择，在第2个下拉列表框中选择ByAttributes选项，在第3个下拉列表框中选择Materialnum单选按钮，在Min,Max,Inc文本框中输入1，在第5个下拉列表框中选择FromFull单选按钮，如下图所示UtilityMenu>Select>Entities命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attachedto选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择FromFull，如下图所示GUI：MainMenu>GeneralPostproc>PlotResults>ContourPlotSelect>Entities命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attachedto选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择FromFull。UI：MainMenu>GeneralPostproc>PlotResults>ContourPlotSelect>Entities命令，在对话框中选择Element选择，在第2个下拉列表框中选择ByAttributes选项，在第3个下拉列表框中选择Materialnum单选按钮，在Min,Max,Inc文本框中输入1，在第5个下拉列表框中选择FromFull单选按钮。UtilityMenu>Select>Entities命令，在对话框中选择Nodes选择，在第2个下拉列表框中选择Attachedto选项，在第3个下拉列表框中选择Elements单选按钮，在第4个选项组中选择FromFull。GUI：MainMenu>GeneralPostproc>ListResults>ReactionSolu命令，在弹出的对话框中保持默认值，单击OK，显示内圆柱体上所有节点的辐射换热量，如图所示。7.命令流/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.5APLOT/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,,,30MSHAPE,0,2DMSHKEY,1AMESH,1APLOTMAT,2AMESH,2ARSYM,X,ALLARSYM,Y,ALLNUMMRG,NODENUMMRG,KPEPLOTALLSEL/AUX12STEF,5.67E-8TOFFSET,273RADOPT,,0.001SPCTEMP,1,20FINISH/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.5SOLVEFINISH/POST1SET,LASTESEL,S,MAT,,1NSLE,SPLNSOL,TEMPESEL,S,MAT,,2NSLE,SPLNSOL,TEMPESEL,S,MAT,,1NSLE,SPRRSOLESEL,S,MAT,,2NSLE,SPRRSOLALLSELFINISH