To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with ...To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.展开更多
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.展开更多
基金the financial support from the National Natural Science Foundation of China(Grant No.51839003)Liaoning Revitalization Talents Program(Grant No.XLYCYSZX 1902)Hubei Key Laboratory for Efficient Utilization and Agglomeration of Metallurgic Mineral Resources(Grant No.2023zy002).
文摘To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.
基金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.