𝐴𝑅1

𝐴𝑀1 −1 𝐴𝑅2
𝐴𝑀2 12 𝐴𝐽1
=
𝐴𝑀3 − 3Τ5 𝐴𝐽2
𝐴𝑀4 15
𝐴𝐽3

𝐴𝑅1
𝐴𝑀1 −1 0 𝐴𝑅2
𝐴𝑀2 12 0 𝐴𝐽1
=
𝐴𝑀3 − 3Τ5 − 4Τ5 𝐴𝐽2
𝐴𝑀4 15 4
𝐴𝐽3

𝐴𝑅1
𝐴𝑀1 −1 0 0 0 0 𝐴𝑅2
𝐴𝑀2 12 0 −1 0 0 𝐴
= 𝐽1
𝐴𝑀3 − 3Τ5 − 4Τ5 0 − 3Τ5 0 𝐴
𝐽2
𝐴𝑀4 15 4 −1 3 −1
𝐴𝐽3
576.0000 72.0000 0 0
72.0000 12.0000 0 0
0 0 41.6667 12.5000
0 0 12.5000 5.0000

1.0e+03 *

1.4910 0.1400 -0.1395 0.1050 -0.0675

0.1400 0.0267 -0.0100 0.0200 -0.0100
-0.1395 -0.0100 0.0170 -0.0075 0.0050
0.1050 0.0200 -0.0075 0.0150 -0.0075
-0.0675 -0.0100 0.0050 -0.0075 0.0050

1491 140 −139.5 105 −67.5

140 26.67 −10 20 −10
−139.5 −10 17 −7.5 5
105 20 −7.5 15 −7.5
−67.5 −10 5 −7.5 5
 1491 140 −139.5 105 −67.5
140 26.67 −10 20 −10
−139.5 −10 17 −7.5 5
105 20 −7.5 15 −7.5
−67.5 −10 5 −7.5 5

2.5992
-0.5208

57.6190
0
4.7619
 Fs remains same as calculated earlier. Independent of loading.

Fully restraint system

Problem which is solved Original problem to be solved
with fix end actions
Problem 2:

Simple flexibility approach:

Post Processing:
Generalized flexibility approach:

Equivalent joint load system

Unassembled member flexibility matrix

1Τ3 1Τ2 𝐿 0 0 0
Τ
𝐿3 1 2 𝐿 1Τ𝐿2 0 0 0
𝐹𝑀 = 0 0 1Τ3 1Τ2 𝐿 0
𝐸𝐼
0 0 1Τ2 𝐿 1Τ𝐿2 0
0 0 0 0 1Τ3
Structural system: Member action system:
Simple stiffness approach:
Generalized stiffness approach:

12𝐸𝐼 6𝐸𝐼
−
𝐿3 𝐿2
6𝐸𝐼 4𝐸𝐼
− 2
𝐿 𝐿
Displacement transformation matrix

0 0 −1 −𝐿 1 0
1 0 0 −1 0 0
𝐵𝑀𝑆 = −𝐿 0 0 0 −1 1
−1 1 0 0 0 0
Computer implementation of stiffness approach for 2D Frames

Input file 1
Nodes x y ∆x ∆y θz Px Py Mz
1 0 0 1 1 1 0 0 0
2 0 L 0 0 0 0 0 0
3 L L 0 0 0 0 0 0
4 L 0 1 1 1 0 0 0

Input file 2

Members Node j Node k EI AE Loading

UDL Point load
1 1 2
2 2 3
3 3 4
Step 1: continued
Step 2: Element stiffness and Equivalent
joint load in local coordinates

𝐴𝐸 𝐴𝐸
0 0 − 0 0
𝐿 𝐿
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼
0 0 − 0
𝐿3 𝐿2 𝐿3 𝐿2 𝑤𝐿 𝑃𝑏
6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼 − −
0 0 − 2 2 𝐿
𝑖 𝐿2 𝐿 𝐿 𝐿
𝑆𝑚 =
𝐴𝐸 𝐴𝐸 𝑤𝐿2 𝑃𝑎𝑏 2
− 0 0 0 0 − − 2
𝐿 𝐿 𝐴𝑖𝑚 = 12 𝐿
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼 0
0 − 3 0 𝑤𝐿 𝑃𝑎
𝐿 𝐿2 𝐿3 𝐿2 − −
6𝐸𝐼 2𝐸𝐼 6𝐸𝐼 4𝐸𝐼 2 𝐿
𝑤𝐿2 𝑃𝑎𝑏 2
0 − 2 0
𝐿 𝐿 𝐿2 𝐿 + 2
12 𝐿
Step 2: Element stiffness and Equivalent joint load in structure/global coordinates

𝑅 𝑇 0 𝛥𝑖
𝛥𝑖𝑠 = 𝑚
0 𝑅𝑇
Step 3: Assembly of stiffness matrix and joint load vector
Step 4: Rearrangement of matrices and analysis
Implementation to grid structures/ 3D Frames/ space trusses
 Work Sheet I
Advanced Structural Analysis (CE6L302) - autumn 2021
Note: Use MATLAB/ Octave for matrix operations.
Problem 1: Define the set of reactions at fixed end A and B of the given beam. The flexural rigidity
of the beam is EI. Use simple flexibility approach.

Fig. 1
Problem 2: The flexural rigidity for both the members in Fig. 2 is same as EI. Complete the
structural analysis using generalized flexibility approach.

Fig. 2 Fig. 3
Problem 3: For the problem given in Fig. 3, calculate external support reactions. Material is same
for all the members.
Problem 4: For the Fig. 4, support A undergoes rotation in clockwise direction by given amount
and support B settles down in the downward direction by Δ amount. Calculate the support reaction
using simple flexibility approach.
 Fig. 4
Problem 5: Calculate the support reactions as well as all member forces in the given pin jointed
member (Fig. 5) using generalized flexibility approach.

Fig. 5