-
砂石等脆性材料由大量的脆性散体颗粒组成,具有松散、多孔、易压缩等特性[1-2],应力波在此类松散脆性材料中的传播[3]通常会呈现显著的衰减规律。散体脆性颗粒作为消波耗能材料已被广泛应用于防护领域,用作冲击、爆炸等的能量吸收[4-5]。因此,深入探讨应力波在散体脆性颗粒中的传播规律和影响因素对防护工程具有重大意义,亦可为防护层的构建提供一定的理论参考。
季顺迎等[6]用石质球体作为冲击物,在散体颗粒底部安装压力传感器,通过散体颗粒底部压力的变化来分析其缓冲性能,研究了颗粒粒径、颗粒形态等对缓冲性能的影响,由于石质球与颗粒接触过程存在不均匀性,因此实验结果的精度不高。赵跃堂等[7]基于平面波加载器,通过在不同深度土层中埋设压力传感器的方法来探究饱和土颗粒在爆炸加载下不同层的压力变化,根据不同土层的压力分布来分析波的传播规律,由于在土颗粒中埋设压力传感器会改变土颗粒的内部结构,因此该方法测得的压力值与真实值存在一定差异。赵凯[8]实验研究了爆炸载荷作用下多层防护层结构对爆炸波的衰减作用,获得了多种分配层下波的衰减规律。魏久淇等[9]基于分离式霍普金森压杆实验,通过透射杆峰值应力与入射杆峰值应力的比值来描述波的衰减,并探究了中高应变率下钙质砂的力学响应。由于钙质砂与透射杆之间存在较大的机械阻抗,直接将透射杆峰值应力作为砂粒端面应力的处理方式对波衰减的准确性有一定影响。Yu等[10]通过改进SHPB装置,将入射波、反射波、透射波相结合来计算试件两端的峰值应力,提出了一种定量计算珊瑚砂试件中应力波衰减规律的新方法,能更准确地得到试件两端面的应力峰值,由于该数据分析方法是基于一维单向应力波传播理论,而应力波在散体颗粒中传播时可能会发生透射和反射,将涉及更复杂的多向波系,因此该方法存在一定误差。
目前,应力波在散体颗粒中传播规律的实验研究尚不成熟,而数值模拟可以避免实验存在的一些技术难题,定性或半定量地分析散体颗粒中波的传播规律。祁原等[11]利用离散元软件PFC2D分析了可破碎颗粒在循环载荷下的耗能特性,并探究了颗粒破碎与系统能量耗散率的关系。郑文[12]、黄俊宇[13]通过PFC3D建立散体颗粒模型,研究了石英砂颗粒在动载荷下的力学响应。迄今为止,数值模拟研究更多关注散体颗粒本构模型的构建和动态力学性能,散体颗粒中应力波的传播规律尚缺乏系统研究。
本研究基于离散元颗粒流软件PFC3D建立典型的石英砂颗粒模型,在保证模型颗粒与实际颗粒等粒径的前提下,通过模型的等比例缩放,模拟侧限条件下石英砂的动态压缩实验,并与实验侧限压缩本构[14-16]进行对比,验证本模型参数的合理性,在此基础上研究冲击载荷下应力波在散体颗粒中的传播规律和影响因素。
-
PFC3D模型是由若干基础颗粒(刚性颗粒)通过一定的黏结键组成的,通过定义基础颗粒间的微观参数来建立与模拟材料具有相同宏观力学性能的离散元模型。由于离散元模型中的微观参数难以通过实验手段直接获得,且尚无有效的理论将其与实际材料的宏观力学特性进行关联,因此目前均通过数值模拟校准的方法建立二者之间的联系[17]。对于散体颗粒,土力学中常采用侧限压缩本构模型[18-20]来表示其本构关系,即通过孔隙比-轴向应力曲线来描述颗粒材料的压缩、变形及破碎过程[21]。
建立如图1所示的簇颗粒模型来表述单个脆性颗粒,每个簇由60~70个基础颗粒组成,由于实际砂颗粒内部存在缺陷和孔隙,删除10%~20%的基础颗粒来模拟实际颗粒的不规则外形和风化腐蚀等对砂砾造成的损伤,通过调整簇中基础颗粒的大小和颗粒间的孔隙以满足实际颗粒尺寸,并保证簇的密度同实际砂粒密度一致。选取石英砂为典型的砂砾代表,其密度为2 636 kg/m3。通过Flat-Joint黏结模型将松散的基础颗粒用黏结键黏结,当基础颗粒间的法向接触力或切向接触力超过其临界值时,黏结键将发生破坏失效。表1为石英砂采用Flat-Joint黏结模型时所对应的微观参数,其中Emod为接触模量,Fj-Kn为法向刚度,Fj-Ks为切向刚度,Fj-Coh为切向黏结强度,Fj-Ten为法向黏结强度,Damp为阻尼系数。
Model Emod/GPa Fj-Kn/GPa Fj-Ks/GPa Fj-Coh/GPa Fj-Ten/GPa Damp Flat-Joint 50 5.25 × 105 9.37 × 105 0.2 0.05 0.15 表 1 Flat-Joint模型微观参数
Table 1. Micro-parameters of the Flat-Joint model
单个石英砂颗粒的粒径分布控制在150~250 μm之间,以42%的自然堆积孔隙率生成散体颗粒,石英砂颗粒内部采用Flat-Joint黏结键黏结,颗粒之间为线性接触。在保证粒径大小同实际一致的情况下对模型进行等比例缩放,生成直径为2.2 mm、高为0.6 mm的圆柱体试样,如图2所示。圆柱体试样的上、下压板层设置为刚性墙体,径向为刚性约束,保持下底板不动,对上压板向下持续施加10 m/s的加载速度,得到图2中的归一化孔隙比-轴向应力曲线,模拟结果与文献[13]中的实验结果整体吻合较好,屈服应力均为10 MPa左右。
-
如图3所示,建立直径为1 mm、高为3.6 mm的圆柱试件,径向为刚性约束,下端面为自由面,上端面向下施加一个速度为50 m/s的冲击加载,由于颗粒间应力分布的离散性较大,在试件内部引入测量球,通过测量球所在空间位置全部粒子的应力平均值来描述应力波幅值,测量球球心所在位置为应力波传播的距离,通过不同位置测量球测得的应力幅值来表征应力波在散体颗粒中的传播规律。图4给出了试件不同位置的应力波波形。整体而言,随着传播距离的增大,应力波峰值逐渐减小,应力波脉宽逐渐发散。
采用观察区域内的应力平均值来表征应力波的传播规律,测量球对应力波传播可能存在一定影响。由于研究内容为轴向的应力波在散体颗粒中的传播规律,因此将测量球的球心固定在圆柱试件的中心轴上,通过改变测量球的直径大小d来研究测量球对应力波传播规律的影响。图5给出了通过不同直径大小的测量球获得的应力波衰减规律。整体来看,所有衰减规律的变化趋势是一致的,但测量球直径越大,应力波衰减越严重。测量球获得的平均应力并非一维应力波,同时涉及了圆柱径向的应力波衰减,故测量球直径越大,横向衰减考虑的也越多。归一化峰值应力
$\bar \sigma = {\sigma _{\rm{s}}}/{\sigma _0}$ ,无量纲化应力波传播距离$\bar L = {L_{\rm{s}}}/D$ ,其中,${\sigma _{\rm{s}}}$ 为当前峰值应力,${\sigma _0}$ 为初始峰值应力,Ls为应力波传播的距离,D为圆柱试件的直径。如果测量球直径过小,由于散体颗粒较大的分散性导致测量区域的应力波动较大;而当选取的测量球直径过大,特别是接近试件直径时,将导致测量的区域靠近侧面固壁端的应力加倍。本研究选取测量球直径为试件直径的0.6倍,此时应力波衰减趋于稳定。从图5中可以看到,峰值应力随传播距离呈指数型衰减,具体可表示为
$\bar \sigma = A + B \times {{\rm{e}}^{ - \beta \times \bar L}}$ 式中:A = 0.20,B = 0.76,β = 2.82。β作为衰减因子,可以定量地描述应力波在散体颗粒中的衰减程度,β越大,说明应力波衰减程度越大。
为更精确地研究应力波的衰减规律,在散体颗粒圆柱中引入更多测量球来观测试件不同位置的应力波形,取各处应力波峰值,归一化后绘于图6(a)。整体而言,归一化峰值应力仍然呈现近似指数型衰减规律,通过对局部数据点拟合可以得出,在传播前期应力波衰减明显,对应的衰减因子β = 5.40;而在传播的中后期,应力波衰减程度趋缓,对应的衰减因子β = 1.16,远小于前期的衰减因子。对照应力波传播过后的颗粒细观结构,如图6(b)所示,其中红色区域代表颗粒间黏结键的断裂。前期应力波衰减显著的区域正好对应于颗粒大量破碎的区域,将其定义为破碎区,中后期衰减因子较小的区域正好对应于颗粒较少甚至没有发生破碎的区域,将其定义为微/未破碎区。由图6可知,颗粒的破碎区主要集中在试件的加载端,而恰好加载端附近区域的峰值应力衰减最为严重,远离加载端的微/未破碎区的颗粒破碎程度显著低于破碎区,对应的应力波峰值应力衰减不明显,趋于稳定,说明颗粒的破碎和应力波衰减具有直接关系。
-
将图3离散元模型中颗粒内部的黏结键强度设置为1 000 MPa(正常情况为50 MPa),在数值模拟过程中黏结键将不发生断裂,即颗粒不破碎,保持其他微观参数不变,在50 m/s的加载速度下进行数值模拟对比。将应力波在不同位置的峰值应力绘于图7,并用式(1)进行拟合,可得不破碎颗粒中的应力波衰减因子β = 1.20,而正常黏结强度下的可破碎颗粒在微/未破碎区中应力波衰减因子β = 1.16,两者相差无几。
图 7 不同区域拟合曲线无量纲传播距离
Figure 7. Dimensionless propagation distance of fitted curves in different regions
将正常黏结强度下的可破碎颗粒在破碎区中的应力波峰值应力进行拟合,得到应力波衰减因子β = 5.40,显著大于不可破碎颗粒的应力波衰减因子(β = 2.92),如图8所示。颗粒破碎过程中,基础颗粒间的黏结键断裂需要消耗能量,并且颗粒破碎使整个散体颗粒系统中的力链更加复杂,由于应力波是单向传播,且颗粒间具有复杂的应力波透反射现象,因此宏观上表现为应力波峰值应力大幅衰减,卸载波形严重弥散。
为深入研究应力波衰减过程中局部衰减因子β与颗粒局部破碎程度的关系,设定基础颗粒间的黏结键强度范围为20~100 MPa,数值模拟得到颗粒在破碎区的局部破碎程度为8%~60%,对应破碎区局部的应力波衰减曲线指数拟合得到β为3.50~6.10。如图9所示,β和颗粒的破碎程度具有显著的正相关性,即β随颗粒破碎程度的增加而增大,近似呈指数关系,说明颗粒破碎是影响散体颗粒中应力波衰减的一个重要因素。
-
应力波在散体颗粒间传播时会发生透射和反射,涉及复杂的多向波系,在弹性散体颗粒中会发生散射现象,即随着传播距离的增大,应力波脉宽会弥散,如图4所示。应力波的衰减程度随着传播距离的增大而逐渐趋缓,除了受到颗粒破碎程度的影响之外,应力波的弥散也会对其产生一定的影响。为避免颗粒破碎带来的干扰,对不可破碎散体颗粒模型(黏结键强度为1 000 MPa)进行不同脉宽的冲击加载。图10给出了相同峰值应力、不同脉宽的应力波随传播距离的衰减规律。由图10可知,加载波形脉宽越大,应力波随传播距离增大所产生的峰值应力衰减越少,弥散也越小,加载初始,应力波具有相近峰值应力而脉宽不同,当传播一定距离之后,应力波具有相近脉宽而峰值应力不同。对不同脉宽的应力波传播过程中的峰值应力利用式(1)进行拟合,得到如图11所示的拟合曲线和对应的衰减因子β。由图11可知,随着应力波脉宽的增大,β从1.46减小到0.61,应力波脉宽亦能对应力波在散体颗粒中的传播产生较大影响。造成该现象的主要原因是:散体颗粒对高频波具有显著的弥散作用,会大幅削弱短脉宽的高频波峰值应力,而对于长脉宽的低频波影响相对较小。
-
图12给出了散体颗粒在不同加载速度(30~110 m/s)下,应力波峰值随传播距离的衰减规律。从图12中可以看出,加载速度越大,β值越大,即应力波的衰减程度随着加载速度的增大而增大,在靠近加载端尤为明显。考虑到颗粒破碎是应力波衰减的重要影响因素,将颗粒内部的黏结键强度设置为1 000 MPa,此时,颗粒不发生破碎,拟合曲线如图13所示。图13结果表明,在颗粒不发生破碎的前提下,在一个较宽泛的加载速度(1~100 m/s)范围内,应力波在散体颗粒中的衰减规律基本不变。
图 12 不同加载速度下可破碎颗粒的峰值应力拟合曲线
Figure 12. Peak stress fitting curves of the crushable particlesat different loading velocities
图 13 不同加载速度下不可破碎颗粒的峰值应力拟合曲线
Figure 13. Peak stress fitting curves of the uncrushable particlesat different loading velocities
将加载速度、衰减因子β、颗粒黏结键的破碎比建立联系,如图14所示。由图14可知,随着加载速度的增大,β呈上升趋势,应力波的衰减程度变大,观察黏结键的破碎比发现,随着加载速度的增大,其破碎比也呈增大趋势,且β与黏结键的破碎比有相似的变化趋势。观察图14中不破碎颗粒的衰减因子,可以看到不同加载速度对不破碎颗粒的应力波衰减影响较小。由此可知,应力波衰减的速度相关性本质上是由散体颗粒的冲击破碎引起的,加载速度越大,颗粒破坏程度越大,进而应力波的衰减程度越大;而在颗粒不破碎的情况下,加载速度增大时,衰减程度变化不明显。
-
利用离散元颗粒流软件PFC3D建立了散体颗粒模型,通过侧限压缩实验验证了模型的准确性。基于该离散元模型,研究了散体颗粒在冲击载荷下应力波的传播规律和主要影响因素。主要结论如下。
(1)应力波在散体颗粒中传播时,其峰值应力呈指数型衰减,且脉宽逐渐弥散。衰减因子β可定量描述应力波的衰减程度,β值越大,应力波的衰减越显著。
(2)散体颗粒在冲击加载后,根据颗粒的微观破坏形貌,可分为破碎区和微/未破碎区,两者的峰值应力也均呈指数衰减。破碎区的衰减程度远大于微/未破碎区,颗粒的破碎是影响散体颗粒中应力波衰减的一个重要因素。
(3)加载脉宽是影响应力波在散体颗粒中衰减的另一个主要因素,加载脉宽越小,应力波峰值衰减越大。应力波在散体颗粒中的传播初期,应力峰值显著衰减,是颗粒破碎和短脉宽高频波的共同作用。
(4)应力波衰减与加载速度显著相关,应力波峰值应力衰减随加载速度的增大而增大。但实际上却是颗粒破碎在宏观上的表象,加载速度越大,颗粒破坏程度越大,进而应力波的衰减程度越大。如果颗粒不发生破碎,加载速度增大时,应力波在散体颗粒中的衰减变化不明显。
应力波在散体颗粒中的传播规律
Attenuation Law of Stress Wave in Granular Particles
-
摘要: 脆性散体颗粒在受到冲击加载时,会对应力波传播产生显著的衰减作用。基于离散元颗粒流软件PFC3D建立了散体颗粒模型,通过不同加载速率下的数值模拟,研究应力波在散体颗粒中的传播规律和影响因素。结果表明:在冲击加载下,散体颗粒中传播的应力波峰值整体呈指数型衰减,随着传播距离的增大,应力波衰减程度逐渐减小,颗粒破碎程度也逐渐减小;应力波在散体颗粒中传播时会发生显著的波形弥散,并且应力波波长越短,传播过程中的衰减越大;应力波衰减的率相关性本质上是由散体颗粒的冲击破碎引起的,加载速度越大,颗粒破坏程度越大,应力波的衰减程度越大,而在颗粒不破碎的情况下,加载速度增大时,应力波的衰减程度变化不明显。Abstract: Brittle loose particles exhibit very complex mechanical behavior during the crushing process and have a significant attenuation effect on the stress wave propagation. In order to explore the attenuation law, this paper builds a brittle loose particle model based on the discrete element software PFC3D, and studies the attenuation of stress wave on the microscopic scale. The results show that: under shock loading, the peak value of the stress wave propagating in the granular particles decays exponentially. As the propagation distance increases, the degree of stress wave attenuation gradually decreases, and the degree of particle fragmentation also decreases. Stress wave propagation in granular particles will cause significant wave dispersion, and the shorter the wavelength of the stress wave, the greater the attenuation during propagation. The rate dependence of the stress wave attenuation is essentially caused by the impact fragmentation of the granular particles. The faster the loading speed, the greater degree of particle damage, and the greater the attenuation of stress wave. When the particle is not broken, the degree of attenuation does not change significantly with the increase of the loading velocity.
-
Key words:
- brittle particles /
- discrete element model /
- degree of fracture /
- stress wave /
- attenuation law
-
表 1 Flat-Joint模型微观参数
Table 1. Micro-parameters of the Flat-Joint model
Model Emod/GPa Fj-Kn/GPa Fj-Ks/GPa Fj-Coh/GPa Fj-Ten/GPa Damp Flat-Joint 50 5.25 × 105 9.37 × 105 0.2 0.05 0.15 -
[1] JAEGER H M, NAGEL S R. Granular solids, liquids, and gases [J]. Reviews of Modern Physics, 1996, 68(4): 1259–1273. doi: 10.1103/RevModPhys.68.1259 [2] 张家铭, 汪稔, 张阳明, 等. 土体颗粒破碎研究进展 [J]. 岩土力学, 2003(S2): 661–665.
ZHANG J M, WANG R, ZHANG Y M, et al. Advance in studies of soil grain crush [J]. Rock and Soil Mechanics, 2003(S2): 661–665.[3] 王礼立. 应力波基础 [M]. 北京: 国防工业出版社, 1985: 45–47.
WANG L L. Foundation of stress waves [M]. Beijing: National Defense Industry Press, 1985: 45–47.[4] 王肖钧. 分层人防工程防护结构对地下爆炸波的阻尼、耗散和导向作用研究 [R]. 合肥: 中国科学技术大学近代力学系, 2008.
WANG X J. Research on the damping, dissipation and guiding effect of layered civil air defense engineering structure on underground explosion waves [R]. Hefei: University of Science and Technology of China, 2008.[5] LI J C, MA G W. Experimental study of stress wave propagation across a filled rock joint [J]. International Journal of Rock Mechanics and Mining Sciences, 2009, 46(3): 471–478. doi: 10.1016/j.ijrmms.2008.11.006 [6] 季顺迎, 李鹏飞, 陈晓东. 冲击荷载下颗粒物质缓冲性能的试验研究 [J]. 物理学报, 2012, 61(18): 184703. doi: 10.7498/aps.61.184703
JI S Y, LI P F, CHEN X D. Experiments on shock-absorbing capacity of granular matter under impact load [J]. Acta Physica Sinica, 2012, 61(18): 184703. doi: 10.7498/aps.61.184703[7] 赵跃堂, 郑守军, 郑大亮, 等. 爆炸波在饱和土介质中传播时压力变化规律的试验研究 [J]. 防灾减灾工程学报, 2004, 24(1): 60–65.
ZHAO Y T, ZHENG S J, ZHENG D L, et al. Experimental investigation on pressure variation during explosion wave propagation in saturated soils [J]. Journal of Disaster Prevention and Mitigation Engineering, 2004, 24(1): 60–65.[8] 赵凯. 分层防护层对爆炸波的衰减和弥散作用研究 [D]. 合肥: 中国科学技术大学, 2007: 25–40.
ZHAO K. The attenuation and dispersion effects on explosive wave of layered protective engineering [D]. Hefei: University of Science and Technology of China, 2007: 25–40.[9] 魏久淇, 吕亚茹, 刘国权, 等. 钙质砂一维冲击响应及吸能特性试验 [J]. 岩土力学, 2019, 40(1): 191–198.
WEI J Q, LÜ Y R, LIU G Q, et al. One-dimensional impact responses and energy absorption of calcareous sand [J]. Rock and Soil Mechanics, 2019, 40(1): 191–198.[10] YU X, CHEN L, FANG Q, et al. Determination of attenuation effects of coral sand on the propagation of impact-induced stress wave [J]. International Journal of Impact Engineering, 2019, 125(1): 63–82. [11] 祁原, 黄俊杰, 陈明祥. 可破碎颗粒体在动力载荷下的耗能特性 [J]. 力学学报, 2015, 47(2): 254–259.
QI Y, HUANG J J, CHEN M X. Energy dissipation characteristics of crushable granules under dynamic excitations [J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(2): 254–259.[12] 郑文. 颗粒物质体系复杂动力学行为研究 [D]. 合肥: 中国科学技术大学, 2013: 17–36.
ZHENG W. A study on complex dynamic properties of granular material [D]. Hefei: University of Science and Technology of China, 2013: 17–36.[13] 黄俊宇. 冲击载荷下脆性颗粒材料多尺度变形破碎特性研究 [D]. 合肥: 中国科学技术大学, 2016: 19–40.
HUANG J Y. Dynamic multiscale deformation behavior and particle-breakage properties of granular materials subjected to impact loading [D]. Hefei: University of Science and Technology of China, 2016: 19–40.[14] 王刚, 叶沁果, 查京京. 珊瑚礁砂砾料力学行为与颗粒破碎的试验研究 [J]. 岩土工程学报, 2018, 40(5): 802–810. doi: 10.11779/CJGE201805004
WANG G, YE Q G, ZHA J J. Experimental study on mechanical behavior and particle crushing of coral sand-gravel fill [J]. Chinese Journal of Geotechnical Engineering, 2018, 40(5): 802–810. doi: 10.11779/CJGE201805004[15] 张科芬, 张升, 滕继东, 等. 颗粒破碎的三维离散元模拟研究 [J]. 岩土力学, 2017, 38(7): 2119–2127.
ZHANG K F, ZHANG S, TENG J D, et al. 3D numerical simulation of particle breaking using discrete element [J]. Rock and Soil Mechanics, 2017, 38(7): 2119–2127.[16] 张家铭, 汪稔, 石祥锋, 等. 侧限条件下钙质砂压缩和破碎特性试验研究 [J]. 岩石力学与工程学报, 2005, 24(18): 3327–3331. doi: 10.3321/j.issn:1000-6915.2005.18.022
ZHANG J M, WANG R, SHI X F, et al. Compression and crushing behavior of calcareous sand under confined compression [J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(18): 3327–3331. doi: 10.3321/j.issn:1000-6915.2005.18.022[17] POTYONDY D O, CUNDALL P A. A bonded-particle model for rock [J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41(8): 1329–1364. doi: 10.1016/j.ijrmms.2004.09.011 [18] ZHANG S, ZHANG F. A thermo-elasto-viscoplastic model for soft sedimentary rock [J]. Soils and Foundations, 2009, 49(4): 583–595. doi: 10.3208/sandf.49.583 [19] WNG J F, LI Y L, GAO Y B, et al. Experimental study on structural properties influencing on compressibility of soft clay [J]. Advanced Materials Research, 2011, 261/262/263: 767–1772. [20] YAN W M, LI X S. Mechanical response of a medium-fine-grained decomposed granite in Hong Kong [J]. Engineering Geology, 2012, 129/130: 1–8. [21] HUANG J Y, XU S L, HU S S. Influence of particle breakage on the dynamic compression responses of brittle granular materials [J]. Mechanics of Materials, 2014, 68: 15–28. doi: 10.1016/j.mechmat.2013.08.002 -