Matrix Methods of Structural Analysis

Why it matters

Matrix Methods of Structural Analysis are crucial for efficiently analyzing complex structures, especially indeterminate ones. These methods allow engineers to handle large systems of equations, making them indispensable in modern structural engineering for designing safe and economical structures.

Key ideas

• Matrix Representation: Structures can be represented using matrices, which simplifies the process of solving equilibrium equations.
• Stiffness Method: This method involves formulating the global stiffness matrix of the structure, which relates nodal displacements to applied forces.
• Flexibility Method: This method uses the flexibility matrix, which relates forces to displacements, and commonly uses redundant forces as unknowns with compatibility equations for indeterminate structures.
• Degrees of Freedom (DOF): The number of independent displacements or rotations that define the system's configuration.
• Global vs Local Coordinates: Transformations between local and global coordinate systems are necessary for assembling the global stiffness matrix.

Formulas

• F = K·d
  • F: Force vector (N)
  • K: Stiffness matrix (N/m)
  • d: Displacement vector (m)
• d = K⁻¹·F
  • d: Displacement vector (m)
  • K⁻¹: Inverse of stiffness matrix (m/N)
  • F: Force vector (N)

Worked example

A simply supported 6 m beam carries a 10 kN downward central load. E = 200 GPa and I = 400×10⁶ mm⁴ = 4×10⁻⁴ m⁴, so EI = 80×10⁶ N·m². Divide it into two 3 m Euler–Bernoulli elements with a node at midspan. End translations are zero; rotations are free.

Use downward displacement v and rotation θ = dv/dx. The free DOFs are [θ_A, v_C, θ_B]. Assembly and support constraints give

(EI/3³) [[36,-18,0],[-18,24,18],[0,18,36]] [θ_A,v_C,θ_B]ᵀ = [0,10000,0]ᵀ.

First and third equations give θ_A = v_C/2 and θ_B = -v_C/2, with the length units carried by matrix coefficients. Substitution gives v_C = 10000/(6EI/27) = 0.0005625 m. Thus θ_A = 0.00028125 rad and θ_B = -0.00028125 rad. The center rotation is zero by symmetry.

Answer: Midspan deflection is 0.5625 mm downward and reactions are 5 kN upward each. The independent check PL³/(48EI) gives the same displacement.

Mixed translation/rotation DOFs give different units in different stiffness entries. Apply constraints before solving the reduced system; the unconstrained global stiffness matrix contains rigid-body modes. Solve linear equations numerically rather than explicitly computing a matrix inverse.

Common mistakes

• Incorrectly assembling the global stiffness matrix due to errors in coordinate transformations.
• Neglecting to apply boundary conditions correctly, leading to incorrect displacement or force calculations.
• Misinterpreting the units, especially when converting between mm and m.

For GATE CE

Questions often involve calculating displacements or forces using the stiffness method. Practice problems on assembling stiffness matrices and applying boundary conditions are essential.

Quick check

1. What is the primary advantage of using matrix methods in structural analysis?
2. How does the stiffness method differ from the flexibility method?
3. What is the significance of degrees of freedom in matrix methods?

Answers: 1. Efficiently solving large systems of equations. 2. Stiffness method uses stiffness matrix; flexibility method uses flexibility matrix. 3. They define the independent displacements or rotations in the system.