Next Article in Journal
CuO Bionanocomposite with Enhanced Stability and Antibacterial Activity against Extended-Spectrum Beta-Lactamase Strains
Next Article in Special Issue
Evaluation of the Behavior of Carbon Short Fiber Reinforced Concrete (CSFRC) Based on a Multi-Sensory Experimental Investigation and a Numerical Multiscale Approach
Previous Article in Journal
Effects of TiAl Alloy as a Binder on Cubic Boron Nitride Composites
Previous Article in Special Issue
Increasing the Fatigue Resistance of Strain-Hardening Cement-Based Composites (SHCC) by Experimental-Virtual Multi-Scale Material Design
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrete Element Modeling and Electron Microscopy Investigation of Fatigue-Induced Microstructural Changes in Ultra-High-Performance Concrete

1
Institute of Materials, Physics and Chemistry of Buildings, Hamburg University of Technology, 21073 Hamburg, Germany
2
Institute of Solids Process Engineering and Particle Technology, Hamburg University of Technology, 21073 Hamburg, Germany
3
Electron Microscopy Unit, Hamburg University of Technology, 21073 Hamburg, Germany
*
Author to whom correspondence should be addressed.
Materials 2021, 14(21), 6337; https://doi.org/10.3390/ma14216337
Submission received: 27 September 2021 / Revised: 15 October 2021 / Accepted: 20 October 2021 / Published: 23 October 2021
(This article belongs to the Special Issue Cyclic Deterioration of Concrete)

Abstract

:
In view of the growing demand for sustainable and lightweight concrete structures, the use of ultra-high-performance concrete (UHPC) is becoming increasingly important. However, fatigue loads occur more frequently in nature than static loads. Despite the impressive mechanical properties of UHPC, a reduced tolerance for cyclic loading is known. For this reason, our paper deals with experimental and numerical investigations regarding the main causes for crack initiation on the meso, micro, and nanoscale. After mechanical fatigue tests, we use both scanning (SEM) and transmission electron microscopy (TEM) to characterize microstructural changes. A new rheological model was developed to apply those changes to the mesoscopic scale. The origins of fatigue damaging can be traced back to a transformation of nanoscale ettringite, resulting in a densification of the surrounding binder matrix. Additionally, a higher content of unhydrated cement clinker in the matrix benefits fatigue resistance. On the mesoscale, stress peaks around aggregate grains expand into the surrounding binder with increasing load cycles and lead to higher degradation.

1. Introduction

Fatigue-induced damage of cementitious structures, similar to concrete foundations of wind turbines or concrete bridges, is more common than static loading in nature. When these structures are subjected to cyclic loading, plastic deformation is observed even though the maximum load level is much smaller than the short-term strength under quasi-static loading conditions. Due to the more frequent use of innovative materials such as ultra-high-performance concrete (UHPC), the focus of research is shifting toward the main crack-initiating factors. UHPC is an almost ideally brittle material with a very dense and high-strength structure. The outstanding mechanical properties are achieved by an optimization on the basis of the maximum packing density theory [1,2] and a very low water to cement ratio (w/c). A high amount of cement, fillers (quartz or limestone powder), and pozzolanic additives such as silica fume lead to the optimized particle size distribution. Due to the low w/c ratio, the degree of cement hydration in UHPC is as low as 30%, so much of the unhydrated cement also acts as physical filler. Furthermore, the pozzolanic reactivity of silica—the additional formation of calcium silicate hydrate (C-S-H)—and its small particle size support a high-strength interfacial transition zone (ITZ) and lead to a lower porosity within the binder [3,4]. As a result, crack propagation from pores or a weaker ITZ cannot be assumed, as it is the case for plain concrete. Recently, investigations regarding crack initiation have concentrated on localizing inhomogeneities, particularly in the case of fatigue loading.
Both scanning (SEM) and transmission electron microscopy (TEM) have been applied to cementitious materials for four decades [5,6,7,8,9], but so far, few if any reports focused on the effects of fatigue on the microstructure of the binder matrix. We aim to investigate these effects, down to their origins on the nanoscale, using SEM as well as scanning TEM on lamella samples prepared using focused ion beam (FIB) techniques. In past research, energy-dispersive X-ray spectroscopy (EDS) applied to cementitious materials has mostly served to determine the local chemical composition of the binder matrix [10,11,12]. However, we additionally employ EDS to detect and identify objects on the nanoscale.
For our numerical investigations, the bonded particle model (BPM) was used. The BPM is an extension of the soft-sphere discrete element method (DEM), where particles can be connected by solid bonds [13,14]. This enables the modeling of porous [15] or coherent materials structures consisting of single or various components. The overall behavior of the material is described by constitutive models for particle interactions and by solid bond models. The possibility to change or remove individual bonds during the simulation allows for the reproduction of small-scale changes or fracture of the material. BPM has been increasingly used by research groups for the numerical analysis of concrete in 2D and 3D simulations during static loading [16,17,18,19,20,21,22], but none of them have applied the BPM to UHPC. Furthermore, only a few models have been developed for cyclic loading [23,24,25].
The main purpose of this article is the illustration of a beneficial application of various experimental and numerical methods to calibrate and concretize a newly derived model approach for fatigue simulations of UHPC. We aim to identify vulnerable elements of the concrete microstructure and to contribute to the development of cementitious building materials with improved fatigue resistance.
A general characterization of components and materials is given in Section 2. Results of SEM and TEM investigations of fatigue-damaged samples are presented in Section 3 and Section 4, respectively. Section 5 introduces the numerical background and proposes a new rheological fatigue model, while Section 6 presents results of fatigue simulations. A short summary and main conclusions are given in Section 7.

2. Materials Characterization

Specimens of UHPC and pure binder were prepared from components given in Table 1. Oxide compositions for cement and silica fume are given in the Supplementary Data (Table S1). Specimens were stripped after 24 h and then stored in water until a minimum age of 56 d, since the compressive strength of UHPC still increases considerably after a standard hydration time of 28 d [26]. This procedure ensures a high degree of hydration and a nearly constant compressive strength. Afterwards, the samples were sawn and ground to cylindrical dimensions of h/d = 180 mm/60 mm.
To investigate the mechanical behavior of UHPC and binder specimens, uniaxial static and cyclic compression tests were carried out. Binder samples have the same composition as the UHPC, excluding the quartz sand.
Uniaxial static compression tests were carried out strain-controlled with a load speed of 0.2 mm/min on three specimens of both materials. Three high-precision inductive displacement transducers were used to capture the longitudinal strain during testing with a sampling rate of 30 Hz. The mean values of the mechanical properties are shown in Table 2. Both the strength and stiffness of the binder samples are lower than for UHPC samples, while the breakage strain is higher. This is related to the missing aggregate grains in the mixture. Nevertheless, the main contributing factor to the high static strength of UHPC and binder is the optimized packing density and low w/c ratio, while the aggregate properties play a secondary role.
After static tests, uniaxial cyclic compression tests were carried out force-controlled with a triangular load regime of 1 Hz. The upper and lower load level ( S u , S l ) corresponded to 80% and 5% of the mean value of the static strength, while the longitudinal strain was measured by three high-precision inductive displacement transducers with a sampling rate of 30 Hz. Figure 1a shows typical stress–strain curves during fatigue testing for both compositions under the mentioned regime, while Figure 1b contains the lifetime values (single values, mean values, and max/min).
The induced fatigue damage leads to an increase of strain and decrease of stiffness for both compositions (Figure 1a). As it is often the case with brittle materials, no failure announcement could be observed, while an abrupt breakage occurred without any post-breakage behavior. Generally, fatigue tests are always subjected to a high scattering caused by material heterogeneity. Nevertheless, a higher fatigue resistance could be observed for binder specimens (Figure 1b).

3. Correlation of Unhydrated Cement Clinker Content and Specimen Lifetime

3.1. Sample Preparation

To find a possible correlation between the lifetime of a specimen and its content of unhydrated cement clinker, we performed SEM investigations and image analysis on sufficiently large sample surfaces to ensure statistical reliability. A similar approach to analyze the hydration degree from segmented images to support numerical models has mostly been applied to μCT data [27], but no reports on the application to SEM data are known so far.
Specimen pieces with a surface of approximately 1 to 2 cm² were cut from suitable fragments of after-failure samples using a low-speed diamond blade saw. The surface was finely ground using diamond-coated sandpaper of four different grits down to 8.4 µm and water cooling. Then, these samples were glued to standard aluminum SEM carriers and coated with a thin electrically conductive layer using carbon evaporation.

3.2. SEM Investigation and Image Processing

Material or atomic number (Z) contrast imaging of backscattered electrons (BSE) with a concentric ring detector was employed to easily discern between aggregate, binder matrix, and unhydrated cement clinker. To survey and analyze large surface areas, individual images with a horizontal field width of ca. 1 mm and a pixel size of below 1 µm were acquired in a serpentine pattern, overlapping by 20% in both the x and y directions. The individual images were fused using a FIJI/ImageJ plugin [28] to create a single large-area image of each sample. These were cropped to exclude specimen edges, covering 8000×8000 pixels (ca. 5.6 × 5.6 mm²) in total.
To improve the image quality and reliability of analysis, several image processing steps were performed:
  • Non-local means denoising (σ = 20) [29,30] to reduce noise, especially in homogeneous parts of the specimen surface, such as aggregate grains;
  • Band-pass filter with lower and upper limits of 3 and 5000 pixels, respectively, to suppress brightness gradients caused by an imperfectly plan-parallel specimen surface;
  • K-means clustering [31,32,33] to segment the image into four individual clusters, corresponding to, from the darkest to brightest: cracks and pore edges, aggregate grains, hydrated binder matrix, and unhydrated cement clinker grains;
  • Gray-level thresholding to select desired components and particle analysis (10 pixels or larger) to exclude negligibly small particles stemming from residual image noise.
Exemplary large-area SEM images, corresponding gray-level histograms, and higher-magnification insets of the marked area are shown in Figure 2.

3.3. Area Calculation and Analysis

The respective area fractions of the individual components in the images were considered a sufficiently accurate approximation of their content in the specimen volume. The relative content x of unhydrated cement clinker in the binder matrix was calculated as a ratio of unhydrated cement clinker area A clinker (taken directly from particle analysis) and binder area A binder (total image area A without aggregate area A aggregate from particle analysis):
x = A clinker A binder = A clinker A A aggregate
We find a positive correlation between the relative content of unhydrated cement clinker in a specimen and its eventual lifetime, the effect being more pronounced for UHPC samples due to their lower total binder content. The corresponding graph is shown in Figure 3. A high content of unhydrated cement clinker, which can be considered a high-strength ceramic material, lowers local w/c ratios and in turn increases local strength.

4. Characterization of Structural Changes on the Nanoscale

The lasting strain introduced with every single load cycle must be connected to irreversible changes in the nanoscale structure of the material. In a recent publication [34], we used TEM for this purpose and observed needle-shaped objects forming in the binder matrix of UHPC already after the early onset of fatigue damage. These objects retain their projected size of ca. 200 nm in length and ca. 30 nm in width but continuously increase in number as fatigue damage progresses.
Due to their increasingly dark contrast with respect to the surrounding material in Z-contrast high-angle annular dark field (HAADF) STEM imaging, we initially concluded these needle-shaped structural changes to be crack precursors, continuously depleting in material while the surrounding matrix is densified.
Upon further investigation, we find that the structural changes (1) either shrink or grow in size upon energy uptake from electron irradiation during TEM investigation and (2) preferentially exist in regions of the binder matrix comparatively rich in aluminum and sulfur and with a low silicon content. A selection of HAADF-STEM image sections after increasing irradiation durations is shown in Figure 4; EDS elemental distribution maps and spectra are displayed in Figure 5.
As these structural changes exhibit dark Z contrast in HAADF-STEM imaging, readily change size and shape upon energy uptake, and preferentially exist in regions rich in aluminum and sulfur, we conclude they are connected to locations of a transformation of ettringite.
Ettringite is a mineral with the nominal chemical formula Ca6Al2(SO4)3(OH)12 · 26 H2O. It is rarely found in nature but plays a crucial role in the early stage of the setting of concrete, forming as a result of the reaction of calcium aluminate and calcium sulfate from Portland cement. A delayed formation of ettringite in hardened concrete can be detrimental to the durability and longevity of the material [35]. Due to its high water content of 26 molecules per formula unit, its average atomic number is considerably lower than for other components of UHPC, resulting in a darker image contrast in HAADF-STEM imaging.
For the same reason, its density is low for a crystalline solid (1.77 g/cm³) and its molar volume is up to four times higher than that of C-S-H phases that constitute the majority of the binder matrix. Ettringite is known to readily recrystallize upon changes in environmental parameters, forming coarse, rigid crystals. It is commonly detectable on fracture surfaces and cracks in damaged concrete samples [36]. A recurring dissolution and crystallization can exert pressure on the surrounding material, leading to solidification and, eventually, the formation of cracks on the nanoscale.
Sample regions containing nanoscale cracks found in after-failure samples exhibit especially high contents of aluminum and sulfur and a very low silicon content, underlining the connection between fatigue damage and nanoscale ettringite transformation. HAADF-STEM images as well as EDS spectra and elemental distribution maps of one such region are displayed in Figure 6.

5. Numerical Background

5.1. Discrete Element Method

The discrete element method is mainly used to discretize granular materials as an assemblage of primary particles, having their own physical properties and motions. Particles interact with each other and with defined boundary walls, with forces between those interactions being governed by contact models. The motion of particles is described by Newton’s second law, and the acceleration is time-integrated to find the new velocities and positions. Initially, DEM was developed for the modeling of systems consisting of spherical particles [37], but nowadays, several extensions of DEM are used for non-spherical particles, the modeling of multiphase flows [38], or chemical reactions. The simulations presented in this contribution have been performed in the open-source DEM simulation framework MUSEN [39]. The latest version of the code can be downloaded from [40]. The calculation of forces and integration of motion is performed on a graphic processor unit (GPU), which reduces the computation time significantly. The current study has been performed on a computer equipped with Intel Xeon Gold 5,118 CPU and NVIDIA A100 GPU.

5.2. Bonded-Particle Model

The bonded particle model is an extension of the DEM, where each pair of the mentioned primary particles can be connected by a solid bond [13,14]. These bonds make it possible to describe a coherent materials structure consisting of single or various components. The overall behavior of the material is described by the interparticle contacts and additional bond models. The possibility to remove individual solid bonds during the simulation allows explicitly simulating the inner fracturing of the material. In this work, a soft-sphere DEM formulation has been used where occurring interparticle or particle–wall overlaps can be seen as local material deformations. Overall, three different types of interactions have been applied in the BPM model:
  • Interparticle contacts: calculation of forces between all primary particles even if they are connected with solid bonds;
  • Particle–wall contacts: interaction between compression walls and primary particles;
  • Solid bonds: forces, moments, and torques acting in single bonds between two particles.

5.3. Contact Models

Due to the almost ideally elastic behavior of UHPC and binder, models with an elastic base were chosen for all three mentioned types of interactions. For the interparticle model, a linear elastic formulation was used, while the Hertz–Mindlin model was chosen for the particle–wall contact. A detailed description of these models can be found in [41,42,43]. With respect to the rheological bond model, a newly derived approach was implemented, which uses elastic and plastic deformations as well as fatigue damaging. Basic equations for the purely elastic sections, which also find application in our model, can be found in [13]. However, the main focus of the following description is on the plastic and fatigue damaging of a bond.
Since concrete is known to fail due to internal transverse tensile stresses, the main damaging mechanisms were integrated into the tensile response of the model. Additionally, to capture fatigue degradation from shear stresses, a fatigue damage was also integrated into shear mode. Note that bonds do not break or undergo any damage during pure compression. Figure 7 shows the rheological models for all three load types.
Generally, the stress acting on a bond is a result of the interaction between the corresponding particles. Solid bonds are modeled as ideally cylindrical objects with an initial diameter d init , which cannot exceed the minimal diameter of the connected particles. Depending on the particle movement or rather the relative velocity v rel , different stresses are acting. In normal direction (compression and tension), the increment of strain Δ ε com , ten and stress Δ σ com , ten is calculated by Equations (2) and (3), while the sum of the increments yields to the current stress value σ current :
  Δ ε com , ten = v rel , n Δ t L init
Δ σ com , ten =   Δ ε     E b  
where v rel , n is the relative velocity of the particles in normal direction, Δ t is the time step, E b is the bond Young’s modulus, and L init is the initial bond length. The used indices for the stress and strain increments represent the respective load type (compression, tension, shear) and the load direction (normal or shear). Once the tensile yield limit σ ten , yield is exceeded, softening occurs, and the bond starts to deform plastically. Then, the increment of the stress Δ σ ten is calculated using Equation (4), where α is the weakening factor between the elastic and plastic section of the tensile bond response:
  Δ σ ten =   Δ ε ten E b α .
In the case of fatigue loading of UHPC, crack propagation occurs on different scales and does not inevitably lead to plastic deformation on the mesoscale. For this reason, it is not sufficient to reproduce fatigue degradation only by a plastic model formulation. Instead, crack growth takes place on all small scales (meso, micro, and nanoscale) simultaneously, which has to be taken into account. To capture this on the mesoscale, a fatigue damage criterion was implemented, where the bond diameter is reduced as soon as fatigue damage occurs. Using this approach, the mentioned micro and nanocracks can be represented without having a bond to break directly and be removed from the simulation. Figure 8 illustrates the connection between small-scale damage and its implementation on the mesoscale bond.
Fatigue damage occurs in the elastic range and thus is separated from plastic deformation, which is important in the case of fatigue mechanisms. Therefore, the bond diameter is determined by the stresses acting on two consecutive simulation time steps. The current and previous stress values σ ten , cur   and σ ten , prev are analyzed and compared to the critical value σ ten , fatigue , which defines the fatigue level in the elastic section:
σ ten , prev > σ ten , fatigue     ;     σ ten , prev < σ ten , yield  
σ ten , cur < σ ten , fatigue .
If the conditions (5) and (6) are fulfilled, then the bond diameter is reduced to a new diameter according to Equation (7), and very small destructions in the supporting material structure are taken into account, as also shown in Figure 8. The larger the value for the damaging factor δ n   , the smaller the represented micro or nanocrack in the bond:
d new = d init   ( 1 10 δ n ) .
To capture fatigue damage under tension and shear load, a similar criterion was additionally implemented in the shear response of the model. Therefore, the bond diameter is linked with the current and previous stress value τ s , cur   and τ s , prev , while τ s , fatigue defines the critical fatigue level:
τ s , prev > τ s , fatigue     ;     τ s , prev < τ max
τ s , current < τ s , fatigue .
If conditions (8) and (9) are fulfilled, then the bond diameter is corrected to a new diameter according to Equation (10). The larger the value for the damaging factor δ s , the smaller the represented micro or nanocrack in the bond:
d new = d init   ( 1 10 δ S ) .
In order to simulate material destruction, total stresses including additional bending and torsional moments are analyzed and compared with specified material properties for normal-tension strength σ ten , max as well as shear strength τ max . If one of the conditions in Equations (11) or (12) is fulfilled, the bond breaks and is removed from the simulation:
τ s , current + M t , b   d new 2 I T , b > τ max
      σ current + M n , b   d new 2 J b > σ t ,   max
where M t , b   and M n , b   are the torsional and bending moments acting on the bond, J b   is the moment of inertia, and I T , b   is the torsional moment of inertia.

5.4. Structural UHPC Model and Calibration of Parameters

With respect to the numerical investigation of the fatigue behavior of UHPC, a structural three-phase UHPC model was used. The model consists of two types of particles (binder and aggregate) and three types of bonds (binder–binder, binder–aggregate, and aggregate–aggregate). These three types of bonds represent the main three phases of concrete (binder, interfacial transition zone, and aggregate) and allow to distinguish and analyze different sections of the modeled UHPC. To consider a realistic non-spherical shape of aggregate, reconstructed 3D models from X-ray computer tomography (µCT) were used for the structural model. The bond formulations presented in Section 5.2 were implemented for the binder and ITZ bonds, while a purely elastic model was used for the aggregate bonds. A detailed procedure of the used algorithms on which the structure generation is based can be found in [41], and only the most significant points will be mentioned here as follows.
Caused by the optimized packing density of the used UHPC composition, a potentially weaker ITZ was not observed in SEM investigations [41]. For this reason, parameters of ITZ and binder bonds were assumed to be equal. This reduced the calibration effort to the remaining phases of binder and aggregate. With respect to the calibration of binder parameters, results of uniaxial compression tests from Table 2 were used for static binder simulations. As cracks were not penetrating the single grains in experiments, only the correct stiffness of aggregates was calibrated, and the bond strength was set to infinity. For the sake of simplicity, the numerical calibration of aggregate parameters was carried out on cylindrical specimens under uniaxial compression. Generally, the model parameters have been calibrated using a gradient-based trial and error principle. The simulations were iteratively repeated using different parameters to minimize the discrepancy between experimental and numerical results. The slope of the stress–strain diagram was used to adjust Young’s modulus of the bonds, while the bond strength parameters were adjusted based on the fracture strength of the macroscopic sample. To estimate an initial parameter set, a linearization-based technique has been used. Here, the strategy based on the direct stiffness method has been applied. The detailed description of this technique is provided in [44].
Finally, using the calibrated parameters for single components such as binder, ITZ, and aggregate, uniaxial static compression of UHPC was simulated and validated with experimental data; see Table 3. Note that the parameters of interparticle friction μ sl have not been included in the calibration strategy and were assumed to be equal between all particles. Figure 9 illustrates the cross-section detail of the structural three-phase model, while Table 4 contains the calibrated model parameters.

6. Fatigue Simulations

Uniaxial Fatigue Simulations

After the calibration of model parameters for single components such as binder, ITZ, and aggregate, uniaxial cyclic compression tests were simulated. Three separately generated samples of UHPC and binder were calculated using the same parameter set and validated with experimental results, respectively. The upper and lower load levels ( S u , S l ) corresponded to 80% and 5% of the static strength of each simulated cylinder. Figure 10a contains two exemplary hysteresis loops obtained from cyclic UHPC and binder simulations, while Figure 10b shows a comparison between the experimental range of concrete fatigue tests and the mean strain-lifetime (S-N) curve of fatigue UHPC simulations.
By comparing Figure 1a and Figure 10a, it can be clearly noticed that our fatigue simulations are in good agreement with the experimental stress–strain hysteresis loops. The simulations show the same macroscopic characteristics, e.g., the increase of strain and decrease of stiffness until fatigue breakage. Additionally, the mean S-N curve obtained from simulations shows a typical S-curve within the range of the experimental data.
Generally, fatigue simulations were carried out for different damaging factors ( δ n ,   δ s   = 3, 6, 8; Table 5). As already mentioned in Section 5.3, the larger the value for δ n and δ s , the smaller the represented micro or nanocrack induced in the bond. However, all simulations showed a higher resistance of binder compared with UHPC, which has been equally shown in fatigue experiments. Furthermore, the different damaging factors can be linked to experimental SEM investigations, where a higher content of unhydrated cement clinker leads to higher fatigue resistance (see Section 3.3). The increased presence of unhydrated cement leads to a local increase of packing density; consequently, transverse strain is blocked through the cement particles, and the induced damage is smaller.
In order to localize the main crack initiation in the investigated UHPC simulations, average stresses on binder particles (at the upper load level) were analyzed at different points of fatigue (N = 1; N/Nf = 0.3, 0.6, 0.9). In the first cycle, particles subjected to high stresses were localized between close-knit grains, while the stress exerted on particles is much lower in regions of binder containing no aggregate (Figure 11; N = 1). Aggregate grains absorb inner stresses and contribute to an increase of macroscopic strength. However, discontinuities caused by the high difference of stiffness induce stress peaks and lead to a higher local degradation of the material.
Furthermore, the stress peaks around the aggregate grains propagate into the surrounding matrix with increasing lifetime (dotted circles), while the less stressed regions without grains do not show any significant changes (squares; Figure 11).
Additionally, average stress distributions on binder particles (in the loading direction) were analyzed at the same fatigue points and load levels (Figure 12a). For higher fatigue stages, the distributions broaden and their maxima drop. This indicates a more inhomogeneous distribution of stress among particles, resulting in the propagation of stress peaks into the surrounding matrix.
Furthermore, the aforementioned stress peaks result in higher stresses within the ITZ bonds and lead to a higher amount of defects during fatigue (Figure 12b). This amount of defects is the ratio of broken bonds compared to its initial number. It can be clearly noticed that the number of broken bonds is dominated by the failure within the ITZ. The largest increase of defects occurs in the first cycles. Afterwards, both curves follow an almost parallel course until breakage. Therefore, the first cycles of fatigue loading play an essential role, as damage and rearrangement within the material structure already occur. Additionally, the binder with its high packing density does not have any buffering effect to reduce stress peaks within the structure, so the damage propagates from the ITZ into the matrix.

7. Conclusions

This paper presents an experimental and numerical analysis of the fatigue behavior of UHPC using both scanning and transmission electron microscopy as well as a newly derived bond formulation. Structural inhomogeneities such as local variations in stiffness or ettringite transformation are possible causes for crack initiation. Similar to these inhomogeneities, deterioration can occur on the meso, micro, and nanoscale.
Based on our results, the following main conclusions may be drawn:
  • Microscopic investigations proved that fatigue damaging also occurs on the nanoscale as a transformation of nanoscale ettringite. Therefore, to reproduce the correct mechanisms of degradation during cyclic simulations, nanoscale changes must be taken into account even in the case of mesoscale modeling. Using sulfate-resistant cements (SR) containing little or no calcium aluminate could possibly prevent a delayed formation of ettringite and slow down the fatigue process.
  • A higher relative content of unhydrated cement clinker in the binder matrix results in higher local strength and consequently leads to a longer specimen lifetime during cyclic loading.
  • Stress peaks caused by a high difference of local stiffness manifest in a higher degradation of the binder directly surrounding aggregate grains during fatigue loading. This promotes a crack initiation in the ITZ and favors a propagation in the surrounding matrix. Using less stiff and even more fine-grained aggregate could minimize local stress peaks and increase fatigue lifetime.
For future investigations, UHPC containing fibers or coarse-grained aggregate could be used to bridge the gap between mesoscale and macroscale in both our experimental and numerical studies.

Supplementary Materials

The following supplementary data are available online at https://www.mdpi.com/article/10.3390/ma14216337/s1. Table S1: Oxide composition of cement and silica fume.

Author Contributions

Conceptualization, S.R. and G.S.; methodology, S.R. and G.S.; software, M.D.; validation, M.D., M.R. and F.S.-D.; formal analysis, S.R., G.S., M.D., M.R. and F.S.-D.; investigation, S.R. and G.S.; data curation, S.R. and G.S.; writing—original draft preparation, S.R. and G.S.; writing—review and editing, S.R., G.S., M.D., M.R. and F.S.-D.; visualization, S.R. and G.S.; supervision, M.D., M.R. and F.S.-D.; project administration, M.D., M.R. and F.S.-D.; funding acquisition, M.D., M.R. and F.S.-D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), grant number 353408149, as part of the priority program “SPP2020: Cyclic damage processes in high-performance concretes in the Experimental Virtual Lab”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We would like to thank our colleagues and partners in the priority program SPP2020 for fruitful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arora, A.; Aguayo, M.; Hansen, H.; Castro, C.; Federspiel, E.; Mobasher, B.; Neithalath, N. Microstructural packing-and rheology-based binder selection and characterization for ultra-high performance concrete (UHPC). Cem. Concr. Res. 2018, 103, 179–190. [Google Scholar] [CrossRef]
  2. Wang, X.; Yu, R.; Song, Q.; Shui, Z.; Liu, Z.; Wu, S.; Hou, D. Optimized design of ultra-high performance concrete (UHPC) with a high wet packing density. Cem. Concr. Res. 2019, 126, 105921. [Google Scholar] [CrossRef]
  3. Alkaysi, M.; El-Tawil, S.; Liu, Z.; Hansen, W. Effects of silica powder and cement type on durability of ultra high performance concrete (UHPC). Cem. Concr. Compos. 2016, 66, 47–56. [Google Scholar] [CrossRef]
  4. Abdulkareem, O.M.; Fraj, A.B.; Bouasker, M.; Khelidj, A. Mixture design and early age investigations of more sustainable UHPC. Constr. Build. Mater. 2018, 163, 235–246. [Google Scholar] [CrossRef] [Green Version]
  5. Groves, G.W. TEM studies of cement hydration. Res. Soc. Symp. Proc. 1987, 85, 3–12. [Google Scholar] [CrossRef]
  6. Richardson, I.G.; Groves, G.W. Microstructure and microanalysis of hardened ordinary Portland cement pastes. J. Mater. Sci. 1993, 28, 265–277. [Google Scholar] [CrossRef]
  7. Diamond, S. Aspects of concrete porosity revisited. Cem. Concr. Res. 1999, 29, 1181–1188. [Google Scholar] [CrossRef]
  8. Gatty, L.; Bonnamy, S.; Felessoufi, A.; Clinard, C.; Richard, P.; van Damme, H. A transmission electron microscopy study of interfaces and matrix homogeneity in ultra-high-performance cement-based materials. J. Mater. Sci. 2001, 36, 4013–4026. [Google Scholar] [CrossRef]
  9. Wenzel, O.; Schwotzer, M.; Müller, E.; Chakravadhanula, V.S.K.; Scherer, T.; Gerdes, A. Investigating the pore structure of the calcium silicate hydrate phase. Mater. Charact. 2017, 133, 133–137. [Google Scholar] [CrossRef]
  10. Diamond, S. The microstructure of cement paste and concrete—A visual primer. Cem. Concr. Compos. 2004, 26, 919–933. [Google Scholar] [CrossRef]
  11. Wei, F.Y.; Lan, X.H.; Lv, Y.N.; Xu, Z. Effect of Pozzolanic Reaction Products on Alkali-silica Reaction. J. Wuhan Univ. Technol. Mater. Sci. Ed. 2006, 21, 168–171. [Google Scholar] [CrossRef]
  12. Sakalli, Y.; Trettin, R. Investigation of C3S hydration mechanism by transmission electron microscope (TEM) with integrated SuperXTM EDS system. J. Microsc. 2017, 267, 81–88. [Google Scholar] [CrossRef]
  13. Potyondy, D.O.; Cundall, P.A. A bonded-particle model for rock. Int. J. Rock Mech. Min. Sci. 2004, 41, 1329–1364. [Google Scholar] [CrossRef]
  14. Dosta, M.; Dale, S.; Antonyuk, S.; Wassgren, C.; Heinrich, S. Numerical and experimental analysis of influence of granule microstructure on its compression breakage. Powder Technol. 2016, 299, 87–97. [Google Scholar] [CrossRef]
  15. Dosta, M.; Bistreck, K.; Skorych, V.; Schneider, G.A. Mesh-free micromechanical modeling of inverse opal structures. Int. J. Mech. Sci. 2021, 204, 106577. [Google Scholar] [CrossRef]
  16. Cusatis, G.; Mencarelli, A.; Pelessone, D.; Baylot, J. Lattice discrete particle model (LDPM) for failure behavior of concrete. II: Calibration and validation. Cem. Concr. Compos 2011, 33, 891–905. [Google Scholar] [CrossRef]
  17. Donze, F.; Magnier, S.; Daudeville, L.; Mariotti, C.; Davenne, L. Numerical study of compressive behavior of concrete at high strain rates. J. Eng. Mech. 1999, 125, 1154–1163. [Google Scholar] [CrossRef]
  18. Hentz, S.; Donze, F.; Daudeville, L. Identification and validation of a discrete element model for concrete. J. Eng. Mech. 2004, 130, 709–719. [Google Scholar] [CrossRef] [Green Version]
  19. Poinard, C.; Piotrowska, E.; Malecot, Y.; Daudeville, L.; Landis, E. Compression triaxial behavior of concrete: The role of the mesostructure by analysis of X-ray tomographic images. Eur. J. Environ. Civ. Eng. 2012, 16, 115–136. [Google Scholar] [CrossRef]
  20. Kim, H.; Buttlar, W.G. Discrete fracture modeling of asphalt concrete. Int. J. Solids Struct. 2009, 46, 2593–2604. [Google Scholar] [CrossRef] [Green Version]
  21. Nitka, M.; Tejchman, J. Modelling of concrete behavior in uniaxial compression and tension with DEM. Granul. Matter 2015, 17, 145–164. [Google Scholar] [CrossRef]
  22. Skarzynski, L.; Nitka, M.; Tejchman, J. Modelling of concrete fracture at aggregate level using FEM and DEM based on X-ray μCT images of internal structure. Eng. Fract. Mech. 2015, 147, 13–35. [Google Scholar] [CrossRef]
  23. Sinaie, S.; Ngo, T.D.; Nguyen, V.P. A discrete element model of concrete for cyclic loading. Comput. Struct. 2018, 196, 173–185. [Google Scholar] [CrossRef]
  24. Nguyen, N.H.T.; Bui, H.H.; Kodikara, J.; Arooran, S.; Darve, F. A discrete element modelling approach for fatigue damage growth in cemented materials. Int. J. Plast. 2018, 112, 68–88. [Google Scholar] [CrossRef]
  25. Song, Z.; Konietzky, H.; Herbst, M. Three-dimensional particle model based numerical simulation on multi-level compressive cyclic loading of concrete. Constr. Build. Mater. 2019, 225, 661–677. [Google Scholar] [CrossRef]
  26. Deutscher, M.; Tran, N.L.; Scheerer, S. Experimental investigation of load-induced temperature development in UHPC subjected to cyclic loading. In Proceedings of the 6th International fib Congress, Krakow, Poland, 27–29 May 2019. [Google Scholar]
  27. Gan, Y.D.; Zhang, H.Z.; Liang, M.F.; Schlangen, E.; van Breugel, K.; Savija, B. A numerical study of fatigue of hardened cement paste at the microscale. Int. J. Fatigue 2021, 151, 106401. [Google Scholar] [CrossRef]
  28. Preibisch, S.; Saalfeld, S.; Tomancak, P. Globally optimal stitching of tiled 3D microscopic image acquisitions. Bioinformatics 2009, 25, 1463–1465. [Google Scholar] [CrossRef]
  29. Buades, A.; Coll, B.; Morel, J.M. Non-Local Means Denoising. Image Process. Line 2011, 1, 208–212. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, G.H.; Yang, L.; Xiong, W.; Li, Y. An improved non-local means filter for color image denoising. Optik 2018, 173, 157–173. [Google Scholar] [CrossRef]
  31. Lloyd, S.P. Least Squares Quantization in PCM. IEEE Trans. Inf. Theory 1982, 28, 129–137. [Google Scholar] [CrossRef]
  32. Kanungo, T.; Mount, D.M.; Netanyahu, N.S.; Piatko, C.D.; Silverman, R.; Wu, A.Y. An Efficient k-Means Clustering Algorithm: Analysis and Implementation. IEEE Trans. Pattern Anal. Mach. Intell. 2002, 24, 881–892. [Google Scholar] [CrossRef]
  33. Chen, L.; Shan, W.B.; Peng, L. Identification of concrete aggregates using K-means clustering and level set method. Structures 2021, 34, 2069–2076. [Google Scholar] [CrossRef]
  34. Schaan, G.; Rybczynski, S.; Schmidt-Döhl, F.; Ritter, M. Fatigue of Concrete Examined on the Nanoscale. TEM Studies of Fatigue-Induced Changes in the Cement Paste of UHPC. Imag. Microsc. 2020, 2, 25–27. [Google Scholar]
  35. Zhang, Y.L.; Pan, Y.F.; Zhang, D.K. A literature review on delayed ettringite formation: Mechanism, affect factors and suppressing methods. Mag. Concr. Res. 2021, 73, 325–342. [Google Scholar] [CrossRef]
  36. Stark, J.; Bollmann, K. Delayed ettringite formation in concrete. Nord. Concr. Res. Publ. 2000, 23, 4–28. [Google Scholar]
  37. Cundall, P.A.; Strack, O.D.L. A discrete numerical model for granular assemblies. Geotechnique 1979, 29, 47–65. [Google Scholar] [CrossRef]
  38. Kieckhefen, P.; Pietsch, S.; Dosta, M.; Heinrich, S. Possibilities and Limits of Computational Fluid Dynamics–Discrete Element Method Simulations in Process Engineering: A Review of Recent Advancements and Future Trends. Annu. Rev. Chem. Biomol. Eng. 2020, 11, 397–422. [Google Scholar] [CrossRef] [Green Version]
  39. Dosta, M.; Skorych, V. MUSEN: An open-source framework for GPU-accelerated DEM simulations. SoftwareX 2020, 12, 100618. [Google Scholar] [CrossRef]
  40. GitHub—Msolids/Musen. Available online: https://github.com/msolids/musen (accessed on 8 October 2021).
  41. Rybczynski, S.; Dosta, M.; Schaan, G.; Ritter, M.; Schmidt-Döhl, F. Numerical study on the mechanical behavior of ultrahigh performance concrete using a three-phase discrete element model. Struct. Concr. 2020. [Google Scholar] [CrossRef]
  42. Mindlin, R.D.; Deresiewicz, H. Elastic spheres in contact under varying oblique forces. J. Appl. Mech. 1953, 20, 327–344. [Google Scholar] [CrossRef]
  43. Tsuji, Y.; Tanaka, T.; Ishida, T. Lagrangian numerical simulation of plug flow of cohesionsless particles in a horizontal pipe. Powder Technol. 1992, 71, 239–250. [Google Scholar] [CrossRef]
  44. Jarolin, K.; Dosta, M. Linearization-based methods for the calibration of bonded-particle models. Comput. Part. Mech. 2020, 8, 511–523. [Google Scholar] [CrossRef]
Figure 1. Results of experimental fatigue tests on UHPC and binder: (a) Typical stress–strain curves during fatigue; (b) Lifetime until fatigue failure (mean and single values).
Figure 1. Results of experimental fatigue tests on UHPC and binder: (a) Typical stress–strain curves during fatigue; (b) Lifetime until fatigue failure (mean and single values).
Materials 14 06337 g001
Figure 2. Steps of SEM image processing and segmentation: (a) Unprocessed image; (b) After non-local means denoising; (c) After band-pass filtering; (d) After k-means clustering segmentation. Scale bars in large-area images and higher-magnification insets are 1 mm and 200 µm, respectively.
Figure 2. Steps of SEM image processing and segmentation: (a) Unprocessed image; (b) After non-local means denoising; (c) After band-pass filtering; (d) After k-means clustering segmentation. Scale bars in large-area images and higher-magnification insets are 1 mm and 200 µm, respectively.
Materials 14 06337 g002aMaterials 14 06337 g002b
Figure 3. Graph of specimen lifetime vs. their relative content of unhydrated cement clinker showing a positive correlation.
Figure 3. Graph of specimen lifetime vs. their relative content of unhydrated cement clinker showing a positive correlation.
Materials 14 06337 g003
Figure 4. Sections of HAADF-STEM images after increasing durations of electron irradiation showing both shrinking and growing structural changes. From left to right: after 6 s; 30 s; 1 min; 2 min; 5 min; 10 min. Scale bar is 200 nm.
Figure 4. Sections of HAADF-STEM images after increasing durations of electron irradiation showing both shrinking and growing structural changes. From left to right: after 6 s; 30 s; 1 min; 2 min; 5 min; 10 min. Scale bar is 200 nm.
Materials 14 06337 g004
Figure 5. Ettringite transformation in the binder matrix after early-stage fatigue: (a) HAADF-STEM image indicating where spectra were extracted; (b) Corresponding EDS spectra; (c) EDS map of aluminum distribution; (d) EDS map of sulfur distribution. Scale bar is 500 nm.
Figure 5. Ettringite transformation in the binder matrix after early-stage fatigue: (a) HAADF-STEM image indicating where spectra were extracted; (b) Corresponding EDS spectra; (c) EDS map of aluminum distribution; (d) EDS map of sulfur distribution. Scale bar is 500 nm.
Materials 14 06337 g005aMaterials 14 06337 g005b
Figure 6. Sample region showing strong ettringite transformation around a nanoscale crack: (a) HAADF-STEM image indicating where spectra were extracted; (b) Corresponding EDS spectra; (c) EDS map of aluminum distribution; (d) EDS map of sulfur distribution. Scale bar is 500 nm.
Figure 6. Sample region showing strong ettringite transformation around a nanoscale crack: (a) HAADF-STEM image indicating where spectra were extracted; (b) Corresponding EDS spectra; (c) EDS map of aluminum distribution; (d) EDS map of sulfur distribution. Scale bar is 500 nm.
Materials 14 06337 g006
Figure 7. Bond response in different load modes: (a) Compression; (b) Tension; (c) Shear.
Figure 7. Bond response in different load modes: (a) Compression; (b) Tension; (c) Shear.
Materials 14 06337 g007
Figure 8. Bond damage caused by a small-scale crack and its implementation: (a) Undamaged bond connection of two primary particles; (b) A fatigue crack manifests in a reduction of the load-bearing cross-cut of a meso bond; (c) Representation of a fatigue crack in the model by a reduced bond diameter.
Figure 8. Bond damage caused by a small-scale crack and its implementation: (a) Undamaged bond connection of two primary particles; (b) A fatigue crack manifests in a reduction of the load-bearing cross-cut of a meso bond; (c) Representation of a fatigue crack in the model by a reduced bond diameter.
Materials 14 06337 g008
Figure 9. Cross-section detail of the structural three-phase UHPC model: (a) Reconstructed 3D models from µCT of aggregate grains; (b) Particle view of binder and aggregate without bonds; (c) Cross-section showing all particles and bonds.
Figure 9. Cross-section detail of the structural three-phase UHPC model: (a) Reconstructed 3D models from µCT of aggregate grains; (b) Particle view of binder and aggregate without bonds; (c) Cross-section showing all particles and bonds.
Materials 14 06337 g009
Figure 10. Macroscopic fatigue simulation results of UHPC using the proposed bond formulation with δ n ,   δ s = 6 : (a) Hysteresis loops obtained from fatigue simulations for UHPC and binder; (b) Comparison of mean S-N curve from UHPC simulations with experimental range.
Figure 10. Macroscopic fatigue simulation results of UHPC using the proposed bond formulation with δ n ,   δ s = 6 : (a) Hysteresis loops obtained from fatigue simulations for UHPC and binder; (b) Comparison of mean S-N curve from UHPC simulations with experimental range.
Materials 14 06337 g010
Figure 11. Visualization of stresses on primary binder particles at the upper load level S u and different fatigue stages; stresses are analyzed in the loading direction.
Figure 11. Visualization of stresses on primary binder particles at the upper load level S u and different fatigue stages; stresses are analyzed in the loading direction.
Materials 14 06337 g011
Figure 12. Evaluation of particle stresses and bond failure of the mesoscale model: (a) Particle stress distributions of binder particles at different fatigue stages; (b) Comparison of ITZ and binder bond failure.
Figure 12. Evaluation of particle stresses and bond failure of the mesoscale model: (a) Particle stress distributions of binder particles at different fatigue stages; (b) Comparison of ITZ and binder bond failure.
Materials 14 06337 g012
Table 1. Composition for UHPC and binder and physical properties of individual components.
Table 1. Composition for UHPC and binder and physical properties of individual components.
ComponentUHPCBinderContent [kg/m3]Vol %Specific Gravity [g/cm3]D50 [μm]
CEM I (52,5)++795.4026.162.75–3.2015.49
Silica fume++168.607.802.20*
Quartz powder++198.407.652.6525.07
Quartz sand (0.125/0.5)+971.0036.582.65266.98
Superplasticizer (PCE)++27.802.701.09-
Water++187.9019.111.00-
* No D50 is available for the used silica fume. Primary grain size is given with 0.1–0.3 μm.
Table 2. Mechanical properties of UHPC and binder obtained from experimental tests.
Table 2. Mechanical properties of UHPC and binder obtained from experimental tests.
MaterialStrength [MPa]Young’s Modulus [GPa]Breakage Strain [‰]
UHPC193.34 ± 3.7746.67 ± 0.944.73 ± 0.26
Binder172.00 ± 1.6336.00 ± 1.635.33 ± 0.29
Table 3. Simulation results obtained from uniaxial compression tests.
Table 3. Simulation results obtained from uniaxial compression tests.
MaterialStrength [MPa]Young’s Modulus [GPa]Breakage Strain [‰]
UHPC198.12 ± 7.9945.49 ± 1.144.95 ± 1.31
Binder169 ± 1.6736.87 ± 0.325.21 ± 0.24
Table 4. Main material parameters of the BPM model after calibration with static experiments.
Table 4. Main material parameters of the BPM model after calibration with static experiments.
ParameterBinderITZAggregate
Particle and bond diameter [mm]0.480.480.48
Number of particles [-]208,125-118,687
Number of solid bonds [-]2,053,425554,0821,031,387
Simulations time step [s]4.00 × 10−94.00 × 10−94.00 × 10−9
Particle Young’s modulus [GPa]87.26-30.00
Bond Young’s modulus [GPa]10.0010.0040.00
Normal and shear stiffness of particles
k n , k s [N/m]
8.47 × 106-8.47 × 106
Bond normal strength [GPa]1.081.08
Bond shear strength [GPa]2.812.81
Weakening factor α [-]0.350.35-
Yield limit in tension σ ten ,   yield [MPa]84.0084.00-
Ratio σ ten ,   fatigue / σ ten ,   yield [-]0.800.80-
Ratio τ s ,   fatigue / τ max [-]0.800.80-
Damaging factors δ s , δ n 3.003.00-
Poisson’s ratio [-]0.190.190.135
Interparticle sliding friction μ sl [-]0.45-0.45
Table 5. Average lifetime of UHPC and binder in fatigue simulations with varied damaging factors δ n , δ s .
Table 5. Average lifetime of UHPC and binder in fatigue simulations with varied damaging factors δ n , δ s .
MaterialAverage N
( δ n ,   δ s = 3 )
Average N
( δ n ,   δ s = 6 )
Average N
( δ n ,   δ s = 8 )
Binder156 ± 29208 ± 22>300
UHPC19 ± 656 ± 14161 ± 12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rybczynski, S.; Schaan, G.; Dosta, M.; Ritter, M.; Schmidt-Döhl, F. Discrete Element Modeling and Electron Microscopy Investigation of Fatigue-Induced Microstructural Changes in Ultra-High-Performance Concrete. Materials 2021, 14, 6337. https://doi.org/10.3390/ma14216337

AMA Style

Rybczynski S, Schaan G, Dosta M, Ritter M, Schmidt-Döhl F. Discrete Element Modeling and Electron Microscopy Investigation of Fatigue-Induced Microstructural Changes in Ultra-High-Performance Concrete. Materials. 2021; 14(21):6337. https://doi.org/10.3390/ma14216337

Chicago/Turabian Style

Rybczynski, Sebastian, Gunnar Schaan, Maksym Dosta, Martin Ritter, and Frank Schmidt-Döhl. 2021. "Discrete Element Modeling and Electron Microscopy Investigation of Fatigue-Induced Microstructural Changes in Ultra-High-Performance Concrete" Materials 14, no. 21: 6337. https://doi.org/10.3390/ma14216337

APA Style

Rybczynski, S., Schaan, G., Dosta, M., Ritter, M., & Schmidt-Döhl, F. (2021). Discrete Element Modeling and Electron Microscopy Investigation of Fatigue-Induced Microstructural Changes in Ultra-High-Performance Concrete. Materials, 14(21), 6337. https://doi.org/10.3390/ma14216337

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