Guyan reduction

In computational mechanics, Guyan reduction, also known as static condensation, is a dimensionality reduction method which reduces the number of degrees of freedom by ignoring the inertial terms of the equilibrium equations and expressing the unloaded degrees of freedom in terms of the loaded degrees of freedom.

Basic concept

The static equilibrium equation can be expressed as:

$\mathbf {K} \mathbf {d} =\mathbf {f}$

where $\mathbf {K}$ is the stiffness matrix, $\mathbf {f}$ the force vector, and $\mathbf {d}$ the displacement vector. The number of the degrees of freedom of the static equilibrium problem is the length of the displacement vector. By partitioning the above system of linear equations with regards to loaded (master) and unloaded (slave) degrees of freedom, the static equilibrium equation may be expressed as:

${\begin{bmatrix}\mathbf {K} _{mm}&\mathbf {K} _{ms}\\\mathbf {K} _{sm}&\mathbf {K} _{ss}\end{bmatrix}}{\begin{Bmatrix}\mathbf {d} _{m}\\\mathbf {d} _{s}\end{Bmatrix}}={\begin{Bmatrix}\mathbf {f} _{m}\\\mathbf {f} _{s}\end{Bmatrix}}$

Focusing on the lower partition of the above system of linear equations, the dependent (slave) degrees of freedom are expressed by the following equation.

$\mathbf {K} _{sm}\mathbf {d} _{m}+\mathbf {K} _{ss}\mathbf {d} _{s}=\mathbf {f} _{s}$

Solving the above equation in terms of the independent (master) degrees of freedom leads to the following dependency relations

$\mathbf {d} _{s}=\mathbf {K} _{ss}^{-1}\mathbf {f} _{s}-\mathbf {K} _{ss}^{-1}\mathbf {K} _{sm}\mathbf {d} _{m}$

Substituting the dependency relations on the upper partition of the static equilibrium problem condenses away the slave degrees of freedom, leading to the following reduced system of linear equations.

$\left[\mathbf {K} _{mm}-\mathbf {K} _{ms}\mathbf {K} _{ss}^{-1}\mathbf {K} _{sm}\right]\mathbf {d} _{m}=\mathbf {f} _{m}-\mathbf {K} _{ms}\mathbf {K} _{ss}^{-1}\mathbf {f} _{s}$

This can be rewritten as:

$\mathbf {K} _{reduced}\mathbf {d} _{m}=\mathbf {f} _{reduced}$

The above system of linear equations is equivalent to the original problem, but expressed in terms of the master's degrees of freedom alone. Thus, the Guyan reduction method results in a reduced system by condensing away the slave degrees of freedom.

Linear transformation

The Guyan reduction can also be expressed as a change of basis which produces a low-dimensional representation of the original space, represented by the master's degrees of freedom. The linear transformation that maps the reduced space onto the full space is expressed as:

${\begin{Bmatrix}\mathbf {d} _{m}\\\mathbf {d} _{s}\end{Bmatrix}}={\begin{bmatrix}\mathbf {I} \\-\mathbf {K} _{ss}^{-1}\mathbf {K} _{sm}\end{bmatrix}}{\begin{Bmatrix}\mathbf {d} _{m}\end{Bmatrix}}={\begin{Bmatrix}\mathbf {T} _{G}\end{Bmatrix}}{\begin{Bmatrix}\mathbf {d} _{m}\end{Bmatrix}}$

where $\mathbf {T} _{G}$ represents the Guyan reduction transformation matrix. Thus, the reduced problem is represented as:

$\mathbf {K} _{G}\mathbf {d} _{m}=\mathbf {f} _{m}$

In the above equation, $\mathbf {K} _{G}$ represents the reduced system of linear equations that's obtained by applying the Guyan reduction transformation on the full system, which is expressed as:

$\mathbf {K} _{G}=\mathbf {T} _{G}^{T}\mathbf {K} \mathbf {T} _{G}$

See also

• Finite element method