Skip to main content

Stiffness matrix: frame with coordinates 1, 2, 3

Short answer

Coordinates are theta_B, sway, theta_C. K=EI[[2,-0.1667,0.6667],[-0.1667,0.2431,-0.375],[0.6667,-0.375,2.3333]]. Solving with equivalent loads gives M_AB=-9.12, M_BA=8.13, M_BC=-3.13, M_CB=14.00, M_CD=-10.00, M_DC=-9.35 kNm.

Generate the stiffness matrix for the frame shown in figure below with respect to coordinates 1, 2 and 3 and use it to analyse the frame if the forces 5 kNm and 4 kN are acting at coordinates 1 and 3 respectively in addition to the external loads as shown in figure. Take EI is constant for all members. [Figure: portal frame; left column 6 m (3 m + 3 m) with 10 kN horizontal at mid-height, fixed base; beam of 3 m with 15 kN/m UDL; right column 4 m, fixed base; coordinates 1 (rotation), 2 (horizontal), 3 (rotation) marked at the top joints.]

Assumptions (figure partly illegible): A(0,0), B(0,6), C(3,6), D(3,2). The left column AB (6 m, 3 m + 3 m) carries the 10 kN horizontal load (to the right) at mid-height; the beam BC is 3 m long with 15 kN/m downward; the right column CD is 4 m. A and D are fixed. EI is the same for all members. Coordinates: 1 = rotation of joint B (clockwise), 2 = horizontal translation (sway) of the beam (to the right), 3 = rotation of joint C (clockwise). The additional forces at coordinates 1 and 3 are taken as a 5 kNm and a 4 kNm clockwise moment (a 'kN' at coordinate 3 cannot act on a rotation, so the 4 is read as kNm).

Stiffness matrix [K][K] (units EI)

Terms: ψAB=Δ/6\psi_{AB}=\Delta/6, ψCD=Δ/4\psi_{CD}=\Delta/4.

  • k11=4EI6+4EI3=2.0000EIk_{11}=\dfrac{4EI}{6}+\dfrac{4EI}{3}=2.0000EI; k13=k31=2EI3=0.6667EIk_{13}=k_{31}=\dfrac{2EI}{3}=0.6667EI; k33=4EI3+4EI4=2.3333EIk_{33}=\dfrac{4EI}{3}+\dfrac{4EI}{4}=2.3333EI
  • k12=−6EI62=−0.1667EIk_{12}=-\dfrac{6EI}{6^2}=-0.1667EI; k32=−6EI42=−0.3750EIk_{32}=-\dfrac{6EI}{4^2}=-0.3750EI; k22=12EI63+12EI43=0.2431EIk_{22}=\dfrac{12EI}{6^3}+\dfrac{12EI}{4^3}=0.2431EI
[K]=EI[2.0000−0.16670.6667−0.16670.2431−0.37500.6667−0.37502.3333][K]=EI\begin{bmatrix}2.0000 & -0.1667 & 0.6667 \\ -0.1667 & 0.2431 & -0.3750 \\ 0.6667 & -0.3750 & 2.3333\end{bmatrix}

Equivalent loads from the member loads (fixed-end effects)

Fixed-end moments: BC ∓15×3212=∓11.25\mp\dfrac{15\times3^2}{12}=\mp11.25 kNm; AB (10 kN at 3 m, horizontal) ∓PL8=∓7.5\mp\dfrac{PL}{8}=\mp7.5 kNm (sign for a horizontal load to the right on a column); the equivalent lateral force at the girder level is the fixed-end shear at B, 10/2=510/2=5 kN.

{F0}={3.7505.000−11.250} (kNm, kN, kNm),{F}={F0}+{504}={8.7505.000−7.250}\{F_0\}=\begin{Bmatrix}3.750\\5.000\\-11.250\end{Bmatrix}\ (\text{kNm, kN, kNm}),\qquad \{F\}=\{F_0\}+\begin{Bmatrix}5\\0\\4\end{Bmatrix}=\begin{Bmatrix}8.750\\5.000\\-7.250\end{Bmatrix}

Analysis

[K]{d}={F}[K]\{d\}=\{F\} gives

EIθB=6.7426,EIΔ=23.1753,EIθC=−1.3090EI\theta_B=6.7426,\quad EI\Delta=23.1753,\quad EI\theta_C=-1.3090

End moments (member equations, kNm, clockwise +)

  • MAB=2EI6(θB−(0.5000)Δ1)−7.50M_{AB}=\frac{2EI}{6}\left(\theta_{B} - (0.5000)\Delta_1\right) - 7.50
  • MBA=2EI6(2θB−(0.5000)Δ1)+7.50M_{BA}=\frac{2EI}{6}\left(2\theta_{B} - (0.5000)\Delta_1\right) + 7.50
  • MBC=2EI3(2θB+θC)−11.25M_{BC}=\frac{2EI}{3}\left(2\theta_{B} + \theta_{C}\right) - 11.25
  • MCB=2EI3(2θC+θB)+11.25M_{CB}=\frac{2EI}{3}\left(2\theta_{C} + \theta_{B}\right) + 11.25
  • MCD=2EI4(2θC−(0.7500)Δ1)M_{CD}=\frac{2EI}{4}\left(2\theta_{C} - (0.7500)\Delta_1\right)
  • MDC=2EI4(θC−(0.7500)Δ1)M_{DC}=\frac{2EI}{4}\left(\theta_{C} - (0.7500)\Delta_1\right)
EndMM (kNm)EndMM (kNm)
AB-9.12BA8.13
BC-3.13CB14.00
CD-10.00DC-9.35

Checks:

  • Joint B: (8.13) + (-3.13) = 5.00 kNm (applied clockwise moment 5.00) ✓
  • Joint C: (14.00) + (-10.00) = 4.00 kNm (applied clockwise moment 4.00) ✓

Bending moment ordinates (kNm)

MemberBM at first endBM at second endMax within member
AB-9.12-8.136.38 at 3.00 m
BC-3.13-14.008.75 at 1.26 m
CD-10.009.35-

Reactions: A: H=−5.16H=-5.16, V=18.88V=18.88, M=9.12M=9.12; D: H=−4.84H=-4.84, V=26.12V=26.12, M=9.35M=9.35.

Answer: [K][K] as above; with the loads given, MAB=−9.12M_{AB}=-9.12, MBA=8.13M_{BA}=8.13, MBC=−3.13M_{BC}=-3.13, MCB=14.00M_{CB}=14.00, MCD=−10.00M_{CD}=-10.00, MDC=−9.35M_{DC}=-9.35 kNm.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Theory of Structures II chapter 3 questions →