MoS<sub>2</sub>/生物质碳复合材料的制备与吸波性能

谢文瀚; 耿浩然; 柳扬; 赵婷婷; 张璇; 董丽杰

doi:10.13801/j.cnki.fhclxb.20210715.001

MoS₂/生物质碳复合材料的制备与吸波性能

武汉理工大学　材料复合新技术国家重点实验室，武汉 430070

基金项目: 武汉理工大学三亚科教创新园开放基金(2020KF0025)

详细信息

通讯作者:
董丽杰，博士，教授，博士生导师，研究方向为先进复合材料　E-mail：dong@whut.edu.cn

中图分类号: TB332
计量
- 文章访问数: 1925
- HTML全文浏览量: 532
- PDF下载量: 102
出版历程
- 收稿日期: 2021-05-11
- 修回日期: 2021-06-13
- 录用日期: 2021-06-29
- 网络出版日期: 2021-07-14
- 刊出日期: 2022-03-22

Preparation and microwave absorbing properties of MoS₂/biomass carbon composite

State Key Laboratory of Advanced Technology for Materials Synthesis and Processing, Wuhan University of Technology, Wuhan 430070, China

摘要

摘要: 为了解决MoS₂吸波材料电导率低的问题，以柚子皮(SP)为原料，采用一锅水热及高温煅烧法制备了MoS₂/生物质碳(BC)复合材料。通过调节初始Mo源、S源的含量来调控该复合材料中MoS₂的含量。微观形貌、结构和电磁参数结果表明，随着复合材料中MoS₂含量的增加，MoS₂在BC表面由片状零散分布转变为花状包覆，MoS₂/BC复合材料的电导率和复介电常数逐渐降低。通过调节MoS₂与BC的比例，实现对MoS₂/BC复合材料电磁参数的有效调控，进而优化其阻抗匹配特性。MoS₂的花状结构有利于促进电磁波的多重反射/散射。同时，MoS₂与BC之间存在丰富界面，有利于促进界面极化，增强MoS₂/BC复合材料对电磁波的衰减能力。所制备的MoS₂/BC-0.8最小反射率损失(RL)值为−40.1 dB，有效吸收带宽为5.9 GHz(11.1~17.0 GHz)。
- 生物质 /
- MoS₂ /
- 柚子皮 /
- 复合材料 /
- 阻抗匹配 /
- 吸波性能
Abstract: In order to solve the problem of low conductivity of MoS₂ absorbing material, the MoS₂/biomass carbon (BC) composite material was prepared by using shaddock peel (SP) as the raw material, using one-pot hydro-thermal and high-temperature calcination methods. The content of MoS₂ in the composite material was adjusted by adjusting the content of the initial Mo source and S source. The results of microscopic morphology, structure and electromagnetic parameters show that with the increase of the MoS₂ content in the composite material, the scattered distribution of MoS₂ on the BC surface changes from flakes to flower-like coatings, and the conductivity and complex permittivity of MoS₂/BC composites gradually decrease. By adjusting the ratio of MoS₂ to BC, the effective control of the electromagnetic parameters of the MoS₂/BC composite material is realized, and its impedance matching characteristics are optimized. The flower-like structure of MoS₂ facilitates the multiple reflection/scattering of electromagnetic waves. At the same time, there are abundant interfaces between flower-like MoS₂ and BC, which is beneficial to promote interface polarization and enhance the attenuation ability of MoS₂/BC composites to electromagnetic waves. The prepared MoS₂/BC-0.8 has a minimum reflectance loss (RL) value of –40.1 dB, and an effective absorption bandwidth of up to 5.9 GHz (11.1-17.0 GHz).
- biomass /
- MoS₂ /
- shaddock peel /
- composites /
- impedance match /
- microwave absorption properties

HTML全文

复合材料具有比强度和比刚度高、可设计性强、疲劳性能好、重量轻、耐腐蚀等优异性能，被广泛应用于航空航天领域^[1]。复合材料薄壁结构承载效率较高，是飞机中常见的结构形式，其中机身侧壁板、梁腹板等结构主要承受剪切载荷的作用。在剪切载荷作用下，机身侧壁板的蒙皮容易局部失稳，已有研究表明，蒙皮失稳后不会立即失效，还能继续承载，表现出较高的承载能力^[2-4]。因此，在复合材料薄壁结构设计中应考虑其剪切稳定性，充分发挥后屈曲承载潜力，这对减轻飞机结构重量、提高结构承载效率和确保飞机服役过程中的安全具有重要工程价值。

目前，国内外学者对复合材料薄壁结构的稳定性研究多集中在压缩载荷工况下^[5-8]，对剪切稳定性的研究相对较少。谭翔飞等^[9]基于蒂尔曼和迪奥模型，进行了加强筋间蒙皮剪切屈曲及后屈曲性能理论分析，得出的理论屈曲载荷与试验屈曲载荷的相对误差为7.2%，后屈曲阶段的理论载荷-应变曲线与试验结果的相对误差在10%以内，并确定了后屈曲角随屈曲比的变化规律。汪厚冰等^[10]对复合材料帽形加筋壁板的剪切屈曲性能进行了试验、理论分析和数值模拟，并指出选择合适的几何初始缺陷系数利用几何非线性分析方法可模拟蒙皮在剪切载荷作用下的屈曲过程。李真等^[11]通过剪切稳定性试验及张力场分析研究了复合材料机身单曲率壁板的剪切失稳及破坏特性。以上研究工作均未考虑冲击损伤，而复合材料结构在使用和维护过程中受到工具坠落、跑道碎石和冰雹等低速冲击的作用容易产生损伤，导致剩余强度急剧下降^[12]，因此冲击损伤的影响不可忽视。

对含冲击损伤复合材料结构的强度预测，一般采用损伤累积法^[13]和软化夹杂法^[14-16]。损伤累积法需要进行冲击及压缩等外载作用的全过程分析，计算规模庞大，收敛困难，不便于工程使用。基于实际冲击损伤，采用软化夹杂法对冲击损伤进行等效是一种可行且高效的方法。Fardin等^[14]将冲击损伤简化为具有不同材料参数的椭圆夹杂物，Debski等^[15]通过减少厚度、降低材料属性将冲击损伤等效，两位学者均高效地计算了复合材料结构冲击后的剩余强度。Li等^[16]基于软化夹杂模型建立了一种边缘冲击损伤的简化方法，快速而准确地预测了加筋板的后屈曲行为。然而，以上软化夹杂模型更适用于冲击后压缩剩余强度预测，对于本文冲击后剪切强度预测，还应考虑剪切非线性效应的影响。

本文以无损伤及含冲击损伤的复合材料层合板为研究对象，基于数字图像相关方法(Digital image correlation，DIC)^[17]对其在剪切载荷作用下的屈曲后屈曲行为进行了实时测量，分析和总结了冲击损伤对屈曲波形和承载能力的影响。随后，采用软化夹杂法将冲击损伤等效简化，将损伤的几何边界信息写入损伤模型中，建立了考虑剪切非线性效应的复合材料渐进损伤模型，并对试验过程进行数值仿真。最后，讨论了数值仿真结果与试验结果的一致性。

1. 试验

1.1 试件

试件材料为碳纤维增强环氧树脂基复合材料CCF300/BA3202，单层名义厚度为0.155 mm，材料属性见表1，参考ASTM标准^[18-22]测得。试件外廓尺寸为370 mm×250 mm(详见图1)，共4件，编号为T1~T4，所有试件铺层顺序为[45/45/0/45/90/0/−45/0]_S，名义厚度为2.48 mm。试件四周用玻璃钢加强，加强片与蒙皮之间用J-116 B胶粘接，加强片厚度为2.5 mm。

表 1 CCF300/BA3202碳纤维增强环氧树脂复合材料性能参数

Table 1. Material properties of CCF300/BA3202 carbon fiber reinforced epoxy composite

Parameter	Value
E₁/GPa	118
E₂/GPa	8.98
G₁₂/GPa	4.21
ν₁₂	0.306
X_T/MPa	1835
X_C/MPa	1296
Y_T/MPa	82.5
Y_C/MPa	240
S_L/MPa	166
G_IC/(N·mm⁻¹)	0.744
G_IIC/(N·mm⁻¹)	1.90
Notes： ${E_1}$ —Longitudinal elastic modulus; ${E_2}$ —Ransverse elastic modulus; ${G_{12}}$ —Shear modulus; ${\nu _{{\text{12}}}}$ —Poisson's ratio; ${X_{\text{T}}}$ —Longitudinal tensile strength; ${X_{\text{C}}}$ —Longitudinal compressive strength; ${Y_{\text{T}}}$ —Transverse tensile strength; ${Y_{\text{C}}}$ —Transverse compressive strength; ${S_{\text{L}}}$ —Shear strength; G_IC—Mode I fracture toughness ; G_IIC—Mode II fracture toughness.

下载: 导出CSV

| 显示表格

1.2 冲击损伤引入

试验前，在T4试件正中心引入了冲击损伤。冲击采用落锤式冲击设备，冲头质量为1.023 kg，端部为半球形，直径为16 mm，冲击能量为8 J。冲击后立即用凹坑深度仪测量，得到的凹坑深度为0.27 mm。随后，采用EPOCH XT超声波探伤仪进行无损检测，并用白色记号笔对损伤进行标识。A扫结果显示：分层损伤尺寸为18 mm×12 mm，如图2所示。

图 1 CCF300/BA3202复合材料层合板尺寸

Figure 1. Dimensions of CCF300/BA3202 composite laminates

下载: 全尺寸图片幻灯片

图 2 T4试件冲击支持方式及分层尺寸

Figure 2. Impact support mode and delamination size of T4 specimen

下载: 全尺寸图片幻灯片

1.3 试验方法

正式试验前，在试件一侧表面喷涂白漆，待白漆干后喷涂黑漆，形成随机的散斑图案。在试件的另一侧表面粘贴花片，花片的0°、45°、90°方向的应变片与试件铺层方向一致。T1~T3试件应变片粘贴位置及编号见图3，通过1^#~4^#花片检验加载对称性，通过中心处的5^#花片判断初始屈曲。T4试件的5^#花片沿着90°方向移动了10 mm。

图 3 无损伤试件(T1~T3)应变片分布及编号

Figure 3. Strain gauges distributions and numbers of undamaged specimens (T1-T3)

下载: 全尺寸图片幻灯片

试验夹具和试件安装如图4所示，采用对角拉伸的方法进行剪切载荷施加。试验在MTS1000 kN试验机上进行，采用道姆公司制造的ARAMIS 3D数字散斑系统测量面外变形，同时采用中国飞机强度研究所自主研发的ST24数据采集系统测量应变。

图 4 CCF300/BA3202复合材料层合板剪切试验方法

Figure 4. Shear test method for CCF300/BA3202 composite laminates

下载: 全尺寸图片幻灯片

2. 数值仿真模型

2.1 失效准则

2.1.1 纤维失效准则

复合材料主要失效模式包括纤维拉伸/压缩失效和基体的拉伸/压缩失效等。Hashin准则^[23]是最常用的失效准则，但Li等^[24]从理论上分析了Hashin纤维拉伸失效判据中包含的剪切应力项可以忽略，并给出了预测纤维拉伸/压缩失效的统一表达式：

${f_{{\text{FF}}}} = \frac{{\left( {{X_{\rm{C}}} - {X_{\rm{T}}}} \right){\sigma _1}}}{{{X_{\text{T}}}{X_{\text{C}}}}} + \frac{{\sigma _1^2}}{{{X_{\text{T}}}{X_{\text{C}}}}}$

(1)

式中： $f_{\mathrm{FF}}$ 为纤维失效函数， $f_{\mathrm{FF}}>1$ 时表征纤维失效； $\sigma_{1}$ 为纤维方向应力； $X_{\mathrm{T}} 、 X_{\mathrm{C}}$ 分别为单向板纤维方向拉伸、压缩强度。

2.1.2 纤维间失效准则

对于纤维间失效，Hashin准则难以反映横向压缩对基体剪切失效的抑制作用。Li等^[24]从Mohr失效面理论的基本假设出发，建立了一种纤维间失效准则，可有效预测横向压缩的增强效应。而且，不同于Puck等^[25]提出的纤维间失效准则，该准则所有系数均可由基本强度及断裂面角度计算给出，准则形式如下：

${f_{{\text{IFF}}}} = \left\{ {\begin{array}{*{20}{c}} {{A_1}{\sigma _{\text{n}}} + {A_2}\sigma _{\text{n}}^2 + {B_2}\tau _{{\text{nl}}}^2 + {C_2}\tau _{{\text{nt}}}^2}&{{\sigma _{\text{n}}} \geqslant 0} \\ {{A_1}{\sigma _{\text{n}}} + {B_2}\tau _{{\text{nl}}}^2 + {C_2}\tau _{{\text{nt}}}^2}&{{\sigma _{\text{n}}} < 0} \end{array}} \right.$

(2)

式中： $f_{\text {IFF }}$ 为纤维间失效函数， $f_{\mathrm{IFF}}>1$ 时表征纤维间失效，即基体失效； ${\sigma }_{\text{n}}、{\tau }_{\text{nt}}、{\tau }_{\text{nl}}$ 为潜在断裂面上的应力(如图5所示)； $A_{1}、A_{2}、B_{2}、 C_{2}$ 为系数，计算公式如下：

$\left\{ \begin{gathered} {\sigma _{\text{n}}} = {\sigma _2}{\cos ^2}{\theta _{{\text{fp}}}} + {\sigma _3}{\sin ^2}{\theta _{{\text{fp}}}} + {\tau _{23}}\sin 2{\theta _{{\text{fp}}}} \\ {\tau _{{\text{nt}}}} = \dfrac{{\left( {{\sigma _3} - {\sigma _2}} \right)\sin 2{\theta _{{\text{fp}}}}}}{2} + {\tau _{23}}\cos 2{\theta _{{\text{fp}}}} \\ {\tau _{{\text{nl}}}} = {\tau _{13}}\sin {\theta _{{\text{fp}}}} + {\tau _{12}}\cos {\theta _{{\text{fp}}}} \\ \end{gathered} \right.$

(3)

$\left\{ \begin{gathered} {A_1} = \dfrac{{{{\sin }^2}{\theta _0} - {{\cos }^2}{\theta _0}}}{{{Y_{\text{C}}}{{\cos }^4}{\theta _0}}} \\ {A_2} = \dfrac{1}{{Y_{\text{T}}^2}}\left[ {1 - {Y_{\text{T}}}\left( {\dfrac{{{{\sin }^2}{\theta _0} - {{\cos }^2}{\theta _0}}}{{{Y_{\text{C}}}{{\cos }^4}{\theta _0}}}} \right)} \right] \\ {B_2} = \dfrac{1}{{S_{\text{L}}^2}} \\ {C_2} = \dfrac{1}{{Y_{\text{C}}^2{{\cos }^4}{\theta _0}}} \\ \end{gathered} \right.$

(4)

式中： ${\sigma }_{2}、{\sigma }_{3}、{\tau }_{12}、{\tau }_{13}、{\tau }_{23}$ 为材料主轴坐标系下的应力分量； ${Y}_{\text{T}}、{Y}_{\text{C}}、{S}_{\text{L}}$ 分别为横向拉伸强度、横向压缩强度和面内剪切强度； $\theta_{0}$ 为纯横向压缩失效时的断裂角，本文参考文献[25]取为53°； $\theta_{\mathrm{fp}}$ 为潜在断裂面角度。

图 5 基体潜在断裂面上的应力分量^[24]

Figure 5. Stress components on potential fracture surface of matrix^[24]

${\sigma _1}$ ,

${\sigma _2}$ ,

${\sigma _3}$ ,

${\tau _{12}}$ ,

${\tau _{13}}$ ,

${\tau _{23}}$ ,

${\tau _{21}}$ ,

${\tau _{31}}$ ,

${\tau _{32}}$ —Stress components in the material's principal axis coordinate system;

${\theta _{{\text{fp}}}}$ —Angle of the potential fracture surface of matrix;

${\sigma _{\text{n}}}$ ,

${\tau _{{\text{nt}}}}$ ,

${\tau _{{\text{nl}}}}$ —Stress components on potential fracture surface of matrix

下载: 全尺寸图片幻灯片

纤维间失效函数 $f_{\mathrm{IFF}}$ 的最大值及相应的潜在断裂面角度 $\theta_{\mathrm{fp}}$ 可通过遍历法^[25]、分区黄金搜索法^[26]、分区抛物线法^[27]求解，综合考虑效率和精度本文采用了分区黄金搜索法。

2.2 剪切非线性模型

在面内纵横剪切试验中发现，CCF300/BA3202表现出了明显的剪切非线性行为。因此，本文采用多项式函数描述其剪切非线性响应。同时，注意到剪应力是剪应变的奇函数，剪切应力-应变关系如下：

${\tau _{12}} = \left\{ \begin{aligned} &{{G_1}{\gamma _{12}} + {G_2}\gamma _{12}^2 + {G_3}\gamma _{12}^3 + \cdots }+{G_i}\gamma _{12}^i \\ &\qquad {{\gamma _{12}} < 0} \\ & {{G_1}{\gamma _{12}} - {G_2}\gamma _{12}^2 + {G_3}\gamma _{12}^3 - \cdots }+{G_i}\gamma _{12}^i\\ &\qquad{{\gamma _{12}} \geqslant 0} \end{aligned} \right.$

(5)

式中： ${\tau }_{12}、{\gamma }_{12}$ 分别为剪切应力和工程剪应变； $G_{i}(i=1,2,3 \cdots)$ 为系数，由试验数据拟合给出。如图6所示，取前4项时已足够精确，各系数分别为：G₁=5.02 GPa、G₂=1.48×10² GPa、G₃=2.16×10³ GPa、G₄=1.18×10⁴ GPa。

图 6 CCF300/BA3202复合材料的剪切非线性行为

Figure 6. Shear nonlinear behavior of CCF300/BA3202

下载: 全尺寸图片幻灯片

由此，等效剪切刚度 $G_{12}^{{\text{eq}}}$ ，计算公式如下：

$G_{12}^{{\text{eq}}} = \left\{ \begin{aligned} & {{G_1} + {G_2}{\gamma _{12}} + {G_3}\gamma _{12}^2 + {G_4}\gamma _{12}^3}&{{\gamma _{12}} < 0} \\ & {{G_1} - {G_2}{\gamma _{12}} + {G_3}\gamma _{12}^2 - {G_4}\gamma _{12}^3}&{{\gamma _{12}} \geqslant 0} \end{aligned} \right.$

(6)

2.3 基于软化夹杂法的损伤演化模型

材料失效后的损伤演化行为主要由断裂韧性控制，参考Linde等^[28]提出的指数型损伤演化模型，并引入单元特征长度以降低模型对网格的依赖性，纤维损伤状态变量d_FF及纤维间损伤状态变量d_IFF定义如下：

$\begin{gathered} {d_{{\text{FF}}}} = \left\{ {\begin{array}{*{20}{c}} {1 - \dfrac{{\exp (X_{\text{T}}^2{L_{\text{C}}}(1 - {f_{{\text{FF}}}})/({E_1}{G_{{\text{FT}}}}))}}{{{f_{{\text{FF}}}}}}}&{{\sigma _1} \geqslant 0} \\ {1 - \dfrac{{\exp (X_{\text{C}}^2{L_{\text{C}}}(1 - {f_{{\text{FF}}}})/({E_1}{G_{{\text{FC}}}}))}}{{{f_{{\text{FF}}}}}}}&{{\sigma _1} < 0} \end{array}} \right. \\ {d_{{\text{IFF}}}} = \left\{ {\begin{array}{*{20}{c}} {1 - \dfrac{{\exp (Y_{\text{T}}^2{L_{\text{C}}}(1 - {f_{{\text{IFF}}}})/({E_2}{G_{{\text{MT}}}}))}}{{{f_{{\text{IFF}}}}}}}&{{\sigma _{\rm{n}}} \geqslant 0} \\ {1 - \dfrac{{\exp (Y_{\text{C}}^2{L_{\text{C}}}(1 - {f_{{\text{IFF}}}})/({E_2}{G_{{\text{MC}}}}))}}{{{f_{{\text{IFF}}}}}}}&{{\sigma _{\rm{n}}} < 0} \end{array}} \right. \end{gathered}$

(7)

式中： $L_{\mathrm{C}}$ 为单元特征长度； $G_{\mathrm{FT}}$ 、 $G_{\mathrm{FC}}$ 为纵向拉伸、压缩断裂韧性，目前还没有公认的标准方法进行测试，参考同类材料的性能^[29-30]， $G_{\mathrm{FT}}$ 、 $G_{\mathrm{FC}}$ 分别取为71.4 N·mm⁻¹、35.6 N·mm⁻¹； $G_{\text {MT }}$ 、 $G_{\mathrm{MC}}$ 为横向拉伸、压缩断裂韧性，近似取为层间I型、II型断裂韧性。

纤维失效后，纵向模量E₁、横向模量E₂、泊松比ν₁₂、剪切模量G₁₂折减为无损材料值的(1−d_FF)倍。纤维间失效后，横向模量E₂、泊松比ν₁₂、剪切模量G₁₂折减为无损材料值的(1−d_IFF)倍。

对于冲击损伤区，采用软化夹杂法进行等效。文献[16, 29]根据损伤的严重程度，将冲击损伤区划分为两个区域，并分别采用不同的系数对初始性能进行折减。本文研究的试件冲击能量小，冲击点处检测未发现纤维大量断裂的情况，为简化模型，整个损伤区采用同一系数η进行性能折减。

上述基于软化夹杂法的损伤演化模型材料性能退化方案如表2所示。

表 2 材料性能退化方案

Table 2. Material performance degradation scheme

Fail mode	No damage zone	Impact damage zone
None	—	$\begin{gathered} {E_1} \to \eta {E_1} \\ {E_2} \to \eta {E_2} \\ G_{12}^{{\text{eq}}} \to \eta G_{12}^{{\text{eq}}} \\ \end{gathered}$
Fiber failure	$\begin{gathered} {E_1} \to (1 - {d_{{\text{FF}}}}){E_1} \\ {E_2} \to (1 - {d_{{\text{FF}}}}){E_2} \\ G_{12}^{{\text{eq}}} \to (1 - {d_{{\text{FF}}}})G_{12}^{{\text{eq}}} \\ {\nu _{12}} \to (1 - {d_{{\text{FF}}}}){\nu _{12}} \\ \end{gathered}$	$\begin{gathered} {E_1} \to (1 - {d_{{\text{FF}}}})\eta {E_1} \\ {E_2} \to (1 - {d_{{\text{FF}}}})\eta {E_2} \\ G_{12}^{{\text{eq}}} \to (1 - {d_{{\text{FF}}}})\eta G_{12}^{{\text{eq}}} \\ {\nu _{12}} \to (1 - {d_{{\text{FF}}}}){\nu _{12}} \\ \end{gathered}$
Inter-fiber failure	$\begin{gathered} {E_2} \to (1 - {d_{{\text{IFF}}}}){E_2} \\ G_{12}^{{\text{eq}}} \to (1 - {d_{{\text{IFF}}}})G_{12}^{{\text{eq}}} \\ {\nu _{12}} \to (1 - {d_{{\text{IFF}}}}){\nu _{12}} \\ \end{gathered}$	$\begin{gathered} {E_2} \to (1 - {d_{{\text{IFF}}}})\eta {E_2} \\ G_{12}^{{\text{eq}}} \to (1 - {d_{{\text{IFF}}}})\eta G_{12}^{{\text{eq}}} \\ {\nu _{12}} \to (1 - {d_{{\text{IFF}}}}){\nu _{12}} \\ \end{gathered}$
Notes: $G_{12}^{{\text{eq}}}$ —Equivalent shear modulus; η—Stiffness reduction factor; ${d_{{\text{FF}}}}$ —Fiber damage state variable; ${d_{{\text{IFF}}}}$ —Matrix damage state variable.

下载: 导出CSV

| 显示表格

2.4 考虑冲击损伤的有限元模型

采用商业有限元软件ABAQUS 6.14对含冲击损伤的层合板进行剪切后屈曲失效分析。

有限元模型中，剪切边框和螺栓进行了如下简化：剪切框架简化为S4 R壳单元(厚度10 mm)，并在8个对角点处施加铰接约束；螺栓简化为B31梁单元(直径8 mm)，分别与试件及上下剪切框架绑定在一起；剪切边框和螺栓材料均为低碳钢。低碳钢的弹性模量为200 GPa、泊松比为0.3。

试件也均采用S4 R壳单元模拟，其中加强区上下各增加一层玻璃钢(厚度2.5 mm)。玻璃钢纵向模量为53.48 GPa、横向模量为17.7 GPa、剪切模量为5.83 GPa、泊松比为0.278。复合材料单向板本构模型由自定义场子程序USDFLD实现，其中场变量为F₁~F₄：F₁表征纤维损伤、F₂表征纤维间损伤、F₃表征剪切非线性行为、F₄区分是否为冲击损伤区，F₄=0时表示非冲击损伤区，F₄=1时表示冲击损伤区。F₁~F₃定义如下：

$\begin{gathered} {F_1} = {d_{{\text{FF}}}} \\ {F_2} = {d_{{\text{IFF}}}} \\ {F_3} = 1 - \frac{{G_{12}^{{\text{eq}}}}}{{{G_{12}}}} \end{gathered}$

(8)

对于含冲击损伤的试件，采用软化夹杂法对冲击损伤进行模拟，并根据单元中心坐标确定单元是否位于冲击损伤区域。试件T4的冲击损伤边界近似于椭圆，且有限元模型原点在试件中心，由此F₄可由下式确定：

${F_4} = \left\{ \begin{aligned} & 1&{\frac{{{x^2}}}{{{a^2}}} + \frac{{{y^2}}}{{{b^2}}} < 1} \\ & 0&{\frac{{{x^2}}}{{{a^2}}} + \frac{{{y^2}}}{{{b^2}}} \geqslant 1} \end{aligned} \right.$

(9)

式中：长半轴a=9 mm；短半轴b=6 mm。

有限元模型网格及边界条件如图7所示：剪切框架网格尺寸为10 mm×10 mm，试件网格尺寸为1 mm×1 mm。边界条件设置如下：在对角坐标系(x'Oy')下，约束y'轴方向的4个主节点的面外平动自由度(U_z'=0)，绕x'、y'轴的转动自由度(UR_x'=UR_y'=0)；x'轴方向的4个主节点仅放开绕z'轴的转动自由度(U_x'=U_y'=U_z'=0，UR_x'=UR_y'=0)。加载时+x'轴方向的2个主节点施加沿x'轴方向4 mm位移(U_x'=4 mm)。

图 7 CCF300/BA3202复合材料层板剪切试验有限元模型

Figure 7. Finite element model for shear test of CCF300/BA3202 composite laminates

MPC—Multi-point constraints

下载: 全尺寸图片幻灯片

通过场变量表征的等效冲击损伤区域如图8所示，本文采用的方法保证了网格的整体质量。

图 8 T4试件等效冲击损伤区域

Figure 8. Equivalent impact damage area of T4 specimen

F₄—Impact damage field variable

下载: 全尺寸图片幻灯片

3. 结果与讨论

3.1 试验结果分析

通过粘贴的花片可以测得0°、45°、90°方向的正应变 ${\varepsilon }_{0}、{\varepsilon }_{45}、{\varepsilon }_{90}$ ，由此可以确定测量点的工程剪应变 ${\gamma _{12}}$ ：

${\gamma _{12}} = {\varepsilon _0} + {\varepsilon _{90}} - 2{\varepsilon _{45}}$

(10)

T1试件的载荷-应变曲线如图9所示，其中虚线为拟合的载荷-应变曲线，根据0~40 kN范围内的载荷及应变线性回归得到。加载到64.6 kN以后应变曲线有明显偏转，由此确定试件屈曲载荷约为64.6 kN。其余试件屈曲载荷确定方法与之类似，试验结果统计如表3所示。

图 9 T1试件载荷-应变曲线

Figure 9. Load-strain curve of T1 specimen

下载: 全尺寸图片幻灯片

表 3 CCF300/BA3202复合材料层板剪切试验结果

Table 3. Shear test results of CCF300/BA3202 composite laminates

Number	Buckling load/kN		Fracture load/kN
Number	Test	Average	Test	Average
T1	64.6	64.2	153.5	156.8
T2	65.0		156.2
T3	63.1		160.7
T4	63.0		141.6

下载: 导出CSV

| 显示表格

不含冲击损伤的试件(T1~T3)屈曲载荷平均值为64.2 kN，破坏载荷平均值为156.8 kN。3个试件的屈曲模态一致，包含3个屈曲半波，沿着垂直于对角拉伸方向分布。破坏模式均为靠近对角拉伸区长边加强片处因应力集中而断裂。无损伤的试件(T1)典型屈曲模态及失效照片如图10所示。

图 10 无损伤的T1试件屈曲及失效时的面外位移云图

Figure 10. Out-of-plane displacement cloud diagram of undamaged T1 specimen during buckling and failure

W—Out of plane displacement

下载: 全尺寸图片幻灯片

含冲击损伤的试件(T4)屈曲载荷测试值为63.0 kN，破坏载荷测试值为141.6 kN，分别比无损伤试件平均值低1.92%、9.69%。如图11所示，T4试件屈曲模态与无损伤试件一致，破坏模式为损伤沿着冲击点向4个角点扩展直至完全失去承载能力，与无损伤试件不同。对比分析可知：试件中心引入冲击损伤后，损伤尺寸较小，整体刚度变化不明显，因此屈曲载荷变化不大。冲击点附近因后屈曲变形导致应力集中，承载能力明显降低。

图 11 含冲击损伤的T4试件屈曲及失效时的面外位移云图

Figure 11. Out-of-plane displacement cloud diagram of T4 specimen with impact damage during buckling and failure

下载: 全尺寸图片幻灯片

3.2 仿真结果与试验对比

试件制造过程中的几何缺陷，特别是翘曲对初始屈曲影响较大，因此数值仿真模型需加以考虑。比较常见的做法是引入一阶屈曲模态作为初始几何缺陷^[31]，几何缺陷因子ξ定义为缺陷幅值与试件考核区厚度之比。汪厚冰等^[10]研究指出几何缺陷因子取值对稳定性影响较大，需进行参数分析确定。

图12为几何缺陷因子和剪切非线性对CCF300/BA3202复合材料层板屈曲载荷、破坏载荷的影响。横坐标表示几何缺陷因子，纵坐标分别表示无损伤试件仿真的屈曲载荷与无损伤试件试验平均屈曲载荷之比、仿真的破坏载荷与无损伤试件试验平均破坏载荷之比。

图 12 几何缺陷因子和剪切非线性对CCF300/BA3202复合材料层合板屈曲载荷、破坏载荷的影响

Figure 12. Effects of geometric defect factor and shear nonlinearity on buckling load and failure load of CCF300/BA3202 composite laminates

下载: 全尺寸图片幻灯片

由图12(a)可以看出，线性剪切响应和非线性剪切响应得到的初始屈曲载荷一致，这是由于初始屈曲时，屈曲应变较小，剪切应力-应变关系接近于线性。随着初始几何缺陷因子增加，分析得到的屈曲载荷随之而减小；初始几何缺陷因子ξ=3%时，与试验结果较接近，相对误差为2.34%。

由图12(b)可以看出，考虑线性剪切响应和非线性剪切响应时的破坏载荷相差较大，采用线性模型误差可达20%以上，而采用非线性模型误差均在5%以内。这是由于临近破坏时，试件剪切应变较高，线性模型已难以反映真实的剪切应力-应变关系。而随着初始几何缺陷因子增加，两种模型分析得到的破坏载荷均无明显变化。

基于以上对无损伤试件的研究结果，对含冲击损伤的试件剪切屈曲后屈曲行为进行分析时，采用非线性剪切响应模型，初始几何缺陷因子取为3%。冲击损伤区刚度折减系数η对承载能力的影响如图13所示，破坏载荷随着η的减小而减小，当η在0.2~0.3之间时与试验结果接近。

图 13 刚度折减系数对CCF300/BA3202复合材料层板破坏载荷的影响

Figure 13. Effects of stiffness reduction factor on failure load of CCF300/BA3202 composite laminates

FEM—Finite element modeling

下载: 全尺寸图片幻灯片

η=0.25时，含冲击损伤的试件屈曲载荷、破坏载荷预测值分别为61 kN、139.8 kN，与试验相比，相对误差分别为−3.17%、−1.27%。仿真的屈曲模态、破坏模式见图14、图15，与图11对比看出，数值仿真的屈曲波形、破坏模式均与试验一致，由此表明了本文考虑剪切非线性效应的软化夹杂模型的有效性。

图 14 T4试件仿真的屈曲模态

Figure 14. Simulated buckling mode of T4 specimen

U₃—Out of plane displacement

下载: 全尺寸图片幻灯片

图 15 T4试件仿真的失效模式

Figure 15. Simulated failure mode of T4 specimen

F₁—Fiber damage field variable; F₂—Matrix damage field variable

下载: 全尺寸图片幻灯片

含冲击损伤试件从初始损伤到最终破坏的过程十分短暂，DIC设备未能捕捉到，因此通过数值仿真结果对其失效过程进行分析。T4试件预测的载荷-位移曲线和损伤扩展过程如图16、图17所示：试件加载到61 kN以后开始屈曲，冲击点附近开始往外鼓包，导致应力集中。加载至132 kN后，冲击点附近出现新增损伤。随着继续加载，损伤沿着4个角点缓慢扩展，达到最大载荷139.8 kN后，损伤急剧扩展并迅速掉载，直到完全破坏。这清晰地反映了冲击损伤对损伤起始和扩展的影响。

图 16 T4试件预测的载荷-位移曲线

Figure 16. Predicted load-displacement curve of T4 specimen

下载: 全尺寸图片幻灯片

图 17 T4试件纤维间损伤扩展过程

Figure 17. Inter-fiber failure propagation process of T4 specimen

P—Load; U—Displacement

下载: 全尺寸图片幻灯片

4. 结论

以无损伤和含冲击损伤的CCF300/BA3202碳纤维增强环氧树脂复合材料层合板为研究对象，基于数字图像相关方法(DIC)进行了试验测试，并对试验过程进行了仿真分析，主要结论如下：

(1) 无损伤试件的屈曲载荷、破坏载荷平均值分别为64.2 kN、156.8 kN；引入冲击损伤后，试件屈曲载荷降为63.0 kN，变化不大；承载能力下降至141.6 kN，降幅为9.69%。

(2) 几何缺陷对复合材料层合板的初始剪切屈曲影响显著，本文通过引入一阶屈曲模态模拟获得了与试验一致的初始屈曲载荷(相对误差仅为2.34%)；后续可将实际测量的翘曲变形引入模型中以更准确地模拟结构的初始屈曲行为。

(3) 剪切非线性效应对层合板剪切承载能力影响较大，线性模型预测的失效载荷与试验测试值误差超过20%，而考虑了剪切非线性响应的仿真模型预测精度可达5%以内。

(4) 本文采用的考虑剪切非线性效应的软化夹杂法可有效地预测含冲击损伤的复合材料层合板的剪切强度，屈曲载荷、破坏载荷误差分别为−3.17%、−1.27%，得到的屈曲波形、破坏模式均与DIC测量结果一致，可为复合材料结构冲击损伤简化和剪切强度预测提供一种思路。

图 1 BC和MoS₂/BC复合材料制备的实验流程图

Figure 1. Experimental flow chart of preparation of BC and MoS₂/BC composites

下载: 全尺寸图片幻灯片

图 2 SP (a) 和BC (b) 的SEM图像

Figure 2. SEM images of SP (a) and BC (b)

下载: 全尺寸图片幻灯片

图 3 MoS₂/BC复合材料的SEM图像 ((a)~(c)) 和TEM图像 ((d)~(f))

Figure 3. SEM ((a)-(c)) and TEM images ((d)-(f)) of MoS₂/BC composites

下载: 全尺寸图片幻灯片

图 4 MoS₂/BC-0.4的TEM图像(a)及EDS图像 ((b)~(d))

Figure 4. TEM (a) and EDS images ((b)-(d)) of MoS₂/BC-0.4

下载: 全尺寸图片幻灯片

图 5 SP、MoS₂/水热碳-0.8 (MoS₂/HTC-0.8) 和MoS₂/BC-0.8的FT-IR图谱

Figure 5. FT-IR images of SP, MoS₂/hydrothermal carbon-0.8 (MoS₂/HTC-0.8) and MoS₂/BC-0.8

下载: 全尺寸图片幻灯片

图 6 BC和MoS₂/BC复合材料的XRD图谱

Figure 6. XRD patterns of BC and MoS₂/BC composites

下载: 全尺寸图片幻灯片

图 7 BC和MoS₂/BC复合材料的XPS图谱:(a) 全谱； (b) Mo3d图谱和S2s图谱；(c) S2p图谱

Figure 7. XPS images of BC and MoS₂/BC composites: (a) Full spectrum; (b) Mo3d spectrum and S2s spectrum; (c) S2p spectrum

下载: 全尺寸图片幻灯片

图 8 BC和MoS₂/BC复合材料的复介电常数实部(a)、复介电常数虚部(b)和介电损耗正切值(c)

Figure 8. Real part of complex permittivity (a), imaginary part of complex permittivity (b) and dielectric loss tangent (c) of BC and MoS₂/BC composites

下载: 全尺寸图片幻灯片

图 9 BC和MoS₂/BC复合材料的Cole-Cole曲线

Figure 9. Cole-Cole curves of BC and MoS₂/BC composites

下载: 全尺寸图片幻灯片

图 10 BC和MoS₂/BC复合材料的反射率损失曲线

Figure 10. Reflection loss curves of BC and MoS₂/BC composites

下载: 全尺寸图片幻灯片

图 11 BC和MoS₂/BC复合材料的阻抗匹配图 (a) 和衰减常数图 (b)

Figure 11. Impedance matching (a) and attenuation constant diagram (b) of BC and MoS₂/BC composites

下载: 全尺寸图片幻灯片

表 1 不同钼酸钠、L-半胱氨酸和柚子皮 (SP) 添加量下制备的生物质碳 (BC)、MoS₂/BC复合材料的参数

Table 1 Parameters of biomass carbon (BC) and MoS₂/BC composites prepared with different additions of sodium molybdate, L-cysteine and shaddock peel (SP)

Samples	SP/g	Na₂MoO₄·2H₂O/g	L-cysteine/g
BC	0.25	0	0
Mo₂S/BC-0.4	0.25	0.10	0.22
Mo₂S/BC-0.6	0.25	0.15	0.33
Mo₂S/BC-0.8	0.25	0.20	0.44

下载: 导出CSV

表 2 BC和MoS₂/BC复合材料的电导率

Table 2 Conductivity of BC and MoS₂/BC composites

Sample	BC	MoS₂/BC-0.4	MoS₂/BC-0.6	MoS₂/BC-0.8
Electrical conductivity/(S·cm⁻¹)	2.53×10⁻²	4.95×10⁻³	2.11×10⁻³	1.13×10⁻⁴

下载: 导出CSV

表 3 MoS₂/BC复合材料与文献报道其他吸波材料性能比较

Table 3 Performance comparison between MoS₂/BC composite and other wave absorbing agents reported in the literature

Absorber	Minimum RL (dB)/ Thickness (mm)	Bandwidth (GHz)/ Thickness (mm)	Ref.
3D MoS₂	−10.3/4.0	0.08/4.0	[11]
3D MoS₂/RGO	−49.5/2.4	5.7/2.3	[11]
MoS₂/GN	−55.3/1.6	5.6/2.2	[12]
CNTs@MoS₂	−55.0/1.5	4.0/1.3	[13]
MoS₂/RGO	−50.9/2.3	5.7/2.0	[36]
MoS₂/BC-0.8	−40.1/3.0	5.9/2.0	This work
Notes: RGO—Reduced graphene oxide; GN—Graphene; CNTs—Carbon nanotubes; RL—Reflectivity loss.

下载: 导出CSV

参考文献(39)

[1]	杨亚楠, 夏龙, 张昕宇, 等. Fe₃O₄@锂铝硅微晶玻璃/还原氧化石墨烯复合材料的制备和吸波性能[J]. 复合材料学报, 2019, 36(11):2651-2664. YANG Yanan, XIA Long, ZHANG Xiyu, et al. Preparation and microwave absorbing properties of Fe₃O₄@lithium aluminum silicate glass ceramic/reduced graphene oxide composite[J]. Acta Materiae Compositae Sinica,2019,36(11):2651-2664(in Chinese).
[2]	马志军, 莽昌烨, 翁兴媛, 等. Zn还原氧化石墨烯(RGO)和ZnO/RGO自组装复合材料的电磁响应行为[J]. 复合材料学报, 2019, 36(7):1776-1786. MA Zhijun, MANG Changyue, WENG Xingyuan, et al. Electromagnetic response behavior of Zn reduced graphene oxide (RGO) and ZnO/RGO self-assembled composites[J]. Acta Materiae Compositae Sinica,2019,36(7):1776-1786(in Chinese).
[3]	LV H, YANG Z, WANG P L, et al. A voltage-boosting strategy enabling a low-frequency, flexible electro-magnetic wave absorption device[J]. Advanced Materials,2018,30(15):1706343. DOI: 10.1002/adma.201706343
[4]	CAO M S, WANG X X, ZHANG M, et al. Variable-tempera-ture electron transport and dipole polarization turning flexible multifunctional microsensor beyond electrical and optical energy[J]. Advanced Materials,2020,32(10):1907156. DOI: 10.1002/adma.201907156
[5]	WANG P, WANG G W, ZHANG J M, et al. Excellent microwave absorbing performance of the sandwich structure absorber Fe@B₂O₃/MoS₂/Fe@B₂O₃ in the Ku-band and X-band[J]. Chemical Engineering Journal,2020,382:122804. DOI: 10.1016/j.cej.2019.122804
[6]	FENG Z, YANG P P, WEN G S, et al. One-step synthesis of MoS₂ nanoparticles with different morphologies for electromagnetic wave absorption[J]. Applied Surface Science,2020,502:144129. DOI: 10.1016/j.apsusc.2019.144129
[7]	NING M Q, JIANG P H, DING W, et al. Phase manipulating toward molybdenum disulfide for optimizing electro-magnetic wave absorbing in gigahertz[J]. Advanced Functional Materials,2021,31(19):2011229. DOI: 10.1002/adfm.202011229
[8]	KUMAR P. Ultrathin 2D nanomaterials for electro-magnetic interference shielding[J]. Advanced Materials Interfaces,2019,6(24):1901454. DOI: 10.1002/admi.201901454
[9]	CAO M S, SHU J C, WANG X X, et al. Electronic structure and electromagnetic properties for 2D electromagnetic functional materials in gigahertz frequency[J]. Annalen Der Physik,2019,531(4):1800390. DOI: 10.1002/andp.201800390
[10]	NING M Q, MAN Q K, TAN G G, et al. Ultrathin MoS₂ nanosheets encapsulated in hollow carbon spheres: A case of a dielectric absorber with optimized impedance for efficient microwave absorption[J]. ACS Applied Materials & Interfaces,2020,12(18):20785-20796.
[11]	DING X, HUANG Y, LI S P, et al. 3D Architecture reduced graphene oxide-MoS₂ composite: Preparation and excellent electromagnetic wave absorption performance[J]. Composites Part A: Applied Science and Manufacturing,2016,90:424-432. DOI: 10.1016/j.compositesa.2016.08.006
[12]	ZHANG D Q, JIA Y X, CHENG J Y, et al. High-performance microwave absorption materials based on MoS₂-graphene isomorphic hetero-structures[J]. Journal of Alloys and Compounds,2018,758:62-71. DOI: 10.1016/j.jallcom.2018.05.130
[13]	WANG R, YANG E Q, QI X S, et al. Constructing and opti-mizing core@shell structure CNTs@MoS₂ nanocomposites as outstanding microwave absorbers[J]. Applied Surface Science,2020,516:146159. DOI: 10.1016/j.apsusc.2020.146159
[14]	ZHANG Y, HUANG Y, ZHANG T F, et al. Broadband and tunable high-performance microwave absorption of an ultralight and highly compressible graphene foam[J]. Advanced Materials,2015,27(12):2049-2053. DOI: 10.1002/adma.201405788
[15]	LIU Y, CHEN Z, ZHANG Y, et al. Broadband and lightweight microwave absorber constructed by in situ growth of hierarchical CoFe₂O₄/reduced graphene oxide porous nanocomposites[J]. ACS Applied Materials & Interfaces,2018,10(16):13860-13868.
[16]	ZHAO J, ZHANG J L, WANG L, et al. Fabrication and investi-gation on ternary heterogeneous MWCNT@TiO₂-C fillers and their silicone rubber wave-absorbing composites[J]. Composites Part A: Applied Science and Manufacturing,2020,129:105714. DOI: 10.1016/j.compositesa.2019.105714
[17]	LIU Y, CHEN Z, XIE W H, et al. In-situ growth and graphiti-zation synthesis of porous Fe₃O₄/carbon fiber composites derived from biomass as lightweight microwave absorber[J]. ACS Sustainable Chemistry & Engineering,2019,7(5):5318-5328.
[18]	ZHOU X F, JIA Z R, FENG A L, et al. Construction of multiple electromagnetic loss mechanism for enhanced electromagnetic absorption performance of fish scale-derived biomass absorber[J]. Composites Part B: Engi-neering,2020,192:107980. DOI: 10.1016/j.compositesb.2020.107980
[19]	TOCMO R, PENA-FRONTERAS J, CALUMBA K F, et al. Valorization of pomelo peel: A review of current utilization, phytochemistry, bioactivities, and mechanisms of action[J]. Comprehensive Reviews in Food Science and Food Safety,2020,19(4):1969-2012. DOI: 10.1111/1541-4337.12561
[20]	LIU J Y, LI H P, ZHANG H S, et al. Three-dimensional hierarchical and interconnected honeycomb-like porous carbon derived from pomelo peel for high performance supercapacitors[J]. Journal of Solid State Chemistry,2018,257:64-71. DOI: 10.1016/j.jssc.2017.07.033
[21]	HU B, WANG K, WU L H, et al. Engineering carbon materials from the hydrothermal carbonization process of biomass[J]. Advanced Materials,2010,22(7):813-828. DOI: 10.1002/adma.200902812
[22]	WANG Q, LI H, CHEN L Q, et al. Monodispersed hard carbon spherules with uniform nanopores[J]. Carbon,2001,39(14):2211-2214. DOI: 10.1016/S0008-6223(01)00040-9
[23]	SEVILLA M, LOTA G, FUERTES A B. Saccharide-based graphitic carbon nanocoils as supports for PtRu nanoparticles for methanol electrooxidation[J]. Journal of Power Sources,2007,171(2):546-551. DOI: 10.1016/j.jpowsour.2007.05.096
[24]	ZHANG Z, TAN J W, GU W H, et al. Cellulose-chitosan framework/polyailine hybrid aerogel toward thermal insulation and microwave absorbing application[J]. Chemical Engineering Journal, 2020, 395: 125190.
[25]	YANG E Q, QI X S, XIE R, et al. Novel "203" type of heterostructured MoS₂-Fe₃O₄-C ternary nanohybrid: Synthesis, and enhanced microwave absorption properties[J]. Applied Surface Science,2018,442:622-629. DOI: 10.1016/j.apsusc.2018.02.175
[26]	SUN X, LI Y. Colloidal carbon spheres and their core/shell structures with noble-metal nanoparticles[J]. Angewandte Chemie International Edition,2004,43(5):597-601. DOI: 10.1002/anie.200352386
[27]	DONG N, HE F Z, XIN J L, et al. A novel one-step hydrothermal method to prepare CoFe₂O₄/graphene-like carbons magnetic separable adsorbent[J]. Materials Research Bulletin,2016,80:186-190. DOI: 10.1016/j.materresbull.2016.04.003
[28]	SEVILLA M, FUERTES A B. Chemical and structural properties of carbonaceous products obtained by hydrothermal carbonization of saccharides[J]. Chemistry-A European Journal,2009,15(16):4195-4203. DOI: 10.1002/chem.200802097
[29]	SEMERCIOZ A S, GOGUS F, ÇELEKLI A, et al. Development of carbonaceous material from grapefruit peel with microwave implemented-low temperature hydrothermal carbonization technique for the adsorption of Cu (II)[J]. Journal of Cleaner Production,2017,165:599-610. DOI: 10.1016/j.jclepro.2017.07.159
[30]	VOLPE M, MESSINEO A, MAKELA M, et al. Reactivity of cellulose during hydrothermal carbonization of lignocellulosic biomass[J]. Fuel Processing Technology,2020,206:106456. DOI: 10.1016/j.fuproc.2020.106456
[31]	ZHAO H Q, CHENG Y, LIU W, et al. Biomass-derived porous carbon-based nanostructures for microwave absorption[J]. Nano-Micro Letters,2019,11(1):81-97. DOI: 10.1007/s40820-019-0312-y
[32]	ZHANG W L, JIANG D G, WANG X X, et al. Growth of polyaniline nanoneedles on MoS₂ nanosheets, tunable electroresponse, and electromagnetic wave attenuation analysis[J]. Journal of Physical Chemistry C,2017,121(9):4989-4998. DOI: 10.1021/acs.jpcc.6b11656
[33]	WANG Y L, YANG S H, WANG H Y, et al. Hollow porous CoNi/C composite nanomaterials derived from MOFs for efficient and lightweight electromagnetic wave absorber[J]. Carbon,2020,167:485-494. DOI: 10.1016/j.carbon.2020.06.014
[34]	ZHAO H Q, CHENG Y, ZHANG Z, et al. Biomass-derived graphene-like porous carbon nanosheets towards ultralight microwave absorption and excellent thermal infrared properties[J]. Carbon,2021,173:501-511. DOI: 10.1016/j.carbon.2020.11.035
[35]	WU Z C, TIAN K, HUANG T, et al. Hierarchically porous carbons derived from biomasses with excellent microwave absorption performance[J]. ACS Applied Materials & Interfaces,2018,10(13):11108-11115.
[36]	WANG Y, DI X C, WU X M, et al. MOF-derived nanoporous carbon/Co/Co₃O₄/CNTs/RGO composite with hierarchi-cal structure as a high-efficiency electromagnetic wave absorber[J]. Journal of Alloys and Compounds,2020,846:156215. DOI: 10.1016/j.jallcom.2020.156215
[37]	SHU R W, ZHANG G Y, ZHANG C, et al. Nitrogen-doping-regulated electromagnetic wave absorption properties of ultralight three-dimensional porous reduced graphene oxide aerogels[J]. Advanced Electronic Materials,2020,7(2):2001001.
[38]	WANG Y F, CHEN D L, YIN X, et al. Hybrid of MoS₂ and reduced graphene oxide: A lightweight and broadband electromagnetic wave absorber[J]. ACS Applied Materials & Interfaces,2015,7(47):26226-26234.
[39]	MO Z C, YANG R L, LU D W, et al. Lightweight, three-dimensional carbon nanotube@TiO₂ sponge with enhanced microwave absorption performance[J]. Carbon,2019,144:433-439. DOI: 10.1016/j.carbon.2018.12.064