Next Article in Journal
Characterization of Thermal Expansion Coefficient of 3D Printing Polymeric Materials Using Fiber Bragg Grating Sensors
Previous Article in Journal
On the Creation of a Material Bond between L-PBF-Manufactured AZ91 and Ti-6Al-4V Components in the Context of Medical Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Elastic Support of Component Glass Panes on Deflection and Stress in Insulating Glass Units—Analytical Model

by
Zbigniew Respondek
Faculty of Civil Engineering, Czestochowa University of Technology, Dąbrowskiego 69, 42-201 Częstochowa, Poland
Materials 2024, 17(18), 4669; https://doi.org/10.3390/ma17184669
Submission received: 11 August 2024 / Revised: 19 September 2024 / Accepted: 21 September 2024 / Published: 23 September 2024
(This article belongs to the Section Construction and Building Materials)

Abstract

:
Insulating glass units (IGUs) are the most common filling for external building envelopes. These elements have many advantages related to the thermal protection of buildings. However, some climatic loads are generated or modified due to the sealed gas cavity between the glass panes. The gas enclosed in the cavities changes its parameters under external load, which affects the operational deflection and stress in an IGU. In most computational models describing this phenomenon, the component panes are assumed to be simply supported on the edge spacer, which is considered a sufficient approximation. This article, which continues previous work, assumes that the component glass panes can be supported elastically at the edges. The parameter describing this connection is rotational stiffness. Based on the theory of linear–elastic plates, coefficients were determined to calculate the change in cavity volume, deflection, and stress in glass panes. Then, the results of calculations of the influence of rotational stiffness and static values in exemplary IGUs of various structures, loaded with changes in atmospheric pressure and wind, are presented. It was found that a feedback loop occurs here. The deflection and stress in elastically supported single panes are lower than in the case of those simply supported. However, the lower susceptibility to deflection of the component panes weakens the gas interaction in the cavity, and the resultant load on these panes increases. The influence of rotational stiffness on the resulting static values may therefore vary. In the analyzed examples, this influence was primarily negative for symmetrical loads and clearly positive for wind loads.

1. Introduction

Insulating glass units (IGUs) are a commonly used filling for windows, glass façades, and other external transparent building envelopes. Their characteristic feature is sealed gas-filled cavities located between the glass component panes of the IGU. The idea behind this structure was to limit heat loss through transparent envelopes. In this context, they are constantly being improved, especially multi-glazed systems [1,2,3] or facades with closed cavities [4]. The tightness of the cavities is achieved by using edge spacers [5], most often with double sealing (Figure 1).
However, it turns out that sealing the cavities in IGUs has a certain side effect. Under the influence of climatic factors, the gas entrapped in the gap changes its parameters, which generates loads and deflections of the component glass panes. For example, an increase in atmospheric pressure or a decrease in gas temperature results in a concave form of deflection of the IGU; the reverse action results in a convex form. This is visually observable through the image distortion reflected from glazing, for example, glass facades [6].
The problem of static values in IGUs (both double- and multi-glazed) is well known. The literature describes relevant analytical models. The most famous and often cited articles are those by Feldmeier [7,8], who assumed that the gas in cavities changes its parameters (pressure, temperature, volume) in accordance with the ideal gas law, while the deflection and stress in glass panes can be estimated based on Kirchhoff’s elastic plate theory. Similar assumptions were made in many other analytical and numerical models [6,9,10,11,12,13]. Another approach to using plate theory, based on Betti’s Analytical Method (BAM) and Green’s functions, was used by Galuppi and Royer–Carfagni [14,15,16]. In the engineering approach, the calculation of deflection and stress in IGUs loaded with climatic factors can be estimated according to the European standard [17]. A comparison of selected analytical methods and a proposal to improve calculation efficiency are presented in [18].
Experimental research in this area is difficult and expensive; static quantities in IGUs depend on many factors, often uncertain and changing over time. Nevertheless, the literature describes such studies, most of which have been conducted recently. Tests of the deflection of IGU components in operational conditions are described in [19] and in simulated pressure or temperature changes in [10,20,21]. In [22], Galuppi et al. experimentally verified the BAM model [14] for double-glazed IGUs, in the range of line-distributed and concentrated loads and also at high temperature. Studies of external load sharing in double- and triple-glazed IGUs are described in [23,24], and tests and numerical analyses on this topic, using IGUs with laminated glass, are presented in [25,26]. In [27], this problem was analyzed numerically, taking into account the effect of large deflections, and in [28] for different shapes of IGUs with varying thicknesses of the component glass. Bedon and Amadio [29,30] developed numerical models of sample IGUs, taking into account all seal elements, using previously experimentally tested parameters of these elements.
The computational practice, confirmed by the experimental studies cited above, reveals some general trends regarding the deflection and stress values in IGUs. In the case of loads acting symmetrically, for example, a change in external atmospheric pressure, interaction of the gas confined in the cavities occurs. Due to the change in the volume of this cavity caused by the deflection of the component glass panes, the gas pressure changes, which partially compensates for the external load. This compensation is greater when the component panes are more susceptible to deflection. As a consequence, the resultant load of the component panes, i.e., the pressure difference between the atmospheric air and the cavity, decreases when the linear dimensions (width, length) of the IGUs increase. It has also been shown [7,8] that for each IGU structure, there is a critical width at which the stress in the component panes is maximum. On the other hand, the use of thicker component panes, as well as increasing the distance between the panes (or increasing the total thickness of the cavities), increases the resultant load [31]. A significant increase in load was also observed when IGUs were used with curved glass [32]. In Cwyl et al. [33], reported a case where a composite spacer was “pulled” into the cavity of a bent IGU, probably due to high underpressure generated under winter conditions.
The situation is different in the case of uniformly distributed loads (e.g., wind action). There is a favorable gas coupling of the IGU component panes. The external load is distributed over all the panes, with large IGU dimensions approximately proportional to the stiffness of the component panes [34].
In most analytical and numerical models, as well as in the standard [17], as a sufficiently accurate approximation, the simple support of the component panes at the edges is assumed, i.e., zero support moments at the connection of the component panes with the edge spacer. It is believed that such an assumption is correct from the point of view of structural safety [14], i.e., taking into account the possibility that elastic support leads to a reduction in the maximum stresses in the IGU, which was shown in experimental studies by Buddenberg et al. [20]. The possibility of elastic support for the component panes was also demonstrated in a study by Stratiy [10]. The IGU model embedded in a window frame and loaded with a pressure difference was tested here. The tests indicated that the deflection of the component slab has intermediate values between those estimated for a simply supported and restrained pane. Similar tests (an IGU was placed in a structure simulating a window frame) described in [21] showed that the measured deflection and stress in the center of the panes differed from the theoretical values for simple support, which may suggest elastic support. Regarding the theoretical analysis of the problem, recently, Galuppi and Zacchei [35] published a theoretical and numerical analysis on the transfer of external load between the panes of double-glazed IGUs, assuming their support with elastic beams, depending on the coefficients describing the bending and torsional stiffness of the support beams.
However, the studies described in the literature do not cover all cases of IGU design and possible variants of climatic loads. It is unquestionable that for elastic support, the deflection and stress in the center of a pane are greater for simple support than for elastic support, but provided that the load on this pane remains constant. Meanwhile, the IGU behaves in a specific way—it is expected that reducing the susceptibility to deflection of the component panes with elastic support increases the resultant load. To estimate possible effects, it is necessary to check the behavior of IGUs with different structures and dimensions. Due to the above, it was noticed that there was a need to supplement the author’s computational model [6,11] with the case of elastic support. The previous article [36] proposed an analytical model enabling the estimation of the resultant load and deflection for IGUs with elastically supported component panes, assuming any rotational stiffness. Therefore, this work continues the analysis presented in [36].
Accordingly, this article aims to complement the computational model by estimating stress in IGUs with elastically supported component panes and then to analyze the influence of rotational stiffness on static values in IGUs with various geometric parameters.
The “Materials and Methods” section presents the proposed analytical model for estimating static quantities in IGUs and its numerical validation. The “Results” section contains an analysis of the influence of rotational stiffness on static values in IGUs of various designs. The “Discussion” section explains some issues regarding the practical application of the proposed model and analyzes what factors may influence the rotational stiffness in real IGUs. The article ends with “Conclusions”.

2. Materials and Methods

2.1. Base Model

As mentioned, the author’s base model for estimating static quantities in climate-loaded IGUs was described in [6,11]. The model assumes conformity of changes in gas parameters in the cavity with the ideal gas law, a linear relationship between the deflection and stress of the component panes, and the simple support of these panes at the edges. The model allows the estimation of static quantities for IGUs with any number of cavities, loaded with pressure and temperature changes and uniform loads per area. The basis of the calculations is to determine the operational gas pressure in individual cavities pop [kPa] and then, based on the load balance, the resultant load q [kPa] for each component pane. The problem was reduced to solving an appropriate quadratic equation (double-glazed IGU) or a system of quadratic equations (multi-glazed IGU). In the simplest case (Figure 2) of a double-glass IGU loaded with pressure changes, temperature changes, and wind, the operating pressure in the cavity is
p op = A 2 B + A 2 B 2 + p 0 v 0 T op B T 0
with
A = ( p a + q z , e x ) α v , ex + p a α v , in v 0
B = α v , ex + α v , in
where
  • p0, T0, v0—initial parameters of gas in the cavity: pressure [kPa], temperature [K] and volume [m3]; parameters of the production process; it is assumed that under these conditions, the resultant load equals 0;
  • pop, Top—operating parameters: pressure [kPa] and temperature [K] in the cavity;
  • pa—current atmospheric pressure [kPa];
  • αv,ex, αv,in—proportionality coefficients describing the change in the volume of the cavity caused by the deflection of a given glass under a unit load per area [m5/kN];
  • qz,ex—wind load [kN/m2].
In the above formulas, the index “ex” denotes the pane on the outside air side, and the index “in” denotes the pane on the room side. If necessary, a sign convention is used—load and deflection values are considered to be positive if they face the interior (Figure 2).
The proportionality coefficients αv were determined on the basis of the deflection function of a simply supported plate [37] according to the following formula:
α v = α v a 6 D       w i t h       D = E d 3 12 1 μ 2
where
  • αv—dimensionless coefficient [−];
  • a—width of the glass pane [m];
  • D—flexural rigidity of glass pane [kNm];
  • E—Young’s modulus of glass [kPa];
  • d—glass pane thickness [m];
  • μ—Poisson’s ratio of glass [−].
After calculating the resultant load q [kPa], the deflection in the center of the glass pane wmax [mm] can be calculated using the following formula:
w m a x = α w q a 4 D · 1000
where
  • αw—dimensionless coefficient [-].
The article [11] also provides dimensionless coefficients for calculating stress in glass panes.

2.2. Deflection and Change of Cavity Volume in IGUs with Elastically Supported Glass Panes

Previous work [36] described an analytical model, taking into account the elastic support of component panes around the perimeter for estimating the resultant load q and maximum deflection wmax in climatically loaded IGUs. The problem was reduced to determining dimensionless coefficients αv and αw for any rotational stiffness of the connection of the glass pane with the spacer (support beam). Further calculations are carried out in accordance with the base model. The analysis of an elastically supported plate is modeled on the rigidly clamped plate solution given by Timoshenko and Woinowsky–Krieger [37], except that different boundary conditions were used.
To determine the coefficients αv and αw, the change in the cavity volume caused by the deflection of one of the component panes was analyzed, assuming that the rotational stiffness of the connection C [kNm/(m·rad)] is known. It was also assumed that the deflection of a pane with dimensions a × b (width × length) is the algebraic sum of the deflection of a simply supported pane loaded with a uniform load q and edge supporting moments (My)y=±b/2 and (Mx)x=±a/2 (Figure 3); see Equations (3)–(6) in [36] based on [37]. The unknown supporting moments were given as follows:
M y y = ± b / 2 = i = 1,3 , 5 ( 1 ) ( i 1 ) / 2 F i cos i π x a
M x x = ± a / 2 = i = 1,3 , 5 ( 1 ) ( i 1 ) / 2 G i cos i π y b
Fi [kNm/m] and Gi [kNm/m] are strongly decreasing numerical sequences whose elements should be determined from the boundary conditions. In this case,
w y = M y C   f o r   y = ± b / 2                 w x = M x C   f o r   x = ± a / 2
The analysis [36] results in the construction and solution of an appropriate system of linear equations in which the first five elements of both sequences are unknown. As a solution, it is convenient to provide dimensionless values of Fi and Gi, such that
F i = F i · q · a 2               G i = G i · q · a 2
For practical reasons, another dimensionless quantity was also defined:
R c = 2 π D a C
Thanks to the use of this parameter, Fi and Gi depend only on the Rc value and the IGU aspect ratio ε = b/a. For the boundary values, C = 0, Rc  ∞ means simple support, and C  ∞, Rc = 0 means clamped fixity.
The Fi and Gi coefficients for selected Rc values and aspect ratios ε = b/a are presented in File S1 (Supplementary Materials). Knowing these values makes it possible to determine the coefficients αw according to Equations (3)–(6) from [36] and then, after integrating these equations, calculate the coefficients αv. Tables with calculated dimensionless coefficients are presented in File S2 (Supplementary Materials).
It should be noted that for the boundary values (Rc  ∞, Rc = 0), the coefficients F′ and G′ (see File S1) describing the supporting moments and the coefficients αw regarding the maximum deflection are consistent (differences do not exceed 0.5%) with tables presented by Timoshenko and Woinowsky–Krieger in [37].

2.3. Stress in Elastically Supported Glass Panes

The stress at the center of the glass pane was estimated using well-known equations from the theory of plates. In general,
σ x = 6 · m x d 2           σ y = 6 · m y d 2
where
m x = D 2 w x 2 + μ 2 w y 2           m y = D 2 w y 2 + μ 2 w x 2
In this case, we use the principle of superposition:
w = w q + w M y + w M x
where
  • wq, wMy, wMx—deflection of a simply supported pane under load q and deflections due to supporting moments My, and Mx [kNm/m].
After integrating the appropriate functions of deflection (see Equations (4)–(6) in [36], based on [37]), for the center of the plate, we obtained
2 w q x 2 x = 0 ,     y = 0 = 4 · q · a 2 D · π 3 i = 1,3 , 5 , 1 i 1 / 2 i 3 1 2 + β i t h β i 2 c h β i
2 w q y 2 x = 0 , y = 0 = 2 · q · a 2 D · π 3 i = 1,3 , 5 , 1 i 1 / 2 i 3 β i t h β i c h β i
2 w M y x 2 x = 0 , y = 0 = 1 D i = 1,3 , 5 . ( 1 ) ( i 1 ) / 2 · F i · β i th β i 2 ch β i
2 w M y y 2 x = 0 , y = 0 = 1 D i = 1,3 , 5 . ( 1 ) ( i 1 ) / 2 · F i · 2 β i th β i 2 ch β i
2 w M x x 2 x = 0 , y = 0 = 1 D i = 1,3 , 5 . ( 1 ) ( i 1 ) / 2 · G i · 2 γ i th γ i 2 ch γ i
2 w M x y 2 x = 0 , y = 0 = 1 D i = 1,3 , 5 . ( 1 ) ( i 1 ) / 2 · G i · γ i th γ i 2 ch γ i
In Equations (14)–(19), auxiliary quantities β i = i π ε 2 and γ i = i π 2 ε were used.
After inserting the appropriate values of partial derivatives into Equation (12) and then into Equation (11), it can be seen that the stress values in the pane center σx,cen [MPa], σy,cen [MPa] can be expressed using dimensionless coefficients kx,cen i ky,cen. Then,
σ x , c e n = k x , c e n 6 · q a 2 d 2 / 1000               σ y , c e n = k y , c e n 6 · q a 2 d 2 / 1000
Dimensionless coefficients for calculating stress for selected Rc and ε = b/a values are presented in File S3 (Supplementary Materials). Calculations were made for two variants of the Poisson’s ratio: μ = 0.2 (according to [38]) and μ = 0.23 (according to [17]).

2.4. Reference IGU Parameters

For the numerical validation of the presented analytical model (Section 2.5) and the analysis of the influence of selected factors on static quantities (Section 3), a reference IGU with the following parameters was adopted: double-glazed IGU, dimensions a × b = 60 × 90 cm2 (ε = 1.5), thickness of glass panes dex = din = 4 mm, cavity thickness s = 16 mm. Initial conditions (state without deformations and stresses): p0 = 100 kPa, T0 = 293.15 K (20 °C). Glass parameters: E = 70 GPa, μ = 0.2.

2.5. Numerical Validation of the Analytical Model

In order to validate the computational model, the results of analytical calculations of a single glass plate (ANA) were compared with numerical calculations (NUM). The following static values were analyzed: maximum deflection wmax, the stress in the plate center σx,cen; σy,cen and edge stress σx,ed; σy,ed.
The numerical model was designed in Autodesk Robot StructuralTM using the “plate” structure type (Figure 4). The computational model was defined as a “shell”. Geometric dimensions and material data were determined according to Section 2.4. Linear supports were assumed on all edges. The support conditions are presented in Table 1. The plate was uniformly loaded with q = 5 kPa in the z direction. Due to the assumed vertical position of the plate, the self-weight load was omitted. The finite element mesh of size 1 × 1 cm2 was created using the simple (Coons) four-node method.
The calculation results are summarized in Table 2 and Table 3. Sample maps of static quantities are presented in File S4 (Supplementary Materials).
Based on the data presented in Table 1 and Table 2, it was found that for the key parameters describing the deflection and stress in the analyzed plate, the analytical results are in good agreement with the numerical ones, which proves the mathematical correctness of the proposed model. A similar agreement was obtained when analyzing other plate thicknesses and dimensions. Significant discrepancies were found only for edge stresses in the case of low C-values (marked gray). This is related to the failure to include in the analytical model the stresses from the reactions preventing the displacement of the plate corners during bending [37]. These additional stresses have a noticeable effect in the case of Rc > 2. Due to these discrepancies, coefficients for calculating edge stress are not given in File S3 for Rc > 2 kNm/(m·rad), and results for C < 2 kNm/(m·rad) are not given in the graphs in the “Results” section.
However, the adopted simplification does not affect the further course of the analysis. In cases of close to simple support, the maximum stress (in terms of absolute value) occurs in the plate center and is much greater than at the edge. With the increase in rotational stiffness, the influence of edge stresses is increasingly stronger. In the analyzed plate, for C ≥ 5 kNm/(m·rad), the stresses σx are greater at the edges than in the center of the plate.

3. Results

3.1. Scope of Analysis

The analysis of the influence of rotational stiffness on static quantities in IGUs was performed in such a way that, based on the reference unit described in Section 2.4, one of the following parameters was changed: width of the glass pane a (0.4–1.2 m), aspect ratio ε (1.0–2.0), glass thickness d (4–10 mm), and cavity thickness s (12–18 mm). Double-glazed IGUs and two types of loading were assumed: symmetrical (change in atmospheric pressure) and asymmetrical (wind pressure). The static quantities analyzed included load reduction r [%], maximum deflection (in the center of the IGU) wmax [mm], and maximum stress (in the IGU center σx,cen [MPa], at the edges σx,ed [MPa]). Stresses σy were not analyzed because they are always smaller than σx (or equal to ε = 1.0). The resulting graphs illustrate the absolute values of the analyzed quantities. To demonstrate the increase in the resultant load with increasing C, the load reduction values are arranged on the ordinate in reverse order.
Based on the results of the tests presented in [10] (described in more detail in Section 4), the range of variation C = 0–5 kNm/(m·rad) was assumed. It was therefore considered that higher C-values are unlikely for real IGUs. To simplify the description of the results, the unit with simply supported panes was designated IGU(C0), while for C = 5 kNm/(m·rad), the symbol IGU(C5) was used.

3.2. Symmetrical Load—Atmospheric Pressure Changes

The load of atmospheric pressure drop of ∆p = 5 kPa was assumed, i.e., from p0 = 100 kPa to pa = 95 kPa. Such pressure fluctuations are possible due to natural changes in weather conditions. It should be noted that similar effects would be caused by an increase in gas temperature in the cavity (for the reference IGU, the equivalent temperature increase is 15.4 K).
The reduction of the symmetrical load rsym was calculated using the following formula:
r s y m = 1 q p · 100 %
where q [kPa] is the resultant load, identical for both glass panes.
This reduction is related to the interaction of the gas enclosed in the cavity. It can be seen that as the sizes of IGUs increase, rsym tends toward 100%. The calculation results are shown in Figure 5, Figure 6, Figure 7 and Figure 8. Spreadsheet printouts for all analyzed cases are presented in File S5 (Supplementary Materials).
As mentioned, the elastic support of the component panes at the edges leads to feedback. With the increase in rotational stiffness, the susceptibility of the panes to deflection decreases, which limits the reduction of the external load, i.e., the resultant load is higher. This is apparent in Figure 5a, Figure 6a, Figure 7a and Figure 8a. This phenomenon can affect the resulting values of the maximum deflection and stress in the IGUs in different ways.
Figure 5b,c indicates that increasing rotational stiffness is disadvantageous for large dimensions of IGUs. For example, for a = 120 cm, the deflection in IGU(C5) is 16.6% greater than in IGU(C0), and the stress in the glass pane center is even 58.7% greater. Additionally, edge stress begins to dominate at C = 2 kNm/(m rad). Fortunately, in large dimensions of IGUs, the stress values are not significant. As is generally known, the greatest stresses occur near the so-called characteristic (critical) widths. For the analyzed IGU, the critical width calculated on the basis of [8] is 37.2 cm. It can be seen that near this value (calculations for a = 40 cm), the stress in the pane center for IGU(C5) is only 2.0% greater than in IGU(C0), and the deflection is 15.0% smaller.
The analysis of different aspect ratios showed that the most advantageous variant for deflection is the square IGU (Figure 6b). The deflection in IGU(C5) is 2.2% smaller than in IGU(C0). The change of ε-value has no significant influence on the resulting dependencies concerning stress (Figure 6c). Edge stresses dominate for C values greater than 4 kNm/(m rad).
Figure 7b,c demonstrates that increasing rotational stiffness is beneficial for thicker glass panes. For example, for d = 10 mm, the deflection in IGU(C5) is 5.5% smaller than in IGU(C0), and the stress in the pane center is 2.9% smaller. The edge stress is significantly smaller. It can also be seen that the stress values are the highest for glass panes with a thickness of 6–8 mm.
Finally, the analysis of different gas cavity thicknesses (Figure 8b,c) details that this parameter does not have a significant effect on the percentage change in static values with increasing C. This is due to the fact that the impact of the cavity thickness on the deflection and stress in the IGUs is approximately linear. So, with increasing thickness, the deflection and stress increase, but the shape of the resulting graphs is similar.

3.3. Asymmetrical Load—Wind Pressure

A wind pressure load of 0.3 kN/m2, acting on the glass pane “ex”, was assumed, which corresponds to wind speed of about 80 km/h [39]. The reduction of this load rwin was calculated using the following formula:
r w i n = 1 q e x q z , e x · 100 %
where qex [kPa] is the resultant load acting on the glass pane “ex”.
This reduction is related to the transmission of part of the load from the glass pane “ex” to the glass pane “in” via the gas cavity. It can be seen that when increasing the IGU dimensions, the reduction of this load tends toward 50%, i.e., the external load is distributed on both panes in half. It should be added that in the case of different thicknesses of the component glass panes, the loads are distributed differently—proportionally to the stiffness of the component glass panes.
The analysis results are presented in Figure 9, Figure 10, Figure 11 and Figure 12. Spreadsheet printouts for all analyzed cases are presented in File S6 (Supplementary Materials).
First, it can be observed that for the load acting only on the “ex” glass pane, the deflection and stress in the pane center increase with increasing the IGU dimensions or decreasing the thickness of the component panes. Therefore, there is no critical width here. Secondly, the load reduction (Figure 9a, Figure 10a, Figure 11a and Figure 12a) is approximately half of that in symmetric loads (Figure 5a, Figure 6a, Figure 7a and Figure 8a). As a result, the increase in the resultant load caused by increasing the rotational stiffness is not as significant. Consequently, the increase in the C-value in each case considered leads to decreased deflection and stress in the pane center.
For different dimensions of IGUs (Figure 9b,c), within the range of the analyzed data, the deflection in IGU(C5) is 38.9 –57.6% smaller than in IGU(C0), and the stress in the pane center drops by 26.7–43.4%. However, it should be noted that at certain C-values, the edge stress starts to dominate (e.g., for a = 10 cm, this limit value of C is about 2.5 kNm/(m·rad).
For different aspect ratios (Figure 10b,c), the deflection in IGU(C5) is 42.1–53.4% smaller than in IGU(C0), and the stress in the pane center drops by 29.5–40.7%. For different pane thicknesses (Figure 11b,c), these values are 8.6–48.4% for deflection and 6.1–35.1% for stress. The change in the gas cavity thickness (Figure 12b,c) has almost no effect on the static values in IGU. The deflection drops by about 48–49%, and the stress in the pane center by about 35–36%. Edge stresses start to dominate at a C greater than about 4 kNm/(m rad).

4. Discussion

Estimating static values in IGUs is a complex issue. The efficiency of the interaction of the gas enclosed in the cavity depends on many factors related to the unit’s design and the assumed load combination. The proposed model introduces an additional factor of influence, namely the possibility of connecting the component panes with a spacer other than a simple support. It was assumed that the parameter describing this connection is the rotational stiffness C. Then, the methodology was developed to determine several dimensionless coefficients and apply them to the basic model describing IGUs with simply supported glass panes. The advantage of the model is the possibility of performing calculations using generally available spreadsheets, e.g., Microsoft ExcelTM (2016). Due to the mathematical complexity of the equations, setting up and testing such a spreadsheet requires quite some effort; however, after completing this task, the results are obtained automatically. A certain disadvantage of the model is the limitation of its application to the range of linear–elastic deflections of the glass pane.
This article analyzed only the most straightforward cases of IGU structures and acting loads to show general trends. However, it is not difficult to analyze any combinations of loads (wind, climatic loads) for multi-glazed IGUs with different stiffness of the component panes. It is also possible to declare different rotational stiffness for the individual panes of the unit.
Taking into account the conducted analysis in the literature, it can be noticed that the issue of rotational stiffness of the connection of component panes with the edge spacer in real IGUs is poorly understood; in particular, there are few experimental studies in this area. Based on original research [20] and earlier studies [40], Buddenberg et al. estimated that the rotational stiffness of the connection should not be greater than 0.2 kNm/(m·rad)—this is a state close to simple support. Also, Kozłowski et al. [32] show agreement of with the measured values of forced overpressure in the cavity with the values calculated for simple support. It can be noted, however, that these tests were conducted using IGU samples that were not embedded in the window frame.
As mentioned, the deflection tests in the IGU embedded in the frame were carried out by Stratiy [10]. Unfortunately, the author did not provide exact numerical values, but presented the results as a graph of the dependence of deflection on the value of forced load (Figure 6 in [10]). In the linear–elastic range, the measured deflection takes intermediate values between the estimated ones for simple support and restraint. Based on approximate values read from that graph, taking into account the parameters given in the article (dimensions 100 × 100 cm2, glass thickness 4 mm), the rotational stiffness was estimated at about 3.5 kNm/(m·rad). For this reason, in Section 3, it was assumed that a slightly higher value is possible for real IGUs in operating conditions.
In the context of the analyzed issue, attention should be paid to the variety of edge spacer designs used [5]. The classic solution is an aluminum or steel profile with a closed cross-section filled with a molecular sieve. Spacers with improved thermal insulation made of plastic with aluminum foil have a similar shape [41]. In most of these solutions, double sealing is used, e.g., polyisobutylene at the contact between the glass pane and the spacer (primary sealant) and a polyurethane or silicone sealant as a filling at the IGU edges (secondary sealant). Composite spacers made of structural foam [42] or a PIB/silicone system [33] are a separate type. These connections have viscoelastic properties [43], i.e., their efficiency depends on factors such as temperature, load duration, number of load repetitions, etc. Also, the aging of sealants may affect their properties [44].
Another factor that can be assumed to affect the rotational stiffness of the connection is the way the IGU is embedded in the frame. Various solutions are used here: rubber and silicone gaskets and bonded glazing systems [45]. Manufacturing and assembly errors may also have an impact—incorrectly selected dimensions of glass fillings mean that they are too shallowly, too deeply, or too tightly embedded in the frame [46]. Therefore, the freedom of rotation of the component glass on the support may be limited by the action of the window frame elements.
Taking the above into account, the author sees a need for experimental studies on the possibility of stiffening the connection of IGU component panes with different spacer structures, also using IGU samples embedded in the window frame. As shown in Section 3, under certain conditions, over-stiffening the connection can lead to increased deflection and stress in the IGUs.

5. Conclusions and Further Work

The paper implements the task of taking into account the possibility of elastic connection of component glass panes with the edge spacer to determine deflection and stress in IGUs. The proposed model is based on the previous analytical solution, which assumes simple support. Based on the Kirchhoff linear–elastic plate theory, appropriate coefficients were determined, which, when used in the basic model, allow estimation of the change in the gas cavity volume, deflection, and stress in IGUs, assuming any rotational stiffness of the connection of component panes with the spacer.
For a single plate loaded uniformly per area, the model was numerically validated for the entire range of C-values (from simply support to clamped fixity). Therefore, it can be helpful for the analysis of plates made of other homogeneous materials, not only glass.
The presented calculation examples concern the influence of rotational stiffness on static quantities for simple cases of climatic loads of double-glazed IGUs: change of atmospheric pressure (symmetrical load) and wind pressure (asymmetrical load). It was found that a kind of feedback loop occurs here—the deflection and stress in the center of a single elastically supported glass pane decrease with increasing rotational stiffness. On the other hand, the reduced susceptibility to deflection of the component panes weakens the interaction of the gas in the cavity, which leads to an increase in the resultant load acting on the component panes. As a result, the influence of rotational stiffness on the resulting static quantities can be different.
This influence was mostly negative for symmetrical loads in the analyzed examples (for C = 0–5 kNm/(m rad)). In some cases, the deflection and stress in the center were over 50% greater than for free support. However, this negative influence was insignificant for dimensions close to critical. A clearly positive influence of the connection stiffening was identified for wind load. In some cases, the static values were reduced by approx. 40 ÷ 50%. More complex load combinations and IGU structural arrangements require individual analysis.
The possibility of an elastic connection of component panes with the edge spacer in IGUs is poorly tested experimentally. Only a few studies are described in the literature. In particular, there is a lack of experiments regarding the possibility of stiffening this connection in IGUs embedded in the window frame, which can be mentioned as a direction for further research.
Of course, the analytical model requires continuous improvement. In the near future, the scope of the model will be extended to cover large deflections. This is very important because, in such cases, the susceptibility of the glass component panes to deflection decreases, which can increase the resultant load. In addition, it is planned to analyze the influence of certain phenomena on the mechanical behavior of IGUs, which have not been previously taken into account in analytical models: firstly, the possibility of thermal deflection of the component panes caused by the temperature difference on both surfaces of a single glass pane; secondly, the behavior of diagonally and horizontally located IGUs under radiative cooling conditions. This involves a temperature drop on the diagonal and horizontal surfaces of building envelopes below the outside air temperature on cloudless nights.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ma17184669/s1, File S1: Coefficients for Calculating Supporting Moments; File S2: Coefficients for Calculating Deflection and Change in Cavity Volume in IGUs; File S3: Coefficients for Calculating Stress in Glass Panes; File S4: Sample Maps of Static Quantities; File S5: Calculations for Symmetrically Loaded IGUs—Spreadsheet Printout; File S6: Calculations for Asymmetrically Loaded IGUs—Spreadsheet Printout.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available upon request.

Conflicts of Interest

The author declares no conflicts of interest.

Nomenclature

Aauxiliary parameter [m3]
awidth (of IGU) [m] or [cm]
Bauxiliary parameters [m5/kN]
blength (of IGU) [m] or [cm]
Crotational stiffness [kNm/(m·rad)]
Dflexural rigidity (of component glass pane) [kNm]
dthickness (of component glass pane) [m] or [mm]
EYoung’s modulus [kPa] or [GPa]
Fnumeric sequence [kNm/m]
Fnumeric sequence [−]
Gnumeric sequence [kNm/m]
Gnumeric sequence [−]
iconsecutive natural number
jconsecutive natural number
kdimensionless stress coefficient [−]
Msupporting moment (at the edges) [kNm/m]
mflexural moment [kNm/m]
ppressure (absolute) [kPa]
qresultant load (of component glass pane) or outside load [kN/m2] or [kPa]
rload reduction [%]
Rcdimensionless quantity [−]
sthickness (of gas cavity) [mm]
Tabsolute temperature (of gas in the cavity) [K]
vvolume (of the cavity) [m3]
wdeflection [m] or [mm]
x-y-zcoordinate system
Greek letters
αproportionality factor, [m5/kN]
αdimensionless coefficient [−]
βauxiliary parameter [−]
γauxiliary parameter [−]
pchange in atmospheric pressure [kPa]
vchange in cavity volume [m3]
εaspect ratio [−]
μPoisson’s ratio [−]
πnumber “pi”
σstress [kPa or MPa]
Subscripts and markings
0initial gas parameters
aatmospheric
cenin the center
edat the edges
exexterior glass pane
IGU(C0)IGU with simply supported panes
IGU(C5)IGU with elastically supported panes C = 5 kNm/(m·rad)
ininterior glass pane
maxmaximum
opoperating gas parameters
qresultant load
R1–R3reference points
S1–S4support beams
symsymmetrical loads
vcavity volume
wdeflection
winwind pressure
x-y-zdirection
zoutside

References

  1. Wakili, K.G.; Rädle, W.; Krammer, A.; Uehlinger, A.; Schüler, A.; Stöckli, T. Ug-value and edge heat loss of triple glazed insulating glass units: A comparison between measured and declared values. J. Build. Eng. 2021, 44, 103031. [Google Scholar] [CrossRef]
  2. Kralj, A.; Drev, M.; Žnidaršič, M.; Černe, B.; Hafner, J.; Jelle, B.P. Investigations of 6-pane glazing: Properties and possibilities. Energy Build. 2019, 190, 61–68. [Google Scholar] [CrossRef]
  3. Starman, B.; Maček, A.; Rus, P.; Obid, Š.; Kralj, A.; Halilovič, M. Primary Seal Deformation in Multipane Glazing Units. Appl. Sci. 2020, 10, 1390. [Google Scholar] [CrossRef]
  4. Veršić, Z.; Binički, M.; Nosil Mešić, M.; Galić, J. Passively Maintained Closed Cavity Façade—Experimental Validation of the Mathematical Thermal Model. Buildings 2023, 13, 2031. [Google Scholar] [CrossRef]
  5. Van Den Bergh, S.; Hart, R.; Jelle, B.P.; Gustavsen, A. Window spacers and edge seals in insulating glass units: A state-of-the-art review and future perspectives. Energy Build. 2013, 58, 263–280. [Google Scholar] [CrossRef]
  6. Respondek, Z. Heat Transfer Through Insulating Glass Units Subjected to Climatic Loads. Materials 2020, 13, 286. [Google Scholar] [CrossRef] [PubMed]
  7. Feldmeier, F. Klimabelastung und Lastverteilung bei Mehrscheiben-Isolierglas. Stahlbau 2006, 75, 467–478. [Google Scholar] [CrossRef]
  8. Feldmeier, F. Bemessung von Dreifach-Isolierglas. Stahlbau Spez. 2011—Glas./Glass Build. 2011, 80, 75–80. [Google Scholar] [CrossRef]
  9. Curcija, C.; Vidanovic, S. Predicting Thermal Transmittance of IGU Subject to Deflection; Lawrence Berkeley National Laboratory, Environmental Energy Technologies Division: Berkeley, CA, USA, 2012.
  10. Stratiy, P. Numerical-and-Analytical 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; Springer: Cham, Switzerland, 2017; Volume 692, pp. 970–979. [Google Scholar]
  11. Respondek, Z.; Kozłowski, M.; Wiśniowski, M. Deflections and Stresses in Rectangular, Circular and Elliptical Insulating Glass Units. Materials 2022, 15, 2427. [Google Scholar] [CrossRef]
  12. 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]
  13. Bedon, C.; Amadio, C. A linear formulation for the ULS design of glass elements under combined loads: Application to IGUs. Glass Struct. Eng. 2018, 3, 289–301. [Google Scholar] [CrossRef]
  14. Galuppi, L.; Royer-Carfagni, G. Betti’s Analytical Method for the load sharing in double glazed units. Compos. Struct. 2020, 235, 111765. [Google Scholar] [CrossRef]
  15. Galuppi, L.; Royer-Carfagi, G. Green’s functions for the load sharing in multiple insulating glazing units. Int. J. Solids Struct. 2020, 206, 412–425. [Google Scholar] [CrossRef]
  16. Galuppi, L. Practical expressions for the design of DGUs. The BAM approach. Eng. Struct. 2020, 221, 110993. [Google Scholar] [CrossRef]
  17. EN 16612:2020; Glass in Building. Determination of the Lateral Load Resistance of Glass Panes by Calculation. CEN: Brussels, Belgium, 2020.
  18. Halilovič, M.; Maček, A.; Mole, N.; Koc, P.; Plešnik, F.; Rus, P.; Žnidaršič, M.; Kralj, A. Accurate determination of the static equilibrium in insulating glass units under climatic loading. J. Build. Eng. 2023, 80, 107955. [Google Scholar] [CrossRef]
  19. Hart, R.; Goudey, H.; Arasteh, D.; Curcija, D.C. Thermal performance impacts of center-of-glass deflections in installed insulating glazing units. Energy Build. 2012, 54, 453–460. [Google Scholar] [CrossRef]
  20. Buddenberg, S.; Hof, P.; Oechsner, M. Climate loads in insulating glass units: Comparison of theory and experimental results. Glass Struct. Eng. 2016, 1, 301–313. [Google Scholar] [CrossRef]
  21. Kozłowski, M.; Respondek, Z.; Wiśniowski, M.; Cornik, D.; Zemła, K. Experimental and Numerical Simulations of Climatic Loads in Insulating Glass Units by Controlled Change of Pressure in the Gap. Appl. Sci. 2023, 13, 1269. [Google Scholar] [CrossRef]
  22. Galuppi, L.; Zacchei, E.; Esteves, M.; Ferrão, J.; Simões, N. Experimental validation of the Betti’s Analytical Method for Double Glass Units. Eng. Struct. 2024, 315, 118468. [Google Scholar] [CrossRef]
  23. McMahon, S.; Norville, H.S.; Morse, S.M. Experimental investigation of load sharing in insulating glass units. J. Archit. Eng. 2018, 24, 04017038. [Google Scholar] [CrossRef]
  24. Huveners, E.M.P.; Van Herwijnen, F.; Soetens, F. Load sharing in insulated double glass units. Heron 2003, 48, 99–122. [Google Scholar]
  25. Zacchei, E.; Simões, N.; Vieira, A.; Esteves, M.; Silva, H. Modelling of layers interactions on the structural behaviour of insulating glasses. Vertical deflection analyses. Case Stud. Constr. Mater. 2023, 18, e02129. [Google Scholar] [CrossRef]
  26. Wang, Z.; Liu, J.; Li, D.; Yang, K.; Chen, M.; Wang, C. Experimental and numerical study on load-bearing performance in triple-glazed insulating glass units. Constr. Build. Mater. 2024, 418, 135385. [Google Scholar] [CrossRef]
  27. Heiskari, J.; Romanoff, J.; Laakso, A.; Ringsberg, J.W. On the thickness determination of rectangular glass panes in insulating glass units considering the load sharing and geometrically nonlinear bending. Thin-Walled Struct. 2022, 171, 108774. [Google Scholar] [CrossRef]
  28. Heiskari, J.; Romanoff, J.; Laakso, A.; Ringsberg, J.W. Influence of the design constraints on the thickness optimization of glass panes to achieve lightweight insulating glass units in cruise ships. Mar. Struct. 2023, 89, 103409. [Google Scholar] [CrossRef]
  29. Bedon, C.; Amadio, C. Mechanical analysis and characterization of IGUs with different silicone sealed spacer connections—Part 1: Experiments. Glass Struct. Eng. 2020, 5, 301–325. [Google Scholar] [CrossRef]
  30. Bedon, C.; Amadio, C. Mechanical analysis and characterization of IGUs with different silicone sealed spacer connections—Part 2: Modelling. Glass Struct. Eng. 2020, 5, 327–346. [Google Scholar] [CrossRef]
  31. Respondek, Z. Influence of Insulated Glass Units Thickness and Weight Reduction on their Functional Properties. Open Eng. 2018, 8, 455–462. [Google Scholar] [CrossRef]
  32. Kozłowski, M.; Respondek, Z.; Wiśniowski, M.; Cornik, D.; Zemła, K. Influence of curvature and geometrical parameters on internal pressure in cylindrical Insulating Glass Units. Thin-Walled Struct. 2023, 188, 110812. [Google Scholar] [CrossRef]
  33. Cwyl, M.; Michalczyk, R.; Wierzbicki, S. Polyisobutylene and Silicone in Warm Edge Glazing Systems—Evaluation of Long-Term Performance. Materials 2021, 14, 3594. [Google Scholar] [CrossRef]
  34. Morse, S.M.; Norville, H.S. Comparison of methods to determine load sharing of insulating glass units for environmental loads. Glass Struct. Eng. 2021, 1, 315–329. [Google Scholar] [CrossRef]
  35. Galuppi, L.; Zacchei, E. Analytical solutions for plates connected by edge beams, under various loading conditions: An application for insulating glass units. Thin-Walled Struct. 2023, 188, 110877. [Google Scholar] [CrossRef]
  36. Respondek, Z. Gas Interaction in Insulating Glass Units in the Case of Elastic Support of Component Glass Panes. Constr. Optim. Energy Potential 2022, 11, 93–101. [Google Scholar] [CrossRef]
  37. Timoshenko, S.; Woinowsky-Krieger, S. Theory of Plates and Shells; McGraw-Hill Book Company: New York, NY, USA; Toronto, ON, Canada; London, UK, 1959. [Google Scholar]
  38. EN 572-1:2012; Glass in Buildings—Basic Soda Lime Silicate Glass Products Part 1: Definitions and General Physical and Mechanical Properties. CEN—European Committee for Standardization: Brussels, Belgium, 2012.
  39. Dynamic Pressure. Available online: https://en.wikipedia.org/wiki/Dynamic_pressure (accessed on 12 June 2024).
  40. Ensslen, F.; Becker, H.-R.; Wittwer, W. Modell für den nachgiebigen Randverbund von Mehrscheiben-Isoliergläsern. Stahlbau 2014, 83, 149–169. [Google Scholar] [CrossRef]
  41. Swisspacer Spacer Bars. Available online: https://www.swisspacer.com/en/products (accessed on 25 June 2024).
  42. Super Spacer® Premium. Available online: https://www.edgetechig.co.uk/our-products/spacer-bars/premium (accessed on 25 June 2024).
  43. Centa, U.G.; Oseli, A.; Mihelčič, M.; Kralj, A.; Žnidaršič, M.; Halilovič, M.; Perše, L.S. Long-Term Viscoelastic Behavior of Polyisobutylene Sealants before and after Thermal Stabilization. Polymers 2024, 16, 22. [Google Scholar] [CrossRef] [PubMed]
  44. Ihara, T.; Gustavsen, A.; Jelle, B.P. Sealant Aging and Its Correlation with Facade Reflectance. Constr. Build. Mater. 2014, 69, 390–402. [Google Scholar] [CrossRef]
  45. Bonded Windows Put Seal on Thermal Envelope. Available online: https://gbr.sika.com/en/construction/building-components/facades/knowledge-articles/bonded-windows-put-seal-on-thermal-envelope.html (accessed on 25 June 2024).
  46. Respondek, Z.; Chęciński, K. Design and assembly errors in the realization of post-and-beam glass facades. Constr. Optim. Energy Potential 2020, 9, 71–78. [Google Scholar] [CrossRef]
Figure 1. Structure of a typical IGU (author’s photo, IGU model from ©Swisspacer).
Figure 1. Structure of a typical IGU (author’s photo, IGU model from ©Swisspacer).
Materials 17 04669 g001
Figure 2. Double-glazed IGU—the symbols.
Figure 2. Double-glazed IGU—the symbols.
Materials 17 04669 g002
Figure 3. Diagram of a set of two glass panes elastically supported on the edges [36].
Figure 3. Diagram of a set of two glass panes elastically supported on the edges [36].
Materials 17 04669 g003
Figure 4. Diagram of the numerical model. S1–S4—support beams; R1—reference point for wmax, σx,cen, and σy,cen; R2—reference point for σx,ed; R3—reference point for σy,ed.
Figure 4. Diagram of the numerical model. S1–S4—support beams; R1—reference point for wmax, σx,cen, and σy,cen; R2—reference point for σx,ed; R3—reference point for σy,ed.
Materials 17 04669 g004
Figure 5. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various widths of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 5. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various widths of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g005aMaterials 17 04669 g005b
Figure 6. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various aspect ratios of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 6. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various aspect ratios of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g006
Figure 7. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various thicknesses of glass panes: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 7. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various thicknesses of glass panes: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g007aMaterials 17 04669 g007b
Figure 8. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various thicknesses of gas cavity: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 8. The influence of rotational stiffness on static quantities in IGUs loaded with a 5 kPa atmospheric pressure decrease for various thicknesses of gas cavity: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g008
Figure 9. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various widths of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 9. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various widths of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g009
Figure 10. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various aspect ratios of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 10. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various aspect ratios of the unit: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g010
Figure 11. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various thicknesses of glass panes: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 11. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various thicknesses of glass panes: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g011
Figure 12. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various thicknesses of gas cavity: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Figure 12. The influence of rotational stiffness on static quantities in IGUs loaded with a wind pressure of 0.3 kN/m2 for various thicknesses of gas cavity: (a) load reduction; (b) maximum deflection; (c) maximum stress—the solid line marks σx,cen, the dotted line marks σx,ed.
Materials 17 04669 g012
Table 1. Support conditions.
Table 1. Support conditions.
SupportRotational Degree of FreedomTranslational Degree of Freedom
Direction xDirection yDirection z
S1 and S2elastic, C (0, ∞)fixedfixed
S3 and S4fixedelastic, C (0, ∞)fixed
Table 2. Model validation. List of coefficients.
Table 2. Model validation. List of coefficients.
C [kNm/(m·rad)]Rcαwkx,cenky,cenkx,edky,ed
00.0077240.0784200.0425790.0000000.000000
0.220.362170.0071930.0743080.040168−0.006987−0.005086
14.072430.0058100.0636010.033856−0.025340−0.018403
22.036220.0048840.0564470.029599−0.037790−0.027422
3.51.163550.0041390.0507110.026148−0.047942−0.034811
50.814490.0037190.0474770.024182−0.053755−0.039096
100.407240.0030820.0425900.021173−0.062691−0.045886
1000.040720.0023010.0366140.017414−0.073807−0.055343
00.0021970.0358170.016907−0.075267−0.056765
Table 3. Model validation. Comparison of analytical and numerical calculations.
Table 3. Model validation. Comparison of analytical and numerical calculations.
C [kNm/(m·rad)]wmax [mm]σx,cen [MPa]σy,cen [MPa]σx,ed [MPa]σy,ed [MPa]
ANANUMdiff.ANANUMdiff.ANANUMdiff.ANANUMdiff.ANANUMdiff.
01.2871.2870.00%5.2935.290−0.06%2.8742.874−0.01%0.0000.022-0.0000.032-
0.21.1991.1990.01%5.0165.013−0.06%2.7112.711−0.01%−0.472−0.4934.60%−0.343−0.3769.52%
10.9680.9680.01%4.2934.290−0.06%2.2852.285−0.01%−1.710−1.7291.10%−1.242−1.2692.19%
20.8140.8140.02%3.8103.808−0.07%1.9981.998−0.01%−2.551−2.5690.69%−1.851−1.8751.28%
3.50.6900.6900.02%3.4233.421−0.07%1.7651.765−0.01%−3.236−3.2540.54%−2.350−2.3710.89%
50.6200.6200.03%3.2053.202−0.07%1.6321.6320.00%−3.628−3.6460.49%−2.639−2.6590.75%
100.5140.5140.06%2.8752.873−0.07%1.4291.4290.01%−4.232−4.2520.48%−3.097−3.1160.59%
1000.3830.3840.08%2.4712.470−0.07%1.1751.1780.24%−4.982−5.0140.64%−3.736−3.7570.58%
→∞0.36600.36640.11%2.4182.416−0.07%1.1411.1450.32%−5.081−5.1160.70%−3.832−3.8570.65%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Respondek, Z. The Influence of Elastic Support of Component Glass Panes on Deflection and Stress in Insulating Glass Units—Analytical Model. Materials 2024, 17, 4669. https://doi.org/10.3390/ma17184669

AMA Style

Respondek Z. The Influence of Elastic Support of Component Glass Panes on Deflection and Stress in Insulating Glass Units—Analytical Model. Materials. 2024; 17(18):4669. https://doi.org/10.3390/ma17184669

Chicago/Turabian Style

Respondek, Zbigniew. 2024. "The Influence of Elastic Support of Component Glass Panes on Deflection and Stress in Insulating Glass Units—Analytical Model" Materials 17, no. 18: 4669. https://doi.org/10.3390/ma17184669

APA Style

Respondek, Z. (2024). The Influence of Elastic Support of Component Glass Panes on Deflection and Stress in Insulating Glass Units—Analytical Model. Materials, 17(18), 4669. https://doi.org/10.3390/ma17184669

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop