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 | 
Numerical Methods, Ordinary Differential Equations (ODEs), Runge-Kutta Methods, Stiff Equations, Computational Efficiency, Single-Degree-of-Freedom (SDOF) Systems, Engineering Education
| [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. | 
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
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
@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}
}
											
										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 -