Advanced Search

Indexed by SCI、CA、РЖ、PA、CSA、ZR、etc .

Baochun ZHOU, Jingtao WANG, Jun WEI. Method of Numerical Modeling for Constitutive Relations of Clay. Journal of Earth Science, 2006, 17(4): 355-360.
Citation: Baochun ZHOU, Jingtao WANG, Jun WEI. Method of Numerical Modeling for Constitutive Relations of Clay. Journal of Earth Science, 2006, 17(4): 355-360.

Method of Numerical Modeling for Constitutive Relations of Clay

Funds:

the Natural Science Foundation of Henan Province 0511045200

the Youth Science Foundation of Xinyang Normal University 20050107

More Information
  • In order to study the method of numerical modeling for constitutive relations of clay, on the basis of the principle of interaction between plastic volumetric strain and plastic generalized shear strain, the two constitutive functionals that include the function of stress path were used as the basic framework of the constitutive model, which are able to demonstrate the dependence of stress path. The two partial differential cross terms appear in the expression of stress-strain increment relation, which are used to demonstrate the interaction between plastic volumetric strain and plastic generalized shear strain. The elasoplastic constitutive models of clay under two kinds of stress paths, CTC and TC, have been constracted using the triaxtal test results. The three basic characteristics of deformation of soils, pressure sensitivity, dilatancy, and dependence of stress path, are well explained using these two models. Using visualization, the three-dimensional surfaces of shear and volume strains in the whole stress field under stress paths of CTC and TC are given. In addition, the two families of shear and volmetric yield Ioei under CTC and TC paths are plotted respectively. By comparing the results of deformation under these two stress paths, it has been found that, there are obvious differences in the strain peaks, the shapes of strain surfaces, and the trends of variation of volumetric yield loci, however both families of shear yield loci are similar. These results demonstrate that the influences of stress path on the constitutive relations of clay are considerably large and not negligible. The nmericul modeling method that can sufficiently reflect the dependence of stress path is superior to the traditional one.

     

  • In jointed and fractured rock masses, it is often difficult to achieve a smooth fracture plane between blasted holes.McKown (1984) indicated that when the angle between rock structural planes and the blast-induced plane is less than 60°, it is difficult to achieve a smooth-blasting surface.From the viewpoint of rock strength theory, Zhang (1988) pointed out that, when the angle between rock structural plane and blast-induced plane is within 25° to 40°, the blast-fractured surface is jagged in the form of "Z" pattern.

    Generally, the factors, which affect blasting results, may be divided into those that can be controlled and those that cannot be controlled.The controllable factors include explosive properties, initiation timing and blast geometry.The uncontrollable factors comprise the rock natural structures, such as joints and fractures, and properties, such as elastic constants, density and strength. Among these, the influence of rock structural planes often contributes to a high degree of variability of blasting results.

    This paper presents a theoretical analysis of rock structural plane influence on blasting results based on elasticity and stress wave propagation theory with an emphasis on smooth blasting techniques.Simulated experimental results, using strain and acoustic emission measurements for verifying the theoretical analysis, are also given.

    In blasting operations, the rock structure conditions shown in Fig. 1 are often encountered, the rock structural plane and blastinduced fracture plane intersect at angle β, where β is the rock structural plane direction with respect to the plane intersecting the boreholes.A line normal to the rock structural plane is given by 'n', and α is the incident angle of the blast-induced stress wave. In other words, α is the angle between the incident wave (blastinduced stress wave) and 'n' (the normal line of structural plane).By geometry, it can be seen that α+β=90°.

    Figure  1.  Relationship between rock structural plane and blast-induced fracture planes.

    To form a smooth fractured plane without wall damage, the tangential explosive stress component, σθ, must be greater than the rock mass dynamic tensile strength of rock and the shear stress component, τα, induced on the rock structural plane must be less than rock joint shear strength as determined by Coulomb's criterion

    (1)

    where σtd is the dynamic tensile strength of the rock joint, while c is the cohesive force on rock structural plane, and ϕ is the friction angle of joint.

    According to elasticity theory (Yin, 1990), normal and shear stress components on the rock structural plane in Fig. 1are determined by the following equations

    (2)

    in which the relationship between the tangential stress component, σθ, and radial stress component, σr, is given by

    (3)

    and ηCs/Cp.Cs and Cp are the propagation speed of shear wave and pressure wave in rock masses, respectively.Suppose that K=1-2η2, then

    (4)

    According to the wave theory (Guo, 1982), the following equations can be established.

    (5)

    In the above equations, E is Young's modulus; ρ is the density, and μ is the Poisson ratio of the rock.Using the relations in equation (5), K can be defined as

    (6)

    By substituting equations (2) into the expression for shear strength (equation(1)) and using equations (4) and (6), replacing α with 90°-β

    (7)

    From equation(7), it can be seen that the blast-induced stress components in the rock mass are functions of Poisson's ratio, μ, cohesion, c, friction angle, φ, and the orientation of the structural plane, β.

    (8)

    Previous research by Wu (1997) shows that the rock structural plane orientation has an influence on the stress wave propagation only when incident angle, α, is within a certain range. Shown in Fig. 1, when the angle β ranges from (90°-α1) to(90°- α2), the rock structural plane affects the propagation of blast-induced stress wave.

    For incident pressure (P) wave with amplitude A1, shown in Fig. 2.

    Figure  2.  Incident pressure wave on rock structure plane.

    The transmission coefficient (A5) and the reflection coefficient (A2) of the P-wave can be calculated by the following (Kolsky, 1983)

    (9)

    Given Poisson's ratio μ=021 and friction angle ϕ=10°, 20°, 25°, and 30° respectively, the computed reflection and transmission coefficients are shown in Fig. 3.From Fig. 3, it can be shown that, regardless of the friction angle value, the transmission coefficient (A5) is minimum while the reflection coefficient (A2) is maximum when the incident angle α is about 60°(β≈30°).T his indicates that the transmission wave energy is a minimum while the reflection wave energy is a maximum when β is about 30°.T he utilization ratio of explosive energy is the lowest under such blasting conditions.In other words, such blasting conditions cannot achieve a smooth fracture surface.

    Figure  3.  Reflection and transmission coefficient with different friction angle ϕ.A2.reflection coefficient; As.transmission coefficient.

    A prerequisite for equation (9) is that the shear stress component on the rock structural plane induced by the incident pressure wave must be larger than or at least equal to the frictional resistance, which can be written as the following form

    (10)

    In Fig. 2, η2=sin2β/sin2α, so equation (10) can be defined as another form.

    (11)

    Plotting values for α and b using equation (11) when Poisson's ratio, μ, takes on the value of 0.25, it can be shown in Fig. 4 that the rock structural plane affects the propagation of blast-induced wave as the incident angle α ranges from 30°to 80° (the friction angle b is about 18°and β is calculated to range from 10°to 60°under this circumstance).

    Figure  4.  Ranges of incident angle α determined by friction angle ϕ.

    Given the theoretical analysis above, it can be stated when the orientation of the rock structural plane and the blast-induced fracture plane are identical, the stresses on the rock structural plane induced by blasting should satisfy the stress-state control equation (7).Rock blasting under this condition will achieve a good smooth surface along the central line of boreholes.If the rock structure angle, β, does not satisfy equation (7) and the joints exhibit a lowly shear strength, the blast-induced surface is usually jagged or in the form of a'Z'.In this case, the rock will be fractured on the structural plane by shear stress, and the tensile stress applied on the structural plane has no effect on the rock breakage.As such, smooth wall fracturing cannot be obtained.

    Blasting experiments in marble and cement mortar were carried out in laboratory conditions using simulated rock structural planes.Marble cores, 50 mm in diameter and 300 mm in length, were prepared with a hole 10 mm in diameter and 60 mm in depth.Shown in Fig. 5, the marble core was cut using a diamond saw to an angle β and glued with a binder to simulate the rock structural plane.An electronic detonator, loaded in the borehole, was used as the explosive source.FoiL-l type strain gages, 3 mm× 15 mm in size, were mounted on both sides of the joint.T he data acquisition system is shown in Fig. 5 and consists of a KD-54super-dynamic amplifier and a CS2092waveform recorder/analyzer.

    Figure  5.  Marble core in strain gage experiment and measurement system.

    Table 1 gives the peak strain results of these experiments as a function of β.Figure 6 shows the relationship between ε2/ε1 and β, showing that the strain ratio, ε2/ε1, is a minimum when β≈30°.T his indicates that a joint orientation of 30°promotes the highest attenuation of explosive energy.In such a case, it will be difficult to achieve a smooth blast-fractured surface.

    Table  1.  RESULTS OF STRAIN MEASUREMENT
     | Show Table
    DownLoad: CSV
    Figure  6.  Relationship between ε2/ε1 and β.

    A second series of experiments were conducted using a cement mortar (cement∶sand∶water=1∶3∶1), to simulate the rock mass, cutting joints similar to the previous experiment.Acoustic emission measurements were made for varying values of β using the set-up shown in Fig. 7.

    Figure  7.  Cement mortar in acoustic emission experiment.

    The results of these tests are given in Table 2.

    Table  2.  EXPERIMENT RESULTS OF ACOUSTIC EMISSION
     | Show Table
    DownLoad: CSV

    Figure 8 shows the relationship between energy per unit sample length, E', and β.It can be seen from Fig. 8 that E' has minimum value when β≈30°, providing similar results to the previous experiments.Typical acoustic waveforms are shown in Fig. 9.

    Figure  8.  Relationship between E' and β.
    Figure  9.  Typical acoustic waveform.

    Both theoretical analysis and simulated experiments in the laboratory show that the influence of rock structural planes contributes to a high degree of variability of blasting results.T he influence of rock structural planes on blasting results depends on the angles between the rock structural planes and the blast-induced fracture planes.It is difficult to achieve better blasting results when the angles between the rock structural planes and the blastinduced fracture planes range from 10° to 60°.Among these angles, 30° is the least favorable angle for achieving a smooth wall. For the angles less than10°and greater than60°, the influence of rock structural planes on blasting results can be ignored.

    ACKNOWLEDGMENT: The authors would like to give special thanks to Dr.Catherine Aimone-Martin, a professor at New Mexico Institute of Mining and Technology, for her help in English writing.
  • Huang, W. X., 1983. Engineering Characteristics of Soil. Hydraulic and Electric Press, Beijing (in Chinese)
    Lade, P. V., 1988. Effects of Voids and Volume Changes on the Behavior of Frictional Materials. International Journal for Numerical and Analy tical Methods in Geomechanics, (12): 351-370
    Lu, P. Y., Chen, S. Y., Xiong, L. Z., et al., 1984. Elastoplastic Constitutive Equation for Hydraulic-Fill Soil of Slurry-Fall Dam. Chinese Journal of Geotechnical Engineering, 6(2): 23-39 (in Chinese with English Abstract)
    Nanjing Hydraulic Research Institute, 1999. Specification of Soil Test. Water Resources and Electric Power Press, Beijing (in Chinese)
    Ren, Q. Y., Wang, J. T., 2005. Numerical Method for Modeling the Constitutive Relationship of Sand under Different Stress Paths. Journal of China University of Geosciences, 16(3): 268-270
    Tarantola, A., 1987. Inverse Problem Theory. Elsevier Science Publishers, Amsterdam.
    Wang, J. T., 2002. Numerical Method in Modeling the Constitutive Relations of Rock and Soil. Journal of Huazhong University of Science and Technology (Urban Science Edition), 19(1): 44-47 (in Chinese with English Abstract)
    Wang, J. T., 2004. Particularity and Unity of the Constitutive Relations for Rock and Soil. Journal of Huazhong University of Science and Technology (Urban Science Edition), 21(4): 5-8 (in Chinese with English Abstract)
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(4)  / Tables(1)

    Article Metrics

    Article views(602) PDF downloads(11) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return