Scale effects on strength of geomaterials, case study: coal
|
|
- Roderick Collins
- 6 years ago
- Views:
Transcription
1 1 ELSEVIER Scale effects on strength of geomaterials, case study: coal Luc Scholtès, Frédéric-Victor Donzé *, Manoj Khanal Earth Science and Resource Engineering, QCAT, CSIRO, Pullenvale, Australia Elsevier use only: Received date here; revised date here; accepted date here Abstract Scale effects on the strength of coal are studied using a Discrete Element Model. The key point of the model is its capability to discriminate between the strictly sample size effect and the Discrete Fracture Network (DFN) density effect on the mechanical response. True triaxial compressive test simulations are carried out to identify their respective roles. The possible bias due to the discretization size distribution of the discrete element model has been investigated in detail by considering low-resolution models. The model is shown to be capable of quantitatively reproducing the dependency of the maximum strength on the size of the sample. This evolution mainly relies on the DFN density. For all given sizes, as long as the DFN density remains constant or if discontinuities are absent, the maximum strength of the sample does not change Elsevier Science. All rights reserved Keywords: Scale effect; Strength properties; Discrete Fracture Network; Fracture density; Discrete Element Method. 1. Introduction Discontinuities are observed at all scale in coal and they greatly affect its properties. For instance, when the density of discontinuities increases, a decrease of the strength is observed. The unconfined compressive strength (UCS) of coal is also known to be dependent on the scale at which it is tested [1, 2, 3, 17]. A non-uniform distribution of discontinuities can also induce en important strength anisotropy within the rock mass [4, 5, 6]. Coal, like rock, is a cohesive frictional geomaterial structure of which usually contains complex networks of discontinuities (cleats). Cleat networks can present very different spatial distributions which depend on the coalification process, ranging from highly organized to disordered configurations. Generally, an organized cleat set can be decomposed into face and orthogonal butt cleats * Corresponding author. Tel.: ; frederic.donze@csiro.au. (see Figure 1). Face cleats are considered as extensional fractures oriented perpendicular to the least compressive principal stress applied on the coal during coalification, whereas butt cleats occur as non persistent fractures, usually perpendicular to face cleats, resulting from the relaxation of the original stress field. Face cleats Butt cleats Figure 1. An illustration of face and butt cleats in coal (after [7]).
2 2 In common bituminous coal, bright coal (vitrain rich) and dull coal (durain rich) are generally the main constituents. In bright coal, cleat spacing occurs at the millimeter scale, whereas it is generally much greater in dull coal (decimeter scale). Based upon the observation that cleating is often confined to the brighter bands in coal, the use of the terms bright and dull infers a measure of volumetric cleat density. The brightness profiles, therefore, enable an assessment to be made of coal seams structure [3]. This variation in density could also explain why dull coal is observed to have a greater strength than bright coal [3]. Master cleats, which are defined as vertically cutting both vitrain and durain layers, have been observed even more widely spaced than the other classes of cleat [4, 7]. In the present study, a numerical model based on the Discrete Element Method (DEM) [8] is used to investigate the contribution of discontinuities on strength properties of a generic bituminous coal. Cleats are represented as discontinuity planes (also referred to here as fracture or joint planes) for which both spatial distribution and mechanical properties are controlled. Simulations of true triaxial compression tests are performed on numerical samples. In order to verify the validity of the method, the predicted results are systematically compared to the experimental ones obtained by Medhurst and Brown on coal [3]. First, the study focuses on intact media, in order to assess the possible effect of discretization, that is, the resolution of the model, on the results. Then, the effect of an increasing density of discontinuities is characterized. Previous relationships between the strength and the sizescale have been proposed. For example, Evans and Pomeroy [21] found that the average compressive strength σ c of unconfined cubes of coal depend on the cube side β length a according to an exponential law σ c = kα, where β is a constant depending on the type of coal. Hoeck and Brown [15] suggested that in the case of an intact rock specimen, its strength can be related to a reference sample strength of 50 mm diameter by the formula σ ( ) c = σ 50 D 50, with D being the section diameter of the sample. For jointed rocks, an empirical formula, using a geological strength index (GSI), has been set up. However, these formulae do not reflect in a quantitative way the overall relationships between the size-scale, the fracture density and the material strength. In this paper, a quantitative investigation on the scale effects on strength is proposed. The size of the sample is increased and the evolution of its strength is studied. The scale dependency of the strength is then discussed. 2. Formulation of the model The algorithm used in the present Discrete (or Distinct) Element Method (DEM) [8], involves two steps. First, based on constitutive laws, interaction forces between Discrete Elements (DE) are computed. Second, Newton second law is applied to determine, for each DE, the resulting acceleration, which is then time integrated to find its new position. This process is repeated until the simulation is finished. All developments and simulations described in this paper have been implemented into the YADE Open DEM platform [9] Constitutive law For small deformations, coal exhibits a linear elastic response characteristic of geomaterials. To reproduce this behavior, linear elastic interaction forces between the discrete elements are sufficient. In the present model, the interaction force F which represents the action of DE a on DE b may be decomposed into a normal force F n and a shear force F s which are related to relative normal and incremental shear displacements through the stiffnesses K n and K s, in the normal and the tangential directions respectively. Once the DE s packing has been generated to form the numerical sample, pairs of initially interacting Discrete Elements are identified within an interaction range γ int, such that, D = γ ( R + R ) (1) eq int a b where D eq is the initial equilibrium distance, R a and R b the radius of elements a and b respectively, and γ int 1 [19]. Contrary to classical DEM which only considers strictly contacting DE (i.e. γ int = 1 ), this feature provides the possibility of controlling the initial number of interacting links (also referred here to as the coordination number), whatever the granular packing, by simply varying γ int before the first timestep of the computation cycle [10, 18]. The normal interaction force is calculated through the local constitutive law (see Figure 2) and can be split into two parts, the compressive and the tensile components. In compression, F n is linear and it is given by, Fn = K n ( D Deq ) (2) where D is the distance between the centroids and K n the normal stiffness.
3 3 F s F n A int C φ b D rupture 1 K n D φ c F n K n ζ 1 TA int F n,max Figure 3. Rupture criterion used in the model. Figure 2. Normal interaction force between DE. In tension, the normal interaction force is also linear using the same stiffness for small deformations. The maximum acceptable tensile force F n,max is defined as a function of the tensile strength T such that, F = TA (3) n,max int where A int = π(min(r a,r b )) 2 is the interacting surface between a and b. After the tensile force reaches its maximum value, the stiffness may be modified by a softening factor ζ to control the energy released by the link breakage: Kn F n = ( D Drupture) (4) ζ When D > D rupture, rupture occurs and interaction forces are set to zero. The shear force F s is calculated by updating its orientation which depends on the orientation of the direction passing through the two centroids of the interacting DE, and by adding the increment of shear force F s as defined by Hart et al. [11]: F = K u (5) s s s with K s being the shear stiffness and u s the relative incremental tangential displacement. The maximum shear force F s,max is characterized by the normal force F n, the cohesion C, the local frictional angle φ b and the local residual frictional angle φ c. The maximum shear force is calculated for a bonded interaction according to, Fs, max = Fn tanϕ b + CAint (6) Only purely frictional new contacts can appear during the simulation, with a maximum shear force defined as, F tanϕ s, max = F n (7) c Since spherical discrete elements are used in this model, the interaction law must include a moment in order to control the rolling occurring between DE during shear displacements. The moment transferred between two interacting DE can be written as, M elast = Σθ K (8) r r where θ r is the relative rotation angle and K r is the rotational stiffness defined by β.k s.r avg [12], where β is a dimensionless coefficient and R avg the average radius of the two interacting DE. An elastic limit is introduced (Figure 4) and when it is reached, the plastic moment M plast defined by M ηf R plast = (9) n avg takes place, η being a dimensionless factor. To model the nonlinear behavior of geomaterials, a modified Mohr-Coulomb model is used (Figure 3).
4 4 n j t j M r M plast K r ηf n R avg n c t c 1 θ r Figure 4. Rolling moment law considered between interacting DE. n j t j Doing so, the sliding process increases and the resulting macroscopic internal friction angle Φ can reach values corresponding to those measured for geomaterials [12] Modeling Discontinuities In order to model discontinuities in the DEM, a specific contact formulation, inspired by the work of Cundall and co-authors on the development of the Synthetic Rock Mass model [13], has been implemented into YADE Open DEM. This joint contact logic is based on the identification of each interaction crossing the plane which represents the discontinuity plane (see Figure 5). Normal and tangential directions ( n r j, t r j ) of these interactions are updated to new orientations ( n r c, t r c ) which depend on the orientation of the joint. With this formulation, it becomes possible to get rid of the interlocking process (due to the spherical geometry of the DE) and of the discretization of the joint surface occurring during shearing displacements [14]. The interactions forces concerned by this modified contact logic use the same formulation as defined in section 2.1, except that they do not transmit any moments. Their mechanical properties (K n, K s, C, T, φ b and φ c ) can be chosen equal to or different from the ones used to represent the intact material. In addition, a dilatation component has been added to the normal force, in order to simulate the aperture of the joint, usually observed with shearing [14]. When sliding occurs between two DE on both sides of a discontinuity, F n is updated such that, F n = K D + u tanψ (10) n s where ψ is the dilatation angle. All the simulations performed in this paper use this joint contact logic when discontinuities are taken into account. (c) Figure 5. Joint contact logic: identification of the interactions crossing the discontinuity plane, DE are colored depending on the side of the plane to which they belong. Contact orientation and displacements induced by shearing displacement: classic contact orientation and modified contact orientation (c). 3. Model setup Within the proposed DEM, intact coal material is represented by a set of spherical DE. In this study, numerical samples are generated using a polydisperse distribution size for the DE in order to avoid any artifacts related to regular packings. A wide range of coal properties can be found in the literature depending on site or testing specificities (see for example [2, 3]). Since the experimental results obtained by Medhurst and Brown [3] provide an often referenced data base, they are used here to setup the mechanical properties of the model (see Table 1). Due to the intrinsic presence of cleats at all scales in coal, it is difficult to assess directly its intact mechanical properties from triaxial tests, even if coal specimens are very small. Measurements made on small dull coal specimens by Medhurst and Brown [3] are thus used to provide reference properties for the intact coal s mechanical properties. Average values of 4 GPa for Young s modulus E, 0.25 for Poisson s ratio ν, 43 for the internal friction angle Φ and 7 MPa for the cohesion C have been estimated by [3] based on the generalized Hoek and Brown criterion parameters [15]. Corresponding local properties used in the DEM to reproduce this generic coal are summarized in Table 2.
5 5 Table 1. Mohr-Coulomb properties for generic coal. Mechanical properties Values Young s modulus E (GPa) 4 Poisson s ratio ν (-) 0.25 Internal friction angle Φ ( ) 48 Cohesion C (MPa) 7 moved vertically at a constant deformation rate while keeping the confining pressure P by adjusting lateral wall positions. Strain and stress variations are computed through wall displacements and resultant forces on walls respectively. ε 1 varies Table 2. Input parameters used to model generic coal. Input parameters Values Normal stiffness K n (MN.m -2 ) Tangential stiffness K s (MN.m -2 ) 9.66 Local friction angle φ= φ b =φc ( ) 40 Tensile strength T (MPa) 11.5 Local cohesion C (MPa) 115 Coefficient of rolling stiffness β (-) 1.0 Coefficient of rolling elastic limit η (-) 0.15 Coefficient of damaging law in release mode (-) 0.0 Damping coefficient η 0.4 Mechanical characteristics of cleats are not straight forward to determine because, in vitrain layers for example, both mineralized and unmineralized cleats can be found in the same seam. A local strength value close to the average tensile strength for coal is thus chosen here (see Table 3). The friction angle is chosen to be equal to 30 whereas cohesion C and tensile strength T are taken to be a tenth of the intact ones. In addition, the dilatation angle ψ is set to 5. Table 3. Input parameters for discontinuities in the DEM. Input parameters for discontinuities Friction angle φ ( ) 30 Values Tensile strength T (MPa) 1.15 Local cohesion C (MPa) 11.5 Dilation angle ψ ( ) 5 4. Triaxial loading tests A series of numerical tests have been performed simulating true triaxial testing apparatus. The loading is applied on the numerical specimen through six rigid frictionless boundary walls with servo-controlled positions depending on the chosen loading path (see Figure 6). The numerical samples are subjected to an isotropic compression up to a predefined confining pressure P. Once this confining pressure is reached, top and bottom plates are σ 2 = σ 3 =P Figure 6. Triaxial testing configuration. σ 2 = σ 3 =P Like the number of elements, the density or the strain rate applied on the specimen need to be selected correctly [12]. The density of the model is chosen to be equal to an average value for coal material (1450 kg.m -3 ), while the strain rateε& 1 is set sufficiently low to limit the inertial contribution of the DE. Quasi-static conditions are systematically verified by checking that dividing ε& 1 by any factor greater than 1 does not produce significant changes in the mechanical response. Since the scope of this article is the scale effect on strength, it is crucial to first identify the possible contribution of the discretization of DE sets on the response of an intact numerical specimen. For this purpose, four different cubical samples with identical volumes V=L 3 are considered. They are made of different numbers of DE with different average diameters d but the distribution size ratio is kept the same, approximately 2 (see Figure 7). Each specimen is generated according to the same procedure, which consists in growing particle s radii iteratively starting from an initial set of non overlapping DE contained within a box of volume V, until the pressure applied on the walls has reached the predefined value P. The four samples are made of a number of DE varying from 400 to , resulting in different L/d ratio, ranging from 7 to 37. Figure 7. Four cubic samples of identical volume V=L 3, composed with discrete elements of different average diameters d, such that L/d = 7, 15, 22, 37.
6 6 The behaviors of the four samples subjected to triaxial compressions under a 200 kpa confining pressure are presented in Figure 8. It can be observed that the stressstrain responses as well as the volumetric deformations are different depending on the discretization size. In addition to the less pronounced variations in the sample response around the peak stress, the strength, as well as Young s modulus, tends to increase when the L/d ratio increases. discontinuities. In that case, the local properties of the fabric may be altered. To get rid of these resolution effects on Young s modulus and the strength when the discretization is low, the local values of the stiffness and tensile strength could be modified. Another alternative without modifying the micromechanical parameters is to control the initial interacting distance for which links between particles are created (see section 2.1). The interaction factor γ int has therefore been modified in order to keep the initial coordination number constant for all samples. The results are presented in Figure 9. (c) (d) Figure 8. The evolution of the differential stress q, volumetric deformation dv/v, porosity (c) and coordination number (d) during triaxial compressions performed for 200 kpa confinement on four cubic samples composed of 400, 3200, and DE respectively. (c) (d) Yet, local interaction laws have been carefully implemented to be size independent and the respective fabric of each sample should be equivalent since they have been generated using the same procedure. One could therefore expect a discretization independency. The porosities and, more particularly, the initial values of the coordination number (average number of interacting links per element) of each sample, are specific to each assembly and thus show an evolution of the internal fabric which contributes to the global response discrepancy. This evolution can be formulated using the ratio between the number of DE which composes the inner part of the sample N int and the number of DE in contact with boundary walls N bound. At the considered discretization levels, N bound /N int varies dramatically (from 0.7 to 0.16). Given an identical volume V, fewer DE induces, proportionally, more frictional contacts with boundaries than cohesive links within the whole specimen, resulting in a lower Young s modulus as well as a weaker strength at the scale of the assembly. Nevertheless, it is worth noting that the differences tend to decrease for higher values of L/d. However, it was of utmost importance in this paper to show the consequences on the mechanical response of using low resolution models. This needs to be taken into consideration when representing highly discontinuous media which may then have a low resolution in between Figure 9. Evolution of the differential stress q, volumetric deformation dv/v, porosity (c) and coordination number (d) during triaxial compressions performed under 200 kpa confinement on four cubic samples composed of 400, 3200, and DE respectively with an identical initial coordination number (d). Contrary to the case where initial coordination numbers were not kept constant (Figure 8), Young s modulus and the differential peak strength converge to the same value whatever the L/d ratio is (except for the low discretization case which is strongly dependent on the spatial distribution of the DE). Thus, the coordination number can be used as a key parameter in DEM simulations when low discretization levels are needed, for example, to keep a low calculation cost. This ratio can counter balance boundary effects when L/d < Contribution of the density of pre-existing discontinuities As the purpose of the present study is to deal with the effect of discontinuities on the overall mechanical properties of cohesive frictional materials, several series of triaxial compression tests have been performed on numerical samples containing a pre-existing set of
7 7 discontinuities (i.e. a Discrete Fracture Network, DFN). The first of these tests concerns the effect of the DFN density on the mechanical properties of a sample with a given size and discretization. In a numerical sample of L = 60 mm with a discretization ratio of L/d=22, different sets of discontinuities have been initially introduced (see section 2.2 for implementation details). In order to represent coal in a realistic but simple manner, persistent discontinuities are are positioned in two directions (bedding planes and face cleats respectively perpendicular to the Y and Z axis), whereas non persistent discontinuities (butt cleats) are placed in the third direction (perpendicular to the X axis) (see Figure 10). A constant ratio between cleat spacing and cleat height is chosen here based upon observations done by Dawson and Esterle [7] on Australian Permian coal. Y X Z Figure 10. Four cubic samples of identical size (L=60 mm) and discretization (L/d=22) with an increase of the density of discontinuities. Four different samples have been built up with an increasing fracture density which can be characterized by the number of cleats per sample length L (from 0 to 4 cleats per length, see Figure 10). Triaxial compressions as described previously have then been performed at several confining pressures (along the vertical Y axis) and model predictions are compared to the Medhurst and Brown results obtained on same size coal samples with different brightness levels [3]. For clarity reasons and since the experimental results have also shown a good agreement with the generalized Hoek and Brown Criterion (HBC) [15], the present numerical results are also compared in terms of this criterion which is defined as: a complement to the HB Criterion. According to the GSI, the strength of rock type materials decreases with the increase of the fracture density or blockiness of the considered rock mass, which is the case here. In addition, it is remarkable that, despite its relative simplicity, the idealized cleat network defined here seems to be relevant to model the change of brightness in coal, at least for the considered case. Table 4. Hoek and Brown Criterion parameters for both experimental [3] and DEM results which show the effects of cleat density. Sample Hoek and Brown parameters σ c (MPa) m s Medhurst and Brown C DEM Intact (no cleat) Medhurst and Brown - C4/C DEM - 1 cleat/l Medhurst and Brown - C3/C2/C DEM - 2 cleats/l Medhurst and Brown - C DEM - 4 cleats/l Figure 12 shows the stress-strain responses as well as the volumetric deformations predicted by the DEM for triaxial compressions performed on all four specimens under a 200 kpa confining pressure. 0.5 mσ 3 1 = σ3 + σc + s σ c σ (11) where σ c is the uniaxial compressive strength of the intact material, and m and s are chosen to fit Medhurst results. The HBC parameter values for both experimental and DEM results are presented in Table 4, whereas the comparison between the experimental data and the numerical results are presented in Figure 11. Numerical results show a decrease of the strength with an increase in the DFN density which is in good agreement with experimental data. This result also agrees with the Geological Strength Index (GSI) approach which is commonly used in rock and mining engineering [15, 16] as Figure 11. Medhurst and Brown [3] experimental results obtained on coal samples with different brightness levels compared to DEM predictions for an increasing fracture density in terms of the generalized Hoek and Brown criterion [15].
8 8 One can note that Young s modulus and the dilatation angle are equally affected by the presence of discontinuities. Both of them tend to decrease with an increase in material cleating, which also confirms experimental observations [3]. The presence of pre-existing discontinuities induces modifications into the fabric of the material which produce dramatic changes at the sample scale. To get more insights on the microscopic mechanisms (i.e. local mechanisms here) induced by cleating, it is possible to compare the deformation and fracturing processes which occur during loading for both, an intact and a discontinuous samples (see Figures 13 and 14). Figure 14. Displacement fields for both intact sample (left) and sample with pre-existing discontinuities (right): global view (top) and Y view (bottom). Figure 12. Effect of cleat density. Differential stress q, volumetric deformation dv/v, obtained during triaxial compressions performed under 200 kpa confinement on four specimens of identical size and discretization but different cleat density. Figure 13. Crack distribution for both intact sample (left) and sample with pre-existing discontinuities (right): global view (top) and Y view (bottom). Obviously, deformations as well as propagation of cracks (i.e. when a local interaction link breaks) are strongly driven by the presence of pre-existing discontinuities which act as planes of weakness inside the material. Due to their lower cohesive properties, the interacting links crossing the discontinuities tend to fail before the ones located in the intact zones of the material. It should also be noted that, as for the intact part of the material, micro-cracking mainly occurs in the direction parallel to the main applied stress (the Y axis here), which is also in agreement with experimental observations made on rocks as long as the confining pressure remains sufficiently low. In order to assess once again, the influence of discretization on DEM predictions, the same series of triaxial compressions on samples have been performed on samples with different L/d ratios while keeping the same networks of discontinuities (Figure 15). As expected, a discretization effect can be observed. The amplitude of the strength decrease is not constant, with a higher strength decrease for the lower L/d values. As noted in section 5, in the case of poor discretization levels, the properties of the intact part of the material can vary dramatically. Here, with the presence of discontinuities, this effect is reported at the scale of each intact cell separated by the discontinuities. Thus, for different L/d, values, each cell has a different fabric, inducing a variation of the global strength of the sample. Thus, attention must be paid when a network of discontinuities is introduced into a DEM and comparisons need to be made between results using different discretization sizes.
9 9 Maximum strength (MPa) Figure 15. Maximum strength values obtained from triaxial compressions test performed on samples with different discretization size, 3200, 10800, 5000 ED (i.e. L/d = 15, 22 and 37 respectively) at 200 kpa confinement. 6. Scale effect on strength Fracture density [m -1 ] A test configuration has been set up to assess whether or not the DEM model can reproduce the scale effect on strength, by comparing its predictions to the Medhurst and Brown results [3] obtained on coal samples of diameters ranging from 60 to 300 mm. Using a given fixed distribution size for the DE (d=8.2 mm), four cubic samples of different edge lengths L i =60, 120, 180 and 300 mm, have been generated (Figure 16), and subjected to triaxial loadings according to the previous procedures. A slight decrease of the strength, which converges to a limit value, can be observed when the sample size increases (see also Figure 20). This slight decrease of the strength is of the same nature as the one observed in section 5. It is due to the low resolution configuration of the smallest sample. Despite this artifact, it can be concluded that the intact material does not show any scale effect regarding its strength, which is confirmed by previous observations [19, 20]. A discrete fracture network (DFN) is now defined and introduced into the numerical samples (see Figure 18). This DFN has the same characteristics as the ones used in section 6 in terms of its spatial distribution (bedding planes, face and butt cleats) and mechanical properties. The difference is that the DFN is kept constant for all the samples, i.e. the space length between the discontinuities is kept constant and equal to the edge length of the smallest sample (60 mm), which is considered to be the reference intact length. Doing so, the numerical simulation procedure is consistent with the Medhurst and Brown study, which considered the smallest sample of 60 mm length to be the intact one (see values in Table 5) when using the HB criterion to quantify scale effects for the considered case. The stress-strain and volumetric curves obtained for the tests performed under a 200 kpa confinement are presented in Figure 19. The corresponding peak strength values are plotted as a function of the sample size in Figure 20. Z Y X Z Y X Figure 16. Four intact cubic samples of different sizes L i=60, 120, 180 and 300 mm, with constant particle size distribution (d=0.82 mm). Results in terms of stress-strain responses and volumetric deformations of the intact material (i.e. no DFN inside) are presented in Figure 17. Figure 17. Differential stress q and volumetric deformation dv/v obtained during triaxial compressions performed for a confinement of 200 kpa, on four intact samples with length L i ranging from 60 to 300 mm. Figure 18. Four cubic samples of different sizes with constant Discrete Element size distribution and constant DFN (i.e. same spacing between discontinuities). Table 5. Hoek and Brown criterion parameters for Moura coal depending on sample size, from [3]. Sample size Hoek and Brown parameters σ c (MPa) m s 60 mm mm mm mm
10 10 Figure 19. Differential stress q and volumetric deformation dv/v evolutions obtained during triaxial compressions performed under 200 kpa confinement on four pre-fractured numerical samples with length L i ranging from 60 to 300 mm. Maximum strength (MPa) Figure 20. Scale effect on the maximum strength of the numerical model. Comparison between intact and pre-fractured samples. In the presence of discontinuities, the decrease in the maximum strength is now well pronounced and it is consistent with the experimental observations provided by Medhurst and Brown and expressed in terms of the generalized Hoek and Brown criterion (see Figure 21). A remarkable aspect concerns the contribution of the DFN density with the sample size. When representing the evolution of this density (computed here as the ratio between the total surface S of the discontinuities and the total volume V of the sample, which is also known to be the P 32 value [19]) versus the length L of the sample, it can be observed that the DFN density converges towards a limit value (see Figure 22). This convergence is associated to the convergence observed for the maximum strength value of the cleated sample (Figure 20). Figure 21. Medhurst and Brown [3] experimental results on coal samples of different sizes compared to DEM predictions for an increasing sample size in terms of the generalized Hoek and Brown criterion [15]. Thus, when the size of the sample increases, as long as the density of the DFN increases, a decrease of the maximum strength is observed. One could extrapolate that given any sample size, as soon as the density of the DFN remains constant, the maximum strength of this sample will no longer change. These results suggest that the scale effect on strength mainly relies on the evolution of the DFN density rather than directly on the size of the sample. Fracture density [m -1 ] Figure 22. Evolution of the density of fractures with the sample length. Comparable numerical results have been recently obtained by Esmaieli et al. [19]. When using disordered DFNs based on observations from a volcanic rock, their model exhibits a decrease of the sample s strength when the size of the sample increases. Depending on the DFN, for the same number of fractures, different densities of fractures (P 30 ) could be also generated. As long as the density of fractures increases, the strength of the sample decreases. An interesting point is that for the same number
11 11 of fractures, different values of P 30 could be obtained. Thus, even if the quantity of fractures increases when the size of the sample increases, a decrease of the strength is obtained only when an increase of the density P 32 occurs. These results fully corroborate our results in terms of scale effect, i.e. the preeminent role of the density of fractures. Based on the 2D method of fractal geometry provided by Zengchao et al. [20] it is suggested that the scale effect on the rock mass strength can be formulated using the fracture quantity, the fractal dimension of fracture quantity distribution and rock mass size. For 3D fracture networks, the 2D information such as the distribution of the fracture length cannot give a unique solution for the total fracture surface. Hence, eventhough the present model provides results comparable to Zengchao et al. [20], one can wonder whether or not their model suffers from the lack of a more complete 3D formulation. 7. Conclusion The effect of discontinuities on strength was numerically investigated using a Discrete Element Method. The model integrates structural effects induced by discontinuities which are known to produce strong behavior changes in rock-like materials, scale dependency above all. Triaxial test simulations were performed on intact samples and samples with pre-existing discontinuities. The results were systematically compared to experimental ones in order to estimate the reliability of the method. Emphasis was put on possible numerical bias. Resolution effects (i.e. Discretization size effects) have been shown for low resolutions which impose special care when discrimination studies and quantitative results are expected. Controlling the coordination number is one of the possible solutions to get rid of this artifact. The study focused on the relationship between the strength and the set of pre-existing discontinuities (DFN) in a numerical sample. Comparisons with experimental data have shown that the model can reproduce a strength dependency on the density of the DFN. The effect of structural features on the fabric of the material can be simulated. It has been shown that, starting from an intact material, increasing the DFN density tends to diminish the strength properties in a way quantitatively comparable to the experimental observations. When increasing the size of the numerical sample, it has been observed that the strength is kept constant when discontinuities are absent or when the DFN density (P 32 ) is not modified. As soon as the DFN density varies, the strength is modified as well. This conclusion is in good agreement with the results recently obtained by Esmaieli et al. [19], where the role of the density of fractures has been explicitly considered in the scale effects on the strength of disordered fractured rock mass. Thus, the formulation of a relationship between the distribution of fractures and the size of the rock mass would need to integrate an accurate value for P 32. Moreover, the contributions of the persistency of the fractures and their interconnectivity would also have to be quantified in this formulation. Acknowledgments Authors would like to acknowledge the Large Open Pit Project for its participation in this work. References [1] Z.T. Beniawski, The effect of specimen size on compressive strength of coal, International Journal of Rock Mechanics and Mining Sciences, 5, pp , [2] N. Deisman, D. Mas Ivars, C. Darcel, R.J. Chalaturnyk, Empirical and numerical approaches for geomechanical characterization of coal seam reservoirs, International Journal of Coal Geology, doi: /j.coal , [3] T.P. Medhurst and E.T. Brown, A study of the mechanical behaviour of coal for pillar design, International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 35(8), pp , [4] G.S. Esterhuizen, The effect of structural discontinuities on coal pillar strength as a basis for improving safety in the design of coal pillar systems, Project report, University of Pretoria, [5] A.B. Swilski, Determination of the anisotropic elastic moduli of coal, International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 21(1), pp. 3-12, [6] S. Okubo, K. Fukui and Qi Qingxin, Uniaxial compression and tension tests of anthracite and loading rate dependence of peak stregth, International Journal of Coal Geology, 68, pp , [7] G.K.W. Dawson and J.S. Esterle, Controls on coal cleat spacing, International Journal of Coal Geology, DOI: /j.coal , [8] F.V., Donzé, V. Richefeu and S.-A. Magnier, Advances in Discrete Element Method applied to Soil, Rock and Concrete Mechanics, in: State of the art of geotechnical engineering, Electronic Journal of Geotechnical Engineering, p. 44, 2009 [9] J. Kozicki and F.V. Donzé, A new open-source software developed for numerical simulations using discrete modeling methods, Computer Methods in Applied Mechanics and Engineering, 197, pp , [10] F. V. Donzé, J. Bouchez, S. A. Magnier, Modeling fractures in rock blasting, International Journal of Rock Mechanics and Mining Sciences, 34, pp [11] R. Hart, P.A. Cundall and J. Lemos, Formulation of threedimensional distinct element model Part II - mechanical calculations of a system composed of many polyhedral blocks. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 25, pp [12] J.P. Plassiard, N. Belheine and F.V. Donzé, A spherical discrete element model: calibration procedure and incremental response, Granular Matter, 11, pp , [13] M. Pierce, D. Mas Ivars, P.A. Cundall and D. Potyondy, A synthetic rock mass model for jointed rock. Proceedings of the first CA-US rock mechanics symposium, Vancouver, BC, pp , [14] D. Mas Ivars, D. Potyondy, M. Pierce and P.A. Cundall, The smoothjoint contact model. Proceedings of the eighth world congress on computational mechanics, fifth European congress on computational
12 12 methods in applied sciences and engineering, paper no. a2735, Venice, Italy, [15] E. Hoek and E.T. Brown, Practical estimates of rock mass strength. Int J Rock Mech Min Sci, 34(8), pp , [16] J.L. Justo, E. Justo and J.M. Azanon, The use of rock mass classification systems to estimate the modulus and strength of jointed rock, Rock Mech Rock Eng, 43(3), , [17] F. Heuze, Scale effects in the determination of rock mass strength and deformability, Rock Mech, 12, , [18] S. Hentz, L. Daudeville and F.V. Donzé, Identification and validation of a discrete element model for concrete, J Eng Mech, 130(6), pp , [19] K. Esmaieli, J. Hadjigeorgiou and M. Grenon, Estimating geometrical and mechanical REV based on synthetic rock mass models at Brunswick Mine, International Journal of Rock Mechanics and Mining Sciences, DOI: /j.ijrmms , online [20] F. Zengchao, Z. Yangsheng and Z. Dong, Investigating the scale effects in strength of fractured rock mass, Chaos, Solitons & Fractals, Volume 41, Issue 5, 15 September 2009, pp , [21] I. Evans and C.D.Pomeroy Strength fracture and workability of coal. London: National Coal Board; 1966.
Discrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationPROGRESSIVE FAILURE MECHANISMS IN JOINTED ROCK: INSIGHT FROM 3D DEM MODELLING
II International Conference on Particle-based Methods Fundamentals and Applications PARTICLES 2011 E. Oñate and D.R.J. Owen (Eds) PROGRESSIVE FAILURE MECHANISMS IN JOINTED ROCK: INSIGHT FROM 3D DEM MODELLING
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 information3D HM-DEM model for Hydro-Fracturing
3D HM-DEM model for Hydro-Fracturing E. Papachristos, F.V. Donzé & B. Chareyre Laboratoire Sols, Solides, Structures, Grenoble, France efthymios.papachristos@3sr-grenoble.fr,, frederic.donze@3srgrenoble.fr,
More informationNumerical Simulations of Triaxial Test with Sand Using DEM
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (2009), No. 3 4, pp. 149 171 IBW PAN, ISSN 1231 3726 Numerical Simulations of Triaxial Test with Sand Using DEM Łukasz Widuliński, Jan
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 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 informationParticle flow simulation of sand under biaxial test
5th International Conference on Civil Engineering and Transportation (ICCET 2015) Particle flow simulation of sand under biaxial test Xiao-li Dong1,2, a *,Wei-hua Zhang1,a 1 Beijing City University, China
More informationA micromechanical approach to describe internal erosion effects in soils
A micromechanical approach to describe internal erosion effects in soils Luc Scholtès, Pierre-Yves Hicher, Luc Sibille To cite this version: Luc Scholtès, Pierre-Yves Hicher, Luc Sibille. A micromechanical
More informationDEM simulation of fracture process of inherently anisotropic rock under Brazilian test condition
Title DEM simulation of fracture process of inherently anisotropic rock under Brazilian test condition Author(s) Kwok, CY; Duan, K Citation The 49th US Rock Mechanics / Geomechanics Symposium, San Francisco,
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 informationModelling dynamic fracture propagation in rock
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part B Faculty of Engineering and Information Sciences 2017 Modelling dynamic fracture propagation in
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 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 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 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 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 informationUsing the Timoshenko Beam Bond Model: Example Problem
Using the Timoshenko Beam Bond Model: Example Problem Authors: Nick J. BROWN John P. MORRISSEY Jin Y. OOI School of Engineering, University of Edinburgh Jian-Fei CHEN School of Planning, Architecture and
More informationThe Influence of Contact Friction on the Breakage Behavior of Brittle Granular Materials using DEM
The Influence of Contact Friction on the Breakage Behavior of Brittle Granular Materials using DEM *Yi-Ming Liu 1) and Hua-Bei Liu 2) 1), 2) School of Civil Engineering and Mechanics, Huazhong University
More informationTable of Contents. Foreword... xiii Introduction... xv
Foreword.... xiii Introduction.... xv Chapter 1. Controllability of Geotechnical Tests and their Relationship to the Instability of Soils... 1 Roberto NOVA 1.1. Introduction... 1 1.2. Load control... 2
More informationMEMORANDUM 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 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 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 informationPLASTICITY FOR CRUSHABLE GRANULAR MATERIALS VIA DEM
Plasticity for crushable granular materials via DEM XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J. Owen, D. Peric and M. Chiumenti
More informationNumerical simulations for isostatic and die compaction of powder by the discrete element method
Numerical simulations for isostatic and die compaction of powder by the discrete element method J.-F. Jerier 1, B. Harthong 1, V. Richefeu 1, D. Imbault 1, F.-V. Donzé 1,1, P. Dorémus 1 a University of
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 informationFluid driven cohesive crack propagation in quasi-brittle materials
Fluid driven cohesive crack propagation in quasi-brittle materials F. Barpi 1, S. Valente 2 Department of Structural and Geotechnical Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129
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 informationGeology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)
Geology 229 Engineering Geology Lecture 5 Engineering Properties of Rocks (West, Ch. 6) Common mechanic properties: Density; Elastic properties: - elastic modulii Outline of this Lecture 1. Uniaxial rock
More informationThe response of an idealized granular material to an incremental shear strain
The response of an idealized granular material to an incremental shear strain Luigi La Ragione and Vanessa Magnanimo Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Bari, Italy E-mail: l.laragione@poliba.it;
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 informationLimit analysis of brick masonry shear walls with openings under later loads by rigid block modeling
Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,
More informationPrediction of torsion shear tests based on results from triaxial compression tests
Prediction of torsion shear tests based on results from triaxial compression tests P.L. Smith 1 and N. Jones *2 1 Catholic University of America, Washington, USA 2 Geo, Lyngby, Denmark * Corresponding
More informationSamantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2
Samantha Ramirez, MSE Stress The intensity of the internal force acting on a specific plane (area) passing through a point. Δ ΔA Δ z Δ 1 2 ΔA Δ x Δ y ΔA is an infinitesimal size area with a uniform force
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 informationROCK MASS CHARACTERISATION IN ENGINEERING PRACTICE
Paul MARINOS NTUA, School of Civil Engineering, 9 Iroon Polytechniou str., Athens, 157 80, Greece, e-mail : marinos@central.ntua.gr ROCK MASS CHARACTERISATION IN ENGINEERING PRACTICE 1. INTRODUCTION The
More informationNORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.
NORMAL STRESS The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. σ = force/area = P/A where σ = the normal stress P = the centric
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 informationNumerical Simulation of Unsaturated Infilled Joints in Shear
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2018 Numerical Simulation of Unsaturated Infilled Joints in Shear Libin Gong University
More informationExperimental study of mechanical and thermal damage in crystalline hard rock
Experimental study of mechanical and thermal damage in crystalline hard rock Mohammad Keshavarz Réunion Technique du CFMR - Thèses en Mécanique des Roches December, 3 nd 2009 1 Overview Introduction Characterization
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 informationStress and fabric in granular material
THEORETICAL & APPLIED MECHANICS LETTERS 3, 22 (23) Stress and fabric in granular material Ching S. Chang,, a) and Yang Liu 2 ) Department of Civil Engineering, University of Massachusetts Amherst, Massachusetts
More informationDiscrete element modeling of self-healing processes in damaged particulate materials
Discrete element modeling of self-healing processes in damaged particulate materials S. Luding 1, A.S.J. Suiker 2, and I. Kadashevich 1 1) Particle Technology, Nanostructured Materials, DelftChemTech,
More informationPowder Technology 205 (2011) Contents lists available at ScienceDirect. Powder Technology. journal homepage:
Powder Technology 25 (211) 15 29 Contents lists available at ScienceDirect Powder Technology journal homepage: www.elsevier.com/locate/powtec Numerical simulation of particle breakage of angular particles
More informationModeling of stress-dependent static and dynamic moduli of weak sandstones
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2011jb009019, 2012 Modeling of stress-dependent static and dynamic moduli of weak sandstones Liming Li 1 and Erling Fjær 1,2 Received 11 November
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
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 informationModule-4. Mechanical Properties of Metals
Module-4 Mechanical Properties of Metals Contents ) Elastic deformation and Plastic deformation ) Interpretation of tensile stress-strain curves 3) Yielding under multi-axial stress, Yield criteria, Macroscopic
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 Rock Engineering
Rock Mechanics and Rock Engineering Overview Rock mechanics is the theoretical and applied science of the mechanical behaviour of rock and rock masses. Rock mechanics deals with the mechanical properties
More informationCh 4a Stress, Strain and Shearing
Ch. 4a - Stress, Strain, Shearing Page 1 Ch 4a Stress, Strain and Shearing Reading Assignment Ch. 4a Lecture Notes Sections 4.1-4.3 (Salgado) Other Materials Handout 4 Homework Assignment 3 Problems 4-13,
More informationShear Rupture of Massive Brittle Rock under Constant Normal Stress and Stiffness Boundary Conditions
Thesis Summary October 2015 Shear Rupture of Massive Brittle Rock under Constant Normal Stress and Stiffness Boundary Conditions Abstract Robert Paul Bewick, Ph.D., P.Eng. Supervisors: Peter K Kaiser (recently
More information1.8 Unconfined Compression Test
1-49 1.8 Unconfined Compression Test - It gives a quick and simple measurement of the undrained strength of cohesive, undisturbed soil specimens. 1) Testing method i) Trimming a sample. Length-diameter
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 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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationConstitutive modelling of fabric anisotropy in sand
Geomechanics from Micro to Macro Soga et al. (Eds) 2015 Taylor & Francis Group, London, ISBN 978-1-138-02707-7 Constitutive modelling of fabric anisotropy in sand Z.W. Gao School of Engineering, University
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 informationDiscrete element investigation of rate effects on the asymptotic behaviour of granular materials
Discrete element investigation of rate effects on the asymptotic behaviour of granular materials David MAŠÍN a,1 and Jan JERMAN a a Faculty of Science, Charles University in Prague, Czech Republic Abstract.
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 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 informationHIGHWALL STABILITY DUE TO GROUND VIBRATIONS FROM BLASTING. Dr. Kyle A. Perry Dr. Kot F. Unrug Kevin Harris Michael Raffaldi
HIGHWALL STABILITY DUE TO GROUND VIBRATIONS FROM BLASTING Dr. Kyle A. Perry Dr. Kot F. Unrug Kevin Harris Michael Raffaldi PURPOSE Objective: Study the effect of blast vibrations on the stability of highwalls
More informationGEO E1050 Finite Element Method Mohr-Coulomb and other constitutive models. Wojciech Sołowski
GEO E050 Finite Element Method Mohr-Coulomb and other constitutive models Wojciech Sołowski To learn today. Reminder elasticity 2. Elastic perfectly plastic theory: concept 3. Specific elastic-perfectly
More information3D Finite Element analysis of stud anchors with large head and embedment depth
3D Finite Element analysis of stud anchors with large head and embedment depth G. Periškić, J. Ožbolt & R. Eligehausen Institute for Construction Materials, University of Stuttgart, Stuttgart, Germany
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
51 Module 4: Lecture 2 on Stress-strain relationship and Shear strength of soils Contents Stress state, Mohr s circle analysis and Pole, Principal stressspace, Stress pathsin p-q space; Mohr-coulomb failure
More informationPROBABILISTIC CALIBRATION OF A DISCRETE PARTICLE MODEL. A Thesis YANBEI ZHANG
PROBABILISTIC CALIBRATION OF A DISCRETE PARTICLE MODEL A Thesis by YANBEI ZHANG Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree
More informationINCREASING RUPTURE PREDICTABILITY FOR ALUMINUM
1 INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM Influence of anisotropy Daniel Riemensperger, Adam Opel AG Paul Du Bois, PDB 2 www.opel.com CONTENT Introduction/motivation Isotropic & anisotropic material
More informationSHEAR STRENGTH OF SOIL UNCONFINED COMPRESSION TEST
SHEAR STRENGTH OF SOIL DEFINITION The shear strength of the soil mass is the internal resistance per unit area that the soil mass can offer to resist failure and sliding along any plane inside it. INTRODUCTION
More informationThree-Dimensional Discrete Element Simulations of Direct Shear Tests
Three-Dimensional Discrete Element Simulations of Direct Shear Tests Catherine O Sullivan Department of Civil and Environmental Engineering, Imperial College London, UK Liang Cui Department of Civil Engineering,
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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationCoupling DEM simulations and Physical Tests to Study the Load - Unload Response of an Ideal Granular Material
Coupling DEM simulations and Physical Tests to Study the Load - Unload Response of an Ideal Granular Material Liang Cui School of Electrical, Electronic and Mechanical Engineering, University College Dublin
More informationStructural behaviour of traditional mortise-and-tenon timber joints
Structural behaviour of traditional mortise-and-tenon timber joints Artur O. Feio 1, Paulo B. Lourenço 2 and José S. Machado 3 1 CCR Construtora S.A., Portugal University Lusíada, Portugal 2 University
More informationAnalysis in Geotechnical Engineering
EOSC433: Geotechnical Engineering Practice & Design Lecture 11: Discontinuum Analysis & the Distinct-Element Method 1 of 45 Erik Eberhardt UBC Geological Engineering EOSC 433 (2017) Analysis in Geotechnical
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 informationTheory at a Glance (for IES, GATE, PSU)
1. Stress and Strain Theory at a Glance (for IES, GATE, PSU) 1.1 Stress () When a material is subjected to an external force, a resisting force is set up within the component. The internal resistance force
More informationCalculation of periodic roof weighting interval in longwall mining using finite element method
Calculation of periodic roof weighting interval in longwall mining using finite element method Navid Hosseini 1, Kamran Goshtasbi 2, Behdeen Oraee-Mirzamani 3 Abstract The state of periodic loading and
More informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More informationChapter 7. Highlights:
Chapter 7 Highlights: 1. Understand the basic concepts of engineering stress and strain, yield strength, tensile strength, Young's(elastic) modulus, ductility, toughness, resilience, true stress and true
More informationDiscrete Element Modelling in Granular Mechanics WUHAN seminars, China, 28th august-2nd september 2017
Discrete Element Modelling in Granular Mechanics WUHAN seminars, China, 28th august-2nd september 2017 F. DARVE, L. SIBILLE, F. DUFOUR Laboratoire Sols, Solides, Structures, Risques, INPG-UGA-CNRS, Grenoble,
More informationCracked concrete structures under cyclic load
Cracked concrete structures under cyclic load Fabrizio Barpi & Silvio Valente Department of Structural and Geotechnical Engineering, Politecnico di Torino, Torino, Italy ABSTRACT: The safety of cracked
More informationPREDICTING STRENGTH OF PILLARS IN FRACTURED ROCK MASS DURING DEVELOPMENT OF APATITE-NEPHELINIC ORES
PREDICTING STRENGTH OF PILLARS IN FRACTURED ROCK MASS DURING DEVELOPMENT OF APATITE-NEPHELINIC ORES Vladimir Leonidovich Trushko, Anatoliy Grigoryevich Protosenya and Pavel Eduardovich Verbilo The Saint-Petersburg
More informationCity, University of London Institutional Repository
City Research Online City, University of London Institutional Repository Citation: Li, Y. Q., Hu, Z., Fang, X. & Fonseca, J. (2015). Analysis of micro characteristics and influence factors of foundation
More informationBifurcation Analysis in Geomechanics
Bifurcation Analysis in Geomechanics I. VARDOULAKIS Department of Engineering Science National Technical University of Athens Greece and J. SULEM Centre d'enseignement et de Recherche en Mecanique des
More informationISMS Paper No Influence of weak planes on rockburst occurrence. Amin Manouchehrian, Ming Cai *
Paper No. 176 ISMS 2016 Influence of weak planes on rockburst occurrence Amin Manouchehrian, Ming Cai * Bharti School of Engineering, Laurentian University, Sudbury, Canada, P3E 2C6 MIRARCO, Laurentian
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
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 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 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 informationComparison of six major intact rock failure criteria using a particle flow approach under true-triaxial stress condition
Geomech. Geophys. Geo-energ. Geo-resour. (1) :3 9 DOI.7/s9-1-3- ORIGINAL ARTICLE Comparison of six major intact rock failure criteria using a particle flow approach under true-triaxial stress condition
More informationEOSC433: Geotechnical Engineering Practice & Design
EOSC433: Geotechnical Engineering Practice & Design Lecture 1: Introduction 1 of 31 Dr. Erik Eberhardt EOSC 433 (Term 2, 2005/06) Overview This course will examine different principles, approaches, and
More informationRock Mechanics Laboratory Tests for Petroleum Applications. Rob Marsden Reservoir Geomechanics Advisor Gatwick
Rock Mechanics Laboratory Tests for Petroleum Applications Rob Marsden Reservoir Geomechanics Advisor Gatwick Summary A wide range of well established and proven laboratory tests are available for petroleum
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 informationModelling Progressive Failure with MPM
Modelling Progressive Failure with MPM A. Yerro, E. Alonso & N. Pinyol Department of Geotechnical Engineering and Geosciences, UPC, Barcelona, Spain ABSTRACT: In this work, the progressive failure phenomenon
More informationThe Stress Variations of Granular Samples in Direct Shear Tests using Discrete Element Method
The Stress Variations of Granular Samples in Direct Shear Tests using Discrete Element Method Hoang Khanh Le 1), *Wen-Chao Huang 2), Yi-De Zeng 3), Jheng-Yu Hsieh 4) and Kun-Che Li 5) 1), 2), 3), 4), 5)
More informationAim of the study Experimental determination of mechanical parameters Local buckling (wrinkling) Failure maps Optimization of sandwich panels
METNET Workshop October 11-12, 2009, Poznań, Poland Experimental and numerical analysis of sandwich metal panels Zbigniew Pozorski, Monika Chuda-Kowalska, Robert Studziński, Andrzej Garstecki Poznan University
More informationForce Transmission Modes of Non-Cohesive and Cohesive Materials at the Critical State
materials Article Force Transmission Modes of Non-Cohesive and Cohesive Materials at the Critical State Ji-Peng Wang ID Building Architecture and Town Planning Department (BATir), Université Libre de Bruxelles,
More informationNumerical Assessment of the Influence of End Conditions on. Constitutive Behavior of Geomaterials
Numerical Assessment of the Influence of End Conditions on Constitutive Behavior of Geomaterials Boris Jeremić 1 and Zhaohui Yang 2 and Stein Sture 3 ABSTRACT In this paper we investigate the behavior
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 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 information