#### Abstract

The integrated thermal protection system (ITPS) is a complicated system that addresses both mechanical and thermal considerations. An M-pattern folded core sandwich panel packed with low-density insulation material provides inherently low mass for a potential ITPS panel. Herein, we identify the most influential geometric parameters and establish a viable, computationally efficient optimization procedure. Variables considered for optimization are geometric dimensions of the ITPS, while temperature and deflection are taken as constraints. A one-dimensional (1D) thermal model based on a modified form of the rule of mixtures was established, while a three-dimensional (3D) model was adopted for linear static analyses. Parametric models were generated to facilitate a design of experiment (DOE) study, and approximate models using radial basis functions were obtained to carry out the optimization process. Sensitivity studies were first conducted to investigate the effect of geometric parameters on the ITPS responses. Then optimizations were performed for both thermal and thermal-mechanical constraints. The results show that the simplified 1D thermal model is able to predict temperature through the ITPS thickness satisfactorily. The combined optimization strategy evidently improves the computational efficiency of the design process showing it can be used for initial design of folded core ITPS.

#### 1. Introduction

Space vehicles encounter severe aerodynamic heating and pressure loads during ascent and reentry. Thermal protection systems (TPS) are therefore required to maintain the underlying structural temperature within acceptable limits [1]. Conventional TPS are not load-bearing components, while an integrated thermal protection system (ITPS) is capable of bearing both mechanical and thermal loads [2, 3]. Both metal foam core and corrugated core have been widely investigated to achieve this goal [3–9]. Metal foam is an attractive proposition as a sandwich core for its desirable combination of mechanical, thermal, acoustic, and impact-absorbing properties [10]. Venkataraman et al. [4] optimized the density distribution to minimize the heat conduction through the metal foam thickness. Zhu et al. [5] evaluated the performance of the metal foam-based ITPS and compared it with the conventional TPS design. Furthermore, an open-cell metal foam core also allows for active cooling of the sandwich structure. Rakow and Waas [2] studied the thermal-mechanical response of actively cooled metal foam sandwich panels for a thermal protection system. It was shown that both temperature and thermal-mechanical deformation were reduced via active cooling. Corrugated core sandwich structures have been proposed for ITPS due to the favourable strength-to-weight ratio and good manufacturability. Bapanapalli et al. [3] studied the corrugated core sandwich panel concept and established an optimization procedure to minimize the weight of the ITPS panel. Sharma et al. [6] homogenized the ITPS panel based on finite element models and constructed multifidelity models to reduce the cost during optimization for a given accuracy. Xie et al. [7, 11] developed a two-step optimization approach to obtain the optimal ITPS dimensions.

Another promising structure suitable for the ITPS concept is a folded core sandwich panel using a three-dimensional core which can be folded from a planar base material along certain geometric lines [12, 13]. It can be made from various sheet materials in a convenient continuous process and can be tailored into different geometries meeting different functional requirements. Besides, the open-cell design of folded core eliminates corrosion issues which is a common problem for metallic honeycomb core structures, especially for those used in the aerospace industry [14].

These advantages have led to several research programmes into folded core structure [15–20]. Wang and Xu [15] experimentally investigated the soundproof characteristics of folded core sandwich plates. Heimbs et al. [16, 17] evaluated the mechanical behaviour of folded structures under compression, shear, and impact loads both experimentally and numerically in a series of studies. Sab and Lebée [18, 19] estimated the bounds for the transverse shear modulus of a folded core by using a unit load method, and later Bending-Gradient plate theory was applied to assess its shear stiffness. Currently, however, most research on folded core structures focuses on its acoustic, mechanical, and impact properties, while little work exists on its capacity for potential thermal protection.

The present study establishes a thermal-mechanical optimization procedure for a folded core sandwich-based ITPS. First, a novel analytical formulation for the effective thermal conductivity is established using a modified rule of mixtures approach, and a comparative study between the homogenized one-dimensional 1-D) and corresponding three-dimensional 3-D) thermal models was performed. Then the simplified 1D thermal model is used to determine the temperature distribution, which is taken as the thermal load in the subsequent linear static analyses of the folded core sandwich panel. Sensitivity studies are then conducted to identify the most influential parameters. Finally, the optimization approach is presented and used to obtain the optimal results of the ITPS.

#### 2. Folded Core TPS

##### 2.1. Geometric Description

The integrated thermal protection system (ITPS) is a sandwich panel with two thin face sheets separated by folded cores, which can be made of homogeneous, isotropic materials such as metals or orthotropic materials such as composite laminates, as shown in Figure 1. The empty space in the folded core is packed with a non-load-bearing insulation material such as SAFFIL [3]. As a single component, this multifunctional structure provides thermal protection as well as load-bearing capabilities.

Typically, a unit cell of M-pattern folded core can be completely described by a set of five independent geometric parameters as shown in Figure 17(a). Using terms in [21], the geometric model is defined by the folded core height* H*, half of the zigzag lines pitch* S*, the saw-tooth lines pitch 2*L*, the zigzag lines amplitude* V*, and spacing between the saw-tooth lines* B*.

##### 2.2. Material Properties

Material selection for different parts of the sandwich panel is based upon previous experience with ARMOR TPS [22] and the corrugated core ITPS [3]. Inconel 718, Ti-6Al-4V titanium alloy, and aluminium 2024 alloy are used for the top-face sheet, folded core, and the bottom-face sheet respectively. SAFFIL foam with a density of 48 kg/m^{3} is used for insulation.

Temperature dependent thermal properties are used in our heat transfer analyses [23, 24], while physical properties are considered constant [7] in linear static and buckling analyses to further simplify the model in this preliminary study. The heat transfer mechanisms through the fibrous insulation involve solid conduction, gas conduction, radiation, and possible natural convection. For simplicity, effective thermal conductivity is used to estimate the combined heat transfer mechanisms in the SAFFIL insulation [25]. All used material properties are found in [23–25].

#### 3. Heat Transfer Analysis

##### 3.1. Thermal Load and Boundary Conditions

Incident heat flux on the vehicle varies with many factors, such as the shape of the vehicle, the trajectory profiles, and different locations on the vehicle. For the current study, incident heat flux during reentry of a reusable launch vehicle (RLV) is used [3]. The outer surface of the top-face sheet is subjected to a transient heat flux as shown in Figure 2. The inner surface of the bottom-face sheet is conservatively assumed to be adiabatic. The initial structural temperature is assumed to be 295 K.

Radiation to the ambient atmosphere is also considered with a constant emissivity of 0.86, which is typical for TPS exterior material and surface treatments. Convection is modelled after the vehicle touches down on the ground. The convective heat transfer coefficient is assumed to be 6.5 W/m^{2}·K [3], and the ambient temperature is assumed to be 295 K.

##### 3.2. Simplified 1D Thermal Model

In order to study the structural performance under temperature environments, it is necessary to obtain the temperature distribution through the panel. Miller and Weaver [26, 27] developed analytical models based on techniques of integral transforms and separation of variables and described the transient temperature distribution through an air-filled box structure and a multilayered plate, respectively. Finite element method (FEM) models were also established for comparison purposes and showed a good fit. In this preliminary study, FEM models are currently adopted using Abaqus [28] to determine the temperature distribution. Due to the periodic characteristic of folded core, it is assumed that there is no lateral heat flow out of the unit cell, which means the heat flux incident on a unit cell only affects this specific unit cell. Although the temperature distribution could possibly be marginally different in a real ITPS panel due to the fact that the ITPS panel is mounted on the stringers and frames which could cause potential lateral heat flow across units, this secondary effect is not taken into consideration in this preliminary design. According to studies by Blosser et al., 1D models can reasonably predict temperatures in TPS systems [25, 29]. In order to reduce the time cost of the transient heat transfer analyses, the folded core ITPS thermal model is further simplified to a 1D model.

In the following the homogenized properties of the folded core are calculated using a modified rule of mixtures approach. The homogenized results for density , specific heat , and thermal conductivity are where is the volume fraction; the subscripts , , and refer to web, SAFFIL, and homogenized core, respectively; the superscripts and represent the parallelepiped part and cuboid part of the core along the direction (see Figures 17(a) and 18). Specifically, (3) shows that thermal conductivity is enhanced by considering the irregular shape of the core. The terms and are inclined angles and is the included angle, as shown in Figure 18. The detailed derivation of the homogenized properties is given in the appendix.

In the case of transient heat transfer analyses, we compare results from the 3D and 1D thermal models. The 8-node brick element (DC3D8) and 2-node link element (DC1D2) were used for the 3D and 1D models, respectively. Mesh convergence studies (Figure 3) for the 3D model show how temperature changes with the number of elements at the midpoint of the bottom-face sheet.

The mesh size with 6195 nodes and 4900 elements was eventually settled upon for the 3D heat transfer model as its temperature results differed only marginally from a finer mesh. Similar studies were also carried out for the 1D model, resulting in 21 nodes and 20 linear link elements of type DC1D2 being used. Figure 4 shows the mesh for the 3D and 1D heat transfer models separately, where Figure 4(a) is the mesh with SAFFIL, while Figure 4(b) refers to the solid portion only, that is, face sheets and web. A single run for each model shows that the 1D model reduces the computational time by nearly 70 percent. Several key locations on the folded core panel are monitored through analyses and are schematically illustrated in Figure 5.

**(a)**

**(b)**

**(c)**

The transient thermal analyses results for a folded core sandwich panel are shown in Figures 6–8, where temperatures versus reentry times at three positions (top, mid, and bottom) are illustrated. Figure 6 shows that on the top-face sheet the 1D model can predict temperature accurately, which is approximately the average value of the four locations’ temperatures in the 3D model. Figure 7 shows that the 1D temperatures are only slightly higher than the 3D model results on the bottom-face sheet, thus showing good validation of the simpler 1D model. Figure 8 shows the temperature variation with time at the mid position of the core. The maximum deviation between the 1D and 3D models occurs when the heat flux is relatively high. This phenomenon is expected to occur due to the existence of the inclined angles. As the 1D prediction lies mostly in the middle of the 3D results and because it reflects the changing tendency of the temperature correctly, it is considered acceptable for use in the preliminary design process. Based upon previous analyses, it can be concluded that the 1D heat transfer model predicts temperature with satisfactory accuracy through the ITPS thickness.

During the transient heat transfer analysis, the maximum temperature of the bottom-face sheet is obtained. At the same time, the temperature distribution through the panel thickness is also extracted for subsequent use in the linear static analyses described in the next section, achieved by using python code from the FEM results file.

#### 4. Linear Static Analyses

In this section, the 3D model used for the linear static analysis is described. As the mechanical properties of the SAFFIL material are negligible when compared to the face sheets and the web, the soft insulation material is not taken into account during structural analyses. In this current study, four unit cells were considered in the folded core panel. To decrease the computational expense, only half of the panel was modelled in the FE linear static analysis. The 3D folded core panel is modelled using the eight-node shell element (S8R). Mesh convergence studies were conducted, with approximately 9715 nodes and 3110 elements used in the linear static model as shown in Figure 9.

The connections between face sheets and the core are accomplished by surface-based tie constraints. Temperatures and aerodynamic pressure loads are applied. Temperature profile data are taken from the 1D heat transfer analysis described in the previous section at the time when the top-face sheet reaches its maximum temperature throughout the whole reentry stage. Then smooth temperature distributions are obtained by fitting a quartic polynomial to the actual temperature data. These polynomials are used to assign temperatures to corresponding linear static models. In this way, we assume that the panel temperatures are uniform throughout their respective lengths and widths, a necessary and limiting consequence of using a 1D thermal model. By considering the temperature variation throughout the panel thickness, this approximation is deemed to be acceptable.

A uniform pressure load of 101325 Pa, which is considered as the maximum pressure that the ITPS may encounter during reentry [3], is applied on the inner surface of the top-face sheet. The model is shown in Figure 10. Symmetric boundary conditions (-axis) are applied to the symmetric edges. For the other three edges, which are the actual panel edges, the bottom-face sheet edges were fixed in vertical displacement (-axis) and all three rotations, and the top-face sheet edges were fixed in all three rotations while allowing all three displacements to occur. In addition, a fixed displacement is applied at the intersection point between two actual edges of the bottom-face sheet to eliminate the rigid body modes.

#### 5. Optimization Study

##### 5.1. Optimization Approach

A combined optimization approach is adopted as shown in Figure 11. Both of the thermal and structural models are generated parametrically in order to carry out the sensitivity study and optimization procedure efficiently. In order to reduce the computational time, simpler, approximate models were established for both objective function and constraints. This was done using radial basis functions (RBF) [30], which are a type of neural network employing a hidden layer of radial units and an output layer of linear units. This neural network is suitable for approximating a wide range of nonlinear spaces.

First, a Latin Hypercube Sampling (LHS) [31] technique was used for the design of experiments (DOE). Variables are normally referred to as factors in a DOE study, while values of variables are known as levels. With the LHS technique, the design space for each factor is uniformly divided, and then these levels are randomly combined to specify sample points defining the design matrix. It provides an efficient method for generating random sample points, which also have a uniform distribution over the entire design space. Then heat transfer and linear static analyses were conducted, respectively, for these samples and corresponding responses obtained.

Afterwards, RBF approximate models were generated based on these data. A set of random points in the complete design space validate its accuracy. These approximate models were continuously updated with additional sample points until necessary accuracy is obtained, which were then used to replace FEM models in the following optimization process.

During optimization, the multi-island genetic algorithm (MIGA) [32] was employed first to locate the approximate global optimal solution. Genetic algorithms (GA) treat complex nonlinear optimizations relatively well. MIGA, which is a further development of GA, divides each population of individuals into several subpopulations called islands, and traditional genetic operations are performed on each island separately. Several individuals are then selected from each island and migrated to different islands periodically. The migration operation maintains the diversity of probable solutions and prevents the premature phenomena. Finally, based on the MIGA solution, sequential quadratic programming was applied to further find the optimum design.

##### 5.2. Sensitivity Study

Sensitivity studies were conducted based on the homogenized 1D transient heat transfer analysis to investigate the contribution of different geometric parameters to the ITPS responses.

Isight [33] is employed to facilitate the conducting of the DOE by Abaqus. Altogether 100 sample points were studied, and then the most influential parameters were obtained through sensitivity analyses. The sensitivity analyses results for maximum bottom-face sheet temperature, areal density of the panel, and maximum top-face sheet deflection are presented using bar graphs in Figures 12–14, blue for positive effect and red for negative effect.

The sensitivity bar graphs plotted in Figure 12 indicate that the folded core height plays a leading role in the inner temperature control. A larger height means thicker insulation material, which blocks most of the heat flux that flows into the panel.* L* also has a remarkable negative effect on the bottom-face sheet temperature; with the increasing, the unit cell density decreases, which leads to a longer heat transfer path. The bottom-face sheet acts as a heat sink component, so thickness of bottom-face sheet is also an influential parameter. A thicker web means more heat is conducted to the bottom-face sheet, so the web thickness adversely affects the bottom-face sheet temperature.

Figure 13 shows the sensitivity indices for the areal density. It is seen that the thickness of the core and face sheets together with the folded core height contributes most to the areal density. In contrast, a larger and mean a larger panel area, which results in a minor areal density.

From Figure 14, the panel dimension and the top-face sheet thickness are seen to be relatively important factors for top-face sheet deflection. Due to the fact that the top-face sheet temperature is close to the radiation equilibrium temperature, the temperature distribution through the panel thickness has little effect on the top-face sheet deflection. Such considerations explain why the sensitivity indices of the top-face sheet deflection with respect to* H*,* t*, , and are relatively low, although they greatly affect the maximum bottom temperature, as shown in Figure 12.

##### 5.3. Thermal Optimization

In this section, only the thermal model is taken into consideration. According to our previous sensitivity study, the spacing has little effect on bottom-face sheet temperature, areal density, and top-face sheet deflection. It is fixed here ( mm) in order to reduce the number of design variables. The ranges of the design variables are the same as those used in the sensitivity study as shown in Table 1.

The mass per unit area of the panel, used as the objective function, is simply given by

The maximum bottom-face sheet temperature, which must be below a certain temperature limit, is taken as the constraint.

The optimization problem is succinctly described by

First, 300 points were selected in the design space using LHS and then used to construct the approximate model for both objective and constraint functions. Another 50 random points were used to evaluate the accuracy of the approximate model. Details are shown in Table 2, where RMSE stands for root mean square error and stands for coefficient of determination. It is shown that the approximations for maximum bottom-face sheet temperature and areal density are of high quality.

Afterwards, the approximate model was used to replace the previous FEM thermal model to carry out the optimization process. The critical parameters of MIGA are as follows: the subpopulation size is 20, the number of islands is three, the number of generations is 50, the rates of crossover and mutation are 1.0 and 0.01, the rate of migration is 0.1, and the migration interval is five.

Three cases are considered, and the optimal results are shown in Table 3. It is seen that and are at their upper bounds to maximize the area of the panel.* V*,* t*, and are at their lower bounds to minimize the proportion of the core web and top-face sheet. The optimal results show the same trends as the sensitivity analysis results. When the ITPS has a less restrictive temperature requirement, is at its lower bound, while increases as the bottom-face sheet temperature limit reduces. As the bottom-face sheet acts as a heat sink component, increases to contain extra heat once reaches its upper bound.

##### 5.4. Thermal-Mechanical Optimization

A large top-face sheet deflection could change the flow around the ITPS, which can lead to a notable increase in aerodynamic heating. On the basis of our previous thermal optimization, a linear static model is now taken into our consideration. The design space and objective function are the same as those used during the thermal optimization, while maximum top-face sheet deflection constraint is imposed as well as the bottom-face sheet temperature limits. Linear static analyses show that extremely high stresses occur locally due to the boundary effects, while the stress of sections away from the boundary lies far within the material strength. Hence, stress is not considered as a constraint.

The optimization problem is simply stated as follows:where refers to the maximum top-face sheet deflection and is the diagonal length of the panel. Here, 400 LHS-obtained points were used first to establish the RBF approximation, whilst 50 random points were used to verify its accuracy. In order to achieve an acceptable error level, 200 additional sample points were sequentially added to the design space. Accuracy evaluation results are shown in Table 4.

Finally, thermal-mechanical optimization was conducted based on this approximate model. The parameter values used in MIGA are the same as those in last section. Optimal results are shown in Table 5. It is seen that* H*, and are mainly influenced by the temperature constraint, similarly to the tendency reflected in the thermal optimization. Top-face sheet thickness is no longer at its lower bound and is dominated by the top-face sheet deflection and also varies slightly with the temperature constraint. From cases –, increases to meet a more challenging deflection constraint, which also means that more heat is conducted into the panel; in order to maintain the bottom-face sheet temperature, increases along with .

The optimized variables for cases and were subsequently input into the finite element analyses. A comparison of the RBF predictions and FEM results is shown in Table 6. The relative error between the RBF and FEM values is less than 5%, which indicates that the approximate model has fairly good accuracy.

Figure 15 illustrates the displacement contour of the ITPS panel in the direction for case . It is seen that the local effects of the pressure loading only appear on the top-face sheet and the maximum deflection occurs at the middle of the top-face sheet.

In the above optimization process, temperature and displacement are taken as constraints, while buckling, which is also likely to cause catastrophic failure, is not considered. Therefore, additional buckling analyses are conducted for case at three different times, that is, times (Time 1, Time 2, and Time 3) when top-face sheet, web middle, and bottom-face sheet reach their respective maximum temperature. The same model as the linear static analysis is used in the buckling analysis. The minimum positive eigenvalues are listed in Table 7 and the corresponding buckling models are shown in Figure 16. It is seen that buckling would probably not occur for the current design of case .

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

From the temperature variations in Figures 6–8, the times when the top-face sheet and web middle reach their maximum temperature are quite close, which leads to a similar buckling mode as shown in Figures 16(a) and 16(b). When the bottom-face sheet temperature reaches its maximum, the temperature distribution through the whole panel thickness is significantly changed compared with the other two. Hence, the web buckles in a different position at Time 3 as shown in Figure 16(c).

#### 6. Conclusions

A novel thermal-mechanical optimization procedure has been developed for the folded core-based integrated thermal protection system (ITPS) concept. In addition, an improved analytical rule of mixtures approach for calculating thermal conductivity considering the inclined shape of the folded core has been presented and a thermal model that homogenizes variations through the structure was developed to reduce the computational cost associated with a 3D structure. Then sensitivity analyses as well as optimization studies were performed to identify key variables of influence. Key conclusions are considered below.

The comparative thermal analyses show that the 1D model has a relatively high accuracy in comparison with 3D models for heat transfer analysis. This approach reduces the computational time by nearly 70 percent, while maintaining sufficiently accurate prediction of temperature distribution. For current material combinations, the bottom-face sheet temperature decreases significantly with an increasing folded core height, saw-tooth lines pitch, and bottom-face sheet thickness. The ideal web thickness is found, not surprisingly, to be as thin as possible in order to reduce the mass and the amount of heat entering the panel. The spacing between the saw-tooth lines has little effect on the responses considered in this paper. The radial basis function approximations are shown to be of high quality and can be used instead of more expensive FEM models to perform optimization studies.

In the current study, only pressure load and thermal load have been considered; however, the ITPS structure is subjected to various kinds of loads during reentry. Buckling constraints and the influence of the folded core topology should be considered in future work.

#### Appendix

#### Derivation of Homogenized Properties of M-Pattern Folded Core

Consider a typical M-pattern folded core unit cell as shown in Figure 17(a). It can be folded continuously from a flat sheet consisting of a repetitive combination of parallelogram and rectangular plates as shown in Figure 17(b). The parameters of the flat sheet can be derived from those of the corresponding folded core as [13, 21]

The included angle , which could vary in the range of 0° to 90°, is related to independent parameters by

As is more frequently used in the following derivations, the relation between the included angle and independent parameters is rearranged as

The effective thermal conductivity of the folded core along the core thickness direction (-direction) can be derived by analysing heat conduction through the half unit cell shown in Figure 18. By considering the inclination of the webs, thermal conductivity along the thickness direction could not be simply homogenized using the general rule of mixtures. Amendments are achieved by introducing inclined angles and for the parallelepiped and cuboid part, respectively. For simplicity, both panels are assumed to be free of process-induced geometrical imperfections and have uniform web thickness.

For the parallelepiped part (see Figure 19(a)), that is, part, the effective thermal conductivity is given bywhere and are the volume fraction of the web and SAFFIL:

**(a)**

**(b)**

Similarly, for the cuboid part (see Figure 19(b)),that is, part, where and are the volume fraction of the web and SAFFIL:

Once the effective thermal conductivity for each part is obtained, the general rule of mixtures is used. The homogenized thermal conductivity of the proposed folded core can then be expressed aswhere and refer to the volume fraction of the parallelepiped and cuboid part, respectively.

By invoking (A.4) and (A.6), (A.8) can be rewritten as

The effective density and specific heat of the homogenized folded core is readily calculated by imposing the general rule of mixtures:where and are the volume fraction of the web and SAFFIL for the folded core unit cell:

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was supported by Funding of Jiangsu Innovation Program for Graduate Education [Grant no. CXLX13_163], the Fundamental Research Funds for the Central Universities [Grant no. NZ2016101], and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). Chen Zhou would also like to thank China Scholarship Council (CSC) for the funding of his research work through University of Bristol.