# Flexibility method

In structural engineering, the flexibility method, also called the method of consistent deformations, is the traditional method for computing member forces and displacements in structural systems. Its modern version formulated in terms of the members' flexibility matrices also has the name the matrix force method due to its use of member forces as the primary unknowns.[1]

## Member flexibility

Flexibility is the inverse of stiffness. For example, consider a spring that has Q and q as, respectively, its force and deformation:

• The spring stiffness relation is Q = k q where k is the spring stiffness.
• Its flexibility relation is q = f Q, where f is the spring flexibility.
• Hence, f = 1/k.

A typical member flexibility relation has the following general form:

${\displaystyle \mathbf {q} ^{m}=\mathbf {f} ^{m}\mathbf {Q} ^{m}+\mathbf {q} ^{om}\qquad \qquad \qquad \mathrm {(1)} }$

where

m = member number m.
${\displaystyle \mathbf {q} ^{m}}$ = vector of member's characteristic deformations.
${\displaystyle \mathbf {f} ^{m}}$ = member flexibility matrix which characterises the member's susceptibility to deform under forces.
${\displaystyle \mathbf {Q} ^{m}}$ = vector of member's independent characteristic forces, which are unknown internal forces. These independent forces give rise to all member-end forces by member equilibrium.
${\displaystyle \mathbf {q} ^{om}}$ = vector of member's characteristic deformations caused by external effects (such as known forces and temperature changes) applied to the isolated, disconnected member (i.e. with ${\displaystyle \mathbf {Q} ^{m}=0}$).

For a system composed of many members interconnected at points called nodes, the members' flexibility relations can be put together into a single matrix equation, dropping the superscript m:

${\displaystyle \mathbf {q} _{M\times 1}=\mathbf {f} _{M\times M}\mathbf {Q} _{M\times 1}+\mathbf {q} _{M\times 1}^{o}\qquad \qquad \qquad \mathrm {(2)} }$

where M is the total number of members' characteristic deformations or forces in the system.

Unlike the matrix stiffness method, where the members' stiffness relations can be readily integrated via nodal equilibrium and compatibility conditions, the present flexibility form of equation (2) poses serious difficulty. With member forces ${\displaystyle \mathbf {Q} _{M\times 1}}$ as the primary unknowns, the number of nodal equilibrium equations is insufficient for solution, in general—unless the system is statically determinate.

## Nodal equilibrium equations

To resolve this difficulty, first we make use of the nodal equilibrium equations in order to reduce the number of independent unknown member forces. The nodal equilibrium equation for the system has the form:

${\displaystyle \mathbf {R} _{N\times 1}=\mathbf {b} _{N\times M}\mathbf {Q} _{M\times 1}+\mathbf {W} _{N\times 1}\qquad \qquad \qquad \mathrm {(3)} }$

where

${\displaystyle \mathbf {R} _{N\times 1}}$: Vector of nodal forces at all N degrees of freedom of the system.
${\displaystyle \mathbf {b} _{N\times M}}$: The resulting nodal equilibrium matrix
${\displaystyle \mathbf {W} _{N\times 1}}$: The vector of forces arising from loading on the members.

In the case of determinate systems, matrix b is square and the solution for Q can be found immediately from (3) provided that the system is stable.

## The primary system

For statically indeterminate systems, M > N, and hence, we can augment (3) with I = M-N equations of the form:

${\displaystyle X_{i}=\alpha Q_{j}+\beta Q_{k}+...\qquad i=1,2,...I\qquad \qquad \mathrm {(4)} }$

The vector X is the so-called vector of redundant forces and I is the degree of statical indeterminacy of the system. We usually choose j, k, ..., ${\displaystyle \alpha }$, and ${\displaystyle \beta }$ such that ${\displaystyle X_{i}}$ is a support reaction or an internal member-end force. With suitable choices of redundant forces, the equation system (3) augmented by (4) can now be solved to obtain:

${\displaystyle \mathbf {Q} _{M\times 1}=\mathbf {B} _{R}\mathbf {R} _{N\times 1}+\mathbf {B} _{X}\mathbf {X} _{I\times 1}+\mathbf {Q} _{v\cdot M\times 1}\qquad \qquad \qquad \mathrm {(5)} }$

Substitution into (2) gives:

${\displaystyle \mathbf {q} _{M\times 1}=\mathbf {f} _{M\times M}{\Big (}\mathbf {B} _{R}\mathbf {R} _{N\times 1}+\mathbf {B} _{X}\mathbf {X} _{I\times 1}+\mathbf {Q} _{v\cdot M\times 1}{\Big )}+\mathbf {q} _{M\times 1}^{o}\qquad \qquad \qquad \mathrm {(6)} }$

Equations (5) and (6) are the solution for the primary system which is the original system that has been rendered statically determinate by cuts that expose the redundant forces ${\displaystyle \mathbf {X} }$. Equation (5) effectively reduces the set of unknown forces to ${\displaystyle \mathbf {X} }$.

## Compatibility equation and solution

Next, we need to set up ${\displaystyle I}$ compatibility equations in order to find ${\displaystyle \mathbf {X} }$. The compatibility equations restore the required continuity at the cut sections by setting the relative displacements ${\displaystyle \mathbf {r} _{X}}$ at the redundants X to zero. That is, using the unit dummy force method:

${\displaystyle \mathbf {r} _{X}=\mathbf {B} _{X}^{T}\mathbf {q} =\mathbf {B} _{X}^{T}{\Big [}\mathbf {f} {\Big (}\mathbf {B} _{R}\mathbf {R} +\mathbf {B} _{X}\mathbf {X} +\mathbf {Q} _{v}{\Big )}+\mathbf {q} ^{o}{\Big ]}=0\qquad \qquad \qquad \mathrm {(7a)} }$
or ${\displaystyle \mathbf {r} _{X}=\mathbf {F} _{XX}\mathbf {X} +\mathbf {r} _{X}^{o}=0\qquad \qquad \qquad \mathrm {(7b)} }$

where

${\displaystyle \mathbf {F} _{XX}=\mathbf {B} _{X}^{T}\mathbf {f} \mathbf {B} _{X}}$
${\displaystyle \mathbf {r} _{X}^{o}=\mathbf {B} _{X}^{T}{\Big [}\mathbf {f} {\Big (}\mathbf {B} _{R}\mathbf {R} +\mathbf {Q} _{v}{\Big )}+\mathbf {q} ^{o}{\Big ]}}$

Equation (7b) can be solved for X, and the member forces are next found from (5) while the nodal displacements can be found by

${\displaystyle \mathbf {r} _{R}=\mathbf {B} _{R}^{T}\mathbf {q} =\mathbf {F} _{RR}\mathbf {R} +\mathbf {r} _{R}^{o}}$

where

${\displaystyle \mathbf {F} _{RR}=\mathbf {B} _{R}^{T}\mathbf {f} \mathbf {B} _{R}}$ is the system flexibility matrix.
${\displaystyle \mathbf {r} _{R}^{o}=\mathbf {B} _{R}^{T}{\Big [}\mathbf {f} {\Big (}\mathbf {B} _{X}\mathbf {X} +\mathbf {Q} _{v}{\Big )}+\mathbf {q} ^{o}{\Big ]}}$

Supports' movements taking place at the redundants can be included in the right-hand-side of equation (7), while supports' movements at other places must be included in ${\displaystyle \mathbf {r} _{X}^{o}}$ and ${\displaystyle \mathbf {r} _{R}^{o}}$ as well.