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.展开更多
Load transformation from the yielding part of the soil to the adjacent part is known as the soil arching effect,which plays an important role in the design of various geotechnical infrastructures.Terzaghi’s trapdoor ...Load transformation from the yielding part of the soil to the adjacent part is known as the soil arching effect,which plays an important role in the design of various geotechnical infrastructures.Terzaghi’s trapdoor test was an importantmilestone in the development of theories on soil arching.The research on earth pressure of the trapdoor problem is presented in this paper using the three-dimensional(3D)discrete element method(DEM).Five 3D trapdoor models with different heights are established by 3DDEMsoftware PFC 3D.The variation of earth pressure on the trapdoor with the downward movement of the trapdoor,the distribution of vertical earth pressure along the horizontal direction,the distribution of vertical earth pressure along the vertical direction,the distribution of lateral earth pressure coefficient along the depth direction,the magnitude and direction of contact force chain are studied,respectively.Related research results show that the earth pressure on the trapdoor decreases rapidly after the downward movement of the trapdoor,and then reaches the minimum earth pressure.After that,the earth’s pressure will rise slightly,and whether this phenomenon occurs depends on the depth ratio.For the bottom soil,due to the stress transfer caused by the soil arching effect,the ratio of earth pressure in the loose area decreases,while the ratio of earth pressure in the stable area increases.With the trapdoor moving down,the vertical earth pressure along the depth in the stable zone is basically consistent with the initial state,which shows an approximate linear distribution.After the trapdoor moves down,the distribution of earth pressure along with the depth in the loose area changes,which is far less than the theoretical value of vertical earth pressure of its self-weight.Because of the compression of the soil on both sides,the lateral earth pressure coefficient of most areas on the central axis of the loose zone is close to the passive earth pressure coefficient Kp.The existence of a‘soil arch’can be observed intuitively from the distribution diagram of the contact force chain in the loose zone.展开更多
The micro-mechanical response of asphalt mixtures was studied using the discrete element method. The discrete element sample of stone mastic asphalt was generated first and the vehicle load was applied to the sample. ...The micro-mechanical response of asphalt mixtures was studied using the discrete element method. The discrete element sample of stone mastic asphalt was generated first and the vehicle load was applied to the sample. A user-written program was coded with the FISH language in PFC3 D to extract the contact forces within the sample and the displacements of the particles. Then, the contact forces within the whole sample, in asphalt mastic, in coarse aggregates and between asphalt mastic and coarse aggregates were investigated. Finally, the movement of the particles in the sample was analyzed. The sample was divided into 15 areas and a figure was drawn to show how the balls move in each area according to the displacements of the balls in each area. The displacements of asphalt mastic balls and coarse aggregates were also analyzed. The experimental results explain how the asphalt mixture bears vehicle load and the potential reasons why the rutting forms from a micro-mechanical view.展开更多
Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were ex...Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.展开更多
Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in...Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical φ numerically calculated is less than the one calculated by use of the limit equilibrium method for the sameC. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.展开更多
The application of a combined finite-discrete element modeling approach to simu late the three-dimensional microscopic compaction behavior of single-layer met al powder system was described. The process was treated as...The application of a combined finite-discrete element modeling approach to simu late the three-dimensional microscopic compaction behavior of single-layer met al powder system was described. The process was treated as a static problem,wit h kinematical component being neglected. Due to ill condition,Cholesky’s metho d failed to solve the system equations,while conjugate gradient method was trie d and yielded good results. Deformation of the particles was examined and compar ed with the results of physical modeling experiments. In both cases,the inner p articles were deformed from sphere to polygonal column,with the edges turning f rom arc to straight line. The edge number of a particle was equal to the number of particles surrounding it. And the experiments show that the ductile metal par ticles can be densified only by their plastic deformation without the occurrence of rearrangement phenomenon.展开更多
The accurate understanding of rockburst mechanism poses a global challenge in the field of rock mechanics.Particularly for strain rockburst,achieving self-initiated static-dynamic state transition is a crucial step in...The accurate understanding of rockburst mechanism poses a global challenge in the field of rock mechanics.Particularly for strain rockburst,achieving self-initiated static-dynamic state transition is a crucial step in the formation of catastrophic events.However,the state transition behavior and its impact on rockburst have not received sufficient attention,and are still poorly understood.Therefore,this study specifically focuses on the state transition behavior,aiming to investigate its abrupt transition process and formation mechanism,and triggering effects on rockburst.To facilitate the study,a novel burst rock-surrounding rock combined laboratory test model is proposed and its effectiveness is validated through experiment verification.Subsequently,corresponding numerical models are established using the three-dimensional(3D)discrete element method(DEM),enabling successful simulation of static brittle failure and rockbursts of varying intensities under quasi-static displacement loading conditions.Moreover,through secondary development,comprehensive recording of the mechanical and energy information pertaining to the combined specimen system and its subsystems is achieved.As a result of numerical investigation studies,the elastic rebound dynamic behavior of the surrounding rock was discovered and identified as the key factor triggering rockburst and controlling its intensity.The impact loading on the burst rock,induced by elastic rebound,directly initiates the dynamic processes of rockburst,serving as the direct cause.Additionally,the transient work and energy convergence towards the burst rock resulting from elastic rebound are recognized as the inherent cause of rockburst.Moreover,it has been observed that a larger extent of surrounding rock leads to a stronger elastic rebound,thereby directly contributing to a more intense rockburst.The findings can provide novel theoretical insights for the exploring of rockburst mechanism and the development of monitoring and prevention techniques.展开更多
To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
The mechanical properties of solid oxide fuel cells(SOFCs)can limit their mechanical stability and lifespan.Understanding the correlation between the microstructure and mechanical properties of porous electrode is ess...The mechanical properties of solid oxide fuel cells(SOFCs)can limit their mechanical stability and lifespan.Understanding the correlation between the microstructure and mechanical properties of porous electrode is essential for enhancing the performance and durability of SOFCs.Accurate prediction of mechanical properties of porous electrode can be achieved by microscale finite element modeling based on three-dimensional(3D)microstructures,which requires expensive 3D tomography techniques and massive computational resources.In this study,we proposed a cost-effective alternative approach to access the mechanical properties of porous electrodes,with the elastic properties of La_(0.6)Sr_(0.4)Co_(0.2)Fe_(0.8)O_(3-δc)athode serving as a case study.Firstly,a stochastic modeling was used to reconstruct 3D microstructures from two-dimensional(2D)cross-sections as an alternative to expensive tomography.Then,the discrete element method(DEM)was used to predict the elastic properties of porous ceramics based on the discretized 3D microstructures reconstructed by stochastic modeling.Based on 2D microstructure and the elastic properties calculated by the DEM modeling of the 3D reconstructed porous microstructures,a convolutional neural network(CNN)based deep learning model was built to predict the elastic properties rapidly from 2D microstructures.The proposed combined framework can be implemented with limited computational resources and provide a basis for rapid prediction of mechanical properties and parameter estimation for multiscale modeling of SOFCs.展开更多
Shield tunnel,composed of several segments,is widely used in urban underground engineering.When the tunnel is under load,relative displacement occurs between adjacent segments.In the past,distributed optical fiber sen...Shield tunnel,composed of several segments,is widely used in urban underground engineering.When the tunnel is under load,relative displacement occurs between adjacent segments.In the past,distributed optical fiber sensing technology was used to perform strain monitoring,but there is an urgent need to determine how to transform strain into displacement.In this study,optical frequency domain reflectometry was applied in laboratory tests.Aiming at the shear process and center settlement process of shield tunnel segments,two kinds of quantitative calculation methods were put forward to carry out a quantitative analysis.Meanwhile,the laboratory test process was simulated numerically utilizing the discrete element numerical analysis method.Optical fiber,an atypical geotechnical material,was innovatively applied for discrete element modeling and numerical simulation.The results show that the measured displacement of the dial gauge,the calculated results of the numerical model,and the displacement quantitatively calculated from the optical fiber data agree with each other in general.The latter two methods can potentially be utilized in engineering application of deformation monitoring at shield tunnel joints,but need to be further calibrated and adjusted in detail.展开更多
The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRE...The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.展开更多
Utilizing the Discrete Element Method,this research studied the stiffness distribution of gap-graded soils by modifying the conventional static method.By acknowledging the inherent particle property disparity between ...Utilizing the Discrete Element Method,this research studied the stiffness distribution of gap-graded soils by modifying the conventional static method.By acknowledging the inherent particle property disparity between coarser and finer particles,this research differentiates the stiffness distribution of gap-graded soils from the perspective of contact and particle types.Results indicate that particle property disparity significantly influence the small-strain stiffness characteristics,consequently altering the overall stiffness distribution in gap-graded soil specimens.Specifically,with the equivalent coarser particle property,an increase in particle Young's modulus of finer particles results in an augmentation of small-strain stiffness values,alongside an increased stiffness distribution contribution from finer particles.Nevertheless,this study reveals that even with a higher particle Young's modulus of finer particles,the proportion of small-strain stiffness transferred by finer particles remains consistently lower than their volume fraction.Furthermore,the proportion of stiffness transferred by finer particles may fall below their contribution to stress transmission.This investigation accentuates the subtle yet significant effects of particle property variations on small strain stiffness and its subsequent distribution,providing a foundation for advancing the significance of particle property disparities in evaluating soil responses.展开更多
Discrete element method(DEM)-based numerical models in the YADE environment are used to simulate the constitutive response of uncemented and bio-cemented sands to investigate the influence of boundary conditions,loadi...Discrete element method(DEM)-based numerical models in the YADE environment are used to simulate the constitutive response of uncemented and bio-cemented sands to investigate the influence of boundary conditions,loading and testing conditions,and material types.Both the classical DEM model and the pore scale finite volume(PFV)-coupled DEM model are used to simulate the response of saturated uncemented and lightly cemented sands with a rigid wall boundary under both drained and undrained triaxial compression.A DEM model with flexible boundaries created using particle facet(PFacet)elements is used to simulate undrained triaxial compression of moderately cemented sands,including the influence of confining stress.The PFacet-based model is used to predict the transition from barreling failure to shear banding when the confining stress or the cementation degree increases.The classical DEM model with cohesive bonds of uniform strength is also used to successfully simulate the uniaxial compression response of a sand with an extremely high degree of cementation.Finally,this paper presents a particle-packing model consisting of multiple solid phases for cemented sands based on the understanding that not all particle types will have the same cohesive properties.This multiple solidphase model is a refinement of the classical DEM model that represents the particle physics more realistically,especially for heterogeneous systems.A preliminary parametric study is carried out considering varying cohesive properties and volume fractions for the different solid phases.展开更多
Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie conditio...Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.展开更多
The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction...The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction behavior between hydraulic fractures(HFs) and gravels. A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method. This paper first investigates the dynamic evolution of HFs in glutenite, then analyzes the influences of various factors such as horizontal stress difference(Δσ), minimum horizontal stress(σh), gravel content(Vg), gravel size(dg), and stiffness ratio of gravel to matrix(Rs) on HF propagation geometries. Results show that penetrating the gravel is the primary HF-gravel interaction behavior, which follows sequential and staggered initiation modes. Bypassing the gravel is the secondary behavior, which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(σH). An offset along the gravel boundary is usually formed while penetrating gravels, and the offsets may cause fracture widths to decrease by 37.8%-84.4%. Even if stress dominates the direction of HF propagation, HFs still tend to deflect within gravels. The deviation angle from σH decreases with rising Δσand increases with the increase of dgand Rs. Additionally, intra-gravel shear HFs(IGS-HFs) are prone to be generated in coarse-grained glutenite under high Δσ, while more gravel-bypassing shear HFs(GBSHFs) tend to be created in argillaceous glutenite with high Rsthan in sandy glutenite with low Rs. The findings above prompt the emergence of a novel HF propagation pattern in glutenite, which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.展开更多
Microstructure topology evolution during severe plastic deformation(SPD)is crucial for understanding and optimising the mechanical properties of metallic materials,though its prediction remains challenging.Herein,we c...Microstructure topology evolution during severe plastic deformation(SPD)is crucial for understanding and optimising the mechanical properties of metallic materials,though its prediction remains challenging.Herein,we combine discrete cell complexes(DCC),a fully discrete algebraic topology model-with finite element analysis(FEA)to simulate and analyse the microstructure topology of pure copper under SPD.Using DCC,we model the evolution of microstructure topology characterised by Betti numbers(β_(0),β_(1),β_(2))and Euler characteristic(χ).This captures key changes in GBNs and topological features within representative volume elements(RVEs)containing several hundred grains during SPD-induced recrystallisation.As SPD cycles increase,high-angle grain boundaries(HAGBs)progressively form.Topological analysis reveals an overall decrease in β_(0)values,indicating fewer isolated HAGB substructures,while β_(2) values show a steady upward trend,highlighting new grain formation.Leveraging DCC-derived RVE topology and FEA-generated plastic strain data,we directly simulate the evolution and spatial distribution of microstructure topology and HAGB fraction in a copper tube undergoing cyclic parallel tube channel angular pressing(PTCAP),a representative SPD technique.Within the tube,the HAGB fraction continuously increases with PTCAP cycles,reflecting the microstructure’s gradual transition from subgrains to fully-formed grains.Analysis of Betti number distribution and evolution reveals the microstructural reconstruction mechanism underpinning this subgrain to grain transition during PTCAP.We further demonstrate the significant influence of spatially non-uniform plastic strain distribution on microstructure reconstruction kinetics.This study demonstrates a feasible approach for simulating microstructure topology evolution of metals processed by cyclic SPD via the integration of DCC and FEA.展开更多
Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensiona...Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensional(3D)benchmark model of a laboratory experiment on graywacke to examine the dynamic hydraulic fracturing process under a polyaxial stress state.In the numerical model,injection pressures after breakdown(postbreakdown)are varied to study the impact on fracture growth.The fluid pressure front and crack front are identified in the numerical model to analyze the dynamic relationship between fluid diffusion and fracture propagation.Following the hydraulic fracturing test,the polyaxial stresses are rotated to investigate the influence of the stress field rotation on the fracture slip behavior and permeability.The results show that fracture propagation guides fluid diffusion under a high postbreakdown injection pressure.The crack front runs ahead of the fluid pressure front.Under a low postbreakdown injection pressure,the fluid pressure front gradually reaches the crack front,and fluid diffusion is the main driving factor of fracture propagation.Under polyaxial stress conditions,fluid injection not only opens tensile fractures but also induces hydroshearing.When the polyaxial stress is rotated,the fracture slip direction of a fully extended fracture is consistent with the shear stress direction.The fracture slip direction of a partly extended fracture is influenced by the increase in shear stress.Normal stress affects the permeability evolution by changing the average mechanical aperture.Shear stress can induce shearing and sliding on the fracture plane,thereby increasing permeability.展开更多
Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres wi...Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.展开更多
Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-int...Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-intensity environment,the dynamic failure evolution and instability mechanism of high-steep bedding slopes are simulated via the discrete element method and shaking table test.The dynamic response characteristics and cumulative failure effects of slopes subjected to continuous ground motion are investigated.The results show that the dynamic response characteristics of slopes under continuous earthquakes are influenced by geological and topographic conditions.Elevation has a distinct impact on both the slope interior and surface,with amplification effects more pronounced on the surface.The weak interlayers have different influences on the dynamic amplification effect of slopes.Weak interlayers have dynamic magnification effects on the slope surface at relative elevations of 0-0.33 and 0.82-1.0 but have weakening effects between 0.33 and 0.82.Moreover,the weak interlayers also have controlling effects on the dynamic instability mode of slopes.The characteristics of intergranular contact failure,fracture propagation,and displacement distribution are analyzed to reveal the dynamic failure evolution and instability mechanism through the discrete-element model.The dynamic instability process of slopes includes three stages:fracture initiation(0-0.2g),fracture expansion(0.2g-0.3g),and sliding instability(0.3g-0.6g).This work can provide a valuable reference for the seismic stability and reinforcement of complex slopes.展开更多
基金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.
基金supports from National Natural Science Foundation of China (NSFC Grant No.52008373)Natural Science Foundation of Zhejiang Province of China (No.Q22E080445)are greatly acknowledged.
文摘Load transformation from the yielding part of the soil to the adjacent part is known as the soil arching effect,which plays an important role in the design of various geotechnical infrastructures.Terzaghi’s trapdoor test was an importantmilestone in the development of theories on soil arching.The research on earth pressure of the trapdoor problem is presented in this paper using the three-dimensional(3D)discrete element method(DEM).Five 3D trapdoor models with different heights are established by 3DDEMsoftware PFC 3D.The variation of earth pressure on the trapdoor with the downward movement of the trapdoor,the distribution of vertical earth pressure along the horizontal direction,the distribution of vertical earth pressure along the vertical direction,the distribution of lateral earth pressure coefficient along the depth direction,the magnitude and direction of contact force chain are studied,respectively.Related research results show that the earth pressure on the trapdoor decreases rapidly after the downward movement of the trapdoor,and then reaches the minimum earth pressure.After that,the earth’s pressure will rise slightly,and whether this phenomenon occurs depends on the depth ratio.For the bottom soil,due to the stress transfer caused by the soil arching effect,the ratio of earth pressure in the loose area decreases,while the ratio of earth pressure in the stable area increases.With the trapdoor moving down,the vertical earth pressure along the depth in the stable zone is basically consistent with the initial state,which shows an approximate linear distribution.After the trapdoor moves down,the distribution of earth pressure along with the depth in the loose area changes,which is far less than the theoretical value of vertical earth pressure of its self-weight.Because of the compression of the soil on both sides,the lateral earth pressure coefficient of most areas on the central axis of the loose zone is close to the passive earth pressure coefficient Kp.The existence of a‘soil arch’can be observed intuitively from the distribution diagram of the contact force chain in the loose zone.
基金Funded by the National Natural Science Foundation of China(Nos.51108237 and 51178112)
文摘The micro-mechanical response of asphalt mixtures was studied using the discrete element method. The discrete element sample of stone mastic asphalt was generated first and the vehicle load was applied to the sample. A user-written program was coded with the FISH language in PFC3 D to extract the contact forces within the sample and the displacements of the particles. Then, the contact forces within the whole sample, in asphalt mastic, in coarse aggregates and between asphalt mastic and coarse aggregates were investigated. Finally, the movement of the particles in the sample was analyzed. The sample was divided into 15 areas and a figure was drawn to show how the balls move in each area according to the displacements of the balls in each area. The displacements of asphalt mastic balls and coarse aggregates were also analyzed. The experimental results explain how the asphalt mixture bears vehicle load and the potential reasons why the rutting forms from a micro-mechanical view.
基金Supported by the China National Science and Technology Major Project(2016ZX05033002,2016ZX05033001).
文摘Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.
文摘Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical φ numerically calculated is less than the one calculated by use of the limit equilibrium method for the sameC. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.
文摘The application of a combined finite-discrete element modeling approach to simu late the three-dimensional microscopic compaction behavior of single-layer met al powder system was described. The process was treated as a static problem,wit h kinematical component being neglected. Due to ill condition,Cholesky’s metho d failed to solve the system equations,while conjugate gradient method was trie d and yielded good results. Deformation of the particles was examined and compar ed with the results of physical modeling experiments. In both cases,the inner p articles were deformed from sphere to polygonal column,with the edges turning f rom arc to straight line. The edge number of a particle was equal to the number of particles surrounding it. And the experiments show that the ductile metal par ticles can be densified only by their plastic deformation without the occurrence of rearrangement phenomenon.
基金supported by the National Natural Science Foundation of China(Grant Nos.U22A20597,and 42142019)the“Unveiling and Commanding”Project of Science and Technology Program of Tibet(Grant No.XZ202303ZY0006G)the Shanghai Peak Plateau Discipline(Class I).
文摘The accurate understanding of rockburst mechanism poses a global challenge in the field of rock mechanics.Particularly for strain rockburst,achieving self-initiated static-dynamic state transition is a crucial step in the formation of catastrophic events.However,the state transition behavior and its impact on rockburst have not received sufficient attention,and are still poorly understood.Therefore,this study specifically focuses on the state transition behavior,aiming to investigate its abrupt transition process and formation mechanism,and triggering effects on rockburst.To facilitate the study,a novel burst rock-surrounding rock combined laboratory test model is proposed and its effectiveness is validated through experiment verification.Subsequently,corresponding numerical models are established using the three-dimensional(3D)discrete element method(DEM),enabling successful simulation of static brittle failure and rockbursts of varying intensities under quasi-static displacement loading conditions.Moreover,through secondary development,comprehensive recording of the mechanical and energy information pertaining to the combined specimen system and its subsystems is achieved.As a result of numerical investigation studies,the elastic rebound dynamic behavior of the surrounding rock was discovered and identified as the key factor triggering rockburst and controlling its intensity.The impact loading on the burst rock,induced by elastic rebound,directly initiates the dynamic processes of rockburst,serving as the direct cause.Additionally,the transient work and energy convergence towards the burst rock resulting from elastic rebound are recognized as the inherent cause of rockburst.Moreover,it has been observed that a larger extent of surrounding rock leads to a stronger elastic rebound,thereby directly contributing to a more intense rockburst.The findings can provide novel theoretical insights for the exploring of rockburst mechanism and the development of monitoring and prevention techniques.
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
基金supported by the National Natural Science Foundation of China(Grant Nos.12172104 and 11932005)the Talent Recruitment Project of Guangdong(2021QN02L892)+3 种基金the Stable Supporting Fund of Shenzhen(GXWD20231130153335002)the Shccig-Qinling Program(SMYJY202300140C)the program of Innovation Team in Universities and Colleges in Guangdong(2021KCXTD006)Development and Reform Commission of Shenzhen(XMHT20220103004).
文摘The mechanical properties of solid oxide fuel cells(SOFCs)can limit their mechanical stability and lifespan.Understanding the correlation between the microstructure and mechanical properties of porous electrode is essential for enhancing the performance and durability of SOFCs.Accurate prediction of mechanical properties of porous electrode can be achieved by microscale finite element modeling based on three-dimensional(3D)microstructures,which requires expensive 3D tomography techniques and massive computational resources.In this study,we proposed a cost-effective alternative approach to access the mechanical properties of porous electrodes,with the elastic properties of La_(0.6)Sr_(0.4)Co_(0.2)Fe_(0.8)O_(3-δc)athode serving as a case study.Firstly,a stochastic modeling was used to reconstruct 3D microstructures from two-dimensional(2D)cross-sections as an alternative to expensive tomography.Then,the discrete element method(DEM)was used to predict the elastic properties of porous ceramics based on the discretized 3D microstructures reconstructed by stochastic modeling.Based on 2D microstructure and the elastic properties calculated by the DEM modeling of the 3D reconstructed porous microstructures,a convolutional neural network(CNN)based deep learning model was built to predict the elastic properties rapidly from 2D microstructures.The proposed combined framework can be implemented with limited computational resources and provide a basis for rapid prediction of mechanical properties and parameter estimation for multiscale modeling of SOFCs.
基金National Natural Science Foundation of China,Grant/Award Numbers:41977218,42222707State Key Laboratory for GeoMechanics and Deep Underground Engineering,Grant/Award Number:SKLGDUEK2117。
文摘Shield tunnel,composed of several segments,is widely used in urban underground engineering.When the tunnel is under load,relative displacement occurs between adjacent segments.In the past,distributed optical fiber sensing technology was used to perform strain monitoring,but there is an urgent need to determine how to transform strain into displacement.In this study,optical frequency domain reflectometry was applied in laboratory tests.Aiming at the shear process and center settlement process of shield tunnel segments,two kinds of quantitative calculation methods were put forward to carry out a quantitative analysis.Meanwhile,the laboratory test process was simulated numerically utilizing the discrete element numerical analysis method.Optical fiber,an atypical geotechnical material,was innovatively applied for discrete element modeling and numerical simulation.The results show that the measured displacement of the dial gauge,the calculated results of the numerical model,and the displacement quantitatively calculated from the optical fiber data agree with each other in general.The latter two methods can potentially be utilized in engineering application of deformation monitoring at shield tunnel joints,but need to be further calibrated and adjusted in detail.
基金funded by the National Natural Science Foundation of China(No.51978048).
文摘The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.
基金Financial supports from the PolyU Distinguished Postdoctoral Fellowship Scheme are highly appreciatedsupported by the National Natural Science Foundation of China (Grant No.52201008)the Fundamental Research Funds for the Central Universities,the State Key Laboratory of Particle Detection and Electronics (Grant No.SKLPDE-KF-202311).
文摘Utilizing the Discrete Element Method,this research studied the stiffness distribution of gap-graded soils by modifying the conventional static method.By acknowledging the inherent particle property disparity between coarser and finer particles,this research differentiates the stiffness distribution of gap-graded soils from the perspective of contact and particle types.Results indicate that particle property disparity significantly influence the small-strain stiffness characteristics,consequently altering the overall stiffness distribution in gap-graded soil specimens.Specifically,with the equivalent coarser particle property,an increase in particle Young's modulus of finer particles results in an augmentation of small-strain stiffness values,alongside an increased stiffness distribution contribution from finer particles.Nevertheless,this study reveals that even with a higher particle Young's modulus of finer particles,the proportion of small-strain stiffness transferred by finer particles remains consistently lower than their volume fraction.Furthermore,the proportion of stiffness transferred by finer particles may fall below their contribution to stress transmission.This investigation accentuates the subtle yet significant effects of particle property variations on small strain stiffness and its subsequent distribution,providing a foundation for advancing the significance of particle property disparities in evaluating soil responses.
基金support for this study from National Science Foundation(NSF)under the Engineering Research Centers(ERC)program,grant EEC-1449501.
文摘Discrete element method(DEM)-based numerical models in the YADE environment are used to simulate the constitutive response of uncemented and bio-cemented sands to investigate the influence of boundary conditions,loading and testing conditions,and material types.Both the classical DEM model and the pore scale finite volume(PFV)-coupled DEM model are used to simulate the response of saturated uncemented and lightly cemented sands with a rigid wall boundary under both drained and undrained triaxial compression.A DEM model with flexible boundaries created using particle facet(PFacet)elements is used to simulate undrained triaxial compression of moderately cemented sands,including the influence of confining stress.The PFacet-based model is used to predict the transition from barreling failure to shear banding when the confining stress or the cementation degree increases.The classical DEM model with cohesive bonds of uniform strength is also used to successfully simulate the uniaxial compression response of a sand with an extremely high degree of cementation.Finally,this paper presents a particle-packing model consisting of multiple solid phases for cemented sands based on the understanding that not all particle types will have the same cohesive properties.This multiple solidphase model is a refinement of the classical DEM model that represents the particle physics more realistically,especially for heterogeneous systems.A preliminary parametric study is carried out considering varying cohesive properties and volume fractions for the different solid phases.
基金a U.S. Federal Railroad Administration (FRA)BAA project,titled “Mitigation of Differential Movement at Railway Transitions for High-Speed Passenger Rail and Joint Passenger/Freight Corridors”the financial support provided by the China Scholarship Council (CSC),which funded Zhongyi Liu’s and Wenjing Li’s time and research efforts for this study
文摘Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.
基金supported by the National Natural Science Foundation of China(Grant No.52304003)the Natural Science Foundation of Sichuan Province(Grant No.2024NSFSC0961)the Postdoctoral Fellowship Program of CPSF(Grant No.GZB20230090).
文摘The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels, making it challenging to perform hydraulic fracturing effectively. To solve this issue, it is essential to study interaction behavior between hydraulic fractures(HFs) and gravels. A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method. This paper first investigates the dynamic evolution of HFs in glutenite, then analyzes the influences of various factors such as horizontal stress difference(Δσ), minimum horizontal stress(σh), gravel content(Vg), gravel size(dg), and stiffness ratio of gravel to matrix(Rs) on HF propagation geometries. Results show that penetrating the gravel is the primary HF-gravel interaction behavior, which follows sequential and staggered initiation modes. Bypassing the gravel is the secondary behavior, which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(σH). An offset along the gravel boundary is usually formed while penetrating gravels, and the offsets may cause fracture widths to decrease by 37.8%-84.4%. Even if stress dominates the direction of HF propagation, HFs still tend to deflect within gravels. The deviation angle from σH decreases with rising Δσand increases with the increase of dgand Rs. Additionally, intra-gravel shear HFs(IGS-HFs) are prone to be generated in coarse-grained glutenite under high Δσ, while more gravel-bypassing shear HFs(GBSHFs) tend to be created in argillaceous glutenite with high Rsthan in sandy glutenite with low Rs. The findings above prompt the emergence of a novel HF propagation pattern in glutenite, which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.
基金support from Outstanding Youth Fund of Jiangsu Province(BK20240077)Key Project(Provincial-Municipal Joint)of Jiangsu Province(BK20243044)+2 种基金Fundamental Research Funds for the Central Universities(NE2024001)National Youth Talents Programof Chinaa project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions.
文摘Microstructure topology evolution during severe plastic deformation(SPD)is crucial for understanding and optimising the mechanical properties of metallic materials,though its prediction remains challenging.Herein,we combine discrete cell complexes(DCC),a fully discrete algebraic topology model-with finite element analysis(FEA)to simulate and analyse the microstructure topology of pure copper under SPD.Using DCC,we model the evolution of microstructure topology characterised by Betti numbers(β_(0),β_(1),β_(2))and Euler characteristic(χ).This captures key changes in GBNs and topological features within representative volume elements(RVEs)containing several hundred grains during SPD-induced recrystallisation.As SPD cycles increase,high-angle grain boundaries(HAGBs)progressively form.Topological analysis reveals an overall decrease in β_(0)values,indicating fewer isolated HAGB substructures,while β_(2) values show a steady upward trend,highlighting new grain formation.Leveraging DCC-derived RVE topology and FEA-generated plastic strain data,we directly simulate the evolution and spatial distribution of microstructure topology and HAGB fraction in a copper tube undergoing cyclic parallel tube channel angular pressing(PTCAP),a representative SPD technique.Within the tube,the HAGB fraction continuously increases with PTCAP cycles,reflecting the microstructure’s gradual transition from subgrains to fully-formed grains.Analysis of Betti number distribution and evolution reveals the microstructural reconstruction mechanism underpinning this subgrain to grain transition during PTCAP.We further demonstrate the significant influence of spatially non-uniform plastic strain distribution on microstructure reconstruction kinetics.This study demonstrates a feasible approach for simulating microstructure topology evolution of metals processed by cyclic SPD via the integration of DCC and FEA.
基金supported by the Knowledge Innovation Program of Wuhan-Basic Research (Grant No.2022010801010159)support from the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES (Contract number VH-NG-1516)supported by the Swedish Radiation Safety Authority (Project SSM2020-2758).
文摘Understanding the hydromechanical behavior and permeability stress sensitivity of hydraulic fractures is fundamental for geotechnical applications associated with fluid injection.This paper presents a three-dimensional(3D)benchmark model of a laboratory experiment on graywacke to examine the dynamic hydraulic fracturing process under a polyaxial stress state.In the numerical model,injection pressures after breakdown(postbreakdown)are varied to study the impact on fracture growth.The fluid pressure front and crack front are identified in the numerical model to analyze the dynamic relationship between fluid diffusion and fracture propagation.Following the hydraulic fracturing test,the polyaxial stresses are rotated to investigate the influence of the stress field rotation on the fracture slip behavior and permeability.The results show that fracture propagation guides fluid diffusion under a high postbreakdown injection pressure.The crack front runs ahead of the fluid pressure front.Under a low postbreakdown injection pressure,the fluid pressure front gradually reaches the crack front,and fluid diffusion is the main driving factor of fracture propagation.Under polyaxial stress conditions,fluid injection not only opens tensile fractures but also induces hydroshearing.When the polyaxial stress is rotated,the fracture slip direction of a fully extended fracture is consistent with the shear stress direction.The fracture slip direction of a partly extended fracture is influenced by the increase in shear stress.Normal stress affects the permeability evolution by changing the average mechanical aperture.Shear stress can induce shearing and sliding on the fracture plane,thereby increasing permeability.
基金the financial support provided by Tianfu Yongxing Laboratory Organized Research Project Funding(No.2023CXXM01)the ARC linkage program(No.LP200100420).
文摘Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.
基金Project(52108361)supported by the National Natural Science Foundation of ChinaProjects(BK20231217,BK20220265)supported by the Basic Research Program of Jiangsu Province,China+5 种基金Project(sklhse-KF-2025-D-02)supported by the Open Research Fund Program of the State Key Laboratory of Hydroscience and Engineering,ChinaProject(2023ZB15)supported by the Independent Research Project of the State Key Laboratory of Subtropical Building and Urban Science,ChinaProject(SKLGME023001)supported by the Key Laboratory of Geomechanics and Geotechnical Engineering Safety,the Chinese Academy of SciencesProject(2025A04J3992)supported by the Basic and Applied Basic Research Project of the Guangzhou Science and Technology Bureau,ChinaProject(SKLGP2022Z015)supported by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project,ChinaProjects(2023YFS0436,2024NSFSC1715)supported by the Science and Technology Department of Sichuan Province,China。
文摘Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-intensity environment,the dynamic failure evolution and instability mechanism of high-steep bedding slopes are simulated via the discrete element method and shaking table test.The dynamic response characteristics and cumulative failure effects of slopes subjected to continuous ground motion are investigated.The results show that the dynamic response characteristics of slopes under continuous earthquakes are influenced by geological and topographic conditions.Elevation has a distinct impact on both the slope interior and surface,with amplification effects more pronounced on the surface.The weak interlayers have different influences on the dynamic amplification effect of slopes.Weak interlayers have dynamic magnification effects on the slope surface at relative elevations of 0-0.33 and 0.82-1.0 but have weakening effects between 0.33 and 0.82.Moreover,the weak interlayers also have controlling effects on the dynamic instability mode of slopes.The characteristics of intergranular contact failure,fracture propagation,and displacement distribution are analyzed to reveal the dynamic failure evolution and instability mechanism through the discrete-element model.The dynamic instability process of slopes includes three stages:fracture initiation(0-0.2g),fracture expansion(0.2g-0.3g),and sliding instability(0.3g-0.6g).This work can provide a valuable reference for the seismic stability and reinforcement of complex slopes.
