A cohesive zone model for dynamic shear faulting based on experimentally inferred constitutive relation and strong motion source parameters

Publisher Website
Google Scholar
Abstract
To understand constitutive behavior near the rupture front during an earthquake source shear failure along a preexisting fault in terms of physics, local breakdown processes near the propagating tip of the slipping zone under mode II crack growth condition have been investigated experimentally and theoretically. A physically reasonable constitutive relation between cohesive stress τ and slip displacementD, τ = (τi− τd)[1 + α log (1 + βD)] exp (−ηD) + τd, is put forward to describe dynamic breakdown processes during earthquake source failure in quantitative terms. In the above equation, τiis the initial shear stress on the verge of slip, τdis the dynamic friction stress, and α, β, and η are constants. This relation is based on the constitutive features during slip failure instabilities revealed in the careful laboratory experiments. These experiments show that the shear stress first increases with ongoing slip during the dynamic breakdown process, the peak stress is attained at a very (usually negligibly) small but nonzero value of the slip displacement, and then the slip‐weakening instability proceeds. The model leads to bounded slip acceleration and stresses at and near the dynamically propagating tip of the slipping zone along the fault in an elastic continuum. The dynamic behavior near the propagating tip of the slipping zone calculated from the theoretical model agrees with those observed during slip failure along the preexisting fault much larger than the cohesive zone. The model predicts that the maximum slip accelerationbe related to the maximum slip velocityand the critical displacementDcby, wherekis a numerical parameter, taking a value ranging from 4.9 to 7.2 according to a value of τi/τp(τpbeing the peak shear stress) in the present model. The model further predicts thatbe expressed in terms ofand the cut off frequencyfmaxsof the power spectral density of the slip acceleration on the fault plane asand thatin terms ofDcandfmaxsas. These theoretical relations agree well with the experimental observations and can explain interrelations between strong motion source parameters for earthquakes. The pulse width of slip acceleration on the fault plane is directly proportional to the timeTcrequired for the crack tip to break down, andfmaxsis inversely proportional toTc.

This publication has 31 references indexed in Scilit: