Structural Optimization and Application Research of Alkali-Activated Slag Ceramsite Compound Insulation Block Based on Finite Element Method

: The research and application of new wall materials have been attracting increasing attention owing to the continuous promotion of sustainable development in the building industry. An alkali-activated slag ceramsite compound insulation block (AASCCIB) is used as the research object. Based on the ﬁnite element method, the effects of different numbers of hole rows and hole ratios on the thermal and mechanical performances of AASCCIBs are analyzed using ANSYS CFX. On this basis, the AASCCIB with the optimal comprehensive performance is determined by a multi-objective optimization analysis. Finally, the improvement effect of the AASCCIB wall on the indoor thermal environment relative to an ordinary block (OB) wall is quantitatively analyzed using ANSYS CFX. The results show that the von Mises equivalent stress and heat transfer coefﬁcient of the AASCCIB decrease with the increase in the hole ratio when the hole shape and number of hole rows are constant. AASCCIB B 1 has the optimal comprehensive performance among six AASCCIBs, with the heat transfer coefﬁcient and average von Mises equivalent stress of 0.446 W/(m 2 · K) and 9.52 MPa, respectively. Compared with the indoor lowest and average temperatures of the building with the OB wall, those of the building with the AASCCIB wall increased by at least 1.39 and 0.82 ◦ C on the winter solstice, respectively. The indoor temperature difference decreased by at least 0.83 ◦ C. In addition, the indoor highest temperature, average temperature, and temperature difference decreased by at least 1.75, 0.79, and 1.89 ◦ C on the summer solstice, respectively.


Introduction
With the progress of science and technology and improvement in living standards, higher requirements for living conditions have been imposed. Instead of limiting the demand for buildings to shelter from wind, rain, and warmth, people emphasize the need for a comfortable, energy-saving, and environmentally friendly living environment [1,2]. However, improvements in indoor thermal comfort are often accompanied by increases in building energy consumption and environmental pollution [3,4]. It is reported that buildings account for 40% of the total global energy consumption, and this proportion is still increasing [5,6]. In addition, approximately 40% of the total human greenhouse gas emission is attributed to the building industry [7,8]. Therefore, how to achieve improvements in indoor comfort without significantly increasing energy consumption and greenhouse gas emissions is a thorny issue that countries must face and solve.
In China, wall materials account for approximately 70% of housing construction materials. Clay bricks are the predominant wall material. However, it is statistically demonstrated that more than 1 billion cubic meters of arable land and 70 million tons of coal are consumed each year by the production of clay bricks [9]. If the growth of clay brick production is uncontrolled, the contradiction between supply and demand of resources, as well as energy, will become more acute. Obviously, it is required to promote the innovation of wall materials to realize the sustainable development of buildings. The self-insulation block is a new type of wall material [10]. It occupies a large share of the construction market because of its good heat insulation and thermal insulation, lightweight, high strength, long service life, simple construction, etc. [11]. Commonly, cement is chosen as the cementitious material for self-insulation blocks [12]. The production and processing of cement, however, consumes large amounts of energy and emits a lot of greenhouse gases. The production of 1 t of cement requires 5000 MJ of energy and emits 1 t of CO 2 [13]. Thus, the cementitious material alternative to cement is sought.
In recent years, alkali-activated slag cementitious material (AASCM) obtained using slag (industrial waste) attracted increasing attention [11]. Compared to the traditional cement production process, the AASCM production process transfers from "two grindings and one burning" to "one grinding". This reduces the energy consumption and greenhouse gas emissions of cementitious materials in the production process [14,15]. Additionally, AASCM has the advantages of high strength, high-temperature resistance, frost resistance, corrosion resistance, etc. [15]. Evidently, it is important to achieve sustainable development of building materials by choosing AASCM to replace traditional cement.
Ceramsite is a typical representative of artificial light aggregates and new energysaving building materials [16]. It is obtained by an industrial solid waste via hightemperature calcination. It has the advantages of recycling waste, environmental protection, convenient production, and low price [17]. It also has a low thermal conductivity, high strength, lightweight, high-temperature resistance, acid and alkali resistance, etc. [18]. Considering these excellent characteristics of ceramsite, it can be used as an aggregate in studies on new wall materials.
To sum up, it is an inevitable choice for wall material development in studies on alkali-activated slag ceramsite self-insulation blocks (AASCSIBs). To improve the thermal performance of the alkali-activated slag ceramsite single self-insulation block, the alkaliactivated slag ceramsite self-insulation hollow block (AASCSIHB) is generally filled with an insulation material. That is, an alkali-activated slag ceramsite compound insulation block (AASCCIB) is prepared. Therefore, conducting research on AASCCIBs will promote the integrated development of energy efficiency and functionality of buildings.
Nevertheless, there are a few reports on AASCCIBs [11,19] which focus on the preparation and performance of the block. In other words, there are very few studies devoted to the structural optimization and application feasibility of AASCCIBs [19], especially those dedicated to the comprehensive optimization of mechanical and thermal performances and the exploration of the thermal insulation effect. Undoubtedly, the above limitations restrict the application and promotion of AASCCIBs.
The most commonly used research approach in the structural optimization and application effect analysis of blocks is numerical simulation [20][21][22]. Various numerical analysis methods are used, including the finite difference method (FDM) [23,24], boundary element method (BEM) [25,26], finite volume method (FVM) [27,28], and finite element method (FEM) [29,30]. The FDM is a dominant numerical method for solving the motion of objects in computational fluids. However, it is challenging to solve the boundary conditions using this method. The BEM transforms the solution of the entire domain into a solution on a regional boundary. This numerical method reduces the computational effort, but is inefficient for the computation of the complex-shape flow field. In addition, the method is generally suitable for only solving homogeneous linear problems. The largest advantage of the FVM is that it provides an accurate integral conservation even with coarse meshes. However, the accuracy of the FVM is only second-order. The FEM is a numerical technique used to obtain approximate solutions to boundary value problems of partial differential equations. The basic idea of the method is to discretize the continuous solution domain, i.e., to divide the continuum into a finite number of tiny blocks with regular shapes. The method is not only applicable to complex geometries and boundary conditions, but also has a high computational accuracy and wide applicability. Furthermore, it has a standard-ized calculation format and can be easily applied programmatically. Thus, the FEM is a powerful, effective, and accurate numerical method.
ANSYS integrates the analyses of structure, fluid, electric, magnetic, and acoustic fields [31]. As far as professional computer-aided engineering software is concerned, ANSYS is the only analysis and design software worldwide that has passed the IS09001 quality certification. Additionally, the advantages of ANSYS are reflected in the broad scope of analysis, powerful coupling analysis functions, and convenient co-simulation platform. ANSYS CFX, a branch of ANSYS, is dedicated to fluid dynamics simulation [31]. Notably, it is the first commercial software in the world to develop and use a fully implicit multi-grid coupled solver technique. Moreover, it has advanced algorithms, rich physical models, and accurate calculation results. Considering the above analysis, ANSYS CFX is recommended for the structural optimization and application effect analysis of AASCCIBs.
In summary, this study aims to optimize the structure and analyze the thermal insulation effect of the AASCCIB using the FEM. The specific work is as follows: • Six types of AASCCIBs with different internal structures are designed based on different numbers of hole rows and hole ratios.

•
Based on the FEM, the thermal and mechanical performances of six AASCCIBs are simulated using ANSYS CFX. Moreover, the AASCCIB with the optimal comprehensive performance is determined though a multi-objective optimization analysis.

•
The improvement effect of the AASCCIB wall on the indoor thermal environment relative to an ordinary block (OB) wall is quantitatively analyzed using ANSYS CFX.

Structural Design of AASCCIBs
Based on the literature [32], the dimensions of the selected blocks were length × width × height = 390 mm × 240 mm × 190 mm. Furthermore, because the thermal resistance of a hollow block with a rectangular hole was larger than those of diamond, square, and circle at the same hole ratio [33], the hole shape was rectangular in this study. On this basis, six types of hollow blocks with different internal structures were designed based on different numbers of hole rows and hole ratios. The diagram of the AASCSIHBs is shown in Figure 1. The distance between rows of the rectangular holes was the same. The distance between the rectangular holes parallel to the length direction of the block was referred to as horizontal rib, and the value was c. The distance between the rectangular holes perpendicular to the width direction of the block was referred to as vertical rib, with a value of f. The length and width of the larger rectangle in each row of the hole were b and e, respectively. The length and width of the smaller rectangle were denoted by d and e, respectively. In addition, a was taken as 25 mm according to the reference [32]. The detailed dimensions of the AASCSIHBs are shown in Table 1. Note: The hole ratio was the ratio of the hole area to the cross-sectional area of the block.

Numerical Simulation of Thermal and Mechanical Performances of AASCCIBs
To improve the thermal performance of the block, the holes of the AASCSIHB were filled with an extruded polystyrene (XPS) foam board, which yielded the AASCCIB. Based on the FEM, the thermal and mechanical performances of the six AASCCIBs were numerically analyzed.

Mathematical Model
The heat transfer calculations for this study were performed on a single block wall. In addition, a single block in the wall was selected for the heat transfer analysis of the wall.
According to Fourier's law and energy conservation equation, the differential equation of heat conduction in the three-dimensional (3D) Cartesian coordinates system was established. In the heat transfer analysis, the heat input and output to the differential cube over time dτ were as follows [34]: The heat input was: The heat output was: The net heat between input and output was as follows: According to Fourier's law [35], By substituting Equation (6) into Equation (5), we obtained Equation (7).
The heat generation of the infinitesimal cube over time dτ was q v d x d y d z · dτ [34]. The increment in the thermodynamic internal energy of the microelement over time dτ was ρc ∂T ∂τ · d x d y · dτ [34].
According to the law of conservation of energy, Considering that, p, c, and λ are constants, as well as there being no internal heat source inside the wall, Equation (9) was simplified to Equation (10). That is, the mathematical model was established.
where T is the temperature, • C. x, y, and z are the 3D Cartesian system coordinates.
The mechanical performance in this study was analyzed using the constitutive model with the von Mises yield criterion. The core concept of this criterion was that when the second invariant (J 2 ) of the stress deviation tensor of a point in a stressed object reached a constant value, the point entered the plastic state.
This criterion can be expressed by the principal stresses [36]: A unidirectional tensile experiment showed that c = σ s 2 3 [36]. According to a pure shear experiment, it was found that c = K 2 [36]. Thus, the von Mises yielding criterion in the principal coordinate system was: In addition, according to the von Mises yield criterion, the material began to yield when the equivalent force reached a constant value, which can be expressed as [37]: where σ is the equivalent stress (Pa), σ i (σ ii ) is the normal stress (Pa), σ s is the yield point, τ ij is the shear stress (Pa), K is the shear yield strength (Pa), The first and second subscripts of the stress components indicated the normal direction of the action plane and direction of the stress action, respectively.

1.
Basic assumptions and geometric model The following assumptions were used before the thermal model was developed: (1) The block performs a one-dimensional heat transfer.
The temperatures on both sides of the block are constant.
The main material of block and the filling material are closely connected. The material performances do not vary with the thermal environment.
Using AASCCIB B 1 as an example, a geometric model for the thermal analysis was developed ( Figure 2). The thermal insulation material (XPS) inside the holes of self-insulation hollow blocks had a small influence on the overall mechanical performances of the compound insulated block. Therefore, the self-insulation hollow block was selected instead of the compound insulation block for the mechanical performance analysis. The geometric model for the mechanical analysis was developed using AASCSIHB A 1 as an example ( Figure 2).

Mesh division
The mesh division of the developed geometric models was conducted. The results are shown in Figure 3.

3.
Material parameter and boundary condition setting The physical parameters of the AASCCIB are listed in Table 2. The block surfaces in contact with air were subjected to convective heat transfer boundary conditions (the third boundary condition). The other surfaces of the block were set with adiabatic boundary conditions. The expressions for the third boundary condition were as follows [38]: − λ 1 ∂t where λ 1 and λ 2 are the thermal conductivities of the hollow block and air layer, respectively (W/(m·K)). h i is the internal surface coefficient of heat transfer with a value of 8.7 W/(m 2 ·K) [39], h e is the external surface coefficient of heat transfer with a value of 23.0 W/(m 2 ·K) [39], t i1 is the indoor air temperature of 281. 15 K [40]. t i2 is the outdoor air temperature of 276.25 K [41], and t e1 and t e2 are the temperatures of the internal and external surfaces of the wall, respectively(K). In the mechanical analysis, a static analysis was used to inform the effect of geometric nonlinearity. All degrees of freedom at the bottom of the finite element model were constrained, and a uniform load (P) of 10 MPa was applied to the upper surface of the model (Figure 4).  In this study, the heat transfer coefficient, equivalent stress, and hole ratio represented the thermal, mechanical, and economic performances, respectively. Thermal, mechanical, and economic performances were included in the comprehensive performance. To obtain the AASCCIB with the optimal comprehensive performance, multi-objective optimization was conducted on the compound insulation blocks with different internal structures. The function expression of the weighted summation method was as follows [42]: where f(x) is the objective function, ω i is the weight coefficient, and f i (x) is the subobjective function.

Calculation of the Weight Coefficient
The weight coefficients of the block were calculated using a hierarchical analysis. According to the relative importance, weights were assigned to each influencing factor (heat transfer coefficient, equivalent stress, and porosity). In addition, a scale of 1 to 9 was used for importance comparisons of pairwise factors ( Table 3). The elements were compared to each other to obtain the judgment matrix A, as shown in Formula (19) [43]: where A is the judgment matrix and a ij is the scale. The maximum eigenvalue of the judgment matrix was calculated, and then a consistency test was performed. The Formulas were as follows [43]: where CI is the consistency index, RI is the average consistency index, λ max is the maximum eigenvalue, n is the order of the judgment matrix, and CR is the consistency ratio. The elements of the matrix A were normalized to obtain the matrix γ. The rows of the matrix γ were summed to obtain the column matrix W. Finally, the matrix W was normalized to obtain the ranking weight vector ω. That is, the percentages of thermal, mechanical, and economic performances that affected the comprehensive performance of the blocks were obtained.

Factor Normalization
The units, sizes, and orders of magnitude of the influencing factors were not consistent, which did not help evaluate the comprehensive performance of the blocks. Therefore, it was necessary to standardize the influencing factors using a normalization method. The calculation Formula for the normalization method was as follows [44]: where y i is the dimensionless constant, and n ∑ i=1 y i is equal to 1 [44].

Influence of AASCCIB Wall on the Indoor Thermal Environment
Through the ANSYS CFX simulation, the improvement effect of the AASCCIB wall on the indoor thermal environment relative to the OB wall was informed.
This study considered rural buildings in Southern Shanxi, China, as the simulation object. The main reason was that Southern Shanxi, China, is located in the zone of hotsummer and cold-winter, where it is hot in summer and wet with cold in winter. Moreover, the region is close to severe-cold and cold zones, where the average daily temperature reaches about −10 • C in winter [45]. Obviously, the indoor thermal environment of rural buildings in this region needs to be improved urgently. Based on the above analysis, this study took rural buildings in Southern Shanxi, China, as the simulation object.

Mathematical Model
According to Section 2.2.1, the 3D unsteady differential equation for heat conduction in the Cartesian coordinate system was as follows: This equation had to be solved in domain Ω. The boundary conditions for this domain were as follows [46]: where ρ is the density (kg/m 3 ), c is the specific heat capacity (W/(kg·K), τ is the time (s), q v is the heat flux (W/m 2 ), T is the temperature ( • C), λ i is the thermal conductivity along the i-axis direction (W/(m·K), n i is the cosine of the outer normal of the boundary along the i-axis direction. q is the heat flux density on the Γ 2 boundary, (W/m 2 ), T is the temperature on the Γ 1 boundary ( • C), T a is the external ambient temperature under natural convection conditions ( • C), T a is the absolute temperature of the boundary layer under forced convection conditions ( • C), and Γ is the whole boundary of domain Ω.

1.
Basic assumptions and geometric models To simplify the model calculation, the following assumptions were used before modeling: (1) The air in each indoor room is considered as a whole. Its heat transfer mode is natural convection heat transfer.
The indoor door is open. The door is half-open.
The influence of indoor human activities and electrical appliances on the indoor temperature is ignored.
We conducted a field study on rural buildings in Ankang (climatic conditions and building characteristics were representative of rural regions in Southern Shanxi [45]) during the period from 2017 to 2019. Based on the field research and reference [40], a representative floor plan for rural buildings in Southern Shanxi, China, was constructed ( Figure 5). Further, the geometric model of the representative building was developed using ANSYS CFX, as shown in Figure 6.

2.
Mesh division In the mesh division, the external surface of the air and internal surface of the external wall, window, door, and roof were mesh-refined. The minimum element size was 0.009 m, and the maximum was 0.967 m. The refined mesh model is illustrated in Figure 7.

Material parameter setting
Based on the field investigation, experimental research, and reference [39], the thermal parameters of the envelope of the building with OB walls and building with AASCCIB walls were obtained (Table 4).

Boundary condition setting
The combined outdoor temperature considered the combined effect of solar radiation and outdoor air temperature on the external envelope of the building. Thus, the comprehensive outdoor temperature was chosen as the boundary condition for the external surface of the wall. Considering that outdoor meteorological parameters generally vary periodically, the comprehensive outdoor temperature was expressed as a sine or cosine function. Further, the comprehensive outdoor temperature was fitted as a periodic sine or cosine function by Fourier series expansion [45]. In addition, the initial temperature of the indoor air body was set to 278.6 K, and the initial temperature of the external surface of the wall was set to 277.1 K in the simulation.
where A o is the zero-order outdoor disturbance ( • C), A n is the amplitude of the external disturbance of the n -th sine wave ( • C), nw is the frequency of the external disturbance of the n -th sine wave, nw = 2π n/T (rad), ϕ n is the initial phase of the external disturbance of the n -th sine wave (rad). T is the period of the function (h), and n is the order of the harmonic. Formula (29) was expanded by the sine, which yielded Formula (30).
Assuming that a 0 2 = A 0 , a n = A n sin ϕ n , and b n = A n cos ϕ n , [a n cos(nωτ) + b n sin(nωτ)] In addition, ϕ n was determined by a n and b n , b n > 0 and a n > 0, ϕ n = arctan a n b n (32) b n > 0 and a n < 0, ϕ n = arctan a n b n + π 2 (33) b n < 0, ϕ n = arctan a n b n + π (34) b n < 0 and a n = a n , ϕ n = π 2 (35) b n < 0 and a n = −a n , ϕ n = 3π 2 Because the external disturbance function of the building envelope was extremely complex, the variation in the external disturbance function was expressed by a series of discrete data with equal intervals.
If the basic period interval (0, T) was divided into N equal parts, the interval only had N discrete point values.
The rectangular superposition summation method was used to obtain the function expressions of a 0 2 , a n , and b n .
where N is the number of measurement points in a period with a value of 24, T is the period, with a value of 24 (h), n is the order of the harmonics with a value of 4. j is the serial number of the sampled values, and ∆τ is the sampling time interval with a value of 1 (h). Therefore, Formula (31) was transformed into Formula (41).
Based on the above analysis, according to the reference [47] and Formula (41), the fourth-order Fourier series expressions for the outdoor comprehensive temperature of different orientations in Ankang region were obtained.

Validation of the model
The correctness of the model was verified to ensure that the simulation results were accurate. The measured data were compared to simulated data. The measured data were the indoor temperatures of a representative rural building in Ankang, Southern Shanxi, China. The results are presented in Figure 8. Each temperature value in the Figure is the average temperature of the room. Figure 8 shows that the measured temperature was highly consistent with the simulation temperature. The R-squared between the measured and simulated temperatures was 0.91. At the same time, it was found that there was a small deviation between the simulated temperature and measured temperature. This deviation mainly originated from the simplifications and assumptions in modeling and accuracy error of the test instrument. Therefore, it was effective to use ANSYS CFX software to develop a finite element model of the building.

Analysis of Thermal and Mechanical Performances of AASCCIBs
Through the simulation by ANSYS CFX, the effects of the different numbers of holes and hole ratios on the thermal and mechanical performances of AASCCIBs were evaluated.

Thermal Performances
Contour plots of temperature and heat flux for the AASCCIB are presented in Figure 9 (for AASCCIB B 1 as an example).  Figure 9a shows that the temperature contour plot of the AASCCIB obeyed the fundamental law of heat transfer. The temperature decreased step-by-step from indoor to outdoor along the Y direction. As shown in Figure 9b, there was a concentration of heat flux at the vertical ribs, which could easily produce the thermal bridge effect. The main reason was that the heat transfer coefficient of the vertical ribs was larger than that of the other parts, which led to a quick heat transfer from the vertical rib of the block.
In addition, the temperature difference and heat flux intensity on both sides of the flat wall of the AASCCIB were obtained by numerical simulations. Furthermore, the heat transfer coefficient of the AASCCIBs was obtained by Formula (50) [48]. The results are listed in Table 5.
where K is the heat transfer coefficient (W/(m 2 ·K)), q is the heat flux intensity, (W/m 2 ), and ∆t is the temperature difference between the internal and external surfaces of the wall (K).  Table 5 shows that the heat transfer coefficient of the AASCCIB decreased with the increase in the hole ratio when the hole type and number of hole rows were constant. This was because the increased hole ratio led to the increase in filler thickness when the hole shape and number of hole rows were constant. This directly increased the block heat transfer hindrance, thereby decreasing the heat transfer coefficient.
In addition, the simulated heat transfer coefficient was compared to the theoretical values to validate the finite element model.
The theoretical value of the average heat transfer coefficient of the AASCCIB was calculated by Formulas (51) to (55) [39]: where R is the thermal resistance of materials ((m 2 ·K)/W), d is the thickness of each layer material (m), λ is the thermal conductivity of each layer material (W/(m·K)), R 0 is the total thermal resistance of the flat wall ((m 2 ·K)/W), R i is the internal surface resistance of heat transfer with the value of 0.11 (m 2 ·K)/W [39], R e is the external surface resistance of heat transfer with the value of 0.04 (m 2 ·K)/W [39], R is the average heat transfer resistance of the combined flat wall ((m 2 ·K)/W), F 0 is the total heat transfer area perpendicular to the direction of heat flux (m 2 ), F 1, F 2, . . . F n is the area of each heat transfer region parallel to the direction of heat flux (m 2 ), R 0,1 , R 0,2, . . . R 0, n is each heat transfer area parallel to the direction of the heat flux ((m 2 ·K)/W), φ is the correction factor with the value of 0.86 [39], and K is the average heat transfer coefficient of the combined wall (W/(m 2 ·K)). A calculation diagram of the average thermal resistance is shown in Figure 10. A division diagram of the heat transfer channels is shown in Figure 11. The average heat transfer coefficient of the AASCCIB was obtained using Formulas (51) to (55).  Based on the above calculation, the average heat transfer coefficient of AASCCIB A 1 was 0.435 W/(m 2 ·K). In addition, using the above calculation method, the average heat transfer coefficients of AASCCIBs A 2 , B 1 , B 2 , C 1 , and C 2 were 0.444, 0.424, 0.448, 0.434, and 0.442 W/(m 2 ·K).
The simulated and calculated values of the heat transfer coefficient were compared, as shown in Figure 12.  Figure 12 shows that the calculated and simulated values of the average heat transfer coefficient were almost equal. The relative error of AASCCIB A 1 was the largest (5.5%). The main reason for this error was the series of assumptions in the modeling. Therefore, the finite element building model obtained using ANSYS CFX was considered effective.

Mechanical Performances
Contour plots of deformation and von Mises equivalent stress for AASCSIHB are presented in Figure 13 (for AASCSIHB A 1 as an example). The simulation results of the von Mises equivalent stress of the AASCSIHBs are listed in Table 6.   Figure 13a presents that a large deformation at the edge of the block due to the principal stresses acted on it. Figure 13b shows that the minimum von Mises equivalent stress appeared in the cross ribs and inner wall of the block. The maximum von Mises equivalent stress appeared at the corners of the block (stress concentration area), and the area where stress concentration occurred was small. In summary, in this study, the maximum von Mises equivalent stress had no decisive influence on the overall bearing capacity and stability of the member. Table 6 indicates that the average von Mises equivalent stress tended to decrease with the increase in the hole ratio when the hole shape and number of hole rows were fixed. The main reason was that the thickness of the middle rib of the block decreased with the increase in the hole ratio when the hole shape and the number of hole rows were constant. This reduced the compressive performance and connectivity between the concrete rib and wall of the block, thereby leading to a decrease in the average von Mises equivalent stress of the block.

Determination of Optimal AASCCIBs
According to Formula (19) and Table 3, judgment matrix A was obtained: Based on Formulas (20) and (21), λ max and CI were 3.092 and 0.0046, respectively. According to reference [49], RI was 0.58. Based on Formula (22), CR was 0.0089. As CR was below 0.1, judgment matrix A satisfied the consistency test.
Matrix A was transformed to obtain the ranking weight vector ω (ω= (0.539 0.297 0.164) T ). Thus, the percentages of thermal, mechanical, and economic performances that affected the comprehensive performance of the AASCCIB were 53.9%, 29.7%, and 16.4%, respectively.
The heat transfer coefficient, Mises equivalent stress, and porosity were transformed to dimensionless using the normalization method. The results are listed in Table 7. The weighted comprehensive value of AASCCIBs was obtained by Formula (18). The results are presented in Figure 14. Note: A smaller heat transfer coefficient indicates better thermal performance. Therefore, the dimensionless value of the heat transfer coefficient was multiplied by −1 [44], to make the monotonicity of the three sub-factor functions of heat transfer coefficient, equivalent stress, and hole ratio consistent.  Figure 14 shows that the weighted comprehensive value of AASCCIB B 1 was largest. Thus, compared with several other AASCCIBs, AASCCIB B 1 had the optimal comprehen-sive performance. Therefore, AASCCIB B 1 was chosen to build the wall in the following modeling. The heat transfer coefficient and average von Mises equivalent stress of AASCCIB B 1 were 0.446 W/(m 2 ·K) and 9.52 MPa, respectively.

Improvement Effect of AASCCIB Wall on the Indoor Thermal Environment
Through the ANSYS CFX simulation, the indoor temperature changes of the building with the OB wall and building with the AASCCIB wall were informed. China is located in the northern hemisphere of the Earth. The day with the shortest sunshine time and the lowest solar altitude angle is the winter solstice. From the perspective of passive heat collection, the amount of solar radiation on this day is the lowest of the year. On the contrary, the summer solstice is the day with the longest sunshine time and the highest solar altitude angle, i.e., the solar radiation on this day is the highest in the year. Therefore, in the simulation analysis of the thermal insulation performance of the wall, the winter solstice (summer solstice) was selected as the date of the outdoor boundary condition. In addition, the indoor temperature of the room with the worst thermal comfort on the first floor (second floor) was used as the temperature of the first floor (second floor) in this study. The simulation results are presented in Figure 15. As shown in Figure 15a, the indoor temperatures of the rooms on the first and second floors of the building with the OB wall were lowest at 6:00 and highest at 14:00. The indoor lowest, highest, and average temperatures of the room on the first floor were 3.10, 5.65, and 4.47 • C, respectively. The indoor lowest, highest, and average temperatures of the room on the second floor were 3.37, 5.63, and 4.61 • C, respectively. For the building with the AASCCIB wall, the indoor temperatures of the rooms on the first and second floors were lowest at 6:00 and highest at 14:00. The indoor lowest, highest, and average temperatures of the room on the first floor were 4.49, 6.21, and 5.31 • C, respectively. The indoor lowest, highest, and average temperatures of the room on the second floor were 4.83, 6.10, and 5.43 • C, respectively. Compared with the indoor lowest and average temperatures of the building with the OB wall, those of the building with AASCCIB walls increased by at least 1.39 and 0.82 • C. In addition, the indoor temperature difference decreased by at least 0.83 • C. Figure 15b shows that the indoor temperatures of the rooms on the first and second floors of the building with the OB wall were lowest at 5:00, 24.06, and 24.64 • C, respectively. The indoor temperatures of the rooms on the first and second floors were highest at 17:00, 29.89, and 29.81 • C, respectively. The average indoor temperatures of the rooms on the first and second floors were 27.37 and 27.54 • C, respectively. For the building with the AASCCIB wall, the indoor temperatures of the rooms on the first and second floors were lowest at 5:00, 24.11, and 24.75 • C, respectively. The indoor temperatures of the rooms on the first and second floors were highest at 15:00, 28.14, and 28.03 • C, respectively. The average indoor temperatures of the rooms on the first and second floors were 26.58 and 26.70 • C, respectively. Compared with the indoor highest temperature of the building with the OB wall, that of the building with the AASCCIB wall decreased by at least 1.75 • C. The indoor average temperature and temperature difference decreased by at least 0.79 and 1.89 • C, respectively.
In summary, the indoor thermal environment of the building with AASCCIB walls was significantly improved compared to that of the building with the OB wall. This was manifested by the increase in the indoor lowest and average temperatures on the winter solstice. The indoor temperature difference decreased. On the summer solstice, the indoor highest temperature, average temperature. and temperature difference decreased.

Conclusions and Outlook
In this study, six types of AASCCIBs with different internal structures were designed. Based on the FEM, the thermal and mechanical performances of the six AASCCIBs were simulated using ANSYS CFX. The AASCCIB with the optimal comprehensive performance was determined through multi-objective optimization. The improvement effect of the AASCCIB wall on the indoor thermal environment relative to the OB wall was quantitatively analyzed using ANSYS CFX. The conclusions of this study can be summarized as follows: 1.
The von Mises equivalent stress and heat transfer coefficient of the AASCCIB decreased with the increase in hole ratio when the hole shape and number of hole rows were constant.

2.
AASCCIB B 1 had the optimal comprehensive performance among the six AASCCIBs. The heat transfer coefficient and average von Mises equivalent stress of AASCCIB B 1 were 0.446 W/(m 2 ·K) and 9.52 MPa, respectively.

3.
Compared with the indoor lowest and average temperatures of the building with the OB wall, those of the building with the AASCCIB wall increased by at least 1.39 and 0.82 • C on the winter solstice, respectively. The indoor temperature difference decreased by at least 0.83 • C. In addition, the indoor highest temperature, average temperature, and temperature difference decreased by at least 1.75, 0.79, and 1.89 • C on the summer solstice, respectively.
This study quantitatively elucidated the improvement effect of the AASCCIB walls on the indoor thermal environment relative to the OB walls. However, the positive effect of the AASCCIB walls on the building energy consumption was unclear. To promote the development of wall materials and building energy efficiency, the reduction effect of the AASCCIB wall on the building energy consumption will be quantitatively analyzed in a follow-up study.
Author Contributions: Funding acquisition, Y.Z.; Methodology, X.F., Y.G. and Q.Z.; Software, X.F., Y.G. and Q.Z.; Writing-original draft, X.F. and Y.G.; Data curation, X.F. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.