Acessibilidade / Reportar erro

Fire behavior of shallow prestressed hollow core slabs from computational modeling

Abstract

Prestressed hollow core slabs are one of the structural systems whose use has increased the most in recent years in Brazil due to its efficiency and versatility. They can be used in many types of structural systems, such as masonry, precast concrete, cast-in-place concrete and steel structures. However, there are few analytical models to evaluate the fire behavior of hollow core slabs. In a simplified way, the fire resistance is evaluated indirectly through the minimum distance of the surface in contact with fire to the reinforcement axis. In this paper, some numerical models in finite element software were developed to analyze the variation of temperature with fire exposure time of shallow hollow core slabs, focusing on the presence of voids in the transversal section of the slab. The 500 °C isotherm method applied to 20 cm high slabs confirmed the Standard Fire Resistance obtained from the tabular method. However, when applied to shallow prestressed hollow core slabs that are 16 cm high, the 500 °C isotherm method indicated that the Standard Fire Resistance of these slabs is lower than values obtained from tabular methods.

Keywords:
shallow hollow-core slabs; fire resistance; precast concrete; computational modeling

Resumo

As lajes alveolares protendidas são um dos sistemas estruturais que mais ganhou força nos últimos anos no Brasil devido a sua eficiência e versatilidade. Podem ser empregadas em vários tipos de sistemas construtivos como: alvenaria estrutural, concreto pré-moldado, concreto moldado in loco e estrutura metálica. Contudo, ainda há poucos modelos de cálculo racionais para avaliar o comportamento ao fogo das lajes alveolares, sendo que, de forma simplificada, sua resistência ao fogo é avaliada de forma indireta pela garantia de uma distância mínima do eixo da armadura até a face do elemento estrutural submetida ao fogo. Neste artigo são apresentados alguns modelos numéricos desenvolvidos em um programa de elementos finitos para analisar a distribuição da temperatura em lajes alveolares de altura reduzida quando em situação de incêndio, com foco na influência da presença dos alvéolos. O método das isotermas de 500 °C aplicado a lajes com 20 cm de altura confirmaram alguns dos Tempos Requeridos de Resistência ao Fogo (TRRF) indicados nos métodos tabulares. Contudo, quando aplicados a lajes alveolares com altura de 16 cm, o método das isotermas de 500 °C indicou que o TRRF dessas lajes é inferior ao indicado nos métodos tabulares.

Palavras-chave:
laje alveolar; incêndio; pré-moldado; modelagem computacional

1. Introduction

In precast construction, the hollow core slab is a structural element whose use has been rapidly increased due to its versatility and agility observed in the execution and assembly of structures. This slab can be used in many types of structural systems: in situ concrete, precast concrete, masonry and even in steel structures. Although the use of hollow core slabs in Brazil has increased significantly, in Europe it is produced in a much higher scale due to the industrialization process and project optimization developed in countries such as Germany, France and England.

The hollow core slab is a widespread system and has been deeply investigated in the last years to determine its behavior in fire conditions. There are many experimental and numerical studies which analyze various situations of hollow core slabs regarding the effects of loading and support of the slabs on the rest of the structure.

In Europe and United States, there are standards and handbooks for the design of structures under fire conditions, such as Eurocode 1 [1[1] EUROPEAN COMMITEE FOR STANDARDIZATION. Eurocode 1 - Actions on structures - Part 1-2: General actions - Actions on structures exposed to fire. 2002.] and ASTM-E119 [2[2] ASTM International. ASTM E119-15 - Standard test methods for fire tests of building construction and materials. U.S.A, 2015.], which deal with the main definitions of fire and show standard fire curves. The Eurocode 2 [3[3] EUROPEAN COMMITTEE FOR STANDARDZATION. Eurocode 2 - Design of concrete structures - Part 1-2: General rules - Structural fire design, 2004.] and the ACI-216.1-M-07 [4[4] AMERICAN CONCRETE INSTITUTE. ACI 216.1M-07 - Code requirements for determining fire resistance of concrete masonry construction assemblies (metric). U.S.A, 2007.] standards present design methods and reduction coefficients for material strength with the elevation of the temperature. The PCI [5[5] PRECAST/PRESTRESSED CONCRETE INSTITUTE. PCI Manual for design of hollowcore slabs. 2. ed. Chicago, 1998.] manual and the BS EN 1168 [6[6] BRITISH STANDARDS. BS EN 1168:2005. Precast concrete products - Hollow core slabs. 2005.] can be used to design prestressed hollow cores slabs in fire conditions.

In Brazil, the fire design of concrete structures is covered in the NBR 15200:2012 [7[7] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).] and NBR 14432:2001 [8[8] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 14432: Fire-resistance requirements for building construction elements - Procedure. Rio de Janeiro, 2001 (in Portuguese).] standards. The NBR 15200:2012 [7[7] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).] standard presents only a tabular method for fire design of concrete structures and does not deal in its tables with the design of hollow core slabs. The recent revision of the Brazilian standard for the design of precast concrete structures, NBR 9062:2017 [9[9] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 9062: Design and execution of precast concrete structures. Rio de Janeiro, 2017 (in Portuguese).], presents some criteria for the design of hollow core slabs in fire conditions, but this criterion is based only on the minimum values for the height of the hollow core slab and for the distance from the strand reinforcement to the face of the structural element subjected to fire.

This article aims to determine the temperature profile in shallow prestressed hollow core slabs under fire conditions by means of computational modeling in finite element software. The presence of the voids in the cross-section was also considered. The developed thermal modeling strategies are validated by comparing the results of a thermo-mechanical numerical model with test results available in the literature of hollow core slabs under fire conditions and subjected to bending. The 500 °C isotherm simplified design method is used to estimate the reduction of bending moment with time of fire exposure of a shallow hollow core slab.

It is important to deepen the knowledge about the behavior of hollow core slabs exposed to fire. There are several studies in the literature about reinforced concrete beams and columns exposed to fire [10-13], as well as steel-concrete composite structures [14-17]. There are also several studies about the behavior of hollow core slabs exposed to fire, either by tests or computational modeling [18-27]. However, most of these studies consider hollow core slabs with high height. The lowest height observed in these studies was 20 cm, and mostly with circular voids. Thus, there is still a gap of knowledge about the hollow core slabs with height smaller than 20 cm with non-circular voids, which are employed in buildings with small spans.

2. Fire design methods of reinforced concrete structures

According to Costa [28[28] COSTA, C.N. Dimensionamento de elementos de concreto armado em situação de incêndio. Tese (Doutorado em Engenharia de Estruturas), Escola Politécnica, Universidade de São Paulo, São Paulo, 2008. doi:10.11606/T.3.2008.tde-04092008-155911.
https://doi.org/10.11606/T.3.2008.tde-04...
], the design of a structure under fire conditions is analyzed by the thermal and mechanical responses of the structure under a transient thermal action. When a structure is analyzed at room temperature, that is, considering only the static service loads, the procedure is simpler due to the well-known structural behavior of the elements (beams, slabs or columns) and its failure mechanisms.

Already in fire conditions, because it is an exceptional action, as well as earthquakes or explosions that introduce dynamic actions on structures, the behavior of structural elements is a hardly known. There is a great influence of the thermal properties of the materials in the way that heat is transferred to the structure. There are also the questions about how the temperature of the structure increases when it is exposed to fire and if the fire is contained in a compartment or not.

According to Costa [28[28] COSTA, C.N. Dimensionamento de elementos de concreto armado em situação de incêndio. Tese (Doutorado em Engenharia de Estruturas), Escola Politécnica, Universidade de São Paulo, São Paulo, 2008. doi:10.11606/T.3.2008.tde-04092008-155911.
https://doi.org/10.11606/T.3.2008.tde-04...
], thermal analysis quantifies the heat transmission of the exposed face to the interior of the element by means of the conduction phenomenon, establishing the profile temperatures in the cross-section. In a more complete analysis, the thermal behavior of a structure is also affected by the radiation and convection phenomena. From the temperature profile of the cross-section, it is possible to determine the mechanical strength of the materials, that is, steel and concrete, at elevated temperatures. When a previous temperature profile is gradually applied to a structural element submitted to mechanical loading, a thermo-mechanical response is obtained.

An adequate fire design of the concrete structure is obtained when the following requirements are fulfilled: the connections between structural elements support thermal actions and the heat cannot propagate rapidly along the structure, which prevent failure of other structural elements which are not directly exposed to fire.

When a structure is designed for fire resistance, it implies that the structure, or a part of the structure or a member fulfills its required functions (load bearing function and/or fire separating function) for a specified load level, for a specified fire exposure and for a specified period of time. The designer must ensure that the structure satisfies the fire-resistance class established in standards for each type of element [29[29] KLEIN JÚNIOR, O. Pilares de concreto armado em situação de incêndio submetidos à flexão normal composta. Dissertação (Mestrado em Engenharia de Estruturas), Escola Politécnica, Universidade de São Paulo, São Paulo, 2011. doi:10.11606/D.3.2011.tde-17082011-154210.
https://doi.org/10.11606/D.3.2011.tde-17...
]. In Brazil, the NBR 14432:2001 standard [8[8] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 14432: Fire-resistance requirements for building construction elements - Procedure. Rio de Janeiro, 2001 (in Portuguese).] defines the desired fire resistance time (TRRF) as the fire resistance of a structural element when exposed to heating according to the standard temperature-time curve. According to Silva [30[30] SILVA, V.P. Projeto de estruturas de concreto em situação conforme a ABNT NBR 15200:2012, São Paulo, Blucher, 1 ed., 2012.], the TRRF does not mean the duration of a fire, or time for dispatch or response of the fire brigade, but rather a design parameter. In other words, the TRRF is an equivalent time of fire exposure, that is, the time of exposure to the standard temperature-time curve that is supposed to have the same heating effect as a real fire in the compartment.

2.1 Tabular method

The tabular methods are the simplest, most traditional and are currently used by designers. Because it is a simple way to determine the fire resistance, these methods involve several restrictions and are very safe. For non-tabulated data, interpolations are required to reach the fire resistance. Application of the tabular method is restricted to usual elements with predefined geometry contained in tables.

The tabular method is based on the principle that the fire resistance of the element increases with the increase of the distance between axis of the reinforcement and the face of concrete exposed to fire due to the reduction of the temperature in the reinforcement. Temperatures tend to decrease in the center of the concrete cross-section due to its low conductivity. Therefore, if the thickness of the element and the distance from the axis of the reinforcement to the face exposed to heat increase, more time will elapse until the structural element reaches failure under mechanical loading. In the tabular method, no further check is required for bending, shear or torsion capacity and anchorage details, and the explosive spalling of concrete is not checked.

The Brazilian standard NBR 15200:2012 [7[7] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).] is based on tabulated data presented in Eurocode 2 [3[3] EUROPEAN COMMITTEE FOR STANDARDZATION. Eurocode 2 - Design of concrete structures - Part 1-2: General rules - Structural fire design, 2004.]. The minimum values of the thicknesses of the structural elements and the axis distance of the reinforcement to the face exposed to fire depend on the structural element type and the TRRF, which is a function of the construction type (residential, commercial, industrial, educational, etc.). Table 1 shows the minimum thickness and axis distance (c1) for solid, simply supported, reinforced concrete slabs for standard fire resistance according to NBR 15200:2012 [7[7] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).].

Table 1
Minimum thickness (h) and axis distance from the reinforcement to the face exposed to fire (c1) for simply supported slabs [7[7] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).]

Tabulated data are based on a reference load level for fire design equal to 70% of the load level for ambient temperature. Therefore, the values shown in tables can be reduced for elements where the load level for fire design is less than 70% of the load level for ambient temperature. Table 2 shows the minimum dimensions for simply supported ribbed slabs, according to Brazilian standard NBR 15200:2012 [7[7] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).]. The values of c1 in Tables 1 and 2 are valid for reinforced concrete and 15 mm must be added for prestressed concrete.

Table 2
Minimum dimensions for simply supported ribbed slabs in a fire [7[7] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).]

The latest version of Brazilian standard NBR 9062:2017 [9[9] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 9062: Design and execution of precast concrete structures. Rio de Janeiro, 2017 (in Portuguese).] shows some recommendations for fire design of hollow core slabs by tabulated data. The minimum axis distance for bending capacity (c1) is expressed as a ratio of factored design moment in fire conditions (MSd) and the design moment capacity (MRd) at ambient temperature. Table 3 shows some combinations between slab thickness (h) and minimum axis distance (c1) for each ratio of MSd/MRd. The slab thickness presented in Table 3 can be considered as the sum of the thickness of the hollow core slab and the thickness of the cast-in-place concrete topping. Moreover, these values are true for hollow core slabs with a void index (ratio between the average thickness and the height of the slab) greater than or equal to 0.5.

Table 3
Minimum slab thickness and axis distance for simply supported hollow core floors [9[9] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 9062: Design and execution of precast concrete structures. Rio de Janeiro, 2017 (in Portuguese).]

2.2 500 °C isotherm method

This method was developed in Sweden and is one of the best known to designers. The simplified calculation method comprises a general reduction of the cross-section size with respect to a heat damaged zone at the concrete surface. The thickness of the damaged concrete is made equal to the average depth of the 500 °C isotherm in the compression zone of the cross-section. Moreover, damage concrete, that is, concrete with temperature in excess 500 °C, is assumed not to contribute to the load bearing capacity of the element. The residual concrete cross-section retains its initial values of strength and modulus of elasticity. Reinforcement steel in zones with a temperature greater than 500 °C is not excluded in the analysis, but its mechanical strength is reduced due to the temperature [21[21] AGUADO, J.V., ESPINOS, A, HOSPITALER, A, ORTEGA, J, ROMERO, L. Influence of reinforcements arrangement in flexure fire behavior of hollow core slabs. Fire Safety Journal, v. 53, p. 72-84, 2012.].

Tests carried out with varied mechanical loading and desired fire resistance time (TRRF) confirmed that the simplification of excluding the cross-section for that with temperature greater than 500°C is valid for normal density concretes and siliceous aggregates [28[28] COSTA, C.N. Dimensionamento de elementos de concreto armado em situação de incêndio. Tese (Doutorado em Engenharia de Estruturas), Escola Politécnica, Universidade de São Paulo, São Paulo, 2008. doi:10.11606/T.3.2008.tde-04092008-155911.
https://doi.org/10.11606/T.3.2008.tde-04...
].

According to Klein Júnior [29[29] KLEIN JÚNIOR, O. Pilares de concreto armado em situação de incêndio submetidos à flexão normal composta. Dissertação (Mestrado em Engenharia de Estruturas), Escola Politécnica, Universidade de São Paulo, São Paulo, 2011. doi:10.11606/D.3.2011.tde-17082011-154210.
https://doi.org/10.11606/D.3.2011.tde-17...
], the 500 °C isotherm method can be applied to any fire curve, since the cross-section data (height, width and position of the reinforcement) and the temperature profile are known. This simplified calculation method can be used for structural elements subjected to simple bending or composite bending.

The basic procedures for this calculation method are: determine the temperature profile in the cross-section of the structure for a specific value of TRRF, reduce the cross-section, excluding the concrete outside the 500 °C isotherm, determine the temperature of reinforcing bars in the tension and compression zone, determine the reduced strength of the steel reinforcement due to temperature, determine the ultimate load bearing capacity with the strength of the reinforcing bars using the conventional calculation method for the reduced cross-section. Finally, the estimated fire resistance should be compared with the fire design loading.

3. Computational modeling

DIANA® 9.6 [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.] is a finite element software package that can be applied to simulate typical structures of civil engineering and others engineering sciences (mechanical, biomechanics, etc.). For thermo-mechanical analysis, the DIANA® 9.6 [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.] software has a specific module to process the structures. This is the junction of the nonlinear structural module with the heat flow analysis. Each module has its own finite element library, set of loadings and constitutive law of materials.

According Rocha [32[32] ROCHA, F.M. Modelos numéricos de vigas mistas de aço e concreto pertencentes a sistemas de pisos mistos de pequena altura em situação de incêndio. Dissertação (Mestrado em Estruturas), Escola de Engenharia de São Carlos, Universidade de São Paulo, São Carlos, 2012. doi:10.11606/D.18.2012.tde-19032012-090525.
https://doi.org/10.11606/D.18.2012.tde-1...
], all data should be considered in only one model, because the DIANA® software will automatically convert the results of thermal analysis into data for structural analysis in coupled thermo-mechanical analyses. The thermal model must consist of elements with linear interpolation, while in the structural model the elements can have quadratic interpolation, the software being responsibility for the compatibilization of the elements.

3.1 Constitutive law of materials

A constitutive law of materials to represent the mechanical behavior of materials at ambient temperature should be defined, and the variation of these properties with temperature should also be defined. In this paper, the reduction of the mechanical properties of concrete and steel with temperature is given in Figure 1, which is similar to that recommended in Eurocode 2 [3[3] EUROPEAN COMMITTEE FOR STANDARDZATION. Eurocode 2 - Design of concrete structures - Part 1-2: General rules - Structural fire design, 2004.]. The steel was considered a perfect elastoplastic material with von Mises yield criterion, both at ambient temperature and at elevated temperature.

Figure 1
Reduction of mechanical properties of concrete and steel with temperature: (a) compressive strength of the concrete (siliceous aggregate), (b) modulus of elasticity of concrete, (c) concrete tensile strength, (d) Mode-I fracture energy of concrete, (e) modulus of elasticity of steel, (f) yield strength of prestressing steel

Figure 2
Predefined uniaxial behavior for total strain model [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]: (a) tension softening; (b) compression

The DIANA® 9.6 software [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.] does not have a multilinear stress-strain relationship for compressive and tensile strength of concrete at elevated temperature, as shown in Eurocode 2 for concrete under compression. Therefore, the Total Strain Crack Model with smeared cracking was used to represent the mechanical behavior of the concrete. In this model, the DIANA® 9.6 software [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.] represents the mechanical behavior of concrete by a predefined uniaxial stress-strain relationship, which depends only on the modulus of elasticity and uniaxial strength. In this paper, the exponential tension softening model shown in Figure 2a for concrete under tension and the parabolic model shown in Figure 2b for concrete under compression were adopted. The Mode-I fracture energy (GF), as well as its variation with the temperature, was adopted from research of Rocha [32[32] ROCHA, F.M. Modelos numéricos de vigas mistas de aço e concreto pertencentes a sistemas de pisos mistos de pequena altura em situação de incêndio. Dissertação (Mestrado em Estruturas), Escola de Engenharia de São Carlos, Universidade de São Paulo, São Carlos, 2012. doi:10.11606/D.18.2012.tde-19032012-090525.
https://doi.org/10.11606/D.18.2012.tde-1...
] and Santos [33[33] SANTOS, D.B.R. Modelagem numérica de lajes mistas de aço e concreto em situação de incêndio. Dissertação (Mestrado em Estruturas), Escola de Engenharia de São Carlos, Universidade de São Paulo, São Carlos, 2014. doi:10.11606/D.18.2014.tde-04062014-153311.
https://doi.org/10.11606/D.18.2014.tde-0...
], shown in Figure 1, which presented good results when compared to experimental results of steel-concrete composite beams. The compressive fracture energy (GC), in all analyses, was taken equal to one hundred times the Mode-I fracture energy [34[34] FEENSTRA, P. H., BORST, R. Aspects of robust computational modelling for plain and reinforced concrete. Heron, , v. 38, n. 4, p. 3-76, 1993.]. Thus, the compressive fracture energy and the uniaxial stress-strain relationship of concrete under compression varied with the temperature according to the variation adopted for the Mode-I fracture energy.

3.2 Finite element mesh

Two types of finite elements are required to perform a heat flow analysis in DIANA® 9.6 software [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]: one to represent the volume of the structural element and another for the regions where the thermal loads are applied, also called boundary conditions. The Q4HT finite element for heat flow analysis of hollow core slabs in two dimensions was used. It is a four-node quadrilateral isoparametric element for general potential flow analysis, shown in Figure 3. It is based on linear interpolation and Gauss integration. For boundary conditions, the B2HT element was used. It is a two-node isoparametric straight boundary element for general potential flow analysis. It is also based on linear interpolation and Gauss integration, shown in Figure 4.

Figure 3
Finite element Q4HT [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]

Figure 4
Finite element B2HT [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]

For three-dimensional thermo-mechanical analyses of the hollow core slab, the CHX60 finite element was used. It is a twenty-node isoparametric solid brick element and is based on quadratic interpolation and Gauss integration (Figure 5). For boundary conditions, the BQ4HT finite element was used. It is a four-node isoparametric quadrilateral element to describe boundaries in three-dimensional general potential flow analysis (Figure 6). It is based on linear interpolation and Gauss integration.

Figure 5
Finite element CHX60 [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]

Figure 6
Finite element BQ4HT [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]

In this thermo-mechanical analysis, the selected boundary condition is based on linear interpolation and it is incompatible with the CHX60 finite element. The DIANA® 9.6 software automatically makes the compatibilization of the CHX60 structural element with the BQ4HT finite element [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]. The temperature on the faces of hollow core slab not exposed to fire was defined as 20 °C.

The prestressing strands were modeled by embedding reinforcement elements, which simulate the perfect bond of bar reinforcement and concrete. The prestress load in the model was applied as initial loading condition.

3.3 Boundary conditions and loading

Two models were developed, i.e., a two-dimensional heat flow analysis and a three-dimensional thermo-mechanical analysis. Convection at boundaries of the continuum model and radiation were considered in models developed in DIANA® 9.6 [31[31] TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.]. The concentrated mechanical loading in the thermo-mechanical model was applied as pressure on the region of load applied.

Boundary elements in heat flow analysis are placed on faces to simulate the heat conduction between the structure and the ambient air. For a hollow core slab, the bottom surface is exposed to fire. The top face and side faces are insulated and do not transfer heat to ambient air. Inside the voids, there is heat flow from its lower (hotter) area to the upper (colder) area through convection and radiation flux of the air inside the voids. Thus, the heat flow along the cross-section of the slab takes into account conduction, as well as the convection and radiation inside the voids.

The thermo-mechanical analysis is performed in three steps. In the first step, the initial stress state due to the prestressing strands is imported as the initial condition. This stress state due prestressing is constant over all the analysis, i.e., it is assumed that there is no prestressing loss due to elevated temperature and all stress loss occurs before exposing the slab to fire. The second step consists in applying the mechanical loading, as percentage of the maximum load capacity of the ambient temperature. This loading is applied incrementally by performing a physical nonlinear analysis and is maintained constant during thermal analysis. Thirdly, the temperature profile is imported from the thermal model, and a transient analysis is evaluated. The time-dependent profile induces thermal strain and degradation of material properties. The calculation is stopped when the hollow core slab is so damaged that convergence cannot be achieved. This methodology is similar to previous research developed and reported in the technical literature [20[20] AGUADO, J. V., ALBERO, V., ESPINOS, A., HOSPITALER, A., ROMERO, M. L. A 3D finite element model for predicting the fire behavior of hollow-core slabs. Engineering Structures, v. 108, p. 12-27, 2016.].

4. Validation of the computational model

Throughout this section, both the heat flow and the thermo-mechanical models developed in DIANA® are validates using a previous experimental result obtained by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] for hollow core slabs exposed to standard fire according to ASTM-E119.

4.1 Geometry

Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] tested and modeled six hollow core slabs subjected to service loading and different fire curves. In this paper, one of the slabs tested by those authors was chosen to validate the computational model. The chosen slab was made with a carbonate coarse aggregate type and was exposed to the standard fire curve defined by ASTM E-119 [2[2] ASTM International. ASTM E119-15 - Standard test methods for fire tests of building construction and materials. U.S.A, 2015.]. The slab had span of 4.0 m, width of 1.20 m, and thickness of 20 cm, with six circular voids and seven prestressing strands of diameter 12.7 mm. Figures 7 and 8 show the cross-section of the hollow core slabs and the mechanical loading scheme, respectively. The hollow core slab was made of concrete with compressive strength of 75 MPa, designed according to the PCI Manual [5[5] PRECAST/PRESTRESSED CONCRETE INSTITUTE. PCI Manual for design of hollowcore slabs. 2. ed. Chicago, 1998.]. The supports were free for rotation. A mechanical loading of 69.4 KN was applied to the slab, which corresponded to 60% of the maximum load capacity of the slab in ambient temperature. Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] also developed numerical models in a finite element software to represent the behavior of this hollow core slab when exposed to fire. Further information can be obtained from the reference [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.].

Figure 7
Cross-section of hollow core slab used in validation of the computational model [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]

Figure 8
Test setup of hollow core slab used in validation of the computational model [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]

4.2 Material properties

Table 4 shows the thermal and physical properties of reinforcing steel and concrete used in computational modeling. Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] reported in their work only the physical properties of concrete and steel at ambient temperature, not reporting its strength reduction with temperature. The reduction of the physical properties of concrete and steel with the temperature, used for computational modeling, are shown in Figure 1.

Table 4
Thermal and physical properties1 of concrete and reinforcing steel

4.3 Finite element mesh and modeling strategies

Three different modeling strategies were performed in DIANA® 9.6 to determine the temperature profile of the slab cross-section. In the first strategy, only the bottom face of the slab was exposed to fire. In the second strategy, the air in the void was included in modeling. In this case, the air in the voids was discretized using the air properties shown in Table 5. Finally, the third strategy considered the heat flow from the lower (hotter) area of the void to the upper (colder) area of the void through convection and radiation flux of the air inside the voids (Figure 9). Figure 10 shows the finite element mesh of the hollow core slab with the discretization of the air in the voids.

Figure 9
Finite element mesh with heat flow in voids, from reference [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]

Figure 10
Mesh of hollow core slab with Q4HT elements with air in the voids included in computational model

Table 5
Thermal properties of the air

The first modeling strategy is typically used in several studies, i.e., the hollow core slab was modeled with voids, but without taking into account the convection and radiation flux inside the voids [35[35] PIERIN, I., SILVA, V. P. Análise térmica de laje hollow core protendida em situação de incêndio. In: 51º Congresso Brasileiro do Concreto, 2009, Curitiba. Anais: Concreto para obras de infraestrutura Sustentáveis. São Paulo: IBRACON, 2009. p. 1-17.]. The second modeling strategy adopt a solid slab, but in the region of voids the thermal properties of the air were used to simulate the heat flow inside the voids. The accurate modeling of hollow core slabs needs to correctly reproduce the thermal exchange through the concrete and the voids. This modeling strategy is compared with the numerical modeling performed by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]. This second modeling strategy is similar to that used by other researchers [27[27] VENANZI, I., BRECCOLOTTI, M., D’ALESSANDRO, A., MATEREZZI, A. L. Fire performance assessment of HPLWC hollow core slabs through full-scale furnace testing. Fire Safety Journal, v. 69, p. 12-22, 2014.].

The third modeling strategy is similar to that performed by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] and Aguado et al. [20[20] AGUADO, J. V., ALBERO, V., ESPINOS, A., HOSPITALER, A., ROMERO, M. L. A 3D finite element model for predicting the fire behavior of hollow-core slabs. Engineering Structures, v. 108, p. 12-27, 2016.]. In this strategy, there are boundary conditions in the lower (hotter) and upper (colder) areas of the void to provide a heat flow inside the voids. However, the temperatures inside the voids was not determined during the test, so Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] proposed the use of a temperature curve to be applied to the upper area of the voids, which differs from the temperature curve on the face exposed to fire. In this paper, this temperature curve was defined from the average of the nodal temperature in lower area of the void. Then, this temperature curve was applied to the upper area of the void to simulate the convection and radiation flux inside the voids. This methodology is similar to that used by Aguado et al. [20[20] AGUADO, J. V., ALBERO, V., ESPINOS, A., HOSPITALER, A., ROMERO, M. L. A 3D finite element model for predicting the fire behavior of hollow-core slabs. Engineering Structures, v. 108, p. 12-27, 2016.].

For the validation of the thermal modeling strategies, a three-dimensional thermo-mechanical modeling was performed (Figure 11). The result was compared with the computational modeling performed by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] in ANSYS software. The main results analyzed are the vertical displacement of the slab and its cracking pattern.

Figure 11
Finite element mesh of the hollow core slab for thermo-mechanical modeling

4.4 Thermal response

Temperature profile along the cross-section of the hollow core slab for the three modeling strategies are shown in Figures 12, 13 and 14. For first strategy, with only the bottom face exposed to fire, the maximum temperatures on the bottom face of the slab were 876 °C and 998 °C after 60 and 120 minutes, respectively (Figure 12). The temperatures on the prestressing strand were 271 °C and 504 °C after 60 and 120 minutes, respectively. Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] obtained in their numerical modeling values of 300 °C after 60 minutes and 500 °C after 120 minutes for the temperature of the prestressing strand. The difference between this modeling strategy and the results of Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] is 9.7% for 60 minutes and 8.4% for 120 minutes of slab exposure to standard fire. In the top face of the slab, not exposed to fire, the temperature obtained from first modeling strategy was 42 °C, a value 82.7% lower than that obtained by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.].

Figure 12
Temperature profile in cross-section of hollow core slab - first strategy: a) 60 minutes; b) 120 minutes

Figure 13
Temperature profile in cross-section of hollow core slab - second strategy: a) 60 minutes; b) 120 minutes

Figure 14
Temperature profile in cross-section of hollow core slab - third strategy: a) 60 minutes; b) 120 minutes

For the second strategy, with the air in voids included in the computational model, the temperature profile is shown in Figure 13. The maximum temperatures on the bottom face of the slab after 60 and 120 minutes of slab exposure to fire were 908 °C and 1000 °C, respectively. These values are greater than the temperatures values obtained from the first modeling strategy. The temperature profile showed a wavy shape along the slab, and the maximum temperatures on the prestressing reinforcement were 273 °C and 494 °C after 60 and 120 minutes, respectively. The difference between this modeling strategy and the results of Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] was 5.9% for 60 minutes and 6.2% for 120 minutes of slab exposure to standard fire. This shows that the heat flow inside the void in fact influences the temperature profile, which implies more elevated temperatures in some regions of the cross-section.

The temperature profile along the cross-section of the hollow core slab for the third modeling strategy is shown in Figure 14. Figure 15 shows the temperature-time curve applied to the upper area of the void to represent the heat flow in the voids. The maximum temperature on the bottom face of the slab was the same as that obtained from the first modeling strategy. However, on prestressing reinforcements the temperature increased to 281 °C after 60 minutes and to 553 °C after 120 minutes. The difference between this modeling strategy and the results of Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] was only 3.1% for 60 minutes, but it was 18.9% greater for 120 minutes of slab exposure to standard fire. Again, this shows that the heat flow inside the void there influences the temperature profile, implying more elevated temperatures on prestressing reinforcements.

Figure 15
Temperature-time curve applied to the upper area of voids in hollow core slab with 20 cm thickness and circular voids

For the top face of the slab, the temperatures obtained from the third modeling strategy are closer to the values obtained by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] in their modeling and test. The maximum temperature obtained from modeling in DIANA® 9.6 was 296 °C, a value only 19.5% greater than the value obtained by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.].

Figure 16 shows the variation of temperature with fire exposure time in some points of the cross-section of the slab for first, second and third modeling strategies. These points were monitored in the test of the slab exposed to standard fire. It is noted that the third modeling strategy best represents the variation of temperature in the inner region of the slab tested by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]. For example, the temperature on the top face of the slab in first strategy after 140 minutes was only 42 °C, while in tests of Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] the temperature in this face reached 248 °C. In the second modeling strategy, the temperature on the top face of the slab was only 49 °C too. That is, in both cases there was a difference about 200 °C for the value obtained in tests by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.].

Figure 16
Variation of temperature with fire exposure time: a) first strategy; b) second strategy; c) third strategy

In the same way, the temperature at the middle point of the cross-section of the hollow core slab in the first strategy was 50% lower than the value obtained by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.], while the value obtained from the second strategy was 93% lower. However, the difference between temperature obtained from the third modeling strategy in DIANA® 9.6 and the value obtained by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] was only 7.8%. In this strategy, however, a more elevated temperature of the prestressing reinforcement after 90 minutes of slab exposure to fire was observed when compared to the values obtained from the modeling and tests of Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]. Therefore, third modeling strategy with the temperature curve shown in Figure 15 is more adequate to represent the inner temperature of the cross-section of the hollow core slab, despite the greater temperature observed in the prestressing reinforcements after 90 minutes.

4.5 Thermo-mechanical response

The thermo-mechanical model evaluated, besides the temperature profile in the cross-section, the vertical displacement and cracking pattern of the slab when exposed to standard fire. The hollow core slab was simply supported and restricted to displacement in the horizontal (or longitudinal) direction of the slab.

The prestress of strand reinforcement used in modeling was 1305 MPa, which gives a force of 13.18 kN in each strand with diameter of 12.7 mm. This stress represents 70% of the yield strength of the prestressing steel and was the value applied in the fabrication of the slab [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]. Mechanical loading was represented by means of two concentrated forces of 69.4 kN each, applied according the test setup shown in Figure 8. This loading corresponds to 60% of the ultimate loading of the slab at ambient temperature.

The variation of the vertical displacement in the middle of the span obtained from thermo-mechanical analysis is shown in Figure 17a. This response is compared to the computational simulation developed by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.], and it is observed that the vertical displacement obtained from DIANA® 9.6 software is greater after 40 minutes. This happened even though the temperature in the prestressing reinforcement obtained from the third thermal modeling strategy is the same as the temperature obtained from Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] up to, approximately, 90 minutes. A possible reason for this divergence may be due to the different constitutive laws of concrete used in this paper and that used by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.].

Figure 17
(a) Variation of displacement with fire exposure time; (b) cracking pattern of slab after 60, 90 and 120 minutes

Despite this difference, a small and similar deflection rate up to 90 minutes from the two computational models is observed, i.e., the deflection rate is approximately 0.6 mm/min from the modeling of Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] and 0.9 mm/min that from the modeling in this paper. The loaded slab was assumed to collapse, according to UNI-EN 1363-1 [27[27] VENANZI, I., BRECCOLOTTI, M., D’ALESSANDRO, A., MATEREZZI, A. L. Fire performance assessment of HPLWC hollow core slabs through full-scale furnace testing. Fire Safety Journal, v. 69, p. 12-22, 2014.], after 90 minutes when the deflection rate reached the limiting value given as follows:

v = d z d t = L 2 9000 d ( m m / min ) (1)

In this expression, L is the span of the slab, and d is the thickness of the slab. For the hollow core slab in this analysis, the limit of the deflection rate is v = 8.9 mm/min. This value was not achieved in any of the computational models, which indicates that this slab would resist up to 140 minutes of standard fire without reaching collapse. The ultimate bending moment of the hollow core slab after fire exposure using a simplified design method was evaluated in next section to confirm this hypothesis.

Figure 17b shows the cracking pattern of the hollow core slab after 60, 90 and 120 minutes of fire exposure, obtained from computational modeling. This cracking pattern is very close to that obtained by Kodur and Shakya [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.]. The cracks in the longitudinal direction started in the test after 20 minutes, as well as in this computational model. The bending cracks started at 60 minutes, and at 75 minutes the shear cracks appear.

4.6 Evaluation of flexural resistance by simplified design method

The flexural resistance of a hollow core slab for a specific time of fire exposure, or residual resistance, can be determined from simplified design methods like the 500 °C isotherm method described in Section 2.2. If the temperature profile in the cross-section, including the temperature of the prestressing reinforcements, is known, the reduced strength of prestressing steel at elevated temperature is determined from Figure 1f, as recommended by Eurocode 2 [3[3] EUROPEAN COMMITTEE FOR STANDARDZATION. Eurocode 2 - Design of concrete structures - Part 1-2: General rules - Structural fire design, 2004.]. For each time, the flexural resistance is determined from the classical theory of reinforced and prestressed concrete considering the rupture domains for the ultimate limit state defined by Brazilian standard NBR 6118:2014 [37[37] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 6118: Design of concrete structures - Procedure. Rio de Janeiro, 2014 (in Portuguese).]. The temperature profile used in this analysis was determined from the first and third modeling strategies described in Section 4.4.

Results obtained from simplified design method are compared with the results obtained from computational modeling developed by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] in Figure 18. It is noted that the ultimate bending moment at ambient temperature for this slab was estimated at 160 kNm by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] and at 179 kNm by the simplified method. Then, the ultimate bending moment estimated by the simplified design method is only 12% greater than value obtained from nonlinear mechanical modeling.

Figure 18
Flexural strength of hollow core slab with time of exposure to standard fire

The resistance of the prestressing reinforcement decreases with the increase of the temperature, which results in a decrease of the flexural resistance of the hollow core slab. Unitary values for partial factors for the strength of steel and concrete are adopted, and the reduction of the compressive strength of the concrete due to long-term effects (creep) is not considered since a fire is an exceptional action. It is noted from Figure 18 that the rate of decrease of the flexural resistance of the hollow core slab estimated by the simplified design method is slightly greater than that predicted by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] for up to 80 minutes of standard fire. After this time, the simplified design method with the temperature profile obtained from the first modeling strategy tends to estimate flexural resistance close to values obtained from the nonlinear mechanical modeling, while the third modeling strategy used to determine the temperature profile tends to provide lower flexural resistance.

From Table 3, this hollow core can satisfy a TRRF of 90 to 120 minutes, depending on the ratio between the factored design moment in fire conditions and the design moment capacity at ambient temperature (MSd/MRd). Since this hollow core slab was loaded by 60% of its ultimate loading at ambient temperature, the MSd/MRd ratio of this slab is equal to 0.6, and it can be concluded from Table 3 that this slab would support 90 minutes of standard fire. Taking MSd= 96 kNm, i.e., 60% of the flexural strength at ambient temperature estimated by Shakya and Kodur [26[26] SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.] and 54% of the flexural strength estimated by the simplified design method, it can be concluded from Figure 18 that the collapse of this slab would be reached between 80 and 90 minutes. This is in agreement with the tabulated data presented in Table 3. There is no influence of the thermal modeling strategy in this analysis, since the temperature in the prestressing reinforcement is the same in all strategies until 90 minutes of fire exposure.

If this hollow core slab is required to support 120 minutes of standard fire without collapse, the MSd/MRd ratio from Table 3 should not be greater than 0.49. That is, the factored design moment in fire conditions of the slab after 120 minutes of exposure to standard fire (MSd) would be 78 kNm. The simplified design method indicates that the slab reaches this flexural strength at 90 to 100 minutes when exposed to standard fire (Figure 18). The greater period of exposure was obtained from the first thermal modeling strategy. This shows that the tabulated data, in this case, overestimate the flexural strength of the slab after 120 minutes of exposure to standard fire.

5. Shallow hollow core slab exposed to fire

This section analyzes the performance of a shallow hollow core slab when exposed to standard fire. The hollow core slab taken for this analysis is 16 cm thick, as commonly used in civil construction with small span or lower live loadings. The geometry of the hollow core slab is shown in Figure 19. The test results of this slab under ambient temperature can be found in reference [36[36] SILVA, R.P.M. Resistência à força cortante de lajes alveolares preenchidas de pequena altura, Goiânia, 2015, Dissertação (mestrado) - Escola de Engenharia Civil e Ambiental, Universidade Federal de Goiás, 201p.]. For thermal analysis, the bottom face of this hollow core slab was exposed to standard fire recommended by ISO 834 [38[38] INTERNATIONAL ORGANIZATION STARDADZATION. ISO 834-1:1999 Fire-resistance tests - Elements of building construction - Part 1: General requirements. Geneva, 1999.]. The temperature-time curve of ISO 834 is similar to the temperature-curve of ASTM E-119 [2[2] ASTM International. ASTM E119-15 - Standard test methods for fire tests of building construction and materials. U.S.A, 2015.] used in the validation of the computational model.

Figure 19
Cross-section of shallow hollow core slab - dimensions in mm [36[36] SILVA, R.P.M. Resistência à força cortante de lajes alveolares preenchidas de pequena altura, Goiânia, 2015, Dissertação (mestrado) - Escola de Engenharia Civil e Ambiental, Universidade Federal de Goiás, 201p.]

5.1 Material properties

The thermal properties of the reinforcing steel and concrete used in this model are the same as those described in Section 4.2, shown in Table 4, while the physical properties are taken from reference [36[36] SILVA, R.P.M. Resistência à força cortante de lajes alveolares preenchidas de pequena altura, Goiânia, 2015, Dissertação (mestrado) - Escola de Engenharia Civil e Ambiental, Universidade Federal de Goiás, 201p.]. The prestressing steel is represented by a perfect elastoplastic material with von Mises yield criterion, both at ambient temperature and at elevated temperature. The yield strength of the steel is 1900 MPa. The Total Strain Crack Model with smeared cracking is used to represent concrete behavior. The exponential curve is used in tension softening (Figure 2a), and the parabolic curve is used for compression behavior (Figure 2b). The strength decrease of the material at elevated temperature is similar to that shown in Figure 1.

5.2 Finite element mesh and modeling strategies

Two different modeling strategies were performed to determine the temperature profile of the slab cross-section: firstly with the temperature-time curve of ISO 834 applied only on the bottom face of the slab, and secondly considering the heat flow from the lower (hotter) area of the void to the upper (colder) area of the void. The temperature-time curve to be applied to upper area of the void to simulate the convection and radiation flux inside the voids is defined from the average of the nodal temperature in the lower area of the void. Figure 20 shows the temperature-time curve applied to the upper area of the void for the shallow hollow core slab.

Figure 20
Temperature-time curve applied to the upper area of voids in shallow hollow core slab

Figure 21 shows the mesh in the cross-section of the slab. The finite elements have an approximate size of 5 mm, resulting in 1266 finite elements. The Q4HT finite element is used in the cross-section of the slab, and the B2HT element is used for boundary conditions. The top face and side faces are insulated and do not transfer heat to ambient air.

Figure 21
Cross-section of the shallow hollow core slab with finite element mesh

5.3 Thermal response

The temperature profile along the cross-section of the shallow hollow core slab after 60 minutes is shown in Figure 22. Figure 23 shows the variation of temperature with fire exposure time in the points A (prestressing strand, with a distance of 32 mm from the bottom face), B (middle thickness of the slab), C (40 mm below the top face) and D (top face). It is noted that the temperature in the prestressing strand for the two modeling strategies is similar until 90 minutes. After this time, it is noted that there is a tendency to more elevated temperatures in modeling strategy with the temperature-time curve applied on upper face of the void. The temperature in points B, C and D in this modeling strategy is much higher, reaching 117 °C on the top face of the slab after 60 minutes of slab exposure to standard fire. At this time, the temperature in the prestressing strand reaches 505 °C. This value is more than that provided by Figure 24 and by the other technical literature [39[39] ASSOCIATION OF MANUFACTURERS OF PRESTRESSED HOLLOW CORE FLOORS - ASSAP. The Hollow Core Floor Design and Applications. 1 ed. Verona: Offset Print Veneta. 2002.], where this value of temperature in the reinforcement was reached only after 120 minutes of slab exposure to fire.

Figure 22
Temperature profile after 60 minutes - shallow hollow core slab: a) first modeling strategy; b) second modeling strategy

Figure 23
Variation of temperature with fire exposure time - shallow hollow core slab

Figure 24
Temperature of reinforcement as a function of the distance from the face exposed to fire [6[6] BRITISH STANDARDS. BS EN 1168:2005. Precast concrete products - Hollow core slabs. 2005.]

Mostly, the temperature profile available in the technical literature was obtained from hollow core slabs with thickness greater than 20 cm. When the thickness of this hollow core slab is increased to 20 cm, the temperature in the prestressing strand reaches 380 °C (Figure 25) after 60 minutes of slab exposure standard fire. This temperature value is similar to that predicted by Figure 24. Then, it can be concluded that besides the distance of the reinforcement to the bottom face of the slab, the thickness of the hollow core has great influence on the variation of reinforcement temperature.

Figure 25
Variation of temperature with fire exposure time - hollow core slab with thickness of 20 cm

5.4 Evaluation of flexural resistance

Initially, the flexural strength at ambient temperature of the shallow hollow core slab was determined from nonlinear mechanical modeling, according to the test setup shown in Figure 26. The bending moment versus deflection of the slab is shown in Figure 27. It is observed that the ultimate bending moment obtained from modeling is 102 kNm, while the ultimate bending moment from the test is 99.6 kNm [36[36] SILVA, R.P.M. Resistência à força cortante de lajes alveolares preenchidas de pequena altura, Goiânia, 2015, Dissertação (mestrado) - Escola de Engenharia Civil e Ambiental, Universidade Federal de Goiás, 201p.]. A difference of only 2.4% was observed, which validates the computational modeling. The ultimate bending moment at ambient temperature estimated by the simplified design method is 110.1 kNm, that is, a difference of 10% from the test value was obtained. This proves the ability of the simplified design method to estimate the flexural strength of the shallow hollow core slab in analysis.

Figure 26
Test setup for mechanical testing of shallow hollow core slab [36[36] SILVA, R.P.M. Resistência à força cortante de lajes alveolares preenchidas de pequena altura, Goiânia, 2015, Dissertação (mestrado) - Escola de Engenharia Civil e Ambiental, Universidade Federal de Goiás, 201p.] - dimensions in cm

Figure 27
Load versus deflection for shallow hollow core slab at ambient temperature [36[36] SILVA, R.P.M. Resistência à força cortante de lajes alveolares preenchidas de pequena altura, Goiânia, 2015, Dissertação (mestrado) - Escola de Engenharia Civil e Ambiental, Universidade Federal de Goiás, 201p.]

Figure 28 shows the decrease in the flexural resistance of the slab with the time of exposure to standard fire obtained from the 500°C isotherm simplified design method, as described in Section 2.2. In this analysis, the first thermal modeling strategy was used, since no influence of the modeling strategy at a temperature of reinforcement until 90 minutes was observed. In this figure, the bending moment of 66.1 kNm is also highlighted, which represents 60% of the ultimate bending moment at ambient temperature.

Figure 28
Flexural strength of shallow hollow core slab with time of exposure to standard fire

According to Table 3, this shallow hollow core slab can satisfy a TRRF of 60 minutes for an MSd/MRd ratio greater than 0.6. However, as can be seen in Figure 28, this MSd/MRd ratio is achieved after 35 minutes of slab exposure to standard fire. After 60 minutes of exposure to standard fire, the flexural strength of the slab is only 29% of the ultimate bending moment at ambient temperature. This result shows that the tabulated data, in this case, greatly overestimate the flexural strength of the slab after 60 minutes of exposure to standard fire.

6. Conclusions

The main conclusions of this research are:

  • From the three thermal modeling strategies analyzed, it was concluded that the best representation of the temperature in the cross-section of the hollow core slab with thickness of 20 cm is given when the temperature-time curve is applied to the upper area of the void to simulate the heat flow into the voids. However, this strategy can overestimate the temperature in the reinforcements for exposure to standard fire greater than 90 minutes.

  • In the absence of a more precise criterion for estimating the temperature-time curve to be applied for voids, the methodology based an average temperature measured in the lower area of the voids can be used to estimate the temperature in concrete. Similar methodologies are recommended in other research available in the literature, but these studies do not show the temperature-time curve used in computational modeling of hollow core slabs for comparison with the curves obtained in this paper.

  • The analysis of temperature-time curves obtained in this paper for voids shows a clear delay in heat flow of the shallow slab with thickness of 16 cm with oblong voids when compared to the hollow core slab with circular voids and thickness of 20 cm. However, this delay did not reduce the final temperature in the upper area of the voids after 140 minutes of exposure to fire. This delay is due to the greater distance between the lower face of the voids and the bottom face of the slab in shallow hollow core slabs. This distance is 35 mm in the shallow hollow core slab and 25 mm in the hollow core slab with circular voids.

  • According to the tabular method presented by Brazilian standard NBR 9062, the 20 cm thick hollow core slab (with c1 = 44 mm) should be able to resist 90 minutes of standard fire when the MSd/MRd is equal to or greater than 0.6 and 120 minutes when this ratio ranges from 0.4 to 0.49. From the 500 °C isotherm simplified design method, it is concluded that this slab resists standard fire without collapse until 80 to 90 minutes for MSd/MRd ratio equal to 0.6 and until 90 to 100 minutes for MSd/MRd ratio equal to 0.49. Thus, the tabulated data are suitable for estimating the fire resistance of the hollow core slab for MSd/MRd ratio equal to 0.6, but the tabulated overestimate it for MSd/MRd ratio equal to 0.49.

  • The thermal modeling of the shallow hollow core slab with thickness of 16 cm and c1= 32 mm showed that the temperature of the prestressing reinforcement is greater than the values predicted in the technical literature, regardless of the thermal strategy adopted. For 60 minutes of exposure to standard fire, a temperature of the reinforcement was reached that would be reached only at 120 minutes according to the technical literature. This behavior is due to the lower thickness of the hollow core slab, since the same slab, but with thickness of 20 cm, has a lower temperature of the reinforcement, similar to the values predicted in the literature.

  • According to the tabulated data of Brazilian standard NBR 9062, a shallow hollow core slab with thickness of 16 cm and c1 = 32 mm should be able to resist 60 minutes of standard fire with MSd/MRd ratio greater than 0.6 without collapse. However, from the 500 °C isotherm simplified design method, the ratio of MSd/MRd equal to 0.6 is reached only after 35 minutes of exposure to standard fire. After 60 minutes of exposure to standard fire, the flexural strength is only 29% of the ultimate bending moment of the slab at ambient temperature. This result shows that the tabulated data greatly overestimate the fire resistance of the shallow hollow core slab. The temperature profiles available in the literature were generally obtained from hollow core slabs with thickness greater than 20 cm, so they underestimate the temperature profile in shallow hollow core slabs.

Results from this research show the importance of using a more refined method for estimation of the fire resistance of hollow core slabs in projects, since other relevant parameters, such as geometry of the voids, thickness of the slab and the distance from voids to the bottom face of the slab, not only the distance from the reinforcement to the slab face exposed to fire, influence the temperature of the reinforcement. The 500 °C isotherm simplified design method associated with a temperature profile obtained from a thermal modeling provides more realistic estimates of fire resistance and should, whenever possible, be adopted for the design of hollow core slabs exposed to fire.

7. Acknowledgment

The authors wish to thank the Brazilian National Council of Research and Development (CNPq) for financing this research (Project number 442716/2014-1). They also wish to thank the FAPEG (Foundation for research support of the state of Goiás) for granting a scholarship.

8. References

  • [1]
    EUROPEAN COMMITEE FOR STANDARDIZATION. Eurocode 1 - Actions on structures - Part 1-2: General actions - Actions on structures exposed to fire. 2002.
  • [2]
    ASTM International. ASTM E119-15 - Standard test methods for fire tests of building construction and materials. U.S.A, 2015.
  • [3]
    EUROPEAN COMMITTEE FOR STANDARDZATION. Eurocode 2 - Design of concrete structures - Part 1-2: General rules - Structural fire design, 2004.
  • [4]
    AMERICAN CONCRETE INSTITUTE. ACI 216.1M-07 - Code requirements for determining fire resistance of concrete masonry construction assemblies (metric). U.S.A, 2007.
  • [5]
    PRECAST/PRESTRESSED CONCRETE INSTITUTE. PCI Manual for design of hollowcore slabs. 2. ed. Chicago, 1998.
  • [6]
    BRITISH STANDARDS. BS EN 1168:2005. Precast concrete products - Hollow core slabs. 2005.
  • [7]
    ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 15200: Fire design of concrete structures. Rio de Janeiro, 2012 (in Portuguese).
  • [8]
    ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 14432: Fire-resistance requirements for building construction elements - Procedure. Rio de Janeiro, 2001 (in Portuguese).
  • [9]
    ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 9062: Design and execution of precast concrete structures. Rio de Janeiro, 2017 (in Portuguese).
  • [10]
    SILVA, V.P. Concrete beams fire design. Enhancement of some recommendations of the Eurocode. Ibracon Structures and Materials Journal, v. 4, n. 2, 277-303, 2011.
  • [11]
    SUAZNABAR, J.S., SILVA, V.P. Combined axial and flexural loads in short reinforced concrete columns in fire: ultimate limit state curves using 500°C isotherm method. Ibracon Structures and Materials Journal, v. 11, n. 7, 163-182, 2018.
  • [12]
    CHEREM, M., SILVA, V.P. Determination of the ULS, for columns with small dimensions, under biaxial bending and symmetrical fire conditions. Ibracon Structures and Materials Journal, v. 10, n. 2, 451-476, 2017.
  • [13]
    SILVA, V.P. Fire design of reinforced concrete columns. An alternative to the tabular method presented by the Brazilian standard NBR 15200:2004. Ibracon Structures and Materials Journal, v. 1, n. 4, 331-392, 2008.
  • [14]
    ROCHA, F. M., MUNAIAR NETO, J. Analysis of slim floor beams in fire: emphasis on the concrete constitutive models. Ibracon Structures and Materials Journal, v. 7, n. 1, 158-177, 2014.
  • [15]
    SIMÕES, Y.S., ROCHA, F. M., MUNAIAR NETO, J. Numerical comparison between the thermo-structural behavior of steel and partially encased steel and concrete composite columns in fire. Ibracon Structures and Materials Journal, v. 11, n. 4, 876-901, 2018.
  • [16]
    PIERIN, I., SILVA, V.P. Fire design of composite ribbed slabs. Ibracon Structures and Materials Journal, v. 7, n. 2, 178-207, 2014.
  • [17]
    RIGOBELLO, R., ALMEIDA, S.J.C., MUNAIAR NETO, J., MALITE, M., SILVA, V.P. The Influence of end conditions on numerical models of cold formed steel and concrete composite beams in fire. Ibracon Structures and Materials Journal, v. 3, n. 3, 357-373, 2010.
  • [18]
    HERTZ, K., GIULIANI, L., SØRENSEN, L. S. Fire resistance of extruded hollow-core slabs. Journal of Structural Fire Engineering, v. 8, p. 324-336, 2017.
  • [19]
    KAKOGIANNIS, D., PASCUALENA, F., REYMEN, B., PYL, L., NDAMBI, J. M., SEGERS, E., LECOMPTE, D., VANTOMME, J., KRAUTHAMMER, T. Blast performance of reinforced concrete hollow core slabs in combination with fire: Numerical and experimental assessment. Fire Safety Journal, v. 57, p. 69-82, 2013.
  • [20]
    AGUADO, J. V., ALBERO, V., ESPINOS, A., HOSPITALER, A., ROMERO, M. L. A 3D finite element model for predicting the fire behavior of hollow-core slabs. Engineering Structures, v. 108, p. 12-27, 2016.
  • [21]
    AGUADO, J.V., ESPINOS, A, HOSPITALER, A, ORTEGA, J, ROMERO, L. Influence of reinforcements arrangement in flexure fire behavior of hollow core slabs. Fire Safety Journal, v. 53, p. 72-84, 2012.
  • [22]
    KODUR, V. K. R., SHAKYA, A. M. Factors governing the shear response of prestressed concrete hollowcore slabs under fire conditions. Fire Safety Journal, v. 88, p. 67-88, 2017.
  • [23]
    KODUR, V. K. R., SHAKYA, A. M. Modelling the response of precast, prestressed concrete hollow-core slabs exposed to fire. PCI Journal, v. 78, p. 78-94, 2014.
  • [24]
    MIN, J., DHAKAL, R. P., MOSS, P. J., BUCHANAN, A. H., ABU, A. K. Modelling the fire resistance of prestressed concrete floors using multi-spring connection elements. Journal of Structural Fire Engineering, v. 3, p. 1-18, 2012.
  • [25]
    SHAKYA, A. M., KODUR, V. K. R. Modelling shear failure in precast prestressed concrete hollowcore slabs under fire conditions. Journal of Structural Engineering, v. 143(9), 2017.
  • [26]
    SHAKYA, A.M, KODUR, V.K.R. Response of precast prestressed concrete hollowcore slabs under fire conditions. Engineering Structures, v. 87, p. 126-138, 2015.
  • [27]
    VENANZI, I., BRECCOLOTTI, M., D’ALESSANDRO, A., MATEREZZI, A. L. Fire performance assessment of HPLWC hollow core slabs through full-scale furnace testing. Fire Safety Journal, v. 69, p. 12-22, 2014.
  • [28]
    COSTA, C.N. Dimensionamento de elementos de concreto armado em situação de incêndio. Tese (Doutorado em Engenharia de Estruturas), Escola Politécnica, Universidade de São Paulo, São Paulo, 2008. doi:10.11606/T.3.2008.tde-04092008-155911.
    » https://doi.org/10.11606/T.3.2008.tde-04092008-155911
  • [29]
    KLEIN JÚNIOR, O. Pilares de concreto armado em situação de incêndio submetidos à flexão normal composta. Dissertação (Mestrado em Engenharia de Estruturas), Escola Politécnica, Universidade de São Paulo, São Paulo, 2011. doi:10.11606/D.3.2011.tde-17082011-154210.
    » https://doi.org/10.11606/D.3.2011.tde-17082011-154210
  • [30]
    SILVA, V.P. Projeto de estruturas de concreto em situação conforme a ABNT NBR 15200:2012, São Paulo, Blucher, 1 ed., 2012.
  • [31]
    TNO. DIANA: Finite Element Analysis - User’s Manual i DIANA. Release 9.6, 1 ed., Netherlands: TNO DIANA, 2015.
  • [32]
    ROCHA, F.M. Modelos numéricos de vigas mistas de aço e concreto pertencentes a sistemas de pisos mistos de pequena altura em situação de incêndio. Dissertação (Mestrado em Estruturas), Escola de Engenharia de São Carlos, Universidade de São Paulo, São Carlos, 2012. doi:10.11606/D.18.2012.tde-19032012-090525.
    » https://doi.org/10.11606/D.18.2012.tde-19032012-090525
  • [33]
    SANTOS, D.B.R. Modelagem numérica de lajes mistas de aço e concreto em situação de incêndio. Dissertação (Mestrado em Estruturas), Escola de Engenharia de São Carlos, Universidade de São Paulo, São Carlos, 2014. doi:10.11606/D.18.2014.tde-04062014-153311.
    » https://doi.org/10.11606/D.18.2014.tde-04062014-153311
  • [34]
    FEENSTRA, P. H., BORST, R. Aspects of robust computational modelling for plain and reinforced concrete. Heron, , v. 38, n. 4, p. 3-76, 1993.
  • [35]
    PIERIN, I., SILVA, V. P. Análise térmica de laje hollow core protendida em situação de incêndio. In: 51º Congresso Brasileiro do Concreto, 2009, Curitiba. Anais: Concreto para obras de infraestrutura Sustentáveis. São Paulo: IBRACON, 2009. p. 1-17.
  • [36]
    SILVA, R.P.M. Resistência à força cortante de lajes alveolares preenchidas de pequena altura, Goiânia, 2015, Dissertação (mestrado) - Escola de Engenharia Civil e Ambiental, Universidade Federal de Goiás, 201p.
  • [37]
    ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS. NBR 6118: Design of concrete structures - Procedure. Rio de Janeiro, 2014 (in Portuguese).
  • [38]
    INTERNATIONAL ORGANIZATION STARDADZATION. ISO 834-1:1999 Fire-resistance tests - Elements of building construction - Part 1: General requirements. Geneva, 1999.
  • [39]
    ASSOCIATION OF MANUFACTURERS OF PRESTRESSED HOLLOW CORE FLOORS - ASSAP. The Hollow Core Floor Design and Applications. 1 ed. Verona: Offset Print Veneta. 2002.

Publication Dates

  • Publication in this collection
    01 June 2020
  • Date of issue
    Mar-Apr 2020

History

  • Received
    21 Dec 2017
  • Accepted
    15 May 2019
  • Published
    26 Mar 2020
IBRACON - Instituto Brasileiro do Concreto Instituto Brasileiro do Concreto (IBRACON), Av. Queiroz Filho, nº 1700 sala 407/408 Torre D, Villa Lobos Office Park, CEP 05319-000, São Paulo, SP - Brasil, Tel. (55 11) 3735-0202, Fax: (55 11) 3733-2190 - São Paulo - SP - Brazil
E-mail: arlene@ibracon.org.br