Optimized Design of Floor Plan and Components of Prefabricated Building with Energy-Cost Effect

: Optimizing building performance and economic benefits through feedback in building design is a hot topic in current academic research. However, few studies on prefabricated buildings have been undertaken in this field. Meanwhile, the methodology used for achieving optimized solutions is still poor. In this paper, genetic algorithms and correlation analysis are employed and two parametric design methods—i.e., the floor plan generation method and the component selection method—are proposed for the modularity of the prefabricated buildings. Taking a typical high-rise building in Tianjin as an example, correlation analyses are performed on the basis of the two proposed methods to enhance the depth of the optimized finding approach. The outcome of this research demonstrates the feasibility of the proposed numerical approach, which can produce the optimized floor plan and construction set under the local conditions. This also reveals that the shape coefficient and window-to-wall ratio are strongly correlated with the energy performance of a building, which can help architects to pursue optimized design solutions in the schematic design process.


Introduction
The issue of global energy consumption is in the spotlight today. The U.S. Information Administration presents a future scenario in which global energy consumption will increase by nearly 50% over the next 30 years or so [1]. Building energy consumption accounts for 21.7% of the national energy consumption in China [2], and there is still work to be done to reduce overall energy consumption in the building industry.
Cost limits, of course, are crucial factors in restricting building energy efficiency [3]. The optimization of energy efficiency without the consideration of cost may result in high incremental costs and the inability to promote applications. As a result, energy usage and cost should be considered throughout the optimization process.
Many earlier researchers have investigated this topic using the genetic algorithm technique. I2n 2002, Caldas et al. used evolutionary algorithms on the DOE2.1E platform to improve the arrangement and size of windows in public buildings to lower the yearly energy consumption [4]. Ferrara et al. optimized near-zero-energy buildings to achieve a low economic cost [5]. Thalfeldt et al. identified the design priorities for cold-climate building facades [6]. However, related research has revealed that diverse optimization targets have conflicting characteristics [7], suggesting that focusing on a single goal is unsuitable in particular optimization schemes. As a result, an increasing number of researchers are attempting to adapt the multi-objective optimum design technique to the purpose of the architecture. With illumination and energy consumption as the optimization targets, Khoroshiltseva et al. employed modeFRONTIER and Daysim to optimize the • Few scholars have applied the synergy of energy consumption and cost to prefabricated buildings.

•
The models of the former studies can mainly be divided into two categories: one is a generic model with similar characteristics to that obtained from our research (it is usually a city building and is used to propose some common optimization conclusions [4,[8][9][10]); the other is generally a specific model, usually for a public building, and the findings primarily relate to the renovation and refurbishment of the building [5,6,11]. However, there are numerous phases in the architectural design process, including conceptual design, preliminary design, and detailed design [16].
In this article, we think that applying the two models to conceptual and preliminary design is most beneficial. • Many of the articles in this area end up focusing on the optimization results, while in practice designers tend to make changes based on these. These articles tend to lack any discussion of which parts need to be changed to have less impact on the optimization results.
Therefore, this paper proposes two design methods based on the genetic algorithm to take building energy consumption into consideration: a floor plan generation method for the conceptual design process and a complement selection method for the preliminary design. The designer can then use the results of the optimization and parameter correlation analysis as a theoretical basis to make further modifications to the computeraided optimization design. The total workflow of the two design methods is shown in Figure 1.

Optimization Method
The genetic algorithm (GA) was utilized as the optimization approach in this work; it has been frequently used in similar papers in recent years [17]. Genetic algorithms, which were researched and proposed by Professor Holland of Michigan University in 1975, are based on the theory of biological evolution and incorporate the evolutionary concepts of reproduction, hybridization, mutation, competition, and selection into the optimization process to achieve global optimization [18]. Figure 2 depicts the optimization concept [19]. This study uses the Galapagos and Octopus plug-ins integrated with the Grasshopper platform. Galapagos is a GA component that comes with the new version of Grasshopper and can perform optimization solutions for a single objective with a simple operation, fast computation, and easy convergence. However, it has the drawback of being able to optimize solutions for only one goal. Developed by the University of Applied Arts Vienna, Austria, and Bollinger + Grohmann Engineering, Germany, the Octopus plug-in is a Grasshopper component that combines Pareto frontier solution sets and GA for the optimization of multiple objectives. Pareto frontier solution sets can provide a basis for analyzing the trade-offs made between design objectives [20].

Optimization Method
The genetic algorithm (GA) was utilized as the optimization approach in this work; it has been frequently used in similar papers in recent years [17]. Genetic algorithms, which were researched and proposed by Professor Holland of Michigan University in 1975, are based on the theory of biological evolution and incorporate the evolutionary concepts of reproduction, hybridization, mutation, competition, and selection into the optimization process to achieve global optimization [18]. Figure 2 depicts the optimization concept [19]. This study uses the Galapagos and Octopus plug-ins integrated with the Grasshopper platform. Galapagos is a GA component that comes with the new version of Grasshopper and can perform optimization solutions for a single objective with a simple operation, fast computation, and easy convergence. However, it has the drawback of being able to optimize solutions for only one goal. Developed by the University of Applied Arts Vienna, Austria, and Bollinger + Grohmann Engineering, Germany, the Octopus plug-in is a Grasshopper component that combines Pareto frontier solution sets and GA for the optimization of multiple objectives. Pareto frontier solution sets can provide a basis for analyzing the trade-offs made between design objectives [20]. Appl The size of the absolute value of the standard regression coefficient, which is the regression coefficient obtained after eliminating the impact between the objective and the  The size of the absolute value of the standard regression coefficient, which is the regression coefficient obtained after eliminating the impact between the objective and the units of the parameters, directly represents the degree of effect of the parameters on the objective [21]. Its regression model can be expressed as Function (1): where β i is the ith parameter's standard regression coefficient, X is the ith parameter's mean, σ X is the ith parameter's standard deviation, ε is a constant, and N is the number of parameters. We need to test the problem of multicollinearity suggested by Frisch in 1934 over the course of the investigation [22]. To detect multicollinearity, a variety of approaches are used, including partial correlation coefficient, tolerance, variance inflation factor (VIF), and conditional index [23]. The variance inflation factor measured by the SPSS software is used to assess the aforementioned problem in this study. A result of greater than one and less than ten generally suggests that the problem does not exist [24].

Pearson Correlation Coefficient
The Pearson correlation coefficient, which is used to estimate the correlation between X and Y variables, can be calculated using Function (2). The coefficient takes values in the range of [−1, 1]: the closer it is to 1, the more likely the two variables are positively correlated; the closer it is to −1, the more likely it is that the two variables are negatively correlated; a value of 0 indicates that the two variables are uncorrelated. The article was followed up with calculations conducted using SPSS to obtain the Pearson correlation coefficient.
where cov(X, Y) is the covariance between X and Y and σ X and σ Y are their respective standard deviations.

Model
In this study, the thresholds of room bays, depths (Table 1), and the laws of arrangement of each functional space ( Figure 3) were summarized from several sets of house types, which were based on the prototype of prefabricated high-rise (one-staircase, two-family) commercial houses with a PC frame shear wall structure in Tianjin. After this, random values in the parameter range were used to form the base floor plan. -----Window-to-wall ratio of each room 0.1-0.9 Interleaved length between each functional space (m) −1.8-1.8 Following that, specific common attributes are assumed to finish the model's development, and their values are provided in Table 2. After the attributes are specified, a complete 3D building model ( Figure 4) can be created, which is the default model used in this study.  Following that, specific common attributes are assumed to finish the model's development, and their values are provided in Table 2. After the attributes are specified, a complete 3D building model ( Figure 4) can be created, which is the default model used in this study.

Parameters of the Floor Plan Generation Method
The control variables required for the floor plan generation method are the windowto-wall ratio of each orientation, the bays and depths of each functional space, and the interleaved length between them. The research object of this paper is the prefabricated building. One of the core tasks to promote the degree of industrialization of the building is to adopt a modular coordinated system to achieve the universality and interchangea-

Parameters of the Floor Plan Generation Method
The control variables required for the floor plan generation method are the windowto-wall ratio of each orientation, the bays and depths of each functional space, and the interleaved length between them. The research object of this paper is the prefabricated building. One of the core tasks to promote the degree of industrialization of the building is to adopt a modular coordinated system to achieve the universality and interchangeability of the components [27]. Therefore, a novel qualification method that considers the modulus of prefabricated buildings is proposed when defining the control variables of space types. This method can reduce the range of control variables based on the modulus to enhance the optimization efficiency and ensure the reasonable size of the prefabricated buildings. The objective function can be seen in Function (3).
where x is the vector of control variables; k i is the step vector of x i control variables; x m is the window-to-wall ratio for each orientation; a m and b m are the minimum and maximum values of the corresponding window-to-wall ratio; x i is the interleaved length of the bay, depth, and spaces, i ∈ {5, 6···n}, m; j is a constant in dimensional transformation which is used to expand the step length of the building modulus in a single transformation, j ∈ N + ; and M is the basic modulus of the building at the dimensional transformation of the prefabricated building (

Parameters of the Component Selection Method
The parameters used in this optimization are the material of the insulation layer, the thickness of the insulation layer, and the construction of the wall components, which can be divided into qualitative and quantitative indexes. The range of values shown Table 3 are determined by summarizing after researching manufacturers. Table 3. Ranges of the component selection method's parameters.

Parameters of Components Ranges
Material of insulation layer {Extruded polystyrene, expanded polystyrene, foamed polyurethane, rock wool} Construction of components {External insulation, internal insulation, sandwich insulation} Thickness of insulation layer in each direction (m) 0-0.5 In the optimization process of the component selection method, the physical properties and cost of each material are considered as attributes, assuming that they are constant throughout the construction phase. The specific attributes are shown in Table 4. In the optimization process, material and construction are qualitative indicators, and there are few desirable types. The optimal solution can be selected using the exhaustive method. On the basis of this solution, the thickness of the insulation material for each orientation is used as the optimization variable, while the total building energy consumption and total cost of insulation material are used as the optimization objectives for the next optimization step. The goal energy consumption in the floor plan generation design technique is confined to the use phase, since most buildings consume around 70% of their total energy over their whole life cycle during their use phase [28]. Lighting energy consumption and equipment energy consumption are not affected by control variables [29]. The annual energy consumption (y ec ) per unit area of the building optimized in this study can be calculated using Function (4): where E h is the annual heating energy consumption of the buildings, kWh; E c is the annual cooling energy consumption of the buildings, kWh; and A is the gross floor area of the buildings, m 2 .

Objective of Component Selection Method
g ec is the same as the expression of y ec in the floor plan generation method. Instead, only the cost of materials used during the construction phase is considered. The changes in the dimensions of the structural material will affect the energy use of buildings. We assume that the structural material cannot be changed during the selection process, which means that the cost of the structural material is constant for the same floor plan. In order to simplify the calculation, the material cost during construction is considered only as the insulation construction cost (g ic ), which can be expressed through Function (5).
where i from 1 to 4 are the four orientations of the buildings; C i is the price of the insulation board used for each facade orientation, CNY/m 3 ; d i is the thickness of the insulation board used in each direction of the external wall, m; and S i is the total area of each facade orientation, m 2 . The optimization objective function can be expressed as Function (6).
where x is the vector of control variables; g is the vector of objectives; and h s (x) is the sth constraint of the vector x, from which the feasible domain B is formed. After completing the above settings, we can obtain the details of the changes in parameters and properties for each phase of the complete workflow, as shown in Figure 5.
where is the vector of control variables; is the vector of objectives; and ℎ ( ) is the sth constraint of the vector , from which the feasible domain is formed. After completing the above settings, we can obtain the details of the changes in parameters and properties for each phase of the complete workflow, as shown in Figure 5:

Results and Discussions
Two design methods were applied on the above case and then the correlation between some parameters or attributes with the energy consumption was analyzed.

Optimization Results
In the floor plan generation method, the specific range of values for window-to-wall ratio, bays, depths, and interleaved length parameters could be determined using Function (3). In this case, = 3. The threshold was adjusted appropriately according to the limitation of ∈ ℤ. Finally, the bays, depths, and interleaved lengths were limited to a small range of values. The Galapagos parameters were set as stated in Table A1 during the optimization process, and convergence was mostly obtained at around 80 iterations, with the optimization ending after around 120 generations. Table 5 shows a comparison of the model before and after the plan optimization with improved energy consumption as the aim. Table A1 in the Appendix A shows the detail of optimization technique used.

Results and Discussion
Two design methods were applied on the above case and then the correlation between some parameters or attributes with the energy consumption was analyzed.

Optimization Results
In the floor plan generation method, the specific range of values for window-to-wall ratio, bays, depths, and interleaved length parameters could be determined using Function (3). In this case, j = 3. The threshold was adjusted appropriately according to the limitation of b i −a i jM ∈ Z. Finally, the bays, depths, and interleaved lengths were limited to a small range of values. The Galapagos parameters were set as stated in Table A1 during the optimization process, and convergence was mostly obtained at around 80 iterations, with the optimization ending after around 120 generations. Table 5 shows a comparison of the model before and after the plan optimization with improved energy consumption as the aim. Table A1 in the Appendix A shows the detail of optimization technique used. It can be seen that after the optimization, the building energy consumption and interleaved length between each functional space are reduced. It can be seen that after the optimization, the building energy consumption and inte leaved length between each functional space are reduced.

Number of Iterations 0 120
Plan shape Energy consumption (kWh/m ) 86. 18 18.11 In the optimization of the component selection method, the model was inherit from the previous optimization step without any modification and then the enumerati method was used to optimize the material and structure. The Octopus parameters we set as stated in Table A3 during the optimization process. The step size of the thickness the optimization process was 0.05 m and the distribution was between 0 and 0.15 m. T curves of the energy consumption and cost for different insulation materials with differe constructions is shown in Figure 6. It can be seen that after the optimization, the building energy consumption and interleaved length between each functional space are reduced.

Number of Iterations 0 120
Plan shape Energy consumption (kWh/m ) 86. 18 18.11 In the optimization of the component selection method, the model was inherited from the previous optimization step without any modification and then the enumeration method was used to optimize the material and structure. The Octopus parameters were set as stated in Table A3 during the optimization process. The step size of the thickness in the optimization process was 0.05 m and the distribution was between 0 and 0.15 m. The curves of the energy consumption and cost for different insulation materials with different constructions is shown in Figure 6.
Energy consumption (kWh/m 2 ) 86. 18 18.11 In the optimization of the component selection method, the model was inherited from the previous optimization step without any modification and then the enumeration method was used to optimize the material and structure. The Octopus parameters were set as stated in Table A3 during the optimization process. The step size of the thickness in the optimization process was 0.05 m and the distribution was between 0 and 0.15 m. The curves of the energy consumption and cost for different insulation materials with different constructions is shown in Figure 6.
(kWh/m ) 86. 18 18.11 In the optimization of the component selection method, the model was inhe from the previous optimization step without any modification and then the enumera method was used to optimize the material and structure. The Octopus parameters w set as stated in Table A3 during the optimization process. The step size of the thickne the optimization process was 0.05 m and the distribution was between 0 and 0.15 m. curves of the energy consumption and cost for different insulation materials with diffe constructions is shown in Figure 6. The energy-cost charts for the four orientations show the same pattern. Rock w board material and polystyrene board material were closer to the 0-coordinate poin The energy-cost charts for the four orientations show the same pattern. Rock wool board material and polystyrene board material were closer to the 0-coordinate point of the coordinate system-i.e., they were better than the other insulation materials in terms of their reducing energy consumption and cost. Considering the cost of fire protection, rock wool board was selected as the optimal material in the next dual-objective optimization. From Figure 6, it can be seen that the two innermost curves are for sandwich insulation, which indicates that the effect of sandwich insulation is better than that of external insulation and internal insulation. Thus, sandwich insulation was chosen as the construction method for the optimization determination.
The thickness was selected as the optimization object, and basic convergence was achieved after 10 iterations. The Pareto frontier solution set (Figure 7) was derived after reaching the maximum number of iterations, and correlation plots between insulation thickness, total energy consumption, and total cost for each orientation were derived (Figure 8).
From the Pareto frontier solution set, it can be seen that the cost kept increasing and the total energy consumption kept decreasing as the total thickness increased within a certain thickness range, in accordance with the objective law.
insulation and internal insulation. Thus, sandwich insulation was chosen as the construction method for the optimization determination.
The thickness was selected as the optimization object, and basic convergence was achieved after 10 iterations. The Pareto frontier solution set (Figure 7) was derived after reaching the maximum number of iterations, and correlation plots between insulation thickness, total energy consumption, and total cost for each orientation were derived (Figure 8).  From the Pareto frontier solution set, it can be seen that the cost kept increasing and the total energy consumption kept decreasing as the total thickness increased within a certain thickness range, in accordance with the objective law. tion method for the optimization determination.
The thickness was selected as the optimization object, and basic convergence was achieved after 10 iterations. The Pareto frontier solution set (Figure 7) was derived after reaching the maximum number of iterations, and correlation plots between insulation thickness, total energy consumption, and total cost for each orientation were derived (Figure 8).  From the Pareto frontier solution set, it can be seen that the cost kept increasing and the total energy consumption kept decreasing as the total thickness increased within a certain thickness range, in accordance with the objective law.

Correlation Analysis of Floor Plan Generation Method
Eleven influencing factors were selected as variables for regression analysis from the parameters that may affect energy consumption-namely, interleaved length, east window wall ratio, south window wall ratio, west window wall ratio, north window wall ratio, east exterior wall area, south exterior wall area, west exterior wall area, north exterior wall area, shape coefficient, and aspect ratio. After exporting 3680 sets of data from the optimization process and eliminating 3 sets of invalid data, we obtained 3677 sets of valid data. The validity of the selected variables was first determined by testing the multicollinearity; if the variance inflation factor (VIF) of each variable was tested to be less than 10 (Table 6), this meant that each variable had a certain degree of independence. A significance test was then performed, yielding an overall p-value of 0.0001 less than 0.05-i.e., the proposed model was valid at a 95% confidence interval. Additionally, the p-value for each variable (Table 7) was less than 0.05-i.e., each variable was significant at a 95% confidence interval. The larger the standard regression coefficient of a variable is, the more important the variable is under the same condition. The variables are ranked in Table 7, and it can be seen that parameters such as the shape coefficient and the window-to-wall ratio are more important than the area of the exterior walls of each orientation-i.e., when making adjustments, the radiant area of the walls and windows of each orientation can be appropriately altered while controlling the shape coefficient and the window-to-wall ratio.

Correlation Analysis of Component Selection Method
In the phase of correlation analysis, only the relationship between the percentage of the thickness of the insulation in each orientation to the total insulation thickness and the total energy consumption is required, given that the types of material and construction have already been determined (Table 8). Table 8. Comparison of the model before and after the generation of building shape.

Thickness of Insulation Layer in Each Direction/Total Thickness Pearson Correlation
East ratio −0.427 West ratio 0.361 South ratio −0.252 North ratio 0.318 A negative correlation could be observed between the insulation thickness in the east/south directions and total energy consumption-i.e., increasing the proportion of the insulation in the east and south directions will decrease the total energy consumption; conversely, increasing the proportion of the insulation in the west and north directions will increase the total energy consumption. From a correlation point of view, in order to reduce the total building energy consumption, the insulation of east-and south-oriented buildings should be appropriately increased and the insulation of west-and north-oriented buildings should be reduced under a certain cost limit.

1.
Simulation-based single-objective or multi-objective optimization can be performed for prefabricated buildings. Unlike traditional buildings, the building modulus and component selection need to be considered in the optimization process. This not only meets the demand for the standardization of prefabricated buildings, but also increases the speed of optimization computation through reducing the number of values taken from parameters.

2.
A novel, modular parametric modeling approach was proposed and applied in the floor plan generation method. After this, the optimal generation of prefabricated highrise buildings in Tianjin was completed based on this method. The correlation between each parameter and energy consumption was also studied, and it was concluded that the shape coefficient and window-to-wall ratio are the main factors affecting the energy consumption of the buildings in Tianjin.

3.
A preliminary component selection method based on computer simulation was proposed-i.e., the component selection for the prefabricated building was mainly carried out to determine the construction of exterior walls, the selection of insulation materials, and the thickness of the insulation layer. By optimizing the generated models, it was finally concluded that sandwich insulation constructions and rock wool board insulation materials should be selected for buildings in Tianjin. According to the correlation analysis, the thickness of the insulation material in the east and south directions should be increased under a certain cost limit in order to reduce the total energy consumption of buildings.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Table A1. Parameter setting of Galapagos.

Project The Maximum Number of Iterations Population Size Multiplier of Initial Boost
Value 120 30 2 Project Proportion of retained elites Crossover ratio value 5% 75%