Numerical Research on Energy Evolution in Granite under Different Confining Pressures Using Otsu ’ s Digital Image Processing and PFC 2 D

Research on energy accumulation and releasing in the rock plays a key role on revealing its failure mechanism. This paper establishes a microscopic structure model of granite using Otsu digital image processing (DIP) technology and particle flow code software (PFC2D). A series of numerical compression tests under different confining pressures were conducted to investigate the macro and micro characteristics of energy evolution in granite. The results showed that the energy evolution of granite is divided into three stages: stable accumulation, slow dissipation, and rapid release. With increasing confining pressure, the strain energy accumulation ratio decreased exponentially and the peak value of strain energy increased linearly. It was found that the energy accumulation speed in the pre-peak stage increased as a linear function, while the energy release speed in the post-peak stage decreased as an exponential function. In addition, the feldspar is the main microstructure which played a major part in accumulating energy in granite. However, the unit mineral energy of mica particles was bigger than that of feldspar and quartz. When subjected to increasing confining pressure, the feldspar’s total energy growth rate was fastest. Meanwhile, the mica’s unit energy growth rate was fastest.


Introduction
Understanding the energy evolution in rock is critical for rock engineering design and assessment and has been one of the key problems facing modern researchers [1][2][3][4][5].The energy evolution in rock is determined by-and is therefore an indicator of-the deformation and failure of the rock [6][7][8][9].In recent years, many scholars have researched the energy accumulation, dissipation, and release in the rock by way of theoretical analysis and laboratory tests.Li et al. [10] researched the energy dissipation and release in rocks during triaxial compression with different loading and unloading paths.Yang et al. [11] studied the blasting response of a jointed rock mass and found that the theoretical criterion based on the minimum strain energy density factor predicted the crack initiation behavior well.Zhang et al. [12] developed a novel analytical solution to calculate the amount of released energy in coal and studied the influence of the coal's joint on energy release.Ning et al. [13] introduced a new energy dissipation method to identify crack initiation, and found that the variation in the pre-peak dissipation ratio was consistent with coal damage.Chen et al. [14] acquired similar results on shale samples.The abovementioned studies mainly focus on the macro characteristics of energy evolution in rocks [15,16].However, a rock is an aggregate of one or more minerals.The different minerals in a rock have significantly different mechanical properties.The complex microstructures thus play an important role in the macroscopic energy evolution and failure of rock.Therefore, it is significant to investigate the energy evolution of the microstructures in rock.The techniques used in existing studies on the micro-mechanical properties of rocks have mainly included the characterization of the heterogeneous distribution [17][18][19], computed tomography (CT) scanning techniques [20][21][22], and digital imaging processing (DIP) technology [23][24][25].DIP usually adopts an image segmentation approach to acquire the image information of microstructures in rock, and the information is then inputted into numerical software to realize the reconstruction of the microstructures.This technique is expected to be widely applied in the research of rock mechanics due to its merits including distinct theory and great applicability to instruments [26].
In this paper, the image information of microstructures in granite was acquired using Otsu's DIP technique.The corresponding microstructural model of granite was reconstructed using particle flow code software (PFC2D).A series of numerical compression tests were carried out under different confining pressures.The macro and micro characteristics of energy evolution in granite were then analyzed in detail.Additionally, the influence of confining pressure on energy accumulation and dissipation was explored.

Acquisition and Characterization of Microstructures in Rock
Rocks are heterogeneous materials composed of a variety of different mineral particles.The same minerals have similar colors, and the difference in color between different minerals is very large.On this basis, color of the mineral surface of rocks is endowed with a corresponding value, and a digital image of the rock comprises pixels with different color values.The color value of each pixel is associated with mineral components.Figure 1 illustrates the change in gray values of the pixels along the measured line on the surface of granite specimen.The abovementioned studies mainly focus on the macro characteristics of energy evolution in rocks [15,16].However, a rock is an aggregate of one or more minerals.The different minerals in a rock have significantly different mechanical properties.The complex microstructures thus play an important role in the macroscopic energy evolution and failure of rock.Therefore, it is significant to investigate the energy evolution of the microstructures in rock.The techniques used in existing studies on the micro-mechanical properties of rocks have mainly included the characterization of the heterogeneous distribution [17][18][19], computed tomography (CT) scanning techniques [20][21][22], and digital imaging processing (DIP) technology [23][24][25].DIP usually adopts an image segmentation approach to acquire the image information of microstructures in rock, and the information is then inputted into numerical software to realize the reconstruction of the microstructures.This technique is expected to be widely applied in the research of rock mechanics due to its merits including distinct theory and great applicability to instruments [26].
In this paper, the image information of microstructures in granite was acquired using Otsu's DIP technique.The corresponding microstructural model of granite was reconstructed using particle flow code software (PFC2D).A series of numerical compression tests were carried out under different confining pressures.The macro and micro characteristics of energy evolution in granite were then analyzed in detail.Additionally, the influence of confining pressure on energy accumulation and dissipation was explored.

Acquisition and Characterization of Microstructures in Rock
Rocks are heterogeneous materials composed of a variety of different mineral particles.The same minerals have similar colors, and the difference in color between different minerals is very large.On this basis, color of the mineral surface of rocks is endowed with a corresponding value, and a digital image of the rock comprises pixels with different color values.The color value of each pixel is associated with mineral components.Figure 1  The gray value was set in the range of 0-255.As shown in Figure 1, the change in gray value corresponds to the mineral composition of the granite.Mica is black with the lowest gray value, feldspar is white with the highest gray value, and quartz is transparent with a medium gray value in DIP.Based on the color of the minerals, the digital image of the rock was segmented and reconstructed to get the characterization image of the microstructures in the rock.
In order to segment the minerals of rock, the gray value range of each mineral needs to be determined.It was assumed that the gray value ranges of mica, quartz, and feldspar in granite were 0-t 1 , t 1 -t 2 , and t 2 -255, respectively.The parameters t 1 and t 2 are the segmentation thresholds of the minerals [27].Otsu's image segmentation method is a conventional image segmentation technique and is often considered for selecting the optimal threshold for image segmentation [28].Otsu's method was adopted to segment the minerals contained in the granite image.The threshold (t 1 , t 2 ) which meets the maximum gray variances σ 2 of all minerals is regarded as an optimal segmentation threshold.The maximum gray variance is as follows: where u is the overall gray value of the granite image, N i is the number of pixels of the different minerals, and u i is the average gray value of the different minerals.
Using Otsu's method, the processed image of the minerals in granite was obtained as shown in Figure 1.The distribution and shapes of minerals in the processed image were well consistent with those in the real image of the granite, and the content of feldspar, quartz and mica in this granite are 66%, 26.1% and 7.9%, respectively.

Microstructural Model Reconstruction Using PFC2D
After obtaining the processed image of microstructures in the granite, a corresponding microstructural model could be established by importing the image information into numerical simulation software.In this paper, using particle flow code software (PFC2D), the microstructural model of granite was reconstructed according to the following steps: (1) Exporting the coordinates and gray values of the pixels.The pixel size of an image is always different from the actual size.Before exporting, there should be a scale conversion.It is supposed that the image's actual size is a × b, and the pixel number is M × N.Each pixel's length l and coordinate (x i , y j ) can be calculated by Equation (2): where s is the scale factor and equals a/M or b/N, accordingly.(2) Importing data into the PFC2D.In this paper, uniaxial compression tests were carried out and the dimension of the sample model was 50 mm × 100 mm.The diameter of the sample's particles ranged from 0.25 to 0.75 mm, and the average value was 0.5 mm which equaled the converted pixel's length l in the processed image (Figure 1).Then, the coordinates and gray values of the pixels were imported into the particle flow model.In the PFC2D model, particles were distributed randomly.A Fish function (a built-in programming langue in PFC2D) was written to find the pixel which was closest to the particle.The particle's color code was set equal to the pixel's gray value.Figure 2 shows an example of incorporating a gray image into the particle flow model.(3) Assigning microscopic parameters to the microstructures.Because the mechanical properties of minerals in rock are different, the respective microscopic parameters should be assigned to the different minerals.Based on the aforementioned particle color codes, the different minerals could be identified and assigned the appropriate parameters.PFC2D usually utilizes a parallel bonded-particle model to simulate rocks [29,30].The parameters, including the friction coefficient μ, contact stiffness kn and ks, parallel bonded stiffness n k and s k , parallel bonded strength σn and σs, and parallel bonded radius λ, need to be supplied as input to the PFC2D.Because these parameters cannot be obtained by laboratory tests directly, a trial-and-error method was used to calibrate the microscopic parameters [31], and the simulated results were compared with those of the laboratory test.The microscopic parameters of the granite in the simulation are displayed in Table 1. Figure 3 shows a comparison of laboratory uniaxial compression test data and results simulated using PFC2D.As shown in Figure 3, the two stress-strain curves agreed well with each other and the failure modes of granite were both split fractures, indicating that the microscopic parameters were selected properly.After calibrating the parameters of the model, the established microstructural model of granite was used to carry out a series of compression tests under different confining pressures to further analyze the macro and micro characteristics of energy evolution in granite.In tests, the lateral and bottom loading plates were used to apply constant confining pressure on the specimens based on servo commands.The confining pressure was set to 0, 10, 20, 30, or 40 MPa.After the confining pressure reached the designated value, axial loading was applied by moving the lower loading plate PFC2D usually utilizes a parallel bonded-particle model to simulate rocks [29,30].The parameters, including the friction coefficient µ, contact stiffness k n and k s , parallel bonded stiffness k n and k s , parallel bonded strength σ n and σ s , and parallel bonded radius λ, need to be supplied as input to the PFC2D.Because these parameters cannot be obtained by laboratory tests directly, a trial-and-error method was used to calibrate the microscopic parameters [31], and the simulated results were compared with those of the laboratory test.The microscopic parameters of the granite in the simulation are displayed in Table 1. Figure 3 shows a comparison of laboratory uniaxial compression test data and results simulated using PFC2D.As shown in Figure 3, the two stress-strain curves agreed well with each other and the failure modes of granite were both split fractures, indicating that the microscopic parameters were selected properly.PFC2D usually utilizes a parallel bonded-particle model to simulate rocks [29,30].The parameters, including the friction coefficient μ, contact stiffness kn and ks, parallel bonded stiffness n k and s k , parallel bonded strength σn and σs, and parallel bonded radius λ, need to be supplied as input to the PFC2D.Because these parameters cannot be obtained by laboratory tests directly, a trial-and-error method was used to calibrate the microscopic parameters [31], and the simulated results were compared with those of the laboratory test.The microscopic parameters of the granite in the simulation are displayed in Table 1. Figure 3 shows a comparison of laboratory uniaxial compression test data and results simulated using PFC2D.As shown in Figure 3, the two stress-strain curves agreed well with each other and the failure modes of granite were both split fractures, indicating that the microscopic parameters were selected properly.After calibrating the parameters of the model, the established microstructural model of granite was used to carry out a series of compression tests under different confining pressures to further analyze the macro and micro characteristics of energy evolution in granite.In tests, the lateral and bottom loading plates were used to apply constant confining pressure on the specimens based on servo commands.The confining pressure was set to 0, 10, 20, 30, or 40 MPa.After the confining pressure reached the designated value, axial loading was applied by moving the lower loading plate reached the designated value, axial loading was applied by moving the lower loading plate until specimen failure.The loading mode was displacement control, and the loading speed was 0.01 mm/s.

The Energy Calculation Method
During the process of the deformation and failure of the granite, the work done by the external load is converted into two parts: one is the elastic strain energy accumulated in the granite, and the other is dissipation energy resulting from friction and movement between mineral grains.The energy can be written as: where E tot is the external input energy, E str is the strain energy, and E dis is the dissipation energy.
In the test, the external input energy was equal the work done by the loading plates' movement on the specimen.The lateral loading plate moved in the negative direction and generated negative work which was regarded as a part of the dissipation energy.Therefore, the external input energy was only the work done by the bottom loading plate, expressed as: where F ε and ∆U ε are the force and displacement increment of the bottom loading plate, respectively, when the specimen strain is ε.
In PFC2D simulation, the strain energy accumulated in the specimen includes two parts [32].One is the contact strain energy E c str stored at all contacts, and the other is the parallel-bond strain energy E pb str stored in parallel bonds, as shown in Figure 4.The strain energy can be expressed as: where In Equations ( 6) and ( 7), F n i , F s i , and M s i are the normal force, shear force, and the moment in the parallel bond i, respectively; F n i and F s i are the normal force and shear force in the contact i, respectively; A i and I i are the area and inertia moment of the bond cross section, respectively; N c is the number of contacts; and N pb is the number of parallel bonds.
until specimen failure.The loading mode was displacement control, and the loading speed was 0.01 mm/s.

The Energy Calculation Method
During the process of the deformation and failure of the granite, the work done by the external load is converted into two parts: one is the elastic strain energy accumulated in the granite, and the other is dissipation energy resulting from friction and movement between mineral grains.The energy can be written as: where Etot is the external input energy, Estr is the strain energy, and Edis is the dissipation energy.
In the test, the external input energy was equal the work done by the loading plates' movement on the specimen.The lateral loading plate moved in the negative direction and generated negative work which was regarded as a part of the dissipation energy.Therefore, the external input energy was only the work done by the bottom loading plate, expressed as: where Fε and ΔUε are the force and displacement increment of the bottom loading plate, respectively, when the specimen strain is ε.
In PFC2D simulation, the strain energy accumulated in the specimen includes two parts [32].One is the contact strain energy c str E stored at all contacts, and the other is the parallel-bond strain energy pb str E stored in parallel bonds, as shown in Figure 4.The strain energy can be expressed as: Where In Equations ( 6) and ( 7   In order to clearly and reasonably describe the process of energy accumulation and release, the following parameters need to be defined as follows: β energy accumulation ratio-equals the ratio of strain energy to total external input energy; v a energy accumulation speed-represents the increasing speed of strain energy in the pre-peak stage; v r energy release speed-represents the decreasing speed of strain energy in the post-peak stage; E s-t mineral energy-represents the total strain energy accumulated in one mineral; E s-u unit mineral energy-represents E s-t divided by the number of mineral particles.

Energy Accumulation and Release
Following the calculation method of energy in the granite model, the external input energy, strain energy, dissipation energy, and energy accumulation ratio β were recorded by compiling the FISH language in PFC2D.Figure 5 shows the energy evolution curves of granite under different confining pressures.As shown in Figure 5, the energy curves for granite could be divided into three stages: stable energy accumulation (Stage I), slow energy dissipation (Stage II), and rapid energy release (Stage III).
In order to clearly and reasonably describe the process of energy accumulation and release, the following parameters need to be defined as follows: β energy accumulation ratio-equals the ratio of strain energy to total external input energy; va energy accumulation speed-represents the increasing speed of strain energy in the pre-peak stage; vr energy release speed-represents the decreasing speed of strain energy in the post-peak stage; Es-t mineral energy-represents the total strain energy accumulated in one mineral; Es-u unit mineral energy-represents Es-t divided by the number of mineral particles.

Energy Accumulation and Release
Following the calculation method of energy in the granite model, the external input energy, strain energy, dissipation energy, and energy accumulation ratio β were recorded by compiling the FISH language in PFC2D.Figure 5 shows the energy evolution curves of granite under different confining pressures.As shown in Figure 5, the energy curves for granite could be divided into three stages: stable energy accumulation (Stage I), slow energy dissipation (Stage II), and rapid energy release (Stage III).In Stage I, due to few cracks having been generated in the granite specimen, only a small amount of external input energy could be dissipated in the form of friction and dislocation movement among granite particles.Most of the external input energy was converted into strain energy stored in the specimen, and the energy accumulation ratio β reached 80% in the initial loading stage.With the continual input of external energy, the energy accumulation ratio β gradually increased and tended to be stable.In Stage II, with increasing load on the specimen, the In Stage I, due to few cracks having been generated in the granite specimen, only a small amount of external input energy could be dissipated in the form of friction and dislocation movement among granite particles.Most of the external input energy was converted into strain energy stored in the specimen, and the energy accumulation ratio β reached 80% in the initial loading stage.With the continual input of external energy, the energy accumulation ratio β gradually increased and tended to be stable.In Stage II, with increasing load on the specimen, the dissipation energy increased gradually and the energy accumulation ratio β decreased distinctly.When the accumulated strain energy reaches the extreme value of energy storage, the excessive energy can accelerate the crack propagation in the specimen until it enters the rapid energy release stage.In Stage III, the accumulated strain energy was released rapidly and the energy accumulation ratio β reduced drastically as a result of many new cracks having been generated in the specimen.Finally, the specimen failed and lost its bearing capacity, but there was still a small amount of residual strain energy because of the confining pressure effect.

Effect of Confining Pressure
To accurately analyze the influence of confining pressure on the energy accumulation and release in the granite specimen, comparisons of the strain energy and accumulation ratio under different confining pressures are drawn in Figures 6 and 7, respectively.With increasing confining pressure, the maximum value of the strain energy accumulation ratio β decreased gradually.Conversely, the peak value of the accumulated strain energy in the specimen increased linearly [33,34].As demonstrated in Figure 8a, the maximum β decreased according to the exponential function β = 8.485e −0.069σ3 + 90.05, and the peak strain energy increased according to the linear function E s = 0.123σ 3 + 1.318.dissipation energy increased gradually and the energy accumulation ratio β decreased distinctly.When the accumulated strain energy reaches the extreme value of energy storage, the excessive energy can accelerate the crack propagation in the specimen until it enters the rapid energy release stage.In Stage III, the accumulated strain energy was released rapidly and the energy accumulation ratio β reduced drastically as a result of many new cracks having been generated in the specimen.Finally, the specimen failed and lost its bearing capacity, but there was still a small amount of residual strain energy because of the confining pressure effect.

Effect of Confining Pressure
To accurately analyze the influence of confining pressure on the energy accumulation and release in the granite specimen, comparisons of the strain energy and accumulation ratio under different confining pressures are drawn in Figures 6 and 7, respectively.With increasing confining pressure, the maximum value of the strain energy accumulation ratio β decreased gradually.Conversely, the peak value of the accumulated strain energy in the specimen increased linearly [33,34].As demonstrated in Figure 8a, the maximum β decreased according to the exponential function β = 8.485e −0.069σ3 + 90.05, and the peak strain energy increased according to the linear function Es = 0.123σ3 + 1.318.dissipation energy increased gradually and the energy accumulation ratio β decreased distinctly.When the accumulated strain energy reaches the extreme value of energy storage, the excessive energy can accelerate the crack propagation in the specimen until it enters the rapid energy release stage.In Stage III, the accumulated strain energy was released rapidly and the energy accumulation ratio β reduced drastically as a result of many new cracks having been generated in the specimen.Finally, the specimen failed and lost its bearing capacity, but there was still a small amount of residual strain energy because of the confining pressure effect.

Effect of Confining Pressure
To accurately analyze the influence of confining pressure on the energy accumulation and release in the granite specimen, comparisons of the strain energy and accumulation ratio under different confining pressures are drawn in Figures 6 and 7, respectively.With increasing confining pressure, the maximum value of the strain energy accumulation ratio β decreased gradually.Conversely, the peak value of the accumulated strain energy in the specimen increased linearly [33,34].As demonstrated in Figure 8a, the maximum β decreased according to the exponential function β = 8.485e −0.069σ3 + 90.05, and the peak strain energy increased according to the linear function Es = 0.123σ3 + 1.318.As shown in Figure 8b, the pre-peak energy accumulation speed in specimens increased in the of the linear function va = 0.01σ3 + 0.314 with increasing confining pressure.In the post-peak stage, the failure status of specimens was converted from brittle failure to ductile failure with increasing confining pressure.The process of strain energy accumulation and release in the specimen occurred repeatedly under the higher confining pressure, resulting in a significant increase of residual energy (Figure 6).Consequently, the energy release speed in the post-peak stage decreased gradually with increasing confining pressure according to the exponential function vr = 73.2e−0.05σ3 .

Energy Distribution of Microscopic Minerals
Due to the differences in the mechanical properties of the minerals in granite, the energy accumulated by different minerals is significantly different.Figure 9 shows a comparison of the total mineral energy and unit mineral energy of mica, feldspar, and quartz under uniaxial loading.In terms of the magnitude of accumulated strain energy in the minerals of granite, the greatest total energy was accumulated in feldspar, followed by quartz and mica.When the granite was subjected to loading, feldspar played a major part in energy accumulation, and the accumulated strain energy in feldspar accounted for approximately 78.5% of the total strain energy.In contrast, the energies accumulated in mica and quartz were less, only accounting for 16.6% and 4.9% of the total strain As shown in Figure 8b, the pre-peak energy accumulation speed in specimens increased in the form of the linear function v a = 0.01σ 3 + 0.314 with increasing confining pressure.In the post-peak stage, the failure status of specimens was converted from brittle failure to ductile failure with increasing confining pressure.The process of strain energy accumulation and release in the specimen occurred repeatedly under the higher confining pressure, resulting in a significant increase of residual energy (Figure 6).Consequently, the energy release speed in the post-peak stage decreased gradually with increasing confining pressure according to the exponential function v r = 73.2e−0.05σ3 .

Energy Distribution of Microscopic Minerals
Due to the differences in the mechanical properties of the minerals in granite, the energy accumulated by different minerals is significantly different.Figure 9 shows a comparison of the total mineral energy and unit mineral energy of mica, feldspar, and quartz under uniaxial loading.In terms of the magnitude of accumulated strain energy in the minerals of granite, the greatest total energy was accumulated in feldspar, followed by quartz and mica.When the granite was subjected to loading, feldspar played a major part in energy accumulation, and the accumulated strain energy in feldspar accounted for approximately 78.5% of the total strain energy.In contrast, the energies accumulated in mica and quartz were less, only accounting for 16.6% and 4.9% of the total strain energy, respectively.Therefore, feldspar was the main microstructure for energy accumulation in the granite.energy, respectively.Therefore, feldspar was the main microstructure for energy accumulation in the granite.The unit mineral energy can be used to reflect the energy accumulation capacity of a particular mineral in rock.As shown in Figure 9, the unit mineral energy of mica was the largest, followed by that of feldspar, and that of quartz was the smallest.It can be seen from Equations ( 5)-( 7) that the strain energy of the particles in the model is inversely proportional to the stiffness.Among the three minerals, the stiffness of mica was the smallest, and that of quartz was largest.In addition, mica particles did not easily cause failure because of the quartz and feldspar particles with higher strength around them.The mica particles generated a larger elastic (recoverable) deformation.Therefore, the unit mineral energy of mica particles was greater than those of both feldspar and quartz.
It should be noted that the results obtained in this section were based on the assumptions of mineral parameters.Though there were deviations between the numerical model and the actual granite at parameters of minerals, the strength and stiffness relationship among the three minerals were in accordance with the actual rock.Therefore, the results had some reference value on energy evolution research of microstructures in rocks.
Figure 10 shows the relationship between the strain energy of microscopic minerals and the confining pressure.With increasing confining pressure, the total energy of the three types of minerals in granite increased linearly at different rates.The growth rate of total energy in feldspar was the fastest, followed by that of quartz, and that of mica was the slowest.Unlike the total mineral energy, the growth rate of the unit energy of mica was the fastest, followed by that of feldspar, and that of quartz was the slowest.The unit mineral energy can be used to reflect the energy accumulation capacity of a particular mineral in rock.As shown in Figure 9, the unit mineral energy of mica was the largest, followed by that of feldspar, and that of quartz was the smallest.It can be seen from Equations ( 5)-( 7) that the strain energy of the particles in the model is inversely proportional to the stiffness.Among the three minerals, the stiffness of mica was the smallest, and that of quartz was largest.In addition, mica particles did not easily cause failure because of the quartz and feldspar particles with higher strength around them.The mica particles generated a larger elastic (recoverable) deformation.Therefore, the unit mineral energy of mica particles was greater than those of both feldspar and quartz.
It should be noted that the results obtained in this section were based on the assumptions of mineral parameters.Though there were deviations between the numerical model and the actual granite at parameters of minerals, the strength and stiffness relationship among the three minerals were in accordance with the actual rock.Therefore, the results had some reference value on energy evolution research of microstructures in rocks.
Figure 10 shows the relationship between the strain energy of microscopic minerals and the confining pressure.With increasing confining pressure, the total energy of the three types of minerals in granite increased linearly at different rates.The growth rate of total energy in feldspar was the fastest, followed by that of quartz, and that of mica was the slowest.Unlike the total mineral energy, the growth rate of the unit energy of mica was the fastest, followed by that of feldspar, and that of quartz was the slowest.energy, respectively.Therefore, feldspar was the main microstructure for energy accumulation in the granite.The unit mineral energy can be used to reflect the energy accumulation capacity of a particular mineral in rock.As shown in Figure 9, the unit mineral energy of mica was the largest, followed by that of feldspar, and that of quartz was the smallest.It can be seen from Equations ( 5)-( 7) that the strain energy of the particles in the model is inversely proportional to the stiffness.Among the three minerals, the stiffness of mica was the smallest, and that of quartz was largest.In addition, mica particles did not easily cause failure because of the quartz and feldspar particles with higher strength around them.The mica particles generated a larger elastic (recoverable) deformation.Therefore, the unit mineral energy of mica particles was greater than those of both feldspar and quartz.
It should be noted that the results obtained in this section were based on the assumptions of mineral parameters.Though there were deviations between the numerical model and the actual granite at parameters of minerals, the strength and stiffness relationship among the three minerals were in accordance with the actual rock.Therefore, the results had some reference value on energy evolution research of microstructures in rocks.
Figure 10 shows the relationship between the strain energy of microscopic minerals and the confining pressure.With increasing confining pressure, the total energy of the three types of minerals in granite increased linearly at different rates.The growth rate of total energy in feldspar was the fastest, followed by that of quartz, and that of mica was the slowest.Unlike the total mineral energy, the growth rate of the unit energy of mica was the fastest, followed by that of feldspar, and that of quartz was the slowest.

Discussion
In this study, energy accumulation and dissipation in granite under different confining pressures were numerically investigated using a microstructure PFC2D model based on DIP.The evolution of input energy, strain energy and dissipation energy obtained by the simulation corresponded with those observed in laboratory tests [13,33], which indicated that the reconstruction and energy analysis methods of the actual microstructural model of rock are reasonable and effective.Due to the difficult operation on cyclic loading and unloading in the post-peak failure stage of rocks, energy evolution in post-peak stage is hard to get by laboratory tests, while the simulation in this paper can get the whole evolution process.Simulation results are supplement and promotion to laboratory test results.
Furtherly, the energy accumulation of different minerals in the granite was researched, and the main microstructures for energy accumulation were found.However, the results cannot be compared with laboratory tests because there is no effective test method to analyze the energy in minerals.The microstructure model in this paper shows great superiority in analyzing the micro-mechanism of energy evolution in rocks.
Although there are many great advantages in actual microstructure reconstruction using the method presented in this paper, there are still some difficulties in the assignment of the micro parameters of different minerals.In the typical homogeneous or heterogeneous model, the set of micro parameters of rock is easily adjusted to provide a good fit to laboratory test data.However, in the actual microstructure model, there are four sets of parameters for different structures which need to be calibrated, which presents a great difficulty [31].Therefore, how to determine suitable micro parameters is a key problem in the application of this microstructure model based on DIP.

Conclusions
This study proposed a method to reconstruct the microstructural model of granite using Otsu's digital image processing technology and PFC2D.A series of compression tests were carried out to investigate the macro and micro characteristics of energy evolution of the granite specimens.The macro-characteristics of energy evolution in other types of granites or other rocks show the similar law.However, the quantificational function relationship of energy in other rocks will be researched in the future work.And the laboratory test methods on measuring energy in minerals need to be further developed.The following conclusions have been drawn from this study.

Discussion
In this study, energy accumulation and dissipation in granite under different confining pressures were numerically investigated using a microstructure PFC2D model based on DIP.The evolution of input energy, strain energy and dissipation energy obtained by the simulation corresponded with those observed in laboratory tests [13,33], which indicated that the reconstruction and energy analysis methods of the actual microstructural model of rock are reasonable and effective.Due to the difficult operation on cyclic loading and unloading in the post-peak failure stage of rocks, energy evolution in post-peak stage is hard to get by laboratory tests, while the simulation in this paper can get the whole evolution process.Simulation results are supplement and promotion to laboratory test results.
Furtherly, the energy accumulation of different minerals in the granite was researched, and the main microstructures for energy accumulation were found.However, the results cannot be compared with laboratory tests because there is no effective test method to analyze the energy in minerals.The microstructure model in this paper shows great superiority in analyzing the micro-mechanism of energy evolution in rocks.
Although there are many great advantages in actual microstructure reconstruction using the method presented in this paper, there are still some difficulties in the assignment of the micro parameters of different minerals.In the typical homogeneous or heterogeneous model, the set of micro parameters of rock is easily adjusted to provide a good fit to laboratory test data.However, in the actual microstructure model, there are four sets of parameters for different structures which need to be calibrated, which presents a great difficulty [31].Therefore, how to determine suitable micro parameters is a key problem in the application of this microstructure model based on DIP.

Conclusions
This study proposed a method to reconstruct the microstructural model of granite using Otsu's digital image processing technology and PFC2D.A series of compression tests were carried out to investigate the macro and micro characteristics of energy evolution of the granite specimens.The macro-characteristics of energy evolution in other types of granites or other rocks show the similar law.However, the quantificational function relationship of energy in other rocks will be researched in the future work.And the laboratory test methods on measuring energy in minerals need to be further developed.The following conclusions have been drawn from this study.

Figure 2 .
Figure 2.An example of incorporating a gray image into the particle flow code (PFC) model.

Figure 3 .
Figure 3.Comparison between numerical simulation and test data.

Figure 2 .
Figure 2.An example of incorporating a gray image into the particle flow code (PFC) model.

Figure 2 .
Figure 2.An example of incorporating a gray image into the particle flow code (PFC) model.

Figure 3 .
Figure 3.Comparison between numerical simulation and test data.

Figure 3 .
Figure 3.Comparison between numerical simulation and test data.After calibrating the parameters of the model, the established microstructural model of granite was used to carry out a series of compression tests under different confining pressures to further analyze the macro and micro characteristics of energy evolution in granite.In tests, the lateral and bottom loading plates were used to apply constant confining pressure on the specimens based on servo commands.The confining pressure was set to 0, 10, 20, 30, or 40 MPa.After the confining pressure

MF
are the normal force, shear force, and the moment in the parallel bond i, respectively; n i F and s i are the normal force and shear force in the contact i, respectively; Ai and Ii are the area and inertia moment of the bond cross section, respectively; Nc is the number of contacts; and Npb is the number of parallel bonds.

Figure 4 .
Figure 4.The calculation model of strain energy in particle flow code software (PFC2D).Figure 4. The calculation model of strain energy in particle flow code software (PFC2D).

Figure 4 .
Figure 4.The calculation model of strain energy in particle flow code software (PFC2D).Figure 4. The calculation model of strain energy in particle flow code software (PFC2D).

Figure 6 .
Figure 6.Curves of strain energy under different confining pressures.

Figure 7 .
Figure 7. Curves of the energy accumulation ratio under different confining pressures.

Figure 6 .
Figure 6.Curves of strain energy under different confining pressures.

Figure 6 .
Figure 6.Curves of strain energy under different confining pressures.

Figure 7 .
Figure 7. Curves of the energy accumulation ratio under different confining pressures.Figure 7. Curves of the energy accumulation ratio under different confining pressures.

Figure 7 .
Figure 7. Curves of the energy accumulation ratio under different confining pressures.Figure 7. Curves of the energy accumulation ratio under different confining pressures.

Figure 8 .
Figure 8. Correlation curves of (a) energy parameters at the peak stress point and (b) energy velocities with different confining pressures.

Figure 8 .
Figure 8. Correlation curves of (a) energy parameters at the peak stress point and (b) energy velocities with different confining pressures.

Figure 9 .
Figure 9.Comparison of strain energy accumulated in the different minerals of granite.

Figure 9 .
Figure 9.Comparison of strain energy accumulated in the different minerals of granite.

Figure 9 .
Figure 9.Comparison of strain energy accumulated in the different minerals of granite.

Figure 10 .
Figure 10.Energy curves of different minerals with increasing confining pressure: (a) total mineral energy and (b) unit mineral energy.

Figure 10 .
Figure 10.Energy curves of different minerals with increasing confining pressure: (a) total mineral energy and (b) unit mineral energy.

Table 1
Parameters of the microstructures in the granite model.

Table 1
Parameters of the microstructures in the granite model.