Size-Dependent Behaviour of Hard Rock Under Triaxial Loading

Understanding size effects is important to rigorously analyse the behaviour of rocks and rock masses at different scales and for different applications. A number of empirical and numerical studies have included size effects on the uniaxial compressive strength of different rocks, but only few have focussed on the triaxial compressive strength. In this study, several triaxial tests were conducted on granite samples at different confinements (from 0.2 to 15 MPa) and sizes (from 30 to 84 mm in diameter). The most relevant strength parameters were recovered including peak and residual strengths, orientation and frictional response of shear bands. Size effects were observed to be less dominant at high confining pressures for both peak and residual strengths. The resulting data were analysed in the context of size-dependent rock strength theories. For peak strength analyses, the unified size effect law (USEL) and the improved unified size law (IUSEL) were used, in which the double trends, ascending descending behaviours were observed. Size-dependent Hoek–Brown modified failure criteria based on USEL and IUSEL were fit to peak strengths exhibiting a good agreement between the models and the laboratory data. The brittle-ductile transition and the frictional behaviour of shear band under triaxial loading were also analysed. A clear brittle behaviour was particularly observed in samples with large-diameters tested at low confinements. Finally, in regards to frictional behaviour, the shear band angle found to be affected by both confinement and sample size. Triaxial tests were conducted on granite samples at different confinements (0.2 to 15 MPa) and sizes (30 to 84 mm in diameter). Size-effect (ascending-descending double trend) was observed on peak strength results. USEL–IUSEL size-effect laws and Hoek-Brown modified failure criterion were used to model the strength behaviour. Brittle-ductile transition and frictional behaviour related to size-effect were analysed. Triaxial tests were conducted on granite samples at different confinements (0.2 to 15 MPa) and sizes (30 to 84 mm in diameter). Size-effect (ascending-descending double trend) was observed on peak strength results. USEL–IUSEL size-effect laws and Hoek-Brown modified failure criterion were used to model the strength behaviour. Brittle-ductile transition and frictional behaviour related to size-effect were analysed.


Introduction
Study of size effect is important in rock engineering particularly for intact rocks where the laboratory data should be extrapolated to the field settings. Size effect has been proven to be an important factor that can influence the mechanical properties of intact rocks. Several studies have included size effect behaviour of intact rocks under various stress conditions. Majority of these have been performed under uniaxial compression (Mogi 1962;Dhir and Sangha 1973;Baecher and Einstein 1981;Hawkins 1998;Thuro et al. 2001a;Yoshinaka et al. 2008;Darlington et al. 2011;Masoumi et al. 2012Masoumi et al. , 2014Masoumi et al. , 2016Masoumi 2013;Roshan et al. 2017;Quiñones et al. 2017), some under point loading (Broch and Franklin 1972;Bieniawski 1975;Brook 1977Brook , 1980Brook , 1985Wijk 1980;Greminger 1982;Hawkins 1998;Thuro et al. 2001b;Masoumi et al. 2012Masoumi et al. , 2018Forbes et al. 2015) and few under tensile testing (Wijk 1980;Andreev 1991a,b;Butenuth 1997;Thuro et al. 2001b;Çanakcia and Pala 2007;Masoumi et al. 2018). Despite extensive studies in this area, few works included the size effect under triaxial condition, focussing mainly on sedimentary rocks (Singh and Huck 1972;Hunt 1973;Medhurst and Brown 1998;Aubertin et al. 2000;Roshan et al. 2017;Masoumi et al. 2017a) and only one study under limited testing conditions was conducted on hard rock by Walton (2018). This could be due to the complexity associated with the triaxial testing on intact brittle (or quasi-brittle) materials. Such a testing is important as most of the design analyses on rock and rock-like materials are performed based on well-known Mohr-Coulomb (Mohr 1900) or Hoek and Brown (H-B) (1997) criteria where the estimation of cohesion (C) and friction angle (ϕ) as well as H-B parameters are essential. It is noteworthy that the brittleness in materials such as rock is relative to the stiffness of them versus stiffness of the loading frame. Therefore, in this study the referral to "brittle" or "ductile" behaviour is indeed associated with the type of loading frame used in the testing program. Our loading frame can be considered as the type with highest stiffness that are commercially available in the market. The need to extrapolate the results obtained from the laboratory to a larger size has led to the development of three main size effect theories for intact rocks. The first is based on the statistical distribution of micro-cracks within an intact rock following Weibull (1939) statistical function. The second theory was proposed by Bazant (1984) which relies on the fracture energy concept. The third one was suggested by Carpinteri et al. (1995) based on multifractal scaling concept. It is noteworthy that a number of empirical size effect models have been developed implicitly based on Weibull (1939) statistical theory (e.g. Hoek and Brown 1980;Brook 1980, 1985and Darlington et al. 2011. From these size effect models, different size-dependent failure criteria have been proposed such as the one by Hoek and Brown (1997) through including an empirical size effect function into their well-known failure criterion. Later, Medhurst and Brown (1998) assessed the predictability of Hoek and Brown sizedependent failure criterion on a number of coal samples with different sizes and at various confining pressures. A comprehensive triaxial size effect criterion was proposed by Aubertin et al. (2000) in which the size effect parameter was included into MSDP (Multiaxial Misses-Schleiner and Drucker-Prager criterion). Masoumi et al. (2017a) suggested a methodology to include the unified size effect law (USEL) (Masoumi et al. 2016) into Hoek and Brown failure criterion and assessed its performance based on the triaxial data obtained from Gosford sandstone.
The genesis of USEL originated from Hawkins (1998) study who conducted an extensive set of laboratory tests on different sedimentary rocks and demonstrated that the uniaxial compressive strengths (UCS) of these rocks at different sizes follow an ascending and then descending trend. At about the same time, Bazant (1997) developed a fractal fracture size effect model to capture such an ascending and then descending behaviour. Masoumi (2013) and Masoumi et al. (2016) extensively investigated such a phenomenon for a particular sandstone, as well as those investigated by Hawkins (1998) and explained that such a strength ascending and then descending behaviour can be associated with the effect of end surface flaws, which are dominant in the samples with small diameters and also the fractal characteristics of rocks. These two physical theories were then supplemented by USEL proposed by Masoumi et al. (2016). Later, such a phenomenon was endorsed by Quiñones et al. (2017) who also reported such a strength ascending and then descending trend on Blanco Mera granite under uniaxial compression as a first study on hard rock samples. In another study, Masoumi et al. (2017a) demonstrated the applicability of strength ascending and then descending trend to Gosford sandstone under triaxial condition, which all were explained based on the two physical theories highlighted earlier including end surface flaws effect and fractality. Zhai et al. (2020) also demonstrated that the strength ascending and descending behaviour is applicable to weak rocks, including Gambier limestone and a synthetic gypsum material with UCS of less than 10 MPa.
The sample size not only affect the elastic and strength properties of an intact rock, it can also impact the overall mechanical behaviour of rock under loading such as transition from brittle to ductile behaviour and formation of shear band particularly under triaxial condition (Roshan et al. 2017). Therefore, it is relevant to advance towards extensive assessment of size effect models for different rocks, in particular those which have been limitedly examined such as hard rocks. To contribute to this endeavour, a dedicated laboratory program was developed to extend the current knowledge gap on size effect behaviour of hard rock under triaxial condition (Arzúa and Alejano 2013;Quiñones et al. 2017;Alejano et al. 2021). As a result, about 100 triaxial compressive tests were performed on Blanco Mera granite at four different diameters including 30, 38, 54 and 84 mm using a set of specially tailored Hoek cells. The triaxial tests were conducted between 0.2 and 15 MPa confining pressures. The resulting shear band orientations at different sizes and confinements were measured and compared along with the transition from brittle to ductile behaviour due to change in size and confining pressure. The strength ascending and then descending trends of granite samples were examined based on USEL and its improved version (IUSEL) and then, their fitability to Hoek and Brown criterion were assessed. Finally, it was concluded that the size-dependent Hoek-Brown failure criterion can accurately capture the strength behaviours of Blanco Mera granite.

Rock Sample Selection and Preparation
The Blanco Mera granite used in this study was supplied from the province of Lugo in Spain. Arzúa and Alejano (2013), Arzúa et al. (2014) and Alejano et al. (2017) performed different mechanical studies on Blanco Mera granite which is a bright white-coloured granite that consists of 35% plagioclase feldspar, 27% alkali feldspar, 20% quartz, 7% muscovite, 5% biotite and other minerals at below 5% including chlorite and sericite. The grain distribution in this granite, to some extent is scattered as shown in Fig. 1a where the size of quartz crystals ranges from 1 to 6 mm and feldspar crystals from 1 to 3 mm along with various mica minerals.
About 100 core samples were prepared at diameters of 30, 38, 54 (NX) and 84 mm (see Fig. 1b and c). The length to diameter ratio of the samples was about 2 complying with ASTM (2014) and ISRM (2007) recommendations which was also suitable to be used in Hoek cell. A number of Hoek cells at different sizes (see Fig. 1d) were designed  2,2.5,5,7.5,10,12.5 and 15 MPa). To account for variability in data, a number of repetitions was performed at different diameters and confining pressures as listed in Table 1.

Testing Procedure
To perform the triaxial tests, a 200 tonne servo-controlled loading frame was used as shown in Fig. 2a. The testing system consisted of a load frame where the Hoek cell was fitted into its central part, the extensometers which were mounted to the Hoek cells, the hydraulic pumps to control the confining pressure, test controller and a computer for data acquisition. Strain gauges were used inside the Hoek cell which were attached to the samples for axial and radial strains measurement as demonstrated in Fig. 2b. The frame movement was also measured axially using Linear Variable Differential Transducer (LVDT). Water was used as a fluid to supply confining pressure and its volumetric changes were measured to report sample volumetric strain which is useful to study the post-failure response.

Test Results
Extensive laboratory data were extracted from the triaxial experiments including peak and residual strengths, elastic limit (derived from gauge measurements) and shear band angles. A typical representation of the data obtained from the tests is shown in Fig. 3. This figure demonstrates the stress-strain curves obtained from the axial and radial strain gauges, the axial strain recorded by LVDTs and the volumetric strain. The peak strength was the maximum stress recorded during the test and the residual strength was calculated as the value of stress towards the end of the experiment where an increase in the strain does not cause any significant (e.g. less than 1 MPa) change in the stress. The elastic limit was calculated from the linear section of the stress-strain curve, using the data collected by the strain gauges where the linear part lies at about 50% peak strength. The average of the resulting peak and residual strengths, elastic limits and shear band angles at different confinements and diameters are presented in Table 2. The resulting stress-strain curves obtained from the sample with 38 mm diameter at various confinements along with the triaxial tests under 0.2 MPa confining pressures on samples with 30, 38, 54 and 84 mm diameters are presented in Fig. 4.
The shear band angle β was defined as the angle between the shear band and the horizontal line as shown in Fig. 5. It is worth remarking the impossibility of measurement of such an angle for almost all samples tested under 0.2 MPa confining pressure due to complex failure patterns which tend to include a number of tensile cracks with no or little shearing as depicted in Fig. 5. Also, all failed samples with 30 mm diameter are shown in Fig. 6. The resulting failure patterns under low confining pressures are mainly axial splitting with a transition to shear band failure mechanism at high confinements.
The peak strengths obtained from various diameters and confining pressures were plotted and compared in Fig. 7a. It is clear that an increase in the confining pressure led to a rise in the peak strengths, where an ascending and then  descending trend is observable at almost all confining pressures due to change in sample diameter which can be associated with end surface flaws effect as well as fractal characteristics of Blanco Mera granite according to Masoumi et al. (2016). It is believed that at 5 MPa confinement, rock heterogeneity was the main factor that such a trend was not clearly visible compared to other confining pressures. From these data, two trends can be deduced. At low confining pressures (less than 5 MPa), there is a clear size effect behaviour as that reported by Masoumi (2013) and Quiñones et al. (2017) where the strength increases from small sizes up to 54 mm diameter and then decreases at larger diameters. At higher confinements (between 5 and 10 MPa), such an ascending and then descending trend is not as clear as those observed at low confinements highlighting Aubertin et al. (2000) argument where an increase in confining pressure can lead to less size effect in rock materials. Such an argument is even more recognisable for samples tested at 12.5 and 15 MPa confining pressures. This behaviour was also reported by Gerogiannopoulos and Brown (1978) as well as Michelis (1987) who tested marble. The role of pre-existing microcracks in tested granite can be considered as the major factor for controlling the size effect at different confining pressures. It is believed that at low confining pressures the pre-existing micro-cracks control the size effect based on any defined theories such as weakest link by Weibull (1939), fracture energy by Bazant (1984) or multifractality by Carpinteri et al. (1995). However, with an increase in confining pressure, these cracks close and thus less size effect behaviour is expected. Such a finding is quite important from practical viewpoint where the data from the laboratory environment is extrapolated to the field setting such as a deep underground rock engineering project.
A size effect analysis was also conducted on the residual strengths obtained from various confinements and diameters. Figure 7b illustrates the size effect behaviour of the resulting residual strengths where only at low confinements (less than 5 MPa), the ascending and then descending trends similar to that observed for the peak strengths are recognisable. The lack of significant size effect at the residual strengths was also reported by Archambault et al. (1993), Roshan et al. (2017), Walton (2018) and Alejano et al. (2021) who argued that due to disintegration of rock samples after failure, the defined size effect theories cannot be applicable to residual state.

Analytical Study
In this section, various aspects of results relevant to size effects on intact rock are analysed. First the peak strengths at different confinements are analysed in the context of different size effect laws to assess their size dependency in regards to confinement effect. Then, a couple of size-dependent Hoek-Brown failure criteria inspired by the work conducted by Masoumi et al. (2017a) are proposed to examine the reliability of such criteria for Blanco Mera granite.

Size Effects Laws and Their Fitability
Unified size effect law (USEL) (Masoumi et al. 2016) and its improved version (IUSEL) (Masoumi et al. 2017b) were used to perform size-dependent analysis on peak strength Fig. 3 Example of stress-strain curves obtained from the triaxial test on a sample with the diameter of 38 mm tested at 12.5 MPa confinement Table 2 Average of the peak and residual strengths, standard deviation and coefficient of variation for peak and residual strengths, elastic limit and shear band angle obtained from tested sam- 1 3 data. The resulting trends were also compared so that the predictability of each model can be assessed. Masoumi et al. (2016) proposed the unified size effect law (USEL) to predict the strength ascending and then descending behaviour of intact rocks, particularly sedimentary rocks. USEL consists of two parts including the strength descending behaviour which was explained based on the fracture energy theory according to Bazant (1984) who proposed size effect law (SEL) and the strength ascending behaviour based on fractal fracture size effect law (FFSEL) introduced by Bazant (1997). A combination of SEL and FFSEL was defined as USEL so that the strength at each size is estimated as the minimum value obtained from SEL and FFSEL as follows:

Unified Size Effect Law
where, σ 0 is the intrinsic or characteristic strength associated with the ascending part, d is the sample size, λ is a dimensional parameter, d 0 is the maximum sample size that can be referred to as a characteristic size and D f is the fractal dimension related to the fracture surface. It is noteworthy that D f = 1 for non-fractal surfaces and D f ≠ 1 for fractal surfaces. In SEL, B is the dimensional parameter and f t is the intrinsic or characteristic strength.
To calibrate the USEL parameters, first the SEL parameters including Bf t and λd 0 are estimated through the curve fitting process. Hence, a similar approach is followed to obtain the FFSEL parameters being σ 0 and D f . It should be noted that for those sizes in which D f = 1, Bf t is equal to σ 0 , meaning that the corresponding data will not show fractal characteristics and thus, the FFSEL becomes the same as SEL. The intersection diameter of SEL and FFSEL (D i ) happens according to: and its corresponding strength can be estimated using the following equation: (2) Masoumi et al. (2017b) proposed an improved version of USEL (IUSEL) based on fractality concept. The main advantage of IUSEL versus USEL is its ability to predict a realistic strength for a sample with an infinite size. Such an improvement was developed based on two fractal models in the same vein as that for USEL including ascending and descending parts. The descending part was developed by Carpinteri et al. (1999) known as "multifractal scaling law" or MFSL. The ascending part was developed by Masoumi et al. (2017b) following the same approach as that implemented by Bazant (1997) for FFSEL leading to "modified multifractal scaling where, f m is a characteristic strength of the ascending zone, d is the sample size, d f is the fractal dimension, l is a material constant with length unit and f c is a characteristic strength associated with the descending zone. The calibration process for IUSEL is consistent with that for USEL where the parameters of descending part are estimated first, followed by the ascending part. The diameter at the intersection of MFSL and MMFSL can be estimated using below equation:

Improved Unified Size Effect Law
and the strength at the intersection diameter or the maximum strength of Eq. (4) is estimated as follows: (4)

Predictability Assessment of USEL and IUSEL on the Triaxial Data with Negligible Confinement
According to Masoumi et al. (2016) and Zhai et al. (2020) three scenarios can be considered for calibration of USEL and IUSEL to the uniaxial compressive strength data or the triaxial data with negligible confinement such as here (0.2 MPa). First, the characteristic diameter with maximum strength is only included in the ascending part, second, it is only included in the descending part and lastly it is included in both ascending and descending parts. All these three scenarios were modelled for both USEL and IUSEL and the resulting parameters are listed in Tables 3  and 4, respectively. Only the ones with the best fit (maximum R 2 ) which are underlined in Tables 3 and 4 are visualised in Fig. 8. It should be noted that for both USEL and IUSEL, the best fits happened when the maximum strength was included in both ascending and descending parts. From Fig. 8, it is clear that predictions by both USEL and IUSEL match well with the experimental data. The major improvement in IUSEL is the realistic prediction of the strength of a sample with an infinite size which is highlighted in Fig. 9 where the strength prediction of a granite with very large size at 0.2 MPa confining pressure by IUSEL tends to a constant value while USEL tends to zero.  Masoumi et al. (2017a) modified the original and widely used Hoek and Brown (1997) failure criterion for intact rocks by incorporating USEL into it leading to the size-dependent

Size-Dependent Hoek-Brown Failure Criterion
Hoek and Brown failure criterion. In the modified criterion, the uniaxial compressive strength (σ ci ) was replaced by USEL resulting in: where, 1 is the axial stress, 3 is the confining pressure and m is the Hoek and Brown constant that controls the curvature of the criterion. The rest of the parameters were defined earlier when USEL was presented. The improved version of USEL was proposed by Masoumi et al. (2017b) and here the IUSEL is incorporated into Hoek and Brown failure criterion in the same vein as that considered by Masoumi et al. (2017a) leading to: where, all the parameters were already defined. For models calibration, a wide range of data sets at different confinements and diameters are required so that the USEL and IUSEL parameters can be calibrated based on the data obtained under uniaxial compression (or in our case under triaxial condition with negligible confining pressure). Hence, m is calibrated for each sample size tested at various confining pressures and then the average of m values obtained from (7) different diameters is used to visualise the three-dimensional graph of axial stress-confining pressure-sample size. As a result, the USEL and IUSEL parameters (see Tables 3 and 4) that were calibrated earlier for the samples with different diameters tested under 0.2 MPa confinement were used to estimate the remaining parameters of sizedependent Hoek and Brown models as listed in Table 5. It is noteworthy that σ ci in Table 5 for USEL and IUSEL is the strength prediction for each diameter based on USEL and IUSEL and the resulting m i for each diameter obtained from the curve fitting.
The size-dependent Hoek-Brown failure criteria based on USEL and IUSEL are presented in Fig. 10.
The resulting values for the parameters of size-dependent Hoek-Brown failure criteria based on USEL and ISUEL found to be very close and thus, the deriving fits from these parameters depicted on Fig. 10 are overlapped. It is noteworthy that the predictions by both models match well with the experimental data. It is also worth mentioning that based on the modelling analysis, it was concluded that m parameter to high extent seem to be size independent as the resulting m values from different diameters are close to each other. Such a size independent characteristic can be further explored for different rocks based on more laboratory triaxial experiments on rocks with various sizes.

Brittle-Ductile Transition and Frictional Behaviours
In this section, the potential size dependency of the brittleductile transition and the frictional behaviour of Blanco Mera granite due to increase in confining pressure are investigated.

Size Effects on Brittle-Ductile Transition
It has been proven that rocks show a transition from brittle to ductile behaviour due to an increase in the confining stress (Mogi 1966;Brady and Brown 2006). Some studies have linked this behaviour to the closure of non-critical microfractures as the confining pressure is increased (Wawersik and Fairhurst 1970;Elliot and Brown 1985). Others have identified the type of micro-failure mechanism at stake being responsible for such a transition (Griffith 1921;Martin and Chandler 1994;Diederichs 2007) where tensile failure associated with micro-cracks propagation and coalescence tends to produce a more brittle response while the shear failure at typically large sizes and high confinements induce more ductile behaviour. Certainly, both reasons can contribute to the transition from brittle to ductile behaviour. While the effects of confining pressure on brittle to ductile transition have been widely studied for different rock types (Edmond and Paterson 1972;Barton 1976;Shah 1997;Walton et al. 2017;Bruning et al. 2019), very limited studies (e.g. Roshan et al. 2017) have analysed the size dependency of such a transition. According to Roshan et al. (2017), the transition from brittle to ductile behaviour can be defined as a non-dimensional transition index (TI) which is the difference between the mean peak and residual strengths at different confining pressures and is normalized with the mean peak strength as follows:  TI can vary from one to zero where for the rocks with high brittleness it tends to one and for those with high ductile behaviour, it tends to zero or negative values. To calculate the TIs for each size and confining pressure, data presented in Table 2 were used. Figure 11 illustrates the TI indices estimated for each sample size and confining pressure confirming that Blanco Mera granite is a very brittle rock. TI reached to one at low confining pressures while it exhibited tendency towards more ductile behaviour as the confining pressure increases.
The size dependency of such a transition is to some extent inconclusive particularly at high confining pressures where the resulting trends mostly overlapped. This can be potentially explained through size independency of the residual strength, which was introduced earlier, and will be discussed in its appropriate context further below.
(9) TI = Mean peak strength − Mean residual strength Mean peak strength Byerlee (1967) suggested that natural faults formed under progressive initial shear failure as many found in nature, showing a behaviour similar to that observed under triaxial experiment of intact rocks. Faults and shear bands typically initiate before reaching the peak strength (Ord et al. 2007), particularly the coalescence process, then they tend to grow to form the macro shear plane during the stress softening part of the stress-strain curve. During the stress softening, rock disintegrates until reaching the full shear failure at the residual state or so-called residual strength (Roshan et al. 2017) following a large shear displacement along the formed shear band. Figure 12a shows the mean shear band angles versus confining pressure at different sizes of Blanco Mera granite extracted from Table 2. As previously noticed, it was not possible to measure the resulting shear bands angle at 0.2 MPa due to axial splitting nature of the failure mode at no or very small confinements. At larger confinements, it is evident that an increase in the confining pressure can lead to decrease in the shear band angle for all sample sizes. Maurer (1965) suggested that this reduction in the shear band angle is mainly associated with the increase in the ductility which was earlier investigated for Blanco Mera granite through TI. Also, Fig. 12b illustrates the variation of shear band angle versus sample sizes at different confining pressures highlighting the size dependency of shear band for Blanco Mera granite. At the 84 mm diameter, the highest values of shear band angles at various confining pressures are observable followed by 54 mm diameter samples and then 38 and 30 mm diameters. The resulting trends of shear band align with those reported by Zhang and Li (2019) for a Norwegian granite and Alejano et al. (2010) for slightly weathered Spanish granites.

Size Effect on Friction Coefficient
The friction coefficient (FC) is defined as the relationship between the shear (τ) and normal (σ n ) stresses acting on the shear plane. The friction coefficient can be therefore expressed as: At high pressures, it is found that the shear and normal stresses during the failure could be approximated by a linear law according to: where A and B are constants and thus, FC could be defined as the slope of linear law presented in Eq. 11 or A. Size effects on friction coefficient have been studied by several researchers including Bandis (1980), Bandis et al. (1981), Schellart (2000), Carpinteri and Paggi (2005), Carpenter et al. (2011Carpenter et al. ( , 2015 and Roshan et al. (2017). These works focussed on assessing the size-dependent behaviour of fractures formation and evolution mainly in structural geology. On the other hand, experimental studies have suggested that the behaviour of natural fractures to some extent is different to that of an intact rock where a fracture starts to occur, i.e. the formation of a shear band (Byerlee 1967). Therefore, the early stages of fault formation tend to be better represented by the formation of shear band in intact rock under triaxial loading where the normal (σ n ) and shear (τ) stresses on the shear plane can be expressed as a function of principal stresses (σ 1 , σ 3 ) and the shear band orientation (β) according to: Estimation of normal and shear stresses for peak and residual strengths are important due to the shear band formation that is initiated before reaching the peak strength and then grows beyond the peak stress all the way to the residual strength (Bèsuelle et al. 2000). As a result, the friction coefficients (FC) can be estimated at both peak and residual strengths which are shown in Fig. 13a and b. It is noteworthy that the average calculated shear and normal stresses for all mean peak and residual strengths were used to develop the trends of FCs at different sizes and confinements. Figure 14 demonstrates the values of the friction coefficient for peak and residual strengths in relation to sample size and based on the line slopes obtained from Fig. 13a and b. It is clear that the peak and residual friction coefficients . 13 a Fitting of shear stress-normal stress Mohr-Coulomb lines for peak strengths to estimate the peak friction coefficients of every sample size group. b Fits of shear stress-normal stress Mohr-Coulomb lines for residual strengths to estimate the residual friction coefficients of every sample size group Fig. 14 Graph presenting peak and residual friction coefficients in relation to sample size. The two branches ascending descending trend is observable both for peak and residual cases tend to follow the similar trends as those observed for strengths with an ascending and then descending zones. They are therefore similar to those previously referred to by USEL and IUSEL trends, meaning the value of the peak and residual friction coefficients increase for small diameters until reaching a maximum value at around 40-50 mm and then decrease at larger diameters.

Conclusions
Extensive triaxial tests were conducted on Blanco Mera granite with samples ranging between 30 and 84 mm and the confinements changing from 0.2 to 15 MPa. The triaxial results at negligible confining pressure (0.2 MPa) were analysed based on unified size effect law (USEL) and its improved version (IUSEL). Also, the size-dependent Hoek and Brown failure criteria based on USEL and IUSEL were used to analyse the triaxial data at different sizes. Lastly, the size dependency of the brittle-ductile transition and the frictional behaviour of shear bands were investigated extensively. The analysis of peak strengths at different confining pressures revealed ascending and then descending size effect behaviour. Such a behaviour is found to be more pronounced at low confining pressures compared to high confinements. It is believed that the closure of micro-cracks at high confining pressures is responsible for less size effect as argued by Aubertin et al. (2000). Similar trends were observed for residual strengths whereby an increase in confining pressure less size effect was recognisable. The USEL and IUSEL models were calibrated to predict the peak strengths resulting in a good agreement between the models and the experimental data. The major advantage of IUSEL over USEL is a realist strength prediction for a sample with an infinite size which was demonstrated for Blanco Mera granite. The sizedependent Hoek and Brown failure criteria were calibrated and their suitable fitability to the triaxial experimental data was highlighted. Also, based on the modelling analysis, it was suggested m in Hoek and Brown failure criterion is size independent, however, additional laboratory experiments particularly on other rock types are recommended to further endorse this suggestion.
The size dependency of brittle-ductile transition in Blanco Mera granite was investigated through a transition index (TI). Such a dimensionless index can quantify the brittle-ductile transition based on peak and residual strengths values. The tested granite, as expected, showed a remarkable brittle behaviour at low confinements and more importantly at large samples. Also, the size dependency of the shear band angle was examined and the resulting trends followed similar size effect trends as those observed for peak strengths at different confining pressures. Finally, the friction coefficients of the resulting shear bands were quantified leading to substantial size-dependent behaviour. Author Contributions MAG-F has developed the general structure of the paper, created the graphic art and did the original writing. LRA designed, found the funding, hired personnel and oversaw the research project and carried out the various general review of the document. XE-V has been in charge of lab work, HM has helped in the implementation of the models used, analysis of results and carried out the general review of the document.
Funding Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. The Spanish Ministry of Science, Innovation & Universities has fund this study as part of the project awarded under Contract Reference No. RTI2018-093563-B-I00, PRE2019-087397, partially financed by means of ERDF funds from the EU.

Availability of Data and Material Not applicable.
Code Availability Not applicable.

Conflict of Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.