首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 832 毫秒
1.
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.  相似文献   

2.
In this paper, a time-dependent infinite element which can be used to simulate transient seepage problems in infinite media is presented. The hydraulic head distribution function of the element has been derived in detail and the property matrices of the element have been well formulated. Since both space and time variables are used in the course of constructing the hydraulic head distribution function of the element, the present infinite element can be referred to as a transient one. Using the present infinite element to model the far field of a system, the mechanism of transient seepage problems in infinite media can be rigorously simulated because the property matrices of the element are evaluated at any time of interest in the analysis. Since explicit expressions can be written for the property matrices of the infinite element, they may be evaluated quite easily and this can be carried out by writing a simple subroutine in a computer program. In order to examine the accuracy and efficiency of the present infinite element, both a one-dimensional (ID) transient seepage problem in a semi-infinite medium and a 2D transient seepage problem in a full plane have been solved using the finite and infinite element technique. It has been demonstrated that the present infinite element is very useful for the numerical simulation of transient seepage problems in infinite media.  相似文献   

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

4.
A time‐domain viscous‐spring transmitting boundary is presented for transient dynamic analysis of saturated poroelastic media with linear elastic and isotropic properties. The u–U formulation of Biot equation in cylindrical coordinate is adopted in the derivation. By this general viscous‐spring boundary, the effective stress and pore fluid pressure on the truncated boundary of the computational area are replaced by a set of continuously distributed spring and dashpot elements, of which the parameters are defined assuming an infinite permeability and considering the two dilatational waves. Numerical examples demonstrate good absorption of both the two cylindrical dilatational waves by the proposed ‘drained’ boundary. For general two‐dimensional wave propagation problems, acceptable accuracy can still be achieved by setting the proposed boundary relatively far away from the scatter. Numerical comparison shows that the results obtained by using this boundary are more accurate for all permeability values than those by the traditional viscous‐spring or viscous boundaries established for u–U formulation. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

5.
A new finite element scheme is proposed, in this paper, for solving two-dimensional wave propagation problems in multilayered soils resting on a rigid base. The multilayered soils are treated as multiple horizontal layers of lateral infinite extension in geometry. Since these horizontal layers can be truncated by two artificially truncated vertical boundaries, two high-order artificial boundary conditions are applied for propagating the incoming waves from the interior domain into the far field of the system. Both the semi-analytical method and the truncated boundary migration procedure are used to derive the high-order artificial boundary conditions, which are comprised of a physically meaningful dashpot and a generalized energy absorber. The main advantage of using the proposed finite element scheme is that the derived artificial boundary condition can be straightforwardly implemented in the finite element analysis, without violating the band/sparse structure of the conventional finite element equation. The related numerical examples have demonstrated that the proposed finite element scheme is of high accuracy in dealing with wave propagation problems in multiple horizontal layers.  相似文献   

6.
Conventional modelling of transport problems for porous media usually assumes that the Darcy flow velocities are steady. In certain practical situations, the flow velocity can exhibit time‐dependency, either due to the transient character of the flow process or time dependency in the boundary conditions associated with potential flow. In this paper, we consider certain one‐ and three‐dimensional problems of the advective transport of a chemical species in a fluid‐saturated porous region. In particular, the advective flow velocity is governed by the piezo‐conduction equation that takes into account the compressibilities of the pore fluid and the porous skeleton. Time‐ and/or mesh‐refining adaptive schemes used in the computational modelling are developed on the basis of a Fourier analysis, which can lead to accurate and optimal solutions for the advective transport problem with time‐ and space‐dependent advective flow velocity distributions. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

7.
Based on the Biot theory, the exact solutions for one‐dimensional transient response of single layer of fluid‐saturated porous media and semi‐infinite media are developed, in which the fluid and solid particles are assumed to be compressible and the inertial, viscous and mechanical couplings are taken into account. First, the control equations in terms of the solid displacement u and a relative displacement w are expressed in matrix form. For problems of single layer under homogeneous boundary conditions, the eigen‐values and the eigen‐functions are obtained by means of the variable separation method, and the displacement vector u is put forward using the searching method. In the case of nonhomogeneous boundary conditions, the boundary conditions are first homogenized, and the displacement field is constructed basing upon the eigen‐functions. Making use of the orthogonality of eigen‐functions, a series of ordinary differential equations with respect to dimensionless time and their corresponding initial conditions are obtained. Those differential equations are solved by the state‐space method, and the series solutions for three typical nonhomogeneous boundary conditions are developed. For semi‐infinite media, the exact solutions in integral form for two kinds of nonhomogeneous boundary conditions are presented by applying the cosine and sine transforms to the basic equations. Finally, three examples are studied to illustrate the validity of the solutions, and to assess the influence of the dynamic permeability coefficient and the fluid inertia to the transient response of porous media. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

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

9.
A set of mapping functions in the form of convergent series for an infinite element, which is capable to include the infinitely distanced constant head boundary condition from the area of disturbance (e.g. pumping), is proposed based on the asymptotic far-field behaviour of typical seepage flow problems. The derived mapping functions have been successfully used in three-dimensional point symmetric, two-dimensional axi-symmetric and one-dimensional unidirectional flow for the fixed head boundary at infinite distance. The result shows excellent agreement with analytical solution. For the first time, the mapping function of an infinite element is presented in the form of a convergent series. The infinite elements are really capable of reducing the cost and efficiency of conventional finite element analysis. Finally, a figure is also proposed to indicate the required size of the near field to obtain accurate drawdown at specified locations based on some calculations for two-dimensional radial flow case.  相似文献   

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

11.
The diffusion equation governs thermal conduction and groundwater flow phenomena. In this paper, we study the two‐dimensional radial propagation of a sinusoidal perturbation radiating from a cylindrical source within an infinite slab of homogeneous material. The solution of this problem has several applications. For instance, it can be used to determine the hydraulic diffusivity of the subsurface based on measurements of the hydraulic head around a vertical well during its development. For thermal problems, it can be used to determine the thermal diffusivity based on measurements of the temperature distribution around a cylindrical heat source generating a sinusoidal power per unit length. In this paper, we present a comprehensive analytical solution of this problem and we compare these solutions with numerical solutions. Two approximate analytical solutions, which can be relevant in practice, are also presented. Finally, we give an upper bound for the survival time of the transient part of the solution and we provide an estimate of the radius of influence of the sinusoidal solicitation. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

12.
A superposition scheme is proposed to obtain a fundamental solution for boundary elements in multi‐layered elastic media. A three‐layered elastic region is obtained by superposing two sets of bonded half‐planes and subtracting one infinite plane. Therefore, the solution for an element in a layered media can be expressed in terms of bonded half‐plane solutions and an infinite‐plane solution. The major advantages of this superposition scheme are: (1) it is unnecessary to introduce elements at the interface, (2) it can be extended to higher‐order element, and (3) it may be applicable to three dimensions easily. The accuracy and performance of the developed model is illustrated by two examples. For the problem of a pressurized two‐dimensional crack within a three‐layered system, the comparison with other numerical results shows the model is quite accurate and efficient. The model is also used for a study of a practical two‐dimensional mining problem in South Africa, i.e. stoping through a dyke with material properties different from the host rock. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

13.
Natural or induced groundwater flow may negatively influence the performance of artificial ground freezing: high water flow velocities can prevent frozen conditions from developing. Reliable models that take into consideration hydraulic mechanisms are then needed to predict the ground freezing development. For forty years, numerous thermo-hydraulic coupled numerical models have been developed. Some of these models have been validated against experimental data but only one has been tested under high water flow velocity conditions. This paper describes a coupled thermo-hydraulic numerical model completely thermodynamically consistent and designed to simulate artificial ground freezing of a saturated and non-deformable porous medium under seepage flow conditions. On some points, less restrictive assumptions than the ones usually used in the literature are considered. As for the constant-porosity assumption, its validity is verified. The model appears to be well validated against analytical solutions and a three-dimensional ground freezing experiment under high seepage flow velocity conditions. It is used to highlight key thermo-hydraulic mechanisms associated with phase change in a porous medium.  相似文献   

14.
A numerical approach is proposed to model the flow in porous media using homogenization theory. The proposed concept involves the analyses of micro‐true flow at pore‐level and macro‐seepage flow at macro‐level. Macro‐seepage and microscopic characteristic flow equations are first derived from the Navier–Stokes equation at low Reynolds number through a two‐scale homogenization method. This homogenization method adopts an asymptotic expansion of velocity and pressure through the micro‐structures of porous media. A slightly compressible condition is introduced to express the characteristic flow through only characteristic velocity. This characteristic flow is then numerically solved using a penalty FEM scheme. Reduced integration technique is introduced for the volumetric term to avoid mesh locking. Finally, the numerical model is examined using two sets of permeability test data on clay and one set of permeability test data on sand. The numerical predictions agree well with the experimental data if constraint water film is considered for clay and two‐dimensional cross‐connection effect is included for sand. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

15.
The formulation of an axi-symmetric infinite element for transient analysis of flow problems in unbounded domain is presented. The theoretical basis as well as the implementation of the element is discussed, and the element decay function is derived using the analytical solution of a one-dimensional axially symmetric configuration. The form of decay within the element is described as a function of both time and space, and thus the hydraulic head distribution in the far field is simulated rigorously. The accuracy and the efficiency of the proposed element are demonstrated through several numerical examples in infinite media. In general, it is shown that using the present infinite element transient flow problems in unbounded domains can be simulated effectively. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

16.
This paper presents a novel analytical solution to the transient, z‐dependent, and asymmetric problem of an infinite wellbore drilled into a fluid‐saturated porous medium. The formulations are based on Biot's linear theory of poroelasticity, in which the dependency of poroelastic field variables to spatial coordinates as well as time domain is considered in the most general form. This gives flexibility to the solution in cases that cannot be analyzed using the conventional plane strain or symmetric models. One such case is when calculating the stress variations around an inclined wellbore where the far‐field stresses are acting over a finite vertical section. The results of our solution to this case with a three‐dimensional state of far‐field stress are used to analyze the stability of inclined wellbores passing through abnormally stressed formations. The presented solution is capable of finding expressions for fundamental solutions with stress or flow boundary conditions at the wellbore. These solutions are here adopted to analyze the pressure disturbances generated by multiprobe formation tester, a standard wireline device that is designed for downhole fluid sampling as well as estimating the directional permeabilities of subsurface earth formations. A comparison with the conventional solution for the relevant pressure diffusion equation indicates that the poroelastic effect is relatively significant in relation to the transient response of the pore pressure. Further, it is shown that the finite dimensions of sink probe would, to a great extent, contribute to the formation's pore pressure variations at its immediate proximity. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

17.
In modeling of many geomechanics problems such as underground openings, soil-foundation structure interaction problems, and in wave propagation problems through semi-infinite soil medium the soil is represented as a region of either infinite or semi-infinite extent. Numerical modeling of such problems using conventional finite elements involves a truncation of the far field in which the infinite boundary is terminated at a finite distance. In these problems, appropriate boundary conditions are introduced to approximate the solution of the infinite or semi-infinite boundaries as closely as possible. However, the task of positioning the finite boundary in conventional finite element discretization and the definition of the boundary and its conditions is very delicate and depends on the modeller's skill and intuition. Moreover, such a choice is influenced by the size of the domain to be discretized. Consequently, the dimensions of the global matrices and the time required for solution of the problem will increase considerably and also selection of the arbitrary location of truncated boundary may lead to erroneous result. In order to over come these problems, mapped infinite elements have been developed by earlier researchers (Simoni and Schrefier, 1987). In the present work the applicability of infinite element technique is examined for different geomechanics problems. A computer program INFEMEP is developed based on the conventional finite element and mapped infinite element technique. It is then validated using selected problems such as strip footing and circular footing. CPU time taken to obtain solutions using finite element approach and infinite element approach was estimated and presented to show the capability of coupled modeling in improving the computational efficiency. Mesh configurations of different sizes were used to explore the enhancement of both computational economy and solution accuracy achieved by incorporation of infinite elements to solve elastic and elasto-plastic problems in semi-infinite/finite domain as applied to geotechnical engineering. © Rapid Science Ltd. 1998  相似文献   

18.
This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

19.
In this study, a simplified analytical closed‐form solution, considering plane strain and axial symmetry conditions, for analysis of a circular pressure tunnel excavated underwater table, is developed. The method accounts for the seepage forces with the steady‐state flow and is based on the generalized effective stress law. To examine the effect of pore pressure variations and also the boundary conditions at the ground surface, the formulations are derived for different directions around the tunnel. The proposed method can be applied for analysis and design of pressure tunnels. Illustrative examples are given to demonstrate the performance of the proposed solution and also to examine the effect of seepage forces on the stability of tunnels. The simplified analytical solution derived in this study is compared with numerical analyses. It is concluded that the classic solutions (Lame's thick‐walled solution), considering the internal pressure as a mechanical load applied to the tunnel surface, are not applicable to pervious media and can result in an unsafe design. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

20.
This paper introduces an exact analytical solution for governing flow equations for one‐dimensional consolidation in unsaturated soil stratum using the techniques of eigenfunction expansion and Laplace transformation. The homogeneous boundary conditions adopted in this study are as follows: (i) a one‐way drainage system of homogenous soils, in which the top surface is considered as permeable to air and water, whereas the base is an impervious bedrock; and (ii) a two‐way drainage system where both soil ends allow free dissipation of pore‐air and pore‐water pressures. In addition, the analytical development adopts initial conditions capturing both uniform and linear distributions of the initial excess pore pressures within the soil stratum. Eigenfunctions and eigenvalues are parts of the general solution and can be obtained based on the proposed boundary conditions. Besides, the Laplace transform method is adopted to solve the first‐order differential equations. Once equations with transformed domain are all obtained, the final solutions, which are proposed to be functions of time and depth, can be achieved by taking an inverse Laplace transform. To verify the proposed solution, two worked examples are provided to present the consolidation characteristics of unsaturated soils based on the proposed method. The validation of the recent results against other existing analytical solutions is graphically demonstrated. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

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

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