The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock mater...The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.展开更多
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 consis...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.展开更多
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 sheari...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.展开更多
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 ...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.展开更多
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 accura...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.展开更多
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 t...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.展开更多
There are relatively few studies on large rotation or deformation by means of the three-dimensional(3D)numerical manifold method(NMM).A new modified symmetric and antisymmetric decomposition(MSAD)theory is developed a...There are relatively few studies on large rotation or deformation by means of the three-dimensional(3D)numerical manifold method(NMM).A new modified symmetric and antisymmetric decomposition(MSAD)theory is developed and implemented into the 3D NMM,eliminating the false-volume expansion and false-rotation strain/stress problems.The Jaumann rate is used to measure the material rotation,and the geometric stiffness built on the Jaumann rate is deduced.The incremental formulas of the MSAD-based 3D NMM and a practical guide on the implementation of the MSAD theory are given in detail and exemplified.The new theory and formulas can be applied to analyze both large rotation and large deformation problems.Based on the hypoelasto-plasticity theory and the unified strength theory,the unified yield criterion with associated flow rule is implemented into the MSAD-based 3D NMM.Several typical examples are studied,showing the advantage and potential of the new MSAD theory and the MSAD-based 3D NMM.展开更多
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 genera...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 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...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.展开更多
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 numer...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.展开更多
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...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.展开更多
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 Ha...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 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 numeric...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.展开更多
Three-dimensional numerical manifold method for unconfined seepage analysis is proposed in this article.By constructing hydraulic potential functions of the manifold element,the element conductivity matrix and the glo...Three-dimensional numerical manifold method for unconfined seepage analysis is proposed in this article.By constructing hydraulic potential functions of the manifold element,the element conductivity matrix and the global simultaneous equations for unconfined seepage analysis are derived in detail.The algorithm of locating the free surface and the formula for seepage forces are also given.Three-dimensional manifold method employs the tetrahedral mathematical meshes to cover the whole material volume.In the iterative process for locating the free surface,the manifold method can achieve an accurate seepage analysis of the saturated domain below the free surface with mathematical meshes unchanged.Since the shape of manifold elements can be arbitrary,the disadvantage of changing the permeability of transitional elements cut by the free surface in the conventional Finite Element Method（FEM） is removed,and the accuracy of locating the free surface can be ensured.Furthermore,the seepage force acting on the transitional elements can be accurately calculated by the simplex integration.Numerical results for a typical example demonstrate the validity of the proposed method.展开更多
A 2nd order numerical manifold method(NMM) based method is developed to simulate the hydraulic fractures propagating process in rock or concrete. The proposed method uses a weak coupling technique to analyze the fluid...A 2nd order numerical manifold method(NMM) based method is developed to simulate the hydraulic fractures propagating process in rock or concrete. The proposed method uses a weak coupling technique to analyze the fluid phase and solid phase. To study the seepage behavior of the fluid phase, all the fractures in solid are identified by a block cutting algorithm and form a flow network. Then the hydraulic heads at crack ends are solved. To study the deformation and destruction of solid phase, the 2-order NMM and sub-region boundary element method are combined to solve the stress-strain field. Crack growth is controlled by the well-accepted criterion, including the tension criterion or Mohr-Coulomb criterion for the initialization of cracks and the maximum circumferential stress theory for crack propagation. Once the crack growth occurs, the seepage and deformation analysis will be resolved in the next simulation step. Such weak coupling analysis will continue until the structure becomes stable or is destructed. Five examples are used to verify the new method. The results demonstrate that the method can solve the SIFs at crack tip and fluid flow in crack network precisely, and the method is effective in simulating the hydraulic facture problem. Besides, the NMM shows great convenience and is of high accuracy in simulating the crack growth problem.展开更多
In order to reach the best numerical properties with the numerical manifold method(NMM),uniform finite element meshes are always favorite while constructing mathematical covers,where all the elements are congruent.In ...In order to reach the best numerical properties with the numerical manifold method(NMM),uniform finite element meshes are always favorite while constructing mathematical covers,where all the elements are congruent.In the presence of steep gradients or strong singularities,in principle,the locally-defined special functions can be added into the NMM space by means of the partition of unity,but they are not available to those complex problems with heterogeneity or nonlinearity,necessitating local refinement on uniform meshes.This is believed to be one of the most important open issues in NMM.In this study multilayer covers are proposed to solve this issue.In addition to the first layer cover which is the global cover and covers the whole problem domain,the second and higher layer covers with smaller elements,called local covers,are used to cover those local regions with steep gradients or strong singularities.The global cover and the local covers have their own partition of unity,and they all participate in the approximation to the solution.Being advantageous over the existing procedures,the proposed approach is easy to deal with any arbitrary-layer hanging nodes with no need to construct super-elements with variable number of edge nodes or introduce the Lagrange multipliers to enforce the continuity between small and big elements.With no limitation to cover layers,meanwhile,the creation of an even error distribution over the whole problem domain is significantly facilitated.Some typical examples with steep gradients or strong singularities are analyzed to demonstrate the capacity of the proposed approach.展开更多
A three-node triangular element fitted to numerical manifold method with continuous nodal stress, called Trig_3-CNS(NMM)element, was recently proposed for linear elastic continuous problems and linear elastic simple c...A three-node triangular element fitted to numerical manifold method with continuous nodal stress, called Trig_3-CNS(NMM)element, was recently proposed for linear elastic continuous problems and linear elastic simple crack problems. The Trig_3-CNS(NMM) element can be considered as a development of both the Trig_3-CNS element and the numerical manifold method(NMM).Inheriting all the advantages of Trig_3-CNS element, calculations using Trig_3-CNS(NMM) element can obtain higher accuracy than Trig_3 element without extra degrees of freedom(DOFs) and yield continuous nodal stress without stress smoothing. Inheriting all the advantages of NMM, Trig_3-CNS(NMM) element can conveniently treat crack problems without deploying conforming mathematical mesh. In this paper,complex problems such as a crucifix crack and a star-shaped crack with many branches are studied to exhibit the advantageous features of the Trig_3-CNS(NMM) element. Numerical results show that the Trig_3-CNS(NMM) element is prominent in modeling complex crack problems.展开更多
In this study, numerical manifold method(NMM) coupled with non-uniform rational B-splines(NURBS) and T-splines in the context of isogeometric analysis is proposed to allow for the treatments of complex geometries and ...In this study, numerical manifold method(NMM) coupled with non-uniform rational B-splines(NURBS) and T-splines in the context of isogeometric analysis is proposed to allow for the treatments of complex geometries and local refinement. Computational formula for a 9-node NMM based on quadratic B-splines is derived. In order to exactly represent some common free-form shapes such as circles, arcs, and ellipsoids, quadratic non-uniform rational B-splines(NURBS) are introduced into NMM. The coordinate transformation based on the basis function of NURBS is established to enable exact integration for the manifold elements containing those shapes. For the case of crack propagation problems where singular fields around crack tips exist, local refinement technique by the application of T-spline discretizations is incorporated into NMM, which facilitates a truly local refinement without extending the entire row of control points. A local refinement strategy for the 4-node mathematical cover mesh based on T-splines and Lagrange interpolation polynomial is proposed. Results from numerical examples show that the 9-node NMM based on NURBS has higher accuracies. The coordinate transformation based on the NURBS basis function improves the accuracy of NMM by exact integration. The local mesh refinement using T-splines reduces the number of degrees of freedom while maintaining calculation accuracy at the same time.展开更多
In this paper,strategies are provided for a powerful numerical manifold method(NMM)with h and p refinement in analyses of elasticity and elasto-plasticity.The new NMM is based on the concept of the independent cover,w...In this paper,strategies are provided for a powerful numerical manifold method(NMM)with h and p refinement in analyses of elasticity and elasto-plasticity.The new NMM is based on the concept of the independent cover,which gets rid of NMM's important defect of rank deficiency when using higher-order local approximation functions.Several techniques are presented.In terms of mesh generation,a relationship between the quadtree structure and the mathematical mesh is established to allow a robust h-refinement.As to the condition number,a scaling based on the physical patch is much better than the classical scaling based on the mathematical patch;an overlapping width of 1%–10%can ensure a good condition number for 2nd,3rd,and 4th order local approximation functions;the small element issue can be overcome after the local approximation on small patch is replaced by that on a regular patch.On numerical accuracy,local approximation using complete polynomials is necessary for the optimal convergence rate.Two issues that may damage the convergence rate should be prevented.The first is to approximate the curved boundary of a higher-order element by overly few straight lines,and the second is excessive overlapping width.Finally,several refinement strategies are verified by numerical examples.展开更多
Numerical simulations of longitudinal wave propagation in a rock bar with microcracks are conducted by using the numerical manifold method which has great advantages in the simulation of discontinuities.Firstly,valida...Numerical simulations of longitudinal wave propagation in a rock bar with microcracks are conducted by using the numerical manifold method which has great advantages in the simulation of discontinuities.Firstly,validation of the numerical manifold method is carried out by simulations of a longitudinal stress wave propagating through intact and cracked rock bars.The behavior of the stress wave traveling in a one-dimensional rock bar with randomly distributed microcracks is subsequently studied.It is revealed that the highly defected rock bar has significant viscoelasticity to the stress wave propagation.Wave attenuation as well as time delay is affected by the length,quantity,specific stiffness of the distributed microcracks as well as the incident stress wave frequency.The storage and loss moduli of the defected rock are also affected by the microcrack properties;however,they are independent of incident stress wave frequency.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.42172312 and 52211540395)support from the Institut Universitaire de France(IUF).
文摘The modeling of crack growth in three-dimensional(3D)space poses significant challenges in rock mechanics due to the complex numerical computation involved in simulating crack propagation and interaction in rock materials.In this study,we present a novel approach that introduces a 3D numerical manifold method(3D-NMM)with a geometric kernel to enhance computational efficiency.Specifically,the maximum tensile stress criterion is adopted as a crack growth criterion to achieve strong discontinuous crack growth,and a local crack tracking algorithm and an angle correction technique are incorporated to address minor limitations of the algorithm in a 3D model.The implementation of the program is carried out in Python,using object-oriented programming in two independent modules:a calculation module and a crack module.Furthermore,we propose feasible improvements to enhance the performance of the algorithm.Finally,we demonstrate the feasibility and effectiveness of the enhanced algorithm in the 3D-NMM using four numerical examples.This study establishes the potential of the 3DNMM,combined with the local tracking algorithm,for accurately modeling 3D crack propagation in brittle rock materials.
基金supported by the Fund of National Engineering and Research Center for Highways in Mountain Area(No.gsgzj-2012-05)the Fundamental Research Funds for the Central Universities of China(No.CDJXS12240003)the Scientific Research Foundation of State Key Laboratory of Coal Mine Disaster Dynamics and Control(No.2011DA105287-MS201213)
文摘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.
基金supported by Laboratory Directed Research and Development(LDRD)funding from Berkeley Labsupported by Open Fund of the State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z017004)。
文摘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.
基金This study is supported by the Youth Innovation Promotion Association of Chinese Academy of Sciences(Grant No.2020327)the National Natural Science Foundation of China(Grant No.51609240).
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘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.
基金Project(51321065,51479191,11672360)supported by the National Natural Science Foundation of China。
文摘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.
基金This research is supported by the National Basic Research Program of China(973 Program,Grant No.2014CB047100)the National Natural Science Foundation of China(Grant Nos.41472289,51179185 and 41807275).
文摘There are relatively few studies on large rotation or deformation by means of the three-dimensional(3D)numerical manifold method(NMM).A new modified symmetric and antisymmetric decomposition(MSAD)theory is developed and implemented into the 3D NMM,eliminating the false-volume expansion and false-rotation strain/stress problems.The Jaumann rate is used to measure the material rotation,and the geometric stiffness built on the Jaumann rate is deduced.The incremental formulas of the MSAD-based 3D NMM and a practical guide on the implementation of the MSAD theory are given in detail and exemplified.The new theory and formulas can be applied to analyze both large rotation and large deformation problems.Based on the hypoelasto-plasticity theory and the unified strength theory,the unified yield criterion with associated flow rule is implemented into the MSAD-based 3D NMM.Several typical examples are studied,showing the advantage and potential of the new MSAD theory and the MSAD-based 3D NMM.
基金supported by the Research Fund for the Doctoral Program of Higher Education of China (20090101120057)the Scientific Research Fund of Zhejiang Provincial Education Department (Y200909163)
文摘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.
基金Supported by the National Natural Science Foundation of China （N0.40574006, N0.40344023）, DGLIGG （L04-02）.
文摘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.
基金supported by the Natural Science Foundation of Shandong Province for Excellent Young and Middle-aged Scientist (2007BS04045 and 2008BS04009)the Natural Science Foundation of Shandong Province(Y2006B24 and Y2008A 11)
文摘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.
文摘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.
文摘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 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.
基金supported by the National Natural Science Foundation of China (Grant Nos. 50725931, 50839004)the Ministry of Education of China for New Century Excellent Talents in University (Grant No. NCET-07-0632)
文摘Three-dimensional numerical manifold method for unconfined seepage analysis is proposed in this article.By constructing hydraulic potential functions of the manifold element,the element conductivity matrix and the global simultaneous equations for unconfined seepage analysis are derived in detail.The algorithm of locating the free surface and the formula for seepage forces are also given.Three-dimensional manifold method employs the tetrahedral mathematical meshes to cover the whole material volume.In the iterative process for locating the free surface,the manifold method can achieve an accurate seepage analysis of the saturated domain below the free surface with mathematical meshes unchanged.Since the shape of manifold elements can be arbitrary,the disadvantage of changing the permeability of transitional elements cut by the free surface in the conventional Finite Element Method（FEM） is removed,and the accuracy of locating the free surface can be ensured.Furthermore,the seepage force acting on the transitional elements can be accurately calculated by the simplex integration.Numerical results for a typical example demonstrate the validity of the proposed method.
基金supported by the National Natural Science Foundation of China(Grant Nos.51439005&51209235)the National Basic Research Program of China("973"Project)(Grant Nos.2013CB035904,2013CB-036406)
文摘A 2nd order numerical manifold method(NMM) based method is developed to simulate the hydraulic fractures propagating process in rock or concrete. The proposed method uses a weak coupling technique to analyze the fluid phase and solid phase. To study the seepage behavior of the fluid phase, all the fractures in solid are identified by a block cutting algorithm and form a flow network. Then the hydraulic heads at crack ends are solved. To study the deformation and destruction of solid phase, the 2-order NMM and sub-region boundary element method are combined to solve the stress-strain field. Crack growth is controlled by the well-accepted criterion, including the tension criterion or Mohr-Coulomb criterion for the initialization of cracks and the maximum circumferential stress theory for crack propagation. Once the crack growth occurs, the seepage and deformation analysis will be resolved in the next simulation step. Such weak coupling analysis will continue until the structure becomes stable or is destructed. Five examples are used to verify the new method. The results demonstrate that the method can solve the SIFs at crack tip and fluid flow in crack network precisely, and the method is effective in simulating the hydraulic facture problem. Besides, the NMM shows great convenience and is of high accuracy in simulating the crack growth problem.
基金supported by the National Basic Research Program of China("973"Project)(Grant Nos.2011CB013505&2014CB047100)the National Natural Science Foundation of China(Grant Nos.11572009&51538001)
文摘In order to reach the best numerical properties with the numerical manifold method(NMM),uniform finite element meshes are always favorite while constructing mathematical covers,where all the elements are congruent.In the presence of steep gradients or strong singularities,in principle,the locally-defined special functions can be added into the NMM space by means of the partition of unity,but they are not available to those complex problems with heterogeneity or nonlinearity,necessitating local refinement on uniform meshes.This is believed to be one of the most important open issues in NMM.In this study multilayer covers are proposed to solve this issue.In addition to the first layer cover which is the global cover and covers the whole problem domain,the second and higher layer covers with smaller elements,called local covers,are used to cover those local regions with steep gradients or strong singularities.The global cover and the local covers have their own partition of unity,and they all participate in the approximation to the solution.Being advantageous over the existing procedures,the proposed approach is easy to deal with any arbitrary-layer hanging nodes with no need to construct super-elements with variable number of edge nodes or introduce the Lagrange multipliers to enforce the continuity between small and big elements.With no limitation to cover layers,meanwhile,the creation of an even error distribution over the whole problem domain is significantly facilitated.Some typical examples with steep gradients or strong singularities are analyzed to demonstrate the capacity of the proposed approach.
基金the National Natural Science Foundation of China(Grant Nos 51609240,11572009&51538001)and the National Basic Research Program of China(Grant No 2014CB047100)
文摘A three-node triangular element fitted to numerical manifold method with continuous nodal stress, called Trig_3-CNS(NMM)element, was recently proposed for linear elastic continuous problems and linear elastic simple crack problems. The Trig_3-CNS(NMM) element can be considered as a development of both the Trig_3-CNS element and the numerical manifold method(NMM).Inheriting all the advantages of Trig_3-CNS element, calculations using Trig_3-CNS(NMM) element can obtain higher accuracy than Trig_3 element without extra degrees of freedom(DOFs) and yield continuous nodal stress without stress smoothing. Inheriting all the advantages of NMM, Trig_3-CNS(NMM) element can conveniently treat crack problems without deploying conforming mathematical mesh. In this paper,complex problems such as a crucifix crack and a star-shaped crack with many branches are studied to exhibit the advantageous features of the Trig_3-CNS(NMM) element. Numerical results show that the Trig_3-CNS(NMM) element is prominent in modeling complex crack problems.
基金supported by the National Basic Research Program of China("973"Project)(Grant No.2014CB047100)the National Natural Science Foundation of China(Grant No.41372316)
文摘In this study, numerical manifold method(NMM) coupled with non-uniform rational B-splines(NURBS) and T-splines in the context of isogeometric analysis is proposed to allow for the treatments of complex geometries and local refinement. Computational formula for a 9-node NMM based on quadratic B-splines is derived. In order to exactly represent some common free-form shapes such as circles, arcs, and ellipsoids, quadratic non-uniform rational B-splines(NURBS) are introduced into NMM. The coordinate transformation based on the basis function of NURBS is established to enable exact integration for the manifold elements containing those shapes. For the case of crack propagation problems where singular fields around crack tips exist, local refinement technique by the application of T-spline discretizations is incorporated into NMM, which facilitates a truly local refinement without extending the entire row of control points. A local refinement strategy for the 4-node mathematical cover mesh based on T-splines and Lagrange interpolation polynomial is proposed. Results from numerical examples show that the 9-node NMM based on NURBS has higher accuracies. The coordinate transformation based on the NURBS basis function improves the accuracy of NMM by exact integration. The local mesh refinement using T-splines reduces the number of degrees of freedom while maintaining calculation accuracy at the same time.
基金supported by the National Natural Science Foundation of China(Grant Nos.52130905 and 52079002)。
文摘In this paper,strategies are provided for a powerful numerical manifold method(NMM)with h and p refinement in analyses of elasticity and elasto-plasticity.The new NMM is based on the concept of the independent cover,which gets rid of NMM's important defect of rank deficiency when using higher-order local approximation functions.Several techniques are presented.In terms of mesh generation,a relationship between the quadtree structure and the mathematical mesh is established to allow a robust h-refinement.As to the condition number,a scaling based on the physical patch is much better than the classical scaling based on the mathematical patch;an overlapping width of 1%–10%can ensure a good condition number for 2nd,3rd,and 4th order local approximation functions;the small element issue can be overcome after the local approximation on small patch is replaced by that on a regular patch.On numerical accuracy,local approximation using complete polynomials is necessary for the optimal convergence rate.Two issues that may damage the convergence rate should be prevented.The first is to approximate the curved boundary of a higher-order element by overly few straight lines,and the second is excessive overlapping width.Finally,several refinement strategies are verified by numerical examples.
文摘Numerical simulations of longitudinal wave propagation in a rock bar with microcracks are conducted by using the numerical manifold method which has great advantages in the simulation of discontinuities.Firstly,validation of the numerical manifold method is carried out by simulations of a longitudinal stress wave propagating through intact and cracked rock bars.The behavior of the stress wave traveling in a one-dimensional rock bar with randomly distributed microcracks is subsequently studied.It is revealed that the highly defected rock bar has significant viscoelasticity to the stress wave propagation.Wave attenuation as well as time delay is affected by the length,quantity,specific stiffness of the distributed microcracks as well as the incident stress wave frequency.The storage and loss moduli of the defected rock are also affected by the microcrack properties;however,they are independent of incident stress wave frequency.
