Elastic Nonlinear Response in Granular Media under Resonance ConditionsX. Jia* and P. A. Johnson**
*Laboratoire de Physique des Matériaux Divisés et des Interfaces, Université de MarnelaVallée, France **Geophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545 AbstractWe are studying the elastic linear and nonlinear behavior of granular media using dynamic wave methods. In the work presented here, our goal is to quantify the elastic nonlinear response by applying both wave resonance and pulsemode measurements. Resonance studies are desirable because they provide the means to easily study amplitude dependencies of elastic nonlinear behavior and thus to characterize the physical nature of the elastic nonlinearity. Pulsemode studies provide information about the time scale of the onset of nonlinear response. This work has implications for a variety of topics, in particular, the in situ nonlinear response of surface sediments and dynamic earthquake triggering. For this work we constructed an experimental cell in which high sensitivity dynamic resonance studies were conducted using granular media under controlled effective pressure. We limit our studies here to bulk modes but have the capability to employ shear waves as well. The granular media are composed of glass beads held under pressure by a piston, while applying resonance waves from transducers as both the excitation and the material probe. The same apparatus is used in pulsemode experiments. We find significant elastic nonlinearity at all effective pressures studied, manifest by the fundamentalmode resonance curves decreasing progressively, at progressively increasing drive level. This is equivalent to progressive material softening with wave amplitude, meaning the wave speed and modulus diminish. The wave dissipation simultaneously increases (Johnson and Sutin, 2004). For example, at 0.11 Mpa effective pressure the observed change in resonance frequency of about 2.6% corresponds to a material bulk modulus decrease of about 5.2%. Strain amplitudes are 107106. Thus, we would predict that surface sediments should have significant elastic nonlinear response beginning at about 106 strain amplitude. Experimental ConfigurationHigh sensitivity dynamic resonance and pulsemode studies were conducted using granular media, glass beads, under controlled effective pressure (see figure below). The glass beads are held under pressure by a piston, while applying resonance waves from transducers as both the excitation and the material probe. The container is closed with two fitted pistons and a normal load is applied to the granular sample across the top piston. Force and displacement are measured directly. Resonance studies are desirable because they provide the means to easily study amplitude dependencies of elastic nonlinear behavior and thus to characterize the physical nature of the elastic nonlinearity. Pulsemode experiments are then conducted to study the time scale of the material elastic nonlinear response. Resonant frequency sweeps with frequencies corresponding to the fundamental bulk mode are conducted with the apparatus. The glass beads used in our experiments are of diameter 0.5 mm, deposited in a duralumin cylinder of diameter 30 mm and height of 15 mm. This corresponds to a granular skeleton acoustic wave velocity of 750 m/s under 50 N of force [0.07 Mpa] (Jia, 2004). The loaded system gives fundamental mode resonances in the audio frequency band at approximately half a wavelength. Planewave generating and detecting transducers of diameter 30 mm are placed on axis at the top and bottom of the cylindrical container in direct contact with the glass beads. The wave signals are detected using a lockin amplifier. Drive frequency is swept from below to above the resonance mode. The resonance frequency at peak amplitude corresponds directly to modulus: where fR is resonance frequency, c is bulk wave speed, M is bulk modulus, m is shear modulus, L is sample length, and r is density. Drive frequency is swept from below to above the resonance mode. A typical frequency sweep is 3 kHz in width with a frequency sampling of 6 Hz. Frequency sweeps are applied at progressively increasing drive voltages to test for nonlineardynamical induced modulus softening. The resonance frequency at peak amplitude corresponds directly to modulus. Material Softening due to Nonlinear Response under Resonance ConditionsFundamentalmode resonancecurves taken under 0.11 MPa effective pressure (80 N force; the pore pressure is 1 atm): Resonance frequency of 17,700 Hz progressively decreases to 17,250 Hz with excitation amplitude. Wave dissipation correspondingly and dramatically increases, seen by the decrease in the normed peakamplitude and the increasing width of the curve. The change in modulus with detected strain amplitude at five effective pressures noted. The Yaxis shows the change in the bulk modulus (DM = M  M0) normalized to the equilibrium (elastically linear) value M0. The righthad Yaxis shows the percent change of the modulus. The Xaxis is the wave strain amplitude. Slow Dynamics in Granular MediaOne of the most intriguing aspects of the dynamical nonlinear behavior of granular media is Slow Dynamics, a phenomenon also observed in rock to a lesser degree (see, for example, Guyer and Johnson, 1999). After the wave disturbance, the material modulus and Q do not immediately return to their original, equilibrium values. Recovery takes hours to days depending on the effective pressure and appears to recover with the log of time. The figure below shows the recovery of the modulus under two different effective pressures. This behavior has never before been observed in granular media. Modulus Decrease with Amplitude in PulseModeThe source is a sinusoidal, onecycle pulse composed of 50 kHz, (duration 20 ms), driven at low (linear elastic) to progressively higher (nonlinear elastic) drive amplitudes. The effective pressure is 0.14 Mpa. Modulus is normalized to its (linearelastic) value. This experiment shows that the modulus decrease is immediate with a pulse, although less significant than under resonance conditions, because wave amplitudes are smaller. Scaling and Physical Nature of the NonlinearityThe scaling relation between the normalized change in frequency Df/f0 and the amplitude provides clues to the physical nature of the elastic nonlinearity, based on existing models. A classical model of solids in one dimension, based on anharmonicity, is described by, where s is stress, e is average strain, and b,d are higher order nonlinear contributions. This model, derived by Landau and Lifshitz (1986) for instance, tells us that the scaling between Df/f0 and the amplitude goes as strain squared. This is not what we observe here, and that is sensible because Landau theory is derived for atomiclevel elastic nonlinearity. Several phenomenological models containing hysteresis in the stressstrain relation (see next column) do describe the scalings we observe, and result in an equation such as, where a describes the hysteresis in the material (see, e.g., Guyer and Johnson (1999). The equation predicts a scaling of Df/f0 and amplitude that is linear, as we (approximately) observe in the granular media below. On the Physical Origin of the BehaviorAs noted, we have a phenomenological description of the elastic behavior. The actual physics has its origin in the graintograin interactions that can be described by HertzMindlin theory (see, e.g., Ostrovsky and Johnson, 2001) for instance. Hertz Mindlin theory describes the the forces between interacting spheres, where hysteretic effects arise from contained fluids, fluid films and shearing during pressure oscillation. Force Chains in Granular MediaThe two figures below show how, under modest pressures, granular media form "force chains" (e.g., Jaeger et al., 1996). Force chains are the network that literally carry the load, and the paths by which sound travels. In the blue figure, the media is placed in torsion, forced from the center. In the second figure, the material is under hydrostatic load; in both cases the bright regions show the chains. The nonlinear response of the material would arise primarily in the force chains in our view, rather than in the entire bulk. This must be rigorously demonstrated however. Summary and ConclusionsThe physics of sand piles is an area of intense research (e.g., Jaeger et al. 1996). We have described dynamical measurements in granular media under resonance and pulsemode conditions. The measurements indicate that granular media is highly nonlinear and hysteretic. Granualr media also exhibit the intriguing phenomenon called slow dynamics, a creeplike phenomenon where, after dynamic wave excitation, the material modulus takes hours to days to return back to its original state. We have also shown that, as effective pressure increases, the nonlinear response decreases. In addition to the intrinsic interest of basic research, regarding the physics of granular media, studies of it's nonlinear behavior have import on understanding soil behavior during strong ground motion and of the behavior of fault gouge, especially in regards to dynamic earthquake triggering. Presentations on these topics are offered during this meeting. Please see: Abstract Reference Number: 9919 Abstract Reference Number: 7055 ReferencesGuyer, R. A., and P. A. Johnson, The astonishing case of mesoscopic elastic nonlinearity, Physics Today, 52, 3035, (1999) . Jia, X., Codalike multiple scattering of elastic waves in dense granular media, Phys. Rev. Lett. 93, 154303154307 (2004). Jaeger, H.M., S. Nagel, and R. P. Behringer, The physics of granular media, Physics Today 49, 3238 (1996). Johnson, P. and A. Sutin, Slow dynamics in diverse solids, J. Acoust. Soc Am., in press (2004). Laudau, L. D. and E. M Lifshitz, Theory of Elasticity, Pergammon New York (1986). Ostrovsky, L. and P. Johnson, Dynamic nonlinear elasticity in geomaterials, Rivista del Nuovo Cimento, 24, 146 (2001). AcknowledgementsWork funded by the DOE Office of Basic Energy Science, Geoscience Research Program (USA) and by the Centre de la Recherche Scientifique (CNRS) [France}. HTML converted by Jeff Simpson.
