張瑾瑜，雷 豹，楊東生，許俊偉，嚴(yán) 旭

(中國運(yùn)載火箭技術(shù)研究院，北京 100076)

0 引言

航天飛行器高溫透波結(jié)構(gòu)是在惡劣使用環(huán)境條件下保護(hù)航天飛行器的通信、遙測(cè)、制導(dǎo)、引爆等系統(tǒng)天線正常工作的一種結(jié)構(gòu)，如天線罩、天線窗等，其材料常選用石英玻璃、石英陶瓷、石英纖維增強(qiáng)二氧化硅、氮化硅、磷酸鹽復(fù)合材料等耐高溫多功能電介質(zhì)材料[1-6]。隨著航天飛行器的不斷發(fā)展，飛行速度不斷提高，面臨的熱環(huán)境也越來越惡劣，而高溫透波結(jié)構(gòu)與其安裝結(jié)構(gòu)的材料熱膨脹系數(shù)存在一定差異，在惡劣熱環(huán)境條件下，它們之間可能產(chǎn)生較大的熱變形，進(jìn)而可能導(dǎo)致高溫透波結(jié)構(gòu)在熱應(yīng)力作用下破壞失效。因此，研究高溫透波結(jié)構(gòu)的熱匹配性能，對(duì)航天飛行器的發(fā)展具有重要意義[7]。

本文針對(duì)一種航天飛行器典型天線窗類高溫透波結(jié)構(gòu)及其安裝形式，選取典型熱環(huán)境條件，采用有限元仿真分析方法，研究高溫透波結(jié)構(gòu)在不同材料、不同尺寸條件下的熱匹配性能。

1 結(jié)構(gòu)及連接形式

航天飛行器典型天線窗及其安裝結(jié)構(gòu)如圖1所示，天線窗通過鋁合金底座和鋁合金壓環(huán)固定連接后，貼壁安裝在航天飛行器金屬殼體內(nèi)壁面上。

圖1 天線窗及其安裝結(jié)構(gòu)示意圖Fig.1 Schematic diagram of antenna window and its installation structure

表1為天線窗結(jié)構(gòu)選用石英玻璃、石英陶瓷、石英纖維增強(qiáng)二氧化硅、氮化硅等不同材料的主要性能參數(shù)[8-11]。

表1 材料性能參數(shù)

2 熱應(yīng)力數(shù)學(xué)模型

結(jié)構(gòu)受熱產(chǎn)生的熱變形、熱應(yīng)力等熱匹配性能主要與材料熱膨脹系數(shù)、彈性模量及溫度變化量有關(guān)[12]。

熱變形量的計(jì)算公式為

ΔL=λLΔt

(1)

式中, ΔL代表結(jié)構(gòu)受熱產(chǎn)生的變形量，λ代表熱膨脹系數(shù)，L代表結(jié)構(gòu)線性尺寸，Δt代表溫度變化量。

假設(shè)兩種材料的溫升相同，則熱應(yīng)力計(jì)算公式為

σ=EcΔt(λc-λs)

(2)

式中,σ代表結(jié)構(gòu)受熱產(chǎn)生的熱應(yīng)力，Ec代表材料彈性模量，Δt代表溫度變化量，λc和λs代表熱匹配結(jié)構(gòu)的材料熱膨脹系數(shù)。

本文主要針對(duì)高溫透波天線窗結(jié)構(gòu)的材料性能、結(jié)構(gòu)尺寸研究其熱匹配性能。

3 有限元分析

3.1 有限元模型建立

1)模型假設(shè)。假設(shè)天線窗與其安裝結(jié)構(gòu)之間無間隙。取非金屬材料表面黑度0.8，金屬材料表面黑度0.3。根據(jù)項(xiàng)目要求，初溫設(shè)置為60℃。考慮熱輻射，波爾茲曼常數(shù)5.67×10-8W/(m2·K4)。針對(duì)金屬殼體等塑性材料，不考慮對(duì)結(jié)果影響較小的倒角、圓角等結(jié)構(gòu)特征。

2)幾何建模。用CATIA軟件建立天線窗及其安裝結(jié)構(gòu)的幾何模型并導(dǎo)入Abaqus，如圖2所示。

圖2 天線窗及其安裝結(jié)構(gòu)的三維模型Fig.2 3D model of antenna window and its installation structure

3)網(wǎng)格劃分。采用力熱耦合四面體單元C3D8T和Abaqus內(nèi)置的力熱耦合分析方法進(jìn)行分析計(jì)算，如圖3所示，共生成94 663個(gè)節(jié)點(diǎn)和67 082個(gè)單元。

圖3 天線窗及其安裝結(jié)構(gòu)的有限元網(wǎng)格Fig.3 Finite element mesh of antenna window and its installation structure

4)材料性能參數(shù)。天線窗結(jié)構(gòu)貼壁安裝在航天飛行器金屬殼體內(nèi)壁面上并與防熱層接觸，金屬殼體材料為ZL114A，防熱層為低密度石英酚醛，天線窗材料性能參數(shù)見表2。

表2 材料性能參數(shù)

5)載荷數(shù)據(jù)。選用典型航天行器熱環(huán)境條件進(jìn)行天線窗結(jié)構(gòu)熱匹配分析，如圖4所示。冷壁熱流最高300 kW/m2，冷壁熱流總加熱量30 000 kJ/m2，恢復(fù)焓最高1 800 kJ/ m2。

圖4 典型天線窗熱環(huán)境條件Fig.4 Typical thermal load curves of antenna window

6)接觸定義及設(shè)置。本文采用的有限元模型為軸對(duì)稱模型，如圖5所示，防熱層與金屬殼體、金屬殼體與壓環(huán)、壓環(huán)與底座之間均定義綁定約束，其余定義接觸，并設(shè)置對(duì)稱邊界，底端固支。

圖5 有限元模型示意圖Fig.5 Finite element analytical model

7)防熱計(jì)算。初溫設(shè)置為60℃，材料表面輻射系數(shù)0.8。提交計(jì)算后，獲得天線窗及其安裝結(jié)構(gòu)的溫度場(chǎng)分布。

8)熱應(yīng)力計(jì)算。將天線窗及其安裝結(jié)構(gòu)的溫度場(chǎng)分布導(dǎo)入計(jì)算模型，作為熱應(yīng)力計(jì)算的輸入。在材料性能參數(shù)中設(shè)置強(qiáng)度數(shù)據(jù)，提交計(jì)算后，獲得天線窗及其安裝結(jié)構(gòu)的熱應(yīng)力分布。

3.2 結(jié)果分析

以石英纖維增強(qiáng)二氧化硅天線窗為例，對(duì)有限元計(jì)算結(jié)果進(jìn)行分析，計(jì)算結(jié)果表明：天線窗熱變形及熱應(yīng)力變化規(guī)律一致，均在末秒達(dá)到最大值，最大變形量為0.3 mm，最大拉應(yīng)力為17.9 MPa，最大壓應(yīng)力為15.2 MPa，位于天線窗側(cè)壁，如圖6和圖7所示。

圖6 天線窗位移分布圖Fig.6 Displacement distribution of antenna window

(a)天線窗最大拉應(yīng)力分布圖

(b)天線窗最大壓應(yīng)力分布圖圖7 天線窗應(yīng)力分布圖Fig.7 Stress distribution of antenna window

對(duì)天線窗、防熱層及金屬殼體溫度隨時(shí)間變化的關(guān)系曲線進(jìn)行分析，如圖8所示。從圖中可以看出：防熱層與天線窗溫度隨時(shí)間變化規(guī)律基本一致，且溫度差值??；金屬殼體與天線窗溫度隨時(shí)間變化規(guī)律差別大，且溫度差值大；金屬殼體末秒溫度達(dá)最高，此時(shí)天線窗熱變形、熱應(yīng)力也達(dá)到最大值。

圖8 溫度隨時(shí)間變化曲線圖Fig.8 Temperature curves with time

分析認(rèn)為，天線窗熱匹配產(chǎn)生的熱變形、熱應(yīng)力主要與金屬殼體受熱膨脹有關(guān)。金屬殼體在末秒溫度達(dá)到最高，由式(1)可知，此時(shí)金屬殼體熱膨脹量也達(dá)最大值，因此,天線窗熱匹配產(chǎn)生的熱變形、熱應(yīng)力在末秒達(dá)到最大值。

3.3 材料性能對(duì)熱匹配的影響

由上文分析可知，天線窗結(jié)構(gòu)熱匹配產(chǎn)生的熱變形與熱應(yīng)力的變化規(guī)律一致，下文主要通過熱應(yīng)力分析材料性能對(duì)熱匹配的影響。在相同熱環(huán)境條件下，保持天線窗結(jié)構(gòu)尺寸參數(shù)不變，對(duì)比分析了選用不同高溫透波材料時(shí)結(jié)構(gòu)熱匹配產(chǎn)生的熱應(yīng)力，仿真計(jì)算結(jié)果得到熱應(yīng)力排序?yàn)椋菏⑻沾?石英玻璃>氮化硅>磷酸鹽復(fù)合材料>石英纖維增強(qiáng)二氧化硅,如表3所示。

表3 不同材料天線窗熱應(yīng)力計(jì)算結(jié)果

分析認(rèn)為,材料彈性模量是影響天線窗熱匹配性能的主要因素。不同材料天線窗結(jié)構(gòu)熱匹配產(chǎn)生的熱應(yīng)力與材料彈性模量變化規(guī)律一致，如圖9所示。由于5種材料的熱膨脹系數(shù)相近，均遠(yuǎn)小于金屬殼體熱膨脹系數(shù)，且熱物理性能參數(shù)相差不大，溫升變化量相當(dāng)，如圖10所示。由式(2)可知，材料彈性模量成為影響結(jié)構(gòu)熱應(yīng)力的主要因素，仿真分析計(jì)算結(jié)果與理論分析一致。

圖9 不同材料熱應(yīng)力及彈性模量變化曲線圖Fig.9 Thermal stress and modulus of elasticity ofdifferent materials

圖10 不同材料天線窗表面溫度隨時(shí)間變化關(guān)系曲線圖Fig.10 Temperature curves of antenna window’s surface with time for different materials

結(jié)合材料抗彎強(qiáng)度，進(jìn)行不同材料天線窗剩余強(qiáng)度系數(shù)分析，剩余強(qiáng)度系數(shù)計(jì)算公式如下

η=σb/σ

(3)

式中，η代表剩余強(qiáng)度系數(shù)，σb代表材料強(qiáng)度極限，σ代表結(jié)構(gòu)產(chǎn)生的最大應(yīng)力。已知天線窗熱匹配最大熱應(yīng)力為拉應(yīng)力，陶瓷材料的抗拉強(qiáng)度為抗彎強(qiáng)度的80%～97%[13]，剩余強(qiáng)度系數(shù)計(jì)算結(jié)果如表4所示。由表4可知，石英纖維增強(qiáng)二氧化硅天線窗結(jié)構(gòu)熱匹配產(chǎn)生的熱應(yīng)力剩余強(qiáng)度系數(shù)大，材料具有優(yōu)越的熱匹配性能。

表4 不同材料天線窗剩余強(qiáng)度系數(shù)

3.4 結(jié)構(gòu)尺寸對(duì)熱匹配的影響

在相同熱環(huán)境、相同材料條件下，對(duì)比分析不同材料尺寸條件下天線窗結(jié)構(gòu)熱匹配產(chǎn)生的熱應(yīng)力，天線窗尺寸組合如表5所示。

表5 天線窗尺寸組合

仿真計(jì)算結(jié)果如表6和表7所示。厚度3 mm氮化硅天線窗，直徑75 mm時(shí)最大拉應(yīng)力為48 MPa，直徑100 mm時(shí)最大拉應(yīng)力為132 MPa；直徑75 mm石英陶瓷天線窗，厚度3 mm時(shí)最大拉應(yīng)力為134 MPa，厚度12 mm時(shí)最大拉應(yīng)力為68 MPa。

表6 不同直徑天線窗熱應(yīng)力計(jì)算結(jié)果(3 mm氮化硅)

表7 不同厚度天線窗熱應(yīng)力計(jì)算結(jié)果(φ75 mm石英陶瓷)

對(duì)計(jì)算結(jié)果進(jìn)行分析后認(rèn)為，天線窗直徑越小、厚度越大，熱應(yīng)力越小，結(jié)構(gòu)熱匹配性能越好，對(duì)實(shí)際工程應(yīng)用具有借鑒意義。

4 結(jié)論

本文針對(duì)一種航天飛行器典型天線窗類高溫透波結(jié)構(gòu)及其安裝形式，通過有限元仿真分析方法，分析了其熱匹配性能及在不同材料、不同尺寸條件下的熱匹配性能，為其在航天飛行器上的應(yīng)用提供了參考，具體結(jié)論如下：

1)天線窗結(jié)構(gòu)熱匹配產(chǎn)生的熱變形、熱應(yīng)力主要與航天飛行器金屬殼體受熱膨脹有關(guān)，與防熱層關(guān)系較小，當(dāng)金屬殼體溫度最高時(shí)天線窗結(jié)構(gòu)熱應(yīng)力達(dá)最大值；

2)在相同熱環(huán)境條件下，不同材料的天線窗熱應(yīng)力排序?yàn)椋菏⑻沾?石英玻璃>氮化硅>磷酸鹽>石英纖維增強(qiáng)二氧化硅，與材料的彈性模量變化規(guī)律一致，結(jié)合材料抗彎強(qiáng)度，分析認(rèn)為石英纖維增強(qiáng)二氧化硅天線窗具有更優(yōu)越的熱匹配性能；

3)天線窗直徑越小，厚度尺寸越大，熱應(yīng)力越小，結(jié)構(gòu)熱匹配性能越好，對(duì)實(shí)際工程應(yīng)用具有借鑒意義。