初始组构各向异性对砂土力学特性及临界状态的影响

采用离散元方法,利用半径扩展法和重力沉积法分别生成具有初始各向同性和各向异性内结构的试样,并开展三轴不排水压缩和拉伸试验,研究不同制样方法产生的初始各向异性对砂土宏微观力学特性及其临界状态的影响。运用组构张量对砂土的各向异性进行量化,分析不同初始组构各向异性对组构张量演化的影响并确定了组构张量的临界值。试验结果表明：初始组构各向异性对试样的剪胀性有重要影响,由于受重力影响形成初始各向异性,其各向异性程度越大、组构方向与加载方向越一致,剪胀性越显著;初始组构各向异性对试样的临界状态没有影响,砂土的组构张量具有唯一的临界状态值。     

An AI-based model for describing cyclic characteristics of granular materials

Modelling cyclic behaviour of granular soils under both drained and undrained conditions with a good performance is still a challenge. This study presents a new way of modelling the cyclic behaviour of granular materials using deep learning. To capture the continuous cyclic behaviour in time dimension, the long short-term memory (LSTM) neural network is adopted, which is characterised by the prediction of sequential data, meaning that it provides a novel means of predicting the continuous behaviour of soils under various loading paths. Synthetic datasets of cyclic loading under drained and undrained conditions generated by an advanced soil constitutive model are first employed to explore an appropriate framework for the LSTM-based model. Then the LSTM-based model is used to estimate the cyclic behaviour of real sands, ie, the Toyoura sand under the undrained condition and the Fontainebleau sand under both undrained and drained conditions. The estimates are compared with actual experimental results, which indicates that the LSTM-based model can simultaneously simulate the cyclic behaviour of sand under both drained and undrained conditions, ie, (a) the cyclic mobility mechanism, the degradation of effective stress and large deformation under the undrained condition, and (b) shear strain accumulation and densification under the drained condition.     

Modelling the pre-failure instabilities of sand

A simple incremental model describing the pre-failure behaviour of granular soils is presented. The model describes both the dry/fully drained and undrained response. It takes into account an initial anisotropy of soil and an initial state defined as either contractive or dilative. A physically sound definition of loading/unloading is assumed, which differs from elasto-plastic approaches. The model is based on extensive empirical data and gives predictions conformable with experimental results. It also describes pre-failure instabilities of granular soils, both dry/fully drained and undrained. The Hill’s criterion was used to examine stability. It was shown that this condition can be formulated either in terms of the effective stresses or by the total stresses. In the extreme cases of either dry/fully drained or undrained conditions, these alternative formulations are equivalent. This is not so in the case of partial drainage of pore water and associated volumetric deformations as well as pore pressure changes. The model describes the pre-failure instabilities well, and additionally allows for analytical derivation of the instability line. It was shown that the second order work, appearing in the Hill’s condition, is equivalent to the entropy source.     

Numerical study of the coupled hydro-mechanical effects in dynamic compaction of saturated granular soils

Dynamic compaction is a widely used method for improvement of loose granular deposits. Its applicability in saturated layers generally considered to be less effective because of the fact that part of the applied energy is absorbed by pore water. Up to now the majority of numerical simulations have focused on the analysis of dynamic compaction in dry/moist soils. In this paper, a fully coupled hydro-mechanical finite element code has been developed and employed to evaluate the dynamic compaction effects on saturated granular soils. After verification of the results by comparing the numerical results with those measured in a real field case of DC treatment in a highway, some sensitivity analyses have been performed to evaluate the effect of water phase on the dimensions of the zone of improvement in the soil beneath the tamper. The results indicate that in the DC process the soil demonstrate two different behaviors. At the very early stage after impact, the soil behaves in an undrained manner and high oscillation of pore pressure occurs. After this phase, consolidation begins during which the pore-water-flow out of the soil mass takes place. The numerical analysis reveals that most of the DC improvement occurs during the undrained phase. The main mechanism responsible for the densification of soil during the undrained phase seems to be the compressibility of pore water. The simulation results indicate that the improvement zone diminishes when the degree of saturation increases.     

Predictive potential of Perzyna viscoplastic modelling for granular geomaterials

This paper reappraises Perzyna-type viscoplasticity for the constitutive modelling of granular geomaterials, with emphasis on the simulation of rate/time effects of different magnitude. An existing elasto-plastic model for sands is first recast into a Perzyna viscoplastic formulation and then calibrated/validated against laboratory test results on Hostun sand from the literature. Notable model features include (1) enhanced definition of the viscous nucleus function and (2) void ratio dependence of stiffness and viscous parameters, to model the pycnotropic behaviour of granular materials with a single set of parameters, uniquely identified against standard creep and triaxial test results. The comparison between experimental data and numerical simulations points out the predicative capability of the developed model and the complexity of defining a unique viscous nucleus function to capture sand behaviour under different loading/initial/boundary and drainage conditions. It is concluded that the unified viscoplastic simulation of both drained and undrained response is particularly challenging within Perzyna's framework and opens to future research in the area. The discussion presented is relevant, for instance, to the simulation of multiphase strain localisation phenomena, such as those associated to slope stability problems in variably saturated soils.     

Undrained expansion of a cylindrical cavity in clays with fabric anisotropy: theoretical solution

This paper presents a novel, exact, semi-analytical solution for the quasi-static undrained expansion of a cylindrical cavity in soft soils with fabric anisotropy. This is the first theoretical solution of the undrained expansion of a cylindrical cavity under plane strain conditions for soft soils with anisotropic behaviour of plastic nature. The solution is rigorously developed in detail, introducing a new stress invariant to deal with the soil fabric. The semi-analytical solution requires numerical evaluation of a system of six first-order ordinary differential equations. The results agree with finite element analyses and show the influence of anisotropic plastic behaviour. The effective stresses at critical state are constant, and they may be analytically related to the undrained shear strength. The initial vertical cross-anisotropy caused by soil deposition changes towards a radial cross-anisotropy after cavity expansion. The analysis of the stress paths shows that proper modelling of anisotropic plastic behaviour involves modelling not only the initial fabric anisotropy but also its evolution with plastic straining.     

Two-phase Material Point Method applied to the study of cone penetration

This paper presents numerical simulations of Cone Penetration Test (CPT) in water-saturated soft soils taking into account pore pressure dissipation during installation. Besides modelling interaction between soil skeleton and pore fluid, the problem involves large soil deformations in the vicinity of the penetrometer, soil–structure interaction, and complex non-linear response of soil. This makes such simulations challenging. Depending on the soil’s permeability and compressibility, undrained, partially drained or drained conditions might occur. Partially drained conditions are commonly encountered in soils such as silts and sand–clay mixtures. However, this is often neglected in CPT interpretation, which may lead to inaccurate estimates of soil properties. This paper aims at improving the understanding of the penetration process in different drainage conditions through advanced numerical analyses. A two-phase Material Point Method is applied to simulate large soil deformations and generation and dissipation of excess pore pressures during penetration. The constitutive behaviour of soil is modelled with the Modified Cam Clay model. Numerical results are compared with experimental data showing good agreement.     

An experimental database for the development,calibration and verification of constitutive models for sand with focus to cyclic loading: part I—tests with monotonic loading and stress cycles

For numerical studies of geotechnical structures under earthquake loading, aiming to examine a possible failure due to liquefaction, using a sophisticated constitutive model for the soil is indispensable. Such a model must adequately describe the material response to a cyclic loading under constant volume (undrained) conditions, amongst others the relaxation of effective stress (pore pressure accumulation) or the effective stress loops repeatedly passed through after a sufficiently large number of cycles (cyclic mobility, stress attractors). The soil behaviour under undrained cyclic loading is manifold, depending on the initial conditions (e.g. density, fabric, effective mean pressure, stress ratio) and the load characteristics (e.g. amplitude of the cycles, application of stress or strain cycles). In order to develop, calibrate and verify a constitutive model with focus to undrained cyclic loading, the data from high-quality laboratory tests comprising a variety of initial conditions and load characteristics are necessary. The purpose of these two companion papers was to provide such database collected for a fine sand. The database consists of numerous undrained cyclic triaxial tests with stress or strain cycles applied to samples consolidated isotropically or anisotropically. Monotonic triaxial tests with drained or undrained conditions have also been performed. Furthermore, drained triaxial, oedometric or isotropic compression tests with several un- and reloading cycles are presented. Part I concentrates on the triaxial tests with monotonic loading or stress cycles. All test data presented herein will be available from the homepage of the first author. As an example of the examination of an existing constitutive model, the experimental data are compared to element test simulations using hypoplasticity with intergranular strain.     

考虑排水条件的软土蠕变特性及模型研究

为了研究珠江三角洲软土在不同的排水条件和应力水平下的蠕变特性,进行了一系列室内三轴排水及不排水蠕变试验。试验结果表明,在排水情况下软土的瞬时变形大,而蠕变变形小,而在不排水情况下蠕变随时间的发展更为显著。在工程中可以提高软土地基的排水性能,来减小蠕变破坏的可能性。引入Singh-Mitchell蠕变方程,建立了珠江三角洲软土排水和不排水的Singh-Mitchell蠕变模型,其中应力-应变关系采用指数函数,应变-时间关系采用幂函数。该模型较好地表达了土体的应力-应变-时间关系。随着排水条件的改变,模型参数的取值有着较大的变化,其原因可归结为,在不同的排水条件下土体的应变-时间关系和应力-应变关系均有很大的不同。     

人工制备结构性土力学特性试验研究

通过对结构性土的研究可以掌握天然土受荷过程中的变形破损过程, 从而为考虑土结构性的结构物设计、地基加固等提供依据。本文尝试了一种简单可行的结构性土样的人工制备方法, 通过在原状土料中加入少量水泥和盐粒以形成颗粒间的胶结作用和大孔隙组构以模拟天然黏土的结构性。然后进行了室内压缩试验和不同胶结强度土样在排水和不排水条件下的三轴剪切试验,分析结构性土样的力学特性,以探讨结构性土的破损过程和变形机理。     

Effects of transport of pore water and material heterogeneity on strain localization of fluid‐saturated gradient‐dependent viscoplastic geomaterial

The purpose of the present paper is to clarify the effects of permeability and initial heterogeneity on the strain localization of fluid‐saturated cohesive soil modelled by a strain gradient‐dependent poro‐viscoplastic constitutive model. The effects of permeability and gradient parameters on the growth rate of the fluctuation were obtained by a linear instability analysis. Deformation behaviour of clay specimens modelled as a viscoplastic model with a second order strain gradient during shear was numerically analysed by a soil–water coupled FEM under both globally undrained and partially drained conditions. It was found that the deformation pattern and the stress–strain curve greatly depend on the permeability, the drainage conditions and the initial non‐homogeneous properties. Copyright © 2005 John Wiley & Sons, Ltd.     

An experimental database for the development,calibration and verification of constitutive models for sand with focus to cyclic loading: part II—tests with strain cycles and combined loading

For numerical studies of geotechnical structures under earthquake loading, aiming to examine a possible failure due to liquefaction, using a sophisticated constitutive model for the soil is indispensable. Such model must adequately describe the material response to a cyclic loading under constant volume (undrained) conditions, amongst others the relaxation of effective stress (pore pressure accumulation) or the effective stress loops repeatedly passed through after a sufficiently large number of cycles (cyclic mobility, stress attractors). The soil behaviour under undrained cyclic loading is manifold, depending on the initial conditions (e.g. density, fabric, effective mean pressure, stress ratio) and the load characteristics (e.g. amplitude of the cycles, application of stress or strain cycles). In order to develop, calibrate and verify a constitutive model with focus to undrained cyclic loading, the data from high-quality laboratory tests comprising a variety of initial conditions and load characteristics are necessary. It is the purpose of these two companion papers to provide such database collected for a fine sand. Part II concentrates on the undrained triaxial tests with strain cycles, where a large range of strain amplitudes has been studied. Furthermore, oedometric and isotropic compression tests as well as drained triaxial tests with un- and reloading cycles are discussed. A combined monotonic and cyclic loading has been also studied in undrained triaxial tests. All test data presented herein will be available from the homepage of the first author. As an example of the examination of an existing constitutive model, the experimental data are compared to element test simulations using hypoplasticity with intergranular strain.     

An anisotropically elastoplastic solution to excavation responses of a circular opening considering three-dimensional strength

Circular opening is commonly encountered in wellbore drilling of petroleum engineering, boring for cast-in situ pile installation, and tunneling excavation. This paper presents a rigorous solution for the elastoplastic responses of the anisotropic soft soil mass around a circular opening excavated under undrained and drained conditions. Both the anisotropic elastoplastic behavior and the 3D strength of the soft clay are incorporated in the present solutions. The well-established anisotropic critical state elastoplastic model S-CLAY1, which can represent the initial fabric anisotropy and stress-induced anisotropy of soft soil, is further modified by the Spatially Mobilized Plane criterion to consider the 3D strength of geomaterials. Then, the investigated problems, excavation of a circular opening under both short-term (undrained) and long-term (drained) conditions, are formulated as a system of first-order differential equations and are solved as initial value problems. The distributions of stress components and anisotropy parameters around the opening, the stress trajectory of a soil particle at the opening wall, as well as the stress–displacement curve at the opening wall are presented to investigate the elastoplastic responses of the opening. Extensive parameters show that the overconsolidated ratio and coefficient of earth pressure at rest (K₀) have remarkable effects on the elastoplastic responses around a circular opening.

     

Stress-strain behaviour of reconstituted illitic clay at different temperatures

Tests on specimens of reconstituted illitic clay have examined the influence of temperature on the mechanical behaviour of clay soils. The program involved consolidation to effective confining pressures up to 1.5 MPa, heating to 100°C, and tests on normally consolidated and overconsolidated specimens with OCR = 2. The tests included isotropic consolidation, undrained triaxial compression with pore water pressure measurement, drained tests along controlled stress paths to investigate yielding behaviour, and undrained tests which involved heating and measurement of the resulting induced pore water pressures. The large strain strength envelope is independent of temperature. However, peak undrained strengths increase with temperature because smaller pore water pressures are generated during shearing. An important contribution from the study is a series of results for the yielding of illitic clay at three different temperatures. For the first time, there is clear evidence of yield loci decreasing in size with increasing temperature. An associated flow rule can be assumed without serious error. The results contribute to the confirmation of a thermal elastic-plastic soil model developed by the authors from cam clay following the addition of a small number of extra assumptions. Depending on the initial stress state, heating under undrained conditions may produce shear failure.     

深圳地区原状花岗岩残积土硬化土模型参数的试验研究

硬化土（hardening soil）模型是岩土工程数值分析中常采用的模型。为研究加载方式和排水条件对原状花岗岩残积土剪切行为的影响,对深圳地区原状花岗岩残积土进行了常规三轴固结排水、常规三轴固结不排水、固结排水侧向卸载和固结不排水侧向卸载4组试验,并由试验结果确定出不同加载方式和排水条件下的硬化土模型参数。结果表明：在常规三轴试验中,试样在剪切过程中表现为先剪缩后剪胀的特性,在侧向卸载试验中,试样始终表现为剪胀的特性;模型参数值与加载路径密切相关,固结不排水侧向卸载试验（CDLU）测得的有效内摩擦角 较常规三轴固结不排水试验（CD）试验大20%,而有效黏聚力 小了46%;CDLU试验得到的三轴压缩试验的参考割线模量 为CD试验的2.9倍,卸载再加载参考割线模量 为CD试验的1.8倍;在相同加载路径下,模型参数值也受排水条件的影响,由CU试验得到的 与CD试验的结果相近,但CU试验得到的 要明显大于CD试验的结果,CU试验得到的 为CD试验的2倍, 为CD试验的3.8倍。因此,在岩土工程数值分析中应根据工程的实际情况确定和选用模型参数值。     

Application of Hysteretic Trends in the Preconsolidation Stress of Unsaturated Soils

This paper involves an evaluation of a relationship describing the evolution in yield stress of unsaturated soils during hydraulic hysteresis, and an application of this relationship in an elasto-plastic framework to predict the compression curves of unsaturated soils under drained (free outflow of air and water with constant suction) or undrained (constant water content with no outflow of water and varying suction) conditions. The yield stress was quantified as the apparent mean effective preconsolidation stress obtained from compression tests reported in the literature on specimens that had experienced different hydraulic paths. It was observed that the preconsolidation stress does not follow a hysteretic path when plotted as a function of matric suction, but does when plotted as a function of the degree of saturation. Accordingly, an existing logarithmic relationship between the preconsolidation stress and matric suction normalized by the air entry suction was found to match the experimental preconsolidation stress results. This same relationship was also able to satisfactorily predict the trends in preconsolidation stress with degree of saturation by substituting the hysteretic soil–water retention curve (SWRC) into the place of the matric suction. The relationship between preconsolidation stress and suction was combined with an elasto-plastic framework to predict the compression curves of soils during drained compression, while the wetting-path relationship between preconsolidation stress and degree of saturation was combined with the framework to predict the compression curves of soils during undrained (constant water content) compression. A good match was obtained with experimental data from the literature, indicating the relevance of considering the hysteretic SWRC and preconsolidation relationships when simulating the behavior of unsaturated soils following different hydro-mechanical paths.     

Non‐isothermal plasticity model for cyclic behaviour of soils

On the one hand, it has been observed that liquefaction‐induced shear deformation of soils accumulates in a cycle‐by‐cycle pattern. On the other hand, it is known that heating could induce plastic hardening. This study deals with the constitutive modelling of the effect that heat may have on the cyclic mechanical properties of cohesive soils, a relatively new area of interest in soil mechanics. In this paper, after a presentation of the thermo‐mechanical framework, a non‐isothermal plasticity cyclic model formulation is presented and discussed. The model calibration is described based on data from laboratory sample tests. It includes numerical simulations of triaxial shear tests at various constant temperatures. Then, the model predictions are compared with experimental results and discussed in the final section. Both drained and undrained loading conditions are considered. The proposed constitutive model shows good ability to capture the characteristic features of behaviour. Copyright © 2007 John Wiley & Sons, Ltd.     

A hypoplastic model for granular soils incorporating anisotropic critical state theory

It is well known that soil is inherently anisotropic and its mechanical behavior is significantly influenced by its fabric anisotropy. Hypoplasticity is increasingly being accepted in the constitutive modeling for soils, in which many salient features, such as nonlinear stress-strain relations, dilatancy, and critical state failure, can be described by a single tensorial equation. However, within the framework of hypoplasticity, modeling fabric anisotropy remains challenging, as the fabric and its evolution are often vaguely assumed without a sound basis. This paper presents a hypoplastic constitutive model for granular soils based on the newly developed anisotropic critical state theory, in which the conditions of fabric anisotropy are concurrently satisfied along with the traditional conditions at the critical state. A deviatoric fabric tensor is introduced into the Gudehus-Bauer hypoplastic model, and a scalar-valued anisotropic state variable signifying the interplay between the fabric and the stress state is used to characterize its impact on the dilatancy and strength of the soils. In addition, fabric evolution during shearing can explicitly be addressed. Modifications have also been undertaken to improve the performance of the undrained response of the model. The anisotropic hypoplastic model can simulate experimental tests for sand under various combinations of principle stress direction, intermediate principal stress (or mode of shearing), soil densities, and confining pressures, and the associated drastic effect of different principal stress orientations in reference to the material axes of anisotropy can be well captured.     

Drained and undrained pullout capacity of a stiff inclusion in a saturated poroelastic matrix

Shear‐lag analysis is used to obtain closed‐form solutions for the problem of a stiff inclusion embedded in a poroelastic soil matrix. The following assumptions are made: the soil matrix and the inclusion are elastic; plane strain conditions apply; and shear stresses at the soil‐inclusion interface follow Coulomb's friction law. Two solutions are obtained, the first one for drained conditions where no excess pore pressures are generated, and the second one for undrained conditions where excess pore pressures are produced and the soil does not change volume during pullout. The solutions are verified by comparing analytical predictions with numerical results obtained using a finite element method. Predictions from the analytical solutions are also compared with results from experiments conducted in a large‐scale pullout box. Both comparisons show good agreement. The analytical solution shows that the pullout capacity in drained and undrained conditions is overall independent of the relative stiffness of the soil and the inclusion. The most important factor controlling the pullout capacity is the coefficient of friction between the soil and the inclusion. Both drained and undrained pullout capacities increase with the coefficient of friction; although the drained capacity shows a proportional increase, it is not so for the undrained capacity. The ratio of undrained to drained pullout capacity is about 0.9 for friction coefficients smaller than 0.2, but can be as small as 0.6 for a coefficient of friction of 1.0. Copyright © 2006 John Wiley & Sons, Ltd.     

Plasticity solutions for soil behaviour around contracting cavities and tunnels

The action of tunnel excavation reduces the in-situ stresses along the excavated circumference and can therefore be simulated by unloading of cavities from the in-situ stress state. Increasing evidence suggests that soil behavior in the plane perpendicular to the tunnel axis can be modelled reasonably by a contracting cylindrical cavity, while movements ahead of an advancing tunnel heading can be better predicted by spherical cavity contraction theory. In the past, solutions for unloading of cavities from in-situ stresses in cohesive-frictional soils have mainly concentrated on the small strain, cylindrical cavity model. Large strain spherical cavity contraction solutions with a non-associated Mohr–Coulomb model do not seem to be widely available for tunnel applications. Also, cavity unloading solutions in undrained clays have been developed only in terms of total stresses with a linear elastic-perfectly plastic soil model. The total stress analyses do not account for the effects of strain hardening/softening, variable soil stiffness, and soil stress history (OCR). The effect of these simplifying assumptions on the predicted soil behavior around tunnels is not known. In this paper, analytical and semi-analytical solutions are presented for unloading of both cylindrical and spherical cavities from in-situ state of stresses under both drained and undrained conditions. The non-associated Mohr-Coulomb model and various critical state theories are used respectively to describe the drained and undrained stress-strain behaviors of the soils. The analytical solutions presented in this paper are developed in terms of large strain formulations. These solutions can be used to serve two main purposes: (1) to provide models for predicting soil behavior around tunnels; (2) to provide valuable benchmark solutions for verifying various numerical methods involving both Mohr–Coulomb and critical state plasticity models. Copyright © 1999 John Wiley & Sons, Ltd.