共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
Four classical geomechanics problems involving semi-infinite linear elastic media have been solved numerically using recently developed mapped infinite elements coupled to finite elements.The effect of the remoteness of the truncated boundary and the location of infinite element coupling on solution accuracy has been studied. The results of conventional analyses using finite elements over a relatively large but restricted region are compared to the coupled analyses. Comparison of the results shows that for the same number of degrees of freedom the performance of the coupled solutions is superior to the conventional approach with respect to accuracy of solution and computational efficiency. Finally, some general guidelines are proposed for the efficient numerical solution of these types of problems using the coupled finite/infinite element approach. 相似文献
3.
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 相似文献
4.
Using the numerical model presented in the first paper of this research,1 a parametric study has been carried out in this paper to investigate the effect of several important parameters on the transient contaminant transport in infinite porous fractured media. From the related numerical results, it has been demonstrated that: (1) transmissive coefficient between the porous block and the fissured network has a significant influence on the value of the concentration but has little effect on the speed of contaminant transport; (2) porosities in the porous block and fissured network have a significant influence on the maximum value of the concentration; (3) average linear velocity of flow has a significant influence on both the concentration distribution and speed of contaminant transport; (4) dispersion coefficient of the medium affects not only the shape of the concentration versus time curve but also the peak value of the concentration. 相似文献
5.
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. 相似文献
6.
Practical civil engineering problems are usually formulated in an infinite half-space domain, and a selected finite domain is required to analyze the dynamic responses of a fluid-saturated porous medium by the finite element method (FEM). Devising a method to deal with the boundaries of the finite domain is the key issue for this open system. In this paper, a two-dimensional spring–dashpot artificial boundary (SDAB) for transient analysis in a fluid-saturated porous media is developed. Based on Biot’s dynamic theory of fluid-saturated porous media, the normal and tangential boundary stress formulae are deduced for out-going cylindrical body waves. The boundary stress is proportional to displacement and velocity, thus continuously distributed dashpots and springs can be placed on the artificial boundaries in the normal and tangential directions to simulate the energy absorption of the infinite media outside of the finite domain for the interior distributed source problems. In this paper, the input seismic motion can be realized by applying an equivalent load on the SDAB for the seismic scattering problems of exterior distributed sources. Numerical examples are given and the analyzed results show that the SDAB and the method of wave motion input have good stability and acceptable accuracy. 相似文献
7.
Abbas El‐Zein 《国际地质力学数值与分析法杂志》2008,32(3):265-287
Finite‐element models of contaminant transport through composite landfill liners require highly refined meshes around the interface between the geomembrane and the clay layer, especially if leakage through holes in the geomembrane is considered. In addition, no general formulation for transport through leaking geomembranes can be found in the literature. The paper develops a general approach to time‐dependent contaminant migration through composite liners with intact or leaking geomembranes. Equations are derived for various combinations of system conditions including Dirichlet and Neumann boundary conditions in the waste, constant mass of contaminants in the waste, steady state or transient transport in the geomembrane, and steady state or transient seepage velocities in the mineral liner. The effect of the geomembrane on transport in the soil is converted into an equivalent boundary condition applicable at the top of the clay layer. Hence, only the media underlying the top geomembrane are explicitly represented in the numerical model, yielding a computationally efficient algorithm. The new formulation is validated in conjunction with finite‐layer, finite‐element and boundary‐element methods, by comparing its predictions to those of more conventional approaches which represent the geomembrane explicitly. The scope of the method is illustrated by modelling a landfill liner with a geomembrane leaking in five locations. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
8.
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. 相似文献
9.
An infinite element is presented to treat wave propagation problems in unbounded saturated porous media. The porous media is modeled by Biot's theory. Conventional finite elements are used to model the near field, whereas infinite elements are used to represent the behavior of the far field. They are constructed in such a way that the Sommerfeld radiation condition is fulfilled, i.e. the waves decay with distance and are not reflected at infinity. To provide the wave information the infinite elements are formulated in Laplace domain. The time domain solution is obtained by using the convolution quadrature method as the inverse Laplace transformation. The temporal behavior of the near field is calculated using standard time integration schemes, e.g. the Newmark method. Finally, the near and far field are combined using a substructure technique for any time step. The accuracy as well as the necessity of the proposed infinite elements, when unbounded domains are considered, is demonstrated by different examples. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
10.
This paper deals with the computational aspects of nonaqueous phase liquid (NAPL) dissolution front instability in two-dimensional
fluid-saturated porous media of finite domains. After the governing equations of an NAPL dissolution system are briefly described,
a combination of the finite element and finite difference methods is proposed to solve these equations. In the proposed numerical
procedure, the finite difference method is used to discretize time, while the finite element method is used to discretize
space. Two benchmark problems, for which either analytical results or previous solutions are available, are used to verify
the proposed numerical procedure. The related simulation results from these two benchmark problems have demonstrated that
the proposed numerical procedure is useful and applicable for simulating the morphological evolution of NAPL dissolution fronts
in two-dimensional fluid-saturated porous media of finite domains. As an application, the proposed numerical procedure has
been used to simulate morphological evolution processes for three kinds of NAPL dissolution fronts in supercritical NAPL dissolution
systems. It has been recognized that: (1) if the Zhao number of an NAPL dissolution system is in the lower range of the supercritical
Zhao numbers, the fundamental mode is predominant; (2) if the Zhao number is in the middle range of the supercritical Zhao
numbers, the (normal) fingering mode is the predominant pattern of the NAPL dissolution front; and (3) if the Zhao number
is in the higher range of the supercritical Zhao numbers, the fractal mode is predominant for the NAPL dissolution front. 相似文献
11.
Modelling of contaminant transport through landfill liners and natural soil deposits is an important area of research activity in geoenvironmental engineering. Conventional mesh‐based numerical methods depend on mesh/grid size and element connectivity and possess some difficulties when dealing with advection‐dominant transport problems. In the present investigation, an attempt has been made to provide a simple but sufficiently accurate methodology for numerical simulation of the two‐dimensional contaminant transport through the saturated homogeneous porous media and landfill liners using element‐free Galerkin method (EFGM). In the EFGM, an approximate solution is constructed entirely in terms of a set of nodes and no characterization of the interrelationship of the nodes is needed. The EFGM employs moving least‐square approximants to approximate the function and uses the Lagrange multiplier method for imposing essential boundary conditions. The results of the EFGM are validated using experimental results. Analytical and finite element solutions are also used to compare the results of the EFGM. In order to test the practical applicability and performance of the EFGM, three case studies of contaminant transport through the landfill liners are presented. A good agreement is obtained between the results of the EFGM and the field investigation data. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
12.
A number of phenomena and processes in geosciences can be summarized by second order partial differential equations. The major numerical methods for their solution include the classical finite difference method and the finite element method newly developed in the last two or three decades. Since 1977 the author has proved that for the Laplace and Poisson equations, these two methods are identical and are different only in the process of formulation. For transient problems, such as heat conduction in the earth and the groundwater and oil-gas unsteady flow in porous media, there are some differences in resulting linear algebraic euqations. In general, two methods give similar results, but when the time step is decreased to some extent, the resulting algebraic equation will be consistent with the anti-heat conduction equation rather than the original heat conduction equation. This is the reason why unrealistic potentials are produced by the finite element method. Such a problem can be overcome by using the 相似文献
13.
无限元是一种有效的人工边界,可用于处理弹性波的传播问题。在传统动力无限元的基础上,提出了一种采用分向插值技术的新型动力无限元,详细地推导了这种无限元的形函数,建立了完全解析形式的刚度矩阵,以提高计算效率,采用该无限元边界,计算了弹性介质中的线源Lamb问题,通过对比解析解答的地基表面位移,验证了该无限元的有效性。算例分析表明,采用此类无限元时,有限元单元边长建议取不超过1/8剪切波波长,网格边界到激励源点的距离宜取5倍剪切波波长。无限单元中的幅值衰减系数对计算结果影响甚微,建议取较小值。 相似文献
14.
D. W. Smith 《国际地质力学数值与分析法杂志》2000,24(8):693-722
Contaminant migration through soil is usually modelled mathematically using the dispersion–advection equation. This type of model finds application when planning the remediation of contaminated land, predicting the movement of polluted groundwater and designing engineered landfills. Usually the analysis assumes that the porous media through which the contaminant migrates is stationary. However, the construction of landfills on clay soils means that the soil beneath the landfill will undergo time‐dependent deformation as the soil consolidates. To date, there are no published data on the effect a deforming porous media may have on contaminant transport beneath a landfill; indeed, there appears to be no theory of contaminant migration through a deforming soil. In this paper, a one‐dimensional theory of contaminant migration through a saturated deforming porous media is developed based on a small and large strain analysis of a consolidating soil and conservation of contaminant mass. By selection of suitable parameters, the new transport equation reduces to the familiar one‐dimensional dispersion–advection equation for a saturated soil with linear, reversible, equilibrium controlled sorption of the contaminant onto the soil skeleton. Analytic solutions to a quasi‐steady‐state contaminant transport problem for a deforming media are presented, and a preliminary assessment made of the potential importance of soil deformation on the results of a contaminant migration analysis. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
15.
16.
A. P. S. Selvadurai 《国际地质力学数值与分析法杂志》2004,28(2):191-208
This paper examines the problem of the advective transport of a contaminant from sources in the shape of either a penny-shaped crack or an elongated needle-shaped cavity located in a porous medium of infinite extent. The advective transport is induced by Darcy flow in the porous medium, where the internal boundary is maintained at a constant potential. The paper presents an approximate analytical solution to this problem, which is deduced from a formulation that models a cavity in the shape of either an oblate or a prolate spheroid. The results also represent one of the few spatially three-dimensional exact analytical solutions for the, albeit linear, hyperbolic problem governing the contaminant transport problem. The paper also presents a canonical proof of uniqueness for advective contaminant transport problems associated with media of infinite extent. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
17.
A finite element method for modeling coupled flow and deformation in porous fractured media 
下载免费PDF全文
下载免费PDF全文 Ahmad Pouya 《国际地质力学数值与分析法杂志》2015,39(16):1836-1852
Modeling the flow in highly fractured porous media by finite element method (FEM) has met two difficulties: mesh generation for fractured domains and a rigorous formulation of the flow problem accounting for fracture/matrix, fracture/fracture, and fracture/boundary fluid mass exchanges. Based on the recent theoretical progress for mass balance conditions in multifractured porous bodies, the governing equations for coupled flow and deformation in these bodies are first established in this paper. A weak formulation for this problem is then established allowing to build a FEM. Taking benefit from recent development of mesh‐generating tools for fractured media, this weak formulation has been implemented in a numerical code and applied to some typical problems of hydromechanical coupling in fractured porous media. It is shown that in this way, the FEM that has proved its efficiency to model hydromechanical phenomena in porous media is extended with all its performances (calculation time, couplings, and nonlinearities) to fractured porous media. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
18.
利用算子分裂迎风均衡格式解对流为主溶质运移问题 总被引：2，自引：0，他引：2
水污染模拟问题是水流问题与溶质运移问题的耦合问题.各种常见的数值解法在以对流为主溶质运移问题的求解中都会遇到困难,如用有限单元法或有限差分法时,会产生数值弥散与过量这两类误差.引入算子分裂迎风均衡格式法求解对流为主的水污染模拟问题,较好地克服了数值弥散和数值解出现振荡问题,该格式具有良好的稳定性、单调性及守恒性特点. 相似文献
19.
非饱和土热湿耦合传输问题的数值解——理论及一维问题的解 总被引：1，自引：0，他引：1
本文介绍了求解非饱和土壤中热量和水分耦合传输问题的一种数值方法——积分有限差分方法(IFDM)。基于菲利普-迪弗瑞斯(Philip-De Vries)多孔介质热湿耦合流动模型,采用积分有限差分方法,编制了求解热湿传输问题的计算程序 HM1,该程序可用于求解各类工程中遇到的多孔介质一维传热传湿问题。文中还给出了计算结果与实测结果相比较的实例。 相似文献
20.
ABAQUS动力无限元人工边界研究 总被引：2，自引：0，他引：2
针对动力场天然无限地基的数值模拟与地震波输入问题进行了一些有意义的研究,评述了现有动力计算常用无限元的优缺点,详细阐述了ABAQUS无限元理论体系框架,并加以改进,提出一种考虑外域地震动影响的ABAQUS动力无限元人工边界。采用等效边界力的叠加原理,对入射波和散射波分开处理,视入射波和散射波在边界上互不影响,将输入地震动转化为作用于有限元无限元交界面上的等效应力的方法来解决外源波的入射问题。算例验证结果表明：内源振动和固定边界会出现失真和扰动现象,同时该计算结果与黏弹性边界的计算结果对比可知,该方法对外行散射波的过滤作用优于黏弹性边界。因此,改进的ABAQUS动力无限元人工边界理论方法有效且具有一定的稳定性。 相似文献