Influence of strain histories on liquefaction resistance of sand

The liquefaction resistance of sand increases with cyclic pre-shearing and pre-shaking as a result of earthquakes if the strain level in the pre-shearing is small. When larger shear strains are imposed, liquefaction resistance decreases. These complicated effects of pre-shearing histories on the liquefaction resistance are investigated in this study through a series of cyclic triaxial tests. Various combinations of cyclic stress amplitude and number of cycles of pre-shearing are examined. The tested sand is Toyoura Sand at 45% relative density, under a confining pressure of 50 kPa. Test results indicate that for the range of shear strain amplitude in pre-shearing smaller than 0.35%, the liquefaction resistance increases with pre-shearing. The increase in the liquefaction resistance depends strongly on the volumetric strain in the pre-shearing, and several effects of the shear stress amplitude and number of cycles can be negligible. Small volumetric strain of the order of 1% doubled the liquefaction resistance. Meanwhile, in the range of shear strain amplitude larger than 0.6%, the liquefaction resistance decreases. The liquefaction resistance decreases as the shear strain amplitude increases. Shear strain amplitude is one of the factors dominating this degrading effect, and the volumetric strain exerts beneficial effects to a certain extent. In this study, another series of tests are conducted to investigate the combined effects of small and large strain amplitude pre-shearing. It is observed that small shear strain pre-shearing cycles subsequent to large shear strain cycles erased the degrading effect of the latter. However, a large shear strain pre-shearing after small strain cycles degrades the beneficial effect of the small shear strain pre-shearing cycles previously applied to the specimens; however, the effects of the former small strain pre-shearing remains.     

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.     

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.     

Linking inherent anisotropy with liquefaction phenomena of granular materials by means of DEM analysis

The liquefaction phenomena of sands have been studied by many researchers to date. Laboratory element tests have revealed key factors that govern liquefaction phenomena, such as relative density, particle size distribution, and grain shape. However, challenges remain in quantifying inherent anisotropy and in evaluating its impact on liquefaction phenomena. This contribution explores the effect of inherent anisotropy on the mechanical response of granular materials using the discrete element method. Samples composed of spherical particles are prepared which have approximately the same void ratio and mean coordination number (CN), but varying degrees of inherent anisotropy in terms of contact normals. Their mechanical responses are compared under drained and undrained triaxial monotonic loading as well as under undrained cyclic loading. The simulation results reveal that cyclic instability followed by liquefaction can be observed for loose samples having a large degree of inherent anisotropy. Since a sample having initial anisotropy tends to deform more in its weaker direction, leading to lower liquefaction resistance, a sample having an isotropic fabric potentially exhibits the greatest liquefaction resistance. Moreover, the effective stress path during undrained cyclic loading is found to follow the instability and failure lines observed for static liquefaction under undrained monotonic loading. From a micromechanical perspective, the recovery of effective stress during liquefaction can be observed when a threshold CN develops along with the evolving induced anisotropy. Realising that the conventional index of the anisotropic degree (a) is not effective when the CN drops to almost zero during cyclic liquefaction, this contribution proposes an alternative index, effective anisotropy (a×CN), with which the evolution of induced anisotropy can be tracked effectively, and common upper and lower bounds can be defined for both undrained monotonic and cyclic loading tests.     

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.     

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.     

Brief note on the influence of shape and percentage of gravel on the shear strength of sand and gravel mixtures

The paper records the influence of the shape and the percentage of gravel on the shear strength/frictional angle of sand and gravel mixtures using direct shear tests. The shear strength is mainly derived from the frictional forces developed due to sliding and interlock; they depend on the maximum particle size and shape, the uniformity coefficient, density and the effective normal stress. As the size of material in a mixture is variable, the shear strength also depends upon the ratio of the specimen diameter to the maximum particle size. In this study, two different shapes of limestone were used, angular and rounded, and the maximum gravel size was 6.3 mm in diameter. Air-dried samples were used in the tests. It is concluded that the shape and percentage of gravel have an important influence on the shear strength properties. Electronic Publication     

Morphological insights into the liquefaction and post-liquefaction response of sands with geotextile inclusions using drained constant volume simple shear tests

The present study aims to explore and bring out morphological insights into the prior-liquefaction, liquefaction, and post-liquefaction response of sands with geotextile inclusions. For this, a series of multi-stage drained constant volume simple shear tests with different cyclic stress ratios (CSR ranging from 0.1125 to 0.225) and different frequencies (f of 0.2 and 1.0 Hz) were carried out on completely dry specimens constituted with granular materials of three distinct grain morphologies (rounded, subrounded, and angular) reinforced with a nonwoven geotextile. The study also consists of morphological quantifications through image analysis algorithms and direct shear tests on sand-geotextile interfaces. Test results revealed that the inclusion of geotextile increased the liquefaction resistance and post-liquefaction shear strength of all the materials, irrespective of their particle morphology. However, the beneficial effects are more in the case of specimens constituted with angular particles. The effect of loading frequency on the response is also established. The interlocking and ploughing tendency of the angular particles leads to the mobilization of the maximum tensile strength of geotextile, which enhances the additional confinement and prevents the lateral movement of particles, thereby providing the maximum benefit.     

Investigation of the ultimate particle size distribution of a carbonate sand

A series of ring shear tests were conducted to investigate the ultimate particle size distribution of a carbonate sand. The tests were carried out under different stress levels, on three types of specimens: 1) uniformly graded specimens made of dry natural sand 2) remoulded specimens of the crushed sand after first shearing to large strains 3) specimens made of natural sand grains but with the same grading as in (2). The first series of tests on type (1), carried out to very large strains, led to apparently stable gradings, distinct for each stress level. Only limited additional particle breakage could be induced by remoulding the specimens after shearing (type (2)) and subjecting them to more shearing. Tests on specimens created at the apparently stable gradings (type (3)) but from the intact sand particles however led to significantly greater breakage. For the three types a stable, fractal grading was achieved. Analyses of the soil particles’ shape showed that the aspect ratio, sphericity and circularity reach a steady value at large strains, in parallel to reaching a stable grading. The mobilized angle of shearing resistance however was not significantly different in the different types of samples, suggesting the final grading dominates the behaviour.     

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.     

Role of particle shape on the shear strength of sand-GCL interfaces under dry and wet conditions

Interface shear strength of geosynthetic clay liners (GCL) with the sand particles is predominantly influenced by the surface characteristics of the GCL, size and shape of the sand particles and their interaction mechanisms. This study brings out the quantitative effects of particle shape on the interaction mechanisms and shear strength of GCL-sand interfaces. Interface direct shear tests are conducted on GCL in contact with a natural sand and a manufactured sand of identical gradation, eliminating the particle size effects. Results showed that manufactured sand provides effective particle-fiber interlocking compared to river sand, due to the favorable shape of its grains. Further, the role of particle shape on the hydration of GCL is investigated through interface shear tests on GCL-sand interfaces at different water contents. Bentonite hydration is found to be less in tests with manufactured sand, leading to better interface shear strength. Grain shape parameters of sands, surface changes related to hydration and particle entrapment in GCL are quantified through image analysis on sands and tested GCL surfaces. It is observed that the manufactured sand provides higher interface shear strength and causes lesser hydration related damages to GCL, owing to its angular particles and low permeability.     

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.     

Effect of particle shape on the response of geogrid-reinforced systems: Insights from 3D discrete element analysis

Understanding soil-geogrid interaction is essential for the analysis and design of reinforced soil systems. Modeling this interaction requires proper consideration for the geogrid geometry and the particulate nature of the backfill soil. This is particularly true when angular soil particles (e.g. crushed limestone) are used as a backfill material. In this study, a three-dimensional (3D) discrete element model that is capable of capturing the response of unconfined and soil-confined geogrid material is developed and used to study the response of crushed limestone reinforced with geogrid and subjected to surface loading. The 3D shape of the crushed limestone is modeled by tracing the surface areas of a typical particle and fitting a number of bonded spheres into the generated surface. Model calibration is performed using triaxial tests to determine the microparameters that allow for the stress-strain behaviour of the backfill material to be replicated. To demonstrate the role of particle shape on the soil-geogrid interaction, the analysis is also performed using spherical particles and the calculated response is compared with that obtained using modeled surfaces. The biaxial geogrid used in this study is also modeled using the discrete element method and the unconfined response is compared with the available index test results. This study suggests that modeling the 3D geogrid geometry is important to accurately capture the geogrid response under both confined and unconfined conditions. Accounting for the particle shape in the analysis can significantly enhance the predicted response of the geogrid-soil system. The modeling approach proposed in this study can be adapted for other reinforced soil applications.     

Particle loss: An initial investigation into size effects and stress-dilatancy

Drained triaxial tests have been performed to explore the effect of particle loss on shearing behaviour and critical states in granular mixtures. The mixtures comprise Leighton Buzzard sand (d₅₀ = 0.8 mm), to which was added 15% by mass of salt particles of different nominal sizes: 0.063 mm, 0.25 mm and 0.5 mm. Shearing behaviours before and after particle loss (by dissolution) were compared. A good fit is observed between the test data and a stress-dilatancy relationship for the post-dissolution tests, highlighting the ability of the stress-dilatancy analysis as a means to interpret the effects of particle loss on shearing. It was noted that critical state strength parameter M is determined by the post-dissolution grading regardless of size of removed particle. However, the duration of contractant volumetric strain increased with the larger removed particles (0.25 mm & 0.5 mm) even when initial specific volumes were virtually identical. It is suggested that a loose volumetric state is reached if the sand particle network is initially disrupted by the amount and/or size of salt particles, which following dissolution results in structural or fabric phenomena that are not reflected in scalar volumetric measures such as specific volume.     

Liquefaction and post-liquefaction resistance of sand reinforced with recycled geofibre

The present study provides an insight into the effect of recycled carpet fibre on the mechanical response of clean sand as backfill material subjected to monotonic loading and cyclic loading as well as post-liquefaction resistance of both unreinforced and carpet fibre reinforced soils. To achieve these goals, a series of multi-stage soil element tests under cyclic loading event resulting in liquefaction followed by undrained monotonic shearing without excess pore water pressure dissipation as well as a series of monotonic undrained shear test is conducted. All the specimens are isotropically consolidated under a constant effective confining stress of 100 kPa by considering the effect of cyclic stress ratio and carpet fibre content ranging from 0.25% to 0.75%. The obtained results revealed the efficiency of carpet fibre inclusion in increasing the secant shear modulus and ductility of clean sand under monotonic shearing without previous loading history. The impact of carpet fibre inclusion on the trend of cyclic excess pore water pressure generation and cyclic stiffness degradation was minimal. However, adding carpet fibre significantly improved both liquefaction and post-liquefaction resistances of clean sand. The liquefaction resistance of clean sand, at a constant 15 loading cycles, improved by 26.3% when the soil was reinforced with 0.75% recycled carpet fibre. In addition, the initial shear modulus of the liquefied specimen significantly increased by adding recycled carpet fibre.     

级配碎石形状量化及其对堆积特性的影响

针对现有离散元模拟多采用规则形状颗粒进行,与颗粒实际形态不符的缺陷。以级配碎石为例,对颗粒形状进行量化分析。采用数字图像处理技术,得到颗粒轮廓,根据轮廓形态特征从多个层次对颗粒形态进行描述,并得到互相独立的扁平度、棱角度形状指标,对指标进行统计分析,发现颗粒形态与粒径无关。最后通过离散元数值模拟探究颗粒真实形态与颗粒堆积孔隙率关系。试验结果表明：在松散情况下,堆积孔隙率随扁平度的增大而增大,而当颗粒近似圆形时,孔隙率随棱角度增大而减小。     

Centrifuge study into the effect of liquefaction extent on permanent settlement and seismic response of shallow foundations

Centrifuge experiments were conducted to investigate how the liquefaction extent affects the seismic and post-seismic settlement of shallow foundations resting on saturated sand. Two rigid foundations with different bearing pressures were placed on the ground surface in a model container. Multiple input motions were applied to achieve different extents of soil liquefaction. The results indicate that foundation settlement can be divided into three distinct phases: (I) during shaking, (II) during the time period after shaking has ceased and before soil reconsolidation in the shallowest layers has taken place, and (III) during soil reconsolidation. Contrary to the free-field ground, most of the total settlement of the foundations occurred before soil reconsolidation, i.e., during Phases I and II. The volumetric strain during these phases was not significant as opposed to the shear strain produced by the foundation surcharge. It was demonstrated that foundation settlement is not necessarily proportional to the liquefied depth of the sand. The extent of the liquefaction in the sand medium mostly affected the post-seismic settlement of the foundations, while the co-seismic settlement was relatively the same for both foundations. The response of the foundations was significantly influenced by the liquefaction extent, whereas the foundations did not experience large accelerations when the soil profile was entirely liquefied. However, the foundations tolerated large settlement under severe liquefaction conditions. The results of this study highlight the role of the liquefaction extent on co-seismic and post-seismic settlement as well as the seismic response of shallow foundations.