CN106126791B - 一种考虑几何不确定性的高超声速机翼气动力/热分析方法 - Google Patents
一种考虑几何不确定性的高超声速机翼气动力/热分析方法 Download PDFInfo
- Publication number
- CN106126791B CN106126791B CN201610438739.XA CN201610438739A CN106126791B CN 106126791 B CN106126791 B CN 106126791B CN 201610438739 A CN201610438739 A CN 201610438739A CN 106126791 B CN106126791 B CN 106126791B
- Authority
- CN
- China
- Prior art keywords
- airfoil
- aerodynamic
- interval
- wing
- hypersonic
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title abstract description 6
- 238000000034 method Methods 0.000 claims abstract description 51
- 230000004044 response Effects 0.000 claims abstract description 50
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 12
- 238000013461 design Methods 0.000 claims description 16
- 238000002076 thermal analysis method Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 7
- 238000003754 machining Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000010438 heat treatment Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000013598 vector Substances 0.000 abstract description 8
- 238000013459 approach Methods 0.000 abstract description 2
- 230000004907 flux Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种考虑几何不确定性的高超声速机翼气动力/热分析方法。根据伯恩斯坦多项式近似模型构造不确定参数的样本点，在给定飞行工况下，利用气动力/热工程算法计算样本点对应机翼外形的气动力/热响应值，同时应用最小二乘法得到多项式近似模型的拟合系数。在此基础上，利用已建立的多项式近似模型求解气动力/热响应值关于各个不确定参数的极大/小值点，并组合形成最大/小值点的向量，最终得到气动力/热响应值的区间上界和区间下界，实现考虑几何不确定性的高超声速机翼气动力/热分析。本发明方法得到的区间边界与蒙特卡洛方法得到的区间边界吻合较好，并且可以实现区间包络，为高超声速机翼外形的总体设计提供了新思路。
Description
技术领域
本发明涉及高超声速飞行器设计技术领域，特别涉及一种考虑几何不确定性的高超声速机翼气动力/热分析方法。
背景技术
气动外形建模方法是高超声速飞行器设计与研制过程中的重要环节。为了提高建模效率和精度，一般需要对设计参数的数目进行控制，并且采用参数化建模方法。以机翼外形的几何建模方法为例，传统的参数化几何建模主要依赖于CAD方法(B样条曲线、NURBS曲线和二次曲线方法)，但该建模方法由于表达参数偏多，不利于概念设计中各学科分析模型的简化。相比而言，基于类函数/形函数转换的CST(Class function and Shape functionTransformation)方法通过修改少量的设计变量，可以生成翼型和机翼几何外形，而且具有快速、准确的拟合效果。
然而，由于加工工艺、制造水平的限制，高超声速飞行器的气动外形存在由加工误差导致的几何不确定性，几何设计参数往往在名义值附近上下波动。对于由加工误差引起的几何不确定性，现有的研究常采用概率方法对不确定因素进行定量化表征，利用标准正态函数描述几何参数的分布特征。但是，当没有足够的数据来验证这些随机变量概率密度的正确性时，概率方法难以可靠地满足精度要求的计算结果。与概率方法相比，非概率区间方法仅需明确不确定参数的分布界限，能够在不确定参数概率密度未知的情况下，利用区间向量对不确定参数进行定量化。非概率区间方法在结构的静、动力特性分析领域已经取得了一定成果，但在高超声速飞行器机翼气动外形建模中还处于起步阶段，相关的研究成果十分有限。同时，当机翼气动外形存在不确定性时，尤其是机翼前缘位置出现的几何波动，会造成机翼气动力/气动热性能存在偏差，这种波动或偏差对飞行安全是不利的。综上所述，亟需发展一种考虑几何不确定性的高超声速机翼气动力/热非概率区间分析方法，以克服传统概率方法对参数大样本容量试验数据的依赖性。
发明内容
本发明要解决的技术问题为：针对高超声速机翼在加工、制造过程中由加工误差所导致的几何不确定性，提出一种考虑几何不确定性的高超声速机翼气动力/热分析方法。该方法在CST参数化建模技术基础上，利用区间向量对机翼几何特征参数进行定量化表征。在给定工况下，利用气动力/热工程算法计算样本点处的气动力/热响应值，通过伯恩斯坦多项式近似模型及最小二乘法建立气动力/热响应值与不确定参数之间的映射关系。在给定区间范围内，利用已建立的多项式近似模型求解气动力/热响应值关于各个不确定参数最大/小值点的向量，计算得到气动力/热响应值的区间上界和区间下界，实现考虑几何不确定性的高超声速机翼气动力/热分析。
本发明解决上述技术问题采用的技术方案为：一种考虑几何不确定性的高超声速机翼气动力/热分析方法，包括以下步骤：
步骤(1)、确定高超声速机翼的平面轮廓参数，包括翼根弦长Cr、翼尖弦长Ct、副翼宽度Ca及机翼半展长L；
步骤(2)、利用类函数/形函数转换方法建立高超声速机翼翼型的参数化表达式，翼型的几何曲线可用下列函数表示：
式中，x/c为翼型弦向的无量纲坐标值，y/c为翼型法向的无量纲坐标值，C(x/c)和S(x/c)分别为类型函数和形状函数，zte/c为翼型后缘点的无量纲坐标值，C(x/c)可表示为：
对于Clark-Ys翼型，取指数N1＝0.5,N2＝1，S(x/c)可表示为：
S(x/c)与翼型前缘半径Rle/c及后缘倾角β满足以下关系：
在形状函数S(x/c)中取n＝3，这样共包含等8个几何设计参数，其中，表示翼型前缘半径，β1和β2表示上、下翼面的后缘倾角，zte/c为翼型后缘点的无量纲坐标值，λ1、λ2、λ1'、λ2'为上、下翼面形状函数多项式的加权系数；
步骤(3)、考虑由加工误差导致翼型几何特征参数存在的波动，将存在不确定性的几何特征参数记为其余参数均视为确定性参数；
步骤(4)、设定翼型几何特征参数α的区间上下界，可表示为：
式中，为高超声速机翼翼型几何特征参数的区间下界，为区间上界；
步骤(5)、利用区间[-1,1]内的n阶伯恩斯坦多项式，构造用于计算多项式系数的样本点，将伯恩斯坦多项式模型中子多项式的极值点t记为：
t＝[t0,t1,…,tr,…,tn-1,tn]T (6)
式中：
步骤(6)、将式(6)所表示的极值点映射至区间参数空间αI内的第s个区间参数αs，则在第s维用于计算该维近似多项式系数的样本点可以表示为：
其中，
式中，m为不确定参数的个数；
步骤(7)、根据式(8)中样本点对应翼型的几何特征参数，建立高超声速机翼的气动外形，进而利用高超声速工程算法计算气动力/热响应值，记为Qs，表示为如下形式：
Qs＝[Qs,0,...,Qs,r,...,Qs,n]T (10)
式中，Qs,r＝Q(αs,r)表示样本点αs,r对应的气动力/热响应值；
步骤(8)、利用最小二乘法计算伯恩斯坦多项式近似模型的待定系数ωs，公式如下：
ωs＝[ωs,0,ωs,1,...,ωs,n]T＝(BTB)-1.BT·Qs (11)
式中，B为伯恩斯坦-范德蒙矩阵，可表示为：
其中，bn,r(t)为伯恩斯坦多项式的子多项式，记为：
步骤(9)、根据步骤(8)中得到的系数ωs，建立气动力/热响应值Q关于第s个不确定参数的伯恩斯坦多项式近似模型Ps(t)，表示如下：
步骤(10)、步骤(9)中Ps(t)的最值点可以通过其导函数的零点与自变量的端点产生，其中导函数零点可通过下式计算：
将零点记为ts＝[ts,1,ts,2,...,ts,n]，根据Ps(t)的自变量范围，可将零点进一步修改为：
其中Im和Re分别表示对应变量的虚部和实部；
步骤(11)、气动力/热响应值Q关于第s个区间参数的最小值点和最大值点可计算为：
式中，Ts＝[ts,-1,1]；
步骤(12)、重复步骤(6)～(11)，以相同的方式计算气动力/热响应值Q关于所有区间参数的最小值点和最大值点，最终可以得到气动力/热响应值在标准区间[-1,1]内的最小值点和最大值点，即：
步骤(13)、将步骤(12)中标准区间内的最小值点和最大值点映射到实际参数空间有：
式中：
步骤(14)、根据αmin和αmax建立其对应的机翼参数化外形，利用高超声速工程算法得到气动力/热响应值Q的区间上下界，可表示为：
式中，Q(αmin)和Q(αmax)分别表示几何特征参数αmin和αmax对应的气动力/热响应值；
步骤(15)、将通过上述方法得到的气动力/热响应值Q的区间上下界与蒙特卡洛方法进行对比，验证所建立方法的有效性。
其中，所述步骤(1)中，机翼平面轮廓参数由表1确定。
表1 机翼平面轮廓参数
其中，所述步骤(2)中，所采用的Clark-Ys标准翼型对应的设计参数由表2确定。
表2 Clark-Ys翼型设计参数
其中，所述步骤(4)中，含不确定性的翼型几何特征参数的区间上下界由表3确定。
表3 翼型几何特征参数的区间界限
其中，所述步骤(7)中，工程算法的计算流程为：首先根据机翼的外形参数建立几何模型，对几何模型进行非结构表面网格划分，然后分别利用切劈法和达黑姆巴克法计算机翼迎风面和背风面的表面压力系数。在此基础上，利用费雷德经验公式计算机翼驻点加热，从而得到机翼表面的最大热流密度、升力系数、阻力系数等三个气动力/热响应值。
其中，所述步骤(15)中，蒙特卡洛方法的计算流程为：在由不确定参数组成的区间范围内，随机选择1000个样本点，逐次计算所有样本点对应的气动力/热响应值，并将响应值中的最大值和最小值作为区间上界和区间下界。
本发明的有益效果是：
本发明利用CST参数化建模方法实现机翼几何外形的快速表征，同时考虑机翼加工过程中由加工误差导致的几何不确定性，建立含区间参数的高超声速机翼气动外形。利用伯恩斯坦多项式近似模型和最小二乘法，建立气动力/热响应值与不确定参数之间的函数关系。通过求解气动力/热响应值关于各个不确定参数的最大/小值点的向量，最终得到气动力/热响应值的区间上界和区间下界，实现考虑几何不确定性的高超声速机翼气动力/热分析。本发明方法得到的区间边界与蒙特卡洛方法得到的区间边界吻合较好，并且可以实现区间包络，为高超声速机翼外形的总体设计提供了新思路。
附图说明
图1为高超声速机翼平面轮廓参数示意图；
图2为Clark-Ys标准翼型示意图；
图3为翼型几何特征参数示意图；
图4为高超声速机翼气动外形示意图；
图5为高超声速机翼表面非结构网格示意图；
图6为机翼表面最大热流密度随马赫数的变化曲线图；
图7为机翼升力系数随马赫数的变化曲线图；
图8为机翼阻力系数随马赫数的变化曲线图；
图9为本发明的方法实现流程图。
具体实施方式
以下将参照附图，对本发明的设计实例进行详细描述。应当理解，所选实例仅为了说明本发明，而不是限制本发明的保护范围。
(1)以类X-37B高超声速飞行器机翼作为研究对象，首先确定高超声速机翼的平面轮廓参数，包括翼根弦长Cr、翼尖弦长Ct、副翼宽度Ca及机翼半展长L。机翼平面轮廓如图1所示，相关参数取值见表4。
表4 机翼平面轮廓参数
(2)以Clark-Ys翼型作为标准翼型(见图2)，利用类型函数/形状函数转换CST方法设定参数化翼型曲线的表达式，其中部分参数的几何含义如图3所示，几何参数由表5给出。
表5 Clark-Ys翼型设计参数
(3)考虑由加工误差导致翼型几何特征参数存在的波动，设定几何特征参数的区间上下界，如表6所示。
表6 翼型几何特征参数的区间界限
(4)利用区间[-1,1]内的n阶伯恩斯坦多项式，构造用于计算多项式系数的样本点，将伯恩斯坦多项式模型中子多项式的极值点t记为：
t＝[t0,t1,…,tr,…,tn-1,tn]T (22)
式中：
本实例中，采用10阶伯恩斯坦多项式，故这里取n＝10。
(5)将式(22)所表示的极值点映射至区间参数空间αI内的第s个区间参数αs，则在第s维用于计算该维近似多项式系数的样本点可以表示为：
其中，
本实例中，共包含3个不确定参数，故这里取m＝3。
(6)根据式(24)中各个样本点αs,r对应的翼型几何特征参数，利用商业软件CATIA的参数化建模功能建立高超声速机翼的气动外形，如图4所示。将气动外形的几何文件导入商业软件ICEM中进行机翼表面非结构网格的自由划分，如图5所示。得到表面网格节点信息后，在给定工况条件下(如表7所示)，分别利用切劈法和达黑姆巴克法计算机翼迎风面和背风面的表面压力系数。在此基础上，利用费雷德经验公式计算机翼驻点加热情况，从而得到机翼表面的最大热流密度qmax、升力系数Cl、阻力系数Cd等三个气动力/热响应值。同时，将气动力/热响应值以向量Qs表示为如下形式：
Qs＝[Qs,0,...,Qs,r,...,Qs,n]T (26)
式中，Qs,r＝Q(αs,r)表示样本点αs,r对应的气动力/热响应值。
表7 飞行工况
注：表7中，H表示飞行高度，T∞表示来流温度，ρ∞表示来流密度，p∞表示压强，alpha表示攻角，Ma表示飞行马赫数。
(7)根据式(22)中极值点分布，计算得到伯恩斯坦-范德蒙矩阵B，如下所示：
其中，bn,r(t)为伯恩斯坦多项式的子多项式，记为：
(8)利用最小二乘法计算伯恩斯坦多项式近似模型的待定系数ωs，公式如下：
ωs＝[ωs,0,ωs,1,...,ωs,n]T＝(BTB)-1.BT·Qs (29)
由此得到气动力/热响应值Q关于第s个不确定参数的伯恩斯坦多项式近似模型Ps(t)，表示如下：
(9)Ps(t)的最值点可以通过其导函数的零点与自变量的端点产生，其中导函数零点可通过下式计算：
将零点记为ts＝[ts,1,ts,2,...,ts,n]，根据Ps(t)的自变量范围，可将零点进一步修改为：
其中Im和Re分别表示对应变量的虚部和实部。
(10)气动力/热响应值Q关于第s个区间参数的最小值点和最大值点可计算为：
式中，Ts＝[ts,-1,1]。
(11)重复步骤(5)～(10)，以相同的方式计算气动力/热响应值Q关于所有区间参数的最小值点和最大值点，最终可以得到气动力/热响应值在标准区间[-1,1]内的最小值点和最大值点，即：
(12)将步骤(11)中的标准区间内的最小值点和最大值点映射到实际参数空间有：
式中：
(13)根据αmin和αmax建立其对应的机翼参数化外形，利用高超声速工程算法得到气动力/热响应值Q的区间上下界，表示为：
式中，Q(αmin)和Q(αmax)分别表示几何参数αmin和αmax对应的气动力/热响应值。
(14)将通过上述方法得到的气动力/热响应值Q的区间上下界与蒙特卡洛方法进行对比，验证所建立方法的有效性。蒙特卡洛方法的计算流程为：在由不确定参数组成的区间范围内，随机选择1000个样本点，逐次计算所有样本点对应的气动力/热响应值，并将响应值中的最大值和最小值作为区间上界和区间下界。
(15)将本发明方法发生得到的机翼表面最大热流密度qmax、升力系数Cl及阻力系数Cd的区间边界与通过蒙特卡洛方法得到的区间边界进行对比，计算结果如表8-10及图6-8所示。从结果对比中可以看出，本发明方法得到的区间边界与蒙特卡洛方法得到的区间边界吻合较好，并且可以包络蒙特卡洛计算结果，从而验证了本发明方法的有效性及可行性。
表8 机翼表面最大热流密度qmax(Kw/m2)的对比
表9 机翼升力系数Cl的对比
表10 机翼升力系数Cd的对比
综上所述，本发明提出了一种考虑几何不确定性的高超声速机翼气动力/热分析方法。该方法充分考虑外形设计中存在的不确定因素，在不确定参数概率密度未知的情况下，利用区间向量实现不确定参数的定量化表征。在给定飞行工况下，根据伯恩斯坦多项式近似模型构造不确定参数的样本点，利用气动力/热工程算法计算样本点处的气动力/热响应值，同时基于最小二乘策略得到多项式近似模型的待定系数。在给定区间范围内，利用已建立的多项式近似模型求解气动力/热响应值关于各个不确定参数的极大/小值点，从而组合形成最大/小值点的向量，并最终得到气动力/热响应值的区间上界和区间下界，实现考虑几何不确定性的高超声速机翼气动力/热分析。与蒙特卡洛方法相比，本发明方法得到的气动力/热响应区间边界与蒙特卡洛方法得到的区间边界误差较小，并且可以包络蒙特卡洛计算结果，从而验证了本方法的有效性及合理性，为高超声速机翼外形的总体设计提供了新思路。
以上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制，其可扩展应用于高超声速飞行器外形设计领域，凡采用等同转换或者等效替换而形成的技术方案，均落在本发明权利保护范围之内。
Claims (6)
1.一种考虑几何不确定性的高超声速机翼气动力/热分析方法，其特征在于实现步骤如下：
步骤(1)、确定高超声速机翼的平面轮廓参数，包括翼根弦长Cr、翼尖弦长Ct、副翼宽度Ca及机翼半展长L；
步骤(2)、利用类函数/形函数转换方法建立高超声速机翼翼型的参数化表达式，翼型的几何曲线用下列函数表示：
式中，x/c为翼型弦向的无量纲坐标值，y/c为翼型法向的无量纲坐标值，C(x/c)和S(x/c)分别为类型函数和形状函数，zte/c为翼型后缘点的无量纲坐标值，C(x/c)表示为：
对于Clark-Ys翼型，取指数N1＝0.5,N2＝1，S(x/c)表示为：
S(x/c)与翼型前缘半径Rle/c及后缘倾角β满足以下关系：
在形状函数S(x/c)中取n＝3，这样共包含等8个几何设计参数，其中，表示翼型前缘半径，β1和β2表示上、下翼面的后缘倾角，zte/c为翼型后缘点的无量纲坐标值，λ1、λ2、λ1'、λ2'为上、下翼面形状函数多项式的加权系数；
步骤(3)、考虑由加工误差导致翼型几何特征参数存在的波动，将存在不确定性的几何特征参数记为其余参数均视为确定性参数；
步骤(4)、设定翼型几何特征参数α的区间上下界，表示为：
式中，为高超声速机翼翼型几何特征参数的区间下界，为区间上界；
步骤(5)、利用区间[-1,1]内的n阶伯恩斯坦多项式，构造用于计算多项式系数的样本点，将伯恩斯坦多项式模型中子多项式的极值点t记为：
t＝[t0,t1,…,tr,…,tn-1,tn]T (6)
式中：
步骤(6)、将式(6)所表示的极值点映射至区间参数空间αI内的第s个区间参数αs，则在第s维用于计算该维近似多项式系数的样本点表示为：
其中，
式中，m为不确定参数的个数；
步骤(7)、根据式(8)中样本点对应翼型的几何特征参数，建立高超声速机翼的气动外形，进而利用高超声速工程算法计算气动力/热响应值，记为Qs，表示为如下形式：
Qs＝[Qs,0,...,Qs,r,...,Qs,n]T (10)
式中，Qs,r＝Q(αs,r)表示样本点αs,r对应的气动力/热响应值；
步骤(8)、利用最小二乘法计算伯恩斯坦多项式近似模型的待定系数ωs，公式如下：
ωs＝[ωs,0,ωs,1,...,ωs,n]T＝(BTB)-1·BT·Qs (11)
式中，B为伯恩斯坦-范德蒙矩阵，表示为：
其中，bn,r(t)为伯恩斯坦多项式的子多项式，记为：
步骤(9)、根据步骤(8)中得到的系数ωs，建立气动力/热响应值Q关于第s个不确定参数的伯恩斯坦多项式近似模型Ps(t)，表示如下：
步骤(10)、步骤(9)中Ps(t)的最值点通过其导函数的零点与自变量的端点产生，其中导函数零点通过下式计算：
将零点记为ts＝[ts,1,ts,2,...,ts,n]，根据Ps(t)的自变量范围，将零点进一步修改为：
其中Im和Re分别表示对应变量的虚部和实部；
步骤(11)、气动力/热响应值Q关于第s个区间参数的最小值点和最大值点计算为：
式中，Ts＝[ts,-1,1]；
步骤(12)、重复步骤(6)～(11)，以相同的方式计算气动力/热响应值Q关于所有区间参数的最小值点和最大值点，最终可以得到气动力/热响应值在标准区间[-1,1]内的最小值点和最大值点，即：
步骤(13)、将步骤(12)中标准区间内的最小值点和最大值点映射到实际参数空间有：
式中：
步骤(14)、根据αmin和αmax建立其对应的机翼参数化外形，利用高超声速工程算法得到气动力/热响应值Q的区间上下界，表示为：
式中，Q(αmin)和Q(αmax)分别表示几何特征参数αmin和αmax对应的气动力/热响应值；
步骤(15)、将通过上述方法得到的气动力/热响应值Q的区间上下界与蒙特卡洛方法进行对比，验证所建立方法的有效性。
2.根据权利要求1所述的一种考虑几何不确定性的高超声速机翼气动力/热分析方法，其特征在于：所述步骤(1)中，
表1机翼平面轮廓参数
机翼平面轮廓参数由表1确定。
3.根据权利要求1所述的一种考虑几何不确定性的高超声速机翼气动力/热分析方法，其特征在于：所述步骤(2)中，
表2Clark-Ys翼型设计参数
所采用的Clark-Ys标准翼型对应的设计参数由表2确定。
4.根据权利要求1所述的一种考虑几何不确定性的高超声速机翼气动力/热分析方法，其特征在于：所述步骤(4)中，
表3翼型几何特征参数的区间界限
含不确定性的翼型几何特征参数的区间上下界由表3确定。
5.根据权利要求1所述的一种考虑几何不确定性的高超声速机翼气动力/热分析方法，其特征在于：所述步骤(7)中，工程算法的计算流程为：首先根据机翼的几何参数建立气动外形，对几何外形进行非结构网格自由划分，然后分别利用切劈法和达黑姆巴克法计算机翼迎风面和背风面的表面压力系数，在此基础上，利用费雷德经验公式计算机翼驻点加热，从而得到机翼表面的最大热流密度、升力系数、阻力系数等三个气动力/热响应值。
6.根据权利要求1所述的一种考虑几何不确定性的高超声速机翼气动力/热分析方法，其特征在于：所述步骤(15)中，蒙特卡洛方法的计算原理为：在由不确定参数组成的区间范围内，随机选择1000个样本点，逐次计算所有样本点对应的气动力/热响应值，并将响应值中的最大值和最小值作为区间上界和区间下界。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610438739.XA CN106126791B (zh) | 2016-06-17 | 2016-06-17 | 一种考虑几何不确定性的高超声速机翼气动力/热分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610438739.XA CN106126791B (zh) | 2016-06-17 | 2016-06-17 | 一种考虑几何不确定性的高超声速机翼气动力/热分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106126791A CN106126791A (zh) | 2016-11-16 |
| CN106126791B true CN106126791B (zh) | 2018-04-24 |
Family
ID=57470574
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610438739.XA Active CN106126791B (zh) | 2016-06-17 | 2016-06-17 | 一种考虑几何不确定性的高超声速机翼气动力/热分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106126791B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106407586B (zh) * | 2016-09-28 | 2019-08-13 | 北京航空航天大学 | 一种基于不确定优化的飞行器结构功能一体化设计平台系统 |
| CN106872195B (zh) * | 2017-01-06 | 2019-04-09 | 北京临近空间飞行器系统工程研究所 | 一种高速飞行器气动热飞行试验数据的关联分析方法 |
| CN109711048B (zh) * | 2018-12-26 | 2020-09-29 | 北京航空航天大学 | 一种考虑气动和结构多类型设计变量的机翼设计方法 |
| CN109969374B (zh) * | 2019-04-09 | 2021-05-18 | 中国空气动力研究与发展中心计算空气动力研究所 | 用于高超声速边界层转捩研究的标模气动布局及设计方法 |
| US11480073B2 (en) | 2020-11-24 | 2022-10-25 | Rolls-Royce Plc | Gas turbine engine nacelle and method of designing same |
| CN112632869B (zh) * | 2020-12-23 | 2022-09-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于网格框架的非结构附面层网格生成方法 |
| CN115320879B (zh) * | 2022-10-14 | 2022-12-09 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种环量控制翼型后缘柯恩达型面设计方法 |
| CN116610905B (zh) * | 2023-07-20 | 2023-09-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于各向异性尺度修正的反距离权重数据插值方法 |
| CN117313414A (zh) * | 2023-10-23 | 2023-12-29 | 哈尔滨工业大学 | 火星不确定环境下航天器气动参数快速估计方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6584413B1 (en) * | 2001-06-01 | 2003-06-24 | Sandia Corporation | Apparatus and system for multivariate spectral analysis |
| CN105046041A (zh) * | 2015-05-18 | 2015-11-11 | 西北工业大学 | 一种压气机可控扩散叶型稳健优化设计方法 |
| CN105574221A (zh) * | 2014-10-11 | 2016-05-11 | 中国航空工业集团公司西安飞机设计研究所 | 一种改进cst翼型参数化方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8488247B2 (en) * | 2008-10-06 | 2013-07-16 | Purdue Research Foundation | System, method and apparatus for modifying the visibility properties of an object |
| US8374821B2 (en) * | 2010-12-22 | 2013-02-12 | Utility Risk Management Corporation, Llc | Thermal powerline rating and clearance analysis using thermal imaging technology |
-
2016
- 2016-06-17 CN CN201610438739.XA patent/CN106126791B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6584413B1 (en) * | 2001-06-01 | 2003-06-24 | Sandia Corporation | Apparatus and system for multivariate spectral analysis |
| CN105574221A (zh) * | 2014-10-11 | 2016-05-11 | 中国航空工业集团公司西安飞机设计研究所 | 一种改进cst翼型参数化方法 |
| CN105046041A (zh) * | 2015-05-18 | 2015-11-11 | 西北工业大学 | 一种压气机可控扩散叶型稳健优化设计方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106126791A (zh) | 2016-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106126791B (zh) | 一种考虑几何不确定性的高超声速机翼气动力/热分析方法 | |
| CN109711048B (zh) | 一种考虑气动和结构多类型设计变量的机翼设计方法 | |
| Ribeiro et al. | An airfoil optimization technique for wind turbines | |
| CN105718634B (zh) | 一种基于非概率区间分析模型的翼型鲁棒优化设计方法 | |
| CN105183996A (zh) | 面元修正与网格预先自适应计算方法 | |
| Narsipur et al. | CFD analysis of multielement airfoils for wind turbines | |
| CN112632869B (zh) | 一种基于网格框架的非结构附面层网格生成方法 | |
| Churchfield et al. | Modeling wind turbine tower and nacelle effects within an actuator line model | |
| CN106126860A (zh) | 一种考虑加工误差的高超声速机翼鲁棒优化设计方法 | |
| Jeong et al. | Optimization of thick wind turbine airfoils using a genetic algorithm | |
| Barrett et al. | Integrated free‐form method for aerostructural optimization of wind turbine blades | |
| Bourgault-Côté et al. | Multi-layer icing methodologies for conservative ice growth | |
| Yuan et al. | Adaptive-surrogate-based robust optimization of transonic natural laminar flow nacelle | |
| Anderson et al. | Geometry parameterization method for multidisciplinary applications | |
| Jin | Numerical simulation of wind turbine wakes based on actuator line method in NEK5000 | |
| Barrett et al. | Airfoil shape design and optimization using multifidelity analysis and embedded inverse design | |
| Sahai et al. | Effect of fineness ratio on minimum-drag shapes in hypersonic flows | |
| CN105260498B (zh) | 一种大型民机机翼变弯度设计方法 | |
| Yang et al. | Response surface technique for static aeroelastic optimization on a high-aspect-ratio wing | |
| Ahmad et al. | High-order accurate CFD/CSD simulation of the UH60 rotor in forward flight | |
| Takenaka et al. | The Application of MDO Technologies to the Design of a High Performance Small Jet Aircraft-Lessons learned and some practical concerns | |
| Lakshminarayan et al. | Fully Automated Surface Mesh Adaptation in Strand Grid Framework | |
| Li et al. | Sensitivity-Based Parameterization for Aerodynamic Shape Global Optimization | |
| Chen et al. | Aerodynamic optimization for a high-lift airfoil/wing configuration | |
| Tang et al. | Aero-structure coupled optimization for high aspect ratio wings using multi-model fusion method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |