非共轴本构模型的数值计算问题

应用角点非共轴本构模型结合有限元方法分析浅基础的受力沉降问题。当非共轴塑性模量较小时,包含非共轴模型子程序的ABAQUS的非线性计算可能不收敛。深入地分析表明,不收敛问题是由于在塑性变形刚开始时,非共轴的塑性应变增量明显大于共轴的塑性应变增量。为了克服数值不收敛问题,对原有的角点非共轴模型加以改进。在改进的模型中,非共轴塑性模量是累加塑性剪应变的函数。计算结果表明,这个改进的非共轴模型提高了非线性计算的收敛性,同时保持了原有非共轴模型的功能。     

非共轴本构模型在地基承载力数值计算中若干影响因素的探讨

土体在剪切变形过程中产生主应力方向的旋转时,主应变增量方向与主应力方向之间存在着非共轴现象,然而传统的弹塑性本构模型未能考虑该现象的影响。通过在屈服面的切线方向增加一项非共轴塑性应变增量,即可实现对非共轴现象的反映。采用显式积分算法和自动分步方法,将非共轴本构模型运用到桶形基础地基承载力问题的有限元计算中,并讨论了流动法则、内摩擦角、膨胀角等因素与非共轴模型的联系。计算结果表明：采用有限元程序默认容许误差时,该本构模型可达到理想的收敛精度,并且,该模型对关联、非关联流动法则均适用。采用共轴模型进行数值计算时,不同流动法则对计算结果的影响可以忽略;采用非共轴模型时,不同流动法则的计算结果之间存在差异。非共轴现象对地基承载力-位移曲线具有软化作用,并且,该软化作用在采用非关联流动法则时变得更加明显     

基于广义位势理论的土的非共轴特性研究

传统的塑性位势理论隐含了应力主方向和塑性应变增量主方向共轴的假定,无法客观地描述主应力轴旋转过程中的非共轴现象。基于广义位势理论提出的拟弹性弹塑性本构模型,把总的塑性应变分解为满足弹性分解准则的拟弹性部分和符合传统塑性理论假设的纯塑性部分,分解后建立的模型更为合理和简便,同时又可以解决土的非共轴问题。通过单剪试验结果的验证表明,基于广义位势理论的拟弹性弹塑性模型的模拟效果较好,传统的弹塑性模型（共轴模型）模拟得到的主应力方向和塑性主应变增量方向保持共轴,而拟弹性弹塑性模型（非共轴模型）的模拟结果则能够合理地描述主应力轴旋转过程中的非共轴特性,结果更符合实际,从而为解决土的非共轴特性问题提供了一种有效的方法。     

Numerical aspects of non-coaxial model implementations

This paper ascertains the reasons for the numerical problems when the yield vertex non-coaxial model is implemented in the finite element analysis to predict the behaviour of complicated geotechnical engineering problems. The numerical problem, reflected in the failure of convergence in the non-linear solutions in ABAQUS, is likely to happen when a smaller non-coaxial plastic modulus is used. It is found that a large non-coaxial influence at the start of elastoplastic loading, which causes the predominance of non-coaxial plastic strain rate over the coaxial plastic strain rate in stress–strain responses, is the reason for the numerical problem. The original yield vertex non-coaxial model is modified to overcome the numerical problem. Instead of a constant non-coaxial plastic modulus in the original yield vertex non-coaxial model, the non-coaxial modulus is made a function of cumulative deviatoric plastic strain. It shows that the modified non-coaxial model retains the functions of the original non-coaxial model. Meanwhile, it overcomes the non-convergence problem when a smaller non-coaxial plastic modulus, representing a larger non-coaxial influence, is used.     

基于临界状态模型的砂土非共轴本构模拟

传统的砂土本构理论隐含了应力和塑性应变率的共轴条件,无法客观描述主应力轴旋转试验中的非共轴现象,并且当密度和围压变化较大时也不适用。基于材料状态相关砂土临界状态概念,将Pietruszczak和Stolle所提出的砂土本构模型进行了改进,并在模型中引入非共轴塑性流动理论来描述非共轴现象。通过对单剪试验和空心圆柱试验进行数值模拟,表明基于临界状态理论的非共轴模型能够合理描述主应力轴旋转过程中砂土的非共轴变形特性     

Elasto-viscoplastic finite element study of the effect of degradation on bearing capacity of footing on clay ground

The bearing capacity of footing has been studied by both conventional and numerical methods by many researchers. However, degradation of the microstructure of material, that is, a change in the microstructure of the soil, has not been adequately taken into account. Degradation of microstructure causes strain softening of materials and it leads to strain localization such as shear bands and slip bands. From an engineering point of view the strain localization is crucial because it is a precursor of failure. In the present study, finite element analyses of the bearing capacity of a shallow foundation on homogeneous and inhomogeneous saturated clay strata have been conducted using an elasto-viscoplastic soil constitutive model of microstructure change. A series of analyses of footing on clay deposit with different microstructure parameters have been carried out. Numerical results show that strain localization can be predicted during the loading of rigid footing on highly structured soil and strain localization affects the footing–soil interaction. The effects of footing roughness on the failure mechanism are also discussed in the study.     

Developing artificial neural network models to predict allowable bearing capacity and elastic settlement of shallow foundation in Sharjah,United Arab Emirates

This research proposes the use of artificial neural network to predict the allowable bearing capacity and elastic settlement of shallow foundation on granular soils in Sharjah, United Arab Emirates. Data obtained from existing soil reports of 600 boreholes were used to train and validate the model. Three parameters (footing width, effective unit weight, and SPT blow count) are considered to have the most significant impact on the magnitude of allowable bearing capacity and elastic settlement of shallow foundations, and thus were used as the model inputs. Throughout the study, depth of footing was limited to 1.5 m below existing ground level and water table depth taken at the level of the footing. Performance comparison of the developed models (in terms of coefficient of determination, root mean square error, and mean absolute error) revealed that the developed artificial neural network models could be effectively used for predicting the allowable bearing capacity and elastic settlement. As such, the developed models can be used at the preliminary stage of estimating the allowable bearing capacity and settlements of shallow foundations on granular soils, instead of the conventional methods.     

Implicit and explicit procedures for the yield vertex non-coaxial theory

The yield vertex non-coaxial model is different from classical elastoplastic models, in that there is an additional plastic strain rate tangential to yield surfaces, as well as the plastic strain rate normal to yield surfaces, when orientations of principal stress change. This feature raises concerns on its finite element implementations. In nonlinear finite element numerical iterations, a large tangential plastic strain rate is likely to make the trial total strain rate direct inside a yield surface, which entails convergence difficulty. Some modifications are introduced on the non-coaxial model itself to make numerical convergence easier in the work published in Yang and Yu (2010) [20]. This paper is an extension of the previous work. Instead of modifying the non-coaxial model itself, this paper concerns the use of finite element explicit procedure, which is suitable for highly discontinuous problems. The simulations of shallow foundation load-settlement responses indicate that the finite element explicit procedure, assisted with a robust and explicit automatic substepping integration scheme of the non-coaxial model, does not encounter numerical difficulty. In addition, the overall trends of implicit and explicit simulations are similar.     

Performance of square footing with structural skirt resting on sand

The paper presents the performance of a square footing with a structural skirt resting on sand and subjected to a vertical load through an experimental study. A series of tests were conducted in a model test tank to evaluate the performance in terms of improvement in bearing capacity and reduction in settlement of a square footing with and without a structural skirt. The results of the study reveal that this type of reinforcement increases the bearing capacity, reduces the settlement and modifies the load settlement behaviour of the footing. The various factors influencing the bearing capacity improvements and settlement reduction using a structural skirt are identified. Skirt factors are proposed which can be introduced into the general ultimate bearing capacity equation for a square footing resting on sand. The predictions made through the modified equation are in reasonable agreement with the experimental results. The bearing capacity of square footing is increased in the range of 11.2 to 70%. The improvement in bearing capacity decreases with the increase in base roughness of the footing. Further, an equation for a settlement reduction factor is proposed which can be used to calculate the settlement of the square footing with structural skirt resting on sand. The settlement reduction factor (SRF, defined as the ratio of settlement of footing with structural skirt to the settlement of footing without structural skirt at a given load) was in the range 0.11 to 1.0 depending on applied load and skirt depth ratio with the use of a structural skirt. The results further reveal that for a given depth of the skirt, the settlement reduction factor decreases with the increase in applied load. The improvement in the bearing capacity and reduction in settlement of a square footing with a structural skirt resting on sand are dependent on the geometrical and structural properties of the skirt, footing, sand characteristics and interface conditions of the sand–skirt–footing system.     

短桩基础桩-土共同作用的原位测试与数值分析

通过对湖沼平原一个6层建筑的短桩原位测试和数值分析,研究了短桩基础的荷载传递与变形特性。试验结果表明,桩基的实际承载力大于计算值,而沉降远小于计算值。桩顶反力和基底土压力的观测结果显示：地基土所分担的荷载较小。对短桩基础的桩-土共同作用进行数值分析,得到了地基土层内的附加应力和分层沉降变化情况。研究结果表明,湖沼平原短桩基础的沉降以桩端以下土层的压缩变形为主。浅部硬土层既作为短桩的持力层,又有利于附加应力的扩散,从而减小了沉降。     

Load‐displacement and bearing capacity of foundations on granular soils using a multi‐surface kinematic constitutive soil model

A finite element approach based on an advanced multi‐surface kinematic constitutive model is used to evaluate the bearing capacity of footings resting on granular soils. Unlike simple elastic‐perfectly plastic models, often applied to granular foundation problems, the present model realistically accounts for stress dependency of the friction angle, strain softening–hardening and non‐associativity. After the model and its implementation into a finite element code are briefly discussed, the numerical difficulty due to the singularity at the footing edge is addressed. The bearing capacity factor N_γ is then calculated for different granular materials. The effect of footing size, shape, relative density and roughness on the ultimate bearing capacity are studied and the computed results compare very favourably with the general experimental trends. In addition, it is shown that the finite element solution can clearly represent counteracting mechanisms of progressive failure which have an important effect on the bearing capacity of granular foundations. Copyright © 2006 John Wiley & Sons, Ltd.     

从各向异性的角度解释和模拟土的非共轴特性

从发挥面的角度出发,分析论证各向异性是引起岩土材料出现非共轴现象的根本原因,得到与材料力学一致的结论。当共轭的两发挥面与沉积面的夹角不相等时,主应力面上将出现塑性应变增量的切向分量,所以塑性应变增量的主方向与应力的主方向非共轴。按照这一结论,对非共轴的数值模拟,也应当根据各向异性本构模型进行。为考虑各向异性影响新近提出的各向异性变换应力法,改变了各应力分量的相对大小,得到的各向异性变换应力张量与真实应力张量的主方向不一致,因此也能反映非共轴。利用各向异性变换应力法,能够在现有的弹塑性本构模型的框架下,描述土的非共轴现象。以各向异性UH模型为例,预测各种加载条件下的非共轴变形,验证了该方法的有效性。     

Stability of eccentrically loaded footings on slopes

The stability of eccentrically loaded strip footings on slopes was investigated using the method of finite element analysis based on the theory of elasto-plasticity. The analysis was done for two different soils involving three levels of slope angle, six footing locations, and two levels of load eccentricity plus central vertical loading. The strip footing analysed was a 3-ft (0.9 m) wide reinforced concrete footing embedded to a depth of 3 ft (0.9 m). The analysis focused on footing settlement, plastic yielding of soil, and ultimate bearing capacity. The results of analysis show that the influence of load eccentricity on footing pressure vs. footing centre settlement is negligibly small. However, the progressive soil yielding and ultimate bearing capacity are greatly affected by load eccentricity. Furthermore, the effect of load eccentricity differs considerably with the load location relative to the footing centre and slope crest. The ultimate bearing capacity for the eccentric load located on the slope side is significantly greater than that for the load located on the other side of the footing centre. For a 2(H): 1(V) slope in silty clay, the effect of slope on footing stability decreases with increasing footing location from slope crest as would be expected, and diminishes when the footing is located from the crest at about 5-times the footing width.     

A macroelement formulation for shallow foundations on cohesive and frictional soils

The scope of this paper is to present a macroelement model for shallow foundations encompassing the majority of combinations of soil and foundation–soil interface conditions that are interesting for practical applications. The basic idea of the formulation is to raise the common assumption that the surface of ultimate loads of the foundation is identified as a yield surface in the space of force parameters which the footing is subjected to. Instead, each non‐linear mechanism participating in the global response of the system is modelled independently and the surface of ultimate loads is retrieved as the combined result of all active mechanisms. This allows formulating each mechanism by respecting its particular characteristics and offers the possibility of activating, modifying or deactivating each mechanism according to the context of application. The model comprises three non‐linear mechanisms: (a) the mechanism of sliding at the soil–footing interface, (b) the mechanism of soil yielding in the vicinity of the footing and (c) the mechanism of uplift as the footing may get detached from the soil. The first two are irreversible and dissipative and are combined within a multi‐mechanism plasticity formulation. The third mechanism is reversible and non‐dissipative. It is reproduced with a phenomenological non‐linear hyperelastic model. The model is validated with respect to the existing results for shallow foundations under quasi‐static loading tests. It is shown that although the ultimate surface of the foundation is not explicitly used in the formulation of the model, the obtained force states by the model are always contained within it. Copyright © 2010 John Wiley & Sons, Ltd.     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

Least‐square support vector machine applied to settlement of shallow foundations on cohesionless soils

This paper examines the potential of least‐square support vector machine (LSVVM) in the prediction of settlement of shallow foundation on cohesionless soil. In LSSVM, Vapnik's ε‐insensitive loss function has been replaced by a cost function that corresponds to a form of ridge regression. The LSSVM involves equality instead of inequality constraints and works with a least‐squares cost function. The five input variables used for the LSSVM for the prediction of settlement are footing width (B), footing length (L), footing net applied pressure (P), average standard penetration test value (N) and footing embedment depth (d). Comparison between LSSVM and some of the traditional interpretation methods are also presented. LSSVM has been used to compute error bar. The results presented in this paper clearly highlight that the LSSVM is a robust tool for prediction of settlement of shallow foundation on cohesionless soil. Copyright © 2008 John Wiley & Sons, Ltd.     

横-竖立体加筋地基中附加应力的分析计算

建立了横-竖立体加筋（H-V筋）地基的有限元模型,通过分析地基中的竖向应力分布、水平向位移分布以及筋-土界面相互作用,发现横-竖立体加筋地基中的竖向应力在筋材下方出现扩散和重分布,并逐渐向土体下部传递,使得土体中整体的应力分布更加均匀;同时,横-竖筋材中的竖筋类似于一个侧壁,其提供的垂直侧向力约束了介于竖筋间的土体,限制了土体的侧向水平位移,使得地基中筋材上部土体的侧向水平位移变小。基于有限元模拟对横-竖立体加筋地基加固机制的认识,将横-竖立体筋视为作用在地基上的一维弹性地基梁,通过弹性地基梁理论,根据弗拉曼解推导求解了横-竖立体加筋地基中任意一点竖向附加应力的计算表达式。将模型计算结果与有限元模拟所得结果进行对比发现两者吻合良好。     

A semi‐empirical method for the liquefaction analysis of offshore foundations

This paper presents a relatively simple method for three‐dimensional liquefaction analysis of granular soil under offshore foundations. In this method, the Mohr–Coulomb model, which defines the elasto–plastic stress–strain relationship under monotonic loading, is modified to accommodate the plastic strains generated by cyclic loading. The effects of cyclic loading, evaluated from the results of laboratory tests on saturated samples of soil, are incorporated into the model. The method is implemented in an efficient finite element program for analyses of three‐dimensional consolidating soil. The practicability of the model is demonstrated by analysis of a typical offshore foundation, and the predictions of the numerical analysis are compared with the observed behaviour of the foundation. Copyright © 2000 John Wiley & Sons, Ltd.     

Constitutive modelling of granular materials using a contact normal-based fabric tensor

This paper presents a fabric tensor-based bounding surface model accounting for anisotropic behaviour (e.g. the dependency of peak strength on loading direction and non-coaxial deformation) of granular materials. This model is developed based on a well-calibrated isotropic bounding surface model. The yield surface is modified by incorporating the back stress which is proportional to a contact normal-based fabric tensor for characterising fabric anisotropy. The evolution law of the fabric tensor, which is dependent on both rates of the stress ratio and the plastic strain, rules that the material fabric tends to align with the loading direction and evolves towards a unique critical state fabric tensor under monotonic shearing. The incorporation of the evolution law leads to a rotational hardening of the yield surface. The anisotropic critical state is assumed to be independent of the initial values of void ratio and fabric tensor. The critical state fabric tensor has the same intermediate stress ratio (i.e. b value) and principal directions as the critical state stress tensor. A non-associated flow rule in the deviatoric plane is adopted, which is able to predict the non-coaxial flow naturally. The stress–strain relation and fabric evolution of model predictions show a satisfactory agreement with DEM simulation results under monotonic shearing with different loading directions. The model is also validated by comparing with laboratory test results of Leighton Buzzard sand and Toyoura sand under various loading paths. The comparison results demonstrate encouraging applicability of the model for predicting the anisotropic behaviour of granular materials.

     

Effects of principal stress rotation on wave-induced soil response in a poro-elastoplastic sandy seabed

In this study, a constitutive model is developed in order to investigate wave–seabed interactions. This model takes into account the impact of principal stress rotation (PSR) and is based on the generalized plasticity theory, in which plastic strain generated by PSR is considered an additional item in the constitutive relationship of soil. The normalized loading direction and plastic flow direction were determined based on the stress tensor invariant. Comparisons between the present model and previous Hollow Cylinder Apparatus tests and geotechnical centrifugal wave tests show good agreement. Numerical results show the effects of PSR on predictions of liquefaction potential due to: (a) the cumulative impact of plastic strain in the seafloor and (b) the buildup of pore pressure. Parametric study shows that the model parameters, including the wave and seabed parameters, have significant effects on the wave-induced soil liquefaction.