[发明专利]弹性体飞行器飞行仿真方法、系统及计算机存储介质有效
申请号: | 201910970180.9 | 申请日: | 2019-10-12 |
公开(公告)号: | CN110874501B | 公开(公告)日: | 2021-08-31 |
发明(设计)人: | 杨磊;张陈安;刘文;李文皓 | 申请(专利权)人: | 中国科学院力学研究所 |
主分类号: | G06F30/15 | 分类号: | G06F30/15;G06F30/28;G06F119/14;G06F113/08 |
代理公司: | 北京维正专利代理有限公司 11508 | 代理人: | 李传亮 |
地址: | 100190 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 弹性体 飞行器 飞行 仿真 方法 系统 计算机 存储 介质 | ||
1.一种弹性体飞行器飞行仿真方法,在仿真过程中,
在前m步中进行CFD/CSD/RBD耦合计算,时间步长采用ts′;其中,N0≤m<N,N为仿真的总步数,N0等于时间推进阶数;
当m<k≤N时,进行多时间尺度的UAM/CFD/CSD/RBD耦合计算,判断k是否为n的整数倍,如果是,则进行CFD/CSD/RBD耦合计算,其中,RBD时间步长采用n·ts,CSD时间步长采用ts;如果否,则进行UAM/CSD耦合计算,时间步长采用ts;其中,k为当前步数,10≤n≤N,ts≤ts′;
完成CFD/CSD/RBD耦合计算或UAM/CSD耦合计算后输出动力学响应参数,进入第k+1步的仿真,直至完成N步的仿真,输出网格、流场及续算数据;
其中,UAM为气动力模型,n为CFD/CSD/RBD耦合计算间隔数。
2.根据权利要求1所述的飞行仿真方法,所述方法还包括,在进行CFD/CSD/RBD耦合计算之前,先根据第k-1步中得到的网格节点坐标对计算网格进行动态变形处理,在进行CFD/CSD/RBD耦合计算之后,存储变形后的网格节点坐标,以为后续计算网格处理做准备。
3.根据权利要求1所述的飞行仿真方法,所述方法还包括,在进行UAM/CSD耦合计算之后,判断k+1是否为n的整数倍,如果是,则存储变形后的网格节点坐标,以为k+1步的计算网格处理做准备。
4.根据权利要求1到3之一所述的飞行仿真方法,所述方法还包括,判断当前步数k是否为Nf的整数倍,如果是,则输出当前的网格、流场及续算数据;其中,Nf为设置的流场输出保存的间隔步数。
5.根据权利要求4所述的飞行仿真方法,所述方法还包括,在进行CFD/CSD/RBD耦合计算或UAM/CSD耦合计算之前,或在进行CFD/CSD/RBD耦合计算或UAM/CSD耦合计算之后判断当前步数k是否为Nf的整数倍。
6.根据权利要求4所述的飞行仿真方法,具体仿真方法步骤包括:
S1,读取主程序设置参数、网格、结构模台参数、飞行力学参数和气动力模型;其中,主程序设置参数中,设置初始时间步长ts′、基准时间步长ts、CFD/CSD/RBD耦合计算间隔n、仿真总步数N、初始求解步数m和流场输出保存的间隔步数Nf的值;
S2,判断当前步数k是否小于等于N,如果是,则进入S3;如果否,则输出网格、流场及续算数据,仿真结束;
S3,判断当前步数k是否小于m或k为n的整数倍,如果是,则进入S4,如果否,则进入S6;
S4,根据第k-1步中得到的网格节点坐标对计算网格进行动态变形处理,进入S5;
S5,进行CFD/CSD/RBD耦合计算并进入S7;
S6,进行UAM/CSD耦合计算后,判断k+1是否为n的整数倍,如果是,则进入S7;如果否,则进入S8;
S7,存储变形后的网格节点坐标并进入S8;
S8,输出动力学响应参数,进入S9;
S9,判断当前步数k是否为Nf的整数倍,如果是,则输出当前的网格、流场及续算数据,进入S10;如果否,则直接进入S10;
S10,对k加1,进入S2。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院力学研究所,未经中国科学院力学研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910970180.9/1.html,转载请声明来源钻瓜专利网。
- 上一篇:微信平台上的电子优惠券系统
- 下一篇:一种生产玉米面粉的方法