On the local stability analysis of the approximate harmonic balance solutions

Periodic response of nonlinear oscillators is usually determined by approximate methods. In the "steady state" type methods, first an approximate solution for the steady state periodic response is determined, and then the local stability of this solution is determined by analyzing the equation of motion linearized about this predicted "solution". An exact stability analysis of this linear variational equation can provide erroneous stability type information about the approximate solutions. It is shown that a consistent stability type information about these solutions can be obtained only when the linearized variational equation is analyzed by approximate methods, and the level of accuracy of this analysis is consistent with that of the approximate solutions. It is demonstrated that these consistent stability results do not imply that the approximate solution is qualitatively correct. It is also shown that the difference between an approximate and the next higher order stability analysis can be used to "guess" the role of higher harmonics in the periodic response. This trial and error procedure can be used to ensure the qualitatively correct and numerically accurate nature of the approximate solutions and the corresponding stability analysis.     

Accurate Higher-Order Analytical Approximate Solutions to Large-Amplitude Oscillating Systems with a General Non-Rational Restoring Force

A new approximate analytical approach for accurate higher-order nonlinear solutions of oscillations with large amplitude is presented in this paper. The oscillatory system is subjected to a non-rational restoring force. This approach is built upon linearization of the governing dynamic equation associated with the method of harmonic balance. Unlike the classical harmonic balance method, simple linear algebraic equations instead of nonlinear algebraic equations are obtained upon linearization prior to harmonic balancing. This approach also explores large parameter regions beyond the classical perturbation methods which in principle are confined to problems with small parameters. It has significant contribution as there exist many nonlinear problems without small parameters. Through some examples in this paper, we establish the general approximate analytical formulas for the exact period and periodic solution which are valid for small as well as large amplitudes of oscillation.     

An equivalent polynomial approximation technique in nonlinear structural dynamics

Summary A new technique is proposed to obtain an approximate probability density for the response of a general nonlinear system under Gaussian white noise excitations. In this new technique, the original nonlinear system is replaced by another equivalent nonlinear system, structured by the polynomial formula, for which the exact solution of stationary probability density function is obtainable. Since the equivalent nonlinear system structured in this paper originates directly from certain classes of real nonlinear mechanical systems, the technique is applied to some very challenging nonlinear systems in order to show its power and efficiency. The calculated results show that applying the technique presented here can yield exact stationary solutions for the nonlinear oscillators. This is obtained by using an energy-dependent system, and for a nonlinearity of a more complex type. A more accurate approximate solution is then available, and is compared with the approximation. Application of the technique is illustrated by examples.     

Geometrically exact curved beam element using internal force field defined in deformed configuration

This paper presents a study on the development of high-performance finite elements for geometrically nonlinear analysis of frame structures with curved members. Based on the geometrically exact beam theory, a highly efficient and accurate mixed finite element is developed. A new approach is proposed for constructing the independent internal force field by including major terms satisfying equilibrium conditions in the deformed configuration. An element-level equilibrium iteration procedure is employed for the condensation of element internal degrees of freedom during the nonlinear solution. Numerical results are presented to demonstrate the excellent performance of the element developed, and it is shown that even when each structural member is modelled with just one element, accurate solutions can still be achieved.     

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.     

A novel construction of homoclinic and heteroclinic orbits in?nonlinear oscillators by a perturbation-incremental method

A novel construction of homoclinic/heteroclinic orbits (HOs) in nonlinear oscillators is presented in this paper. An accurate analytical solution of a HO for small perturbation can be obtained in terms of trigonometric functions. An advantage of the present construction is that it gives an accurate approximate solution of a HO for large parametric value in relatively few harmonic terms while other analytical methods such as the Lindstedt?CPoincaré method and the multiple scales method fail to do so.     

Analysis of a laminar boundary layer flow over a flat plate with injection or suction

An analysis is performed to study a laminar boundary layer flow over a porous flat plate with injection or suction imposed at the wall. The basic equations of this problem are reduced to a system of nonlinear ordinary differential equations by means of appropriate transformations. These equations are solved analytically by the optimal homotopy asymptotic method (OHAM), and the solutions are compared with the numerical solution (NS). The effect of uniform suction/injection on the heat transfer and velocity profile is discussed. A constant surface temperature in thermal boundary conditions is used for the horizontal flat plate.     

Darcy-Forchheimer flow with nonlinear mixed convection

An analysis of the mixed convective flow of viscous fluids induced by a nonlinear inclined stretching surface is addressed. Heat and mass transfer phenomena are analyzed with additional effects of heat generation/absorption and activation energy, respectively. The nonlinear Darcy-Forchheimer relation is deliberated. The dimensionless problem is obtained through appropriate transformations. Convergent series solutions are obtained by utilizing an optimal homotopic analysis method (OHAM). Graphs depicting the consequence of influential variables on physical quantities are presented. Enhancement in the velocity is observed through the local mixed convection parameter while an opposite trend of the concentration field is noted for the chemical reaction rate parameter.     

Nonlinear vibration of stringer shell by means of extended Hamiltonian approach

In this paper, the nonlinear free vibration of a stringer shell is studied. The mathematical model of the string shell, which is the most convenient for frequency analysis, is considered. Due to the geometrical properties of the vibrating shell, strong nonlinearities are evident. Approximate analytical expressions for the nonlinear vibration are provided by introducing the extended version of the Hamiltonian approach. The method suggested in the paper gives the approximate solution for the differential equation with dissipative term for which the Lagrangian exists. The aim of this study is to provide engineers and designers with an easy method for determining the shell nonlinear vibration frequency and nonlinear behavior. The effects of different parameters on the ratio of nonlinear to linear natural frequency of shells are studied. This analytical representation gives excellent approximations to the numerical solutions for the whole range of the oscillation amplitude, reducing the respective error of the angular frequency in comparison with the Hamiltonian approach. This study shows that a first-order approximation of the Hamiltonian approach leads to highly accurate solutions that are valid for a wide range of vibration amplitudes.     

Remarks on the Perturbation Methods in Solving the Second-Order Delay Differential Equations

The paper presents a study on the validity of perturbation methods, suchas the method of multiple scales, the Lindstedt–Poincaré method and soon, in seeking for the periodic motions of the delayed dynamic systemsthrough an example of a Duffing oscillator with delayed velocityfeedback. An important observation in the paper is that the method ofmultiple scales, which has been widely used in nonlinear dynamics, worksonly for the approximate solutions of the first two orders, and givesrise to a paradox for the third-order approximate solutions of delaydifferential equations. The same problem appears when theLindstedt–Poincaré method is implemented to find the third-orderapproximation of periodic solutions for delay differential equations,though it is effective in seeking for any order approximation ofperiodic solutions for nonlinear ordinary differential equations. Apossible explanation to the paradox is given by the results obtained byusing the method of harmonic balance. The paper also indicates thatthese perturbation methods, despite of some shortcomings, are stilleffective in analyzing the dynamics of a delayed dynamic system sincethe approximate solutions of the first two orders already enable one togain an insight into the primary dynamics of the system.     

非线性转子-轴承系统的周期解及近似解析表达式

通过对普通打靶方法进行改造提出一种确定非线性系统周期轨道及周期的新型打靶算法。首先通过改变系统的时间尺度，将非线性系统周期轨道的周期显式地出现在非线性系统的系统方程中，然后对传统打靶法进行改造，将周期也作为一个参数一起参与打靶法的迭代过程，迭代过程包含对周期轨道和周期的求解，迭代过程中的增量通过优化方法选择，从而能迅速确定出系统的周期轨道及其周期。应用所求的结果结合谐波平衡方法求得了非线性系统的周期轨道的近似解析表达式，理论上通过增加谐波的阶数任何精度的周期解都可以得到。最后将该方法应用于非线性转子轴承系统，求出了在某些参数下转子的周期解及其近似解析表达式，通过与四阶Runge-Kutta数值积分结果比较，验证了方法的有效性，计算结果对于转子系统运动的定量控制有重要理论指导意义。     

A Modified Perturbation Technique Depending Upon an Artificial Parameter

In this paper, a modified perturbation method is proposed to search for analytical solutions of nonlinear oscillators without possible small parameters. An artificial perturbation equation is carefully constructed by embedding an artificial parameter, which is used as expanding parameter. It reveals that various traditional perturbation techniques can be powerfully applied in this theory. Some examples, such as the Duffing equation and the van der Pol equation, are given here to illustrate its effectiveness and convenience. The results show that the obtained approximate solutions are uniformly valid on the whole solution domain, and they are suitable not only for weak nonlinear systems, but also for strongly nonlinear systems. In applying the new method, some special techniques have been emphasized for different problems.     

Approximate symmetry group classification for a nonlinear fractional filtration equation of diffusion-wave type

A one-dimensional nonlinear fractional filtration equation with the Riemann–Liouville time-fractional derivative is proposed for modeling fluid flow through a porous medium. This equation is derived under an assumption that the fluid has a fractional equation of state in which the fluid density depends on the time-fractional derivative of pressure. The obtained equation belongs to the diffusion-wave type of equations. A case when the order of fractional differentiation is close to an integer number is considered, and a small parameter is introduced into the fractional filtration equation under consideration. An expansion of the Riemann–Liouville time-fractional derivative into the series with respect to this small parameter is obtained. Using this expansion, a first-order approximation of the derived fractional filtration equation is performed. Next, the problem of approximate Lie point symmetry group classification for this approximate nonlinear filtration equation with a small parameter is studied. It is shown that approximate symmetry groups admitted by different realizations of the approximate filtration equation have much more dimensions than the corresponding exact Lie point symmetry groups admitted by unperturbed fractional diffusion-wave equations. Obtained classification results permit to construct approximate invariant solutions for the considered nonlinear time-fractional filtration equations.     

An equivalent nonlinearization method for strongly nonlinear oscillations

An equivalent nonlinearization method is proposed for the study of certain kinds of strongly nonlinear oscillators. This method is to express the nonlinear restored force of an oscillatory system by a polynomial of degree two or three such that the asymptotic solutions can be derived in terms of elliptic functions. The least squares method is used to determine the coefficients of approximate polynomials. The advantage of present method is that it is valid for relatively large oscillations. As an application, a strongly nonlinear oscillator with slowly varying parameters resulted from free-electron laser is studied in detail. Comparisons are made with other methods to assess the accuracy of the present method.     

Accurate analytical approximate solutions to general strong nonlinear oscillators

A new approach is presented for establishing the analytical approximate solutions to general strong nonlinear conservative single-degree-of-freedom systems. Introducing two odd nonlinear oscillators from the original general nonlinear oscillator and utilizing the analytical approximate solutions to odd nonlinear oscillators proposed by the authors, we construct the analytical approximate solutions to the original general nonlinear oscillator. These analytical approximate solutions are valid for small as well as large oscillation amplitudes. Two examples are presented to illustrate the great accuracy and simplicity of the new approach.     

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.     

Fully-coupled nonlinear analysis of single lap adhesive joints

This paper presents novel closed-form and accurate solutions for the edge moment factor and adhesive stresses for single lap adhesive bonded joints. In the present analysis of single lap joints, both large deflections of adherends and adhesive shear and peel strains are taken into account in the formulation of two sets of nonlinear governing equations for both longitudinal and transverse deflections of adherends. Closed-form solutions for the edge moment factor and the adhesive stresses are obtained by solving the two sets of fully-coupled nonlinear governing equations. Simplified and accurate formula for the edge moment factor is also derived via an approximation process. A comprehensive numerical validation was conducted by comparing the present solutions and those developed by Goland and Reissner, Hart-Smith and Oplinger with the results of nonlinear finite element analyses. Numerical results demonstrate that the present solutions for the edge moment factor (including the simplified formula) and the adhesive stresses appear to be the best as they agree extremely well with the finite element analysis results for all ranges of material and geometrical parameters.     

Generalized Strongly Nonlinear Quasi-Complementarity Problems

ＧＥＮＥＲＡＬＩＺＥＤＳＴＲＯＮＧＬＹＮＯＮＬＩＮＥＡＲＱＵＡＳＩ－ＣＯＭＰＬＥＭＥＮＴＡＲＩＴＹＰＲＯＢＬＥＭＳＬｉＨｏｎｇ－ｍｅｉ（李红梅）ＤｉｎｇＸｉｅ－ｐｉｎｇ（丁协平）（ＳｉｃｈｕａｎＮｏｒｍａｌＵｎｉｖｅｒｓｉｔｙ），Ｃｈｅｎｇｄｕ（Ｒｅ...     

Hyperbolic conservation laws with umbilic degeneracy,I

In this paper a compactness framework for approximate solutions to nonlinear hyperbolic systems with umbilic degeneracy is established by combining techniques of compensated compactness with some classical methods, and by a detailed analysis of a highly singular equation of Euler-Poisson-Darboux type. Then this framework is successfully applied to prove the convergence of the viscosity method and to prove the existence of global entropy solutions for the Cauchy problem with large initial data for a canonical class of the systems with quadratic flux form.     

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.