Stiffness Method (Notes - Beams, Frames and Truss) PDF
Stiffness Method (Notes - Beams, Frames and Truss) PDF
Outline
Stiffness Method
➢ Analysis of Beams
➢ Analysis of Frames
➢ Analysis of Truss
Chapter 2: Stiffness Method Topic 2.0 : Outline
Timber Frame
Chapter 2: Stiffness Method Topic 2.0 : Learning Outcomes
(CO2PO10)
Chapter 2: Stiffness Method Topic 2.3 : Truss Analysis
q N AE 1 − 1 d N
q = − 1 1 d
F L F
q = k'd Where,
AE 1 − 1
q = member end forces along
k'= local axis-x
L − 1 1 k’ = member stiffness matrix
d = displacement along local
axis-x
Chapter 2: Stiffness Method Topic 2.3 : Truss Analysis
Chapter 2 Stiffness Method Topic 2.3 : Truss Analysis
DN x
d N C x C y 0 0 DN y
d = 0 0 C C D
F x y
Fx
DF
y
d = TD
Chapter 2: Stiffness Method Topic 2.3: Truss Analysis
TT is use to transforms local forces q acting at the ends of the member into global for components Q
C x qF
Fx
QN x = q N C x ; QN y = q N C y QF 0 C y
y
QFx = qF C x ; QFy = qF C y
Q = TTq
Chapter 2: Stiffness Method Topic 2.3 : Truss Analysis
Substitute 1 into
Q = TTq
Q = T T k ' TD or Q = kD
Where;
k = T T k'T
Performing the matrix operations yields:-
Cx 2 CxC y − Cx
2
− CxC y
2
− CxC y − C y
2
AE C x C y Cy
k=
L − Cx − CxC y CxC y
2 2
Cx
− C x C y − C y
2 2
CxC y C y
Chapter 2 : Stiffness Method Topic 2.3 : Truss Analysis
xF − x N
C x = cos x =
L
xF − x N
=
( xF − x N ) 2 + ( y F − y N ) 2
yF − yN
C y = cos y =
L
yF − yN
=
( xF − x N ) 2 + ( y F − y N ) 2
Chapter 2 Stiffness Method Topic 2.3 : Truss Analysis
xF − x N
C x = cos x =
L
yF − yN
C y = cos y =
(XN , YN) L
x
N L = Length of member
Example 1 : Calculate the value of Cx and Cy by the coordinate system for the
truss shown in Figure Q1
Solution:
Figure Q1
Chapter 2 Stiffness Method Topic 2.3 : Truss Analysis
3−0 3−0
(3 , 4)
Cx = =1 Cx = = 0.6
3 5
0−0 4−0
Cy = =0 Cy = = 0.8
3 5
(0 , 0) (3 , 0)
0 − (−3) 0 − (−3)
(0 , 0) Cx = =1 Cx = = 0.6
3 5
− 4 − (−4) 0 − (−4)
Cy = =0 Cy = = 0.8
3 5
(-3 , -4) (0 , -4)
Chapter 2 : Stiffness Method Topic 2.3 : Truss Analysis
Member 1 , L=3m
x = 0
C x = cos x = cos 0 = 1
N F
1 C y = sin x = cos 0 = 0
Member 1 , L=5m
F
x = 53.13
C x = cos x = cos 53.13 = 0.6
2
x = 53.13 C y = sin x = sin 53.13 = 0.8
N
Chapter 2: Stiffness Method Topic 2.3 : Truss Analysis
Qk = K11 Du + K12 Dk ; Qu = K 21 Du + K 22 Dk
Often Dk = 0, since support are not displaced;
Qk = K11 Du ; Qu = K 21Du
Du = K11 Qk
−1
Member forces,
q = k ' TD
Chapter 2 : Stiffness Method Topic 2.3 : Truss Analysis
DNx
qN AE 1 − 1 C x C x 0 0 DNy
q = L − 1 1 0 0 C C x DFx
F
x
DFy
Since qN = -qF for equilibrium, only one of the forces has to be found
DNx
D
qF =
AE
L
− Cx − C y Cx C y
Ny
DFx
DFy
Chapter 2: Stiffness Method Topic 2.3 : Truss Analysis
Ny 0 = AET y
(QFx ) 0 − Cx
(QFy ) 0 − C y
(QNx ) 0 Cx
(Q ) C Qo for Fabrication Error
Ny 0 = AEL y
(QFx ) 0 L − Cx
Fy 0
(Q ) −
y
C
Chapter 2 :Stiffness Method Topic 2.3 : Truss Analysis
Temperature change,
DNx
D
qF =
AE
L
− Cx − C y Cx Cy Ny
DFx
− AET
DFy
Fabrication Error,
DNx
D
AEL
qF =
AE
L
− Cx − C y Cx Cy Ny
DFx
−
L
DFy
1
Example 2:
A pin-jointed plane truss is shown in Figure Q2. Member 4 was made 0.01 too short before it was fitted into
place. Assume axial rigidity AE is constant for all members, and A = 2000 mm2 and E = 210GPa.
a) Establish the element global stiffness matrices. (10 marks)
b) Construct the stiffness matrix of the structure. (5 marks)
c) Determine the unknown nodal displacements at 3. (4 marks)
d) Determine the axial forces in members 1, 2 and 3. (6 marks)
75 kN
5
3 4
50 kN
Solution:
3 5m
2
1
1 2
10 m
Figure Q2
Chapter 2 : Stiffness Method Topic 2.3 : Truss Analysis
1
(0,0) 2
1 (10,0)
10 m 4
5 7
Chapter 2: Stiffness Method Topic 2.3 : Truss Analysis
0 − 0.1
7
0.1 0
5 (10,0)
N 1 5
F
0 0 0 0
k1 = AE
6
Cx = 0
2 0 0 0 0 5
0 0.2 0 − 0.2
k 2 = AE
6
5−0 0
Cy = =1 0 0 0 1
6 5
0 − 0.2 0 0.2 2
5
N
(0,0)
Chapter 2 : Stiffness Method Topic 2.3 : Truss Analysis
Member 3, L = 5m 2
4 Member 4, L = 11.18m
(0,5)
F
3 1
N
(10,5)
4
10 − 10 3
Cx = =0 8
5
7
(10,0)
5−0 8 F
Cy = =1
5 7 10 − 0 0−5
Cx = = 0.894 Cy = = −0.447
N
(10,0)
11.18 11.18
7 8 3 4 1 2 7 8
0 − 0.2 0.036 − 0.018 − 0.036 0.018
8
0 0.2 4
Chapter 2: Stiffness Method Topic 2.3 : Truss Analysis
8
10 − 0 5−5
Cx = =1 Cy = =0 (Q1 )O 0.894
10 7
10 (Q ) − 0.447 (10,0)
= AE ( −0.01) F
Combine:
2 O
Question 2(c)
385.832 71820 − 15120 D1
=
− 242. 916 − 15120 91560 D 2
D1 4.988(10−3 )m
= −3
D 2 − 1.829(10 )m
DNx
D
AEL
− Cx − C y Cx C y −
Question 2(d) AE
qF =
Ny
(CO2PO10)
1
Kinematic indeterminacy
➢ If consider effects of both bending & shear
- Each node on the beam can have 2 degrees of freedom, vertical displacement & rotation
➢ These displacement will be identified by code numbers
➢ The lowest code numbers will be used to identify the unknown displacement & the highest numbers
are used to identify the known displacement
q = kd
Chapter 2: Stiffness Method Topic 2.1 : Beam Analysis
Chapter 2 : Stiffness Method Topic 2.1 : Beam Analysis
Local Global
12 EI 6 EI 12 EI 6 EI
+ L3 + L2 − L3 + L2
6 EI 4 EI 6 EI 2 EI
+ 2 + − 2 +
k'= k = L L L L
− 12 EI − 6 EI + 12 EI − 6 EI
L3 L2 L3 L2
6 EI 2 EI 6 EI 4 EI
+ 2 + − 2 +
L L L L
1
Q = KD
Actual loading
Qk = K11Du + K12 Dk
Qu = K 21Du + K 22 Dk
Member Forces, q
Fixed-end element loading
q = kd + q0 of joints
Example 3 : Calculate the intermediate loadings for the beam shown in Figure Q3
5 1 3 7
20 kN/m
6 2 4 8
Figure Q3
1
I 2 3I
A B C 3 D
2m
3m 3m
5 1
(Q5 ) O 0
Element 1 6 2 (Q ) 0
6 O =
A
I
(Q1 ) O 0
1m 2 O 0
(Q )
Chapter 2: Stiffness Method Topic 2.1 : Beam Analysis
wL/2 wL/2
Element 2 Element 3
w kN/m 1 3 7
3
20 kN/m
20 kN/m 8
2 4 4
A B
L
wL2/12 wL2/12 C D
B I C 3I
wL 2m
2 3m
RA wL2
FEM
AB (Q1 ) O 20kN (Q3 )O 30kN
= 12 (Q ) 6.67kNm (Q ) 15kNm
RB wL
2 2 O =
4 O =
FEM BA
(Q3 )O 20kN (Q7 )O 30kN
wL
2
− 12 4 O − 6.67kNm
(Q ) 8 O − 15kNm
(Q )
Chapter 2 : Stiffness Method Topic 2.1 : Beam Analysis
q = kd + q0
12EI 6 EI 12EI 6 EI
+ L3 + 2
L
− 3
L
+ 2
L d (q )
qNY 6 EI 4 EI 6 EI 2 EI NY NY 0
q + 2 + − 2 + d (q )
NZ = L L L L NZ + NZ 0
qFY 12EI 6 EI 12EI 6 EI d FY (qFY )0
− L3 − 2 + 3 − 2
L L L d (qFZ )0
qFZ 6 EI 2 EI 6 EI 4 EI FZ
+ 2 + − 2 +
L L L L
1
Example 4 :
A continuous beam ABCD is subjected to the loads as shown in Figure Q4 and the flexural rigidity, EI, of
the member is as indicated. When the loads are acting, support B yields vertically downwards by 60/EI.
a) Establish the element stiffness matrices and the equivalent element nodal force vectors in the global
coordinate system. (12 marks)
b) Form the global stiffness relations corresponding to the unconstrained degrees of freedom and solve
for the non-zero degrees of freedom. (8 marks)
c) Determine the axial forces for element 1 only. (5
marks)
100 kN
30 kN/m
20 kN/m
A D
1 E 2 3
3m 3m 8m 6m
Solution :
Chapter 2 : Stiffness Method Topic 2.1: Beam Analysis
Assembling k1 , k2 and k3
Chapter 2 : Stiffness Method Topic 2.1 : Beam Analysis
Q = KD + Qo
Chapter 2 : Stiffness Method Topic 2.1 : Beam Analysis
Question 4 (b)
Chapter 2: Stiffness Method Topic 2.1 : Beam Analysis
q = kd + q0
q5 0.1111 0.3333 − 0.1111 0.3333 0 50
q 0.3333 1.3333 − 0.3333 0.6667 75
6 = EI 0
+
q4 − 0.1111 − 0.3333 0.1111 − 0.3333 − 60 / EI 50
q2 0.3333 0.6667 − 0.3333 1.3333 − 22.35 / EI − 75
q5 49.22kN
q 80.10kNm
6 =
q4 50.78kN
q2 − 84.80kNm
Chapter 2: Stiffness Method Topic 2.0 : Learning Outcomes
(CO2PO10)
1
The resulting 6 load-displacement relations for the member can be expressed in matrix form or can
be written in the form of : -
q = k'd
Displacement
transformation
matrix
Cx = cos x ; C y = cos y
Chapter 2: Stiffness Method Topic 2.2 : Plane Frame Analysis
Load-displacement relations
1
d = TD
d Nx' C x C y 0 0 0 0 DNx
d
Ny ' − C y C x 0 0 0 0 DNy
d ' 0 0 1 0 0 0 DNz
=Nz
d Fx' 0 0 0 Cx C y 0 DFx
d ' 0 0 0 − C y Cx 0 DFy
Fy
d Fz' 0 0 0 0 0 1 DFz
QNx C x − C y 0 0 0 0 q Nx'
Q C q
Ny y C x 0 0 0 0 Ny '
QNz 0 0 1 0 0 0 q Nz '
= q
QFx 0 0 0 Cx − C y 0 Fx'
QFy 0 0 0 C y Cx 0 qFy'
QFz 0 0 0 0 0 1 qFz'
Frame-member global stiffness matrix
q = k ' TD + qo
Q = T T k'TD + Qo
Q = kD Where k = TTk’T
Chapter 2 : Stiffness Method Topic 2.2 : Plane Frame Analysis
q = k ' TD + qo
1
Example 5 :
A rigid-jointed frame ABC is subjected to moment and horizontal load at B as shown in Figure Q5. Support A
rotates in counterclockwise direction by 80/EI. All the members have the following properties:
EI = 50,000 kNm2
EA = 1600,000 kN
45 kNm
60 kN C
B 2 Solution :
1
6m Figure Q5
4m
Chapter 2 : Stiffness Method Topic 2.2: Plane Frame Analysis
Kinematic Indeterminacy
2 8
If consider effects of axial loads qNx’. qFx’,
3
45kNm 9 shear loads qNy’, qFy’ and bending
moment qNz’, qFz’
1 7
60kN B 2 C
d.o.f = 3 The code numbers are specified with
d.o.f = 0
numbers assigned first to the
unconstrained degrees of freedom.
1
Degree of freedom is 3
Support A rotates in
counterclockwise
direction by 80/EI
5
6
A 4
d.o.f = 0
Chapter 2: Stiffness Method Topic 2.2 : Plane Frame Analysis
Question 5 (b)
Chapter 2 : Stiffness Method Topic 2.2 : Plane Frame Analysis
Q8 0 -9.4 - 0 0 0 0 9.4 - 0 0 8
18.8 18.8
Q9 0 18.8 25.0 0 0 0 0 - 50.0 0 0 9
18.8
Chapter 2: Stiffness Method Topic 2.2 : Plane Frame Analysis
Support Reaction
Q4 -2.78 0 8.3 0 0 0
-0.1639
Q5 0 -267 8.3 8.3 0 0
Q6 33.3 1.6x10- 0
8.3 0 16.7 3
= 103 0.0595 + 103 +
Q7 -400 0 0 0 0 0
Q8 0 -9.4 -18.8 0 0 0
-0.8567
Q9 0 18.8 25.0 0 0 0
Q4 -6.65kN
Q5 -22.94kN
Q6 -15.61kNm The moment reaction at C, Q9 = -20.30 kNm
=
Q7 65.56kN
Q8 15.55kN
Q9 -20.30kNm