Parametric excitations and lock-in of flexible hydrofoils in two-phase flows

This work introduces a reduced-order method to study the parametric excitations and lock-in of flexible hydrofoils caused by unsteady two-phase (cavitating) flow. The reduced-order method is based on a 1-DOF structural model coupled with a van der Pol wake oscillator with empirically derived relations for the variation in lift, cavity-length, and cavity-shedding frequency as a function of a non-dimensional cavitation parameter. The results are compared with several available data from both numerical simulations and experimental measurements. The frequency content of both the predicted and measured vibrations suggested that, in addition to the primary cavity-shedding frequency and the hydrofoil natural frequencies, unsteady two-phase flows may excite additional modulated frequencies due to time-varying fluid-added mass effects. The results show that these frequency modulations might cause the flexible hydrofoil to undergo higher-order resonances, as well as parametric resonances. While the maximum deformations for the primary and higher-order resonances were observed to damp out, parametric resonances might persist even with realistic fluid damping coefficients (4–12%). It was observed that with higher effective foil flexibility, the cavity-shedding frequencies may be significantly modified from the rigid foil trends, and may instead lock-in with the system natural frequencies.     

带非线性弹簧的两自由度准零刚度隔振系统力传递率特性

基于多自由度系统中的反共振特性,分别在传统线性隔振系统的上、下两层引入非线性倾斜弹簧负刚度机构,构成两自由度准零刚度隔振器。通过静态特性分析,推导出系统满足零刚度条件时,各参数之间的关系,分析了力学参数及结构参数对系统刚度特性的影响。建立两自由度准零刚度隔振系统的非线性动力学方程,利用平均法求解,推导出力传递率表达式,结合数值分析方法,探讨系统在不同的上、下层隔振器阻尼比、竖直刚度比及质量比情况下的力传递率特性,并与单自由度准零刚度隔振系统及线性斜弹簧两自由度准零刚度隔振系统进行对比研究。结果表明：当结构参数 （即：倾斜弹簧处于静平衡位置的长度与倾斜弹簧原长的比值）较小且倾斜弹簧为软化弹簧时,可在平衡位置附近获得较小的系统刚度及较大的低刚度区间;通过选择适当的上、下层隔振器阻尼比、竖直刚度比与质量比,可减小系统的起始隔振频率,增宽隔振频带,加快系统力传递率在特定频段内的衰减速率,改善系统的低频隔振性能。     

Frequency lock-in is caused by coupled-mode flutter

The mechanism underlying the lock-in of frequencies in flow-induced vibrations is analysed using elementary linear dynamics. Considering the case of lock-in in vortex-induced vibrations (VIV), we use a standard wake oscillator model, as in previous studies, but in its simplest form where all nonlinear terms and all dissipative terms are neglected. The stability of the resulting linear system is analysed, and a range of coupled-mode flutter is found. In this range, the frequency of the most unstable mode is found to deviate from the Strouhal law when the frequency of the wake oscillator approaches that of the free cylinder motion. Simultaneously the growth rate resulting from coupled-mode flutter increases, which would lead to higher vibration amplitudes. The extent of the range of lock-in is then compared with experimental data, showing a good agreement. It is therefore stated that the lock-in phenomenon, such as in VIV, is a particular case of linear coupled-mode flutter.     

低雷诺数串列双方柱流致振动质量比效应的数值研究

为澄清串列双方柱流致振动的质量比效应, 采用数值模拟方法, 在雷诺数为150时, 研究了质量比($m^{\ast }=3$, 10, 20)对下游方柱振动响应特性的影响规律, 分析了下游方柱尾流模态的演变过程, 探讨了导致下游方柱振动的流固耦合机制. 结果表明: 质量比对下游方柱的流致振动有重要影响, 低质量比($m^{\ast }=3$)时下游方柱的振动响应更为复杂, 随着折减速度的增大, 下游方柱并未出现传统“锁定”现象(即振动频率比$f_{y}$/$f_{\rm n} \approx1$的锁定), 而发生了“弱锁定”现象(即$f_{y}/f_{\rm n}<1$的锁定); 随着质量比的增加($m^{\ast }=10$和20), “弱锁定”现象消失, 而出现传统“锁定”现象, 且下游方柱横流向最大振幅减小. 质量比对串列双方柱的柱心间距有明显影响, 低质量比($m^{\ast }=3$)时的柱间距在振动锁定区内会急剧减小, 而较高质量比($m^{\ast }=10$和20)下的柱间距则变化不大. 此外, 质量比对串列双方柱的尾流模态和流固耦合机制也有显著影响, 其中低质量比($m^{\ast }=3$)下的情况更为多样.      

Free v/s forced vibrations of a cylinder at low Re

Free vibrations of a circular cylinder of low non-dimensional mass are investigated at low Reynolds numbers. Computations are carried out for 5% blockage. Lock-in is observed for a range of Re and is accompanied with hysteresis at both lower as well as higher Re ends of the synchronisation/lock-in region. It is well known that the lock-in regime for free vibrations depends on the non-dimensional mass of the oscillator. The results from the present computations are compared with the data for forced vibrations from Koopmann (Journal of Fluid Mechanics, 28, 501–512, 1967) on a Y _max/D vs. f* plot, where Y _max is the maximum oscillation amplitude and f* is the ratio of cylinder vibration frequency to the vortex shedding frequency for a stationary cylinder. Good agreement is observed for the critical amplitude needed for onset of synchronisation between the forced and free vibrations. The results from the free vibrations are compared to the predictions from the linear oscillator model by assuming that the forces on the cylinder are unaffected as a result of vibrations. It is found that, for low mass oscillators, the modification of vortex shedding frequency and lift coefficient due to cylinder oscillations leads to the enhancement of the lock-in regime.     

Piezoelectric energy harvesting from concurrent vortex-induced vibrations and base excitations

We investigate the potential of using a piezoelectric energy harvester to concurrently harness energy from base excitations and vortex-induced vibrations. The harvester consists of a multilayered piezoelectric cantilever beam with a circular cylinder tip mass attached to its free end which is placed in a uniform air flow and subjected to direct harmonic excitations. We model the fluctuating lift coefficient by a van der Pol wake oscillator. The Euler–Lagrange principle and the Galerkin procedure are used to derive a nonlinear distributed-parameter model for a harvester under a combination of vibratory base excitations and vortex-induced vibrations. Linear and nonlinear analyses are performed to investigate the effects of the electrical load resistance, wind speed, and base acceleration on the coupled frequency, electromechanical damping, and performance of the harvester. It is demonstrated that, when the wind speed is in the pre- or post-synchronization regions, its associated electromechanical damping is increased and hence a reduction in the harvested power is obtained. When the wind speed is in the lock-in or synchronization region, the results show that there is a significant improvement in the level of the harvested power which can attain 150 % compared to using two separate harvesters. The results also show that an increase of the base acceleration results in a reduction in the vortex-induced vibrations effects, an increase of the difference between the resonant excitation frequency and the pull-out frequency, and a significant effects associated with the quenching phenomenon.     

Modeling and analysis of piezoaeroelastic energy harvesters

This work investigates the influence of structural and aerodynamic nonlinearities on the dynamic behavior of a piezoaeroelastic system. The system is composed of a rigid airfoil supported by nonlinear torsional and flexural springs in the pitch and plunge motions, respectively, with a piezoelectric coupling attached to the plunge degree of freedom. The analysis shows that the effect of the electrical load resistance on the flutter speed is negligible in comparison to the effects of the linear spring coefficients. The effects of aerodynamic nonlinearities and nonlinear plunge and pitch spring coefficients on the system’s stability near the bifurcation are determined from the nonlinear normal form. This is useful to characterize the effects of different parameters on the system’s output and ensure that subcritical or “catastrophic” bifurcation does not take place. Numerical solutions of the coupled equations for two different configurations are then performed to determine the effects of varying the load resistance and the nonlinear spring coefficients on the limit-cycle oscillations (LCO) in the pitch and plunge motions, the voltage output and the harvested power.     

低雷诺数下弹性圆柱体涡激振动及影响参数分析

利用Fluent软件数值求解不可压缩粘性流体的N-S方程,研究均匀来流Re=200时弹性圆柱体的涡激振动.圆柱体运动简化为质量-弹簧-阻尼系统,将Newmark-β方法代码写入用户自定义函数(UDF)求解运动方程,柱体与流体间的非线性耦合作用通过动网格技术实现.详细分析了涡激力系数、柱体位移特征值和尾流涡结构随频率比的变化关系,获得"相位开关"、"拍"等现象.考虑流向振动对横向振动影响时,圆柱体最大横向振幅为0.65倍直径.当固定频率比,而质量比或折合阻尼增大时,圆柱体流向与横向振动均呈非线性衰减趋势,但增大质量比对流向平均位移的偏离起到更好的控制效果.     

On energy pumping, synchronization and beat phenomenon in a nonideal structure coupled to an essentially nonlinear oscillator

In this paper, we investigated the effectiveness of the linear electromechanical vibration absorber (LEVA) and a nonlinear electromechanical vibration absorber (NEVA) in the vibration attenuation for nonideal structures (NIS). This electromechanical damping device consists of an electrical system coupled magnetically to a mechanical structure under a nonideal excitation. An analysis of the effects of the parameters of coupling and of nonlinear coefficients with increasing of constant torque of the DC motor is presented.     

Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

两种典型螺栓连接结构的刚度特性比较研究

针对法兰对接和径向套接两种典型舱段螺栓连接形式,基于有限元静刚度计算、动力学简化建模、冲击响应特性分析及结构冲击试验,系统研究了两种连接形式的轴向刚度特性及其对动力学冲击响应的影响.有限元静刚度分析揭示了法兰对接与径向套接的轴向拉压刚度分别是非对称的和对称的,而法兰对接的平均刚度更大.之后,为两种连接形式建立了统一的动力学模型,证明非对称的轴向拉压刚度导致结构在受到横向载荷作用时会产生附加的耦合轴向振动,并且利用高精度幂函数拟合刚度跳变,得到耦合轴向振动频率是弯曲振动频率的二倍的结论.最后,通过冲击动力学试验证明了法兰对接存在二倍频的耦合轴向振动,而径向套接则不存在该耦合振动.径向套接虽然一阶频率较低,但阻尼效果更好.     

Limit cycle oscillation suppression of 2-DOF airfoil using nonlinear energy sink

This paper presents a novel mechanical attachment, i.e., nonlinear energy sink （NES）, for suppressing the limit cycle oscillation （LCO） of an airfoil. The dynamic responses of a two-degree-of-freedom （2-DOF） airfoil coupled with an NES are studied with the harmonic balance method. Different structure parameters of the NES, i.e., mass ratio between the NES and airfoil, NES offset, NES damping, and nonlinear stiffness in the NES, are chosen for studying the effect of the LCO suppression on an aeroelastic system with a supercritical Hopf bifurcation or subcritical Hopf bifurcation, respectively. The results show that the structural parameters of the NES have different influence on the supercritical Hopf bifurcation system and the subcritical Hopf bifurcation system.     

Numerical and experimental studies on the hydrodynamic damping of a zero-thrust propeller

Fluid added mass and damping are significant parameters when predicting the dynamic response of a submerged structure. The hydrodynamic damping of underwater rotating machinery is numerically and experimentally investigated by a zero-thrust propeller in this paper. The lifting surface method(LSM) combined with forced vibration was introduced as the numerical method to compute the corresponding unsteady thrust, while the experimental method of measuring added damping was accomplished by a propeller undergoing rotation combined heave motion. Results of the theoretical method are in good agreement with the experimental results before cavitation occurs, as cavitation is regarded to weaken the unsteady response of the propeller partly. The calculation results also show that both the frequency ratio(vibration frequency divided by rotation frequency) and the blade angle have a significant influence on the hydrodynamic damping. Therefore, the effect of blade angle on hydrodynamic damping should be considered during the design phase.     

Spatially localized vibrations in a rotor subjected to flutter

The current push toward lightweight structures in aerospace and aeronautical engineering is leading to slender design airfoils, which are more likely to undergo large deformation, hence experiencing geometrical nonlinearities. The problem of vibration localization in a rotor constituted by N coupled airfoils with plunge and pitch degrees of freedom subjected to flutter instability is considered. For a single airfoil, it is shown that depending on the system parameters, multiple static and dynamic equilibria coexist which may be a fixed point, a limit cycle, or irregular motion. By elastically coupling N airfoils, a simplified rotor model is obtained. The nonlinear dynamical response of the rotor is studied via time integration with particular attention to the emergence of localized vibrating solutions, which have been classified introducing a localization coefficient. Finally, the concept of basin stability is exploited to ascertain the likelihood of the system to converge to a certain localized state as a function of the airstream velocity. We found that homogeneous and slightly localized states are more likely to appear with respect to strongly localized states.

     

Primary resonance of fractional-order van der Pol oscillator

In this paper the primary resonance of van der Pol (VDP) oscillator with fractional-order derivative is studied analytically and numerically. At first the approximately analytical solution is obtained by the averaging method, and it is found that the fractional-order derivative could affect the dynamical properties of VDP oscillator, which is characterized by the equivalent damping coefficient and the equivalent stiffness coefficient. Moreover, the amplitude–frequency equation for steady-state solution is established, and the corresponding stability condition is also presented based on Lyapunov theory. Then, the comparisons of several different amplitude–frequency curves by the approximately analytical solution and the numerical integration are fulfilled, and the results certify the correctness and satisfactory precision of the approximately analytical solution. At last, the effects of the two fractional parameters, i.e., the fractional coefficient and the fractional order, on the amplitude–frequency curves are investigated for some typical excitation amplitudes, which are different from the traditional integer-order VDP oscillator.     

Lock-in of the vortex-induced vibrations of a long tensioned beam in shear flow

The occurrence of lock-in, defined as the local synchronization between the vortex shedding frequency and the cross-flow structural vibration frequency, is investigated in the case of a tensioned beam of length to diameter ratio 200, free to move in both the in-line and cross-flow directions, and immersed in a linear shear current. Direct numerical simulation is employed at three Reynolds numbers, from 110 to 1100, so as to include the transition to turbulence in the wake. The Reynolds number influences the response amplitudes, but in all cases we observed similar fluid–structure interaction mechanisms, resulting in high-wavenumber vortex-induced vibrations consisting of a mixture of standing and traveling wave patterns.Lock-in occurs in the high oncoming velocity region, over at least 30% of the cylinder length. In the case of multi-frequency response, at any given spanwise location lock-in is principally established at one of the excited vibration frequencies, usually the locally predominant one. The spanwise patterns of the force and added mass coefficients exhibit different behaviors within the lock-in versus the non-lock-in region. The spanwise zones where the flow provides energy to excite the structural vibrations are located mainly within the lock-in region, while the flow damps the structural vibrations in the non-lock-in region.     

Dynamic design of a nonlinear energy sink with NiTiNOL-steel wire ropes based on nonlinear output frequency response functions

A novel vibration isolation device called the nonlinear energy sink(NES)with Ni Ti NOL-steel wire ropes(Ni Ti-ST) is applied to a whole-spacecraft system. The Ni Ti-ST is used to describe the damping of the NES, which is coupled with the modified Bouc-Wen model of hysteresis. The NES with Ni Ti-ST vibration reduction principle uses the irreversibility of targeted energy transfer(TET) to concentrate the energy locally on the nonlinear oscillator, and then dissipates it through damping in the NES with Ni Ti-ST.The generalized vibration transmissibility, obtained by the root mean square treatment of the harmonic response of the nonlinear output frequency response functions(NOFRFs),is first used as the evaluation index to analyze the whole-spacecraft system in the future.An optimization analysis of the impact of system responses is performed using different parameters of NES with Ni Ti-ST based on the transmissibility of NOFRFs. Finally, the effects of vibration suppression by varying the parameters of Ni Ti-ST are analyzed from the perspective of energy absorption. The results indicate that NES with Ni Ti-ST can reduce excessive vibration of the whole-spacecraft system, without changing its natural frequency. Moreover, the NES with Ni Ti-ST can be directly used in practical engineering applications.     

COUPLED ANALYSIS FOR THE HARVESTING STRUCTURE AND THE MODULATING CIRCUIT IN A PIEZOELECTRIC BIMORPH ENERGY HARVESTER

The authors analyze a piezoelectric energy harvester as an electro-mechanically coupled system. The energy harvester consists of a piezoelectric bimorph with a concentrated mass attached at one end, called the harvesting structure, an electric circuit for energy storage, and a rectifier that converts the AC output of the harvesting structure into a DC input for the storage circuit. The piezoelectric bimorph is assumed to be driven into flexural vibration by an ambient acoustic source to convert the mechanical energies into electric energies. The analysis indicates that the performance of this harvester, measured by the power density, is characterized by three important non-dimensional parameters, i.e., the non-dimensional inductance of the storage circuit, the non-dimensional aspect ratio （length/thickness） and the non-dimensional end mass of the harvesting structure. The numerical results show that： （1） the power density can be optimized by varying the non-dimensional inductance for each fixed non-dimensional aspect ratio with a fixed non-dimensional end mass; and （2） for a fixed non-dimensional inductance, the power density is maximized if the non-dimensional aspect ratio and the non-dimensional end mass are so chosen that the harvesting structure, consisting of both the piezoelectric bimorph and the end mass attached, resonates at the frequency of the ambient acoustic source.     

含机械调频式动力吸振器的准零刚度隔振系统隔振性能分析

基于动力吸振器原理,在单自由度准零刚度隔振器基础上耦合可调频动力吸振器构成两自由度隔振系统。首先,对动力吸振器工作原理进行理论分析并提出其力学模型;其次,通过静力学分析,推导出系统满足零刚度条件时,各参数间的关系并分析其对系统刚度特性的影响;然后,建立两自由度隔振系统非线性动力学方程,利用谐波平衡法进行幅频响应解析分析,得到力传递率表达式;最后,数值分析动力吸振器阻尼、刚度、质量、激励力幅值和弹簧片有效长度对力传递率的影响规律,并与单自由度准零刚度隔振系统及两自由度线性隔振系统对比分析。结果表明：通过选择适当的动力吸振器参数不仅可以减小系统的起始隔振频率,增宽隔振频带,且还能加快系统力传递率在特定频段内的衰减速率,改善系统的低频隔振性能,实现激励频率的可适应性。     

AN INTEGRATED NONLINEAR APPROACH FOR TURBOMACHINERY FORCED RESPONSE PREDICTION. PART II: CASE STUDIES

This paper discusses the application of an advanced turbomachinery forced prediction model to two representative turbomachinery cases: an HP turbine and a rig fan. The approach is based on an integrated nonlinear multi-passage analysis that includes both the stator and the rotor blade-rows. The numerical model has advanced features such as nonlinear friction damping for turbine blades, tip gap flows and blade vibratory motion. A series of inviscid and viscous computations were performed for an HP turbine with 36 stator and 92 rotor blades. The peak-to-peak maximum displacements were predicted with and without root friction dampers and the findings were compared with available experimental data. Good agreement was observed in most cases. It was found that most of the unsteady forcing was due to the potential effects. In a second phase of the turbine work, the response to low engine-order excitation was predicted using a multi-row whole-annulus model. The stator assembly was assumed to have blades with varying throat widths and the magnitude of the unsteady aerodynamic forcing was found to increase with increasing scatter in throat width variation. A second forced response study was conducted for a rig fan with 15 variable inlet guide vanes (VIGVs) and 20 rotor blades. For a 30° VIGV opening, a good match was observed between the predicted and measured wake profiles. Similarly, the measured and predicted rotor blade vibration levels were also found to be in good agreement. It is concluded that the proposed methodology can be applied, with reasonable confidence, to the study of industrial cases.