[发明专利]一种大包线强机动飞行器动力学高精度仿真方法有效
申请号: | 201610847855.7 | 申请日: | 2016-09-23 |
公开(公告)号: | CN106484969B | 公开(公告)日: | 2019-08-09 |
发明(设计)人: | 李争学;刘刚;韩鹏鑫;李杰奇;晏资湘;郭金花;王宁宇;刘峰;张广春;王飞;蔡巧言;朱永贵;韩威;曾凡文 | 申请(专利权)人: | 中国运载火箭技术研究院 |
主分类号: | G06F17/50 | 分类号: | G06F17/50 |
代理公司: | 中国航天科技专利中心 11009 | 代理人: | 臧春喜 |
地址: | 100076 北京*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 一种大包线强机动飞行器动力学高精度仿真方法,属于飞行器动力学与控制领域。该方法首先建立了大包线强机动飞行器的刚体弹性耦合动力学高阶模型,建模过程仅考虑小幅线性振动假设,充分考虑了飞行器大包线强机动飞行过程中显著且快时变的气动力和力矩、发动机推力和力矩、重力、姿态强机动、气动舵快速运动与结构弹性振动之间的相互耦合影响,模型包含了全面的高阶非线性项,刚体运动与弹性振动之间的耦合影响项、姿态机动和气动舵快速运动对弹性振动和姿态的影响项。因此,该高阶模型能够真实反映大包线强机动飞行器真实状态,利用该高阶模型进行仿真分析,可以用于验证飞行器设计的合理性,验证结果较传统模型更准确可靠。 | ||
搜索关键词: | 一种 大包 机动 飞行器 动力学 高精度 仿真 方法 | ||
【主权项】:
1.一种大包线强机动飞行器动力学高精度仿真方法,其特征在于包括如下步骤:(1)建立大包线强机动飞行器的刚体弹性耦合动力学高阶模型,所述刚体弹性耦合动力学高阶模型包括地面坐标系下的飞行器平动动力学方程:
飞行器本体坐标系下的飞行器转动动力学方程:
飞行器本体坐标系下的飞行器弹性振动动力学方程:
Abg为地面坐标系到本体坐标系的坐标转换矩阵;Abi为惯性坐标系到本体坐标系的坐标转换矩阵;Abrk为气动舵k本体坐标系到飞行器本体坐标系的坐标转换矩阵;Abv为速度坐标系到本体坐标系的坐标转换矩阵;Aei为惯性坐标系到地球坐标系的坐标转换矩阵;Agb为本体坐标系到地面坐标系的坐标转换矩阵;Age为地球坐标系到地面坐标系的坐标转换矩阵;Agrk为气动舵k本体坐标系到地面坐标系的坐标转换矩阵;Aog为地面坐标系到弹道坐标系的坐标转换矩阵;Arkb为飞行器本体坐标系到气动舵k本体坐标系的坐标转换矩阵;Avo为弹道坐标系到速度坐标系的坐标转换矩阵;dbg为b系原点相对g系原点的矢径在b系中的表示;ddg为b系原点相对g系原点的矢径在g系中的表示;dgi为g系原点相对i系原点的矢径在i系中的表示;di为b系原点相对i系原点的矢径在i系中的表示;drk为气动舵k本体系原点在b系中的未变形矢径;Fbg为飞行器弹性振动对飞行器转动的耦合系数矩阵;Ffrk为飞行器弹性振动对气动舵k转动的耦合系数矩阵;g为重力加速度;Jbb为飞行器相对b系的转动惯量矩阵;Jrk为气动舵k相对本体坐标系的转动惯量矩阵;M为飞行器的总质量;Mrk为气动舵k的质量;
为飞行器弹性振动的等效质量矩阵;![]()
Pbg为飞行器转动对飞行器平动的耦合系数矩阵;Pf为飞行器弹性振动对飞行器平动的耦合系数矩阵;Prk为气动舵k转动对飞行器平动的耦合系数矩阵;Qb为作用于飞行器的气动力矩和推力产生的力矩;Qbc为发动机推力产生的力矩在飞行器本体坐标系下的表示;Qbq为作用于飞行器上的气动力矩在飞行器本体坐标系下的表示;Qdc为发动机推力在飞行器本体坐标系下的表示;Qdq为气动力在飞行器本体坐标系下的表示;Qrk为气动舵的驱动力矩和作用于气动舵上的干扰力矩对弹性振动产生的广义力;
为气动分布力对飞行器弹性振动模态坐标的广义力;qdq为气动分布力在飞行器体坐标系上的表示;Rbg_rk为气动舵k转动对飞行器转动的耦合系数矩阵;rb为飞行器机身上任一质点dm相对i系的矢径;rrk为气动舵k上任一质点dm相对i系的矢径;
为飞行器本体的模态振型函数矩阵;
为
在飞行器本体与发动机连接点处的值;
为
在飞行器本体与气动舵k连接点处的值;ηb为飞行器弹性振动的模态坐标列向量;ρb为飞行器机身上任一质点dm相对b系坐标原点的未变形矢径;ρrk为气动舵k上任一质点dm相对本体系的矢径;ρrk0为气动舵k的质心在其本体系中的表示;ωbg_b为b系相对g系的角速度在b系中的表示;ωbi_b=ωbg_b+Abiωgi_i为b系相对i系的角速度在b系中的表示;ωgi_i为g系相对i系的角速度在i系中的表示;ωrkb_rk为rk系相对b系的角速度在rk系中的表示;ωrkg_rk=ωrkb_rk+Arkbωbg_b为rk系相对g系的角速度在rk系中的表示;ωrki_rk=ωrkb_rk+Arkb(ωbg_b+Abiωgi_i)为rk系相对i系的角速度在rk系中的表示;Λb为飞行器弹性振动的特征频率对角阵;ξb为飞行器弹性振动的粘性阻尼对角阵;i系为惯性坐标系,b系为飞行器本体坐标系,rk系为气动舵k本体坐标系,g系为地面坐标系;(2)利用步骤(1)得到的高阶模型,根据飞行器位置、速度、姿态、角速度、弹性振动模态坐标以及模态坐标速度的初始值以及除这些量以外的所有其他参数的设计值进行飞行器刚体弹性耦合动力学仿真,得到飞行器位置、速度、姿态、角速度以及弹性振动模态坐标随时间变化的关系曲线,用于从动力学角度验证飞行器设计的合理性。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国运载火箭技术研究院,未经中国运载火箭技术研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201610847855.7/,转载请声明来源钻瓜专利网。