一种考虑惯性耦合的柔性飞行动力学建模与分析方法与流程

1.本发明属于飞行器设计技术领域，具体涉及一种考虑惯性耦合的柔性飞行动力学建模与分析方法，可以用于采用细长体布局的高速飞行器、大型运输机、超大展弦比无人机和微小型等柔性飞行器飞行动力学建模和动力学分析及后续的控制设计。


背景技术：

2.随着空天科学技术的不断进步，飞行器速度、高度、航时的记录不断突破，创新气动布局不断涌新，轻质材料广泛使用，对飞行器设计工程实践和学科发展提出新的要求。采用机身刚体假设的传统飞行力学建模无法满足结构柔性飞行动力学建模与控制设计问题的需要。柔性飞行器飞行动力学建模是实现动力学分析和控制设计的基础。对于柔性飞行器，飞行动力学模型需要反映飞行器大范围运动，又能反映结构的动力学行为，同时还需要能够充分利用已有的分析方式和可获取的气动、结构等参数的数据源。


技术实现要素：

3.本发明针对柔性飞行器动力学建模问题，利用准坐标系下的拉格朗日方程建立了柔性飞行器飞行动力学模型，其中能够考虑惯性耦合，方便整理成用于分析和仿真的形式，分析采用平均轴系假设和忽略惯性耦合后，建模方法之间的退化形式，以及刚体运动和弹性运动之间的耦合机理。具体的技术方案为：
4.一种考虑惯性耦合的柔性飞行动力学建模与分析方法，包括以下步骤：
5.(1)采用以下假设：不考虑地球曲率和自转；飞机结构描述为一系列的集中质量，相应的集中质量和惯性矩为mi和ji；变形满足线弹性理论和小变形假设；为了方便结构模型处理，认为能够获得满足正交化假设的结构模态；
6.(2)定义以下坐标系：1.惯性系(xe,ye,ze)：与大地固连，根据平面大地假设可认为是惯性系，坐标原点oe；2.体轴系(xb,yb,zb)：未变形状态的机体坐标系，原点ob取在飞行器未变形时的质心上，xb指向头部为正，yb垂直于xb，方向指向机身右侧，zb按右手定则确定；3.当地坐标系(xi,yi,zi)：在未变形机体上的每个节点，相对体轴系原点的位置向量为ri，坐标轴与体轴系平行。
7.(3)建立飞行器运动学关系方程。考虑未变形状态，与刚体飞行器固联的体轴系在惯性系的位置向量r
0e
＝(x
0e
,y
0e
,z
0e
)
t
，用欧拉角φ,θ,ψ表达从惯性坐标系到体轴系的姿态变换θ＝(φ,θ,ψ)
t
，速度和角速度在体轴系的表达为vb＝(u,v,w)
t
，ωb＝(p,q,r)
t
，描述飞行器质点运动和姿态运动的运动学关系方程为
[0008][0009][0010]
其中，φ、θ、ψ分别为滚转角、俯仰角和偏航角；p、q、r分别为滚转角速度、俯仰角速度和偏航角速度；u、v、w分别为飞行器速度在体轴系三个轴上的投影；
[0011][0012][0013]
体轴系原点位置对于每个变形后集中质量位置在体轴系中的表达为
[0014][0015]
其中r0为惯性系中体轴系坐标原点的位置向量在体轴系的表达，di为当地坐标系中的变形，si为集中质量在当地坐标系中的位置向量，为扭转变形。t表示变换矩阵，描述位置向量si从未变形位置到变形后位置的扭转。在线性假设下，对其泰勒展开取一阶
[0016][0017]
其中表示的反对称矩阵。
[0018]
集中质量位移ri：
[0019][0020]
集中质量速度
[0021][0022]
集中质量的转动速度ωi由体轴系的转动ωb和变形带来的转动共同组成：
[0023][0024]
节点的弹性位移可以表达为平移和转动，根据假设条件，可以满足正交条件的自由-自由边界条件下的模态表达，设弹性变形的广义坐标为ηe，节点弹性位移可以表达为：
[0025][0026]
其中，di为节点的运动变形，为节点的扭转变形，和分别为广义坐标ηe中平动和转动自由度对应的广义坐标模态；
[0027]
(4)建立飞行器动力学关系方程。采用准坐标系下表达的拉格朗日方程，建立动力学方程：
[0028]
包含平动自由度方程：
[0029][0030]
其中pg＝[f m]
t
为节点上的载荷，包含非保守外力。非保守力包含气动力和其他外力，其中气动力可以采用准定常方式获得气动力表达。为广义模态坐标下的广义模
态，为
[0031]
转动自由度方程：
[0032][0033]
在计算分析中，h,j,每一步都需要重新计算。对于变形后的飞行器，惯性张量展开为
[0034][0035]
将表达式中的每一项编号表示为：
[0036][0037]
a1不随时间变化，来自坐标轴偏移。
[0038]
有
[0039][0040][0041][0042][0043][0044][0045][0046]
e为单位向量，e1＝[1 0 0]
t
，e2＝[0 1 0]
t
，e3＝[0 0 1]
t
；
[0047]
有
[0048]
[0049][0050][0051][0052][0053][0054][0055]
其中，
[0056][0057]ai1
＝s
iz
φ
giyer-s
iy
φ
gizer
[0058]ai2
＝s
ix
φ
gizer-s
iz
φ
gixer
[0059]ai3
＝s
iy
φ
gixer-s
ix
φ
giyer
[0060]
有
[0061][0062][0063][0064][0065][0066]
a5＝a4
t
[0067]
有
[0068][0069]
[0070][0071][0072][0073]
有
[0074][0075][0076][0077][0078][0079]
a8＝a6
t
[0080]
a9＝a7
t
[0081]
定义
[0082]
惯性张量j的表达式整理为：
[0083][0084]
j对时间求导：
[0085][0086]
类似的，动量矩h可以用以下公式进行计算：
[0087][0088]
表达式中的各项编号，可以写为：有
[0089][0090][0091]
[0092][0093][0094]
其中，
[0095][0096][0097][0098][0099][0100][0101][0102][0103][0104][0105][0106][0107][0108][0109][0110][0111]
特别的，为方便模型的求解，需要对系统拉格朗日方程的原始表达形式进行解释：
[0112]
[0113][0114][0115]
其中l为拉格朗日算子，l＝e
k-e
p
，ek、e
p
分别为动能和势能，q
t
，qr，qe为非保守外力，满足：虚功
[0116]
动能利用下式进行计算：
[0117][0118]
对于动能项ek求解，需要计算每个集中质量的平动动能e
kt
和转动动能e
kr
：
[0119][0120]
其中，ri为集中质量的位置矢量，ji为集中质量的转动惯量；
[0121]
结合实用平均轴系假设，求得平动动能：
[0122][0123]
其中平动动能公式最后两项体现了刚性运动与弹性变形之间的耦合。
[0124]
转动动能：
[0125][0126]
令总动能可以求解、整理为
[0127][0128]
利用广义坐标和模态表示变形，第二项和第五项整理为
[0129][0130]
因此，总动能表达式为
[0131][0132]
其中，ri为集中质量的位置矢量，ji为集中质量的转动惯量；
[0133]
对于势能项e
p
计算，势能包含重力势能和弹性势能，利用线弹性假设，全机势能表达式为
[0134][0135]
其中，t
eb
＝t
be-1
，ge＝[0 0 g]，g为当地重力加速度，k
gg
为刚度阵，ug为弹性变形；
[0136]
对于非保守外力求解，主要计算气动力、推力和其他外力，其中气动力表达为剖面升力其中为动压，αs为剖面迎角，机翼升力由剖面升力沿展向积分得到，即为方便计算气动力会转化为升力系数的形式，为该领域基本问题，在此不赘述。
[0137]
考虑每个节点上的载荷p
g，i
＝[f
i mi]
t
，非保守力和力矩的虚功用下式计算：
[0138][0139]
其中δa为体轴系的虚角位移。也可以利用每个节点的虚位移和刚体运动模态表达
[0140][0141]
其中广义力：
[0142][0143][0144][0145]
建立飞行器弹性关系方程。通过拉格朗日变量对广义变形坐标求导，并在模型中加入对称形式的阻尼弹性关系方程可以整理为：
[0146][0147]
其中m
ee
为广义质量矩阵，k
ee
为广义刚度矩阵。
[0148]
本方法得到的飞行器弹性关系方程为：
[0149][0150]
其中，ηe为弹性变形的广义坐标，ej为单位向量，来源于转动惯量j的表达式各项，有
[0151][0152][0153][0154][0155][0156][0157][0158][0159][0160][0161][0162][0163][0164][0165][0166][0167][0168]
其中，
[0169]
来源于动量矩表达式中的各项，有：
[0170][0171][0172][0173]
(5)利用数值方法对方程其进行求解，计算本发明所述建模方式在初始条件下的时域响应。本方法建立考虑惯性耦合的柔性飞行器飞行动力学方程总结如下：
[0174]
运动学关系方程：
[0175][0176][0177]
动力学关系方程：
[0178][0179][0180][0181][0182]
进一步，通过数值方法计算模型在特定初始条件作用下运动参数的时域响应，分析飞行器模型的稳定性。
[0183]
本方法利用准坐标系描述的拉格朗日方程推导了柔性飞行器飞行动力学方程，模型采用的假设条件和结构模型便于利用已有的气动参数和结构有限元模型，方程中包含惯性耦合项，如果采用更为苛刻的简化条件，能够获得模型的进一步简化，便于不同模型精度下的分析，适应不同的分析和控制设计需求。如不考虑刚体运动和弹性运动惯性耦合，可以对动力学模型进行简化为常规动力学方程形式：
[0184]
[0185][0186]
附图说明
[0187]
图1为飞行动力学建模流程示意图；
[0188]
图2为实施例结构及坐标系示意图；
[0189]
图3为实施例节点1在体轴系内的位移时域响应；
[0190]
图4为实施例在节点1在体轴系内的位移的时域响应；
[0191]
图5为实施例全机角速度时域响应曲线；
[0192]
图6为实施例节点1转角的时域响应。
具体实施方式
[0193]
下面结合附图和实施例对本发明的考虑惯性耦合的柔性飞行器的飞行动力学建模与分析方法进一步详细说明。
[0194]
具体步骤如图1所述，本实施例针对虑极端简化的飞机模型(图2)，将机身和发动机考虑为欧拉梁连接中部和两端的集中质量，对其进行柔性飞行动力学建模。
[0195]
第一步，采用以下假设：不考虑地球曲率和自转；飞机结构描述为一系列的集中质量，相应的集中质量和惯性矩为mi和ji；变形满足线弹性理论和小变形假设；为了方便结构模型处理，认为能够获得满足正交化假设的结构模态；
[0196]
第二步，对飞行器进行坐标系定义，体轴系原点位于中心处，初始位置与惯性系重合，确定飞行器简化模型模型各参数为
[0197]
m1＝1kg,m2＝2kg,m3＝1kg,j1＝8
×
10-4
kg.m2,j2＝2.5
×
10-3
kg.m2,j3＝8
×
10-4
kg.m2[0198]
a＝1
×
10-5
m2,e＝2
×
109n/m2,l0＝1m,i＝1
×
10-8
m-4
[0199]
其中，m1和m3为吊舱质量，j1和j3为吊舱转动惯量，m2为机身质量，j2为机身转动惯量，a为机翼截面积，e为机翼刚度，i为机翼截面惯性矩。
[0200]
第三步，建立模型的运动学关系方程。当前模型在仅考虑横向平面运动前提下，运动学关系方程简化为
[0201]
第四步，建立模型的动力学关系方程，弹性变形以向量描述为
[0202][0203]
将位移用广义坐标形式表达为μg＝φ
ge
ηe[0204]
其中φ
ge
为欧拉梁模型在自由-自由边界条件下获得的振型。
[0205]
分析得到模型动力学关系方程如下：
[0206]
[0207]
其中，分析外力矩m
x
在短时间内的作用，设m
x
为
[0208][0209]
建立模型的弹性关系方程，模型中加入对称形式的阻尼，即
[0210]
飞行器的弹性关系方程可以表达为：
[0211][0212]
第五步，通过积分计算本发明所述建模方式在给定力矩作用下运动参数的时域响应。如图3为节点1的位置在体轴系的变化，图4为节点1在y轴的位移随时间的变化。结果显示，在体轴系内观察节点1的位移，可以发现，惯性耦合后的结果在滚转力矩的作用下存在y方向的位移，即伸长。图5所示为全机角速度的时域响应，图6所示为节点1的时域响应。
[0213]
可见对于研究飞机刚体运动为主的飞行动力学与控制问题，飞行器飞行过程中的载荷分析和计算以及相关的气动弹性问题时，考虑惯性耦合可以提供精度更高的分析结果尤其是存在较为明显的集中质量和做类似高速滚转等特殊机动的情况。