留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于非线性气动力的失速颤振计算与试验研究

戴玉婷 严慧 王林鹏

戴玉婷, 严慧, 王林鹏. 基于非线性气动力的失速颤振计算与试验研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.03.0141
引用本文: 戴玉婷, 严慧, 王林鹏. 基于非线性气动力的失速颤振计算与试验研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.03.0141
Yu-ting DAI, Hui YAN, Lin-peng WANG. CALCULATION AND EXPERIMENTAL STUDY OF STALL FLUTTER BASED ON NONLINEAR AERODYNAMICS[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.03.0141
Citation: Yu-ting DAI, Hui YAN, Lin-peng WANG. CALCULATION AND EXPERIMENTAL STUDY OF STALL FLUTTER BASED ON NONLINEAR AERODYNAMICS[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.03.0141

基于非线性气动力的失速颤振计算与试验研究

doi: 10.6052/j.issn.1000-4750.2019.03.0141
基金项目: 国家自然科学基金(11672018)
详细信息
    作者简介:

    严 慧(1996−),女,安徽人,硕士生,主要从事气动弹性力学研究(E-mail: hui1805423@buaa.edu.cn)

    王林鹏(1987−),男,辽宁人,博士生,主要从事气动弹性力学研究(E-mail: wanglinpeng@buaa.edu.cn)

    通讯作者: 戴玉婷(1985−),女,湖北人,副教授,主要从事气动弹性力学研究(E-mail: yutingdai@buaa.edu.cn)
  • 中图分类号: V215.34

CALCULATION AND EXPERIMENTAL STUDY OF STALL FLUTTER BASED ON NONLINEAR AERODYNAMICS

  • 摘要: 飞行器大攻角飞行过程中的动态失速会导致结构自激扭转或俯仰运动,造成非线性失速颤振现象,直接影响飞行器飞行安全与结构安全。该文对标准Leishman-Beddoes非线性非定常气动力模型进行马赫数修正,使其适用于低速不可压情形的动态失速气动力计算,然后基于二元翼段气动弹性模型,采用Newmark时域推进方法进行工程失速颤振计算。依据计算结果设计并完成了二元翼段失速颤振风洞试验。试验结果表明,多数试验状态,基于L-B模型的失速颤振计算结果与试验结果均吻合较好。结果验证了修正的L-B模型可以用来进行低速大展弦比平直翼段翼型的失速颤振工程分析与极限环振荡评估,同时,失速颤振速度与极限环幅值受初始攻角的影响很大。
  • 图  1  二元翼段气动弹性系统

    Figure  1.  Two-dimensional airfoil aeroelastic system

    图  2  攻角为$\alpha = {15^ \circ }{\rm{ + }}{10^ \circ }\sin \omega t$时翼型的法向力与力矩系数

    Figure  2.  Normal force and moment coefficient of airfoil when $\alpha = {15^ \circ }{\rm{ + }}{10^ \circ }\sin \omega t$

    图  3  气动弹性系统的俯仰分岔图和极限环振荡曲线图

    Figure  3.  Pitch bifurcation diagram and limit cycle oscillation graph of aeroelastic system

    图  4  翼段气弹系统风洞试验模型

    Figure  4.  Wind tunnel test model for wing airframe system

    图  5  初始攻角16度不同风速下翼段失速颤振试验曲线

    Figure  5.  Wing stall flutter test curve of different wind speed at initial angles of 16 degrees

    图  6  不同初始攻角下翼段俯仰角随风速变化的分岔曲线

    Figure  6.  Bifurcation curve of wing pitch angle with wind speed under different initial angles of attack

    图  7  初始攻角16度试验和数值模拟计算曲线对比

    Figure  7.  Test and numerical simulation calculation comparison curve at initial attack angle of 16 degrees

    图  8  初始攻角16度风速10 m/s试验和数值振荡频率曲线对比

    Figure  8.  Comparison curve of test and numerical oscillation frequency with initial attack angle of 16 degrees and wind speed of 10 m/s

    图  9  不同初始攻角下翼段俯仰角随风速变化的试验和计算分岔曲线对比

    Figure  9.  Comparison of test and calculated bifurcation curves of wing pitch angle with wind speed under different initial angles of attack

    表  1  试验重要参数

    Table  1.   Important parameters of the test

    参数值参数值
    C0.3 ml0.9 m
    m16.69 kg/mρ1.225 kg/m3
    Iθ0.31 kgm2Kh30.5 N/mm
    Kθ13.1 Nm/radPitch axis0.115 m
    下载: 导出CSV
  • [1] 吴长波, 周拜豪, 陆庆飞, 崔海涛. 小展弦比压气机转子叶片颤振研究[J]. 航空动力学报, 2016, 31(4): 1000 − 1005.

    Wu Changbo, Zhou Baihao, Lu Qingfei, Cui Haitao. Study on rotor blade flutter of compressor with small aspect ratio [J]. Journal of Aerospace Power, 2016, 31(4): 1000 − 1005. (in Chinese)
    [2] 刘廷瑞, 于子晴. 风力机叶片失速非线性颤振伺服气弹智能控制[J]. 中南大学学报(自然科学版), 2016, 47(10): 3562 − 3569.

    Liu Tingrui, Yu Ziqing. Wind turbine blade stall nonlinear flutter servo aeroelastic intelligent control [J]. Journal of Central South University (natural science edition), 2016, 47(10): 3562 − 3569. (in Chinese)
    [3] Sheng, Galbraith, Mcd R A, F. N. Return from airfoil stall during ramp-down pitching motions [J]. Journal of Aircraft, 2007, 44(6).
    [4] 谢丹, 徐敏. 基于特征正交分解的三维壁板非线性颤振分析[J]. 工程力学, 2015, 32(1): 1 − 9.

    Xie Dan, Xu Min. Three-dimensional panel nonlinear flutter analysis based on proper orthogonal decomposition method [J]. Engineering Mechanics, 2015, 32(1): 1 − 9. (in Chinese)
    [5] Lawrence, W., Carr, And, M. S., Chandrasekhara. Compressibility effects on dynamic stall[J].
    [6] John, A., Ekaterinaris, And, Max, F., Platzer. Computational prediction of airfoil dynamic stall[J].
    [7] Guilmineau E, Queutey P. Numerical study of dynamic stall on several airfoil sections [J]. Aiaa Journal, 2015, 37(1): 128 − 130. doi:  10.2514/2.677
    [8] Geissler W, Dietz G, Mai H. Dynamic stall on a supercritical airfoil [J]. Aerospace Science & Technology, 2015, 9(5): 390 − 399.
    [9] Spentzos A, Barakos G N, Badcock K J, Richards B E, Wernert P, Schreck S, Raffel M. Investigation of three-dimensional dynamic stall using computational fluid dynamics [J]. Aiaa Journal, 2015, 43(5): 1023 − 1033. doi:  10.2514/1.8830
    [10] Dimitriadis G, Li J. Bifurcation behavior of airfoil undergoing stall flutter oscillations in low-speed wind tunnel [J]. Aiaa Journal, 2009, 47(47): 2577 − 2596.
    [11] Yabili S, Smith M J, Dimitriadis G. Unsteady Navier-Stokes simulation of low-Reynolds stall flutter[C]//Aiaa Aerospace Sciences Meeting and Exhibit, 2012.
    [12] Sheng W, Galbraith R a M, Coton F N. A modified dynamic stall model for low Mach numbers [J]. Journal of Solar Energy Engineering, 2007, 130(3): 653 − 653.
    [13] Galbraith R. Return from Airfoil Stall During Ramp-Down Pitching Motions [J]. Journal of Aircraft, 2007, 44(6).
    [14] 任勇生, 刘廷瑞. 具有结构阻尼的复合材料薄壁梁的动力失速非线性颤振特性[J]. 振动与冲击, 2013, 32(18): 146 − 152. doi:  10.3969/j.issn.1000-3835.2013.18.028

    Ren Yongsheng, Liu Tingrui. Nonlinear stall flutter characteristics of composite thin-walled beams with structural damping [J]. Journal of Vibration and Shock, 2013, 32(18): 146 − 152. (in Chinese) doi:  10.3969/j.issn.1000-3835.2013.18.028
    [15] Costes M, Richez F, Pape A L, Gavériaux R. Numerical investigation of three-dimensional effects during dynamic stall [J]. Aerospace Science & Technology, 2015, 47(3): 216 − 237.
    [16] Chen X, Zha G C, Yang M T. Numerical simulation of 3-D wing flutter with fully coupled fluid-structural interaction [J]. Computers & Fluids, 2007, 36(5): 856 − 867.
    [17] Zanotti A, Nilifard R, Gibertini G, Guardone A, Quaranta G. Assessment of 2D/3D numerical modeling for deep dynamic stall experiments [J]. Journal of Fluids & Structures, 2007, 51: 97 − 115.
    [18] Naughton J, Strike J, Hind M, Magstadt A, Babbitt A. Measurements of dynamic stall on the DU wind turbine airfoil series [C]//AHS Forum, 2013.
    [19] Bhat S S, Govardhan R N. Stall flutter of NACA 0012 airfoil at low Reynolds numbers [J]. Journal of Fluids & Structures, 2013, 41(August 2013): 166 − 174.
    [20] Uruba V. Near wake dynamics around a vibrating airfoil by means of PIV and Oscillation Pattern Decomposition at Reynolds number of 65000 [J]. Journal of Fluids & Structures, 2015, 55: 372 − 383.
    [21] Šidlof P, Vlček V, Štěpán M. Experimental investigation of flow-induced vibration of a pitch-plunge NACA 0015 airfoil under deep dynamic stall [J]. Journal of Fluids and Structures, 2016, 67: 48 − 59. doi:  10.1016/j.jfluidstructs.2016.08.011
    [22] Song S, Qinghua Z, Zhang C, Xianping N. Airfoil aeroelastic flutter analysis based on modified Leishman-Beddoes model at low Mach number [J]. Chinese Journal of Aeronautics, 2011, 24(5): 550 − 557. doi:  10.1016/S1000-9361(11)60064-2
  • [1] 任若松, 梁新华, 刘小兵, 马文勇, 刘庆宽.  准流线型桥梁断面气动力特性的雷诺数效应研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2019.04.S023
    [2] 周佳豪, 马文勇, 黄伯城.  临界雷诺数区光滑圆柱振动与气动力研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2018.05.S050
    [3] 梁洪超, 楼文娟, 丁浩, 卞荣.  非线性振型结构HFFB试验模态力计算方法及不确定性分析 . 工程力学, doi: 10.6052/j.issn.1000-4750.2018.01.0008
    [4] 张景钰, 张明金, 李永乐, 房忱, 向活跃.  高速铁路路堤-路堑过渡段复杂风场和列车气动效应风洞试验研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2017.09.0749
    [5] 李尚斌, 林永峰, 樊枫.  倾转旋翼气动特性风洞试验与数值模拟研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2017.07.0561
    [6] 郑云飞, 刘庆宽, 刘小兵, 马文勇.  端部状态对斜拉索节段模型气动特性的影响 . 工程力学, doi: 10.6052/j.issn.1000-4750.2016.03.S036
    [7] 苏洋, 李永乐, 彭栋, 陈科宇, 孟新利.  公路风屏障流场特性及自身风荷载的足尺模型风洞试验研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2016.08.0588
    [8] 徐昕宇, 李永乐, 廖海黎, 任森.  双层桥面桁架梁三塔悬索桥颤振性能优化风洞试验 . 工程力学, doi: 10.6052/j.issn.1000-4750.2015.12.0953
    [9] 张志田, 张显雄, 陈政清.  桥梁气动力CFD模拟中湍流模型的应用现状 . 工程力学, doi: 10.6052/j.issn.1000-4750.2015.10.0821
    [10] 周晅毅, 祖公博, 顾明.  TTU标准模型表面风压大涡模拟及风洞试验的对比研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2014.06.0504
    [11] 李永乐, 徐昕宇, 郭建明, 向活跃, 陈克坚.  六线双层铁路钢桁桥车桥系统气动特性风洞试验研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2014.09.0741
    [12] 谢长川, 胡锐, 王斐, 刘燚, 常楠.  大展弦比柔性机翼气动弹性风洞模型设计与试验验证 . 工程力学, doi: 10.6052/j.issn.1000-4750.2015.04.0254
    [13] 马文勇, 袁欣欣, 张晓斌, 尉耀元, 刘庆宽.  圆形断面在35k~330k雷诺数范围的气动力特性研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2014.05.S022
    [14] 向活跃, 李永乐, 胡喆, 廖海黎.  桥上孔隙式风屏障缩尺模型模拟方法的风洞试验 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.05.0321
    [15] 刘晓燕, 杨 超, 吴志刚.  适用于气动弹性的气动力自适应辨识方法 . 工程力学,
    [16] 李永乐, 王 涛, 廖海黎.  斜拉桥并列拉索尾流驰振风洞试验研究 . 工程力学,
    [17] 梁枢果, 邹良浩, 郭必武.  基于刚性模型测压风洞试验的 武汉国际证券大厦三维风致响应分析 . 工程力学,
    [18] 张 伟, 葛耀君.  速度场演化结构非定常气动力 . 工程力学,
    [19] 吕 斌, 吴志刚, 杨 超.  T型尾翼颤振特性分析方法 . 工程力学,
    [20] 何 辉, 张伟伟, 叶正寅.  基于动态近似边界条件的气动弹性数值模拟 . 工程力学,
  • 加载中
计量
  • 文章访问数:  29
  • HTML全文浏览量:  22
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-30
  • 修回日期:  2020-04-13
  • 网络出版日期:  2020-06-02

基于非线性气动力的失速颤振计算与试验研究

doi: 10.6052/j.issn.1000-4750.2019.03.0141
    基金项目:  国家自然科学基金(11672018)
    作者简介:

    严 慧(1996−),女,安徽人,硕士生,主要从事气动弹性力学研究(E-mail: hui1805423@buaa.edu.cn)

    王林鹏(1987−),男,辽宁人,博士生,主要从事气动弹性力学研究(E-mail: wanglinpeng@buaa.edu.cn)

    通讯作者: 戴玉婷(1985−),女,湖北人,副教授,主要从事气动弹性力学研究(E-mail: yutingdai@buaa.edu.cn)
  • 中图分类号: V215.34

摘要: 飞行器大攻角飞行过程中的动态失速会导致结构自激扭转或俯仰运动,造成非线性失速颤振现象,直接影响飞行器飞行安全与结构安全。该文对标准Leishman-Beddoes非线性非定常气动力模型进行马赫数修正,使其适用于低速不可压情形的动态失速气动力计算,然后基于二元翼段气动弹性模型,采用Newmark时域推进方法进行工程失速颤振计算。依据计算结果设计并完成了二元翼段失速颤振风洞试验。试验结果表明,多数试验状态,基于L-B模型的失速颤振计算结果与试验结果均吻合较好。结果验证了修正的L-B模型可以用来进行低速大展弦比平直翼段翼型的失速颤振工程分析与极限环振荡评估,同时,失速颤振速度与极限环幅值受初始攻角的影响很大。

English Abstract

戴玉婷, 严慧, 王林鹏. 基于非线性气动力的失速颤振计算与试验研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.03.0141
引用本文: 戴玉婷, 严慧, 王林鹏. 基于非线性气动力的失速颤振计算与试验研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.03.0141
Yu-ting DAI, Hui YAN, Lin-peng WANG. CALCULATION AND EXPERIMENTAL STUDY OF STALL FLUTTER BASED ON NONLINEAR AERODYNAMICS[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.03.0141
Citation: Yu-ting DAI, Hui YAN, Lin-peng WANG. CALCULATION AND EXPERIMENTAL STUDY OF STALL FLUTTER BASED ON NONLINEAR AERODYNAMICS[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.03.0141
  • 自1980年以后,科学界对动态失速和失速颤振的仿真计算有很大兴趣,各种计算模型不断发展。在过去几十年中,许多研究人员通过实验,数值模拟等方式[1-4]开展了研究。目前人们普遍认为没有通用的动态失速模型,即使是先进的CFD计算[5- 9]也不能提供一致的结果。2009年,Dimitriadis和Li[10]对NACA0012翼型的俯仰和沉浮自由度做了一些试验研究。证明沉浮运动对失速颤振的发生没有影响,并得出失速颤振呈现出不同的分支的结论。2012年,SachaYabili[11]等利用自己制定的求解器对二元翼段低雷诺数下的动态失速和失速颤振进行了CFD计算,并与Li的试验结果做了对比验证。

    对于工程应用而言,半经验动态失速模型在直升机和风力涡轮机空气动力学中仍然普遍存在。在这些半经验动态失速模型中,Leishman-Beddoes模型是最受欢迎的模型之一,主要适用于0.3 Ma-0.8 Ma。2007年,W.Sheng等[12]提出了一种改进的动态失速模型,该模型适用于0.3 Ma以下气动力计算。由Sheng等人[13]开发的一种新的失速发生标准已被用来取代基于Evans-Mort相关性的原始标准。任勇生等[14]研究了复材薄壁梁非线性失速颤振特性,得出结构阻尼在抑制颤振、增强结构稳定性方面效果明显。

    失速颤振试验相关研究也从动态失速[15]、试验规律、影响因素[16]等方面有了新的进展。随着计算机,激光和相关技术的发展,二维和三维粒子图像测速(PIV)技术成功地用于测量翼型周围的流场速度[17-18],分析非定常气动力产生机理。2013年Bhat和Govardhan[19]使用PIV和测压元件测量了NACA 0012俯仰振荡翼型周围流场和非定常气动力。2015年Uruba[20]通过振荡模式分解(OPD)方法分析了俯仰振荡的NACA 0012翼型周围的PIV速度场。2016年Petr等人[21]展示了一种自振荡俯仰的翼型设计方法,该模型适用于中等雷诺数下流动引起的振动。

    在Sheng[12]的修改模型基础上,本文对标准Leishman-Beddoes非线性非定常气动力模型参数进行了修正,建立二维翼段气动弹性模型进行失速颤振计算,并根据计算结果设计并完成了失速颤振风洞试验,进一步验证了修正的L-B模型可以用来进行低速翼型的失速颤振工程分析与极限环振荡评估。

    • 原始L-B模型适用的速度范围在0.3 Ma-0.8 Ma之间,为了使该非定常气动力模型适用于低速情况,本文在Beddoes和Sheng[3, 11]研究基础上,充分考虑翼型表面涡对气动载荷的影响,在低马赫数条件下对L-B模型进行修正,使其适用于低速条件的非定常气动力计算。

      翼型的附面层的气流扰动需要一段时间才能发展成一定强度的涡,从而使前缘压力超过临界值而发生失速。为了将其计入计算中,本文借鉴原模型中延迟的方法,对${C'_{\rm{N}}}$进行一阶延迟得到二次延迟的前缘压力系数${C''_{\rm{N}}}$,其延迟计算可由下式表示:

      $${C''_{\rm{N}}} = {C'_{\rm{N}}} - {D_{{b_n}}}\quad\qquad\qquad\qquad\qquad\qquad\qquad$$ (1)
      $${D_{{b_n}}} = {D_{{b_{n - 1}}}}\exp \left(\frac{{ - \Delta S}}{{{T_b}}}\right) + ({C'_{{{\rm{N}}_n}}} - {C'_{{{\rm{N}}_{n - 1}}}})\exp \left(\frac{{ - \Delta S}}{{2{T_b}}}\right)$$ (2)

      式中${D_{{b_n}}}$为中间变量,${T_b}$为时间延迟常数,$S$为半弦长下涡移动的距离。

      在上升沿失速过程中,涡的形成和分离靠近前缘位置,在低马赫数时附面层在分离的时候仍有部分附着在翼型的上表面。因此,其附着效应将为上表面带来一个附加升力,称之为“超调”。超调升力是低马赫数下特有的量,该部分的升力系数可由下式计算[22]

      $$\Delta C_{\rm{N}}^{\rm{v}} = {B_1}(f'' - f){V_{\rm{x}}}$$ (3)

      式中${B_1}{\rm{ = }}1.0$为与翼型有关的参量;$f$为弦向分离点的位置;$f''$为初始模型中对于分离点$f$的一阶修正;${V_{\rm{x}}}$为由涡产生的法向力系数的形函数。

      $${V_{\rm{x}}} = \left\{ {\begin{array}{*{20}{c}} {{{\sin }^{3/2}}\left( {\dfrac{{\pi \tau }}{{2{T_{\rm{v}}}}}} \right),\begin{array}{*{20}{c}} {}&{0 < \tau < {T_{\rm{v}}}} \end{array}} \\ {{{\cos }^2}\left( {\dfrac{{\pi (\tau - {T_{\rm{v}}})}}{{{T_{\rm{vl}}}}}} \right),\begin{array}{*{20}{c}} {}&{\tau > {T_{\rm{v}}}} \end{array}} \end{array}} \right.$$ (4)

      式中$\tau $为无量纲时间;${T_{\rm{v}}}{\rm{ = }}7.2$为沿弦向传播的时间常数;${T_{\rm{vl}}}{\rm{ = }}4.5$为涡传播距离的时间常数。

      翼面上由附加效应产生的附加力矩系数$\Delta C_{\rm{m}}^{\rm{v}}$可由下式计算:

      $$\Delta C_{\rm{m}}^{\rm{v}} = {B_2}\left( {1 - \cos \frac{{\pi {\tau _{\rm{v}}}}}{{{T_{\rm{vl}}}}}} \right)\Delta C_{\rm{N}}^{\rm{v}}$$ (5)

      式中${B_2}{\rm{ = }}0.32$为翼型相关参量。${\tau _{\rm{v}}}$为涡传递的无量纲时间参数。

      在下降途中,气流重新附着到翼型表面时,会受到涡的影响使得附着延迟,进而实际气动力要小于未做延迟处理的模型所得值。该差值为下降沿的超调值。${\alpha _{\min 0}}{\rm{ = }}16.57$°表示下降沿气流再次附着阶段的起始角度,${T_{\rm{r }}}{\rm{ = }} 7.02$为气流再次附着阶段涡运动时间常数,当$\alpha < {\alpha _{\min 0}}$时气流再次进入附着阶段,翼型上表面涡运动停止,气流由分离进入完全附着阶段。在下降沿中将超调过程开始时对应的攻角设为${\alpha _{\min 0}}$,则超调过程与超调开始对应的攻角之间的关系见下式:

      $${\alpha _{\min }} = {\alpha _{\min 0}} + {T_{\rm{r}}}q$$ (6)

      式中${T_{\rm{r}}}$为重新附着过程中的时间延迟常数,$q$为减缩频率。

      则回程过程中,在下降沿中法向力系数的超调值可写成:

      $$\Delta C_{\rm{nr}}^{\rm{v}} = {B_1}(f'' - f){V_{\rm{xr}}}$$ (7)

      式中${V_{\rm{xr}}}$为回程中基于翼型的法向力系数的形函数,其计算式为

      $${V_{\rm{xr}}} = \left\{ {\begin{array}{*{20}{c}} {{{\sin }^{3/2}}\left( {\dfrac{{\pi {\tau _{\rm{r}}}}}{{2{T_{\rm{r}}}}}} \right),\begin{array}{*{20}{c}} {}&{0 < {\tau _{\rm{r}}} < {T_{\rm{r}}}} \end{array}} \\ {{{\cos }^2}\left( {\dfrac{{\pi ({\tau _{\rm{r}}} - {T_{\rm{r}}})}}{{2{T_{\rm{r}}}}}} \right),\begin{array}{*{20}{c}} {}&{{\tau _{\rm{r}}} > {T_{\rm{r}}}} \end{array}} \end{array}} \right.$$ (8)

      式中${\tau _{\rm{r}}}$为回程的无量纲时间。

      对应的力矩系数:

      $$\Delta C_{\rm{mr}}^{\rm{v}} = {B_2}\left( {1 - \cos \frac{{\pi {\tau _{\rm{v}}}}}{{{T_{{\rm{vl}}}}}}} \right)\Delta C_{{\rm{Nr}}}^{\rm{v}}$$ (9)

      低速修正模型最终的法向力系数和力矩系数:

      $${C_{\rm{N}}}(t) = C_{\rm{N}}^{\rm{f}}(t) + C_{\rm{N}}^{\rm{v}}(t) + \Delta C_{\rm{N}}^{\rm{v}} + \Delta C_{\rm{nr}}^{\rm{v}}\;\;$$ (10)
      $${C_{\rm{M}}}(t) = C_{\rm{M}}^{\rm{f}}(t) + C_{\rm{M}}^{\rm{v}}(t) + \Delta C_{\rm{m}}^{\rm{v}} + \Delta C_{\rm{mr}}^{\rm{v}}$$ (11)

      其中$\Delta C_{\rm{N}}^{\rm{v}}$$\Delta C_{\rm{nr}}^{\rm{v}}$表示去程和回程涡升力的超调值,这是低速情况与亚声速情况L-B模型的不同之处。

    • 图1为二维翼段气动弹性系统,该系统具有俯仰和沉浮自由度。综合1.1节非定常气动力计算方法,则该系统的气动弹性方程可以写为:

      $$\begin{split}\left[ {\begin{array}{*{20}{c}} {{Q_{\rm{h}}}} \\ {{Q_{{\text{θ}}} }} \end{array}} \right] =& \left[ {\begin{array}{*{20}{c}} m&{{S_{\text{θ}} }} \\ {{S_{\text{θ}} }}&{{I_{\text{θ}} }} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\ddot h} \\ {\ddot \theta } \end{array}} \right] + \\&\left[ {\begin{array}{*{20}{c}} {{C_{\rm{h}}}}&0 \\ 0&{{C_{\text{θ}} }} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\dot h} \\ {\dot \theta } \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{K_{\rm{h}}}}&0 \\ 0&{{K_{\text{θ}} }} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} h \\ \theta \end{array}} \right]\end{split}$$ (12)

      式中$\theta $为翼型的俯仰角;$h$为翼型的垂向位移;$m$为翼型单位长度质量;${K_{\rm{h}}}$${K_{{\text{θ}}}} $分别为沉浮和俯仰弹簧刚度;${C_{\rm{h}}}$${C_{\text{θ}} }$分别为沉浮和俯仰阻尼系数;${S_{\text{θ}} }$为质量静距;${I_{\text{θ}} }$为翼型相对1/4弦线处的极惯性矩;${Q_{\rm{h}}}$${Q_{\text{θ}} }$分别为外部气动力和外部气动力矩;其表达式在1.1节中已求得,具体如下:

      图  1  二元翼段气动弹性系统

      Figure 1.  Two-dimensional airfoil aeroelastic system

      $$\left[ {\begin{array}{*{20}{c}} {{Q_{\rm{h}}}} \\ {{Q_{\text{θ}} }} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} { - \dfrac{1}{2}\rho {V^2}{C_{\rm{L}}}cl} \\ {\dfrac{1}{2}\rho {V^2}{C_{\rm{M}}}{c^2}l + \dfrac{1}{2}\rho {V^2}{C_{\rm{L}}}cl{x_{\rm{h}}}} \end{array}} \right]$$ (13)

      式中的CLCM分别为升力系数和气动力矩系数;c为弦长;l为展长;xh为弹性轴距1/4弦线的距离;ρ为空气密度;V为来流速度。

      L-B模型中的攻角α和俯仰减缩角速率q可由下式进行计算:

      $$\left[ {\begin{array}{*{20}{c}} \alpha \\ q \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\theta + \dfrac{{\dot h}}{V}} \\ {\dfrac{{\dot \theta c}}{V}} \end{array}} \right]$$ (14)

      将方程(3)改写为时域状态空间方程形式,可得到:

      $${{A}}{\dot{ x}} = {{Bx}} + {{C}}$$ (15)

      式中

      $$ \begin{array}{l} {{{{x}} = }}\left[ \begin{array}{l} {\dot h} \\ {\dot \theta } \\ h \\ \theta \\ \end{array} \right], \;\;{{A}} = \left[ {\begin{array}{*{20}{c}} m&{{S_{\text{θ}} }}&0&0 \\ {{S_{\text{θ}} }}&{{I_{\text{θ}} }}&0&0 \\ 0&0&1&0 \\ 0&0&0&1 \end{array}} \right],\\ {{{{B}} = }}\left[ {\begin{array}{*{20}{c}} { - {C_{\rm{h}}}}&0&{ - {K_{\rm{h}}}}&0 \\ 0&{ - {C_{\text{θ}} }}&0&{ - {K_{\text{θ}} }} \\ 1&0&0&0 \\ 0&1&0&0 \end{array}} \right], \;\;{{{{C}} = }}\left[ \begin{array}{l} {Q_{\rm{h}}} \\ {Q_{\text{θ}} } \\ 0 \\ 0 \\ \end{array} \right] \end{array} $$

      当给定初始条件$x$,利用显式或者隐式积分可由求得其时域响应:

      $$\frac{{\delta {{{x}}^n}}}{{\delta t}} = {{{{ A}}}^{ - 1}}\left( {{{B}}{{{x}}^n} + {{C}}} \right)$$ (16)

      在本文中,使用四阶定步长龙格库塔法求解该状态空间方程。

    • 为验证低马赫数下的修正模型,仍然选用NACA0012翼型。并将计算结果与原始模型和试验以及Sheng的模型的结果[3]做了对比。验证中来流速度为0.1马赫。对比图中给出了深度失速条件(攻角$\alpha = {15^ \circ }{\rm{ + }}{10^ \circ }\sin \omega t$)下修正L-B模型的计算结果与试验对比,对比结果如图2所示。

      图  2  攻角为$\alpha = {15^ \circ }{\rm{ + }}{10^ \circ }\sin \omega t$时翼型的法向力与力矩系数

      Figure 2.  Normal force and moment coefficient of airfoil when $\alpha = {15^ \circ }{\rm{ + }}{10^ \circ }\sin \omega t$

      从图中可以看出,在低马赫数下,原始${\rm{L - B}}$模型计算数据较试验数据有较大差异,主要体现在超调位置,原始模型和试验数据重合性较差。而修正之后的${\rm{L - B}}$模型在超调位置上与试验数据吻合的较好。说明修正之后的${\rm{L - B}}$模型有较好的精度,能够用来计算翼段在低速情况下的气动力特性。模型中俯仰力矩系数还受${B_2}$, ${T_{\rm{vl}}}$等因素影响,所以力矩对的不是特别好。

    • 本文选用参考文献[10]中的试验条件和试验结果作为数值算例的参数,来验证失速颤振计算结果的正确性。文献展示了在低马赫数下基于NACA0012翼型的失速振荡试验。该试验的重要参数见下表1

      本文计算了二维翼段失速颤振的俯仰分岔图和极限环振荡曲线,并与试验[10]的计算结果进行了对比,对比如图3所示。

      表 1  试验重要参数

      Table 1.  Important parameters of the test

      参数值参数值
      C0.3 ml0.9 m
      m16.69 kg/mρ1.225 kg/m3
      Iθ0.31 kgm2Kh30.5 N/mm
      Kθ13.1 Nm/radPitch axis0.115 m

      图  3  气动弹性系统的俯仰分岔图和极限环振荡曲线图

      Figure 3.  Pitch bifurcation diagram and limit cycle oscillation graph of aeroelastic system

      从分岔图中可以看出,当风速为12.8 m/s,翼段开始周期性俯仰振荡,即开始分岔。当风速较小时,尤其当风速V低于15 m/s时,计算数据与试验数据吻合较好,但是在风速V大于15 m/s时,俯仰角振荡幅值明显高于试验值。且从试验数据来看,随着风速的增加俯仰振荡幅值未有明显变化。这与结构中沉浮刚度和俯仰刚度可能存在非线性或阻尼非线性等因素有关。

    • 为了探究失速颤振现象的形成机理,根据风洞的参数,设计并制作了二维翼段失速颤振的风洞试验模型,进行了多种状态下的风洞试验。失速颤振试验在北京航空航天大学D4开口风洞中完成,开口段尺寸为1.5 m。开口风速范围小于70 m/s。

      翼段翼型采用与仿真条件相同的NACA0012翼型,弦长为0.36 m,翼展为1 m。翼段翼肋之间相距0.125 m。另外,在翼弦向20%和70%处分别布置前梁和后梁。厚度和展长能保证最大瞬时风洞堵塞因子小于5%,满足试验要求。在前梁两端处安装加速度传感器,转轴处安装角位移计。

      翼段质量为2.363 kg,转轴位于弦长的37.4%位置处,重心处于弦长的39.04%位置,绕转轴的转动惯量${I_{\text{θ}} }$$0.135\;{\rm{kg}}{{\rm{m}}^{\rm{2}}}$。支撑系统为转轴滑轨系统。图4为试验装置模型,其中,图4a为试验系统的原理图,图4b为试验系统的实物图。图4a中,1为八个线性弹簧,2为四只初始攻角调节器,3为两组滑块系统,4为沉浮刚度调节的线性弹簧,5为支撑系统,6为两组转轴摇杆机构,7为角位移传感器,8为翼段转轴,9为加速度传感器,10为翼段。最终翼段沉浮刚度${k_{\rm{h}}}$$4896.595\;{\rm{kg/m}}$,俯仰刚度${k_{\text{θ}} }$$135.28\;{\rm{Nm/rad}}$。角位移传感器的精度为0.1度。

      图  4  翼段气弹系统风洞试验模型

      Figure 4.  Wind tunnel test model for wing airframe system

      试验目的是通过给定初始攻角,测量翼段发生颤振的速度,以及发生颤振后翼段的俯仰振幅。在试验过程中,通过调节拉杆两端连接器(图4a中2)的长度,由角度传感器测量结果,将初始攻角调节到要求的位置。然后给定一系列稳定风速,测量该初始攻角和给定速度下的角加速度和角位移,判断失速颤振特性。

    • 试验就四个初始攻角进行了测试,依次为7.5度、10度、13度、16度。试验风速从6 m/s到20 m/s。图5为16度初始攻角时,在翼段转轴处的角度传感器在不同风速下测得角位移动态曲线图。从下列图中可以看出在风速为6 m/s时,曲线几乎平直,幅值近乎为0。风速增加时,角度幅值信号随风速的增加而增加,且呈等幅震荡形式。将不同风速下的俯仰角振荡曲线幅值提取出来,则可生成攻角为16度时翼段俯仰角度幅值随风速变化的分岔图,分岔图见图6图6为不同初始攻角下翼段俯仰角随风速变化的分岔曲线图,从图中可以看出,在初始攻角为16度时,翼段发生俯仰自由度等幅振荡的风速为6.5 m/s,即分岔风速为6.5 m/s。在初始攻角一定的情况下,俯仰角的幅值随风速的增加而增加。另外,当初始攻角增大时,俯仰角的分岔速度明显降低。

      图  5  初始攻角16度不同风速下翼段失速颤振试验曲线

      Figure 5.  Wing stall flutter test curve of different wind speed at initial angles of 16 degrees

      图  6  不同初始攻角下翼段俯仰角随风速变化的分岔曲线

      Figure 6.  Bifurcation curve of wing pitch angle with wind speed under different initial angles of attack

      根据第2节二维气动弹性计算模型与试验模型参数,对试验模型失速颤振特性进行了数值模拟。模型的阻尼大小按照经验取0.1倍的刚度。图7为初始攻角为16度时风洞试验与计算模型在不同风速下的对比图。其中图7a为风速8 m/s时俯仰角响应,图7b为风速10 m/s时的俯仰角响应。从图中可以看出,仿真计算振幅和试风洞试验振幅吻合度较好。在风速较小时,受风速调节控制系统、开口外界环境影响和支撑支座的影响,气动弹性时域响应有一定的波动。图8为初始攻角16度、风速10 m/s试验和仿真的振荡频率曲线对比,从图中可以看出,振荡频率吻合较好。图9为不同初始攻角下,试验测得的翼段俯仰角等幅振荡的幅值随风速变化的分岔曲线与数值计算的分岔曲线的对比图。从图中可以看出,除了在初始攻角为7.5度时稍有偏差,数值计算的结果和试验测得结果吻合较好。尤其是初始攻角较大时,吻合度较高。说明修正后的L-B模型能够较好应用于大攻角失速颤振分析。

      图  7  初始攻角16度试验和数值模拟计算曲线对比

      Figure 7.  Test and numerical simulation calculation comparison curve at initial attack angle of 16 degrees

      图  8  初始攻角16度风速10 m/s试验和数值振荡频率曲线对比

      Figure 8.  Comparison curve of test and numerical oscillation frequency with initial attack angle of 16 degrees and wind speed of 10 m/s

      图  9  不同初始攻角下翼段俯仰角随风速变化的试验和计算分岔曲线对比

      Figure 9.  Comparison of test and calculated bifurcation curves of wing pitch angle with wind speed under different initial angles of attack

    • (1) 修正的L-B模型可以用于二元翼段或低速大展弦比飞机大迎角失速时的非定常气动力分析。但需要指出的是,L-B模型中超调时间参数对结果影响较大,应仔细选择。

      (2) 设计了二元翼段大攻角失速颤振试验,试验表明初始攻角越大,失速颤振分岔速度越低,极限环振荡幅值越大。风洞试验结果与非线性失速颤振理论计算结果吻合较好。说明基于修正L-B失速气动力模型具有工程上可以接受的精度。

目录

    /

    返回文章
    返回