Laboratory Study for Influence of Clay Content (CC) on Wave-Induced Liquefaction in Marine Sediments     

B. Liu 《Marine Georesources & Geotechnology》2016,34(3):280-292

In this artice, the influence of clay content on the wave-induced liquefaction in marine sediments was reported. The one-dimensional (1-D) equipment was setup with a vertical cylinder and 1.8–m–thick clayey sandy deposit and 0.2–m–thick water above the deposit. Unlike the previous experimental study for a single soil layer, this study used sand-kaolin mixtures, sand-illite mixtures, and sand-bentonite mixtures as the experimental samples. A series of experiments with 3,000 wave cycles in each test were conducted under numerous wave and soil conditions, which allowed us to examine the influence of clay content (CC) on wave-induced liquefaction in marine sediments. The experimental results showed that the clayey sandy deposit will become prone to liquefaction with the increase of CC when CC is less than a critical value, which depends on the type of clay. However, when CC is greater than the critical value, liquefaction depth will decrease as CC increases. Furthermore, when the CC value reaches a certain level, liquefaction will not occur. For example, no liquefaction occurs when CC ≥ 33% for both kaolin-sand and illite-sand mixtures and CC ≥ 16.36% for bentonite-sand mixtures.  相似文献   

Cavity expansion analyses of crushable granular materials with state‐dependent dilatancy     

M. J. Jiang  Y. G. Sun 《国际地质力学数值与分析法杂志》2012,36(6):723-742

Crushability is one of the important behaviors of granular materials particularly under high stress states, and affects both the deformability and strength of the materials that are in essence associated with state‐dependent dilatancy. In this presentation, first, a new critical state model is proposed to take into account the three different modes of compressive deformation of crushable granular materials, i.e. particle rearrangement, particle crushing and pseudo‐elastic deformation. Second, the governing equations for cavity expansion in crushable granulates are introduced, in which the state‐dependent dilatancy as well as the bounding surface plasticity model are used. Then, the procedure to obtain semi‐analytical solutions to cavity expansion in the material is described in detail, in which a commercial differential equation solver is employed. Finally, cavity expansion analyses are carried out on Toyoura sand, a well‐documented granular material, to demonstrate the effects of crushability and state‐dependent dilatancy. The study shows that particle crushing does occur at both high stress and critical states and affects the stress fields and the deformation behavior of the material surrounding the cavity in association with state‐dependent dilatancy. This leads to conclusion that particle crushing and state‐dependent dilatancy have to be taken into account when cavity expansion theory is used to interpret cone penetration tests and pressuremeter tests. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Explicit integration of a porosity‐dependent hydro‐mechanical model for unsaturated soils          下载免费PDF全文

Yue Zhang  Annan Zhou 《国际地质力学数值与分析法杂志》2016,40(17):2353-2382

An adaptive substepping explicit integration scheme is developed for a porosity‐dependent hydro‐mechanical model for unsaturated soils. The model is referred to as the modified σ –Θ model in this paper, which features the employment of the subloading surface plasticity and the stress–saturation approach. On numerical aspects, convex/nonconvex subloading surfaces in the σ –Θ space may result in incorrect loading–unloading decisions during the integration. A new loading–unloading decision method is developed here to solve the problem and then embedded into the explicit integration scheme for the modified σ –Θ model. In addition, to enhance the accuracy of the explicit integration, local errors from both hydraulic and mechanical components are included in the error control for each substep. A drift correction method is also developed to ensure the state point lies on the subloading surface in the σ –Θ space within a set error level. The performance of the loading–unloading decision method for the modified σ –Θ model is discussed through comparing it with the conventional loading–unloading decision method. The importance of involving the hydraulic component in the error control is also demonstrated. The accuracy and efficiency of the proposed adaptive substepping explicit integration scheme for the modified p–Θ model are also studied via several numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Bearing capacity of strip foundations reinforced with geogrid sheets by using upper bound finite‐element limit analysis     

Jyant Kumar  Jagdish Prasad Sahoo 《国际地质力学数值与分析法杂志》2013,37(18):3258-3277

By using the upper bound finite‐elements limit analysis, with an inclusion of single and two horizontal layers of reinforcements, the ultimate bearing capacity has been computed for a rigid strip footing placed over (i) fully granular, (ii) cohesive‐frictional, and (iii) fully cohesive soils. It is assumed that (i) the reinforcements are structurally strong so that no axial tension failure can occur, (ii) the reinforcement sheets have negligible resistance to bending, and (iii) the shear failure can take place between the reinforcement and soil mass. It is expected that the different approximations on which the analysis has been based would generally remain applicable for reinforcements in the form of geogrid sheets. A method has been proposed to incorporate the effect of the reinforcement in the analysis. The efficiency factors, η_c and η_γ, to be multiplied with N_c and N_γ , for finding the bearing capacity of reinforced foundations, have been established. The results have been obtained (i) for different values of ? in case of fully granular and cohesive‐frictional soils, and (ii) for different rates at which the cohesion increases with depth for a fully cohesive soil. The optimum positions of the reinforcements' layers have also been determined. The effect of the reinforcements' length on the results has also been analyzed. As compared to cohesive soils, the granular soils, especially with higher values of ?, cause a much greater increase in the bearing capacity. The results compare reasonably well with the available theoretical and experimental data from literature. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Stress‐dependent hardening and failure surfaces of dry sand     

W. Ehlers  O. Avci 《国际地质力学数值与分析法杂志》2013,37(8):787-809

During several triaxial compression experiments on plastic hardening, softening, and failure properties of dense sand specimens, it was found on various stress paths that the size of the failure surface was not constant. Instead, it changed depending on the current state of hydrostatic pressure. This finding is in contrast to the standard opinion consisting of the fact that the failure surface remains constant, once it has been reached during an experiment or in situ. In general, the behaviour of cohesionless granular‐material‐like sand is somehow characterised in between fluid and solid, where the solid behaviour results from the angle of internal friction and the confining pressure. Although the friction angle is an intrinsic material property, the confining pressure varies with the boundary conditions, thus defining different solid properties like plastic hardening, softening, and also failure. Based on our findings, it was the goal of the present contribution to introduce an improved setting for the plastic strain hardening and softening behaviour including the newly found yield properties at the limit state. For the identification of the material parameters, a complete triaxial experimental analysis of the tested sand is given. The overall elasto‐plasticity concept is validated by numerical computations of several laboratory foundation‐ and slope‐failure experiments. The performance of the proposed approach is compared with the standard concept of a constant failure surface, where the corresponding yield surfaces are understood as contours of equivalent plastic work or plastic strain. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Micromechanical viscoelasto‐plastic models and finite element implementation for rate‐independent and rate‐dependent permanent deformation of stone‐based materials     

Qingli Dai 《国际地质力学数值与分析法杂志》2010,34(13):1321-1345

This paper presents parallel and serial viscoelasto‐plastic models to simulate the rate‐independent and the rate‐dependent permanent deformation of stone‐based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto‐plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement‐based incremental formulation for the viscoelastic part and an elastic predictor—plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto‐plastic model employed a viscoelastic predictor—plastic corrector algorithm. The stone‐based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto‐plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto‐plastic matrix models. The simulation results showed that the serial viscoelasto‐plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto‐plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate‐independent and rate‐dependent viscoelasto‐plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto‐plastic behaviors of the stone‐based materials. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Description of deformation process in inherently anisotropic granular materials     

S. Pietruszczak  P. Guo 《国际地质力学数值与分析法杂志》2013,37(5):478-490

The main focus in this work is on modeling of mechanical response of granular materials that display inherent anisotropy. Both the experimental and numerical investigations are described. First, the results of direct shear as well as drained/undrained triaxial tests that involve crushed limestone with elongated angular‐shaped particles are reviewed. Afterward, a mathematical framework is presented for modeling of elastic/ inelastic deformation that incorporates the multi‐laminate approach. The deformation is monitored on a set of randomly oriented planes, and the formulation incorporates the thickness of the shear band that is associated with sliding/separation process. A systematic procedure for identification of material functions/ parameters is outlined that is based on the results of direct shear tests, and the framework is later applied to simulate the behavior under triaxial conditions. The results of numerical simulations are compared with the experimental data. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Garnet crystal plasticity in the continental crust,new example from south Madagascar     

J.‐E. MARTELAT  K. MALAMOUD  P. CORDIER  B. RANDRIANASOLO  K. SCHULMANN  J.‐M. LARDEAUX 《Journal of Metamorphic Geology》2012,30(4):435-452

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.  相似文献   

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation     

H.H.T. Tran  K.C. Toh  K.K. Phoon 《国际地质力学数值与分析法杂志》2013,37(17):2972-2986

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Explicit integration of bounding surface model for the analysis of earthquake soil liquefaction     

Konstantinos I. Andrianopoulos  Achilleas G. Papadimitriou  George D. Bouckovalas 《国际地质力学数值与分析法杂志》2010,34(15):1586-1614

This paper presents a new plasticity model developed for the simulation of monotonic and cyclic loading of non‐cohesive soils and its implementation to the commercial finite‐difference code FLAC, using its User‐Defined‐Model (UDM) capability. The new model incorporates the framework of Critical State Soil Mechanics, while it relies upon bounding surface plasticity with a vanished elastic region to simulate the non‐linear soil response. Stress integration of constitutive relations is performed using a recently proposed explicit scheme with automatic error control and substepping, which so far has been employed in the literature only for constitutive models aiming at monotonic loading. The overall accuracy of this scheme is evaluated at element level by simulating cyclic loading along complex stress paths and by using iso‐error maps for paths involving change of the Lode angle. The performance of the new constitutive model and its stress integration scheme in complex boundary value problems involving earthquake‐induced liquefaction is evaluated, in terms of accuracy and computational cost, via a number of parametric analyses inspired by the successful simulation of the VELACS centrifuge Model Test No. 2 studying the lateral spreading response of a liquefied sand layer. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献