共查询到20条相似文献，搜索用时 31 毫秒
1.
二维扩散输移问题的一种新的有限体积算法 总被引：1，自引：0，他引：1
从空间离散格式上比较了有限元格式与有限体积格式的异同处,证明有限体积格式是有限元格式的一种特例,有限元格式也是一种守恒型格式。并根据有限元格式与有限体积格式的比较,提出了一种新的高精度、稳定性好的有限体积格式。通过模拟非定常的纯对流方程和二维对流扩散方程,证明了有限体积格式的优点。 相似文献
2.
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. 相似文献
3.
锚索与锚杆联合锚固支护岩坡的有限元分析 总被引：1，自引：0，他引：1
针对广东某高速公路74 m高的7级挖方边坡,采用三维弹塑性有限元数值计算方法分析山区路堑开挖岩质边坡的锚固支护的受力情况。分别计算了坡体无锚和3种不同的锚固方案。通过坡体无锚方案的计算确定坡体的滑动面并为锚固方案的设计提供依据。3种锚固方案中的预应力锚索和锚杆按不同方式进行了组合。所建立的锚固坡体三维有限元分析模型对锚固坡体进行了合理简化并节约了计算量。计算中锚索的预应力采用预先施加应变的方式施加。计算分析表明,边坡加锚后可显著减少边坡塑性区。实际工程中采用预应力锚索和锚杆交错布设方案进行边坡锚固支护可较好地达到安全与经济并重的边坡处治原则。现场监测结果表明,锚固边坡有限元分析合理可行。 相似文献
4.
浅水湖泊风生流的迎风有限元数值模型研究 总被引：7，自引：1，他引：6
从控制方程组出发,以Galerkin有限元法为基础,引入沿流线加权的权函数(沿流线加权的迎风有限元法)和选择性集中系数矩阵方法,推导、建立了一适合于浅水湖泊风生流计算的二维迎风有限元数值模型。并以太湖为例,对模型作了检验,分析了均匀、定常风场持续作用下太湖风生流场的形态特征。 相似文献
5.
It is well known that the trial process for seeking the safety factor in the shear strength reduction finite element method (SSRFEM) is quite expensive, particularly for large 3-D slope stability analyses. The search algorithm for the safety factor is crucial to the entire solution process for the shear strength reduction finite element method, but few studies have attempted to exploit it. Among search algorithms, the commonly used bracketing and bisection search has not been fully optimised. Consequently, to improve the search scheme for the safety factor associated with the shear strength reduction finite element method, two strategies are suggested. First, a generalised bisection search algorithm is proposed to reduce the possibility of encountering non-convergence from a statistical point of view. To further improve the efficiency, a new two-grid scheme, characterised by a coarse mesh search and followed by a fine mesh search, is developed. Based on the drained or undrained analyses of the 3-D slope examples, the new search algorithm can markedly outperform the commonly used bisection search algorithms based on a single finite element mesh. 相似文献
6.
This paper presents a time integration scheme capable of simulating blast loading of relatively high frequency on porous media, using coarse meshes. The scheme is based on the partition of unity finite element method. The discontinuity is imposed on the velocity field, while the displacement field is kept continuous. The velocity discontinuity is postulated to occur in the time domain. The developed time integration scheme is unconditionally stable and has controllable numerical dissipation in the high frequency range. An important feature of the time scheme is that it allows for controlling the numerical damping in a consistent way. The time scheme has been implemented in combination with Biot’s theory of wave propagation in saturated porous media. Numerical examples have demonstrated that the proposed time scheme is, in addition to being accurate and stable, highly effective for coarse meshes. This makes the developed scheme suitable for large scale finite element analysis. 相似文献
7.
Slope stability has been a long-standing problem in geotechnical practice. The difficulty in analyzing slope stability may be attributed to geological complexity and versatile environments, but the limitations of analysis methods should not be ignored. In this work, the conventional limit equilibrium methods and the FE stress-based limit equilibrium methods with special emphasis on their limitations are reviewed, and the advantages of the finite element method with the capability of strength reduction are pointed out and illustrated by using some slope examples. Particularly, a flowchart of decision tree on selecting single-grid or two-grid search scheme is proposed for the shear strength reduction finite element method of slope stability analysis. Based on some slope examples, it is demonstrated that the automatic selection or switch between the single-grid scheme and the two-grid scheme can be remarkably superior to the single mesh scheme from the perspectives of computational cost and software development. 相似文献
8.
利用算子分裂迎风均衡格式解对流为主溶质运移问题 总被引：2，自引：0，他引：2
水污染模拟问题是水流问题与溶质运移问题的耦合问题.各种常见的数值解法在以对流为主溶质运移问题的求解中都会遇到困难,如用有限单元法或有限差分法时,会产生数值弥散与过量这两类误差.引入算子分裂迎风均衡格式法求解对流为主的水污染模拟问题,较好地克服了数值弥散和数值解出现振荡问题,该格式具有良好的稳定性、单调性及守恒性特点. 相似文献
9.
储层流固耦合的数学模型和非线性有限元方程 总被引：2，自引：0，他引：2
根据饱和多孔介质固体骨架的平衡方程和多孔介质中流体的连续性方程,建立了储层流固耦合数学模型。模型中引入了Jaumann应力速率公式描述多孔介质固体骨架的大变形效应,并考虑了地应力、初始孔隙压力、初始流体密度和初始孔隙度对耦合模型的影响。基于与微分方程等价的加权余量公式,在空间域采用有限元离散,对时间域进行隐式差分格式离散,导出了以单元节点位移和单元节点孔隙压力为未知量的储层流固耦合的非线性有限元增量方程。该模型在石油工程中有广泛的应用,为储层流固耦合的数值模拟奠定了理论基础。 相似文献
10.
One proposed geothermal heat extraction scheme relies on water circulation in a large vertical crack created by hydraulic fracturing in a hot dry impermeable rock mass. Water flow, heat convection and crack opening widths are analysed by finite elements. Governing field equations of the problem are first set up rigorously and then various small terms are identified and neglected, retaining the effects of pressure gradient, buoyancy, velocity head (kinetic energy) and head loss due to viscous friction in the water flow equation, and the effects of heat convection in water and heat conduction in rock in the heat transfer equation. The finite element scheme for water flow is based on a variational principle that is typical for diffusion problems, and for heat transfer it is based on the method of least-square residuals. The system of differential equations is highly non-linear. The non-linear terms and coefficients are treated in the fiaite element analysis as constant; the finite element analysisof, the steady-state pressures, fluxes and temperatures is then iterated, evaluating all non-linear terms and coefficients on the basis of the solution obtained in the previous iteration. Numerically calculated fields at various times after the start ofcooling are presented. They indicate some features favourable for the geothermal scheme, such as formation of eddy currents, and downward flux of water toward hotter rock. However, other important questions would have to be solved to gain full understanding, of this proposed geothermal scheme. 相似文献
11.
An improved search strategy for the critical slip surface using finite element stress fields 总被引：17，自引：0，他引：17
The finite element method can be used to advantage in slope stability problems. This paper proposes a technique to search for the critical slip surface as well as to define and calculate the factor of safety for the slope, when the finite element method is used to model its formation. First, stresses are estimated at each Gaussian point from the finite element analysis. Then, the global stress smoothing method is applied to get a continuous stress field. Based on this stress field, the factor of safety is calculated for a specified slip surface by a stress integration scheme. An improved search strategy is proposed for a noncircular critical surface which starts with a search method for a circular critical surface. During the search process, points defining a trial slip surface can freely move in the finite element mesh subject to some kinematical constraints. This method can be applied to both the limit equilibrium method and the finite element method. Effects of the slope stress history and soil parameters on the resulting critical surface are investigated. 相似文献
12.
13.
This paper is concerned with numerical methods for the modeling of flow and transport of contaminant in porous media. The
numerical methods feature the mixed finite element method over triangles as a solver to the Darcy flow equation and a conservative
finite volume scheme for the concentration equation. The convective term is approximated with a Godunov scheme over the dual
finite volume mesh, whereas the diffusion–dispersion term is discretized by piecewise linear conforming triangular finite
elements. It is shown that the scheme satisfies a discrete maximum principle. Numerical examples demonstrate the effectiveness
of the methodology for a coupled system that includes an elliptic equation and a diffusion–convection–reaction equation arising
when modeling flow and transport in heterogeneous porous media. The proposed scheme is robust, conservative, efficient, and
stable, as confirmed by numerical simulations.
相似文献
14.
The equations governing the dynamic behavior of saturated porous media as well as a finite element spatial discretization of these equations are summarized. A three-parameter time integration scheme called the Hilber–Hughes–Taylor α-method is used together with a predictor/multi-corrector algorithm, instead of the widely used Newmark's method, to integrate the spatially discrete finite element equations. The new time integration scheme possess quadratic accuracy and desirable numerical damping characteristics. The proposed numerical solution and bounding surface plasticity theory to describe the constitutive behaviour of soil have been implemented as the computer code DYSAC2. Predictions made by DYSAC2 code are verified using dynamic centrifuge test results for a clay embankment. Importance of initial state of a soil on its dynamic behaviour is demonstrated. 相似文献
15.
A space-time discontinuous Galerkin finite element method is proposed and applied to a convection-dominant single-phase flow
problem in porous media. The numerical scheme is based on a coupled space-time finite element discretization allowing for
discontinuous approximations in space and in time. The continuities on the element interfaces are weakly enforced by the flux
treatments, so that no extra penalty factor has to be determined. The resulting space-time formulation possesses the advantage
of capturing the steep concentration front with sharp gradients efficiently. The stability and reliability of the proposed
approach is demonstrated by numerical experiments.
The author is grateful to the DFG (German Science Foundation—Deutsche Forschungsgemeinschaft) for the financial support under
the grant number Di 430/4-2. 相似文献
16.
An unconditionally stable, fully explicit and highly precise multiple timescale finite element modeling scheme is described for a fully coupled hydro-mechanical (FCHM) analysis of saturated poroelastic media. The finite element method (FEM) is used for the discretization of the FCHM differential equation in the space domain. Direct integration is performed based on the precise time step integration method (PTSIM) for the time derivatives. Two configurations for the proposed scheme are constructed (abbreviated as PTSIM-f1 and -f2, respectively). The stability and convergence of the PTSIM-f1 and -f2 are proved using a matrix-based spectral analysis in the time domain. It is demonstrated that the explicit scheme proposed in this paper is unconditionally stable and independent of the time-step size. The algorithmic error estimation results indicate that the numerical modeling performed using PTSIM-f1 and -f2 in the time domain match the computer precision. Theoretically, the algorithmic error is caused by only the mesh discretization. Therefore, the proposed modeling scheme is a semi-analytical scheme. The applicability and accuracy of the proposed scheme are examined using sample calculations. By comparing with the analytical solutions, it is indicated that the modeling results have significant advantages over the standard FEM in terms of precision and computational efficiency for large timescales. 相似文献
17.
The yield vertex non-coaxial model is different from classical elastoplastic models, in that there is an additional plastic strain rate tangential to yield surfaces, as well as the plastic strain rate normal to yield surfaces, when orientations of principal stress change. This feature raises concerns on its finite element implementations. In nonlinear finite element numerical iterations, a large tangential plastic strain rate is likely to make the trial total strain rate direct inside a yield surface, which entails convergence difficulty. Some modifications are introduced on the non-coaxial model itself to make numerical convergence easier in the work published in Yang and Yu (2010) [20]. This paper is an extension of the previous work. Instead of modifying the non-coaxial model itself, this paper concerns the use of finite element explicit procedure, which is suitable for highly discontinuous problems. The simulations of shallow foundation load-settlement responses indicate that the finite element explicit procedure, assisted with a robust and explicit automatic substepping integration scheme of the non-coaxial model, does not encounter numerical difficulty. In addition, the overall trends of implicit and explicit simulations are similar. 相似文献
18.
In finite element calculations the constitutive model plays a key role. The evaluation of the stress response of the constitutive relation for a given strain increment, which is a time integration in the case of models of the rate type, is a typical sub task in such calculations. Adaptive behaviour of the time integration is essential to assure numerical stability and to control the accuracy of the solution. An adaptive second order semi-implicit method is developed in this paper. Its numerical behaviour is compared with an adaptive second order explicit scheme. The two proposed methods control the local error and guarantee numerical stability of the time integration. We include several numerical geotechnical element tests using hypoplasticity with intergranular strain. The element tests simulate the behaviour of a finite element method based on the displacement formulation. 相似文献
19.
George Zyvoloski 《国际地质力学数值与分析法杂志》1983,7(1):75-86
Two finite element algorithms suitable for long term simulation of geothermal reservoirs are presented. Both methods use a diagonal mass matrix and a Newton iteration scheme. The first scheme solves the 2N unsymmetric algebraic equations resulting from the finite element discretization of the equations governing the flow of heat and mass in porous media by using a banded equation solver. The second method, suitable for problems in which the transmissibility terms are small compared to the accumulation terms, reduces the set of N equations for the Newton corrections to a symmetric system. Comparison with finite difference schemes indicates that the proposed algorithms are competitive with existing methods. 相似文献
20.
Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method 
下载免费PDF全文
下载免费PDF全文 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. 相似文献