[发明专利]一种固定翼飞机阵风响应建模分析方法有效
申请号: | 201811551488.1 | 申请日: | 2018-12-18 |
公开(公告)号: | CN109840349B | 公开(公告)日: | 2023-02-10 |
发明(设计)人: | 龚喜盈;王丹;李道春;田茂江;向锦武;张琳;王安文;李凡凡 | 申请(专利权)人: | 西安爱生技术集团公司;北京航空航天大学 |
主分类号: | G06F30/15 | 分类号: | G06F30/15;G06F30/23;G06F17/16;G06F30/17 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 710000 陕西省西安市高*** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 固定 飞机 阵风 响应 建模 分析 方法 | ||
1.一种固定翼飞机阵风响应建模分析方法,具体步骤如下:
1)针对带被动阵风减缓装置的固定翼飞机,进行有限元分析,提取飞机结构的前几阶模态,这些模态代表飞机结构的动力学响应,有效降低结构模型阶数;
2)分析翼尖装置的结构模态,使用正交分解的方法将翼尖装置的结构模态和机翼的结构模态解耦,得到翼尖装置的转动模态;
3)定义单位阶跃输入为:
用多项式分别拟合各阶模态的变形,将变形量插值到气动网格节点上,通过动网格方法,使用CFD计算所述结构各模态单位阶跃输入下的各时间步的非定常广义气动力,记为y(n);
4)用一阶的Volterra级数表示非定常广义气动力其中h(k),k=1,2,…,K是Volterra级数的核函数,由下式得到:
5)构造Hankel矩阵:
其中,n为大于0的整数,r和s用来确定矩阵的大小,r表示该矩阵的行数,s表示列数,该矩阵为对称矩阵,r与s满足n+s+r-2≤K;
6)对上述Hankel矩阵进行奇异值分解得到构造出广义非定常气动力的状态方程:
其中,xa是该状态空间的状态变量,u为该状态空间的输入,该输入为结构的形变量,即广义位移,Fa为系统的输出量,即为广义气动力;
其中,与定义如下:
Ip与Im是维度为p与m的单位方阵,0p与0m是维度为p与m的零方阵,与的维度分别是p×rp,m×sm;
7)引入1-cos阵风模型,飞机飞行动力学模型,结构动力学模型,并与步骤6)中的广义非定常气动力模型耦合,得到带有被动阵风减缓装置的固定翼飞机阵风响应模型,使用该阵风响应模型进行阵风响应分析,计算各阶结构模态在阵风作用下的响应。
2.根据权利要求1所述的建模分析方法,所述步骤7)中的飞机的结构动力学模型为:
其中,M1为飞机的广义质量,C1为飞机的广义阻尼,K1为飞机的广义刚度,M2为被动翼尖广义质量,C2为飞机的广义阻尼,K2为飞机的广义刚度,A1为舵面的气动力矩阵,A12为飞机运动作用在被动翼尖气动力矩阵,A2为被动翼尖运动的气动力矩阵,xd为结构广义坐标,为结构广义坐标对时间求导,为结构广义坐标对时间的二阶导。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安爱生技术集团公司;北京航空航天大学,未经西安爱生技术集团公司;北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201811551488.1/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种三重加速的拓扑优化方法
- 下一篇:一种电力系统动态仿真方法及系统