CN112904898B - 旋转弹箭非定常气动响应特性评估方法和系统 - Google Patents
旋转弹箭非定常气动响应特性评估方法和系统 Download PDFInfo
- Publication number
- CN112904898B CN112904898B CN202110120733.9A CN202110120733A CN112904898B CN 112904898 B CN112904898 B CN 112904898B CN 202110120733 A CN202110120733 A CN 202110120733A CN 112904898 B CN112904898 B CN 112904898B
- Authority
- CN
- China
- Prior art keywords
- moment
- damping
- angle
- alpha
- static
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000004044 response Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000013016 damping Methods 0.000 claims abstract description 123
- 230000003068 static effect Effects 0.000 claims abstract description 99
- 238000004458 analytical method Methods 0.000 claims abstract description 90
- 230000033001 locomotion Effects 0.000 claims abstract description 15
- 238000013461 design Methods 0.000 claims abstract description 14
- 230000008878 coupling Effects 0.000 claims description 12
- 238000010168 coupling process Methods 0.000 claims description 12
- 238000005859 coupling reaction Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 9
- 230000006641 stabilisation Effects 0.000 claims description 6
- 238000011105 stabilization Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 4
- 239000011248 coating agent Substances 0.000 claims description 2
- 238000000576 coating method Methods 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 7
- 238000012360 testing method Methods 0.000 abstract description 4
- 230000002159 abnormal effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000000295 complement effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005312 nonlinear dynamic Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000006880 cross-coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/107—Simultaneous control of position or course in three dimensions specially adapted for missiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/08—Aerodynamic models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Fluid Mechanics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种旋转弹箭非定常气动响应特性评估方法和系统,包括:步骤1:进行非定常气动力建模;步骤2:进行气动非线性分析,包括静稳定力矩非线性分析、Magnus力矩非线性分析和阻尼力矩非线性分析;步骤3:进行动态稳定性分析,包括平衡点附近线性化分析和基于非线性模型的α‑β相平面流型分析。本发明解决了传统线性化无法充分考虑非定常气动特性、进而评估失效的难题,可直接用于指导旋转弹箭气动设计,避免飞行试验中气动非定常效应导致异常运动现象的风险。
Description
技术领域
本发明涉及弹箭飞行技术领域,具体地,涉及一种旋转弹箭非定常气动响应特性评估方法和系统。
背景技术
常规导弹的动态响应特性采用小扰动线性化方法进行分析,弹体动力学可简化为二阶系统,由时间常数和阻尼系数即可描述弹体的动态响应情况。类似的方法也被用于旋转弹箭动态响应特性分析。但是,由于弹体旋转会引起显著的气动非定常效应,气动力和力矩非线性度很高,采用线性化方法通常难以准确描述旋转弹箭弹体的实际响应情况。尤其是对于新型的高性能旋转弹箭,设计静稳定度小、阻尼非线性严重,叠加旋转弹箭特有的Magnus效应,使得其气动响应特性相当复杂。在以往的飞行试验中,旋转弹箭会出现不收敛的锥形运动,导致脱靶量增大,试验失败。按照传统方法难以给出有说服力的解释。
根据非定常响应特性对气动设计结果进行评估是旋转弹箭总体设计的重要内容,直接影响弹箭的实际性能。因此,需要发展一种基于非定常气动预测的旋转弹箭动态响应评估方法,用于指导总体设计,避免极限环等非期望运动发生,从而提高弹箭的飞行性能。
专利文献CN107220403A(申请号:CN201710259360.7)公开了一种飞行器弹性模态的控制关联建模方法,包括如下步骤:建立飞行器弹性结构模态的机理模型,获得弹性模态的动态响应数据,基于梁单元确立弹性模态的代理模型形式,采用系统辨识方法整定系数,获得弹性模态的控制关联模型,进而基于拟合优度验证所构建代理模型的有效性,分析弹性模态对于飞行器本体开环和闭环性能的影响,以及本体特性与弹性模态之间的关联关系,找到弹性模态对于系统总体以及控制系统设计的影响规律。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种旋转弹箭非定常气动响应特性评估方法和系统。
根据本发明提供的旋转弹箭非定常气动响应特性评估方法,包括:
步骤1:进行非定常气动力建模,模型的横侧向力矩由静稳定性力矩、Magnus力矩、阻尼力矩和控制力矩组成,静稳定力矩和Magnus力矩的多项式拟合次数不小于2次,阻尼力矩的多项式拟合次数不小于3次;
步骤2:进行气动非线性分析,包括静稳定力矩非线性分析、Magnus力矩非线性分析和阻尼力矩非线性分析,根据静稳定力矩非线性确定最大过载能力,根据Magnus力矩非线性确定静态耦合度,根据阻尼力矩非线性确定极限环状态;
步骤3:进行动态稳定性分析,包括平衡点附近线性化分析和基于非线性模型的α-β相平面流型分析,根据线性化分析获得动态稳定的参数设计边界,根据α-β相平面流型分析获得动态稳定的初值稳定域。
优选的,非定常气动力模型具有如下形式:
其中,Fy和Fz分别表示横侧向力,My和Mz分别表示横侧向力矩;α和β分别表示准攻角和准侧滑角;ωY1和ωZ1分别表示偏航角速度和俯仰角速度;δY1和δZ1分别表示偏航舵偏角和俯仰舵偏角;Q表示动压;S表示特征面积;L表示特征长度;V表示来流速度;表示法向力系数;表示Magnus力系数;表示舵偏诱导法向力系数;表示静稳定力矩相关系数;表示Magnus力矩系相关数;表示阻尼力矩相关系数;表示控制力矩系数;αφ表示全攻角,且αφ=arccos(cosαcosβ)。
优选的,静稳定力矩非线性分析包括静稳定、零攻角附近静不稳定和完全静不稳定;
静稳定:
动态响应攻角收敛到零点;
零攻角附近静不稳定,即存在αequ,使得动态响应攻角收敛至αequ;
完全静不稳定:
动态响应攻角发散;
对于静稳定和零攻角附近静不稳定的情况,获得最大配平攻角,得到弹体的最大过载:
其中δmax表示最大舵偏角,αmax表示最大舵偏角对应的平衡攻角,Fmax表示最大法向力,最大法向力与弹体重力的比即为最大过载;
对于完全静不稳定的情况,最大过载由失速攻角大小决定。
优选的,Magnus力矩非线性分析包括力矩是否随攻角变化发生符号改变和Magnus力矩与静稳定力矩的最大比值;
Magnus力矩符号改变则锥形运动的方向改变,Magnus力矩与静稳定力矩的比值越大,耦合度越高,对系统的平衡点影响越大。
优选的,阻尼力矩非线性分析包括负阻尼、负阻尼→正阻尼和负阻尼→正阻尼→负阻尼;
负阻尼:
使得系统收敛至平衡点位置;
负阻尼→正阻尼:
在负阻尼区域,弹体的收敛特性与静态分析一致,在正阻尼区域,弹体呈现发散状态;
负阻尼→正阻尼→负阻尼:
当αcirc1和αcirc2的值比较接近时,弹体呈现负阻尼状态,出现一个类极限环的动态收敛区域,当αcirc1和αcirc2的值相差较大时,正阻尼区域起主导作用,导致弹体呈现发散状态。
根据本发明提供的旋转弹箭非定常气动响应特性评估系统,包括:
模块M1:进行非定常气动力建模,模型的横侧向力矩由静稳定性力矩、Magnus力矩、阻尼力矩和控制力矩组成,静稳定力矩和Magnus力矩的多项式拟合次数不小于2次,阻尼力矩的多项式拟合次数不小于3次;
模块M2:进行气动非线性分析,包括静稳定力矩非线性分析、Magnus力矩非线性分析和阻尼力矩非线性分析,根据静稳定力矩非线性确定最大过载能力,根据Magnus力矩非线性确定静态耦合度,根据阻尼力矩非线性确定极限环状态;
模块M3:进行动态稳定性分析,包括平衡点附近线性化分析和基于非线性模型的α-β相平面流型分析,根据线性化分析获得动态稳定的参数设计边界,根据α-β相平面流型分析获得动态稳定的初值稳定域。
优选的,非定常气动力模型具有如下形式:
其中,Fy和Fz分别表示横侧向力,My和Mz分别表示横侧向力矩;α和β分别表示准攻角和准侧滑角;ωY1和ωZ1分别表示偏航角速度和俯仰角速度;δY1和δZ1分别表示偏航舵偏角和俯仰舵偏角;Q表示动压;S表示特征面积;L表示特征长度;V表示来流速度;表示法向力系数;表示Magnus力系数;表示舵偏诱导法向力系数;表示静稳定力矩相关系数;表示Magnus力矩系相关数;表示阻尼力矩相关系数;表示控制力矩系数;αφ表示全攻角,且αφ=arccos(cosαcosβ)。
优选的,静稳定力矩非线性分析包括静稳定、零攻角附近静不稳定和完全静不稳定;
静稳定:
动态响应攻角收敛到零点;
零攻角附近静不稳定,即存在αequ,使得动态响应攻角收敛至αequ;
完全静不稳定:
动态响应攻角发散;
对于静稳定和零攻角附近静不稳定的情况,获得最大配平攻角,得到弹体的最大过载:
其中δmax表示最大舵偏角,αmax表示最大舵偏角对应的平衡攻角,Fmax表示最大法向力,最大法向力与弹体重力的比即为最大过载;
对于完全静不稳定的情况,最大过载由失速攻角大小决定。
优选的,Magnus力矩非线性分析包括力矩是否随攻角变化发生符号改变和Magnus力矩与静稳定力矩的最大比值;
Magnus力矩符号改变则锥形运动的方向改变,Magnus力矩与静稳定力矩的比值越大,耦合度越高,对系统的平衡点影响越大。
优选的,阻尼力矩非线性分析包括负阻尼、负阻尼→正阻尼和负阻尼→正阻尼→负阻尼;
负阻尼:
使得系统收敛至平衡点位置;
负阻尼→正阻尼:
在负阻尼区域,弹体的收敛特性与静态分析一致,在正阻尼区域,弹体呈现发散状态;
负阻尼→正阻尼→负阻尼:
当αcirc1和αcirc2的值比较接近时,弹体呈现负阻尼状态,出现一个类极限环的动态收敛区域,当αcirc1和αcirc2的值相差较大时,正阻尼区域起主导作用,导致弹体呈现发散状态。
与现有技术相比,本发明具有如下的有益效果:
(1)本发明的旋转弹箭非定常气动响应特性评估方法,从非定常气动力建模出发,结合物理实际,给出了横侧向气动力矩多项式拟合的阶次约束,并通过具体分析,得到了非线性气动力的影响情况。这些分析可初步确定旋转弹箭的过载能力、横侧向耦合强度、攻角-侧滑角平面响应流型的基本情况。根据气动力建模和分析,开展系统的动态稳定性分析,进而得到旋转弹箭气动响应的稳定性评估结果。三个步骤相辅相成,解决了传统线性化无法充分考虑非定常气动特性、进而评估失效的难题;
(2)本发明的旋转弹箭非定常气动响应特性评估方法,可直接用于指导旋转弹箭气动设计,避免飞行试验中气动非定常效应导致异常运动现象的风险。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1是本发明较佳实施例的旋转导弹稳定性边界计算方法流程图;
图2是(α,β)平面运动流型Magnus力矩约为静稳定力矩10%的示意图;
图3是(α,β)平面运动流型Magnus力矩约为静稳定力矩50%的示意图;
图4是(α,β)平面运动流型小攻角负阻尼、大攻角正阻尼情况的示意图;
图5是(α,β)平面运动流型阻尼随攻角增加呈负、正、负情况的示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
实施例:
本发明的旋转弹箭非定常气动响应特性评估方法包括非定常气动力建模、气动非线性分析、动态稳定性分析三个主要步骤。
步骤1、非定常气动力建模的横侧向力矩由静稳定性力矩、Magnus力矩、阻尼力矩和控制力矩等组成,静稳定力矩和Magnus力矩的多项式拟合次数不小于2次,阻尼力矩的多项式拟合次数不小于3次。
步骤2、气动非线性分析包括静稳定力矩非线性分析、Magnus力矩非线性分析和阻尼力矩非线性分析,根据静稳定力矩非线性确定最大过载能力,根据Magnus力矩非线性确定静态耦合度,根据阻尼力矩非线性确定极限环状态。
步骤3、动态稳定性分析包含平衡点附近线性化分析和基于非线性模型的α-β相平面流型分析;线性化分析获得动态稳定的参数设计边界,α-β相平面流型分析获得动态稳定的初值稳定域。
如图1所示,非定常气动力模型为气动非线性分析和动态稳定性分析的基础,气动非线性分析可初步确定旋转弹箭的过载能力、横侧向耦合强度、攻角-侧滑角平面响应流型的基本情况,同时也作为动态稳定性分析输入。三个步骤相辅相成,解决了传统线性化无法充分考虑非定常气动特性、进而评估失效的难题。
本发明的旋转弹箭非定常气动响应特性评估方法,非定常气动力模型具有如下形式:
其中,Fy和Fz分别表示横侧向力,My和Mz分别表示横侧向力矩;α和β分别表示准攻角和准侧滑角;ωY1和ωZ1分别表示偏航角速度和俯仰角速度;δY1和δZ1分别表示偏航舵偏角和俯仰舵偏角;Q表示动压;S表示特征面积;L表示特征长度;V表示来流速度;表示法向力系数;表示Magnus力系数;表示舵偏诱导法向力系数;表示静稳定力矩相关系数;表示Magnus力矩系相关数;表示阻尼力矩相关系数;表示控制力矩系数;αφ表示全攻角,且αφ=arccos(cosαcosβ)。
非定常气动力的横侧向力包含法向力、Magnus力和舵偏诱导法向力等,横侧向力矩包含静稳定性力矩、Magnus力矩、阻尼力矩和控制力矩等。静稳定力矩和Magnus力矩的多项式拟合次数不小于2次,阻尼力矩的多项式拟合次数不小于3次,即和的n不小于1,的n不小于2。
本发明的旋转弹箭非定常气动响应特性评估方法,气动非线性分析包括静稳定力矩非线性分析、Magnus力矩非线性分析和阻尼力矩非线性分析。
静稳定力矩非线性分析可分为三种情况,分别为静稳定、零攻角附近静不稳定和完全静不稳定。
静稳定:
动态响应攻角收敛到零点。
零攻角附近静不稳定,即存在αequ,使得动态响应攻角收敛至αequ。
完全静不稳定:
动态响应攻角发散。
对于静稳定和零攻角附近静不稳定的情况,可以获得最大配平攻角,得到弹体的最大过载。
其中δmax表示最大舵偏角,αmax表示最大舵偏角对应的平衡攻角,Fmax表示最大法向力,最大法向力与弹体重力的比即为最大过载
对于完全静不稳定的情况,最大过载由失速攻角大小决定。
Magnus力矩非线性分析包括力矩是否随攻角变化发生符号改变和Magnus力矩与静稳定力矩的最大比值。
Magnus力矩作用的结果有两个方面,一是横侧向耦合,使得弹体运动呈锥形运动形式;二是通过交叉耦合作用影响系统的平衡点分布。Magnus力矩符号改变则锥形运动的方向改变,Magnus力矩与静稳定力矩的比值越大,耦合度越高,对系统的平衡点影响越大。
图2和图3为零攻角附近静不稳定设计旋转弹箭α-β相平面流型情况,图2为Magnus力矩约为静稳定力矩10%的情况,极限环5.13°,图3为Magnus力矩约为静稳定力矩50%的情况,极限环6.22°。不考虑Magnus力矩作用的理论平衡攻角为5°,可见Magnus力矩越大对平衡点位置影响越大。
阻尼力矩非线性分析可分为三种情况,分别为负阻尼、负阻尼→正阻尼和负阻尼→正阻尼→负阻尼。
负阻尼:
可使得系统收敛至平衡点位置。
负阻尼→正阻尼:
在负阻尼区域,弹体的收敛特性与静态分析一致,在正阻尼区域,弹体通常会呈现发散状态。
负阻尼→正阻尼→负阻尼:
当αcirc1和αcirc2的值比较接近时,弹体主要呈现负阻尼状态,但会出现一个类极限环的动态收敛区域,当αcirc1和αcirc2的值相差较大时,正阻尼区域会其主导作用,通常将导致弹体呈现发散状态。
图4为小攻角负阻尼、大攻角正阻尼的情况,在图中,初值αφ<4.535时,系统呈现收敛状态,如两条实线所示;当初值αφ≥4.535时,系统呈现发散状态,如两条虚线所示。
图5为阻尼随攻角增加呈负、正、负的情况,在图中,初值αφ较小时,系统呈现收敛状态,如虚线所示;当初值αφ增加后,系统呈现极限环响应状态,如实线所示。
本发明的旋转弹箭非定常气动响应特性评估方法,动态稳定性分析基于旋转弹箭非线性动力学模型开展,包含平衡点附近线性化分析和基于非线性模型的α-β相平面流型分析。
平衡点附近线性化首先求解系统的平衡点,采用一阶Lyapunov线性化方法得到平衡点附近的线性模型,然后根据线性理论获得系统稳定条件,进而得到主导气动参数的设计边界。α-β相平面流型分析直接基于非线性动力学模型,采用数值方法获得动态稳定的初值稳定域。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以,本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件,而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构;也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (6)
1.一种旋转弹箭非定常气动响应特性评估方法,其特征在于,包括:
步骤1:进行非定常气动力建模,模型的横侧向力矩由静稳定性力矩、Magnus力矩、阻尼力矩和控制力矩组成,静稳定力矩和Magnus力矩的多项式拟合次数不小于2次,阻尼力矩的多项式拟合次数不小于3次;
步骤2:进行气动非线性分析,包括静稳定力矩非线性分析、Magnus力矩非线性分析和阻尼力矩非线性分析,根据静稳定力矩非线性确定最大过载能力,根据Magnus力矩非线性确定静态耦合度,根据阻尼力矩非线性确定极限环状态;
步骤3:进行动态稳定性分析,包括平衡点附近线性化分析和基于非线性模型的α-β相平面流型分析,根据线性化分析获得动态稳定的参数设计边界,根据α-β相平面流型分析获得动态稳定的初值稳定域;
非定常气动力模型具有如下形式:
其中,Fy和Fz分别表示横侧向力,My和Mz分别表示横侧向力矩;α和β分别表示准攻角和准侧滑角;ωY1和ωZ1分别表示偏航角速度和俯仰角速度;δY1和δZ1分别表示偏航舵偏角和俯仰舵偏角;Q表示动压;S表示特征面积;L表示特征长度;V表示来流速度;表示法向力系数;表示Magnus力系数;表示舵偏诱导法向力系数;表示静稳定力矩相关系数;表示Magnus力矩系相关数;表示阻尼力矩相关系数;表示控制力矩系数;αφ表示全攻角,且αφ=arccos(cosαcosβ);
静稳定力矩非线性分析包括静稳定、零攻角附近静不稳定和完全静不稳定;
静稳定:
动态响应攻角收敛到零点;
零攻角附近静不稳定,即存在αequ,使得动态响应攻角收敛至αequ;
完全静不稳定:
动态响应攻角发散;
对于静稳定和零攻角附近静不稳定的情况,获得最大配平攻角,得到弹体的最大过载:
其中δmax表示最大舵偏角,αmax表示最大舵偏角对应的平衡攻角,Fmax表示最大法向力,最大法向力与弹体重力的比即为最大过载;
对于完全静不稳定的情况,最大过载由失速攻角大小决定。
2.根据权利要求1所述的旋转弹箭非定常气动响应特性评估方法,其特征在于,Magnus力矩非线性分析包括力矩是否随攻角变化发生符号改变和Magnus力矩与静稳定力矩的最大比值;
Magnus力矩符号改变则锥形运动的方向改变,Magnus力矩与静稳定力矩的比值越大,耦合度越高,对系统的平衡点影响越大。
4.一种旋转弹箭非定常气动响应特性评估系统,其特征在于,包括:
模块M1:进行非定常气动力建模,模型的横侧向力矩由静稳定性力矩、Magnus力矩、阻尼力矩和控制力矩组成,静稳定力矩和Magnus力矩的多项式拟合次数不小于2次,阻尼力矩的多项式拟合次数不小于3次;
模块M2:进行气动非线性分析,包括静稳定力矩非线性分析、Magnus力矩非线性分析和阻尼力矩非线性分析,根据静稳定力矩非线性确定最大过载能力,根据Magnus力矩非线性确定静态耦合度,根据阻尼力矩非线性确定极限环状态;
模块M3:进行动态稳定性分析,包括平衡点附近线性化分析和基于非线性模型的α-β相平面流型分析,根据线性化分析获得动态稳定的参数设计边界,根据α-β相平面流型分析获得动态稳定的初值稳定域;
非定常气动力模型具有如下形式:
其中,Fy和Fz分别表示横侧向力,My和Mz分别表示横侧向力矩;α和β分别表示准攻角和准侧滑角;ωY1和ωZ1分别表示偏航角速度和俯仰角速度;δY1和δZ1分别表示偏航舵偏角和俯仰舵偏角;Q表示动压;S表示特征面积;L表示特征长度;V表示来流速度;表示法向力系数;表示Magnus力系数;表示舵偏诱导法向力系数;表示静稳定力矩相关系数;表示Magnus力矩系相关数;表示阻尼力矩相关系数;表示控制力矩系数;αφ表示全攻角,且αφ=arccos(cosαcosβ);
静稳定力矩非线性分析包括静稳定、零攻角附近静不稳定和完全静不稳定;
静稳定:
动态响应攻角收敛到零点;
零攻角附近静不稳定,即存在αequ,使得动态响应攻角收敛至αequ;
完全静不稳定:
动态响应攻角发散;
对于静稳定和零攻角附近静不稳定的情况,获得最大配平攻角,得到弹体的最大过载:
其中δmax表示最大舵偏角,αmax表示最大舵偏角对应的平衡攻角,Fmax表示最大法向力,最大法向力与弹体重力的比即为最大过载;
对于完全静不稳定的情况,最大过载由失速攻角大小决定。
5.根据权利要求4所述的旋转弹箭非定常气动响应特性评估系统,其特征在于,Magnus力矩非线性分析包括力矩是否随攻角变化发生符号改变和Magnus力矩与静稳定力矩的最大比值;
Magnus力矩符号改变则锥形运动的方向改变,Magnus力矩与静稳定力矩的比值越大,耦合度越高,对系统的平衡点影响越大。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110120733.9A CN112904898B (zh) | 2021-01-28 | 2021-01-28 | 旋转弹箭非定常气动响应特性评估方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110120733.9A CN112904898B (zh) | 2021-01-28 | 2021-01-28 | 旋转弹箭非定常气动响应特性评估方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112904898A CN112904898A (zh) | 2021-06-04 |
CN112904898B true CN112904898B (zh) | 2022-11-29 |
Family
ID=76119950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110120733.9A Active CN112904898B (zh) | 2021-01-28 | 2021-01-28 | 旋转弹箭非定常气动响应特性评估方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112904898B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114896830B (zh) * | 2022-07-14 | 2022-11-08 | 中国空气动力研究与发展中心计算空气动力研究所 | 导弹非线性非定常气动力微分方程模型辨识方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5379968A (en) * | 1993-12-29 | 1995-01-10 | Raytheon Company | Modular aerodynamic gyrodynamic intelligent controlled projectile and method of operating same |
CN109540452A (zh) * | 2018-11-21 | 2019-03-29 | 中国航空工业集团公司沈阳空气动力研究所 | 一种旋转弹箭三自由度角运动模拟试验装置 |
CN110750836A (zh) * | 2019-10-15 | 2020-02-04 | 北京电子工程总体研究所 | 一种基于稳态锥动的飞行器俯仰阻尼力矩计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111323198B (zh) * | 2020-03-10 | 2022-01-04 | 中国航天空气动力技术研究院 | 一种船尾外形旋转弹马格努斯力滚转导数获取方法 |
-
2021
- 2021-01-28 CN CN202110120733.9A patent/CN112904898B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5379968A (en) * | 1993-12-29 | 1995-01-10 | Raytheon Company | Modular aerodynamic gyrodynamic intelligent controlled projectile and method of operating same |
CN109540452A (zh) * | 2018-11-21 | 2019-03-29 | 中国航空工业集团公司沈阳空气动力研究所 | 一种旋转弹箭三自由度角运动模拟试验装置 |
CN110750836A (zh) * | 2019-10-15 | 2020-02-04 | 北京电子工程总体研究所 | 一种基于稳态锥动的飞行器俯仰阻尼力矩计算方法 |
Non-Patent Citations (2)
Title |
---|
杨永强 等.旋转导弹非线性动力学建模方法研究.《上海航天》.2017,第34卷 * |
钱龙 等.固定鸭舵修正弹非线性与非定常气动数值模拟.《弹道学报》.2020,第32卷(第1期), * |
Also Published As
Publication number | Publication date |
---|---|
CN112904898A (zh) | 2021-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106078742B (zh) | 一种针对带有输出约束的柔性机械臂的振动控制方法 | |
Tewari et al. | Aeroservoelasticity | |
CN113806871B (zh) | 一种考虑结构非线性的柔性飞行动力学建模方法 | |
CN109062055A (zh) | 一种基于Back-stepping鲁棒自适应动态面的近空间飞行器控制系统 | |
Andrievsky et al. | Simple adaptive control for airfoil flutter suppression | |
CN112904898B (zh) | 旋转弹箭非定常气动响应特性评估方法和系统 | |
An et al. | Neural adaptive control of air-breathing hypersonic vehicles robust to actuator dynamics | |
CN111159812A (zh) | 一种吸气式高超声速无人机损伤特性分析方法 | |
Ahmadian et al. | A novel L1 adaptive-hybrid control with guaranteed stability for a class of uncertain nonlinear systems: A case study on SA330 Puma | |
Karpel | Sensitivity derivatives of flutter characteristics and stability margins for aeroservoelastic design | |
Li et al. | Adaptive maneuver load alleviation via recurrent neural networks | |
Kazemi-Lari et al. | Flexural-torsional flutter analysis of a deep cantilever beam subjected to a partially distributed lateral force | |
Emami et al. | Robustness investigation of a ducted-fan aerial vehicle control, using linear, adaptive, and model predictive controllers | |
Tokunaga et al. | Flight evaluation of fault-tolerant control system using simple adaptive control method | |
Roy et al. | Robust altitude control of an unmanned autonomous helicopter using backstepping | |
Gregory | Dynamic inversion to control large flexible transport aircraft | |
CN112015196A (zh) | 一种姿控系统限幅值设计方法、存储介质及服务器 | |
Kasumov | Modeling of helicopter behavior dynamics under the controlling action of the stabilization and control system | |
Ju et al. | New Improved g Method for Flutter Solution | |
Yang et al. | Adaptive regulation of hypersonic vehicle systems with partial nonlinear parametrization | |
Antequera et al. | A helicopter control based on eigenstructure assignment | |
Pirzadeh et al. | Design of generalized predictive control for the stabilizing loop from a two-axis gimbal seeker, considering cross-coupling in between two channels | |
Wang et al. | Virtual flight simulation of the basic finner projectile based on fuzz control | |
CN112527007B (zh) | 抑制大迎角引发机翼摇滚的直接自适应模糊逻辑控制方法 | |
Xu et al. | MRAC control with prior model knowledge for asymmetric damaged aircraft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |