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.     

Uniform difference scheme for a singularly perturbed linear 2nd order hyperbolic problem with zeroth order reduced equation

In this paper a singularly perturbed linear second order hyperbolic problem with zeroth order reduced equation is discussed. Firstly, an energy inequality of the solution and an estimate of the remainder term of the asymptotic solution are given. Then an exponentially fitted difference scheme is developed in an equidistant mesh. Finally, uniform convergence in small parameter is proved in the sense of discrete energy norm.     

用DQM求解成层地基一维非线性固结问题

由于多层地基的一维非线性固结问题求解的复杂性,其解析解很难求得。本文基于Davis和Raymond一维非线性固结理论,利用DQM(Differential Quadrature Method)导了初始有效应力沿深度变化、任意边界条件、任意荷载作用下成层地基一维非线性固结的统一表达式,求得了孔压、有效应力和平均固结度的解答。通过解的收敛性分析讨论了DQM解的有效性。由于DQM解对于固结间题各种复杂条件具有统一的矩阵表达式,更便于编程计算和工程应用。最后,用本文解答对三层地基一维非线性固结问题进行了讨论。     

饱和轴对称黏弹性体的求解

应用Laplace变换和Hankel变换求解Blot固结方程,得出了饱和半无限空间黏弹性体在均布荷载作用下的应力、位移和超孔隙水压力的显式解答。进行Kelvin模型黏弹性算例分析,从而验证了本文理论的准确性和有效性。计算表明：超孔隙水压力随着深度的增长先增加后减小,体积应变在荷载作用初期有明显的增加;随着时间的增长,体积应变的增加速度放缓。     

One-dimensional consolidation of layered soils with impeded boundaries under time-dependent loadings

On the basis of Terzaghi‘s one-dimensional consolidation theory, the variation of effective stress ratio in layered saturated soils with impeded boundaries under timedependent loading was studied. By the method of Laplace transform, the solution was presented. Influences of different kinds of cyclic loadings and impeded boundaries conditions were discussed. Through numerical inversion of Laplace transform, useful illustrations were given considering several common time-dependent loadings. Pervious or impervious boundary condition is just the special case of the problem considered here. Compared with average index method, the results from the method illustrated are more accurate.     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

A multi-scale theory of swelling porous media: I. Application to one-dimensional consolidation

A theory is developed which describes flow in multi-scale, saturated swelling media. To upscale information, both the hybrid theory of mixtures and the homogenization technique are employed. In particular, a model is formulated in which vicinal water (water adsorbed to the solid phase) is treated as a separate phase from bulk (non-vicinal) water. A new form of Darcy's law governing the flow of both vicinal and bulk water is derived which involves an interaction potential to account for the swelling nature of the system. The theory is applied to the classical one-dimensional consolidation problem of Terzaghi and to verify Low's empirical, exponential, swelling result for clay at the macroscale.     

A new analytical solution to axisymmetric Blot's consolidation of a finite soil layer

A new analytical method is presented to study the axisymmetric Blot's consolidation of a finite soil layer. Starting from the governing equations of axisymmetric Blot's consolidation, and based on the property of Laplace transform, the relation of basic variables for a point of a finite soil layer is established between the ground surface (z= 0) and the depth z in the Laplace and Hankel transform domains. Combined with the boundary conditions of the finite soil layer, the analytical solution of any point in the transform domain can be obtained. The actual solution in the physical domain can be obtained by inverse Laplace and Hankel transforms. A numerical analysis for the axisymmetric consolidation of a finite soil layer is carried out.     

Analytical solution for bending beam subject to lateral force with different modulus

A bending beam, subjected to state of plane stress, was chosen to investigate. The determination of the neutral surface of the structure was made, and the calculating formulas of neutral axis, normal stress, shear stress and displacement were derived. It is concluded that, for the elastic bending beam with different tension-compression modulus in the condition of complex stress, the position of the neutral axis is not related with the shear stress, and the analytical solution can be derived by normal stress used as a criterion,improving the multiple cyclic method which determines the position of neutral point by the principal stress. Meanwhile, a comparison is made between the results of the analytical solution and those calculated from the classic mechanics theory, assuming the tension modulus is equal to the compression modulus, and those from the finite element method (FEM) numerical solution. The comparison shows that the analytical solution considers well the effects caused by the condition of different tension and compression modulus. Finally, a calculation correction of the structure with different modulus is proposed to optimize the structure.     

Fluid-solid coupling mathematical model of contaminant transport in unsaturated zone and its asymptotical solution

IntroductionThetransportofcontaminantsinunsaturatedzonehascausedmuchattention .Inearly1960s,contaminationproblemsofsoilandgroundwaterhadbeenstudiedathomeandabroad[1].Andinrecentyears ,thetransformationandtransportationofcontaminantshavebeendeeplystudiedinthefieldsofhydrogeology ,petroleumengineering ,environmentalengineeringandsoon[2 ,3].Somecontaminanttransportmodelshavebeenpresentedsofar.Forexample ,Paker[4 ]etal.presentedaconstitutivemodelgoverningparametersofwater,gasandcontaminantswhenth…     

Analytical solution for functionally graded anisotropic cantilever beam subjected to linearly distributed load

The bending problem of a functionally graded anisotropic cantilever beam subjected to a linearly distributed load is investigated.The analysis is based on the exact elasticity equations for the plane stress problem.The stress function is introduced and assumed in the form of a polynomial of the longitudinal coordinate.The expressions for stress components are then educed from the stress function by simple differentiation. The stress function is determined from the compatibility equation as well as the bound- ary conditions by a skilful deduction.The analytical solution is compared with FEM calculation,indicating a good agreement.     

Analytical solution of a simply supported piezoelectric beam subjected to a uniformly distributed loading

IntroductionIntelligentstructureisakindofnewstructuremodelsandhasbeenreceivedmuchattentioninrecentyears.Usedaspiezoelectricsensorsandactuators,thiskindofintelligentstructureshasmanyadvantagessuchasthepromptnessofresponseandtheconvenienceforthesignalcontrol.Soitisfollowedwithinterestbothinthetheoreticalstudyandintheengineeringapplications[1].Forexample ,theresearchoffundamentalsolution[2 ,3]andellipsoidalinclusion[4 ]forthree_dimensionalpiezoelectricmaterialbyuseofFouriertransformation ;thestu…     

STRUCTURAL DAMAGE MODEL OF UNSATURATED EXPANSIVE SOIL AND ITS APPLICATION IN MULTI-FIELD COUPLE ANALYSIS ON EXPANSIVE SOIL SLOPE

Focused on the sensitivity to climate change and the special mechanical characteristics of undisturbed expansive soil, an elastc-plastic damage constitutive model was proposed based on the mechanics of unsaturated soil and the mechanics of damage. Undisturbed expansive soil was considered as a compound of non-damaged part and damaged part. The behavior of the non-damaged part was described using non-linear constitutive model of unsaturated soil. The property of the damaged part was described using a damage evolution equation and two yield surfaces, i.e., loading yield (LY) and shear yield (SY). Furthermore, a consolidation model for unsaturated undisturbed expansive soil was established and a FEM program named UESEPDC was designed. Numerical analysis on solid-liquid-gas tri-phases and multi-field couple problem was conducted for four stages and fields of stress, displacement, pore water pressure, pore air pressure, water content, suction, and the damage region as well as plastic region in an expansive soil slope were obtained.     

Advances in modeling of water in the unsaturated zone

This paper reviews recent advances in analytical and numerical solution of problems of water flow through rigid soils in the unsaturated zone. The Richards model remains the most widely accepted and fertile framework for water flow analyses. More general formulations are reserved for the analysis of problems involving macroporosity, thermal effects, and air pressure effects. Many exact and approximate solutions have been derived for particular boundary value problems of homogeneous soils using methods such as quasi-linear analysis, Green-Ampt analysis, perturbation, and the kinematic wave approximation. Numerical simulators have become bigger and more accurate due to improvements in the areas of nonlinear solution procedures, mass conservation, computational efficiency, and computer hardware. Problems of natural heterogeneity have been addressed primarily through application of various stochastic methods to the Richards model. The stochastic formulations generally refute the concept of simple equivalent homogeneous properties, but do themselves offer a certain limited potential for a predictive capability.     

Application of dispersion theory to time domain reflectometry in soils

With time domain reflectometry (TDR) two dispersive parameters, the dielectric constant, , and the electrical conductivity, can be measured. Both parameters are nonlinear functions of the volume fractions in soil. Because the volume function of water ( _w) can change widely in the same soil, empirical equations have been derived to describe these relations. In this paper, a theoretical model is proposed based upon the theory of dispersive behaviour. This is compared with the empirical equations. The agreement between the empirical and theoretical aproaches was highly significant: the ( _w) relation of Topp et al. had a coefficient of determination r ² = 0.996 and the (_u) relation of Smith and Tice, for the unfrozen water content, _u, at temperatures below 0°C, had an r ² = 0.997. To obtain ( _w) relations, calibration measurements were performed on two soils: Caledon sand and Guelph silt loam. For both soils, an r ² = 0.983 was obtained between the theoretical model and the measured values. The correct relations are especially important at low water contents, where the interaction between water molecules and soil particles is strong.     

An Analytical Solution to the One-Dimensional Solute Advection-Dispersion Equation in Multi-Layer Porous Media

An analytical solution to the one-dimensional solute advection-dispersion equation in multi-layer porous media is derived using a generalized integral transform method. The solution was derived under conditions of steady-state flow and arbitrary initial and inlet boundary conditions. The results obtained by this solution agree well with the results obtained by numerically inverting Laplace transform-generated solutions previously published in the literature. The analytical solution presented in this paper provides more flexibility with regard to the inlet conditions. The numerical evaluation of eigenvalues and matrix exponentials required in this solution technique can be accurately and efficiently computed using the sign-count method and eigenvalue evaluation methods commonly available. The illustrative calculations presented herein have shown how an analytical solution can provide insight into contaminant distribution and breakthrough in transport through well defined layered column systems. We also note that the method described here is readily adaptable to two and three-dimensional transport problems.     

State space solution to 3D multilayered elastic soils based on order reduction method

Starting with the governing equations in terms of displacements of 3D elastic media, the solutions to displacement components and their first derivatives are obtained by the application of a double Fourier transform and an order reduction method based on the Cayley-Hamilton theorem. Combining the solutions and the constitutive equations which connect the displacements and stresses, the transfer matrix of a single soil layer is acquired. Then, the state space solution to multilayered elastic soils is further obtained by introducing the boundary conditions and continuity conditions between adjacent soil layers. The numerical analysis based on the present theory is carried out, and the vertical displacements of multilayered foundation with a weak and a hard underlying stratums are compared and discussed.     

Exact analytical solution to three-dimensional phase change heat transfer problems in biological tissues subject to freezing

Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.     

A General Formulation of Stress Distribution in Cylinders Subjected to Non-Uniform External Pressure

This work presents a two-dimensional stress analysis for elastic solid cylinders subjected to combined loading. The loading is generally formed with a number of concentrated and partially distributed forces all applied radially on the outer surface. The distributed forces cause pressures with non-uniform intensity along the circumferential direction. The cylinder is assumed to be long so that a state of plane-strain is valid. To obtain the stress distribution for the problem of partially distributed forces a new approach is followed first introduced in this paper. It is based on the expressions formed after using the theory of simple radial stress distribution when point-forces are applied on the cylinder and leads to the solution after direct integration. The total stresses due to both concentrated and distributed forces are obtained using the method of superposition. Apart from its simplified formulation, this general solution is always preferable since it proved to have a great advantage. As a result of not containing Fourier series, it eliminates some problems of convergence of the series at the boundaries that appear due to the Gibbs phenomena when the boundary conditions are a discontinuous function. Numerical results are presented for some interesting cases of loading conditions. This revised version was published online in August 2006 with corrections to the Cover Date.