Time integration of non-linear dynamic equations by means of a direct variational method

Abstract

Non-linear dynamic problems governed by ordinary (ODE) or partial differential equations (PDE) are herein approached by means of an alternative methodology. A generalized solution named WEM by the authors and previously developed for boundary value problems, is applied to linear and non-linear equations. A simple transformation after selecting an arbitrary interval of interest T allows using WEM in initial conditions problems and others with both initial and boundary conditions. When dealing with the time variable, the methodology may be seen as a time integration scheme. The application of WEM leads to arbitrary precision results. It is shown that it lacks neither numerical damping nor chaos which were found to be present with the application of some of the time integration schemes most commonly used in standard finite element codes (e.g., methods of central difference, Newmark, Wilson-θ, and so on). Illustrations include the solution of two non-linear ODEs which govern the dynamics of a single-degree-of-freedom (s.d.o.f.) system of a mass and a spring with two different non-linear laws (cubic and hyperbolic tangent respectively). The third example is the application of WEM to the dynamic problem of a beam with non-linear springs at its ends and subjected to a dynamic load varying both in space and time, even with discontinuities, governed by a PDE. To handle systems of non-linear equations iterative algorithms are employed. The convergence of the iteration is achieved by taking n partitions of T. However, the values of T/n are, in general, several times larger than the usual Δt in other time integration techniques. The maximum error (measured as a percentage of the energy) is calculated for the first example and it is shown that WEM yields an acceptable level of errors even when rather large time steps are used.