Primary Seal Deformation in Multipane Glazing Units

Featured Application: The presented topic is related to the development of a new six-pane facade glazing unit Q-Air that has been recently successfully put on the market (http: // www.reﬂex.si / en / q-air). During R&D, the herein presented methodology for the durability assessment has been developed based on theoretical mechanics principles and used as a tool in a design phase for di ﬀ erent types and sizes of glazing units. Abstract: Quadruple glazing has become a high-end standard in the ﬁeld of sealed insulating glass units. With more than three glass panes a set of internal heat-related technical problems emerges. Durability, being the most pressing problem, requires careful management of the primary sealant deformation, especially under summer environmental conditions. Namely, the role of the primary seal in insulating glass units is to protect their humidity-sensitive, low-emissivity coatings from moisture and the quality of such protection is a key factor in durability. In this paper, a new methodology is proposed for a feasibility assessment of a new multipane insulating glass unit design, where proper design enables avoidance of excessive strain on the primary seal. It focuses on the calculation of the primary sealant strains, which are nonhomogeneous and multiaxial according to di ﬀ erent loading conditions. This approach leads to analytical expressions that enable convenient identiﬁcation of the critical location on the primary sealant. Finally, feasibility is assessed with the proposed methodology for the newly developed highly insulative six-pane facade unit by means of a comparison of the calculated strain state with the anticipated allowable strains, based on technical practice.


Introduction
In 1865, Stetson patented the first double glass unit (DGU) as an improvement of window glass that effectively reduces heat loss and improves sound protection [1]. The first insulating glass unit (IGU) was single-sealed and was soon improved by a desiccant in order to protect the interior side of the glass panes from condensation. In the late 1970s, edge seal design was improved in order to increase durability. Double-sealed IGUs were introduced, which have a primary seal between the spacers and glass, and a secondary seal above the spacers. The primary sealant is usually polyisobutylene (PIB) and, due to its thermoplastic mechanical behaviour, it can be described as a non-structural component. Its main role is to protect the insulating chambers against water vapour permeation. The secondary sealants, i.e., polysulphides, polyurethanes, or silicones [2,3], represent a mechanical restraint that

Issues and Limitations of IGU Durability
The major limitation regarding the service life of an IGU is vapour condensation fogging in the insulating chambers, i.e., on the internal glass surfaces. The durability of units is therefore related to the water vapour permeability of the edge seal, which was experimentally examined in [5], where time to fogging was predicted based on measurements of the dew point temperature of the cavity gas fill in IGUs under several ambient conditions. The moisture content in the insulating chambers was experimentally examined in [6] and can be assessed using standards for IGU durability testing [7]. A novel experimental method for diffusion gas loss in IGUs was presented in [8].
Multipane glazing units are attractive due to high energy efficiency of those buildings where they are built-in [9]. However, in multipane glazing arrangements (i.e., in glazing with more than three glass panes) a further restriction stems from the maximum permitted operating temperature of the PIB primary sealant. In two-pane and three-pane arrangements an 80 • C temperature limit [10] is difficult to exceed; however, exceeding the maximum permitted PIB temperature is easy in multipane arrangements [10].
The water vapour permeability of the dual-edge seal is almost exclusively determined by the permeability of the primary sealant because the permeability of the secondary sealant is significantly lower [11]. However, the diffusion process, as a main component of the gas permeability in the primary sealant, is exponentially dependent on temperature [12], hence resistance to atmospheric loads is the lowest during the summer. Not only a change in the material properties, but also the deformation caused by climatic loads plays a crucial role in the longevity of any IGU. Increased temperature causes a pressure increase in the insulating chambers. The pressure deflects the glass panes and the consequent change in the pressure and chamber volume is dictated by the ideal gas law. Thin glass panes enable a higher chamber volume change due to lower stiffness, which compensates for the pressure increase and lowers the pressure load on the edge seal. However, high edge seal displacements are not acceptable from a durability point of view and also due to the potential risk of glass failure [13]. Therefore, glass panes normally feature significant stiffness, which results in a pressure increase in the insulating chambers that further stretches the edge seal. Such stretching has a direct influence on moisture permeation because the deformed effective diffusion cross-section of the primary seal increases and the diffusion length decreases.
A double-seal system between two glass panes is made of two layers perpendicular to the glass: in the first one there is a spacer and a primary sealant, and the second layer enveloping the first layer represents a secondary seal. From a mechanical point of view, there are two parallel elements subjected to the same displacement, but the strains in the two sealants differ significantly. While the strain in the secondary seal is nearly homogeneous, the primary seal has to compensate for the vast majority of the displacement due to the underlying spacer stiffness. The described phenomenon is experimentally well documented together with the development of a new type of spacer in [14] and confirmed by measurements of the gas concentration decrease in the IGU's insulating chamber [15]. Therefore, in a typical IGU, the strain on the primary sealant is significantly higher than the corresponding strain on the secondary sealant.
Furthermore, the material properties of both sealants also differ significantly. The primary sealant is a thermoplastic-like liquid and is not employed as a structural component. It stretches and compresses as it is forced by the surrounding structural components. As long as the elastic deformation remains small, cyclic straining is mainly reversible. Under cyclic loading, a pumping action is observed [16] and the excessive opening between the spacer and the glass can overstretch the primary sealant.
The role of the secondary sealant, which is a structural elastomeric material, is to restrain the extension of the edge seal because the extension of the secondary seal determines how much the primary Appl. Sci. 2020, 10,1390 3 of 20 seal is permitted to deform under the expected load. As a matter of concept, once the load is reduced again, elastic recovery of the secondary sealant forces the primary sealant into its original position.
The elastic recovery of the secondary sealant is temperature-dependent and if it is not able to recover, the primary seal remains permanently deformed until negative pressure compresses the primary seal [11]. The described phenomenon increases the effective cross-section for the diffusion of moisture. Both seals are also water-sensitive and must be protected from both rain and condensed water. Any secondary seal swelling directly affects the primary seal stretching and must be avoided [17].

State of the Art
It has been recognised that IGU durability is mainly dependent on the primary seal [18], hence this aspect has been extensively investigated [19]. A thorough state-of-the-art review of a large number of primary seal systems and comparison of the performance thereof is presented by Bergh et al. [2]. Song et al. [20,21] addressed the influence of spacer material on the inside surface condensation. In [22] Wolf describes the impact of the spacer rigidity and thermal expansion coefficient (CTE) on the PIB seal deformation under various loadings. It was found that a stainless-steel spacer had less effect on the change in the PIB cross-section due to the lower CTE, while an aluminium spacer had the largest effect due to the high CTE. In order to homogenise the strain in the primary sealant, a new warm edge system, i.e., Ködispace, was developed [14]. In such a double glass unit (DGU) there is a primary seal only between two glass panes, while in multipane glazing there are more interacting systems, which significantly complicates IGU design.
Recently, multipane IGUs have appeared on the market; Grynning et al. analysed and compared their performance in [23]. The authors exposed pressing issues that are specific to multipane IGUs. Namely, in multipane glazing the intermediate chambers are highly insulated by the other chambers, therefore high cavity temperatures due to solar energy absorption can pose a problem for the robustness of such units. There are several layers of PIB and their mechanical behaviour can differ significantly from one layer to another and also along the circumference. The design requirements for such a complex structure were partially disclosed in [10].
A numerical approach is one of the possibilities to gain insight into the mechanical behaviour of an IGU. Hagl [24] performed compressive tests of edge seal spacers and built a numerical model for the sake of comparison. Stewart et al. [25] used the finite element method (FEM) to analyse the thermal displacements of the spacer frame and glass panes. However, to obtain the PIB response, the entire IGU should be exactly modelled, which requires very dense mesh and expensive analyses.

Study Aims
A determination of the strain field of all PIB layers seems an almost impossible task for several reasons. The most obvious is the PIB mechanical and rheological nonlinearity, which results in loading history-dependent behaviour. Even if a specific cyclic loading over the lifetime of the unit were somehow known, there are too many geometrical and material uncertainties that make a calculation of the exact PIB shape over several years highly unreliable. In order to obtain some guidelines for multipane IGU design, a comparative approach to assessing PIB longevity is undertaken herein.
From the IGU durability point of view, the longevity of the primary sealant is one of the major concerns. Based on the assumption that primary sealant longevity is closely related to primary sealant deformation, this paper deals with primary seal strain under different loading conditions. There is no intention herein to derive the primary sealant strain state after years of loading, but to develop an alternative methodology that provides a feasibility assessment for any new IGU design.
The methodology is based on the idea that primary sealant strains can be classified according to type of strain (normal, shear) and the intensity thereof. The type of strain can be obtained by analysing each loading case, i.e., the bulging of the IGU due to temperature increase in the chambers, temperature expansion of IGU components, dead weight, wind, etc. Each deformation mode of the seal, presented in Section 2, and the combined effect of their influences are analysed, as presented in Section 3. The most important theoretical contribution of the paper is given in Section 4, where analytical expressions for primary seal strains are derived for the most significant deformation mode. Once the strains are classified, the intensity needs to be assessed. Since there exist no allowable values in this regard, they have to be determined from technical practice. Here we rely on the fact that there already exist older IGUs that have been on the market for a long time and their durability has been proven through real application in civil engineering. Analysing such older IGUs by the methodology presented in Section 5.1 gives allowable strains for the new IGU, which is analysed using the same methodology in Section 5.2. Here it is assumed that if the primary sealant withstood a certain level of strain of a specific type in the older IGU, then the primary sealant in the new IGU is capable of withstanding the same strain level. The presented methodology has been implemented in the development of a new six-pane glass unit Q-Air, which has been successfully put on the market, as briefly presented in Section 5.3.
The main aim of this work is herein presented derivation of analytical strain field calculation, which represents a shortcut in order to avoid time-consuming FEM analyses of primary sealant behaviour, thus final analytical equations serve as a useful tool in a design phase for fast calculations and feasibility assessment of each particular panel configuration. A feasibility of assumptions used in derivation is confirmed by comparison of analytical and FEM results.

Methodology
Since there are no general criteria regarding the sealant system strain-stress state for determining the serviceable IGU lifetime, it can be deduced from the abovementioned literature where PIB strain is the limiting factor. The strain state is not equal in all PIB layers and varies along the circumference of the IGU as well. It is reasonable to assume that the location of maximal strain with proportional loading from the initial state represents a critical point also for cyclic loading and thus is the decisive factor in high-performance IGU lifespan. Figure 1 shows a schematic representation of different strain modes that occur under loading conditions in a single-chamber double-glass unit (DGU). The climatic conditions, including solar radiation, influence the temperature of the glass panes and gas within the chamber. The glass panes' dimensions vary vs. the spacer due to temperature change, and the elongation differences between the panes and spacer cause shear strains, as presented by modes A and B of the elementary PIB volume. The heat-related pressure increase in the chamber generates the extension of the PIB, shown as mode D. Furthermore, the edge rotates due to the consequent glass pane bulging, which causes shear and normal strain on the inner and outer sides of the PIB, shown as mode C.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 20 Section 3. The most important theoretical contribution of the paper is given in Section 4, where analytical expressions for primary seal strains are derived for the most significant deformation mode. Once the strains are classified, the intensity needs to be assessed. Since there exist no allowable values in this regard, they have to be determined from technical practice. Here we rely on the fact that there already exist older IGUs that have been on the market for a long time and their durability has been proven through real application in civil engineering. Analysing such older IGUs by the methodology presented in Section 5.1 gives allowable strains for the new IGU, which is analysed using the same methodology in Section 5.2. Here it is assumed that if the primary sealant withstood a certain level of strain of a specific type in the older IGU, then the primary sealant in the new IGU is capable of withstanding the same strain level. The presented methodology has been implemented in the development of a new six-pane glass unit Q-Air, which has been successfully put on the market, as briefly presented in Section 5.3.
The main aim of this work is herein presented derivation of analytical strain field calculation, which represents a shortcut in order to avoid time-consuming FEM analyses of primary sealant behaviour, thus final analytical equations serve as a useful tool in a design phase for fast calculations and feasibility assessment of each particular panel configuration. A feasibility of assumptions used in derivation is confirmed by comparison of analytical and FEM results.

Methodology
Since there are no general criteria regarding the sealant system strain-stress state for determining the serviceable IGU lifetime, it can be deduced from the abovementioned literature where PIB strain is the limiting factor. The strain state is not equal in all PIB layers and varies along the circumference of the IGU as well. It is reasonable to assume that the location of maximal strain with proportional loading from the initial state represents a critical point also for cyclic loading and thus is the decisive factor in high-performance IGU lifespan. Figure 1 shows a schematic representation of different strain modes that occur under loading conditions in a single-chamber double-glass unit (DGU). The climatic conditions, including solar radiation, influence the temperature of the glass panes and gas within the chamber. The glass panes' dimensions vary vs. the spacer due to temperature change, and the elongation differences between the panes and spacer cause shear strains, as presented by modes A and B of the elementary PIB volume. The heat-related pressure increase in the chamber generates the extension of the PIB, shown as mode D. Furthermore, the edge rotates due to the consequent glass pane bulging, which causes shear and normal strain on the inner and outer sides of the PIB, shown as mode C.  To summarize, in the PIB there is a complex strain state with the non-zero strain components being γ xz , γ xy , ε xx , and ε yy . It should be noted that strains due to other loading types, such as self-weight or wind-load, are of the same type and can be algebraically added to these strains. In order to estimate the critical strain state, the following approach is proposed: First, we determine the assumed quasi-steady-state temperature field of the analysed IGU structure. According to the prescribed climatic conditions (including solar radiation-induced internal heating), individual pane temperatures are determined. With all glass-pane temperatures known, the temperatures of other structural elements may be determined using linear interpolation, which is an approximation of the alternative 3D heat transfer analysis. The temperature field serves as the input data for the subsequent thermo-mechanical analysis.

2.
A 3D finite element method (FEM) thermo-mechanical IGU model is built, where the mechanical influence of the PIB on the global mechanical response of the structure is neglected. Although the same assumption was also applied by Stewart et al. in [25], PIB notably contributes to the rigidity of the seal system when cold. However, due to the exponential growth of vapour permeation with increasing temperature, the primary longevity-related issue is the summer, when the PIB is also soft and its mechanical rigidity really becomes negligible. Additional mechanical loads, e.g., wind load and self-weight, can be optionally added. The analyses give a reliable strain-stress state for all structural components, except for the omitted PIB. However, the results include the displacements of the glass panes and spacers to which the PIB is attached. These displacements are the input data for the identification of the type and severity of the PIB deformation modes, as presented in Figure 1.

3.
The extracted, numerically determined, displacements represent the boundary conditions for the shear and normal strain calculations of the PIB for all deformation modes separately.

4.
In accordance with the superposition principle, the analytically calculated strains for all deformation modes are summed into a strain tensor for each point of the PIB along the unit's circumference and for all PIB layers. Lastly, an equivalent strain calculation serves for the final comparative evaluation.
Each step will be described in detail in the next section, where the approach is applied to a six-pane Q-Air glazing unit [10] with a specific large-unit design, 10-18-4-18-4-18-4-18-4-18-66.4, cf. Figure 2. This consists of 10 PIB layers designed to be 0.3 mm thick. Their precise inclusion in the full 3D FEM model would be too taxing due to the extremely fine mesh required to adequately model the strain field [26].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 20 To summarize, in the PIB there is a complex strain state with the non-zero strain components being , , , and . It should be noted that strains due to other loading types, such as self-weight or wind-load, are of the same type and can be algebraically added to these strains. In order to estimate the critical strain state, the following approach is proposed: 1. First, we determine the assumed quasi-steady-state temperature field of the analysed IGU structure. According to the prescribed climatic conditions (including solar radiation-induced internal heating), individual pane temperatures are determined. With all glass-pane temperatures known, the temperatures of other structural elements may be determined using linear interpolation, which is an approximation of the alternative 3D heat transfer analysis. The temperature field serves as the input data for the subsequent thermo-mechanical analysis. 2. A 3D finite element method (FEM) thermo-mechanical IGU model is built, where the mechanical influence of the PIB on the global mechanical response of the structure is neglected. Although the same assumption was also applied by Stewart et al. in [25], PIB notably contributes to the rigidity of the seal system when cold. However, due to the exponential growth of vapour permeation with increasing temperature, the primary longevity-related issue is the summer, when the PIB is also soft and its mechanical rigidity really becomes negligible. Additional mechanical loads, e.g., wind load and self-weight, can be optionally added. The analyses give a reliable strain-stress state for all structural components, except for the omitted PIB. However, the results include the displacements of the glass panes and spacers to which the PIB is attached. These displacements are the input data for the identification of the type and severity of the PIB deformation modes, as presented in Figure 1. 3. The extracted, numerically determined, displacements represent the boundary conditions for the shear and normal strain calculations of the PIB for all deformation modes separately. 4. In accordance with the superposition principle, the analytically calculated strains for all deformation modes are summed into a strain tensor for each point of the PIB along the unit's circumference and for all PIB layers. Lastly, an equivalent strain calculation serves for the final comparative evaluation.
Each step will be described in detail in the next section, where the approach is applied to a six-pane Q-Air glazing unit [10] with a specific large-unit design, 10-18-4-18-4-18-4-18-4-18-66.4, cf. Figure 2. This consists of 10 PIB layers designed to be 0.3 mm thick. Their precise inclusion in the full 3D FEM model would be too taxing due to the extremely fine mesh required to adequately model the strain field [26].

Individual Pane Temperature Calculations
Pane temperatures for the six-pane IGU are calculated using LBNL (Lawrence Berkeley National Laboratory) Window 7.4 software [27]. Glass pane temperatures for 90% argon-filled sealed glass blocks with 18 mm spacers are calculated under the following climatic conditions: 40 • C exterior temperature, stagnant air, 24 • C interior temperature, and 783 W/m 2 solar radiation. The method was experimentally verified in [10].
The temperatures of the other components were estimated using linear interpolation, as shown in Figure 3. The gas temperature in each individual cavity is set to the average of both neighbouring glass panes. The temperature field of the unit is used as the input data for the subsequent mechanical analysis.

Individual Pane Temperature Calculations
Pane temperatures for the six-pane IGU are calculated using LBNL (Lawrence Berkeley National Laboratory) Window 7.4 software [27]. Glass pane temperatures for 90% argon-filled sealed glass blocks with 18 mm spacers are calculated under the following climatic conditions: 40 °C exterior temperature, stagnant air, 24 °C interior temperature, and 783 W/m 2 solar radiation. The method was experimentally verified in [10].
The temperatures of the other components were estimated using linear interpolation, as shown in Figure 3. The gas temperature in each individual cavity is set to the average of both neighbouring glass panes. The temperature field of the unit is used as the input data for the subsequent mechanical analysis.

Figure 3.
Temperature distribution in the structural elements. The recommended upper limit is the hybrid spacer and PIB primary sealant useful upper temperature limit [10].

IGU Mechanical Analysis
Thermal loading is the main cause of long-term strains in IGUs with more than three glass panes. The temperature increase in the chambers leads to increased cavity pressure, which affects the gas volume change in the insulation core. The temperature, pressure, and volume in each insulation chamber follow the ideal gas law so as to achieve a static equilibrium. The modelled gas is isothermal, but pressure and volume at equilibrium are affected by the deflection of the glass panes, and the numerical model calculates the equilibrium state. The shape of the unit defines the volume of the given chamber and also recursively affects the pressure within the chamber. It should also be noted that the argon-filled chambers in the investigated six-pane configuration are pneumatically interconnected by 3 mm holes in the intermediate glass panes.
In the FEM model, 20-node quadratic stress/displacement continuum elements are used for the glass panes and sealants, and 8-node doubly-curved thick stress/displacement shell elements are used for the spacers. Additionally, a hydrostatic fluid volume element interconnects the temperature, pressure, and volume of each insulation chamber according to the ideal gas law. The boundary conditions consider an installed IGU as part of a curtain wall facade. A half-model is analysed due to symmetry. The model considers thermal expansion and loading due to cavity pressure. Other load types, e.g., external pressure representing wind loads and glass self-weight, are also considered in the model. The key FEM results needed for further detailed sealant analysis are the displacements of the structural elements attached to the PIB.

Primary Sealant Mechanical Behaviour
The global FEM model described above gives a general overview of the structural integrity of the product, i.e., the deflection of all panes and their stresses and the strains on the spacers and secondary sealant. It is of major importance for further analysis that this model gives the displacements of all points of the glass panes and spacers that are in contact with the primary

IGU Mechanical Analysis
Thermal loading is the main cause of long-term strains in IGUs with more than three glass panes. The temperature increase in the chambers leads to increased cavity pressure, which affects the gas volume change in the insulation core. The temperature, pressure, and volume in each insulation chamber follow the ideal gas law so as to achieve a static equilibrium. The modelled gas is isothermal, but pressure and volume at equilibrium are affected by the deflection of the glass panes, and the numerical model calculates the equilibrium state. The shape of the unit defines the volume of the given chamber and also recursively affects the pressure within the chamber. It should also be noted that the argon-filled chambers in the investigated six-pane configuration are pneumatically interconnected by 3 mm holes in the intermediate glass panes.
In the FEM model, 20-node quadratic stress/displacement continuum elements are used for the glass panes and sealants, and 8-node doubly-curved thick stress/displacement shell elements are used for the spacers. Additionally, a hydrostatic fluid volume element interconnects the temperature, pressure, and volume of each insulation chamber according to the ideal gas law. The boundary conditions consider an installed IGU as part of a curtain wall facade. A half-model is analysed due to symmetry. The model considers thermal expansion and loading due to cavity pressure. Other load types, e.g., external pressure representing wind loads and glass self-weight, are also considered in the model. The key FEM results needed for further detailed sealant analysis are the displacements of the structural elements attached to the PIB.

Primary Sealant Mechanical Behaviour
The global FEM model described above gives a general overview of the structural integrity of the product, i.e., the deflection of all panes and their stresses and the strains on the spacers and secondary sealant. It is of major importance for further analysis that this model gives the displacements of all  The primary sealant could be analysed as a FEM submodel by prescribing displacements as boundary conditions, but it turns out that such an approach is computationally expensive and has several convergence problems due to the PIB material properties where PIB could be treated as an incompressible solid of very low modulus. For practical reasons, an alternative approach that considers the combined effects of separate deformation modes, with each of them representing a separate task based on the type of loading, was chosen.
Mode A: Shear strains occur mainly due to the thermal expansion difference between the glass pane and spacer in the longitudinal (z) direction. In cases where a unit is double-symmetric, as presented in Figure 1, is zero on the symmetry planes and is the highest in the corners. In general, however, self-weight affects shear strains too. We assume that is location-independent at the PIB cross-section and is dependent on the location on the unit's circumference only; therefore, ( , , ) = ( ) can be calculated from: here h is the thickness of the PIB layer. Mode B: Shear strains occur due to the thermal expansion difference between the pane and the spacer in the transverse (y) direction and due to the increased cavity pressure in the chamber that loads the spacer from the inner side. Simple shear is assumed for convenience, which means that is the same over the entire PIB cross-section, hence ( , , ) = ( ). In reality, a spacer made of thin-walled material of low rigidity deforms locally and the spacer wall is no longer precisely parallel with the glass pane; therefore, the average displacements of both contact surfaces with the PIB are needed for the following calculation of the shear strain: Mode C: Glass panes bulge under increased cavity pressure; therefore, the glass pane is no longer parallel with the surface of the spacer. The inclination of the PIB layer can be described by the shear strain and normal strain . The maximal values of both strains on the PIB cross-section can be approximated with: , Mode D: Increased pressure in the chambers moves the glass panes apart and the role of the secondary sealant is to capture this effect. However, the primary sealant follows the deformations of the secondary sealant, causing the thinning of the primary sealant layer. Obviously, when this deformation mode is established, the PIB layer thins, meaning that the effective cross-section for moisture diffusion increases for the duration of the highest summer climatic loadings and is thus of the highest interest. Moreover, the instantaneously deformed state would not be that critical if the material behaviour of the primary sealant were not viscoplastic. Low elastic recovery causes a The primary sealant could be analysed as a FEM submodel by prescribing displacements as boundary conditions, but it turns out that such an approach is computationally expensive and has several convergence problems due to the PIB material properties where PIB could be treated as an incompressible solid of very low modulus. For practical reasons, an alternative approach that considers the combined effects of separate deformation modes, with each of them representing a separate task based on the type of loading, was chosen.
Mode A: Shear strains γ xz occur mainly due to the thermal expansion difference between the glass pane and spacer in the longitudinal (z) direction. In cases where a unit is double-symmetric, as presented in Figure 1, γ xz is zero on the symmetry planes and is the highest in the corners. In general, however, self-weight affects shear strains too. We assume that γ xz is location-independent at the PIB cross-section and is dependent on the location on the unit's circumference only; therefore, γ xz (x, y, z) = γ xz (z) can be calculated from: here h is the thickness of the PIB layer. Mode B: Shear strains γ xy occur due to the thermal expansion difference between the pane and the spacer in the transverse (y) direction and due to the increased cavity pressure in the chamber that loads the spacer from the inner side. Simple shear is assumed for convenience, which means that γ xy is the same over the entire PIB cross-section, hence γ xy (x, y, z) = γ xy (z). In reality, a spacer made of thin-walled material of low rigidity deforms locally and the spacer wall is no longer precisely parallel with the glass pane; therefore, the average displacements of both contact surfaces with the PIB are needed for the following calculation of the shear strain: Mode C: Glass panes bulge under increased cavity pressure; therefore, the glass pane is no longer parallel with the surface of the spacer. The inclination of the PIB layer can be described by the shear strain γ xy and normal strain ε xx . The maximal values of both strains on the PIB cross-section can be approximated with: Mode D: Increased pressure in the chambers moves the glass panes apart and the role of the secondary sealant is to capture this effect. However, the primary sealant follows the deformations of the secondary sealant, causing the thinning of the primary sealant layer. Obviously, when this deformation mode is established, the PIB layer thins, meaning that the effective cross-section for moisture diffusion increases for the duration of the highest summer climatic loadings and is thus of the highest interest. Moreover, the instantaneously deformed state would not be that critical if the material behaviour of Appl. Sci. 2020, 10, 1390 8 of 20 the primary sealant were not viscoplastic. Low elastic recovery causes a pumping effect during cyclic loading and the loss of sealant integrity; therefore, it is of primary importance that the strain remains as low as possible.
The primary sealant deformation mode D cross-section consists of three strains, ε xx , ε yy , and γ xy . The multiaxial strain state is non-homogeneous and its determination is nontrivial. As will be shown in the results of a case study, mode D is the most influential among all modes, thus the derivation of analytical expressions for all strain components for mode D is crucial and is presented in the next section. The final expression shows that maximal values for the mode D strain field can be expressed as: Total and maximal strain: Once all components of the analysed strain modes are calculated, the same kind of components are combined to get the strain tensor for the multiaxial strain state.
Strains at different locations are compared via a magnitude of maximal shear strain, calculated as:

Derivation of Analytical Expressions for Deformation Mode D
Deformation mode D considers opening or closing the PIB layer between the glass pane and spacer ( Figure 5). Relative displacements in the normal direction cause the thinning or thickening of the PIB layer. As shown in Figure 1, the deformation pattern of the PIB cross-section includes normal strains ε xx (x, y, z), ε yy (x, y, z), and shear strains γ xy (x, y, z). Determination of all these strain field distributions will be based on the following analytical approach.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 20 pumping effect during cyclic loading and the loss of sealant integrity; therefore, it is of primary importance that the strain remains as low as possible. The primary sealant deformation mode D cross-section consists of three strains, , , and . The multiaxial strain state is non-homogeneous and its determination is nontrivial. As will be shown in the results of a case study, mode D is the most influential among all modes, thus the derivation of analytical expressions for all strain components for mode D is crucial and is presented in the next section. The final expression shows that maximal values for the mode D strain field can be expressed as: Total and maximal strain: Once all components of the analysed strain modes are calculated, the same kind of components are combined to get the strain tensor for the multiaxial strain state.
Strains at different locations are compared via a magnitude of maximal shear strain, calculated as:

Derivation of Analytical Expressions for Deformation Mode D
Deformation mode D considers opening or closing the PIB layer between the glass pane and spacer ( Figure 5). Relative displacements in the normal direction cause the thinning or thickening of the PIB layer. As shown in Figure 1, the deformation pattern of the PIB cross-section includes normal strains ( , , ), ( , , ), and shear strains ( , , ). Determination of all these strain field distributions will be based on the following analytical approach. The first question defines the shape of the PIB when the distance between both surfaces (glass and spacer) increases/decreases. If the PIB is considered to be a liquid medium, the answer is given The first question defines the shape of the PIB when the distance between both surfaces (glass and spacer) increases/decreases. If the PIB is considered to be a liquid medium, the answer is given Appl. Sci. 2020, 10, 1390 9 of 20 by Honschoten et al. [28], who examined the capillary liquid bridge between solid surfaces. Liquid takes a form that minimizes its outer surface area due to surface tension, q sur f , which holds a curved differential patch of a doubly-curved liquid surface in equilibrium. Both principal radii of the curvature are related to the Young-Laplace equation, which also considers the pressure difference between the liquid pressure just below the surface and the gas pressure just above the surface, ∆p.
In the absence of gravitational and electrostatic forces, the pressure gradient inside a fluid is zero. For plane strain conditions, one of the principal radii is infinite, meaning that a cross-sectional curve is a circular arc ( Figure 6).
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 20 by Honschoten et al. [28], who examined the capillary liquid bridge between solid surfaces. Liquid takes a form that minimizes its outer surface area due to surface tension, , which holds a curved differential patch of a doubly-curved liquid surface in equilibrium. Both principal radii of the curvature are related to the Young-Laplace equation, which also considers the pressure difference between the liquid pressure just below the surface and the gas pressure just above the surface, ∆ .
In the absence of gravitational and electrostatic forces, the pressure gradient inside a fluid is zero. For plane strain conditions, one of the principal radii is infinite, meaning that a cross-sectional curve is a circular arc ( Figure 6). For the sake of brevity, let us assume an initially rectangular domain of PIB with dimensions 2 × 2 , which is stretched by glass panes in the -direction. Due to symmetry, only a quarter of the domain is presented in Figure 7. All points of the PIB at = are attached on a glass pane, hence their displacements follow the pane, i.e., they move horizontally by ∆ . The free surface at = forms a circular arc ( , ) in a deformed state. We also assume that all points at the same height also form a circular arc ( , ) with radius ( ), which follows the equation: where ( ) determines the location of the circle centre for each arc along y and is determined from the condition that each circle must lead through the point ( + Δ , ) = . Each circular arc is then defined as where (location in a deformed state) is assumed to be linearly scaled from the initial domain as: For the sake of brevity, let us assume an initially rectangular domain of PIB with dimensions 2a × 2b, which is stretched by glass panes in the x-direction. Due to symmetry, only a quarter of the domain is presented in Figure 7. All points of the PIB at x = a are attached on a glass pane, hence their displacements follow the pane, i.e., they move horizontally by ∆a. The free surface at y = b forms a circular arc Y(X, b) in a deformed state. We also assume that all points at the same height y also form a circular arc Y(X, y) with radius R(y), which follows the equation: where B(y) determines the location of the circle centre for each arc along y and is determined from the condition that each circle must lead through the point Y(a + ∆a, y) = y. Each circular arc is then defined as where X (location in a deformed state) is assumed to be linearly scaled from the initial domain as: If we consider the incompressibility of the PIB, the radius ( ) is determined from the volume preservation, i.e., the initial volume (or area under plane strain conditions) 0 ( ) = is the same as the deformed volume (area) ( ), defined as: Considering 0 ( ) = ( ), it follows an implicit nonlinear equation for K: However, K cannot be extracted explicitly, but since ≪ 1, it is admissible to perform a Taylor series expansion for some functions in the above expression. Fourth-order expansions would give √1 − 2 ≈ − 3 2 and Arcsin( ) ≈ + 3 6 , which gives a linear expression for K and consequently it can be extracted as: Which finally gives the linear approximation of the radius ( ): Insertion of ( ) into Equation (11) gives the function of all circular arcs along with the height of domain ( , ), which is presented in Figure 7: If we consider the incompressibility of the PIB, the radius R(y) is determined from the volume preservation, i.e., the initial volume (or area under plane strain conditions) V 0 (y) = ay is the same as the deformed volume (area) V(y), defined as: Considering V 0 (y) = V(y), it follows an implicit nonlinear equation for K: However, K cannot be extracted explicitly, but since K 1, it is admissible to perform a Taylor series expansion for some functions in the above expression. Fourth-order expansions would give K √ 1 − K 2 ≈ K − K 3 2 and Arcsin(K) ≈ K + K 3 6 , which gives a linear expression for K and consequently it can be extracted as: Which finally gives the linear approximation of the radius R(y) : Insertion of R(y) into Equation (11) gives the function of all circular arcs along with the height of domain Y(X, y), which is presented in Figure 7: A point initially located in (x, y) moves to the displaced coordinates (X, Y), hence the vertical displacement is given by v(x, y) = Y(X(x), y) − y, which finally gives the vertical component of the displacement field: v(x, y) = (a + ∆a) Considering ∆a a, we can approximate the vertical displacement by means of a Taylor series expansion, which, after rearrangement, gives: Now strains ε yy (x, y) can be calculated as: The next task is to determine the horizontal displacement component u(x, y). For incompressible solids no volume change occurs, hence, according to small strain theory, the trace of a strain tensor is zero: The elaborated case can be considered as a plane strain problem, where ε zz (x, y) = 0, therefore the absolute values of the strains are the same for the whole domain ε xx (x, y) = −ε yy (x, y). Horizontal displacement is now calculated from the definition: Considering the symmetry boundary condition u(0, y) = 0, it follows: All displacement and strain fields are presented in Figure 8.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 11 of 20 Considering Δ ≪ , we can approximate the vertical displacement by means of a Taylor series expansion, which, after rearrangement, gives: Now strains ( , ) can be calculated as: The next task is to determine the horizontal displacement component ( , ). For incompressible solids no volume change occurs, hence, according to small strain theory, the trace of a strain tensor is zero: The elaborated case can be considered as a plane strain problem, where ( , ) = 0, therefore the absolute values of the strains are the same for the whole domain ( , ) = − ( , ). Horizontal displacement is now calculated from the definition: Considering the symmetry boundary condition (0, ) = 0, it follows: All displacement and strain fields are presented in Figure 8. If the obtained solution is checked with regard to boundary conditions, it is evident that all conditions are fulfilled, except that the shear strains on the free (upper) edge ( , ) are not zero, as they should be, since the edge is traction-free. It is reasonable to assume that the solution is a fair approximation of the far-field strains and displacements, but in the vicinity of the upper edge, the solution has to be upgraded.
For such an upgrade, it is necessary to define the region of influence that the free edge affects. To answer this issue, let us consider the deformation pattern on the free edge. Initially, the regular pattern has to be perpendicular in all free surface points in a deformed state (see Figure 9). The points that are initially at any arbitrary = remain a vertical = distance from the If the obtained solution is checked with regard to boundary conditions, it is evident that all conditions are fulfilled, except that the shear strains on the free (upper) edge Y(x, b) are not zero, as they should be, since the edge is traction-free. It is reasonable to assume that the solution is a fair approximation of the far-field strains and displacements, but in the vicinity of the upper edge, the solution has to be upgraded.
For such an upgrade, it is necessary to define the region of influence that the free edge affects. To answer this issue, let us consider the deformation pattern on the free edge. Initially, the regular pattern has to be perpendicular in all free surface points in a deformed state (see Figure 9). The points that are initially at any arbitrary x = const remain a vertical X = const distance from the surface, but tend to incline in the upper region to be perpendicular with the surface. To capture the behavior in the upper part, two functions are defined. The first one is f bound (x), which divides both regions, and the second one is f corr (y), which will help to define the correction of the horizontal displacement field in the upper part, i.e., y ≥ f bound (x). It turns out that both functions can be defined as quadratic functions, f bound (x) = a 1 + a 2 x + a 3 x 2 and f corr (y) = a 4 + a 5 y + a 6 y 2 (see Figure 9). The points that are close to the right edge cannot have a wide region of correction. Although on the surface the condition of perpendicularity must be fulfilled, below the surface the points have to reach a vertical position not far from the surface; the limit case is exactly the right edge, which is vertical by definition, hence the correction from the inclined upper surface to the vertical position has to be theoretically made at a single point, i.e., in the upper right corner.  ′ (0) = 0: Assuming the symmetry of the influential region is reasonable.  (0) = − 7 : This equation defines the range of the function on the symmetry plane, but according to the new unknown 7 , another condition is needed. The free surface far from the right edge is traction-free (and consequently has no strains) on the free surface, hence let us consider (0, ) = 0. However, the strains must take into consideration the new displacement field, which already includes the correction: The equation system is linear and all seven unknowns, 1 to 7 , can be obtained analytically. Using known functions, the upgraded displacement field can be defined as: There are six parameters for both functions, which are defined according to the following assumptions: • f corr ( f bound (x)) = 0: The function f corr will not correct the far-field solution of the horizontal displacement of points lying on the boundary.
The perpendicularity of the parabola f corr and circular arc Y(x, b) on the surface. • f corr ( f bound (x)) = 0: The smooth transition of the deformed curve, which was initially x = const between both regions. • f bound (a) = b: The points that are close to the right edge cannot have a wide region of correction. Although on the surface the condition of perpendicularity must be fulfilled, below the surface the points have to reach a vertical position not far from the surface; the limit case is exactly the right edge, which is vertical by definition, hence the correction from the inclined upper surface to the vertical position has to be theoretically made at a single point, i.e., in the upper right corner. • f corr (0) = 0: Assuming the symmetry of the influential region is reasonable.
This equation defines the range of the function f bound on the symmetry plane, but according to the new unknown a 7 , another condition is needed. The free surface far from the right edge is traction-free (and consequently has no strains) on the free surface, hence let us consider ε xx (0, b) = 0. However, the strains must take into consideration the new displacement field, which already includes the correction: The equation system is linear and all seven unknowns, a 1 to a 7 , can be obtained analytically. Using known functions, the upgraded displacement field can be defined as: ; ; Once u ALL (x, y) is known, v ALL (x, y) is derived analogically to Equation (22). Without the loss of generality, the solution can be converted into a more convenient form using non-dimensional variables ξ = x a , η = y b and ε 0 = ∆a a . The final solution of the displacement field is: ; ; ; ; where f bound (ξ) = 1 − a b 2 1 − ξ 2 . The displacement and strain fields are presented in Figure 10 for a = 4mm, b = 10mm, and ∆a = 0.1mm.   As a final result, we can also extract extreme values for all strains. Extreme values are in the rightmost and uppermost corner, i.e., in the location (ξ = η = 1): Verification: In the presented derivation, there are a moderate number of assumptions and approximations. The first check of the final solution correctness is to verify the fulfilment of the boundary conditions. The displacement field Equations (26) and (27) give u ALL (0, η) = 0, v ALL (ξ, 0) = 0 and u ALL (1, η) = aε 0 = ∆a, which confirms that all kinematic boundary conditions are exactly met.
The second verification is related to the entire displacement and strain fields. Analytical expressions will be compared with the results of a numerical analysis performed with the finite element method using the Abaqus FEM program [29].
In FEM, a domain 4 mm × 10 mm has been discretized into 4000 plane strain 8-node biquadratic hybrid elements (Abaqus notation CPE8H). The left and lower edges have been symmetrically supported, and on the right edge vertical displacements have been restrained, while the horizontal displacements are defined as ∆a =0.1 mm. The upper edge is traction free. The results are presented in Figure 10 and a comparison with the analytical displacement and strain field confirms that the derived expressions, which have a quite manageable form, are a fair approximation of the strain field calculations of the presented problem.

Reference Case
The described methodology gives the strain state in all layers of the primary sealant in a multipane IGU, but there exist no standardised allowable values to compare this strain state with. In order to obtain allowable values, we will rely on engineering practice. SIKA Facade Systems [30] are well known for their long use in civil engineering, which proves the longevity of such units. From a mechanical point of view, the behaviour of SIKA double glass units (DGUs) can be analysed by means of the same methodology as described above, and such a calculation will be called a reference case.
The calculated primary sealant strains in this reference case are conservatively treated as the maximal strains that the primary sealant is capable of withstanding, which represents the allowable strain value for the primary sealant of any other arbitrary IGU configuration. It is reasonable to assume that a new IGU configuration is feasible if its primary sealant strains do not exceed the values from the reference case.
For this reason, we chose a classic 1 m high and 0.5 m wide double glass unit 6-12-6 at 20 kPa internal overpressure for analysis using the same methodology. Figure 11 presents the maximal shear strain on the primary sealant along the edge of both layers for one-quarter of the model. The strains are equal in both layers due to the symmetrical unit configuration 6-12-6 and the presence of thermal loads only. The analysis reveals that the pressure change in mode D has by far the greatest impact on the total strain. The maximal shear strain distribution has two extreme values, the first in the middle of the longest edge due to the expansion of the glass panes, and the second in the corner due to the compression of the deformed, i.e., doubly-curved, glass panes.
shear strain on the primary sealant along the edge of both layers for one-quarter of the model. The strains are equal in both layers due to the symmetrical unit configuration 6-12-6 and the presence of thermal loads only. The analysis reveals that the pressure change in mode D has by far the greatest impact on the total strain. The maximal shear strain distribution has two extreme values, the first in the middle of the longest edge due to the expansion of the glass panes, and the second in the corner due to the compression of the deformed, i.e., doubly-curved, glass panes. The analysis of the reference DGU conservatively considers only climatic loading and no other external loads, e.g., self-weight, wind loads, etc. In this specific case, the absolute maximal value is The analysis of the reference DGU conservatively considers only climatic loading and no other external loads, e.g., self-weight, wind loads, etc. In this specific case, the absolute maximal value is reached in the middle of the longest edge. Since the longevity of the DGU is verified through engineering practice, the values from Figure 11 will be called the reference values and the maximum will be used as the allowable strain in the feasibility assessment of a newly designed multipane IGU.
In the next subsection, the multipane IGU strain will be compared with the allowable strain. If the calculated primary sealant strain in all locations of the IGU is below the allowable strain, the longevity of the analysed IGU is confirmed.

Analysis of the Multipane IGU
The same methodology as for the DGU is also used for the above-described six-pane IGU. The analysis considers climate loading, self-weight, and wind load. Wind load is modelled as a uniform pressure of 3.6 kPa, which is equivalent to a wind gust of about 250 km/h.
The six-pane IGU has 10 layers of primary sealant. The maximal shear strain distribution for all the layers is presented in Figure 12. The analysis is performed in three consecutive steps: in the first one only self-weight is considered, in the second step thermal loading is added, and finally, wind load is applied.
The same methodology as for the DGU is also used for the above-described six-pane IGU. The analysis considers climate loading, self-weight, and wind load. Wind load is modelled as a uniform pressure of 3.6 kPa, which is equivalent to a wind gust of about 250 km/h. The six-pane IGU has 10 layers of primary sealant. The maximal shear strain distribution for all the layers is presented in Figure 12. The analysis is performed in three consecutive steps: in the first one only self-weight is considered, in the second step thermal loading is added, and finally, wind load is applied.  The results reveal that shear strains caused by self-weight are classified as mode A on the vertical edges and mode B on the horizontal edges and are significantly lower in comparison to the other loadings. In the final step, the thermal loads are applied, which cause all four modes, but it is evident that mode D is by far the most dominant deformation mode due to its large contribution to the magnitude of the strain components ( Figure 13). The same observation is recorded for the wind load, which means that the transverse pressure loading caused either by chamber temperature increase or wind loads has a much more important effect on the primary sealant strains than shear The results reveal that shear strains caused by self-weight are classified as mode A on the vertical edges and mode B on the horizontal edges and are significantly lower in comparison to the other loadings. In the final step, the thermal loads are applied, which cause all four modes, but it is evident that mode D is by far the most dominant deformation mode due to its large contribution to the magnitude of the strain components ( Figure 13). The same observation is recorded for the wind load, which means that the transverse pressure loading caused either by chamber temperature increase or wind loads has a much more important effect on the primary sealant strains than shear strains due to the self-weight or glass temperature difference.

6.
7. 8. 9. 10. The results reveal that shear strains caused by self-weight are classified as mode A on the vertical edges and mode B on the horizontal edges and are significantly lower in comparison to the other loadings. In the final step, the thermal loads are applied, which cause all four modes, but it is evident that mode D is by far the most dominant deformation mode due to its large contribution to the magnitude of the strain components ( Figure 13). The same observation is recorded for the wind load, which means that the transverse pressure loading caused either by chamber temperature increase or wind loads has a much more important effect on the primary sealant strains than shear strains due to the self-weight or glass temperature difference. Figure 13. The contribution of deformation modes to the strain components in the primary sealant for different load cases (a) self-weight, (b) self-weight + temperature, (c) self-weight + temperature + wind; location on the edge: Point A (cf. Figure 12).
The most important result of this analysis is that the maximal strains are below the allowable strain. The maximal shear strains in the six-pane Q-Air glazing unit also take into consideration Figure 13. The contribution of deformation modes to the strain components in the primary sealant for different load cases (a) self-weight, (b) self-weight + temperature, (c) self-weight + temperature + wind; location on the edge: Point A (cf. Figure 12).
The most important result of this analysis is that the maximal strains are below the allowable strain. The maximal shear strains in the six-pane Q-Air glazing unit also take into consideration self-weight and wind loads, but they are still below the allowable strain, which was calculated using climatic loading of 20 kPa internal over-pressure, as described in the previous section. It can be also noticed that the strains in the middle of the edge and in the corner are in the same range and that there is no especially highly loaded region, which can become critical during the lifetime of the unit. From these results, the value of such an analysis for the design phase can be recognised. The described methodology can provide answers to the composition of a multipane IGU where strains in all primary sealant layers are sufficiently low.

Practical Aspects
Quadruple glazing is a cost-efficient solution. Quadruple glazing offers glazing with U w values as low as 0.4 W/m 2 K, which is sufficient to abandon otherwise common external, modulated sun-shading and trading for the solar control coating. Abandoning of external sun-shading is sufficient to offset for the cost of the extra glass pane whose add-on cost is at about 10€/m 2 , so payback is immediate as savings are direct on the up-front basis when compared with triple-pane glazing with external sun-shading. Cost benefits don't end here. In some cases, it is possible to abandon central heating for the building and make the building so-called "zero heating". A further advantage is marked acceptance. In Norway, for example, there are many triple glazed buildings now in their third or fourth decade-ready to be refurbished.
Four to six-pane Q-Air units were first produced in 2010 and have subsequently been field-tested since 2012. Industrial property for handling sealed chamber insulating gas pressure variations and for handling temperature effects for intermediate panes and sealants was obtained. Their first commercial application was in 2014. Six-pane Q-Air glazing units were used in a major renovation of the Wergelandsveien 7 building in Oslo (Figure 14), where the old triple-pane glazing with thermal transmittance of glass U g -value 2.2 W/(m 2 K) was replaced with the new system with a U g -value of 0.24 W/(m 2 K). A total of 550 small fixed glass, window-size, six-pane units were installed in 2015.
In 2016 a further 20 large (jumbo) size six-pane units were installed on the ground floor. The presented methodology has been used for durability assessment in a design phase of the Q-Air. abandon central heating for the building and make the building so-called "zero heating". A further advantage is marked acceptance. In Norway, for example, there are many triple glazed buildings now in their third or fourth decade-ready to be refurbished.
Four to six-pane Q-Air units were first produced in 2010 and have subsequently been field-tested since 2012. Industrial property for handling sealed chamber insulating gas pressure variations and for handling temperature effects for intermediate panes and sealants was obtained. Their first commercial application was in 2014. Six-pane Q-Air glazing units were used in a major renovation of the Wergelandsveien 7 building in Oslo (Figure 14), where the old triple-pane glazing with thermal transmittance of glass Ug-value 2.2 W/(m 2 K) was replaced with the new system with a Ug-value of 0.24 W/(m 2 K). A total of 550 small fixed glass, window-size, six-pane units were installed in 2015. In 2016 a further 20 large (jumbo) size six-pane units were installed on the ground floor. The presented methodology has been used for durability assessment in a design phase of the Q-Air. Figure 14. The Wergelandsveien 7 building in Oslo before (left) and after renovation (right) [10].

Conclusions
Primary sealant strain is of great significance to insulating glass unit (IGU) lifetime. After cyclic loading with excessive strain, the seal functionality of the insulation chambers can be prematurely compromised. Excessive thinning of the primary sealant increases moisture permeation and may provoke premature fogging of the glass panes.

Conclusions
Primary sealant strain is of great significance to insulating glass unit (IGU) lifetime. After cyclic loading with excessive strain, the seal functionality of the insulation chambers can be prematurely compromised. Excessive thinning of the primary sealant increases moisture permeation and may provoke premature fogging of the glass panes.
The presented new methodology assesses the strain state in all primary sealant layers. The strain state is based on the displacement of the surrounding structural elements, determined by a global FEM model. Primary sealant strains are then calculated using derived analytical expressions. The methodology considers the combination of four deformation modes arising from different loading scenarios. The most influential and of the highest interest is the extension of the sealant (i.e., mode D), which is also the most challenging to characterize due to its multiaxial strain state. Such behaviour cannot be described with simple engineering expressions. The main theoretical contribution of the paper is thus a derivation of analytical expressions for a strain field estimation for deformation mode D, which is the most influential among all the other modes.
The analyses reveal that an increase in pressure in the insulating chambers following a temperature increase represents by far the most significant influence on the sealant strain state; therefore, pressure variations in the insulating chambers should be a major concern in the design phase of any multipane IGU. This means that self-weight and differences in linear thermal expansion coefficients play only a minor role in multipane IGU durability.
Our method is built on the assumption of proportional loading from the initial state. However, during the lifetime of any IGU very complex loading occurs, i.e., cyclic weather conditions (summer/winter), stochastic external loads (wind), etc. Since the behaviour of incompressible primary sealant can be described as perfectly plastic, loading history plays a crucial role. It would seem impossible to determine realistic sealant strains during an unknown loading history, thus the presented methodology is based on a comparison of the maximal strain in two glass units, where the first one is an already established reference maximum from civil engineering practice and the second is the newly designed unit. If the peak strain in the investigated unit does not exceed the allowable strain, it is deemed that the unit is validated and ready to be applied in practice. The developed methodology thus represents a useful tool for IGU durability assessment.