Zhan Zhifa, Wang Jilin, Du Jiahu, Wang Bibo, Yang Zhigang, Huang Xiaolin
China Renewable Energy Engineering Institute Limited Company, Andingmenwai Street No. 57A, Dongcheng District, Beijing, 100011, China.
Department of Civil Engineering, School of Human Settlements and Civil Engineering, Xi'an Jiaotong University, Xianning West Road No. 28, Xi'an, 710049, Shaanxi, China.
Sci Rep. 2025 Jul 1;15(1):22445. doi: 10.1038/s41598-025-02322-2.
This study investigates the synergistic effects of crystal habit and confining pressure on the compressive mechanical behavior of crystalline rocks using a grain-based model implemented in the Universal Distinct Element Code. Synthetic granite samples with euhedral, subhedral, and anhedral microstructures were simulated by adjusting grain boundary irregularity, quantified by fractal dimension (D), which increased from euhedral to anhedral. Deformation and failure were depicted using a compression-hardening contact model and a cohesion-weakening-friction-strengthening material model. Results show that Young's modulus and Poisson's ratio nonlinearly increased with confining pressure, insensitive to crystal habit. Under low confining pressure, compressive stress-strain curves exhibited significant dispersion during crack damage due to crystal habit variations. However, pre-peak curves increasingly overlapped at higher pressures, while post-peak curves remained dispersed. Peak compressive strength initially decreased and then increased with D under low pressure, but this trend diminished as pressure rose. Strength followed the Hoek-Brown criterion, with envelopes dispersing at low pressure but converging under higher pressure. Hoek-Brown parameters first increased and then decreased with D. Tensile stress magnitude and the number of concentration zones increased with increasing D but were suppressed by higher pressure, reducing their influence. As confining pressure increased, tensile stress concentration effects were inhibited, while grain contact shear strength increased, eventually matching grain material strength, diminishing the role of crystal habit and making confining pressure dominant. Additionally, with increasing D and confining pressure, intergranular cracks decreased in number and length, while intragranular damage zones became denser, highlighting the interplay between microstructure heterogeneity and stress state in controlling failure mechanisms.
本研究使用通用离散元代码中实现的基于颗粒的模型,研究晶体习性和围压对晶体岩石压缩力学行为的协同作用。通过调整由分形维数(D)量化的晶界不规则性来模拟具有自形、半自形和他形微观结构的合成花岗岩样品,分形维数从自形到他形逐渐增加。使用压缩硬化接触模型和内聚力弱化-摩擦强化材料模型来描述变形和破坏。结果表明,杨氏模量和泊松比随围压非线性增加,对晶体习性不敏感。在低围压下,由于晶体习性变化,压缩应力-应变曲线在裂纹损伤期间表现出显著的离散性。然而,在较高压力下,峰值前曲线越来越重叠,而峰值后曲线仍然离散。在低压下,峰值抗压强度最初随D减小然后增加,但随着压力升高,这种趋势减弱。强度遵循Hoek-Brown准则,在低压下包络线分散,但在较高压力下收敛。Hoek-Brown参数首先随D增加然后减小。拉伸应力大小和集中区数量随D增加而增加,但受到较高压力的抑制,从而降低了它们的影响。随着围压增加,拉伸应力集中效应受到抑制,而颗粒接触剪切强度增加,最终与颗粒材料强度匹配,削弱了晶体习性的作用,使围压起主导作用。此外,随着D和围压增加,晶间裂纹数量和长度减少,而晶内损伤区变得更密集,突出了微观结构非均质性和应力状态在控制破坏机制中的相互作用。