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.

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 Methods:A Review and Recent Developments

The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to complex symbolic operations.Many variants have been developed to improve the HB method,among which the time domain HB-like methods are regarded as crucial improvements because of their fast computation and simple derivation.So far,there are two problems remaining to be addressed.i)A dozen of different versions of HB-like methods,in frequency domain or time domain or in hybrid,have been developed;unfortunately,misclassification pervades among them due to the unclear borderlines of different methods.ii)The time domain HB-like methods suffer from non-physical solutions,which have been shown to be caused by aliasing(mixture of the high-order into the low-order harmonics).Although a series of dealiasing techniques have been developed over the past two decades,the mechanism of aliasing and the final solution to dealiasing are still not well known to the academic community.This paper aims to provide a comprehensive review of the development of HB-like methods and enunciate their principal differences.In particular,the time domain methods are emphasized with the famous aliasing phenomenon clearly addressed.

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.

Nonlinear vibration analysis of a circular composite plate harvester via harmonic balance

Alumped parameter transversevibration model of a composite plate harvester is analyzed via harmonic balance approaches. The harvester is mainly composed of a piezoelectriccircular composite clamped by two steel rings and a proof mass on the plate.The lumped parameter model is a 1.5 degree-of-freedom strongly nonlinear system with a higher order polynomial stiffness. Aharmonic balance approach is developed to analyze the system, and the resulting algebraic equations are numerically solved by adopting an arc-length continuation technique. Anincremental harmonic balance approach is also developedfor the lumped parameter model. The two approaches yieldthe same results.The amplitude-frequency responses produced by the harmonic balance approach are validated by the numericalintegrations and the experimental data. The investigation reveals that there coexist hardening and softening characteristics in the amplitude-frequency response curves under sufficiently large excitations. The harvester with thecoexistenceof hardening and softening nonlinearitiescan outperform not only linear energy harvesters but also typical hardening nonlinear energy harvesters.

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.

Period-doubling bifurcation in two-stage power factor correction converters using the method of incremental harmonic balance and Floquet theory

In this paper, period-doubling bifurcation in a two-stage power factor correction converter is analyzed by using the method of incremental harmonic balance （IHB） and Floquet theory. A two-stage power factor correction converter typically employs a cascade configuration of a pre-regulator boost power factor correction converter with average current mode control to achieve a near unity power factor and a tightly regulated post-regulator DC-DC Buck converter with voltage feedback control to regulate the output voltage. Based on the assumption that the tightly regulated postregulator DC-DC Buck converter is represented as a constant power sink and some other assumptions, the simplified model of the two-stage power factor correction converter is derived and its approximate periodic solution is calculated by the method of IHB. And then, the stability of the system is investigated by using Floquet theory and the stable boundaries are presented on the selected parameter spaces. Finally, some experimental results are given to confirm the effectiveness of the theoretical analysis.

Nonlinear vibration analysis of a circular micro-plate in two-sided NEMS/MEMS capacitive system by using harmonic balance method

In this study, forced nonlinear vibration of a circular micro-plate under two-sided electrostatic, two-sided Casimir and external harmonic forces is investigated analytically. For this purpose, at first, von Karman plate theory including geometrical nonlinearity is used to obtain the deflection of the micro-plate. Galerkin decomposition method is then employed, and nonlinear ordinary differential equations (ODEs) of motion are determined. A harmonic balance method (HBM) is applied to equations and analytical relation for nonlineaT frequency response (F-R) curves are derived for two categories (including and neglecting Casimir force) separately. The analytical results for three cases:(1) semi-linear vibration;(2) weakly nonlinear vibration;(3) highly non linear vibration, are validated by comparing with the numerical solutio ns. After validation, the effects of the voltage and Casimir force on the natural frequency of two-sided capacitor system are investigated. It is shown that by assuming Casimir force in small gap distances, reduction of the natural frequency is considerable. The influences of the applied voltage, damping, micro-plate thickness and Casimir force on the frequency response curves have been presented too. The results of this study can be useful for modeling circular parallel-plates in nano /microelectromechanical transducers such as microphones and pressure sensors.

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.

Modeling, Steady-State Analysis of a SEPIC dc-dc Converter Based on Switching Function and Harmonic Balance Technique

The paper presents modeling approach of a Single Ended Primary Inductance Converter (SEPIC) system. The complete model derivation of the SEPIC converter system has been presented in different modes of operation. Steady state and small signal analysis was carried out on the converter dynamic equations using the method of Harmonic balance Technique. The steady state variables and their respective ripple quantities obtained were plotted against duty ratio D. The results obtained for a supply input voltage of 60 volts to the converter at a duty ratio of D = 0.8 , compares well with simulation results.

Analysing chaos in fractional-order systems with the harmonic balance method

In this paper, the fractional-order Genesio-Tesi system showing chaotic behaviours is introduced, and the corresponding one in an integer-order form is studied intensively. Based on the harmonic balance principle, which is widely used in the frequency analysis of nonlinear control systems, a theoretical approach is used to investigate the conditions of system parameters under which this fractional-order system can give rise to a chaotic attractor. Finally, the numerical simulation is used to verify the validity of the theoretical results.

Period-One Rotating Solutions of Horizontally Excited Pendulum Based on Iterative Harmonic Balance

In this study, the iterative harmonic balance method was used to develop analytical solutions of period-one rotations of a pendulum driven horizontally by harmonic excitations. The performance of the method was evaluated by two criteria, one based on the system energy error and the other based on the global residual error. As a comparison, analytical solutions based on the multi-scale method were also considered. Numerical solutions obtained from the Dormand-Prince method (ODE45 in MATLAB&#169) were used as the baseline for evaluation. It was found that under lower-level excitations, the multi-scale method performed better than the iterative method. At higher-level excitations, however, the performance of the iterative method was noticeably more accurate.

Improvement of Harmonic Balance Using Jacobian Elliptic Functions

We propose a method for finding approximate analytic solutions to autonomous single degree-of-freedom nonlinear oscillator equations. It consists of the harmonic balance with linearization in which Jacobian elliptic functions are used instead of circular trigonometric functions. We show that a simple change of independent variable followed by a careful choice of the form of anharmonic solution enable to obtain highly accurate approximate solutions. In particular our examples show that the proposed method is as easy to use as existing harmonic balance based methods and yet provides substantially greater accuracy.

AN IMPROVEMENT TO THE HARMONIC-BALANCE TECHNIQUE FOR ANALYZING NONLINEAR CIRCUITS

In this paper,an improved method,which is based on the harmonic-balance tech-nique,for speeding up the convergence iteration procedure is presented.As an example,it isapplied to the nonlinear analysis of microwave mixer,in which a simple time domain to fre-quency domain transformation technique for the mixer diode is firstly introduced.The use of thistechnique saves computing time dramatically.And the analyzed results are in good agreementwith the results published.The improved method can also be extended to the CAA of nonlinearcircuits with multi-devices.

Joint eigenvalue estimation by balanced simultaneous Schur decomposition

The problem of joint eigenvalue estimation for the non-defective commuting set of matrices A is addressed. A procedure revealing the joint eigenstructure by simultaneous diagonalization of. A with simultaneous Schur decomposition （SSD） and balance procedure alternately is proposed for performance considerations and also for overcoming the convergence difficulties of previous methods based only on simultaneous Schur form and unitary transformations, it is shown that the SSD procedure can be well incorporated with the balancing algorithm in a pingpong manner, i. e., each optimizes a cost function and at the same time serves as an acceleration procedure for the other. Under mild assumptions, the convergence of the two cost functions alternately optimized, i. e., the norm of A and the norm of the left-lower part of A is proved. Numerical experiments are conducted in a multi-dimensional harmonic retrieval application and suggest that the presented method converges considerably faster than the methods based on only unitary transformation for matrices which are not near to normality.

THE LARGE-SIGNAL DESIGN FOR FET FUNDAMENTAL AND HARMONIC OSCILLATORS

A Large-signal model for GaAs FET is derived based on its small-signal S parame-ters and DC characteristics. The harmonic balance algorithm is applied to analyze and optimizethe FET fundamental and harmonic oscillators, and the values of steady current are obtained.In the solving process, a simplified CAD approach is used to obtain the parameters of matchingnetwork when the output power is maximum. Finally, a fundamental oscillator and a harmonicoscillator of Q-band are fabricated. The measurements show that the theoretical analysis andexperimental results are in good agreement.

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.

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.