Next Article in Journal
Linear Contrails Detection, Tracking and Matching with Aircraft Using Geostationary Satellite and Air Traffic Data
Next Article in Special Issue
Dynamic Analysis of a Large Deployable Space Truss Structure Considering Semi-Rigid Joints
Previous Article in Journal
Velocity Mapping of an H2 − O2 Exhaust Jet in Air by Means of Schlieren Image Velocimetry (SIV)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stress Characteristics and Structural Optimization of Spacecraft Multilayer Insulation Components

1
Mechanical Electrical Engineering School, Beijing Information Science & Technology University, Beijing 100192, China
2
Department of Mechanical Engineering, State Key Laboratory of Tribology, Tsinghua University, Beijing 100084, China
3
Beijing Institute of Spacecraft Environment Engineering, Beijing Key Laboratory of Environment & Reliability Test Technology for Aerospace Mechanical & Electrical Products, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(7), 577; https://doi.org/10.3390/aerospace10070577
Submission received: 20 May 2023 / Revised: 20 June 2023 / Accepted: 21 June 2023 / Published: 21 June 2023

Abstract

:
Multilayer insulation (MLI) components are important parts of spacecraft, in which thin films are the core elements. The film thicknesses are generally small, between 10 and 30 μm, to reduce the weight of the spacecraft. During the launch of a spacecraft, there is a rapid drop in the internal pressure, which causes the internal gas to flow out rapidly through the component. The resulting fluid force may cause film damage and failure, thereby directly affecting the normal operation of the spacecraft. Therefore, the mechanical characteristics of the thin films under rapid decompression conditions were investigated during this study. Considering the effects of the flow field stress distributions on the films during the rapid decompression, a fluid–structure interaction (FSI) model for a component was first constructed. The results show that the stress is largest in the outlet film, which is the part of the overall structure most vulnerable to failure. Furthermore, the effects of the structural parameters of the component on the stresses in the different film layers were analyzed using the orthogonal experimental method. The results show that the film thickness had the largest influence, followed by the film hole diameter, the number of component layers, and finally the staggered hole distance. Finally, a structurally optimized design scheme for the component is proposed based on a parameter range analysis with respect to the maximum film stress. After optimization, the maximum stress of the thin film decreased by 97.6%. This research has practical engineering value for the structural design and optimization of MLI components.

1. Introduction

Multilayer insulation (MLI) components are the key pieces of equipment installed on the outer surfaces of spacecraft. They are used for heat insulation and rapid decompression. Failure of these components can seriously affect the normal operation of the spacecraft [1,2,3,4]. These components are subjected to gas friction during spacecraft launching. Additionally, the external pressure of a spacecraft rapidly drops during launching, which causes the internal gas to rapidly flow out through the component due to the rapid drop in the external ambient pressure. These components are also simultaneously affected by random spacecraft vibrations, resulting in the vibration of their structures [5,6,7]. The components are composed of nylon nets and stacks of thin insulation films. The material properties, film density, and other structural parameters of the components can affect the thermal performance of the components, and improving the thermal insulation properties of the components is an important research focus in this field [8,9,10,11]. Therefore, the structural design of components operating in complex conditions is particularly important.
Bapat et al. [12] conducted an experimental investigation using different combinations of spacer and shield materials and compared the results with those of a numerical model. The results showed that the combination of an aluminized 12 μm thick Mylar layer and a 76.2 μm thick glass fabric layer was the best for MLI applications. Wei et al. [13] experimentally investigated the structural and shape effects of a perforated MLI blanket. A new calorimeter with spherical top and bottom surfaces was compared with a cylindrical calorimeter; the results indicated that 13% excess heat flux was transferred in the new calorimeter. This occurred because of the “edge effect” formed by the MLI deformation. Thus, it was concluded that the shape and structure affect MLI performance. Zhou et al. [14] studied the effects of the interlayer vacuum degree, the interlayer compaction pressure, and the layer density on the performance of a high-vacuum multilayer insulation component.
However, the structural parameters of a component impact not only its thermal performance but also its internal flow field. Emphasizing the thermal performance of a component while ignoring its mechanical properties is not conducive to obtaining an integral optimized design for the component. The first automated transfer vehicle (ATV), which was named Jules Verne after the French writer, was launched on 9 March 2008. Immediately after the launch, direct observations of the Jules Verne ATV using cameras onboard the international space station (ISS) indicated that some of its MLI components (primarily those located on the pressurized cylinder) appeared to have partially detached from the structure, as shown in Figure 1.
Ren et al. [16] conducted research on the failure mode of an MLI component in a low-altitude environment. The results showed that a strong gas flow shearing force in the powered phase during the launch and a bulging force caused by decompression after entering the vacuum state were the primary factors that led to component damage. The results also indicated that considering the thermal performance of the component while ignoring the breathability could easily cause damage to the component. Woodward et al. [17] designed and tested a multilayer thermal insulation system that also provides debris and micrometeorite damage detection. Li et al. [18] studied the transient pressure on each film of a spacecraft during rapid decompression by establishing a computational fluid dynamics (CFD) model of the MLI component. The environmental adaptability and reliability of the component were evaluated, and it was revealed that the component and brush seal structures had some common points. Huang et al. [19,20,21] adopted a three-dimensional slice model to solve for the complex structures associated with multiple brush wires using CFD calculations. These studies provided important reference values for the study discussed in this paper. To date, relatively little research has been performed regarding stress analyses of MLI components, but other fluid–structure interaction (FSI) calculation methods adopted under similar working conditions were used as references when conducting this research. Boustani et al. [22] simulated the descent process of the ASPIRE SR01 supersonic parachute using an FSI model. Yang et al. [23] also solved for the flow process of a flexible buffer liquid inside a solid rocket motor using an FSI model.
With the rapid development of the finite element technique, the FSI calculation method has been applied in an increasing number of fields [24,25,26]. The FSI method can capture the impact of flow field stresses on solid deformation more accurately than can the CFD method, and it can perform the complex calculations associated with the impact of a solid deformation on the flow field.
The MLI component of a spacecraft deforms due to the gas flow during the rapid decompression process. However, to reduce the weight of the spacecraft, the MLI film thickness is usually only a few microns to tens of microns, thereby causing the films to be prone to damage and failure when they are subjected to strong forces. To prevent damage to the MLI component during rapid liftoff, a combination of nylon nets and thin films is used in the structural design of the MLI component, where the nylon nets prevent the film from producing large deformations and thus being destroyed, and the hole design in the film quickly releases the internal fluid pressure to prevent the film from producing large deformations. However, most existing studies focus on the thermal resistance effect of MLI components, while neglecting to explore the mechanical properties, especially the size of the holes, the distribution of the holes, and the thickness of the thin film. Therefore, the structural design of this component is a problem that cannot be ignored. The stress problem associated with spacecraft MLI components was investigated during this study. The internal flow field and film stress were first studied using an FSI model under different structural parameter values. Then, the internal film stress characteristics were further explored. The results of this study provide a basis for spacecraft MLI component structural designs.

2. FSI Model of an MLI Component

2.1. FSI Model

An MLI component of a spacecraft has a three-dimensional structure similar to that of a plane, and it has a large developed area. It has a multilayer structure formed by stacking nylon nets and thin films. The films are made of aluminized polyimide materials with good thermal insulation properties. The film thicknesses are only 10–30 μm and their lengths and widths are generally a few hundred millimeters. The films have holes to allow gas to flow through them. The space between two film layers is filled with a nylon net with a thickness of approximately 200 μm, and gas can flow through the nylon net. The component is composed of nylon nets and thin films, and the gas flow through the nylon nets cannot be ignored. Therefore, significant mesh generation and calculation challenges arise when modeling an MLI component.
To solve these problems, a typically sized FSI model was established for the component during this study so FSI calculations could be conducted. The five-unit MLI component (which included five nylon nets and six films) shown in Figure 2 was taken as an example.
The FSI component model developed during this study consists of a fluid model and a solid model. The solid model was formed by stacking nylon nets and thin films, while the gaps between the nylon nets and the thin films represent the fluid domain. The nylon nets were modeled using a single strand of nylon rope, and circular nodes were adopted at the intersections to facilitate mesh division and improve the mesh quality. The holes in the films are the key channels that connect the fluid domains of the adjacent nylon nets. The holes in adjacent film layers are staggered, which means that the holes in adjacent film layers have different positions.
During the process of assembling the nylon nets and films, there was randomness in the placement of the film holes and the nylon nets. Therefore, some reasonable simplifying assumptions were made for this model when analyzing the impacts of different structural parameters on the film stress [19,20,21]:
(1) There is only one film hole in each diamond-shaped hole of a nylon net, and it is positioned in the middle of the diamond structure. This assumption was made to study the effect of the distance between the staggered holes of two adjacent film layers on the flow field.
(2) Each film layer has only one film hole. There is a certain distance between the staggered holes in two adjacent film layers, and the holes are arranged periodically. This means that the hole position in each odd film is the same and that the hole position in each even film is also the same.

2.2. Calculation Method

FSI calculations were performed for the five-unit MLI component described above. The FSI simulation method fully considers the interactions between the fluid and the solid in complex models. The FSI model can apply the internal fluid pressure to the solid surface to obtain solid deformation and stress results that are more accurate than those produced by CFD calculations.
Due to the rapid drop in the external pressure and the rapid outflow of the internal gas, the films of the component deform as a result of the gas flow during launching. Based on this process, the calculation method shown in Figure 3 was used to determine the flow field changes and the film stresses in the component. The fluid model and the solid model were first established and their meshes were generated. The interface between the fluid model and the solid model was set as the FSI surface, which forms the key data transmission bridge for the FSI process. Next, the pressure distribution on the FSI surface was obtained from finite element calculations of the fluid model. Finally, static calculations were performed for the solid model. During the calculation process, the stress distribution and magnitude on each side of each film layer were imported from the CFD model results to obtain the deformation and stress results for the films subjected to the fluid pressure.

2.3. CFD Mathematical Model

The controlling equation for the internal flow field of the MLI component is described next. The gas volume force was ignored, and the fluid was considered to be a compressible gas. The continuity equation, the momentum equation, and the energy equation can be respectively expressed as follows:
ρ t + ρ v = 0 ,
t ρ v + ρ v v = p + τ ¯ + S M ,
t ρ h t o t p t + ρ v h t o t = λ T + v τ e + v S M + S E ,
where ρ represents the gas density, t is the time, v is the velocity vector, p is the pressure, λ is the effective thermal conductivity, τ represents the stress tensor, T is the temperature, htot is the total enthalpy, and SM and SE are the generalized source terms for the mass and thermal energy (internal energy), respectively.
τ can be expressed as follows:
τ = μ ( v + ( v ) T 2 3 δ v ) .
δ, which represents the Kronecker Delta function, can be defined as follows:
δ = 1 0 0 0 1 0 0 0 1 .
The relationship between htot and the static enthalpy, hstat, can be expressed as follows:
h t o t = h s t a t + 1 2 v v , h t a t = h s t a t ( p , T ) .
The equation of state for an ideal gas can be expressed as follows:
ρ = p / R T e m / M w ,
where Mw represents the molar mass constant of the gas, R is the ideal gas constant, and Tem is the room temperature, which is equal to 300 K.
The gas accelerates in each film hole of the component, and the entire flow process is an entropy increment process. The maximum Reynolds number of the gas exceeds 9000. Therefore, the k–ε turbulence model, the turbulent kinetic energy equation, k, and the energy dissipation equation, ε, can be used together as follows:
t ρ k + ρ k v = μ + μ t σ k k + G k + G b ρ ε Y M + S k ,
t ρ ε + ρ ε v = μ + μ t σ ε ε + C 1 ε ε k G k + C 3 ε G b C 2 ε ρ ε 2 k + S ε ,
where G k represents the turbulent kinetic energy generated by the average velocity gradient, Gb represents the turbulent kinetic energy generated by the buoyancy, YM represents the turbulence dissipation rate generated by the turbulent expansion in the compressible turbulence, and μ t = ρ C μ k 2 / ε is the turbulent viscosity. C1ε, C2ε, C3ε, and Cμ are empirical constants with values as listed in Table 1, σk and σε are the Prandtl constants of the turbulent kinetic energy equation, k, and the energy dissipation equation, ε, respectively, with values as listed in Table 1, and Sk and Sε are the source terms.
A transient solver and a compressible ideal gas were used in the numerical solution process. The SIMPLEC algorithm, which is based on the SIMPLE algorithm but adds a pressure correction relaxation factor, was adopted to achieve better convergence during the calculation process. In the numerical solution process, the convergence accuracy was set to 1 × 10−6 to ensure the stability of the results and the accuracy of the calculations.

2.4. Mesh Generation

The mesh generation was challenging because the fluid model structure had a small thickness value but large length and width values. To solve this problem, a polygonal mesh was used to mesh all the fluid domains, as shown in Figure 4; this polygonal mesh produced a better mesh quality with fewer cells. Establishing circular nodes at the nylon net intersection points enables a smoother fluid mesh transition, which improves the mesh quality and the computational accuracy.
Mesh independence verification was conducted for the fluid model, and mesh refinement was performed in some regions to reduce the computational error in key regions of the model. The mesh independence verification results are shown in Figure 5. The results show that the mesh quality tended to stabilize when there were more than 8 million cells (the minimum cell size was 0.02 mm), thereby indicating that the CFD model calculations had reached mesh independence. Based on this result, a minimum cell size of 0.02 mm was used for the CFD calculations in this study.

2.5. Boundary Conditions

During launching, the external pressure of the spacecraft gradually decreases. Therefore, the fluid model outlet boundary was set as a variable pressure outlet with a pressure that varied over time. The fluid model boundary conditions are shown in Figure 6. The model only had an outlet for gas discharge, and the inlet was the wall boundary. The model was surrounded by symmetry boundaries, and the initial pressure of each fluid domain was set to an absolute pressure of 0.1 MPa.
The relationship between the external pressure (the outlet pressure of the fluid model) and time during launching is shown in Figure 7. The outlet pressure gradually decreases with time, and the stress on each surface in the fluid model reaches a maximum value when the pressure decreases most rapidly. Therefore, the moment when the films are subjected to the maximum pressure was selected for the FSI analysis. This was completed to obtain the stress characteristics of the films under maximum fluid pressure conditions.
By solving the CFD model, an outlet leakage–time curve was obtained, as shown in Figure 8. Figure 8 shows that the outlet leakage reaches a maximum near the 58th second. Based on this result, it was assumed that the films in the CFD model experience the maximum fluid pressure at the 58th second. Therefore, the fluid pressure results of the FSI model were extracted at the 58th second for the subsequent solid finite element analyses.
In practice, buckles are used to fasten the MLI components to the exterior of a spacecraft. Therefore, fixed supports were added around the films and at the end points of the nylon nets in the solid model, as shown in Figure 9.

3. Results and Discussion

3.1. Fluid Model Result Analysis

The structural parameters of the MLI component, as well as the corresponding symbols used to represent them in this paper, include the film hole diameter (D), the staggered hole distance (L), the number of component layers (N), and the film thickness (c). The five-unit MLI component used as an example in this study had the following structural parameters: D = 1 mm, L = 4 mm, N = 5, and c = 0.01 mm. The structural model shows that a fluid domain existed between every two film layers. Due to the small fluid domain thickness, the drop in the fluid pressure in the thickness direction could be ignored; therefore, the pressure in the fluid domain between two films was the same throughout that domain. The pressure contours in the flow field on one side of each of the six film layers are shown in Figure 10, where Ni represents the film layer number as counted from the inlet. Figure 10 shows that the nylon net forms a barrier to the gas flow and that the fluid pressure exhibits significant stratification at the nylon net. Meanwhile, the gas is not evenly distributed throughout the fluid domain of the nylon net; this result is primarily reflected in the relatively high pressure at the inlet side of the film hole, as well as in the observation that the fluid pressure gradually decreases from the inlet side to the outlet side of the hole in each film layer.
The pressure distribution on each film layer can be obtained from Figure 10. The fluid pressure differences between the films can be obtained by calculating the differences between the average fluid pressures on the surfaces of both sides of each film. To further study the effects of the structural parameters on the flow field pressure differences between the films, the flow field pressure differences calculated by CFD models with different structural parameter values were obtained by the single-factor variable method, as shown in Figure 11. The figure shows that the fluid pressures on the different films are not equal and that the fluid pressure gradually increases from the first layer to the last layer. This phenomenon occurs because the gas inside the component can only be discharged through the outlet. A film that is closer to the outlet experiences less gas resistance, and thus, it has a larger internal pressure difference.
Figure 11 shows that the pressure difference of the last film layer is always the largest. The effects of the structural parameters on the flow field pressure difference are described in the following six points.
(1) The film hole diameter, the staggered hole distance, and the film layer sequence are all inversely proportional to the maximum pressure difference. The film thickness, however, is positively correlated with the maximum pressure difference.
(2) The parameters that affect the maximum film stress, in descending order according to their degrees of influence, are the film hole diameter, the number of film layers, the staggered hole distance, and the film thickness.
(3) Of all the parameters investigated, the film hole diameter has the greatest impact on the maximum pressure difference. In the component, the hole diameter directly affects the flow of the internal gas, and the throttling effect of the thin film hole is crucial to prevent leakage. A smaller hole diameter causes a greater resistance to the gas flowing through the film hole, and gas that flows more fully into the nylon net fluid domain between the films causes more fluid kinetic energy to be converted into pressure energy in the fluid domain.
(4) As the number of the film layers increases, the number of film-layer nylon net sequences increases. This increase results in an increase in the internal gas resistance and a decrease in the flow velocity, which enables better conversion of gas kinetic energy into pressure energy in the nylon net fluid domain. Therefore, the influence of the number of film layers on the maximum pressure difference is also significant.
(5) The staggered hole distance changes the pressure difference by changing the gas resistance in each nylon net layer, which has a limited impact on the pressure distribution.
(6) The influence of the film thickness on the maximum pressure difference is minimal. Because the film thickness is very small, the time it takes for the gas to flow through the film hole is short; therefore, there is not enough time for energy conversion or dissipation.
Therefore, through finite element analyses of components with different structural parameters, the impacts of the structural parameters on the pressure difference of each film layer were obtained. These results lay a foundation for the structural design of MLI components. The pressure difference acting on a film can intuitively reflect the mechanical characteristics of the film. However, the component is composed of stacked films and nylon nets, and the nylon nets also impact the mechanical characteristics of the films. Therefore, finite element calculations were also performed for the solid component model to further analyze the film stresses and deformations, as discussed in the following section.

3.2. Solid Model Result Analysis

The solid model was formed by stacking films and nylon nets while leaving gaps between each layer of films and nets during the initial assembly. When a film deforms due to the influence of the fluid pressure, it may come into contact with a nylon net. Therefore, the films and the nylon nets were set as frictional contacts in the finite element model. The five-unit MLI component model, which was used as an example, had solid model structural parameters of D = 1 mm, L = 4 mm, N = 5, and c = 0.01 mm. Deformation nephograms for the films were obtained by introducing the fluid pressure onto both sides of the films, as shown in Figure 12.
The results show that the maximum film deformations occur in the middle regions of the films due to the fixed constraints around the film edges. Additionally, the overall maximum deformation occurs in the sixth film layer. It can also be concluded that the deformations of the first three film layers are relatively small and that these films do not come into contact with the nylon nets; therefore, the nylon nets do not significantly deform. The fourth and fifth film layers, however, come into contact with the nylon nets, and the nylon nets do significantly deform; their deformation nephograms are shown in Figure 13. The maximum nylon net deformations occur in the middle regions of the nets (the reverse direction of the Z-axis in Figure 13 is the direction of the gas flow). Due to the small nylon net deformation values, the deformation results are amplified 50 times in Figure 13 for better observation. They indicate that the nylon nets experience Z-direction deformations in their middle regions, which are caused by the Z-direction deformations of the films and the contact between the films and the nets.
Simultaneous lateral displacements also occur inside the nylon nets; these are primarily reflected in the X-direction motion in Figure 13. This occurs because the gas enters a nylon net through the hole in the upper film layer and flows out through the hole in the next film layer. This flow process generates lateral gas pressure on the nylon nets, resulting in deformations. Because the contact stress is greater than the internal stress, the deformations of the nylon nets in the Z-direction are greater than those in the X-direction.
Figure 12 and Figure 13 show the solid deformation trends for the FSI model. However, the research presented above is still limited for systematic analyses of structural parameters. To further study the effects of the structural parameters on the maximum deformation of each film layer, the single-factor variable method was used with the FSI model to calculate the deformation for different structural parameter values. The maximum deformation of each film layer was thus obtained, as shown in Figure 14.
Figure 14a,b shows the effects of different film hole diameters and different staggered hole distances on the maximum deformation of each film layer, respectively. The results show that the deformations of the fifth film layer are the same despite changes in the parameter values. This result indicates that the film has come into contact with the nylon net, which limits further deformation of the film and reduces the risks of film damage and failure caused by excessive deformation. The sixth film layer exhibits the largest deformation because it is subjected to the maximum fluid pressure and there is no nylon net to limit its deformation. Figure 14c,d depicts the effects of the different numbers of film layers and different thicknesses on the film deformation, respectively. The results show that the largest deformation occurs in the film at the outlet. As the film thickness increases, the film strength significantly increases and the maximum film deformation significantly decreases.
The fluid pressure was introduced onto both sides of the solid film models to obtain the film stress distributions, and the results are shown in Figure 15. The maximum film stresses are generated at the edges of the films, indicating that the fluid pressure generated by the gas flowing out from between the films is the key cause of film failure.
Figure 16 depicts the impacts of the different structural parameters on the maximum stress in each film layer. The results show that the film at the outlet has the maximum stress value because this is where the maximum fluid pressure occurs. Therefore, the film at the outlet is the part of the MLI component that is most vulnerable to failure. Increasing the film thickness significantly increases the film strength and reduces the maximum film stress.
The effects of the structural parameters on the deformations and stress distributions of the films were obtained by analyzing the results of the FSI model. The single-factor variable analysis method revealed that the film thickness is the most critical factor affecting the maximum film stress; it indicated that increasing the film thickness can effectively reduce the film stress. The impacts of the structural parameters on the film stress of the component are complex, and a method of applying the rules discussed above when optimizing the design of a component has great significance. Therefore, an optimal design scheme for the component is proposed next.

3.3. Structure Optimization Analysis

The effects of the MLI component’s structural parameters on the maximum film stress are complex, and strict restrictions are imposed on the maximum weights of components used on spacecraft. Therefore, the film thickness should be as small as possible so that the weight requirements can be met. According to the weight limitations and the actual manufacturing scenarios of components, the design limitations were obtained for a certain type of component. The parameter values were selected based on these design limitations, as shown in Table 2.
The process used to analyze the parameter effects on the maximum film stress using an orthogonal experimental method is described next. First, the correct orthogonal analysis table was selected according to the number of factors to be analyzed and the typical parameters to be calculated. Second, through a variance analysis of typical experimental data, the effects of the parameters on the maximum film stresses were obtained. Finally, the accuracy of the analysis results was verified using statistical methods. The letter j is used to represent the serial number of the corresponding factor, where j = 1, 2, 3, and 4 represent the number of component layers, the film thickness, the film hole diameter, and the staggered hole distance, respectively. The letter k is used to indicate the serial number of the corresponding factor level, and k = 1, 2, and 3 represent the values selected for factors, as shown in Table 2. An L9(34) orthogonal table, which has three levels and four factors, was selected, and nine corresponding numerical experiments were conducted. The results of the nine experiments are shown in Table 3.
The maximum film stress was the key evaluation index used to analyze the component film failure. The maximum film stress was calculated for all the numerical experiments listed in Table 3. The influence trends of the structural parameters with respect to the maximum film stress are different. A range analysis was used to analyze the influence degrees of the different structural parameters on the maximum film stress. A larger range indicated a greater impact of the corresponding structural parameter on the maximum film stress. The definition of the range, Rj, of factor j is as follows:
R j = max K j k ¯ min K j k ¯ j / k = 1 , 2 , 3 , 4 ,
where K j k ¯ represents the mean value of all the maximum film stress values when factor j is at level k, while k = 1, 2, and 3 represent the values selected for the factors in Table 2.
The effects of the different factors on the maximum film stress were obtained according to the ranges of the orthogonal experiment results. The parameter order, ranked from the largest to the smallest according to the parameter influence on the maximum film stress, is the film thickness, the film hole diameter, the number of component layers, and the staggered hole distance. To reduce the maximum stresses in the thin films, a structurally optimized design scheme for the component is proposed: D = 1.2 mm, L = 7.75 mm, N = 15, and c = 30 µm. This scheme is illustrated in Figure 17.
According to the results presented in Section 3.2, the film at the outlet has the maximum stress and is the most prone to failure. The maximum deformation of the film at the outlet of the five-unit MLI component (D = 1 mm, L = 4 mm, N = 5, and c = 0.01 mm) is 0.549 mm, and the maximum stress is 4.17 MPa. According to the modeling and calculations that were conducted based on the optimized design scheme, the maximum film deformation is 0.0017 mm, and the maximum stress is 0.098 MPa. Therefore, the optimized design produced significant decreases in the deformation and stress values.

4. Conclusions

The mechanical characteristics of the MLI component of a spacecraft during rapid decompression were investigated in this study. An FSI model for the component was established, and the internal flow field stress trends, the film deformation trends, and the film stress distribution trends were obtained for the component through FSI calculations and analyses. Then, a structurally optimized design scheme was proposed for the component, which provides an important basis for the design of spacecraft MLI components. This study produced five primary conclusions:
(1) An FSI model for the component was established. By reasonably simplifying the model and solving the modeling and mesh generation problems for the complex structure, simulation calculation and analysis methods were proposed for the component.
(2) Finite element calculation results for the flow field indicate that the film stress gradually increases in successive layers from the inlet to the outlet for any structural parameter values. Therefore, the film stress is the largest in the outlet film, which causes this film to be the part of the component that is most prone to damage and failure. The film hole diameter, the staggered hole distance, and the film layer sequence are all inversely proportional to the maximum film pressure difference, while the film thickness is positively correlated with the maximum film pressure difference. The structural parameter with the greatest influence on the maximum film pressure difference is the film hole diameter; it is followed by the number of component layers, the staggered hole distance, and finally the film thickness.
(3) The calculation results of the solid FSI model indicate that the film thickness is the factor that most critically affects the maximum film stress and that increasing the film thickness can effectively reduce the film stress. It was also verified that the nylon nets can effectively limit the film deformation, thereby better protecting the films.
(4) An orthogonal experiment was used to analyze the effects of typical component structural parameter values on the maximum film stress. The results show that the parameter with the largest influence is the film thickness; it is followed by the film hole diameter, the number of component layers, and finally the staggered hole distance. Based on this result, a structurally optimized design scheme was proposed for the component, for which the deformation and stress values decreased significantly.
(5) The method presented in this paper is generally applicable to the study of the film stresses in MLI components of spacecraft. FSI model calculations can be conducted, orthogonal experiments can be used to analyze the effects of the structural parameters on the stresses, and optimal design schemes can be proposed to avoid component failure. The influence of various structural parameters on the maximum stress of the thin film presented in this paper can be directly applied to the structural design of the MLI components of spacecraft. Under the condition that the heat insulation performance of the MLI components meets the requirements, the failure risk of the film can be effectively reduced by appropriately increasing the film thickness and the film hole diameter. However, only the effect of the fluid pressure on the component films under rapid decompression conditions during spacecraft launching was investigated in this study. The interference of external factors, such as random vibrations during spacecraft launching, was not considered. The results of this study do not reflect the overall film stress situation. In the future, the mechanical characteristics of MLI component films will be investigated under combined fluid and random vibration effects.

Author Contributions

Conceptualization, W.S., Y.L. and S.H.; methodology, W.S. and Y.L.; software, Y.L. and Q.A.; validation, W.S., S.H., F.L. and Q.A.; formal analysis, W.S. and Y.L.; investigation, W.S. and Y.L.; resources, S.H. and F.L.; data curation, Y.L. and Q.A.; writing—original draft preparation, W.S. and Y.L.; writing—review and editing, S.H., F.L. and Q.A.; visualization, W.S. and Q.A.; supervision, S.H. and F.L.; project administration, Y.L. and S.H.; funding acquisition, Y.L. and S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (No. 52075043), and the Beijing Information Science & Technology University Research Foundation (No. 2023XJJ03).

Data Availability Statement

All relevant data are within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, Y.; Ning, X.W.; Su, S. Experiment research on heat-transfer capability of the multi-layer insulation blankets under different vacuum degrees. Sci. Sin. Technol. 2015, 45, 263–267. [Google Scholar] [CrossRef] [Green Version]
  2. Yang, D.; Li, M.; Shang, K. Development of Thermal Insulation Materials Technology for Spacesuit. J. Aeronaut. Mater. 2016, 36, 87–96. [Google Scholar] [CrossRef]
  3. Chen, J. Spacesuit Engineering; National Defense Industry Press: Beijing, China, 2004; pp. 206–215. [Google Scholar]
  4. Keller, C.W.; Cunnington, G.R.; Glassford, A.P. Thermal Performance of Multilayer Insulations Gas Evacuation Characteristics of Three Selected Multilayer Insulation Composites; NASA: Sunnyvale, CA, USA, 1974; N74–N22; p. 564.
  5. Jiang, J. Multilayer Insulation Materials and Their Application to Spacecrafts. Aerosp. Mater. Technol. 2000, 2000, 17–25. [Google Scholar] [CrossRef]
  6. Abramov, I.P.; Pozdnyakov, S.S.; Severin, G.I. Main Problems of the Russian Orlan-M space suit utilization for evas on the iss. Acta Astronaut. 2001, 48, 265–273. [Google Scholar] [CrossRef] [PubMed]
  7. Gilmore, D.G. Spacecraft Thermal Control Handbook, Volume I: Fundamental Technologies; The Aerospace Press: El Segundo, CA, USA, 2002. [Google Scholar]
  8. Liu, Z.; Li, Y.; Xie, F. Thermal performance of foam/MLI for cryogenic liquid hydrogen tank during the ascent and on orbit period. Appl. Therm. Eng. 2016, 98, 430–439. [Google Scholar] [CrossRef]
  9. Johnson, W.L. Launch ascent testing of a representative Altair ascent stage methane tank. Cryogenics 2012, 52, 278–282. [Google Scholar] [CrossRef]
  10. Riesco, M.E.; Mclean, C.H.; Mills, G.L. Venting and High Vacuum Performance of Low Density Multilayer Insulation. In Proceedings of the Advances Cryogenic Engineering Conference (CEC), Tucson, AZ, USA, 9 April 2010. [Google Scholar] [CrossRef]
  11. Fesmire, J.E.; Johnson, W.L.; Meneghelli, B.J.; Coffman, B.E. Cylindrical Boiloff Calorimeters for Testing of Thermal Insulation Systems. In Proceedings of the Advances Cryogenic Engineering Conference (CEC), Tucson, AZ, USA, 28 June–2 July 2015. [Google Scholar] [CrossRef]
  12. Bapat, S.L.; Narayankhedkar, K.G.; Lukose, T.P. Performance prediction of multilayer insulation. Cryogenics 1990, 30, 700–710. [Google Scholar] [CrossRef]
  13. Wei, W.; Li, X.; Wang, R.; Li, Y. Effects of structure and shape on thermal performance of Perforated Multi-Layer Insulation Blankets. Appl. Therm. Eng. 2008, 29, 1264–1266. [Google Scholar] [CrossRef]
  14. Zhou, C.; Wang, R.; Lu, X.; Gu, A. The main influencing factors of high vacuum multilayer insulation. Vac. Cryog. 1998, 1998, 7–10. [Google Scholar]
  15. Ferrero, A.; Palestro, R.; Tavera, S. Qualification Campaign on ATV Enhanced Design MLI Blankets. In Proceedings of the 40th International Conference on Environmental Systems, Barcelona, Spain, 11–15 July 2010. [Google Scholar]
  16. Ren, H.; Zhao, X. The failure mechanism of multilayer insulation (MLI) in low-altitude environment and its protection. Spacecr. Environ. Eng. 2008, 25, 523–525. [Google Scholar] [CrossRef]
  17. Woodward, S.E.; Taylor, B.D.; Jones, T.W. A Method to have Multi-Layer Thermal Insulation Provide Damage Detection. In Proceedings of the 48th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Honolulu, HI, USA, 23–26 April 2007. [Google Scholar]
  18. Li, F.; Huang, S.; Qi, Z.; Hao, J. Fluid Simulation and Reliability Evaluation for the Spacecraft Multilayer Insulation in the Rapid Pressure Relief Process of the Powered Phase. In Proceedings of the 2022 Global Reliability and Prognostics and Health Management Conference, Yantai, China, 2 March 2022. [Google Scholar] [CrossRef]
  19. Liu, Y.; Suo, S.; Guo, F.; Huang, M.; Sun, W.; Tan, B.; Huang, S.; Li, F. Stress uniformity of spacecraft multilayer thermal insulation components. J. Tsinghua Univ. Sci. Technol. 2023, 26, 1–10. [Google Scholar] [CrossRef]
  20. Yan, S.; Hu, F.; Huang, S. Leakage and inter-stage imbalance analysis of two-stage brush seals structure based on CFD and neural network. Aeroengine 2022, 48, 70–75. [Google Scholar] [CrossRef]
  21. Huang, S.; Suo, S.; Li, Y.; Yang, J.; Liu, S.; Wang, Y. Flows in brush seals based on a 2-D staggered tube bundle model. J. Tsinghua Univ. Sci. Technol. 2016, 56, 160–166. [Google Scholar] [CrossRef]
  22. Boustani, J.; Cadieux, F.; Kenway, G.K. FSI simulations of the ASPIRE SR01 supersonic parachute flight test. Aerosp. Sci. Technol. 2022, 2022, 126. [Google Scholar] [CrossRef]
  23. Yang, H.Q.; West, J.; Harris, R. Coupled FSI Analysis of Solid Rocket Motor with Flexible Inhibitors. J. Spacecr. Rocket. 2018, 55, 303–314. [Google Scholar] [CrossRef]
  24. Diosady, L.; Murman, S.; Carton, C.; Wiart, D. A higher-order space-time finite-element method for moving-body and FSI problems. In Proceedings of the Tenth International Conference on Computational Fluid Dynamics, Barcelona, Spain, 9–13 July 2018. [Google Scholar]
  25. Sances, D.J.; Gangadharan, S.N.; Sudermann, J.E. CFD Fuel Slosh Modeling of FSI in Spacecraft Propellant Tanks with Diaphragms. In Proceedings of the 51st AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Orlando, FL, USA, 12–15 April 2010. [Google Scholar]
  26. Brehm, C.; Barad, M.F.; Kiris, C.C. An Immersed Boundary Method for Solving the Compressible Navier-Stokes Equations with FSI. In Proceedings of the 34th AIAA Applied Aerodynamics Conference, Washington, DC, USA, 13–17 June 2016. [Google Scholar] [CrossRef] [Green Version]
Figure 1. A problem with the multilayer insulation (MLI) components on the automated transfer vehicle [15].
Figure 1. A problem with the multilayer insulation (MLI) components on the automated transfer vehicle [15].
Aerospace 10 00577 g001
Figure 2. Schematic of the FSI model.
Figure 2. Schematic of the FSI model.
Aerospace 10 00577 g002
Figure 3. Schematic of the FSI model calculation method. Different colors represent the result value, with red being the largest value and blue the smallest. The specific result values are shown in the following section.
Figure 3. Schematic of the FSI model calculation method. Different colors represent the result value, with red being the largest value and blue the smallest. The specific result values are shown in the following section.
Aerospace 10 00577 g003
Figure 4. Schematic of the polygonal mesh for the global model and the nylon net edge region.
Figure 4. Schematic of the polygonal mesh for the global model and the nylon net edge region.
Aerospace 10 00577 g004
Figure 5. Mesh independence verification results.
Figure 5. Mesh independence verification results.
Aerospace 10 00577 g005
Figure 6. Schematic of the boundary conditions of the fluid model.
Figure 6. Schematic of the boundary conditions of the fluid model.
Aerospace 10 00577 g006
Figure 7. Outlet pressure–time curve.
Figure 7. Outlet pressure–time curve.
Aerospace 10 00577 g007
Figure 8. Outlet leakage–time curve. The leakage reaches a maximum at the 58th second.
Figure 8. Outlet leakage–time curve. The leakage reaches a maximum at the 58th second.
Aerospace 10 00577 g008
Figure 9. Solid model boundary conditions: (a) film boundary and (b) nylon net boundary.
Figure 9. Solid model boundary conditions: (a) film boundary and (b) nylon net boundary.
Aerospace 10 00577 g009
Figure 10. Pressure contours in the flow field: (a) N1, (b) N2, (c) N3, (d) N4, (e) N5, and (f) N6. Different shapes represent different pressure distribution.
Figure 10. Pressure contours in the flow field: (a) N1, (b) N2, (c) N3, (d) N4, (e) N5, and (f) N6. Different shapes represent different pressure distribution.
Aerospace 10 00577 g010
Figure 11. Effects of the structural parameters on the flow field pressure difference: (a) film hole diameter effect, (b) staggered hole distance effect, (c) film layer sequence effect, and (d) film thickness effect.
Figure 11. Effects of the structural parameters on the flow field pressure difference: (a) film hole diameter effect, (b) staggered hole distance effect, (c) film layer sequence effect, and (d) film thickness effect.
Aerospace 10 00577 g011
Figure 12. Deformation nephograms of the films: (a) N1, (b) N2, (c) N3, (d) N4, (e) N5, and (f) N6. Different shapes represent different deformation distribution. The maximum deformation areas are marked.
Figure 12. Deformation nephograms of the films: (a) N1, (b) N2, (c) N3, (d) N4, (e) N5, and (f) N6. Different shapes represent different deformation distribution. The maximum deformation areas are marked.
Aerospace 10 00577 g012
Figure 13. Deformation nephograms for the nylon nets (unit: mm): (a) primary deformation of M4, (b) deformation of M4 in the Z-direction, (c) deformation of M4 in the X-direction, (d) primary deformation of M5, (e) deformation of M5 in the Z-direction, and (f) deformation of M5 in the X-direction. Different shapes represent different deformation distribution. The maximum areas are marked.
Figure 13. Deformation nephograms for the nylon nets (unit: mm): (a) primary deformation of M4, (b) deformation of M4 in the Z-direction, (c) deformation of M4 in the X-direction, (d) primary deformation of M5, (e) deformation of M5 in the Z-direction, and (f) deformation of M5 in the X-direction. Different shapes represent different deformation distribution. The maximum areas are marked.
Aerospace 10 00577 g013
Figure 14. Effects of the structural parameters on the film deformations: (a) film hole diameter effect, (b) staggered hole distance effect, (c) film layer sequence effect, and (d) film thickness effect.
Figure 14. Effects of the structural parameters on the film deformations: (a) film hole diameter effect, (b) staggered hole distance effect, (c) film layer sequence effect, and (d) film thickness effect.
Aerospace 10 00577 g014
Figure 15. Stress distributions of the films (unit: MPa): (a) N1, (b) N2, (c) N3, (d) N4, (e) N5, and (f) N6. Different shapes represent different stress distribution. The maximum stress areas are marked.
Figure 15. Stress distributions of the films (unit: MPa): (a) N1, (b) N2, (c) N3, (d) N4, (e) N5, and (f) N6. Different shapes represent different stress distribution. The maximum stress areas are marked.
Aerospace 10 00577 g015
Figure 16. Effects of the structural parameters on the film stress: (a) film hole diameter effect, (b) staggered hole distance effect, (c) film layer sequence effect, and (d) film thickness effect.
Figure 16. Effects of the structural parameters on the film stress: (a) film hole diameter effect, (b) staggered hole distance effect, (c) film layer sequence effect, and (d) film thickness effect.
Aerospace 10 00577 g016
Figure 17. Schematic of the optimized design scheme for the MLI component.
Figure 17. Schematic of the optimized design scheme for the MLI component.
Aerospace 10 00577 g017
Table 1. Empirical constants and Prandtl constants.
Table 1. Empirical constants and Prandtl constants.
C1εC2εC3εCμσkσε
1.441.921.30.0911.3
Table 2. Orthogonal parameter scheme that includes four factors and three levels.
Table 2. Orthogonal parameter scheme that includes four factors and three levels.
FactorLevel
123
1Number of component layers51015
2Film thickness (10−3 mm)102030
3Film hole diameter (mm)0.811.2
4Staggered hole distance (mm)467.75
Table 3. Results of the orthogonal experiment.
Table 3. Results of the orthogonal experiment.
ModelNumber of Component LayersFilm Thickness (10−3 mm)Film Hole Diameter (mm)Staggered Hole Distance (mm)
15100.84
252016
35301.27.75
4101017.75
510201.24
610300.86
715101.26
815200.87.75
9153014
K j 1 ¯ 1.7173.3612.32.007
K j 2 ¯ 1.5280.7651.1671.34
K j 3 ¯ 1.3140.4331.0921.213
Rj0.4032.9281.2080.794
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, W.; Liu, Y.; An, Q.; Huang, S.; Li, F. Stress Characteristics and Structural Optimization of Spacecraft Multilayer Insulation Components. Aerospace 2023, 10, 577. https://doi.org/10.3390/aerospace10070577

AMA Style

Sun W, Liu Y, An Q, Huang S, Li F. Stress Characteristics and Structural Optimization of Spacecraft Multilayer Insulation Components. Aerospace. 2023; 10(7):577. https://doi.org/10.3390/aerospace10070577

Chicago/Turabian Style

Sun, Weiwei, Yue Liu, Qi An, Shouqing Huang, and Fangyong Li. 2023. "Stress Characteristics and Structural Optimization of Spacecraft Multilayer Insulation Components" Aerospace 10, no. 7: 577. https://doi.org/10.3390/aerospace10070577

APA Style

Sun, W., Liu, Y., An, Q., Huang, S., & Li, F. (2023). Stress Characteristics and Structural Optimization of Spacecraft Multilayer Insulation Components. Aerospace, 10(7), 577. https://doi.org/10.3390/aerospace10070577

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop