Veined Rock Performance under Uniaxial and Triaxial Compression Using Calibrated Finite Element Numerical Models
Abstract
:1. Introduction
2. Calibration Methodology
- Unveined UCS and TCS numerical calibration and validation of the FEM model configuration using LdB granite:
- Model configuration: external boundary; material boundaries defining dimensions of specimen and loading platens; position of x- and y-pinned and x- or y-restrained external and material boundary nodes and elements are defined to simulate physical laboratory conditions.
- Sensitivity analysis to model settings: Mesh element density, displacement rate, and stiffness properties of the platen–specimen joint element are varied.
- Comparison of numerical results to physical laboratory test data: Model settings are selected based on their relative influence on the model and whether the numeric output matches target output.
- Unveined and veined UCS calibration and TCS validation using CMET mafic wallrock with single sulfide veins:
- Calibration of unveined specimen: Strength and stiffness properties of laboratory UCS test data compared to UCS model results.
- Incorporation of vein: Vein is incorporated into UCS model, first represented only by the intact vein material, followed by including joint elements at vein–wallrock contacts.
- Sensitivity analysis of inputs: Material properties of all components of the UCS models defined including wallrock and vein materials, and joint elements at vein–wallrock contacts. The solution is indeterminate, so each parameter cannot be individually calibrated. Therefore, a sensitivity analysis is employed to evaluate the effect of each unknown variable on peak strength of the specimen.
- Calibration of veined models to laboratory results: Inputs are adjusted based on the sensitivity analysis results to achieve a set of calibrated input parameters that produce UCS peak strength results that match the physical laboratory test data.
- Sensitivity analysis of model parameters on outputs: Calibrated parameters are systematically varied to ensure one parameter alone is not heavily influencing numeric output.
- Validation: TCS models are used to validate the output Mohr–Coulomb strength parameters against reported laboratory test data.
3. Model Configuration
4. Calibration of Lac du Bonnet Granite Models
4.1. Definition of Calibration Target
4.2. Calibration of Mesh Density
4.3. Calibration of Applied Displacement Rate
4.4. Stiffness and Strength of Platen-Specimen Interface
4.5. Validation of LdB Granite Models
5. Calibration of CMET Models
5.1. Numerical Input Parameters
- For all materials, there is zero residual cohesion (cr) or residual tensile strength (σtr) after yield as yielding in RS2 represents shear or tensile failure of the material at a mesh node. Although RS2 cannot simulate true detachment of nodes because a FEM model is fundamentally a continuum [26], it can be assumed that these values drop to zero after yield. Based on guidance from Li and Bahrani [30], values of 0.1 MPa were assigned to cr and σtr in all three materials to avoid numerical convergence errors.
- Residual friction angle (ϕr) is equal to the calibrated peak friction angle (ϕp), following guidance from Markus [29] and confirmed by a sensitivity analysis in this study comparing peak strength results, which showed that changing ϕr by ±10° resulted in a maximum difference in peak strength of just 3% when ϕr > ϕp and a 0% difference when ϕp > ϕr.
- For the vein–wallrock contact joint element, the Kn:Ks ratio was maintained at 1:1 because the contact is intact [34].
5.2. Published Test Data of CMET Unit
5.3. Numerical Sensitivity Analysis
5.3.1. Selection of Base Case Inputs
5.3.2. Sensitivity Analysis: First Iteration
5.3.3. Sensitivity Analysis: Second Iteration
5.3.4. Sensitivity Analysis: Third Iteration
- (i)
- The vein–wallrock joint elements control specimen failure when using the base case inputs. This is evidenced when cp(v) is increased, as the peak strength of the specimen does not change. This indicates that it has exceeded the strength of the joint element and therefore increases in material strength cannot improve the peak strength output. Additionally, when cp(j) increases, the peak strength slightly increases. This suggests that while the strength of the joint element may have the most influence on the specimen failure, it is not so significant to cause major increases in peak strength given the critical orientation of the vein in these models.
- (ii)
- Kn(j), Ks(j), and E(v) influence the specimen peak strength. When stiffness is increased, the rate at which stress accumulates in the specimen increases and the allowable deformation in the specimen before yield decreases. However, it is important that system stiffness and particularly Kn and Ks are high enough to enable adequate stress transfer through the vein and into the wallrock material below the vein to achieve a realistic numerical simulation.
5.4. Calibration Results
5.4.1. Validation of Unveined CMET Models
5.4.2. Validation of Veined CMET Models
6. Numerical TCS Experiments on Veined CMET
6.1. Modelling Program
6.2. Numerical Experiment Results
6.2.1. Specimen Elastic Response
6.2.2. Failure Types
- (i)
- Type A failures: Failure primarily occurs through both the vein material and vein–wallrock contact joint element. Minor material yield in the wallrock may occur near the ends of the specimen.
- (ii)
- Type B failures: Failure primarily occurs through the vein material but not the vein–wallrock contact joint element. Minor material yield in the wallrock may occur near the vein–wallrock contact.
- (iii)
- Type C failures: Failure occurs in wide bands of predominantly shear through the wallrock material on one or both sides of the vein. Some failure may occur in the vein material, but no failure occurs in the vein–wallrock contact joint element.
6.2.3. Effects of Vein Orientation and Thickness
6.2.4. Intact Vein Shear Strength
7. Discussion
7.1. Model Evaluation
7.2. Calibration Results
7.3. Numerical UCS and TCS Experiments
- (1)
- this study shows peak strength of the specimen is reduced even by the presence of non-critically oriented weakening veins (55° < α < 90°);
- (2)
- this study shows the vein orientation resulting in the weakest peak strength ranges from 20–30°, depending on the magnitude of confining stress (σ3).
8. Conclusions
Author Contributions
Funding
Data Availability Statement
Acknowledgments
Conflicts of Interest
Abbreviations
2D | 2-dimensional |
3D | 3-dimensional |
c | Cohesion (MPa) |
CMET | Complejo Máfico El Teniente/El Teniente Mafic Complex |
E | Young’s modulus (MPa) |
E2D | Young’s modulus for plane strain analysis (MPa) |
E3D | Young’s modulus measured from physical laboratory test (MPa) |
FEM | Finite element method |
ISRM | International Society for Rock Mechanics and Rock Engineering |
Kn | Joint normal stiffness (MPa/m) |
Ks | Joint shear stiffness (MPa/m) |
LdB | Lac du Bonnet granite |
TCS | Triaxial compressive strength |
UCS | Uniaxial (or Unconfined) compressive strength |
α | Angle between core axis and vein in diamond drill core |
γ | Unit weight (MN/m3) |
μ | Coefficient of friction |
ν | Poisson’s ratio |
ν2D | Poisson’s ratio for plane strain analysis |
ν3D | Poisson’s ratio measured from physical laboratory test |
ϕ | Friction angle (degrees, °) |
σ1 | Major principal stress (MPa) |
σ3 | Minor principal stress (MPa) |
σc | Peak compressive strength (MPa) |
σn | Normal stress (MPa) |
τ | Shear stress (MPa) |
References
- Day, J.J. Brittle overbreak prediction in deep excavations for hydrothermally altered and heterogeneous rockmasses. Bull. Eng. Geol. Environ. 2020, 79, 1041–1060. [Google Scholar] [CrossRef]
- Hajiabdolmajid, V.; Kaiser, P.K.; Martin, C.D. Modelling brittle failure of rock. Int. J. Rock Mech. Min. Sci. 2002, 39, 731–741. [Google Scholar] [CrossRef]
- Day, J.J. How to incorporate variability of rockmass structures into equivalent continuum numerical models using Composite Geological Strength Index. In Proceedings of the Rocscience International Conference; The Evolution of Geotech—25 Years of Innovation; CRC Press/Balkema: Leiden, The Netherlands, 2021; pp. 110–116. [Google Scholar] [CrossRef]
- Hoek, E.; Carter, T.G.; Diederichs, M.S. Quantification of Geological Strength Index chart. In Proceedings of the 47th US Rock Mechanics/Geomechanics Symposium, San Francisco, CA, USA, 23–26 June 2013; American Rock Mechanics Association: Alexandria, VA, USA, 2013; pp. 1757–1764. [Google Scholar]
- Bieniawski, Z.T. Engineering Rock Mass Classifications: A Complete Manual for Engineers and Geologists in Mining, Civil, and Petroleum Engineering; Wiley: New York, NY, USA, 1989. [Google Scholar]
- Brzovic, A.; Villaescusa, E. Rock mass characterization and assessment of block-forming geological discontinuities during caving of primary copper ore at the El Teniente mine, Chile. Int. J. Rock Mech. Min. Sci. 2007, 44, 565–583. [Google Scholar] [CrossRef]
- Clark, M.D.; Day, J.J. Mineralogical and sample selection implications for geomechanical properties of intact heterogeneous and veined rocks from the Legacy skarn deposit. Eng. Geol. 2021, 285, 106067. [Google Scholar] [CrossRef]
- Turichshev, A.; Hadjigeorgiou, J. Quantifying the effects of vein mineralogy, thickness, and orientation on the strength of intact veined rock. Eng. Geol. 2017, 226, 199–207. [Google Scholar] [CrossRef]
- Hoek, E.; Carranza-Torres, C.; Corkum, B. Hoek-Brown failure criterion—2002 edition. In Proceedings of the North American Rock Mechanics Symposium—Tunnelling Association of Canada Conference, Toronto, ON, Canada, 7–10 July 2002; pp. 267–273. Available online: https://www.rocscience.com/assets/resources/learning/hoek/Hoek-Brown-Failure-Criterion-2002.pdf (accessed on 2 April 2020).
- Hoek, E.; Brown, E.T. The Hoek-Brown failure criterion ant GSI—2018 edition. J. Rock Mech. Geotech. Eng. 2019, 11, 445–463. [Google Scholar] [CrossRef]
- Day, J.J.; Diederichs, M.S.; Hutchinson, D.J. Optimization of structural contact stiffness and strength for discrete simulation of progressive failure of healed structure. Geomech. Tunnelbau 2015, 8, 414–420. [Google Scholar] [CrossRef]
- Sanipour, S.; Bahrani, N.; Corkum, A. Simulation of brittle failure around Canada’s mine-by experiment tunnel using 2D continuum-based Voronoi tessellated models. Rock Mech. Rock Eng. 2022, 55, 6387–6408. [Google Scholar] [CrossRef]
- Goodman, R.E.; Taylor, R.L.; Brekke, L. A model for the mechanics of jointed rock. J. Soil Mech. Found. Div. 1968, 94, 637–659. [Google Scholar] [CrossRef]
- Coulomb, C.A. Esai Sur Une Application res degeles des Maximum et Minimis des Stalique Relations a L’Architecture. Mémoires de Mathematique de l’Academie Royale de Science Paris 1776, 7, 343–382. [Google Scholar]
- Parry, R.H.G. Mohr Circles, Stress Paths and Geotechnics, 2nd ed.; Spoon Press: London, UK; New York, NY, USA, 2004. [Google Scholar]
- Barton, N.; Bandis, S. Review of predictive capabilities of JRC-JCS model in engineering practice. In Proceedings of the International Symposium on Rock Joints (Barton & Stephansson eds.), Loen, Norway, 4–6 June 1990 ; pp. 603–610. [Google Scholar]
- Rocscience. RS2: 2D Finite Element Analysis, v11.003. 2020. Available online: https://www.rocscience.com/ (accessed on 14 March 2020).
- Lajtai, E.Z. Shear strength of weakness planes in rock. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1969, 6, 499–515. [Google Scholar] [CrossRef]
- Bustos, N.; van Sint Jan, M.; Seguel, J.; Cavieres, P. Factors to analyse strength of mineral infilled rock mass discontinuities. J. Min. Eng. Res. 2019, 1, 166–174. [Google Scholar] [CrossRef]
- Rudderham, G.A. Geomechanics of Hydrothermal Veins: Insights from Laboratory Direct Shear and Numerical Triaxial Experiments. Master’s Thesis, Queen’s University, Kingston, ON, Canada, 2022. Available online: https://qspace.library.queensu.ca/handle/1974/31376 (accessed on 9 January 2023).
- Martin, C.D. The Strength Massive Lac du Bonnet Granite around Underground Openings. Ph.D. Thesis, University of Manitoba, Winnipeg, MB, Canada, 1993. Available online: http://hdl.handle.net/1993/9785 (accessed on 12 April 2022).
- Labeid, M.T.A. Improvement of Laboratory Protocols for Discontinuity Geomechanical Characterization and Investigation of the Effect of Saturation on Granite Strength. Master’s Thesis, Queen’s University, Kingston, ON, Canada, 2019. Available online: https://qspace.library.queensu.ca/handle/1974/26635 (accessed on 12 April 2022).
- de los Santos Valderrama, C.G. Effecto de la Mineralogîa, Alteración y Geometría en la Resistencia Mecánica de las Vetillas, mena El Teniente, Región del Libertador Bernardo O’Higgins, Chile. Master’s Thesis, Universidad de Concepción, Concepción, Chile, 2011. [Google Scholar]
- ISRM. International society for rock mechanics commission on standardization of laboratory and field tests: Suggested methods for determining the strength of rock materials in triaxial compression. Int. J. Rock Mech. Min. Sci. 1978, 15, 47–51. [Google Scholar] [CrossRef]
- Bieniawski, Z.T.; Bernede, M.J. International Society for Rock Mechanics Suggested methods for determining the uniaxial compressive strength and deformability of rock materials: Part 1. Suggested method for determination of the uniaxial compressive strength of rock materials. Int. J. Rock Mech. Min. Sci. 1979, 16, 137. [Google Scholar] [CrossRef]
- Jing, L. A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering. Int. J. Rock Mech. Min. Sci. 2003, 40, 283–353. [Google Scholar] [CrossRef]
- Askeland, D.R.; Fulay, P.P.; Wright, W.J. The Science and Engineering of Materials, 6th ed.; Cengage Learning: Stamford, CT, USA, 2011. [Google Scholar]
- Knowles, K.M. The plane strain Young’s modulus in cubic materials. J. Elast. 2017, 128, 147–173. [Google Scholar] [CrossRef]
- Markus, S. Stepwise Approach to Verification of the Combined Finite-Discrete Element Method for Modelling Instability around Tunnels in Brittle Rock. Master’s Thesis, Queen’s University, Kingston, ON, Canada, 2019. Available online: https://qspace.library.queensu.ca/handle/1974/27437 (accessed on 12 April 2022).
- Li, Y.; Bahrani, N. A continuum grain-based model for intact and granulated Wombeyan marble. Comput. Geotech. 2021, 129, 103872. [Google Scholar] [CrossRef]
- Farahmand, K. Characterization of Rockmass Properties and Excavation Damage Zone (EDZ) Using a Synthetic Rock Mass (SRM) Approach. Ph.D. Thesis, Queen’s University, Kingston, ON, Canada, 2017. Available online: https://qspace.library.queensu.ca/handle/1974/22775 (accessed on 12 April 2022).
- Rabbat, B.G.; Russell, H.G. Friction coefficient of steel on concrete or grout. J. Struct. Eng. 1985, 111, 505–515. [Google Scholar] [CrossRef]
- Read, J.; Stacey, P. Guidelines for Open Pit Slope Design; CSIRO, CRC Press/Balkema: Leiden, The Netherlands, 2009. [Google Scholar]
- Day, J.J.; Diederichs, M.S.; Hutchinson, D.J. Component and system deformation properties of complex rockmasses with healed structure. In Proceedings of the 48th US Rock Mechanics/Geomechanics Symposium, Minneapolis, MN, USA, 1–4 June 2014; American Rock Mechanics Association: Alexandria, VA, USA, 2014. 11p. [Google Scholar]
- Day, J.J.; Diederichs, M.S.; Hutchinson, D.J. New direct shear testing protocols and analyses for fractures and healed intrablock rockmass discontinuities. Eng. Geol. 2017, 229, 53–72. [Google Scholar] [CrossRef]
- Simmons, G.; Birch, F. Elastic Constants of Pyrite. J. Appl. Phys. 1963, 34, 2736–2738. [Google Scholar] [CrossRef]
- Birch, F. Compressibility: Elastic constants. In Handbook of physical constants, Geological Society of America Memoir; Clark, S.P., Jr., Ed.; Geological Society of America: New York, NY, USA, 1966; Volume 97, pp. 97–173. [Google Scholar] [CrossRef]
- LePage, Y.; Rodgers, J.R. Ab initio elasticity of FeS2 pyrite from 0 to 135 GPa. Phys. Chem. Miner. 2005, 32, 564–567. [Google Scholar] [CrossRef]
- Prasad, S.C.; Wooster, W.A. The Elasticity of Iron Pyrites, FeS2. Acta Crystallogr. 1956, 9, 169–173. [Google Scholar] [CrossRef]
- Whitaker, M.L.; Liu, W.; Wang, L.; Li, B. Acoustic velocities and elastic properties of pyrite (FeS2) to 9.6 GPa. J. Earth Sci. 2010, 21, 792–800. [Google Scholar] [CrossRef]
- Lama, R.D.; Vutukuri, V.S. Handbook on Mechanical Properties of Rocks; Trans Tech Publications: Clausthal, Germany, 1978; Volume 2. [Google Scholar]
- Jaeger, J.C.; Cook, N.G.W. Fundamentals of Rock Mechanics, 3rd ed.; Chapman and Hall: London, UK, 1979. [Google Scholar]
- Hudson, J.A.; Harrison, J.P. Engineering Rock Mechanics: An Introduction to the Principles, 1st ed.; Elsevier Science & Technology: Amsterdam, The Netherland, 1997. [Google Scholar]
- Patton, F.D. Multiple modes of shear failure in rocks. In Proceedings of the 1st Congress of the International Society for Rock Mechanics, Lisbon, Portugal, 25 September–1 October 1966; Volume 1, pp. 509–513. [Google Scholar]
- Puzrin, A.M. Constitutive Modelling in Geomechanics: Introduction; Springer: Berlin/Heidelber, Germany, 2012. [Google Scholar] [CrossRef]
- Sillitoe, R.H. Porphyry copper systems. Econ. Geol. 2010, 105, 3–41. [Google Scholar] [CrossRef]
Parameter (Units) | Value |
---|---|
Unit weight, γ (MN/m3) | 0.077 |
Young’s modulus, E (GPa) | 200 |
Poisson’s ratio, ν | 0.27 |
Parameter (Units) | Physical Laboratory Test Data | Numerical Input Value | Calibrated Numerical Output |
---|---|---|---|
Young’s modulus, E (GPa) | 69.54 *; 69 ± 5.8 † | 64.83 § | 69.4 |
Poisson’s ratio, ν | 0.36 *; 0.26 ± 0.04 † | 0.26 § | 0.35 |
Peak uniaxial compressive strength, σc (MPa) | 228 *; 200 ± 22 † | - | 228.5 |
Tensile strength, σt (MPa) | 9.3 ± 1.3 † | 9.3 | 8.9 ** |
Cohesion, c (MPa) | 30 †,‡ | 30.4 | 29.0 |
Friction angle, ϕ (°) | 59 †,‡ | 62.2 | 61.6 |
Unit weight, γ (MN/m3) | N/A | 0.027 | - |
Base Parameter (Units) | CMET Wallrock Material (w) | Pyrite Vein Material (v) | Vein–Wallrock Contact Joint Element (j) |
---|---|---|---|
Unit weight (MN/m3) | γ(w) = 0.027 | γ(v) = 0.047 | - |
Young’s modulus (GPa) | E(w) = 58.5 | E(v) = 125 | - |
Poisson’s ratio | ν(w) = 0.14 | ν(v) = 0.17 | - |
Peak friction angle (°) | ϕp(w) = 32 | ϕp(v) = 48 | ϕp(j) = 42 |
Residual friction angle (°) | ϕr(w) = 32 | ϕr(v) = 48 | ϕr(j) = 42 |
Peak cohesion (MPa) | cp(w) = 43 | cp(v) = 26 | cp(j) = 22 |
Residual cohesion (MPa) | cr(w) = 0.1 | cr(v) = 0.1 | cr(j) = 0.1 |
Peak tensile strength (MPa) | σtp(w) = 14 | σtp(v) = 7.5 | σtp(j) = 5 |
Residual tensile strength (MPa) | σtr(w) = 0.1 | σtr(v) = 0.1 | σtr(j) = 0.1 |
Normal stiffness (GPa/m) | - | - | Kn(j) = 1750 |
Shear stiffness (GPa/m) | - | - | Ks(j) = 1750 |
Parameter (Units) | VDA-04 | VDA-06 | VDA-09 |
---|---|---|---|
Vein primary mineral % | Pyrite 45% | Chalcopyrite 67% | Pyrite 60% |
Vein secondary mineral % | Chalcopyrite 27% | Quartz 22% | Chalcopyrite 24% |
Vein tertiary mineral % | Quartz 24% | Pyrite 6% | Quartz 8% |
Vein orientation, α (°) | 22 | 20 | 30 |
Vein thickness (mm) | 4 | 1.5 | 6 |
Wallrock lithology | CMET | CMET | CMET |
σ3 (MPa) | 5 | 0 | 15 |
σ1peak (MPa) | 91 | 85 | 131 |
Failure type | Through vein | Through vein | Through vein |
Angle of rupture relative to core axis (°) | 19 | 22 | 28 |
τpeak (MPa) | 29.98 | 27.22 | 50.13 |
σn(peak) (MPa) | 17.11 | 9.91 | 43.94 |
Model Component | Input Parameter | Assumed Values | Value from Literature | Data from |
---|---|---|---|---|
Vein | E(v) (GPa) | 240 | [36]: 231.1 {14}; [37]: 235 {10}; [38,39]: 306.5 {2}; [40]: 262.8 {1} (Format: average value {# of specimens}) | |
ν(v) | 0.17 | [38,41] | ||
ϕp(v) (°) | 47 | [41] | ||
ϕr(v) (°) | Equal to ϕp(v) | |||
cp(v) (MPa) | 4.7 | [41] | ||
cr(v) (MPa) | 0.1 | |||
σtp(v) (MPa) | 2 | [41] | ||
σtr(v) (MPa) | 0.1 | |||
Vein–wallrock contact joint element | Kn(j) (MPa/m) | Equal to Ks(j) | Following guidance from [34,35] | |
Ks(j) (MPa/m) | Equal to Kn(j) | Following guidance from [34,35] | ||
ϕp(j) (°) | 37 | Interpreted from [23] | ||
ϕr(j) (°) | Equal to ϕp(j) | |||
cp(j) (MPa) | 18 | Interpreted from [23] | ||
cr(j) (MPa) | 0.1 | |||
σtp(j) (MPa) | 5 | Interpreted from [23] | ||
σtr(j) (MPa) | 0.1 |
Parameter | Minimum | Intermediate Low | 1st Iteration Base Case | Intermediate High | Maximum | |
---|---|---|---|---|---|---|
Vein material | E(v) (GPa) | 100 | 150 | 205 | 240 | 275 |
ν(v) | 0.13 | 0.15 | 0.17 | 0.19 | 0.21 | |
ϕp(v) (°) | 37 | 42 | 47 | 52 | 57 | |
cp(v) (MPa) | 1 | 5 | 10 | 15 | 20 | |
σtp(v) (MPa) | 1 | 2.5 | 5 | 7.5 | 10 | |
Wallrock (CMET) | E(w) (GPa) | 40 | 50 | 58.5 | 70 | 85 |
ν(w) | 0.1 | 0.12 | 0.14 | 0.16 | 0.18 | |
ϕp(w) (°) | 22 | 27 | 32 | 37 | 42 | |
cp(w) (MPa) | 23 | 33 | 43 | 50 | 60 | |
σtp(w) (MPa) | 8 | 11 | 14 | 17 | 21 |
Parameter | 1st Iteration Base Case | 2nd Iteration Base Case | Minimum | Intermediate | Maximum | |
---|---|---|---|---|---|---|
Vein–wallrock contact joint element | Kn (GPa/m) | - | 1500 | 2000 | 2500 | 5000 |
Ks (GPa/m) | - | 1500 | 2000 | 2500 | 5000 | |
σtp(j) (MPa) | - | 2.5 | 1 | 5 | - | |
ϕp(j) (°) | - | 37 | 32 | 42 | 47 | |
cp(j) (MPa) | - | 18 | 15 | 21 | 24 | |
Vein material | E(v) (GPa) | 205 | 125 | 100 | 150 | 175 |
ν(v) | 0.17 | 0.17 | - | - | - | |
ϕp(v) (°) | 47 | 2.5 | 1 | 5 | - | |
cp(v) (MPa) | 10 | 47 | 42 | 52 | - | |
σtp(v) (MPa) | 5 | 20 | 15 | 25 | - | |
Wallrock (CMET) | E(w) (GPa) | 58.5 | 58.5 | - | - | - |
ν(w) | 0.14 | 0.14 | - | - | - | |
ϕp(w) (°) | 32 | 32 | - | - | - | |
cp(w) (MPa) | 43 | 43 | - | - | - | |
σtp(w) (MPa) | 14 | 14 | - | - | - |
Parameter | Minimum | 3rd Iteration Base Case | Maximum | |
---|---|---|---|---|
Vein–wallrock contact joint element | Kn (GPa/m) | 1500 | 1750 | 2000 |
Ks (GPa/m) | 1500 | 1750 | 2000 | |
σtp(j) (MPa) | 2.5 | 5 | 7.5 | |
ϕp(j) (°) | 37 | 42 | 47 | |
cp(j) (MPa) | 17 | 22 | 27 | |
Vein material | E(v) (GPa) | 100 | 125 | 150 |
ν(v) | - | 0.17 | - | |
ϕp(v) (°) | 5 | 7.5 | 10 | |
cp(v) (MPa) | 43 | 48 | 53 | |
σtp(v) (MPa) | 21 | 26 | 31 | |
Wallrock (CMET) | E(w) (GPa) | - | 58.5 | - |
ν(w) | - | 0.14 | - | |
ϕp(w) (°) | - | 32 | - | |
cp(w) (MPa) | - | 43 | - | |
σtp(w) (MPa) | - | 14 | - |
Confining Stress, σ3 (MPa) | Vein Thickness (mm) | Vein Orientation, α (°) | Increments of Vein Orientation between Models (°) |
---|---|---|---|
0 | 4 | 10–80 | 10 |
5 | 4 | 10–80 | 10 |
10 | 1, 4, and 8 | 5–85 | 5 |
15 | 4 | 10–80 | 10 |
Comparison Item | Data from Jaeger and Cook [42] | Numerical UCS/TCS Experiment Results from This Study |
---|---|---|
Discontinuity type | Foliation | Single hydrothermal vein |
Range of discontinuity orientation where failure occurs on discontinuity, i.e., critical angles (α, °) | 10–52 | 5–55 |
Weakest orientation (α, °) | 30 | 20–30 (for different σ3) |
Influence of discontinuity outside critical angles | No influence | Decreases specimen peak strength |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Rudderham, G.A.; Day, J.J. Veined Rock Performance under Uniaxial and Triaxial Compression Using Calibrated Finite Element Numerical Models. Geotechnics 2023, 3, 1219-1250. https://doi.org/10.3390/geotechnics3040067
Rudderham GA, Day JJ. Veined Rock Performance under Uniaxial and Triaxial Compression Using Calibrated Finite Element Numerical Models. Geotechnics. 2023; 3(4):1219-1250. https://doi.org/10.3390/geotechnics3040067
Chicago/Turabian StyleRudderham, Gisèle A., and Jennifer J. Day. 2023. "Veined Rock Performance under Uniaxial and Triaxial Compression Using Calibrated Finite Element Numerical Models" Geotechnics 3, no. 4: 1219-1250. https://doi.org/10.3390/geotechnics3040067
APA StyleRudderham, G. A., & Day, J. J. (2023). Veined Rock Performance under Uniaxial and Triaxial Compression Using Calibrated Finite Element Numerical Models. Geotechnics, 3(4), 1219-1250. https://doi.org/10.3390/geotechnics3040067