A nonlinear programming analysis of unconfined steady-state seepage

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.     

饱和模型非定常渗流的数值求解方法

提出了饱和模型非稳定流的有限元算法，并对自由面流量补给项提出一种简单的处理方法，从而避免了在迭代中求解自由面的位置，最后给出了计算实例。     

Unconfined seepage analysis in earth dams using smoothed fixed grid finite element method

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.     

Kinematic modelling of shear band localization using discrete finite elements

Modelling shear band is an important problem in analysing failure of earth structures in soil mechanics. Shear banding is the result of localization of deformation in soil masses. Most finite element schemes are unable to model discrete shear band formation and propagation due to the difficulties in modelling strain and displacement discontinuities. In this paper, a framework to generate shear band elements automatically and continuously is developed. The propagating shear band is modelled using discrete shear band elements by splitting the original finite element mesh. The location or orientation of the shear band is not predetermined in the original finite element mesh. Based on the elasto‐perfect plasticity with an associated flow rule, empirical bifurcation and location criteria are proposed which make band propagation as realistic as possible. Using the Mohr–Coulomb material model, various results from numerical simulations of biaxial tests and passive earth pressure problems have shown that the proposed framework is able to display actual patterns of shear banding in geomaterials. In the numerical examples, the occurrence of multiple shear bands in biaxial test and in the passive earth pressure problem is confirmed by field and laboratory observations. The effects of mesh density and mesh alignment on the shear band patterns and limit loads are also investigated. Copyright © 2003 John Wiley & Sons, Ltd.     

基于二次场二维起伏地形MT有限元数值模拟

通过计算二次场来进行二维大地电磁数值模拟;导出了二维大地电磁二次场的微分方程,利用有限单元法来解微分方程;对矩形网格进行对角线的二次剖分,更容易且真实地模拟起伏地形。对几个典型模型进行了试算,与前人总场法的计算结果做了比较,两者视电阻率曲线一致,证明本文算法是正确的;通过2个简单的算例说明复杂地表下2种极化模式的MT观测资料都有明显的异常,视电阻率在TM模式下比TE模式更易受地形影响,TE模式下视电阻率曲线形态与地形呈"正相关",TM模式下反之。     

基于互补理论的扩展有限元接触问题实现

在处理含有裂隙这类不连续问题时,常规有限元方法需要对裂隙尖端部位进行局部网格加密,当裂隙扩展时还需要进行网络重构。基于单位分解思想的扩展有限元成功解决了常规有限元难以处理的裂隙类不连续问题。在研究复合裂隙时,通常需要考虑裂隙的接触问题。基于互补理论,建立了裂隙面上相对位移和接触力的互补方程,并采用牛顿法求解,无需开闭迭代,且能够快速收敛。最后,对含裂隙平板进行受压数值试验,计算结果表明,基于互补理论的扩展有限元接触算法能够有效地阻止裂隙两端网格的相互嵌入,且获得裂隙面上的应力分布与实际一致。     

岩体裂隙网络非稳定渗流分析与数值模拟

针对裂隙岩体的非稳定渗流问题,通过将Darcy定理扩展到包含干区的整个裂隙网络区域,并令潜在溢出边界条件为Signorini型互补边界条件,将湿区上的非稳定渗流问题转化为全域上的一个新的初边值问题。为降低试探函数选取的难度,建立与定义在整个裂隙网络区域上的偏微分方程（PDE）提法等价的抛物型变分不等式（PVI）提法,并给出裂隙网络非稳定渗流分析的有限元数值分析格式和迭代算法,与砂槽模型试验数据的对比分析,验证其有效性。最后,将文中发展的计算方法应用到含复杂裂隙网络的边坡非稳定渗流分析,计算结果很好地反映出边坡内部自由面随库水降落的变化规律,并能准确地描述裂隙网络内部渗流运动特征及流量分布的不均匀性。     

Slipping strip analysis of reinforced earth

Reinforced earth in plane strain is idealized as a homogeneous material with the strips attached to the elastic soil matirx by a conceptual shear zone. A ‘no-slip’ finite element model is derived by assigning a large shear modulus to the shear zone. Relaxation of this modulus using a tangential stiffness algorithm in conjunction with a Mohr–Coulomb strip-slip criterion allows slipping to be simulated. The finite element formulation is validated and the finite element discretization assessed by comparisons against exact solutions for a simple test problem. An idealized reinforced earth wall example is used to demonstrate the feasibility of the method and to answer the question: ‘is slipping significant?’ The method is shown to be potentially useful, and slipping is shown to be significant.     

基于CAD地形图的初始地应力场反演快速建模程序开发

基于Visual C#2008.NET开发环境,利用AutoCAD ActiveX技术,对AutoCAD2008进行二次开发,开发出快速建立初始地应力场反演计算三维有限元模型的面向对象程序,通过该程序可以根据需要自动读取一定范围内的网格节点地面高程,能够为有限元模型变地面线处理直接利用,避免了常规人工读取CAD地面高程的繁琐性和误差较大的问题,极大地提高了工作效率和准确性,快速便捷。可以快速自动建立用于初始地应力反演计算的有限元网格模型,该模型能够较好地反映地形貌地的起伏状况,为有限元计算提供了良好的基础条件。     

模拟三维裂纹问题的扩展有限元法

扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其计算网格与不连续面相互独立,因此模拟移动不连续面时无需对网格进行重新剖分。给出了模拟三维裂纹问题的扩展有限元法。在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性。用两个水平集函数表示裂纹。采用线性互补法求解裂纹面非线性接触条件,不需要迭代,提高了计算效率。采用两点位移外推法计算裂纹前缘应力强度因子。给出了3个三维弹性静力问题算例,其结果显示了所提方法能获得高精度的应力强度因子,并能有效地处理裂纹面间的接触问题,同时表明扩展有限元结合线性互补法求解不连续问题具有较好的前景。     

A finite volume approach with local adaptation scheme for the simulation of free surface flow in porous media

We present a method for solving steady‐state flow with a free surface in porous media. This method is based on a finite volume approach and is halfway between a fixed and an adaptive mesh method, taking advantage of both approaches: computational efficiency and localization accuracy. Most of the mesh remains fixed during the iterative process, while the cells in contact with the free surface (free surface cells) are being reshaped. Based on this idea, we developed two methods. In the first one, only the volumes of the free surface cells are adapted. In the second one, the computational nodes of the free surface cells are relocated exactly at the free surface. Both adaptations are designed for a better application of the free surface boundary conditions. Implementation details are given on a regular finite volume mesh for the case of homogeneous and heterogeneous rectangular dams in 2D and 3D. Accuracy and convergence properties of the proposed approach are demonstrated by comparison with an analytical solution and with existing references. Copyright © 2011 John Wiley & Sons, Ltd.     

地下厂房岩锚梁锚杆涂沥青段有限元模拟研究与应用

基于AutoCAD和OpenGL有限元可视化建模,利用Fortran语言对模型文件中的单元信息进行改写,提出了一种在地下厂房有限元模型中生成沥青单元网格的方法,显著减少了模型单元数与节点数,提高了建模与计算效率。采用隐式杆单元法模拟涂沥青锚杆,隐式柱单元法模拟普通锚杆,推导了其刚度矩阵及有限元迭代计算公式,并应用于地下厂房岩锚梁锚杆涂沥青段的有限元分析中。结果表明,吊车梁附近围岩破坏区减少,围岩深部锚杆应力增大,说明围岩受力特性得以改善,围岩和吊车梁整体安全度增大,符合工程实践经验。所提出的计算和单元再分方法可操作性和实用性大,可为类似工程有限元分析提供参考。     

裂隙网络无压渗流分析的初流量法

借鉴初流量法的思想,引入初流速来抵消在裂隙网络干区实际不存在的流速,将Darcy定理扩展到整个区域,定义了在整个区域上的非线性边值问题,并将潜在溢出面边界条件归纳为Signorini型边界条件,建律了等价的变分不等式提法。通过结合连续型的Heaviside函数,给出了基于变分不等式的初流量有限元迭代算法,研制开发了相应的计算程序,通过两个典型算例验证了本文算法在求解复杂裂隙网络渗流自由面的有效性和鲁棒性。     

Use of the finite element method in modelling a static earth pressure problem

Results obtained with the finite element method are greatly affected by the overall dimensions of the mesh and the number and size of the elements used. Proper choice of the boundary conditions is also necessary to obtain an accurate estimate of the unknowns of the problem. This paper discusses some factors that influence the results of a finite element idealization of the problem of earth pressure behind a gravity wall with a dry, cohesionless backfill.     

Coupled computations for an elastic‐perfectly plastic soil using adaptive mesh refinement

An adaptive mesh refinement algorithm has been developed for non‐linear computations in geomechanics, based on a smoothed stress–strain finite element formulation. This uses estimates of error in the incremental shear strain invariant to guide the regeneration of unstructured meshes at regular intervals during loading. Following each mesh‐update, no re‐analysis of previous increments with the new mesh is necessary. Algorithm performance has been investigated by analysing a passive earth pressure problem using a linear elastic‐perfectly plastic Mohr–Coulomb soil model. Perfectly drained behaviour has been considered, as have partially drained situations using hydromechanical coupling, while undrained behaviour has been approximated using time steps close to zero. In all cases, mesh adaptivity has been successful in capturing regions of high strain gradient. The results have been compared with analytical solutions. Accurate computations of limit load and shear band orientation have been obtained for a wide range of material dilation angles. Copyright © 2000 John Wiley & Sons, Ltd.     

多尺度有限单元法求解非均质多孔介质中的三维地下水流问题

应用多尺度有限单元法模拟非均质多孔介质中的三维地下水流问题。与传统有限单元法相比，多尺度有限单元法的基函数具有能反映单元内参数变化的优点，所以这种方法能在大尺度上抓住解的小尺度特征获得较精确的解。在介绍多尺度有限单元法求解非均质多孔介质中三维地下水流问题的基本原理之后，对参数水平方向渐变垂直方向突变的非均质多孔介质中的三维地下水流和Borden实验场的三维地下水流分别用多尺度有限单元法和传统等参有限单元法进行了计算，结果表明在模拟高度非均质多孔介质中的三维地下水流问题时，多尺度有限单元法比传统有限单元法有效，既节省计算量又有较高的精度；在模拟非均质性弱的多孔介质中的三维地下水流问题时，多尺度有限单元法虽然也能在大尺度上获得较为精确的解，但效果不明显。     

俯斜式混凝土重力挡土墙强震反应数值模拟

首先,介绍了基于OpenSees独立开发的一套用于挡土墙-土地震反应相互作用有限元分析计算软件RW_2DPS.据此建立了俯斜式混凝土重力挡土墙-土强震相互作用有限元模型.模型中,引入非线性有限元计算方法,选用多屈服面弹塑性本构模型模拟砂土的动力属性,应用零长度接触单元模拟墙与土体之间的接触特性,且采用一致耗能阻尼边界与速度边界条件.最后,输入随机地震动,进行挡土墙-土强震反应分析,并重点探讨墙背地震土压力和水平地震惯性力沿挡土墙高度分布规律.结果表明,墙背动土压力峰值出现在距挡土墙底约1/3墙高处;挡土墙背加速度具有放大效应,加速度峰值出现在挡土墙顶部;不同地震动作用下,加速度放大系数沿墙高分布规律不同,动土压力沿墙高变化规律基本一致.     

Algorithm of Coupled Normal Stress and Fluid Flow in Fractured Rock Mass by the Composite Element Method

This paper presents a composite element algorithm of coupled normal stress and fluid flow process for fractured rock mass, developed from the composite element method (CEM). The coupled relation between the fracture flow and normal stress makes use of the “filled model”, which examines the asperities in the fracture as a layer of granular medium having high porosity and being clipped by the two parallel plates. The existence of fractures is not considered in the mesh generation, but it will be considered explicitly in the mapped composite element. The coupled normal stress and fluid flow process has been simulated by applying a cross iterative algorithm between the two fields. The proposed algorithm considers not only the flow through the fractures, but also the flow exchange between fractures and the surrounding rock blocks. In addition, it can be used for both the filled and non-filled fractures. The verification of the proposed algorithm has been conducted through the illustration of three examples by comparison with the conventional finite element method (FEM), from which the advantages and reliability of the proposed algorithm have been shown clearly.