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.展开更多
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.展开更多
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.展开更多
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.展开更多
Structure-type rockbursts frequently occur in deep tunnels,with structural planes and stress conditions being critical factors in their formation.In this study,we utilized specially developed analogous materials that ...Structure-type rockbursts frequently occur in deep tunnels,with structural planes and stress conditions being critical factors in their formation.In this study,we utilized specially developed analogous materials that exhibit the high brittleness and strength characteristics of deep hard rock to construct physical models representing different types of structural planes,including composite,exposed,non-exposed,and throughgoing structural planes.Physical simulation experiments were conducted on structuretype rockbursts in deep horseshoe-shaped tunnels,focusing on strain differentiation characteristics,critical triggering conditions,critical crack opening displacement,the incubation process,the reduction effects of structural planes on failure intensity,and formation mechanisms.These experiments were complemented by acoustic and optical monitoring,as well as discrete element numerical simulations,to provide a comprehensive analysis.The results revealed that the most significant strain heterogeneity in the surrounding rock occurs at the tip of the structural plane along the tunnel's minimum principal stress direction,driven by the combined effects of tensile and shear forces.We quantitatively determined the critical stress and strain conditions for structure-type rockbursts and evaluated the intensity of rockbursts induced by different structural planes using critical crack opening displacement(COD)values,the uniformity coefficient,and the curvature coefficient.Analysis of acoustic emission events,including frequency,amplitude,and b-value,indicated that the macro-fracture process is governed by both the principal stress differential and the characteristics of the structural plane.Furthermore,using the bearing capacity reduction coefficient,we found that exposed structural planes have the most significant weakening effect on rock mass strength,followed by non-exposed and throughgoing structural planes.The analysis of average frequency(AF)and rise angle(RA)parameters revealed a close correlation between the failure modes of structure-type rockbursts,the rock mass structure,and the stress levels.These findings provide critical theoretical support for the prediction and prevention of structure-type rockburst disasters.展开更多
沥青混合料是由沥青和粗细集料等粒状单元共同组成的多相复合材料，目前现有的沥青路面结构设计理论均将路面结构当作均匀连续、各向同性的弹性体进行计算分析，这与沥青混合料自身属于颗粒体非连续微细观结构的实际情况不相符。采用离...沥青混合料是由沥青和粗细集料等粒状单元共同组成的多相复合材料，目前现有的沥青路面结构设计理论均将路面结构当作均匀连续、各向同性的弹性体进行计算分析，这与沥青混合料自身属于颗粒体非连续微细观结构的实际情况不相符。采用离散元方法，运用二维颗粒流程序（PFC2D），对比分析了两种颗粒生成方法的特点及其适用条件，以开级配大粒径沥青碎石混合料（Open—graded Large Stone Asphalt Mixes，OLSM）OLSM-25推荐级配的上限为1#级配，建立了1#级配OLSM-25混合料的二维离散元模型，验证了二维离散元方法用于沥青混合料颗粒体材料数值模拟的可行性，弥补了有限元方法的不足之处，为OLSM-25的细观结构分析开辟了新的解决思路和途径。展开更多
基金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 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 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 (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.
基金supported by the National Natural Science Foundation of China(Grant Nos.42307241 and 42107211)the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project(Grant No.SKLGP2022Z008).
文摘Structure-type rockbursts frequently occur in deep tunnels,with structural planes and stress conditions being critical factors in their formation.In this study,we utilized specially developed analogous materials that exhibit the high brittleness and strength characteristics of deep hard rock to construct physical models representing different types of structural planes,including composite,exposed,non-exposed,and throughgoing structural planes.Physical simulation experiments were conducted on structuretype rockbursts in deep horseshoe-shaped tunnels,focusing on strain differentiation characteristics,critical triggering conditions,critical crack opening displacement,the incubation process,the reduction effects of structural planes on failure intensity,and formation mechanisms.These experiments were complemented by acoustic and optical monitoring,as well as discrete element numerical simulations,to provide a comprehensive analysis.The results revealed that the most significant strain heterogeneity in the surrounding rock occurs at the tip of the structural plane along the tunnel's minimum principal stress direction,driven by the combined effects of tensile and shear forces.We quantitatively determined the critical stress and strain conditions for structure-type rockbursts and evaluated the intensity of rockbursts induced by different structural planes using critical crack opening displacement(COD)values,the uniformity coefficient,and the curvature coefficient.Analysis of acoustic emission events,including frequency,amplitude,and b-value,indicated that the macro-fracture process is governed by both the principal stress differential and the characteristics of the structural plane.Furthermore,using the bearing capacity reduction coefficient,we found that exposed structural planes have the most significant weakening effect on rock mass strength,followed by non-exposed and throughgoing structural planes.The analysis of average frequency(AF)and rise angle(RA)parameters revealed a close correlation between the failure modes of structure-type rockbursts,the rock mass structure,and the stress levels.These findings provide critical theoretical support for the prediction and prevention of structure-type rockburst disasters.
文摘沥青混合料是由沥青和粗细集料等粒状单元共同组成的多相复合材料，目前现有的沥青路面结构设计理论均将路面结构当作均匀连续、各向同性的弹性体进行计算分析，这与沥青混合料自身属于颗粒体非连续微细观结构的实际情况不相符。采用离散元方法，运用二维颗粒流程序（PFC2D），对比分析了两种颗粒生成方法的特点及其适用条件，以开级配大粒径沥青碎石混合料（Open—graded Large Stone Asphalt Mixes，OLSM）OLSM-25推荐级配的上限为1#级配，建立了1#级配OLSM-25混合料的二维离散元模型，验证了二维离散元方法用于沥青混合料颗粒体材料数值模拟的可行性，弥补了有限元方法的不足之处，为OLSM-25的细观结构分析开辟了新的解决思路和途径。
