CN105631125A - 基于降阶模型的气动-热-结构耦合分析方法 - Google Patents
基于降阶模型的气动-热-结构耦合分析方法 Download PDFInfo
- Publication number
- CN105631125A CN105631125A CN201511001064.4A CN201511001064A CN105631125A CN 105631125 A CN105631125 A CN 105631125A CN 201511001064 A CN201511001064 A CN 201511001064A CN 105631125 A CN105631125 A CN 105631125A
- Authority
- CN
- China
- Prior art keywords
- reduced
- thermal
- aerodynamic
- model
- design
- 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.)
- Pending
Links
Classifications
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种基于降阶模型的气动-热-结构耦合分析方法,包括:建立飞行器分析模型;建立初始气动热降阶模型;验证所建立的气动热降阶模型的总体平均相对误差;建立符合设计精度要求的气动热降阶模型;获得飞行器在全弹道飞行过程中,结构内部温度场的变化情况;在热载荷作用下的飞行器颤振分析。本发明采用降阶模型方法以及加点方法,快速、准确的获得有关物理量降阶模型,提高气动热结构分析效率,能够在飞行器设计过程中提供高效、精确的分析结果,节约了产品设计成本,提高了设计效率,缩短了产品设计周期。
Description
技术领域
本发明属于航空航天领域,具体而言,涉及一种基于降阶模型的气动-热-结构耦合分析方法。
背景技术
高超声速飞行器通常是指能以不小于5马赫速度飞行,以冲压发动机为主要动力,并能实现在临近空间内的长距离飞行的飞行器。高超声速技术是航空宇航技术的一个重要分支。在运输、国防等领域有广阔的应用前景。高超声速飞行器由于涉及到高超声速流动、推进***、飞行器结构以及控制***的复杂的耦合关系,对其进行建模与分析是一项艰巨的任务。而且,随着新型材料以及薄壁结构在现代高超声速飞行器上的广泛应用,高超声速飞行器的气动热弹性问题日益严重。因此准确的对高超声速飞行器的气动热结构问题进行准确的分析在高超声速飞行器设计中显得十分重要。高超声速流动中的气动弹性问题与其他气动弹性问题存在显著的区别,包括严重的气动加热效应、真实气体效应和激波边界层干扰等,这些都给高超声速飞行器的气动热结构分析与设计带来巨大挑战。传统基于简单几何假设以及无粘流假设的各类工程算法无法对高超声速流动给出足够准确预测结果,而基于连续介质力学的计算流体力学虽然可以对高超声速流动中的一系列复杂问题给出比较准确的解答,但由于其计算量巨大,在应用于需要反复修改、迭代的工程设计、优化问题中时,其计算量往往是不可接受的。而基于有限组准确结果来建立的高超声速流动降阶模型可以在计算的准确性与计算效率直接取得一个很好的平衡,适用于实际工程中的反复迭代的设计过程。已有的气动热结构分析方法主要是基于工程算法或是计算流体力学方法。前者虽然计算效率高,但由于无法对高超声速流动中的复杂现象进行准确描述,故基于工程算法的气动热结构分析方法的计算精度受到明显的限制;基于计算流体力学的高超声速气动热结构分析方法,虽然可以给出准确的结果,却受限于其巨大的计算量,在实际工程应用中,面对大量的迭代、优化过程,难以广泛应用。
发明内容
本发明旨在提供一种基于降阶模型的气动-热-结构耦合分析方法。
本发明主要为了解决现有技术中高超声速气动热结构在全弹道计算中计算效率与计算精度相矛盾的问题。本发明通用性强,进而改善高超声速飞行器气动热结构分析方式,提高设计效率,降低设计成本,能够满足高超声速飞行器***多学科优化设计需求。
本发明的目的是通过下述技术方案实现的:
步骤1,根据飞行器设计要求,建立飞行器分析模型,确定设计变量和设计空间B0=[xlb,xub],确定初始样本点数量N0,确定降阶模型相对误差允许阈值Ee;
步骤2,运用试验设计方法(DoE)获得设计空间的样本点X0,进而通过实验或计算的方法获得各样本点的响应值Y0;建立初始气动热降阶模型;
步骤3,采用留一交叉验证法验证上一步所建立的气动热降阶模型的总体平均相对误差E0;
步骤4,建立符合设计精度要求的气动热降阶模型。若气动热降阶模型的总体平均相对误差E0小于相对误差允许阈值Ee,则此气动热降阶模型即为最终符合要求的降阶模型,此时的总体平均相对误差E0即为最终降阶模型总体平均相对误差Efinal。若此气动热降阶模型的总体平均相对误差E0大于相对误差允许阈值Ee,则通过降阶模型加点算法,获得新的样本点Xi add。之后得到新增样本点Xi add处的响应值Yi add。将此时全部的样本点(即原有的X0与所有新增的Xi add)Xall定义为新的初始样本点X’ 0,全部的样本点响应值(即原有的Y0与所有新增的Yi add)Yall定义为新的初始样本点Y’0,返回步骤3;
步骤5,获得飞行器在全弹道飞行过程中,结构内部温度场的变化情况。采用上一步得到的气动热降阶模型,获得飞行器全弹道飞行过程中的气动加热结果,以此为边界条件,进行飞行器结构的瞬态热传导计算。同时,由于考虑到高超声速飞行器在实际飞行时的气动加热随飞行状态变化强烈,又考虑到计算量的权衡,在时间离散上采用自适应时间步长的方法,计算得到全弹道各离散时刻点处的飞行器结构温度场分布;
步骤6,在热载荷作用下的飞行器颤振分析。根据上一步得到的全弹道上一系列离散时刻点处的飞行器结构温度场分布,通过计算各温度场相应的飞行器结构热模态。将飞行器的气动-弹性方程投影到模态域,将各阶模态对应的广义坐标随时间变化关系假设为周期振荡形式,采用气动弹性问题的频域解法得到飞行器在该时刻的颤振特性。
有益效果:
本发明采用降阶模型方法以及加点方法,快速、准确的获得有关物理量降阶模型,提高气动热结构分析效率,能够在飞行器设计过程(例如修改飞行弹道或飞行器结构内部参数)中提供高效、精确的分析结果,节约了产品设计成本;提高了设计、优化效率,缩短了产品设计周期。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为基于降阶模型的气动热结构分析方法;
图2为典型升力面几何外形及其有限元模型;
图3为典型飞行弹道;
图4为典型升力面计算流体网格;
图5为典型时刻升力面结构温度分布;
图6为t=10.8s时刻的升力面前六阶热模态;
图7为典型升力面颤振速度随时间变化规律;
图8为t=10.8s时刻升力面颤振分析中的V-ω图;
图9为t=10.8s时刻升力面颤振分析中的V-g图。
具体实施方式
为了更好的说明本发明的目的与优点,下面通过高超声速典型升力面气动热分析实例,结合附图与表格对本发明做进一步说明。
实施例1
基于降阶模型的气动热结构分析方法,如图1所示,具体步骤如下:
步骤1,以气动热弹性分析中常用的典型升力面F-104机翼为研究对象,机翼结构几何模型及有限元模型图2所示。有限元模型中所有单元均采用壳单元,材料为钛合金。针对典型的高超声速飞行器升力面的气动热结构分析中的气动热降阶模型,变量确定为随飞行弹道变化的飞行马赫数、飞行高度和飞行攻角,设计空间B如表1所示。
表1设计变量和设计空间
初始样本点个数取30个,降阶模型相对误差允许阈值Ee取5%,飞行弹道如图3所示。
步骤2,本实例中气动热分析中的响应值为机翼表面温度,通过CFD商业软件CFD-Fastran计算得到,升力面计算流体网格模型见图4。
在设计空间中采用拉丁超立方抽样方法生成30个初始样本点(即飞行工况)。并通过CFD-Fastran计算得到在各样本点对应工况下的升力面温度分布,即升力面上流体网格节点处的温度值。
步骤3,采用留一交叉验证法得到初始降阶模型的误差。根据N0个初始样本点(第一次执行此步计算时,N0=30),采用Kriging模型作为降阶模型算法,将每个初始样本点作为测试工况,用其他N0-1个样本点及其对应的升力面温度分布来构造Kriging模型,并得到该Kriging模型在所对应的测试工况下升力面的表面温度分布的预测值,并计算该预测值与CFD-Fastran计算得到的升力面温度分布的相对误差。这样可构造ns个Kriging模型(第一次执行此步计算时,ns=30),并且得到ns个预测值与响应值的相对误差,并计算总体平均相对误差E0。
步骤4,若Kriging模型的总体平均相对误差E0小于相对误差允许阈值5%,则此Kriging模型即为符合要求的高超声速气动热降阶模型,即可用于之后的的气动热结构分析中。若Kriging模型的总体平均相对误差E0大于相对误差允许阈值5%,则采用基于模糊聚类的加点算法通过模糊聚类方法获得降阶模型在设计空间中的误差分布较大区域,之后在其中重点增加新的样本点Xi add,之后得到新增样本点Xi add处的响应值Yi add。将此时全部的样本点(即原有的X0与所有新增的Xi add)Xall定义为新的初始样本点X’ 0,全部的样本点响应值(即原有的Y0与所有新增的Yi add)Yall定义为新的初始样本点响应值Y’0,返回步骤3。
步骤5,根据飞行弹道,获得升力面在全弹道飞行过程中,结构内部温度场的变化情况。采用上一步得到的气动热降阶模型,获得在全弹道飞行过程中的升力面表面温度随时间变化关系,通过插值,获得有限元模型表面节点处的温度随时间变化关系。以此为边界条件,采用有限元软件对升力面结构进行瞬态热传导计算。得到全弹道各离散时刻点处的飞行器结构温度场分布,典型时刻升力面结构温度分布如图5所示。
步骤6,在热载荷作用下的飞行器颤振分析。根据上一步得到的全弹道上一系列离散时刻点处的飞行器结构温度场分布,通过有限元方法计算各时刻温度场相应的飞行器结构热模态。不失一般性,选取飞行过程中五个时间点,给出升力面结构的前六阶模态频率,如表2所示。
表2升力面前六阶模态频率变化表
t=10.8s时刻的升力面前六阶热模态如图6所示。根据热模态计算得到各时刻的前六阶模态数据,将飞行器的气动-弹性方程投影到模态域,得到模态坐标下的颤振运动方程,其中,非定常气动力由统一升力面理论得到,分析所需的大气密度、压强、温度等数据根据各时刻弹道飞行高度在标准大气表中插值得到。然后通过频域实现升力面结构的频域颤振计算,得到全弹道飞行过程中各时刻升力面的颤振特性。升力面颤振速度随时间变化规律如图7所示。不失一般性,t=10.8s时刻升力面颤振分析中的V-ω图和V-g图分别如图8和图9所示。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于降阶模型的气动-热-结构耦合分析方法,其特征在于:所述方法包括:
步骤1:建立飞行器分析模型;
步骤2:建立初始气动热降阶模型;
步骤3:验证步骤2所建立的气动热降阶模型的总体平均相对误差;
步骤4:建立符合设计精度要求的气动热降阶模型;
步骤5:获得飞行器在全弹道飞行过程中,结构内部温度场的变化情况;
步骤6:在热载荷作用下的飞行器颤振分析。
2.根据权利要求1所述的方法,其特征在于,所述步骤1包括:
确定设计变量,设计空间,初始样本点数量和降阶模型相对误差允许阈值Ee。
3.根据权利要求2所述的方法,其特征在于,所述步骤2包括:
采用试验设计方法(DoE)获得设计空间的样本点X0,然后通过实验或计算的方法获得各样本点的响应值Y0。
4.根据权利要求3所述的方法,其特征在于,所述步骤3包括:
采用留一交叉验证法验证所建立的初始气动热降阶模型的总体平均相对误差E0。
5.根据权利要求4所述的方法,其特征在于,所述步骤4包括:
比较相对误差允许阈值Ee与相对误差E0,若E0<Ee,则降阶模型符合要求,E0即为最终降阶模型总体平均相对误差Efinal;若E0>Ee,则通过降阶模型加点算法,重新获取初始样本后返回所述步骤3。
6.根据权利要求5所述的方法,其特征在于,所述步骤5包括:
步骤5,采用自适应时间步长的方法计算全弹道各离散时刻点处的飞行器结构温度场分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511001064.4A CN105631125A (zh) | 2015-12-28 | 2015-12-28 | 基于降阶模型的气动-热-结构耦合分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511001064.4A CN105631125A (zh) | 2015-12-28 | 2015-12-28 | 基于降阶模型的气动-热-结构耦合分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105631125A true CN105631125A (zh) | 2016-06-01 |
Family
ID=56046054
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511001064.4A Pending CN105631125A (zh) | 2015-12-28 | 2015-12-28 | 基于降阶模型的气动-热-结构耦合分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631125A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106872195A (zh) * | 2017-01-06 | 2017-06-20 | 北京临近空间飞行器***工程研究所 | 一种高速飞行器气动热飞行试验数据的关联分析方法 |
CN107609307A (zh) * | 2017-10-10 | 2018-01-19 | 北京理工大学 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
CN107766620A (zh) * | 2017-09-27 | 2018-03-06 | 北京理工大学 | 一种基于降阶模型的气动‑热‑结构优化方法 |
CN108170948A (zh) * | 2017-12-27 | 2018-06-15 | 西安电子科技大学 | 高超音速飞行目标流场模型与电磁模型耦合方法 |
CN111767685A (zh) * | 2020-07-08 | 2020-10-13 | 中国飞机强度研究所 | 一种热颤振特性试验方法 |
CN111859534A (zh) * | 2020-06-19 | 2020-10-30 | 空气动力学国家重点实验室 | 热气动弹性分析适用的热固耦合结构动力学降阶模型方法 |
CN111881629A (zh) * | 2020-06-19 | 2020-11-03 | 西北工业大学 | 一种气动热-结构热传导耦合非线性降阶模型方法 |
CN112560167A (zh) * | 2020-11-10 | 2021-03-26 | 北京航空航天大学 | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 |
CN113110100A (zh) * | 2021-03-29 | 2021-07-13 | 北京空天技术研究所 | 适用于飞行器的模型降阶分析方法 |
CN114818550A (zh) * | 2022-06-30 | 2022-07-29 | 中国飞机强度研究所 | 一种飞机振动试验中时变气动载荷地面等效模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN103902782A (zh) * | 2014-04-11 | 2014-07-02 | 北京理工大学 | 基于pod和代理模型的高超声速气动热模型降阶方法 |
CN105138766A (zh) * | 2015-08-21 | 2015-12-09 | 北京理工大学 | 基于模糊聚类的高超声速气动热降阶模型的加点方法 |
-
2015
- 2015-12-28 CN CN201511001064.4A patent/CN105631125A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN103902782A (zh) * | 2014-04-11 | 2014-07-02 | 北京理工大学 | 基于pod和代理模型的高超声速气动热模型降阶方法 |
CN105138766A (zh) * | 2015-08-21 | 2015-12-09 | 北京理工大学 | 基于模糊聚类的高超声速气动热降阶模型的加点方法 |
Non-Patent Citations (2)
Title |
---|
李昱霖 等: "高超声速飞行器气动热结构耦合问题研究", 《飞航导弹》 * |
陈鑫 等: "基于本征正交分解和代理模型的高超声速气动热模型降阶研究", 《航空学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106872195A (zh) * | 2017-01-06 | 2017-06-20 | 北京临近空间飞行器***工程研究所 | 一种高速飞行器气动热飞行试验数据的关联分析方法 |
CN106872195B (zh) * | 2017-01-06 | 2019-04-09 | 北京临近空间飞行器***工程研究所 | 一种高速飞行器气动热飞行试验数据的关联分析方法 |
CN107766620A (zh) * | 2017-09-27 | 2018-03-06 | 北京理工大学 | 一种基于降阶模型的气动‑热‑结构优化方法 |
CN107609307A (zh) * | 2017-10-10 | 2018-01-19 | 北京理工大学 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
CN107609307B (zh) * | 2017-10-10 | 2019-09-10 | 北京理工大学 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
CN108170948A (zh) * | 2017-12-27 | 2018-06-15 | 西安电子科技大学 | 高超音速飞行目标流场模型与电磁模型耦合方法 |
CN108170948B (zh) * | 2017-12-27 | 2019-06-14 | 西安电子科技大学 | 高超音速飞行目标流场模型与电磁模型耦合方法 |
CN111859534A (zh) * | 2020-06-19 | 2020-10-30 | 空气动力学国家重点实验室 | 热气动弹性分析适用的热固耦合结构动力学降阶模型方法 |
CN111881629A (zh) * | 2020-06-19 | 2020-11-03 | 西北工业大学 | 一种气动热-结构热传导耦合非线性降阶模型方法 |
CN111881629B (zh) * | 2020-06-19 | 2022-11-08 | 西北工业大学 | 一种气动热-结构热传导耦合非线性降阶模型方法 |
CN111859534B (zh) * | 2020-06-19 | 2023-12-26 | 空气动力学国家重点实验室 | 热气动弹性分析适用的热固耦合结构动力学降阶模型方法 |
CN111767685A (zh) * | 2020-07-08 | 2020-10-13 | 中国飞机强度研究所 | 一种热颤振特性试验方法 |
CN111767685B (zh) * | 2020-07-08 | 2022-05-06 | 中国飞机强度研究所 | 一种热颤振特性试验方法 |
CN112560167A (zh) * | 2020-11-10 | 2021-03-26 | 北京航空航天大学 | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 |
CN112560167B (zh) * | 2020-11-10 | 2023-03-14 | 北京航空航天大学 | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 |
CN113110100A (zh) * | 2021-03-29 | 2021-07-13 | 北京空天技术研究所 | 适用于飞行器的模型降阶分析方法 |
CN114818550A (zh) * | 2022-06-30 | 2022-07-29 | 中国飞机强度研究所 | 一种飞机振动试验中时变气动载荷地面等效模拟方法 |
CN114818550B (zh) * | 2022-06-30 | 2022-09-13 | 中国飞机强度研究所 | 一种飞机振动试验中时变气动载荷地面等效模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105631125A (zh) | 基于降阶模型的气动-热-结构耦合分析方法 | |
CN111079228B (zh) | 一种基于流场预测的气动外形优化方法 | |
CN105183996A (zh) | 面元修正与网格预先自适应计算方法 | |
CN105718634A (zh) | 一种基于非概率区间分析模型的翼型鲁棒优化设计方法 | |
Chang et al. | A parallel implicit hole-cutting method based on background mesh for unstructured Chimera grid | |
CN105138766B (zh) | 基于模糊聚类的高超声速气动热降阶模型的加点方法 | |
CN107766620A (zh) | 一种基于降阶模型的气动‑热‑结构优化方法 | |
Deng et al. | Overall performance analysis–oriented aerodynamic configuration optimization design for hypersonic vehicles | |
Wang et al. | An adaptive sampling strategy for construction of surrogate aerodynamic model | |
Yang et al. | Static aeroelastic modeling and rapid analysis of wings in transonic flow | |
Cummings et al. | Supersonic, turbulent flow computation and drag optimization for axisymmetric afterbodies | |
Bisson et al. | Adjoint-based aerodynamic optimization framework | |
Zheng et al. | Flutter analyses of complete aircraft based on hybrid grids and parallel computing | |
Duong et al. | Fluid–structure coupling in time domain for dynamic stall using purely lagrangian vortex method | |
Fang et al. | Design optimization of unsteady airfoils with continuous adjoint method | |
Lyu et al. | Strategies for solving high-fidelity aerodynamic shape optimization problems | |
Takenaka et al. | The Application of MDO Technologies to the Design of a High Performance Small Jet Aircraft-Lessons learned and some practical concerns | |
Wang et al. | Fundamentals of Computational Fluid Dynamics | |
Nejati et al. | Application of the adjoint optimisation of shock control bump for ONERA-M6 wing | |
Vedam et al. | Evaluation of Gradient and Curvature-Based Adaptive Mesh Refinement for Viscous Transonic Flows | |
Zuo et al. | Efficient aero-structural design optimization: Coupling based on reverse iteration of structural model | |
Peng et al. | Rapid Aerodynamic Shape Optimization With Payload Size Constraints for Hypersonic Vehicle | |
Azab et al. | High order aerodynamic optimization using new hybrid sequential quadratic programing-particle swarm intelligence technique | |
Hamdi et al. | Link-Wise Artificial Compressibility Method for attached and separated flows past airfoils | |
Katzenmeier et al. | Correction technique for quality improvement of doublet lattice unsteady loads by introducing CFD small disturbance aerodynamics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160601 |