Global residue harmonic balance method to periodic solutions of a class of strongly nonlinear oscillators

A new approach, namely the global residue harmonic balance method, was advanced to determine the accurate analytical approximate periodic solution of a class of strongly nonlinear oscillators. A class of nonlinear jerk equation containing velocity-cubed and velocity times displacements-squared was taken as a typical example. Unlike other harmonic balance methods, all the former residual errors are introduced in the present approximation to improve the accuracy. Comparison of the result obtained using this approach with the exact one and simplicity and efficiency of the proposed procedure. The method can be easily extended to other strongly nonlinear oscillators.     

Periodic wave solutions of coupled integrable dispersionless equations by residue harmonic balance

We introduce the residue harmonic balance method to generate periodic solutions for nonlinear evolution equations. A PDE is firstly transformed into an associated ODE by a wave transformation. The higher-order approximations to the angular frequency and periodic solution of the ODE are obtained analytically. To improve the accuracy of approximate solutions, the unbalanced residues appearing in harmonic balance procedure are iteratively considered by introducing an order parameter to keep track of the various orders of approximations and by solving linear equations. Finally, the periodic solutions of PDEs result. The proposed method has the advantage that the periodic solutions are represented by Fourier functions rather than the sophisticated implicit functions as appearing in most methods.     

Forward residue harmonic balance for autonomous and non-autonomous systems with fractional derivative damping

Both the autonomous and non-autonomous systems with fractional derivative damping are investigated by the harmonic balance method in which the residue resulting from the truncated Fourier series is reduced iteratively. The first approximation using a few Fourier terms is obtained by solving a set of nonlinear algebraic equations. The unbalanced residues due to Fourier truncation are considered iteratively by solving linear algebraic equations to improve the accuracy and increase the number of Fourier terms of the solutions successively. Multiple solutions, representing the occurrences of jump phenomena, supercritical pitchfork bifurcation and symmetry breaking phenomena are predicted analytically. The interactions of the excitation frequency, the fractional order, amplitude, phase angle and the frequency amplitude response are examined. The forward residue harmonic balance method is presented to obtain the analytical approximations to the angular frequency and limit cycle for fractional order van der Pol oscillator. Numerical results reveal that the method is very effective for obtaining approximate solutions of nonlinear systems having fractional order derivatives.     

Nonlinear analysis of a parametrically excited beam with intermediate support by using Multi-dimensional incremental harmonic balance method

In this paper, a nonlinear Euler-Bernoulli beam under a concentrated harmonic excitation with intermediate nonlinear support is investigated. Continuous expression for the kinetic energy, potential energy and dissipation function are constructed. An energy method based on the Lagrange equation combined with the Galerkin truncation is used for discretizing the governing equation. The Multi-dimensional incremental harmonic balance method (MIHBM) is derived, and the comparisons between the numerical results and the approximate analytical solutions based on the MIHBM verify the excellent accuracy of the MIHBM. The steady state dynamic of the beam is investigated by MIHBM. In order to investigate the energy transmission and understand the vibration response of the Euler-Bernoulli beam, the effects of the key parameters on the dynamic behaviors are studied and discussed, individually. The results show that the amplitude-frequency curves exhibits softening nonlinear behavior in the super-harmonic resonance region, and near resonant region the hardening nonlinear behavior is observed depending on the different parameters. Nonlinear dynamic analysis, such as bifurcation, 3-D frequency spectrum, waveform, frequency spectrum, phase diagram and Poincaré map, are also presented in order to study the influences of the key parameters on the vibration behaviors for the beam in a more accurate manner. In addition, the path to chaotic motion is observed to be through a sequence of the periodic motion and quasi-periodic motion.     

Study of non-linear dynamic behavior of open cracked functionally graded Timoshenko beam under forced excitation using harmonic balance method in conjunction with an iterative technique

The nonlinear modeling and subsequent dynamic analysis of cracked Timoshenko beams with functionally graded material (FGM) properties is investigated for the first time using harmonic balance method followed by an iterative technique. Crack is assumed to be open throughout. During modeling, nonlinear strain–displacement relation is considered. Rotational spring model is adopted in order to model the open cracks. Energy formulations are established using Timoshenko beam theory. Nonlinear governing differential equations of motion are derived using Lagrange's equation. In order to incorporate the influence of higher order harmonics, harmonic balance method is employed. This reduces the governing differential equations into nonlinear set of algebraic equations. These equations are solved using two different iterative techniques. Methodology is computationally easier and efficient as well. This is observed that although assumption of simple harmonic motion (SHM) simplifies the problem, it yields to erroneous results at higher amplitude of motion. However, accuracy of the solution is improved considerably when the contribution of higher order harmonic terms are considered in the analysis. Results are compared with the available results, which confirm the validity of the methodology. Subsequent to that a parametric study on influence of forcing term, material indices and crack parameters on large amplitude vibration of Timoshenko beams is performed for two different boundary conditions.     

Nonlinear analysis of buckling,free vibration and dynamic stability for the piezoelectric functionally graded beams in thermal environment

Employing Euler–Bernoulli beam theory and the physical neutral surface concept, the nonlinear governing equation for the functionally graded material beam with two clamped ends and surface-bonded piezoelectric actuators is derived by the Hamilton’s principle. The thermo-piezoelectric buckling, nonlinear free vibration and dynamic stability for the piezoelectric functionally graded beams, subjected to one-dimensional steady heat conduction in the thickness direction, are studied. The critical buckling loads for the beam are obtained by the existing methods in the analysis of thermo-piezoelectric buckling. The Galerkin’s procedure and elliptic function are adopted to obtain the analytical solution of the nonlinear free vibration, and the incremental harmonic balance method is applied to obtain the principle unstable regions of the piezoelectric functionally graded beam. In the numerical examples, the good agreements between the present results and existing solutions verify the validity and accuracy of the present analysis and solving method. Simultaneously, validation of the results achieved by rule of mixture against those obtained via the Mori–Tanaka scheme is carried out, and excellent agreements are reported. The effects of the thermal load, electric load, and thermal properties of the constituent materials on the thermo-piezoelectric buckling, nonlinear free vibration, and dynamic stability of the piezoelectric functionally graded beam are discussed, and some meaningful conclusions have been drawn.     

Bifurcating periodic solution for a class of first-order nonlinear delay differential equation after Hopf bifurcation

In this paper, we predict the accurate bifurcating periodic solution for a general class of first-order nonlinear delay differential equation with reflectional symmetry by constructing an approximate technique, named residue harmonic balance. This technique combines the features of the homotopy concept with harmonic balance which leads to easy computation and gives accurate prediction on the periodic solution to the desired accuracy. The zeroth-order solution using just one Fourier term is applied by solving a set of nonlinear algebraic equations containing the delay term. The unbalanced residues due to Fourier truncation are considered iteratively by solving linear equations to improve the accuracy and increase the number of Fourier terms of the solutions successively. It is shown that the solutions are valid for a wide range of variation of the parameters by two examples. The second-order approximations of the periodic solutions are found to be in excellent agreement with those obtained by direct numerical integration. Moreover, the residue harmonic balance method works not only in determining the amplitude but also the frequency of the bifurcating periodic solution. The method can be easily extended to other delay differential equations.     

The iterative homotopy harmonic balance method for conservative Helmholtz-Duffing oscillators

An approach that the iterative homotopy harmonic balance method which incorporates salient features of both the parameter-expansion and the harmonic balance is presented to solve conservative Helmholtz-Duffing oscillators. Since the behaviors of the solutions in the positive and negative directions are quite different, the asymmetric equation is separated into two auxiliary equations. The auxiliary equations are solved by proposed method. The results show it works very well for the whole range of initial amplitudes in a variety of cases, and the excellent agreement of the approximate periods and periodic solutions with exact ones have been demonstrated and discussed. And, the proposed method is very simple in its principle and has great potential to be applied to other nonlinear oscillators.     

Elliptic harmonic balance method for two degree-of-freedom self-excited oscillators

Elliptic harmonic balance (EHB) method as an analytical method is widely used for strongly non-linear oscillators with a single degree-of-freedom (DOF). The oscillation equations are expressed by elliptic functions, and then the expressions are expanded as harmonics of elliptic type while only the first harmonic is retained. To the best of our knowledge, however, it seems that the EHB method has not been found applications in two or multiple DOFs systems. One possible reason is that the EHB method may cause a difficult problem that the number of equations obtained by harmonic balancing is not equal to that of unknowns. To this end, in the present paper, the EHB method is therefore extended to study a class of strongly self-excited oscillators with two DOFs. Prior to harmonic balancing, an additional equation is established to tackle the aforementioned problem. Illustrative examples show that solutions of limit cycles obtained by the proposed method are in good agreement with the numerical solutions.     

Exact front, soliton and hole solutions for a modified complex Ginzburg-Landau equation from the harmonic balance method

A novel approach of using harmonic balance (HB) method is presented to find front, soliton and hole solutions of a modified complex Ginzburg-Landau equation. Three families of exact solutions are obtained, one of which contains two parameters while the others one parameter. The HB method is an efficient technique in finding limit cycles of dynamical systems. In this paper, the method is extended to obtain homoclinic/heteroclinic orbits and then coherent structures. It provides a systematic approach as various methods may be needed to obtain these families of solutions. As limit cycles with arbitrary value of bifurcation parameter can be found through parametric continuation, this approach can be extended further to find analytic solution of complex quintic Ginzburg-Landau equation in terms of Fourier series.     

Application of the harmonic balance method to ground moling machines operating in periodic regimes

A new system for ground moling has been patented by the University of Aberdeen and licensed world-wide. This new system is based on vibro-impact dynamics and offers significant advantages over existing systems in terms of penetrative capability and reduced soil disturbance. This paper describes current research into the mathematical modelling of the system. Periodic response is required to achieve the optimal penetrating conditions for the ground moling process, as this results in reduced soil penetration resistance. Therefore, there is a practical need for a robust and efficient methodology to calculate periodic responses for a wide range of operational parameters. Due to the structural complexity of a real vibro-impact moling system, the dynamic response of an idealised impact oscillator has been investigated in the first instance. This paper presents a detailed study of periodic responses of the impact oscillator under harmonic forcing using the alternating frequency-time harmonic balance method. Recommendations of how to effectively adapt the alternating frequency-time harmonic balance method for a stiff impacting system are given.     

Relationship between the homotopy analysis method and harmonic balance method

This paper presents a study of the relationship between the homotopy analysis method (HAM) and harmonic balance (HB) method. The HAM is employed to obtain periodic solutions of conservative oscillators and limit cycles of self-excited systems, respectively. Different from the usual procedures in the existing literature, the HAM is modified by retaining a given number of harmonics in higher-order approximations. It is proved that as long as the solution given by the modified HAM is convergent, it converges to one HB solution. The Duffing equation, the van der Pol equation and the flutter equation of a two-dimensional airfoil are taken as illustrations to validate the attained results.     

Incremental harmonic balance method for nonlinear flutter of an airfoil with uncertain-but-bounded parameters

The limit cycle oscillation of a two-dimensional airfoil with parameter variability in an incompressible flow is investigated using the incremental harmonic balance (IHB) method. The variable parameters, such as the wind speed, the cubic plunge and pitch stiffness coefficients, are modeled as either bounded uncertain or stochastic parameters. In the solution process of the IHB method, the bounded parameters are considered as an active increment. Taking all values over the considered bounded regions of the active parameters provides us with a series of IHB solutions of limit cycle oscillations of the airfoil. With the aid of the attained solutions, the bounds and some statistical properties of the limit cycle oscillations are determined and compared with Monte Carlo simulation (MCS) results. Numerical examples show that the proposed approach is valid and effective for analyzing strongly nonlinear vibration problems with bounded uncertainties.     

Analysis and modeling of a flexible rectangular cantilever plate

Flexible plate structures with large deflection and rotation are commonly used structures in engineering. How to analyze and solve the cantilever plate with large deflection and rotation is still an unsolved problem. In this paper, a general nonlinear flexible rectangular cantilever plate considering large deflection and rotation angle is modeled, solved and analyzed. Hamilton's principle is applied to obtain the nonlinear differential dynamic equations and boundary conditions by introducing a coordinate transformation between the Cartesian coordinate system and the deformed local coordinate system. Stress function relating to in-plane force resultants and shear forces is given for the first time for complex coupling equations caused by coordinate transformation. The nonlinear equations and the solving method are validated by experiments. Then, harmonic balance method is adopted to get the nonlinear frequency-response curves, which shows strong hardening spring characteristic of this system. Runge–Kutta methods are used to reveal complex nonlinear behaviors such as 5 super-harmonic resonance, bifurcations and chaos for general nonlinear flexible rectangular cantilever plate.     

Newton–harmonic balancing approach for accurate solutions to nonlinear cubic–quintic Duffing oscillators

This paper presents a new approach for solving accurate approximate analytical higher-order solutions for strong nonlinear Duffing oscillators with cubic–quintic nonlinear restoring force. The system is conservative and with odd nonlinearity. The new approach couples Newton’s method with harmonic balancing. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton’s method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of nonlinear algebraic equations without analytical solution. Using the approach, accurate higher-order approximate analytical expressions for period and periodic solution are established. These approximate solutions are valid for small as well as large amplitudes of oscillation. In addition, it is not restricted to the presence of a small parameter such as in the classical perturbation method. Illustrative examples are presented to verify accuracy and explicitness of the approximate solutions. The effect of strong quintic nonlinearity on accuracy as compared to cubic nonlinearity is also discussed.     

The extended homogeneous balance method and its applications for a class of nonlinear evolution equations

The extended homogeneous balance method with the aid of computer algebraic system Maple, is proposed for seeking the travelling wave solutions for a class of nonlinear evolution equations, in which the homogeneous balance method is applied to solve the Riccati equation and the reduced nonlinear evolution equations, respectively. Many new exact travelling wave solutions are successfully obtained. The method is straightforward and concise, and it can be also applied to other nonlinear evolution equations.     

Oscillatory region and asymptotic solution of fractional van der Pol oscillator via residue harmonic balance technique

In this paper, a powerfully analytical technique is proposed for predicting and generating the steady state solution of the fractional differential system based on the method of harmonic balance. The zeroth-order approximation using just one Fourier term is applied to predict the parametric function for the boundary between oscillatory and non-oscillatory regions of the fractional van der Pol oscillator. The unbalanced residues due to Fourier truncation are considered iteratively by solving linear algebraic equations to improve the accuracy of the solutions successively. The highly accurate solutions to the angular frequency and limit cycle of fractional van der Pol oscillator are obtained and compared. The results reveal that the technique described in this paper is very effective and simple for obtaining asymptotic solution of nonlinear system having fractional order derivative.     

基于时域配点法的充液储箱系统多模态方程求解

采用时域配点法研究了充液储箱系统多模态方程的稳态周期解.在模型求解过程中,利用牛顿迭代法求解了配点法得到的非线性代数方程组,而牛顿迭代的初值来自谐波平衡法求解得到的低阶谐波近似.数值仿真结果验证了时域配点法的有效性,并验证以二倍激励频率为基频的第二模态的假设形式更为有效.最终通过对比谐波系数数量级提出一种更为简洁有效的模态表达形式.     

Multi-harmonic measurements and numerical simulations of nonlinear vibrations of a beam with non-ideal boundary conditions

This study presents a direct comparison of measured and predicted nonlinear vibrations of a clamped–clamped steel beam with non-ideal boundary conditions. A multi-harmonic comparison of simulations with measurements is performed in the vicinity of the primary resonance. First of all, a nonlinear analytical model of the beam is developed taking into account non-ideal boundary conditions. Three simulation methods are implemented to investigate the nonlinear behavior of the clamped–clamped beam. The method of multiple scales is used to compute an analytical expression of the frequency response which enables an easy updating of the model. Then, two numerical methods, the Harmonic Balance Method and a time-integration method with shooting algorithm, are employed and compared one with each other. The Harmonic Balance Method enables to simulate the vibrational stationary response of a nonlinear system projected on several harmonics. This study then proposes a method to compare numerical simulations with measurements of all these harmonics. A signal analysis tool is developed to extract the system harmonics’ frequency responses from the temporal signal of a swept sine experiment. An evolutionary updating algorithm (Covariance Matrix Adaptation Evolution Strategy), coupled with highly selective filters is used to identify both fundamental frequency and harmonic amplitudes in the temporal signal, at every moment. This tool enables to extract the harmonic amplitudes of the output signal as well as the input signal. The input of the Harmonic Balance Method can then be either an ideal mono-harmonic signal or a multi-harmonic experimental signal. Finally, the present work focuses on the comparison of experimental and simulated results. From experimental output harmonics and numerical simulations, it is shown that it is possible to distinguish the nonlinearities of the clamped–clamped beam and the effect of the non-ideal input signal.     

Observing nonlinear stochastic resonance with piecewise constant driving forces by the method of moments

The original method of moments confined within linear response theory is improved to calculate the nonlinear dynamic response of the standard noisy bistable stochastic systems in the general response sense by proposing a different operating technique. Especially, the proposed technique is simple and efficient to be used to the cases where the driving forces are not harmonics. Using the piecewise constant driving forces for demonstration, our comparative analysis shows that the long time ensemble average and the first three harmonic susceptibilities calculated by the proposed technique are of high accuracy. The dependence of the spectral amplification parameters at the first three harmonics on the noise intensity is also investigated, and the analysis to the resonant curves suggests a possible way to induce the even-order harmonic stochastic resonance.