A novel adaptive harmonic balance method with an asymptotic harmonic selection

The harmonic balance method(HBM)is one of the most widely used methods in solving nonlinear vibration problems,and its accuracy and computational efficiency largely depend on the number of the harmonics selected.The adaptive harmonic balance(AHB)method is an improved HBM method.This paper presents a modified AHB method with the asymptotic harmonic selection(AHS)procedure.This new harmonic selection procedure selects harmonics from the frequency spectra of nonlinear terms instead of estimating the contribution of each harmonic to the whole nonlinear response,by which the additional calculation is avoided.A modified continuation method is proposed to deal with the variable size of nonlinear algebraic equations at different values of path parameters,and then all solution branches of the amplitude-frequency response are obtained.Numerical experiments are carried out to verify the performance of the AHB-AHS method.Five typical nonlinear dynamic equations with different types of nonlinearities and excitations are chosen as the illustrative examples.Compared with the classical HBM and Runge-Kutta methods,the proposed AHB-AHS method is of higher accuracy and better convergence.The AHB-AHS method proposed in this paper has the potential to investigate the nonlinear vibrations of complex high-dimensional nonlinear systems.

Harmonic balance method with alternating frequency/time domain technique for nonlinear dynamical system with fractional exponential

Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time （HB-AFT） do- main technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.

Improving convergence of incremental harmonic balance method using homotopy analysis method

We have deduced incremental harmonic balance an iteration scheme in the （IHB） method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.

Application of the Material Balance Method in Paleoelevation Recovery: A Case Study of the Longmen Mountains Foreland Basin on the Eastern Margin of the Tibetan Plateau

We applied the material balance principle of the denudation volume and sedimentary flux to study the denudation-accumulation system between the Longmen Mountains （Mts.） and the foreland basin. The amount of sediment in each sedimentation stage of the basin was estimated to obtain the denudation volume, erosion thickness and deposit thickness since the Late Triassic Epoch, to enable us to recover the paleoelevation of the provenance and the sedimentary area. The results show the following： （1） Since the Late Triassic Epoch, the elevation of the surface of the Longmen Mts. has uplifted from 0 m to 2751 m, and the crust of the Longmen Mts. has uplifted by 9.8 km. Approximately 72% of the materials introduced have been denuded from the mountains. （2） It is difficult to recover the paleoelevation of each stage of the Longmen Mts. foreland basin quantitatively by the present-day techniques and data. （3） The formation of the Longmen Mts. foreland basin consisted of three stages of thrust belt tectonic load and three stages of thrust belt erosional unload. During tectonic loading stages （Late Triassic Epoch, Late Jurassic-Early Cretaceous, Late Cretaceous-Miocene）, the average elevation of Longmen Mts. was lower （approximately 700-1700 m）. During erosional unloading stages （Early and Middle Jurassic, Middle Cretaceous and Jiaguan, Late Cenozoic）, the average elevation of Longmen Mts. was high at approximately 2000-2800m.

A Modified Homogeneous Balance Method and Its Applications

A modified homogeneous balance method is proposed by improving some key steps in the homogeneousbalance method.Bilinear equations of some nonlinear evolution equations are derived by using the modified homogeneousbalance method.Generalized Boussinesq equation,KP equation,and mKdV equation are chosen as examples to illustrateour method.This approach is also applicable to a large variety of nonlinear evolution equations.

INCREMENTAL HARMONIC BALANCE METHOD FOR AIRFOIL FLUTTER WITH MULTIPLE STRONG NONLINEARITIES

The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.

Primary resonance of fractional-order Duffing–van der Pol oscillator by harmonic balance method

The dynamical properties of fractional-order Duffing–van der Pol oscillator are studied, and the amplitude–frequency response equation of primary resonance is obtained by the harmonic balance method. The stability condition for steady-state solution is obtained based on Lyapunov theory. The comparison of the approximate analytical results with the numerical results is fulfilled, and the approximations obtained are in good agreement with the numerical solutions. The bifurcations of primary resonance for system parameters are analyzed. The results show that the harmonic balance method is effective and convenient for solving this problem, and it provides a reference for the dynamical analysis of similar nonlinear systems.

Homogenous Balance Method and Exact Analytical Solutions for Whitham-Broer-Kaup Equations in Shallow Water

Based on the homogenous balance method and with the help of mathematica, the Backlund transformation and the transfer heat equation are derived. Analyzing the heat-transfer equation, the multiple soliton solutions and other exact analytical solution for Whitham-Broer-Kaup equations(WBK) are derived. These solutions contain Fan's, Xie's and Yan's results and other new types of analytical solutions, such as rational function solutions and periodic solutions. The method can also be applied to solve more nonlinear differential equations.

Exact Solutions for the Generalized KdV Equation： Modified Homogeneous Balance Method

In this paper, we present a solution methodology to obtain exact solutions of some nonlinear evolution equation by modifying the homogeneous balance method. Based on the modified homogeneous balance method, several kinds of exact（new） solutions of the generalized KdV equation are obtained.

A semi-analytical multi-harmonic balance method on full-3D contact model for dynamic analysis of dry friction systems

Dry friction damping structures are widely-used in aero-engines to mitigate vibration.The nonlinear nature of friction and the two-dimensional in-plane motion on the contact interface bring challenges to accurately and efficiently predict the forced response of frictionally damped structures.The state-of-the-art Multi-Harmonic Balance Method(MHBM)on quasi-3D contact model in engineering cannot precisely capture the kinematics on the friction interface although the efficiency is high.The full-3D contact model can describe the constitutive relationship of the interface in a more accurate manner;however,the efficiency and convergence are not guaranteed for large-scale models.In this paper,a semi-analytical MHBM on full-3D contact model is proposed.The original Trajectory Tracking Method(TTM)for evaluating the contact force is reformulated to make the calculation more concise and the derivation of the Analytical Jacobian Matrix(AJM)feasible.Based on the chain rule of derivation,the AJM which is the core to upgrade the performance is deduced.Through a shrouded blade finite element model,the accuracy and efficiency of the proposed method are compared with both the MHBM on full-3D contact model with numerical Jacobian matrix and the MHBM on quasi-3D contact model with AJM.The results show that the AJM improves significantly the efficiency of the MHBM on full-3D contact model.The time cost of the proposed method is in the same order of magnitude as that of the MHBM on quasi-3D contact model.We also confirm that the full-3D contact model is necessary for the dynamic analyses of shrouded blades.If one uses the quasi-3D model,the estimation relative error of damping can even reach 31.8%in some cases.In addition,the AJM also brings benefits for stability analysis.It is highly recommended that engineers use the MHBM on full-3D contact model for the dynamic analysis and design of shrouded blades.

Harmonic balance simulation of the influence of component uniformity and reliability on the performance of a Josephson traveling wave parametric amplifier

A Josephson traveling wave parametric amplifier(JTWPA),which is a quantum-limited amplifier with high gain and large bandwidth,is the core device of large-scale measurement and control systems for quantum computing.A typical JTWPA consists of thousands of Josephson junctions connected in series to form a transmission line and hundreds of shunt LC resonators periodically loaded along the line for phase matching.Because the variation of these capacitors and inductors can be detrimental to their high-frequency characteristics,the fabrication of a JTWPA typically necessitates precise processing equipment.To guide the fabrication process and further improve the design for manufacturability,it is necessary to understand how each electronic component affects the amplifier.In this paper,we use the harmonic balance method to conduct a comprehensive study on the impact of nonuniformity and fabrication yield of the electronic components on the performance of a JTWPA.The results provide insightful and scientific guidance for device design and fabrication processes.

Determination of initial cable force of cantilever casting concrete arch bridge using stress balance and influence matrix methods

Cantilever casting concrete arch bridge using form traveller has a broad application prospect.However,it is difficult to obtain reasonable initial cable force in construction stage.In this study,stress balance and influence matrix methods were developed to determine the initial cable force of cantilever casting concrete arch bridge.The stress balance equation and influence matrix of arch rib critical section were established,and the buckle cable force range was determined by the allowable stress of arch rib critical section.Then a group of buckle cable forces were selected and substituted into the stress balance equation,and the reasonable initial buckle cable force was determined through iteration.Based on the principle of force balance,the initial anchor cable force was determined.In an engineering application example,it is shown that the stress balance and influence matrix methods for the determination of initial cable force are feasible and reliable.The initial cable forces of arch rib segments only need to be adjusted once in the corresponding construction process,which improves the working efficiency and reduces the construction risk.It is found that the methods have great advantages for determining initial cable force in cantilever casting construction process of concrete arch bridge.

Periodic response analysis of a misaligned rotor system by harmonic balance method with alternating frequency/time domain technique

A dynamic model is established for an offset-disc rotor system with a mechanical gear coupling, which takes into consideration the nonlinear restoring force of rotor support and the effect of coupling misalignment. Periodic solutions are obtained through harmonic balance method with alternating frequency/time domain(HB-AFT) technique, and then compared with the results of numerical simulation. Good agreement confirms the feasibility of HB-AFT scheme. Moreover, the Floquet theory is adopted to analyze motion stability of the system when rotor runs at different speed intervals. A simple strategy to determine the monodromy matrix is introduced and two ways towards unstability are found for periodic solutions: the period doubling bifurcation and the secondary Hopf bifurcation. The results obtained will contribute to the global response analysis and dynamic optimal design of rotor systems.

Revisiting the space-time gradient method:A time-clocking perspective, high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric GaussSeidel/Gauss-Seidel(LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-countreduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.

Solution of nonlinear cubic-quintic Duffing oscillators using He’s Energy Balance Method

In this study, He＇s Energy Balance Method （EBM） was applied to solve strong nonlinear Duffing oscillators with cubic-quintic nonlinear restoring force. The complete EBM solution procedure of the cubic-quintic Duffing oscillator equation is presented. For illustration of effectiveness and convenience of the EBM, different cases of cubic-quintic Duffing oscillator with different parameters of α,β and y were compared with the exact solution. We found that the solutions were valid for small as well as large amplitudes of oscillation. The results show that the EBM is very convenient and precise, so it can be widely appli- cable in engineering and other sciences.

Almost Periodic Solutions by the Harmonic Balance Method

We consider non-autonomous ordinary differential equations in two cases.One is the one-dimensional case that admits a condition of hyperbolicity,and the other one is the higher-dimensional case that admits an exponential dichotomy.For differential equations of this kind,we give a rigorous treatment of the Harmonic Balance Method to look for almost periodic solutions.

Numerical Analysis of Balanced Methods for the Impulsive Stochastic Differential Equations

Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong convergence can preserve the mean square stability with the sufficiently small stepsize. Weak variants and their mean square stability are also considered. Several numerical experiments are given for illustration and show that the fully implicit methods are superior to those of the explicit methods in terms of mean-square stabilities for relatively large stepsizes especially.

General Modified Split-Step Balanced Methods for Stiff Stochastic Differential Equations

A class of general modified split-step balanced methods proposed in the paper can be applied to solve stiff stochastic differential systems with m-dimensional multiplicative noise. Compared to some other already reported split-step balanced methods, the drift increment function of the methods can be taken from any chosen ane-step ordinary differential equations （ODEs） solver. The schemes is proved to be strong convergent with order one. For the mean-square stability analysis, the investigation is confined to two cases. Some numerical experiments are reported to testify the performance and the effectiveness of the methods.

A vertical track nonlinear energy sink

Eliminating the effects of gravity and designing nonlinear energy sinks(NESs)that suppress vibration in the vertical direction is a challenging task with numerous damping requirements.In this paper,the dynamic design of a vertical track nonlinear energy sink(VTNES)with zero linear stiffness in the vertical direction is proposed and realized for the first time.The motion differential equations of the VTNES coupled with a linear oscillator(LO)are established.With the strong nonlinearity considered of the VTNES,the steady-state response of the system is analyzed with the harmonic balance method(HBM),and the accuracy of the HBM is verified numerically.On this basis,the VTNES prototype is manufactured,and its nonlinear stiffness is identified.The damping effect and dynamic characteristics of the VTNES are studied theoretically and experimentally.The results show that the VTNES has better damping effects when strong modulation responses(SMRs)occur.Moreover,even for small-amplitude vibration,the VTNES also has a good vibration suppression effect.To sum up,in order to suppress the vertical vibration,an NES is designed and developed,which can suppress the vertical vibration within certain ranges of the resonance frequency and the vibration intensity.

时域谐波平衡求解器中非物理解出现根源探析

The time domain harmonic balance method is an attractive reduced order method of analyzing unsteady flow for turbomachines. However, the method can admit non-physical solutions. Non-physical solutions were encountered from a three-blade-row compressor configuration in a time domain harmonic balance analysis. This paper aims to investigate the root cause of the non-physical solutions. The investigation involves several strategies, which include increasing the number of harmonics, increasing the number of time instants, including scattered modes,including the rotor-rotor interaction, and the use of a new method-the approximate time domain nonlinear harmonic method. Numerical analyses pertinent to each strategy are presented to reveal the root cause of the non-physical solution. It is found that the nonlinear interaction of unsteady flow components with different fundamental frequencies is the cause of the non-physical solution. The non-physical solution can be eliminated by incorporating extra scattered modes or using the approximate time domain nonlinear harmonic method.