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

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.

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.

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.

The Reduced Space Method for Calculating the Periodic Solution of Nonlinear Systems

A hybrid method combined the reduced Sequential Quadratic Programming(SQP)method with the harmonic balance method has been developed to analyze the characteristics of mode localization and internal resonance of nonlinear bladed disks.With the aid of harmonic balance method,the nonlinear equality constraints for the constrained optimization problem are constructed.The reduced SQP method is then utilized to deal with the original constrained optimization problem.Applying the null space decomposition technique to the harmonic balance algebraic equations results in the vanishing of the nonlinear equality constraints and a simple optimization problem involving only upper and lower bound constraints on the optimization variables is formed and solved.Finally,numerical results are given for several test examples to validity the proposed method.The efficiency of the solution method to trace the family of energy dependent nonlinear modes is illustrated.The localization nonlinear normal modes of bladed disks related to various types of internal resonances are explored.

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.

混联Ⅱ型惯容非线性能量阱的动力学特性研究

分别使用非线性恢复力、非线性阻尼替代惯容减振系统中的线性恢复力、线性阻尼,并考虑摩擦力的影响,提出了混联Ⅱ型惯容非线性能量阱.建立了主系统的动力学方程,利用谐波平衡法求解系统在简谐激励下的幅频响应曲线.采用弧长算法和数值法相结合的方法研究了系统的惯质比、非线性阻尼、非线性刚度和摩擦力单个参数对其减振性能的影响.发现非线性刚度和非线性阻尼数值的增大会使峰值先减小后增大,不同的是,前者幅频响应曲线逐渐向右上方向弯曲,后者产生峰值的位置向低频段转移.分析了惯质比、非线性阻尼、非线性刚度3种参数两两组合下对系统减振效果的影响.研究表明,在激励幅值为0.005 m时,惯质比和阻尼同时变化减振效果最好:当ε=0.1时,系统主结构位移峰值的最小值约为0.01 m;而在参数ε=0.001时,整体取值范围内其最大值约为0.061 m;当惯质比取得最佳值0.1时,非线性阻尼和非线性刚度κ_(21)的取值范围变大.在摩擦力的作用下,系统的最大幅值都有不同程度的增加.上述研究可为振动系统减振的研究提供参考.

面向光栅刻划机的负刚度结构隔振性能优化

大面积光栅刻划机实验室存在的复杂振动因素,会对它的运行及刻画精度造成影响,针对光栅刻划机设计一种正负刚度并联的被动隔振结构,改善仅有空气弹簧时对低频隔振无效的现象,拓宽隔振频带。为模拟隔振器隔振性能,建立系统的动力学模型和动力学方程,使用仿真软件MATLAB进行仿真,为磁致负刚度弹簧进行尺寸结构设计,对比并联磁致负刚度弹簧前后的隔振效果。仿真结果表明采用空气弹簧和磁致负刚度结构并联后能够有效的隔离低频振动,并且磁吸致负刚度和磁斥致负刚度进行组合能够有效的改善单一磁吸致负刚度或磁斥致负刚度结构的非线性。

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.

直齿锥齿轮系统动态特性参数解域界结构

为研究含间隙直齿锥齿轮系统的幅值跳跃、齿面冲击特性,用谐波平衡法建立系统动力学方程。在频率与间隙、时变刚度、静态综合误差和负载构建的参数平面内用Broyden拟牛顿法和伪弧长延续算法获得系统参数解域界结构,探索系统幅值跳跃、多值解和啮合冲击特性对参数的敏感性。结果表明在系统啮频和轴频主共振区均出现幅值跳跃和多值解现象,在啮频处动态特性更为复杂丰富。小间隙下存在严重的齿面冲击现象,当间隙b>0.98时系统的跳跃与齿面冲击均趋于稳定,系统动态特性对时变啮合刚度不敏感,在高速轻载、较大齿面误差激励下系统的非线性跳跃和齿面冲击均加剧。参数解域界结构为锥齿轮结构设计提供数据支持。

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.

重构谐波平衡法及其求解复杂非线性问题应用

谐波平衡法是求解非线性动力学系统周期解的最常用方法,但对非线性项进行高阶近似需要庞杂的公式推导,限制了该方法的超高精度解算.通过对频域非线性量的时域等价重构,提出了重构谐波平衡法(RHB法),解决了经典谐波平衡法超高阶次计算难题.然而,上述两种方法均要求动力学系统为多项式型非线性,且无法直接用来求解非线性系统的拟周期解.针对上述问题,文章提出一种将RHB法和复杂非线性系统等价重铸法相结合的计算方法,首先将一般非线性问题无损重铸为多项式型非线性系统,然后用RHB法进行高精度求解;针对拟周期响应求解问题,提出基于“补频”思想的RHB方法,通过基频的优化筛选,实现拟周期响应的快速精准捕捉.选取非线性单摆、相对论谐振子和非线性耦合非对称摆等典型系统进行仿真计算,仿真结果表明,所提出的RHB-重铸法在解非多项式型非线性系统的稳态响应时精度保持为10^(-12)量级,达计算机精度,远超现有方法水平.补频RHB法则实现了对拟周期问题的高效解算,拓宽了方法对真实物理响应的求解范围.