首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Ai  Zhi Yong  Ye  Zi  Song  Xiaoyu  Wang  Lu Jun 《Acta Geotechnica》2019,14(4):1143-1160

We develop a new numerical model based on a precise integration method to investigate the coupled thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source. To obtain the relational matrices of the extended precise integration method, we first convert the governing equations of the problem into ordinary differential matrix equations through the Laplace–Hankel transform. Then, the cylindrical heat source is divided into a series of plane heat sources, and the plane temperature load term is added to the state vector between layer elements. By combining the layer elements, we build a layered transversely isotropic numerical model containing a cylindrical heat source in the transformed domain. Finally, we solve the model in the transformed domain and obtain the solution of the problem in the real domain through the Laplace–Hankel transform inversion. The accuracy of this method is verified by comparing the solutions with the results of the analytical method and the finite element method. Then, we study the influence of the anisotropy of thermal parameters, the embedded depth, the length/radius ratio, the type of heat source and the stratification of the medium on the thermo-mechanical coupled performance.

  相似文献   

2.
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.  相似文献   

3.
Axisymmetric consolidation of a poroelastic soil layer with a compressible fluid constituent induced by groundwater drawdown was studied based on Biot’s axisymmetric consolidation theory. Laplace and Hankel transforms were employed to solve the governing equation. Explicit analytical solutions are obtained in the Laplace–Hankel transform domain when groundwater drawdown is induced by a constant pumping well. Based on the solutions, numerical computations were performed to study the influences of the compressibility of the fluid constituent on the consolidation behavior of the soil layer.  相似文献   

4.
This paper presents the analytical layer‐element method to analyze the consolidation of saturated multi‐layered soils caused by a point sink by considering the anisotropy of permeability. Starting from the governing equations of the problem, the solutions of displacements and stresses for a single soil layer are obtained in the Laplace–Hankel transformed domain. Then, the analytical layer‐element method is utilized to further derive the solutions for the saturated multi‐layered soils in the transformed domain by combining with the boundary conditions of the soil system and continuity conditions between adjacent layers. The actual solutions in the physical domain can be acquired by the inversion of Laplace–Hankel transform. Numerical results are carried out to show the accuracy and stability of the proposed method and evaluate the influence of sink depth and anisotropic permeability on excess pore pressure and surface settlement. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

5.
With the aid of integral transform techniques, this paper presents an extended precise integration solution for thermal consolidation problems of a multilayered porous thermo-elastic medium with anisotropic thermal diffusivity and permeability due to a heat source. From the fundamental governing equations, ordinary differential equations are derived by employing Laplace–Hankel transforms. By applying the extended precise integration method, equations in the transformed domain can be solved, and the actual solutions are further obtained by adopting a numerical inverse transformation. The accuracy and feasibility of the proposed theory is demonstrated by contrastive analysis with existing studies. Finally, several examples are carried out to investigate the influence of heat source’s type, axial distance, burial depth of heat source, ratio of thermo-permeability, permeability anisotropy, thermal diffusivity anisotropy and stratification on the thermal consolidation process.  相似文献   

6.
After the application of a Laplace–Hankel transform, the governing equations of Biot’s consolidation are solved analytically by using the eigenvalue approach. Then the analytical layer-element of a single soil layer can be obtained in the transformed domain by synthesizing the generalized displacements and stresses, which are both expressed by six arbitrary constants. The elements of the analytical layer-element only contain negative exponential functions, which leads to a considerable improvement in computation efficiency and stability. The global stiffness matrix equation of multilayered soils is further obtained by assembling the interrelated layer-elements, and the actual solution is achieved by numerical inversion of the Laplace–Hankel transform after the solution in the transformed domain is obtained. Numerical examples are given to demonstrate the accuracy of this method and to study the influence of the layered soil properties and time history on the consolidation behavior.  相似文献   

7.
白冰 《岩土力学》2011,32(10):2901-2906
基于考虑热渗效应和等温热流效应的热-水-力耦合的线性热弹性固结控制方程,建立无限长空心圆柱饱和多孔介质热固结问题的一种理论求解方法。该方法先给出Laplace变换域上的解,然后,利用Stehfest法求其数值逆变换。该理论解考虑了空心圆柱体内、外透水界面随时间变化的外力和温度荷载耦合作用过程。最后,通过一算例分析了饱和多孔介质的热固结特征,给出其温度、孔压、位移和应力的演化规律  相似文献   

8.
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.  相似文献   

9.
This paper presents a stable and efficient method for calculating the transient solution of layered saturated media subjected to impulsive loadings by means of the analytical layer element method. Starting with the field equations based on Biot's linear theory for porous, fluid‐saturated media, and the seepage continuity equation, an analytical layer element for a single layer is established by applying Laplace‐Hankel integral transform. The global stiffness matrix in the transform domain for a layered saturated half‐space subjected to a transient circular patch loading is obtained by assembling the layer elements of each layer. The displacements in the time domain are derived by Laplace‐Hankel inverse transform of the global stiffness matrix. Numerical examples are conducted to verify the accuracy of the method and to demonstrate the influences of type of transient loading, buried depth of loading, permeability, and stratification of materials on the transient response of the multilayered saturated poroelastic media.  相似文献   

10.
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.  相似文献   

11.
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.  相似文献   

12.
The governing differential equations of unsaturated soils considering the thermo‐poro‐mechanical behaviour consist of equilibrium, moisture air and heat transfer equations. In this paper at first, following some necessary simplifications, the thermal three‐dimensional fundamental solution for an unsaturated deformable porous medium with linear elastic behaviour in Laplace transform domain is presented. Subsequently, the closed‐form time domain fundamental solutions are derived by analytical inversion of the Laplace transform domain solutions. Then a set of numerical results are presented, which demonstrate the accuracies and some salient features of the derived analytical transient fundamental solutions. Finally, the closed‐form time domain fundamental solution will be verified mathematically by comparison with the previously introduced corresponding fundamental solution. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

13.
基于Terzaghi一维固结理论,分析了考虑半透水边界条件的分数阶导数黏弹性饱和土层在随时间变化的任意荷载作用下一维固结问题。首先,应用Laplace变换联立求解饱和土层一维固结微分方程和分数阶Kelvin-Voigt黏弹性本构方程,推导出有效应力和沉降在Laplace变换域内的解析解,采用Crump方法进行Laplace逆变换,得到了时间域内的半解析解。然后将本文得到的半解析解分别退化为半透水边界条件下基于黏弹性假设的一维固结半解析解和双面透水边界条件下基于分数阶黏弹性假设的一维固结半解析解,结果与已有文献的半解析解相同,验证了本研究所提出解的可靠性。最后通过算例分别考察了半透水边界参数、分数阶黏弹性模型参数和荷载参数对饱和土层固结沉降的影响。研究表明,半透水边界条件参数、分数阶次与黏滞系数主要影响饱和土层固结的发展快慢,而饱和土层的最终沉降量主要受到土层压缩模量的影响;另外,饱和土层的固结规律与外荷载变化规律一致。  相似文献   

14.
This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

15.
Based on Fredlund’s one-dimensional consolidation equation for unsaturated soil, Darcy’s law and Fick’s law, a semi-analytical solution was presented to the free drainage well with a finite thickness under application of uniform vertical loading and the boundary of the top and bottom surfaces impermeable to water and air. According to the polar governing equations of water and air phases and the boundary and initial conditions, the excess pore-air and pore-water pressures and the soil layer settlement in the Laplace transformed domain are obtained by performing the Laplace transform and utilizing the Bessel functions. Crump’s method is used to perform the inversion of Laplace transform in order to obtain numerical solutions in the real time domain. Finally, a typical example is given to illustrate the changes in the excess pore-air and pore-water pressures and soil layer settlement with time factor at different ratios of air–water permeability coefficient and/or different distances from the well.  相似文献   

16.
An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

17.
This paper presents a theoretical approach to analyse coupled, linear thermoporoelastic fields in a saturated porous medium under radial and spherical symmetry. The governing equations account for compressibility and thermal expansion of constituents, heat sink due to thermal dilatation of water and thermal expansion of the medium, and thermodynamically coupled heatwater flow. It has been reported in the literature that thermodynamically coupled heat–water flows known as thermo-osmosis and thermal filtration have the potential to significantly alter the flow fields in clay-rich barriers in the near field of a underground waste containment scheme. This study presents a mathematical model and examines the effects of thermo-osmosis and thermal-filtration on coupled consolidation fields in a porous medium with a cavity. Analytical solutions of the governing equations are presented in the Laplace transform space. A numerical inversion scheme is used to obtain the time-domain solutions for a cylindrical cavity in a homogeneous or a non-homogeneous medium. A closed form time-domain solution is presented for a spherical cavity in a homogeneous medium. Selected numerical solutions for homogeneous and non-homogeneous media show a significant increase in pore pressure and displacements due to the presence of thermodynamically coupled flows and a negligible influence on temperature. © 1998 John Wiley & Sons, Ltd.  相似文献   

18.
In this paper, a coupling method between finite element and analytical layer‐elements is utilized to analyze the time‐dependent behavior of a plate of any shape and finite rigidity resting on layered saturated soils. Based on the integral transform techniques together with the aid of an order reduction method, an analytical layer‐element solution is derived from the governing equations for three‐dimensional Biot consolidation with respect to a Cartesian coordinate system and then extended to be the fundamental solution for the layered saturated soil under a point load. The Mindlin plate is modeled by eight‐noded isoparametric elements. The governing equations of the interaction between soil and plate in the Laplace‐Fourier transformed domain are deduced by referring to the coupling theory of FEM/BEM, and the final solution is obtained by applying numerical inversion. Numerical examples concerned with the time‐dependent response of a plate are performed to demonstrate the influence of soil and plate properties on the interaction process. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

19.
多层横观各向同性地基轴对称固结的传递矩阵解   总被引:2,自引:1,他引:1  
艾智勇  王全胜 《岩土力学》2009,30(4):921-925
从横观各向同性地基轴对称Biot固结的基本方程出发,通过关于t的Laplace变换和关于r的Hankel变换,得到关于z的一阶常微分方程组。然后,对变换域内的基本未知量进行线性化处理,建立了变换域内的基本状态变量在z = 0处和任意深度处z的显式关系。利用传递矩阵法,结合层间连续性条件和边界条件,得到了多层横观各向同性地基的Biot固结轴对称问题的解答。该解答能避免随着层数增加而需要求解大型方程组的困难,明显地提高了计算效率。  相似文献   

20.
汪磊  李林忠  徐永福  夏小和  孙德安 《岩土力学》2018,39(11):4142-4148
为描述饱和土体的流变特性,引入分数阶导数Kelvin-Voigt黏弹性模型,采用解析方法对半透水边界下的分数阶黏弹性饱和土一维固结特性进行了研究。分别对骤加恒载下饱和土一维固结微分方程和分数阶Kelvin-Voigt黏弹性本构方程进行Laplace变换,并联立求解得到了双边半透水边界条件下分数阶黏弹性饱和土在Laplace变换域内的解析表达式。通过Crump方法实现Laplace数值反演,得到时间域内的半解析解。将所得到的解分别退化为分数阶黏弹性饱和土一维固结半解析解和双边半透水黏弹性饱和土一维固结半解析解,结果与已有文献半解析解相同,验证了提出的双边半透水边界条件下分数阶黏弹性饱和土一维固结解的可靠性。通过算例考察了半透水边界条件和分数阶黏弹性饱和土参数对一维固结特性的影响。研究表明,双边半透水边界下分数阶黏弹性饱和土一维固结发展过程与半透水边界条件、分数阶次和黏滞系数有关,且土体的压缩模量对饱和土一维固结最终沉降量有显著影响。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号