[发明专利]一种高超声速飞行器热气动弹性特性的耦合计算方法有效
申请号: | 201710621122.6 | 申请日: | 2017-07-27 |
公开(公告)号: | CN107368661B | 公开(公告)日: | 2020-06-05 |
发明(设计)人: | 刘磊;杨肖锋;曾磊;代光月;肖光明;魏东;杜雁霞;桂业伟;张昊元 | 申请(专利权)人: | 中国空气动力研究与发展中心计算空气动力研究所 |
主分类号: | G06F30/20 | 分类号: | G06F30/20;G06F119/08 |
代理公司: | 成都九鼎天元知识产权代理有限公司 51214 | 代理人: | 邓世燕 |
地址: | 621052 四*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 高超 声速 飞行器 热气 弹性 特性 耦合 计算方法 | ||
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用于更新流场网格。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国空气动力研究与发展中心计算空气动力研究所,未经中国空气动力研究与发展中心计算空气动力研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201710621122.6/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种预测周期性点阵材料屈服面的分析方法
- 下一篇:一种封闭式安全高效取暖器