谐波平衡法与复平均法计算强非线性系统稳态响应的区别

近些年，很多学者致力于利用非线性增强振动响应减少的效果或者能量采集器的效率。因而非线性系统的响应值需要从理论计算方面更准确地预测。另外，根据学者已取得的研究成就，非线性能量汇(NES)中存在的立方刚度非线性可以将结构中宽频域的振动能量传递至非线性振子部分。文章将一种由NES和压电能量采集器组成的NES-piezo装置与两自由度主结构耦合连接，系统受谐和激励作用。文章采用谐波平衡法和复平均法分别推导了系统稳态响应，参照数值结果，对比两种近似解析方法在求解强非线性系统稳态响应时的异同。计算结果表明，系统体现较弱非线性时，二者计算结果差异很小；当系统体现强非线性时，复平均法不能准确地呈现系统高阶响应，提高阶数的谐波平衡法能更准确地表示系统响应值。基于谐波平衡法和数值算法，讨论NES-piezo装置对于系统宽频域减振的影响。与仅加入非线性能量汇情况对比，结果表明NES-piezo装置不会恶化宽频域减振效果，并且在第一阶共振频率附近，可以稍微提高结构减振效率。另外，计算结果也表明，采用恰当的NES-piezo装置可实现宽频域范围的结构减振和压电能量采集一体化。此项研究工作为研究不同情形强非线性系统的响应提供了理论方法的指导。另外，研究结果也为宽频域范围的结构减振和压电能量采集一体化提供了理论依据。     

Responses of Duffing oscillator with fractional damping and random phase

This paper aims to investigate dynamic responses of stochastic Duffing oscillator with fractional-order damping term, where random excitation is modeled as a harmonic function with random phase. Combining with Lindstedt–Poincaré (L–P) method and the multiple-scale approach, we propose a new technique to theoretically derive the second-order approximate solution of the stochastic fractional Duffing oscillator. Later, the frequency–amplitude response equation in deterministic case and the first- and second-order steady-state moments for the steady state in stochastic case are presented analytically. We also carry out numerical simulations to verify the effectiveness of the proposed method with good agreement. Stochastic jump and bifurcation can be found in the figures of random responses, and then we apply Monte Carlo simulations directly to obtain the probability density functions and time response diagrams to find the stochastic jump and bifurcation. The results intuitively show that the intensity of the noise can lead to stochastic jump and bifurcation.     

Nonlinear vibration of a rotating laminated composite circular cylindrical shell: traveling wave vibration

In this paper, the large-amplitude (geometrically nonlinear) vibrations of rotating, laminated composite circular cylindrical shells subjected to radial harmonic excitation in the neighborhood of the lowest resonances are investigated. Nonlinearities due to large-amplitude shell motion are considered using the Donnell’s nonlinear shallow-shell theory, with account taken of the effect of viscous structure damping. The dynamic Young’s modulus which varies with vibrational frequency of the laminated composite shell is considered. An improved nonlinear model, which needs not to introduce the Airy stress function, is employed to study the nonlinear forced vibrations of the present shells. The system is discretized by Galerkin’s method while a model involving two degrees of freedom, allowing for the traveling wave response of the shell, is adopted. The method of harmonic balance is applied to study the forced vibration responses of the two-degrees-of-freedom system. The stability of analytical steady-state solutions is analyzed. Results obtained with analytical method are compared with numerical simulation. The agreement between them bespeaks the validity of the method developed in this paper. The effects of rotating speed and some other parameters on the nonlinear dynamic response of the system are also investigated.     

Constrained parameter-splitting perturbation method for the improved solutions to the nonlinear vibrations of Euler–Bernoulli cantilevers

In this paper, a new method named constrained parameter-splitting perturbation method for improving the solutions obtained from the parameter-splitting perturbation method is proposed for solving the problems in some extremal cases, such as the strongly nonlinear vibration of an Euler–Bernoulli cantilever. The proposed method takes the advantages of both the perturbation method and the harmonic balance method. The idea is that the solution obtained by the parameter-splitting perturbation method is substituted into the equation of motion and then the accumulative error of the equation is minimized for determining the unknown splitting parameters under the constraints constructed under the frame of harmonic balance method. The forced vibration of an oscillator with cubic geometric nonlinearity and inertia nonlinearity and the forced vibration of a planar microcantilever beam with a lumped tip mass are studied as examples to reveal the efficacy of the proposed method. The inspection of the steady-state response including its stability is conducted by means of comparing the frequency-response curves obtained by the proposed method with those obtained by the numerical continuation method and harmonic balance method, respectively, to show the efficacy and the advantages of the proposed method. Meanwhile, the nonlinear ordering effect on the solutions of the proposed method is also studied by comparing the results obtained by using different nonlinear orderings in the systems. In the last, we found through convergence examinations that it is necessary to have corrections to the erroneous solution which are obtained by harmonic balance method and Floquet theory in stability analysis.

     

Analytical solutions for resonant response of suspended cables subjected to external excitation

In this study, two analytical methods are applied to study the primary resonances response of suspended cables subjected to external excitation. We choose four different sag-to-span ratios and the first two modes to investigate the differences in nonlinear responses obtained with analytical methods. First, we summarize the equations of motion by applying the Hamilton’s principle and quasi-static assumption, and then these equations are discretized by the Galerkin procedure. Second, the multiple-scale method and homotopy analysis method are adopted to obtain the approximate solutions. Moreover, numerical integrations are introduced in order to verify the obtained approximate results. The numerical results show that frequency response curves obtained by different analytical methods show different quantitative predictions in some cases of motion, modes, and particular sag-to-span ratios. Finally, the differences in displacement fields and axial tension forces are compared and analyzed.     

Performance analysis of a nonlinear inerter-based vibration isolator with inerter embedded in a linkage mechanism

This study presents an inerter-based nonlinear vibration isolator with geometrical nonlinearity created by configuring an inerter in a diamond-shaped linkage mechanism. The isolation performance of the proposed nonlinear isolator subjected to force or base-motion excitations is investigated. Both analytical and alternating frequency-time harmonic balance methods as well as numerical integration method are used to obtain the dynamic response. Beneficial performance of the nonlinear isolator is demonstrated by various performance indices including the force and displacement transmissibility as well as power flow variables. It is found that the use of the nonlinear inerter in the isolator can shift and bend the peaks of the transmissibility and time-averaged power flow to the low-frequency range, creating a larger frequency band of effective vibration isolation. It is also shown that the inertance-to-mass ratio and the initial distance of the nonlinear inerter can be effectively tailored to achieve reduced transmissibility and power transmission at interested frequencies. Anti-resonant peaks appear at specific frequency, creating near-zero energy transmission and significantly reducing vibration transmission to a base structure on which the proposed isolator is mounted.

     

静载荷对夹层圆板超谐波共振的影响

研究静载荷作用下夹层圆板的超谐波共振问题．基于Hoff型夹层板理论,给出了静载荷作用下夹层圆板的非线性动力学方程．应用Galerkin法推导了静载荷作用下夹层圆板的轴对称非线性振动方程．运用多尺度法分别对系统的三次超谐波问题和二次超谐波问题进行了求解,并依据Lyapunov稳定性理论得到了系统稳态运动的稳定性判据．通过算例,得到了周边简支约束下夹层圆板三次超谐波共振和二次超谐波共振的幅频响应曲线图、振幅-静载荷响应曲线图、振幅-激励力幅值响应曲线图;研究了不同参数对系统振幅的影响规律,并对解的稳定性进行了分析．     

Subharmonic resonance of a clamped-clamped buckled beam with 1:1 internal resonance under base harmonic excitations

The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated. The nonlinear partial integrodifferential equation of the motion of the buckled beam with both quadratic and cubic nonlinearities is given by using Hamilton's principle. A set of second-order nonlinear ordinary differential equations are obtained by spatial discretization with the Galerkin method. A high-dimensional model of the buckled beam is derived, concerning nonlinear coupling. The incremental harmonic balance (IHB) method is used to achieve the periodic solutions of the high-dimensional model of the buckled beam to observe the nonlinear frequency response curve and the nonlinear amplitude response curve, and the Floquet theory is used to analyze the stability of the periodic solutions. Attention is focused on the subharmonic resonance caused by the internal resonance as the excitation frequency near twice of the first natural frequency of the buckled beam with/without the antisymmetric modes being excited. Bifurcations including the saddle-node, Hopf, perioddoubling, and symmetry-breaking bifurcations are observed. Furthermore, quasi-periodic motion is observed by using the fourth-order Runge-Kutta method, which results from the Hopf bifurcation of the response of the buckled beam with the anti-symmetric modes being excited.     

Primary resonance of coupled cantilevers subjected to magnetic interaction

Based on a distributed-parameter model, the forced vibration of a cantilever pair excited by a sinusoidal base movement is analyzed. Two cantilevers are coupled at their free ends by a linear spring. A nonlinear concentrated magnetic force acts on the tip of one cantilever, serving at the nonlinear boundary condition of the continuous model. The magnetic force is modeled as a fractional function, strongly dependent on the distance between two magnets. Via the method of multiple scales, the primary resonance is analyzed for all modes. A second-order approximate solution and its stability condition are analytically captured. It is revealed that the frequency–response curves are sensitive to the distance between the two magnets. The curve may exhibit the hardening-type, softening-type or linear behavior due to the existence of the quadratic nonlinearity. The outcomes are supported by the numerical simulations very well.     

发电机组转子轴系扭振模化系统的三重共振

建立了发电机组转子轴系扭振模化系统的数学模型，此模型具有平方非线性并受简谐激励作用．研究了系统的固有频率存在双重内共振关系ω３≈２ω２，ω２≈２ω１且Ω≈ω２时的三重共振问题．文中应用非线性振动的改进平均法，求得了系统三重共振的一次近似解，对三重共振的定常解进行了理论分析和数值计算，并进行了奇异性分析，文中指出三重共振解具有双饱和特性，对二种主要的理论分析和数值计算结果进行了实验验证，实验结果与理论结果相符．     

Tuning methods for tuned inerter dampers coupled to nonlinear primary systems

This study develops displacement- and kinetic energy-based tuning methods for the design of the tuned inerter dampers (TIDs) coupled to both linear and nonlinear primary systems. For the linear primary system, the design of the TID is obtained analytically. The steady-state frequency–response relationship of the nonlinear primary system with a softening or hardening stiffness nonlinearity is obtained using the harmonic balance (HB) method. Analytical and numerical tuning approaches based on HB results are proposed for optimal designs of the TID to achieve equal peaks in the response curves of the displacement and the kinetic energy of the primary system. Via the developed approaches, the optimal stiffness of the TID can be obtained according to the stiffness nonlinearity of the primary system and the inertance of the absorber. Unlike the linear primary oscillator case, for a nonlinear primary oscillator the shape of its resonant peaks is mainly affected by the damping ratio of the TID, while the peak values depend more on the stiffness ratio. The proposed designs are shown to be effective in a wide range of stiffness nonlinearities and inertances. This study demonstrates the benefits of using inerters in vibration suppression devices, and the adopted methods are directly applicable for nonlinear systems with different types of nonlinearities.

     

A symbolic algorithm for the automatic computation of multitone-input harmonic balance equations for nonlinear systems

A symbolic algorithm is developed for the automatic generation of harmonic balance equations for multitone input for a class of nonlinear differential systems with polynomial nonlinearities. Generalized expressions are derived for the construction of balance equations for a defined multitone signal form. Procedures are described for determining combinations for a given output frequency from the desired set obtained from box truncated spectra and their permutations to automate symbolic algorithm. An application of method is demonstrated using the well-known Duffing–Van der Pol equation. Then the obtained analytical results are compared with numerical simulations to show the accuracy of the approach. The computation times for both the numerical solutions of equations versus the number of frequency components and the symbolic generation of the equations versus the power of nonlinearity are also investigated.     

Complex dynamics of a harmonically excited structure coupled with a nonlinear energy sink

Nonlinear behaviors are investigated for a structure coupled with a nonlinear energy sink. The structure is linear and subject to a harmonic excitation, modeled as a forced single-degree-of-freedom oscillator. The nonlinear energy sink is modeled as an oscillator consisting of a mass,a nonlinear spring, and a linear damper. Based on the numerical solutions, global bifurcation diagrams are presented to reveal the coexistence of periodic and chaotic motions for varying nonlinear energy sink mass and stiffness. Chaos is numerically identified via phase trajectories, power spectra,and Poincaré maps. Amplitude-frequency response curves are predicted by the method of harmonic balance for periodic steady-state responses. Their stabilities are analyzed.The Hopf bifurcation and the saddle-node bifurcation are determined. The investigation demonstrates that a nonlinear energy sink may create dynamic complexity.     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

Stability analysis of duffing oscillator with time delayed and/or fractional derivatives

The periodic motions of the fractional order and/or delayed nonlinear systems are investigated in the frequency domain using a harmonic balance method with the analytical gradients of the nonlinear quality constraints and the sensitivity information of the Fourier coefficients can also obtained. The properties of fractional order derivatives and trigonometric functions are utilized to construct the fractional order derivatives, delayed and product operational matrices. The operational matrices are used to derive the analytical formulae of nonlinear systems of algebraic equations. The stability of periodic solutions for the delayed nonlinear systems is identified by an eigenvalue analysis of quasi-polynomials characteristic equations. Sensitivity analysis is performed to study the influence of the structural parameters on the system responses. Finally, three numerical examples are presented to illustrate the validity and feasibility of the developed method. It is concluded that the proposed methodology has the potential to facilitate highly efficient optimization, as well as sensitivity and uncertainty analysis of nonlinear systems with fractional derivatives and/or time delayed.     

Nonlinear dynamic analysis of a quarter vehicle system with external periodic excitation

In this paper, a nonlinear dynamic model of a quarter vehicle with nonlinear spring and damping is established. The dynamic characteristics of the vehicle system with external periodic excitation are theoretically investigated by the incremental harmonic balance method and Newmark method, and the accuracy of the incremental harmonic balance method is verified by comparing with the result of Newmark method. The influences of the damping coefficient, excitation amplitude and excitation frequency on the dynamic responses are analyzed. The results show that the vibration behaviors of the vehicle system can be control by adjusting appropriately system parameters with the damping coefficient, excitation amplitude and excitation frequency. The multi-valued properties, spur-harmonic response and hardening type nonlinear behavior are revealed in the presented amplitude-frequency curves. With the changing parameters, the transformation of chaotic motion, quasi-periodic motion and periodic motion is also observed. The conclusions can provide some available evidences for the design and improvement of the vehicle system.     

Application of a modified Lindstedt–Poincaré method in coupled TDOF systems with quadratic nonlinearity and a constant external excitation

An analytical approach is developed for the nonlinear oscillation of a conservative, two-degree-of-freedom (TDOF) mass-spring system with serial combined linear–nonlinear stiffness excited by a constant external force. The main idea of the proposed approach lies in two categories, the first one is the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation. Another is the treatment a quadratic nonlinear oscillator (QNO) by the modified Lindstedt–Poincaré (L-P) method presented recently by the authors. The first-order and second-order analytical approximations for the modified L-P method are established for the QNOs with satisfactory results. After solving the nonlinear differential equation, the displacements of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, the modified L-P method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and classical harmonic balance methods. Two examples of nonlinear TDOF mass-spring systems excited by a constant external force are selected and the approximate solutions are verified with the exact solutions derived from the Jacobi elliptic function and also the numerical fourth-order Runge–Kutta solutions.     

Nonlinear dynamics of an inclined beam subjected to a moving load

In this paper, the nonlinear dynamic response of an inclined pinned-pinned beam with a constant cross section, finite length subjected to a concentrated vertical force traveling with a constant velocity is investigated. The study is focused on the mode summation method and also on frequency analysis of the governing PDEs equations of motion. Furthermore, the steady-state response is studied by applying the multiple scales method. The nonlinear response of the beam is obtained by solving two coupled nonlinear PDEs governing equations of planar motion for both longitudinal and transverse oscillations of the beam. The dynamic magnification factor and normalized time histories of mid-pint of the beam are obtained for various load velocity ratios and the outcome results have been illustrated and compared to the results with those obtained from traditional linear solution. The appropriate parametric study considering the effects of the linear viscous damping, the velocity of the traveling load, beam inclination angle under zero or nonzero axial load are carried out to capture the influence of the effect of large deflections caused by stretching effects due to the beam’s immovable ends. It was seen that quadratic nonlinearity renders the softening effect on the dynamic response of the beam under the act of traveling load. Also in the case where the object leaves the inclined beam, its planar motion path is derived and the targeting accuracy is investigated and compared with those from the rigid solution assumption. Moreover, the stability analysis of steady-state response for the modes equations having quadratic nonlinearity was carried out and it was observed from the frequency response curves that for the considered parameters in the case of internal-external primary resonance, both saturation phenomenon and jump phenomenon can be predicted for the longitudinal excitation.     

Evolution of the double-jumping in pipes conveying fluid flowing at the supercritical speed

Non-linear vibration of viscoelastic pipes conveying fluid around curved equilibrium due to the supercritical flow is investigated with the emphasis on steady-state response in external and internal resonances. The governing equation, a non-linear integro-partial-differential equation, is truncated into a perturbed gyroscopic system via the Galerkin method. The method of multiple scales is applied to establish the solvability condition in the first primary resonance and the 2:1 internal resonance. The approximate analytical expressions are derived for the frequency–amplitude curves of the steady-state responses. The stabilities of the steady-state responses are determined. The generation and the vanishing of a double-jumping phenomenon on the frequency–amplitude curves are examined. The analytical results are supported by the numerical integration results.     

Nonlinear dynamic analysis of cantilevered piezoelectric energy harvesters under simultaneous parametric and external excitations

The nonlinear dynamics of cantilevered piezoelectric beams is investigated under simultaneous parametric and external excitations. The beam is composed of a substrate and two piezoelectric layers and assumed as an Euler–Bernoulli model with inextensible deformation. A nonlinear distributed parameter model of cantilevered piezoelectric energy harvesters is proposed using the generalized Hamilton's principle. The proposed model includes geometric and inertia nonlinearity, but neglects the material nonlinearity. Using the Galerkin decomposition method and harmonic balance method, analytical expressions of the frequency–response curves are presented when the first bending mode of the beam plays a dominant role. Using these expressions, we investigate the effects of the damping, load resistance, electromechanical coupling, and excitation amplitude on the frequency–response curves. We also study the difference between the nonlinear lumped-parameter and distributedparameter model for predicting the performance of the energy harvesting system. Only in the case of parametric excitation, we demonstrate that the energy harvesting system has an initiation excitation threshold below which no energy can be harvested.We also illustrate that the damping and load resistance affect the initiation excitation threshold.