Stress-strain response of hydrate-bearing sands: Numerical study using discrete element method simulations

Gas hydrate is a crystalline solid found within marine and subpermafrost sediments. While the presence of hydrates can have a profound effect on sediment properties, the stress-strain behavior of hydrate-bearing sediments is poorly understood due to inherent limitations in laboratory testing. In this study, we use numerical simulations to improve our understanding of the mechanical behavior of hydrate-bearing sands. The hydrate mass is simulated as either small randomly distributed bonded grains or as "ripened hydrate" forming patchy saturation, whereby sediment clusters with 100% pore-filled hydrate saturation are distributed within a hydrate-free sediment. Simulation results reveal that reduced sand porosity and higher hydrate saturation cause an increase in stiffness, strength, and dilative tendency, and the critical state line shifts toward higher void ratio and higher shear strength. In particular, the critical state friction angle increases in sands with patchy saturation, while the apparent cohesion is affected the most when the hydrate mass is distributed in pores. Sediments with patchy hydrate distribution exhibit a slightly lower strength than sediments with randomly distributed hydrate. Finally, hydrate dissociation under drained conditions leads to volume contraction and/or stress relaxation, and pronounced shear strains can develop if the hydrate-bearing sand is subjected to deviatoric loading during dissociation.

The efficiency of compensation grouting in sands

Physical model experiments on compensation grouting in sands were performed in two different setups (Cambridge and Delft). The effect of water-cement (w/c) ratio, bentonite content (b.c.) and injection rate on compensation efficiency was investigated. Results show a considerable drop in compensation efficiency resulted from reducing the soil density. Injection in dense sand (R.D. = 93%) resulted in efficiencies between 40-90%, whereas injection in medium-dense sand (R.D. = 60-75%) yielded in reduced efficiencies between 10-40%. When the w/c ratio increased from 0.5 to 1.5 for a given density (R.D. = 93%) and the b.c. of 4%, the compensation efficiency value decreased. Typical efficiencies were between 60% and 40-50% for w/c ratios of 0.5 and 1.5, respectively. The values of compensation and grout efficiencies were almost equal, suggesting that pressure filtration happens mainly during injection. Increasing the b.c. improved the compensation efficiency. When a higher b.c. of 12% to 14% was used, typical compensation efficiencies in dense sand were 78 and 90% for w/c ratios of 1.5 and 1.8 respectively. © 2012 Taylor & Francis Group.

Stiffness of sands through a laboratory test database

Deformations of sandy soils around geotechnical structures generally involve strains in the range small (0·01%) to medium (0·5%). In this strain range the soil exhibits non-linear stress-strain behaviour, which should be incorporated in any deformation analysis. In order to capture the possible variability in the non-linear behaviour of various sands, a database was constructed including the secant shear modulus degradation curves of 454 tests from the literature. By obtaining a unique S-shaped curve of shear modulus degradation, a modified hyperbolic relationship was fitted. The three curve-fitting parameters are: an elastic threshold strain γe, up to which the elastic shear modulus is effectively constant at G0; a reference strain γr, defined as the shear strain at which the secant modulus has reduced to 0·5G0; and a curvature parameter a, which controls the rate of modulus reduction. The two characteristic strains γe and γr were found to vary with sand type (i.e. uniformity coefficient), soil state (i.e. void ratio, relative density) and mean effective stress. The new empirical expression for shear modulus reduction G/G0 is shown to make predictions that are accurate within a factor of 1·13 for one standard deviation of random error, as determined from 3860 data points. The initial elastic shear modulus, G0, should always be measured if possible, but a new empirical relation is shown to provide estimates within a factor of 1·6 for one standard deviation of random error, as determined from 379 tests. The new expressions for non-linear deformation are easy to apply in practice, and should be useful in the analysis of geotechnical structures under static loading.

BEHAVIOUR OF FIXED CANTILEVER WALLS SUBJECT TO LATERAL SHAKING.

A programme of research on the seismic behaviour of retaining walls has been under way at Cambridge since 1981. Centrifuge tests have presently been conducted both on cantilever walls and isolated mass walls, retaining dry sands of varying grading and density. This paper is devoted to the modelling of fixed-base cantilever walls retaining Leighton Buzzard (14/25) sand of relative density 99% with a horizontal surface level with the crest of the wall. The base of the centrifuge container was used to fix the walls, and to provide a rigid lower boundary for the sand. No attempt was made to inhibit the propagation of compression waves from the side of the container opposite the inside face of the model wall. The detailed analysis of dynamic deflections and bending moments was made difficult by the anelastic nature of reinforced concrete, and the difficulty of measuring bending strains thereon. A supplementary programme of well-instrumented tests on Dural walls of similar stiffness, including the modelling of models, was therefore carried out. Refs.

Measurement of coefficient of consolidation during reconsolidation of liquefied sand

Saturated sands particularly at low relative density commonly exhibit rises in excess pore pressure when subjected to earthquake loading. The excess pore pressure can approach a maximum value, limited by the initial vertical effective stress. After the completion of earthquake shaking, these excess pore pressures dissipate according to the consolidation equation, which can be solved to produce a Fourier series solution. It will be shown by manipulation of this Fourier series that excess pore pressure traces provide a method for back-calculation of coefficient of consolidation Cv. This method is validated against dissipation curves generated using known values of C v and seen to be more accurate in the middle of the layer. The method is then applied to data recorded in centrifuge tests to evaluate Cv throughout the reconsolidation process following liquefaction conditions. C v is seen to fit better as a function of excess pore pressure ratio than effective stress for the stress levels considered. For the soil investigated, Cv is about three times smaller at excess pore pressure ratio of 0.9 compared to excess pore pressure ratio of 0. Copyright © 1996-2011 ASTM.

3D FE Analyses of Buried Pipeline with Elbows Subjected to Lateral Loading

This study investigates the interaction between soil and pipeline in sand subjected to lateral ground displacements with emphasis on the peak force exerted to a bended elbow-pipe. A series of three-dimensional (3D) finite-element (FE) analyses were performed in both opening and closing modes of the elbow section for different initial pipe bending angles. To model the mechanical behavior of sands, two soil models were adopted: Mohr-Coulomb and Nor-Sand soil model. Investigations also included the effects of pipe embedment depth and soil density. Results show that the opening mode exhibits higher ultimate forces and greater localized deformations than the closing mode. Nondimensional charts that account for pipeline location, bending angle, and soil density are developed. Soil-spring pipeline analyses of an elbow-pipe were performed using modified F-δ soil-spring models based on the 3D FE results and were compared to the findings of conventional spring model analyses using the standard two-dimensional soil-spring model. Results show that the pipe strain does not change in the closing mode case. However, in the opening mode case, the pipe strain computed by the modified analysis is larger than that by the conventional analysis and the difference is more pronounced when the pipe stiffness is stiffer. © 2011 American Society of Civil Engineers.