Variational inequality-based particle discontinuous deformation analysis

Orienting the circular and rigid particle medium, the variational inequality-based discontinuous deformation analysis (DDA) is established. In the proposed DDA, the global stiffness matrix, the penalty parameters, and the open-close iteration are successfully avoided. The contact constraint is transferred into the problem of variational or quasi-variational inequalities. And explicit variational expression on the contact force is firstly established. To speed up the rate of solving contact force, on the basis of the two-stage prediction-correction method, we design a compatibility iteration algorithm (PPC-CI). The C++ code is developed in multicore environment through the open multi-processing (OpenMP) in order to take advantage of the parallelizable features of the new DDA. Numerical tests suggest that the presented DDA is effective and promising.     

A generalized rigid particle contact model for fracture analysis

Rigid particle models taking directly into consideration the physical mechanisms and the influence of the material meso‐structure have recently been developed for fracture studies of quasi‐brittle material such as concrete. The formulation of a generalized contact model for rigid particle simulations is presented in which the contact discretization is a model parameter. The contact model performance for different discretizations is evaluated for uniaxial tensile tests, for uniaxial compression tests and for a notched beam in mode I. Copyright © 2005 John Wiley & Sons, Ltd.     

On the Implementation of augmented Lagrangian method in the two‐dimensional discontinuous deformation Analysis

The discontinuous deformation analysis (DDA) is a discontinuum‐based method, which employs a penalty method to represent the contact between blocks. The penalty method is easy to be implemented in the program, but the contact constraint is only approximately satisfied. Penetrations between contacting blocks are unavoidable even if the penalty value is very large. To improve the contact precision in the DDA, an augmented Lagrangian method is introduced, which can make use of advantages of both the Lagrangian multiplier method and the penalty method. This paper provides a detailed implementation of the augmented Lagrangian method in the DDA program and compares it with the standard DDA on the computational efficiency and contact precision. Copyright © 2013 John Wiley & Sons, Ltd.     

接触摩擦问题的数值模拟

无网格伽辽金法(EFGM)可脱离单元的概念,特别适合岩体裂纹面的接触摩擦分析。基于EFGM,在裂纹面引入罚参数,通过迭代计算,得到裂纹面真实的应力状态,从而模拟闭合裂纹的粘接、滑移和张开行为,数值结果表明该方法是合理可行的。     

基于线性互补的非连续变形分析

非连续变形分析(DDA)方法是一种新的用来分析块体系统运动和变形的非连续介质数值计算方法。研究的核心工作是致力于对现有DDA接触问题处理方法的改进。DDA主要采用罚函数法和Lagrange乘子法处理接触问题,合理设定罚参数很困难,此外,因开闭迭代而引起的刚度矩阵的不连续变化也会导致收敛方面的困难。为避免引入罚参数及传统意义上的开闭迭代,用混合线性互补模型(LCDDA)对DDA方法进行了重新描述。在此基础上,综合基于非光滑分析的Newton法的局部平方收敛和最速下降法的全局线性收敛的优势,提出求解LCDDA模型的有效算法。根据上述思想及理论研究成果编制了完整的计算程序,算例计算结果证明了方法的精度及可行性。     

A Model of Point-to-Face Contact for Three-Dimensional Discontinuous Deformation Analysis

Summary The key to three-dimensional discontinuous deformation analysis (3D DDA) is a rigorous contact theory that governs the interaction of many three-dimensional blocks. This theory must provide algorithms to judge contact types and locations and the appropriate state of each contact, which can be open, sliding or locked. This paper presents a point-to-face contact model, which forms a part of the contact theory, to be used in 3D DDA. Normal spring, shear spring and frictional force submatrices are derived by vector analysis and the penalty method. Also given are the open-close iteration criteria and operations performed for different changes in contact state. Sliding at a contact can occur in any direction parallel to the contact face, as opposed to one of two directions in two-dimensional DDA. This point-to-face contact model has been implemented into a 3D DDA computer program, and numerical results from several test cases demonstrate the validity of the model and the capability of the program.     

Unified displacement boundary constraint formulation for discontinuous deformation analysis (DDA)

Displacement boundary constraints in discontinuous deformation analysis (DDA) are applied using stiff penalty springs. A co‐ordinate‐free formulation for displacement boundary constraints is presented here for DDA, which unifies previous derivations for points of fixity, and for points constrained to induce or prohibit block motion in specified directions as a function of location or time. Examples for each type of constraint are used to illustrate the behaviour of the algorithm and provide a link with previous formulations for each case. The new, unified formulation has five benefits: (1) simple to express algorithmically; (2) easy to program and verify; (3) penalty values in different directions may be chosen to allow fixed points, lines, curves or planes; (4) formulation works for 2D and 3D; (5) displacement constraint may be a function of time or location or both. Feedback in the algorithm may induce internal resonance in homogeneously deformable discrete elements used in DDA, and resonance in block‐to‐block contact interactions. Consequently, high mass problems with insufficient damping may suffer from excessive ‘vibrational hammering’, inducing physically implausible behaviour such as elastic rebound. Copyright © 2005 John Wiley & Sons, Ltd.     

A new model of edge-to-edge contact for three dimensional discontinuous deformation analysis

A robust contact theory can be regarded as key to three dimensional discontinuous deformation analyses (3D-DDA). Not only must this theory provide an efficient algorithm to judge the type and location of contacts but also it must be able to present comprehensive formulations for every kind of contact (open, sliding and locked contact). There are six types of contact in three dimensional discontinuous deformation analyses (vertex-to-vertex, vertex-to-edge, vertex-to-face, edge-to-edge, edge-to-face and face-to-face) that can be converted to vertex-to-face and edge-to-edge contacts. This paper presents a new model of edge-to-edge contact to three dimensional discontinuous deformation analyses (3D-DDA). This new model considers both kinds of edge-to-edge contact (cross-over and parallel edge-to-edge contact) and presents a criterion for inter-penetration. Sub matrices of normal and shear spring and friction force are derived by geometrical analysis and penalty method. This new model is implemented in a 3D-DDA computer programme, and the example results demonstrate the validity of the model.     

Numerical modeling of hydraulic fracture propagation,closure and reopening using XFEM with application to in‐situ stress estimation

In this paper, a fully coupled model is developed for numerical modeling of hydraulic fracturing in partially saturated weak porous formations using the extended finite element method, which provides an effective means to simulate the coupled hydro‐mechanical processes occurring during hydraulic fracturing. The developed model is for short fractures where plane strain assumptions are valid. The propagation of the hydraulic fracture is governed by the cohesive crack model, which accounts for crack closure and reopening. The developed model allows for fluid flow within the open part of the crack and crack face contact resulting from fracture closure. To prevent the unphysical crack face interpenetration during the closing mode, the crack face contact or self‐contact condition is enforced using the penalty method. Along the open part of the crack, the leakage flux through the crack faces is obtained directly as a part of the solution without introducing any simplifying assumption. If the crack undergoes the closing mode, zero leakage flux condition is imposed along the contact zone. An application of the developed model is shown in numerical modeling of pump‐in/shut‐in test. It is illustrated that the developed model is able to capture the salient features bottomhole pressure/time records exhibit and can extract the confining stress perpendicular to the direction of the hydraulic fracture propagation from the fracture closure pressure. Copyright © 2016 John Wiley & Sons, Ltd.     

A new contact algorithm in the material point method for geotechnical simulations

Contact between stiff structural elements and soil is encountered in many applications in geotechnical engineering. Modelling of such contact is challenging as it often involves impact that would lead to large deformation and failure of the soil. The Material Point Method (MPM) is a mesh‐free method that has been applied to simulate such phenomena. However, the frictional contact algorithm commonly used in MPM only supports Coulomb friction and cannot model fully or partially rough contact conditions in terms of geotechnical engineering. Moreover, because of very different stiffness of contacting materials, the contact force predicted by the previous frictional contact algorithms usually suffers from severe oscillation when applied in structure–soil interaction. This paper presents a new contact algorithm, termed Geo‐contact, designed for geotechnical engineering. In Geo‐contact, a penalty function is incorporated to reduce the oscillation in contact computation, and a limited shear stress is specified along the contact interface. The proposed Geo‐contact algorithm has been implemented to simulate smooth, partially rough and rough contact in typical large deformation penetration problems. The resistance–displacement curves obtained using the Geo‐contact are compared with analytical solutions of limit analysis and large deformation finite element results to verify the accuracy and robustness of the proposed contact algorithm. Copyright © 2014 John Wiley & Sons, Ltd.     

Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method

A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Parametric study of the smooth‐joint contact model on the mechanical behavior of jointed rock

The smooth‐joint contact model based on distinct element method has been widely used to represent discontinuity in the simulation of fractured rock mass, but there is rare efficient guidance for the selection of proper parameters of smooth‐joint contact model, which is the basement for using this model properly. In this paper, the effect of smooth joint parameters on the macroscopic properties and failure mechanism of jointed rock under triaxial compression test is investigated. The numerical results reveal that the friction coefficient of smooth joint plays a dominant role in controlling mechanical behaviors. The stiffness of smooth joint has a relative small influence on the mechanical behaviors. Poisson ratio decreases with the reduction of normal stiffness but increases with the reduction of shear stiffness. The reduction of smooth joint strength, which is determined by normal strength, cohesion, and friction angle of smooth joint, contributes to the breakage of bonded smooth joint and ultimately decreases the strength of the specimen. We proposed a detailed calibration process for smooth‐joint contact model according to the relationship between smooth‐joint parameters and mechanical properties. By following this process, the numerical results are validated against corresponding experimental results and good agreement between them can be found in stress‐strain curves and failure modes of different joint orientations. Further analyses from the microperspective are performed by looking at transmission of contact force, the nature and distribution of microcracks, and the particle displacement to show the failure process and failure modes.     

Numerical evaluation of soft inter‐slab joint in concrete‐faced rockfill dam with dual mortar finite element method

Recently constructed concrete‐faced rockfill dams (CFRDs) often use soft inter‐slab joints to prevent axial compression‐induced extrusion damage in the concrete face. Due to the complexity of the multibody contact and the lack of information on the actual behavior of soft joints, it is highly challenging to numerically assess the effect of soft joints in CFRDs. In this paper, we present a numerical approach for the three‐dimensional modeling of CFRDs with hard and soft joints. A dual mortar finite element method with Lagrange multiplier is developed to treat the multibody contact in hard joints with impenetrability condition. The soft joint slab‐filler‐slab contact system is modeled using an equivalent contact interface approach, where the soft contact constraints are imposed using a perturbed Lagrange formulation. Through a series of laboratory tests, the mechanical behavior of soft joint is investigated. An extrusion model for the soft joint is presented and implemented in the dual mortar finite element method. The proposed numerical method is applied to the three‐dimensional analysis of Tianshengqiao‐1 CFRD. Despite the complex multibody contact and strong material and geometry nonlinearities in the CFRD, the proposed method is stable and capable of capturing salient characteristics of the CFRD. Numerical results show that in Tianshengqiao‐1, the employment of soft joints can effectively reduce the axial compression stress, thus greatly alleviating the risk of extrusion damage in the concrete face.     

Non-rigid disk-based DDA with a new contact model

In this paper, a new disk-based DDA formulation is presented. In the original disk-based DDA, disks are considered to be rigid and the penalty method is used to enforce disk contact constraints. In order to improve the accuracy of the disk-based DDA, new formulations of stiffness and force matrices for non-rigid disks using a new efficient contact model are presented in this paper. Blocks are considered deformable without need to do more computations for contact detection. In the proposed contact model, disk–disk and disk–boundary contacts are transformed into the form of point-to-line contacts and normal spring, shear spring and frictional force sub-matrices are derived by vector analysis. The penalty method is quite simple to implement, but has some major disadvantages. In the presented contact model, not only the simplicity of the penalty method is retained but also the limitations are overcome by using the augmented Lagrangian method. Moreover, unlike the contact model used in the original disk-based DDA, reference line can be obtained directly by using only coordinates of disk centers and their radii, and no more computations are needed. The validity and capability of the new disk-based DDA formulation are demonstrated by several illustrative examples.     

A development of the discontinuous deformation analysis for rock fall analysis

Discontinuous deformation analysis (DDA), a discrete numerical analysis method, is used to simulate the behaviour of falling rock by applying a linear displacement function in the computations. However, when a block rotates, this linear function causes a change in block size called the free expansion phenomenon. In addition, this free expansion results in contact identification problems when the rotating blocks are close to each other. To solve this problem of misjudgment and to obtain a more precise simulation of the falling rock, a new method called Post‐Contact Adjustment Method has been developed and applied to the program. The basic procedure of this new method can be divided into three stages: using the linear displacement function to generate the global matrix, introducing the non‐linear displacement function to the contact identification, and applying it to update the co‐ordinates of block vertices. This new method can be easily applied to the original DDA program, demonstrating better contact identification and size conservation results for falling rock problems than the original program. Copyright © 2005 John Wiley & Sons, Ltd.     

重,磁界面计算的广义反演方法

介绍一种重、磁界面计算的广义反演方法。该方法的主要特点是引入了模型参数的协方差矩阵使得反演过程能够利用模型参数的先验信息而加快收敛速度。文中首先给出了失代公式，并给出了详细算法。其次，通过理论模型的计算对方法作了测试。根据这些测试结果，重点讨论了Ｃｍ在反演中的作用。最后给出了一个应用实例，并与频率域方法作了比较。     

u‐p semi‐Lagrangian reproducing kernel formulation for landslide modeling

This paper presents a u‐p (displacement‐pressure) semi‐Lagrangian reproducing kernel (RK) formulation to effectively analyze landslide processes. The semi‐Lagrangian RK approximation is constructed based on Lagrangian discretization points with fixed kernel supports in the current configuration. As a result, it tracks state variables at discretization points while allowing extreme deformation and material separation that is beyond the capability of Lagrangian formulations. The u‐p formulation following Biot theory is incorporated into the formulation to describe poromechanics of saturated geomaterials. In addition, a stabilized nodal integration method to ensure stability of the domain integration and kernel contact algorithms to model contact between bodies are introduced in the u‐p semi‐Lagrangian RK formulation. The proposed method is verified with several numerical examples and validated with an experimental result and the field data of an actual landslide.     

Three‐dimensional numerical analysis of concrete‐faced rockfill dam using dual‐mortar finite element method with mixed tangential contact constraints

Concrete‐faced rockfill dam (CFRD) is a popular alternative to traditional dam types in the last two decades. The modelling of CFRD involves complex multi‐body contact and strong geometry and material nonlinearities. We present a numerical approach for the modelling of CFRDs in this paper. Based on the dual‐mortar finite element method, the presented approach considers different parts of rockfill and all concrete slabs as independent deformable continuum. The multi‐body contacts are modelled using Lagrange multipliers with a weak form segment‐to‐segment contact strategy. To alleviate instability induced by strong geometry nonlinearity in the slab–slab contact, we propose a mixed type of constraints for the tangential contact. A general transformation scheme is introduced to simplify the implementation of contact constraints. Three‐dimensional analysis of Tianshengqiao‐1 CFRD is performed. The nonlinear and time‐dependent deformation of the rockfill is considered. We study the influence of the rockfill deformation on the reliability of the concrete face. Three major concerns of the face, that is, the axial compression, the slab–slab separation and the face‐rockfill separation, are discussed in detail. The numerical results are compared with data from in‐situ observation. Copyright © 2016 John Wiley & Sons, Ltd.