Skip to main content

Stiffness matrix: continuous beam BMD

Short answer

Using theta_B and theta_C, solve EI[[1.0714,0.25],[0.25,0.5]]{theta}={-9.694,-55}. Final moments are M_AB=-17.07, M_BA=45.44, M_BC=-45.44, M_CB=0 kNm; maximum sagging in BC is 51.53 kNm at 5.38 m from B.

Analyse the continuous beam and draw bending moment diagram which is loaded as shown in figure below. Use the stiffness matrix method. [Figure: beam ABC, A fixed; 10 kN at 2 m and 30 kN at 5 m from A; B support at 7 m (2 m + 3 m + 2 m); BC = 8 m (4 m + 4 m) with 15 kN/m UDL on the right 4 m; C roller.]

Data: A is fixed; AB = 7 m with 10 kN at 2 m and 30 kN at 5 m from A; B is a roller support; BC = 8 m with 15 kN/m on the right-hand 4 m (from 4 m to 8 m from B); C is a roller (free to rotate). EI is constant.

Unknown displacements (coordinates)

1 = θB\theta_B, 2 = θC\theta_C (clockwise +). θA=0\theta_A=0 and there are no support settlements.

Fixed-end moments (clockwise +, kNm)

  • AB (point loads): load 10 kN (a=2a=2, b=5b=5): ∓Pab2L2, ±Pa2bL2\mp\dfrac{Pab^2}{L^2},\ \pm\dfrac{Pa^2b}{L^2} = −10.204-10.204, +4.082+4.082; load 30 kN (a=5a=5, b=2b=2): −12.245-12.245, +30.612+30.612. Total: FEMAB=−22.449FEM_{AB}=-22.449, FEMBA=34.694FEM_{BA}=34.694.
  • BC (partial UDL, 15 kN/m over 4 m, from the right end): FEMBC=−25.000FEM_{BC}=-25.000, FEMCB=55.000FEM_{CB}=55.000 (integration of the point-load formulas over the loaded length).

Member stiffness equations

  • MAB=2EI7(θB)−22.45M_{AB}=\frac{2EI}{7}\left(\theta_{B}\right) - 22.45
  • MBA=2EI7(2θB)+34.69M_{BA}=\frac{2EI}{7}\left(2\theta_{B}\right) + 34.69
  • MBC=2EI8(2θB+θC)−25.00M_{BC}=\frac{2EI}{8}\left(2\theta_{B} + \theta_{C}\right) - 25.00
  • MCB=2EI8(2θC+θB)+55.00M_{CB}=\frac{2EI}{8}\left(2\theta_{C} + \theta_{B}\right) + 55.00

Stiffness matrix equation ([K]{d}={P}[K]\{d\}=\{P\})

The roller C gives MCB=0M_{CB}=0 and the joint B gives MBA+MBC=0M_{BA}+M_{BC}=0:

EI[1.07140.25000.25000.5000]{θBθC}={−9.694−55.000}EI\begin{bmatrix}1.0714 & 0.2500 \\ 0.2500 & 0.5000\end{bmatrix}\begin{Bmatrix}\theta_{B} \\ \theta_{C}\end{Bmatrix}=\begin{Bmatrix}-9.694 \\ -55.000\end{Bmatrix}

K11=4EI7+4EI8=1.0714EIK_{11}=\dfrac{4EI}{7}+\dfrac{4EI}{8}=1.0714EI, K12=K21=2EI8=0.2500EIK_{12}=K_{21}=\dfrac{2EI}{8}=0.2500EI, K22=4EI8=0.5000EIK_{22}=\dfrac{4EI}{8}=0.5000EI.

EIθB=18.8140,EIθC=−119.4070EI\theta_{B}=18.8140,\quad EI\theta_{C}=-119.4070

Final end moments (kNm)

EndMM (kNm)EndMM (kNm)
AB-17.07BA45.44
BC-45.44CB0.00

Checks: Joint B: 45.44+(−45.44)=0.0045.44+(-45.44)=0.00 ✓; roller C: MCB=0.00M_{CB}=0.00 ✓.

Bending moment diagram (sagging +, kNm)

MemberBM at first endBM at second endMax within member
AB-17.07-45.4411.23 at 5.00 m
BC-45.440.0051.53 at 5.38 m
BM under the 10 kN load: 6.25; under the 30 kN load: 11.23 kNm.

Reactions: A: H=0.00H=0.00, V=11.66V=11.66, M=17.07M=17.07; B: H=0.00H=0.00, V=49.02V=49.02, M=0.00M=0.00; C: H=0.00H=0.00, V=39.32V=39.32, M=0.00M=0.00.

Answer: MAB=−17.07M_{AB}=-17.07, MBA=45.44M_{BA}=45.44, MBC=−45.44M_{BC}=-45.44 kNm, MCB=0M_{CB}=0; maximum sagging in BC = 51.53 kNm at 5.38 m from B.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Theory of Structures II chapter 3 questions →