CN107368661B - 一种高超声速飞行器热气动弹性特性的耦合计算方法 - Google Patents
一种高超声速飞行器热气动弹性特性的耦合计算方法 Download PDFInfo
- Publication number
- CN107368661B CN107368661B CN201710621122.6A CN201710621122A CN107368661B CN 107368661 B CN107368661 B CN 107368661B CN 201710621122 A CN201710621122 A CN 201710621122A CN 107368661 B CN107368661 B CN 107368661B
- Authority
- CN
- China
- Prior art keywords
- coupling
- flow field
- aircraft
- temperature
- field
- 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
Images
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
-
- 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)
- 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)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种高超声速飞行器热气动弹性特性的耦合计算方法,从气动力、气动热、结构传热和结构应力/变形等物理场特征时间出发,在考虑现有计算资源和不降低耦合计算精度的前提下,有效减小了耦合分析方法的计算量,可用于高超声速飞行器实际结构的热气动弹性分析。本发明通过监控结构温度场的变化情况实现了耦合时间步长的动态调整,在有效保证耦合计算精度的情况下,大幅度提升耦合计算效率这一难题。该方法可有效实现高超声速飞行器整机结构或部件的热气动弹性特性分析;同时,对同样涉及飞行器流‑热‑固耦合计算问题也具备求解能力,譬如气动热与传热耦合问题、结构热安全性评估问题等。
Description
技术领域
本发明涉及一种高超声速飞行器热气动弹性特性的耦合计算方法。
背景技术
新一代高超声速飞行器受高速度、高机动性及长航时的任务需求影响,多采用细长体、升力体或乘波体构型。同时,由于质量限制具有较大的结构柔度。在长时间的气动力、气动热的综合作用下,温升引起的结构应力、刚度和模态的显著变化,对飞行器气动弹性特性带来了较大影响,由此形成的热气动弹性问题给飞行器热安全带来了前所未有的挑战。针对真实飞行器结构的热气动弹性特性预测问题,国内外均无试验研究能力,主要依靠数值分析开展研究。且受制于气动力/热/结构多场耦合的复杂性,目前尚无开展真实飞行器结构热气动弹性计算分析的能力。
现有流-热-固耦合计算方法通常采用紧耦合方案,耦合时间根据各物理场数值求解方法中最小的一个决定。此时,开展热气动弹性问题计算需要花费相当长的时间(计算时间通常以月为单位计算)。
发明内容
为了克服现有技术的上述缺点,本发明提出了一种高超声速飞行器热气动弹性特性的耦合计算方法,通过合理的调整耦合时间间隔,大幅度降低总计算量,计算时间可以以天为单位计算,使得热气动弹性问题计算具有实用价值。
本发明解决其技术问题所采用的技术方案是:一种高超声速飞行器热气动弹性特性的耦合计算方法,包括如下步骤:
步骤一、对流场和温度场进行初始化;
步骤二、对ti时刻绕飞行器的流体域进行稳态流场计算,获得稳态流场近飞行器壁面的热流密度qw,i和壁面压力pw,i;
步骤三、将稳态流场的壁面压力pw,i和壁面热流密度qw,i插值传递给飞行器结构表面,分别作为结构热应力/热变形计算的力边界载荷条件和结构热传导计算的热边界载荷条件;
步骤四、基于施加的热边界载荷,自ti时刻开始在飞行器结构内进行瞬态热传导计算,并监控全飞行器结构的温升变化情况,当结构温升变化达到预设的最大温升阈值△Tmax时,停止热传导计算并将当前时刻记录为ti+1,同时保存ti+1时刻的结构温度场和结构壁面温度Tw,i+1;
步骤五、计算ti+1时刻的结构位移场和动力响应特性;
步骤六、将ti+1时刻的结构壁面温度Tw,i+1和壁面位移量Uw,i+1传递给流场;
步骤七、返回步骤二进入下一个时间步长继续进行计算,直至耦合推进时间到达预设的结束时间tend。
与现有技术相比,本发明的积极效果是:本发明从气动力、气动热、结构传热和结构应力/变形等物理场特征时间出发,在考虑现有计算资源的基础上、在不降低耦合计算精度的前提下,有效减小了耦合分析方法的计算量,可用于高超声速飞行器实际结构的热气动弹性分析。本发明通过监控结构温度场的变化情况实现了耦合时间步长的动态调整,在有效保证耦合计算精度的情况下,大幅度提升耦合计算效率这一难题。该方法可有效实现高超声速飞行器整机结构或部件的热气动弹性特性分析;同时,对同样涉及飞行器流-热-固耦合计算问题也具备求解能力,譬如气动热与传热耦合问题、结构热安全性评估问题等。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1为本发明的流程图。
具体实施方式
一种高超声速飞行器热气动弹性特性的耦合计算方法,如图1所示,包括如下步骤:
步骤一:在初始t0时刻,通过给定飞行器结构温度Tinitial初始化结构温度场,并给定结构位移约束条件;通过物理场间数据传递方法将结构壁面温度传递给流场作为流场计算的边界条件,并根据飞行状态初始化流场。
步骤二:基于温度边界条件和流动状态,对绕飞行器的流体域进行稳态流场计算,获得稳态流场近飞行器壁面的热流密度qw和壁面压力pw。
步骤三:通过场间数据传递方法将稳态流场的壁面压力pw和壁面热流密度qw插值传递给飞行器结构表面,分别作为结构热应力或热变形计算的力边界载荷条件和结构热传导计算的热边界载荷条件。
步骤四:基于施加的热边界载荷,在飞行器结构内进行瞬态热传导计算,并监控全飞行器结构的温升变化情况。当结构温升变化达到预设的最大温升阈值△Tmax时,停止热传导计算并确认此时所处时刻t1,保存t1时刻的结构温度场和结构壁面温度Tw。在热气动弹性问题中,该阈值△Tmax取10K可有效平衡耦合计算精度与计算效率。
步骤五:时间推进从t0到t1,基于施加的力边界载荷和飞行器结构t1时刻的瞬态温度场,在结构内进行结构热应力/热变形计算和结构动力响应计算,获得t1时刻的结构位移场和动力响应特性;
步骤六:通过场间数据传递方法将t1时刻的结构壁面温度Tw和壁面位移量Uw传递给流场,其中壁面温度Tw作为流场计算的边界条件,而壁面位移量Uw则用于更新流场网格;
步骤七:至此完成一个流-热-固耦合计算时间步长的计算,然后进入下一个耦合计算时间步长继续进行计算,直至耦合推进时间到达预设的结束时间tend。
Claims (6)
1.一种高超声速飞行器热气动弹性特性的耦合计算方法,其特征在于:包括如下步骤:
步骤一、对流场和温度场进行初始化;
步骤二、对ti时刻绕飞行器的流体域进行稳态流场计算,获得稳态流场近飞行器壁面的热流密度qw,i和壁面压力pw,i;
步骤三、将稳态流场的壁面压力pw,i和壁面热流密度qw,i插值传递给飞行器结构表面,分别作为结构热应力/热变形计算的力边界载荷条件和结构热传导计算的热边界载荷条件;
步骤四、基于施加的热边界载荷,自ti时刻开始在飞行器结构内进行瞬态热传导计算,并监控全飞行器结构的温升变化情况,当结构温升变化达到预设的最大温升阈值△Tmax时,停止热传导计算并将当前时刻记录为ti+1,同时保存ti+1时刻的结构温度场和结构壁面温度Tw,i+1;
步骤五、计算ti+1时刻的结构位移场和动力响应特性;
步骤六、将ti+1时刻的结构壁面温度Tw,i+1和壁面位移量Uw,i+1传递给流场;
步骤七、返回步骤二进入下一个时间步长继续进行计算,直至耦合推进时间到达预设的结束时间tend。
2.根据权利要求1所述的一种高超声速飞行器热气动弹性特性的耦合计算方法,其特征在于:步骤一所述对流场进行初始化是指通过物理场间数据传递方法将结构壁面温度传递给流场作为流场计算的边界条件,并根据飞行状态对流场进行初始化。
3.根据权利要求1所述的一种高超声速飞行器热气动弹性特性的耦合计算方法,其特征在于:步骤一所述对温度场进行初始化是指通过给定飞行器结构温度Tinitial对结构温度场进行初始化,同时给定结构位移约束条件。
4.根据权利要求1所述的一种高超声速飞行器热气动弹性特性的耦合计算方法,其特征在于:步骤三所述将稳态流场的壁面压力pw,i和壁面热流密度qw,i插值传递给飞行器结构表面是通过场间数据传递方法实现的。
5.根据权利要求1所述的一种高超声速飞行器热气动弹性特性的耦合计算方法,其特征在于:所述△Tmax取10K。
6.根据权利要求1所述的一种高超声速飞行器热气动弹性特性的耦合计算方法,其特征在于:步骤六所述将ti+1时刻的结构壁面温度Tw,i+1和壁面位移量Uw,i+1传递给流场是通过场间数据传递方法实现的,其中壁面温度Tw,i+1作为流场计算的边界条件,壁面位移量Uw,i+1用于更新流场网格。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710621122.6A CN107368661B (zh) | 2017-07-27 | 2017-07-27 | 一种高超声速飞行器热气动弹性特性的耦合计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710621122.6A CN107368661B (zh) | 2017-07-27 | 2017-07-27 | 一种高超声速飞行器热气动弹性特性的耦合计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107368661A CN107368661A (zh) | 2017-11-21 |
CN107368661B true CN107368661B (zh) | 2020-06-05 |
Family
ID=60307275
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710621122.6A Active CN107368661B (zh) | 2017-07-27 | 2017-07-27 | 一种高超声速飞行器热气动弹性特性的耦合计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107368661B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109918808B (zh) * | 2019-03-13 | 2023-10-20 | 北京宇航系统工程研究所 | 一种气热弹三场耦合仿真分析方法 |
CN111859531B (zh) * | 2020-06-15 | 2023-10-13 | 空气动力学国家重点实验室 | 一种求解化学非平衡流动化学反应源项的预估修正方法 |
CN112329148A (zh) * | 2020-11-18 | 2021-02-05 | 西北工业大学 | 一种高超声速飞行器热防护结构优化方法及系统 |
CN112651156B (zh) * | 2020-12-21 | 2022-03-04 | 中国科学院力学研究所 | 一种高速来流条件下激光毁伤效应的多场耦合分析方法 |
CN112989485B (zh) * | 2021-01-18 | 2022-04-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于热壁修正的沿弹道热流插值方法 |
CN114626313B (zh) * | 2022-03-04 | 2023-05-16 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种可解析时变热响应的高速气动热cfd求解方法 |
CN114429052B (zh) * | 2022-04-06 | 2022-06-21 | 中国空气动力研究与发展中心计算空气动力研究所 | 复合材料广义热弹性耦合问题计算方法 |
CN116720264B (zh) * | 2023-08-04 | 2023-10-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 考虑气动力/热累积变形反向几何预置的气动布局方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN105095603A (zh) * | 2015-09-09 | 2015-11-25 | 哈尔滨工业大学 | 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法 |
CN106844914A (zh) * | 2017-01-09 | 2017-06-13 | 西北工业大学 | 一种空天飞行器机翼振动响应的快速仿真方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0840397A (ja) * | 1994-07-28 | 1996-02-13 | Mitsubishi Heavy Ind Ltd | 高速飛行体 |
CN102012953B (zh) * | 2010-11-04 | 2013-05-08 | 西北工业大学 | Cfd/csd耦合求解非线性气动弹性仿真方法 |
WO2017106048A1 (en) * | 2015-12-16 | 2017-06-22 | University Of Notre Dame Du Lac | Methods and apparatus for pulsed-dc dielectric barrier discharge plasma actuator and circuit |
-
2017
- 2017-07-27 CN CN201710621122.6A patent/CN107368661B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077259A (zh) * | 2011-10-26 | 2013-05-01 | 上海机电工程研究所 | 高超声速导弹多场耦合动力学一体化仿真分析方法 |
CN105095603A (zh) * | 2015-09-09 | 2015-11-25 | 哈尔滨工业大学 | 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法 |
CN106844914A (zh) * | 2017-01-09 | 2017-06-13 | 西北工业大学 | 一种空天飞行器机翼振动响应的快速仿真方法 |
Non-Patent Citations (2)
Title |
---|
高超声速翼面的流固热多场耦合分析;徐飞;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20160715;第2016年卷(第7期);正文第3.2.1节 * |
高超声速进气道前缘流场-热-结构耦合分析;张胜涛 等;《空气动力学学报》;20170630;第35卷(第3期);正文第1.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN107368661A (zh) | 2017-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107368661B (zh) | 一种高超声速飞行器热气动弹性特性的耦合计算方法 | |
Wolfe et al. | CFD calculations of S809 aerodynamic characteristics | |
WO2019071763A1 (zh) | 一种考虑管道动态特性的电-热耦合系统状态估计方法 | |
CN106919731B (zh) | 一种用于不同风向角的风电机组尾流确定方法 | |
CN109058562A (zh) | 一种高速开关阀复合pwm控制方法 | |
CN106980734B (zh) | 一种考虑材料和接触非线性的rpv密封性能计算方法 | |
CN110705066B (zh) | 基于投影积分的气电耦合园区综合能源系统动态仿真方法 | |
CN115422654B (zh) | 基于cfd/csd技术的跨/超声速飞行器高效高精度颤振时域分析方法 | |
CN107871058A (zh) | 电热联合系统的潮流计算方法、装置、设备及存储介质 | |
CN101240775A (zh) | 风电机组液压变桨距系统的控制方法 | |
Tabatabaei et al. | Wind Turbine Aerodynamic Modeling in Icing Condition: Three-Dimensional RANS-CFD Versus Blade Element Momentum Method | |
CN110737197B (zh) | 一种基于无模型的柴油机转速自适应自抗扰控制方法 | |
CN103941724A (zh) | 一种低故障诊断精度要求的长时延网络控制系统容错控制方法 | |
CN102033975B (zh) | 一种换流阀用饱和电抗器热力性能设计方法 | |
CN117311135A (zh) | 一种双余度eha驱动桨扇发动机变距系统的建模方法 | |
CN104749959A (zh) | 一种基于广义滑模估计器的机组变桨容错控制方法 | |
CN112084680A (zh) | 一种基于dqn算法的能源互联网优化策略方法 | |
CN109635370B (zh) | 开裂式阻力方向舵静气动弹性特性分析方法 | |
CN111031613B (zh) | 一种复杂对流环境下的功率器件主动热控制方法 | |
Dai et al. | Modeling and analysis of electrical heating system based on entransy dissipation-based thermal resistance theory | |
CN201826286U (zh) | 闸门槽热管融冰装置 | |
CN113459083B (zh) | 一种事件触发下的机械臂自适应固定时间控制方法及系统 | |
CN203798530U (zh) | 冰风洞风速管 | |
CN110705135A (zh) | 一种用于混合三通管道的高周疲劳寿命估算方法 | |
CN203270657U (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 | ||
GR01 | Patent grant |