Processing math: 58%

弱可见光催化活性La-TiO2-还原氧化石墨烯填充聚偏氟乙烯共混膜的构建

黄露露, 张艳玲, 王挺, 吴礼光, 董春颖

黄露露, 张艳玲, 王挺, 等. 弱可见光催化活性La-TiO2-还原氧化石墨烯填充聚偏氟乙烯共混膜的构建[J]. 复合材料学报, 2020, 37(11): 2749-2759. DOI: 10.13801/j.cnki.fhclxb.20200318.001
引用本文: 黄露露, 张艳玲, 王挺, 等. 弱可见光催化活性La-TiO2-还原氧化石墨烯填充聚偏氟乙烯共混膜的构建[J]. 复合材料学报, 2020, 37(11): 2749-2759. DOI: 10.13801/j.cnki.fhclxb.20200318.001
HUANG Lulu, ZHANG Yanling, WANG Ting, et al. Fabrication of polyvinylidene fluoride blending membranes filled by La-TiO2-reduced graphene oxide with photocatalytic activity[J]. Acta Materiae Compositae Sinica, 2020, 37(11): 2749-2759. DOI: 10.13801/j.cnki.fhclxb.20200318.001
Citation: HUANG Lulu, ZHANG Yanling, WANG Ting, et al. Fabrication of polyvinylidene fluoride blending membranes filled by La-TiO2-reduced graphene oxide with photocatalytic activity[J]. Acta Materiae Compositae Sinica, 2020, 37(11): 2749-2759. DOI: 10.13801/j.cnki.fhclxb.20200318.001

弱可见光催化活性La-TiO2-还原氧化石墨烯填充聚偏氟乙烯共混膜的构建

基金项目: 国家自然科学基金(21776250);浙江省自然科学基金(LY20B060001; LY19B060004)
详细信息
    通讯作者:

    王挺,博士,副研究员,研究方向为功能膜材料及其应用 E-mail:zjwtwaiting@hotmail.com

    吴礼光,博士,教授,研究方向为功能膜材料及其应用 E-mail:wulg64@hotmail.com

  • 中图分类号: TQ051.893

Fabrication of polyvinylidene fluoride blending membranes filled by La-TiO2-reduced graphene oxide with photocatalytic activity

  • 摘要: 为提高聚偏氟乙烯(PVDF)超滤膜的通量及抗污染性能,首先利用吸附相反应技术耦合乙醇热处理制备La掺杂TiO2-还原氧化石墨烯(La-TiO2-RGO),再将其与PVDF共混制备La-TiO2-RGO/PVDF抗污染超滤膜。结果表明,均匀分散于PVDF高分子中表面亲水的La-TiO2-RGO增多,La-TiO2-RGO/PVDF共混膜的水通量和抗污染性能也显著提升。当La-TiO2-RGO/PVDF共混膜中出现团聚体,则会削弱其膜通量和抗污染性。在La-TiO2-RGO填充量(与PVDF质量比)为2.0%时,La-TiO2-RGO/PVDF共混膜具有最优纯水通量。La-TiO2-RGO/PVDF共混膜最高纯水通量可达171.5 L·m−2·h−1,是PVDF膜的5倍以上,其通量衰减速率也明显低于PVDF膜。另外,由于La-TiO2-RGO具有可见光催化活性,被污染后的La-TiO2-RGO/PVDF共混膜经过光照处理后用水清洗,其膜通量恢复率较直接用水清洗后的通量恢复率大幅提高;热处理温度升高,La-TiO2-RGO弱可见光活性增强,光照后La-TiO2-RGO/PVDF共混膜通量恢复率变大。但过高热处理温度抑制了La-TiO2-RGO中Ti3+形成,且削弱其光活性,La-TiO2-RGO/PVDF共混膜通量恢复率反而下降;对于La-TiO2-RGO填充量为2.0%的La-TiO2-RGO/PVDF共混膜,被污染后分别采用直接水清洗、仅光照处理2 h、先光照处理2 h后水清洗的膜通量恢复率分别为79.28%、52.42%、90.01%。
    Abstract: To improve the flux and antifouling performance of polyvinylidene fluoride (PVDF) ultrafiltration membrane, La doped TiO2-reduced graphene (La-TiO2-RGO) was first synthesized by the adsorption phase reaction coupled with alcohol solvothermal reduction processes. Followed that, La-TiO2-RGO was blended with PVDF to fabricate a La-TiO2-RGO/PVDF ultrafiltration membrane with high anti-fouling performance. The results show that the water flux and antifouling performance of the La-TiO2-RGO/PVDF blending membrane increase, when more well-distributed La-TiO2-RGO with hydrophilic groups added. The aggregations generate in the La-TiO2-RGO/PVDF blending membranes will depress their water flux and antifouling performance. When the loading content (mass ratio to PVDF) of La-TiO2-RGO is 2.0%, the La-TiO2-RGO/PVDF blending membrane has the best pure water flux. The optimum pure water flux of the La-TiO2-RGO/PVDF blending membrane reaches 171.5 L·m−2·h−1, which is 5 times as high as that of the PVDF membrane. And the flux decay rate of the La-TiO2-RGO/PVDF blending membrane is significantly lower than that of the PVDF membrane. The flux recovery of the contaminated La-TiO2-RGO/PVDF blending membrane treated by illumination then washing is obviously higher than that just treated by washing, due to the addition of La-TiO2-RGO with photocatalytic activity. The increase in the solvothermal temperature enhances the photocatalytic activity of La-TiO2-RGO, thus improving the flux recovery rate of the La-TiO2-RGO/PVDF membranes after light irradiation. However, too high solvothermal temperature inhibites the formation of Ti3+ in La-TiO2-RGO, which weakens its photoactivity and decreases the flux recovery rate of the resulting La-TiO2-RGO/PVDF blending membranes. For the La-TiO2-RGO/PVDF blending membrane filled with La-TiO2-RGO loading of 2.0%, the flux recovery rates of contaminated membrane are 79.28%, 52.42% and 90.01%, respectively, after washing, illumination for 2 h and illumination for 2 h then washing.
  • 纺织复合材料相比传统金属材料拥有较高的比强度与比模量,且抗冲击性优异,有着较强的可设计性,被广泛应用于航空航天领域[1]。2.5D机织复合材料是纺织复合材料的一种[2],其经纱与纬纱之间形成联锁,在厚度方向上以一定角度进行交织,使其相比二维纺织复合材料具有更好的整体性和层间抗剪切性能。而相比三维正交复合材料,2.5D复合材料具有更低廉的制造成本和更短的生产周期,并且更易于制备回转构件,例如头锥、壳体等复杂构件,具有非常广泛的应用前景,得到了国内外学者的密切关注与研究。

    针对2.5D机织复合材料的工程应用和结构设计,需对其真实力学性能和破坏机制有明确的认识。力学试验是获取材料力学性能最直接有效的方法之一。Jian等[3]通过试验研究了2.5D机织复合材料在不同温度下的力学响应、损伤演化和失效行为特征,发现随着温度的升高,材料的弹性模量和强度均呈下降趋势,但纤维束的强度对温度和加载方向均不敏感。Teng等[4]采用数字图像相关法开展了2.5D机织复合材料的面内拉伸试验研究,有效表征了应变的周期性分布特征,发现平均应变的偏差随计算面积的增加而减小,提出计算区域的大小应大于单位区域大小的两倍,且当面积大小是材料单胞大小的整数倍时,可以获得更准确的结果。Yuan等[5]制备了浅交弯联、浅交直联、深交联三种2.5D机织复合材料试样,并分别进行了低速冲击试验,试验结果对比显示三者中浅交弯联结构抗冲击性最佳。Younes等[6]基于三维Tsai-Wu强度准则,针对2.5D机织复合材料代表性体积单元的损伤容限和弹性模量进行了优化设计,并根据试验结果进行了验证。

    解析模型是一种可以快速预测纺织复合材料力学性能的方法,其优点在于计算效率高及弹性性能预测基本准确,且相比于试验成本极低。杨振宇等[7]采用刚度平均法和有限元法分别预测了2.5D编织结构的等效力学性能,结果显示理论解与有限元解一致性较高。Hallal等[8]提出了一种2.5D机织复合材料的改进型解析模型—“三阶段均一化(3SHM)”,并通过与有限元结果的对比证明了模型的有效性与泛用性。Byun等[9]通过确定结构的结构参数与几何参数建立了一种2.5D复合材料理论模型,使用均一化方法和张量变换来预测弹性常数与纤维体积分数,并分析了各种参数对其有效弹性性能的影响。Chen等[10]则提出了一种根据2.5D机织复合材料的微观几何构型参数和结构特性来表达弹性常数的半解析方法,并分析了几何变量对2.5D机织复合材料力学性能的影响。

    作为一种三维纺织结构,2.5D机织复合材料的损伤模式相比二维复合材料更加复杂多样;单纯通过力学试验及断面显微观测等方式难以观察到材料内部的损伤,而通过CT扫描等方式则无法对渐进损伤过程进行实时观测。另一方面,由于其几何构型的复杂性,2.5D复合材料的解析模型的建模极为复杂,且在强度预测和失效行为分析方面仍具有较大的局限性。因此,对2.5D机织复合材料的渐进损伤机制的研究,需借助计算机数值模拟方法进行辅助分析。

    有限元模拟是一种可获取材料内部应力和损伤分布特征的较有效的数值模拟方法,近年来被广泛应用于各种复合材料的等效力学性能预测和渐进损伤行为分析。Dong等[11]建立了一种可以反映真实纤维束轮廓及走势的2.5D机织复合材料的有限元模型,并对模型进行了验证,发现模拟与试验结果相一致。Liu等[12]提出了一种可模拟渐进失效行为的多尺度非均匀有限元建模方法,分析了2.5D机织复合材料在单轴拉伸下的力学行为,成功地预测了该材料的弹性性能、断裂强度和渐进损伤过程。Zhang等[13]通过CT扫描获取2.5D机织复合材料的几何参数,进而建立了一种细观体素网格有限元模型,准确预测了材料的弹性模量、破坏强度及失效模式。Lu等[14]建立了考虑纤维损伤、基体开裂和界面损伤的2.5D角互锁编织复合材料的细观单胞模型,准确地预测了其拉伸应力-应变响应及损伤演化规律,并发现界面损伤对拉伸强度影响很大。

    关于2.5D机织复合材料损伤失效行为的已有研究,大多以拉伸为主,相比之下对其压缩力学性能及失效模式的研究还较匮乏。卢子兴等[15]开展了2.5D机织复合材料的准静态压缩试验,并建立了单层有限元模型模拟了材料的压缩失效过程,但其并未分析简化模型与真实试样的尺寸及结构差异对模拟结果造成的影响。因此本文针对2.5D机织复合材料的压缩力学性能开展了试验测试,并建立了其细观尺度完整单胞的有限元模型,结合复合材料三维失效准则和损伤演化模型,采用ABAQUS/EXPLICIT求解器及VUMAT自定义子程序,模拟了2.5D机织复合材料在经向和纬向压缩载荷下的力学响应,并对材料的压缩损伤失效机制进行了系统的总结分析。

    2.5D机织复合材料使用增强纤维型号为Toray T800的碳纤维,基体材料为Cytec CYCOM® PR520,其力学性能参数分别来自Xing等[16]和Jia等[17]的报道。复合材料试样及结构示意图如图1所示,其最小体积代表性单胞(RVE)结构的尺寸为13.68 mm×9.74 mm×8.54 mm,并以此尺寸切割制备压缩试样。准静态压缩试验按照标准GB/T 1448—2005[18],分别对材料单胞试样的经向和纬向进行压缩试验;压缩试验采用DNS100电子万能试验机(长春试验机研究所),经向和纬向加载速度分别为0.82 mm/min和0.58 mm/min,以此来保证0.001/s的准静态加载应变率。两种工况各重复三次,试验得到的应力-应变曲线重复性较好。因此,本文选取了其中一条代表性曲线进行后续失效行为分析和模型验证。加载过程中,使用配备有远心镜头的V711高速摄像机(美国,Vision Research Inc)进行同步拍摄,并通过DIC技术计算得到试样的真实应变,试验装置系统见图2

    图  1  2.5D机织复合材料试样及结构示意图
    Figure  1.  Sample and structure diagram of 2.5D woven composites
    RVE—Representative volume element
    图  2  准静态试验系统及设备
    Figure  2.  Quasi-static test system and equipment

    细观有限元模型通过在几何上区分复合材料纤维束和纯基体区,可以很好地反映复合材料织物的真实纺织构型;通过分别对纤维束和基体赋予不同的材料属性与失效模型,细观有限元模型非常适用于分析纺织复合材料内部的应力分布和损伤演化规律。本文建立的2.5D机织复合材料细观有限元模型是基于复合材料单胞的真实几何构型之上建立的,且为全六面体网格,不易出现应力集中与单元畸变,计算效率高;使用区分损伤形式的三维失效准则即渐进损伤模型,可以较大程度模拟材料的真实局部变形、应力-应变分布及各种损伤的起始与扩展。

    2.5D角联锁结构是2.5D机织复合材料中最具代表性的结构之一,主要包含浅交直联、浅交弯联和深交联几种结构形式。本文研究的材料属于经向增强浅交直联形式,共有6层经纱层和7层纬纱层。在此结构中,接结经纱层间形成周期性联锁交织,纬纱近似呈现直线状,与经纱共同作用,提高结构厚度方向的力学性能;填充经纱则有效提高了材料整体的纤维体积分数,使纤维的性能得以充分利用。

    本文通过四个步骤建立材料单胞的细观几何模型:(1) 纤维束几何尺寸测量;(2) 纤维束截面选型;(3) 理想单胞几何模型构建;(4) 单胞几何模型修正。

    本文采用计算机图像测定法[19]计算复合材料纤维束的几何尺寸参数:通过体视显微镜分别拍摄试样垂直于经纱和纬纱的两个截面,作为含有经纱和纬纱横截面的一组切片。通过在图像处理软件Photoshop 2018中设置适当的灰度区间,将不同材料区域大致区分开,再通过人工调整完善灰度分区,如图3所示,最后利用下式计算单根纤维束横截面积S

    S=L1L2×M1n×M2 (1)

    其中:L1L2分别是切面的真实长度和宽度;M1为纤维束横截面区域像素点个数;M2为切片像素点总数;n为切片所包含具有完整横截面纤维束的数目。

    模型中经纱和纬纱的纤维束截面形状采用椭圆截面,其中接结经纱与填充经纱尺寸及纤维体积分数相同。椭圆长短轴的真实长度ab测量方法与纤维束截面测量方法类似,计算式如下:

    a=xL1N1 (2)
    b=yL2N1 (3)

    其中: xy分别是其长短轴的像素长度;N1是切面长边的像素长度。对多个纤维束椭圆截面进行计算,最终取平均值作为模型的输入几何参数。

    图  3  计算机图像测定法处理2.5D机织复合材料试样横截面
    Figure  3.  Computer image measuring method for processing sample section of 2.5D woven composites

    纤维束模型横截面的几何参数见表1,在此基础上,结合编织结构的纱线走势,利用三维织物建模软件(Texgen),可建立2.5D机织复合材料单胞模型,如图4(a)所示。可以观察到,此时的几何模型尺寸与真实试样有一定的偏差,且相邻纤维束之间存在一定程度的侵入。这是由于2.5D机织复合材料的制造工艺特点,成型后织物的实际结构与理论结构会有一定的差异,主要表现在纬纱和填充经纱不再为直线形状,接结经纱的波动角会有所改变。因此,需对理想几何模型进行一些修正调整。通过在Texgen中对经纱和纬纱中节点的面外方向坐标进行修改,使单胞模型的厚度及纤维波动角与真实试样一致,最终得到的几何模型如图4(b)所示。

    表  1  纤维束模型截面几何参数
    Table  1.  Section geometric parameters of fiber tow model
    Width/mmHeight/mmSection power
    Warp yarn 1.21 0.47 1
    Weft yarn 3.15 0.67 0.8
    下载: 导出CSV 
    | 显示表格
    图  4  2.5D机织复合材料几何模型
    Figure  4.  2.5D woven composite geometric model

    需要说明的是,由于真实试样中纤维束之间的相互挤压导致了部分纤维束截面的非对称变形和扭转,而本文建立的几何模型无法精确反映这些几何特征。本文的几何模型修正,仍是基于理想模型的近似几何修正,并结合纤维束体积分数的调整,实现纤维总体积分数与真实试样相一致的目的。

    基于单胞的细观几何模型,通过Texgen软件内置功能可生成三维体素有限元网格,见图5。Zhang等[13]在其关于2.5D复合材料拉伸的论文中研究了大、中、小三种尺寸的体素网格对2.5D复合材料力学性能的影响,结果表明,模拟值随网格尺寸降低而逐渐收敛,且中、小网格模拟值与试验相当接近,并基于此结论在研究2.5D复合材料压缩性能时直接使用了中等网格,得到了较好的结果。本文模型的网格相对尺寸正是基于中网格尺寸所选取的,因此认为计算结果收敛。单胞模型单元总数量为80×64×98,所有单元类型均为八节点实体缩减积分单元(C3D8R)。由于压缩试验试样为单个单胞,因此有限元模型无需使用周期性边界条件,通过建立两个刚性平台模型,定义平台和试样的接触,设置相应的加载速度模拟试验平台的匀速压缩加载。

    图  5  2.5D机织复合材料最小体积代表性单胞(RVE)有限元模型网格划分
    Figure  5.  Meshing of 2.5D woven composite representative volume element (RVE) mesoscopic finite element model

    基体材料PR520是一种高韧性的环氧树脂,其压缩力学响应可以用弹塑性本构来简化表征,并添加了基于位移的失效准则;基于Goldberg等[20]的试验数据,对其弹性和塑性响应进行了定义,对其在压缩、剪切、拉伸三种不同应力三轴度下的断裂应变分别设定为0.07、0.1、0.05。

    2.5D机织复合材料的纤维束由碳纤维丝束和环氧树脂组成,在不考虑纤维束波动的前提下,可将其视为横观各向同性材料。纤维束轴向为1方向,垂直于纤维轴向为2和3方向,其刚度和强度参数可以用Huang[21]的桥联模型计算得到。其中,单向复合材料强度预测公式如下[22]

    轴向拉伸强度:

    F1t=(VfEf11+VmEm)Ef11σft (4)

    轴向压缩强度:

    F1c=(VfEf11+VmEm)Ef11σfc (5)

    横向拉伸强度:

    F2t=(Vf+0.4Vm)Ef22+0.6VmEm(0.4Ef22+0.6Em)Kt22σmt (6)

    横向压缩强度:

    F2c=(Vf+0.4Vm)Ef22+0.6VmEm(0.4Ef22+0.6Em)Kc22σmc (7)

    面内剪切强度:

    Fls=(Vf+αVm)Gf12+(1α)VmGmαGf12+(1α)Gmσms (8)

    式中:VfVm为纤维和基体体积分数;Ef11Ef22Gf12为纤维丝的纵向、横向和剪切模量;Em为基体的弹性模量;σftσfcσmtσmcσms分别为纤维丝的拉伸、压缩强度、基体的拉伸、压缩、剪切强度;Kt22Kc22分别为基体的横向拉伸应力集中系数和横向压缩应力集中系数,α为桥联参数,一般可在0.4到0.5之间取值,本文取0.45。表2为T800/PR520复合材料经纱和纬纱的刚度和强度性能参数。

    表  2  T800碳纤维/PR520环氧树脂纤维束的刚度与强度性能参数
    Table  2.  Stiffness and strength properties of T800 carbon fiber/PR520 epoxy resin fiber tows
    Warp yarnWeft yarn
    Fiber volume fraction Vf /% 81 64
    Elastic modulus E11/GPa 238.96 169.26
    Elastic modulus E22(E33)/GPa 14.59 9.30
    Shear modulus G12(G13)/GPa 7.63 3.62
    Shear modulus G23/GPa 4.97 2.93
    Tensile strength F1t/MPa 4 462 3160
    Compression strength F1c/MPa 1 747 1 237
    Tensile strength F2t/MPa 59 58
    Compression strength F2c/MPa 119 130
    Shear strength Fls/MPa 107 93
    下载: 导出CSV 
    | 显示表格

    本文采用的纤维束失效准则基于Hashin[23]和Hou等[24]提出的失效判据改进得到:

    纤维拉伸失效(σ11>0)

    eft=(σ11F1t)2+a(σ12+σ13Fls)2 (9)

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

    {e_{{\rm{fc}}}} = {\left( {\frac{{{\sigma _{{\rm{11}}}}}}{{{F_{{\rm{1c}}}}}}} \right)^2} \geqslant 1 (10)

    基体拉伸失效\left( {{\sigma _{{\rm{22}}}} > 0} \right)

    {e_{{\rm{mt}}}} = {\left( {\frac{{{\sigma _{{\rm{22}}}}}}{{{F_{{\rm{2t}}}}}}} \right)^2} + {\left( {\frac{{{\sigma _{{\rm{12}}}}}}{{{F_{{\rm{ls}}}}}}} \right)^2} + {\left( {\frac{{{\sigma _{{\rm{23}}}}}}{{{F_{{\rm{ts}}}}}}} \right)^2} \geqslant 1 (11)

    基体压缩失效\left( {{\sigma _{{\rm{22}}}} < 0} \right)

    {e_{{\rm{mc}}}} = \frac{1}{4}{\left( {\frac{{{\sigma _{{\rm{22}}}}}}{{{F_{{\rm{ls}}}}}}} \right)^2} + \frac{{F_{{\rm{2c}}}^{\rm{2}}{\rm{ }}{\sigma _{{\rm{22}}}}}}{{4{F_{{\rm{ls}}}}^2{F_{{\rm{2c}}}}}} - {\left( {\frac{{{\sigma _{{\rm{22}}}}}}{{{F_{{\rm{2c}}}}}}} \right)^2} + {\left( {\frac{{{\sigma _{{\rm{12}}}}}}{{{F_{{\rm{ts}}}}}}} \right)^2} \geqslant 1 (12)

    式中,eftefcemtemc为表征材料失效的变量,当其小于1时代表材料未失效,大于等于1时则表示材料失效。在材料发生失效前,材料表现为线弹性行为。式(9)中系数a为剪切应力对纤维拉伸失效的影响因子。

    一旦某种损伤满足其失效判据,继续加载则会造成材料刚度下降,Zhao等[25]和刘鹏等[26]在文章里对其进行了详细的解释。在该模型的损伤演化表达式中,通过引入特征元素长度(实体单元体积的立方根),可计算单元的断裂能GI

    {G_{\rm{I}}} = \frac{1}{2} \sigma _{{\rm{eq}}}^{\rm{f}}\varepsilon _{{\rm{eq}}}^{\rm{f}}{l_{\rm{c}}} (13)

    其中:\sigma _{{\rm{eq}}}^{\rm{f}}\varepsilon _{{\rm{eq}}}^{\rm{f}}分别为等效破坏强度和等效破坏应变;lc是单元特征长度。基于断裂能和失效时的等效应力可得到相应的等效失效位移,而单元的损伤变量的演化则可基于等效位移的线性关系得到:

    {d_i} = \frac{{\delta _{{{i}},{\rm{eq}}}^{\rm{f}}\left( {{\delta _{i,{\rm{eq}}}} - \delta _{i,{\rm{eq}}}^0} \right)}}{{{\delta _{i,{\rm{eq}}}}\left( {\delta _{i,{\rm{eq}}}^{\rm{f}} - \delta _{i,{\rm{eq}}}^0} \right)}} (14)

    其中:i=ft、fc、mt、mc;\delta _{i,{\rm{eq}}}^{\rm{f}}表示相应失效模式下的等效失效位移;\delta _{i,{\rm{eq}}}^{\rm{0}}表示相应失效模式的等效损伤起始位移;\delta _{i,{\rm{eq}}}^{}则表示相应失效模式的实时等效位移。

    损伤起始后,复合材料性能的渐进退化可通过含损伤变量的二阶对称张量矩阵进行描述,其中di=0表示材料完好,di=1意味材料已完全失效,详细的参数退化方法参见文献[26]。

    以上包含纤维束失效判据和损伤演化的三维渐进损伤模型可以通过ABAQUS用户子程序VUMAT实现,并结合ABAQUS/EXPLICIT进行求解。

    采用所建立的2.5D机织复合材料细观有限元模型,开展了其准静态压缩仿真分析,结合试验结果,对材料内部的渐进损伤行为与失效机制进行了分析与讨论。

    图6为2.5D机织复合材料单胞试样经向和纬向的压缩破坏形貌。由图6(a)可以看出,经向压缩的主要破坏模式包括经纱的横向开裂、纤维束和界面的分层破坏、纬纱边缘部分的压溃损伤及基体开裂,其中经纱的开裂一般随机发生在经纱波动角最大的位置附近,而分层破坏在经纱和纬纱的交织区域最明显。由图6(b)可以看出,纬向压缩的主要破坏模式为纬纱的压溃、纬纱纤维束的断裂及基体开裂;裂纹首先出现在个别纬纱中,随后扩展到纯基体区域,并继续发展至相邻的纬纱中,最终形成宏观剪切裂纹。值得关注的是,经纱损伤只出现在部分经纱中,而剩余经纱则保持了较高的完整性。

    图  6  2.5D机织复合材料RVE试样准静态压缩破坏形貌
    Figure  6.  Quasi-static compression failure morphology of 2.5D woven composite RVE specimen

    图7(a)图7(b)分别为有限元模拟的2.5D机织复合材料单胞试样经向和纬向压缩应力-应变曲线与试验结果的对比。其中应力是通过加载面的总支反力除以试样加载端面积得到,应变为加载方向的位移除以试样初始长度。通过计算在0.1%~0.2%应变范围内应力-应变曲线的斜率,可以得到材料的弹性模量。破坏应力与破坏应变分别为曲线峰值点对应的应力与应变值。在线弹性阶段,并未观察到明显损伤,仿真曲线与试验曲线吻合较好。表3为弹性模量、破坏应力和破坏应变的试验与模拟值对比。

    图  7  2.5D机织复合材料在准静态压缩下有限元模拟与试验应力-应变曲线对比
    Figure  7.  Comparison of quasi-static compression stress-strain curves between finite element modeling and experiment of 2.5D woven composites
    表  3  T800/PR520 2.5D机织复合材料弹性模量、破坏应力、破坏应变的仿真预测和试验测试结果对比
    Table  3.  Comparison of numerical predicted and experimental measured elastic modulus, failure stress and failure strain of T800/PR520 2.5D woven composite
    Warp directionWeft direction
    Elastic modulus/GPaStrength/MPaFailure strain/%Elastic modulus/GPaStrength/MPaFailure strain/%
    Simulation 59.3 300.0 0.50 53.1 268.8 0.52
    Experiment 56.5 308.6 0.51 48.6 211.5 0.49
    Error/% 4.7 2.8 2.0 8.5 21.3 5.7
    下载: 导出CSV 
    | 显示表格

    结合图7表3可见,经向压缩模拟和仿真结果吻合度很好,对弹性模量、破坏应力、破坏应变的预测误差均在5%以内;纬向压缩仿真结果误差相对较大,其中弹性模量和破坏应变的误差较低,分别为8.5%和5.7%,但是破坏强度的误差较大,超过了20%。由图7(b)可见,试验曲线在高应变区间 (0.3%~0.5%)有明显的非线性行为,而仿真曲线未能有效地模拟该非线性响应。这种非线性行为主要是由于制造加工所产生的纤维束扭转,进而可能导致纤维束的提前屈曲失效。然而,由于本文细观模型存在一定局限性,未能重构纤维束的局部扭结行为,难以精确反映其所产生的非线性应力-应变响应,造成了模拟与试验曲线的差异。总体来说,本文建立的细观单胞有限元模型具有较高的准确性,可以用于2.5D机织复合材料的失效模式和渐进损伤行为分析。

    作为一种纤维增强复合材料,2.5D机织复合材料中纤维束的损伤状态直接影响着复合材料的实时力学响应和最终失效。基于验证后的细观有限元模型,可对不同加载状态下纤维束的损伤模式、损伤起始和失效过程进行系统的分析。

    图8为2.5D机织复合材料在0.41%和0.5%经向压缩应变下,纤维束和基体的损伤云图。其中,数值0代表没有损伤,数值1代表完全损伤。图9为复合材料出现损伤前(压缩应变= 0.39%)的应力分布云图。由图8图9可见,填充经纱承受了较高的轴向压缩载荷,并于0.41%应变时最先出现轴向压缩损伤。这主要有两方面的原因:从几何构型和力的传递来看,接结经纱穿插于不同的纬纱层之间,具有较大的纤维束波动角,在上升和下降波动处会产生方向相反的弯矩,使经纱波动角最小的位置外侧产生明显的轴向拉伸应力,内侧产生轴向压缩应力;同时,由于经纱轴向模量远大于横向模量,因而波动较小的填充经纱在经向的模量大于接结经纱,因此轴向填充经纱起主要承力作用。

    图  8  2.5D机织复合材料经向压缩渐进损伤过程
    Figure  8.  Progressive damage process of 2.5D woven composite under warp direction compression
    图  9  2.5D机织复合材料在经向压缩载荷下纤维束的应力分布云图(ε=0.39%)
    Figure  9.  Stress distribution contours of fiber bundles for 2.5D woven composite under warp compression (ε=0.39%)

    此外,材料经向压缩会引起一定程度的纬向变形,该变形会导致经纱和基体及相邻经纱之间的挤压,从而产生横向压缩应力。随着该应力的增加,最终导致在接结经纱波动角最大位置附近出现横向压缩损伤,卢子兴等[15]在其研究中也对这一现象有所描述。在图6(a)中可以很明显地观察到,多根接结经纱断裂的位置几乎都位于波动角最大处,虽然由于试样的微观缺陷与工艺误差,断面分布具有随机性,在厚度方向并不是完全一致,但这仍然很好地证明了模拟结果的有效性。另一方面,接结经纱受到了与之相邻的纬纱的横向挤压,其反作用力在纬纱表面引起较高的横向压缩应力,从而使纬纱在受挤压位置发生局部横向压缩损伤。

    随着压缩载荷的提高,当应变达到0.51%时,纤维束的损伤程度有显著上升。填充经纱中的局部轴向压缩损伤基本扩展至整个纤维束横截面,直接导致了填充经纱的断裂及应力-应变曲线的卸载;接结经纱上的横向压缩损伤也有所提高,但并没有达到完全失效的程度,因此个别接结经纱可能出现断裂,但大部分仍然可以继续承载,这也是经向压缩应力-应变曲线从峰值载荷下降后仍可以保持一定载荷的原因。在达到破坏应变时,纬纱中的横向压缩损伤由表面局部向纤维束内部扩展,具体表现为纬纱发生面外翘曲及横向压溃损伤。总体来说,2.5D机织复合材料在经向准静态压缩下的渐进损伤过程较复杂,其失效是由填充经纱的轴向压缩失效、接结经纱的横向压缩失效及纬纱的横向压溃共同作用而导致的。

    图10为2.5D机织复合材料在纬向压缩载荷下纤维束和基体的渐进损伤过程。图11为纤维束在纬向压缩应变0.35%时的应力分布云图。可见,在纬向压缩条件下,与压头接触的外侧经纱表现出较高的横向压缩应力,而接结经纱中的轴向应力集中在其最大波动角附近;纬纱承受了最大比例的压缩载荷,该载荷导致的轴向压缩应力主要分布于纬纱纤维最小波动角处的受压一侧。在全局应变达到0.38%时,导致沿纬纱分布的横向压缩损伤,此时纤维束和基体几乎没有其他损伤,材料整体压缩刚度基本保持不变。

    图  10  2.5D机织复合材料纬向压缩渐进损伤过程
    Figure  10.  Progressive damage process of 2.5D woven composite under weft direction compression
    图  11  2.5D机织复合材料纬向压缩下的纤维束应力分布云图(ε=0.35%)
    Figure  11.  Yarns stress distribution contours of 2.5D woven composite under weft compression (ε=0.35%)

    随着压缩载荷的持续增加,在应变达到0.52%时,外侧经纱的损伤逐步扩展至整根经纱,内部经纱之间由于相互挤压而发生了明显的横向压缩损伤。另一方面,纬纱承受了大部分的压缩载荷,在其纤维束波动角最大的位置分解为轴向和横向的应力分量,并以轴向压缩损伤和横向压缩损伤的混合形式表现出来,导致了纬纱的断裂和压缩应力-应变曲线的卸载。此外,纬纱的损伤随着纤维束裂纹也扩展至基体中,并贯穿相邻的纬纱,形成明显的裂纹带,这一现象与试验的断裂模式图(图6(b))基本吻合。总体来说,2.5D机织复合材料在纬向准静态压缩下的失效是由纬纱的轴向压缩失效和横向压缩失效主导,并伴随着经纱的横向压缩损伤与基体失效。

    针对2.5D机织复合材料单胞试样分别开展了沿经向和纬向的准静态压缩试验,获得了相应的压缩应力-应变响应与材料破坏模式。另一方面,建立了带有三维失效准则和损伤演化模型的细观单胞有限元模型,并针对2.5D机织复合材料的压缩响应进行了数值模拟,基于模拟结果讨论了其损伤分布与失效机制。得到的主要结论如下。

    (1) 所建立的细观有限元模型能够准确预测2.5D机织复合材料单胞在压缩载荷下的应力-应变曲线,证明了有限元模型的有效性。

    (2) 经向压缩载荷下,2.5D机织复合材料的失效主要是由填充经纱的轴向压缩失效、接结经纱的横向开裂、纤维束间的界面分层破坏及纬纱的横向压溃共同作用而导致的。

    (3) 纬向压缩载荷下,2.5D机织复合材料的失效主要是由纬纱的压溃和纬纱纤维束的断裂主导,并伴随有经纱的横向压缩失效与基体开裂的发生。

  • 图  1   不同温度热处理后La-TiO2-RGO的HRTEM图像(内插图为TEM图像)

    Figure  1.   HRTEM images of La-TiO2-RGO treated by solvothermal reduction under different temperatures (Inset are TEM images)

    图  2   不同温度热处理后La-TiO2-RGO的XRD图谱

    Figure  2.   XRD patterns of La-TiO2-RGO treated by solvothermal reduction under different temperatures

    图  3   不同温度热处理后La-TiO2-RGO的FTIR图谱

    Figure  3.   FTIR spectra of La-TiO2-RGO treated by solvothermal reduction under different temperatures

    图  4   不同温度热处理后La-TiO2-RGO的XPS图谱

    Figure  4.   XPS spectra of La-TiO2-RGO treated by solvothermal reduction under different temperatures

    图  5   不同温度热处理后La-TiO2-RGO可见光激发光电流响应曲线

    Figure  5.   Visible light excitation photocurrent response curves of La-TiO2-RGO treated by solvothermal reduction under different temperatures

    图  6   PVDF膜和La-TiO2-RGO/PVDF共混膜表面的SEM图像

    Figure  6.   SEM images of surface of PVDF membrane and La-TiO2-RGO/PVDF blending membranes

    图  7   PVDF膜和La-TiO2-RGO/PVDF共混膜断面的SEM图像

    Figure  7.   SEM images of cross-section of PVDF membrane and La-TiO2-RGO/PVDF blending membranes

    图  8   光照对PVDF膜和La-TiO2-RGO(170)/PVDF-2混合膜通量衰减过程的影响

    Figure  8.   Effect of irradiation on flux decay process of PVDF membrane and La-TiO2-RGO(170)/PVDF-2 blending membrane

    图  9   PVDF膜和及La-TiO2-RGO/PVDF共混膜的抗污染性能

    Figure  9.   Anti-fouling performance of PVDF membrane and La-TiO2-RGO/PVDF blending membranes

    表  1   La-TiO2-还原氧化石墨烯(RGO)/聚偏氟乙烯(PVDF)共混膜各组分含量

    Table  1   Contents of components of La-TiO2-reduced graphene oxide(RGO)/ polyvinylidene fluoride(PVDF) blending membranes

    MembraneLa-TiO2-RGOVDMAc/
    mL
    Mass of
    PVDF/g
    Mass of
    PEG2000/g
    Mass ratio of La-TiO2-
    RGO to PVDF/%
    Emax/
    MPa
    PVDF 50 7.5 2.5 0 75.3
    La-TiO2-RGO(160)/PVDF-1 La-TiO2-RGO(160) 50 7.5 2.5 1.0 86.7
    La-TiO2-RGO(160)/PVDF-2 La-TiO2-RGO(160) 50 7.5 2.5 2.0 98.9
    La-TiO2-RGO(160)/PVDF-3 La-TiO2-RGO(160) 50 7.5 2.5 3.0 88.1
    La-TiO2-RGO(170)/PVDF-1 La-TiO2-RGO(170) 50 7.5 2.5 1.0 87.2
    La-TiO2-RGO(170)/PVDF-2 La-TiO2-RGO(170) 50 7.5 2.5 2.0 101.1
    La-TiO2-RGO(170)/PVDF-3 La-TiO2-RGO(170) 50 7.5 2.5 3.0 89.5
    Notes: Emax—Maximum tensile strength of membrane; La-TiO2-RGO(160)/PVDF-(1/2/3)—Blending membrane in which mass ratio of La-TiO2-RGO(160) to PVDF is 1%, 2%, 3%, respectively, La-TiO2-RGO(170)/PVDF-(1/2/3) is also understood like this; VDMAc—Volume of N,N-dimethyl acetamide.
    下载: 导出CSV

    表  2   PVDF膜和La-TiO2-RGO/PVDF共混膜的孔隙结构、表面亲水角和膜通量

    Table  2   Porosities, mean pore sizes, water contact angles and flux of PVDF membrane and La-TiO2-RGO/PVDF blending membranes

    MembranePorosity/%Mean pore size/nmWater contact angle/(°)Jw/(L·m−2·h−1)Jp/(L·m−2·h−1)
    PVDF 30.2 35.1 91.3 37.9 7.4
    La-TiO2-RGO(160)/PVDF-1 61.3 51.5 66.7 145.3 31.5
    La-TiO2-RGO(160)/PVDF-2 71.1 62.3 65.3 169.2 39.6
    La-TiO2-RGO(160)/PVDF-3 68.8 49.7 61.2 119.1 22.4
    La-TiO2-RGO(170)/PVDF-1 65.3 57.2 56.4 149.2 33.4
    La-TiO2-RGO(170)/PVDF-2 70.7 65.3 55.8 171.5 42.8
    La-TiO2-RGO(170)/PVDF-3 69.7 59.7 58.9 155.1 38.7
    Notes: Jw—Water flux; Jp—Bull serum albumin (BSA) flux.
    下载: 导出CSV

    表  3   PVDF膜和La-TiO2-RGO/PVDF共混膜的抗污染性能

    Table  3   Anti-fouling performance of PVDF membrane and La-TiO2-RGO/PVDF blending membranes

    MembraneR/%RE/%REr/%
    PVDF 32.65 14.19 34.86
    La-TiO2-RGO(160)/PVDF-1 77.40 34.70 82.22
    La-TiO2-RGO(160)/PVDF-2 79.79 48.47 84.51
    La-TiO2-RGO(160)/PVDF-3 63.69 33.55 72.09
    La-TiO2-RGO(170)/PVDF-1 78.59 43.78 86.00
    La-TiO2-RGO(170)/PVDF-2 79.28 52.42 90.01
    La-TiO2-RGO(170)/PVDF-3 66.12 40.87 72.81
    Notes: R—Flux recovery rate of contaminated membrane after washing; RE—Flux recovery rate after illumination for 2 h; REr—Flux recovery rate after illumination for 2 h then washing.
    下载: 导出CSV
  • [1] 朱振亚, 白成玲, 王磊, 等. 氧化石墨烯-氨基酰化酶/聚偏氟乙烯复合膜的制备及特性[J]. 复合材料学报, 2019, 36(11):2495-2501.

    ZHU Zhenya, BAI Chengling, WANG Lei, etal. Preparation and characteristic of graphene oxide-acylase/poly(vinylidenefluoride) composite membrane[J]. Acta Materiae Compositae Sinica,2019,36(11):2495-2501(in Chinese).

    [2] 朱振亚, 白成玲, 王磊, 等. 磺化氧化石墨烯/聚砜复合膜的制备及抗污染性能[J]. 复合材料学报, 2019, 36(11):2515-2521.

    ZHU Zhenya, BAI Chengling, WANG Lei, et al. Preparation and antifouling property of sulfonated graphene oxide/polysulfone composite membrane[J]. Acta Materiae Compositae Sinica,2019,36(11):2515-2521(in Chinese).

    [3] 冯雪婷, 杨盛, 文晨, 等. Ag2CO3@PVDF/氧化石墨烯超滤膜及其分离性能[J]. 化工学报, 2017, 68(5):2169-2176.

    FENG Xueting, YANG Sheng, WEN Chen, et al. Ag2CO3@PVDF/GO ultrafiltration membrane for water purification[J]. CIESC Journal,2017,68(5):2169-2176(in Chinese).

    [4]

    VENAULT A, CHOU Y N, WANG Y H, et al. A combined polymerization and self-assembling process for the fouling mitigation of PVDF membranes[J]. Journal of Membrane Science,2018,547:134-145. DOI: 10.1016/j.memsci.2017.10.040

    [5]

    SUN C G, FENG X S. Enhancing the performance of PVDF membranes by hydrophilic surface modification via amine treatment[J]. Separation and Purification Technology,2017,185:94-102.

    [6]

    OTITOJU T A, AHMAD A L, OOI B S. Polyvinylidene fluoride (PVDF) membrane for oil rejection from oily wastewater: A performance review[J]. Journal of Water Process Engineering,2016,14:41-59. DOI: 10.1016/j.jwpe.2016.10.011

    [7]

    MA N, CAO J J, LI H Y, et al. Surface grafting of zwitterionic and PEGylated cross-linked polymers toward PVDF membranes with ultralow protein adsorption[J]. Polymer,2019,167:1-12.

    [8]

    RANA D, MATSUURA T. Surface modifications for antifouling membranes[J]. Chemical Reviews,2010,110(4):2448-2471. DOI: 10.1021/cr800208y

    [9]

    FARAHANIA M H D A, VATANPOUR V. A comprehensive study on the performance and antifouling enhancement of the PVDF mixed matrix membranes by embedding different nanoparticulates: Clay, functionalized carbon nanotube, SiO2 and TiO2[J]. Separation and Purification Technology,2018,197:372-381. DOI: 10.1016/j.seppur.2018.01.031

    [10]

    BET-MOUSHOUL E, MANSOURANAH Y, FARHADI K, et al. TiO2 nanocomposite based polymeric membranes: A review on performance improvement for various applications in chemical engineering processes[J]. Chemical Engineering Journal,2016,283:29-46.

    [11]

    CHEN F T, SHI X X, CHEN X B, et al. Preparation and characterization of amphiphilic copolymer PVDF-g-PMABS and its application in improving hydrophilicity and protein fouling resistance of PVDF membrane[J]. Applied Surface Science,2018,427:787-797. DOI: 10.1016/j.apsusc.2017.08.096

    [12]

    YUAN H K, REN J. Preparation of poly(vinylidene fluoride) (PVDF)/acetalyzed poly(vinyl alcohol) ultrafiltration membrane with the enhanced hydrophilicity and the anti-fouling property[J]. Chemical Engineering Research and Design,2017,121:348-359. DOI: 10.1016/j.cherd.2017.03.023

    [13]

    LIU F, XU Y Y, ZHU B K, et al. Preparation of hydrophilic and fouling resistant poly(vinylidene fluoride) hollow fiber membranes[J]. Journal of Membrane Science,2009,345(1-2):331-339. DOI: 10.1016/j.memsci.2009.09.020

    [14] 李妍, 周晓吉, 沈舒苏, 等. 一种两亲性共聚物的合成及其对 PVDF膜的改性研究[J]. 膜科学与技术, 2016, 36(6):70-77.

    LI Yan, ZHOU Xiaoji, SHEN Shusu, et al. Synthesis of a new type of amphiphilic copolymer and its effects on the properties of the modified poly-vinylidene fluoride membrane[J]. Membrane Science and Technology,2016,36(6):70-77(in Chinese).

    [15]

    HEGAB H M, ZOU L. Graphene oxide-assisted membranes: Fabrication and potential applications in desalination and water purification[J]. Journal of Membrane Science,2015,484:95-106. DOI: 10.1016/j.memsci.2015.03.011

    [16]

    ONG C S, GOH P S, LAU W J, et al. Nanomaterials for biofouling and scaling mitigation of thin film composite membrane: A review[J]. Desalination,2016,393:2-15. DOI: 10.1016/j.desal.2016.01.007

    [17]

    XU Z W, WU T F, SHI J, et al. Photocatalytic antifouling PVDF ultrafiltration membranes based on synergy of graphene oxide and TiO2 for water treatment[J]. Journal of Membrane Science,2016,520:281-293. DOI: 10.1016/j.memsci.2016.07.060

    [18]

    DAMODAR R A, YOU S J, CHOU H H. Study theself cleaning, antibacterial and photocatalytic properties of TiO2 entrapped PVDF membranes[J]. Journal of Hazardous Materials,2009,172(2-3):1321-1328. DOI: 10.1016/j.jhazmat.2009.07.139

    [19]

    MOGHADAM M T, LESAGE G, MOHAMMADI T, et al. Improved antifouling properties of TiO2/PVDF nanocomposite membranes in UV-coupled ultrafiltration[J]. Journal of Applied Polymer Science,2015,132(21):41731.

    [20]

    LIAO C J, YU P, ZHAO J Q, et al. Preparation and characterization of NaY/PVDF hybrid ultrafiltration membranes containing silver ions as antibacterial materials[J]. Desalination,2011,272(1-3):59-65.

    [21]

    WU L G, ZHANG X Y, WANG T, et al. Enhanced performance of polyvinylidene fluoride ultrafiltration membranes by incorporating TiO2/graphene oxide[J]. Chemical Engineering Research and Design,2019,141:492-501.

    [22] 许智勇, 李冰蕊, 潘家豪, 等. TiO2复合催化剂弱光催化降解模拟海水中苯酚及其催化活性的影响[J]. 环境科学学报, 2017, 37(12):4593-4601.

    XU Zhiyong, LI Bingrui, PAN Jiahao, et al. Photodegradation of phenol in artificial seawater by TiO2 composite catalysts under weak UV irradiation[J]. Acta Scientiae Circumstantiae,2017,37(12):4593-4601(in Chinese).

    [23]

    LI J F, XU Z L, YANG H. Microporous polyethersulfone membranes prepared under the combined precipitation conditions with non-solvent additives[J]. Polymers for Advanced Technologies,2008,19(4):251-257. DOI: 10.1002/pat.982

    [24]

    FENG C S, SHI B L, LI G M, et al. Preparation and properties of microporous membrane from poly(vinylidene fluoride-co-tetrafluoroethylene) (F2.4) for membrane distillation[J]. Journal of Membrane Science,2004,237(1-2):15-24. DOI: 10.1016/j.memsci.2004.02.007

    [25]

    WANG T, ZHANG Y L, PAN J H, et al. Hydrothermal reduction of commercial P25 photocatalysts to expand their visible-light response and enhance their performance for photodegrading phenol in high-salinity wastewater[J]. Applied Surface Science,2019,480:896-904. DOI: 10.1016/j.apsusc.2019.03.052

    [26]

    VATANPOUR V, MADAENI S S, MORADIAN R, et al. Fabrication and characterization of novel antifouling nanofiltration membrane prepared from oxidized multiwalled carbon nanotube/polyethersulfone nanocomposite[J]. Journal of Membrane Science,2011,375(1-2):284-294.

    [27]

    SILVA T L S, MORALES-TORRES S, FIGUEIREDO J L, et al. Multi-walled carbon nanotube/PVDF blended membranes with sponge- and finger-like pores for direct contact membrane distillation[J]. Desalination,2015,357:233-245. DOI: 10.1016/j.desal.2014.11.025

  • 期刊类型引用(7)

    1. 唐宇哲,傅梅珍,刘恒. UHPC-NC组合梁抗弯承载能力影响因素分析. 公路工程. 2025(01): 41-49+120 . 百度学术
    2. 金清平,陈川,黄桥连,高永红,雷万雄,尹少龙. 超高性能混凝土中纤维作用效应研究综述. 水利与建筑工程学报. 2024(02): 55-65 . 百度学术
    3. 乔险涛,喻鹏,陈喜健,周佳乐,张连杰,李健朋. 非晶合金纤维和钢纤维单/混掺增强混凝土开裂及断裂性能研究. 硅酸盐通报. 2024(06): 2137-2148 . 百度学术
    4. 王冲,刘国磊,何国志,李公成,马秋峰,刘树龙,李庆国. 全尾砂固废基充填材料胶凝特性与纤维强化效应研究. 矿业研究与开发. 2024(06): 112-120 . 百度学术
    5. 刘志峰,卢涛,包杰,王祥雨,黄相又. 无粘结预应力UHPC梁剪切性能试验研究. 世界桥梁. 2024(05): 69-75 . 百度学术
    6. 梁林,王秋维,史庆轩. 粗骨料UHPC单轴拉伸性能及损伤本构模型. 建筑材料学报. 2024(09): 806-815 . 百度学术
    7. 郭旭,张品乐,胡静,查叶铭,陶忠. 钢-PVA混杂纤维水泥基力学性能与工程经济评价. 材料导报. 2024(S2): 179-184 . 百度学术

    其他类型引用(4)

图(9)  /  表(3)
计量
  • 文章访问数:  2045
  • HTML全文浏览量:  699
  • PDF下载量:  45
  • 被引次数: 11
出版历程
  • 收稿日期:  2019-12-22
  • 录用日期:  2020-02-26
  • 网络出版日期:  2020-03-18
  • 刊出日期:  2020-11-14

目录

/

返回文章
返回