2. 中国航空工业集团公司 成都飞机设计研究所, 四川 成都 610091;
3. 南京航空航天大学 机械结构力学及控制国家重点实验室, 江苏 南京 210016
航天飞机轨道器在发射和再入阶段受到气动加热的作用[1-3], 为保证内部机体结构在可承受的温度范围内, 需要在机体表面附加热防护系统(TPS)[4-6], 即热防护系统起到阻隔外部热量进入到内部结构的作用, 其平面形状通常为矩形, 宽度大约为50~300 mm, TPS通常由涂层、防热瓦、及SIP构成。TPS最上面为高辐射率的涂层, 辐射率通常高于0.8, 其主要作用是将大部分气动热辐射到大气中, 对结构传热性能影响较小; 下层是防热瓦, 它是一种低导热系数的材料, 是TPS隔热的最主要部件, 其厚度一般在20~100 mm范围内; 防热瓦通过应变隔离垫与机体粘接在一起, 而应变隔离垫主要作用是保证防热瓦与机体之间的变形协调, 其厚度往往在1~8 mm之间, 且主要由nomex纤维构成。
TPS除了抵御外部气动加热外, 还受到外部气动力的作用, 在这些机械载荷作用下SIP将产生应力、应变响应, 若SIP的应力或应变超出其允许的最大值时SIP将发生破坏, 造成防热瓦与机体表面的分离与TPS的失效, 轨道器将发生灾难性的事故, 故TPS设计时必须进行SIP的力学响应分析, 其前提条件是需要获得SIP的力学性能参数。
目前应变隔离垫力学特性的研究还非常少, Sawyer等[7-8]通过试验研究了室温下SIP的力学特性, 获得了在拉伸和剪切载荷作用下SIP的应力-应变曲线和破坏应力, 取得了一些重要的研究成果。但是前人的研究环境为室温, 由于气动加热的作用SIP通常处于100~350℃温度环境下, 故需要研究高温对SIP力学特性的影响, 获得高温下SIP的力学性能参数, 并为后续SIP的力学分析提供输入条件, 这对SIP的强度以及TPS的完整性至关重要。
为了研究SIP在高温下的力学特性, 本文设计了SIP的平拉、平压和平剪试验件以及试验夹具, 在常温(23±2℃)和300℃高温环境下进行了拉伸、压缩和剪切试验, 获得了试验件的载荷-位移曲线、弹性模量及破坏应变, 并根据以上试验结果进行了高温和室温环境下SIP的应力分布规律研究。
1 SIP高温力学性能试验 1.1 试验件SIP是TPS的重要组成部分,如图 1所示。本文测试用的SIP是通过针刺技术制备而成的nomex织物, 为了研究其在高温下的基本力学性能, 本文设计了3类试验件, 分为拉伸、压缩和剪切试验件, 试验件尺寸和实物图分别如图 2和图 3所示。为了研究温度对SIP的力学特性的影响规律, 考虑了常温(23±2℃)和300℃高温2种试验环境。拉伸试验件上、下立方体铁块通过中间SIP粘接在一起; 压缩试验件仅为SIP; 而剪切试验件上、下铁板通过中间SIP粘接在一起, 构成单搭接的型式。图 2中的试验件尺寸L1=L2=60 mm, L3=120 mm及L4=60 mm, 此外SIP厚度H=2 mm。由此可知共有6组试验, 每组试验有5件试验件, 试验件共30件, 试验件编号规则为:T(C, S)-N(H)-1(2, 3, …, 5)。其中T代表拉伸试验, C代表压缩试验, S代表剪切试验; N代表常温环境, H代表 300℃高温环境; 1(2, 3, …, 5)为试验件号。
1.2 试验方法本文的试验环境分常温(23±2℃)和300℃高温两类, 试验采用MTS370.25单轴拉、压试验机。针对高温试验, 将试验机夹头和试验件置于高温环境箱中以达到加热试验件到指定温度的目的。此外高温试验还需要在试验结束后向环境箱注入液氮以达到快速冷却环境箱中的温度的目的。
根据拉伸试验标准ASTM Standard C297/C297M-04[9]、压缩试验标准ASTM Standard C365-03[10]和剪切试验标准ASTM Standard C273/C273M-07a[11], 拉伸、压缩和剪切试验件的加载情况如图 4所示, 通过试验机的测力传感器和位移传感器可以测得试验过程中载荷-位移曲线。
在拉伸和压缩试验过程中载荷P可表示为:
(1) |
式中, K为载荷-位移曲线的斜率, δ为位移传感器测试的位移, σ为SIP应力, S为SIP面积。又因为应力-应变关系:
(2) |
式中, E为SIP的拉伸或压缩弹性模量, ε为SIP应变, t0为SIP厚度。由(1)式和(2)式可得SIP的拉伸和压缩弹性模量的计算式:
(3) |
而在剪切试验过程中载荷P可表示为:
(4) |
式中, d为位移传感器测试的剪切试验件受载后的相对位移, τ为SIP剪应力。又因为应力-应变关系:
(5) |
式中, G为粘接层剪切弹性模量, γ为粘接层剪切应变。由(4)式和(5)式可得SIP的剪切弹性模量的计算式:
(6) |
当获得试验过程中的载荷-位移曲线的斜率K后, 可以通过(3)式和(6)式计算SIP的拉伸模量、压缩模量和剪切模量。由于拉伸和剪切试验中SIP会破坏, 而一般情况下压缩试验中SIP并不会破坏, 载荷会持续上升, 最后由夹具承载。故针对拉伸和剪切试验, 选取50%最大载荷与原点连线的斜率作为K值, 而针对压缩试验, 选取初始近似线性段的斜率最为K值。
2 试验结果 2.1 载荷-位移曲线首先考虑压缩试验, 高温下典型的载荷-位移曲线如图 5所示。
可以观察到载荷-位移曲线呈现2个阶段:第1阶段曲线的近似直线的斜率可视为SIP的刚度, 而随后曲线并没有出现载荷下降的现象, 且试验结束后进行SIP的检查, 并未发现明显的破坏。此外本文的SIP致密性较高, 而非传统的多孔介质材料, 故试验后SIP的体积变化也很小。综上所述可得试验过程中SIP并未发生明显破坏, 载荷持续上升, 最后主要由上、下夹具承受载荷, 即图中的第2阶段。
针对拉伸和剪切试验, 高温下典型的载荷-位移曲线如图 6所示, 同样可以观察到载荷-位移曲线呈现2个阶段, 第1阶段曲线为载荷上升阶段, 其近似直线的斜率可视为SIP的刚度, 而载荷-位移曲线因为试验件的破坏而造成的载荷下降现象, 即第2阶段的曲线。
2.2 高温下SIP的力学性能参数本文的拉伸、压缩和剪切试验测得了试验件的载荷-位移曲线, 将载荷-位移曲线的斜率K带入(3)式和(6)式即可获的SIP的拉伸、压缩和剪切模量, 此外还测得了SIP的破坏应变, 如表 1所示。从表 1中的结果数据可以得到以下结论:
试验件编号 | 弹性模量 | 破坏应变 | |||||
均值/MPa | 标准差/MPa | 变异系数 | 均值 | 标准差 | 变异系数 | ||
T-N | 4.96 | 0.93 | 0.188 | 0.69 | 0.04 | 0.058 | |
T-H | 4.48 | 0.54 | 0.121 | 0.61 | 0.06 | 0.098 | |
C-N | 2.49 | 0.16 | 0.064 | 无 | 无 | 无 | |
C-H | 2.26 | 0.24 | 0.106 | 无 | 无 | 无 | |
S-N | 2.29 | 0.29 | 0.127 | 1.64 | 0.08 | 0.049 | |
S-H | 2.09 | 0.11 | 0.053 | 1.43 | 0.10 | 0.070 |
1) 相比常温, 在300℃高温下SIP的弹性模量均有所下降, 其中SIP拉伸、压缩和剪切弹性模量分别下降了9.68%, 9.24%和8.73%, 即高温引起了SIP材料的软化, 导致了其刚度的下降;
2) 相比常温, 在300℃高温下SIP的拉伸和剪切破坏应变分别下降了11.59%和12.81%, 即高温导致了SIP破坏应变的下降;
3) SIP拉伸弹性模量明显大于其压缩和剪切弹性模量, 而压缩和剪切弹性模量却相近。此外SIP的剪切破坏应变明显大于其拉伸破坏应变。
SIP的以上试验结果可作为TPS强度设计时SIP的输入参数, 考虑高温下SIP刚度和强度的下降以及各向异性特征能获得更精确的SIP强度分析结果, 提高TPS的完整性。
3 SIP高温应力分布规律航天飞机轨道器再入过程中受到气动吸力和表面摩擦力的作用, 此外在进行TPS的安装时需要在防热瓦外表面施加压力以完成与机体表面的粘接过程。本文研究了TPS在外部气动吸力、安装压力以及剪切力作用下SIP的应力分布规律, 如图 7所示。
分析模型自上至下依次为隔热瓦、SIP和机体结构。其中TPS宽度W=150 mm, 防热瓦厚度H1=45 mm, SIP厚度H2=2 mm, 机体为铝合金材料, 厚度H3=1.5 mm, 弹性模量E1=72 GPa, 泊松比为ν1=0.3。防热瓦弹性模量E2=70 MPa, 泊松比为ν2=0.25, SIP弹性模量采用表 1中的试验结果。根据以上参数建立了二维有限元模型, 并且采用四边形平面应变单元进行模拟, 由于本文研究的是SIP的应力分布规律, 故假设气动吸力、安装压力以及剪切力均为0.01 MPa, 且机体底部固定。
分析获得了在气动吸力、安装压力和剪切力作用下的SIP正应力及剪切应力的分布情况, 如图 8所示, 从图中可获得以下结论:
1) 在局部吸力、安装压力和剪切力作用下, SIP应力呈现中间低边缘高的趋势, 峰值出现在在边缘处, 这是由于泊松效应造成的, 在粘接结构中广泛存在此现象;
2) 在局部吸力和安装压力作用下, SIP应力主要为沿厚度方向正应力, 而在剪切力作用下, SIP不仅产生剪切应力, 还会产生沿厚度方向正应力, 并且正应力一端受拉, 一端受压;
3) SIP边缘是最薄弱的地方, 最容易发生破坏, 应给予重点考察和足够的重视。
4 结论1) 在常温(23±2℃)和300℃高温2种试验环境下进行了SIP的拉伸、压缩和剪切试验, 获得了试验件的载荷-位移曲线、弹性模量及破坏应变。
2) 相比常温下, 在300℃高温下SIP的拉伸、压缩和剪切弹性模量分别下降了9.68%, 9.24%和8.73%, 在300℃高温下SIP的拉伸和剪切破坏应变分别下降了11.59%和12.81%, 即高温引起了SIP刚度和破坏应变的下降。
3) 在局部吸力、安装压力和剪切力作用下, SIP应力呈中间低边缘高的趋势, 即SIP的边缘处是高应力区域, 是危险区域。
[1] | Knight D, Chazot O, Austin J, et al. Assessment of Predictive Capabilities for Aerodynamic Heating in Hypersonic Flow[J]. Progress in Aerospace Sciences, 2017, 90: 39-53. DOI:10.1016/j.paerosci.2017.02.001 |
[2] | Li J, Chen H, Zhang S, et al. On the Response of Coaxial Surface Thermocouples for Transient Aerodynamic Heating Measurements[J]. Experimental Thermal and Fluid Science, 2017, 86: 141-148. DOI:10.1016/j.expthermflusci.2017.04.011 |
[3] | Vasil'evskii S A, Gordeev A N, Kolesnikov A F. Local Modeling of the Aerodynamic Heating of the Blunt Body Surface in Subsonic High-Enthalpy Air Flow. Theory and Experiment on a High-Frequency Plasmatron[J]. Fluid Dynamics, 2017, 52(1): 158-164. DOI:10.1134/S001546281701015X |
[4] | Huang J, Li P, Yao W X. Thermal Protection System Gap Analysis Using a Loosely Coupled Fluid-Structural Thermal Numerical Method[J]. Acta Astronautica, 2018, 146: 368-377. DOI:10.1016/j.actaastro.2018.02.047 |
[5] |
黄杰, 姚卫星, 陈炎, 等. 热防护系统分区协调耦合推进方法[J]. 宇航学报, 2018, 39(1): 27-34.
Huang Jie, Yao Weixing, Chen Yan, et al. Division Coordinating Coupled Marching Method on Thermal Protection System[J]. Journal of Astronautics, 2018, 39(1): 27-34. (in Chinese) |
[6] | Huang J, Yao W X, Li P. Uncertainty Dynamic Theoretical Analysis on Ceramic Thermal Protection System Using Perturbation Method[J]. Acta Astronautica, 2018, 148: 41-47. DOI:10.1016/j.actaastro.2018.04.027 |
[7] | Sawyer J W, Rummler D R. Room Temperature Mechanical Properties of Shuttle Thermal Protection System Materials[R]. NASA TM-81786, 1980 |
[8] | Sawyer J W, Waters W A. Room Temperature Shear Properties of the Strain Isolator Pad for the Shuttle Thermal Protection System[R]. NASA TM-81900, 1981 |
[9] | American Society for Testing and Materials, Standard Test Method for Flatwise Tensile Strength of Sandwich Constructions[S]. C297/C297M-04, 2004 |
[10] | American Society for Testing and Materials, Standard Test Method for Flatwise Compressive Properties of Sandwich Cores[S]. C365-03, 2003 |
[11] | American Society for Testing and Materials, Standard Test Method for Shear Properties of Sandwich Core Materials[S]. C273/C273M-07a, 2007 |
2. AVIC Chengdu Aircraft Design and Research Institute, Chengdu 610091, China;
3. State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China