Demystifying Matrix Methods of Structural Analysis: Flexibility vs. Stiffness

Matrix methods are structural engineering tools designed to solve complex frames, beams, and trusses efficiently by converting continuous structural behavior into linear algebraic systems of the form $[A]\{x\} = \{b\}$.

1. Fundamental Overview: Flexibility vs. Stiffness Method

Both formulation strategies solve structural systems by satisfying equilibrium, compatibility, and constitutive (force-displacement) relationships simultaneously.

Flexibility Method (Force / Compatibility Method)

• Primary Focus: Works by selecting redundant joint forces/moments as primary unknowns.
• Formulation Basis: Built on structural flexibility (compliance).
• Core Governing Equation: $\{\delta\} = [f]\{P\} + \{\delta^0\}$
• Variables: $[f]$ is the Flexibility Matrix, $\{P\}$ is the Nodal Load Vector, $\{\delta^0\}$ represents initial displacements (such as support settlements or thermal strain), and $\{\delta\}$ is the Total Displacement Vector.

Stiffness Method (Displacement / Equilibrium Method)

• Primary Focus: Uses unknown joint displacements (degrees of freedom) as primary variables.
• Formulation Basis: Built on structural stiffness.
• Core Governing Equation: $[K]\{\delta\} = \{P\} + \{P^0\}$
• Variables: $[K]$ is the System Stiffness Matrix, $\{\delta\}$ is the Nodal Displacement Vector, $\{P\}$ represents External Nodal Loads, and $\{P^0\}$ accounts for Equivalent Fixed-End Forces.

2. Head-to-Head Comparison Matrix

Aspect	Flexibility Method (Force Method)	Stiffness Method (Displacement Method)
Primary Unknowns	Redundant Forces / Moments ($R$)	Nodal Displacements / Rotations ($\delta, \theta$)
Governing Matrix	Flexibility Matrix $[f]$	Stiffness Matrix $[K]$
Matrix Size	Equal to Degree of Static Indeterminacy ($D_s$)	Equal to Degree of Kinematic Indeterminacy ($D_k$)
Symmetry	Symmetric ($f_{ij} = f_{ji}$)	Symmetric ($K_{ij} = K_{ji}$)
Computational Efficiency	Less efficient for large structures	Highly efficient; easily automated for software
Best Suited For	Small to medium, low-$D_s$ indeterminate systems	Large, complex, highly indeterminate structures
Application Domain	Theoretical analysis, hand calculations	Commercial Finite Element Analysis (FEA) software

Useful Trick: Deciding Which Method to Use

To pick the fastest route by hand, compare the degrees of indeterminacy:

• If $D_s < D_k \rightarrow$ Use the Flexibility Method (fewer redundant force unknowns).
• If $D_k < D_s \rightarrow$ Use the Stiffness Method (fewer joint displacement unknowns).
• For computer programming $\rightarrow$ Always use the Stiffness Method. Its boundary-condition assembly process is systematic and easily scales in computer code.

3. Fundamental Types of Displacements & Formulation Coefficients

Flexibility ($f$) is the displacement produced by a unit force ($f = \frac{\delta}{P}$), whereas Stiffness ($k$) is the force required to produce a unit displacement ($k = \frac{P}{\delta}$). Notice that $k = \frac{1}{f}$.

Deformation Type	Displacement Formula ($\delta$ or $\theta$)	Flexibility Coefficient ($f$)	Stiffness Coefficient ($k$)
1. Axial Displacement	$\delta = \frac{PL}{AE}$	$f = \frac{L}{AE}$	$k = \frac{AE}{L}$
2. Shear Displacement	$\delta = \frac{PL}{k_{s}GA}$	$f = \frac{L}{k_{s}GA}$	$k = \frac{k_{s}GA}{L}$
3. Bending (Cantilever End Load)	$\delta = \frac{PL^3}{3EI}$	$f = \frac{L^3}{3EI}$	$k = \frac{3EI}{L^3}$
Bending (Simply Supported Center Load)	$\delta = \frac{PL^3}{48EI}$	$f = \frac{L^3}{48EI}$	$k = \frac{48EI}{L^3}$
4. Torsional Displacement	$\theta = \frac{TL}{GJ}$	$f = \frac{L}{GJ}$	$k = \frac{GJ}{L}$

Note on 2D Plane Frame Joint Displacements: Every unconstrained node in a 2D frame possesses 3 degrees of freedom: horizontal translation ($u$), vertical translation ($v$), and rotation ($\theta$).

4. Worked Examples

Example A: Flexibility Method (Propped Cantilever Beam)

Problem: Analyze a propped cantilever beam of length $L$ carrying a uniform load $w$, with constant $EI$. Find the vertical reaction at the prop support $B$.

Step 1: Determine $D_s$ & Select Redundant
Fixed support $A$ (3 reactions) + Prop $B$ (1 reaction) = 4 reactions. Equations of equilibrium = 3.
$D_s = 4 - 3 = 1$. Select prop vertical reaction $R_B$ as the redundant force ($P_1$).

Step 2: Primary Released Structure
Remove prop $B$ to form a determinate cantilever beam subjected to uniform load $w$.
Downward displacement at end $B$ due to load $w$:
$$\delta_1^0 = \frac{wL^4}{8EI}$$

Step 3: Flexibility Coefficient
Apply a unit upward load $P_1 = 1$ at $B$. The upward deflection at $B$ is:
$$f_{11} = \frac{L^3}{3EI}$$

Step 4: Compatibility Equation & Solve
Total displacement at unyielding support $B$ must equal zero:
$$\delta_1 = \delta_1^0 + f_{11} R_B = 0$$ $$-\frac{wL^4}{8EI} + \left(\frac{L^3}{3EI}\right) R_B = 0 \implies R_B = \frac{3}{8}wL$$

Example B: Stiffness Method (2D Frame Node Assembly)

Problem: Consider a portal frame with fixed bases at Nodes 3 and 4, and rigid joints at Nodes 1 and 2 subjected to a point load $P$. Set up the system unknown vector and solve for displacements.

Step 1: Identify Nodal Degrees of Freedom ($D_k$)
• Node 1 (Unconstrained): Displacements $\{u_1, v_1, \theta_1\}$
• Node 2 (Unconstrained): Displacements $\{u_2, v_2, \theta_2\}$
• Nodes 3 & 4 (Fixed): $\{u_3=0, v_3=0, \theta_3=0\}$ and $\{u_4=0, v_4=0, \theta_4=0\}$
System Nodal Displacement Vector: $\{\delta\} = [u_1, v_1, \theta_1, u_2, v_2, \theta_2]^T$

Step 2: Assemble System Stiffness Matrix $[K]$ & Load Vector $\{P\}$
Form member stiffness matrices in global coordinates and map them into the global system matrix $[K]_{6 \times 6}$.
Equilibrium equation: $[K]\{\delta\} = \{P\}$

Step 3: Solve for Displacements and Member Actions
Solve the matrix system: $\{\delta\} = [K]^{-1}\{P\}$
Once joint displacements $\{\delta\}$ are determined, back-calculate member end actions (bending moments, shear forces, axial forces) and support reactions using individual element transformation matrices.

5. Key Takeaways for Students

• Matrix methods convert continuous structural physics into solvable linear algebraic equations.
• Stiffness method dominates industry software due to its standardized assembly and scalability for large structures.
• Accurate modeling of joint displacements is essential for safe structural design.

Yogendra Gopal Borse

Civil Engineer | Assistant Engineer Grade-I, Maharashtra PWD

B.Tech (Civil) from VJTI Mumbai. Experienced in bridge design, road works, estimation, project monitoring and digital engineering tools. Creator of YogiPWD – practical technical resources for civil engineers.

Read full About →