Effects of small and large shear histories on multiple liquefaction properties of sand with initial static shear

It has been reported that soils belonging to slope grounds show different types of liquefaction behavior than those belonging to horizontal grounds. Some research has also revealed that liquefaction histories can significantly affect the shear behavior of sandy soils. However, the combined effects of the slope angle and the magnitude of past shear histories on the liquefaction properties of soils have not been studied comprehensively. Based on this background, several multiple liquefaction tests with initial static shear were conducted on Toyoura sand. In each of these tests, a single specimen was sheared several times up to small or large double amplitude shear strain under a constant volume condition using a specially designed stacked-ring shear apparatus. The behavior of the Toyoura sand observed in these tests was discussed considering various perspectives, such as the increase in relative density, the induced anisotropy, the change in liquefaction resistance, and the shear strain accumulation. The findings of this study established that shear histories of smaller magnitude had relatively less influence on densification and induced anisotropy than those of larger magnitude. Moreover, shear histories of smaller magnitude also resulted in the relatively higher liquefaction resistance of sand specimens against the next cyclic shear, while the opposite trend was observed in the case of specimens subjected to shear histories of larger magnitude. Finally, shear strain accumulated less easily in tests with small shear histories than in those with large shear histories.     

Effects of specimen density and initial cyclic loading history on correlation between shear wave velocity and liquefaction resistance of Toyoura sand

Some previous studies have shown a good correlation between the shear wave velocity, V_s, and the cyclic resistance ratio, CRR. Recently, however, a V_s-based liquefaction assessment method has become an alternative and supplementary method to the conventional N_SPT-based method. It is known that the CRR is influenced not only by the specimen density, but also by the soil fabric. Unfortunately, there are concerns that different combinations of the effects of the specimen density and the soil fabric may generate different relations between V_s and the CRR even if the tested specimens are of the same soil material. In the current study, a series of V_s measurements and undrained cyclic triaxial tests is performed on Toyoura sand specimens with different soil fabrics for three different specimen densities. The fabric of the specimens is varied by applying initial cyclic loading. The results of the V_s measurements indicate that the V_s of the specimen is affected by the initial cyclic loading histories, and the results of the undrained cyclic triaxial tests show that there is a good correlation between V_s and the CRR. However, the correlation varies depending on the specimen density even when the tested material is Toyoura sand only. In other words, the soil-type specific correlation between V_s and the CRR depends on the specimen density. Therefore, the results indicate that both V_s and the specimen density are necessary parameters for an accurate assessment of the CRR.     

Microscopic investigation into liquefaction resistance of pre-sheared sand: Effects of particle shape and initial anisotropy

Liquefaction resistance of sand can be either increased or reduced due to an undrained cyclic pre-shearing depending on the degree of pre-shearing, which hinders a better prediction of liquefaction potential to be established. The mechanism of such changes in liquefaction resistance has been poorly understood. This contribution aims to gain micromechanical insights into pre-shearing effects on liquefaction resistance of sand using discrete element method (DEM) simulations. In particular, effects of particle shape and initial anisotropy on liquefaction resistance are investigated. The simulation results from samples consisting of non-spherical particles with an initial anisotropy can qualitatively capture the mechanical responses observed in equivalent laboratory experiments. The samples which yielded a qualitative agreement with the laboratory results are further analyzed micromechanically, and the relationships between liquefaction resistance and some microscopic parameters before cyclic loading are discussed. Microscopic analyses reveal that mean mechanical coordination number is well correlated with liquefaction resistance, whereas liquefaction resistance is less sensitive to anisotropy in particle orientation induced by pre-shearing.     

Effects of induced anisotropy on multiple liquefaction properties of sand with initial static shear

The multiple liquefaction phenomenon has been attracting the attention of more and more researchers and engineers since the 2010–2011 Christchurch Earthquakes and the 2011 Great East Japan Earthquake. However, little has been known about the multiple liquefaction properties of sloped grounds. In this study, therefore, multiple liquefaction tests that consider the initial static shear stress, which have never been conducted before, were carried out with a special designed apparatus, the stacked-ring shear apparatus. A series of multiple liquefaction tests revealed that induced anisotropy, which is weak against the loading opposite to the direction of the initial static shear stress, was produced by the liquefaction of a sloped ground. As a result, a significant decrease in the liquefaction resistance during the next cyclic of shearing occurred. The indicators which influenced the magnitude of anisotropy were also discussed from the perspective of the reconsolidation procedures, the magnitude of initial static shear stress, and the type of ending of the previous liquefaction stage. In addition, it turned out that in the multiple liquefaction test with a larger initial static shear stress, the re-liquefaction resistance was higher because the shear stress opposite to the direction of the initial static shear stress, causing large negative dilatancy due to anisotropy, was smaller in that test.     

Aging effects on liquefaction resistance of sand estimated from laboratory investigation

Liquefaction resistance is known to increase concomitantly with the increase in time after construction or sedimentation. Nevertheless, the mechanisms of its aging effect on liquefaction have not been completely elucidated. To clarify the mechanisms of aging in sandy soils, the liquefaction resistance (CRR), initial and secant shear moduli (G₀ and G_sec), and laboratory penetration resistance of long-term consolidated sand specimens were examined using cyclic undrained triaxial tests, local small strain (LSS) tests equipped with bender elements (BEs), and penetration index tests, respectively. Based on the existing reports, the CRR was inferred from G₀, G_sec, and the laboratory penetration resistance. In the case of Toyoura sand of D_r = 40%, the CRR increased by about 14% with a 360-day consolidated specimen in the cyclic undrained triaxial tests. However, increases in the CRR evaluated from G₀ and the laboratory penetration resistance were nothing and only 2%, respectively. G_sec started to degrade at greater shear strain in the long-term consolidated specimens. An increase in the CRR, evaluated from the G_sec of 0.01% shear strain, had a much better agreement with that obtained from the cyclic undrained triaxial tests.     

Insight into excess pore pressure generation leading to liquefaction of sand with stress history under saturated and unsaturated conditions

Assessments of the liquefaction resistance of clean sand still involve considerable uncertainties, which are a current research topic in the field of soil liquefaction. The factors considered and discussed in this study include the loading history, degree of saturation, and partial drainage. The effects of each of these factors on pore pressure generation and liquefaction resistance have been studied for decades in the laboratory, and empirical relationships have been derived. In this paper, an attempt is made to explain these effects using the unique index of volumetric strain. A pore pressure generation model is developed which is similar to that of Martin et al. (1975), but based on stress-controlled triaxial tests. The model is verified through comparisons of its results with those of laboratory tests. It is confirmed that the plastic volumetric strain that has accumulated in sand, either by drained or undrained cyclic loading, dominates the increase in the liquefaction resistance of the sand. However, the plastic volumetric strain caused by overconsolidation is less effective in reducing the volumetric strain potential for subsequent cyclic shearing, thus enhancing its resistance to liquefaction. The model provides a better understanding of the physical processes leading to the liquefaction of saturated and unsaturated sand with and without stress history.     

A unified approach to link small-strain shear modulus and liquefaction resistance of pumiceous sand

Natural pumiceous (NP) sands containing pumice particles, a type of volcanic soil, are commonly found in the central part of the North Island in New Zealand. The pumice particles are highly crushable, compressible, lightweight and angular, making engineering assessment of their properties problematic. In this paper, several series of bender element and undrained cyclic triaxial tests were performed on reconstituted and undisturbed NP sands to determine their small-strain shear modulus (G_max) and cyclic resistance ratio (CRR). Furthermore, similar tests were also conducted on normal hard-grained sands (e.g., Toyoura sand) for the purpose of comparison. The results showed that the NP sands have considerably lower G_max compared to normal sands, resulting in their higher deformability during the initial stages of the cyclic loading test. The high angularity of NP sands play an important role toward the end of the cyclic loading and contributed to their higher CRR. Next, the ratio of CRR/G_max for each sample was correlated to a level of strain denoted as cyclic yield strain (ε_ay), which was found to be significantly dependent on the percentages of pumice particles present in the natural soils. On the other hand, the ε_ay was found to be less sensitive to the consolidation stress (σ′_c) and the relative density (Dr) of the materials. For example, over different values of σ′_c and Dr, NP sands have substantially higher values of cyclic yield strain due to their lower G_max and higher CRR when compared with those of ordinary sands.     

Evaluation of post-liquefaction volumetric strain of reconstituted samples based on soil compressibility

In this study, we performed undrained cyclic–torsional–shear tests on clean and silty sands to examine post-liquefaction volumetric strain. Accordingly, we proposed a novel density index—ultimate volumetric strain—to evaluate the density of clean sand and silty soils; the proposed density index is based on soil compressibility. We discussed the test results with respect to relative density, skeleton void ratio, and ultimate volumetric strain. Consequently, we confirmed that ultimate volumetric strain showed good correlation with post-liquefaction volumetric strain irrespective of soil type, whereas values for relative density and skeleton void ratio were discordant. Further, we observed that, compared to maximum shear strain, accumulated shear strain was a better index for representing the effect of loading history. From these results therefore, we proposed a chart to evaluate post-liquefaction volumetric strain based on ultimate volumetric strain and accumulated shear strain.     

Direct strain evaluation method for laboratory-based pillar performance

Pillar stability is one of important aspects for underground mines.Generally,the stability of the pillars is evaluated empirically based on case studies and site-specific rock mass conditions in mines.Nevertheless the empirical approach applicability can sometimes be constrained.The numerical-based approaches are potentially more useful as parametric studies can be undertaken and,if calibrated,can be more representative.Both empirical and numerical approaches are dependent on the strength evaluation of the pillars while the strain developing in the pillars is seldom taken into consideration.In this paper,gypsum and sandstone samples were tested in laboratory with different width-to-height ratios(W/H)to adapt the strain evaluation method to the laboratory-based pillars.A correlation was then developed between the strain and the width-to-height ratio for pillar monitoring purposes.Based on the results,a flowchart was created to conduct back analysis for the existing pillars to evaluate their stability and design new pillars,considering the strain analysis of the existing pillars with the W/H ratios modelled.     

Numerical study on liquefaction characteristics of granular materials under Rayleigh-wave strain conditions using 3D DEM

In liquefaction analyses, liquefaction is conventionally assumed to originate from the vertical propagation of shear waves. However, some field and theoretical evidence has demonstrated that the risk of liquefaction may be induced or increased by surface waves. In this study, the liquefaction characteristics of K₀-consolidated granular materials under Rayleigh-wave strain conditions, ideal deformation conditions under the assumption of constant volume, were investigated by the three-dimensional discrete element method (3D DEM). The results indicate that Rayleigh-wave strain conditions combine pure and simple shear modes. As the ratio of the shear strain amplitude to the normal strain amplitude (RSN) increases from 1 to +∞, granular materials tend to have a slower liquefaction rate and higher liquefaction resistance; however, the difference in the undrained responses is negligible when the RSN is less than 1. The distribution of the magnitude and orientation of the contact forces also varies with the RSN, while it is similar when the RSN is less than 1. The degradation of the skeleton structure and the evolution of the structural anisotropy accelerate the liquefaction of granular materials. Moreover, in situations with the same accumulated equivalent strain per cycle, the Rayleigh-wave strain condition with a low RSN value will make granular materials more vulnerable to liquefaction compared with Love- and SH-wave strain conditions.     

Deformation and cyclic resistance of sand in large-strain undrained torsional shear tests with initial static shear stress

This paper presents the findings from an experimental study focusing on the undrained cyclic behavior of sand in the presence of initial static shear stress. A series of undrained cyclic torsional shear tests was performed on saturated air-pluviated Toyoura sand specimens up to single amplitude shear strain ($\gamma$_SA) exceeding 50%. Two types of cyclic loading conditions, namely, stress reversal (SR) and stress non-reversal (SNR), were employed by changing the amplitude of the combined initial static shear and cyclic shear stresses. The tests covered a broad range of initial states in terms of relative density (Dr = 20–74%) and the initial static shear stress ratio (α = 0–0.30). The following five distinct modes of deformation were identified from the tests based on the density state, the transient undrained peak shear stress, and the combined cyclic and static shear stresses: 1) static liquefaction, 2) cyclic liquefaction, 3) cyclic mobility, 4) shear deformation failure, and 5) limited deformation. Of these, cyclic liquefaction and static liquefaction are the most critical. They occur in very loose sand (D_r ≤ 24%) under SR and SNR, respectively, and are characterized by abrupt flow-type shear deformation. Cyclic mobility occurs under SR in loose to dense sand with D_r ≥ 24%. Contrarily, shear deformation failure typically occurs under SNR in sand with 24 < D_r < 65%, and limited deformation may take place in dense sand with D_r ≥ 65%. In this paper, a stress-void ratio-based predictive method is proposed to identify the likely mode of deformation/failure in sand under undrained shear loading with static shear. Furthermore, the cyclic resistance is evaluated at three different levels of $\gamma$_SA (i.e., small, $\gamma$_SA = 3%; moderate, $\gamma$_SA = 7.5%; and large, $\gamma$_SA = 20%). The results show that, independent of the density state, the cyclic resistance continuously decreases with an increase in α at the small $\gamma$_SA level, while it first decreases and then increases for both loose and dense sand at the moderate and large $\gamma$_SA levels.     

Influence of inherent anisotropy on the seismic behavior of liquefiable sandy level ground

Inherent anisotropy is a crucial aspect to consider for an improved understanding of the strength and deformation characteristics of granular materials. It has been the focus of intense investigation since the mid-1960s. However, inherent anisotropy’s influence on ground seismic responses, such as liquefaction, has not been extensively studied. In this paper, inherent anisotropy’s influence on ground seismic responses is examined through a series of dynamic centrifuge model tests on liquefiable level sand deposits. During the model setup, five different deposition angles (0, 30, 45, 60, and 90 degrees) were achieved using a specially designed rigid container. The models were exposed to tapered sinusoidal input accelerations and the recorded results were fully investigated. It was found that deposition angle-caused inherent anisotropy significantly influenced the excess pore pressure responses during the shaking and dissipation phases. The amount of excess pore pressure build-up and the high excess pore pressure duration increased with the deposition angle, while the dissipation rate decreased as the deposition angle increased. The inherent anisotropy also influenced liquefaction-induced ground settlement, with volumetric strain increasing along with the deposition angle. With respect to response acceleration, inherent anisotropy’s effects depended on the amount of excess pore pressure build-up (i.e., degree of liquefaction). In view of these results, it was concluded that a sandy ground, deposited at a higher angle (i.e., closer to 90 degrees), is more susceptible to liquefaction and that inherent anisotropy’s influence should be considered when evaluating the liquefaction potential and performing effective stress analyses.     

Characterizing the effects of fines on the liquefaction resistance of silty sands

Recent earthquakes in New Zealand and Japan indicate that evaluating the liquefaction potential of silty sands remains an area of difficulty and uncertainty in geotechnical engineering. This paper presents a comprehensive experimental study along with analysis and interpretation in the framework of critical state soil mechanics, with the aim to address the complicated effects of fines. Two series of sand-silt mixtures, formed by mixing two different base sands with the same type of non-plastic silt, were tested under a range of packing density, confining pressure and silt content, and a unified correlation was established between the cyclic resistance and the state parameter that collectively accounts for the effect of packing density and confining pressure. The proposed correlation is independent of packing density, confining pressure, fines content and base sand, and allows prediction of the cyclic resistance of silty sands under different states. Furthermore, the mechanism of the fines-content induced reduction of cyclic resistance and the mechanism of the base-sand effect observed from the tests are elaborated in the sound theoretical context. The present study suggests that the critical state soil mechanics is a rational and appealing framework for liquefaction analysis of both clean and silty sands.     

Extended equivalent linear model (X-ELM) to assess liquefaction triggering: Application to the City of Urayasu during the 2011 Tohoku earthquake

This paper proposes a hybrid iterative approach to evaluating the liquefaction potential and settlements when liquefaction is not triggered for a saturated sand deposit subjected to seismic motion. The proposed method is a combination of an equivalent linear dynamic analysis and an empirical pore pressure build-up model. The two concepts are combined to account for the increase of pore-water pressure in the dynamic response. This extended equivalent linear model, or X-ELM, introduces a parameter $\chi$ identified by means of an extensive set of computations and observations. The approach has been validated in the case of Urayasu City during the 2011 Great East Japan Earthquake. The X-ELM enables the assessment of the triggering of the soil liquefaction of 12 representative soil profiles, namely, four profiles representing non-liquefied grounds in the old town of Urayasu and eight profiles representing the liquefaction-prone grounds in the reclaimed area of the city. The X-ELM computed areas of triggered liquefaction or non-liquefaction are found to be in good agreement with qualitative observations and in-situ measured responses.     

Using a state-dependent constitutive model in strain wedge method for laterally loaded piles in sand

Strain Wedge Method is an effective technique for the analysis of laterally loaded flexible piles in sand. In this method, the three-dimensional soil-pile interaction problem is reduced to the problem of one-dimensional beam on elastic foundation. Various modifications have been applied to the conventional Strain Wedge Method to promote its predictive capacity for different problems. Nevertheless, in the existing Strain Wedge Methods, because the impact of soil state on the stress-strain-strength behavior of granular soils is ignored, multiple calibrations for a wide domain of soil densities are required. To overcome this drawback, a simple but versatile state-dependent constitutive model is implemented in the Strain Wedge Method. It is shown that the modified Strain Wedge Method can effectively simulate the response of laterally loaded model piles embedded in sands of loose, medium, and dense states using a single set of parameters.     

“振冲碎石桩＋强夯”彻底消除砂层液化的施工技术

本文结合福建泰山石化仓储项目石油化工品罐区工程,主要介绍采用“振冲碎石桩＋强夯”的技术措施对地基进行处理,在提高地基承载力的同时,重点解决砂层在地震作用下的液化问题。     

Geo-mechanical behavior of clay soils stabilized at ambient temperature with fly-ash geopolymer-incorporated granulated slag

Geopolymer is a cementitious material that can replace ordinary Portland cement in several geotechnical engineering applications, such as soil stabilization, with the advantages of much lower harmful emissions and energy consumption. This paper presents a rigorous evaluation of the geo-mechanical behavior of different types of clay soils treated with geopolymer, including the influence of soil characteristics and mineralogy. Two natural clay soils in addition to a commercially available kaolin clay were used for this investigation. Laboratory experiments were performed including unconfined compressive strength (UCS) and consolidated undrained (CU) triaxial compression tests under different confining pressures. The UCS and triaxial tests indicated that the addition of geopolymer considerably increased the yield stress and initial stiffness of all examined clays. With the increase of geopolymer content, the stress–strain behavior of treated clays was found to develop progressively from ductile response into a post-peak brittle fashion. The CU tests also demonstrated that the addition of geopolymer changed the initial characteristics of remolded clays from quasi-over-consolidated to heavily over-consolidated, rendering high yield surface and more effective shear strength parameters (i.e., cohesion and friction angle). Moreover, although the overall qualitative stress–strain and stress path responses of the clays were similar, significant quantitative differences were observed, particularly in terms of the attainable yield strength, stiffness, and shear strength. These differences can be attributed mainly to the heterogeneity associated with the soil mineralogy and the corresponding differences in the interaction between the clay/non-clay minerals and geopolymer.     

A novel 2D-3D conversion method for calculating maximum strain of geosynthetic reinforcement in pile-supported embankments

For design of a geosynthetic-reinforced pile-supported (GRPS) embankment over soft soil, the methods used to calculate strains in geosynthetic reinforcement at a vertical stress were mostly developed based on a plane-strain or two-dimensional (2-D) condition or a strip between two pile caps. These 2-D-based methods cannot accurately predict the strain of geosynthetic reinforcement under a three-dimensional (3-D) condition. In this paper, a series of numerical models were established to compare the maximum strains and vertical deflections (also called sags) of geosynthetic reinforcement under the 2-D and 3-D conditions, considering the following influence factors: soil support, cap shape and pattern, and a cushion layer between cap and reinforcement. The numerical results show that the maximum strain in the geosynthetic reinforcement decreased with an increase of the modulus of subgrade reaction. The 2-D model underestimated the maximum strain and sag in the geosynthetic reinforcement as compared with the 3-D model. The cap shape and pattern had significant influences on the maximum strains in the geosynthetic reinforcements. An empirical method involving the geometric factors of cap shape and pattern, and the soil support was developed to convert the calculated strains of geosynthetic reinforcement in piled embankments under the 2-D condition to those under the 3-D condition and verified through a comparison with the results in the literature.     

Modelling the degradation of penetration resistance during cyclic T-bar tests in a Gulf of Mexico clay

Installation and repetitive movement during the operation of a pipeline causes remolding and softening of the surrounding soil. Similar effects can occur around foundations during their operating life. The degradation of the undrained shear strength of soft clays near the seabed is a critical component of the design of subsea facilities, including pipelines and shallow foundations. This paper presents a two-stage strength degradation model based on the framework developed by Hodder et al. (2010) for repeated vertical movement of a cylindrical object (T-bar penetrometer) embedded in a soft Gulf of Mexico (GOM) clay. The model is compared with the behavior observed in four box cores sampled in the Gulf of Mexico, in which 8 cyclic T-bar penetrometer tests were performed. By varying the parameters in the strength degradation model, its applicability for GOM clay is examined, and the effects of the location of cycling within the box core and parameter variability between box cores relative to the average parameters are explored.     

Shaking table tests on wall-type gravel and rubber drains as a liquefaction countermeasure in silty sand

In this paper, the seismic performance of wall-type gravel and rubber drains as a liquefaction countermeasure in silty sand has been addressed using a series of 1 g shaking table tests. In these tests, the liquefaction resistance of silty sand was studied by changing soil relative density, silt content, number of walls, and materials of drainage walls. In order to evaluate shear wave propagation, the generation and dissipation of pore water pressure (PWP), and the ground surface settlement, various accelerometers, displacement and PWP transducers were placed. The obtained results indicated that an increase in the relative density and the number of wall-type gravel and rubber drains in liquefiable silty sand reduces the settlement, liquefaction-induced deformations, as well as the excess PWP and eventually improves the liquefaction resistance. Totally, it can be noted that the reinforced silty sand with wall-type gravel drains revealed less excess PWP and settlement than that with wall-type rubber drains.