3D dynamic response of a transversely isotropic multilayered medium subjected to a moving load

The dynamic problem of a transversely isotropic multilayered medium is reducible to quasi‐static problem by introducing a moving system that travels synchronously with the load. Based on the governing equations in the moving system, the ordinary differential equations in the Fourier transformed domain are deduced. An extended precise integration method is adopted to solve the ordinary differential equations, and the solution in the physical domain is recovered by the inverse Fourier transform. Numerical examples are performed to verify the accuracy of the presented method and to analyze the influence of material properties and the load characteristic.     

层状横观各向同性地基平面应变问题的扩展精细积分解

利用扩展精细积分法求解横观各向同性地基的平面应变问题。扩展精细积分法具有高精度和较高计算效率的特点,是求解微分方程的有效方法,相比于解析法可以节省大量理论推导工作量。从直角坐标下弹性力学控制方程出发,推导出Fourier变换域内地基的常微分矩阵方程;之后对地基微层元进行消元合并,进一步得到荷载作用在地基内部时层状地基的扩展精细积解。与已有文献的对比验证了方法的精确性,并分析了横观各向同性参数、层状性质和荷载作用点对计算结果的影响。结果表明：土体竖向位移随着横观各向同性参数n的增大而减小,而随着横观各向同性参数m的增大而增大;荷载作用点 的变化只对作用点以上的土体有影响,而上层土体的模量对竖向位移计算结果的影响更为显著,土体成层性对沉降的影响要比对竖向应力的影响更为显著。     

横观各向同性层状地基动应力响应的求解与分析

动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。     

Analytical layer element solutions for deformations of transversely isotropic multilayered elastic media under nonaxisymmetric loading

This paper presents the analytical layer element solutions for deformations of transversely isotropic elastic media subjected to nonaxisymmetric loading at an arbitrary depth. The state vectors for the nonaxisymmetric problem are deduced through the substitution of the Hu Hai‐chang solutions into the basic equations for the transversely isotropic elastic media. From the state vectors, the analytical layer element of a single layer is obtained in the Hankel transformed domain. The analytical layer element is an exact and symmetric stiffness matrix whose elements are without positive exponential functions, which can not only simplify the calculation but also improve the stability of computation. On the basis of the continuity conditions between adjacent layers, the global stiffness matrix is obtained by assembling the interrelated layer elements. The solutions for the multilayered elastic media in the transformed domain are obtained by solving the algebraic equation of the global stiffness matrix, which satisfies the boundary conditions. The actual solutions in the physical domain are further obtained by inverting the Hankel transform. Finally, some cases are analyzed to verify the solutions and evaluate the influences of the transversely isotropic character and stratified character of the media on the load–displacement responses. The numerical results show that the variations of the elastic properties between layers have a great effect on the displacements of the multilayered media. Copyright © 2014 John Wiley & Sons, Ltd.     

Extended precise integration method for axisymmetric thermo‐elastic problem in transversely isotropic material

This paper presents a numerical solution for the analysis of the axisymmetric thermo‐elastic problem in transversely isotropic material due to a buried heat source by means of extended precise integral method. By virtue of the Laplace–Hankel transform applied into the basic governing equations, an ordinary differential matrix equation is achieved, which describes the relationship between the generalized stresses and displacements in transformed domain. An extended precise integration method is introduced to solve the aforementioned matrix equation, and the actual solution in the physical domain is acquired by inverting the Laplace–Hankel transform. Numerical examples are carried out to demonstrate the accuracy of the proposed method and elucidate the influence of the character of transverse isotropy, the anisotropy of linear expansion coefficient, the anisotropy of thermal diffusivity, and medium's stratification on the thermo‐elastic response. Copyright © 2015 John Wiley & Sons, Ltd.     

Static analysis of vertically loaded pile group in multilayered transversely isotropic soils

In this paper, a coupling approach is presented to study the static responses of vertically loaded pile group embedded in multilayered transversely isotropic soils. The individual pile in pile group is modeled by the finite element method, while the analytical layer-element method is applied to represent the soil's behavior. Then, the interaction equation of piles and soils is obtained by considering the force equilibrium and displacement compatibility conditions and solved by a FORTRAN program. The results computed by the proposed approach compare favorably with those from some existing solutions and field test. Some typical parametric analysis cases are investigated to study the effect of soil anisotropy, pile stiffness ratio, and pile spacing on the behavior of vertically loaded pile group.     

三维横观各向同性成层地基的传递矩阵解

通过解耦变换推导出三维直角坐标系下横观各向同性地基的非耦合状态方程;利用双重Fourier变换以及Cayley-Hamilton定理得到了单层地基的传递矩阵;结合边界条件和层间连续条件进而得其传递矩阵解。编制了相应程序并进行了数值计算与分析,结果表明：数值结果与已有文献结果十分吻合,地基的横观各向同性性质与成层性质对受荷地基中竖向位移和应力的影响较为显著。     

列车荷载下上覆弹性层横观各向同性饱和地基的动力响应

基于Biot多孔弹性介质的波动理论,构造了轨道、道碴、枕木及弹性层的横观各向同性饱和地基在列车荷载下的动力计算模型。利用Fourier变换,得到了列车荷载作用下横观各向同性软土地基上弹性层动力响应的解析结果。利用离散Fourier逆变换得到数值计算结果,分析了荷载速度、地基的各向异性参数、弹性层刚度系数及厚度对位移和孔压响应的影响。分析结果表明：弹性层对控制地基振动作用显著,地表振动幅值随荷载速度的增加而增大,软土的横向弹性模量对地表振动及土中孔压有较大影响。     

多层横观各向同性地基轴对称固结的传递矩阵解

从横观各向同性地基轴对称Biot固结的基本方程出发,通过关于t的Laplace变换和关于r的Hankel变换,得到关于z的一阶常微分方程组。然后,对变换域内的基本未知量进行线性化处理,建立了变换域内的基本状态变量在z = 0处和任意深度处z的显式关系。利用传递矩阵法,结合层间连续性条件和边界条件,得到了多层横观各向同性地基的Biot固结轴对称问题的解答。该解答能避免随着层数增加而需要求解大型方程组的困难,明显地提高了计算效率。     

Triaxial behaviour of transversely isotropic materials: application to sedimentary rocks

Failure and long‐term behaviour of oriented solids are studied. Transversely isotropic materials are considered and a mathematical formulation that respect the material symmetry is developed and applied to model the triaxial behaviour of sedimentary rocks. Two failure criteria and a viscoplastic constitutive model that describe, respectively, triaxial failure and triaxial creep tests are presented and discussed. The application of the developed models to describe the mechanical behaviour of Tournemire shale shows that theoretical predictions are in good 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 © 2007 John Wiley & Sons, Ltd.     

On the strength of transversely isotropic rocks

Accurate prediction of strength in rocks with distinct bedding planes requires knowledge of the bedding plane orientation relative to the load direction. Thermal softening adds complexity to the problem since it is known to have significant influence on the strength and strain localization properties of rocks. In this paper, we use a recently proposed thermoplastic constitutive model appropriate for rocks exhibiting transverse isotropy in both the elastic and plastic responses to predict their strength and strain localization properties. Recognizing that laboratory‐derived strengths can be influenced by material and geometric inhomogeneities of the rock samples, we consider both stress‐point and boundary‐value problem simulations of rock strength behavior. Both plane strain and 3D loading conditions are considered. Results of the simulations of the strength of a natural Tournemire shale and a synthetic transversely isotropic rock suggest that the mechanical model can reproduce the general U‐shaped variation of rock strength with bedding plane orientation quite well. We show that this variation could depend on many factors, including the stress loading condition (plane strain versus 3D), degree of anisotropy, temperature, shear‐induced dilation versus shear‐induced compaction, specimen imperfections, and boundary restraints.     

横观各向同性地基上刚性矩形基础分析

对横观各向同性地基上刚性矩形板进行了求解。首先,利用表面受矩形均布荷载作用下的层状横观各向同性地基的位移解答,获得地基的柔度矩阵;然后,通过刚性矩形基础与层状横观各向同性地基的协调条件,建立刚性矩形基础与横观各向同性地基共同作用的方程,进而求得基础的地基反力。通过编制相应的程序,确定了合理的网格划分值;最后,进行算例分析,分析了地基横观各向同性性质、矩形刚性基础的长宽比以及地基分层性对地基反力的影响。分析表明：以上3种因素对地基反力有重要影响。     

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.     

列车荷载下轨道系统-层状横观各向同性饱和地基动力响应

基于Biot波动理论,构建列车荷载-轨道系统-双层状横观各向同性饱和地基模型,将模型分为上覆路轨系统和地层系统。对上覆路轨系统和地层系统处理,并利用双重Fourier变换技术,在变换域中将横观各向同性饱和地基动力响应的求解简化为求解一个6阶控制方程的特征值问题,进而得到了列车荷载作用下双层横观各向同性饱和地基力响应的解析结果。利用离散Fourier逆变换得到数值计算结果,重点分析了上下土层的刚度和泊松比对位移和孔隙水压力和剪切应力响应的影响,结果表明,上、下土层刚度差异对地基动力响应有较大影响,土层各向异性参数中模量的影响较泊松比大。计算结果可为软土路基加固深度的确定提供理论依据。     

有限元和有限层元研究横向同性土对建筑物沉降的影响

从常规三轴压缩试验得到横观各向同性土层的力学参数,综合9种来自于不同国家和地区的横观各向同性土层的土样参数,研究分析了横观各向同性土层参数对建筑物沉降的影响。作为典型例子,对上海一幢剪力墙结构的箱形基础与横观各向同性土进行共同作用分析,并与各向同性土层上的分析结果进行比较,结果表明：横观各向同性土层参数对沉降的影响有些不同,尤其是不排水加载条件确定的土层参数影响的结果大不相同。     

Fundamental solutions for a fluid‐saturated,transversely isotropic,poroelastic solid

The fundamental solutions were obtained for step‐like point forces acting in three orthogonal directions and an instantaneous fluid point source in a fluid‐saturated, porous, infinite solid of transversely isotropic elasticity and permeability. After expressing the governing equations in the form of matrix in the Laplace space, we employed Kupradze's method together with the triple Fourier transforms. This method reduces the simultaneous partial differential equations with respect to three displacement components and a pore fluid pressure to a differential equation in terms of only one potential scalar function, which can be operationally solved in the transformed space. After the Laplace inversion of the potential, the residue theorem was applied to its Fourier inverse transform with respect to one of the transformation variables. The Fourier transforms with respect to two other variables were rewritten into the Hankel transforms. Copyright © 2002 John Wiley & Sons, Ltd.     

The axisymmetric consolidation of a semi‐infinite transversely isotropic saturated soil

This paper presents an exact analytical solution to fully coupled axisymmetric consolidation of a semi‐infinite, transversely isotropic saturated soil subjected to a uniform circular loading at the ground surface. The analysis is under the framework of Biot's general theory of consolidation. First, the governing equations of consolidation are transformed into a set of equivalent partial differential equations with the introduction of two auxiliary variables. These partial differential equations are then solved using Hankel–Laplace integral transforms. Once solutions in the transformed domain have been obtained, the actual solutions in the physical domain for displacements and stress components of the solid matrix, pore‐water pressure and fluid discharge can be finally obtained by direct numerical inversion. The accuracy of the numerical solutions developed is confirmed by comparison with an existing exact solution for an isotropic and saturated soil that is a special case of the more general problem addressed. Numerical analyses are also presented to investigate the influence of the degree of material anisotropy on the consolidation settlement. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

横观各向同性含液饱和多孔介质黏弹性瞬态反应

在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。