Acta Polytechnica CTU Proceedings doi:10.14311/APP.2017.7.0058 Acta Polytechnica CTU Proceedings 7:58–63, 2017 © Czech Technical University in Prague, 2017 available online at http://ojs.cvut.cz/ojs/index.php/app EXPERIMENTAL, NUMERICAL AND ANALYTICAL INVESTIGATIONS OF WIND-INDUCED NET PRESSURES FOR INDUSTRIAL BUILDINGS WITH ENVELOPE POROSITIES Vanessa Saubke∗, Rüdiger Höffer Ruhr-Universität Bochum, Universitätsstraße 150, 44801 Bochum, Germany ∗ corresponding author: vanessa.saubke@rub.de Abstract. The magnitude and the spatial distribution of wind-induced net pressures (external and internal) on buildings are frequently discussed among research communities and construction industries. This paper deals with this topic based on a case study about an industrial building in Denmark, which was damaged due to the wind impact during a storm when a large part of the roof covering was blown off. In order to detect the reason for the damage the wind-induced loads were studied by i) wind tunnel experiments on the external pressures due to different wind directions, ii) analytical investigations of internal pressure due to envelope porosities and planned openings and iii) numerical analyses for the internal and the external pressure. The Reynolds averaged Navier-Stokes (RANS) method is employed to build a numerical model. The experimental, analytical and numerical results are compared with the indicated characteristic loads from the Eurocode [1]. Keywords: net pressure, internal pressure, industrial building, Eurocode. 1. Introduction The storm "Christian" caused a lot of significant dam- ages in Northern Europe in October 2013, including also failure of supporting structures of buildings. Espe- cially in Denmark high wind velocities were measured. According to the Danish meteorological service local wind peaks of 53.5 m/s were listed, which were the highest values ever measured in Denmark. The wind velocity distribution for the maximum values for this day is given in Figure 1 (b). The diagram in Figure 1 (c) shows the data of the two measurement stations nearby the building which is considered here. The peaks of the wind velocity at 1 p.m. and at 3 p.m. are clearly visible. Figure 1. (a) Cardinal points (Google Earth [2], modified), (b) distribution of the peak wind velocity in Denmark during hurricane "Christian" in October 2013 (Danmark Meteorologiske Institut [3], modified) and (c) measurements of the wind velocities at stations near to the building. At 3.15 p.m. the industrial building was damaged due to the wind impact. A large part of the roof structure breaks off and leads to an extended hole in the roof construction (see Fig. 3 (a)). The debris were blown over the building and damaged also the leeward building (see Fig. 2). The wind direction was west-southwest (for the cardinal points see Fig. 1 (a)). Figure 2. Pictures from a security camera taken during the storm in 2013, the roof failure is visible. The damage started at the edge of the roof where the foil peeled off. This led to a deformation of the parapet wall (see Fig. 3 (d)). Approximately 7.5 % of the roof construction was completely blown off, additional 6 % of the roof foil was torn off and the trapezoidal sheets were folded and buckled (see Fig. 3 (b)). In total more than 10 % of the roof construction was destroyed. Futhermore this caused a damage at the sprinkler system and the leaking water destroyed a big part of the stored goods. The impact of the wind load was huge so that also the girders of the supporting roof construction were deformed (see Fig. 3 (c)). 58 http://dx.doi.org/10.14311/APP.2017.7.0058 http://ojs.cvut.cz/ojs/index.php/app vol. 7/2017 Exp., num. and anal. investigations of wind-induced net pressures Figure 3. Qualitative description of typical damages: (a) Partically opened roof, (b) buckled sheets, (c) bent girder and (d) deformed parapet. The considered building is located in the east part of Denmark. It is situated at the west of a complex of three storage buildings of the same type. Each building is about 135 m long, 55 m wide and 38 m tall. The shelfs for the stored goods are part of the supporting steel construction which consists of a regular framework of stiffening frames (see Fig. 4). This main structure is covered by a sandwich-element facade. Figure 4. Example of a warehouse building under construction (Google Maps [2]). The roof consists of supporting girders, trapezoidal sheets, mineral insulating panels and a roof foil. The sheets are connected to the supporting girders by screws and the foil is fixed by spikes. For the sandwich- element facade a connection type with unsealed joints was used. After Lange and Rädel [4] this specific joint construction can lead to a leakage of the building walls. For the consideration of this leakage an effective joint was estimated as 3.5 mm (for investigations on the leakage of joint see Kuhnhenne [5]). With this information the general porosity φ, the relation of the area of the openings to the area of the surfaces, was calculated: φ = ∑ A′∑ A = 0.0026 = 0.26%. (1) At the southern facade the gates for the delivery of the goods are located behind an antecedent flat hall. If the gates are not closed the total area of openings is 61.43 m2. 2. Materials and Methods 2.1. Wind loads after 1991-1-4 After the Eurocode external and internal pressures due to wind impact are taken into account [1]. A superposition of both leads to the net pressure. The definition of the signs of this values are given in Figure 5. Figure 5. Three typical cases for the internal pressure (Cook [6], modified). The code defines the wind load we,i as multipli- cation of the peak velocity pressure and a pressure coefficient: we = qp · cpe and wi = qp · cpi. (2) The peak velocity pressure qp is influenced by different factors like the terrain category, the building height and the wind zone. The external pressure coefficient cpe is influenced by the shape of a building and the internal pressure coefficient cpi by the distribution of the openings over the four faces of the building. In the following investigations pressure coefficients are considered and compared. For this purpose the external pressure coefficients for the roof of this building are determined in Table 2. Figure 6. Distribution of the pressure coefficients for a flat roof with (a) wind direction on the short side and (b) wind direction on the long side. Regarding the flow around a rectangular cube of a certain height (see Hucho [7]) the code estimation of 59 Vanessa Saubke, Rüdiger Höffer Acta Polytechnica CTU Proceedings the distribution of the pressure coefficients is given in Figure 6. The distribution depends on the dimensions of the building and has to be estimated for every building and every wind direction individually due to the regulations. The cpe,1 values should be used for loaded areas of less than 1 m2 and the cpe,10 values for loaded areas of more than 10 m2 with the possibility of a logarithmic interpolatin for intermediate values. The internal pressure can be calculated according to two different regulations: (1.) For homogeneously distributed openings the pres- sure coefficient can be determined from a diagram by calculating the opening ratio µ, which is the re- lation from the area of the openings in the suction region divided by the area of all openings: µ = ∑ area of openings with cpe ≤ 0∑ area of all openings (3) (2.) A Face with an area of openings that is more than twice the area of the sum of the openings in the remaining faces is considerd as a dominant face. • The internal pressure coefficient can be calculated with the following equation: cpi = 0.75 · cpe. (4) • If the area is more than three times the remaining openings then it is calculated with the following equation: cpi = 0.9 · cpe. (5) The considered building has firstly a general porosity because of the leakage due to the facade-elements. Secondly it has large openings for the delivery of the stored goods at the short face where the open flat building is located. It is not known yet if openings were open or closed during the storm. Therefore, three different situations for two wind directions are consid- ered. These different cases are explained in Table 1. The results for the internal pressure coefficients for these six cases are listed in the Table 3. Case Windward Face Porosity Delivery Gates 1 X 2 short X X 3 X 4 X 5 long X X 6 X Table 1. Description of the six load cases. 2.2. Experimental tests Wind tunnel tests were performed at the Boundary Layer Wind Tunnel at Ruhr-Universität Bochum to build up a reference data-base of wind-induced loads on the structure for specifically chosen spots of the building envelope, especially on the parts damaged during the storm event. The focus of the measure- ments was set on the roof construction which failed (see Fig. 7 (a)). The effect of the parapet wall was ne- glected due to the small size in relation to the building. This leads to results on the safe side. Figure 7. (a) Positions of the pressure measurement points and (b) the principal set-up of the model. The tests were performed in a geometrical scale of 1:300 for different wind directions from southwest to northwest. The building was investigated as a (i) single standing, isolated building and (ii) in a group arrangement with the neighboring buildings (see Fig. 7 (b) for the set-up). The evaluated cpe values repre- sent the 78 % fractiles of the analyzed measurement data. The pressures are measured over a point-like averaging area which represents an area smaller than 1 m2 in full scale. The minimum values from the flow direction on the long face and of all measured directions for every zone are listed in Table 2. Figure 8. (a) The minimal cpe values of the flow direction west (on long face) and (b) of all measured flow directions. 2.3. Numerical simulations For the numerical investigations a CFD method ap- plying finite volumes and the RANS method with the standard k-�-model is used. The applied inlet velocity profile was introduced by Richards and Hoxey [8] in 1993 and is commonly used for this method. Effects of building interferences were neglected by running all simulations with a single building. Figure 9. (a) Idealization of the joints and (b) ideal- ization of the openings. 60 vol. 7/2017 Exp., num. and anal. investigations of wind-induced net pressures Because a simulation of all explicit openings and joints would lead to an extremely high number of volume elements the openings were idealized. The general porosity of the joints was considered with four explicit idealized joints distributed over the height (see Fig. 9, left). The openings of the delivery gates were also idealized as one large opening with the same area which is located in the center of gravity of all gates (see Fig. 9, right). Figure 10. Numerical results for the internal pressure coefficients cpi: (a) case 1, (b) case 2, (c) case 3, (d) case 4, (e) case 5 and (f) case 6. Figure 11. Numerical results of the external pressure coefficients of the roof for (a) flow on the short side and (b) flow on the long side. Geometrical models and meshes were generated for all six cases (see Table 1) with approximately three million volume elements per case. The distribution of the external pressure coefficients is pictured in Figure 11. The maximum values for every zone are listed in Table 2. The internal pressure coefficients for the six different cases are plotted in Figure 10 and numbered in Table 3. 2.4. Analytical calculations An analytical investigation is performed for the es- timation of the internal pressure. As Cook [6] 1990 already described, the internal pressure pi depends on the external pressure pe at the openings A′. It is defined as the balance between the inflow and outflow of the building (see also Fig. 5 (c)): inflow∑ ( A′ (pe − pi) 1 2 ) = outflow∑ ( A′ (pi − pe) 1 2 ) . (6) The equation cannot be solved directly but iteratively with an estimated value for the internal pressure. The left-hand and right-hand side of Equation 6 are not fixed during the iteration process, because an outflow can become an inflow and vice versa. Figure 12. Distribution of the external pressure coefficients from a wind tunnel test for (a) the wind direction on the short face and (b) the wind direction on the long face (TPU Aerodynamic Database [9], modified). To calculate the internal pressure analytically a dis- tribution of the external pressure is needed. Because the performed wind tunnel tests are putting the focus on the damaged area of the roof construction, com- plete roof and wall distributions of external pressures have not been measured and are missing for analytical calculations of internal pressures. Therefore mean val- ues of external pressure coefficients from wind tunnel tests of a building with similar proportional dimen- sions from the TPU Aerodynamic Database [9] were used (see Fig. 12). Because the building from the database was smaller than the original one, the area of the openings was scaled down whereas the opening ratio µ stayed unchanged. In this scale the distribu- tions of the pressure coefficients were evaluated and the area of the different pressure coefficients calcu- lated. In order to maintain the comparability the same idealized opening of the delivery gates like those for the numerical investigations in Section 2.3 were used. For the joints a general porosity of 0.26 % was considered. The internal pressure was calculated for all six cases as given in Table 1 after Cook’s method, [6]. The results are listed in Table 3. 61 Vanessa Saubke, Rüdiger Höffer Acta Polytechnica CTU Proceedings 3. Results The results of the wind tunnel investigations corre- spond relatively good to the required values of the Eurocode, for example the Eurocode requires a value of cpe = 2.2 for the zone F and the wind tunnel results leads to the highest value of cpe = 2.14 considering all investigated directions. The results are compared to the cpe,1 values although the area of the measure- ment points in the experiments correspond to an area of smaller 1 m2 in full scale. The normative values of zone F and G suit to the minimum values of the testing results for all directions while the normative value of zone H fits better to the averaged value. But as one can see in Figure 8 the values of this zone have a certain variation, they decrease in direction to the leeward side of the building. Therefore, it is obvious that a design for the maximum value would not be economic, but one should be aware of the possibility that higher values can occur. Zone F G H Eurocode1) -2.2 -1.8 -1.2 Windtunnel test (long face) -1.64 -1.34 -1.43 (-1.192)) Windtunnel test (all directions) -2.14 -1.84 -1.63 (-1.292)) Numerical simulation (Case 1-3) -1.19 -1.11 -1.01 Numercial simulation (Case 4-6) -1.18 -1.13 -1.03 Table 2. Comparison of the external pressure coeffi- cients cpe for the roof construction. 1)cpe,1 values in order to maintain comparability 2)averaged value of zone H The cpe values of the wind tunnel tests for the wind on the long face can be compared to the numer- ical results of Cases 4 to 6. It should also be taken into account that the RANS method does not lead to peak values and deliver average value only. The distribution of the external pressure coefficient shows a satisfactory similarity to the experimentally deter- mined distribution from the windtunnel tests. But RANS underestimates the external pressure in the separation area because the turbulence is modelled mathematically and local instationary vortices are neglected. The values of the numerical results are between 72 % to 84 % of the values of the wind tunnel test. The investigations of the internal pressure coefficients lead to different results for the six cases. The mag- nitude of the values are similar, in most of the cases, like Case 1, 3 and 6, the Eurocode overestimates the values and is on the safe side, but in Case 5 it under- estimates the internal pressure in comparison to the analytical solution. Case 1 2 3 4 5 6 Windward face short long Porosity X X X X Delivery gates X X X X Eurocode -0.30 0.34 0.64 -0.03 -0.33 -0.83 Numerical simulation -0.21 -0.12 0.52 -0.24 -0.31 -0.63 Analytical calculation -0.20 0.34 0.59 -0.23 -0.48 -0.62 Table 3. Comparison of the internal pressure coeffi- cients cpi for the six considered load cases. In Case 2 the numerical analysis lead to a complete different value than the results of the analytical solu- tion, which seems to fit perfectly to the value of the Eurocode. In a further investigation, e.g. with the LES-method, the numerical results should be checked, but it is sure, that this is connected to the position of the opening. If the position of the opening in the analytical calculation is changed to a higher position in the middle of the face (see Fig. 13), the value increases from cpi = 34 to cpi = 0.45. Figure 13. Changed position of the opening at the windward face (TPU Aerodynamic Database [9], mod- ified). Also in some additional numerical investigations on a 1m x 1m x 1m cube, where the position of the openings over the height was changed, the same behaviour was determined (see Fig. 14). This shows that the internal pressure reacts quite sensitive to the position of the openings, especially for openings in the windward side. Figure 14. Investigation on a 1m x 1m x 1m cube due to the position of the opening: (a) Opening at the bottom, (b) opening in the middle and (c) opening at the top. In Case 4 the Eurocode diagram leads to a value of nearly zero (with µ = 0.63 in Fig. 15) in contrast to the numerical and analytical solution, which leads to good corresponding results of -0.24 and -0.23. In 62 vol. 7/2017 Exp., num. and anal. investigations of wind-induced net pressures the old German code near to zero values were avoided. In the range of the diagram where the value leads to zero one has to consider two values (0.2 and -0.3) on the safe side (see the red square in Fig. 15 (b)). Figure 15. Diagram for the cpi values: (a) Eurocode [1] and (b) the previous code DIN 1055 [10]. 4. Conclusions The results show, that the regulation of the internal pressure for a dominant face is conservative for di- mensioning a building (see regulation (2.) in Section 2.1). The regulation for homogeneously distributed openings however seems to falsify the internal pressure (see regulation (1.) in Section 2.1). As well for very homogeneously distributed openings, like a general porosity, the Regulation (1.) does not work properly, like in Case 4. Also in Case 5 the value is underesti- mated in comparison to the analytical value. The normative DIN EN 1991-1-4 [1] offers the oppor- tunity to use the fixed values of +0.2 and -0.3, if it is not possible or not justifiable to estimate a µ. These values seems to represent an underestimation of the effects of the internal pressures as one can understand from the results of the investigations for cpi ranging from +0.59 to -0.63. This apparent load underesti- mation through the quoted simplified code regulation does not necessarily lead to damages of the structure as it can be assumed to be based on experiences and sound practice requirements which are included into code regulations. However, such configuration can be insecure in connection with a certain systematical underestimation of external pressures after specific design regulation for trapezoidal sheets. For example, the German DIN 18807-3 [11] considers the higher pressures at the edge of the roof (zone F and G) just for the design of the screw connections; for the de- sign of the sheets the high suction in these zones is neglected, which leads to an underestimation of the external pressure at the edge of the roof. Together with the simplified Eurocode regulation as described above this can result into a safety relevant lack of net design pressures. The research on the topic of the internal pressure should be continued to validate the results. To per- form some additional wind tunnel tests on buildings with different proportions in combination with numer- ical simulations would be useful. For the numerical simulations some further methods should be used to improve the results. The aim would be to find a sta- ble model that leads to good results for the different proportions so that it can be used as a tool to check and improve the regulation of the code. List of symbols A Area of the faces [m2] A′ Area of the openings [m2] cpe External pressure coefficient cpi Internal pressure coefficient pe External pressure [N/mm2] pi Internal pressure [N/mm2] we External pressure (Eurocode) [N/mm2] wi Internal pressure (Eurocode) [N/mm2] µ Opening ratio φ General porosity References [1] E. C. for Standardization. EN 1991-1-4 - Eurocode 1: Actions on structures - Part 1-4: General actions - Wind actions. 2010. [2] Google Maps. Panoramio, 2007. [3] Danmark Meteorologiske Institut (DMI). Orkan satte Danmarksrekord i vind - se kortene, 2013. [4] J. Lange, F. Rädel. Die Fugendichtigkeit von Sandwichelementkonstruktionen - Wasser- und Luftdichtigkeit in Längsfugen und Fensteranschlussfugen. In Festschrift zum 60. Geburtstag von Gerhard Hanswille, Heft 20. Institut für Konstruktiven Ingenieurbau, Bergische Universität Wuppertal, 2011. [5] M. Kuhnhenne. Energetische Qualität von Gebäudehüllen in Stahl-Sandwichbauweise. Dissertation, Rheinisch-Westfälische Technische Hochschule Aachen, Fakultät für Bauingenieurwesen, 2009. [6] N. J. Cook. The designer’s guide to wind loading of building structures - Part 2: Static structures. Butterworths, 1990. [7] W.-H. Hucho. Aerodynamik der stumpfen Körper. Vieweg+Teubner, 2nd edn., 2011. [8] P. Richards, R. Hoxey. Appropriate boundary conditions for computational wind engineering models using the k-� turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46-47:145–153, 1993. [9] T. P. U. Aerodynamic Database. Database of isolated low-rise building without eaves. http://www.wind.arch.t-kougei.ac.jp/info_ center/windpressure/lowrise/mainpage.html. [10] DIN 1055-4 - Action on structures - Part 4: Wind loads. DIN Deutsches Institut für Normung e.V., 2005. [11] DIN 18807-3 - Trapezoidal sheeting in buildings; steel trapezoidal sheeting; strength analysis, structural design. DIN Deutsches Institut für Normung e.V., 1987. 63 http://www.wind.arch.t-kougei.ac.jp/info_center/windpressure/lowrise/mainpage.html http://www.wind.arch.t-kougei.ac.jp/info_center/windpressure/lowrise/mainpage.html Acta Polytechnica CTU Proceedings 7:58–63, 2017 1 Introduction 2 Materials and Methods 2.1 Wind loads after 1991-1-4 2.2 Experimental tests 2.3 Numerical simulations 2.4 Analytical calculations 3 Results 4 Conclusions List of symbols References