一种适用于饱和黏土循环动力分析边界面塑性模型的隐式积分算法

修正了传统隐式回映算法,建立了适用于饱和黏土循环动力分析的边界面塑性模型的完全隐式积分格式。该模型基于无弹性域概念和临界状态理论,采用各向同性、运动硬化准则、旋转的边界面,并引入表征土体结构损伤和重塑程度的损伤变量以反映循环载荷作用下饱和黏土的各向异性、刚度、强度软化及塑性变形累积等特征。针对等压固结 和偏压固结 的饱和高岭黏土的不排水三轴试验进行模拟,采用不同的应变增量步长进行计算,并与试验数据对比,结果表明,修正隐式回映算法应用于该类边界面模型的合理性、积分格式的精确性和稳定性;另外,结合有限元软件自动时间步长的增量迭代解法,对饱和黏土应力控制的不排水动三轴试验进行预测,结果表明,修正的适用于该边界面的塑性模型隐式回映算法可以得到比较合理的数值分析结果,能够反映饱和黏土的循环刚度的退化和强度的弱化等动力特性。     

Implicit and explicit integration schemes in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay

Two integration algorithms, namely the implicit return mapping and explicit sub-stepping schemes, are adopted in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay and are implemented into finite element code. The model is a representative of a series of bounding surface models that have typical characteristics, including isotropic and kinematic hardening rules and a rotational bounding surface to capture complex but important cyclic behaviours of soils, such as cyclic shakedown and degradation. However, there is no explicit current yield surface in the model to which the conventional implicit algorithm returns the stress state back or the sub-stepping integration corrects the drift of the stress state. Hence, necessary modifications have been made for both of the integration schemes. First, the image stress point is mapped or corrected to the bounding surface instead of mapping back or correcting the stress state to the yield surface. Second, the unloading–loading criterion is checked to determine the image stress point rather than checking the yield criterion after giving the trial stress state in a conventional way. Comparative studies on the accuracy, stability and efficiency of the two integration schemes are conducted not only at the element level but also in solving boundary value problems of monotonic and cyclic bearing behaviours of rigid footings on saturated clay. For smaller strain increments, there is no significant difference in the accuracy between the two integration schemes, but the explicit integration shows a higher efficiency and accuracy. For relatively larger increments, the implicit return mapping algorithm presents good accuracy and more robustness, while the sub-stepping algorithm shows deteriorating accuracy and suffers the convergence problem. With the tolerance used in the present model, the bearing capacity of the rigid footing predicted by the return mapping algorithm is closer to the available analytical and numerical solutions, while the bearing capacity predicted by the sub-stepping algorithm shows a marginal increase.     

Implicit integration under mixed controls of a breakage model for unsaturated crushable soils

This paper discusses a series of stress point algorithms for a breakage model for unsaturated granular soils. Such model is characterized by highly nonlinear coupling terms introduced by breakage‐dependent hydro‐mechanical energy potentials. To integrate accurately and efficiently its constitutive equations, specific algorithms have been formulated using a backward Euler scheme. In particular, because implementation and verification of unsaturated soil models often require the use of mixed controls, the incorporation of various hydro‐mechanical conditions has been tackled. First, it is shown that the degree of saturation can be replaced with suction in the constitutive equations through a partial Legendre transformation of the energy potentials, thus changing the thermomechanical state variables and enabling a straightforward implementation of a different control mode. Then, to accommodate more complex control scenarios without redefining the energy potentials, a hybrid strategy has been used, combining the return mapping scheme with linearized constraints. It is shown that this linearization strategy guarantees similar levels of accuracy compared with a conventional strain–suction‐controlled implicit integration. In addition, it is shown that the use of linearized constraints offers the possibility to use the same framework to integrate a variety of control conditions (e.g., net stress and/or water‐content control). The convergence profiles indicate that both schemes preserve the advantages of implicit integration, that is, asymptotic quadratic convergence and unconditional stability. Finally, the performance of the two implicit schemes has been compared with that of an explicit algorithm with automatic sub‐stepping and error control, showing that for the selected breakage model, implicit integration leads to a significant reduction of the computational cost. Such features support the use of the proposed hybrid scheme also in other modeling contexts, especially when strongly nonlinear models have to be implemented and/or validated by using non‐standard hydro‐mechanical control conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of liquefaction behaviour of granular materials using Discrete Element Method

In this paper, numerical simulation of 3-dimensional assemblies of 1000 polydisperse sphere particles using Discrete Element Method (DEM) is used to study the liquefaction behaviour of granular materials. Numerical simulations of cyclic triaxial shear tests under undrained conditions are performed at different confining pressures under constant strain amplitude. Results obtained in these numerical simulations indicate that with increase in confining pressure there is an increase in liquefaction resistance.     

Investigation of monotonic and cyclic behavior of sand using a bounding surface plasticity model

Developing the pore water pressures in loose to medium sands below the water table may lead to liquefaction during earthquakes. The simulation of liquefaction (cyclic mobility and flow liquefaction) in sandy soils is one of the major challenges in constitutive modeling of soils. This paper presents the simulation of sand behavior using a critical state bounding surface plasticity model (Dafalias and Manzari’s model, 2004) during monotonic and cyclic loading. The drained, undrained, and cyclic triaxial tests were simulated using Dafalias and Manzari’s model. The simulation results showed that the model predicts behavior of sand, reasonably well. Also, for CSR?<?0.2, number of cycles for liquefaction is significantly increased. The residual strength of Babolsar sand is produced when it is deformed to an axial strain of 20 to 25%.     

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.     

Laboratory Investigation of the Resistance of Tailings and Natural Sediments to Cyclic Loading

A number of cyclic triaxial tests were carried out on mine tailings and natural sediment samples under undrained conditions to investigate their resistance to cyclic loading. The tests were performed on more than 100 samples with a cyclic shear stress ratio ranging from 0.10 to 0.40 under varying void ratio and the same confining pressure. It was observed that the axial strain and excess pore water pressure increased with the number of loading cycles while the effective stress decreased with increasing number of loading cycles. The liquefaction resistance of the tailings was also observed to be higher than that of natural soils with similar particle size distribution, void ratio and plasticity index. It was observed that the influence of specific gravity on the cyclic strength of mine tailings is significant. The results showed that the cyclic resistance of the tailings was not strongly influenced by plasticity index for low plasticity tailings. A boundary relationship between void ratio and normalized cyclic resistance ratio was established based on the results.     

纤维加筋砂土抗液化试验与数值模拟

砂土液化是导致重大地震灾害的主要原因之一。本研究探讨了天然纤维加筋砂土在循环荷载作用下的抗液化性能。在不排水条件下,对具有不同纤维含量的加筋砂土试样进行了一系列循环三轴试验,研究了饱和砂土的液化特性以及循环剪应变幅值、纤维含量对饱和砂土抗液化性能的影响。此外,通过模拟已完成的循环三轴试验,建立了二维有限元数值模型,并对具有不同纤维含量的加筋砂土进行了参数标定。研究结果表明:(1)增加循环剪应变幅值将促进超孔压累积,使得滞回曲线斜率和平均有效应力减小速度加快;(2)纤维的存在能够减缓超孔压的累积,随着纤维含量增加,加筋砂土抗液化能力得到明显提高;(3)标定后的本构模型参数能可靠地用于模拟纤维加筋砂土的液化响应。研究结果为饱和砂土抗液化问题与纤维加筋砂土的数值模拟提供了有价值的参考。     

Microscale Energy Dissipation Mechanisms in Cyclically-Loaded Granular Soils

This paper utilizes the Discrete Element Method to characterize energy dissipation mechanisms in cyclically loaded soils based on micromechanical considerations. Computational simulations of consolidated undrained cyclic triaxial tests were conducted at various relative densities and were subjected to cyclic loading of different frequencies and shear strain amplitudes. The different components of microscale energies were monitored during the course of the simulations and characterized into input and dissipated energies. A comparison is made between the dissipated energy computed from microscopic energy components and macroscopic energy calculated based on the area of the deviator stress-axial strain loops. These energies are then used to obtain the specific damping capacity defined as the ratio of dissipated energy during one cycle to the maximum stored elastic energy during the same cycle. The conducted simulations highlight the importance of calculating actual stored energy in the system as opposed to approximating it to be that calculated as the triangular area under the secant modulus. Finally, a series of simulations that resulted in liquefaction are discussed, and the amount of energy dissipated to liquefaction is examined based on these results.     

Implicit numerical integration for a kinematic hardening soil plasticity model

Soil models based on kinematic hardening together with elements of bounding surface plasticity, provide a means of introducing some memory of recent history and stiffness variation in the predicted response of soils. Such models provide an improvement on simple elasto‐plastic models in describing soil behaviour under non‐monotonic loading. Routine use of such models requires robust numerical integration schemes. Explicit integration of highly non‐linear models requires extremely small steps in order to guarantee convergence. Here, a fully implicit scheme is presented for a simple kinematic hardening extension of the Cam clay soil model. The algorithm is based on the operator split methodology and the implicit Euler backward integration scheme is proposed to integrate the rate form of the constitutive relations. This algorithm maintains a quadratic rate of asymptotic convergence when used with a Newton–Raphson iterative procedure. Various strain‐driven axisymmetric triaxial paths are simulated in order to demonstrate the efficiency and good performance of the proposed algorithm. Copyright © 2001 John Wiley & Sons, Ltd.     

Shortest link method for contact detection in discrete element method

With the increasing demand for discrete element simulations with larger number of particles and more realistic particle geometries, the need for efficient contact detection algorithms is more evident. To date, the class of common plane (CP) methods is among the most effective and widely used contact detection algorithms in discrete element simulations of polygonal and polyhedral particles. This paper introduces a new approach to obtain the CP by employing a newly introduced concept of ‘shortest link’. Among all the possible line segments that connect any point on the surface of particle A to any point on the surface of particle B, the one with the shortest length defines the shortest link between the two particles. The perpendicular bisector plane of the shortest link fulfils all the conditions of a CP, suggesting that CP can be obtained by seeking the shortest link. A new algorithm, called shortest link method (SLM), is proposed to obtain the shortest link and subsequently the CP between any two polyhedral particles. Comparison of the analysis time between SLM and previously introduced algorithms demonstrate that SLM results in a substantial speed up for polyhedral particles contact detection. Copyright © 2006 John Wiley & Sons, Ltd.     

粗粒土边界面塑性模型及其积分算法

粗粒土的强度、变形特性对土石坝、边坡和路基等工程的安全性与稳定性有着至关重要的影响。针对粗粒土在复杂应力状态下的强度和变形特性,在边界面塑性理论和临界状态理论框架下,通过引入状态参数和动态临界状态线建立了粗粒土状态相关边界面塑性模型。模型不仅能够模拟粗粒土的应变硬化和体积收缩行为,还能描述应变软化和体积膨胀特性。基于ABAQUS的二次开发平台,结合带误差控制的改进欧拉积分算法编写了边界面塑性模型的UMAT子程序。通过设置不同的应变增量步和误差容许值,对改进欧拉积分算法的精确性和收敛性进行了分析。对不同密实状态和围压下粗粒土三轴排水剪切试验进行了模拟,验证了带误差控制的改进欧拉积分算法应用于粗粒土边界面塑性模型的合理性,为进一步工程应用奠定了基础。     

A unified plasticity model for large post-liquefaction shear deformation of sand

Based on previous experimental findings and theoretical developments, this paper presents the formulation and numerical algorithms of a novel constitutive model for sand with special considerations for cyclic behaviour and accumulation of large post-liquefaction shear deformation. Appropriate formulation for three volumetric strain components enables the model to accurately predict loading and load reversal behaviour of sand, fully capturing the features of cyclic mobility. Compliance with the volumetric compatibility condition, along with reversible and irreversible dilatancy, allows for physically based simulation of the generation and accumulation of shear strain at zero effective stress after initial liquefaction. A state parameter was incorporated for compatibility with critical state soil mechanics, enabling the unified simulation of sand at various densities and confining pressures with a same set of parameters. The determination methods for the 14 model parameters are outlined in the paper. The model was implemented into the open source finite-element framework OpenSees using a cutting-plane stress integration scheme with substepping. The potentials of the model and its numerical implementation were explored via simulations of classical drained and undrained triaxial experiments, undrained cyclic torsional experiments, and a dynamic centrifuge experiment on a single pile in liquefiable soil. The results showed the model’s great capabilities in simulating small to large deformation in the pre- to post-liquefaction regime of sand.     

随机地震荷载作用下饱和粉土的液化特性

为了考虑随机地震荷载作用下饱和粉土的液化特性，研究了如何通过合理的地震响应分析确定饱和粉土层受到的随机地震荷载，进一步探讨了随机地震荷载作用下的动三轴试验技术，以及利用此种试验分析饱和粉土液化特性的方法。研究结果为形成判别饱和粉土层地震液化可能性的新方法提供了依据。     

Mechanics of origin of flow liquefaction instability under proportional strain triaxial compression

We define a flow liquefaction potential for determining flow liquefaction susceptibility during proportional strain triaxial compression. The flow liquefaction potential is a function of inconsistency between the natural dilative tendency of the soil and the imposed dilatancy during proportional strain triaxial compression. It helps us analyze why given the right conditions, a loose soil that contracts during drained triaxial compression and liquefies under undrained triaxial compression may be stable under proportional strain triaxial compression. Conversely, we also use the flow liquefaction potential to analyze why a dense soil that dilates during drained triaxial compression and is stable under undrained triaxial compression may liquefy under proportional strain triaxial compression. The undrained loose case is a special case of proportional strain triaxial compression under which a soil can liquefy. The central objective of this paper was to investigate the origins of flow liquefaction instability. Hence, we also analyze stress evolution during proportional strain triaxial compression and discuss the mechanics of the test leading up to flow liquefaction instability. We arrive at a necessary precursor for instability, which can serve as a warning sign for flow liquefaction instability, while the soil is still stable. The precursor is not a condition of sufficiency and should also not be confused with the onset of instability itself. The same loading must be applied continuously to induce flow liquefaction instability. The current progress is encouraging and facilitates a deeper understanding of origin of flow liquefaction instabilities.     

Improved numerical algorithms for frictional contact in pile penetration analysis

This paper presents a numerical formulation for frictional contact problems associated with pile penetration. The frictional contact at the soil–pile interface is formulated using the theory of hardening/softening plasticity, so that advanced models for the interface can be dealt with. A smooth discretisation of the pile surface is proposed using Bézier polynomials. An automatic load stepping scheme is proposed, which features an error control algorithm and automatic subincrementation of the load increments. The numerical algorithms are then used to analyse the installation process of pushed-in axial piles. It is shown that the smooth discretisation of the pile surface is effective in reducing the oscillation in the predicted pile resistances and the automatic load stepping scheme outperforms the classical Newton–Raphson scheme for this type of problem.     

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.     

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.     

Numerical modeling of changes in anisotropy during liquefaction using a generalized constitutive model

A constitutive model with rotation hardening was generalized from the triaxial compression state to a general stress state. With the generalized model, numerical simulations of cyclic and monotonic undrained triaxial tests were conducted to reproduce the phenomenon of continuous, orderly and rapid changes in anisotropy during liquefaction. The simulated results demonstrated that when sand enters the liquefaction stage, the yield surface in the stress space rotates quickly, causing continuous and rapid changes in anisotropy. Through comparison of the simulated and experimental results, the generalized constitutive model was shown to be able to capture the fundamental behaviors of sand demonstrated by the experimental data, and the rotational hardening rule adopted in the generalized model was proven suitable for describing the continuous, orderly and rapid changes in anisotropy that occur during liquefaction.