To reveal the rock burst mechanism,the stress and failure characteristics of coal-rock strata under different advancing speeds of mining working face were explored by theoretical analysis,simulation,and engineering mo...To reveal the rock burst mechanism,the stress and failure characteristics of coal-rock strata under different advancing speeds of mining working face were explored by theoretical analysis,simulation,and engineering monitoring.The relationship between energy accumulation and release was analyzed,and a reasonable mining speed according to specific projects was recommended.The theoretical analysis shows that as the mining speed increases from 4 to 15 m/d,the rheological coefficient of coal mass ranges from 0.9 to 0.4,and the elastic energy of coal mass accumulation varies from 100 to 900 kJ.Based on the simulation,there is a critical advancing speed,the iteration numbers of simulation are less than 15,000 per mining 10 m coal seam,the overburden structure is obvious,the abutment pressure in coal mass is large,and the accumulated energy is large,which is easy to cause strong rock burst.When the iteration number is greater than 15,000,the static force of coal mass increases slightly,but there is no obvious rock burst.Based on engineering monitoring,the mining speed of a mine is less than 8 m/d,and the periodic weighting distance is about 17 m;as the mining speed is greater than 10 m/d,and the periodic weighting distance is greater than 20 m;as the mining speed is 3-8 m/d,and the range of high stress in surrounding rock is 48 m;as the advancing speed is 8-12 m/d,and the high-stress range in surrounding rock is 80 m.Moreover,as the mining speed is less than 8 cut cycles,the micro seismic energy is less than 10,000 J;as the mining speed is 12 cut cycles,the microseismic energy is about 20,000 J.In summary,the advancing speed is positively correlated with the micro seismic event;as the mining speed increases,the accumulated elastic energy of surrounding rock is greater,which is easy to cause rock burst.The comprehensive analysis indicates the daily advance speed of the mine is not more than 12 cut cycles.展开更多
Ground anchor drilling is a promising technology for investigating the mechanical properties and environmental variability of lunar regolith in low-gravity environments,with minimal demands for reactive cutting.This s...Ground anchor drilling is a promising technology for investigating the mechanical properties and environmental variability of lunar regolith in low-gravity environments,with minimal demands for reactive cutting.This study explores the interaction behavior during ground anchor drilling of lunar regolith by employing a coupled approach that integrates the Material Point Method(MPM)and the ContinuouseDiscontinuous Element Method(CDEM),considering the interactions among numerous particles and blocks.The numerical parameters are calibrated based on experimental penetration resistance data of lunar regolith simulant.The numerical approach effectively captures key mechanical properties of the simulant,such as particle flow and scattering patterns,anchor penetration effects,and disturbance-related ultimate bearing characteristics.Additionally,this study examines the influence of inter-particle friction and compactness on penetration resistance.By combining the Golden Section Search Method(GSSM)with ground anchor drilling simulations,an inverse analysis model for penetration resistance is developed,allowing for the determination of mechanical parameters of the lunar regolith simulant.The feasibility of this parameter inversion method is verified,providing valuable insights for engineering applications in lunar exploration and construction.展开更多
Precipitation strengthening is a pivotal mechanism for enhancing the mechanical properties of low-density alloys.A detailed analysis of microstructural evolution during thermal processing is imperative to thoroughly u...Precipitation strengthening is a pivotal mechanism for enhancing the mechanical properties of low-density alloys.A detailed analysis of microstructural evolution during thermal processing is imperative to thoroughly understand its strengthening behavior.This study employed the Bähr D805L quenching dilatometer system to study the formation,evolution,and impact on the contribution of nano-precipitates on the mechanical behavior of Fe-21Mn-10Al-5Ni-C(nominal composition)low-density alloy during continuous cooling.The study unveiled the precipitation mechanism of nano-particles within the austenite(γ)matrix at cooling rates in the range of 40-0.1℃·s^(-1).Moreover,the addition of Ni in Fe-21Mn-10Al-5Ni-C low-density alloy enhances the atomic size factor,promoting alloy spinodal decomposition and ordering.During slow cooling,B2 phases precipitate along grain boundaries,accompanied by the formation of a precipitation-free zone(PFZ)near the boundaries and the dissolution of some later nucleated small particles.These phenomena are a primary mechanism that suppresses the precipitation of B2 phases within theγmatrix.展开更多
Polymerization-induced self-assembly(PISA)has become one of the most versatile approaches for scalable preparation of linear block copolymer nanoparticles with various morphologies.However,the controlled introduction ...Polymerization-induced self-assembly(PISA)has become one of the most versatile approaches for scalable preparation of linear block copolymer nanoparticles with various morphologies.However,the controlled introduction of branching into the core-forming block and the effect on the morphologies of block copolymer nanoparticles under PISA conditions have rarely been explored.Herein,a series of multifunctional macromolecular chain transfer agents(macro-CTAs)were first synthesized by a two-step green light-activated photoiniferter polymerization using two types of chain transfer monomers(CTMs).These macro-CTAs were then used to mediate reversible addition-fragmentation chain transfer(RAFT)dispersion polymerization of styrene(St)to prepare block copolymers with different core-forming block structures and the assemblies.The effect of the core-forming block structure on the morphology of block copolymer nanoparticles was investigated in detail.Transmission electron microscopy(TEM)analysis indicated that the brush-like core-forming block structure facilitated the formation of higher-order morphologies,while the branched core-forming block structure favored the formation of lower-order morphologies.Moreover,it was found that using macroCTAs with a shorter length also promoted the formation of higher-order morphologies.Finally,structures of block copolymers and the assemblies were further controlled by changing the structure of macro-CTA or using a binary mixture of two different macro-CTAs.We expect that this work not only sheds light on the synthesis of block copolymer nanoparticles but also provide important mechanistic insights into PISA of nonlinear block copolymers.展开更多
Continuum-based discrete element method(CDEM)is an explicit numerical method used for simulation of progressive failure of geological body.To improve the efficiency of contact detection and simplify the calculation st...Continuum-based discrete element method(CDEM)is an explicit numerical method used for simulation of progressive failure of geological body.To improve the efficiency of contact detection and simplify the calculation steps for contact forces,semi-spring and semi-edge are introduced in calculation.Semispring is derived from block vertex,and formed by indenting the block vertex into each face(24semisprings for a hexahedral element).The formation process of semi-edge is the same as that of semi-spring(24semi-edges for a hexahedral element).Based on the semi-springs and semi-edges,a new type of combined contact model is presented.According to this model,six contact types could be reduced to two,i.e.the semi-spring target face contact and semi-edge target edge contact.By the combined model,the contact force could be calculated directly(the information of contact type is not necessary),and the failure judgment could be executed in a straightforward way(each semi-spring and semi-edge own their characteristic areas).The algorithm has been successfully programmed in C++program.Some simple numerical cases are presented to show the validity and accuracy of this model.Finally,the failure mode,sliding distance and critical friction angle of Jiweishan landslide are studied with the combined model.展开更多
One of the main issues in environmental hydraulics is pollutant release from sediments.For instance,the strong affinity between phosphorus and sediment permits most of the phosphorus to be adsorbed on the surface of t...One of the main issues in environmental hydraulics is pollutant release from sediments.For instance,the strong affinity between phosphorus and sediment permits most of the phosphorus to be adsorbed on the surface of the sediment particles in rivers or lakes.Post sediment resuspension,phosphorus is desorbed from the sediment to the overlying water.The release of phosphorus from the resuspended sediment is an important process in the secondary pollution of water.Herein,a coupled mechanical model of the overlying water,sediment,and pollutant was established based on the experimentally gathered data.Two types of sediment with different adsorption and desorption characteristics were selected to simulate the process of sediment resuspension and phosphorus release under different hydrodynamic conditions.The simulation results were subsequently used to analyze the relationship between the flow field characteristics and phosphorus concentration,from which the relationships between velocity,particle volume fraction,turbulent kinetic energy,total phosphorus concentration,desorbed phosphorus concentration,and time were elucidated.Based on the results,phosphorus is rapidly released into the overlying water from the resuspended sediment,and it reaches a peak value in a short duration.Unlike the release process of non-adsorption pollutants,hydrodynamic conditions and sediment properties play a crucial role in the phosphorus release process.The turbulent kinetic energy rapidly increases with the flow velocity,whereas the desorbed phosphorus concentration exhibits a certain relationship with the particle volume fraction and turbulent kinetic energy.In particular,the turbulent kinetic energy increases the desorbed phosphorus concentration per unit time.Additionally,the time taken by the total phosphorus concentration to attain its peak value is closely related to the characteristics of the flow field,whereas the amount of phosphorus is closely related to sediment properties.Post sediment resuspension,the release of phosphorus shows the characteristics of a centralized and massive release,which suggests that the total phosphorus concentration in the overlying water would change in a short duration and cause secondary pollution in the water environment.展开更多
The aboveground biomass(AGB)of shrubs and small trees is the main component for the productivity and carbon storage of understory vegetation in subtropical secondary forests.However,few allometric models exist to accu...The aboveground biomass(AGB)of shrubs and small trees is the main component for the productivity and carbon storage of understory vegetation in subtropical secondary forests.However,few allometric models exist to accurately evaluate understory biomass.To estimate the AGB of five common shrub(diameter at base<5 cm,<5 m high)and one small tree species(<8 m high,trees’s seedling),206 individuals were harvested and species-specific and multi-species allometric models developed based on four predictors,height(H),stem diameter(D),crown area(Ca),and wood density(ρ).As expected,the six species possessed greater biomass in their stems compared with branches,with the lowest biomass in the leaves.Species-specific allometric models that employed stem diameter and the combined variables of D~2H andρDH as predictors accurately estimated the components and total AGB,with R^(2) values from 0.602 and 0.971.A multi-species shrub allometric model revealed that wood density×diameter×height(ρDH)was the best predictor,with R^(2) values ranging from between 0.81 and 0.89 for the components and total AGB,respectively.These results indicated that height(H)and diameter(D)were effective predictors for the models to estimate the AGB of the six species,and the introduction of wood density(ρ)improved their accuracy.The optimal models selected in this study could be applied to estimate the biomass of shrubs and small trees in subtropical regions.展开更多
This paper focuses on the relationship between the microstructure and tensile properties of Fe-Mn-Al-C low-density high-strength steel processes by hot-rolling and air-cooling process. The microstructure analysis reve...This paper focuses on the relationship between the microstructure and tensile properties of Fe-Mn-Al-C low-density high-strength steel processes by hot-rolling and air-cooling process. The microstructure analysis reveals that the combination of hot-rolling and air-cooling results in the formation of heterogeneous structures comprising different-sized γ and B_(2) phases in the low-density steel with the addition of nickel (Ni). The addition of Ni promotes the formation of the B_(2) phase and induces the pinning of B_(2) phase particles at the γ grain boundaries. This pinning effect effectively hinders the growth of the γ grains, leading to grain refinement. The tensile test results demonstrate that LDS-5Ni (low-density steel, LDS) exhibits excellent high strength and ductility combination, e.g., a tensile strength of 1535 MPa, yield strength of 1482 MPa, and elongation of 23.3%. These remarkable mechanical properties are primarily attributed to the combined strengthening contributions of grain refinement and duplex nano-sized second-phase precipitation hardening.展开更多
To simulate the progressive failure of slope,a block particle coupled model is introduced.Particle oriented cell mapping(POCM)algorithm is used to enhance the search efficiency,and particlepoint,particle-edge,particle...To simulate the progressive failure of slope,a block particle coupled model is introduced.Particle oriented cell mapping(POCM)algorithm is used to enhance the search efficiency,and particlepoint,particle-edge,particle-face contact detecting method is adopted to establish contact pair between particles and blocks precisely.Strain softening Mohr Coulomb model with tensile cutoff is adopted for blocks,and brittle Mohr Coulomb model is used for particles.The particle-block replacement approach is used to describe the fracture and fragmentation process of continuum media.Once the cohesion or tensile strength of one block reaches zero,the block will be deleted,and particles are generated at the same place with all information inherited from the deleted block.Some numerical cases related to landslides demonstrate the precision and rationality of the coupled model.展开更多
Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle in...Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle induced areas.It is a scenario of strong coupled thermomechanical process involving discontinuity behaviours of rocks.In this context,a numerical model was proposed to investigate the thermal cracking of rocks,in a framework of the continuous-discontinuous element method(CDEM)for efficiently capturing the initiation and propagation of multiple cracks.A simplex integration strategy was adopted to account for the influences of temperature-dependent material properties.Several benchmark tests were considered and the obtained results were compared with analytical solutions and numerical results from the literature.The results show that the fracture degree of the cases when considering temperature-dependent material parameters had 10%differences approximately compared with the cases with constant parameters.展开更多
基金supported by the National Natural Science Foundation of China(52174109)Program for Innovative Research Team(in Science and Technology)in University of Henan Province(22IRTSTHN005)+1 种基金Key Research and Development Project of Henan Province(242102240029)Key Research Project of Institutions of Higher Education in Henan Province(24A580001).
文摘To reveal the rock burst mechanism,the stress and failure characteristics of coal-rock strata under different advancing speeds of mining working face were explored by theoretical analysis,simulation,and engineering monitoring.The relationship between energy accumulation and release was analyzed,and a reasonable mining speed according to specific projects was recommended.The theoretical analysis shows that as the mining speed increases from 4 to 15 m/d,the rheological coefficient of coal mass ranges from 0.9 to 0.4,and the elastic energy of coal mass accumulation varies from 100 to 900 kJ.Based on the simulation,there is a critical advancing speed,the iteration numbers of simulation are less than 15,000 per mining 10 m coal seam,the overburden structure is obvious,the abutment pressure in coal mass is large,and the accumulated energy is large,which is easy to cause strong rock burst.When the iteration number is greater than 15,000,the static force of coal mass increases slightly,but there is no obvious rock burst.Based on engineering monitoring,the mining speed of a mine is less than 8 m/d,and the periodic weighting distance is about 17 m;as the mining speed is greater than 10 m/d,and the periodic weighting distance is greater than 20 m;as the mining speed is 3-8 m/d,and the range of high stress in surrounding rock is 48 m;as the advancing speed is 8-12 m/d,and the high-stress range in surrounding rock is 80 m.Moreover,as the mining speed is less than 8 cut cycles,the micro seismic energy is less than 10,000 J;as the mining speed is 12 cut cycles,the microseismic energy is about 20,000 J.In summary,the advancing speed is positively correlated with the micro seismic event;as the mining speed increases,the accumulated elastic energy of surrounding rock is greater,which is easy to cause rock burst.The comprehensive analysis indicates the daily advance speed of the mine is not more than 12 cut cycles.
基金financial support from the National Natural Science Foundation of China(Grant Nos.52178324,12102059,and 12472207).
文摘Ground anchor drilling is a promising technology for investigating the mechanical properties and environmental variability of lunar regolith in low-gravity environments,with minimal demands for reactive cutting.This study explores the interaction behavior during ground anchor drilling of lunar regolith by employing a coupled approach that integrates the Material Point Method(MPM)and the ContinuouseDiscontinuous Element Method(CDEM),considering the interactions among numerous particles and blocks.The numerical parameters are calibrated based on experimental penetration resistance data of lunar regolith simulant.The numerical approach effectively captures key mechanical properties of the simulant,such as particle flow and scattering patterns,anchor penetration effects,and disturbance-related ultimate bearing characteristics.Additionally,this study examines the influence of inter-particle friction and compactness on penetration resistance.By combining the Golden Section Search Method(GSSM)with ground anchor drilling simulations,an inverse analysis model for penetration resistance is developed,allowing for the determination of mechanical parameters of the lunar regolith simulant.The feasibility of this parameter inversion method is verified,providing valuable insights for engineering applications in lunar exploration and construction.
基金supported by the CNPC Scientific Research and Technology Development Project(Nos.2021ZZ03,2023ZZ11,2022DQ03-02 and 2020B-4020).
文摘Precipitation strengthening is a pivotal mechanism for enhancing the mechanical properties of low-density alloys.A detailed analysis of microstructural evolution during thermal processing is imperative to thoroughly understand its strengthening behavior.This study employed the Bähr D805L quenching dilatometer system to study the formation,evolution,and impact on the contribution of nano-precipitates on the mechanical behavior of Fe-21Mn-10Al-5Ni-C(nominal composition)low-density alloy during continuous cooling.The study unveiled the precipitation mechanism of nano-particles within the austenite(γ)matrix at cooling rates in the range of 40-0.1℃·s^(-1).Moreover,the addition of Ni in Fe-21Mn-10Al-5Ni-C low-density alloy enhances the atomic size factor,promoting alloy spinodal decomposition and ordering.During slow cooling,B2 phases precipitate along grain boundaries,accompanied by the formation of a precipitation-free zone(PFZ)near the boundaries and the dissolution of some later nucleated small particles.These phenomena are a primary mechanism that suppresses the precipitation of B2 phases within theγmatrix.
基金financially supported by the National Natural Science Foundation of China(Nos.22171055 and 52222301)the Guangdong Natural Science Foundation for Distinguished Young Scholar(No.2022B1515020078)the Science and Technology Program of Guangzhou(No.2024A04J2821)。
文摘Polymerization-induced self-assembly(PISA)has become one of the most versatile approaches for scalable preparation of linear block copolymer nanoparticles with various morphologies.However,the controlled introduction of branching into the core-forming block and the effect on the morphologies of block copolymer nanoparticles under PISA conditions have rarely been explored.Herein,a series of multifunctional macromolecular chain transfer agents(macro-CTAs)were first synthesized by a two-step green light-activated photoiniferter polymerization using two types of chain transfer monomers(CTMs).These macro-CTAs were then used to mediate reversible addition-fragmentation chain transfer(RAFT)dispersion polymerization of styrene(St)to prepare block copolymers with different core-forming block structures and the assemblies.The effect of the core-forming block structure on the morphology of block copolymer nanoparticles was investigated in detail.Transmission electron microscopy(TEM)analysis indicated that the brush-like core-forming block structure facilitated the formation of higher-order morphologies,while the branched core-forming block structure favored the formation of lower-order morphologies.Moreover,it was found that using macroCTAs with a shorter length also promoted the formation of higher-order morphologies.Finally,structures of block copolymers and the assemblies were further controlled by changing the structure of macro-CTA or using a binary mixture of two different macro-CTAs.We expect that this work not only sheds light on the synthesis of block copolymer nanoparticles but also provide important mechanistic insights into PISA of nonlinear block copolymers.
基金the National Basic Research Program of the Ministry of Science and Technology of China (Grant No. 2010CB731506)the National Key Technology Research and Development Program of the Ministry of Science and Technology of China (Grant No. 2012BAK10B01)the Youth Science Fund of National Natural Science Foundation of China (Grant No. 11302230)
文摘Continuum-based discrete element method(CDEM)is an explicit numerical method used for simulation of progressive failure of geological body.To improve the efficiency of contact detection and simplify the calculation steps for contact forces,semi-spring and semi-edge are introduced in calculation.Semispring is derived from block vertex,and formed by indenting the block vertex into each face(24semisprings for a hexahedral element).The formation process of semi-edge is the same as that of semi-spring(24semi-edges for a hexahedral element).Based on the semi-springs and semi-edges,a new type of combined contact model is presented.According to this model,six contact types could be reduced to two,i.e.the semi-spring target face contact and semi-edge target edge contact.By the combined model,the contact force could be calculated directly(the information of contact type is not necessary),and the failure judgment could be executed in a straightforward way(each semi-spring and semi-edge own their characteristic areas).The algorithm has been successfully programmed in C++program.Some simple numerical cases are presented to show the validity and accuracy of this model.Finally,the failure mode,sliding distance and critical friction angle of Jiweishan landslide are studied with the combined model.
基金This work was financially supported by the Strategic Priority Research Program of the National Key R&D Program of China(Grants 2018YFC1505500 and 2018YFC1505504)the Chinese Academy of Science(Grant XDB10030303)and the National Natural Science Foundation of China(NSFC)(Grants 11802313 and 11872117).
文摘One of the main issues in environmental hydraulics is pollutant release from sediments.For instance,the strong affinity between phosphorus and sediment permits most of the phosphorus to be adsorbed on the surface of the sediment particles in rivers or lakes.Post sediment resuspension,phosphorus is desorbed from the sediment to the overlying water.The release of phosphorus from the resuspended sediment is an important process in the secondary pollution of water.Herein,a coupled mechanical model of the overlying water,sediment,and pollutant was established based on the experimentally gathered data.Two types of sediment with different adsorption and desorption characteristics were selected to simulate the process of sediment resuspension and phosphorus release under different hydrodynamic conditions.The simulation results were subsequently used to analyze the relationship between the flow field characteristics and phosphorus concentration,from which the relationships between velocity,particle volume fraction,turbulent kinetic energy,total phosphorus concentration,desorbed phosphorus concentration,and time were elucidated.Based on the results,phosphorus is rapidly released into the overlying water from the resuspended sediment,and it reaches a peak value in a short duration.Unlike the release process of non-adsorption pollutants,hydrodynamic conditions and sediment properties play a crucial role in the phosphorus release process.The turbulent kinetic energy rapidly increases with the flow velocity,whereas the desorbed phosphorus concentration exhibits a certain relationship with the particle volume fraction and turbulent kinetic energy.In particular,the turbulent kinetic energy increases the desorbed phosphorus concentration per unit time.Additionally,the time taken by the total phosphorus concentration to attain its peak value is closely related to the characteristics of the flow field,whereas the amount of phosphorus is closely related to sediment properties.Post sediment resuspension,the release of phosphorus shows the characteristics of a centralized and massive release,which suggests that the total phosphorus concentration in the overlying water would change in a short duration and cause secondary pollution in the water environment.
基金supported by the Special Major Science and Technology Project of Anhui Province(S202103b06020066)the 2020 Annual Graduate Innovation Fund of Anhui Agricultural University(2020YSJ-21)。
文摘The aboveground biomass(AGB)of shrubs and small trees is the main component for the productivity and carbon storage of understory vegetation in subtropical secondary forests.However,few allometric models exist to accurately evaluate understory biomass.To estimate the AGB of five common shrub(diameter at base<5 cm,<5 m high)and one small tree species(<8 m high,trees’s seedling),206 individuals were harvested and species-specific and multi-species allometric models developed based on four predictors,height(H),stem diameter(D),crown area(Ca),and wood density(ρ).As expected,the six species possessed greater biomass in their stems compared with branches,with the lowest biomass in the leaves.Species-specific allometric models that employed stem diameter and the combined variables of D~2H andρDH as predictors accurately estimated the components and total AGB,with R^(2) values from 0.602 and 0.971.A multi-species shrub allometric model revealed that wood density×diameter×height(ρDH)was the best predictor,with R^(2) values ranging from between 0.81 and 0.89 for the components and total AGB,respectively.These results indicated that height(H)and diameter(D)were effective predictors for the models to estimate the AGB of the six species,and the introduction of wood density(ρ)improved their accuracy.The optimal models selected in this study could be applied to estimate the biomass of shrubs and small trees in subtropical regions.
基金funding from the National Key Technologies Research and Development Program of China(No.2021YFB3703500)the National Natural Science Foundation of China(No.51771014).
文摘This paper focuses on the relationship between the microstructure and tensile properties of Fe-Mn-Al-C low-density high-strength steel processes by hot-rolling and air-cooling process. The microstructure analysis reveals that the combination of hot-rolling and air-cooling results in the formation of heterogeneous structures comprising different-sized γ and B_(2) phases in the low-density steel with the addition of nickel (Ni). The addition of Ni promotes the formation of the B_(2) phase and induces the pinning of B_(2) phase particles at the γ grain boundaries. This pinning effect effectively hinders the growth of the γ grains, leading to grain refinement. The tensile test results demonstrate that LDS-5Ni (low-density steel, LDS) exhibits excellent high strength and ductility combination, e.g., a tensile strength of 1535 MPa, yield strength of 1482 MPa, and elongation of 23.3%. These remarkable mechanical properties are primarily attributed to the combined strengthening contributions of grain refinement and duplex nano-sized second-phase precipitation hardening.
基金supported by the National Key Research and Development Project of China,the Ministry of Science and Technology of China(Grant 2018YFC1505504)。
文摘To simulate the progressive failure of slope,a block particle coupled model is introduced.Particle oriented cell mapping(POCM)algorithm is used to enhance the search efficiency,and particlepoint,particle-edge,particle-face contact detecting method is adopted to establish contact pair between particles and blocks precisely.Strain softening Mohr Coulomb model with tensile cutoff is adopted for blocks,and brittle Mohr Coulomb model is used for particles.The particle-block replacement approach is used to describe the fracture and fragmentation process of continuum media.Once the cohesion or tensile strength of one block reaches zero,the block will be deleted,and particles are generated at the same place with all information inherited from the deleted block.Some numerical cases related to landslides demonstrate the precision and rationality of the coupled model.
基金the financial support from the Natural Science Foundation of Hebei Province(Grant No.E2020050012)the National Natural Science Foundation of China(NSFC)(Grant No.52178324)the National Key Research and Development Project of China,the Ministry of Science and Technology of China(Grant No.2018YFC1505504).
文摘Thermal cracking of rocks can significantly affect the durability of underground structures in engineering practices such as geothermal energy extraction,storage of nuclear waste and tunnelling in freezeethaw cycle induced areas.It is a scenario of strong coupled thermomechanical process involving discontinuity behaviours of rocks.In this context,a numerical model was proposed to investigate the thermal cracking of rocks,in a framework of the continuous-discontinuous element method(CDEM)for efficiently capturing the initiation and propagation of multiple cracks.A simplex integration strategy was adopted to account for the influences of temperature-dependent material properties.Several benchmark tests were considered and the obtained results were compared with analytical solutions and numerical results from the literature.The results show that the fracture degree of the cases when considering temperature-dependent material parameters had 10%differences approximately compared with the cases with constant parameters.
