搜索
第六章热辐射分析6.1热辐射的定义热辐射是一种通过电磁波传递热能的方式。电磁波以光的速度进行传递，而能量传递与辐射物体之间的介质无关。热辐射只在电磁波的频谱中占小部分的带宽。由于辐射产生的热流与物体表面的绝对温度的四次方成正比，因此热辐射有限元分析是高度非线性的。物体表面的辐射遵循Stefan-Boltzmann定律：式中：—物体表面的绝对温度；—Stefan-Boltzmann常数，英制为0.119×10-10BTU/hr-in-R，公制为5.67×10-86.2基本概念下面是对辐射分析中用到的一些术语的定义：黑体黑体被定义为在任意温度下，吸收并发射最大的辐射能的物体；通常的物体为“灰体”，即ε1；在某些情况下，辐射率（黑度）随温度变化；辐射率（黑度）物体表面的辐射率（黑度）定义为物体表面辐射的热量与黑体在同一表面辐射热量之比。式中：－辐射率（黑度）－物体表面辐射热量－黑体在同一表面辐射热量形状系数形状系数用于计算两个面之间的辐射热交换，在ANSYS中，可以用隐藏/非隐藏的方法计算2维和三维问题，或者用半立方的方法来计算3维问题。表面I与表面J之间的形状系数为：形状系数是关于表面面积、面的取向及面间距离的函数；由于能量守恒，所以：根据相互原理：由辐射矩阵计算的形状系数为：式中：－单元法向与单元I,J连线的角度－单元I,J重心的距离有限单元模型的表面被处理为单元面积dAI及dAJ，然后进行数字积分。辐射对在辐射问题中，辐射对由一些相互之间存在辐射的面组成，可以是开放的或是闭合的。在ANSYS中，可以定义多个辐射对，它们相互之间也可以存在辐射ANSYS使用辐射对来计算一个辐射对中各面间的形状系数；每一个开放的辐射对都可以定义自己的环境温度，或是向周围环境辐射的空间节点。Radiosity求解器当所有面上的温度已知时，Radiosity求解器方法通过计算每一个面上的辐射热流来得到辐射体之间的热交换。而面上的热流为接下来的热传导分析提供了有限元模型的边界条件。重复上面的过程，就会由于新的时间步或者新的迭代循环会得到新的热流边界条件，从而计算出新的温度分布。在计算中使用的每个表面的温度必须是均匀的，这样才能满足辐射模型的条件。6.3分析热辐射问题针对不同的情况ANSYS为热辐射分析提供了四种方法。热辐射线单元（LINK31），模拟两节点间（或多对节点）间辐射；表面效应单元（SURF151及SURF152），模拟点对面（线）的辐射；利用AUX12生成辐射矩阵，模拟更一般的面与面（或线与线）的辐射(只有ANSYS/MultiphysicsANSYS/Mechanical和ANSYS/Professional这些产品提供辐射矩阵生成器)；Radiosity求解器方法，求解二维、三维面与面之间的热辐射，该方法对所有含温度自由度的二维和三维单元都适用。(只有ANSYS/Multiphysics，ANSYS/Mechanical和ANSYS/Professional这些产品提供Radiosity求解器)可以将上面四种辐射方法中的任何一种用于稳态或瞬态热分析中。辐射是一种非线性现象，因此需要进行平衡迭代来得到收敛解。6.4节点间的热辐射非线性线单元LINK31用于计算两节点间或多对节点间的简单辐射热传递，节点的位置是任意的，可作为其它单元的节点。LINK31需要定义如下数据：材料属性：EMIS辐射率（可以随温度变化）实常数：AREA（Ai）（有效辐射面积）FORMF（Fij）（形状系数）SBCONST（Stefan-Boltzman常数）有关LINK31的使用实例，请参考《ANSYS校验手册》：VM106RadiantenergyemissionVM107Thermocoupleradiation6.5点与面间的热辐射应用表面效应单元可以方便地计算点与面间的辐射，包括2D的SURF151及3D的SURF152单元：首先在实体单元的辐射表面覆盖一层表面效应单元；单元关键选项KEYOPT(9)激活这些单元的热辐射分析功能；如果设置KEYOPT(9)＝1（缺省值），则可在单元实常数中定义形状系数；如果设置KEYOPT(9)＝2或3，则程序基于单元表面的法向与附加节点的位置关系、考虑余弦效应，计算形状系数。使用表面单元进行热辐射分析的实例，请参考《ANSYS校验手册》：VM192Coolingofabilletbyradiation6.6AUX12―辐射矩阵生成器只有ANSYS/Multiphysics，ANSYS/Mechanical，ANSYS/Thermal等产品提供AUX12，用于计算多个辐射面之间的辐射热传递。AUX12生成辐射面间的形状系数矩阵，并将此矩阵作为超单元用于热分析。分析模型可包含隐藏面或部分隐藏面，以及能吸收辐射能量的“空间节点”。使用AUX12计算面与面之间的辐射可分为以下三个步骤：定义辐射面生成辐射矩阵将辐射矩阵用于热分析6.6.1定义辐射面分析类型在2D有限元模型的辐射边上覆盖一层LINK32单元，或在3D有限元模型的辐射面上覆盖一层SHELL57单元，步骤如下：在前处理中创建热分析模型。由于辐射表面不支持对称条件，包含辐射的模型就无法利用几何上的对称性，因此必须建立完整的分析模型。二维和三维的辐射面示意图6-1如下：在辐射面上覆盖一层SHELL57（3D）或LINK32（2D）单元，如图6-2所示。最好的方法是先选择辐射面的节点，然后用下面的方法创建面上的单元：命令：ESURFGUI:MainMenuPreprocessorCreateElementsSurfEffectExtraNodeMainMenuPreprocessorCreateElementsSurfEffectNoextraNode在执行上述命令之前，要确定首先激活了相应的单元类型及材料属性，如果不同辐射面的辐射率不同，建议用材料编号来区分这些面。还要注意SHELL57或LINK32的节点一定要与相应实体单元对应节点编号重合。否则计算结果将是错误的。所覆盖的SHELL57或LINK32单元的方向是非常重要的。AUX12假设辐射方向是SHELL57单元坐标系的正Z向或LINK32单元坐标系的正Y向，因此必须正确地划分覆盖层的网格以使辐射面相对。单元的方向是由节点的排列顺序决定的，如下图所示:图6-3覆盖单元的方向显示单元辐射方向的方法如下：命令：/PSYMB,ESYS,1GUI:UtilityMenuPlotCtrlsSymbols,将ESYSElementCoordinate设置为ON。定义一个空间节点，用于吸收没有被模型中其它辐射面吸收的辐射能量。这个节点的位置是任意的，对于一个开放系统通常需要空间节点，而对于封闭系统则不得设置空间节点6.6.2生成辐射矩阵计算辐射矩阵可按如下的步骤完成：1．进入Aux12命令：/AUX12GUI：MainMenuRadiation2.选择构成辐射面的节点和单元。较简便的方法是根据单元属性选择（如单元类型）选择单元，然后选择所有Attachedto单元的节点（同时应注意将空间节点也选择进来）:命令：ESEL,S,TYPE和NSELGII:UtilityMenuSelectEntities3.确定所分析的模型是3D还是2D命令：GEOMGUI:MainMenuRadiationMatrixOtherSettingsAUX12采用不同的算法分别计算2D和3D模型的形状系数，AUX12默认为3D。2D可以是平面的（NDIV＝0），也可以是轴对称的（NDIV0），缺省为平面的。轴对称模型在内部展成3D，NDIV是圆周方向分割数量。例如NDIV＝10，则每段为36度。4．确定辐射率（缺省为1.0）：命令：EMISGUI:MainMenuRadiationMatrixEmmisivities5.定义Stefan-Baltzmann常数（缺省情况下，英制单位为0.119×10-10Btu/hr-in2-R4，国际单位制为5.67×10－8W/m2K4）。命令：STEFGUI：MainMenuRadiationMatrixOtherSettings6.确定计算形状系数的方法。命令：VTYPEGUI:MainMenuRadiationMatrixWriteMatrix可选择非隐藏或隐藏两种方法之一，非隐藏方法计算每个单元对其它单元的形状系数，无论两单元之间是否有阻碍；隐藏方法（默认）首先用“隐藏线”算法确定两单元之间是否“可见”，如果目标单元与辐射单元的辐射方向指向对方，而且设有其它单元阻碍，则它们是“可见”的，形状系数按如下方法计算：每一个辐射单元被封闭成一个半径为单位值的半球（3D）或半圆（2D）；所有的目标单元向这个半球或半圆投影；一定数量（默认为20）的射线由辐射单元面投向半球或半圆。这样，形状系数就是投到投影面上的射线数量与辐射面发出的射线的数量之比，通常设定的射线数量越多，形状系数的精度越高。可以通过设定VTYPE命令的变量NZONE或上述的菜单来设定射线数量7．如果有必要（例如开放系统），应指定空间节点：命令：SPACEGUI:MainmenuRaditionMatrixOtherSettings8．将辐射矩阵写到文件jobname.sub中,如果想要写更多的辐射矩阵，为不同的矩阵指定不同的文件名：命令：WRITEGUI:MainMenuRadiationMatrixWriteMatrix如果需要打印出辐射矩阵，应在执行Write命令之前执行命令：mprint,1。9．选择所有节点及单元命令：ALLSELGUI:UtilityMenuSelectEverything现在就已经将辐射矩阵作为一个超单元写入到一个文件中了。6.6.3使用辐射矩阵进行热分析生成了辐射矩阵之后，重新进入前处理器，定义作为超单元的辐射矩阵。步骤如下：1．重新进入前处理器，选择单元MATRIX50（超单元），并设置单元Keyoption为热辐射分析。命令：/PREP7GUI：MainMenuPreprocessor2．设置缺省单元类型为超单元命令：TYPEGUI:MainMenuPreprocessorCreateElementElemAttributes3．读入辐射超单元矩阵命令：SEGUI:MainMenuReprocessorCreateElements–Superelements-From.SUB4．不选择或删除用于生成辐射矩阵的SHELL57或LINK32单元,因为在热分析中已经不再需要了。命令：EDELEGUI:MainMenuPreproccssorModelingDeleteElements5．进入求解器，定义空间节点的热边界条件，空间节点的典型热边界为温度（环境温度），也可能是热流率。边界条件应能够反映被模拟的环境的真实情况。命令：D,FGUI:MainMenuSolutionLoads-Apply…6．其它步骤与普通热分析相同6.7使用空间节点的几点建议尽管模拟热辐射并不总是需要定义空间节点，但使用或不使用空间节点可能会明显影响计算精度，分析中请注意有关空间节点的如下几点：6.7.1对于非隐藏方法用非隐藏方法计算形状系数，不对空间节点做特别的考虑，也可以得到系统足够精确的解。通常对于封闭系统不应定义空间节点；而对于开放系统应当定义。只有当开放系统中含有灰体（辐射率小于1）时，才必须定义一个空间节点，以保证计算精度。6.7.1对于隐藏方法AUX12中形状系数计算的精度会影响到空间节点的辐射计算，由于计算的误差在空间节点上累积，在封闭或接近封闭系统中空间节点形状系数的相对误差会过大。使用隐藏方法时，可能会需要增大计算形状系数时的射线数量，并细化网格，以便得到更精确的形状系数。如果上述方法不能实施，可考虑如下建议：对于封闭系统，即所有的辐射面形成一个封闭空间，不向外界辐射，不要使用空间节点。如果问题的实质允许只模拟辐射面间辐射（忽略向空