We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensiti...We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.展开更多
The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The bas...The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The basic parameters of sandstone are determined through a series of static and dynamic tests,including uniaxial compression,Brazilian disc,triaxial compression under varying confining pressures,hydrostatic compression,and dynamic compression and tensile tests with a split Hopkinson pressure bar.Based on the sandstone test results from this study and previous research,a step-by-step procedure for parameter calibration is outlined,which accounts for the categories of the strength surface,equation of state(EOS),strain rate effect,and damage.The calibrated parameters are verified through numerical tests that correspond to the experimental loading conditions.Consistency between numerical results and experimental data indicates the precision and reliability of the calibrated parameters.The methodology presented in this study is scientifically sound,straightforward,and essential for improving the TSID model.Furthermore,it has the potential to contribute to other rock constitutive models,particularly new user-defined models.展开更多
Water effects on the mechanical properties of rocks have been extensively investigated through experiments and numerical models.However,few studies have established a comprehensive link between the microscopic mechani...Water effects on the mechanical properties of rocks have been extensively investigated through experiments and numerical models.However,few studies have established a comprehensive link between the microscopic mechanisms of water-related micro-crack and the constitutive behaviors of rocks.In this work,we shall propose an extended micromechanical-based plastic damage model for understanding weakening effect induced by the presence of water between micro-crack’s surfaces on quasi-brittle rocks,based on the Mori-Tanaka homogenization and irreversible thermodynamics framework.Regarding the physical mechanism,water strengthens micro-crack propagation,which induces damage evolution during the pre-and post-stage,and weakens the elastic effective properties of rock matrix.After proposing a special calibration procedure for the determination of model parameters based on the laboratory compression tests,the proposed micromechanical-based model is verified by comparing the model predictions to the experimental results.The model effectively captures the mechanical behaviors of quasibrittle rocks subjected to the weakening effects of water.展开更多
This paper presents a new finite element model updating method for estimating structural parameters and detecting structural damage location and severity based on the structural responses(output-only data).The method ...This paper presents a new finite element model updating method for estimating structural parameters and detecting structural damage location and severity based on the structural responses(output-only data).The method uses the sensitivity relation of transmissibility data through a least-squares algorithm and appropriate normalization of the extracted equations.The proposed transmissibility-based sensitivity equation produces a more significant number of equations than the sensitivity equations based on the frequency response function(FRF),which can estimate the structural parameters with higher accuracy.The abilities of the proposed method are assessed by using numerical data of a two-story two-bay frame model and a plate structure model.In evaluating different damage cases,the number,location,and stiffness reduction of the damaged elements and the severity of the simulated damage have been accurately identified.The reliability and stability of the presented method against measurement and modeling errors are examined using error-contaminated data.The parameter estimation results prove the method’s capabilities as an accurate model updating algorithm.展开更多
To investigate the specific creep behavior of ultra-deep buried salt during oil and gas exploitation,a set of triaxial creep experiments was conducted at elevated temperatures with constant axial pressure and unloadin...To investigate the specific creep behavior of ultra-deep buried salt during oil and gas exploitation,a set of triaxial creep experiments was conducted at elevated temperatures with constant axial pressure and unloading confining pressure conditions.Experimental results show that the salt sample deforms more significantly with the increase of applied temperature and deviatoric loading.The accelerated creep phase is not occurring until the applied temperature reaches 130℃,and higher temperature is beneficial to the occurrence of accelerated creep.To describe the specific creep behavior,a novel three-dimensional(3D)creep constitutive model is developed that incorporates the thermal and mechanical variables into mechanical elements.Subsequently,the standard particle swarm optimization(SPSO)method is adopted to fit the experimental data,and the sensibility of key model parameters is analyzed to further illustrate the model function.As a result,the model can accurately predict the creep behavior of salt under the coupled thermo-mechanical effect in deep-buried condition.Based on the research results,the creep mechanical behavior of wellbore shrinkage is predicted in deep drilling projects crossing salt layer,which has practical implications for deep rock mechanics problems.展开更多
In order to evaluate the performance of deep geological disposal of radioactive waste,an underground research laboratory(URL)was constructed by Andra in the Callovo-Oxfordian(COx)claystone formation at the Meuse/Haute...In order to evaluate the performance of deep geological disposal of radioactive waste,an underground research laboratory(URL)was constructed by Andra in the Callovo-Oxfordian(COx)claystone formation at the Meuse/Haute-Marne(MHM).The construction of URL induced the excavation damage of host formations,and the ventilation in the galleries desaturated the host formation close to the gallery wall.Moreover,it is expected that the mechanical behaviour of COx claystone is time-dependent.This study presents a constitutive model developed to describe the viscoplastic behaviour of unsaturated and damaged COx claystone.In this model,the unsaturation effect is considered by adopting the Bishop effective stress and the van Genuchten(VG)water retention model.In terms of the viscoplastic behaviour,the nonstationary flow surface(NSFS)theory for unsaturated soils is used with consideration of the coupled effects of strain rate and suction on the yield stress.A progressive hardening law is adopted.Meanwhile,a non-associated flow rule is used,which is similar to that in Barcelona basic model(BBM).In addition,to describe the damage effect induced by suction change and viscoplastic loading,a damage function is defined based on the crack volume proportion.This damage function contains two variables:unsaturated effective stress and viscoplastic volumetric strain,with the related parameters determined based on the mercury intrusion porosimetry(MIP)tests.For the model validation,different tests on COx claystone under different loading paths are simulated.Comparisons between experimental and simulated results indicated that the present model is able to well describe the viscoplastic behaviour of damaged COx claystone,including swelling/shrinkage,triaxial extension and compression,and triaxial creep.展开更多
The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element metho...The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.展开更多
The damage evolution process of non-penetrating cracks often causes some unexpected engineering disasters.Gypsum specimens containing non-penetrating crack(s)are used to study the damage evolution and characteristics ...The damage evolution process of non-penetrating cracks often causes some unexpected engineering disasters.Gypsum specimens containing non-penetrating crack(s)are used to study the damage evolution and characteristics under cyclic loading.The results show that under cyclic loading,the relationship between the number of non-penetrating crack(s)and the characteristic parameters(cyclic number,peak stress,peak strain,failure stress,and failure strain)of the pre-cracked specimens can be represented by a decreasing linear function.The damage evolution equation is fitted by calibrating the accumulative plastic strain for each cycle,and the damage constitutive equation is proposed by the concept of effective stress.Additionally,non-penetrating cracks are more likely to cause uneven stress distribution,damage accumulation,and local failure of specimen.The local failure can change the stress distribution and relieve the inhibition of non-penetrating crack extension and eventually cause a dramatic destruction of the specimen.Therefore,the evolution process caused by non-penetrating cracks can be regarded as one of the important reasons for inducing rockburst.These results are expected to improve the understanding of the process of spalling formation and rockburst and can be used to analyze the stability of rocks or rock structures.展开更多
Introducing Neutral Polymeric bonding agents(NPBA) into the Nitrate Ester Plasticized Polyether(NEPE)propellant could improve the adhesion between filler/matrix interface, thereby contributing to the development of ne...Introducing Neutral Polymeric bonding agents(NPBA) into the Nitrate Ester Plasticized Polyether(NEPE)propellant could improve the adhesion between filler/matrix interface, thereby contributing to the development of new generations of the NEPE propellant with better mechanical properties. Therefore,understanding the effects of NPBA on the deformation and damage evolution of the NEPE propellant is fundamental to material design and applications. This paper studies the uniaxial tensile and stress relaxation responses of the NEPE propellant with different amounts of NPBA. The damage evolution in terms of interface debonding is further investigated using a cohesive-zone model(CZM). Experimental results show that the initial modulus and strength of the NEPE propellant increase with the increasing amount of NPBA while the elongation decreases. Meanwhile, the relaxation rate slows down and a higher long-term equilibrium modulus is reached. Experimental and numerical analyses indicate that interface debonding and crack propagation along filler-matrix interface are the dominant damage mechanism for the samples with a low amount of NPBA, while damage localization and crack advancement through the matrix are predominant for the ones with a high amount of NPBA. Finally, crosslinking density tests and simulation results also show that the effect of the bonding agent is interfacial rather than due to the overall crosslinking density change of the binder.展开更多
After the excavation of the roadway,the original stress balance is destroyed,resulting in the redistribution of stress and the formation of an excavation damaged zone(EDZ)around the roadway.The thickness of EDZ is the...After the excavation of the roadway,the original stress balance is destroyed,resulting in the redistribution of stress and the formation of an excavation damaged zone(EDZ)around the roadway.The thickness of EDZ is the key basis for roadway stability discrimination and support structure design,and it is of great engineering significance to accurately predict the thickness of EDZ.Considering the advantages of machine learning(ML)in dealing with high-dimensional,nonlinear problems,a hybrid prediction model based on the random forest(RF)algorithm is developed in this paper.The model used the dragonfly algorithm(DA)to optimize two hyperparameters in RF,namely mtry and ntree,and used mean absolute error(MAE),rootmean square error(RMSE),determination coefficient(R^(2)),and variance accounted for(VAF)to evaluatemodel prediction performance.A database containing 217 sets of data was collected,with embedding depth(ED),drift span(DS),surrounding rock mass strength(RMS),joint index(JI)as input variables,and the excavation damaged zone thickness(EDZT)as output variable.In addition,four classic models,back propagation neural network(BPNN),extreme learning machine(ELM),radial basis function network(RBF),and RF were compared with the DA-RF model.The results showed that the DARF mold had the best prediction performance(training set:MAE=0.1036,RMSE=0.1514,R^(2)=0.9577,VAF=94.2645;test set:MAE=0.1115,RMSE=0.1417,R^(2)=0.9423,VAF=94.0836).The results of the sensitivity analysis showed that the relative importance of each input variable was DS,ED,RMS,and JI from low to high.展开更多
As one of the most common occurring geological landforms in deep rock formations, the dynamic mechanical properties of layered composite rock bodies under impact loading have been widely studied by scholars. To study ...As one of the most common occurring geological landforms in deep rock formations, the dynamic mechanical properties of layered composite rock bodies under impact loading have been widely studied by scholars. To study the dynamic properties of soft and hard composite rocks with different thickness ratios, this paper utilizes cement, quartz sand and gypsum powder to construct soft and hard composite rock specimens and utilizes a combination of indoor tests, numerical calculations, and theoretical analyses to investigate the mechanical properties of soft and hard composite rock bodies. The test results reveal that:(1) When the proportion of hard rock increases from 20% to 50%, the strength of the combined rock body increases by 69.14 MPa and 87 MPa when the hard rock face and soft rock face are loaded, respectively;however, when the proportion of hard rock is the same, the compressive strength of the hard rock face impact is 9%-17% greater than that of the soft rock face impact;(2) When a specimen of soft and hard combined rock body is subjected to impact loading, the damage mode involves mixed tension and shear damage, and the cracks generally first appear at the ends of the specimen, then develop on the laminar surface from the impact surface, and finally end in the overall damage of the soft rock part. The development rate and the total number of cracks in the same specimen when the hard rock face is impacted are significantly greater than those when the soft rock face is impacted;(3) By introducing Weibull’s statistical strength theory to establish the damage variables of soft-hard combined rock bodies, combined with the DP strength criterion, the damage model and the Kelvin body are concatenated to obtain a statistical damage constitutive model, which can better fit the full stress-strain curve of soft-hard combined rock body specimens under a single impact load.展开更多
Due to differences in the properties of composition materials and construction techniques,unreinforced masonry is characterized by low strength,anisotropy,nonuniformity,and low ductility.In order to accurately simulat...Due to differences in the properties of composition materials and construction techniques,unreinforced masonry is characterized by low strength,anisotropy,nonuniformity,and low ductility.In order to accurately simulate the mechanical behavior of unreinforced brick masonry walls under static and dynamic loads,a new elastoplastic damage constitutive model was proposed and the corresponding subroutine was developed based on the concrete material constitutive model.In the proposed constitutive model,the Rankine strength theory and the Drucker-Prager strength theory were used to define the tensile and compressive yield surface function of materials,respectively.Moreover,the stress updating algorithm was modified to consider the tensile plastic permanent deformation of masonry materials.To verify the accuracy of the proposed constitutive model,numerical simulations of the brick masonry under monotonic and cyclic uniaxial tension and compression loads were carried out.Comparisons among the numerical and theoretical and experimental results show that the proposed model can properly reflect the masonry material mechanical properties.Furthermore,the numerical models of four pieces of masonry walls with different mortar strengths were established.Low cyclic loadings were applied and the results show that the proposed constitutive model can properly simulate the wall shear failure characteristics,and the force-displacement hysteretic curves obtained by numerical simulation are in good agreement with the tests.Overall,the proposed elastic-plastic damage constitutive model can simulate the nonlinear behavior of unreinforced brick masonry walls very well,and can be used to predict the structural response of masonry walls.展开更多
Given the difficulty in accurately evaluating the fatigue performance of large composite wind turbine blades(referred to as blades),this paper takes the main beam structure of the blade with a rectangular cross-sectio...Given the difficulty in accurately evaluating the fatigue performance of large composite wind turbine blades(referred to as blades),this paper takes the main beam structure of the blade with a rectangular cross-sectionas the simulation object and establishes a composite laminate rectangular beam structure that simultaneouslyincludes the flange,web,and adhesive layer,referred to as the blade main beam sub-structure specimen,throughthe definition of blade sub-structures.This paper examines the progressive damage evolution law of the compositelaminate rectangular beam utilizing an improved 3D Hashin failure criterion,cohesive zone model,B-K failurecriterion,and computer simulation technology.Under static loading,the layup angle of the anti-shear web hasa close relationship with the static load-carrying capacity of the composite laminate rectangular beam;under fatigueloading,the fatigue damage will first occur in the lower flange adhesive area of the whole composite laminaterectangular beam and ultimately result in the fracture failure of the entire structure.These results provide a theoreticalreference and foundation for evaluating and predicting the fatigue performance of the blade main beamstructure and even the full-size blade.展开更多
Block piezoelectric ceramics(PZTs)are often used in impact igniters to provide activation energy for electric initiators.Under the action of strong impact stress,PZTs release electric energy accompanied by crack initi...Block piezoelectric ceramics(PZTs)are often used in impact igniters to provide activation energy for electric initiators.Under the action of strong impact stress,PZTs release electric energy accompanied by crack initiation,propagation and crushing.At present,the electrical output performance of PZTs in projectile is usually calculated by quasi-static piezoelectric equation without considering the dynamic effect caused by strong impact and the influence of crack propagation on material properties.So the ignition parameters are always not accurately predicted.To tackle this,a PZT dynamic damage constitutive model considering crack propagation is established based on the dynamic impact test and the crack propagation theory of brittle materials.The model is then embedded into the ABAQUS subroutine and used to simulate the electromechanical response of the impact igniter during the impact of a small caliber projectile on the target.Meanwhile,the experiments of projectile with impact igniter impact on the target are carried out.The comparison between experimental and numerical simulation results show that the established dynamic damage model can effectively predict the dynamic electromechanical response of PZTs in the missile service environment.展开更多
Aiming at the problem of temperature-mechanics-chemical(T-M-C)action encountered by rocks in underground engineering,sandstone was selected as the sample for acid corrosion treatment at pH 1,3,5 and 7,the acid corrosi...Aiming at the problem of temperature-mechanics-chemical(T-M-C)action encountered by rocks in underground engineering,sandstone was selected as the sample for acid corrosion treatment at pH 1,3,5 and 7,the acid corrosion treated samples were then subjected to high-temperature experiments at 25,300,600,and 900℃,and triaxial compression experiments were conducted in the laboratory.The experimental results show that the superposition of chemical damage and thermal damage has a significant impact on the quality,wave velocity,porosity and compression failure characteristics of the rock.Based on the Lemaitre strain equivalent hypothesis theory,the damage degree of rock material was described by introducing damage variables,and the spatial mobilized plane(SMP)criterion was adopted.The damage constitutive model can well reflect the stress-strain characteristics of the rock triaxial compression process,which verified the rationality and reliability of the model parameters.The experiment and constitutive model analyzed the change law of mechanical properties of rock after chemical corrosion and high temperature thermal damage,which had certain practical significance for rock engineering construction.展开更多
Non aqueous reactive polymer materials produced by the reaction of isocyanate and polyol have been widely used in infrastructure construction,which may be subjected to explosion loads during complex service conditions...Non aqueous reactive polymer materials produced by the reaction of isocyanate and polyol have been widely used in infrastructure construction,which may be subjected to explosion loads during complex service conditions.The blast response of composite materials is a crucial aspect for applications in engineering structures potentially subjected to extreme loadings.In this work,damage caused to rebar reinforced polymer slabs by surface explosive charges was studied experimentally and numerically.A total of 6 field tests were carried out to investigate the performances of the failure modes of rebar reinforced polymer slabs under contact and near-field explosions.The influence of explosive quantity(10-40 g)and stand-off distances(0-20 cm)at the damage modes were studied.The results show that the failure modes of rebar reinforced polymer slabs under near-field explosion mainly were bending and surface spalling,while under the impact of contact explosion,the failure modes were craters of the top surface,spalling of the bottom surface,and middle perforation.Furthermore,a detailed fully coupled model was developed and validated with the test data.The influences of explosive quantity and slab thickness on rebar reinforced polymer slabs under contact explosion were studied.Based on this,the calculation formula between breach diameter,explosive quantity,and slab thickness is fitted.展开更多
Over the past century,the safety of dams has gradually attracted attention from all parties.Research on the dynamic response and damage evolution of dams under extreme loads is the basis of dam safety issues.In recent...Over the past century,the safety of dams has gradually attracted attention from all parties.Research on the dynamic response and damage evolution of dams under extreme loads is the basis of dam safety issues.In recent decades,scholars have studied the responses of dams under earthquake loads,but there is still much room for improvement in experimental and theoretical research on small probability loads such as explosions.In this paper,a 50-m-high concrete gravity dam is used as a prototype dam,and a water explosion model test of a 2.5-m-high concrete gravity dam is designed.The water pressure and the acceleration response of the dam body in the test are analysed.The pressure characteristics and dynamic response of the dam body are assessed.Taking the dam damage test as an example,a numerical model of concrete gravity dam damage is established,and the damage evolution of the dam body is analysed.By combining experiments and numerical simulations,the damage characteristics of the dam body under the action of different charge water explosions are clarified.The integrity of the dam body is well maintained under the action of a small-quantity water explosion,and the dynamic response of the dam body is mainly caused by the shock wave.Both the shock wave and the bubble pulsation cause the dam body to accelerate,and the peak acceleration of the dam body under the action of the bubble pulsation is only one percent of the peak acceleration of the dam body under the action of the shock wave.When subjected to explosions in large quantities of water,the dam body is seriously damaged.Under the action of a shock wave,the dam body produces a secondary acceleration response,which is generated by an internal interaction after the dam body is damaged.The damage evolution process of the dam body under the action of a large-scale water explosion is analysed,and it is found that the shock wave pressure of the water explosion causes local damage to the dam body facing the explosion.After the peak value of the shock wave,the impulse continues to act on the dam body,causing cumulative damage and damage inside the dam body.展开更多
The interaction between soil and marine structures like submarine pipeline/pipe pile/suction caisson is a complicated geotechnical mechanism process.In this study,the interface is discretized into multiple mesoscopic ...The interaction between soil and marine structures like submarine pipeline/pipe pile/suction caisson is a complicated geotechnical mechanism process.In this study,the interface is discretized into multiple mesoscopic contact elements that are damaged randomly throughout the shearing process due to the natural heterogeneity.The evolution equation of damage variable is developed based on the Weibull function,which is able to cover a rather wide range of distribution shapes by only two parameters,making it applicable for varying scenarios.Accordingly,a statistical damage model is established by incorporating Mohr–Coulomb strength criterion,in which the interfacial residual strength is considered whereby the strain softening behavior can be described.A concept of“semi-softening”characteristic point on shear stress–displacement curve is proposed for effectively modeling the evolution of strain softening.Finally,a series of ring shear tests of the interfaces between fine sea sand and smooth/rough steel surfaces are conducted.The predicted results using the proposed model are compared with experimental data of this study as well as some results from existing literature,indicating that the model has a good performance in modeling the progressive failure and strain softening behavior for various types of soil–structure interfaces.展开更多
The creep characteristics of rock under high crustal stress are of important influence on the long‐term stability of deep rock engineering.To study the creep characteristics and engineering application of sandstone u...The creep characteristics of rock under high crustal stress are of important influence on the long‐term stability of deep rock engineering.To study the creep characteristics and engineering application of sandstone under high crustal stress,this study constructed nonlinear creep damage(NCD)constitutive mode based on the triaxial graded loading‒unloading creep test of sandstone in the Yuezhishan Tunnel.A numerical NCD constitutive model and a breakable lining(BL)model were developed based on FLAC3D and then applied to the stability analysis of the Yuezhishan Tunnel.Based on the creep test results of sandstone,a power function of creep rate and stress level was constructed,by which the long‐term strength was solved.The results show that the long‐term strength of the red sandstone based on the related function of the steady‐state creep rate and stress level is close to the measured stress value in engineering.The NCD model considering damage factors reflects the instantaneous and viscoelastic plasticity deformation characteristics of the red sandstone.The numerical NCD constitutive model and the BL model can reflect surrounding rock deformation characteristics and lining failure characteristics in practical engineering.The research results provide theoretical references for long‐term stability analysis of rock engineering and the deformation control of surrounding rock under high crustal stress.展开更多
Damage statistical mechanics model of horizontal section height in the top caving was constructed in the paper. The influence factors including supporting pressure, dip angle and characteristic of coal on horizontal s...Damage statistical mechanics model of horizontal section height in the top caving was constructed in the paper. The influence factors including supporting pressure, dip angle and characteristic of coal on horizontal section height were analyzed as well. By terms of the practice project analysis, the horizontal section height increases with the increase of dip angle β and thickness of coal seam M. Dip angle of coal seam β has tremendous impact on horizontal section height, while thickness of coal seam M has slight impact. When thickness of coal seam is below 10m, horizontal section height increases sharply. While thickness exceeds 15m, it is not major factor influencing on horizontal section height any long.展开更多
基金funding received by a grant from the Natural Sciences and Engineering Research Council of Canada(NSERC)(Grant No.CRDPJ 469057e14).
文摘We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.
基金funded by the National Natural Science Foundation of China(Grant No.12272247)National Key Project(Grant No.GJXM92579)Major Research and Development Project of Metallurgical Corporation of China Ltd.in the Non-Steel Field(Grant No.2021-5).
文摘The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The basic parameters of sandstone are determined through a series of static and dynamic tests,including uniaxial compression,Brazilian disc,triaxial compression under varying confining pressures,hydrostatic compression,and dynamic compression and tensile tests with a split Hopkinson pressure bar.Based on the sandstone test results from this study and previous research,a step-by-step procedure for parameter calibration is outlined,which accounts for the categories of the strength surface,equation of state(EOS),strain rate effect,and damage.The calibrated parameters are verified through numerical tests that correspond to the experimental loading conditions.Consistency between numerical results and experimental data indicates the precision and reliability of the calibrated parameters.The methodology presented in this study is scientifically sound,straightforward,and essential for improving the TSID model.Furthermore,it has the potential to contribute to other rock constitutive models,particularly new user-defined models.
基金financially supported by the National Natural Science Foundation of China(Nos.42001053 and 42277147)the General Scientific Research Fund of Zhejiang Provincial Education Department(No.Y202352363)the University Natural Science Foundation of Jiangsu Province(No.23KJD130001)。
文摘Water effects on the mechanical properties of rocks have been extensively investigated through experiments and numerical models.However,few studies have established a comprehensive link between the microscopic mechanisms of water-related micro-crack and the constitutive behaviors of rocks.In this work,we shall propose an extended micromechanical-based plastic damage model for understanding weakening effect induced by the presence of water between micro-crack’s surfaces on quasi-brittle rocks,based on the Mori-Tanaka homogenization and irreversible thermodynamics framework.Regarding the physical mechanism,water strengthens micro-crack propagation,which induces damage evolution during the pre-and post-stage,and weakens the elastic effective properties of rock matrix.After proposing a special calibration procedure for the determination of model parameters based on the laboratory compression tests,the proposed micromechanical-based model is verified by comparing the model predictions to the experimental results.The model effectively captures the mechanical behaviors of quasibrittle rocks subjected to the weakening effects of water.
文摘This paper presents a new finite element model updating method for estimating structural parameters and detecting structural damage location and severity based on the structural responses(output-only data).The method uses the sensitivity relation of transmissibility data through a least-squares algorithm and appropriate normalization of the extracted equations.The proposed transmissibility-based sensitivity equation produces a more significant number of equations than the sensitivity equations based on the frequency response function(FRF),which can estimate the structural parameters with higher accuracy.The abilities of the proposed method are assessed by using numerical data of a two-story two-bay frame model and a plate structure model.In evaluating different damage cases,the number,location,and stiffness reduction of the damaged elements and the severity of the simulated damage have been accurately identified.The reliability and stability of the presented method against measurement and modeling errors are examined using error-contaminated data.The parameter estimation results prove the method’s capabilities as an accurate model updating algorithm.
基金This research was financially supported by the Scientific and technological research projects in Sichuan province(Grant Nos.2022YFSY0007 and 2021YFH0010)the National Scientific Science Foundation of China(Grant No.U20A20266).
文摘To investigate the specific creep behavior of ultra-deep buried salt during oil and gas exploitation,a set of triaxial creep experiments was conducted at elevated temperatures with constant axial pressure and unloading confining pressure conditions.Experimental results show that the salt sample deforms more significantly with the increase of applied temperature and deviatoric loading.The accelerated creep phase is not occurring until the applied temperature reaches 130℃,and higher temperature is beneficial to the occurrence of accelerated creep.To describe the specific creep behavior,a novel three-dimensional(3D)creep constitutive model is developed that incorporates the thermal and mechanical variables into mechanical elements.Subsequently,the standard particle swarm optimization(SPSO)method is adopted to fit the experimental data,and the sensibility of key model parameters is analyzed to further illustrate the model function.As a result,the model can accurately predict the creep behavior of salt under the coupled thermo-mechanical effect in deep-buried condition.Based on the research results,the creep mechanical behavior of wellbore shrinkage is predicted in deep drilling projects crossing salt layer,which has practical implications for deep rock mechanics problems.
文摘In order to evaluate the performance of deep geological disposal of radioactive waste,an underground research laboratory(URL)was constructed by Andra in the Callovo-Oxfordian(COx)claystone formation at the Meuse/Haute-Marne(MHM).The construction of URL induced the excavation damage of host formations,and the ventilation in the galleries desaturated the host formation close to the gallery wall.Moreover,it is expected that the mechanical behaviour of COx claystone is time-dependent.This study presents a constitutive model developed to describe the viscoplastic behaviour of unsaturated and damaged COx claystone.In this model,the unsaturation effect is considered by adopting the Bishop effective stress and the van Genuchten(VG)water retention model.In terms of the viscoplastic behaviour,the nonstationary flow surface(NSFS)theory for unsaturated soils is used with consideration of the coupled effects of strain rate and suction on the yield stress.A progressive hardening law is adopted.Meanwhile,a non-associated flow rule is used,which is similar to that in Barcelona basic model(BBM).In addition,to describe the damage effect induced by suction change and viscoplastic loading,a damage function is defined based on the crack volume proportion.This damage function contains two variables:unsaturated effective stress and viscoplastic volumetric strain,with the related parameters determined based on the mercury intrusion porosimetry(MIP)tests.For the model validation,different tests on COx claystone under different loading paths are simulated.Comparisons between experimental and simulated results indicated that the present model is able to well describe the viscoplastic behaviour of damaged COx claystone,including swelling/shrinkage,triaxial extension and compression,and triaxial creep.
基金funded by the Research Project Supported by Shanxi Scholarship Council of China(2022-067)the Opening Project of State Key Laboratory of Explosion Science and Technology(Beijing Institute of Technology)(KFJJ22-14M).
文摘The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.
基金supported by the National Natural Science Foundation of China(Nos.52204092 and 52274203).
文摘The damage evolution process of non-penetrating cracks often causes some unexpected engineering disasters.Gypsum specimens containing non-penetrating crack(s)are used to study the damage evolution and characteristics under cyclic loading.The results show that under cyclic loading,the relationship between the number of non-penetrating crack(s)and the characteristic parameters(cyclic number,peak stress,peak strain,failure stress,and failure strain)of the pre-cracked specimens can be represented by a decreasing linear function.The damage evolution equation is fitted by calibrating the accumulative plastic strain for each cycle,and the damage constitutive equation is proposed by the concept of effective stress.Additionally,non-penetrating cracks are more likely to cause uneven stress distribution,damage accumulation,and local failure of specimen.The local failure can change the stress distribution and relieve the inhibition of non-penetrating crack extension and eventually cause a dramatic destruction of the specimen.Therefore,the evolution process caused by non-penetrating cracks can be regarded as one of the important reasons for inducing rockburst.These results are expected to improve the understanding of the process of spalling formation and rockburst and can be used to analyze the stability of rocks or rock structures.
基金National Natural Science Foundation of China(U22B20131)for supporting this project.
文摘Introducing Neutral Polymeric bonding agents(NPBA) into the Nitrate Ester Plasticized Polyether(NEPE)propellant could improve the adhesion between filler/matrix interface, thereby contributing to the development of new generations of the NEPE propellant with better mechanical properties. Therefore,understanding the effects of NPBA on the deformation and damage evolution of the NEPE propellant is fundamental to material design and applications. This paper studies the uniaxial tensile and stress relaxation responses of the NEPE propellant with different amounts of NPBA. The damage evolution in terms of interface debonding is further investigated using a cohesive-zone model(CZM). Experimental results show that the initial modulus and strength of the NEPE propellant increase with the increasing amount of NPBA while the elongation decreases. Meanwhile, the relaxation rate slows down and a higher long-term equilibrium modulus is reached. Experimental and numerical analyses indicate that interface debonding and crack propagation along filler-matrix interface are the dominant damage mechanism for the samples with a low amount of NPBA, while damage localization and crack advancement through the matrix are predominant for the ones with a high amount of NPBA. Finally, crosslinking density tests and simulation results also show that the effect of the bonding agent is interfacial rather than due to the overall crosslinking density change of the binder.
基金funded by the National Science Foundation of China(42177164)the Distinguished Youth Science Foundation of Hunan Province of China(2022JJ10073)the Innovation-Driven Project of Central South University(2020CX040).
文摘After the excavation of the roadway,the original stress balance is destroyed,resulting in the redistribution of stress and the formation of an excavation damaged zone(EDZ)around the roadway.The thickness of EDZ is the key basis for roadway stability discrimination and support structure design,and it is of great engineering significance to accurately predict the thickness of EDZ.Considering the advantages of machine learning(ML)in dealing with high-dimensional,nonlinear problems,a hybrid prediction model based on the random forest(RF)algorithm is developed in this paper.The model used the dragonfly algorithm(DA)to optimize two hyperparameters in RF,namely mtry and ntree,and used mean absolute error(MAE),rootmean square error(RMSE),determination coefficient(R^(2)),and variance accounted for(VAF)to evaluatemodel prediction performance.A database containing 217 sets of data was collected,with embedding depth(ED),drift span(DS),surrounding rock mass strength(RMS),joint index(JI)as input variables,and the excavation damaged zone thickness(EDZT)as output variable.In addition,four classic models,back propagation neural network(BPNN),extreme learning machine(ELM),radial basis function network(RBF),and RF were compared with the DA-RF model.The results showed that the DARF mold had the best prediction performance(training set:MAE=0.1036,RMSE=0.1514,R^(2)=0.9577,VAF=94.2645;test set:MAE=0.1115,RMSE=0.1417,R^(2)=0.9423,VAF=94.0836).The results of the sensitivity analysis showed that the relative importance of each input variable was DS,ED,RMS,and JI from low to high.
基金supported by the Xi’an Key Laboratory of Geotechnical and Underground Engineering Open Fund Project (XKLGUEKF20-03)the Natural Science Basic Research Program of Shaanxi Province General Project-Youth Project(2024JC-YBQN-0258)。
文摘As one of the most common occurring geological landforms in deep rock formations, the dynamic mechanical properties of layered composite rock bodies under impact loading have been widely studied by scholars. To study the dynamic properties of soft and hard composite rocks with different thickness ratios, this paper utilizes cement, quartz sand and gypsum powder to construct soft and hard composite rock specimens and utilizes a combination of indoor tests, numerical calculations, and theoretical analyses to investigate the mechanical properties of soft and hard composite rock bodies. The test results reveal that:(1) When the proportion of hard rock increases from 20% to 50%, the strength of the combined rock body increases by 69.14 MPa and 87 MPa when the hard rock face and soft rock face are loaded, respectively;however, when the proportion of hard rock is the same, the compressive strength of the hard rock face impact is 9%-17% greater than that of the soft rock face impact;(2) When a specimen of soft and hard combined rock body is subjected to impact loading, the damage mode involves mixed tension and shear damage, and the cracks generally first appear at the ends of the specimen, then develop on the laminar surface from the impact surface, and finally end in the overall damage of the soft rock part. The development rate and the total number of cracks in the same specimen when the hard rock face is impacted are significantly greater than those when the soft rock face is impacted;(3) By introducing Weibull’s statistical strength theory to establish the damage variables of soft-hard combined rock bodies, combined with the DP strength criterion, the damage model and the Kelvin body are concatenated to obtain a statistical damage constitutive model, which can better fit the full stress-strain curve of soft-hard combined rock body specimens under a single impact load.
基金National Key Research and Development Program of China under Grant Nos.2018YFC1504400 and 2019YFC1509301Natural Science Foundation of China under Grant No.52078471Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration under Grant No.2019EEEVL0402。
文摘Due to differences in the properties of composition materials and construction techniques,unreinforced masonry is characterized by low strength,anisotropy,nonuniformity,and low ductility.In order to accurately simulate the mechanical behavior of unreinforced brick masonry walls under static and dynamic loads,a new elastoplastic damage constitutive model was proposed and the corresponding subroutine was developed based on the concrete material constitutive model.In the proposed constitutive model,the Rankine strength theory and the Drucker-Prager strength theory were used to define the tensile and compressive yield surface function of materials,respectively.Moreover,the stress updating algorithm was modified to consider the tensile plastic permanent deformation of masonry materials.To verify the accuracy of the proposed constitutive model,numerical simulations of the brick masonry under monotonic and cyclic uniaxial tension and compression loads were carried out.Comparisons among the numerical and theoretical and experimental results show that the proposed model can properly reflect the masonry material mechanical properties.Furthermore,the numerical models of four pieces of masonry walls with different mortar strengths were established.Low cyclic loadings were applied and the results show that the proposed constitutive model can properly simulate the wall shear failure characteristics,and the force-displacement hysteretic curves obtained by numerical simulation are in good agreement with the tests.Overall,the proposed elastic-plastic damage constitutive model can simulate the nonlinear behavior of unreinforced brick masonry walls very well,and can be used to predict the structural response of masonry walls.
基金the Science and Technology Programs of Gansu Province(Grant Nos.21JR1RA248,23YFGA0050)the Young Scholars Science Foundation of Lanzhou Jiaotong University(Grant Nos.2020039,2020017)+2 种基金the Special Funds for Guiding Local Scientific and Technological Development by the Central Government(Grant No.22ZY1QA005)the National Natural Science Foundation of China(Grant No.72361019)the Gansu Provincial Outstanding Graduate Students Innovation Star Program(Grant No.2023CXZX-574).
文摘Given the difficulty in accurately evaluating the fatigue performance of large composite wind turbine blades(referred to as blades),this paper takes the main beam structure of the blade with a rectangular cross-sectionas the simulation object and establishes a composite laminate rectangular beam structure that simultaneouslyincludes the flange,web,and adhesive layer,referred to as the blade main beam sub-structure specimen,throughthe definition of blade sub-structures.This paper examines the progressive damage evolution law of the compositelaminate rectangular beam utilizing an improved 3D Hashin failure criterion,cohesive zone model,B-K failurecriterion,and computer simulation technology.Under static loading,the layup angle of the anti-shear web hasa close relationship with the static load-carrying capacity of the composite laminate rectangular beam;under fatigueloading,the fatigue damage will first occur in the lower flange adhesive area of the whole composite laminaterectangular beam and ultimately result in the fracture failure of the entire structure.These results provide a theoreticalreference and foundation for evaluating and predicting the fatigue performance of the blade main beamstructure and even the full-size blade.
基金supported by the National Natural Science Foundation of China(Grant No.12172232)the project of Key Laboratory of Impact and Safety Engineering(Ningbo University,China)+1 种基金Ministry of Education(CJ202206)supported by the scientific research support plan of introducing high-level talents from Shenyang Ligong University。
文摘Block piezoelectric ceramics(PZTs)are often used in impact igniters to provide activation energy for electric initiators.Under the action of strong impact stress,PZTs release electric energy accompanied by crack initiation,propagation and crushing.At present,the electrical output performance of PZTs in projectile is usually calculated by quasi-static piezoelectric equation without considering the dynamic effect caused by strong impact and the influence of crack propagation on material properties.So the ignition parameters are always not accurately predicted.To tackle this,a PZT dynamic damage constitutive model considering crack propagation is established based on the dynamic impact test and the crack propagation theory of brittle materials.The model is then embedded into the ABAQUS subroutine and used to simulate the electromechanical response of the impact igniter during the impact of a small caliber projectile on the target.Meanwhile,the experiments of projectile with impact igniter impact on the target are carried out.The comparison between experimental and numerical simulation results show that the established dynamic damage model can effectively predict the dynamic electromechanical response of PZTs in the missile service environment.
文摘Aiming at the problem of temperature-mechanics-chemical(T-M-C)action encountered by rocks in underground engineering,sandstone was selected as the sample for acid corrosion treatment at pH 1,3,5 and 7,the acid corrosion treated samples were then subjected to high-temperature experiments at 25,300,600,and 900℃,and triaxial compression experiments were conducted in the laboratory.The experimental results show that the superposition of chemical damage and thermal damage has a significant impact on the quality,wave velocity,porosity and compression failure characteristics of the rock.Based on the Lemaitre strain equivalent hypothesis theory,the damage degree of rock material was described by introducing damage variables,and the spatial mobilized plane(SMP)criterion was adopted.The damage constitutive model can well reflect the stress-strain characteristics of the rock triaxial compression process,which verified the rationality and reliability of the model parameters.The experiment and constitutive model analyzed the change law of mechanical properties of rock after chemical corrosion and high temperature thermal damage,which had certain practical significance for rock engineering construction.
基金supported by the National Natural Science Foundation of China(Grant Nos.52009126,51939008)Foundation of Hubei Key Laboratory of Blasting Engineering(Grant No.BL202104)First-class Project Special Funding of Yellow River Laboratory(No.YRL22IR08)。
文摘Non aqueous reactive polymer materials produced by the reaction of isocyanate and polyol have been widely used in infrastructure construction,which may be subjected to explosion loads during complex service conditions.The blast response of composite materials is a crucial aspect for applications in engineering structures potentially subjected to extreme loadings.In this work,damage caused to rebar reinforced polymer slabs by surface explosive charges was studied experimentally and numerically.A total of 6 field tests were carried out to investigate the performances of the failure modes of rebar reinforced polymer slabs under contact and near-field explosions.The influence of explosive quantity(10-40 g)and stand-off distances(0-20 cm)at the damage modes were studied.The results show that the failure modes of rebar reinforced polymer slabs under near-field explosion mainly were bending and surface spalling,while under the impact of contact explosion,the failure modes were craters of the top surface,spalling of the bottom surface,and middle perforation.Furthermore,a detailed fully coupled model was developed and validated with the test data.The influences of explosive quantity and slab thickness on rebar reinforced polymer slabs under contact explosion were studied.Based on this,the calculation formula between breach diameter,explosive quantity,and slab thickness is fitted.
文摘Over the past century,the safety of dams has gradually attracted attention from all parties.Research on the dynamic response and damage evolution of dams under extreme loads is the basis of dam safety issues.In recent decades,scholars have studied the responses of dams under earthquake loads,but there is still much room for improvement in experimental and theoretical research on small probability loads such as explosions.In this paper,a 50-m-high concrete gravity dam is used as a prototype dam,and a water explosion model test of a 2.5-m-high concrete gravity dam is designed.The water pressure and the acceleration response of the dam body in the test are analysed.The pressure characteristics and dynamic response of the dam body are assessed.Taking the dam damage test as an example,a numerical model of concrete gravity dam damage is established,and the damage evolution of the dam body is analysed.By combining experiments and numerical simulations,the damage characteristics of the dam body under the action of different charge water explosions are clarified.The integrity of the dam body is well maintained under the action of a small-quantity water explosion,and the dynamic response of the dam body is mainly caused by the shock wave.Both the shock wave and the bubble pulsation cause the dam body to accelerate,and the peak acceleration of the dam body under the action of the bubble pulsation is only one percent of the peak acceleration of the dam body under the action of the shock wave.When subjected to explosions in large quantities of water,the dam body is seriously damaged.Under the action of a shock wave,the dam body produces a secondary acceleration response,which is generated by an internal interaction after the dam body is damaged.The damage evolution process of the dam body under the action of a large-scale water explosion is analysed,and it is found that the shock wave pressure of the water explosion causes local damage to the dam body facing the explosion.After the peak value of the shock wave,the impulse continues to act on the dam body,causing cumulative damage and damage inside the dam body.
基金financially supported by the China Postdoctoral Science Foundation(Grant No.2023M732997)the National Natural Science Foundation of China(Grant Nos.51890912,52008268)Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering,Hohai University(Grant No.2023007)。
文摘The interaction between soil and marine structures like submarine pipeline/pipe pile/suction caisson is a complicated geotechnical mechanism process.In this study,the interface is discretized into multiple mesoscopic contact elements that are damaged randomly throughout the shearing process due to the natural heterogeneity.The evolution equation of damage variable is developed based on the Weibull function,which is able to cover a rather wide range of distribution shapes by only two parameters,making it applicable for varying scenarios.Accordingly,a statistical damage model is established by incorporating Mohr–Coulomb strength criterion,in which the interfacial residual strength is considered whereby the strain softening behavior can be described.A concept of“semi-softening”characteristic point on shear stress–displacement curve is proposed for effectively modeling the evolution of strain softening.Finally,a series of ring shear tests of the interfaces between fine sea sand and smooth/rough steel surfaces are conducted.The predicted results using the proposed model are compared with experimental data of this study as well as some results from existing literature,indicating that the model has a good performance in modeling the progressive failure and strain softening behavior for various types of soil–structure interfaces.
基金National Science and Technology Major Project,Grant/Award Number:2017YFC1503102National Natural Science Foundation of China,Grant/Award Number:51704144。
文摘The creep characteristics of rock under high crustal stress are of important influence on the long‐term stability of deep rock engineering.To study the creep characteristics and engineering application of sandstone under high crustal stress,this study constructed nonlinear creep damage(NCD)constitutive mode based on the triaxial graded loading‒unloading creep test of sandstone in the Yuezhishan Tunnel.A numerical NCD constitutive model and a breakable lining(BL)model were developed based on FLAC3D and then applied to the stability analysis of the Yuezhishan Tunnel.Based on the creep test results of sandstone,a power function of creep rate and stress level was constructed,by which the long‐term strength was solved.The results show that the long‐term strength of the red sandstone based on the related function of the steady‐state creep rate and stress level is close to the measured stress value in engineering.The NCD model considering damage factors reflects the instantaneous and viscoelastic plasticity deformation characteristics of the red sandstone.The numerical NCD constitutive model and the BL model can reflect surrounding rock deformation characteristics and lining failure characteristics in practical engineering.The research results provide theoretical references for long‐term stability analysis of rock engineering and the deformation control of surrounding rock under high crustal stress.
基金This work was financially supported by the National Natural Science fund of China (No.50274058).
文摘Damage statistical mechanics model of horizontal section height in the top caving was constructed in the paper. The influence factors including supporting pressure, dip angle and characteristic of coal on horizontal section height were analyzed as well. By terms of the practice project analysis, the horizontal section height increases with the increase of dip angle β and thickness of coal seam M. Dip angle of coal seam β has tremendous impact on horizontal section height, while thickness of coal seam M has slight impact. When thickness of coal seam is below 10m, horizontal section height increases sharply. While thickness exceeds 15m, it is not major factor influencing on horizontal section height any long.