更全的杂志信息网

基于CFD方法的主动襟翼控制旋翼翼型涡特性研究

更新时间:2009-03-28

降低直升机振动水平及噪声大小是直升机优化设计中关注的主要问题。近年来,旋翼主动襟翼控制(Active flap control,AFC) [1-2]技术凭借良好的减振[1]作用及降噪[3-4]特性逐渐成为直升机界的研究热点。

AFC技术最早由Lemnios和Smith在使用伺服襟翼降低旋翼振动载荷[5]的研究中被提出。马里兰大学的Milgram在文献[2]中将其发展应用于降低旋翼气动载荷。而后,Schmitz等人将其用于旋翼厚度噪声的降噪研究[3],Aoyama又将其发展应用于降低桨-涡干扰噪声[4]。在国内,夏品奇等也对AFC旋翼的载荷控制[6]与气动特性[7]进行了相关的研究。

图1给出了AFC旋翼桨叶的示意图,与常规直升机旋翼桨叶相比,AFC旋翼桨叶不仅具有复杂的外形和结构,且其后缘小翼在伴随桨叶运动的同时还具有相对独立的偏转运动,使得流场非定常变化,难以精确模拟。另外,在网格生成方面,要求能准确处理主桨叶和后缘小翼间缝隙处的网格。此外,考虑到旋翼本身的挥舞、摆振及变距,这将使得后缘小翼的运动相当复杂,运动嵌套的实现也变得更为棘手。其复杂性使得目前针对AFC旋翼气动特性的研究并不多见,而详细探究AFC旋翼翼型参数对其尾涡产生的影响,并能用于指导AFC旋翼设计的工作则更为少见。

  

图1 AFC旋翼桨叶示意图Fig.1 Schematic of AFC rotor blade

针对上述问题,本文采用基于Euler方程与Jameson中心格式的数值模拟求解方法,编写了可用于非结构多重嵌套条件下的嵌套网格处理程序,并结合基于坐标系变换的运动描述方式,提出可用于AFC旋翼数值模拟的运动嵌套网格方法。然后建立了适于AFC旋翼翼型复杂运动条件下的流场数值模拟方法。应用所建立的方法,深入地开展了AFC旋翼翼型参数对翼型涡量影响的计算研究,并得出了一些新的结论。

技术的飞速发展已经开始改变人类历史几千年的文明惯性和生活惯性,且这种改变将加速到来。人工智能的社会效应,既可能极善,也可能极恶。前者是大同或天堂;后者是新奴隶社会或人类毁灭,且不可逆。如何对待这种不确定性?在找到预防和控制其负面效应的手段之前,人类是否应该轻易地发展人工智能?

1 计算方法和模型

1.1 流场控制方程

针对AFC旋翼流场的复杂性,并考虑计算效率,本文以欧拉方程为主控方程,其微分方程形式为

 

(1)

为验证CFD方法对AFC翼型的计算准确性,选取类似的“30P30N”多段翼型[10],图4给出了Ma=0.2、迎角8.5°下“30P30N”翼型后缘襟翼部分对应的压力系数计算结果与文献[10]实验值的对比。由图可见,其吻合程度是合理的,可用于后文的分析。

恒山景区部分宗教活动场所和宗教人员的行为恶劣,有违法设立“功德香”、“烧高香”等敛财行为。虽然国家和地方虽出台了一系列相关的法律、法规、条例,也先后颁布一些修订后的旅游市场的价格法律法规,但旅游价格管理立法仍稍显滞后,使得打击旅游价格违法行为的实时性有所降低,不能有效地约束和制约经营者,导致了景区混乱局面的出现。在景区管理方面,没有做到明确分工,各司其能。春季干燥,是森林火灾的多发季节,应加强监管和宣传,自2006年恒山森林失火惊动央视后,后期虽在努力弥补和监督,但仍收效甚微。

W=

F=

1.1 一般资料 选取自2010—2017年北部战区总医院行原位心脏移植术的34例患者为研究对象。按心脏移植术后是否出现过排斥反应分为排斥组(n=12)及无排斥组(n=22)。排斥组中,男性10例,女性2例;年龄45~63岁,平均年龄(56±6)岁。无排斥组中,男性16例,女性6例;年龄30~66岁,平均年龄(50±10)岁。两组患者一般资料比较,差异无统计学意义(P>0.05),具有可比性。本研究经医院伦理委员会批准,患者均签署知情同意书。

对照组ALT(36.54±11.24)U/L、AST(62.54±12.12)U/L、TBIL(30.24±12.27)umol/L,实验组ALT(30.64±10.65)U/L、AST(58.48±13.14)U/L、TBIL(20.64±12.16)umol/L(p<0.05),见表1。

H=

(2)

式中:ρp分别为流体密度和压强;uvw为速度在3个方向上的分量。

空间离散采用Jameson中心格式,其表达式为

①苏轼《卜算子》(缺月挂疏桐):双调44字,上阕4句2仄韵,下阕4句2仄韵。句式:5575。5575。

 

(3)

1.2 嵌套网格处理方法

本文构建的嵌套网格方法与常规的inverse-map或hole-map方法不同,这里先进行网格贡献单元搜索以确定是否有贡献单元及贡献单元所在位置,再根据网格分类准则参数来决定网格所属分类,所得分类单元的边界即为嵌套网格边界。该嵌套网格方法能处理相互嵌套的网格块,因此非常适合AFC旋翼主桨叶与后缘小翼间的网格嵌套。

关于网格分类准则,考虑到在流场计算中网格尺度对流场求解的精度和收敛效果有很大影响,网格尺度越小,求解精度越高、收敛性越好,故本文使用Löhner[8]的物面距和网格尺度组合参数s作为判据进行分类,s定义为

s=dphq

(4)

式中:d为物面距,h为网格尺度,通过对指数pq的调节可以改变物面距和网格尺度的权重。本文选取p=1,q=1,这样既考虑了物面距影响,又考虑了网格尺度影响。

对于贡献单元搜索,本文采用相邻单元搜索算法与交流数字树(Alter-nating digital tree, ADT)搜索算法[9]相结合的方法。对网格单元i,其在网格块b中贡献单元搜索算法如下:

(1) 取目标坐标点为单元i的形心。

能入选特种部队的人开始都感觉自己特有本事,但不是每名队员都能完全承受训练,陆续有人调离、转到后勤甚至退伍。开始1000人的新兵队伍,几个月后只剩下不到500人的作战人员。

(3) 连线当前单元ic的形心与目标坐标点,判断该线段与单元边的相交情况,如果无相交边,则当前单元为网格单元i的贡献单元,并结束搜索;否则,继续执行步骤(4)。

(4) 判断相交边的类型,如果为中间面,即相邻单元存在,则以该单元为新的出发单元,跳转执行步骤(3);如果为外边界(通常为压力远场),则该单元保留为活动单元,搜索结束;如果该边为内边界(通常为物面边界),则继续执行步骤(5)。

(5) 采用ADT算法,判断步骤(3)中得到的线段与内边界相交情况,如果相交边(或面)个数为奇数,则目标坐标点位于物面内,即单元i在内边界内,并标记为非活动单元;如果为偶数,则取距目标坐标点最近的内边界,并以其相邻单元作为新的出发单元,跳转到步骤(3)并继续。

数值实践表明,上述嵌套网格方法适应范围广、效率高,并具有良好的鲁棒性。

1.3 运动描述方式

旋翼桨叶运动除绕旋转轴旋转外,还具有挥舞、摆振和变距运动。特别地,对包含后缘小翼的AFC旋翼,该小翼还将绕其安装轴运动。上述运动的合成使得描述小翼对地面运动变得十分复杂。

本文提出一种基于坐标系变换的运动描述方式,其基本思路是:通过描述一个坐标系相对另一个坐标系的运动,使得复杂运动能经过多次简单描述后转换为对地面运动的表示。该方法分为建立装配关系以及对运动进行描述两大部分,见图2。

(1) 处理装配关系。对两个已知坐标系O1O2,处理装配关系时先固定参考点,可在空间中任取点b,记其在坐标系O1O2中的坐标为b1b2。然后,固定参考方向,即过参考点b建立新坐标系Ob1Ob2,得坐标系基向量在O1O2中的向量组nb1nb2。注意到点b1b2虽在两坐标系中具有不同坐标值,但实际上为同一位置点。于是,对空间中任意点P,类似的有坐标P1P2,若已知b1nb1b2nb2P2,便可求得P1

 

(5)

(2) 在装配的基础上,增加运动描述项。设时刻t时坐标系Ob1Ob2变为运动后的P1P2成为定义运动平动分量为函数向量T(t),转动分量为函数矩阵R(t)。由于P2相对坐标系Ob2静止,即于是在已知R(t)的条件下,可解得P1 ′为

 

(6)

实践表明,描述两个原本重合的坐标系间的相对运动容易实现。对于AFC后缘小翼,只需依次建立后缘小翼坐标系-主桨叶坐标系-扭转坐标系-摆振坐标系-挥舞坐标系-桨毂坐标系-地面坐标系,便可完整地描述具有挥摆扭运动的AFC旋翼后缘小翼对地面运动方式。

12月18日至19日,浙江省自然资源厅党组书记、厅长黄志平一行先后到绍兴、宁波、台州、温州等地调研。绍兴市、宁波市、台州市、温州市政府分管领导分别陪同调研。

  

图2 坐标系变换示意图Fig.2 Schematic of coordinate transformation

2 算例验证

2.1 运动嵌套网格验证

图3给出了某时刻下AFC后缘小翼附近的嵌套网格,图中后缘部分网格与主旋翼部分相互嵌套,网格嵌套情况符合预期要求。

2.2 “30P30N”多段翼型算例验证

式中:W为守恒变量;FGH为对流通量,其具体表达式如下

  

图3 AFC后缘小翼附近的嵌套网格Fig.3 Overset grid near AFC rotor trailing edge

  

图4 “30P30N”翼型后缘压力系数计算结果Fig.4 Calculation results of pressure coefficient at trailing edge of ″30P30N″ airfoil compared with experimental results[10]

2.3 C-T旋翼算例验证

图5给出了桨尖马赫数0.439下C-T旋翼在0.96R处的翼型表面压力系数分布与实验值[11]的对比,结果吻合良好,说明对于典型旋翼工作条件下本文运动嵌套数值模拟方法能有效地模拟桨叶表面压强分布。

  

图5 C-T旋翼表面压力系数计算结果Fig.5 Calculation results of C-T rotor surface pressure coefficient compared with experimental results[11]

3 AFC旋翼翼型数值模拟

3.1 AFC旋翼翼型数值模拟条件

  

图6 不同时刻下后缘小翼附近涡分布图Fig.6 Vortex near trailing edge at different time

与典型固定翼多段翼型相比,旋翼AFC后缘小翼工作环境存在较大差别。固定翼后缘襟翼多用于起飞降落阶段,来流速度通常在Ma=0.1~0.2之间,且襟翼转动缓慢,可视为定常或准定常流场;但后者位于靠近桨尖部分,其周围流场呈现明显的非定常特性,两者有着显著差别,故需为AFC后缘小翼开展单独的分析工作。本文中AFC旋翼数值模拟的基准条件如表1,2所示。

 

表1 AFC旋翼主桨叶参数

 

Tab.1 Parameters of AFC rotor main blade

  

参数数值翼型NACA0012旋翼半径R/m1弦长c/m0.1桨尖马赫数Ma0.65总距角/(°)6前进比0

 

表2 AFC后缘小翼参数

 

Tab.2 Parameters of AFC trailing edge

  

参数数值弦长0.15c与主桨叶间距0.1%c分布位置0.85R~0.95R

3.2 AFC翼型涡的特性

AFC旋翼对气动特性影响主要通过尾部所产生的涡来实现。图6给出了后缘小翼某一摆动周期内其附近的涡量图。从图中可以观察到,后缘小翼是AFC旋翼涡的主要产生部件,后缘上下表面均可从后缘小翼与主桨叶间缝隙处产生涡。其中,上表面涡沿顺时针方向旋转,而下表面为逆时针。由于总距角的存在,使得下表面涡的强度总体上较上表面强。此外,从图6中还可以观察到,在后缘小翼向下摆动过程中,由于小翼对下表面涡的阻塞,使得两个涡发生了合并,表现出非定常特性。

G=

综上所述,可以得出结论:AFC旋翼翼型涡主要从后缘小翼与主桨叶间缝隙产生,上下表面涡旋转方向相反,并呈现显著的非定常特性。

其中COD、TN、TP作为主要污染物去除指标,其每年污染物减少量为:COD:438t/a;TN:291t/a;TP:18.3t/a;HN3-H:109.5t/a。

由于串线的每一步都是沿着正六边形的一条边进行运动,所以将每一条边赋值为(见图6).假设每圈运动起点坐标为,则有上行线和下行线在第圈运动到第步的坐标表达式,即当时,有;当时,有;当时,有;当时,有;当时,有;当时,有.

4 AFC翼型参数化特性分析

为探究影响涡产生的因素,并为日后开展AFC旋翼设计、运动控制方式设计等提供参考依据,有必要进行AFC翼型参数对涡的影响研究。

4.1 桨尖速度对涡的影响

在第3节的基准桨叶参数基础上,改变旋翼桨尖速度进行算例计算。图7给出了旋翼桨尖马赫数从Ma=0.4开始,以Ma=0.05步进、后缘小翼同一摆动角度下的涡量分布图。可以观测到,随着旋翼桨尖速度的增大,涡产生速度加快,且单个涡强度增大。

  

图7 不同桨尖马赫数下后缘小翼处涡量分布Fig.7 Vorticity distribution of trailing edge with different tip Mach numbers

4.2 总距角对涡的影响

图8给出了旋翼总距从2.5°开始,以2.5°步进、后缘小翼同一摆动角度下的涡量分布图。可观察到,随着旋翼总距角的增加,单个涡强度逐渐增大,且在迎角过大时主桨叶也将产生一定量的涡。

(2) 在网格块b中,初始化搜索的出发单元作为当前单元ic

  

图8 不同总距角下后缘小翼处涡量分布Fig.8 Vorticity distribution of trailing edge with different collective pitch angles

4.3 后缘小翼摆动频率对涡的影响

  

图9 不同后缘小翼摆动频率下涡量分布Fig.9 Vorticity distribution of trailing edge with different oscillating frequencies

AFC旋翼摆动周期通常为桨叶片数k与旋翼转速Ω乘积的整倍数。图9给出了后缘小翼摆动角速度2Ω,4Ω,6Ω,8Ω时后缘小翼同一摆动角度下的涡量分布图。图中不同摆动频率下观察到的涡数目几乎相同,即随着后缘小翼摆动频率的增加,单位时间内产生涡的数目也增多。另外,还可看出单个涡的强度略有减小。

4.4 后缘小翼摆动幅度对涡的影响

图10给出了后缘小翼摆动幅度从6°开始,以3°步进、后缘小翼同一摆动角度下的涡量分布图。对比可看出,随着摆动幅度的增加,涡的强度将明显增大。

  

图10 不同后缘小翼摆动幅度下涡量分布Fig.10 Vorticity distribution of trailing edge with different swing amplitudes

4.5 后缘小翼缝隙宽度对涡的影响

根据第3节涡的历程可知,后缘小翼与主桨叶间缝隙是绝大部分涡的形成区域,为此,本节分析缝隙宽度改变对涡的影响。图11给出了缝隙宽度在0.5%c~3%c,以0.5%c步进、后缘小翼同一摆动角度下的涡量分布图。可以观测到,从0.5%c开始,随着缝隙的增大,单个涡强度逐渐增大,但缝隙宽度3%c下单个涡的强度却较2.5%c时下降,即在一定范围内,增大缝隙宽度,单个涡的强度将增大。

  

图11 不同后缘小翼与主桨叶间间距下涡量分布Fig.11 Vorticity distribution of trailing edge with different gap widths

4.6 后缘小翼长度对涡的影响

图12给出了后缘小翼弦长在0.1c~2c,以0.025c步进、后缘小翼同一摆动角度下的涡量分布图。可见,随着后缘小翼弦长的增大,涡逐渐由连续的贴体涡转变为分离涡,单个涡强度逐渐增加,但涡产生频率逐渐减少。

  

图12 不同后缘小翼长度下涡量分布Fig.12 Vorticity distribution of trailing edge with different chord lengths

5 结 论

(1) AFC旋翼涡主要从后缘小翼与主桨叶间缝隙产生,对本文算例,涡在上表面顺时针旋转,而下表面旋向相反为逆时针,呈现显著非定常现象。

1.1 研究对象 本研究以上海市某地区失去独生子女的父母作为研究对象。纳入标准:(1)没有领养意愿及行为,已经丧失再生育能力的夫妇;(2)年龄≥50 岁;(3)失去独生子女 1年以上;(4)能独立完成问卷调查。排除认知障碍及患有重大精神疾病者。

本文建立的运动嵌套数值模拟方法满足AFC旋翼翼型的数值模拟要求。应用该方法对AFC旋翼翼型后缘小翼涡的数值模拟可得到如下结论:

(2) 提高桨尖马赫数、增加后缘小翼摆动频率能加快翼型涡产生速度;而提高桨尖马赫数、增大后缘小翼摆动幅度及增加后缘小翼长度能增大单个涡强度。

3、多家保险公司并进发展。2012年,我市由阳光农业保险公司一家开展试点,经过几年的发展,到2016年,形成了阳光互助农业保险公司、阳光农业保险公司、中国人保财险、平安保险公司等多家保险公司齐头并进,共同发展的良好势头,为保险部门在创新农业保险产品,提高服务质量等方面,起到了积极的促进作用。

(3) 在一定范围内,增加总距角能提高AFC旋翼翼型单个涡强度;增大后缘小翼与主桨叶缝隙间距也能提高单个涡强度,但间距过大时涡强度反而减少。

参考文献:

[1] FRIEDMANN P P, MILLOTT T A. Vibration reduction in rotorcraft using active control: A comparison of various approaches [J]. AIAA Journal of Guidance, Control, and Dynamics, 1995, 18(4):664-673.

[2] MILGRAM J H. A comprehensive aeroelastic analysis of helicopter main rotors with trailing edge flaps for vibration reduction [C] // 47th AIAA Aerospace Sciences Meeting and Including.[S.l.]:AIAA,1997:9736601.

[3] SCHMITZ F H, GOPALAN G. High-speed impulsive helicopter noise reduction possibilities through on-blade acoustic control [C] // Proceedings of the International Forum on Rotorcraft Multidisciplinary Technology.Seoul, Korea:[s.n.], 2007.

[4] AOYAMA T. Comparison of noise reduction effect between AFC and conventional IBC by moving overlapped grid method [C] // 27th AIAA Aeroacoustics Conference.[S.l.]:AIAA, 2006:8-10.

[5] MILLOTT T,FRIEDMANN P. Vibration reduction in helicopter rotors using an active control surface located on the blade[C]//33rd Structures, Structural Dynamics and Materials Conference. Dallas,TX,USA:[s.n.],1992.

坐在父母身边,父亲对母亲说:“你看老三这头发多好,油黑油黑的。”我冲二老一笑:“你们的老花眼哪能看到我冒出的白头发。”饭桌上,母亲一头白发,父亲一头稀发,我一头黑发,凑到一起大笑起来。墙上的老相框里,头发乌黑浓密的父亲、编着麻花长辫的母亲,一直面露他俩结婚时的喜悦笑容瞧着这个家,瞧了五十多年。

[6] 王荣,夏品奇. 多片后缘小翼对直升机旋翼桨叶动态失速及桨毂振动载荷的控制 [J]. 航空学报, 2013, 34(5):1083-1091.

WANG Rong, XIA Pinqi. Control of helicopter rotor blade dynamic stall and hub vibration by multiple trailing edge flaps [J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(5):1083-1091.

[7] 马奕扬,招启军,赵国庆. 基于后缘小翼的旋翼翼型动态失速控制分析 [J]. 航空学报, 2017, 38(3):127-137.

MA Yiyang, ZHAO Qijun, ZHAO Guoqing. Dynamic stall control of rotor airfoil via trailing-edge flap [J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(3):127-137.

[8] LÖHNER R, SHAROV D, LUO H. Overlapping unstructured grids [C]// 39th AIAA Aerospace Sciences Meeting and Exhibit.[S.l.]:AIAA, 2001:8-11.

[9] 田书玲. 基于非结构网格方法的重叠网格算法研究 [D].南京:南京航空航天大学, 2008.

TIAN Shuling. Investigation of overset unstructured grids algorithm [D]. Nanjing:Nanjing University of Aeronautics and Astronautics, 2008.

[10] PASCIONI K, CATTAFESTA L N, CHOUDHARI M M.An experimental investigation of the 30P30N multi-element high-lift airfoil [C]// 35th AIAA Aeroacoustics Conference.[S.l.]:AIAA, 2014, 204(1):1-13.

[11] CARADONNA F X, TUNG C. Experimental and analytical studies of a model helicopter rotor in hover [R]. NASA TM-81232, 1981.

2010年全国开展经济普查活动,作为公司统计,马成贵责无旁贷地担负起了此项任务。经济普查工作庞杂烦琐,需要大量的时间和精力,要对公司辖区内所有企业、个体工商户等经营户进行排查、摸底、登记、造册,要登记清楚每个辖区内的经营户成立的时间、注册资本、性质、经营范围等。由于公司所属煤矿分散在乌鲁木齐、托克逊和吐鲁番等辖区内,他要先走访了解各个煤矿辖区内的经营户有多少家,然后再去工商局调档进行登记。他没有抱怨工作忙、工作量大,而是默默地投入繁杂的工作中。没有帮手,他经常加班,双休日几乎没有休息过。在他的努力下,按时完成了各项工作。

 
胡志远,徐国华,史勇杰
《南京航空航天大学学报》2018年第02期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号