Analytical approximation to large-amplitude oscillation of a non-linear conservative system

This paper deals with non-linear oscillation of a conservative system having inertia and static non-linearities. By combining the linearization of the governing equation with the method of harmonic balance, we establish analytical approximate solutions for the non-linear oscillations of the system. Unlike the classical harmonic balance method, linearization is performed prior to proceeding with harmonic balancing, thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations. Hence, we are able to establish analytical approximate formulas for the exact frequency and periodic solution. These analytical approximate formulas show excellent agreement with the exact solutions, and are valid for small as well as large amplitudes of oscillation.     

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.     

An analytical approximate technique for a class of strongly non-linear oscillators

An analytical approximate technique for large amplitude oscillations of a class of conservative single degree-of-freedom systems with odd non-linearity is proposed. The method incorporates salient features of both Newton's method and the harmonic balance method. 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 non-linear algebraic equations without analytical solution. With carefully constructed iterations, only a few iterations can provide very accurate analytical approximate solutions for the whole range of oscillation amplitude beyond the domain of possible solution by the conventional perturbation methods or harmonic balance method. Three examples including cubic-quintic Duffing oscillators are presented to illustrate the usefulness and effectiveness of the proposed technique.     

A New Method for Approximate Analytical Solutions to Nonlinear Oscillations of Nonnatural Systems

This paper deals with nonlinear oscillations of a conservative,nonnatural, single-degree-of-freedom system with odd nonlinearity. Bycombining the linearization of the governing equation with the method ofharmonic balance, we establish approximate analytical solutions for thenonlinear oscillations of the system. Unlike the classical harmonicbalance method, the linearization is performed prior to proceeding withharmonic balancing thus resulting in linear algebraic equations insteadof nonlinear algebraic equations. Hence, we are able to establish theapproximate analytical formulas for the exact period and periodicsolution. These approximate solutions are valid for small as well aslarge amplitudes of oscillation. Two examples are presented toillustrate that the proposed formulas can give excellent approximateresults.     

On the harmonic balance with linearization for asymmetric single degree of freedom non-linear oscillators

This paper deals with analytical approximation of non-linear oscillations of conservative asymmetric single degree of freedom systems, using the method of harmonic balance with linearization. This technique which consists of linearizing the governing equations prior to harmonic balance permits us to avoid solving complicated non-linear algebraic equations. But it could be applied only to symmetric oscillations for which it proves to be very simple and effective. This restriction is due to the fact that the method requires an appropriate initial approximate solution as input. Such a solution could not be readily identified for nonsymmetric oscillations, contrary the symmetric case where the fundamental harmonic works well. For these nonsymmetric oscillations, we propose in this paper to consider an initial approximation which consists of a small bias plus the fundamental harmonic. By expanding the corresponding harmonic balance equations respectively to first and second order in the bias, we are able to easily determine the bias and thus the required initial approximate solution that yields consistent solution at higher order. We use three examples to illustrate the proposed approach and reveal its simplicity and its very good convergence.     

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.     

Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

A Frequency Method for Predicting Limit Cycle Bifurcations

The paper studies the bifurcations of limit cycles in a rather general class of nonlinear dynamic systems. Relying on the classical harmonic balance approach as applied in control engineering neat frequency conditions for such bifurcations are derived. These results, approximate in nature, make clear the structural mechanism of the considered phenomena and can be applied to predict the occurrence of bifurcations as a function of system parameters. The application to several examples of different complexity shows the simplicity and accuracy of the proposed method for solving complicated problems of nonlinear dynamics.     

Construction of Basis Functions and Their Application to Boundary-Value Problems of Mechanics of Continuous Media

A unified method for constructing basis (eigen) functions is proposed to solve problems of mechanics of continuous media, problems of cubature and quadrature, and problems of approximation of hypersurfaces. Numerical-analytical methods are described, which allow obtaining approximate solutions of internal and external boundary-value problems of mechanics of continuous media of a certain class (both linear and nonlinear). The method is based on decomposition of the sought solutions of the considered partial differential equations into series in basis functions. An algorithm is presented for linearization of partial differential equations and reduction of nonlinear boundary-value problems, which are reduced to systems of linear algebraic equations with respect to unknown coefficients without using traditional methods of linearization.     

电动力学中衰减机械系统的非线性振荡方法

衰减机械系统的非线性振荡可用来研究长约瑟夫逊结的电动力学方程式，而这方程式等同于弱衰减机械系统的非线性振荡。本文应用的方法是将控制方程线性化及结合谐波平衡法（线性谐波平衡法）而产生色散关系，再把平均法应用在弱非线性的耗散系统中得到非常准确的瞬变反应。在此提出的方法不仅考虑能量耗散，而且利用简单的线性代数等式关系来代替冗长及复杂的分析近似解。     

一种新的求解非线性系统周期解方法

本文利用切比雪夫多项式的若干良好性质,对非自治强非线性动力系统进行分析。将状态矢量在主周期上先展开谐波级数的形式,再将各谐波展开为切比雪夫多项式的形式,从而将求周期解的问题转变为非线性代数方程组的求解问题,得出一种可以方便、迅速地获得近似周期解的解析方法。此方法不依赖于小参数假设,可以用于分析强非线性问题和高维问题,而且对参数激励系统同样有效。以Duffing系统周期解的计算为例,通过与标准谐波平衡方法和四阶Runge-Kutta数值积分结果比较,说明此方法的有效性。     

Euler杆大挠度屈曲的解析逼近

研究Euler杆大挠度屈曲问题，将控制方程的线性化与谐波平衡法组合起来，分别建立以杆端转角形式表示的屈曲荷载及最大挠度的解析逼近公式，这些公式既适用于小变形又适用于大变形。     

Approximate analytical solution in slow-fast system based on modified multi-scale method

A simple, yet accurate modified multi-scale method (MMSM) for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed. This method depends on the classical multi-scale method (MSM) and the method of variation of parameters. Assuming that the forced excitation is a constant, one could easily obtain the approximate analytical solution of the simplified system based on the traditional MSM. Then, this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation. To certify the correctness and precision of the proposed analytical method, the van der Pol system with two scales subject to slowly periodic excitation is investigated; this system presents rich dynamical phenomena such as spiking (SP), spiking-quiescence (SP-QS), and quiescence (QS) responses. The approximate analytical expressions of the three types of responses are given by the MMSM, and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method (HBM). The results obtained by the present method are considerably better than those obtained by traditional methods, quantitatively and qualitatively, particularly when the excitation frequency is far less than the natural frequency of the system.     

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.     

Adaptive neural DSC for stochastic nonlinear constrained systems under arbitrary switchings

This paper focuses on the primary resonance analysis of a dual-rotor system having two rotor unbalance excitations of different rotating speeds and being connected by an inter-shaft ball bearing. Due to the complex excitation condition and the complicated nonlinear bearing forces of the inter-shaft bearing, the general analytical methods, e.g., the multiple scales method or the harmonic balance method, are failed to give the theoretical solutions. Thus, the harmonic balance–alternating frequency/time domain (HB–AFT) method is formulated to deal with this problem. The basic idea of the method is using the inverse discrete Fourier transform and the discrete Fourier transform, instead of the direct analytical relationship between the supposed solutions of the system and the nonlinear forces, to construct the harmonic expressions of the nonlinear forces, which is the so-called alternating frequency/time domain technique. By using the HB–AFT method, therefore, a Newton– Raphson iteration procedure can be performed to demonstrate the approximate solutions of the system. Accordingly, the frequency responses of the system affected by some critical parameters, such as rotating speed ratio, unbalances of both the inner and outer rotors, and clearance of the inter-shaft bearing, are analyzed, respectively. A variety of phenomena including double resonance peaks, biperiodic and quasi-periodic behaviors, and resonance hysteresis phenomenon are obtained, which are discussed in detail through diagrams for separated frequency responses with different frequency components and rotors’ orbits for different combinations of system parameters. Most prominently, for a relatively small unbalance of rotor as well as a relatively large clearance of the inter-shaft bearing, the resonance hysteresis phenomena are more obvious. The results obtained are also compared with the direct numerical simulation results, and the comparisons show good agreements. In addition, the methodology formulated in this paper is a general approach, which can be applied to other engineering systems with multi-frequency excitations.     

Approximate analytical solution in slow-fast system based on modi?ed multi-scale method

A simple, yet accurate modi?ed multi-scale method(MMSM) for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed. This method depends on the classical multi-scale method(MSM) and the method of variation of parameters. Assuming that the forced excitation is a constant, one could easily obtain the approximate analytical solution of the simpli?ed system based on the traditional MSM. Then, this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation. To certify the correctness and precision of the proposed analytical method, the van der Pol system with two scales subject to slowly periodic excitation is investigated; this system presents rich dynamical phenomena such as spiking(SP),spiking-quiescence(SP-QS), and quiescence(QS) responses. The approximate analytical expressions of the three types of responses are given by the MMSM, and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method(HBM). The results obtained by the present method are considerably better than those obtained by traditional methods,quantitatively and qualitatively, particularly when the excitation frequency is far less than the natural frequency of the system.     

建立在同伦基础之上的一种非线性分析方法(2)——单自由度守恒系统周期的通用公式

本文以非线性的单自由度守恒系统为例,描述了一种新的非线性分析方法。此方法建立在光滑映射的基础上,不依赖于小参数,从而特别适合于求解不具有小参数的强非线性问题。本文所给近似公式的精确度对参数值的变化不敏感,对任意大小的参数都是一致有效的。     

Forced response of quadratic nonlinear oscillator: comparison of various approaches

The primary resonances of a quadratic nonlinear system under weak and strong external excitations are investigated with the emphasis on the comparison of different analytical approximate approaches. The forced vibration of snap-through mechanism is treated as a quadratic nonlinear oscillator. The Lindstedt-Poincaré method, the multiple-scale method, the averaging method, and the harmonic balance method are used to determine the amplitude-frequency response relationships of the steady-state responses. It is demonstrated that the zeroth-order harmonic components should be accounted in the application of the harmonic balance method. The analytical approximations are compared with the numerical integrations in terms of the frequency response curves and the phase portraits. Supported by the numerical results, the harmonic balance method predicts that the quadratic nonlinearity bends the frequency response curves to the left. If the excitation amplitude is a second-order small quantity of the bookkeeping parameter, the steady-state responses predicted by the second-order approximation of the LindstedtPoincaré method and the multiple-scale method agree qualitatively with the numerical results. It is demonstrated that the quadratic nonlinear system implies softening type nonlinearity for any quadratic nonlinear coefficients.     

Localized Basis Function Method for Computing Limit Cycle Oscillations

An alternate approach to the standard harmonic balance method (based on Fourier transforms) is proposed. The proposed method begins with an idea similar to the harmonic balance method, i.e. to transform the initial set of differential equations of the dynamics to a set of discrete algebraic equations. However, as distinct from previous harmonic balance techniques, the proposed method uses a set of basis functions which are localized in time and are not necessarily sinusoidal. Also as distinct from previous harmonic balance methods, the algebraic equations obtained after the transformation of the differential equations of the dynamics are solved in the time domain rather than the frequency domain. Numerical examples are provided to demonstrate the performance of the method for autonomous and forced dynamics of a Van der Pol oscillator.     

一种单元谐波平衡法

基于有限元离散,对于工程中的非线性响应问题,提出一种单元谐波平衡法．与常规的谐波平衡法不同,本文将谐波平衡方程建立在有限元素上,从而兼顾了有限元素法和常规谐波平衡法两大优势．有限元技术的应用能使得求解问题的范围扩大到复杂工程结构,而谐波平衡概念的使用将使得含有复杂变形和复杂本构关系的动力学响应问题得到有效解决．所提方法能适用于工程结构中具有复杂非线性关系的动力学响应问题．由于谐波平衡法的实施依赖于谐波系数方程及其切线刚度矩阵的解析推导,尽管已经局限到有限元素上,但对于较为复杂一些的本构关系,推导仍非易事．为解决这些问题,放弃通常对于变形梯度和应变张量所作的向量假设,而是从连续介质力学中基本的几何关系入手,提出一种矩阵分解形式．通过利用张量的内蕴导数定义以及关于迹函数的有关性质,给出应力增量的一种新的表现形式．当它与变形梯度的矩阵分解相结合时,使得切线刚度矩阵的导出变得十分简单,而且所得计算形式也比通常紧凑和方便许多．