Skip to main content

Stiffness matrix: frame, 25 kN, 100 kN

Short answer

Using the stiffness matrix method with joint rotations and sway as unknowns, the equilibrium equations are solved to give the final end moments: M_AB = -122.33 kNm, M_BA = 110.27 kNm, M_BC = -110.27 kNm, M_CB = 174.70 kNm, M_CD = -174.70 kNm, M_DC = -169.28 kNm.

Analyse the frame shown in figure below by using the stiffness matrix method and draw bending moment diagram. [Figure: frame; left column 2I (3 m + 7 m) with 100 kN horizontal; 25 kN vertical at top-left; beam (2I) of 8 m with 30 kN/m UDL; right column I of 5 m; fixed bases.]

Assumptions: A(0,0), B(0,10), C(8,10), D(8,5). Left column AB (2I) is 10 m high (3 m + 7 m) with the 100 kN horizontal force (to the right) acting 7 m above A; beam BC (2I) is 8 m and carries 30 kN/m (the '30 kNm' of the figure is read as a UDL of 30 kN/m); right column CD (I) is 5 m; A and D are fixed. The 25 kN vertical load at B only produces axial force in AB (it acts through the column axis) and no bending moment. Axial deformation is neglected.

Unknowns

θB,θC\theta_B,\theta_C and the sway Δ\Delta: ψAB=Δ/10\psi_{AB}=\Delta/10, ψCD=Δ/5\psi_{CD}=\Delta/5.

Fixed-end moments (clockwise +)

  • AB (100 kN at a=7a=7 m from A, b=3b=3 m): FEMAB=−63.00FEM_{AB}=-63.00, FEMBA=147.00FEM_{BA}=147.00 kNm (from Pab2/L2=100×7×9/100Pab^2/L^2=100\times7\times9/100 and Pa2b/L2=100×49×3/100Pa^2b/L^2=100\times49\times3/100)
  • BC: ∓30×8212=∓160\mp\dfrac{30\times8^2}{12}=\mp160 kNm

Stiffness equations

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

Equations of equilibrium

Joint B: MBA+MBC=0M_{BA}+M_{BC}=0; joint C: MCB+MCD=0M_{CB}+M_{CD}=0; sway (∑H=0\sum H=0, with HA=MAB+MBA−100×310H_A=\dfrac{M_{AB}+M_{BA}-100\times3}{10} and HD=MCD+MDC5H_D=\dfrac{M_{CD}+M_{DC}}{5}): MAB+MBA10+MCD+MDC5+100×710=0\dfrac{M_{AB}+M_{BA}}{10}+\dfrac{M_{CD}+M_{DC}}{5}+\dfrac{100\times7}{10}=0. When the fixed-end moments are moved to the right-hand side, the load term of this row becomes the fixed-end reaction of the 100 kN load at B, Pa2(L+2b)L3=100×49×161000=78.40\dfrac{Pa^2(L+2b)}{L^3}=\dfrac{100\times49\times16}{1000}=78.40 kN (=70+147−6310)\left(=70+\dfrac{147-63}{10}\right).

EI[1.80000.5000−0.12000.50001.8000−0.2400−0.1200−0.24000.1200]{θBθCΔ1}={13.000−160.00078.400}EI\begin{bmatrix}1.8000 & 0.5000 & -0.1200 \\ 0.5000 & 1.8000 & -0.2400 \\ -0.1200 & -0.2400 & 0.1200\end{bmatrix}\begin{Bmatrix}\theta_{B} \\ \theta_{C} \\ \Delta_{1}\end{Bmatrix}=\begin{Bmatrix}13.000 \\ -160.000 \\ 78.400\end{Bmatrix} EIθB=56.5023,EIθC=−13.5535,EIΔ1=682.7287EI\theta_{B}=56.5023,\quad EI\theta_{C}=-13.5535,\quad EI\Delta_{1}=682.7287

Final end moments (kNm, clockwise +)

EndMM (kNm)EndMM (kNm)
AB-122.33BA110.27
BC-110.27CB174.70
CD-174.70DC-169.28

Checks:

  • Joint B: (110.27) + (-110.27) = 0.00 kNm (applied clockwise moment 0.00) ✓
  • Joint C: (174.70) + (-174.70) = 0.00 kNm (applied clockwise moment 0.00) ✓

Bending moment diagram ordinates (sagging/inside tension +, kNm)

MemberBM at first endBM at second endMax within member
AB-122.33-110.2796.11 at 7.00 m
BC-110.27-174.7098.59 at 3.73 m
CD-174.70169.28-

Reactions: A: H=−31.21H=-31.21, V=136.95V=136.95, M=122.33M=122.33; D: H=−68.79H=-68.79, V=128.05V=128.05, M=169.28M=169.28.

Answer: MAB=−122.33M_{AB}=-122.33, MBA=110.27M_{BA}=110.27, MBC=−110.27M_{BC}=-110.27, MCB=174.70M_{CB}=174.70, MCD=−174.70M_{CD}=-174.70, MDC=−169.28M_{DC}=-169.28 kNm.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Theory of Structures II chapter 3 questions →