共查询到20条相似文献，搜索用时 15 毫秒
1.
A new method to calculate the transport of dissolved species in aquifers is presented. This approach is an extension of the stream tubes which are used for flow computation. The flux tubes defined here are conservative for solutes, but not for water mass. The flux tubes are first defined in a general domain and then calculated in a two-dimensional uniform flow field. The tubes?? computation is based on a parametric solution. The method is extended further in order to deal with heterogeneous media. A particle-tracking algorithm is used where the deviation of the flux-tube boundaries due to dispersion is accounted for. The approximate solution obtained by this approach is compared to classical numerical solutions given by a finite difference code (RT3D) and a finite element code (FEFLOW). This comparison was performed for several test cases with increasing complexity. The differences between the flux-tube approach and the other methods always remain small, even regarding mass conservation. The major advantage of the flux-tube approach is the ability to reach a solution quickly, as the method is hundreds to thousands of times faster than classical finite difference or finite element models. 相似文献
2.
One dimensional solutions for the classic critical upward seepage gradient/quick condition and the time rate of consolidation problems are obtained using coupled routines for the finite volume method (FVM) and discrete element method (DEM), and the results compared with the analytical solutions. The two phase flow in a system composed of fluid and solid is simulated with the fluid phase modeled by solving the averaged Navier–Stokes equation using the FVM and the solid phase is modeled using the DEM. A framework is described for the coupling of two open source computer codes: YADE-OpenDEM for the discrete element method and OpenFOAM for the computational fluid dynamics. The particle–fluid interaction is quantified using a semi-empirical relationship proposed by Ergun [12]. The two classical verification problems are used to explore issues encountered when using coupled flow DEM codes, namely, the appropriate time step size for both the fluid and mechanical solution processes, the choice of the viscous damping coefficient, and the number of solid particles per finite fluid volume. 相似文献
3.
二维扩散输移问题的一种新的有限体积算法 总被引：1，自引：0，他引：1
从空间离散格式上比较了有限元格式与有限体积格式的异同处,证明有限体积格式是有限元格式的一种特例,有限元格式也是一种守恒型格式。并根据有限元格式与有限体积格式的比较,提出了一种新的高精度、稳定性好的有限体积格式。通过模拟非定常的纯对流方程和二维对流扩散方程,证明了有限体积格式的优点。 相似文献
4.
H. R. Thomas 《国际地质力学数值与分析法杂志》1988,12(1):31-44
The subject of this paper is the nonlinear analysis of heat and moisture transfer in partly saturated soil. In particular, an extension of the previous work of the author to two-dimensional applications is presented. Problems of greater complexity and engineering significance can hence be solved. The theoretical basis of the problem is first outlined before the two-dimensional formulation of the numerical solution is given. The finite element method is employed to give the spatial distributions of the variables, with a finite difference scheme being used to predict the temporal variations. The method allows the nonlinear nature of the soil parameters to be modelled, predictions being given of the coupled transient simultaneous transfer of heat and volumetric moisture content. The two-dimensional nature of the work is illustrated for the case of surface evaporative soil moisture losses from sloping ground. Subsoil conditions of a uniform deposit of loam are considered with values of representative material parameters assumed after a literature search. The coupled nature of heat and mass transfer is illustrated by means of an analysis of the moisture migration patterns that occur. Two-dimensional moisture flow, out of and into the soil, is shown to take place. 相似文献
5.
Coupled theory of mixtures for clayey soils 总被引：4，自引：0，他引：4
In this work, elasto-plastic coupled equations are formulated in order to describe the time-dependent deformation of saturated cohesive soils (two-phase state). Formulation of these equations is based on the principle of virtual work and the theory of mixtures for inelastic porous media. The theory of mixtures for a linear elastic porous skeleton was first developed by Biot (Theory of elasticity and consolidation for a porous anisotropic solid, Journal of Applied Physics, 1955, 26, 188–185). An extension of Biot's theory into a nonlinear inelastic media was performed by Prevost (Mechanics of continuous porous media, International Journal of Engineering Science, 1980, 18, 787–800). The saturated soil is considered as a mixture of two deformable media, the solid grains and the water. Each medium is regarded as a continuum and follows its own motion. The flow of pore-water through the voids is assumed to follow Darcy's law. The coupled equations are developed for large deformations with finite strains in an updated Lagrangian reference frame. The coupled behavior of the two-phase materials (soil-water state) is implemented in a finite element program. A modified Cam-clay model is adopted and implemented in the finite element program in order to describe the plastic behavior of clayey soils. Penetration of a piezocone penetrometer in soil is numerically simulated and implemented into a finite element program. The piezocone penetrometer is assumed to be infinitely stiff. The continuous penetration of the cone is simulated by applying an incremental vertical movement of the cone tip boundary. Results of the finite element numerical simulation are compared with experimental measurements conducted at Louisiana State University using the calibration chamber. The numerical simulation is carried out for two cases. In the first case, the interface friction between the soil and the piezocone penetrometer is neglected. In the second case, interface friction is assumed between the soil and the piezocone. The results of the numerical simulations are compared with experimental laboratory measurements. 相似文献
6.
7.
滑坡稳定性分析与安全系数取值研究 总被引：5，自引：0，他引：5
在地质灾害中,滑坡灾害分布最广、发生频率最高、危害最大,是我国地质灾害防治的主要对象。滑坡研究的主要任务之一是稳定性分析。论文针对三板溪水电站进水口滑坡和水布垭水电站大岩淌滑坡,分别采用刚体极限平衡法、块体单元法、有限单元法对稳定安全系数和条间推力分布曲线进行分析和比较,研究各种方法的差异和计算精度,以及各种方法稳定安全系数的取值标准,希望为水利水电工程边坡设计规范有关有限单元法、块体单元法条目的编写和相应的允许安全系数取值标准的制定提供一定的参考。研究结果表明:①块体单元法计算的安全系数与刚体极限平衡法计算的安全系数可以采用相同的取值标准,而有限单元法计算的安全系数取值标准可适当降低;②块体单元法计算的条间水平推力、条间剪切力与刚体极限平衡法(尤其是Sarma法)计算的条间水平推力、条间剪切力非常相近,有限单元法计算的条间水平推力、条间剪切力较小;③分析方法、力学模型与参数取值应与安全系数取值标准相匹配。 相似文献
8.
9.
By following a simple set of rules, a flow net can be manually constructed to obtain a graphical solution to the Laplace equation for simple two‐dimensional (2‐D) flow problems. With the advent of numerical solutions such as the finite difference and finite element methods, it is more common to generate a flow net automatically using the nodal head and flow values output by the computer program. Two methods have been published for automatically generating flow nets from finite element solutions: the stream‐function method and the particle‐tracking method. The stream‐function method works well for many cases, but it does not work for problems with holes in the mesh or internal sources or sinks. The particle‐tracking method works for all cases, but previously published algorithms that utilize this method do not result in the properly sized flow channels. A new approach is presented in this paper that is a hybrid of the stream‐function and particle‐tracking approaches. This method works for all cases and generates the properly sized flow channels. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
10.
11.
The paper presents a computational algorithm designed for efficient modelling of apparent resistivity over complex geological
structures, using finite element method. The algorithm can be used to study variations of apparent resistivities using any
electrode configuration at any point on the earth’s surface, not necessarily regular. A Schlumberger apparent resistivity
sounding curve over a buried anticline, is presented here as an example and compared with the corresponding analytical curve,
to demonstrate the correctness of the FEM algorithm.
The various potential derivatives required for the computation of apparent resistivities evaluated through different electrode
configurations have been obtained by calculating the ‘influence coefficients’ using reciprocal theorems, an approach successfully
applied in structural engineering. In essence, a set of self balancing nodal currents, obtained from the appropriate derivative(s)
of the shape functions of the elements contributing to the point of observation, is applied as the load vector.
The resulting quantities corresponding to the potential distribution in traditional finite element method, then, turn out
to be the potential derivatives at the point of observation for different positions of the current electrodes. These are known
as influence coefficients.
The continuum nature of the domain beyond the region of interest has been modelled by using ‘infinite elements’ across which
the potential is assumed to decay exponentially. 相似文献
12.
多尺度有限单元法求解非均质多孔介质中的三维地下水流问题 总被引：7，自引：0，他引：7
应用多尺度有限单元法模拟非均质多孔介质中的三维地下水流问题。与传统有限单元法相比，多尺度有限单元法的基函数具有能反映单元内参数变化的优点，所以这种方法能在大尺度上抓住解的小尺度特征获得较精确的解。在介绍多尺度有限单元法求解非均质多孔介质中三维地下水流问题的基本原理之后，对参数水平方向渐变垂直方向突变的非均质多孔介质中的三维地下水流和Borden实验场的三维地下水流分别用多尺度有限单元法和传统等参有限单元法进行了计算，结果表明在模拟高度非均质多孔介质中的三维地下水流问题时，多尺度有限单元法比传统有限单元法有效，既节省计算量又有较高的精度；在模拟非均质性弱的多孔介质中的三维地下水流问题时，多尺度有限单元法虽然也能在大尺度上获得较为精确的解，但效果不明显。 相似文献
13.
Savithru Jayasinghe David L. Darmofal Marshall C. Galbraith Nicholas K. Burgess Steven R. Allmaras 《Computational Geosciences》2018,22(2):527-542
This paper analyzes the adjoint equations and boundary conditions for porous media flow models, specifically the Buckley-Leverett equation, and the compressible two-phase flow equations in mass conservation form. An adjoint analysis of a general scalar hyperbolic conservation law whose primal solutions include a shock jump is initially presented, and the results are later specialized to the Buckley-Leverett equation. The non-convexity of the Buckley-Leverett flux function results in adjoint characteristics that are parallel to the shock front upstream of the shock and emerge from the shock front downstream of the shock. Thus, in contrast to the behavior of Burgers’ equation where the adjoint is continuous at a shock, the Buckley-Leverett adjoint, in general, contains a discontinuous jump across the shock. Discrete adjoint solutions from space-time discontinuous Galerkin finite element approximations of the Buckley-Leverett equation are shown to be consistent with the derived closed-form analytical solutions. Furthermore, a general result relating the adjoint equations for different (though equivalent) primal equations is used to relate the two-phase flow adjoints to the Buckley-Leverett adjoint. Adjoint solutions from space-time discontinuous Galerkin finite element approximations of the two-phase flow equations are observed to obey this relationship. 相似文献
14.
This paper adopts an upper bound procedure using the cell-based smoothed finite element method (CS-FEM) to estimate the seismic bearing capacity of shallow strip footings, focussing on seismic soil-structure interactions. In simulations, soil behaviour is assumed as the Mohr–Coulomb material, and increment of plasticity deformation obeys the associated flow rule. The first step of the numerical procedure involves approximating the kinematically admissible displacement fields using the cell-based smoothed finite element method, while the second relates to the establishments of the optimization problem as the conic programming. The inclusion of seismic conditions in the simulations was made using the pseudo-static approach. Initially, three seismic bearing capacity factors were resolved for both smooth and rough foundations by including horizontal and vertical inertia forces caused by the soil weight, the superstructure and the surcharge in the analyses. All seismic bearing capacity components obtained are in excellent agreement with those obtained using the method of characteristics and other finite element analyses. Subsequently, the reduction coefficients that correlate static and seismic bearing capacity factors were computed to facilitate the seismic design of the foundation.
相似文献15.
将地基视为多层各向同性弹性体,对Euler-Bernoulli梁进行有限单元离散分析,对地基-梁接触面采用边界积分法求解,根据地基-梁接触面的竖向位移协调和光滑接触条件,应用有限元与边界元耦合的方法推导出各向同性成层弹性地基上的Euler-Bernoulli梁的半数值半解析解。基于文中理论,编制了相应的程序,通过与现有文献对比,验证了文中理论的正确性,对比分析了分层地基与等效均匀地基两种模型上的弹性地基梁。研究结果表明,分层地基与等效均匀地基两种模型上弹性地基梁性状差异较大,实际弹性地基梁计算中应采用分层地基模型。 相似文献
16.
为解决城区渗流场中缝隙多、区域的大尺度与物体的小尺寸相差了几个数量级,有限元分析上存在的困难,根据缝隙流动的解析解,经适当地简化,将缝隙附近的渗流场分为缝外区和缝内区两部分,将缝外区流动简化为一个半径为半缝宽的井流,缝内区流动简化为一个均匀流。引入附加阻力系数,考虑缝隙出入口附近的由于过水断面突然变化而引的局部水头损失,并得到附加阻力系数的表达式。利用该式,结合渗流场中井点水头的修正公式,得到了能有效模拟渗流场中建筑物间缝隙的修正线单元公式。利用该线单元编制了相关程序模拟城区渗流场中的缝隙,使得缝隙出入口处的网格尺寸为缝隙两边的建筑物边长1/3～1/4时就能得到较为准确的结果,避免了区域性渗流场中小尺寸物体的网格划分问题,通过算例验证了该公式的精度和边界适应性。 相似文献
17.
The analysis of elasto-plastic boundary value problems using the finite element method involves many discretizations. These lead to the problem of yield surface drift in which the stress state predicted at the end of an elasto-plastic increment of loading does not lie on the current yield surface. As such discrepancies are comulative it is important to ensure that the stresses are corrected back to the yield surface during each increment of loading. In this paper five methods of accounting for this drift are examined. These involve correcting the stresses by projecting back along the plastic flow, the total strain increment and the accumulated effective stress direction. In addition a ‘correct’, method which accounts for the changes in elastic strains which accompany any stress correction is considered. This method is theoretically more sound than the other approximate approaches. All five methods have been used in finite element analyses of the stress changes that occur adjacent to a single pile installed in a uniform deposit of soil on pile loading. The soil was assumed to be normally consolidated and was modelled using a form of modified Cam Clay. Comparison of these results with an analysis, in which yield surface drift was negligible indicated that only the ‘correct’ method and the method involving projecting back along the plastic flow direction give accurate predictions. Substantiai errors occur if the other methods of correcting for yield surface drift are employed. It is recommended that the ‘correct’ method be adopted for finite element calculations. 相似文献
18.
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. 相似文献
19.
为解决城区渗流场中缝隙多、区域的大尺度与物体的小尺寸相差了几个数量级，有限元分析上存在的困难，根据缝隙流动的解析解，经适当地简化，将缝隙附近的渗流场分为缝外区和缝内区两部分，将缝外区流动简化为一个半径为半缝宽的井流，缝内区流动简化为一个均匀流。引入附加阻力系数，考虑缝隙出入口附近的由于过水断面突然变化而引的局部水头损失，并得到附加阻力系数的表达式。利用该式，结合渗流场中井点水头的修正公式，得到了能有效模拟渗流场中建筑物间缝隙的修正线单元公式。利用该线单元编制了相关程序模拟城区渗流场中的缝隙，使得缝隙出入口处的网格尺寸为缝隙两边的建筑物边长1/3～1/4时就能得到较为准确的结果，避免了区域性渗流场中小尺寸物体的网格划分问题，通过算例验证了该公式的精度和边界适应性。 相似文献
20.
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. 相似文献