No Arabic abstract
This article focuses on liquefaction of saturated granular soils, triggered by earthquakes. Liquefaction is definedhere as the transition from a rigid state, in which the granular soil layer supports structures placed on its surface, toa fluidlike state, in which structures placed initially on the surface sink to their isostatic depth within the granularlayer.We suggest a simple theoretical model for soil liquefaction and show that buoyancy caused by the presence ofwater inside a granular medium has a dramatic influence on the stability of an intruder resting at the surface of themedium.We confirm this hypothesis by comparison with laboratory experiments and discrete-element numericalsimulations. The external excitation representing ground motion during earthquakes is simulated via horizontalsinusoidal oscillations of controlled frequency and amplitude. In the experiments, we use particles only slightlydenser than water, which as predicted theoretically increases the effect of liquefaction and allows clear depth-of-sinkingmeasurements. In the simulations, a micromechanical model simulates grains using molecular dynamicswith friction between neighbors. The effect of the fluid is captured by taking into account buoyancy effects onthe grains when they are immersed. We show that the motion of an intruder inside a granular medium is mainlydependent on the peak acceleration of the ground motion and establish a phase diagram for the conditions underwhich liquefaction happens, depending on the soil bulk density, friction properties, presence of water, and peak acceleration of the imposed large-scale soil vibrations.We establish that in liquefaction conditions, most cases relaxtoward an equilibrium position following an exponential in time.We also show that the equilibrium position itself,for most liquefaction regimes, corresponds to the isostatic equilibrium of the intruder inside a medium of effectivedensity. The characteristic time to relaxation is shown to be essentially a function of the peak ground velocity.
Soil liquefaction is a significant natural hazard associated with earthquakes. Some of its devastating effects include tilting and sinking of buildings and bridges, and destruction of pipelines. Conventional geotechnical engineering practice assumes liquefaction occurs via shear-driven compaction and consequent elevation of pore pressure. This assumption guides construction for seismically hazardous locations, yet evidence suggests that liquefaction strikes also under currently unpredicted conditions. Here we show, using theory, simulations and experiments, another mechanism for liquefaction in saturated soils, without necessitating high pore fluid pressure or special soils, whereby seismically triggered liquefaction is controlled by buoyancy forces. This new mechanism supplements the conventional pore pressure mechanism, enlarges the window of conditions under which liquefaction is predicted to occur, and may explain previously not understood cases such as liquefaction in well-compacted soils, under drained conditions, repeated liquefaction cases, and the basics of sinking in quicksand. These results may greatly impact hazard assessment and mitigation in seismically active areas.
During an earthquake, part of the released elastic strain energy is dissipated within the slip zone by frictional and fracturing processes, the rest being radiated away via elastic waves. Frictional heating thus plays a crucial role in the energy budget of earthquakes, but, to date, it cannot be resolved by seismological data. Here we investigate the dynamics of laboratory earthquakes by measuring frictional heat dissipated during the propagation of shear instabilities at typical seismogenic depth stress conditions. We perform, for the first time, the full energy budget of earthquake rupture and demonstrate that increasing the radiation efficiency, i.e. the ratio of energy radiated away via elastic waves compared to that dissipated locally, increases with increasing thermal - frictional - weakening. Using an in-situ carbon thermometer, we map frictional heating temperature heterogeneities - heat asperities - on the fault surface. Combining our microstructural, temperature and mechanical observations, we show that an increase in fault strength corresponds to a transition from a weak fault with multiple strong asperities, but little overall radiation, to a highly radiative fault, which behaves as a single strong asperity.
Soil has been recognized as an indirect driver of global warming by regulating atmospheric greenhouse gases. However, in view of the higher heat capacity and CO2 concentration in soil than those in atmosphere, the direct contributions of soil to greenhouse effect may be non-ignorable. Through field manipulation of CO2 concentration both in soil and atmosphere, we demonstrated that the soil-retained heat and its slow transmission process within soil may cause slower heat leaking from the earth. Furthermore, soil air temperature was non-linearly affected by soil CO2 concentration with the highest value under 7500 ppm CO2. This study indicates that the soil and soil CO2, together with atmospheric CO2, play indispensable roles in fueling the greenhouse effect. We proposed that anthropogenic changes in soils should be focused in understanding drivers of the globe warming.
The analysis of strong motion recordings in structures is crucial to understand the damaging process during earthquakes. A very precise time-frequency representation, the reassigned smoothed pseudo-Wigner-Ville method, allowed us to follow the variation of the Millikan Library (California) and the Grenoble City Hall building (France) resonance frequencies during earthquakes. Under strong motions, a quick frequency drop, attributed to damage of the soil-structure system, followed by a slower increase is found. However, in the case of weak earthquakes, we show that frequency variations come from the ground motion spectrum and cannot be interpreted in terms of change of the soil-structure system.
Desiccation cracking in clayey soils occurs when they lose moisture, leading to an increase in their compressibility and hydraulic conductivity and hence significant reduction of soil strength. The prediction of desiccation cracking in soils is challenging due to the lack of insights into the complex coupled hydro-mechanical process at the grain scale. In this paper, a new hybrid discrete-continuum numerical framework, capable of capturing hydro-mechanical behaviour of soil at both grain and macro scales, is proposed for predicting desiccation cracking in clayey soil. In this framework, a soil layer is represented by an assembly of DEM particles, each occupies an equivalent continuum space and carries physical properties governing unsaturated flow. These particles move freely in the computational space following the discrete element method (DEM), while their contact network and the continuum mixture theory are used to model the unsaturated flow. The dependence of particle-to-particle contact behaviour on water content is represented by a cohesive-frictional contact model, whose material properties are governed by the water content. In parallel with the theoretical development is a series of experiments on 3D soil desiccation cracking to determine essential properties and provide data for the validation of mechanical and physical behaviour. Very good agreement in both physical behaviour (e.g. evolution of water content) and mechanical behaviour (e.g. occurrence and development of cracks, and distribution of compressive and tensile strains) demonstrates that the proposed framework is capable of capturing the hydro-mechanical behaviour of soil during desiccation. The capability of the proposed framework facilitates numerical experiments for insights into the hydro-mechanical behaviour of unsaturated soils that have not been possible before.