The Hermite–Obreshkov–Padé (HOP) procedure is an implicit method for the numerical solution of a system of ordinary differential equations (ODEs) applicable to stiff dynamical systems. This procedure applies an Obreshkov condition to multiple derivatives of the system state vector, both at the start and end of a time step in the numerical solution. That condition is shown to be satisfied by the Hermite interpolating polynomial that matches the state vector and its derivatives, also at the start and end of a time step. The Hermite polynomial, in turn, can be specified in terms of the system state and its derivatives at the start of a step together with a collection of free parameters. Adjusting these free parameters to minimize magnitudes of the ODE residual and its derivatives at the end of a step serves as a proxy for matching the system state and its derivatives. A high-order Taylor expansion at the start of a time step interval models the residual and its derivatives over the entire interval. A variant of this procedure adjusts those parameters to match integrals of the system state over the duration of that interval. This is done by minimizing magnitudes of integrals of the ODE residual calculated from the extrapolating Taylor-series expansion, a process that avoids the need to determine integration constants for multiple integrals of the state. This alternative method eliminates the calculation of high-order derivatives of the system state and hence avoids loss in accuracy from floating-point round off. Numerical performance is evaluated on a dynamically unbalanced constant-velocity (CV) coupling having a high spring rate constraining shaft deflection.

## References

1.
Hairer
,
E.
, and
Wanner
,
G.
,
1991
,
Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems
,
Springer
,
Berlin
.
2.
Gear
,
C.
,
1971
, “
Simultaneous Numerical Solution of Differential-Algebraic Equations
,”
IEEE Trans. Circuit Theory
,
18
(
1
), pp.
89
95
.10.1109/TCT.1971.1083221
3.
Negrut
,
D.
,
Haug
,
E. J.
, and
German
,
H. C.
,
2003
, “
An Implicit Runge–Kutta Method for Integration of Differential Algebraic Equations of Multibody Dynamics
,”
Multibody Syst. Dyn.
,
9
(
2
), pp.
121
142
.10.1023/A:1022506312444
4.
Eriksson
,
K.
,
Johnson
,
C.
, and
Logg
,
A.
,
2004
, “
Explicit Time-Stepping for Stiff ODEs
,”
SIAM J. Sci. Comput. (USA)
,
25
(
4
), pp.
1142
1157
.10.1137/S1064827502409626
5.
Gear
,
C. W.
, and
Kevrekidis
,
I. G.
,
2003
, “
Projective Methods for Stiff Differential Equations: Problems With Gaps in Their Eigenvalue Spectrum
,”
SIAM J. Sci. Comput. (USA)
,
24
(
4
), pp.
1091
1106
.10.1137/S1064827501388157
6.
Milenkovic
,
P.
, “
Numerical Solution of Stiff Multi-Body Dynamic Systems Based on Kinematic Derivatives
,”
ASME J. Dyn. Syst., Meas., Control
(submitted).
7.
Griewank
,
A.
,
1996
, “
ODE Solving via Automatic Differentiation and Rational Prediction
,”
Numerical Analysis 1995
, Vol.
344
(Pitman Research Notes in Mathematics Series),
D. F.
Griffiths
and
G. A.
Watson
, eds.,
Longman
,
Harlow, UK
, pp.
36
56
.
8.
Corliss
,
G.
,
Griewank
,
A.
, and
Henneberger
,
P.
,
1997
, “
High-Order Stiff ODE Solvers via Automatic Differentiation and Rational Prediction
,”
Numerical Analysis and Its Applications
, Vol.
1196
(Lecture Notes in Computer Science),
L.
Vulkov
,
J.
Wasniewski
, and
P.
Yalamov
, eds.,
Springer
,
Berlin
, pp.
114
125
.
9.
,
E.
,
Nakhla
,
M.
, and
Achar
,
R.
,
2009
, “
A-Stable and L-Stable High-Order Integration Methods for Solving Stiff Differential Equations
,”
IEEE Trans. Comput.-Aided Des.
,
28
(
9
), pp.
1359
1372
10.
Zhou
,
Y.
,
,
E.
, and
Nakhla
,
M. S.
,
2012
, “
Structural Characterization and Efficient Implementation Techniques For-Stable High-Order Integration Methods
,”
IEEE Trans. Comput.-Aided Des.
,
31
(
1
), pp.
101
108
11.
Milenkovic
,
P.
,
2011
, “
Solution of the Forward Dynamics of a Single-Loop Linkage Using Power Series
,”
ASME J. Dyn. Syst., Meas., Control
,
133
(
6
), p.
061002
.10.1115/1.4004766
12.
Griewank
,
A.
, and
Walther
,
A.
,
2008
,
Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation
,
Society for Industrial and Applied Mathematics (SIAM)
,
.
13.
Griewank
,
A.
, and
Walther
,
A.
,
2006
, “
On the Efficient Generation of Taylor Expansions for DAE Solutions by Automatic Differentiation
,”
Applied Parallel Computing, State of the Art in Scientific Computing
,
J.
Dongarra
,
K.
, and
J.
Wasniewski
, eds.,
Springer
,
Berlin
, pp.
1089
1098
.
14.
Milenkovic
,
P.
,
2012
, “
Series Solution for Finite Displacement of Single-Loop Spatial Linkages
,”
ASME J. Mech. Rob.
,
4
(
2
), p.
021016
.10.1115/1.4006193
15.
Scholz
,
H. E.
,
1998
, “
The Examination of Nonlinear Stability and Solvability of the Algebraic Equations for the Implicit Taylor Series Method
,”
Appl. Numer. Math.
,
28
(
2
), pp.
439
458
.10.1016/S0168-9274(98)00059-2
16.
Wall
,
H. S.
,
1948
,
Analytic Theory of Continued Fractions
,
D. Van Nostrand
,
New York
.
17.
Barrio
,
R.
,
2005
, “
Performance of the Taylor Series Method for ODEs/DAEs
,”
Appl. Math. Comput.
,
163
(
2
), pp.
525
545
.10.1016/j.amc.2004.02.015
18.
Wilf
,
H. S.
,
1994
,
Generatingfunctionology
,
,
San Diego, CA
.
19.
Milenkovic
,
P.
,
2013
, “
Projective Constraint Stabilization for a Power Series Forward Dynamics Solver
,”
ASME J. Dyn. Syst., Meas., Control
,
135
(
3
), p.
031004
.10.1115/1.4023212