[发明专利]考虑气动面曲面效应和法向运动的非定常气动力计算方法有效
申请号: | 201910206454.7 | 申请日: | 2019-03-19 |
公开(公告)号: | CN110162823B | 公开(公告)日: | 2020-12-08 |
发明(设计)人: | 刘燚;刘凯;操小龙;何海波 | 申请(专利权)人: | 北京机电工程研究所 |
主分类号: | G06F30/20 | 分类号: | G06F30/20;G06F30/15;G06F119/14 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 100074 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 考虑 气动 曲面 效应 运动 非定常 气动力 计算方法 | ||
1.一种考虑气动面曲面效应和法向运动的非定常气动力计算方法,其特征在于,所述非定常气动力计算方法包括:
对飞行器的变形后气动面进行网格划分;
对飞行器进行结构动力学分析以获取飞行器的结构弹性模态;
从所述飞行器的结构弹性模态中提取与变形后气动面相关的结构弹性模态并插值到变形后气动面的气动网格上;
在变形后气动面的网格上布置偶极子基本解,求解考虑气动面变形的核函数,根据所述核函数求解气动力影响系数矩阵;考虑气动面变形的核函数为λ1=xi-xj,其中,(xi,yi,zi)是接收点坐标,(xj,yj,zj)是扰动点坐标,ni是(xi,yi,zi)处的翼面法线方向,nj是(xj,yj,zj)处的翼面法线方向,ω为结构振动圆频率,U∞为远方来流速度,R'为接收点与扰动点之间的距离,a∞为远前方来流的声速,M∞为来流马赫数;
计算变形后气动面的气动网格的当地法向量;
根据变形后气动面的结构弹性模态和当地法向量求解变形后气动面的法向模态,所述变形后气动面的法向模态根据变形后气动面的结构弹性模态与所述当地法向量相乘来获取;根据变形后气动面的法向模态求解法向运动边界条件;
根据所述法向运动边界条件和所述气动力影响系数矩阵求解变形后气动面网格上的非定常气动力;
根据变形后气动面网格上的非定常气动力和变形后气动面的法向模态以获取广义化的曲面非定常气动力,所述广义化的曲面非定常气动力根据变形后气动面网格上的非定常气动力与所述变形后气动面的法向模态相乘来获取。
2.根据权利要求1所述的考虑气动面曲面效应和法向运动的非定常气动力计算方法,其特征在于,所述法向运动边界条件为其中,w为网格无量纲法向下洗速度,k为减缩频率,b为参考弦长,f为所述变形后气动面的法向模态,n'为曲面法向,为实部,为虚部。
3.根据权利要求2所述的考虑气动面曲面效应和法向运动的非定常气动力计算方法,其特征在于,变形后气动面网格上的非定常气动力Δcp可根据w=DΔcp来获取,其中,w为网格无量纲法向下洗速度,D为气动力影响系数矩阵,Δcp为气动网格的非定常气动力组成的列向量。
4.根据权利要求3所述的考虑气动面曲面效应和法向运动的非定常气动力计算方法,其特征在于,所述气动力影响系数矩阵D可根据来求解,其中,Δxj为第j个网格的中剖面长度,lj为第j个网格的过1/4弦点的展长,χj为网格1/4弦线后掠角,Kij为核函数。
5.根据权利要求1至4中任一项所述的考虑气动面曲面效应和法向运动的非定常气动力计算方法,其特征在于,所述变形后气动面的弦向网格至少大于5,展向网格可根据网格单元长细比进行确定。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京机电工程研究所,未经北京机电工程研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910206454.7/1.html,转载请声明来源钻瓜专利网。