Pressure shock fronts formed by ultra-fast shear cracks in viscoelastic materials

Tensile cracks can shatter classical speed limits

Brittle materials fail by means of rapid cracks. Classical fracture mechanics describes the motion of tensile cracks that dissipate released elastic energy within a point-like zone at their tips. Within this framework, a "classical" tensile crack cannot exceed the Rayleigh wave speed, [Formula: see text]. Using brittle neo-hookean materials, we experimentally demonstrate the existence of "supershear" tensile cracks that exceed shear wave speeds, [Formula: see text]. Supershear cracks smoothly accelerate beyond [Formula: see text], to speeds that could approach dilatation wave speeds. Supershear dynamics are governed by different principles than those guiding "classical" cracks; this fracture mode is excited at critical (material dependent) applied strains. This nonclassical mode of tensile fracture represents a fundamental shift in our understanding of the fracture process.

The role of pore fluids in supershear earthquake ruptures

The intensity and damage potential of earthquakes are linked to the speed at which rupture propagates along sliding crustal faults. Most earthquakes are sub-Rayleigh, with ruptures that are slower than the surface Rayleigh waves. In supershear earthquakes, ruptures are faster than the shear waves, leading to sharp pressure concentrations and larger intensities compared with the more common sub-Rayleigh ones. Despite significant theoretical and experimental advances over the past two decades, the geological and geomechanical controls on rupture speed transitions remain poorly understood. Here we propose that pore fluids play an important role in explaining earthquake rupture speed: the pore pressure may increase sharply at the compressional front during rupture propagation, promoting shear failure ahead of the rupture front and accelerating its propagation into the supershear range. We characterize the transition from sub-Rayleigh to supershear rupture in fluid-saturated rock, and show that the proposed poroelastic weakening mechanism may be a controlling factor for intersonic earthquake ruptures.

Theoretical estimates of the parameters of longitudinal undular bores in polymethylmethacrylate bars based on their measured initial speeds

We study the evolution of the longitudinal release wave that is generated by induced tensile fracture as it propagates through solid rectangular polymethylmethacrylate (PMMA) bars of different constant cross-section. High-speed multi-point photoelasticity is used to register the strain wave at three distances from the fracture site in each experiment. In all cases, oscillations develop at the bottom of the release wave that exhibit the qualitative features of an undular bore. The pre-strain, post-strain, strain rate of the release wave and the cross-section dimensions determine the evolution of the oscillations. From the wave speed and strain rate close to the fracture site, we estimate the strain rate of the release wave as well as the growth of the amplitude and duration of the leading oscillation away from the fracture site by using formulae derived from the simple analytical solution of the linearized Gardner equation (linearized near the pre-strain level at fracture). Our estimates are then compared to experimental data, where good agreements of these three parameters are found between the predictions of the model and the experimental observations. Thus, we developed an approach to estimating the key characteristics of the developing unsteady undular bore based on the measured initial speeds of the longitudinal and shear waves. This does not require a prior knowledge of the elastic moduli for the conditions of the experiments, which in PMMA are known to be strain rate dependent.

Structure Characterization and Impact Effect of Al-Cu Graded Materials Prepared by Tape Casting

With the need of developing new materials, exploring new phenomenon, and discovering new mechanisms under extreme conditions, the response of materials to high-pressure compression attract more attention. However, the high-pressure state deviating from the Hugoniot line is difficult to realize by conventional experiments. Gas gun launching graded materials could reach the state. In our work, the corresponding Al-Cu composites and graded materials are prepared by tape casting and hot-pressing sintering. The microstructure and the acoustic impedance of the corresponding Al-Cu composites are analyzed to explain the impact behavior of Al-Cu graded materials. Computed tomographic testing and three-dimension surface profilometry machine results demonstrated well-graded structure and parallelism of the graded material. Al-Cu GMs with good parallelism are used to impact the Al-LiF target at 2.3 km/s using a two-stage light-gas gun, with an initial shock impact of 20.6 GPa and ramping until 27.2 GPa, deviating from the Hugoniot line.

Dynamic rupture initiation and propagation in a fluid-injection laboratory setup with diagnostics across multiple temporal scales

Fluids are known to trigger a broad range of slip events, from slow, creeping transients to dynamic earthquake ruptures. Yet, the detailed mechanics underlying these processes and the conditions leading to different rupture behaviors are not well understood. Here, we use a laboratory earthquake setup, capable of injecting pressurized fluids, to compare the rupture behavior for different rates of fluid injection, slow (megapascals per hour) versus fast (megapascals per second). We find that for the fast injection rates, dynamic ruptures are triggered at lower pressure levels and over spatial scales much smaller than the quasistatic theoretical estimates of nucleation sizes, suggesting that such fast injection rates constitute dynamic loading. In contrast, the relatively slow injection rates result in gradual nucleation processes, with the fluid spreading along the interface and causing stress changes consistent with gradually accelerating slow slip. The resulting dynamic ruptures propagating over wetted interfaces exhibit dynamic stress drops almost twice as large as those over the dry interfaces. These results suggest the need to take into account the rate of the pore-pressure increase when considering nucleation processes and motivate further investigation on how friction properties depend on the presence of fluids.

Evolution of dynamic shear strength of frictional interfaces during rapid normal stress variations

Pressure shear plate impact tests have revealed that when normal stress changes rapidly enough, the frictional shear resistance is no longer proportional to the normal stress but rather evolves with slip gradually. Motivated by these findings, we focus on characterizing the dynamic shear strength of frictional interfaces subject to rapid variations in normal stress. To study this problem, we use laboratory experiments featuring dynamic shear cracks interacting with a free surface and resulting in pronounced and rapid normal stress variations. As dynamic cracks tend to propagate close to the wave speeds of the material, capturing their behavior poses the metrological challenge of resolving displacements on the order of microns over timescales microseconds. Here we present our novel approach to quantify the full-field behavior of dynamic shear ruptures and the evolution of friction during sudden variations in normal stress, based on ultrahighspeed photography (at 1-2 million frames/sec) combined with digital image correlation. Our measurements allow us to capture the evolution of dynamic shear cracks during these short transients and enable us to decode the nature of dynamic friction.

Enhanced Digital Image Correlation Analysis of Ruptures with Enforced Traction Continuity Conditions Across Interfaces

Accurate measurements of displacements around opening or interfacial shear cracks (shear ruptures) are challenging when digital image correlation (DIC) is used to quantify strain and stress fields around such cracks. This study presents an algorithm to locally adjust the displacements computed by DIC near frictional interfaces of shear ruptures, in order for the local stress fields to satisfy the continuity of tractions across the interface. In the algorithm, the stresses near the interface are extrapolated by local polynomials that are constructed using a constrained inversion. This inversion is such that the traction continuity (TC) conditions are satisfied at the interface while simultaneously matching the displacements produced by the DIC solution at the pixels closest to the center of the subset, where the DIC fields are more accurate. We apply the algorithm to displacement fields of experimental shear ruptures obtained using a local DIC approach and show that the algorithm produces the desired continuous traction field across the interface. The experimental data are also used to examine the sensitivity of the algorithm against different geometrical parameters related to construction of the polynomials in order to avoid artifacts in the stress field.