一种横观各向同性岩体蠕变模型

层状岩体由于层理面或者结构面的存在,在力学上具有横观各向同性的特点,现有的各向同性蠕变模型难以全面反映横观各向同性岩体的蠕变力学特性。为构建能够反映横观各向同性岩体的三维蠕变模型,以能反映瞬时应变、减速蠕变和稳态蠕变特征的黏弹性Burgers模型为基础,基于常泊松比假定,在三维各向同性蠕变本构方程的基础上,按照算子替换的方法,将横观各向同性柔度矩阵代替各向同性柔度矩阵,并考虑了平行和垂直层理方向岩体蠕变力学行为的差异性,推导了横观各向同性岩体的三维蠕变本构方程。根据本构方程的特点,提出了根据平行和垂直方向岩体蠕变试验结果进行三维蠕变本构模型中蠕变参数的辨识方法。将提出的模型应用于三轴蠕变试验参数辨识,从而获得了一套完整的三维蠕变参数,并与试验结果进行对比分析,从而验证了所提模型与方法的合理性和有效性。进一步,指出了传统蠕变试验设计方案的局限性,给出了横观各向同性材料蠕变试验设计建议。研究成果为研究岩体三维蠕变机制提供了新思路,可对岩体蠕变试验设计提供相应的科研支撑。     

Mathematical model for isotropic compression behaviour of cemented soil cured under stress

A mathematical model is proposed to simulate isotropic compression test results obtained for specimens of silty sand mixed with different cement contents and cured at different void ratios. The ability of the model to represent the isotropic test results of specimens with different cement contents cured at high and low void ratios under high confining stresses was successfully checked. The physical meanings of the parameters established by the model were also evaluated. Finally, the applicability of the model to other artificially cemented soils was demonstrated.     

温度作用下岩石的热黏弹塑性模型研究

在西原体模型的基础上,通过引入弹簧的热膨胀系数,将温度效应引入到模型中,推导了单轴应力状态下微分形式的热粘塑性本构方程。通过蠕变曲线和卸载曲线的分析,考虑了温度效应的西原体模型当应力水平较低时为稳定蠕变,存在瞬时弹性、弹性后效和由温度引起的粘性流动;当应力水平较高时为不稳定蠕变,存在瞬时弹性、弹性后效和由温度和应力共同引起的粘性流动,该模型较全面地反映了岩石在温度作用下的蠕变性质。     

A simple model for isotropic non-linear elasticity of frictional materials

A simple isotropic nonlinear elastic model has been developed to describe the behaviour of granular materials immediately after stress reversal. Its three parameters can be measured in the same drained triaxial test.     

A 3D non-orthogonal elastoplastic constitutive model for transversely isotropic soil

Acta Geotechnica - The significant influence of transverse isotropy on mechanical behaviours of soil has long been recognised, but it is still difficult to capture well under 3D stress conditions....     

弹性矩形板下横观各向同性多层地基分析

利用弹性矩形板与多层地基表面的竖向位移协调条件与光滑接触条件,由横观各向同性多层地基应力与位移非耦合的传递矩阵解,推导出弹性矩形板下竖向应力和位移的解析解。在此基础上,编制了相应的程序,并进行了数值计算。计算结果表明：矩形板刚度对板底竖向位移及板中心下的竖向应力有着较为显著的影响;板底竖向位移及板中心下的竖向应力随着板刚度的增加而减小,相同荷载作用下横观各向同性地基与均匀各向同性地基模型的计算结果差异较大,实际工程中很有必要采用更符合土体性质的横观各向同性地基模型     

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

考虑温度影响的冻土蠕变本构模型

冰作为冻土的基本组成部分,对冻土蠕变的加速蠕变阶段有重要影响。温度通过影响冻土中冰的冻结与融化过程,及其黏塑性流动,引起冻土结构的强化与弱化,从而成为了决定冻土蠕变力学行为的关键因素之一。同时,外部应力也造成冻土的强化与弱化,影响着冻土的蠕变。通过引入硬化因子与损伤因子来考虑温度、应力造成的冻土材料强化与弱化。硬化因子H代表了蠕变过程中强化效应的大小,而损伤因子D则代表了由弱化效应造成的冻土材料相关参数的折减比例,进而提出了适用于冻土的改进西原蠕变本构模型。该模型预测值与试验数据的比较表明：改进的模型不仅能较好地描述初始蠕变阶段、稳定蠕变阶段,而且相比于传统模型,能更好地描述冻土加速蠕变阶段,具有合理性与一定的实用性。     

考虑重力侵蚀影响的分布式土壤侵蚀模型

黄土高原地区的土壤侵蚀具有水力侵蚀和重力侵蚀相伴发生的特点,在大多数针对黄土高原的侵蚀模型研究中未考虑重力侵蚀的影响,使得模型的模拟精度较差。针对黄土高原的侵蚀特点,采用量化影响重力侵蚀发生的主要因素,确定重力侵蚀发生的具体沟道栅格单元的方法,从而考虑了重力侵蚀对产输沙过程的影响。所建立的分布式土壤侵蚀模型以逐网格汇流的水文模型为基础,采用逐网格侵蚀输沙的模拟方式,能够模拟上方来水来沙对侵蚀输沙的影响。经小理河流域实测资料验证,模型具有一定的计算精度。     

Can unbiased source be retrieved from anisotropic waveforms by using an isotropic model of the medium?

Anisotropy is frequently present in geological structures, but usually neglected when source parameters are determined through waveform inversion. Due to the coupling of propagation and source effects in the seismic waveforms, such neglect of anisotropy will lead to an error in the retrieved source. The distortion of the mechanism of a double-couple point source located in an anisotropic medium is investigated when inverting waveforms using isotropic Green's functions. The anisotropic medium is considered to be transversely isotropic with six levels of anisotropy ranging from a fairly weak to rather strong anisotropy, up to about 24% in P waves and 11% in S waves. Inversions are based on either only direct P waves or both direct P and S waves. Two different algorithms are employed: the direct parametrization (DIRPAR, a nonlinear algorithm) and the indirect parametrization (INPAR, a hybrid scheme including linear and nonlinear steps) of the source. The orientation of the double-couple mechanism appears to be robustly retrieved. The inclination of the resulting nodal planes is very small, within 10° and 20° from the original solution, even for the highest degree of anisotropy. However, the neglect of anisotropy results in the presence of spurious isotropic and compensated linear-vector dipole (CLVD) components in the moment tensor (MT). This questions the reliability of non-double-couple components reported for numerous earthquakes.     

中主应力系数影响下的冻结砂土损伤本构模型

人工冻结法是饱水砂层开挖过程中常用的止水和临时支护方法,通过冻土损伤特性研究为冻土力学特性和冻结体稳定性分析奠定基础。为研究冻结砂土的损伤力学特性,在-5℃下进行了不同中主应力系数的冻结砂土三维室内试验。从冻土微元破坏服从Weibull随机分布的特点出发,将Drucker-Prager强度准则作为冻土微元统计分布变量,利用应变等价性假说,建立了三维应力状态下冻结砂土损伤本构模型;在此基础上,讨论模型参数F₀与m和中主应力系数的关系,对模型参数进行合理修正,建立中主应力系数影响下的冻结砂土损伤本构模型,并与试验结果进行对比。分析结果表明：参数F₀和m随着中主应力系数的增大呈现先减小后增大的趋势;参数F₀反映了冻结砂土的强度特性,参数m代表了冻结砂土的延性及脆性特征,考虑中主应力系数影响的冻结砂土损伤本构模型能很好地模拟冻结砂土应力-应变全过程曲线。研究成果为人工冻结法工程设计提供一定的理论依据。移动阅读     

On the effects of deformation induced anisotropy in isotropic materials

Effects of recoverable deformation induced anisotropy in the elastic stiffness of isotropic materials are described. In isotropic materials, thermodynamics predicts coupling of hydrostatic and deviatoric responses. It is shown that the coupling of the two responses is more significant than previously recognized in the literature. Properly accounting for the coupling of hydrostatic and deviatoric responses requires re‐evaluating elastic materials characterization data, allowing for the coupled response. The result is an apparent decrease in the pressure sensitivity of the elastic shear modulus. The decrease in the pressure sensitivity of the shear modulus leads to stress paths that are more tangential to the yield surface in stress space, resulting in an increase in predicted elastic strain at each step of an elastic–plastic stress update. Consequently, predicted plastic strains and, in particular, volumetric plastic strains, are smaller than if recoverable deformation induced anisotropy had been neglected. The result is an associated plasticity model, which appears to be non‐associated. Copyright © 2012 John Wiley & Sons, Ltd.     

黄土的湿载结构性本构模型研究

土结构性是土结构空间排列几何特征和联结特征所表现出的力学性质,是影响土力学特性的重要内在因素。它与土的密度、粒度、湿度相联系,能够综合反映土的强度、变形等力学特性的变化和发展。从综合反映增湿导致土结构内在主动变化和加载变形导致土结构被动变化两个方面出发,沿用土综合结构势的应力比结构性参数,首先建立压缩、剪切及增湿作用下黄土的结构性演化规律;其次,通过压剪湿作用下黄土的应力-应变特性、屈服破坏特性及其与结构性之间的关系,揭示结构性变化对土力学性质的影响;进而,建立了反映压剪湿作用下黄土结构性发展演化的湿载结构性本构模型。最后通过与试验数据、修正剑桥模型的比较分析,初步验证该模型的准确性。     

A mesoscopic model for the behaviour of concrete under high confinement

When impact loaded, concrete is submitted to high triaxial stresses. The experimental response of concrete under quasi‐static triaxial compression is studied using a triaxial press capable of applying a mean pressure greater than 1 GPa on cylindrical samples measuring 7 cm in diameter and 14 cm high. A numerical analysis of these previous experiments is performed herein at a mesoscopic scale. Concrete is modelled as a biphasic material consisting of a mortar (cement paste and fine aggregates) and roughly spherical aggregates (with a diameter exceeding 2 mm) whose characteristics are applied on a regular cubic finite element mesh. A damage‐plasticity model is then used to model the behaviour of mortar. An identification of model parameters on mortar samples and the subsequent comparison between numerical and experimental tests will be presented for hydrostatic and triaxial compression. Copyright © 2009 John Wiley & Sons, Ltd.     

A model for the irradiative coloration of smoky feldspar and the inhibiting influence of water

Radiation-induced smoky color and associatedelectron paramagnetic resonance (EPR) signals develop only in potassium feldspar (KAlSi₃O₈) free of structurally bound molecular water. Fluid inclusion water does not influence coloration. The integrated intensity of each of the four bands (11,600, 16,200, 19,100, and 27,200 cm^?1) in the optical absorption spectra are linearly correlated with the doubly-integrated intensity of a broad, asymmetric first derivative atg _eff=2.027 in EPR spectra. In microcline, the EPR pattern is resolved into an asymmetric six-line pattern atg _eff=2.024 and a single derivative atg _eff=2.009 which, based on analogy to alkali-silicate glass, are due respectively to [SiO₄/K⁺]²⁺ and a hole shared between two nonbonding oxygens on Si. We propose that structural water inhibits formation of smoky centers in feldspar by releasing atomic hydrogen during irradiation which destroys centers while diffusing towards a stable site.     

An approach to model the mechanical behavior of transversely isotropic materials

The present work proposes an approach to adapt existing isotropic models to transversely isotropic materials. The main idea is to introduce equivalence relations between the real material and a fictitious isotropic one on which one can take all the advantages of the well‐established isotropic theory. Two applications of this approach are presented here: a failure criterion and a damage model that takes into account the load‐induced anisotropy. In both cases, theoretical predictions are in agreement with the experimental data. In the present paper, the developed approach is applied to sedimentary rock materials; nevertheless, it can be generalized to any material that exhibits transverse isotropy. Copyright © 2015 John Wiley & Sons, Ltd.     

A constitutive model for coal-fouled ballast capturing the effects of particle degradation

Rail tracks undergo degradation owing to particle breakage and fouling of ballast by various fines including coal and subgrade soil. As the ballast becomes fouled, its strength and drainage capacity are compromised, sometimes resulting in differential settlement and reduced track stability. This paper demonstrates a continuum mechanics based framework to evaluate the detrimental effect of fines on the strength, deformation and degradation of coal-fouled ballast under monotonic loading. An elastoplastic constitutive model that considers the effect of fines content and energy consumption associated with particle breakage during shearing is presented. This multiphase constitutive model is developed within a critical state framework based on a kinematic-type yield locus and a modified stress-dilatancy approach. A general formulation for the rate of ballast breakage and coal particle breakage during triaxial shearing is presented and incorporated into the plastic flow rule to accurately predict the stress–strain response of coal-fouled ballast at various confining pressures. The behaviour of ballast at various levels of fouling is analysed and validated by experimental data.     

量子地震模型:对汶川地震的解释

从大量的地震参数分析可知,在由同一主震引发的所有的余震发生的位置都不同。在发生地震之前,地震的震中位置都是一些相互独立、互不连续且发生了弹性应变的单元个体,笔者称这些独立不连续的弹性应变单元称为应变量子。综合前人研究成果,笔者建立了一个地震模型。运用该模型对2008年5月12日发生在中国汶川地震作了解释:由于在龙门山断裂带周围形成了许多的应变量子,这些应变量子的形成阻碍了龙门山断裂的运动,产生滑移亏损,从而造成在地震前测量到的龙门山断裂带速度场变化很不显著。2008年5月12日14时28分时,在龙门山断裂周围已形成应变量子中的其中一应变量子最先达到它的最大储能,释放它所储存的应变能,引发了汶川Ms8.0地震。此次地震产生的地震波引发了龙门山断裂周围地壳应力的重分布,使得其他应变量子提前达到最大储能,释放出能量,于是触发成千上万次余震。此外,我们或许可以通过观测断裂滑移速率的变化情况来预测断裂周围应变量子的形成,从而来预测该断裂是否存在潜在地震的可能。     

A transversely isotropic thermo-poroelastic model for claystone: parameter identification and application to a 3D underground structure

A transversely isotropic thermo-poroelastic constitutive law is developed and implemented in the finite element code Code_Aster (EDF, France). It is then validated using an analytic solution for an inclined borehole in a transversely isotropic medium. A strategy for identifying the parameters of the transversely isotropic thermo-poroelastic model based on an inverse method is proposed on the basis of different laboratory tests. To demonstrate the efficiency and applicability of the model, it is then applied in a three-dimensional numerical model of an underground structure in a parameter sensitivity study. The results of the modelling highlight the importance of accounting for anisotropic phenomena when determining the dimensions of underground facilities. The whole approach is presented in the paper, from model development to application to 3D numerical modelling to an engineering case study.     

A numerical model for the transient analysis of offshore foundations under cyclic loading

A comprehensive numerical model for the analysis of offshore foundations under a general transient loading is presented here. The theoretical basis of the model lies on the Swansea formulation of Biot’s equations of dynamic poroelasticity combined with a constitutive model that reproduces key aspects of cyclic soil behaviour in the frame of the theory of generalised plasticity. On the practical side, the adoption of appropriate finite element formulations may prevent the appearance of spurious numerical instabilities of the pore pressure field. In this respect, the use of a coupled enhanced-strain element is here proposed. On the other hand, the practicality of the presented model depends ultimately on its computational efficiency. Some practical recommendations concerning the solution strategies, the matrix storage/handling procedures and the parallel multi-processor computation are here provided. Finally, the performance of the model with a benchmark study case and its practical application to analyse the soil–structure interaction of an offshore monopile under a realistic transient storm loading are discussed.