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.

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.

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.

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.

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.

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.

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.

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.

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.

Adjacent mode resonance of a hydraulic pipe system consisting of parallel pipes coupled at middle points

The coupling vibration of a hydraulic pipe system consisting of two pipes is studied.The pipes are installed in parallel and fixed at their ends,and are restrained by clips to one bracket at their middle points.The pipe subjected to the basement excitation at the left end is named as the active pipe,while the pipe without excitation is called the passive pipe.The clips between the two pipes are the bridge for the vibration energy.The adjacent natural frequencies will enhance the vibration coupling.The governing equation of the coupled system is deduced by the generalized Hamilton principle,and is discretized to the modal space.The modal correction is used during the discretization.The investigation on the natural characters indicates that the adjacent natural frequencies can be adjusted by the stiffness of the two clips and bracket.The harmonic balance method(HBM)is used to study the responses in the adjacent natural frequency region.The results show that the vibration energy transmits from the active pipe to the passive pipe swimmingly via the clips together with a flexible bracket,while the locations of them are not node points.The adjacent natural frequencies may arouse wide resonance curves with two peaks for both pipes.The stiffness of the clip and bracket can release the vibration coupling.It is suggested that the stiffness of the clip on the passive pipe should be weak and the bracket should be strong enough.In this way,the vibration energy is reflected by the almost rigid bracket,and is hard to transfer to the passive pipe via a soft clip.The best choice is to set the clips at the pipe node points.The current work gives some suggestions for weakening the coupled vibration during the dynamic design of a coupled hydraulic pipe system.

A novel way for vibration control of FGM fluid-conveying pipes via NiTiNOL-steel wire rope

In this study,a coupling model of fluid-conveying pipes made of functionally graded materials(FGMs)with NiTiNOL-steel(NiTi-ST)for vibration absorption is investigated.The vibration responses of the FGM fluid-conveying pipe with NiTi-ST are studied by the Galerkin truncation method(GTM)and harmonic balance method(HBM).The harmonic balance solutions and the numerical results are consistent.Also,the linearized stability of the structure is determined.The effects of the structure parameters on the absorption performance are also studied.The results show that the NiTi-ST is an effective means of vibration absorption.Furthermore,in studying the effect of the NiTi-ST,a closed detached response(CDR)is first observed.It is noteworthy that the CDR may dramatically change the vibration amplitude and that the parameters of the NiTi-ST may determine the emergence or disappearance of the CDR.This vibration absorption device can be extended to offer more general vibration control in engineering applications.

Nonlinear dynamic analysis of cantilevered piezoelectric energy harvesters under simultaneous parametric and external excitations

The nonlinear dynamics of cantilevered piezoelectric beams is investigated under simultaneous parametric and external excitations. The beam is composed of a substrate and two piezoelectric layers and assumed as an Euler-Bernoulli model with inextensible deformation. A nonlinear distributed parameter model of cantilevered piezoelectric energy harvesters is proposed using the generalized Hamilton＇s principle. The proposed model includes geometric and inertia nonlinearity, but neglects the material nonlinearity. Using the Galerkin decomposition method and harmonic balance method, analytical expressions of the frequency-response curves are presented when the first bending mode of the beam plays a dominant role. Using these expressions, we investigate the effects of the damping, load resistance, electromechanical coupling, and excitation amplitude on the frequency-response curves. We also study the difference between the nonlinear lumped-parameter and distributed- parameter model for predicting the performance of the energy harvesting system. Only in the case of parametric excitation, we demonstrate that the energy harvesting system has an initiation excitation threshold below which no energy can be harvested. We also illustrate that the damping and load resistance affect the initiation excitation threshold.

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.

Nonlinear vibration of functionally graded graphene platelet-reinforced composite truncated conical shell using first-order shear deformation theory

In this study,the first-order shear deformation theory(FSDT)is used to establish a nonlinear dynamic model for a conical shell truncated by a functionally graded graphene platelet-reinforced composite(FG-GPLRC).The vibration analyses of the FG-GPLRC truncated conical shell are presented.Considering the graphene platelets(GPLs)of the FG-GPLRC truncated conical shell with three different distribution patterns,the modified Halpin-Tsai model is used to calculate the effective Young’s modulus.Hamilton’s principle,the FSDT,and the von-Karman type nonlinear geometric relationships are used to derive a system of partial differential governing equations of the FG-GPLRC truncated conical shell.The Galerkin method is used to obtain the ordinary differential equations of the truncated conical shell.Then,the analytical nonlinear frequencies of the FG-GPLRC truncated conical shell are solved by the harmonic balance method.The effects of the weight fraction and distribution pattern of the GPLs,the ratio of the length to the radius as well as the ratio of the radius to the thickness of the FG-GPLRC truncated conical shell on the nonlinear natural frequency characteristics are discussed.This study culminates in the discovery of the periodic motion and chaotic motion of the FG-GPLRC truncated conical shell.

NES cell

A broadband adaptive vibration control strategy with high reliability and flexible versatility is proposed.The high vibration damping performance of nonlinear energy sink(NES)has attracted attention.However,targeted energy transfer may cause severe vibration of NES.Besides,it is difficult to realize pure nonlinear stiffness without the linear part.As a result,the reliability of NES is not high.The low reliability of NES has hindered its application in engineering.In addition,the performance of NES depends on its mass ratio of the primary system,and NES lacks versatility for different vibration systems.Therefore,this paper proposes the concept of NES cell.The advantages of the adaptive vibration control of NES are applied to cellular NES.By applying a large number of NES cells in parallel,the reliability of NES and its versatility to complex vibration structures are improved.An elastic beam is used as the primary vibration structure,and a limited NES is used as the cell.The relationship between the vibration suppression effect of NES cells and the number of NES cell is studied.In addition,the effect of the distribution of NES cells on the multi-mode resonance suppression of the beam is also studied.In summary,the mode of the primary structure can be efficiently controlled by a large number of lightweight NES cell.Therefore,the lightweight NES cell is flexible for vibration control of complex structures.In addition,it improves the reliability of NES applications.Therefore,the distributed application of NES cells proposed in this paper is a valuable vibration suppression strategy.