Chapter 18 Matrix Analysis of Beams and Frames by The Direct Stiffness Method
Chapter 18 Matrix Analysis of Beams and Frames by The Direct Stiffness Method
Chapter 18 Matrix Analysis of Beams and Frames by The Direct Stiffness Method
HH
+ H
~ , u
18:1 Introducti.on
In Chapter 17 we discussed the analysis of trusses using the direct stiff
ness method. In this chapter we extend the method to structures in which
loads. may be applied to joints as well as to members between joints, and
induce both axial forces and shears and moments. Whereas in the case of.
trusses we had to consider only joint displacements as unknowns in set
tingup the equilibrium equations, for frarneswemtlsf adrljoint rotations.
Consequently, a total of three' equations of equilibrium, two for forees and
one for mQm~t, can be written for each joint in a plane frame.
Even thou,'gh the analysis of a plane frame using the direct stiffness
method involves three displacement components per joint (8, Il", Il), we
can often reduce the number of equations to be solved by neglecting the
change in length of the members. In typical beams or frames, this sim
plification introduces little error in the results.
In the analysis of any structure using the stiffness method, the value
of any quantity (for example, shear, moment, or displacement) is obtained
from the sum of two parts. The first part is obtained from the analysis of
a restrained structure in which all the joints are restrained against move
ment. The moments induced at the ends of each member are fixedend
moments. This procedure is similar to that used in the moment distribu- '.
tion method in Chapter 13. After the net restraining forces are computed
and the signs reversed at each joirtt, these restraining forces are applied
to the original structure in the second part of the analysis to determine
the effect induced by joint displacements,
The superposition of forces and. displacements from two parts can be
explained using as an example theframe in Figure I8.la. This frame is
composed, of two members connected by a rigid joint at B. Under the
..."' .....
684
Chapter 18
I
/ !'
t 1
t,
I
/'
t
Ab
\1
\1
MD=Mb+M;)
Milc+MBC
~
.. '.~
Mb
MilCB
Milc
MeB
MBA = MBA
I
AD
MBA
r:.======~::::;::!f
+
il
MAB
/
(a)
(c)
(b)
loading shown, the structure will defonn and develop shears, moments,
and axial loads in both members. Because of the changes in length
induced by the axial forces, joint B will experience, in addition to a rota
tion 8B , small displacements in the x and y directions. Since these dis
placements are small and do not appreciably affect the member forces, we
neglect thePl .. With this simplification we can analyze the frame as having
only one degree of kinematic indeterminacy (Le., the rotation of joint B).
. In the first part of the analysis, which we. designate as the restrained
condition, .we introduce a rotational restraint (an imaginary c1amp)at
joint R (see Fig. 18.1h). The annition of the clamp transfonns the stmc
ture into two fixed-end beams. The analysis of these beams can be read
ily carried out using tables (e.g., see Table 12.5). The deflected shape and
the corresponding moment diagrams (directly under the sketch of the
frame) are shown in Figure 18.lb. Forces and displacements associated
with this case are superscripted with a prime.
Since the counterclockwise moment M applied by the clamp at B
does not exist in the original structure, we must eliminate its effect. We
'a-,...... _
Section 18.2
do this in the second part of the analysis by solving for the rQtation OR of
joint B produced by an applied moment that is equal in magnitude but
opposite in sense to the moment applied by the clamp. The moments and
displacements in the members for the second part of analysis are super
scripted with a double prime, as shown in Figure 18.1c. The final results,
shown in Figure 18.la, follow by direct superposition of the cases in Fig
ure 18.lb and c.
We note that not only are the final moments obtained by adding the
values in the restrained case to those produced by the joint rotation OR'
but also any other force or displacement can be obtained in the same
manner. For example, the deflection directly under the load IlD equals the
sum of the corresponding deflections at D in Figure 18.1 band c, that is,
IlD = Ilh
+ 1l'D
..... . .................................................
:1i;2J)lrSt;~ct~~~ stiff~;~;M~t~i~
685
686
Chapter 18
18.3
~.iI
In this section we derive the member stiffness matrix for an individual flex
ural element using only joint rotations as degrees of freedom. The 2 X 2
matrix that relates moments and rotations at the ends of the member is
important beca,Qse it can be used directly in the solution of many practi
cal problems, such as continuous beams and braced frames where joint
translations are prevented. Furthermore, it is a basic item in the deriva
tion of the more general 4 X 4 member stiffness matrix to be presented
.in Section 18.4 .
. Figure 18.2 shows beam of length L with cnd moments M j and Mj
. As a sign convention the end rotations OJ and OJ are positive when clock
wise and negative when counterclockwise. Similarly, clockwise end
moments are also positive, and counterclockwise moments are negative.
To highlight the fact that the derivation to follow is independent of the
member orientation, the axis of the element is drawn with an arbitrary
inclination a.
.
In matrix notation, the relationship between the end moments and the
resulting end rotations can be written as
chord
(18.2)
Figure 18.2: End rotations produced by member
end moments.
........... -
.......... -
... .....
.;.
Section 18.3
and
2EI
MJ = -(0
L
'
687
+ 20.)
J
(18.4)
[~;]= 2:/[~ ~J
[:;] ,
(185)
k is
- = 2E/[2
k
IJ
(18.6)
1 2
EXAMPLE 18.1
Using the direct stiffness method, analyze the frame shown in Figure 18.3a.
The change in length of the members may be neglected. The frame con
1
L
W=
Solution
trates the positive direction (clockwise) selected for the rotational degree
of freedom at joint B.
1----- L - - - i
(a)
....... ..- -
..-.-.:i.. .......
688
Chapter 18
0.125 PL
degree of freedom
Ih~_
0.125 PL
E
18
,
,I
A
(b)
_
1
I
I
.)0.25PL
0.125 PL
~)0.125PL
~
0.25 PL
(d)
(c)
into two fixed-end beams (Fig. 18.3c). The fixed-end moments (see Fig.
~~2 ;; _
3: (~~) ;;
,PL
MBA = -MAB =
PL
4
(18.7)
... '(18.8)'
M' = _PL
Be
8
,
MeB
(18.9)
,PL
= -MBc = 8
(18.10)
Figure 18.3c shows the fixed-end moments and the deflected shape
of the restrained frame. To illustrate the calculation of the restraining
moment Mt> a free-body diagram of joint B.is also shown in the lower
right comer of Figure I8.3c. For clarity, shears acting on the joint are
omitted. From the requirement of rotational eqUilibrium of the joint
('2:MB = 0) we obtain
.
PL
PL
M
+-+
48
(18.11)
S
In this I-degree of freedom problem, the value of M t with its sign
reversed is the only element in the restraining force vector F (see Eq.
Section 18.3
689
18.1). Figure 18.3d shows the moment diagrams for the members in the
restrained structure.
[MAB]
MBA
2El[2
L 1
2El[2
L 1
'i
+
~I~J= +
(e)
[4EIJ
L,-:-=~~-..;::-...,......,""'c
2EI
Kl1
4El
L
=-+
4EI
8EI
=L
L
'I
\ .
a~J~ [Ul]
and
[M BC ]
MeB
4EI I
L I
2EI
(f)
(18.12)
In this problem the value given by Equation 18.12 is the only element of
the stiffness matrix K. The moment diagrams for the members corre
sponding to the condition OB = 1 rad are shown in Figure 18.3/
Ka =F
8EL 0 __ PL
L
B 8
(18.13)
........... -
..
690
Chapter 18
MAB -
/I
2EI Q
L UB
4EI
__
2EI
L
MeB
I/
(18.15)
PL
PL
32
-16
(18.16)
PL
16
(18.17)
__
PL
(18.18)
= - 32
The double prime indicates that these moments are associated with the
joint displacement condition.
= MAE + MAE = - PL
4
I
MAB
/I
(PL)
- 32
9PL
=-3'2
3PL
=-
16
0.142PL
r""
+ M Bc =
+ MCB
M Bc
M CB
= MCB
M Bc
."
O.S3L
~""'-'--">.
9PL
/I
PL
8
PL
8
_ 3PL
+ (PL)
-16 --16
(PL)'"
+ -
32
3PL
32
32
(g)
'"
';
Section 18.3
691
EXAMPLE 18.2
,Solution
inspection of the structure indicates that the degree of kinematic
'indeterminacy is 3. The positive directions selected for the 3 degrees of
freedom (rotations at joints B, C, and D) are shown with curved arrows
in Figure 18.4b.
! An
The fixed-end
moments induced in the restrained structure by the applied loads are cal
culated using the formulas in Figure 12.5. Figure 18.4c shows the moment
diagram for the restrained condition and the free-body diagrams of the
joints that are used to calculate the forces in the restraints. Considering
moment eqUilibrium, we compute the restraining moments as follows:
Joint B:
Joint C:
Joint D:
+ 100
M[
150
M2
\.-150+M3
M3
Ml
-100
+ M2 +
clarity);
= -100 kipft
= -'-50 kipft .
150 kipft
20 kips
100]
.
50 kipft
[
-150
4.5 kips/ft
~IIW
,,___ / . . . c~
_/~~~
.
. _. . . . .
(18.19)
J?
"iW~
150
150
150
20'
I'C,M2 1'[/<13
EI =constant
(~)("
(a)
100
150
(c)
(b)
.e'. . . . . _
150
~)
c:=::::J
20'--1--
150
692
Chapter 18
The forces
at the ends of the members resulting from the introduction of unit dis
placements at each one of the degrees of freedom are shown in Figure
18.4d to f. The elements of the structure stiffness matrix are readily cal
culated from the free-body diagrams of the joints. Summing moments,
we calculate from Figure 18.4d:
-0.2EI
O.lEI
-0.05EI
+ Kll = 0
and
+ K21
:::::
and
K21
= 0.05E1
K31 :::::
and
K31
81 =1 rad
/C~K21
(~
0.05EI.
(d)
a
0.2El + K22 ::::: a
and
K12 = 0.05EI
and
K22
= 0.3EI
+ K32 = 0
and
K32
= O.lEI
-0.05EI
-O.1El -
-O.lEI
+ K12
I'C~K22
/"B-<!2
Figure 18.4: (e) Stiffness coefficients produced
by a unit rotation of joint C with joints Band D
restrained;
~)
(~)(~
O.1EI 0.2El O.lEl
0.051
(e)
I'D~K32
Section 18.3
693
= 0
and
Kl3
-0. lEI
+ K23
= 0
and
K 23
-0.2EI
+ K33
and
= O.lEI
K33 = 0.2EI
I lowing
0.3
EI [ ~.05
0.05
0.3
0.1
~.1]
(18.20)
0.2
(n
Figure 18.4: .
Stiffness coefficie~ts produced
by a unit rotation of joint D with joints Band C
restrained;
0.lEIO.2EI
(f)
1
]
100]
-1~~
(UUl)
1[
25.8.6]
448.3
_
8f)I]
2
f)3
EI -974.1
(18.22)
The
moments produced by the actual Joint rotations are determined by multi
plying the moments produced by the unit displacements (see Fig. 18.4d
............ ,
694
Chapter 18
+ 92 (0.05E/) + 9 3 (0)
MBC
9 1 (0. lEI)
48.3 kipft
(18.23)
MCD
(18.24)
MEC] = 2El [2
[MCB
40 1
1] 1 [258.6]
2 El 448.3
[48.3]
57.8
(18.25)
(18.26)
(18.27)
150
(g)
Section 18.4
695
157.8
(h)
in Local Coordinates
.. ..... ';;.
.. .....
,.,;;.
696
Chapter 18
\,
y'
(a)
(b)
and
2EI
..
.
L (20 i + OJ - 3tf;)
(18.28)
2EI
MJ = -(20.
+ O - 3'/')
L
J.''I' .
(18.29)
I: :. j -I:::..,
tf; = - - " - -
(18.30)
2EI(
. 3
M=
I
L 20+0+
I
)
L
(18.32)
1:::...I
(18.33)
(18.34)
........... -
Section 18.4
(3
~
L2 L'l )
6
vJ = -2El
L -().
L + -().
L J + -A
L2
I
j.
697
. (18.35)
[Z}
Vi
2EI
L
Vi
3
L
L
3
3'
L
L
6
6
L2
L2
3
3
L
3
L
3
L
6
3
L - L2
t
(18.36)
I
I.
y
OJ
Ili
(a)
The 4 X 4 member stiffness matrix can also be derived using the basic
approach of introducing unit displacements at each one of the degrees of
freedom. The external forces, at the DOF, required to satisfy eqUilibrium
in each defomled configuration are the elements of the member stiffness
matrix in the column corresponding to that DOE Refer to Figure 18.6 for
the following derivations.
(}j=
It:l~;:::;;;::l:~
k' - 4El
11 L
k' _ 2El
21 -
k' _ 6EI
31 L2
k' __ 6EI
41 . L2
(18.37)
k' _ 6EI
32- L2
(c)
= 1 rad)
k' _ 2E1
12 L
k' _ 4El
22- L
k' _ 6El
32 L2.
(18.38)
r
I
\
698
Chapter 18
I
I
[2 1] [1]
[1]
2EI
~
= 6EI
L 1 2L 1
U 1
(18.39)
The end moments and con-esponding shears (calculated from statics) are
depicted in Figure 18.6d; again we have
12EI
' _ 6EI
k23 -
' _ 6EI
k13-
(18.40)
12EI
k' _ _ 12EI
k' __ 6EI
L2
L3
34 -
(18.41)
k'
2El
L
3
L
3
L
3
L
3
L
3
L
3
L
6
L2
3
L
3
L
6
L2
(18.42)
6
L2
Section 18.4
699
T[:;]
[()iC]
()Jc
(18.43)
Al
Aj
where T is the transformation matrix and the subscript c has been added
to distinguish between rotations measured with respect to the chord and
rotations with respect to the local axis x' .
The elements of the transformation matrix T can be obtained with the
aid of Figure 18.7. From there we have
0ic = ()i ()jc
= OJ -
if!
if!
(18.44)
(18.45)
'"
.,K"~......."..-:r~{lj
{ljr
A - A
t/J=
(18.30)
f).
Ojc
= ()j
Aj
+----!.
(18.46)
Aj
(18.47)
(18.48)
700
Chapter 18
From Section 17.7 we know that if two sets of coordinates are. geo
metrically related; then if the stiffness matrix is known in one set of coor
dinates, it can be transformed to the other by the following operation:
k' = TTkT
(18.49)
where k is the 2 x 2 rotational stiffness matrix. CEq. 18.6) and k' is the
4 X 4 member stiffness matrix in local coordinates. Substituting the T
matrix in Equation 18.48 and the rotational stiffness matrix of Equation
18.6 for k, we get
k'
o
1
L
1
1
1
L
I
The multiplication of the matrices shown above yields the same beam
element stiffness matrix as derived previously arid presented as Equation
18.42; the verification is left as an exercise for the reader.
EXAMPLE18.3
Analyze the plane frame shown in Figure I8.8a. The frame is made up
of two colUmnS of moment of inertia I, rigidly connected to a horizontal
beam whose moment of inertia is 31. The structure .Sllpports a concen
trated load of 80 kips acting horizontally to the right at the niidheight of
column AB. Neglect the deformations due to axial forces.
.
Solution
Because axial deformations are neglected, joints Band C do not move
vertically but have the same horizontal displacement. In Figure 18.8b we
use arrows to show the positive sense of the three independent joint dis
placement components. We now apply the five-step solution procedure
utilized in the preceding examples.
2
3
"I
8'
80 kips
(b)
(a)
........
--
Section 18.4
The 4
701
+ Ml
-160.0
Mj
M2 = 0
= 0
+ F3
40.0
= 160.0
F3 = -40.0
F [-16~.0]
(18.50)
40.0
where forcis are in kips and moments are in kipft.
,
..
'
The
deformed confIgurations, corresponding to unit displacements at each
degree of freedom, are shown in Figure 18.8d. The moments at the end of
c ,
160 BT
80 kips
160
.....
Ml
I
\
\......Jf
.,
M2
~
160
-+F3
(c)
..
702
Chapter 18
+ Kll = 0
-0.2EI + K21 = 0
0.0234EI + K31 = 0
-O.4El - 0.25El + K22 = 0
0.0234EI + K32 = 0
-0.25EI - O.4El
or
KIl
0.65EI
or
K21
= 0.20El
or
K31
or
= -0.0234E1
K22 = 0.65EI
or
K32
-0.0234El
The elements of the third row of the structure stiffness matrix are eval
uated by introducing a unit horizontal displacement at the top of the frame
K;~ . . .
\ Us
O.4EI
0.2El
~
0.'25EI
B
K12
MEl
0.25E[
)
: [Ii
.......
0;2EI
_ ..K22
0.25El
0.0234E[
L13 = 1
hC
K 13
I
-----~.-~
I 0.0234EI
'
I 0.0234EI
/1/
I
0.0234EI
23
'Us
"'-'
'-'
0.0234El
B
0.0234EI
0.0029EI
0.0029E[
,I 0.0234EI
....."'.-
_ ....K
........
D
(d)
..."'.- -
Section 18.4
703
(DOF 3). The forces in the members are calculated as follows. From Fig
ure 18.8d we see that for this condition member BC remains undeformed,
thus having no moments or shears. The columns, members AB and DC,
are subjected to the deformation pattern given by
[~}[~l
where the subscripts i and) are used to designate the near and the far joints,
respectively. Notice that by defining the columns as going from A to B
and from D to C, both local y axes are in accordance with the previously
established sign convention, directed to the right, thus making the dis
placement A = 1 positive.
The moments and shears in each column are obtained by substituting
the displacements shown above into Equation 18.36, that is,
rJ~
Vi
VJ
2EI
L
3
L
3
L
3
L
3
L
3
L
3
L
3
L
3
L
6
6
'L2 L2
6
L2
6
L2
Substituting L = 16 ft gives
[M]
[-0.0234]
M
-0.0234
= EI
j
Vi
-0.0029
0.0029
+ K33 = 0
or
K33
= 0.0058EI
0.65
0.20
K = EI
0.20 _'
0.65
[
-0.0234 -0.0234
'-0.0234]
-0.0234
0.0058
.":;;' 4IIioa
704
Chapter 18
EI
0.65
-O.0234][(JI ]
0.20
0.65
-0.0234 (J2
0.20
[
-0.0234 -0.0234 0.0058
60 3
-160.0]
0.0
[
40.0
Solving yields
(Jl ]
1 [-57.0]
O
= - 298.6
[ 602
EI 7793.2
3
[0]
(Jo
1 -57.0
(JA]
=
[ 60A
60B
EI
0.0
7793.2
(JO]
(Je =
[ 60B
1
. EI
[-57,0]
298.6
0
0 .
60e
::]
[ 60D
60e
1[
EI
29~'6]
0
7793.2
I!
Section 18.5
705
196.9
36.9
163.6
80 kips
---jIooo
189.&
349.8
145.3
(f)
(e)
in Local Coordinates
While virtually all members in real structures are subjectto both axial and
flexural deformations, it is often possible to obtain accurate solutions
using analytical models in which only one deformation mode (flexural or
axial) is considered. For example, as we showed in Chapter 17, the
analysis of trusses can be carried out using a member stiffness matrix
that relates axial loads and deformations; bending effects, although pres
ent (since real joints do not behave as frictioriless pins, and the dead weight
of a member produces moment), are negligible. In other structures, such as
beams and frames treated in the previous sections of this chapter, often the
axial deformations have a negligible effect, and the analysis can be carried
out considering bending deformations orily. When it is necessary to include
both deformation components, in this section we derive a member stiff
ness matrix in local coordinates that will allow us to consider both axial
and bending effects simultaneously.
When bending and axial deformations are considered, each joint has
3 degrees of freedom; thus the order of the membcrstiffncss matrixis 6.
Figure 18.9 shows the positive direction of the degrees of freedom (joint
displacements) in local coordinates; notice that the sign convention for
end rotations and transverse displacements (degrees of freedom 1 through
4) is identical to that previously used in the derivation of the member stiff
ness matrix given by Equation 18.36. The displacements in the axial direc
tion (degrees of freedom 5 and 6) are positive in the direction of the posi
tive x/axis, which, as stated previously, runsfrom the near to the far joint.
706
Chapter 18
kin
k62
k63
k64 =
(18.51)
k46 =
(18.53)
Notice that the coefficients in Equations 18.51 and 18.53 satisfy sym
metry (Betti's law).
Organizing all the stiffness coefficients in a matrix, we obtain the 6 X 6
member stiffness matrix in local coordinates as
DOF:
4EI
2EI
6EI
L2
4EI
6El
L2
L
=
6EI
k'
2EI
6EI
6EI
L'I.
L2
6EI
6EI
6El
12El -12El
12El
12EI
- L2
---
L3
IF
L2
0
(18.54)
AE
AE
AE
AE
........ . .
Section 18.5
707
in
Analyze the frame Figure I8.lOa, considering both axial and flexural
deformations. The flexural. and axial stiffnesses EI and AE are the same .
for both members and equal24 X 106 kip'in2 and 0.72 X 106 kips, respec
tively. The structure supports a concentrated load of 40 kips that acts ver
tically down at the center of span Be.
EXAMPLE 18.4
40 kips
Solution
With axial elongations considered, the structure has 3 degrees of kine
matic indeterminacy, as shown in Figure I8.10b. The five-step solution
procedure follows:
+ 20.0
+ 250.0 = 0
Y2
M3
=0
=0
Xl
or
l-..l
or
50'-----1
(a)
f:::. . . . ..
=0
Y2 = -20.0
or
I-- 30'
Reversing the sign of these restraining forces to construct the force vec
tor F gives
[ 0]
(b)
(18.55)
20.0
3000.0
40 kips
(c)
...,... .a- _
708
Chapter 18
The stiff
ness matrices in local coordinates for members AB and BC are identical
because their properties are the same. Substituting the numerical values
for EI, AE, and the length L, which is 600 in., into Equation 18.54 gives
k' = 102
1600
800
4
-4
0
0
800
1600
4
-4
0
0
-4
-4
-0.0133
0.0133
0
0
4
4
0.0133
-0.0133
0
0
0
0
0
0
12
-12
0
0
0
(18.56)
0
-12
12
(d)
A
B
dA
dB
OA
DB
0
0
0
0.8
0
0.6
(18,57)
0
0
0
0
1
0
(18.58)
(}e
dB
de
8B'
De
Mj
Vi
\'i
Fj
Fj
-320.0
-320.0
-1.064
1.064
-720.0
720.0
(18.59)
Section 18.5
709
M/
Mj
VI
\j
F,
0
0
0
0
1200.0
-1200.0
(18.60) .
+ (720
X 0.8)
=0
+ 320.0 = 0
(1.067 X 0.6)
if.31
=0
or Kll = 1632.85
or
K21
= -575.36
or
K31
= -320.0
In Figure 18.10e we show the deformed configuration for a unit displacement at degree of freedom 2. Proceeding as before, we find the mem
ber deformations. For member AB,
0
0.4
0
0
(18.61)
0.6
0
8A
-0.8
88
OB
llA
llB "'"
K22
, ..:-',K33
K I2 __
240
p t
)400
~ 0.8 1.33
960
(e)
tlB
0
0
08
(}c
tlc
=
0
88
8e
0
0
(18.62)
.:0. .....
71 0
Chapter 18
Mi
Mj
Vi
"i
Fi
Fj
-240.0
-240.0
-0.8
0.8
960.0
-960.0
(18.63)
400.0
400.0
1.333
-1.333
0
0
(18.64)
Mi
Mj
Vi
"i
Fi
FJ
Given the internal member forces, the external forces required for
equilibrium at the degrees of freedom are readily found; referring to the
free-body diagram of joint B in Figure I8.10e, we calculate the follow
ing stiffness coefficients:
K12
K22
=a
or
K 12
:::::
=0
400 = 0
or
K22
= 769.81
or
K32
+ 240
-575.36
160.0
AB
0
SA
0
SB
and for member BC,
eB
(Jc
AB
Ac
6B
6c
........ .-
..
(18.66)
Section 18.5
Cuunlillnh:~
711
8000
160,000
400
-400
Vi
Vj
Fj
Fj
(18.67)
F;
160,000
80,000
400
==
-400
Fj
Mi
Mj
Vi
Vj
(18.68)
K23 + 400
K33
+ 400
0.8 = 0
and
Ku = -320
and
K 23 = 160
160,000 = 0
and
K33
X 0.6 - 400
160,000
DOF3
320,000
K 23
.1632.85
K = -575.36
[
-320.0
(18.69)
-320.0][L\1]
160.0 L\,
320,000.0
03
[ 0]
20.0
3000.0
K13 __ ~33t )
'
400
160,000\ ) 400
160000
'
(f)
,~
=1
(18.70)
L\l] =
L\2
83
[0.014 ]
0.0345
0.00937
(18.71)
[continues on next page]
712
Chapter 18
+ (0.0345
X -0.8)
+ (0.00937
X (0.0)
-0.0192
=0
e8
= 0.00937
.il. A
a8
5A
8B
X -0.8)
= -0.0192
eB = 0.00937
fie = 0
aB = 0.0345
ac =
5B
= 0.014
5e
=0
M"AB
MBA
V"AB
=
V"BA
F"AB
F"BA
736.98
1486.71
3.706
-3.706
23.04
-23.04
.........
--
........
(18.72)
Section 18.6
MilBe
Men
V"Be
V"eB
F"Be
F"eB
713
126.1
1513.29
763.54
3.79
-3.79
16.80
16.80
(18.73)
/l12l9
"""'"<J
63.6
moments
(kip.ft)
61.4
The results given by Equations 18.72 and 18.73 are plotted in Figure
I8.lOg. Note that the units of moment in the figure are kipfeet.
+16.8
-23.04
+ tension
axial forces
(kips)
281.3
~
:;':;.
::ik~' .'...
123.9
(g)
'. :.'.
.
313.6
+16.8
V~Ion
a~:~:rm
~.
-23.04 .
(kips)
,,' :'>/J
(II)
H .......................U.H
n . . . . . . . . . . . . . . u
u . . . . ..
.
. in Global Coordinates
...........
714
Chapter 18
(a)
TTk'T
where k' is the member stiffness matrix in local coordinates (Eq. 18.54),
k is the member stiffness matrix in global coordinates, and T is the trans
fomlation matrix. The T matrix is formed from the geometric relationships
that exist between the local and the global coordinates. In matrix form
(b)
= TA
(18.75)
where 5 and A are the vectors of local and global joint displacements,
respectively.
Refer to Figure 18.I1a and b for the member ij expressed in the local
and global coordinate systems, respectively. Note that the components of
translation are different at each end, but the rotation is identical. The
relationship between the local displacement vector 5 and the global dis
placement vector A is established as follows. Figure I8.11e and d shows
the displacement components in the local coordinate system produced by
global displacements Aix and Aiy, at joint i, respectively. From the figure,
(c)
(d)
(18.74)
(COS)(AiY)
(18.76)
(18.77)
(18.78)
(18.79)
........
Section 18.6
(j)
(e)
715
Together with two identity equations for joint rotations (Oi = 0i and OJ
OJ), the relationship between Band 4. is
0
0
OJ
Ai
Aj
0;
8j
1
0
$
0
C
0 0 0
c -s 0
0 0 0
0
0
0
0
0
c
0
-$
0
0
0
$
0
1
0
0
0
0
Ah
Aiy
OJ
(18.80)
Aj ..
A jy
OJ
where 8 = sin cp, c = cos cp, and the 6 X 6 matrix is the transformation.
matrix T.
From Equation 18.74, the member stiffness matrix in global coordi
nates is
TTk'T
o0
s 0
o0 c 0
1 000
000 s
000 c
o1 0 0
c
-8
0
0
0
0
0
0
0
c
-8
4EI
L
2EI
L
6EI
2EI
L
4EI
L
6EI
6EI
6EI
L2
6EI
L2
6EI
L2
-12EI
L3
-----0-
12El
L3
AE
L
AE
L
AE
L
AE
L
12EI
L3
12EI
6EI
L2
6EI
a:-... . - _
...
. -.....
.-
0 0 1 0
0 0 o0
s c o0
0 0 o s
C -$ o 0
0 0 o C
0 0
0 1
0 0
c 0
0 0
-8 0
716
Chapter 18
NC2+p~
k = EI
L
(Nc 2 + PS2)
se(N - P)
-Qs
Ne 2 +Ps?
sc(-N+P)
Ns 2 + Pel
Qs
Qc
4
Symmetric about main diagonal
-sc(-N+P)
Qs
-(Ns2 + Pe2)
Qc
-Qe
.2
-Qs
se(-N+ P)
2
2
-Qe
Ns + Pe
4
.-.,
........."'.."'...............t
~n
(18.81)
u ~ H.;o u _
k = [kJ;, kJ;,F]
k}\'N
kF
(18.82)
EXAMPLE
18.5
Using the direct stiffness method, assemble the structure stiffness matrix
for the frame in Figure 18.10a.
Solution
Figure 18.12a illustrates the structure and identifies the degrees of freedom.
Note that the degrees of freedom are numbered in the order x, y, z and are
shown in the sense of the positive direction of the global axes; this order is
necessary to take advantage of the special form of Equation 18.82.
Since the frame considered has three joints, the total number of inde
pendent joint displacement components, before any supports are intro
duced, is 9. Figure 18.12b shows the stiffness matrices for the two mem
bers (expressed in the format of Eq. 18.82), properly located within the
Section 18.7
s s
717
kN I
kNF
r;:; rn
C
2
zr: x
~-l
2
kNF
3
S
global
axes
kFN
structure stiffness
matrix
kFN2
kl
k}+k~
(1 R.R~)
0.72
N =
gL2 = ~
=
6002
= 24.0 == 0.03 m-
33 33 X 10- 6 in- 2
Q ==
~L = ~
= 001 in- 1
600
.
E1
24.0 X 106
600
L ;:::
40,000 kipin
'.-:0.'- _
432.85-575.36
k} = . - 575.36
768.48
[
-320
....,.240
-320
-240
160,000
l'
(18.84)
718
Chapter 18
I'
.
and
k~
[1200
0
0
0
1.33
400
400
160,000
0]
(18.85)
1632.85
K = -575.36
[
-320
-575.36
-320]
769.81
160
320,000
160
Ans.
(18.86)
Summary
"", .
K.1 = F
719
Problems
16'
10 kips
P = 20 kips
8'
--*,-'-15'
6'
P1 B.1
r-
IO'
(El)beams = 2(El)columns
D
1---20'
6'~
P18.3
.1.
1.
35'
D
P18.2
,I.
,I.
E1 = constant
P1B.S
........
720
Chapter 18.
!..-- 10 / -------1
P18.7
....... u .... u ...................................... u . . . . . . . . . . . . . . . n.u ...... u ............... u ............ n ............. u ...... u ... u ....... n ....................... UUH ....... U
......... .-
u . . . . . . . ". . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
i'
I