CN113919081B - 一种考虑惯性耦合的柔性飞行动力学建模与分析方法 - Google Patents
一种考虑惯性耦合的柔性飞行动力学建模与分析方法 Download PDFInfo
- Publication number
- CN113919081B CN113919081B CN202111345980.5A CN202111345980A CN113919081B CN 113919081 B CN113919081 B CN 113919081B CN 202111345980 A CN202111345980 A CN 202111345980A CN 113919081 B CN113919081 B CN 113919081B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- equation
- expression
- inertial
- generalized
- 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
- 230000008878 coupling Effects 0.000 title claims abstract description 24
- 238000010168 coupling process Methods 0.000 title claims abstract description 24
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 24
- 238000004458 analytical method Methods 0.000 title abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 18
- 230000014509 gene expression Effects 0.000 claims description 24
- 238000005381 potential energy Methods 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000005489 elastic deformation Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 238000010205 computational analysis Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000009471 action Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 239000012141 concentrate Substances 0.000 description 3
- 238000005096 rolling process Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑惯性耦合的柔性飞行动力学建模与分析方法,从准坐标系下表达的拉格朗日方程出发,建立能够反映惯性耦合的柔性飞行器飞行动力学关系方程;针对忽略惯性耦合的柔性飞行器动力学模型,利用准定常气动力模型建立了与采用刚体假设的飞行动力学模型类似的柔性飞行动力学模型,保证了一定计算精度条件下,提高计算效率。本发明的方法可以用于采用细长体布局的高速飞行器、大型运输机和超大展弦比无人机飞行动力学建模及后续动力学分析,从而为相关设计提供指导思想,且具有一定的工程应用价值。
Description
技术领域
本发明属于飞行器设计技术领域,具体涉及一种考虑惯性耦合的柔性飞行动力学建模与分析方法,可以用于采用细长体布局的高速飞行器、大型运输机、超大展弦比无人机和微小型等柔性飞行器飞行动力学建模和动力学分析及后续的控制设计。
背景技术
随着空天科学技术的不断进步,飞行器速度、高度、航时的记录不断突破,创新气动布局不断涌新,轻质材料广泛使用,对飞行器设计工程实践和学科发展提出新的要求。采用机身刚体假设的传统飞行力学建模无法满足结构柔性飞行动力学建模与控制设计问题的需要。柔性飞行器飞行动力学建模是实现动力学分析和控制设计的基础。对于柔性飞行器,飞行动力学模型需要反映飞行器大范围运动,又能反映结构的动力学行为,同时还需要能够充分利用已有的分析方式和可获取的气动、结构等参数的数据源。
发明内容
本发明针对柔性飞行器动力学建模问题,利用准坐标系下的拉格朗日方程建立了柔性飞行器飞行动力学模型,其中能够考虑惯性耦合,方便整理成用于分析和仿真的形式,分析采用平均轴系假设和忽略惯性耦合后,建模方法之间的退化形式,以及刚体运动和弹性运动之间的耦合机理。具体的技术方案为:
一种考虑惯性耦合的柔性飞行动力学建模与分析方法,包括以下步骤:
(1)采用以下假设:不考虑地球曲率和自转;飞机结构描述为一系列的集中质量,相应的集中质量和惯性矩为mi和Ji;变形满足线弹性理论和小变形假设;为了方便结构模型处理,认为能够获得满足正交化假设的结构模态;
(2)定义以下坐标系:1.惯性系(xe,ye,ze):与大地固连,根据平面大地假设可认为是惯性系,坐标原点oe;2.体轴系(xb,yb,zb):未变形状态的机体坐标系,原点ob取在飞行器未变形时的质心上,xb指向头部为正,yb垂直于xb,方向指向机身右侧,zb按右手定则确定;3.当地坐标系(xi,yi,zi):在未变形机体上的每个节点,相对体轴系原点的位置向量为ri,坐标轴与体轴系平行。
(3)建立飞行器运动学关系方程。考虑未变形状态,与刚体飞行器固联的体轴系在惯性系的位置向量R0e=(x0e,y0e,z0e)T,用欧拉角φ,θ,ψ表达从惯性坐标系到体轴系的姿态变换Θ=(φ,θ,ψ)T,速度和角速度在体轴系的表达为Vb=(u,v,w)T,Ωb=(p,q,r)T,描述飞行器质点运动和姿态运动的运动学关系方程为
其中,φ、θ、ψ分别为滚转角、俯仰角和偏航角;p、q、r分别为滚转角速度、俯仰角速度和偏航角速度;u、v、w分别为飞行器速度在体轴系三个轴上的投影;
体轴系原点位置对于每个变形后集中质量位置在体轴系中的表达为
其中R0为惯性系中体轴系坐标原点的位置向量在体轴系的表达,di为当地坐标系中的变形,si为集中质量在当地坐标系中的位置向量,为扭转变形。T表示变换矩阵,描述位置向量si从未变形位置到变形后位置的扭转。在线性假设下,对其泰勒展开取一阶
其中表示/>的反对称矩阵。
集中质量位移Ri:
集中质量速度
集中质量的转动速度Ωi由体轴系的转动Ωb和变形带来的转动共同组成:
节点的弹性位移可以表达为平移和转动,根据假设条件,可以满足正交条件的自由-自由边界条件下的模态表达,设弹性变形的广义坐标为ηE,节点弹性位移可以表达为:
其中,di为节点的运动变形,为节点的扭转变形,/>和/>分别为广义坐标ηE中平动和转动自由度对应的广义坐标模态;
(4)建立飞行器动力学关系方程。采用准坐标系下表达的拉格朗日方程,建立动力学方程:
包含平动自由度方程:
其中Pg=[F M]T为节点上的载荷,包含非保守外力。非保守力包含气动力和其他外力,其中气动力可以采用准定常方式获得气动力表达。为广义模态坐标下的广义模态,为/>
转动自由度方程:
在计算分析中,h,J,/>每一步都需要重新计算。对于变形后的飞行器,惯性张量展开为
将表达式中的每一项编号表示为:
A1不随时间变化,来自坐标轴偏移。
有
e为单位向量,e1=[1 0 0]T,e2=[0 1 0]T,e3=[0 0 1]T;
有
其中,
ai1=sizΦgiyEr-siyΦgizEr
ai2=sixΦgizEr-sizΦgixEr
ai3=SiyΦgixEr-SixΦgiyEr
有
A5=A4T
有
有
A8=A6T
A9=A7T
定义
惯性张量J的表达式整理为:
J对时间求导:
类似的,动量矩h可以用以下公式进行计算:
表达式中的各项编号,可以写为:有
其中,
特别的,为方便模型的求解,需要对系统拉格朗日方程的原始表达形式进行解释:
其中L为拉格朗日算子,L=Ek-Ep,Ek、Ep分别为动能和势能,Qt,Qr,QE为非保守外力,满足:虚功
动能利用下式进行计算:
对于动能项Ek求解,需要计算每个集中质量的平动动能Ekt和转动动能Ekr:
其中,Ri为集中质量的位置矢量,Ji为集中质量的转动惯量;
结合实用平均轴系假设,求得平动动能:
其中平动动能公式最后两项体现了刚性运动与弹性变形之间的耦合。
转动动能:
令总动能可以求解、整理为
利用广义坐标和模态表示变形,第二项和第五项整理为
因此,总动能表达式为
其中,Ri为集中质量的位置矢量,Ji为集中质量的转动惯量;
对于势能项Ep计算,势能包含重力势能和弹性势能,利用线弹性假设,全机势能表达式为
其中,Teb=Tbe -1,ge=[0 0 g],g为当地重力加速度,Kgg为刚度阵,ug为弹性变形;
对于非保守外力求解,主要计算气动力、推力和其他外力,其中气动力表达为剖面升力其中/>为动压,αs为剖面迎角,机翼升力由剖面升力沿展向积分得到,即/>为方便计算气动力会转化为升力系数的形式,为该领域基本问题,在此不赘述。
考虑每个节点上的载荷Pg,i=[Fi Mi]T,非保守力和力矩的虚功用下式计算:
其中δa为体轴系的虚角位移。也可以利用每个节点的虚位移和刚体运动模态表达
其中广义力:
建立飞行器弹性关系方程。通过拉格朗日变量对广义变形坐标求导,并在模型中加入对称形式的阻尼弹性关系方程可以整理为:
其中MEE为广义质量矩阵,KEE为广义刚度矩阵。
本方法得到的飞行器弹性关系方程为:
其中,ηE为弹性变形的广义坐标,ej为单位向量,来源于转动惯量J的表达式各项,有
其中,
来源于动量矩表达式中的各项,有:
(5)利用数值方法对方程其进行求解,计算本发明所述建模方式在初始条件下的时域响应。本方法建立考虑惯性耦合的柔性飞行器飞行动力学方程总结如下:
运动学关系方程:
动力学关系方程:
进一步,通过数值方法计算模型在特定初始条件作用下运动参数的时域响应,分析飞行器模型的稳定性。
本方法利用准坐标系描述的拉格朗日方程推导了柔性飞行器飞行动力学方程,模型采用的假设条件和结构模型便于利用已有的气动参数和结构有限元模型,方程中包含惯性耦合项,如果采用更为苛刻的简化条件,能够获得模型的进一步简化,便于不同模型精度下的分析,适应不同的分析和控制设计需求。如不考虑刚体运动和弹性运动惯性耦合,可以对动力学模型进行简化为常规动力学方程形式:
/>
附图说明
图1为飞行动力学建模流程示意图;
图2为实施例结构及坐标系示意图;
图3为实施例节点1在体轴系内的位移时域响应;
图4为实施例在节点1在体轴系内的位移的时域响应;
图5为实施例全机角速度时域响应曲线;
图6为实施例节点1转角的时域响应。
具体实施方式
下面结合附图和实施例对本发明的考虑惯性耦合的柔性飞行器的飞行动力学建模与分析方法进一步详细说明。
具体步骤如图1所述,本实施例针对虑极端简化的飞机模型(图2),将机身和发动机考虑为欧拉梁连接中部和两端的集中质量,对其进行柔性飞行动力学建模。
第一步,采用以下假设:不考虑地球曲率和自转;飞机结构描述为一系列的集中质量,相应的集中质量和惯性矩为mi和Ji;变形满足线弹性理论和小变形假设;为了方便结构模型处理,认为能够获得满足正交化假设的结构模态;
第二步,对飞行器进行坐标系定义,体轴系原点位于中心处,初始位置与惯性系重合,确定飞行器简化模型模型各参数为
m1=1kg,m2=2kg,m3=1kg,J1=8×10-4kg.m2,J2=2.5×10-3kg.m2,J3=8×10- 4kg.m2
A=1×10-5m2,E=2×109N/m2,l0=1m,I=1×10-8m-4
其中,m1和m3为吊舱质量,J1和J3为吊舱转动惯量,m2为机身质量,J2为机身转动惯量,A为机翼截面积,E为机翼刚度,I为机翼截面惯性矩。
第三步,建立模型的运动学关系方程。当前模型在仅考虑横向平面运动前提下,运动学关系方程简化为
第四步,建立模型的动力学关系方程,弹性变形以向量描述为
将位移用广义坐标形式表达为μg=ΦgEηE
其中ΦgE为欧拉梁模型在自由-自由边界条件下获得的振型。
分析得到模型动力学关系方程如下:
/>
其中,分析外力矩Mx在短时间内的作用,设Mx为
建立模型的弹性关系方程,模型中加入对称形式的阻尼,即
飞行器的弹性关系方程可以表达为:
第五步,通过积分计算本发明所述建模方式在给定力矩作用下运动参数的时域响应。如图3为节点1的位置在体轴系的变化,图4为节点1在y轴的位移随时间的变化。结果显示,在体轴系内观察节点1的位移,可以发现,惯性耦合后的结果在滚转力矩的作用下存在y方向的位移,即伸长。图5所示为全机角速度的时域响应,图6所示为节点1的时域响应。
可见对于研究飞机刚体运动为主的飞行动力学与控制问题,飞行器飞行过程中的载荷分析和计算以及相关的气动弹性问题时,考虑惯性耦合可以提供精度更高的分析结果尤其是存在较为明显的集中质量和做类似高速滚转等特殊机动的情况。
Claims (4)
1.一种考虑惯性耦合的柔性飞行动力学建模与分析方法,其特征在于,包括以下步骤:
(1)对建模的条件进行假设;
(2)对飞行器进行坐标系定义;
(3)建立飞行器运动学关系方程;
(4)建立飞行器动力学关系方程和弹性关系方程;
(5)利用数值方法对方程其进行求解,计算模型在初始条件下的时域响应;
所述步骤(4)具体包括:
采用准坐标系下表达的拉格朗日方程,建立动力学方程:
包含平动自由度方程:
其中Pg=[F M]T为节点上的载荷,包含非保守外力,非保守力包含气动力和其他外力,其中气动力采用准定常方式获得气动力表达,为广义模态坐标下的广义模态,为
转动自由度方程:
在计算分析中,每一步都需要重新计算,对于变形后的飞行器,惯性张量展开为
将表达式中的每一项编号表示为:
A1不随时间变化,来自坐标轴偏移;
有
e为单位向量,e1=[1 0 0]T,e2=[0 1 0]T,e3=[0 0 1]T;
有
其中,
ai1=sizΦgiyEr-siyΦgizEr
ai2=sixΦgizEr-sizΦgixEr
ai3=siyΦgixEr-sixΦgiyEr
有
A5=A4T
有
有
A8=A6T
A9=A7T
定义惯性张量J的表达式整理为:
J对时间求导:
动量矩h用以下公式进行计算:
表达式中的各项编号写为:有
其中,
对系统拉格朗日方程的原始表达形式进行解释:
其中L为拉格朗日算子,L=Ek-Ep,Ek、Ep分别为动能和势能,Qt,Qr,QE为非保守外力,满足:虚功
动能利用下式进行计算:
其中,对于动能Ek求解,需要计算每个集中质量的平动动能Ekt和转动动能Ekr:
结合实用平均轴系假设,求得平动动能:
其中平动动能公式最后两项体现了刚性运动与弹性变形之间的耦合;
转动动能:
令总动能求解、整理为
利用广义坐标和模态表示变形,第二项和第五项整理为
因此,总动能表达式为
对于势能项Ep计算,势能包含重力势能和弹性势能,利用线弹性假设,全机势能表达式为
对于非保守外力求解,计算气动力、推力和其他外力,其中气动力表达为剖面升力其中/>为动压,αs为剖面迎角,机翼升力由剖面升力沿展向积分得到,即/>将气动力转化为升力系数的形式;
考虑每个节点上的载荷Pg,i=[Fi Mi]T,非保守力和力矩的虚功用下式计算:
其中δa为体轴系的虚角位移;
建立飞行器弹性关系方程,通过拉格朗日变量对广义变形坐标求导,并在模型中加入对称形式的阻尼弹性关系方程整理为:
其中MEE为广义质量矩阵,KEE为广义刚度矩阵;
得到的飞行器弹性关系方程为:
其中,ηE为弹性变形的广义坐标,ej为单位向量,来源于转动惯量J的表达式各项,有
ai1=sizΦgiyEr-siyΦgizEr
ai2=sixΦgizEr-sizΦgixEr
其中,ai3=siyΦgixEr-sixΦgiyEr
来源于动量矩表达式中的各项,有:
2.根据权利要求1所述的一种考虑惯性耦合的柔性飞行动力学建模与分析方法,其特征在于,步骤(1)采用以下假设:
不考虑地球曲率和自转;飞机结构描述为一系列的集中质量,相应的集中质量和惯性矩为mi和Ji;变形满足线弹性理论和小变形假设;为了方便结构模型处理,认为能够获得满足正交化假设的结构模态。
3.根据权利要求1所述的一种考虑惯性耦合的柔性飞行动力学建模与分析方法,其特征在于,步骤(2)定义以下坐标系:
惯性系(xe,ye,ze):与大地固连,根据平面大地假设可认为是惯性系,坐标原点oe;
体轴系(xb,yb,zb):未变形状态的机体坐标系,原点ob取在飞行器未变形时的质心上,xb指向头部为正,yb垂直于xb,方向指向机身右侧,zb按右手定则确定;
当地坐标系(xi,yi,zi):在未变形机体上的每个节点,相对体轴系原点的位置向量为ri,坐标轴与体轴系平行。
4.根据权利要求1所述的一种考虑惯性耦合的柔性飞行动力学建模与分析方法,其特征在于,通过准坐标系下表达的拉格朗日方程,建立的飞行器动力学模型中包括动力学关系和弹性关系方程;动力学方程中包含惯性耦合项
其中为飞行器惯性张量J对时间的导数,Ωb为飞行器在体轴系的角速度,h为飞行器动量,/>为飞行器动量对时间的导数;弹性关系方程中包含惯性耦合项中,ηE和/>分别为弹性变形的广义坐标和广义坐标对时间的导数,ej为单位向量,有e1=[1 0 0]T,e2=[01 0]T,e3=[0 0 1]T,/>来源于飞行器转动惯量的表达式中的各项,/>来源于飞行器动量矩表达式中的各项,/>为飞行器在体轴系的角加速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111345980.5A CN113919081B (zh) | 2021-11-15 | 2021-11-15 | 一种考虑惯性耦合的柔性飞行动力学建模与分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111345980.5A CN113919081B (zh) | 2021-11-15 | 2021-11-15 | 一种考虑惯性耦合的柔性飞行动力学建模与分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113919081A CN113919081A (zh) | 2022-01-11 |
CN113919081B true CN113919081B (zh) | 2023-09-19 |
Family
ID=79246624
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111345980.5A Active CN113919081B (zh) | 2021-11-15 | 2021-11-15 | 一种考虑惯性耦合的柔性飞行动力学建模与分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113919081B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116611166B (zh) * | 2023-05-23 | 2024-03-12 | 南京航空航天大学 | 一种考虑初始误差的飞机部件多层级装配偏差分析方法 |
CN117852306B (zh) * | 2024-03-04 | 2024-05-24 | 西北工业大学 | 基于相对动量矩定理的变体高速飞行器动力学建模方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106484969A (zh) * | 2016-09-23 | 2017-03-08 | 中国运载火箭技术研究院 | 一种大包线强机动飞行器动力学高精度仿真方法 |
CN108427322A (zh) * | 2018-03-16 | 2018-08-21 | 南京航空航天大学 | 一种大柔性飞行器基于在线辨识的建模方法 |
CN108959796A (zh) * | 2018-07-18 | 2018-12-07 | 哈尔滨工业大学 | 一种大惯量旋转载荷卫星的刚柔磁耦合动力学建模方法 |
CN111931292A (zh) * | 2020-08-07 | 2020-11-13 | 北京航空航天大学 | 一种翼尖铰接组合式飞行平台飞行力学建模方法 |
CN112364571A (zh) * | 2020-10-09 | 2021-02-12 | 天津大学 | 大型复杂耦合航天器动力学模型建模方法 |
CN113111430A (zh) * | 2021-03-06 | 2021-07-13 | 北京航空航天大学 | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 |
CN113158528A (zh) * | 2021-05-14 | 2021-07-23 | 南京航空航天大学 | 一种空间充气展开结构的动力学建模方法及系统 |
CN113191099A (zh) * | 2021-04-28 | 2021-07-30 | 北京航空航天大学 | 一种考虑结冰影响的无人机动力学建模方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8380473B2 (en) * | 2009-06-13 | 2013-02-19 | Eric T. Falangas | Method of modeling dynamic characteristics of a flight vehicle |
CN108972558B (zh) * | 2018-08-16 | 2020-02-21 | 居鹤华 | 一种基于轴不变量的多轴机器人动力学建模方法 |
-
2021
- 2021-11-15 CN CN202111345980.5A patent/CN113919081B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106484969A (zh) * | 2016-09-23 | 2017-03-08 | 中国运载火箭技术研究院 | 一种大包线强机动飞行器动力学高精度仿真方法 |
CN108427322A (zh) * | 2018-03-16 | 2018-08-21 | 南京航空航天大学 | 一种大柔性飞行器基于在线辨识的建模方法 |
CN108959796A (zh) * | 2018-07-18 | 2018-12-07 | 哈尔滨工业大学 | 一种大惯量旋转载荷卫星的刚柔磁耦合动力学建模方法 |
CN111931292A (zh) * | 2020-08-07 | 2020-11-13 | 北京航空航天大学 | 一种翼尖铰接组合式飞行平台飞行力学建模方法 |
CN112364571A (zh) * | 2020-10-09 | 2021-02-12 | 天津大学 | 大型复杂耦合航天器动力学模型建模方法 |
CN113111430A (zh) * | 2021-03-06 | 2021-07-13 | 北京航空航天大学 | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 |
CN113191099A (zh) * | 2021-04-28 | 2021-07-30 | 北京航空航天大学 | 一种考虑结冰影响的无人机动力学建模方法 |
CN113158528A (zh) * | 2021-05-14 | 2021-07-23 | 南京航空航天大学 | 一种空间充气展开结构的动力学建模方法及系统 |
Non-Patent Citations (4)
Title |
---|
Liang Xu等.Fight dynamics modeling and analysis of high aspect ratio flexible aircraft.《IEEE Xplore》.2018,全文. * |
于扬等.柔性翼飞行器动力学建模与循环求解仿真方法.《北京理工大学学报》.2017,第37卷(第12期),第0-2章. * |
向锦武等.长航时无人机关键技术研究进展.《哈尔滨工业大学学报》.2020,第52卷(第6期),全文. * |
郭东 ; 徐敏 ; 陈士橹 ; 钱宇 ; .准坐标系下的弹性飞行器飞行动力学建模.系统仿真学报.2010,(第11期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113919081A (zh) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109460596B (zh) | 一种飞翼无人机非线性载荷计算方法 | |
CN113919081B (zh) | 一种考虑惯性耦合的柔性飞行动力学建模与分析方法 | |
BUTTRILL et al. | Nonlinear simulation of a flexible aircraft in maneuvering flight | |
Ashraf et al. | Dynamic modeling of the airship with Matlab using geometrical aerodynamic parameters | |
Zhao et al. | Multibody dynamic approach of flight dynamics and nonlinear aeroelasticity of flexible aircraft | |
CN113868771B (zh) | 一种考虑结构和气动非线性的飞行动力学建模方法 | |
Otsuka et al. | Aeroelastic deployable wing simulation considering rotation hinge joint based on flexible multibody dynamics | |
Selitrennik et al. | Computational aeroelastic simulation of rapidly morphing air vehicles | |
Scarlett et al. | Multibody dynamic aeroelastic simulation of a folding wing aircraft | |
Wang et al. | Modeling and simulation of mass-actuated flexible aircraft for roll control | |
Su et al. | Nonlinear aeroelasticity of highly flexible joined-wing aircraft using unsteady vortex-lattice method | |
Liu et al. | Continuous dynamic simulation for morphing wing aeroelasticity | |
Ting et al. | Development of an Integrated Nonlinear Aeroservoelastic Flight Dynamic Model of the Truss-Braced Wing Aircraft | |
Natella et al. | The effects of aeroelastic tailoring on flight dynamic stability | |
CN115729264A (zh) | 一种基于柔性自适应翼梢小翼的变稳隐身飞机控制方法 | |
Riso et al. | Coupled flight dynamics and aeroelasticity of very flexible aircraft based on commercial finite element solvers | |
Mardanpour et al. | Effect of engine placement on aeroelastic trim and stability of flying wing aircraft | |
Zhao et al. | Modeling and analysis of landing collision dynamics for a shipborne helicopter | |
Noll et al. | Active flexible wing program | |
Otsuka et al. | Parametric studies for the aeroelastic analysis of multibody wings | |
Kier et al. | Development of aircraft flight loads analysis models with uncertainties for pre-design studies | |
Qiu et al. | Improving aileron effectiveness based on changing the position of aileron connectors | |
Nguyen et al. | Nonlinear Aeroelasticity of Flexible Wing Structure Coupled with Aircraft Flight Dynamics | |
Tsushima et al. | Structural and aerodynamic models for aeroelastic analysis of corrugated morphing wings | |
Sabri et al. | Aeroelastic behaviour of a flexible morphing wing design for unmanned aerial vehicle |
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 |