CES521 - 3 - Stifeness Method (Truss)
CES521 - 3 - Stifeness Method (Truss)
ECS 468
INDETERMINATE STRUCTURES
Topic 2
Stiffness Method
(Truss)
AE AE
q' 'N = − dF ; q' 'F = dF
L L
q N = AE 1 − 1 d N
q F L − 1 1 d F
q = k'd eqn 14.3
AE 1 − 1 This matrix, k’ is called the member stiffness
k'= eqn 14.4 matrix
L − 1 1 qN = – qF
CES 521 - INDETERMINATE STRUCTURES Slide 6 of 49 Topic 2 – Stiffness Method (Truss)
Displacement & Force Transformation matrices
(Pulau Pinang)
λ x = cos θ x ; λ y = cos θ y
Consider member NF of the truss as shown
in Figure (b), the coordinates of N & F are
(xN, yN) and (xF, yF) respectively Figure (b)
xF − xN
λ x = cos θ x =
L
xF − xN
= eqn 14.5
(xF − xN )2 + ( yF − y N )2
yF − yN
λ y = cos θ y =
L
yF − yN
= eqn 14.6
(xF − xN )2 + ( yF − y N )2
d N = D N x cos θ x + D N y cos θ y
d F = D Fx cos θ x + D F y cos θ y
d N = D N x λ x + D N y λ y ; d F = D Fx λ x + D F y λ y
In matrix form,
DNx
d N = λ x λ y 0 0 DN y
eqn 14.7
d F 0 0 λ x λ y D Fx
D F y
d = TD eqn 14.8
Q Fx = q F cos θ x ; Q F y = q F cos θ y
Q N x λ x 0
Q N λ y 0 q N
y= 0 eqn 14.10
Q Fx λ x q F
0 λ y
F y
Q
– In this case, TT transforms the 2 local forces q acting at the ends of the
member into 4 global force components Q
– This force transformation matrix is the transpose of the displacement
transformation matrix
If we subtitute equation 14.8 into equation 14.3, we can determine the member’s
forces q in terms of the global displacement D at its end points.
q = k’TD ---------- eqn. 14.13
Subtitute this equation into equation 14.11 yields the final result:
Q = TTk’TD
Once all the member stiffness matrices are formed in the global coordinates, it
becomes necessary to assemble them in the proper order so that the stiffness
matrix K for the entire truss can be found
This is done by designating the rows & columns of the matrix by the 4 code
numbers used to identify the 2 global degrees of freedom that can occur at
each end of the member
The structure stiffness matrix will then have an order that will be equal to the
highest code number assigned to the truss since this rep the total no. of degree
of freedom for the structure
This method of assembling the member matrices to form the structure stiffness
matrix will now be demonstrated by numerical e.g.
This process is somewhat tedious when performed by hand but is rather easy
to program on computer
Determine the structure stiffness matrix for the 2 member truss shown in Figure
below. Take AE as constant.
4m
3m
2 kN
Step 1.
Determine degree of freedom (DOF).
a. Joint 1 and 3 – DOF = 0
b. Joint 2 – DOF = 2
Step 2.
Unknown displacement (DOF>0) need to code
numbered first. If the joints have two unknown
displacement, start with the horizontal
displacement first. The number followed by those at
joints having known displacement (DOF=0) .
Step 3.
Locate the origin of the global coordinate system
2 kN
(it can be located at any point).
The members are identified arbitrarily & arrows
are written along 2 members to identify the near
& far ends of each member
Step 4.
XF − XN YF − Y N
Determine λx and λy for every member. λx = ; λy =
L L
Member 1 Member 2
Member 1 Member 2
Step 6.
Corresponding elements of those 2 matrices are added algebraically to form the structure
stiffness matrix. As the member global stiffness matrices are generated, they are placed
directly into their respective element positions in the K matrix
The global force components Q acting on the truss can then be related to its
global displacements D using;
Q = KD ------ eqn. 14.17
Qk K 11 K 12 Du
= ------- eqn. 14.18
Qu K 21 K 22 Dk
Often Dk = 0 since the supports are not displaced, eqn. 14.19 becomes ;
Qk = K11(Du) ------- eqn. 14.21
The member forces can be determined using equation 14.13. Expanding this
equation yields;
DNx
q N AE 1 - 1 λx λy 0 0 DNy
Qu =
qF L - 1 1 0 0 λx λ y DFx
DFy
Since with qN = -qF for equilibrium; Only one need to be calculated, lets
determine qF .
D Nx
AE D Ny
qF =
L
[
- λx - λy λx λy ]
D Fx ------- eqn. 14.23
D Fy
Determine the forces in each member for the 2 member truss shown in Figure
below. Take AE as constant.
4m
3m
2 kN
Step 1.
Determine degree of freedom (DOF).
a. Joint 1 and 3 – DOF = 0
b. Joint 2 – DOF = 2
Step 2.
Unknown displacement (DOF>0) need to code
numbered first. If the joints have two unknown
displacement, start with the horizontal
displacement first. The number followed by those at
joints having known displacement (DOF=0) .
Step 3.
2 kN
Locate the origin of the global coordinate
system (it can be located at any point). By inspection, it is seen that the
The members are identified arbitrarily & known external displacement are
arrows are written along 2 members to D3=D4=D5=D6=0
identify the near & far ends of each member Also, the known external loads are
Q1=0, Q2=-2kN
Step 4.
XF − XN YF − Y N
Determine λx and λy for every member. λx = ; λy =
L L
Member 1 Member 2
Member 1 Member 2
Step 6.
Corresponding elements of those 2 matrices are added algebraically to form the structure
stiffness matrix. As the member global stiffness matrices are generated, they are placed
directly into their respective element positions in the K matrix
Step 7.
Writing equation, Q = KD.
Partition between Known and Unknown Forces / Displacement
Qk K 11 K 12 Du
= ------- eqn. 14.18
Qu K 21 K 22 Dk
Step 8.
From eqn. 14.19, Qk = K11(Du) + K12 (Dk), using Matrix multiplication, solve to determine
the value of D1 and D2.
Step 9.
From eqn. 14.20, Qu = K21(Du) + K22 (Dk), using Matrix multiplication, solve to determine
the value of Q3, Q4, Q5, and Q6,.
Q3 − 0.333 0 0 Q3 = − 1 .5 kN
D
Q
4 = AE 0 0
1 + 0 Q4 = 0 kN
Q5 − 0.072 − 0.096 D2 0 Q5 = 1 .5 kN
Q6 = 2 .0 kN
Q6 − 0.096 − 0.128 0
Step 10.
Calculate the member forces using eqn. 14.23.
D Nx
AE D
λy
Ny
qF =
L
[
- λx - λy λx ]
D Fx ------- eqn. 14.23
D Fy
For member 1,
λ x = 1, λ y = 0 , L = 3m ⇒ q1 = − 1 .5 kN
For member 2,
λ x = 0 .6, λ y = 0 .8, L = 5 m ⇒ q 2 = 2 .5 kN
Step 1. - The fixed end forces necessary to prevent movement of the nodes as
caused by temperature or fabrication are calculated
Step 2. - Equal but opposite forces are placed on the truss at the nodes & the
displacement of the nodes are calculated using the matrix analysis
Step 3. - The actual forces in the members & the reactions on the truss are
determined by superposing these 2 results
If the member is too short, then ∆L becomes negative & these forces will reverse.
In global coordinates, these forces are
(Q Nx )0 λx
(Q Ny )0 AE∆ L λ y
=
(Q Fx )0 L - λx ------- eqn. 14.26
(Q Fy )0 - λ y
CES 521 - INDETERMINATE STRUCTURES Slide 31 of 49 Topic 2 – Stiffness Method (Truss)
Thermal Changes & Fabrication Errors
(Pulau Pinang)
With the truss subjected to applied forces, temperature changes and fabrication
errors, the initial force-displacement relationship for the truss then becomes:
Member forces are obtained by superposition. The far end member forces
becomes;
D Nx
AE D Ny
qF =
L
- λx [ - λy λx λy ]
DFx
+ (q F )0 ------- eqn. 14.30
DFy
AE∆L
(q F )0 = − (Fabricati on Error)
L
(q F )0 = − AEα∆T (Thermal Changes)
Step 1.
Determine degree of freedom (DOF).
a. Joint 1, 3 and 4 – DOF = 0
b. Joint 2 – DOF = 2
Step 2.
Unknown displacement (DOF>0) need to code
numbered first. If the joints have two unknown
displacement, start with the horizontal
displacement first. The number followed by those at
joints having known displacement (DOF=0) .
Step 3.
Locate the origin of the global coordinate
system (it can be located at any point). By inspection, it is seen that the
The members are identified arbitrarily & known external displacement are
arrows are written along 2 members to D3=D4=D5=D6=D7=D8 = 0
identify the near & far ends of each member Also, the known external loads are
Q1 and Q2 =-0
Step 4.
XF − XN YF − Y N
Determine λx and λy for every member. λx = ; λy =
L L
Member 1
4−4 3−0
λx = = 0; λy = =1
3 3
Member 2
0−4 0−3
λx = = −0.8; λy = = −0.6
5 5
Member 3
4 −0 3−3
λx = = 1; λy = =0
4 4
Step 5.
Determine member Global Stiffness Matrix, k, for every member.
Member 1 Member 3
3 4 1 2 7 8 1 2
0 0 0 0 3 0.25 0 − 0.25 0 7
0 0.333 0 − 0.333 4 0 0 0 0 8
k1 = AE
0 0 1 k 3 = AE
0 0 − 0.25 0 0.25 0 1
0 − 0.333 0 0.333 2
0 0 0 0 2
Member 2
1 2 5 6
0.128 0.096 − 0.128 − 0.096 1
0.096 0.072 − 0.096 − 0.072 2
k 2 = AE
− 0.128 − 0.096 0.128 0.096 5
− 0.096 − 0.072 0.096 0.072 6
CES 521 - INDETERMINATE STRUCTURES Slide 37 of 49 Topic 2 – Stiffness Method (Truss)
Example 3
(Pulau Pinang)
Step 6.
Corresponding elements of those 3 matrices are added algebraically to form the structure
stiffness matrix. As the member global stiffness matrices are generated, they are placed
directly into their respective element positions in the K matrix
1 2 3 4 5 6 7 8
0.378 0.096 0 0 − 0.128 − 0.096 − 0.25 0 1
0.096 0.405 0 − 0.333 − 0.096 − 0.072 0 0 2
0 0 0 0 0 0 0 0 3
0 − 0.333 0 0.333 0 0 0 0 4
k = AE
− 0.128 − 0.096 0 0 0.128 0.096 0 0 5
− 0.096 − 0.072 0 0 0.096 0.072 0 0 6
− 0.025 0 0 0 0 0 0.25 0 7
0 0 0 0 0 0 0 0 8
Step 7.
Since the member is short, then ∆L = -0.01m. Applying equation 14.26 to member 2.
Step 8.
Writing equation, Q = KD + Q0.
Partition between Known and Unknown Forces / Displacement
Step 9.
From eqn. 14.28, Qk = K11(Du) + K12 (Dk) + (Qk)0, using Matrix multiplication, solve to
determine the value of D1 and D2.
D1 = − 0 .003704 m
0 0.378 0.096 D1 0 0.0016
= AE + + AE
0 0.096 0.405 D2 0 0.0012 D 2 = − 0 .002084 m
Step 10.
From eqn. 14.29, Qu = K21(Du) + K22 (Dk) + (Qu)0, using Matrix multiplication, solve to
determine the value of Q3 Q4, Q5, Q6, Q7 and Q8.
Q3 0 0 0 0
Q
4 0 − 0.333 0
0
Q5 − 0.128 − 0.096 D1 0 - 0.0016
= AE
+ + AE
Q
6 − 0.096 − 0.072 D
2 0 - 0.0012
Q − 0.025 0 0 0
7
8
Q 0 0
0 0
Step 11.
Calculate the member forces using eqn. 14.30.
D Nx
AE D Ny
qF =
L
- λx [ - λy λx λy ]
DFx
+ (q F )0 ------- eqn. 14.30
DFy
λ x = 0, λ y = 1, L = 3m, AE = 8(10 3 )kN
q 1 = − 5.56kN
λ x = − 0.8, λ y = − 0.6, L = 5m, AE = 8(10 3 )kN
q 2 = 9.26kN
CES 521 - INDETERMINATE STRUCTURES Slide 42 of 49 Topic 2 – Stiffness Method (Truss)
Example 4
(Pulau Pinang)
1. Determine the stiffness matrix K for the truss. Take A = 1000 mm2, E = 200 Gpa.
2. Determine the horizontal displacement of joint 1 and the force in member 2. Take A =
1000 mm2, E = 200 Gpa.
3. Determine the force in member 2 if its temperature is increased by 550C. Take A =
1000 mm2, E = 200 Gpa, α = 11.7 (10-6)/0C
2 kN
4m
4m 3m
A plane truss ABCD is subjected to uniformly distributed load and point loads as shown in
Figure 5. It carries points loads 30 kN, 40 kN and 35 kN at joints A, B and C respectively. The
upper chord ABC has uniformly distributed loads of 3 kN/m acted on it. The truss is pinned at
A and C, while D is on roller. Find all the support reactions
30 kN 40 kN 35 kN
3 kN/m
A C
B
BD = 3000 mm2
Figure 5
Review
Problems
A plane truss ABC is pinned at A and roller supported at B Figure P1(a). There is a horizontal
load of 50 kN acting at B of the truss. Compounding the external load, member AC is
fabricated 15 mm too long. The numbering of members and the nodal degrees of freedom
(DOF) are shown in Figure P1(b). Given modulus of elasticity, E = 200 kN/mm2 and constant
cross sectional area, A = 400 mm2.
Figure P2
Hibbeler R.C, (2012), “Structural Analysis, 8th Edition in S.I. Units”, Pearson,
Singapore
Hibbeler R.C, (2009), “Structural Analysis, 7th Edition in S.I. Units”, Pearson,
Singapore
Hibbeler R.C, (2006), “Structural Analysis, 6th Edition in S.I. Units”, Pearson,
Singapore
UiTM Structural Division (2003), “Basic Structural Analysis”, Cerdik
Publications Sdn. Bhd.
UiTM Structural Division (2003), “Basic Structural Mechanics”, Cerdik
Publications Sdn. Bhd.
Goh, C. H. et. al. (2012), “Problem and Solution on Stiffness Method for
Structures”, UiTM Press.