Welcome to Francis Academic Press

Academic Journal of Computing & Information Science, 2023, 6(2); doi: 10.25236/AJCIS.2023.060216.

PFC2D numerical simulation of sandstone damage evolution based on MATLAB image processing


Yingying Zhu

Corresponding Author:
Yingying Zhu

Shandong University of Science and Technology, Qingdao, China


In this paper, mesoscopic parameters of sandstone are obtained by SEM and image processing, and mechanical parameters of sandstone are obtained by uniaxial compression experiment, based on which the PFC2D model of red sandstone is established. Through numerical simulation of uniaxial compression, crack propagation of sandstone was observed, and parameters such as stress-strain curve and number of bond fractures were obtained. Fractal dimension of sandstone was calculated by MATLAB, and the relationship between damage degree, fractal dimension and crack propagation was analyzed. The results show that the crack development of sandstone can be divided into four stages: elastic deformation, crack initiation, crack propagation and penetration. The damage degree corresponds to each stage of fracture development, the minimum is in elastic deformation stage, and the maximum is in penetrating stage. The fractal dimension changes differently in each stage, and shows an increasing trend in the germination stage. In the expansion stage, it shows a downward trend. At the through-through stage, there is a substantial growth trend. Therefore, the box dimension can be used to indirectly evaluate the internal damage and study the damage evolution characteristics, which provides a new idea for studying the damage evolution of sandstone and has important significance for evaluating the stability of engineering rock mass.


MATLAB image processing; PFC2D; fractal dimension; damage degree; red sandstone

Cite This Paper

Yingying Zhu. PFC2D numerical simulation of sandstone damage evolution based on MATLAB image processing. Academic Journal of Computing & Information Science (2023), Vol. 6, Issue 2: 120-126. https://doi.org/10.25236/AJCIS.2023.060216.


[1] B. Kang, C. Yanan, L. Jian, “The unloading failure characteristics of shale under different water absorption time using the PFC numerical method,” Rock and Soil Mechanics, Vol.41, pp. 355-367, July 2020. 

[2] Y. Shengqi, S. Bowen, T. Wenling, “ Discrete element simulation of the mechanical properties of shale with different bedding inclinations under conventional triaxial compression,” Journal of Engineering Science, Vol.44, pp. 430-439, March 2022.

[3] H. Yanhua, Y. Shengqi, “Discrete element study on strength and failure behavior of jointed sandstone with two sets of cross-joint,” Journal of China Coal Society , Vol. 40, pp. 76-84, June 2015.

[4] D. Qinghai, H. Shanxiang, X. Yongqing, “Particle flow simulation of uniaxial compression failure of prefabricated cracked rock,” Hydropower energy science, Vol.35, pp. 95-98, November 2017

[5] D. Haigang, S. Jianwei, X. Jun, “Effect of weak acid and surfactant on coal pore and fractal characteristics,” Mining Safety and Environmental Protection, Vol.48, pp.12-16, June 2021.

[6] X. Xingjun, “Study on pore distribution characteristics of pervious concrete based on ct scan test,” Bulletin of the Chinese Ceramic Society, Vol.38, pp.3670-3674, November, 2019.

[7] G. Shuang, Z. Yixin, W. Zhen, “Effect of bedding on the fractal characteristics of dynamic crack propagation in coal rocks,” Journal of China Coal Society, Vol.46, pp.2574-2582, August 2021.

[8] L. Zezheng, L. Jiangteng, “Fracture fractal and crack propagation rate of layered slate,” Mining and Metallurgical Engineering, Vol.37, pp.15-18, October 2017.

[9] D. Ziwei, L. Xiaofei, T. Qingbao, “Study on correlation between fractal characteristics of pore distribution and strength of sandstone particles,” Chinese Journal of Rock Mechanics and Engineering, Vol. 39, pp.1787-1796, September 2020.

[10] Z. Fei, L. Haibo, Z. Qingchun, “Fractal features study of rock-like material damage based on gray correlation of digital images,” Rock and Soil Mechanics, Vol.33, pp.731-738, March 2012.

[11] Z. Yu, W. Shunchuan, X. Xueliang, “Particle flow analysis of acoustic emission characteristics during rock failure process,” Chinese Journal of Rock Mechanics and Engineering, 2013, Vol.32, pp. 951-959, May 2013.

[12] P. Launeau, P.-Y. F. Robin, “Fabric analysis using the intercept method,” Tectonophysics, Vol.267, pp.91-119, 1996.

[13] S. Hongliang, W. Baoren, “Bend Measuring Method Based on MATLAB Image Processing,” Journal of Shandong University of Science and Technology (Natural Science), Vol.36, pp. 83-90, June 2017.

[14] R. Haoyu, L. Guichen, L. Dongxu, “Particle flow simulation of mechanical properties of high stress rock under the influence of stress path,” Chinese Journal of Mining and Safety Engineering, Vol.39, pp. 163-173, January 2022.

[15] H. Yanhua, Y.Shengqi, “Particle flow analysis on crack coalescence behavior of sandstone specimen containing three pre-existing fissures under uniaxial compression,” Chinese Journal of Applied Basic and Engineering Sciences, Vol.24, pp.1232-1247, December 2016

[16] H. Dongmei, C. Xikun, L. Xiaofei, F. Yingxial, “The Fractal Dimension of Rock Crack Under Uniaxial Compression”. Journal of Shandong University of Science and Technology (Natural Science), Vol.33, pp. 58-62, April 2014.

[17] L. Xuelian, L. Xueying, L. Xinchao, “Experimental on indirect tensile test for asphalt mixture based on discrete element simulation and fractal dimension,” Journal of Chang'an University (Natural Science Edition), Vol.39, pp. 39-48, November 2019

[18] D. Zhengding, X. Shuai, Z. Jianrong, “Rate correlation and deformation of damage evolution of non-penetrating fractured rock masse,” Explosion and Impact, Vol.39, pp. 118-129, August 2019.

[19] G. Mudan, Y. Fenfang, “Research on the strength damage law of rock mass containing non-coalescence structural planes,” Jiangsu Architecture, pp. 88-91, August 2018