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 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.

Incremental harmonic balance method for periodic forced oscillation of a dielectric elastomer balloon

Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.

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.

Spectral Method for Three-Dimensional Nonlinear Klein-Gordon Equation by Using Generalized Laguerre and Spherical Harmonic Functions

In this paper,a generalized Laguerre-spherical harmonic spectral method is proposed for the Cauchy problem of three-dimensional nonlinear Klein-Gordon equation. The goal is to make the numerical solutions to preserve the same conservation as that for the exact solution.The stability and convergence of the proposed scheme are proved.Numerical results demonstrate the efficiency of this approach.We also establish some basic results on the generalized Laguerre-spherical harmonic orthogonal approximation,which play an important role in spectral methods for various problems defined on the whole space and unbounded domains with spherical geometry.

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.

Fully discrete Jacobi-spherical harmonic spectral method for Navier-Stokes equations

A fully discrete Jacobi-spherical harmonic spectral method is provided for the Navier-Stokes equations in a ball. Its stability and convergence are proved. Numerical results show efficiency of this approach. The proposed method is also applicable to other problems in spherical geometry.

APPLICATION OF HARMONIC ANALYSIS METHOD TO RESEARCH ON ROTOR AIRLOADS

According to the rotor vortex theory,the rotor circulation and the rotor induced velocity are developed into Fourier series.The circulation distribution along blade spanwise is expressed in terms of segment-by-segment linear functions.In consequence the induced velocity equations and the circulation equations are derived.The engineering application of the rotor vortex theory is provided.Then the induced velocity and its harmonic components are obtained to provide a quantitative basis for the vortex model.For calculating each harmonic component of the induced velocity a simplified method is put forward which considers the effects of each order circulation with neglecting those of higher order.The method saves the computer time and is of significant benefit.

Harmonics analysis of the ITER poloidal field converter based on a piecewise method

Poloidal field（PF） converters provide controlled DC voltage and current to PF coils. The many harmonics generated by the PF converter flow into the power grid and seriously affect power systems and electric equipment. Due to the complexity of the system, the traditional integral operation in Fourier analysis is complicated and inaccurate. This paper presents a piecewise method to calculate the harmonics of the ITER PF converter. The relationship between the grid input current and the DC output current of the ITER PF converter is deduced. The grid current is decomposed into the sum of some simple functions. By calculating simple function harmonics based on the piecewise method, the harmonics of the PF converter under different operation modes are obtained.In order to examine the validity of the method, a simulation model is established based on Matlab/Simulink and a relevant experiment is implemented in the ITER PF integration test platform.Comparative results are given. The calculated results are found to be consistent with simulation and experiment. The piecewise method is proved correct and valid for calculating the system harmonics.

A harmonic-wave bio-thermal method for continuous monitoring skin thermal conductivity and capillary perfusion rate

The revelation of thermal energy exchange mechanism of human body is challenging yet worthwhile,because it can clearly explain the changes in human symptoms and health status.Understanding,the heat transfer of the skin is significant because the skin is the foremost organ for the energy exchange between the human body and the environment.In order to diagnose the physiological conditions of human skin without causing any damage,it is necessary to use a non-invasive measurement technique by means of a conformal flexible sensor.The harmonic method can minimize the thermal-induced injury to the skin due to its low heat generating properties.A novel type of computational theory assessing skin thermal conductivity,blood perfusion rate of capillaries in the dermis,and superficial subcutaneous tissues was formed by combining the multi-medium thermal diffusion model and the bio-thermal model(Pennes equation).The skins of the hand back of six healthy subjects were measured.It was found that the results revealed no consistent changes in thermal conductivity were observed across genders and ages.The measured blood perfusion rates were within the range of human capillary flow.It was found that female subjects had a higher perfusion rate range(0.0058-0.0061 s^(-1))than male subjects(0.0032-0.0049 s^(-1)),which is consistent with invasive medical studies about the gender difference in blood flow rates and stimulated effects in relaxation situations.

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.

Odd–even harmonic emission from asymmetric molecules:Identifying the mechanism

We study odd–even high-order harmonic generation（HHG） from oriented asymmetric molecules He H2＋numerically and analytically. The variational method is used to improve the analytical description of the ground-state wave function for the asymmetric system, with which the ground-state-continuum-state transition dipole is evaluated. The comparison between the odd–even HHG spectra and the improved dipoles allows us to identify and clarify the complex generation mechanism of odd–even harmonics from asymmetric molecules, providing deep insights into the relation between the odd–even HHG and the asymmetric molecular orbital.

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.

VISCO—ELASTIC SYSTEMS UNDER BOTH DETERMINISTIC HARMONIC AND RANDOM EXCITATION

The response of visco_elastic system to combined deterministic harmonic and random excitation was investigated. The method of harmonic balance and the method of stochastic averaging were used to determine the response of the system. The theoretical analysis was verified by numerical results. Theoretical analyses and numerical simulations show that when the intensity of the random excitation increase, the nontrivial steady state solution may change from a limit cycle to a diffused limit cycle. Under some conditions the system may have two steady state solutions and jumps may exist.

1/3 SUBHARMONIC SOLUTION OF ELLIPTICAL SANDWICH PLATES

The problem of nonlinear forced oscillations for elliptical sandwich plates is dealt with. Based on the governing equations expressed in terms of five displacement components, the nonlinear dynamic equation of an elliptical sandwich plate under a harmonic force is derived. A superpositive-iterative harmonic balance (SIHB) method is presented for the steady-state analysis of strongly nonlinear oscillators. In a periodic oscillation, the periodic solutions can be expressed in the form of basic harmonics and bifurcate harmonics. Thus, an oscillation system which is described as a second order ordinary differential equation, can be expressed as fundamental differential equation with fundamental harmonics and incremental differential equation with derived harmonics. The 1/3 subharmonic solution of an elliptical sandwich plate is investigated by using the methods of SIHB. The SIHB method is compared with the numerical integration method. Finally, asymptotical stability of the 1/3 subharmonic oscillations is inspected.

Professor Ruixia Pei's Clinical Experience in Treating Globus Hystericus with"Harmonizing Method"

To summarize professor Pei9s experience in treating globus hystericus.Methods:Learn from your teacher.Results:Professor Pei had been suffering from the disease for more than 30 years and had her unique views on the treatment of globus hystericus,which are mostly for emotional dysfunction,liver Qi is not comfortable,functioning of Qi is not adjusted,Yin and Yang imbalance.The basic principle of treatment is to grasp the core pathogenesis,take harmony as the method and balance as the duration,harmonize Qi,regulate Yin and Yang,use drugs to disperse the liver and rectify Qi,as well as auxiliary treatment with products to promote blood circulation,remove blood stasis,dryness and dampness,and clear heat.

Clinical Experience of"'Harmonizing MethodM based on Usage of Minor Bupleurum Decoction

Minor Bupleurum Decoction is a common TCM compound and a classic prescription that can best reflect the"harmonizing method',of TCM.But there are many differences on the usage and dosage of Minor Bupleurum Decoction.Combined with the clinical treatment experience,this paper summarizes the experience of Minor Bupleurum Decoction in the clinical treatment of diseases.It is hoped that relevant suggestions can be provided for clinicians to improve the clinical therapeutic effect of Minor Bupleurum Decoction.