Research Article | | Peer-Reviewed

Benchmarking Numerical Solvers for Damped Single-Degree-of-Freedom Vibration Systems: Technical Evaluation and Educational Insights

Received: 4 October 2025     Accepted: 14 October 2025     Published: 30 October 2025
Views:       Downloads:
Abstract

This study presents a systematic benchmarking of numerical methods for solving ordinary differential equations (ODEs) applied to damped single-degree-of-freedom (SDOF) vibration systems. Ten solvers—including Runge–Kutta variants, Adams–Bashforth–Moulton, Rosenbrock, and Backward Differentiation Formula (BDF)—were evaluated under both non-stiff and stiff conditions by varying mass, damping, and stiffness parameters. Analytical solutions were used as references to quantify global error, convergence behavior, and computational efficiency. High-order adaptive solvers, such as Verner’s and Runge–Kutta 7/8, consistently achieved the highest accuracy, reducing global error by up to 15% compared with the classical Runge–Kutta (ODE45) method. Implicit methods, including Rosenbrock and BDF, demonstrated superior stability in stiff and highly damped cases. In contrast, low-order approaches, particularly the Trapezoidal rule, exhibited the largest errors, exceeding 30% in oscillatory regimes. The results confirm that solver performance is problem-dependent, emphasizing that no single algorithm is universally optimal. Beyond the technical contributions, this study introduces a pedagogical framework that allows engineering students to visualize solver trade-offs, quantify numerical accuracy, and interpret computational efficiency. The educational integration strengthens conceptual understanding of numerical methods and supports data-driven solver selection in vibration analysis and related engineering applications.

Published in American Journal of Civil Engineering (Volume 13, Issue 5)
DOI 10.11648/j.ajce.20251305.15
Page(s) 304-312
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2025. Published by Science Publishing Group

Keywords

Numerical Methods, Ordinary Differential Equations (ODEs), Runge-Kutta Methods, Stiff Equations, Computational Efficiency, Single-Degree-of-Freedom (SDOF) Systems, Engineering Education

References
[1] Städter, P., Schälte, Y., Schmiester, L., Weindl, D., & Hasenauer, J. (2021). Benchmarking of numerical integration methods for ODE models of biological systems. Scientific Reports, 11, 2696.
[2] Owolabi, K. M., & Olayinka, O. A. (2023). Time‑accurate and highly‑stable explicit peer methods for stiff ODEs. Applied Mathematics and Computation, 449, 127903.
[3] Deka, P. J., & Einkemmer, L. (2021). Efficient adaptive step size control for exponential integrators. Journal of Computational Physics, 435, 110239.
[4] Gaudreault, S., Rainwater, G., & Tokman, M. (2018). KIOPS: A fast adaptive Krylov subspace solver for exponential integrators. Journal of Computational Physics, 372, 236–255.
[5] Choi, B., Bathe, K.-J., & Noh, G. (2022). Time‑splitting ratio in the ρ∞‑Bathe time integration method for higher‑order accuracy in structural dynamics and heat transfer. Computers & Structures, 270, 106814.
[6] Alhayki, E., & Dettmer, W. (2024). New implicit time integration schemes for structural dynamics combining high‑frequency damping and high second‑order accuracy. Computers & Structures, 303, 107587.
[7] Malakiyeh, M. M., Shojaee, S., Hamzehei‑Javaran, S., & Bathe, K.-J. (2023). The explicit β1/β2‑Bathe time integration method. Computers & Structures, 286, 107092.
[8] Bathe, K.-J. (2019). The Bathe time integration method with controllable spectral radius: The ρ∞‑Bathe method. Computers & Structures, 212, 289–298.
[9] Chung, J., & Hulbert, G. M. (1993). A time integration algorithm for structural dynamics with improved numerical dissipation: The generalized‑α method. Journal of Applied Mechanics, 60(2), 371–375.
[10] Jansen, K. E., Whiting, C. H., & Hulbert, G. M. (2000). A generalized‑α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method. Computer Methods in Applied Mechanics and Engineering, 190(3-4), 305–319.
[11] Cox, S. M., & Matthews, P. C. (2002). Exponential time differencing for stiff systems. Journal of Computational Physics, 176(2), 430–455.
[12] Hosea, M. E., & Shampine, L. F. (1996). Analysis and implementation of TR‑BDF2. Applied Numerical Mathematics, 20(1-2), 21–37.
[13] Simo, J. C., & Tarnow, N. (1992). The discrete energy‑momentum method: Conserving algorithms for nonlinear elastodynamics. Zeitschrift für Angewandte Mathematik und Physik, 43, 757–792.
[14] Bathe, K.-J., & Noh, G. (2012). Insight into an implicit time integration scheme for structural dynamics. Computers & Structures, 98–99, 1–6.
[15] Lopez, S. (2020). A predictor–corrector time integration algorithm for dynamic analysis of nonlinear systems. Nonlinear Dynamics, 102, 829–850.
[16] Lázaro, M., Rubio, H., & Morata, I. (2022). Damping perturbation‑based time integration asymptotic method for structural dynamics. International Journal of Applied Mechanics, 14(8), 2250022.
[17] Yang, C., Yang, B., Zhu, T., & Xiao, S. (2017). Comparison and assessment of time‑integration algorithms for nonlinear vibration systems. Journal of Central South University, 24, 1090–1097.
[18] Kim, W., & Park, T. (2020). A comparative study of implicit and explicit composite time integration schemes. International Journal of Structural Stability and Dynamics, 20(12), 2041003.
[19] Lopez, S. (2024). An explicit time integration method based on the Verlet scheme with improved characteristics in numerical dispersion. Journal of Engineering Mechanics, 150(6), 04024063.
Cite This Article
  • APA Style

    Cannon, J., Zirakian, T. (2025). Benchmarking Numerical Solvers for Damped Single-Degree-of-Freedom Vibration Systems: Technical Evaluation and Educational Insights. American Journal of Civil Engineering, 13(5), 304-312. https://doi.org/10.11648/j.ajce.20251305.15

    Copy | Download

    ACS Style

    Cannon, J.; Zirakian, T. Benchmarking Numerical Solvers for Damped Single-Degree-of-Freedom Vibration Systems: Technical Evaluation and Educational Insights. Am. J. Civ. Eng. 2025, 13(5), 304-312. doi: 10.11648/j.ajce.20251305.15

    Copy | Download

    AMA Style

    Cannon J, Zirakian T. Benchmarking Numerical Solvers for Damped Single-Degree-of-Freedom Vibration Systems: Technical Evaluation and Educational Insights. Am J Civ Eng. 2025;13(5):304-312. doi: 10.11648/j.ajce.20251305.15

    Copy | Download

  • @article{10.11648/j.ajce.20251305.15,
      author = {John Cannon and Tadeh Zirakian},
      title = {Benchmarking Numerical Solvers for Damped Single-Degree-of-Freedom Vibration Systems: Technical Evaluation and Educational Insights
    },
      journal = {American Journal of Civil Engineering},
      volume = {13},
      number = {5},
      pages = {304-312},
      doi = {10.11648/j.ajce.20251305.15},
      url = {https://doi.org/10.11648/j.ajce.20251305.15},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.ajce.20251305.15},
      abstract = {This study presents a systematic benchmarking of numerical methods for solving ordinary differential equations (ODEs) applied to damped single-degree-of-freedom (SDOF) vibration systems. Ten solvers—including Runge–Kutta variants, Adams–Bashforth–Moulton, Rosenbrock, and Backward Differentiation Formula (BDF)—were evaluated under both non-stiff and stiff conditions by varying mass, damping, and stiffness parameters. Analytical solutions were used as references to quantify global error, convergence behavior, and computational efficiency. High-order adaptive solvers, such as Verner’s and Runge–Kutta 7/8, consistently achieved the highest accuracy, reducing global error by up to 15% compared with the classical Runge–Kutta (ODE45) method. Implicit methods, including Rosenbrock and BDF, demonstrated superior stability in stiff and highly damped cases. In contrast, low-order approaches, particularly the Trapezoidal rule, exhibited the largest errors, exceeding 30% in oscillatory regimes. The results confirm that solver performance is problem-dependent, emphasizing that no single algorithm is universally optimal. Beyond the technical contributions, this study introduces a pedagogical framework that allows engineering students to visualize solver trade-offs, quantify numerical accuracy, and interpret computational efficiency. The educational integration strengthens conceptual understanding of numerical methods and supports data-driven solver selection in vibration analysis and related engineering applications.
    },
     year = {2025}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Benchmarking Numerical Solvers for Damped Single-Degree-of-Freedom Vibration Systems: Technical Evaluation and Educational Insights
    
    AU  - John Cannon
    AU  - Tadeh Zirakian
    Y1  - 2025/10/30
    PY  - 2025
    N1  - https://doi.org/10.11648/j.ajce.20251305.15
    DO  - 10.11648/j.ajce.20251305.15
    T2  - American Journal of Civil Engineering
    JF  - American Journal of Civil Engineering
    JO  - American Journal of Civil Engineering
    SP  - 304
    EP  - 312
    PB  - Science Publishing Group
    SN  - 2330-8737
    UR  - https://doi.org/10.11648/j.ajce.20251305.15
    AB  - This study presents a systematic benchmarking of numerical methods for solving ordinary differential equations (ODEs) applied to damped single-degree-of-freedom (SDOF) vibration systems. Ten solvers—including Runge–Kutta variants, Adams–Bashforth–Moulton, Rosenbrock, and Backward Differentiation Formula (BDF)—were evaluated under both non-stiff and stiff conditions by varying mass, damping, and stiffness parameters. Analytical solutions were used as references to quantify global error, convergence behavior, and computational efficiency. High-order adaptive solvers, such as Verner’s and Runge–Kutta 7/8, consistently achieved the highest accuracy, reducing global error by up to 15% compared with the classical Runge–Kutta (ODE45) method. Implicit methods, including Rosenbrock and BDF, demonstrated superior stability in stiff and highly damped cases. In contrast, low-order approaches, particularly the Trapezoidal rule, exhibited the largest errors, exceeding 30% in oscillatory regimes. The results confirm that solver performance is problem-dependent, emphasizing that no single algorithm is universally optimal. Beyond the technical contributions, this study introduces a pedagogical framework that allows engineering students to visualize solver trade-offs, quantify numerical accuracy, and interpret computational efficiency. The educational integration strengthens conceptual understanding of numerical methods and supports data-driven solver selection in vibration analysis and related engineering applications.
    
    VL  - 13
    IS  - 5
    ER  - 

    Copy | Download

Author Information
  • Sections