3.2.6.7. Hilber-Hughes-Taylor Method

integrator HHT $alpha <$gamma $beta>

This command is used to construct a Hilber-Hughes-Taylor (HHT) integration object. This is an implicit method that allows for energy dissipation and second order accuracy (which is not possible with the regular Newmark method). Depending on choices of input parameters, the method can be unconditionally stable.

Argument

Type

Description

$alpha

float

\(\alpha\) factor

$gamma

float

\(\gamma\) factor

$beta

float

\(\beta\) factor

Example

  1. Tcl Code

integrator HHT 0.9
  1. Python Code

integrator('HHT', 0.9)

Note

\(\alpha\) is defined differently that in the paper, we use \(\alpha = \alpha_{HHT} - 1\) where \(\alpha_{HHT}\) is that used in the paper.

  • Like Newmark and all the implicit schemes, the unconditional stability of this method applies to linear problems. There are no results showing stability of this method over the wide range of nonlinear problems that potentially exist. Experience indicates that the time step for implicit schemes in nonlinear situations can be much greater than those for explicit schemes.

  • \(\alpha = 1.0\) corresponds to the Newmark method.

  • \(\alpha\) should be between 0.67 and 1.0. The smaller the \(\alpha\) the greater the numerical damping.

  • \(\gamma\) and \(\beta\) are optional. The default values ensure the method is second order accurate and unconditionally stable when \(\alpha\) is \(\tfrac{2}{3} <= \alpha <= 1.0\). The defaults are: \(\beta = \frac{(2 - \alpha)^2}{4}\) and \(\gamma = \frac{3}{2} - \alpha\)

3.2.6.7.1. Theory

The HHT method (sometimes called the :math:` alpha` method) is a one step implicit method for solving the transient problem which attempts to increase the amount of numerical damping present without degrading the order of accuracy. In the HHT method, the same Newmark approximations are used:

\[U_{t+\Delta t} = U_t + \Delta t \dot U_t + [(0.5 - \beta) \Delta t^2] \ddot U_t + [\beta \Delta t^2] \ddot U_{t+\Delta t}\]
\[\dot U_{t+\Delta t} = \dot U_t + [(1-\gamma)\Delta t] \ddot U_t + [\gamma \Delta t ] \ddot U_{t+\Delta t}\]

but the time-discrete momentum equation is modified:

\[R_{t + \alpha \Delta t} = F_{t+\Delta t}^{ext} - M \ddot U_{t + \Delta t} - C \dot U_{t+\alpha \Delta t} - F^{int}(U_{t + \alpha \Delta t})\]

`

where the displacements and velocities at the intermediate point are given by:

\[U_{t+ \alpha \Delta t} = (1 - \alpha) U_t + \alpha U_{t + \Delta t}\]
\[\dot U_{t+\alpha \Delta t} = (1-\alpha) \dot U_t + \alpha \dot U_{t + \Delta t}\]

Following the methods outlined for Newmarks method, loinearization of the nonlinear momentum equation results in the following linear equations:

\[K_{t+\Delta t}^{*i} d U_{t+\Delta t}^{i+1} = R_{t+\Delta t}^i\]

where

\[K_{t+\Delta t}^{*i} = \alpha K_t + \frac{\alpha \gamma}{\beta \Delta t} C_t + \frac{1}{\beta \Delta t^2} M\]

and

\[R_{t+\Delta t}^i = F_{t + \Delta t}^{ext} - F(U_{t + \alpha \Delta t}^{i-1})^{int} - C \dot U_{t+\alpha \Delta t}^{i-1} - M \ddot U_{t+ \Delta t}^{i-1}\]

The linear equations are used to solve for \(U_{t+\alpha \Delta t}, \dot U_{t + \alpha \Delta t} \ddot U_{t+\Delta t}\). Once convergence has been achieved the displacements and velocities at time \(t + \Delta t\) can be computed.