Skip to main content

Stiffness matrix: frame, end moments

Short answer

The structure stiffness matrix is EI[2.1333 0.6667 -0.2400; 0.6667 2.9333 -0.4800; -0.2400 -0.4800 0.2880] for {θ_B, θ_C, Δ}. End moments: M_AB = -53.27, M_BA = -36.76, M_BC = 36.76, M_CB = 103.45, M_CD = -93.45, M_DC = -116.51 kNm.

Generate the stiffness matrix of the frame loaded as shown in figure below. Also determine the end moments considering stiffness equations of each member. [Figure: portal frame ABCD; beam BC (2I) with 50 kN at mid-span (3 m + 3 m) and 10 kNm moment at C; column AB (I) 5 m, column CD (2I) 5 m; 60 kN horizontal at B; A and D fixed.]

Assumptions: columns AB = 5 m (I) and CD = 5 m (2I), girder BC = 6 m (2I); A and D fixed. The 60 kN horizontal load at B acts to the right, the 50 kN load acts downward at mid-span of BC and the 10 kNm moment at C is clockwise. Axial deformation is neglected.

Coordinates

1 = θB\theta_B, 2 = θC\theta_C (clockwise +), 3 = Δ\Delta (sway to the right). ψAB=ψCD=Δ/5\psi_{AB}=\psi_{CD}=\Delta/5.

Fixed-end moments

FEMBC=−PL8=−50×68=−37.5FEM_{BC}=-\dfrac{PL}{8}=-\dfrac{50\times6}{8}=-37.5 kNm, FEMCB=+37.5FEM_{CB}=+37.5 kNm. Other members unloaded.

Stiffness equations of each member

  • MAB=2EI5(θB−(0.6000)Δ1)M_{AB}=\frac{2EI}{5}\left(\theta_{B} - (0.6000)\Delta_1\right)
  • MBA=2EI5(2θB−(0.6000)Δ1)M_{BA}=\frac{2EI}{5}\left(2\theta_{B} - (0.6000)\Delta_1\right)
  • MBC=2(2EI)6(2θB+θC)−37.50M_{BC}=\frac{2(2EI)}{6}\left(2\theta_{B} + \theta_{C}\right) - 37.50
  • MCB=2(2EI)6(2θC+θB)+37.50M_{CB}=\frac{2(2EI)}{6}\left(2\theta_{C} + \theta_{B}\right) + 37.50
  • MCD=2(2EI)5(2θC−(0.6000)Δ1)M_{CD}=\frac{2(2EI)}{5}\left(2\theta_{C} - (0.6000)\Delta_1\right)
  • MDC=2(2EI)5(θC−(0.6000)Δ1)M_{DC}=\frac{2(2EI)}{5}\left(\theta_{C} - (0.6000)\Delta_1\right)

Stiffness matrix of the frame

Joint B: K11=4EI5+4(2EI)6K_{11}=\dfrac{4EI}{5}+\dfrac{4(2EI)}{6}, joint C: K22=4(2EI)6+4(2EI)5K_{22}=\dfrac{4(2EI)}{6}+\dfrac{4(2EI)}{5}, coupling K12=2(2EI)6K_{12}=\dfrac{2(2EI)}{6}, sway terms K13=−6EI52K_{13}=-\dfrac{6EI}{5^2}, K23=−6(2EI)52K_{23}=-\dfrac{6(2EI)}{5^2}, K33=12EI53+12(2EI)53K_{33}=\dfrac{12EI}{5^3}+\dfrac{12(2EI)}{5^3}.

[K]=EI[2.13330.6667−0.24000.66672.9333−0.4800−0.2400−0.48000.2880][K]=EI\begin{bmatrix}2.1333 & 0.6667 & -0.2400 \\ 0.6667 & 2.9333 & -0.4800 \\ -0.2400 & -0.4800 & 0.2880\end{bmatrix}

Load vector (applied loads minus fixed-end effects): {F}={37.500, −27.500, 60.000}T\{F\}=\{37.500,\ -27.500,\ 60.000\}^T (rows 1 and 2: −∑FEM-\sum FEM plus the applied moment at B and C; row 3: the 60 kN lateral force).

EI[2.13330.6667−0.24000.66672.9333−0.4800−0.2400−0.48000.2880]{θBθCΔ1}={37.500−27.50060.000}EI\begin{bmatrix}2.1333 & 0.6667 & -0.2400 \\ 0.6667 & 2.9333 & -0.4800 \\ -0.2400 & -0.4800 & 0.2880\end{bmatrix}\begin{Bmatrix}\theta_{B} \\ \theta_{C} \\ \Delta_{1}\end{Bmatrix}=\begin{Bmatrix}37.500 \\ -27.500 \\ 60.000\end{Bmatrix} EIθB=41.2829,EIθC=28.8240,EIΔ1=290.7758EI\theta_{B}=41.2829,\quad EI\theta_{C}=28.8240,\quad EI\Delta_{1}=290.7758

End moments using the member stiffness equations (kNm, clockwise +)

EndMM (kNm)EndMM (kNm)
AB-53.27BA-36.76
BC36.76CB103.45
CD-93.45DC-116.51

Checks:

  • Joint B: (-36.76) + (36.76) = 0.00 kNm (applied clockwise moment 0.00) ✓
  • Joint C: (103.45) + (-93.45) = 10.00 kNm (applied clockwise moment 10.00) ✓
  • Storey shear: −90.03+(−209.97)5+60=0.00\dfrac{-90.03+(-209.97)}{5}+60=0.00 ✓
MemberBM at first endBM at second endMax within member
AB-53.2736.76-
BC36.76-103.4541.65 at 3.00 m
CD-93.45116.51-

Reactions: A: H=−18.01H=-18.01, V=1.63V=1.63, M=53.27M=53.27; D: H=−41.99H=-41.99, V=48.37V=48.37, M=116.51M=116.51.

Answer: MAB=−53.27M_{AB}=-53.27, MBA=−36.76M_{BA}=-36.76, MBC=36.76M_{BC}=36.76, MCB=103.45M_{CB}=103.45, MCD=−93.45M_{CD}=-93.45, MDC=−116.51M_{DC}=-116.51 kNm.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Theory of Structures II chapter 3 questions →