Post-liquefaction undrained monotonic behaviour of sands: experiments and DEM simulations

This paper presents the results of a series of servo-controlled cyclic triaxial tests and numerical simulations using the three- dimensional discrete element method (DEM) on post-liquefaction undrained monotonic strength of granular materials. In a first test series,undrained monotonic tests were carried out after dissipating the excess pore water pressure developed during liquefaction. The influence of different parameters such as amplitude of axial strain,relative density and confining pressure prior to liquefaction on the post-liquefaction undrained response have been investigated.The results obtained highlight an insignificant influence of amplitude of axial strain, confining pressure and a significant influence of relative density on the post-liquefaction undrained monotonic stress-strain response.In the second series, undrained monotonic tests were carried out on similar triaxial samples without dissipating the excess pore water pressure developed during liquefaction. The results highlight that the amplitude of axial strain prior to liquefaction has a significant influence on the post-liquefaction undrained monotonic response.In addition,DEM simulations have been carried out on an assembly of spheres to simulate post-liquefaction behaviour.The simulations were very similar to the experiments with an objective to understand the behaviour of monotonic strength of liquefied samples from the grain scale. The numerical simulations using DEM have captured qualitatively all the features of the post-liquefaction undrained monotonic response in a manner similar to that of the experiments.In addition,a detailed study on the evolution of micromechanical parameters such as the average coordination number and induced anisotropic coefficients has been reported during the post-liquefaction undrained monotonic loading.

Seismic horizontal pullout capacity of vertical anchors in sands

The problem of finding the horizontal pullout capacity of vertical anchors embedded in sands with the inclusion of pseudostatic horizontal earthquake body forces, was tackled in this note. The analysis was carried out using an upper bound limit analysis, with the consideration of two different collapse mechanisms: bilinear and composite logarithmic spiral rupture surfaces. The results are presented in nondimensional form to find the pullout resistance with changes in earthquake acceleration for different combinations of embedment ratio of the anchor (lambda), friction angle of the soil (phi), and the anchor-soil interface wall friction angle (delta). The pullout resistance decreases quite substantially with increases in the magnitude of the earthquake acceleration. For values of delta up to about 0.25-0.5phi, the bilinear and composite logarithmic spiral rupture surfaces gave almost identical answers, whereas for higher values of delta, the choice of the logarithmic spiral provides significantly smaller pullout resistance. The results compare favorably with the existing theoretical data.

Damping of Sands for Varying Saturation

A series of resonant column tests have been performed in the torsional mode of vibration to assess the effect of saturation, starting from the near dry state to the fully saturated state, on the damping ratio of sands corresponding to the threshold strain level. Tests were carried out on three different gradations of sand for various combinations of relative density and effective confining pressure. For fine sands, a certain optimum degree of saturation exists at which the damping ratio minimizes; it is known that a decrease in Sr from a fully saturated state leads to a continuous increase in the matric suction. With an increase in the relative density, the optimum degree of saturation for fine sand increases marginally from 1.38 to 1.49%, but does not show any dependency on the effective confining pressure. In contrast, the minimum values of the damping ratio for medium and coarse sands are found to always correspond to the near dry state. The damping ratio decreases continuously with increases in relative density and effective confining pressure. The threshold strain level has been found to decrease continuously with increases in relative density and effective confining pressure. (C) 2013 American Society of Civil Engineers.

Shaking table tests to investigate the influence of various factors on the liquefaction resistance of sands

This paper presents the shaking table studies to investigate the factors that influence the liquefaction resistance of sand. A uniaxial shaking table with a perspex model container was used for the model tests, and saturated sand beds were prepared using wet pluviation method. The models were subjected to horizontal base shaking, and the variation of pore water pressure was measured. Three series of tests varying the acceleration and frequency of base shaking and density of the soil were carried out on sand beds simulating free field condition. Liquefaction was visualized in some model tests, which was also established through pore water pressure ratios. Effective stress was calculated at the point of pore water pressure measurement, and the number of cycles required to liquefy the sand bed were estimated and matched with visual observations. It was observed that there was a gradual variation in pore water pressure with change in base acceleration at a given frequency of shaking. The variation in pore water pressure is not significant for the range of frequency used in the tests. The frequency of base shaking at which the sand starts to liquefy when the sand bed is subjected to any specific base acceleration depends on the density of sand, and it was observed that the sand does not liquefy at any other frequency less than this. A substantial improvement in liquefaction resistance of the sand was observed with the increase in soil density, inferring that soil densification is a simple technique that can be applied to increase the liquefaction resistance.

Experimental studies on the influence of intermediate principal stress and inclination on the mechanical behaviour of angular sands

A comprehensive experimental study has been made on angular sand to investigate various aspects of mechanical behavior. A hollow cylinder torsion testing apparatus is used in this program to apply a range of stress conditions on this angular quartzitic fine sand under monotonic drained shear. The effect of the magnitude and inclination of the principal stresses on an element of sand is studied through these experiments. This magnitude and inclination of the principal stresses are presented as an ``ensemble measure of fabric in sands''. This ensemble measure of fabric in the sands evolves through the shearing process, and reaches the final state, which indeed has a unique fabric. The sand shows significant variation in strength with changing inclination of the principal stresses. The locus of the final stress state in principal stress space is also mapped from these series of experiments. Additional aspects of non-coaxiality, a benchmarking exercise with a few constitutive models is presented here. This experimental approach albeit indirect shows that a unique state which is dependent on the fabric, density and confining stress exists. This suite of experiments provides a well-controlled data set for a clear understanding on the mechanical behavior of sands.

Seismic passive earth pressures in soils

Seismic passive earth pressure coefficients were computed by the method of limit equilibrium using a pseudostatic approach for seismic forces. Composite curved rupture surfaces were considered in the analysis. While earlier studies using this type of analysis were mainly for sands, seismic passive earth pressure coefficients were obtained in the present study considering the effects of cohesion, surcharge, and own weight. The minimum seismic passive force was obtained by adding the individual minimum values of these components and the validity of the principle of superposition was examined. Other parameters considered in the analysis were wall batter angle, ground surface slope, soil friction angle, wall friction angle, wall adhesion to soil cohesion ratio, and horizontal and vertical seismic accelerations. The seismic earth pressure coefficients were found to be highly sensitive to the seismic acceleration coefficients both in the horizontal and vertical directions. Results of the study are presented in the form of figures and tables. Comparisons of the proposed method with available theories in the seismic case are also presented.

Effect of relative density and confining pressure on Poisson ratio from bender and extender elements tests

By using the bender and extender elements tests, together with measurements of the travel times of shear (S) and primary (P) waves, the variation of Poisson ratio (nu) was determined for dry sands with respect to changes in relative densities and effective confining pressures (sigma(3)). The tests were performed for three different ranges of particle sizes. The magnitude of the Poisson ratio decreases invariably with an increase in both the relative density and the effective confining pressure. The effect of the confining pressure on the Poisson ratio was found to become relatively more significant for fine-grained sand as compared with the coarse-grained sand. For a given material, at a particular value of sigma(3), the magnitude of the Poisson ratio decreases, almost in a linear fashion, with an increase in the value of maximum shear modulus (G(max)). The two widely used correlations in literature, providing the relationships among G(max), void ratio (e) and effective confining pressure (sigma(3)), applicable for angular granular materials, were found to compare reasonably well with the present experimental data for the fine- and medium-grained sands. However, for the coarse-grained sand, these correlations tend to overestimate the values of G(max).

Support vector machine for evaluating seismic-liquefaction potential using shear wave velocity

The use of the shear wave velocity data as a field index for evaluating the liquefaction potential of sands is receiving increased attention because both shear wave velocity and liquefaction resistance are similarly influenced by many of the same factors such as void ratio, state of stress, stress history and geologic age. In this paper, the potential of support vector machine (SVM) based classification approach has been used to assess the liquefaction potential from actual shear wave velocity data. In this approach, an approximate implementation of a structural risk minimization (SRM) induction principle is done, which aims at minimizing a bound on the generalization error of a model rather than minimizing only the mean square error over the data set. Here SVM has been used as a classification tool to predict liquefaction potential of a soil based on shear wave velocity. The dataset consists the information of soil characteristics such as effective vertical stress (sigma'(v0)), soil type, shear wave velocity (V-s) and earthquake parameters such as peak horizontal acceleration (a(max)) and earthquake magnitude (M). Out of the available 186 datasets, 130 are considered for training and remaining 56 are used for testing the model. The study indicated that SVM can successfully model the complex relationship between seismic parameters, soil parameters and the liquefaction potential. In the model based on soil characteristics, the input parameters used are sigma'(v0), soil type. V-s, a(max) and M. In the other model based on shear wave velocity alone uses V-s, a(max) and M as input parameters. In this paper, it has been demonstrated that Vs alone can be used to predict the liquefaction potential of a soil using a support vector machine model. (C) 2010 Elsevier B.V. All rights reserved.

Seepage effects on sand-bed channels

ExperimentS were conducted to find the effects of seepage on flow over a sand bed in a straight rectangular flume under two conditions: (1) When the channel bed is plane, horizontal, and nontransporting; and (2) when the bed is transporting at a constant sediment concentration. Effects of both injection and suction, caused by seepage flow into and out of the channel bed, are studied for condition 1; and only suction effects are studied for condition 2, Three sands, sizes 0.34 mm, 0.53 mm, and 0.80 mm, are used in the study. It is found that seepage can cause an increase or decrease in the bed shear stress relative to no seepage for the two conditions. The change in bed shear stress depends on the relative magnitudes of the bed shear stress and the critical shear stress of particles under the no-seepage condition, sediment concentration, and the seepage rate. Quantitative relationships giving the ratio of bed shear stresses with and without seepage are presented for both conditions of the bed. A procedure to estimate the changes in bed shear stress, friction factor, Manning's n, and stream power due to seepage for known initial conditions of the channel and the amount of applied seepage is presented

Triaxial A-value versus swell or collapse for compacted soil - Discussion

A discussion of a technical note with the aforementioned title by Day and Marsh, published in this journal (Volume 121, Number 7, July 1995), is presented. Discussers Robinson and Allam assert that the authors' application of the pore-pressure parameter A to predict and quantify swell or collapse of compacted soils is hard to use because the authors visualize the collapse-swell phenomenon to occur in compacted soils broadly classified as sands and clays. The literature demonstrates that mineralogy has an important role in the volume change behavior of fine-grained soils. Robinson and Allam state that the A-value measurements may not completely predict the type of volume change anticipated in compacted soils on soaking without soil clay mineralogy details. Discussion is followed by closure from the authors.

Vertical Uplift Resistance of Pipes Buried in Sand

By incorporating the variation of peak soil friction angle (phi) with mean principal stress (sigma(m)), the effect of pipe diameter (D) on the vertical uplift resistance of a long horizontal pipeline embedded in sand has been investigated. The analysis has been performed by using the lower bound finite-element limit analysis in combination with nonlinear optimization. Three well-defined phi versus sigma(m) curves reported from literature for different sands have been used. It is observed that for a given embedment ratio, with an increase in pipe diameter, the magnitude of the uplift factor (F-gamma) reduces quite significantly, which indicates the importance of considering scale effects while designing buried pipe lines. The scale effects have been found to become even more substantial with an increase in the embedment ratio. The analysis compares well with various theoretical results reported from literature. On the other hand, as compared to available centrifuge test results, the present analysis has been found to provide quite a higher magnitude of the uplift resistance when the theoretical prediction is based on peak soil friction angle. However, if the theoretical analysis is performed by using the friction angle that accounts for the progressive shear failure, the difference between the theoretical and centrifuge test results decreases quite significantly.(C) 2013 American Society of Civil Engineers.

Uplift Resistance of Long Pipelines in the Presence of Seismic Forces

The vertical uplift resistance of long pipes buried in sands and subjected to pseudostatic seismic forces has been computed by using the lower-bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The soil mass is assumed to follow the Mohr-Coulomb failure criterion and an associated flow rule. The failure load is expressed in the form of a nondimensional uplift factor F-gamma. The variation of F-gamma is plotted as a function of the embedment ratio of the pipe, horizontal seismic acceleration coefficient (k(h)), and soil friction angle (phi). The magnitude of F-gamma is found to decrease continuously with an increase in the horizontal seismic acceleration coefficient. The reduction in the uplift resistance becomes quite significant, especially for greater values of embedment ratios and lower values of friction angle. The predicted uplift resistance was found to compare well with the existing results reported from the literature. (C) 2014 American Society of Civil Engineers.

Effect of Groundwater Seepage on Uplift Resistance of Buried Pipelines

The uplift resistance of pipelines buried in sands, in the presence of inclined groundwater flow, considering both upward and downward flow directions, has been determined by using the lower bound finite elements limit analysis in conjunction with nonlinear optimization. A correction factor (f (gamma) ), which needs to be multiplied with the uplift factor (F (gamma) ), has been computed to account for groundwater seepage. The variation of f (gamma) has been obtained as a function of i(gamma (w) /gamma (sub) ) for different horizontal inclinations (theta) of groundwater flow; where i = absolute magnitude of hydraulic gradient along the direction of flow, gamma (w) is the unit weight of water and gamma (sub) is the submerged unit weight of soil mass. For a given magnitude of i, there exists a certain critical value of theta for which the magnitude of f (gamma) becomes the minimum. An example has also been presented to illustrate the application of the results obtained for designing pipelines in presence of groundwater seepage.

Uplift Resistance of Interfering Pipelines Buried in Sand

The vertical uplift resistance of interfering pipelines buried in sands has been computed using the lower-bound limit analysis in conjunction with finite elements and nonlinear optimization. The soil mass is assumed to follow the Mohr-Coulomb failure criterion and an associated flow rule. It is specified that all the pipes fail simultaneously at the same magnitude of the failure load. For different clear spacing (S) between the pipes, the magnitude of the efficiency factor (xi(gamma)) is determined. Because of pipes' interference, with a reduction in the spacing between the pipelines, the magnitude of xi(gamma) is found to decrease continuously. The results were found to compare quite well with the available data from literature for horizontal strip anchors. (C) 2015 American Society of Civil Engineers.