论文部分内容阅读
随着技术的发展航天器有向着高热流密度发展的趋势,如何在高热流密度、热不平衡、热不均匀的状态下精确计算出航天器的温度场成为了一个不可回避的问题。本文首先通过对目前国内外的航天器的热控现状,包括热控措施、设计方法及热传导方程的求解方法的研究分析,最终选用精细积分方法来解决高热流密度下瞬态温度场计算的问题。精细积分方法首先被广泛应用于结构动力学方程,之后又发展到传热学中。在传热学中求解常微分方程通常用的是有限差分法,然而有限差分法会带来一定的计算误差和存在刚性问题,与之相比,精细积分方法可以达到计算机步长范围内的计算精度、不受计算步长的限制、无条件稳定且不会出现刚度问题。详细阐述精细积分的原理及其在传热方面的应用与改进,分别讨论了指数矩阵的算法、线性齐次和非齐次热传导方程、非线性热传导方程的精细积分;对于非齐次方程的求解,本文的思路是利用叠加原理,线性常微分方程的解可以表达成通解加特解,然后分别求通解和特解的精细积分。通解的积分可达到计算机步长范围内的精度,特解的精细积分是将积分域化成2N份并对其进行积分,这里的非齐次项是多项式形式,将积分求和转化成一个递推过程,只需要N次矩阵的乘法就完成了,这样避免了矩阵的求逆且具有和通解的精细积分具有同样的精度和高效计算效率;对于非线性热传导的情况,非线性项是温度的函数,温度又是时间的函数,本文采用拟合手段和利用预测-校正算法将非线性项转化成关于时间的多项式形式,然后用非齐次方程的特解精细积分来计算,同样具有高精度性和计算效率。对于各种离散格式,只要能生成其标准的格式,精细积分就可以与它相结合生成高精度的求解格式,之后通过多个算例来验证了其有效性及高精度性。但是到目前为止还没有将精细积分应用于航天器的热分析中,本文以美国Thermal Desktop软件为依托,作为前处理建模手段,利用其输出各个参数文件并基于上述精细积分原理,成功的设计出了一款精细积分瞬态温度场求解器,首次将精细积分算法应用于整星在轨运行的热分析中。通过算例证明了在高热流密度、热不平衡、热不均匀的情况下其高精度性和高效的计算效率,并应用该精细积分求解器对分布式SAR做了整星的在轨运行瞬态热分析。当然,该精细积分方法可以和其他有限元、有限差分等软件相结合,构成相应的求解器,这些工作有待于下一步研究。