Extension of the Zone Method of Eurocode 2 for reinforced concrete columns subjected to standard fire

Purpose – The purpose of this paper is to develop a method suitable for the design of reinforced concrete columns subjected to a standard fire. Design/methodology/approach – The Zone Method – a ’simplified calculation method” included in Eurocode 2 – has been developed by Hertz as a manual calculation scheme for the check of fire resistance of concrete sections. The basic idea is to disregard the thermal strains and to calculate the resistance of a cross-section by reducing the concrete cross-section by a “damaged zone”. It is assumed that all fibers can reach their ultimate, temperature dependent strength. Therefore, it is a plastic concept; the information on the state of strain is lost. The calculation of curvatures and deflections is thus only possible by making further assumptions. Extensions of the zone method toward a general calculation method, suitable for the implementation in commercial design software and using the temperature dependent stress–strain curves of the Advanced Calculation Method, have been developed in Germany. The extension by Cyllok and Achenbach is presented in detail. The necessary assumptions of the Zone Method are reviewed, and an improved proposal for the consideration of the reinforcement in this extended Zone Method is presented. Findings – The principles and assumptions of the Zone Method proposed by Hertz can be validated. Originality/value – An extension of the Zone Method suitable for the implementation in design software is proposed.


Introduction
The Zone Method for the check of fire resistance of reinforced concrete cross-sections is part of EN 1992EN -1-2 (2004b)).It was developed by Hertz (1985) as a scheme for a manual calculation.The concept is to reduce the effective dimensions of the cross-section to simulate the mean deterioration of the concrete.It is assumed that the remaining cross-section can be modeled with a constant concrete strength.Therefore, simplifications and assumptions are necessary: thermal strains are disregarded; the resistance of the cross-section is calculated with a reduced cross-section and simplified material properties; and the ultimate buckling load of a column is determined with Rankine's formula.
The description of the Zone Method in EN 1992-1-2 (2004b) is incomplete, especially for concrete columns.It is not possible to calculate the fire resistance without background information, as given by Hertz (1985).Hence, the use of the Zone Method in Germany is only allowed by building authorities, if the comments and interpretations by Zilch et al. (2010) or Cyllok and Achenbach (2009) are considered.Both papers are listed as non-contradictory complementary information (NCI) in the German National Annex (DIN EN 1992-1-2/NA, 2010).
These interpretations develop the Zone Method to a general method similar to EN 1992EN -1-1, 5.8.6 (2004a) ) with a reduced cross-section and temperature dependent material properties.Zilch, Müller and Reitmayer propose to consider the effect of thermal strains by shifting the stress-strain curves for steel and concrete.Cyllok and Achenbach suggest to reduce the strength of the compressed reinforcement.
The intention of both advancements is to establish an extended Zone Method, which can be easily implemented in commercial software for the design of concrete columns.In contrast to the Advanced Calculation Method, the Zone Method does not need a detailed temperature distribution of the cross-section for the mechanical analysis, as the concrete zone is modeled with constant temperature and material properties.This makes it easy to adjust commercial software for the design of the fire resistance, a coupling of the thermal and mechanical analysis of the cross-section is not required.For example, the temperatures for the centroid and the reinforcement at a given fire resistance can be taken from figures displayed in EN 1992-1-2 (2004b), Annex A.
The benefit of an extended Zone Method is its suitability for the "hot" design of reinforced concrete columns.This method follows the same principles which are used for the design at room temperature, mainly a concrete cross-section with constant material properties.The purpose of an extended Zone Method is the determination of the required area of reinforcement.
For a calculation with the Advanced Calculation Method the area of reinforcement must be known a priori.The result of the calculation is a fire resistance time for a given heating curve or an ultimate load for a fixed time step of a given heating curve.So it is suitable for the recalculation of laboratory tests or existing structures.
In Section 2, the basic assumptions and principles of the Zone Method are presented.Section 3 shows the extensions by Cyllok and Achenbach and gives information on recalculation of fire tests.In Section 4, the assumptions of Section 2 are reviewed.A proposal for the extension of the Zone Method to a general calculation method is given in Section 5. Finally, one worked example is presented in Section 6. Hertz (1985) assumes that the thermal strains of heated concrete cross-sections can be recovered by transient thermal strains.The explicit creep model by Anderberg and

Extension of the Zone Method of Eurocode 2
Thelandersson (1976) is used by Hertz for proof.The transient thermal strains tr are described here as follows: where ϭ stress; f ck ϭ strength of the concrete at 20°C; and th ϭ thermal strain.At a stress level ||/f ck Ն 0.43 thermal strains can be fully recovered by transient thermal strains.Therefore, Hertz concludes that the thermal strains can be disregarded.
In Figure 1, a concrete wall, heated on both sides, is displayed.It is assumed that the isotherms are parallel to the surface, and that all fibers achieve their ultimate strength.The resistance of the concrete wall per unit length N R,c,fi is as follows: where h ϭ height of the wall; k c () ϭ f c, /f ck ; and f c, ϭ concrete strength at temperature .The parameter k c,m denotes the mean strength of the concrete: The reduced mean strength of the concrete is represented by a minimized cross-section with an equivalent resistance.The strength at the center M is used for reference: where M ϭ temperature at the center M.
The relative height f of the reduced cross-section is calculated by: The width of the "damaged" zone a z,f is given by: The resistance of the reinforcement N R,s,fi is calculated using the strength at 0.2 per cent proof strain: where A s ϭ total area of reinforcement; k s,0.2% ϭ f s,0.2%, /f yk ; f s,0.2%; ϭ reinforcement strength at temperature for 0.2 per cent proof strain; and f yk ϭ reinforcement yield strength at 20°C.It must be mentioned that Hertz (1985) uses simplified equations for the reduction factors k c () and k s,0.2% ().The factors given in EN 1992EN -1-2 (2004b) ) are not based on simplifications -they are derived from the stress-strain curves supplied for the Advanced Calculation Method.
Hertz states that the peak strain c1, corresponding to f c, can be approximated by Խ c1,Խ Ϸ 3.5%/k c (). Reducing the strength of the concrete by k c () and increasing the peak strain c1, by 1/k c () leads to the following: where E c ϭ modulus of elasticity at 20°C.Under the assumption that all fibers are under constant compression, the mean stiffness of the concrete wall EI c,fi can be calculated by: The reduced mean stiffness is represented by an equivalent cross-section with the material properties at the center M: where M ϭ temperature at the center M.
The relative height of the cross-section with equivalent stiffness EI is given by: The "damaged" zone a z,EI for a cross-section of equivalent stiffness is given by:

Extension of the Zone Method of Eurocode 2
The calculation of equation ( 11) can be laborious and is not suitable for a manual calculation scheme.Therefore, Hertz proposes a simplified empiric equation: Equation ( 12) can be rewritten using equations ( 13) and ( 5) as follows: According to Hertz, the reduction of the strength of the reinforcement is proportional to the reduction of the modulus of elasticity: where E s ϭ modulus of elasticity at 20°C.For a centrically loaded wall with pin ended conditions, the Euler buckling load for concrete N b,c,fi and steel N b,s,fi can be calculated as follows: where l col ϭ length of the wall; EI c,fi ϭ bending stiffness according to equation (10); z s ϭ lever arm of reinforcement; and A s ϭ total area of reinforcement.
The ultimate load N u,fi of the column according to Hertz is given by Rankine's formula: Where N R,c,fi and N b,c,fi are both calculated using a z,EI as "damaged" zone.

Extension by Cyllok and Achenbach
Richter and Zehfuß (1999) compared the results of the Zone Method with those of the Advanced Calculation Method for bending members.The calculated moment resistance of the Zone Method has always been lower compared to the Advanced Calculation Method.Therefore, the use of the Zone Method for bending members is allowed according to the German National Annex to EN 1992-1-2 (DIN EN 1992-1-2/NA, 2010).The use of the Zone Method for concrete compression members has been discussed controversially in Germany.The description of the Zone Method for columns in the Eurocode (EN 1992(EN -1-2,2004b) is incomplete, which led to misunderstandings and misinterpretations.In particular, it has not been considered that Rankine's formula -as given in equation ( 18) -is used for the calculation of the ultimate load in the Zone Method.JSFE 7,2 Cyllok and Achenbach (2009) examined three interpretations of the Zone Method and made a proposal for an extension.The aim and purpose of this extension was the development toward a general method according to EN 1992-1-1 (2004a) with only slight modifications of the Zone Method as described in Section 2. All research done for bending members, e.g. from Richter and Zehfuß (1999), should still be valid for an extended Zone Method.
The proposed calculation scheme for a square concrete column is represented in Figure 2, left hand side.The assumptions are as follows: • Thermal stresses are neglected.
• The cross-section is reduced by a z,EI .
• The concrete is represented with constant material properties from the Advanced Calculation Method at the temperature M .
• The peak strain of the concrete Խ c1,Խ is at least 3.5‰.
• The reinforcement is considered using the material properties from the Advanced Calculation Method.• The strength of the compressed reinforcement is reduced by ⌬ ϭ 2‰ • E s ( s ).
The Assumptions 1 to 5 are in accordance with Hertz because the material properties for simple calculation methods supplied in the Eurocode (EN 1992(EN -1-2 (2004b))) are derived from the stress-strain curves for the Advanced Calculation Method.The effect of Assumption 6 compared to the Advanced Calculation Method has been examined by Cyllok and Achenbach (2009) for three columns.The shapes of the calculated moment-curvature diagrams are comparable to those of the Advanced Calculation Method.The ultimate load for a given fire resistance deviates ϩ20/Ϫ18 per cent from the Advanced Calculation Method.
Laboratory tests carried out at TU Braunschweig and documented in the thesis of Haß (1986) have been recalculated by Cyllok and Achenbach (2011) using this extended Zone Method.The mean ratio of calculated time to failure to experimental time to failure ϭ t cal /t exp is 0.904 with a standard deviation ϭ 0.218.Haß (1986) also recalculates these 47 columns with an Advanced Calculation Method and obtains a mean of ϭ 0.966 and standard deviation of ϭ 0.207.Hence, the results from Cyllok and Achenbach (2011) are comparable and the proposed method can be regarded as safe and accurate.
The proposed extension of the Zone Method by Cyllok and Achenbach ( 2009) is mentioned as NCI in the German National Annex (DIN EN 1992-1-2/NA, 2010) to  EN 1992EN -1-2 (2004b) ) and can be used for the design of columns subjected to fire.But the reduction of the compression strength is chosen empirically, which has been criticized by Gelien (2011).Also, the exact limits of this approximation compared to the Advanced Calculation Method have not been determined yet.Furthermore, the Zone Method itself makes many assumptions, which have not been proved with the material laws given in EN 1992-1-2 (2004b).

Proof of the assumptions
The basic ideas of the Zone Method are that thermal strains can be disregarded and that each fiber of the cross-section can achieve its ultimate strength.This means, it is a plastic concept.It cannot be proven analytically -it can be validated by comparing the results with the Advanced Calculation Method.Verification is only possible by the recalculation of laboratory tests.
Therefore, the stress distribution in a wall, considering different material models, is examined.The material models from EN 1992-1-2 (2004b), Anderberg and Thelandersson (1976) and Gernay and Franssen (2012) are used for the Advanced Calculation Method.In the Eurocode model (EN 1992(EN -1-2 (2004b))), the effects of transient thermal strains are included in the stress-strain relationship for concrete (implicit creep).Anderberg and Thelandersson (1976) use an empirical creep function, which has been developed from their own laboratory tests.Gernay and Franssen (2012) also use an empirical creep function.It is based on the material parameters given in ENV 1992-1-2 (1995).Both latter material models use an explicit creep function, which is able to represent stress and temperature dependent stress redistributions.
The stress distribution is displayed in Figure 3 for an axial force per unit length of Ϫ 0.35 • N Rd and Ϫ 0.70 • N Rd with: where A c,net ϭ net concrete cross-section; ␣ ϭ 0.85; f ck ϭ concrete strength at 20°C; ␥ c ϭ 1.5;A s,tot ϭ total area of reinforcement; f yk ϭ steel yield strength at 20°C; and ␥ s ϭ 1.15.The effect of the stress redistribution is evident for the explicit creep models.As expected, it is more distinctive for higher compression forces.Hence, the assumption of Hertz seems to be reasonable, though it cannot be proven in general.It must also be mentioned that the reinforcing steel hinders the transient thermal strains of the concrete fibers, which has not been fully considered by Hertz.
Hertz states that the axial resistance of a cross-section can be calculated with a reduced concrete cross-section [equation ( 4)] and the 0.2 per cent proof strength of the reinforcement [equation ( 7)].The axial resistance of walls using the method of Hertz and the Advanced Calculation Method with the material models from Eurocode, Anderberg and Thelandersson, Gernay and Franssen is calculated for the parameters given in Table I.The net concrete cross-section has been considered in all calculations.The height of the "damaged" zone for the Zone Method has been calculated with equation ( 14) because the reduction of stiffness is decisive for compressed walls.
The results are displayed in Figure 4.For each examined cross-section, the result for the Advanced Calculation Method using material models from Eurocode is plotted on the y-axis.The result of the other models is plotted on the x-axis.All results are above the bisecting line, i.e. the Eurocode model gives the highest axial resistance.The results from the Zone Method are close to the results of the Advanced Calculation Method with the Eurocode material model.Therefore, the calculation of the axial resistance with the Zone Method can be considered as validated for walls.
The peak strains c1, corresponding to f c, are plotted in Figure 5 for the considered material models.The approximation Խ c1,Խ Ϸ 3.5‰/k c () matches the peak strains of the explicit creep models up to a temperature Ϸ 700°C.It is obvious that the peak strains of the implicit creep model of the Eurocode are much larger.
The proposed simplification given in equation ( 8) can be validated by comparison with the explicit creep models.The ratio k Ec () ϭ E c,fi ()/E c,fi (20°C) is pictured in Figure 6 for the considered material models.The proposed approach k Ec () Ϸ (k c ()) 2 is closer to the results from Gernay and Franssen, where the model from Anderson and Thelandersson The calculation of the "damaged" zone a z,EI should also lead to comparable results for each of the considered explicit creep models.The results for a z,f [equation ( 6)], a z,EI [equation ( 12)] and the related simplification given in equation ( 14) is calculated for walls with a height between 10 and 60 cm.The results are displayed in Figure 7.The evaluation of equations ( 12) and ( 14) leads to almost identical results.For h Ն 30 cm the calculated a z,EI is nearly constant.Equations ( 8)-( 12) are also evaluated for the modulus of elasticity for the models from Anderberg and Thelandersson and Gernay and Franssen: the simplification given in equation ( 8) is replaced by E c, of the corresponding model.Both models lead to larger "damaged" zones, the results seem not to be constant for "thick" cross-sections.
The derivation of the damaged zone a z,EI is only valid under the assumptions that the isotherms are parallel to the surface, and all fibers are under constant compression.These assumptions are violated for columns heated on four sides because the isotherms are not parallel.They are even more violated for columns with a large eccentricity because the stress distribution is not constant for the cross-section.For columns, proof can only be given by recalculation of laboratory tests and by the comparison with the Advanced Calculation Method of calculated moment-curvature diagrams.It must be mentioned that the tested columns (Haß, 1986) cover only cross-sections 20 ϫ 20 cm (14 tests), 30 ϫ 30 cm (32 tests) and 30 ϫ 40 cm (1 test).Interpreting Figure 7 within the mainly tested range 20 cm Յ h Յ 30 cm confirms that the simplified calculation of a z,EI given in equation ( 14) is of sufficient accuracy for a simplified calculation method.

Proposal for an extension
As shown in Section 4, the axial resistance of a wall calculated with the Zone Method is close to results calculated by the Advanced Calculation Method, if the resistance of the reinforcement is calculated using the 0.2 per cent proof strength.If the Zone Method for manual calculation is extended toward a general method, suitable for a calculation by computer, the ultimate strength of the concrete f c, and reinforcing steel f s, ϭ k s,0.2% • f yk are replaced by the stress curves given by the Advanced Calculation Method.The new proposal is to reduce the strength of the reinforcement by s .The resistance of the reinforcement is calculated by: where A s ϭ area of reinforcement; s () ϭ reduction factor for reinforcing steel; and s (,) ϭ stress of reinforcing steel.
The parameter s can be interpreted as a reduction of the area of reinforcement, analogous to the factors f and EI for the concrete cross-section.This concept was also proposed by Holmberg and Anderberg (1993) for the 500°C isotherm method.They proposed to reduce the area by 50 per cent, if only four corner bars are present.For more reinforced columns, the corner bars are neglected.
Regarding the results for the axial strength given in Section 4, the limitation to the 0.2 per cent proof strength gives conservative results compared to the Advanced

91
Extension of the Zone Method of Eurocode 2 Calculation Method.The results for the ratio f s,0.2% ()/f s,y (), where f s,y () ϭ yield strength of reinforcement, are displayed in Figure 8 for cold worked steel.
Cyllok and Achenbach proposed to reduce the compression strength of the reinforcing steel by ⌬ ϭ 2% • E s ( s ).This approach has been proofed by the recalculation of laboratory fire test, as explained in Section 3. The reduced compression strength s,red and the yield strength f sy,red are given by: Խ s,red (, ) Խ ϭ Խ s (, ) Խ Ϫ0.002 • E s (), Յ 0 (21) where s (,) ϭ stress of reinforcing steel; E s () ϭ tangent modulus of elasticity of reinforcing steel; and f sy () ϭ reinforcement yield strength at temperature .
The ratio Խ f sy,red () Խ / Խ f sy () Խ is also plotted in Figure 8.It is obvious that the results are conservative for temperatures colder than 400°C: at room temperature only 20 per cent of the yield strength is considered.For temperatures above 400°C, the relative strength is almost constant with a value of 50 per cent of the yield strength at elevated temperatures.
It is proposed to use the following reduction factor s () for the reinforcing steel: where ϭ temperature of the reinforcement (°C).
The reduction factor s () is also displayed in Figure 8.The effect of hindered thermal expansions and corresponding stress redistributions does not occur for reinforcing steel in tension.Therefore, the reinforcing steel is considered to act with its full strength, as written in equation (23a).
• The concrete is represented with a constant temperature M using the stress-strain curves of the Advanced Calculation Method.
• The peak strain of the concrete Խ c1,Խ is at least 3.5‰.
• The stress-strain curves of the Advanced Calculation Method are used for the reinforcement.• The strength of the compressed reinforcement is reduced by s (), as given in equation ( 23b).
The proposed modification for the consideration of the effect of hindered thermal strains does only concern Step 6. Steps 1 to 5 are equal to the scheme proposed by Cyllok and Achenbach, as described in Section 3. The proposed function s () for reduction of the compressed reinforcement is close to the approach by Cyllok and Achenbach, as displayed in Figure 8.Therefore, it can be supposed that the calibration will lead to comparable statistical key data, as given in Section 3.But reliable data can only be obtained by the recalculation of laboratory tests.

Worked example
The column displayed in Figure 9 has been tested at TU Braunschweig (Haß, 1987) as number 46 of the test series.The diameter of the stirrup is not given by Haß (1987).It is assumed to be 8 mm, in accordance with previous test series (Haß and Klingsch, 1980).The yield strength f yk of the reinforcement is considered with 526 MPa.The concrete strength ␤ W,t has been measured with 200 mm cubes at the age of test.The equivalent strength f ck is calculated by equation ( 24) with the factors proposed by Schnell and Loch (2009): where k 150 ϭ strength of 150/200 mm cubes; k cyl ϭ strength of cylinders/cubes; and k cure ϭ strength of wet cured/dry cured concrete.
The moisture is considered with 3 per cent.The thermal conductivity is taken into account with its lower limit (EN 1992(EN -1-2, 2004b)).The temperature distribution of the cross-section is calculated with a finite difference scheme for every time step ⌬t ϭ 1 min.The moment-curvature diagram of the cross-section is also calculated for each time step.For a given curvature, the strain at the centroid is varied until the inner forces N R,fi of the cross-section are equal to the applied axial force N E,fi to satisfy equilibrium.
The transfer matrix method with second-order effects (Petersen, 1982) is used for the determination of the state of strain.The column is divided in ten sections with constant stiffness, which is derived from the moment-curvature diagram.The failure of the column occurs when no equilibrium can be found.
The column failed in the laboratory at t f,exp ϭ 50 min.The calculated time t f,cal for the Advanced Calculation Method is 45 min (Ϫ10 per cent) and 51 min (ϩ2 per cent) for the extended Zone Method.The moment-curvature diagrams for both methods are displayed in Figure 10 for t f ϭ 40 min and t f ϭ 50 min.The curve for the acting moment M E,fi is plotted for reference.It is calculated with a nominal curvature, assuming a parabolic distribution: where ϭ curvature; N E,fi ϭ acting axial force; e 0 ϭ eccentricity; and l col ϭ buckling length of column.
The shape of the moment-curvature curve of the extended Zone Method is comparable to the curve of the Advanced Calculation Method.In this example, the cross-section behavior of the Advanced Calculation Method is "softer" compared to the extended Zone Method.It is also obvious, that the line of the acting moment M E,fi and the curve of cross-section resistance M R,f do not have a clear intersection.This means that the results of calculations close to the ultimate limit state are dependent on the used calculation method and the chosen tolerance levels.

Conclusion
The basic ideas and assumptions of Hertz can be confirmed by the comparison with the results from explicit creep models.The manual calculation scheme can be extended toward a general method, suitable for the implementation in commercial design software.For the "hot design", the material laws of the Advanced Calculation Method can be used, if the strength of the compressed reinforcement is reduced.One "rough" approach has been proposed by Cyllok and Achenbach (2009).This approach has been checked by the recalculation of laboratory tests (Cyllok and Achenbach, 2011).
This first approach is enhanced by introducing the factor s , which can be interpreted as a reduction of the area of reinforcement.In one worked example, taken from the tests carried out at TU Braunschweig, the results are close to the Advanced Calculation Method.
The recalculation of further laboratory tests is needed to calibrate the proposed extended Zone Method.The proof of equivalent structural safety, compared to the Advanced Calculation Method, can only be given by reliability methods.
Figure 1.Concrete wall heated on both sides, cross-section (left), temperature distribution (center) and distribution of concrete strength (right) Figure 2. Square cross-section for the Zone Method by Cyllok and Achenbach (left) and proposal for the extended Zone Method (right)

Table I .
lower modulus of elasticity.As expected, the implicit material model of the Eurocode gives the lowest results. gives