Semi-analytical solutions to one-dimensional consolidation for unsaturated soils with semi-permeable drainage boundary

The semi-analytical solutions to Fredlund and Hasan's one-dimensional(1 D)consolidation for unsaturated soils with a semi-permeable drainage boundary are presented. Two variables are introduced to transform the two coupled governing equations of pore-water and pore-air pressures into an equivalent set of partial differential equations(PDFs), which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pressures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are conducted on the pore-air and pore-water pressures at different ratios(the air permeability coefficient to the water permeability coefficient) and depths.     

Semi-analytical solution to one-dimensional consolidation in unsaturated soils

In this paper, a series of semi-analytical solutions to one-dimensional consolidation in unsaturated soils are obtained. The air governing equation by Fredlund for unsaturated soils consolidation is simplified. By applying the Laplace transform and the Cayley-Hamilton theorem to the simplified governing equations of water and air, Darcy＇s law, and Fick＇s law, the transfer function between the state vectors at top and at any depth is then constructed. Finally, by the boundary conditions, the excess pore-water pressure, the excess pore-air pressure, and the soil settlement are obtained under several kinds of boundary conditions with the large-area uniform instantaneous loading. By the Crump method, the inverse Laplace transform is performed, and the semi-analytical solutions to the excess pore-water pressure, the excess pore-air pressure, and the soils settlement are obtained in the time domain. In the case of one surface which is permeable to air and water, comparisons between the semi-analytical solutions and the analytical solutions indicate that the semi-analytical solutions are correct. In the case of one surface which is permeable to air but impermeable to water, comparisons between the semi-analytical solutions and the results of the finite difference method are made, indicating that the semi-analytical solution is also correct.     

Analytical solution to one-dimensional consolidation in unsaturated soils

This paper presents an analytical solution of the one-dimensional consolidation in unsaturated soil with a finite thickness under vertical loading and confinements in the lateral directions. The boundary contains the top surface permeable to water and air and the bottom impermeable to water and air. The analytical solution is for Fredlund＇s one-dimensional consolidation equation in unsaturated soils. The transfer relationship between the state vectors at top surface and any depth is obtained by using the Laplace transform and Cayley-Hamilton mathematical methods to the governing equations of water and air, Darcy＇s law and Fick＇s law. Excess pore-air pressure, excess pore-water pressure and settlement in the Laplace-transformed domain are obtained by using the Laplace transform with the initial conditions and boundary conditions. By performing inverse Laplace transforms, the analytical solutions are obtained in the time domain. A typical example illustrates the consolidation characteristics of unsaturated soil from an- alytical results. Finally, comparisons between the analytical solutions and results of the finite difference method indicate that the analytical solution is correct.     

Exact solutions of multi-term fractional diffusion-wave equations with Robin type boundary conditions

General exact solutions in terms of wavelet expansion are obtained for multi- term time-fractional diffusion-wave equations with Robin type boundary conditions. By proposing a new method of integral transform for solving boundary value problems, such fractional partial differential equations are converted into time-fractional ordinary differ- ential equations, which are further reduced to algebraic equations by using the Laplace transform. Then, with a wavelet-based exact formula of Laplace inversion, the resulting exact solutions in the Laplace transform domain are reversed to the time-space domain. Three examples of wave-diffusion problems are given to validate the proposed analytical method.     

Exact solutions of multi-term fractional difusion-wave equations with Robin type boundary conditions

General exact solutions in terms of wavelet expansion are obtained for multiterm time-fractional difusion-wave equations with Robin type boundary conditions. By proposing a new method of integral transform for solving boundary value problems, such fractional partial diferential equations are converted into time-fractional ordinary diferential equations, which are further reduced to algebraic equations by using the Laplace transform. Then, with a wavelet-based exact formula of Laplace inversion, the resulting exact solutions in the Laplace transform domain are reversed to the time-space domain. Three examples of wave-difusion problems are given to validate the proposed analytical method.     

Generalized solutions of transient thermal shock problem with bounded boundaries

In this work, the generalized thermoelastic solutions with bounded boundaries for the transient shock problem are proposed by an asymptotic method. The governing equations are taken in the context of the generalized thermoelasticity with one relaxation time (L–S theory). The general solutions for any set of boundary conditions are obtained in the physical domain by the Laplace transform techniques. The corresponding asymptotic solutions for a thin plate with finite thickness, subjected to different sudden temperature rises in its two boundaries, are obtained by means of the limit theorem of Laplace transform. In the context of these asymptotic solutions, two specific problems with different boundary conditions have been conducted. The distributions of displacement, temperature and stresses, as well as the propagations, intersections and reflections of two elastic waves, named as thermoelastic wave and thermal wave separately, are obtained and plotted. These results are agreed with the results obtained in the existing literatures.     

动态载荷下功能梯度复合材料的圆币形裂纹问题

研究了动态载荷下功能梯度材料中的圆币形裂纹问题．假设材料为横观各向同性,并且含有多个垂直于厚度方向的裂纹,材料参数沿轴向（与裂纹面垂直的方向）为变化的,沿该方向将材料划分为许多单层,各单层材料参数为常数,利用Hankel变换祛,在Laplace域内推导出了控制问题的对偶积分方程组．利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子和能量释放率．研究了含两个裂纹的功能梯度接头结构,分析了材料非均匀性参数对应力强度因子和能量释放率的影响．     

Dynamic response for functionally graded materials with penny-shaped cracks

This paper provides a method for studying the penny-shaped cracks configuration in functionally graded material(FGM) structures subjected to dynamic or steady loading. It is assumed that the FGMs are transversely isotropic and all the material properties only depend on the axial coordinatez. In the analysis, the elastic region is treated as a number of layers. The material properties are taken to be constants for each layer. By utilizing the Laplace transform and Hankel transform technique, the general solutions for the layers are derived. The dual integral equations are then obtained by introducing the mechanical boundary and layer interface conditions via the flexibility/stiffness matrix approach. The stress intensity factors are computed by solving dual integral equations numerically in Laplace transform domain. The solution in time domain is obtained by utilizing numerical Laplace inverse. The main advantage of the present model is its ability for treating multiple crack configurations in FGMs with arbitrarily distributed and continuously varied material properties by dividing the FGMs into a number of layers with the properties of each layer slightly different from one another. This work was supported by Failure Mechanics Laboratory of State Education Commission and the Post-doctor Research Fund of China.     

考虑在地基固结和流变的上部结构和地基基础的共同作用问题

本文提出一种能同时考虑地基的瞬时变形,团结变形和流变的地基基础与上部结构共同作用的一种新的计算方法,在这种方法中,土的计算考虑了Biot固结理论和流变,上部结构利用子结构法进行凝聚可得边界刚度矩阵和边界力,利用位移协讯条件则可以得到共同作用的控制方程,对上述方程进行Laplace变换,可以得到Laplace变换域内的控制方程,在Laplace变换域内对上述方程进行数值求解并进行相应的Laplace逆变换则可以得到时间境内任意时刻的解。     

Dynamic Analysis of an Infinite Cylindrical Hole in a Saturated Poroelastic Medium

In this paper, the dynamic response of an infinite cylindrical hole embedded in a porous medium and subjected to an axisymmetric ring load is investigated. Two scalar potentials and two vector potentials are introduced to decouple the governing equations of Biot’s theory. By taking a Fourier transform with respect to time and the axial coordinate, we derive general solutions for the potentials, displacements, stresses and pore pressures in the frequency-wave-number domain. Using the general solutions and a set of boundary conditions applied at the hole surface, the frequency-wave-number domain solutions for the proposed problem are determined. Numerical inversion of the Fourier transform with respect to the axial wave number yields the frequency domain solutions, while a double inverse Fourier transform with respect to frequency as well as the axial wave number generates the time-space domain solution. The numerical results of this paper indicate that the dynamic response of a porous medium surrounding an infinite hole is dependant upon many factors including the parameters of the porous media, the location of receivers, the boundary conditions along the hole surface as well as the load characteristics.     

Transient analysis of viscoelastic helical rods subject to time-dependent loads

The transient analysis of viscoelastic helical rods subject to time-dependent loads are examined in the Laplace domain. The governing equations for naturally twisted and curved spatial rods obtained using the Timoshenko beam theory are rewritten for cylindrical helical rods. The curvature of the rod axis, effect of rotary inertia and, shear and axial deformations are considered in the formulation. The material of the rod is assumed to be homogeneous, isotropic and linear viscoelastic. The viscoelastic constitutive equations are written in the Boltzmann–Volterra form. Ordinary differential equations in canonical form obtained in the Laplace domain are solved numerically using the complementary functions method to calculate the dynamic stiffness matrix of the problem. The solutions obtained are transformed to the real space using an appropriate numerical inverse Laplace transform method. Numerical results for quasi-static and dynamic response of viscoelastic models are presented in the form of graphics.     

A new method for solving Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system

A new method is developed to solve Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system. Based on the governing equations of Biot＇s consolidation and the technique of Laplace transform, Fourier expansions and Hankel transform with respect to time t, coordinate θ and coordinate r, respectively, a relationship of displacements, stresses, excess pore water pressure and flux is established between the ground surface （z = 0） and an arbitrary depth z in the Laplace and Hankel transform domain. By referring to proper boundary conditions of the finite soil layer, the solutions for displacements, stresses, excess pore water pressure and flux of any point in the transform domain can be obtained. The actual solutions in the physical domain can be acquired by inverting the Laplace and the Hankel transforms.     

Thermoelastic transient response of an infinitely long annular multilayered cylinder

Thermoelastic transient response of multilayered annular cylinders of infinite lengths subjected to known temperature at traction-free inner and outer surfaces are considered. A method based on the Laplace transformation and finite difference method has been developed to analyze the thermoelasticity problem. Using the Laplace transform with respect to time, the general solutions of the governing equation are obtained in transform domain. The solution is obtained by using the matrix similarity transformation and inverse Laplace transform. Solutions for the temperature and thermal stress distributions in a transient state were obtained. It was found that the temperature distribution, the displacement and the thermal stresses change slightly as time increases. There is no limit of number of annular layers of the cylinder in the presented computational procedures.     

Dynamic analysis of beams on viscoelastic foundation

This study is intended to analyze dynamic behavior of beams on Pasternak-type viscoelastic foundation subjected to time-dependent loads. The Timoshenko beam theory is adopted in the derivation of the governing equation. Ordinary differential equations in scalar form obtained in the Laplace domain are solved numerically using the complementary functions method to calculate exactly the dynamic stiffness matrix of the problem. The solutions obtained are transformed to the real space using the Durbin's numerical inverse Laplace transform method. The dynamic response of beams on viscoelastic foundation is analyzed through various examples.     

Two-dimensional transient thermo-hydro-mechanical fundamental solutions of multiphase porous media in frequency and time domains

In this paper, the closed form two-dimensional fundamental solutions for a non-isothermal unsaturated deformable porous medium have been derived for a symmetric polar domain in both Laplace transform and time domains. The governing differential equations of the non-isothermal unsaturated soil consist of equilibrium, moisture, air and heat transfer equations including the suction effect, temperature effect and dissolved air in water. The derived fundamental solution has been verified mathematically by comparison with the previously presented corresponding fundamental solutions in three limiting cases including the steady-state thermo-hydro-mechanical, steady-state hydro-mechanical and elastostatic fundamental solutions. Also these 2D kernel functions are tested in comparison with a finite element method (FEM).     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

Nonsteady flow of groundwater in aquifer system with viscoelastic properties

The transient flow behavior of groundwater in aquifer-aquitard system with viscoelastic properties is studied. On the basis of previous works (Hantush, Neuman, Brutsaert, Corapcioglu), the new partial differential-integral equations are derived. The well-known equations (Hantush, Brutsaert) are the special cases of the new equations. The new equations describe the flow of a slightly compressible groundwater in layers with viscoelastic properties.Analytical solutions of the partial differential-integral equations are obtained by using the method of Laplace transform. The viscoelastic properties enhance the heterogeneities of elastic aquifer system which have delay and feed qualities. The agreements between the numericl inversion results of Laplace transform and the analytical solutions are good. The formulae predict the transient flow behavior of groundwater in this heterogeneous layers.     

时域边界元法分析撞水响应

基于势流理论，考虑流场的可压缩性，首先利用积分变换导得了势流问题的一个动力学倒易定理，在此基础上，进而求得问题对应的时空边界积分方程，然后通过对边界和时间轴同时离散，建立了一组有递推形式的时间边界元方程最后结合液面条件和物体运动方程耦全求解得到了刚体的撞水响应。     

层状饱和土Biot固结问题状态空间法

针对饱和多孔介质空间非轴对Biot固结问题，引入状态变量，构造了两组相比独立的状态变量方程，利用Fourier级数和Laplace-Hankel变换，将状态变量方程转换为两组一阶常微分方程组，提出了均质饱和多孔介质空间非轴对称Biot固结问题的传递矩阵，得到以状态变量和传递矩阵乘积的形式表示的均质饱和多孔介质空间非轴对称Biot固结问题的解，利用层间完全接触的条件，可得到N层饱和多孔介质空间非轴对称Biot固结问题的一般解析表达式，文中考虑几种不同的边界条件，分析了两个算例，数值结果表明该方法具有较高的计算精度和良好的计算稳定性。     

理想导体中具有两个松弛时间的电磁热弹性波

研究处于均布磁场中的理想导体的二维电磁热弹性耦合问题，引入势函数使控制方程转化为3个偏微分方程．运用Laplace变换和Fourier变换得到该问题在变换域内的精确表达式，再通过级数展开和Laplace逆变换法求得在时间较短时的逆变换，得到时间-空间域内问题的解．运用此方法研究了表面受到热冲击的半无限空间问题．给出了电磁热弹性波、膨胀波和横向波传播的速度，并通过数值计算，给出了各个场量的分布图．所得结论与已有的结论一致．