YZX有限元方法及CAE软件软件有限元方法及机械电子及信息系统研究所耿海鹏haipenggeng@126.comANSYS/MultiphysicsANSYS/LS-DYNA?LS-ANSYS/Emag?ANSYS/FLOTRAN?ANSYS/Mechanical?由NordriDesign提供www.nordridesign.comANSYS/Thermal?ANSYS/Structural?ANSYS/LinearPlus?热辐射的定义热辐射是一种通过电磁波传递热能的方式。电磁波以光的速度进行传递，而能量传递与辐射物体之间的介质无关。热辐射只在电磁波的频谱中占小部分的带宽。由于辐射产生的热流与物体表面的绝对温度的四次方成正比，因此热辐射有限元分析是高度非线性的。物体表面的辐射遵循Stefan-Boltzmann定律：q/A=σT4式中：T—物体表面的绝对温度；σ—Stefan-Boltzmann常数，英制为0.119×10-10BTU/hr-in-R，公制为5.67×10-8基本概念下面是对辐射分析中用到的一些术语的定义：黑体黑体被定义为在任意温度下，吸收并发射最大的辐射能的物体；通常的物体为“灰体”，即ε<1；在某些情况下，辐射率（黑度）随温度变化；辐射率ε黑度）辐射率（黑度）物体表面的辐射率（黑度）定义为物体表面辐射的热量与黑体在同一表面辐射热量之比。ε=q/qb式中：ε－辐射率（黑度）q－物体表面辐射热量qb－黑体在同一表面辐射热量基本概念(conti)基本概念?形状系数形状系数用于计算两个面之间的辐射热交换，在ANSYS中，可以用隐藏/非隐藏的方法计算2维和三维问题，或者用半立方的方法来计算3维问题。表面I与表面J之间的形状系数为：形状系数是关于表面面积、面的取向及面间距离的函数；由于能量守恒，所以：根据相互原理由辐射矩阵计算的形状系数为：式中:θ－单元法向与单元I,J连线的角度r－单元I,J重心的距离有限单元模型的表面被处理为单元面积dAI及dAJ，然后进行数字积分。基本概念(conti)基本概念辐射对在辐射问题中，辐射对由一些相互之间存在辐射的面组成，可以是开放的或是闭合的。在ANSYS中，可以定义多个辐射对，它们相互之间也可以存在辐射ANSYS使用辐射对来计算一个辐射对中各面间的形状系数；每一个开放的辐射对都可以定义自己的环境温度，或是向周围环境辐射的空间节点。Radiosity求解器当所有面上的温度已知时，Radiosity求解器方法通过计算每一个面上的辐射热流来得到辐射体之间的热交换。而面上的热流为接下来的热传导分析提供了有限元模型的边界条件。重复上面的过程，就会由于新的时间步或者新的迭代循环会得到新的热流边界条件，从而计算出新的温度分布。在计算中使用的每个表面的温度必须是均匀的，这样才能满足辐射模型的条件。分析热辐射问题针对不同的情况ANSYS为热辐射分析提供了四种方法。热辐射线单元（LINK31），），模拟两节点间（或多对节点）间辐射；热辐射线单元（LINK31），表面效应单元（SURF151及SURF152），模拟点对面（线）的辐射；表面效应单元利用AUX12生成辐射矩阵，模拟更一般的面与面（或线与线）的辐射(只有ANSYS/Multiphysics,ANSYS/Mechanical和ANSYS/Professional这些产品提供辐射矩阵生成器)；Radiosity求解器方法，求解二维、三维面与面之间的热辐射，该方法对所有含温度自由度的二维和三维单元都适用。(只有Multiphysics，ANSYS/Mechanical和ANSYS/Professional这些产品提供Radiosity求解器)可以将上面四种辐射方法中的任何一种用于稳态或瞬态热分析中。辐射是一种非线性现象，因此需要进行平衡迭代来得到收敛解。求解方法1：应用节点间单元的热辐射问题求解方法：应用节点间单元的求解热辐射热辐射非线性线单元LINK31用于计算两节点间或多对节点间的简单辐射热传递，节点的位置是任意的，可作为其它单元的节点。LINK31需要定义如下数据：材料属性：EMIS辐射率（可以随温度变化）实常数：AREA（Ai）（有效辐射面积）FORMF（Fij）（形状系数）SBCONST（Stefan-Boltzman常数）有关LINK31的使用实例，请参考《ANSYS校验手册》：VM106RadiantenergyemissionVM106VM107ThermocoupleradiationVM107黑体辐射如图所示为一黑体,表面积为1m2,形状系数和辐射率均为1,温度T为2000℃，环境温度Ta为0°C,求黑体的辐射热流率