The Scientific World Journal
Volume 2020 (2020), Article ID 2978257, 20 pages
https://doi.org/10.1155/2020/2978257

Research Article

Modeling and Numerical Analysis in 3D of Anisotropic and Nonlinear Mechanical Behavior of Tournemire Argillite under High Temperatures and Dynamic Loading

Foguieng Wembe Marius,¹ Mambou Ngueyep Luc Leroy,2,3 and Ngapgue François1,4

¹Unité de Recherche de Mécanique et de Modélisation des Systèmes Physiques (UR-2MSP), Department of Physics, Dschang School of Science and Technology, University of Dschang, P.O. Box 67, Dschang, Cameroon
²Laboratory of Material Sciences, Department of Physics, Faculty of Science, University of Yaoundé 1, P.O. Box 812, Yaoundé, Cameroon
³Department of Mine Mineral Processing and Environment, School of Geology and Mining Engineering, University of Ngaoundéré, P.O. Box 115, Meiganga, Cameroon
⁴Laboratory of Industrial and Systems Engineering Environment (LISIE), Department of Civil Engineering, Fotso Victor Institute of Technology, Dschang School of Science and Technology, University of Dschang, P.O. Box 134, Bandjoun, Cameroon

Correspondence should be addressed to Mambou Ngueyep Luc Leroy; mamboulucleroy@gmail.com

Received 4 February 2020; Revised 30 May 2020; Accepted 1 June 2020; Published 23 June 2020

Academic Editor: Silvia Scaglione

Copyright © 2020 Foguieng Wembe Marius et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This work proposes a model that takes into account the anisotropy of material with its inhomogeneity and geometrical and material nonlinearities. According to Newton’s second law, the investigations were carried out on the simultaneous effects of mechanical load and thermal treatment on the Tournemire argillite material. The finite difference method was used for the numerical resolution of the problem by the MATLAB 2015a software in order to determine the peak stress and strain of argillite as a function of material nonlinearity and demonstrated the inhomogeneity parameter Ω. The critical temperature from which the material damage was pronounced is 500°C. Indeed, above this temperature, the loss of rigidity of argillite reduced significantly the mechanical performance of this rock. Therefore, after 2.9 min, the stress reduction in X or Y direction was 75.5% with a peak stress value of 2500 MPa, whereas in Z direction, the stress reduction was 74.1% with a peak stress value of 1998 MPa. Meanwhile, knowing that the material inhomogeneity was between 2995 and 3256.010, there was an increase in peak stress of about 75%. However, the influence of the material nonlinearity was almost negligible. Thus, the geometrical nonlinearity allows having the maximal constant strain of about 1.25 in the direction of the applied dynamic mechanical force.

1. Introduction

High temperature has a great influence on the microstructure of rocks as reported in several literature investigations [1–3]. However, under high temperature, the microstructure of some rocks can completely change, while it may be rearranged at certain temperature. Thermal processing of rocks involves many microcracks and microvoids and intensifies granular cleavage planes which have significant effect on physical properties of rocks such as porosity, density, and permeability [4–6].

Consequently at certain levels of temperature, the mechanical properties of rocks such as elastic modulus, Poisson’s ratio, thermal expansion coefficient, tensile, and compressive strength decrease [7–18].

From the experimental point of view, previous elucidated works and others [19–21] have been done in this literature. However, few theoretical studies had been done and numerical results in this field were obtained. For instance, Nubissie et al. [22] investigated the dynamic behavior of a wooden beam under mechanical loading and fire. The authors established one model that takes into account the variation of physicomechanical parameters of wooden beam as a function of temperature. At the end, the authors concluded that the prediction of time to the structural failure is necessary for safety consideration. Ndoukouo et al. [23] studied the dynamics of fire-exposed steel beam under mechanical load and showed that an increase in the deflection versus time grows, while the bending moment presents a nonmonotonic behavior under a sinusoidal load. In the case of rocks, Mambou et al. [24] investigated numerically the mechanical properties of a granite rock specimen subjected to uniaxial loading and fire to analyze the internal stress and strain. From this investigation, they reported that beyond three minutes of exposure to thermal load, the mechanical energy required to fragment its rock specimen reduced up to 80%. Once more, in our previous work [25], we investigated the theoretical behavior of mechanical properties of sandstone rock specimen at high temperatures. At low stress, due to the closure of microcracks and changes in mechanical properties of the rock, we have introduced the material nonlinearity in the established model and showed the loss of rigidity of this sandstone. The same work also reported that 450°C is the critical temperature required to damage the physical and mechanical properties. Recently, Wang et al. [26] have studied the effects of treatment temperature and strain rate on the mechanical behaviors of granite samples. They used a statistical damage constitutive model for the rock based on the Weibull distribution to characterize the entire stress-strain response during rock failure. Finally, they showed that at high temperature, the enhancement effect of the strain rate on dynamic compressive strength is permanent. All these cited works did not deal with the anisotropy and geometric nonlinear behavior of rocks. From our best knowledge, the model of mechanical response of rock under high temperature taking into account the effect of the anisotropy and geometric nonlinear behavior of argillite rocks is not reported in the literature. One of the fundamental behaviors of rock specimen which is not explored in high temperature is known as anisotropy.

In this regard, the main objective of this work was to establish a model that takes into account the anisotropy of argillite with inhomogeneity and geometrical and material nonlinearities. The other objectives of this study include determining the mechanical behavior of rock specimen in 3D when subjected to high temperature and mechanical load, evaluating the peak stress and strain as a function of temperature and time with the corresponding inhomogeneity.

2. Modeling of Argillite Rock Specimen under Uniaxial Mechanical and Thermal Load

Figure 1 depicts a rock specimen subjected to thermomechanical load according to the ISRM norm.

Figure 1: Rock specimen under thermal load (Fth) and dynamic mechanical load (Fm).

By applying Newton’s second law to the model presented in Figure 1, we obtained the following equation (1) with and V the bulk density and volume, respectively:

where is the sum of force due to inertia, is the sum of internal forces due to internal stresses, is the sum of forces due to gravity, and represents the vector of displacement of rocks particle.

Following X, Y, and Z directions, equation (1) becomes

where represent the spatial coordinates of rock; σxy = σyx, σxz = σzx, and σyz = σzy, σxx, σyy, and σzz are the components of stress tensor in the X, Y, and Z directions; and intensity of gravity.

In laboratory, the anisotropy is usually investigated by the standard testing practices such as uniaxial compressive strength, triaxial test, and direct shear strength. In order to evaluate the anisotropic behavior of this rock, a uniaxial compression test should be carried out on the specimen on X, Y, and Z directions and then the elastic parameters in different directions were determined. In this case, tangential stresses should be equal to zero () and equation (2) becomes

The stress-strain relation given by Hooke’s law applied on rock mechanics is as follows:

To describe the elastic response of a transverse isotropic material, five independent elastic constants are necessary E₁, E₂, υ₁, υ₂ and G₁₂. If the isotropy plane is the XY plane, the parameters (E₁, υ₁₂) are determined from the uniaxial compression tests carried out in the plane. On the other hand, the parameters (E₂, υ₂₃) are determined by the tests carried out in the direction perpendicular to the plane.

From this assumption and as reported by Masri et al. [27], argillite presents transverse isotropic behavior, with plane XY as the symmetric plane, then Poisson’s coefficients, and Young’s modulus in the 3 directions, respectively. Equation (4) becomes

2.1. Mechanical Effect

Due to their polymineral constitution, most of the rocks contain voids and microcracks generally occupied by gases, water, and inclusions. If we assume that argillite presents a material nonlinearity behavior at the lower state of stresses characterized by the closing of the microcracks and the variation of the mechanical properties of rocks, then we can modify stress-strain relation as reported by Inserra et al. [28], considering a second order approximation in the form of

If we assume that a rock material has an inhomogeneity as found in the functionally graded materials, we can express the Young E modulus according to a power law [29] given by the following relations:

where represents the inhomogeneity parameter, and , the nominal Young’s modulus.

Then, similarly, we can express the coefficient of thermal expansion, the density, and Poisson coefficient as in the following equation:

where is the nominal thermal expansion, is the nominal density, and is the nominal Poisson’s coefficient.

In addition, the general equation of strain-displacement is given by equation (9) in X, Y, or Z directions:

The terms (∂u/∂x)², (/∂y)², and (/∂z)² represent geometric nonlinearity at the high state of stresses.

Many experimental works show that high temperature and mechanical loading affect physical and mechanical properties of rock. Considering the rock matrix with its geological history, the behavior of anisotropy may not be the same at high temperature. It is recognized as the dissimilarity response under mechanical or physical effects [30]. In light to this, equation (3) becomes

2.2. Thermal Effect (Fth)

In this section, we suppose that thermal excitation (Fth) is due to fire. For modeling of the fire effect, we use the mathematical formula of the ISO 834 fire as in ref. [31]: because it is a conventional fire which is used to have resistance tests for the numerical modeling of structures exposed to fire. The international standard time-temperature curve of the ISO 834 fire is defined as in ref. [31]:

where (min) is the time, and θ₀ = 20°C represents the room temperature; the thermal stress is calculated as in

where , the thermal strain, and , the variation of thermal expansion coefficient.

Finally, thermal force can be expressed as in

where d represents the length of specimen in X, Y, and Z directions.

By combining equations (6)–(11) and (13), we obtain the following equations governing the displacements of the rock:

By referring to the experimental works of Masri et al. [27], the coefficient of thermal expansion, density, Young’s modulus, and Poisson’s ratio were obtained by the following equations:

3. Numerical Analysis of Argillite Rock Specimen under Uniaxial Mechanical Load and Thermal Load on X or Y and Z Directions

To solve equations (14a)–(14c), initial conditions and boundary conditions could be well defined. We assume that the both ends of specimen are free and obtained equation (16) using uniaxial load in each direction of the specimen:

Initial condition is obtained by solving equation (16) without external forces and nonlinearities, and thus equation (16) becomes

Thus, the general solution assuming both ends of specimen are free is as follows:

where ɷnx = ɷny = nπ/d.(E₀₁/ρ)1/2 and ɷnz = nπ/d.(E₀₂/ρ)1/2, with n = 1, 2, 3, …

Initial condition of our model is

We use the centered discretization scheme for numerical approach respectively of the first and second spatial derivatives and the second temporal derivative as follows:

Physical and mechanical parameters of argillite used in this subsection are adopted from the experimental works of Masri et al. [27] and the numerical work of Mambou et al. [24]:

4. Numerical Analysis of Argillite under Uniaxial Mechanical Load and Thermal Load in X or Y and Z Directions

The analysis in this part is done at the center of rock specimen (x = y = z = 0.025). For each uniaxial compression applied on X, Y, and Z directions, we plot the evolution of the internal stress and internal strain as a function of temperature and time. These temperatures vary from 20°C to 1120°C. We determine numerically the peak stress and peak strain for each direction, and then, the inhomogeneous parameter () which characterized the rock specimen for each temperature. The effect of nonlinearity parameter beta () was studied by taking values 0.001, 0.01, and 0.1 as in ref. [25].

4.1. Case of Argillite in Which Failure Occurred at 100°C

Figures 2 and 3 show the evolution of internal stress and internal strain versus temperature and time of the rock specimen subjected to mechanical load and fire for inhomogeneity parameter  = 2027.

(a)

(b)

(c)

(d)

Figure 2: Internal stress and strain of the specimen subjected to dynamic mechanical loading in X-axis and fire (for  = 2027). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy and strain_zz with time.

(a)

(b)

(c)

(d)

Figure 3: Internal stress and strain of the specimen subjected to dynamic mechanical loading in Z-axis and fire (for  = 2192.5). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

Figures 2(a) and 2(b) show the evolution of internal stress and internal strain versus temperature and time, respectively, of the rock specimen subjected to mechanical load and fire for inhomogeneity parameter and material nonlinearity . From these figures, considering different values of = 0.1, 0.01, and 0.001, we have observed the same evolution in X direction. But we have noted that for  = 0.1, we have observed peak stress at θ = 100°C and the corresponding time was t = 0.087 min. These figures have presented in the X or Y direction. Consequently, we can conclude that this internal stress with  = 0.1 dominated on the other stress which have equals to 0.01 or 0.001.

From Figures 2(c) and 2(d), we have in the X or Y direction the maximum strain and the peak strain at the same temperature and time. Then, the peak strain .

Figures 3(a) and 3(b) show that in the Z direction at 100°C and 0.087 min, the peak stress for beta = 10⁻¹ and , which is the maximum amplitude value of internal stress. In the case of strain, we have noted that and . The reduction of stress in Z direction is about 5.8% and in X or Y direction 4% at 100°C. In Figures 3(c) and 3(d), the peak strain value is 0.42 approximately in X or Y direction.

4.2. Case of Rock Specimen in Which Failure Occurred at 300°C

The curves in Figures 4 and 5 show the evolution of internal stress and internal strain versus temperature and time of the rock specimen subjected to mechanical load and fire for inhomogeneity parameter  = 2654 and  = 2800.40, respectively.

(a)

(b)

(c)

(d)

Figure 4: Internal stress and strain of the specimen subjected to dynamic mechanical loading in X-axis and fire (for  = 2654). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

(a)

(b)

(c)

(d)

Figure 5: Internal stress and strain of specimen submitted to dynamic mechanical loading in Z-axis and fire (for  = 2800.40). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

Considering in Figures 4(a) and 4(b) that different values of beta equal to 0.1, 0.01, and 0.001, we have the same behavior in X direction in terms of stress or strain.

In Figures 4(c) and 4(d), we have in the X or Y direction the peak strain at 300°C and 0.69 min. Then, the peak strain in the Z direction. We have noted that for beta = 10⁻¹, we have peak stress at θ = 300°C and the corresponding time is t = 0.69 min in the X direction with maximum stress.

Figures 5(a) and 5(b) show that in the Z direction at 300°C and 0.69 min, the peak stress for beta = 10⁻¹. In the case of strain, we noted that and . At 300°C, the reduction of is about 57.8%, and in X or Y direction, we have 49.5%. This percentage implies the beginning of damage of rock material. Less than 1 min (0.69 min) of exposure to fire at 100°C, we can conclude as in ref. [25] that the peak stress of this rock is reduced of about 47.42% in X or Y direction and 55.55% in Z direction.

4.3. Case of Rock Specimen in Which Failure Occurred at 500°C

Figures 6 and 7 show the behavior of internal stress and internal strain versus temperature and time of the rock specimen subjected to mechanical load and fire for inhomogeneity parameter  = 2995 and  = 3122.5, respectively.

(a)

(b)

(c)

(d)

Figure 6: Internal stress and strain of the specimen subjected to dynamic mechanical loading in X-axis and fire (for  = 2995). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

(a)

(b)

(c)

(d)

Figure 7: Internal stress and strain of the specimen subjected to dynamic mechanical loading in Z-axis and fire (for  = 3122.5). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

Considering in Figures 6(a) and 6(b) that different values of beta equal to 10⁻¹, 10⁻², and 10⁻³, we have the same evolution in Z direction.

In Figures 6(c) and 6(d), we have in the X or Y direction the peak strain at 500°C and 2.9 min. We have noted that for beta = 10⁻¹, peak stress at θ = 500°C and corresponding time is t = 2.9 min in the X direction.

Figures 7(a) and 7(b) show that at 500°C with associated time 2.9 min, the peak stress for beta = 10⁻¹ in the Z direction. In the case of strain, we noted that and . At 500°C, the reduction of stress is about 75.5% in the X or Y direction, and 74.1% in the Z direction. At this temperature, damage of material is very pronounced. Consequently, this temperature will be considered as critical temperature. This result is similar to that in refs. [24, 25].

4.4. Case of Rock Specimen in Which Failure Occurred at 600°C

Figure 8 shows the behavior of internal stress and internal strain versus temperature and time of the rock specimen subjected to mechanical load and fire for inhomogeneity parameter  = 3132.

(a)

(b)

(c)

(d)

Figure 8: Internal stress and strain of the specimen subjected to dynamic mechanical loading in X-axis and fire (for  = 3132). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

Figures 8(a) and 8(b) show the evolution of internal stress and internal strain versus temperature and time, respectively, of the rock specimen subjected to mechanical load and fire for inhomogeneity parameter and material nonlinearity beta. Considering in the same figures that different values of beta equal to 10⁻¹, 10⁻², and 10⁻³, we have the same evolution. But we have noted that for beta = 10⁻¹, we have peak stress at θ = 600°C and the corresponding time t = 5.9 min as in ref. [23] in which the temperature t = 5.87 min. We have also noted that the decrease in the peak stress values is approximately 67.71% in X or Y direction and 84.02% in Z direction compared with those values at 100°C.

In Figures 8(c) and 8(d), in the X or Y direction, the maximum strain and peak strain are the same at the same temperature and time. Then, the peak strain .

Figures 9(a) and 9(b) show in Z direction at 600°C and 5.9 min, the peak stress for beta = 10⁻¹. In the case of strain, we have noted that and . The reduction of stress in Z direction is about 85% and in X direction 78.8% at 600°C. This result clearly shows the significant effect of fire on the Tournemire argillite. These results are similar to those obtained in ref. [24].

(a)

(b)

(c)

(d)

Figure 9: Internal stress and strain of the specimen subjected to dynamic mechanical loading in Z-axis and fire (for  = 3250.50). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

4.5. Case of Rock Material in Which Failure Occurred at 700°C

Figure 10 plots the behavior of internal stress and internal strain versus temperature and time of the argillite rock specimen subjected to mechanical load and fire for inhomogeneity parameter = 3256.010.

(a)

(b)

(c)

(d)

Figure 10: Internal stress and strain of the specimen subjected to dynamic mechanical loading in X-axis and fire (for  = 3256.010). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

In Figures 10(a) and 10(b), the peak stress is in X or Y direction. In X direction, in the case of strain, we noted that and in Figures 10(c) and 10(d).

Figures 11(a) and 11(b) show at θ = 700°C and t = 10.8 min, the peak stress for beta = 10⁻¹ in Z direction. In the case of strain, we have observed that and .

(a)

(b)

(c)

(d)

Figure 11: Internal stress and strain of the specimen subjected to dynamic mechanical loading in Z-axis and fire (for  = 3368.17). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

Finally, at θ = 700°C in Z direction, the reduction of stress is about 92.3%, whereas in X or Y direction, it is 67.43%, suggesting that the material is damaged.

4.6. Case of Rock Specimen in Which Failure Occurred at 900°C

Figures 12 and 13 show the behavior of internal stress and internal strain versus temperature and time of the rock specimen subjected to mechanical load and fire for inhomogeneity parameter  = 3477.050.

(a)

(b)

(c)

(d)

Figure 12: Internal stress and strain of the specimen subjected to dynamic mechanical loading in X-axis and fire (for  = 3477.050). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

(a)

(b)

(c)

(d)

Figure 13: Internal stress and strain of the specimen subjected to dynamic mechanical loading in Z-axis and fire (for  = 3573). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

In Figures 12(a) and 12(b), we have the same evolution of stress for different beta values. In Figures 12(c) and 12(d), in the X or Y direction, the peak strain at 900°C and 44.5 min. We noted that for beta = 10⁻¹, the peak stress at θ = 900°C and the corresponding time is t = 44.5 min in the X or Y direction. At this temperature, internal stress tends to be zero.

Figures 13(a) and 13(b) show that at θ = 900°C and t = 44.5 min, the peak stress for beta = 10⁻¹ in Z direction. In the case of strain, we noted that and . At 900°C, in Z direction, the reduction of stress is about 98%, and in X direction, it is 97%.

4.7. Case of Rock Material in Which Failure Occurred at 1120°C

In Figure 14, we have noted the same behavior of stress for different beta values. In Figures 14(c) and 14(d), in the X or Y direction, the peak strain at 1120°C and 167 min. We have noted that for beta = 10⁻¹, the peak stress at θ = 1120°C and the corresponding time is t = 167 min in the X direction.

(a)

(b)

(c)

(d)

Figure 14: Internal stress and strain of the specimen subjected to dynamic mechanical loading in X-axis and fire (for  = 3672). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

Figures 15(a) and 15(b) show that at 1120°C and 199.5 min, the peak stress for beta = 0.1 in Z direction. In the case of strain, we noted that and .

(a)

(b)

(c)

(d)

Figure 15: Internal stress and strain of the specimen subjected to dynamic mechanical loading in Z-axis and fire (for  = 3772.0). (a) Uniaxial stress with temperature; (b) uniaxial stress with time; (c) strain_xx, strain_yy, and strain_zz with temperature; and (d) strain_xx, strain_yy, and strain_zz with time.

In X direction, in the case of strain, we noted that and . When θ = 1120°C, in Z direction, the reduction of stress is about 99.25%, whereas in X or Y direction, it is 99%, suggesting that the material is more damaged in its totality.

In general, Figures 5(c)–15(c) (in temperature) and 5(d)–15(d) (in time) show the same strain in X or Y direction when dynamic mechanical load is applied in Z direction. We obtained the peak strain value of . In Z direction, the maximum strain is . We can observe that the strain is higher in the direction at which dynamic mechanical load is applied than other directions. The observation results are same when dynamic mechanical force is applied in X or Y direction. In this case, the maximum constant strain is in X direction. In Z direction, we obtained . The maximal strain is constant in the direction at which the force is applied, but the strain increases gradually in another direction.

Many experimental works have been done for various temperature ranges from room temperature (20°C) to 1200°C. This work proposed the peak stress and peak strain which are obtained with associated material inhomogeneity in X or Y direction (Table 1) and in Z direction (Table 2), at certain temperatures.

Table 1: Inhomogeneity and peak stress at different temperatures in X or Y direction.

Temperature (°C)
100
300
500
600
700
900
1120

Inhomogeneity
2027.0
2654.0
2995.0
2150.0
995.0
300.0
101.39

Peak stress on X or Y direction (MPa)
9700
5100
2500
3132
3256.010
3477.050
3672

Table 2: Inhomogeneity and peak stress at different temperatures in Z direction.

Temperature (°C)
100
300
500
600
700
900
1120

Inhomogeneity
2192.5
2800.4
3122.5
3250.5
3368.17
3573.0
3772.0

Peak stress on Z direction (MPa)
7200
3200
1998
1150
595
155
50

As in Figure 16(a), we have noted that the evolution of peak stress in X or Y and Z direction decreased similarly for inhomogeneity below 3000. When inhomogeneity is up to 3000, the peak stress is independent of direction X or Z. In general, Figures 16(a) and 16(b) show that the peak stress decreases with inhomogeneity and temperature, respectively. For inhomogeneity greater than 3500, the peak stress in general tends to be zero. We have also noted that the peak stress in Z direction is lower than the peak stress in X or Y direction. Consequently, we can conclude that material has a great rigidity in X or Y direction. Moreover, it is noted that for inhomogeneity  = 3125 peak stress has the same value approximatively. Then, the corresponding temperature was between 500°C and 573.5°C. In this interval, many experimental works presented a phase transition from quartz to [7, 8, 32] around 573°C.

(a)

(b)

Figure 16: Evolution of the peak stress: (a) versus inhomogeneity; (b) versus temperature.

This model can predict the fire resistance of argillite rock compared with those results in other element structures such as wood beam [22], steel beam [23], granite rock [24], and sandstone rock [25], which are most adequate for engineering project in order to perform the safety time and avoid the damage of structure

5. Conclusion

In this work, the mechanical behavior in terms of internal stress and internal strain of anisotropic Tournemire argillite under high temperature and dynamic loading was investigated. As results, the peak stress, peak strain, and inhomogeneity parameters were predicted at different temperatures. In general, the internal stress with temperature and time decreases and tends to zero, while at the same time and temperature, internal strain increases. Inhomogeneity of rocks like argillite has a great influence on its mechanical properties as peak stress and peak strain. The damage temperature and time at which material will be destroying strongly depend on these parameters. The geometrical nonlinearity allows having the maximal constant strain of about 1.25 in the direction of the applied mechanical force. We recorded 500°C as a critical temperature at which damage of material was pronounced after 2.9 min. The reduction of stress was 75.5% in X direction with a peak stress value of 2500 MPa and 74.1% in Z direction with a peak stress value 1998 MPa. Consequently, the energy of damage of argillite is reduced at 75% with inhomogeneity included in interval 2995 to 3256.010. However, the material nonlinearity has a negligible influence on the material.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

O. Alm, “The influence of micro crack density on the elastic and fracture mechanical properties of stripa granite,” Physics of the Earth and Planetary Interiors, vol. 40, pp. 61–179, 1985.

L. Y. Zhang, X. B. Mao, and W. T. Lu, “Experimental research on mechanical properties of sandstone at high temperature,” Journal of Mining, Safety Engineering, vol. 24, no. 3, pp. 293–297, 2007.

A. K. Verma, M. S. Saini, T. N. Singh, A. Dutt, and R. K. Bajpai, “Effect of excavation stages on stress and pore pressure changes for an underground nuclear repository,” Arabian Journal of Geosciences, vol. 6, no. 3, pp. 635–645, 2013.

N. N. Bajpai, B. Mahanta, P. G. Ranjith, and T. N. Singh, “Effects of thermal treatment on physic-morphological properties of Indian fine-grained sandstone,” Bulletin of Engineering Geology and the Environment, vol. 78, no. 2, pp. 883–897, 2017.

L. He, Q. Yin, and H. Jing, “Laboratory investigation of granite permeability after high-temperature exposure,” Processes, vol. 6, no. 4, pp. 36–14, 2018.

X. Meng, W. Liu, and T. Meng, “Experimental investigation of thermal cracking and permeability evolution of granite with varying initial damage under high temperature and triaxial compression,” Advances in Materials Science and Engineering, vol. 2018, Article ID 8759740, 9 pages, 2018.

F. E. Heuze, “High temperature mechanical, physical and thermal properties of granitic rocks-a review,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 20, no. 1, pp. 3–10, 1983.

F. Homand-Etienne, “Mechanical behavior of rocks as a function of temperature,” Earth Sciences, vol. 46, p. 261, 1986.

A. Sygała, M. Bukowska, and T. Janoszek, “High temperature versus geomechanical parameters of selected rocks—the present state of research,” Journal of Sustainable Mining, vol. 12, no. 4, pp. 45–51, 2014.

X. Mao, L. Zhang, R. Liu, and M. Dan, “Mechanical and thermal damage properties of sandstone at high temperatures,” Electronic Journal of Geotechnical Engineering, vol. 19, pp. 3137–3150, 2014.

R. Qiu-Hua, W. Zhi, X. Hai-Feng, and X. Qiang, “Experimental study of mechanical properties of sandstone at high temperatures,” Journal of Central South University of Technology, vol. 14, no. 1, pp. 478–483, 2007.

L. Zhang, X. Mao, R. Liu, X. Guo, and D. Ma, “The mechanical properties of mudstone at high temperatures: an experimental study,” Rock Mechanics and Rock Engineering, vol. 47, no. 4, pp. 1479–1484, 2014.

L. T. Yang, T. C. Ekneligoda, D. Wanatowski, A. M. Marshall, and L. R. Stace, “Experimental study on the mechanical properties of sandstone at high temperature,” in Proceedings of the 7th International Conference on Mining Science and Technology, Xuzhou, China, 2015.

P. K. Gautam, A. K. Verma, M. K. Jha, K. Sarkar, T. N. Singh, and R. K. Bajpai, “Study of strain rate and thermal damage of Dholpur sandstone at elevated temperature,” Rock Mechanics and Rock Engineering, vol. 49, no. 9, pp. 3805–3815, 2016.

S.-Q. Yang, W.-L. Tian, and G. P. Ranjith, “Failure mechanical behavior of Australian Strathbogie granite at high temperatures: insights from particle flow modeling,” Energies, vol. 10, no. 756, pp. 1–19, 2017.

H. Tian, G. Mei, G.-S. Jiang, and Y. Qin, “High-temperature influence on mechanical properties of diorite,” Rock Mechanics and Rock Engineering, vol. 50, no. 6, pp. 1661–1666, 2017.

T. D. Rathnaweera, P. G. Ranjith, X. Gu et al., “Experimental investigation of thermomechanical behaviour of clay-rich sandstone at extreme temperatures followed by cooling treatments,” International Journal of Rock Mechanics and Mining Sciences, vol. 107, pp. 208–223, 2018.

X. Wu, Z. Huang, H. Song et al., “Variations of physical and mechanical properties of heated granite after rapid cooling with liquid nitrogen,” Rock Mechanics and Rock Engineering, vol. 52, no. 7, pp. 2123–2139, 2019.

Z. Peng and Z.-J. Feng, “Thermal deformation of granite under different temperature and pressure pathways,” Advances in Materials Science and Engineering, vol. 2019, Article ID 7869804, 8 pages, 2019.

D. Healy, E. N. Timms, and A. M. Pearce, “The variation and visualisation of elastic anisotropy in rock forming minerals,” Solids Earth, vol. 11, no. 2, pp. 259–286, 2019.

Y. N. L. Wong, Y. Zhang, and Z. Wu, “Rock strengthening or weakening upon heating in the mild temperature range?” Engineering Geology, vol. 272, Article ID 105619, pp. 1–32, 2020.

A. Nubissie, N. A. Ndoukouo, and P. Woafo, “Dynamical behaviour of a wooden beam under mechanical loading and fire,” Materials and Design, vol. 32, pp. 1331–1336, 2010.

A. N. Ndoukouo, A. Nubissie, and P. Woafo, “On the dynamics of fire-exposed steel beam under mechanical load,” Journal of Constructional Steel Research, vol. 67, no. 12, pp. 1864–1871, 2011.

N. L. L. Mambou, J. Ndop, and J.-M. Ndjaka, “Modeling and numerical analysis of granite rock specimen under mechanical loading and fire,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 7, pp. 105–108, 2014.

N. L. L. Mambou, J. Ndop, and J. M. B. Ndjaka, “Theoretical investigations of mechanical properties of sandstone rock specimen at high temperatures,” Journal of Mining Science, vol. 50, no. 1, pp. 69–80, 2015.

Z. L. Wang, H. Shi, and J. G. Wang, “Mechanical behavior and damage constitutive model of granite under coupling of temperature and dynamic loading,” Rock Mechanics and Rock Engineering, vol. 51, no. 10, pp. 3045–3059, 2018.

M. Masri, M. Sibai, J. F. Shao, and M. Mainguy, “Experimental investigation of the effect of temperature on the mechanical behavior of Tournemire shale,” International Journal of Rock Mechanics and Mining Sciences, vol. 70, pp. 185–191, 2014.

C. Inserra, S. Biway, and C. Youqing, “Non linear ultrasonic characterization of thermally damaged westerly granite,” Japanese Journal of Applied Physics, vol. 48, pp. 101–106, 2009.

A. N. Eraslan and T. Akis, “On the plane strain and plane stress solutions of functionally graded rotating solid shaft and solid disk problems,” Acta Mechanica, vol. 181, no. 1–2, pp. 43–63, 2006.

H. A. Abbass, Z. Mohamed, and S. F. Yasir, “A review of methods, techniques and approaches on investigation of rock anisotropy,” AIP Conference Proceedings, vol. 1, pp. 2–9, 2018.

J. Seputro, “Effect of support conditions on steel beams exposed to fire,” School of Engineering, University of Canterbury, Christchurch, New Zealand, 2001, Master thesis in Fire Engineering.

E. F. Homand and J. P. Troalen, “Behaviour of granites and limestones subjected to slow and homogeneous temperature changes,” Engineering Geology, vol. 20, no. 3, pp. 219–233, 1984.