Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research...Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.展开更多
The particle morphological properties,such as sphericity,concavity and convexity,of a granular assembly significantly affect its macroscopic and microscopic compressive behaviors under isotropic loading condition.Howe...The particle morphological properties,such as sphericity,concavity and convexity,of a granular assembly significantly affect its macroscopic and microscopic compressive behaviors under isotropic loading condition.However,limited studies on investigating the microscopic behavior of the granular assembly with real particle shapes under isotropic compression were reported.In this study,X-ray computed tomography(mCT)and discrete element modeling(DEM)were utilized to investigate isotropic compression behavior of the granular assembly with regard to the particle morphological properties,such as particle sphericity,concavity and interparticle frictions.The mCT was first used to extract the particle morphological parameters and then the DEM was utilized to numerically investigate the influences of the particle morphological properties on the isotropic compression behavior.The image reconstruction from mCT images indicated that the presented particle quantification algorithm was robust,and the presented microscopic analysis via the DEM simulation demonstrated that the particle surface concavity significantly affected the isotropic compression behavior.The observations of the particle connectivity and local void ratio distribution also provided insights into the granular assembly under isotropic compression.Results found that the particle concavity and interparticle friction influenced the most of the isotropic compression behavior of the granular assemblies.展开更多
Seafloor irregularities influence rupture behavior along the subducting slab and in the overriding plate,thus affecting earthquake cycles.Whether seafloor irregularities increase the likelihood of large earthquakes in...Seafloor irregularities influence rupture behavior along the subducting slab and in the overriding plate,thus affecting earthquake cycles.Whether seafloor irregularities increase the likelihood of large earthquakes in a subduction zone remains contested,partially due to focus put either on fault development or on rupture pattern.Here,we simulate a subducting slab with a seafloor irregularity and the resulting deformation pattern of the overriding plate using the discrete element method.Our simulations illustrate the rupture along three major fault systems:megathrust,splay and backthrust faults.Our results show different rupture dimensions of earthquake events varying from tens to ca.140 km.Our results suggest that the recurrence interval of megathrust events with rupture length of ca.100 km is ca.140 years,which is overall comparable to the paleoseismic records at the Mentawai area of the Sumatran zone.We further propose the coseismic slip amounts decrease and interseismic slip amounts increase from the surface downwards gradually.We conclude that the presence of seafloor irregularities significantly affects rupture events along the slab as well as fault patterns in the overriding plate.展开更多
To investigate migration and evolution rules of coarse aggregates in the static compaction process, an algorithm of generating digital coarse aggregates that can reflect real morphology（ such as shape, size and fract...To investigate migration and evolution rules of coarse aggregates in the static compaction process, an algorithm of generating digital coarse aggregates that can reflect real morphology（ such as shape, size and fracture surface） of aggregate particles, is represented by polyhedral particles based on the discrete element method（ DEM）. A digital specimen comprised of aggregates and air voids is developed. In addition,a static compaction model consisting of a digital specimen and three plates is constructed and a series of evaluation indices such as mean contact force σMCF, wall stress in direction of zcoordinate σWSZZ, porosity and coordination numbers are presented to investigate the motion rules of coarse aggregates at different compaction displacements of 7. 5, 15 and 30 mm. The three-dimensional static compaction model is also verified with laboratory measurements. The results indicate that the compaction displacements are positively related to σMCF and σWSZZ, which increase gradually with the increase in iterative steps. When the compaction proceeds, the digital specimen porosity decreases, but the coordination number increases. The variation ranges of these four indices are different at different compaction displacements. This study provides a method to analyze the compaction mechanism of particle materials such as asphalt mixture and graded broken stone.展开更多
The acoustic emission （AE） features in rock fracture are simulated numerically with discrete element model （DEM）. The specimen is constructed by using spherical particles bonded via the parallel bond model. As a r...The acoustic emission （AE） features in rock fracture are simulated numerically with discrete element model （DEM）. The specimen is constructed by using spherical particles bonded via the parallel bond model. As a result of the heterogeneity in rock specimen, the failure criterion of bonded particle is coupled by the shear and tensile strengths, which follow a normal probability distribution. The Kaiser effect is simulated in the fracture process, for a cubic rock specimen under uniaxial compression with a constant rate. The AE number is estimated with breakages of bonded particles using a pair of parameters, in the temporal and spatial scale, respectively. It is found that the AE numbers and the elastic energy release curves coincide. The range for the Kaiser effect from the AE number and the elastic energy release are the same. Furthermore, the frequency-magnitude relation of the AE number shows that the value of B determined with DEM is consistent with the experimental data.展开更多
Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensiona...Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensional(3D)benchmark model of a laboratory experiment on graywacke to examine the dynamic hydraulic fracturing process under a polyaxial stress state.In the numerical model,injection pressures after breakdown(postbreakdown)are varied to study the impact on fracture growth.The fluid pressure front and crack front are identified in the numerical model to analyze the dynamic relationship between fluid diffusion and fracture propagation.Following the hydraulic fracturing test,the polyaxial stresses are rotated to investigate the influence of the stress field rotation on the fracture slip behavior and permeability.The results show that fracture propagation guides fluid diffusion under a high postbreakdown injection pressure.The crack front runs ahead of the fluid pressure front.Under a low postbreakdown injection pressure,the fluid pressure front gradually reaches the crack front,and fluid diffusion is the main driving factor of fracture propagation.Under polyaxial stress conditions,fluid injection not only opens tensile fractures but also induces hydroshearing.When the polyaxial stress is rotated,the fracture slip direction of a fully extended fracture is consistent with the shear stress direction.The fracture slip direction of a partly extended fracture is influenced by the increase in shear stress.Normal stress affects the permeability evolution by changing the average mechanical aperture.Shear stress can induce shearing and sliding on the fracture plane,thereby increasing permeability.展开更多
Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft...Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.展开更多
Bulking characteristics of gangue are of great significance for the stability of goafs in mining overburden in the caving zones.In this paper,a particle discrete element method with clusters to represent gangue was ad...Bulking characteristics of gangue are of great significance for the stability of goafs in mining overburden in the caving zones.In this paper,a particle discrete element method with clusters to represent gangue was adopted to explore the bulking coefficient time effect of the broken rock in the caving zone under three-dimensional triaxial compression condition.The phenomena of stress corrosion,deformation,and failure of rock blocks were simulated in the numerical model.Meanwhile,a new criterion of rock fragments damage was put forward.It was found that the broken rock has obvious viscoelastic properties.A new equation based on the Burgers creep model was proposed to predict the bulking coefficient of broken rock.A deformation characteristic parameter of the prediction equation was analyzed,which can be set as a fixed value in the mid-and long-term prediction of the bulking coefficient.There are quadratic function relationships between the deformation characteristic parameter value and Talbot gradation index,axial pressure and confining pressure.展开更多
The lateral resistance of sleeper plays an important role in ensuring the stability of a railway track, which may change in the operation of railway, due to the fouling in the ballast bed. In this work, discrete eleme...The lateral resistance of sleeper plays an important role in ensuring the stability of a railway track, which may change in the operation of railway, due to the fouling in the ballast bed. In this work, discrete element method was adopted to investigate the effect of fouling on the lateral resistance of sleeper. The shape information of ballast was captured by method of three-dimensional vision reconstruction. In order to calibrate the mechanical parameters and verify the models, a lateral resistance field test was carried out by using a custom-made device. The contact force distributions in the different parts of sleeper as well as the interaction between ballast and sleeper were discussed in depth. The results show that fouling of ballast bed evidently reduces the lateral resistance of sleeper and the decreasing degree is also related to the fouled position of ballast bed, in the order of shoulder > bottom > side.Therefore, the effect of fouling, especially the fouling in the ballast shoulder, on the lateral resistance of sleeper, should be taken into account in ballast track maintenance work.展开更多
The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment ...The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment interaction simulation was carried out,the accuracy of the particle interaction parameters will directly affect the simulation results.This study proposed a method to systematically calibrate the interaction parameters between deep-sea sediment and grouser through the combination of experiment and simulation.The uniaxial compression test and macro adhesion test and corresponding discrete element numerical simulation were carried out,modifying the contact parameters until the simulation results are close to the experimental results.Then the micro-parameters of the JKR adhesion contact model were back calibrated with the test results,and the contact parameters between soil particle-soil particle and soil particle-metal are calibrated.Besides,the adhesion test shows that the adhesion forces were ranked in the order of 5052<STi80<TA2<TC4 under the same surface roughness,which indicates the aluminum alloy 5052 has the best anti-adhesion performance.The relationship between surface adhesion force and microscopic contact parameters was studied by discrete element numerical simulation,and the result shows that the coefficient of static friction and the coefficient of rolling friction has little effect on adhesion force.While it is mainly affected by the coefficient of restitution and surface energy,the surface adhesion force tends to decrease with the increase of the coefficient of restitution and increase with the growth of surface energy.The obtained parameters of soil particle to soil particle and soil particle to metal affecting the adhesion can contribute to the design optimization for the grouser of mining machines to decrease surface adhesion and enhance its movability and mining efficiency.展开更多
The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fractur...The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.展开更多
Wellbore instability is a key problem restricting efficient production of coal-bed methane. In order to perform thorough and systematic research regarding coal-bed wellbore stability problems, a new discrete element m...Wellbore instability is a key problem restricting efficient production of coal-bed methane. In order to perform thorough and systematic research regarding coal-bed wellbore stability problems, a new discrete element model which fully considers the features of cleat coal-beds is established based on the Kirsch equation. With this model, the safe pipe tripping speed, drilling fluid density window and coal- bed collapse/fracture pressure are determined; in addition, the relationships between pipe tripping speed and pipe size, cleat size, etc. and wellbore stability are analyzed in the coal-bed drilling and pipe tripping processes. The case studies show the following results： the wellbore collapses （collapse pressure： 4.33 MPa） or fractures （fracture pressure： 12.7 MPa） in certain directions as a result of swab or surge pressure when the pipe tripping speed is higher than a certain value; the cleat face size has a great influence on wellbore stability, and if the drilling fluid pressure is too low, the wellbore is prone to collapse when the ratio of the face cleat size to butt cleat size is reduced; however, if the drilling fluid pressure is high enough, the butt cleat size has no influence on the wellbore fracture; the factors influencing coal-bed stability include the movement length, pipe size, borehole size.展开更多
To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 ...To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 simulations of one-dimensional compression tests on coarse-grained sand using the three-dimensional(3D)discrete element method(DEM)were conducted to construct a database.In this process,the positions of the particles were randomly altered,and the particle assemblages changed.Interestingly,besides confirming the influence of particle size distribution parameters,the stress-strain curves differed despite an identical gradation size statistic when the particle position varied.Subsequently,the obtained data were partitioned into training,validation,and testing datasets at a 7:2:1 ratio.To convert the DEM model into a multi-dimensional matrix that computers can recognize,the 3D DEM models were first sliced to extract multi-layer two-dimensional(2D)cross-sectional data.Redundant information was then eliminated via gray processing,and the data were stacked to form a new 3D matrix representing the granular soil’s fabric.Subsequently,utilizing the Python language and Pytorch framework,a 3D convolutional neural networks(CNNs)model was developed to establish the relationship between the constrained modulus obtained from DEM simulations and the soil’s fabric.The mean squared error(MSE)function was utilized to assess the loss value during the training process.When the learning rate(LR)fell within the range of 10-5e10-1,and the batch sizes(BSs)were 4,8,16,32,and 64,the loss value stabilized after 100 training epochs in the training and validation dataset.For BS?32 and LR?10-3,the loss reached a minimum.In the testing set,a comparative evaluation of the predicted constrained modulus from the 3D CNNs versus the simulated modulus obtained via DEM reveals a minimum mean absolute percentage error(MAPE)of 4.43%under the optimized condition,demonstrating the accuracy of this approach.Thus,by combining DEM and CNNs,the variation of soil’s mechanical characteristics related to its random fabric would be efficiently evaluated by directly tracking the particle assemblages.展开更多
This paper presents the technological relevance of a concurrent algorithm-based discrete element modelling (DEM) system, HADES. This new system is the successor of SPACE that is limited to spherical grains only. It ca...This paper presents the technological relevance of a concurrent algorithm-based discrete element modelling (DEM) system, HADES. This new system is the successor of SPACE that is limited to spherical grains only. It can realistically simulate the packing of arbitrary-shaped particles up to the fully compacted state. Generation of families of such particles, i.e., generally representing aggregate of fluvial origin and crushed rock, respectively, and the forming way of particulate structure are described. Similarly shaped particles are proposed for simulation of cement paste because of conformity with experimental results obtained by the X-ray tomography method. Technologically relevant territories inside and outside concrete technology are presently explored in this efficient, reliable, and economic way. Some results obtained by this DEM approach are presented.展开更多
An inverse method for parameters identification of discrete element model combined with experiment is proposed.The inverse problem of parameter identification is transmitted to solve an optimization problem by minimiz...An inverse method for parameters identification of discrete element model combined with experiment is proposed.The inverse problem of parameter identification is transmitted to solve an optimization problem by minimizing the distance between the numerical calculations and experiment responses.In this method,the discrete element method is employed as numerical calculator for the forward problem.Then,the orthogonal experiment design with range analysis was used to carry out parameters sensitivity analysis.In addition,to improve the computational efficiency,the approximate model technique is used to replace the actual computational model.The intergeneration projection genetic algorithm(IP-GA)is employed as the optimization algorithm.Consequently,the parameters of the discrete element model are determined.To verify the effectiveness and accuracy of the inverse results,the comparisons of shape deviation experiments with discrete element simulations are provided.It indicates that the effective and reliable discrete element model parameters can be quickly obtained through several sets of experimental data.Hence,this inverse method can be applied more widely to determine the parameters of discrete element model for other materials.展开更多
To improve the application of discrete element models(DEM)to the design of agricultural crushers,in this study a new highly accurate model is elaborated.The model takes into account the fiber structure,porous nature o...To improve the application of discrete element models(DEM)to the design of agricultural crushers,in this study a new highly accurate model is elaborated.The model takes into account the fiber structure,porous nature of the material and the leaf sheath coating structure.Dedicated experimental tests are conducted to determine the required“intrinsic”and basic contact parameters of the considered banana straw materials.A large number of bonding parameters are examined in relation to the particle aggregation model in order to characterize different actual banana straws.Using the particle surface energy contact model,the viscosity characteristics of the crushed material are determined together with the related stacking angle(considered as the main response factor).Through single factor experiment analysis,it is found that when the surface energy is 0.9 J·m-2,the relative error between simulations and physical experiments is 5.288%.展开更多
The high-temperature creep behavior of asphalt mixture was investigated based on micromechanical modeling and virtual test by using three-dimensional discrete element method（DEM）. A user-defined micromechanical mode...The high-temperature creep behavior of asphalt mixture was investigated based on micromechanical modeling and virtual test by using three-dimensional discrete element method（DEM）. A user-defined micromechanical model of asphalt mixture was established after analyzing the irregular shape and gradation of coarse aggregates, the viscoelastic property of asphalt mastic, and the random distribution of air voids within the asphalt mixture. Virtual uniaxial static creep test at 60 ℃ was conducted by using Particle Flow Code in three dimensions（PFC3D） and was validated by laboratory test. Based on virtual creep test, the micromechanical characteristics between aggregates, within asphalt mastic, and between aggregate and asphalt mastic were analyzed for the asphalt mixture. It is proved that the virtual test based on the micromechanical model can efficiently predict the creep deformation of asphalt mixture. And the high-temperature behavior of asphalt mixture was characterized from micromechanical perspective.展开更多
The issues of seepage in calcareous sand foundations and backfillshave a potentially detrimental effect on the stability and safety of superstructures.Simplifying calcareous sand grains as spheres or ellipsoids in num...The issues of seepage in calcareous sand foundations and backfillshave a potentially detrimental effect on the stability and safety of superstructures.Simplifying calcareous sand grains as spheres or ellipsoids in numerical simulations may lead to significantinaccuracies.In this paper,we present a novel intelligence framework based on a deep convolutional generative adversarial network(DCGAN).A DCGAN model was trained using a training dataset comprising 11,625 real particles for the random generation of three-dimensional calcareous sand particles.Subsequently,3800 realistic calcareous sand particles with intra-particle voids were generated.Generative fidelityand validity of the DCGAN model were well verifiedby the consistency of the statistical values of nine morphological parameters of both the training dataset and the generated dataset.Digital calcareous sand columns were obtained through gravitational deposition simulation of the generated particles.Directional seepage simulations were conducted,and the vertical permeability values of the sand columns were found to be in accordance with the objective law.The results demonstrate the potential of the proposed framework for stochastic modeling and multi-scale simulation of the seepage behaviors in calcareous sand foundations and backfills.展开更多
Coriolis effects,encompassing the dilative,compressive,and deflective manifestations,constitute pivotal considerations in the centrifugal modelling of high-speed granular run-out processes.Notably,under the deflective...Coriolis effects,encompassing the dilative,compressive,and deflective manifestations,constitute pivotal considerations in the centrifugal modelling of high-speed granular run-out processes.Notably,under the deflective Coriolis condition,the velocity component parallel to the rotational axis exerts no influence on the magnitude of Coriolis acceleration.This circumstance implies a potential mitigation of the Coriolis force's deflective impact.Regrettably,extant investigations predominantly emphasize the dilative and compressive Coriolis effects,largely neglecting the pragmatic import of the deflective Coriolis condition.In pursuit of this gap,a series of discrete element method(DEM)simulations have been conducted to scrutinize the feasibility of centrifugal modelling for dry granular run-out processes under deflective Coriolis conditions.The findings concerning the deflective Coriolis effect reveal a consistent rise in the run-out distance by 2%–16%,a modest increase in bulk flow velocity of under 4%,and a slight elevation in average flow depth by no more than 25%.These alterations display smaller dependence on the specific testing conditions due to the granular flow undergoing dual deflections in opposing directions.This underscores the significance and utility of the deflective Coriolis condition.Notably,the anticipated reduction in error in predicting the final run-out distance is substantial,potentially reaching a 150%improvement compared to predictions made under the dilative and compressive Coriolis conditions.Therefore,the deflective Coriolis condition is advised when the final run-out distance of the granular flow is the main concern.To mitigate the impact of Coriolis acceleration,a greater initial height of the granular column is recommended,with a height/width ratio exceeding 1,as the basal friction of the granular material plays a crucial role in mitigating the deflective Coriolis effect.For more transverse-uniform flow properties,the width of the granular column should be as large as possible.展开更多
Facing the high demand for faster and heavier freight trains in Australia,researchers and practitioners are endeavouring to develop more innovative and resilient ballasted tracks.In recent years,many studies have been...Facing the high demand for faster and heavier freight trains in Australia,researchers and practitioners are endeavouring to develop more innovative and resilient ballasted tracks.In recent years,many studies have been conducted by the researchers from Transport Research Centre at the University of Technology Sydney(TRC-UTS)to examine the feasibility of incorporating recycled tyre/rubber into rail tracks.This paper reviews three innovative applications using recycled rubber products such as(1)a synthetic energy-absorbing layer for railway subballast using a composite of rubber crumbs and mining byproducts,(2)using rubber intermixed ballast stratum to replace conventional ballast,and(3)installing recycled rubber mat to mitigate ballast degradation under the impact loading.Comprehensive laboratory and field tests as well as numerical modelling have been conducted to examine the performance of rail tracks incorporating these innovative inclusions.The laboratory and field test results and numerical modelling reveal that incorporating these rubber products could increase the energy-absorbing capacity of the track,and mitigate the ballast breakage and settlement significantly,hence increasing the track stability.The research outcomes will facilitate a better understanding of the performance of ballast tracks incorporating these resilient waste tyre materials while promoting more economical and environmentally sustainable tracks for greater passenger comfort and increased safety.展开更多
基金support of the National Natural Science Foundation of China(51604275)the Key Laboratory of Urban Under Ground Engineering of Ministry of Education(TUE2018-01)+1 种基金Yue Qi Young Scholar Project of China University of Mining&Technology,Beijingthe Fundamental Research Funds for the Central Universities(2016QL02).
文摘Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.
基金the Universidad Nacional de San Agustín(UNSA)through the joint Center for Mining Sustainability with the Colorado School of Mines is highly acknowledged.
文摘The particle morphological properties,such as sphericity,concavity and convexity,of a granular assembly significantly affect its macroscopic and microscopic compressive behaviors under isotropic loading condition.However,limited studies on investigating the microscopic behavior of the granular assembly with real particle shapes under isotropic compression were reported.In this study,X-ray computed tomography(mCT)and discrete element modeling(DEM)were utilized to investigate isotropic compression behavior of the granular assembly with regard to the particle morphological properties,such as particle sphericity,concavity and interparticle frictions.The mCT was first used to extract the particle morphological parameters and then the DEM was utilized to numerically investigate the influences of the particle morphological properties on the isotropic compression behavior.The image reconstruction from mCT images indicated that the presented particle quantification algorithm was robust,and the presented microscopic analysis via the DEM simulation demonstrated that the particle surface concavity significantly affected the isotropic compression behavior.The observations of the particle connectivity and local void ratio distribution also provided insights into the granular assembly under isotropic compression.Results found that the particle concavity and interparticle friction influenced the most of the isotropic compression behavior of the granular assemblies.
基金This work is Earth Observatory of Singapore contribution(No.M4430217.B50.706022)the Ministry of Science and Technology(Grant Nos.MOST 109-2116-M-008-029-MY3,MOST 110-2124-M-002-008,and MOST 110-2634-F-008-008)。
文摘Seafloor irregularities influence rupture behavior along the subducting slab and in the overriding plate,thus affecting earthquake cycles.Whether seafloor irregularities increase the likelihood of large earthquakes in a subduction zone remains contested,partially due to focus put either on fault development or on rupture pattern.Here,we simulate a subducting slab with a seafloor irregularity and the resulting deformation pattern of the overriding plate using the discrete element method.Our simulations illustrate the rupture along three major fault systems:megathrust,splay and backthrust faults.Our results show different rupture dimensions of earthquake events varying from tens to ca.140 km.Our results suggest that the recurrence interval of megathrust events with rupture length of ca.100 km is ca.140 years,which is overall comparable to the paleoseismic records at the Mentawai area of the Sumatran zone.We further propose the coseismic slip amounts decrease and interseismic slip amounts increase from the surface downwards gradually.We conclude that the presence of seafloor irregularities significantly affects rupture events along the slab as well as fault patterns in the overriding plate.
基金The National Natural Science Foundation of China(No.51108081)
文摘To investigate migration and evolution rules of coarse aggregates in the static compaction process, an algorithm of generating digital coarse aggregates that can reflect real morphology（ such as shape, size and fracture surface） of aggregate particles, is represented by polyhedral particles based on the discrete element method（ DEM）. A digital specimen comprised of aggregates and air voids is developed. In addition,a static compaction model consisting of a digital specimen and three plates is constructed and a series of evaluation indices such as mean contact force σMCF, wall stress in direction of zcoordinate σWSZZ, porosity and coordination numbers are presented to investigate the motion rules of coarse aggregates at different compaction displacements of 7. 5, 15 and 30 mm. The three-dimensional static compaction model is also verified with laboratory measurements. The results indicate that the compaction displacements are positively related to σMCF and σWSZZ, which increase gradually with the increase in iterative steps. When the compaction proceeds, the digital specimen porosity decreases, but the coordination number increases. The variation ranges of these four indices are different at different compaction displacements. This study provides a method to analyze the compaction mechanism of particle materials such as asphalt mixture and graded broken stone.
基金supported by the National Basic Research Program of China (2010CB731502)
文摘The acoustic emission （AE） features in rock fracture are simulated numerically with discrete element model （DEM）. The specimen is constructed by using spherical particles bonded via the parallel bond model. As a result of the heterogeneity in rock specimen, the failure criterion of bonded particle is coupled by the shear and tensile strengths, which follow a normal probability distribution. The Kaiser effect is simulated in the fracture process, for a cubic rock specimen under uniaxial compression with a constant rate. The AE number is estimated with breakages of bonded particles using a pair of parameters, in the temporal and spatial scale, respectively. It is found that the AE numbers and the elastic energy release curves coincide. The range for the Kaiser effect from the AE number and the elastic energy release are the same. Furthermore, the frequency-magnitude relation of the AE number shows that the value of B determined with DEM is consistent with the experimental data.
基金supported by the Knowledge Innovation Program of Wuhan-Basic Research (Grant No.2022010801010159)support from the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES (Contract number VH-NG-1516)supported by the Swedish Radiation Safety Authority (Project SSM2020-2758).
文摘Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensional(3D)benchmark model of a laboratory experiment on graywacke to examine the dynamic hydraulic fracturing process under a polyaxial stress state.In the numerical model,injection pressures after breakdown(postbreakdown)are varied to study the impact on fracture growth.The fluid pressure front and crack front are identified in the numerical model to analyze the dynamic relationship between fluid diffusion and fracture propagation.Following the hydraulic fracturing test,the polyaxial stresses are rotated to investigate the influence of the stress field rotation on the fracture slip behavior and permeability.The results show that fracture propagation guides fluid diffusion under a high postbreakdown injection pressure.The crack front runs ahead of the fluid pressure front.Under a low postbreakdown injection pressure,the fluid pressure front gradually reaches the crack front,and fluid diffusion is the main driving factor of fracture propagation.Under polyaxial stress conditions,fluid injection not only opens tensile fractures but also induces hydroshearing.When the polyaxial stress is rotated,the fracture slip direction of a fully extended fracture is consistent with the shear stress direction.The fracture slip direction of a partly extended fracture is influenced by the increase in shear stress.Normal stress affects the permeability evolution by changing the average mechanical aperture.Shear stress can induce shearing and sliding on the fracture plane,thereby increasing permeability.
基金Special Fund of Marine Commonweal Industry under contact Nos 201105016 and 201205007supported by National Marine Environment Forecasting Centrethe National Natural Science Foundation of China under contact No.41176012
文摘Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.
基金This work was supported by the National Natural Science Foundation of China,NSFC(Nos.U1803118 and 51974296)and the China Scholarship Council(CSC)(award to Fanfei Meng for PhD period at Kyushu University).
文摘Bulking characteristics of gangue are of great significance for the stability of goafs in mining overburden in the caving zones.In this paper,a particle discrete element method with clusters to represent gangue was adopted to explore the bulking coefficient time effect of the broken rock in the caving zone under three-dimensional triaxial compression condition.The phenomena of stress corrosion,deformation,and failure of rock blocks were simulated in the numerical model.Meanwhile,a new criterion of rock fragments damage was put forward.It was found that the broken rock has obvious viscoelastic properties.A new equation based on the Burgers creep model was proposed to predict the bulking coefficient of broken rock.A deformation characteristic parameter of the prediction equation was analyzed,which can be set as a fixed value in the mid-and long-term prediction of the bulking coefficient.There are quadratic function relationships between the deformation characteristic parameter value and Talbot gradation index,axial pressure and confining pressure.
基金Projects(U1234211,61472029,51208034)supported by the National Natural Science Foundation of China
文摘The lateral resistance of sleeper plays an important role in ensuring the stability of a railway track, which may change in the operation of railway, due to the fouling in the ballast bed. In this work, discrete element method was adopted to investigate the effect of fouling on the lateral resistance of sleeper. The shape information of ballast was captured by method of three-dimensional vision reconstruction. In order to calibrate the mechanical parameters and verify the models, a lateral resistance field test was carried out by using a custom-made device. The contact force distributions in the different parts of sleeper as well as the interaction between ballast and sleeper were discussed in depth. The results show that fouling of ballast bed evidently reduces the lateral resistance of sleeper and the decreasing degree is also related to the fouled position of ballast bed, in the order of shoulder > bottom > side.Therefore, the effect of fouling, especially the fouling in the ballast shoulder, on the lateral resistance of sleeper, should be taken into account in ballast track maintenance work.
基金Project(12072309)supported by the National Natural Science Foundation of ChinaProject(19B546)supported by the Education Department Foundation of Hunan Province,ChinaProject(2019RS1059)supported by the Hunan Innovative Province Construction Project,China。
文摘The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment interaction simulation was carried out,the accuracy of the particle interaction parameters will directly affect the simulation results.This study proposed a method to systematically calibrate the interaction parameters between deep-sea sediment and grouser through the combination of experiment and simulation.The uniaxial compression test and macro adhesion test and corresponding discrete element numerical simulation were carried out,modifying the contact parameters until the simulation results are close to the experimental results.Then the micro-parameters of the JKR adhesion contact model were back calibrated with the test results,and the contact parameters between soil particle-soil particle and soil particle-metal are calibrated.Besides,the adhesion test shows that the adhesion forces were ranked in the order of 5052<STi80<TA2<TC4 under the same surface roughness,which indicates the aluminum alloy 5052 has the best anti-adhesion performance.The relationship between surface adhesion force and microscopic contact parameters was studied by discrete element numerical simulation,and the result shows that the coefficient of static friction and the coefficient of rolling friction has little effect on adhesion force.While it is mainly affected by the coefficient of restitution and surface energy,the surface adhesion force tends to decrease with the increase of the coefficient of restitution and increase with the growth of surface energy.The obtained parameters of soil particle to soil particle and soil particle to metal affecting the adhesion can contribute to the design optimization for the grouser of mining machines to decrease surface adhesion and enhance its movability and mining efficiency.
文摘The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.
文摘Wellbore instability is a key problem restricting efficient production of coal-bed methane. In order to perform thorough and systematic research regarding coal-bed wellbore stability problems, a new discrete element model which fully considers the features of cleat coal-beds is established based on the Kirsch equation. With this model, the safe pipe tripping speed, drilling fluid density window and coal- bed collapse/fracture pressure are determined; in addition, the relationships between pipe tripping speed and pipe size, cleat size, etc. and wellbore stability are analyzed in the coal-bed drilling and pipe tripping processes. The case studies show the following results： the wellbore collapses （collapse pressure： 4.33 MPa） or fractures （fracture pressure： 12.7 MPa） in certain directions as a result of swab or surge pressure when the pipe tripping speed is higher than a certain value; the cleat face size has a great influence on wellbore stability, and if the drilling fluid pressure is too low, the wellbore is prone to collapse when the ratio of the face cleat size to butt cleat size is reduced; however, if the drilling fluid pressure is high enough, the butt cleat size has no influence on the wellbore fracture; the factors influencing coal-bed stability include the movement length, pipe size, borehole size.
基金supported by the National Key R&D Program of China (Grant No.2022YFC3003401)the National Natural Science Foundation of China (Grant Nos.42041006 and 42377137).
文摘To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 simulations of one-dimensional compression tests on coarse-grained sand using the three-dimensional(3D)discrete element method(DEM)were conducted to construct a database.In this process,the positions of the particles were randomly altered,and the particle assemblages changed.Interestingly,besides confirming the influence of particle size distribution parameters,the stress-strain curves differed despite an identical gradation size statistic when the particle position varied.Subsequently,the obtained data were partitioned into training,validation,and testing datasets at a 7:2:1 ratio.To convert the DEM model into a multi-dimensional matrix that computers can recognize,the 3D DEM models were first sliced to extract multi-layer two-dimensional(2D)cross-sectional data.Redundant information was then eliminated via gray processing,and the data were stacked to form a new 3D matrix representing the granular soil’s fabric.Subsequently,utilizing the Python language and Pytorch framework,a 3D convolutional neural networks(CNNs)model was developed to establish the relationship between the constrained modulus obtained from DEM simulations and the soil’s fabric.The mean squared error(MSE)function was utilized to assess the loss value during the training process.When the learning rate(LR)fell within the range of 10-5e10-1,and the batch sizes(BSs)were 4,8,16,32,and 64,the loss value stabilized after 100 training epochs in the training and validation dataset.For BS?32 and LR?10-3,the loss reached a minimum.In the testing set,a comparative evaluation of the predicted constrained modulus from the 3D CNNs versus the simulated modulus obtained via DEM reveals a minimum mean absolute percentage error(MAPE)of 4.43%under the optimized condition,demonstrating the accuracy of this approach.Thus,by combining DEM and CNNs,the variation of soil’s mechanical characteristics related to its random fabric would be efficiently evaluated by directly tracking the particle assemblages.
文摘This paper presents the technological relevance of a concurrent algorithm-based discrete element modelling (DEM) system, HADES. This new system is the successor of SPACE that is limited to spherical grains only. It can realistically simulate the packing of arbitrary-shaped particles up to the fully compacted state. Generation of families of such particles, i.e., generally representing aggregate of fluvial origin and crushed rock, respectively, and the forming way of particulate structure are described. Similarly shaped particles are proposed for simulation of cement paste because of conformity with experimental results obtained by the X-ray tomography method. Technologically relevant territories inside and outside concrete technology are presently explored in this efficient, reliable, and economic way. Some results obtained by this DEM approach are presented.
基金supported by the National Natural Science Foundation of China(11602212)the Natural Science Foundation of Hunan Province of China(2018JJ3509)supported by the National Natural Science Foundation of China(51605409,11802258,51775468).
文摘An inverse method for parameters identification of discrete element model combined with experiment is proposed.The inverse problem of parameter identification is transmitted to solve an optimization problem by minimizing the distance between the numerical calculations and experiment responses.In this method,the discrete element method is employed as numerical calculator for the forward problem.Then,the orthogonal experiment design with range analysis was used to carry out parameters sensitivity analysis.In addition,to improve the computational efficiency,the approximate model technique is used to replace the actual computational model.The intergeneration projection genetic algorithm(IP-GA)is employed as the optimization algorithm.Consequently,the parameters of the discrete element model are determined.To verify the effectiveness and accuracy of the inverse results,the comparisons of shape deviation experiments with discrete element simulations are provided.It indicates that the effective and reliable discrete element model parameters can be quickly obtained through several sets of experimental data.Hence,this inverse method can be applied more widely to determine the parameters of discrete element model for other materials.
基金Youth Fund of National Natural Science Foundation of China,61801175,JJ,https://www.nsfc.gov.cn/.
文摘To improve the application of discrete element models(DEM)to the design of agricultural crushers,in this study a new highly accurate model is elaborated.The model takes into account the fiber structure,porous nature of the material and the leaf sheath coating structure.Dedicated experimental tests are conducted to determine the required“intrinsic”and basic contact parameters of the considered banana straw materials.A large number of bonding parameters are examined in relation to the particle aggregation model in order to characterize different actual banana straws.Using the particle surface energy contact model,the viscosity characteristics of the crushed material are determined together with the related stacking angle(considered as the main response factor).Through single factor experiment analysis,it is found that when the surface energy is 0.9 J·m-2,the relative error between simulations and physical experiments is 5.288%.
基金Funded by the National Natural Science Foundation of China(No.51378006)the Huoyingdong Foundation of China(No.141076)+1 种基金the Fundamental Research Funds for the Central Universities(No.2242015R30027)the Natural Science Foundation of Jiangsu Province(BK20161421 and BK20140109)
文摘The high-temperature creep behavior of asphalt mixture was investigated based on micromechanical modeling and virtual test by using three-dimensional discrete element method（DEM）. A user-defined micromechanical model of asphalt mixture was established after analyzing the irregular shape and gradation of coarse aggregates, the viscoelastic property of asphalt mastic, and the random distribution of air voids within the asphalt mixture. Virtual uniaxial static creep test at 60 ℃ was conducted by using Particle Flow Code in three dimensions（PFC3D） and was validated by laboratory test. Based on virtual creep test, the micromechanical characteristics between aggregates, within asphalt mastic, and between aggregate and asphalt mastic were analyzed for the asphalt mixture. It is proved that the virtual test based on the micromechanical model can efficiently predict the creep deformation of asphalt mixture. And the high-temperature behavior of asphalt mixture was characterized from micromechanical perspective.
基金financially supported by the National Natural Science Foundation of China(Grant No.42077232)the National Natural Science Foundation for Excellent Young Scholars of China(Grant No.52222110)the Fundamental Research Funds for the Central Universities(Grant No.14380229).
文摘The issues of seepage in calcareous sand foundations and backfillshave a potentially detrimental effect on the stability and safety of superstructures.Simplifying calcareous sand grains as spheres or ellipsoids in numerical simulations may lead to significantinaccuracies.In this paper,we present a novel intelligence framework based on a deep convolutional generative adversarial network(DCGAN).A DCGAN model was trained using a training dataset comprising 11,625 real particles for the random generation of three-dimensional calcareous sand particles.Subsequently,3800 realistic calcareous sand particles with intra-particle voids were generated.Generative fidelityand validity of the DCGAN model were well verifiedby the consistency of the statistical values of nine morphological parameters of both the training dataset and the generated dataset.Digital calcareous sand columns were obtained through gravitational deposition simulation of the generated particles.Directional seepage simulations were conducted,and the vertical permeability values of the sand columns were found to be in accordance with the objective law.The results demonstrate the potential of the proposed framework for stochastic modeling and multi-scale simulation of the seepage behaviors in calcareous sand foundations and backfills.
基金supported by the National Natural Science Foundation of China(Grant Nos.42120104008 and 42307214)the Postdoctoral Fellowship Program of CPSF(Grant No.GZB20230620).
文摘Coriolis effects,encompassing the dilative,compressive,and deflective manifestations,constitute pivotal considerations in the centrifugal modelling of high-speed granular run-out processes.Notably,under the deflective Coriolis condition,the velocity component parallel to the rotational axis exerts no influence on the magnitude of Coriolis acceleration.This circumstance implies a potential mitigation of the Coriolis force's deflective impact.Regrettably,extant investigations predominantly emphasize the dilative and compressive Coriolis effects,largely neglecting the pragmatic import of the deflective Coriolis condition.In pursuit of this gap,a series of discrete element method(DEM)simulations have been conducted to scrutinize the feasibility of centrifugal modelling for dry granular run-out processes under deflective Coriolis conditions.The findings concerning the deflective Coriolis effect reveal a consistent rise in the run-out distance by 2%–16%,a modest increase in bulk flow velocity of under 4%,and a slight elevation in average flow depth by no more than 25%.These alterations display smaller dependence on the specific testing conditions due to the granular flow undergoing dual deflections in opposing directions.This underscores the significance and utility of the deflective Coriolis condition.Notably,the anticipated reduction in error in predicting the final run-out distance is substantial,potentially reaching a 150%improvement compared to predictions made under the dilative and compressive Coriolis conditions.Therefore,the deflective Coriolis condition is advised when the final run-out distance of the granular flow is the main concern.To mitigate the impact of Coriolis acceleration,a greater initial height of the granular column is recommended,with a height/width ratio exceeding 1,as the basal friction of the granular material plays a crucial role in mitigating the deflective Coriolis effect.For more transverse-uniform flow properties,the width of the granular column should be as large as possible.
基金financial support from the Australian Research Council for ARCLP200200915 and ARCDP220102862financial and technical support from industry partners including Sydney Trains,SMEC Australia Pty.
文摘Facing the high demand for faster and heavier freight trains in Australia,researchers and practitioners are endeavouring to develop more innovative and resilient ballasted tracks.In recent years,many studies have been conducted by the researchers from Transport Research Centre at the University of Technology Sydney(TRC-UTS)to examine the feasibility of incorporating recycled tyre/rubber into rail tracks.This paper reviews three innovative applications using recycled rubber products such as(1)a synthetic energy-absorbing layer for railway subballast using a composite of rubber crumbs and mining byproducts,(2)using rubber intermixed ballast stratum to replace conventional ballast,and(3)installing recycled rubber mat to mitigate ballast degradation under the impact loading.Comprehensive laboratory and field tests as well as numerical modelling have been conducted to examine the performance of rail tracks incorporating these innovative inclusions.The laboratory and field test results and numerical modelling reveal that incorporating these rubber products could increase the energy-absorbing capacity of the track,and mitigate the ballast breakage and settlement significantly,hence increasing the track stability.The research outcomes will facilitate a better understanding of the performance of ballast tracks incorporating these resilient waste tyre materials while promoting more economical and environmentally sustainable tracks for greater passenger comfort and increased safety.
