留言板

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

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

抑制柱体结构涡激振动的非线性能量阱减振装置优化设计

陈东阳 顾超杰 朱卫军 李迺璐 杨俊伟 芮筱亭

陈东阳, 顾超杰, 朱卫军, 李迺璐, 杨俊伟, 芮筱亭. 抑制柱体结构涡激振动的非线性能量阱减振装置优化设计[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.10.0613
引用本文: 陈东阳, 顾超杰, 朱卫军, 李迺璐, 杨俊伟, 芮筱亭. 抑制柱体结构涡激振动的非线性能量阱减振装置优化设计[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.10.0613
Dong-yang CHEN, Chao-jie GU, Wei-jun ZHU, Nai-lu LI, Jun-wei YANG, Xiao-ting RUI. AN OPTIMUM DESIGN OF NONLINEAR ENERGY SINK VIBRATION ABSORBING DEVICE FOR SUPPRESSING VORTEX-INDUCED VIBRATION OF CYLINDRICAL STRUCTURES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.10.0613
Citation: Dong-yang CHEN, Chao-jie GU, Wei-jun ZHU, Nai-lu LI, Jun-wei YANG, Xiao-ting RUI. AN OPTIMUM DESIGN OF NONLINEAR ENERGY SINK VIBRATION ABSORBING DEVICE FOR SUPPRESSING VORTEX-INDUCED VIBRATION OF CYLINDRICAL STRUCTURES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.10.0613

抑制柱体结构涡激振动的非线性能量阱减振装置优化设计

doi: 10.6052/j.issn.1000-4750.2019.10.0613
基金项目: 江苏省自然科学基金项目(BK20190871);国家自然科学基金项目(11672261);江苏省研究生科研与实践创新计划项目(KYCX19_2104)
详细信息
    作者简介:

    顾超杰(1997−),男,无锡市人,硕士生,主要从事流固耦合动力学研究(E-mail: 727713838@qq.com)

    朱卫军(1976−),男,南通市人,教授,博士,博导,主要从事空气动力学研究(E-mail: wjzhu@yzu.edu.cn)

    李迺璐(1986−),女,扬州市人,副教授,博士,硕导,主要从事风力机叶片动力学与控制研究(E-mail: nlli@yzu.edu.cn)

    杨俊伟(1990−),男,扬州市人,讲师,博士生,主要从事空气动力学研究(E-mail: yangjunwei@yzu.edu.cn)

    芮筱亭(1956−),男,镇江市人,中国科学院院士,博士,南京理工大学发射动力学研究所所长,主要从事多体系统动力学研究(E-mail: ruixt@163.net)

    通讯作者: 陈东阳(1987−),男,连云港市人,讲师,博士,硕导,主要从事流固耦合动力学研究(E-mail: cdy_1988@sina.cn)
  • 中图分类号: O353.1

AN OPTIMUM DESIGN OF NONLINEAR ENERGY SINK VIBRATION ABSORBING DEVICE FOR SUPPRESSING VORTEX-INDUCED VIBRATION OF CYLINDRICAL STRUCTURES

  • 摘要: 为了抑制柱体结构的涡激振动(Vortex induced vibration, VIV)问题,该文基于Van der Pol尾流振子模型、结构动力学理论和优化算法,建立了非线性能量阱(Nonlinear energy sink, NES)减振装置优化设计仿真模型。采用Van der Pol尾流振子模型,将其与二维弹性支撑柱体结构动力学方程结合,建立二维弹性支撑柱体单自由度涡激振动预报模型,并与文献实验数据对比,验证了模型的准确性;在模型中嵌入NES控制模型,完成NES作用下的单自由度柱体VIV预测模型,并将其与优化算法相结合,通过优化算法改变非线性能量阱的参数,设计出对柱体涡激振动抑制效果满足设计要求的非线性能量阱;通过仿真分析优化后的NES结构参数对柱体涡激振动的抑制效果和机理。案例优化结果表明,在$U_{\rm{r}} = 5.5$时,NES作用下的柱体振幅减小了66.39%,减振效果明显。同时,该文采用优化方法获得的NES结构参数符合实际生产要求,确保了实物生产设计的可能性。
  • 图  1  物理模型示意图

    Figure  1.  Physical model

    图  2  优化流程图

    Figure  2.  Optimized flow chart

    图  3  基于Van der Pol尾流振子模型的VIV计算结果

    Figure  3.  The VIV simulation results based on Van der Pol wake oscillator model

    图  4  设计速度下的计算结果

    Figure  4.  Calculation results at design speed

    图  5  振动响应对比图

    Figure  5.  Comparisons of vibration response

    图  6  柱体涡激振动二维相图

    Figure  6.  Two-dimensional phase diagram of VIV of the cylinder

  • [1] Zhang W, Li X, Ye Z, et al. Mechanism of frequency lock-in in vortex-induced vibrations at low Reynolds numbers [J]. Journal of Fluid Mechanics, 2015, 783: 72 − 102. doi:  10.1017/jfm.2015.548
    [2] Khalak A, Williamson C H K. Fluid forces and dynamics of a hydroelastic structure with very low mass and damping [J]. Journal of Fluids and Structures, 1997, 11(8): 973 − 982. doi:  10.1006/jfls.1997.0110
    [3] Khan N B, Ibrahim Z, Nguyen L T T, et al. Numerical investigation of the vortex-induced vibration of an elastically mounted circular cylinder at high Reynolds number (Re= 104) and low mass ratio using the RANS code [J]. PloS one, 2017, 12(10): e0185832.
    [4] 唐友刚, 邵卫东, 张杰, 等. 深海顶张力立管参激-涡激耦合振动响应分析[J]. 工程力学, 2013, 30(5): 282 − 286.

    Tang Yougang, Shao Weidong, Zhang Jie, et al. Dynamic response analysis for coupled parametric vibration and vortex-induced vibration of top-tensioned riser in deep-sea [J]. Engineering Mechanics, 2013, 30(5): 282 − 286. (in Chinese)
    [5] 周佳豪, 马文勇, 黄伯城. 临界雷诺数区光滑圆柱振动与气动力研究[J]. 工程力学, 2019, 36(增刊1): 306 − 310.

    Zhou Jiahao, Ma Wenyong, Huang Bocheng. Study on vibration and aerodynamic forces of a smooth cylinder in the critical Reynolds number range [J]. Engineering Mechanics, 2019, 36(Suppl1): 306 − 310. (in Chinese)
    [6] 喻梅. 大跨度桥梁颤振及涡激振动主动控制 [D]. 四川: 西南交通大学, 2013.

    Yu Mei. Active control for flutter and vortex-induced vibration of long-span bridges [D]. Sichuan: Southwest Jiaotong University, 2013. (in Chinese)
    [7] Liu M, Zhao M, Lu L, et al. Numerical Investigation of Vortex-Induced Vibration of Circular Cylinder with Multiple Control Rods at Low Reynolds Number [C]//Symposium on Fluid-Structure-Sound Interactions and Control. Springer, Singapore, 2017: 217 − 221.
    [8] 徐万海, 栾英森, 余杨, 覃雯琪. 大倾角附螺旋列板倾斜圆柱涡激振动抑制分析[J]. 工程力学, 2018, 35(1): 236 − 245.

    Xu Wanhai, Luan Yingsen, Yu Yang, Tan Wenqi. Vortex-induced vibration suppression for an inclined flexible cylinder at large inclination angle [J]. Engineering Mechanics, 2018, 35(1): 236 − 245. (in Chinese)
    [9] Bichiou Y, Hajj M R, Nayfeh A H. Effectiveness of a nonlinear energy sink in the control of an aeroelastic system [J]. Nonlinear Dynamics, 2016, 86(4): 2161 − 2177. doi:  10.1007/s11071-016-2922-y
    [10] Ahmadabadi Z N. Nonlinear energy transfer from an engine crankshaft to an essentially nonlinear attachment [J]. Journal of Sound and Vibration, 2019, 443: 139 − 154. doi:  10.1016/j.jsv.2018.11.040
    [11] Blanchard A, Bergman L A, Vakakis A F. Vortex-induced vibration of a linearly sprung cylinder with an internal rotational nonlinear energy sink in turbulent flow [J]. Nonlinear Dynamics, 2020, 99(1): 593 − 609.
    [12] Blanchard A B, Gendelman O V, Bergman L A, et al. Capture into slow invariant manifold in the fluid-structure dynamics of a sprung cylinder with a nonlinear rotator [J]. Journal of Fluids & Structures, 2016, 63: 155 − 173.
    [13] Blanchard A, Bergman L A, Vakakis A F. Targeted energy transfer in laminar vortex-induced vibration of a sprung cylinder with a nonlinear dissipative rotator [J]. Physica D Nonlinear Phenomena, 2017, 350: 26 − 44. doi:  10.1016/j.physd.2017.03.003
    [14] Tumkur R K R, Calderer R, Masud A, et al. Computational study of vortex-induced vibration of a sprung rigid circular cylinder with a strongly nonlinear internal attachment [J]. Journal of Fluids and Structures, 2013, 40: 214 − 232. doi:  10.1016/j.jfluidstructs.2013.03.008
    [15] Mehmood A, Nayfeh A H, Hajj M R. Effects of a non-linear energy sink (NES) on vortex-induced vibrations of a circular cylinder [J]. Nonlinear Dynamics, 2014, 77(3): 667 − 680. doi:  10.1007/s11071-014-1329-x
    [16] Dai H L, Abdelkefi A, Wang L. Vortex-induced vibrations mitigation through a nonlinear energy sink [J]. Communications in Nonlinear Science and Numerical Simulation, 2017, 42: 22 − 36. doi:  10.1016/j.cnsns.2016.05.014
    [17] Chen D Y, Abbas L K, Wang G P, Rui X T, Marzocca P. Numerical study of flow-induced vibrations of cylinders under the action of nonlinear energy sinks (NESs) [J]. Nonlinear Dynamics, 2018, 94(2): 925 − 957.
    [18] Tumkur R K R, Domany E, Gendelman O V, et al. Reduced-order model for laminar vortex-induced vibration of a rigid circular cylinder with an internal nonlinear absorber [J]. Communications in Nonlinear Science and Numerical Simulation, 2013, 18(7): 1916 − 1930. doi:  10.1016/j.cnsns.2012.11.028
    [19] Facchinetti M L, Langre E D, Biolley F. Vortex-induced travelling waves along a cable [J]. European Journal of Mechanics, B/Fluids, 2004, 23(1): 199 − 208. doi:  10.1016/j.euromechflu.2003.04.004
    [20] Facchinetti M L, De Langre E, Biolley F. Coupling of structure and wake oscillators in vortex-induced vibrations [J]. Journal of Fluids and Structures, 2004, 19(2): 123 − 140. doi:  10.1016/j.jfluidstructs.2003.12.004
    [21] 葛继科, 邱玉辉, 吴春明, 蒲国林. 遗传算法研究综述[J]. 计算机应用研究, 2008(10): 2911 − 2916. doi:  10.3969/j.issn.1001-3695.2008.10.008

    Ge Jike, Qiu Yuhui, Wu Chunming, Pu Guolin. Summary of genetic algorithms research [J]. Application Research of Computers, 2008(10): 2911 − 2916. (in Chinese) doi:  10.3969/j.issn.1001-3695.2008.10.008
    [22] Holland J H. Concerning efficient adaptive systems [M]// Yovits M C, eds. Self-Organizing Systems, 1962: 215 − 230.
    [23] Holland J H. Adaptation in natural and artificial systems [M]. Cambridge, MA: MIT Press, 1992.
    [24] Bagley J D. The behavior of adaptive systems which employ genetic and correlation algorithms : technical report [R]. University of Michigan, United States, 1967.
    [25] Stappenbelt B, Lalji F, Tan G. Low mass ratio vortex-induced motion [C]. The 16th Australasian Fluid Mechanics Conference, Australasian, 2007: 1491 − 1497.
    [26] 陈东阳, 肖清, 顾超杰, 王国平, 芮筱亭. 柱体结构涡激振动数值计算[J]. 振动与冲击, 2020 (录用)

    Chen Dongyang, Xiao Qing, Gu Chaojie, Wang Guoping, Rui Xiaoting. Numerical calculation of vortex-induced vibration of a cylinder structure[J]. Journal of Vibration and Shock, 2020 (录用). (in Chinese)
  • [1] 颜欣桐, 徐龙河.  基于遗传算法的钢筋混凝土框架-剪力墙结构失效模式多目标优化 . 工程力学, doi: 10.6052/j.issn.1000-4750.2016.12.0981
    [2] 蒲利东, 罗务揆, 严泽洲.  气动伺服弹性系统结构陷幅滤波器优化设计 . 工程力学, doi: 10.6052/j.issn.1000-4750.2016.12.0971
    [3] 许坤, 葛耀君.  基于尾流振子模型的桥梁节段至实桥涡激共振振幅转换关系 . 工程力学, doi: 10.6052/j.issn.1000-4750.2015.07.0588
    [4] 陈鑫, 李爱群, 徐庆阳, 张志强.  基于遗传算法的高耸结构环形TLD满意优化设计 . 工程力学, doi: 10.6052/j.issn.1000-4750.2014.10.0918
    [5] 喻梅, 廖海黎, 李明水, 马存明, 贾宏宇, 罗楠.  基于经验线性涡激力的桥梁涡激振动TMD控制 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.08.0630
    [6] 秦伟, 康庄, 宋儒鑫.  基于离散点涡的双自由度涡激振动尾流振子模型研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2011.05.0264
    [7] 薛晓敏, 孙 清, 张 陵, 伍晓红.  基于遗传算法策略的含时滞结构振动主动控制研究 . 工程力学,
    [8] 陈伟民, 张立武, 李 敏.  采用改进尾流振子模型的柔性海洋立管的涡激振动响应分析 . 工程力学,
    [9] 刘书田, 林哲祺.  考虑作动器联接方式的结构形状控制优化 . 工程力学,
    [10] 朱南海, 赵晓华.  基于遗传算法的Tikhonov正则参数优化计算 . 工程力学,
    [11] 黄冀卓, 王 湛.  基于遗传算法的离散型结构拓扑优化设计 . 工程力学,
    [12] 燕乐纬, 蹇开林, 黄晓刚.  基于广义遗传算法的结构动力响应优化 . 工程力学,
    [13] 张连振, 黄 侨, 王潮海.  基于多目标遗传算法的传感器优化布点研究 . 工程力学,
    [14] 李守巨, 刘迎曦.  改进遗传算法在非线性热传导参数识别中的应用 . 工程力学,
    [15] 黄维平, 刘娟, 李华军.  基于遗传算法的传感器优化配置 . 工程力学,
    [16] 修英姝, 崔德刚.  复合材料层合板稳定性的铺层优化设计 . 工程力学,
    [17] 郭惠勇, 蒋健, 张陵.  改进遗传算法在MRFD半主动控制系统优化配置中的应用 . 工程力学,
    [18] 郭鹏飞, 韩英仕.  离散变量结构优化设计的拟满应力遗传算法 . 工程力学,
    [19] 钟登华, 梅传书, 韩圣章.  基于复合形算子的基础支护桩优化设计智能算法研究 . 工程力学,
    [20] 谢楠, 陈英俊.  遗传算法的改进策略及其在桥梁抗震优化设计中的应用效果 . 工程力学,
  • 加载中
计量
  • 文章访问数:  90
  • HTML全文浏览量:  37
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-27
  • 修回日期:  2020-03-19
  • 网络出版日期:  2020-06-02

抑制柱体结构涡激振动的非线性能量阱减振装置优化设计

doi: 10.6052/j.issn.1000-4750.2019.10.0613
    基金项目:  江苏省自然科学基金项目(BK20190871);国家自然科学基金项目(11672261);江苏省研究生科研与实践创新计划项目(KYCX19_2104)
    作者简介:

    顾超杰(1997−),男,无锡市人,硕士生,主要从事流固耦合动力学研究(E-mail: 727713838@qq.com)

    朱卫军(1976−),男,南通市人,教授,博士,博导,主要从事空气动力学研究(E-mail: wjzhu@yzu.edu.cn)

    李迺璐(1986−),女,扬州市人,副教授,博士,硕导,主要从事风力机叶片动力学与控制研究(E-mail: nlli@yzu.edu.cn)

    杨俊伟(1990−),男,扬州市人,讲师,博士生,主要从事空气动力学研究(E-mail: yangjunwei@yzu.edu.cn)

    芮筱亭(1956−),男,镇江市人,中国科学院院士,博士,南京理工大学发射动力学研究所所长,主要从事多体系统动力学研究(E-mail: ruixt@163.net)

    通讯作者: 陈东阳(1987−),男,连云港市人,讲师,博士,硕导,主要从事流固耦合动力学研究(E-mail: cdy_1988@sina.cn)
  • 中图分类号: O353.1

摘要: 为了抑制柱体结构的涡激振动(Vortex induced vibration, VIV)问题,该文基于Van der Pol尾流振子模型、结构动力学理论和优化算法,建立了非线性能量阱(Nonlinear energy sink, NES)减振装置优化设计仿真模型。采用Van der Pol尾流振子模型,将其与二维弹性支撑柱体结构动力学方程结合,建立二维弹性支撑柱体单自由度涡激振动预报模型,并与文献实验数据对比,验证了模型的准确性;在模型中嵌入NES控制模型,完成NES作用下的单自由度柱体VIV预测模型,并将其与优化算法相结合,通过优化算法改变非线性能量阱的参数,设计出对柱体涡激振动抑制效果满足设计要求的非线性能量阱;通过仿真分析优化后的NES结构参数对柱体涡激振动的抑制效果和机理。案例优化结果表明,在$U_{\rm{r}} = 5.5$时,NES作用下的柱体振幅减小了66.39%,减振效果明显。同时,该文采用优化方法获得的NES结构参数符合实际生产要求,确保了实物生产设计的可能性。

English Abstract

陈东阳, 顾超杰, 朱卫军, 李迺璐, 杨俊伟, 芮筱亭. 抑制柱体结构涡激振动的非线性能量阱减振装置优化设计[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.10.0613
引用本文: 陈东阳, 顾超杰, 朱卫军, 李迺璐, 杨俊伟, 芮筱亭. 抑制柱体结构涡激振动的非线性能量阱减振装置优化设计[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.10.0613
Dong-yang CHEN, Chao-jie GU, Wei-jun ZHU, Nai-lu LI, Jun-wei YANG, Xiao-ting RUI. AN OPTIMUM DESIGN OF NONLINEAR ENERGY SINK VIBRATION ABSORBING DEVICE FOR SUPPRESSING VORTEX-INDUCED VIBRATION OF CYLINDRICAL STRUCTURES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.10.0613
Citation: Dong-yang CHEN, Chao-jie GU, Wei-jun ZHU, Nai-lu LI, Jun-wei YANG, Xiao-ting RUI. AN OPTIMUM DESIGN OF NONLINEAR ENERGY SINK VIBRATION ABSORBING DEVICE FOR SUPPRESSING VORTEX-INDUCED VIBRATION OF CYLINDRICAL STRUCTURES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.10.0613
  • 结构流致振动是有着广泛工程应用背景的流固耦合问题,在航空工程、海洋工程等领域,当流体流过钝体时,钝体下游流场持续产生和脱落旋涡,从而导致结构受到周期变化的流体力作用,使结构以一定的频率和振幅振动,即为涡激振动(vortex-induced vibration, VIV)[1-5]。在流体作用下,当结构振动频率接近结构固有频率时,结构将产生涡激共振,此时振幅将远大于正常情况,并且共振时会发生“频率”锁定现象,在一定来流速度内,结构振动频率都不会发生改变。涡激共振的发生,结构横向上将产生大幅度的振动,受到周期性的疲劳应力作用,从而可能导致钝体结构的疲劳失效,甚至造成对结构的破坏。若结构长期处于频率锁定状态,将大大减少结构的疲劳寿命。

    为了减弱涡激振动对高柔性柱体结构疲劳寿命的影响,通常采用主动控制和被动控制来抑制其涡激振动。主动控制[6]通过输入外部能量来扰动流场,如采用声波激励、敲击振动等方式,将旋涡脱落控制在合适情况下,但是该类方法技术含量高,成本也更昂贵。而被动控制通常采用外加控制柱[7]、螺旋列板[8]等扰流装置扰乱流场结构或干扰旋涡的形成与泄放。但这些扰流装置往往会使阻力增大,并且还会引发其他形式的振动。非线性能量阱(nonlinear energy sink, NES)是能实现定向能量传递的具有立方非线性的吸振器[9-12]。NES通过与柱体共振,将来自柱体系统的振动能量定向传递到吸振器的振子上,并通过阻尼消耗能量,从而达到减小涡激目的[13-14]。NES不仅具有宽频吸振特性,而且该装置内置于柱体内部,不改变柱体的形状,不会增加额外阻力,成本低。

    采用NES来抑制涡激振动时,不同参数的NES对柱体振动的抑制效果不同。文献[15-17]通常对比不同NES参数下的涡激振动响应,寻找振动响应与参数间的规律。Mehmood等[15]使用计算流体力学(computational fluid dynamics, CFD)方法耦合结构动力学建立预测模型,研究NES的质量比对低雷诺数柱体结构涡激振动的抑制效果。Dai等[16]采用尾流振子模型研究了NES作用下的柱体结构频率耦合机理。本人也在前期工作中[17]研究不同NES参数对中等雷诺数柱体结构涡激振动的减振效果。但在实际设计中,人为选取NES的参数,通常设计出的NES往往并不能达到最优的抑制效果,同时效率较低。因此本文为了避免漫无目的地选取NES参数,建立了用于柱体结构涡激振动抑制的NES减振装置优化设计仿真模型,使用优化算法快速设计出合适的NES来有效抑制柱体结构的涡激振动,为实物设计提供参考设计目标。

    • 物理模型如图1所示:图1(a)中为二维平面内包含NES的单自由度柱体模型,把柱体结构看成质量-弹簧-尼系统,嵌入的NES也同样看作质量块,和柱体横向上有阻尼器和立方非线性弹簧连接;图1(b)为三维柱体涡激振动模型,来流经过柱体后,对柱体横向上产生振动。

      图  1  物理模型示意图

      Figure 1.  Physical model

      根据该物理模型可得到NES作用下单自由度柱体运动的控制方程[18]为:

      $$\begin{split}& (m - m_{\rm {nes}})\ddot y_1 + c\dot y_1 + ky_1 + c_{\rm {nes}}(\dot y_1 - \dot y_2) + \\& \qquad k_{\rm {nes}}{(y_1 - y_2)^3} = F_L(t) \end{split} $$ (1)
      $$m_{\rm {nes}}\ddot y_2 + c_{\rm {nes}}(\dot y_2 - \dot y_1) + k_{\rm {nes}}{(y_2 - y_1)^3} = 0$$ (2)

      式中:$m$为圆柱体的质量;$c$为结构阻尼系数;$k$为结构刚度系数;$y_1$$\dot y_1$$\ddot y_1$分别为柱体的横向位移、速度、加速度;$y_2$$\dot y_2$$\ddot y_2$分别为NES的横向位移、速度、加速度;$F_L$为作用在柱体横向上的流体力;$m_{\rm {nes}}$$c_{\rm {nes}}$$k_{\rm {nes}}$分别是NES的质量、阻比系数和刚度系数。

      涡激升力$F_L(t)$的半经验计算公式[16]如下:

      $$F_L(t) = - \frac{1}{2}C_{\rm D}\rho _{\rm f}DU\dot y_1 + \frac{1}{4}C_{\rm L0}\rho _{\rm f}D{U^2}q(t)$$ (3)

      式中:$C_{\rm D}$为阻力系数;$C_{\rm L0}$为稳定横向升力系数,都由实验和经验确定;$\rho $为流体密度;D为柱体特征长度;U为来流风速;$q(t)$为折减升力系数,它描述了邻近尾流的状况,并且能用如下Van der Pol等式[19-20]表达:

      $$\ddot q + \lambda \omega _{\rm s}({q^2} - 1)\dot q + \omega {_{\rm s}^2}q = \frac{P}{D}\ddot y_1$$ (4)

      式中:${\textit{ω}} _{\rm s} = 2\pi \cdot St \cdot U/D$为涡脱频率;$St$为Strouhal数;$\lambda $$P$等系数也同样为由实验和经验确定;

      根据式(1)~式(4)联立方程,建立柱体VIV模型:

      $$\begin{split}& (1 - \beta )\ddot y_1 + 2\varsigma \omega _0\dot y_1 + \omega {_0^2}y_1 + 2\varsigma _{\rm {nes}}\omega _0(\dot y_1 - \dot y_2) + \\&\qquad v\frac{{\omega {_0^2}}}{{{D^2}}}{(y1 - y2)^3} = - \frac{1}{{2m}}C_{\rm D}\rho _{\rm f}DU\dot y_1 +\\&\qquad \frac{1}{{4m}}C_{\rm L0}\rho _{\rm f}D{U^2}q(t) \end{split} \quad$$ (5)
      $$\ddot q + \lambda \omega _{\rm s}({q^2} - 1)\dot q + \omega {_{\rm s}^2}q = \frac{P}{D}\ddot y_1\qquad\qquad\qquad\quad$$ (6)
      $$\beta \ddot y_2 + 2\varsigma _{\rm {nes}}\omega _0(\dot y_2 - \dot y_1) + v\frac{{\omega {_0^2}}}{{{D^2}}}{(y_2 - y_1)^3} = 0\qquad$$ (7)

      式中:$k$为柱体的结构刚度;$k_{\rm {nes}}$为NES的结构刚度;$c$为柱体的结构阻尼;$c_{\rm {nes}}$为NES的结构阻尼;$\omega _0$为柱体的固有频率;$\varsigma $为柱体的阻尼比;$\varsigma _{\rm {nes}}$为NES的阻尼比;$\beta = {{m_{\rm {nes}}} / {m_{\rm cylinder}}}$$\gamma = ({{k_{\rm {nes}} \cdot {D^2})} / }$$k$$\xi = {{\varsigma _{\rm {nes}}} / \varsigma }$$\beta $$\gamma $$\xi $分别为NES与柱体结构的无量纲质量之比、刚度之比和阻尼之比。

      通过龙格库塔法求解该方程组,可得到柱体的振动响应。求解时,初始值给一个小扰动,设初始值$y_1 = 0.000001$$y_2 = \dot y_1 = \dot y_2 = q = \dot q = 0$,时间步长取$0.001\;{\rm{ s}}$

    • 使用优化算法与VIV仿真模型相结合,以NES对柱体涡激振动的抑制要求作为优化目标,对NES参数进行优化设计,建立NES减振装置优化设计仿真模型。

    • 优化模块中,以NES的三个参数(无量纲质量之比$\beta $、无量纲阻尼之比$\xi $和无量纲刚度之比$\gamma $)为设计参数。这三个参数是设计NES的重要参数,它们决定了NES对柱体涡激振动的抑制效果,也是优化过程中的设计变量。

      涡激共振的发生,结构横向上将产生大幅度的振动,受到周期性的疲劳应力作用,从而可能导致钝体结构的疲劳失效,甚至造成对结构的破坏。若结构长期处于频率锁定状态,将大大减少结构的疲劳寿命。为了减小涡激振动对柱体结构的疲劳损伤,通常设计减振装置以限制设计来流速度范围内的最大振幅。本文将VIV预测模型计算出的设计工况下柱体稳定振动时的最大横向振幅为优化目标函数,使柱体横向振幅在一定折减速度内产生的最大值控制在最小值(即设定风速段内各振幅值均方差最小),达到限制涡激振动的目的。

      考虑到NES的工程实际应用,本优化设计中对NES的无量纲质量之比、阻尼之比和刚度之比进行约束。工程实际中,需要NES质量越小越好:一方面减小柱体结构的承重;另一方面便于安装在柱体结构内部。所以对NES质量之比做出约束,令$0 < \beta < 0.3$。同理,工程上也难以生产出过高阻尼和过高或者过柔刚度的NES,所以对NES的无量纲阻尼之比和刚度之比做出限制:$0 < \xi < 2$$0 < \gamma < 3$。这样使NES在有效抑制横向振动的同时,也符合实际情况。为防止该风速段内某一风速下出现极大振幅造成结构损坏,对该风速段内出现的最大振幅进行限制,使NES作用下的柱体在设定风速段内出现的最大无量纲振幅都在0.2以下,即$\left(\dfrac{{y_1}}{D}\right)\max < 0.2$。同时为了防止放置在柱体结构内部的NES振动与柱体内壁发生碰撞,限制$\dfrac{{\left| {{y_1}{\rm{ - }}{y_2}} \right|}}{D} < 0.5$

    • 优化流程主要是基于柱体VIV模型对一定设计参数下的柱体振动响应进行模拟,在此基础上优化模块对设计参数进行优化使柱体振动响应接近目标条件,最后判断NES对柱体涡激振动的抑制效果是否满足输出条件。满足后输出设计参数,若不满足则将优化后的设计参数再代入第一步循环,直到输出满足条件的NES参数。优化流程如图2所示。

      图  2  优化流程图

      Figure 2.  Optimized flow chart

    • 本文采用遗传算法完成对设计参数的优化。遗传算法[21]是一种模拟自然进化过程寻找最优解的方法。20世纪70年代,遗传算法由密切根大学Holland[22-23]最先开始研究,并运用于模拟生物学的进化过程。之后,由Bagley[24]正式提出“遗传算法”一词,并发表出世界上第一篇关于遗传算法应用的论文。遗传算法通过对变量进行编码,能将实际问题转化为对编好代码进行的处理,方便解决了变量的离散性问题。遗传算法使用目标函数本身建立优化方向,不需要大量的求导求逆等运算,也能为优化结果设置约束条件,因此常常被用来求解多目标的无功优化问题。本文将遗传算法与能快速计算出结果的经验模型相结合,建立了一个NES减振装置优化设计模型,运用遗传算法对设计参数进行全局寻优,找到最接近优化目标满足条件的NES参数。

      遗传算法模型的主要求解步骤如下:

      1)对NES参数进行编码;

      2)初始化种群,确定起始搜索点的初始种群数据;

      3)根据优化目标振幅建立适应度函数,以适应度的大小判定个体的优劣,振幅越小,个体越优,遗传机率越大(求解适应度函数即求解VIV预测模型);

      4)对种群进行选择,交叉,变异操作,以迭代生成下一代遗传群体;

      5)判断遗传群体是否满足收敛条件和约束条件,若满足则结束运算,得出最优解;否则重复第4)步。

    • 为了验证Van der Pol尾流振子模型能够以一定的精度预测柱体的涡激振动响应,本文使用与Stappenbelt等[25]实验中相同的实验参数,计算得到无NES作用下的圆柱振动响应,对比实验数据验证模型。取文献中两组参数作为Van der Pol尾流振子模型的验证工况:第一组柱体直径$D = $$ 0.0554\;{\rm{ m}}$,柱体阻尼比$\varsigma = 0.0056$,圆频率$ \omega _0 =$$ 7.486\;{\rm{ rad/s}}$,质量比${m^*} = 6.54$;第二组柱体直径$D = 0.0554\;{\rm{ m}}$,柱体阻尼比$\varsigma = 0.0057$,圆频率$\omega _0 = 6.0276\;{\rm{ rad/s}}$,质量比${m^*} = 10.63$

      本文模型与Stappenbelt实验数据和CFD计算结果对比如图3所示。图3(a)为不同折减速度下柱体的振幅分布,从图中可看出,高质量比情况下,模型模拟结果在下端分支表现出与实验结果更为贴合。Van der Pol尾流振子模型与实验结果有一致的变化趋势,所以Van der Pol尾流振子模型可以一定程度上模拟出柱体的涡激振动。但由于经验公式的限制性,当$U{\rm _r}$在3~6时,模型与实验和CFD数据都有一些差距。Van der Pol尾流振子模型中有较少的经验参数,仿真的结果取决于经验参数的选取。如果经验参数并不是十分适用于该工况,就会与实验结果产生较大的偏差,而通常又无法给一个非常准又普适的经验参数。Dai等[16]同样采用Van der Pol尾流振子模型仿真柱体VIV,观察得到:该模型在频率锁定区间,基本上能捕捉到最大振幅。图3(b)为频率比随约化速度变化图,从图中可以观察到频率“锁定”现象。当$U{\rm _r}$在4~8时,柱体的振动频率与固有频率之比接近于1,同时对应于图3(a),振幅在$U{\rm _r} =$$ 4 \sim 8$时也显著大于其他折减速度下。

      图  3  基于Van der Pol尾流振子模型的VIV计算结果

      Figure 3.  The VIV simulation results based on Van der Pol wake oscillator model

      文中CFD方法计算柱体涡激振动不仅需要使用动网格,而且建立VIV响应预测模型困难。在本人的前期工作文献[26]中已建立柱体结构的涡激振动高保真仿真模型。从图3(a)中可看出,基于CFD和嵌套网格技术的涡激振动模型具有更高的仿真精度,但是计算效率远远小于Van der Pol尾流振子模型,不适用于与优化算法结合进行设计。文献[26]中同样的单核CPU计算条件下,采用建立的CFD二维弹性支撑柱体涡激振动仿真系统计算一个工况需要20 h左右,而采用Van der Pol尾流振子模型仅需要10 s左右的时间。因此,用于描述涡激升力的Van der Pol尾流振子模型计算效率高,也基本可以捕捉到柱体的涡激振动特性,可应用于本文中建立NES优化模型。

    • 运用上述建立的优化设计仿真模型,为质量$m = 15.708$ kg、刚度$k = 2530.1$ N/m、阻尼比$\varsigma = 0.0013$、圆柱直径$D = 0.02$ m、固有频率${f_n} = 2$ Hz的圆柱结构设计NES,减小其在折减速度$U{\rm _r} = 5 \sim 7$时的横向振动。最终优化获得NES参数$\beta = 0.1316$$\xi = 1.1356$$\gamma = 1.6496$,求得在此NES作用下折减速度$U{\rm _r} = 5 \sim 7$内的柱体涡激振动最大振幅为$\dfrac{{y_{1\max} }}{D} = 0.1553$,此时折减速度$U{\rm _r} = 6$

      对该NES作用下的柱体涡激振动情况进行验证,在该NES作用下,设计折减速度内的计算结果如图4所示。图4(a)为最大振幅随折减速度变化曲线图,从图中可看出,无NES情况下,折减速度$U{\rm _r} = 5.5$时,振幅接近0.5并达到最大值。NES作用情况下,在折减速度$U{\rm _r} = 3 \sim 4.5$$U{\rm _r} = 6.5 \sim 7$时,振幅都几乎为0;而在折减速度$U{\rm _r} = 4.5 \sim 6.5$区间,柱体横向振幅显著增加,并在$U{\rm _r} = 6$时,振幅达到最大值。在折减速度$U{\rm _r} = 6$时,柱体产生最大振幅,但此时${{y_{1\max} } / D}$依然在0.2以内,说明在该折减速度范围内柱体的涡激振动振幅都较小,满足设计减振要求。其中增加了文献[15]中NES参数下的振幅随折减速度变化曲线,$\beta = 0.05$$\xi = 0.8$$\gamma = 0.8$,文中发现在这几个参数下NES对涡激振动的抑制效果良好,故用于文中进行对比。从图中可看到,该NES作用下柱体振幅略小于无NES情况下,可见减振效果尚未达到目标。图4(b)为频率比随折减速度变化曲线对比图。从图中可以看出,单个柱体在折减速度$U{\rm _r} = 4.5 \sim 5.5$时,发生了频率锁定现象,对应于图4(a)中振幅也达到较大值;而加了NES以后,频率比在$U{\rm _r} = 4.5 \sim 5.5$这一区间依旧继续上升,避免了涡激共振的产生,以此达到减振效果。由此可知,通过本方法设计完成的NES对柱体涡激振动具有良好的抑制作用,可根据本模型设计方法应用于柱体减振装置设计。

      图  4  设计速度下的计算结果

      Figure 4.  Calculation results at design speed

      观察到优化NES与文献[15]中NES主要是$\beta $的变化,所以同样将文献[15]中NES的$\beta $修改为0.1和0.15进行对比,得到如图4(c)所示的最大振幅随折减速度变化曲线图。从图中可以看出,随着$\beta $的增大,NES对柱体涡激振动的抑制效果更好。优化NES的效果要优于文献[15]中NES$\beta = 0.15$的情况,同时$\beta $也小于0.15,证明了本文建立优化模型的可用性。

      柱体有无NES作用的振动响应对比如图5所示。图5(a)为振动位移对比图,振幅小的曲线为在设计NES作用下柱体产生最大振幅来流速度下$U{\rm _r} = 6$的振动位移,振幅大的曲线则是单个柱体产生最大振幅来流速度下$U{\rm _r} = 5.5$的涡激振动位移图。从图中可以看出,无NES情况下的柱体最大无量纲振幅将达到0.5,远大于NES作用情况下,说明在该折减速度范围内下柱体的涡激振动都得到了较好的抑制,满足了设计要求。图5(b)则是对应于图5(a)中振动位移曲线的功率密度谱,从图中可以看出,无NES情况下频谱曲线在2 Hz处到达峰值,对应于柱体的固有频率处,代表在该折减速度下正发生着涡激共振;红线为NES作用情况下,频谱曲线产生了更多的波动,但频谱曲线峰值避开了柱体固有频率2 Hz,避免了频率锁定的发生,这也是NES能抑制涡激共振的主要原因。

      图  5  振动响应对比图

      Figure 5.  Comparisons of vibration response

      $U{\rm _r} = 5.5$时,柱体有无NES情况下涡激振动二维相图如图6所示,图中横坐标为柱体的横向位移,纵坐标为柱体横向振动速度,外部曲线为无NES情况下的柱体振动相图,内部曲线为NES作用下柱体振动相图。可以从图中看出,无NES情况下的曲线相轨被限制在极限环上,此时柱体发生等幅振动;而NES作用下的相轨则在一定范围内波动,说明在NES作用下柱体横向振动变得不规律,但NES作用下的曲线最大半径远小于另一条曲线,说明此时NES对柱体振动起到了限制作用。

      图  6  柱体涡激振动二维相图

      Figure 6.  Two-dimensional phase diagram of VIV of the cylinder

    • 本文基于Van der Pol尾流振子模型、结构动力学和遗传算法,引入NES被动控制装置,建立了用于抑制柱体结构涡激振动的NES减振装置优化设计仿真模型。应用本方法针对一实例进行验证,对比分析了优化设计前后的柱体的振动响应结果,并研究了NES的减振机理。

      (1)用于描述涡激升力的Van der Pol尾流振子模型,相较于计算流体力学方法,计算效率高,可以用于快速捕捉到柱体的涡激振动特性。由于其高计算效率,使之能与遗传算法相结合,建立了以模型预测结果为优化函数的高效优化模型。

      (2)优化后的NES能使柱体振动错开频率锁定区域,避免产生涡激共振,从而有效抑制柱体横向振动。在$U{\rm _r} = 5.5$时候,NES作用下的柱体振幅最大减小了66.39%,减振效果明显。

      (3)不可压缩定常流动下,柱体涡激振动通常为等幅振动;而NES有效抑制柱体振动情况下,柱体往往在一定振幅范围内不规则振动,振动频率避开频率锁定区间。

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回