Skip to main content

Flexibility matrix method: portal frame

Short answer

Use the flexibility method with redundants at the hinge D. Compatibility gives X1 = 42.23 kN (left) and X2 = 104.55 kN (upward), so HD = 42.23 kN left and VD = 104.55 kN upward.

Analyse the portal frame ABCD shown in the figure below using the flexibility matrix method. Take EI=12000 kN/m2EI = 12000\ \text{kN/m}^2. [Figure: frame with column AB of height 4 m (2 m + 2 m) with 50 kN horizontal at mid-height, beam BC of 8 m carrying 25 kN/m UDL, column CD of height 3 m; A fixed, D supported.]

Reading of the figure: A (0, 0) fixed; column AB is 4 m high with a 50 kN horizontal load (towards the right) at mid-height E; beam BC is 8 m long with 25 kN/m; column CD is 3 m high, so D is 1 m above the level of A, and D is hinged. EI=12000EI = 12000 kNm2^2 is constant.

r=3+2=5r = 3 + 2 = 5, so Ds=2D_s = 2. Release the hinge at D and take the reactions at D as the redundants: X1X_1 = horizontal (→\rightarrow), X2X_2 = vertical (↑\uparrow). The primary structure is a cantilever fixed at A. In the flexibility method the compatibility condition is [F]{X}+{Δ0}=0[F]\{X\} + \{\Delta_0\} = 0.

Moments in the primary structure (sign: + tension on the inner face)

  • Loads only (M0M_0): at A, −(50×2+25×8×4)=−900-(50\times2 + 25\times8\times4) = -900; at E, −800-800; at B, −800-800; the beam falls parabolically to 0 at C.
  • Unit X1X_1 (m1m_1): m1=m_1 = lever arm of the unit horizontal force; at C and along BC it is 33, at B it is 33, at A it is −1-1.
  • Unit X2X_2 (m2m_2): at B it is 88, falling linearly to 0 at C; along AB it is 88.

Flexibility matrix and load vector

f11=90.333EIf12=128.000EIf22=426.667EIΔ10=−9566.7EIΔ20=−39200.0EI\begin{aligned} f_{11} &= \frac{90.333}{EI} \\ f_{12} &= \frac{128.000}{EI} \\ f_{22} &= \frac{426.667}{EI} \\ \Delta_{10} &= \frac{-9566.7}{EI} \\ \Delta_{20} &= \frac{-39200.0}{EI} \end{aligned}

In units of m/kN with EI=12000EI = 12000:

[F]=[0.007530.010670.010670.03556],{Δ0}={−0.79722−3.26667} m[F] = \begin{bmatrix} 0.00753 & 0.01067 \\ 0.01067 & 0.03556 \end{bmatrix},\qquad \{\Delta_0\} = \begin{Bmatrix} -0.79722 \\ -3.26667 \end{Bmatrix}\ \text{m}

Compatibility equations

90.333 X1+128.000 X2=9566.7128.000 X1+426.667 X2=39200.0\begin{aligned} 90.333\,X_1 + 128.000\,X_2 &= 9566.7 \\ 128.000\,X_1 + 426.667\,X_2 &= 39200.0 \end{aligned} {X}=−[F]−1{Δ0} ⇒ X1=−42.23 kN,X2=104.55 kN\{X\} = -[F]^{-1}\{\Delta_0\} \ \Rightarrow\ X_1 = -42.23\ \text{kN},\quad X_2 = 104.55\ \text{kN}

Final bending moments (kNm)

SectionM0M_0m1m_1m2m_2M=M0+∑XmM=M_0+\sum X m
A-900.00-1.008.00-21.41
E-800.001.008.00-5.87
B-800.003.008.00-90.34
mid BC-200.003.004.0091.48
C0.003.000.00-126.70
D0.000.000.000.00

Maximum sagging moment in BC: 91.89 kNm at 3.82 m from B.

Reactions

  • At D (hinge): HD=42.23H_D = 42.23 kN (←\leftarrow), VD=104.55V_D = 104.55 kN (↑\uparrow)
  • At A: HA=7.77H_A = 7.77 kN (←\leftarrow), VA=95.45V_A = 95.45 kN (↑\uparrow), MA=21.41M_A = 21.41 kNm (anticlockwise)

Check: ΣV=200.00=25×8=200\Sigma V = 200.00 = 25 \times 8 = 200 kN; ΣH=50−7.77−42.23=0\Sigma H = 50 - 7.77 - 42.23 = 0.

Answer: X1=HD=42.23X_1 = H_D = 42.23 kN towards the left and X2=VD=104.55X_2 = V_D = 104.55 kN upward.

All Theory of Structures II chapter 2 questions →