Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the d...Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the deformation properties of complex fractured rock mass.It comprehends both the advantages of the discrete fracture network model and the equivalent continuum model to capture the features of discontinuities explicitly while reducing computational intensity.The complex fracture network is stochastically split into a number of subfracture networks according to the domain,length or angle.An analytical solution is derived to infer theoretically the relationship between the elastic moduli of the original complex fractured rock mass and the split subfractured rock masses by introducing a correction term based on the deformation superposition principle.Numerical simulations are conducted to determine the elastic moduli of split subfractured rock masses using universal distinct element code(UDEC),while the elastic modulus of the original model is estimated based on the currently proposed analytical relationship.The results show that the estimation accuracy with the current domainbased splitting model is far superior compared to those with the other two splitting models.Thus,the estimation method of elastic modulus of complex fractured rock mass based on domain splitting mode of fracture network is identified as the multi-domain equivalent method proposed in this paper.The reliability of this method is evaluated,and its high computational efficiency is demonstrated through exemplification with regard to different geometric configurations for stochastically artificial discrete fracture network.The proposed multi-domain equivalent method constructs the theoretical framework except for the regression analysis hypothesis compared to the density-reduced model equivalent method.展开更多
To solve the problems of divergence,low accuracy and project application of membrane wrinkling analysis,an analysis method of zero shear modulus and equivalent stiffness was proposed.This method is an improvement to t...To solve the problems of divergence,low accuracy and project application of membrane wrinkling analysis,an analysis method of zero shear modulus and equivalent stiffness was proposed.This method is an improvement to the previous method (Method I) of local coordinate transposition and stiffness equivalence.The new method is derived and the feasibility is theoretically proved.A small-scale membrane structure is analyzed by the two methods,and the results show that the computational efficiency of the new method (Method II) is approximately 23 times that of Method I.When Method II is applied to a large-scale membrane stadium structure,it is found that this new method can quickly make the second principal stress of one way wrinkled elements zero,and make the two principal stresses of two-way wrinkled elements zero as well.It could attain the correct load responses right after the appearance of wrinkled elements,which indicates that Method II can be applied to wrinkling analysis of large-scale membrane structures.展开更多
Methods that can efficiently model the effects of rock joints on rock mass behavior can be beneficial in rock engineering. The suitability of equivalent rock mass(ERM) technique based upon particle methods is investig...Methods that can efficiently model the effects of rock joints on rock mass behavior can be beneficial in rock engineering. The suitability of equivalent rock mass(ERM) technique based upon particle methods is investigated. The ERM methodology is first validated by comparing calculated and experimental data of lab triaxial compression test on a set of cylindrical rock mass samples, each containing a single joint oriented in various dip angles. The simulated results are then used to study the stress-strain nonlinearity and failure mechanism as a function of the joint dip angle and confining stress. The anisotropy and size effects are also investigated by using multi-scale cubic ERM models subjected to triaxial compression test. The deformation and failure behavior are found to be influenced by joint degradation, the micro-crack formation in the intact rock, the interaction between two joints, and the interactions of micro-cracks and joints.展开更多
In the longitudinal seismic deformation method for shield tunnels,one of the most commonly used is the longitudinal equivalent stiffness beam model(LES)for simulating the mechanical behavior of the lining.In this mode...In the longitudinal seismic deformation method for shield tunnels,one of the most commonly used is the longitudinal equivalent stiffness beam model(LES)for simulating the mechanical behavior of the lining.In this model,axial deformation and bending deformation are independent,so the equivalent stiffness is a constant value.However,the actual situation is that axial deformation and bending deformation occur simultaneously,which is not considered in LES.At present,we are not clear about the effect on the calculation results when axial deformation and bending deformation occur simultaneously.Therefore,in this paper,we improve the traditional LES by taking the relative deformation as a load and considering the coordinated deformation of axial and bending degrees of freedom.This improved model is called DNLES,and its neutral axis equations are an explicit expression.Then,we propose an iterative algorithm to solve the calculation model of the DNLES-based longitudinal seismic deformation method.Through a calculation example,we find that the internal forces based on LES are notably underestimated than those of DNLES in the compression bending zone,while are overestimated in the tension bending zone.When considering the combined effect,the maximum bending moment reached 13.7 times that of the LES model,and the axial pressure and tension were about 1.14 and 0.96 times,respectively.Further analysis reveals the coordinated deformation process in the axial and bending directions of the shield tunnel,which leads to a consequent change in equivalent stiffness.This explains why,in the longitudinal seismic deformation method,the traditional LES may result in unreasonable calculation results.展开更多
An equivalent continuum method is developed to analyze the effective stiffness of three-dimensional stretching dominated lattice materials. The strength and three-dimensional plastic yield surfaces are calculated for ...An equivalent continuum method is developed to analyze the effective stiffness of three-dimensional stretching dominated lattice materials. The strength and three-dimensional plastic yield surfaces are calculated for the equivalent continuum. A yielding model is formulated and compared with the results of other models. The bedding-in effect is considered to include the compliance of the lattice joints. The predicted stiffness and strength are in good agreement with the experimental data, validating the present model in the prediction of the mechanical properties of stretching dominated lattice structures.展开更多
The rock mass consists of rock blocks and structural planes,which can reduce its integrity and strength.Therefore,accurately obtaining the characteristics of the rock mass structural plane is a prerequisite for evalua...The rock mass consists of rock blocks and structural planes,which can reduce its integrity and strength.Therefore,accurately obtaining the characteristics of the rock mass structural plane is a prerequisite for evaluating stability and designing supports in underground engineering.Currently,there are no effective testing methods for the characteristic parameters of the rock mass structural plane in underground engineering.The paper presents the digital drilling technology as a new testing method of rock mass structural planes.Flawed rock specimens with cracks of varying widths and angles were used to simulate the rock mass structural planes,and the multifunctional rock mass digital drilling test system was employed to carry out the digital drilling tests.The analysis focuses on the variation laws of drilling parameters,such as drilling pressure and drilling torque,affected by the characteristics of prefabricated cracks,and clarifies the degradation mechanism of rock equivalent compressive strength.Additionally,an identification model for the characteristic parameters of rock mass structural planes during drilling is established.The test results indicate that the average difference of the characteristics of prefabricated cracks identified by the equivalent compressive strength is 2.45°and 0.82 mm,respectively.The identification model while drilling is verified to be correct due to the high identification accuracy.Based on this,a method for testing the characteristic parameters of the surrounding rock structural plane while drilling is proposed.The research offers a theoretical and methodological foundation for precise in situ identification of structural planes of the surrounding rock in underground engineering.展开更多
This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock m...This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.展开更多
Quantitative damage identification of surrounding rock is important to assess the current condition and residual strength of underground tunnels.In this work,an underground tunnel model with marble-like cementitious m...Quantitative damage identification of surrounding rock is important to assess the current condition and residual strength of underground tunnels.In this work,an underground tunnel model with marble-like cementitious materials was first fabricated using the three-dimensional(3D)printing technique and then loaded to simulate its failure mode in the laboratory.Lead zirconate titanate piezoelectric(PZT)transducers were embedded in the surrounding rock around the tunnel in the process of 3D printing.A 3D monitoring network was formed to locate damage areas and evaluate damage extent during loading.Results show that as the load increased,main cracks firstly appeared above the tunnel roof and below the floor,and then they coalesced into the tunnel boundary.Finally,the tunnel model was broken into several parts.The resonant frequency and the peak of the conductance signature firstly shifted rightwards with loading due to the sealing of microcracks,and then shifted backwards after new cracks appeared.An overall increase in the root-mean-square deviation(RMSD)calculated from conductance signatures of all the PZT transducers was observed as the load(damage)increased.Damage-dependent equivalent stiffness parameters(ESPs)were calculated from the real and imaginary signatures of each PZT at different damage states.Satisfactory agreement between equivalent and experimental ESP values was achieved.Also,the relationship between the change of the ESP and the residual strength was obtained.The method paves the way for damage identification and residual strength estimation of other 3D printed structures in civil engineering.展开更多
This paper proposes a new technique which introduces the high-order single-step-β method(HSM)into the experimental study on the substructure pseudo-dynamic testing(SPDT).The technique is based on the proposed concept...This paper proposes a new technique which introduces the high-order single-step-β method(HSM)into the experimental study on the substructure pseudo-dynamic testing(SPDT).The technique is based on the proposed concept of equivalent shear stiffness which can meet the requirement of the HSM algorithm.A study is done to theoretically validate the technique by the numerical analysis of two-storey shear building structure,in comparison of the proposed substructure pseudo-dynamic testing algorithm with the central difference method(CDM).Then,a full-scale SPDT model,the three-storey frame-supported reinforced concrete short-limb masonry shear wall structure,is designed and tested to simulate the seismic response of the corresponding six-storey structure and verify the proposed force control HSM technique.Meanwhile,the techniques of both stiffness correction and force control are suggested to control algorithmic error,control error and measurement error.The results indicate that the force control HSM can be used in the full-scale multi-degree-of-freedom(MDOF)substructure pseudo-dynamic testing before descent segment of structure restoring force properties.展开更多
基金financial support by the National Natural Science Foundation of China(Grant Nos.52008152,U1965204,52061160367,U2067203 and 52008153)Natural Science Foundation of Hebei Province of China(Grant No.E2021202087)Hebei Department of Human Resource(Grant No.E2020050015)。
文摘Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the deformation properties of complex fractured rock mass.It comprehends both the advantages of the discrete fracture network model and the equivalent continuum model to capture the features of discontinuities explicitly while reducing computational intensity.The complex fracture network is stochastically split into a number of subfracture networks according to the domain,length or angle.An analytical solution is derived to infer theoretically the relationship between the elastic moduli of the original complex fractured rock mass and the split subfractured rock masses by introducing a correction term based on the deformation superposition principle.Numerical simulations are conducted to determine the elastic moduli of split subfractured rock masses using universal distinct element code(UDEC),while the elastic modulus of the original model is estimated based on the currently proposed analytical relationship.The results show that the estimation accuracy with the current domainbased splitting model is far superior compared to those with the other two splitting models.Thus,the estimation method of elastic modulus of complex fractured rock mass based on domain splitting mode of fracture network is identified as the multi-domain equivalent method proposed in this paper.The reliability of this method is evaluated,and its high computational efficiency is demonstrated through exemplification with regard to different geometric configurations for stochastically artificial discrete fracture network.The proposed multi-domain equivalent method constructs the theoretical framework except for the regression analysis hypothesis compared to the density-reduced model equivalent method.
基金Project(020940) supported by the Natural Science Foundation of Guangdong Province,China
文摘To solve the problems of divergence,low accuracy and project application of membrane wrinkling analysis,an analysis method of zero shear modulus and equivalent stiffness was proposed.This method is an improvement to the previous method (Method I) of local coordinate transposition and stiffness equivalence.The new method is derived and the feasibility is theoretically proved.A small-scale membrane structure is analyzed by the two methods,and the results show that the computational efficiency of the new method (Method II) is approximately 23 times that of Method I.When Method II is applied to a large-scale membrane stadium structure,it is found that this new method can quickly make the second principal stress of one way wrinkled elements zero,and make the two principal stresses of two-way wrinkled elements zero as well.It could attain the correct load responses right after the appearance of wrinkled elements,which indicates that Method II can be applied to wrinkling analysis of large-scale membrane structures.
基金Projects(51074014,51174014) supported by the National Natural Science Foundation of China
文摘Methods that can efficiently model the effects of rock joints on rock mass behavior can be beneficial in rock engineering. The suitability of equivalent rock mass(ERM) technique based upon particle methods is investigated. The ERM methodology is first validated by comparing calculated and experimental data of lab triaxial compression test on a set of cylindrical rock mass samples, each containing a single joint oriented in various dip angles. The simulated results are then used to study the stress-strain nonlinearity and failure mechanism as a function of the joint dip angle and confining stress. The anisotropy and size effects are also investigated by using multi-scale cubic ERM models subjected to triaxial compression test. The deformation and failure behavior are found to be influenced by joint degradation, the micro-crack formation in the intact rock, the interaction between two joints, and the interactions of micro-cracks and joints.
基金the National Natural Science Foundation of China(Grant Nos.52130808 and 51878566)National Key R&D Program of China(Key Projects for International Science and Technology Innovation Cooperation between Governments,Grant No.2022YFE0104300).
文摘In the longitudinal seismic deformation method for shield tunnels,one of the most commonly used is the longitudinal equivalent stiffness beam model(LES)for simulating the mechanical behavior of the lining.In this model,axial deformation and bending deformation are independent,so the equivalent stiffness is a constant value.However,the actual situation is that axial deformation and bending deformation occur simultaneously,which is not considered in LES.At present,we are not clear about the effect on the calculation results when axial deformation and bending deformation occur simultaneously.Therefore,in this paper,we improve the traditional LES by taking the relative deformation as a load and considering the coordinated deformation of axial and bending degrees of freedom.This improved model is called DNLES,and its neutral axis equations are an explicit expression.Then,we propose an iterative algorithm to solve the calculation model of the DNLES-based longitudinal seismic deformation method.Through a calculation example,we find that the internal forces based on LES are notably underestimated than those of DNLES in the compression bending zone,while are overestimated in the tension bending zone.When considering the combined effect,the maximum bending moment reached 13.7 times that of the LES model,and the axial pressure and tension were about 1.14 and 0.96 times,respectively.Further analysis reveals the coordinated deformation process in the axial and bending directions of the shield tunnel,which leads to a consequent change in equivalent stiffness.This explains why,in the longitudinal seismic deformation method,the traditional LES may result in unreasonable calculation results.
基金Project supported by the Key Project of Chinese Ministry of Education (No.106015).
文摘An equivalent continuum method is developed to analyze the effective stiffness of three-dimensional stretching dominated lattice materials. The strength and three-dimensional plastic yield surfaces are calculated for the equivalent continuum. A yielding model is formulated and compared with the results of other models. The bedding-in effect is considered to include the compliance of the lattice joints. The predicted stiffness and strength are in good agreement with the experimental data, validating the present model in the prediction of the mechanical properties of stretching dominated lattice structures.
基金supported by the National Key Research and Development Program of China(Grant No.2023YFC2907600)the National Natural Science Foundation of China(Grant Nos.42277174 and 52204260).
文摘The rock mass consists of rock blocks and structural planes,which can reduce its integrity and strength.Therefore,accurately obtaining the characteristics of the rock mass structural plane is a prerequisite for evaluating stability and designing supports in underground engineering.Currently,there are no effective testing methods for the characteristic parameters of the rock mass structural plane in underground engineering.The paper presents the digital drilling technology as a new testing method of rock mass structural planes.Flawed rock specimens with cracks of varying widths and angles were used to simulate the rock mass structural planes,and the multifunctional rock mass digital drilling test system was employed to carry out the digital drilling tests.The analysis focuses on the variation laws of drilling parameters,such as drilling pressure and drilling torque,affected by the characteristics of prefabricated cracks,and clarifies the degradation mechanism of rock equivalent compressive strength.Additionally,an identification model for the characteristic parameters of rock mass structural planes during drilling is established.The test results indicate that the average difference of the characteristics of prefabricated cracks identified by the equivalent compressive strength is 2.45°and 0.82 mm,respectively.The identification model while drilling is verified to be correct due to the high identification accuracy.Based on this,a method for testing the characteristic parameters of the surrounding rock structural plane while drilling is proposed.The research offers a theoretical and methodological foundation for precise in situ identification of structural planes of the surrounding rock in underground engineering.
基金supported by the National Natural Science Youth Foundation of China(Grant No.51309101)the Henan Province Major Scientific and Technological Projects(Grant No.172102210372)the Cooperative Project of Production,Teaching and Research in Henan Province(Grant No.18210700031)
文摘This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.
基金The study is financially supported by the National Major Research Instrument Development Project of the National Natural Science Foundation of China(Grant No.51627812)the National Natural Science Foundation of China(Grant No.52078181)the Natural Science Foundation of Hebei Province,China(Grant No.E2019202484)。
文摘Quantitative damage identification of surrounding rock is important to assess the current condition and residual strength of underground tunnels.In this work,an underground tunnel model with marble-like cementitious materials was first fabricated using the three-dimensional(3D)printing technique and then loaded to simulate its failure mode in the laboratory.Lead zirconate titanate piezoelectric(PZT)transducers were embedded in the surrounding rock around the tunnel in the process of 3D printing.A 3D monitoring network was formed to locate damage areas and evaluate damage extent during loading.Results show that as the load increased,main cracks firstly appeared above the tunnel roof and below the floor,and then they coalesced into the tunnel boundary.Finally,the tunnel model was broken into several parts.The resonant frequency and the peak of the conductance signature firstly shifted rightwards with loading due to the sealing of microcracks,and then shifted backwards after new cracks appeared.An overall increase in the root-mean-square deviation(RMSD)calculated from conductance signatures of all the PZT transducers was observed as the load(damage)increased.Damage-dependent equivalent stiffness parameters(ESPs)were calculated from the real and imaginary signatures of each PZT at different damage states.Satisfactory agreement between equivalent and experimental ESP values was achieved.Also,the relationship between the change of the ESP and the residual strength was obtained.The method paves the way for damage identification and residual strength estimation of other 3D printed structures in civil engineering.
基金Sponsored by the National Natural Science Foundation of China(Grant No.50508012)
文摘This paper proposes a new technique which introduces the high-order single-step-β method(HSM)into the experimental study on the substructure pseudo-dynamic testing(SPDT).The technique is based on the proposed concept of equivalent shear stiffness which can meet the requirement of the HSM algorithm.A study is done to theoretically validate the technique by the numerical analysis of two-storey shear building structure,in comparison of the proposed substructure pseudo-dynamic testing algorithm with the central difference method(CDM).Then,a full-scale SPDT model,the three-storey frame-supported reinforced concrete short-limb masonry shear wall structure,is designed and tested to simulate the seismic response of the corresponding six-storey structure and verify the proposed force control HSM technique.Meanwhile,the techniques of both stiffness correction and force control are suggested to control algorithmic error,control error and measurement error.The results indicate that the force control HSM can be used in the full-scale multi-degree-of-freedom(MDOF)substructure pseudo-dynamic testing before descent segment of structure restoring force properties.