共查询到20条相似文献，搜索用时 21 毫秒
1.
2.
岩体裂隙网络非稳定渗流分析与数值模拟 总被引：1，自引：0，他引：1
针对裂隙岩体的非稳定渗流问题,通过将Darcy定理扩展到包含干区的整个裂隙网络区域,并令潜在溢出边界条件为Signorini型互补边界条件,将湿区上的非稳定渗流问题转化为全域上的一个新的初边值问题。为降低试探函数选取的难度,建立与定义在整个裂隙网络区域上的偏微分方程（PDE）提法等价的抛物型变分不等式（PVI）提法,并给出裂隙网络非稳定渗流分析的有限元数值分析格式和迭代算法,与砂槽模型试验数据的对比分析,验证其有效性。最后,将文中发展的计算方法应用到含复杂裂隙网络的边坡非稳定渗流分析,计算结果很好地反映出边坡内部自由面随库水降落的变化规律,并能准确地描述裂隙网络内部渗流运动特征及流量分布的不均匀性。 相似文献
3.
A finite element procedure is developed to accurately locate the free surface of unconfined seepage flow through porous media. The free surface is taken as the boundary between wet and dry soils, with flow in the saturated region characterized by Darcy's law. The method involves equations and meshing which are fully consistent with a general formulation for geotechnical engineering problems involving simultaneous solution of pore fluid pressures and soil skeleton displacements. Accuracy and versatility of the proposed procedure are demonstrated by solving various unconfined seepage flow problems through earth structures. Free surfaces and flownets are presented for the calculated flow fields. 相似文献
4.
Farhang Daneshmand 《国际地质力学数值与分析法杂志》2012,36(6):780-797
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. 相似文献
5.
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. 相似文献
6.
针对无压渗流问题,在传导矩阵调整法的基础上,对穿越自由面的单元用复合单元来处理,每个复合单元具有两套结点水头,一套用于饱和区的水头插值;另一套则用于无水区。通过变分原理,推出了控制方程并在程序中实现。应用复合单元法可以较好地解决自由面的插值问题,消除了穿过自由面单元的水上部分对自由面附近结点水头的贡献,同时对逸出面边界条件的处理也更加方便和准确。算例表明了该方法的有效性。 相似文献
7.
Investigation of a hybrid polygonal finite element formulation for confined and unconfined seepage 
下载免费PDF全文
下载免费PDF全文 This study presents a formulation for field problems using hybrid polygonal finite elements, taking steady state seepage through a porous material as the focus. We make comparisons with a conventional finite element formulation based on a single primary variable, focussing on the advantages of the hybrid formulation in terms of flux field accuracy and extension to convex polygonal shaped elements. For the unconfined case, we adopt a head dependent hydraulic conductivity that does not require remeshing. The performance of the hybrid polygonal element formulation is demonstrated through a series of numerical examples. The results show a sensitivity of the location of the free surface in unconfined seepage to mesh configuration for hybrid quadrilateral meshes with various aspect ratios, but not for hybrid polygonal meshes with various orientations and irregularity. Examination of the free surface location results for several conforming shape function options shows an insensitivity to choice of interpolation function, provided that it conforms with the assumptions in the formulation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
8.
A finite element, variable mesh analysis of unconfined steady-state seepage problems is presented based on a nonlinear programming algorithm. It is shown that the minimization of an objective function which merely represents a measure of the total flux leaving or entering the mesh at the free surface nodes (except those that belong also to pervious boundaries) does not permit a unique definition of the free surface geometry. This problem, which is apparently related to the numerical instabilities often met when using variable mesh approaches, can be eliminated by adding to the objective function a term representing a sort of overall ‘regularity’ condition for the shape of the free surface. The modified solution procedure turns out to be stable and able to provide meaningful results for practical problems even when rather coarse meshes are adopted. 相似文献
9.
Grid adaptive methods combined with domain adaptation are discussed for two-dimensional seepage flow problems with free boundaries through porous media. Examples of grid and domain adaptive methods are presented to demonstrate several ways to predict grids and shapes of free boundaries using an iterative scheme. Finally, the combined adaptive methods are applied to obtain smooth non-oscillatory shape of a free boundary of seepage flow through non-homogeneous porous media. 相似文献
10.
11.
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. 相似文献
12.
岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。 相似文献
13.
B‐bar based algorithm applied to meshfree numerical schemes to solve unconfined seepage problems through porous media 下载免费PDF全文
下载免费PDF全文 The object of this work is to establish a meshfree framework for solving coupled, steady and transient problems for unconfined seepage through porous media. The Biot's equations are formulated in displacements (or u − w) assuming an elastic solid skeleton. The free surface location and its evolution in time are obtained by interpolation of pore water pressures throughout the domain. Shape functions based on the principle of local maximum entropy are chosen for the meshfree approximation schemes. In order to avoid the locking involved in the fluid phase of the porous media, a B‐bar based algorithm is devised to compute the average volumetric strain in a patch composed of various integration points. The efficiency of such an implementation for one phase problems is shown through the Benchmark problem, Cook's membrane loaded by a distributive shear load. The proposed methodology is firstly applied to various classical examples in unconfined steady seepage problems through earth dams, then to the dynamic consolidation of a soil column. The results obtained for both problems are quite satisfactory and demonstrate the feasibility of the proposed method in solving coupled problems in porous media. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
14.
Noboru Kikuchi 《国际地质力学数值与分析法杂志》1977,1(3):283-297
The theory of variational inequalities enables us to formulate and solve free boundary problems in fixed domains, while most other methods assume the position of the unknown domain in solving the problem. Here the problem of seepage flow through a rectangular dam with a free boundary is formulated as a vertical inequality following the ideas of Baiocchi. In order to demonstrate the essential ideas of extending the domain of the solution of problems with free boundaries, the problem of the deflection, of a string on a rigid support is first examined. Next, variational inequalities are derived which are associated with several cases of seepage problems. An approximation theory, including a priori error estimates, is developed using finite element methods, and an associated numerical scheme is given. It is shown that for linear and quadratic finite element methods, the rates of convergence are 0(h) and 0(h1.25-δ), 0 < δ < 0.25, respectively, if the permeability is constant. 相似文献
15.
Stochastic analysis of free surface flow through earth dams 总被引：1，自引：0，他引：1
This research investigated the unconfined flow through dams. The hydraulic conductivity was modeled as spatially random field following lognormal distribution. Results showed that the seepage flow produced from the stochastic solution was smaller than its deterministic value. In addition, the free surface was observed to exit at a point lower than that obtained from the deterministic solution. When the hydraulic conductivity was strongly correlated in the horizontal direction than the vertical direction, the flow through the dam has markedly increased. It is suggested that it may not be necessary to construct a core in dams made from soils that exhibit high degree of variability. 相似文献
16.
By extending Darcy's law to the dry domain above the free surface and specifying the boundary condition on the potential seepage surfaces as Signorini's type, a partial differential equation (PDE) defined in the entire domain of interest is formulated for non‐steady seepage flow problems with free surfaces. A new parabolic variational inequality (PVI) formulation equivalent to the PDE formulation is then proposed, in which the flux part of the complementary condition of Signorini's type in the PDE formulation is transformed into natural boundary condition. Consequently, the singularity at the seepage points is eliminated and the difficulty in selecting the trial functions is significantly reduced. By introducing an adaptive penalized Heaviside function in the finite element analysis, the numerical stability of the discrete PVI formulation is well guaranteed. The proposed approach is validated by the existing laboratory tests with sudden rise and dropdown of water heads, and then applied to capture the non‐steady seepage flow behaviors in a homogeneous rectangular dam with five drainage tunnels during a linear dropdown of upstream water head. The non‐steady seepage flow in the surrounding rocks of the underground powerhouse in the Shuibuya Hydropower Project is further modeled, in which a complex seepage control system is involved. Comparisons with the in situ monitoring data show that the calculation results well illustrate the non‐steady seepage flow process during impounding and the operation of the reservoir as well as the seepage control effects of the drainage hole arrays and drainage tunnels. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
17.
含自由面的无压渗流问题本质上是一类非线性自由边值问题,固定网格的结点虚流量法在全域范围内不断扣除虚域流量贡献,从而使该问题得到求解。它具有网格依赖性小、出逸点收敛快等优点,但其内在理论基础尚未被完全揭示。通过引入互补型约束条件建立了结点虚流量法和Signorini型变分不等式提法的等价性桥梁,在此基础上引入过渡区放大系数 对自由面判别准则进行优化,并以砂槽模型试验为例进行验证。对比结果表明,优化后算法数值稳定性更好,计算结果与试验数据吻合度更高。研究成果为超大规模网格的渗控结构优化设计提供了有效分析手段。 相似文献
18.
提出了求解有自由面渗流问题的三维数值流形方法,通过构造任意形状流形单元的水头函数,推导了流形单元的渗透矩阵和无压渗流分析的总体控制方程,并给出了自由面的迭代求解策略和渗透体积力的计算方法。典型算例的数值分析表明,该方法采用数学网格覆盖整个材料区域,在自由面的迭代求解过程中数学网格保持不变,只考虑自由面以下渗流区的介质,只对自由面以下的流形单元形成总体渗透矩阵,具有精度高、收敛速度快、编程简单等优点,而且能够通过单纯形积分精确计算被自由面穿越单元的渗透作用力,因此,特别适用于有自由面渗流问题的模拟。 相似文献
19.
剖面二维裂隙网络渗流计算方法 总被引：6，自引：2，他引：6
王恩志 《水文地质工程地质》1993,20(4):27-29
本文针对剖面裂隙网络渗流问题,采用全区域不变网络分析方法,并引用初流量概念,化非线性分析为线性分析,成功的解决了岩体剖面上有自由面的裂隙网络渗流计算问题。 相似文献
20.
面板堆石坝垂直缝破坏下三维渗流场有限元模拟 总被引：1，自引：0，他引：1
采用改进节点虚流量法求解无压稳定渗流场,并引入无厚度的裂缝模型对破坏的垂直缝渗流行为进行模拟,得到面板堆石坝裂缝渗流问题的有限元分析方法,并编制Fortran程序。以某混凝土面板堆石坝为例,计算了面板单一垂直缝破坏和多条垂直缝破坏条件下的三维渗流场,得到不同条件下渗流场的水头分布、浸润线以及渗漏量,系统分析了面板堆石坝在垂直缝破坏条件下的稳定渗流场规律和特点。结果表明,该方法能对渗流逸出点和浸润线进行准确定位,还能很好地模拟面板垂直缝破坏对坝区渗流场的影响,可以为面板堆石坝的接缝设计提供参考。 相似文献