Numerical Study on Wave Effect of the Frictional Interface
-
摘要: 界面摩擦是一种普遍的自然现象。基于摩擦的界面微接触断裂机制,采用线弹性本构关系和D-P破坏准则,建立了包含三角形微凸起的二维界面摩擦模型,采用有限元分析对入射波和摩擦界面的相互作用进行数值模拟。结果表明:在主动加载的微过程中,界面存在显著的应力波动及精细结构特征,波阵面在界面近区域内的演化具有对称扩散性,应力扰动作用于界面微凸起可诱发其断裂,从而以断裂面为中心形成纵波、横波和界面波结构。一个有趣的现象是,在加载的瞬间,界面几乎同步产生了微应力扰动,以纵波形式向基体内传播,更多比较算例和分析证实该扰动产生的物理机制同作用在界面的整体重力微调整有关。该工作揭示了摩擦早期的界面波动效应及其微断裂机制,有望为地震预测提供新的有效途径,从而实现将地震预测时间提前。Abstract: Interface friction is a common natural phenomenon. Based on the micro-contact fracture mechanism of friction, a two-dimensional interface friction model including a triangular micro bulge is established with linear elastic constitutive relationship and D-P failure criterion. The early dynamic behavior of the interface under transient loading is numerically calculated and analyzed by the finite element simulation method. The research shows that in the micro process of loading, there exist significant stress fluctuations and fine structure characteristics at frictional interfaces. The evolution of the wavefront in the near region of the interface has symmetrical diffusion. The interaction of the incoming stress disturbance and the micro bulge will induce the fracture of the bulge, resulting in a three-wave profile centered on the fracture surface: longitudinal wave, transverse wave, and interface wave. A new interesting phenomenon is that at the moment of loading, a micro stress disturbance is generated synchronously from the interface and propagates to the substrate in the form of longitudinal waves. More comparative examples and analysis show that the mechanism of this disturbance is related to the overall gravity micro-adjustment acting on the interface. This work reveals the early wave effect of interface friction and its micro fracture mechanism, which is expected to provide an effective way for earthquake prediction and to advance the earthquake prediction time.
-
Key words:
- wave mechanics /
- interface friction /
- simulation /
- fracture /
- earthquake prediction
-
界面摩擦是材料的固有属性之一[1],普遍存在于各种材料和结构中。摩擦界面的宏观运动遵循经典摩擦定律。然而,若探究摩擦的物理机制,需在微细层面研究界面接触的结构特征和界面运动的微过程。研究界面静/动转化的传统方法一般通过静摩擦系数(fs)和动摩擦系数(fd)的变化来描述,然而这种描述过于粗糙,无法刻画界面运动状态转化的细节和机制。事实上,界面由静到动变化的瞬间,界面存在显著的波动现象和波的精细结构变化[2-4]。弄清这些波动规律,探索波动的起源和物理机制,将有助于深入理解摩擦的微观本质,拓展对界面摩擦动力学行为的认知,从而为工程中构件的摩擦设计提供理论支撑,为自然界各种摩擦运动(如地震等)的预判提供新的波动方法。
目前,对界面摩擦行为已有较多研究,大多集中在界面的摩擦特征,研究主要包括两个方向:(1) 探索各种物理过程中不同尺度下界面摩擦系数、摩擦运动的规律[5-11];(2) 关注特定尺度界面的摩擦运动对更大尺度物理过程的影响,如界面粗糙度对撞击或侵彻的影响[12-17]、纤维/基体界面强度对整体复合材料性能的影响[18]、颗粒接触摩擦对颗粒群剪切流动变形的影响[19-20]以及摩擦运动对地震的影响[21-23]等。上述研究工作涉及的时间尺度相对较长,贯穿整个相关物理过程。因此,不论是研究界面摩擦运动规律还是其影响,都没有过多关注界面摩擦运动早期由于摩擦运动导致的界面波动现象。事实上,不存在理想的准静态荷载,扰动的外荷载可诱发界面波动的发生和传播,这些波动信号携带了界面的实际微接触几何、物理特征和摩擦微运动信息。实验上实时检测该波动信号,理论上建立波结构和界面接触微结构的关系,从而分析其传播及演化规律,将有助于全过程理解界面的动态摩擦力学行为。
实验方面,直接观察摩擦界面上的波动现象是极具挑战性的,目前只有少数实验做到了这一点。Pyrak-Nolte 等[24]首先在断裂面上观察到一类新的弹性界面波,这类界面波包含两类,即快波(fast waves)和慢波(slow waves),其产生机制同界面的断裂刚度有关。Xia等[25]在有机玻璃摩擦界面滑动过程中,通过光学方法直接观测到界面上的亚瑞利波和超剪切波。Rubinstein等[2-3]也开展了有机玻璃准静态剪切实验,通过光学方法精细测量了界面微滑移过程中实际接触面的几何形态变化,结果表明,在界面由静到动的微转化瞬间,接触界面产生并传播了3种应力扰动,即快瑞利波、超剪切波和慢波。Ferrer等[26]实时测量了界面摩擦运动瞬时来自于界面的声辐射信号,认为这些信号幅值和数量的变化可以用来界定界面由静到动转化的特征时间节点。上述实验揭示了一个重要现象:在摩擦界面由静到动转化的瞬时,界面的运动和变形行为是个微动态过程,并伴随显著的波动特征。
理论方面,Braun等[27]建立了弹簧-滑块模型,采用接触刚度综合描述界面的接触特征,结果表明,在界面运动尚未转为整体摩擦滑移之前,由于界面的波动作用,界面初始未受扰动区已经转化为高度非均匀应力状态。Svetlizky等[28]提出了粗糙界面有限断裂的理论模型,认为界面精细波结构的产生与界面上大量实际微接触群被瞬时剪断事件有关,基于模型模拟计算的波结构与实验中的慢波[2-3]吻合较好,但与其他两种波有一定差异。Bartolomeo等[29]建立了界面波结构特征和界面运动转化的联系,认为弄清该波结构特征和传播规律将有助于理解界面摩擦运动的机制。Kammer等[30]建立了断裂能与界面摩擦强度的关系模型,其三维有限元模拟计算表明,界面存在亚瑞利波和超声波,这两类波的产生和转化取决于两种尺度的相互作用,即界面黏性区域宽度和非均匀区特征尺度。事实上,各种尺度的界面均是粗糙不平的,这种粗糙性在更小尺度上表现为界面的实际自然接触具有更精细的微尺度起伏性,这些起伏的结构特征和演化规律必然深刻影响界面的摩擦强度和动力学行为[1, 5]。界面摩擦运动诱发的界面波动本质上来源于界面实际接触中的几何结构变化,也必然蕴含了这些接触的特征信息。然而,困难在于难以对界面上波的精细结构进行定量描述和解释,从而建立界面动力学行为和波结构的联系。
本工作主要研究界面微接触断裂事件引起的波动效应,为此建立了简单界面接触模型,以模拟界面微滑动过程中微接触的断裂行为以及界面上相应的波结构特征。该模型将展示摩擦滑动过程中界面微凸起断裂的演变过程及其相应的波结构,以揭示波的精细结构特征和产生机理。
1. 平面摩擦滑动模型
采用ABAQUS软件建立了平面摩擦滑动模型,如图1所示。上基体(Part-1)、下基体(Part-2)的宽度均为10 mm,高度分别为40和10 mm。为表征界面的粗糙性,宏观上赋予界面摩擦系数,细观上在界面中心设计了边长为0.1 mm的等边三角形微接触凸起。考虑到摩擦界面的地震背景,摩擦界面的样本材料选择弹脆性玄武岩体,材料的破坏采用D-P准则描述,相关弹性参数和断裂参数见表1。应力脉冲荷载作用下,摩擦界面的动力学行为是一个微过程的波传播问题,为展现波传播的细节和精细结构,单元网格的划分非常小,尤其是在微凸起近区域,最小网格尺寸5 μm,时间步长10−9 s。
t=0时刻在上部滑块左侧施加瞬态荷载(σ*)10 MPa,下滑块底面为固定约束,其他边界设置为无反射类型。如图1所示,为更好地揭示波动过程和应力波的精细结构,选取两类单元:(1) 在界面上以微凸起为中心,上下、左右对称选取8组单元,形成4个区域,通过这些单元的应力历程展现界面上的波动效应和精细波结构;(2) 在微凸起的正上方,旨在揭示基体内的波结构。
表 1 计算材料参数Table 1. Material parameters for calculationDensity/
(kg·m−3)Elastic modulus/GPa Shear modulus/GPa P wave velocity/(m·s–1) S wave velocity/(m·s–1) Friction coefficient Internal friction angle/(°) 2300 62.8 24.1 5225 3237 0.1 44 Expansion angle/(°) Hardening coefficient Fracture strain Tensile
strength/MPaCohesion strength/MPa Shear stress ratio Absolute plastic strain 0 6.98 0.0075 3.5 8 0.33 0 2. 计算结果及分析
2.1 界面上的波结构特征
图2为上界面单元的3种应力波结构特征。从上到下依次为σ11、σ22和σ12,从左至右依次为Zone-1、Zone-2、Zone-3和Zone-4。σ11和σ22分别为图1中x和y方向的正应力,σ12为单元的剪应力。图2中σ11的分布演化表明,在3 μs微过程中,σ11的变化呈现出显著的3波结构特征。最先到达的波速度约为5200 m/s,幅值为10 MPa,显然该波来自于加载端面的平面纵波;后继两个波的平均速度约为5200和3000 m/s,可能是纵波和Rayleigh界面波。由于Zone-2和Zone-3、Zone-1和Zone-4均关于微凸起对称分布,从波的到达时间来看,这两个扰动来自于微凸起的断裂,即微凸起断裂瞬时,沿接触界面将传播纵波和界面波,速度约为0.93Cs。从σ22和σ12波结构来看,来自于加载端面的纵波同样会引起界面单元竖直方向应力分量σ22和剪应力σ12的变化,幅值强度分别为1.22和0.12 MPa,但二者跳跃变化的机制不同,σ22的变化源于泊松效应,而σ12的微扰动与界面的微摩擦相对滑动有关。另外,也可以看到,微接触凸起断裂诱发的纵波和界面波不会引起Part-1部分界面单元σ22应力分量变化,但Rayleigh波会导致该类界面单元应力σ12的扰动。值得注意的是,在σ22的波形结构中,左端面纵波未到之前已存在一个微扰动(图2和图3中σ22应力波动结构图中灰色矩形框),并且这些微扰动信号的幅值一致,时间同步(到达时约为0.3 μs)。对比相应的σ11和σ12的波结构发现,并不存在类似波形,说明该扰动并不是沿着界面传播的,而是起源于界面,进一步的物理机制需要结合微凸起的断裂过程来深入分析。
图3为下界面单元的3种应力波结构特征。在图3所示σ11的波动信号中,也呈现类似的三波结构:其中后到达的两个波同界面上方对应单元一致,均来自于界面微凸起断裂产生的纵波和界面波;与上部对称单元相比,最先到达的波的速度和到达时间是一致的,说明两个波同源,但其幅值差别较大,这是因为上部单元位于Part-1上,承受主动荷载,而下部单元在Part-2上,其σ11的变化主要在于界面的摩擦作用。对比图2和图3的波形结构,可以看出,除个别点外,界面上下对称单元的剪应力σ12和垂直应力σ22的扰动变化几乎是完全一致的,这源于扰动起源的一致性和所取单元位置的对称性。另外,在σ22的变化过程中,相同时刻主动加载纵波未到达之前也出现了相似的微扰信号。
综合图2和图3展示的界面波动信号可以看出,在施加载荷的瞬间,界面存在显著的波动现象和精细的波结构,这种精细结构特征源自于加载波和微凸起的相互作用造成的微接触断裂。
2.2 应力波与界面微缺陷相互作用过程
图4以应力云图的形式给出了加载弹性波的传播图像及与界面微凸起的相互作用过程。由于图4主要展示了波传播过程及波的精细结构特征,因此并没有绘出反映应力大小的标度柱状图,这样使得图4中波结构更简明清晰。t=0时,加载弹性波包含3个应力扰动,即σ11、σ22和σ12以平面纵波的形式向基体内传播。图4中,t=0.905 μs时,加载纵波传播至界面微凸起位置,并与之相互作用,由于应力集中引起微凸起断裂,从而以断裂点为中心,形成球面P波和S波向上、下基体内传播,同时沿着界面形成Rayleigh波,图4(b)中σ22的变化清晰展示了该过程。对于界面上的单元,只能感应到P波和Rayleigh波,因此在图2和图3中,在主动加载波后,σ11为双波结构(P波和Rayleigh波),而σ22和σ12只有Rayleigh波,即界面上的单元无法反映由于微凸起断裂而形成的S波。图4(b)中,自t=0.311 μs开始,可以清晰看到两个P波自界面分别向上、下基体中传播,如图2和图3中σ22信号中的灰色区域所示。剪应力σ12的扰动波阵面的演化如图4(c)所示。整体来看,剪应力σ12的扰动阵面是以界面为中心呈近似对称平面圆锥形,随着主平面纵波σ11阵面一起传播,在传播过程中,圆锥面逐渐向上、下基体扩大。
2.3 微接触的断裂过程及波动效应
如图4所示,在应力波演化过程中的一个重要事件是界面上微接触发生断裂,但图4并没有清晰地描绘该过程中的断裂过程以及相应波的精细结构。图5通过σ22应力云图给出了微凸起断裂的微过程。图5中,当t=1.272 μs时,由于应力集中,三角形微凸起的左下角首先开始起裂,接着裂纹向下部基体和沿凸起根部传播(t=1.562 μs和t=1.604 μs)。与此同时,微凸起的右下角位置也开始出现裂纹(t=1.660 μs),两个角裂纹沿着界面运动(t=1.724 μs),最终微凸起从根部剪断(t=1.740 μs)。在该过程中,主要以微凸起的左、右下角为中心形成新的应力扰动。
微凸起断裂形成新的次声波结构如图6所示。在t=1.780 μs时的应力云图中,可以清晰地观察到纵波、横波和界面波三波结构。通过图2可知,该界面波主要引起界面σ11和σ12的变化,其传播速度约为0.93Cs。显然,只有上、下基体内的单元能感应到纵波和横波扰动。为证实这一点,提取微接触的正上方单元的应力扰动变化信息,如图7所示。图7中,当断裂引起的球面纵波传播至捡取单元时,主要引起σ22和σ12的较大变化,σ11有微小扰动,而横波主要体现在σ12的变化上。因此,从σ12的扰动信号中可以清晰地观察到凸起断裂形成的双波结构,其波速分别为5225和3239 m/s。
3. 讨 论
有限元模拟结果表明,在施加载荷的瞬间,摩擦界面存在显著的波动效应。界面上波的精细结构主要包含3部分:主动加载脉冲、界面微凸起断裂引起的球面扰动和起源于界面的微扰动。基于传统的固体波动理论和特征线分析方法[31],可以在x-t平面上绘出前两种波的传播及结构特征,如图8所示。由于模拟对象处于近似一维应变状态,因此初始时刻的加载弹性波包含3个应力扰动,σ11、σ22和σ12独立向基体内传播。图9给出了t=0.315 μs时3个应力波动的云图。宏观上,3个应力扰动的前沿阵面传播速度是一致的,但由于摩擦界面效应,可以清晰地看到,界面近区域波前沿阵面的形状发生了显著变化。具体而言,σ11在波阵面的右下尖端有应力集中带并延伸至下基体,σ22在界面近区域形成关于界面对称的三角锥,剪应力σ12的波阵面整体上是平面纵波,但在界面近区域呈现出较大的对称圆锥,这种近似圆锥形的剪应力阵面同Rubino等[23]的实验结果是一致的,但其形成机制需要进一步深入研究。
图2和图3的σ22波结构中还存在一个精细的微小扰动。图4(b)所示的σ22应力云图表明,该扰动在t=0.300 μs开始在界面上出现,此时左端面的加载纵波传播了约1.56 mm,这对于介质的内部单元(尤其是近界面区域)来说,意味着波尚未传至时已产生微应力扰动。该扰动并没有沿界面传播,而是起源于界面且向上、下基体材料内传播。为进一步验证这一点,在微凸起正上方取一系列单元,其σ22扰动信号如图10所示,对灰色框所示的信号放大,从而清晰地观察到一个速度约为5220 m/s的微扰动信号自界面处向上传播,显然该波是起源于界面的平面纵波。
基于上述分析形成了初步认识,该波起源于界面,是平面纵波,但产生的物理机制仍不清楚。调整界面摩擦系数(μ=0,0.3,0.5)对该波的形成没有影响(图11(a)),推测该波可能是微凸起引起的,但改变微凸起的大小(三角形微接触边长a=0.1,0.5,1.0 mm)和形状甚至将微凸起去掉后,发现该波仍然存在(图11(b))。基于对该波到达时的分析,推测该波可能与上部基体的重力作用于界面有关,为此进一步设计了两个模型:一个是将界面去掉,上下为一个整体,此时该波不存在;另一个是将原模型顺时针旋转90°,即重力不再直接作用于界面,此时该波亦不存在。上述两种模型说明该波的形成可能与作用在界面上的重力调整有关。如图12(a)所示,根据惠更斯原理,在外部荷载作用下,界面上的每个点都是一个新的波源,产生球面扰动,这些球面波系的前沿波阵面形成一个包络面,该包络面的传播速度与球面P波一致,形式上表现为平面纵波的特征。图12(b)所示的数值模拟结果清楚地展现了界面上球面波系与材料内向上下对称方向传播的新纵波扰动的关系。但这里仍然有些基本问题有待商榷,如:界面上单元最初的扰动从哪里来,是否与重力作用有关系,该波出现的时间是否固定,同哪些因素有关。
该波的存在给了我们深刻的启示。目前,地震预测系统建立在地震P波和S波之间的速度差上,当地震发生时,基于本研究中的模型结果,在球面P波抵达地表探测系统之前,应有重力扰动引起的纵波存在,若能捡取、分离、明确该信号扰动,将有助于将地震预报的时间提前。为了验证其可行性,进行了宏观大尺度地震模拟,如图13所示。显然,先于地震P波(图13(b)中蓝色虚线部分)的纵波扰动(图13(b)中红色虚线部分)是存在的,这方面的实验和进一步的数值模拟研究将在后续工作中逐步开展。
4. 结 论
建立了简单的界面摩擦模型,通过数值模拟计算分析外荷载作用瞬时界面上的波动效应,得到如下主要结论:
(1) 在加载瞬时,界面运动的微过程中界面上存在清晰的波的精细结构;
(2) 界面上波的精细结构与界面微接触的断裂有关,断裂将产生纵波、横波和界面波;
(3) 在界面微凸起断裂之前,由于重力扰动的作用,自界面产生微扰动,该微扰动以平面纵波的形式向基体内传播。
客观上,界面的粗糙度和起伏度要复杂得多,但本研究通过一个简单的三角形微凸起模型揭示了一个重要规律,即界面的摩擦动力学行为与界面的粗糙起伏引起的断裂密切相关,这种断裂事件将以波动的形式在基体和界面上传播。当考虑界面的实际粗糙度时,相应波结构产生的变化方面的研究将在后续工作中继续开展。另外,应力扰动σ22和σ12的传播过程中波阵面形状的变化机制也需要进一步深入分析,给出合理的物理解释。
-
表 1 计算材料参数
Table 1. Material parameters for calculation
Density/
(kg·m−3)Elastic modulus/GPa Shear modulus/GPa P wave velocity/(m·s–1) S wave velocity/(m·s–1) Friction coefficient Internal friction angle/(°) 2300 62.8 24.1 5225 3237 0.1 44 Expansion angle/(°) Hardening coefficient Fracture strain Tensile
strength/MPaCohesion strength/MPa Shear stress ratio Absolute plastic strain 0 6.98 0.0075 3.5 8 0.33 0 -
[1] 许金泉. 界面力学 [M]. 北京: 科学出版社, 2006.XU J Q. The mechanics of interface [M]. Beijing: Science Press, 2006. [2] RUBINSTEIN S M, COHEN G, FINEBERG J. Detachment fronts and the onset of dynamic friction [J]. Nature, 2004, 430(7003): 1005–1009. doi: 10.1038/nature02830 [3] RUBINSTEIN S M, COHEN G, FINEBERG J. Dynamics of precursors to frictional sliding [J]. Physical Review Letters, 2007, 98(22): 226103. doi: 10.1103/PhysRevLett.98.226103 [4] BEN-DAVID O, COHEN G, FINEBERG J. The dynamics of the onset of frictional slip [J]. Science, 2010, 330(6001): 211–214. doi: 10.1126/science.1194777 [5] ZHU Y D, ZHENG Z J, ZHANG Y L, et al. Adhesion of elastic wavy surfaces: interface strengthening/weakening and mode transition mechanisms [J]. Journal of the Mechanics and Physics of Solids, 2021, 151: 104402. doi: 10.1016/j.jmps.2021.104402 [6] PERSSON B N J. Sliding friction: physical principles and applications [M]. 2nd ed. Berlin: Springer, 2000. [7] TA W R, QIU S M, WANG Y L, et al. Volumetric contact theory to electrical contact between random rough surfaces [J]. Tribology International, 2021, 160: 107007. doi: 10.1016/j.triboint.2021.107007 [8] GERDE E, MARDER M. Friction and fracture [J]. Nature, 2001, 413(6853): 285–288. doi: 10.1038/35095018 [9] BAUMBERGER T, BERTHOUD P, CAROLI C. Physical analysis of the state- and rate-dependent friction law. Ⅱ. dynamic friction [J]. Physical Review B, 1999, 60(6): 3928–3939. doi: 10.1103/PhysRevB.60.3928 [10] BRAUN O M, MANINI N, TOSATTI E. Size scaling of static friction [J]. Physical Review Letters, 2013, 110(8): 085503. doi: 10.1103/PhysRevLett.110.085503 [11] BARRAS F, AGHABABAEI R, MOLINARI J F. Onset of sliding across scales: how the contact topography impacts frictional strength [J]. Physical Review Materials, 2021, 5(2): 023605. doi: 10.1103/PHYSREVMATERIALS.5.023605 [12] SHAO R L, WAHLE M, ZIMMERMANN M. A model for the dynamic friction behaviour of rubber-like materials [J]. Tribology International, 2021, 164: 107220. doi: 10.1016/j.triboint.2021.107220 [13] 张磊, 王文帅, 苗春贺, 等. 花岗岩粗糙表面动摩擦形态演化 [J]. 高压物理学报, 2021, 35(3): 031201. doi: 10.11858/gywlxb.20200640ZHANG L, WANG W S, MIAO C H, et al. Rough surface morphology of granite subjected to dynamic friction [J]. Chinese Journal of High Pressure Physics, 2021, 35(3): 031201. doi: 10.11858/gywlxb.20200640 [14] BERMAN N, COHEN G, FINEBERG J. Dynamics and properties of the cohesive zone in rapid fracture and friction [J]. Physical Review Letters, 2020, 125(12): 125503. doi: 10.1103/PhysRevLett.125.125503 [15] WANG P F, JIANG H B, XU S L, et al. Dynamic plastic instability of ring-shaped aluminum alloy with different interface behaviors [J]. International Journal of Impact Engineering, 2021, 155: 103898. doi: 10.1016/j.ijimpeng.2021.103898 [16] 赵剑衡, 孙承纬, 段祝平, 等. 玻璃样品表面对失效波萌生的影响 [J]. 力学学报, 2001, 33(6): 834–838. doi: 10.3321/j.issn:0459-1879.2001.06.014ZHAO J H, SUN C W, DUAN Z P, et al. Effect of impacted surface of K9 glass sample on formation of failure wave [J]. Acta Mechanica Sinica, 2001, 33(6): 834–838. doi: 10.3321/j.issn:0459-1879.2001.06.014 [17] 刘均伟, 张先锋, 刘闯, 等. 考虑摩擦因数变化的弹体高速侵彻混凝土质量侵蚀模型研究 [J]. 爆炸与冲击, 2021, 41(8): 083301. doi: 10.11883/bzycj-2020-0250LIU J W, ZHANG X F, LIU C, et al. Study on mass erosion model of projectile penetrating concrete at high speed considering variation of friction coefficient [J]. Explosion and Shock Waves, 2021, 41(8): 083301. doi: 10.11883/bzycj-2020-0250 [18] POCHIRAJU K V, TANDON G P, PAGANO N J. Analyses of single fiber pushout considering interfacial friction and adhesion [J]. Journal of the Mechanics and Physics of Solids, 2001, 49(10): 2307–2338. doi: 10.1016/S0022-5096(01)00045-X [19] 王蕉, 楚锡华. 冲击载荷下颗粒材料临边界区域的波动行为及变形特征分析 [J]. 力学学报, 2021, 53(9): 2395–2403. doi: 10.6052/0459-1879-21-242WANG J, CHU X H. Analysis of wave behavior and deformation characteristics of granular materials in pro-border zone under impact load [J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(9): 2395–2403. doi: 10.6052/0459-1879-21-242 [20] ZHENG W, ZHANG S Y, XU N. Jamming of packings of frictionless particles with and without shear [J]. Chinese Physics B, 2018, 27(6): 066102. doi: 10.1088/1674-1056/27/6/066102 [21] SCHOLZ C H. Earthquakes and friction laws [J]. Nature, 1998, 391(6662): 37–42. doi: 10.1038/34097 [22] KANAMORI H, ANDERSON D L, HEATON T H. Frictional melting during the rupture of the 1994 Bolivian earthquake [J]. Science, 1998, 279(5352): 839–842. doi: 10.1126/science.279.5352.839 [23] RUBINO V, ROSAKIS A J, LAPUSTA N. Understanding dynamic friction through spontaneously evolving laboratory earthquakes [J]. Nature Communications, 2017, 8: 15991. doi: 10.1038/ncomms15991 [24] PYRAK-NOLTE L J, XU J P, HALEY G M. Elastic interface waves propagating in a fracture [J]. Physical Review Letters, 1992, 68(24): 3650–3653. doi: 10.1103/PhysRevLett.68.3650 [25] XIA K W, ROSAKIS A J, KANAMORI H. Laboratory earthquakes: the sub-Rayleigh-to-supershear rupture transition [J]. Science, 2004, 303(5665): 1859–1861. doi: 10.1126/science.1094022 [26] FERRER C, SALAS F, PASCUAL M, et al. Discrete acoustic emission waves during stick-slip friction between steel samples [J]. Tribology International, 2010, 43(1/2): 1–6. doi: 10.1016/j.triboint.2009.02.009 [27] BRAUN O M, BAREL I, URBAKH M. Dynamics of transition from static to kinetic friction [J]. Physical Review Letters, 2009, 103(19): 194301. doi: 10.1103/PhysRevLett.103.194301 [28] SVETLIZKY I, FINEBERG J. Classical shear cracks drive the onset of dry frictional motion [J]. Nature, 2014, 509(7499): 205–208. doi: 10.1038/nature13202 [29] DI BARTOLOMEO M, MASSI F, BAILLET L, et al. Wave and rupture propagation at frictional bimaterial sliding interfaces: from local to global dynamics, from stick-slip to continuous sliding [J]. Tribology International, 2012, 52: 117–131. doi: 10.1016/j.triboint.2012.03.008 [30] KAMMER D S, MUÑOZ D P, MOLINARI J F. Length scale of interface heterogeneities selects propagation mechanism of frictional slip fronts [J]. Journal of the Mechanics and Physics of Solids, 2016, 88: 23–34. doi: 10.1016/j.jmps.2015.12.014 [31] 李永池. 波动力学 [M]. 合肥: 中国科学技术大学出版社, 2015.LI Y C. Wave mechanics [M]. Hefei: University of Science and Technology of China Press, 2015. -