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

复合材料固化工艺的直热模具温度场均匀性分析

陆阳, 齐俊伟, 肖军, 石甲琪

陆阳, 齐俊伟, 肖军, 等. 复合材料固化工艺的直热模具温度场均匀性分析[J]. 复合材料学报, 2021, 38(9): 2954-2963. DOI: 10.13801/j.cnki.fhclxb.20201222.003
引用本文: 陆阳, 齐俊伟, 肖军, 等. 复合材料固化工艺的直热模具温度场均匀性分析[J]. 复合材料学报, 2021, 38(9): 2954-2963. DOI: 10.13801/j.cnki.fhclxb.20201222.003
LU Yang, QI Junwei, XIAO Jun, et al. Analysis on the uniformity of temperature field of direct heating mold for composite material curing process[J]. Acta Materiae Compositae Sinica, 2021, 38(9): 2954-2963. DOI: 10.13801/j.cnki.fhclxb.20201222.003
Citation: LU Yang, QI Junwei, XIAO Jun, et al. Analysis on the uniformity of temperature field of direct heating mold for composite material curing process[J]. Acta Materiae Compositae Sinica, 2021, 38(9): 2954-2963. DOI: 10.13801/j.cnki.fhclxb.20201222.003

复合材料固化工艺的直热模具温度场均匀性分析

基金项目: “高档数控机床与基础制造装备”科技重大专项(2017ZX04009001)
详细信息
    通讯作者:

    齐俊伟,硕士,高级工程师,硕士生导师,研究方向为先进复合材料自动化成型技术  E-mail:qijunwei@nuaa.edu.cn

  • 中图分类号: TB332

Analysis on the uniformity of temperature field of direct heating mold for composite material curing process

  • 摘要: 针对复合材料固化成型工艺的直热模具温度场均匀性进行了研究。建立模具温度场和复合材料固化反应温度场的耦合传热学模型,并对该模型进行有限元建模仿真分析。针对影响模具表面温度均匀性的主要因素,即电加热管的间距和功率,设计正交试验优化,优化后模具表面最大温差为3.5℃,达到行业标准。此外,对影响温度场均匀性的其他因素,即加热管与模具的接触热阻、复合材料层合板厚度进行了探讨,接触热阻的存在使得模具表面最大温差达到7.24℃,模具加热到指定温度多用时800 s,降低了效率。研究层合板对模具温度均匀性的影响时发现未加入复合材料时模具表面最大温差为4.44℃,加入层合板耦合后最大温差为3.5℃;厚度为毫米级时,层合板对直热模具表面温度均匀性影响不大。
    Abstract: The uniformity of the temperature field of the direct heating mold of the composite material curing molding process was studied. A coupled heat transfer model of mold temperature field and composite material curing reaction temperature field was established, and finite element modeling and simulation analysis was performed on the model. In view of the main factors that affect the uniformity of the mold surface temperature, which are the distance and power of the electric heating tube, the orthogonal experiment was designed to optimize, and the maximum temperature difference of the mold surface after optimization was 3.5℃, which reached the industry standard. In addition, other factors that affect the uniformity of the temperature field, the thermal resistance of the contact between the heating tube and the mold, and the thickness of the composite laminate were discussed. The existence of contact thermal resistance makes the maximum temperature difference of the mold surface reach 7.24℃, and it takes 800 s to heat the mold to the specified temperature, which reduces the efficiency. When studying the effect of laminates on mold temperature uniformity, it is found that the maximum temperature difference on the mold surface without composite materials is 4.44℃; the maximum temperature difference after coupling with laminates is 3.5℃, and the thickness of the laminate is millimeter-level. The uniformity has little effect.
  • 聚对苯二甲酸乙二醇酯(PET)是全球用量最大的高分子材料之一,其纤维制品俗称涤纶,是全球第一大化纤品种[1]。但PET属于易燃材料,且其燃烧时还伴有严重的熔滴现象,极易导致火灾蔓延和二次伤害,这为PET的应用带来了极大的安全隐患,也使其在军事、工业等诸多领域的应用受到限制[2]。目前PET阻燃改性的主流方法是引入含磷阻燃剂,但含磷聚酯的阻燃性主要是通过熔体滴落带走燃烧区域的热量和火焰来实现的[3-4],未能解决PET的熔滴问题。因此,同时赋予PET阻燃和抗熔滴特性是目前PET阻燃改性面临的一大难点。

    相关研究表明,提高聚合物高温下的熔体黏度和炭化能力是实现阻燃和抗滴落的关键[5-6]。但由于PET分子链的线性结构,其在高温下具有较低的熔融黏度,这一特性在赋予其可纺性的同时,也导致其在燃烧过程中表现出严重的熔滴行为,并难以形成连续的炭层[7]。在PET中引入可交联结构单元是提高熔体黏度的一个有效途径[6, 8]。例如,Wu等[9]合成了一种芳香族席夫碱5-(亚苄基-氨基)-间苯二甲酸二甲酯,并将其用作PET的共聚单体;研究结果显示,芳族席夫碱可以在PET的熔融温度和分解温度之间形成稳定的交联网络,并在燃烧中进一步转变为致密的炭层,因而使得PET共聚酯显示出优异的自熄和抗滴落行为。但尽管上述共聚酯表现出较为理想的阻燃性和抗滴落性能,问题仍然存在。除了工艺复杂、成本高、工业化困难之外,共聚法往往会破坏分子链的规整性,从而损害PET的力学性能和可纺性。

    近年来,碳基阻燃剂作为一种新型阻燃剂表现出巨大的发展潜力。研究表明,碳纳米管[10]、富勒烯[11]、纳米炭黑[12]、石墨烯[13]、碳微球[14]等碳基阻燃剂在改善聚合物的成炭质量、降低热释放速率、提高热稳定性等方面表现突出。与传统阻燃剂相比,往往少量碳基阻燃剂的引入即可显著提高聚合物的阻燃性,除此之外,引入碳基阻燃剂还能不同程度地改善聚合物的力学、热学以及电学等性能[15]。其中,碳纳米球(Carbon nanospheres,CNSs)具有无卤环保、粒径小、热稳定性高等优势,能满足PET高温加工和熔融纺丝的理论要求,但迄今为止,将碳基阻燃剂应用于PET的研究较少。

    基于此,本文以CNSs为基体,在其表面接枝芳香席夫碱4-苯亚甲基氨基苯酚(4-phenyl-methyleneamino-phenol,BA)制备了一种新型碳纳米球基复合阻燃剂(CNSs-BA),旨在将CNSs和BA的优势有机结合,使阻燃剂同时具备“智能自交联”特性,从而在燃烧时在PET基体中形成三维交联网络,进而改善熔滴。重点研究了CNSs-BA/PET复合材料的阻燃性及其阻燃机制。

    蒸馏水,使用XY-ZL-20型蒸馏水器自制。碳纳米球(CNSs),纯度99.99%,宁波金雷纳米材料科技有限公司;30%过氧化氢,优级纯,上海沃凯生物技术有限公司;硫酸,分析纯,华东医药股份有限公司;二氯甲烷(DCM),分析纯,华东医药股份有限公司;4-苯亚甲基氨基苯酚(BA),纯度98%,东京化成工业株式会社;4-二甲氨基吡啶(DMAP),纯度99%,阿拉丁试剂(上海)有限公司;无水乙醇,分析纯,嘉兴市甬宏化工有限公司;PET切片,半消光型SD500,中国石化仪征化纤股份有限公司。

    首先采用酸化法[16]制备羧基化碳纳米球(CNSs-COOH);在DCM中加入适量氯化钙,常温振荡12 h以去除其中的水分,使用时用吸管吸取上层液体用。在三口烧瓶中加入一定量的无水DCM (50~100 mL)作为反应溶剂。搅拌状态下加入3 g CNSs-COOH,再加入1 g BA和DMAP (作为碱性催化剂,DMAP的用量为BA用量的1/10),升温至30℃,之后在搅拌状态下反应2 h,反应完成后抽滤除去液体。将反应所得固体依次用乙醇和蒸馏水洗涤,100℃下干燥6 h后研磨均匀,即得CNSs-BA阻燃剂,CNSs-BA的制备示意图见图1

    图  1  碳纳米球(CNSs)-芳香席夫碱4-苯亚甲基氨基苯酚(BA)制备示意图
    Figure  1.  Schematic diagram of the preparation of carbon nanospheres (CNSs)-Schiff bases 4-phenyl-methyleneamino-phenol (BA)

    首先将纯PET切片和制备的CNSs-BA阻燃剂分别在120℃下真空干燥12 h后冷却至室温备用,然后分别按CNSs-BA占PET的质量分数为0.5wt%、1.0wt%、2.0wt%、3.0wt% 的比例将CNSs-BA与PET切片混合均匀后喂入BP-8188型转矩流变仪(东莞市宝品精密仪器有限公司)中,转矩流变仪各区的温度依次为150、255、273、275℃,转速为35~60 r/min,熔体依次经过熔融共混、挤出、切粒工序得到CNSs-BA/PET复合材料。

    CNSs-BA阻燃剂的测试与表征:用JSM-6510LA型场发射扫描电镜(SEM,日本电子株式会社)和EM420型透射电子显微镜(赛默飞世尔科技)观察阻燃剂的微观形貌,加速电压3 kV 和20 kV。用Perkin Elmer Frontier型傅里叶变换红外光谱仪(FTIR),溴化钾压片法分析测定阻燃剂表面性质和化学结构,光谱记录范围4000~400 cm−1。用Perkin Elmer TG4000型热重分析仪(TG),在N2气氛下测试阻燃剂的热稳定性,气体流速20 mL/min,程序设定为:30℃恒温1 min后以20℃/min的升温速率升温至800℃。

    CNSs-BA/PET复合材料的测试与表征:用TM606 数显氧指数测试仪(青岛睿新杰仪器有限公司),按照GB/T 2406.2—2009[17]测试PET及其阻燃复合材料的极限氧指数(LOI),样条尺寸为120 mm×6.5 mm×3 mm。用CZF-5水平垂直燃烧仪(沧州冀路试验仪器有限公司),按照 GB/T 2408—2008[18]判定PET及其阻燃复合材料的UL-94垂直燃烧等级,样品尺寸为130 mm×13 mm×3 mm。用C-1087型锥形量热仪(英国FTT),按照ISO 5660-1: 2015[19]测试PET及其阻燃复合材料的燃烧热释放(HRR)等参数,样品尺寸:100 mm×100 mm×3 mm,辐射照度50 kW/m2。用Perkin Elmer TG4000型热重分析仪(TG),在N2气氛下测试PET及其阻燃复合材料(阻燃剂含量2wt%)的热稳定性,气体流速20 mL/min,程序设定为:30℃恒温1 min后以20℃/min的升温速率升温至800℃。用Netzsch STA449F5型同步热分析仪(TG-DSC)研究PET及其阻燃复合材料(阻燃剂含量2.0wt%)的交联行为,氩气做保护气,空气气氛,气体流速20 mL/min,升温速率为10℃/min,测试温度范围为30~800℃。用气相Agilent 6980N色谱仪,Agilent 5975质谱仪,采用HP-5MS色谱柱对PET及其阻燃复合材料(阻燃剂含量2.0wt%)进行裂解-气相色谱-质谱联用(Py-GC-MS)测试。裂解条件:裂解温度750℃,时间20 s,升温速率200℃/s。色谱条件:柱温在50℃保持5 min,然后以10℃/min升温至260℃,在260℃保持10 min;进样温度220℃,传输温度280℃,He做载气,流量1.0 mL/min;裂解产物经色谱柱分离后进入质谱仪,电子能量为70 eV。

    残炭的测试与表征:采用SEM观察PET及其阻燃复合材料(阻燃剂含量2.0wt%)燃烧后残炭的形貌,并用其配套的EDS能谱设备对残炭进行元素分析。采用TG在N2气氛下测试残炭的热稳定性,气体流速20 mL/min,程序设定为:30℃恒温1 min后以20℃/min的升温速率升温至800℃。

    图2为原CNSs (图2(a))和CNSs-BA (图2(b))的SEM、TEM和EDS能谱图。可知:CNSs和CNSs-BA二者均呈规则的球形颗粒状。不同的是,原始CNSs表面光滑,平均粒径约45 nm。而经BA接枝后,CNSs-BA的表面变得粗糙,平均粒径增大到50 nm左右,由EDS谱图可知:纯CNSs中的主要成分为C元素,与CNSs相比,CNSs-BA的表面增加了N元素,源自其表面接枝的BA中的氨基。

    图  2  CNSs (a)和CNSs-BA (b)的SEM图像和EDS能谱图
    Figure  2.  SEM and EDS images of CNSs (a) and CNSs-BA (b)

    图3是CNSs和CNSs-BA的红外图谱。对比CNSs和CNSs-BA的FTIR曲线可知,在CNSs-BA的红外曲线中,3330 cm−13380 cm−1处对应N—H的伸缩振动峰,2925 cm−12850 cm−1处为亚甲基的伸缩振动峰,1450 cm−1处的特征峰是芳环骨架的伸缩振动峰,1269 cm−1处的特征峰是酯基C(O)—O的伸缩振动峰,1045 cm−1处的特征峰是芳环上1, 4位取代的振动峰,以上特征峰源自CMSs表面接枝的BA。

    图  3  CNSs和CNSs-BA的红外图谱
    Figure  3.  Infrared spectra of CNSs和CNSs-BA

    图4为CNSs和CNSs-BA的TG曲线。可知:纯CNSs的初始分解温度(Tonset,定义为热失重5wt%时的温度)大于800℃,经BA接枝后Tonset降低到483.1℃,该温度远高于PET的加工温度和热分解温度。纯CNSs和CNSs-BA的最高热分解速率的温度(Tmax)分别为158.4℃和177.4℃,说明CNSs经BA接枝后热分解速率减慢。CNSs在30~800℃之间表现出3个较为明显的失重阶段,338.7℃之前对应CNSs中的少量结晶水和无定形碳的分解,338.7~544.1℃之间对应CNSs主体的热分解,544.1℃之后对应热分解产物的再分解。而CNSs-BA可划分为4个失重阶段:前两个阶段分别对应结晶水、无定形碳、小分子产物的分解以及CNSs-BA主体的分解,值得注意的是,544.1℃之后,CNSs-BA出现一个较为明显的失重峰,而CNSs的DTG曲线上并无该峰,该失重峰的出现证明CNSs表面接枝的BA在第二阶段(PET的熔融温度和分解温度之间)形成了一个较为稳定的交联网络结构,这将十分有助于燃烧时保护炭层的形成。

    图  4  CNSs和CNSs-BA的TG (a)和DTG (b)曲线
    Figure  4.  TG (a) and DTG (b) curves of CNSs and CNSs-BA
    Tonset—Initial weightlessness temperature; Tmax—Maximum weightlessness temperature

    表1是CNSs-BA/PET复合材料的LOI和UL-94垂直燃烧测试结果。可知,与CNSs/PET相比,CNSs-BA/PET复合材料的LOI进一步提高,二者LOI规律变化一致,即随着阻燃剂含量的增大,LOI先提高后降低,当CNSs-BA含量为2.0wt%时,CNSs-BA/PET的LOI指数达到最大值28.1%,此时与纯PET相比,CNSs-BA/PET的LOI提高了33.8%。UL-94垂直燃烧测试结果表明,CNSs-BA/PET复合材料的抗熔滴性能较CNSs/PET也有明显提高,两次施加火焰后的余焰时间明显缩短,当CNSs-BA的添加量超过2.0wt%时,CNSs-BA/PET复合材料的阻燃等级可达到V-0级。

    表  1  复合材料的极限氧指数(LOI)和UL-94垂直燃烧测试结果
    Table  1.  Limiting oxygen index (LOI) and UL-94 vertical burning test results of composites
    Sample Flame retardant content/wt% LOI/% UL-94 vertical combustion test results
    t1/s t2/s t3/s Ignite cotton? Rate
    PET 21.0 Burn out Yes NR
    CNSs/PET 0.5 23.2 2.6 2.5 0 Yes V-2
    1.0 25.0 2.4 2.4 0 Yes V-2
    2.0 26.2 2.4 2.8 0 Yes V-2
    3.0 24.6 3.1 2.2 0 Yes V-2
    CNSs-BA/PET 0.5 24.0 1.5 2.3 0 Yes V-2
    1.0 26.9 1.2 2.1 0 Yes V-2
    2.0 28.1 0.5 2.2 0 No V-0
    3.0 27.5 0.6 1.9 0 No V-0
    Notes: PET—Polyethylene terephthalate; t1—Afterglow time after the first application of flame; t2—Afterglow time after the second application of flame; t2—Afterglow time; NR—No rate.
    下载: 导出CSV 
    | 显示表格

    锥形量热仪测试结果见图5表2。热释放速率(HRR)是表征材料火灾危险性的主要依据。结合图5表2可知,纯PET被点燃后热释放速率急剧增大,其峰值热释放速率(pk-HRR)为810.45 kW/m2,总热释放(THR)为150.27 MJ/m2。与之相比,CNSs-BA/PET的THR与之接近,但HRR曲线却明显变平缓。值得注意的是,CNSs-BA/PET的HRR曲线表现出两个明显的热释放阶段,即在热释放速率达到峰值之后又出现了一个较为平缓的放热平台(当CNSs-BA含量为0.5wt%时表现为放热峰),这意味着CNSs-BA/PET在燃烧过程中的热释放受到抑制,这是由燃烧时炭层的形成或可燃气体减少导致的[20]。此外,表2表明,与CNSs/PET相比,CNSs-BA/PET复合材料的pk-HRR进一步降低。当CNSs-BA 含量为2.0wt%时,CNSs-BA/PET的pk-HRR最小,为435 kW/m2,该值与相同阻燃剂含量的CNSs/PET相比降低了7.6%,较纯PET降低了46.3%,说明CNSs经BA接枝后对PET的燃烧抑制作用进一步增强,阻燃效果进一步提高。

    图  5  CNSs-BA/PET的热释放速率(HRR) (a)和总热释放(THR) (b)曲线
    Figure  5.  Curves of heat release rate (HRR) (a) and total heat release (THR) (b) of CNSs-BA/PET
    表  2  复合材料的锥形量热仪测试数据
    Table  2.  Data of cone calorimeter test of composites
    Sample FR content/wt% TTI/s Time to pk-HRR/s pk-HRR/(kW·m−2) THR/(MJ·m−2)
    PET 0 47 104 810.45 150.27
    CNSs/PET 0.5 44 34 528.96 151.64
    1 34 34 503.44 148.03
    2 40 41 470.72 146.06
    3 30 29 501.49 143.19
    CNSs-BA/PET 0.5 35 39 485.54 146.54
    1 31 55 469.98 156.04
    2 34 39 435.00 146.54
    3 30 39 466.05 156.69
    Notes: TTI—Time to ignition; pk-HRR—Peak heat release rate; FR—Flame retardant.
    下载: 导出CSV 
    | 显示表格

    为了研究阻燃剂的引入对PET的热降解行为的影响,对PET、CNSs/PET和CNSs-BA/PET在氮气气氛下的TG-DTG曲线作了对比分析,并计算了CNSs/PET和CNSs-BA/PET在500℃时残炭量的理论值,如图6表3所示。由图6可知:在氮气气氛下,PET、CNSs/PET和CNSs-BA/PET三者的TG曲线和DTG曲线基本重合,说明加入少量(2.0wt%)的CNSs和CNSs-BA均不会对PET的无氧降解行为造成明显影响。由表4可知,PET、CNSs/PET和CNSs-BA/PET三者的TonsetTmax均较为接近,但三者在高温(500℃)下的残余质量有所不同。经计算发现CNSs/PET和CNSs-BA/PET二者在高温下残炭量的实际值(CR500℃,exp)均大于理论值(CR500℃,cal),这说明阻燃剂CNSs和CNSs-BA对PET均有促进成炭作用。其中,CNSs-BA/PET在500℃下残炭量的实际值与理论值的差值(∆CR500℃)大于CNSs/PET,这说明CNSs经BA接枝后对PET的促进成炭作用加强。聚合物在高温下形成的残炭越多,燃烧时发生热分解的部分就越少[21],这便是阻燃复合材料热释放速率降低的主要原因之一。

    图  6  CNSs、CNSs-BA以及PET、CNSs/PET、CNSs-BA/PET的TG (a)和DTG (b)曲线
    Figure  6.  TG (a) and DTG (b) curves of CNSs, CNSs-BA, PET, CNSs/PET and CNSs-BA/PET
    表  3  CNSs、CNSs-BA以及PET、CNSs/PET、CNSs-BA/PET在氮气气氛下的TG-DTG数据
    Table  3.  TG-DTG data of CNSs, CNSs-BA, PET, CNSs/PET and CNSs-BA/PET under nitrogen atmosphere
    Sample Tonset/℃ Tmax/℃ CR500℃/% CR500℃/%c
    exp.a/cal.b
    CNSs >800 96.88/—
    CNSs-BA 476.4 94.92/—
    PET 379.1 419.1 10.09/—
    CNSs/PET 380.1 421.4 13.97/11.52 2.45
    CNSs-BA/PET 382.0 420.1 15.68/11.79 3.89
    Notes: a CR500℃,exp. is the experimental value of char residue; b CR500℃,cal. is the calculated value of char residue; cCR500℃=CR500℃,exp.CR500℃,cal..
    下载: 导出CSV 
    | 显示表格
    表  4  PET、CNSs/PET和CNSs-BA/PET在空气气氛下的TG-DTG数据
    Table  4.  TG-DTG data of PET, CNSs/PET and CNSs-BA/PET under air atmosphere
    Sample Tonset/℃ Tmax-1/℃ Tmax-2/℃
    PET 397.3 433.4 585.1
    CNSs/PET 359.1 438.4 567.5
    CNSs-BA/PET 391.0 439.7 563.3
    Notes: Tmax-1—Maximum weightlessness temperature in the first stage; Tmax-2—Maximum weightlessness temperature of the second stage.
    下载: 导出CSV 
    | 显示表格

    对纯PET、CNSs/PET和CNSs-BA/PET锥形量热仪测试后的残炭做了SEM和TG分析以进一步研究阻燃机制。

    炭层的形貌结构和稳定性对于提高聚合物的阻燃性能至关重要,有效的炭层可通过阻止聚合物内部与可燃气体、氧气的接触来实现阻燃目的。图7为纯PET、CNSs/PET和CNSs-BA/PET炭层的SEM图像。可见,纯PET燃烧生成的炭层稀薄空且松散,表面存在大量气体逸出形成的气孔,显然这种形貌的炭层无法形成有效的屏障作用。与纯PET相比,CNSs/PET的炭层的致密性明显提高,气孔明显变小,意味着炭层有效性的提高。值得注意的是,与CNSs/PET相比,CNSs-BA/PET炭层的致密性和连续性得到了进一步改善,表面气孔也明显变少和变小,另外还存在大量鼓起的未破裂气泡,这种形貌的炭层在燃烧时一方面能有效地阻隔热量的传递,另一方面还能有效地阻隔PET燃烧降解生成的气态可燃物的逸出,起到隔热和隔氧的作用[22-23]。除此之外,CNSs-BA受热分解生成的CO2、氨气、氮气等难燃性气体能够稀释燃烧区域可燃气体的浓度,抑制燃烧的发展,这便是CNSs-BA/PET阻燃性提高的重要原因。

    图  7  PET (a)、CNSs/PET (b)和CNSs-BA/PET (c)残炭的SEM图像
    Figure  7.  SEM images of the char residues of PET (a), CNSs/PET (b)和CNSs-BA/PET (c)

    图8是PET、CNSs/PET和CNSs-BA/PET炭层的TG曲线。可知,纯PET炭层的Tonset较低,为215.08℃,其中100℃前失重为4.36wt%,这主要是由于纯PET的炭层结构松散、孔洞较多,容易吸收水分和储存小分子气体所致,其800℃时的残余质量为87.2wt%。在整个升温过程中,纯PET的炭层表现出3个失重阶段,第一个失重阶段发生在100℃之前,主要对应炭层中贮存的水分以及气态小分子的降解;第二个失重阶段发生在100~530℃之间,对应炭层主体部分的降解;第三个失重阶段发生在530℃之后,对应炭层热降解产物的再降解。与之相比,CNSs/PET炭层的Tonset提高到615.37℃,800℃时的残余质量提高到89.1wt%,这主要是由于CNSs/PET的炭层的致密性提高所致,其TG曲线基本保持了纯PET炭层的3个失重阶段。与PET和CNSs/PET的炭层相比,CNSs-BA/PET炭层的Tonset提高到800℃以上,意味着炭层在燃烧时能耐受更高的温度,从而更持久有效地起到凝聚相阻燃作用。值得注意的是,其炭层在热分解过程中只有一个较为明显的失重平台,并未像PET和的CNSs/PET的炭层一样经历3个失重阶段,说明阻燃剂CNSs-BA能使PET燃烧生成结构稳定的炭层,该炭层在燃烧过程中能耐受较高的火焰温度,从而对内部的基体起到持久有效的保护作用。

    图  8  残炭的TG曲线
    Figure  8.  TG curves of the char residues

    聚合物的交联直接影响其热性能、流变性、成炭性、熔滴和自熄行为,并有助于聚合物的芳香化或炭化[24]图9是PET、CNSs/PET和CNSs-BA/PET在热氧降解过程中的TG-DSC曲线,相关数据见表4。由图9可以看出,PET、CNSs/PET和CNSs-BA/PET在空气中均有两个失重阶段,说明PET及其复合材料发生的是两步降解反应[25]。第一个失重阶段是PET的主要失重阶段,发生在360~470℃之间。第二个失重阶段发生在470~590℃之间,该阶段对应第一个降解阶段生成的降解产物的进一步降解。值得注意的是,纯PET的Tonset为397.3℃,而PET的燃点通常在420℃左右,这说明PET在燃烧之前,首先会发生一定程度的降解并生成一些可燃性的气体或挥发性产物,以此来维持燃烧的进行。与纯PET相比,在第一个失重阶段,CNSs/PET和CNSs-BA/PET的热失重曲线稍向低温方向移动,但二者在第一个降解阶段结束时的剩余质量却大于PET,且该阶段的最大失重率所对应温度(Tmax-1)大于PET,说明阻燃剂的存在使PET的主体降解提前,但在该阶段却重组生成了热稳定性较高的物质。DSC曲线表明,CNSs-BA/PET在熔融峰和分解峰之间出现了明显的放热峰,该峰是PET的交联峰[26],而在PET和CNSs/PET的DSC曲线上交联峰却不明显,这说明CNSs经BA接枝改性后促进了PET的交联,这是由于阻燃剂表面接枝的芳香族席夫碱(BA)可以在PET的熔融温度和分解温度之间形成稳定的交联网络。

    图  9  PET、CNSs/PET和CNSs-BA/PET的TG-DSC曲线
    Figure  9.  TG-DSC curves of PET, CNSs/PET and CNSs-BA/PET

    裂解-气相色谱-质谱联用(Py-GC-MS)是目前研究聚合物高温裂解产物的常用方法[27]。为了研究阻燃剂的引入对PET的热裂解行为及其高温裂解产物的影响,对PET、CNSs/PET和CNSs-BA/PET做了Py-GC-MS分析,三者的高温裂解产物对比见表5。可知,与纯PET的裂解产物相比,CNSs/PET和CNSs-BA/PET的裂解产物中都包含更多的杂环、稠环、共轭芳环类化合物,这些裂解产物具有较高的热稳定性,是难燃性的保护炭层形成的物质基础[28]。而CNSs-BA/PET的裂解产物中出现了诸如二甲基胺、偶氮苯等含氮产物,这是由于阻燃剂CNSs-BA表面接枝的苯亚甲基氨基苯酚所致。另外,与PET和CNSs/PET相比,CNSs-BA/PET的裂解产物中菲、萘、苊、芴等稠环芳烃类以及联苯类产物明显增多,佐证了CNSs-BA促进了PET降解过程中的交联,该交联一方面通过增大熔体黏度改善了熔滴现象,另一方面提高了炭层的致密性和热稳定性,这就是CNSs-BA/PET阻燃性和抗熔滴性提高的主要原因。

    表  5  PET、CNSs/PET和CNSs-BA/PET裂解产物
    Table  5.  Pyrolysis products of PET, CNSs/PET and CNSs-BA/PET
    Pyrolysis products which found only in PET Tetrahydropyran; 2,2-dimethylpropanal; 4,8,12-trimethyl-tridecanoic acid methyl ester; 2,2-dimethoxybutane; 2-methyl-1,5-hexadien-3-yne; 1,6-heptadiyne; p-xylene; Decane; Methyl benzoate; Dodecylethyl ketone; 1-(3-methylphenyl)benzyl(2-methyl-1-methylenepropylidene); 4-methylphenyl-1-pentyn-3-ol phenol; Dimethyl 1,3-benzenedicarboxylate; Vinylmethyl terephthalate; Diphenylacetylene; Biphenyl-4-ylacetophenone; 1-(5,5-dimethyl-1,3-dioxocyclohexan-2-ylidene)-2-(N-ethylbenzothiazol-2-ylidene)-ethanes; Phthalic acid 4-formylphenyl ester; o-tertiaryl tricyclic [8.2.2.2(4,7)]hexadeca-2,4,6,8,10,12,13,15-octene; 4-(diethylaminomethyl)-2,5-dimethylphenol
    Pyrolysis products which found only in CNSs/PET Phenol; 1,2-dihydro-indene; 1-(4-methylphenyl)-ethanone; Stilbene; 1H-cyclopropyl[l]phenanthrene; Dihydro-p-terphenyl; 1-naphthol; Fluorene-9-methanol; 2-ethyl-1,1'-biphenyl; 1,1-diphenylethene; 4-(2-benzoyl-5-phenyl-3-thienyl)-1,2-dihydrophenanthrene; 2-phenylnaphthalenyl benzoate; 1,1-dihydro-2-phenylnaphthalenyl benzoate; 3-chlorobenzylnonyl; 1-(2,5-dimethylphenyl)ethanone; 1-(2,5-dimethylphenethyl) ethanone; Dimethyl-1H-indene; Diethylmalonic acid; 3-chlorobenzylnonyl ester
    Pyrolysis products which found only in CNSs-BA/PET 1,5-hexadiyne; Dimethylamine; Nitrous oxide; 1,1'-(1,4-phenylene)bis-acetophenone; 2-methylindene; Azobenzene; Benzene; (1-methyl-2-cyclopropen-1-yl)-2-methylindene; Stilbene; Ethylketone; 1-(3,4-dimethylphenyl); 1-(4-methylphenyl); 1-ethenyl-4-methylbenzene; Dibenzofuran; 2-naphthol; 4-hydroxy-1,2,3,4-tetrahydrophenanthrene; 9,10-dihydrophenanthrene; Benzopropiophenone; Fluorene; 4-vinylbiphenyl; 1,2,3,4-tetrahydrofil; 9,10-dihydrofil; 4-vinylbiphenyl; 1,4-vinylbiphenyl; Phenylacetone; 1,3,5-cycloheptatriene; 2-phenylnaphthalene; 1-acrylbenzene; 2-methylnaphthalene; 4-(2-benzoyl-5-phenyl-3-thienyl)-methylbenzoic acid; 1,3-dimethyl-1H-indene; Tricyclohexen-8-ol; Hexaethylcyclohexane; 9-phenyl-9-fluorenol; Ethylene oxide; Methoxyphenyltricyclohexadecen-5-ylmethanol; 4-benzylbiphenyl; Tritylbenzene; 9-phenylanthracene; 3-(1-phenylethoxy)-3H-isobenzofuran-1-one; 4-phenyl-3,4-dihydroisoquinoline; Oxetane; 2-phenyl; 3-phenylethynyl; Tetraphenyl; 1-[4-(2-phenylethenyl)phenyl]-ethanone; Acenaphthene; 1,2,3,5-tetraisopropyl-cyclohexane; 6,9-dimethoxy-phenazine-1-carboxylic acid; [1,1'-biphenyl]-4-yl-phenylmethanone; 1,1':4',1''-3'-methyltriphenylene
    Pyrolysis products which found both in PET and CNSs/PET Acetophenone; Benzoic acid; Biphenyl; 2-methyl-1,1'-biphenyl; 1,1'-(1,4-phenylene)bisacetophenone; p-terphenyl
    Pyrolysis products which found both in PET and CNSs-BA/PET Styrene; Acetophenone; Benzoic acid; Biphenyl; 2-ethyl-1,1'-biphenyl; Benzophenone; 9H-fluoren-9-one; p-terphenyl
    Pyrolysis products which found both in CNSs/PET and CNSs-BA/PET Benzene; Biphenyl; Acetophenone; Naphthalene; Toluene; Phenanthrene; Indene; 6,6-diphenylfulvene; p-terphenyl; Methylstyrene; Biphenylacetophenone; 4-ethylbiphenyl; Diphenylmethane
    Pyrolysis products found in PET, CNS/PET and CNSs-BA/PET Acetophenone; Benzoic acid; Biphenyl; p-terphenyl
    下载: 导出CSV 
    | 显示表格

    图10为PET阻燃复合材料的抗拉强度和断裂伸长率图。可知,随着阻燃剂含量的增加,CNSs/PET和CNSs-BA/PET复合材料的抗拉强度和断裂伸长率均呈下降趋势。尤其是当阻燃剂含量超过2.0wt%时,PET复合材料的抗拉强度和断裂伸长率大幅度下降。这是由于高含量的阻燃剂在PET基体中形成了较大的团聚体,破坏了基体的连续性,阻碍了应力的传递所致,后续研究中应重点关注材料力学性能的改善。

    图  10  PET复合材料的抗拉强度(a)和断裂伸长率(b)
    Figure  10.  Tensile strength (a) and elongation at break (b) of PET composites

    (1)为同时改善聚对苯二甲酸乙二醇酯(PET)的阻燃性和抗熔滴性,在碳纳米球表面接枝4-苯亚甲基氨基苯酚制备了一种新型碳基复合阻燃剂(CNSs-BA)。CNSs-BA为粒径约50 nm的球形颗粒,热稳定性良好。

    (2) CNSs-BA的引入可显著提高PET的阻燃性和抗熔滴性。当CNSs-BA添加量为2.0wt%时,CNSs-BA/PET复合材料的极限氧指数(LOI)从PET的21.0%提高至28.1%,阻燃等级达到V-0级,热释放速率峰值降低了46.3%。

    (3) CNSs-BA/PET表现出典型的凝聚相阻燃机制。CNSs-BA的引入能促进PET成炭,CNSs-BA/PET的高温残炭量(CR500℃)比PET提高了55.4%,且成炭量的实际值大于理论值。与纯PET的炭层相比,CNSs-BA/PET燃烧生成的炭层的致密性、连续性以及热稳定性都显著提高。这是由于CNSs-BA的引入促进了PET的高温交联,使其高温降解生成了更多的难燃性焦炭物质。

    (4)本文为碳基阻燃剂的发展提供了重要理论补充,对开发无卤、阻燃、抗熔滴的PET材料具有一定的指导意义。

  • 图  1   复合材料直热模具综合传热示意图

    Figure  1.   Schematic diagram of comprehensive heat transfer of composite material direct heating mold

    图  2   加热管电阻随温度变化的曲线

    Figure  2.   Curve of heating tube resistance with temperature

    图  3   12层铺层的M21层合板真实放热功率

    Figure  3.   True heat release power of 12-layer M21 laminate

    图  4   复合材料层合板热导率计算串联模型和并联模型

    Figure  4.   Thermal conductivity calculation series model and parallel model of composite laminate

    图  5   复合材料直热模具有限元模型

    Figure  5.   Finite element model of composite material direct heating mold

    图  6   复合材料固化直热模具尺寸和加热管排布间距示意图

    Figure  6.   Schematic diagram of mold size and heating tube arrangement spacing of composite material curing direct heat mold

    图  7   复合材料固化直热模具表面测温点排布

    Figure  7.   Arrangement of temperature measuring points on mold surface of composite material curing direct heat mold

    图  8   复合材料固化直热模具各测温点温度随时间变化情况

    Figure  8.   Temperature of each temperature measurement point changes with time of composite material curing direct heat mold

    图  9   初始加热管设置的模具表面温度云图

    Figure  9.   Mold surface temperature cloud map of initial heating tube setting

    图  10   更换边缘加热管间距和功率后模具表面温度云图

    Figure  10.   Mold surface temperature cloud map after changing the edge heating tube spacing and power

    图  11   横截面a、b、c上模具表面相应的温度值

    Figure  11.   Corresponding temperature value of the mold surface on the cross section a, b, c

    图  12   正交试验优化后模具表面温度云图

    Figure  12.   Mold surface temperature cloud map after orthogonal experiment optimization

    图  13   正交试验优化后横截面a、b、c上模具表面相应的温度值

    Figure  13.   Corresponding temperature value of the mold surface on cross-section a, b, c after optimization by orthogonal test

    图  14   考虑接触热阻时模具表面测温点随时间变化情况

    Figure  14.   Temperature measurement point of the mold surface changes with time when considering the contact thermal resistance

    图  15   未涂抹导热硅脂和涂抹导热硅脂时模具表面温度云图

    Figure  15.   Mold surface temperature cloud map when thermal grease is not applied and when thermal grease is applied

    图  16   无层合板和有层合板且厚度不同时模具表面温度云图

    Figure  16.   Mold surface temperature cloud map when there is no laminate and when there is a laminate and the thickness is different

    表  1   干空气热物理性质

    Table  1   Thermal physical properties of dry air

    Temperature/℃λ/(102W(mK)1)v/(106m2s1)Pr
    20 2.59 15.06 0.703
    60 2.90 18.97 0.696
    90 3.13 22.10 0.690
    Notes: λ—Thermal conductivity of dry air; v—Acceleration of dry air; Pr—Prandtl Number.
    下载: 导出CSV

    表  2   M21预浸料的基本参数

    Table  2   Basic parameters of M21 prepreg

    ρc/(kg·m−3)Cc/(J·(kg·K)−1)hc/mmVc/wt%kc/(W·(m·K)−1)
    1580 926 0.184 Vf=59.2
    Vr=40.8
    λP=2.56
    λT=1.48
    Notes: ρc—Density of M21 prepreg; Cc—Specific heat capacity of M21 prepreg; hc—Thickness of every piece of prepreg; Vc—Volume fraction; Vf—Volume fraction of carbon fiber; Vr—Volume fraction of epoxy resin; kc—Thermal conductivity of M21 prepreg; λP—Thermal conductivity of the prepreg along the fiber direction; λT—Thermal conductivity of the prepreg perpendicular to the fiber direction.
    下载: 导出CSV

    表  3   结构钢的物性参数

    Table  3   Physical parameters of structural steel

    ρm/
    (kg·m−3)
    Cm/
    (J·(kg·K)−1)
    km/(W·(m·K)−1)
    20℃120℃205℃
    7850 434 60.4 55.9 52.1
    Notes: ρm—Density of structural steel; Cm—Specific heat capacity of structural steel; km—Thermal conductivity of structural steel.
    下载: 导出CSV

    表  4   加热管具体参数

    Table  4   Heating tube specific parameters

    L1/mmL2/mmL3/mmP1/WP2/WP3/W
    336640300300300
    下载: 导出CSV

    表  5   正交实验前加热管设置

    Table  5   Heating tube setting before orthogonal experiment

    L1/mmL2/mmL3/mmP1/WP2/WP3/W
    336625300300400
    下载: 导出CSV

    表  6   正交试验因素水平

    Table  6   Orthogonal test factor level

    LevelFactor
    L1/mmL2/mmP1/WP2/W
    13570350250
    23966300300
    34362250350
    下载: 导出CSV

    表  7   正交试验方案及结果

    Table  7   Orthogonal test plan and results

    Test numberFactorMaximum temperature
    difference/℃
    L1/mmL2/mmP1/WP2/W
    1 1 1 1 1 6
    2 1 2 2 2 5.11
    3 1 3 3 3 4.82
    4 2 1 2 3 4.65
    5 2 2 3 1 3.51
    6 2 3 1 2 6.39
    7 3 1 3 2 4.51
    8 3 2 1 3 5.46
    9 3 3 2 1 3.95
    下载: 导出CSV

    表  8   正交实验极差分析

    Table  8   Orthogonal experiment range analysis

    VariableFactor
    L1/mmL2/mmP1/WP2/W
    K1 15.93 15.16 17.85 13.46
    K2 14.55 14.08 13.71 16.01
    K3 13.92 15.16 12.84 14.93
    Excellent level 3 2 3 1
    R 2.01 1.08 5.01 2.55
    Order C>D>A>B
    Notes: Ki(i=1,2,3)—Sum of the i-th level test indicators of each factor, R—Extremely poor, that is, the difference between the maximum value and the minimum value of Ki in each column.
    下载: 导出CSV
  • [1] 沈军, 谢怀勤. 先进复合材料在航空航天领域的研发与应用[J]. 材料科学与工艺, 2008(5):737-740. DOI: 10.3969/j.issn.1005-0299.2008.05.036

    SHENG Jun, XIE Huaiqin. Research, development and application of advanced composite materials in the aerospace field[J]. Materials Science and Technology,2008(5):737-740(in Chinese). DOI: 10.3969/j.issn.1005-0299.2008.05.036

    [2] 朱晋生, 王卓, 欧峰. 先进复合材料在航空航天领域的应用[J]. 新技术新工艺, 2012(10):76-79. DOI: 10.3969/j.issn.1003-5311.2012.10.025

    ZHU Jinshen, WANG Zhuo, OU Feng. Application of advanced composite materials in aerospace[J]. New Technology & New Process,2012(10):76-79(in Chinese). DOI: 10.3969/j.issn.1003-5311.2012.10.025

    [3] 杜善义. 先进复合材料与航空航天[J]. 复合材料学报, 2007(1):1-12. DOI: 10.3321/j.issn:1000-3851.2007.01.001

    DU Shanyi. Advanced composite materials and aerospace[J]. Acta Materiae Compositae Sinica,2007(1):1-12(in Chinese). DOI: 10.3321/j.issn:1000-3851.2007.01.001

    [4]

    MAIR R I. Advanced composite structures research in Australia[J]. Composite Structures,2002,57(1-4):3-10. DOI: 10.1016/S0263-8223(02)00058-2

    [5] 韩远飞, 孙相龙, 邱培坤. 颗粒增强钛基复合材料先进加工技术研究与进展[J]. 复合材料学报, 2017, 34(8):1625-1635.

    HAN Yuanfei, SUN Xianglong, QIU Peikun. Research and development of processing technology on particulate reinforced titanium matrix composites[J]. Acta Materiae Compositae Sinica,2017,34(8):1625-1635(in Chinese).

    [6] 晏冬秀, 刘卫平, 黄钢华. 复合材料热压罐成型模具设计研究[J]. 航空制造技术, 2012(7):49-52. DOI: 10.3969/j.issn.1671-833X.2012.07.006

    YAN Dongxiu, LIU Weiping, HUANG Ganghua. Design and research of molding die for composite material autoclave[J]. Aeronautical Manufacturing Technology,2012(7):49-52(in Chinese). DOI: 10.3969/j.issn.1671-833X.2012.07.006

    [7] 唐爱霞. 拉挤玻璃钢的成型工艺及市场开发前景[J]. 玻璃钢/复合材料, 1996(2):52.

    TANG Aixia. The forming process and market development prospect of pultruded FRP[J]. Fiber Reinforced Plastics/Composite,1996(2):52(in Chinese).

    [8] 赵洪凯, 钱春香, 周效谅. 拉挤工艺成型连续纤维增强热塑性FRP的性能研究[J]. 化学建材, 2006(1):23-25.

    ZHAO Hongkai, QIAN Chunxiang, ZHOU Xiaoliang. Study on properties of continuous fiber reinforced thermoplastic FRP formed by pultrusion process[J]. Green Building,2006(1):23-25(in Chinese).

    [9] 齐俊伟, 宋伟, 肖军. C_f/E预浸料先进拉挤固化工艺[J]. 宇航材料工艺, 2011, 41(5):44-49. DOI: 10.3969/j.issn.1007-2330.2011.05.011

    QI Junwei, SONG Wei, XIAO Jun. Advanced pultrusion curing process of C_f/E prepreg[J]. Aerospace Materials and Technology,2011,41(5):44-49(in Chinese). DOI: 10.3969/j.issn.1007-2330.2011.05.011

    [10] 王健, 郑学丰, 付昌云. 碳纤维/环氧树脂复合材料-铝合金层合板深拉成型特性[J]. 复合材料学报, 2019, 36(12):2786-2794.

    WANG Jian, ZHENG Xuefeng, FU Changyun. Deep drawing characteristics of carbon fiber/epoxy resin compositealuminum alloy laminates[J]. Acta Materiae Compositae Sinica,2019,36(12):2786-2794(in Chinese).

    [11] 花蕾蕾, 安鲁陵, 匡海华. 复合材料构件热压罐成型模具温度均匀性分析[J]. 南京航空航天大学学报, 2019, 51(3):357-365.

    HUA Leilei, AN Luling, KUANG Haihua. Analysis of temperature field uniformity of composite mold in autoclave[J]. Journal of Nanjing University of Aeronautics & Astronautics,2019,51(3):357-365(in Chinese).

    [12] 王宇, 齐俊伟, 王跃全. 多区间热压对拉挤制件成型质量的影响[J]. 玻璃钢/复合材料, 2019(10):10-18.

    WANG Yu, QI Junwei, WANG Yuequan. The influence of multi-zone hot pressing on the forming quality of pultruded parts[J]. Fiber Reinforced Plastics/Composite,2019(10):10-18(in Chinese).

    [13] 齐俊伟, 李勇, 肖军. 先进拉挤成形技术及其在大飞机复合材料结构中的应用[J]. 航空制造技术, 2011(15):58-60. DOI: 10.3969/j.issn.1671-833X.2011.15.012

    QI Junwei, LI Yong, XIAO Jun. Advanced pultrusion technology and its application in the composite structure of large aircraft[J]. Aeronautical Manufacturing Technology,2011(15):58-60(in Chinese). DOI: 10.3969/j.issn.1671-833X.2011.15.012

    [14]

    CHEN S C, JONG W R, CHANG J A. Dynamic mold surface temperature control using induction heating and its effects on the surface appearance of weld line[J]. Journal of Applied Polymer Science,2006,101(2):1174-1180. DOI: 10.1002/app.24070

    [15]

    CHEN S C, PENG H S, CHANG J A. Simulations and verifications of induction heating on a mold plate[J]. International Communications in Heat and Mass Transfer,2004,31(7):971-980. DOI: 10.1016/j.icheatmasstransfer.2004.05.007

    [16]

    CHEN S C, JONG W R, CHANG Y J. Rapid mold temperature variation for assisting the micro injection of high aspect ratio micro-feature parts using induction heating technology[J]. Journal of Micromechanics and Microengineering,2006,16(9):1783-1791. DOI: 10.1088/0960-1317/16/9/005

    [17]

    HUANG M S, HUANG Y L. Effect of multi-layered induction coils on efficiency and uniformity of surface heating[J]. International Journal of Heat and Mass Transfer,2010,53(11-12):2414-2423. DOI: 10.1016/j.ijheatmasstransfer.2010.01.042

    [18] 戴春晖. 大型复合材料模具内置循环水加热系统设计研究[D]. 长沙: 国防科学技术大学, 2007.

    DAI Chunhui. Design and research on the built-in circulating water heating system of large composite material mould[D]. Changsha: National University of Defense Technology, 2007(in Chinese).

    [19] 史红瑞, 李文斌. 喷管模压料成型过程温度控制系统的研究[J]. 机械工程与自动化, 2008(4):129-131, 134. DOI: 10.3969/j.issn.1672-6413.2008.04.048

    SHI Hongrui, LI Wenbing. Research on temperature control system of nozzle molding material[J]. Mechanical Engineering & Automation,2008(4):129-131, 134(in Chinese). DOI: 10.3969/j.issn.1672-6413.2008.04.048

    [20] 李金国, 蒋宁, 高增梁. 反应成型模具耦合温度场数值模拟与试验研究[J]. 机械工程学报, 2014, 50(8):73-80. DOI: 10.3901/JME.2014.08.073

    LI Jingguo, JIANG Ning, GAO Zengliang. Numerical simulation and experimental research on coupled temperature field of reaction molding die[J]. Journal of Mechanical Engineering,2014,50(8):73-80(in Chinese). DOI: 10.3901/JME.2014.08.073

    [21] 姜碧羽, 齐俊伟, 刘小林. 预浸料先进拉挤成型的固化传热过程数值模拟[J]. 复合材料学报, 2020, 37(6):1496-1504.

    JIANG Biyu, QI Junwei, LIU Xiaolin. Numerical simulation of curing and heat transfer process of prepreg in advanced pultrusion[J]. Acta Materiae Compositae Sinica,2020,37(6):1496-1504(in Chinese).

    [22]

    PILLING M, YATES B, BLACK M. The thermal conductivity of carbon fibre-reinforced composites[J]. Journal of Materials Science,1979(14):1326-1338.

    [23]

    SPRINGER S G, TSAI S W. Thermal conductivity of unidirectional materials[J]. Journal of Composite Materials,1967(1):166-173.

    [24]

    ADAMS D, DONER D. Longitudinal shear loading of a unidirectional composite[J]. Journal of Composite Materials,1967(1):4-17.

    [25]

    SAFONOV A A, CARLONE P, AKHATOV I. Mathematical simulation of pultrusion processes: A review[J]. Compo-site Structures,2018(184):153-177.

图(16)  /  表(8)
计量
  • 文章访问数:  1218
  • HTML全文浏览量:  513
  • PDF下载量:  112
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-21
  • 录用日期:  2020-12-13
  • 网络出版日期:  2020-12-22
  • 刊出日期:  2021-08-31

目录

/

返回文章
返回