Visualization of material stiffness in geomechanics analysis

This paper presents novel visualization techniques to simplify representation of the fourth‐order material stiffness tensor as a set of three‐dimensional geometric objects. Stiffness visualization aids in understanding the complex stiffness characteristics of highly non‐linear constitutive models including modelled material anisotropy and loading path dependent stiffness variation. Stiffness visualization is relevant for understanding the relationship of material stiffness to global behaviour in the analysis of a boundary value problem. The spherical pulse stiffness visualization method, developed in the acoustics field, is extended to visualize stiffness of geomaterials using three three‐dimensional objects. This method is limited to relatively simple constitutive models with symmetric stiffness matrices insensitive to loading magnitude and direction. A strain dependent stiffness visualization method is developed that allows the examination of material stiffness for a range of loading directions and is suitable for highly non‐linear and path dependent material models. The proposed stiffness visualization can be represented as 3‐D, 2‐D and 1‐D objects. The visualization technique is used to represent material stiffness and its evolution during simulated soil laboratory tests and deep excavation construction. Copyright © 2005 John Wiley & Sons, Ltd.     

A middle surface concept (MSC) model for saturated sands in general stress space

An elastoplastic constitutive model is proposed for saturated sands in general stress space using the middle surface concept (MSC). In MSC, different features of stress–strain response of a material are divided into different pseudo‐yield surfaces. The true‐yield surface representing the true response is established by using various links between the yield surfaces. In this MSC sand model, several well‐known features of sand response are represented by three different pseudo‐yield surfaces, which are developed in a simple and straightforward way. These features include the critical state behaviour, the effects of state parameter, unloading and reloading plastic deformation, the influence of fabric anisotropy, and phase transformation line related behaviour. Finally, the model predictions and test results are compared for two different types of sands under a variety of loading conditions and good comparisons are obtained. The application of MSC to saturated sand modelling shows the versatility of MSC as a general concept for modelling stress–strain response of materials. Copyright © 2006 John Wiley & Sons, Ltd.     

Finite element formulation and algorithms for unsaturated soils. Part I: Theory

This paper presents a complete finite‐element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress–strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time‐stepping scheme which automatically adjusts the time‐step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non‐linearity caused by suction‐dependent plastic yielding, suction‐dependent degree of saturation, and saturation‐dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress–strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction. Copyright © 2003 John Wiley & Sons, Ltd.     

Probabilistic seismic response and reliability assessment of isolated bridges

Bridge seismic isolation strategy is based on the reduction of shear forces transmitted from the superstructure to the piers by two means: shifting natural period and earthquake input energy reduction by dissipation concentrated in protection devices. In this paper, a stochastic analysis of a simple isolated bridge model for different bridge and device parameters is conducted to assess the efficiency of this seismic protection strategy. To achieve this aim, a simple nonlinear softening constitutive law is adopted to model a wide range of isolation devices, characterized by only three essential mechanical parameters. As a consequence of the random nature of seismic motion, a probabilistic analysis is carried out and the time modulated Kanai-Tajimi stochastic process is adopted to represent the seismic action. The response covariance in the state space is obtained by solving the Lyapunov equation for a stochastic linearized system. After a sensitivity analysis, the failure probability referred to extreme displacement and the mean value of dissipated energy are assessed by using the introduced stochastic indices of seismic bridge protection efficiency. A parametric analysis for protective devices with different mechanical parameters is developed for a proper selection of parameters of isolation devices under different situations.     

往复荷载作用下锈蚀钢筋混凝土柱黏结滑移模型及数值模拟研究

为实现地震作用下锈蚀钢筋混凝土柱精细化数值模拟分析,基于已有研究成果建立往复荷载作用下锈蚀钢筋与混凝土间的黏结滑移本构模型:结合课题组前期试验结果,采用ABAQUS有限元分析软件对建立的黏结滑移本构模型进行有效性验证,通过对数值计算结果与试验结果之间误差分析,进一步对黏结滑移模型中的摩擦黏结应力系数和退化系数进行修正,最终建立更为合理的锈蚀钢筋与混凝土间黏结滑移本构模型。通过数值计算结果与试验结果的再次比较,验证修正后黏结滑移本构模型的有效性。结果表明:修正后的锈蚀钢筋与混凝土间黏结滑移模型可更好地反映往复荷载作用下锈蚀钢筋混凝土柱的滞回性能。该成果可为地震作用下锈蚀钢筋混凝土结构的数值分析计算提供理论参考。     

冻土弹塑性本构模型研究现状

土体的本构模型是分析计算土工结构变形规律的关键。自剑桥模型提出以来,对于土体本构模型的研究,各国学者在不同环境条件下发展研究了适用于相应试验条件下的各类模型。随着寒区经济的发展,关于冻土力学性质的研究日益增多,许多学者借鉴常规融土的研究方法建立了不同条件下的冻土弹塑性本构模型。为了进一步理清冻土变形行为的特征,完善适用于冻土弹塑性行为的本构模拟理论和方法,作者总结和分析了各类冻土弹塑性本构模型的理论基础、建模方法、参数确定方法等内容,对进一步发展可准确描述冻土复杂力学行为的本构模型具有指导意义。     

含水量对高含冰量冻结砂土屈服面的影响

以往关于冻土屈服特性的研究很少考虑含水量的影响,但实际工程中经常遇到的是变含水量情况,因此本文通过-5.0℃条件下高含冰量冻结砂土的一系列三轴压缩试验,系统地研究了含水量对屈服特性的影响,并且由此建立了带有含水量参数的偏应力-剪应变增量型关系式。试验结果表明：随含水量的变化,应力-应变曲线类型有明显变化,即在不同的含水量区间,含水量对硬化规律有不同的影响特性,因此为了使屈服函数的形式更加简单和提高拟合准确度,对于屈服函数的具体形式应该根据塑性剪应变（硬化参数）和含水量不同区间进行分别确定;当塑性剪应变较小时（0.00%~1.00%）,随含水量的增大,偏应力逐渐增加,而当塑性剪应变较大时（大于1.00%）,随含水量的增大,偏应力有一个先减小后增大的变化过程,并且42.0%可以作为含水量转折点;通过分段硬化原则建立的带有含水量参数的屈服函数与偏应力-剪应变的增量型关系式的模拟值与试验点基本相吻合,这说明得到的屈服函数与偏应力-剪应变的增量型关系式可以用于不同含水量条件下屈服面和偏应力-剪应变曲线的预测。     

水平分层胶结充填体损伤本构模型及强度准则

以矿山分层充填为背景,针对水平分层充填体,提出了初始分层损伤、荷载损伤与总损伤的概念;利用损伤力学理论,考虑分层效应,建立了充填体损伤演化及本构模型;基于全微分方法构建了考虑分层效应的充填体强度准则。开展了不同分层充填体力学特性试验。研究结果表明：所建立的分层充填体损伤本构模型及强度准则理论曲线与试验曲线高度吻合,验证了模型的正确性;初始分层损伤随分层数增多,呈现 的二次多项式递增,分层效应与荷载耦合作用劣化了充填体总损伤;充填体总损伤率呈先增大后减小趋势。在峰值之前,分层数越多,则总损伤率越大;而峰值之后,分层数越多,则总损伤率越小,且总损伤率最大值随分层数增多而增大。     

考虑颗粒破碎效应的堆石料静动力本构模型

为合理反映颗粒破碎对堆石料力学特性的影响,基于试验结果分析,得出了堆石料在压缩和剪切作用下的颗粒破碎特性规律。通过引入压缩破碎和剪切破碎的相关参数,借鉴已有本构模型的合理定义,吸收临界状态理论和边界面理论的优点,发展了考虑颗粒破碎和状态相关的堆石料静动力统一弹塑性本构模型,并阐述了模型参数的确定方法。该模型不仅能够反映堆石料在静力荷载作用下的低压剪胀、高压剪缩、应变软化和硬化等特性,还能够反映在循环荷载作用下应力-应变的滞回特性和残余变形的累积效应。最后为验证模型的合理性,分别对堆石料的静力三轴和循环三轴试验进行了数值模拟预测,结果表明：模型预测与试验数据吻合良好,所提出的本构模型能够合理地描述颗粒破碎对堆石料静动力变形特性的影响。