Laboratory Investigation on Assessing Liquefaction Resistance of Sandy Soils by Shear Wave Velocity

Shear wave velocity (Vs) offers engineers a promising alternative tool to evaluate liquefaction resistance of sandy soils, and the lack of sufficient in-situ databases makes controlled laboratory study very important. In this study, semitheoretical considerations were first given based on review of previous liquefaction studies, which predicted a possible relationship between laboratory cyclic resistance ratio (CRRtx) and Vs normalized with respect to the minimum void ratio, confining stress and exponent n of Hardin equation. Undrained cyclic triaxial tests were then performed on three reconstituted sands with Vs measured by bender elements, which verified this soil-type-dependent relationship. Further investigation on similar laboratory studies resulted in a database of 291 sets of data from 34 types of sandy soils, based on which the correlation between liquefaction resistance and Vs was established statistically and further converted to equivalent field conditions with well-defined parameters, revealing that CRR will vary proportionally with (Vs1)4. Detailed comparisons with Vs-based site-specific investigations show that the present lower-bound CRR–Vs1 curve is a reliable prediction especially for sites with higher CSR or Vs1. The framework of liquefaction assessment based on the present laboratory study is proposed for engineering practice.     

Influence of Nonplastic Fines on Shear Wave Velocity-Based Assessment of Liquefaction

Many false positives (no liquefaction detected when the normalized shear wave velocity-cyclic stress ratio (Vs1-CSR) combination indicated that it should have been) are observed in the database used in the simplified liquefaction assessment procedure based on shear wave velocity. Two possible reasons for false positives are the presence of a thick surface layer of nonliquefiable soil and the effects of fines on cyclic shear resistance (CRR) and Vs1. About 67% of the false positives that could not have been caused by an overlying thick surface layer are associated with silty sands with less than 35% fines. The effects of fines on the liquefaction resistance of silty sands and on the shear wave velocity are analyzed. Theoretical CRRfield?versus?Vs1 curves for silty sands containing 0 to 15% nonplastic fines are established. They show that the theoretical CRR-Vs1 correlations for silty sands with 5 to 15% nonplastic fines are all located to the far left of the semi-empirical curves that separate liquefaction from no-liquefaction zones in the simplified liquefaction potential assessment procedures. The results suggest the currently used shear wave velocity-based liquefaction potential curves may be overly conservative when applied to sands containing nonplastic fines.     

Liquefaction Resistance of Clean and Nonplastic Silty Sands Based on Cone Penetration Resistance

Liquefaction of granular soil deposits is one of the major causes of loss resulting from earthquakes. The accuracy in the assessment of the likelihood of liquefaction at a site affects the safety and economy of the design. In this paper, curves of cyclic resistance ratio (CRR) versus cone penetration test (CPT) stress-normalized cone resistance qc1 are developed from a combination of analysis and laboratory testing. The approach consists of two steps: (1) determination of the CRR as a function of relative density from cyclic triaxial tests performed on samples isotropically consolidated to 100 kPa; and (2) estimation of the stress-normalized cone resistance qc1 for the relative densities at which the soil liquefaction tests were performed. A well-tested penetration resistance analysis based on cavity expansion analysis was used to calculate qc1 for the various soil densities. A set of 64 cyclic triaxial tests were performed on specimens of Ottawa sand with nonplastic silt content in the range of 0–15% by weight, and relative densities from loose to dense for each gradation, to establish the relationship of the CRR to the soil state and fines content. The resulting (CRR)7.5-qc1 relationship for clean sand is consistent with widely accepted empirical relationships. The (CRR)7.5-qc1 relationships for the silty sands depend on the relative effect of silt content on the CRR and qc1. It is shown that the cone resistance increases at a higher rate with increasing silt content than does liquefaction resistance, shifting the (CRR)7.5-qc1 curves to the right. The (CRR)7.5-qc1 curves proposed for both clean and silty sands are consistent with field observations.     

Accounting for Soil Aging When Assessing Liquefaction Potential

It has been recognized that liquefaction resistance of sand increases with age due to processes such as cementation at particle contacts and increasing frictional resistance resulting from particle rearrangement and interlocking. As such, the currently available empirical correlations derived from liquefaction of young Holocene sand deposits, and used to determine liquefaction resistance of sand deposits from in situ soil indices [standard penetration test (SPT), cone penetration test (CPT), shear wave velocity test (Vs)], are not applicable for old sand deposits. To overcome this limitation, a methodology was developed to account for the effect of aging on the liquefaction resistance of old sand deposits. The methodology is based upon the currently existing empirical boundary curves for Holocene age soils and utilizes correction factors presented in the literature that comprise the effect of aging on the in situ soil indices as well as on the field cyclic strength (CRR). This paper describes how to combine currently recorded SPT, CPT, and Vs values with corresponding CRR values derived for aged soil deposits to generate new empirical boundary curves for aged soils. The method is illustrated using existing geotechnical data from four sites in the South Carolina Coastal Plain (SCCP) where sand boils associated with prehistoric earthquakes have been found. These sites involve sand deposits that are 200,000?to?450,000?years in age. This work shows that accounting for aging of soils in the SCCP yields less conservative results regarding the current liquefaction potential than when age is not considered. The modified boundary curves indicate that old sand deposits are more resistant to liquefaction than indicated by the existing empirical curves and can be used to evaluate the liquefaction potential at a specific site directly from the current in situ properties of the soil.     

Updated Liquefaction Resistance Correction Factors for Aged Sands

Data from over 30 sites in 5 countries are analyzed to develop updated factors for correcting liquefaction resistance for aged sand deposits. Results of cyclic laboratory tests on relatively undisturbed and reconstituted specimens suggest an increase in the correction factors of 0.12 per log cycle of time and an average reference age of 2 days for the reconstitute specimens. Laboratory and field test results combined with cyclic resistance ratio (CRR) charts suggest an increase in the correction factors of 0.13 per log cycle of time and an average reference age of 23 years. A reference age of 23 years seems appropriate for the commonly used CRR charts derived from field liquefaction and no liquefaction case history data. Because age of natural deposits is often difficult to accurately determine, a relationship between measured to estimated shear-wave velocity ratio (MEVR) and liquefaction resistance correction factor is also derived directly from the compiled data. This new MEVR-liquefaction resistance correction factor relationship is not as sensitive to MEVR as in the relationship derived indirectly in a previous paper.     

Coupled Use of Cone Tip Resistance and Small Strain Shear Modulus to Assess Liquefaction Potential

Resistance against earthquake-related liquefaction is usually assessed using relationships between an index of soil strength such as normalized cone tip resistance and the cyclic resistance ratio (CRR) developed from observed field performance. The alternative approach based on laboratory testing is rarely used, mainly because of the apprehension that laboratory results may not reflect field behavior since the quality of laboratory data is often compromised by sampling disturbance. In this study, a database of laboratory data obtained mainly from cyclic testing of frozen (undisturbed) samples and in situ index measurements from near sampling locations comprised of cone tip resistance, qc, and shear wave velocity, Vs, have been assembled. These data indicate that neither normalized cone tip resistance nor normalized shear wave velocity individually correlate well with laboratory-measured CRR. However, the ratio of qc to the small strain shear modulus, G0, relates reasonably with CRR via separate correlations depending on geologic age. The derived qc/G0-CRR relationships were also found to be consistent with earthquake field-performance case histories.     

Use of Artificial Neural Networks in the Prediction of Liquefaction Resistance of Sands

A backpropagation artificial neural network (ANN) model has been developed to predict the liquefaction cyclic resistance ratio (CRR) of sands using data from several laboratory studies involving undrained cyclic triaxial and cyclic simple shear testing. The model was verified using data that was not used for training as well as a set of independent data available from laboratory cyclic shear tests on another soil. The observed agreement between the predictions and the measured CRR values indicate that the model is capable of effectively capturing the liquefaction resistance of a number of sands under varying initial conditions. The predicted CRR values are mostly sensitive to the variations in relative density thus confirming the ability of the model to mimic the dominant dependence of liquefaction susceptibility on soil density already known from field and experimental observations. Although it is common to use mechanics-based approaches to understand fundamental soil response, the results clearly demonstrate that non-mechanistic ANN modeling also has a strong potential in the prediction of complex phenomena such as liquefaction resistance.     

CPT-Based Evaluation of Liquefaction and Lateral Spreading in Centrifuge

Two series of centrifuge model tests were conducted using Nevada sand. Four saturated models placed in a mildly inclined laminar box and simulating a 6-m-thick deposit were shaken inducing liquefaction effects and lateral spreading. The sand was deposited at a relative density, Dr = 45 or 75%; two of the 45% models were subjected to overconsolidation or preshaking. The second series involved in-flight measurements of static cone tip penetration resistance, qc, simulating the standard cone penetration test (CPT) 36-mm cone. Values of qc increased with Dr, overconsolidation, and preshaking. A normalized resistance, qc1N, was assigned to each of the four liquefaction/lateral spreading models. Increases in Dr, overconsolidation, and preshaking decreased liquefaction and ground deformation, but relative density alone captured these effects rather poorly. Conversely, qc1N predicted extremely well the liquefaction and lateral spreading response of the four models, confirming Seed’s hypothesis to explain the success of penetration-based seismic liquefaction charts. The depth to liquefaction measured in the four centrifuge models is consistent with the field CPT liquefaction chart.     

Liquefaction Resistance of Sandy Soils under Partially Drained Condition

In order to simulate the effect of drainage on soils adjacent to gravel drains that are installed as countermeasure against liquefaction, several series of cyclic triaxial tests were performed on saturated sands under partially drained conditions. The condition of partial drainage under cyclic loading was simulated in the laboratory using triaxial testing equipment installed with a drainage control valve to precisely regulate the volume of water being drained from test specimens. Effects of both drainage conditions and loading frequencies on cyclic response were incorporated through the coefficient of drainage effect, α*. Experimental results showed that for sand exhibiting strain softening, the partially drained response was controlled by the critical effective stress ratio while for sand showing strain hardening behavior, the controlling factor was the phase transformation stress ratio. Moreover, test results indicated that the minimum liquefaction resistance under partially drained conditions can be used as a parameter to describe the liquefaction resistance of sands improved by the gravel drain method. From these results, a simplified procedure for designing gravel drains based on the factor of safety (FL) concept was proposed.     

Correlation between Cyclic Resistance and Shear-Wave Velocity for Providence Silts

As an alternative to a field-based liquefaction resistance approach, cyclic triaxial tests with bender elements were used to develop a new correlation between cyclic resistance ratio (CRR) and overburden stress-corrected shear-wave velocity (VS1) for two nonplastic silts obtained from Providence, Rhode Island. Samples of natural nonplastic silt were recovered by block sampling and from geotechnical borings/split-spoon sampling. The data show that the correlation is independent of the soils’ stress history as well as the method used to prepare the silt for cyclic testing. The laboratory results indicate that using the existing field-based CRR-VS1 correlations will significantly overestimate the cyclic resistance of the Providence silts. The strong dependency of the CRR-VS1 curves on soil type also suggests the necessity of developing silt-specific liquefaction resistance curves from laboratory cyclic tests performed on reconstituted samples.     

Saturation and Preloading Effects on the Cyclic Behavior of Sand

In order to study pore water pressure response and liquefaction characteristics of sand, which has previously experienced liquefaction, two series of cyclic triaxial tests were run on medium dense sand specimens. In the first test series the influence of the soil saturation under undrained cyclic loading has been studied. It summarizes results of cyclic triaxial tests performed on Hostun-RF sand at various values of the Skempton’s pore-pressure coefficient. Analysis of experimental results gives valuable insights on the effect of soil saturation on sand response to undrained cyclic paths. In the second series of tests, the preloading influence on the resistance to the sands liquefaction has been realized on samples at various histories of loading. It was found that a large preloading induces a reduction of the resistance of sands to liquefaction.     

Characterization of Liquefaction Susceptibility of Sands by Means of Extreme Void Ratios and/or Void Ratio Range

The liquefaction susceptibility of various graded fine to medium saturated sands are evaluated by stress controlled cyclic triaxial laboratory tests. Cyclic triaxial tests are performed on reconstituted specimens having global relative density of 60%. In all cyclic triaxial tests, loading pattern is selected as a sinusoidal wave form with 1.0 Hz frequency and effective consolidation pressure is chosen as 100 kPa. Liquefaction resistance is defined as the required cyclic stress ratio causing initial liquefaction in 10 cycles during the cyclic triaxial test. The results are used to draw conclusions on the effect of the extreme void ratios and void ratio range on the liquefaction resistance of various graded sands.     

Liquefaction-Induced Softening of Load Transfer between Pile Groups and Laterally Spreading Crusts

Laterally spreading nonliquefied crusts can exert large loads on pile foundations causing major damage to structures. While monotonic load tests of pile caps indicate that full passive resistance may be mobilized by displacements on the order of 1–7% of the pile cap height, dynamic centrifuge model tests show that much larger relative displacements may be required to mobilize the full passive load from a laterally spreading crust onto a pile group. The centrifuge models contained six-pile groups embedded in a gently sloping soil profile with a nonliquefied crust over liquefiable loose sand over dense sand. The nonliquefied crust layer spread downslope on top of the liquefied sand layer, and failed in the passive mode against the pile foundations. The dynamic trace of lateral load versus relative displacement between the “free-field” crust and pile cap is nonlinear and hysteretic, and depends on the cyclic mobility of the underlying liquefiable sand, ground motion characteristics, and cyclic degradation and cracking of the nonliquefied crust. Analytical models are derived to explain a mechanism by which liquefaction of the underlying sand layer causes the soil-to-pile-cap interaction stresses to be distributed through a larger zone of influence in the crust, thereby contributing to the softer load transfer behavior. The analytical models distinguish between structural loading and lateral spreading conditions. Load transfer relations obtained from the two analytical models reasonably envelope the responses observed in the centrifuge tests.     

Degree of Saturation and Liquefaction Resistances of Sand Improved with Sand Compaction Pile

Sand compaction pile (SCP) is a ground improvement technique extensively used to ameliorate liquefaction resistance of loose sand deposits. This paper discusses results of laboratory tests on high-quality undisturbed samples obtained by the in situ freezing method at six sites where foundation soils had been improved with SCP. Inspection of samples revealed that the improved ground was desaturated during the ground improvement. Degree of saturation (Sr) was lower than 77% for the sand piles and 91% for the improved sand layers, while Sr was approximately 100% for improved clayey and silty soils. A good correlation was found between Sr and 5% diameter of the soil; the larger 5% diameter of soils (D5), the lower the degree of saturation. It appeared that the variation of Sr with D5 for soils within a month after the ground improvement work was quite similar in trend to that after more than several years. Degree of saturation of soils after several years was noticeably, but not significantly, higher as compared with that shortly after ground improvement, indicating longevity of air bubbles injected in the improved soil. Undrained cyclic shear tests were also carried out on saturated and unsaturated specimens and effects of desaturation on undrained cyclic shear strength were studied. The test results were summarized in a form of liquefaction resistance with reference to normalized standard penetration test N-value.     

Effects of Nonplastic Fines on the Liquefaction Resistance of Sands

A laboratory parametric study utilizing cyclic triaxial tests was performed to clarify the effects of nonplastic fines on the liquefaction susceptibility of sands. Studies previously published in the literature have reported what appear to be conflicting results as to the effects of silt content on the liquefaction susceptibility of sandy soils. The current study has shown that if the soil structure is composed of silt particles contained within a sand matrix, the resistance to liquefaction of the soil is controlled by the relative density of the soil and is independent of the silt content of the soil. For soils whose structure is composed of sand particles suspended within a silt matrix, the resistance to liquefaction is again controlled by the relative density of the soil, but is lower than for soils with sand-dominated matrices at similar relative densities. In this case, the resistance to liquefaction is essentially independent of the amount and type of sand. These findings suggest the need for further evaluation of the effects of nonplastic fines content upon penetration resistance, and the manner in which this relationship affects the simplified methods currently used in engineering practice to evaluate the liquefaction resistance of silty soils.     

Development of Dissipation Model of Excess Pore Pressure in Liquefied Sandy Ground

Recently, several studies on the dissipation of excess pore pressure in liquefied sandy grounds have been reported on evaluation of postliquefaction behavior of structures. To further contribute to the understanding of this complicated dynamic phenomenon, centrifuge tests were performed in this research to analyze the liquefaction behavior of level saturated sandy grounds. The test results showed that the excess pore pressure in the liquefied sand was dissipated by the combined process of the solidification of the sand grains and the consolidation of the solidified layer. Based on the test results, a nonlinear model for the solidified layer thickness versus time, i.e., the solidification velocity, was developed. A new dissipation model was also developed by combining the nonlinear solidification model with Scott’s theory to improve the prediction for the time history of excess pore pressure. In addition, a method for evaluating the input parameters in the dissipation model was proposed by relating the parameters to the particle size and the relative density of soils. The proposed dissipation model properly estimated the dissipation of excess pore pressure with time observed in the centrifuge experiments.     

Method for Estimating Dynamic Compaction Effect on Sand

This paper proposes a method for estimating the degree and depth of improvement resulting from dynamic compaction on sand. The method is based on extensive two-dimensional finite element analyses, which are benchmarked using centrifuge model results. By appropriately normalizing the results from a wide range of soil properties, initial conditions and operating parameters such as momentum per blow, energy per blow, and number of blows, the results for numerous cases can be summarized into a relatively small number of plots, which can be used for predictive purposes. A preliminary assessment of the reliability of these curves was conducted using results from a centrifuge and a field study. The results of this comparative exercise indicate that the predictions using the method are reasonably realistic.     

DSC Model for Soil and Interface Including Liquefaction and Prediction of Centrifuge Test

Realistic predictions of dynamic soil–structure interaction problems require appropriate constitutive models for the characterization of soils and interfaces. This paper presents a unified model based on the disturbed state concept (DSC). The parameters for the models for the Nevada sand, and sand–metal interface are obtained based on available triaxial test data on the sand and interfaces. The predicted stress–strain–pore water pressure behavior for the sand using the DSC model is compared with the test data. In addition, a finite element procedure with the DSC model, based on the generalized Biot’s theory, is used to predict the measured responses for a pile (aluminum) sand foundation problem obtained by using the centrifuge test. The predictions compared very well with measured pore water pressures. The DSC model is used to identify microstructural instability leading to liquefaction. A procedure is proposed to apply the proposed method for analysis and design for dynamic response and liquefaction.     

Observed Seismic Lateral Resistance of Liquefying Sand

This paper presents results from a study of the dynamic response of pile foundations in liquefying sand during seismic loading. The study included a series of dynamic centrifuge tests of pile-supported structures and the back-calculation of time histories for the lateral resistance p and relative displacement y between a pile and the free-field soil. Details of the centrifuge experiments and the procedures used to back-calculate p and y time histories are described. The back-calculated p-y time histories provide a concise representation of the experimental results and can be compared to the equivalent p-y behavior predicted by soil-pile interaction analysis methods. The observed p-y behavior provides insight into the mechanisms of soil-pile interaction in liquefying sand, showing characteristics that are consistent with the undrained cyclic loading behavior of saturated sand, including the effects of relative density, cyclic degradation, pore-pressure generation, prior displacement (strain) history, and phase transformation behavior.     

Dynamic and Static Load Testing of Model Piles Driven into Dense Sand

Systematic studies of pile response during static and dynamic load tests are generally too expensive to conduct in the field, and at model scale may be limited by scaling effects and the ability to obtain accurate stress waves. This paper describes such a study, conducted at model scale in a geotechnical centrifuge, for piles driven into dense sand. Adverse scaling effects were minimized by the use of extremely fine sand (silica flour), and accurate stress waves were obtained using high-frequency data logging, together with a Hopkinson bar arrangement for the measurement of pile-head velocity. The overall aim of the study was to compare dynamic and static test data, for open- and closed-ended piles driven into dense sand, for a range of delivered hammer energies. Open- and close-ended model piles were driven into dense sand and statically load tested at different penetrations, without stopping the centrifuge. Stress-wave data were collected, during continuous driving and from single blows immediately prior to the commencement of static load testing. An assessment of the accuracy of the mobilized soil resistance estimated from dynamic testing, using different levels of analysis, has been made by direct comparison with the static load test data. Particular emphasis has been placed on the performance of the dynamic analyses in light of different driving conditions and delivered hammer energy. From measurements of the delivered hammer energy, the efficiency of the centrifuge pile driving system was also assessed.