• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

长杆高速侵彻的特征控制参量分析

尹志勇 陈小伟

尹志勇, 陈小伟. 长杆高速侵彻的特征控制参量分析[J]. 爆炸与冲击, 2021, 41(2): 023302. doi: 10.11883/bzycj-2020-0057
引用本文: 尹志勇, 陈小伟. 长杆高速侵彻的特征控制参量分析[J]. 爆炸与冲击, 2021, 41(2): 023302. doi: 10.11883/bzycj-2020-0057
YIN Zhiyong, CHEN Xiaowei. Analysis of characteristic control parameters of long-rod penetration[J]. Explosion And Shock Waves, 2021, 41(2): 023302. doi: 10.11883/bzycj-2020-0057
Citation: YIN Zhiyong, CHEN Xiaowei. Analysis of characteristic control parameters of long-rod penetration[J]. Explosion And Shock Waves, 2021, 41(2): 023302. doi: 10.11883/bzycj-2020-0057

长杆高速侵彻的特征控制参量分析

doi: 10.11883/bzycj-2020-0057
基金项目: 国家自然科学基金(11627901,11872118)
详细信息
    作者简介:

    尹志勇(1998- ),男,硕士研究生,3120200140@bit.edu.cn

    通讯作者:

    陈小伟(1967- ),男,博士,教授,chenxiaoweintu@bit.edu.cn

  • 中图分类号: O389

Analysis of characteristic control parameters of long-rod penetration

  • 摘要: 针对理想长杆侵彻,通过对长杆侵彻Alekseevskii-Tate模型近似解进行分析,指出单一的无量纲速度衰减系数α(deceleration index)不足以完全表征长杆高速侵彻的准定常阶段。在此基础上,重新定义了2个无量纲特征参量:Johnson破坏数ΦJp和特征时间系数β,2个参量之间的关系为α=β/ΦJp。分析表明,ΦJpβ(或αβ)可实现对长杆高速侵彻准定常阶段的弹尾速度的完全表征;若再引入长杆弹相对临界速度vc*,则可完全表征长杆侵彻的准定常阶段。此外,还证明了α能够判定侵彻过程偏离定常状态的程度,并指出通过确定ΦJpβ(或αβ),可针对攻防需求对长杆弹侵彻设计进行指导。
  • 长杆弹是一类由钨合金和贫铀合金等高密度金属制成的大长径比(L/D10)动能武器。不同于刚性弹,高速下长杆弹体与靶体相互作用时,作用面上压力远高于材料强度,弹靶发生严重质量侵蚀。由于侵彻机理的独特性和军事应用需求,长杆高速侵彻问题已成为穿甲侵彻领域的研究热点。

    长杆高速侵彻(半流体侵彻)的速度范围大致为1.5~3.0 km/s,不同弹靶材料对应的侵彻速度范围存在差异。在超过3.0 km/s(不同弹靶材料有差异)的超高速碰撞下,将发生完全的流体侵彻,弹靶强度影响可以忽略,同时需要考虑冲击波和可压缩性等因素。在低于1.5 km/s的低速撞击下,长杆将以刚性弹的方式侵彻或无法侵彻(形成界面击溃)。在临界速度范围内,对应不同的弹靶强度关系,将发生两类典型的侵彻模式转变[1]

    与其他工程研究领域类似,长杆高速侵彻的研究包括实验研究、数值模拟和理论分析3个方面。根据实验与模拟所得到的结果进行抽象与分析,建立能够反映实验与模拟中典型现象的物理模型,一直是长杆高速侵彻研究中的重点与难点。长杆高速侵彻的理论模型经历了从早期的流体动力学理论[2]到经典的Alekseevskii-Tate模型[3-5]再到更复杂模型的发展。Alekseevskii-Tate模型是长杆高速侵彻最经典且仍最常用的理论分析模型。

    长杆高速侵彻通常包括初始瞬态、主要侵彻、次级侵彻和靶体回弹等4个阶段。其中主要侵彻阶段对侵彻影响最显著,其最大特征就是准定常,弹靶作用力与速度不会发生显著变化,又称这个阶段为准定常侵彻阶段。包括Alekseevskii-Tate模型在内的几乎所有理论模型都只针对该阶段进行分析,也即通常用准定常侵彻阶段代表长杆高速侵彻的全过程[6]。为了描述侵彻过程准定常阶段与定常状态之间的差异,Jiao等[7-8]首先对Alekseevskii-Tate模型进行近似处理,得到显式的近似理论解析解,并在此基础上定义了无量纲速度衰减系数α。Jiao等[7]认为α可反映侵彻过程中弹尾速度的总衰减程度,可作为判断长杆高速侵彻不同状态的依据。但事实上,仅由速度衰减程度α不足以完全表征长杆高速侵彻的准定常侵彻阶段,本文中将对此作进一步分析。

    Alekseevskii-Tate模型[3-5]是最经典最常用的长杆高速侵彻理论模型。在不可压流体动力学模型[2]的基础上,Alekseevskii-Tate模型在Bernoulli方程中加入弹靶强度项进行修正。其控制方程组如下:

    12ρp(vu)2+Yp=12ρtu2+Rt
    (1)
    ρpldvdt=Yp
    (2)
    dldt=(vu)
    (3)
    dpdt=u
    (4)

    式中:ρpρt分别为弹材和靶材的密度,Yp为弹体强度,Rt为靶体侵彻阻力,pl分别为侵彻深度和弹体剩余长度。需特别指出的是,Alekseevskii-Tate模型是针对理想长杆侵彻,即不考虑弹靶的可压缩性,而弹靶材料的其他特性如应变硬化、应变率硬化和热软化效应等,将非常笼统地体现在弹靶的强度项中。即弹靶的强度项YpRt基本上是由经验确定的。因此,本文中的分析也将在该前提下进行。

    由于方程组(1)~(4)的非线性,无法由Alekseevskii-Tate模型直接显式给出相关物理量在侵彻过程中的变化规律和影响因素。针对其准定常侵彻的特点,Jiao等[8]对Alekseevskii-Tate模型中弹体剩余长度的对数表达式ln(l/L)对时间进行线性近似,获得两组显式的近似理论解,其中一组近似解中速度为常数,这里给出非常数速度的近似解1为:

    vv0=1+2ˉμΦJpKln(1K2ˉμtτ)
    (5)
    uv0=11+μ[1+2ˉμΦJpKln(1K2ˉμtτ)]12μv2c
    (6)
    l=Lμv01+μ[(1+1+μ2μ2v2c2ˉμΦJpK)t+2ˉμΦJpKln(1K2ˉμtτ)(t2ˉμKτ)]
    (7)
    p=v01+μ[(11+μ2μv2c2ˉμΦJpK)t+2ˉμΦJpKln(1K2ˉμtτ)(t2ˉμKτ)]
    (8)

    式中:vc=vc/v0为长杆弹相对临界速度,vc=2|RtYp|/ρp为临界弹体速度,μ=ρt/ρpˉμ=(1μ2)/μ为与弹靶密度比相关的2个无量纲参数,Johnson破坏数ΦJp=ρpv02/Yp为表征结构破坏程度的无量纲参量,τ=L/v0为特征时间(Lv0分别为初始弹长和初始弹速),K为无量纲系数,其取值为:

    K=1μ+ˉμμv2c2+(1μ+ˉμμv2c2)24ˉμ(1μ)ΦJp
    (9)

    在较高的撞击速度下,式(9)中的v2c1/ΦJp将变为可以忽略的小量,式(9)可以近似简化为:

    K˜K=2(1μ)
    (10)

    将式(10)代回上述近似解,可以完全地推导得到Walters等[9]提出的一阶摄动解,即一阶摄动解是近似解在高速下的特殊形式。

    近似解给出了瞬时速度关于时间的显示代数表达,因此可以进一步分析长杆高速侵彻过程中的速度衰减。在此基础上,Jiao等[7]通过对长杆侵彻速度-时间图进行分析,定义了无量纲速度衰减系数α。以表1中一组长杆高速侵彻典型算例[10]为例,图1为弹尾速度在侵彻过程中的变化情况。如图1所示,在无量纲化瞬时弹尾速度曲线上作A(0,1)点的切线并与直线BDtv0/L=Tv0/L=2ˉμ/K,其中T为图1中直线BD所对应的时刻)交于点C,Jiao等[7]定义无量纲速度衰减系数α为:

    表  1  长杆侵彻设计工况中的相关参数
    Table  1.  Related parameters of long-rod penetration design
    v0/(km·s−1)L/mm ρp/(g·cm−3) ρt/(g·cm−3)Yp/GPaRt/GPa
    1.581.717.47.82.04.94
    下载: 导出CSV 
    | 显示表格
    图  1  长杆高速侵彻过程无量纲化速度衰减示意图
    Figure  1.  Schematic of the deceleration of dimensionless tail velocity during long-rod penetration process
    α=BCBD=2ˉμKΦJp
    (11)

    Jiao等[7]认为线段BCBD的比值α能反映弹尾速度的总衰减程度:当α=0时,瞬时弹尾速度在侵彻过程中不发生变化,此即为定常侵彻;当α0时,侵彻过程偏离定常侵彻,α的值实际上刻画了该偏离程度。因此,焦文俊等[1]认为:无量纲速度衰减系数α能够判定侵彻过程偏离定常侵彻状态的程度,可以表征长杆高速侵彻过程。

    α代入式(5),可得:

    vv0=1+αln(1K2ˉμtτ)
    (12)

    类似地,可对式(6)及侵彻速度和长杆侵蚀长度等进行简化,这里不再赘述。由式(12)可知,在控制无量纲弹尾速度v/v0与无量纲时间t/τ的关系中,除无量纲速度衰减系数α之外,还有无量纲系数K/2ˉμ。显然,当α一定时,系数Kˉμ均可发生变化。因此仅由α无法控制长杆高速侵彻的速度变化。

    前文已定义Johnson破坏数ΦJp=ρpv02/YpΦJp值的大小可表征撞击事件的严重程度。在此之外,引入无量纲特征时间系数 β=2ˉμ/K β实质是长杆高速侵彻过程的无量纲时间。则近似解1的式(5)~(8)可简化为:

    vv0=1+βΦJpln(11βtτ)
    (13)
    uv0=1(1+μ)[1+βΦJpln(11βtτ)]12μv2c
    (14)
    lL=1μt(1+μ)τ[(1+1+μ2μ2v2cβΦJp)+βΦJpln(11βtτ)(1βτt)]
    (15)
    pL=t(1+μ)τ[(11+μ2μv2cβΦJp)+βΦJpln(11βtτ)(1βτt)]
    (16)

    由式(13)可知,与无量纲速度衰减系数α的单一表达不同,Johnson破坏数ΦJp和特征时间系数 β共同控制并决定着长杆高速侵彻的弹尾速度变化。

    由式(14)可知,在控制无量纲侵彻速度u/v0与无量纲时间t/τ的关系中,除无量纲系数ΦJp β外,还有无量纲系数1/(1+μ)v2c/(2μ)。由式(10)可知,在较高的撞击速度下,μ β存在关系式:

    μ=1β1
    (17)

    式(17)说明 μ可由 β单一表达。

    由于无量纲系数ΦJp β均不含靶体强度项,因此无法对长杆弹相对临界速度vc进行表征。本文中主要讨论通过控制无量纲系数ΦJp β对弹尾速度v进行描述,因此在此对vc不做过多讨论。事实上,可将长杆弹相对临界速度vc作为第3个无量纲数引入,其与ΦJp β可以共同控制长杆侵彻过程侵彻速度变化。

    类似地,也可对式(15)~(16)进行分析,通过vcΦJp β这3个无量纲数,可以实现对长杆剩余长度(或弹体侵蚀长度)和侵彻深度的完全表征。

    事实上,ΦJp β在无量纲化瞬时弹尾速度曲线中均有明确的物理意义。对式(13)求导,可得:

    d(v/v0)d(t/τ)=1ΦJp1(1/β)(t/τ)1
    (18)

    t=0时,有:

    d(v/v0)d(t/τ)|tτ=0=1ΦJp
    (19)

    由式(19)可知,|1/ΦJp|反映了无量纲化瞬时弹尾速度曲线在初始时刻的斜率,即ΦJp控制着无量纲瞬时弹尾速度曲线在初始时刻的衰减速率,图1|1/ΦJp|即为线段BCAB的比值。由式(18)可知,当t/τ=β时,无量纲弹尾速度曲线斜率为正无穷,即图1D点坐标为 β=2ˉμ/K,因此 β控制着长杆侵彻过程完整的无量纲时间。

    简言之,Johnson破坏数ΦJp控制着弹尾速度曲线初始斜率,特征时间系数 β控制着侵彻时间,因此ΦJp β可以实现对长杆高速侵彻过程中弹尾速度的完全表征。而Jiao等[7]定义的速度衰减程度α图1中为线段BCBD的比值,无法控制侵彻时间,因此仅由α无法完全表征长杆高速侵彻过程。

    实际上,在特征时间系数 β(OD)确定的情况下,α β能够反映弹尾速度曲线初始斜率(BC/OD=α/β)。由式(11)中α=β/ΦJp也可以看出,即长杆侵彻过程弹尾速度既可通过ΦJp β完全表征,也可通过α β完全表征。若考虑弹体侵彻速度以及长杆剩余长度(或弹体侵蚀长度)和侵彻深度,则还需引入长杆弹相对临界速度vc,就可完全表征长杆侵彻的准定常阶段。

    Jiao等[7]通过无量纲速度衰减系数α判定不同侵彻状态的依据为:α值能反映瞬时速度曲线与坐标轴所围成区域的形状特征,α值越趋近于0,形状越接近矩形,代表长杆侵彻过程越接近定常侵彻。但该判断依据仅从观察得到,缺少严格的证明,且仅通过观察很难判断形状接近于矩形的程度。

    为了更准确地描述不同工况下长杆高速侵彻过程偏离定常侵彻的程度,定义长杆高速侵彻的弹尾速度衰减程度为:实际侵彻过程无量纲弹尾速度-时间曲线与坐标轴围成面积与该工况下定常侵彻无量纲弹尾速度-时间曲线与坐标轴围成面积的比值,如图1所示,即扇形OAD面积与矩形OABD的比值。

    下面计算任意工况下长杆高速侵彻的弹尾速度衰减程度。实际侵彻过程中,无量纲弹尾速度-时间曲线与坐标轴围成面积,即为式(5)在t/τ=0t/τ=T/τ=2ˉμ/K的积分:

    S1=2ˉμK(12ˉμKΦJp)
    (20)

    定常侵彻过程中,无量纲的弹尾速度-时间曲线与坐标轴围成面积为:

    S2=2ˉμK
    (21)

    因此,长杆高速侵彻的弹尾速度衰减程度为:

    S1S2=12ˉμKΦJp=1α
    (22)

    这里仍定义2ˉμ/KΦJpα,与Jiao等[7]通过几何关系推导不同,该无量纲系数α通过数学运算推导得到。

    由式(22)及上述推导过程可知,无量纲速度衰减系数α可作为判断长杆高速侵彻不同侵彻状态的依据:当α=0时,S1/S2=1,此时为定常侵彻;当α0时,S1/S21,侵彻过程偏离定常侵彻,且α值越趋近于0,S1/S2越趋近于1代表侵彻过程越接近定常侵彻。

    下面采用表1中算例分析v0ρpρtYpRt对无量纲速度衰减系数α以及Johnson破坏数ΦJp和特征时间系数 β的影响。由前述分析可知,两个无量纲参数ΦJp β(或α β),可完全表征长杆侵彻过程的弹尾速度。若进一步引入长杆弹相对临界速度vc,就可完全表征长杆侵彻的准定常阶段。因此,通过对三个无量纲控制参量的分析,可确定弹靶的相关参数,进而针对攻防需求对长杆弹侵彻设计进行指导。为简单起见,这里仅讨论无量纲参数ΦJp β的影响。

    由Jiao等[7]的分析可知,无量纲速度衰减系数α与初始撞击速度、弹靶密度以及靶体强度负相关,而与弹体强度正相关,且靶体材料性质对α的影响显著小于弹体材料性质。

    ΦJp=ρpv02/Yp可知,无量纲数ΦJp与弹体密度、初始撞击速度成正比,与弹体强度成反比。

    图2 β随相关参数的变化情况,从图2中可以看出, β与初始撞击速度、弹体密度以及弹体强度正相关,而与靶体密度、靶体强度负相关。且从图2(b)(c)中可以看出,当弹、靶密度或强度在相同范围变化时,无量纲数 β的变化相近,即弹、靶的强度和密度分别对 β影响程度相当。

    图  2  无量纲系数 βv0 ρp ρtYpRt的变化情况
    Figure  2.  Variation of dimensionless coefficient  β with v0,  ρp,  ρt, Yp and Rt

    在后续实验中,可在确定两个无量纲控制参数ΦJp β(或α β)的前提下,通过改变v0ρpρtYpRt等参数,设计需要的工况。由式(10)可知,在较高初始撞击速度下,K可用˜K=2(1μ)代替。以表1中算例为例,初始撞击速度分别为1.5、2.25和3 km/s时,对应的˜K/K分别为85.41%、93.33%和96.22%,因此可用˜K代替K,此时α˜α=Ypρ1/2p(ρ1/2p+ρ1/2t)v20 β˜β=ρp/ρt+1,即可忽略靶体强度Rt对长杆高速侵彻过程速度衰减的影响。

    为进一步说明无量纲参数ΦJp β以及无量纲速度衰减系数α的影响因素及相互关系,这里根据Jiao等[7]的参数设计出9组具有不同侵彻状态的工况。相关参数列于表2,侵彻过程中各组工况的瞬时弹尾速度变化情况如图3所示。

    表  2  设计工况中相关参数
    Table  2.  Summary of parameters in the designed cases
    工况v0/
    (km·s−1)
    L/
    mm
     ρp/
    (g·cm−3)
     ρt/
    (g·cm−3)
    Yp/GPaα/%ΦJp β
    11.510019.009.002.011.4821.3752.453
    21.5 5019.009.002.011.4821.3752.453
    33.010019.009.002.02.8785.5002.453
    43.010013.009.002.03.7658.5002.202
    53.010013.005.002.04.4758.5002.612
    61.510019.009.001.05.7442.7502.453
    71.510023.755.562.011.4826.7193.067
    81.510019.0012.562.211.4819.4322.230
    91.510019.009.00002.453
    下载: 导出CSV 
    | 显示表格
    图  3  不同工况下瞬时弹尾速度在侵彻过程中的变化
    Figure  3.  Variation of instant tail velocities during penetration process in different cases

    图3中,工况1和工况2两条曲线完全重合,说明弹体长度对速度衰减程度无影响;对比工况1与工况3可看出,侵彻过程在较高的撞击速度下趋于近定常状态;对比工况3~5可看出,速度衰减程度与ρpρt正相关;对比工况1与工况6可以看出,速度衰减程度与Yp负相关。

    特别需要指出的是,对比工况1与工况7~8可看出,在α不变的前提下,通过选取合适的参数可以改变ΦJp β,即在无量纲速度衰减系数不变的情况下,可改变无量纲弹尾速度初始衰减速率和无量纲衰减时间。这也进一步表明:无量纲速度衰减系数α只能判定侵彻过程偏离定常侵彻的程度,无法完全表征长杆高速侵彻过程。而无量纲数ΦJp β(或α β)通过分别控制初始衰减速率和衰减时间,完全决定着长杆高速侵彻中弹尾速度的变化。

    对于Yp=0的极端情况(工况9),侵彻过程为定常状态。

    对长杆高速侵彻过程速度衰减进行了详细分析。通过对文献[7]中提出的Alekseevskii-Tate模型近似解进行讨论,明确指出无量纲速度衰减系数α不足以完全表征长杆高速侵彻过程。重新提出Johnson破坏数ΦJp和特征时间系数 β(或α β),理论分析其共同控制并完全表征着长杆高速侵彻的弹尾速度变化;若进一步引入长杆弹相对临界速度vc,3个无量纲数就可完全表征长杆侵彻的准定常阶段。通过数学推导,严格证明无量纲速度衰减系数α可作为判定侵彻过程偏离定常侵彻状态程度的依据。结合参数讨论v0ρpρtYpRt等对ΦJp β以及α的影响,指出通过确定无量纲数ΦJp β(或α β),针对攻防需求可对长杆弹侵彻设计进行指导。

  • 图  1  长杆高速侵彻过程无量纲化速度衰减示意图

    Figure  1.  Schematic of the deceleration of dimensionless tail velocity during long-rod penetration process

    图  2  无量纲系数 βv0 ρp ρtYpRt的变化情况

    Figure  2.  Variation of dimensionless coefficient  β with v0,  ρp,  ρt, Yp and Rt

    图  3  不同工况下瞬时弹尾速度在侵彻过程中的变化

    Figure  3.  Variation of instant tail velocities during penetration process in different cases

    表  1  长杆侵彻设计工况中的相关参数

    Table  1.   Related parameters of long-rod penetration design

    v0/(km·s−1)L/mm ρp/(g·cm−3) ρt/(g·cm−3)Yp/GPaRt/GPa
    1.581.717.47.82.04.94
    下载: 导出CSV

    表  2  设计工况中相关参数

    Table  2.   Summary of parameters in the designed cases

    工况v0/
    (km·s−1)
    L/
    mm
     ρp/
    (g·cm−3)
     ρt/
    (g·cm−3)
    Yp/GPaα/%ΦJp β
    11.510019.009.002.011.4821.3752.453
    21.5 5019.009.002.011.4821.3752.453
    33.010019.009.002.02.8785.5002.453
    43.010013.009.002.03.7658.5002.202
    53.010013.005.002.04.4758.5002.612
    61.510019.009.001.05.7442.7502.453
    71.510023.755.562.011.4826.7193.067
    81.510019.0012.562.211.4819.4322.230
    91.510019.009.00002.453
    下载: 导出CSV
  • [1] 焦文俊, 陈小伟. 长杆高速侵彻问题研究进展 [J]. 力学进展, 2019, 49(1): 312–391. DOI: 10.6052/1000-0992-17-021.

    JIAO W J, CHEN X W. Review on long-rod penetration at hypervelocity [J]. Advances in Mechanics, 2019, 49(1): 312–391. DOI: 10.6052/1000-0992-17-021.
    [2] ALLEN W A, ROGERS J W. Penetration of a rod into a semi-infinite target [J]. Journal of the Franklin Institute, 1961, 272(4): 275–284. DOI: 10.1016/0016-0032(61)90559-2.
    [3] ALEKSEEVSKII V P. Penetration of a rod into a target at high velocity [J]. Combustion, Explosion and Shock Waves, 1966, 2(2): 63–66. DOI: 10.1007/BF00749237.
    [4] TATE A. A theory for the deceleration of long rods after impact [J]. Journal of the Mechanics and Physics of Solids, 1967, 15(6): 387–399. DOI: 10.1016/0022-5096(67)90010-5.
    [5] TATE A. Further results in the theory of long rod penetration [J]. Journal of the Mechanics and Physics of Solids, 1969, 17(3): 141–150. DOI: 10.1016/0022-5096(69)90028-3.
    [6] 焦文俊. 长杆高速侵彻的1D和2D理论模型 [D]. 合肥: 中国科学技术大学, 2019.
    [7] JIAO W J, CHEN X W. Analysis of the velocity relationship and deceleration of long-rod penetration [J]. Acta Mechanica Sinica, 2019, 35(4): 852–865. DOI: 10.1007/s10409-019-00862-1.
    [8] JIAO W J, CHEN X W. Approximate solutions of the Alekseevskii-Tate model of long-rod penetration [J]. Acta Mechanica Sinica, 2018, 34(2): 334–348. DOI: 10.1007/s10409-017-0672-9.
    [9] WALTERS W, WILLIAMS C, NORMANDIA M. An explicit solution of the Alekseevski-Tate penetration equations [J]. International Journal of Impact Engineering, 2006, 33(1−12): 837–846. DOI: 10.1016/j.ijimpeng.2006.09.057.
    [10] ANDERSON Jr C E, WALKER J D. An examination of long-rod penetration [J]. International Journal of Impact Engineering, 1991, 11(4): 481–501. DOI: 10.1016/0734-743X(91)90015-8.
  • 期刊类型引用(1)

    1. Ziqi Zhao,Tong Li,Donglin Sheng,Jian Chen,Amin Yan,Yan Chen,Haiying Wang,Xiaowei Chen,Lanhong Dai. Machine learning optimization strategy of shaped charge liner structure based on jet penetration efficiency. Defence Technology. 2024(09): 23-41 . 必应学术

    其他类型引用(0)

  • 加载中
图(3) / 表(2)
计量
  • 文章访问数:  918
  • HTML全文浏览量:  433
  • PDF下载量:  99
  • 被引次数: 1
出版历程
  • 收稿日期:  2020-03-05
  • 修回日期:  2020-06-30
  • 网络出版日期:  2021-02-02
  • 刊出日期:  2021-02-05

目录

/

返回文章
返回