Incremental-iterative model for time-variant analysis of SFRC subjected to flexural fatigue

Fatigue behaviour of plain (PC) and steel fibre reinforced concrete (SFRC) is of growing interest in concrete engineering since structural reliability often depends on the concrete’s damage state. The current paper deals with its investigation. Starting from a new and universal SN-approach for PC and SFRC based on stress- and material-dependent ductility and a cycle-dependent strain evolution under centric pulsating loads an isotropic and time-dependent material damage parameter is derived. In the framework of the elasto-plastic damage theory, wherein PC and SFRC are idealized homogeneously on macroscopic level, this damage parameter in conjunction with the established envelope concept enables to compute time-dependent stiffness and strength degradations as well as increasing plastic strains. Additionally, an assessment of the total number of cycles to failure of specimens subjected to multi-staged cyclic loading is permitted by a specific damage accumulation procedure. Decisive parameters like fibre type, length, orientation, dosages and bond on the axial static and cyclic material response are covered. To compute numbers of cycles to failure, deflections and stress redistributions of members exposed to flexural fatigue, the findings are integrated into a plastic hinge model. This accounts for localization of damage after cracking in a discontinuity region typical for both, PC and SFRC. Cycle- and crack-width-dependent deflections are obtained for valid states of equilibriums between external loads and internal stress resultants. To assess macroscopic cracking a strain criterion serves. In the numerical simulation, time-increments are performed according to Lemaitre’s jump-in-cycles procedure. For verification test results from literature are recalculated. Theoretical and experimental data are in good accordance on average proving a comparable fatigue behaviour of PC and SFRC in principle. However, fibre’s impact on fatigue life and deformation capacities of concrete is ambivalent, depending on specific material and load characteristics.

Abstract Fatigue behaviour of plain (PC) and steel fibre reinforced concrete (SFRC) is of growing interest in concrete engineering since structural reliability often depends on the concrete's damage state. The current paper deals with its investigation. Starting from a new and universal SN-approach for PC and SFRC based on stress-and material-dependent ductility and a cycle-dependent strain evolution under centric pulsating loads an isotropic and time-dependent material damage parameter is derived. In the framework of the elasto-plastic damage theory, wherein PC and SFRC are idealized homogeneously on macroscopic level, this damage parameter in conjunction with the established envelope concept enables to compute time-dependent stiffness and strength degradations as well as increasing plastic strains. Additionally, an assessment of the total number of cycles to failure of specimens subjected to multi-staged cyclic loading is permitted by a specific damage accumulation procedure. Decisive parameters like fibre type, length, orientation, dosages and bond on the axial static and cyclic material response are covered. To compute numbers of cycles to failure, deflections and stress redistributions of members exposed to flexural fatigue, the findings are integrated into a plastic hinge model. This accounts for localization of damage after cracking in a discontinuity region typical for both, PC and SFRC. Cycle-and crack-width-dependent deflections are obtained for valid states of equilibriums between external loads and internal stress resultants. To assess macroscopic cracking a strain criterion serves. In the numerical simulation, time-increments are performed according to Lemaitre's jump-in-cycles procedure. For verification test results from literature are recalculated. Theoretical and experimental data are in good accordance on average proving a comparable fatigue behaviour of PC and SFRC in principle. However, fibre's impact on fatigue life and deformation capacities of concrete is ambivalent, depending on specific material and load characteristics. Today, steel fibre reinforced concrete (SFRC) is widely used in structural engineering since mechanically anchored fibres are suited to reduce crackwidths and spacings and enhance structural integrity significantly compared to plain concrete (PC). The governing but also interacting parameters on macroscopic level are fibre type, bond, dosage and orientation as well as concrete's composition [1][2][3]. Usually, SFRC exhibits a strong ductile but softening axial material behaviour after cracking due to progressive fibre pullout in static loading. For specific applications strain hardening behaviour can also be ensured by employing high dosages of well anchored fibres. Those fibres of elaborated shape ensure strong bond behaviour with concrete's matrix and often fail with fibre tears. Nevertheless, Young's modulus and concrete strength are not significantly enhanced compared to PC, neither when subjected to compression nor tension.
In case of cyclic loading, internal degradation mechanisms of PC and SFRC are similar and can be compared to the ones in monotonic loading [4][5][6]. Most results available in literature agree that static load deformation curves form an envelope of deformations at failure observed in constant amplitude tests [7][8][9][10][11]. However, this observation seems to hold true in case of centric rather than eccentric loading only, since e.g. the authors in [12] attribute measured deviations from the envelope to a size-dependent fracture process zone in case of flexural fatigue mainly. Actually, repeated eccentric loading enables a redistribution of stresses on cross-sectional level that leads to more reversals than observed in axial tests [13]. Concurrently, the stress redistribution sets higher demands on numerical simulation of concrete's material behaviour being met e.g. by plastic hinge models.
To assess the number of cycles to failure (N f ) in constant amplitude, low and high cyclic loading, SN or so called Wöhler-curves are well established (e.g. [14,15]). For these, the stress level S stat , the stress ratio R according to Eq. (1) as well as the frequency F and the waveform of the load history are considered to be the governing parameters of fatigue life.
Consensus is a greater number of cycles endurable decreasing S stat or increasing R or F, respectively (e.g. [16,17]). In comparison, alternating stresses, also called stress reversals (R \ 0), and rectangular waveforms cause failure after lower numbers of cycles than pulsating stresses (R [ 0) or sinusoidal and triangular waveforms [18]. However, an influence of the loading type (uni-or multiaxial, centric or eccentric) on the fatigue life obtained from SN-curves can nearly be eliminated relating the maximum fatigue strength r f max to the corresponding static one f c/t . Accordingly, in case of flexural fatigue, the stress-level S stat is calculated from the ratio of maximum flexural stress and the flexural tensile strength of concrete. Likewise centrically induced compressive (subscript c) or tensile (subscript t) stresses are related to the material's axial strength and thus treated conceptually analogue [19]. Doing so, effects of concrete's composition by means of grading curve, pore volume, water-cement-ratio and moisture conditions are included sufficiently, too [20].
In contrast to these generally accepted governing factors, an impact of fibres on the fatigue behaviour of concrete is discussed controversially among experts. It is reported to depend at least on fibre type and loading conditions, e.g. [5,[21][22][23]. Anyway, to employ the established SN-curves for PC containing variable amounts of steel fibres as well, modifications become necessary.
In Lee and Barr [24] an overview of recent achievements in research on fatigue behaviour of PC and SFRC in the framework of Brite-EuRam [25] is provided. From this it is concluded that the threephased process of progressive internal structural changes (interpreted as damage) until fatigue failure is similar for PC and SFRC. The increase of strains with cycles runs generally affine to well-known s-shaped creep curves. But, fibres might accelerate or decelerate degradation in certain circumstances which is traced back to a dual impact they might have. On the one hand fibres increase the pore and initial micro-crack densities, on the other hand they are able to bridge micro-cracks and retard their growth [26,27]. To accurately capture this dual impact of fibres on the fatigue life of concrete, a performance based parameter shall be identified and considered additional to the loading conditions as usual.
In the remainder the fatigue behaviour of PC and SFRC subjected to centric and eccentric loads is illuminated. For this purpose, equations for the axial and time-invariant material response of concrete are presented. Subsequently, a cycle dependent material damage parameter is derived and verified using experimental data from the literature. The new approach is enhanced, transferred to members subjected to bending and finally verified using flexural fatigue tests from the literature.

Uniaxial fatigue of PC and SFRC
Starting from a macroscopic description of static and uniaxial material response in tension according to Strack [28] a new SN-a approach is proposed. To provide clearly arranged equations, the derivations are limited to tension subsequently. However, the procedure can easily be extended to compressive loading conditions using the formulae presented in [29] for instance. The introduced new ductility index a in conjunction with stress level S stat and ratio R as well as frequency F and waveform (a F-t ) are used to derive consistent Wöhlercurves valid for PC and SFRC. Then, s-shaped cyclic creep curves in accordance with the envelope concept are presented taking account for dual fibre effect. In the framework of elasto-plastic damage theory, in which PC and SFRC are assumed to be homogenous materials on macroscopic level, time-dependent degradations of stiffness and strength can be assessed. Additionally, an approach for damage accumulation in case of multistaged centric cyclic loading is proposed.

Uniaxial material response in tension
The phenomenological modelling with respect to [28] superimposes effects of fibre geometry, dosage and bond behaviour with effects related to PC matrix and fracture energy. Up to the tensile strength f t linearelastic material behaviour according to Eq. (2) is assumed.
Material softening after cracking is captured by a mechanically based stress crack-mouth-opening relation r t (w). It additively superimposes two stress contingents. The first one is of exponential type and covers the PC matrix contribution with respect to the fracture energy G f .
A multi-factored approach models the second one (cf. Eq. (4)) employing two efficiency factors j I and j II assigned to the steel fibres effect that considers crack-width-dependent pullout and bond conditions. j I gathers the fibres' stress-state with respect to length l f , diameter d f , dosage V f and orientation g u as well as the mean bond stress s b . Since it is dimensionless (stress is related to the concrete's tensile strength) j I lies between the limits of zero and one, while zero marks PC behaviour. j II covers shape-bond and degradation effects with increasing crack-width by an associated fibre-(l h ) and pullout-length (l peak ). Thereby, a pull-out length of l peak = 0.3 -0.7 is typical for mechanically anchored steel fibres, while j II = 1.0 indicates exclusive frictional but no shape resistance, as with straight fibres [28,29].
with: Figure 1 recalls the complex load transfer mechanisms of discrete fibres bridging crack edges. Postcracking tensile forces can be covered as long as the crack-width is smaller than the total bond length l b , which is equal to half of the fibre's length l f .
Hence, an integrative stress approach is favored here which smears stresses over a crack's length. In Fig. 1 the development of g bond with increasing crack width w according to Eq. (4) is displayed too. It reaches its maximum at 1.0 if the crack width equals l peak and decreases to j II at a crack width of l h , which itself depends on the fibre-type.

Modified SN-a approach considering ductility of PC and SFRC
A stress-and material-dependent ductility index a independent from size effects and the loading rate requires pooling of the static uniaxial stress-strain and stress crack-mouth-opening relations given in Eqs. (2) and (3). Coupled, both form a continuous material model for SFRC covering all intermediate material states up to failure. The conceptual idea behind is to transfer the well-known envelope concept of PC for application with SFRC which says: On an equivalent stress level, the strains at failure e fmax in a constant amplitude fatigue test and the ones induced by static loading according to the stress-strain curve's decreasing branch coincide (c.f. Fig. 2). Therefore, the stress crack-mouth-opening relation (Eq. 3) based on the ''Fictitious Crack Model'' by Hillerborg [30] is formally converted into a stressstrain one obeying the ''Crack Band Model'' by Bazant and Oh [31]. However, and in contrast to [31], the reference length is not associated with the finite η bond (w) 1 size of the fracture process zone as a multiple of the maximum grain diameter but generally set to w cr,E = 1.0 mm instead. Here, the choice of a constant reference length serves to form a unique envelope of stress-strain and stress crack-mouth-opening relation to assess material's ductility. The outcome is not intended to be used for structural analysis directly, as proposed with the original one according to the ''Crack Band Model''. Figure 2 presents the holistic stress-strain relation derived. The dimensionless index a(r f max ) is then defined according to Eq. (5) as the ratio of elastic to total deformation energy integrated from zero to e fmax . It implies that (steel fibre reinforced) concrete's deformation energy per unit area at cyclic loading correlates to the integral of the stress-strain relation up to failure strain at monotonic loading. Thus, materials enhanced capacity to absorb energy due to the strengthening mechanism of fibres results in enhanced values a. Considering balance of energy, damping effects within hysteretic loops can be neglected since a does not quantify correlations of dissipated energy and fatigue lifetime directly. In fact, it rather serves to classify different concretes based on their deformation energy capacities, interpreted as ductility.
In the framework of the envelope concept, e fmax directly follows from r f max and thus, the influence of fibre type, bond, dosage and orientation on the static and cyclic behaviour is treated analogously. Per definition the ductility index ranges from 0 B a B 1.0, while large values of a mark enhanced ductility features under static and cyclic loading and vice versa.
A sophisticated derivation and detailed verification is contained in Heek and Mark [32]. The ductility index is used to derive SN-a curves semi-empirically: Regression coefficients b 0 ¼ ffiffiffiffiffiffiffi ffi 103 p and b 1 ¼ b 2 ¼ 0:5 Á p therein have been determined via response surface methodology (RSM), evaluating a database of more than 560 fatigue test results from the literature. As suggested in [33] the effect of load frequency and history of the waveform on the logarithmic number of cycles to failure is generally assigned to the ratedependent concrete strength, which leads to a modified effective stress ratio S dyn to replace S stat .
In tension and compression, the rate-dependent strength of PC and SFRC (f c/t,dyn ) in Eq. (7) can be derived using the mean stress rate of cyclic loading.
with a F-t = 1.0 in case of rectangular waveform, a F-t = 2.0 in case of sinusoidal and triangular waveform.
To exclude exponential strength gains at high stress rates, Eq. (8) is limited to the quasi-static domain in dependence on Model Code [1]. Nevertheless, such a restriction is not a limit for most structures of practical relevance. In case of higher stress rates, beyond the scope of this paper, detrimental effects due to appreciable temperature gradients would have to be considered, too [34].
Exemplarily, Fig. 3 presents a subset of experimental fatigue data from the 560 tests already used to fit Eq. (6). Since tensile tests are generally not available yet, the focus is set on compressive and flexural tests of PC and SFRC specimens. In case of flexural fatigue data, inverse analysis using the plastic hinge model presented in Sect. 3 Obviously, the experiments are focused on higher load levels of practical relevance associated with 0.65 \ S dyn \ 0.95. By reason of enhanced numbers of data available, fitting and prediction are statistically more appropriate there as indicated by confidence bounds of reduced extend. Due to the inherent scatter in all fatigue data, serial mean values per load level are calculated and indicated in red. They prove a good accordance with the proposed mean trend functions. For comparison, two representative parameter sets are evaluated employing Eq. (6) that differ in a only. A high ductility index is representative for SFRC and a low one for PC by trend. Thereby, the approach yields enhanced numbers of cycles to failure with increasing capacities of material's deformation energy, expressed by a. This goes along with further test results in [37], which proof a strong correlation between fatigue life and dissipated energy during constant amplitude centric loading of PC and SFRC. Correspondingly, the authors in [26] also observed that the amount of cumulative energy during cyclic bending tests is noticeably related to the maximum number of cycles.

.1 has been
Obviously, the new SN-a approach is able to consider dual fibre effect. It distinguishes between fibres being helpful, since they bridge cracks, or not, since they initiate cracking by means of the stress level S dyn . While at high stress levels additional fibres usually do not influence the number of cycles to failure significantly, their crack bridging ability leads to higher lifetimes at lower stress levels. The approach has originally been derived for the range of R [ 0 but can be extended to stress reversals (R \ 0) as well. Up to now, it has been verified for normal strength (fibre reinforced) concrete with softening material response after cracking in the range of -0.3 B R B 0.5, 1 B F B 20 Hz and 0.55 B S dyn B 0.95 [32].

Evolutions of strains and damage
Progressive damage accumulation in load cycles can be captured by a concrete damage plasticity model. To account for the inelastic material behaviour of concrete, isotropic damaged elasticity and isotropic plasticity concepts are combined herein. Fundamental for the elasto-plastic damage model is an additive split of strains in elastic and plastic portions (subscripts el and pl) according to e.g. [38] and Eq. (9), which leads to the stress-strain relation wherein d t denotes a scalar damage parameter in its common limits of 0 and 1. Employing d t , damaging effects of concrete cracking are modelled by a reduced   Fig. 3 Verification of SN-a approach illustrating the influence of a on N f Young's modulus E c fat . Assuming a constant ratio b t of plastic to inelastic strains the damage parameter has been found according to Eq. (11), time-invariant at first.
Agreeing to results of (compliance) tests reported in the literature, e.g. [12] a constant and time-invariant value of b t = 0.50 is convenient in numerical simulations for both materials, PC and SFRC. Physically interpreted, half of the inelastic strains remain plastic. By contrast, Mark [39] recommends b c = 0.70 for PC in compression.
To elaborate the damage parameter to a cycle dependency, the evolution of strains in hysteretic loops is analyzed next. As reported in [18] and [24] for instance, constant amplitude loading causes a threestaged evolution of strains in case of PC and SFRC similar to creep curves. The first part of retarded crack growth involves the weak regions in the concrete matrix, e.g. due to initial shrinkage cracks, and is termed flaw initiation. The second stage is characterized by inherent flaws progressively growing up to a critical size and exhibits accelerated crack growth. For this stage, several authors like in [7,8] provide evidence for a strong correlation between the rate of crack-growth and the number of cycles to failure. Consequently, the same parameters have to be decisive to describe the evolution of strains that influence fatigue life. Finally, the third stage ends with rupture of the specimen due to accumulations and conjunctions of micro-to macro-cracks and a formation of fracture surfaces. To describe the three-staged cycledependent evolution of strains for both, PC and SFRC, a composite creep function is derived subsequently, incorporating S dyn , R and a [cf. Eq. (6)].
The creep function e c/t (n) in Eq. (12) is formulated free of dimensions in the timescale of related numbers of cycles to failure 0 B n = N/N f B 1. Thus, an impact of specific material characteristics (e.g. strain at failure e fmax ) on the shape of cyclic creep function (especially at phase transitions) could be disregarded with intent to grant a wide field of application. Three individual summands J i (n) are superimposed to trace the strains with time. Analogue, the first phase of retarded crack-growth is described with a rational function followed by a linear one, in line with constant crack-growth, and a final power function indicating accelerated crack-growth up to failure. At the intermediate but load-type dependent points of transition, indicated with n 1 and n 2 , the strain increments De 1 and De 2 with respect to the initial strains e f0 induced by the maximum fatigue stress r f max (cf. Fig. 2) fix the strain evolution curve. The type of the strain increment's function results from regression analysis of experimental data from the literature (cf. Fig. 4). Thereby, a cubic relation of De i and S dyn is postulated to ensure an s-shaped creep function in either case. The best fit coefficients (b, c, k, f) are contained in Table 1.
In Fig. 4 calculated and documented strains from uniaxial fatigue tests performed with PC and SFRC in compression and tension are compared supporting the accuracy of the developed approach. Obviously, all stages of de-and acceleration of strains are reflected accurately.
The cycle-dependent material damage parameter d t (n) can now be expressed analogue to Eq. (11) with respect to time-variant strains e c/t (n) by substituting r t by r f max . Of course, physically, d t (n) must be of monotonically increasing type.

Damage accumulation in variable amplitude loading
To consider multi-staged loading with variable stress amplitudes, continuity of the damage parameter d t (n) is postulated similar to the recommendation in [43] and the experimental verification in [44]. Figure 5 exemplarily illustrates the damage progress related to two successive constant amplitude stress levels S i , while S 1 [ S 2 . Separately analysed, S 1 would cause a damage level at failure d t,1 (S 1 ) according to Eq. (13) which is lower than d t,2 (S 2 ). Here, in a two staged procedure, unloading is assumed after N 1 load cycles (n 1 (S 1 )) while the damage level remains fixed d t,1 (n 1 ) = d t,2 (n 1 ) = d t,1/2 as postulated. According to its original progress curve d t,2 (n, S 2 ) the damage level d t,2 (n 1 ) theoretically corresponds to n 2 (S 2 ). Since N 1 / N f 2 = n 2 (S 2 ) the total number of cycles to failure N f 2 with respect to Eq. (6) has to be modified (N f ;mod 2 ) to ensure N 1 / N f ;mod 2 = n 2 (S 2 ). Thus, in case of unloading in a two-staged variable amplitude scenario, greater numbers of cycles to failure are obtained compared to the original number of cycles to failure which matches physical expectation. Moreover, subsequent monotonic strength tests of run-outs often exhibit higher residual strength f fat t than expected from plain monotonic testing which is usually attributed to consolidation of mortar [7,18]. In the continuous damage accumulation approach proposed in Fig. 5 the small range of d t,2 (n) \ d t,1 (n) must be excluded. Therein, no interactions in case of multi-staged loadings are postulated, which leads to unmodified numbers of cycles to failure N f 1 and N f 2 , respectively. Otherwise, increased numbers of cycle to failure N f 1 (S 1 ) would result, if load cycles of lower amplitudes S 2 were applied previously to S 1 .
In Fig. 6 predictions of total numbers of cycles to failure in a two staged loading procedure are compared to experimental data [45] and a best-fit curve presented therein. On the abscissa, numbers of cycles applied on the first stress level are presented, while the ordinate illustrates residual numbers of cycles to failure (N 2,res ) on the second one. Good accordance is shown except the first third of the unloading branch. However, the well-established Palmgren and Miner rule [46,47] Fig. 5 Continuity of damage during multi-staged axial cyclic loading way unable to distinguish between loading and unloading as well as sequence effects at all.

Incremental-iterative model for time-variant structural analysis
To analyse (fibre reinforced) concrete structures subjected to flexural fatigue, the axial and timedependent material model is implemented into a plastic hinge model originally proposed in [28]. It is briefly summarized in its key aspects first and then modified to account for cyclic loading.

Plastic hinge model for monotonic static loading
Based on parametric non-linear finite element analysis, Strack [28] developed a plastic hinge model taking account for distinct localization of damage after cracking, typical for PC and SFRC with softening material response [1][2][3]. Concrete cracking causes a discontinuity region of limited extent (plastic hinge) characterized by a nonlinear strain distribution over the height (h). The hinge's rotation and the neutral axis are obtained with respect to averaged deflections of the tensile and compressive fibres within the region, respectively (cf. Fig. 7). Equilibrium is checked on cross-sectional level employing lamellae techniques for integration of axial material response. It delivers the crack-width dependent load bearing capacity of the whole structure with respect to the hinge's rotation angle u. Asides the plastic hinge, sections remain plane, e.g. a linear strain distribution over the crosssections height is maintained according to Bernoulli's theory. Thus, in transition the strains e ucr are directly obtained from linear elasticity theory employing bending (M) and resistance (W ucr ) moments according to Eq. (14), in which the bending moment can be induced by three-point (e = 0) or four-point bending (e [ 0).
In contrast to this, the averaged tensile strains e t in Eq. (15) are obtained at the plastic hinge by the sum of  Fig. 7 Plastic hinge model according to [28] elastic strains and maximum crack width w smeared over the hinge's width s = 2 h according to Fig. (7) and [28].
Thereby, a linear crack opening over the crosssection's height is assumed. Cracking occurs, if the induced strains exceed the cracking strain e t,cr . The mean strain e el results from integrating the edge tensile strains e t,u along the plastic hinge. Accordingly, the nonlinear compressive strains e c,o summarize to e c . It mainly depends on a crack-width dependent contraction of the compressive zone n(w) and iteratively follows from Eq. (16).
with: e c;max ¼ e ucr j jþ e c j jÀ e ucr j j arctanð1=gÞÁ g 3 þg ð Þ À g 2 and Thus, edge-rotations of the plastic hinge define the rotation angle h which itself consists of elastic and plastic portions (subscripts el and pl) according to Eq. (17).
with: h el ¼ 4 Á e el j j Likewise, the vertical deflections d of the beam result from superposition of elastic and plastic parts in case of equilibrium of outer loads and inner stress resultants. For statically determinate bending members of length L the crack-width dependent plastic deflections can be read from Eq. (18).
All computations are carried out in a deformation controlled way. Beam's geometry, load pattern and material behaviour by means of the stress-strainrelations as well as one concrete edge strain are provided. Its counterpart on the other side is iteratively computed employing structural equilibrium conditions. From a permissible cross-sectional strain distribution load-size and the beam's deflection are obtained.
Step-wise repetition with continuously increased strain finally delivers the entire loaddeflection curve.

Modified plastic hinge model for fatigue analysis
Subsequently, the presented approach is enhanced to simulate cyclic behaviour. For convenience, the belonging algorithm for incremental-iterative fatigue analysis has been implemented in Matlab Ò . In general, the upper equations remain unchanged. However by contrast, computations are now performed load-and time-controlled. Thus, geometry, load-pattern and load-size-the latter one uniquely defined by means of its cyclic parameters S dyn (S stat , F, a F-t ) \ 1.0 and R [ 0-are initially preset, while mean edge strains at top and bottom of the plastic hinge are computed iteratively. They replace the two unknowns load-size and edge strain of the original approach. Furthermore, it is assumed that damage localization in cyclic loading occurs similar to the monotonic case, as experimental results gained by acoustic emission analysis advice [48].
Corresponding to the number of load-cycles applied, the (degraded) material behaviour is accounted for lamella-wise. To evaluate material's degradation, the load history is incremented using a jump-in-cycles procedure with respect to [49]. At first, in all increments the lamella stresses are assumed equivalent to the ones obtained from the preceding step in equilibrium. These are always related to the degraded strength according to the envelopeconcept to obtain the current stress-level S dyn . Ductility index a and current stress ratio R follow analogue, while the time-invariant stress amplitude Dr = r f max À r f min causes variable stress ratios R over time, c.f. [50]. Employing Eq. (6) gives theoretical numbers of cycles to failure per lamella that are related to the total number of load cycles being applied. The dimensionless number of cycles to failure n along with the strain state and the degraded stress-strain-relation of the preceding time-step consistently enable to update the damage parameter d t (n) according to Eq. (13). Thereby, strains at failure e fmax correspond to the initial value of r f max at the outermost fibre of the hinge (cf. Fig. 2). A progressively increasing damage parameter in case of stress redistributions is granted by considering the degraded concrete strength to compute S dyn as outlined above. Both, damage parameter and stress-strain relation are updated lamella-wise. That way, the total stress and strain states of the plastic hinge are iterated meanwhile equilibrium is preserved. If equilibrium cannot be reached anymore, the belonging time-step is associated with fatigue failure. It should be noted that it is utmost important to choose sufficiently small time-increments, since all are sought to initially have constant stresses here. Thus, PFANNER [43] suggests not to keep constant increments at all but to adapt their sizes according to the expected threephase deformation course. Correspondingly, the timeincrements in phases one (flaw initiation) and three (accelerated crack growth) should be relatively small. Figure 8 simplifies the total approach for fatigue analysis by means of a flowchart. Figure 9 shows the calculated evolution of stresses, damage and ductility index for PC and SFRC over the normalized cross-section's height (0 B z/h B 1.0) in case of constant amplitude flexural fatigue with time (N/N f ), respectively. While the first row contains information on all intermediate steps but separately for PC and SFRC, the second row contrasts both materials but restricted to initial and final states only. Obviously, cyclic loading causes significant stress redistributions with time resulting in progressive unloading of the outermost lamellae of the plastic hinge. This comes along with a contraction of the compressive zone, as experimentally verified in e.g. [23], and a steadily growing crack process zone (subpictures 9a vs. 9b). The evolution of the crack process zone with time differs for PC and SFRC so that the ductility enhancing impact of fibres due to crack bridging leads to larger numbers of cycles to failure N f SFRC [ N f PC here. Compared at fracture, the total damage state d t in case of SFRC is greater than the one for PC (subpicture 9c).
Thereby, only the initial stress-strain relation in tension is assumed to be affected by fibers here. This approach is not restricting when investigating the flexural fatigue response of PC and SFRC, since the degradation process in the compressive zone is generally subordinate to the one in tension due to lower stress levels. It particularly becomes obvious by considering the time dependent evolution of a in compression (subpicture 9d). There, ductility indices are almost identical for concrete with and without steel fibres and approximately constant over time. Small differences are exclusively attributed to the enhanced stress redistribution capacity of SFRC. By contrast, cycle dependent unloading of the outermost lamellae in tension cause progressively increasing ductility indices that significantly differ with fibre's effect.

Verification
To demonstrate practical relevance and to verify the predicted influence of steel fibres on structural behaviour by means of cyclic crack growth and deflections, experimental data from the literature are recalculated subsequently. The four point standard bending tests in [21] enable to analyze the correlation between deflections at mid-span and fibre dosage. Thereby, four different fibre amounts were investigated in total. As an input for the fatigue model, the uniaxial static response has been derived via inverse analysis of the documented static bending tests using the original plastic hinge model as summarized in Sect. 3.1. Main results are recapitulated in Fig. 10.
Both, fatigue model and test data exhibit threestaged deformation characteristics, in which an increase in fibre content, here primarily expressed by j I , causes an increase of ultimate deflections and initial input geometrie, loading conditions, initial material response output data stress and strains, theoretical N f per lamella acc. to eq.   Test result: logN f = 4.7 / 5.2 / 5.9 / 6.1 Fig. 10 Recalculation of test-data taken from [21] specimens comprise beams of different geometry in order to investigate size effects on flexural fatigue performance of PC and SFRC. During testing, a constant deformation rate was imposed until the load dropped to approximately 95 % of the peak load. Then constant amplitude cyclic loading started. Such a specific load history cannot be adapted by the proposed computational model directly, in which not preloaded specimens can be analysed exclusively. Thus, experimental and theoretical values of N f cannot be expected comparable. But actually, fatigue dependent crack growth is the main issue here, so the theoretical numbers of cycles to failure are assumed to be equal to the experimental ones (cf. Fig. 11). Then, experimental and theoretical initiation and growth of fatigue induced cracking are in good accordance on average. Actually, maximum crack width at testing are larger throughout, due to a progressive and irregular reduction of maximum fatigue stresses until the minimum fatigue stress is reached. By contrast, the theoretical model considers just constant amplitude loading for simplification here, which has been the case over a wide range of fatigue life during testing as well. In contrast to experimental data, size effects do not emerge in the recalculation employing the proposed model.

Conclusions
Experimental results published in the literature on the static and cyclic behaviour of PC and SFRC prove a general comparability on macroscale. Subjected to cyclic loads, three-phased deformation and damage evolution characteristics are typically observed. A significant fibre contribution occurs after cracking of concrete and increases both, post-cracking tensile strength as well as ductility, especially when subjected to tension or bending. With respect to the number of cycles to failure, results at constant amplitude testing highlight a potential contrary impact of fibres. Dependent on stress level, load type and fibre characteristics, increased or decreased numbers of cycles to failure are stated. By means of a newly proposed stress and material-dependent ductility index a-that relies on the envelope concept-the dual impact of fibres can sufficiently be covered and integrated into an SN-a approach valid for (fibre reinforced) normal strength concrete with softening material behaviour after cracking. An experimentally proven strong correlation between strain evolution, especially regarding the strain rate in the second phase of fatigue life, and log N f advises to characterize both with the same set of parameters. A consistent integration of the derived cyclic creep curves into an elasto-plastic damage model permits to assess stiffness and strength degradations as well as to monitor the evolution of plastic strains over the entire fatigue lifetime.
A close to reality damage accumulation covering multi-staged load sequences gets possible, claiming continuity of the damage progress. It yields higher numbers of cycles to failure in case of unloading and vice versa.
Inclusion of the fatigue damage model in a modified plastic hinge approach enables to deal with flexural fatigue. It has been proven numerically that -small -medium -large Fig. 11 Recalculation of test-data published in [12] additional steel fibres help to improve to redistribute stresses on cross-sectional level contrasted to PC. Cycle-dependent it comes along with a contradiction of the compressive zone and an increased deformation capacity. Since in PC brittle cracking characterizes macroscopic fatigue failure steel fibres are able to bridge cracks and transfer stresses even after cracking and thus lead to higher numbers of cycles to failure. Obviously, the gains depend on fibre type and stresslevel. In general, cracking starts earlier in case of SFRC than with PC, while at failure, higher damage indices are computed. Thus, the findings confirm the experimentally observed duality of fibers' impact on concrete, which is initiation and bridging of cracks. Which of both dominates in a specific practical case can only be judged considering the interactive contributions of fibre type, orientation and dosage as well as concrete strength and fatigue load by means of S dyn (S stat , F, a F-t ) carefully. The newly proposed fatigue model contributes in this framework.