Skip to main content

Stiffness matrix: frame with 100 kN, end reactions

Short answer

The frame's 4×4 stiffness matrix is assembled from member stiffness equations ignoring axial deformation. Solving for joint rotations and sway gives end moments M_BA=-19.63, M_BC=19.63, M_CB=53.27, M_CD=-53.27, M_DC=-140.19 kNm and horizontal reactions H_A=3.27 kN, H_D=96.73 kN to the left.

Generate the stiffness matrix for the frame shown and determine the end reactions at the support. [Figure: portal frame ABCD; 100 kN horizontal at B; beam BC (EI) of 6 m; columns EI, AB 6 m high and CD 2 m high; A hinged, D fixed.]

Assumptions: A(0,0) is a hinge, B(0,6) and C(6,6) are the top corners (AB = 6 m, BC = 6 m, CD = 2 m high, so D is 4 m above A and fixed). EI is the same for all members; the 100 kN load at B acts to the right. Axial deformations are ignored.

Degrees of freedom (coordinates)

  1. θA\theta_A (hinge rotation) 2. θB\theta_B 3. θC\theta_C 4. Δ\Delta (horizontal sway of the girder). Chord rotations: ψAB=Δ/6\psi_{AB}=\Delta/6, ψCD=Δ/2\psi_{CD}=\Delta/2, ψBC=0\psi_{BC}=0.

Member stiffness equations (Mij=2EIL(2θi+θj−3ψ)M_{ij}=\tfrac{2EI}{L}(2\theta_i+\theta_j-3\psi), no member loads)

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

Equilibrium equations

  • Hinge A: MAB=0M_{AB}=0
  • Joint B: MBA+MBC=0M_{BA}+M_{BC}=0
  • Joint C: MCB+MCD=0M_{CB}+M_{CD}=0
  • Storey shear: MAB+MBA6+MCD+MDC2+100=0\dfrac{M_{AB}+M_{BA}}{6}+\dfrac{M_{CD}+M_{DC}}{2}+100=0

Stiffness matrix of the frame [K]{d}={F}[K]\{d\}=\{F\}

EI[0.66670.33330.0000−0.16670.33331.33330.3333−0.16670.00000.33332.6667−1.5000−0.1667−0.1667−1.50001.5556]{θAθBθCΔ1}={0.0000.0000.000100.000}EI\begin{bmatrix}0.6667 & 0.3333 & 0.0000 & -0.1667 \\ 0.3333 & 1.3333 & 0.3333 & -0.1667 \\ 0.0000 & 0.3333 & 2.6667 & -1.5000 \\ -0.1667 & -0.1667 & -1.5000 & 1.5556\end{bmatrix}\begin{Bmatrix}\theta_{A} \\ \theta_{B} \\ \theta_{C} \\ \Delta_{1}\end{Bmatrix}=\begin{Bmatrix}0.000 \\ 0.000 \\ 0.000 \\ 100.000\end{Bmatrix}

The matrix is symmetric. The only non-zero load is the 100 kN at the sway coordinate.

Solution

EIθA=44.8598,EIθB=−14.0187,EIθC=86.9159,EIΔ1=151.4019EI\theta_{A}=44.8598,\quad EI\theta_{B}=-14.0187,\quad EI\theta_{C}=86.9159,\quad EI\Delta_{1}=151.4019

End moments (kNm, clockwise on the member end +)

EndMM (kNm)EndMM (kNm)
AB0.00BA-19.63
BC19.63CB53.27
CD-53.27DC-140.19

Checks: MAB=0.00M_{AB}=0.00 ✓; Joint B: −19.63+(19.63)=0.00-19.63+(19.63)=0.00 ✓; Joint C: 53.27+(−53.27)=0.0053.27+(-53.27)=0.00 ✓.

End reactions

Column shears: HA=MAB+MBA6H_A=\dfrac{M_{AB}+M_{BA}}{6} and HD=MCD+MDC2H_D=\dfrac{M_{CD}+M_{DC}}{2} (magnitudes), and the vertical reactions from the girder moments.

  • Support A: HA=3.27H_A=3.27 kN (to the left), VA=−12.15V_A=-12.15 kN, MA=0M_A=0
  • Support D: HD=96.73H_D=96.73 kN (to the left), VD=12.15V_D=12.15 kN, MD=140.19M_D=140.19 kNm (anticlockwise +)
  • Check: ∑H=100−3.27−96.73=0.00\sum H=100-3.27-96.73=0.00 ✓, ∑V=0.00\sum V=0.00 ✓
MemberBM at first endBM at second endMax within member
AB0.0019.63-
BC19.63-53.27-
CD-53.27140.19-

Answer: EIΔ=151.402EI\Delta=151.402; MBA=−19.63M_{BA}=-19.63, MBC=19.63M_{BC}=19.63, MCB=53.27M_{CB}=53.27, MCD=−53.27M_{CD}=-53.27, MDC=−140.19M_{DC}=-140.19 kNm; reactions HA=3.27H_A=3.27 kN, HD=96.73H_D=96.73 kN (both to the left).

✓ Answer checked by Claude Opus 5.5 (high effort)

All Theory of Structures II chapter 3 questions →