《万方数据-数字化期刊群》全文上网期刊
CNKI《中国学术期刊(网络版)》全文收录期刊
《中文科技期刊数据库》(维普网)全文收录期刊
超星期刊域出版平台、博看网全文收录期刊
日本科学技术振兴机构数据库收录

留言板

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

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

空气舵气动力‒脉动压力‒结构耦合响应分析

郭静 张忠 谷迎松 王建民 李海波 高博

郭静, 张忠, 谷迎松, 等. 空气舵气动力‒脉动压力‒结构耦合响应分析[J]. 航天器环境工程, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004
引用本文: 郭静, 张忠, 谷迎松, 等. 空气舵气动力‒脉动压力‒结构耦合响应分析[J]. 航天器环境工程, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004
GUO J, ZHANG Z, GU Y S, et al. Analysis of aerodynamic force- fluctuating pressure- structure coupled response for air rudder[J]. Spacecraft Environment Engineering, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004
Citation: GUO J, ZHANG Z, GU Y S, et al. Analysis of aerodynamic force- fluctuating pressure- structure coupled response for air rudder[J]. Spacecraft Environment Engineering, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004

空气舵气动力‒脉动压力‒结构耦合响应分析

doi: 10.12126/see.2020.03.004
基金项目: 国家自然科学基金项目(编号:11502024;11502023;11402028;51805036)
详细信息
    作者简介:

    郭 静(1983—),女,博士学位,高级工程师,主要从事热气动弹性仿真与试验技术研究。E-mail:guojing2662632@126.com

  • 中图分类号: TB115.7; V414.8

Analysis of aerodynamic force- fluctuating pressure- structure coupled response for air rudder

  • 摘要: 气动力、脉动压力、结构振动相互作用,组成复杂的多场耦合系统,给动力学分析带来极大的挑战。文章基于PCL和DMAP语言自主研发了气动力–脉动压力–结构耦合响应分析软件,以复合材料空气舵为研究对象,建立其有限元模型,并开展模态分析;进而,建立基于Van Dyke修正活塞理论的气动模型,基于模态法分析了气动力–脉动压力–结构三者耦合的空气舵响应,并与不考虑气动力效应的非耦合结构响应进行对比,探究了气动力耦合效应对空气舵响应的影响规律。结果表明,气动弹性效应能使得空气舵振动响应从随机振动变为发散极限环振荡形式的高阶运动,显著改变脉动压力响应谱。可以预测,结构声疲劳分析中必须考虑气动弹性效应。
  • 图  1  复合材料空气舵有限元模型

    Figure  1.  FEM of air rudder made of composite materials

    图  2  复合材料空气舵固有频率和模态振型

    Figure  2.  Natural frequencies and modal shapes of composite air rudder

    图  3  气动力–脉动压力–结构耦合响应分析软件界面

    Figure  3.  Interface of software for analyzing aerodynamic force- fluctuating pressure- structure coupled response

    图  4  脉动压力时域曲线

    Figure  4.  Time domain curve of the fluctuating pressure

    图  5  耦合与非耦合加速度频域响应对比

    Figure  5.  Comparison of coupled and uncoupled acceleration frequency response

    图  6  不同飞行速度下空气舵耦合加速度频域响应

    Figure  6.  Coupled acceleration frequency response of the air rudder at different flight velocities

    图  7  空气舵耦合加速度RMS值随飞行速度变化

    Figure  7.  Variation of coupled acceleration RMS of the air rudder against flight velocities

  • [1] ROBERTSON J E. Prediction of in flight fluctuating pressure environments including protuberance induced flow: N71-36677[R]. Wyle Laboratories-Research Staff Report WR71-10
    [2] 黄寿康. 流体力学弹道载荷环境[M]. 北京: 宇航出版社, 1991: 65-81
    [3] 龙万花, 陈伟芳, 宋松和. 旋成体跨音速脉动压力环境分析与预测[J]. 国防科技大学学报, 2004, 26(1): 17-20 doi:  10.3969/j.issn.1001-2486.2004.01.005

    LONG W H, CHEN W F, SONG S H. Analysis and prediction of the fluctuating pressure induced by rotated aircraft at transonic Mach numbers[J]. Journal of National University of Defense Technology, 2004, 26(1): 17-20 doi:  10.3969/j.issn.1001-2486.2004.01.005
    [4] 刘振皓, 任方, 王骁峰, 等. 旋成体模型仪器舱脉动压力空间相关特性研究[J]. 宇航学报, 2016, 37(12): 1299-1305

    LIU Z H, REN F, WANG X F, et al. Study of fluctuating pressure spatial correlation characteristics of instrument cabin model in a rotated body[J]. Journal of Astronautics, 2016, 37(12): 1299-1305
    [5] 陈予恕, 钟顺, 郭虎伦. 高超声速空间飞行器中若干非线性动力学问题进展[J]. 科学技术与工程, 2012, 12(25): 6401-6409 doi:  10.3969/j.issn.1671-1815.2012.25.029

    CHEN Y S, ZHONG S, GUO H L. Advances on dynamic problems of hypersonic space vehicle[J]. Science Technology and Engineering, 2012, 12(25): 6401-6409 doi:  10.3969/j.issn.1671-1815.2012.25.029
    [6] 辛健强, 屈强, 许小静, 等. 高超声速飞行器舵面热气动弹性分析方法研究[J]. 强度与环境, 2015, 42(3): 1-9 doi:  10.3969/j.issn.1006-3919.2015.03.001

    XIN J Q, QU Q, XU X J, et al. Aerothermoelastic analysis method of hypersonic vehicle rudder[J]. Structure & Environment Engineering, 2015, 42(3): 1-9 doi:  10.3969/j.issn.1006-3919.2015.03.001
    [7] 郭静, 吴振强, 张伟, 等. 热噪声复合环境下飞行器结构动响应预示技术研究进展[J]. 强度与环境, 2014, 41(6): 1-10 doi:  10.3969/j.issn.1006-3919.2014.06.001

    GUO J, WU Z Q, ZHANG W, et al. Review of study on prediction of dynamical response for aerocraft structure under a combination of thermal and acoustic environments[J]. Structure & Environment Engineering, 2014, 41(6): 1-10 doi:  10.3969/j.issn.1006-3919.2014.06.001
    [8] 刘时秀, 张卫红, 呙道军, 等. 飞行器宽带声振响应预示技术研究[J]. 导弹与航天运载技术, 2011(5): 50-53 doi:  10.3969/j.issn.1004-7182.2011.05.013

    LIU S X, ZHANG W H, GUO D J, et al. Study on wide-band vibro-acoustic response prediction for the aerocraft[J]. Missiles and Space Vehicles, 2011(5): 50-53 doi:  10.3969/j.issn.1004-7182.2011.05.013
    [9] MEI C. A finite element approach for nonlinear panel flutter[J]. AIAA Journal, 1977, 15(8): 1107-1110 doi:  10.2514/3.60760
    [10] 刘超峰, 李海东, 杨炳渊. 身‒翼‒舵组合体高超声速变攻角颤振计算[J]. 宇航学报, 2011, 32(8): 1663-1669 doi:  10.3873/j.issn.1000-1328.2011.08.002

    LIU C F, LI H D, YANG B Y. Flutter analysis of body-wing-rudder configurations at hypersonic velocities and variable attack-angles[J]. Journal of Astronautics, 2011, 32(8): 1663-1669 doi:  10.3873/j.issn.1000-1328.2011.08.002
    [11] 杨超, 许赟, 谢长川. 高超声速飞行器气动弹性力学研究综述[J]. 航空学报, 2010, 31(1): 1-11

    YANG C, XU Y, XIE C C. Review of studies on aeroelasticity of hypersonic vehicles[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(1): 1-11
    [12] 郭静, 张忠, 王建民, 等. 基于活塞理论的壁板非线性气动弹性分析[J]. 强度与环境, 2018, 45(5): 7-12

    GUO J, ZHANG Z, WANG J M, et al. Study on nonlinear aeroelastic response of panel structure based on classical and local piston theories[J]. Structure & Environment Engineering, 2018, 45(5): 7-12
    [13] 贺尔铭, 陈兵, 张忠, 等. 考虑后屈曲的2D C/SiC复合材料板热振动分析[J]. 强度与环境, 2016, 43(4): 49-56 doi:  10.3969/j.issn.1006-3919.2016.04.008

    HE E M, CHEN B, ZHANG Z, et al. Modal analysis of 2D C/SiC composite panels in thermal environment considering post buckling[J]. Structure & Environment Engineering, 2016, 43(4): 49-56 doi:  10.3969/j.issn.1006-3919.2016.04.008
  • [1] 朱世清, 连荫俊, 郝永波, 李博, 臧立彬.  大范围弱刚性运动模拟系统的振动特性分析与振动抑制方法 . 航天器环境工程, 2023, 40(5): 455-460. doi: 10.12126/see.2023093
    [2] 范为, 王缅, 邢夏宁, 姚延风.  基于最小曲率半径的软波导安全裕度校核方法 . 航天器环境工程, 2022, 39(2): 205-209. doi: 10.12126/see.2022.02.014
    [3] 刘明辉, 刘闯, 向树红, 韩晓健, 周春燕, 邹雍.  模拟自由边界对卫星动力学特性的影响 . 航天器环境工程, 2022, 39(3): 225-229. doi: 10.12126/see.2022.03.001
    [4] 李青, 逯运通, 赵会光.  多舱段航天器振动基频分配速算方法 . 航天器环境工程, 2020, 37(1): 25-30. doi: 10.12126/see.2020.01.004
    [5] 李健, 胡明, 孟少华, 张立建, 高兴文, 张丹丹.  扇形太阳电池阵折展机构模态及谐响应分析 . 航天器环境工程, 2020, 37(5): 435-439. doi: 10.12126/see.2020.05.003
    [6] 谢伟华, 尹家聪, 林勇文, 董锴, 陈曦, 李正举.  卫星星箭耦合力学分析模型二次缩聚方法 . 航天器环境工程, 2019, 36(4): 330-334. doi: 10.12126/see.2019.04.005
    [7] 陈涛, 郝建滨.  基于三维有限元模型的弹翼胶接结构胶层应力分析 . 航天器环境工程, 2018, 35(3): 228-232. doi: 10.12126/see.2018.03.005
    [8] 赵星, 张翔.  三单元立方体卫星的结构特性分析 . 航天器环境工程, 2018, 35(2): 165-169. doi: 10.12126/see.2018.02.011
    [9] 张群, 毕京丹, 王英诚, 卫国, 胡彦平.  小卫星声-振联合试验耦合响应分析及预估 . 航天器环境工程, 2017, 34(5): 495-499. doi: 10.12126/see.2017.05.007
    [10] 刘明辉, 王剑, 高海洋, 岳志勇, 张俊刚.  大型振动台结构设计中的有限元分析与应用 . 航天器环境工程, 2016, 33(3): 257-261. doi: 10.12126/see.2016.03.005
    [11] 秦玉灵1, 韩增尧2, 邹元杰1, 王建民3基于粒子群优化算法的整星有限元模型修正方法 . 航天器环境工程, 2015, 32(3): 243-246. doi: 10.12126/see.2015.03.003
    [12] 刘明辉, 王剑, 焦安超, 韩晓健.  一种卫星在轨自由边界条件模拟方法 . 航天器环境工程, 2014, 31(1): 19-24. doi: 10.12126/see.2014.01.004
    [13] 吕海炜, 李映辉, 李中华, 李亮.  超声速气流下黏弹性夹层壁板非线性颤振分析 . 航天器环境工程, 2013, 30(1): 40-48.
    [14] 陈洪, 刘谦, 周宏亮, 欧峰.  加速度-气动力综合试验系统研制 . 航天器环境工程, 2013, 30(5): 544-549.
    [15] 任青梅, 杨志斌, 成 竹, 蒋军亮.  气动加热与热响应耦合分析及试验研究 . 航天器环境工程, 2012, 29(3): 304-307.
    [16] 秦玉灵, 孔宪仁, 罗文波.  基于统计分析的有限元模型修正 . 航天器环境工程, 2011, 28(3): 239-245. doi: 10.12126/see.2011.03.006
    [17] 蒋华兵, 李春丽, 陈强洪.  再入飞行器脉动压力环境特性分析 . 航天器环境工程, 2010, 27(3): 378-382.
    [18] 韩 潇, 祁 妍.  大型真空容器结构设计中的有限元分析与应用 . 航天器环境工程, 2009, 26(2): 150-153.
    [19] 应用有限元法对某热试验支架进行热-结构耦合变形分析 . 航天器环境工程, 2005, 22(6): 327-329.
    [20] KM5A真空容器屈曲稳定有限元分析 . 航天器环境工程, 2004, 21(3): 23-28.
  • 加载中
图(7)
计量
  • PDF下载量:  35
  • 文章访问数:  278
  • HTML全文浏览量:  58
文章相关
  • 中图分类号:  TB115.7; V414.8
  • 收稿日期:  2019-08-28
  • 修回日期:  2020-05-12
  • 网络出版日期:  2020-06-17
  • 刊出日期:  2020-06-01

空气舵气动力‒脉动压力‒结构耦合响应分析

doi: 10.12126/see.2020.03.004
    基金项目:  国家自然科学基金项目(编号:11502024;11502023;11402028;51805036)
    作者简介:

    郭 静(1983—),女,博士学位,高级工程师,主要从事热气动弹性仿真与试验技术研究。E-mail:guojing2662632@126.com

  • 中图分类号: TB115.7; V414.8

摘要: 气动力、脉动压力、结构振动相互作用,组成复杂的多场耦合系统,给动力学分析带来极大的挑战。文章基于PCL和DMAP语言自主研发了气动力–脉动压力–结构耦合响应分析软件,以复合材料空气舵为研究对象,建立其有限元模型,并开展模态分析;进而,建立基于Van Dyke修正活塞理论的气动模型,基于模态法分析了气动力–脉动压力–结构三者耦合的空气舵响应,并与不考虑气动力效应的非耦合结构响应进行对比,探究了气动力耦合效应对空气舵响应的影响规律。结果表明,气动弹性效应能使得空气舵振动响应从随机振动变为发散极限环振荡形式的高阶运动,显著改变脉动压力响应谱。可以预测,结构声疲劳分析中必须考虑气动弹性效应。

English Abstract

郭静, 张忠, 谷迎松, 等. 空气舵气动力‒脉动压力‒结构耦合响应分析[J]. 航天器环境工程, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004
引用本文: 郭静, 张忠, 谷迎松, 等. 空气舵气动力‒脉动压力‒结构耦合响应分析[J]. 航天器环境工程, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004
GUO J, ZHANG Z, GU Y S, et al. Analysis of aerodynamic force- fluctuating pressure- structure coupled response for air rudder[J]. Spacecraft Environment Engineering, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004
Citation: GUO J, ZHANG Z, GU Y S, et al. Analysis of aerodynamic force- fluctuating pressure- structure coupled response for air rudder[J]. Spacecraft Environment Engineering, 2020, 37(3): 235-239 doi:  10.12126/see.2020.03.004
    • 航天飞行器在大气中飞行时,气流与结构之间的相互作用会导致其表面气流出现复杂的流动形式,包括湍流、分离和激波振荡等[1],从而使得飞行器遭受严酷的脉动压力载荷,这种脉动压力亦被称为“气动噪声”[2]。脉动压力频带较宽,通常覆盖了低频、中频和高频,总声压级为130~170 dB[2-3]。强烈的脉动压力载荷对飞行器结构和内部仪器设备而言构成恶劣的振动与噪声环境,造成结构的疲劳损坏或仪器设备失效甚至任务失败。因此,预测脉动压力引起的飞行器结构振动响应,对于飞行器载荷环境设计及保证飞行可靠性具有重要意义[4]

      另一方面,航天飞行器以高马赫数飞行时将承受严酷的气动力载荷。气动力与结构惯性力、弹性力相互耦合,引发了气动弹性问题[5-6]。严重的耦合效应将可能导致飞行器结构或全机发生颤振,造成飞行器性能下降甚至结构破坏,危及飞行安全。

      气动力、脉动压力、结构振动相互作用,互为因果,组成复杂的多场耦合系统,给动力学分析带来极大挑战。目前,国内外现有的耦合响应分析方法主要包括两类:一类是不考虑气动力与结构的耦合作用,仅考虑脉动压力作用下结构的动响应,即为声振响应预示问题[7-8];另一类是不考虑脉动压力的作用,仅考虑气动力与结构的耦合效应,即为气动弹性问题[9-11]。有关气动力–脉动压力–结构多场耦合环境下的动力学响应分析研究相对较少。

      本文基于PCL和DMAP语言自主研发了气动力–脉动压力–结构耦合响应分析软件,以复合材料空气舵为研究对象,建立其有限元模型,并开展模态分析。进而建立基于Van Dyke修正活塞理论的气动模型,基于模态法分析气动力–脉动压力–结构三者耦合的空气舵响应,并与不考虑气动力耦合效应的脉动压力作用下空气舵结构的非耦合响应进行对比分析,以探究气动力耦合效应对空气舵响应的影响规律。

    • 空气舵有限元模型如图1所示。复合材料空气舵的蒙皮和腹板的材料为石墨–环氧树脂,其材料力学性能和层合板的铺叠方式参见文献[12],其厚度均为1 mm,采用壳单元模拟;筋条材料为钛合金,用以加强蒙皮和翼肋,采用梁单元模拟;筋条交叉骨架点上布置集中质量单元,每个为1 kg,用以模拟结构的质量分布;边界条件为根部部分结点简支[13]。分析获得了该模型前10阶固有频率和模态振型,如图2所示。

      图  1  复合材料空气舵有限元模型

      Figure 1.  FEM of air rudder made of composite materials

      图  2  复合材料空气舵固有频率和模态振型

      Figure 2.  Natural frequencies and modal shapes of composite air rudder

    • 矩阵形式表示的气动力–脉动压力–结构耦合振动方程为

      $${{M}}\ddot {{u}} + {{C}}\dot {{u}} + {{Ku}} = {{f}} + {{A}}\left( {{{u}}\left( t \right)} \right),$$

      式中:M为质量矩阵;C为阻尼矩阵;K为刚度矩阵;A(u)为气动力向量[12]f为脉动压力。A(u)由Van Dyke修正活塞理论求得。若同时考虑式(1)右端两项,则为气动力–脉动压力–结构三者耦合的空气舵响应。若仅考虑式(1)右端第一项,则为不考虑气动力耦合效应的脉动压力作用下空气舵结构的非耦合响应。

      本文充分利用MSC.Nastran软件中基于模态法的瞬态响应求解模块SOL112,采用DMAP语言编写程序,进行二次开发,增加基于活塞理论的高超声速非定常气动力计算、脉动压力施加、考虑气动力耦合效应的耦合响应分析模块与不考虑气动力耦合效应的非耦合响应分析模块,开发了气动力–脉动压力–结构耦合响应分析软件,实现结构耦合与非耦合响应分析,并采用PCL语言在Patran软件中设计了嵌入式前后置处理界面。气动力–脉动压力–结构耦合响应分析软件界面如图3所示。

      图  3  气动力–脉动压力–结构耦合响应分析软件界面

      Figure 3.  Interface of software for analyzing aerodynamic force- fluctuating pressure- structure coupled response

    • 上述空气舵的脉动压力时域曲线如图4所示,为频域范围50~500 Hz的平直谱,采用开发的耦合响应分析软件,基于模态法和活塞理论,计算了气动力–脉动压力–结构三者耦合的空气舵响应和不考虑气动力耦合效应的脉动压力作用下空气舵结构的非耦合响应。

      图  4  脉动压力时域曲线

      Figure 4.  Time domain curve of the fluctuating pressure

      针对马赫数2.0、空气密度1.225 kg/m3、飞行速度2400 m/s的分析工况,本文主要关心舵面整体的动态响应,结合响应求解的频率范围,选取了前10阶模态,积分步长为0.000 1 s,模态阻尼比为0.01,计算得到了空气舵耦合与非耦合加速度频域响应,对比结果如图5所示。由图可知,空气舵耦合与非耦合加速度响应的RMS值分别为1.48g和1.61g。研究表明,在一定飞行速度下,气动力在响应分析中表现为阻尼效应,使得考虑气动力耦合效应的空气舵结构响应比非耦合的响应偏小。

      图  5  耦合与非耦合加速度频域响应对比

      Figure 5.  Comparison of coupled and uncoupled acceleration frequency response

      为了探究气动力耦合效应对空气舵响应的影响规律,进一步计算得到了飞行速度分别为2400、2700、3000、3300、3450和3550 m/s时空气舵耦合加速度频域响应以及加速度RMS值随飞行速度的变化曲线,分别如图6图7所示。

      图  6  不同飞行速度下空气舵耦合加速度频域响应

      Figure 6.  Coupled acceleration frequency response of the air rudder at different flight velocities

      图  7  空气舵耦合加速度RMS值随飞行速度变化

      Figure 7.  Variation of coupled acceleration RMS of the air rudder against flight velocities

      分析表明:空气舵耦合响应随着飞行速度的增加不断增大;与非耦合响应相比,飞行速度低时,由于远离颤振速度,舵面气动弹性系统的阻尼较大,空气舵耦合响应主要为脉动压力产生的强迫振动,气动力表现为阻尼效应,使得耦合响应比非耦合响应偏小;飞行速度增大到某个值时,耦合响应与非耦合响应大小相等;随着飞行速度继续增大,接近颤振速度时,系统接近失稳,空气舵耦合响应主要变为自激振动,导致响应迅速增大,直至颤振发生,临界颤振速度为3600 m/s。

    • 本文自主研发了气动力–脉动压力–结构耦合响应分析软件,分析了复合材料空气舵气动力–脉动压力–结构耦合响应,并与不考虑气动力耦合效应的脉动压力作用下空气舵结构的非耦合响应进行对比,探究了气动力耦合效应对空气舵响应的影响规律。研究结果表明:气动力耦合效应能够显著影响脉动压力作用下结构的响应特性;飞行速度低时,耦合响应主要为脉动压力产生的强迫振动,气动力在响应分析中表现为阻尼效应;随着飞行速度继续增大,耦合响应特性由强迫振动变为自激振动,发生极限环振荡,响应发散。

      本文研究为航天飞行器典型结构的合理设计和多场耦合响应分析提供了技术支撑。气动弹性效应能够改变结构耦合响应特性的机理有待深入研究。

参考文献 (13)

目录

    /

    返回文章
    返回