Static friction coefficient depends on the external pressure and block shape due to precursor slip

Amontons’ law states that the maximum static friction force on a solid object is proportional to the loading force and is independent of the apparent contact area. This law indicates that the static friction coefficient does not depend on the external pressure or object shape. Here, we numerically investigate the sliding motion of a 3D viscoelastic block on a rigid substrate using the finite element method (FEM). The macroscopic static friction coefficient decreases with an increase in the external pressure, length, or width of the object, which contradicts Amontons’ law. Precursor slip occurs in the 2D interface between the block and substrate before bulk sliding. The decrease in the macroscopic static friction coefficient is scaled by the critical area of the precursor slip. A theoretical analysis of the simplified models reveals that bulk sliding results from the instability of the quasi-static precursor slip caused by velocity-weakening local friction. We also show that the critical slip area determines the macroscopic static friction coefficient, which explains the results of the FEM simulation.


Introduction
A friction force prevents the relative sliding motion between two objects in contact.Friction plays a crucial role in various situations, such as the contact surface between the ground and the sole of a shoe, brakes and bearings in machines, and tectonic plates that cause earthquakes.Many studies on friction have been conducted, but the elucidation of the fundamental mechanism of friction is essential for science and technology [1][2][3][4][5][6][7] .
Amontons' law states that the maximum static friction force on a solid object is independent of the apparent contact area and proportional to the load [1][2][3][4][5][6][7] .This law has been taught in high school physics textbooks and is believed to hold true for diverse systems.When the friction force obeys Amontons' law, the friction coefficient, which is the ratio of the friction force to the loading force, does not depend on the pressure, size, or object shape.On a rough frictional interface with numerous asperities, only a tiny fraction of the surfaces forms junctions, the so-called real contact points.Amontons' law is explained by the proportionality of the total area of real contact points to the loading force [1][2][3][4][5][6][7][8][9] .
The above explanation for the origin of Amontons' law implicitly assumes uniformity of the stress field.Therefore, Amontons' law is not expected to hold if a macroscopic deformation exists.In fact, recent numerical studies have reported the breakdown of Amontons' law in macroscopic viscoelastic objects 10,11 , revealing that it is related to local quasi-static precursor slips before the onset of bulk sliding owing to non-uniform deformation [10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28] .The relationship between precursor slips and the breakdown of Amontons' law has been confirmed previously in an experiment with an acrylic glass block 12 .However, previous studies have only investigated systems with a 1D frictional interface.Friction usually occurs in 2D interfaces of 3D objects.However, it is not clear whether the results in previous studies apply to more realistic 3D systems.
In this study, we numerically investigate the sliding motion of a 3D viscoelastic object on a rigid substrate using the finite element method (FEM).The macroscopic static friction coefficient decreases with an increase in the pressure or size of the object.The precursor slip propagates in a 2D frictional interface.Bulk sliding occurs when the area of the precursor slip reaches a critical value, which determines the macroscopic static friction coefficient.An analysis of the simplified models reveals that the instability of the precursor slip leads to bulk sliding.

3D FEM simulation
We numerically investigate a viscoelastic block on a rigid substrate with width W , length L, and height H along the x-, y-, and z-axes, respectively, as shown in Fig. 1 (see Methods for details).The area of the frictional interface is denoted by A 0 = LW .The density, Young's modulus, and Poisson's ratio of the block are denoted by ρ, E, and ν, respectively.The dissipation in the block is characterized by two viscosity coefficients: η 1 and η 2 .We assume that Amontons' law holds locally at the interface arXiv:2211.01066v2[cond-mat.mtrl-sci]20 Jan 2023 between the block and the rigid substrate (z = 0), and the magnitude of the local frictional stress, σ (fric) (x, y) in the interface is locally determined as where p(x, y) is the bottom pressure, and µ(v) is the friction coefficient, which depends on the magnitude of the local slip velocity v(x, y) when v(x, y) = 0 29 .Here, µ(v) is characterized by the characteristic velocity of velocity-weakening friction v c and the local static and kinetic friction coefficients denoted by µ S and µ K (see Methods).The rigid rod quasi-statically pushes the center of the side surface along the y direction.The macroscopic friction force F T is measured as the force on the rigid rod in the y direction.The loading force applied to the top of the block is given by F N = P ext A 0 with the external pressure to the top surface P ext .

Rigid rod
Viscoelastic block

Rigid substrate
Friction force   4) and ( 6).(c) Macroscopic static friction coefficient µ M against P ext for various W /H values with L/H = 1.The thin solid lines represent the analytical results with α B = 0.2 given by equations ( 6) and (11).The dotted and dashed lines represent µ S and µ K , respectively.The ratio F T /F N is plotted against the displacement of the rigid rod U for L/H = 1, W /H = 2, and P ext /E = 0.006 in Fig. 2 (a).First, F T /F N increases linearly with U.After obtaining a maximum value lower than µ S , F T /F N rapidly decreases to a value close to µ K .This rapid drop is associated with bulk sliding.The significant drop after the linear increase periodically repeats itself.This periodic behavior corresponds to the stick-slip motion of the object.The maximum value of F T /F N represents the macroscopic static friction coefficient, µ M .Figures 2 (b) and (c) display the macroscopic static friction coefficient µ M against pressure P ext for various L/H and W /H values, respectively.The magnitude of µ M decreases with increasing P ext , which is qualitatively consistent with the results for a system with a 1D friction interface in ref. 10.The previous study reported the size dependence of µ M while maintaining the aspect ratio L/H = 2 10 , whereas Figs. 2 (b) and (c) demonstrate that the friction coefficient µ M also decreases with increasing aspect ratios L/H and W /H.These results indicate that Amontons' law breaks down in systems with 2D interfaces.
Figure 3 (a) shows the spatial distribution of the slip region with nonzero slip velocity in the frictional interface at z = 0 for U = U 1 ,U 2 ,U 3 , and U 4 shown in Fig. 2 (a).Here, we choose U 1 /L = 50 × 10 −3 , U 2 /L = 56 × 10 −3 , U 3 /L = 61.38 × 10 −3 , and U 4 /L = 62.71 × 10 −3 , which corresponds to the stationary stick-slip region.See Methods for the definition of the slip region.In Fig. 3 (a), the local precursor slip starts from the region under the rigid rod for U = U 1 .As U increases (U 2 and U 3 ), the region expands gradually.After U = U 3 , the entire area slips with v > v c , resulting in bulk sliding.Note that the slip occurs almost along the y direction.Figure 3 (b) shows the area of precursor slip A normalized by the area of frictional interface A 0 against displacement U. First, the area of the precursor slip increases gradually with displacement U. When the area A reaches the critical area A c just before bulk sliding (dotted line), the propagation speed of the area suddenly increases.Owing to rapid propagation, A reaches A 0 and then returns to 0. We demonstrate the normalized critical area A c /A 0 against pressure P ext in Figs. 3 (c) and (d) for various L/H values with W /H = 1 and for various W /H values with L/H = 1, respectively.The normalized critical area A c /A 0 decreases as P ext , L/H, or W /H increases.This decrease is similar to that of µ M in Figs. 2 (b) and (c), respectively.
In Fig. 4, we present the macroscopic friction coefficient µ M against the normalized critical area A c /A 0 for various L/H and W /H values.The macroscopic friction coefficient µ M for different L/H and W /H values approximately collapses onto a master curve, which indicates a linear increase in µ M with A c /A 0 .The minimum value close to A c /A 0 = 0 is almost equal to µ K , whereas the maximum value at A c /A 0 = 1 is equal to µ S .Figure 5 shows the spatial distribution of the ratio σ (fric) /p in the frictional interface for L/H = 1, W /H = 2, and P ext /E = 0.006 at U = U 1 ,U 2 ,U 3 , and U 4 .It should be noted that the direction of the frictional stress is almost opposite to the driving direction, that is, the y direction.In the no-slip region, the local static friction can take any value for 0 < σ (fric) /p < µ S .Before the onset of precursor slip, that is, just after bulk sliding, σ (fric) /p takes a value almost equal to µ K , the local kinetic friction coefficient, in the entire interface, as explained below.At U = U 1 , σ (fric) /p reaches the local static friction coefficient, µ S , near the rigid rod at (x/H, y/H) = (1, 0).As the displacement U increases to U 2 and U 3 , the area with σ (fric) /p µ S gradually increases.The region of σ (fric) /p µ S coincides with the local precursor slip region in Fig. 3 (a).Except for the slip region, σ (fric) /p remains approximately at µ K .Immediately after U 3 , bulk sliding with v > v c occurs, and the fast slip leads to σ (fric) /p = µ K at U 4 .Bulk sliding rapidly decelerates, and the slip velocity v decreases to 0, when σ (fric) /p increases to µ S in the frictional interface.However, the internal deformation is not able to follow the rapid change, and the ratio of static frictional stress to bottom pressure finally returns to σ (fric) /p µ K after bulk sliding.Consequently, σ (fric) /p is almost equal to µ K after bulk sliding.The macroscopic static friction coefficient µ M is approximately expressed by the average of σ (fric) /p over the entire frictional interface at U 3 immediately before bulk sliding.This result explains the dependence of µ M on A c /A 0 shown in Fig. 4, where µ M approaches µ S for A c /A 0 = 1.

Analysis based on simplified models
To theoretically analyze the numerical results, we employ two simplified models, which explain the dependence of µ M on L/H and W /H (see Supplementary Note online for details).

Model for large L/H
To discuss the behavior of increasing L/H while maintaining W /H = 1, we employ a 1D effective model, as shown in Fig. 6 (a).
The slip region propagates along the y direction, as shown in the Supplementary Note and Supplementary Video S1.Therefore, in this model, the degrees of freedom in the z and x directions are neglected by assuming W /H L/H, and the deformation is characterized only by the y-dependent displacement in the y direction, u y (y,t), at the interface z = 0. We also assume a uniform bottom pressure P ext .The equation of motion is given by where üy and uy are the second-and first-order time derivatives of the displacement, respectively.Here, α A represents the effect of the block thickness and is treated as a fitting parameter.The normal stress σ yy is given by with the elastic constant in the plane stress state by considering the block as a thin plate (see Methods).The quasi-static solution u (a) y (y) of equation ( 3) with üy = uy = 0 is obtained analytically, where the precursor slip area A increases with U (see Supplementary Note).A linear stability analysis reveals that the quasi-static solution becomes unstable, and bulk sliding occurs when A reaches the critical area A c owing to the competition between velocity-weakening friction and viscosity.The critical area A c satisfies (see Supplementary Note).For A c /A 0 1, equation (4) yields This equation indicates that the normalized critical area A c /A 0 decreases as L/H or P ext increases, which is consistent with the FEM results shown in Fig. 3 (c).We plot A c /A 0 obtained from equation ( 4) as thin solid lines in Fig. 3 (c) by choosing the fitting parameter α A = 0.2 to match the results of the FEM simulations.The analytical results semi-quantitatively reproduce the numerical results except for L/H = 1.
The quasi-static solution u This is consistent with the FEM simulations, as shown by the solid line in Fig. 4. For A c /A 0 1, substituting equation ( 5) into this equation, we obtain

5/11
This equation indicates that the macroscopic static friction coefficient µ M decreases as P ext or L/H increases.We plot µ M given by equations ( 4) and ( 6) as thin solid lines in Fig. 2 (b), which semi-quantitatively reproduces the results of the FEM simulations except for L/H = 1.
In a previous study 10 , µ M is obtained analytically as for A c /A 0 1 in a system with small L/H.Here, α is the fitting parameter.The power-law exponents in equation ( 8) for the dependence on P ext and L/H differ from those in equation (7).The present model assumes L/H 1, which results in a uniform bottom pressure, as shown in the Supplementary Note.For a small L/H, the bottom pressure increases along the driving direction owing to the torque effect 10,21 , and the analytical results deviate from those of FEM due to the non-uniform pressure as shown in Figs. 2 (b) and 3 (c), which leads to different exponents from those in the present study.

Model for large W /H
To discuss the behavior of increasing W /H while maintaining L/H = 1, we employ a 1D effective model, as shown in Fig. 6 (b).For W /H 1, the slip region propagates along the x direction, as shown in Supplementary Note and Supplementary Video S2.Hence, in this model, we characterize the deformation only by the x-dependent displacement in the y direction, u y (x,t), at the interface z = 0 by assuming L/H W /H. We also assume the rod size is sufficiently small and negligible.The equation of motion is given by Here, α B represents the effect of the block thickness and is treated as a fitting parameter.The shear stress σ xy is given by with the elastic constant E 2 = E/ {2(1 + ν)} and the viscous constant η 1 /2 (see Methods).
The quasi-static solution u (a) y (x) is also obtained analytically, where the precursor slip area A increases with the value of U (see Supplementary Note).The linear stability analysis reveals that the precursor slip becomes unstable, and bulk sliding occurs when A reaches the critical area A c satisfying For A c /A 0 1, this equation yields The power-law exponents for the pressure and aspect ratio are the same as those in equation ( 5).This equation indicates that A c /A 0 decreases as P ext or W /H increases.We plot A c /A 0 given by equation ( 11) as thin solid lines in Fig. 3 (d), which semi-quantitatively reproduces the results of the FEM analysis by choosing α B = 0.2 except for W /H = 1.For small W /H, the size of the rod and the y-dependence of the displacement become relevant, which leads to the deviation between the numerical and theoretical results.The macroscopic static friction coefficient µ M is given by equation (6).For A c /A 0 1, substituting equation ( 12) into equation ( 6), we obtain The macroscopic static friction coefficient µ M decreases as P ext or W /H increases.The thin solid lines shown in Fig. 2 (c) are given by equations ( 6) and ( 11), and they semi-quantitatively reproduce the results of the FEM simulations except for W /H = 1.

Discussion
In this study, we numerically investigate the sliding motion of a 3D viscoelastic object using the FEM.The critical area of the precursor slip and macroscopic static friction coefficient decrease with an increase in the external pressure, length, or width of the object.The analysis based on the simplified models reveals that the stability condition determines the critical area of the precursor slip owing to the competition between the velocity-weakening friction and viscosity.The analysis explains the dependence of macroscopic static friction in the FEM simulations.
In a previous study 10 , the aspect ratio of the system is fixed at L/H = 2 to investigate the size and load dependences of the precursor slip and the breakdown of Amontons' law.For L/H = 2, the nonuniformity of the bottom pressure is remarkable, which is considered to be the origin of the precursor slip and the breakdown of Amontons' law.However, the present results with various aspect ratios show that the nonuniformity of shear stress also causes these behaviors without non-uniform pressure.Although the model considered in the previous study reproduces the results of systems with a smaller L/H better, the simplified model in this study is more appropriate for systems with a large L/H (see Supplementary Note).
The parameters for the FEM simulations employed here are those of a virtual material, and different from those of poly methyl methacrylate (PMMA) employed in experiments 12,15 .We choose them to compare our results with the 2D simulations of previous studies 10 and to reduce the computational load (see Methods).It also should be noted that the driving rod employed in experiments is hard but has finite stiffness, which is different from the rigid rod used in this study.The effect of the finite stiffness of the driving rod is considered to be small because it is taken into account as a deformation of the viscoelastic block around the driving point.In addition, we have ignored the aging effect 7 in the local friction model because a previous experiment using PMMA 15 indicates that the time scale of the aging is larger than that of the stick of the macroscopic stick-slip motion.The difference in the parameters, the driving method, and the local friction model may affect our results.However, FEM simulations employing similar parameters semi-quantitatively reproduce the external pressure dependence of the macroscopic static friction coefficient obtained in the experiment using PMMA 12 .The dependence of the macroscopic static friction coefficient on the aspect ratio for PMMA is also considered to be consistent with our present results.The dependence on material parameters, the driving methods, and the local friction model will be investigated in future work.
The dependence of the static and kinetic friction coefficients on the pressure or block shape has been studied in experiments using rubber blocks [30][31][32] .The results of these experiments are partially consistent with ours, but there is a difference in the dependence on the aspect ratio.In these experiments, the methods to change the aspect ratio and drive the block differ from those used in this study.For the rubber block, the local Amontons' law used in this study may not be applicable because the real contact area can become comparable to the apparent contact area, which contradicts the assumption of the Amontons' law.We need further investigations to determine the origin of the difference.
Recent numerical simulations of spring-block models have shown that the friction coefficient changes with the geometric pattern of the frictional interface [33][34][35][36] .However, our results indicate that an object shape can also control the macroscopic static friction coefficient.This might lead to new insights into methods for controlling friction in various objects, including shoe soles and tires.
Precursor slip has been investigated experimentally for the sliding motion of PMMA blocks based on fracture mechanics 26,[37][38][39][40][41] .Such a precursor slip is related to pre-earthquakes that occur a few days or months before a major earthquake [42][43][44] , which are studied using frictional spring-block models 45 .However, these studies have focused on 1D frictional interfaces or discrete models, which differ from 2D friction interfaces in more realistic systems.Our results for a 3D system with a 2D interface will provide new insights into the precursor slip observed in realistic situations.

Setting of system
The equation of motion for a viscoelastic body is given by with displacement u u u, stress σ σ σ , and second-order time derivative ü u u of displacement.The stress σ σ σ is given by the sum of the elastic stress σ σ σ (E) obeying Hooke's law and the viscous stress σ σ σ (V) , which is proportional to the strain rate.We assume that the viscoelastic body is isotropic.The elastic stress tensor σ (E) i j is given by with the Kronecker delta δ i j and the strain tensor ε i j .The viscous stress tensor σ (V) i j is given by with the strain rate tensor εi j 46 .The boundary conditions for the top surface at z = H are σ zz = −P ext and σ zx = σ zy = 0.At the free surface for x = 0,W or y = 0, L, we assume σ σ σ • n n n = 0 0 0 with the normal vector n n n of the surface.The boundary conditions at the contact surface with a rigid rod (y = 0) are given by σ yx = σ yz = 0 and uy = V rod , where uy is the velocity in the y direction and V rod is the velocity of rigid rod.At the bottom of the block (z = 0) in contact with a rigid substrate, the bottom pressure p = −σ zz is determined such that the displacement u z in the z direction is 0. However, the bottom pressure is limited to p ≥ 0. The region of the bottom surface with u z > 0 and p = 0 becomes a free surface with σ σ σ • n n n = 0 0 0. The boundary condition in the tangential direction at the bottom with p > 0 is given by with the tangential stress vector t t t(x, y) = (σ zx , σ zy ), local slip velocity vector v v v(x, y) = ( ux , uy ), velocity ux in the x direction, and velocity uy in the y direction.The direction of the frictional stress is opposite to that of the local slip velocity.Frictional stress is defined as σ (fric) (x, y) = |t t t|.The slip velocity is defined as v(x, y) = |v v v(x, y)|.
The frictional stress σ (fric) is given by equation ( 1).In the case v(x, y) = 0, the frictional stress is balanced with the local shear stress, where the maximum magnitude of the former is given by µ S p(x, y).The local friction coefficient µ(v) linearly decreases from µ S to µ K for 0 < v ≤ v c and µ K for v > v c .Amontons' law is expected to hold locally if the local region considered in the frictional interface contains a sufficiently large number of real contact points and has negligibly small spatial variations in internal stress 8,9,47 .
To treat static friction in the numerical simulation, we introduce a small velocity scale v e .The local friction coefficient µ(v) is given by We consider the limit v e → +0.The region with 0 ≤ v ≤ v e corresponds to static friction.The slip area A is defined as the region with v > v e .

Details of 3D FEM simulation
The viscoelastic block is divided into cubes with length ∆x consisting of six tetrahedra.The displacements and velocities within each element are approximated using a linear interpolation.We choose the characteristic velocity v e /V rod = 2.5 × 10 −2 such that v e /V rod 1 is satisfied.In the FEM simulations, we select ∆x/H = 1/40, ∆t/(H ρ/E) ≈ 10 −6 , and V rod ρ/E = 2 × 10 −5 .
We have confirmed that the numerical results do not change, even if we use smaller values.First, we apply an external uniform pressure P ext to the top surface and relax the system to an equilibrium state.After relaxation, the center of the side surface (x, y, z) = (W /2, 0, H/2) is pushed along the y direction by a rigid rod from time t = 0 with a sufficiently slow speed V rod .The displacement of the rigid rod is denoted by U(t) = V rod t.The length of one side of a rigid square rod is 0.1H, and the height of its center from the bottom is 0.5H.

Details of analysis based on simplified models
Model for large L/H: The second term on the right-hand side of equation ( 2) represents local friction.Here, we assume a constant bottom pressure given by P ext , which is verified in the FEM simulations for L/H 1 as shown in the Supplementary Note and Supplementary Video S1.The local friction coefficient µ is expressed as a function of v = | uy |.Note that 0 ≤ µ ≤ µ S when v = 0.The boundary conditions are ∂ u y (y = L,t)/∂ y = 0 and u y (y = 0,t) = U(t).In our analysis, we set the origin of U immediately after the bulk sliding and assume that the ratio of the frictional stress to P ext is equal to µ K at U = 0. Model for large W /H: The second term on the right-hand side of equation ( 9) represents the friction.The bottom pressure is almost independent of x in the FEM simulations, as shown in the Supplementary Note and Supplementary Video S2.Therefore, we assume a constant bottom pressure given by P ext .The boundary conditions are ∂ u y (|x| = W /2,t)/∂ x = 0 and u y (x = 0,t) = U(t).

Parameters
The parameters for the viscoelastic object are chosen as ν = 0.34, η 1 /(H √ ρE) = 2, and η 2 /η 1 = 1, whereas we set the parameters for the friction as µ S = 0.38, µ K = 0.1, and v c ρ/E = 3.4 × 10 −4 , following previous FEM simulations 10 .These values are different from those adopted for the experiment using PMMA 10,12 .The parameters for the PMMA blocks 12  Therefore, in the simplified model with L/H 1, we neglect the x dependence.The spatial distribution of pressure p at the bottom is shown in Fig. S1 (c).The pressure p is approximately equal to the pressure P ext at the top surface of the entire bottom plane.Based on this numerical result, we assume that the bottom pressure is uniform in the simplified model for L/H 1.The spatial distributions in the frictional interface for W /H = 4 with L/H = 1 are shown in Fig. S2, corresponding to Supplementary Video S2.Figures S2 (a) and (b) demonstrate the spatial distributions of the slip region and the ratio σ (fric) /p, respectively.The regions for precursor slip and σ (fric) /p ≈ µ S propagate from the center at x/H = 2, which is different from those for L/H 1 shown in Fig. S1.The characteristic behavior is the extension of the slip region along the x direction, as shown in Supplementary Video S2.Therefore, in the model for W /H 1, we neglect the dependence on y.The spatial distribution of pressure p at the bottom is shown in Fig. S2 (c).The variation in p is not strong.Therefore, we assume a constant pressure in the simplified model for W /H 1.
0 1 x/H x/H x/H  x/H x/H x/H Using the same procedure as that in Sec. 2, the equation for the critical area A c is obtained as For A c /A 0 1, A c /A 0 is expressed as where the exponents for P ext and L/H differ from those in equation ( 5) in the main text.The difference in the exponents results from the different assumptions for the bottom pressure adopted in both studies, i.e., uniform pressure in the present study and non-uniform pressure in the previous study.Note that the macroscopic friction coefficient µ M in this model is also given by equation (6) in the main text.In Fig. S3, we show µ M against P ext for various lengths, L/H with W /H = 1, obtained from the FEM simulations in the present study.The analytical results in ref. 1 given by equations (S43) and (6), and in the present study given by equations (S17) and ( 6) are also shown.The fitting parameter α is determined such that the analytical result matches the FEM results for L/H = 2 and P ext /E > 0.01. Figure S3 shows that the analytical results in the previous study better reproduce the results of the FEM analysis for L/H ≤ 2. This is because the nonuniformity of the bottom pressure is significant for L/H ≤ 2. However, the analytical results deviate from the FEM analysis for L/H ≥ 4, where the bottom pressure is almost uniform, as shown in Fig. S1 (c) and Supplementary Video S1 (c).Instead, the analytical results of the present study agree with the FEM analysis for L/H ≥ 4, as shown in Fig. 3 and Fig. S3.

Figure 1 .
Figure 1.Schematic of a 3D viscoelastic block on a fixed rigid substrate.

Figure 2 .
Figure 2. (a) Ratio F T /F N against the displacement of the rigid rod U for L/H = 1, W /H = 2, and P ext /E = 0.006.The red horizontal line represents the macroscopic static friction coefficient µ M .(b) Macroscopic static friction coefficient µ M against pressure P ext for various L/H values with W /H = 1.The thin solid lines represent the analytical results with α A = 0.2 given by equations (4) and (6).(c) Macroscopic static friction coefficient µ M against P ext for various W /H values with L/H = 1.The thin solid lines represent the analytical results with α B = 0.2 given by equations (6) and(11).The dotted and dashed lines represent µ S and µ K , respectively.

Figure 3 .
Figure 3. (a) Spatial distribution of the slip region in the frictional interface at U = U 1 ,U 2 ,U 3 , and U 4 for L/H = 1, W /H = 2, and P ext /E = 0.006.The yellow area represents the slip region.The rigid rod is pushing the block at (x/H, y/H) = (1, 0).(b) Normalized precursor slip area A/A 0 against displacement U.The dotted line represents the normalized critical area A c /A 0 .(c) Normalized critical area A c /A 0 against pressure P ext for various L/H values with W /H = 1.The thin solid lines represent the analytical results with α A = 0.2 given by equation (4).(d) Normalized critical area A c /A 0 against P ext for various W /H values with L/H = 1.The thin solid lines represent the analytical results with α B = 0.2 given by equation (11).

Figure 4 .Figure 5 .
Figure 4. Macroscopic static friction coefficient µ M against the normalized critical area A c /A 0 for various L/H and W /H values.The solid line represents the analytical result given by equation (6).The dotted and dashed lines represent µ S and µ K , respectively.

8 / 11 1
are estimated as L/H = 5, W /H = 0.25, P ext /E ≈ 3 × 10 −4 , ν = 0.4, µ S = 1.2, and µ K = 0.2, and much smaller v c ρ/E and η 1 /(H √ ρE) are used in the previous study 10 .Spatial distributions of quantities in the frictional interface This section presents the spatial distributions of the quantities in the frictional interface for the 3D FEM simulations, which supports the assumption used to derive the simplified models in the main text.The spatial distributions for L/H = 4 and W /H = 1 are shown in Fig. S1, corresponding to Supplementary Video S1.Figures S1 (a) and (b) demonstrate the spatial distributions of the slip region and the ratio σ (fric) /p, respectively.Both distributions exhibit a slight dependency on x.

Figure S1 .
Figure S1.Spatial distributions of different quantities in the frictional interface for P ext /E = 0.012, L/H = 4, W /H = 1, and U/L = 2.6 × 10 −2 .The rigid rod pushes the block at (x/H, y/H) = (0.5, 0).(a) Spatial distribution of the slip region, which is represented by the yellow area.(b) Spatial distribution of the ratio σ (fric) /p.(c) Spatial distribution of p.