A new finite element formulation for one-dimensional analysis of elastic-plastic materials

It is well accepted that severe numerical difficulties arise when using the conventional finite element displacement method to analyse incompressible, or nearly incompressible, solids. These effects are caused by the kinematic constraints imposed on the nodal velocities by the constant volume condition. In elastic-plastic analysis, these effects are due to a conflict between the plastic flow rule and the finite element discretization. Although several methods have been proposed to cope with this problem, none has been based on the appropriate choice of displacement interpolation to minimise the constraints. In this paper, a new displacement interpolation, which is able to reduce the imposed constraints, is adopted. Comparisons of the results with those from a conventional linear displacement interpolation are made for predictions of cylindrical and spherical cavity expansion limit pressures in elastic-plastic solids. This study suggests that the proposed displacement interpolation is preferable to the conventional one in the elastic-plastic finite element analysis of one dimensional-axisymmetric problems which involve nearly incompressible material behaviour.     

A note on finite element simulations of pile driving

The analyses of pile driving using solid, axisymmetric finite elements reported in Reference 1 are extended to different soil conditions and hammer characteristics. Correlation with ‘pile analyser’ procedures is made.     

Variational principles and finite element simulations for thermo-elastic consolidation

In this paper, a general variational principle for the initial boundary value problem of quasi-static thermoelastic consolidation is developed by assuming infinitesimal deformation and an incompressible fluid flowing through a linearly elastic solid. By manipulating the coupling operators, an extended form of the variational pronciple is derved. The associated finite element formulation based on this principle is presented and numerical applications for plane strain thermo-elastic consolidation are revealed.     

Fast 2.5D finite element simulations of borehole resistivity measurements

We develop a rapid 2.5-dimensional (2.5D) finite element method for simulation of borehole resistivity measurements in transversely isotropic (TI) media. The method combines arbitrary high-order \(H^{1}\)- and \(\mathbf {H}\)(curl)-conforming spatial discretizations. It solves problems where material properties remain constant along one spatial direction, over which we consider a Fourier series expansion and each Fourier mode is solved independently. We propose a novel a priori method to construct quasi-optimal discretizations in physical and Fourier space. This construction is based on examining the analytical (fundamental) solution of the 2.5D formulation over multiple homogeneous spaces and assuming that some of its properties still hold for the 2.5D problem over a spatially heterogeneous formation. In addition, a simple parallelization scheme over multiple measurement positions provides efficient scalability. Our method yields accurate borehole logging simulations for realistic synthetic examples, delivering simulations of borehole resistivity measurements at a rate faster than 0.05 s per measurement location along the well trajectory on a 96-core computer.     

A stabilized extended finite element framework for hydraulic fracturing simulations

We present a stabilized extended finite element formulation to simulate the hydraulic fracturing process in an elasto‐plastic medium. The fracture propagation process is governed by a cohesive fracture model, where a trilinear traction‐separation law is used to describe normal contact, cohesion and strength softening on the fracture face. Fluid flow inside the fracture channel is governed by the lubrication equation, and the flow rate is related to the fluid pressure gradient by the ‘cubic’ law. Fluid leak off happens only in the normal direction and is assumed to be governed by the Carter's leak‐off model. We propose a ‘local’ U‐P (displacement‐pressure) formulation to discretize the fluid‐solid coupled system, where volume shape functions are used to interpolate the fluid pressure field on the fracture face. The ‘local’ U‐P approach is compatible with the extended finite element framework, and a separate mesh is not required to describe the fluid flow. The coupled system of equations is solved iteratively by the standard Newton‐Raphson method. We identify instability issues associated with the fluid flow inside the fracture channel, and use the polynomial pressure projection method to reduce the pressure oscillations resulting from the instability. Numerical examples demonstrate that the proposed framework is effective in modeling 3D hydraulic fracture propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

基于刚体平动运动单元的上限有限元研究

刚性块体极限分析上限法常应用于岩土工程稳定性研究,然而应用时需假定刚性块体破坏模式并递推繁琐的几何关系。为此,提出一种适应性更广的基于非线性规划模型的刚体平动运动单元上限有限元法,并解决了其优化模型初始值的确定问题。通过引入有限单元思想,将计算区域离散成刚体单元,同时以单元速度和节点坐标作为决策变量,由上限定理建立非线性规划模型获得上限解。利用编制的上限有限元程序进行边坡和浅埋隧道稳定性算例验证,表明运动单元上限有限元法能调整速度间断线至较优方位,所得破坏模式特征鲜明,上限解精度高,可广泛应用于边坡、隧道等稳定性分析研究。     

A stabilized assumed deformation gradient finite element formulation for strongly coupled poromechanical simulations at finite strain

An adaptively stabilized finite element scheme is proposed for a strongly coupled hydro‐mechanical problem in fluid‐infiltrating porous solids at finite strain. We first present the derivation of the poromechanics model via mixture theory in large deformation. By exploiting assumed deformation gradient techniques, we develop a numerical procedure capable of simultaneously curing the multiple‐locking phenomena related to shear failure, incompressibility imposed by pore fluid, and/or incompressible solid skeleton and produce solutions that satisfy the inf‐sup condition. The template‐based generic programming and automatic differentiation (AD) techniques used to implement the stabilized model are also highlighted. Finally, numerical examples are given to show the versatility and efficiency of this model. Copyright © 2013 John Wiley & Sons, Ltd.     

Surge simulations for 1999 Orissa super cyclone using a finite element model

The Orissa super cyclone which crossed the Orissa coastal region near Paradip on October 29, 1999 proved to be disastrous. The strong winds, torrential rains with heavy rainfall and high storm surge associated with the cyclone caused havoc that resulted in the death of thousands of people, cattle and extensive damage to agricultural land, paddy crop, transmission lines, power supply, roads and buildings. In the present study, a fine resolution finite element model is used to simulate surges due to this super cyclone. The model results are compared with observed surges available from post-storm survey reports. Comparison of results show that they are in good agreement with the observed surges, and the finite element model can be used for real time surge forecasts.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

基于粗细网格的有限元并行分析方法

并行计算己成为求解大规模岩土工程问题的一种强大趋势。探讨了粗细网格与预处理共轭梯度法结合的并行有限元算法。从多重网格刚度矩阵推得有效的预处理子。该算法在工作站机群上实现。用地基处理时土体强夯的数值模拟分析进行了数值测试,对其并行性能进行了详细分析。计算结果表明：该算法具有良好的并行加速比和效率,是一种有效的并行算法。     

Modeling discontinuous rock mass based on smoothed finite element method

This paper aims at developing a method for modeling rock mass with preexisting multiple discontinuities within the framework of the smoothed finite element method (SFEM). The discontinuity is simulated by an interface element with zero thickness, the stiffness matrix of which are derived explicitly based on the SFEM. An elastic damage constitutive relation with residual strength is introduced in order to describe the nonlinear mechanical behavior of the discontinuities. The computation codes of the present method were developed. The present method has been verified to be a sound approach for modeling discontinuous rock mass, inheriting the advantages of the SFEM.     

铜王公路2#滑坡侵蚀弹塑性有限元数值模拟

本研究采用铜王公路2#滑坡侵蚀实例,定性分析了该滑坡侵蚀的诱发因素、变形特征、复活的原因;并在研究该滑坡侵蚀地貌特征、岩土体组成、水文地质特征后,建立了地质模型和计算模型;采用弹塑性有限元法对斜坡的初始应力场特征进行模拟研究。计算成果及分析表明,铜王公路2#滑坡处于非稳定状态。最后,对该滑坡滑动后可能造成的灾害程度及侵蚀量进行评估。      

有限单元法分析边坡稳定的若干问题研究

系统介绍了有限单元法分析边坡稳定的基本原理,对安全指标与判断标准两个重要方面进行了详细论述,并采用面向对象技术编制了相应的可视化分析软件AFEAS。结合一个边坡稳定工程实例,对比了有限单元法和刚体极限平衡法的分析结果。分析结果证明了所给的安全指标与判断标准是可行的,并表明了有限单元法分析边坡稳定的优越性。     

A basal slip model for Lagrangian finite element simulations of 3D landslides

A Lagrangian numerical approach for the simulation of rapid landslide runouts is presented and discussed. The simulation approach is based on the so‐called Particle Finite Element Method. The moving soil mass is assumed to obey a rigid‐viscoplastic, non‐dilatant Drucker–Prager constitutive law, which is cast in the form of a regularized, pressure‐sensitive Bingham model. Unlike in classical formulations of computational fluid mechanics, where no‐slip boundary conditions are assumed, basal slip boundary conditions are introduced to account for the specific nature of the landslide‐basal surface interface. The basal slip conditions are formulated in the form of modified Navier boundary conditions, with a pressure‐sensitive threshold. A special mixed Eulerian–Lagrangian formulation is used for the elements on the basal interface to accommodate the new slip conditions into the Particle Finite Element Method framework. To avoid inconsistencies in the presence of complex shapes of the basal surface, the no‐flux condition through the basal surface is relaxed using a penalty approach. The proposed model is validated by simulating both laboratory tests and a real large‐scale problem, and the critical role of the basal slip is elucidated. Copyright © 2016 John Wiley & Sons, Ltd.     

A dynamic large deformation finite element method based on mesh regeneration

In this paper, a large deformation finite element (LDFE) approach termed ‘remeshing and interpolation technique with small strain (RITSS)’ is extended from static to dynamic soil-structure interaction applications. In addition, a technique termed ‘element addition’ is developed to improve the computational efficiency of both static and dynamic LDFE analyses that involve moving boundaries. The RITSS approach is based on frequent mesh generation to avoid element distortion. In dynamic RITSS, the field variables mapped from the old to the new mesh involve not only the stresses and material properties, but also the nodal velocities and accelerations. Using the element addition technique, new soil elements are attached to the domain boundaries periodically when the soil near the boundaries becomes affected by large displacements of the structure. The procedures of this Abaqus-based dynamic LDFE analysis and element addition technique are detailed, and the robustness of the techniques is validated and assessed through three example analyses: penetration of a flat footing into a half-space and movement of rigid and deformable landslides down slopes.     

基于ABAQUS平台的水力裂缝扩展有限元模拟研究

随着扩展有限元理论的深入研究,利用扩展有限元方法模拟水力压裂具有了一定的可操作性。相比于常规有限元方法,XFEM方法具有计算结果精度高和计算量小的优点。但是,如何模拟射孔孔眼、如何模拟流体与岩石相互作用以及分析水力裂缝的扩展规律仍然是难题。以研究水力压裂裂缝扩展规律为目的,建立了岩石多孔介质应力平衡方程、流体渗流连续性方程和边界条件。通过有限元离散化方法对耦合方程矩阵进行处理。通过富集函数定义初始裂缝（射孔孔眼）,选择最大主应力及损伤变量D分别作为裂缝起裂和扩展判定准则,利用水平集方法模拟水力裂缝扩展过程。数值模拟结果显示：增加射孔方位角、压裂液排量和减小水平地应力差,起裂压力上升;黏度对起裂压力无明显影响。增加射孔方位角、压裂液排量、黏度和减小水平地应力差值有助于裂缝宽度的增加。增加水平地应力差值、压裂液排量和减小射孔方位角以及压裂液黏度有助于裂缝长度增加,反之亦然。基于ABAQUS的水力裂缝扩展有限元法可对不同井型和诸多储层物性参数及压裂施工参数进行分析,且裂缝形态逼真,裂缝面凹凸程度清晰,结果准确。此研究可作为一种简便有效研究水力压裂裂缝扩展规律的方法为油田水力压裂设计与施工提供参考与依据。     

基于有限元数值计算的爆破震动强度分析

爆破震动强度的预测一直是爆破震动危害控制的前提和基础,也是其研究的一个重要课题。基于有限元数值模拟方法的基础上对反映爆破震动强度的重要指标--质点振动速度的计算方法进行了探讨,并对空气和水两种不同间隔装药结构下的震动强度计算结果进行了对比分析。数值模拟结果表明,在爆心距不太大的情况下,计算结果与实测结果吻合较好,研究结果对爆破设计与施工有一定的指导意义。     

基于方位离散线性化的上限原理有限元法

上限原理有限元法不仅可以得到边坡的安全系数,还可以给出临界滑动面,且具有比极限平衡法更严谨的理论基础,因此,拥有更广阔的应用前景。针对传统的上限有限元法不能考虑强度各向异性的问题,提出了一种新的摩尔-库仑屈服面线性化方法。该方法在对方位角离散化的基础上,建立了线性化的方位离散塑性流动约束方程,丰富了基于线性规划的上限法理论。两个算例结果表明：该方法可以稳定地从极限解的上方收敛;且对边坡进行稳定性分析,若忽略了边坡的强度各向异性,则会高估边坡的稳定性,得到较大的安全系数。     

基于非结构化有限元的三维井地电阻率法约束反演

电磁探测反演是典型的不适定问题,易造成反演结果的多解性,不适定性是反演自身固有的特征,没有求解的附加信息这一本质困难是很难克服的,解决该问题的有效方法是研究约束反演。本文采用目前较为主流的高斯牛顿—共轭梯度法(GN-CG),在反演目标函数中直接施加约束条件,将介质电阻率的取值范围作为先验信息和约束条件以外点罚函数法的方式引入到反演目标函数中,与常规三维电阻率反演目标函数相比,增加了不等式约束项的目标函数,理论上可以压制反演的多解性。通过多种理论模型的测试结果表明,本文基于不等式约束的三维井地电阻率反演算法有效地改善了反演结果的精度,以惩罚函数法施加不等式约束条件的方式是现实可行及有效的。     

基于三维非线性有限元的边坡稳定分析方法

刚体极限平衡法不能反映岩体中实际的应力分布,而基于有限元的强度折减系数法在判断收敛性方面存在一些问题。为了解决这些问题,采用多重网格法,分别建立用于有限元计算的结构网格和用于计算滑面稳定安全系数的滑面网格,可以方便地获得任意滑面或滑块的稳定安全系数,从而将非线性有限元和极限平衡分析结合起来。为了提高计算规模和计算精度,采用有限元并行计算程序TFINE.Pfem进行计算,分析了网格密度对计算结果精度的影响,并应用于锦屏高边坡的稳定分析中。与刚体极限平衡法结果的对比分析表明,由于考虑了计算过程中的非线性应力调整,该方法的计算结果比刚体极限平衡法偏大,而且更符合实际情况。