In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. ...In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.展开更多
With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their conta...With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.展开更多
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 user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂F...A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.展开更多
The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated ...The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated by discrete element method(DEM) in three dimensions under variant rotation speeds,filling degree,based on the background of induration process of iron ore pellets.The influences of the mentioned factors on the maximum thickness of the active layer and the average velocity of particles have been investigated.The average velocity of particles increases with Froude number following the power function over a wide range,and the maximum thickness rises with increasing rotation speed in a way of logarithm.The influence of the filling degree f on the maximum thickness exhibits a good linearity under two classic regimes,but the increasing of the average velocity of the active layer is limited at f=0.4.This basic research highlights the impact of the active layer within rotary kilns,and lays a good foundation for the further investigation in mixing and heat transfer within the particle bed inside rotary kilns.展开更多
Methods that can efficiently model the effects of rock joints on rock mass behavior can be beneficial in rock engineering. The suitability of equivalent rock mass(ERM) technique based upon particle methods is investig...Methods that can efficiently model the effects of rock joints on rock mass behavior can be beneficial in rock engineering. The suitability of equivalent rock mass(ERM) technique based upon particle methods is investigated. The ERM methodology is first validated by comparing calculated and experimental data of lab triaxial compression test on a set of cylindrical rock mass samples, each containing a single joint oriented in various dip angles. The simulated results are then used to study the stress-strain nonlinearity and failure mechanism as a function of the joint dip angle and confining stress. The anisotropy and size effects are also investigated by using multi-scale cubic ERM models subjected to triaxial compression test. The deformation and failure behavior are found to be influenced by joint degradation, the micro-crack formation in the intact rock, the interaction between two joints, and the interactions of micro-cracks and joints.展开更多
Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a d...Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a deep circular tunnel in a rock mass with multiple weakness planes using a 2D combined finite element method/discrete element method(FEM/DEM).Conventional triaxial compression tests were performed on typical hard rock(marble)specimens under a range of confinement stress conditions to validate the rationale and accuracy of the proposed numerical approach.Parametric analysis was subsequently conducted to investigate the influence of inclination angle,and length on the crack propagation behavior,failure mode,energy evolution,and displacement distribution of the surrounding rock.The results show that the inclination angle strongly affects tunnel stability,and the failure intensity and damage range increase with increasing inclination angle and then decrease.The dynamic disasters are more likely with increasing weak plane length.Shearing and sliding along multiple weak planes are also consistently accompanied by kinetic energy fluctuations and surges after unloading,which implies a potentially violent dynamic response around a deeply-buried tunnel.Interactions between slabbing and shearing near the excavation boundaries are also discussed.The results presented here provide important insight into deep tunnel failure in hard rock influenced by both unloading disturbance and tectonic activation.展开更多
In order to study the interaction between various fouling particles and ballast,a multi-layer and multi-scale discrete element model(DEM)including the sleeper,ballast bed and the surface layer of subgrade was develope...In order to study the interaction between various fouling particles and ballast,a multi-layer and multi-scale discrete element model(DEM)including the sleeper,ballast bed and the surface layer of subgrade was developed.Two typical fouling particles,the hard particles(sand)and soft ones(coal fines),are considered.A support stiffness test of the ballast bed under various fouling conditions was conducted to calibrate the microscopic parameters of the contact model.With the model,the influence of fouling particles on the mechanical behavior and deformation of the ballast bed was analyzed from macro and micro perspectives.The results show that the increase in the strength of the fouling particles enlarges the stiffness of the ballast bed.Hard particles increase the uniformity coefficient of the contact force bondγof ballast by 50.4%.Fouling particles increase the average stress in the subgrade,soft particles by 2 kPa and hard particles by 1 kPa.Hard particles can reduce the elasticity,plastic deformation and energy dissipation in the track structure.As the fouling particle changes from hard to soft,the proportion of the settlement in ballast bed increases to 40.5%and surface layer of swbgrade settlement decreases to 59.5%.Thus,the influence of fouling particles should be considered carefully in railway design and maintenance.展开更多
For a deeper understanding of the deformation failure behavior of jointed rock, numerical compression simulations are carried out on a rock specimen containing non-persistent joints under confining pressure with the b...For a deeper understanding of the deformation failure behavior of jointed rock, numerical compression simulations are carried out on a rock specimen containing non-persistent joints under confining pressure with the bondedparticle model. The microscopic parameters which can reflect the macroscopic mechanical properties and failure behavior of artificial jointed specimens are firstly calibrated. Then, the influence of joint inclination and confining pressure on stress-strain curves, crack patterns, and contact force distributions of jointed rock are investigated. The simulation results show that both the compressive strength and elastic modulus of the specimens increase with increasing confining pressure, and these two mechanical parameters decrease first and then increase with the increase of joints inclination. The sensitivity of strength and elastic modulus to confining pressure is not the same in different joints inclinations, which has the least impact on specimens with α=90°. Under low confining pressure, the failure modes are controlled by the joint inclination. As the confining pressure increased, the initiation and propagation of tensile crack are gradually inhibited, and the failure mode is transferred from tensile failure to shear-compression failure. Finally, the reinforcement effect of prestressed bolt support on engineering fractured rock mass is discussed.展开更多
This study aims to reveal the macroscopic permanent deformation(PD)behavior and the internal structural evolution of construction and demolition waste(CDW)under loading.Firstly,the initial matric suction of CDW was me...This study aims to reveal the macroscopic permanent deformation(PD)behavior and the internal structural evolution of construction and demolition waste(CDW)under loading.Firstly,the initial matric suction of CDW was measured by the filter paper method.Secondly,the PD of CDW with different humidity and stress states was investigated by repeated load triaxial tests,and a comprehensive prediction model was established.Finally,the discrete element method was performed to analyze the internal structural evolution of CDW during deformation.These results showed that the VAN-GENUCHTEN model could describe the soil-water characteristic curve of CDW well.The PD increases with the increase of the deviator stress and the number of cyclic loading,but the opposite trend was observed when the initial matric suction and confining pressure increased.The proposed model in this study provides a satisfactory prediction of PD.The discrete element method could accurately simulate the macroscopic PD of CDW,and the shear force,interlock force and sliding content increase with the increase of deviator stress during the deformation.The research could provide useful reference for the deformation stability analysis of CDW under cyclic loading.展开更多
The macro mechanical properties of materials with characteristics of large scale and complicated structural composition can be analyzed through its reconstructed meso-structures.In this work,the meso-structures of tal...The macro mechanical properties of materials with characteristics of large scale and complicated structural composition can be analyzed through its reconstructed meso-structures.In this work,the meso-structures of talus deposits that widely exist in the hydro-power engineering in the southwest of China were first reconstructed by small particles according to the in-situ photographs based on the self-adaptive PCNN digital image processing,and then numerical direct shear tests were carried out for studying the mechanical properties of talus deposits.Results indicate that the reconstructed meso-structures of talus deposits are more consistent with the actual situation because the self-adaptive PCNN digital image processing has a higher discrimination in the details of soil-rock segmentation.The existence and random distribution of rock blocks make the initial shear stiffness,the peak strength and the residual strength higher than those of the "pure soil" with particle size less than 1.25 cm apparently,but reduce the displacements required for the talus deposits reaching its peak shear strength.The increase of rock proportion causes a significant improvement in the internal friction angle of talus deposit,which to a certain degree leads to the characteristics of shear stress-displacement curves having a changing trend from the plastic strain softening deformation to the nonlinear strain hardening deformation,while an unconspicuous increase in cohesion.The uncertainty and heterogeneity of rock distributions cause the differences of rock proportion within shear zone,leading to a relatively strong fluctuation in peak strengths during the shear process,while movement features of rock blocks,such as translation,rotation and crossing,expand the scope of shear zone,increase the required shear force,and also directly lead to the misjudgment that the lower shear strength is obtained from the samples with high rock proportion.That,however,just explains the reason why the shear strength gained from a small amount of indoor test data is not consistent with engineering practice.展开更多
A micromechanical investigation on simple shear behavior of dense granular assemblies was carried out by discrete element method.Three series of numerical tests were performed to examine the effects of initial porosit...A micromechanical investigation on simple shear behavior of dense granular assemblies was carried out by discrete element method.Three series of numerical tests were performed to examine the effects of initial porosity,vertical stress and particle shape on simple shear behavior of the samples,respectively.It was found that during simple shear the directions of principal stress and principal strain increment rotate differently with shear strain level.The non-coaxiality between the two directions decreases with strain level and may greatly affect the shear behavior of the assemblies,especially their peak friction angles.The numerical modelling also reveals that the rotation of the principal direction of fabric anisotropy lags behind that of the major principal stress direction during simple shear,which is described as fabric hyteresis effect.The degrees of fabric and interparticle contact force anisotropies increase as particle angularity increases,whereas the orientations of these anisotropies have not been significantly influenced by particle shape.An extended stress–dilatancy relationship based on ROWE-DAVIS framework was proposed to consider the non-coaxiality effect under principal stress rotation.The model was validated by present numerical results as well as some published physical test and numerical modelled data.展开更多
The discrete element method was used to investigate the microscopic characteristics of granular materials under simple shear loading conditions. A series of simple tests on photo-elastic materials were used as a bench...The discrete element method was used to investigate the microscopic characteristics of granular materials under simple shear loading conditions. A series of simple tests on photo-elastic materials were used as a benchmark. With respect to the original experimental observations, average micro-variables such as the shear stress, shear strain and the volumetric dilatancy were extracted to illustrate the performance of the DEM simulation. The change of anisotropic density distributions of contact normals and contact forces was demonstrated during the course of simple shear. On the basis of microscopic characteristics, an analytical approach was further used to explore the macroscopic behaviors involving anisotropic shear strength and anisotropic stress-dilatancy. This results show that under simple shear loading, anisotropic shear strength arises primarily due to the difference between principal directions of the stress and the fabric. In addition, non-coaxiality, referring to the difference between principal directions of the strain rate and the stress, generates less stress-dilatancy. In particular, the anisotropic hardening and anisotropic stress-dilatancy will reduce to the isotropic hardening and the classical Taylor’s stress-dilatancy under proportional loading.展开更多
The layout of the buckets for tunnel boring machine(TBM)directly affects the muck removal efficiency of cutterhead during excavation.In order to improve the muck removal performance for TBM,the optimal design of bucke...The layout of the buckets for tunnel boring machine(TBM)directly affects the muck removal efficiency of cutterhead during excavation.In order to improve the muck removal performance for TBM,the optimal design of bucket layout was investigated.The whole muck transfer process was simulated by discrete-element method(DEM),including the muck falling,colliding,pilling up,shoveling and transferring into the hopper.The muck model was established based on size distribution analysis of muck samples from the water-supply tunnel project in Jilin Province,China.Then,the influence of the bucket number and the interval angle between buckets on muck removal performance was investigated.The results indicated that,as the number of buckets increased from four to eight,the removed muck increased by 29%and the residual volume decreased by 40.5%,and the process became steadier.Different interval angles between buckets were corresponding to different removed muck irregularly,but the residual muck number increased generally with the angles.The optimal layout of buckets for the cutterhead in this tunnel project was obtained based on the simulation results,and the muck removal performance of the TBM was verified by the actual data in the engineering construction.展开更多
In order to exactly provide scientific basis for pressure dynamic balance control of working chamber of earth pressure balance shield (EPBS),study on optimal arrangement of pressure measurement points in working chamb...In order to exactly provide scientific basis for pressure dynamic balance control of working chamber of earth pressure balance shield (EPBS),study on optimal arrangement of pressure measurement points in working chamber was conducted. Based on mathematical description of optimal arrangement for pressure measurement points,fuzzy clustering analysis and discriminant analysis were used to divide pressure regions of nodes on bulkhead. Finally,the selection method of optimal measurement points was proposed,and by selecting d6.28 m EPBS as study object,the case study was conducted. By contrast,based on optimal arrangement scheme of pressure measurement points,through adopting weighted algorithm,the absolute error mean of equivalent pressure of working chamber is the smallest. In addition,pressure curve of optimal arrangement points presents parabola,and it can show the state of pressure distribution on bulkhead truly. It is concluded that the optimal arrangement method of pressure measurement points in working chamber is effective and feasible,and the method can provide basis for realizing high precision pressure control of EPBS.展开更多
Mechanical behaviors of granular materials are complicated and greatly influenced by the particle shape.Current,some composite approaches have been proposed for realistic particle shape modelling within discrete eleme...Mechanical behaviors of granular materials are complicated and greatly influenced by the particle shape.Current,some composite approaches have been proposed for realistic particle shape modelling within discrete element method(DEM),while they cannot give a good representation to the shape and mass properties of a real particle.In this work,a novel algorithm is developed to model an arbitrary particle using a cluster of non-overlapping disks.The algorithm mainly consists of two components:boundary filling and domain filling.In the boundary filling,some disks are placed along the boundary for a precise representation of the particle shape,and some more disks are placed in the domain to give an approximation to the mass properties of the particle in the domain filling.Besides,a simple method is proposed to correct the mass properties of a cluster after domain filling and reduce the number of the disks in a cluster for lower computational load.Moreover,it is another great merit of the algorithm that a cluster generated by the algorithm can be used to simulate the particle breakage because of no overlaps between the disks in a cluster.Finally,several examples are used to show the robust performance of the algorithm.A current FORTRAN version of the algorithm is available by contacting the author.展开更多
Soil-rock mixture(S-RM)is a widely distributed geotechnical medium composed of "soil" and "rock block" different both in size and strength. Internal rock blocks form special and variable meso-struc...Soil-rock mixture(S-RM)is a widely distributed geotechnical medium composed of "soil" and "rock block" different both in size and strength. Internal rock blocks form special and variable meso-structural characteristics of S-RM. The objective of this work was to study the control mechanism of meso-structural characteristics on mechanical properties of S-RM. For S-RM containing randomly generated polygonal rock blocks, a series of biaxial tests based on DEM were conducted. On the basis of research on the effects of rock blocks' breakability and sample lateral boundary type(rigid, flexible) on macroscopic mechanical behavior of S-RM, an expanded Mohr-Coulomb criterion in power function form was proposed to represent the strength envelop. At the mesoscopic level, the variations of meso-structure such as rotation of rock block, and the formation mechanism and evolution process of the shear band during tests were investigated. The results show that for S-RM with a high content of rock block, translation, rotating and breakage of rock blocks have crucial effects on mechanical behavior of S-RM. The formation and location of the shear band inside S-RM sample are also controlled by breakability and arrangement of rock blocks.展开更多
Adding a moving baffle to the drum is a new way to enhance the motion and mixing of particles in rotating drums.To obtain its influence on binary particles,horizontal rotating drums provided with a moving baffle were ...Adding a moving baffle to the drum is a new way to enhance the motion and mixing of particles in rotating drums.To obtain its influence on binary particles,horizontal rotating drums provided with a moving baffle were investigated by discrete element method(DEM).AtΩ=15 r/min,increasing the length of moving baffle can increase the fluctuation amplitude of average particle velocity.AtΩ=60 r/min,the influence of the moving baffle on the average velocity fluctuation tends to be more random.At both rotational speeds,the moving baffle causes the average particle velocity to fluctuate more sharply.The moving baffle can enhance particle mixing.AtΩ=15 r/min,the moving baffle with length ofδ=1/3 can best enhance particle mixing.However,atΩ=60 r/min,only the moving baffle with a specific length(δ=1/4)can enhance mixing.This basic research has a positive reference value for the application of the moving baffle in industry.展开更多
Irregular columnar jointed structure is a primary irregular columnar morphological tensile fracture. In order to study the geometric features of irregular columnar joints and the new problems in geotechnical engineeri...Irregular columnar jointed structure is a primary irregular columnar morphological tensile fracture. In order to study the geometric features of irregular columnar joints and the new problems in geotechnical engineering, hydraulic and hydropower engineering caused by columnar jointed basaltic mass, Voronoi graph from geometry was introduced to simulate the irregular columnar jointed basaltic mass at Baihetan hydraulic station. Discrete element software UDEC was used to simulate the whole process of rigid bearing plate test. Anisotropic constitutive of columnar joints was adopted to analyze the stress diffusion of rock mass at dam base of Baihetan. The results show that, the compaction property and hysteresis effect are well simulated based on discrete element simulation of Voronoi joint structure by UDEC. Four stages of cyclic loading and unloading process are imaged clearly. The results from in situ rigid bearing plate tests are explicated and the stress diffusion rule of anisotropic body is affected by structure surface. The elements in the stress state of 4-5 MPa are the most, about more than 35% of the total. Appropriate constitutive must be proposed to columnar jointed rock mass with different styles. It has important significance to realize the nonlinear mechanical behavior of irregular columnar jointed basaltic mass.展开更多
The landslide disaster caused by the argillaceous interlayer not only destroys buildings,cultivated land,and roads but also seriously endangers human life and safety.This study concerns the mineral composition of sele...The landslide disaster caused by the argillaceous interlayer not only destroys buildings,cultivated land,and roads but also seriously endangers human life and safety.This study concerns the mineral composition of selected argillaceous interlayer and their strength characteristics.To study the mineral composition of argillaceous interlayers,8 kinds of samples in the southern Jiangsu region of China were analyzed utilizing X-ray diffraction(XRD).The repeated direct shear strength tests(RDST)were carried out on the undisturbed specimens of the argillaceous interlayer.The results show that the argillaceous interlayer with high content of kaolinite shows ductile failure mode,which means that there is no obvious residual strength in the shear process.The arrangement of mineral particles on the shear surface of the specimens after different shear displacements was observed under the scanning electron microscope(SEM).It was observed that mineral particles on the shear surface showed a more directional arrangement with the increase of shear displacement.Furthermore,the influence of shear direction on the argillaceous interlayer with completely oriented mineral particles was studied through numerical experiments with four shear strength mechanisms proposition proposed.The influence of the mineral arrangement on the action occasion and magnitude of dilatancy component of shear strength is clarified in the shear mechanism.展开更多
文摘In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.
基金Project(U1234211)supported by the National Natural Science Foundation of ChinaProject(2013G009-B)supported by China Railway Corporation
文摘With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.
基金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(51208178)supported by the National Natural Science Foundation of ChinaProject(2012M520991)supported by China Postdoctoral Science Foundation
文摘A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.
基金Project(FRF-AS-10-0058) supported by the Fundamental Research Funds for the Central Universities,China
文摘The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated by discrete element method(DEM) in three dimensions under variant rotation speeds,filling degree,based on the background of induration process of iron ore pellets.The influences of the mentioned factors on the maximum thickness of the active layer and the average velocity of particles have been investigated.The average velocity of particles increases with Froude number following the power function over a wide range,and the maximum thickness rises with increasing rotation speed in a way of logarithm.The influence of the filling degree f on the maximum thickness exhibits a good linearity under two classic regimes,but the increasing of the average velocity of the active layer is limited at f=0.4.This basic research highlights the impact of the active layer within rotary kilns,and lays a good foundation for the further investigation in mixing and heat transfer within the particle bed inside rotary kilns.
基金Projects(51074014,51174014) supported by the National Natural Science Foundation of China
文摘Methods that can efficiently model the effects of rock joints on rock mass behavior can be beneficial in rock engineering. The suitability of equivalent rock mass(ERM) technique based upon particle methods is investigated. The ERM methodology is first validated by comparing calculated and experimental data of lab triaxial compression test on a set of cylindrical rock mass samples, each containing a single joint oriented in various dip angles. The simulated results are then used to study the stress-strain nonlinearity and failure mechanism as a function of the joint dip angle and confining stress. The anisotropy and size effects are also investigated by using multi-scale cubic ERM models subjected to triaxial compression test. The deformation and failure behavior are found to be influenced by joint degradation, the micro-crack formation in the intact rock, the interaction between two joints, and the interactions of micro-cracks and joints.
基金Projects(52004143,51774194)supported by the National Natural Science Foundation of ChinaProject(2020M670781)supported by the China Postdoctoral Science Foundation+2 种基金Project(SKLGDUEK2021)supported by the State Key Laboratory for GeoMechanics and Deep Underground Engineering,ChinaProject(U1806208)supported by the NSFC-Shandong Joint Fund,ChinaProject(2018GSF117023)supported by the Key Research and Development Program of Shandong Province,China。
文摘Natural geological structures in rock(e.g.,joints,weakness planes,defects)play a vital role in the stability of tunnels and underground operations during construction.We investigated the failure characteristics of a deep circular tunnel in a rock mass with multiple weakness planes using a 2D combined finite element method/discrete element method(FEM/DEM).Conventional triaxial compression tests were performed on typical hard rock(marble)specimens under a range of confinement stress conditions to validate the rationale and accuracy of the proposed numerical approach.Parametric analysis was subsequently conducted to investigate the influence of inclination angle,and length on the crack propagation behavior,failure mode,energy evolution,and displacement distribution of the surrounding rock.The results show that the inclination angle strongly affects tunnel stability,and the failure intensity and damage range increase with increasing inclination angle and then decrease.The dynamic disasters are more likely with increasing weak plane length.Shearing and sliding along multiple weak planes are also consistently accompanied by kinetic energy fluctuations and surges after unloading,which implies a potentially violent dynamic response around a deeply-buried tunnel.Interactions between slabbing and shearing near the excavation boundaries are also discussed.The results presented here provide important insight into deep tunnel failure in hard rock influenced by both unloading disturbance and tectonic activation.
基金Project(51978045) supported by the National Natural Science Foundation of ChinaProject([2017]7) supported by Shenshuo Science and Technology Development Project,China。
文摘In order to study the interaction between various fouling particles and ballast,a multi-layer and multi-scale discrete element model(DEM)including the sleeper,ballast bed and the surface layer of subgrade was developed.Two typical fouling particles,the hard particles(sand)and soft ones(coal fines),are considered.A support stiffness test of the ballast bed under various fouling conditions was conducted to calibrate the microscopic parameters of the contact model.With the model,the influence of fouling particles on the mechanical behavior and deformation of the ballast bed was analyzed from macro and micro perspectives.The results show that the increase in the strength of the fouling particles enlarges the stiffness of the ballast bed.Hard particles increase the uniformity coefficient of the contact force bondγof ballast by 50.4%.Fouling particles increase the average stress in the subgrade,soft particles by 2 kPa and hard particles by 1 kPa.Hard particles can reduce the elasticity,plastic deformation and energy dissipation in the track structure.As the fouling particle changes from hard to soft,the proportion of the settlement in ballast bed increases to 40.5%and surface layer of swbgrade settlement decreases to 59.5%.Thus,the influence of fouling particles should be considered carefully in railway design and maintenance.
基金Projects(52004145,51904164)supported by the National Natural Science Foundation of ChinaProject(ZR2020QE119)supported by the Natural Science Foundation of Shandong Province,ChinaProject(SICGM202107)supported by the Open Fund of the Key Laboratory of Mining Disaster Prevention and Control,China。
文摘For a deeper understanding of the deformation failure behavior of jointed rock, numerical compression simulations are carried out on a rock specimen containing non-persistent joints under confining pressure with the bondedparticle model. The microscopic parameters which can reflect the macroscopic mechanical properties and failure behavior of artificial jointed specimens are firstly calibrated. Then, the influence of joint inclination and confining pressure on stress-strain curves, crack patterns, and contact force distributions of jointed rock are investigated. The simulation results show that both the compressive strength and elastic modulus of the specimens increase with increasing confining pressure, and these two mechanical parameters decrease first and then increase with the increase of joints inclination. The sensitivity of strength and elastic modulus to confining pressure is not the same in different joints inclinations, which has the least impact on specimens with α=90°. Under low confining pressure, the failure modes are controlled by the joint inclination. As the confining pressure increased, the initiation and propagation of tensile crack are gradually inhibited, and the failure mode is transferred from tensile failure to shear-compression failure. Finally, the reinforcement effect of prestressed bolt support on engineering fractured rock mass is discussed.
基金Project(52025085)supported by the National Science Fund for Distinguished Young Scholars,ChinaProjects(51927814,51878078)supported by the National Natural Science Foundation of China+3 种基金Project(2018-025)supported by the Training Program for High-level Technical Personnel in Transportation Industry,ChinaProject(CTKY-PTRC 2018-003)supported by the Design Theory,Method and Demonstration of Durability Asphalt Pavement Based on Heavy-duty Traffic Conditions in Shanghai Area,ChinaProject(2020RC4048)supported by the Science and Technology Innovation Program of Hunan Province,ChinaProject(SJCX202001)supported by the Construction Project for Graduate Students of Changsha University of Science&Technology,China。
文摘This study aims to reveal the macroscopic permanent deformation(PD)behavior and the internal structural evolution of construction and demolition waste(CDW)under loading.Firstly,the initial matric suction of CDW was measured by the filter paper method.Secondly,the PD of CDW with different humidity and stress states was investigated by repeated load triaxial tests,and a comprehensive prediction model was established.Finally,the discrete element method was performed to analyze the internal structural evolution of CDW during deformation.These results showed that the VAN-GENUCHTEN model could describe the soil-water characteristic curve of CDW well.The PD increases with the increase of the deviator stress and the number of cyclic loading,but the opposite trend was observed when the initial matric suction and confining pressure increased.The proposed model in this study provides a satisfactory prediction of PD.The discrete element method could accurately simulate the macroscopic PD of CDW,and the shear force,interlock force and sliding content increase with the increase of deviator stress during the deformation.The research could provide useful reference for the deformation stability analysis of CDW under cyclic loading.
基金Project(2013BAB06B00) supported by the National Key Technology R&D Programof ChinaProject(2011CB013504) supported by the National Basic Research Program of ChinaProject(50911130366) supported by the National Natural Science Foundation of China
文摘The macro mechanical properties of materials with characteristics of large scale and complicated structural composition can be analyzed through its reconstructed meso-structures.In this work,the meso-structures of talus deposits that widely exist in the hydro-power engineering in the southwest of China were first reconstructed by small particles according to the in-situ photographs based on the self-adaptive PCNN digital image processing,and then numerical direct shear tests were carried out for studying the mechanical properties of talus deposits.Results indicate that the reconstructed meso-structures of talus deposits are more consistent with the actual situation because the self-adaptive PCNN digital image processing has a higher discrimination in the details of soil-rock segmentation.The existence and random distribution of rock blocks make the initial shear stiffness,the peak strength and the residual strength higher than those of the "pure soil" with particle size less than 1.25 cm apparently,but reduce the displacements required for the talus deposits reaching its peak shear strength.The increase of rock proportion causes a significant improvement in the internal friction angle of talus deposit,which to a certain degree leads to the characteristics of shear stress-displacement curves having a changing trend from the plastic strain softening deformation to the nonlinear strain hardening deformation,while an unconspicuous increase in cohesion.The uncertainty and heterogeneity of rock distributions cause the differences of rock proportion within shear zone,leading to a relatively strong fluctuation in peak strengths during the shear process,while movement features of rock blocks,such as translation,rotation and crossing,expand the scope of shear zone,increase the required shear force,and also directly lead to the misjudgment that the lower shear strength is obtained from the samples with high rock proportion.That,however,just explains the reason why the shear strength gained from a small amount of indoor test data is not consistent with engineering practice.
基金Projects(50909057,51208294,41372319)supported by the National Natural Science Foundation of ChinaProject(15ZZ081)supported by Innovation Program of Shanghai Municipal Education Commission,ChinaProject(20131129)supported by Innovation Program of Shanghai Postgraduate Education,China
文摘A micromechanical investigation on simple shear behavior of dense granular assemblies was carried out by discrete element method.Three series of numerical tests were performed to examine the effects of initial porosity,vertical stress and particle shape on simple shear behavior of the samples,respectively.It was found that during simple shear the directions of principal stress and principal strain increment rotate differently with shear strain level.The non-coaxiality between the two directions decreases with strain level and may greatly affect the shear behavior of the assemblies,especially their peak friction angles.The numerical modelling also reveals that the rotation of the principal direction of fabric anisotropy lags behind that of the major principal stress direction during simple shear,which is described as fabric hyteresis effect.The degrees of fabric and interparticle contact force anisotropies increase as particle angularity increases,whereas the orientations of these anisotropies have not been significantly influenced by particle shape.An extended stress–dilatancy relationship based on ROWE-DAVIS framework was proposed to consider the non-coaxiality effect under principal stress rotation.The model was validated by present numerical results as well as some published physical test and numerical modelled data.
基金Projects(10972159, 41272291, 51238009) supported by the National Natural Science Foundation of China Project supported by the Fundamental Research Funds of the central Universities
文摘The discrete element method was used to investigate the microscopic characteristics of granular materials under simple shear loading conditions. A series of simple tests on photo-elastic materials were used as a benchmark. With respect to the original experimental observations, average micro-variables such as the shear stress, shear strain and the volumetric dilatancy were extracted to illustrate the performance of the DEM simulation. The change of anisotropic density distributions of contact normals and contact forces was demonstrated during the course of simple shear. On the basis of microscopic characteristics, an analytical approach was further used to explore the macroscopic behaviors involving anisotropic shear strength and anisotropic stress-dilatancy. This results show that under simple shear loading, anisotropic shear strength arises primarily due to the difference between principal directions of the stress and the fabric. In addition, non-coaxiality, referring to the difference between principal directions of the strain rate and the stress, generates less stress-dilatancy. In particular, the anisotropic hardening and anisotropic stress-dilatancy will reduce to the isotropic hardening and the classical Taylor’s stress-dilatancy under proportional loading.
基金Project(51475478)supported by the National Natural Science Foundation of ChinaProject(2012AA041801)supported by the National High Technology Research and Development Program of China+1 种基金Project(2014FJ1002)supported by the Science and Technology Major Project of Hunan Province,ChinaProject(2013CB035401)supported by the National Basic Research Program of China。
文摘The layout of the buckets for tunnel boring machine(TBM)directly affects the muck removal efficiency of cutterhead during excavation.In order to improve the muck removal performance for TBM,the optimal design of bucket layout was investigated.The whole muck transfer process was simulated by discrete-element method(DEM),including the muck falling,colliding,pilling up,shoveling and transferring into the hopper.The muck model was established based on size distribution analysis of muck samples from the water-supply tunnel project in Jilin Province,China.Then,the influence of the bucket number and the interval angle between buckets on muck removal performance was investigated.The results indicated that,as the number of buckets increased from four to eight,the removed muck increased by 29%and the residual volume decreased by 40.5%,and the process became steadier.Different interval angles between buckets were corresponding to different removed muck irregularly,but the residual muck number increased generally with the angles.The optimal layout of buckets for the cutterhead in this tunnel project was obtained based on the simulation results,and the muck removal performance of the TBM was verified by the actual data in the engineering construction.
基金Project(2007CB714006) supported by the National Basic Research Program of China
文摘In order to exactly provide scientific basis for pressure dynamic balance control of working chamber of earth pressure balance shield (EPBS),study on optimal arrangement of pressure measurement points in working chamber was conducted. Based on mathematical description of optimal arrangement for pressure measurement points,fuzzy clustering analysis and discriminant analysis were used to divide pressure regions of nodes on bulkhead. Finally,the selection method of optimal measurement points was proposed,and by selecting d6.28 m EPBS as study object,the case study was conducted. By contrast,based on optimal arrangement scheme of pressure measurement points,through adopting weighted algorithm,the absolute error mean of equivalent pressure of working chamber is the smallest. In addition,pressure curve of optimal arrangement points presents parabola,and it can show the state of pressure distribution on bulkhead truly. It is concluded that the optimal arrangement method of pressure measurement points in working chamber is effective and feasible,and the method can provide basis for realizing high precision pressure control of EPBS.
基金Project(2011CB013504)supported by the National Basic Research Program(973 Program)of ChinaProject(2013BAB06B01)supported by Key Projects in the National Science&Technology Pillar Program during the Twelfth Five-year Plan Period,China+1 种基金Projects(51309089,51479049)supported by National Natural Science Foundation of ChinaProject(487237)supported by Natural Sciences and Engineering Research Council of Canada
文摘Mechanical behaviors of granular materials are complicated and greatly influenced by the particle shape.Current,some composite approaches have been proposed for realistic particle shape modelling within discrete element method(DEM),while they cannot give a good representation to the shape and mass properties of a real particle.In this work,a novel algorithm is developed to model an arbitrary particle using a cluster of non-overlapping disks.The algorithm mainly consists of two components:boundary filling and domain filling.In the boundary filling,some disks are placed along the boundary for a precise representation of the particle shape,and some more disks are placed in the domain to give an approximation to the mass properties of the particle in the domain filling.Besides,a simple method is proposed to correct the mass properties of a cluster after domain filling and reduce the number of the disks in a cluster for lower computational load.Moreover,it is another great merit of the algorithm that a cluster generated by the algorithm can be used to simulate the particle breakage because of no overlaps between the disks in a cluster.Finally,several examples are used to show the robust performance of the algorithm.A current FORTRAN version of the algorithm is available by contacting the author.
基金Projects(51323014,51479095)supported by the National Natural Science Foundation of ChinaProject(20111081125)supported by Independent Research Plan of Tsinghua University,ChinaProject(2013-KY-4)supported by the State Key Laboratory of Hydroscience and Engineering Project,China
文摘Soil-rock mixture(S-RM)is a widely distributed geotechnical medium composed of "soil" and "rock block" different both in size and strength. Internal rock blocks form special and variable meso-structural characteristics of S-RM. The objective of this work was to study the control mechanism of meso-structural characteristics on mechanical properties of S-RM. For S-RM containing randomly generated polygonal rock blocks, a series of biaxial tests based on DEM were conducted. On the basis of research on the effects of rock blocks' breakability and sample lateral boundary type(rigid, flexible) on macroscopic mechanical behavior of S-RM, an expanded Mohr-Coulomb criterion in power function form was proposed to represent the strength envelop. At the mesoscopic level, the variations of meso-structure such as rotation of rock block, and the formation mechanism and evolution process of the shear band during tests were investigated. The results show that for S-RM with a high content of rock block, translation, rotating and breakage of rock blocks have crucial effects on mechanical behavior of S-RM. The formation and location of the shear band inside S-RM sample are also controlled by breakability and arrangement of rock blocks.
基金Project(51676032)supported by the National Natural Science Foundation of ChinaProject(IRT_17R19)supported by the Program for Changjiang Scholars and Innovative Research Team in University,China
文摘Adding a moving baffle to the drum is a new way to enhance the motion and mixing of particles in rotating drums.To obtain its influence on binary particles,horizontal rotating drums provided with a moving baffle were investigated by discrete element method(DEM).AtΩ=15 r/min,increasing the length of moving baffle can increase the fluctuation amplitude of average particle velocity.AtΩ=60 r/min,the influence of the moving baffle on the average velocity fluctuation tends to be more random.At both rotational speeds,the moving baffle causes the average particle velocity to fluctuate more sharply.The moving baffle can enhance particle mixing.AtΩ=15 r/min,the moving baffle with length ofδ=1/3 can best enhance particle mixing.However,atΩ=60 r/min,only the moving baffle with a specific length(δ=1/4)can enhance mixing.This basic research has a positive reference value for the application of the moving baffle in industry.
基金Projects(50979030, 50911130366) supported by the National Natural Science Foundation of ChinaProject(2009B14014) supported by the Fundamental Research Funds for the Central Universities in China
文摘Irregular columnar jointed structure is a primary irregular columnar morphological tensile fracture. In order to study the geometric features of irregular columnar joints and the new problems in geotechnical engineering, hydraulic and hydropower engineering caused by columnar jointed basaltic mass, Voronoi graph from geometry was introduced to simulate the irregular columnar jointed basaltic mass at Baihetan hydraulic station. Discrete element software UDEC was used to simulate the whole process of rigid bearing plate test. Anisotropic constitutive of columnar joints was adopted to analyze the stress diffusion of rock mass at dam base of Baihetan. The results show that, the compaction property and hysteresis effect are well simulated based on discrete element simulation of Voronoi joint structure by UDEC. Four stages of cyclic loading and unloading process are imaged clearly. The results from in situ rigid bearing plate tests are explicated and the stress diffusion rule of anisotropic body is affected by structure surface. The elements in the stress state of 4-5 MPa are the most, about more than 35% of the total. Appropriate constitutive must be proposed to columnar jointed rock mass with different styles. It has important significance to realize the nonlinear mechanical behavior of irregular columnar jointed basaltic mass.
基金Project(41672258) supported by the National Natural Science Foundation of ChinaProject(2018045) supported by the Land and Resources Science&Technology Project of Jiangsu Province,China。
文摘The landslide disaster caused by the argillaceous interlayer not only destroys buildings,cultivated land,and roads but also seriously endangers human life and safety.This study concerns the mineral composition of selected argillaceous interlayer and their strength characteristics.To study the mineral composition of argillaceous interlayers,8 kinds of samples in the southern Jiangsu region of China were analyzed utilizing X-ray diffraction(XRD).The repeated direct shear strength tests(RDST)were carried out on the undisturbed specimens of the argillaceous interlayer.The results show that the argillaceous interlayer with high content of kaolinite shows ductile failure mode,which means that there is no obvious residual strength in the shear process.The arrangement of mineral particles on the shear surface of the specimens after different shear displacements was observed under the scanning electron microscope(SEM).It was observed that mineral particles on the shear surface showed a more directional arrangement with the increase of shear displacement.Furthermore,the influence of shear direction on the argillaceous interlayer with completely oriented mineral particles was studied through numerical experiments with four shear strength mechanisms proposition proposed.The influence of the mineral arrangement on the action occasion and magnitude of dilatancy component of shear strength is clarified in the shear mechanism.
