Effects of pretwist and presetting on coupled bending vibrations of rotating thin-walled composite beams

A refined dynamic theory of rotating blades modeled as anisotropic composite thin-walled beams, experiencing the flapping-lagging-transverse shear coupling is presented. The structural model encompasses a number of non-standard features, such as anisotropy and transverse shear, pretwist and presetting angles, the presence of a rigid hub on which the beam is mounted, and the rotatory inertia. The developed theory and the methodology used to determine the eigenfrequency characteristics are validated against the results available in the literature, and new results emphasizing the influence played by the ply-angle, pretwist and presetting, coupled with that of the rotating speed on blade free vibration characteristics are supplied, and pertinent conclusions are outlined.     

NONLINEAR DYNAMIC RESPONSE OF PRE-DEFORMED BLADE WITH VARIABLE ROTATIONAL SPEED UNDER 2:1 INTERNAL RESONANCE$^{\bf 1)}$

在工程实际中,涡轮机叶片的转速在很多应用场景下不是一个定常值,比如发动机在启动、变速、停机等工况下,转子输入与输出功率失衡,伴随产生扭振,产生速度脉冲. 另外,由于服役环境、安装误差等因素会引起叶片在所难免的预变形. 本文主要研究预变形叶片,在变转速条件下的非线性动力学行为. 考虑叶片转速由一定常转速和一简谐变化的微小扰动叠加而成. 应用拉格朗日原理得到变转速叶片的动力学控制方程,并采用假设模态法将偏微分方程转为常微分方程,通过引入无量纲,使方程更具有一般性. 运用多尺度方法求解了该参激振动系统,得到了在 2:1 内共振情形下的平均方程,进而获得系统的稳态响应. 详细研究温度梯度、阻尼以及转速扰动幅值等系统参数对叶片动力学响应的影响规律,同时考察了立方项在 2:1 内共振下对方程的影响. 对原动力方程进行正向、反向扫频积分来观察其跳跃现象,并对解析解进行验证. 结果发现参数的变化对叶片均有不同程度影响,在 2:1 内共振下立方项对系统响应的影响很小,解析解与数值解吻合很好.     

2:1内共振条件下变转速预变形叶片的非线性动力学响应

在工程实际中,涡轮机叶片的转速在很多应用场景下不是一个定常值,比如发动机在启动、变速、停机等工况下,转子输入与输出功率失衡,伴随产生扭振,产生速度脉冲. 另外,由于服役环境、安装误差等因素会引起叶片在所难免的预变形. 本文主要研究预变形叶片,在变转速条件下的非线性动力学行为. 考虑叶片转速由一定常转速和一简谐变化的微小扰动叠加而成. 应用拉格朗日原理得到变转速叶片的动力学控制方程,并采用假设模态法将偏微分方程转为常微分方程,通过引入无量纲,使方程更具有一般性. 运用多尺度方法求解了该参激振动系统,得到了在 2:1 内共振情形下的平均方程,进而获得系统的稳态响应. 详细研究温度梯度、阻尼以及转速扰动幅值等系统参数对叶片动力学响应的影响规律,同时考察了立方项在 2:1 内共振下对方程的影响. 对原动力方程进行正向、反向扫频积分来观察其跳跃现象,并对解析解进行验证. 结果发现参数的变化对叶片均有不同程度影响,在 2:1 内共振下立方项对系统响应的影响很小,解析解与数值解吻合很好.      

Non-linear dynamic response of a rotating radial Timoshenko beam with periodic pulse loading at the free-end

Consideration is given to the dynamic response of a Timoshenko beam under repeated pulse loading. Starting with the basic dynamical equations for a rotating radial cantilever Timoshenko beam clamped at the hub in a centrifugal force field, a system of equations are derived for coupled axial and lateral motions which includes the transverse shear and rotary inertia effects, as well. The hyperbolic wave equation governing the axial motion is coupled with the flexural wave equation governing the lateral motion of the beam through the velocity-dependent skew-symmetric Coriolis force terms. In the analytical formulation, Rayleigh-Ritz method with a set of sinusoidal displacement shape functions is used to determine stiffness, mass and gyroscopic matrices of the system. The tip of the rotating beam is subjected to a periodic pulse load due to local rubbing against the outer case introducing Coulomb friction in the system. Transient response of the beam with the tip deforming due to rub is discussed in terms of the frequency shift and non-linear dynamic response of the rotating beam. Numerical results are presented for this vibro-impact problem of hard rub with varying coefficients of friction and the contact-load time. The effects of beam tip rub forces transmitted through the system are considered to analyze the conditions for dynamic stability of a rotating blade with intermittent rub.     

风力机叶片非线性挥舞分析

将风力机叶片简化为绕轮毂旋转的变截面Euler-Bernoulli悬臂梁,基于Greenberg公式给出非线性气动力,建立叶片挥舞振动非线性控制方程.由于变截面梁的弯曲刚度和线密度是沿梁轴线变化的函数,无法给出模态函数解析式,论文提出使用假设模态法计算的模态函数,作为基函数对控制方程进行Galerkin截断,通过将挥舞振动分解为静态位移和动态扰动合成,对其进行动态响应分析,同时讨论了叶轮转速、风速和旋转位置对振动特性的影响.研究表明:(1)叶轮转速对叶片挥舞特性影响显著,风速和叶片转角对振动特性影响很小.(2)静态位移随风速增加而增大,大体上成线性关系,气动阻尼随风速增加而减小.(3)风速较低时,非线性挥舞振动表现为衰减振动,随着风速增加,振动由衰减振动演化为周期运动,再由周期运动演化为拟周期运动.     

NONLINEAR DYNAMICS OF AXIALLY ACCELERATING VISCOELASTIC BEAMS BASED ON DIFFERENTIAL QUADRATURE

This paper investigates nonlinear dynamical behaviors in transverse motion of an axially accelerating viscoelastic beam via the differential quadrature method. The governing equation, a nonlinear partial-differential equation, is derived from the viscoelastic constitution relation using the material derivative. The differential quadrature scheme is developed to solve numerically the governing equation. Based on the numerical solutions, the nonlinear dynamical behaviors are identified by use of the Poincare map and the phase portrait. The bifurcation diagrams are presented in the case that the mean axial speed and the amplitude of the speed fluctuation are respectively varied while other parameters are fixed. The Lyapunov exponent and the initial value sensitivity of the different points of the beam, calculated from the time series based on the numerical solutions, are used to indicate periodic motions or chaotic motions occurring in the transverse motion of the axially accelerating viscoelastic beam.     

Aerothermoelastic behavior of supersonic rotating thin-walled beams made of functionally graded materials

In this study, a thin-walled beam made of functionally graded material (FGM) which is used as rotating blades in turbomachinery under aerothermoelastic loading is investigated. The governing equations, which are based on first-order shear deformation theory, include the effects of the presetting angle, the secondary warping, temperature gradient through the wall thickness of the beam and also the rotational speed. Moreover, quasi-steady aerodynamic pressure loadings are determined using first-order piston theory, and steady beam surface temperature is obtained from gas dynamics theory. Then, the blade partial differential equations are transformed into a set of ordinary differential equations using the extended Galerkin method. Finally, having solved the resulting structural–fluid–thermal eigenvalue system of equations, the effects of Mach number and geometric parameters on natural frequencies are presented. The results demonstrate that the natural frequencies decrease under aerothermoelastic loading at high Mach numbers.     

Studies on Bifurcation and Chaos of a String-Beam Coupled System with Two Degrees-of-Freedom

In this paper, research on nonlinear dynamic behavior of a string-beam coupled system subjected to parametric and external excitations is presented. The governing equations of motion are obtained for the nonlinear transverse vibrations of the string-beam coupled system. The Galerkin's method is employed to simplify the governing equations to a set of ordinary differential equations with two degrees-of-freedom. The case of 1:2 internal resonance between the modes of the beam and string, principal parametric resonance for the beam, and primary resonance for the string is considered. The method of multiple scales is utilized to analyze the nonlinear responses of the string-beam coupled system. Based on the averaged equation obtained here, the techniques of phase portrait, waveform, and Poincare map are applied to analyze the periodic and chaotic motions. It is found from numerical simulations that there are obvious jumping phenomena in the resonant response–frequency curves. It is indicated from the phase portrait and Poincare map that period-4, period-2, and periodic solutions and chaotic motions occur in the transverse nonlinear vibrations of the string-beam coupled system under certain conditions. An erratum to this article is available at .     

Nonlinear vibrations of a composite laminated cantilever rectangular plate with one-to-one internal resonance

The nonlinear vibrations of a composite laminated cantilever rectangular plate subjected to the in-plane and transversal excitations are investigated in this paper. Based on the Reddy??s third-order plate theory and the von Karman type equations for the geometric nonlinearity, the nonlinear partial differential governing equations of motion for the composite laminated cantilever rectangular plate are established by using the Hamilton??s principle. The Galerkin approach is used to transform the nonlinear partial differential governing equations of motion into a two degree-of-freedom nonlinear system under combined parametric and forcing excitations. The case of foundational parametric resonance and 1:1 internal resonance is taken into account. The method of multiple scales is utilized to obtain the four-dimensional averaged equation. The numerical method is used to find the periodic and chaotic motions of the composite laminated cantilever rectangular plate. It is found that the chaotic responses are sensitive to the changing of the forcing excitations and the damping coefficient. The influence of the forcing excitation and the damping coefficient on the bifurcations and chaotic behaviors of the composite laminated cantilever rectangular plate is investigated numerically. The frequency-response curves of the first-order and the second-order modes show that there exists the soft-spring type characteristic for the first-order and the second-order modes.     

Nonlinear energy harvesting with dual resonant zones based on rotating system

An electromagnetic nonlinear energy harvester(NEH)based on a rotating system is proposed,of which the host system rotates at a constant speed and vibrates harmonically in the vertical direction.This kind of device exhibits several resonant phenomena due to the combinations of the rotating and the vibration frequencies of the host system as well as the cubic nonlinearity of the NEH.The governing equation of motion for the NEH is derived,and the dynamic responses and output power are investigated with the multiple scale method under the 1:1 primary and 2:1 superharmonic resonant conditions.The effects of system parameters including the nondimensional external frequency,the rotating speed,and the nonlinear stiffness on the responses of free vibration for the system are studied.The results of the primary resonance show that the responses exhibit not only the resonant characteristics but also the nonlinear dynamic characteristics such as the saddle-node(SN)bifurcation.The coexistence of multiple solutions and the varying trends of responses are verified with the direct numerical simulation.Moreover,the effects of system parameters on the average output power are investigated.The results of the analyses on the two resonant conditions indicate that the large power can be harvested in two resonant frequency bands.The effect of resonance on the output power is dominant for the 2:1 superharmonic resonance.Moreover,the results also show that introducing the nonlinearity can increase the value of the output power in large frequency bands and induce the occurence of new frequency bands to harvest the large power.The efficiency of the harvested power could be improved by the combined effects of the resonance as well as the nonlinearity of the NEH device.Suitable parameter conditions could help optimize the power harvesting in design.     

Dynamic analysis of a simply supported beam resting on a nonlinear elastic foundation under compressive axial load using nonlinear normal modes techniques under three-to-one internal resonance condition

In this paper, the Nonlinear Normal Modes (NNMs) analysis for the case of three-to-one (3:1) internal resonance of a slender simply supported beam in presence of compressive axial load resting on a nonlinear elastic foundation is studied. Using the Euler?CBernoulli beam model, the governing nonlinear PDE of the beam??s transverse vibration and also its associated boundary conditions are extracted. These nonlinear motion equation and boundary condition relations are solved simultaneously using four different approximate-analytical solution techniques, namely the method of Multiple Time Scales, the method of Normal Forms, the method of Shaw and Pierre, and the method of King and Vakakis. The obtained results at this stage using four different methods which are all in time?Cspace domain are compared and it is concluded that all the methods result in a similar answer for the amplitude part of the transverse vibration. At the next step, the nonlinear normal modes are obtained. Furthermore, the effect of axial compressive force in the dynamic analysis of such a beam is studied. Finally, under three-to-one-internal resonance condition the NNMs of the beam and the steady-state stability analysis are performed. Then the effect of changing the values of different parameters on the beam??s dynamic response is also considered. Moreover, 3-D plots of stability analysis in the steady-state condition and the beam??s amplitude frequency response curves are presented.     

Global Bifurcations and Chaotic Dynamics in Nonlinear Nonplanar Oscillations of a Parametrically Excited Cantilever Beam

This paper presents the analysis of the global bifurcations and chaotic dynamics for the nonlinear nonplanar oscillations of a cantilever beam subjected to a harmonic axial excitation and transverse excitations at the free end. The governing nonlinear equations of nonplanar motion with parametric and external excitations are obtained. The Galerkin procedure is applied to the partial differential governing equation to obtain a two-degree-of-freedom nonlinear system with parametric and forcing excitations. The resonant case considered here is 2:1 internal resonance, principal parametric resonance-1/2 subharmonic resonance for the in-plane mode and fundamental parametric resonance–primary resonance for the out-of-plane mode. The parametrically and externally excited system is transformed to the averaged equations by using the method of multiple scales. From the averaged equation obtained here, the theory of normal form is applied to find the explicit formulas of normal forms associated with a double zero and a pair of pure imaginary eigenvalues. Based on the normal form obtained above, a global perturbation method is utilized to analyze the global bifurcations and chaotic dynamics in the nonlinear nonplanar oscillations of the cantilever beam. The global bifurcation analysis indicates that there exist the heteroclinic bifurcations and the Silnikov type single-pulse homoclinic orbit in the averaged equation for the nonlinear nonplanar oscillations of the cantilever beam. These results show that the chaotic motions can occur in the nonlinear nonplanar oscillations of the cantilever beam. Numerical simulations verify the analytical predictions.     

Nonlinear analysis of the dynamics of articulated composite rotor blades

The general nonlinear intrinsic equations of motion of an elastic composite beam are solved in order to obtain the elasto-dynamic response of a rotating articulated blade. The solution utilizes the linear Variational-Asymptotic Method (VAM) cross-sectional analysis, together with an improved damped nonlinear model for the rigid-body motion analysis of helicopter blades in coupled flap and lead-lag motions. The explicit (direct) integration algorithm implements the perturbation method in order to solve the transient form of the nonlinear intrinsic differential equations of motion and obtain the elasto-dynamic behavior of an accelerating composite blade. The specific problem considered is an accelerating articulated helicopter blade of which its motion is analyzed since it starts rotating from rest until it reaches the steady-state condition. It is observed that the steady-state solution obtained by this method compares very well with other available solutions. The resulting simulation code is a powerful tool for analyzing the nonlinear response of composite rotor blades; and for serving the ultimate aim of efficient noise and vibration control in helicopters.     

Nonlinear aeroelastic stability analysis of wind turbine blade with bending–bending–twist coupling

In this study, the nonlinear aeroelastic stability of wind turbine blade with bending–bending–twist coupling has been investigated for composite thin-walled structure with pretwist angle. The aerodynamic model used here is the differential dynamic stall nonlinear ONERA model. The nonlinear aeroelastic equations are reduced to ordinary equations by Galerkin method, with the aerodynamic force decomposition by strip theory. The nonlinear resulting equations are solved by a time-marching approach, and are linearized by small perturbation about the equilibrium point. The nonlinear aeroelastic stability characteristics are investigated through eigenvalue analysis, nonlinear time domain response, and linearized time domain response.     

气动载荷下旋转运动圆板磁弹性主共振分析

基于Kirchhoff薄板理论与哈密顿原理,建立旋转运动导电圆板的磁-气动弹性非线性动力学方程．根据电磁场基本原理得到旋转运动圆板所受电磁力表达式,同时采用一种简化的气动模型以描述作用于板上的气动载荷．基于贝塞尔函数形式振型函数的选取,应用伽辽金法得到旋转圆板的磁气动弹性轴对称非线性振动微分方程．应用多尺度法推导出主共振下系统的幅频响应方程,并依据Lyapunov方法得到系统稳态运动稳定性判据．通过算例,得到周边夹之约束下圆板主共振的幅频特性曲线图,以及振幅随磁感应强度和激励力幅值的变化曲线图;阐述了不同参数对系统共振幅值的影响规律,并对解的稳定性进行了分析．     

时变张力作用下轴向运动梁的分岔与混沌

轴向运动结构的横向参激振动一直是非线性动力学领域的研究热点之一. 目前研究较多的是轴向速度摄动的动力学模型,参数激励由速度的简谐波动产生. 但在工程应用中,存在轴向张力波动的运动结构较为广泛,而针对轴向张力摄动的模型研究较少. 本文研究了时变张力作用下轴向变速运动黏弹性梁的分岔与混沌. 考虑随着时间周期性变化的轴向张力,计入线性黏性阻尼,采用Kelvin模型的黏弹性本构关系,给出了梁横向非线性 振动的积分--偏微分控制方程. 首先应用四阶Galerkin截断方法将控制方程离散化,然后采用四阶Runge-Kutta方法计算系统的数值解,进而确定其动力学行为. 基于梁中点的横向位移和速度的数值结果,仿真了梁沿平均轴速、张力摄动幅值、张力摄动频率以及黏弹性系数变化的倍周期分岔与混 沌运动,并且通过计算系统的最大李雅普诺夫指数来识别其混沌行为. 结果表明:较小的平均轴速有助于梁的周期运动,梁在临界速度附近容易发生倍周期分岔与混沌行为. 随着张力摄动幅值的增大,梁的振动幅值的混沌区间不断增大. 较小的黏弹性系数和张力摄动频率更容易使梁发生混沌运动. 最后,给出时程图、频谱图、相图以及Poincaré 映射图来确定梁的混沌运动.      

Thermo-mechanical nonlinear dynamics of a buckled axially moving beam

The thermo-mechanical nonlinear dynamics of a buckled axially moving beam is numerically investigated, with special consideration to the case with a three-to-one internal resonance between the first two modes. The equation of motion of the system traveling at a constant axial speed is obtained using Hamilton??s principle. A closed form solution is developed for the post-buckling configuration for the system with an axial speed beyond the first instability. The equation of motion over the buckled state is obtained for the forced system. The equation is reduced into a set of nonlinear ordinary differential equations via the Galerkin method. This set is solved using the pseudo-arclength continuation technique to examine the frequency response curves and direct-time integration to construct bifurcation diagrams of Poincaré maps. The vibration characteristics of the system at points of interest in the parameter space are presented in the form of time histories, phase-plane portraits, and Poincaré sections.     

Complicated nonlinear responses of a simply supported FGM rectangular plate under combined parametric and external excitations

In this paper, we use the asymptotic perturbation method based on the Fourier expansion and the temporal rescaling to investigate the nonlinear oscillations and chaotic dynamics of a simply supported rectangular plate made of functionally graded materials (FGMs) subjected to a through-thickness temperature field together with parametric and external excitations. Material properties are assumed to be temperature-dependent. Based on the Reddy’s third-order plate theory, the governing equations of motion for the plate are derived using the Hamilton’s principle. The Galerkin procedure is employed to obtain a two-degree-of-freedom nonlinear system including the quadratic and cubic nonlinear terms. The resonant case considered here is 1:2 internal resonance, principal parametric resonance-1/2 subharmonic resonance. Based on the averaged equation in polar coordinate form, the stability of steady state solutions is analyzed. The phase portrait, waveform and Poincaré map are used to analyze the periodic and chaotic motions of the FGM rectangular plate. It is found that the FGM rectangular plate exhibits the chaotic motions under certain circumstances. It is seen that the nonlinear dynamic responses of the FGM rectangular plate are more sensitive to transverse excitation. The excitation force can be used as a controlling factor which can change the response of the FGM rectangular plate from periodic motion to the chaotic motion.     

确定性载荷作用下Timoshenko薄壁梁的弯扭耦合动力响应

建立了一种普遍的解析理论用于求解确定性载荷作用下Timoshenko薄壁梁的弯扭耦合动力响应。首先通过直接求解单对称均匀Timoshenko薄壁梁单元弯扭耦合振动的运动偏微分方程，给出了计算其自由振动的精确方法，并导出了Timoshenko弯扭耦合薄壁梁自由振动主模态的正交条件。然后利用简正模态法研究了确定性载荷作用下单对称Timoshenko薄壁梁的弯扭耦合动力响应，该弯扭耦合梁所受到的荷载可以是集中载荷或沿着梁长度分布的分布载荷。最后假定确定性载荷是谐波变化的，得到了各种激励下封闭形式的解，并对动力弯曲位移和扭转位移的数值结果进行了讨论。     

耦合变形对大范围运动柔性梁动力学建模的影响

柔性梁在作大范围空间运动时,产生弯曲和扭转变形,这些变形的相互耦合形成了梁在纵向以及横向位移的二次耦合变量。本文考虑了变形产生的几何非线性效应对运动柔性梁的影响,在其三个方向的变形中均考虑了二次耦合变量,利用弹性旋转矩阵建立了准确的几何非线性变形方程,通过Lagrange方程导出系统的动力学方程。仿真结果表明,在大范围运动情况下,仅在纵向变形中计及了变形二次耦合量的一次动力学模型,与考虑了完全几何非线性变形的模型具有一定的差异。