Exploiting Sparsity in SDP Relaxation for Harmonic Balance Method

In general, harmonic balance problems are extremely nonconvex and difficult to solve. A convex relaxation in the form of semidefinite programming has attracted a lot of attention recently, as it finds a global solution with high accuracy without the need for initial values. However, the computational cost of solving large-scale optimization poses a major challenge for the application in many real-world practical cases. This work proposes a heuristic optimization approach to find the Fourier coefficients of harmonic balance problems. The structural sparsity in the Harmonic Balance problem is exploited to improve numerical tractability and efficiency at the cost of adding smaller-sized semidefinite constraints in the problem formulation. After exploiting sparsity, the simulation results show that the size of the largest semidefinite constraint and the number of decision variables are greatly reduced. In addition, the computation speed shows an improvement rate of 3 or 7.5 times for larger instance problems and a reduction in memory occupation. Moreover, the proposed formulation can also solve nonlinear circuits with nonpolynomial nonlinearities with high accuracy.

[1]  S. Pearson.  Moments , 2020, Narrative inquiry in bioethics.

[2]  Malte Krack,et al.  Harmonic Balance for Nonlinear Vibration Problems , 2020, Mathematical Engineering.

[3]  C. Lim,et al.  Dynamic Response and Stability Analysis with Newton Harmonic Balance Method for Nonlinear Oscillating Dielectric Elastomer Balloons , 2018, International Journal of Structural Stability and Dynamics.

[4]  Nathan A. Wukie,et al.  Unsteady Turbomachinery Simulations Using Harmonic Balance on a Discontinuous Galerkin Discretization , 2018, Volume 2C: Turbomachinery.

[5]  Daniel K. Molzahn,et al.  Lasserre Hierarchy for Large Scale Polynomial Optimization in Real and Complex Variables , 2017, SIAM J. Optim..

[6]  L. Xie,et al.  Bifurcation tracking by Harmonic Balance Method for performance tuning of nonlinear dynamical systems , 2017 .

[7]  Federico Bizzarri,et al.  Harmonic Balance Based on Two-Step Galerkin Method , 2016, IEEE Transactions on Circuits and Systems I: Regular Papers.

[8]  Luc Masset,et al.  The harmonic balance method for bifurcation analysis of large-scale nonlinear mechanical systems , 2015, 1604.05621.

[9]  Takashi Tanaka,et al.  LQG Control With Minimum Directed Information: Semidefinite Programming Approach , 2015, IEEE Transactions on Automatic Control.

[10]  Martin S. Andersen,et al.  Chordal Graphs and Semidefinite Optimization , 2015, Found. Trends Optim..

[11]  Simao Marques,et al.  Prediction of Transonic Limit-Cycle Oscillations Using an Aeroelastic Harmonic Balance Method , 2015 .

[12]  J. Lasserre.  An Introduction to Polynomial and Semi-Algebraic Optimization , 2015 .

[13]  Earl H. Dowell,et al.  Harmonic balance methods applied to computational fluid dynamics problems , 2013 .

[14]  Andrew Y. T. Leung,et al.  Residue harmonic balance analysis for the damped Duffing resonator driven by a van der Pol oscillator , 2012 .

[15]  Fabrice Thouverez,et al.  On a new harmonic selection technique for harmonic balance method , 2012 .

[16]  R. Jabr.  Exploiting Sparsity in SDP Relaxations of the OPF Problem , 2012, IEEE Transactions on Power Systems.

[17]  Chenjie Gu,et al.  QLMOR: A Projection-Based Nonlinear Model Order Reduction Approach Using Quadratic-Linear Representation of Nonlinear Systems , 2011, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[18]  Martin S. Andersen,et al.  Linear matrix inequalities with chordal sparsity patterns and applications to robust quadratic optimization , 2010, 2010 IEEE International Symposium on Computer-Aided Control System Design.

[19]  J. Lasserre.  Moments, Positive Polynomials And Their Applications , 2009 .

[20]  Masakazu Kojima,et al.  Exploiting Sparsity in SDP Relaxation for Sensor Network Localization , 2009, SIAM J. Optim..

[21]  Masakazu Kojima,et al.  Sparsity in sums of squares of polynomials , 2005, Math. Program..

[22]  R. W. Menzies,et al.  Stability Domain Calculations of Period-1 Ferroresonance in a Nonlinear Resonant Circuit , 2002, IEEE Power Engineering Review.

[23]  Jacob K. White,et al.  A trajectory piecewise-linear approach to model order reduction and fast simulation of nonlinear circuits and micromachined devices , 2001, IEEE/ACM International Conference on Computer Aided Design. ICCAD 2001. IEEE/ACM Digest of Technical Papers (Cat. No.01CH37281).

[24]  Kazuo Murota,et al.  Exploiting Sparsity in Semidefinite Programming via Matrix Completion I: General Framework , 2000, SIAM J. Optim..

[25]  S. H. A. Chen,et al.  Application of the incremental harmonic balance method to cubic non-linearity systems , 1990 .

[26]  Charles R. Johnson,et al.  Positive definite completions of partial Hermitian matrices , 1984 .

[27]  Cheng H. Yang,et al.  A Semidefinite Programming Approach for Harmonic Balance Method , 2019, IEEE Access.

[28]  Robin Deits,et al.  Computing Large Convex Regions of Obstacle-Free Space Through Semidefinite Programming , 2014, WAFR.

[29]  F. H. Ling,et al.  An Alternating Frequency/Time Domain Method for Calculating the Steady-State Response of Nonlinear Dynamic Systems , 2007 .

[30]  M. Kojima,et al.  Sums of Squares and Semidefinite Program Relaxations for Polynomial Optimization Problems with Structured Sparsity , 2006, SIAM J. Optim..

[31]  J. Lofberg,et al.  YALMIP : a toolbox for modeling and optimization in MATLAB , 2004, 2004 IEEE International Conference on Robotics and Automation (IEEE Cat. No.04CH37508).

[32]  Daniel Liberzon,et al.  Common Lyapunov functions for families of commuting nonlinear systems , 2005, Syst. Control. Lett..