Numerical analysis of prestressed hollow core slabs under long term loading

Análise numérica de lajes alveolares protendidas submetidas a carregamentos de longa duração

Numerical analysis of prestressed hollow core slabs under long term loading

Análise numérica de lajes alveolares protendidas submetidas a carregamentos de longa duração

S. R. Pereira; J. M. Calixto; T. P. Bortone

Escola de Engenharia - Universidade Federal de Minas Gerais. Belo Horizonte, MG. Brasil 31270-090. [e-mail], [e-mail], [e-mail]

ABSTRACT

This paper presents a numerical analysis of prestressed hollow core slabs under long term loading. The model considers the time dependence of material and rheological properties in order to predict the actual stage of displacements, strains and stresses. It also takes into account load changes. For the analysis, each slab is divided in a finite number of bar elements, in which the cross section is described in concrete elements, parallel to the flexural axis, and prestressed steel elements. For the results evaluation, the effective concrete area is considered. The numerical results are compared with experimental tests performed on two series of prestressed hollow core slabs. Each series had a different geometry, rate and distribution of prestressing strands. Mid-span displacements were evaluated up to 127 days after initial loading. Good correlation was achieved with both series at and below the service load level.

Keywords: numerical analysis, prestressed concrete, hollow core slabs, long term loading.

RESUMO

Este artigo apresenta os resultados de uma análise numérica de lajes alveolares protendidas extrudadas submetidas a carregamentos de longa duração. O modelo numérico considera os efeitos do tempo nas propriedades mecânicas dos materiais bem como nas propriedades reológicas, de modo a prever flechas, deformações e tensões. O modelo leva em conta também variações no carregamento. Para a análise, cada laje é subdividida em um número finito de elementos de barra cuja seção transversal é discretizada em elementos de concreto paralelos ao eixo de flexão, e elementos de armadura passiva. Para avaliação dos resultados é considerada a área líquida de concreto. Os resultados do modelo numérico foram comparados com valores medidos em ensaios de duas séries distintas de lajes alveolares protendidas. Cada série de lajes possuía diferentes geometrias, taxa e distribuição de cordoalhas. Deslocamentos no meio do vão de cada laje foram avaliados até a idade de 127 dias após o carregamento inicial. Boa correlação foi obtida com o modelo numérico para cargas menores ou iguais as de serviço.

Palavras-chave: análise numérica, lajes alveolares protendidas, carregamentos de longa duração.

1. Introduction

A hollow core slab is a precast, prestressing concrete member with continuous voids provided to reduce weight and cost. They are primarily used as a floor deck system in residential and commercial buildings as well as in parking structures because they are economical, have good fire resistance and sound insulation properties, and are capable of spanning long distances with relatively small depths. Hollow core slabs can make use of prestressing strands, which allow slabs with depths between 150 and 260 mm to span over 9 meters.

When used in buildings, several hollow core slabs are placed next to each other to form a continuous floor system. The small gap that is left between each slab is usually filled with a non-shrink grout. To give the floor a smooth finished surface, a topping slab overlay is poured on the top surface of the hollow core slabs. This topping slab is typically 5 cm deep. The voids in a hollow core slab may be used for electrical or mechanical runs. For example, routing of a lighting circuit through the cores can allow fixtures in an exposed slab ceiling without unsightly surface mounted conduit.

Structurally, a hollow core slab provides the efficiency of a prestressed member for load capacity, span range, and deflection control. In addition, a basic diaphragm is provided for resisting lateral loads by the grouted slab assembly provided proper connections and details exist (PCI [1]).

The basic manufacturing method currently in use for the production of hollow core slabs is a dry cast or extrusion system where a very low slump concrete is forced through the machine. The cores are formed with augers or tubes with the concrete being compacted around the cores.

In this scenario, the objective of this paper is to present a numerical analysis of prestressed hollow core slabs under long term loading. The modelling considers not only the load history of hollow core slabs but also the effects of time in the mechanical properties of the concrete. The slabs are divided in a number of elements each having different material properties as well as load histories. The results of the numerical analysis are compared to experimental tests performed on two series of prestressed hollow core slabs. Each series had a different geometry and rate and distribution of strands as well as material properties. Loads and mid-span displacements were evaluated up to 127 days after initial loading.

2. The numerical model

2.1 Assumptions

In the development of the numerical modelling [2], the following assumptions are introduced:

1. The hollow core slab, the applied loads and the deformations lie in a plane; the plane of loads is a plane of symmetry for the slab.

2. The slab is slender, that is, its length is much larger than its lateral dimensions.

3. Transverse and longitudinal displacements are infinitesimal.

4. Only normal strains parallel to the axis of the slab are considered.

5. The geometry of the slab can vary with respect to time as well along the length.

6. The material properties, in each section, can be different in each element.

7. The effect of shrinkage, creep and cracking as well as the evolution of the mechanical properties of the concrete with time are taken into account.

8. The losses due to steel relaxation in the strands are also considered.

9. The applied loads vary along time.

2.2 Formulation

In the derivation of the model, the slab is divided in a number of sections, each of one composed of elements. These elements can have different geometrical and material properties as well load histories. It is assumed that at a generic time t_u the total strain in a concrete element i is known and given by:

The total strain in this case includes the effects of the applied loads as well creep and shrinkage. These effects correspond to the first, second and third terms respectively of the right hand side in the equation above. The expression for the shrinkage strain and creep coefficient as well as the concrete tangent modulus of elasticity at time τ are the ones presented in NBR 6118 [3]. If a generic time interval (t_u+1 – t₀) is divided in a number of time steps, the total strain of each concrete element i in each section of the slab becomes:

From equation (2) it can be seen that each concrete element can have different material and load history at every time step. It is important to point out that the creep effects are considered to the time step (t_u+1) and p(t_u+1,t_u+1) = 0.

In equation (2) the incremental stress in each concrete element Δσ_ci (t_u+1,t_u) can be explicitly determined if one rewrites the equation after splitting the second summation. Thus:

With this incremental stress, the final stress in each concrete element at time t_u+1 can be obtained as:

This final stress, in each time step, has to be compared to the concrete strength. For case of concrete in tension the limit is defined by the concrete tensile strength f_ctm. If this tensile strength is exceeded, the contribution of this element is neglected.

For the stress evaluation in each element of a prestressed strand is necessary to compute the losses of stress due to the combination of the concrete shrinkage and creep effects as well as the steel relaxation. According to Trevino and Ghali [4], the stress loss due to the steel relaxation in a strand element i is given by:

where:

Thus, the stress loss in a strand element i including the concrete shrinkage and creep effects can be evaluated by:

with

In this case, Δσ_pi(c+s+r) should be determined by iterative process. For the first attempt, consider ξ_r = 1

Then, the stress and strain at a strand element i at time t_u+1 are given by:

For each strand element i, the maximum tensile stress is equal to the steel yield strength.

The internal axial force and bending moment in each cross section are determined as follows:

The equilibrium between external and internal forces in every cross-section of the hollow core slab is given by:

The correct state of strain and the corresponding equilibrium position of the slab at each time step are obtained when both equations (13) and (14) are satisfied at every cross-section. If this is not achieved a new strain state is introduced and the procedure is repeated.

To ensure convergence of the iterative search process for the state of strain which satisfies the equations (13) and (14), one must consider the compressive strength of concrete and steel stress unlimited. The value of these stresses should be checked at the end of the analysis. Stresses levels above usual values for service state level indicate insufficient cross-section dimensions.

After the equilibrium conditions are satisfied at each cross section, the hollow core slab deflected shape can be determined. With the strain of each concrete (equation 2) and strand (equation 10) element in every section, the curvature in each section along the slab is calculated. In this case the transverse displacement curvature relationship is given by:

The rotations are obtained by integrating the above equation. Thus:

The transverse displacements are calculated then by:

Since each slab has been divided in sections and considering the number of the first section equal to 1, the above equation is replaced by:

It is important to point out that equation (18) has to satisfy the slab cinematic conditions.

3. Experimental program

In order to verify the numerical model an experimental program was carried on. It consisted on the testing of two series of hollow core slabs. Each series had three specimens. In series 1, the hollow core slabs were 1245 mm in width and 260 mm in depth with constant hollow cores of 199 mm in diameter. They were 10 m long and pre-tensioned with 4 seven-wire low relaxation strands (12.7 mm in nominal diameter) and 6 seven-wire low relaxation strands (11.1 mm in nominal diameter). Both strands were placed with a 30 mm bottom cover. Slabs of series 2 had the same width and length as of series 1 but with a total depth of 210 mm; the hollow cores were 152 mm in diameter. They had 8 seven-wire low relaxation strands of 12.7 mm in nominal diameter placed with a 30 mm bottom cover. Both series were pre-tensioned to 75% of the ultimate strength fpu = 1900 MPa, before the concrete was cast using dry mix extrusion procedures. The applied prestressing force was transferred to the concrete 20 hours after casting. The details of the slabs cross-sections are shown in figure 1.

The concrete was produced using Brazilian type V cement and limestone as coarse aggregate. Its concrete compressive strength was evaluated employing the impact hammer. The impact hammer is a practical method for determining the concrete strength of slabs cast using dry mix extrusion procedures and considerably more precise than test cylinders. This is because the compaction of the machine cannot be accurately duplicated in making the cylinders. The concrete strength was evaluated at the age of 30 days for slabs of series 1 and 35 days for series 2. These ages correspond to the initial testing time in which a uniform distributed live load was applied to each slab. The average concrete compressive strength at this age was 40.6 and 28.3 MPa for series 1 and 2 respectively. Both of these values correspond to an average of 9 readings taken with the impact hammer in the horizontal position on the side of each slab. The correlation to compressive cylinder strengths was made with the use of curves printed on a plate attached to the instrument. This compressive strength was the only concrete mechanical property measured throughout the study.

The hollow core slab tests were conducted in an open area at Precon Industrial in Pedro Leopoldo, Minas Gerais. The test set-up is shown in figure 2. The specimens were subjected to a uniformly distributed load on a simple span. The clear span between the supports was 9.9 m Cement bags, with a nominal weight of 0.5 kN each, were used for the uniformly distributed load (figure 2b). A dial-gage placed at mid-span was used to measure the deflection.

The uniform distributed live load was applied incrementally after each slab had been placed over supports. During the first four days, loads were incrementally applied daily. Each live load increment corresponded to 1.52 and 1.01 kN/m for series 1 and 2 respectively. These loads were kept for approximately 114 days. After that period additional loads were applied up to each slab series service live load (8.875 and 6.65 kN/m for series 1 and 2 respectively). Before each live load increment and during the period of 127 days, mid-span deflections were measured. These readings were always taken early in the morning to avoid temperature effects.

The temperature and relative humidity were measured daily (also early in the morning) during the testing period. Their average values were calculated and correspond to 25°C and 40% respectively.

4. Comparative study

A comparison of the numerical and experimental results is presented herein. The mid-span deflection for each series was selected for this study since it represents the overall behavior the hollow core slabs. It is important to point out that the comparison is carried out during a 4-month loading period in which the values of the loads were at or below service level.

In the numerical model, the concrete properties including shrinkage and creep effects were derived from the compressive strength measured at the initial testing time and from the temperature and relative humidity average values shown before. For a proper evaluation of the needed parameters, the formulation prescribed in NBR 6118 sections 8.2.8, 12.3.3, A.2.2.3 and A.2.3.2 was used. For each slab, 10 cross-sections along the span were analyzed at each time step with 50 concrete elements in every cross-section. The number of strand elements in every cross-section corresponded to the number of strands for each slab series: 10 for series 1 and 8 for series 2.

The results of the comparative study are presented both in tabular and graphical form. Table 1 shows the results for slabs of Series 1 while the values correspondent to slabs of Series 2 are in table 2. The comparison between the actual and predicted overall behavior can be quantified based on the model error y defined as the ratio between the evaluated mid-span deflection and the measured value. Values of y (last columns of tables 1 and 2) smaller than the unity indicate a stiffer behavior predicted by the numerical model. Each graph, on the other hand, shows an overall visual comparison.

The comparative study for slabs of Series 1 is well represented in figure 3 which shows the load versus the midspan deflection along 119 days for specimen A1. These results reveal a less stiff behavior predicted by the numerical model throughout the load spectrum investigated. The overall comparison for the slabs of this Series is made based on the statistical analyses of y which include its average (m), the standard deviation (s) as well as its coefficient of variation (COV). The average of y is used as a measure of the conservativeness of the numerical modeling and the coefficient of variation is taken as an indication of its accuracy. These values are presented in table 3. The analysis of these results shows that the model predicts conservatively and accurately the measured midspan deflections with a model error average of 1.14 and COV of 12.6 % for all the slabs of this series.

The load versus the midspan deflection along 126 days for Series 2 (slab B2) is presented in Figure 4. It can be seen that the numerical model predicts well the observed behavior. The statistical analysis (table 3) for this series also shows that the numerical procedure evaluates conservatively the measured midspan deflections. On the other hand the model was less accurate in this case since the COV was equal to 25 % for all the slabs of this series.

The overall analysis of the comparative study indicates the conservative bias of the numerical procedure: model error average equals to 1.10 for all the slabs investigated. With respect to accuracy, the modeling also shows good results with an overall COV of 20 %.

The midspan deflection, temperature and relative humidity measurements were taken daily at the same hour (early in the morning) during the investigation campaign. This is an important aspect of the research and it may be one of the reasons for the good correlation achieved between the experimental and modeling results.

5. Conclusions

A numerical model for the analysis of prestressed hollow core slabs under long term loading has been presented. The model considers explicitly the geometrical changes in cross section as well as the time dependence of the loads and of material properties in order to predict the actual stage of displacements, strains and stresses. The numerical results were compared with experimental tests performed on two series of prestressed hollow core slabs. Each series had a different geometry, rate and distribution of prestressing strands as well as material properties. Loads and mid-span displacements were evaluated up to 127 days after initial loading. The overall analysis of the comparative study indicates the conservative bias of the numerical procedure: model error average equals to 1.10 for all the slabs investigated. With respect to accuracy, the modeling scheme also shows good results with an overall COV of 20 %.

6. Acknowledgements

The authors would like to thank Universidade Federal de Minas Gerais (UFMG) and PRECON Industrial for the financial and infrastructural support.

7. References

8. Notation

The following symbols are used in this paper:

ε_ci (t_u, t_o) = total strain in a concrete element iin time t_u;

ε_ci (t_u+1, t_o) = total strain in a concrete element i in time t_u+1;

ε_cci (t, t_o) = creep strain in a concrete element i in time t;

ε_csi (t, t_o) = shrinkage strain in a concrete element i in time t;

ε_pi (t_u+1, t_p) = strain in a strand element i after initial losses;

= total strain in a strand element i in time t_u+1;

E_c(τ) = concrete tangent modulus of elasticity at time τ (NBR 6118 sections 12.3.3 and 8.2.8);

E_c28 = concrete tangent modulus of elasticity at time 28 days (NBR 6118 section 8.2.8);

E_s = steel modulus of elasticity ;

f_ctm = average concrete tensile strength;

φ (t_u,τ) = creep coefficient (NBR 6118 section A.2.2.3);

ε_cs∞ = maximum shrinkage strain (NBR 6118 section A.2.3.2);

t₀ = initial time in which the stress state changes in a cross section;

t_p = time of the prestressing force transfer to the concrete slab;

t_s = time corresponding to the end of the curing period;

β_s(t_u) = function that describes the development of shrinkage strain along the time; (NBR 6118 section A.2.3.2);

Δσ_ci(t_j,t_j-1) = incremental stress in a concrete element i between times t_jand t_j-1;

Δσ_pri(t,t₀) = stress loss in a strand element i due to steel relaxation between times t and t₀;

λ = stress intensity ratio in a prestressed strand;

η_∞ = steel relaxation coefficient;

Δσ_pr∞ = maximum stress loss in a strand due to steel relaxation;

Ω = error coefficient of iteration process;

ξ_r = steel relaxation correction coefficient;

Δσ_pi(c+s+r) = stress loss in a strand element i due to the combination of concrete shrinkage and creep effects as well as steel relaxation;

A_ci = area of the concrete element i;

A_pi = area of the strand element i;

Y_ci = distance from the center of gravity of concrete element i to the top fiber of the slab cross-section;

Y_pi = distance from steel element i to the top fiber of the slab cross-section;

α_pi = inclination angle of strand element i in a slab cross-section;

M_R (t_u+1, t₀) = internal bending moment at time t_u+1;

N_R (t_u+1, t₀) = internal axial force at time t_u+1;

N_EXT = external axial force in a slab cross-section;

M_EXT = external bending moment in a slab cross-section;

r (t_u+1) = radius of curvature of a generic cross-section at time t_u+1;

r_j (t_u+1) = radius of curvature of cross-section j at time ;

θ(t_u+1) = rotation of a generic cross-section at time ;

θ_j (t_u+1) = rotation of cross-section j, at a distance j₁ from the origin, at time t_u+1;

y (t_u+1) = transverse displacement of a generic cross-section at time t_u+1;

y_j (t_u+1) = transverse displacement of cross-section j, at a distance j₁ from the origin, at time t_u+1;

Δℓ _{i, i+1} = longitudinal distance between cross sections i and i+1.

Received: 24 Feb 2013

Accepted: 23 Jul 2013

Available Online: 12 Aug 2013

• [01] PRECAST/PRESTRESSED CONCRETE INSTITUTE, PCI Manual for the Design of Hollow Core Slabs - 2nd edition; Chicago, Illinois, 1998.
• [02] PEREIRA, S. S. R. Análise do Comportamento Reológico de Pontes em Balanços Sucessivos, Doctoral Thesis, COPPE - Universidade Federal do Rio de Janeiro, February, 1999.
• [3] ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS, NBR 6118 – Projeto de Estruturas de Concreto - Procedimento, Rio de Janeiro, 2007.
• [4] TREVINO, J. and GHALI, A. Relaxation of Steel in Prestressed Concrete, Prestressed Concrete Institute Journal, vol. 30, n. 5, September-October, 1985, p. 82-94.

Publication Dates

• Publication in this collection
  27 Aug 2013
• Date of issue
  Aug 2013

History

• Received
  24 Feb 2013
• Accepted
  23 July 2013