Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Additively Manufactured Continuous Processing Reactor System for Producing Liquid-Based Pharmaceutical Substances
Previous Article in Journal
Theory of Refraction, Ray–Wave Tilt, Hidden Momentum, and Apparent Topological Phases in Isotropy-Broken Materials Based on Electromagnetism of Moving Media
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dynamic Damage Constitutive Model of Rock-like Materials Based on Elastic Tensile Strain

Northwest Institute of Nuclear Technology, Xi’an 710024, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(16), 6852; https://doi.org/10.3390/app14166852
Submission received: 17 June 2024 / Revised: 27 July 2024 / Accepted: 2 August 2024 / Published: 6 August 2024

Abstract

:
To accurately characterize the damage of rock-like materials under simultaneous or alternating tensile and compressive loading, a dynamic damage constitutive model for rock-like materials based on elastic tensile strain is developed by integrating the classical compressive plastic damage model and the tensile elastic damage model. The model is based on the Holmquist–Johnson–Cook (HJC) and Kuszmaul (KUS) models, categorizing the element stress state into tensile and compressive states through positive and negative elastic volumetric strain. It utilizes elastic tensile strain to enhance the calculation method for tensile cracks, determining the tensile strength of the principal direction based on the contribution rate of tensile principal stress for uniaxial/multiaxial loading. Additionally, it establishes a maximum elastic tensile strain rate function to rectify the model’s effect on the tensile strain rate. Through the LS-DYNA subroutine development, the model proficiently delineates the distribution of ring-shaped cracks on the frontal side and strip-shaped cracks on the rear side of the reinforced concrete slab subjected to impact loading. Numerical simulations demonstrate that the model provides more accurate damage prediction results for stress conditions involving simultaneous or alternating compression and tension, offering valuable insights for damage analysis in engineering blasting or impact penetration.

1. Introduction

The dynamic response process of damage and fracture in rock-like materials under blast loading is highly complex, and the study of the dynamic constitutive model of these materials, including rock and concrete, represents a significant research focus. Numerous experimental studies have demonstrated that rock-like materials exhibit significant tensile brittleness, the tensile–compressive strength difference effect, the damage effect, the confining pressure effect, and the strain rate effect. To accurately characterize the dynamic mechanical properties of rock-like materials, various classical dynamic damage constitutive models have been proposed in the industry. These can be broadly categorized into two groups: one being the compressive-shear plastic damage model, exemplified by the Holmquist–Johnson–Cook (HJC) model [1], the Riedel–Hiermaier–Thoma (RHT) model [2], and the Karagozian & Case (K&C) model [3], primarily employed in applications such as impact penetration; the other is the tensile statistical damage model, represented by the Grady–Kipp (GK) model [4], the Taylor–Chen–Kuszmaul (TCK) model [5], and the Kuszmaul (KUS) model [6], mainly utilized in fields like engineering blasting. In recent years, several scholars have conducted extensive research and made significant improvements regarding the applicability of classical models in specific engineering calculations. Xiong et al. [7,8] contended that the HJC model demonstrates superior performance in shear-compression damage description, computational efficiency, and convenience compared to the RHT model and the K&C model, and refined its parameter determination method. Du et al. [9] conducted an analysis and concluded that the HJC model offers a more precise depiction of damage and maximum deflection in comparison to the experimental results. Jin [10] integrated the TCK model for tensile damage with the RHT model for shear damage to formulate a comprehensive damage model applicable to concrete impact or penetration scenarios; the calculated results indicate that the model is capable of predicting crater formation, perforation, and post-perforation scabbing on the inner face of the tunnel and the back face of the target during penetration into concrete or reinforced concrete materials, and the residual velocity of the penetrator aligns well with experimental data. Wu et al. [11] developed a novel damage model, utilizing the TCK model and the HJC model, to characterize the explosive damage effects in concrete; the numerical simulation results depict the failure of reinforced concrete and the formation of a blasting crater, specifically highlighting the compression of the central blasting cavity and the spalling due to free surface tension, which aligns with both experimental and theoretical findings. Qiang et al. [12] refined the treatment of the compression–tension interface and introduced the TCK-HJC composite constitutive model to characterize the penetration and perforation behavior of thin plates; the calculated results indicate that the dispersion of failed concrete fragments under two distinct penetration angles closely aligns with the experimental findings. Hu et al. [13] compared five typical blast damage models and selected the KUS tensile model with superior computational result while also incorporating the Rheological–Dynamical Analogy (RDA) compression damage model to develop a novel combined tensile–compression damage model for characterizing rock blast damage, the numerical simulation results demonstrate the model’s capability to accurately depict the blasting-induced damage in rock. Chen et al. [14] and Xie et al. [15] both enhanced the equivalent tensile fracture damage model proposed by Yang et al. [16] and developed their constitutive models for tensile–compressive damage. The former’s model correlates compressive damage with compressive plastic work and the history of tensile damage, while the latter’s model still associates compressive damage with equivalent tensile strain. These studies demonstrate that directly integrating a compressive damage model and a tensile damage model into a new framework is a scientifically sound, efficient, and practical approach based on classical models.
Building on the aforementioned concepts, this study endeavors to integrate the widely acknowledged HJC model and KUS model to formulate a composite damage constitutive model known as the HJC-KUS (HK) composite damage model. This new model combines the attributes of the HJC model for characterizing distortion/volumetric plastic damage under shear and compression, along with those of the KUS model for describing microcrack extension damage under volumetric tensile deformation. In the course of this investigation, several typical issues are identified in the direct composite model (HK0), including the influence of compression plastic history on tensile strength, excessive suppression of tensile damage under volumetric tensile conditions, and discrepancies between uniaxial/multiaxial tensile strength and experimental observations. Consequently, targeted improvements are implemented in this study to obtain the HK model, and the HK model is integrated into LS-DYNA for further validation.

2. The Construction of the HK Model

2.1. Damage Evolution

2.1.1. Evolution of Compressive Damage

Holmquist et al. [1] proposed that under strong compressive loading, the accumulation of volumetric plastic strain and equivalent (shear) plastic strain leads to material strength degradation or failure. The compressive damage Dc is defined as in Equation (1).
D c = Δ ε p + Δ u p ε f p + u f p
ε f p + u f p = D 1 ( p * + T * ) D 2
In these equations, Δ ε p and Δ u p represent the increments in equivalent plastic strain and volumetric plastic strain at the current time step, D1 and D2 are material parameters, p * = p / f c is the normalized pressure (p denotes hydrostatic pressure, fc represents uniaxial compressive strength), and T * = T / f c is the normalized maximum tensile hydrostatic pressure (T denotes the maximum tensile hydrostatic pressure).

2.1.2. Evolution of Tensile Damage

Kuszmaul et al. [6] proposed that when a material is subjected to volumetric tensile loading, internal cracks will be activated, and the relationship between the tensile damage Dt, the Poisson’s ratio ν , crack density Cd and its evolution rate can be expressed as in Equations (3)–(5).
D t = 16 9 × 1 ν ¯ 2 1 2 ν ¯ C d
ν ¯ = ν ( 1 16 9 C d )
C ˙ d = 5 2 k m K I C ρ 0 V μ ˙ max 2 μ m 1 μ ˙ ( 1 D t )
In these equations, k and m represent material parameters, μ denotes the volumetric strain acting in tension (non-negative), and μ ˙ max is the maximum volumetric strain rate. The symbol ρ 0 represents density, V stands for sound velocity, and KIC signifies the first-type fracture toughness. The crack density Cd represents the ratio of the volume of the cracked region to the total volume, and it increases with more severe damage.

2.1.3. Issues and Resolutions

Referring to the constitutive coupling method [12], the zero value of elastic volumetric strain is utilized as the classification criterion for compressive and tensile states, enabling direct coupling of the compressive damage evolution of the HJC model with the tensile damage evolution of the KUS model. This results in the formulation of a damage constitutive model (HK0 model) expressed by Equations (1)–(5). Nevertheless, there are numerous issues associated with directly employing this model for numerical calculations. To enhance its applicability, this paper addresses the following problems to develop the final HK model.
  • The impact of accumulated compressive plasticity.
In the context of explosion and layered fracture, compressive loading precedes tensile loading. If the material undergoes plastic deformation under compression during this process, a discrepancy may arise between the zero value of hydrostatic pressure and the zero value of volumetric strain when transitioning from compression to tension:
p = 0 μ e = 0 , μ p < 0 , μ = μ e + μ p < 0
In this equation, μ e and μ p represent the elastic and plastic strains of the material, respectively. Subsequently, the strains are increased as follows:
μ e = μ p > 0 μ = μ e + μ p = 0 p < 0
During the loading process from p = 0 and μ < 0 to p < 0 and μ = 0, the tensile stress will continue to increase based on elastic calculations without accounting for damage. The damage evolution relationship is unable to cover all stress states, ultimately resulting in a significantly higher calculated tensile strength than the actual/input strength.
Since the sign of hydrostatic pressure is consistent with the sign of the elastic body strain in the elastic state, this study exclusively employs the elastic component of the volumetric strain to assess tensile damage, aiming to avoid the impact of compressive plastic loading history on both tensile damage and strength.
2.
Conservative predictions of tensile damage.
Owing to the evident tensile brittle characteristics of rock-like materials, their tensile strength is approximately one order of magnitude lower than their compressive strength. Under compression with p > 0, the material may still undergo tensile failure, as observed in Brazilian split testing, where an approximate relationship is as shown in Equation (8) [17].
σ y min = 3 σ x max
In this equation, σ y min represents the minimum compressive stress in the direction of the load, and σ x max represents the maximum tensile stress in the perpendicular direction of the load. Currently, the hydrostatic pressure of the adjacent unit near the specimen center is positive, and the volumetric strain of the adjacent unit near the specimen center is negative. According to Equation (3), it exclusively computes crack density and tensile damage when μ > 0. Despite Equation (3) indicating that no tensile damage will occur at the center of the specimen under Brazilian split loading, in reality, tensile failure typically occurs at this location. Therefore, Equation (3) yields an excessively conservative prediction for tensile damage.
Yang et al. [16] posited that the equivalent tensile strain is the primary cause of material damage, and utilizing this measure to calculate damage could broaden the conditions under which material damage occurs. The equivalent tensile strain θ is defined as depicted in Equation (9), as described by the authors.
θ = i = 1 3 ( ε i + ε i ) / 2
In this equation, ε i denotes the strain in the three principal directions. While directly substituting the equivalent tensile strain θ for the volumetric strain μ in Equation (2) can broaden the conditions for calculating tensile damage, it results in excessively lenient criteria. In the elastic state, the strain in the principal direction comprises both the strain induced by tensile stress and that caused by the Poisson’s effect as shown in Equation (10).
ε i e = 1 E σ i + ν E ( σ i σ k k ) , i = 1 , 2 , 3
In this equation, ε i e denotes the elastic strain in the three principal directions, σ i represents the stress in these directions, and E stands for Young’s modulus. In an elastic state, the Poisson effect causes tensile and compressive strains in one principal direction to produce strains in the other principal directions. This causes Equation (9) to include an unnecessary Poisson effect that influences in its calculated value, as observed in the case of the uniaxial compression condition where the following relationships exist in Equation (11).
ε 1 = ε 2 = ν ε 3 > ε 3 μ = ( 1 2 ν ) ε 3 < 0 θ = 2 ν ε 3 > 0
In this scenario, substituting the equivalent tensile strain θ for the volumetric strain μ in Equation (5) will result in the calculation of tensile cracking, contradicting the expectation that the material should experience compressive damage rather than tensile fracture under uniaxial compression. Utilizing equivalent tensile strain for the calculation of tensile crack density may lead to an overestimation in damage assessment.
This study posits that the elastic strain resulting from tensile stress in the principal direction of the strain is responsible for causing tensile damage in materials, rather than the additional strain induced by the Poisson effect. In accordance with the method of equivalent tensile strain, the elastic tensile strain ϕ e is defined as in Equation (12).
ϕ e = i = 1 3 σ i + σ i / 2 E
Ultimately, the substitution of elastic tensile strain for volumetric strain is employed to enhance the crack density calculation method, thereby expanding the scope of tensile damage assessment. The crack density evolution rate function is shown in Equation (13).
C ˙ d = 5 2 k m K I C ρ 0 V ϕ ˙ max e 2 ϕ e m 1 ϕ ˙ e ( 1 D t )
3.
Prediction of anomalies in multi-axial tensile strength.
A substantial body of experimental evidence suggests that the tensile strength of rock-like materials remains consistent under uniaxial, biaxial, and triaxial loading conditions [18]. However, most existing models fail to account for this behavior; for instance, the HJC model predicts a tensile strength with ft/3 = 2ftt/3 = fttt, while the KUS model also exhibits similar discrepancies in predicting tensile strength under uniaxial/multiaxial loading conditions. It is imperative to consider the evolution of tensile damage in different principal directions to enhance the characterization of tensile strength within the model. Therefore, the principal stress crack density contribution factor λ i and the increment in crack density in the principal direction C ˙ d i are introduced, as shown in Equations (14) and (15).
λ i = σ i + σ i i = 1 3 σ i + σ i , i = 1 , 2 , 3
C ˙ d i = C ˙ d λ i m , i = 1 , 2 , 3
By combining Equations (3) and (4), the tensile damage Dti in various principal directions can be determined. The tensile damage accumulates separately in the three principal directions, and the maximum value is taken as the final tensile damage:
D t = max D t 1 , D t 2 , D t 3
4.
The impact of tensile strain rate.
Experimental results indicate a strain rate dependency in the tensile strength of rock-like materials [19], with a weaker effect at low strain rates and a more pronounced impact at high strain rates. To accurately characterize the strain rate effect of materials, an auxiliary strain rate ϕ ˙ 0 is introduced to establish the lower limit of the strain rate parameter ϕ ˙ max e , thus achieving a segmented depiction of the tensile strain rate effect. ϕ ˙ max e is shown in Equation (17).
ϕ ˙ max e = ϕ ˙ max e ϕ ˙ max e ϕ ˙ 0 ϕ ˙ 0 ϕ ˙ max e < ϕ ˙ 0
Therefore, the damage calculation expressions of the HK model are as follows:
D t i = 16 9 × 1 - ν ¯ 2 1 - 2 ν ¯ C d i ϕ e > 0 D t = max D t 1 , D t 2 , D t 3 D c = Δ ε p + Δ u p ε f p + u f p μ e 0
The damage variable D of the material is defined as the maximum of Dt and Dc, taking into consideration the most adverse effects of damage.

2.2. Constitutive Relationship

2.2.1. Compressive Constitutive Relationship

The HJC model is used to describe the elastic–plastic response of materials under compression, and the distortion behavior of the material is represented by an equivalent strength function, as shown in Equation (19).
σ y * = [ A ( 1 D c ) + B p * n ] ( 1 + C ln ε ˙ * )
In this equation, σ y * represents the normalized equivalent intensity, and ε ˙ * denotes the dimensionless equivalent strain rate, while A, B, and n are parameters defining the intensity limit surface, and C stands for the strain rate coefficient. When the equivalent stress surpasses the equivalent strength, the material undergoes plastic deformation, and the plastic strain is determined using the corresponding flow law.
The volumetric deformation behavior of materials is characterized by a three-stage state equation, as shown in Equation (20).
p = K μ e e l a s t i c   s t a g e p c r u s h + K t r a n ( μ μ c r u s h ) t r a n s i t i o n a l   s t a g e K 1 μ ¯ + K 2 μ ¯ 2 + K 3 μ ¯ 3 c o m p a c t i o n   s t a g e
In this equation, p c r u s h and μ c r u s h represent the pressure and volumetric strain at which the cavity initiates collapse, K denotes the bulk modulus, Ktran represents the bulk modulus of the transitional stage, while K1, K2, and K3 are material constants, and μ ¯ is the modified volumetric strain.

2.2.2. Tensile Constitutive Relationship

According to the KUS model, the tensile state is described by a linear elastic damage constitutive equation, as shown in Equation (21), which causes a decrease in modulus due to tensile damage.
p = K 1 D μ e S i j = 2 G 1 D e i j
In this equation, K and G are the initial bulk modulus and shear modulus, respectively, and e i j is the deviatoric strain, while S i j is the deviatoric stress.
Equations (18)–(21) constitute the theoretical basis of the HK model.

3. The HK Model Testing

The HK model theory is implemented in the commercial software LS-DYNA to conduct simple and comprehensive loading numerical calculations aiming to validate its rationality.

3.1. Material Parameter Determination

The parameters for the compression stage can be obtained by referring to the parameter fitting method of the HJC model [7]. The parameters for the stretching stage include KIC, k, and m. KIC can be derived from experimental data. The parameter m for granite is commonly assigned as 7 [5,13], whereas for concrete, it is typically set at 6 [10,11,12]. The parameter k is fitted based on the uniaxial tensile strength of the material. By employing the aforementioned approach, the material parameters for the HK model of high-strength steel fiber concrete are derived and presented in Table 1.

3.2. Simple Loading Numerical Calculations

3.2.1. Numerical Calculation of a Single-Element Model

The numerical calculation of a single-element model is hypothetical. A single-element model is established, and numerical simulations are performed using the HK and HK0 models for pre-compression followed by tensile loading, uniaxial and triaxial equal loading, as well as tensile loading at different strain rates. The material parameters are detailed in Table 1, with corresponding results presented in Figure 1. In Figure 1a, the initiation of damage occurs earlier in the HK model as compared to the HK0 model. The hydrostatic tensile strength calculated by the HK model aligns with the input value of parameter T, while that computed by the HK0 model significantly exceeds the specified strength. In Figure 1b, the uniaxial and triaxial tensile strengths determined by the HK model are identical, whereas substantial disparities exist in those calculated by the HK0 model. In Figure 1c, under loading conditions with a strain rate greater than the reference, both material strengths computed by the HK and HK0 models exhibit rapid escalation with an increasing strain rate. Conversely, under loading conditions with a strain rate lower than the reference, while material strength remains constant for the HK model, it declines sharply for the HK0 model.
The computational results indicate the following: (1) When subjected to plastic compression followed by reverse stretching, the HK model effectively accounts for the influence of accumulated plasticity on damage calculation. Additionally, it synchronously calculates tensile damage and stress without encountering situations where the evolution relationship of damage fails to cover tensile stress as observed in the HK0 model calculations, as shown in Figure 1a. (2) In comparison to the tensile strength predicted by the HK0 model for uniaxial/triaxial tensile tests, the tensile strength predicted by the HK model is found to be consistent, aligning with established experimental findings, as shown in Figure 1b. (3) The HK model can reasonably describe the strain rate effect of the segmented form of the material, and the tensile strength obtained under low strain rate loading is more reasonable than that obtained using the HK0 model, as shown in Figure 1c.

3.2.2. Numerical Computation of Split Test

The numerical computation of the split test is a hypothetical calculation. The split test is a classical experimental technique used for indirectly measuring the tensile strength of rock-like materials. In this test, the units near the center of the specimen experience both tensile and compressive stresses simultaneously. Performing split calculations enables the analysis of the damage behavior of the specimen under combined tensile and compressive loading conditions. The specimen is sized to be Φ5.04 cm × 2.4 cm, and the displacement load is applied to the specimen through the opposite motion of two steel plates, with a loading strain rate of 5/s. The numerical computations are conducted using the HK and HK0 constitutive models, with the material parameters detailed in Table 1.
The numerical results from the split loading simulations (Figure 2a,b) indicate that in the HK model, the damage distribution is primarily concentrated at the loading location and in the vertical direction, exhibiting a central crack consistent with the typical experimental failure pattern (Figure 2c). In contrast, for the HK0 model, damage distribution is solely focused on the loading location without any damage occurring at the center of the specimen, thus failing to simulate tensile failure behavior.

3.3. Comprehensive Loading Numerical Calculation

Wang et al. studied the damage results of reinforced concrete slabs through explosion experiments [20]. In the experiment, the dimensions of the reinforced concrete slab were 1.25 m × 1.25 m × 0.05 m, with a single layer of bi-directional φ6@75 mm reinforcement. The slab was constrained by clamp plates on both sides. The explosive was located directly above the center of the slab, 50 cm away from the slab, and the explosive charge was 640 g. Based on the experimental conditions, a numerical simulation was conducted to analyze the impact of an explosion on a reinforced concrete slab. The complete model was utilized for the numerical computation, with fully constrained vertical and rotational displacements on both sides of the slab and a joint connection between the steel bar and the concrete. A Lagrange mesh with a grid size of 0.5 cm was applied for the reinforced concrete slab, while a Euler mesh with a grid size of 1 cm was used for explosives and air. Additionally, material interaction was implemented using the structured–arbitrary Lagrange Euler method. The material parameters of the concrete are detailed in Table 2.
The numerical calculations for the front and back surfaces of the thin plate are shown in Figure 3. In the experimental findings [20], a small crack is observed at the front center of the slab, along with a visible curved crack running along the symmetry centerline and several circular cracks on the surface. At the rear of the slab, there are pull-off pits, circular cracks, and radial cracks extending outward from the center. Compared with the experimental findings, the computed results of the HK model demonstrate the formation of ring-shaped cracks, partial vertical fractures, and damage at the central position on the front surface of the slab. Additionally, damage is observed on the back surface, including a central area affected, ring-shaped cracks, and cracks extending outward from the center. The simulation effectively characterizes the damage morphology of both surfaces of the slab.
Du et al. compared the damage calculation results of the HJC model, RHT model, and K&C model based on the experiment [9]. On the front side of the slab, the damage calculation results of the HJC model indicate a lateral damage zone caused by the central curved deformation, and an approximately circular damage area at the center with a total damage value of less than 0.3, while the damage calculation results of the RHT model show lateral bending damage at the center and simulated radial damage cracks around it. Meanwhile, for the K&C model, extensive surface damage is observed. On the rear side of the plate, the simulation results of the HJC model do not exhibit any significant damage phenomena different from the other models, while the simulation results of the RHT model show cross-crack patterns and radial cracks on the surface, and the simulation results of the K&C model also show extensive surface damage. The degree of similarity of the damage distribution of the HK model to the experimental results is significantly superior to that attained by employing the HJC model, RHT model, and K&C model.

4. Conclusions

This study adopts the constitutive coupling method to establish the HK model capable of describing rock-like material damage evolution under impact loading, based on the HJC compressive damage model and KUS tensile damage model. Specific improvements have been made to the model to enhance its applicability, including:
  • Replacing volumetric strain with elastic tensile strain to enhance the method for calculating tensile crack density, thus overcoming the influence of plastic accumulation on tensile stress calculation and expanding the range of damage assessment under tensile loading,
  • Calculating the crack density on the principal direction based on the contribution rate of the principal stress to unify the tensile strength under uniaxial and multiaxial loading conditions,
  • Introducing a maximum elastic tensile strain rate function to make the description of the material strain rate effect more accurate.
Based on this, the HK model is imported into the finite element software and verified through calculations for both simple and combined scenarios. Some advantages are demonstrated by the computational results. The HK model overcame some inherent limitations of creating the HK0 model by directly combining the HJC and KUS models, thus avoiding overly conservative estimates of tensile damage, which helps reduce safety risks related to the design. The HK model has significantly higher accuracy than the HJC, RHT, and K&C models in predicting damage in complex loading conditions involving tensile and compressive loads. It is suggested that the HK model proposed in this study holds practical utility for predicting damage in engineering applications such as blasting or impact penetration.
However, the HK model established in this paper only utilized elastic tensile strain for calculating tensile damage. Therefore, the model is unable to depict the shear failure process under a tensile state. Subsequent enhancements will be implemented to refine the calculation method for shear damage under a tensile state, thereby improving the model’s capacity to characterize damage in rock-like materials. Concurrently, verification experiments on various rock-like materials will be conducted to further validate the rationality of the model.

Author Contributions

Conceptualization, X.Z. and Y.X.; methodology, X.Z., Y.Z., W.W. and F.Z.; validation, X.Z., Y.X. and L.W.; resources, X.Z. and Y.Z.; data curation, X.Z. and L.W.; writing—original draft preparation, X.Z.; writing—review and editing, X.Z., Y.X., W.W. and F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in this study are included in this article; further inquiries can be directed to the corresponding authors.

Acknowledgments

The authors would like to thank Li-peng Hu of the Northwest Institute of Nuclear Technology for his assistance in the split experiment.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Holmquist, T.; Johnson, G.; Cook, W. A computational constitutive model for concrete subjected to large strains, high strain rates and high pressures. In Proceedings of the 14th International Symposium on Ballistics, Quebec, QC, Canada, 26–29 September 1993; pp. 591–600. [Google Scholar]
  2. Riedel, W.; Thoma, K.; Hiermaier, S. Penetration of reinforced concrete by BETA-B-500 numerical analysis using a new macroscopic concrete model for hydrocodes. In Proceeding of the 9th International Symposium, Interaction of the Effects of Munitions with Structures, Berlin, Germany, 3–7 May 1999; pp. 315–322. [Google Scholar]
  3. Malvar, L.; Crawford, J.; Wesevich, J. A plasticity concrete material model for DYNA3D. Int. J. Impact Eng. 1997, 19, 847–873. [Google Scholar] [CrossRef]
  4. Grady, D.; Kipp, M. Continuum modelling of explosive fracture in oil shale. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1980, 17, 147–157. [Google Scholar] [CrossRef]
  5. Taylor, L.; Chen, E.; Kuszmaul, J. Microcrack-induced damage accumulation in brittle rock under dynamic loading. Comput. Methods Appl. Mech. Eng. 1986, 55, 301–320. [Google Scholar] [CrossRef]
  6. Kuszmaul, J. A new constitution model for fragmentation of rock under dynamic loading. In Proceedings of the 2nd International Symposium on Rock Fragmentation by Blasting, Keystone, CO, Canada, 23−26 August 1987; pp. 412–423. [Google Scholar]
  7. Xiong, Y.B.; Chen, J.J.; Hu, Y.L. Study on the key parameters of the Johnson-Holmquist constitutive model for concrete. Eng. Mech. 2012, 29, 121–127. [Google Scholar]
  8. Xiong, Y.B.; Peng, L.; Wang, L.Y. Numerical Simulation on Damage Effects of Underground Structure of Steel Pipeline-Concrete-Rock Under Explosive Loading. Mod. Appl. Phys. 2014, 5, 51–58. [Google Scholar]
  9. Du, C.; Song, S.; Zhang, J.P. Comparative study on three concrete constitutive models under blast loading. J. Ordnance Equip. Eng. 2022, 43, 49–56. [Google Scholar]
  10. Jin, Q.K. Dynamic Damage and Failure Model for Concrete Materials. J. Mil. Eng. 2006, 27, 10–13. [Google Scholar]
  11. Wu, H.J.; Huang, F.L.; Fu, Y.S. Numerical Simulation of Reinforced Concrete Breakage Under Internal Blast Loading. Trans. Beijing Inst. Technol. 2007, 27, 200–204. [Google Scholar]
  12. Qiang, H.F.; Sun, X.Y.; Wang, G. Numerical Simulation of Penetration in Concrete Sheet Based on SPH Method. Chin. J. High Press. Phys. 2019, 33, 67–76. [Google Scholar]
  13. Hu, Y.G.; Lu, W.B.; Chen, M. Comparison and improvement of blasting damage models for rock. Rock Soil Mech. 2012, 33, 3278–3284. [Google Scholar]
  14. Chen, J.H.; Zhang, J.S.; Li, X.P. Model of rock blasting-induced damage considering integrity of rock mass and its application. Chin. J. Geotech. Eng. 2016, 38, 857–866. [Google Scholar]
  15. Xie, F.J.; Zhang, J.S.; Chen, J.H. Dynamic damage model of rock under impact loads of compression and tension. J. Cent. South Univ. Sci. Technol. 2019, 50, 420–427. [Google Scholar]
  16. Yang, R.; Brwden, W.; Katsabanis, P. A new constitutive model for blast damage. Int. J. Rock Mech. Min. Sci. 1996, 33, 245–254. [Google Scholar] [CrossRef]
  17. Chen, G.; Jiang, Y.C.; He, P. Experimental study on static and dynamic tensile properties of granite. In Proceedings of the 6th National Symposium on Experimental Technology in Explosives Dynamics, Changsha, China, 23 July 2010; pp. 394–398. [Google Scholar]
  18. Dong, M.L. The Fundamentals of Nonlinear Mechanics of Concrete, 1st ed.; China Architecture & Building Press: Beijing, China, 1997; pp. 7–12. [Google Scholar]
  19. Gong, F.Q.; Lu, D.H.; Li, X.B. Experimental research of sandstone dynamic strength criterion under different strain rates. Rock Soil Mech. 2013, 34, 2433–2441. [Google Scholar]
  20. Wang, W.; Zhang, D.; Lu, F.Y. Experimental study on scaling the explosion resistance of a one-way square reinforced concrete slab under a close-in blast loading. Int. J. Impact Eng. 2012, 49, 158–164. [Google Scholar] [CrossRef]
Figure 1. Numerical calculation results of two models under different loading methods: (a) the evolution curve of hydrostatic pressure and damage under compress then tensile loading with volumetric strain (the tensile loading stage); (b) stress and damage evolution curves under uniaxial/triaxial tensile loading; (c) the curve illustrating the variation in tensile strength with strain rate.
Figure 1. Numerical calculation results of two models under different loading methods: (a) the evolution curve of hydrostatic pressure and damage under compress then tensile loading with volumetric strain (the tensile loading stage); (b) stress and damage evolution curves under uniaxial/triaxial tensile loading; (c) the curve illustrating the variation in tensile strength with strain rate.
Applsci 14 06852 g001
Figure 2. Results of damage in split specimen test: (a) damage cloud map of the HK model; (b) damage cloud map of the HK0 model; (c) typical experimental damage result.
Figure 2. Results of damage in split specimen test: (a) damage cloud map of the HK model; (b) damage cloud map of the HK0 model; (c) typical experimental damage result.
Applsci 14 06852 g002
Figure 3. The damage results of the thin plate after impact loading: (a) front surface; (b) back surface.
Figure 3. The damage results of the thin plate after impact loading: (a) front surface; (b) back surface.
Applsci 14 06852 g003
Table 1. Parameters of high-strength steel fiber concrete [7,10,11,12].
Table 1. Parameters of high-strength steel fiber concrete [7,10,11,12].
G/GPafC/MPaT/MPaCABnSmax
17.3131.211.40.0070.541.530.8210
ε f m i n D1D2PC/MPaUCPL/MPaULK1/GPa
0.020.041.092.40.0048000.185
K2/GPaK3/GPaKIC/(MN·cm−3/2)k/(cm−3)m
−1712082.747 × 10−43.0 × 10146
Table 2. Parameters of concrete [9,10,11,12].
Table 2. Parameters of concrete [9,10,11,12].
G/GPafC/MPaT/MPaCABnSmax
16.739.54.20.0070.791.60.617.0
ε f m i n D1D2PC/MPaUCPL/MPaULK1/GPa
0.010.0351.023.30.001586000.15817.4
K2/GPaK3/GPaKIC/(MN·cm−3/2)k/(cm−3)m
38.829.82.747 × 10−45.55 × 10166
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zou, X.; Xiong, Y.; Wang, L.; Zhou, Y.; Wang, W.; Zhong, F. A Dynamic Damage Constitutive Model of Rock-like Materials Based on Elastic Tensile Strain. Appl. Sci. 2024, 14, 6852. https://doi.org/10.3390/app14166852

AMA Style

Zou X, Xiong Y, Wang L, Zhou Y, Wang W, Zhong F. A Dynamic Damage Constitutive Model of Rock-like Materials Based on Elastic Tensile Strain. Applied Sciences. 2024; 14(16):6852. https://doi.org/10.3390/app14166852

Chicago/Turabian Style

Zou, Xuan, Yibo Xiong, Leiyuan Wang, You Zhou, Wanpeng Wang, and Fangping Zhong. 2024. "A Dynamic Damage Constitutive Model of Rock-like Materials Based on Elastic Tensile Strain" Applied Sciences 14, no. 16: 6852. https://doi.org/10.3390/app14166852

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop