1.

Derive flexibility & stiffness matrix

2
1

To find flexibility matrix

[f] =
[ f 11 f 12
f 21 f 22 ]
1
1 To get I column
1

F21
F11

1
+

The conjugate beam with load

1/EI

B1
A1
RA RB
1 1

ε M FB=0
1
ε F y =0

1 1 2 1
( *L* ) L +RA1*L = 0 ½ *L* +RA1 + RB1 = 0
2 EI 3 EI
−l −l
R A1 = (↓) RB1 = (↓)
3 EI 6 EI

Therefore slope at A = SF at A1
 −l
-1 =
l
3 EI
f11 =
3 EI

Slope at B = SF at B1

l
1=
6 EI
l
F21 =
6 EI

1
2 To get II column

A 1
B

f12 F22

1
+

The conjugate beam

1
EI

l
6 EI l
3 EI
Slope at A = SF at A1

−l
-1 =
6 EI
l
F21 =
6 EI

Slope at B = SF at B1

l
F22 =
3 EI
 l
= 1
3 EI

[ ]
l l
3 EI 6 EI
[f] =
l l
6 EI 3 EI

To find stiffness matrix:-

1
1 2
EI

[K] =
[ K 11 K 12
K 21 K 22 ]
To get
K I column
11

K 21

1
K 11 +¿ K 21

conjugate beam
K 11 /EI
B1
A1 K 21 / EI

Deflection at A = B.M at A1 = 0

1 2
-½ ( K 11 /EI )L* L – (1/2 K 11 /EI L) L= 0
3 3

k 11=2 k 21

ε F y =0
 K 11 1 K 21
R A1 + ½ l− l=0
EI 2 EI

k 11 l K 21
R A1 + - l=0
2 EI 2 EI

k 11 l
R A1 + (↓)
4 EI

Slope at A = SF @ A1

−k 11 l
-1 =
4 EI

−2 EI
K21 =
l

Fixed end moments

w
1. A B

l
MFBA
MFAB

2 2
−w l wl
MFAB = MFBA =
12 12

P
A B

MFAB
MFBA
L/2 L/2

−pl pl
MFAB = MFBA =
8 8

3. p
A B

MFAB
MFBA
a
 b

2 2
−wa b +w a b
MFAB = 2 MFBA = 2
l l
A B
M
4.

Mb 3b Ma 3a
MFAB = [2− ] MFBA = [2− ]
l l l l
MFAB MFBA
5.

A
B

MFAB MFBA
l

2 2
−w l wl
MFAB = MFBA =
30 20
A w
6. B

a C b
MFAB MFBA

−w 2 2 3
MFAB = 2 [12ab(b+c) + 6cb +4c (a+b) + c ]
12l

w 2 2 3
MFBA = 2 [12ab(b+c) + 6ca +4c (a+b) + c ]
12l

7.A B
w

MFAB MFBA
 5 5
MFAB = −w l 2∗¿ 96 ¿ MFBA = −w l 2∗¿ 96 ¿

2.Repeat the problem (1) by choosing moment at A & B as redundant

10 0KN 10 KN/m
B
A
2I C I

5m 5m 10m

Redundant:

1
3
4

2
1x 3x 4x
2x

[F] {P} = {Δ}

[F] = [A]T [Fx] [A]

1.To find [A]:

1 5KN-m 10* ½ = 5KN-m
1 A
C B

RA = ½ RB= ½ ½ ½ ½ ½

{}
0
5
ε F y =0 A1 =
−5 3
0

RA + RB + 1 = 0

1
C
A
RA = - ½

2 0.5KN-m 1/20*10 = 0.5

A
B C B

1/20 1/2
1/20

{}
0
0.5
ε M A =0 A1 =
0.5
0

-RB *20+ 1 = 0

RB = +1/20

2. 0.5 0.5
1 1

1/20 1/20 1/20 1/20 1/20 1/20

{ }
0
−0.5
A3 =
0.5
0
1/20*10=0.
0.5 5 1
1

1/20 1/20
ε M A =0

{ }
0
−0.5
-RB*20 + 1 =0 A3 =
0.5
0

RB = 1/20
 [ ]
0 0 1 0
5 0.5 −0.5 0.5
[A] =
−5 0.5 0.5 −0.5
0 0 0 1

From previous problem

[ ]
1.67 −0.833 0 0
1 −0.833 1.67 0 0
[Fx] =
EI 0 0 3.33 −1.67
0 0 −1.67 3.33

[F] = [A]T [Fx] [A]

Consider [Fx] [A]

[ ][ ]
1.67 −0.833 0 0 0 0 1 0
1 −0.833 1.67 0 0 5 0.5 −0.5 0.5
[Fx][A] =
EI 0 0 3.33 −1.67 −5 0.5 0.5 −0.5
0 0 −1.67 3.33 0 0 0 1

[ ]
−4.165 −0.4165 2.08 −0.416
1 8.35 0.835 −1.665 0.835
EI −16.65 1.665 1.665 −3.33
8.35 −0.835 −0.835 4.16

[ ][ ]
0 5 −5 0 −4.165 −0.4165 2.08 −0.416
1 0 0.5 0.5 0 8.35 0.835 −1.665 0.835
[F] =
EI 1 −0.5 0.5 0 −16.65 1.665 1.665 −3.33
0 0.5 −0.5 1 8.35 −0.835 −0.835 4.16

[ ]
125 −4.15 −16.66 20.82
1 −4.15 1.25 0 −1.25
[F] =
EI −16.66 0 3.74 −2.5
20.85 −1.25 −2.5 6.24

We know that

[F]{P} = {Δ}

[ ]
125 −4.15 −16.66 20.82
1 −4.15 1.25 0 −1.25
{ P }={ Δ }
EI −16.66 0 3.74 −2.5
20.85 −1.25 −2.5 6.24
To get force vector

{ }{ }{ }
P1 0 50+50
P2 30
= − 125−83.33
P3 MA −125
P4 MB 83.33

{ }{ }
P1 −100
P2 −11.67
=
P3 M A +125
P4 M B +83.33

To find MA & MB

[ ]{ } { }
125 −4.15 −16.66 20.82 −100 D1
1 −4.15 1.25 0 −1.25 −11.67
= D2
EI −16.66 0 3.74 −2.5 M A +125 0
20.85 −1.25 −2.5 6.24 M B + 83.33 0

[−16.66
20.85
0
]{ }[ ]{
−100 + 3.74 −2.49 M A + 125 = 0
−1.25 −11.67 −2.5 6.24 M B−83.33 0 } {}
1
{ }[
1666 + 3.74 −2.49 M A + 125 = 0
EI −2070.41 −2.5 6.24 M B−83.33 0 ]{ } {}
[−2.5
3.74
]{
−2.49 M A +125
6.24 M B −83.33 }{
= −1666
2070.41 }
{ M A +125
M B−83.33 }{
= −304.30
209.10 }
To find force:

{}
P1
P2 = P
{ FEM }+ [ A ] [P]
P3
P4

{ }[ ]{ }
−125 0 0 1 0 −100
125 −5 0.5 −0.5 0.5 −11.67
= +
−83.33 5 0.5 0.5 −0.5 −304.3
83.33 0 0 0 1 209.10
 { }{ }
−125 −304.30
125 + −248.1
=
−83.33 236.43
83.33 209.7

{ }{ }
P1 −430
P2 −123.10
=
P3 153.10
P4 293.03

Analyze the France by flexibility method & draw BMD

B 2I
HC
3m
6
10 KN/m Im
VC

A
DOIm= 5-3 = 2

In this case co-ordinates 3 Is acting only at C. Therefore no need of finding reaction at B

1 2
3
HC

1
VC Released structure
4

3x 4x

To find [A] Local Co-ordiantes

1x
1. 1 C 1 0 0
B

0 1
A
 {}
−1
−1
A1 =
0
0 1 1 1

2.

0 1

{}
−1
1
A2 =
−1
1

3.

{}
−6
0
A3 =
0
0
3.

3m
 {}
3
−3
A4 =
3
0

[ ]
−1 −1 −6 3
1 1 0 −3
[A] =
0 −1 0 3
0 1 0 0

To find [Fx]

[ ]
l −l
1 3 6
[f] =
EI −l l
6 3

[ ]
6 −6
[f]1 =
1 3
EI −6
6
6
=
1 2 −1
EI −1 2 [ ]
6 3

[ ]
3 −3
[f]2 =
1 3
2 EI −3
6
3
=
1 0.5
EI −0.25 [
−0.25
0.5 ]
6 3

x
[F ] =
[ [−12 −12 ]
0 [−0.25
0.5
0

−0.25
0.5 ] ]
To find [F] = [A]T [Fx] [A]

[ ][ ]
2 −1 0 0 −1 −1 −6 3
1 −1 2 0 0 1 1 0 −3
[F] =
EI 0 0 0.5 −0.25 0 −1 0 3
0 0 −0.25 0.5 0 1 0 0

[ ]
−3 −3 −12 9
3 3 6 −9 1
0 −0.75 0 1.5 EI
0 0.75 0 −0.75
 [ ]
6 6 18 −18
6 7.5 18 −20.25 1
[F] =
18 18 72 −54 EI
−18 −20.25 −54 58.5

To find Δ:

[F] {P} = {Δ}

{ }{ }{ }{ }
P1 0 30 −30
P2 30 0
= − 0 =
P3 H C 0 H C
P4 V C 0 V C

[ ]{ } { }
6 6 18 −18 −30 D1
1 6 7.5 18 −20.25 0
= D2
EI 18 18 72 −54 HC 0
−18 −20.25 −54 58.5 VC 0

[−18
18 18
]{ } [
−20.25 0 −54 58.5 V C]{ } { }
−30 + 72 −54 H C = 0
0

{−540
540 }+
[ 72
−54 ]{ } { }
−54 H C
58.4 V C
=0
0

[−54
72
]{ } { }
−54 H C
58.4 V C
= 540
−540

{ }{ }
HC
VC
= 1.843
−7.5

{P} = {PFEM} + {A} {P}

{ }[ ]{ }
−30 −1 −1 −6 3 −30
30 + 1 1 0 −3 0
=
0 0 −1 0 3 1.843
0 0 1 0 0 −7.5

{ }{ } { }
−30 −3.558 −33.55
30 + −7.5 = 22.5
=
0 −22.5 −22.5
0 0 0