Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Civil-Comp Proceedings
ISSN 1759-3433
CCP: 65
Edited by: B.H.V. Topping
Paper X.6

Automatic Time Stepping Algorithms for Implicit Numerical Simulations of Non-Linear Dynamics

L. Noels, L. Stainier, J.P. Ponthot and J. Bonini

Laboratoire de Techniques Aéronautiques et Spatiales, Université de Liège, Liège, Belgium

Full Bibliographic Reference for this paper
L. Noels, L. Stainier, J.P. Ponthot, J. Bonini, "Automatic Time Stepping Algorithms for Implicit Numerical Simulations of Non-Linear Dynamics", in B.H.V. Topping, (Editor), "Finite Elements: Techniques and Developments", Civil-Comp Press, Edinburgh, UK, pp 263-272, 2000. doi:10.4203/ccp.65.10.6
Variable step strategies are specially well suited to deal with problems characterized by high non-linearity and contact/impact. Both phenomena are typical of dynamic simulations of the interactions between a turbine blade and its casing, the most dramatic example being blade loss. Constant step size strategies do not give satisfactory answer for this kind of problems, since it is very difficult, if not impossible, for the user to find an appropriate time step that does not lead to divergence nor generate extremely costly computations. An automatic time stepping algorithm is proposed, which takes into account the recent history of accelerations in the bodies under consideration. More precisely, the adaptation algorithm is based on estimators of the integration error of the differential dynamic balance equations. This allows for adaptation of the time step to capture correctly the transient phenomena, with characteristic times which can range from relatively long (in regime) to very short (blade loss), thus ensuring precision while keeping the computation cost to a minimum. Additionally, the proposed algorithm automatically takes decisions regarding the necessity of updating the tangent matrix or stopping the iterations, further reducing the computational cost. As an illustration of the capabilities of this algorithm, several numerical simulations of both academic and industrial (the contact/impact between a turbine blade and the casing) problems will be presented.

purchase the full-text of this paper (price £20)

go to the previous paper
go to the next paper
return to the table of contents
return to the book description
purchase this book (price £80 +P&P)