A thermo-mechanical model for SFRC beams or slabs at elevated temperatures

The bearing capacity of steel fibre reinforced concrete (SFRC) at elevated temperatures is the subject of significant ongoing research, as the effect of steel fibres on concrete performance at higher temperatures is poorly understood. On one hand, steel fibres increase the average thermal conductivity of the concrete cross section and lead to greater heating within a concrete structural member during fire exposure, and on the other, fibres reduce crack widths and prevent excessive spalling. The former effect negatively impacts SFRC performance at high temperature, whereas the latter effect protects the inner structure from direct fire exposure. Additionally, the decreasing strength of steel at higher temperatures can result in the sudden failure of fibre or traditionally reinforced concrete. Within this contribution, a coupled thermo-mechanical model is developed in order to investigate the influence of steel fibres on the thermal loading of concrete. The effect of fibres on heat transmission within concrete, the length of time for which concrete can sustain thermal loads, and on the bending stiffness of reinforced concrete beams or slabs is investigated. The heat transfer process is modelled using Fourier’s partial differential equation of transient heat conduction. A modified plastic hinge model and moment–curvature relations are used to describe stress-dependent deformations. Thus, two alternative approaches are used to adequately track the localisation of damage for single cracks and for distributed and multiple cracking. Thermo-mechanical coupling is achieved by means of temperature-dependent stress–strain relations. These are derived for SFRC based on experimental data from the literature. Experiments are performed in which concrete slabs reinforced with variable amounts of fibres and rebar are exposed to combined thermal and mechanical loadings. The results of these experiments are used to validate the proposed model. The measured and predicted results agree well and indicate that steel fibres have a positive effect on the fire resistance of structures, assuming additional rebar is provided to prevent crack localisation. Additionally, it is shown that temperature fields within concrete remain almost unaffected by variations in fibre content.


List of symbols
Thermal conductivity q Density r, r c , r t , r s1 Stress (general, concrete in compression and tension, in longitudinal rebar) r SB Stefan-Boltzmann constant u, u el , u pl Rotation angle (overall, elastic and plastic portions)

Introduction
Steel fibre reinforced concrete (SFRC) is widely used in civil engineering as it considerably enhances the reliability of reinforced concrete structures. Although many different fibre types are available, mechanically anchored steel macro fibres have been proven to be especially suitable for retarding crack growth in concrete after crack formation. The governing parameters that influence the post-cracking tensile strength of SFRC are fibre type, dosage, orientation and the fibre-concrete bond conditions [1,2]. At subcritical fibre contents, SFRC exhibits a softening material response after cracking. Low-fibre-content reinforcement is therefore often used in conjunction with traditional rebar to ensure the structural integrity of reinforced concrete members on a cross sectional level. In contrast, supercritical fibre contents induce material stiffening after concrete cracking. In common ultimate limit state (ULS) design [3], the necessary rebar reinforcement ratios can be significantly reduced if steel fibres are also included as a component of the reinforcement scheme [4,5]. However, in the case of fire resistant design, the postcracking tensile strength of SFRC is not taken into account by any current or pre-normative standard [6]. It is conservatively assumed that only rebar carries tensile loads at elevated temperatures. This leads to difficulties in determining the necessary residual bearing capacity of a reinforced concrete member [7], as the contribution of the fibres to the strength of the remaining cross-section after loading is neglected. Significant strength reserves may therefore be overlooked in the standard design process. The thermal conductivity and the temperature-variant stress-strain response of SFRC are therefore an object of interest in current SFRC research [8][9][10][11][12][13][14][15][16].
The impact of fibres on heat transmission is a controversial subject among experts. Various sources report that the presence of steel fibres accelerates heating rates within concrete [17], whereas other sources claim that fibres retard the heating process [18]. The enhanced conductivity of steel compared to that of concrete may result in varying temperature distributions within the composite material, which may lead to local damage within a given cross section. Conversely, if cracks develop within a concrete member during heating, steel fibres prevent excessive concrete spalling and thus protect the inner structure from direct damage under sustained fire exposure. This protective effect is, however, only present as long as the heated fibres are able to sufficiently transfer tensile forces across a crack. If the fibres themselves are heated to such an extent that their tensile strength becomes negligible, the fibres cannot fulfil their intended purpose.
In order to predict SFRC behaviour under thermal loading, it is essential to properly analyse the interaction between the thermal and mechanical behaviour of the material. Within this contribution, a coupled thermo-mechanical model for SFRC is presented. The model is used to investigate both the heat propagation as a function of the fibres' volume contents and the overall load-bearing capacity of SFRC members subjected to bending while exposed to fire. Effects of local concrete spalling or time-dependent behaviour, such as creep or viscosity changes, are not considered explicitly. Large-scale fire tests conducted on concrete slabs reinforced with both steel fibres and traditional reinforcement are used to verify the results. The tested slabs are heated according to a uniform temperature time curve according to Eurocode 1-1-2 [19].

Modelling the thermo-mechanical behaviour
The weakly coupled model is derived in three steps. At first, the equation of heat transfer in solid bodies is presented. Second, the bending behaviour of SFRC elements with or without rebar under mechanical loadings is introduced and third, the two strain fields are coupled by overlapping portions and introducing temperature-dependent stress-strain laws. Since stress-states are assumed not to influence heat transmission, the heat flux is derived a priori. Throughout, damage induced material loss due to spalling is neglected so that the initial geometry is preserved.

Heat transmission
The model is based upon the principles of continuous composite body heat transmission. For the sake of simplicity, fire exposure or down-cooling to ambient conditions are assumed to occur at the outer boundaries from convective and radiative portions only while inner heat sources are generally neglected.
The temperature field #(x,t) in the solid depends on the spatial location x and time t, which changes progressively due to temperature exposure and the corresponding temperature gradients entailing a heat flow in the gradient's direction. Considering heat balance in an incremental volume element Fourier's partial differential equation of transient heat transmission [20,21] yields: The solid is assumed to be homogeneous and isotropic, which also holds true for any arbitrary, infinitesimal volume therein. Consequently, the thermal diffusivity a simplifies to a scalar, denoted by a ¼ k=ðq c p Þ, the ratio of the thermal conductivity k to the product of material density q and specific heat capacity c p .
The heat flow volume Q [J] of a solid is computed from its mass m, its specific heat capacity and temperature T [K] relative to absolute zero at 0 K = -273.15°C.
Two volume fractions with a contact area A exchange heat flow volume increments dQ per time increment dt, which is described by the heat flux density _ q.
Since heat flux and the field of temperature gradients are proportionally coupled, Cartesian components can be derived in the three spatial directions x, y and z.
Separation into x, y and z components benefits the numerical solution of heat transfer by discretising space and time domains [22,23] using consistent rectangular volume elements and time steps Dt, respectively. Figure 1 illustrates such a spatial discretisation of a reinforced beam with boundary contact elements simulating ambient or fire conditions. The arrangement is chosen to effectively derive transient temperature propagation in spreadsheet environments [24][25][26][27]. Material elements are idealised by specific cells in the spreadsheet surrounded by others to model boundary conditions of convection and radiation.
Starting with an initial temperature field #(x,t) at a time t, the current distribution at the end of a time step t ? Dt is obtained by explicit time integration. Temperatures are assumed to remain stepwise constant in Dt while the heat flow changes. For convenience derivations are restricted to two-dimensional cases in the remainder.
Cross-sections are discretised by rectangular elements. Each element i is characterised by an associated area A i = x i z i , its average temperature # i at a time t and its thermal diffusivity a i . Internal elements are surrounded by g = 4 solid elements. By contrast, boundary elements exhibit one or two free surfaces (r) that reduce the number of the inner surfaces to g = 4r. An alteration of heat flow volume (DQ i / Dt) results from the sum of up to four densities _ q i along the inner or boundary contact lengths l i,j and l i,r .
Thereby, the directed heat flow density between two differently heated internal elements i and k depends on the element lengths s i , s k , the thermal conductivities k i , k k and the temperature difference between # i and # k .
Boundary conditions, i.e. heat sources or cooling at the surfaces, are treated similarly. For this purpose the corresponding component is replaced by a fictitious transfer quantity from or to the outside of the body.
Convective and radiative portions are modelled by an equivalent temperature # ID depending on a fictitious heat-transfer coefficient a K [28] and the surface temperature # 0 [22,29].
In case of fire exposure, the radiative heat flow _ q rad becomes a function of the grey body emission coefficients of the solid e m and the surrounding fumes e f [21]. r SB denotes the Boltzmann constant.
The radiation temperature # A and the surface temperature # 0 are both related to absolute zero. For the sake of simplicity, # A can be set to the discharge temperature of the fumes, # A = h [28] and # 0 to the corresponding element's one # 0 = # i . Hence, the temperature field at the end of a time step reads: To ensure numerical stability the time step size Dt [22,23] has to be limited. For two-dimensional cases it depends on thermal parameters (c p , q, k) and element sizes (x i , z i ): To reduce the computational demand the check might be restricted to the smallest element of a discretisation.

SFRC beams or slabs under mechanical loading
Elevated temperatures considerably affect the mechanical material response of plain concrete (PC), SFRC and steel rebar. Stresses are a function of strains and both depend on temperature. Introducing a separation of strains into thermal (superscript th) and purely stress-dependent (superscript r) portions according to [7] gives: Thus, mechanical and thermal analyses can be conducted independently.
A softening material effect due to elevated temperatures can be considered in mechanical analyses only. In this case e th purely arises from temperature differences D# and the coefficients of thermal expansion a T (#).
Mechanical modelling considers bending moments M and axial forces N only, shear deformations are neglected. On cross-sectional level strains are distributed linearly over the cross-section's height h. Moreover, the bond between concrete and rebar is assumed perfect and slip free. Uniaxial stress-strain laws describe the mechanical response of PC, SFRC and rebar which depend on the current temperature of the material and thus change with time. Applying equilibrium conditions for a two dimensional case yields: introducing stresses r c of PC or SFRC, the crosssectional area A, n layers of reinforcement A s with stresses r s and the vertical coordinate z. Figure 2 illustrates the situation for a standard beam with one reinforcement layer at the bottom. A nonlinear stressstrain law for SFRC is applied considering the beneficial effect of fibres in a smeared way by additional tensile bearing capacities.
For an iterative numerical solution of Eqs. (13,14), the cross-sectional height is partitioned into equally spaced layers, corresponding to the discretisation of the thermal analysis to allow direct superposition of thermal and mechanical portions according to Eq. (12). The appropriate strain plane is derived using optimisation methods [24,30] that robustly and efficiently solve the nonlinear equilibrium problem by minimisation of residuals between outer forces N and M and corresponding stress integrals [24].
To facilitate a transfer of the proposed calculation method to the longitudinal direction of beams or slabs it has to be distinguished between single (localised) and multiple cracking.

Multiple cracking
Bernoulli's theorem of thin and transversely rigid beams enables the derivation of temperature-variant curvatures j(#) which depend on the strains in equilibrium at the outermost fibres of a cross-section. These are generally located at the bottom e c,bot and at the top e c,top (Fig. 2) Since both strains are induced by stresses, j reflects the ratio of bending moment (M) to bending stiffness (E c I) in non-cracked as well as in cracked states. Thus, an individual moment-curvature relation has to be derived for each cross-section being considered. The principle of virtual work and Simpson's rule are employed to numerically integrate the curvature j(x) and the virtual bending moment MðxÞ along x resulting in the vertical deflections d(x). Of course, a suited number of integration points along the length l has to be provided to adequately cover stiffness reductions due to cracking. Furthermore, tension stiffening might be approximated by modified stress-strain relations of rebar [5].

Localisation due to single cracking
The softening post-cracking response of SFRC beams containing subcritical amounts of steel fibres but no additional rebar causes a localisation of deflections at single cracks [2,3,32]. Within the fracture process zone a discontinuity, often denoted a plastic hinge, forms that governs the beam's deflections by its rotation. Rotation depends on the strain distribution in the region that develops nonlinearly over the width of the hinge. To transfer localisation to a quasi-continuous approach, Strack [33] developed a plastic hinge model. It is based on average strains of the tensile and compressive fibres in the discontinuity region, while equilibrium is checked at the cracked cross-section only, Fig. 3. In this regard, uniaxial material functions by means of stress-crack-width relations r(w) according to [32,33] are transferred to stress-strain ones r(e) and vice versa introducing a characteristic length l ch according to the notion of a crack band [34]. Assuming linear elastic stress-strain behaviour in the noncracked state, the post cracking strains consist of elastic (e el ) and plastic (e pl ) portions: In case of common macro fibres the characteristic length l ch corresponds to about twice the fibre's length and thus to an almost constant value. Subsequently, l ch is set equal to 140 mm based on [2]. This value approximately refers to fibre lengths of 50-80 mm and-at least for thin structures in bending-resembles l ch & h, as advised in [3].
Strack's plastic hinge model was originally proposed for SFRC subjected to short term loadings at normal temperatures [33], but allows for generalisations in case of fatigue [35] or elevated temperatures. Elevated temperature requires customisation of the material behaviour only. Moreover, the calculation procedure is more reliable when load instead of deformation controlled. While the geometry of the beam, its loading pattern and temperature-dependent material behaviour are prescribed for each discrete time step, the strains at the border of the hinge are iteratively computed. Subsequently, the basic model is , left: structural system as well as geometry and deflection characteristics; right: strain and stress distributions briefly summarised in its key aspects respecting temperature dependencies and axial forces for the sake of simplicity. Presuming a stress distribution at 45°from the compressed to the tensile edge, the longitudinal extent of the plastic hinge equals about twice the crosssection's height, s = 2 h [33]. In non-cracked regions (subscript ncr), strains e ncr are obtained based on linear elastic theory.
Here, W ncr denotes the sectional bending resistance and M the bending moment, which can be induced by three (e = 0) or four-point (e [ 0) bending from the force F. Average tensile ( e c; bot ) and compressive ( e c; top ) strains at the edges describe the hinge's overall deformation. They depend on the deflections in equilibrium and result from the tensile and the compressive zone of the cracked cross-section, respectively. Within the tensile zone an artificial strain e c;bot ¼ e el þ e w is introduced by summing up average elastic strains ( e el ) and the maximum crack width w smeared over the hinge's width [33]. Thereby, a linear crack opening over the cross-section's height is assumed.
Cracking occurs if the induced strains exceed e t,cr , which equals the ratio of concrete's tensile strength f t to its Young's modulus E c . The average compressive strains follow iteratively considering the height of the compressive zone x c = n h and exhibit a distinct nonlinear course as indicated by e c,top [33] in Fig. 3 The same characteristic applies to the corresponding deflections d that are superimposed analogously. with and Therein, I denotes the geometrical moment of inertia of the non-cracked cross-section and l the total length of an assumed beam section.

Coupling
Thermo-mechanical coupling is achieved via temperature-dependent stress-strain relations. Starting from initial conditions with a constant temperature state in the solid (e.g. # 0 = 20°C at t = t 0 ) the temperature distribution at the end of a time step t ? Dt follows according to Sect. 2.1. It affects the stress-strain relations in each lamella of the discretised body and thus governs sectional equilibrium (Sect. 2.2). Since material's stiffness and strength decrease progressively with increasing temperature, strains have to increase to ensure equilibrium of forces. If equilibrium cannot be achieved anymore, the structure fails. Figure 4 recapitulates this basic calculation scheme.
To obtain total deflections d, the corresponding stress (d r ) and temperature (d th ) induced portions are additively superimposed. Since distributed cracking is expected for SFRC structures reinforced with additional rebar and Eq. (16) serves to derive d r , it must be replaced by Eq. (24) in case of single localised cracking.
d th depends on the thermal strains e th , which usually turn out highly nonlinear distributed over the height with maxima on the heated face. Since e th equals to the product of temperature difference and thermal expansion coefficient-Eq. (12)-the temperature gradient can be calculated instead to obtain d th . For convenience, the temperature field over the cross section's height is partitioned into constant, linear and residual nonlinear parts [36]. Only the linear part-the temperature gradient D#-causes vertical deflections. The others induce axial deformations or constraints. The gradient D# reads or when discretised with i = 1 to m lamellae A i over the cross-section's height Employing the principle of virtual work, thermal deflections are obtained integrating the product of thermal curvatures j th = (a T D#)/h and a virtual bending moment MðxÞ along the effectively heated length l eff,h of a specimen [31].

Thermal parameters of SFRC
The conductivity of SFRC k SFRC is proposed to be computed from a weighted average of volumes of concrete and steel. Thus, a homogenous composite is assumed which holds true, if the representative length of specimens and discretisation considerably exceeds the internal lengths of fibres and aggregates. Introducing the fibre volume content V f and the conductivities k c and k s of concrete and steel, k SFRC reads: While literature on the rule of mixture is generally ample (e.g. [8,38,39]), an increased fibre volume content V f in Eq. (30) yields enhanced conductivities and vice versa. By contrast, the fibre's impact on concrete's density, specific heat capacity and heattransfer as well as thermal expansion coefficients, is considered subordinate [40] and thus neglected. Pure concrete characteristics are assumed.

Stress-strain relations of SFRC
Temperature-dependent stress-strain relations of concrete in compression [7] are adopted for SFRC, too, since they do not remarkably differ up to the softening branch [10] as it is with normal temperatures [1][2][3]. By contrast, a more specific behaviour is proposed in tension. The discrete crack-bridging ability of fibres in fire is captured in a smeared way introducing a multilinear stress-strain relation r t (e t ,#). Based on the material response at a reference temperature of # = 20°C, softening with increasing temperatures is modelled by a global, multiplicative softening parameter k t f (#). It represents the ratio of concrete's degraded to initial strength and is defined in the range of 0-1. To assess the temperature-dependent stiffness a temperature-variant Young's modulus of concrete E c (#) is adopted [7]. It represents an artificial numerical value to also integrate the effects of high-temperature creeping [36].
Linear-elastic behaviour is assumed up to the tensile strength f t . A bilinear function is used to model the fibre's post-cracking bearing capacity employing two residual strength values f 1 f and f 2 f at fixed strains e 1 and e 2 , respectively. The concept of two strain values with associated strengths is chosen to comply with the established procedure of material classifications via standard bending tests in [1][2][3]. The strains e 1 = 3.5% and e 2 = 25% correspond to crack widths of w 1 & 0.45 mm and w 2 & 3.15 mm employing l ch-= 140 mm according to [2]. While the strength at w 1 is associated with the serviceability limit state (SLS), the strength at w 2 corresponds to ULS conditions [1,4].
The definition of temperature-invariant, constant strain values e 1 and e 2 in Eq. (31) leads to a conservative underestimation since [7,37] report on enlarged strain boundaries for PC and rebar under sustained fire exposure. Figure 5 illustrates the derived stress-strain relations in tension and compression. It should be noted that discontinuity in the plotted temperaturedependent stress-strain curves is virtually due to intended cuts in the strain scales highlighted by dashdotted lines. Fictive heat-transfer coefficient [28] Non temperature-exposed concrete surface 4 To consider radiative portions at the non temperature-exposed concrete surface 9 Temperature-exposed concrete surface 25 a T (#) Thermal expansion coefficient a T = e th /D# Concrete Acc. to EC 2-1-2, 3.3.1 (1), [7] Steel Acc. to EC 2-1-2, 3.4 (1), [7] Experimental data from literature [9,[11][12][13][14][15][16] as well as closed formulas already suggested for PC [7] and SFRC [40] give advice to capture the concrete's temperature-dependent strength independent of a material's composition, the crack width and the mechanical loading factor (SLS or ULS). Then, the relative softening parameter k t f (#) relies on temperatures only.
It is proposed as a conservative, trilinear function according to Eq. (32) evaluating literature data, shown in Fig. 6 (dashed line). Softening starts, when # approximately exceeds 150°C and continues up to # = 700°C. Beyond, no bearing capacity is assumed. As intended, Eq. (32) accurately reflects test results by trend but still conservatively [40].
Since test data reported in [41][42][43] confirms that bond between rebar and concrete is not significantly influenced by steel fibres, neither at normal nor at elevated temperatures, no additional bond conditions are specified for the stress analyses.

Verification to large-scale fire tests
For verification purposes large-scale fire tests are carried out in cooperation of the Technical University of Kaiserslautern, the Ruhr-University of Bochum (both Germany) and NV Bekaert SA (Belgium). All specimens exhibit the same geometry and concrete strength class C35/45, but are reinforced with different ratios of steel fibres to rebar. The general idea is to vary the ratio in such a way, that reinforcement type changes but bending capacities at ambient temperatures remain almost constant.

Experiments
In total five slabs are exposed to combined thermomechanical actions. Table 2  Eq. (32) Fig. 6 Experimental data indicates the course of the relative softening parameter k t f regarding applied temperatures calculated bending capacity for ambient conditions at midspan that nearly coincides for P1, P2 and P5. To achieve maximum homogeneously dispersed fibres in the concrete, the mix contains a great amount of cement paste and fines. Thereby, the particle-size distribution of the quartzitic aggregates is continuous with a maximum grain diameter of 16 mm. Two types of Dramix steel fibres with one or two-times hooked ends and aspect ratios of 65 are used in dosages of V f = 0.5 and 1.0 vol%, respectively. The mean postcracking tensile strengths f 1 f and f 2 f of the fibres have been determined before fire testing on small scale reference beams according to [2]. Despite a strong inherent scatter of the post-cracking strength [35], the fibre contents at hand are always subcritical: f . However, the distinct calculative bearing contributions of fibres regarding M t0 emphasise fibres' governing impact even when combined with considerable amounts of conventional rebar. For experiments, rectangular mats of cold-formed rebar, steel class N [7], are used.
To ensure moisture contents of practical relevance, i.e. of about 4 vol% at testing, the specimens are conditioned for at least three months after concreting. For the experiments, an enclosed fire room is created exposing the single-span slabs on the exterior walls of the furnace and providing insulation elsewhere. Two half shells serve as line supports. Thus, the effective length of the specimens follows to l eff ¼ 400 À 16 ¼ 384 cm, whereas only the inner part of 300 cm is directly exposed to fire, (Fig. 7). Hence, the effectively heated length l eff,h lies between l eff and 300 cm. Motivated by thermal images (Fig. 7, right) it is assumed to the balanced average of l eff,h = (384 ? 300)/2 = 342 cm.
Mechanical loadings are applied by means of steel plates on top of the slabs. They induce a bending moment of about 17 kNm at mid-span. Heating is applied to the bottom-side according to a standardised temperature time curve [19] for 120 min at maximum with a fume temperature of As failure criterions, a maximum deflection d max or rather a maximum deflection rate _ d max apply [44] that both depend on the static height d of the slab (d & 15 cm) and its length.
Vertical displacements are continuously recorded at the quarter points of the specimens while temperature data is collected at 25 discrete measurement points uniformly distributed over the slab's volume. Panel thermometers monitor heating of the furnace by gas powered burners.

Fire resistances and failures
The purely fibre reinforced samples collapsed early after 7 (P4) and 41 (P5) minutes of fire exposure. A single and widely opened crack characterised failure and was accompanied by intensive spalling at the bottom face. This can be traced back to two effects: first, a scattering post-cracking tensile strength of SFRC and second, the statically determined single span system that does not allow stresses to redistribute and induces local cracks at the weakest section. By contrast, P1, P2 and P3

Temperature fields
For numerical analysis of heat transmission the slabs are discretised in the x-z-plane while the transverse direction is neglected. Thus, a constant temperature over the width is supposed as documented by another thermographic image in Fig. 7, right, too. In spreadsheet analysis [23][24][25][26], the specimen and its boundaries are idealised by a set of perpendicular cells, similar to Fig. 1. Each cell is attributed either to an element made from concrete, rebar or to environmental conditions, i.e. exposed to fire or to ambient air. Individual element sizes depend on the material and the position within the solid. Here, the slabs are subdivided into 40 elements lengthwise and 25 elements over the height, adaptively reducing their sizes to consider rebar and concrete cover. Reinforcing mats are represented by two layers of elements superimposed upon each other. The lower layer models bars in longitudinal direction, the upper one those across. Up to 120 min at maximum, the transient temperature fields are continuously and explicitly updated. Element sizes and assigned thermal parameters limit the size of a time step to Dt = 1 s regarding Eq. (11).
Subsequently, the numerical results-particularly with regard to the smeared parameter of thermal conductivity k SFRC -are validated by experimental data. On the left, Fig. 8 displays general temperature distributions after 15, 60 and 120 min of heating. For the same points in time, calculated temperature distributions over the section's height (coloured dashed lines) are contrasted to experimental ones on the right. By trend all dots fall on a straight line (perfect match) whereas different fibre volume contents of V f = 0, 0.5 and 1% can be distinguished. In detail, temperatures are assessed at the top and bottom side as well as at five additional points across the thickness. As the inner sensors could not be fixed directly to the formwork, they are more sensitive to slipping during concreting. Deviations from the sensors' reference positions might have caused some scatter in the results. As expected, numerical and experimental results document clearly nonlinear temperature gradients over the height. Generally, the temperature increase lasts on the longer the fire burns, whereas the trend mitigates with less distance to the fire-exposed side. While maximum temperatures occur at the fire facing side, minimum values are found on the opposite face.
Theoretically, increasing the fibre volume contents goes along with enhanced body temperatures due to accelerated heat flux. The dashed lines in Fig. 8, right reveal this aspect by certain offsets that rise with time. However and despite the scatter inherent to measured temperatures, no significant influence of the fibre amount (V f B 1.0 vol%) on the concrete's heating rate becomes obvious during the experiments. Moreover, recalculations employing unmodified thermal parameters of concrete with respect to [7] (k SFRC = k c according to [40]) match to the test data best. Assuming V f in correspondence to its real physical values even slightly reduces the accuracy of recalculations.

Transient bearing capacities
Based on the transient temperature fields, temperatureand stress-dependent deflections are obtained. Since   Table 2 or supplemented according to standard code provisions of [45].

Cross-sectional level
To capture the total impact of heating on the bearing behaviour, the cross-sectional response is analysed first. For sections with rebar (P1 to P3) exposed to fire for 120 min and a maximum surface temperature of 1000°C moment-curvature relations are derived. In Fig. 9 they are contrasted to the ones obtained initially for a uniform temperature of h = 20°C. As expected, maximum bending capacities decrease for all three fibre contents (P1, P2, P3) with increased temperatures. However, this is still true after heating. Similar to additional bearing capacities provided by rebar, fibres contribute to the overall resistance even when exposed to fire. Although strength and stiffness are on a significant lower level then. Most surprising fibres remain beneficial even if-from the model's point of view-almost all performance within the concrete cover at temperatures higher than 700°C is lost.

Structural level
All five slab tests are finally recalculated and compared to the experimental results. To achieve transient evaluations, the mid-span deflections d are contrasted over a continuously increasing fume temperature h in Fig. 10  While temperature effects dominate at that point, stress induced portions become important.

Conclusions
Within this contribution, a thermo-mechanical model is proposed to investigate transient bearing capacities of SFRC slabs with and without additional rebar exposed to elevated temperatures. The temperaturedependent tensile strength of the SFRC is derived based on test results from the literature. The accuracy of the coupled model is validated using results obtained from large-scale fire tests performed on reinforced concrete slabs. Both tests and numerical data are, on average, in good agreement. Based on the presented research, the following conclusions can be drawn: • Concrete's heating rate is not significantly enhanced by steel fibres in dosages up to about 1.0% in volume. The thermal parameters describing heat capacity, density and conductivity as well as the thermal expansion coefficient of PC can be adopted for SFRC as a first order estimate. • The temperature-dependent but strain-independent softening of fibre's post-cracking tensile strength can be reliably predicted by means of a bilinear function, as originally proposed in [40]. This function assumes undisturbed conditions until the solid temperature reaches 150°C and assumes a full loss of bearing capacity at 700°C. • SFRC without rebar might extensively spall when subjected to fire. Brittle failure was observed in the single-span systems very early after fire exposure, i.e. after 7 and 41 min, due to single cracking. In both cases, the radiation temperature of the fumes did not exceed approximately 600°C. • SFRC containing additional rebar exhibits multiple cracking, pronounced deflections and reduced spalling. The ratio of specimen's length to its maximum vertical deflection reaches about d = l/ 30. Thus, a ductile material behaviour prevails up to about 1000°C of outer exposure. • The combination of mechanically anchored steel fibres and rebar results in enhanced fire resistance durations and higher bending stiffness compared to reinforced concrete. However, the necessary residual strength values corresponding to fire resistance durations of 90 min (R90) could be met in both cases. [°C] [