Next Article in Journal
Temporary Wettability Tuning of PCL/PDMS Micro Pattern Using the Plasma Treatments
Next Article in Special Issue
Prediction of Corrosive Fatigue Life of Submarine Pipelines of API 5L X56 Steel Materials
Previous Article in Journal
Adaptivity in Bayesian Inverse Finite Element Problems: Learning and Simultaneous Control of Discretisation and Sampling Errors
Previous Article in Special Issue
Grain-Size Distribution Effects on the Attenuation of Laser-Generated Ultrasound in α-Titanium Alloy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Dynamic Mechanical Properties of Concrete under Uniaxial Compression

1
Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, Beijing 100124, China
2
Department of Civil Engineering, Faculty of Engineering, Fayoum University, 63514 Fayoum, Egypt
3
College of Water Resources and Civil Engineering, China Agricultural University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(4), 643; https://doi.org/10.3390/ma12040643
Submission received: 29 January 2019 / Revised: 14 February 2019 / Accepted: 15 February 2019 / Published: 20 February 2019
(This article belongs to the Special Issue Deformation, Fatigue and Fracture of Materials)

Abstract

:
Based on the base force element method (BFEM), the dynamic mechanical behavior of concrete under uniaxial compression loading at different strain rates is investigated. The concrete can be considered as a three-phase composite material composed of aggregate, cement mortar, and interfacial transition zone (ITZ) on the meso-level. A two-dimensional random aggregate model is generated by the Monte Carlo method. A multi-linear two-dimensional damage model is applied to describe the damage properties of each phase in the concrete. The strain-softening behavior, strain-rate effect, and failure patterns of the concrete are studied. The numerical results find that the peaks of compressive stress and compressive strain of concrete show the rate-sensitivity in various degrees under different strain rates. The calculated results of the dynamic enhancement factors are in a good agreement with the formula given by the Comité Euro-International du Béton (CEB) and other experimental results. The failure diagram of the specimen clearly describes the compressive failure process of the concrete specimen. This failure’s characteristics are similar to the experimental results.

1. Introduction

In practice, the concrete used in infrastructure is usually subjected to dynamic loading, including impact loading and sustained loading. The previous obtained results show that the mechanical properties and damage characteristics of concrete under dynamic loading are very different from those under static loading [1,2]. Therefore, it is of interest to investigate the dynamic behavior of concrete under dynamic loading at high strain rates.
Hitherto, much research based on the traditional experimental test has been conducted. Abram [3] was the first to find out the compressive strength of concrete and show the strain rate sensitivity under dynamic loading. Bischoff et al. [4] concluded that the strain rate plays a significant role in both the dynamic ultimate strength and the dynamic deformation behavior of plain concrete at high strain rates. Ross et al. [5] and John et al. [6] respectively carried out Split-Hopkinson pressure bar (SHPB) tests to study the dynamic behavior of concrete under tensile stress. Malvar et al. [7] and Williams [8] respectively conducted literature reviews to investigate the strain-rate effects on the dynamic increase factor (DIF) of strength and elastic moduli; thus, the empirical equations were given correspondingly. In addition, the mechanical properties of the interfacial transition zone (ITZ) and the effect of the ITZ on the properties of rock-concrete materials were investigated [9]. Erzar et al. [10] developed direct tensile tests, spalling tests, and edge-on impact tests at high strain rates. The influence of the ITZ between the aggregates and the cement mortar on the uniaxial dynamic tensile strength of concrete was also studied. The SHPB apparatus was used to study the mechanical properties of cement-based materials. Deformation and stress distribution in the specimens are non-uniform due to the composite microstructure of the materials. The strength and failure patterns were also studied in both quasi-static and dynamic loadings [11,12,13,14,15,16,17]. Li et al. [18] used both cylinder and cube concrete specimens to conduct compressive tests and the effects of the specimen shape and size on concrete strength subjected to different loading rates were investigated. Jurowski and Grzeszczyk [19] concluded that the stabilized dynamic elastic modulus of concrete is proportional to the initial static elastic modulus, and the coefficient of proportionality is affected by the type of aggregates. Using various inner diameters of specimens, Zhang et al. [20] investigated the inertial effect on the tensile strength of concrete materials under dynamic loading.
In recent years, numerical analysis with computer modeling has become an efficient method to study the properties of concrete. On the basis of a damage parameter, Simo and Ju [21,22] proposed a continuum isotropic model and an anisotropic elastoplastic-damage model, and the results are consistent with the existing experimental data. The mechanical properties of small-eccentric loaded reinforced concrete (RC) columns and the dynamic behavior of steel fiber-reinforced concrete (SFRC) beams under impact loading were investigated using Abaqus software [23,24]. Furthermore, finite element theory is one of the major methods used to analyze the mechanical properties of cement-based composite material [25,26,27,28]. In the meso-scale, concrete can be deemed as a three-phase composite which comprises aggregate and cement mortar, with an ITZ between both. Regarding concrete studies, several scholars around the world have emphasized the relationship between microstructure and macroscopic mechanical properties of concrete under dynamic loading with a numerical method.
Georgin et al. [29] exploited a viscous plastic model to simulate the SHPB test and studied the influence of inertial force and strain rates on concrete dynamic behavior. Snozzi et al. [30,31] and Gatuingt et al. [32,33] proposed a computational model to investigate the mechanical properties of concrete, composed of aggregate and mortar paste matrix, under dynamic loading of tension and compression. Also, Park et al. [34] analyzed the influence of impact loading at high strain rates on concrete-like materials using a dynamic finite element simulation. Zhou et al. [35,36] and Hao et al. [37,38] adopted numerical methods to analyze the influence of the ITZ on the dynamic failure patterns of the concrete, which is considered as a three-phased composite consisting of aggregate, mortar, and ITZ. Additionally, with the energy theory and the micro prestressed-solidification theory applied, Cusatis [39] utilized the developed Confinement Shear Lattice (CSL) model to analyze the effect of strain rates on the strength and failure behavior of concrete. However, Wu et al. [40] established a new experimental method of numerical simulation to identify the rate sensitive to the concrete dynamic tensile behavior. In this regard, in order to simulate the dynamic behavior of concrete under tensile, Zhou et al. [41] presented a two-dimensional meso-scale finite element model, validated by comparing with the experimental data from spall tests. Chen et al. [42] formulated a new dynamic compressive constitutive model applied to investigations of strain-rate effects and damage effects within the specimen.
To date, the relationship between the failure mechanism of concrete and the stress distribution at high strain rate is not clear; however, the mesoscopic components remarkably affect the macroscopic mechanical properties of plain concrete. Therefore, on the meso-level, it is important to know the dynamic behavior of the concrete.
In this paper, according to the investigations of concrete-like materials, under static loading with the base force element method (BFEM) proposed by Peng et al. [43,44,45], a dynamic base force element model is developed. In addition, the failure process of concrete under uniaxial dynamic loading at high strain rates is simulated. Strain-softening behavior and the failure process of modeled concrete under uniaxial compressive loading is also investigated. Finally, the effects of different loading velocities on compressive strength and the stress distribution in the specimens are studied.

2. Establishment of the Dynamic Base Force Element Equilibrium Equation

Generally, the behavior of real physical structures subjected to loadings or displacements is almost dynamic. According to the D’Alembert’s principle, considering the actions of inertial force and damping force, the balance equation of a single-degree-of-freedom system can be obtained as follows:
[ M ] { u ¨ ( t ) } + [ C ] { u ˙ ( t ) } + [ K d ] { u ( t ) } = { P ( t ) }
in which { u ¨ ( t ) } , { u ˙ ( t ) } , and { u ( t ) } respectively indicate the acceleration vector quantity, the velocity vector quantity, and the displacement vector quantity of each node in structure; [ M ] means global mass matrix; [ C ] means the damping matrix; [ K d ] is the global stiffness matrix based on the base force as proposed in the paper; and { P ( t ) } is the dynamic load array of the structure.
Furthermore, the full variable form of the Newmark-β method is adopted to solve the dynamic equilibrium equation in the study.

2.1. Base Force Element Stiffness Matrix

According to the BFEM, which is based on the potential energy principle, a plane triangular element matrix K I J [46], expressed by base force, can be obtained as follows (shown in Figure 1):
K I J = E 2 A ( 1 + ν ) [ 2 ν 1 2 ν m I m J + m I J U + m J m I ]
where E , ν , and A mean the Young’s modulus, the Poisson’s ratio, and the area of an element, respectively; U is the unit tensor; u I , u J , u K represent the displacements of the triangular element vertex.
For plane strain problems, x and y represent Cartesian coordinate system, and the element matrix K I J can be expressed as follows:
[ K I J ] e = E 2 A ( 1 + ν ) [ 2 2 ν 1 2 ν m x I m x J + m y I m y J 2 ν 1 2 ν m x I m y J + m y I m x J 2 ν 1 2 ν m y I m x J + m x I m y J 2 2 ν 1 2 ν m y I m y J + m x I m x J ]
where m I J = m I m J , m I and m J can be described as follows and are shown in Figure 2:
m I = m i I e i = 1 2 ( L I J n I J + L K I n K I ) = 1 2 ( L I J n i I J e i + L K I n i K I e i ) = 1 2 ( L I J n i I J + L K I n i K I ) e i ,
m J = m i J e i = 1 2 ( L J K n J K + L I J n I J ) = 1 2 ( L J K n i J K e i + L I J n i I J e i ) = 1 2 ( L J K n i J K + L I J n i I J ) e i ,
{ m x I m y I } = 1 2 ( L I J { n x I J n y I J } + L L I { n x K I n y K I } ) ,
{ m x J m y J } = 1 2 ( L J K { n x J K n y J K } + L I J { n x I J n y I J } ) .
in which I , J , K represent the vertexes of the triangular element; L I J , L K I , L J K are the lengths of the element boundary lines; n I J , n K I , n J K are the normal vectors of the element boundary lines. Average strain components, which could take the place of the real strain in the case of small deformation, are given as follows:
ε ¯ x = 1 A I = 1 n ( u I x m x I ) ,
ε ¯ y = 1 A I = 1 n ( u I y m y I ) ,
γ ¯ x y = 1 A I = 1 n ( u I x m y I + u I y m x I ) .

2.2. Damping Matrix and Mass Matrix

In the finite element dynamic analysis, lumped mass matrix, one of the most common forms of the element mass matrix, is applied to simplify the calculation and reduce the storage space. According to the equivalent static principle, the mass of the element can be evenly distributed among three vertexes of a triangular element, according to the hypothesis that there is no interaction between the inertial forces of each vertex.
In dynamic response problems it is generally assumed that the viscous damping force is proportional to the velocity of particle motion.
Therefore, the effect of damping in the dynamic structural analysis at a high loading rate cannot be ignored. As a widely used orthogonal damping model, the Rayleigh damping model can be expressed as follows:
[ C ] = α [ M ] + β [ K d ] ,
α = 2 ζ ω 1 ω 2 ω 1 + ω 2 ,
β = 2 ζ ω 1 + ω 2 ,
[ M ] e = ρ b T A 3 g [ 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 ] ,
in which α and β are the mean coefficients of proportionality, which can be calculated from Equations (6) and (7); ζ means the damping ratio; ω 1 and ω 2 are the first and the second order angular frequencies of the specimen.
It should be noted that, in the process of the meso-mechanical failure of concrete, there are few research works on damping models. The damping theory is a methodology which aims to describe the microscopic mechanism of damping in a macroscopic way. For the general structural analysis, Rayleigh damping is usually used to approximately describe the damping characteristics. Normally, the damping ratio of the engineering structure ranges from 0.01 to 0.1. The average value 0.05 is applied in this paper.
In this paper, to simplify the model, the influence of damage development on damping matrices is not considered. Consequently, the initial stiffness matrix of the specimen is adopted in the present dynamic analysis and the damping ratio is regarded as a constant. It is assumed that the mass matrix is independent from the damage state due to the conservation of mass.

3. Meso-Structure of Concrete

It is well known that concrete is a multi-phase heterogeneous brittle material. The mechanical properties of concrete are determined by components in its mixture. In this paper, concrete is described as a three-phase composite material at the meso-level, composed of coarse aggregates, cement mortar, and an ITZ between both phases. With the spherical aggregate applied, the microstructure in concrete can be depicted as in Figure 3. In the model only coarse aggregates with a particle size greater than 5 mm are represented clearly, while the other smaller aggregates are mixed up in the cement paste matrix.

3.1. Random Aggregate Model

In this paper, in order to simulate the real distribution of aggregate as possible, the Monte Carlo method is used to establish the two-dimensional random aggregate model. Based on the Fuller three-dimensional aggregate gradation curve, Walraven [47] developed the two-dimensional cross-section aggregate gradation curve. The cumulative probability of aggregate particles satisfied the condition of D < D 0 (herein D means the diameter of the aggregate) is calculated using the following formula:
P c ( D D 0 ) = P k [ 1.065 ( D 0 D max ) 1 / 2 0.053 ( D 0 D max ) 4 0.012 ( D 0 D max ) 6 0.0045 ( D 0 D max ) 8 + 0.0025 ( D 0 D max ) 10 ]
in which P k represents the percentage of the volume of aggregate in the whole concrete specimen; P c represents the cumulative probability of aggregate particles with a size smaller than D 0 ; D max is the maximum diameter of the aggregate particles.
The number of aggregate particles meeting the requirement of D 1 < D < D 2 can be calculated by the following equation:
n = [ P c ( D < D 2 ) P c ( D < D 1 ) ] × A / A i
where A is the cross-section area of the concrete specimen and A i is the area of the representative particle diameter of the aggregate.
According to the theory of Fuller’s maximum density curve, three representative diameters (10 mm, 20 mm, and 32.5 mm) of aggregates are selected to calculate number of aggregate particles. The numbers are, respectively, 66, 11, and 3, as obtained from Equation (20). The aggregates calculated are put into the two-dimensional region, the same size as concrete specimen, by using the Monte Carlo method with three sets of different random numbers. Each aggregate should be in line with the boundary conditions and not overlap with the other aggregates before it is put on. The modeled aggregates are shown in Figure 4.

3.2. Mesh Generation Method and Element Attributes

In the finite element analysis, two common methods are adopted to deal with the thickness of ITZ in order to simplify the calculation. One considers the thickness of the ITZ from 0.5 mm to 2 mm [48,49], and the other does not consider the thickness of the ITZ [50,51,52,53]. In this paper, the dynamic behavior of concrete is studied by using the former method and the size of the mesh element is deemed as the ITZ thickness to simplify the mesoscopic model.
The finite element mesh is shadowed on the random aggregate model. The different mechanical properties are assigned to the corresponding elements. A linear elastic triangular finite element grid is applied in the paper, shown in Figure 5. The type of element is determined by the position of element nodes. The element could be deemed as an aggregate (or cement mortar) element when the three nodes of the element are all located at the aggregate (or cement mortar) region; otherwise, the element is the ITZ when the element is located in the both the aggregate and cement mortar.

4. Dynamic Behavior for Concrete Meso-Components

4.1. Concrete Dynamic Damaged Model

As known, the mechanical properties of microscopic components have a great influence on the fracture damage behavior of concrete. Therefore, the heterogeneity of concrete should be taken into account in this numerical simulation. In this paper, concrete can be treated as a three-phase composite material composed of coarse aggregate, cement mortar, and the ITZ. Each phase material in concrete is assumed to be homogeneous and isotropic. The constitutive relation is presented in Figure 6. The maximum principal stress criterion is applied as the failure criteria in the study. The reduction of the elastic modulus E of the material in the damage processes can be expressed as follows:
E = E 0 ( 1 D )
in which E 0 is the initial modulus of elasticity and D is the damage factor, defined as follows:
D t = { 0 ε max ε t 0 1 ε t 0 ε max + ε max ε t 0 η t ε t 0 ε t 0 ε t 0 ε max ( 1 α ) ε t 0 < ε max η t ε t 0 1 α ξ t η t ε max η t ε t 0 ε max + α ε t 0 ε max η t ε t 0 < ε max ξ t ε t 0 1 ε max > ξ t ε t 0
D c = { 1 β γ ε max λ ε c 0 1 1 β 1 λ ε max λ ε c o ε max β ε c o ε max λ ε c 0 < ε max ε c 0 1 1 γ 1 η c ε max ε c o ε max ε c o ε max ε c 0 < ε max η c ε c 0 1 γ ε c 0 ε max η c ε c 0 < ε max ξ c ε c 0 1 ε max > ξ c ε c 0
in which ε 0 means the principle strain; η denotes the residual strain coefficient; ξ represents the ultimate strain coefficient.
In Figure 6, f c and f t stand for compressive and tensile strength; furthermore, the subscripts t and c respectively symbolize tensile and compressive features. The material parameters are given in Table 1, which is obtained after calculations.

4.2. Dynamic Increase Factor (DIF) for Concrete

In practice, damage patterns and mechanical properties of concrete under dynamic loading present different forms, which is called the strain-rate effect and is characterized by the dynamic increase factor (DIF). The dynamic increase factor for the compressive strength is recommended by Comité Euro-International du Béton [54] as the following:
D I F c = f c d / f c = { ( ε ˙ ε ˙ s ) 1.026 α ε ˙ 30 s 1 γ ( ε ˙ ε ˙ s ) 1 / 3 ε ˙ > 30 s 1
in which f c d is the dynamic uniaxial compressive strength and f c is the quasi-static uniaxial compressive strength; ε ˙ is the quasi-dynamic strain rate and ε ˙ s is the quasi-static strain rate; f c s is the quasi-static uniaxial compression strength and f c 0 = 10   MPa . It should be noted that α = 1 / ( 5 + 9 f c s / f c 0 ) and lg γ = 6.156 α 2 . As known from the empirical formula, the value of D I F c is 1 when the strain rate ε ˙ is 30 × 10 6 s 1 , which is called the quasi-static load modal. In this paper, the minimum strain rate 10 3 s 1 is deemed as the quasi-static strain rate for comparison. According to the previous study, the strain-rate sensitivity of concrete characteristics, such as the Poisson’s ratio, elasticity modulus, energy dissipation capacity, and so forth, are much lower than the tensile and compressive strength of concrete [55].

5. Numerical Examples and Results

5.1. Boundary Conditions and Loading Model

In this section, the dynamic test of concrete subjected to uniaxial compressive loads is simulated with the BFEM. A standard concrete specimen, with the dimensions 150 mm × 150 mm, is chosen to conduct the test. The particle diameter of the coarse aggregates ranges from 5 mm to 40 mm. The boundary conditions are shown in Figure 7. The bottom and top surfaces of the specimen are restricted only in the y-direction; other surfaces of the specimen are free in all directions. The influence of friction between the specimen and the loaded end on the compressive strength of concrete is ignored.
Continuous and uniform vertical displacement loading is used on the specimens. The direction of vertical displacement loading is parallel to the y-axis. Dynamic displacement step-load is applied in this study and the velocity of displacement loading is well controlled by the duration of each load-step. Vertical strain rate ε ˙ m under different loading velocity can be calculated as the following:
ε ˙ m = 2 v h
where v is the constant loading velocity and h is the height of the specimen. The loading curves of concrete specimens are shown in Figure 8.

5.2. Dynamic Failure Behavior of Concrete under Uniaxial Compressive Stress

Based on the meso-mechanical numerical simulation proposed in the present study, the specimen can be subdivided into triangular finite element meshes where the mesh size is 1 mm. With eight groups with different loading rates applied, dynamic uniaxial compression tests are carried out on the concrete specimens. The corresponding macroscopic nominal strain rates are 10 3 / s , 10 2 / s , 10 1 / s , 1 / s , 10 / s , 30 / s , 80 / s , and 100 / s .
The stress–strain curves of concrete specimens are shown in Figure 9. The peak values of strain and stress are shown in Figure 10 at the different strain-rates. The damage process of concrete specimens under different strain-rate compression levels is shown in Figure 11. The failure patterns at different rates are illustrated in Figure 12. The distribution of the maximum principal stress in the damaged concrete is shown in Figure 13. It is to be noted that ε ˙ means the strain rate, ε is the strain, and σ is the compressive stress. In Figure 13, the positive numbers represent the elements under tensile stress and the negative numbers represent the elements are under compressive stress.
The stress–strain curves are similar to the static at low strain rates from 0.01 / s to 0.1 / s and the compressive strengths present a small enhancement. However, at the high strain rates from 1 / s to 100 / s , the curves slow down gradually as the strain rate increases and the stress peak and strain peak increase significantly, and especially the ratio of dynamic compressive strength at the rate of 100 / s to the static compressive strength is about 2.67. It is also found that the increasing trend of the strain is more prominent than that of the stress from Figure 10.
The meso-cracks first occur and spread in the ITZ, or the weak link, in cement mortar. Finally, a region is formed, filled with a concentrative zonal crack, which results in the destruction of concrete. At the high strain rates from 10 / s to 100 / s , the number of cracks greatly increases and the cracks are diffused through the whole concrete specimen—a few cracks even pass through the aggregate area. Coalescent cracks as a diffusion state can be seen in Figure 12. In addition, from Figure 9 and Figure 12, there is also a minor effect of the distribution of aggregates on the compressive stress. However, the distribution of aggregates can significantly affect the initiation and propagation of cracks.
In Figure 14, it can be observed that the curve of the present result goes up precipitously at the strain rate of 0.1 / s . The results obtained show good agreement with the Comité Euro-International du Béton (CEB) standard and other experimental data at the strain rates from 10 / s to 100 / s . At the strain rates from 0.01 / s to 10 / s , the DIF is lower than the CEB standard.

6. Discussion

The BFEM for dynamic analysis of concrete introduced in the present study verifies an efficient numerical simulation method to investigate the damage mechanism of concrete. The strain-rate effect and the failure characteristics are explored using two-dimensional models. Several conclusions can be obtained as follows:
(1) The present results prove that the failure process of concrete under dynamic compression is simulated well. With the increase of stress, the material enters the nonlinear stage and the stress–strain curves show a nonlinear increase relationship, simultaneously.
(2) The variation tendency of the DIF of concrete at different strain rates is consistent with the available experimental test data and the CEB empirical formula.
(3) At high strain rates, cracks increase and, in a diffusion state, some of the elements present fracture damage and more energy is released, which could enhance the dynamic strength of the concrete.
(4) The strain peaks of concrete present rate-sensitivity under different strain rates similar to stress peaks.
For future work, a three-dimensional model of dynamic problems could be addressed by the BFEM and other conditions could be considered; for instance, considering the shape of aggregates, complicated boundary constraints, and so on, so as to simulate real concrete as much as possible.

Author Contributions

Y.P. made the revision and the final approval of the article. Q.W. carried out the numerical experiments, interpretation of results and article. L.Y. provided the technical support of the computer software program and the comprehensive review. M.M.A.K. made the check and revision of language and grammar of the article. H.P. carried out the data collection and the data processing.

Funding

This research received no external funding.

Acknowledgments

This work is supported by National Science Foundation of China (10972015, 11172015), the Beijing Natural Science Foundation (8162008), and the Pre-exploration Project of Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology (USDE201404).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bischoff, P.H.; Perry, S.H. Compressive behavior of concrete at high strain rates. Mater. Struct. 1991, 24, 425–450. [Google Scholar] [CrossRef]
  2. Evans, R.H. Effect of rate of loading on the mechanical properties of some materials. J. Ice 1942, 18, 296–306. [Google Scholar] [CrossRef]
  3. Abrams, D.A. Effect of rate of application of load on the compressive strength of concrete. ASTM J. 1917, 17, 70–78. [Google Scholar]
  4. Bischoff, P.H.; Perry, S.H. Impact behavior of plain concrete loaded in uniaxial compression. J. Eng. Mech. 1995, 25, 685–693. [Google Scholar] [CrossRef]
  5. Ross, C.A.; Tedesco, J.W.; Kuennen, S.T. Effects of strain rate on concrete strength. Aci Mater. J. 1995, 92, 37–47. [Google Scholar]
  6. John, R.; Antoun, T.; Rajendran, A.M. Effect of strain rate and size on tensile strength of concrete. Shock Compression Condensed Matter 1992, 1994, 501–504. [Google Scholar]
  7. Malvar, L.J.; Crawford, J.E. Dynamic Increase Factors for Concrete. In Proceedings of the Twenty-Eighth DDESB Seminar, Orlando, FL, USA, 18–20 August 1998; pp. 1–17. [Google Scholar]
  8. Williams, M.S. Modeling of local impact effects on plain and reinforced concrete. Aci Struc. J. 1994, 91, 178–187. [Google Scholar]
  9. Kishen, J.M.C. Fracture of Rock-Concrete Interfaces: Laboratory Tests and Applications. ACI Struct. J. 2004, 101, 325–331. [Google Scholar]
  10. Erzar, B.; Forquin, P. Experiments and mesoscopic modelling of dynamic testing of concrete. Mech. Mater 2011, 43, 505–527. [Google Scholar] [CrossRef]
  11. Grote, D.L.; Park, S.W.; Zhou, M. Experimental characterization of the dynamic failure behavior of mortar under impact loading. J. Appl. Phys. 2001, 89, 2115–2123. [Google Scholar] [CrossRef]
  12. Grote, D.L.; Park, S.W.; Zhou, M. Dynamic behavior of concrete at high strain rates and pressures: I. experimental characterization. Int. J. Impact Eng. 2001, 25, 869–886. [Google Scholar] [CrossRef]
  13. Jin, X.C.; Hou, C.; Fan, X. Quasi-static and dynamic experimental studies on the tensile strength and failure pattern of concrete and mortar discs. Sci Rep. 2017, 7, 15305. [Google Scholar] [CrossRef] [PubMed]
  14. Guo, Y.B.; Gao, G.F.; Jing, L. Response of high-strength concrete to dynamic compressive loading. Int. J. Impact Eng. 2017, 108, 114–135. [Google Scholar] [CrossRef]
  15. Levi-Hevroni, D.; Kochavi, E.; Kofman, B. Experimental and numerical investigation on the dynamic increase factor of tensile strength in concrete. Int. J. Impact Eng. 2018, 114, 93–104. [Google Scholar] [CrossRef]
  16. Zhou, J.; Chen, X.; Wu, L. Influence of free water content on the compressive mechanical behavior of cement mortar under high strain rate. Sadhana 2011, 36, 357–369. [Google Scholar] [CrossRef]
  17. Quast, M.; Curbach, M. Concrete under biaxial dynamic compressive loading. Procedia Eng. 2017, 210, 24–31. [Google Scholar] [CrossRef]
  18. Li, M.; Hao, H.; Shi, Y. Specimen shape and size effects on the concrete compressive strength under static and dynamic tests. Constr. Build. Mater. 2018, 161, 84–93. [Google Scholar] [CrossRef]
  19. Jurowski, K.; Grzeszczyk, S. Influence of Selected Factors on the Relationship between the Dynamic Elastic Modulus and Compressive Strength of Concrete. Materials 2018, 11, 477. [Google Scholar] [CrossRef]
  20. Zhang, S.; Lu, Y.; Jiang, X. Inertial effect on concrete-like materials under dynamic direct tension. Int. J. Prot. Struct. 2018, 9, 377–396. [Google Scholar] [CrossRef]
  21. Simo, J.C.; Ju, J.W. Strain- and stress-based continuum damage model—II. Comput aspects. Int. J. Solids Struct. 1989, 23, 841–869. [Google Scholar] [CrossRef]
  22. Simo, J.C.; Ju, J.W. Strain- and stress-based continuum damage models—I. Formulation. Int. J. Solids Struct. 1987, 23, 821–840. [Google Scholar] [CrossRef]
  23. Jin, L.; Zhang, R.; Dou, G.; Xu, J.; Du, X. Experimental and numerical study of reinforced concrete beams with steel fibers subjected to impact loading. Int. J. Damage Mech. 2018, 27, 1058–1083. [Google Scholar] [CrossRef]
  24. Jin, L.; Ding, Z.; Li, D.; Du, X. Experimental and numerical investigations on the size effect of moderate high-strength reinforced concrete columns under small-eccentric compression. Int. J. Damage Mech. 2018, 27, 657–685. [Google Scholar] [CrossRef]
  25. Li, W.S.; Wu, J.Y. A consistent and efficient localized damage model for concrete. Int. J. Damage Mech. 2018, 27, 541–567. [Google Scholar] [CrossRef]
  26. Wardeh, M.A.; Toutanji, H.A. Parameter estimation of an anisotropic damage model for concrete using genetic algorithms. Int. J. Damage Mech. 2017, 26, 801–825. [Google Scholar] [CrossRef]
  27. Yoo, D.Y.; Kang, S.T.; Banthia, N.; Yoon, Y.S. Nonlinear finite element analysis of ultra-high-performance fiber-reinforced concrete beams. Int. J. Damage Mech. 2017, 26, 735–757. [Google Scholar] [CrossRef]
  28. Liang, J.; Ren, X.; Li, J. A competitive mechanism driven damage-plasticity model for fatigue behavior of concrete. Int. J. Damage Mech. 2016, 25, 377–399. [Google Scholar] [CrossRef]
  29. Georgin, J.F.; Reynouard, J.M. Modeling of structures subjected to impact: concrete behavior under high strain rate. Cem. Concr. Compos. 2003, 25, 131–143. [Google Scholar] [CrossRef]
  30. Snozzi, L.; Caballero, A.; Molinari, J.F. Influence of the meso-structure in dynamic fracture simulation of concrete under tensile loading. Cem. Concr. Res. 2011, 41, 1130–1142. [Google Scholar] [CrossRef]
  31. Snozzi, L.; Gatuingt, F.; Molinari, J.F. A meso-mechanical model for concrete under dynamic tensile and compressive loading. Int. J. Fract. 2012, 178, 179–194. [Google Scholar] [CrossRef]
  32. Gatuingt, F.; Snozzi, L.; Molinari, J.F. Numerical determination of the tensile response and the dissipated fracture energy of concrete: role of the mesostructure and influence of the loading rate. Int. J. Numer. Anal. Methods Geomech. 2013, 37, 3112–3130. [Google Scholar] [CrossRef]
  33. Gatuingt, F.; Snozzi, L.; Molinari, J.F. Determination of the Dynamic Tensile Response and Dissipated Fracture Energy of Concrete with a Cohesive Element Model. International Conference on Numerical Modeling Strategies for Sustainable Concrete Structures—SSCS 2012, June 2012, Aix en Provence, France. 2017. Available online: https://hal.archives-ouvertes.fr/hal-01624529 (accessed on 15 January 2019).
  34. Park, S.W.; Xia, Q.; Zhou, M. Dynamic behavior of concrete at high strain rates and pressures: II. numerical simulation. Int. J. Impact Eng. 2001, 25, 887–910. [Google Scholar] [CrossRef]
  35. Zhou, X.Q.; Hao, H. Mesoscale modelling of concrete tensile failure mechanism at high strain rates. Comput. Struct. 2008, 86, 2013–2026. [Google Scholar] [CrossRef]
  36. Zhou, X.Q.; Hao, H. Modelling of compressive behavior of concrete-like materials at high strain rate. Int. J. Solids Struct. 2008, 45, 4648–4661. [Google Scholar] [CrossRef]
  37. Hao, H.; Hao, Y.; Li, Z.X. A numerical study of factors influencing high-speed impact tests of concrete material properties. In Proceedings of the 8th International Conference on Shock & Impact Loads on Structures, Adelaide, Australia, 2–4 December 2009; pp. 37–52. [Google Scholar]
  38. Hao, H.; Zhou, X.Q. Concrete material model for high rate dynamic analysis. In Proceedings of the Seventh International Conference on shock & impact loads on structures, Beijing, China, 17–19 October 2007; pp. 753–768. [Google Scholar]
  39. Cusatis, G. Strain-rate effects on concrete behavior. Int. J. Impact Eng. 2011, 38, 162–170. [Google Scholar] [CrossRef]
  40. Wu, H.; Zhang, Q.; Huang, F. Experimental and numerical investigation on the dynamic tensile strength of concrete. Int. J. Impact Eng. 2005, 32, 605–617. [Google Scholar] [CrossRef]
  41. Zhou, W.; Tang, L.; Liu, X. Mesoscopic simulation of the dynamic tensile behavior of concrete based on a rate-dependent cohesive model. Int. J. Impact Eng. 2016, 95, 165–175. [Google Scholar] [CrossRef]
  42. Chen, X.; Wu, S.; Zhou, J. Experimental and modeling study of dynamic mechanical properties of cement paste, mortar and concrete. Constr. Build. Mater. 2013, 47, 419–430. [Google Scholar] [CrossRef]
  43. Peng, Y.J.; Liu, Y.H.; Pu, J.W.; Zhang, L. Application of base force element method to mesomechanics analysis for recycled aggregate concrete. Math Probl. Eng. 2013, 2013, 1–8. [Google Scholar] [CrossRef]
  44. Peng, Y.J.; Pu, J.W. Micromechanical investigation on size effect of tensile strength for recycled aggregate concrete using BFEM. Int. J. Mech. Mater. Des. 2016, 12, 525–538. [Google Scholar] [CrossRef]
  45. Peng, Y.J.; Chu, H.; Pu, J.W. Numerical simulation of recycled concrete using convex aggregate model and base force element method. Adv. Mater. Sci. Eng. 2016, 2016, 1–10. [Google Scholar] [CrossRef]
  46. Peng, Y.J.; Dong, Z.L.; Peng, B.; Liu, Y.H. Base force element method (BFEM) on potential energy principle for elasticity problems. Int. J. Mech. Mater. Des. 2011, 7, 245–251. [Google Scholar] [CrossRef]
  47. Walraven, J.C.; Reinhardt, H.W. Theory and experiments on the mechanical behavior of cracks in plain and reinforced concrete subjected to shear loading. Heron. 1991, 26, 26–35. [Google Scholar]
  48. Du, X.L.; Jin, L.; Ma, G. A meso-scale analysis method for the simulation of nonlinear damage and failure behavior of reinforced concrete members. Int. J. Damage Mech. 2013, 22, 878–904. [Google Scholar] [CrossRef]
  49. Jayasuriya, A.; Adams, M.P.; Bandelt, M.J. Understanding variability in recycled aggregate concrete mechanical properties through numerical simulation and statistical evaluation. Constr. Build. Mater. 2018, 178, 301–312. [Google Scholar] [CrossRef]
  50. Nilsen, A.U.; Monteiro, P.J.M. Concrete: A three phase material. Cem. Concr. Res. 1993, 23, 147–151. [Google Scholar] [CrossRef]
  51. Agioutantis, Z.G.; Stiakakis, C.; Kleftakis, S. Numerical simulation of the mechanical behavior of epoxy based mortars under compressive loads. Comput. Struct. 2002, 80, 2071–2084. [Google Scholar] [CrossRef]
  52. Fakhari Tehrani, F.; Absi, J.; Allou, F. Heterogeneous numerical modeling of asphalt concrete through use of a biphasic approach: Porous matrix/inclusions. Comput. Mater. Sci. 2013, 69, 186–196. [Google Scholar] [CrossRef]
  53. Du, X.L.; Jin, L.; Ma, G. Numerical simulation of dynamic tensile-failure of concrete at meso-scale. Int. J. Impact Eng. 2014, 66, 5–17. [Google Scholar] [CrossRef]
  54. Comité Euro-International du Béton. CEB-FIP Model Code 1990; Redwood Books: Trowbridge, Wiltshire, UK, 1993. [Google Scholar]
  55. Dilger, W.H.; Koch, R.; Kowalczyk, R. Ductility of plain and confined concrete under different strain rates. ACI J. 1984, 81, 73–81. [Google Scholar]
Figure 1. A triangular element.
Figure 1. A triangular element.
Materials 12 00643 g001
Figure 2. Construction of a triangular element stiffness.
Figure 2. Construction of a triangular element stiffness.
Materials 12 00643 g002
Figure 3. Three-phase sphere model of concrete.
Figure 3. Three-phase sphere model of concrete.
Materials 12 00643 g003
Figure 4. Random aggregate models of concrete. (a) Specimen 1; (b) specimen 2; (c) specimen 3.
Figure 4. Random aggregate models of concrete. (a) Specimen 1; (b) specimen 2; (c) specimen 3.
Materials 12 00643 g004
Figure 5. The diagram of mesh generation.
Figure 5. The diagram of mesh generation.
Materials 12 00643 g005
Figure 6. Mechanical constitutive model of materials.
Figure 6. Mechanical constitutive model of materials.
Materials 12 00643 g006
Figure 7. Loading model of compression test.
Figure 7. Loading model of compression test.
Materials 12 00643 g007
Figure 8. Loading curve of dynamic uniaxial compression test.
Figure 8. Loading curve of dynamic uniaxial compression test.
Materials 12 00643 g008
Figure 9. Stress–strain curve of concrete under dynamic uniaxial compressive stress at different strain rates.
Figure 9. Stress–strain curve of concrete under dynamic uniaxial compressive stress at different strain rates.
Materials 12 00643 g009
Figure 10. The comparison of the peak values of strain and stress.
Figure 10. The comparison of the peak values of strain and stress.
Materials 12 00643 g010
Figure 11. The damage process of concrete specimen 1 under dynamic uniaxial compressive stress; (a) ε ˙ = 0.001   s 1 ; (b) ε ˙ = 1   s 1 ; (c) ε ˙ = 10   s 1 ; (d) ε ˙ = 100   s 1 .
Figure 11. The damage process of concrete specimen 1 under dynamic uniaxial compressive stress; (a) ε ˙ = 0.001   s 1 ; (b) ε ˙ = 1   s 1 ; (c) ε ˙ = 10   s 1 ; (d) ε ˙ = 100   s 1 .
Materials 12 00643 g011aMaterials 12 00643 g011b
Figure 12. The failure pattern of concrete specimens under dynamic uniaxial compressive stress; (a) specimen 1; (b) specimen 2; (c) specimen 3.
Figure 12. The failure pattern of concrete specimens under dynamic uniaxial compressive stress; (a) specimen 1; (b) specimen 2; (c) specimen 3.
Materials 12 00643 g012
Figure 13. The distribution of maximum principal stress at different strain rates of specimen 1.
Figure 13. The distribution of maximum principal stress at different strain rates of specimen 1.
Materials 12 00643 g013
Figure 14. The effect of strain rate on compressive strength of plain concrete.
Figure 14. The effect of strain rate on compressive strength of plain concrete.
Materials 12 00643 g014
Table 1. Mechanical parameters of materials.
Table 1. Mechanical parameters of materials.
Mechanical ParametersCement MortarInterfacial Transition Zone (ITZ)Aggregate
Density ρ (kg/m3)210017002700
Poisson’s ratio ν 0.220.20.16
Strength (tensile/compressive) σ (MPa)3.2/322.5/257/70
λ 0.250.250.80
β 0.850.650.90
γ 0.350.350.35
α 0.30.30.3
η t / η c 4/43/35/5
ξ t / ξ c 10/1010/1010/10

Share and Cite

MDPI and ACS Style

Peng, Y.; Wang, Q.; Ying, L.; Kamel, M.M.A.; Peng, H. Numerical Simulation of Dynamic Mechanical Properties of Concrete under Uniaxial Compression. Materials 2019, 12, 643. https://doi.org/10.3390/ma12040643

AMA Style

Peng Y, Wang Q, Ying L, Kamel MMA, Peng H. Numerical Simulation of Dynamic Mechanical Properties of Concrete under Uniaxial Compression. Materials. 2019; 12(4):643. https://doi.org/10.3390/ma12040643

Chicago/Turabian Style

Peng, Yijiang, Qing Wang, Liping Ying, Mahmoud M. A. Kamel, and Hongtao Peng. 2019. "Numerical Simulation of Dynamic Mechanical Properties of Concrete under Uniaxial Compression" Materials 12, no. 4: 643. https://doi.org/10.3390/ma12040643

APA Style

Peng, Y., Wang, Q., Ying, L., Kamel, M. M. A., & Peng, H. (2019). Numerical Simulation of Dynamic Mechanical Properties of Concrete under Uniaxial Compression. Materials, 12(4), 643. https://doi.org/10.3390/ma12040643

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