利用DMAP重用热分析中的辐射交换系数矩阵.docxVIP

PAGE 1 PAGE 1 利用DMAP重用热分析中的辐射交换系数矩阵  在卫星热设计、热分析中，有大量的计算时间用在辐射交换系数矩阵的重复计算上。本文分析了Nastran中辐射交换系数矩阵的计算过程，并给出了重用的DMAP代码，大大节约了辐射热分析的计算时间。  1引言  在卫星、星载设备的热分析中，往往同样一套模型，需要计算多个状态。有瞬态热分析，也有稳态热分析。卫星在境内、境外的工作模式不同导致热功耗不同，是不同的状态。在热控设计阶段，不同的热管设定温度对星载设备的影响，也需要进行多个状态的分析。此外，利用热平衡试验结果，对热分析模型中的热阻等重要参数进行修正，也需要进行多轮计算。  对于有较多表面参与辐射热交换的复杂结构，每次计算首先都要花费大量的时间在辐射交换系数矩阵的计算上。但是，辐射交换系数矩阵只与辐射表面单元的数量、单元大小、单元间的距离、角度有关，对于许多只是修改了热流、温度边界条件的工况，或者只是修改热阻、材料热力学参数的工况，大量的计算时间花在重复计算辐射交换系数矩阵上了。假如可以重复利用第一轮计算中产生的辐射交换系数矩阵，将大大节约后续多轮分析的计算时间。  2辐射交换系数分析  2.1辐射视角系数 图1视角系数的定义 辐射视角系数定义如图1，从表面1向表明2传递的辐射热流密度q12为  其中T1和T2分别为表面1和表面2的温度，ε，σ分别为Stefan-Boltzmann常数和辐射效率。f12为视角系数，积分如下：  2.2.f06文件、pch文件中的辐射交换系数  在Nastran的bdf输入文件中，RADCAV卡片的PRTPCH域填写0，将在f06文件中打印视角系数计算的相关结果：  ***VIEWFACTORMODULE***OUTPUTDATA***CAVITYID=1***  ELEMENTTOELEMENTVIEWFACTORSC*PARTIAL  SURF-ISURF-JAREA-IAI*FIJFIJERRORSHADINGERRORSCALE  其中FIJ列的计算结果，即公式(1.2)中的f（可以用简洁的有解析解的例子进行确认）。在RADCAV卡片的PRTPCH域填写0后，同时还会产生一个pch文件。该文件是个文本文件，内容是RADMTX和RADLST卡片。通过查阅QuickReferenceGuide，RADMTX卡片定义的就是辐射交换系数矩阵。与f06对比可知，RADMTX卡片中的FIJ域对应的是f06文件中的AI*FIJ列，即辐射交换系数。  F06文件中的FIJ为辐射视角系数fij；RADMTX中的FIJ为辐射交换系数Ai*Fij。  在卫星热设计、热分析中，有大量的计算时间用在辐射交换系数矩阵的重复计算上。本文分析了Nastran中辐射交换系数矩阵的计算过程，并给出了重用的DMAP代码，大大节约了辐射热分析的计算时间。  1引言  在卫星、星载设备的热分析中，往往同样一套模型，需要计算多个状态。有瞬态热分析，也有稳态热分析。卫星在境内、境外的工作模式不同导致热功耗不同，是不同的状态。在热控设计阶段，不同的热管设定温度对星载设备的影响，也需要进行多个状态的分析。此外，利用热平衡试验结果，对热分析模型中的热阻等重要参数进行修正，也需要进行多轮计算。  对于有较多表面参与辐射热交换的复杂结构，每次计算首先都要花费大量的时间在辐射交换系数矩阵的计算上。但是，辐射交换系数矩阵只与辐射表面单元的数量、单元大小、单元间的距离、角度有关，对于许多只是修改了热流、温度边界条件的工况，或者只是修改热阻、材料热力学参数的工况，大量的计算时间花在重复计算辐射交换系数矩阵上了。假如可以重复利用第一轮计算中产生的辐射交换系数矩阵，将大大节约后续多轮分析的计算时间。  2辐射交换系数分析  2.1辐射视角系数 图1视角系数的定义 辐射视角系数定义如图1，从表面1向表明2传递的辐射热流密度q12为  其中T1和T2分别为表面1和表面2的温度，ε，σ分别为Stefan-Boltzmann常数和辐射效率。f12为视角系数，积分如下：  2.2.f06文件、pch文件中的辐射交换系数  在Nastran的bdf输入文件中，RADCAV卡片的PRTPCH域填写0，将在f06文件中打印视角系数计算的相关结果：  ***VIEWFACTORMODULE***OUTPUTDATA***CAVITYID=1***  ELEMENTTOELEMENTVIEWFACTORSC*PARTIAL