Sensitivity analysis of fracture energies for the combined finite-discrete element method (FDEM)
Penghai Deng,Quansheng Liu,Xing Huang,Yin Bo,Qi Liu,Weiwei Li
DOI: https://doi.org/10.1016/j.engfracmech.2021.107793
IF: 5.4
2021-06-01
Engineering Fracture Mechanics
Abstract:<p>The combined finite-discrete element method (FDEM) has been widely used in the numerical study of rock mechanics and geotechnical engineering. Selection of the correct parameter values is a prerequisite to ensure the accuracy of the simulation results. However, current calibration procedures are relatively cumbersome, and the parameters are strongly dependent on the element size. To eliminate the dependence of parameter values on the element size, this paper proposes a new constitutive model based on the stress-strain relationship. In addition, uniaxial compression and direct tension simulations are used to investigate the sensitivities of type II and type I fracture energies, i.e., <em>G</em><sub>II</sub> and <em>G</em><sub>I</sub>, including the influence of the <em>G</em><sub>I</sub> (or <em>G</em><sub>II</sub>) value on the simulation results of direct tension (or uniaxial compression) and the influence of different parameters (such as Young's modulus <em>E</em>, Poisson's ratio <em>ν</em>, cohesion <em>c</em>, internal friction angle <em>φ<sub>i</sub></em>, tensile strength <em>f<sub>t</sub></em> and element size <em>h</em>) on the selection of the values of <em>G</em><sub>II</sub> and <em>G</em><sub>I</sub>. The following results are obtained. <a class="workspace-trigger" href="#e0005">(1)</a> As <em>G</em><sub>I</sub> or <em>G</em><sub>II</sub> increases, the strength of the rock sample increases, and the plastic characteristics are more obvious, but the failure mode tends to be stable. <a class="workspace-trigger" href="#e0010">(2)</a> As the <em>E</em> value increases, the value of <em>G</em><sub>I</sub> or <em>G</em><sub>II</sub> decreases as a power function. As the value of <em>f<sub>t</sub></em> or <em>c</em> increases, the value of <em>G</em><sub>I</sub> or <em>G</em><sub>II</sub> increases as a power function. However, other parameters, including <em>ν</em>, <em>φ<sub>i</sub></em> and <em>h</em>, have no effect on the selection of the values of <em>G</em><sub>II</sub> and <em>G</em><sub>I</sub>. <a class="workspace-trigger" href="#e0015">(3)</a> Actual uniaxial compression and triaxial compression simulations prove that the calibrated <em>G</em><sub>II</sub> and <em>G</em><sub>I</sub> values are reasonable. In this study, only two fracture energies, <em>G</em><sub>II</sub> and <em>G</em><sub>I</sub>, must be calibrated; other macroparameters can be obtained through laboratory tests, and other microparameters can be obtained from previous empirical or theoretical values, which greatly improves the calibration efficiency of FDEM input parameters.</p>
mechanics