The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete elem...The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.展开更多
A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)t...A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.展开更多
Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model sl...Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.展开更多
This study presents the first step of a research project that aims at using a three-dimensional （3D） hybridfinite-discrete element method （FDEM） to investigate the development of an excavation damaged zone（EDZ） ...This study presents the first step of a research project that aims at using a three-dimensional （3D） hybridfinite-discrete element method （FDEM） to investigate the development of an excavation damaged zone（EDZ） around tunnels in a clay shale formation known as Opalinus Clay. The 3D FDEM was first calibratedagainst standard laboratory experiments, including Brazilian disc test and uniaxial compression test. Theeffect of increasing confining pressure on the mechanical response and fracture propagation of the rockwas quantified under triaxial compression tests. Polyaxial （or true triaxial） simulations highlighted theeffect of the intermediate principal stress （s2） on fracture directions in the model： as the intermediateprincipal stress increased, fractures tended to align in the direction parallel to the plane defined by themajor and intermediate principal stresses. The peak strength was also shown to vary with changing s2. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
We simulated three-dimensional heat transfer inside a horizontal rotating drum using the discrete element method and a thermal conduction model.The aim was to determine the effect of end-wall heating on thermal behavi...We simulated three-dimensional heat transfer inside a horizontal rotating drum using the discrete element method and a thermal conduction model.The aim was to determine the effect of end-wall heating on thermal behavior of a granular bed.The simulation showed that the end-wall heating significantly affects the axial temperature profile of the bed,particularly when the length-to-diameter ratio is low.Particles near the wall heated faster and became more thermally uniform than those in the center of the drum.The region affected by the end heating gradually increased over time.Increasing the rotation speed enhanced the heat conduction rate,and increasing the fill level reduced the mean temperature and thermal uniformity of the granular bed.Heat transfer was also simulated for drums with different length-to-diameter ratios.展开更多
Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in...Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine,according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999.The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope.If the discrete element method adopts the same assumption as the limit equilibrium method,the results of two methods are in good agreement;while if the assumption is not adopted in the discrete element method,the criticalφnumerically calculated is less than the one calculated by use of the limit equilibrium method for the sameC.Thus,from an engineering point of view,the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.展开更多
To address the problem that granular compound fertilizer is prone to agglomeration during mechanized direct seeding of oilseed rape in the middle and lower reaches of the Yangtze River,which causes clogging of the fer...To address the problem that granular compound fertilizer is prone to agglomeration during mechanized direct seeding of oilseed rape in the middle and lower reaches of the Yangtze River,which causes clogging of the fertilizer discharger and leads to a reduction in the uniformity and stability of fertilizer discharge,research on the crushing mechanism of caking compound fertilizer was performed.Considering that it is difficult to measure the bonding force between caking fertilizer particles directly,a simulation model of caking composite fertilizer was established with the bonding model in EDEM discrete element software.To decrease error between the simulation and physical test results,the normal contact stiffness,tangential contact stiffness,critical normal stress,critical tangential stress,bonding radius,and other parameters of the bonding model of caking composite fertilizer were calibrated.The three-dimensional structure of the caking composite fertilizer was obtained via three-dimensional scanning,the critical crushing displacement and critical crushing force of the caking composite fertilizer were measured via compression testing with a mass spectrometer,and the optimal parameter combination of the bonding model was determined via EDEM discrete element simulation of the Plackett-Burman test,steepest ascent test,and Box-Behnken test.The results of the simulated compression tests under the optimal parameter combination show that the relative errors of the critical crushing displacement and critical crushing force with respect to the physical test results were 0.296%and 0.343%,respectively.Using the crushing rate of caking compound fertilizer as an evaluation index,the feasibility of the calibrated parameters was verified for a four-head spiral two-row fertilizer discharger installed in a direct seeding machine for oilseed rape.The results show that the relative errors of the caking fertilizer crushing rates from the simulation relative to those of the bench and field tests were 5.81%and 5.06%,respectively,indicating that the calibration parameters of the discrete element model were accurate and could be used for parameter analysis of caking fertilizer with a discrete element model.These results can provide a reference for the structural optimization of fertilizer discharger crushing of caking fertilizer of direct seeding machine for oilseed rape.展开更多
This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The ...This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The cohesive elements with softening traction-separation relations and damage initiation and evolution laws are embedded between solid elements in regions of interest in the initial mesh to model potential cracks. The initial mesh can consist of tetrahedrons, wedges, bricks or a mixture of these elements. Neither remeshing nor objective crack propagation criteria are needed. Four examples of concrete specimens, including a wedgesplitting test, a notched beam under torsion, a pull-out test of an anchored cylinder and a notched beam under impact, were modelled and analysed. The simulated crack propagation processes and load-displacement curves agreed well with test results or other numerical simulations for all the examples using initial meshes with reasonable densities. Making use of Abaqus＇s rich pre/post- processing functionalities and powerful standard/explicit solvers, the developed method offers a practical tool for engineering analysts to model complex 3D fracture problems.展开更多
In the past decade, numerical modelling has been increasingly used for simulating the mechanical behaviour of naturally fractured rock masses. In this paper, we introduce new algorithms for spatial and temporal analys...In the past decade, numerical modelling has been increasingly used for simulating the mechanical behaviour of naturally fractured rock masses. In this paper, we introduce new algorithms for spatial and temporal analyses of newly generated fractures and blocks using an integrated discrete fracture network (DFN)-finite-discrete element method (FDEM) (DFN-FDEM) modelling approach. A fracture line calculator and analysis technique (i.e. discrete element method (DEM) fracture analysis, DEMFA) calculates the geometrical aspects of induced fractures using a dilation criterion. The resultant two-dimensional (2D) blocks are then identified and characterised using a graph structure. Block tracking trees allow track of newly generated blocks across timesteps and to analyse progressive breakage of these blocks into smaller blocks. Fracture statistics (number and total length of initial and induced fractures) are then related to the block forming processes to investigate damage evolution. The combination of various proposed methodologies together across various stages of modelling processes provides new insights to investigate the dependency of structure's resistance on the initial fracture configuration.展开更多
基金Project(51378006) supported by National Natural Science Foundation of ChinaProject(141076) supported by Huoyingdong Foundation of the Ministry of Education of China+1 种基金Project(2242015R30027) supported by Excellent Young Teacher Program of Southeast University,ChinaProject(BK20140109) supported by the Natural Science Foundation of Jiangsu Province,China
文摘The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.
基金Project supported by the National Natural Science Foundation of China(Nos.59978005 and 10232024)the National Distinguished Youth Fund of China(No.10025212).
文摘A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.
基金The authors wish to thank National Key R&D Program of China(Grant No.2022YFC308100)the National Nature Science Foundation of China(Grant Nos.42107172 and 42072303)for financial support.
文摘Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.
文摘This study presents the first step of a research project that aims at using a three-dimensional （3D） hybridfinite-discrete element method （FDEM） to investigate the development of an excavation damaged zone（EDZ） around tunnels in a clay shale formation known as Opalinus Clay. The 3D FDEM was first calibratedagainst standard laboratory experiments, including Brazilian disc test and uniaxial compression test. Theeffect of increasing confining pressure on the mechanical response and fracture propagation of the rockwas quantified under triaxial compression tests. Polyaxial （or true triaxial） simulations highlighted theeffect of the intermediate principal stress （s2） on fracture directions in the model： as the intermediateprincipal stress increased, fractures tended to align in the direction parallel to the plane defined by themajor and intermediate principal stresses. The peak strength was also shown to vary with changing s2. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
文摘We simulated three-dimensional heat transfer inside a horizontal rotating drum using the discrete element method and a thermal conduction model.The aim was to determine the effect of end-wall heating on thermal behavior of a granular bed.The simulation showed that the end-wall heating significantly affects the axial temperature profile of the bed,particularly when the length-to-diameter ratio is low.Particles near the wall heated faster and became more thermally uniform than those in the center of the drum.The region affected by the end heating gradually increased over time.Increasing the rotation speed enhanced the heat conduction rate,and increasing the fill level reduced the mean temperature and thermal uniformity of the granular bed.Heat transfer was also simulated for drums with different length-to-diameter ratios.
基金supported by the National 973 Project(Grant No.2002CB412703)the Important AspecProject of the Chinese Academy of Sciences(Grant No.KJCX2_SW_L1).
文摘Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine,according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999.The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope.If the discrete element method adopts the same assumption as the limit equilibrium method,the results of two methods are in good agreement;while if the assumption is not adopted in the discrete element method,the criticalφnumerically calculated is less than the one calculated by use of the limit equilibrium method for the sameC.Thus,from an engineering point of view,the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.
基金supported by the China Agricultural Research System of MOF and MARA(NO:CARS-12).
文摘To address the problem that granular compound fertilizer is prone to agglomeration during mechanized direct seeding of oilseed rape in the middle and lower reaches of the Yangtze River,which causes clogging of the fertilizer discharger and leads to a reduction in the uniformity and stability of fertilizer discharge,research on the crushing mechanism of caking compound fertilizer was performed.Considering that it is difficult to measure the bonding force between caking fertilizer particles directly,a simulation model of caking composite fertilizer was established with the bonding model in EDEM discrete element software.To decrease error between the simulation and physical test results,the normal contact stiffness,tangential contact stiffness,critical normal stress,critical tangential stress,bonding radius,and other parameters of the bonding model of caking composite fertilizer were calibrated.The three-dimensional structure of the caking composite fertilizer was obtained via three-dimensional scanning,the critical crushing displacement and critical crushing force of the caking composite fertilizer were measured via compression testing with a mass spectrometer,and the optimal parameter combination of the bonding model was determined via EDEM discrete element simulation of the Plackett-Burman test,steepest ascent test,and Box-Behnken test.The results of the simulated compression tests under the optimal parameter combination show that the relative errors of the critical crushing displacement and critical crushing force with respect to the physical test results were 0.296%and 0.343%,respectively.Using the crushing rate of caking compound fertilizer as an evaluation index,the feasibility of the calibrated parameters was verified for a four-head spiral two-row fertilizer discharger installed in a direct seeding machine for oilseed rape.The results show that the relative errors of the caking fertilizer crushing rates from the simulation relative to those of the bench and field tests were 5.81%and 5.06%,respectively,indicating that the calibration parameters of the discrete element model were accurate and could be used for parameter analysis of caking fertilizer with a discrete element model.These results can provide a reference for the structural optimization of fertilizer discharger crushing of caking fertilizer of direct seeding machine for oilseed rape.
基金supported by EPSRC UK(No.EP/F00656X/1)Xiangting Su's one-year visit to the University of Liverpoosupported by the China Scholarship Council and the National Natural Science Foundation of China(No.50579081).
文摘This study proposes an algorithm of embedding cohesive elements in Abaqus and develops the computer code to model 3D complex cragk propagation in quasi-brittle materials in a relatively easy and efficient manner. The cohesive elements with softening traction-separation relations and damage initiation and evolution laws are embedded between solid elements in regions of interest in the initial mesh to model potential cracks. The initial mesh can consist of tetrahedrons, wedges, bricks or a mixture of these elements. Neither remeshing nor objective crack propagation criteria are needed. Four examples of concrete specimens, including a wedgesplitting test, a notched beam under torsion, a pull-out test of an anchored cylinder and a notched beam under impact, were modelled and analysed. The simulated crack propagation processes and load-displacement curves agreed well with test results or other numerical simulations for all the examples using initial meshes with reasonable densities. Making use of Abaqus＇s rich pre/post- processing functionalities and powerful standard/explicit solvers, the developed method offers a practical tool for engineering analysts to model complex 3D fracture problems.
文摘In the past decade, numerical modelling has been increasingly used for simulating the mechanical behaviour of naturally fractured rock masses. In this paper, we introduce new algorithms for spatial and temporal analyses of newly generated fractures and blocks using an integrated discrete fracture network (DFN)-finite-discrete element method (FDEM) (DFN-FDEM) modelling approach. A fracture line calculator and analysis technique (i.e. discrete element method (DEM) fracture analysis, DEMFA) calculates the geometrical aspects of induced fractures using a dilation criterion. The resultant two-dimensional (2D) blocks are then identified and characterised using a graph structure. Block tracking trees allow track of newly generated blocks across timesteps and to analyse progressive breakage of these blocks into smaller blocks. Fracture statistics (number and total length of initial and induced fractures) are then related to the block forming processes to investigate damage evolution. The combination of various proposed methodologies together across various stages of modelling processes provides new insights to investigate the dependency of structure's resistance on the initial fracture configuration.
