No Arabic abstract
The interest of the mining industry on the assessment of tailings static liquefaction has exacerbated after recent failures of upstream-raised tailings storage facilities (TSF). Standard practices to evaluate global stability of TSFs entail the use of limit equilibrium analyses considering peak and residual undrained shear strengths; thus, neglecting the work input required to drive the softening process that leads to progressive failure of susceptible tailings. This paper presents a simplified procedure to evaluate the static liquefaction triggering of upstream-raised TSFs by means of finite element models employing the well-known Hardening Soil model with small-strain stiffness (HSS). A calibration methodology is proposed to overcome the model limitation of not being implemented in a critical state framework, focusing on the stiffness parameters that control the rate of shear-induced plastic volumetric strains. A real TSF is modelled in Plaxis 2D to evaluate its vulnerability to liquefy due to an undrained lateral spreading at the foundation. Results show that minor movements near the toe induce the material into a strain-softening regime that leads to a progressive failure towards the structure crest.
Recent failures of upstream-raised tailings storage facilities (TSF) raised con-cerns on the future use of these dams. While being cost-effective, they entail higher risks than conventional dams, as stability largely relies on the strength of tailings, which are loose and normally-consolidated materials that may exhibit strain-softening during un-drained loading. Current design practice involves limit equilibrium analyses adopting a fully-softened shear strength; while being conservative, this practice neglects the work input required to start the softening process that leads to progressive failure. This paper describes the calibration and application of the NGI-ADPSoft constitutive model to evaluate the potential of static liquefaction of an upstream-raised TSF and provides an indirect measure of resilience. The constitutive model incorporates undrained shear strength anisotropy and a mesh-independent anisotropic post-peak strain softening. The calibration is performed using laboratory testing, including anisotropically-consolidated triaxial compression tests and direct simple shear tests. The peak and residual undrained shear strengths are validated by statistical interpretation of the available CPTu data. It is shown that this numerical exercise is useful to verify the robustness of the TSF design.
A likely source of earthquake clustering is static stress transfer between individual events. Previous attempts to quantify the role of static stress for earthquake triggering generally considered only the stress changes caused by large events, and often discarded data uncertainties. We conducted a robust two-fold empirical test of the static stress change hypothesis by accounting for all events of magnitude M>=2.5 and their location and focal mechanism uncertainties provided by catalogs for Southern California between 1981 and 2010, first after resolving the focal plane ambiguity and second after randomly choosing one of the two nodal planes. For both cases, we find compelling evidence supporting the static triggering with stronger evidence after resolving the focal plane ambiguity above significantly small (about 10 Pa) but consistently observed stress thresholds. The evidence for the static triggering hypothesis is robust with respect to the choice of the friction coefficient, Skemptons coefficient and magnitude threshold. Weak correlations between the Coulomb Index (fraction of earthquakes that received positive Coulomb stress change) and the coefficient of friction indicate that the role of normal stress in triggering is rather limited. Last but not the least, we determined that the characteristic time for the loss of the stress change memory of a single event is nearly independent of the amplitude of the Coulomb stress change and varies between ~95 and ~180 days implying that forecasts based on static stress changes will have poor predictive skills beyond times that are larger than a few hundred days on average.
Time-history deformation analyses of upstream-raised tailings dams use seismic records as input data. Such records must be representative of the in-situ seismicity in terms of a wide range of intensity measures (IMs) including peak ground acceleration (PGA), Arias intensity (AI), cumulative absolute velocity (CAV), source-to-site distance, duration, among others. No single IM is a sufficient descriptor of a given seismic demand (e.g. crest settlement) because different records, all of them compliant with any IM, can produce a very wide range of results from insignificant damage to global failure. The use of brute force, where hundreds of seismic records compliant with a set of IMs are employed, has proven to be a reasonable workaround of this limitation, at least able to produce a probabilistic density function of demand indicators. This procedure, however, requires a large number of runs, and is therefore expensive and time-consuming. Analyses can be optimized if an a priori simple tool is used to predict which seismic records would yield a given demand, thus obtaining estimations with much fewer runs. In order to perform a more precise selection, a semi-analytical screening procedure is presented in this paper. The procedure makes use of the spectral properties of the seismic record, considering only the intensity of the frequency content which is not filtered by the dam to obtain an a priori estimate of demand, expressed in this case in terms of displacements. The tool is validated using analytical and numerical models that prove insensitivity to the constitutive model used in the analysis, and is applied to a large tailings dam subjected to strong earthquakes.
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.
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.