Preprint
Article

Numerical Approach to Determine the Resistance of Threaded Anchors in UHPFRC

Altmetrics

Downloads

53

Views

28

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

26 September 2024

Posted:

27 September 2024

You are already at the latest version

Alerts
Abstract
In this study, the structural behavior of a short threaded anchor with a 20 mm diameter and an embedment length of 50 mm (2.5Ø) in UHPFRC is investigated using non-linear Finite Element models. UHPFRC is assumed to exhibit tensile strain-hardening behavior, with tensile strengths of 7 MPa and 11 MPa, respectively. The modelled anchor was subjected to a continuously increasing uniaxial pull-out force. The results indicate that the fracture mechanism of threaded anchors in UHPFRC is primarily characterized by the formation of a tensile membrane within the UHPFRC, which acts as the main resisting element against the pull-out force. Additionally, the influence of UHPFRC's tensile properties on the pull-out behavior and ultimate resistance of the threaded anchors was determined.
Keywords: 
Subject: Engineering  -   Architecture, Building and Construction

1. Introduction

1.1. Industrial Demands and Research Motivation

The use of Ultra-High Performance Fibre Reinforced Cementitious Composites (UHPFRC) in structural engineering, particularly for constructing thin structural elements such of façades, bridges, and highly stressed fixed railway tracks, has revealed a need for anchorage systems with short embedment lengths. Current commercial anchorage systems were originally designed for concrete applications, with a minimum length of 80 mm, which is non-applicable for the constructive design of UHPFRC structural elements with thicknesses between 30 mm and 100 mm.
A recent application in the domain of fixed railway tracks [1] involved replacing damaged concrete monoblocs of 24 cm height by 11.9 cm-thin UHPFRC pedestals after only 40 years of operation. A comparison of project variants showed that, due to its relatively high mechanical strength and durability, UHPFRC can provide significantly more effective engineering solutions, both economically and technically, compared to traditional reinforced concrete. One of the challenges in this project was to define rail fasteners that were compatible with the entire fastening system. These fasteners had to be 20 mm in diameter but required a short embedment length. Additionally, the pull-out resistance of the anchorage in UHPFRC had to be determined. This project highlighted the need to further investigate the pull-out resistance of short anchors to optimize UHPFRC-based fixed railway track systems.
Moreover, the mentioned minimal anchor embedment length of 80 mm in concrete presupposes anchoring within the reinforcement layer near the surface. The pull-out resistance of such anchors in concrete has been well-studied through numerous laboratory tests and investigations, resulting in four main fracture modes [3]: (a) Pull-out of a concrete cone, (b) Fracture of the anchor, (c) Pull-out of the anchor, and (d) Combined fracture of the concrete cone and anchor pull-out.
Since UHPFRC has significantly different properties from concrete, particularly in terms of tensile strength, fatigue resistance, and ductile deformation behavior, the fracture modes shown in Figure 1 for concrete cannot be directly adopted for UHPFRC.

1.2. Research Objective

The main objective of this research is to investigate the structural behavior of UHPFRC under quasi-static, continuously increasing forces applied to short threaded anchors. This will be achieved by using non-linear finite element analysis as the first step in the investigation process.
To the best of the author’s knowledge, no study is currently available on the structural behavior and pull-out resistance of UHPFRC, with tensile properties as shown in Figure 1, for anchors with short embedment lengths subjected to pull-out loading.

2. Materials and Methods

2.1. UHPFRC Material Properties

The Ultra-High Performance Fibre Reinforced Cementitious Composites (UHPFRC) is composed of a cementitious matrix of optimized high packing density, which is strengthened by a high amount of short, slender steel fibers. It has a compressive strength of more than 150 MPa and a tensile strength of more than 12 MPa while developing pronounced tensile strain-hardening and strain-softening behavior [6]. A further relevant property of UHPFRC is its extremely low permeability for water, leading to a highly durable material.
The UHPFRC response in tension (Figure 2) and compression is fundamentally different. Under direct compression, the material remains quasi-elastic almost until reaching the compressive strength. Under direct tension, three domains of material behavior are observed: (1) the elastic domain until the elastic limit stress fUte is reached, (2) beyond, quasi-linear strain hardening occurs when matrix discontinuities form, until the tensile strength fUtu is attained, and (3) gradual softening of a single discrete crack in the weakest cross-section of the specimen, finally leading to material fracture. The fracture zone (3) is a transitional regime of a so-called fictitious crack, developing from a continuum to a real (stress-free) crack.
In addition, the fatigue endurance limit has been determined to be at a maximum fatigue stress level of 70% of the elastic limit stress for UHPFRC subjected in the elastic domain, at about 55 – 65% of the elastic limit stress for UHPFRC subjected in the strain hardening domain and 45% in the strain softening domain.

2.2. Methodology

To determine the failure mechanism of short threaded anchors embedded in UHPFRC and their pull-out resistance, a numerical approach was found to be efficient. The tensile and compressive properties of UHPFRC had already been tested through many uniaxial tensile and bending tests [2], [5]. These tests were used to calibrate the simulated material properties for the numerical model, which was subjected to a continuously increasing uniaxial pull-out force.
The numerical analysis was conducted using the finite element software ATENA[12]. The suitability of this program for the intended purpose was confirmed by the first author’s prior experience in investigating damage of load-bearing structures made of reinforced concrete, including the Gotschna road tunnel and the Adler railway tunnel in Switzerland. Additionally, ATENA was used for designing rehabilitation measures for the fixed railway track at Zurich Airport station using UHPFRC [1,2]. In the damage investigations of the mentioned tunnels, monitoring was conducted, which validated the results of the numerical analysis.
In the frame of this research, the numerical modeling was carried out as follows:
  • Simulation of previously conducted 4-point bending tests, following [4] (pp. 4-5), using UHPFRC properties from test B4 (p. 6). The numerical model was also tested using UHPFRC properties of specimen T1-2 from [3] (p. 5 Table 2, p. 9 ) and specimen TIII 000-1 from [5] (p. 19 Table 3.2, p. 94). The simulated 4-point bending test served as an auxiliary model.
  • Calibration of the finite element mesh using the auxiliary model.
  • Simulation of the numerical model to investigate pull-out resistance.
  • Sensitivity analysis and calibration of the finite element mesh size.
The geometry of the steel anchor was slightly modified for numerical purposes; the flanks of the threaded anchor were modeled with the same surface, thread angle, and pitch as a B500B rebar with a diameter of 20 mm, but without a helix angle.

2. Literature Review

Previous laboratory tests were mainly conducted on fiber-reinforced concretes, high performance fiber-reinforced concretes, or different types of UHPFRC with low steel fiber content, usually 0.1 % – 1.6 %. In addition, epoxy glue and other synthetic materials were used as bonding materials for anchors. In some cases, UHPFRC itself was used as a bonding material for post-installed anchors in concrete [4]. All reported laboratory tests failed by anchor pull-out (Figure 1(c)), fracture of the anchor (Figure 1 (b)), or a combined fracture (Figure 1 (d)), except the particular type of facade anchors discussed below. Thus, the pull-out resistance of UHPFRC could not be directly determined, nor could the structural behavior be investigated based on these laboratory tests.
Satoh et al. [4] conducted pull-out tests on post-installed adhesive anchors in concrete, which were fixed in a borehole of 30 mm diameter using UHPFRC and analyzed the results using a finite element model. The anchors with 25 mm, 50 mm, and 75 mm embedment lengths were tested. The UHPFRC, with a compressive strength of 100 MPa and tensile strength of 10 MPa, had a fiber content of only 1.57 vol. %. The tests revealed that the pull-out of a concrete cone (similar to Figure 1 (a)) occurs in the case of the shortest anchor with an embedment length of 25 mm, and the combined fracture mode was observed in the case of 50 mm and 75 mm embedment lengths (Figure 1 (d)). This research showed that using UHPFRC instead of the grouting mortar increases the bond between anchor and concrete and thus increases the pull-out resistance of anchors fixed in concrete. In addition, the finite element analysis [4] showed that critical cracks appeared in concrete and led to the formation of the fracture cone.
Kalthoff et al. [5] investigated the pull-out behavior of threaded anchors with a relatively short pitch in UHPC (ultra-high performance concrete) with 0.62% to 1.86% fiber content and concrete. The pull-out resistance in the tested UHPC was up to 53% higher than in concrete. Moreover, it was pointed out that using short anchors in concrete is unreliable and thus not recommended. A good bond was achieved just in the case of anchorage in UHPC. Furthermore, the same fracture mechanism for UHPC and concrete was initially assumed. However, no specific investigation of fracture mode was conducted, and thus, the results could not be used for the interpretation of the fracture mechanisms observed in this work.
Fehling and Lorenz [6] investigated the bond behavior of rebars in UHPC. They performed three series of pull-out tests on rebars with diameters ds = 12 mm and 14 mm anchored in UHPC with 0.1% to 1.5% steel fiber contents, in which different parameters were varied:
  • The first test series investigated the influence of UHPC cover thickness on the fracture mode. The pull-out force was applied to the rebars with embedment lengths ranging from 2ds to 12ds anchored in the UHPC specimens with the UHPC cover of ds to 2.5ds. It was found that the increase of UHPC cover leads to the disproportionately lower increase of the average bond strength at maximum force. Furthermore, splitting of UHPC was observed in the case of rebars with an embedment length of 8ds anchored in UHPC of 1% fiber content and UHPC-cover smaller than 2.5 ds.
  • The effect of conventional steel reinforcement bars deployed in the transverse direction was investigated in the second series. The pull-out force was applied to the rebars with an embedment length of 2ds to 5ds anchored in the specimens with a UHPC cover of 1ds to 2.5ds. The results show that transverse reinforcement bars do not increase the pull-out resistance in the tested case.
  • The influence of fiber content on bond strength was investigated in the third series, which was performed on the UHPC specimens of the length of 780 mm, width of 300mm, and height of 360 mm with the anchored rebars with embedment lengths of 3.5ds and 5ds. The steel fiber content was 0.5, 1 and 1.5%. The type of fibers and the heat treatment were the same as in the second test series. The UHPC cone fracture type was observed in the case of UHPC cover larger than 2ds and embedment length of 3.5ds to 5ds. In the case of a very large (not exactly defined) UHPC cover, a combination of UHPC-cone and splitting cracking was observed.
The tests reveal that UHPC with fiber content up to 1.5% behaves similarly to concrete. Thereby, the UHPC cone fracture in analogy to Figure 1a occurs just in the case of the shortest embedment length, which was of sufficient size to allow for fracture cone formation. Moreover, on the photo showing the fracture mode (Figure 7 in [6]), radial cracks are identified, which indicates that the tensile strength of UHPC in the tangential direction around the rebars was reached.
Grzesiak et al. [7] investigated facade anchors fixed in high-performance polypropylene fibre reinforced concrete (HPFRC). The stiff polypropylene fibres with a diameter of 0.70 mm and a length of 30 mm are three times larger than usual steel microfibers (for instance, the steel microfibers in a commercially available UHPFRC-mix are 13 mm long and have a diameter of 0.175 mm). Two types of facade anchors with an embedment length of 25 and 26 mm were used. In both cases, the fracture mode was a pulled-out HPFRC cone, which is very similar to the typical fracture cone in concrete, according to Figure 1a. Thereby, the inclination between the fracture surface and the specimen surface was about 35°. In the case of the adjustable suspended tension anchor, the pull-out resistance of the HPFRC cone increased up to 250% with increasing polypropylene fibre content from 10 to 35 kg/m3. The pull-out tests for the adjustable spacer bolt showed no significant improvement in the pull-out resistance due to increased fiber content.
The research findings indicate that increasing steel microfibre content increases the pull-out resistance. Generally, a better bond can be achieved between a threaded steel anchor and UHPFRC than in concrete. With increasing steel microfibre content, the fracture mechanism changes from a pure fracture cone, according to Figure 1a, to a combined fracture mechanism, according to Figure 1d. Furthermore, in UHPFRC with polypropylene fibers, just as in a steel fiber-reinforced concrete, an increase in the pull-out resistance of the fracture cone is observed. The fibers that traverse the fracture surface are anchored on both faces of the crack, transmitting tensile stress.

3. Numerical Simulation

3.1. Aim

Current non-linear FE programs allow for accurate modeling of the properties of cementitious building materials and structural elements to simulate the structural behavior up to maximum resistance and into the post-peak domain. These FE models are often used for detailed analysis of laboratory test results and investigation of the mechanical behavior of structures, considering specific loading conditions. Such finite element simulations presuppose that the cementitious building material is idealized as a continuum with a predefined constitutive stress-strain law.
In the following, the structural behavior of UHPFRC near a short, threaded anchor subjected to uniaxial quasi-static loading is investigated by means of non-linear finite element analyses. Moreover, it will be demonstrated that numerical simulation is an efficient approach for designing laboratory tests or in-situ monitoring of structures.

3.2. Numerical Model

The presented numerical model corresponds to a common and frequently encountered application scenario for increasing the load-bearing capacity of a reinforced concrete (RC) slab due to the application of a UHPFRC layer and the need to fix a tensile element on the top of a slab. To the existing RC element, a 50 mm thin layer of UHPFRC is added. An alternative application case to the building of the numerical model could also be a model that represents a thin UHPFRC structural part, but, in this case, an additional boundary condition would be needed to obtain compliance of this thin structural part under the pull-out force.
Further, as known from numerous laboratory tests and realized projects, a monolithic bond is obtained between the concrete and the UHPFRC layer. Thus, the numerical model assumes a monolithic bond between both building materials and that no cracking or fracture occurs in the concrete layer below the anchor. Based on this, assigning only UHPFRC properties to the macro-elements representing the UHPFRC part of the composite slab is sufficient, making the numerical analysis easier and more efficient. Figure 3a presents the concept of the numerical model. A rectangular element ① of 0.5mx0.5mx0.5m is cut out from a large composite slab, consisting of a concrete part with a thickness of 0.45 m and a UHPFRC layer of 50 mm, which are monolithically connected. A short, threaded anchor of 20 mm diameter ②, with an embedment length of 50 mm and six force-transmitting flanks, is fixed in the UHPFRC layer during casting. The size of the anchor flank surface corresponds to a common rebar. Figure 3b shows the numerical model using the finite element program ATENA [8]. The anchor is modeled with six 5 mm wide flanks distributed evenly along the embedded part according to Figure 4. Because of the non-orthogonal shape, the tetrahedra mesh with an element size of 3 mm was assigned to the anchor. The simulated slab part is subdivided into three areas. The inner orthogonal macro element around the anchor with a size of 0.1mx0.1mx0.1m (Figure 3b) has the same mesh type and size as the anchor. The macro-element anchor and the macro-element UHPFRC are monolithically connected to simulate the UHPFRC – steel bond behavior. The outer upper area of the UHPFRC element is subdivided into eight orthogonal macro elements, to which the hexahedra mesh type of the element size of 23 mm was assigned. The lower area of the UHPFRC element serves only for self-weight and is statically irrelevant, as only the upper area will be affected, as shown below. This area was subdivided into nine orthogonal macro elements with hexahedra mesh of 100 mm element size for geometrical reasons. All macro elements representing UHPFRC are monolithically connected. The size of finite elements in the numerical model was iteratively determined, calibrated using the auxiliary model according to 2.2 and tested to ensure force transmission from the anchor to the UHPFRC over the surfaces of the flanks.
Note that the size of the representative slab part (Figure 3a) has been chosen so that the horizontal supports ④ are placed outside the area with UHPFRC – anchor interaction. The self-weight of this representative slab part is greater than the maximum pull-out force, such that the reaction of the vertical support ③ is always in compression. In addition, since the focus is to analyze the structural behavior of UHPFRC and how the pull-out force is introduced in UHPFRC structural parts, the numerical simulation of anchor slip in UHPFRC is not performed, since the term slip refers to the limited abrupt displacement of the anchor, during which no force is introduced into the UHPFRC. A sensitivity analysis was performed to justify this assumption. For this, the interface elements between the anchor and the UHPFRC macro-elements were inserted into the original model, and the horizontal stiffness (friction) was varied. This sensitivity analysis showed that no slip was observed even at very low friction in the case of this specific numerical model, and as a result, the compressive force introduced in UHPFRC above the anchor flanks remains the same during each loading step in the numerical models with and without interfaces. The only difference between both models obviously was the computing time, which was longer in the case of the model with interfaces. This sensitivity analysis does, however, not claim that a slip between UHPFRC and the anchor could not occur in real structural elements.
Further, Figure 4 shows the concept of force introduction from the threaded anchor into the UHPFRC bulk material, as considered in the numerical model. The applied pull-out force causes tensile stress in the anchor, which leads to compressive stress acting on the contact surfaces between the anchor flanks and the UHPFRC. The resulting inclined upward compressive stress fields in UHPFRC lead to the formation of a horizontal near surface tensile stress field in UHPFRC, as will be discussed below. The maximum compressive force that can be introduced into the UHPFRC over a flank is limited by the compressive strength of UHPFRC.

3.3. Modeling of UHPFRC

The theoretical background of the numerical modeling of UHPFRC and the programming aspects are described in detail in [13]. The following is the summary of the relevant steps of numerical computation in ATENA. Table 1 summarizes the mechanical properties of two UHPFRC products used in this study. Both commercially available UHPFRC products were already used in previous research for laboratory tests. UHPFRC 7 corresponds to TIII 000-1 according to [5] (p. 94) and UHPFRC 11 corresponds to T1-2 according to [3] (p.9).
The mechanical properties of UHPFRC were modeled using a fracture-plastic material model (CC3DNonLinCementitious2User [13]), in which the constitutive models for tensile and compressive behavior are combined. The fracture model is based on the smeared crack model that presupposes the Rankine failure criterion. A pointwise defined tensile softening law close to a linear function is used (Figure 5). The plasticity model is based on the Menétrey-Willam failure surface. The fracture and plasticity models are combined by the strain decomposition method after De Borst [13,14] and the return mapping algorithm is used to determine the elastic predictor. Due to strain decomposition after De Borst [14] into elastic εije, plastic εijp, and fracture εijf components (εij= εije+ εijp + εijf, I,j=1,2,3,4), the stress equivalence in both models is guaranteed. The new stress state is determined according to the formula σijn= σijn-1+Eijklεkl- Δεklp - Δεklf), where the Δεklp and Δεklf are increments of plastic and fracture strain. As mentioned above, the Rankine criterion Fif= σijt-fti≤0 is considered in the numerical simulation of the softening behavior [13], where σijt is the tensile stress calculated according to the formula: σijt = σijn-1+Eijklεkle calc+ Δεklp calc + Δεklf calc), fti is the tensile strength in the principal stress direction i, and k is the failure surface. If the calculated tensile stress σiit is higher than the tensile strength fti, then, in the next iteration step n, the tensile stress is recalculated such that the final tensile stress satisfies the condition: Fif= σiin-fti =0.
Assuming that the increment of fracture strain Δεklf is normal to the failure surface k, it follows that Δεklf=Δλ·∂Fkf/∂σij= Δλ·δik, where Δλ is the increment of the fracturing multiplier. By the substitution of the tensile stress σiin= σijt- Eijkl·Δεklf and the increment of the fracture strain Δεklf into Fif= σijt-fti =0, the fracturing multiplier is determined as follows: Δλ= σkkt-ftk/Ekkkk. Further, the crack opening wk max is determined by the sum of the total value of the fracture strain εkkf and the current increment of the fracturing strain Δεkkf = Δλ·δkk, multiplied by the characteristic length Lt: wk max =Ltkkf+ Δεkkf). It should be noted that, in the present study, the characteristic length Lt is defined as the projection of the geometric size of the finite element on the axis in the direction of the principal strain and corresponds to the characteristic crack1 bandwidth [11], as shown in Figure 3b.
Further, Figure 5 presents the tensile law of UHPFRC, modeled pointwise in the elastic, strain-hardening, and strain-softening domains. Until the localization of strain (εtloc), the stress-strain curve corresponds to the stress-strain constitutive law. After exceeding the localized strain εtloc, the material law is adjusted to the actual crack band Lt [16] as follows: εfl current= εtloc+( εfl- εtloc)/Lt/Ltch, thereby, Ltch is the measuring base used to determine the strain value.
Note that the mentioned crack band model was initially developed for concrete. Since concrete is a cementitious material with low tensile strength and without strain-hardening domain, tensile strain softening occurs immediately after reaching the tensile strength when the fictitious crack starts to develop. Unlike concrete, UHPFRC shows matrix discontinuities and no real crack until reaching the deformation w (εfl=lf/2) corresponding to half of the maximum fiber length lf. Therefore, the crack band model can be interpreted as a discontinuity band model in the case of UHPFRC.

4. Numerical Simulation of the Quasi-Static Structural Behavior of UHPFRC Substrate

4.1. Modeling of Quasi-Static Loading

The quasi-static pull-out force is modeled as a gradually increasing uniaxial force with 5kN per loading step acting at the top of the steel anchor. To simulate the post-peak domain, the arc-length solution method was used [13]. Note that for studying this model, knowing the exact force introduced into the UHPFRC body is crucial. Therefore, the introduced force in the UHPFRC body is monitored on the level of the top UHPFRC surface, as shown in Figure 7. This modeling results in pull-out force–displacement curves shown in Figure 6. The finite element simulation ran until no more change of the introduced pull-out force was observed in the post-peak domain. The numerical simulation was interrupted manually when at constant pull-out force the displacement at point O (Figure 7) was increasing which was interpreted as pulling out of the threaded steel bar from the UHPFRC substrate.

5. Results

5.1. Pull-Out Force-Displacement Curves

To investigate the effect of increasing the tensile strength of UHPFRC on the force-displacement behavior and the fracture mode, the force–displacement diagram was normalized. Figure 7 presents the normalized pull-out force–displacement curves obtained for UHPFRC11 and UHPFRC7. Thereby, the maximum vertical displacement is assumed to be the displacement at which no more force change was registered; respectively, the threaded anchor will be pulled out from the UHPFRC substrate. It can be seen that the ultimate force Fu of the anchorage was achieved at 16% of the maximum vertical displacement for UHPFRC11 and 30% of the maximum displacement for UHPFRC7.
Analysis of the absolute magnitudes of force and displacement (Figure 6) reveals that the computed ultimate force of the anchorage is similar, i.e. 234 kN in UHPFRC11 and 227 kN in UHPFRC7, reached at the same vertical displacement of 0.7 mm in both cases. However, the post-peak domain is different. In UHPFRC7, the post-peak domain shows relatively little deformation while the force decreases to 70% Fu with increasing vertical deformation thereafter at constant force, indicating pulling-out mechanism. In the post-peak domain of UHPFRC11, the force decreases to 33% Fu, then only the vertical deformation near the anchor increases.
In both cases, no real crack occurs in the tensile membrane during the entire loading process. The strain increase in the membrane leads to a widening of the opening between UHPFRC and the anchor. This means that the ultimate force of the anchorage does not correspond to the ultimate resistance of the UHPFRC membrane. Moreover, the numerically simulated pull-out behavior shows three domains with different slopes, corresponding to three stages of the pull-out process. Therefore, for the detailed analysis of the numerical results, the entire loading process is subdivided into three phases: Phase I: loading up to 30% Fu, Phase II: loading up to 75% Fu, and Phase III: loading up to 100% Fu and post-peak. Using plots of principal strains at each loading phase and in the post-peak domain, the tension and compression fields and, thus, the fracture process are analyzed as follows.

5.2. Pull-Out and Fracture Process

Figure 9a,b show plots of maximum principal strain at the end of phases I, II, and when achieving 100%Fu. These plots are supplemented with red and blue arrows to clarify the development of tension and compression fields. The red arrows highlight the tensile strain propagation in the radial direction, and the blue arrows show the direction of the compression field extension extracted from the plots of minimum principal strain. Positive values of principal strains correspond to tension and negative values to compression.
In phase I, which extends from 0 kN to 30% Fu, the pull-out process starts to propagate in the upper part of the UHPFRC layer around the threaded anchor within an approximately 5 mm thin horizontal circular tension strain field (εrad). Additionally, an inclined tensile strain field (εcone) develops perpendicular to the UHPFRC cone forms.
These strain fields reveal that the pull-out force introduced via the threaded anchor produces initially conical compression () and circular horizontal tension fields (εrad) in the near surface UHPFRC layer. If no fracture of the steel anchor occurs, and compression (σcompession(-ε)≤fUc) and tension (σtensionrad) ≤fUtu) are transmitted in the UHPFRC substrate, the anchorage is functional. Thereby, the UHPFRC next to the threaded anchor enters into the strain-softening domain, already at F=10-15% Fu. UHPFRC in the softening domain in the state at the end of phase I (30%Fu) is shown with the orange color for UHPFRC7 and the purple color for UHPFRC11 according to the scale on the right of Figure 9a,b. Further, the inclined tensile strain field (εcone) extends over the entire UHPFRC region that comprises the threaded anchor and the UHPFRC cone, including the mentioned circular horizontal tensile (εrad) and conical compression () strain fields. In addition, phase I is characterized by the quasi-linear path on the force-displacement diagram indicating elastic behavior.
During phase II, extending up to 75% Fu, the horizontal circular tensile strain field (εrad) increases further and encompasses the UHPFRC cone (lcone) completely, as depicted by the plots of principal strains.
Moreover, with increasing circular tension strain field (εrad , lmembrane), the height of the UHPFRC cone decreases (lcone → 0). And at the end of phase II, the anchorage mechanism consists of alone the circular tension field of thickness (lmembrane) which is equal to the embedment length of the threaded anchor lembedment. From a static viewpoint, the circular tension strain field (εrad) behaves like a tensile membrane. On the force–displacement diagram, phase II is characterized by non-linear force-displacement behavior between 30% Fu and 75% Fu.
Phase III begins with the fully formed tensile membrane of thickness lmembrane ≈1.0- 1.1 lembedment as the main resisting part against the pull-out force. Since the pull-out force acting through the threaded anchor leads to compression at the contact surfaces between the anchor flanks and UHPFRC (Figure 4), the inclined conical compression strain field (, Figure 8a,b) is also observed within the tensile membrane and causes additional increasing tensile strain in the upper third of the tensile membrane.

5.3. Failure Process within the Tensile Membrane

Next, the failure process within the tensile membrane is analyzed. Figure 9 presents schematically the strain (εrad) and stress (ft) distribution within the tensile membrane. Figure 12 and Figure 13 show the state of the strain distribution in the radial direction ( εrad) in the annular circular tensile membrane every 54 kN until Fu is reached. The step of 54kN was chosen iteratively to monitor the development of the tensile membrane. The dashed arrows indicate the radius of the tensile membrane that corresponds to the loading given on the dashed arrows. According to the legend below the diagrams, the color curves present the strain distribution within the tensile membrane in the mentioned loading steps. From these diagrams, the development of the tensile membrane and the strain distribution during the entire loading process can be depicted. The diagram shows that the maximum radius (Rmembrane) is 250 mm, and the outer part of the membrane, corresponding to 60% of the membrane radius, remains in the elastic domain during the entire loading process. UHPFRC next to the anchor enters into the softening domain, still within the loading phase I according to Figure 7. The comparison of Figure 10 and Figure 11 reveals that the tensile membrane in UHPFRC11 is bell-shaped and more pronounced in the inner part than in UHPFRC7, whereas the tensile membrane in UHPFRC7 is less pronounced and shows the shape of an arch.
Further, the diagrams shown in Figure 12 provide the main parameters, i.e., the radius of the tensile membrane Rmembrane and the thickness of the membrane lmembrane during the loading process. During loading phases I and II, the tensile membrane forms, and the rate of increase in the radial direction is up to 5 times higher than in the vertical direction.
The tensile stress remains constant within a limited part of the membrane with a constant diameter where failure occurs. However, more volume is activated within the substrate around the anchor where the tensile membrane develops, and the tensile stress is distributed such that the outer part of the tensile membrane, corresponding to 60% of the radius, remains elastic. Thus, the membrane maintains force balance during the entire loading process.

5.3. Features of the Failure Mode

This study aims principally to investigate the structural behavior of the anchorage of common threaded steel bars in UHPFRC with given tensile properties as subjected to pull-out force using finite element analysis. The investigation of the structural behavior of UHPFRC was possible because the numerical modeling allows to assign sufficiently high tensile strength to the steel to avoid anchor fracture such that the failure mechanism in UHPFRC becomes the determinant failure mechanism of the system.
The numerical simulations show that the introduction of a pull-out force in the UHPFRC body via a short, threaded anchor leads to the formation of a tensile membrane in the near surface UHPFRC with a thickness 1 – 1,1 lembedment. The outer part of this tensile membrane, which corresponds to about 60% of the membrane radius, remains in the elastic domain during the entire pull-out process while the thickness and radius of the tensile membrane increase. Thus, the fracture process zone is localized in the inner part of the annular circular tensile membrane of UHPFRC, corresponding to about 35% of the membrane radius. In the radial direction, the tensile membrane develops up to 5 times faster than the thickness in the vertical direction. At about 85% of the ultimate force Fu, the tensile membrane is completely formed with the thickness corresponding to the anchor embedment length lembedment and radius of about 5 lembedment. The stress transfer in the completely formed tensile membrane occurs further via compressive stresses acting on the contact surface between anchor flanks and UHPFRC, forming compression stress fields and increasing tensile stresses in the upper third of the tensile membrane. Finally, anchorage failure occurs due to the widening of the inner part of the annular circular UHPFRC tensile membrane around the anchor and, as a result, pulling-out the anchor from the UHPFRC substrate. Thereby, no real cracks occur.
In summary, this study has revealed that the pull-out behavior of threaded anchors in a UHPFRC substrate fundamentally differs from the well-known pull-out behavior in concrete, according to Figure 1. For short, threaded anchors in UHPFRC, only two fracture modes are possible: (1) Fu,anchor > Fu,membrane : pulling-out of the threaded anchor from the UHPFRC substrate (fracture mode (1) Figure 13a) , or (2) Fu,anchor < Fu,membrane: fracture of the anchor (fracture mode (2) Figure 13b).
Another feature of the structural behavior of the anchorage system is the plastic deformation in terms of failure indication:
  • In the case of the tensile resistance of the anchor being lower than the pull-out resistance of the UHPFRC substrate (Fu,anchor < Fu,membrane), the plastic deformation of the threaded anchor would occur. This is the relevant case for the design of anchors in structural applications in engineering practice.
  • In the case of the tensile resistance of the anchor being higher than the pull-out resistance of the UHPFRC substrate (Fu,anchor > Fu,membrane), the tensile membrane deforms within the inner 35% of the tensile membrane radius Rmembrane around the anchor. Comparing the vertical deformations in UHPFRC11 and UHPFRC7 at the instant of anchorage failure, the vertical deformation of the annular circular tensile membrane with higher tensile properties was found to be more pronounced and localized around the threaded anchor. However, the vertical deformations remain very small compared to the fracture deformation of the anchor in steel. These vertical deformations are only 2mm for UHPFRC7 and 4.5mm for UHPFRC11. However, this failure mode can only occur if the threaded anchor of diameter 20 mm is made of high-strength steel (ft ≥ 700MPa).
The tensile stress remains constant within a limited part of the membrane with a constant diameter where failure occurs. However, more volume is activated within the substrate around the anchor where the tensile membrane develops, and the tensile stress is distributed such that the outer part of the tensile membrane, corresponding to 60% of the radius, remains elastic. Thus, the membrane maintains force balance during the entire loading process.

6. Conclusions

  • This study shows that FE modeling is an efficient approach to discover possible failure mechanisms before laboratory tests are designed to focus on the research-relevant case. FE modeling allows to vary building material properties such as to avoid failure modes that are not research-relevant, i.e., fracture of the steel anchor.
  • Due to the UHPFRC’s relatively high tensile strength and tensile deformation capacity, significant tensile stress fields form a tensile membrane in the UHPFRC substrate near the anchor when subjected to increasing quasi-static pull-out force.
  • With sufficient substrate size (i.e., the distance from the anchor to the element edge should be longer than five times the embedment length of the anchor), the tensile membrane develops both in horizontal and vertical directions. The tensile membrane’s maximum extension was about five times the membrane thickness in the horizontal radial direction.
  • The tensile membrane’s outer part, corresponding to 60% Rmembrane, remained in the elastic domain during the loading process until pulling-out of the threaded anchor.
  • No fracture of the UHPFRC substrate was observed; the failure of the anchorage occurred due to the widening of the inner part of the UHPFRC annular circular tensile membrane and pulling out of the threaded anchor.

Author Contributions

Conceptualization, A.H. and E.B.; methodology, A.H, E.B.; validation, A.H., E.B; formal analysis, A.H.; investigation, A.H.; writing—original draft preparation, A.H.; writing—review and editing, E.B.; supervision, E.B.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This article is a revised and expanded version of the conference paper titled “Numerical Simulation of the Structural Behavior of Threaded Anchors in UHPFRC”, presented at the 3rd International Interactive Symposium on Ultra-High Performance Concrete, Wilmington, Delaware, USA, June 4-7, 2023 [17], and “Failure Mechanism of Threaded Anchors in UHPFRC as Simulated by Finite Element Analysis”, which will be presented at the 4th International Symposium on Ultra-High Performance Fibre-Reinforced Concrete (UHPFRC), Menton, France, October 21-23, 2024 [18].

Conflicts of Interest

The authors declare no conflicts of interest.
1
The definition of "crack" in the case of UHPFRC means multiple fine cracking of the cementitious matrix of UHPFRC. Tensile stress is transmitted by the steel fibres. In this paper, this type of cracking is called “matrix discontinuity”.

References

  1. Hochuli A, Brühwiler E, Trouillet J.-G., "Improvement of fixed railway tracks using fatigue-resistant cementitious fibre reinforced composite material, Proceedings of The Fifth International Conference on Railway Technology: Research, Development and Maintenance, ", Montpellier, France. [CrossRef]
  2. Hochuli A, Trouillet J.-G., Brühwiler E, Schulz O, Hoffmann M, "Renewal of a fixed railway track using UHPFRC – part 2", Bautechnik, 2022, Volume 99, pp 10-17. [CrossRef]
  3. Shen X, Brühwiler E, “Influence of local fiber distribution on tensile behavior of strain hardening UHPFRC using NDT and DIC”. Cement and Concrete Research, 2020, Volume 132,. [CrossRef]
  4. Shen X, Brühwiler E, Peng W “Biaxial flexural response of Strain-Hardening UHPFRC circular slab elements”. Construction and Building Materials, 2020, Volume 255,. [CrossRef]
  5. Oesterlee C, Structural Response of Reinforced UHPFRC and RC Composite Members, doctoral thesis, Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Thèse N°4848(2010).
  6. Makita T, Brühwiler E, "Damage models for UHPFRC and R-UHPFRC tensile fatigue behaviour, "Engineering Structures, Bd. 90, pp. 61-70, 2015. [CrossRef]
  7. Eligehausen R, A. Cook R, Appl J, "Behavior and Design of Adhesive Bonded Anchors, "ACI Structural Journal, vol. 103, Nr. 6, pp. 882-832, 2006. [CrossRef]
  8. Satoh A, Sakagami Y, Mitarai S, "Finite Element Analysis on Pull-out Behavior of Post-installed Adhesive Anchor Filled with UHPFRC," Proceedings of the AFGC-ACI-fib-RILEM International Conference on Ultra-High Performance Fibre-Reinforced Concrete, Montpellier, France, 2017.
  9. Kalthoff M, Raupach M, "Pull - out behaviour of threaded anchors in fibre reinforced ordinary concrete and UHPC for machine tool constructions," Journal of Building Engineering, Vols. 33, no 2021, no. 2021, 2020. [CrossRef]
  10. Fehling E, Lorenz P, "characterisation of rebars anchorage in UHPC, "in RILEM-fib-AFGC Int. Symposium on Ultra-High Performance Fibre-Reinforced Concrete, Marseille, France, 2013.
  11. Grzesiak S., Pahn M., Schultz-Cornelius M., Susanne Bies N., "Influence of Different Fiber Dosages on the Behaviour of Façade Anchors in High-Performance Concrete," MDPI CivilEng, 2021. [CrossRef]
  12. Červenka Consulting s.r.o., FEM Programm ATENA (Vers. 5.7.0), Praha, 2021.
  13. Červenka V, Jendele L, Červenka J, ATENA Program Documentation Part 1, Theory, Prague, 2021.
  14. De Borst R, "Non-linear analysis of frictional materials, "TU Delft, Delft, 1986, http://resolver.tudelft.nl/uuid:d52ac7f8-d56f-4e0c-ab14-455ed02fcffe.
  15. Cervenka J, Cervenka V, Laserna S, "On crack band model in finite element analysis of concrete fracture in engineering practice, "Engineering Fracture Mechanics, Bd. 197, pp. 27-47, 2018.
  16. Pukl R, Novák D, Sajdlová T, Lehký D, Červenka J, Červenka V, "Probabilistic design of fibre concrete structures," in IOP Conference Series: Materials Science and Engineering, Prague, Czech Republic, 2017.
  17. Hochuli, A., Brühwiler, E., “Numerical Simulation of the Structural Behavior of Threaded Anchors in UHPFRC”, 3th International Interactive Symposium on Ultra-High Performance Concrete 2023, Wilmington, Delaware, USA, June 4-7, 2023.
  18. Hochuli, A., Brühwiler, E., “Failure Mechanism of Threaded Anchors in UHPFRC as Simulated by Finite Element Analysis”, 4th International Symposium on Ultra-High Performance Fibre-Reinforced Concrete (UHPFRC), Menton, France, October 21-23, 2024.
Figure 1. Schematic presentation of the fracture modes of cast-in-place threaded steel anchors in concrete or cementitious material with low tensile strength according to [7].
Figure 1. Schematic presentation of the fracture modes of cast-in-place threaded steel anchors in concrete or cementitious material with low tensile strength according to [7].
Preprints 119436 g001
Figure 2. Tensile response of UHPFRC.
Figure 2. Tensile response of UHPFRC.
Preprints 119436 g002
Figure 3. (a) Schematic illustration of the set-up; (b) numerical model and illustration of the formation of discontinuities and propagation of a fictitious crack in finite element.
Figure 3. (a) Schematic illustration of the set-up; (b) numerical model and illustration of the formation of discontinuities and propagation of a fictitious crack in finite element.
Preprints 119436 g003
Figure 4. Stress transfer from steel anchor to UHPFRC along the contact surface of the anchor rib with the UHPFRC.
Figure 4. Stress transfer from steel anchor to UHPFRC along the contact surface of the anchor rib with the UHPFRC.
Preprints 119436 g004
Figure 5. Assumed tensile stress-strain curves based on the tensile properties of the two UHPFRC products according to Table 1.
Figure 5. Assumed tensile stress-strain curves based on the tensile properties of the two UHPFRC products according to Table 1.
Preprints 119436 g005
Figure 6. Force-Displacement Diagram UHPFRC7 and UHPFRC11.
Figure 6. Force-Displacement Diagram UHPFRC7 and UHPFRC11.
Preprints 119436 g006
Figure 7. Normalized pull-out force–displacement diagrams for pulling the threaded anchor Ø 20 mm lembedment = 50 mm out of two different UHPFRC types.
Figure 7. Normalized pull-out force–displacement diagrams for pulling the threaded anchor Ø 20 mm lembedment = 50 mm out of two different UHPFRC types.
Preprints 119436 g007
Figure 8. (a) Development of the membrane in the UHPFRC7 (Output ATENA principal strain max) (b) Development of the membrane in the UHPFRC11 (Output ATENA principal strain max).
Figure 8. (a) Development of the membrane in the UHPFRC7 (Output ATENA principal strain max) (b) Development of the membrane in the UHPFRC11 (Output ATENA principal strain max).
Preprints 119436 g008
Figure 9. Schematic representation of the stress and strain distribution within the tensile membrane.
Figure 9. Schematic representation of the stress and strain distribution within the tensile membrane.
Preprints 119436 g009
Figure 10. Strain-distribution in the tensile membrane for UHPFRC11: (a) strain distribution 0-44 0/00, (b) strain distribution 0-4 0/00 .
Figure 10. Strain-distribution in the tensile membrane for UHPFRC11: (a) strain distribution 0-44 0/00, (b) strain distribution 0-4 0/00 .
Preprints 119436 g010
Figure 11. Strain-distribution in the tensile membrane for UHPFRC7: (a) strain distribution 0-14 0/00, (b) strain distribution 0-1 0/00 .
Figure 11. Strain-distribution in the tensile membrane for UHPFRC7: (a) strain distribution 0-14 0/00, (b) strain distribution 0-1 0/00 .
Preprints 119436 g011
Figure 12. Radius (Ra) and thickness (Th) of the membrane, with increasing pull-out force.
Figure 12. Radius (Ra) and thickness (Th) of the membrane, with increasing pull-out force.
Preprints 119436 g012
Figure 13. Failure modes of the threaded anchors in UHPFRC.
Figure 13. Failure modes of the threaded anchors in UHPFRC.
Preprints 119436 g013
Table 1. Tensile properties of the two used UHPFRC products.
Table 1. Tensile properties of the two used UHPFRC products.
Tensile properties of UHPFRC
according to [3] (p.5 Table 2, p. 9) and [5] (p. 19 Table 3.2, p.94)
UHPFRC11 UHPFRC7
Tensile strength fUtu [MPa] 11.5 7.1
Elastic limit fUte [MPa] 9.8 6.9
Strain in UHPFRC at the tensile strength εUtu [‰] 2.29 1.5
Strain in UHPFRC at the elastic tensile strength εUte [‰] 0.22 0.15
Fiber length lf [mm] 13 13
Fiber diameter df [mm] 0.175 0.16
fiber content Vol. [%] 3.8 3.0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated