Equivalent continuum analyses of jointed rockmass: Some case studies
|
|
- Bertram Chester Ross
- 5 years ago
- Views:
Transcription
1 Volume 5, Number, March 9, pp.9-5 [INVITED PAPER] Equivalent continuum analyses of jointed rockmass: Some case studies T. G. SITHARAM Professor, Department of Civil Engineering, Indian Institute of Science, Bangalore-56, India. Received 6 8; accepted 9 ABSTRACT The paper presents summary of the work carried out by the author on the modelling of jointed rock mass with some field applications. The equivalent continuum analyses approach presented attempts to use statistical relations, which are simple and obtained after analyzing a large data from the literature on laboratory test results of jointed rock masses. Systematic investigations were done including laboratory experiments to develop the methodologies to determine the equivalent material properties of rock mass and their stress-strain behaviour, using a hyperbolic approach. Present study covers the development of equivalent continuum model for rock mass, implementation of the model in FLACD for -dimensional applications and subsequently verification leading to real field application involving jointed rocks. The model was rigorously validated by simulating jointed rock specimens. Element tests were conducted for both uniaxial and triaxial cases and then compared with the respective experimental results. The numerical test program includes laboratory tested cylindrical rock specimens of different rock types, from plaster of Paris representing soft rock to granite representing very hard rock. The results of the equivalent continuum modelling were also compared with explicit modelling results where joints were incorporated in the model as interfaces. Several case studies have been presented with the developed model. Keywords: Jointed rock mass, equivalent continuum analyses, modeling, field applications, explicit modeling of joints, power house caverns, large slopes. BACKGROUND Rock is distinguished from other engineering materials by the presence of inherent discontinuities such as joints, bedding planes and faults that control its behavior. Hence, the prediction of the response of rocks and rock masses derives largely from their discontinuous and variable nature. Reliable characterization of the strength and deformation behavior of jointed rocks is important for the safe and economic design of civil and mining structures such as arch dams, bridge piers, tunnels, slopes and large caverns. Realistic evaluation of the shear strength and deformation characteristics presents formidable theoretical and experimental difficulties due to the complex geometry and behavior of jointed rock. A numerical approach to treating the rock mass with equivalent material properties for obtaining the overall response has been advocated in recent years (Sitharam, 7; Sitharam et al., ; Sitharam and Latha, ). Several numerical methods have been developed by various researchers to model jointed rock masses using various techniques. Singh (97) has presented continuum characterization methods for jointed rock masses and expressions were presented to estimate the elastic moduli of the equivalent continuum anisotropic rock mass. Zienkiewicz et al. (977) has used the equivalent continuum approach, referred to as a multi-laminate model to simulate a discontinuous rock mass. Gerrard (98) has also used an equivalent continuum approach by expressing the compliance of an element as the sum of the compliances of the intact rock and that of the individual joint sets. Cai and Horii (99) proposed a constitutive model that presents the effects of density, orientation and connectivity of joints, as well as the property of the joints themselves. The constitutive equation is derived from the relation between the average stress and strain over a representative volume, which consists of many fractures, on the basis of micromechanics. This analysis is valid for small displacements before the failure of the joints. In the other micromechanics based model of Yoshida and Horii (), the jointed rock mass is replaced with an equivalent continuum body whose constitutive equation is obtained from the relation between the average stress and strain over a representative volume element. The constitutive equations directly reject the orientation and spacing of the dominant joints and can incorporate strong induced anisotropy. In the crack tensor model of Oda et al. (99) the compliance tensor of the jointed rock mass is given as the sum of the elastic compliance tensor of the base rock and that corresponding to the crack deformation. The crack tensor is defined in terms of the size of the joints and components of the unit vector normal to the joints. The damage tensor model developed by Kawamoto et al. (988) is based on the concept that the effective cross section of the material is reduced due to the damage. When the number of joints in the rock mass is many and it is not possible to obtain information about all of them, thus, it is not possible to deal with each joint individually. It is necessary to replace the JCRM All rights reserved.
2 4 T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp.9-5 jointed rock mass with an equivalent continuum body for analysis with an appropriate associated constitutive model. However, the constitutive models which are discussed above and many others, are very complex and need much input data from experimental or field testing in order to carry out the analysis. So, there is a need for a simple technique where the equivalent continuum method can capture sufficiently well the behavior of a jointed rock mass using minimal input from the field or from tests and experiments. Considering the inherently inhomogeneous nature of rock masses, the practical equivalent continuum approach attempts to use statistical relations, which are simple and obtained after analyzing a large set of data from the literature on laboratory test results relating to jointed rock masses. The material properties are represented by a set of relations which express the properties of the jointed rock mass as a function of a joint factor and the properties of the intact rock. The joint factor for a given jointed rock is estimated based on the joint fabric. The tangent elastic modulus of the intact rock is represented by a confining stress-dependent hyperbolic relation. The jointed rock has been represented as an equivalent continuum using the statistical relations arrived at from the analysis of the experimental data of Roy (99), Arora (987), Yaji (984), Brown and Trollope (97), VidyaBhushan Maji (7) and Einstein and Hirschfield (97). These relations express the tangent modulus of the jointed medium as a function of the joint factor and the tangent modulus of the intact rock. The results have been presented in the form of stress strain curves for the jointed rocks and compared with the experimental results. These results are also compared with explicit modeling of rock mass where the joints are modeled using interface element. The developed model called as Practical Equivalent Continuum Model has been applied to the analyses of well-documented field case studies of large rock excavations and slopes.. EMPIRICAL RELATIONSHIPS An effort has been made to arrive at empirical relations, which express the strength and deformation of jointed rock as a function of intact rock properties and a joint factor. These relations are determined by statistical analysis of a large amount of experimental data from Brown and Trollope (97), Einstein and Hirschfeld (97), Yaji (984), Arora (987), VidyaBhushan Maji (7), and Roy (984). The experimental data covers a wide range of rocklike material and rocks namely plaster of Paris, different kinds of sandstone, granite and gypsum plaster with paralled and unparalled joints with different joint fabric for different confining pressures. Thus, it is hoped that the developed equations for jointed rock mass may have a reasonably wide applicability. Based on the statistical analysis, the uniaxial compressive strength and elastic modulus obtained from uniaxial compressive tests and triaxial tests of jointed rock at different confining pressures are expressed as the function of the joint factor and intact rock properties. So, knowing the intact rock properties and the joint factor, the jointed rock properties can be estimated. Statistical analyses were carried out to arrive at possible empirical relations for the tangent modulus at different confinements and uniaxial compressive strength. A large amount of experimental data of uniaxial compressive strength ratio and elastic modulus ratio versus joint factor of the jointed rock specimens was digitally filtered to reduce the scatter in the data. Linear and nonlinear relationships between the uniaxial compressive strength, tangent elastic modulus at different confinements, and joint factor have been arrived at by using least-squares fitting for linear relationships and Lorentzian minimization for nonlinear relationships. Least-squares minimization assumes that the x values are accurately determined and that an error exists only in the dependant variable y. The errors are assumed to map a Gaussian profile and are normally distributed. Lorenzian minimization is very robust when the data is noisy and also converges quite rapidly. The correlation coefficient and the standard error of the fitted relationship is defined below. The correlation coefficient, r of a relationship fitted to x, y data is expressed as n i r [ y ( y ) ] n i i ( y y) i p i where y i = y value for a given x; y data and the (y p ) i = value of y computed using the relationship fitted; and y = mean of y values. The standard error of the relation fitted is computed as S () n i [ yi ( yp) i] () n In the present analysis the variable y is either uniaxial compressive strength ratio (σ cr ) or modulus ratio (E r ) and variable x is joint factor J f. Both linear and nonlinear relations are obtained for σ cr and E r as a function of a joint factor J f using the procedure described above. The experimental data is compiled in the form of a joint factor J f, the uniaxial compressive strength ratio σ cr, and elastic modulus ratio E r. The experimental data used in the analysis covers a wide range of rocks like plaster of Paris, different kinds of sandstone, granite, and gypsum plaster with filled and unfilled joints with different joint fabric for different confining pressures. The experimental data is digitally filtered to reduce the scatter in the data for a better fit. The values of correlation coefficient for the fitted equations are given for the digitally filtered data set and the standard error of the fit is calculated (Sridevi, ) with respect to the actual experimental data. Rank of the relationship fitted is in the order of minimum standard error. The relation for which the standard error is minimum is ranked as and the next relation is ranked,, and 4 based on the increase in the standard error. Based on the statistical analysis of the data, empirical relationships for the uniaxial compressive strength ratio as a function of joint factor (J f ) are derived. Figure shows the plot of the uniaxial compressive strength ratio versus the joint factor of the experimental data with the empirical relationship fitted for the equation with rank. The elastic modulus ratio from unconfined compressive strength test, from triaxial test, with different confining pressures is considered in the analysis. Linear and nonlinear relationships between the elastic modulus ratio at the three different confining pressures and the joint factor are derived using the experimental data. Much more unconfined compressive strength test data is available compared to triaxial test data. Hence the relationships arrived for elastic modulus using unconfined compressive test data are more reliable than those obtained using triaxial test data. Sample plot of the experimental data with the fitted relationship is shown in Figure. All the 4
3 Ratio of Elastic modulus at = Uniaxial compressive strength ratio T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp statistical relations arrived in the analysis are for < J f < 8. For intact rock, i.e., for J f =, σ cr and E r =. (σ cr ) Plaster of Paris Jamrani sandstone Arora (987) Agra sandstone Plaster of Paris Kota sandstone Yaji (984) Granite Plaster of Paris Singh and Dev (988) Kota sandstone Sharma (989) Einstein and Hirschfeld (97) Brown and Trollope (97) Brown (97) Present experimental data Joint Factor Figure. Uniaxial compressive test data and the fitted relation between J f and σ cr. Data from Arora (987), Yaji (984), Einstein and Hirschfeld(97), Brown and Trollope (97), Brown (97), Roy (99). (E r ) σ Plaster of Paris Agra sandstone Jamrani sandstone Plaster of Paris Kota sandstone Granite Einstein and Hirschfeld (97) Brown and Trollope (97) Brown (97) Joint Factor J f cr e.64 Arora (987) Yaji (984) E r e.86 J f Figure. Experimental data and the fitted relation between J f and E r. Data from Arora (987), Yaji (984), Einstein and Hirschfeld (97), Brown and Trollope (97), Brown (97), and Roy (99). Extensive laboratory testing of intact and jointed specimens of different grades of plaster of Paris, sandstone and granite in uniaxial and triaxial compression revealed that the important factors which influence the strength and modulus values of the jointed rock are (i) joint frequency, J n, (ii) joint orientation, β, with respect to the major principal stress direction and (iii) joint strength. The joint factor (Ramamurthy, 99) for a given jointed rock is estimated using the following equation: J J n f () nr. where, J n is number of joints per meter depth, n is the inclination parameter depending on the orientation of the joint β, r is the roughness or joint strength parameter depending on the joint condition. The value of n is obtained by taking the ratio of log (strength reduction) at β = 9 o to log (strength reduction) at the desired value of β. The joint strength parameter r is obtained from a shear test along the joint and is given as r = τ j /σ nj where τ j is the shear strength along the joint and σ nj is the normal stress on the joint. This inclination parameter is independent of joint frequency. The values of n are given for various orientation angles and the joint strength parameter r is given for various uniaxial compressive strengths of intact rock and they are presented by Ramamurthy (99), based on extensive laboratory testing of rocks. So, knowing the intact rock properties and the joint factor, the jointed rock properties can be estimated using the following groups of empirical/ statistical relations. Group : The jointed rock is modeled using the statistical relations given as below. cj J.4.89exp f cr (4) 6. ci E j J E a b f r exp (5) E c i These relationships arrived are based on the statistical analysis of a large amount of experimental data. E r is the tangent modulus ratio, E j is the tangent modulus of jointed rock, E i is the tangent modulus of intact rock, J f is the joint factor, a, b, and c are statistical constants. The values of these constants are given in Table for,. and 5. MPa confining pressures. For the confining pressures other than listed in the Table the values of a, b, and c have been linearly interpolated or extrapolated based on the value of a, b, and c for 5. MPa confining pressure. Table. Values of empirical constant a, b, c for different confining pressure. Confining Pressure, MPa Value of a in Eq. (5) Value of b in Eq. (5) Value of c in Eq. (5) Group : The jointed rock is modeled using the following empirical relations given by Ramamurthy (99). cj c r = =exp-.8 Jf (6) ci Ej Er = J E σ - =exp-.5 f σ i Tangent elastic modulus of jointed rock for any other σ is derived from the tangent elastic modulus of jointed rock at = using the formula given below. Ej Ej ( ) cj -exp -. Equation (7) is valid for =. For different confining pressures, the elastic modulus of jointed rock is calculated using Equation (8), where j is obtained from Equation (6) and E j at = is obtained from Equation (7). These relations were given by Ramamurthy and Arora (994) based on laboratory studies on numerous artificial joints. (7) (8) 4
4 4 T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp.9-5 Group : The jointed rock is modeled by a set of statistical relations after Sitharam (7) as given below. cj c r = =exp-.65j (9) f ci i Ej Er = =exp a J E f () These relations are arrived at Equations (9) and () by fitting the equations of the same form as Ramamurthy (99) to the experimental data using statistical analysis. The value of the empirical constant a is given in tabular form Table for different confining pressures. For confining pressure other than those listed in the Table the value of a can be interpolated or extrapolated based on the value of a for to 7. MPa confining pressure. Table. Value of empirical constant a for different confining pressure. Confining Pressure, MPa Value of a HYPERBOLIC STRESS-STRAIN RELATIONSHIPS The constitutive relations for the intact and jointed rock used in the equivalent continuum analysis are described in detail in this section. The stress strain behavior of rocks over a wide range of stress field is nonlinear and dependent upon confining pressure. The nonlinear elastic confining stressdependent model following a hyperbolic relation proposed by Duncan and Chang (97) is used in the present study. The material behavior of intact rock is modeled using the following non-linear relation: d( ) E ( ) E R E ti r f ti d( a) ( ) f ar f Eti f () n E K P P () ti a a ( ) f R f () ( ) ult Where, E t is the tangent modulus computed at 5 percent of failure stress for the intact rock, E ti is the initial tangent modulus, is the major principal stress, is the confining pressure, a is axial strain, K is the modulus number, n is modulus exponent and P a is atmospheric pressure. ( ) f is the failure stress and ( ) ult is the asymptotic value of stress. The value of modulus exponent K and modulus number n are determined from the plots of tangent modulus versus confining pressure for the intact rock. For brittle failure in rocks the failure stress and the asymptotic value of stress are almost same so R f =. For heavily jointed rocks which have a similar stress-strain relationship of soils, the asymptotic value of stress and failure stress are slightly different and so R f is taken as.9. For the linear stress-strain relation an R f =. is taken. 4. EQUIVALENT CONTINUUM MODEL First implementation of this equivalent continuum model has been done with finite element (Sitharam, ), and the numerical model has been developed from an existing finite element code for a non-linear soil structure interaction program (NLSSRIP). As mentioned, the jointed rock is represented as an equivalent continuum whose properties have been derived from the intact rock properties and a joint factor based on the joint fabric. Nonlinearity in the finite element analysis has been incorporated in the form of material non-linearity of both the intact and jointed rock. The incremental method is used for the solution of the non-linear problem by the finite element method. The load is increased in a series of steps or increments. Each increment is analyzed twice: the first time using the moduli values for the elements based on the stresses at the beginning of the increment; and the second time using moduli value based on the average stresses during the increment. The changes in stress and strain of the elements and changes in the nodal displacements during each increment are added to the values at the beginning of the increment. At the beginning of each new increment of loading, an appropriate modulus value is selected for each element on the basis of the values of stress or strain in that element. The non-linear stress strain behavior is approximated by a series of straight lines. The displacements at each increment of loading are accumulated to give the total displacement at any stage of loading, and the incremental process is repeated until the total load is applied. The principal advantages of this procedure are its complete generality and its ability to provide a relatively complete description of the load deformation behavior. Initial stresses may be readily accounted for as the tangent modulus is expressed in terms of the stresses only. In the equivalent continuum approach, the discontinuous rock body is modeled using four-node quadrilateral elements, the properties of each element being defined in terms of some combination of the properties of the intact rock and those of the joints. Non-linearity in the finite element analysis has been incorporated in the form of material non-linearity of both the intact and jointed rock. During loading, if a rock element is found to fail in shear, this is noted, but no changes are effected, and the element is allowed to follow the hyperbolic relation as before, in keeping with the non-linear elastic formulation of the problem. During analysis, it has been found that some elements which fail at one stage recoup at a subsequent higher loading stage. If the elements fail in tension, they are assigned very small values of the elastic modulus for the subsequent loads. In the case of explicit representation of joints, the intact rock is represented by four-node quadrilateral elements and the joints are represented by a two dimensional gap and friction elements. This element is a two-node non-linear interface element used to model node-to-node contact between two bodies with or without friction. It is represented by a pair of coupled non-linear orthogonal springs in the normal and tangential directions to the interface, which are assumed to be very stiff, relative (with stiffness, K n and K t ) to the bodies they are attached to. The Coulomb law is used for friction. Frictionless contact may be modeled by specifying a zero coefficient of friction. The element may assume open or 4
5 T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp a closed status on relative displacement in the normal direction. The closed status may be sticking or sliding depending on whether the friction limit μ f n is reached, where μ is the coefficient of friction and f n is the normal compressive force in the gap. Later the model has also been incorporated in the commercial finite difference code First Langrangian Analysis of Continua FLAC (Sitharam and Madhavi Latha, ). A FISH function has been written to incorporate joint factor model with Duncan-Chang non-linear hyperbolic relationships in FLAC. The work verifies the validity of proposed model for different field case studies, namely two large power station caverns, one in Japan and the other in Himalayas and Kiirunavara mine in Sweden. Sequential excavation was simulated in the analysis by assigning null model available in FLAC to the excavated rock mass in each stage. The settlement and failure observations reported from field studies for these different cases were compared with the predicted observations from the numerical analysis in this study. The results of numerical modeling applied to these different cases are systematically analyzed to investigate the efficiency of the numerical model in estimating the deformations and stress distribution around the excavations. Results indicated that the model is capable of predicting the settlements and failure observations made in field fairly well. Results from this study confirmed the effectiveness of the practical equivalent continuum approach and the joint factor model used together for solving various problems involving excavations in jointed rocks. Recently, the applicability of the model has been extended by incorporating the model in the FLAC-D (Sitharam, 5). Initially, it has been validated with simple element tests i.e by simulating -D triaxial testing of rock specimen and is compared with experimental results as well as explicit modeling. Later, the applicability of the model for field problems is investigated by undertaking numerical modeling of the Nathpa Jhakri powerhouse cavern, India, considering -dimensional geometry and stresses along with stages of excavation. 4. Validation of the model 4.. Equivalent continuum model For validation purposes, the finite element analysis has been carried out for single-, multiple- and block jointed rocks (Figure ) of sandstone, granite, Agra sandstone and gypsum plaster using the proposed model. The jointed rock properties are expressed using Equations (6)-(8). The properties of the intact rocks used for the simulation are presented in Table. The results have been plotted in the form of stress strain curves and compared with the experimental results. The stress strain curves compare well with the experimental results for single-, multiple- and block jointed specimens with different confining pressures and joint inclination angles. Some sample stress strain plots for multiple-jointed specimens of Agra Sandstones is shown in Figure 4 and block-jointed specimens of gypsum plaster for different confining pressures shown in Figure 5 along with the experimental results. It can be seen from Figures 4 and 5 that the equivalent continuum model provides a good approximation of the jointed rock behavior. 4.. Comparison of the equivalent continuum model with explicit modeling of joints The results obtained using explicit modeling of joints for single-jointed specimens of sandstone and granite and multiple-jointed specimens of Agra sandstone are compared with the results obtained for the same specimens using equivalent continuum analysis. The comparison of stress strain curves for single-jointed specimens of sandstone is given in Figure 7. The experimental results of Yaji (984) and Arora (987) are also plotted in the figure for comparison. Though only a few comparison plots are shown here, the trend of the results is consistent for different rocks, different joint inclination angle and confining pressure. Figure. Jointed rock specimens and corresponding finite element. Table. Properties of intact rocks used for the numerical modeling. Property Jamrani sandstone Agra Sandstone Gypsum Plaster Mass density, kn/m UCS, MN/m 7 Mod. No. K Mod.Exp. n Cohesion, MN/m Friction, o Classification Hard Rock Very hard rock Soft Rock 4. Validation of the FLAC model 4.. Equivalent continuum model Element tests of jointed rock mass with one to four joints with different orientation and subjected to different confining 4
6 Deviatoric stress (MPa) Deviatoric Stress(MPa) Deviatoric stress (MPa) 44 T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp.9-5 pressure have been simulated using the proposed equivalent continuum model. To validate the results of equivalent continuum model, it is compared with the actual laboratory Experimental, joint Experimental, joints Experiemntal, joints Corresponding FLACD-ECM results Corresponding FLACD-explicit results One joint Two joints Three joints 4.. Explicit modeling using FLAC Explicit modeling of jointed rock sample has also been done to know the efficiency of Numerical modeling with FLAC-D. Introduction of interfaces in cylindrical sample FLAC-D is complicated job and creating more than 4-5 interfaces is not recommended. Here, up to four joints have been tried. Figure 6 shows the cylindrical sample with one and three joints as interfaces. The material properties for explicit modeling, assumed to be elastic which is given in Table Axial strain (%) Figure 4. Stress strain plot for different confining pressures (experimental data after Arora, 987). m Block separation to Interfaces introduce interface Figure 6. Cylindrical samples with interfaces σ =.8 MPa σ =6.9 MPa σ =.4 MPa σ =.4 MPa Table 4. Properties of rocks for explicit modeling. Properties Sandstone Agra Sandstone Density, kn/m.5.7 UCS, MPa 7 Cohesion, MPa. 9. Friction, degrees Elas. Mod., GPa 5. Classification Hard rock Very Hard rock Axial Strain (%) Figure 5. Stress strain plot for different confining pressures (experimental data after Brown and Trollope, 97) Experimental, σ =5MPa Experimental, σ =.5MPa Experimental, σ =MPa FLACD Equivalent continuum FLACD Explcit σ =5MPa σ =.5MPa experiment results and also by explicit modeling of joints in rock. Equivalent continuum analysis has been carried out for two types of rocks namely Jamrani sandstone and Agra sandstone for both intact and jointed specimens (Arora, 987; Yaji, 984). The properties of the intact rocks used for the simulation are presented in Table 6. The jointed rock properties are expressed using Equations (), () and (). The stress-strain curves for each load increment has been studied for different values of confining pressure for both intact and jointed specimens with different orientations of joints. The experimental values (Arora, 987) then compared with the numerical results. However, testing was performed in three different confining pressures namely,.5 and 5 MPa. Various Plot for deviatoric stress vs axial strain, at respective confining pressure, for both experimental as well as numerical testing results are shown in Figure 7. It can be seen from the results that the equivalent continuum model developed for the jointed rockmass matches well with the actual experimental test results Axial strain (%) σ =MPa Figure 7. Comparison of numerical and experimental results of sandstone with single joint at 6 degree inclination for,.5 and 5 MPa confining pressure. 5. FIELD CASE STUDIES 5. Finite Element analysis of Shiobara power house cavern The equivalent continuum model developed has been applied for the analysis of a large-scale cavern in jointed rock mass for the Shiobara power station in Japan (Horii et al., 999). This case study was selected for modeling using the equivalent continuum model because the complete field 44
7 Displacement (cm) Displacement (cm) Displacement (cm) T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp BI BI BI5 BI7 BI9 m 5m 8m m 4.7MPa.5MPa 5 BI BI BI4 BI6 BI8 Figure 8 Figure 9 Figure 8. Cross section of cavern and location of displacement transducer Figure 9. Finite element mesh for the cavern and the surrounding rock. the value of the joint inclination parameter n =.5 and the joint strength parameter r is based on the uniaxial compressive strength of the rock and is taken as r =.9. The cavern, along with the surrounding rock, has been analyzed using the above approach at the completion of the whole excavation. The results of the finite element analysis have been plotted in terms of relative displacement at the completion of the whole excavation along the each measurement line. The measured values of relative displacement given by Horii et al. (999 and ) are also plotted for comparison. The relative displacement versus distance from the cavern wall along the measurement lines BI, BI, BI7 and BI8 are as shown in Figures () to () along with the measured values. information, including displacement data, was available for possible analysis and comparison with the proposed model. The Shiobara power station cavern constructed by the Tokyo Electric Company is a large cavern for a pumped storage power station with a maximum output of 9 MW. The rock mass surrounding the cavern was characterized mainly as rhyolite consisting of platy and columnar joints. The cavern was located at a depth of m below the ground level. The cavern (Figure 8) measured 8 m in width, 5 m in height and 6 m in length. The amount of jointed rock mass excavated due to opening of the cavern was estimated to about 9, m. The three in situ principal stresses were recorded as 5.,.9 and.8 MPa. The reported average intact rock compressive strength and elastic modulus were 8. MPa and 4. GPa, respectively. The elastic modulus of the jointed rock mass was in the range of 5 GPa and the strength parameters measured for the jointed rock mass were c = MPa and = 45 o. Multi-point bore hole extensometer (MPBX) data were available at several locations along different measurement lines BI to BI9 (Figure 8) around the cavern for a possible comparison with our model. The cross section of the cavern along with the location of the MPBXs is as shown in Figure 8. The measurement lines along which the field displacements are available are also shown in Figure 8. The jointed rock mass surrounding the cavern has been analyzed by the finite element method using the proposed equivalent continuum approach. Equivalent material properties for jointed rock were modeled using the statistical relations given in Equations (4) and (5). The values for a, b, and c constants used here are the same as the ones listed in Table. Figure 9 shows the finite element mesh ( m m) used for the above problem (Sitharam, ). The problem was analyzed with the initial stresses existing in the surrounding rock. The joint properties, the dip angle and the average spacing of the three dominant joints present in the surrounding rock are given in Table 5. The joint factor for the above problem has been estimated as follows. The number of joints per meter depth in the surrounding rock which have an influence on the cavern determine the joint frequency value J n as 5. Since the most critical joint has a joint orientation of o, Table 5. Properties of joint sets for Shiobara power house cavern. Joint set Dip angle spacing Joint frequency (J n) Joint inclination parameter (n) Joint strength parameter (r) Joint factor (J f) J n/(nr) I 6 R cm II 6 L cm III L cm Observed Values (after Horii et.al) Eq Continuum resullts 5 5 Distance from Cavern wall (m) Figure. Measured (after Horii et al., 999 and ) and calculated relative displacements along the measurement line BI. 5 4 Observed Values (after Horii et.al) Eq Continuum resullts 5 5 Distance from Cavern wall (m) Figure. Measured (after Horii et al., 999 and ) and calculated relative displacements along the measurement line BI 5 4 Observed Values (after Horii et.al) Eq Continuum resullts 5 5 Distance from Cavern wall (m) Figure. Measured (after Horii et al., 999 and ) and calculated relative displacements along the measurement line BI7. 45
8 Displacement (cm) 46 T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp Observed Values (after Horii et.al) Eq Continuum resullts minor horizontal virgin stress (Perpendicular to the ore body) in MPa and z m is the mine level in meters (z-axis positive upwards). 5 5 Distance from Cavern wall (m) Figure. Measured (after Horii et al., 999 and ) and calculated relative displacements along the measurement line BI8. 5. FLAC analysis of Kiirunavaara mine,sweden. The Kiirunavaara mine, which is 4 m long, with an average width of 9 m, is located 44 kilometers north of the Arctic Circle in the city of Kiruna in North Sweden and currently produces Mt of magnetic iron ore annually. The mine was initially mined as open pit, starting at the top of Kiirunavaara mountain. Underground mining using sublevel caving is started in 95s. The mine strikes nearly northsouth and dips, on average 6 to the east. The ore body is relatively strong surrounded by competent quartz porphyry on the hanging wall and syenite porphyry on the footwall. The rock mass is relatively jointed, with two to three joint sets occurring in 8 different structural domains of the mine. One joint set is oriented roughly parallel to the ore body as the other two strike obliquely to the ore body. All dip fairly steeply (6 to 9). A representative section of the mine, namely section Y is selected for numerical modeling. The geometry of the mine at this section showing different mining levels is given in Figure 4. The mine has experienced large-scale stability problems. Because of sub level caving, the hanging wall continuously caves in as the ore is mined. Due to this, large subsidence of the ground surface is observed. Instabilities and large-scale failures were observed both in hanging wall and footwall of the mine. Signs of failure were first observed in the footwall in the year 985, with more widespread cracking underground in year 989. The locations where first cracks were observed in 985 are mapped by Sjoberg (999) and are presented in Figure 5. Failure surface can be drawn connecting these failure observations, which can be compared using numerical simulations from present study. Qualitative and quantitative estimation of this instability is very important as the city of the Kiruna and the railroads are located on the hanging wall side of the mine. Several researchers attempted to predict the instabilities using different approaches. Sjoberg (999) simulated the failure using numerical model incorporating Hoek and Brown failure criterion, which involves exact evaluation of various parameters for the equivalent continuum model. The input parameters are selected for the numerical analysis from the observations reported by Lupo (996) as given in Table 6. The input for the virgin stress distribution is taken from the average stress components presented by Sjoberg (999) as: Table 6. Input parameters for numerical modeling of Kiirunavaara mine (after Lupo, 996). Parameter Hanging Foot wall wall Orebody Uniaxial Compressive strength, MPa 4 5 Elastic Modulus, MPa Poisson's ratio..7.5 Density, kg/m Cohesion, MPa..6.6 Angle of internal friction Results from the numerical analysis showed that the type of failure occurred in this model are primarily shear failure. In continuum type of model, since the displacements are continuous, we can not see the development of cracks as observed in the field. Failure can only be observed from the concentration of shear strains in the model. The path of concentrated shear strains represents the failure surface in the model. The failure thus simulated by the numerical model using practical equivalent continuum approach can be compared with the failure observations in the field. Shear failure was observed in the footwall of the model while excavating for the mining level of -586 m, agreeing with the field observations as reported by Sjoberg (999). Typical failure surface for the mining step of 586 m is presented in Figure 5. It can be observed from the figure that failure develops along a curved surface within the foot wall. The figure is zoomed between the excavation steps -65 m to -586 m so as to clearly visualize the failure surface. The shear strains are concentrated near the slope, at a depth of 55 m, which is comparable with the field observations presented in Figure 5. The predicted failure surface matched fairly well with the field observations. Ground Surface -6 m Open pit -9 m - m Hanging Wall -65 m m m 984 Foot Wall -586 m m m 99-7 m Ore body Figure 4. FLAC model for Kiirunavaara mine showing different mining levels. v = -.7z m -.6, H = -.6z m -.6, h = -.4z m -.4. Where v is the vertical virgin stress, H is the major horizontal virgin stress (Parallel to the ore body), h is the 46
9 T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp overburden is taken into consideration by applying equivalent amount of pressure at the top. Analysis procedure was similar to FLAC-D as mentioned above accept that here the - dimensinal geometry has been considered and the displacements values measured were consistent with the actual experimented values. Figure 5. Shear strains for FLAC model for 586 m mining level. 5. FLAC Analysis: Naptha Jhakri Power House Cavern Nathpa-Jhakri hydropower project in the state of Himachal Pradesh in India involved a major opening for powerhouse cavern of dimensions 6 m m 49 m (length width height) at a depth of 6.5 m below the ground level (Figure 6). The opening is located in the left bank at about 5 m from the Sutlej river. The rock mass around the cavern is mainly quartz mica schist with joints and discontinuities (Bhasin 996). The measured in-situ stresses are 4.7 and 6.4 MPa in E-W and N-S and directions respectively and 5.89 MPa in vertical direction. The properties of intact rock and joints are given in Tables 7 and 8 respectively (Varadarajan, ). The cavern is excavated in stages. The sequence of excavation with the elevation of each excavation step and the locations of installation of multi point borehole extensometers (MPBX) for displacement measurements are shown in Figure 7. Since the cavern is symmetric, only half of the cavern is analyzed (Sitharam and Madhavi Latha, ). The finite difference grid used for the analysis is of size m 45 m with rectangular zones. The excavation steps are simulated in the numerical analysis and the locations of the installation of extensometers are identified for obtaining the displacements for comparison with the measured displacements from instrumentation of the cavern in field. The variation of displacements with time is also obtained from numerical analysis by solving for equilibrium after each excavation step. Comparison of the observed and predicted deformations along the measurement line at different locations for various excavation levels after the completion of excavation is presented in Table 9. In all the cases, the minimum value of deformation corresponds to the measuring point of the MPBX positioned close to the cavern wall and the maximum value corresponds to the farthest measuring point from the wall. It can be observed from Table 9 that the numerical model is efficient in predicting the deformations around the cavern wall. The variation of displacement with time for a particular excavation step and with progress of excavation steps is well represented in the analysis. Analysis of the cavern has also been extended to FLAC- D for better understanding of stress and deformational behavior of the cavern. For the -D numerical simulation, the grid selected is shown in Figure 7. As shown in Figure 7, only 4 m overburden is taken into model, rest m Figure 6. Excavation sequence and locations of extensometers for Nathpa-Jhakri powerhouse cavern. Table 7. Properties of intact rock for Nathpa-Jhakri powerhouse cavern. Modulus number K 45 Modulus exponent n.4 Cohesion, MPa 6.8 Poisson ratio.9 Angle of internal friction, degrees Table 8. Properties of joints for Nathpa-Jhakri powerhouse cavern. Frequency of joints J n 5 Orientation of critical joint 55 Inclination parameter n.8 Joint strength parameter r.6 Joint factor J f.9 47
10 89m 45m 48 T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp.9-5 Table 9 Measured and predicted deformations at for Nathpa- Jhakri cavern. Excavation Deformation along the Location of level line (mm) From Stage To MPBX EL, at EL, Observed Predicted EL, m m m (A) (B) (B) -. to (C) (D) (E) Figure 7. D Grid of the power house cavern with FLAC-D (All dimensions are in meters). 5.4 Seismic Stability Analysis of A River Abutment Slope The Chenab river forms about 6 m deep gorge in a V- shaped valley in the area between Bakkal and Kauri villages (Figure 8). The railway link between Katra and Laole section in Jammu and Kashmir has been planned to have a steel arch bridge on Chenab river with around 95 m span at a height of 6 m. The rocks present at the bridge site are heavily jointed. The stability analysis of the right abutment slope at Kauri side of Chenab river between Katra and Laole, in Jammu and Kashmir, India was simulated using FLAC as a plane strain problem. Seismic stability analyses have been done using pseudostatic approach and also by applying complete time history of a real earthquake. The assessment of slope stability of large slopes in jointed rocks requires the use of non-linear material models. Here, the concept of joint factor (J f ) along with hyperbolic model has been used for the study. The subsurface essentially consists of dolomitic limestone with different degrees of weathering and fracturing. The rock mass ratings for surface weathered and fractured dolomitic limestone are given in Table. The intact rock properties of the right abutment of the slope are presented in Table. Calculated equivalent rock mass properties, cohesion.785 MPa, friction angle o, Young s modulus 4.4 GPa, and Poisson s ratio.5. The Hoek and Brown Parameter m for broken rock.59 where as the s value calculated for RMR of 4 is.7. Table. The ratings for surface weathered and fractured dolomitic limestone (NIRM Report 4). Parameter Fractured Surface weathered dolomitic dolomitic limestone limestone RQD 49 Q 6..5 RMR 48 5 GSI 4 7 J v 6 Table. Intact rock properties for the slope at right abutment (NIRM Report 4). Property Average value Density, kg/m 76 Young s modulus, GPa 65 Poisson s ratio.5 UCS, MPa 5 Cohesion, MPa Friction, degrees.76 Hoek & Brown parameter, m.5 The slope is simulated here, with FLAC in plane strain using small strain mode. Relatively high discretization with approximately 6 zones for a slope height of 59 m is considered in the analyses. Horizontal displacements are fixed for nodes along the left and right boundaries while both horizontal and vertical displacements are fixed along the bottom boundary. 6m Kauri side Right abutment Bakkal side Left abutment 8m 69m Figure 9 Numerical mesh of the slope using FLAC. Figure 8. Slope at the Chenab river site. Slope stability analysis was done with the shear strength reduction technique, where the simulations are run for a series of trial factor of safety F trial with c and adjusted according to the equations: 48
11 Acceleration (g) T. G. SITHARAM / International Journal of the JCRM vol.5 (9) pp c trial trialc (4) F.5.5 trial tan trial tan F (5).5.5 Corrected transverse component of acceleration The value F trial at which slope fails is found efficiently in FLAC using bracketing and bisection method. The static factor of safety for the slope with rock mass properties of cohesion.785 MPa and friction o is found to be.86. Analyses represent the effects of earthquake shaking by pseudostatic accelerations that produce inertial forces, which act through the centroid of the failure mass. The magnitude of the horizontal pseudo static force is a W F h h k W g h (6) Where, a h is horizontal pseudostatic accelerations, k h is dimensionless horizontal pseudostatic coefficients, and W is the weight of the failure mass. The horizontal pseudostatic forces are assumed to act in directions that produce positive driving moments. With the values of cohesion.785 MPa and friction o, the factor of safety for the peak a h =. was found.. Figure shows the maximum shear strain rate for the slope along with velocity vectors. Failure surface is almost circular from top of the slope to the toe. Figure shows the plasticity indicator plot of the slope. They reveal those zones in which the stresses satisfy the yield criterion. A failure indicated by the contiguous line of active plastic zones that join two surfaces. The yielded zone due to shear can be seen in the Figure and as observed, it also follows a circular/non-circular failure surface. Plastic Zone In tension Elastic Zone Figure Figure Figure. Critical failure surface. Figure. Plasticity indicator plot depicting the yielded zones. To carry out the full dynamic analysis with the real earthquake records, the corrected transverse component of acceleration time history was applied (Figure ) at the base of the slope. Free field boundary has been used to minimize the wave reflection. It has been observed that the overall effect of continuous inertial forces may lead to accumulation of the displacement of a particular section of the slope. Once the applied ground motions generated due to inertia forces have ceased, no further deformation has occurred, as there is no marked loss in the strength of rock. The variation in displacements along the slope values were recorded and studied for different Rayleigh damping values Time (s) Figure The corrected transverse component of acceleration time history of Uttarkashi earthquake, Oct,99 :5 IST, used in the study (Data from Dept. of Earthquake Engg., IIT Roorkee). 6. CONCLUDING REMARKS The paper presents the summary of the work on jointed rock mass along with development of equivalent continuum model, validation and field case studies. The statistical equations for representing the jointed rock mass properties are very simple and give a fair estimate of jointed rock properties in the absence of reliable experimental data. The accuracy of the estimation of the jointed rock properties depends upon the estimation of joint factor. Joint factor is very important for the description of jointed rock and to find the correlation between the parameters and the geological data. Since the database covers a wide range of rock properties the statistical relationships arrived more or less give a good estimate of uni-axial compressive strength and elastic modulus for all rock types. Jointed rock has been successfully modeled as an equivalent continuum whose properties represent the properties of the jointed rock. It can be seen from the results that the equivalent continuum model developed works well for single-, multiple- and block-jointed rock with different joint fabric and joint orientation under a wide range of confining pressures. In the equivalent continuum model, the jointed rock properties are expressed as a simple function of intact rock properties and joint properties. As the joints are not modeled separately, the analysis becomes simpler as highly complex joint fabrics can be modeled using a simple finite element mesh / finite difference (FLAC) unlike explicit modeling. The major advantage of explicit modeling of joints is that the mode of failure in the jointed rock mass can be reasonably estimated and the zones most susceptible to failure can be identified which is not possible in the equivalent continuum modeling. Explicit modeling of joints is efficient only when the jointed rock mass has few major joints. Equivalent continuum analysis can be applied to a single-jointed rock mass through to heavily jointed rock mass effectively without compromising the accuracy of the results. Moreover, the input data required for the equivalent continuum model are minimal. It can be inferred from the results that the equivalent continuum model developed provides a reasonable estimate of rock mass behavior in the absence of detailed experimental data. The only input data required for the analysis are the properties of intact rock and the joint properties for estimating the joint factor. It can be concluded from the analysis of various case studies that the 49
MEMORANDUM SUBJECT: CERTIFICATE IN ROCK MECHANICS PAPER 1 : THEORY SUBJECT CODE: COMRMC MODERATOR: H YILMAZ EXAMINATION DATE: OCTOBER 2017 TIME:
MEMORANDUM SUBJECT: CERTIFICATE IN ROCK MECHANICS PAPER 1 : THEORY EXAMINER: WM BESTER SUBJECT CODE: COMRMC EXAMINATION DATE: OCTOBER 2017 TIME: MODERATOR: H YILMAZ TOTAL MARKS: [100] PASS MARK: (60%)
More informationModule 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression
FAILURE CRITERIA OF ROCK AND ROCK MASSES Contents 5.1 Failure in rocks 5.1.1 Hydrostatic compression 5.1.2 Deviatoric compression 5.1.3 Effect of confining pressure 5.2 Failure modes in rocks 5.3 Complete
More informationBehaviour of Blast-Induced Damaged Zone Around Underground Excavations in Hard Rock Mass Problem statement Objectives
Behaviour of Blast-Induced Damaged Zone Around Underground Excavations in Hard Rock Mass Problem statement Blast-induced damaged zone can affect the affect stability and performance of tunnel. But, we
More informationON THE FACE STABILITY OF TUNNELS IN WEAK ROCKS
33 rd 33 Annual rd Annual General General Conference conference of the Canadian of the Canadian Society for Society Civil Engineering for Civil Engineering 33 e Congrès général annuel de la Société canadienne
More informationSeismic analysis of horseshoe tunnels under dynamic loads due to earthquakes
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2010 Seismic analysis of horseshoe tunnels under dynamic loads due to earthquakes Navid
More informationAN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS
AN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS S. Kontoe, L. Pelecanos & D.M. Potts ABSTRACT: Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified
More informationIntroduction and Background
Introduction and Background Itasca Consulting Group, Inc. (Itasca) has been participating in the geomechanical design of the underground 118-Zone at the Capstone Minto Mine (Minto) in the Yukon, in northwestern
More informationThe effect of stope inclination and wall rock roughness on backfill free face stability
The effect of stope inclination and wall rock roughness on backfill free face stability Dirige, A. P. E., McNearny, R. L., and Thompson, D. S. Montana Tech of the University of Montana, Butte, Montana,
More informationEffect of buttress on reduction of rock slope sliding along geological boundary
Paper No. 20 ISMS 2016 Effect of buttress on reduction of rock slope sliding along geological boundary Ryota MORIYA *, Daisuke FUKUDA, Jun-ichi KODAMA, Yoshiaki FUJII Faculty of Engineering, Hokkaido University,
More informationROCK MASS PROPERTIES FOR TUNNELLING
ROCK MASS PROPERTIES FOR TUNNELLING Robert Bertuzzi 2 nd November 2017 1 Driver Estimating the strength and deformation characteristics of a rock mass for tunnel design is generally based on empiricism
More informationFurther Research into Methods of Analysing the October 2000 Stability of Deep Open Pit Mines EXECUTIVE SUMMARY
EXECUTIVE SUMMARY This report presents the results of a program of further research into the use of a combined approach of numerical and centrifuge modeling in assessing the stability of deep open pit
More informationFailure and Failure Theories for Anisotropic Rocks
17th international Mining Congress and Exhibition of Turkey- IMCET 2001, 2001, ISBN 975-395-417-4 Failure and Failure Theories for Anisotropic Rocks E. Yaşar Department of Mining Engineering, Çukurova
More informationPLANES OF WEAKNESS IN ROCKS, ROCK FRCTURES AND FRACTURED ROCK. Contents
PLANES OF WEAKNESS IN ROCKS, ROCK FRCTURES AND FRACTURED ROCK Contents 7.1 Introduction 7.2 Studies On Jointed Rock Mass 7.2.1 Joint Intensity 7.2.2 Orientation Of Joints 7.2.3 Joint Roughness/Joint Strength
More informationGROUND RESPONSE AND SUPPORT MEASURES FOR PIR PANJAL TUNNEL IN THE HIMALAYAS
Ground IGC 2009, Response Guntur, and INDIA Support Measures for Pir Panjal Tunnel in the Himalayas GROUND RESPONSE AND SUPPORT MEASURES FOR PIR PANJAL TUNNEL IN THE HIMALAYAS K.S. Rao Professor, Department
More informationPillar strength estimates for foliated and inclined pillars in schistose material
Pillar strength estimates for foliated and inclined pillars in schistose material L.J. Lorig Itasca Consulting Group, Inc., Minneapolis, MN, USA A. Cabrera Itasca S.A., Santiago, Chile ABSTRACT: Pillar
More informationSession 3: Geology and Rock Mechanics Fundamentals
Session 3: Geology and Rock Mechanics Fundamentals Geotechnical Engineering Appreciation Course (Jointly organised by IES Academy and GeoSS) Dr Zhou Yingxin, Senior Principal Engineer, DSTA Adjuct Associate
More informationApplication of a transversely isotropic brittle rock mass model in roof support design
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2012 Application of a transversely isotropic brittle rock mass model in roof support
More informationNumerical Modelling of Dynamic Earth Force Transmission to Underground Structures
Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures N. Kodama Waseda Institute for Advanced Study, Waseda University, Japan K. Komiya Chiba Institute of Technology, Japan
More informationNumerical modelling for estimation of first weighting distance in longwall coal mining - A case study
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2012 Numerical modelling for estimation of first weighting distance in longwall coal
More informationMining. Slope stability analysis at highway BR-153 using numerical models. Mineração. Abstract. 1. Introduction
Mining Mineração http://dx.doi.org/10.1590/0370-44672015690040 Ricardo Hundelshaussen Rubio Engenheiro Industrial / Doutorando Universidade Federal do Rio Grande do Sul - UFRS Departamento de Engenharia
More informationFlin Flon Mining Belt
EOSC433: Geotechnical Engineering Practice & Design Lecture 7: Stress Analysis around Underground Openings 1 of 40 Erik Eberhardt UBC Geological Engineering EOSC 433 (2007) Flin Flon Mining Belt Since
More informationNumerical investigation of EDZ development around a deep polymetallic ore mine
Paper No. 198 ISMS 2016 Numerical investigation of EDZ development around a deep polymetallic ore mine Mountaka Souley a *, Marwan Al Heib a, Vincent Renaud a a INERIS, c/o Ecole des Mines de Nancy, Campus
More informationHaulage Drift Stability Analysis- A Sensitivity Approach
Haulage Drift Stability Analysis- A Sensitivity Approach W. Abdellah University of Assiut, Assiut, Egypt ABSTRACT Haulage drifts are the primary access to the mining blocks of an ore body in a multi-level
More informationStability Assessment of a Heavily Jointed Rock Slope using Limit Equilibrium and Finite Element Methods
Indian Geotechnical Conference 2017 GeoNEst 14-16 December 2017, IIT Guwahati, India Stability Assessment of a Heavily Jointed Rock Slope using Limit Equilibrium and Finite Element Methods Aswathi CK Amalesh
More informationModule 9 : Foundation on rocks
LECTURE 32 9.3 BEARING CAPCITY contd... 9.3.2 Safe bearing pressure There are different methods are available to determine the safe bearing pressure on rocks. However, the applicability of different methods
More informationAnalysis of Controlling Parameters for Shear behavior of Rock Joints with FLAC3D
International Journal of Structural and Civil Engineering Research Vol. 5, No. 1, February 2016 Analysis of Controlling Parameters for Shear behavior of Rock Joints with FLAC3D Prasoon Tiwari and Hakan
More informationModel tests and FE-modelling of dynamic soil-structure interaction
Shock and Vibration 19 (2012) 1061 1069 1061 DOI 10.3233/SAV-2012-0712 IOS Press Model tests and FE-modelling of dynamic soil-structure interaction N. Kodama a, * and K. Komiya b a Waseda Institute for
More informationSOIL MODELS: SAFETY FACTORS AND SETTLEMENTS
PERIODICA POLYTECHNICA SER. CIV. ENG. VOL. 48, NO. 1 2, PP. 53 63 (2004) SOIL MODELS: SAFETY FACTORS AND SETTLEMENTS Gabriella VARGA and Zoltán CZAP Geotechnical Department Budapest University of Technology
More informationDynamic Analysis Contents - 1
Dynamic Analysis Contents - 1 TABLE OF CONTENTS 1 DYNAMIC ANALYSIS 1.1 Overview... 1-1 1.2 Relation to Equivalent-Linear Methods... 1-2 1.2.1 Characteristics of the Equivalent-Linear Method... 1-2 1.2.2
More informationTHE VOUSSOIR BEAM REACTION CURVE
THE VOUSSOIR BEAM REACTION CURVE Yossef H. Hatzor Ben-Gurion University, Department of Geological and Environmental Sciences Beer-Sheva, Israel, 84105 ABSTRACT: The influence of joint spacing (s) on the
More informationReciprocal of the initial shear stiffness of the interface K si under initial loading; reciprocal of the initial tangent modulus E i of the soil
Appendix F Notation a b B C c C k C N C s C u C wt C θ D r D 1 D 2 D 10 D 30 Reciprocal of the initial shear stiffness of the interface K si under initial loading; reciprocal of the initial tangent modulus
More informationUnderground Excavation Design Classification
Underground Excavation Design Underground Excavation Design Classification Alfred H. Zettler alfred.zettler@gmx.at Rock Quality Designation Measurement and calculation of RQD Rock Quality Designation index
More informationInfluence of foliation on excavation stability at Rampura Agucha underground mine
Recent Advances in Rock Engineering (RARE 2016) Influence of foliation on excavation stability at Rampura Agucha underground mine P Yadav, A Panda, M Sonam, B Banerjee, S Parihar, and DC Paneri Geotechnical
More informationGround Support in Mining and Underground Construction
Ground Support in Mining and Underground Construction Proceedings of the Fifth International Symposium on Ground Support 28-30 September 2004, Perth, Western Australia Edited by Ernesto Villaescusa Yves
More informationRock Material. Chapter 3 ROCK MATERIAL HOMOGENEITY AND INHOMOGENEITY CLASSIFICATION OF ROCK MATERIAL
Chapter 3 Rock Material In all things of nature there is something of the marvelous. Aristotle ROCK MATERIAL The term rock material refers to the intact rock within the framework of discontinuities. In
More informationRock Mechanics and Seismology Laboratory
CAEE 211 Geology Laboratory 6 Lab Date: 29 July 2016 Rock Mechanics and Seismology Laboratory Due Date: 5 August 2016 Attendance (based on signing): 30 points Submitting Student Name: Members of laboratory
More informationAnalysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis
Analysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis R.E. Hammah, T. Yacoub, B. Corkum & F. Wibowo Rocscience Inc., Toronto, Canada J.H. Curran Department of Civil Engineering
More informationExcavation method in Goushfill mine
International Journal of Engineering and Technology, 2 (3) (2013) 225-229 Science Publishing Corporation www.sciencepubco.com/index.php/ijet Excavation method in Goushfill mine Masoud Cheraghi Seifabad*,
More informationModule 6: Stresses around underground openings. 6.2 STRESSES AROUND UNDERGROUND OPENING contd.
LECTURE 0 6. STRESSES AROUND UNDERGROUND OPENING contd. CASE : When σ x = 0 For σ x = 0, the maximum tangential stress is three times the applied stress and occurs at the boundary on the X-axis that is
More informationA circular tunnel in a Mohr-Coulomb medium with an overlying fault
MAP3D VERIFICATION EXAMPLE 9 A circular tunnel in a Mohr-Coulomb medium with an overlying fault 1 Description This example involves calculating the stresses and displacements on a fault overlying a 5 m
More informationGround support modelling involving large ground deformation: Simulation of field observations Part 1
Ground Support 2016 E. Nordlund, T.H. Jones and A. Eitzenberger (eds) Ground support modelling involving large ground deformation: Simulation of field observations Part 1 D.Saiang, Luleå University of
More informationEDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d
Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi
More informationConceptual Numerical Modeling of Large-Scale Footwall Behavior at the Kiirunavaara Mine, and Implications for Deformation Monitoring
Rock Mech Rock Eng (2016) 49:943 960 DOI 10.1007/s00603-015-0750-x ORIGINAL PAPER Conceptual Numerical Modeling of Large-Scale Footwall Behavior at the Kiirunavaara Mine, and Implications for Deformation
More informationEffect of intermediate principal stresses on compressive strength of Phra Wihan sandstone
Rock Mechanics, Fuenkajorn & Phien-wej (eds) 211. ISBN 978 974 533 636 Effect of intermediate principal stresses on compressive strength of Phra Wihan sandstone T. Pobwandee & K. Fuenkajorn Geomechanics
More informationBorehole Camera And Extensometers To Study Hanging Wall Stability Case Study Using Voussoir beam - Cuiabá Mine
Rock Mechanics for Natural Resources and Infrastructure ISRM Specialized Conference 09-13 September, Goiania, Brazil CBMR/ABMS and ISRM, 2014 Borehole Camera And Extensometers To Study Hanging Wall Stability
More informationAn introduction to the Rock Mass index (RMi) and its applications
Reference: A. Palmström, www.rockmass.net An introduction to the Rock Mass index (RMi) and its applications by Arild Palmström, Ph.D. 1 Introduction Construction materials commonly used in civil engineering
More informationChapter 4 STRESSES AROUND BOREHOLE:BOREHOLE INSTABILITY CRITERIA
Chapter 4 STRESSES AROUND BOREHOLE:BOREHOLE INSTABILITY CRITERIA 4.1 NEAR WELLBORE STRESS- STATE Before drilling, rock stress is described by the in-situ stresses and effective stresses. As the hole is
More informationThree-dimensional settlement analysis of a primary crusher station at a copper mine in Chile
Three-dimensional settlement analysis of a primary crusher station at a copper mine in Chile B. Méndez Rizzo Associates Chile S.A., Santiago, Chile D. Rivera Rizzo Associates Inc., Pittsburgh, PA, USA
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
More informationSeismic Analyses of Concrete Gravity Dam with 3D Full Dam Model
Seismic Analyses of Concrete Gravity Dam with 3D Full Dam Model Haibo Wang, Deyu Li & Huichen Yang China Institute of Water Resources and Hydropower Research, Beijing, China SUMMARY: Seismic analyses of
More informationSTRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION TESTS
Studia Geotechnica et Mechanica, Vol. 37, No. 1, 2015 DOI: 10.1515/sgem-2015-0003 STRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION
More informationPRACTICAL THREE-DIMENSIONAL EFFECTIVE STRESS ANALYSIS CONSIDERING CYCLIC MOBILITY BEHAVIOR
PRACTICAL THREE-DIMENSIONAL EFFECTIVE STRESS ANALYSIS CONSIDERING CYCLIC MOBILITY BEHAVIOR Hiroyuki Yoshida 1, Kohji Tokimatsu 2, Tatsuya Sugiyama 3 and Tadahiko Shiomi 4 1 Member, Arch. & Struct. Eng.
More informationArticle Discrete Element Modelling of Footwall Rock Mass Damage Induced by Sub-Level Caving at the Kiirunavaara Mine
Article Discrete Element Modelling of Footwall Rock Mass Damage Induced by Sub-Level Caving at the Kiirunavaara Mine Mikael Svartsjaern * and David Saiang Department of Civil, Environmental and Natural
More informationCritical Borehole Orientations Rock Mechanics Aspects
Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds
More informationCase study: Geotechnical design of bored piles in rock for the Kawarau Falls Bridge
Case study: Geotechnical design of bored piles in rock for the Kawarau Falls Bridge Proc. 20 th NZGS Geotechnical Symposium. Eds. GJ Alexander & CY Chin, Napier Case study: Geotechnical design of bored
More informationGeotechnical & Mining Engineering Services
Geotechnical & Mining Engineering Services Southwest Research Institute San Antonio, Texas A s an independent, nonprofit research and development organization, Southwest Research Institute (SwRI ) uses
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationNumerical models on anisotropy of rocks
NGM 206 Reykjavik Proceedings of the 7 th Nordic Geotechnical Meeting Challenges in Nordic Geotechnic 25 th 28 th of May Numerical models on anisotropy of rocks Henok M. Kassa Statoil ASA, Norway, hmka@statoil.com
More informationCentrifuge Shaking Table Tests and FEM Analyses of RC Pile Foundation and Underground Structure
Centrifuge Shaking Table s and FEM Analyses of RC Pile Foundation and Underground Structure Kenji Yonezawa Obayashi Corporation, Tokyo, Japan. Takuya Anabuki Obayashi Corporation, Tokyo, Japan. Shunichi
More informationNumerical Approach to Predict the Strength of St. Peter Sandstone Pillars acted upon by Vertical Loads A case study at Clayton, IA, USA.
IOSR Journal of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vol. 05, Issue 01 (January. 2015), V2 PP 36-41 www.iosrjen.org Numerical Approach to Predict the Strength of St. Peter Sandstone
More informationIn situ fracturing mechanics stress measurements to improve underground quarry stability analyses
In situ fracturing mechanics stress measurements to improve underground quarry stability analyses Anna M. Ferrero, Maria R. Migliazza, Andrea Segalini University of Parma, Italy Gian P. Giani University
More informationSeismic stability safety evaluation of gravity dam with shear strength reduction method
Water Science and Engineering, 2009, 2(2): 52-60 doi:10.3882/j.issn.1674-2370.2009.02.006 http://kkb.hhu.edu.cn e-mail: wse@hhu.edu.cn Seismic stability safety evaluation of gravity dam with shear strength
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More information16. Mining-induced surface subsidence
16. Mining-induced surface subsidence 16.1 Types and effects of mining-induced subsidence Subsidence - Lowering of the ground surface following underground extraction of an orebody. - Types: continuous
More informationDynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models
Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global
More informationNumerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach
Numerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach S. Stefanizzi GEODATA SpA, Turin, Italy G. Barla Department of Structural and Geotechnical Engineering,
More informationRock slope failure along non persistent joints insights from fracture mechanics approach
Rock slope failure along non persistent joints insights from fracture mechanics approach Louis N.Y. Wong PhD(MIT), BSc(HKU) Assistant Professor and Assistant Chair (Academic) Nanyang Technological University,
More informationSome Aspects on the Design of Near Surface. Tunnels - Theory and Practice. Thomas Marcher¹, Taner Aydogmus²
Mitteilungen für Ingenieurgeologie und Geomechanik Band 10 6 th Colloquium Rock Mechanics Theory and Practice p. 179-188 Wien 2013 Some Aspects on the Design of Near Surface Tunnels - Theory and Practice
More informationMethods of Interpreting Ground Stress Based on Underground Stress Measurements and Numerical Modelling
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2006 Methods of Interpreting Ground Stress Based on Underground Stress Measurements and
More information1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION 1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION
Module 6 Lecture 40 Evaluation of Soil Settlement - 6 Topics 1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION 1.5.1 Definition of Stress Path 1.5. Stress and Strain Path for Consolidated Undrained Undrained
More informationSuccessful Construction of a Complex 3D Excavation Using 2D and 3D Modelling
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2015 Successful Construction of a Complex 3D Excavation Using 2D and 3D Modelling Yvette
More informationCURRENT METHODOLOGY AT THE BUREAU OF RECLAMATION FOR THE NONLINEAR ANALYSES OF ARCH DAMS USING EXPLICIT FINITE ELEMENT TECHNIQUES
CURRENT METHODOLOGY AT THE BUREAU OF RECLAMATION FOR THE NONLINEAR ANALYSES OF ARCH DAMS USING EXPLICIT FINITE ELEMENT TECHNIQUES Barbara Mills-Bria, P.E., 1 Larry Nuss, P.E. 2 and Dr. Anil Chopra 2 1
More informationVerification of the Hyperbolic Soil Model by Triaxial Test Simulations
1 Introduction Verification of the Hyperbolic Soil Model by Triaxial Test Simulations This example simulates a series of triaxial tests that can be used to verify that the Hyperbolic constitutive model
More informationTectonics. Lecture 12 Earthquake Faulting GNH7/GG09/GEOL4002 EARTHQUAKE SEISMOLOGY AND EARTHQUAKE HAZARD
Tectonics Lecture 12 Earthquake Faulting Plane strain 3 Strain occurs only in a plane. In the third direction strain is zero. 1 ε 2 = 0 3 2 Assumption of plane strain for faulting e.g., reverse fault:
More informationCOMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM
COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM Sujan MALLA 1 ABSTRACT The seismic safety of the 147 m high Gigerwald arch dam in Switzerland was assessed for
More informationRecent Research on EPS Geofoam Seismic Buffers. Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada
Recent Research on EPS Geofoam Seismic Buffers Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada What is a wall (SEISMIC) buffer? A compressible inclusion placed between
More informationEXAMINATION PAPER MEMORANDUM
EXAMINATION PAPER MEMORANDUM SUBJECT: CERTIFICATE IN ROCK MECHANICS PAPER 3.1 : HARD ROCK TABULAR EXAMINER: PJ LE ROUX SUBJECT CODE: COMRMC EXAMINATION DATE: MAY 2015 TIME: MODERATOR: WM BESTER TOTAL MARKS:
More informationSimulation of the cutting action of a single PDC cutter using DEM
Petroleum and Mineral Resources 143 Simulation of the cutting action of a single PDC cutter using DEM B. Joodi, M. Sarmadivaleh, V. Rasouli & A. Nabipour Department of Petroleum Engineering, Curtin University,
More informationDeformability Modulus of Jointed Rocks, Limitation of Empirical Methods and Introducing a New Analytical Approach
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2016 Deformability Modulus of Jointed Rocks, Limitation of Empirical Methods and Introducing
More informationTHEME A. Analysis of the elastic behaviour of La Aceña arch-gravity dam
THEME A Analysis of the elastic behaviour of La Aceña arch-gravity dam Gjorgi KOKALANOV, Professor, Faculty of Civil Eng., Skopje, Republic of Macedonia Ljubomir TANČEV, Professor, Faculty of Civil Eng.,
More informationCHAPTER 8 ANALYSES OF THE LATERAL LOAD TESTS AT THE ROUTE 351 BRIDGE
CHAPTER ANALYSES OF THE LATERAL LOAD TESTS AT THE ROUTE 351 BRIDGE.1 INTRODUCTION An important objective of this research is to determine whether accurate analyses of the lateral load-deflection behavior
More informationTIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL
TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq
More informationSeismic Response Analysis of Structure Supported by Piles Subjected to Very Large Earthquake Based on 3D-FEM
Seismic Response Analysis of Structure Supported by Piles Subjected to Very Large Earthquake Based on 3D-FEM *Hisatoshi Kashiwa 1) and Yuji Miyamoto 2) 1), 2) Dept. of Architectural Engineering Division
More informationDynamic Response of EPS Blocks /soil Sandwiched Wall/embankment
Proc. of Second China-Japan Joint Symposium on Recent Development of Theory and Practice in Geotechnology, Hong Kong, China Dynamic Response of EPS Blocks /soil Sandwiched Wall/embankment J. C. Chai 1
More informationFinite difference modelling in underground coal mine roadway
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2017 Finite difference modelling in underground coal mine roadway Ali Akbar Sahebi University
More informationModule 9 : Foundation on rocks. Content
FOUNDATION ON ROCKS Content 9.1 INTRODUCTION 9.2 FOUNDATION TYPES ON ROCKS 9.3 BEARING CAPCITY- SHALLOW FOUNDATION 9.3.1 Ultimate bearing capacity 9.3.2 Safe bearing pressure 9.3.3 Estimation of bearing
More informationDERIVATIVE OF STRESS STRAIN, DEVIATORIC STRESS AND UNDRAINED COHESION MODELS BASED ON SOIL MODULUS OF COHESIVE SOILS
International Journal of Civil Engineering and Technology (IJCIET) Volume 6, Issue 7, Jul 2015, pp. 34-43, Article ID: IJCIET_06_07_005 Available online at http://www.iaeme.com/ijciet/issues.asp?jtypeijciet&vtype=6&itype=7
More informationSYLLABUS AND REFERENCES FOR THE STRATA CONTROL CERTIFICATE. METALLIFEROUS MINING OPTION Updated November 1998
CHAMBER OF MINES OF SOUTH AFRICA SYLLABUS AND REFERENCES FOR THE STRATA CONTROL CERTIFICATE METALLIFEROUS MINING OPTION Updated November 1998 1 PART 1 : THEORY 1.1 Basic principles of rock engineering
More informationINFLUENCE OF LOCAL PERTURBATION ON REGIONAL STRESS AND ITS IMPACT ON THE DESIGN OF MAJOR UNDERGROUND STRUCTURE IN HYDROELECTRIC PROJECT
INFLUENCE OF LOCAL PERTURBATION ON REGIONAL STRESS AND ITS IMPACT ON THE DESIGN OF MAJOR UNDERGROUND STRUCTURE IN HYDROELECTRIC PROJECT *D.S. Subrahmanyam National Institute of Rock Mechanics, Bangalore
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More informationA Document on ASSESSMENT OF ENGINEERING PROPERTIES OF ROCK-MASS IN CIVIL ENGINEERING APPLICATIONS. Compiled by Mahendra Singh
A Document on ASSESSMENT OF ENGINEERING PROPERTIES OF ROCK-MASS IN CIVIL ENGINEERING APPLICATIONS Compiled by Mahendra Singh Contributors: T. Ramamurthy Bhawani Singh M.N. Viladkar K. S. Rao T.G. Sitharam
More informationWeak Rock - Controlling Ground Deformations
EOSC 547: Tunnelling & Underground Design Topic 7: Ground Characteristic & Support Reaction Curves 1 of 35 Tunnelling Grad Class (2014) Dr. Erik Eberhardt Weak Rock - Controlling Ground Deformations To
More informationDETERMINING THE STRESS PATTERN IN THE HH RAILROAD TIES DUE TO DYNAMIC LOADS 1
PERIODICA POLYTECHNICA SER. CIV. ENG. VOL. 46, NO. 1, PP. 125 148 (2002) DETERMINING THE STRESS PATTERN IN THE HH RAILROAD TIES DUE TO DYNAMIC LOADS 1 Nándor LIEGNER Department of Highway and Railway Engineering
More informationDetermination of stope geometry in jointed rock mass at Pongkor Underground Gold Mine
Volume 5, Number 2, April 2009, pp.63-68 [TECHNICAL NOTES] Determination of stope geometry in jointed rock mass at Pongkor Underground Gold Mine Budi SULISTIANTO *, M. Safrudin SULAIMAN **, Ridho Kresna
More informationInterpretation of Pile Integrity Test (PIT) Results
Annual Transactions of IESL, pp. 78-84, 26 The Institution of Engineers, Sri Lanka Interpretation of Pile Integrity Test (PIT) Results H. S. Thilakasiri Abstract: A defect present in a pile will severely
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationNumerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation
Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation P. Sushma Ph D Scholar, Earthquake Engineering Research Center, IIIT Hyderabad, Gachbowli,
More informationTWO DIMENSIONAL MODELING AND STABILITY ANALYSIS OF SLOPES OVERLAYING TO SHAHID RAGAEE POWER PLANT
4 th International Conference on Earthquake Geotechnical Engineering June 25-28, 2007 Paper No. 1637 TWO DIMENSIONAL MODELING AND STABILITY ANALYSIS OF SLOPES OVERLAYING TO SHAHID RAGAEE POWER PLANT Mohammad
More informationEUROCODE EN SEISMIC DESIGN OF BRIDGES
Brussels, 18-20 February 2008 Dissemination of information workshop 1 EUROCODE EN1998-2 SEISMIC DESIGN OF BRIDGES Basil Kolias Basic Requirements Brussels, 18-20 February 2008 Dissemination of information
More information1 Introduction. Abstract
Abstract This paper presents a three-dimensional numerical model for analysing via finite element method (FEM) the mechanized tunneling in urban areas. The numerical model is meant to represent the typical
More information