Predicting the probability distribution of Martian rocks mechanical property based on microscale rock mechanical experiments and accurate grain-based modeling

The exploration of Mars would heavily rely on Martian rocks mechanics and engineering technology.As the mechanical property of Martian rocks is uncertain,it is of utmost importance to predict the probability distribution of Martian rocks mechanical property for the success of Mars exploration.In this paper,a fast and accurate probability distribution method for predicting the macroscale elastic modulus of Martian rocks was proposed by integrating the microscale rock mechanical experiments(micro-RME),accurate grain-based modeling(AGBM)and upscaling methods based on reliability principles.Firstly,the microstructure of NWA12564 Martian sample and elastic modulus of each mineral were obtained by micro-RME with TESCAN integrated mineral analyzer(TIMA)and nanoindentation.The best probability distribution function of the minerals was determined by Kolmogorov-Smirnov(K-S)test.Secondly,based on best distribution function of each mineral,the Monte Carlo simulations(MCS)and upscaling methods were implemented to obtain the probability distribution of upscaled elastic modulus.Thirdly,the correlation between the upscaled elastic modulus and macroscale elastic modulus obtained by AGBM was established.The accurate probability distribution of the macroscale elastic modulus was obtained by this correlation relationship.The proposed method can predict the probability distribution of Martian rocks mechanical property with any size and shape samples.

Influence of grain size on strength of polymineralic crystalline rock: New insights from DEM grain-based modeling

Grain size effect on rock strength is a topic of great interest in geotechnical engineering.A consensus obtained from earlier laboratory tests is that rock strength generally decreases with the increase of grain size for both silicate and carbonate rocks;however,some recent numerical results conflict with such laboratory test results.To address this intriguing issue,the effect of grain size on strength of polymineralic crystalline rock with low porosity is investigated numerically using the grain-based modeling(GBM)approach in discrete element method(DEM)by interpreting micro-cracking process in response to loading.In agreement with some previous DEM simulation results,the simulated rock strength is found to increase with increasing grain size for both homogeneous and heterogeneous models,even when the number of assembled disks in one mineral grain changes.The mechanism of strength increase with increasing grain size is mainly associated with the number of assembled smooth-joint contacts along grain interfaces and the generation of grain boundary cracks in response to loading.The grain interfaces significantly weaken the integrity of the rock model,which is similar to effects of inherent defects in real rock.As the grain size increases,fewer grain interfaces are built in the model and the rock strength becomes much higher.Hence,by solely changing the mineral grain size in a model,the mechanism of grain size effect as observed in laboratory tests cannot be replicated.To address this issue,a method of degradation of grain boundary strength parameters is used to mimic the possible mechanism of grain size effect.The simulated strength using the method becomes comparable with those obtained from laboratory tests when the heterogeneity in the rock is considered.Degradation of grain boundary parameters with increasing grain size provides a plausible explanation for the grain size effect on rock strength.

3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing

A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Calibration of coupled hydro-mechanical properties of grain-based model for simulating fracture process and associated pore pressure evolution in excavation damage zone around deep tunnels

The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.

Deterioration of equivalent thermal conductivity of granite subjected to heating-cooling treatment

Understanding the thermal conductivity of granite is critical for many geological and deep engineering applications.The heated granite was subjected to air-,water-,and liquid nitrogen(LN2-)coolings in this context.The transient hot-wire technique was used to determine the equivalent thermal conductivity(ETC)of the granite before and after treatment.The deterioration mechanism of ETC is analyzed from the meso-perspective.Finally,the numerical model is used to quantitatively study the impact of cooling rate on the microcrack propagation and heat conduction characteristics of granite.The results show that the ETC of granite is not only related to the heating temperature,but also affected by the cooling rate.The ETC of granite decreases nonlinearly with increasing heating temperature.A faster cooling rate causes a greater decrease in ETC at the same heating temperature.The higher the heating temperature,the stronger the influence of cooling rate on ETC.The main explanation for the decrease in ETC of granite is the increase in porosity and microcrack density produced by the formation and propagation of pore structure and microcracks during heating and cooling.Further analysis displays that the damage of granite at the heating stage is induced by the difference in thermal expansion and elastic properties of mineral particles.At the cooling stage,the faster cooling rate causes a higher temperature gradient,which in turn produces greater thermal stress.As a result,it not only causes new cracks in the granite,but also aggravates the damage at the heating stage,which induces a further decrease in the heat conduction performance of granite,and this scenario is more obvious at higher temperatures.

Effect of the mineral spatial distribution heterogeneity on the tensile strength of granite:Insights from PFC3D-GBM numerical analysis

The mechanical characteristics of crystalline rocks are affected by the heterogeneity of the spatial distribution of minerals.In this paper,a novel three-dimensional(3D)grain-based model(GBM)based on particle flow code(PFC),i.e.PFC3D-GBM,is proposed.This model can accomplish the grouping of mineral grains at the 3D scale and then filling them.Then,the effect of the position distribution,geometric size,and volume composite of mineral grains on the cracking behaviour and macroscopic properties of granite are examined by conducting Brazilian splitting tests.The numerical results show that when an external load is applied to a sample,force chains will form around each contact,and the orientation distribution of the force chains is uniform,which is independent of the external load level.Furthermore,the number of high-strength force chains is proportional to the external load level,and the main orientation distribution is consistent with the external loading direction.The main orientation of the cracks is vertical to that of the high-strength force chains.The geometric size of the mineral grains controls the mechanical behaviours.As the average grain size increases,the number of transgranular contacts with higher bonding strength in the region connecting both loading points increases.The number of high-strength force chains increases,leading to an increase in the stress concentration value required for the macroscopic failure of the sample.Due to the highest bonding strength,the generation of transgranular cracks in quartz requires a higher concentrated stress value.With increasing volume composition of quartz,the number of transgranular cracks in quartz distributed in the region connecting both loading points increases,which requires many high-strength force chains.The load level rises,leading to an increase in the tensile strength of the numerical sample.

Influence of inter-grain cementation stiffness on the effective elastic properties of porous Bentheim sandstone

Effective elastic properties of porous media are known to be significantly influenced by porosity.In this paper,we investigated the influence of another critical factor,the inter-grain cementation stiffness,on the effective elastic properties of a granular porous rock(Bentheim sandstone)using an advanced numerical workflow with realistic rock microstructure and a theoretical model.First,the disparity between the experimentally tested elastic properties of Bentheim sandstone and the effective elastic properties predicted by empirical equations was analysed.Then,a micro-computed tomography(CT)-scan based approach was implemented with digital imaging software AVIZO to construct the 3D(three-dimensional)realistic microstructure of Bentheim sandstone.The microstructural model was imported to a mechanics solver based on the 3D finite element model with inter-grain boundaries modelled by cohesive elements.Loading simulations were run to test the effective elastic properties for different shear and normal intergrain cementation stiffness.Finally,a relation between the macroscale Young’s modulus and inter-grain cementation stiffness was derived with a theoretical model which can also account for porosity explicitly.Both the numerical and theoretical results indicate the influence of the inter-grain cementation stiffness,on the effective elastic properties is significant for porous sandstone.The calibrated normal and shear stiffnesses at the inter-grain boundaries are 1.2×10^(5) and 4×10^(4) GPa/m,respectively.

Numerical Study on the Mechanical Behavior of Hydraulic Fractures under Stress Disturbance,based on a Multi Pb-GBM Method

Water inrush is a common disaster in submarine mining.The key to preventing this disaster is to restrict the water-conducting fissure propagation induced by mining from forming a flow channel to communicate with the overlying aquifer.The mechanical behavior of hydraulic fractures under stress disturbance lies at the heart of the problem.Hence,the multiple parallel bond-grain-based model(multi Pb-GBM)is introduced to explore the hydraulic fracture evolution law of crystalline granite under the influence of stress disturbance.The results show that:hydraulic fracturing under stress disturbance is clearly affected by the stress wave frequency;the higher the frequency is,the more difficult it is to crack,but the crack propagation speed is faster after crack initiation;the propagation direction of a crack is deflected towards the propagation direction of the stress wave and the crack dip angle is controlled by the maximum principal stress;the internal crystal boundary of the same mineral is the most stable one among the three contacts,the contact boundary between different minerals being the most fragile one.This research not only has a practical application value to seabed mining engineering,but also has important theoretical significance in enriching deep rock mechanics theory.