Processing math: 100%

装药缺陷对浇注装药快烤响应特性的影响

梁明阳 智小琦 于永利 肖游

李胜华, 李金良, 王旭之. 具有碳纤维材料环的超高压模具[J]. 高压物理学报, 2021, 35(1): 013301. doi: 10.11858/gywlxb.20200580
引用本文: 梁明阳, 智小琦, 于永利, 肖游. 装药缺陷对浇注装药快烤响应特性的影响[J]. 高压物理学报, 2025, 39(4): 045101. doi: 10.11858/gywlxb.20240893
LI Shenghua, LI Jinliang, WANG Xuzhi. Ultra-High Pressure Die with Carbon Fiber Ring[J]. Chinese Journal of High Pressure Physics, 2021, 35(1): 013301. doi: 10.11858/gywlxb.20200580
Citation: LIANG Mingyang, ZHI Xiaoqi, YU Yongli, XIAO You. Effect of Charge Defects on the Fast Cook-off Response Characteristics of Cast PBX Explosive Charge[J]. Chinese Journal of High Pressure Physics, 2025, 39(4): 045101. doi: 10.11858/gywlxb.20240893

装药缺陷对浇注装药快烤响应特性的影响

doi: 10.11858/gywlxb.20240893
详细信息
    作者简介:

    梁明阳(1998-),男,硕士研究生,主要从事不敏感弹药研究. E-mail:1522775234@qq.com

    通讯作者:

    智小琦(1963-),女,博士,教授,主要从事战斗部毁伤技术及不敏感弹药研究. E-mail:zxq4060@sina.com

  • 中图分类号: O521.9; TJ55

Effect of Charge Defects on the Fast Cook-off Response Characteristics of Cast PBX Explosive Charge

  • 摘要: 为研究炸药装药缺陷对快速烤燃响应特性的影响,对Ⅰ型烤燃弹(无缺陷装药)和Ⅱ型烤燃弹(含缺陷装药)进行了快速烤燃试验。Ⅱ型烤燃弹的响应时间(128 s)短于Ⅰ型烤燃弹(132 s),且在5 m处的最大冲击波超压峰值(62.7 kPa)高于Ⅰ型烤燃弹(12.5 kPa)。试验结果表明,点火后Ⅱ型烤燃弹的响应较无缺陷的Ⅰ型烤燃弹更为剧烈,但二者的响应等级仍为燃烧反应。在此基础上,通过Fluent软件建立了池火与烤燃试件相互耦合的计算模型,模拟了试件在火焰中的受热情况。研究发现,缺陷越靠近装药表面,缺陷处的局部温度越高,但对装药的响应时间无明显影响。

     

  • 在当今世界,知识密集、技术密集的高新技术产业快速发展,涉及新材料技术、信息技术、机电一体化技术、航空航天技术等重点领域。超硬材料(如金刚石、立方氮化硼等)因具有摩擦系数低、天然状态的硬度高、耐磨性好、导热系数高等一系列优异的性能而成为新材料队伍中的一支新军[1]。超高压设备是人工合成超硬材料必不可少的基础设备。世界上用于生产超硬材料的压机主要有两种类型:两面顶压机和铰链式六面顶压机,如图1[2]所示。六面顶压机在我国应用广泛,为超硬材料产业发展作出了重要贡献,然而其极限腔体压力较低,难以合成高品质金刚石。这也导致我国人造金刚石虽然产量占全球总产量的90%[3],但是品质却处于中、低档水平。鉴于此,我国对两面顶超高压模具的研究一直持续进行。

    图  1  两种超高压装置
    Figure  1.  Two ultra-high pressure devices

    近年来,两面顶超高压模具的发展方向主要有两个:一是改进内部压砧结构,二是改变外部预紧方式。对于后者,大多采用多层钢环、钢丝缠绕、钢带缠绕3种超高压模具预紧形式,然而它们却存在钢环直径过大、钢丝和钢带容易高温蠕变及松弛等问题[4]。碳纤维具有高比强度、高比模量等优异的材料性能,最初被用于制造大型飞机的结构材料、缠绕压力容器等。从碳纤维复合材料的材料性能上看,将其运用到模具预紧以改进预紧方式将会是一种不错的选择。在这方面,文献[5]中虽有分析,但未对进一步应用进行深入探讨,为此本研究拟对具有碳纤维复合材料预紧的两面顶超高压模具进行研究,以期为其设计和制造提供一定的科学参考。

    本研究设计的具有复合材料环的超高压模具从结构上看属于年轮式两面顶超高压装置。压缸是模具的主要承压部件,在工作状态下所受的周向拉应力是压缸破坏的主要原因。为避免压缸处于周向拉应力环境,基于“未裂先分”的思想,对压缸进行切分[6],在周向范围内分成离散的多块。这种方式能够最大限度地减小压缸周向拉应力,降低压缸失效的危险程度。为提高模具的承压极限,使用可承受高于钢环周向抗拉强度的碳纤维复合材料环代替模具最外层环。在模具大型化过程中,压缸体积增大带来的烧结困难可以通过对压缸剖分来解决,而支撑环的大直径、大厚度问题则可以通过使用复合材料环来解决。在实际生产中,通常在最外侧支撑环外增加一层钢环,以保护模具,同时避免因断裂破坏而对操作人员造成伤害。而以缠绕方式制作的复合材料环的破坏是由内而外循序渐进的,即使最内层纤维层部分断裂,压力也会分配给其余的纤维层,直到整个纤维层破坏,因此可不再设置保护环。图2显示了具有碳纤维复合材料环的超高压模具,其中图2(a)为模具结构示意图,图2(b)为模具施加材质之后的渲染图。

    图  2  具有碳纤维复合材料环的超高压模具
    Figure  2.  Ultra-high pressure mold with carbon fiber composite ring

    本研究中具有碳纤维复合材料环的超高压模具是在全钢环模具的基础上设计的。先从压缸的6块剖分式全钢环模具结构入手,内腔直径设计为40 mm,取压缸内外径比为1∶4.5,压腔的径高比为1∶1,模具整体高径比为2∶1,锥角为40°,则压缸的外径为180 mm,压腔高度h为40 mm,模具的总高度H为80 mm。

    图3给出了模具结构参数。以压腔压力为6 GPa作为工作条件,根据倍压原理可以得到最内层支撑环受到的内压约为1333 MPa。根据多层厚壁圆筒的最佳设计原则,可计算出支撑环所需的层数和各层尺寸,即

    图  3  超高压模具的几何尺寸
    Figure  3.  Geometry of ultra-high pressure mold
    pint=n[τ](ω2/n1)ω2/n
    (1)
    δi=2RiE(pipi1)(k2i+1k2i+11+1k2i1)
    (2)

    式中:pint为最大工作内压,MPa;n为支撑环层数;ω为模具的总径比;δi为第i层支撑环的过盈量,mm;Ri为第i层支撑环的内径,mm;E为支撑环的弹性模量,MPa;pi为第i层支撑环所受的合成压力,MPa;ki为第i层支撑环的径比。根据式(1),支撑环层数取5,各层径比均为1.224,相应的各层支撑环半径分别为110.00、134.80、165.00、202.00和247.00 mm。最外层环所受接触压力为266.67 MPa。由式(2)计算得到支撑环相邻过盈量由内而外分别为0.400、0.279、0.340、0.419和0.513 mm。

    本研究将一层复合材料环应用于全钢环模具上,替代在预紧状态下提供266.67 MPa预紧压力的、厚度为45 mm、内径为202 mm的最外层钢环。复合材料环采用与全钢环等厚度设计,即复合材料环的厚度为45 mm,以与全钢环模具压缸获得相同预紧压力为目标,复合材料环与内层钢环的过盈量为参数,基于Workbench的参数化设计模块,得到最佳过盈量为0.734 mm。

    超高压模具中压缸的材料选择YG8硬质合金,外层钢环选择超高强度钢材料45CrNiMoVA,主要性能参数列于表1[7],其中:ρ0为密度,μ为泊松比,σc为剪切强度。

    表  1  压缸和支撑环的材料参数[7]
    Table  1.  Material parameters of cylinder and support ring[7]
    Materialρ0/(g·cm−3)E/GPaμFailure strength/GPaσc/GPa
    YG814.605780.226.23.25
    45CrNiMoVA7.832100.291.60.80
    下载: 导出CSV 
    | 显示表格

    碳纤维复合材料环采用以T1000碳纤维为增强体、双马来酰亚胺树脂为基体制作而成的T1000/BM1复合材料,其性能参数如表2[8]所示,其中:G为剪切模量,ft为抗拉强度,fc为抗压强度,下标xyz表示xyz方向,下标xyyzxz表示xyyzxz面。

    表  2  T1000/BM1复合材料参数[8]
    Table  2.  Parameters of T1000/BM1 composite material[8]
    Ex/GPaEy/GPaEz/GPaμxyμyzμxz
    231.413.213.20.280.250.28
    Gxy/GPaGyz/GPaGxz/GPaft,x/MPaft,y/MPaft,z/MPa
    9.175.009.171731.066.566.5
    fc,x/MPafc,y/MPafc,z/MPaσc,xy/MPaσc,yz/MPaσc,xz/MPa
    −1379.0−268.2−268.2133.8100.0133.8
    下载: 导出CSV 
    | 显示表格

    考虑到超高压模具具有对称性,为节约计算资源,提高仿真实验效率,仅建立模具的1/12模型进行分析。各接触对之间均选用摩擦接触和增广拉格朗日算法,其中:压缸的各剖分块之间的摩擦因数为0.10,压缸与钢环之间、各层钢环之间、钢环与复合材料环之间的摩擦因数均为0.15。模具实际工作时,内壁可简化为均布压力,压缸锥面上的正压力p(s)与压缸内壁的均布压力存在一定的数学关系,其表达式[7, 9]

    p(s)=p0exp(τst)
    (3)

    式中:p0为压缸内壁受到的均布压力,MPa;τ为传压介质的内摩擦系数;t为锥面上密封介质的厚度,mm;s为传压介质到压腔端部边缘的距离,mm。复合材料环实体模型如图4所示,模具1/12模型的边界条件如图5所示。

    图  4  复合材料环实体模型
    Figure  4.  Model of composite ring
    图  5  模型的边界条件
    Figure  5.  Boundary condition of model

    压缸是模具中最重要的结构,其材料特性表现为承受压应力的能力较强而承受拉应力的能力较弱。图6显示了压缸预紧状态和工作状态时的周向应力分布。预紧状态下,压缸周向应力范围为−1357.65~−3065.67 MPa,表现为压应力,且随着半径的增大而逐渐减小,内壁处压应力最大。工作状态下,压缸受到工作载荷和预紧力的共同作用,压缸周向应力范围为−818.80~519.49 MPa,靠近压缸内壁部分因预紧力不足以完全抵消加压产生的周向拉应力,从而表现为拉应力,但是这部分拉应力在材料的可承受范围之内;此外,压缸块内壁中间位置的拉应力最大,两边位置因与相邻压缸块处于摩擦状态,应力状态大幅改善,可见剖分之后压缸周向拉应力的改善效果明显。

    图  6  压缸的周向应力分布
    Figure  6.  Circumferential stress distribution of pressure cylinder

    图7显示了压缸处于预紧状态和工作状态时的最大剪应力分布。预紧状态下,压缸最大剪应力的范围为703.69~1534.30 MPa,整体上看,最大剪应力随半径的增大而逐渐减小,最大值位于压缸内壁与锥面相交处,最小值位于外壁。工作状态下,压缸最大剪应力为194.89~3253.63 MPa,由内而外呈逐渐减小的分布规律,最大值位于压缸内壁且几乎达到压缸的抗剪强度,说明压缸的最危险位置处于压缸内壁面。

    图  7  压缸的最大剪应力分布
    Figure  7.  Maximum shear stress distribution of pressure cylinder

    图8显示了压缸在预紧状态和工作状态下的等效应力分布。预紧状态下,压缸整体受到的等效应力范围为1311.58~2996.61 MPa;工作状态下,等效应力范围为342.60~ 6057.87 MPa。从应力云图上看,无论是预紧状态还是工作状态,等效应力分布情况与最大剪应力分布情况十分相似,均随着半径的增大而逐渐减小。压缸内壁的等效应力最大,但是压缸块之间接触面处的等效应力相对内壁面有所降低。

    图  8  压缸的等效应力分布
    Figure  8.  Equivalent stress distribution of pressure cylinder

    综合以上分析,压缸的剖分对于减小压缸周向拉应力和压缸内壁等效应力都是有益的。另外,在工作状态下,压缸周向应力分布很不均匀,靠近内壁处为拉应力,靠近外壁处为压应力,是引起压缸破坏失效的主要原因。压缸最大剪应力和等效应力由内而外的变化幅度较大,外层材料没有得到充分利用。

    为了充分发挥复合材料层的性能,在超高压模具的外层增加一层保护钢环,其1/4模型如图9所示,对加保护环的模具进行应力分析。虽然外加保护环时模具的应力基本没有变化,但是却可以提高模具的安全系数,避免模具破坏造成不可预料的危险。

    图  9  外加保护层的超高压模具1/4模型
    Figure  9.  1/4 model of ultra-high pressure mold with protective layer

    除了对压缸进行应力分析外,还需要对最外层的复合材料环进行失效分析。复合材料的失效判别准则选用应用较广泛的蔡-吴准则。逆储备系数c是判断复合材料是否安全的常数。当c > 1时,可认为材料失效破坏;当c < 1时,材料安全,且越接近1,安全裕度越小。本研究以蔡-吴准则为计算依据,提取整个复合材料环的逆储备系数。图10显示了复合材料环处于预紧状态和工作状态时的逆储备系数分布。可见,逆储备系数由内而外逐渐递减,最内层最先发生破坏。两种状态下,逆储备系数的最大值分别为0.949和0.987,说明即使在工作状态下,复合材料环的最危险层(靠近内壁面的第1层铺层)仍满足强度要求。

    图  10  复合材料环的逆储备系数分布
    Figure  10.  Distribution diagram of inverse reserve coefficient of composite ring

    本研究以等剪应力原则设计并计算了模具尺寸,其中第4层与第5层支撑环之间的设计接触压力为266.67 MPa。图11给出了工作状态下第4层和第5层支撑环的接触面压力。可见,接触面上边缘处的接触压力较大,最大值为346.81 MPa,最小值为223.27 MPa,平均接触压力为279.28 MPa,与设计压力接近,说明复合材料环能够提供设计所需的预紧压力。

    图  11  第4层与第5层支撑环的接触面压力
    Figure  11.  Contact surface pressure of the fourth and fifth layers of support rings

    以具有碳纤维复合材料环的超高压模具为研究对象,通过使用碳纤维复合材料环代替全钢环模具的最外层钢制支撑环,建立了一种以钢环加复合材料环共同对压缸实施预紧的预紧方式,避免了大直径钢环难以制造加工的问题。计算结果表明:与全钢环模具最外层支撑环相比,等厚度的复合材料环在模具预紧状态和工作状态下均满足结构强度要求,且应力分布更加合理;将压缸剖分后,能够大大减小压缸的周向拉应力、最大剪应力和等效应力。

  • 图  烤燃系统压力测点俯视图

    Figure  1.  Top view of pressure points in the cook-off system

    图  快速烤燃示意图

    Figure  2.  Schematic diagram of fast cook-off

    图  Ⅱ型烤燃弹预制缺陷CT图

    Figure  3.  CT images of the pre-fabricated defect of TypeⅡcook-off bomb

    图  温度-时间历史曲线

    Figure  4.  Temperature-time history curves

    图  两型弹响应时的状态

    Figure  5.  During-response condition of two types of cook-off bombs

    图  两型弹响应后的状态

    Figure  6.  Post-response condition of two types of cook-off bombs

    图  烤燃模拟试件模型

    Figure  7.  Simulation models of Cook-off bombs

    图  计算域模型

    Figure  8.  Computational domain model

    图  池火模拟场景

    Figure  9.  Pool fire simulation scenario

    图  10  池火截面温度云图

    Figure  10.  Temperature of the cross-section of pool fire

    图  11  Ⅰ型烤燃弹火焰中测点温度随时间变化曲线

    Figure  11.  Temperature-time curves of gauging points in flame of type Ⅰ cook-off bomb

    图  12  不同时刻药柱表面温度云图(Ⅰ型烤燃弹)

    Figure  12.  Temperature contour map of the charge surface at different times (Type Ⅰ cook-off bomb)

    图  13  点火时刻药柱表面温度云图(Ⅰ型烤燃弹)

    Figure  13.  Temperature contour map of the propellant surface at the moment of ignition (Type Ⅰ cook-off bomb)

    图  14  Ⅱ型烤燃弹火焰中测点温度随时间变化曲线

    Figure  14.  Temperature-time curves of gauging points in flame of type Ⅱ cook-off bomb

    图  15  点火时刻药柱表面温度云图(Ⅱ型烤燃弹)

    Figure  15.  Temperature contour map of the propellant surface at the moment of ignition (Type Ⅱ cook-off bomb)

    图  16  点火前缺陷和药柱截面温度云图

    Figure  16.  Temperature contour map at the defect and propellant cross-section before ignition

    图  17  区域A缺陷位置示意图

    Figure  17.  Defect location diagram in area A

    图  18  区域B缺陷位置示意图

    Figure  18.  Defect location diagram in area B

    图  19  响应前不同位置缺陷及药柱截面温度云图

    Figure  19.  Temperature contour maps of defects at different positions and charge cross-sections before ignition

    图  20  缺陷与药柱的最高温度随时间变化曲线

    Figure  20.  Time-dependent curves of the maximum temperature of defects and explosive charges

    表  1  冲击波超压峰值

    Table  1.   Peak overpressure of shock wave

    Cook-off
    bomb
    Monitor point
    location
    Distance/m Overpressure
    peak/kPa
    Cook-off
    bomb
    Monitor point
    location
    Distance/m Overpressure
    peak/kPa
    TypeⅠ Axial 5 5.9 TypeⅡ Axial 5 47.5
    8 5.4 8 15.9
    10 4.8 10 11.8
    Radial 5 12.5 Radial 5 62.7
    8 9.4 8 32.4
    10 8.9 10 3.2
    下载: 导出CSV

    表  2  材料的物性参数

    Table  2.   Physical parameters of the material

    Material ρ/(kg·m−3) C/(J·kg−1·K−1) λ/(W·m−1·K−1)
    Steel 7 850 480 43
    Charge 1 809 1 012 0.49
    Air 1.225 1 004 0.024 2
    下载: 导出CSV

    表  3  反应动力学参数

    Table  3.   Reaction kinetic parameters of explosives

    E/(J·mol−1)A/s−1Q/(J·kg−1)R/(J·mol−1·K−1)
    1158202.030×10142.293×1068.314
    下载: 导出CSV
  • [1] 江鹏. 烤燃作用下PBX炸药的热分解及细观热损伤研究 [D]. 北京: 北京理工大学, 2016.

    JIANG P. Thermal decomposition and microstructure study of PBX explosive subjected to cook-off [D]. Beijing: Beijing Institute of Technology, 2016.
    [2] 黄开书, 江涛, 吴欣欣, 等. 提高浇注PBX炸药装药品质技术途径 [J]. 兵工自动化, 2023, 42(10): 60–62. doi: 10.7690/bgzdh.2023.10.014

    HUANG K S, JIANG T, WU X X, et al. Technical approach to improve charging quality of cast PBX explosive [J]. Ordnance Industry Automation, 2023, 42(10): 60–62. doi: 10.7690/bgzdh.2023.10.014
    [3] CUENOT B, SHUM-KIVAN F, BLANCHARD S. The thickened flame approach for non-premixed combustion: principles and implications for turbulent combustion modeling [J]. Combustion and Flame, 2022, 239: 111702. doi: 10.1016/j.combustflame.2021.111702
    [4] WEI R C, WANG B W, HE Q Z, et al. An experimental study on aviation kerosene pool fire flame and thermal radiation [J]. Journal of Physics: Conference Series, 2022, 2263(1): 012018. doi: 10.1088/1742-6596/2263/1/012018
    [5] 王帅, 智小琦, 贾秋琳, 等. 基于mass_flux法的炸药火烧试验与数值仿真 [J]. 兵器装备工程学报, 2020, 41(8): 1–6. doi: 10.11809/bqzbgcxb2020.08.001

    WANG S, ZHI X Q, JIA Q L, et al. Experiment and numerical simulation of explosive firing based on mass_flux method [J]. Journal of Ordnance Equipment Engineering, 2020, 41(8): 1–6. doi: 10.11809/bqzbgcxb2020.08.001
    [6] 肖游, 智小琦, 王琦, 等. 火焰特征量对快速烤燃的影响 [J]. 高压物理学报, 2022, 36(5): 055201. doi: 10.11858/gywlxb.20220557

    XIAO Y, ZHI X Q, WANG Q, et al. Influence of flame characteristics on fast cook-off [J]. Chinese Journal of High Pressure Physics, 2022, 36(5): 055201. doi: 10.11858/gywlxb.20220557
    [7] 王洪波, 王旗华, 卢永刚, 等. PBX炸药细观孔洞缺陷对其冲击点火特性的影响 [J]. 火炸药学报, 2015, 38(5): 31–36. doi: 10.14077/j.issn.1007-7812.2015.05.006

    WANG H B, WANG Q H, LU Y G, et al. Effect of Meso-defect holes on the shock-to-ignition characteristics of PBX explosives [J]. Chinese Journal of Explosives & Propellants, 2015, 38(5): 31–36. doi: 10.14077/j.issn.1007-7812.2015.05.006
    [8] 周婷婷, 楼建锋. 含孔洞炸药晶体HMX冲击响应的分子动力学模拟 [J]. 含能材料, 2024, 32(1): 65–75. doi: 10.11943/CJEM2023167

    ZHOU T T, LOU J F. Molecular dynamic studies on the shock responses of energetic crystal HMX with cylindrical voids [J]. Chinese Journal of Energetic Materials, 2024, 32(1): 65–75. doi: 10.11943/CJEM2023167
    [9] GROSS M L, HEDMAN T D, MEREDITH K V. Considerations for fast cook-off simulations [J]. Propellants, Explosives, Pyrotechnics, 2016, 41(6): 1036–1043. doi: 10.1002/prep.201500253
    [10] 王志富. 有装药缺陷的凝聚相炸药烤燃的数值模拟 [D]. 淮南: 安徽理工大学, 2020.

    WANG Z F. Numerical simulation of the burning of defective condensed phase explosives [D]. Huainan: Anhui University of Science and Technology, 2020.
    [11] PETERS N. Turbulent combustion [M]. Cambridge: Cambridge University Press, 2000.
    [12] Ansys Inc. Fluent A.: ANSYS fluent theory guide [Z]. USA, 2013, 15317: 724–746.
    [13] ATTAR A A, POURMAHDIAN M, ANVARIPOUR B. Experimental study and CFD simulation of pool fires [J]. International Journal of Computer Applications, 2013, 70(11): 9–15. doi: 10.5120/12004-5790
    [14] 任玉新, 陈海昕. 计算流体力学基础 [M]. 北京: 清华大学出版社, 2006.

    REN Y X, CHEN H X. Fundamentals of computational fluid dynamics [M]. Beijing: Tsinghua University Press, 2006.
    [15] 付淑青. 国产航空煤油多组分替代燃料的燃烧反应机理研究 [D]. 沈阳: 沈阳航空航天大学, 2020.

    FU S Q. Study on combustion reaction mechanism of domestic aviation kerosene multi-component surrogate fuels [D]. Shenyang: Shenyang Aerospace University, 2020.
    [16] 黄生洪, 徐胜利, 刘小勇. 煤油超燃冲压发动机两相流场数值研究(Ⅲ)煤油在超燃流场中的多步化学反应特征 [J]. 推进技术, 2005, 26(2): 101–105. doi: 10.3321/j.issn:1001-4055.2005.02.002

    HUANG S H, XU S L, LIU X Y. Combustion flow of kerosene-fueled scramjet with 3D cavity (Ⅲ) multi-step chemistry characteristics of kerosene [J]. Journal of Propulsion Technology, 2005, 26(2): 101–105. doi: 10.3321/j.issn:1001-4055.2005.02.002
    [17] KUNDU K, PENKO P, YANG S. Reduced reaction mechanisms for numerical calculations in combustion of hydrocarbon fuels [C]//36th AIAA Aerospace Sciences Meeting and Exhibit. Reno: AIAA, 1998.
    [18] 曾娇. 开放空间航空煤油池火燃烧数值模拟 [D]. 哈尔滨: 哈尔滨工程大学, 2016.

    ZENG J. Numerical simulation of aviation fuel pool fire in open air [D]. Harbin: Harbin Engineering University, 2016.
    [19] 傅献彩. 普通高等教育十五国家级规划教材: 物理化学(上) [M]. 北京: 高等教育出版社, 2005: 99–103.
  • 加载中
图(20) / 表(3)
计量
  • 文章访问数:  51
  • HTML全文浏览量:  22
  • PDF下载量:  14
出版历程
  • 收稿日期:  2024-09-24
  • 修回日期:  2024-10-22
  • 网络出版日期:  2025-04-01
  • 刊出日期:  2025-04-05

目录

/

返回文章
返回