Open Access Library Journal
Vol.03 No.12(2016), Article ID:72665,20 pages
10.4236/oalib.1103214

Evaluating the Triaxial Characteristics of Ultamafic Rocks from Central Greece Using the Physical, Dynamic and Mechanical Properties

Konstantinos Diamantis1*, George Exarhakos2, George Migiros1, Efthimios Gartzos1

1Department of Sciences, Division of Geological Science and Atmospheric Environment, Laboratory of Mineralogy-Geology, Agricultural University of Athens, Athens, Greece

2Department of Civil Engineering, Piraeus University of Applied Sciences, Athens, Greece

Copyright © 2016 by authors and Open Access Library Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: November 8, 2016; Accepted: November 21, 2016; Published: December 9, 2016

ABSTRACT

Due to the fact that the estimation of the triaxial parameters is very important in any geotechnical project, many attempts have become with the aim of their determination. In spite of its relative simplicity, triaxial tests are time consuming, expensive, and require a large number of well prepared (regularly shaped) rock specimens and suitable testing procedures. For this reason, they cannot be easily determined in rocks which are thinly bedded, highly fractured, weak, and so indirect, fast, practical, easy and economical ways should be used such as the correlation of them with some other characteristics which are calculated easily. The aim of this study is to apply correlation analysis to investigate the relationships between physical, dynamic and mechanical characteristics and triaxial parameters for ultramafic rocks. Thus, sixteen samples, taken from the western part of Othrysmt (ten samples) and the Kallidromomt (six samples, central Greece), were tested and the relations among the properties were described by simple regression analyses. The study reveals strong negative, logarithmic correlations between the effective porosity and triaxial parameters (c, φ). Positive, linear relationships are also indicated between wave velocities, apparent cohesion and friction angle, while the increase of serpentinization percentage causes decrease of the c, φ. Both dry unit weight and Schmidt Hammer Values are logarithmically affected with c, φ. The point load strength index and apparent cohesion and friction angle are strongly correlated by both logarithmic and linear functions, but logarithmic trends present higher determination coefficients.

Subject Areas:

Civil Engineering

Keywords:

Triaxial Tests, Apparent Cohesion, Friction Angle, Ultramafic, Correlation, Greece

1. Introduction

The construction of geotechnical engineering projects like foundations, rock slopes and underground structures is based on the estimation of the rock mass properties. The determination of the triaxial parameters is very important because the rock-mass properties are estimated through them. Despite of its relative simplicity, triaxial tests are time consuming, expensive, and require a large number of well prepared (regularly shaped) rock specimens and suitabletesting procedures. Some rocks are thinly bedded, highly fractured, weak, present alterations such as ultramafic rocks with the result they are usually not suitable for preparing specimens and the determination of triaxial parameters is very difficult. On the other hand the fact that ultramafic rocks are met in many areas all over the world and especially south-eastern Europe has as a result many engineering geology projects to be constructed on/in them [1] . For the above-mentioned reasons, indirect, reliable, quick, sample and economical ways are necessary to estimate the triaxial parameters for the ultramafic rocks.

Ultramafic rocks are members of ophiolitic suite rocks which represent remnants of the Earth’s oceanic crust and upper mantle. Basic, hypabyssal and extrusive rocks cover ultramafic rocks. This succession is idealized and in most cases some members may be absent. This research focuses on ultramafic rocks. They include particular types of geological formation with both petrographic variety (harzburgites, lherzolites, plagioclastic peridotites, dunites, etc.) and structural complexity due to tectonic deformation and alteration. This ocean-floor metamorphism (serpentinization) lead to a modification of their petrographic characteristics (serpentinized varieties of them). Thus, they are from massive strong to weak rocks (peridotites and serpentinites, [1] - [8] ).

Within the framework of the present study, sixteen ultramafic rocks (six peridotites and ten serpentinites) were taken from the Kallidromomt and the western part of Othrysmt (Central Greece), the dry unit weight (γd), the effective porosity (ne), the compressional wave velocity (Vp), the shear wave velocity (Vs), the Schmidt Hammer Values (SHV), the point load strength index (Is50) were determined and presented. The Uniaxial Compressive Strength (σci) and the material constant (mi) were estimated by the Hoek and Brown failure criterion and the triaxial characteristics (apparent cohesion, c and friction angle φ) were presented.

The main objective of this research is was to determine the predictability of apparent cohesion (c) and friction angle (φ) of ultramafic rocks (serpentinites and peridotites) with a simple, fast, practical and economical way through the physical, dynamic and mechanical properties at the preliminary site investigation stage. The results were statistically analyzed using the method of least-squares regression. The relationships among these parameters were described by the best fit equations and the highest correlation coefficient, in each relationship, was also determined. The same samples were subjected to petrographic studies with the aim of describing the main mineralogical composition and the serpentinization percentage of ultramafic rocks.

2. Location and Geology of the Study Area

The study area is to be found in the central part of Greece (Figure 1(a)) and specifically

Figure 1. (a) Study areas and ophiolitic complexes in Greece; (b) Geologigal map and sampling points of the study area at Othrys Mt; (c) Geologigal map and sampling points of the study area at Kallidromo Mt.

at western part of Othrys Mt. (near the villages of Moschokarya, Mega Isoma, LoutraKaitsas and Domokos, Figure 1(b)) and at the eastern-southeastern side of Kallidromo Mt (between the villages of Kallidromo, Reginio and Modio, Figure 1(c)). It belongs to the Sub-Pelagonian geotectonical unit, which is part of the Internal Hellenides and is mainly comprised of Alpine formations ( [9] [10] [11] [12] ). They can be distinguished into:

1) A carbonate sequence of Triassic-Jurassic age which constitutes the basement of the area.

2) A tectonic nappe, mainly ophiolitic (volcano-sedimentary formations, basaltic lavas, basic rocks and ultramafic masses, [11] ).

3) An unconformable sequence of Cretaceous limestones which passes upward to flysch.

The present study research focuses on ophiolitic sequence and mainly on ultramafic masses. They are mainly represented by the peridotites whose degree of serpentinization varies.

3. Methods and Results

3.1. Sample Preparation and Test Procedures

The data pertains to ultramafic rocks sampled from sixteen sites in the study areas (eleven from Othrys and six from Kallidromo, Figure 1). The samples were taken from the surface and cylindrical rock specimens for laboratory testing were prepared according [13] requirements. To reduce the uncertainty regarding the influence of the sample size on the measured properties and especially on strength [14] [15] , their diameters range between 53 and 55 mm, the length to diameter ratio was between 2.0 and 2.5 and the edges of the specimens were cut parallel and smooth [13] . After a macroscopical inspection, seventy-four isotropic, homogeneous, unweathered (or slightly weathered) and free of visible joints ultramafic specimens were considered. The dry unit weight and the effective porosity were determined in accordance with [16] . The sound velocities (Vp, Vs) were calculated from the travel time with the application of ultrasonic compression wave pulses to the samples in accordance with test designations [17] . Moreover, point load strength index values were determined according to [18] requirements (in sixteen specimens, ten serpentinites and six peridotites) and a Type L Schmidt Hammer was used for the determination of Schmidt Hammer Values. Ten points on each specimen were taken and the mean values were calculated. The Schmidt hammer tests were carried out according to [16] . In addition, the triaxial tests were determined using a Hoek cell according to [19] . Three or four different confining pressures (σ3) applied during the triaxial tests (fifty-eight specimens). Range of confining pressures of 0 < σ3< 0.5σci, proposed by [14] , was used. The specimens were first subjected to the required confining pressure and then the axial load was applied until the specimen failed.

Thin sections, from samples of the sixteen areas, were prepared and examined under the polarizing microscope, at the Mineralogy-Geology Laboratory of the Agricultural University of Athens, with the aim of describing the main mineralogical composition and the serpentinization percentage of ultramafic rocks. Fractures created by the tests, did not follow internal discontinuities and were always fresh. The tests were carried out in dry conditions for a better relation of the results.

3.2. Petrographic Properties of Ultramafic Rocks

The ultramafic rocks include rock types (peridotites) with large variability (harzburgites, lherzolites, plagioclasticperidotites, dunites, etc.), which due to their tectonic deformation and a low grade metamorphism (serpentinization, [20] ) lead to a modification of their petrographic characteristics (serpentinites). Serpentinisation is the transformation of ferromagnesian minerals, olivine in particular, to serpentine, which is a lattice mineral of either fibrous or laminar form. This alteration allows the understanding of the dramatic changes in the behavior of ultramafic rocks, reduces their strength and generally the physical, dynamic and mechanical characteristics [1] - [8] .

In this paper, sixteen thin sections were prepared and examined under a polarizing microscope. The studied rock samples are unaltered or slightly serpentinized (serpentinization < 30% by volume) ultramafic rock types, characterized as peridotites and partially to completely serpentinizedperidotites (serpentinization > 70% by volume) called as serpentinites. The investigated peridotites mainly comprised of olivines, ortho and clino pyroxenes which are characterised as parent-primary minerals (Pm > 74% by volume, Table 1).

When the peridotites are influenced by the serpentinization, they transformed into Serpentinites. Serpentinites are mainly composed of secondary minerals (70% - 92%, by volume, Table 1). Serpentine is the major mineral phase (62% - 84% by volume, Table1, [1] ). Also, some other secondary minerals are chlorite, talc and actinolite comprising up to 13%. Thus, the percentage of the remaining parent rock minerals (olivine, pyroxene), is less than 30%. The mean values of serpentinization for peridotites and serpentinites are 9.17% and 79.90% respectively, while the standard deviations are 8.75 and 7.02 respectively (Table 2).

3.3. Physical, Dynamic and Mechanical Properties

The dry unit weight (γd) and effective porosity (ne) were determinedfor 6 peridotites and 10 serpentinites samples and their values are listed in Table 3. The γd and ne for peridotites range from 31.67 to 33.17 KN/m3 and from 0.07 to 0.19 respectively. On the other hand, the values of γd and ne for serpentinites fluctuate between 24.95 and 26.79 KN/m3 and from 0.39 to 3.59 respectively. The serpentinites present lower dry unit weight values and higher effective porosity than peridotites because they are mainly composed of secondary minerals (as it is afore-mentioned) which have lower specific gravity and more voids.

The compressional wave velocity (Vp) and shear wave velocity (Vs), which were determined as above-mentioned in chapter 3.1, are given in Table 3 as well. The values of Vp for peridotites and serpentinites fluctuate between 7412 and 7991 m/sec and from 4955 to 5645 m/sec respectively, while the mean values are 7781 m/sec (S.D. 211) and 5361

Table 1. Primary minerals, secondary minerals and serpentinization percentage for ultramafic rocks.

KP: Peridotite of Kallidromo, OP: Peridotite of Othrys, KS: Serpentinite of Kallidromo, OS: Serpentinite of Othrys.

Table 2. Statistical analysis of the serpentinization percentage for ultramafic rocks.

m/sec (S.D. 225) respectively. As far as Vs values are concerned, they vary between 4171 and 4539 m/sec for peridotites and from 2464 to 3044 m/sec for serpentinites. The mean values for serpentinites and peridotites are 2795 and 4392 m/sec respectively. According to [2] [3] and [5] , both unit weight and velocities decrease with increasing of serpentinization.

Moreover, the Schmidt Hammer Values (SHV) range between 61.25 and 64.69 for peridotites (Table 3), present mean value 63.75 and their standard deviation is 1.24. On

Table 3. Physical, dynamic and mechanical property values for ultramafic rocks.

the other hand, in serpentinites the SHV vary from 49.52 to 54.49, the mean value is 52.17 and their standard deviation is 1.61 (Table 4).

In this research, only diametrical point load tests were carried out on the samples. The point load strength index (Is50) (referred to a standard size 50 mm) values are listed in Table 3. The point load strength index values range from 2.01 to 4.91 MPa in serpentinites and between 7.27 and 10.63 MPa in peridotites. The mean values are 3.41 (S.D. 0.91) and 9.54 MPa (S.D. 1.42) respectively (Table 4).

In this research, the triaxial tests were carried out in fifty-eight specimens (three confining pressures at least in each sample ranging from 3 to 18 MPa, Table 5) and the material constant (mi) and the σci were estimated by the Hoek and Brown failure criterion (Table 5). For this criterion [14] [21] [22] [23] [24] , the maximum and minimum total stresses at failure (σ1 and σ3) for intact rock are correlated by the following equation:

(1)

where,

σ1, σ3, the maximum and minimum (confining) total stresses at failure for intact rock,

Table 4. Statistical analysis of the physical, dynamic and mechanical properties for ultramafic rocks.

σci, the uniaxial compressive strength of intact rock,

mi, the material constant which depends on the properties of intact rock.

Finally the apparent cohesion (c) and the friction angle (φ) for intact rock were calculated using triaxial compressive strength test data (Table 5). The confining pressures were plotted against axial stresses with the σ3 as abscissa and the σ1 as ordinates (Figure 2 and Figure 3). The best linear lines passing through the points are fitted using the least squares technique. The regression functions of this line for peridotites and serpentinites are respectively:

(2)

and

(3)

Furthermore, c and φ were calculated using the following equations:

Table 5. The maximum and minimum (confining) total stresses, the material constant, the friction angle, the apparent cohesion and the estimated uniaxial compressive Strength values for ultramafic rocks.

Figure 2. The confining pressures (σ3) against axial stresses (σ1) for Peridotites.

Figure 3. The confining pressures (σ3) against axial stresses (σ1) for Serpentinites.

(4)

(5)

where

a, the gradient of the Equations (2) and (3).

4. Correlation Analysis

Least squares regression analysis was applied in order to describe the relationships among triaxial characteristics (c, φ), physical, dynamic, mechanical properties and the petrographic data of the ultramafic rocks. The equation of the best-fit line and the determination coefficient (R2) were determined for each regression. Techniques from Excel 2003 software (Analysis ToolPak program) was used to process the data.

Correlation between Triaxial Parameters and Physical, Dynamic and Mechanical Properties

In this study, an attempt to correlate c and φ with dry unit weight is respectively presented in Figure 4 and Figure 5. The logarithmic equation seems to fit better (R2 = 0.84) the dry unit weight with the apparent cohesion while the relationship between φ and γd is better described by the linear equation (R2 = 0.92). The regression lines representing the best fit between effective porosity and c, φ are logarithmic (Figure 6 and Figure 7). The apparent cohesion and the friction angle using the effective porosity values can be expressed by the empirical equations:

(6)

(7)

Moreover, the plots of the apparent cohesion values as a function of sound velocities (Vp, Vs) values are shown in Figure 8. The best fit relationships are found to be linear and the determination coefficients are 0.83 and 0.84 respectively. The correlations between the friction angle and the wave velocities were also determined. A positive

Figure 4. Relationship between apparent cohesion and dry unit weight for ultramafic rocks.

Figure 5. Correlation between friction angle and dry unit weight for ultramafic rocks.

relationship was found between the φ and the Vp, Vs (Figure 9) indicating that the friction angle increase with increasing wave velocities. The logarithmic model gives high R2 values, but better correlation is described by the linear trends (Figure 9).

As it is illustrated in Figure 10, Figure 11 the friction angle and apparent cohesion present visible correlations with the degree of serpentinization (β). The c, φ decrease linearly with the increase of β. The relationship between φ and β presents higher

Figure 6. Relation between apparent cohesion and effective porosity for ultramafic rocks.

Figure 7. Logarithmic function between friction angle and effective porosity for ultramafic rocks.

determination coefficient (R2 = 0.93) than that between c and β (R2 = 0.86).

In addition, Figure 12 presents the plot of the apparent cohesion against Schmidt Hammer Values (SHV). The best fit equation is logarithmic and the relation between them presents high determination coefficient (R2 = 0.88). The plots of the friction angle as a function of the SHV are shown in Figure 13. The best fit equation is logarithmic with high determination coefficient (R2 = 0.95).

As shown in Table 6, linear and logarithmic functions can describe the relationships between c and Is50 and between φ and Is50. The logarithmic equations seem to fit better

Figure 8. Relationship of wave velocities with apparent cohesion for ultramafic rocks.

Figure 9. Correlation between wave velocities and friction angle for ultramafic rocks.

the above-mentioned properties than the linear equations, exhibiting R-square values 0.87 and 0.93 respectively. The logarithmic functions are illustrated in Figure 14 and Figure 15.

5. Conclusions

The estimation of the triaxial parameters is considered to be the most important component in any engineering geology project because the rock-mass properties are calculated through them. Despite of their simple, fast and easy determination, they require a

Figure 10. Relationship between degree of serpentinization and apparent cohesion for ultramafic rocks.

Figure 11. Correlation between degree of serpentinization and friction angle for ultramafic rocks.

large number of well prepared (regularly shaped) rock specimens. But ultramafic rocks (especially serpentinites) are usually not suitable for preparing specimens. Thus, the immediate determination of the triaxial characteristics is usually difficult for these rocks. For this reason, this study mainly attempts to develop empirical equations between triaxial parameters and physico-mechanical properties.

Figure 12. Relationship of Schmidt hammer values with apparent cohesion for ultramafic rocks.

Figure 13. Variation of Schmidt hammer values against friction angle for ultramafic rocks.

In this paper ultramafic samples, taken from sixteen sites of central Greece (western part of Othrysmt and the Kallidromomt), were tested in laboratory and the triaxial properties were predicted through the wave velocities by simple regression analysis.

The research demonstrates that the friction angle (φ) exhibits strong linear correlation with the dry unit weight (γd) (R2 = 0.92) while the logarithmic equation seems to correlate better (R2 = 0.84) the dry unit weight with the apparent cohesion (c).

Significant negative logarithmic relationships exist between the triaxial parameters

Figure 14. Correlation between point load strength index and friction angle for ultramafic rocks.

Figure 15. Relationship between point load strength index and apparent cohesion for ultramafic rocks.

and the effective porosity (ne). The relation between φ and ne presents higher determination coefficient (R2 = 0.96) than that between c and ne (R2 = 0.91).

The wave velocities are positively correlated with c and φ and provide high determination coefficient. In particular, φ presents higher correlation with Vp, Vs (R2 = 0.91, R2 = 0.92 respectively) than c (R2 = 0.83, R2 = 0.84 respectively).

The effect of alteration on the physico-mechanical properties of ultramafic rocks can be characterized quantitatively by the serpentinization percentage (β). Strong inverse linear

Table 6. Regression equations and determination coefficients (R2) for ultramafic rocks.

relationships exist between the β and c, φ. The relationship between φ and β presents higher determination coefficient (R2 = 0.93) than that between c and β (R2 = 0.86).

The relations between c, φ and Schmidt hammer values (SHV) are also expressed by logarithmic functions. The SHV shows higher correlation (R2 = 0.95) with the friction angle than that with apparent cohesion (R2 = 0.88).

Both the c and φ exhibit negative relationships with the point load strength index (Is50). Τhe best fit trends are logarithmic and φ presents higher correlation with Is50 (R2 = 0.93) than that c with Is50 (R2 = 0.87).

Some observed deviations may be due to the petrographic variety, the different serpentinization degree, the structural complexity and the internally imprinted tectonic deformation of studied rocks.

All empirical methods evaluated in this study can be used for an assessment of the triaxial characteristics of ultramafic rocks, especially for serpentinites whose preparation in standard size cores, is quite difficult. However, it is commonly known that the prediction equations derived by different researches are dependent on rock types, quality, test conditions and different tectonic settings.

Acknowledgements

This study was funded by the State Scholarship Foundation of Greece (I.K.Y). The authors would also like to express their thanks to the Public Works Central Laboratory of Greece (KEDE) and to Ass. Professor Anastasios Tsagalidis for his help in petrography.

Cite this paper

Diamantis, K., Exarhakos, G., Migiros, G. and Gartzos, E. (2016) Evaluating the Triaxial Charac- teristics of Ultamafic Rocks from Central Greece Using the Physical, Dynamic and Mechanical Properties. Open Access Library Journal, 3: e3214. http://dx.doi.org/10.4236/oalib.1103214

References

  1. 1. Diamantis, K. (2010) Engineering Geological Properties of the Ultrabasic Rocks in Othrys and Kallidromo Mountains (Central Greece). PhD Thesis, Agricultural University of Athens, Athens, 386.

  2. 2. Christensen, N.I. (1966) Shear-Wave Vetocities in Metamorphic Rocks at Pressures to 10 Kbar. Journal of Geophysical Research, 71, 3549-3556.
    https://doi.org/10.1029/JZ071i014p03549

  3. 3. Christensen, N.I. (2004) Serpentinites, Peridotites and Seismology. International Geology Review, 46, 795-816.
    https://doi.org/10.2747/0020-6814.46.9.795

  4. 4. Ramana, Y.V., Gogte, B.S. and Sarma, K. (1986) Physical Properties of Indus Ophiolites from Kashmir Himalaya. Physics of the Earth and Planetary Interiors, 43, 104-122.
    https://doi.org/10.1016/0031-9201(86)90079-8

  5. 5. Escartin, J., Hirth, G. and Evans, B. (2001) Strength of Slightly Serpentinized Peridotites: Implications for the Tectonics of Oceanic Lithosphere. Geological Society of America, 29, 1023-1026.
    https://doi.org/10.1130/0091-7613(2001)029<1023:sosspi>2.0.co;2

  6. 6. Marinos, P., Hoek, E. and Marinos, V. (2006) Variability of the Engineering Properties of Rock Masses Quantified by the Geological Strength Index: The Case of Ophiolites with Special Emphasis on Tunnelling. Bulletin of Engineering Geology and the Environment, 65, 129-142.
    https://doi.org/10.1007/s10064-005-0018-x

  7. 7. Ozsoy, E.A., Yilmaz, G. and Arman, H. (2010) Physical, Mechanical and Mineralogical Properties of Ophiolitic Rocks at the Yakakayi Dam Site, Eskisehir, Turkey. Scientific Research and Essays, 5, 2579-2587.

  8. 8. Diamantis, K., Gartzos, E. and Migiros, G. (2014) Influence of Petrographic Characteristics on Physico-Mechanical Properties of Ultrabasic Rocks from Central Greece. Bulletin of Engineering Geology and the Environment, 73, 1273-1292.
    https://doi.org/10.1007/s10064-014-0584-x

  9. 9. Marinos, G. (1974) Geology of Orthrys and Issues on Its Ophiolites. Ann Geol d Pays Jieneiiiciue.st26, University of Athens, Athens, 118-148.

  10. 10. Mountrakis, D., Sapountzis, E., Kilias, A., Eleftheriadis, G. and Christofides, G. (1983) Paleogeographic Conditions in the Western Pelagonian Margin in Greece during the Initial Rifting of the Continental Area. Canadian Journal of Earth Sciences, 20, 1673-1681.
    https://doi.org/10.1139/e83-158

  11. 11. Katsikatsos, G., Migiros, G., Triantaphyllis, M. and Mettos, A. (1986) Geological Structure of Internal Hellenides (E. Thessaly-SW Macedonia. Euboea-Attica-Northern Cyclades Islands and Lesvos). IGME, Geol. and Geoph. Res., Special Issue, 191-212.

  12. 12. Migiros, G. (1990) The Lithostratigraphic-Tectonic Structure of Othris (Central Greece). Bulletin of the Geological Society of Greece, 26, 107-120.

  13. 13. ASTM (2001) Standard Practices for Preparing Rock Core Specimens and Determining Dimensional and shape Tolerances. American Society for Testing and Materials, D4543.

  14. 14. Hoek, E. and Brown, E.T. (1980) Underground Excavations in Rock. Inst. Min. 537 Metall, London.

  15. 15. Hawkins, A.B. (1998) Aspects of Rock Strength. Bulletin of Engineering Geology and the Environment, 57, 17-30.
    https://doi.org/10.1007/s100640050017

  16. 16. ISRM (2007) The Complete ISRM Suggested Methods for Rock Characterization, Testing and Monitoring: 1974-2006. Suggested Methods Prepared by the Commission on Testing Methods, International Society for Rock Mechanics, Compilation Arranged by the ISRM Turkish National Group Ankara, Turkey, 628 p.

  17. 17. ASTM (1983) Test Methods for Ultra Violet Velocities Determination. American Society for Testing and Materials, D2845.

  18. 18. ASTM (2005) Standard Test Method for the Determination of the Point Load Strength Index of Rock. American Society for Testing and Materials, D5731.

  19. 19. ISRM (1981) Rock Characterization Testing and Monitoring. Brown, E., Ed., Pergamon Press, Oxford, 211 p.

  20. 20. Foucault, A. and Rault, J.F. (1995) Dictionnaire de geologie. Masson, Paris.

  21. 21. Hoek, E. and Brown, E.T. (1988) The Hoek-Brown Failure Criterion a 1988 Update. Proceedings of the 15th Canadian Rock Mechanics Symposium, Toronto, 3-4 October 1988, 31-38.

  22. 22. Hoek, E., Wood, D. and Shah, S. (1992) A Modified Hoek-Brown Criterion for Jointed Rock Masses. Rock Characterization: ISRM Symposium, Chester, 14-17 September 1992, 209-213.

  23. 23. Hoek, E., Kaiser, P.K. and Bawden, W.F. (1995) Support of Underground Excavations in Hard Rock. AA Balkema, Rotterdam.

  24. 24. Hoek, E., Carranza-Torres, C.T. and Corkum, B. (2002) Hoek-Brown Failure Criterion— 2002 Edition. Proceedings of the 5th North American Rock Mechanics Symposium, Toronto, 7-10 July 2002, 267-273.