García-Alonso, Fernando, Reyes, José Antonio, Cortés-Molina, Mónica
An Algorithm for the Numerical Integration of Perturbed and Damped Second-Order ODE Systems
García-Alonso F, Reyes JA, Cortés-Molina M. An Algorithm for the Numerical Integration of Perturbed and Damped Second-Order ODE Systems. Mathematics. 2020; 8(11):2028. https://doi.org/10.3390/math8112028
URI: http://hdl.handle.net/10045/110475
DOI: 10.3390/math8112028
ISSN: 2227-7390
Abstract: 
A new method of numerical integration for a perturbed and damped systems of linear second-order differential equations is presented. This new method, under certain conditions, integrates, without truncation error, the IVPs (initial value problems) of the type: x’’(t) + Ax’(t) + Cx(t) = εF(x(t), t), x(0) = x0, x’(0) = x’0, t ∈ [a, b] = I, which appear in structural dynamics, astrodynamics, and other fields of physics and engineering. In this article, a succession of real functions is constructed with values in the algebra of m x m matrices. Their properties are studied and we express the solution of the proposed IVP through a serial expansion of the same, whose coefficients are calculated by means of recurrences involving the perturbation function. This expression of the solution is used for the construction of the new numerical method. Three problems are solved by means of the new series method; we contrast the results obtained with the exact solution of the problem and with its first integral. In the first problem, a quasi-periodic orbit is integrated; in the second, a problem of structural dynamics associated with an earthquake is studied; in the third, an equatorial satellite problem when the perturbation comes from zonal harmonics J2 is solved. The good behavior of the series method is shown by comparing the results obtained against other integrators.
Keywords:Numerical solutions of second-order ODE systems, Perturbed and damped, Quasi-periodic orbit, Earthquake, Equatorial satellite
MDPI
info:eu-repo/semantics/article