REVIEW 3 major objections 6 minor 137 references
pyHB: an open-source automatic-differentiation-enhanced semi-analytical solver for nonlinear dynamics
T0 review · 3 major / 6 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read A solver that needs only the system equations, not hand-derived Jacobians, maps full periodic-response landscapes—stable and unstable—of million-scale vibration problems.
desk verdict A real and well-tested HB tool whose headline numbers come from localized nonlinearities; the 'general/one-stop' wrapper overclaims that scope, but the core is worth serious engagement. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the localized-nonlinearity decomposition f_nl = Ξ f_nl^(r)(Θx'', Θx', Θx, τ, ω) with q,r << n. It is what lets automatic differentiation operate on a low-dimensional computational graph; it also makes the nonlinear Jacobian a sparse scatter of a small dense tensor. Around it, three pieces do the heavy lifting: precomputed Fourier-basis tensors for the nonlinear Jacobian update; a Kronecker-structured linear Jacobian with sparse LU factorization; and a blocked solve of the arc-length-augmented linear system that reuses one sparse factorization per continuation step. Periodic-orbit stability is computed by propagating a monodromy matrix via an implicit trapezoidal in
What would settle it
Run pyHB on a system with n=2000 where the nonlinear force couples every DOF to every other (e.g., a fully populated cubic stiffness matrix) with m=50 harmonics. If per-iteration time or GPU memory scales with n rather than staying near 0.44 s and 250 MB, the central scalability claim fails. Reproducing the reported 2000-DOF beam benchmark and measuring total wall time for a full branch sweep would also confirm or refute the several-hundred-fold speedup.
Extended reading notes
Core claim
The central claim is that the entire harmonic balance workflow can be driven by user-defined equations alone. The key move is rewriting the nonlinear force as a selection-matrix product that restricts it to a few degrees of freedom, so automatic differentiation builds the required tangent matrices for this reduced force, and the Jacobian for the full system is assembled from precomputed linear blocks plus a sparse nonlinear contribution obtained by FFT-based harmonic projection. The paper reports that this makes a 2000-DOF discretized beam with 50 harmonics (about 202,000 unknowns) solvable at roughly 0.44 s per continuation point, several hundred times faster than stepping through ten perio
Load-bearing premise
The paper assumes nonlinear forces act on only a few degrees of freedom (local contact, joint, or support), so the automatic-differentiation graph stays small; if a nonlinearity couples many degrees of freedom at once, the claimed speed and flat memory would not hold.
Editorial extensions
If this is right
- Users can change a nonlinear force model by editing one function, without re-deriving tangent stiffness or damping matrices, so design iterations in isolators, harvesters, and rotors become fast.
- Full stable and unstable branches become routinely available for high-dimensional finite-element models, allowing prediction of jump phenomena, hysteresis, and hidden resonances that time integration misses.
- The memory footprint of automatic-differentiation-enhanced harmonic balance no longer grows with the number of degrees of freedom or harmonics in the observed cases, making single-GPU analysis of 2000-DOF systems with 50 harmonics feasible.
- A standard benchmark platform emerges: the same modular pipeline reproduces low-DOF strongly nonlinear isolator responses and high-DOF aeroengine and beam responses, facilitating comparison of future harmonic balance methods.
- The automatically computed tangent blocks could feed other Newton variants, such as quasi-Newton or matrix-free implementations, as preconditioners, extending the approach toward systems with millions of unknowns.
Reading between the lines
- The localization assumption (Eq. 47) is untested for globally coupled nonlinearities; if a system has dense nonlinear coupling (e.g., a fully populated cubic stiffness matrix), the automatic-differentiation graph would grow with the global DOF count and the flat GPU-memory plateau would break. A direct test would run pyHB on an all-to-all coupled 2000-DOF chain and watch memory.
- The reported speedup is against a 10-period time integration that already starts from the harmonic balance solution; for a full operating-domain sweep the time-integration baseline would likely be far slower, so the practical speedup may exceed the paper's stated number, but a fair comparison would report total wall-clock per complete branch, not just per continuation point.
- The modular architecture suggests pyHB's localized Jacobian could serve as a preconditioner in matrix-free Newton-Krylov implementations, which the paper itself lists as future work but does not test.
- The implicit-trapezoidal monodromy integration requires inverting a 2n-by-2n matrix per time segment; whether that remains efficient at lower sparsity or higher DOF counts is untested and could become a bottleneck for stability analysis of very large systems.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents pyHB, an open-source Python implementation of the Harmonic Balance (HB) method for periodic response analysis of nonlinear dynamical systems. The central methodological contribution is the use of PyTorch automatic differentiation (AD) to compute the Jacobian of the nonlinear force, applied only to a reduced, localized representation of that force (Eq. (47), Sections 3.4 and 3.5). The authors supplement this with FFT-based residual projection, sparse assembly of the linear Jacobian part, a blocked solver for the arc-length continuation system, weighted arc-length scaling, and Floquet stability analysis via an implicit trapezoidal integrator. Four numerical examples are presented: a quasi-zero-stiffness isolator, a piezoelectric-magnetic energy harvester, a 284-DOF aeroengine dual-rotor model, and a 2000-DOF Bernoulli beam with a localized nonlinear spring. The results are verified against explicit Runge-Kutta or Newmark-β time integration for stable branches, and unstable branches are classified via computed Floquet multipliers. The headline performance claim is that the AD-enhanced solver reaches about 0.44 s per continuation point at 202,000 HB unknowns with 637.8 MB additional RAM and 243.5 MB GPU memory, while HB-AD (the authors' prior method) runs out of memory.
Significance. If the claims hold, pyHB would be a useful and reproducible engineering tool: it removes the need for manual derivation of nonlinear-force Jacobians in HB analysis, provides a modular workflow from model definition through stability assessment, and demonstrates good scalability on a genuinely high-dimensional example (2000 DOFs, 202,000 HB unknowns). The comparison with the authors' previous HB-AD implementation is informative and shows a substantive improvement in GPU memory footprint. The open-source release, the inclusion of executable examples, and the verification against independent time-integration baselines are concrete strengths that support the credibility of the reported results. The main significance risk is that the generality and scalability claims are conditioned on a localization assumption that is not stated as a limitation and is not tested in the distributed-nonlinearity regime.
major comments (3)
- [§3.4, Eq. (47), and §5/Table 1] The scalability and memory claims rest on the localization assumption f_nl = Ξ f_nl^(r)(Θx'', Θx', Θx, τ, ω) with q, r << n. The paper states in §3.5 that the AD computational graph 'depends mainly on the local dimensions q and r', and Table 1 shows GPU memory staying near 243 MB precisely because the nonlinearity is localized. For distributed nonlinearities with q, r = O(n), the reduced block H_nl^(r) and the AD graph scale with n, so the claimed 'controllable GPU memory' and near-linear scaling would not hold. This is not flagged as a scope limitation in §6; the abstract's 'general user-defined nonlinear systems' and 'one-stop benchmark platform' therefore overstate the demonstrated capability. A caveat and a small distributed-nonlinearity experiment or complexity estimate are needed.
- [§3.5 and §4.3 (Eqs. (62)-(63))] The paper advertises support for non-smooth nonlinear forces (gaps, piecewise stiffness, Heaviside functions) and asserts that PyTorch's subgradients provide a 'robust solution'. However, for Newton-HB iterations the relevant object is a consistent tangent of the projected residual; for discontinuous or non-differentiable forces a subgradient is not automatically a valid Newton derivative. The aeroengine example uses a Heaviside function and a fractional exponent, yet no convergence diagnostics, mesh-refinement study, or comparison of subgradient-based Jacobians against a smoothed approximation is reported. As it stands, the claim that non-smooth forces are handled 'without additional manual processing' is not supported by evidence. At minimum, the paper should report behavior at switching points and discuss conditions under which the AD subgradient yields a convergent HB iteration.
- [§4.4 and §5, Table 1] The performance comparison with Newmark-β is presented as evidence of a 'several-hundred-fold speedup'. The comparison is per continuation point versus a 10-period Newmark trajectory initialized directly from the HB solution. This is a legitimate pointwise cost comparison, but it does not compare full workflows (e.g., obtaining a complete response curve with continuation versus a swept time-integration run), and the text in §4.4 even contrasts '15 s for the complete pyHB curve' with '>80 s for 10 periods at one ω'. The abstract's 'several-hundred-fold' number comes from the per-point ratio, which should be stated explicitly. I recommend clarifying the comparison basis in the text so that readers do not infer a total-workflow speedup of that magnitude.
minor comments (6)
- [Eq. (17)] In the expression for ̄C(τ), the term ∂f_nl/∂x'' should presumably be ∂f_nl/∂x'. As written, the perturbation damping term contains a derivative with respect to x'', which appears twice and is inconsistent with the preceding linearization.
- [Throughout] Several typographical errors should be corrected: 'baed' (p. 10), 'Simlarly' (p. 8), 'discribed' (p. 7), 'anbitrary' (p. 9), 'utlized' (p. 17), 'omega' (p. 17), and inconsistent spellings of 'preprocessing/preparation'.
- [§3.6] The weighted arc-length scaling factors q_s and ω_s are introduced as user choices with no guidance or sensitivity study. Since the continuation behavior and even convergence can depend on these factors, an automatic or recommended scaling procedure would strengthen reproducibility.
- [§3.9 and §4.1] The frequency resolution parameter F_r is described as enabling subharmonic analysis, and Example 4.1 uses F_r = 1/3. The paper should explain more concretely how F_r, the period T, and the number of FFT samples N1 interact, and whether the harmonic basis is always commensurate.
- [§5] The comparison with COCO and NLvib cites the earlier computational environment in reference [128]. This is acceptable as a reference point, but the sentence 'the computational cost of both methods exceeds 24 hours' should be accompanied by the original hardware/software conditions to avoid implying a same-machine benchmark.
- [§1 and §4.3] The statement that 'the explicit Runge-Kutta method adopted in the previous examples is unstable for this system' is imprecise: DOP853 is not globally unstable; rather it would require excessively small steps for this stiff (contact-bearing) problem. A more careful wording would avoid an unfair characterization of explicit integrators.
Circularity Check
No significant circularity: the HB derivations are self-contained and the key claims are validated against independent time-integration baselines.
full rationale
The paper's central derivation chain is not circular. The HB residual and Jacobian construction in Eqs. (10)-(13), (27)-(54) follows standard harmonic-balance algebra: the linear part is precomputed from the system matrices, and the nonlinear part is formed by FFT-based projection, tensor contraction, and AD of the user-supplied local force f_nl^(r). No target result is inserted into the derivation as an assumption; the AD procedure (Alg. 1) computes derivatives of the user-defined nonlinear force rather than fitting the response. The validation examples compare pyHB against independent Runge-Kutta and Newmark-β time integration, and the reported speedups and memory figures are measured benchmarks, not quantities implied by the formulation. The unstable branches, which time integration cannot produce, are identified by Floquet analysis, so agreement on stable branches is not forced by construction. Self-citations, notably the HB-AD comparison with reference [70], are used as empirical baselines rerun under identical conditions; this is code-reproduced evidence, not an appeal to an unverified prior result. The localized-nonlinearity premise in Eq. (47) is a stated modeling assumption that limits the claimed generality, but it is an input assumption rather than a circular step: the paper does not define the efficiency result in terms of that premise, and the derivation remains valid for the systems satisfying it. Overall, no load-bearing step reduces to its own inputs or to a self-citation chain, so the appropriate circularity score is 0.
Assumptions & free parameters
free parameters (2)
- Arc-length scaling factors q_s and ω_s
- Frequency resolution F_r
assumptions (4)
- domain assumption Nonlinear forces are localized to a small subset of DOFs (q,r << n) as in Eq. (47)
- ad hoc to paper PyTorch's subgradients provide valid Jacobians for non-smooth nonlinear forces
- domain assumption The implicit trapezoidal integrator gives a faithful monodromy matrix for Floquet analysis
- domain assumption The system is periodically excited and the steady-state response is periodic with commensurate frequencies
Cite this review
Pith. "Pith review of pyHB: an open-source automatic-differentiation-enhanced semi-analytical solver for nonlinear dynamics." pith.science (2026). https://pith.science/paper/C24N3O2X
@misc{pith2026260717577,
author = {Pith},
title = {Pith review of: pyHB: an open-source automatic-differentiation-enhanced semi-analytical solver for nonlinear dynamics},
year = {2026},
howpublished = {\url{https://pith.science/paper/C24N3O2X}},
note = {Machine review of arXiv:2607.17577}
}
abstract
The Harmonic Balance (HB) method is widely used to compute and analyze the periodic responses of nonlinear systems. However, its application to high-dimensional complex systems is limited by the burden of handling the partial derivatives of the nonlinearities. This work presents pyHB, an open-source, automatic-differentiation-enhanced semi-analytical framework that integrates the complete HB workflow for general user-defined nonlinear systems. The proposed formulation exploits localized nonlinearities and applies PyTorch-based automatic differentiation (AD) only to the reduced nonlinear force, thereby avoiding the need for user-supplied derivatives of the nonlinear force and maintaining controllable GPU memory usage. Weighted arc-length continuation, sparse matrix assembly, a blocked solution strategy for the augmented continuation equations, and Floquet-based stability analysis are incorporated within a modular architecture that separates model definition from reusable numerical procedures. Hence, pyHB can provide a complete landscape of the nonlinear system's periodic response based solely on the user-defined dynamical equations. Four examples, including a quasi-zero-stiffness isolator, a nonlinear piezoelectric energy harvester, a 284 degrees of freedom (DOFs) aeroengine model, and a 2000 DOFs Bernoulli beam, demonstrate the ability of pyHB to trace stable and unstable solution branches and capture subharmonic resonance, combination resonance, and mixed-order electromechanical responses. Notably, in the Bernoulli beam example with 202000 HB unknowns, the AD-enhanced solver requires approximately 0.44s per continuation point, achieving several-hundred-fold speedup compared to the Newmark-$\beta$ method and remaining 637.8MB of additional RAM and 243.5MB of GPU memory. The proposed pyHB provides a general, one-stop benchmark platform for HB-based nonlinear dynamics analysis.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
S. Mei, C. Caprani, D. Cantero, Dynamic response of a combined beam and spring-mass system under moving load: A closed-form solution using higher-order ordinary differential equations, Structures 87 (2026) 111521
2026
-
[2]
Zhang, Y
Y. Zhang, Y. Teng, H. Feng, W. Zhang, Nonlinear vibration analysis and experimental research on rotating pre-twisted tc4 blades under combined transverse aerodynamic loads and tip clearance airflow, Thin-Walled Structures 228 (2026) 115115
2026
-
[3]
Theodosiou, K
C. Theodosiou, K. Sikelis, S. Natsiavas, Periodic steady state response of large scale mechanical models with local nonlinearities, International Journal of Solids and Structures 46 (2009) 3565–3576
2009
-
[4]
Su, Y.-F
R.-X. Su, Y.-F. Guo, Y.-P. Zhao, Steady-state dynamical behavior analysis of an improved fitzhugh-nagumo neuron system based on the generalized cell mapping method, Chinese Journal of Physics 103 (2026) 1–17
2026
-
[5]
Y. Dong, H. Hu, L. Wang, A comprehensive study on the coupled multi-mode vibrations of cylindrical shells, Mechanical Systems and Signal Processing 169 (2022) 108730
2022
-
[6]
Karam, J
M. Karam, J. C. Sutherland, T. Saad, Low-cost Runge-Kutta integrators for incompressible flow simulations, Journal of Computational Physics 443 (2021) 110518
2021
-
[7]
K. Wang, C. Tang, D. Zhao, X. Qian, S. Jiang, Simulation of transient heat conduction in granular materials using the numerical manifold method with a high-accuracy explicit time integration scheme, Journal of Computational Physics 564 (2026) 115162
2026
-
[8]
Sokolov, P
A. Sokolov, P. Barkhayev, S. Turek, Numerical study of the rbf-fd parallel-in-time contour integration method for convectiondiffusion equations, Journal of Computational and Applied Mathematics 475 (2026) 117017
2026
Show all 137 references
-
[9]
Zhang, J
Q. Zhang, J. Zhang, Z. zhong Sun, Optimal convergence rate of the explicit euler method for convection-diffusion equations, Applied Mathematics Letters 131 (2022) 108048
2022
-
[10]
improved euler
O. Hanna, New explicit and implicit "improved euler" methods for the integration of ordinary differential equations, Computers & Chemical Engineering 12 (1988) 1083–1086
1988
-
[11]
Faleichik, A
B. Faleichik, A. Moisa, Explicit runge-kutta-chebyshev methods of second order with monotonic stability polynomial, Journal of Computational and Applied Mathematics 476 (2026) 117061
2026
-
[12]
L. Mei, Y. Chen, Explicit multistep method for the numerical solution of rlw equation, Applied Mathematics and Computation 218 (2012) 9547–9554
2012
-
[13]
Chang, L
Z. Chang, L. Hou, R. Lin, Y. Chen, P. Masarati, Monolithic approaches to transient thermo-mechanical interaction in nonlinear rotor systems, International Journal of Mechanical Sciences 290 (2025) 110066
2025
-
[14]
X. Qin, J. Yu, Z. Jiang, L. Huang, C. Yan, Explicit strong stability preserving second derivative multistep methods for the euler and navier- stokes equations, Computers & Fluids 268 (2024) 106089
2024
-
[15]
El Abbassi, A
A. El Abbassi, A. Bichri, Y. Achenani, L. Souinida, A. Tilioua, Control of a nonlinear dynamic system by backstepping: Numerical simulation and modelling using electronic circuits, Procedia Computer Science 265 (2025) 728–733
2025
-
[16]
J. Li, K. Yu, X. Li, An identical second-order single step explicit integration algorithm with dissipation control for structural dynamics, International Journal for Numerical Methods in Engineering 122 (2021) 1089–1132
2021
-
[17]
T. Liu, C. Zhao, Q. Li, L. Zhang, An efficient backward euler time-integration method for nonlinear dynamic analysis of structures, Computers & Structures 106-107 (2012) 20–28
2012
-
[18]
Y. Gui, R. Du, Second-order semi-implicit imex runge-kutta methods for micromagnetic simulations in antiferromagnetics, Communications in Nonlinear Science and Numerical Simulation 161 (2026) 110088
2026
-
[19]
Z. Li, H. Du, J. Si, Z. Wang, W. Xiong, Rigid-flexible coupling dynamic modeling and validation of a helicopter rescue simulator based on an inverted Stewart platform, Ocean Engineering 319 (2025) 120164
2025
-
[20]
B. Liu, Y. Ji, Nonlinear dynamics of three-dimensional curved geometrically exact beams by a quadrature element formulation, Nonlinear Dynamics 112 (2024) 14925–14958
2024
-
[21]
Chowdhury, S
S. Chowdhury, S. Adhikari, Nonlinear stiffened inertial amplifier tuned mass friction dampers, Soil Dynamics and Earthquake Engineering 191 (2025) 109264
2025
-
[22]
R. Chen, Y. Lü, R. Hu, Y. Zhang, G. Liu, Z. Xie, C. Li, Nonlinear dynamic characteristics of flexible bolted rotor system supported by journal bearings considering time varying bending stiffness of flange interface, Mechanical Systems and Signal Processing 245 (2026) 113858
2026
-
[23]
F. Wang, J. Gao, L. liang Xu, Y. fang Zheng, C. ping Chen, Analysis of nonlinear dynamic response of magneto-electro-elastic laminated beams with weak interface bonding, Communications in Nonlinear Science and Numerical Simulation 162 (2026) 110403
2026
-
[24]
X. Yin, Y. Quan, L. Wu, T. Kaiersaer, Z. Huang, A 3d vehicle-bridge interaction framework integrating energy-conserving hamilton’s principle and stabilized newmark- method, Advances in Engineering Software 211 (2025) 104022
2025
-
[25]
Piterskaya, M
A. Piterskaya, M. Mortensen, A comparison of implicit-explicit runge-kutta time integration schemes in numerical solvers based on the galerkin and petrov-galerkin spectral methods for two-dimensional magneto-hydrodynamic problems, Results in Applied Mathematics 29 (2026) 100683
2026
-
[26]
X. Yin, Y. Quan, Y. Zhou, Z. Huang, J. Liu, H. Tang, P. Wei, A substructure method with coupled explicit-implicit strategy for long-span wave-bridge fsi analysis, Ocean Engineering 363 (2026) 126751
2026
-
[27]
Zhang, Y
J. Zhang, Y. Zhang, S. Chen, G. Hong, A hybrid implicit-explicit time integration for stiff chemically reacting flows based on adaptive component-splitting method, Computers & Fluids 307 (2026) 106977
2026
-
[28]
Cardona, T
A. Cardona, T. Coune, A. Lerusse, M. Geradin, A multiharmonic method for non-linear vibration analysis, International Journal for Numerical Methods in Engineering 37 (1994) 1593–1608
1994
-
[29]
Q. Ke, H. Wang, Z. Chen, J. Li, Y. Lin, Quantitative analysis of limit cycles in two-stroke oscillators with exponential functions based on Perturbation Incremental Method, Chaos, Solitons & Fractals 188 (2024) 115549
2024
-
[30]
X. Su, H. Kang, W. Zhang, T. Hu, Y. Cong, T. Liu, Nonlinear coupled vibrations of a resonant two-cable system with springs under the support motion, Applied Mathematical Modelling 161 (2027) 117138. Yuhong Jin et al.: Preprint submitted to Elsevier Page 28 of 32 pyHB: an open-...
2027
-
[31]
Y. Kang, Z. Zhang, X. Zhang, Z. Zhang, Q. Liu, Investigation of a tuned clutch inerter mass damper for vibration reduction: An ANN-based harmonic balance method and parameter optimization, Applied Mathematical Modelling 149 (2026) 116321
2026
-
[32]
Chang, L
Z. Chang, L. Hou, Y. Chen, Investigation on the 1:2 internal resonance of an FGM blade, Nonlinear Dynamics 107 (2022) 1937–1964
2022
-
[33]
Z. Li, G. Li, S. Lan, H. Wang, Y. Hu, Study on superharmonic-parametric combined resonance and multistability of translating piezoelectric sandwich plate, Thin-Walled Structures 231 (2026) 115370
2026
-
[34]
Z. Wu, Y. Zhang, G. Yao, 3/2 superharmonic resonance and 1/2 subharmonic resonance of functionally graded carbon nanotube reinforced composite beams, Composite Structures 241 (2020) 112056
2020
-
[35]
R. Lin, L. Hou, S. Zhong, Y. Chen, Nonlinear vibration and stability analysis of a dual-disk rotor-bearing system under multiple frequency excitations, Nonlinear Dynamics 112 (2024) 12815–12846
2024
-
[36]
Y. Li, J. Huang, F. Liao, W. Zhu, Coupled-mode subharmonic resonance of a piecewise-linear gear transmission system with 1:2 internal resonance, Mechanical Systems and Signal Processing 240 (2025) 113383
2025
-
[37]
M. Xu, B. Cai, C. Li, H. Zhang, Z. Liu, D. He, Y. Zhang, Dynamic characteristics and reliability analysis of ball screw feed system on a lathe, Mechanism and Machine Theory 150 (2020) 103890
2020
-
[38]
X.-D. Yang, W. Zhang, Nonlinear dynamics of axially moving beam with coupled longitudinal–transversal vibrations, Nonlinear Dynamics 78 (2014) 2547–2556
2014
-
[39]
S. Wang, L. Hou, Q. Meng, G. Cui, X. Wang, Three-magnet-ring quasi-zero stiffness isolator for low-frequency vibration isolation, International Journal of Mechanical System Dynamics 4 (2024) 153–170
2024
-
[40]
M. A. AL-Shudeifat, E. A. Butcher, C. R. Stern, General harmonic balance solution of a cracked rotor-bearing-disk system for harmonic and sub-harmonic analysis: Analytical and experimental approach, International Journal of Engineering Science 48 (2010) 921–935
2010
-
[41]
Mohamed, N
H. Mohamed, N. Brödling, F. Duddeck, Autoencoder-based reduction of harmonic balance Fourier coefficients for uncertainty analysis of multi-DOF gear systems, Results in Engineering 28 (2025) 107902
2025
-
[42]
S. L. Lau, Y. K. Cheung, Amplitude Incremental Variational Principle for Nonlinear Vibration of Elastic Systems, Journal of Applied Mechanics 48 (1981) 959–964
1981
-
[43]
Y. B. Kim, S. T. Noah, Stability and Bifurcation Analysis of Oscillators With Piecewise-Linear Characteristics: A General Approach, Journal of Applied Mechanics 58 (1991) 545–553
1991
-
[44]
Chopra, J
I. Chopra, J. Dugundji, Non-linear dynamic response of a wind turbine blade, Journal of Sound and Vibration 63 (1979) 265–286
1979
-
[45]
Z. Yan, H. Dai, Q. Wang, S. N. Atluri, Harmonic Balance Methods: A Review and Recent Developments, Computer Modeling in Engineering & Sciences 137 (2023) 1419–1459
2023
-
[46]
Zhang, K
K. Zhang, K. Lu, S. Chai, H. Cheng, C. Fu, D. Guo, Dynamic modeling and parameter sensitivity analysis of AUV by using the POD method and the HB-AFT method, Ocean Engineering 293 (2024) 116693
2024
-
[47]
Litewka, R
P. Litewka, R. Lewandowski, Non-linear harmonic vibrations of laminate plates with viscoelastic layers described by extended fractional models, Thin-Walled Structures 228 (2026) 115091
2026
-
[48]
Chang, L
Z. Chang, L. Hou, P. Masarati, R. Lin, Z. Li, Y. Chen, Modeling and nonlinear analysis of a coupled thermo-mechanical dual-rotor system, Nonlinear Dynamics 112 (2024) 17811–17842
2024
-
[49]
R. Lin, L. Hou, S. Dun, Y. Cai, C. Sun, Y. Chen, Synchronous impact phenomenon of a high-dimension complex nonlinear dual-rotor system subjected to multi-frequency excitations, Science China Technological Sciences 66 (2023) 1757–1768
2023
-
[50]
Q. Meng, L. Hou, R. Lin, Y. Chen, N. A. Saeed, A. Fouly, E. Awwad, On a quasi-zero stiffness vibration isolator with multiple zero stiffness points for mass load deviation, Applied Mathematical Modelling 145 (2025) 116112
2025
-
[51]
H. Lei, H. Ren, P. Zhou, F. Guo, Y. Wang, Embedding intelligent excitation adaptability for resonance-suppressed QZS vibration isolation, International Journal of Mechanical Sciences 307 (2025) 110931
2025
-
[52]
W. Wu, Y. Wang, S. Xiao, Q. Wu, Y. Feng, D. Cao, H. Han, Rotational-inertia-enhanced quasi-zero-stiffness chiral metamaterial for high- performance low-frequency vibration isolation, International Journal of Mechanical Sciences 324 (2026) 111787
2026
-
[53]
Wei, Y.-Z
L.-S. Wei, Y.-Z. Wang, Y.-S. Wang, Nonreciprocal transmission of nonlinear elastic wave metamaterials by incremental harmonic balance method, International Journal of Mechanical Sciences 173 (2020) 105433
2020
-
[54]
Y. Chen, B. Gong, Z. Zheng, On the subcritical period doubling of a non-smooth network system by incremental harmonic balance method, International Journal of Non-Linear Mechanics 148 (2023) 104256
2023
-
[55]
R. Jin, Q. Liu, Y. Ma, X. Shen, Y. Wang, H. Qian, D. Jiang, Physics-informed harmonic balance identification of high-dimensional structures with nonlinear stiffness and damping, International Journal of Non-Linear Mechanics 181 (2026) 105282
2026
-
[56]
H. Dai, X. Yue, J. Yuan, S. N. Atluri, A time domain collocation method for studying the aeroelasticity of a two dimensional airfoil with a structural nonlinearity, Journal of Computational Physics 270 (2014) 214–237
2014
-
[57]
H. Dai, Z. Yan, X. Wang, X. Yue, S. N. Atluri, Collocation-based harmonic balance framework for highly accurate periodic solution of nonlinear dynamical system, International Journal for Numerical Methods in Engineering 124 (2023) 458–481
2023
-
[58]
Liang, C
H. Liang, C. Zang, X. Wei, H. Wang, An integrated approach for response prediction in large-scale rotor-bearing system with local nonlinear joints based on FRF-based harmonic balance method, Journal of Sound and Vibration 583 (2024) 118450
2024
-
[59]
N. Chen, S. Cao, Y. Hou, An indirect harmonic balance method based on frequency response functions simplification for periodical response analysis of local nonlinearity systems, Computers & Structures 310 (2025) 107663
2025
-
[60]
Colaïtis, A
Y. Colaïtis, A. Batailly, The harmonic balance method with arc-length continuation in blade-tip/casing contact problems, Journal of Sound and Vibration 502 (2021) 116070
2021
-
[61]
Y. Li, J. Huang, W. Zhu, A generalized incremental harmonic balance method by combining a data-driven framework for initial value selection of strongly nonlinear dynamic systems, International Journal of Non-Linear Mechanics 169 (2025) 104951
2025
-
[62]
Y. L. Li, J. L. Huang, W. D. Zhu, An enhanced incremental harmonic balance method to improve the computational efficiency and convergence for systems with non-polynomial nonlinearities, Nonlinear Dynamics 113 (2025) 8265–8294. Yuhong Jin et al.: Preprint submitted to Elsevier Pa...
2025
-
[63]
Q. Wu, Y. Wang, D. Cao, Improving the Solution Procedure of Incremental Harmonic Balance Method for Multi-Degree-of-Freedom Self-Excited Vibration Systems, International Journal of Applied Mechanics 17 (2025) 2450125
2025
-
[64]
R. Ju, W. Fan, W. D. Zhu, An efficient Galerkin averaging-incremental harmonic balance method for nonlinear dynamic analysis of rigid multibody systems governed by differential–algebraic equations, Nonlinear Dynamics 105 (2021) 475–498
2021
-
[65]
X. Zhou, A. Heya, T. Inoue, A constraint elimination-based harmonic balance method for multibody systems, International Journal of Mechanical Sciences 317 (2026) 111490
2026
-
[66]
Huang, T
J. Huang, T. Wang, W. Zhu, An incremental harmonic balance method with two time-scales for quasi-periodic responses of a van der pol-mathieu equation, International Journal of Non-Linear Mechanics 135 (2021) 103767
2021
-
[67]
Y. Li, J. Huang, W. Zhu, A two-time-scale convergence-enhanced and efficient incremental harmonic balance method for solving quasi- periodic responses of nonlinear systems, Communications in Nonlinear Science and Numerical Simulation 152 (2026) 109042
2026
-
[68]
html, 2025
Ansys Inc., Ansys help, https://ansyshelp.ansys.com/public/account/secured?returnurl=/Views/Secured/main_page. html, 2025. Ansys documentation
2025
-
[69]
Malte, G
K. Malte, G. Johann, Nlvib-a matlab tool for nonlinear vibration problems., https://www.ila.uni-stuttgart.de/nlvib, 2019
2019
-
[70]
Y. Chen, Y. Jin, R. Lin, Y. Jiang, X. Mei, L. Hou, Y. Wang, N. T. Yong, A. Guo, Harmonic balance-automatic differentiation method: A practical nonlinear dynamics solver, International Journal of Mechanical Sciences 312 (2026) 111192
2026
-
[71]
Y. Chen, L. Hou, R. Lin, Y. Wang, N. A. Saeed, Y. Chen, Combination resonances of a dual-rotor-bearing-casing system, Nonlinear Dynamics 112 (2024) 4063–4083
2024
-
[72]
H. Wang, S. Zhou, X. Liang, T. Ding, T. Zhou, C. Dong, K. M. Lim, Dynamic Characteristics of Horizontal Rotor–Stator Nonlinear Interaction Under Additional Acceleration, AIAA Journal (2026) 1–20
2026
-
[73]
Arango Montoya, O
J. Arango Montoya, O. Chiello, J.-J. Sinou, R. Tufano, A Harmonic Balance Method with contact condensation for the frequency-domain computation of self-sustained nonlinear vibration related to railway curve squeal, Journal of Sound and Vibration 626 (2026) 119628
2026
-
[74]
G.-Y. Lee, K. Ahn, A computational framework based on the harmonic balance method for nonlinear bulk-flow analysis, Tribology International 214 (2026) 111336
2026
-
[75]
Zhang, H
H. Zhang, H. Xu, S. Wang, S. Zhou, Theoretical analysis of coupled thermo-electric-elastic rotational piezoelectric energy harvesters based on Green’s function and harmonic balance method, Applied Mathematical Modelling 139 (2025) 115815
2025
-
[76]
Saadatmand, J
M. Saadatmand, J. Kook, Nonlinear vibration analysis of circular multilayer graphene-based NEMS sensors using harmonic balance and pseudo-arclength continuation methods, International Journal of Non-Linear Mechanics 178 (2025) 105187
2025
-
[77]
D. Zhou, Z. Ji, J. Wu, E. Shen, T. Guo, Z. Lu, Multi-Passage harmonic balance method for flutter prediction in multi-row configurations, Aerospace Science and Technology 150 (2024) 109212
2024
-
[78]
Mohamed, N
H. Mohamed, N. Brödling, F. Duddeck, Multi-fidelity harmonic balance method: A non-intrusive approach, Mechanism and Machine Theory 217 (2025) 106251
2025
-
[79]
X. F. Wang, W. D. Zhu, A modified incremental harmonic balance method based on the fast Fourier transform and Broyden’s method, Nonlinear Dynamics 81 (2015) 981–989
2015
-
[80]
Huang, W
J. Huang, W. Zhou, W. Zhu, Quasi-periodic motions of high-dimensional nonlinear models of a translating beam with a stationary load subsystem under harmonic boundary excitation, Journal of Sound and Vibration 462 (2019) 114870
2019
-
[81]
T.-C. Yuan, J. Yang, L.-Q. Chen, A harmonic balance approach with alternating frequency/time domain progress for piezoelectric mechanical systems, Mechanical Systems and Signal Processing 120 (2019) 274–289
2019
-
[82]
Y. Sun, Z. Yin, W.-B. Shangguan, S. Rakheja, A computationally efficient method for dynamic response analyses of accessory drive systems using harmonic balance method together with alternating frequency/time domain technique, Mechanical Systems and Signal Processing 213 (2024) 111307
2024
-
[83]
Y. Chen, L. Hou, G. Chen, H. Song, R. Lin, Y. Jin, Y. Chen, Nonlinear dynamics analysis of a dual-rotor-bearing-casing system based on a modified HB-AFT method, Mechanical Systems and Signal Processing 185 (2023) 109805
2023
-
[84]
R. Ju, W. Fan, W. D. Zhu, Comparison Between the Incremental Harmonic Balance Method and Alternating Frequency/Time-Domain Method, Journal of Vibration and Acoustics 143 (2021) 024501
2021
-
[85]
R. Lin, L. Hou, Y. Chen, Y. Jin, N. A. Saeed, Y. Chen, A novel adaptive harmonic balance method with an asymptotic harmonic selection, Applied Mathematics and Mechanics 44 (2023) 1887–1910
2023
-
[86]
Speksnijder, U
A. Speksnijder, U. Karacadagli, H. Seyffert, A. Grammatikopoulos, Application of the harmonic balance method for ship-cargo interaction with intermittent contact nonlinearities, Journal of Sound and Vibration 601 (2025) 118925
2025
-
[87]
T. Li, Y. Cai, Modeling and dynamic analysis of belt drive systems using a harmonic balance and alternating frequency/time domain method, Engineering Computations 43 (2026) 452–482
2026
-
[88]
Zhang, W
J. Zhang, W. Zhang, Y. Zhang, Nonlinear vibrations of porous-hyperelastic cylindrical shell under harmonic force using harmonic balance and pseudo-arc length continuation methods, Thin-Walled Structures 198 (2024) 111767
2024
-
[89]
Woiwode, N
L. Woiwode, N. N. Balaji, J. Kappauf, F. Tubita, L. Guillot, C. Vergez, B. Cochelin, A. Grolet, M. Krack, Comparison of two algorithms for Harmonic Balance and path continuation, Mechanical Systems and Signal Processing 136 (2020) 106503
2020
-
[90]
Lee, Y.-H
G.-Y. Lee, Y.-H. Park, A proper generalized decomposition-based harmonic balance method with arc-length continuation for nonlinear frequency response analysis, Computers & Structures 275 (2023) 106913
2023
-
[91]
Q. Wu, Y. Wang, D. Cao, Nonlinear dynamic analysis of high aspect ratio wings via IHB method, Nonlinear Dynamics 113 (2025) 16225– 16244
2025
-
[92]
F. Guo, H. Ren, Y. Yao, W. Fan, R. Ju, P. Zhou, Semi-analytical Hopf bifurcation analysis of self-excited hunting motion in railway vehicles, Applied Mathematical Modelling (2026) 117144
2026
-
[93]
B. Zhou, F. Thouverez, D. Lenoir, A variable-coefficient harmonic balance method for the prediction of quasi-periodic response in nonlinear systems, Mechanical Systems and Signal Processing 64–65 (2015) 233–244. Yuhong Jin et al.: Preprint submitted to Elsevier Page 30 of 32 pyH...
2015
-
[94]
Huang, T
J. Huang, T. Wang, W. Zhu, An incremental harmonic balance method with two time-scales for quasi-periodic responses of a Van der Pol–Mathieu equation, International Journal of Non-Linear Mechanics 135 (2021) 103767
2021
-
[95]
Klinger, F
T. Klinger, F. Greiner, A. Rohde, A. Piel, M. E. Koepke, Van der Pol behavior of relaxation oscillations in a periodically driven thermionic discharge, Physical Review E 52 (1995) 4316–4327
1995
-
[96]
Seydel, Practical Bifurcation and Stability Analysis, volume 5 of Interdisciplinary Applied Mathematics, Springer New York, New York, NY, 2010
R. Seydel, Practical Bifurcation and Stability Analysis, volume 5 of Interdisciplinary Applied Mathematics, Springer New York, New York, NY, 2010. doi:10.1007/978-1-4419-1740-9
2010 doi
-
[97]
E. J. Doedel, R. C. Paenroth, A. R. Champneys, T. F. Fairgrieve, Auto 2000 : Continuation and bifurcation software for ordinary differential equations (with homcont), 1997. URL: https://api.semanticscholar.org/CorpusID:117201224
2000
-
[98]
Taege, H
S. Taege, H. Schiefhauer, A. Ali, B. Müller, Floquet theory based stability analysis of low-speed sensorless control, IFAC-PapersOnLine 59 (2025) 680–685
2025
-
[99]
Y. Wang, C. Cai, C. Liao, Z. Hu, L. Zhang, X. Sun, X. Zhong, Q. Li, Aeroelastic stability analysis of large-scale wind turbine blades under different operating conditions based on system identification and floquet theory, Energy 326 (2025) 136279
2025
-
[100]
Riahi, M
M. Riahi, M. Hayani Choujaa, S. Wang, Flow reversal and neimarksacker bifurcations in zero-mean time-modulated temperatures in a hele-shaw cell: Floquet theory and wentzelkramersbrillouin analysis, International Communications in Heat and Mass Transfer 169 (2025) 109507
2025
-
[101]
Ricci, P
R. Ricci, P. Pennacchi, Discussion of the dynamic stability of a multi-degree-of-freedom rotor system affected by a transverse crack, Mechanism and Machine Theory 58 (2012) 82–100
2012
-
[102]
C. S. Hsu, W. Cheng, Applications of the theory of impulsive parametric excitation and new treatments of general parametric excitation problems, Journal of Applied Mechanics 40 (1973) 78–86
1973
-
[103]
Cheung, S
Y. Cheung, S. Chen, S. Lau, Application of the incremental harmonic balance method to cubic non-linearity systems, Journal of Sound and Vibration 140 (1990) 273–286
1990
-
[104]
D. Yang, J. Huang, W. Zhu, Nonlinear dynamics of a misaligned and eccentric rotor system with three base motions, Applied Mathematical Modelling 148 (2025) 116240
2025
-
[105]
Jiang, K
W. Jiang, K. Liu, X. Yuan, H. Cao, J. Shi, Q. Qin, Nonlinear dynamics of rotor-support-casing system with support looseness fault, International Journal of Mechanical Sciences 300 (2025) 110482
2025
-
[106]
Y. Yao, X. Huang, X. Yang, L. Jézéquel, Nonlinear dynamics analysis of rotor systems based on extended nonlinear modal theory with gyroscopic coupling, Chinese Journal of Aeronautics (2026) 104180
2026
-
[107]
Baglioni, F
S. Baglioni, F. Cianetti, C. Braccesi, D. M. De Micheli, Multibody modelling of N DOF robot arm assigned to milling manufacturing: Dynamic analysis and position errors evaluation, Journal of Mechanical Science and Technology 30 (2016) 405–420
2016
-
[108]
L. Cao, Y. Chen, L. Wang, L. Salles, Z. Zheng, Perturbation Function Iteration Method: A New Framework for Solving Periodic Solutions of Non-linear and Non-smooth Systems, 2025. arXiv:2510.23071
2025
-
[109]
Deshpande, S
S. Deshpande, S. Mehta, G. N. Jazar, Optimization of secondary suspension of piecewise linear vibration isolation systems, International Journal of Mechanical Sciences 48 (2006) 341–377
2006
-
[110]
R. D. O. Teloli, S. Da Silva, A new way for harmonic probing of hysteretic systems through nonlinear smooth operators, Mechanical Systems and Signal Processing 121 (2019) 856–875
2019
-
[111]
S. Wang, Y. Zhang, W. Guo, T. Pi, X. Li, Vibration analysis of nonlinear damping systems by the discrete incremental harmonic balance method, Nonlinear Dynamics 111 (2023) 2009–2028
2023
-
[112]
S. Wang, L. Hua, C. Yang, X. Han, Z. Su, Applications of incremental harmonic balance method combined with equivalent piecewise linearization on vibrations of nonlinear stiffness systems, Journal of Sound and Vibration 441 (2019) 111–125
2019
-
[113]
R. E. Wengert, A simple automatic derivative evaluation program, Communications of the ACM 7 (1964) 463–464
1964
-
[114]
J. W. Hanson, J. S. Caviness, C. Joseph, Analytic differentiation by computer, Communications of the ACM 5 (1962) 349–355
1962
-
[115]
A. G. Baydin, B. A. Pearlmutter, A. A. Radul, J. M. Siskind, Automatic differentiation in machine learning: a survey, Journal of Machine Learning Research 18 (2017) 5595–5637
2017
-
[116]
Paszke, S
A. Paszke, S. Gross, F. Massa, A. Lerer, J. Bradbury, G. Chanan, T. Killeen, Z. Lin, N. Gimelshein, L. Antiga, A. Desmaison, A. Köpf, E. Yang, Z. DeVito, M. Raison, A. Tejani, S. Chilamkurthy, B. Steiner, L. Fang, J. Bai, S. Chintala, Pytorch: An imperative style, high-perform...
2019 arXiv
-
[117]
PyTorch 2.12 documentation, last updated Jan
PyTorch Contributors, Autograd mechanics, https://docs.pytorch.org/docs/2.12/notes/autograd.html, 2026. PyTorch 2.12 documentation, last updated Jan. 6, 2026; accessed Jul. 4, 2026
2026
-
[118]
Y. Chen, L. Hou, R. Lin, W. Toh, T. Ng, Y. Chen, A general and efficient harmonic balance method for nonlinear dynamic simulation, International Journal of Mechanical Sciences 276 (2024) 109388
2024
-
[119]
Riks, An incremental approach to the solution of snapping and buckling problems, International Journal of Solids and Structures 15 (1979) 529–551
E. Riks, An incremental approach to the solution of snapping and buckling problems, International Journal of Solids and Structures 15 (1979) 529–551
1979
-
[120]
snap-through
M. Crisfield, A fast incremental/iterative solution procedure that handles “snap-through”, Computers & Structures 13 (1981) 55–62
1981
-
[121]
E. J. Doedel, A. R. Champneys, F. Dercole, T. F. Fairgrieve, Y. A. Kuznetsov, B. E. Oldeman, R. C. Paffenroth, B. Sandstede, X. Wang, C. Zhang, AUTO-07P: Continuation and Bifurcation Software for Ordinary Differential Equations, Software manual, Concordia University, Montreal, C...
2007
-
[122]
Peletan, S
L. Peletan, S. Baguet, M. Torkhani, G. Jacquet-Richardet, A comparison of stability computational methods for periodic solution of nonlinear problems with application to rotordynamics, Nonlinear Dynamics 72 (2013) 671–682
2013
-
[123]
G. H. Gaonkar, D. S. Simha Prasad, D. Sastry, On Computing Floquet Transition Matrices of Rotorcraft, Journal of the American Helicopter Society 26 (1981) 56–61
1981
-
[124]
J. Fu, H. Miao, J. Hong, C. Li, X. Qu, Y. Ma, Y. Wang, Adaptive optimization design framework for vibration transfer suppression in aero-engine dual-rotor systems with shared bearing frame, Aerospace Science and Technology 176 (2026) 112475. Yuhong Jin et al.: Preprint submitt...
2026
-
[125]
H. Su, W. Zhou, L. Yan, N. Qiu, B. Gao, Theoretical modeling and experimental validation of a dual-disk rotor system with coupled spatial crack and bearing outer race defects, Mechanical Systems and Signal Processing 255 (2026) 114471
2026
-
[126]
P. Gao, Y. Chen, L. Hou, Nonlinear thermal behaviors of the inter-shaft bearing in a dual-rotor system subjected to the dynamic load, Nonlinear Dynamics 101 (2020) 191–209
2020
-
[127]
L. Tang, Y. Yang, A nonlinear piezoelectric energy harvester with magnetic oscillator, Applied Physics Letters 101 (2012) 094102
2012
-
[128]
Ponsioen, S
S. Ponsioen, S. Jain, G. Haller, Model reduction to spectral submanifolds and forced-response calculation in high-dimensional mechanical systems, Journal of Sound and Vibration 488 (2020) 115640
2020
-
[129]
Dankowicz, F
H. Dankowicz, F. Schilder, Recipes for Continuation, Society for Industrial and Applied Mathematics, Philadelphia, PA, 2013. doi:10.1137/ 1.9781611972573
2013
-
[130]
X. F. Wang, W. D. Zhu, Dynamic Analysis of an Automotive Belt-Drive System With a Noncircular Sprocket by a Modified Incremental Harmonic Balance Method, Journal of Vibration and Acoustics 139 (2017) 011009
2017
-
[131]
Veerse, Efficient iterative time preconditioners for harmonic balance rf circuit simulation, in: ICCAD-2003
F. Veerse, Efficient iterative time preconditioners for harmonic balance rf circuit simulation, in: ICCAD-2003. International Conference on Computer Aided Design (IEEE Cat. No.03CH37486), 2003, pp. 251–254. doi: 10.1109/ICCAD.2003.159698
2003
-
[132]
Lindblad, N
D. Lindblad, N. Andersson, Convergence Acceleration of the Harmonic Balance Method Using a Time-Level Preconditioner, AIAA Journal 58 (2020) 4908–4922
2020
-
[133]
C. Tan, Y. Su, F. Yang, X. Zeng, New time-domain preconditioners for hb jacobian of rf circuits, in: 2025 62nd ACM/IEEE Design Automation Conference (DAC), 2025, pp. 1–7. doi: 10.1109/DAC63849.2025.11132845
2025
-
[134]
R. J. Kuether, A. Steyer, Large-scale harmonic balance simulations with Krylov subspace and preconditioner recycling, Nonlinear Dynamics 112 (2024) 3377–3398
2024
-
[135]
S. Jain, T. Breunung, G. Haller, Fast computation of steady-state response for high-degree-of-freedom nonlinear systems, Nonlinear Dynamics 97 (2019) 313–341
2019
-
[136]
R. Ju, S. M. Yang, H. Ren, W. Fan, R. C. Ni, P. Gu, Steady-State Rotary Periodic Solutions of Rigid and Flexible Mechanisms With Large Spatial Rotations Using the Incremental Harmonic Balance Method for Differential-Algebraic Equations, Journal of Computational and Nonlinear Dy...
2024
-
[137]
X. Zhou, A. Heya, T. Inoue, A constraint elimination-based harmonic balance method for multibody systems, International Journal of Mechanical Sciences 317 (2026) 111490. Yuhong Jin et al.: Preprint submitted to Elsevier Page 32 of 32
2026
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.