CN115809515A - Optimal design method for multilayer heat insulation structure of high-speed aircraft - Google Patents

Optimal design method for multilayer heat insulation structure of high-speed aircraft Download PDF

Info

Publication number
CN115809515A
CN115809515A CN202310085526.3A CN202310085526A CN115809515A CN 115809515 A CN115809515 A CN 115809515A CN 202310085526 A CN202310085526 A CN 202310085526A CN 115809515 A CN115809515 A CN 115809515A
Authority
CN
China
Prior art keywords
heat insulation
layer
insulation structure
optimization
temperature
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310085526.3A
Other languages
Chinese (zh)
Other versions
CN115809515B (en
Inventor
郝东
王琪
苏杰
刘建霞
丁智坚
武龙
贺元元
王�锋
余婧
杨大伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Aerospace Technology of China Aerodynamics Research and Development Center
Original Assignee
Institute of Aerospace Technology of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Aerospace Technology of China Aerodynamics Research and Development Center filed Critical Institute of Aerospace Technology of China Aerodynamics Research and Development Center
Priority to CN202310085526.3A priority Critical patent/CN115809515B/en
Publication of CN115809515A publication Critical patent/CN115809515A/en
Application granted granted Critical
Publication of CN115809515B publication Critical patent/CN115809515B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Laminated Bodies (AREA)

Abstract

The invention discloses an optimal design method of a high-speed aircraft multi-layer heat insulation structure, which relates to the technical field of heat insulation, wherein the high-speed aircraft heat insulation structure is divided into multiple layers of materials, and design variables, namely a selectable material database and the geometric dimensions of each layer of material, are determined; determining an optimized objective function heat insulation efficiency-mass ratio and setting constraint conditions; the optimization objective function and the constraint condition are gathered, and an optimization model is obtained according to the design variables; and solving the optimization model by adopting an intelligent optimization algorithm to obtain optimized heat insulation material selection and the geometric dimension of each layer of heat insulation material. The method introduces the material selection, the geometric dimension and the temperature change characteristic of the thermal conductivity of the material into the optimization design model to obtain the optimized material selection and the geometric dimension of the material, thereby better meeting the requirements of practical application.

Description

Optimal design method for multilayer heat insulation structure of high-speed aircraft
Technical Field
The invention relates to the technical field of heat insulation, in particular to a method for optimally designing a multilayer heat insulation structure of a high-speed aircraft.
Background
The heat insulation structure is a key factor for protecting the high-speed aircraft to safely work and complete various tasks, and the effective heat insulation structure can meet the basic requirements of heat prevention and heat insulation of the high-speed aircraft, so that the high-speed aircraft is prevented from being damaged by an extreme aerodynamic thermal environment. With the development and upgrade of the national high-speed aircraft technology, the traditional heat insulation structure is difficult to meet the requirements of low cost, light weight, integration of prevention and heat insulation and the like, so that the design of the multilayer heat insulation structure is widely concerned by various countries. The schematic diagram of the multilayer heat insulation structure of the high-speed aircraft is shown in fig. 1, the high-temperature resistant layer of the multilayer heat insulation structure is heated in a forced convection heat exchange manner in an extreme high-temperature thermal environment borne by the high-speed aircraft in flight, and meanwhile, the burning resistant layer performs radiation heat dissipation to the outer space to remove part of heat. The main heat in the remaining heat is blocked by the insulating layer, and a small part of heat is absorbed and stored by the material of the burning-resistant layer and is delivered to the base layer, and the base layer structure and the internal air perform natural convection heat exchange. The quality of the heat insulation performance of the multilayer heat insulation structure mainly depends on the material parameters of the heat insulation layer and the thickness of each layer, and the equivalent mechanical property of the whole structure is mainly determined by the material parameters of the base layer, so that the base layer is sometimes called as a bearing layer according to the function of the base layer.
At present, the optimized design of the high-speed aircraft heat insulation structure mostly considers that the thermal conductivity of materials is constant, the thermal conductivity of the materials is sensitive to temperature, the heat transfer performance changes violently along with the temperature, the thermal conductivity of the existing aircraft heat insulation material does not consider the condition along with the temperature change, the temperature change characteristic of the thermal conductivity of the heat insulation material influences the heat insulation effect, and the designed high-speed aircraft multilayer heat insulation structure does not consider the influence of the temperature on the thermal conductivity and has larger difference with the actual heat insulation performance; in addition, when the design optimization of the heat insulation structure is carried out, the current material selection is determined, the material preference is not considered, and the material selection and the size variable are not simultaneously taken into the design optimization of the heat insulation structure.
The existing optimization design method cannot simultaneously solve the problems of material optimization and structure optimization, and has large limitation.
Disclosure of Invention
The invention aims to: aiming at the optimal design of the multilayer heat insulation structure of the high-speed aircraft, a material selection design variable and a material heat conductivity temperature variation characteristic are introduced into an optimal design model, so that the problem that the optimization design method of the existing multilayer heat insulation structure of the high-speed aircraft does not consider the material optimization and the influence of the material heat conductivity temperature variation is solved.
The technical scheme adopted by the invention is as follows:
the invention relates to an optimal design method of a high-speed aircraft multi-layer heat insulation structure, the heat insulation structure of the high-speed aircraft is divided into a plurality of layers of materials,
determining design variables of the heat insulation structure, namely a selectable material database and the geometric dimensions of materials of each layer;
determining an optimized objective function heat insulation mass ratio, and setting constraint conditions;
the optimization objective function and the constraint condition are gathered, and an optimization model is obtained according to the design variables;
and solving the optimization model by adopting an intelligent optimization algorithm to obtain optimized design variables, namely optimized material selection and the geometric dimension of each layer of heat insulation material.
Further, the multilayer material specifically comprises an outer layer material, an intermediate layer material and an inner layer material, and the selectable material data are respectively
Figure SMS_1
Figure SMS_2
And
Figure SMS_3
then the material database can be selected as
Figure SMS_4
The geometric dimensions of the three layers of materials are respectively the thickness of the outer layer of material
Figure SMS_5
Thickness of the interlayer Material
Figure SMS_6
Thickness of inner layer material
Figure SMS_7
Further, the optimization objective function thermal insulation efficiency-to-mass ratio, that is, the mass ratio between the thermal insulation efficiency and the thermal insulation structure, is given by the formula:
Figure SMS_8
the heat insulation efficiency is the effective working time length
Figure SMS_9
The thermal insulation rate is as follows:
Figure SMS_10
wherein ,
Figure SMS_11
is the temperature of the heat source and is,
Figure SMS_12
inner wall temperature response;
the heat insulation structure comprises the following components in parts by mass:
Figure SMS_13
wherein ,
Figure SMS_14
the density of the material of the outer layer,
Figure SMS_15
the density of the material of the intermediate layer,
Figure SMS_16
is the density of the material of the inner layer,
according to the concept of the maximum heat insulation efficiency-to-mass ratio, the obtained optimization objective function is as follows:
Figure SMS_17
further, the constraint conditions are specifically:
the inner layer of the heat insulation structure has the temperature of the inner surface wall at one side of the inner cavity
Figure SMS_18
Time of day temperature
Figure SMS_19
Is less than or equal to
Figure SMS_20
The thickness dimension of each layer of material satisfies
Figure SMS_22
Figure SMS_25
Figure SMS_27
, wherein ,
Figure SMS_23
Figure SMS_26
Figure SMS_28
the upper limit of the thickness dimension of each layer of material,
Figure SMS_29
Figure SMS_21
Figure SMS_24
the lower limit of the thickness dimension of each layer of material;
equivalent modulus
Figure SMS_30
The equivalent modulus is defined as:
Figure SMS_31
wherein ,
Figure SMS_32
is the minimum value of the modulus body allowed,
Figure SMS_33
is the modulus of the material of the outer layer,
Figure SMS_34
is the modulus of the material of the intermediate layer,
Figure SMS_35
is the modulus of the inner layer material;
material selection belongs to the Material data set
Figure SMS_36
Furthermore, the temperature of the inner layer of the heat insulation structure is lower than that of the inner surface wall of one side of the internal cavity
Figure SMS_37
Time of day temperature
Figure SMS_38
It is necessary to solve through a heat transfer equation,
the heat transfer equation:
Figure SMS_39
the boundary conditions are as follows:
Figure SMS_40
and
Figure SMS_41
wherein ,
Figure SMS_43
in order to obtain the thermal conductivity of the material,
Figure SMS_47
the specific heat capacity is the specific heat capacity,
Figure SMS_51
it is the temperature that is set for the purpose,
Figure SMS_44
as a matter of time, the time is,
Figure SMS_48
is a coordinate of the space, and is,
Figure SMS_52
the thermal conductivity of the material of the inner layer,
Figure SMS_54
the heat exchange coefficient between the inner layer structure and the air,
Figure SMS_42
is the thermal conductivity of the material of the outer layer,
Figure SMS_46
the heat exchange coefficient between the outer layer structure and the air,
Figure SMS_50
is the stefin-boltzmann constant,
Figure SMS_53
it is the temperature of the air inside that,
Figure SMS_45
the temperature of the outer wall surface of the outer layer,
Figure SMS_49
either the outside air temperature or the heat source temperature.
Further, the optimization objective function and the constraint condition are aggregated to obtain an optimization design mathematical model:
Figure SMS_55
the design variables are defined as:
Figure SMS_56
the optimization model is defined as:
Figure SMS_57
furthermore, the design variable is a mixed type of coexistence of discrete variable and continuous variable, which is a mixed design variable optimization problem, and an intelligent optimization algorithm is adopted to solve the optimization model to obtain optimized design variables, namely, the material selection variable M and the geometric dimension of the three layers of materials
Figure SMS_58
Figure SMS_59
And
Figure SMS_60
the intelligent optimization algorithm is a genetic algorithm or a particle swarm algorithm.
In summary, due to the adoption of the technical scheme, the invention has the beneficial effects that:
the invention relates to an optimal design method of a high-speed aircraft multi-layer heat insulation structure, aiming at the optimal design of the multi-layer heat insulation structure of the high-speed aircraft, the material selection, the geometric dimension and the temperature variation characteristic of the material heat conductivity are introduced into an optimal design model, the problem that the material optimization and the temperature variation influence of the material heat conductivity are not considered in the existing optimal design method of the high-speed aircraft multi-layer heat insulation structure is solved, the optimized material selection and the optimized material geometric dimension are obtained, and the requirements of practical application are better met.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and that for those skilled in the art, other relevant drawings can be obtained according to the drawings without inventive effort, wherein:
FIG. 1 is a schematic view of a high speed aircraft multi-layer insulation structure;
FIG. 2 is a flow chart of the optimization method of the present invention;
FIG. 3 is a diagram of the solution optimized using the Monte Carlo method;
FIG. 4 is a result of a solution using PSO optimization.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the detailed description and specific examples, while indicating the preferred embodiment of the invention, are intended for purposes of illustration only and are not intended to limit the scope of the invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
It is to be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
As shown in fig. 2, the invention is an optimized design method of a multilayer heat insulation structure of a high-speed aircraft, the heat insulation structure of the high-speed aircraft is divided into multiple layers of materials,
determining design variables of the heat insulation structure, namely a selectable material database and the geometric dimensions of materials of each layer;
specifically, as shown in fig. 1, the multilayer material specifically includes an outer layer material, an intermediate layer material and an inner layer material, and the selectable material data are respectively
Figure SMS_61
Figure SMS_62
And
Figure SMS_63
then the material database can be selected as
Figure SMS_64
The geometric dimensions of the three layers are respectively the thickness of the outer layer material
Figure SMS_65
Thickness of the interlayer Material
Figure SMS_66
Thickness of inner layer material
Figure SMS_67
Determining an optimized objective function heat insulation mass ratio, and setting constraint conditions;
specifically, the optimization objective function is a heat insulation efficiency/heat insulation structure mass ratio, and the formula is as follows:
Figure SMS_68
the heat insulation efficiency is the effective working time length
Figure SMS_69
The thermal insulation rate was:
Figure SMS_70
wherein ,
Figure SMS_71
is the temperature of the heat source,
Figure SMS_72
inner wall temperature response;
the quality of the heat insulation structure is as follows:
Figure SMS_73
wherein ,
Figure SMS_74
the density of the material of the outer layer,
Figure SMS_75
the density of the material of the intermediate layer,
Figure SMS_76
is the density of the material of the inner layer,
according to the concept of the maximum heat insulation efficiency-to-mass ratio, the obtained optimization objective function is as follows:
Figure SMS_77
the optimization objective function is the ratio of quality to efficiency, so the smaller the optimization result of the optimization objective function is, the better the optimization result is.
Specifically, the constraint condition is specifically:
heat insulation structureThe temperature of the inner layer is within the range of the inner surface wall of one side of the internal cavity
Figure SMS_78
Time of day temperature
Figure SMS_79
Is less than or equal to
Figure SMS_80
The thickness dimension of each layer of material satisfies
Figure SMS_82
Figure SMS_86
Figure SMS_88
, wherein ,
Figure SMS_83
Figure SMS_85
Figure SMS_87
the upper limit of the thickness dimension of each layer of material,
Figure SMS_89
Figure SMS_81
Figure SMS_84
the lower limit of the thickness dimension of each layer of material;
equivalent modulus
Figure SMS_90
The equivalent modulus is defined as:
Figure SMS_91
wherein ,
Figure SMS_92
is the minimum value of the modulus body allowed,
Figure SMS_93
is the modulus of the material of the outer layer,
Figure SMS_94
is the modulus of the material of the intermediate layer,
Figure SMS_95
is the modulus of the inner layer material;
material selection belongs to the material data set
Figure SMS_96
Specifically, the inner layer of the heat insulation structure has the inner surface wall temperature towards one side of the internal cavity body
Figure SMS_97
Time of day temperature
Figure SMS_98
It needs to be solved by the heat transfer equation,
the heat transfer equation:
Figure SMS_99
the boundary conditions are as follows:
Figure SMS_100
and
Figure SMS_101
wherein ,
Figure SMS_104
in order to obtain the thermal conductivity of the material,
Figure SMS_107
the specific heat capacity is the specific heat capacity,
Figure SMS_111
is the temperature of the liquid to be treated,
Figure SMS_103
as a matter of time, the time is,
Figure SMS_109
is a coordinate of the space, and is,
Figure SMS_112
the thermal conductivity of the material of the inner layer,
Figure SMS_114
the heat exchange coefficient between the inner layer structure and the air,
Figure SMS_102
is the thermal conductivity of the material of the outer layer,
Figure SMS_106
the heat exchange coefficient between the outer layer structure and the air,
Figure SMS_110
is the stefin-boltzmann constant,
Figure SMS_113
the inside air temperature (background temperature in figure 2),
Figure SMS_105
the temperature of the outer wall surface of the outer layer,
Figure SMS_108
either the outside air temperature or the heat source temperature.
The partial differential equation (heat transfer equation) is solved by adopting a finite difference implicit format, and the spatial region sequence is represented by i, and the time sequence is represented by j. The space grid line and the time grid line (i, j) represent a node of a time-space, and the state of the j moment at the space position i replaces the derivative by finite difference, so that the partial differential equation is converted into a differential equation:
Figure SMS_115
using taylor expansion, the finite difference expression of its derivative is:
Figure SMS_116
Figure SMS_117
the above two formulas are known:
Figure SMS_118
the difference in partial differentiation with respect to time is:
Figure SMS_119
the above results are substituted into a heat transfer differential equation, and an iterative recursion format in a discrete form can be obtained:
Figure SMS_120
recur in time by 1, i.e. j becomes j +1, have
Figure SMS_121
Then there are:
Figure SMS_122
the above formula is a finite difference iterative recurrence formula.
The optimization objective function and the constraint condition are gathered, and an optimization model is obtained according to the design variables;
specifically, the optimization objective function and the constraint condition are aggregated to obtain an optimization design mathematical model:
Figure SMS_123
the design variables are defined as:
Figure SMS_124
the optimization model is defined as:
Figure SMS_125
and solving the optimization model by adopting an intelligent optimization algorithm to obtain optimized design variables, namely optimized material selection and the geometric dimension of each layer of heat insulation material.
Specifically, the design variable is a mixed type of coexistence of a discrete variable and a continuous variable, and is a mixed design variable optimization problem, and an intelligent optimization algorithm is adopted to solve an optimization model to obtain an optimized design variable, namely a material selection variable M and the geometric dimension of three layers of materials
Figure SMS_126
Figure SMS_127
And
Figure SMS_128
the intelligent optimization algorithm is a genetic algorithm or a particle swarm algorithm, but is not limited to such an optimization algorithm.
Example of the implementation
The optimization method is verified and analyzed by using an example. The parameter values in the optimization model are:
Figure SMS_129
=0.5mm、
Figure SMS_134
=3mm、
Figure SMS_138
=5mm、
Figure SMS_130
=50mm、
Figure SMS_133
=2mm、
Figure SMS_137
=5mm,
Figure SMS_140
=10Gpa,
Figure SMS_131
=1200℃,
Figure SMS_135
=20℃,
Figure SMS_139
=100 ℃, heat exchange coefficient of outer layer structure and air
Figure SMS_141
= 230W/(m 2. DEG C), and the heat exchange coefficient of the inner layer structure and the air
Figure SMS_132
=10 W/(m2·℃),
Figure SMS_136
=300s。
Outer layer materials data set table 1:
TABLE 1 ablation resistance layer (outer layer) Material parameters (outer layer)
Name of Material Density (kg)·m-3) Specific heat capacity (J. (kg. K) -1) Thermal conductivity coefficient (W. (m. K) -1) Young's modulus (GPa)
C/C composite material 1600 713 1.6 22.32
C/SiC composite material 2100 1420 5 49.8
ZrB 2-based composite material 6300 430 53 490
Inner layer material data set:
TABLE 2 bearing layer (inner layer) Material parameters
Name of Material Density (kg. M-3)) Specific heat capacity (J. (kg. K) -1) Thermal conductivity coefficient (W. (m. K) -1) Young's modulus (GPa)
Aluminum alloy 7075 2800 962 129.4 71
GH4099 8470 624 18.9 175
Titanium alloy TC4 4440 659 9.13 115.85
Interlayer material dataset:
TABLE 3 thermal barrier layer (interlayer) Material parameters
Name of Material Density (kg. M-3) Specific heat capacity (J. (kg. K) -1) Thermal conductivity coefficient (W. (m. K) -1) Young's modulus (GPa)
Carbon aerogels (carbon aerogels) 811 870 3.521505E-7T2+7.847416E-4T+1.104313 2.84
SiO2 aerogel 140 1050 1.135885E-10T3-2.256084E-8T2+1.640566E-5T+1.481429E-2 0.010
SiC aerogel 9.7 579 5.269242E-11T3-2.221319E-8T2+1.622211E-5T+2.377993E-2 0.013
And (4) optimizing and solving by adopting a Monte Carlo method to obtain 200 feasible solution optimization results, wherein the optimization results are shown in figure 3.
With PSO optimization, the population number is 50, the number of iterations is 200, and the optimization result is shown in fig. 4.
The PSO optimization method has better effect than the Monte Carlo method, and the PSO optimization method is adopted in the invention.
According to the optimization result of the PSO, the optimal value of the optimization objective function is 10.9167, and the corresponding design variables are as follows: x = (2,2,1,5.4070e-04, 0.0120, 0.0020), namely that the outer layer material is selected from C/SiC composite material, the intermediate layer material is selected from carbon aerogel, the inner layer material is selected from GH4099, the thickness of the outer layer is 5.4070e-01mm, the thickness of the intermediate layer is 12mm, and the thickness of the inner layer is 2mm.
The above description is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be made by those skilled in the art without inventive work within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope defined by the claims.

Claims (7)

1. A method for optimally designing a multilayer heat insulation structure of a high-speed aircraft is characterized in that the heat insulation structure of the high-speed aircraft is divided into multiple layers of materials,
determining design variables of the heat insulation structure, namely a selectable material database and the geometric dimensions of materials of each layer;
determining an optimized objective function heat insulation mass ratio, and setting constraint conditions;
the optimization objective function and the constraint condition are aggregated, and an optimization model is obtained according to design variables;
and solving the optimization model by adopting an intelligent optimization algorithm to obtain optimized design variables, namely optimized material selection and the geometric dimension of each layer of heat insulation material.
2. The method for optimally designing the multilayer heat insulation structure of the high-speed aircraft according to claim 1, wherein the method comprises the following steps of: the multilayer material specifically comprises an outer layer material, an intermediate layer material and an inner layer material, and the selectable material data are respectively
Figure QLYQS_1
Figure QLYQS_2
And
Figure QLYQS_3
then, thenThe database of selectable materials is
Figure QLYQS_4
The geometric dimensions of the three layers are respectively the thickness of the outer layer material
Figure QLYQS_5
Thickness of the interlayer Material
Figure QLYQS_6
Thickness of inner layer material
Figure QLYQS_7
3. The method for optimally designing the multilayer heat insulation structure of the high-speed aircraft according to claim 2, wherein the optimal objective function is the heat insulation efficiency-to-mass ratio of the heat insulation structure, and the formula is as follows:
Figure QLYQS_8
the heat insulation efficiency is the effective working time length
Figure QLYQS_9
The thermal insulation rate was:
Figure QLYQS_10
wherein ,
Figure QLYQS_11
is the temperature of the heat source,
Figure QLYQS_12
inner wall temperature response;
the heat insulation structure comprises the following components in parts by mass:
Figure QLYQS_13
wherein ,
Figure QLYQS_14
the density of the material of the outer layer,
Figure QLYQS_15
the density of the material of the intermediate layer,
Figure QLYQS_16
is the density of the material of the inner layer,
according to the concept of the maximum heat insulation efficiency-to-mass ratio, the obtained optimization objective function is as follows:
Figure QLYQS_17
4. the method for optimally designing the multilayer heat insulation structure of the high-speed aircraft according to claim 3, wherein the constraint conditions are specifically as follows:
the inner layer of the heat insulation structure has the temperature of the inner surface wall at one side of the inner cavity
Figure QLYQS_18
Time of day temperature
Figure QLYQS_19
Is less than or equal to
Figure QLYQS_20
The thickness dimension of each layer of material satisfies
Figure QLYQS_23
Figure QLYQS_26
Figure QLYQS_28
, wherein ,
Figure QLYQS_22
Figure QLYQS_25
Figure QLYQS_27
the upper limit of the thickness dimension of each layer of material,
Figure QLYQS_29
Figure QLYQS_21
Figure QLYQS_24
the lower limit of the thickness dimension of each layer of material;
equivalent modulus
Figure QLYQS_30
The equivalent modulus is defined as:
Figure QLYQS_31
wherein ,
Figure QLYQS_32
is the minimum value of the modulus body allowed,
Figure QLYQS_33
is the modulus of the material of the outer layer,
Figure QLYQS_34
is the modulus of the material of the intermediate layer,
Figure QLYQS_35
is the modulus of the inner layer material;
material selection falls within material dataCollection
Figure QLYQS_36
5. The method for optimally designing the multilayer heat insulation structure of the high-speed aircraft according to claim 4, wherein the temperature of the inner surface wall of the inner layer of the heat insulation structure on one side of the inner cavity is within the range of the temperature of the inner surface wall on one side of the inner cavity
Figure QLYQS_37
Time of day temperature
Figure QLYQS_38
It is necessary to solve through a heat transfer equation,
the heat transfer equation:
Figure QLYQS_39
the boundary conditions are as follows:
Figure QLYQS_40
and
Figure QLYQS_41
wherein ,
Figure QLYQS_44
in order to obtain the thermal conductivity of the material,
Figure QLYQS_49
the specific heat capacity is the specific heat capacity,
Figure QLYQS_52
is the temperature of the liquid to be treated,
Figure QLYQS_45
in the form of a time, the time,
Figure QLYQS_48
is a coordinate of the space, and is,
Figure QLYQS_51
the thermal conductivity of the material of the inner layer,
Figure QLYQS_54
the heat exchange coefficient between the inner layer structure and the air,
Figure QLYQS_42
is the thermal conductivity of the material of the outer layer,
Figure QLYQS_46
the heat exchange coefficient between the outer layer structure and the air,
Figure QLYQS_50
is the stefin-boltzmann constant,
Figure QLYQS_53
it is the temperature of the air inside that,
Figure QLYQS_43
the temperature of the outer wall surface of the outer layer,
Figure QLYQS_47
either the outside air temperature or the heat source temperature.
6. The method for optimally designing the multilayer heat insulation structure of the high-speed aircraft according to claim 4, wherein the optimization objective function and the constraint conditions are combined to obtain an optimal design mathematical model:
Figure QLYQS_55
the design variables are defined as:
Figure QLYQS_56
the optimization model is defined as:
Figure QLYQS_57
7. the method for optimally designing the multilayer heat insulation structure of the high-speed aircraft according to claim 6, wherein the method comprises the following steps of: the design variable is a mixed type of coexistence of discrete variable and continuous variable, is a mixed design variable optimization problem, and adopts an intelligent optimization algorithm to carry out optimization model solution to obtain optimized design variables, namely a material selection variable M and the geometric dimension of three layers of materials
Figure QLYQS_58
Figure QLYQS_59
And
Figure QLYQS_60
the intelligent optimization algorithm is a genetic algorithm or a particle swarm algorithm.
CN202310085526.3A 2023-02-09 2023-02-09 Multi-layer heat insulation structure optimization design method for high-speed aircraft Active CN115809515B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310085526.3A CN115809515B (en) 2023-02-09 2023-02-09 Multi-layer heat insulation structure optimization design method for high-speed aircraft

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310085526.3A CN115809515B (en) 2023-02-09 2023-02-09 Multi-layer heat insulation structure optimization design method for high-speed aircraft

Publications (2)

Publication Number Publication Date
CN115809515A true CN115809515A (en) 2023-03-17
CN115809515B CN115809515B (en) 2023-05-02

Family

ID=85487728

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310085526.3A Active CN115809515B (en) 2023-02-09 2023-02-09 Multi-layer heat insulation structure optimization design method for high-speed aircraft

Country Status (1)

Country Link
CN (1) CN115809515B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116052820A (en) * 2023-03-22 2023-05-02 中国空气动力研究与发展中心计算空气动力研究所 Material thermal performance evaluation method, device, equipment and readable storage medium
CN117634313A (en) * 2023-12-15 2024-03-01 浙江阿斯克建材科技股份有限公司 Heat insulation material modularized design method and system

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105843996A (en) * 2016-03-18 2016-08-10 中国飞机强度研究所 Optimization method of heat-proof structure
CN106021734A (en) * 2016-05-23 2016-10-12 北京航空航天大学 Multi-objective reliability optimum design method for forehold thermal protection system of reusable launch vehicle
WO2016190176A1 (en) * 2015-05-22 2016-12-01 旭硝子株式会社 Layered heat insulator having through hole, and heat insulating structure
WO2017072688A1 (en) * 2015-10-28 2017-05-04 Termokir Industries (1980) Ltd. Wall profile generator for structures
CN106874588A (en) * 2017-02-08 2017-06-20 北京航空航天大学 A kind of analysis of multilayer thermal protection system non-probabilistic uncertainty and Optimization Design based on experimental design
CN106971015A (en) * 2016-01-13 2017-07-21 南京航空航天大学 A kind of re-entry space vehicle thermal protection system is quickly analyzed and design method
CN212685893U (en) * 2020-07-02 2021-03-12 中南大学 Hypersonic aircraft prevents thermal-insulated integral structure
CN113502144A (en) * 2021-08-17 2021-10-15 北京理工大学 Directional heat conduction and insulation material and preparation method thereof
CN114741915A (en) * 2022-03-17 2022-07-12 上海机电工程研究所 Optimization method and system for thermal protection scheme of heat-proof and heat-insulation integrated structure

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016190176A1 (en) * 2015-05-22 2016-12-01 旭硝子株式会社 Layered heat insulator having through hole, and heat insulating structure
WO2017072688A1 (en) * 2015-10-28 2017-05-04 Termokir Industries (1980) Ltd. Wall profile generator for structures
CN106971015A (en) * 2016-01-13 2017-07-21 南京航空航天大学 A kind of re-entry space vehicle thermal protection system is quickly analyzed and design method
CN105843996A (en) * 2016-03-18 2016-08-10 中国飞机强度研究所 Optimization method of heat-proof structure
CN106021734A (en) * 2016-05-23 2016-10-12 北京航空航天大学 Multi-objective reliability optimum design method for forehold thermal protection system of reusable launch vehicle
CN106874588A (en) * 2017-02-08 2017-06-20 北京航空航天大学 A kind of analysis of multilayer thermal protection system non-probabilistic uncertainty and Optimization Design based on experimental design
CN212685893U (en) * 2020-07-02 2021-03-12 中南大学 Hypersonic aircraft prevents thermal-insulated integral structure
CN113502144A (en) * 2021-08-17 2021-10-15 北京理工大学 Directional heat conduction and insulation material and preparation method thereof
CN114741915A (en) * 2022-03-17 2022-07-12 上海机电工程研究所 Optimization method and system for thermal protection scheme of heat-proof and heat-insulation integrated structure

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
李世斌等: "高速飞行器组合式热防护系统研究进展" *
杨明等: "航天器多层绝热材料的非稳态传热分析" *
柴峻: "临近空间飞行器多层隔热结构传热分析及优化设计" *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116052820A (en) * 2023-03-22 2023-05-02 中国空气动力研究与发展中心计算空气动力研究所 Material thermal performance evaluation method, device, equipment and readable storage medium
CN117634313A (en) * 2023-12-15 2024-03-01 浙江阿斯克建材科技股份有限公司 Heat insulation material modularized design method and system

Also Published As

Publication number Publication date
CN115809515B (en) 2023-05-02

Similar Documents

Publication Publication Date Title
CN115809515A (en) Optimal design method for multilayer heat insulation structure of high-speed aircraft
CN109960878B (en) Hypersonic aircraft-oriented overall active/passive thermal protection system coupling design method
Qian et al. Topology optimization of a coupled thermal-fluid system under a tangential thermal gradient constraint
CN111967089B (en) Irregular three-dimensional lattice structure heat insulation performance optimization design method based on gradual change cell size
CN114547790B (en) Calculation method for evaluating heat insulation performance of complex multilayer heat protection structure
CN106570252B (en) A kind of thermal protection system design method based on probabilistic technique
CN114781219A (en) Phase change heat sink design method based on topology optimization
Li et al. Characteristics and optimation of heat pipe radiator for space nuclear propulsion spacecraft
Mahdi et al. Modeling the thermal behavior of solar salt in electrical resistance heaters for the application in PV-CSP hybrid power plants
CN117113567A (en) Design method of self-adaptive plate heat exchanger
CN112836399B (en) Nonlinear contact thermal resistance analysis solving method based on finite element algorithm
CN109992812A (en) A kind of Heat Transfer Calculation of cylinder thermal insulation layer
CN112728297B (en) Device for improving heat insulation effect based on additive manufacturing
CN110059372A (en) A kind of objective design method of the shell-and-tube heat exchanger based on differential evolution algorithm
Bulut et al. Multi-objective optimization of liquid cooling system for a twelve-cell battery module
CN112131801A (en) Thermal design method for box type energy storage system
JP2004311885A (en) Heat sink and shape calculating method for the same
Mekki et al. Voxel-Based Topology Optimization of Heat Exchanger Fins
CN115577652A (en) Gas shell type cooling structure of turbine box and design method thereof
CN114757125B (en) Self-organizing heat sink structure design method based on diffusion back diffusion system
CN110600074B (en) Heat spreading device for heat conduction and heat radiation
Chen et al. Effect of low wick permeability on transient and steady-state performance of heat pipes
Yun et al. New and advanced materials and technologies in ultralow-energy buildings
Giancarli et al. Conceptual design of a high temperature water-cooled divertor for a fusion power reactor
CN117709169B (en) Determination method of composite furnace lining temperature rise curve and application thereof

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant