Newmark-beta method

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,

$M{\ddot {u}}+C{\dot {u}}+f^{\textrm {int}}(u)=f^{\textrm {ext}}\,$

here $M$ is the mass matrix, $C$ is the damping matrix, $f^{\textrm {int}}$ and $f^{\textrm {ext}}$ are internal force per unit displacement and external forces, respectively.

Using the extended mean value theorem, the Newmark- $\beta$ 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 }\,

where

{\ddot {u}}_{\gamma }=(1-\gamma ){\ddot {u}}_{n}+\gamma {\ddot {u}}_{n+1}~~~~0\leq \gamma \leq 1

therefore

{\dot {u}}_{n+1}={\dot {u}}_{n}+(1-\gamma )\Delta t~{\ddot {u}}_{n}+\gamma \Delta t~{\ddot {u}}_{n+1}.

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 }

where again

{\ddot {u}}_{\beta }=(1-2\beta ){\ddot {u}}_{n}+2\beta {\ddot {u}}_{n+1}~~~~0\leq 2\beta \leq 1

The discretized structural equation becomes

${\begin{aligned}&{\dot {u}}_{n+1}={\dot {u}}_{n}+(1-\gamma )\Delta t~{\ddot {u}}_{n}+\gamma \Delta t~{\ddot {u}}_{n+1}\\&u_{n+1}=u_{n}+\Delta t~{\dot {u}}_{n}+{\frac {\Delta t^{2}}{2}}\left((1-2\beta ){\ddot {u}}_{n}+2\beta {\ddot {u}}_{n+1}\right)\\&M{\ddot {u}}_{n+1}+C{\dot {u}}_{n+1}+f^{\textrm {int}}(u_{n+1})=f_{n+1}^{\textrm {ext}}\,\end{aligned}}$

Explicit central difference scheme is obtained by setting $\gamma =0.5$ and $\beta =0$

Average constant acceleration (Middle point rule) is obtained by setting $\gamma =0.5$ and $\beta =0.25$

Stability Analysis

A time-integration scheme is said to be stable if there exists an integration time-step $\Delta t_{0}>0$ so that for any $\Delta t\in (0,\Delta t_{0}]$ , a finite variation of the state vector $q_{n}$ at time $t_{n}$ induces only a non-increasing variation of the state-vector $q_{n+1}$ calculated at a subsequent time $t_{n+1}$ . Assume the time-integration scheme is

$q_{n+1}=A(\Delta t)q_{n}+g_{n+1}(\Delta t)$

The linear stability is equivalent to $\rho (A(\Delta t))\leq 1$ , here $\rho (A(\Delta t))$ is the spectral radius of the update matrix $A(\Delta t)$ .

For the linear structural equation

$M{\ddot {u}}+C{\dot {u}}+Ku=f^{\textrm {ext}}\,$

here $K$ is the stiffness matrix. Let $q_{n}=[{\dot {u}}_{n},u_{n}]$ , the update matrix is $A=H_{1}^{-1}H_{0}$ , and

${\begin{aligned}H_{1}={\begin{bmatrix}M+\gamma \Delta tC&\gamma \Delta tK\\\beta \Delta t^{2}C&M+\beta \Delta t^{2}K\end{bmatrix}}\qquad H_{0}={\begin{bmatrix}M-(1-\gamma )\Delta tC&-(1-\gamma )\Delta tK\\-({\frac {1}{2}}-\beta )\Delta t^{2}C+\Delta tM&M-({\frac {1}{2}}-\beta )\Delta t^{2}K\end{bmatrix}}\end{aligned}}$

For undamped case ( $C=0$ ), the update matrix can be decoupled by introducing the eigenmodes $u=e^{i\omega _{i}t}x_{i}$ of the structural system, which are solved by the generalized eigenvalue problem

$\omega ^{2}Mx=Kx\,$

For each eigenmode, the update matrix becomes

${\begin{aligned}H_{1}={\begin{bmatrix}1&\gamma \Delta t\omega _{i}^{2}\\0&1+\beta \Delta t^{2}\omega _{i}^{2}\end{bmatrix}}\qquad H_{0}={\begin{bmatrix}1&-(1-\gamma )\Delta t\omega _{i}^{2}\\\Delta t&1-({\frac {1}{2}}-\beta )\Delta t^{2}\omega _{i}^{2}\end{bmatrix}}\end{aligned}}$

The characteristic equation of the update matrix is

$\lambda ^{2}-\left(2-(\gamma +{\frac {1}{2}})\eta _{i}^{2}\right)\lambda +1-(\gamma -{\frac {1}{2}})\eta _{i}^{2}=0\,\qquad \eta _{i}^{2}={\frac {\omega _{i}^{2}\Delta t^{2}}{1+\beta \omega _{i}^{2}\Delta t^{2}}}$

As for the stability, we have

Explicit central difference scheme ( $\gamma =0.5$ and $\beta =0$ ) is stable when $\omega \Delta t\leq 2$ .

Average constant acceleration (Middle point rule) ( $\gamma =0.5$ and $\beta =0.25$ ) is unconditionally stable.

References

Newmark, Nathan M. (1959), "A method of computation for structural dynamics", Journal of the Engineering Mechanics Division, 85 (EM3) (3): 67–94, doi:10.1061/JMCEA3.0000098

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[1] Newmark, Nathan M. (1959), "A method of computation for structural dynamics", Journal of the Engineering Mechanics Division, 85 (EM3) (3): 67–94, doi:10.1061/JMCEA3.0000098

Numerical methods for integration
First-order methods	Euler method Backward Euler Semi-implicit Euler Exponential Euler
Second-order methods	Verlet integration Velocity Verlet Trapezoidal rule Beeman's algorithm Midpoint method Heun's method Newmark-beta method Leapfrog integration
Higher-order methods	Exponential integrator Runge–Kutta methods List of Runge–Kutta methods Linear multistep method General linear methods Backward differentiation formula Yoshida Gauss–Legendre method
Theory	Symplectic integrator