Next Article in Journal
Numerical Strength Analysis of Laser-Welded Differential Housing and Gear Considering Residual Stress
Next Article in Special Issue
Genetic Algorithm Optimization of Beams in Terms of Maximizing Gaps between Adjacent Frequencies
Previous Article in Journal
Triple-Band and Ultra-Broadband Switchable Terahertz Meta-Material Absorbers Based on the Hybrid Structures of Vanadium Dioxide and Metallic Patterned Resonators
Previous Article in Special Issue
Nonlinear Vibration of Double-Walled Carbon Nanotubes Subjected to Mechanical Impact and Embedded on Winkler–Pasternak Foundation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis, Optimization, and Multi-Criteria Design of Vacuum Insulated Glass Composite Panels

1
Department of Structural Mechanics, Lodz University of Technology, Politechniki 6, 93-590 Lodz, Poland
2
Department of Concrete Structures, Lodz University of Technology, Politechniki 6, 93-590 Lodz, Poland
*
Author to whom correspondence should be addressed.
Materials 2023, 16(13), 4722; https://doi.org/10.3390/ma16134722
Submission received: 9 June 2023 / Revised: 26 June 2023 / Accepted: 27 June 2023 / Published: 29 June 2023
(This article belongs to the Special Issue Modelling of Microstructured Media)

Abstract

:
The subject of this study is Vacuum Insulated Glass (VIG) panels, which consist of two glass panes with an evacuated space and evenly distributed micro-support pillars between them. The deflection of panes towards the centre of the structure caused by atmospheric pressure is a mechanical problem that occurs in this type of structure. The aim of this study was to extend previous research on the optimal arrangement of support pillars in terms of eigenfrequencies and dynamics to include aesthetic aspects. Using Abaqus/CAE v2017 software, a large number of numerical models were created and subjected to a comprehensive multi-criteria analysis. Fractal analysis was employed to automatically assess the aesthetics of the proposed solutions. The study presents theoretical solutions that could be implemented in industrial production. The presented study shows that it is possible to effectively extend the criteria for optimizing the arrangement of pillars with new design criteria. Most studies focus on pillar placement, amount, or shape in terms of panes thermal or mechanical properties. Due to the increasing number of VIG panels applications in places exposed to external vibrations, other design criteria for VIG panels are also required and are provided by the following study.

1. Introduction

1.1. VIG Panels

Nowadays, energy management is an important problem in the world economy. Its price has risen sharply in recent years. It is caused by the need to move away from fossil fuels and the dwindling resources of the planet. The recent geopolitical problems are an additional factor deepening the upward trend in energy prices. Therefore, any initiative aimed at saving energy is a priority for the economy.
The subject of this work is VIG (Vacuum Insulating Glass) panels. They represent the latest technological achievement that is already implemented in many countries. The VIG technology consists of two high-strength glass panes with a vacuum between them (Figure 1). In order to counteract the atmospheric pressure and the contact of the panels, they are connected to micro-support pillars—elements usually made of steel with a diameter of 0.3–0.5 mm [1]. Moreover, the panes are sealed around the edges. Consequently, this design provides excellent thermal properties to VIG panels. Their U-value can be as low as 0.2–0.4 W/m2K [2].

1.2. Fractals in Engineering

In recent years, scientists and engineers have been looking for inspiration to create new, original structures with better properties. Among the many interesting trends, the attempts to apply structures called fractals in various fields of technology are remarkable [3,4,5]. These are mathematical structures discovered and described by Mandelbrot in the 1970s. Their representation in nature is significant [6]. They are characterized by a sophisticated, self-similar structure, the similarity scale of which is described by a parameter called fractal dimension. Over the next years, fractals will be subject of the extensive research. The individual section of mathematics, such as fractal geometry, was created. The wide possibilities of fractal applications were also analysed, both to create new technical structures and to use the so-called fractal analysis to describe the features of existing structures [7]. Fractals application in architecture and civil engineering [8] deserves special attention—particularly in the dynamics of plate structures [9]. Fractals have been seen in architecture across many eras and cultures.
The richness and variety of fractal forms naturally occurs in nature but also in the architecture of many cultures. This leads to natural questions of the aesthetic value of fractals. The proportions related to the similarity of fractal structures provides the observer an aesthetic feeling which has already been the subject of scientific research [10]. The mathematical description of fractal structures gives the engineers possibility to quantify the aesthetic criterion that can be used for generative design of building objects [11].

1.3. Motivation

Thermal properties of VIG panels [12,13], strength properties of support-pillars and their static properties [14,15] were subject of numerous studies. However, in recent years VIG technology has been analysed also in terms of its acoustic properties [16].
An interesting trend in the recent research on VIG technology are studies on its dynamic properties. In [17], the authors, using an advanced FEM model, analysed a number of plates with different geometries and compared their dynamic properties. By determining the value of the natural frequency of these structures, it is possible to determine the possibility of using VIG plates in buildings exposed to external dynamic factors, such as traffic or industrial vibrations. Another work worth mentioning is [18], in which the authors made an attempt to optimize pillar distribution in terms of the dynamic properties of the entire structure. They used an analytical model based on the analysis of classical equations of motion of thin plates elastically connected with connectors. Then, using genetic algorithms and applying dynamic optimization criteria, they obtained non-trivial topological solutions for the structure of VIG plates.
The motivation to undertake this work was the willingness to summarize and also continue the research contained in [18]. It is worth emphasizing that authors use a completely different theoretical model of VIG plates. Combining, comparing and analysing advanced FEM models with a simple analytical approach allows the selection of models of this type of structure that is easier to apply in the future.
Another motivation to start this work were the conclusions from the work [18]. The authors, after analysing the dynamic properties of many VIG-type structures, concluded that the optimization used criteria require further refinement and extension, as they slightly affect the distribution of pillars. Hence, the need to extend the optimization criterion with new aspects, which is undoubtedly the fractal analysis. By determining the correlation coefficient of a given structure with the mathematical structure of the fractal, it is possible to evaluate the aesthetic properties of the arrangement of pillars in the VIG panels. This would enable a more comprehensive approach to the optimization of these types of structures.

1.4. Purpose of the Work

The first aim of this work is to compare the accuracy of various numerical models of VIG plates in terms of dynamics. Demonstrated models are simple and easier to apply. They provide correct and accurate results in relation to advanced three-dimensional structure models. They allow for faster computations, which is important in the case of genetic calculations and optimization.
The second important goal of the work is the use of multi-criteria optimization to design the arrangement of connectors in VIG panels. By extending the objective function with a parameter that determines the degree of similarity of a given structure to a fractal, a simple quantitative parameter characterizing the aesthetic properties of panels was obtained. The relationship between the aesthetic properties of the given structure and its fractal properties has been investigated by other authors for various problems [10]. The distribution of pillars, which in the case of structural elements intended for use in exposed places such as window partitions, is of key importance for the comfort of users. To determine a parameter characterizing the degree of similarity in the arrangement of pillars on the VIG plate, and thus to investigate its aesthetic properties, the mass-radius method was used [19].

2. Numerical Modelling

2.1. Geometry and Support Conditions

The main issue in selecting the geometry and support conditions was the need to enable compatibility between all analysed modelling methods, which differ significantly in their approach. Only by providing the same support conditions obtained results were comparable to each other. Ultimately, simple support conditions were assumed on all edges of the plate (i.e., all displacements were blocked and all rotations were released at all nodes—Figure 2).
In the case of the three-dimensional FEM model, the displacements were blocked only in the internal panes of both plates (Figure 3). In the case of blocking the displacements along the entire thickness of the glass elements, restraint would have appeared on a given edge [17].
Support of the middle panes of both glass plates was not considered as this would significantly increase the stiffness of the entire model. The thickness of the glass plates is much greater (6–8 mm) than the vacuum thickness (0.3 mm). This makes the arm of horizontal support forces much (about 20 times) greater than the arm of corresponding forces in the simplified model.

2.2. Simple FEM Model

A simple FEM model was made using the commercial engineering software Autodesk Robot Structural Analysis Professional v2022. A shell-member model was built, in which the glass plates were represented as flat shell elements simply supported in their middle panes. Pillars, represented as bar elements with six degrees of freedom in each node, were placed between plates (Figure 4). The division of the regular FEM mesh of both plates was adjusted so the mesh stands out at the point of the pillar’s connection.

2.3. 3D FEM Model

The model implemented in ABAQUS based on three-dimensional C3D8R elements: first order, reduced integration element with hourglass control activated was selected as the most advanced 3D model using FEM [17] (Figure 5).
The interpolation function for C3D8R element shown in Figure 5 is given by the equation:
u = N I g , h , r u I   sum   on   I .
where I denote the node of the element. The last four vectors Γ α I , α = 1, 2, 3, 4 are the hourglass base vectors. The gradient matrix B I is defined by integrating over the element:
N I g , h , r = 1 8 Σ I + 1 4 g Λ 1 I + 1 4 h Λ 2 I + 1 4 r Λ 3 I + 1 2 h r Γ 1 I + 1 2 g r Γ 2 I + 1 2 g h Γ 3 I + 1 2 g h r Γ 4 I ,
The shape functions are the same as for the C3D8 element and can be found in [20] (2):
B i I = 1 V e l V e l   N i I g , h , r d V e l ,   N i I g , h , r = N I x i ,
where V e l is volume of the element and I = 1, 2, 3. In the centroidal strain formulation the gradient matrix is simply given:
B i I = N i I 0 , 0 , 0 .
Considering the above, it may be seen that centroidal strain formulation reduces the amount of effort required to compute the gradient matrix. In ABAQUS the artificial stiffness method and the artificial damping method showed in [21] are used to control the hourglass modes in these elements. In the work [22], the effectiveness of C3D8R elements with hourglass control was compared to other types of elements.

2.4. Analytical Model

The analytical model included in this study was based on the assumptions made in [18]. Let Ox1x2x3 be an orthogonal Cartesian coordinate system. Two plates of thicknesses hi, i = 1, 2 made of the same material (density ρ, Young’s modulus E and Poisson number ν) were considered. There are n springs, each of stiffness kj, located in points with coordinates x 1 , j , x 2 , j ,   j = 1 , , n are located between those plates (Figure 6).
Due to the small thickness of the glass plates and their slight deformations, the classical theory of thin plates, known as the Kirchhoff—Love theory, was applied. According to the studies conducted by the authors in [15], the results obtained using the linear and classical theory for frequency values are fully realistic. According to the theory of thin plates, the displacement takes the form of:
w α x = w α 0 x 1 , x 2 x 3 w 0 x α = w α 0 x 3 w , α 0 ;             α = 1 , 2 w 3 x = w 0 x 1 , x 2
strain-displacement relations are presented as:
ε α β = 1 2 w α , β 0 + w β , α 0 + w , α 0 w , β 0 x 3 w , α 0 ε α 3 = 0 ;   ε 33 = 0 .
The action functional of the problem is given by:
A = t 1 t 2 W K d t ,
where the strain energy is equal to:
W = 1 2 A i = 1 2 α β w i D i α β γ ω γ ω w i d A + 1 2 j = 1 n k j w 2 x 1 , j , x 2 , j , t w 1 x 1 , j , x 2 , j , t 2 ,
and kinetic energy is given by:
K = 1 2 A i = 1 2 μ i w ˙ i 2 d A .
Let w i = w i x 1 , x 2 , t be the transverse deflection; overdot stands for the derivative with respect to time t = t 0 , t 1 . The coefficients D α β γ δ i represent the plate stiffness tensor elements and are given by the following formula:
D α β γ ω i = D 0 i v δ α β δ γ ω + 1 v 2 δ α γ δ β ω + δ α ω δ β γ ,
where the plate flexural rigidity is equal to:
D 0 i = E h i 3 12 1 ± v 2 .
The derivatives are given as follows:
α w i w i x α ,     α w i 2 w i x α x β ,       w ˙ i = w i t .
Coefficient α β stands for Kronecker delta. The indices α , β , γ , ω are related to the Ox1x2 coordinate system and run through 1, 2.

2.5. Natural and Forced Vibrations

In this study, the dynamic properties of VIG plates were investigated in terms of natural vibrations in which various approaches to numerical modelling were compared and in terms of forced vibrations. The eigenproblem can be expressed in a standard form:
K ω 2 M u = 0 ,
where K, M and u stand for matrix of stiffness, matrix of mass and vector of displacement, respectively. The displacement vector is equal to:
u = a 1 1 a 1 2   a 1 m a 2 1 a 2 1   a 2 m .
Elements of matrices K and M are, respectively, given by following formulas:
K i , j = 2 W u i u j ,
M i , j = 2 K u i u j .
Considering forced damped vibrations the system of equations is assumed as:
Ku + C u ˙ + M u ¨ = Q   cos Ω t .
where K, C, M stand for stiffness, damping and mass matrices, respectively, Q is load, u is displacement vector Ω is forcing frequency and t is time. The damping matrix has been defined so that the damping could be considered as slight:
C = α MK 1 2 ,
where α is assumed equal to 10−2 [18]. The formula describing a uniformly distributed load is assumed as:
  Q x 1 , x 2 , t = j = 1 m q j φ j x 1 , x 2 c o s Ω t ,
where qj are known coefficients dependent on the loading case.
Application of the multimodal approach leads to the following solution:
u = a C c o s Ω t + a S s i n Ω t .
Coefficients aC and aS are defined as:
a C = K Ω 2 M + Ω 2 C K Ω 2 M 1 1 Q ,
a S = Ω K ω 2 M 1 C K Ω 2 M + Ω 2 C K Ω 2 M 1 1 Q .

2.6. Simple Genetic Algorithm (SGA)

SGA (simple genetic algorithm) was selected as the tool for topology optimization of the plate. A single individual (VIG plate) was encoded by genotype in the form of an array with values 1 or 0 representing the existence or absence of a pilaster at a given position in the plate [18]. An example of coded pillars arrangement is shown in Figure 7.
Other SGA parameters, such as the selection of an appropriate method of interbreeding of individuals taking into account the characteristics of the problem, are described in detail in [9].

3. Optimization Criteria and Fractal Analysis

In this paper, several optimization criteria of various types have been defined. The first criterion used was the dynamic one, defined analogously to the work [18], i.e., as the average of the ratio of the maximum deflection of the lower and upper plate in the frequency range from 220 Hz to 440 Hz:
f = a v e r a g e m a x w 2 m a x w 1 .
The function (23) for a solid plate (i.e., a plate with evenly spaced pillars) is not less than 0.985. However, as shown in [18], limiting ourselves to this criterion gives only a small possibility of topological optimization of the joint arrangement.
Fractal analysis was used as another criterion for advanced optimization of pillar distribution. It determines the degree of similarity of a given structure to the theoretical fractal structure. There are many methods of determining the so-called fractal dimension of the structure. Due to the point nature of the pillars distribution decision to use an algorithm for the objective function calculation based on the simplified mass-radius method described in [19] was made. This algorithm was implemented in two steps:
  • The so-called spinning radius R (Figure 8) varying from 0 to the maximum value depending on the size of the given plate, and for each R, the number of pilasters in total M(R) was counted in circles with centres in all pillars of a given spinning radius R;
  • Because for an ideal fractal structure there is a relationship [19]:
M R ~ R D ,
where D is the so-called fractal dimension.
In order to investigate the degree of similarity between VIG plate structure and the ideal fractal a logarithmic approximation of M(R) data series was used. Then, r coefficient was determined which defines the degree of adaptation the logarithmic curve to the data described in the previous step. Classical formulas for the least squares approximation were used:
r = n i = 1 n x i y i i = 1 n x i i = 1 n y i n i = 1 n x i 2 i = 1 n x i 2 n i = 1 n y i 2 i = 1 n y i 2 ,
where x i = log R and y i = log M R .
The correlation coefficient r obtained in (25) was used directly as an objective function. Its values are in the range from 0 to 1.
Equal weights and linear combination of the above objective functions were used for the multi-criteria analysis:
w = 0.5 f + 0.5 r .
The results obtained solely on the basis of the criterion defined by the Formula (21) were shown and analysed in [18]. This paper focuses on the results obtained for the criteria defined by the Formulas (23) and (24).
Figure 8. The idea of calculation of the fractal dimension of a structure based on the spinning radius [19].
Figure 8. The idea of calculation of the fractal dimension of a structure based on the spinning radius [19].
Materials 16 04722 g008

4. Calculation Results

4.1. Comparison of Various Numerical Models

Comparative numerical calculations were performed for two pillars systems. In both cases, a plate with a glass pane thickness of h = 7.5 mm, vacuum width hv = 0.3 mm and plate dimensions 2.0 m × 1.0 m was used. Glass with Young’s modulus E = 72 GPa, Poisson’s ratio ν = 0.33 and density ρ = 2500 kg/m3 was assumed. Support-pillars are 0.3 mm height and diameter. The material for pillars is steel, with Young modulus Es = 200 GPa, the mass density ρs = 7850 kg/m3, and the Poisson’s ratio νs = 0.3. The following pillars systems were used (Figure 9):
For the analysis using FEM, in addition to the given material data, the following approaches were adopted:
  • Simple FEM model—flat plate quadrilateral finite elements with a relatively sparse division of 50 mm.
  • The 3D FEM model—three-dimensional finite elements were used. The glass panels were divided with a mesh size of 2.5 mm (0.5 mm near the support pillars), while the support pillars themselves were divided more accurately.
In the case 0375, the following results of the first several natural frequencies were obtained (Table 1).
The 0375X was analysed as the second case, which differed in more pillars in the centre of the slab. In the case of the 0375X, the following results of the first sixteen natural frequencies were obtained (Table 2).
When comparing the values of individual frequencies of natural vibrations, the modes of natural vibrations obtained in all models were analysed. They can be divided into two groups: vibrations in phase and vibrations in counter-phase. In the former case, both plates swing in the same direction at the same time (Figure 10a–d), in the latter case, plates swing in the opposite direction (Figure 10e–h).
It is worth noting that in the case of vibrations in the phase, the pilasters have no significant effect on the shape and the result of the natural vibrations values because both panes vibrate identically—the pillars do not undergo major deformations. However, in the case of vibrations in the counter-phase, the pillars undergo significant deformations, and their parameters and modelling method are of key importance for the final result. In general, the obtained results are in line with expectations. In most cases, the results obtained from all three models are consistent with each other, with differences not exceeding 5–6%.
However, in both cases, both 0375 and 0375X, the results obtained for models I and II (analytical model and FEM model) in the case of the lowest vibration frequencies significantly differ from the values obtained for model III (FEM 3D)—even by 24% (marked green in Table 1 and Table 2). The results for the 3D model are significantly higher. This is due to a completely different method of supporting these models, which is of greatest importance precisely for the first forms and frequency of vibrations. As the thicknesses of the glass plates is much greater than the thickness of the vacuum, a realistic 3D model is much stiffer than the corresponding 2D model. It was observed after analysing the fragment of the first mode of natural vibrations located in the support zone (Figure 11).
The form of deformation, in the case of 3D FEM model (Figure 10a), resembles a fixed scheme more than a simply supported scheme (Figure 10b).
For some cases of counter-phase vibrations large discrepancies in the results between the models were obtained, i.e., 3 (Figure 10e) and 10 (Figure 10g). Method of modelling and meshing of pillars may have had an influence on them. This problem was analysed in [17].
The results from the 0375 and 0375X models were also compared. It can be clearly seen that, in the case of the 0375X, the use of pillars in the centre of the slab completely eliminated the 3rd form of vibrations in the 0375 model (Figure 10e, Table 1 and Table 2).

4.2. Multi-Criteria Optimization

The rest of the numerical calculations for a number of cases which varied the dimensions of the VIG plate were performed (length L1, width L2 and the thickness of the glass plate h). Other geometrical and material parameters were assumed as in the previous analysis. The list of case numbers along with the corresponding combinations of dimensions are summarized in Table 3.
The analytical model of VIG plate and optimization using previously described genetic algorithms were used for the calculations. The following pillars distributions inside the VIG plate were obtained (Figure 12).
For square plates, the optimal forms are point-symmetrical with rotational symmetry of 90° angle (cases 0125 and 0175). As the disproportion between plate length and width increases, the optimal forms develop toward axisymmetric systems. An interesting relationship for the cases with a thicker glass plate (0175 and 0275) may be noticed. Due to the greater stiffness of the individual glass panes and the consequent different degree of energy transfer between the plates correlated with the objective function (21), the optimal forms aim at the ones shown in Figure 13.
In the case of using the fractal objective function defined by the Formula (23), the optimal solutions become more complicated. It is clear that trivial solutions are rejected. In the case of a square plate (0125), it is characterized by quasi-rotational symmetry with 90° angle. However, in the case of rectangular plates (0225 and 0325) the evolution of forms in relation to those shown in Figure 13 towards more differentiated solutions, but with the preservation of some areas of a linear arrangement of pillars.

5. Summary and Conclusions

As a result of the numerical calculations and analyses the following conclusions were formed:
  • Natural frequencies analysis of VIG panels requires 3D model application, the results obtained from simplified models are significantly underestimated;
  • VIG panels are characterized by two types of vibration: in phase (both glass panes bend in the same direction) and in counter-phase (both glass panes bend in opposite directions), pillars, their geometry and modelling method are of key importance for vibrations in counter-phase;
  • It is possible to effectively extend the criteria for optimizing the arrangement of connectors with new design criteria;
  • Fractal analysis can be a tool for VIG panels design.
In further research, the authors intend to focus on refining the programming genetic algorithm and determining the fractal size of the structure, so that it is possible to significantly tighten the pilaster mesh, as well perform a detailed analysis of pillars parameters on the values and forms of natural vibrations.

Author Contributions

Conceptualization, D.K., Ł.D. and A.W.; methodology, Ł.D. and A.W.; software, I.K., D.K., S.K., M.R., Ł.D. and A.W.; validation, I.K., D.K., S.K., M.R., Ł.D. and A.W.; formal analysis, I.K., D.K., S.K., M.R., Ł.D. and A.W.; investigation I.K., D.K., S.K., M.R., Ł.D. and A.W.; resources, I.K., D.K., S.K., M.R., Ł.D. and A.W.; data curation, I.K., D.K., S.K., M.R., Ł.D. and A.W.; writing—original draft preparation, I.K. and A.W.; writing—review and editing, I.K. and A.W.; visualization, I.K., D.K., S.K., M.R., Ł.D. and A.W.; supervision, M.R., Ł.D. and A.W.; project administration, I.K., D.K., S.K., M.R., Ł.D. and A.W.; funding acquisition, I.K., D.K., S.K., M.R., Ł.D. and A.W. All authors have read and agreed to the published version of the manuscript.

Funding

The contribution is founded by Polish Minister of Education and Science under contract no. SKN/SP/496315/2021. This research was supported in part by PL-Grid Infrastructure.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the result.

References

  1. Collins, R.E.; Simko, T.M. Current Status of the Science and Technology of Vacuum Glazing. Sol. Energy 1998, 62, 189–213. [Google Scholar] [CrossRef]
  2. Cho, S.; Kim, S.-H.; Eom, J.-Y. The Evaluation of Thermal Performance of Vacuum Glazing by Composition and the Pillar Arrangement through Test Method of Thermal Resistance. J. Korean Sol. Energy Soc. 2015, 35, 61–68. [Google Scholar] [CrossRef] [Green Version]
  3. Lévy-Véhel, J.; Lutton, E. Fractals in Engineering: New Trends in Theory and Applications; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  4. Wu, Y.; Shao, Z.; Wang, F. Study on Wood Fracture Parallel to the Grains Based on Fractal Geometry. Int. J. Fract. 2012, 176, 163–169. [Google Scholar] [CrossRef]
  5. Hu, W.; Liu, Y.; Li, S. Characterizing Mode I Fracture Behaviors of Wood Using Compact Tension in Selected System Crack Propagation. Forests 2021, 12, 1369. [Google Scholar] [CrossRef]
  6. Mandelbrot, B.B. The Fractal Geometry of Nature/Revised and Enlarged Edition; WH Free. Co.: New York, NY, USA, 1983. [Google Scholar]
  7. Jiang, S.F.; Su, Y. Fractal Theory and Its Application in Civil Engineering. Eng. Mech. 2009, 26, 148–152. [Google Scholar]
  8. Vrdoljak, A.; Miletić, K. Principles of Fractal Geometry and Applications in Architecture and Civil Engineering. e-Zb. Electron. Collect. Pap. Fac. Civ. Eng. 2019, 9, 40–52. [Google Scholar]
  9. Wirowski, A. The Hybrid Method of Finding the Natural Frequency Spectrum for Quasi-Fractal Plate Bands. J. Appl. Math. Comput. Mech. 2015, 14, 153–165. [Google Scholar] [CrossRef] [Green Version]
  10. Coleman, R. Fractal Analysis of Stealthy Pathfinding Aesthetics. Int. J. Comput. Games Technol. 2009, 2009, 670459. [Google Scholar] [CrossRef]
  11. Rian, I.M.; Sassone, M. Fractal-Based Generative Design of Structural Trusses Using Iterated Function System. Int. J. Space Struct. 2014, 29, 181–203. [Google Scholar] [CrossRef]
  12. Park, J.; Oh, M.; Lee, C. Thermal Performance Optimization and Experimental Evaluation of Vacuum-Glazed Windows Manufactured via the in-Vacuum Method. Energies 2019, 12, 3634. [Google Scholar] [CrossRef] [Green Version]
  13. Hu, D.; Li, Y.; Liu, C.; Li, Y. Analysis for the Heat Transfer of Fully Tempered Vacuum Glazing Based on the Thermal Resistance Model and Finite Element Model. Adv. Mech. Eng. 2018, 10, 1687814018795985. [Google Scholar] [CrossRef] [Green Version]
  14. Simko, T.M. Heat Transfer Processes and Stresses in Vacuum Glazing. Ph.D. Thesis, The University of Sydney, Camperdown, NSW, Australia, 1996. [Google Scholar]
  15. Xi, X.; Shi, Y.; Shan, X.; Zhang, Y.; Shen, H.; Zhang, R. Mechanical Properties of Tempered Vacuum Glazing with Continuous Vacant Support Pillars. Vacuum 2021, 188, 110165. [Google Scholar] [CrossRef]
  16. Ashmore, N.; Cabrera, D.; Kocer, C. Acoustic Properties of Vacuum Insulating 2015. Available online: https://www.acoustics.asn.au/conference_proceedings/AAS2015/papers/p74.pdf (accessed on 8 June 2023).
  17. Kowalczyk, I.; Kozanecki, D.; Krasoń, S.; Rabenda, M. Computational Modelling of VIG Plates Using FEM: Static and Dynamic Analysis. Materials 2022, 15, 1467. [Google Scholar] [CrossRef] [PubMed]
  18. Kowalczyk, I.; Domagalski, Ł.; Wirowski, A. Placement Optimization of Micro-Support Pillars in Vacuum Glazing Using SGA in Terms of Dynamics Criteria. In Proceedings of the Proceedings M2D2021: 9th International Conference on Mechanics and Materials in Design, Funchal, Portugal, 26–30 June 2022; p. 19109. [Google Scholar]
  19. Omiotek, Z. The Use of the Fractal Dimension for Analysis of the Contour of Objects. Inform. Autom. Pomiary W Gospod. I Ochr. Środowiska 2012, 2, 8–11. [Google Scholar]
  20. Abaqus ABAQUS Version 6.6 Documentation. Available online: https://classes.engineering.wustl.edu/2009/spring/mase5513/abaqus/docs/v6.6/index.html (accessed on 10 February 2023).
  21. Xiao, B.; Wang, W.; Zhang, X.; Long, G.; Fan, J.; Chen, H.; Deng, L. A Novel Fractal Solution for Permeability and Kozeny-Carman Constant of Fibrous Porous Media Made up of Solid Particles and Porous Fibers. Powder Technol. 2019, 349, 92–98. [Google Scholar] [CrossRef]
  22. Cho, S. Analysis of the Performance of Vacuum Glazing in Office Buildings in Korea: Simulation and Experimental Studies. Sustainability 2017, 9, 936. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic diagram of a vacuum glazing.
Figure 1. Schematic diagram of a vacuum glazing.
Materials 16 04722 g001
Figure 2. VIG plate computational model—geometry.
Figure 2. VIG plate computational model—geometry.
Materials 16 04722 g002
Figure 3. The method of support used in the FEM 3D model.
Figure 3. The method of support used in the FEM 3D model.
Materials 16 04722 g003
Figure 4. A Shell-member model created in Robot Structural Analysis Professional.
Figure 4. A Shell-member model created in Robot Structural Analysis Professional.
Materials 16 04722 g004
Figure 5. C3D8R element designations. Numbers stand for nodes.
Figure 5. C3D8R element designations. Numbers stand for nodes.
Materials 16 04722 g005
Figure 6. Computational model of the plate.
Figure 6. Computational model of the plate.
Materials 16 04722 g006
Figure 7. Example of chromosome coding: (a) schematic plate with pillars (black point—pillar occurs, grey point—pillar lacking), (b) part of chromosome representing an individual plate.
Figure 7. Example of chromosome coding: (a) schematic plate with pillars (black point—pillar occurs, grey point—pillar lacking), (b) part of chromosome representing an individual plate.
Materials 16 04722 g007
Figure 9. Applied pillars arrangement in two calculation examples: (a) case number 0375, (b) case number 0375X. Grey dots stand for the lack of the pillars and black stand for pillars.
Figure 9. Applied pillars arrangement in two calculation examples: (a) case number 0375, (b) case number 0375X. Grey dots stand for the lack of the pillars and black stand for pillars.
Materials 16 04722 g009
Figure 10. Selected z-normalized wane vibration modes for plate 0375: (ad) vibration in phase, (eh) vibration in counter phase; (a) 1st mode shape, (b) 2nd mode shape, (c) 4th mode shape, (d) 11th mode shape, (e) 3rd mode shape, (f) 6th mode shape, (g) 10th mode shape, (h) 14th mode shape.
Figure 10. Selected z-normalized wane vibration modes for plate 0375: (ad) vibration in phase, (eh) vibration in counter phase; (a) 1st mode shape, (b) 2nd mode shape, (c) 4th mode shape, (d) 11th mode shape, (e) 3rd mode shape, (f) 6th mode shape, (g) 10th mode shape, (h) 14th mode shape.
Materials 16 04722 g010
Figure 11. Fragment of the first z-normalized wane mode of natural vibrations for case 0375 in two calculation models (a) 3D FEM model; (b) simple FEM model.
Figure 11. Fragment of the first z-normalized wane mode of natural vibrations for case 0375 in two calculation models (a) 3D FEM model; (b) simple FEM model.
Materials 16 04722 g011
Figure 12. Optimal pillars distributions for the mixed criterion for different cases of plate geometry. Grey dots stand for the lack of the pillars and black stand for pillars.
Figure 12. Optimal pillars distributions for the mixed criterion for different cases of plate geometry. Grey dots stand for the lack of the pillars and black stand for pillars.
Materials 16 04722 g012
Figure 13. Optimal pillars distributions for the fractal criterion for different cases of plate geometry. Grey dots stand for the lack of the pillars and black stand for pillars.
Figure 13. Optimal pillars distributions for the fractal criterion for different cases of plate geometry. Grey dots stand for the lack of the pillars and black stand for pillars.
Materials 16 04722 g013
Table 1. The first eighteen eigenfrequencies calculated in three different numerical models for the case number 0375. The frequencies that are most influenced by the problem of mapping with the method of supporting the 3D model are marked in green. Vibrations in counter-phase are marked in blue.
Table 1. The first eighteen eigenfrequencies calculated in three different numerical models for the case number 0375. The frequencies that are most influenced by the problem of mapping with the method of supporting the 3D model are marked in green. Vibrations in counter-phase are marked in blue.
Shape
Number
I. Analytical,
(Hz)
II. Simple FEM,
(Hz)
III. 3D FEM,
(Hz)
(I–II)/II,
(%)
(II–III)/III,
(%)
Figure Number,
[-]
124.2023.9331.56−12410a
238.7238.2043.53−11210b
358.1839.2656.30−483010e
462.9161.9765.38−2510c
582.2781.3977.94−1−4-
683.5385.4983.992−210f
7--92.66---
896.7995.2494.43−2−1-
996.7995.3495.92−21-
1094.1693.59107.46−11310g
11120.99118.59119.04−2010d
12139.23138.00130.76−1−6-
13--136.29---
14140.35139.41139.48−1010h
15145.59151.14149.694−1-
16154.86155.23149.840−4-
17169.41177.29171.854−3-
18--176.64---
Table 2. The first sixteen eigenfrequencies calculated in three different numerical models for the case number 0375X. The frequencies that are most influenced by the problem of mapping with the method of supporting the 3D model are marked in green. Vibrations in counter-phase are marked in blue.
Table 2. The first sixteen eigenfrequencies calculated in three different numerical models for the case number 0375X. The frequencies that are most influenced by the problem of mapping with the method of supporting the 3D model are marked in green. Vibrations in counter-phase are marked in blue.
Shape
Number
I. Analytical,
(Hz)
II. Simple FEM,
(Hz)
III. 3D FEM,
(Hz)
(I–II)/II,
(%)
(II–III)/III,
(%)
124.2023.9331.65−124
238.7238.2044.36−114
362.9161.9766.02−26
482.2781.3981.31−10
589.69-89.13--
690.32-89.71--
796.7995.2494.80--
896.7995.3496.50−21
9120.99118.60119.78−21
10140.35138.01131.57−2−5
11141.79138.10138.01−30
12--138.166−9
13141.85151,15150.27--
14154.86153.73---
15-158.87---
16179.06177.29172.81−1−3
Table 3. The list of case numbers with the corresponding combinations of dimensions.
Table 3. The list of case numbers with the corresponding combinations of dimensions.
Case NumberL1 (m)L2 (m)h (mm)
01251.001.002.50
01501.001.005.00
01751.001.007.50
02251.501.002.50
02501.501.005.00
02751.501.007.50
03252.001.002.50
03502.001.005.00
03752.001.007.50
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

Kowalczyk, I.; Kozanecki, D.; Krasoń, S.; Rabenda, M.; Domagalski, Ł.; Wirowski, A. Numerical Analysis, Optimization, and Multi-Criteria Design of Vacuum Insulated Glass Composite Panels. Materials 2023, 16, 4722. https://doi.org/10.3390/ma16134722

AMA Style

Kowalczyk I, Kozanecki D, Krasoń S, Rabenda M, Domagalski Ł, Wirowski A. Numerical Analysis, Optimization, and Multi-Criteria Design of Vacuum Insulated Glass Composite Panels. Materials. 2023; 16(13):4722. https://doi.org/10.3390/ma16134722

Chicago/Turabian Style

Kowalczyk, Izabela, Damian Kozanecki, Sylwia Krasoń, Martyna Rabenda, Łukasz Domagalski, and Artur Wirowski. 2023. "Numerical Analysis, Optimization, and Multi-Criteria Design of Vacuum Insulated Glass Composite Panels" Materials 16, no. 13: 4722. https://doi.org/10.3390/ma16134722

APA Style

Kowalczyk, I., Kozanecki, D., Krasoń, S., Rabenda, M., Domagalski, Ł., & Wirowski, A. (2023). Numerical Analysis, Optimization, and Multi-Criteria Design of Vacuum Insulated Glass Composite Panels. Materials, 16(13), 4722. https://doi.org/10.3390/ma16134722

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