The three-dimensional numerical manifold method(3DNMM) method is further enriched to simulate wave propagation across homogeneous/jointed rock masses. For the purpose of minimizing negative effects from artificial bou...The three-dimensional numerical manifold method(3DNMM) method is further enriched to simulate wave propagation across homogeneous/jointed rock masses. For the purpose of minimizing negative effects from artificial boundaries, a viscous nonreflecting boundary, which can effectively absorb the energy of a wave, is firstly adopted to enrich 3DNMM. Then, to simulate the elastic recovery property of an infinite problem domain, a viscoelastic boundary, which is developed from the viscous nonreflecting boundary, is further adopted to enrich 3DNMM. Finally, to eliminate the noise caused by scattered waves, a force input method which can input the incident wave correctly is incorporated into 3DNMM. Five typical numerical tests on P/S-wave propagation across jointed/homogeneous rock masses are conducted to validate the enriched 3DNMM. Numerical results indicate that wave propagation problems within homogeneous and jointed rock masses can be correctly and reliably modeled with the enriched 3DNMM.展开更多
We present a numerically stable one-point quadrature rule for the stiffness matrix and mass matrix of the three-dimensional numerical manifold method(3D NMM).The rule simplifies the integration over irregularly shaped...We present a numerically stable one-point quadrature rule for the stiffness matrix and mass matrix of the three-dimensional numerical manifold method(3D NMM).The rule simplifies the integration over irregularly shaped manifold elements and overcomes locking issues,and it does not cause spurious modes in modal analysis.The essential idea is to transfer the integral over a manifold element to a few moments to the element center,thereby deriving a one-point integration rule by the moments and making modifications to avoid locking issues.For the stiffness matrix,after the virtual work is decomposed into moments,higher-order moments are modified to overcome locking issues in nearly incompressible and bending-dominated conditions.For the mass matrix,the consistent and lumped types are derived by moments.In particular,the lumped type has the clear advantage of simplicity.The proposed method is naturally suitable for 3D NMM meshes automatically generated from a regular grid.Numerical tests justify the accuracy improvements and the stability of the proposed procedure.展开更多
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.展开更多
基金supported by the Youth Innovation Promotion Association CAS(Grant No. 2020327)the National Natural Science Foundation of China(Grant Nos. 12202024, 52130905, 12272393, and 12072357)。
文摘The three-dimensional numerical manifold method(3DNMM) method is further enriched to simulate wave propagation across homogeneous/jointed rock masses. For the purpose of minimizing negative effects from artificial boundaries, a viscous nonreflecting boundary, which can effectively absorb the energy of a wave, is firstly adopted to enrich 3DNMM. Then, to simulate the elastic recovery property of an infinite problem domain, a viscoelastic boundary, which is developed from the viscous nonreflecting boundary, is further adopted to enrich 3DNMM. Finally, to eliminate the noise caused by scattered waves, a force input method which can input the incident wave correctly is incorporated into 3DNMM. Five typical numerical tests on P/S-wave propagation across jointed/homogeneous rock masses are conducted to validate the enriched 3DNMM. Numerical results indicate that wave propagation problems within homogeneous and jointed rock masses can be correctly and reliably modeled with the enriched 3DNMM.
基金supported by the National Natural Science Foundation of China(Grant Nos.42302331,52130905 and 52079002)。
文摘We present a numerically stable one-point quadrature rule for the stiffness matrix and mass matrix of the three-dimensional numerical manifold method(3D NMM).The rule simplifies the integration over irregularly shaped manifold elements and overcomes locking issues,and it does not cause spurious modes in modal analysis.The essential idea is to transfer the integral over a manifold element to a few moments to the element center,thereby deriving a one-point integration rule by the moments and making modifications to avoid locking issues.For the stiffness matrix,after the virtual work is decomposed into moments,higher-order moments are modified to overcome locking issues in nearly incompressible and bending-dominated conditions.For the mass matrix,the consistent and lumped types are derived by moments.In particular,the lumped type has the clear advantage of simplicity.The proposed method is naturally suitable for 3D NMM meshes automatically generated from a regular grid.Numerical tests justify the accuracy improvements and the stability of the proposed procedure.
基金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.