Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

2D numerical manifold method based on quartic uniform B-spline interpolation and its application in thin plate bending

A new numerical manifold （NMM） method is derived on the basis of quartic uniform B-spline interpolation. The analysis shows that the new interpolation function possesses higher-order continuity and polynomial consistency compared with the conven- tional NMM. The stiffness matrix of the new element is well-conditioned. The proposed method is applied for the numerical example of thin plate bending. Based on the prin- ciple of minimum potential energy, the manifold matrices and equilibrium equation are deduced. Numerical results reveal that the NMM has high interpolation accuracy and rapid convergence for the global cover function and its higher-order partial derivatives.

Numerical manifold method modeling of coupled processes in fractured geological media at multiple scales

The greatest challenges of rigorously modeling coupled hydro-mechanical(HM)processes in fractured geological media at different scales are associated with computational geometry.These challenges include dynamic shearing and opening of intersecting fractures at discrete fracture scales as a result of coupled processes,and contact alteration along rough fracture surfaces that triggers structural and physical changes of fractures at micro-asperity scale.In this paper,these challenges are tackled by developing a comprehensive modeling approach for coupled processes in fractured geological media based on numerical manifold method(NMM)at multiple scales.Based on their distinct geometric features,fractures are categorized into three different scales:dominant fracture,discrete fracture,and discontinuum asperity scales.Here the scale is relative,that of the fracture relative to that of the research interest or domain.Different geometric representations of fractures at different scales are used,and different governing equations and constitutive relationships are applied.For dominant fractures,a finite thickness zone model is developed to treat a fracture as a porous nonlinear domain.Nonlinear fracture mechanical behavior is accurately modeled with an implicit approach based on strain energy.For discrete fractures,a zero-dimensional model was developed for analyzing fluid flow and mechanics in fractures that are geometrically treated as boundaries of the rock matrix.With the zero-dimensional model,these fractures can be modeled with arbitrary orientations and intersections.They can be fluid conduits or seals,and can be open,bonded or sliding.For the discontinuum asperity scale,the geometry of rough fracture surfaces is explicitly represented and contacts involving dynamic alteration of contacts among asperities are rigorously calculated.Using this approach,fracture alteration caused by deformation,re-arrangement and sliding of rough surfaces can be captured.Our comprehensive model is able to handle the computational challenges with accurate representation of intersections and shearing of fractures at the discrete fracture scale and rigorously treats contacts along rough fracture surfaces at the discontinuum asperity scale.With future development of three-dimensional(3D)geometric representation of discrete fracture networks in porous rock and contacts among multi-body systems,this model is promising as a basis of 3D fully coupled analysis of fractures at multiple scales,for advancing understanding and optimizing energy recovery and storage in fractured geological media.

Searching for critical slip surfaces of slopes using stress fields by numerical manifold method

This study first reviews the numerical manifold method(NMM)which possesses some advantages over the traditional limit equilibrium methods(LEMs)in calculating the factors of safety(Fs)of the slopes.Then,with regard to a trial slip surface(TSS),associated stress fields reproduced by NMM as well as the enhanced limit equilibrium method are combined to compute Fs.In order to search for the potential critical slip surface(CSS),the MAX-MIN ant colony optimization algorithm(MMACOA),one of the best performing algorithms for some optimization problems,is adopted.Procedures to obtain Fs in conjunction with the potential CSS are described.Finally,the proposed numerical model and traditional methods are compared with stability analysis of three typical slopes.The numerical results show that Fs and CSSs of the slopes can be accurately calculated with the proposed model.

A linearly-independent higher-order extended numerical manifold method and its application to multiple crack growth simulation

The numerical manifold method(NMM)can be viewed as an inherent continuous-discontinuous numerical method,which is based on two cover systems including mathematical and physical covers.Higher-order NMM that adopts higher-order polynomials as its local approximations generally shows higher precision than zero-order NMM whose local approximations are constants.Therefore,higherorder NMM will be an excellent choice for crack propagation problem which requires higher stress accuracy.In addition,it is crucial to improve the stress accuracy around the crack tip for determining the direction of crack growth according to the maximum circumferential stress criterion in fracture mechanics.Thus,some other enriched local approximations are introduced to model the stress singularity at the crack tip.Generally,higher-order NMM,especially first-order NMM wherein local approximations are first-order polynomials,has the linear dependence problems as other partition of unit(PUM)based numerical methods does.To overcome this problem,an extended NMM is developed based on a new local approximation derived from the triangular plate element in the finite element method(FEM),which has no linear dependence issue.Meanwhile,the stresses at the nodes of mathematical mesh(the nodal stresses in FEM)are continuous and the degrees of freedom defined on the physical patches are physically meaningful.Next,the extended NMM is employed to solve multiple crack propagation problems.It shows that the fracture mechanics requirement and mechanical equilibrium can be satisfied by the trial-and-error method and the adjustment of the load multiplier in the process of crack propagation.Four numerical examples are illustrated to verify the feasibility of the proposed extended NMM.The numerical examples indicate that the crack growths simulated by the extended NMM are in good accordance with the reference solutions.Thus the effectiveness and correctness of the developed NMM have been validated.

An overview of particle-based numerical manifold method and its application to dynamic rock fracturing

This review summarizes the development of particle-based numerical manifold method(PNMM)and its applications to rock dynamics.The fundamental principle of numerical manifold method(NMM)is first briefly introduced.Then,the history of the newly developed PNMM is given.Basic idea of PNMM and its simulation procedure are presented.Considering that PNMM could be regarded as an NMM-based model,a comparison of PNMM and NMM is discussed from several points of view in this paper.Besides,accomplished applications of PNMM to the dynamic rock fracturing are also reviewed.Finally,some recommendations are provided for the future work of PNMM.

Numerical Method Based on Compatible Manifold Element for Thin Plate Bending

The typical quadrangular and triangular elements for thin plate bending based on Kirchhoff assumptions are the non- conforming elements with low computational accuracy and limitative application range in fmite element method（FEM）. Some compatible elements can be developed by the means of supplementing correction functions, increasing nodes in element or on the boundaries, expanding nodal degrees of freedom（DOF）, etc, but these elements are inconvenient to apply in practice for the high calculation complexity. In this paper, in order to overcome the defects of thin plate bending finite element, numerical manifold method（NMM） was introduced to solve thin plate bending deformation problem. Rectangular mesh was adopted as mathematical mesh to form f＇mite element cover system, and then 16-cover manifold element was proposed. Numerical manifold formulas were constructed on the basis of minimum potential energy principle, displacement boundary conditions are implemented by penalty function method, and all the element matrixes were derived in details. The 16-cover element has a simple calculation process for employing only the transverse displacement cover DOFs as the basic unknown variables, and has been proved to meet the requirements of completeness and full compatibility. As an application, the presented 16-cover element has been used to analyze bending deformation of square thin plate under different loads and boundary conditions, and the results show that numerical manifold method with compatible element, compared with finite element method, can improve computational accuracy and convergence greatly.

A generalized cover renewal strategy for multiple crack propagation in two-dimensional numerical manifold method

Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references.In this study,a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems.The proposed algorithm updates the cover system with a bottom-up process:1)identification of fractured manifold elements according to the previous and latest crack tip position;and 2)local topological update of the manifold elements,physical patches,block boundary loops,and non-persistent joint loops according to the scenario classification of the propagating crack.The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis.Three crack propagation examples show that the proposed algorithm performs well in updating the cover system.This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.

Modeling rock failure using the numerical manifold method followed bythe discontinuous deformation analysis

A complete rock failure process usually involves opening/sliding of preexisting discontinuities as well as frac- turing in intact rock bridges to form persistent failure sur- faces and subsequent motions of the generated rock blocks. The recently developed numerical manifold method （NMM） has potential for modelling such a complete failure process. However, the NMM suffers one limitation, i.e., unexpected material domain area change occurs in rotation modelling. This problem can not be easily solved because the rigid body rotation is not represented explicitly in the NMM. The discontinuous deformation analysis （DDA） is specially de- veloped for modelling discrete block systems. The rotation- induced material area change in the DDA modelling can be avoided conveniently because the rigid body rotation is represented in an explicit form. In this paper, a transition technique is proposed and implemented to convert a NMMmodelling to a DDA modelling so as to simulate a complete rock failure process entirely by means of the two methods, in which the NMM is adopted to model the early fracturing as well as the transition from continua to discontinua, while the DDA is adopted to model the subsequent motion of the generated rock blocks. Such a numerical approach also im- proves the simulation efficiency greatly as compared with a complete NMM modelling approach. The fracturing of a rock slab with pre-existing non-persistent joints located on a slope crest and the induced rockfall process are simulated. The validity of the modelling transition from the NMM to the DDA is verified and the applicability of the proposed nu- merical approach is investigated.

THE VARIATIONAL PRINCIPLES AND APPLICATION OF NONLINEAR NUMERICAL MANIFOLD METHOD

The physical-cover-oriented variational principle of nonlinear numerical manifold method (NNMM) for the analysis of plastical problems is put forward according to the displacement model and the characters of numerical manifold method (NMM). The theoretical calculating formulations and the controlling equation of NNMM are derived. As an example, the plate with a hole in the center is calculated and the results show that the solution precision and efficiency of NNMM are agreeable.

Incompatible numerical manifold method for fracture problems

The incompatible numerical manifold method （INMM） is based on the finite cover approximation theory, which provides a unified framework for problems dealing with continuum and discontinuities. The incompatible numerical manifold method employs two cover systems as follows. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the weighted functions, and the physical cover system describes geometry of the domain and the discontinuous surfaces therein. In INMM, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the process of displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the INMM is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the INMM, the analytical solution is used at the tip of a crack, and thus the cover displacement functions are extended with higher precision and computational efficiency. Some numerical examples are given.

THEORETICAL STUDY OF THREE-DIMENSIONAL NUMERICAL MANIFOLD METHOD

The three-dimensional numerical manifold method（NMM） is studied on the basis of two-dimensional numerical manifold method. The three-dimensional cover displacement function is studied. The mechanical analysis and Hammer integral method of three-dimensional numerical manifold method are put forward. The stiffness matrix of three-dimensional manifold element is derived and the dissection rules are given. The theoretical system and the numerical realizing method of three-dimensional numerical manifold method are systematically studied. As an example, the cantilever with load on the end is calculated, and the results show that the precision and efficiency are agreeable.

THE VARIATIONAL PRINCIPLE AND APPLICATION OF NUMERICAL MANIFOLD METHOD

The physical-cover-oriented variational principle of numerical manifold method (NMM) for the analysis of linear elastic static problems was put forward according to the displacement model and the characters of numerical manifold method. ne theoretical calculating formulations and the controlling equation of NMM were derived. As an example, the plate with a hole in the center is calculated and the results show that the solution precision and efficiency of NMM are agreeable.

Numerical Manifold Method and Its Application in the Study of Crustal Movements in the Sichuan-Yunnan Area

The numerical manifold method （NMM） can calculate the movements and deformations of structures or materials. Both the finite element method （FEM） for continua and the discontinuous deformation analysis （DDA） for block systems are special cases of NMM. NMM has separate mathematical covers and physical meshes： the mathematical covers define only fine or rough approximations; as the real material boundary, the physical mesh defines the integration fields. The mathematical covers are triangle units; the physical mesh includes the fault boundaries, joints, blocks and interfaces of different crust zones on the basis of a geological tectonic background. Aiming at the complex problem of continuous and discontinuous deformation across the Chinese continent, the numerical manifold method （NMM） is brought in to study crustal movement of the Stchuan-Yunnan area. Based on the GPS velocity field in the Sichuan-Yunnan area, a crustal strain and stress field is simulated and analyzed. Moreover, results show that the NMM is a more suitable method than DDA in simulating the movement of the Sichuan-Yunnan area. Finally, a kind of mechanism of crustal motion in the Sichuan-Yunnan area is discussed in the paper.

Numerical simulation of crack propagation by manifold method

In this paper, by means of the maximum circle tensile stress on curve of constant ω and stress intensity factors by a path independent contour integral method, trajectories of maxed mode crack propagation are simulated through numerical manifold method. The crack propagation is traced dynamically by modifying the neighboring connection between the crack-top and nodes within elements in the calculating process. This method has the advantages such as less modified area, easiness of programming, high realizability and so on. Then a single sharp nicked specimen is used to verified the numerical result. It is shown that the provided method is reasonable and effective.

NUMERICAL MANIFOLD METHOD AND ITS APPLICATION IN UNDERGROUND OPENINGS

A brief introduction is made for the Numerical Manifold Method and its analysing process in rock mechanics. Some aspects of the manifold method are improved in implementing process according to the practice of excavating underground openings. Corresponding formulas are given and a computer program of the Numerical Manifold Method has been completed in this paper.

耦合渗流-变形的数值流形法裂隙岩质 边坡稳定性分析

岩石裂隙渗流问题作为影响岩石边坡失稳的重要因素吸引了越来越多研究者采用不同的数值方法模拟其力学行为。提出了一种基于渗流-变形耦合模型的数值流形方法,并相应地讨论了岩体裂隙网络中地下水的流动、渗透压力和岩石变形的耦合效应。基于最小能量原理,分别推导渗流条件下单元边界流体压力和系统的整体平衡方程,以及块体接触算法求解任意岩石裂缝的局部安全系数。通过一系列验证实例对所提出的数模模型进行检验,验证了该数值模型的鲁棒性和有效性。利用该方法模拟了一个岩石边坡由于渗流作用而发生的坍塌事故,再现了岩石边坡的破坏过程。模拟结果表明,过大的水力压力导致垂直裂缝张开,增加岩体变形量,最终导致边坡破坏。研究表明,该方法在模拟大规模工程问题方面具有很大的潜力。

酶底物法测定污泥中粪大肠菌群的方法研究及比对

为了提高污泥微生物学指标粪大肠菌群的检测效率,验证采用酶底物法检测该指标的可行性及适应性,开展方法研究及比对试验。采用酶底物法和多管发酵法对污泥中粪大肠菌群进行同步测定,检测结果采用SPSS软件进行数据处理及配对t检验,全国6家实验室201组比对结果显示:两组数据P=0.141(>0.05),相关系数r=0.964,表明两种检测方法相关性较好,没有统计学意义上的显著性差异,具有一致性。3家实验室的6个实验组分别对低、中、高三个浓度的实际样品以及有证标准物质配制的样品进行了6次重复测定:实验室内相对标准偏差范围分别为4.5%~13.7%、1.6%~6.7%、1.4%~3.6%和1.1%~4.5%;实验室间相对标准偏差分别为8.7%、5.2%、5.3%、3.1%;重复性限依次为0.44、0.39、0.43和0.27;再现性限依次为0.44、0.49、0.81和0.31;实验室内的相对误差范围为0.91%~2.8%;相对误差最终值为1.9%±1.5%。通过试验验证,酶底物法测定污泥中粪大肠菌群的方法精密度、准确度符合要求,确认方法检出限为10MPN/g。多管发酵法、滤膜法分别是现行检测污泥中粪大肠菌群的国标方法和行标方法,酶底物法与前两者相比,具有操作简单、准确度高、不需要确证试验、检测效率高等特点。酶底物法可作为测定污泥中粪大肠菌群指标的补充方法。

Fundamentals and progress of the manifold method based on independent covers

Aiming to solve mesh generation,computational stability,accuracy control,and other problems encountered with existing numerical methods,such as the finite element method and the finite volume method,a new numerical computational method for continuum mechanics,namely the manifold method based on independent covers(MMIC),is proposed based on the concept of mathematical manifolds,to form partitioned series solutions of partial differential equations.As partitions,the cover meshes have the characteristics of arbitrary shape,arbitrary connection,and arbitrary refinement.They are expected to fundamentally solve the mesh generation problem and can also simulate the precise geometric boundaries of the CAD model and strictly impose boundary conditions.In the selection of series solutions,local analytical solutions(such as series solutions at crack tips and series solutions in infinite domains)or proper forms of complete series can be used to reflect the local or global characteristics of the physical field to accelerate convergence.Various applications are presented.A new method of beam,plate,and shell analysis is proposed.The deformation characteristics of beams,plates,and shells are simulated with polynomial series of suitable forms,and the analysis of curved beams and shells with accurate geometric representation is realized.For the static elastic analysis of two-dimensional structures,a mesh splitting algorithm is proposed,and h-p version adaptive analysis is carried out with error estimation.Thus,automatic computation integrated with CAD is attempted.Adaptive analysis is also attempted for the solution of differential equations of fluids.For the one-dimensional convection-diffusion equation and Burgers equation,calculation results with high precision are obtained in strong convection and shock wave simulations,avoiding nonphysical oscillations.And solving the two-dimensional incompressible Navier-Stokes equation is also attempted.The series solution formula is used to obtain the physical quantity of interest of the material at a space point to eliminate the convection terms.Thus,geometrically nonlinear problems can be analyzed in fixed meshes,and a new method of free surface tracking is proposed.

Three-dimensional numerical manifold method for heat conduction problems with a simplex integral on the boundary

The three-dimensional numerical manifold method(3D-NMM),which is based on the derivation of Galerkin's variation,is a powerful calculation tool that uses two cover systems.The 3D-NMM can be used to handle continue-discontinue problems and extend to THM coupling.In this study,we extended the 3D-NMM to simulate both steady-state and transient heat conduction problems.The modelling was carried out using the raster methods(RSM).For the system equation,a variational method was employed to drive the discrete equations,and the crucial boundary conditions were solved using the penalty method.To solve the boundary integral problem,the face integral of scalar fields and two-dimensional simplex integration were used to accurately describe the integral on polygonal boundaries.Several numerical examples were used to verify the results of 3D steady-state and transient heat-conduction problems.The numerical results indicated that the 3D-NMM is effective for handling 3D both steadystate and transient heat conduction problems with high solution accuracy.