Author: Zbigniew Respondek
Source: Materials 2020, 13(2), 286
https://doi.org/10.3390/ma13020286
(This article belongs to the Special Issue Energy in Construction and Building Materials)
Abstract
One of the structural elements used in the construction of insulating glass units (IGUs) are tight gaps filled with gas, the purpose of which is to improve the thermal properties of glazing in buildings. Natural changes in weather parameters: atmospheric pressure, temperature, and wind influence the gas pressure changes in the gaps and, consequently, the resultant loads and deflections of the component glass panes of a unit. In low temperature conditions and when the atmospheric pressure increases, the component glass panes may have a concave form of deflection, so that the thickness of the gaps in such loaded glazing may be less than its nominal thickness.
The paper analyses the effect of reducing this thickness in winter conditions on the design heat loss through insulating glass units. For this purpose, deflections of glass in sample units were determined and on this basis the thickness of the gaps under operating conditions was estimated. Next, the thermal transmittance and density of heatflow rate determined for gaps of nominal thickness and of thickness reduced under load were compared. It was shown that taking into account the influence of climatic loads may, under certain conditions, result in an increase in the calculated heat loss through IGUs. This happens when the gaps do not transfer heat by convection, i.e., in a linear range of changes in thermal transmittance.
For example, for currently manufactured tripleglazed IGUs in conditions of “mild winter”, the calculated heat losses can increase to 5%, and for doubleglazed IGUs with 10–14 mm gaps this ratio is about 4.6%. In other cases—e.g., large thickness of the gaps in a unit, large reduction in outside temperature—convention appears in the gaps. Then reducing the thickness of the gaps does not worsen the thermal insulation of the glazing. This effect should be taken into account when designing IGUs. It was also found that the wind load does not significantly affect the thickness of the gaps.
1. Introduction
Generally used in the building industry as a filling of windows or glass facades, insulating glass units (IGUs) consist of two or more component glass panes, connected at the edges with a glass spacer. The space between the component glass panes forms a tight gap filled with gas. In order to improve the thermal performance of the building partition constructed this way, the gap is filled with gas with lower thermal conductivity than air, most often argon. Further improvement of the performance is achieved by the use of component glass panes with a lowE coating—such a coating must be located on the side of the gap because it corrodes quickly when exposed to weather conditions. The tightness of the gap in insulating glass units is therefore a necessary factor to maintain good thermal insulation of transparent glazing [1,2].
Tight gaps, however, determine some specific properties of IGUs in the context of environmental loads transfer and the associated deformation of structural elements. The gap is filled with gas in the production process of the unit, therefore the gas in the gap has some initial parameters of pressure, temperature, and volume. Under operating conditions, an insulating glass unit is exposed to climatic loads which generate loads and deflections of the component glass panes due to the pressure difference between the gap and the environment. For example, an increase in atmospheric pressure or a decrease in the gas temperature in the gap results in a concave form of deflection of the panes (Figure 1a) and the opposite changes of these parameters in a convex form (Figure 1b). The magnitude of the under or overpressure in the gap depends not only on the value of climatic loads, but also on the structure of IGU. In general, it increases with reduced IGU dimensions (width × length), increased thickness of the gas gap, and increased thickness of the component glass plates. How the pressure difference affects the deflections in IGU will be presented later in this article.
In the case of wind exposure (Figure 1c), the tightness of the gap has a positive effect on the load distribution in an IGU. Due to changes in gas pressure in the gaps, the external load is partly transferred to the other panes of the unit.
The deflections of the glass described above result in deformation of the image viewed in the light reflected from the glass in windows or on glass facades (Figure 2). It is important that under conditions of low air temperatures, i.e., during the heating season, insulating glass units tend to take the concave form of deflection. The result is a reduction in the thickness of the gas space—especially in the central part of the glazing, where the component glass panes are closest to each other—which makes it possible to reduce the thermal insulation of the IGU.
The aim of the analysis carried out in the paper was to determine the effect of taking into account the reduction in thickness of gasfilled gaps in insulating glass units in winter conditions on the calculated heat losses through these partitions. The analysis was made for example for double and tripleglazed IGUs. A detailed numerical quantification of this phenomenon was carried out for various IGU constructions.
In the literature, studies describing previous research in this area can be found: Barnier and Bourret [3] analyzed the effect of plate curvature in IGUs on the thermal transmittance (Ugvalue). The authors determined the Ugvalue for IGUs with variable gap thickness (limited by the surfaces of deflected panes), considering the average gap thickness in the loaded IGU as reliable. The authors stated that this assumption becomes reasonable when plate curvature is small and it is certainly acceptable in the conduction regime, where the convective movement is not significant. This article provides the results of sample calculations for double and tripleglazed units under winter conditions. It was found that taking into account the plate curvature increases the calculated Ugvalue from 4.4% to 5.8%. Calculations were also made to account for changes in weather conditions (typical meteorological year) for Montreal and Toulouse. The results indicate that Ug may vary up to 5% above and 10% below the yearly average.
Hart et al. [4] analyzed the Ugvalue calculated from real deflections of double and tripleglazed units, measured in summer and winter at several locations in the USA. It was found that a 20 °C temperature difference reduces thermal performance by 4.6% for doubleglazed IGUs and by 3.6% for tripleglazed IGUs.
Penkova et al. [5] presented examples of numerical analysis and experimental research regarding both parameters related to heat flow and climate loads. However, no detailed analysis of the change in thermal transmittance related to the deflections in the IGUs was carried out.
Thermal imaging photographs illustrating a decrease in thermal insulation in the central part of the glazing were published as well [6]. Examples are presented where the temperature in the central part of the glazing is 1–3 °C higher than the average on its surface (in images from the outside). An example of an IGU in which the component panes came into contact due to climatic loads is also presented.
2. Methodology for the Calculation of Static Quantities in IGUs
The methods of calculation of static quantities in doubleglazed IGUs loaded with climatic factors are described in the literature. Mention may be made here of analytical models presented in papers [7,8,9,10] and numerical models allowing for consideration of the possibility of nonlinear deflections of component glass [11,12]. The results of calculations presented in this paper were obtained using the author’s analytical model proposed in the article [13], which allows to calculate the load and deflection of component glass panes in units with any number of tight gaps.
The basis for the calculation of static quantities in IGUs is the assumption that the gas in the gaps meets the ideal gas equation:
where:
 p₀, T₀, v₀—initial gap gas parameters: pressure [kPa], temperature [K], volume [m³], obtained in the production process,
 pop, Top, vop—operating parameters—analogously.
It is also assumed that the glass panes are simply supported at the edges and that the linear dependence of deflection w [m] of the component glass pane on its resultant surface load q [kN/m²] is assumed. The latter assumption is considered to be sufficiently accurate if the deflection is not greater than the thickness of the glass [14]. The deflection function of a simply supported single pane of the a [m] width and b [m] length, subjected to the q [kN/m²] load, placed centrally in the xy coordinate system, can be recorded as [15]:
with
D [kNm] is the flexural rigidity of glass pane:
where:
 d—is the glass pane thickness [m],
 E—is the Young’s modulus of glass [kPa],
 μ—is the Poisson’s ratio [].
Change in gap volume ∆v [m³] resulting from the deflection of one of the limiting glass pane may be determined by integration of Equation (2):
After the relevant calculations have been made:
where:
 α′v—is the dimensionless coefficient dependent on the b/a ratio (Table 1) [],
 αv—is the proportionality factor, [m⁵/kN].
Any change in climatic conditions (atmospheric pressure, temperature, wind) results in a change in the gas pressure in the gaps which affects the resultant operating load of each of the component glass panes. For each gap of an IGU it is possible to formulate the equation of state:
where:
Σ∆v—is the change in gap volume caused by deflection of both panes limiting it [m³].
As already mentioned in the article, double and tripleglazed IGUs were analyzed. In the remaining part, the parameters of the individual component glass panes and gaps were marked with appropriate indices (Figure 3). It is also assumed that loads and deflections are positive if they face the interior, i.e., from left to right as in Figure 3.
Taking into account the adopted markings and conventions, Equation (8) for a doubleglazed IGU can be presented in the form:
with
where:
 pa—current atmospheric pressure [kPa],
 qz,ex, qz,in—load per area from outer factors, primarily wind [kN/m²], almost always qz,in = 0.
After the relevant transitions have been made:
with
Equation (11) has one solution giving nonnegative results:
In the case of a tripleglazed unit, a system of quadratic equations should be solved:
This system has no analytical solution, but it can be solved numerically by iteration.
After calculating the operating pressure pop for each of the gaps, the resultant loading q for each of the component glass panes can be determined:

for a doubleglazed IGU
 for a tripleglazed IGU
Deflection wc [mm] in the center of the glass pane can be determined by the formula:
where:
α′w—is the dimensionless coefficient dependent on the b/a ratio (Table 1) [].
However, the average deflection of the component glass panes wₘ [mm] was determined from the formula:
Table 1. Coefficients for calculating volume change and deflection for simply supported glass pane.
3. Materials and Methods
Thermal transmittance Ug [W/(m²·K)] of IGUs was calculated on the basis of the methodology described in standard [16], and heat losses were expressed by density of heatflow rate Φ [W/m²] from the formula:
where:
tᵢ, tₑ—are the internal and external air temperature [°C].
The heat flow through an insulating glass unit is complex—through conduction, convection, and radiation. The thermal resistance of gasfilled gaps Rs [(m²·K)/W] has the greatest influence on the Uvalue. For each gap:
with
where:
 hr—is the thermal conductance by radiation [W/(m²·K)],
 hg—is the thermal conductance of gas (by conduction and convection) [W/(m²·K)].
 λg—is the thermal conductivity of gas [W/(m·K)],
 s—is the gas gap thickness [m],
 Nu—is the Nusselt number [],
 σ—is the Boltzmann constant 5.6693 × 10⁻⁸ W/(m²·K⁴)
 Tm—is the average temperature of both surfaces delimiting the gap [K],
 εsur1, εsur2—are the emissivity of surfaces delimiting the gap [].
It is particularly important whether convection occurs in the gaps. In the case of narrow gaps (Nu < 1) it is assumed that convection does not occur—thermal insulation of the gap increases linearly with its thickness. If a certain gap thickness limit (for Nu = 1) is exceeded, the effect of convection is taken into account. In this nonlinear range (for Nu > 1) thermal insulation of the IGU does not improve. The value of this thickness limit depends on many factors (see also [17,18]), first of all on:
 the type of gas; the calculation was based on the use of argon,
 location in the structure; the calculations assume a horizontal position, in units situated horizontally or diagonally convection increases.
 increasing the temperature difference on the surfaces of the glass panes limiting the gap affects the increase in convection,
 convection also increases when the average gas temperature in the gap increases.
The thermal resistance of the gaps is primarily influenced by the use of lowemission glass. Glass without coating has a standard coefficient of emission of ε = 0.837. Application of lowemission coating reduces the emissivity of the plate surface, which results in a significant reduction of heat transfer by radiation. Currently, in Central and Northern Europe, IGUs are most often produced, in which each gap is adjacent to one coated surface and one without coating (Figure 3). This solution is most often used in units currently produced in Central and Northern Europe. The values εsur1 = 0.837 and εsur2 = 0.04 were used in the calculations.
Glass conducts heat well, therefore the thickness of the glass panes has no significant effect on the Ugvalue. Physical parameters of argon were adopted on the basis of the standard [17].
Of course, thermal insulation is also affected by thermal surface resistance at the external side (Re [(m²·K)/W]) and at the internal side of a window (Ri [(m²·K)/W]). They depend primarily on the positioning of the window in the structure and the velocity of air (a short analysis on this subject is presented in Chapter 5). The calculations assumed Ri = 0.13 (m²·K)/W (vertical position) and Re = 0.04 (m²·K)/W (for wind velocity V = 4 m/s). These are often accepted comparative values, also in the standard [16].
The calculations according to the adopted model require the use of numerical methods, because we encounter several interdependent values here. For example, the temperature values of gas and glass surfaces depend on the temperature distribution in the crosssection of the IGU. This distribution depends on the resulting thermal resistance values. The results of calculations were obtained by iteration after building the appropriate spreadsheet, assuming the steady state of heat transfer.
Figure 4 shows the effect of gap thickness s [mm] on the design Ugvalue for double and tripleglazed IGUs, with glass thickness d = 4 mm, assuming ti = 20 °C and in two variants of the outside air temperature te = 0 °C and te = −20 °C. The dashed line was used to determine the limits of gap thickness at which Nu = 1.
Figure 4 shows that at low temperatures te the limit thickness decreases. It can also be stated that in the case of tripleglazed IGUs, the difference in temperature in the gap is smaller, and the thickness of the boundary increases.
4. IGUs Under Pressure and Temperature Changes—Presentation and Discussion of Test Results
An analysis of the influence of climate loads on heat loss through IGUs under winter conditions was carried out for sample units with dimensions 0.7 × 1.4 m. Glass material parameters were adopted according to the standard [19]: E = 70 GPa, μ = 0.2.
It was also assumed that the following initial parameters were obtained in the argonfilled gaps during the production process T₀ = 20 °C = 293.15 K, p₀ = 100 kPa. In these conditions, the component glass panes are flat.
Two variants of the temperature drop load were used.
Variant 1. Reduced temperature conditions: ti = 20 °C, te = −20 °C; the gas temperature in each gap was calculated for each case based on the temperature distribution in the particular IGU: for a doubleglazed IGU Top1 = −2.37 to −2.25 °C, for tripleglazed IGU Top1 = −10.09 to −9.66 °C, Top2 = 7.60 to 7.97 °C.
Variant 2. Conditions for a “mild winter”: ti = 20 °C, te = 0 °C; gas temperature: for a doubleglazed IGU Top1 = 8.80 to 9.01 °C, for tripleglazed IGU Top1 = 4.97 to 5.08 °C, Top2 = 13.66 to 14.10 °C.
First, the effect of varying glass thickness on the gap width in the loaded set was investigated. IGUs with 16 and 12 mm gap thickness were analyzed in various combinations of 3, 4, and 6 mm thick panes. It was assumed that IGUs are only loaded with the temperature drop, as in variant 1, i.e., the current atmospheric pressure pa = p₀ = 100 kPa. The results of the calculations are presented in Table 2.
Table 2. Static quantities and gap thicknesses in IGUs—under reduced temperature conditions (Variant 1).
The resultant loading qex and qin (absolute value) illustrates the underpressure in the gaps in relation to atmospheric pressure. The parameter q₁₋₂ is the difference in operating pressure between the gaps in a tripleglazed IGU. From Equations (18) and (19) the extreme deflection (in the center of the pane) wc and the average deflection wm (the wₘ values are given between parentheses) were calculated for each pane.
On the basis of these deflections, the minimum gap thickness in the center of the IGU sc [mm] and the average gap thickness sₘ[mm] were calculated.
On the basis of the calculations presented in Table 2, it was found that the calculated values of sc and sₘ for analyzed IGUs with gaps of the same nominal thickness do not differ much from each other. This is despite the fact that the deflection of component glass panes varies considerably. The effect of gas interactions in tight gaps can be seen here. Rigid panes are less susceptible to deflection, but the external load is less compensated for by the gas pressure in the gap. After changing the thickness of all the panes in a unit, the load changes, but the deflections are similar. Therefore, when one of the glass panes changes to a stiffer one, the absolute load values of component glass panes increase, although their algebraic sum for each IGU is equal to 0. After such conversion, the less rigid panes deflect more because they are exposed to a higher loading—for this reason, a loaded IGU has approximately constant volume of gaps, despite the change in thickness of the component panes.
To identify the extent of the phenomenon described above, another example was solved. Figure 5 shows the influence of IGU width (at a constant ratio b/a = 2) on the maximum deflection of component panes wc in doubleglazed units at 3, 4 i 6 mm thick panes and 16 mm thick gap. It was assumed that IGU is loaded only by a change in atmospheric pressure by ∆p = pa − po = 3 kPa. This means that the current atmospheric pressure is pa = 103 kPa. It can be added here that the results of calculations of static quantities are not very sensitive to the value of po, and to a significant extent to ∆p. This means that if we assumed, for example, po = 950 kPa i pa = 980 kPa, the results would be almost identical.
Figure 5 shows that greater diversity of deflections for units of different component panes thickness occurs in the case of smaller IGU sizes. Then, however, the deflection values are smaller and it can be expected that changes in the gap thickness are also small.
In the context of the above, further analysis was carried out for IGUs with the same thickness of component glass panes d = 4 mm and different gap thicknesses s were assumed. Table 3 presents calculations of wc, wm, sc, and sm values for units loaded with temperature change as in variant 1 and simultaneously operating loading with external atmospheric pressure increase of ∆p = 3 kPa. These are particularly unfavorable operating conditions in the context of reducing the thickness of the gaps. An analogous calculation was carried out for variant 2 (Table 4).
Table 3. Static quantities and gap thicknesses in IGUs under reduced temperature conditions (Variant 1) and atmospheric pressure increase by ∆p = 3 kPa.
Table 4. Static quantities and gap thicknesses in IGUs under conditions for a “mild winter” (Variant 2) and atmospheric pressure increase by ∆p = 3 kPa.
Based on the above data, Table 5 (for variant 1) and Table 6 (for variant 2) present the results of calculations of the thermal transmittance U and the density of heatflow rate Φ [W/m²]:
Table 5. Quantities describing heat losses by IGUs under conditions of reduced temperature (Variant 1) and atmospheric pressure increase by ∆p = 3 kPa.
Table 6. Quantities describing heat losses by IGUs under conditions for a “mild winter” (Variant 2) and atmospheric pressure increase by ∆p = 3 kPa.
 Ug, Φg—describe heat loss without taking into account the curvature of the panes, calculated for the nominal thickness of the gaps,
 Uc, Φc—describe possible local heat loss near the IGU center, i.e., where the distance between the panes is the smallest, calculated for the thickness of the gaps sc,
 Um, Φₘ—describe the average heat loss through the IGU, calculated for the thickness of the gaps sₘ.
Finally, the percentage increase in the calculated quantities is also presented (∆Φc, ∆Φm) for units of nominal gap thickness.
The data presented in Table 5 and Table 6 and Figure 4 indicate that the reduction in the thickness of the gaps of insulating glass units due to their deflection under a drop in gas temperature and a rise in atmospheric pressure may result in an increase in design heat losses in relation to the calculations without taking into account the curvature of the panes. The increase in heat loss occurs in the linear range of the Ugvalue change, i.e., when the conditions inside the gap lead to Nu < 1. It is different when the Ugvalue changes in the nonlinear range (Nu > 1). Heat losses do not increase. Then, the reduction of gap thickness can lead to a slight decrease in the calculated Nu value, which translates into a slight reduction in the calculated heat losses.
In this context, it is preferable to design IGUs such that it has Nu > 1 with a certain margin based on glazing deflections. However, this task should be approached with great caution, taking into account local climate conditions. It is necessary to check if the thickening of the gaps between the panes will not lead to excessive overpressure during the summer, due to the heating of gas in the gaps.
One more feature of the described phenomenon should be noted. In the linear range of changes in the Ugvalue, the indices ∆Φc and ∆Φm almost do not depend on the thickness of the gaps. This is due to the fact, as additional calculations have shown, that the relationship between the thickness of IGU gaps and static quantities (resultant loading of component glass and their deflections) is also linear.
For many years, doubleglazed IGUs dominated the market. Currently, due to the need to save energy, in Central and Northern Europe, tripleglazed IGU 4164164 is the most commonly produced and sold glazing for windows. Figure 6 presents an analysis illustrating the dependence of the percentage change in the calculated heat loss ∆Φm for these units on their width (at a constant ratio b/a = 2), under different external temperature conditions te. Simultaneous pressure increase ∆p = 3 kPa was assumed. Other data was used as in previous examples.
It was found that the described effect is important for the currently sold glazing in “mild winter” conditions, i.e., when the outside temperatures fluctuate within between −5 °C and 5 °C. For IGU width above 0.7m, the ratio ∆Φm changes from 3.9% to 5.0%. These values are characteristic of the average temperature during the winter months in many places around the world.
5. Notes on IGUs Wind Load
Wind pressure or suction are also factors that cause deflection of the component glass panes in an IGU. As already mentioned, the wind velocity pressure acts directly only on the outer pane, but due to the change in the gas pressure in the gaps, the resultant load is distributed over all the panes of the unit. Table 7 shows the resultant loads and deflections in sample unit’s surface loaded with 0.3 kN/m², which approximately corresponds [20] to a pressure of wind with velocity V of approx. 80 km/h (22.2 m/s).
Table 7. Static quantities in IGUs loaded with wind pressure of 0.3 kN/m².
Table 7 demonstrates that in the majority of units the deflections of component glass have similar values. Greater variations may occur when thicker panes are used, but the deflection values are small. It can therefore be concluded that the change in the thickness of the gaps due to wind load is small and has no noticeable effect on heat loss by IGUs.
Wind velocity has an indirect effect on heat loss. It is a factor influencing external thermal surface resistance on the outside, which translates into the Ugvalue. Graphic illustration of this effect is shown in Figure 7. The calculations were made for units with gap thickness of 16 mm. It can be noted that in the case of tripleglazed IGUs, the effect of wind velocity is negligible.
6. Conclusions
One of the factors influencing thermal transmittance Ug of insulating glass units is the thickness of gasfilled tight gaps. It is assumed in the calculation procedures that this thickness is not dependent on temporary changes in climatic factors. The thickness is variable under real operating conditions. In winter conditions in particular, IGU component glass panes take a concave form of deflection, which reduces the thickness of the gaps. This effect increases if the atmospheric pressure increases at the same time.
Based on the example calculations carried out, it has been shown that the increase in the calculated heat losses associated with the reduction of the gap thickness occurs when the conditions in the gap lead to Nu < 1, i.e., when the thermal transmittance of the gas layer is linearly dependent on its thickness. Heat losses can then increase to about 4.6% for doubleglazed IGUs and to about 5% for tripleglazed ones, for external air temperature te = 0 °C. These values almost do not depend on the nominal thickness of the gaps, which results from the linear dependence of static quantities in an IGU on this thickness. Under certain conditions, heat losses calculated according to standard procedures may therefore be underestimated.
It is different in the nonlinear range of the Ugvalue change (Nu > 1), i.e., when the outside temperature drops significantly or the gaps are thick enough. The thermal performance of glazing does not deteriorate. It is therefore advantageous to design IGUs so that Nu > 1, but it is necessary to take into account local climatic conditions and analyze loads that may also occur during the summer period.
In the case of the most commonly sold tripleglazed units 4164164 heat losses may be underestimated when the outside temperatures fluctuate between −5 °C and 5 °C. For large IGU dimensions, the ∆Φm index totals then from 3.9% to 5.0%.
It was also shown that the effect of wind load on gap thickness change is negligible in the context of heat loss estimation.
References
 Huang, S.; Wang, Z.; Xu, J.; Lu, D.; Yuan, T. Determination of optical constants of functional layer of online LowE glass based on the Drude theory. Thin Solid Films 2008, 516, 3179–3183. [Google Scholar] [CrossRef]
 Solovyev, A.A.; Rabotkin, S.V.; Kovsharov, N.F. Polymer films with multilayer lowE coatings. Mater. Sci. Semicond. Process. 2015, 38, 373–380. [Google Scholar] [CrossRef]
 Bernier, M.; Bourret, B. Effects of Glass Plate Curvature on the UFactor of Sealed Insulated Glazing Units. ASHRAE Trans. 1997, 103, 4038. [Google Scholar]
 Hart, R.; Goudey, H.; Arasteh, D.; Curcija, D.C. Thermal performance impacts of centerofglass deflections in installed insulating glazing units. Energy Build. 2012, 54, 453–460. [Google Scholar] [CrossRef][Green Version]
 Penkova, N.; Krumov, K.; Zashkova, L.; Kassabov, I. Heat transfer and climatic loads at insulating glass units in window systems. Int. J. Adv. Sci. Eng. Technol. 2017, 5, 22–28. [Google Scholar]
 Derwiński, W. Okna w kamerze termowizyjnej. Świat Szkła 2014, 12, 24–26. [Google Scholar]
 Solvason, K.R. Pressures and Stresses in Sealed Double Glazing Units; Technical Paper No. 423; Division of Building Research, National Research Council Canada: Ottawa, ON, Canada, 1974.
 Feldmeier, F. Klimabelastung und Lastverteilung bei MehrscheibenIsolierglas. Stahlbau 2006, 75, 467–478. [Google Scholar] [CrossRef]
 Curcija, C.; Vidanovic, S. Predicting Thermal Transmittance of IGU Subject to Deflection; Lawrence Berkeley National Laboratory, Environmental Energy Technologies Division: Berkeley, CA, USA, 2012.
 Feldman, M.; Kaspar, R.; Abeln, B.; Gessler, A.; Langosch, K.; Beyer, J.; Schneider, J.; Schula, S.; Siebert, G.; Haese, A.; et al. Guidance for European Structural Design of Glass Components; Support to the Implementation, Harmonization and Further Development of the Eurocodes; Luxembourg Publications Office of the European Union: Luxembourg, 2014.
 Stratiy, P. NumericalandAnalytical Method of Estimation Insulated Glass Unit Deformations Caused by Climate Loads. In International Scientific Conference Energy Management of Municipal Transportation Facilities and Transport EMMFT 2017; Murgul, V., Popovic, Z., Eds.; Advances in Intelligent Systems and Computing 2017; Springer International Publishing: New York, NY, USA, 2017; Volume 692, pp. 970–979. [Google Scholar]
 Velchev, D.; Ivanov, I.V. A finite element for insulating glass units. In Challenging Glass 4 & COST Action TU0905 Final Conference 2014; Louter, C., Bos, F., Belis, J., Lebet, J.P., Eds.; Taylor & Francis Group: London, UK, 2014; pp. 311–318. [Google Scholar]
 Respondek, Z. Rozkład obciążeń środowiskowych w wielokomorowej szybie zespolonej. Constr. Optim. Energy Potential 2017, 1, 105–110. [Google Scholar] [CrossRef]
 Klindt, L.B.; Klein, W. Glas als Baustoff: Eigenschaften, Anwendung, Bemessung; Verlagsgesellschaft R. Müller: KölnBraunsfeld, Germany, 1997. [Google Scholar]
 Timishenko, S.; WoinowskyKrieger, S. Theory of Plates and Shells; McGrawHill Book Company: New York, NY, USA; Toronto, ON, Canada; London, UK, 1959. [Google Scholar]
 EN 673:2011. Glass in Building—Determination of Thermal Transmittance (U Value)—Calculation Method; European Committee for Standardization (CEN): Brussels, Belgium, 2011. [Google Scholar]
 Aydın, O. Conjugate heat transfer analysis of double pane windows. Build. Environ. 2006, 41, 109–116. [Google Scholar] [CrossRef]
 Arıcı, M.; Kan, M. An investigation of flow and conjugate heat transfer in multiple pane windows with respect to gap width, emissivity and gas filling. Renew. Energy 2015, 75, 249–256. [Google Scholar] [CrossRef]
 EN 5721:2004. Glass in Building—Basic Soda Lime Silicate Glass Products—Part 1: Definitions and General Physical and Mechanical Properties; British Standards Institution: London, UK, 2004. [Google Scholar]
 Dynamic Pressure. Available online: https://en.wikipedia.org/wiki/Dynamic_pressure (accessed on 7 September 2019).