一种快速计算旋翼刚体挥舞稳态周期响应的方法及系统技术方案

本发明专利技术涉及一种快速计算旋翼刚体挥舞稳态周期响应的方法及系统，属于旋翼挥舞响应计算技术领域，利用旋翼区域网格迭代计算刚体挥舞角、挥舞速度和流场速度等，再根据叶素理论计算旋翼的气动力；利用气动力计算旋翼动量源项；将旋翼动量源项添加至纳维

【技术实现步骤摘要】
一种快速计算旋翼刚体挥舞稳态周期响应的方法及系统


[0001]本专利技术涉及旋翼挥舞响应计算
，特别是涉及一种快速计算旋翼刚体挥舞稳态周期响应的方法及系统。

技术介绍

[0002]对于常规直升机、共轴双旋翼直升机、自转旋翼机、倾转旋翼机等旋翼飞行器，其旋翼存在周期性的挥舞响应。在前飞等旋翼气动力不平衡的状态下，气动力大的一侧旋翼将受迫向上挥舞从而使迎角减小，气动力也随之减小，气动力小的一侧旋翼将受迫向下挥舞从而增大迎角，气动力也随之增大，旋翼的挥舞响应可大幅减少旋翼飞行器的倾翻力矩，对其稳定性和操纵性作用明显，同时旋翼的挥舞响应特性也对旋翼飞行器的气动特性与性能、动力学特性以及飞行性能等有重要影响，因此，在设计与应用阶段如何快速、准确地计算和评估旋翼的挥舞响应特性是开展气动、动力学、飞行力学等特性分析的前提。旋翼挥舞响应计算属于流固耦合仿真领域，流体和结构双向耦合的特性导致旋翼挥舞响应计算方法复杂、计算量极大且收敛慢，如何提高计算效率、减小计算时长一直是旋翼挥舞响应计算的热点和难点问题，然而现有的快速计算方法不能全面的考虑除旋翼外的其他气动部件对流场的影响。

技术实现思路

[0003]本专利技术的目的是提供一种快速计算旋翼刚体挥舞稳态周期响应的方法及系统，将旋翼动量源项和旋翼刚体挥舞响应松耦合以快速计算旋翼稳态刚体挥舞响应，在全面考虑气动部件对流场影响的前提下有效提高了计算效率。
[0004]为实现上述目的，本专利技术提供了如下方案：
[0005]一种快速计算旋翼刚体挥舞稳态周期响应的方法，所述方法包括：
[0006]获取包含旋翼的一个任意空间作为旋翼CFD流场计算域；
[0007]对所述旋翼CFD流场计算域进行网格划分，得到旋翼区域网格；
[0008]在桨毂坐标系下利用旋翼位于各方位角处的刚体尖端挥舞位移X
n
计算旋翼位于各方位角处的刚体挥舞角β
flapn
；
[0009]在桨毂坐标系下利用旋翼的尖端挥舞速度V
n
计算旋翼的各叶素微段的挥舞速度V
flapn
；所述叶素微段为旋翼沿展向的一截；
[0010]利用惯性坐标系下所述旋翼区域网格处的x方向流场速度V
Gxn
、y方向流场速度V
Gyn
和z方向流场速度V
Gzn
获取桨毂坐标系下旋翼的各叶素微段位于各方位角处的第一法向流场速度Vn
planeofdiskn
、第一径向流场速度Vr
planeofdiskn
和第一切向流场速度Vt
planeofdiskn
；所述x方向流场速度、y方向流场速度和z方向流场速度分别为惯性坐标系的x方向的流场速度、y方向的流场速度和z方向的流场速度；所述第一法向流场速度、第一径向流场速度和第一切向流场速度分别为桨毂坐标系的法向方向的流场速度、径向方向的流场速度和切向方向的流场速度；
[0011]利用桨毂坐标系下的所述旋翼刚体挥舞角β
flapn
、所述旋翼的各叶素微段位于各方位角处的第一法向流场速度Vn
planeofdiskn
和第一径向流场速度Vr
planeofdiskn
，计算叶素微段坐标系下旋翼的各叶素微段位于各方位角处的第二法向流场速度Vn
bladen
和第二切向流场速度Vt
bladen
；所述第二法向流场速度和第二切向流场速度分别为叶素微段坐标系的法向方向的流场速度和切向方向的流场速度；
[0012]在叶素微段坐标系下利用旋翼的各叶素微段的挥舞速度V
flapn
和旋翼的各叶素微段位于各方位角处的第二法向流场速度Vn
bladen
计算流场相对于叶素微段的法向速度Vn
bladebladen
；
[0013]获取旋翼总距操纵输入φ0、旋翼纵向周期变距操纵输入A以及旋翼横向周期变距操纵输入B，并计算旋翼桨距角；
[0014]在叶素微段坐标系下，利用旋翼的各叶素微段位于各方位角处的第二切向流场速度Vt
bladen
、流场相对于叶素微段的法向速度Vn
bladebladen
和所述旋翼桨距角，根据叶素理论计算旋翼的各叶素微段位于各方位角处的气动力Fareo
n
；
[0015]利用所述气动力Fareo
n
计算旋翼位于各方位角处的旋翼动量源项Fsource
n
；
[0016]将所述旋翼动量源项Fsource
n
添加至纳维
‑
斯托克斯方程的动量源项中，并根据所述纳维
‑
斯托克斯方程利用计算流体力学的方法更新所述旋翼区域网格处在惯性坐标系下的x方向流场速度、y方向流场速度和z方向流场速度；
[0017]利用所述气动力Fareo
n
以松耦合的方式更新旋翼刚体挥舞响应的尖端挥舞加速度尖端挥舞速度和尖端挥舞位移并返回步骤“在桨毂坐标系下利用旋翼位于各方位角处的刚体尖端挥舞位移X
n
计算旋翼位于各方位角处的刚体挥舞角β
flapn”；
[0018]直到相邻两次迭代的尖端挥舞位移的差值小于预设误差值，停止迭代。
[0019]本专利技术还提供了一种快速计算旋翼刚体挥舞稳态周期响应的系统，所述系统包括：
[0020]计算域选取单元，用于获取包含旋翼的一个任意空间作为旋翼CFD流场计算域；
[0021]网格划分单元，用于对所述旋翼CFD流场计算域进行网格划分，得到旋翼区域网格；
[0022]刚体挥舞角计算单元，用于在桨毂坐标系下利用旋翼位于各方位角处的刚体尖端挥舞位移X
n
计算旋翼位于各方位角处的刚体挥舞角β
flapn
；
[0023]挥舞速度计算单元，用于在桨毂坐标系下利用旋翼的尖端挥舞速度V
n
计算旋翼的各叶素微段的挥舞速度V
flapn
；所述叶素微段为旋翼沿展向的一截；
[0024]流场速度获取单元，用于利用惯性坐标系下所述旋翼区域网格处的x方向流场速度V
Gxn
、y方向流场速度V
Gyn
和z方向流场速度V
Gzn
获取桨毂坐标系下旋翼的各叶素微段位于各方位角处的第一法向流场速度Vn
planeofdiskn
、第一径向流场速度Vr
planeofdiskn
和第一切向流场速度Vt
planeofdiskn
；所述x方向流场速度、y方向流场速度和z方向流场速度分别为惯性坐标系的x方向的流场速度、y方向的流场速度和z方向的流场速度；所述第一法向流场速度、第一径向流场速度和第一切向流场速度分别为桨毂坐标系的法向方向的流场速度、径向方向的流场速度和切向方向的流场本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种快速计算旋翼刚体挥舞稳态周期响应的方法，其特征在于，所述方法包括：获取包含旋翼的一个任意空间作为旋翼CFD流场计算域；对所述旋翼CFD流场计算域进行网格划分，得到旋翼区域网格；在桨毂坐标系下利用旋翼位于各方位角处的刚体尖端挥舞位移X
n
计算旋翼位于各方位角处的刚体挥舞角β
flapn
；在桨毂坐标系下利用旋翼的尖端挥舞速度V
n
计算旋翼的各叶素微段的挥舞速度V
flapn
；所述叶素微段为旋翼沿展向的一截；利用惯性坐标系下所述旋翼区域网格处的x方向流场速度V
Gxn
、y方向流场速度V
Gyn
和z方向流场速度V
Gzn
获取桨毂坐标系下旋翼的各叶素微段位于各方位角处的第一法向流场速度Vn
planeofdiskn
、第一径向流场速度Vr
planeofdiskn
和第一切向流场速度Vt
planeofdiskn
；所述x方向流场速度、y方向流场速度和y方向流场速度分别为惯性坐标系的x方向的流场速度、y方向的流场速度和z方向的流场速度；所述第一法向流场速度、第一径向流场速度和第一切向流场速度分别为桨毂坐标系的法向方向的流场速度、径向方向的流场速度和切向方向的流场速度；利用桨毂坐标系下的所述旋翼刚体挥舞角β
flapn
、所述旋翼的各叶素微段位于各方位角处的第一法向流场速度Vn
planeofdiskn
和第一径向流场速度Vr
planeofdiskn
，计算叶素微段坐标系下旋翼的各叶素微段位于各方位角处的第二法向流场速度Vn
bladen
和第二切向流场速度Vt
bladen
；所述第二法向流场速度和第二切向流场速度分别为叶素微段坐标系的法向方向的流场速度和切向方向的流场速度；在叶素微段坐标系下利用旋翼的各叶素微段的挥舞速度V
flapn
和旋翼的各叶素微段位于各方位角处的第二法向流场速度Vn
bladen
计算流场相对于叶素微段的法向速度Vn
bladebladen
；获取旋翼总距操纵输入φ0、旋翼纵向周期变距操纵输入A以及旋翼横向周期变距操纵输入B，并计算旋翼桨距角；在叶素微段坐标系下，利用旋翼的各叶素微段位于各方位角处的第二切向流场速度Vt
bladen
、流场相对于叶素微段的法向速度Vn
bladebladen
和所述旋翼桨距角，根据叶素理论计算旋翼的各叶素微段位于各方位角处的气动力Fareo
n
；利用所述气动力Fareo
n
计算旋翼位于各方位角处的旋翼动量源项Fsource
n
；将所述旋翼动量源项Fsource
n
添加至纳维
‑
斯托克斯方程的动量源项中，并根据所述纳维
‑
斯托克斯方程利用计算流体力学的方法更新所述旋翼区域网格处在惯性坐标系下的x方向流场速度、y方向流场速度和z方向流场速度；利用所述气动力Fareo
n
以松耦合的方式更新旋翼刚体挥舞响应的尖端挥舞加速度尖端挥舞速度V
in+1
和尖端挥舞位移并返回步骤“在桨毂坐标系下利用旋翼位于各方位角处的刚体尖端挥舞位移X
n
计算旋翼位于各方位角处的刚体挥舞角β
flapn”；直到相邻两次迭代的尖端挥舞位移的差值小于预设误差值，停止迭代。2.根据权利要求1所述的快速计算旋翼刚体挥舞稳态周期响应的方法，其特征在于，在所述对所述旋翼CFD流场计算域进行网格划分，得到旋翼区域网格之后，还包括：对所述旋翼区域网格中每个网格处的x方向流场速度、y方向流场速度、z方向流场速度
和流场密度赋初值；对旋翼位于各方位角处的刚体尖端挥舞位移和尖端挥舞速度赋初值。3.根据权利要求1所述的快速计算旋翼刚体挥舞稳态周期响应的方法，其特征在于，所述对所述旋翼CFD流场计算域进行网格划分，得到旋翼区域网格，具体包括：对所述旋翼CFD流场计算域进行网格划分，得到网格计算域；在所述网格计算域内使用旋翼动量源方法模拟旋翼，得到旋翼区域网格。4.根据权利要求1所述的快速计算旋翼刚体挥舞稳态周期响应的方法，其特征在于，所述在桨毂坐标系下利用旋翼位于各方位角处的刚体尖端挥舞位移X
n
计算旋翼位于各方位角处的刚体挥舞角β
flapn
的计算公式为：β
flapn
＝X
n
/R其中，β
flapn
表示第n次迭代中桨毂坐标系下旋翼位于各方位角处的刚体挥舞角，X
n
表示第n次迭代中桨毂坐标系下旋翼位于各方位角处的刚体尖端挥舞位移，R表示旋翼半径。5.根据权利要求1所述的快速计算旋翼刚体挥舞稳态周期响应的方法，其特征在于，所述在桨毂坐标系下利用旋翼的尖端挥舞速度V
n
计算旋翼的各叶素微段的挥舞速度V
flapn
的计算公式为：V
flapn
＝V
n
×
r/R其中，V
flapn
表示第n次迭代中桨毂坐标系下旋翼的各叶素微段的挥舞速度，V
n
表示第n次迭代中的桨毂坐标系下的尖端挥舞速度，r表示叶素微段坐标系中心所在的半径，R表示旋翼半径。6.根据权利要求1所述的快速计算旋翼刚体挥舞稳态周期响应的方法，其特征在于，桨毂坐标系下旋翼的各叶素微段位于各方位角处的第一法向流场速度Vn
planeofdiskn
的计算公式为：Vn
planeofdiskn
＝V
Gzn
；桨毂坐标系下旋翼的各叶素微段位于各方位角处的第一径向流场速度Vr
planeofdiskn
的计算公式为：Vr
planeofdiskn
＝V
Gxn
×
sin(θ
n
)+V
Gyn
×
cos(θ
n
)；其中，θ
n
为第n迭代步旋翼方位角；桨毂坐标系下旋翼的各叶素微段位于各方位角处的第一切向流场速度Vt
planeofdiskn
的计算公式为：Vt
planeofdiskn
＝V
Gxn
×
cos(θ
n
)
‑
V
Gyn
×
sin(θ
n
)。7.根据权利要求1所述的快速计算旋翼刚体挥舞稳态周期响应的方法，其特征在于，所述叶素微段坐标系下旋翼...

【专利技术属性】
技术研发人员：招启军，王宏亮，赵国庆，王博，陈希，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：