[发明专利]一种直升机旋翼/尾桨气动干扰数值仿真方法有效
申请号: | 201510078995.8 | 申请日: | 2015-02-13 |
公开(公告)号: | CN104881510B | 公开(公告)日: | 2018-06-05 |
发明(设计)人: | 徐国华;史勇杰;叶舟;樊枫 | 申请(专利权)人: | 南京航空航天大学 |
主分类号: | G06F17/50 | 分类号: | G06F17/50 |
代理公司: | 江苏圣典律师事务所 32237 | 代理人: | 贺翔 |
地址: | 210016 江*** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开的直升机旋翼/尾桨气动干扰数值仿真方法,根据直升机相关参数建立了一种准确的旋翼/尾桨气动干扰的数值仿真模型,通过该数值仿真模型在各种状态下对旋翼/尾桨干扰流场进行CFD计算,并对比分析计算结果。本发明的数值仿真方法建立的旋翼/尾桨气动干扰的数值仿真模型在结果上更接近于实际的干扰结果。将不同的旋翼/尾桨参数输入到计算输入文件中,并进行计算,得出的干扰流场结果,就可以模拟直升机旋翼/尾桨在不同飞行状态的干扰情况,并避免了传统涡系方法中对旋翼/尾桨气动干扰进行计算而带来的误差,可直接服务于直升机设计、飞行性能计算以及飞行品质分析。 | ||
搜索关键词: | 尾桨 气动干扰 旋翼 数值仿真模型 直升机旋翼 数值仿真 干扰流场 直升机设计 参数输入 对比分析 飞行品质 飞行性能 飞行状态 输入文件 相关参数 涡系 直升机 分析 服务 | ||
【主权项】:
一种直升机旋翼/尾桨气动干扰数值仿真方法,其特征在于:包括以下步骤:1)根据常规直升机的构型参数,确定旋翼和尾桨的展弦比以及相对位置,从而建立用于旋翼/尾桨气动干扰数值模拟的计算模型;2)生成相应的旋翼和尾桨网格;3)分别在悬停状态、低速前飞状态、中速前飞状态和高速前飞状态下,针对底向前和底向后尾桨、高位和低位尾桨,对旋翼/尾桨干扰流场进行CFD计算;所述CFD计算包括以下步骤:301)读入直升机计算参数、旋翼和尾桨桨叶贴体网格以及背景网格;302)对直升机旋翼和尾桨桨叶网格以及背景网格进行前处理;303)对旋翼和尾桨网格在背景网格中进行挖洞和贡献单元搜索,并将两者的嵌套关系进行储存;304)进行流场初始化;305)更新直升机旋翼和尾桨桨叶网格坐标;306)更新直升机旋翼和尾桨桨叶网格和背景网格之间的嵌套关系;307)进行Lynx直升机旋翼和尾桨桨叶网格上的流场计算,包括无粘通量、湍流模型、粘性通量计算以及LU‑SGS隐式时间推进,并更新背景网格上对应的洞边界单元的流场信息;在旋翼/尾桨干扰流场计算中采用N‑S方程作为主控方程:式中,W为守恒变量,F(W)和G(W)分别表示无粘通量和粘性通量;Vol表示任意单元体,dVol表示单元体Vol的体积,表示Vol的边界,dS表示上的面元;使用MUSCL格式对单元内流场变量进行重构,并采用低耗散的Roe格式进行空间离散,根据Roe格式,网格面上的无粘通量计算公式为: 式中,(·)ij表示网格面ij上的值,R和L分别表示网格面的左右两侧;选用Venkatakrishnan限制器:其中,Δl,min=Wmin‑Wi,Δl,max=Wmax‑Wi,而Wmin和Wmax为单元体i及其所有相邻单元中最小和最大流体变量值;l2为与网格尺度有关的控制变量;Δ2定义为单元体i的流体变量沿格心到面心矢量的方向梯度;采用双时间法模拟干扰流场的非定常效应,引入伪时间,并对伪时间和物理时间导数分别进行差分处理,可得到: 其中,m为伪时间步,n为物理时间步;Δτ为伪时间增量,Δt为物理时间增量;在双时间法中的伪时间步引入LU‑SGS隐式时间;无矩阵形式的LU‑SGS格式,具体计算方法可分为两个步骤:3071)向前扫略: 3072)向后扫略:
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于南京航空航天大学,未经南京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201510078995.8/,转载请声明来源钻瓜专利网。
- 上一篇:自校准无源互调测试仪
- 下一篇:暖奶器用外凸式电加热锅