不同ZrC含量的(C/C)/SiC-ZrC复合材料的抗烧蚀性能

王玲玲; 闫联生; 郭春园; 张宏亮; 王坤杰

doi:10.13801/j.cnki.fhclxb.20200110.003

不同ZrC含量的(C/C)/SiC-ZrC复合材料的抗烧蚀性能

1.
西安航天复合材料研究所，西安 710025
2.
空军装备部驻西安第八军事代表室，西安 710065

基金项目: 国家自然科学基金(21676163)

详细信息

通讯作者:
王玲玲，硕士，高级工程师，研究方向为碳纤维增强陶瓷基复合材料性能　E-mail：1985wllzpz@sina.com

中图分类号: TB332
计量
- 文章访问数: 1229
- HTML全文浏览量: 608
- PDF下载量: 82
出版历程
- 收稿日期: 2019-10-22
- 录用日期: 2019-12-13
- 网络出版日期: 2020-01-09
- 刊出日期: 2020-09-14

Anti-ablative property of (C/C)/SiC-ZrC composites with different ZrC content

1.
Xi’an Aerospace Composites Research Institute, Xi’an 710025, China
2.
The Eighth military Representative Room of air force equipment department in Xi’an, Xi’an 710065, China

HTML全文

摘要

摘要: 采用先驱体转化(PIP)法制备了不同ZrC含量的(C/C)/SiC-ZrC复合材料，考察了ZrC含量对复合材料微观结构和抗烧蚀性能的影响。结果表明，氧乙炔烧蚀600 s后，(C/C)/SiC复合材料表面疏松，出现了较大的烧蚀凹坑；而(C/C)/SiC-ZrC复合材料表面相对较致密，被白色氧化物质覆盖，烧蚀率均有所降低。在较低的ZrC含量下，(C/C)/SiC-ZrC复合材料表面形成ZrO₂-SiO₂二元共熔体系氧化膜，有效抑制氧化性气氛向复合材料内部渗透，同时氧化物不断熔化和挥发，降低了复合材料烧蚀表面的温度；而当ZrC体积分数为12.4vol%时，在烧蚀过程中(C/C)/SiC-ZrC复合材料表面能形成一个ZrO₂外层/SiO₂内层的双层结构保护膜，ZrO₂是一种优异的热障材料，且导热系数较低，使烧蚀过程中烧蚀区域热扩散降低，因此(C/C)/SiC-ZrC复合材料表现为较高的表面温度，但双层氧化膜阻挡有氧气氛进一步进入复合材料内部，使复合材料表现出优异的抗烧蚀性能。
- (C/C)/SiC-ZrC /
- 复相陶瓷 /
- 抗烧蚀性能 /
- 表面温度 /
- 烧蚀机制
Abstract: The (C/C)/SiC-ZrC composites with different ZrC content were fabricated by precursor infiltration and pyrolysis (PIP) process, and effects of the ZrC content on the microstructures and anti-ablative properties of the (C/C)/SiC-ZrC composites were investigated. The results show that after oxyacetylene ablation 600 s, the surface of (C/C)/SiC composites are loose, and a rather large ablation pit appeared. Whereas the surface of (C/C)/SiC-ZrC composites are compact relatively and coverd by white oxide. Compared with (C/C)/SiC composites, ablation rate of (C/C)/SiC-ZrC composites with different ZrC content decreases. When the ZrC content is low, ZrO₂-SiO₂ binary eutectic systems a oxide film is formed on the surface, which can effectively prevent the infiltration of oxidation atmosphere, and also the oxide melt and volatilize continuously, which can reduce temperature of ablation surface of the (C/C)/SiC-ZrC composites. When the ZrC volume fraction rises to 12.4vol%, outer ZrO₂ and inner SiO₂ are formed by the oxidation of ZrC and SiC, respectively. ZrO₂ is a excellent heat barrier material with low heat conduction coefficient, which can alleviate the heat diffusion of ablation area during ablation, leading the higher surface temperature of the (C/C)/SiC-ZrC composites, but two dense oxidation film can prevent the corrosion of oxidation atmosphere, resulting in a good ablation resistance of (C/C)/SiC-ZrC composites.
- (C/C)/SiC-ZrC /
- multiphase ceramics /
- anti-ablative property /
- surface temperature /
- ablation mechnism

HTML全文

蜂窝夹层结构具有比模量大、隔音、隔热、吸能强、成本低等优点，在航空航天、交通运输、土木建筑、海洋船舶等许多领域都有广泛的应用^[1]。碳纤维增强树脂复合材料（CFRP）蒙皮-铝蜂窝夹层板主要通过蜂窝芯和蒙皮的变形、损伤及内能形式将冲击能量加以耗散，从而达到吸能的目的^[2-3]。CFRP板和铝蜂窝在生产和使用过程中都容易受到低速冲击的影响，而且蜂窝夹层板在低速冲击下的损伤一般是目视不可见的，但这种损伤却能够明显降低零部件的寿命和强度^[4-5]。因此，研究结构参数对CFRP蒙皮-铝蜂窝夹层结构低速冲击性能的影响很有必要。

国内外有大量关于蜂窝夹层板在不同载荷下力学性能的研究^[6-11]。Zhang等^[12]通过摆锤冲击试验和数值模拟方法研究了铝蜂窝夹层结构在低速冲击下的动力学特性，并得出结论：蜂窝芯是主要的吸能构件。Chen等^[13]通过落锤式低速冲击试验与数值模拟方法验证了试验与仿真结果的一致性，并探究了CFRP面板-Nomex蜂窝夹层结构的损伤机制。A. Riccio等^[14]结合落锤冲击试验和仿真模型探究了亚麻天然纤维蜂窝的吸能效果，试验与仿真结果吻合较好。Sun等^[15]通过解析法研究了金属蒙皮蜂窝夹层板在低速冲击下的动态响应，通过与已有的试验结果对比，证明了其解析法的可靠性。Vincenzo Crupi等^[16]结合试验方法和数值模型，研究了玻璃纤维增强树脂复合材料（GFRP）面板对低速冲击下铝蜂窝结构动力学性能的影响，发现GFRP面板是增强蜂窝夹层板能量吸收和承载能力的重要构件。Inés Ivañez等^[17]基于材料用户子程序建立了CFRP面板-铝蜂窝夹层梁低速冲击有限元模型，通过多次试验对比验证了仿真模型的有效性，并得出结论：蜂窝夹层结构受低速冲击时，冲击能量较大时，蒙皮吸能较多，冲击能量较小时，蜂窝芯吸能较多。

目前对不同材料蜂窝夹层板在低速冲击下动态响应的研究有很多^[18-22]，但CFRP蒙皮-铝蜂窝夹层板在低速冲击下的损伤模式相对较复杂，关于建立该类型夹层板在低速冲击下的精细化仿真模型，以探究其结构参数对抗低速冲击性能的影响以及各组件损伤对吸能占比的研究比较少，研究结构参数对该类型蜂窝夹层板低速冲击性能的影响，对其在实际工程中的应用及结构参数优化有重要的指导意义。本文使用半球头式落锤冲击试验平台和高速摄像分析平台，进行了蜂窝芯单元边长对CFRP蒙皮-铝蜂窝夹层板抗低速冲击性能影响的试验探究，根据试验结果验证了所建仿真模型的准确性。利用该数值模型进一步探究了蜂窝芯高度、蒙皮厚度和蜂窝芯壁厚等结构参数对蜂窝夹层板低速冲击性能的影响。

1. 试验

1.1 试验材料及设备

基于D7766/D7766M—11标准^[23]，搭建了如图1所示的CFRP蒙皮-铝蜂窝夹层板落锤冲击试验平台，平台主要构成设备包括：XBL-300型落锤式低速冲击试验机及数据采集装置，厂家为长春科新试验仪器有限公司。其中，控制平台、落锤、夹具和防二次冲击装置是试验机的主要构件。使用的i-speed 2系列高速摄像机及配备的ProAnalyst数码分析软件由Xcitex公司提供，如图2所示。

图 1 落锤式低速冲击试验平台

Figure 1. Low-speed impact test system of drop hammer

下载: 全尺寸图片幻灯片

图 2 高速摄像分析平台

Figure 2. High-speed camera analysis platform

下载: 全尺寸图片幻灯片

试件蒙皮均采用由威海光威复合材料有限公司提供的CFRP层合板，CFRP蒙皮是型号为T300/7901的碳纤维增强环氧树脂复合材料，铺层顺序为(45/-45/0)_2s，单层厚度为0.1 mm，单块蒙皮长L=100 mm，宽W=100 mm，厚T_f =1.2 mm；铝蜂窝采用Al 3003-H19铝箔粘接而成，蜂窝芯壁厚T_c=0.06 mm，铝芯高H_c=20 mm，正六边形蜂窝芯单元边长L_c=2 mm、4 mm、6 mm；蜂窝芯与蒙皮之间采用LJM-170胶膜进行粘接；夹层板总体尺寸为100 mm×100 mm×22.4 mm。图3为蜂窝夹层板的几何构型示意图及试件，表1为T300/7901层合板的材料参数^[24]，表2为Al 3003-H19铝箔的材料参数。

表 1 CFRP蒙皮的材料参数

Table 1. Mechanical properties of CFRP laminates

T300/7901	Value	Adhesive	Value
${E_1}{\rm{/MPa}}$	125 000	$G_{\rm{n}}^{\rm{C}}{\rm{/(N}} \cdot {\rm{m}}{{\rm{m}}^{{\rm{ - 1}}}}{\rm{)}}$	0.52
${E_2},{E_3}{\rm{/MPa}}$	11 300	$G_{\rm{s}}^{\rm{C}},G_{\rm{t}}^{\rm{C}}{\rm{/(N}} \cdot {\rm{m}}{{\rm{m}}^{ - {\rm{1}}}}{\rm{)}}$	0.92
${G_{12}},{G_{13}}{\rm{/MPa}}$	5 430	${\sigma _{{\rm{n}},\max }}{\rm{/MPa}}$	50
${G_{23}}{\rm{/MPa}}$	3 979	${\sigma _{{\rm{s}},\max }}{\rm{/MPa}}$	94
ν₁₂, ν₂₃	0.3	${\sigma _{{\rm{t}},\max }}{\rm{/MPa}}$	94
ν₂₃	0.42	${K_{\rm{n}}}{\rm{/(N}} \cdot {\rm{m}}{{\rm{m}}^{ - {\rm{3}}}}{\rm{)}}$	100 000
${X_{\rm{T}}}{\rm{/MPa}}$	2 000	${K_{\rm{s} } },{K_{\rm{t} } }{\rm{/(N} } \cdot {\rm{m} }{ {\rm{m} }^{ - {\rm{3} } } }{\rm{)} }$	100 000
${X_{\rm{C}}}{\rm{/MPa}}$	1 100
${Y_{\rm{T}}},{Z_{\rm{T}}}{\rm{/MPa}}$	80
${Y_{\rm{C}}},{Z_{\rm{T}}}{\rm{/MPa}}$	280
$S{\rm{/MPa}}$	120
Notes: E_i(i=1,2,3) is Young’s modulus in i direction; G_ij(i,j=1,2,3) is shear modulus in i-j plane; ν_ij(i,j=1,2,3) is Poisson’s ratio in i-j plane; X_t/X_c and Y_t/Y_c are the tensile/compressive strengths in 1 and 2 directions, respectively; S is shear strength; $G_{\rm{n}}^{\rm{C}}$ is toughness in tension; $G_{\rm{s}}^{\rm{C}}$ and $G_{\rm{t}}^{\rm{C}}$ are toughness components in shear; σ_n,max is maximum nominal stress of normal-only mode; σ_s,max and σ_t,max are maximum nominal stress in 1 and 2 directions, respectively; K_n is stiffness in tension; K_s and K_t are stiffness components in shear.

下载: 导出CSV

| 显示表格

表 2 Al 3003-H19铝箔材料参数

Table 2. Material properties of Al3003-H19 aluminum alloy foil

Property	Value
Density/(kg·m⁻³)	2 730
Young’s modulus/GPa	70
Poisson’s ratio	0.3
Yield strength/MPa	183

下载: 导出CSV

| 显示表格

图 3 碳纤维增强树脂复合材料（CFRP）蒙皮-铝蜂窝夹层板几何构型示意图及试件

Figure 3. Schematic of geometric configurations and specimens of aluminum honeycomb sandwich plates with carbon fiber reinforced polymer (CFRP) face sheets

下载: 全尺寸图片幻灯片

1.2 试验过程

在进行落锤冲击试验时，用夹具将试件固定，设置好落锤下降高度后进行冲击试验。落锤完成冲击时，防二次冲击装置会立刻夹紧落锤。冲击时压力传感器采集的冲击力数据由数据采集装置记录并保存；同时通过高速摄像机记录冲击过程，以便使用ProAnalyst数码分析软件获得落锤在冲击过程中的运动数据。冲击过程中的冲击试验参数：落锤质量m=2.5 kg，半球形冲头直径为25 mm，落锤下降高度h=0.48 m，理论冲击能量为E=11.76 J。由于落锤在下降过程中存在摩擦，通过高速摄像平台测得落锤接触夹层板时初始冲击速度为3.0 m/s，即实际冲击能量为11.25 J。

2. 数值计算模型

2.1 CFRP层合板本构模型

2.1.1 CFRP层合板层内渐进损伤模型

本文引用文献[24-25]中基于渐进损伤力学（CDM）方法对受损的CFRP蒙皮进行刚度退化，此方法定义了损伤因子D_f，其值随着材料损伤的扩展而逐渐增大。层合板损伤后利用基于应力描述的三维Hashin准则进行损伤判定：

纤维拉伸失效 $\left( {{\sigma _{11}} \geqslant 0} \right)$ ：

${\left( {\frac{{{\sigma _{11}}}}{{{X_{\rm{T}}}}}} \right)^2} + \frac{1}{{{S^2}}}\left( {\sigma _{12}^2 + \sigma _{13}^2} \right) \geqslant 1$

(1)

纤维压缩失效 $\left( {{\sigma _{11}} < 0} \right)$ ：

${\left( {\frac{{{\sigma _{11}}}}{{{X_{\rm{C}}}}}} \right)^2} \geqslant 1$

(2)

基体开裂失效 $\left( {{\sigma _{22}} + {\sigma _{33}} \geqslant 0} \right)$ ：

$\frac{1}{{Y_{\rm{T}}^2}}{\left( {{\sigma _{22}} + {\sigma _{33}}} \right)^2} + \frac{1}{{{S^2}}}\left( {\sigma _{23}^2 - {\sigma _{22}}{\sigma _{23}}} \right) + \frac{1}{{{S^2}}}\left( {\sigma _{12}^2 + \sigma _{13}^2} \right) \geqslant 1$

(3)

基体压溃失效 $\left( {{\sigma _{22}} + {\sigma _{33}} < 0} \right)$ ：

$\begin{split} \frac{1}{{{Y_{\rm{C}}}}}\left[ {{{\left( {\frac{{{Y_{\rm{C}}}}}{{2S}}} \right)}^2} - 1} \right]\left( {{\sigma _{22}} + {\sigma _{33}}} \right) + \frac{1}{{4{S^2}}}{\left( {{\sigma _{22}} + {\sigma _{33}}} \right)^2} + \\ \frac{1}{{{S^2}}}\left( {\sigma _{23}^2 - {\sigma _{22}}{\sigma _{23}}} \right) + \frac{1}{{{S^2}}}\left( {\sigma _{12}^2 + \sigma _{13}^2} \right) \geqslant 1 \end{split}$

(4)

其中，σ_ij（i，j=1,2,3）为CFRP面板单元的各个方向应力张量。

2.1.2 CFRP层合板层间胶层损伤模型

基于内聚力（Cohesive）单元的双线性本构模型，建立了CFRP层合板胶层模型以模拟CFRP层合板的分层损伤。内聚力单元的双线性本构模型如图4所示^[26-28]，σ轴表示单元的应力，δ横轴表示单元的相对分离位移。当δ<δ₀时，内聚力单元未破坏，其应力和相对位移呈线性关系，刚度为K，即O→A段；当δ₀≤δ<δ_max时，内聚力单元发生部分破坏，卸载时单元进入线性软化阶段，如图4中虚线部分，此时刚度为(1-D_f)K，D_f为损伤因子，即O→A→B→O段，如重加载，初始阶段内聚力单元的σ与δ的关系即对应O→B段；当δ≥δ_max时，内聚力单元完全失效。

图 4 内聚力单元的双线性本构模型

Figure 4. Bilinear constitutive model of cohesive element

下载: 全尺寸图片幻灯片

A点表示材料开始出现损伤，其对应的损伤准则为

${\left( {\frac{{{{\left( {{t_{\rm{n}}}} \right)}^2}}}{{t_{\rm{n}}^0}}} \right)^2} + {\left( {\frac{{{{\left( {{t_{\rm{s}}}} \right)}^2}}}{{t_{\rm{s}}^0}}} \right)^2} + {\left( {\frac{{{{\left( {{t_{\rm{t}}}} \right)}^2}}}{{t_{\rm{t}}^0}}} \right)^2} \geqslant 1$

(5)

其中， $t_{\rm{n}}^0$ 、 $t_{\rm{s}}^0$ 、 $t_{\rm{t}}^0$ 分别表示内聚力单元的法向和两个切向强度。AB段表示内聚力单元的损伤扩展，本文选择与试验吻合度较高的基于能量释放率的B-K准则^[24]：

$G_{\rm{n}}^{\rm{C}} + \left( {G_{\rm{s}}^{\rm{C}} - G_{\rm{n}}^{\rm{C}}} \right){\left( {\frac{{{G_{\rm{S}}}}}{{{G_{\rm{T}}}}}} \right)^\eta } = {G^{\rm{C}}}$

(6)

其中：G_S=G_s+G_t，G_T=G_n+G_S， $G_{\rm{s}}^{\rm{C}} = G_{\rm{t}}^{\rm{C}}$ ；G_i（i=n,s,t）为Cohesive单元的应变能释放率； $G_i^{\rm{C}}$ 为单元的临界应变能释放率；G^C对应图4中△OAC的面积，即单元开裂对应的总临界能释放率；η为单元材料修正系数。

2.2 CFRP蒙皮-铝蜂窝夹层板低速冲击有限元模型

基于CFRP层合板本构模型和用户材料子程序VUMAT，在ABAQUS有限元分析软件上建立蜂窝夹层板在低速冲击下的精细化仿真模型，如图5所示，蜂窝夹层板结构参数如表3所示。其中，层合板采用C3D8R单元，层间胶层采用COH3D8单元，铝蜂窝芯采用S4R壳单元；夹具和落锤均设置为刚体，夹具限制所有自由度，落锤仅保留冲击方向的自由度；假设蒙皮与蜂窝芯始终保持粘接，其接触类型设置为“Tie”，冲头与蒙皮采用“面-面接触”；为了节约计算成本，仅对临近冲击区域蒙皮进行网格细化，细化区域网格大小为1.5 mm×1.5 mm，其余区域网格大小为3 mm×3 mm，蜂窝芯网格大小为1 mm×2 mm。

表 3 CFRP蒙皮-铝蜂窝夹层板的结构参数

Table 3. Structural parameters of aluminum honeycomb sandwich plate with CFRP face sheets

Label	Cell side length ${L_{\rm{c}}}$ /mm	Core height ${H_{\rm{c}}}$ /mm	Facesheet thickness ${T_{\rm{f}}}$ /mm	Cell wall thickness ${T_{\rm{c}}}$ /mm
A1	2	20	1.2	0.06
A2	4	20	1.2	0.06
A3	6	20	1.2	0.06
B1	4	10	1.2	0.06
B2	4	20	1.2	0.06
B3	4	30	1.2	0.06
C1	4	20	1.2	0.06
C2	4	20	1.8	0.06
C3	4	20	2.4	0.06
D1	4	20	1.2	0.06
D2	4	20	1.2	0.10
D3	4	20	1.2	0.14

下载: 导出CSV

| 显示表格

图 5 CFRP蒙皮-铝蜂窝夹层板的低速冲击有限元模型

Figure 5. Finite element model of aluminum honeycomb sandwich plate with CFRP face sheets for low-velocity impact simulation

下载: 全尺寸图片幻灯片

3. 模型验证

为了证明数值模拟方法所建有限元模型的准确性，以本文的试验参数为参照，将数值模拟计算结果与试验结果进行对比。在保证其他参数相同的条件下，蜂窝芯单元边长L_c=2 mm、4 mm、6 mm的CFRP蒙皮-铝蜂窝夹层板在质量m=2.5 kg、下降高度h=0.46 m的半球头落锤的冲击下（仿真中设置落锤初始速度为3 m/s），仿真与试验的接触力历程对比如图6所示，通过高速摄像平台采集分析了落锤的速度历程，仿真与试验速度历程对比如图7所示。

图 6 试验与仿真中CFRP蒙皮-铝蜂窝夹层板的接触力历程对比

Figure 6. Comparison of experimental and numerical contact force historiesof aluminum honeycomb sandwich plate with CFRP face sheets

下载: 全尺寸图片幻灯片

图 7 试验与仿真中落锤的速度历程对比

Figure 7. Comparison of experimental and numerical velocity histories of the impactor

下载: 全尺寸图片幻灯片

通过对比图6和图7的仿真和试验数据可知：在相同能量的冲击作用下，三种蜂窝芯单元边长夹层板在仿真中的接触力峰值分别为4 037 N、3 666 N和2 854 N，试验d接触力峰值分别为3 799 N、3 476 N和2 811 N，仿真和试验接触力峰值的偏差分别为6.24%、5.47%和1.53%；仿真和试验的接触时长的差值分别为0.72 ms、0.40 ms和0.11 ms。试验与仿真在接触力历程和速度历程的总体趋势方面具有良好的一致性。

将A组所有受冲击后的蜂窝夹层板由冲击区域的中心纵向切开，会发现受冲击位置的上蒙皮与芯层之间的胶层脱粘，蜂窝芯存在折叠和屈曲，上蒙皮的下表面均有部分基体开裂的现象。试验与仿真的蜂窝芯损伤深度对比如表4所示，不同蜂窝芯单元边长的夹层板在相同的冲击能量下得到了三种不同的变形深度，A组试验中受低速冲击区域蜂窝芯变形深度分别为2.15 mm、2.77 mm和4.32 mm，仿真中的蜂窝芯变形深度δ分别为2.20 mm、3.03 mm和4.10 mm。仿真与试验得到夹层板的蜂窝芯变形量之间的偏差分别为2.33%、9.39%和5.09%。由于本文设置的低速冲击能量较小，下蒙皮未发生任何形式的损伤，故不再在文中赘述其受载情况。

表 4 试验与仿真中蜂窝芯损伤深度δ的对比

Table 4. Comparison of experimental and numerical deformations δ of honeycomb cores

Label	Top face sheet	Cross-section
Label	Top face sheet	Experiment	Simulation-Honeycomb core
A1
A2
A3

下载: 导出CSV

| 显示表格

综上所述，所搭建的低速冲击有限元仿真模型可以较为准确地预测低速冲击下CFRP蒙皮-铝蜂窝夹层板的动力学响应。

图8为三种不同蜂窝芯边长的CFRP蒙皮-铝蜂窝夹层板在质量m=2.5 kg、初始冲击速度为3 m/s的半球头落锤的冲击下，有限元模型预测的A组落锤动能历程，图9为预测的接触力-位移关系。可知，相同冲击能量下的蜂窝芯边长L_c分别为2 mm、4 mm、6 mm的CFRP蒙皮-铝蜂窝夹层板的最大变形深度依次增大，分别为5.26 mm、6.25 mm、7.40 mm。

图 8 A组CFRP 蒙皮-铝蜂窝夹层板试件落锤的动能历程对比

Figure 8. Comparison of kinetic energy histories of the impactor for group A aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

图 9 A组CFRP 蒙皮-铝蜂窝夹层板试件落锤的接触力-位移对比

Figure 9. Comparison of contact force-central displacement histories of the impactor for group A aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

通过分析可知：随着蜂窝芯单元边长的增大，夹层板受冲击时，落锤动能减小速率逐渐降低，接触时长逐渐增加，接触力峰值随之减小，夹层板受冲击位置的变形深度随之增大。可见，增大蜂窝芯单元的边长会降低蜂窝夹层板的刚度，但其吸能效果会随之增强。

4. 结构参数对CFRP蒙皮-铝蜂窝夹层板吸能的影响

4.1 蜂窝芯高度对CFRP蒙皮-铝蜂窝夹层板吸能的影响

为了确定在低速冲击载荷下，蜂窝芯高度H_c对CFRP蒙皮-铝蜂窝夹层板吸能特性的影响，在仿真中使用了H_c=10 mm、20 mm、30 mm三种不同高度的蜂窝芯，蜂窝芯单元边长L_c=4 mm，蒙皮厚度T_f=1.2 mm，蜂窝芯壁厚T_c=0.06 mm，其余条件保持不变。图10为该工况下落锤的接触力历程对比，图11为落锤的动能历程对比，图12为落锤的接触力-位移关系对比。

图 10 B组CFRP 蒙皮-铝蜂窝夹层板试件落锤的接触力历程对比

Figure 10. Comparison of contact force histories of the impactor for group B aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

图 11 B组CFRP 蒙皮-铝蜂窝夹层板试件落锤的动能历程对比

Figure 11. Comparison of kinetic energy histories of the impactor for group B aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

图 12 B组CFRP 蒙皮-铝蜂窝夹层板试件落锤的接触力-位移对比

Figure 12. Comparison of contact force-central displacement histories of the impactor for group B aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

由图10~12可知：在蜂窝芯高度H_c=10 mm、20 mm、30 mm的该组仿真中，三种蜂窝夹层板在低速冲击载荷下，落锤的接触力峰值分别为3852N、3666N、3674N，接触时长分别为5.81 ms、5.85 ms、5.82 ms，最大变形深度分别为6.25 mm、6.26 mm、6.38 mm，且落锤的接触力历程、动能历程和接触力-位移的图形几乎保持一致。B组仿真中预测的蜂窝夹层板冲击损伤如表5所示。可知，蜂窝芯受冲击区域的最终变形量δ分别为2.97 mm、3.03 mm和3.08 mm，其变形量相差较小。由该组仿真结果可见，蜂窝芯高度的变化对于CFRP蒙皮-铝蜂窝夹层板的刚度和抗低速冲击性能影响不大。

表 5 仿真预测的B组CFRP蒙皮-铝蜂窝夹层板冲击损伤

Table 5. Simulation predicted impact damage of group B aluminum honeycomb sandwich plate with CFRP face sheets

Label	Cross-section
Label	Honeycomb sandwich plate	Honeycomb core
B1
B2
B3

下载: 导出CSV

| 显示表格

4.2 蒙皮厚度对CFRP蒙皮-铝蜂窝夹层板吸能的影响

为了确定蒙皮厚度T_f对低速冲击载荷下CFRP蒙皮-铝蜂窝夹层板吸能特性的影响，在仿真中使用了T_f=1.2 mm、1.8 mm、2.4 mm三种不同厚度的CFRP蒙皮，蜂窝芯单元边长为4 mm，蜂窝芯高度H_c=20 mm，蜂窝芯壁厚T_c=0.06 mm，其余条件保持不变。图13为该工况下落锤的接触力历程对比，图14为落锤的动能历程对比，图15为落锤的接触力-位移关系对比。

图 13 C组CFRP 蒙皮-铝蜂窝夹层板试件落锤的接触力历程对比

Figure 13. Comparison of contact force histories of the impactor for group C aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

图 14 C组CFRP 蒙皮-铝蜂窝夹层板试件落锤的动能历程对比

Figure 14. Comparison of kinetic energy histories of the impactor for group C aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

图 15 C组CFRP 蒙皮-铝蜂窝夹层板试件落锤的接触力-位移对比

Figure 15. Comparison of contact force-central displacement histories of the impactor for group C aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

由图13~15可知：在蒙皮厚度T_f=1.2 mm、1.8 mm、2.4 mm的该组仿真中，三种蜂窝夹层板在相同冲击能量的低速冲击载荷下，落锤的接触力峰值分别为3666N、3923N、4314N，接触时长分别为5.85 ms、5.13 ms、4.28 ms，最大变形深度分别为6.26 mm、5.14 mm、4.39 mm。C组仿真中预测的蜂窝夹层板冲击损伤如表6所示。可知，蜂窝芯受冲击区域的最终变形量δ分别为3.03 mm、2.41 mm和2.13 mm。随着蜂窝夹层板蒙皮厚度的增大，冲击时，落锤动能的减小速率随之增大，接触时间随之减少，接触力的峰值随之增大，夹层板受冲击区域的变形量随之减小。可见，蒙皮厚度的增加会使蜂窝夹层板的刚度增大，吸能效果随之降低。

表 6 仿真预测的C组CFRP蒙皮-铝蜂窝夹层板冲击损伤

Table 6. Simulation predicted impact damage of group C aluminum honeycomb sandwich plate with CFRP face sheets

Label	Cross-section
Label	Honeycomb sandwich plate	Honeycomb core
C1
C2
C3

下载: 导出CSV

| 显示表格

4.3 蜂窝芯壁厚对CFRP蒙皮-铝蜂窝夹层板吸能的影响

为了探究低速冲击载荷下，蜂窝芯壁厚T_c对CFRP蒙皮-铝蜂窝夹层板吸能特性的影响，在仿真中使用了T_c=0.06 mm、0.10 mm、0.14 mm三种不同壁厚的蜂窝芯，蜂窝芯单元边长为4 mm，蜂窝芯高度H_c=20 mm，蒙皮厚度T_f=1.2 mm，其余条件保持不变。图16为该工况下落锤的接触力历程对比，图17为落锤的动能历程对比，图18为落锤的接触力-位移关系对比。

图 16 DCFRP 蒙皮-铝蜂窝夹层板组试件落锤的接触力历程对比

Figure 16. Comparison of contact force histories of the impactor for group D aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

图 17 DCFRP 蒙皮-铝蜂窝夹层板组试件落锤的动能历程对比

Figure 17. Comparison of kinetic energy histories of the impactor for group D specimens

下载: 全尺寸图片幻灯片

图 18 D组CFRP 蒙皮-铝蜂窝夹层板试件落锤的接触力-位移对比

Figure 18. Comparison of contact force-central displacement histories of the impactor for group D aluminum honeycomb sandwich plate with CFRP face sheets specimens

下载: 全尺寸图片幻灯片

根据图16~18可知：在蜂窝芯壁厚T_c=0.06 mm、0.10 mm、0.14 mm的该组仿真中，三种蜂窝夹层板在相同冲击能量的低速冲击载荷下，落锤的接触力峰值分别为3666N、4153N、4454N，接触时长分别为5.85 ms、5.31 ms、4.99 ms，最大变形深度分别为6.26 mm、5.50 mm、5.28 mm。仿真预测的D组蜂窝夹层板冲击损伤如表7所示。可知，蜂窝芯受冲击区域的最终变形量δ分别为3.03 mm、2.83 mm和2.77 mm。随着蜂窝夹层板蜂窝芯壁厚的增大，在低速冲击下，落锤动能的减小速率随之增大，接触时间随之减少，落锤与夹层板之间的接触力峰值随之增大，夹层板受冲击区域的变形量随之减小。可见，增加蜂窝芯的壁厚会使蜂窝夹层板的刚度增大，抗低速冲击的性能增强。

表 7 仿真预测的D组CFRP蒙皮-铝蜂窝夹层板冲击损伤

Table 7. Simulation predicted impact damage of group D aluminum honeycomb sandwich plate with CFRP face sheets

Label	Cross-section
Label	Honeycomb sandwich plate	Honeycomb core
D1
D2
D3

下载: 导出CSV

| 显示表格

4.4 各组件吸能占比

表8为不同结构参数的CFRP蒙皮-铝蜂窝夹层板在相同能量的低速冲击下，各组件的吸能占比情况。可知：在本文设置的较低能量级冲击载荷下，CFRP蒙皮-铝蜂窝夹层板的吸能占比均大于60%，足以证明蜂窝夹层结构具有较好的冲击吸能特性；在该规格的冲击条件下，CFRP蒙皮-铝蜂窝夹层板最主要的吸能组件是铝蜂窝芯，上蒙皮吸能次之，下蒙皮吸能最少。由A组夹层板吸能占比分析可知，随着蜂窝芯单元边长的增大，蜂窝芯吸能占比逐渐增大，上蒙皮吸能占比逐渐减小，夹层板总吸能逐渐增大；由B组数据分析可知，随着蜂窝芯高度的增大，蜂窝芯吸能占比逐渐增大，上蒙皮吸能逐渐减少，夹层板总吸能逐渐增大；由C组数据分析可知，随着蒙皮厚度的增大，蜂窝芯层吸能占比逐渐减小，上蒙皮吸能占比逐渐增大，夹层板总吸能逐渐增大；由D组数据分析可知，随着蜂窝芯壁厚的增大，蜂窝芯吸能占比逐渐减小，上蒙皮吸能占比逐渐增大，夹层板总吸能逐渐增大。综上所述，在本文所设置的载荷条件下，增大蜂窝芯单元边长、蜂窝芯高度、蒙皮厚度和蜂窝芯壁厚，均能增加蜂窝夹层板的吸能量。然而吸能效果需结合碰撞时的接触时长、接触力峰值、变形损伤深度和吸能量综合考虑，故需结合图表才能得到合理的结论。

表 8 不同结构参数的CFRP蒙皮-铝蜂窝夹层板组件吸能占比

Table 8. Absorbed energy rates by components of aluminum honeycomb sandwich plate with CFRP face sheets with different structural parameters %

Parameters	Cell side length ${L_{\rm{c}}}$			Core height ${H_{\rm{c}}}$			Facesheet thickness ${T_{\rm{f}}}$			Cell wall thickness ${T_{\rm{c}}}$
Parameters	A1	A2	A3	B1	B2	B3	C1	C2	C3	D1	D2	D3
Upper plate	20.43	17.29	16.49	19.53	17.29	16.04	17.29	26.97	31.72	17.29	20.88	30.02
Honeycomb	45.16	50.01	58.60	38.45	50.01	57.35	50.01	43.20	39.78	50.01	45.61	41.75
Lower plate	1.88	1.70	2.51	2.15	1.70	1.70	1.70	1.08	0.72	1.70	2.69	3.32
Impactor	32.53	31.00	22.40	39.87	31.00	24.91	31.00	28.05	27.78	31.00	30.82	24.91
Sandwich plate	67.37	69.00	77.60	60.13	69.00	75.09	69.00	71.95	72.22	69.00	69.18	75.09

下载: 导出CSV

| 显示表格

5. 结论

建立了碳纤维增强树脂复合材料（CFRP）蒙皮-铝蜂窝夹层板精细化低速冲击仿真模型，通过对比不同蜂窝芯边长的夹层板在相同冲击能量下的仿真与试验结果，证明了试验与仿真结果具有较好的一致性。进而利用该数值模型探究了蜂窝芯高度、蒙皮厚度和蜂窝芯壁厚等结构参数对于蜂窝夹层板低速冲击性能的影响。

（1）随着蜂窝芯单元边长的增大，夹层板受冲击时，落锤动能减小速率逐渐降低，接触时长逐渐增加，接触力峰值随之减小，夹层板受冲击位置的变形深度随之增大。可见，蜂窝芯单元边长尺寸的增大会降低蜂窝夹层板的刚度，但其吸能效果会随之增强。

（2）蜂窝夹层板的芯层高度对夹层板的刚度及抗低速冲击性能影响较小。

（3）增大蜂窝夹层板的蒙皮厚度，可以提高夹层板的刚度，但会降低夹层板的吸能效果。

（4）增大蜂窝芯壁厚，会提高蜂窝夹层板的刚度和抗低速冲击性能。

（5）在本文的低速冲击条件下，CFRP蒙皮-铝蜂窝夹层板的蜂窝芯是主要的吸能组件，上蒙皮吸能占比次之，下蒙皮吸能最少。

图 1 (C/C)/SiC-ZrC复合材料的XRD图谱

Figure 1. XRD patterns of (C/C)/SiC-ZrC composites

下载: 全尺寸图片幻灯片

图 2 (C/C)/SiC-ZrC复合材料的SEM图像((a)~(c))及EDS图谱((d)~(f))

Figure 2. SEM images ((a)–(c)) and EDS spectra((d)–(f)) of (C/C)/SiC-ZrC composites

下载: 全尺寸图片幻灯片

图 3 (C/C)/SiC和(C/C)/SiC-ZrC复合材料烧蚀后的照片

Figure 3. Photographs of (C/C)/SiC and (C/C)/SiC-ZrC composites after ablation((a) (C/C)/SiC; (b) (C/C)/SiC-ZrC-1^#; (c) (C/C)/SiC-ZrC-2^#; (d) (C/C)/SiC-ZrC-3^#)

下载: 全尺寸图片幻灯片

图 4 (C/C)/SiC-ZrC-1^#复合材料烧蚀中心的SEM图像

Figure 4. SEM images of ablation left of (C/C)/SiC-ZrC-1^# composites((a) glassy mixture; (b) surface fibers)

下载: 全尺寸图片幻灯片

图 5 (C/C)/SiC-ZrC-2^#复合材料烧蚀中心的SEM图像

Figure 5. SEM images of ablation left of (C/C)/SiC-ZrC-2^# composites((a) glassy mixture; (b) surface fibers)

下载: 全尺寸图片幻灯片

图 6 (C/C)/SiC-ZrC-3^#复合材料烧蚀中心的SEM图像

Figure 6. SEM images of ablation left of (C/C)/SiC-ZrC-3^# composites((a) glassy mixture; (b) surface fibers)

下载: 全尺寸图片幻灯片

图 7 图4~图6中(C/C)/SiC-ZrC复合材料的两种氧化物的EDS图谱

Figure 7. EDS spectra of two oxidates of (C/C)/SiC-ZrC composites in Fig.4–Fig.6((a) A^#; (b) B^#)

下载: 全尺寸图片幻灯片

图 8 (C/C)/SiC和(C/C)/SiC-ZrC复合材料烧蚀中心的表面温度

Figure 8. Surface temperature of ablation left of (C/C)/SiC and (C/C)/SiC-ZrC composites

下载: 全尺寸图片幻灯片

表 1 (C/C)/SiC-ZrC复合材料各组分的体积分数

Table 1 Volume fraction of different components of (C/C)/SiC-ZrC composites

Sample	Density/(g·cm^–³)	Volume fraction/vol%				Porosity/%
Sample	Density/(g·cm^–³)	V_f	V_C	V_SiC	V_ZrC	Porosity/%
(C/C)/SiC	1.820	24.5	38.9	28.8	0	7.83
(C/C)/SiC-ZrC-1^#	1.861	24.3	38.9	27.8	4.02	4.98
(C/C)/SiC-ZrC-2^#	1.870	24.0	38.9	25.1	7.22	4.78
(C/C)/SiC-ZrC-3^#	1.950	23.5	38.9	21.3	12.40	3.90
Notes: V_f—Volume fraction of carbon fiber; V_C—Volume fraction of pyrolytic carbon; V_SiC—Volume fraction of SiC; V_ZrC—Volume fraction of ZrC fraction of SiC; V_ZrC—Volume fraction of ZrC.

下载: 导出CSV

表 2 (C/C)/SiC和(C/C)/SiC-ZrC复合材料的线烧蚀率、质量烧蚀率及烧蚀深度

Table 2 Linear ablation rate, mass ablation rate and ablation depth of (C/C)/SiC and (C/C)/SiC-ZrC composites

Sample	Ablation time/s	Linear ablation rate/(mm·s^–1)	Mass ablation rate/(g·s^–1)	Ablation depth/mm
(C/C)/SiC	600	0.0047	0.0033	2.79
(C/C)/SiC-ZrC-1^#	600	0.0027	0.0019	1.61
(C/C)/SiC-ZrC-2^#	600	0.0019	0.0016	1.14
(C/C)/SiC-ZrC-3^#	600	0.0012	0.0004	0.72

下载: 导出CSV

参考文献(24)

[1]	刘小冲, 成来飞, 张立同, 等. C/SiC复合材料在空间环境中的性能研究进展[J]. 材料导报, 2013, 27(9):127-130. DOI: 10.3969/j.issn.1005-023X.2013.09.028 LIU Xiaochong, CHENG Laifei, ZHANG Litong, et al. Research progress on properties of the C/SiC composites under the space environment[J]. Materials Review,2013,27(9):127-130(in Chinese). DOI: 10.3969/j.issn.1005-023X.2013.09.028
[2]	张立同, 成来飞. 连续纤维增韧陶瓷基复合材料可持续发展战略探讨[J]. 复合材料学报, 2007, 24(2):1-6. DOI: 10.3321/j.issn:1000-3851.2007.02.001 ZHANG Litong, CHENG Laifei. Discussion on strategies of sustainable development of continous fiber reinforced ceramic matrix composites[J]. Acta Materiae Compositae Sinica,2007,24(2):1-6(in Chinese). DOI: 10.3321/j.issn:1000-3851.2007.02.001
[3]	TRIANTOU K I, MERGIA K, PEREZ B, et al. Thermal shock performance of carbon-bonded carbon fiber composite and ceramic matrix composite joints for thermal protection re-entry applications[J]. Composites Part B: Engineering,2017,111:270-278. DOI: 10.1016/j.compositesb.2016.12.020
[4]	JAYASEELAN D D, XIN Y, VANDEPERRE L, et al. Development of multilayered thermal protection system (TPS) for aerospace applications[J]. Composites Part B: Engineering,2015,79:392-405. DOI: 10.1016/j.compositesb.2015.04.017
[5]	PAUL A, VENUGOPAL S, BINNER J G P, et al. UHTC carbon fibre composites: Preparation, oxyacetylene torch testing and characterization[J]. Journal of the European Ceramic Society,2013,33(2):423-432. DOI: 10.1016/j.jeurceramsoc.2012.08.018
[6]	ZHAO L Y, JIA D C, DUAN X M, et al. Oxidation of ZrC-30vol% SiC composite in air from low to ultrahigh temperature[J]. Journal of the European Ceramic Society,2012,32(4):947-954. DOI: 10.1016/j.jeurceramsoc.2011.10.024
[7]	LI K Z, XIE J, FU Q G, et al. Effects of porous C/C density on the densification behavior and ablation property of C/C-ZrC-SiC composites[J]. Carbon,2013,57:161-168. DOI: 10.1016/j.carbon.2013.01.059
[8]	武海棠, 魏玺, 于守泉, 等. 整体抗氧化C/C-ZrC-SiC复合材料的超高温烧蚀性能研究[J]. 无机材料学报, 2011, 26(8):852-856. WU Haitang, WEI Xi, YU Shouquan, et al. Ablation performances of multi-phased C/C-ZrC-SiC ultra-high temperature composites[J]. Journal of Inorganic Materials,2011,26(8):852-856(in Chinese).
[9]	WANG Z, DONG S M, DING Y S, et al. Mechanial poperties and microstructures of C<sub>f</sub>/SiC-ZrC composites using T700SC carbon fibers as reinforcements[J]. Ceramics International,2011,37(3):695-700. DOI: 10.1016/j.ceramint.2010.09.048
[10]	LIU C X, CAO L X, CHEN J X, et al. Microstructure and ablation behavior of SiC coated C/C-SiC-ZrC composites prepared by a hybrid infiltration process[J]. Carbon,2013,65:196-205. DOI: 10.1016/j.carbon.2013.08.014
[11]	LI Y, CHEN S A, MA X, et al. Effects of high-tempenrature annealing on the microstructure and properties of C/SiC-ZrC composites[J]. Ceramic International,2016,42(4):5171-5176. DOI: 10.1016/j.ceramint.2015.12.039
[12]	XIE J, LI K Z, LI H J, et al. Ablation behavior and mechanism of C/C-ZrC-SiC composites under an oxyacetylene torch at 3 000℃[J]. Ceramics International,2013,39(4):4171-4178. DOI: 10.1016/j.ceramint.2012.10.273
[13]	王玲玲, 嵇阿琳, 崔红, 等. ZrC改性C/C-SiC复合材料的力学和抗烧蚀性能[J]. 复合材料学报, 2016, 33(2):373-378. WANG Lingling, JI Alin, CUI Hong, et al. Mechanical and anti-ablative preformances of C/C-SiC composites modified by ZrC[J]. Acta Materiae Compositae Sinica,2016,33(2):373-378(in Chinese).
[14]	杨星, 崔红, 闫联生, 等. 陶瓷前驱体配比对C/C-ZrC-SiC复合材料烧蚀性能的影响[J]. 材料导报, 2015, 29(2):52-56. YANG Xing, CUI Hong, YAN Liansheng, et al. Effects of ceramic precursor ratio on ablation properties of C/C-ZrC-SiC composites[J]. Materials Review,2015,29(2):52-56(in Chinese).
[15]	杨星, 崔红, 闫联生, 等. ZrC改性C/C-SiC复合材料的烧蚀性能研究[J]. 化工新型材料, 2014, 42(11):169-171. YANG Xing, CUI Hong, YAN Liansheng, et al. Investigation on ablation performance of ZrC modified C/C-SiC composites[J]. New Chemical Materials,2014,42(11):169-171(in Chinese).
[16]	国防科学技术工业委员会. 烧蚀材料烧蚀试验方法: GJB 323A—96[S]. 北京: 国防科工委军标出版社, 1997. National Defense Science and Technology Commission. Test methods for ablation of ablators: CJB 323A—96[S]. Beijing: National Defense Science and Technology Commission Military Publishing House, 1997(in Chinese).
[17]	庄磊, 付前刚, 李贺军, 等. 聚合物浸渍裂解法制备C/C-ZrC-SiC复合材料的氧化行为及抗烧蚀性能研究[J]. 中国材料进展, 2015, 34(6):425-431. ZHUANG Lei, FU Qiangang, LI Hejun, et al. Oxidation behavior and ablation resistance of C/C-ZrC-SiC composites prepared by polymer infiltration and pyrolysis[J]. Materials China,2015,34(6):425-431(in Chinese).
[18]	王其坤, 胡海峰, 陈朝辉, 等. 先驱体转化法制备2D C/SiC-ZrB<sub>2</sub>复合材料及其性能[J]. 复合材料学报, 2009, 26(1):108-112. DOI: 10.13801/j.cnki.fhclxb.2009.01.030 WANG Qikun, HU Haifeng, CHEN Zhaohui, et al. Perparation and property of 2D C/SiC-ZrB<sub>2</sub> composites via precursor infiltration and pyrolysis[J]. Acta Materiae Compositae Sinica,2009,26(1):108-112(in Chinese). DOI: 10.13801/j.cnki.fhclxb.2009.01.030
[19]	童长青, 成来飞, 刘永胜, 等. 2D C/SiC-ZrB<sub>2</sub>复合材料的烧蚀性能[J]. 航空材料学报, 2012, 32(2):69-72. DOI: 10.3969/j.issn.1005-5053.2012.2.013 TONG Changqing, CHENG Laifei, LIU Yongsheng, et al. Ablation properties of 2D C/SiC-ZrB<sub>2</sub> composite[J]. Journal of Aeronautical Materials,2012,32(2):69-72(in Chinese). DOI: 10.3969/j.issn.1005-5053.2012.2.013
[20]	FENG B, LI H J, ZHANG Y L, et al. Effect of SiC/ZrC ratio on the mechanical and ablation properties of C/C-SiC-ZrC composites[J]. Corrosion Science,2014,82:27-35.
[21]	王少雷, 李红, 任慕苏, 等. ZrC-SiC-C/C复合材料的制备及其烧蚀性能[J]. 复合材料学报, 2017, 34(5):1040-1047. WANG Shaolei, LI Hong, REN Musu, et al. Fabrication and ablation performances of ZrC-SiC-C/C composites[J]. Acta Materiae Compositae sinica,2017,34(5):1040-1047(in Chinese).
[22]	LI Z Q, LI H J, ZHANG S Y, et al. Microstructures and ablation properties of C/C-SiC-ZrC composites prepared using C/C skeletons with various densities[J]. Ceramic International,2013,39(7):8173-8181. DOI: 10.1016/j.ceramint.2013.03.093
[23]	解静, 李克智, 付前刚, 等. 聚合物浸渍裂解法制备C/C-ZrC-SiC-ZrB<sub>2</sub>复合材料及其性能研究[J]. 无机材料学报, 2013, 28(6):605-610. DOI: 10.3724/SP.J.1077.2013.12450 XIE Jing, LI Kezhi, FU Qiangang, et al. Preparation and properties of C/C-ZrC-SiC-ZrB<sub>2</sub> composites via polymer infiltration and pyrolysis[J]. Journal of Inorganic Materials,2013,28(6):605-610(in Chinese). DOI: 10.3724/SP.J.1077.2013.12450
[24]	张强, 赵景鹏, 崔红, 等. C/C-ZrB<sub>2</sub>-SiC 复合材料的制备及其性能[J]. 硅酸盐学报, 2016, 44(7):1003-1007. ZHANG Qiang, ZHAO Jingpeng, CUI Hong, et al. Preparation and properties of C/C-ZrB<sub>2</sub>-SiC composites[J]. Journal of the Chinese Ceramic Society,2016,44(7):1003-1007(in Chinese).