波浪作用下饱和砂质海床土体与管线相互作用规律研究

针对波浪引起的饱和砂质海床土体和管线相互作用问题,将Biot动力固结理论与笔者课题组提出的砂土液化变形弹塑性本构模型相耦合,较为合理地再现了简谐波浪作用下较浅饱和砂质海床中管线周围可液化海床土体的超静孔隙水压力瞬态累积变化规律与液化过程。数值计算结果与Sumer等的试验规律一致。结果表明：由于管线的存在,改变了饱和砂质海床液化区域的空间分布。液化首先由管线下部土体开始产生,随着波浪荷载的持续作用,液化区域沿着管线外壁向上演化;同时海床表层土体产生液化并向深层发展,最终管线周围土体都发生液化,这是导致空管上浮的主要原因。当饱和砂质海床中存在管线时,管线附近海床土体液化深度明显变深。超静孔压累积和渗透力变化的耦合作用是导致饱和砂质海床土体产生液化的原因。与将海床土体视为饱和弹性多孔介质相比,可考虑液化全过程的弹塑性动力分析能更为合理地揭示实际波浪作用下饱和砂质海床土体的渗流场和应力场的瞬态时空演变规律。     

非线性波浪荷载作用下海上风机管桩基础周围海床液化研究

为研究近浅海领域内海上风机大直径单桩基础周围的砂质海床在非线性荷载作用下的瞬时液化的稳定性,在OpenFOAM平台上建立波浪-单桩-海床三维数值模型(WSSI模型),采用olaFoam求解器求解RANS方程模拟波浪的非线性运动,将Biot方程作为海床模型的控制方程研究单桩周围海床的液化规律。对以往的数值解析结果和实验结果进行对比分析,验证本文建立三维波浪模型和海床模型的准确性和有效性。将建立的波浪-管桩-海床相互作用三维数值模型应用到实际工况中,对数值模型在实际工程中的海床响应进行分析评估。研究结果表明,单桩周围海床的孔隙水压力场和有效应力场在水深方向发生较大变化,在水平方向变化较小。在该工况下海床的最大液化深度可达到10 m,且在垂直于波浪运动方向的海床更容易发生液化,为实际工程预测海床液化提供一定的技术支持。     

近海欠密实砂质海床内波致渐进液化特征研究

利用一个经过广泛验证的数值模型FSSI-CAS 2D为计算工具,采用砂土的高级本构模型Pastor-Zienkiewicz-Mark III (PZIII) 描述海床砂土的动态力学行为,定量研究松散海床地基土在波浪作用下,其内部的液化过程和特征,以加深对波致海床液化特征、性质的认识。计算结果分析表明,开发的耦合数值模型FSSI-CAS 2D能够很好地捕捉到波浪作用下欠密实海床的动力响应特征,以及海床内的累积液化过程等一些列的非线性物理现象。研究表明,波浪导致的松砂海床液化是一个渐进过程,海床表面首先液化,并逐渐向下扩展。     

Breaking wave‐induced response and instability of seabed around caisson breakwater

Breaking‐wave‐induced dynamic response and instability of seabed around a caisson breakwater are investigated. A seabed‐rubble‐breakwater system is modeled using finite elements. The impact response of the porous seabed and rubble foundation is assumed to be governed by the coupled Biot equations, and three possible formulations are considered with respect to the inclusion of inertial terms. The response is presented in terms of shear stress and pore pressure distributions at three locations underneath the breakwater. The effect of seabed and wave parameters and the inertial terms on the impact response is investigated through parametric studies. Analyses show that usually partly dynamic formulation yields the largest response amplitudes as compared to the fully dynamic formulation, which is the most complete form. The instability of seabed and rubble mound as a result of instantaneous liquefaction is also studied. Breaking wave‐induced pressures in some cases are found to cause liquefaction in the rubble and the seabed. The effect of some parameters on the instability is found to be significant. Copyright © 2011 John Wiley & Sons, Ltd.     

Analytical definition of collapse surface in multiaxial space as a criterion for flow liquefaction occurrence

Prediction of flow liquefaction occurrence in multiaxial stress space is presented in this paper. The closed-form relation analytically attained by means of Dafalias-Manzari constitutive equations gives the flow stress ratio as a function of model parameters, state parameter, direction of undrained loading and especially initial condition of consolidation. These proposed flow condition can predict the occurrence of flow liquefaction initiated by any loadings whether the soil is consolidated isotropically or anisotropically. Suggested predictions are compared with the results of triaxial tests performed on Hostun RF sand. The results show that increasing consolidation stress ratio leads to flow stress ratio increase.     

An elastoplastic constitutive model for frozen saline coarse sandy soil undergoing particle breakage

The mechanical property of frozen saline sandy soil is complicated due to its complex components and sensitivity to salt content and confining pressure. Thus, a series of triaxial compression tests were carried out on sandy samples with different Na₂SO₄ contents under different confining pressures to explore the effects of particle breakage, pressure melting, shear dilation and strain softening or hardening. The test results indicate that the stress–strain curves exhibit strain softening/hardening phenomena when the confining pressures are below or above 6 MPa, respectively. A shear dilation phenomenon was observed in the loading process. With increasing confining pressure, the strength firstly increases and then decreases. By taking into consideration the changes between the grain size distributions before and after triaxial compression tests, a failure strength line incorporating the influences of both particle breakage and pressure melting is proposed. In order to describe the deformation characteristics of frozen saline sandy soil, an elastoplastic incremental constitutive model is established based on the test results. The proposed model considers the plastic compressive, plastic shear and breakage mechanisms by adopting the non-associated flow rule. The breakage mechanism can be reflected by an index related to the initial, current and ultimate grain size distributions. The hardening parameters corresponding to compressive and shear mechanisms consider the influence of particle breakage. Then the effect of particle breakage on both the stress–strain and volumetric strain curves is analyzed. The calculated results fit well with the test results, indicating that the developed constitutive model can well describe the mechanical and deformation features of frozen saline sandy soil under various stress levels and stress paths. In addition, the strain softening/hardening, contraction, high dilation and particle breakage can be well captured.

     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

Modeling cyclic shearing of sands in the semifluidized state

Liquefaction is associated with the loss of mean effective stress and increase of the pore water pressure in saturated granular materials due to their contractive tendency under cyclic shear loading. The loss of mean effective stress is linked to loss of grain contacts, bringing the granular material to a “semifluidized state” and leading to development and accumulation of large cyclic shear strains. Constitutive modeling of the cyclic stress-strain response in earthquake-induced liquefaction and post-liquefaction is complex and yet very important for stress-deformation and performance-based analysis of sand deposits. A new state internal variable named strain liquefaction factor is introduced that evolves at low mean effective stresses, and its constitutive role is to reduce the plastic shear stiffness and dilatancy while maintaining the same plastic volumetric strain rate in the semifluidized state. This new constitutive ingredient is added to an existing critical state compatible, bounding surface plasticity reference model, that is well established for constitutive modeling of cyclic response of sands in the pre-liquefaction state. The roles of the key components of the proposed formulation are examined in a series of sensitivity analyses. Their combined effects in improving the performance of the reference model are examined by simulating undrained cyclic simple shear tests on Ottawa sand, with focus on reproducing the increasing shear strain amplitude as well as its saturation in the post-liquefaction response.     

Predicting flow liquefaction,a constitutive model approach

In this paper, flow liquefaction criterion for contractive loose sands is analytically extracted based on the fundamental definition of flow liquefaction. In order to obtain the closed form of this criterion, Dafalias–Manzari constitutive model is employed; so the stress ratio at the onset of flow liquefaction is presented as a function of model parameters, state parameter and void ratio. Flow liquefaction line, as a graphical form of suggested criterion in stress space, shows that the peak points of undrained stress paths with same void ratios are not necessarily in a straight line. In order to validate the reliability of proposed flow liquefaction line to predict the onset of instability, it has been compared with the results of experimental tests performed on Toyoura, Ottawa and Leighton Buzzard sands. The verification results show that the present criterion can satisfactorily predict the onset of flow liquefaction in monotonic and cyclic undrained tests of saturated sands as well as the structural collapse in constant deviatoric stress tests of loose dry sands.     

波浪作用下黄河三角洲粉质土海床动力响应分析

应用饱和土动力固结理论和Pastor-Zienkiewicz III动力本构关系,对波浪作用下黄河三角洲粉质土海床的动力响应特征进行了有限元分析,应用总超孔压准则对海床进行了液化判别,并将计算结果与现场观测资料进行对比。结果表明：波浪导致的海床超孔压由瞬态孔压和累积孔压两部分组成;相比均质海床,拥有表面硬层的海床瞬态孔压沿深度衰减更快,累积孔压在表层增长速度更大;不同波浪条件下,瞬态孔压值及其变化趋势较为一致,累积孔压则具有较大的不同。年平均波浪条件下海床不会发生液化;5 a和50 a一遇极端波浪条件下,考虑三维效应和具有表面硬层的海床更容易液化,最大液化深度在海床表面以下2～3 m范围内。计算所得的海床最大液化深度与实测的黄河三角洲海底灾害地貌深度有较好的一致性,表明了文中方法的有效性和合理性。     

考虑土体结构性的弹塑性硬化模型

在对人工制备结构性土样等应力比压缩试验结果分析基础上,确定出结构性土体初始屈服面形状和土体初始屈服后塑性应变增量的方向,推导出结构性土体屈服函数的表达式;硬化参数采用塑性功的函数,根据三轴排水剪切试验结果确定出土体的硬化规律。由此构建能反映土体结构性的弹塑性硬化本构模型,并用试验进行了验证。本文提出一种基于试验的建模方法,不依赖经典塑性力学理论的正交流动规则,并建立可考虑土体结构性影响的本构模型,对土体结构性研究具有借鉴意义。     

A modification to dense sand dynamic simulation capability of Pastor–Zienkiewicz–Chan model

A modification to the nonlinear Pastor–Zienkiewicz–Chan (PZC) constitutive model without any change in the number of model parameters is introduced in order to simulate stiffness degradation of dense sands at dynamic loading. The PZC model is based on generalized plasticity and was verified by good prediction of liquefaction and undrained behavior of saturated sand. The PZC is a robust model that can predict drained dynamic behavior of sands, especially stiffness increase in loose sand at reloading of dynamic loading. Yet, this model does not show stiffness degradation of dense sand at reloading. The modification is made through modifying the stress memory factor, H _DM, which is multiplied by the plastic modulus, H _L. This modification does not influence reloading behavior of loose sand. The modified PZC model is verified via results of drained cyclic tests. Two cyclic triaxial tests on loose and dense specimens, along with two cyclic plane strain tests on dense sand are utilized for validation. The model simulation shows that the modified PZC model is able to predict the stiffness degradation of dense sand at reloading well.     

Discrete element modeling of static liquefaction of shield tunnel face in saturated sand

The failure mechanism and limit support pressure are essential factors for the stability assessment of tunnel face during shield tunneling. Previous studies were mostly based on the assumption that soil fails when the plastic limit is reached. The static liquefaction of saturated sand at pre-failure state often largely reduces shear strength; a stable tunnel face could be destabilized and even collapses abruptly. This paper aims to explore the static liquefaction failure of tunnel face in saturated sands by 3D discrete element simulation. The particle shape is considered by using clump composed of two identical spheres, and the micro-material parameters are calibrated by fitting against triaxial tests. The tunnel face is initially supported by the earth pressure at rest, and then the pressure is reduced by applying a displacement on the free surface until soil body completely collapses. The second-order work during the unloading process is calculated, and a negative value denotes the initiation of static liquefaction. The distribution of negative second-order work point shows the initiation and propagation of static liquefaction. The ultimate failure mode is shown to be composed of a sliding wedge and an overlying chimney, and the chimney is confined to a local area for soil arching. The experienced stress path of tunnel face verifies the existence of static liquefaction instability, indicating the stability analysis should consider static liquefaction rather than simply using the conventional plastic limit state analysis.

     

Modified Structured Cam Clay: A generalised critical state model for destructured,naturally structured and artificially structured clays

This paper presents a generalised constitutive model for destructured, naturally structured and artificially structured clays that extends the Structured Cam Clay (SCC) model. This model is designated as “Modified Structured Cam Clay (MSCC) model”. The influence of structure and destructuring on the mechanical behaviour of clay can be explained by the change in the modified effective stress, which is the sum of the current mean effective stress and the additional mean effective stress due to structure (structure strength). The presence of structure increases the modified mean effective stress and yield surface, enhancing the cohesion, peak strength and stiffness. The destructuring begins when the stress state is on the virgin yield surface. After the failure (peak strength) state, the abrupt destructuring occurs as the soil–cementation structure is crushed; hence the strain softening. The soil structure is completely removed at the critical state when the yield surface becomes identical to the destructured surface. The destructuring law is proposed based on this premise. In the MSCC model, the yield function is the same shape as that of the Modified Cam Clay (MCC) model. A plastic potential is introduced so as to account for the influence of structure on the plastic strain direction for both hardening and softening behaviours. The required model parameters are divided into those describing destructured properties and those describing structured properties. All the parameters have physical meaning and can be simply determined from the conventional triaxial tests. Thus, the MSCC model is a useful tool for geotechnical practitioners. The capability of the model is verified by the test results of destructured, natural structured and artificially structured clays.     

Constitutive analysis of the mechanical anisotropy of Opalinus Clay

This paper aims to analyse the anisotropic features of behaviour of Opalinus Clay using the theory of plastic multi-mechanisms. The results of triaxial tests conducted with different load levels and directions showed that the mechanical behaviour of this shale is cross-anisotropic. The stiffer samples are those in which the loading direction is parallel to the bedding plane. This indicates that the preconsolidation stress depends on the orientation of the load with respect to the fabric of Opalinus Clay. It is proposed to interpret the observed cross-anisotropy with an elastoplastic model based on four plastic strain mechanisms that may be successively mobilised depending on the loading direction. The predicted stress–strain responses vary according to the directions of the space as a result of the hardening process, depending on the number of plastic strain mechanisms that have been mobilised. The numerical predictions show overall good agreement with the experimental data in terms of deviatoric stress versus axial strain, demonstrating that multi-mechanism plasticity is a suitable constitutive tool for the interpretation of the mechanical anisotropy of this shale.     

A two-mechanism soil–structure interface model for three-dimensional cyclic loading

A three-dimensional (3D) soil–structure interface model is proposed within the two-mechanism constitutive theory and bounding surface theory originally established for soils. The proposed model has two main characteristics: first, the model is formulated based on two different and superposed deformation mechanisms. The first mechanism is due to the stress ratio increment, and the second is due to the normal stress increment. Each mechanism induces a shear strain component and a normal strain component. The proposed model can be reduced to the conventional single-mechanism interface model. Second, the plastic modulus and stress dilatancy are defined using the bounding surface theory. The plastic flow rule under cyclic loading is modified and assumed to be dependent on both the stress state of the mapping point and the stress reversal loading direction. The proposed model was validated against the available 3D interface tests and was found to satisfactorily reflect the salient features of the interfaces under monotonic and cyclic loading paths with different normal boundaries. The problem in which the elastic normal stiffness in conventional single-mechanism interface models is often underestimated to enhance the simulation performance under varying normal stress conditions is solved by incorporating the second mechanism. And the effect of the second mechanism on the modeling behavior is discussed. The modified plastic flow direction accurately simulates the 3D cyclic shear response, and the difference between the model simulation and test result increases with the number of cycles by use of the plastic flow direction defined in conventional bounding surface theory.     

A kinematic hardening constitutive model for sands: the multiaxial formulation

This paper explores the possibility of using well-accepted concepts—Mohr-Coulomb-like strength criterion, critical state, existence of a small strain elastic region, hyperbolic relationship for representing global plastic stress–strain behaviour, dependence of strength on state parameter and flow rules derived from the Cam-Clay Model—to represent the general multiaxial stress–strain behaviour of granular materials over the full range of void ratios and stress level (neglecting grain crushing). The result is a simple model based on bounding surface and kinematic hardening plasticity, which is based on a single set of constitutive parameters, namely two for the elastic behaviour plus eight for the plastic behaviour, which all have a clear and easily understandable physical meaning. In order to assist the convenience of the numerical implementation, the model is defined in a ‘normalized’ stress space in which the stress–strain behaviour does not undergo any strain softening and so certain potential numerical difficulties are avoided. In the first part the multiaxial formulation of the model is described in detail, using appropriate mixed invariants, which rationally combine stress history and stress. The model simulations are compared with some experimental results for tests on granular soils along stress paths lying outside the triaxial plane over a wide range of densities and mean stresses, using constitutive parameters calibrated using triaxial tests. Furthermore, the study is extended to the analysis of the effects induced by the different shapes of the yield and bounding surfaces, revealing the different role played by the size and the curvature of the bounding surface on the simulated behaviour of completely stress- and partly strain-driven tests. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical modelling of the effects of foundation scour on the response of a bridge pier

Foundation scour can have a detrimental effect on the performance of bridge piers, inducing a significant reduction of the lateral capacity of the footing and accumulation of permanent settlement and rotation. Although the hydraulic processes responsible for foundation scour are nowadays well known, predicting their mechanical consequences is still challenging. Indeed, its impact on the failure mechanisms developing around the foundation has not been fully investigated. In this paper, numerical simulations are performed to study the vertical and lateral response of a scoured bridge pier founded on a cylindrical caisson foundation embedded in a layer of dense sand. The sand stress–strain behaviour is reproduced by employing the Severn-Trent model. The constitutive model is firstly calibrated on a set of soil element tests, including drained and undrained monotonic triaxial tests and resonant column tests. The calibration procedure is implemented considering the stress and strain nonuniformities within the samples, by simulating the laboratory tests as boundary value problems. The numerical model is then validated against the results of centrifuge tests. The results of the simulations are in good agreement with the experimental results in terms of foundation capacity and settlement accumulation. Moreover, the model can predict the effects of local and general scour. The numerical analyses also highlight the impact of scouring on the failure mechanisms, revealing that the soil resistance depends on the hydraulic scenario.

     

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.     

脉冲地震动中竖向加速度对场地液化的影响

为研究近断层脉冲地震动中竖向加速度对砂土场地液化的影响,基于有限元平台OpenSees开发的边界面塑性本构模型,建立了动单剪单元试验模型和饱和砂土三维有限元模型。选取台湾Chi-Chi地震中10条具有速度脉冲特性的地震波,对比分析了水平双向脉冲波与三向脉冲波作用下土柱竖向位移、循环应力比、孔压比及等效循环周数的差异性,继而明确了脉冲地震动中竖向加速度对砂土液化的影响规律。研究表明,三向脉冲地震波中竖向加速度分量对场地永久位移值影响较小,但使永久位移的发展持时明显增大;土柱循环应力比受竖向地震动影响较小,因此分析脉冲地震动对场地剪切特性的影响时,可将三向脉冲地震动简化为水平双向地震动;考虑竖向地震动的三向脉冲地震波引起的孔压比变化幅度较大,孔压消散时间较长;三向脉冲地震波对应的等效循环周数较大,地震动发展持时长,可认为竖向加速度对场地液化有促进作用。