Unconfined seepage analysis in earth dams using smoothed fixed grid finite element method

One major difficulty in seepage analyses is finding the position of phreatic surface which is unknown at the beginning of solution and must be determined in an iterative process. The objective of the present study is to develop a novel non‐boundary‐fitted mesh finite‐element method capable of solving the unconfined seepage problem in domains with arbitrary geometry and continuously varied permeability. A new non‐boundary‐fitted finite element method named as smoothed fixed grid finite element method (SFGFEM) is used to simplify the solution of variable domain problem of unconfined seepage. The gradient smoothing technique, in which the area integrals are transformed into the line integrals around edges of smoothing cells, is used to obtain the element matrices. The solution process starts with an initial guess for the unknown boundary and SFGFEM is used to approximate the field variable. The boundary shape is then modified to eventually satisfy nonlinear boundary condition in an iterative process. Some numerical examples are solved to evaluate the applicability of the proposed method and the results are compared with those available in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Energy‐work‐based numerical manifold seepage analysis with an efficient scheme to locate the phreatic surface

A major challenge in seepage analysis is to locate the phreatic surface in an unconfined aquifer. The phreatic surface is unknown and assumed as a discontinuity separating the seepage domain into dry and wet parts, thus should be determined iteratively with special schemes. In this study, we systematically developed a new numerical manifold method (NMM) model for unconfined seepage analysis. The NMM is a general numerical method for modeling continuous and discontinuous deformation in a unified mathematical form. The novelty of our NMM model is rooted in the NMM two‐cover‐mesh system: the mathematical covers are fixed and the physical covers are adjusted with iterations to account for the discontinuity feature of the phreatic surface. We developed an energy‐work seepage model, which accommodates flexible approaches for boundary conditions and provides a form consistent with that in mechanical analysis with clarified physical meaning of the potential energy. In the framework of this energy‐work seepage model, we proposed a physical concept model (a pipe model) for constructing the penalty function used in the penalty method to uniformly deal with Dirichlet, Neumann, and material boundaries. The new NMM model was applied to study four example problems of unconfined seepage with varying geometric shape, boundary conditions, and material domains. The comparison of our simulation results to those of existing numerical models for these examples indicates that our NMM model can achieve a high accuracy and faster convergence speed with relatively coarse meshes. This NMM seepage model will be a key component of our future coupled hydro‐mechanical NMM model. Copyright © 2014 John Wiley & Sons, Ltd.     

A high‐frequency open boundary for transient seepage analyses of semi‐infinite layers by extending the scaled boundary finite element method

A high‐frequency open boundary has been developed for the transient seepage analyses of semi‐infinite layers with a constant depth. The scaled boundary finite element equation of pore water pressure is formulated first in the frequency domain. With the eigenvalue problem, the equation can be decoupled into modal equations whose modal dynamic permeability equation can be determined. The continued fraction technique is adopted to formulate the continued fraction solution in the frequency domain. All constants in the solution are determined recursively at the high‐frequency limit. By introducing auxiliary variables and the continued fraction solution to the relationship between the prescribed seepage flow and the pore water pressure in the frequency domain, the open boundary condition is obtained. After transformed to the time domain, the open boundary condition is expressed as a system of fractional differential equations. No convolution integral is required. The accuracy of the analysis results increases with the increasing orders of continued fraction. Copyright © 2015 John Wiley & Sons, Ltd.     

A local artificial boundary for transient seepage problems with unsteady boundary conditions in unbounded domains

Numerical analysis of transient seepage in unbounded domains with unsteady boundary conditions requires a more sophisticated artificial boundary approach to deal with the infinite character of the domain. To that end, a local artificial boundary is established by simplifying a global artificial boundary. The global artificial boundary conditions (ABCs) at the truncated boundary are derived from analytical solutions for one‐dimensional axisymmetric diffusion problems. By applying Laplace transforms and introducing some specially defined auxiliary variables, the global ABCs are simplified to local ABCs to significantly enhance the computational efficiency. The proposed local ABCs are implemented in a finite element computer program so that the solutions to various seepage problems can be calculated. The proposed approach is first verified by the computation of a one‐dimensional radial flow problem and then tentatively applied to more general two‐dimensional cylindrical problems and planar problems. The solutions obtained using the local ABCs are compared with those obtained using a large element mesh and using a previously proposed local boundary. This comparison demonstrates the satisfactory performance and obvious superiority of the newly established boundary to the other local boundary. Copyright © 2017 John Wiley & Sons, Ltd.     

A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows

The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Semi‐analytical elastostatic analysis of unbounded two‐dimensional domains

Unbounded plane stress and plane strain domains subjected to static loading undergo infinite displacements, even when the zero displacement boundary condition at infinity is enforced. However, the stress and strain fields are well behaved, and are of practical interest. This causes significant difficulty when analysis is attempted using displacement‐based numerical methods, such as the finite‐element method. To circumvent this difficulty problems of this nature are often changed subtly before analysis to limit the displacements to finite values. Such a process is unsatisfactory, as it distorts the solution in some way, and may lead to a stiffness matrix that is nearly singular. In this paper, the semi‐analytical scaled boundary finite‐element method is extended to permit the analysis of such problems without requiring any modification of the problem itself. This is possible because the governing differential equations are solved analytically in the radial direction. The displacement solutions so obtained include an infinite component, but relative motion between any two points in the unbounded domain can be computed accurately. No small arbitrary constants are introduced, no arbitrary truncation of the domain is performed, and no ill‐conditioned matrices are inverted. Copyright © 2002 John Wiley & Sons, Ltd.     

Semi‐analytical far field model for three‐dimensional finite‐element analysis

A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd.     

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

杂交边界点法求解带自由面的渗流问题

将杂交边界点法与迭代法相结合,求解有自由面的渗流问题。杂交边界点法基于杂交位移变分原理和移动最小二乘近似,利用基本解插值域内的场函数,而边界上的变量则用移动最小二乘近似,是一种纯边界类型的无网格方法。利用该方法只需在边界上布点而不需要划分任何网格的特性,先假定自由面的初始位置,再进行迭代求解。数值算例表明,该方法精度较高、计算量较小,适合于求解各种具有自由面的渗流问题。     

One‐dimensional analysis of a poroelastic medium during freezing

A solution to the problem of freezing of a poroelastic material is derived and analysed in the case of one‐dimensional deformation. The solution is sought within the framework of thermo‐poroelasticity, with specific account of the behaviour of freezing materials. The governing equations of the problem can be combined into a pair of coupled partial differential equations for the temperature and the fluid pressure, with particular forms in the freezing and the unfrozen regions. In the freezing region, the equations are highly non‐linear, partly due to the dependence of thermal and hydraulic properties on water saturation, which varies with temperature. Consequently, the solution is obtained through numerical methods, with special attention to the propagation of the freezing front boundary. The response to one‐dimensional freezing is illustrated for the case of cement paste. Finally, the influence on the solution of varying selected parameters is analysed, such as the temperature boundary conditions, the parameters characterizing the geometry of the porous system, the ratio of fluid and thermal diffusivities, and the rate of cooling applied at the freezing end. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamic responses of unsaturated half‐space soil to a moving harmonic rectangular load

The problem of the dynamic responses of a semi‐infinite unsaturated poroelastic medium subjected to a moving rectangular load is investigated analytical/numerically. The dynamic governing equations are obtained with consideration of the compressibility of solid grain and pore fluid, inertial coupling, and viscous drag as well as capillary pressure in the unsaturated soil, and they can be easily degraded to the complete Biot's theory. Using the Fourier transform, the general solution for the equations is derived in the transformed domain, and then a corresponding boundary value problem is formulated. By introducing fast Fourier transform algorithm, the unsaturated soil vertical displacements, effective stresses, and pore pressures induced by moving load are computed, and some of the calculated results are compared with those for the degenerated solution of saturated soils and confirmed. The influences of the saturation, the load speed, and excitation frequency on the response of the unsaturated half‐space soil are investigated. The numerical results reveal that the effects of these parameters on the dynamic response of the unsaturated soil are significant.     

Explanation of locking of four‐node plane element by considering it as elastic Dirichlet‐type boundary value problem

The aim of this study is to arrive at a better understanding of the phenomenon of locking of low‐order compatible displacement type of finite elements in particular for the hour‐glass mode of the plane four‐node element and dilative materials. To this end the properties of finite elements are investigated in an analytical way, where a finite element is considered as a plane boundary value problem with prescribed boundary displacement (Dirichlet problem). In this paper for the sake of simplicity the simplest possible linear comparison solid, namely isotropic linear elasticity, is applied, although recognizing fully that for a dilative material elasto‐plasticity would be more realistic. From the study described in this paper it is concluded that locking of the four‐node element is not due to any particular numerical formulation of this compatible finite element since, even the analytical solution suffers from this problem. The locking of this element is not related to incompressibility of the material either as the analytical solution shows locking to occur at a parameter set which differs significantly from the one in case of incompressibility. It is shown that locking is a consequence of the combination of the dilative material behaviour and the compatible displacement type of boundary conditions, which leads to infinite isotropic stresses in the element. These infinite isotropic stresses occur at the limit of uniqueness of the solution, which for this element is shown to occur outside the parameter range of the sufficiency of uniqueness. Copyright © 2000 John Wiley & Sons, Ltd.     

An efficient finite–discrete element method for quasi‐static nonlinear soil–structure interaction problems

An efficient finite–discrete element method applicable for the analysis of quasi‐static nonlinear soil–structure interaction problems involving large deformations in three‐dimensional space was presented in this paper. The present method differs from previous approaches in that the use of very fine mesh and small time steps was not needed to stabilize the calculation. The domain involving the large displacement was modeled using discrete elements, whereas the rest of the domain was modeled using finite elements. Forces acting on the discrete and finite elements were related by introducing interface elements at the boundary of the two domains. To improve the stability of the developed method, we used explicit time integration with different damping schemes applied to each domain to relax the system and to reach stability condition. With appropriate damping schemes, a relatively coarse finite element mesh can be used, resulting in significant savings in the computation time. The proposed algorithm was validated using three different benchmark problems, and the numerical results were compared with existing analytical and numerical solutions. The algorithm performance in solving practical soil–structure interaction problems was also investigated by simulating a large‐scale soft ground tunneling problem involving soil loss near an existing lining. Copyright © 2011 John Wiley & Sons, Ltd.     

One‐dimensional consolidation with a threshold gradient: a Stefan problem with rate‐dependent latent heat

During the process of one‐dimensional consolidation with a threshold gradient, the seepage front moves downward gradually, and the problem is indicated as a Stefan problem. The novel feature in this Stefan problem is a latent heat that varies inversely with the rate of the moving boundary. An exact solution for the external load that increases in proportion to the square root of time is constructed using the similarity transformation technique. Computational examples concerning the effect of different parameters on the motion of the seepage front are presented. The exact solution provides a worthwhile benchmark for verifying the accuracy of numerical and approximate methods. Copyright © 2013 John Wiley & Sons, Ltd.     

昆仑山隧道冻融特征分析

从分析昆仑山隧道渗漏水地形特征入手,得出沟谷地形是多年冻土隧道渗漏水的原因,重新定义移动边界概念,提出产生移动边界特征的充要条件,为进一步认识多年冻土提供了理论依据。引入移动边界计算模型,分析了计算参数的取值原则,采用有限元技术,利用该模型计算和预测冻土的融化深度,并将计算结果与工程实际进行比较,最后得出昆仑山隧道冻融区的变化规律,为隧道病害治理提供了可靠途径。     

A moving boundary solution for solidification of lava lake and magma intrusion in the presence of time-varying contact temperature

During the solidification of a lava lake heat is released convectively from the top surface as well as conductively into the country rock from the base, leading to non-uniform solidification. The upper solidified layer grows at a faster rate than the lower solidified layer. Similarly, solidification of magma intrusion within the crust is also non-uniform due to the presence of thermal gradient in the crust. Available analytical solution for solidification of a melt layer assumes only symmetric cooling about the centre of the layer. In the present work a moving boundary solution for thermal evolution and non-uniform solidification of a melt layer incorporating time-varying contact temperature conditions at both of its boundaries is developed. The solution is obtained by using the Fourier spectral approach in the space domain and a modified finite difference scheme in the time domain, and is validated with available analytical solutions for simple cases and a semi-analytical solution for the case involving temperature gradient in the country rock. This solution can be used to analyse solidification of lava lakes and magma intrusions experiencing time-dependent temperature variation at their contacts with the country rock.     

Theoretical and numerical study of the steady‐state flow through finite fractured porous media

This paper investigates the two‐dimensional flow problem through an anisotropic porous medium containing several intersecting curved fractures. First, the governing equations of steady‐state fluid flow in a fractured porous body are summarized. The flow follows Darcy's law in matrix and Poiseuille's law in fractures. An infinite transversal permeability is considered for the fractures. A multi‐region boundary element method is used to derive a general pressure solution as a function of discharge through the fractures and the pressure and the normal flux on the domain boundary. The obtained solution fully accounts for the interaction and the intersection between fractures. A numerical procedure based on collocation method is presented to compute the unknowns on the boundaries and on the fractures. The numerical solution is validated by comparing with finite element solution or the results obtained for an infinite matrix. Pressure fields in the matrix are illustrated for domains containing several interconnected fractures, and mass balance at the intersection points is also checked. Copyright © 2013 John Wiley & Sons, Ltd.     

Fractal approach applied to calculate hour‐glass eigenvalue of elastic square 4‐node quad

It is shown that the property of the scale invariance of the eigenvalues and eigenmodes of a finite element can be used as a basis to calculate good approximations to the analytical magnitudes of eigenvalues. This requires the subdivision of the element into a mesh of small elements with the same shape as the large element, the enforcement of the modal boundary displacements of the large element to the mesh of small elements and finally the application of the conditions of both the nodal equilibrium and the equality of the nodal work at both scales. Due to the self‐similarity of the elements at all scales the authors propose to call this method the fractal approach. The method is applied to calculate the hour‐glass eigenvalue of a plane square 4‐node quad for isotropic linear elastic material. The resulting hour‐glass eigenvalue is shown to be a good approximation of the analytical magnitude as derived in a companion paper. Copyright © 2000 John Wiley & Sons, Ltd.     

An operator‐split ALE model for large deformation analysis of geomaterials

Analysis of large deformation of geomaterials subjected to time‐varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator‐split arbitrary Lagrangian–Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid–fluid coupling and strong material non‐linearity. Each time step of the operator‐split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one‐dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd.     

An artificial boundary approach for unbounded transient seepage problems

A new artificial boundary approach for transient seepage problems in unbounded domain is presented. The artificial boundary condition at the truncated boundary is derived from the analytical solutions for transient seepage problems in one dimension, including solutions, respectively, for flow in one‐dimensional infinite space and for radial flow in an infinite layer, and then it is tentatively applied for some two dimensional problems in addition to the one‐dimensional problems mentioned above. The boundary conditions derived relate the time‐dependent boundary flux with the time derivative of the hydraulic head at the truncated boundary, which makes the implementation much easier compared with the infinite element method. The accuracy and efficiency of the artificial boundary are validated by several numerical examples, which shows that the proposed boundary can give very good results for one‐dimensional transient seepage problems, as expected, whereas reasonable results can be also obtained for two‐dimensional problems, such as two‐dimensional axisymmetric flow and flow in an infinite plane. Copyright © 2014 John Wiley & Sons, Ltd.