Determination of Constitutive Parameters of Crystalline Limestone Based on Improved RHT Model

Tian, Hao-Fan; Bao, Tai; Li, Zhuo; Peng, Hao; You, Shuai; Xiao, Shi-li

doi:https://doi.org/10.1155/2022/3794898

Advances in Materials Science and Engineering

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Mechanical Properties of Foundation Materials Under Multiple Loads

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 3794898 | https://doi.org/10.1155/2022/3794898

Determination of Constitutive Parameters of Crystalline Limestone Based on Improved RHT Model

Hao-Fan Tian,¹Tai Bao,¹Zhuo Li,²Hao Peng,¹Shuai You,³and Shi-li Xiao¹

Academic Editor: Yonghong Wang

Received31 Dec 2021

Accepted23 Mar 2022

Published13 Apr 2022

Abstract

This study improves the RHT constitutive model based on some shortcomings, and the main improvements work are as follows: (1) hyperbolic function is used to make up for the deficiency of tensile strain rate enhancement factor of RHT constitutive model at high strain rate; (2) considering the constant change of hydrostatic pressure, the initial damage value of rock is introduced into the constitutive model, and the improved damage accumulation expression is derived based on this value. (3) In the equation of residual failure surface, the influence of the Lode angle factor is considered, the RHT constitutive parameters of crystalline limestone are determined by the SHPB experiment, and the difference in rock failure forms before and after the correction of constitutive parameters is simulated by ANSYS/LS-DYNA software. The results show that the improved constitutive parameters can effectively represent the yield strength surface and residual strength surface of the rock, which verifies the rationality of the improved RHT constitutive model, and also shows that it is an effective and simple method to obtain RHT constitutive parameters by SHPB experiment.

1. Introduction

Compared with the HJC [1] and K & C [2] model, the RHT model [3] fully considers the compression effect, strain rate effect, and damage accumulation on the failure intensity of the rock under explosive impact. The elastic limit surface equation, failure surface equation, and residual strength surface equation related to pressure are embedded in this model, which is used to describe the variation law of yield strength, failure strength, and residual strength of concrete and brittle materials under high strain rate. With the rapid development of numerical simulation technology, the model is widely used in the numerical simulation of high strain rates such as explosion impact. When Tu, Hansson, and Skoglund [4–8] et al. applied the RHT model, they found that the model was incomplete for tensile damage description, so the model was modified for tensile and high strain rate effects. Zhang ruoqi and Pavlovic A [9–14] et al. adjusted and improved the failure surface equation and residual surface equation in RHT constitutive parameters based on concrete experiments and verified them by AUTODYN element simulation, all of which achieved good results. A total of 34 constitutive parameters of the RHT model need to be determined. Most scholars can only refer to Riedel’s [3, 15–17] research results because of its complexity when applying the RHT model. Li [18] used marble as an example to explain in detail the method for determining the constitutive parameters of RHT. Xie [19, 20] and Wang [21–23] et al. studied the constitutive parameters of rock under the action of high ground stress and successfully applied them to the numerical simulation analysis of rock crack propagation, cut hole blasting, and cyclic blasting. For the convenience of parameter determination, Xu and Ye[24], Liu et al. [25], and Li et al. [26–28] used an orthogonal experimental method and optimized Latin hypercube design method to analyze the sensitivity of the parameters. The analysis results show that the strength parameters of failure surface and residual surface have a sensitive influence on the fracture morphology and damage accumulation of rock. Therefore, some parameters that have little influence on rock failure can be neglected when determining rock constitutive parameters, which makes the RHT model more convenient to use.

The study of the above scholars for the RHT constitutive model mainly includes the following: the modification of the constitutive equation, the sensitivity analysis of the model parameters, and the method of determining the constitutive parameters. In the improvement of the constitutive equation, most scholars use the linear interpolation method for the residual strength of the rock in the damage softening stage, while the initial damage of rock has not been considered. In fact, for rocks in the karst area, due to the development of karst phenomena, the rock has certain initial damage, which will reduce the strength of the rock. In addition, when rock is subjected to dynamic loads such as explosion, vibration, and impact, the hydrostatic pressure on its failure surface often changes with the change of load, so it is necessary to update the expression of damage accumulation, considering that many complex experiments are usually needed in the process of determining RHT constitutive parameters. Based on this, the residual surface equation of the RHT constitutive model is improved, and the constitutive parameters of crystalline limestone are determined by the Hopkinson bar impact test.

2. RHT Constitutive Model

2.1. Failure Surface Equation

In the RHT model, the equivalent stress intensity of the failure surface is a function of normalized pressure , Lode angle , and strain rate .

In formula (1), the in the upper right corner of the formula indicates the normalized uniaxial compressive strength (unit: ) of the material, and is the compression meridian strength.

In formula (2), is the failure surface constant; is the failure surface index; and is the normalized layer crack strength. is the strain rate dynamic enhancement factor:

In formula (3), is uniaxial tensile strength (unit: )；；； is the current strain rate; is a compression strain rate index; and is the tensile strain rate index.

is the Lode angle factor: it is a function of the Lode angle and the tension-compression meridian ratio .

In formula (4), ； and are deviatoric stress invariants; and is the meridian ratio of tension and compression, and in order to consider the effect of pressure on strength, it is expressed as follows:

In formula (5), is the initial meridian ratio of tension and compression, and is the material constant.

2.2. Elastic Ultimate Surface Equation

In the RHT model, the elastic limit surface equation is derived from the failure surface equation:

In formula (6), is the elastic scaling function and is the cap function, which can effectively solve the problem that materials are always in the elastic state without yielding under high hydrostatic pressure. The expressions for , are as follows:

In equation (7), is the uniaxial tensile elastic ultimate stress, is the uniaxial compressive elastic ultimate stress, is the compressive yield surface parameter, and is the tensile yield surface parameter.

In equation (8), is the pressure at which the void material begins to crush.

2.3. Residual Strength Surface Equation

When the equivalent stress strength of the material is greater than the failure stress strength, the damage of the material begins to accumulate and enters the damage softening phase. Damage variable is the ratio of cumulative equivalent plastic strain increment to final failure equivalent plastic strain. The plastic strain of the material at ultimate failure is as follows:

In formula (9), is the minimum plastic strain at the time of material damage, and is the destruction of the cutoff pressure. and are material parameters. In fact, when the pressure of the material exceeds the yield stress, the friction and confining pressure still exist in the broken part of the material, which makes it have the ability to resist shear. Therefore, the residual strength surface equation is introduced into the RHT model.

In formula (10), is the residual stress intensity parameter and is the residual stress intensity index. The equivalent stresses between the failure strength surface and the residual strength surface are obtained by linear interpolation.

3. Improvement of the RHT Model

3.1. Correction of Strain Rate Enhancement Factor

In the RHT model, the dynamic strain rate enhancement factor (DIF) is sensitive to the high strain rate region of materials in tensile, which is different from most impact tests. Qi and Qian [29] research suggests the following: the strength of the material is not always infinitely enhanced with the increase in strain rate, when the strain rate is in the high strain rate region, and if the strain rate continues to increase, the strength of the material will slowly increase. Malvar and Ross [30] proposed a modified CEB model, which is consistent with the results of most experiments. However, the model still cannot show the characteristics of the explosion impact problem under a high strain rate. Based on this, this study uses the hyperbolic function (tanh) proposed by Gebbeken and Greulich [31] to describe the tensile strain rate enhancement factor of rock materials under a high strain rate. This function divides the strength change of rock into low strain rate region, medium strain rate region, and high strain rate region. As shown in Figure 1, the modified formula is as follows:

In formula (11), is the reference strain rate. According to the data of Tedesco and Ross [32] experiments, the fitting formula (11) gives , , , and

3.2. Model Modification Considering Initial Damage

In the RHT model, the initial damage value of brittle materials such as rock is 0 by default, but in practical engineering, the rock will be affected by different degrees of excavation disturbance, weathering, and groundwater dissolution, which leads to the initial damage value of rock not being zero. Therefore, it is necessary to introduce the initial damage value to modify the damage accumulation expression of the RHT model, and the improved damage value is as follows:

In formula (12), is equivalent plastic strain increment, is the initial injury of the rock mass, is the wave velocity of the rock mass, and is the wave velocity of the rock test piece.

3.3. Correction of Residual Strength Surfaces

The Lode angle factor is considered in the failure surface equation, which can transform the compression meridian into the tension meridian. Drawing on this idea, the rudder angle factor can be taken into account in the calculation formula of residual stress. In the RHT model, the equivalent stress strength between the failure stress surface and the residual stress surface is obtained by linear interpolation, which is inconsistent with the actual stress state in the rock mass, because the hydrostatic pressure in this interval is usually constantly changed. Therefore, the normalized hydrostatic pressure in equation (10) can be corrected to the normalized hydrostatic pressure corresponding to the current stress state surface, and finally, the corrected residual surface strength can be obtained as follows:

When the strain exceeds the strain corresponding to the yield stress, it can be assumed that there is a linear relationship between the strain increment and the plastic strain increment.

From formula (14), it can be inferred that

In equation (15), is the strain corresponding to the residual strength and is the strain corresponding to the failure strength when a brittle material such as rock loses its load-bearing capacity, and the new equivalent stress intensity between the failure stress surface and the residual stress surface can be obtained by linear interpolation:

4. Crystalline Limestone RHT Parameter Determination

4.1. Parameter Determination of the Standard RHT Model

4.1.1. Determination of Static Load Mechanic Parameters

The rock sample used in the test is crystalline limestone in a mining area in Guizhou province. The size of the specimen is 50 × 50 mm cylinder, the aspect ratio is 1.0, and core taking, cutting, and polishing of rock samples are carried out by core-taking machine, cutting machine, double-face grinding machine, and sandpaper. The nonparallelism and nonperpendicularity of rock specimens are limited below 0.02 mm in order to make the subsequent tests go on smoothly. The rock sample is shown in Figure 2.

(a)

(b)

The and of crystalline limestone can be obtained by conventional uniaxial compressive test, Brazilian splitting test, and rock wave velocity test. The rock density and initial porosity were measured by the balance method and the mass method. The mechanical parameters of crystalline limestone are shown in Table 1:

4.1.2. State Equation Parameter Determination

Brittle materials such as rock and concrete contain a large number of voids in the interior. When the material is under a strong dynamic load, the material will be subjected to the combined action of shear stress and high hydrostatic pressure, which makes the mechanical response of the material more complex. Based on this, Herrman [33] puts forward an equation of state considering the internal voids of brittle materials in 1969, which is called equation , and its compression process is shown in Figure 3:

From Figure 3, it can be seen that when the pressure value is , the initial density of the brittle material is , and when the pressure linearly increases to , then the density of the material is compacted to . When the pressure increases in a nonlinear way until the internal voids of the rock are completely compressed, the corresponding pressure is and compaction density is . Generally, the dense relationship is expressed in Mie-Greisen [34] form as follows:

In the above formula, , , and are the Rankine-Hugoniot polynomial coefficients. is the sound velocity in rock mass when the pressure is zero, and is an empirical parameter. Calcium is the main component of crystalline limestone. Available through Meyers [35] literature , , , , and are the parameters of the equation of state. According to the literature [36], the pressure at the beginning of crushing of the void can be obtained as . Generally, the pressure at the time of crushing of the material is and the compression index is

4.1.3. Determination of RHT Constitutive Parameters

Hoek and Brown [37] empirical formula is used to estimate rock strength under different confining pressures, and its expression is as follows:

and in equation (22) are the maximum and minimum principal stresses at the failure of the rock mass, respectively, reflects the soft and hard degrees of rock mass, , and reflects the degree of fragmentation of the rock mass and is taken as . , and . The mechanical parameters of the rock mass under different loading conditions can be obtained according to equation (22), as shown in Table 2.

The Mohr-Coulomb criterion was chosen to calculate the shear strength and the cohesion of the rock mass.

in equation (23), and is the internal friction angle of rock mass, which is determined by Mohr-Coulomb criterion:

The data brought into Table 2 are calculated as and .

When the rock is under quasi-static loading conditions, , , and the strength of the rock under and confining pressures is selected to determine the failure surface parameters. The equations for the compression meridian are as follows:

According to equation (25), the failure surface parameters and can be obtained. So far, most of the parameters in the constitutive model have been determined. The determination of parameters and is complicated, and the sensitivity analysis of Li et al. [26, 38] to the parameters of the RHT model shows that it has little influence on the calculation results. Therefore, the value determined by the Riedel experiment can be used. The RHT constitutive parameters of crystalline limestone are shown in Table 3.

4.2. Parameter Determination of the Modified RHT Model

4.2.1. SHPB Experiment

Hopkinson pressure bar test was carried out, and the average impact pressures of the projectile in the experiment were 0.2, 0.3, 0.4, 0.5, and 0.6 MPa, respectively. Based on the theory of one-dimensional elastic waves and the assumption of stress uniformity, the three-wave method is used to calculate the stress-strain and the average strain rates of the specimen based on the experimentally data measured, the strain rates corresponding to different shock pressures are 30.37 s⁻¹, 49.64 s⁻¹, 87.05 s⁻¹, 124.75 s⁻¹, and 138.18 s⁻¹ respectively, and the resulting stress-strain curve is shown in Figure 4.

In order to determine the elastic ultimate strength, failure strength, and residual strength of rock, it is necessary to smooth the stress-strain curves obtained from experiments. Elastic ultimate strength is the point at which the slope of the stress-strain curve begins to change, failure strength is the point at which the maximum stress value is on the stress-strain curve, and residual strength is the first turning point after the peak value of the stress-strain curve. The stress in the direction of the compression bar is , the corresponding strain is , and the strain rate is . The strength indexes are shown in Table 4.

It can be seen from Figure 3 that when the impact pressure is 0.2 MPa, the corresponding yield strength is 61.02 MPa, and with the gradual increase in impact pressure, the failure strength of the rock also gradually increases. When the impact pressure is 0.6 MPa, the failure strength of rock reaches 148.56 MPa, which is 2.43 times higher than the yield strength when the impact pressure is 0.2 MPa and 3.36 times higher than the uniaxial compressive strength. It can be seen that with the increase in impact pressure, the yield strength of the rock is also gradually improved.

Based on the experimental results in Table 4, it is possible to plot the relationship between and , the curve equation is fitted by the modified residual surface equation of equation (13), and decision factor . The fitting curve is shown in Figure 4:

From Figure 5 we can obtain and , and are calculated by taking corresponding to 0.2 MPa and 0.6 MPa, respectively, and the rest of the parameters are the same as those taken in Table 3.

5. Validation of the Modified RHT Model

5.1. Numerical Simulation

In order to verify the correctness of the parameters of the improved RHT model, the numerical model is established by ANSYS/LS-DYNA software according to the experimental process of 1 : 1, in which the incident rod is 2.0 m long, the transmission rod is 1.5 m long, and the diameter of the rod is 0.05 m. The keyword LOAD_SEGMENT_SET is used to load the sinusoidal stress wave on the end face of the incident rod, which not only ensures the consistency of the waveform but also plays a role in simplifying the model. The model is meshed into 255,400 cells and 274,373 nodes using SOLID 164 cells. The contact between the bar and the specimen is set as surface-to-surface erosion contact, the contact stiffness is the default value, and the dynamic and static friction coefficients are zero. The finite element analysis model is shown in Figure 6. The incident stress wave curve is selected when the impact pressure is 0.2 MPa for loading. When , the peak value of incident stress is 134.26 MPa, and the stress wave curve is shown in Figure 7.

(a)

(b)

5.2. Analysis of Results

In order to verify the rationality of the improved RHT model, it is necessary to compare the difference between the rock failure patterns before and after the modification of constitutive parameters in numerical simulation and the experimental results, and the keyword^∗MAT_ADD-EROSION in LS-DYNA is used to control the failure of the element in order to simulate the crack propagation in rock. The results of numerical simulation before correction and the experimental impact effect are shown in Figure 8.

As can be seen from Figure 8, the stress-strain curves obtained by numerical simulation show that the elastic ultimate strength is 6.98 MPa higher than the experimental value, and the yield strength is 1.45 MPA higher than the experimental value. At the same time, it is not difficult to find that the strain value corresponding to the elastic ultimate strength in the numerical simulation and test results is 5 × 10⁻⁴ more than the test value. When the strength exceeds the yield value, the residual strength surface described by the numerical simulation results is different from the experimental results. Although the surface crack propagation trend of numerical simulation is similar to the experimental results, only a few elements fail on the side of rock samples. This is because the unmodified constitutive parameters are conservative in describing the yield strength surface and residual strength surface, which leads to the failure mode of rock not fully displayed. In order to make the simulation results more consistent with the failure mode of rock, the improved constitutive parameters are used for calculation, and the results are shown in Figure 9.

From Figure 9, it can be found that the description of elastic ultimate strength surface, yield strength surface, and residual failure surface by the improved RHT model is basically consistent with the test results. This is because the modified RHT model considers the constant change of hydrostatic pressure in the failure surface equation and the influence of the Lode angle factor in the residual strength surface equation, which makes the constitutive parameters of the RHT model more reasonable to describe the failure strength surface and the residual strength surface, which also verifies the correctness of the modified RHT model.

6. Conclusions

In this study, the shortcomings of the RHT constitutive model are studied, and the conclusions are as follows:(1)It is practical to use the hyperbolic function to describe the tensile strain rate enhancement factor of rock at high strain rates.(2)The introduction of the initial damage variable makes the expression of damage variable in the RHT constitutive model more complete.(3)Considering the influence of the Lode angle factor and hydrostatic pressure, the updated damage variable can better reflect the actual failure mode of the rock mass.(4)It is a simple and effective method to determine the constitutive parameters of the RHT model by SHPB experiment.

Data Availability

All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest. Tian Hao-fan and Bao Tai equally contributed to this manuscript.

Acknowledgments

The author’s research work was supported by the National Natural Science Foundation of China (no. 51664007).

References

G. R. Johnson and W. H. Cook, “A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures,” Engineering Fracture Mechanics, vol. 21, pp. 541–548, 1983.
View at: Google Scholar
L. J. Malvar, J. E. Crawford, J. W. Wesevich, and D. Simons, “A plasticity concrete material model for DYNA3D,” International Journal of Impact Engineering, vol. 19, no. 9-10, pp. 847–873, 1997.
View at: Publisher Site | Google Scholar
W. Riedel, K. Thoma, and S. Hiermaier, “Penetration of reinforced concrete by BETA-B-500 numerical analysis using a new macroscopic concrete model for hydrocodes,” in Proceedings of the 9th International Symposium, Interaction of the Effects of Munitions with Structures, pp. 315–322, IBMAC, Berlin, Germany, 1999.
View at: Google Scholar
P. Skoglund, “Simulation of concrete penetration in 2D and 3D with the RHT material model,” Swedish Defense Research Agency, 2002.
View at: Google Scholar
Z. Tu and Y. Lu, “Evaluation of typical concrete material models used in hydrocodes for high dynamic response simulations,” International Journal of Impact Engineering, vol. 36, no. 1, pp. 132–146, 2009.
View at: Publisher Site | Google Scholar
Z. Tu and Y. Lu, “Modifications of RHT material model for improved numerical simulation of dynamic response of concrete,” International Journal of Impact Engineering, vol. 37, no. 10, pp. 1072–1082, 2010.
View at: Publisher Site | Google Scholar
T. L. Ling, Y. T. Wang, and D. S. Liu, “Application of modified RHT model in numerical simulation of rock blast response,” Journal of Coal, vol. 43, no. S2, pp. 434–442, 2018.
View at: Google Scholar
C. P. Kuang and X. K. Yuan, “Analysis and simulation of strength parameters of RHT concrete principal structure model,” Quarterly Journal of Mechanics, vol. 33, no. 1, pp. 158–163, 2012.
View at: Google Scholar
M. Abdel-Kader, “Modified settings of concrete parameters in RHT model for predicting the response of concrete panels to impact,” International Journal of Impact Engineering, vol. 132, pp. 103312.1–103312.18, 2019.
View at: Publisher Site | Google Scholar
A. Pavlovic, C. Fragassa, and A. Disic, “Comparative numerical and experimental study of projectile impact on reinforced concrete,” Composites Part B: Engineering, vol. 108, pp. 122–130, 2017.
View at: Publisher Site | Google Scholar
D. S. Shi, G. Y. Shi, and J. Y. Chen, “Improvement and numerical fitting of RHT model,” Journal of Ningbo University (Science and Technology Edition), vol. 27, no. 03, pp. 93–96, 2014.
View at: Google Scholar
Y. Wang, Research on Rock Blasting Fragmentation Model Based on RHT Ontology, China University of Mining and Technology, Beijing, China, 2015.
Z. Wang, Y. Ni, and J. Cao, “Progress in the study of dynamic mechanical properties of concrete under impact loading,” Explosion and Impact, vol. 2005, no. 6, pp. 519–527, 2005.
View at: Google Scholar
R. Zhang, Y. Ding, and W. Tang, “Failure strength parameters of concrete HJC and RHT intrinsic structure models,” Journal of High Pressure Physics, vol. 25, no. 1, pp. 15–22, 2011.
View at: Google Scholar
A. By, A. Zl, and A. Yc, “Mechanical and microstructural properties of recycling granite residual soil reinforced with glass fiber and liquid-modified polyvinyl alcohol polymer,” Chemosphere, 2021.
View at: Google Scholar
L. Wang, F. Guo, H. Yang, Y Wang, and S Tang, “Comparison of fly ash, PVA fiber, MGO and shrinkage-reducing admixture on the frost resistance of face slab concrete via pore structural and fractal analysis,” Fractals, vol. 29, Article ID 2140002, 2020.
View at: Google Scholar
L. Wang, R. Luo, W. Zhang, M. Jin, and S. Tang, “Effects of fineness and content of phosphorus slag on cement hydration, permeability, pore structure and fractal dimension of concrete,” Fractals, vol. 29, Article ID 2140004, 2020.
View at: Publisher Site | Google Scholar
H. Li, Research on the Theory of Rock RHT Model and the Method of Determining the Main Parameters, China University of Mining and Technology, Beijing, China, 2016.
L. X. Xie, W. B. Lu, Q. B. Zhang, Q. H. Jiang, M. Chen, and J. Zhao, “Analysis of damage mechanisms and optimization of cut blasting design under high in-situ stresses,” Tunnelling and Underground Space Technology, vol. 66, pp. 19–33, 2017.
View at: Publisher Site | Google Scholar
C. Pan, L-X. Xie, X. Li, K. Liu, P.-F. Gao, and L.-G. Tian, “Numerical investigation of the effect of eccentric decoupled charge structure on blasting-induced rock damage,” Journal of Central South University, vol. 29, pp. 1–16, 2022.
View at: Google Scholar
Z. Wang, H. Wang, J. Wang, and N. Tian, “Finite element analyses of constitutive models performance in the simulation of blast-induced rock cracks,” Computers and Geotechnics, vol. 135, p. 135, 2021.
View at: Publisher Site | Google Scholar
B. Yuan, Z. Li, Z. Zhao, H. Ni, Z. Su, and Z. Li, “Experimental study of displacement field of layered soils surrounding laterally loaded pile based on transparent soil,” Journal of Soils and Sediments, vol. 21, no. 4, pp. 1–12, 2021.
View at: Publisher Site | Google Scholar
B. Yuan, M. Sun, Y. Wang, L. Zhai, Q. Luo, and X. Zhang, “Full 3D displacement measuring system for 3D displacement field of soil around a laterally loaded pile in transparent soil,” International Journal of Geomechanics, vol. 2019, no. 5, Article ID 04019028, 2019.
View at: Google Scholar
C. Xu and G. Ye, “Application of orthogonal experimental design for sensitivity analysis of numerical model parameters,” Hydrogeology & Engineering Geology, vol. 4, no. 01, pp. 95–97, 2004.
View at: Google Scholar
D. Liu, N. Liu, and T. Gao, “Parameter sensitivity study of RHT model by applying orthogonal test method,” Journal of Beijing University of Technology, vol. 39, no. 06, pp. 558–564, 2019.
View at: Google Scholar
H. Li, Y. Chen, and D. Liu, “Study on the sensitivity of main parameters and determination method of rock RHT model,” Journal of Beijing University of Technology, vol. 38, no. 08, pp. 779–785, 2018.
View at: Google Scholar
Z. Nie, Y. Peng, and R. Chen, “Parameter sensitivity analysis of RHT model for rock-like materials under intrusion conditions,” Vibration and Shock, vol. 40, no. 14, pp. 108–116, 2021.
View at: Google Scholar
X. Jian, “Sensitivity analysis of RHT model parameters under the effect of explosion,” Naval Electronics Engineering, vol. 39, no. 04, pp. 111–113+122, 2019.
View at: Google Scholar
C. Qi and Q. Qian, “Physical mechanism of strain rate dependence of dynamic strength of brittle materials such as rocks,” Journal of Rock Mechanics and Engineering, vol. 3, no. 02, pp. 177–181, 2003.
View at: Google Scholar
L. J. Malvar and C. A. Ross, “Review of strain rate effects for concrete in tension,” Materials Journal, vol. 95, no. 6, 1998.
View at: Google Scholar
N. Gebbeken and S. Greulich, “A new material model for SFRC under high dynamic loadings,” in Proceedings of the 1 Ith International Symposium on Interaction of the Effects of Munitions with Structures, Mannheim, Germany, 2003.
View at: Google Scholar
J. W. Tedesco, J. C. Powell, and C. A. Ross, “A strain-rate-dependent concrete material model for ADINA,” Computers & Structures, vol. 64, no. 5, 1997.
View at: Publisher Site | Google Scholar
W. Herrmann, “Constitutive equation for the dynamic compaction of ductile porous materials,” Journal of Applied Physics, vol. 40, no. 6, pp. 2490–2499, 1969.
View at: Publisher Site | Google Scholar
X. Zhang, The Point-of-Matter Method, Point-of-matter method, 2013.
M. A. Meyers, Dynamic Behavior of Materials, Dynamic Behavior of Materials, New Jersey, NJ, USA, 1994.
H. G. R. A. Johnson, “Computational constitutive model for glass subjected to large strains, high strain rates and high pressures,” Journal of Applied Mechanics, vol. 78, no. 5, Article ID 051003, 2011.
View at: Google Scholar
E. Hoek and E. T. Brown, Underground Excavations in Rock, Underground Excavations in Rock, London, UK, 1980.
M. Yu, L. He, and L. Song, “Twin shear stress theory and its generalization,” Ence in China Ser A, vol. 10, 1985.
View at: Google Scholar

Copyright

Copyright © 2022 Hao-Fan Tian 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.

PDF Download Citation

Download other formats

Order printed copies