Restricted Pad scheme in computational structural dynamics
Restricted Pad schemes where the denominator is a complete power p of a linear function have been examined analytically by Nrsett. Here this method is evaluated with paying special attention to numerical aspects when solving typical finite-element systems in structural dynamics with n degrees of freedom and with essential properties like sparsety and symmetry. This paper presents an efficient implicit restricted time-stepping scheme with a final coefficient matrix of order n and a leading error term proportional to h4. This is a significant improvement compared with the most popular Pad- -Newmark scheme with matrix order n, too, but error term h3 on the one hand and the well established Pad- -discontinuous Hamiltonian approach with error term h4 but matrix order 2n on the other hand. The classical -scheme and the new restricted -scheme are compared in a modal manner and by assessing the correlation between local error estimator and global error.
pad schemes computational structural dynamics error estimation time-finite elements time-step control modal assessment