Numerical Laplace inversion in problems of elastodynamics: Comparison of four algorithms

0480700 - ÚT 2018 RIV GB eng J - Journal Article
Adámek, V. - Valeš, František - Červ, Jan
Numerical Laplace inversion in problems of elastodynamics: Comparison of four algorithms.
Advances in Engineering Software. Roč. 113, November (2017), s. 120-129. ISSN 0965-9978. E-ISSN 1873-5339
R&D Projects: GA ČR(CZ) GAP101/12/2315
Institutional support: RVO:61388998
Keywords : inverse Laplace transform * numerical algorithm * wave propagation * multi-precision computation * Maple code
OECD category: Acoustics
Impact factor: 3.198, year: 2017

The objective of this work is to find a suitable algorithm for numerical Laplace inversion which could be used for effective and precise solution of elastodynamic problems. For this purpose, the capabilities of four algorithms are studied using three transforms resulted from analytical solutions of longitudinal waves in a thin rod, flexural waves in a thin beam and plane waves in a strip. In particular, the Gaver–Stehfest algorithm, the Gaver–Wynn’s rho algorithm, the Fixed-Talbot algorithm and the FFT algorithm combined with Wynn’s epsilon accelerator are tested. The codes written in Maple 16 employing multi- precision computations are presented for each method. Given the results obtained, the last mentioned algorithm proves to be the best. It is most efficient and it gives results of reasonable accuracy nearly for all tested times ranging from 3 ×10 −7 s to 3 ×10 3 s .
Permanent Link: http://hdl.handle.net/11104/0276813