CN114462137A - 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法 - Google Patents

一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法 Download PDF

Info

Publication number
CN114462137A
CN114462137A CN202111318071.2A CN202111318071A CN114462137A CN 114462137 A CN114462137 A CN 114462137A CN 202111318071 A CN202111318071 A CN 202111318071A CN 114462137 A CN114462137 A CN 114462137A
Authority
CN
China
Prior art keywords
mrve
leading edge
dimensional orthogonal
woven
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
CN202111318071.2A
Other languages
English (en)
Other versions
CN114462137B (zh
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202111318071.2A priority Critical patent/CN114462137B/zh
Publication of CN114462137A publication Critical patent/CN114462137A/zh
Priority to PCT/CN2023/070104 priority patent/WO2023078475A1/zh
Application granted granted Critical
Publication of CN114462137B publication Critical patent/CN114462137B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/26—Composites
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Woven Fabrics (AREA)

Abstract

本发明公开了一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，属于工程热物理技术领域，可以在不增加冷却措施的情况下有效降低前缘温度。本发明首先建立多变量线性回归模型，寻找前缘的理论优化导热配置，使得前缘部件温度得到有效降低；其次建立了适用于三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式，用于指导三维正交机织陶瓷基复合材料的介观机织结构；最后将多变量线性回归模型得到的优化结果与通用公式直接结合，设计出不同区域所需的介观机织结构，实现了前缘宏观温度优化与三维正交机织陶瓷基复合材料介观机织结构的协同设计。

Description

一种基于三维正交机织复合材料的高超音速飞行器前缘热防 护设计方法
技术领域
本发明属于工程热物理技术领域，特别涉及一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法。
背景技术
随着高超声速飞行器以高马赫数飞行时，动能在前缘驻点处转化为内能，使头部区域温度急剧升高(Bull JD,Rasky DJ,Tran HK,et al.Material Response of DiborideMatrix Composites to Low Pressure Simulated Hypersonic Flows.NASA CP-32351994；Part 2(5):653-673.)。由于高超声速飞行器存在严重的气动加热，飞行器前缘的材料必须满足热防护要求。与传统材料相比，陶瓷基复合材料(CMC)具有优异的耐温性和高温力学性能，在高超声速应用中具有显著优势。然而，CMC在超高温(>2000K)环境中的氧化使其力学性能迅速下降。因此，热防护设计对高超音速带CMC前缘的安全稳定工作具有重要意义。
为了实现前沿的热防护设计，研究人员开发了多种基于CMC的热防护结构。在这些热防护结构中，一种是被动热防护结构，以夹层结构为代表，以实现隔热。另一种是主动热保护结构，其表现为引入冷却剂(Ferrari L,Barbato M,Esser B,et al.Sandwichstructured ceramic matrix composites with periodic cellular ceramic cores:anactive cooled thermal protection for space vehicles.Composite Structures2016；154:61-68)。然而，夹层结构的复杂结构给材料制造和构件重量控制带来了挑战。冷却剂的引入将消耗高超声速飞行器的额外能量，因为非空气冷却***需要携带外部冷却剂，这对高超声速飞行器来说是一个负担。因此，迫切需要在不引入额外挑战的情况下实现热防护设计。
除了上述宏观热结构外，CMC中的纤维结构对CMC部件的温度分布有重大影响。一方面，CMC内部的纤维可以编织成不同的结构，这直接决定了CMC的导热性。另一方面，具有不同导热性的各种机织纤维可用于CMC。特别是，一些纤维的超高导热系数大于600W/(mK)。先前的研究(Glass DE.Ceramic Matrix Composite(CMC)Thermal ProtectionSystems(TPS)and Hot Structures for Hypersonic Vehicles.15th AIAAInternational Space Planes and Hypersonic Systems and TechnologiesConference:AIAA；2008)已经证明，具有高导热纤维的合理编织结构将有助于前缘结构散热，避免部件过热，从而避免对高超音速飞行器造成额外负担。因此，采用这种方法进行热防护设计是有利的。
目前大多数的热防护纤维结构设计都是根据停滞点的最高温度直接考虑整个前缘的机织结构。所设计的编织结构只考虑了驻点的最高温度，可以满足热防护设计的要求。然而，基于这种设计，需要大量昂贵的高导热纤维，这最终会增加前沿组件的制造成本。同时，前缘热流分布不均匀，这使得设计对驻点处的散热贡献很大，但对其他区域的散热贡献过大。因此，考虑前缘热流分布不均匀的特点，建立合理的机织结构设计，对提高热防护设计的效率具有重要意义。
为了找到合理的前缘机织结构设计，需要确定具有不同机织结构的CMC的热导率值。目前，有一种方法可以估算三维正交编织CMC的热导率([18]
Lee SE,Yoo JS,Kang JH,et al.Prediction of the thermal conductivitiesof four-axial non-woven composites.Composite Structures 2009；89(2):262-269.)，但它适用于特定的三维正交编织结构。如果机织结构发生变化，则很难通过调整导热系数计算来获得相应结构的导热系数。因此，迫切需要建立一种有效的通用三维正交机织CMC导热系数计算方法。
发明内容
本发明针对高超音速前缘部件热防护的需求，考虑到目前热防护设计中冷却剂的引入将消耗高超声速飞行器的额外能，本发明提出了一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，该方法通过三维正交机织陶瓷基复合材料前缘部件的宏观温度场优化和介观机织结构的协同设计来降低了前缘部件的温度，可以在不增加冷却措施的情况下有效降低前缘温度。首先，建立多变量线性回归模型，寻找前缘的理论优化导热配置，使得前缘部件温度得到有效降低。其次，建立了适用于三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式，用于指导三维正交机织陶瓷基复合材料的介观机织结构。最后，将多变量线性回归模型得到的优化结果与通用公式直接结合，设计出不同区域所需的介观机织结构，实现了前沿宏观温度优化与三维正交机织陶瓷基复合材料介观机织结构的协同设计。
为实现上述目的，本发明采用的技术方案为：
一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，包括以下步骤：
步骤一：将前缘模型导入Comsol Multiphysics软件，对前缘模型进行各向异性热物性参数赋值并进行网格划分，对前缘模型外表面给热流密度边界条件，底部给定温边界条件开展温度场的有限元计算，获得初始条件下前缘部件的温度场以及热流密度场；
步骤二：根据步骤一中获取的温度场以及热流密度场特征，将前缘结构划分为不同的区域，得到前缘热疏导模型；
步骤三：将步骤二得到的前缘热疏导模型导入Comsol Multiphysics软件，对所述前缘热疏导模型不同区域进行热物性参数赋值，并重复步骤一的操作，同时使得热物性参数在一定范围内变化，即可得到一组前缘温度与对应温度下区域热物性之间的数据样本；
步骤四：根据步骤三得到的前缘温度与各区域热物性之间的数据样本，通过线性回归的方法，建立起前缘温度多变量线性回归优化模型，通过该模型可以得到在一定导热系数范围内，每个区域导热系数该如何配置，才能使得前缘部件温度最低；
步骤五：根据三维正交机织材料的介观机织结构特征，建立一个适用于三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式；
步骤六：根据步骤四中所述多变量线性回归优化模型得到的最优导热系数配置，几何步骤五的导热系数通用计算公式，即可得到不同区域的材料机织介观结构，实现前缘热防护设计和材料介观结构的一体化协同设计，以降低前缘部件的温度；
所述步骤四中，前缘温度多变量线性回归优化模型如下式：
Figure BDA0003342531670000031
Figure BDA0003342531670000032
ε＝(ε1…εn)T
β＝(β0 β1…βm)T＝(β0…βxx-i βyy-i βzz-i…)T,i＝1,2,…,6
其中T表示前缘最高温度，β表示各变量的回归系数，ε为残差，En为n介单位矩阵，X表示变量值矩阵，σ2为平方差。
表1各变量回归系数值
Figure BDA0003342531670000041
所述步骤五中，三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式，如下式：
Figure BDA0003342531670000042
vmrve,m＝1-vmrve,f
λmrve,a＝λmvmrve,m+λf，avmrve,f
Figure BDA0003342531670000043
λp1,x＝λmrve,a
λp2,x＝λmrve,t
λp3,x＝λm
Figure BDA0003342531670000044
λp1,y＝λmrve,t
λp2,y＝λmrve,a
λp3,y＝λm
Figure BDA0003342531670000051
Figure BDA0003342531670000052
λp1,z＝λmrve,t
λp2,z＝λmrve,t
λp3,z＝λm
Figure BDA0003342531670000053
Figure BDA0003342531670000054
λp4,x＝λmrve,a
λp5,x＝λmrve,t
λp6,x＝λm
Figure BDA0003342531670000055
λp4,y＝λmrve,t
λp5,y＝λmrve,t
λp6,y＝λm
Figure BDA0003342531670000056
Figure BDA0003342531670000057
λp4,z＝λmrve,t
λp5,z＝λmrve,a
λp6,z＝λm
Figure BDA0003342531670000061
Figure BDA0003342531670000062
Figure BDA0003342531670000063
Figure BDA0003342531670000064
Figure BDA0003342531670000065
其中：D为纤维直径，L为代表单元的长度。代表单元X方向上的纤维数目为：Nxy×Nxz,代表单元Y方向上的纤维数目为：Nyx×Nyz,代表单元Z方向上的纤维数目为：Nzx×Nzy，vmrve,f表示单胞单元中纤维的体积分数，vmrve,m表示单胞单元中基体的体积分数，λmrve,a表示单胞轴向的等效导热系数，vmrve,m表示单胞的横向等效导热系数，λm表示纤维的基体的导热系数，λp1,x、λp2,x、λp3,x、λp1,y、λp2,y、λp3,y、λp1,z、λp2,z、λp3,z分别表示典型代表单元中层I区域p1、p2、p3沿X、Y、Z3个方向的等效导热系数，λI,x、λI,y、λI,z分别表示典型代表单元中层I沿X、Y、Z3个方向的等效导热系数，λp4,x、λp5,x、λp5,x、λp4,y、λp5,y、λp6,y、λp4,z、λp5,z、λp6,z分别表示典型代表单元中层II区域p4、p5、p6沿X、Y、Z3个方向的等效导热系数，λII,x、λII,y、λII,z分别表示典型代表单元中层II沿X、Y、Z3个方向的等效导热系数，λRVE,x为三维正交材料沿X方向的导热系数，λRVE,y为三维正交材料沿Y方向的导热系数，λRVE,z为三维正交材料沿Z方向的导热系数。
有益效果：本发明提供了一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法；首先建立多变量线性回归模型，寻找前缘的理论优化导热配置，使得前缘部件温度得到有效降低；其次建立了适用于三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式，用于指导三维正交机织陶瓷基复合材料的介观机织结构；最后将多变量线性回归模型得到的优化结果与通用公式直接结合，设计出不同区域所需的介观机织结构，实现了前沿宏观温度优化与三维正交机织陶瓷基复合材料介观机织结构的协同设计。本发明通过三维正交机织陶瓷基复合材料前缘部件的宏观温度场优化和介观机织结构的协同设计来降低了前缘部件的温度，能够可以在不增加冷却措施的情况下有效降低前缘温度，在热防护设计的同时实现了对前缘部件材料的设计。
附图说明
图1为前缘模型及边界条件示意图；
图2为前缘未温度场优化前材料结构以及前缘分区示意图。其中：(a)前缘温度场未优化前材料结构示意图，(b)前缘分区示意图；
图3为前缘温度场示意图，(a)前缘未优化温度场示意图，(b)前缘优化后温度场示意图；
图4为多元线性回归温度场优化示意图；
图5为三维机织复合材料机织结构示意图；
图6为优化后前缘各区域机织结构示意图；
图7为优化前后前缘特征线上温度分布曲线。
具体实施方式
下面结合附图和具体实施例对本发明作更进一步的说明。
实施例1
以某高超音速飞行器前缘为例阐述基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，如图1所示，前缘结构的几何特征以及计算边界条件；前缘主要由两部分组成：激波针和肩部；激波针驻点区域为半球形，直径Ds＝50mm，冲击波针的长度Ls＝75mm，前缘的长度Ll＝300mm，前缘宽度Wl＝250mm，前缘高度Hl＝85mm，驻点区域具有最强的气动加热，热通量5000KW/m2，肩部区域的热通量为3000KW/m2，前缘底部设置为恒温边界，T＝293.15K，其他表面的热通量为500W/m2，同时，考虑了前缘表面对环境的热辐射，表面发射率ε为0.9，环境温度为T∞是293.15K。
未优化前，整个前缘结构都使用的同一种三维正交机织复合材料，如图2(a)所示。在该机制结构下材料3个方向上的导热系数λx＝λν＝λz＝176W/(m·K)。前缘温度场计算主要采用Comsol Multiphysics软件中的稳态固体传热模块，通过计算可以得到未优化前，前缘温度场分布，如图3(a)所示，前缘的头部和两肩出现高温，尤其在头部区域温度高达2800K。
根据未优化前的前缘温度场分布，在前缘温度场优化时将前缘划分为如图2(b)所示的6个区域，在优化时每个区域的导热系数λx、λν、λz可以在100W/(m·K)达150W/(m·K)之间变化，将导热系数变量带入到如下的多变量线性回归温度优化模型中，即可得到每个区域的最优导热系数配置，在该配置下前缘具有最低温度。
Figure BDA0003342531670000081
Figure BDA0003342531670000082
ε＝(ε1…εn)T
β＝(β0 β1…βm)T＝(β0…βxx-i βyy-i βzz-i…)T,i＝1,2,…,6
表1各变量回归系数值
Figure BDA0003342531670000083
其中T表示前缘最高温度，β表示各变量的回归系数，ε为残差，En为n介单位矩阵，X表示变量值矩阵，σ2为平方差。
多变量线性回归温度优化模型下每个区域导热系数变化与前缘温度的关系如图4所示，多变量线性回归温度优化模型给出的每个区域导热系数配置如表2所示：
表2每个区域导热系数配置
Figure BDA0003342531670000091
三维正交编织复合材料的机织结构示意图，如图5所示。材料由三个方向上的机织纤维和基体材料组成。为了制造出满足各个区域所需导热系数的三维正交复合材料，则需要调整材料三个方向上的纤维数目比。将优化模型得到的导热系数带入到本发明建立的适用于三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式，如下式：即可得到每个区域材料在三个编织方向上的纤维数目，如表2所示。
Figure BDA0003342531670000092
vmrve,m＝1-vmrve,f
λmrve,a＝λmvmrve,m+λf，avmrve,f
Figure BDA0003342531670000093
λp1,x＝λmrve,a
λp2,x＝λmrve,t
λp3,x＝λm
Figure BDA0003342531670000094
λp1,y＝λmrve,t
λp2,y＝λmrve,a
λp3,y＝λm
Figure BDA0003342531670000095
Figure BDA0003342531670000096
λp1,z＝λmrve,t
λp2,z＝λmrve,t
λp3,z＝λm
Figure BDA0003342531670000101
Figure BDA0003342531670000102
λp4,x＝λmrve,a
λp5,x＝λmrve,t
λp6,x＝λm
Figure BDA0003342531670000103
λp4,y＝λmrve,t
λp5,y＝λmrve,t
λp6,y＝λm
Figure BDA0003342531670000104
Figure BDA0003342531670000105
λp4,z＝λmrve,t
λp5,z＝λmrve,a
λp6,z＝λm
Figure BDA0003342531670000106
Figure BDA0003342531670000107
Figure BDA0003342531670000108
Figure BDA0003342531670000109
Figure BDA00033425316700001010
其中：λRVE,x为三维正交材料沿X方向的导热系数，λRVE,y为三维正交材料沿Y方向的导热系数，λRVE,z为三维正交材料沿Z方向的导热系数。
综上即实现前缘宏观温度优化与三维正交机织陶瓷基复合材料介观机织结构的协同设计，整个前缘在传热优化后不同区域的材料机织结构如图6所示，其中区域1为前缘激波针的半球区，区域2为前缘激波针的圆柱区域，区域3为激波针与肩部的连接区域，区域4为前缘两侧的肩部区域，区域5为两肩部区域的中间区域，区域6为前缘的底部区域。从图7前缘特征线上的温度分布曲线可以看出，基于本发明方法可以在不引入冷却措施的情况下有效降低前缘部件的温度，尤其在头部区域，最多可以减低前缘约300K左右的温度。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (4)

1.一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，其特征在于，包括以下步骤：
步骤一：将前缘模型导入Comsol Multiphysics软件，对前缘模型进行各向异性热物性参数赋值并进行网格划分，对前缘模型外表面给热流密度边界条件，底部给定温边界条件开展温度场的有限元计算，获得初始条件下前缘部件的温度场以及热流密度场；
步骤二：根据步骤一中获取的温度场以及热流密度场特征，将前缘结构划分为不同的区域，得到前缘热疏导模型；
步骤三：将步骤二得到的前缘热疏导模型导入Comsol Multiphysics软件，对所述前缘热疏导模型不同区域进行热物性参数赋值，并重复步骤一的操作，同时使得热物性参数在一定范围内变化，即可得到一组前缘温度与对应温度下区域热物性之间的数据样本；
步骤四：根据步骤三得到的前缘温度与各区域热物性之间的数据样本，通过线性回归的方法，建立起前缘温度多变量线性回归优化模型，通过该模型可以得到在一定导热系数范围内，每个区域导热系数该如何配置，才能使得前缘部件温度最低；
步骤五：根据三维正交机织材料的介观机织结构特征，建立一个适用于三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式；
步骤六：根据步骤四中所述多变量线性回归优化模型得到的最优导热系数配置，几何步骤五的导热系数通用计算公式，即可得到不同区域的材料机织介观结构，实现前缘热防护设计和材料介观结构的一体化协同设计，以降低前缘部件的温度。
2.根据权利要求1所述的基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，其特征在于，步骤四中前缘温度多变量线性回归优化模型如下式：
Figure FDA0003342531660000011
Figure FDA0003342531660000012
ε＝(ε1…εn)T
β＝(β0 β1…βm)T＝(β0…βxx-i βyy-i βzz-i…)T,i＝1,2,…,6
其中T表示前缘最高温度，β表示各变量的回归系数，ε为残差，En为n介单位矩阵，X表示变量值矩阵，σ2为平方差。
3.根据权利要求2所述的基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，其特征在于，所述各变量的回归系数值为：
Figure FDA0003342531660000021
4.根据权利要求1所述的基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法，其特征在于，步骤五中三维正交机织陶瓷基复合材料不同机织结构各向异性导热系数计算的通用计算公式，如下式：
Figure FDA0003342531660000022
vmrve,m＝1-vmrve,f
λmrve,a＝λmvmrve,m+λf，avmrve,f
Figure FDA0003342531660000023
λp1,x＝λmrve,a
λp2,x＝λmrve,t
λp3,x＝λm
Figure FDA0003342531660000024
λp1,y＝λmrve,t
λp2,y＝λmrve,a
λp3,y＝λm
Figure FDA0003342531660000031
Figure FDA0003342531660000032
λp1,z＝λmrve,t
λp2,z＝λmrve,t
λp3,z＝λm
Figure FDA0003342531660000033
Figure FDA0003342531660000034
λp4,x＝λmrve,a
λp5,x＝λmrve,t
λp6,x＝λm
Figure FDA0003342531660000035
λp4,y＝λmrve,t
λp5,y＝λmrve,t
λp6,y＝λm
Figure FDA0003342531660000036
Figure FDA0003342531660000037
λp4,z＝λmrve,t
λp5,z＝λmrve,a
λp6,z＝λm
Figure FDA0003342531660000041
Figure FDA0003342531660000042
Figure FDA0003342531660000043
Figure FDA0003342531660000044
Figure FDA0003342531660000045
其中：D为纤维直径，L为代表单元的长度。代表单元X方向上的纤维数目为：Nxy×Nxz,代表单元Y方向上的纤维数目为：Nyx×Nyz,代表单元Z方向上的纤维数目为：Nzx×Nzy，vmrve,f表示单胞单元中纤维的体积分数，vmrve,m表示单胞单元中基体的体积分数，λmrve,a表示单胞轴向的等效导热系数，vmrve,m表示单胞的横向等效导热系数，λm表示纤维的基体的导热系数，λp1,x、λp2,x、λp3,x、λp1,y、λp2,y、λp3,y、λp1,z、λp2,z、λp3,z分别表示典型代表单元中层I区域p1、p2、p3沿X、Y、Z3个方向的等效导热系数，λI,x、λI,y、λI,z分别表示典型代表单元中层I沿X、Y、Z3个方向的等效导热系数，λp4,x、λp5,x、λp5,x、λp4,y、λp5,y、λp6,y、λp4,z、λp5,z、λp6,z分别表示典型代表单元中层II区域p4、p5、p6沿X、Y、Z3个方向的等效导热系数，λII,x、λII,y、λII,z分别表示典型代表单元中层II沿X、Y、Z3个方向的等效导热系数，λRVE,x为三维正交材料沿X方向的导热系数，λRVE,y为三维正交材料沿Y方向的导热系数，λRVE,z为三维正交材料沿Z方向的导热系数。
CN202111318071.2A 2021-11-08 2021-11-08 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法 Active CN114462137B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202111318071.2A CN114462137B (zh) 2021-11-08 2021-11-08 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法
PCT/CN2023/070104 WO2023078475A1 (zh) 2021-11-08 2023-01-03 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111318071.2A CN114462137B (zh) 2021-11-08 2021-11-08 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法

Publications (2)

Publication Number Publication Date
CN114462137A true CN114462137A (zh) 2022-05-10
CN114462137B CN114462137B (zh) 2024-07-23

Family

ID=81406732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111318071.2A Active CN114462137B (zh) 2021-11-08 2021-11-08 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法

Country Status (2)

Country Link
CN (1) CN114462137B (zh)
WO (1) WO2023078475A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023078475A1 (zh) * 2021-11-08 2023-05-11 南京航空航天大学 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106250575A (zh) * 2016-05-19 2016-12-21 南京航空航天大学 一种考虑纤维走向的编织复合材料涡轮叶片温度场计算方法
CN107792392A (zh) * 2017-09-26 2018-03-13 北京航天长征飞行器研究所 一种飞行器控制舵尖化前缘主动式热疏导试验***及方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020211012A1 (zh) * 2019-04-17 2020-10-22 大连理工大学 一种面向混杂纤维复合材料板壳结构的快速协同优化方法
CN110083900B (zh) * 2019-04-17 2023-02-10 大连理工大学 一种面向混杂纤维复合材料板壳结构的快速协同优化方法
CN110909495A (zh) * 2019-10-31 2020-03-24 南京航空航天大学 基于全尺寸微观结构模型的编织cmc材料薄壁构件等效导热系数预估方法
CN112149236A (zh) * 2020-10-12 2020-12-29 南京航空航天大学 考虑微观编织结构及纤维束方向变化的陶瓷基复合材料涡轮叶片热分析方法
CN114462137B (zh) * 2021-11-08 2024-07-23 南京航空航天大学 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106250575A (zh) * 2016-05-19 2016-12-21 南京航空航天大学 一种考虑纤维走向的编织复合材料涡轮叶片温度场计算方法
CN107792392A (zh) * 2017-09-26 2018-03-13 北京航天长征飞行器研究所 一种飞行器控制舵尖化前缘主动式热疏导试验***及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENWEI ZHAO等: ""Multi-scale collaborative design method for macroscopic thermal optimization and mesoscopic woven structure of hypersonic vehicle’s TOCMC leading edge"", 《CHINESE JOURNAL OF AERONAUTICS》, vol. 37, no. 4, 30 April 2024 (2024-04-30), pages 524 - 541, XP087493952, DOI: 10.1016/j.cja.2024.01.003 *
穆军武;边天涯;唐斐;董志兴;: "高超声速飞行器热防护***分析与数值计算", 科技导报, no. 05, 13 March 2015 (2015-03-13) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023078475A1 (zh) * 2021-11-08 2023-05-11 南京航空航天大学 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法

Also Published As

Publication number Publication date
CN114462137B (zh) 2024-07-23
WO2023078475A1 (zh) 2023-05-11

Similar Documents

Publication Publication Date Title
Shittu et al. High performance and thermal stress analysis of a segmented annular thermoelectric generator
WO2023078475A1 (zh) 一种基于三维正交机织复合材料的高超音速飞行器前缘热防护设计方法
CN108101568B (zh) 一种高导热碳/碳复合材料及其制备方法
CN104553102A (zh) 超高温梯度隔热材料及其制备方法
WO2009061492A1 (en) Composite material compositions, arrangements and methods having enhanced thermal conductivity behavior
CN104876580B (zh) 一种轻质高导热碳基材料的制备方法
CN112918025B (zh) 一种烧蚀型防隔热吸波一体化复合材料及其制备方法
CN109960878A (zh) 一种面向高超声速飞行器总体的主/被动热防护***耦合设计方法
Chen et al. Design and heat transfer analysis of a compound multi-layer insulations for use in high temperature cylinder thermal protection systems
CN112149236A (zh) 考虑微观编织结构及纤维束方向变化的陶瓷基复合材料涡轮叶片热分析方法
CN112851313A (zh) 一种高温热敏电阻材料及其微波制备方法
CA2046560C (en) Cooled refractory structure and manufacturing process therefor
Hemrick et al. Woven graphite fiber structures for use in ultra-light weight heat exchangers
CN112251906A (zh) 一种适用于陶瓷基复合材料构件的气膜冷却孔编织结构及其制作方法
GB2340595A (en) Composite heat exchanger
CN112668113A (zh) 一种复合材料多尺度防热优化方法
CN116976007A (zh) 一种基于变换热学的飞行器定向疏导热防护层设计方法
CN112728297A (zh) 一种基于増材制造提高隔热效果的装置
EP2298716A2 (en) Methods of rapidly densifying complex-shaped, asymmetrical porous structures
CN104080309A (zh) 具有散热效果的复合材料纤维布及其制造方法
EP1254346A1 (en) Composite heat exchanger having strengthened joints
CN109152298B (zh) 一种轻质温控装置及其制造方法
Wang et al. Lightweight chopped carbon fiber/carbon composites with low thermal conductivity fabricated by vacuum filtration method
CN111692612A (zh) 一种三维机织航空火焰稳定器及其制备方法
CN109736972B (zh) 火箭底部热防护面板、液体火箭底部热防护***及方法

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