Search results

Accurate presentation of the rock microstructure is critical to the grain-scale analysis of rock deformation and failure in numerical modelling. 3D granite microstructure modelling has only been used in limited studies with the mineral pattern often remaining poorly constructed. In this study, the authors developed a new approach for generating 2D and 3D granite microstructure models from a 2D image by combining a heterogeneous material reconstruction method (simulated annealing method) with Voronoi tessellation.

More specifically, the stochastic information in the 2D image is first extracted using the two-point correlation function (TPCF). Then an initial 2D or 3D Voronoi diagram with a random distribution of the minerals is generated and optimised using a simulated annealing method until the corresponding TPCF is consistent with that in the 2D image. The generated microstructure model accurately inherits the stochastic information (e.g. volume fraction and mineral pattern) from the 2D image. Lastly, the authors compared the topological characteristics and mechanical properties of the 2D and 3D reconstructed microstructure models with the model obtained by direct mapping from the 2D image of a real rock sample.

The good agreements between the mapped and reconstructed models indicate the accuracy of the reconstructed microstructure models on topological characteristics and mechanical properties.

The newly developed reconstruction method successfully transfers the mineral pattern from a granite sample into the 2D and 3D Voronoi-based microstructure models ready for use in grain-scale modelling.

Evaluating mechanical properties of simply made samples by 3D printing technology at nanoscale provides a clear path to better understand larger-scale responses of complex natural rocks. Therefore, to realize the similarity between synthetically manufactured materials and natural geomaterials, this study focused on nanoscale mechanical characterization of a 3D printed object with only two constituent components (gypsum powder and infiltrant).

The study method includes nanoindentation technique combined with numerical simulation via discrete element method (DEM).

Four typical load-displacement curves were identified from nanoindentation of total test points indicating a typical elastic-plastic behavior of the 3D printed gypsum rock sample. Mechanical parameters such as Young’s modulus and hardness were calculated by energy-based methods and a positive correlation was observed. The infiltrant was found to considerably be responsible for the majority of the sample nano-mechanical behavior rather than the gypsum particles, thus expected to control macroscale properties. This was decided from deconvolution and clustering of elastic modulus data. Particle flow modeling in DEM was used to simulate the nanoindentation process in a porous media yielding rock-alike mechanical behavior.

The results show a matching load-displacement response between experimental and simulation results, which verified the credibility of simulation modeling for mechanical behavior of 3D printed gypsum rock at nanoscale. Finally, differential effective medium theory was used to upscale the nanoindentation results to the macroscale mechanical properties, which provided an insight into the geomechanical modeling at multiscale.

Discrete element methods are numerical procedures for simulating the complete behaviour of systems of discrete, interacting bodies. Three important aspects of discrete element programs are examined: (1) the representation of contacts; (2) the representation of solid material; and (3) the scheme used to detect and revise the set of contacts. A proposal is made to define what constitutes a discrete element program, and four classes of such programs are described: the distinct element method, modal methods, discontinuous deformation analysis and the momentum‐exchange method. Several applications and examples are presented, and a list is given of suggestions for future developments.

The correlations between mechanical behaviour, tensile strength, and rock parameters of metasedimentary rock samples in Karak, Pahang’s New Austrian Tunnelling Method (NATM) were statistically evaluated from the rock mechanic laboratory works at the selected sections around 2,000 m of the tunnel (named as NATM-1). According to a statistical analysis, lithotypes, geological structures, and region geology have a significant impact on the mechanical behaviour of the metasedimentary rock. In the Brazilian test, the fracture behaviour of the disc specimens was highly related to the reliability and precision of the experimental data by validations of methods. In this work, the impact of different loading methods and rock lithotypes on the failure mechanism of Brazilian discs was examined utilising five different metasedimentary rock types and three different loading methods. During the loading operation, the strain and displacement fields of the specimens were recorded and evaluated using a computerised strain gauge system. The rock types, according to experimental data, have a significant impact on the peak load and deformation properties of Brazilian discs. With the method below, tensile strength point of a disc specimen is clearly regulated by the material stiffness and tensile–compression ratio. Seismic occurrences have had a substantial impact on changing the rock and exerting forces that may affect its mechanical characteristics as well as its vulnerability to weathering effects or discontinuities. As a result, the goal of this study is to look into the connection between rock mechanics and metasedimentary rock stress analysis in NATM-1, Karak, Pahang.

This paper aims to present a digital image processing (DIP)-based discrete element method (DEM) for the analysis of heterogeneous geomaterials. Taking a soil and rock mixture as an example, the direct shear test is used to illustrate the application of this method. The numerical result is validated by the laboratory experiment and implies its feasibility in the analysis of heterogeneous geomaterials.

This method has two major steps. Based on a modification of the connected-component labeling algorithm, a novel vectorization method, which can transform the digital photos to vectorized geometry automatically, is proposed first. Then, a simple yet effective method for the generation of heterogeneous DEM models is presented using the simulation of simplicity technique.

DIP-DEM method is a feasible approach for the analysis of mechanical behavior of heterogeneous material. For soil and rock mixtures (SRM), the horizantal deformation at peak shear point becomes larger with the normal stress. Compared with pure soil, the rock aggregates mainly improve the friction angle of SRM.

As a universal method taking advantage of both DIP and DEM, this method has broad application prospects in related fields.

Hydrofracturing technology has been widely used in tight oil and gas reservoir exploitation, and the fracture network formed by fracturing is crucial to determining the resources recovery rate. Due to the complexity of fracture network induced by the random morphology and type of fluid-driven fractures, controlling and optimising its mechanisms is challenging. This paper aims to study the types of multiscale mode I/II fractures, the fluid-driven propagation of multiscale tensile and shear fractures need to be studied.

A dual bilinear cohesive zone model (CZM) based on energy evolution was introduced to detect the initiation and propagation of fluid-driven tensile and shear fractures. The model overcomes the limitations of classical linear fracture mechanics, such as the stress singularity at the fracture tip, and considers the important role of fracture surface behaviour in the shear activation. The bilinear cohesive criterion based on the energy evolution criterion can reflect the formation mechanism of complex fracture networks objectively and accurately. Considering the hydro-mechanical (HM) coupling and leak-off effects, the combined finite element-discrete element-finite volume approach was introduced and implemented successfully, and the results showed that the models considering HM coupling and leak-off effects could form a more complex fracture network. The multiscale (laboratory- and engineering-scale) Mode I/II fractures can be simulated in hydrofracturing process.

Based on the proposed method, the accuracy and applicability of the algorithm were verified by comparing the analytical solution of KGD and PKN models. The effects of different in situ stresses and flow rates on the dynamic propagation of hydraulic fractures at laboratory and engineering scales were investigated. when the ratio of in situ stress is small, the fracture propagation direction is not affected, and the fracture morphology is a cross-type fracture. When the ratio of in situ stress is relatively large, the propagation direction of the fracture is affected by the maximum in situ stress, and it is more inclined to propagate along the direction of the maximum in situ stress, forming double wing-type fractures. Hydrofracturing tensile and shear fractures were identified, and the distribution and number of each type were obtained. There are fewer hydraulic shear fractures than tensile fractures, and shear fractures appear in the initial stage of fracture propagation and then propagate and distribute around the perforation.

The proposed dual bilinear CZM is effective for simulating the types of Mode I/II fractures and seizing the fluid-driven propagation of multiscale tensile and shear fractures. Practical fracturing process involves the multi-type and multiscale fluid-driven fracture propagation. This study introduces general fluid-driven fracture propagation, which can be extended to the fracture propagation analysis of potential fluid fracturing, such as other liquids or supercritical gases.

To study fracture characteristics of jointed rock masses under blasting load, the RFPA2D analysis software for dynamic fracture of rocks based on the finite element method and statistical damage theory was used.

On this basis, this research simulated the fracture process of rock masses in blasting with different joint geometrical characteristics and mainly analysed the influences of distance from joints to blasting holes, the length of joints, the number of joints and joint angle on fracture of rock masses.

The calculation results show that with the constant increase of the distance from joints to blasting holes, the influences of joints on blasting effects of rock masses gradually reduced. Rock masses with long joints experienced more serious damages than those with short joints. Damages obviously increased with the changing from rock masses without joints to rock masses with joints, and when there were three joints, the further increase of the number of joints had unobvious changes on blasting effects of rock masses. Joints showed significant guidance effect on the propagation of cracks in blasting: promoting propagation of main vertical cracks deflecting to the ends of joints.

The research results are expected to provide some theoretical bases in practical application of engineering blasting.

The purpose of the paper was the application of computational fluid dynamics (CFD) techniques in fluid flow using Maxwell’s equation for partial slip modelling, estimating the flow parameters, and selecting tangential momentum accommodation coefficient (TMAC) for tight rock samples in permeability calculations.

The paper presents a numerical analysis of fluid flow in a low-porosity rock sample by using CFD. Modelling results allowed to determine mass flow rates in a rock sample and to calculate permeability values using a modified Darcy’s equation. Three-dimensional (3D) geometrical model of rock sample generated using computed X-ray tomography was used in the analysis. Steady-state calculations were carried out for defined boundary conditions in the form of pressure drop. The simulations were applied taking into account the slip phenomenon described by Maxwell’s slip model and TMAC.

Values of permeability were calculated for different values of TMAC, which vary from 0 to 1. Results in the form of gas mass flow rates were compared with the measured value of permeability for rock sample, which confirmed the high accuracy of the presented model.

Calculations of fluid flow in porous media using CFD can be used to determine rock samples’ permeability. In slip flow regime, Maxwell’s slip model can be applied and the empirical value of TMAC can be properly estimated.

This paper presents the usage of CFD, Maxwell’s equation for partial slip modelling, in fluid flow mechanism for tight rock samples. 3D geometric models were generated using created pre-processor (poROSE software) and applied in the raw form for simulation.

Continua and discontinua coexist in natural rock materials. This paper aims to present an improved approach for addressing the mechanical response of rock masses based on the combined finite-discrete element method (FDEM) proposed by Munjiza.

Several algorithms have been programmed in the new approach. The algorithms include (1) a simpler and more efficient algorithm to calculate the contact force; (2) An algorithm for tangential contact force closer to the actual physical process; (3) a plastic yielding criterion (e.g. Mohr-Coulomb) to modify the elastic stress for fitting the mechanical behavior of elastoplastic materials; and (4) a complete code for the mechanical calculation to be implemented in Matrix Laboratory (MATLAB).

Three case studies, including two standard laboratory experiments (uniaxial compression and Brazilian split test) and one engineering-scale anti-dip slop model, are presented to illustrate the feasibility of the Y-Mat code and its ability to deal with multi-scale rock mechanics problems. The results, including the progressive failure process, failure mode and trajectory of each case, are acceptable compared to other corresponding studies. It is shown that, the code is capable of modeling geotechnical and geological engineering problems.

This article gives an improved FDEM-based numerical calculation code. And, feasibility of the code is verified through three cases. It can effectively solve the geotechnical and geological engineering problems.

The purpose of this study is to evaluate how accurately a 3D printer could manufacture basic porous models. Geoscience research is evolving toward numerical prediction of porous rock properties, but laboratory tests are still considered a standard practice. 3D printing digital designs of porous models (proxies) is a way to bridge the gap between these two realms of inquiry.

Digital designs of simple porous models have been 3D-printed on an inkjet-style (polyjet) 3D printer. Porosity and pore-throat size distribution of proxies have been measured with helium porosimetry, mercury porosimetry and computed tomography (CT) image analysis. Laboratory results on proxies have been compared with properties calculated on digital designs and CT images.

Bulk volume of proxies was by 0.6-6.7 per cent lower than digital designs. 3D-printed porosity increased from 0.2 to 1.9 per cent compared to digital designs (0-1.3 per cent). 3D-printed pore throats were thinner than designed by 10-31 per cent.

Incomplete removal of support material from pores yielded inaccurate property measurements. The external envelope of proxies has been 3D-printed at higher accuracy than pores.

Characterization of these simple models improves understanding of how more complex rock models can be 3D-printed accurately and how both destructive (mercury porosimetry) and non-destructive (CT and helium porosimetry) methods can be used to characterize porous models.

Validation of 3D-printed porous models using a suite of destructive and non-destructive methods is novel.