Microwave-assisted rock-breaking technology,as a novel hybrid approach,is anticipated to facilitate the efficient excavation of complex rock formations.It is therefore crucial to understand the damage and failure mech...Microwave-assisted rock-breaking technology,as a novel hybrid approach,is anticipated to facilitate the efficient excavation of complex rock formations.It is therefore crucial to understand the damage and failure mechanisms of rocks that have been subjected to irradiation.In this study,uniaxial compression experiments were conducted on granite specimens after 1.4 kW microwave irradiation for varying durations.Furthermore,a numerical method was proposed to solve electromagnetic-thermal-mechanical coupling problems by integrating finite and discrete elements.The results demonstrated a differential temperature distribution(high temperature in the middle and low-temperature areas at the ends)in the granite specimens under microwave irradiation,which resulted in a notable reduction in their physical and mechanical properties.As the duration of irradiation increased,the rate of heating and the extent of strength reduction both diminished,while the morphology and distribution of cracks at ultimate failure became increasingly complex.The numerical method effectively addresses the simulation challenges associated with the electromagnetic selective heating of granite containing multiple polar minerals under microwave irradiation.This approach accounted for the non-uniform thermal expansion of the minerals and provided a comprehensive model of damage progression under compression.展开更多
Rockbursting in deep tunnelling is a complex phenomenon posing significant challenges both at the design and construction stages of an underground excavation within hard rock masses and under high in situ stresses. Wh...Rockbursting in deep tunnelling is a complex phenomenon posing significant challenges both at the design and construction stages of an underground excavation within hard rock masses and under high in situ stresses. While local experience, field monitoring, and informed data-rich analysis are some of the tools commonly used to manage the hazards and the associated risks, advanced numerical techniques based on discontinuum modelling have also shown potential in assisting in the assessment of rockbursting. In this study, the hybrid finite-discrete element method(FDEM) is employed to investigate the failure and fracturing processes, and the mechanisms of energy storage and rapid release resulting in bursting, as well as to assess its utility as part of the design process of underground excavations.Following the calibration of the numerical model to simulate a deep excavation in a hard, massive rock mass, discrete fracture network(DFN) geometries are integrated into the model in order to examine the impact of rock structure on rockbursting under high in situ stresses. The obtained analysis results not only highlight the importance of explicitly simulating pre-existing joints within the model, as they affect the mobilised failure mechanisms and the intensity of strain bursting phenomena, but also show how the employed joint network geometry, the field stress conditions, and their interaction influence the extent and depth of the excavation induced damage. Furthermore, a rigorous analysis of the mass and velocity of the ejected rock blocks and comparison of the obtained data with well-established semi-empirical approaches demonstrate the potential of the method to provide realistic estimates of the kinetic energy released during bursting for determining the energy support demand.展开更多
Mechanical cutting provides one of the most flexible and environmentally friendly excavation methods.It has attracted numerous efforts to model the rock chipping and fragmentation process,especially using the explicit...Mechanical cutting provides one of the most flexible and environmentally friendly excavation methods.It has attracted numerous efforts to model the rock chipping and fragmentation process,especially using the explicit finite element method(FEM) and bonded particle model(BPM),in order to improve cutting efficiency.This study investigates the application of a general-purpose graphic-processing-unit parallelised hybrid finite-discrete element method(FDEM) which enjoys the advantages of both explicit FEM and BPM,in modelling the rock chipping and fragmentation process in the rock scratch test of mechanical rock cutting.The input parameters of FDEM are determined through a calibration procedure of modelling conventional Brazilian tensile and uniaxial compressive tests of limestone,A series of scratch tests with various cutting velocities,cutter rake angles and cutting depths is then modelled using FDEM with calibrated input parameters.A few cycles of cutter/rock interactions,including their engagement and detachment process,are modelled for each case,which is conducted for the first time to the best knowledge of the authors,thanks to the general purpose graphic processing units(GPGPU) parallelisation.The failure mechanism,cutting force,chipping morphology and effect of various factors on them are discussed on the basis of the modelled results.Finally,it is concluded that GPGPU-parallelised FDEM provides a powerful tool to further study rock cutting and improve cutting efficiencies since it can explicitly capture different fracture mechanisms contributing to the rock chipping as well as chip formation and the separation process in mechanical cutting.Moreover,it is concluded that chipping is mostly owed to the mix-mode Ⅰ-Ⅱ fracture in all cases although mode Ⅱ cracks and mode Ⅰ cracks are the dominant failures in rock cutting with shallow and deep cutting depths,respectively.The chip morphology is found to be a function of cutter velocdty,cutting depth and cutter rake angle.展开更多
Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed ...Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.展开更多
Due to the long construction life,improper design methods,brittle material properties and poor construction techniques,most existing masonry structures do not perform well during earthquakes.The retrofitting method us...Due to the long construction life,improper design methods,brittle material properties and poor construction techniques,most existing masonry structures do not perform well during earthquakes.The retrofitting method using an external steel-meshed mortar layer is widely used to retrofit existing masonry buildings.Assessing the seismic performance of masonry walls reinforced by an external steel-meshed mortar layer reasonably and effectively is a difficult subject in the research field of masonry structures.Based on the combined finite-discrete elements method,the numerical models of retrofitted brick walls with four different masonry mortar strengths by an external mortar layer are established.The shear strength of mortar and the contact between the retrofitted mortar layer and the brick blocks are discussed in detail.The failure patterns and load-displacement curves of the retrofitted brick walls were obtained by applying low cycle reciprocating loads to the numerical model,and the bearing capacity and the failure mechanism of the retrofitted walls were obtained by comparing the failure patterns,ultimate bearing capacity,deformability and other aspects with the tests.This study provides a basis for improving the seismic strengthening design method of masonry structures and helps to better assess the seismic performance of masonry structures after retrofitting.展开更多
This paper presents the application of a hybrid finite-discrete element method to study blast-induceddamage in circular tunnels. An extensive database of field tests of underground explosions above tunnelsis used for ...This paper presents the application of a hybrid finite-discrete element method to study blast-induceddamage in circular tunnels. An extensive database of field tests of underground explosions above tunnelsis used for calibrating and validating the proposed numerical method; the numerical results areshown to be in good agreement with published data for large-scale physical experiments. The method isthen used to investigate the influence of rock strength properties on tunnel durability to withstand blastloads. The presented analysis considers blast damage in tunnels excavated through relatively weak（sandstone） and strong （granite） rock materials. It was found that higher rock strength will increase thetunnel resistance to the load on one hand, but decrease attenuation on the other hand. Thus, undercertain conditions, results for weak and strong rock masses are similar. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid ...Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Soil desiccation cracking is a common phenomenon on the earth surface.Numerical modeling is an effective approach to study the desiccation cracking mechanism of soil.This work develops a novel 3D moisture diffusion di...Soil desiccation cracking is a common phenomenon on the earth surface.Numerical modeling is an effective approach to study the desiccation cracking mechanism of soil.This work develops a novel 3D moisture diffusion discrete model that is capable of dynamically assessing the effect of cracking on moisture diffusion and allowing moisture to be discontinuous on both sides of the cracks.Then,the parametric analysis of the moisture exchange coefficient in the 3D moisture diffusion discrete model is carried out for moisture diffusion in continuous media,and the selection criterion of the moisture exchange coefficient for the unbroken cohesive element is given.Subsequently,an example of moisture migration in a medium with one crack is provided to illustrate the crack hindering effect on moisture migration.Finally,combining the 3D moisture diffusion discrete model with the finite-discrete element method(FDEM),the moisture diffusion-fracture coupling model is built to study the desiccation cracking in a strip soil and the crack pattern of a rectangular soil.The evolution of crack area and volume with moisture content is quantitatively analyzed.The modeling number and average width of cracks in the strip soil show a good consistency with the experimental results,and the crack pattern of the rectangular soil matches well with the existing numerical results,validating the coupled moisture diffusion-fracture model.Additionally,the parametric study of soil desiccation cracking is performed.The developed model offers a powerful tool for exploring soil desiccation cracking.展开更多
Based on the finite-discrete element method,a three-dimensional numerical model for axial impact rock breaking was established and validated.A computational method for energy conversion during impact rock breaking was...Based on the finite-discrete element method,a three-dimensional numerical model for axial impact rock breaking was established and validated.A computational method for energy conversion during impact rock breaking was proposed,and the effects of conical tooth forward rake angle,rock temperature,and impact velocity on rock breaking characteristics and energy transfer laws were analyzed.The results show that during single impact rock breaking with conical tooth bits,merely 7.52%to 12.51%of the energy is utilized for rock breaking,while a significant 57.26%to 78.10%is dissipated as frictional loss.An insufficient forward rake angle increases tooth penetration depth and frictional loss,whereas an excessive forward rake angle reduces penetration capability,causing bit rebound and greater energy absorption by the drill rod.Thus,an optimal forward rake angle exists.Regarding environmental factors,high temperatures significantly enhance impact-induced rock breaking.Thermal damage from high temperatures reduces rock strength and inhibits its energy absorption.Finally,higher impact velocities intensify rock damage,yet excessively high velocities increase frictional loss and reduce the proportion of energy absorbed by the rock,thereby failing to substantially improve rock breaking efficiency.An optimal impact velocity exists.展开更多
For rapid and cost-effective hammer drilling,accurate prediction of rock impact response is crucial for designing optimal bits and maximising rock fragmentation.Current design optimisation workflows combine numerical ...For rapid and cost-effective hammer drilling,accurate prediction of rock impact response is crucial for designing optimal bits and maximising rock fragmentation.Current design optimisation workflows combine numerical simulations and experiments but often require numerous iterations to pinpoint the optimal design.Although physics-based models can potentially reduce experimental expenses,their significant computational demands present challenges when simulating the complex fragmentation dynamics during drill bit-rock interactions.This study introduces a data-driven artificial intelligence(AI)model,employing a multilayer perceptron(MLP)as a surrogate.The model leverages the hybrid finitediscrete element model(FDEM)as a powerful method in rock fracture mechanics to generate a sufficiently large training dataset.An automated workflow has been developed for generating the training data,comprising a pipeline that includes pre-processing,solving,and post-processing modules.Subsequently,the AI models were integrated into an optimisation framework alongside uncertainty quantification to demonstrate their potential in enhancing drilling efficiency through optimised bit design and operations.The MLP exhibits high accuracy in predicting key parameters,including rebound velocity,total crack length,quantities of fragments with different sizes and maximum contact force between rock and insert.Notably,this approach achieves real-time prediction compared to the 5-7 min simulation times of FDEM.Integrating this data-driven model into a design framework enables rapid assessment of different bit designs under various operational conditions.More broadly,this approach has the potential to impact other applications,such as digital twins,serving as a forward and inverse model for predicting rock type and optimising drilling performance.展开更多
Particle morphology is critical in affecting the crushing behavior of rockfill materials.In contrast,most current single particle simulations lack satisfactory morphology accuracy,and the resulting crushing modes devi...Particle morphology is critical in affecting the crushing behavior of rockfill materials.In contrast,most current single particle simulations lack satisfactory morphology accuracy,and the resulting crushing modes deviate from observations to some extent.Therefore,we reconstruct the real particle morphology with the spherical harmonic(SH)method and employ the finite-discrete element method(FDEM)to simulate the one-dimensional(1D)compressive crushing process of basalt particles commonly used in rockfill.The influences of four main morphological parameters,i.e.sphericity,aspect ratio,roundness,and convexity,on the single particle strength and the crushing modes are discussed.The results show that with the SH degree set to 15 and a mesh number of 20,480,the FDEM models of reconstructed particles achieve sufficient morphology accuracy and high computational efficiency.Based on the model,the simulation results demonstrate that the aspect ratio has the most significant impact on single particle strength,followed by sphericity.In contrast,roundness and convexity have a weaker effect than the above two parameters.Also,it is revealed that single particle strength decreases with increasing aspect ratio and sphericity,while it increases with higher roundness and convexity.Furthermore,aspect ratio significantly changes the initial crushing position,sphericity dominates post-crushing fragment size and quantity,and roundness mainly affects post-crushing morphology.The model results have been employed in establishing a support vector regression(SVR)-based predicted model,exhibiting good predictive performance and advantages for the optimization of rockfill particles in engineering.展开更多
Frost heave in water-bearing rock masses poses significant threats to geotechnical engineering.This paper developed a novel three-dimensional(3D)frost model,based on the combined finite-discrete element method(FDEM),t...Frost heave in water-bearing rock masses poses significant threats to geotechnical engineering.This paper developed a novel three-dimensional(3D)frost model,based on the combined finite-discrete element method(FDEM),to investigate the frost heave process in rock masses where thermal transfer,water migration,water-ice phase transition(ice growth)and ice-rock interaction are explicitly simulated.The proposed model is first validated against existing experimental and analytical solutions,and further applied to investigate path-dependent frost heave behavior under various freezing conditions.Results show that freezing direction plays a vital role in the dynamic ice growth and ice-rock interaction,thus affecting the frost heave behavior.In the top-down freezing regime,ice plugs form first at the crack's top surface,sealing the crack and preventing water migration,which can amplify ice pressure.Parametric studies,including rock Young's modulus,ice-rock friction,and rock hydraulic conductivity,further reveal that the temporal aspects of ice development and rock mechanical response strongly affect ice-rock interaction and hence the frost heave mechanism.Furthermore,some typical phenomena(e.g.water/ice extrusion and frost cracking)can also be well captured in this model.This novel numerical framework sheds new light on frost heave behavior and enriches our understanding of frost heave mechanisms and ice-rock interaction processes within cold environment engineering projects.展开更多
The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly ...The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.展开更多
This paper presents the development of a coupled modeling approach to simulate cryogenic thermo-hydro-mechanical(THM)processes associated with a freezing medium,which is then implemented in the combined finite-discret...This paper presents the development of a coupled modeling approach to simulate cryogenic thermo-hydro-mechanical(THM)processes associated with a freezing medium,which is then implemented in the combined finite-discrete element method code(FDEM)for multi-physics simulation.The governing equations are deduced based on energy and mass conservation,and static equilibrium equations,considering water/ice phase change,where the strong couplings between multi-fields are supplemented by critical coupling parameters(e.g.unfrozen water content,permeability,and thermal conductivity).The proposed model is validated against laboratory and field experiments.Results show that the cryogenic THM model can well predict the evolution of strongly coupled processes observed in frozen media(e.g.heat transfer,water migration,and frost heave deformation),while also capturing,as emergent properties of the model,important phenomena(e.g.latent heat,cryogenic suction,ice expansion and distinct three-zone distribution)caused by water/ice phase change at laboratory and field scales,which are difficult to be all revealed by existing THM models.The novel modeling framework presents a gateway to further understanding and predicting the multi-physical coupling behavior of frozen media in cold regions.展开更多
A new thermomechanical(TM)coupled finite-discrete element method(FDEM)model,incorporating heat conduction,thermal cracking,and contact heat transfer,has been proposed for both continuous and discontinuous geomaterials...A new thermomechanical(TM)coupled finite-discrete element method(FDEM)model,incorporating heat conduction,thermal cracking,and contact heat transfer,has been proposed for both continuous and discontinuous geomaterials.This model incorporates a heat conduction model that can accurately calculate the thermal field in continuousediscontinuous transition processes within a finite element framework.A modified contact heat transfer model is also included,which accounts for the entire contact area of discrete bodies.To align with the finite strain theory utilized in the FDEM mechanics module,the TM coupling module in the model is based on the multiplicative decomposition of the deformation gradient.The proposed model has been applied to various scenarios,including heat conduction in both continuous and discontinuous media during transient states,thermal-induced strain and stress,and thermal cracking conditions.The thermal field calculation model and the TM coupling model have been validated by comparing the numerical results with experiment findings and analytical solutions.These numerical cases demonstrate the reliability of the proposed model convincingly,making it suitable for use across a wide range of continuous and discontinuous media.展开更多
This paper introduces a deep learning workflow to predict phase distributions within complex geometries during two-phase capillary-dominated drainage.We utilize subsamples from Computerized Tomography(CT)images of roc...This paper introduces a deep learning workflow to predict phase distributions within complex geometries during two-phase capillary-dominated drainage.We utilize subsamples from Computerized Tomography(CT)images of rocks and incorporate pixel size,interfacial tension,contact angle,and pressure as inputs.First,an efficient morphology-based simulator creates a diverse dataset of phase distributions.Then,two commonly used convolutional and recurrent neural networks are explored and their deficiencies are highlighted,particularly in capturing phase connectivity.Subsequently,we develop a Higher-Dimensional Vision Transformer(HD-ViT)that drains pores solely based on their size,with phase connectivity enforced as a post-processing step.This enables inference for images of varying sizes,resolutions,and inlet-outlet setup.After training on a massive dataset of over 9.5 million instances,HD-ViT achieves excellent performance.We demonstrate the accuracy and speed advantage of the model on new and larger sandstone and carbonate images.We further evaluate HD-ViT against experimental fluid distribution images and the corresponding Lattice-Boltzmann simulations,producing similar outcomes in a matter of seconds.In the end,we train and validate a 3D version of the model.展开更多
Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-ex...Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-existing joints in the damage evolution around the underground opening is of critical importance as they govern the fracturing mechanisms and influence the brittle responses of these hard rock masses under highly anisotropic in situ stresses.In this study,the main focus is the impact of joint network geometry,joint strength and applied field stresses on the rock mass behaviours and the evolution of excavation induced damage due to the loss of confinement as a tunnel face advances.Analysis of such a phenomenon was conducted using the finite-discrete element method(FDEM).The numerical model is initially calibrated in order to match the behaviour of the fracture-free,massive Lac du Bonnet granite during the excavation of the Underground Research Laboratory(URL)Test Tunnel,Canada.The influence of the pre-existing joints on the rock mass response during excavation is investigated by integrating discrete fracture networks(DFNs)of various characteristics into the numerical models under varying in situ stresses.The numerical results obtained highlight the significance of the pre-existing joints on the reduction of in situ rock mass strength and its capacity for extension with both factors controlling the brittle response of the material.Furthermore,the impact of spatial distribution of natural joints on the stability of an underground excavation is discussed,as well as the potentially minor influence of joint strength on the stress induced damage within joint systems of a non-persistent nature under specific conditions.Additionally,the in situ stress-joint network interaction is examined,revealing the complex fracturing mechanisms that may lead to uncontrolled fracture propagation that compromises the overall stability of an underground excavation.展开更多
The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fractur...The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.展开更多
A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevo...A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.展开更多
Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage deve...Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method （FDEM） code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone （EDZ） around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
基金funded by the Postgraduate Research and Practice Innovation Program of Jiangsu Province(Grant No.KYCX23_2744)the Fundamental Research Funds for the Central Universities(Grant No.2023XSCX051)the Graduate Innovation Program of China University of Mining and Technology(Grant No.2023WLKXJ182).
文摘Microwave-assisted rock-breaking technology,as a novel hybrid approach,is anticipated to facilitate the efficient excavation of complex rock formations.It is therefore crucial to understand the damage and failure mechanisms of rocks that have been subjected to irradiation.In this study,uniaxial compression experiments were conducted on granite specimens after 1.4 kW microwave irradiation for varying durations.Furthermore,a numerical method was proposed to solve electromagnetic-thermal-mechanical coupling problems by integrating finite and discrete elements.The results demonstrated a differential temperature distribution(high temperature in the middle and low-temperature areas at the ends)in the granite specimens under microwave irradiation,which resulted in a notable reduction in their physical and mechanical properties.As the duration of irradiation increased,the rate of heating and the extent of strength reduction both diminished,while the morphology and distribution of cracks at ultimate failure became increasingly complex.The numerical method effectively addresses the simulation challenges associated with the electromagnetic selective heating of granite containing multiple polar minerals under microwave irradiation.This approach accounted for the non-uniform thermal expansion of the minerals and provided a comprehensive model of damage progression under compression.
文摘Rockbursting in deep tunnelling is a complex phenomenon posing significant challenges both at the design and construction stages of an underground excavation within hard rock masses and under high in situ stresses. While local experience, field monitoring, and informed data-rich analysis are some of the tools commonly used to manage the hazards and the associated risks, advanced numerical techniques based on discontinuum modelling have also shown potential in assisting in the assessment of rockbursting. In this study, the hybrid finite-discrete element method(FDEM) is employed to investigate the failure and fracturing processes, and the mechanisms of energy storage and rapid release resulting in bursting, as well as to assess its utility as part of the design process of underground excavations.Following the calibration of the numerical model to simulate a deep excavation in a hard, massive rock mass, discrete fracture network(DFN) geometries are integrated into the model in order to examine the impact of rock structure on rockbursting under high in situ stresses. The obtained analysis results not only highlight the importance of explicitly simulating pre-existing joints within the model, as they affect the mobilised failure mechanisms and the intensity of strain bursting phenomena, but also show how the employed joint network geometry, the field stress conditions, and their interaction influence the extent and depth of the excavation induced damage. Furthermore, a rigorous analysis of the mass and velocity of the ejected rock blocks and comparison of the obtained data with well-established semi-empirical approaches demonstrate the potential of the method to provide realistic estimates of the kinetic energy released during bursting for determining the energy support demand.
基金the support of CSIRO and the Australia-Japan Foundation(Grant No.17/20470)supported by the Japan Society for the Promotion of Science KAKENHI(Grant No.JP18K14165)for Grant-in-Aid for Young Scientists。
文摘Mechanical cutting provides one of the most flexible and environmentally friendly excavation methods.It has attracted numerous efforts to model the rock chipping and fragmentation process,especially using the explicit finite element method(FEM) and bonded particle model(BPM),in order to improve cutting efficiency.This study investigates the application of a general-purpose graphic-processing-unit parallelised hybrid finite-discrete element method(FDEM) which enjoys the advantages of both explicit FEM and BPM,in modelling the rock chipping and fragmentation process in the rock scratch test of mechanical rock cutting.The input parameters of FDEM are determined through a calibration procedure of modelling conventional Brazilian tensile and uniaxial compressive tests of limestone,A series of scratch tests with various cutting velocities,cutter rake angles and cutting depths is then modelled using FDEM with calibrated input parameters.A few cycles of cutter/rock interactions,including their engagement and detachment process,are modelled for each case,which is conducted for the first time to the best knowledge of the authors,thanks to the general purpose graphic processing units(GPGPU) parallelisation.The failure mechanism,cutting force,chipping morphology and effect of various factors on them are discussed on the basis of the modelled results.Finally,it is concluded that GPGPU-parallelised FDEM provides a powerful tool to further study rock cutting and improve cutting efficiencies since it can explicitly capture different fracture mechanisms contributing to the rock chipping as well as chip formation and the separation process in mechanical cutting.Moreover,it is concluded that chipping is mostly owed to the mix-mode Ⅰ-Ⅱ fracture in all cases although mode Ⅱ cracks and mode Ⅰ cracks are the dominant failures in rock cutting with shallow and deep cutting depths,respectively.The chip morphology is found to be a function of cutter velocdty,cutting depth and cutter rake angle.
文摘Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.
基金National Key Research and Development Program of China under Grant Nos. 2018YFC1504400 and 2019YFC1509301Natural Science Foundation of China under Grant No. 52078471Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration under Grant No. 19EEEVL0402
文摘Due to the long construction life,improper design methods,brittle material properties and poor construction techniques,most existing masonry structures do not perform well during earthquakes.The retrofitting method using an external steel-meshed mortar layer is widely used to retrofit existing masonry buildings.Assessing the seismic performance of masonry walls reinforced by an external steel-meshed mortar layer reasonably and effectively is a difficult subject in the research field of masonry structures.Based on the combined finite-discrete elements method,the numerical models of retrofitted brick walls with four different masonry mortar strengths by an external mortar layer are established.The shear strength of mortar and the contact between the retrofitted mortar layer and the brick blocks are discussed in detail.The failure patterns and load-displacement curves of the retrofitted brick walls were obtained by applying low cycle reciprocating loads to the numerical model,and the bearing capacity and the failure mechanism of the retrofitted walls were obtained by comparing the failure patterns,ultimate bearing capacity,deformability and other aspects with the tests.This study provides a basis for improving the seismic strengthening design method of masonry structures and helps to better assess the seismic performance of masonry structures after retrofitting.
文摘This paper presents the application of a hybrid finite-discrete element method to study blast-induceddamage in circular tunnels. An extensive database of field tests of underground explosions above tunnelsis used for calibrating and validating the proposed numerical method; the numerical results areshown to be in good agreement with published data for large-scale physical experiments. The method isthen used to investigate the influence of rock strength properties on tunnel durability to withstand blastloads. The presented analysis considers blast damage in tunnels excavated through relatively weak（sandstone） and strong （granite） rock materials. It was found that higher rock strength will increase thetunnel resistance to the load on one hand, but decrease attenuation on the other hand. Thus, undercertain conditions, results for weak and strong rock masses are similar. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金supported by the Natural Sciences and Engineering Research Council of Canada through Discovery Grant 341275 (G. Grasselli) and Engage EGP 461019-13
文摘Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金supported by the State Key Laboratory of Intelligent Construction and Healthy Operation and Maintenance of Deep Underground Engineering(Grant No.SKLGDUEK2206)National Natural Science Foundation of China(Grant No.11872340).
文摘Soil desiccation cracking is a common phenomenon on the earth surface.Numerical modeling is an effective approach to study the desiccation cracking mechanism of soil.This work develops a novel 3D moisture diffusion discrete model that is capable of dynamically assessing the effect of cracking on moisture diffusion and allowing moisture to be discontinuous on both sides of the cracks.Then,the parametric analysis of the moisture exchange coefficient in the 3D moisture diffusion discrete model is carried out for moisture diffusion in continuous media,and the selection criterion of the moisture exchange coefficient for the unbroken cohesive element is given.Subsequently,an example of moisture migration in a medium with one crack is provided to illustrate the crack hindering effect on moisture migration.Finally,combining the 3D moisture diffusion discrete model with the finite-discrete element method(FDEM),the moisture diffusion-fracture coupling model is built to study the desiccation cracking in a strip soil and the crack pattern of a rectangular soil.The evolution of crack area and volume with moisture content is quantitatively analyzed.The modeling number and average width of cracks in the strip soil show a good consistency with the experimental results,and the crack pattern of the rectangular soil matches well with the existing numerical results,validating the coupled moisture diffusion-fracture model.Additionally,the parametric study of soil desiccation cracking is performed.The developed model offers a powerful tool for exploring soil desiccation cracking.
基金Supported by Major Instrument Project of National Natural Science Foundation of China(52327803)Major Project of National Natural Science Foundation of China(52192622).
文摘Based on the finite-discrete element method,a three-dimensional numerical model for axial impact rock breaking was established and validated.A computational method for energy conversion during impact rock breaking was proposed,and the effects of conical tooth forward rake angle,rock temperature,and impact velocity on rock breaking characteristics and energy transfer laws were analyzed.The results show that during single impact rock breaking with conical tooth bits,merely 7.52%to 12.51%of the energy is utilized for rock breaking,while a significant 57.26%to 78.10%is dissipated as frictional loss.An insufficient forward rake angle increases tooth penetration depth and frictional loss,whereas an excessive forward rake angle reduces penetration capability,causing bit rebound and greater energy absorption by the drill rod.Thus,an optimal forward rake angle exists.Regarding environmental factors,high temperatures significantly enhance impact-induced rock breaking.Thermal damage from high temperatures reduces rock strength and inhibits its energy absorption.Finally,higher impact velocities intensify rock damage,yet excessively high velocities increase frictional loss and reduce the proportion of energy absorbed by the rock,thereby failing to substantially improve rock breaking efficiency.An optimal impact velocity exists.
基金supported by the European Union's Horizon 2020 research and innovation programme under grant agreement No.101006752(ORCHYD project).
文摘For rapid and cost-effective hammer drilling,accurate prediction of rock impact response is crucial for designing optimal bits and maximising rock fragmentation.Current design optimisation workflows combine numerical simulations and experiments but often require numerous iterations to pinpoint the optimal design.Although physics-based models can potentially reduce experimental expenses,their significant computational demands present challenges when simulating the complex fragmentation dynamics during drill bit-rock interactions.This study introduces a data-driven artificial intelligence(AI)model,employing a multilayer perceptron(MLP)as a surrogate.The model leverages the hybrid finitediscrete element model(FDEM)as a powerful method in rock fracture mechanics to generate a sufficiently large training dataset.An automated workflow has been developed for generating the training data,comprising a pipeline that includes pre-processing,solving,and post-processing modules.Subsequently,the AI models were integrated into an optimisation framework alongside uncertainty quantification to demonstrate their potential in enhancing drilling efficiency through optimised bit design and operations.The MLP exhibits high accuracy in predicting key parameters,including rebound velocity,total crack length,quantities of fragments with different sizes and maximum contact force between rock and insert.Notably,this approach achieves real-time prediction compared to the 5-7 min simulation times of FDEM.Integrating this data-driven model into a design framework enables rapid assessment of different bit designs under various operational conditions.More broadly,this approach has the potential to impact other applications,such as digital twins,serving as a forward and inverse model for predicting rock type and optimising drilling performance.
基金financial support to this study from the National Natural Science Foundation of China,NSFC(Grant No.52278367).
文摘Particle morphology is critical in affecting the crushing behavior of rockfill materials.In contrast,most current single particle simulations lack satisfactory morphology accuracy,and the resulting crushing modes deviate from observations to some extent.Therefore,we reconstruct the real particle morphology with the spherical harmonic(SH)method and employ the finite-discrete element method(FDEM)to simulate the one-dimensional(1D)compressive crushing process of basalt particles commonly used in rockfill.The influences of four main morphological parameters,i.e.sphericity,aspect ratio,roundness,and convexity,on the single particle strength and the crushing modes are discussed.The results show that with the SH degree set to 15 and a mesh number of 20,480,the FDEM models of reconstructed particles achieve sufficient morphology accuracy and high computational efficiency.Based on the model,the simulation results demonstrate that the aspect ratio has the most significant impact on single particle strength,followed by sphericity.In contrast,roundness and convexity have a weaker effect than the above two parameters.Also,it is revealed that single particle strength decreases with increasing aspect ratio and sphericity,while it increases with higher roundness and convexity.Furthermore,aspect ratio significantly changes the initial crushing position,sphericity dominates post-crushing fragment size and quantity,and roundness mainly affects post-crushing morphology.The model results have been employed in establishing a support vector regression(SVR)-based predicted model,exhibiting good predictive performance and advantages for the optimization of rockfill particles in engineering.
基金supported by the Natural Sciences and Engineering Research Council of Canada(Grant Nos.Discovery 341275,and CRDPJ 543894-19)NSERC/Energi Simulation Industrial Research Chair programState Key Laboratory of Geohazard Prevention and Geoenvironment Protection Open Fund(Grant No.SKLGP2024K001).
文摘Frost heave in water-bearing rock masses poses significant threats to geotechnical engineering.This paper developed a novel three-dimensional(3D)frost model,based on the combined finite-discrete element method(FDEM),to investigate the frost heave process in rock masses where thermal transfer,water migration,water-ice phase transition(ice growth)and ice-rock interaction are explicitly simulated.The proposed model is first validated against existing experimental and analytical solutions,and further applied to investigate path-dependent frost heave behavior under various freezing conditions.Results show that freezing direction plays a vital role in the dynamic ice growth and ice-rock interaction,thus affecting the frost heave behavior.In the top-down freezing regime,ice plugs form first at the crack's top surface,sealing the crack and preventing water migration,which can amplify ice pressure.Parametric studies,including rock Young's modulus,ice-rock friction,and rock hydraulic conductivity,further reveal that the temporal aspects of ice development and rock mechanical response strongly affect ice-rock interaction and hence the frost heave mechanism.Furthermore,some typical phenomena(e.g.water/ice extrusion and frost cracking)can also be well captured in this model.This novel numerical framework sheds new light on frost heave behavior and enriches our understanding of frost heave mechanisms and ice-rock interaction processes within cold environment engineering projects.
基金support from the Natural Science Foundation of China(Grant Nos.41941018,U21A20153,42177140).
文摘The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.
基金supported by the Natural Sciences and Engineering Research Council of Canada (NSERC)Discovery Grants 341275,NSERC CRDPJ 543894-19,and NSERC/Energi Simulation Industrial Research Chair programfunding he received from Lassonde International Graduate Scholarship in Mining at the University of Toronto+1 种基金supported by the FCE Start-up Fund for New Recruits at the Hong Kong Polytechnic University (P0034042)the Early Career Scheme and the General Research Fund Scheme of the Research Grants Council of the Hong Kong SAR,China (Project Nos.PolyU 25220021 and PolyU 15227222).
文摘This paper presents the development of a coupled modeling approach to simulate cryogenic thermo-hydro-mechanical(THM)processes associated with a freezing medium,which is then implemented in the combined finite-discrete element method code(FDEM)for multi-physics simulation.The governing equations are deduced based on energy and mass conservation,and static equilibrium equations,considering water/ice phase change,where the strong couplings between multi-fields are supplemented by critical coupling parameters(e.g.unfrozen water content,permeability,and thermal conductivity).The proposed model is validated against laboratory and field experiments.Results show that the cryogenic THM model can well predict the evolution of strongly coupled processes observed in frozen media(e.g.heat transfer,water migration,and frost heave deformation),while also capturing,as emergent properties of the model,important phenomena(e.g.latent heat,cryogenic suction,ice expansion and distinct three-zone distribution)caused by water/ice phase change at laboratory and field scales,which are difficult to be all revealed by existing THM models.The novel modeling framework presents a gateway to further understanding and predicting the multi-physical coupling behavior of frozen media in cold regions.
基金supported by the Research Grants Council of Hong Kong (General Research Fund Project Nos.17200721 and 17202423)the National Natural Science Foundation of China (Grant No.42377149).
文摘A new thermomechanical(TM)coupled finite-discrete element method(FDEM)model,incorporating heat conduction,thermal cracking,and contact heat transfer,has been proposed for both continuous and discontinuous geomaterials.This model incorporates a heat conduction model that can accurately calculate the thermal field in continuousediscontinuous transition processes within a finite element framework.A modified contact heat transfer model is also included,which accounts for the entire contact area of discrete bodies.To align with the finite strain theory utilized in the FDEM mechanics module,the TM coupling module in the model is based on the multiplicative decomposition of the deformation gradient.The proposed model has been applied to various scenarios,including heat conduction in both continuous and discontinuous media during transient states,thermal-induced strain and stress,and thermal cracking conditions.The thermal field calculation model and the TM coupling model have been validated by comparing the numerical results with experiment findings and analytical solutions.These numerical cases demonstrate the reliability of the proposed model convincingly,making it suitable for use across a wide range of continuous and discontinuous media.
基金supported by the International Cooperation Programme of Chengdu City(No.2020-GH02-00023-HZ)。
文摘This paper introduces a deep learning workflow to predict phase distributions within complex geometries during two-phase capillary-dominated drainage.We utilize subsamples from Computerized Tomography(CT)images of rocks and incorporate pixel size,interfacial tension,contact angle,and pressure as inputs.First,an efficient morphology-based simulator creates a diverse dataset of phase distributions.Then,two commonly used convolutional and recurrent neural networks are explored and their deficiencies are highlighted,particularly in capturing phase connectivity.Subsequently,we develop a Higher-Dimensional Vision Transformer(HD-ViT)that drains pores solely based on their size,with phase connectivity enforced as a post-processing step.This enables inference for images of varying sizes,resolutions,and inlet-outlet setup.After training on a massive dataset of over 9.5 million instances,HD-ViT achieves excellent performance.We demonstrate the accuracy and speed advantage of the model on new and larger sandstone and carbonate images.We further evaluate HD-ViT against experimental fluid distribution images and the corresponding Lattice-Boltzmann simulations,producing similar outcomes in a matter of seconds.In the end,we train and validate a 3D version of the model.
基金the Natural Sciences and Engineering Research Council of Canadathe Ministry of National Defensethe RMC Green Team for providing the funding and the resources
文摘Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-existing joints in the damage evolution around the underground opening is of critical importance as they govern the fracturing mechanisms and influence the brittle responses of these hard rock masses under highly anisotropic in situ stresses.In this study,the main focus is the impact of joint network geometry,joint strength and applied field stresses on the rock mass behaviours and the evolution of excavation induced damage due to the loss of confinement as a tunnel face advances.Analysis of such a phenomenon was conducted using the finite-discrete element method(FDEM).The numerical model is initially calibrated in order to match the behaviour of the fracture-free,massive Lac du Bonnet granite during the excavation of the Underground Research Laboratory(URL)Test Tunnel,Canada.The influence of the pre-existing joints on the rock mass response during excavation is investigated by integrating discrete fracture networks(DFNs)of various characteristics into the numerical models under varying in situ stresses.The numerical results obtained highlight the significance of the pre-existing joints on the reduction of in situ rock mass strength and its capacity for extension with both factors controlling the brittle response of the material.Furthermore,the impact of spatial distribution of natural joints on the stability of an underground excavation is discussed,as well as the potentially minor influence of joint strength on the stress induced damage within joint systems of a non-persistent nature under specific conditions.Additionally,the in situ stress-joint network interaction is examined,revealing the complex fracturing mechanisms that may lead to uncontrolled fracture propagation that compromises the overall stability of an underground excavation.
文摘The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.
基金Supported by Marie Curie International Incoming Fellowship (No. PIIF-GA-2009-253453)
文摘A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.
基金supported by the Natural Science and Engineering Research Council (NSERC) of Canada in the form of discovery grant No. 341275the Swiss National Cooperative for the Disposal of Radioactive Waste (NAGRA)
文摘Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method （FDEM） code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone （EDZ） around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
