Rubble mound breakwaters, a prevalent type of sloping breakwater structure, are extensively employed in port and coastal infrastructure projects. Under soft soil foundation conditions, the process of squeezing silt by...Rubble mound breakwaters, a prevalent type of sloping breakwater structure, are extensively employed in port and coastal infrastructure projects. Under soft soil foundation conditions, the process of squeezing silt by riprap is implemented to enhance bearing capacity through soft soil replacement and compaction. However, predicting the depth law of squeezing silt by riprap and understanding its mechanism remain significant engineering design challenges.This study employs particle flow code(PFC) based on the discrete element method to simulate the squeezing silt process by riprap, examining variations in depth law under different geological conditions and its mechanical characteristics.Through calibration of the PFC model's meso-parameters via macro-experiments, the study analyzes the effects of riprap size, drop height, and soft soil properties on the depth law of squeezing silt. Findings demonstrate that riprap drop height and soft soil thickness substantially influence the depth, while appropriate calibration of meso-parameters enhances simulation accuracy. This research contributes theoretical and practical guidance for optimizing rubble mound breakwater design, understanding squeezing silt mechanisms, construction practices, and riprap quantity estimation.展开更多
Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass...Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass is very important to elucidate the permeability evolution mechanisms, which is difficult to realize in laboratory, but easy to be achieved in numerical simulations. Therefore, the particle flow code （PFC）, a discrete element method, is used to simulate permeability behaviors of rock materials in this study. Owe to the limitation of the existed solid-fluid coupling algorithm in PFC, an improved flow-coupling algorithm is presented to better reflect the preferential flow in rock fractures. The comparative analysis is conducted between original and improved algorithm when simulating rock permeability evolution during triaxial compression, showing that the improved algorithm can better describe the experimental phenomenon. Furthermore, the evolution of pore pressure and flow rate distribution during the flow process are analyzed by using the improved algorithm. It is concluded that during the steady flow process in the fractured specimen, the pore pressure and flow rate both prefer transmitting through the fractures rather than rock matrix. Based on the results, fractures are divided into the following three types： I） fractures link to both the inlet and outlet, II） fractures only link to the inlet, and III） fractures only link to the outlet. The type I fracture is always the preferential propagating path for both the pore pressure and flow rate. For type II fractures, the pore pressure increases and then becomes steady. However, the flow rate increases first and begins to decrease after the flow reaches the stop end of the fracture and finally vanishes. There is no obvious pore pressure or flow rate concentration within type III fractures.展开更多
This study presents a calibration process of three-dimensional particle flow code(PFC3D)simulation of intact and fissured granite samples.First,laboratory stressestrain response from triaxial testing of intact and fis...This study presents a calibration process of three-dimensional particle flow code(PFC3D)simulation of intact and fissured granite samples.First,laboratory stressestrain response from triaxial testing of intact and fissured granite samples is recalled.Then,PFC3D is introduced,with focus on the bonded particle models(BPM).After that,we present previous studies where intact rock is simulated by means of flatjoint approaches,and how improved accuracy was gained with the help of parametric studies.Then,models of the pre-fissured rock specimens were generated,including modeled fissures in the form of“smooth joint”type contacts.Finally,triaxial testing simulations of 1 t 2 and 2 t 3 jointed rock specimens were performed.Results show that both elastic behavior and the peak strength levels are closely matched,without any additional fine tuning of micro-mechanical parameters.Concerning the postfailure behavior,models reproduce the trends of decreasing dilation with increasing confinement and plasticity.However,the dilation values simulated are larger than those observed in practice.This is attributed to the difficulty in modeling some phenomena of fissured rock behaviors,such as rock piece corner crushing with dust production and interactions between newly formed shear bands or axial splitting cracks with pre-existing joints.展开更多
Based on the principle of 3D particle flow code,a numerical landslide run-out model is presented to simulate the failure process of the Zhenggang landslide(in southwestern China) under the effect of water after a rain...Based on the principle of 3D particle flow code,a numerical landslide run-out model is presented to simulate the failure process of the Zhenggang landslide(in southwestern China) under the effect of water after a rainfall.The relationship between the micro-mechanical parameters and the macro-shear strength of the grain material is determined through numerical calibrations.Then the rainfall effect is considered in numerical simulations and rain-induced sliding processes are performed,which help us to discuss the mechanism of deformation and failure of this landslide together with field observations.It shows the Zhenggang landslide would most likely be activated in Zone I and would gain momentum in Zone II.In order to prevent the potential disaster,a tailing dam is advised to be designed about 175 m downstream from the current landslide boundary of Zone II.Verified by field observations,the presented landslide model can reflect the failure mechanism after rainfall.It can also provide a method to predict the potential disaster and draft disaster prevention measures.展开更多
Polyurethane polymer grouting materials were studied with conventional triaxial tests via the particle flow code in two dimensions（PFC^（2D）） method, and the simulation results agreed with the experimental data. Th...Polyurethane polymer grouting materials were studied with conventional triaxial tests via the particle flow code in two dimensions（PFC^（2D）） method, and the simulation results agreed with the experimental data. The particle flow code method can simulate the mechanical properties of the polymer. The triaxial cyclic loading tests of the polymer material under different confining pressures were carried out via PFC^（2D） to analyze its mechanical performance. The PFC^（2D） simulation results show that the value of the elastic modulus of the polymer decreases slowly at first and fluctuated within a narrow range near the value of the peak strength; the cumulative plastic strain increases slowly at first and then increases rapidly; the peak strength and elastic modulus of polymer increase with the confining pressure; the PFC^（2D） method can be used to quantitatively evaluate the damage behavior of the polymer material and estimate the fatigue life of the materials under fatigue load based on the number and the location of micro-cracks. Thus, the PFC^（2D） method is an effective tool to study polymers.展开更多
Several special mechanical properties,such as dilatancy and compressibility,of cemented paste backfill(CPB)are controlled by its internal microstructure and evolution.The mesoscopic structure changes of CPB during the...Several special mechanical properties,such as dilatancy and compressibility,of cemented paste backfill(CPB)are controlled by its internal microstructure and evolution.The mesoscopic structure changes of CPB during the development process were investigated.On the basis of the scanning electron microscopy(SEM)and mechanical test results of CPB,the particle size information of CPB was extracted,and a two-dimensional particle flow code(PFC)model of CPB was established to analyze the evolution rule of mesoscopic parameters during CPB development.The embedded FISH language in PFC was used to develop a program for establishing a PFC model on the basis of the SEM results.The mesoscopic parameters of CPB samples at different curing times,such as coordination number(C_(n)),contact force chain,and rose diagram,were obtained by recording and loading and used to analyze the intrinsic relationship between mesoscopic parameter variations and macroscopic mechanical response during CPB development.It is of considerable significance to establish the physical model of CPB using the PFC to reveal the mesoscopic structure of CPB.展开更多
A two-dimensionM discrete element code, particle flow code （PFC2D）, is employed to investigate foundations reinforced with horizontal-vertical （H-V） inclusions. The initial states and loading processes of both unr...A two-dimensionM discrete element code, particle flow code （PFC2D）, is employed to investigate foundations reinforced with horizontal-vertical （H-V） inclusions. The initial states and loading processes of both unreinforced and H-V reinforced foundations are simulated by PFC2D method. The interface between particles and reinforcements, and the reinforcement mechanism of the H-V reinforced foundations are studied through stress distribution graphs, displacement vector graphs and contact force graphs. The simulation results demonstrate that the vertical elements of the H-V reinforcement keep the particles from being displaced under the applied load. The H-V reinforcement can distribute the load uniformly over a wider area, thereby improving the bearing capacity of soil foundation.展开更多
The underground or open-pit methods are used for the extraction of mineral resources,each of which is divided into different categories.Coal is one of the mineral resources,which is exploited either by the surface or ...The underground or open-pit methods are used for the extraction of mineral resources,each of which is divided into different categories.Coal is one of the mineral resources,which is exploited either by the surface or the underground methods.The long-wall mining is one of the methods for the underground coal mining.In this method,which is a mechanized one,some machines such as the shearer or plow are used for the mining.The coal mine in Parvadeh,Tabas is a mechanized mine that is extracted by the long-wall mining.The modeling with particle flow code software was used in this mine for the evaluation of plow performance using the coal specifications.In this regard,the sample was first calibrated by sampling from the Parvadeh coal mine and performing the uniaxial and Brazilian tests on the model.Then,the modeling was done by constructing the model and using the variables such as the clearance angle and the linear velocity of the plow.After making 28 models for the plow,the best model of the plow was selected based on the maximum force applied to the machine in the X direction.Finally,the results of this study showed that the best plow performance is for a model with the clearance angle of zero and the linear velocity of 9 mm/min,and the maximum force applied to this model is equal to 39,000 kN in the X direction.展开更多
This study investigated the effects of weathering depth and thickness on the failure mechanisms of rock samples through experimental and numerical methods.The first configuration involved conducting artificial weather...This study investigated the effects of weathering depth and thickness on the failure mechanisms of rock samples through experimental and numerical methods.The first configuration involved conducting artificial weathering on limestone using the freezing and thawing(F-T)for 40 cycles.The mechanical parameters of the samples were measured at the end of the 40th cycle.In the second configuration,a series of specimens underwent salt crystallization(S-C)tests for 20 cycles.Experimental results were validated using discrete element method(DEM).Next,the weathered limestone model with dimensions of 108 mm54 mm were prepared.The weathering layers were tested at four different thicknesses(i.e.2.5 mm,5 mm,7.5 mm,and 10 mm)and three different positions(at the surface,5 mm under the rock surface,and 10 mm under the rock surface).According to the results,weathering depth and thickness have a considerable effect on the failure process.The results also showed a correlation between the values of compressive strength and failure mechanisms associated with the weathering layer.The numerical results revealed that the tension crack was the dominant factor.Additionally,with increasing weathering thickness,Young's modulus,crack initiation stress,and final strength decreased in constant weathering depth.The results also demonstrated that the failure progress of the numerical models was similar to that observed in the laboratory.展开更多
基金financially supported by the Fundamental Research Funds for the Central Universities (Grant Nos.B200202087 and B200204032)the National Natural Science Foundation of China (Grant No.51609071)。
文摘Rubble mound breakwaters, a prevalent type of sloping breakwater structure, are extensively employed in port and coastal infrastructure projects. Under soft soil foundation conditions, the process of squeezing silt by riprap is implemented to enhance bearing capacity through soft soil replacement and compaction. However, predicting the depth law of squeezing silt by riprap and understanding its mechanism remain significant engineering design challenges.This study employs particle flow code(PFC) based on the discrete element method to simulate the squeezing silt process by riprap, examining variations in depth law under different geological conditions and its mechanical characteristics.Through calibration of the PFC model's meso-parameters via macro-experiments, the study analyzes the effects of riprap size, drop height, and soft soil properties on the depth law of squeezing silt. Findings demonstrate that riprap drop height and soft soil thickness substantially influence the depth, while appropriate calibration of meso-parameters enhances simulation accuracy. This research contributes theoretical and practical guidance for optimizing rubble mound breakwater design, understanding squeezing silt mechanisms, construction practices, and riprap quantity estimation.
基金Project（BK20150005） supported by the Natural Science Foundation of Jiangsu Province for Distinguished Young Scholars, China Project（2015XKZD05） supported by the Fundamental Research Funds for the Central Universities, China
文摘Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass is very important to elucidate the permeability evolution mechanisms, which is difficult to realize in laboratory, but easy to be achieved in numerical simulations. Therefore, the particle flow code （PFC）, a discrete element method, is used to simulate permeability behaviors of rock materials in this study. Owe to the limitation of the existed solid-fluid coupling algorithm in PFC, an improved flow-coupling algorithm is presented to better reflect the preferential flow in rock fractures. The comparative analysis is conducted between original and improved algorithm when simulating rock permeability evolution during triaxial compression, showing that the improved algorithm can better describe the experimental phenomenon. Furthermore, the evolution of pore pressure and flow rate distribution during the flow process are analyzed by using the improved algorithm. It is concluded that during the steady flow process in the fractured specimen, the pore pressure and flow rate both prefer transmitting through the fractures rather than rock matrix. Based on the results, fractures are divided into the following three types： I） fractures link to both the inlet and outlet, II） fractures only link to the inlet, and III） fractures only link to the outlet. The type I fracture is always the preferential propagating path for both the pore pressure and flow rate. For type II fractures, the pore pressure increases and then becomes steady. However, the flow rate increases first and begins to decrease after the flow reaches the stop end of the fracture and finally vanishes. There is no obvious pore pressure or flow rate concentration within type III fractures.
基金The University of Vigo is acknowledged for financing part of the first author’s PhD studiesthe Spanish Ministry of Economy and Competitiveness for funding of the project‘Deepening on the behaviour of rock masses:Scale effects on the stressestrain response of fissured rock samples with particular emphasis on post-failure’,awarded under Contract Reference No.RTI2018-093563-B-I00partially financed by means of European Regional Development Funds from the European Union(EU)。
文摘This study presents a calibration process of three-dimensional particle flow code(PFC3D)simulation of intact and fissured granite samples.First,laboratory stressestrain response from triaxial testing of intact and fissured granite samples is recalled.Then,PFC3D is introduced,with focus on the bonded particle models(BPM).After that,we present previous studies where intact rock is simulated by means of flatjoint approaches,and how improved accuracy was gained with the help of parametric studies.Then,models of the pre-fissured rock specimens were generated,including modeled fissures in the form of“smooth joint”type contacts.Finally,triaxial testing simulations of 1 t 2 and 2 t 3 jointed rock specimens were performed.Results show that both elastic behavior and the peak strength levels are closely matched,without any additional fine tuning of micro-mechanical parameters.Concerning the postfailure behavior,models reproduce the trends of decreasing dilation with increasing confinement and plasticity.However,the dilation values simulated are larger than those observed in practice.This is attributed to the difficulty in modeling some phenomena of fissured rock behaviors,such as rock piece corner crushing with dust production and interactions between newly formed shear bands or axial splitting cracks with pre-existing joints.
基金The National Natural Science Foundation of China(Grants Nos.51309089,41472272 and 11202063)the National Key Technology R&D Program(Grant No.2013BAB06B00)+1 种基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20130846)the Fundamental Research Funds for the Central Universities(Grant No.2015B06014)supported this work
文摘Based on the principle of 3D particle flow code,a numerical landslide run-out model is presented to simulate the failure process of the Zhenggang landslide(in southwestern China) under the effect of water after a rainfall.The relationship between the micro-mechanical parameters and the macro-shear strength of the grain material is determined through numerical calibrations.Then the rainfall effect is considered in numerical simulations and rain-induced sliding processes are performed,which help us to discuss the mechanism of deformation and failure of this landslide together with field observations.It shows the Zhenggang landslide would most likely be activated in Zone I and would gain momentum in Zone II.In order to prevent the potential disaster,a tailing dam is advised to be designed about 175 m downstream from the current landslide boundary of Zone II.Verified by field observations,the presented landslide model can reflect the failure mechanism after rainfall.It can also provide a method to predict the potential disaster and draft disaster prevention measures.
基金the National Key R&D Program of China(No.2017YFC0405002)
文摘Polyurethane polymer grouting materials were studied with conventional triaxial tests via the particle flow code in two dimensions（PFC^（2D）） method, and the simulation results agreed with the experimental data. The particle flow code method can simulate the mechanical properties of the polymer. The triaxial cyclic loading tests of the polymer material under different confining pressures were carried out via PFC^（2D） to analyze its mechanical performance. The PFC^（2D） simulation results show that the value of the elastic modulus of the polymer decreases slowly at first and fluctuated within a narrow range near the value of the peak strength; the cumulative plastic strain increases slowly at first and then increases rapidly; the peak strength and elastic modulus of polymer increase with the confining pressure; the PFC^（2D） method can be used to quantitatively evaluate the damage behavior of the polymer material and estimate the fatigue life of the materials under fatigue load based on the number and the location of micro-cracks. Thus, the PFC^（2D） method is an effective tool to study polymers.
基金financially supported by the National Natural Science Foundation of China(Nos.51874229,52074212,51674188,51504182,51404191,and 51405381)the Natural Science Basic Research Plan of Shaanxi Province of China(Nos.2015JQ5187,2018JQ5183,and 2018JM5161)+3 种基金the Scientific Research Program funded by the Shaanxi Education Department(No.15JK1466)the China Postdoctoral Science Foundation(No.2015M582685)the Outstanding Youth Science Fund of Xi’an University of Science and Technology(No.2018YQ2-01)supported by the National Research Council of Science&and Technology(NST)grant by the Korea Korean government(MSIP)(No.CRC-16-38502-KICT)。
文摘Several special mechanical properties,such as dilatancy and compressibility,of cemented paste backfill(CPB)are controlled by its internal microstructure and evolution.The mesoscopic structure changes of CPB during the development process were investigated.On the basis of the scanning electron microscopy(SEM)and mechanical test results of CPB,the particle size information of CPB was extracted,and a two-dimensional particle flow code(PFC)model of CPB was established to analyze the evolution rule of mesoscopic parameters during CPB development.The embedded FISH language in PFC was used to develop a program for establishing a PFC model on the basis of the SEM results.The mesoscopic parameters of CPB samples at different curing times,such as coordination number(C_(n)),contact force chain,and rose diagram,were obtained by recording and loading and used to analyze the intrinsic relationship between mesoscopic parameter variations and macroscopic mechanical response during CPB development.It is of considerable significance to establish the physical model of CPB using the PFC to reveal the mesoscopic structure of CPB.
基金the National Natural Science Foundation of China (Nos. 41202215 and 40972192)the Innovation Program of Shanghai Municipal Education Commission (No. 11ZZ88)the Innovation Foundation of Shanghai University
文摘A two-dimensionM discrete element code, particle flow code （PFC2D）, is employed to investigate foundations reinforced with horizontal-vertical （H-V） inclusions. The initial states and loading processes of both unreinforced and H-V reinforced foundations are simulated by PFC2D method. The interface between particles and reinforcements, and the reinforcement mechanism of the H-V reinforced foundations are studied through stress distribution graphs, displacement vector graphs and contact force graphs. The simulation results demonstrate that the vertical elements of the H-V reinforcement keep the particles from being displaced under the applied load. The H-V reinforcement can distribute the load uniformly over a wider area, thereby improving the bearing capacity of soil foundation.
文摘The underground or open-pit methods are used for the extraction of mineral resources,each of which is divided into different categories.Coal is one of the mineral resources,which is exploited either by the surface or the underground methods.The long-wall mining is one of the methods for the underground coal mining.In this method,which is a mechanized one,some machines such as the shearer or plow are used for the mining.The coal mine in Parvadeh,Tabas is a mechanized mine that is extracted by the long-wall mining.The modeling with particle flow code software was used in this mine for the evaluation of plow performance using the coal specifications.In this regard,the sample was first calibrated by sampling from the Parvadeh coal mine and performing the uniaxial and Brazilian tests on the model.Then,the modeling was done by constructing the model and using the variables such as the clearance angle and the linear velocity of the plow.After making 28 models for the plow,the best model of the plow was selected based on the maximum force applied to the machine in the X direction.Finally,the results of this study showed that the best plow performance is for a model with the clearance angle of zero and the linear velocity of 9 mm/min,and the maximum force applied to this model is equal to 39,000 kN in the X direction.
文摘This study investigated the effects of weathering depth and thickness on the failure mechanisms of rock samples through experimental and numerical methods.The first configuration involved conducting artificial weathering on limestone using the freezing and thawing(F-T)for 40 cycles.The mechanical parameters of the samples were measured at the end of the 40th cycle.In the second configuration,a series of specimens underwent salt crystallization(S-C)tests for 20 cycles.Experimental results were validated using discrete element method(DEM).Next,the weathered limestone model with dimensions of 108 mm54 mm were prepared.The weathering layers were tested at four different thicknesses(i.e.2.5 mm,5 mm,7.5 mm,and 10 mm)and three different positions(at the surface,5 mm under the rock surface,and 10 mm under the rock surface).According to the results,weathering depth and thickness have a considerable effect on the failure process.The results also showed a correlation between the values of compressive strength and failure mechanisms associated with the weathering layer.The numerical results revealed that the tension crack was the dominant factor.Additionally,with increasing weathering thickness,Young's modulus,crack initiation stress,and final strength decreased in constant weathering depth.The results also demonstrated that the failure progress of the numerical models was similar to that observed in the laboratory.
