CN115809515A - 一种高速飞行器多层隔热结构优化设计方法 - Google Patents
一种高速飞行器多层隔热结构优化设计方法 Download PDFInfo
- 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
Links
- 238000009413 insulation Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000000463 material Substances 0.000 claims abstract description 126
- 238000005457 optimization Methods 0.000 claims abstract description 67
- 239000012774 insulation material Substances 0.000 claims abstract description 7
- 239000010410 layer Substances 0.000 claims description 109
- 239000011229 interlayer Substances 0.000 claims description 5
- 230000002068 genetic effect Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 2
- 238000013480 data collection Methods 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 3
- 239000002131 composite material Substances 0.000 description 4
- 239000004966 Carbon aerogel Substances 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 239000004964 aerogel Substances 0.000 description 2
- 239000011204 carbon fibre-reinforced silicon carbide Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052681 coesite Inorganic materials 0.000 description 1
- 229910052906 cristobalite Inorganic materials 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 235000012239 silicon dioxide Nutrition 0.000 description 1
- 229910052682 stishovite Inorganic materials 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 229910052905 tridymite Inorganic materials 0.000 description 1
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Laminated Bodies (AREA)
Abstract
本发明公开了一种高速飞行器多层隔热结构优化设计方法,涉及隔热技术领域,本发明的高速飞行器隔热结构分为多层材料,确定设计变量,即可选材料数据库和各层材料的几何尺寸;确定优化目标函数隔热效质比,并设置约束条件;将优化目标函数、约束条件集合起来,根据设计变量,得到优化模型;采用智能优化算法求解优化模型,得到优化的隔热材料选择和各层隔热材料的几何尺寸。本发明将材料选择、几何尺寸和材料热导率温变特性引入优化设计模型之中,得到优化的材料选择和材料几何尺寸,更符合实际应用的需求。
Description
技术领域
本发明涉及隔热技术领域,特别是一种高速飞行器多层隔热结构优化设计方法。
背景技术
隔热结构是保护高速飞行器安全工作和完成各项任务的关键因素,有效的隔热结构可以满足高速飞行器防热和隔热的基本需求,使其免于极端气动热环境的伤害。随着国家高速飞行器技术的发展和升级,传统的隔热结构难以满足低成本、轻量化、防/隔热一体化等方面的要求,因此多层隔热结构的设计受到各国广泛关注。高速飞行器多层隔热结构示意图如图1所示,高速飞行器在飞行中所承受的极端高温热环境通过强制对流换热方式对多层隔热结构的耐高温层进行加热,同时耐烧灼层又向外部空间进行辐射散热,排走部分热量。余下热量中的主要热量被隔热层阻挡,小部分热量由耐烧灼层自身材料吸收存储和递给基层,基层结构与内部空气进行自然对流换热。多层隔热结构隔热性能的优劣主要取决于隔热层材料参数和各层厚度尺寸,整体结构的等效力学性能主要由基层材料参数决定,因此基层有时候根据其功能又称为承力层。
在高速飞行器多层隔热结构中,每一层的材料往往具有多种材料可供选择,目前,对于高速飞行器隔热结构的优化设计大多考虑的是材料热导率是常数,材料热导率对温度敏感,传热性能随温度变化剧烈,现有的飞行器隔热材料热导率没有考虑随温度变化的情况,而隔热材料热导率的温变特性对隔热效果是有影响的,没有考虑温度对热导率影响而设计得到的高速飞行器多层隔热结构往往与实际的隔热性能差距较大;此外,在进行隔热结构设计优化时,目前的材料选型是确定的,没有考虑材料优选的情况,没有将材料选择和尺寸变量同时纳入隔热结构优化设计变量。
现在的优化设计方法不能同时解决材料优选和结构优化的问题,局限性较大。
发明内容
本发明的目的在于:提出一种高速飞行器多层隔热结构优化设计方法,针对高速飞行器的多层隔热类结构优化设计,将材料选择设计变量和材料热导率温变特性引入优化设计模型之中,解决现有高速飞行器多层隔热结构优化设计方法没有考虑材料优选和材料热导率温变影响的问题。
本发明采用的技术方案如下:
本发明是一种高速飞行器多层隔热结构优化设计方法,所述高速飞行器的隔热结构分为多层材料,
确定隔热结构的设计变量,即可选材料数据库和各层材料的几何尺寸;
确定优化目标函数隔热效质比,并设置约束条件;
将优化目标函数、约束条件集合起来,根据设计变量,得到优化模型;
采用智能优化算法求解优化模型,得到优化的设计变量,即优化的材料选择和各层隔热材料的几何尺寸。
进一步的,所述优化目标函数隔热效质比,即隔热效率和隔热结构的质量比,公式为:
所述隔热结构的质量为:
根据最大隔热效质比的概念,得到优化目标函数为:
进一步的,所述约束条件具体为:
所述传热方程:
边界条件为:
其中,为材料热导率,为比热容,为温度,为时间,为空间坐标,为内层材料热导率,为内层结构与空气换热系数,为外层材料热导率,为外层结构与空气的换热系数,为史蒂芬-玻尔兹曼常数,为内部空气温度,为外层外壁面温度,为外部空气温度或热源温度。
进一步的,将所述优化目标函数、约束条件集合起来,得到优化设计数学模型:
设计变量定义为:
则优化模型定义为:
进一步的,所述设计变量为离散变量和连续变量共存的混合型,是一种混合设计变量式优化问题,采用智能优化算法进行优化模型求解,得到优化的设计变量,即材料选择变量M和三层材料的几何尺寸、和,所述智能优化算法为遗传算法或粒子群算法。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
本发明是一种高速飞行器多层隔热结构优化设计方法,针对高速飞行器的多层隔热类结构优化设计,将材料选择、几何尺寸和材料热导率温变特性引入优化设计模型之中,解决现有高速飞行器多层隔热结构优化设计方法没有考虑材料优选和材料热导率温变影响的问题,得到优化的材料选择和材料几何尺寸,更符合实际应用的需求。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图,其中:
图1是高速飞行器多层隔热结构示意图;
图2 是本发明优化方法流程图;
图3 是采用蒙特卡洛方法优化求解结果;
图4 是采用PSO优化求解结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
如图2所示,本发明是一种高速飞行器多层隔热结构优化设计方法,所述高速飞行器的隔热结构分为多层材料,
确定隔热结构的设计变量,即可选材料数据库和各层材料的几何尺寸;
确定优化目标函数隔热效质比,并设置约束条件;
具体的,所述优化目标函数隔热效质比,即隔热效率和隔热结构的质量比,公式为:
隔热结构的质量为:
根据最大隔热效质比的概念,得到优化目标函数为:
该优化目标函数为质量与效率的比,所以该优化目标函数优化结果越小越好。
具体的,所述约束条件具体为:
所述传热方程:
边界条件为:
其中,为材料热导率,为比热容,为温度,为时间,为空间坐标,为内层材料热导率,为内层结构与空气换热系数,为外层材料热导率,为外层结构与空气的换热系数,为史蒂芬-玻尔兹曼常数,为内部空气温度(图2中的背景温度),为外层外壁面温度,为外部空气温度或热源温度。
采用有限差分隐式格式求解上述偏微分方程(传热方程),将空间区域序列用i表示,时间序列用j表示。空间网格线和时间网格线(i,j)代表了时间-空间的一个节点,及空间位置i处j时刻的状态,用有限差分代替微商,将偏微分方程转化为差分方程:
采用泰勒展开,其导数的有限差分表达式为:
有上面两个式子可知:
关于时间的偏微分的差分为:
将上面的结果带入传热微分方程,可得离散形式的迭代递推格式:
时间上递推1即j变为j+1,有
则有:
上式即为有限差分迭代递推公式。
将优化目标函数、约束条件集合起来,根据设计变量,得到优化模型;
具体的,将所述优化目标函数、约束条件集合起来,得到优化设计数学模型:
设计变量定义为:
则优化模型定义为:
采用智能优化算法求解优化模型,得到优化的设计变量,即优化的材料选择和各层隔热材料的几何尺寸。
具体的,所述设计变量为离散变量和连续变量共存的混合型,是一种混合设计变量式优化问题,采用智能优化算法进行优化模型求解,得到优化的设计变量,即材料选择变量M和三层材料的几何尺寸、和,所述智能优化算法为遗传算法或粒子群算法,但是不仅限于此种优化算法。
实施案例
采用一个算例对本优化方法进行验证和分析。优化模型中的参数取值为:=0.5mm、=3mm、=5mm、=50mm、=2mm、=5mm,=10Gpa,=1200℃,=20℃,=100℃,外层结构与空气的换热系数=230W/(m2·℃),内层结构与空气换热系数=10 W/(m2·℃),=300s。
外层材料数据集如表1:
表1,耐烧蚀层(外层)材料参数(外层)
材料名称 | 密度 (kg·m-3) | 比热容 (J·(kg·K)-1) | 导热系数 (W·(m·K)-1) | 杨氏模量 (GPa) |
C/C复合材料 | 1600 | 713 | 1.6 | 22.32 |
C/SiC复合材料 | 2100 | 1420 | 5 | 49.8 |
ZrB2基复合材料 | 6300 | 430 | 53 | 490 |
内层材料数据集:
表格2,承力层(内层)材料参数
材料名称 | 密度 (kg·m-3) | 比热容 (J·(kg·K)-1) | 导热系数 (W·(m·K)-1) | 杨氏模量 (GPa) |
铝合金7075 | 2800 | 962 | 129.4 | 71 |
GH4099 | 8470 | 624 | 18.9 | 175 |
钛合金TC4 | 4440 | 659 | 9.13 | 115.85 |
中间层材料数据集:
表格3,隔热层(中间层)材料参数
材料名称 | 密度 (kg·m-3) | 比热容 (J·(kg·K)-1) | 导热系数 (W·(m·K)-1) | 杨氏模量 (GPa) |
碳气凝胶(carbon aerogels) | 811 | 870 | 3.521505E-7T2+7.847416E-4T+1.104313 | 2.84 |
SiO2气凝胶 | 140 | 1050 | 1.135885E-10T3-2.256084E-8T2+1.640566E-5T+1.481429E-2 | 0.010 |
SiC气凝胶 | 9.7 | 579 | 5.269242E-11T3-2.221319E-8T2+1.622211E-5T+2.377993E-2 | 0.013 |
采用蒙特卡洛方法优化求解,得到200个可行解优化结果,优化结果如图3所示。
采用PSO优化,种群数量50,迭代次数200,优化结果如图4所示。
PSO优化方法效果优于蒙特卡洛方法,本发明采用PSO优化方法。
根据PSO的优化结果,优化目标函数最优值为10.9167,对应的设计变量为:x=(2,2,1,5.4070e-04,0.0120,0.0020),即外层材料选择C/SiC复合材料,中间层材料选择碳气凝胶,内层材料选择GH4099,外层厚度5.4070e-01mm,中间层厚度12mm,内层厚度2mm。
以上所述,仅为本发明的优选实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域的技术人员在本发明所揭露的技术范围内,可不经过创造性劳动想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书所限定的保护范围为准。
Claims (7)
1.一种高速飞行器多层隔热结构优化设计方法,其特征在于,所述高速飞行器的隔热结构分为多层材料,
确定隔热结构的设计变量,即可选材料数据库和各层材料的几何尺寸;
确定优化目标函数隔热效质比,并设置约束条件;
将优化目标函数、约束条件集合起来,根据设计变量,得到优化模型;
采用智能优化算法求解优化模型,得到优化的设计变量,即优化的材料选择和各层隔热材料的几何尺寸。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310085526.3A CN115809515B (zh) | 2023-02-09 | 2023-02-09 | 一种高速飞行器多层隔热结构优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310085526.3A CN115809515B (zh) | 2023-02-09 | 2023-02-09 | 一种高速飞行器多层隔热结构优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115809515A true CN115809515A (zh) | 2023-03-17 |
CN115809515B CN115809515B (zh) | 2023-05-02 |
Family
ID=85487728
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310085526.3A Active CN115809515B (zh) | 2023-02-09 | 2023-02-09 | 一种高速飞行器多层隔热结构优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115809515B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116052820A (zh) * | 2023-03-22 | 2023-05-02 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种材料热性能评估方法、装置、设备及可读存储介质 |
CN117634313A (zh) * | 2023-12-15 | 2024-03-01 | 浙江阿斯克建材科技股份有限公司 | 一种绝热材料模块化设计方法和*** |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105843996A (zh) * | 2016-03-18 | 2016-08-10 | 中国飞机强度研究所 | 一种防热结构优化方法 |
CN106021734A (zh) * | 2016-05-23 | 2016-10-12 | 北京航空航天大学 | 一种可重复使用飞行器前舱热防护***多目标可靠性优化设计方法 |
WO2016190176A1 (ja) * | 2015-05-22 | 2016-12-01 | 旭硝子株式会社 | 貫通穴付き積層断熱体および断熱構造 |
WO2017072688A1 (en) * | 2015-10-28 | 2017-05-04 | Termokir Industries (1980) Ltd. | Wall profile generator for structures |
CN106874588A (zh) * | 2017-02-08 | 2017-06-20 | 北京航空航天大学 | 一种基于试验设计的多层热防护***非概率不确定性分析和优化设计方法 |
CN106971015A (zh) * | 2016-01-13 | 2017-07-21 | 南京航空航天大学 | 一种空天飞行器热防护***快速分析与设计方法 |
CN212685893U (zh) * | 2020-07-02 | 2021-03-12 | 中南大学 | 高超声速飞行器防隔热一体化结构 |
CN113502144A (zh) * | 2021-08-17 | 2021-10-15 | 北京理工大学 | 一种定向导热隔热材料及其制备方法 |
CN114741915A (zh) * | 2022-03-17 | 2022-07-12 | 上海机电工程研究所 | 针对防隔热一体化结构热防护方案的优化方法及*** |
-
2023
- 2023-02-09 CN CN202310085526.3A patent/CN115809515B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016190176A1 (ja) * | 2015-05-22 | 2016-12-01 | 旭硝子株式会社 | 貫通穴付き積層断熱体および断熱構造 |
WO2017072688A1 (en) * | 2015-10-28 | 2017-05-04 | Termokir Industries (1980) Ltd. | Wall profile generator for structures |
CN106971015A (zh) * | 2016-01-13 | 2017-07-21 | 南京航空航天大学 | 一种空天飞行器热防护***快速分析与设计方法 |
CN105843996A (zh) * | 2016-03-18 | 2016-08-10 | 中国飞机强度研究所 | 一种防热结构优化方法 |
CN106021734A (zh) * | 2016-05-23 | 2016-10-12 | 北京航空航天大学 | 一种可重复使用飞行器前舱热防护***多目标可靠性优化设计方法 |
CN106874588A (zh) * | 2017-02-08 | 2017-06-20 | 北京航空航天大学 | 一种基于试验设计的多层热防护***非概率不确定性分析和优化设计方法 |
CN212685893U (zh) * | 2020-07-02 | 2021-03-12 | 中南大学 | 高超声速飞行器防隔热一体化结构 |
CN113502144A (zh) * | 2021-08-17 | 2021-10-15 | 北京理工大学 | 一种定向导热隔热材料及其制备方法 |
CN114741915A (zh) * | 2022-03-17 | 2022-07-12 | 上海机电工程研究所 | 针对防隔热一体化结构热防护方案的优化方法及*** |
Non-Patent Citations (3)
Title |
---|
李世斌等: "高速飞行器组合式热防护***研究进展" * |
杨明等: "航天器多层绝热材料的非稳态传热分析" * |
柴峻: "临近空间飞行器多层隔热结构传热分析及优化设计" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116052820A (zh) * | 2023-03-22 | 2023-05-02 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种材料热性能评估方法、装置、设备及可读存储介质 |
CN117634313A (zh) * | 2023-12-15 | 2024-03-01 | 浙江阿斯克建材科技股份有限公司 | 一种绝热材料模块化设计方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN115809515B (zh) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115809515A (zh) | 一种高速飞行器多层隔热结构优化设计方法 | |
de Risi et al. | Modelling and optimization of transparent parabolic trough collector based on gas-phase nanofluids | |
CN109960878B (zh) | 一种面向高超声速飞行器总体的主/被动热防护***耦合设计方法 | |
Qian et al. | Topology optimization of a coupled thermal-fluid system under a tangential thermal gradient constraint | |
CN111159939B (zh) | 一种肋片构型拓扑优化***及方法 | |
CN114547790B (zh) | 用于评估复杂多层热防护结构隔热性能的计算方法 | |
CN106570252B (zh) | 一种基于概率技术的热防护***设计方法 | |
Modrek et al. | An optimization case study to design additively manufacturable porous heat sinks based on triply periodic minimal surface (TPMS) lattices | |
CN111967089A (zh) | 一种基于渐变胞元尺寸的异型三维点阵结构隔热性能优化设计方法 | |
CN114781219A (zh) | 一种基于拓扑优化的相变热沉设计方法 | |
Li et al. | Characteristics and optimation of heat pipe radiator for space nuclear propulsion spacecraft | |
CN115577566A (zh) | 多层防热结构接续烧蚀时的处理方法、装置、设备和介质 | |
CN112836399B (zh) | 一种基于有限元算法的非线性接触热阻热分析求解方法 | |
CN109992812A (zh) | 一种圆筒隔热层的传热计算方法 | |
CN112728297B (zh) | 一种基于増材制造提高隔热效果的装置 | |
CN110059372A (zh) | 一种基于差分进化算法的管壳式换热器的多目标设计方法 | |
CN112131801A (zh) | 一种箱式储能***热设计方法 | |
WO2023078475A1 (zh) | 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法 | |
CN117113567A (zh) | 一种自适应板式换热器设计方法 | |
Mekki et al. | Voxel-Based Topology Optimization of Heat Exchanger Fins | |
Ma et al. | Optimization design of the grate cooler based on the power flow method and genetic algorithms | |
Chen et al. | Effect of low wick permeability on transient and steady-state performance of heat pipes | |
CN110600073B (zh) | 一种针对热辐射的热透明装置 | |
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 |
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 |