Transcript Slide 1

CHAPTER 3
DIRECT STIFFNESS METHOD FOR TRUSSES:
3.1 INTRODUCTION
In the previous chapter the procedure for obtaining the structure
stiffness matrix was discussed. The structure stiffness matrix was
established by the following equation:
[K] = [T]T [kc] [T]
-----------(2.21)
However if a large and complicated structure is to be analyzed
and if more force components are to be included for an element then
the size of composite stiffness matrix [kc] and deformation
transformation matrix [T] will be increased. Therefore the procedure
outlined in the preceding chapter for formation of structure stiffness
matrix appears to be inefficient, furthermore this procedure is not
suitable to automatic computation.
1
In this chapter an alternative procedure called the “direct stiffness
method” is introduced. This procedure provides the basis for most
computer programs to analyze structures. In this method each
individual element is treated as a structure and structure stiffness
matrix is obtained for this element using the relationship:
[K]m = [T]Tm[k]m [T]m
-----------(3.1)
Where
[K]m =
[T]m =
individual
[k]m =
Structure stiffness matrix of an individual element.
Deformation Transformation matrix of an
element.
Member Stiffness matrix of an individual element.
Total structure stiffness matrix can be obtained by superimposing the
structure stiffness matrices of the individual elements..
2
As all members of a truss are not in the same direction i.e.
inclination of the longitudinal axes of the elements varies, therefore
stiffness matrices are to be transformed from element coordinate
system to structure or global coordinate system. When the matrices for
all the truss elements have been formed then adding or combining
together the stiffness matrices of the individual elements can generate
the structure stiffness matrix [K] for the entire structure, because of
these considerations two systems of coordinates are required.
i) Local or member or element coordinate system:
In this coordinate system x-axis is collinear with the longitudinal
axis of the element or member. Element stiffness is calculated with
respect to this axis. This system is illustrated in figure 3.1
3
Element 'A'
Element 'B'
Local or member
or element axes
for element 'A'
Local or member
or element axes for
element B
A
Structure or
globel axes
B
Figure 3.1
4
ii)
Structure or global coordinate system:
A single coordinate system for the entire structure is chosen, with
respect to which stiffness of all elements must be written.
3.2 PROCEDURE FOR THE FORMATION OF TOTAL STRUCTURE
STIFFNESS MATRIX FOR AN ELEMENT USING DIRECT STIFFNESS
METHOD:
Following is the procedure for the formation of structure
stiffness matrix:
i)
Formation of the element stiffness matrix using equation 2.16.
k m
 1
AE 

L 
 1
 1


1
----------- (2.16)
ii)
Formation of the deformation transformation
matrix [T] for a single element:
[]m = [T]m []m
-----------(3.2)
5
where
[]m= Element or member deformation matrix.
[]m= Structure deformation matrix of an element or member
[T]m= Element or member deformation transformation matrix.
(iii)Formation of a structure stiffness matrix [K] or an element using the
relation
[K]m = [T]Tm[k]m [T]m -----------(3.1)
The above mentioned procedure is discussed in detail in the subsequent
discussion.
3.2.1
The formation of element stiffness matrix in local co-ordinates:
It has already been discussed in the previous chapter. However it
is to be noted that for horizontal members the structure stiffness matrix
and element stiffness matrix are identical because both member
coordinate systems and structure coordinate systems are identical. But
for inclined members deformation transformation matrices are to be
used because member coordinate system and structure coordinate
system are different; therefore their structure stiffness matrix and
6
element stiffness matrix will also be different.
3.2.2
The formation of deformation transformation matrix:
As the main difference between the previously discussed method and direct
stiffness method is the formation of the deformation transformation matrix.
In this article deformation transformation matrix for a single element will be
derived.
Before the formation of deformation transformation matrix following
conventions are to be established in order to identify joints, members,
element and structure deformations.
1) The member is assigned a direction. An arrow is written along the member,
with its head directed to the far end of the member.
2) i, j, k, and l, are the x and y structure deformations at near and far (tail
& head) ends of the member as shown in figure 3.2. These are positive in the
right and upward direction.
3) r and s are the element deformations at near and far (tail & head) ends of
the member as shown in figure 3.2.
4)The member axis (x-axis of member coordinate system) makes an angle x,
with the x-axis of the structure coordinate system as shown in figure 3.2.
5)The member axis (x-axis of member coordinate system) makes an angle y,
with the y-axis of the structure co-ordinate system as can be seen from figure
3.2.
6) The cosines of these angles are used in subsequent discussion Letters l and m7
represent these respectively.
l
s
Far
end
y
Near
end
j
r
Member
direction
x
i
k
j
r
Near
end
(a)
i
y

Member
direction
k
Far
end
Figure 3.2
(b)
s
x
8
 = Cos x and m = Cos y
 = Cos x =
X 2  X1

L
m = Cos y =
X 2  X1
 X 2  X 1 2  Y2  Y1 2
Y2  Y1

L
Y2  Y1
 X 2  X 1 2  Y2  Y1 2
The algebraic signs of θ‘s will be automatically accounted for the members
which are oriented in other quadrants of X-Y plan.
In order to form deformation transformation matrix, once again consider
the member of a truss shown in figure 3.1.
9
s=0
s=0

r

=1
j


r
c
i
l

x

os
j
x
=1
i
i
y
i
r
j
(a)
x

s
o
j c

y=
m

r

(b)
10
s
y=
m
s
y
k

s
o
c
k


s
k
1
y
x=
x
=1
l
1
s
s

k=1
co
1
x
r=0
r=0
(c)
(d)
Figure 3.3
11
Following four cases are considered to form deformation transformation matrix.
Case-1: Introduction of horizontal deformation to the structure i = 1, while far
end of the member is hinged (restrained against movement). From the geometry
of figure 3.3(a)
r = i Cos x = 1. Cos x = Cos x = l
r =l
---------- (3.3)
s = 0
---------- (3.4)
Case-2: Introduction of vertical deformation to the structure j = 1, while near
end of the member is hinged (restrained against movement). From the geometry of
figure 3.3 (b)
r = j Sin x = j Cos y = 1. Cos y = m
r = m
---------- (3.5)
s = 0
---------- (3.6)
Case-3: Introduction of horizontal deformation to the structure k = 1, while
near end of the member is hinged (restrained against movement).From the
geometry of figure 3.3(c)
s = k Cos x = 1. Cos x = Cos x = l
r = 0
---------- (3.7)
s = l
---------- (3.8)
12
Case-4: Introduction of vertical deformation to the structure L = 1, while
near end is hinged
From the geometry of figure 3.3(d)
L = 1
r = 0
---------- (3.9)
s =  L
Cos y = m
---------- (3.10)
Effect of four structure deformations and two member deformations can be
written as
r = i Cos x + j Cos y + k .0 +  L .0
s = i .0 + j .0 + k Cos x +  L Cos y
r = l.i + m.j + 0.k + 0.  L
s = 0.i + 0.j + l.k + m.  l
----------- (3.11)
----------- (3.12)
Arranging equations 3.11 and 3.12 in matrix from
 i 
 r    m 0 0    j  ----------- (3.13)
 s   0 0  m    k 
 
 L 
13
Comparing this equation with the following equation
[]m = [T]m []m
---------------(3.2)
After comparing equation (3.14) and (3.2) following deformation
transformation matrix is obtained.
T m

m
0
0


0 0  m 
---------- (3.14)
This matrix [T]m transforms four structure deformation (i, j, k, L )
into two element deformation (r and s).
3.2.3
Formation of structure stiffness matrix:
Structure Stiffness matrix of an individual member is first to be
transformed from member to structure coordinates. This can be done by
using equation 3.1.
[K] = [T]Tm [k]m [T]m
--------(3.1)
l
T m  
0
m
0
0
l
0
m 
----- (3.15)
14
k m
K m
AE  1  1

1
L  1
------(2.16)
i  0
i j k L


j m 0   1  1  m 0 0 
 



k  0    1 1  0 0  m


L  0 m
K m
i  0
i
j
k
L


j m 0  AE  
m    m
 



L
k 0  

m 
   m


L  0 m
i
K m
j
k
L
 2
m   2  m i


2
2 j
m
 m  m 
AE  m



L    2  m  2
m  k


2
2


m

m

m
m

 L
--------------------------------(3.16)
15
3.2.4 ASSEMBLING OF THE INDIVIDUAL STRUCTURE ELEMENT
STIFFNESS MATRICES TO FORM TOTAL STRUCTURE
STIFFNESS MATRIX
Combining the stiffness matrix of the individual members can generate
the stiffness matrix of a structure However the combining process
should be carried out by identifying the truss joint so that matrix
elements associated at particular member stiffness matrices are
combined. The procedure of formation of structure stiffness matrix is
as follows:
Step 1: Form the individual element stiffness matrices for each member.
Step 2: Form a square matrix, whose order should be equal to that of
structure deformations.
Step 3: Place the elements of each individual element stiffness matrix
framed into structure in the corresponding rows and columns of
structure stiffness matrix of step-2.
Step 4: If more than one element are to be placed in the same location of
the structure stiffness matrix then those elements will be added.
16
Following figure shows the above mentioned process.
Structure
stiffness matrix
1
2
3
Structure
stiffness
matrix for
member-1
1
2
2
3
Structure
stiffness
matrix for
member-2
17
1
 k1
K 1  
  k1
2
 k1  1
 k 1  2
K 2
2
 k 2

 k 2
3
 k2  2
 k 2  3
K   K 1  K 2
1
 k1
K    k1
 0
2
 k1
k1  k 2
 k2
3
0 1
 k 2  2
k 2  3
In the above matrix the element in the second row and
second column is k1+k2 where k1 is for member 1 and k2
is for member 2.This is because both k1 and k2 have
structure coordinates 2-2.
18
Illustrative Examples Regarding the Formation of [K] Matrix:
Example 3.1:
following
Form the structure stiffness matrix for the
truss by direct stiffness method.
12
11.5
2
2
1
(x2,y2) = (L/2,2L/3)
2L/3
1
(x1,y1) = (L/2,2L/3)
1
2
6
5
(x1,y1) = (0,0)
6
(a)
L/2
4
(c)
2
3
(x2,y2) = (L,0)
1
4
1
3
5
(x2,y2) = (0,0)
3
(x1,y1) = (L,0)
(Free body diagram indicating structure forces
and coordinates)
L/2
6
5
(b)
3
2
4
3
(Structure forces and deformations)
19
Member
Length Near
End
X 1 Y1
1
5L/6
0
2
5L/6
L/
Far End
X2, Y2
0
i
j k L l=Cosx =
X2-X1/
Length
L/2 2L/3 5 6 1 2
m=Cosy
Y2-Y1/
Length
0.6
0.8
2L/3
L
0
1 2 3 4
0.6
-0.8
0
0
0
3 4 5 6
-1
0
2
3
L
L
According to eq. (3.16) we get
i
l

AE lm
K   
L l
 lm
2
m
2
j
lm
m2
 lm
 m2
k
l
 l 2  lm  i
 lm  m 2 j

l2
lm k
lm
m2 
l
20
5
6
1
2
0.432
0.576  0.432  0.576 5

6
AE 0.576
0
.
768

0
.
576

0
.
768
K 1   0.432  0.576 0.432 0.576  1
L
 0.576  0.768 0.576 0.768  2
1
2
3
4
0.432  0.576  0.432 0.576 1

2
AE  0.576 0.768
0
.
576

0
.
768
K 2   0.432 0.576 0.432  0.576 3
L
 0.576  0.768  0.576 0.768  4
3
4
5
6
1.00 0  1.00 0 3

0
K 3  AE 
L  1.00
 0
0
0
0
1.00
0
0

0 4
0 5
0 6
21
1
0.864
0.000


AE  0.432
K  
L  0.576
 0.432
 0.576
2
0.000
1.536
0.576
3
4
5
6
 0.432 0.576  0.432  0.576 1

0.576  0.768  0.576  0.768 2
1.432
 0.576  1.000
 0.768  0.576
 0.576  1.000
0.768
0.000
0.000
1.432
 0.768
0.000
0.576
0.000
0.000  3
0.000  4
0.576  5
0.768  6
22
Example 3.2:
Form the structure stiffness matrix of the following truss.
P
2
3
1
2
(L,L)
1
L
2
(L,L)
1
1
L
3
8
2

4
Structure to be analysed
y
7
(0,0)
1
x
(L,L)
2
5
L
2
3
1
(2L,0)
(L,0)
Free body diagram showing structure forces and
coordinates
7
4
6
8
6
5
3
Structure forces and deformations
23
Member Near End
Far End Length =Cosx m=Cosy
i
j
k L
X1
Y1
X2
Y2
1
0
0
L
L
2L
0.707
0.707
7
8
1
2
2
L
L
L
0
L
0
-1
1
2
5
6
3
L
L
2L
0
2L
0.707
-0.707
1
2
3
4

AE  lm
L  l
 lm
l2
K m

2
 lm 
lm
l 2
m2
 lm  m 2 
 lm
 m2
l2
lm
lm
m2


So
24
7
0.5
K 1  AE 0.5
2L .5
.5
1 2
0 0
AE 0 1
K 2 
L 0 0
0 1
5 6
0 0 1

0 1 2
0 0 5
0 1 6
8
0.5
0.5
.5
.5
1
.5
.5
0.5
0.5
2
.5 7

.5 8
0.5 1
0.5 2
1
K 3 
AE
2L
 .5
 .5
 .5.5

2
3
 .5  .5
.5
.5
.5
.5
 .5  .5
4

 .5

 .5

.5
1
.5
2
3
4
Now the matrix can be written as [k]2:
1
2
0
0
K 2  AE 0 1.414
0
2L 0
0 1.414
5
6
0
0 1

0 1.414 2
0
0 5
0 1.414 6
So [k] = [k]1 + [k]2 + [k]3
7
8
1
2
 .5 .5 .5 .5 7
K 1  AE .5 .5 .5 .5 8
2L .5 .5 .5 .5 1
.5 .5 .5 .5 2
25
1
2
3
4 5
6
7
0 .5 0.5 0
0  .5
 1
 0 2.414 .5 .5 0 1.414 .5
.5
.5 .5 0
0
0
.5
.5 .5
.5 0
0
0
AE  .5
K  
0
0
0 0
0
0
2L  0
 0 1.414 0 0 0 1.414 0
.5
0
0 0
0
.5
.5
.5
.5
0
0
0
0
.5
8
.5
1
.5
2
0 3
0 4
0 5
0 6
7
.5
.5
8
26
3.4 Force transformation matrix:
The axial force (wm) in the members of a truss can be
found by using the relationship between member forces
and member deformations [equation (2.15)] and between
element and structure deformations [equation (3.1)]
[w]m= [k]m []m ------------ (2.15)
[]m = [T]m []m ------------ (3.3)
Substituting value of []m from equation 3.1 into equation
2.15
[w]m = [k]m [T]m []m ------------ (3.18)
Substituting value of [k]m and [T]m from equation 2.16 and
3.15 respectively
 i 
 
AE  1 1  
m



m

wm  1 1     m  m   j 

  k 
L 
 
 L
27
 i 
 
 wr  AE  
m



m

 j
 ws   L    m 
m   
 
k
 L 
In above Equations, wr= Force at near end
See figure below:
ws = Force at far end
As wr = -ws so:
AE
ws    l
L
m l
 i 
 
m  j 
 k 
 
 l 
-------------- (3.19)
Sign conventions
wr
ws
28
If
wr is positive then member is in compression.
ws is negative then member is in compression.
wr is negative then member is in tension.
ws is positive then member is in tension.
3.5 ANALYSIS OF TRUSSES USING DIRECT STIFFNESS METHOD:
Basing on the derivations in the preceding sections of this chapter
a truss can be completely analyzed. The analysis comprises of
determining.
i)
Joint deformations.
ii) Support reactions.
iii) Internal member forces.
As the first step in the analysis is the determination of unknown joint
deformation. Using the equations can do this.
[W] = [K] []
The matrices [W], [K] and [] can be divided in submatrices in the
following form
Wk 
Wu 
 K11
 K 21
K12 
K 22 
 u 
 k 
----------- (3.20)
29
Where
Wk =
Wu =
u =
k =
conditions.
Known values of loads at joints.
Unknown support reaction.
Unknown joint deformation.
Known deformations, generally zero due to support
K11, K12, K21, K22 are the sub-matrices of [K]
Expanding equation 3.20
Wk = K11 u + K12 k ----------- (3.21)
Wu = K21 u + K22 k ----------- (3.22)
If the supports do not move, then k = 0 therefore equation 3.21
& 3.22 can be written as
Wk = K11 u
----------- (3.23)
Wu = K21 u
----------- (3.24)
By pre-multiplying equation 3.23 by [K11]-1 following equation is obtained
[K11]-1 [WK]=[K11]-1 [K11][U]
[u] = [K11]-1 [Wk]
----------- (3.25)
Substituting value u from equation 3.25 into equation 3.24
[Wu] = [K21] [K11]-1 [Wk] ----------- (3.26)
30
Using equation 3.25, 3.26 and 3.19, joint deformations, support reactions
and internal member forces can be determined respectively.
As this method does not depend upon degree of indeterminacy so it can
be used for determinate as well as indeterminate structures.
Using the basic concepts as discussed in the previous pages, following
are the necessary steps for the analysis of the truss using stiffness
method.
1-Identify the separate elements of the structure numerically and specify
near end and far end of the member by directing an arrow along the
length of the member with head directed to the far end as shown in the
fig. 3.2
2-Establish the x,y structure co-ordinate system. Origin be located at one
of the joints. Identify all nodal co-ordinates by numbers and specify
two different numbers for each joint (one for x and one for y). First
number the joints with unknown displacements.
3-Form structure stiffness matrix for each element using equation 3.16.
4-Form the total structure stiffness matrix by superposition of the element
stiffness matrices.
5-Get values of unknown displacements using equation 3.25.
6-Determine support reaction using equation 3.26.
7-Compute element or member forces using equation 3.19.
31
3.6- Illustrative Examples Regarding Complete Analysis of Trusses:
Example 3.3: Solve truss in example 3.1 to find member forces.
12
11.5
2
2
1
(x2,y2) = (L/2,2L/3)
2L/3
1
(x1,y1) = (L/2,2L/3)
1
2
6
5
(x1,y1) = (0,0)
6
(a)
L/2
4
(c)
2
3
(x2,y2) = (L,0)
1
4
1
3
5
(x2,y2) = (0,0)
3
(x1,y1) = (L,0)
(Free body diagram indicating structure forces and
coordinates)
(b)
4
6
5
2
3
3
(Structure forces and deformations)
32
The [K] matrix for the truss as formed in example 3.1 is:
1
 0.864
0.000

 0.432
K  AE 
L  0.576
 0.432
 0.576
2
0.000
1.536
3
4
5
6
 0.432 0.576  0.432  0.576 1

0.576  0.768  0.576  0.768 2
0.576
1.432
 0.576 1.000
 0.768  0.576
 0.576 1.000
0.768
0.000
0.000
1.432
 0.768
0.000
0.576
0.000

0.000  3
0.000  4
0.576  5
0.768  6
Using the relation [W] = [K] [] we get
Wk   K11
W   K
 u   21
K12 
K 22 
 u 
 
 k
33
1
2
3
4
5
6
W1 
0.864
0.000  0.432 0.576  0.432  0.576  1 

W 
 
0
.
000
1
.
536
0
.
576

0
.
768

0
.
576

0
.
768

  2
 2
W3  AE  0.432 0.576
1.432  0.576  1.000 0.000    3 
 
 
W
L
 0.576  0.768  0.576 0.768 0.000 0.000   4 
 4
W5 
 0.432  0.576  1.000 0.000
1.432
0.576    5 

 
 
W
 6 
 0.576  0.768 0.000 0.000 0.576 0.768   6 
From equation (3.25) [u] = [K11]-1 [Wk]
1
 1 
L 0.864 0.000 W1 
  
0.000 1.536 W 
AE

  2
 2
 
 1


 2 


L 1.157 0.000 11.5


 
0
.
000
0
.
651

12

AE 
 

 1  L  13.31 
   AE 7.812 


 2
34
Using the relation:
[Wu] = [K21] [u]
W3   0.432 0.576 
W  
 13.310
0
.
576

0
.
768

 4  

W5   0.432  0.576  7.812
  

W

0
.
576

0
.
768
 6  

W3  10.25
W  

13.66
 4  

W5   1.25 
  

W
 6   1.66 
Now the member forces can be found using the relation:
35
 i 
 
AE
wr    m    m  j 
L
 k 
 
 L
 0.00 
 0.00 
6 AE
 L  2.09
w1  
0.6 0.8 0.6 0.8 
 13.31  AE
5L


 7.812
kips
 13.31 
 7.812
6 AE
 L  17.09
w2  
0.6 0.8 0.6 0.8 
 0.00  AE
5L


 0.00 
0.00
0.00
AE
w3   1.0 0.0 1.0 0.0   L  0.00
0.00 AE
L


0.00
kips
kips
36
Final Results & sketch
12k
11.5k
kip
2.0
9k
09
ips
17.
13.66k
1.66k
1.25k
0 kips
10.25k
37
Example 3.4:
Solve the truss in example 3.2 to find member forces.
P
2
3
1
2
(L,L)
1
L
2
(L,L)
1
1
2
y
7
(0,0)
L
3
8
1
x
4
(L,L)
(Structure to be analysed)
2
3
2
5
L
1
(2L,0)
(L,0)
(Free body diagram showing structure forces and
coordinates)
7
4
6
8
6
5
3
(Structure forces and deformations)
38
The [K] matrix for the truss as formed in example 3.2 is as follows:
1
1
0
  .5
K   AE  .5
2L  0
0
  .5
.5
2
3
4
5
0
0
0
0
0
0
1.414
0
0
0
1.414
0
.5
.5
0
0
0
0
0
0
0
0
.5
.5
0
2.414
.5
.5
6
7
8
.5 0.5 0
0
.5 .5 1
.5 .5 0 1.414 .5 .5 2
.5 .5 0
0
0
0 3
.5 .5 0
0
0
0 4
5
0 6
.5  7
.5 
8
0
Now as W = k
39
1
2
3
4
1
0
 .5 0.5
W1 

W 
 2
 0 2.414 .5  .5
W3 
 .5
.5
.5  .5

 
 .5
 .5 .5
W4   AE  .5
W5 
0
0
0
2L  0
 
W6 
 0  1.414 0 0
W 
 .5  .5 0 0
 7
W8 
 .5  .5 0 0
5
0
0
0
0
0
0
0
0
6
7
8
0
 .5  .5

 1.414  .5  .5
0
0
0 
0
0
0 
0
0
0 
1.414
0
0 
0
.5
.5 
0
.5
.5 
 1 
 
 2
 3 
 
 4 
 5 
 
 6 
 
 7
  8 
k12   u 
k 22   k 
Wk   k11
W   k
 u   21
As k = 0 or 3, 4, 5, 6, 7,8, all are zero due to supports and using
[Wk] = [k11] [u] we get
0 
W1 
AE 1

W 
0 2.414
2
L
 2


Which gives
1 
 
 2
40
W1 = 0
W2 = P
Also using [Wu] = [k21] [u] we have
0   1 
0
AE 1
 P 
0 2.414  
2
L
 

 2 
1
0  0
 1 
2 L 1
  
0 2.414  P
AE

  
 2
Or
2L
1 2.414

.

AE 2.414  0
0  0 
 P
1
 
Or
41
1 
2L

 
 2  2.414 AE
0
 P
 
Or
2 
2 LP
2.414 AE
Or
 2  0.5858
Now
PL
AE
Wu = [K21] [u]
.5 
W3 
.5
W 
 .5


.
5
4
 


W5 
0  1 
AE  0
 

 
W
0

1414
.
2L 
 6
  2 
W7 
.5
.5 
 


W

.
5

.
5
 8 


W3 
W 
 4
W5 
1



W
 6  2.414
W7 


W8 
.5
 .5

 0

 0
.5

.5

.5 
0  0 

1.414   P 
.5 

.5 
.5
42
W3 
 0.5 P 
W 
 0.5 P 
4
 


W5 
1  0 
 


W

1414
.
P
2
.
414
 6


W7 
 0.5 P 
 


W8 
 0.5 P 
W3 = 0.207 P
W4 = -0.207 P
W5 = 0
W6 = -0.5858 P
W7 = -0.207 P
W8 = -0.207 P
Now as
 wr  AE  l
 w   L  l

 s
m l
m l
So
ws 
AE
 l
L
m l
i 
 m  j 
m   k 
 
l 
 i  7
  8
m  j 
 k  1
 
l  2
43
So, for member # 1:
AE
 0.707  0.707 0.707 0.707
ws  w1 
L
Or
ws= 0.707 x 0.5858 P
w1= 0.414 P

0


0


0

0.5858P







L 
AE 

Similarly, for member # 2:
ws  w2 
Or
Or
AE
0
L
1 0
0

1
PL 

.
5858
2


AE
 1

5
0


0

 6
w2 = 1 x 0.5858 P
w2 = 0.5858 P
44
For member # 3:
 i  1
  2
AE
 l  m l m  j 
Ws 
 k  3
L
 
l  4
0


 .5858PL 
AE
 .707 .707 .707  .707  AE 
w3 
L
0




0
Or
w3= 0.707 x 0.5858 P
w3= 0.414 P
45
P
0.
0.207P
0.
41
4P
0.586P
4P
1
4
0.207P
0.207P
0.59P
0.207P
w1 = 0.414 P
w2 = 0.586 P
w3 = 0.414 P
46
Example 3.5
Analyze the truss shown in the figure.
8'
W6,6
W3,3
W5 ,5
3
4
W4,4
5
2
6'
W7,7
8'
W8,8
1
5k
W1,1
W2,2
Structure forces and
displacements
2k
Truss to be analysed
A and E are constant for each member.
47
Member
Length
l
m
i
j
k
l
1
10
-0.6
-0.8
7
8
1
2
2
11.34
-0.7071
0.7071
1
2
3
4
3
8
1
0
3
4
5
6
4
6
1
0
5
6
7
8
5
8
0
1
1
2
5
6
Using the properties given in the above table we can find the structures
stiffness matrices for each element as follows.
48
w5 , 5
w6 , 6
3
w7 , 7
w8 , 8
w10 , 10
w4 , 4
4
w1 , 1
5
1
2
w9 , 9
w3 , 3
w2 , 2
Element forces and displacements
49
1
8
7
2
0.048  0.036  0.048 7
 0.036
 0.048
0.064  0.048  0.064 8

k1  
0.048  1
  0.036  0.048 0.036


064
.
0
048
.
0
064
.
0

048
.
0

2

4
3
2
1
 0.044  0.044
  0.044 0.044
k2  
  0.044 0.044

 0.044  0.044
5
4
3
 0.125

0

k3 
  0.125

0

0
0
0
0
 0.044
0.044
0.044
 0.044
0.044  1
 0.044 2

 0.044 3

0.044  4
6
0.125 0 3
0 4
0

0.125 0 5

0 6
0
50
5
6
 0.167

0

k4 
  0.167

0

3
0
0
k5  
0

0
0
0
0
0
4
0
0.125
0
7
8
 0.167 0
0
0

0.167 0

0
0
7
0
0
0
 0.125 0
5
6
7
8
8
0

 0.125

0


0.125 
3
4
7
8
Using equation [K] = [K1]+[K2]+[K3]+[K4]+[K5] we get the structure
stiffness matrix of (8x8) dimensions. Partitioning this matrix with respect
to known and unknown deformations we get [K11] and [K12] portions as
follows.
51
1
1
2
 0.004
 0.044

 0
K 12   
 0
 0.036

  0.048
 0.08 0.004 1
K11   
2
0
.
004
0
.
233


2
0.044  3
 0.044 4
0 5

 0.125 6
 0.048 7

 0.064 8
Using equation [Du]=[K11]–1[Wk] we get the
1
 1   0.08 0.004  5   62.983
   0.004 0.233  2   9.665
   

 2 
Putting the above value in equation [Wu]=[K21][Du]
52
W1   0.004
W   0.044
 2 
W3   0
 
W4   0
W5   0.036
  
W6   0.048
0.044 
 0.044
0 

 0.125
 0.048

 0.064
 3.197
 3.197 


 62.983  0 
 9.665   1.208 

 

  1.803


 2.405
Element forces can be calculated using equation as follows.
53
 0 
 0 
1
  3.006 kip
0.6 0.8  0.6  0.8 
w1 

62.983
10


 9.665
 62.983
 9.665
1
  4.540 kip
0.707  0.707  0.707 0.707 
w2 
 0 
11.314


 0 
1.208 K
3.197 K
0K
 62.983
 9.665
1
  1.208 kip
w5  0  1 0 1 
 0 
8


 0 
4k
ip
s
-3.
00
6k
4.5
1.803 K
ips
3.197 K
0 
0 
1
w 4   1 0 1 0    0 kip
0 
6
 
0 
2.405 K
1.208 kips
0 
0 
1
w3   1 0 1 0    0 kip
0 
8
 
0 
5K
2K
54