Loading [MathJax]/jax/output/SVG/jax.js

连续碳纤维3D打印圆形增强蜂窝的面内压缩性能

孟云聪, 周光明, 蔡登安

孟云聪, 周光明, 蔡登安. 连续碳纤维3D打印圆形增强蜂窝的面内压缩性能[J]. 复合材料学报, 2024, 41(4): 1776-1787. DOI: 10.13801/j.cnki.fhclxb.20230721.001
引用本文: 孟云聪, 周光明, 蔡登安. 连续碳纤维3D打印圆形增强蜂窝的面内压缩性能[J]. 复合材料学报, 2024, 41(4): 1776-1787. DOI: 10.13801/j.cnki.fhclxb.20230721.001
MENG Yuncong, ZHOU Guangming, CAI Deng'an. In-plane compression properties of 3D printed continuous carbon fiber circular improved honeycomb[J]. Acta Materiae Compositae Sinica, 2024, 41(4): 1776-1787. DOI: 10.13801/j.cnki.fhclxb.20230721.001
Citation: MENG Yuncong, ZHOU Guangming, CAI Deng'an. In-plane compression properties of 3D printed continuous carbon fiber circular improved honeycomb[J]. Acta Materiae Compositae Sinica, 2024, 41(4): 1776-1787. DOI: 10.13801/j.cnki.fhclxb.20230721.001

连续碳纤维3D打印圆形增强蜂窝的面内压缩性能

详细信息
    通讯作者:

    周光明,博士,教授,博士生导师,研究方向为先进复合材料的结构设计 E-mail: zhougm@nuaa.edu.cn

  • 中图分类号: TB332

In-plane compression properties of 3D printed continuous carbon fiber circular improved honeycomb

  • 摘要: 为提高圆形蜂窝(CH)的抗压缩性能和吸能性能,以CH结构为基础,在横向和竖向上增设树叶形支撑,提出了单向增强圆形蜂窝(SEH)和双向增强圆形蜂窝(DEH)两种改进型蜂窝。以碳纤维(CF)作为增强体,聚乳酸(PLA)为基体,使用连续纤维3D打印技术制造了试验件,并规划结构内部CF的成型路径,同时设置PLA对照组。通过准静态压缩试验研究各蜂窝的面内压缩性能、吸能特性和结构的变形失效模式。结果表明:CF增强后的DEH-CF相较CH-CF在比吸能上提升167.63%。CH、SEH和DEH在采用CF增强后,比吸能相比PLA对照组分别提高43.37%、63.17%和161.58%,平均压缩力分别提高51.72%、61.81%和96.09%。研究发现,CF增强结构内部的纤维路径规划会影响结构的刚度和变形行为,采用“支撑一体化成型”路径的DEH-CF在压缩时,其结构动态泊松比保持在PLA对照组33.36%以下。

     

    Abstract: To improve the compression resistance and energy absorption performance of circular honeycomb (CH), two improved honeycomb, single enhanced circular honeycomb (SEH) and double enhanced circular honeycomb (DEH) were designed on the basis of CH structure, and leaf shaped supports were added horizontally and vertically. Using carbon fiber (CF) as reinforcement and polylactic acid (PLA) as matrix, continuous fiber 3D printing technology was used to manufacture test parts, and the forming path of CF bundle inside the structure was designed, PLA control group was set. The in-plane compression properties, energy absorption characteristics and deformation failure modes of the honeycomb structures were investigated by quasi-static compression tests. The results show that the specific energy absorption (SEA) of CF enhanced DEH-CF is improved by 167.63% compared with CH-CF. The SEA are increased by 43.37%, 63.17% and 161.58% and mean crushing force are increased by 51.72%, 61.81% and 96.09% compared with the PLA control group, respectively. The results indicate that the fiber path planning inside the CF reinforced structure would affect the stiffness and deformation behavior of the structure. The dynamic Poisson's ratio of the DEH-CF using the "strut integrated molding path" during compression remains 33.36% lower than that of the PLA control group.

     

  • 近年来,碳化硅纤维(SiCf)由于电性能可调范围广、耐高温(≥1600℃)、高强高模等优点,由其增强的SiCf/SiC复合材料已成为新型空天高温隐身部件的理想候选材料[1-4]。目前,用于SiCf/SiC吸波复合材料的细观特征主要为二维层合结构。然而,该结构不仅吸波频带窄且可设计性差,还存在易分层和面外承载性能低问题,难以满足新型装备的结构-功能一体化苛刻服役要求[5-7]。相比之下,2.5D机织结构通过面内和面外不同电磁参数纤维同时排布,可实现高承载和优异吸波设计,为高可靠、多功能、高承载高温隐身结构部件提供了有效途径[8-10]

    近年来,学者们将吸波特性与三维纺织结构相结合,开展了诸多实验工作。Mu等[11]首先对SiCf进行脱碳处理,然后采用先驱体浸渍裂解(PIP)工艺制备了2.5D SiCf/SiC复合材料,测得其最大反射损耗可达−23.5 dB,有效吸波带宽可扩展至5.9 GHz。Han等[12]采用PIP工艺制备了三维四向和三维五向编织SiCf/SiC复合材料,并比较了两种结构在X波段的吸波性能。结果表明,由于三维五向结构中轴向纱的存在,致使三维五向编织 SiCf/SiC 复合材料呈现出更优异的微波吸收性能。吕丽华等[13]和Zhang等[14]以羰基铁粉 (CIP) 和炭黑 (CB) 为吸波剂,制备了蜂窝状三维机织结构型吸波复合材料。结果显示,该材料在高频 15.48 GHz 下的吸波损耗为−15.06 dB,且有效吸收带宽大于 3 GHz。Chen等[15]采用纬编工艺制备了菱形针织柔性吸波复合材料,测得其有效吸波带宽可达13 GHz。李丹丹[16]研究了三向正交、经向增强角联锁和层层角联锁结构碳纤维复合材料的吸波性能。结果表明,三向正交结构的吸波性能最好,层层角联锁结构的吸波性能相对最差,但3种结构吸波性能差异不明显。综上,目前所选用的三维纺织结构仍局限于传统的几种类型,导致吸波性能调节幅度有限。同时,受限于昂贵成本,有关2.5D 机织SiCf /SiC 复合材料的吸波性能研究尚处于探索阶段。

    除此之外,电磁数值仿真也是研究三维纺织结构材料吸波性能的重要手段[17-19]。苏钦城等[20]和杨飙[21]利用基于有限积分法的CST软件,对二维平纹织物吸波性能进行了数值分析。孙娜等[22]将二维机织结构织物简化为网格模型,利用CST仿真分析了30~1500 MHz 频段内织物的屏蔽效能,并通过实验验证了仿真的准确性。Yin等[23]设计制备了一种具有双向角联锁结构的3D电磁吸波机织织物,并采用实验与仿真相结合的方法研究了吸波性能。结果表明,仿真与实验结果基本一致,该3D织物具有吸波频率宽、吸波性能好(−20 dB)的优点。Wang等[24]设计制备了双层蜂窝状三维编织吸波复合材料,将CST电磁模拟与实验分析相结合揭示了其电磁波吸收性能。结果表明,仿真和实验得到的反射峰位置分别为10.78 GHz和9.28 GHz,最大反射损耗分别为−24.35 dB和−23.5 dB,二者吻合良好。然而,围绕纺织复合材料的电磁仿真,学者们将细观结构模型进行了大幅简化,未准确体现纱线的交织特性,这也导致其预测数据的准确性和可靠性不足。

    本文选用国产SiC纤维,采用PIP工艺制备了2.5D机织SiCf/SiC复合材料,利用弓形法开展反射损耗测试,获取了其吸收峰位置和有效吸波带宽。在此基础上,建立了基于X射线计算机断层扫描(Micro-CT)图像的全厚度细观模型,利用CST电磁仿真软件开展了数值分析,并与上述实验结果进行了对比验证。之后,进一步揭示了2.5D机织SiCf/SiC复合材料的吸波机制,分析了几何特性、电磁参数、电磁波电场极化方向和入射角度对材料吸波性能的影响规律,为具有优异吸波特性的三维纺织超材料结构设计和制备提供了方法和数据支撑。

    图1给出了2.5D机织SiCf/SiC复合材料的制备工艺流程。首先,选用规格为0.5 K的国产SiC纤维(福建立亚新材有限公司商用的cansas3301 SiC纤维),并利用多综眼织机,织造了2.5D机织预制体;其次,以BCl3和NH3作为原料气(山东德芯特种气体有限责任公司),以Ar (北京普莱克斯实用气体有限公司)作为保护气,H2 (北京普莱克斯实用气体有限公司)作为催化气,通过物理气相沉积法制备了BN界面层[25];然后,以聚碳硅烷(PCS,中国科学院化学研究所)浆料为先驱体,采用PIP法制备了2.5D机织SiCf/SiC复合材料粗坯,浸渍-裂解循环为8~10次,直至质量增加<1.5wt%。最后,经过近净尺寸加工和再致密化等后处理得到2.5D机织SiCf/SiC复合材料样件。制备各阶段实物如图2所示,材料样件规格参数如表1所示。

    参照标准《雷达吸波材料反射率测试方法》(GJB 2038—94)[26],利用弓形法对2.5D机织SiCf/SiC复合材料的反射损耗进行测试。测试系统如图3所示:(1) 测试之前,需在弓形架上方安装好对称的发射天线和接收天线,并将2.5D机织SiCf/SiC复合材料测试试样放置在样品台上;(2) 开始测试后,电磁波由发射天线发射,经样品台反射到接收天线并传输到矢量网络分析仪上,建立起相应特性频率的参考功率电平,测出接收到的反射电磁波功率;(3) 用金属板取代测试试样,获取其反射电磁波的功率。上述金属板与测试试样接收的功率电平之差,即为反射率。本实验在中国航空技术研究院复合材料中心完成,测试试样尺寸为40 mm×40 mm×3 mm,测试频率范围为8~18 GHz。

    图  1  先驱体浸渍裂解(PIP)法制备2.5D机织SiCf/SiC复合材料工艺流程图
    Figure  1.  Process flow chart of 2.5D woven SiCf/SiC composites prepared by precursor infiltration and pyrolysis (PIP)
    PCS—Polycarbosilane
    图  2  (a) SiC纤维预制体;(b) 带BN界面层的SiC纤维预制体;(c) 2.5D机织SiCf/SiC复合材料粗坯;(d) 2.5D机织SiCf/SiC复合材料样件
    Figure  2.  (a) SiC fiber preform; (b) SiC fiber preform with BN interface layer; (c) Rough blank of 2.5D woven SiCf/SiC composites; (d) Sample of 2.5D woven SiCf/SiC composites
    表  1  2.5D机织SiCf/SiC复合材料规格参数
    Table  1.  Specification parameters of 2.5D woven SiCf/SiC composites
    Prefabricated structure Dimension/
    mm3
    Warp density/
    (yarn·cm−1)
    Weft density/
    (yarn·cm−1)
    Yarn fineness/tex Number of
    layer
    Volume
    fraction/vol%
    2.5D woven angle interlock 40×40×3 9 3.5 185 8 45
    下载: 导出CSV 
    | 显示表格

    使用德国Diondo-d2微焦点CT扫描系统扫描2.5D机织SiCf/SiC复合材料,扫描规格为20 mm×20 mm×3 mm。所得复合材料经向截面扫描形态和纬向截面扫描形态,分别如图4(a)图4(b)所示。可知,其纱线截面可近似看作椭圆形或扁六边形。统计分析扫描结果,得到纱线截面几何参数如图4(c)图4(d)所示。

    图  3  弓形法测2.5D机织SiCf/SiC复合材料反射损耗系统示意图
    Figure  3.  System schematic diagram of measuring reflection loss of 2.5D woven SiCf/SiC composites by bow method

    基于图4所示纱线截面几何参数和表1所示复合材料规格参数,在SOLIDWORKS 2018对2.5D机织结构预制体进行全厚度细观模型重构。如图5所示,分别建立了以椭圆为纱线截面和以扁六边形为纱线截面的2.5D机织 SiCf/SiC 复合材料模型。仿真发现,两种不同截面的模型对仿真结果几乎没有影响。但从仿真速度来看,以扁六边形为截面所建模型远快于以椭圆为截面所建模型。这是由于在CST电磁仿真软件中,当网格无法分辨较精细曲面结构时,就会出现CST网格划分短路现象,称为CC短路现象(Critical cell),需要进行全局网格优化或者局部加密网格来解决问题。这导致以椭圆为截面所建模型划分网格数量要高于以六边形为截面所建模型划分网格数量。因此,后续研究以扁六边形为纱线截面所建模型为准。

    图  4  2.5D机织SiCf/SiC复合材料的X射线计算机断层扫描(Micro-CT)扫描结果:(a) 经向截面扫描形态;(b) 纬向截面扫描形态;(c) 椭圆截面几何参数;(d) 扁六边形截面几何参数
    Figure  4.  X-ray computed tomography (Micro-CT) scanning results of 2.5D woven SiCf/SiC composites: (a) Scanning morphology of warp cross-section; (b) Scanning morphology of weft cross-section; (c) Geometric parameters of ellipse cross-section; (d) Geometric parameters of flat hexagonal cross-section
    图  5  2.5D机织SiCf/SiC复合材料模型:(a) 以椭圆为纱线截面;(b) 以扁六边形为纱线截面
    Figure  5.  2.5D woven SiCf/SiC composites model: (a) Cross-section of yarn is ellipse; (b) Cross-section of yarn is flat hexagonal

    对于吸波材料,一般将A(ω)、R(ω)和T(ω)分别定义为吸收率、反射率和透射率。在CST Studio Suite 2020中,反射率R(ω)=|S11|2,透射率T(ω)=|S21|2,此时吸收率A(ω)如下式所示[27]

    A(ω)=1R(ω)T(ω)=1|S11|2|S21|2 (1)

    其中:S11为反射系数;S21为透射系数。在仿真过程中,将边界条件中的Z轴负方向的Zmin设置为电边界,这可以防止入射电磁波传输透过底板,因此透射率为0,此时吸收率A(ω)如下式所示:

    A(ω)=1R(ω)=1|S11|2 (2)

    合格的吸波材料要求吸收率尽可能高,反射率尽可能低,而反射率由等效相对阻抗决定,如下式所示:

    R(ω)=Zeffη0Zeff+η0 (3)

    将式(3)代入式(2),此时吸收率A(ω)如下式所示:

    A(ω)=2η0Re(Zeff)+Im(Zeff)+η0 (4)

    由式(4)可知,若要吸收率A(ω)接近1,则吸波材料的等效阻抗实部Re(Zeff)应接近空气阻抗η0,而虚部Im(Zeff)应接近0。由S参数反演法[28-29]计算等效相对阻抗如下式所示:

    Zeff=(1+S11)2S212(1S11)2S212 (5)

    图6示出了利用 CST 微波工作室分析2.5D 机织SiCf/SiC复合材料吸波性能的仿真流程。在CST Studio Suite 2020中,仿真频率设置为1~18 GHz、边界条件Zmin设置为电边界(无透射)、网格划分为四面体自适应性网格划分、求解器设置为频域求解器、材料参数为介电常数和损耗角正切(波导法测得,由中国航空技术研究院复合材料中心测试提供)。运行仿真后在1D Results中查看反射损耗曲线,在2D/3D Results中查看场分布图。仿真得到反射损耗曲线S11S21后,在Post-processing模块的Result Templates Tools中键入式(5),重新运行仿真,在tables中查看等效电磁参数曲线。

    图  6  CST 微波工作室仿真流程图
    Figure  6.  Simulation flow chart of CST microwave studio

    图7对比分析了实验测试与仿真计算所得的反射损耗曲线。从图中可以看出,室温下实验测得的吸收峰出现在9.3 GHz处,吸收峰处反射损耗为−17 dB,反射损耗小于−10 dB的有效吸波带宽为3 GHz。而对应的仿真结果显示吸收峰出现在10.1 GHz处,吸收峰处反射损耗为−16 dB,有效吸波带宽也为3 GHz。且在高温800℃下,2.5D 机织SiCf/SiC复合材料仍然具有一定的吸波性能。显然,在8~18 GHz频率范围内,实验与仿真结果吻合良好,证明了所建模型的合理性和仿真结果的准确性,这为后续利用仿真模拟进行吸波机制分析的可靠性提供了有力支撑。但也不可否认的是,二者在吸收峰处略有差异,其中主要原因包括:(1) 仿真发生在理想条件下,而实验会受到材料性能差异、测量误差和环境等条件的真实影响;(2) 所建模型没有充分考虑材料内部孔隙和纱线挤压问题,这些因素会在后续研究中展开分析。

    图  7  2.5D机织SiCf/SiC复合材料反射损耗的仿真结果与实验结果对比
    Figure  7.  Comparison between simulation results and experimental results of reflection loss of 2.5D woven SiCf/SiC composites model

    图8给出了2.5D 机织SiCf/SiC复合材料的等效电磁参数曲线。由图可知,在10.1 GHz处,Re(Zeff)和Im(Zeff)分别为1和0。此时,复合材料表面阻抗与空气阻抗近似相等,这能保证电磁波尽可能多的进入材料内部被吸收。而在10.1 GHz两侧,Re(Zeff)和Im(Zeff)分别越来越远离1和0,这导致材料表面阻抗与空气阻抗失配,复合材料产生了对电磁波的强反射作用。这也验证了图7所示反射损耗曲线中,吸收峰出现在10.1 GHz附近,吸收峰向两边的反射损耗越来越小,复合材料吸波性能越来越差。

    图  8  2.5D机织SiCf/SiC复合材料的等效电磁参数曲线
    Figure  8.  Curves of equivalent electromagnetic parameters of 2.5D woven SiCf/SiC composites
    Re(Zeff)—Real part of the equivalent impedance; Im(Zeff)—Imaginary part of the equivalent impedance

    为了进一步研究2.5D机织SiCf/SiC复合材料与电磁波之间的相互作用,揭示其吸波机制, 开展了基于CST电磁仿真的场分布图研究。在频率1~18 GHz之间,分别仿真了频率两端1 GHz、18 GHz、中间点的频率9.5 GHz以及谐振点频率10.1 GHz共4个频点的电磁场分布。

    图9(a)给出了电磁波在2.5D机织SiCf/SiC复合材料和空气中传输的电场分布图。从图中可以看出,当电磁波频率为1 GHz时,电场强度沿Y轴负方向分布,且在整个仿真空间内电场强度方向无变化;当电磁波频率升高至9.5 GHz时,吸波材料上方空气中场强沿Y轴负方向分布,而吸波材料内部场强转为沿Y轴正方向分布;当电磁波频率继续上升至18 GHz时,电场强度的方向发生了多次转变。这与苏钦城[30]的研究结果一致,即不同频率下的电场方向随着电磁波频率的变化而变化,且电磁波频率越高,电场方向改变越快。除此之外,电磁波频率为1 GHz时的最大电场强度为2938 V/m,随着电磁波频率升高至9.5 GHz,最大电场强度迅速升高至29523 V/m,谐振频率点10.1 GHz处的最大电磁强度为29361 V/m。而随着电磁波频率继续升高,最大电场强度开始逐渐下降,当电磁波频率为18 GHz时,最大电场强度下降至16669 V/m。可以看出,本文所制备的2.5D机织SiCf/SiC复合材料主要通过电损耗的方式吸收电磁波,且电场强度越大,损耗能力越强。又由于最大电场强度发生在9.5 GHz,而不是谐振频率10.1 GHz处,可以判断该复合材料对电磁波的损耗作用除了电损耗之外,还存在其他的损耗类型。

    图  9  2.5D机织SiCf/SiC复合材料的场分布图:(a) 电场(E-Field)分布图;(b) 磁场(H-Field)分布图;(c) 表面电流密度分布图;(d) 能量损耗密度分布图
    Figure  9.  Field distribution map of 2.5D woven SiCf/SiC composites: (a) Electric field (E-Field) distribution map; (b) Magnetic field (H-Field) distribution map; (c) Current density distribution map; (d) Power loss density distribution map

    图9(b)给出了电磁波在2.5D机织SiCf/SiC复合材料和空气中传输的磁场分布图。由图可知,当电磁波频率为1 GHz时,整个仿真空间内磁场强度沿X轴负方向分布。而随着电磁波频率的升高,磁场方向逐渐发生变化,且频率越高变化越快。频率为1 GHz时的最大磁场强度为34.9 A/m,随着电磁波频率升高至9.5 GHz时,最大磁场强度升高至55.7 A/m,谐振频率点10.1 GHz处的最大磁场强度为53.6 A/m。而随着电磁波频率继续升高,最大磁场强度开始逐渐下降,当电磁波频率为18 GHz时,最大磁场强度下降至31 A/m。可以看出,在1~18 GHz整个频段内,电磁波形成的磁场强度较弱,2.5D机织SiCf/SiC复合材料由磁损耗贡献的吸波效果较小。这是由于SiC纤维本身磁性较弱,难以在复合材料中形成强磁场。

    图9(c)给出了2.5D机织SiCf/SiC复合材料内部的电流密度分布图。可以看出,电磁波频率为1 GHz时,材料内部电流均沿Y轴负方向分布,这与1 GHz处的电场方向一致。同样的,当电磁波频率为9.5 GHz、10.1 GHz和 18 GHz时,吸波材料内部电流方向均与对应频点处的电场方向一致,沿Y轴正方向分布。这是由于在电场的作用下,大量自由电子沿电场方向发生定向移动产生了电流。当电磁波频率为1 GHz时,材料内部电流密度仅为393 A/m2。随着电磁波频率升高至9.5 GHz时,电流密度迅速升高至31909 A/m2,谐振频率点10.1 GHz处的电流密度为32690 A/m2,此时2.5D机织SiCf/SiC复合材料吸波性能达到最佳。

    图9(d)给出了1~18 GHz间的能量损耗密度分布图。由图可知,当电磁波频率为1 GHz时,能量损耗密度为3.03×105 W/m3,此时复合材料反射损耗绝对值为全频段最小的0.004 dB;随着电磁波频率的上升,能量损耗密度逐渐增大,在谐振点频率10.1 GHz处达到最大值2.66×108 W/m3,此时对应的反射损耗绝对值为全频段最大的16.2 GHz;随着电磁波频率继续升高,能量损耗密度逐渐减小,在18 GHz时能量损耗密度减小为1.44×108 W/m3,此时对应的反射损耗绝对值也减小至3.5 dB。因此,电磁波在2.5D机织SiCf/SiC复合材料内部以能量的形式进行传播,在传播过程中,复合材料将利用各种损耗机制将其转化为热量耗散掉。

    为了更完整地表达2.5D机织SiCf/SiC复合材料的吸波机制,图10给出了电磁波在材料内部的传输路径。由图可知,入射电磁波在接触到复合材料表面时,一部分直接从表面反射回去,一部分穿过表面透波层进入复合材料内部的损耗层。在损耗层中,一方面,由于SiC纤维固有的材料属性,入射波在材料内部会发生电损耗而被吸收;另一方面,采用PIP法制备的2.5D机织SiCf/SiC复合材料内部存在较多随机分布的孔隙(Micro-CT扫描测得该2.5D机织SiCf/SiC复合材料样件内部孔隙率为10.7%),这些孔隙大大延长了电磁波在材料内部的传播路径,并使电磁波在材料内部发生多次散射,这进一步提高了其损耗能力。

    图  10  2.5D机织SiCf/SiC复合材料的吸波机制示意图
    Figure  10.  Schematic diagram of microwave absorbing mechanism of 2.5D woven SiCf/SiC composites

    在上述数值模型与实验结果对比分析基础上,本部分利用有限积分法和CST Studio Suite 2020仿真软件,着重分析几何特性、纤维电磁参数、电磁波电场极化方向和入射角度等对2.5D 机织SiCf/SiC复合材料吸波性能的影响规律,以期为耐高温结构吸波复合材料的选材与结构拓展提供设计依据。

    图11(a)给出了厚度d为1.5 mm、3 mm和6 mm时2.5D机织SiCf/SiC复合材料的反射损耗曲线。可以看出,随着厚度的增加,有效吸波带宽无变化,但材料吸收峰明显向低频移动。同时,最大反射损耗随着厚度增加逐渐增大。且当厚度为6 mm时,在低频和高频处都出现了一个明显的吸收峰。这可以用1/4波长理论[31]来解释,厚度的增加延长了电磁波在材料内部的传播路径,使电磁波尽多被吸收。因此,材料厚度主要影响吸收峰的位置和个数,且在允许的范围内,厚度越大对吸波效能越有利。

    图  11  材料几何特性对2.5D机织SiCf/SiC复合材料吸波性能的影响:(a) 厚度d;(b) 纬密;(c) 纤维层数
    Figure  11.  Influence of geometric properties of materials on microwave absorbing properties of 2.5D woven SiCf/SiC composites: (a) Thickness d; (b) Weft density; (c) Number of fiber layers

    图11(b)给出了纬密为2束/cm、3.5束/cm和5束/cm时2.5D机织SiCf/SiC复合材料的反射损耗曲线。由图可知,随着纬密的增大,材料吸收峰逐渐向低频移动,材料对电磁波的损耗能力逐渐增强,但有效吸收带宽无明显改善。这是由于纬密越大,高介电的SiC纤维所占体积分数越大,复合材料对电磁波的损耗能力越强。同样的,仿真分析经密对材料吸波性能的影响,结果同纬密,这与苏钦城等[20]的研究结论一致。

    图11(c)给出了当厚度为3 mm,纤维层数为6层、8层和10层时2.5D机织SiCf/SiC复合材料的反射损耗曲线。可以看出,在一定厚度下,纤维层数越多,材料对电磁波的损耗能力越强,且吸收峰逐渐向低频移动,这与莫红松等[32]的研究结果一致。这是由于纤维层数越多,电磁波进入材料内部后在各层之间的散射次数越多,电磁波能发生多次损耗。

    图12(a)给出了介电常数实部(ε')变化对2.5D机织SiCf/SiC复合材料反射损耗的影响规律。从图中可以看出,随着纤维ε'的增大,材料吸收峰明显向低频移动,最大反射损耗逐渐增大,但有效吸波带宽无明显变化。因此,纤维ε'主要影响2.5D机织SiCf/SiC复合材料的吸收峰位置和最大反射损耗值,而对有效吸波带宽影响较小。兰州大学应用磁学研究所研究表明[33],吸收峰的出现需要两种机制共同作用,一是材料具有较强的电磁损耗能力;二是吸波材料表面与电磁波产生相位相消作用。随着ε'的增大,材料特征阻抗与空气阻抗之间具有良好匹配特性的频点向低频移动,因此吸收峰逐渐向低频移动。

    图  12  (a) 纤维介电常数实部ε'对吸波性能的影响;(b) 纤维介电损耗角正切tanδ对吸波性能的影响;(c) 材料厚度为3 mm时纤维混杂层数对吸波性能的影响;(d) 材料厚度为6 mm时纤维混杂层数对吸波性能的影响
    Figure  12.  (a) Influence of the fiber real part of permittivity ε' on the absorption performance; (b) Influence of fiber dielectric loss angle tangent tanδ on the absorption performance; (c) Influence of fiber hybrid layers on the absorption properties when the material thickness is 3 mm; (d) Influence of fiber hybrid layers on the absorption properties when the material thickness is 6 mm

    图12(b)给出了介电损耗角正切(tanδ)变化对2.5D机织SiCf/SiC复合材料反射损耗的影响。可以看出,随着纤维tanδ的增大,材料有效吸波带宽和最大反射损耗都逐渐增大,但吸收峰位置无明显变化。不难得知,纤维tanδ主要影响2.5D机织SiCf/SiC复合材料的有效吸波带宽和最大反射损耗,而对吸收峰位置几乎没有影响。

    为了进一步研究纤维电磁参数对2.5D机织SiCf/SiC复合材料吸波性能的影响,开展了不同电磁参数纤维的层间混杂分析。根据阻抗渐变原则,预制体上层应使用低介电常数纤维来提高材料和空气间的阻抗匹配,而在下层应使用高介电常数纤维提高电磁波损耗能力。按照上述混编思路,图12(c)给出了材料厚度为3 mm时,混杂不同层低介电常数SiC纤维对2.5D机织SiCf/SiC复合材料吸波性能的影响(左下角小图为混编思路示意图,黄色(上)为低介电匹配层,绿色(下)为高介电损耗层)。由图可知,随着混杂层数的增多,有效吸波带宽略有拓宽,但吸收峰逐渐向高频移动,且最大反射损耗逐渐减少。因此,当材料整体厚度较小时,不一定都是低介电常数纤维在上层作匹配层,高介电常数纤维在下层作损耗层。这与刘海韬等[34]的研究一致:阻抗渐变原则中的阻抗是指各层材料的输入阻抗,而输入阻抗不仅和材料电磁参数有关,还与材料厚度有关。为了进一步分析材料厚度对不同电磁参数纤维混杂效果的影响,图12(d)给出了材料厚度为6 mm时,混杂不同层低介电常数SiC纤维对2.5D机织SiCf/SiC复合材料吸波性能的影响。可以看出,3条反射损耗曲线都有两个明显的吸收峰。随着混杂层数的增多,两个吸收峰都略向高频移动,第一个峰对应的反射损耗逐渐减小,而第二个峰对应的反射损耗逐渐增大。整体而言,复合材料有效吸波带宽随着混杂层数的增多而逐渐增大。这说明只有当材料厚度较大时,不同电磁参数纤维混杂的方式才满足特征阻抗渐变原则,材料整体厚度越大,越有将低介电常数纤维放在上层的趋势。

    考虑到实际应用中,入射电磁波的电场极化方向是未知的,因此吸波材料的极化稳定性研究具有重要意义。图13(a)给出了不同极化角度下2.5D机织SiCf/SiC复合材料的反射损耗曲线。可以看出,当电磁波垂直入射到材料表面时,随着极化角度的增大,最大反射损耗逐渐减小,且吸收峰略向高频移动,不利于低频的高效吸收。因此本文所述的2.5D机织SiCf/SiC复合材料不具备很好的极化稳定性,这与材料的结构对称性有关。

    图  13  极化角度(phi) (a)和入射角度(θ) (b)对2.5D机织SiCf/SiC复合材料吸波性能的影响
    Figure  13.  Influence of polarization angle (phi) (a) and incidence angle (θ) (b) on the microwave absorption properties of 2.5D woven SiCf/SiC composites

    图13(b)给出了不同电磁波入射角度下复合材料的反射损耗曲线。由图可知,在0°~60°之间,最大反射损耗随着入射角度的增大而增大,这是由于材料的感应电流密度在此范围内随入射角的变大而增加,且电流密度分布不均匀,产生了极化效应[30]。而当入射角度>60°时,无有效吸波带宽,这说明本文所研究的2.5D机织SiCf/SiC复合材料具有良好的宽角度入射特性。

    (1) 实验结果显示,2.5D机织SiCf/SiC复合材料吸收峰出现在9.3 GHz处,吸收峰处的反射损耗为−17 dB,反射损耗小于−10 dB的有效吸波带宽为3 GHz。对应的仿真结果显示,吸收峰出现在10.1 GHz处,吸收峰处的反射损耗为−16 dB,有效吸波带宽也为3 GHz。整体来看,实验与仿真结果吻合良好,证明了所建模型的合理性和仿真结果的准确性。二者在吸收峰处的差异一方面是由于实验环境和理想化仿真环境不同导致的,另一方面是由于所建模型没有充分考虑材料内部孔隙和纱线挤压问题。后续研究将进一步精细化仿真模型。

    (2) 2.5D机织SiCf/SiC复合材料优异的吸波性能是材料特性和结构设计协同作用的结果。一方面,SiC纤维固有的材料属性使电磁波在材料内部以电损耗的形式被吸收;另一方面,2.5D机织SiCf/SiC复合材料内部随机分布的孔隙引起的多次反射和多次散射行为也是一个重要的吸波机制。关于材料内部孔隙分布和孔隙率对吸波性能的影响将在下一步研究中展开。

    (3) 相比经纬密和纤维层数,厚度对2.5D机织SiCf/SiC复合材料吸波性能的影响更加显著。且当厚度较大时,通过不同介电常数纤维混杂来提升材料吸波性能的效果更好。除此之外,纤维介电常数实部对材料吸收峰的位置和有效吸波带宽都有明显影响,而纤维介电损耗角正切只对吸收峰深度有影响。由于该复合材料内部结构没有明显的对称性,因此不具备良好的极化特性。电磁波入射角度的变化会引起表面电流密度极化效应,该复合材料具有良好的宽角度入射特性。这些规律为吸波复合材料的选材和结构设计提供了可靠支撑。

  • 图  1   (a) 连续纤维3D打印设备;(b) 加工示意图

    Figure  1.   (a) Continuous fiber 3D printing equipment; (b) Processing diagram

    图  2   (a) 圆形蜂窝(CH)-碳纤维(CF)试件;(b) 单向增强圆形蜂窝(SEH)-CF试件;(c) 双向增强圆形蜂窝(DEH)-CF试件

    Figure  2.   (a) Circular honeycomb (CH)-carbon fiber (CF) specimen; (b) Single enhanced circular honeycomb (SEH)-CF specimen; (c) Double enhanced circular honeycomb (DEH)-CF specimen

    图  3   (a) DEH单胞示意图;(b) SEH单胞成型轨迹;(c) CH成型轨迹

    Figure  3.   (a) Diagram of DEH monocell; (b) SEH monocell forming path; (c) CH forming path

    r2—Radius of leaf support; ; ϕ—Radian of leaf support curve; L, H—Width, height of honeycombs

    图  4   DEH成型轨迹

    Figure  4.   Forming path of DEH

    图  5   采用“支撑一体化成型”路径打印的DEH结构

    Figure  5.   DEH printed by the path of "strut integrated molding"

    图  6   实验设备及试件安装

    Figure  6.   Experimental equipment and installation of specimen

    图  7   不同蜂窝的面内压缩名义应力-应变曲线:(a) CH;(b) SEH;(c) DEH

    Figure  7.   In-plane compression stress-strain curves of different honeycombs: (a) CH; (b) SEH; (c) DEH

    图  8   各蜂窝的抗压强度

    Figure  8.   Compressive strength of honeycombs

    图  9   变形过程:(a) CH-PLA;(b) CH-CF

    Figure  9.   Deformation process: (a) CH-PLA; (b) CH-CF

    εy—Nominal strain

    图  10   变形过程:(a) SEH-PLA;(b) SEH-CF

    Figure  10.   Deformation process: (a) SEH-PLA; (b) SEH-CF

    图  11   变形过程:(a) DEH-PLA;(b) DEH-CF

    Figure  11.   Deformation process: (a) DEH-PLA; (b) DEH-CF

    图  12   蜂窝破坏模式:(a) CH-PLA;(b) CH-CF;(c) SEH-PLA;(d) SEH-CF;(e) DEH-PLA;(f) DEH-CF

    Figure  12.   Failure mode of honeycombs: (a) CH-PLA; (b) CH-CF; (c) SEH-PLA; (d) SEH-CF; (e) DEH-PLA; (f) DEH-CF

    图  13   DEH-PLA和DEH-CF侧向视图

    Figure  13.   Side view of DEH-PLA and DEH-CF

    图  14   DEH蜂窝的泊松比-名义应变曲线

    Figure  14.   Poisson's ratio-nominal strain curves of DEH

    图  15   DEH-PLA变形过程

    Figure  15.   Deformation process of DEH-PLA

    A1-A3—Left end point of DEH; B1-B3—Right end point of DEH

    图  16   DEH-CF变形过程

    Figure  16.   Deformation process of DEH-CF

    表  1   蜂窝结构设计参数

    Table  1   Geometric parameters of honeycomb

    r1/mmr2/mmϕ/(°)L/mmH/mmq/mm
    1014.1490606020
    Notes: r1—Radius of honeycomb outer circle; q—Thickness of honeycombs.
    下载: 导出CSV

    表  2   打印成型参数

    Table  2   Parameters of printing process

    tz/mmV1/(mm·min–1)T1/℃n1tH/mm
    0.4130200501.7
    Notes: tz—Thickness of layer; V1—Speed of moulding; T1—Nozzle temperature; n1—Number of layers; tH—Honey-comb wall and strut thickness.
    下载: 导出CSV

    表  3   蜂窝试验件参数

    Table  3   Parameters of honeycomb specimens

    TypeL/mmH/mmq/mmm/gρ
    CH-PLA61612018.10.24
    CH-CF61612018.80.24
    SEH-PLA61612031.10.40
    SEH-CF61612032.10.40
    DEH-PLA61612043.30.57
    DEH-CF61612044.90.57
    Notes: m—Mass of specimen; ρ—Relative density; PLA—Poly-lactic acid.
    下载: 导出CSV

    表  4   各蜂窝结构的吸能参数

    Table  4   Energy absorption parameters of each honeycomb

    TypeεdS/(J·g−1)(SCFSPLA)/SPLAM/kN(MCFMPLA)/MPLA
    CH-PLA0.76 3.62 1.45
    CH-CF0.74 5.19 43.37% 2.2051.72%
    SEH-PLA0.59 3.53 3.09
    SEH-CF0.61 5.76 63.17% 5.0061.81%
    DEH-PLA0.42 5.31 8.95
    DEH-CF0.5913.89161.58%17.5596.09%
    Notes: εd—Densification strain; S—Specific energy absorption; M—Mean crushing force.
    下载: 导出CSV
  • [1]

    WANG Z G. Recent advances in novel metallic honeycomb structure[J]. Composites Part B: Engineering,2019,166:731-741. DOI: 10.1016/j.compositesb.2019.02.011

    [2]

    FENG G Z, LI S, XIAO L J, et al. Energy absorption performance of honeycombs with curved cell walls under quasi-static compression[J]. International Journal of Mechanical Sciences, 2021, 210: 106746.

    [3]

    QI C, JIANG F, YANG S. Advanced honeycomb designs for improving mechanical properties: A review[J]. Composites Part B: Engineering,2021,227:109393. DOI: 10.1016/j.compositesb.2021.109393

    [4]

    NOVAK N, HOKAMOTO K, VESENJAK M, et al. Mechanical behaviour of auxetic cellular structures built from inverted tetrapods at high strain rates[J]. International Journal of Impact Engineering,2018,122:83-90. DOI: 10.1016/j.ijimpeng.2018.08.001

    [5] 沈振峰, 张新春, 白江畔, 等. 负泊松比内凹环形蜂窝结构的冲击响应特性研究[J]. 振动与冲击, 2020, 39(18):89-95, 117.

    SHEN Zhenfeng, ZHANG Xinchun, BAI Jiangpan, et al. Dynamic response characteristics of re-entrant circular honeycombs with negative Poisson's ratio[J]. Journal of Vibration and Shock,2020,39(18):89-95, 117(in Chinese).

    [6] 周星驰, 唐振刚, 周徐斌, 等. CFRP圆形胞元蜂窝芯层面外剪切模量[J]. 复合材料学报, 2018, 35(10):2777-2785. DOI: 10.13801/j.cnki.fhclxb.20180202.001

    ZHOU Xingchi, TANG Zhengang, ZHOU Xubin, et al. External shear modulus of CFRP circular cell honeycomb[J]. Acta Materiae Compositae Sinica,2018,35(10):2777-2785(in Chinese). DOI: 10.13801/j.cnki.fhclxb.20180202.001

    [7] 章娅菲, 闵世威, 王海涛, 等. TPU蜂窝结构及吸能特性研究[J]. 塑性工程学报, 2023, 30(3):113-122. DOI: 10.3969/j.issn.1007-2012.2023.03.016

    ZHANG Yafei, MIN Shiwei, WANG Haitao, et al. Study on compression and energy absorption characteristics of TPU honeycomb structure[J]. Journal of Plasticity Engineering,2023,30(3):113-122(in Chinese). DOI: 10.3969/j.issn.1007-2012.2023.03.016

    [8] 孙靖, 王旭琴, 柳玉文, 等. 熔融沉积成型连续碳纤维增强尼龙蜂窝芯材的压缩特性[J]. 上海航天, 2020, 37(3):32-37.

    SUN Jing, WANG Xuqin, LIU Yuwen, et al. Compressive behavior of continuous carbon fiber reinforced nylon honeycomb cores prepared by fused deposition modeling[J]. Aerospace Shanghai,2020,37(3):32-37(in Chinese).

    [9] 杨来侠, 刘波, 刘腾飞, 等. 3D打印连续纤维增强聚碳酸酯复合材料预浸丝制备与性能[J]. 复合材料学报, 2023, 40(10): 5654-5665.

    YANG Laixia, LIU Bo, LIU Tengfei, et al. Preparation and properties of 3D printing continuous fiber reinforced polycarbonate composite prepreg filaments[J]. Acta Materiae Compositae Sinica, 2023, 40(10): 5654-5665(in Chinese).

    [10]

    QUAN C, HAN B, HOU Z H, et al. 3D printed continuous fiber reinforced composite auxetic honeycomb structures[J]. Composites Part B: Engineering, 2020, 187: 107858.

    [11]

    CHENG Y Y, LI J J, QIAN X P, et al. 3D printed recoverable honeycomb composites reinforced by continuous carbon fibers[J]. Composite Structures,2021,268:113974. DOI: 10.1016/j.compstruct.2021.113974

    [12]

    ZENG C, LIU L W, BIAN W, et al. Compression behavior and energy absorption of 3D printed continuous fiber reinforced composite honeycomb structures with shape memory effects[J]. Additive Manufacturing, 2021, 38: 101842.

    [13] 张亚男. 连续芳纶纤维增强TPU复合材料蜂窝结构设计与性能研究[D]. 济南: 山东大学, 2019.

    ZHANG Yanan. Study on the design and properties of honeycomb structure of continuous PPTA fiber reinforced TPU composites[D]. Jinan: Shandong University, 2019(in Chinese).

    [14] 胡锦顺, 林永水, 陈威, 等. 改进星形蜂窝结构面内动力学响应及能量吸收特性研究[J]. 振动与冲击, 2022, 41(23):119-128. DOI: 10.13465/j.cnki.jvs.2022.23.015

    HU Jinshun, LIN Yongshui, CHEN Wei, et al. In-plane dynamic response and energy absorption characteristics of improved star-shape honeycomb structure[J]. Journal of Vibration and Shock,2022,41(23):119-128(in Chinese). DOI: 10.13465/j.cnki.jvs.2022.23.015

    [15] 杨玉平, 张中伟, 李玮洁, 等. 碳/碳蜂窝制备工艺及压缩与剪切行为[J]. 复合材料学报, 2023, 40(12): 6639-6648.

    YANG Yuping, ZHANG Zhongwei, LI Weijie, et al. Preparation process, compression and shear behavior of carbon/carbon honeycomb[J]. Acta Materiae Compositae Sinica, 2023, 40(12): 6639-6648(in Chinese).

    [16]

    YANG X F, XI X L, PAN Q F, et al. In-plane dynamic crushing of a novel circular-celled honeycomb nested with petal-shaped mesostructure.[J]. Composite Structures,2019,226:111219. DOI: 10.1016/j.compstruct.2019.111219

    [17]

    LIN T C, CHEN T J, HUANG J S. In-plane elastic constants and strengths of circular cell honeycombs[J]. Composites Science and Technology,2012,72(12):1380-1386. DOI: 10.1016/j.compscitech.2012.05.009

    [18]

    ORUGANTI R K, GHOSH A K. FEM analysis of transverse creep in honeycomb structures[J]. Acta Materialia,2008,56(4):726-735. DOI: 10.1016/j.actamat.2007.10.019

    [19] 张帆, 魏培祥, 赵圆圆, 等. 基于体素的机械臂连续碳纤维3D打印路径规划[J]. 计算机集成制造系统, 2023, 29(5):1517-1527.

    ZHANG Fan, WEI Peixiang, ZHAO Yuanyuan, et al. Voxel-based path planning method for continuous carbon fiber 3D printing with robot arm[J]. Computer Integrated Manufacturing Systems,2023,29(5):1517-1527(in Chinese).

    [20] 韩宁达. 连续纤维复合材料构件三维打印工艺及机电系统设计[D]. 南京: 南京师范大学, 2019.

    HAN Ningda. 3D printing technology and electromechanical system design of continuous carbon fiber composite components[D]. Nanjing: Nanjing Normal University, 2019(in Chinese).

    [21] 秦若森, 孙守政, 韩振宇, 等. 3D打印连续纤维增强热塑性复合材料成型质量的研究进展[J]. 材料导报, 2022, 36(17):200-208.

    QIN Ruosen, SUN Shouzheng, HAN Zhengyu, et al. 3D printing for continuous fiber-reinforced thermoplastic composites: A review on molding quality[J]. Materials Reports,2022,36(17):200-208(in Chinese).

    [22] 田小永, 张亚园, 刘腾飞, 等. 连续碳纤维增强尼龙复合材料预浸丝制备与3D打印性能研究[J]. 航空制造技术, 2021, 64(15):24-33.

    TIAN Xiaoyong, ZHANG Yayuan, LIU Tengfei, et al. Prepreg preparation and 3D printing of continuous carbon fiber reinforced nylon composite[J]. Aeronautical Manufacturing Technology,2021,64(15):24-33(in Chinese).

    [23]

    MATSUZAKI R, UEDA M, NAMIKI M, et al. Three-dimensional printing of continuous-fiber composites by in-nozzle impregnation[J]. Scientific Reports,2016,6:23058. DOI: 10.1038/srep23058

    [24]

    WU H X, ZHANG X C, LIU Y. In-plane crushing behavior of density graded cross-circular honeycombs with zero Poisson's ratio[J]. Thin-Walled Structures,2020,151:106767. DOI: 10.1016/j.tws.2020.106767

    [25] 中国国家标准化管理委员会. 高聚物多孔材料 压缩应力应变特性的测试: GB/T 18942—2003[S]. 北京: 中国标准出版社, 2003.

    Standardization Administration of the People's Republic of China. Flexible cellular polymeric materials—Determination of stress-strain characteristics in compression: GB/T 18942—2003[S]. Beijing: China Standards Press, 2003.

    [26] 魏路路, 余强, 赵轩, 等. 内凹-反手性蜂窝结构的面内动态压溃性能研究[J]. 振动与冲击, 2021, 40(4):261-269. DOI: 10.13465/j.cnki.jvs.2021.04.036

    WEI Lulu, YU Qiang, ZHAO Xuan, et al. In-plane dynamic crushing characteristics of re-entrant anti-trichiral honeycomb[J]. Journal of Vibration and Shock,2021,40(4):261-269(in Chinese). DOI: 10.13465/j.cnki.jvs.2021.04.036

  • 期刊类型引用(1)

    1. 邱海鹏,王晓猛,陈明伟,王岭,徐俊杰,赵马娟,张典堂,宋雨晴. 层合结构吸波陶瓷基复合材料设计及性能研究. 航空制造技术. 2024(Z2): 34-42 . 百度学术

    其他类型引用(0)

  • 目的 

    圆形蜂窝结构受面内压缩时,在抗压强度和吸能性能上相较其它基础构型蜂窝有更好的性能。当前基于复杂构型蜂窝的性能研究多使用熔融式三维打印技术制造试验件,但由于纯树脂的力学性能较低,导致试件性能不佳。本文使用连续纤维三维打印技术制造出蜂窝试验件并开展面内压缩试验,探究圆形蜂窝在新材料体系和新构型下的抗压和吸能性能。

    方法 

    在试验件制造方面,使用连续纤维三维打印制造出三种构型的蜂窝试验件,同时设置聚乳酸(PLA)三维打印试验件对照组。三维打印设备基于预浸渍工艺自行研制,使用了直径为1mm的纤维体积含量为8.3vol%的碳纤维增强聚乳酸基(CF/PLA)丝材。以圆形蜂窝(CH)构型为基础,在横向和竖向上增设树叶形支撑,设计了单向增强圆形蜂窝(SEH)和双向增强圆形蜂窝(DEH)两种改进型蜂窝。讨论并规划了以上三种构型蜂窝碳纤维增强试验件的纤维成型路线。在试验方面,使用MTS Model E45型试验机对蜂窝结构进行准静态压缩试验,获取了试件面内压缩时的载荷位移曲线,并对试验数据进行处理,得到比吸能和平均压缩力等参数。并借助光学显微镜对试件的断口形貌进行了观察。

    结果 

    对试验件的面内准静态压缩试验结果表明,同种构型蜂窝在采用碳纤维增强材料后,重量增加约4%的同时,各项性能提升明显:在抗压强度上,CH、SEH和DEH结构在采用CF增强后,相比PLA对照组分别提高83%、43%和40%,比吸能分别提高43%、63%和162%,平均压缩力分别提高52%、62%和96%。在圆形蜂窝基础上提出的改进型蜂窝DEH拥有最高的抗压能力和吸能特性,DEH-CF的抗压强度和比吸能分别为12.3Mpa和13.9J/g。研究发现,连续碳纤维三维打印结构内部的纤维路径规划会影响结构的刚度和变形行为,采用“支撑一体化成型”路径的DEH-CF在压缩时,其结构动态泊松比保持在PLA对照组35%以下。

    结论 

    连续碳纤维三维打印技术制造的蜂窝试件相较对照组在抗压和吸能性能上提升明显,结构中的碳纤维束阻止了树脂基体中裂纹的扩展,防止结构的侧壁和支撑在受压时过早断裂,有效提升了结构的抗压缩和吸能性能。基于圆形蜂窝构型设计的SEH和DEH结构在准静态压缩下,应力应变曲线表现出了多胞材料的典型应力响应特征,即出现弹性阶段、平台阶段和密实化阶段。在连续纤维三维打印成型中,结构内部的纤维路径设计可影响结构的刚度和变形模式。

图(16)  /  表(4)
计量
  • 文章访问数:  929
  • HTML全文浏览量:  507
  • PDF下载量:  90
  • 被引次数: 1
出版历程
  • 收稿日期:  2023-05-25
  • 修回日期:  2023-06-25
  • 录用日期:  2023-07-05
  • 网络出版日期:  2023-07-23
  • 刊出日期:  2024-03-31

目录

/

返回文章
返回