The Newmark-beta method is a method of numerical integration used to solve certain differential equations. It is widely used in numerical evaluation of the dynamic response of structures and solids such as in finite element analysis to model dynamic systems. The method is named after Nathan M. Newmark,[1] former Professor of Civil Engineering at the University of Illinois at Urbana–Champaign, who developed it in 1959 for use in structural dynamics. The semi-discretized structural equation is a second order ordinary differential equation system,
here
is the mass matrix,
is the damping matrix,
and
are internal and external forces.
Using the extended mean value theorem, the Newmark-
method states that the first time derivative (velocity in the equation of motion) can be solved as,
![\dot{u}_{n+1}=\dot{u}_n+ \Delta t~\ddot{u}_\gamma \,](https://wikimedia.org/api/rest_v1/media/math/render/svg/929d0f404ef505c8258cdc5f1ca78473aed95038)
where
![\ddot{u}_\gamma = (1 - \gamma)\ddot{u}_n + \gamma \ddot{u}_{n+1}~~~~0\leq \gamma \leq 1](https://wikimedia.org/api/rest_v1/media/math/render/svg/096f403fea2353b20bea317ee10a0481452794e6)
therefore
![\dot{u}_{n+1}=\dot{u}_n + (1 - \gamma) \Delta t~\ddot{u}_n + \gamma \Delta t~\ddot{u}_{n+1}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/0835ab6572026430bceeacdfa79aac4b81e049d4)
Because acceleration also varies with time, however, the extended mean value theorem must also be extended to the second time derivative to obtain the correct displacement. Thus,
![u_{n+1}=u_n + \Delta t~\dot{u}_n+\begin{matrix} \frac 1 2 \end{matrix} \Delta t^2~\ddot{u}_\beta](https://wikimedia.org/api/rest_v1/media/math/render/svg/a3971502feffa6bf78faadac970daef0a2bbf7a3)
where again
![\ddot{u}_\beta = (1 - 2\beta)\ddot{u}_n + 2\beta\ddot{u}_{n+1}~~~~0\leq 2\beta\leq 1](https://wikimedia.org/api/rest_v1/media/math/render/svg/7715190cae1f94bfb9aacec77e2957a74f747df3)
The discretized structural equation becomes
Explicit central difference scheme is obtained by setting
and
Average constant acceleration (Middle point rule) is obtained by setting
and
Stability Analysis[]
A time-integration scheme is said to be stable if there exists an integration time-step
so that for any
, a finite variation of the state vector
at time
induces only a non-increasing variation of the state-vector
calculated at a subsequent time
. Assume the time-integration scheme is
The linear stability is equivalent to
, here
is the spectral radius of the update matrix
.
For the linear structural equation
here
is the stiffness matrix. Let
, the update matrix is
, and
For undamped case (
), the update matrix can be decoupled by introducing the eigenmodes
of the structural system, which are solved by the generalized eigenvalue problem
For each eigenmode, the update matrix becomes
The characteristic equation of the update matrix is
As for the stability, we have
Explicit central difference scheme (
and
) is stable when
.
Average constant acceleration (Middle point rule) (
and
) is unconditionally stable.
References[]
- ^ Newmark, Nathan M. (1959), "A method of computation for structural dynamics", Journal of the Engineering Mechanics Division, 85 (EM3): 67–94
|
---|
First-order methods | |
---|
Second-order methods | |
---|
Higher-order methods | |
---|
Theory | |
---|