ANSYS热辐射解析.docx

第六章热辐射分析6・1热辐射的定义热辐射是一种通过电磁波传递热能的方式电磁波以光的速度进行传递，而 能量传递与辐射物体之间的介质无关热辐射只在电磁波的频谱中占小部分的带 宽由于辐射产生的热流与物体表面的绝对温度的四次方成正比，因此热辐射有 限元分析是高度非线性的物体表面的辐射遵循Stefan-Boltzmann定律：式中：于一物体表面的绝对温度；—Stefan-Boltzmann 常数，英制为 O.119x1O-ioBTU/hr-in-R，公制为 5.67x10-8 6.2基本概念下面是对辐射分析中用到的一些术语的定义：黑体黑体被定义为在任意温度下，吸收并发射最大的辐射能的物体；通常的物体为“灰体”，即区1；在某些情况下，辐射率（黑度）随温度变化；辐射率（黑度）物体表面的辐射率（黑度）定义为物体表面辐射的热量与黑体在同一表面辐 射热量之比式中：£—辐射率（黑度）—物体表面辐射热量芻一黑体在同一表面辐射热量形状系数形状系数用于计算两个面之间的辐射热交换，在ANSYS中，可以用隐藏/非 隐藏的方法计算2维和三维问题，或者用半立方的方法来计算3维问题表面I与表面J之间的形状系数为：F 一表面J接受的由表面I发出的辐射能D 由表面I发出的辐肘能形状系数是关于表面面积、面的取向及面间距离的函数；由于能量守恒，所以：FI1 + FI2 + ...+r11 = 1根据相互原理：直局-血局由辐射矩阵计算的形状系数为：式中：9—单元法向与单元I,J连线的角度f —单元IJ重心的距离有限单元模型的表面被处理为单元面积dA,及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校验手册》：VM106Radiant energy emissionVM107Thermocouple radiation6・5点与面间的热辐射应用表面效应单元可以方便地计算点与面间的辐射，包括2D的SURF151及 3D 的 SURF152 单元：首先在实体单元的辐射表面覆盖一层表面效应单元；单元关键选项KEY0PT(9)激活这些单元的热辐射分析功能；如果设置KEY0PT(9) = 1 (缺省值)，则可在单元实常数中定义形状系数；如果设置KEYOPT(9)=2或3,则程序基于单元表面的法向与附加节点的位置 关系、考虑余弦效应，计算形状系数。

使用表面单元进行热辐射分析的实例，请参考《ANSYS校验手册》：VM192Cooling of a billet by radiation6・6AUX12-辐射矩阵生成器只有 ANSYS/Multiphysics, ANSYS/Mechanical, ANSYS/Thermal 等产品提供 AUX12, 用于计算多个辐射面之间的辐射热传递AUX12生成辐射面间的形状系数矩阵， 并将此矩阵作为超单元用于热分析分析模型可包含隐藏面或部分隐藏面，以及 能吸收辐射能量的“空间节点”使用AUX12计算面与面之间的辐射可分为以下三 个步骤：定义辐射面生成辐射矩阵将辐射矩阵用于热分析6.6.1定义辐射面分析类型庄阳 iited tor danty：-图I二维和三维的辐射面 i® ◎辐射百上覆盖的单元在2D有限元模型的辐射边上覆盖一层LINK32单元，或在3D有限元模型的辐 射面上覆盖一层SHELL57单元，步骤如下：在前处理中创建热分析模型由于辐射表面不支持对称条件，包含辐射的模 型就无法利用几何上的对称性，因此必须建立完整的分析模型二维和三维的辐 射面示意图6-1如下：在辐射面上覆盖一层SHELL57（3D）或LINK32（2D）单元，如图6-2所示。

最好的方法是先选择辐射面的节点，然后用下面的方法创建面上的单元：命令：ESURFGUI:Main Menu>Preprocessor>Create>Elements>Surf Effect>Extra Node Main Menu>Preprocessor>Create>Elements>Surf Effect>No extra Node在执行上述命令之前，要确定首先激活了相应的单元类型及材料属性，如果 不同辐射面的辐射率不同，建议用材料编号来区分这些面还要注意SHELL57 或LINK32的节点一定要与相应实体单元对应节点编号重合否则计算结果将是 错误的所覆盖的SHELL57或LINK32单元的方向是非常重要的AUX12假设辐射方向 是SHELL57单元坐标系的正Z向或LINK32单元坐标系的正Y向，因此必须正确 地划分覆盖层的网格以使辐射面相对单元的方向是由节点的排列顺序决定的， 如下图所示：图6-3覆盖单元的方向Viewing 耳 wao ar：nQSbrf^i显示单元辐射方向的方法如下：命令：/PSYMB,ESYS,1GUI: Utility Menu>PlotCtrls>Symbols, 将 ESYS Element Coordinate 设置为 ON。

定义一个空间节点，用于吸收没有被模型中其它辐射面吸收的辐射能量这 个节点的位置是任意的，对于一个开放系统通常需要空间节点，而对于封闭系统 则不得设置空间节点6.6.2 生成辐射矩阵计算辐射矩阵可按如下的步骤完成：1．进入 Aux12命令： /AUX12GUI： Main Menu>Radiation2. 选择构成辐射面的节点和单元较简便的方法是根据单元属性选择（如单 元类型）选择单元，然后选择所有 Attached to 单元的节点（同时应注意将空间节 点也选择进来）:命令：ESEL,S,TYPE 和 NSELGII:Utility Menu>Select>Entities3.确定所分析的模型是3D还是2D命令： GEOMGUI:Main Menu>Radiation Matrix>Other SettingsAUX12采用不同的算法分别计算2D和3D模型的形状系数,AUX12默认为3D 2D可以是平面的（NDIV=O），也可以是轴对称的（NDIV >0），缺省为平面的 轴对称模型在内部展成3D, NDIV是圆周方向分割数量例如NDIV=10，则每段 为 36 度4．确定辐射率（缺省为 1.0）：命令： EMISGUI:Main Menu>Radiation Matrix>Emmisivities5. 定义 Stefan-Baltzmann 常数（缺省情况下，英制单位为 O.119x1O-io Btu/hr-in2-R4, 国际单位制为5.67x10-8W/m2K4）。

命令： STEFGUI： Main Menu>Radiation Matrix>Other Settings6. 确定计算形状系数的方法命令： VTYPEGUI:Main Menu>Radiation Matrix>Write Matrix可选择非隐藏或隐藏两种方法之一，非隐藏方法计算每个单元对其它单元的 形状系数，无论两单元之间是否有阻碍；隐藏方法（默认）首先用“隐藏线”算法 确定两单元之间是否“可见”，如果目标单元与辐射单元的辐射方向指向对方，而 且设有其它单元阻碍，则它们是“可见”的，形状系数按如下方法计算：每一个辐射单元被封闭成一个半径为单位值的半球（3D）或半圆（2D）；所有的目标单元向这个半球或半圆投影；一定数量（默认为20）的射线由辐射单元面投向半球或半圆这样，形状系 数就是投到投影面上的射线数量与辐射面发出的射线的数量之比，通常设定的射 线数量越多，形状系数的精度越高可以通过设定VTYPE命令的变量NZONE或 上述的菜单来设定射线数量7．如果有必要（例如开放系统），应指定空间节点：命令：SPACEGUI:Main menu>Radition>Matrix>Other Settings8•将辐射矩阵写到文件jobname.sub中如果想要写更多的辐射矩阵，为不同 的矩阵指定不同的文件名：命令： WRITEGUI: Main Menu>Radiation Matri仪〉Write Matrix如果需要打印出辐射矩阵，应在执行Write命令之前执行命令：mprint,1。

9.选择所有节点及单元命令： ALLSELGUI:Utility Menu>Select>Everything现在就已经将辐射矩阵作为一个超单元写入到一个文件中了6.6.3 使用辐射矩阵进行热分析生成了辐射矩阵之后，重新进入前处理器，定义作为超单元的辐射矩阵步 骤如下：1. 重新进入前处理器，选择单元MATRIX50 （超单元），并设置单元Keyoption 为热辐射分析命令： /PREP7GUI： Main Menu>Preprocessor2. 设置缺省单元类型为超单元命令： TYPEGUI:Main Menu>Preprocessor>Create>Element>Elem Attributes3. 读入辐射超单元矩阵命令： SEGUI:Main Menu>Reprocessor>Create>Elements-Super elements-From .SUB4. 不选择或删除用于生成辐射矩阵的SHELL57或LINK32单元,因为在热分析 中已经不再需要了命令： EDELEGUI:Main Menu>Prepr。