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

爆炸加载下装配垫片对金属柱壳膨胀断裂影响研究

张世文 金山 陈艳 郭昭亮 但加坤 刘明涛 汤铁钢

张世文, 金山, 陈艳, 郭昭亮, 但加坤, 刘明涛, 汤铁钢. 爆炸加载下装配垫片对金属柱壳膨胀断裂影响研究[J]. 爆炸与冲击, 2022, 42(8): 083102. doi: 10.11883/bzycj-2021-0456
引用本文: 张世文, 金山, 陈艳, 郭昭亮, 但加坤, 刘明涛, 汤铁钢. 爆炸加载下装配垫片对金属柱壳膨胀断裂影响研究[J]. 爆炸与冲击, 2022, 42(8): 083102. doi: 10.11883/bzycj-2021-0456
ZHANG Shiwen, JIN Shan, CHEN Yan, GUO Zhaoliang, DAN Jiakun, LIU Mingtao, TANG Tiegang. Influence of a cushion on dynamic expansion and fracture of an explosively-driven metallic cylinder[J]. Explosion And Shock Waves, 2022, 42(8): 083102. doi: 10.11883/bzycj-2021-0456
Citation: ZHANG Shiwen, JIN Shan, CHEN Yan, GUO Zhaoliang, DAN Jiakun, LIU Mingtao, TANG Tiegang. Influence of a cushion on dynamic expansion and fracture of an explosively-driven metallic cylinder[J]. Explosion And Shock Waves, 2022, 42(8): 083102. doi: 10.11883/bzycj-2021-0456

爆炸加载下装配垫片对金属柱壳膨胀断裂影响研究

doi: 10.11883/bzycj-2021-0456
基金项目: 国家自然科学基金重点项目(11932018);国家自然科学基金面上项目(12072332)
详细信息
    作者简介:

    张世文(1971- ),男,博士,副研究员,zhangswxueshu@163.com

    通讯作者:

    金 山(1981- ),男,学士,高级工程师,jin9835@163.com

  • 中图分类号: O346.1

Influence of a cushion on dynamic expansion and fracture of an explosively-driven metallic cylinder

  • 摘要: 采用多普勒光纤探针测速技术(Doppler pins system,DPS,又称全光纤位移干涉测速技术)和高速摄影技术,研究装配垫片对金属柱壳膨胀断裂的影响,获得了有无垫片对应柱壳外表面位置的速度曲线和垫片对柱壳膨胀断裂影响明显的高速摄影图像。实验结果表明:与无垫片区域相比,垫片区域的柱壳外表面经历了先凸起后内凹的过程,导致垫片对应柱壳的径向运动位移发生反复错位,最终低于无垫片区域约0.34 mm,该位移差可能导致柱壳发生径向剪切断裂;实验结果还表明,在垫片与间隙交界处两侧(沿垫片方向约7.5°、沿间隙方向约9°)处各增加了一条裂纹,该断裂模式既不同于环向拉伸断裂,也不同于45°的剪切断裂,而是由垫片/间隙边界产生的两束稀疏应力波传到柱壳外表面引起的扰动影响所致,这个新的断裂模式与柱壳材料的动态力学性能密切相关。数值模拟结果表明,装配垫片对柱壳断裂机制影响不仅包含该处附加的质量效应,还应考虑炸药通过垫片后作用在柱壳上的冲击加载幅值变化、冲击加载时序与其他部位不同步的差异,以及垫片/间隙交界处引起的表面波传播对柱壳断裂模式的后续发展行为的影响。
  • 爆炸加载下金属柱壳的膨胀断裂已有了大量研究:第1类为炸药与金属柱壳最终破裂速度的关系[1-4],关注破片终速度和毁伤效应;第2类为金属柱壳的断裂模式(拉伸断裂、纯剪切断裂和拉剪混合型断裂等) [5-9],认为不同的断裂模式发展进程不同,进而影响到柱壳的破裂时间,关注炸药与柱壳材料的影响[8-10]、不同壁厚柱壳断裂模式差异[11]和柱壳发生单旋剪切断裂模式的诱导因素[12]等。另外,目前,因材料制备、加工和工程设计的不同,关注集中在柱壳微缺陷[13]、表面微加工痕迹[14]和装配状态等对柱壳断裂的影响规律。而对微缺陷、表面塑性层和装配薄垫片对柱壳断裂的微小影响,传统实验数据不易精确判断,因此对实验的方法和精度的要求更高。

    装配垫片是轴对称结构各部件对心的重要组成部分。在工程装配过程中,炸药与金属柱壳之间不可避免会产生装配垫片或空气的间隙,这显然影响炸药对金属柱壳爆炸驱动的同步性,造成金属柱壳膨胀断裂的差异。与垫片对应的柱壳到底是外凸还是内凹于其他部位,垫片对柱壳膨胀断裂究竟如何影响,需通过精细实验予以判断。李涛等[15]开展了在平面爆轰加载下垫片与间隙对单层平面飞片运动影响的实验,指出相比空气间隙,垫片导致金属飞片的速度起跳较晚、起跳速度较高、末速度较低。该实验采用炸药正向驱动飞层,与滑移爆轰驱动柱壳有一定的差别。

    本文中,开展炸药滑移加载下装配垫片的柱壳膨胀断裂实验,通过DPS测速阵列、高速摄影获得有无垫片时柱壳膨胀断裂的速度曲线、动态演化图像,研究垫片对柱壳膨胀断裂的影响,并分析影响柱壳断裂的主要因素。

    实验装置如图1(a)所示:柱壳材料为45钢,高120.02 mm,外径59.10 mm,内径51.08 mm,壁厚4.01 mm;主炸药为PBX柱状炸药,直径50.01 mm,高100.00 mm;传爆药直径20.02 mm,高4.98 mm,中心点滑移爆轰。采用厚0.496 mm、宽20.05 mm、长118.20 mm的两条透明胶带垫片,沿轴向对称贴在主炸药180°两侧。测点布局如图1(b)所示:DPS探头测点1~5对应有垫片部位,测点6~10对应无垫片部位,相邻测点夹角3°,测点1~10角度分别为13.5°~13.5°,距主炸药上端面43.00 mm。测点11~13分别距主炸药上端面23.00、53.00和63.02 mm,角度均为0°(即垫片边缘位置)。

    图  1  柱壳实验装置及测点布局
    Figure  1.  The cylindrical shell experimental device and its measuring point layout

    图2所示,高速摄影拍摄角度与测速位置对称,用炸药爆轰冲击氩气发光作为前照明光源,白纱布为背景屏增加对比度。转镜式分幅相机转速为1.2×105 r/min,图像的时间间隔为1.0 μs,可获得40幅图像,获取从表面裂纹萌生到产物严重泄漏的膨胀断裂全过程。

    图  2  速度测量和高速摄影的现场布局
    Figure  2.  Field layout of velocity measurement and high-speed photography
    2.1.1   激光测速整体趋势

    图3为测点速度曲线,在相同高度测点1~10的速度曲线基本重合,测点11~13的速度曲线反映了主炸药在滑移爆轰作用下不同柱壳位置外表面的膨胀过程。速度曲线特征明显,包含弹性前驱波ue和冲击波usw,以及柱壳在炸药驱动下的持续加速过程。不同位置的弹性前驱波幅值几乎相同,约70 m/s;冲击波幅值略有差异,约1 100 m/s。45钢密度ρ=7.85 g/cm3,声速c0=3.896 km/s,λ=1.611,据p=ρDu=ρ(c0+λu)u,可得45钢弹性极限为1.086 GPa(45钢泊松比为0.3,可得屈服应力为620.6 MPa),炸药作用于柱壳的冲击幅值为20.65 GPa。并且,尽管弹性前驱波ue基本不变,但冲击波usw变化趋势明显,测点11、5或6、12和13冲击波幅值分别为1 091.49、1 047.74、1 030.01和1 014.63 m/s,越靠近起爆端冲击波越强。这说明,通常意义的点起爆柱状炸药产生滑移爆轰时定常段并非严格意义的定常段,在精密实验中需特别注意。

    图  3  测点的速度曲线
    Figure  3.  Velocity curves of measuring points

    图4为点起爆炸药对不同位置柱壳驱动的示意图,越靠近起爆端,冲击波波阵面与柱壳夹角θ越小,冲击压力越大,柱壳外表面的第一速度幅值越大。当柱壳外表面离起爆端距离较远时,冲击波波阵面与柱壳角度变化不大,可近似为定常段,采用平面透镜起爆方式比点起爆方式可更好获得滑移爆轰定常段实验数据。这为后面的二维数值模拟中炸药采用同步起爆提供依据。

    图  4  点起爆柱状炸药驱动柱壳运动的示意图
    Figure  4.  The diagram of cylindrical shell movement driven by point-initiated cylindrical explosive
    2.1.2   垫片对速度曲线的影响

    测点1偏离正中−13.5°,对应位置偏离正中11.78 mm,约在垫片(宽20 mm)中部位置,其速度具有装配垫片对柱壳影响的典型特征。测点10为无垫片柱壳位置。因测点5~6为有无垫片交界面处,在前期有可能有微小差别,但随着时间的推移,两者将会趋于一致。

    为了更清楚地观察有无垫片对柱壳膨胀断裂影响,选择DPS阵列中测点1和10、测点5和6的速度曲线进行对比,如图5(a)所示。测点1和10相比,有垫片时第一峰值速度大(1054.84−1020.54)34.30 m/s,后期无垫片的速度高于有垫片部位的,在17.47 μs时速度差可达(1637.84−1594.08)43.76 m/s。对4条曲线积分获得位移,以测点1径向位移为基点,如图5(b)所示。测点10先低于测点1最大为0.05 mm,后在17.85 μs高于测点1为0.34 mm。两点之间初始弧长为13.9 mm,膨胀变形50%后弧长约20 mm,估算平均径向剪应变为1.7%。测点6与测点5相比,先低0.02 mm再高0.02 mm,径向错位变化量为0.04 mm,径向剪应变变化量为1.8%,两者均远高于45钢材料塑性剪应变,测点5和6之间的径向剪应变反向可能导致裂纹的萌生、发展。

    图  5  有无垫片时柱壳膨胀速度和位移差曲线
    Figure  5.  Expansion velocity and displacement difference curves of the cylindrical shell with or without a cushion

    以实验装置雷管起爆为零时,时标雷管设置为装置起爆后5 μs。由高速摄影图像可知,时标雷管在第7幅时开始发光,可推算传爆药起爆在2 μs时,即第2幅对应传爆药起爆。以此类推,第19幅对应传爆药起爆17 μs时的柱壳膨胀断裂图像。

    图6为17、19、21和23 μs时的柱壳膨胀图像。垫片处柱壳明显不同于其他位置,特别在19、21和23 μs时,明显有垫片位置的痕迹,局部有少量产物溢出。这与DPS阵列速度曲线结果相互印证。

    图  6  垫片对柱壳膨胀断裂的影响
    Figure  6.  Influence of cushion on cylindrical shell expansion fracture

    图7为25和27 μs时的柱壳膨胀图像,爆轰产物从垫片与间隙交界处少量溢出。由柱壳演化过程,因测点5和6相距较近,速度差异远不如测点1和10的差异大。这是由于,在膨胀过程中,45钢柱壳环向有拉伸应力,该拉伸应力抑制了垫片交界处引起的径向位移差,即径向具有一定的调节功能,速度较高的测点6和速度较低的测点5相互牵制而趋于一致。当测点1和10的速度差沿环向传播过来,测点5和6处无法承受径向方向的剪应力而发生剪切断裂。该断裂模式既不同于传统的环向拉伸断裂,也不同于45°角剪切断裂。由速度曲线和高速摄影图像可见,测点5和6、测点1和10处在冲击波7.5 μs时已有速度差,测点5和6处发生断裂所需时间约为8 μs。这与材料的屈服应力和延伸率有关,屈服应力越大,延伸率越大,垫片处发生径向剪切断裂越晚。

    图  7  垫片对柱壳膨胀断裂爆轰产物的影响
    Figure  7.  Influence of cushion on detonation products of cylindrical shell expansion fracture

    图67中尺标距推算,垫片处柱壳明显裂纹的破片宽度小于20 mm。为了更深入理解垫片对柱壳膨胀断裂的影响,进一步分析速度原始波形(见图8)。测点3、6和9的速度在后期出现分叉,其余测点未出现分叉情况,可推断在测点3、6和9处出现了裂纹。对应的每条破片环向张角为9°,如按原始柱壳尺寸(外径59 mm)计算宽度为4.63 mm,若按半径膨胀15 mm(由60 mm膨胀至90 mm)计算,破片宽度为7.07 mm,低于垫片宽度。

    图  8  激光测速谱图
    Figure  8.  Original spectra of laser velocity measurement

    因实验数据有限,再对上述实验进行数值模拟,研究垫片对柱壳的影响。一方面,确认冲击波经过垫片和间隙的时序及柱壳速度的差异;另一方面,分析垫片影响柱壳膨胀断裂的主要因素。

    直接采用三维计算模型,计算量大,且垫片、柱壳和间隙之间三维接触计算难度较大。因此,数值模拟采用两种简化二维计算模型,一种为轴对称计算模型,另一种为平面应变计算模型。两种计算模型结果均表明,垫片处的速度因附加质量导致最终速度小于无垫片处的。这里,重点关注冲击波通过垫片与间隙后对柱壳膨胀速度的影响细节。

    在轴对称模型中,将传爆药与主炸药合二为一,如图9所示。上端中心点起爆,炸药与柱壳之间一个有垫片,一个没有垫片,垫片(透明胶带)密度约1 g/cm3,垫片材料参数暂用有机玻璃的代替。分析离主炸药上端面43 mm处的金属柱壳内界面单元和外界面单元的压力曲线。由计算结果可知,爆轰冲击波通过垫片时压力明显衰减,最大压力比无垫片的冲击波压力小,但冲击波到达柱壳内表面早于无垫片。在空气隙扩散时炸药产生的爆轰产物冲击压力衰减斜率远大于有垫片的,到达柱壳外表面时冲击压力反而低于有垫片的。这使有垫片柱壳外表面的第1幅值大于无垫片的。柱壳外表面垫片和间隙的速度第1幅值和最终速度规律与实验测速一致。

    图  9  轴对称计算模型和柱壳内外表面的冲击压力曲线
    Figure  9.  Axisymmetric calculation model and impact pressure curves on the inner and outer surfaces of the cylindrical shell

    图10为平面应变模型,由点起爆柱状炸药驱动柱壳运动性质(见图4),采用所有炸药零时同步起爆方式,关注垫片与间隙交界处垫片、柱壳外表面运动规律。比较偏离垫片和间隙交界处0°、7°、9°和18°处柱壳外表面速度曲线,在垫片方向,越靠近交界处,速度第1幅值越小。这说明,柱壳外表面速度受到交界处影响,并向垫片传播环向表面波。

    图  10  平面应变计算模型、垫片变形和柱壳外表面速度曲线
    Figure  10.  Plane strain calculation model, cushion deformation and velocity curves of cylindrical shell outer surface

    图11为垫片/间隙对应柱壳外表面速度曲线和环向表面波传播示意图,垫片中间区域对应的柱壳外表面测点1速度起跳幅值高于边界处测点5的,末速度低于测点5的。同样,柱壳外表面测点10速度起跳幅值低于测点6的,末速度高于测点6的。结合数值模拟和实验结果,当垫片足够宽或测点离交界处较远时,柱壳外表面第1速度幅值基本一样,不受交界处影响,计算中9°和18°处(对应单元为5911和5921)第1幅值差别很小,实验中测点1~3第1速度幅值无法区分大小,同样的规律呈现在间隙方向。垫片/间隙交界处对柱壳膨胀断裂的影响引入两个扰动,交界处质量突变,可能造成实验中测点6速度分叉;交界处发射的两束应力波沿着垫片和间隙方向传播至柱壳外表面,可能造成测点3和9速度分叉。应力波传播角度和影响范围由垫片(包括垫片厚度、垫片横向移动,见图10中垫片变形)和柱壳(包括柱壳材料声速、柱壳厚度)决定。通常垫片与柱壳的质量比越大,垫片屈服应力越大,交界处垫片边界横向滑动越小,交界处质量梯度差异越大,越可能在交界处首先断裂。柱壳厚度越小,在偏离交界处引起的断裂破片宽度越小。

    图  11  垫片/间隙对应柱壳外表面速度曲线和环向表面波传播示意图
    Figure  11.  Velocity curves of the outer surface of cylindrical shell and cushion/gap and propagation of annular surface waves

    开展了装配垫片对柱壳膨胀断裂影响实验。根据DPS阵列速度曲线和高速摄影图像,获得了如下结论。

    (1) 炸药与金属柱壳的垫片改变了施加在金属柱壳的冲击压力的时序、冲击幅值等。与无垫片相比,有垫片处柱壳外表面速度启动快,速度第1峰值大,在0.5 mm透明胶垫片作用下,启动时间差可达0.16 μs,速度差可达43 m/s。

    (2) 装配垫片严重影响柱壳膨胀断裂过程,由于附加质量效应,装配垫片处的末速度低于柱壳其余部位的,结合垫片启动不一致,使垫片柱壳部位在膨胀过程中先外凸、后内凹于无垫片对应位置,交界处对应位置径向方向发生二次错位,容易产生径向剪切断裂。

    (3) 垫片边界的横向移动和环向表面波具有抑制径向剪切断裂的趋势。

    (4) 在垫片/间隙交界处的两边(沿垫片方向约7.5°和沿间隙方向9°)裂纹来自于交界处发射的两束稀疏波和环向表面波的传播影响。

  • 图  1  柱壳实验装置及测点布局

    Figure  1.  The cylindrical shell experimental device and its measuring point layout

    图  2  速度测量和高速摄影的现场布局

    Figure  2.  Field layout of velocity measurement and high-speed photography

    图  3  测点的速度曲线

    Figure  3.  Velocity curves of measuring points

    图  4  点起爆柱状炸药驱动柱壳运动的示意图

    Figure  4.  The diagram of cylindrical shell movement driven by point-initiated cylindrical explosive

    图  5  有无垫片时柱壳膨胀速度和位移差曲线

    Figure  5.  Expansion velocity and displacement difference curves of the cylindrical shell with or without a cushion

    图  6  垫片对柱壳膨胀断裂的影响

    Figure  6.  Influence of cushion on cylindrical shell expansion fracture

    图  7  垫片对柱壳膨胀断裂爆轰产物的影响

    Figure  7.  Influence of cushion on detonation products of cylindrical shell expansion fracture

    图  8  激光测速谱图

    Figure  8.  Original spectra of laser velocity measurement

    图  9  轴对称计算模型和柱壳内外表面的冲击压力曲线

    Figure  9.  Axisymmetric calculation model and impact pressure curves on the inner and outer surfaces of the cylindrical shell

    图  10  平面应变计算模型、垫片变形和柱壳外表面速度曲线

    Figure  10.  Plane strain calculation model, cushion deformation and velocity curves of cylindrical shell outer surface

    图  11  垫片/间隙对应柱壳外表面速度曲线和环向表面波传播示意图

    Figure  11.  Velocity curves of the outer surface of cylindrical shell and cushion/gap and propagation of annular surface waves

  • [1] MOTT N F. Fragmentation of rings and shells: a theory of the fragmentation of shells and bombs [M]. Berlin, Germany: Springer Press, 2006: 243−294.
    [2] MOTT N F. Fragmentation of shell cases [J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1947, 189(1018): 300–308. DOI: 10.1098/rspa.1947.0042.
    [3] GRADY D. Investigation of explosively driven fragmentation of metals: two-dimensional fracture and fragmentation of metal shells [R]. Livermore, USA: Lawrence Livermore National Laboratory, 2003.
    [4] GURNEY R W. The initial velocities of fragments from bombs, shells and grenades [R]. Aberdeen, UK: Ballistic Research Laboratory, 1943.
    [5] HIROE T, FUJIWARA K, HATA H, et al. Explosively driven expansion and fragmentation behavior for cylinders, spheres and rings of 304 stainless steel [J]. Materials Science Forum, 2010, 638: 1035–1040.
    [6] HIROE T, FUJIWARA K, HATA H, et al. Deformation and fragmentation behaviour of exploded metal cylinders and the effects of wall materials, configuration, explosive energy and initiated locations [J]. International Journal of Impact Engineering, 2008, 35(12): 1578–1586. DOI: 10.1016/j.ijimpeng.2008.07.002.
    [7] 禹富有, 董新龙, 俞鑫炉, 等. 不同填塞装药下金属柱壳断裂特性的实验研究 [J]. 兵工学报, 2019, 40(7): 1418–1424. DOI: 10.3969/j.issn.1000-1093.2019.07.011.

    YU F Y, DONG X L, YU X L, et al. Fracture characteristics of metal cylinder shells with different charges [J]. Acta Armamentarii, 2019, 40(7): 1418–1424. DOI: 10.3969/j.issn.1000-1093.2019.07.011.
    [8] 汤铁钢, 李庆忠, 孙学林, 等. 45钢柱壳膨胀断裂的应变率效应 [J]. 爆炸与冲击, 2006, 26(2): 129–133. DOI: 10.11883/1001-1455(2006)02-0129-05.

    TANG T G, LI Q Z, SUN X L, et al. Strain-rate effects of expanding fracture of 45 steel cylinder shells driven by detonation [J]. Explosion and Shock Waves, 2006, 26(2): 129–133. DOI: 10.11883/1001-1455(2006)02-0129-05.
    [9] 胡八一, 董庆东, 韩长生, 等. TC4钛合金自然破片的引燃机理分析 [J]. 爆炸与冲击, 1995, 15(3): 254–258.

    HU B Y, DONG Q D, HAN C S, et al. Analysis of the firing mechanics for TI-6AL-4V natural fragments [J]. Explosion and Shock Waves, 1995, 15(3): 254–258.
    [10] 沈飞, 王辉, 屈可朋, 等. 不同晶粒度无氧铜管在爆轰加载下的膨胀及断裂特性 [J]. 爆炸与冲击, 2020, 40(2): 022201. DOI: 10.11883/bzycj-2019-0063.

    SHEN F, WANG H, QU K P, et al. Expansion and fracture characteristics of oxygen-free copper tubes with different grain sizes under detonation loading [J]. Explosion and Shock Waves, 2020, 40(2): 022201. DOI: 10.11883/bzycj-2019-0063.
    [11] 汤铁钢, 谷岩, 李庆忠, 等. 爆轰加载下金属柱壳膨胀破裂过程研究 [J]. 爆炸与冲击, 2003, 23(6): 529–533.

    TANG T G, GU Y, LI Q Z, et al. Expanding fracture of steel cylinder shell by detonation driving [J]. Explosion and Shock Waves, 2003, 23(6): 529–533.
    [12] 胡海波, 汤铁钢, 胡八一, 等. 金属柱壳在爆炸加载断裂中的单旋现象 [J]. 爆炸与冲击, 2004, 24(2): 97–107.

    HU H B, TANG T G, HU B Y, et al. A study of uniform shear bands orientation selection tendency on explosively loaded cylindrical shells [J]. Explosion and Shock Waves, 2004, 24(2): 97–107.
    [13] 金山, 张世文, 龙建华. 缺陷对圆管膨胀断裂影响的实验研究 [J]. 高压物理学报, 2011, 25(2): 188–192. DOI: 10.11858/gywlxb.2011.02.017.

    JIN S, ZHANG S W, LONG J H. Experimental study on the influences of defects on expanding fracture of a metal cylinder [J]. Chinese Journal of High Pressure Physics, 2011, 25(2): 188–192. DOI: 10.11858/gywlxb.2011.02.017.
    [14] LIU M T, GUO Z L, FAN C, et al. Modeling spontaneous shear bands evolution in thick-walled cylinders subjected to external high-strain-rate loading [J]. International Journal of Solids and Structures, 2016, 97: 336–354. DOI: 10.1016/j.ijsolstr.2016.07.014.
    [15] 李涛, 刘明涛, 王晓燕, 等. 装配垫层与间隙对爆轰加载下金属飞片运动特征的影响 [J]. 高压物理学报, 2018, 32(4): 044202. DOI: 10.11858/gywlxb.20170576.

    LI T, LIU M T, WANG X Y, et al. Effects of explosive device with foam cushion and air clearance on kinetic characteristic of steel flyer under detonation loading [J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 044202. DOI: 10.11858/gywlxb.20170576.
  • 期刊类型引用(0)

    其他类型引用(1)

  • 加载中
图(11)
计量
  • 文章访问数:  540
  • HTML全文浏览量:  171
  • PDF下载量:  92
  • 被引次数: 1
出版历程
  • 收稿日期:  2021-11-04
  • 修回日期:  2022-03-15
  • 网络出版日期:  2022-03-29
  • 刊出日期:  2022-09-09

目录

/

返回文章
返回