冲击荷载对无烟煤微观孔隙影响的分形研究

王以贤 梁为民

孙远翔, 秦嘉鸿, 陈岩武, 王成. 基于PIV实验研究带破口双层圆柱结构附近气泡的动力学特性[J]. 高压物理学报, 2024, 38(6): 064104. doi: 10.11858/gywlxb.20240756
引用本文: 王以贤, 梁为民. 冲击荷载对无烟煤微观孔隙影响的分形研究[J]. 高压物理学报, 2020, 34(5): 054203. doi: 10.11858/gywlxb.20200528
SUN Yuanxiang, QIN Jiahong, CHEN Yanwu, WANG Cheng. Bubble Dynamics Characteristics Near Double-Layer Cylindrical Structures with a Hole Based on PIV Experiments[J]. Chinese Journal of High Pressure Physics, 2024, 38(6): 064104. doi: 10.11858/gywlxb.20240756
Citation: WANG Yixian, LIANG Weimin. Fractal Study on Influence of Impact Load on Microscopic Pore of Anthracite[J]. Chinese Journal of High Pressure Physics, 2020, 34(5): 054203. doi: 10.11858/gywlxb.20200528

冲击荷载对无烟煤微观孔隙影响的分形研究

doi: 10.11858/gywlxb.20200528
基金项目: 国家自然科学基金(41772163)
详细信息
    作者简介:

    王以贤(1981–),男,博士研究生,讲师,主要从事爆破工程及岩土工程研究.E-mail:56653648@qq.com

    通讯作者:

    梁为民(1967–),男,博士,教授,主要从事爆破工程及岩土工程研究.E-mail:liangwm@hpu.edu.cn

  • 中图分类号: O347.2

Fractal Study on Influence of Impact Load on Microscopic Pore of Anthracite

  • 摘要: 为了揭示冲击荷载对无烟煤微观孔隙结构的影响规律,利用分离式霍普金森压杆(SHPB)冲击实验系统模拟了冲击应力在不同衰减过程中的冲击波和应力波,基于冲击前后压汞实验及低温液氮实验测试数据,运用分形理论,研究了赵固二矿不同方向煤体(与层理方向分别呈垂直、平行、45°斜交)冲击前后孔隙结构的分形特征。结果表明:对于渗流孔,冲击荷载提高了瓦斯的渗流与运移速度,对于吸附孔,冲击荷载减小了吸附孔的吸附能力,促进了瓦斯的解吸;分形维数具有明显的冲击方向性,且吸附孔的分形维数明显小于渗流孔;不同方向无烟煤的最佳冲击荷载不同,垂直和斜交层理方向的最佳冲击荷载为51.80 MPa,平行层理方向的最佳冲击荷载为28.46 MPa。研究结果可为冲击荷载促进瓦斯抽采机理的探讨提供参考。

     

  • 目前,潜艇是世界各国使用最多的大型水下航行体,潜艇具有隐蔽性和机动性,其制海范围比水面舰船要大,同时潜艇还具备较好的攻击性能,使其成为海战中的重要组成部分。潜艇的典型外观结构形式为圆柱壳形,我国潜艇借鉴了俄罗斯风格,采用双壳体结构。双壳体潜艇与单壳体潜艇的最大区别在于压载水舱的布置方式,双壳体潜艇的压载水舱包裹在潜艇耐压结构外侧,形成了非耐压壳包裹耐压壳体的双层结构。对于此种潜艇的双层结构,当冲击波载荷将外侧的非耐压壳体破坏后,内部的耐压壳体会受到气泡脉动载荷和射流载荷的冲击。王玉等[1]研究了3种聚能战斗部对典型潜艇目标的毁伤效应,发现3种战斗部的毁伤均表现为潜艇壳体小尺寸穿孔和大面积凹陷变形。贺铭等[2]对近场水下爆炸气泡与双层破口结构的相互作用机理进行了研究,发现水下爆炸冲击波首先使外层板出现破口,随后气泡通过破口进入双层板间舱室,进而对内底板造成毁伤。

    气泡载荷对双壳体潜艇耐压壳的作用过程本质上是近场水下爆炸气泡与双层圆柱破口结构的相互作用问题。目前,对于水下爆炸冲击波作用于外板毁伤效应的研究较多[37],但往往忽略了后续气泡与破口边界的耦合问题。气泡与具有初始破口的结构的相互作用问题较为复杂,需要考虑气泡与不完整边界的耦合效应,探索气泡对结构的毁伤规律[8]。Quah等[9]利用高速相机同时捕捉气泡动力学的侧视图像和前视图像,发现气泡形成射流的方向取决于双层板间的无量纲间隙宽度。Chen等[10]建立了一个带有圆孔的双层板简化模型,探究圆孔直径、两板间距以及爆距对气泡运动的影响。陈岩武等[11]进行了电火花气泡与双层结构模型实验,基于实验结果建立了实尺度的舰船双层底舱段模型,模拟结果表明,水下炸药在距离舰船底部结构2.5倍装药半径爆炸时,其毁伤威力最大。

    本研究将基于双层潜艇结构建立带破口的双层圆柱模型,利用电火花装置产生爆炸脉动气泡,观测气泡与破口之间的相互作用,利用粒子图像测速(particle image velocimetry,PIV)技术测量气泡周围流场速度以及气泡溃灭后产生的水射流速度,并测量圆柱内壁面处的压力载荷,探究爆距和破口直径对气泡脉动过程的影响。

    依据典型双壳体潜艇结构设计的双层圆柱模型如图1所示,外壳外径为70 mm,内壳内径为40 mm,外壳和内壳均为5 mm厚的透明亚克力板制成,以保证高速摄像机能够清晰捕捉到内外壳间的气泡脉动图像。不同破口直径的圆柱模型尺寸见表1。外壳底部预制直径dW分别为10、20和30 mm的破口,探究不同破口直径对气泡脉动的影响;内壁面预制直径dN为6 mm的破口,用于放置PCB压力传感器,测量内壁面受到的压力载荷。

    图  1  双层圆柱模型(单位:mm)
    Figure  1.  Double layered cylindrical model (Unit:mm)
    表  1  双层圆柱模型的尺寸以及参数
    Table  1.  Dimensions and parameters of the double layered cylindrical model
    Model No.Outer radius/mmInner radius/mmPlate thickness/mmdW/mmdN/mm
    170455106
    270455206
    370455306
    下载: 导出CSV 
    | 显示表格

    实验在长、宽、高均为60 cm的透明水箱内进行,水深约为50 cm,实验装置布置如图2所示。双层圆柱模型浸没于水面下方,实验中主要观察内外壳间的气泡脉动情况,内壳体环境对实验的影响较小,因此模型两端不封口,整体浸没于水中。搭接直径为0.3 mm的细铜丝制作电火花气泡发生装置,用于代替传统TNT炸药水下爆炸产生的气泡,定义铜丝交点到圆柱模型破口的距离为爆距DM,该装置产生的气泡内部为高压气体,避免了爆炸产物对气泡观测的影响,同时使得气泡脉动以及周围流场情况测量更为准确。实验中使用的电火花气泡发生装置产生的气泡最大直径不超过6 cm,水箱直径为60 cm,大于10倍气泡直径,因此水箱侧面和底部边界对气泡脉动的影响可忽略不计。

    图  2  实验装置
    Figure  2.  Experimental device

    为测量水射流速度,引入PIV技术对圆柱模型底部水下爆炸时的流场进行观测及测量。在待测流场中均匀布满粒径为10~60 μm的玻璃示踪粒子,玻璃示踪粒子为硅酸钠、二氧化硅和碳酸钠等烧结的混合物,密度为1.05 g/cm3,折射率为1.33%,适用于1~50 m/s流速范围内的流体速度测量实验。激光器产生的激光通过扩束器及透镜组后转变为厚度2~3 mm的均匀片光源覆盖待测流场,示踪粒子在激光照射后折射一部分激光进入CCD相机,通过高速摄像机连续采集时间间隔为Δt的目标场粒子位置图像,得到不同时刻的瞬态粒子场,测速系统如图3所示。

    图  3  PIV原理示意图
    Figure  3.  Schematic diagram of PIV method

    在PIV lab后处理软件中对比t时刻和t+Δt时刻2帧图像中示踪粒子的位置,利用2个瞬态粒子场的空间信息得出某一瞬时时刻t的速度向量场,重复此步骤即可得到完整气泡脉动周期的流场速度和水射流速度。PIV测量实验的设备参数:使用10 W连续激光器产生532 nm激光作为光源,利用高速摄像机(型号:Photron SA4)拍摄清晰的气泡运动画面,相机帧率为10000 幅/秒,分辨率为768×768,捕捉时间间隔为0.1 ms的清晰流场图像。对高速摄像机拍摄的画面进行标定,调整相机的位置和焦距,使拍摄画面尺寸与实际尺寸的比例为1∶1。

    1.4.1   自由场放电气泡脉动参照组实验

    为了确定自由场条件下气泡的动力学行为,首先进行了3组自由场放电气泡脉动实验作为参照组,控制电火花装置放电电压均为500 V。参照组实验中,气泡的最大半径分别为4.95、5.01和5.00 cm,记为工况1~工况3,气泡脉动周期均为4.8 ms。如图4所示,以工况1为例进行分析,放电时铜丝交点处短路燃烧放热,使得周围液体汽化,并在内部高压驱动下迅速膨胀形成电火花气泡;经过2.4 ms后气泡膨胀到最大,此时气泡的最大半径为4.95 cm;气泡达到最大体积后,在内外压力差的作用下开始收缩,4.8 ms时被压缩到最小体积,完成一个气泡脉动周期。由于细铜线的直径远小于气泡直径,因此其对气泡脉动过程的影响可以忽略不计。

    图  4  参照组实验(工况1)的气泡脉动过程
    Figure  4.  Bubble pulsation process of the reference group experiment (Case 1)
    1.4.2   压力载荷实验

    为了研究双层圆柱外壳结构对内壳壁面受到的压力载荷的影响,分别在自由场条件下和双层圆柱结构附近进行压力载荷实验。在自由场条件下,利用PCB自由场传感器,分别测量4、5、6、7、8 cm爆距处的压力数据,记为工况F-1~工况F-5。自由场中的爆距DF为铜丝交点(爆心)到传感器的水平距离。将传感器垂直放置于铜丝交点附近,连接电荷转换器和示波器,自由场传感器灵敏度为145 mV/MPa,将示波器所示电压换算为压力数值(压力单位为MPa)。

    对模型2进行爆距为1、2、3、4、5 cm的电火花压力载荷实验,记为工况M1~工况M5测量内壁面的冲击波峰值压力、气泡脉动压力和气泡射流载荷压力。为了使实验工作具有普遍意义,将爆距D和破口直径dW 2个变量进行无量纲化,定义爆距参数γ和破口参数λ分别为

    γ=D/db
    (1)
    λ=dW/db
    (2)

    式中:db为自由场条件下气泡脉动的最大直径,取参照组工况3的实验结果5.00 cm进行计算。工况M1~M5的γ分别为0.2、0.4、0.6、0.8和1.0。

    考虑到破口到内壁面的距离为3 cm,因此对比自由场与模型2的实验压力数据时,选用的模型2的工况爆距比自由场条件下的爆距小3 cm,例如,将模型2爆距为1 cm的实验工况与自由场条件下爆距为4 cm的实验工况进行对比。

    1.4.3   不同破口参数电火花爆炸实验

    为了研究不同破口参数对气泡脉动的影响,对3个模型分别进行爆距参数γ为0、0.1、0.2、0.3、0.4、0.6和0.8的水下电火花爆炸实验,实验工况如表2所示。

    表  2  模型1~模型3的实验工况
    Table  2.  Experimental cases for Model 1–Model 3
    Model Case λ γ
    1 M1-1 0.2 0
    M1-2 0.2 0.1
    M1-3 0.2 0.2
    M1-4 0.2 0.3
    M1-5 0.2 0.4
    M1-6 0.2 0.6
    M1-7 0.2 0.8
    2 M2-1 0.4 0
    M2-2 0.4 0.1
    M2-3 0.4 0.2
    M2-4 0.4 0.3
    M2-5 0.4 0.4
    M2-6 0.4 0.6
    M2-7 0.4 0.8
    3 M3-1 0.6 0
    M3-2 0.6 0.1
    M3-3 0.6 0.2
    M3-4 0.6 0.3
    M3-5 0.6 0.4
    M3-6 0.6 0.6
    M3-7 0.6 0.8
    下载: 导出CSV 
    | 显示表格

    为了验证PIV实验测量的准确性,通过LS-DYNA软件对自由场下的第3参照组进行模拟,水下爆炸气泡的最大半径

    Rm=3.56(WH+10.33)1/3
    (3)

    式中:W为装药质量,单位kg;H为炸药所处位置水深,单位m。由式(3)计算可知,在自由场条件下,0.0037 gTNT产生的水下爆炸气泡最大直径为5.01 cm,与实验中电火花装置产生的气泡最大直径十分接近,因此,选取该当量TNT作为爆炸源进行模拟。

    由于自由场实验工况具有对称性,为了提高计算效率,建立图5所示的1/4有限元模型,模拟水域大小为30 cm×30 cm×60 cm,炸药位于30 cm水深处,采用渐变网格对模型进行划分,爆心区域网格尺寸为0.45 mm,最远端网格尺寸为10 mm,炸药采用*MAT_HIGH_EXPLOSIVE_BURN材料模型和JWL爆轰产物状态方程描述,即

    图  5  有限元模型
    Figure  5.  Finite element model
    pe=A(1ωR1V)eR1V+B(1ωR2V)eR2V+ωEZV
    (4)

    式中:pe为爆轰产物的压力;V为爆轰产物的相对体积;EZ为炸药单位体积内能;ABR1R2ω为表征炸药特性的参数,取值见表3。表中:DCJ为CJ爆速,pCJ为CJ爆轰压力。

    表  3  炸药的材料模型及状态方程参数
    Table  3.  Material model and equation of state parameter of the explosive
    ρ0/(kgm3) DCJ/(ms1) pCJ/GPa EZ/GPa A/GPa B/GPa R1 R2 ω
    1630 6930 21.0 7.17 373.7 3.74 4.15 0.9 0.35
    下载: 导出CSV 
    | 显示表格

    水介质选用*MAT_NULL材料模型和*EOS_LINEAR_POLYNOMIAL多项式状态方程描述

    pw=C0+C1μ+C2μ2+C3μ3+(C4+C5μ+C6μ2)EC
    (5)

    式中:pw为压力,C0C6为常数,EC为水的单位体积内能。具体参数取值见表4

    表  4  水的材料模型及状态方程参数
    Table  4.  Material model and equation of state parameter of water
    ρ0/(kgm3) C0/GPa C1/GPa C2/GPa C3/GPa C4 C5 C6 EC
    1010 1.01×104 2.036 8.432 8.014 0.4934 1.3937 0 2.05×104
    下载: 导出CSV 
    | 显示表格

    由于自由场条件下不会产生气泡射流,因此对第1次脉动周期结束时的气泡收缩速度进行分析,第3参照组气泡的最大收缩速度如图6所示,选取气泡最大直径db、气泡脉动周期Tb和气泡最大收缩速度(vb)3个关键参数进行对比(表5),模拟计算结果与实验结果的最大误差为12.1%,说明本实验中采用的PIV测速系统的测量结果准确。

    图  6  实验工况3气泡最大收缩速度测量结果
    Figure  6.  Measurement results of maximum bubble contraction velocity of Case 3
    表  5  实验与数值模拟结果对比
    Table  5.  Comparison between experimental and simulation results
    db Tb vb
    Exp./cm Sim./cm Error/% Exp./ms Sim./ms Error/% Exp./(m·s−1) Sim./(m·s−1) Error/%
    5.00 4.99 0.2 4.80 4.42 8.5 43.118 49.069 12.1
    下载: 导出CSV 
    | 显示表格

    在自由场条件下,电火花水下爆炸气泡能量主要集中于第1个气泡脉动周期,第2、3个脉动周期中能量较少,且第2、3个脉动周期中的气泡直径与第1个脉动周期相比大幅减小,因此研究重点放在气泡的第1个脉动周期。自由场条件下压力测试实验结果如表6图7所示。第1个压力峰值为冲击波压力峰值,由铜丝通电短路燃烧瞬间产生,第2个压力峰值为气泡二次脉动压力峰值。气泡结束第1个脉动周期后,气泡体积缩小到最小,其内部压力高于周围环境的平衡压力,于是进行二次膨胀,产生气泡二次脉动压力。在自由场条件下,冲击波压力峰值和气泡二次脉动压力峰值均随着DF的增大而减小。

    表  6  自由场压力测量实验结果
    Table  6.  Experimental results of free field pressure measurement
    CaseDF/cmShock wave pressure/MPaPressure of bubble secondary wave/MPa
    F-140.5737.692
    F-250.4545.818
    F-360.3724.772
    F-470.2913.798
    F-580.2612.482
    下载: 导出CSV 
    | 显示表格
    图  7  自由场条件下不同工况的压力时程曲线
    Figure  7.  Pressure-time history curves of different cases under free field conditions

    与自由场条件下相比,模型2的气泡脉动情况较为复杂。表7列出了模型2的压力实验结果。当周围存在圆柱壁面结构时,气泡脉动会受到结构的影响,在膨胀阶段受到结构表面的轻微排斥,而在坍塌阶段被结构强烈吸引,气泡二次和三次脉动时释放的气泡能占比增大。与自由场条件下相比,当周围存在壁面结构时,气泡二次和三次脉动时的最大直径会增大。对比自由场电火花实验结果与模型2的压力实验结果可以看出,圆柱结构壁面对冲击波峰值压力的影响不大,主要影响气泡脉动压力和气泡射流载荷,如图8所示。

    表  7  模型2的压力实验结果
    Table  7.  Pressure experiment results of Model 2
    Case D/cm γ Shock wave
    pressure/MPa
    Pressure of secondary
    wave/MPa
    Pressure of triple
    wave/MPa
    Bubble jet
    load/MPa
    M-1 1 0.2 0.594 3.575 1.025
    M-2 2 0.4 0.466 2.085 1.297
    M-3 3 0.6 0.381 6.325 2.338 0.404
    M-4 4 0.8 0.284 7.317 3.032
    M-5 5 1.0 0.267 2.332 0.372
    下载: 导出CSV 
    | 显示表格
    图  8  自由场和模型2的压力-爆距曲线对比
    Figure  8.  Comparison of pressure-explosion distance curves between the free field and Model 2

    图9(a)~图9(b)所示,模型2中,当γ为0.2和0.4时,气泡二次脉动压力小于自由场条件下的二次峰值压力,且气泡溃灭后产生了气泡射流,对双层圆柱内壁面施加了射流载荷。结合图10所示的高速摄影机拍摄的脉动过程,分析发现:当γ=0.2 时,电火花装置起爆后首先产生的第1个峰值压力为冲击波压力;随后气泡开始膨胀进入第1个气泡脉动周期,由于气泡距离破口较近,气泡膨胀时发生了“腔吸效应”,气泡被分割成“舱内”气泡和“舱外”气泡,部分气泡能量随着“舱内”气泡的溃灭而释放,而“舱外”气泡体积相对于自由场中的气泡小,在5.0 ms时收缩到最小体积,产生的气泡二次脉动压力相较于自由场条件下也较小;受破口的限制,在7.1 ms,第2个气泡脉动周期未完成时气泡溃灭形成了射流,产生第3个峰值压力,即气泡射流载荷。

    图  9  不同爆距下模型2的压力时程曲线
    Figure  9.  Pressure-time history curves of Model 2 with different γ
    图  10  γ = 0.2时模型2的气泡脉动过程
    Figure  10.  Bubble pulsation process of Model 2 with γ = 0.2

    图11所示,当γ=0.6时,由于起爆点距离破口较远,破口对气泡形态的影响减弱,在2.5 ms时,气泡膨胀到最大体积,气泡在第1个脉动周期中保持完整的球形。但受壁面结构的限制,更多气泡能量在第2个脉动周期中释放,4.8 ms时产生的气泡二次脉动压力峰值大于自由场条件下的气泡二次脉动压力峰值。气泡上方仍有空间,允许气泡完成第2个脉动周期,气泡在7.0 ms时膨胀到第二次脉动周期中的最大体积,在7.5 ms时产生了气泡三次脉动压力,9.5 ms时溃灭形成射流产生气泡射流载荷,如图9(c)所示。

    图  11  γ = 0.6时模型2的气泡脉动过程
    Figure  11.  Bubble pulsation process of Model 2 with γ = 0.6

    γ=0.8 时,气泡脉动过程与γ=0.6的工况相近,此工况条件下气泡二次脉动压力峰值最大,产生了气泡三次脉动压力波,但内壁面没有检测到气泡射流载荷,如图9(d)所示。当γ=1.0时,气泡二次脉动峰值压力小于自由场条件下的二次压力波载荷,但存在气泡三次压力波,如图9(e)所示。此时,圆筒结构对气泡脉动的影响较小,模型2的压力实验结果与自由场接近。

    2.2.1   气泡脉动实验结果及分析

    不同γ下3个模型的水下电火花爆炸实验结果见表8表10。利用PIV lab软件对高速摄像机拍摄的图像进行后处理,得到气泡射流方向和速度。由于气泡溃灭生成射流时的速度在极短时间内变化较大,高速相机难以捕捉图像,因此选取拍摄到的气泡收缩至最小时的图像进行分析,取计算得到的最大速度为瞬时射流速度。通过测量水射流产生之后到达内壁面的行进距离以及花费的时间,计算得到速度,将其记为平均射流速度。表8表10中的射流速度均指朝向圆筒内壁即向上的射流速度。

    表  8  模型1的气泡脉动实验结果
    Table  8.  Bubble pulsation experiment results of Model 1
    Case γ Maximum bubble
    diameter/cm
    Bubble pulsation
    period/ms
    Jet direction Instantaneous jet
    velocity/(ms1)
    Mean jet velocity/
    (ms1)
    M1-1 0 4.50 4.0 Upward and downward 34.75 19.55
    M1-2 0.1 4.95 4.5 Upward and downward 41.44 25.56
    M1-3 0.2 5.08 4.9 Upward and downward 38.56 17.49
    M1-4 0.3 4.95 5.0 Upward and downward 38.88 12.68
    M1-5 0.4 5.02 5.1 Upward and downward 41.68 12.01
    M1-6 0.6 5.01 5.0 Upward 24.51 7.85
    M1-7 0.8 5.02 4.9
    下载: 导出CSV 
    | 显示表格
    表  9  模型2的气泡脉动实验结果
    Table  9.  Bubble pulsation experiment results of Model 2
    Case γ Maximum bubble
    diameter/cm
    Bubble pulsation
    period/ms
    Jet direction Instantaneous jet
    velocity/(ms1)
    Mean jet velocity/
    (ms1)
    M2-1 0 4.97 4.4 Upward and downward 44.74 25.08
    M2-2 0.1 4.95 4.8 Upward and downward 42.18 23.03
    M2-3 0.2 5.15 5.0 Upward and downward 44.31 17.88
    M2-4 0.3 5.05 4.9 Upward and downward 42.01 12.94
    M2-5 0.4 5.00 4.7 Upward 40.13 10.42
    M2-6 0.6 5.19 5.1 Upward 5.61
    M2-7 0.8 5.00 5.0
    下载: 导出CSV 
    | 显示表格
    表  10  模型3的气泡脉动实验结果
    Table  10.  Bubble pulsation experiment results of Model 3
    Case γ Maximum bubble
    diameter/cm
    Bubble pulsation
    period/ms
    Jet direction Instantaneous jet
    velocity/(ms1)
    Mean jet velocity/
    (ms1)
    M3-1 0 5.37 4.8 Upward 40.91 30.49
    M3-2 0.1 4.90 4.9 Upward 43.17 26.92
    M3-3 0.2 4.75 5.0 Upward 42.79 16.37
    M3-4 0.3 5.01 5.0 Upward 46.51 17.69
    M3-5 0.4 4.95 4.8 Upward 28.76 9.38
    M3-6 0.6 4.98 4.9 Upward 2.88
    M3-7 0.8 5.00 5.0
    下载: 导出CSV 
    | 显示表格

    爆距参数γ和破口参数λ是影响气泡脉动和气泡射流的重要参数。当γ0.6时,双层圆柱结构对气泡脉动的影响逐渐减小,气泡第1次脉动周期和气泡最大直径逐渐趋近于自由场条件下的脉动周期和最大直径。选取模型3(γ=0.6)作为典型工况进行分析,结果如图12所示。初始时刻铜丝通电,气泡随之产生并膨胀,2.5 ms时膨胀至最大,4.9 ms时收缩至最小,随后气泡进行二次脉动,二次脉动周期结束后气泡溃灭,9.5 ms时产生了向上的射流。

    图  12  γ = 0.6时模型3的气泡脉动过程
    Figure  12.  Bubble pulsation process of Model 3 with γ = 0.6

    γ0.4时,气泡在第1次脉动周期内会受到破口边界的阻碍作用,气泡分为舱内气泡和舱外气泡,如图13所示。舱内气泡还受到破口大小的影响,破口越小,对舱内气泡膨胀的限制越大。模型1破口直径最小,舱内气泡的运动范围也最小。

    图  13  γ = 0.2时不同模型破口附近的气泡脉动形态
    Figure  13.  Bubble pulsation patterns near different models with γ = 0.2

    图14所示,当γ=0 时,舱外气泡脉动会受到破口参数的影响。模型1的破口参数λ=0.2,舱外气泡膨胀受到破口的限制较大,舱外气泡的最大直径最小;模型2的破口参数λ=0.4,舱外气泡受破口的影响较小,舱外气泡的最大直径与自由场气泡的最大直径接近;模型3的破口参数λ=0.6,舱外气泡的最大直径大于自由场气泡的最大直径。

    图  14  γ = 0时不同模型破口附近的气泡脉动形态
    Figure  14.  Bubble pulsation patterns near different models with γ = 0
    2.2.2   气泡射流实验结果及分析

    爆距参数γ决定了是否能够产生有效射流以及射流速度,如图15所示,选取工况M3-6和M3-7作为典型工况进行分析,可以看出:γ0.8 时,射流无法到达内壁面,属于无效射流;当γ=0.6 时,气泡溃灭后形成的射流在到达内壁面时速度已经大幅衰减,毁伤效果减弱。

    图  15  工况M3-6和工况M3-7气泡的射流形态
    Figure  15.  Bubble jet pattern of Case M3-6 and Case M3-7

    PIV测量的不同γ下气泡的射流速度流场结果见表8表10。当γ0.6 时,射流的成因可以用著名的Bjerknes效应来解释。在Bjerknes力作用下,气泡会被水中壁面结构吸引而向壁面结构方向运动,由于气泡在圆柱模型下方生成,因此在Bjerknes力、浮力、惯性力等因素的影响下, 气泡脉动时会向上朝圆柱模型运动,气泡下表面往上收缩,沿着气泡坍塌方向形成高速射流。当γ=0.8 时,模型1~模型3附近均不产生有效射流。当γ0.4λ0.4γ0.3λ=0.6时,受破口的限制,气泡一次脉动后就溃灭形成射流,最大射流速度在34.75~46.51 m/s范围内。

    几种典型工况下气泡的射流速度测量结果如图16所示。工况M1-1中,气泡的瞬时射流速度最小,这是由于起爆点距离破口最近且破口直径最小,气泡脉动受到的限制最大。工况M3-5中,当γ=0.4λ=0.6时,气泡距离结构较远,且破口直径较大,气泡上方有充足的空间可以进行二次脉动,气泡完成二次脉动周期后溃灭形成气泡射流,二次脉动释放了部分能量,此时气泡的瞬时射流速度减小至28.76 m/s,如图16(d)所示。如图16(e)所示,当γ=0.6时,模型1~模型3破口附近气泡均进入二次脉动周期,测得工况M1-6的瞬时射流速度为24.51 m/s,但模型2、模型3破口附近在二次脉动周期时气泡运动的范围过小,软件无法分析,无瞬时射流速度。

    图  16  几种典型工况下瞬时射流速度的PIV结果
    Figure  16.  Jet velocities obtaied by PIV for several typical cases

    γ较小时,λ决定了射流的方向。当λ0.6 时,舱内气泡受到破口的限制作用较小,与舱外气泡同时脉动并溃灭,形成向上的射流。当λ0.4 时,出现向上、向下及其他不同方向的水射流,这是由于舱内气泡受到破口的限制作用较大,其收缩速度大于舱外气泡,舱内气泡先溃灭形成向下的射流,击穿舱外气泡,随后舱外气泡继续收缩,最终形成向上的射流,如图17所示。不同工况下射流的形成条件如图18所示。模型1的λ=0.2,当γ0.4时,有向上、向下射流产生;当0.4<γ0.6 时,只产生向上射流;γ0.8 时,不产生有效射流。模型2的λ=0.4,当γ0.3时,有向上、向下射流产生;当0.3<γ0.6 时,只产生了向上射流;当γ0.8时,不产生有效射流。模型3的λ=0.6,当爆距参数γ0.6时,只产生向上的有效射流;当γ0.8 时,不产生有效射流。从表8表10中的数据可以看出,平均射流速度随着γ的增大而减小,这是因为随着射流向前运动,射流速度不断衰减。

    图  17  γ = 0.2时不同模型破口附近的气泡射流形态
    Figure  17.  Bubble jet pattern near different models with γ = 0.2
    图  18  不同模型破口附近射流的形成条件
    Figure  18.  Conditions of jet formation near different models

    基于PIV技术对带破口双层圆柱结构附近气泡脉动和气泡载荷特性开展了实验研究,探究了爆距参数和破口参数对气泡脉动和流场特征的影响,得出如下结论。

    (1) 带破口的双层圆柱结构对水下爆炸冲击波峰值压力的影响较小,主要影响气泡脉动压力和射流载荷压力。当γ=0.2, γ=0.4时,气泡的二次压力波相对于自由场工况小,且产生了射流载荷;当γ=0.6, γ=0.8时,气泡的二次压力波相对于自由场工况大,且存在三次压力波;当γ=1.0时,压力结果与自由场结果相似。

    (2) 爆距参数γ和破口参数λ是影响气泡脉动和气泡射流的2个重要参数。当γ0.4 时,气泡脉动的第1个周期与自由场大致相同;当γ<0.4 时,脉动受到破口边界的阻碍作用,产生的气泡分为舱内气泡和舱外气泡;当γ=0 时,舱外气泡的脉动受破口参数λ的影响;当λ=0.6 时,舱外气泡的最大直径大于自由场条件下气泡的最大直径;当λ0.2 时,舱外气泡的膨胀受到破口的限制较大,舱外”气泡的最大直径小于自由场气泡的最大直径。

    (3) 采用PIV技术,获得了不同时刻气泡周围水介质的速度矢量图,测量得到不同γ下双层圆柱破口附近的瞬时最大水射流速度为24.51~46.51 m/s,平均水射流速度(2.88~30.49 m/s)随着爆距的增大而减小。

    (4) 爆距参数γ决定了是否产生有效射流。当γ0.6 时,气泡溃灭产生的射流能够到达双层圆柱内壁面,内壁面受到气泡射流压力的影响;当γ0.8 时,气泡射流无法到达圆筒内壁面,属于无效射流。

    (5) 当γ0.4 时,破口参数λ决定了射流的方向。当λ0.6 时,气泡溃灭只产生向上的射流;当λ0.4 时,舱内气泡先溃灭形成向下的射流,射流击穿舱外气泡后,舱外气泡继续收缩,最终形成向上的射流。

  • 图  煤样图片及尺寸

    Figure  1.  Picture and size of the coal sample

    图  实验方案

    Figure  2.  Experimental scheme

    图  垂直层理方向渗流孔分形维数的对比

    Figure  3.  Fractal dimension of the seepage hole in the vertical bedding direction

    图  平行层理方向渗流孔分形维数的对比

    Figure  4.  Fractal dimension of the seepage hole in the parallel bedding direction

    图  斜交层理方向渗流孔分形维数的对比

    Figure  5.  Fractal dimension of the seepage hole in the oblique bedding direction

    图  垂直层理方向吸附孔分形维数的对比

    Figure  6.  Fractal dimension of the adsorption hole in the vertical bedding direction

    图  平行层理方向吸附孔分形维数的对比

    Figure  7.  Fractal dimension of the adsorption hole in the parallel bedding direction

    图  斜交层理方向吸附孔分形维数的对比

    Figure  8.  Fractal dimension of the adsorption hole in the oblique bedding direction

    图  不同方向煤体渗流孔与吸附孔分形维数减小率随冲击荷载的变化

    Figure  9.  Fractal dimension reduction rates of the seepage hole and adsorption hole vary with the impact load

    图  10  不同方向煤体的分形维数与冲击荷载之间的关系

    Figure  10.  Relationship between the fractal dimension and the impact load of coal bodies in different directions

    表  1  煤样基本参数

    Table  1.   Basic parameters of coal samples

    No.ρ/(g·cm−3)fc/MPa CL/(km·s−1)R0,max/%Mass fraction/%
    AshVolatile
    component
    Fixed carbonExiniteVitriniteInertinite
    C01.42317.291.3813.325.626.0583.981981
    P01.46115.131.852
    X01.42211.651.653
    下载: 导出CSV

    表  2  实验数据表

    Table  2.   Impact test data

    DirectionpI/MPapL/MPaNo.Size/mm × mmDsˉDsR2sDxˉDxR2xdˉd
    Vertical
    bedding
    00YC12.992.980.972.752.780.95
    YC22.990.992.780.96
    YC32.950.982.810.98
    0.1028.46C1 49.75 × 50.332.852.850.952.312.230.950.140.14
    C249.33 × 49.892.920.962.210.930.16
    C349.34 × 50.422.790.952.170.960.12
    0.1532.59C449.42 × 49.572.732.730.992.192.190.950.180.16
    C549.66 × 49.932.780.972.230.960.12
    C649.42 × 50.842.690.942.150.950.18
    0.2041.43C749.39 × 49.932.632.682.052.130.960.180.17
    C849.42 × 49.762.710.982.210.970.17
    C949.45 × 50.282.690.992.130.980.15
    0.3051.80C1049.52 × 50.472.562.510.992.052.050.990.170.19
    C1149.36 × 49.782.640.962.090.960.19
    C1249.27 × 50.272.340.942.010.960.22
    0.5058.70C1349.38 × 50.472.742.730.922.112.120.920.330.44
    C1449.45 × 50.212.790.942.170.930.57
    C1549.49 × 49.442.670.952.070.980.42
    Parallel
    bedding
    00YP12.842.870.992.892.860.95
    YP22.910.962.860.94
    YP32.870.972.830.96
    0.1028.46P149.29 × 50.112.572.570.992.032.030.980.130.13
    P249.36 × 50.372.650.982.050.950.11
    P349.47 × 49.882.500.992.020.960.15
    0.1532.59P449.31 × 49.322.582.580.972.152.151.000.310.25
    P549.33 × 50.282.610.982.190.950.22
    P649.31 × 49.162.560.952.110.970.22
    Parallel
    bedding
    0.2041.43P749.33 × 50.332.642.640.922.202.190.950.180.28
    P849.60 × 49.852.680.932.140.960.22
    P949.31 × 50.112.610.952.240.980.44
    0.3051.80P1049.31 × 50.402.812.810.942.212.180.950.370.34
    P1149.05 × 49.872.850.952.160.970.30
    P1249.42 × 49.872.780.942.180.960.35
    0.5058.70P1349.40 × 49.932.892.810.952.822.730.970.540.50
    P1449.62 × 49.962.810.992.730.960.49
    P1549.32 × 49.792.730.952.630.960.47
    45° oblique bedding00YX12.812.810.942.772.780.95
    YX22.850.952.710.95
    YX32.780.972.850.96
    0.1028.46X149.49 × 50.032.752.750.992.382.380.970.170.16
    X249.27 × 49.632.810.982.280.950.19
    X349.31 × 49.302.680.972.470.960.12
    0.1532.59X449.24 × 49.932.672.650.992.392.310.980.290.24
    X549.27 × 49.332.750.952.310.960.21
    X649.27 × 49.822.540.982.230.950.22
    0.2041.43X749.25 × 49.902.762.700.952.082.120.940.200.22
    X849.60 × 49.382.710.922.170.920.21
    X949.29 × 49.052.640.932.120.930.25
    0.3051.80X1049.40 × 50.282.472.550.982.012.060.930.210.18
    X1149.38 × 50.882.610.962.070.920.17
    X1249.71 × 50.022.580.972.110.950.16
    0.5058.70X1349.66 × 50.442.792.730.982.262.280.960.630.71
    X1449.73 × 49.912.750.972.310.940.72
    X1549.46 × 50.052.640.962.280.930.78
    下载: 导出CSV
  • [1] 高为, 易同生, 金军, 等. 黔西地区煤样孔隙综合分形特征及对孔渗性的影响 [J]. 煤炭学报, 2017, 42(5): 1258–1265.

    GAO W, YI T S, JIN J, et al. Pore integrated fractal characteristics of coal sample in western Guizhou and its impact to porosity and permeability [J]. Journal of China Coal Society, 2017, 42(5): 1258–1265.
    [2] 秦勇, 邱爱慈, 张永民. 高聚能重复强脉冲波煤储层增渗新技术试验与探索 [J]. 煤炭科学技术, 2014, 42(6): 1–7.

    QIN Y, QIU A C, ZHANG Y M. Experiment and discovery on permeability improved technology of coal reservoir based on repeated strong pulse waves of high energy accumulation [J]. Coal Science and Technology, 2014, 42(6): 1–7.
    [3] WEI G Y, SHAN Z Y, ZHANG Z M. Research on hydraulic slotting technology controlling coal-gas outbursts [J]. Journal of Coal Science and Engineering, 2008, 14(1): 67–72. doi: 10.1007/s12404-008-0014-4
    [4] 姜永东, 李业, 崔悦震, 等. 声场作用下煤储层渗透性试验研究 [J]. 煤炭学报, 2017, 42(S1): 154–159.

    JIANG Y D, LI Y, CUI Y Z, et al. Experimental study on characteristics of coal reservoir permeability under acoustic wave [J]. Journal of China Coal Society, 2017, 42(S1): 154–159.
    [5] 赵丽娟, 秦勇. 超声波作用对改善煤储层渗透性的试验分析 [J]. 天然气地球科学, 2014, 25(5): 747–752.

    ZHAO L J, QIN Y. Experiment on improving the permeability of coal reservoir under ultrasound [J]. Natural Gas Geoscience, 2014, 25(5): 747–752.
    [6] FORRESTAL M J, FREW D J, HANCHAK S J, et al. Penetration of grout and concrete targets with ogive-nose steel projectiles [J]. International Journal of Impacting Engineering, 1996, 18(5): 465–476. doi: 10.1016/0734-743X(95)00048-F
    [7] 穆朝民, 王海露, 黄文尧, 等. 高瓦斯低透气性煤体定向聚能爆破增透机制 [J]. 岩土力学, 2013, 34(9): 2496–2500.

    MU C M, WANG H L, HUANG W Y, et al. Increasing permeability mechanism using directional cumulative blasting in coal seams with high concentration of gas and low permeability [J]. Rock and Siol Mechanics, 2013, 34(9): 2496–2500.
    [8] 李恒乐, 秦勇, 张永民, 等. 重复脉冲强冲击波对肥煤结构影响的实验研究 [J]. 煤炭学报, 2015, 40(4): 915–921.

    LI H L, QIN Y, ZHANG Y M, et al. Experimental study on the effect of strong repetitive pulse shockwave on the pore structure of fat coal [J]. Journal of China Coal Society, 2015, 40(4): 915–921.
    [9] CAI Y D, LIU D M, PAN Z J, et al. Pore structure and its impact on CH4 adsorption capacity and flow capability of bituminous and subbituminous coals from Northeast China [J]. Fuel, 2013, 103: 258–268. doi: 10.1016/j.fuel.2012.06.055
    [10] MANDELBROT B B. The fractal geometry of nature [M]. San Francisco: Freeman, 1982: 35.
    [11] BIRD N, DÍAZ M C, SAA A, et al. Fractal and multifractal analysis of pore-scale images of soil [J]. Journal of Hydrology, 2006, 322(1/2/3/4): 211–219.
    [12] DATHE A, EINS S, NIEMEYER J, et al. The surface fractal dimension of the soil-pore interface as measured by image analysis [J]. Geoderma, 2001, 103(1/2): 203–229.
    [13] HUANG S J, YU Y C, LEE T Y, et al. Correlations and characterization of porous solids by fractal dimension and porosity [J]. Physica A: Statistical Mechanics and Its Applications, 1999, 274(3/4): 419–432.
    [14] HILDGEN P, NEKKA F, HILDGEN F, et al. Macroporosity measurement by fractal analysis [J]. Physica A: Statistical Mechanics and Its Applications, 1997, 234(3/4): 593–603.
    [15] YANG F, NING Z F, LIU H Q. Fractal characteristics of shales from a shale gas reservoir in the Sichuan Basin, China [J]. Fuel, 2014, 115: 378–384. doi: 10.1016/j.fuel.2013.07.040
    [16] LIU X J, XIONG J, LIANG L X. Investigation of pore structure and fractal characteristics of organic-rich Yanchang formation shale in central China by nitrogen adsorption/desorption analysis [J]. Journal of Natural Gas Science and Engineering, 2015, 22(1): 62–72.
    [17] 朱汉卿, 贾爱林, 位云生, 等. 蜀南地区五峰-龙马溪组页岩微观孔隙结构及分形特 [J]. 科学技术与工程, 2018, 18(10): 12–19. doi: 10.3969/j.issn.1671-1815.2018.10.003

    ZHU H Q, JIA A L, WEI Y S, et al. Microscopic pore structure and fractal characteristics of Wufeng-Longmaxi shale, south Sichuan [J]. Science Technology and Engineering, 2018, 18(10): 12–19. doi: 10.3969/j.issn.1671-1815.2018.10.003
    [18] GARBACZ J K. Fractal description of partially mobile single gas adsorption on energetically homogeneous solid adsorbent [J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 1998, 143(1): 95–101.
    [19] LEE G J, PYUN S Ⅰ, RHEE C K. Characterization of geometric and structural properties of pore surfaces of reactivated microporous carbons based upon image analysis and gas adsorption [J]. Microporous and Mesoporous Materials, 2006, 93(1/2/3): 217–225.
    [20] YAO Y B, LIU D M, TANG D Z, et al. Fractal characterization of adsorption-pores of coals from North China: an investigation on CH4 adsorption capacity of coals [J]. International Journal of Coal Geology, 2008, 73(1): 27–42. doi: 10.1016/j.coal.2007.07.003
    [21] XU L J, ZHANG D J, XIAN X F. Fractal dimensions of coals and cokes [J]. Journal of Colloid and Interface Science, 1997, 190(2): 357–359. doi: 10.1006/jcis.1997.4885
    [22] 宋昱, 姜波, 李凤丽, 等. 低-中煤级构造煤纳米孔分形模型适用性及分形特征 [J]. 地球科学, 2018, 43(5): 1611–1622.

    SONG Y, JIANG B, LI F L, et al. Applicability of fractal models and nanopores’ fractal characteristics for low-middle rank tectonic deformed coals [J]. Earth Science, 2018, 43(5): 1611–1622.
    [23] 王登科, 刘淑敏, 魏建平, 等. 冲击载荷作用下煤的破坏特性试验研究 [J]. 采矿与安全工程学报, 2017, 34(3): 594–600.

    WANG D K, LIU S M, WEI J P, et al. The failure characteristics of coal under impact load in laboratory [J]. Journal of Mining & Safety Engineering, 2017, 34(3): 594–600.
    [24] 刘晓辉, 张茹, 刘建锋. 不同应变率下煤岩冲击动力试验研究 [J]. 煤炭学报, 2012, 37(9): 1528–1534.

    LIU X H, ZHANG R, LIU J F. Dynamic test study of coal rock under different strain rates [J]. Journal of China Coal Society, 2012, 37(9): 1528–1534.
    [25] XODOT B B. 煤与瓦斯突出 [M]. 宋世钊, 王佑安, 译. 北京: 中国工业出版社, 1966: 310–318.
    [26] YUAN J, TAKEDA N, WAAS A M. A note on data processing in the split Hopkinson pressure bar tests [J]. Experimental Techniques, 1998, 22(5): 21–24. doi: 10.1111/j.1747-1567.1998.tb02317.x
    [27] SHAN R L, JIANG Y S, LI B Q. Obtaining dynamic complete stress-strain curves for rock using the split Hopkinson pressure bar technique [J]. International Journal of Rock Mechanics and Mining Sciences, 2000, 37(6): 983–992. doi: 10.1016/S1365-1609(00)00031-9
    [28] 陈俊宇, 裴向军, 杜瑞锋, 等. 冲击载荷作用下砂岩的动力学特性及能耗规律 [J]. 科学技术与工程, 2019, 19(31): 304–310. doi: 10.3969/j.issn.1671-1815.2019.31.045

    CHEN J Y, PEI X J, DU R F, et al. Dynamic characteristics and energy consumption of sandstone under impact loading [J]. Science Technology and Engineering, 2019, 19(31): 304–310. doi: 10.3969/j.issn.1671-1815.2019.31.045
    [29] KOLSKY H. An investigation of the mechanical properties of materials at very high rates of loading [J]. Proceedings of the Physical Society (Section B), 1949, 62: 676–700. doi: 10.1088/0370-1301/62/11/302
    [30] LUNDBERG B. A split Hopkinson bar study of energy absorption in dynamic rock fragmentation [J]. International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts, 1976, 13(6): 187–197.
    [31] 赵迪斐, 郭英海, WANG G, 等. 基于分形建模的高煤级煤孔隙结构特征量化表征-以阳泉矿区山西组煤样为例 [J]. 东北石油大学学报, 2019, 43(3): 53–67. doi: 10.3969/j.issn.2095-4107.2019.03.006

    ZHAO D F, GUO Y H, WANG G, et al. Quantitative characterization of pore structure characteristics of high quality coal based on fractal modeling: taking coal samples from Shanxi Formation in Yangquan Mining Area as an example [J]. Journal of Northeast Petroleum University, 2019, 43(3): 53–67. doi: 10.3969/j.issn.2095-4107.2019.03.006
    [32] RIGBY S P. Predicting surface diffusivities of molecules from equilibrium adsorption isotherms [J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2005, 262(1/2/3): 139–149.
    [33] 杨峰, 宁正福, 张世栋, 等. 基于氮气吸附实验的页岩孔隙结构表征 [J]. 天然气工业, 2013, 33(4): 135–140. doi: 10.3787/j.issn.1000-0976.2013.04.025

    YANG F, NING Z F, ZHANG S D, et al. Characterization of pore structures in shales through nitrogen adsorption experiment [J]. Natural Gas Industry, 2013, 33(4): 135–140. doi: 10.3787/j.issn.1000-0976.2013.04.025
    [34] 吉小峰. 煤中纳米孔隙发育特征及其对气体运移的控制机理研究 [D]. 焦作: 河南理工大学, 2018.

    JI X F. Development characteristics of nanopores in coal and its controlling mechanism on gas migration [D]. Jiaozuo: Henan Polytechnic University, 2018.
    [35] YAO Y B, LIU D M, TANG D Z, et al. Fractal characterization of seepage-pores of coals from China: an investigation on permeability of coals [J]. Computers & Geosciences, 2009, 35(6): 1159–1166.
    [36] 宋晓夏, 唐跃刚, 李伟, 等. 中梁山南矿构造煤吸附孔分形特征 [J]. 煤炭学报, 2013, 38(1): 134–139.

    SONG X X, TANG Y G, LI W, et al. Fractal characteristics of adsorption pores of tectonic coal from Zhongliangshan southern coalmine [J]. Journal of China Coal Society, 2013, 38(1): 134–139.
    [37] 林海飞, 刘静波, 严敏, 等. CO2/CH4在煤储层中扩散规律的分子动力学模型 [J]. 中国安全生产科学技术, 2017, 13(1): 84–89.

    LIN H F, LIU J B, YAN M, et al. Molecular dynamics simulation on diffusion rules of CO2/CH4 in coal reservoir [J]. Journal of Safety Science and Technology, 2017, 13(1): 84–89.
    [38] 董怡静, 韩雨桢, 候泉林, 等. 煤变形产气的力化学机理探讨 [J]. 煤炭学报, 2017, 42(4): 942–949.

    DONG Y J, HAN Y Z, HOU Q L, et al. Mechanochemistry mechanism of gas generation during coal deformation [J]. Journal of China Coal Society, 2017, 42(4): 942–949.
    [39] 刘运通, 高文学. 爆炸荷载下岩石损伤的数值模拟研究 [J]. 岩石力学与工程学报, 2001, 20(6): 789–792. doi: 10.3321/j.issn:1000-6915.2001.06.007

    LIU Y T, GAO W X. Numerical simulation on rock damage under explosion loading [J]. Chinese Journal of Rock Mechanics and Engineering, 2001, 20(6): 789–792. doi: 10.3321/j.issn:1000-6915.2001.06.007
    [40] 唐红梅, 周云涛, 廖云平. 地下工程施工爆破围岩损伤分区研究 [J]. 振动与冲击, 2015, 34(23): 202–206.

    TANG H M, ZHOU Y T, LIAO Y P. Damage zone of surrounding rock of underground engineering under construction blasting [J]. Journal of Vibration and Shock, 2015, 34(23): 202–206.
  • 加载中
图(10) / 表(2)
计量
  • 文章访问数:  6595
  • HTML全文浏览量:  2432
  • PDF下载量:  20
出版历程
  • 收稿日期:  2020-03-12
  • 修回日期:  2020-03-26

目录

/

返回文章
返回