冷壁-热壁热流的结构瞬态温度场分析

开发背景 

飞行器在大气层内做超声速或高超声速飞行时,由于飞行器与环境气体之间存在巨大动能差以及气体黏性的作用,环境气体会对飞行器表面产生剧烈的加热作用,这种现象称为气动加热。以6马赫飞行的高超声速飞行器的前端天线罩锥部的瞬时热流密度可高达1.2MW/㎡,驻点温度将超过1200°C。航天飞机穿越大气层时其机体、机翼、垂尾等大部分区域的温度在750°C~1450°C之间,前锥端部和进气道等部位甚至会出现接近1800°C的驻点温度。如此极端恶劣的高温热环境条件,使得高声速飞行器材料和结构的热防护和热强度问题成为事关研制成败的关键问题。

飞行器设计初期进行气动热环境评估中,通常以“冷壁热流”载荷作为设计输入。冷壁热流假设结构表面在受气动加热时表面温度恒定于热力学温度零度不变,它剥离了材料响应及其与环境的耦合作用,单纯依靠气动热计算得出。

热壁热流是考虑了结构热响应的结构净吸收热流,通过气动加热(冷壁热流、压强、恢复焓等)和结构瞬态温度场实时耦合计算得出,其中外壁表面温度会随着加热时间的改变而不断变化,形成所谓“热壁”,而这种不断变化的冷壁边界条件以及复杂耦合关系,使得热壁热流的计算变得很困难。

 

开发过程 

为解决飞行器的气动热环境仿真和热防护方案评估,优解未来团队在自主研发的新一代国产工业软件OptFuture上线了冷壁-热壁热流高效实时仿真功能,对标某国外商业CAE软件计算结果,OptFuture计算结果与“预期值”基本完全吻合,足可验证冷壁-热壁热流功能的正确性和可靠性。

 

验证案例 

分别对20mm x 20mm x 2mm和20mm x 20mm x 5mm的石英盖板进行冷壁热流分析。石英盖板初始温度25℃,环境温度23℃,冷壁面同时存在对环境辐射和冷壁热流边界,其中冷壁表面发射率为0.4,如图1。冷壁压力、焓值和热流密度随时间变化如图2所示。

图1 石英盖板一侧表面的冷壁热流和环境辐射示意

图2 压力、焓值和热流密度随温度的变化曲线
 
3.1 温度云图计算结果

在OptFuture软件进行瞬态热仿真分析,求解得到的2mm石英盖板在500秒的温度云图(沿着板厚方向)如图3所示,分析图1可知,在面内材料和边界条件都是均匀的,在厚度方向受冷壁热流和环境辐射影响,因此沿着板厚方向将出现温度梯度,图3的结果刚好符合该问题的物理规律。图4给出了5mm石英盖板在500秒的温度云图,同样沿着板厚方向,规律与2mm盖板类似。

图3 2mm薄板温度分布

 

图4 5mm薄板温度分布

 

3.2 与国外商业软件对比

将国外商业CAE软件计算所得冷壁中心位置温度记为 “ 国外商软T1”,热壁面中心位置温度记为 “ 国外商软T2”,将OptFuture计算所得冷壁中心位置温度记为“OptFuture T1”,热壁面中心位置温度记为“OptFuture T2”,图5给出了2mm石英盖板的温升曲线,图6给出了5mm石英盖板的温升曲线。从图5和图6可以看出,OptFuture计算所得温升曲线与国外商业CAE软件计算结果高度吻合。

图5 2mm薄板冷壁热流计算结果对比

图6 5mm薄板冷壁热流计算结果对比

结论 

以国外商业CAE软件冷壁热流程序计算结果作为对标值,OptFuture软件的计算结果误差小于0.1%,验证了OptFuture软件所开发的冷壁热流功能的正确性与可靠性
error: