Next Article in Journal
Analysis of Behaviors of the Railway Subgrade with a New Waterproof Seal Layer
Next Article in Special Issue
Compressive Fatigue Investigation on High-Strength and Ultra-High-Strength Concrete within the SPP 2020
Previous Article in Journal
Selective Electron Beam Melting (SEBM) of Pure Tungsten: Metallurgical Defects, Microstructure, Texture and Mechanical Properties
Previous Article in Special Issue
Influence of Moisture Content and Wet Environment on the Fatigue Behaviour of High-Strength Concrete
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Damage Evolution of Steel Fibre-Reinforced High-Performance Concrete in Low-Cycle Flexural Fatigue: Numerical Modeling and Experimental Validation

1
Institute for Structural Engineering, Faculty of Architecture and Civil Engineering, Bergische Universität Wuppertal, Pauluskirchstraße 11, 42285 Wuppertal, Germany
2
Institute of Mechanics, Department Civil Engineering, Faculty of Engineering, University of Duisburg-Essen, Universitätsstraße 15, 45141 Essen, Germany
*
Authors to whom correspondence should be addressed.
Materials 2022, 15(3), 1179; https://doi.org/10.3390/ma15031179
Submission received: 12 November 2021 / Revised: 14 January 2022 / Accepted: 21 January 2022 / Published: 3 February 2022
(This article belongs to the Special Issue Cyclic Deterioration of Concrete)

Abstract

:
This contribution aims to analyze the deterioration behaviour of steel fibre-reinforced high-performance concrete (HPC) in both experiments as well as numerical simulations. For this purpose, flexural tensile tests are carried out on beams with different fibre contents and suitable damage indicators are established to describe and calibrate the damage behaviour numerically using a phase-field model approach. In addition to conventional measurement methods, the tests are equipped with acoustic emission sensors in order to obtain a more precise picture of crack evolution by observing acoustic events. It is shown that, in addition to classical damage indicators, such as stiffness degradation and absorbed energy, various acoustic indicators, such as the acoustic energy of individual crack events, can also provide information about the damage progress. For the efficient numerical analysis of the overall material behaviour of fibre-reinforced HPC, a phenomenological material model is developed. The data obtained in the experiments are used to calibrate and validate the numerical model for the simulation of three-point bending beam tests. To verify the efficiency of the presented numerical model, the numerical results are compared with the experimental data, e.g., load-CMOD curves and the degradation of residual stiffness.

1. Introduction

Steel fibre-reinforced concrete (SFRC), high-performance concrete (HPC) with steel fibres, as well as ultra-high-performance concrete (UHPC) with fibres have been objects of research worldwide and have been widely applied throughout the last decades [1,2]. The advantages of steel fibres as reinforcement in concrete include increasing the load-bearing capacity of cracked concrete, reducing crack width to increase durability, and increasing energy absorption capacity and ductility under cyclic loading. Therefore, steel fibre-reinforced concrete has been applied in several areas. For example, in precast segments for tunnel linings, e.g., for the Brightwater tunnels or the Crossrail project in London [3,4,5]. Meanwhile, a state-of-the-art report for using steel fibres in precast tunnel segments has been published [6]. Another application of SFRC are joints in reinforced beam–column structures, which may be subjected to earthquakes [7,8,9,10]. In this case, steel fibres are added to improve the energy dissipation of the joints.
The load-bearing capacity and crack width-reduction capability of SFRC is also used for slabs on ground, such as industrial floorings, in transportation infrastructure [11,12], or in applications for field-cast ultra-high-performance concrete joints [13,14] as well as SFRC for concrete bridge decks or bridge deck overlays [13,15,16,17,18]. In the latter applications, additional cyclic loadings have to be considered. High-cycle loadings exceeding 10 7 load cycles are expected, for example, for the support structures of offshore wind energy turbines or high-performance materials for grouting joints of precast segments of wind towers. Here, especially ultra-high-performance concrete with and without fibres is often used in the so called grouted connections [19,20].
In addition, the increasing compressive strength of HPC and UHPC is expected to decrease the dead weight of concrete structures, increasing their susceptibility to dynamic loads. Therefore, the performance of SFRC in terms of cyclic behaviour subjected to very high numbers of load cycles has to be closely investigated. Since the properties of SFRC differ depending on geometry, tensile strength, aspect ratio and amount of steel fibres, it is not possible to test every single concrete mix with respect to its high-cycle performance.
In order to gain deeper understanding of the composite material SFRC under high-cyclic loading, the German Science Foundation (DFG) has launched a comprehensive priority program (SPP2020). The project described here is part of this priority program. The overall goal of this project is to develop an efficient experimental–numerical approach that will allow researchers and engineers to combine short-term tests, such as low-cycle flexural fatigue tests, with numerical models to evaluate the performance of high-performance SFRCs under high-cycle loading without the need for time-consuming high-cycle fatigue tests. To calibrate the numerical models, the evolution of damage in the SFRC has to be known. Therefore, different damage indicators such as the evolution of consumed energy per load cycle, the evolution of stiffness or the evolution of acoustic emission events are described, which can be used to calibrate damage functions in the numerical model. A comparable approach was described by [21] for polymer fibre-reinforced concrete and by de Smedt [22] and Boulekbache et al. [23] for normal-strength fibre-reinforced concrete. In the former two contributions, a sectional analysis approach was applied as numerical models.
This paper describes specific, low-cycle, flexural tests on steel fibre-reinforced concrete with three different amounts of steel fibres, the development of various damage indicators including acoustic emission, and a numerical phase-field approach to fracture. Finally, the ability of the model to calculate the load crack-opening behaviour and the damage evolution thereof is shown.
The phase-field approach for fracture has recently become popular because of its prediction capability for complex crack propagation, see [24,25,26,27,28]. The works on the phase-field models for ductile fracture (see [29,30,31,32,33]) are followed for the development of an elastoplastic phase-field model. In our recent work using the small-strain elastoplastic phase-field model for fractures, see [34,35,36], the single steel-fibre pullout tests are elaborated experimentally and numerically to set up a framework for the study of fibre’s influence on the overall material behaviour of HPC. In this contribution the formulations presented in [37] are further developed to capture the numerical response of three-point bending beam tests for steel fibre-reinforced HPC.
This paper, as well as the underlying research effort, is divided into two main areas. Within the scope of the first area, the material behaviour is investigated experimentally with the aim of determining suitable indicators to describe damage development. For this purpose, as well as for calibrating and validating the numerical model, these indicators are needed to describe the development of the deterioration in the post-cracking phase [38,39,40,41,42,43]. For the determination of these damage indicators as well as for characterizing low- and high-cycle fatigue behaviour, unloading–reloading cycles are necessary. The second area is focused on the development of a phenomenological model for the prediction of the complex elastoplastic material behaviour of steel fibre-reinforced HPC under fatigue.

2. Materials and Experimental Methods

2.1. Concrete Mixtures and Fibres

All tests described here were performed on high-performance concrete bending beams of dimension of 700 × 150 × 150 mm 3 . The concrete composition used was developed as a reference mixture within the German priority programme SPP 2020; the exact component breakdown can be found in Table 1. As listed there, it is based on a Portland cement of type CEM I 52.5R. As aggregate, quartz sand is used in grain fractions from 0 to 2 mm and basalt gravel for those from 2 to 8 mm in grain size. Furthermore, a superplasticizer and stabilizer are added for better fresh concrete workability. The w–c ratio of the mixture is 0.352. In addition to the tests described in this contribution, comparable sets of tests are performed using ultra-high-performance concrete (UHPC) with high-strength, short, steel fibres.
The mechanical properties of the concrete mix shown in Table 2 were determined on standard cylindrical specimens. Stripping was performed one day after casting and testing started 28 days after production. All specimens were stored in water until testing.
The experiments were carried out on specimens containing 0/23/57/115 kg/m 3 hook-end macro steel fibres of type 3D 65/60 with a tensile strength of 1150 MPa, kindly provided by Bekaert. These correspond to contents of 0/0.3/0.75/1.5% by volume, respectively. A total of 12 bending beam specimens are investigated in the test series described here, with three specimens per examined fibre set.

2.2. Experimental Setup and Loading

All specimens have been stored under water for 28 days until testing. In order to investigate damage development and, at the same time, to ensure a basic comparability between tests, three-point bending tests are carried out on notched specimens according to EN 14 651, with additional periodic unloading–reloading cycles. The test setup is shown in Figure 1.
In the tests accompanied by acoustic emission (AE) recording, piezoelectric AE-sensors type VS150-K by Vallen Systeme are attached to both sides of the beam along the longitudinal axis, in two planes containing four sensors each. The first plane is intended to provide the most accurate two-dimensional resolution possible of the primary crack propagation over a short distance of 20 mm from the notch of the beam. The second plane of sensors is intended to quantify the size of the three-dimensional crack process zone over a longer distance of 100 mm to the other side of the notch. The sensors are connected to a Vallen AMSY-6 system through pre-amplifiers with a gain level of 34 dB. The most relevant calibration settings for the system are summarized in Table 3.
As mentioned earlier, the experiments are performed with periodic unloading–reloading cycles. The cycles at the beginning of the test (CMOD values of up to 0.05 mm ) are carried out at crack-opening increments of 0.01 mm in order to make the crack initiation traceable as precisely as possible. For larger crack openings, the individual crack opening target value is doubled for each, following the unloading cycle, namely CMOD values of 0.1/0.2/0.4/…/3.2 mm. All tests were stopped at a CMOD of 4.0 mm. All phases with an increasing crack mouth opening displacement CMOD (loading phases) were controlled using the CMOD. All phases with a decreasing CMOD (unloading phases) were force-controlled in order to be able stop at a given lower load limit automatically.

2.3. Damage Indicators

The aim of this contribution is to experimentally quantify and numerically simulate not only the envelope of the load-CMOD curve but also to quantify the ongoing deterioration within the fibre-modified high-performance concrete. For this purpose, different damage indicators can be derived from the experiments. As stated above, unloading curves are necessary to differentiate between elastic crack openings, plastic crack openings and to further determine residual stiffnesses, especially after the occurrence of the first crack. In the literature [22,23,39,41,42,43] different indicators are proposed to quantify damage development throughout a low-cycle or high-cycle fatigue test. Mechanical indicators that are often mentioned are the plastic part of the crack-opening, the residual stiffness (sometimes referred to as compliance) and the energy that is consumed to open the crack. With regard to the damage-effective energy, different approaches must be distinguished. For example, in static tests, an integral over the load-CMOD curve can be used to determine the energy consumed [39,41]. For force-controlled, high-cycle fatigue tests, the damage energy can be interpreted as the portions of the load-CMOD curve enclosed by the hysteresis loops, as shown in Bode et al. [46]. A second approach for such tests is to interpret the trapezoid enclosed by the upper and lower load levels and initial and residual stiffness as the damage energy, an approach taken, for example, in Gebuhr et al. [41], De Smedt et al. [43]. The latter approach is transferable to fatigue tests with a low number of load cycles, such as those investigated in this paper, since residual stiffnesses can also be determined here, thus enabling better comparability with high-cycle fatigue tests.
In this contribution, three damage indicators are derived and displayed from the experiments and subsequently used in the numerical model. In Figure 2 these indicators are sketched.
For all indicators, a reference frame has to be set. In this contribution, the total CMOD for each unloading cycle (CMOD tot , i ) is applied, defined as the intersection of the unloading and consecutive reloading branch of the load-CMOD curve of each individual load cycle. The CMOD tot at the beginning of each unloading cycle is defined according to CMOD tot , i = CMOD pl , i + CMOD el , i , whereby the plastic CMOD pl , i is additionally used as a damage indicator. Another well-known damage indicator is the gradient module of the i t h unloading–reloading cycle, the residual stiffness s i , calculated as
s i = Δ F i CMOD tot , i CMOD pl , i ,
where Δ F describes the load difference between the CMOD tot , i - and CMOD pl , i -assigned points of the curve. In a second step this stiffness can be transformed to a dimensionless damage value, q i , as used in the numerical model, by normalizing it using the initial stiffness s 0 , i.e.,
q i = s 0 s i s 0 .
The dimensionless damage parameter q i is then used for the calibration of the numerical model to formulate the degradation function g ( q i ) . The last indicator is the plastic or damage energy, W pl , i . It is defined as the integral of the load-CMOD curve up to CMOD tot , i ; reduced by the elastic energy at this point, it reads
W pl , i = 0 CMOD tot , i F d CMOD CMOD el , i * F max , i 2 .
For the calibration of the numerical model of the load-CMOD curve, plastic crack opening CMOD pl , i and the residual stiffness are of primary importance, as the calculation of the plastic energy is based on them. The challenge for the numerical model is to correctly predict not only the named indicators but also their development with an increasing fibre content.
Apart from the aforementioned mechanical damage indicators, an analysis of the acoustic phenomena can also be used to assess and describe the damage process. The analysis of the data is divided into different levels, all of which revolve around the analysis of the so called hits, i.e., the signal shape after a predefined threshold value has been exceeded. The first level considered here is the direct evaluation of the time-domain characteristics for hits recorded by the individual microphones, such as the maximum amplitude and overall energy of the hit, as depicted in Figure 3. These values can be used to determine the general activity in the test specimen, as shown e.g., in De Smedt et al. [43], Nguyen [47].

2.4. Numerical Model

For the numerical solution, we used the small-strain, elastoplastic phase-field model for fracture in reinforced high-performance concrete, see Pise et al. [35], Storm et al. [36], Pise et al. [37]. In the context of small strain, the displacement field is u and the phase-field parameter is q 0 , 1 , which serves as a damage parameter, see Miehe et al. [25]. An additive structure of the macroscopic stored energy function per unit volume ψ is introduced as
ψ = v HPC ψ HPC ε , ε p , HPC , q , q , α HPC + v F ψ F ε , M , e p , F , α F with ε e , HPC : = ε ε p , HPC ,
where the variables α HPC and α F are the equivalent plastic strains in the HPC phase and in the fibre phase, respectively. Elastic ε e , HPC and plastic ε p , HPC strains in the HPC phase are the parts of the total strain tensor ε that can be expressed using the symmetric displacement gradient, s u , as ε = s u = 1 2 u + T u . It is observed in the experiments that during the failure of steel fibre-reinforced HPC the fracture mainly occurs only in the HPC phase. On that account, the phase-field parameter q [ 0 , 1 ] is considered only in the HPC phase, representing the damage therein. Its gradient is denoted by q . The symbols v HPC and v F describe the volume fraction of the HPC phase and fibre phase, respectively. The are conserved by the constraint v HPC = 1 v F .
The energy function ψ HPC describes the mechanical behaviour of the HPC phase, cf. [32], reads
ψ HPC = g q ψ 0 e + , HPC + ψ 0 p , HPC ψ c , HPC + ψ 0 e , HPC + ψ c , HPC + 2 ψ c ζ l 1 2 l q 2 + l 2 | | q | | 2 .
The specific critical fracture energy ψ c , HPC > 0 controls the crack threshold and the vanishing length scale parameter l 0 achieves a sharp crack interface. The parameter ζ controls the post-critical stress softening due to the evolution of the fracture and g q , denoting a degradation function, depends on the phase-field parameter q. A positive part ψ 0 e + , HPC ( ε e , HPC ) and a negative part ψ 0 e , HPC ( ε e , HPC ) of the reference elastic energy function ψ 0 e , HPC ( ε e , HPC ) as proposed in Amor et al. [48], read
ψ 0 e + , HPC ( ε e , HPC ) = κ tr [ ε e , HPC ] + 2 / 2 + μ | | dev ε e , HPC | | 2 and ψ 0 e , HPC ( ε e , HPC ) = κ tr [ ε e , HPC ] 2 / 2 .
This formulation ensures that the evolution of damage only occurs due to the volume expansion of material. Therein, μ and κ are the Lam e ´ coefficients of the HPC phase and Macaulay’s notation describes the operation ± = 1 / 2 ( ± | | ) . The considered plastic energy, ψ p , HPC , reads
ψ p , HPC α HPC = g q ψ 0 p , HPC with ψ 0 p , HPC = y 0 HPC α HPC + 1 2 h HPC ( α HPC ) 2 ,
which depends on the equivalent plastic strain α HPC , yield stress y 0 HPC and hardening parameter h HPC for the HPC phase, respectively.
The energy function for the fibre phase ψ F characterizing the embedded steel fibre in a preferred fibre direction a with | | a | | = 1 reads
ψ F ε , M , e p , F , α F = ψ e , F ε , M , e p , F + ψ p , F α F where e e , F = e F e p , F with e F = ε : M ,
where α F is the equivalent plastic strain for the fibre phase. The total strain tensor e F of the fibre phase along the preferred fibre directions can be decomposed into elastic e e , F and plastic e p , F parts. Structural tensor M can be formulated using the preferred fibre direction a as M : = a a , cf. [49]. With this the one-dimensional elastoplasticity problem representing the three-dimensional elastoplasticity for the embedded steel fibre, oriented in a preferred fibre direction, a is formulated. Therein, an elastic ψ e , F and a plastic ψ p , F energy function are considered as
ψ e , F ε , M , e p , F = 1 2 E F e F e p , F 2 and ψ p , F α F = y 0 F α F + 1 2 h F ( α F ) 2 ,
where E F , y 0 F and h F are the elastic moduli, initial yield stress and hardening parameter for the fibre phase, respectively. The stress tensor additively, composed by the stress tensor σ HPC for the HPC phase and σ F for the fibre, phase reads
σ : = v HPC σ HPC + v F σ F .
The stress tensor σ F for the fibre phase is defined as
σ F : = σ F a a with σ F = e e , F ψ F = E F e F e p , F .
A stress tensor for the HPC phase σ HPC can be obtained by
σ HPC : = ε e , HPC ψ HPC = g q κ HPC tr ε e , HPC + I + 2 μ HPC dev ε e , HPC σ 0 + , HPC + κ HPC tr ε e , HPC I σ 0 , HPC ,
where σ 0 + , HPC and σ 0 + , HPC are the positive and the negative part of the effective stress tensor σ 0 HPC for the HPC phase. Here, the degradation function g q is a data-driven function that depends on the evolution of the phase-field parameter q. Two different data sets are considered for the degradation function, g + q , in tension, and g q , in compression. These degradation functions are calibrated using the experimental data for the three-point bending beam test for pure HPC at low cycle. For the computation of the governing equation for the phase-field parameter q, the degradation function g q = 1 q 2 is considered. It is calculated from the derivative of the macroscopic stored energy function ψ with respect to the phase-field parameter q, which results in
q l 2 Div q 1 q H HPC = 0 ,
where the maximum value of the local history field H , which ensures the irreversibility of the evolution of damage, is given by
H HPC : = max s [ 0 , t ] ζ v HPC [ ψ 0 e + , HPC + ψ 0 p , HPC ] ψ c , HPC 1 .
The unique behaviour of concrete in tension and in compression is captured by using different parameters for the critical fracture energy in tension ψ t c , HPC and in compression ψ c c , HPC , for details see Storm et al. [36]. They are differentiated from each other by the sign of the first invariant of the stress tensor tr σ HPC for the HPC phase. Two different yield criteria are considered for the description of non-linear behaviour of the HPC and fibre phases. This gives the flexibility to use different yield criteria to predict the elastoplastic behaviour of the reinforcement, as well as of the concrete materials. The non-associative Drucker–Prager yield criteria is used for the HPC phase, which is capable of capturing the asymmetric tension–compression property of concrete. A potential function of Drucker–Prager yield criteria, ϕ p HPC , is defined by cf. [50,51],
ϕ p HPC = 1 2 | | dev σ 0 HPC | | β p tr σ 0 HPC κ p HPC with κ p HPC : = α HPC ψ 0 p , HPC = y 0 HPC + h HPC α HPC ,
and the plastic potential function, ϕ n HPC , is considered as
ϕ n HPC = 1 2 | | dev σ 0 HPC | | β n tr σ 0 HPC ,
where β p and β n are the material parameters. Note that these potentials are the functions of the effective stress tensor σ 0 HPC for the HPC phase. The equation describing the evolution of plastic strains ε ˙ p , HPC , equivalent plastic strains α ˙ HPC and expressed in terms of the incremental plastic consistency parameter λ p , HPC , and the Kuhn–Tucker conditions for the HPC phase read
ε ˙ p , HPC = λ p , HPC ϕ n HPC σ 0 HPC , α ˙ HPC = λ p , HPC and ϕ p HPC 0 , λ p , HPC 0 , λ p , HPC ϕ p HPC = 0 .
To describe the non-linear behaviour along the preferred fibre direction, the one-dimensional von Mises yield criterion is used, i.e.,
ϕ F σ F , κ p F = | σ F | κ p F with κ p F : = α F ψ p , F = y 0 F + h F α F .
The evolution equation of the plastic strain e ˙ p , F and the equivalent plastic strain α F using Kuhn–Tucker conditions to describe the loading and unloading conditions for the fibre phase are
e ˙ p , F = λ p , F ϕ F σ F , α ˙ F = λ p , F and ϕ F 0 , λ p , F 0 , λ p , F ϕ F = 0 ,
where λ p , F denotes the incremental plastic consistency parameter for the fibre phase. The system of equations contains the balance of linear momentum using Equation (10) and the governing equation for the phase-field parameter, i.e., Equation (13), in domain B along with boundary conditions σ · n = t and q · n = 0 on the surface of the domain B . For the analysis of failure in reinforced HPC, the presented numerical model is implemented in the framework of the finite element method and solved using the incrementally decoupled updates using a staggered scheme, cf. [25,30].

3. Experimental and Numerical Results

3.1. Load-CMOD Curves

The load-CMOD curves displayed in Figure 4 show the effects of fibre reinforcement on the ductility and load-bearing capacity. All specimens tested show similar curves up to a CMOD of approx. 0.015 mm. From this point, plastic deformation starts, as the crack initiates. The mixture without fibres (black line) has only very limited ability to bridge this crack opening. Failure occurs nearly immediately after the maximum load has been reached. The mixture with a fibre content of 23 kg/m 3 and 0.3 vol.-% (green line) shows typical subcritical behaviour. The forces transferred over the crack decrease continuously after the flexural strength is reached. The blue line represents the behaviour with a fibre content of 57 kg/m 3 respectively 0.5 vol.-%. This line shows typical critical behaviour, which means that force transfer can be maintained by the fibres after the crack opening starts and the flexural strength is preserved, even at higher crack openings. In this curve it can also be seen that it is important to distinguish between the load that can be absorbed at first crack, which is reached when the force CMOD curve deviates from a linear elastic slope, and the maximum load that can be absorbed. In the case of strain-hardening behaviour as shown here, this force is much higher than the force at first crack, meaning that the high flexural strength is attributable to the high number of fibres and their load bearing capacity. This difference is more pronounced if the fibre content is increased to 115 kg/m 3 or 1.5% by volume. Here, the first crack still occurs at about 25 kN, but flexural strength is reached at about 58 kN at CMOD tot , i of about 0.6 mm .

3.2. Mechanical Damage Indicators

As prior research has shown [39,41,43], one feasible indicator for the overall deterioration process is the evolution of the residual stiffnesses, i.e., the gradient modules of the individual load cycles. Figure 5 depicts the development of residual stiffnesses of concretes containing different fibre contents. To get a better resolution of the first unloading cycles up to a CMOD of 0.06 mm, which is the failure point of the fibre-free specimens, a detail view is added up for this CMOD, which will be given for all following figures depicting damage indicators.
Since it holds true for all beams that the first unloading cycle initiated at a CMOD tot of 0.01 mm is regraded as the linear elastic part of the load-CMOD curve, the first stiffness determined can be understood as the initial stiffness of the uncracked material. Since the fibres in the specimen are not activated at this point, all stiffnesses are close to each other, on average, at approx. 1300 kN/mm. As stated above, the specimens start to develop plastic deformations at a CMOD of approx. 0.015 mm. While the crack opening itself leads to a decrease in stiffness for all fibre contents, the steepness of this decrease is determined by the respective fibre content, with the fibre-free specimens degrading fastest. Furthermore, a higher fibre quantity seems to lead to a more pronounced scatter in the curve. Except for the fibre-free specimens, which fail at a crack-mouth opening of about 0.06 mm, the fibre content-dependent degradation rate continues up to crack openings of about 2 mm. From this point onward, the supercritical fibre contents of above 57 kg/m 3 are, again, closer together, at a stiffness level of about 100 kN/mm. Beams with a subcritical fibre content of 23 kg/m 3 achieve residual stiffnesses of about 75 kN / mm at these crack openings.
The development of the plastic crack opening CMOD pl at the lower reversal point of each load cycle as shown in Figure 6 also indicates a clear influence of the fibre content on the behaviour of the specimens starting shortly after the initiation of the damage to the concrete matrix. This becomes evident from a CMOD tot of 0.03 mm, above which the fibre-free specimens begin to show significantly larger plastic crack openings than the beams with fibres. The plastic deformations of the remaining beams first start to scatter significantly more up to this value. From a CMOD tot of 0.4 mm the specimens plastic crack openings start to order according to their fibre contents, with the beams with the 115 kg / m 3 showing the smallest values. This trend increases in the rest of the test, up to 3.2 mm total crack opening, where these beams show an average CMOD pl of 2.6 mm. The beams with 23 kg/m 3 fibres are at 2.99 mm at this point, i.e., on average, with approx. at 12% larger plastic crack openings. The difference in these quantities is primarily due to the greater residual forces that can be transferred by the material as the fibre quantity increases. Since a convergence of the stiffnesses between the different fibre contents takes place with larger crack openings, the plastic deformation is correspondingly smaller. Overall, it is also noticeable that the plastic parts of the crack openings in the tests on the beams with a fibre content of 115 kg / m 3 start to scatter more than with the lower fibre contents.
Figure 7 shows the progression of the absorbed damage energy. Here, the increasing absorption capacity of the beams with increasing fibre contents becomes clear. After the beams again show a behaviour, as with the aforementioned damage indicators, up to a crack opening of 0.3 mm—including the progressive fanning out into a wider dispersion of the calculated energies—from this point onwards, they again arrange themselves according to their fibre content. This process continues up to a crack opening of 0.8 mm. From this point on, the energy for an increase in crack opening rises approximately proportionally to the fibre content.

3.3. Acoustic Emission Measurements

Figure 8 shows an exemplary experiment on a bending beam with 57 kg/m 3 fibres, which corresponds to 0.75% by volume. Here, both the progression of CMOD as well as load are given over time. These progressions are now compared to the the data collected with the aid of the AE-sensors. Firstly, the AE-data is given via the relative number of raw hits detected after each point in time in relation to the total number of recorded hits over the course of the test. Secondly, the results are normalized using the respective energy of each of those hits to calculate cumulative energy.
The results are divided into distinct phases, starting with the pre-crack phase up to a CMOD of about 0.02 mm at about 800 s. Here, the crack initiation phase starts—mechanically defined by the loss of stiffness and beginning of plastic deformation, as described above. In the visualization above, this loss of stiffness can be seen in the gradient change of the residual force marked with a red circle. Following this point, micro-cracks begin to accumulate into a macro-crack, resulting into non-linearly increasing CMOD values. In this phase, the number of hits registered rises significantly. Interestingly, the beginning of the increase in total hits corresponds to the first significant increase in overall level of applied force after the beginning of the crack initiation phase. This effect however seems to have nearly no impact on the relative cumulative acoustic energy registered. While the absolute count of hits surpasses 60% of the total number of hits registered over the course of the experiment at the 3800 s mark of the test, the cumulative energy of the hits stays well below 5% of the total energy of hits registered by the microphones. The energy first starts to rise shortly after the maximum force is reached at a CMOD of approximately 0.8 mm, again labeled with a red circle in the curve. This point marks the beginning of a macro-crack-opening phase. As the acoustic events start to emit more energy, the overall load-bearing capacity of the beam degrades. This indicates a disparity between the first activation of the fibre and the minimum fibre pull-out necessary to mechanically activate the fibre.

3.4. Simulations of Three-Point Bending Tests at Low Cycle for Reinforced HPC

In this section, the simulation results of three-point bending beam tests at low cycle are reported. For these simulations, the mechanical properties of HPC are taken from the cumulative results of various partners within SPP 2020 program, for comparison, [34,44,45] (see Table 2). The material parameters and interpolated values of the data-driven degradation functions for tension g + ( q ) and for compression g ( q ) are calibrated simultaneously for the uniaxial cyclic tension test, uniaxial cyclic compression test and three-point bending test for pure HPC at low cycle. The mechanical properties of HPC and the calibrated parameters used for the simulation are listed in Table 4.
The calibrated values for the interpolation of degradation functions for HPC in tension g + ( q ) and compression g ( q ) , with respect to corresponding values of the phase-field parameter, are listed in Table 5.
The geometry and the boundary conditions of the three-point bending beam test are shown in Figure 1a. All the dimensions shown in Figure 1a are taken as per the European Standard EN 14651. For all simulations, a beam of reinforced HPC with fibre content of 23 kg/mm 2 , i.e., 0.3 vol.-%, is considered. The three-dimensional state is considered by a plane stress approximation. Therefore the x-y-cross section of the beam is discretized by a finite element mesh with 3673 eight-noded hexahedral elements and unit-length thickness. A Maximum element size of 0.5 mm is taken for elements in the area surrounding the notch at the middle of the beam. As shown in Figure 1a, load is applied by a vertical displacement boundary condition, u t , on the top surface of the beam, and the reaction forces of the constrained nodes at this tip are computed. To calculate the crack mouth-opening displacement (CMOD) the relative horizontal displacement of the opposite corners of the notch at the bottom of the beam is used, see Figure 1a. The distribution of stresses in the x-direction, σ x (in GPa), the equivalent, plastic strain. α HPC , the phase-field parameter, q, at CMOD = 0.018 mm are shown in Figure 9a,c,e, and at CMOD = 0.108 mm in Figure 9b,d,f, respectively.
In Figure 10a, the experimental and numerical load-CMOD curves are plotted. Therein, a scatter band of the experimental data and a representative experimental curve are shown in the gray area and by the black dotted curve, respectively. The numerical load-CMOD curve for pure HPC—the blue curve in Figure 10a—and reinforced HPC for a preferred fibre direction along x-axis—green curve in Figure 10a— represent the lower and upper bound for the load-CMOD curves for different combinations of preferred fiber direction.
The load-CMOD curve using five preferred fibre orientations—red curve in Figure 10a— between the angles of −10° to 10° (see Figure 10b) and using 24 preferred fibre orientations—cyan curve in Figure 10a—between the angles of −90° to 90° (see Figure 10c) are plotted in Figure 10a. For the analysis, load-CMOD curves for the experiment and for the simulation using five preferred fibre orientations are compared in Figure 10d. The development of the residual stiffness of the beam for simulation using five preferred fibre orientations is calculated for the each loading cycle, using the procedure explained in Figure 2b. For the calculation, the value of load at midpoint of the unloading and consequent reloading branch is calculated (see top point of dotted red vertical line in Figure 10d). The calculated residual stiffness-CMOD curves for the experiment and for the simulation using five preferred fibre orientations are plotted in the black and red circles, respectively, in Figure 10e. These calculated values of residual stiffness are used to get the interpolated residual stiffness-CMOD experimental, i.e., dotted black curve, and numerical curves, i.e., solid red curve. These interpolated residual stiffness–CMOD curves in Figure 10e are compared for the validation of the capabilities of the presented numerical model.

4. Discussion

Steel fibres significantly influence the behaviour of the otherwise very brittle HPC. While pure HPC fails after minimal crack openings, specimens with higher fibre contents can still transfer forces exceeding the initial crack load significantly, even at crack openings of more than 3 mm. As shown in Figure 5 and Figure 7, the plastic CMOD, residual stiffness and damage energy are suitable indicators to not only describe the deterioration process but also to differentiate between different fibre contents. Especially, the progression of the residual stiffness shows significantly higher values and a smoother transition to the opening of the macro-crack, especially in the crack-initiation phase from a CMOD of 0.02 mm for high fibre contents. The influence of the fibre quantity is also evident in the progress of the plastic fractions of the crack openings as well as the consumed energy. The results and development of damage indicators displayed comply well with studies in which similar loading–unloading schemes were used. For example, de Smedt [22] shows similar damage evolutions for two different types and lesser contents of hooked-end steel fibres in normal strength concrete. The tests performed here, therefore, show a generalizability to greater fibre contents and higher compressive strengths of the concrete. Results published by Boulekbache et al. [23] regarding residual stiffness and consumed energy are confirmed. The trends for all mentioned damage indicators as well as fibre content are comparable. Their results furthermore show positive effects of an increasing aspect ratio of the fibres as well as compressive strength of the concrete. However, the absolute values presented here are not directly comparable to Boulekbache et al. [23], since four-point bending tests are used there, resulting in multiple cracks in the tensile zone. In contrast, in the tests presented here, a single primary crack is induced in the center of the beam due to the weakened cross-section.
Nevertheless, scattering of the damage indicators is noticeable. However, it requires further investigation as to whether this is directly attributable to the amount of fibres, e.g., due to a more inhomogeneous fibre distribution resulting from comparatively poorer fresh concrete properties with increasing fibre quantities. The analysis of the acoustic emission provides additional insights into the damage behaviour. The comparison between the development of number of hits to the development of the energy of the acoustic events seems suitable to make the phases of crack initiation and macro-crack opening differentiable. The crack initiation is characterized by high numbers of hits with low energy content, whereas macro-crack opening is characterized by high energies at those acoustic events. This distinction between a crack-initiation phase and a crack-opening phase has first been introduced by Schorn [52] and Kopp [53]. De Smedt [22] reports that the evolution of acoustic emission can be used as alternative damage indicator. This finding can be supported here.
A phenomenological material model incorporating phase-field approach for fracture representing the damage behaviour of reinforced HPC at low cycle is presented. To describe the nonlinear features of HPC the Drucker–Prager plasticity model is formulated for the HPC phase. Additionally, a one-dimensional von Mises plasticity model is used to capture the non-linearity along preferred fibre direction. The Drucker–Prager plasticity model for the HPC phase, combined with different values of the critical fracture energy in the phase–field model, empower the presented model to calibrate the different parameters for concrete in tension and compression. The comparison of experimental and numerical results of three-point bending beam tests at low cycle shows the efficiency of the presented model. It is able to reproduce the main characteristics of the experimental results using measured mechanical properties and calibrated material parameters. For the validation, the degradation of the residual stiffness calculated using the numerical results and experimental data are compared. It can be observed that the model is not only capable of predicting the load-CMOD curves in terms of the load envelop, but also to efficiently predict the deterioration development in terms of residual stiffness and plastic CMOD during unloading-reloading cycles.

5. Conclusions

In this contribution a combined experimental–numerical approach is outlined aiming at the assessment of the high-cycle performance of steel fibre-reinforced HPC and UHPC. As a first step, low-cycle flexural fatigue tests on concrete with different fibre contents are conducted and different damage indicators, including acoustic emission analysis, are derived. Simultaneously, a phase-field approach for the fracture of steel fibre-reinforced HPC is presented.
From the experimental work, it is found that the evolution of stiffness and the plastic fraction of crack-mouth opening are damage indicators that show clear characteristics depending on the fibre content and are suitable for the calibration of numerical models. The absorbed damage energy is highly selective with respect to different fibre contents, but its applicability for calibration is numerically complex. The acoustic emission method allows the differentiation of different fracture processes, such as crack initiation and crack opening.
Furthermore it is shown, that the phenomenological material model incorporating a phase field approach for fracture in the pure HPC can be extended to predict fracture in steel fibre-reinforced HPC. The main feature is that, in contrast to incorporating the real spatial fibre distribution, the formulation using a structural tensor based on effective overall fibre distributions reduces the effort of their consideration in finite element discretization. This enables the feasible study of the influence of the different orientations and distribution of the embedded fibres on the overall material behaviour of HPC, e.g., load-CMOD curves. Additionally, the material model is not only capable of predicting the load-CMOD curve but also to predict damage evolution in terms the residual stiffness and the plastic CMOD. To study the capabilities, the presented model will be further developed for higher fibre contents that provide supercritical behaviour. Likewise, the modeling of higher cycle numbers is in progress and a comparison of damage development in low-cycle and high-cycle fatigue will be the future steps. Initial results, see [41], promise that the damage evolution is comparable.

Author Contributions

Methodology (experimental aspect), G.G. and S.A.; methodology (modelling aspects), M.P., D.B. and J.S.; writing—original draft preparation, M.P. and G.G.; writing—review and editing, S.A. and D.B.; supervision, project administration and funding acquisition, S.A., D.B. and J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—353513049 (AN1113/2-2, BR5278/2-2, SCHR570/32-2) within the DFG Priority Programme 2020 “Cyclic deterioration of high-performance concrete in an experimental-virtual lab”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Computing time provided on the supercomputer magnitUDE, funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—263348352 (INST 20876/209-1 FUGG, INST 20876/243-1 FUGG), is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reinhardt, H.; Naaman, A. (Eds.) High Performance Fiber Reinforced Cement Composites; RILEM Publications SARL: London, UK, 1996. [Google Scholar]
  2. Naaman, A.E. Half a Century of Progress Leading to Ultra-High Performance Fiber Reinforced Concrete: Part 1-Overall Review; RILEM Publications SARL: Rio de Janeiro, Brazil, 2012. [Google Scholar]
  3. Tiberti, G.; Chiriotti, E.; Plizzari, G.A. Twenty Years of FRC Tunnel Segments Practice: Lessons Learnt and Proposed Design Principles; International Tunnelling Association (ITA): Montreal, QC, Canada, 2016. [Google Scholar]
  4. Coughlan, D.; Diez, R.; Comins, J.; Stärk, A. Crossrail project: Use of sprayed concrete tunnel linings on London’s Elizabeth line. Proc. Inst. Civ. Eng.- Civ. Eng. 2017, 170, 39–46. [Google Scholar] [CrossRef] [Green Version]
  5. Abbas, S.; Soliman, A.; Nehdi, M. Structural Behaviour of Ultra-High Performance Fibre Reinforced Concrete Tunnel Lining Segments. In Proceedings of the FRC 2014 Joint ACI-fib International Workshop Fibre Reinforced Concrete: from Design to Structural Applications, International Federation for Structural Concrete, Lausanne, France, 24–25 July 2014. [Google Scholar]
  6. Blasini, G.; Kodra, L.; Bragard, C.; Jackson, P.; de la Fuente, A.; Pannetier, T.; Dehn, F.; Plizzari, G.A.; Dobashi, H.; Rinaldi, Z.; et al. fib Bulletin 83. Precast Tunnel Segments in Fibre-Reinforced Concrete; International Federation for Structural Concrete: Lausanne, France, 2017. [Google Scholar]
  7. Aghniaey, N. Behaviour of Self Consolidating Steel Fiber Reinforced Concrete Beams under Reversed Cyclic Loading. M.Sc. Thesis, University of Ottawa, Ottawa, ON, Canada, March 2013. [Google Scholar]
  8. Jiuru, T.; Chaobin, H.; Kaijian, Y.; Yongcheng, Y. Seismic Behavior and Shear Strength of Framed Joint Using Steel-Fiber Reinforced Concrete. J. Struct. Eng. 1992, 118, 341–358. [Google Scholar] [CrossRef]
  9. Abbass, A.; Abid, S.; Özakça, M. Experimental Investigation on the Effect of Steel Fibers on the Flexural Behavior and Ductility of High-Strength Concrete Hollow Beams. Adv. Civ. Eng. 2019, 2019, 1–13. [Google Scholar] [CrossRef]
  10. Tsonos, A.D. Steel fiber high-strength reinforced concrete: A new solution for earthquake strengthening of old R/C structures. WIT Trans. Built Environ. 2009, 104, 153–164. [Google Scholar] [CrossRef] [Green Version]
  11. Sørensen, E.; Westhof, L.; Yde, E.; Serednicki, A. Fatigue Life of High Performance Grout for Wind Turbine Grouted Connection in Wet or Dry Environment; NORSK BETONGFORENING: Oslo, Norway, 2011. [Google Scholar]
  12. Technical Report 34 (TR34): Concrete Industrial Ground Floors; Technical Report; The Concrete Society: Surrey, UK, 2016.
  13. Birely, A.; Park, P.; McMahon, J.; Shi, X.; Rew, Y. Fiber Reinforced Concrete for Improved Performance of Transportation Infrastructure; Technical Report; Federal Highway Administration: Washington, DC, USA, 21 June 2018. [Google Scholar]
  14. Greybeal, B. Behavior of Field-Cast Ultra-High Performance Bridge Deck Connections Under Cyclic and Static Structural Loading; Federal Highway Administration: Washington, DC, USA, 2010. [Google Scholar]
  15. Perry, V.; Weiss, G. Innovative Field Cast UHPC Joints for Precast Bridge Decks. Design, Prototype Testing and Projects. In Designing and Building with UHPFRC; John Wiley & Sons, Ltd: Hoboken, NJ, USA, 2011; Chapter 28; pp. 421–436. [Google Scholar] [CrossRef]
  16. McMahon, J.A.; Birely, A.C. Experimental Performance of Steel Fiber Reinforced Concrete Bridge Deck. J. Bridge Eng. 2018, 23, 04018074. [Google Scholar] [CrossRef]
  17. Akhnoukh, A.K.; Buckhalter, C. Ultra-high-performance concrete: Constituents, mechanical properties, applications and current challenges. Case Stud. Constr. Mater. 2021, 15, e00559. [Google Scholar] [CrossRef]
  18. Graybeal, B. Design and Construction of Field-Cast UHPC Connections; Technical Report; Federal Highway Administration: Washington, DC, USA, 2014. [Google Scholar] [CrossRef]
  19. Aboubakr, A.A.M. Behaviour Study of Grouted Connection for Offshore Wind Turbine Structures with Brittle Cement Based Grouts. Ph.D. Thesis, Fachbereich Bauingenieur- und Umweltingenieurwesen, Universität Kassel, Kassel, Germany, 2020. [Google Scholar]
  20. Schaumann, P.; Raba, A.; Bechtel, A. Fatigue behaviour of grouted connections at different ambient conditions and loading scenarios. Energy Procedia 2017, 137, 196–203. [Google Scholar] [CrossRef]
  21. Vrijdaghs, R. Creep of Synthetic Fiber Reinforced Concrete: A Multi-Scale and Two Phased Approach. Ph.D. Thesis, KU Leuven, Leuven, Belgium, 13 February 2019. [Google Scholar]
  22. De Smedt, M. Cyclic Behaviour of Steel Fibre Reinforced Concrete with Acoustic Emission-Based Multi-Scale Analysis. Ph.D. Thesis, KU Leuven, Leuven, Belgium, 22 December 2021. [Google Scholar]
  23. Boulekbache, B.; Hamrat, M.; Chemrouk, M.; Amziane, S. Flexural behaviour of steel fibre-reinforced concrete under cyclic loading. Constr. Build. Mater. 2016, 126, 253–262. [Google Scholar] [CrossRef]
  24. Bourdin, B.; Francfort, G.; Marigo, J. The variational approach to fracture. J. Elast. 2008, 91, 5–148. [Google Scholar] [CrossRef]
  25. Miehe, C.; Hofacker, M.; Welschinger, F. A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits. Comput. Methods. Appl. Mech. Eng. 2010, 199, 2765–2778. [Google Scholar] [CrossRef]
  26. Miehe, C.; Welschinger, F.; Hofacker, M. Thermodynamically consistent phase-field models of fracture: Variational principles and multi-field FE implementations. Int. J. Numer. Methods. Eng. 2010, 83, 1273–1311. [Google Scholar] [CrossRef]
  27. Kuhn, K.; Müller, R. A continuum phase field model for fracture. Eng. Frac.t Mech. 2010, 77, 3625–3634. [Google Scholar] [CrossRef]
  28. Ambati, M.; Gerasimov, T.; De Lorenzis, L. A review on phase-field models of brittle fracture and a new fast hybrid formulation. Comput. Mech. 2015, 55, 383–405. [Google Scholar] [CrossRef]
  29. Kuhn, C.; Noll, T.; Müller, R. On phase field modeling of ductile fracture. GAMM-Mitteilungen 2016, 39, 35–54. [Google Scholar] [CrossRef]
  30. Ambati, M.; Gerasimov, T.; De Lorenzis, L. Phase-field modeling of ductile fracture. Comput. Mech. 2015, 55, 1017–1040. [Google Scholar] [CrossRef]
  31. Borden, M.; Hughes, T.; Landis, C.; Anvari, A.; Lee, I. A phase-field formulation for fracture in ductile materials: Finite deformation balance law derivation, plastic degradation, and stress triaxiality effects. Int. J. Numer. Methods. Eng. 2016, 312, 130–166. [Google Scholar] [CrossRef] [Green Version]
  32. Miehe, C.; Aldakheel, F.; Raina, A. Phase field modeling of ductile fracture at finite strains: A variational gradient-extended plasticity-damage theory. Int. J. Plast 2016, 84, 1–32. [Google Scholar] [CrossRef]
  33. Alessi, R.; Ambati, M.; Gerasimov, T.; Vidoli, S.; De Lorenzis, L. Comparison of Phase-Field Models of Fracture Coupled with Plasticity. Comput. Methods. Appl. Sci. 2018, 46, 1–21. [Google Scholar]
  34. Gebuhr, G.; Pise, M.; Sarhil, M.; Anders, S.; Brands, D.; Schröder, J. Analysis and evaluation of single fibre pull-out behaviour of hooked steel fibres embedded in high-performance concrete for calibration of numerical models. Struct. Concr. 2019, 20, 1254–1264. [Google Scholar] [CrossRef]
  35. Pise, M.; Brands, D.; Sarhil, M.; Schröder, J.; Gebuhr, G.; Anders, S. Numerical calibration of elastoplastic phase-field modeling of fracture for experimental pullout tests of single steel fibres embedded in high-performance concrete. In Advances in Engineering Materials, Structures and Systems: Innovations, Mechanics and Applications, Proceedings of the 7th International Conference on Structural Engineering, Mechanics and Computation (SEMC 2019), Cape Town, South Africa, 2–4 September 2019; Zingoni, A., Ed.; Taylor & Francis Group: London, UK, 2019. [Google Scholar]
  36. Storm, J.; Pise, M.; Brands, D.; Schröder, J.; Kaliske, M. A comparative study of micro-mechanical models for fiber pullout behavior of reinforced high-performance concrete. Eng. Fract. Mech. 2021, 243, 107506. [Google Scholar] [CrossRef]
  37. Pise, M.; Brands, D.; Schröder, J.; Gebuhr, G.; Anders, S. Macroscopic model for steel fiber-reinforced high-performance concrete based on unit cell calculations. Proc. Appl. Math. Mech. 2021, 21, e202100180. [Google Scholar] [CrossRef]
  38. Barros, J.A.; Cunha, V.M.; Ribeiro, A.F.; Antunes, J. Post-cracking behaviour of steel fiber-reinforced concrete. Mater. Struct. 2005, 38, 47–56. [Google Scholar] [CrossRef]
  39. Germano, F.; Tiberti, G.; Plizzari, G. Post-peak fatigue performance of steel fiber-reinforced concrete under flexure. Mater. Struct. 2016, 49, 4229–4245. [Google Scholar] [CrossRef]
  40. Zhang, J. Fatigue Fracture of Fibre Reinforced Concrete. Ph.D. Thesis, Technische Universität von Dänemark, Kongens Lyngby, Denmark, 1998. [Google Scholar]
  41. Gebuhr, G.; Pise, M.; Sarhil, M.; Anders, S.; Brands, D.; Schröder, S. Deterioration development of steel fiber-reinforced high-performance concrete in low cycle fatigue. In Advances in Engineering Materials, Structures and Systems: Innovations, Mechanics and Applications, Proceedings of the 7th International Conference on Structural Engineering, Mechanics and Computation (SEMC 2019), Cape Town, South Africa, 2–4 September 2019; Zingoni, A., Ed.; Taylor & Francis Group: London, UK, 2019. [Google Scholar]
  42. Subramaniam, K.V.; Popovics, J.; Shah, S. Fatigue fracture of concrete subjected to biaxial stresses in the tensile ct region. J. Eng. Mech. 2002, 128, 668–676. [Google Scholar] [CrossRef]
  43. De Smedt, M.; Vrijdaghs, R.; Van Steen, C.; Verstrynge, E.; Vandewalle, L. Damage analysis in steel fibre-reinforced concrete under monotonic and cyclic bending by means of acoustic emission monitoring. Cem. Concr. Compos. 2020, 114, 103765. [Google Scholar] [CrossRef]
  44. Leusmann, T.; Basutkar, G.; Lunardelli, M. Characterizing the 3D mesostructure of high-performance concrete with the help of computed tomography. In Proceedings of the International Conference on Sustainable Materials, Systems and Structures (SMSS 2019)—Novel Methods for Characterization of Materials and Structures, Rovinj, Croatia, 18–22 March 2019; Gabrijel, I., Grosse, C., Skazlic, M., Eds.; RILEM: Paris, France, 2019. [Google Scholar]
  45. Scheiden, T.; Oneschkow, N. Influence of coarse aggregate type on the damage mechanism in high-strength concrete under compressive fatigue loading. Struct. Concr. 2019, 20, 1212–1219. [Google Scholar] [CrossRef]
  46. Bode, M.; Marx, S.; Vogel, A.; Völker, C. Dissipationsenergie bei Ermüdungsversuchen an Betonprobekörpern. Beton- Stahlbetonbau 2019, 114, 548–556. [Google Scholar] [CrossRef]
  47. Nguyen, T. Damage classification in concrete, mortar and cement-paste beams under bending by Acoustic Emission Technique (AET). Mater. Sci. Eng. 2020, 869, 072007. [Google Scholar] [CrossRef]
  48. Amor, H.; Marigo, J.J.; Maurini, C. Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments. Comput. Mech. 2009, 57, 1209–1229. [Google Scholar] [CrossRef]
  49. Boehler, J.P. Applications of Tensor Functions in Solid Mechanics; Springer: Vienna, Austria, 1987. [Google Scholar]
  50. Drucker, D.C.; Prager, W. Soil mechanics and plastic analysis or limit design. Q. Appl. Math. 1952, 10, 157–165. [Google Scholar] [CrossRef] [Green Version]
  51. Kossa, A. Analytical strain solution for the Drucker-Prager elastoplasticity model with linear isotropic hardening. Period. Polytech. Mech. Eng. 2012, 56, 27–31. [Google Scholar] [CrossRef] [Green Version]
  52. Schorn, H.; Rode, U. Numerical simulation of crack propagation from microcracking to fracture. Cem. Concr. Res. 1991, 13, 87–94. [Google Scholar] [CrossRef]
  53. Kopp, S. Proof, Assessment and Quantification of Structural Changes in Concrete with Energetic Parameters (German Title: Nachweis, Beurteilung und Quantifizierung von Gefügeveränderungen im Beton mit Energetischen Kenngrößen). Ph.D. Thesis, Ruhr-Universität Bochum, Bochum, Germany, 1998. [Google Scholar]
Figure 1. Three-point bending beam test: (a) geometry for the experimental setup and boundary value problem and (b) detailed photo of setup.
Figure 1. Three-point bending beam test: (a) geometry for the experimental setup and boundary value problem and (b) detailed photo of setup.
Materials 15 01179 g001
Figure 2. Sketches for the damage indicators in quasi-static tests: (a) CMOD plastic (CMOD pl , i ), (b) initial ( s 0 ) and residual ( s i ) stiffness and (c) damage energy ( W pl , i ).
Figure 2. Sketches for the damage indicators in quasi-static tests: (a) CMOD plastic (CMOD pl , i ), (b) initial ( s 0 ) and residual ( s i ) stiffness and (c) damage energy ( W pl , i ).
Materials 15 01179 g002
Figure 3. Schematic waveform of a signal.
Figure 3. Schematic waveform of a signal.
Materials 15 01179 g003
Figure 4. Exemplary load-CMOD curves: (a) zoomed view of (b) plot for different fibre contents.
Figure 4. Exemplary load-CMOD curves: (a) zoomed view of (b) plot for different fibre contents.
Materials 15 01179 g004
Figure 5. Development of residual stiffness: (a) zoomed view of (b) plot for different fibre contents.
Figure 5. Development of residual stiffness: (a) zoomed view of (b) plot for different fibre contents.
Materials 15 01179 g005
Figure 6. Plastic deformation after unloading: (a) zoomed view of (b) plot for different fibre contents.
Figure 6. Plastic deformation after unloading: (a) zoomed view of (b) plot for different fibre contents.
Materials 15 01179 g006
Figure 7. Absorbed damaging energy W pl , i : (a) zoomed view of (b) plot for different fibre contents.
Figure 7. Absorbed damaging energy W pl , i : (a) zoomed view of (b) plot for different fibre contents.
Materials 15 01179 g007
Figure 8. Phases of deterioration, progression of (a) CMOD, (b) load and (c) registered hits (total and weighted by respective energy) over time.
Figure 8. Phases of deterioration, progression of (a) CMOD, (b) load and (c) registered hits (total and weighted by respective energy) over time.
Materials 15 01179 g008
Figure 9. Three-point bending beam test at low cycle for fiber-reinforced HPC: distribution of stresses σ x (in GPa) in horizontal direction (x-axis), the equivalent plastic strain α HPC for HPC phase, the phase-field parameter q for HPC phase in (a,c,e) at CMOD = 0.018 mm and (b,d,f) at CMOD = 0.108 mm, respectively.
Figure 9. Three-point bending beam test at low cycle for fiber-reinforced HPC: distribution of stresses σ x (in GPa) in horizontal direction (x-axis), the equivalent plastic strain α HPC for HPC phase, the phase-field parameter q for HPC phase in (a,c,e) at CMOD = 0.018 mm and (b,d,f) at CMOD = 0.108 mm, respectively.
Materials 15 01179 g009
Figure 10. Three-point bending beam test for fiber-reinforced HPC: (a) load-CMOD diagramm for experimental data and simulations using (b) five preferred orientations between the angles of −10° to 10° and (c) 24 preferred orientations between the angles of −90° to 90°. Comparison of experimental and simulated (d) load-CMOD diagram and (e) calculated and interpolated residual stiffness-CMOD diagram.
Figure 10. Three-point bending beam test for fiber-reinforced HPC: (a) load-CMOD diagramm for experimental data and simulations using (b) five preferred orientations between the angles of −10° to 10° and (c) 24 preferred orientations between the angles of −90° to 90°. Comparison of experimental and simulated (d) load-CMOD diagram and (e) calculated and interpolated residual stiffness-CMOD diagram.
Materials 15 01179 g010
Table 1. Composition of the used concrete mixtures, data from [34].
Table 1. Composition of the used concrete mixtures, data from [34].
   IngredientsQuantity in kg/m 3 , HPC
   cement: CEM I 52,5R500
   fine sand75
   sand 0/2850
   basalt 2/5350
   basalt 5/8570
   silica fume570
   superplasticizer5
   stabilizer3
   water176
   steel fibres0/23/57/115
   number of specimens3/3/3/3
Table 2. Mechanical properties of HPC, adapted from [34], compare [44,45].
Table 2. Mechanical properties of HPC, adapted from [34], compare [44,45].
   PropertyMean ValueUnit
   tensile strength f t 5.7MPa
   compressive strength f c 112MPa
   Young’s modulus E39.976GPa
   Poisson’s ratio μ 0.192
Table 3. Settings for the measurements of acoustic emission.
Table 3. Settings for the measurements of acoustic emission.
   ParameterValue
   duration discrimination time0.4 ms
   rearm time0.4 ms
   threshold30.1 dB
   filter95 kHz–850 kHz
   gain34 dB
Table 4. Measured mechanical properties, see Table 2, and calibrated material parameters for HPC.
Table 4. Measured mechanical properties, see Table 2, and calibrated material parameters for HPC.
E HPC ν HPC f t f c ψ t c , HPC ψ c c , HPC y 0 HPC β p β n h HPC l ζ E F v F y 0 F h F
GPaMPaMPaMPaMPaMPammGPaMPaMPa
HPC39.9760.1925.71124.2 × 10 4 0.136.20.50.1213,0001412100.003660130
Table 5. Calibrated values for the interpolation of degradation functions for HPC in tension and compression.
Table 5. Calibrated values for the interpolation of degradation functions for HPC in tension and compression.
q00.10.20.30.40.50.60.70.80.91
g + ( q ) 11110.930.850.740.60.390.180.003
g ( q ) 10.8520.7140.5850.4640.3530.2520.1640.0890.030.003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gebuhr, G.; Pise, M.; Anders, S.; Brands, D.; Schröder, J. Damage Evolution of Steel Fibre-Reinforced High-Performance Concrete in Low-Cycle Flexural Fatigue: Numerical Modeling and Experimental Validation. Materials 2022, 15, 1179. https://doi.org/10.3390/ma15031179

AMA Style

Gebuhr G, Pise M, Anders S, Brands D, Schröder J. Damage Evolution of Steel Fibre-Reinforced High-Performance Concrete in Low-Cycle Flexural Fatigue: Numerical Modeling and Experimental Validation. Materials. 2022; 15(3):1179. https://doi.org/10.3390/ma15031179

Chicago/Turabian Style

Gebuhr, Gregor, Mangesh Pise, Steffen Anders, Dominik Brands, and Jörg Schröder. 2022. "Damage Evolution of Steel Fibre-Reinforced High-Performance Concrete in Low-Cycle Flexural Fatigue: Numerical Modeling and Experimental Validation" Materials 15, no. 3: 1179. https://doi.org/10.3390/ma15031179

APA Style

Gebuhr, G., Pise, M., Anders, S., Brands, D., & Schröder, J. (2022). Damage Evolution of Steel Fibre-Reinforced High-Performance Concrete in Low-Cycle Flexural Fatigue: Numerical Modeling and Experimental Validation. Materials, 15(3), 1179. https://doi.org/10.3390/ma15031179

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