STRUCTURAL-Analysis
STRUCTURAL-Analysis
COM
BASIC CONCEPTS
AND
CONVENTIONAL METHODS
OF
STUCTURAL ANALYSIS
(LECTURE NOTES)
ACKNOWLEDGEMENTS
I realize the profound truth that He who created all things inert as well as live is
GOD. But many things in this world are created through knowledge of some living
beings and these living beings are groomed by their teachers and through their own effort
of self study and practice.
Those who are gifted by God are exceptions and those who are gifted by their
teachers are lucky. I am thankful to God as He has been grateful to me much more than I
deserve and to my all teachers from school level to University heights as I am product of
their efforts and guidance and hence this work.
It is not out of place to mention some names that did excellent job of teaching and
showed path to learning, teaching and research. Mr. Sainani of Baba Thakurdas Higher
Secondary School, Lucknow, India did excellent job of teaching of mathematics at school
level. Mr. Tewarson of Lucknow Christian College, Lucknow, India was an excellent
teacher of mathematics at college level. Professor Ramamurty taught theory of Structures
excellently at I.I.T. Kharagpur, India. Visiting Professor Paul Andersen from U.S.A.
demonstrated the techniques of teaching through his very well prepared lectures and
course material on Structural Mechanics, Soil Mechanics and his consultancy practice.
Visiting Professor Gerald Picket taught Advanced Theory of Elasticity, Plates and shells.
It was a rare opportunity to be their student during Master of Engineering Course of
Calcutta University at Bengal Engineering College, Howrah, India.
The present work is the result of inspirations of my teachers and also any future
work that I may accomplish.
Last but not the least, I express my gratitude to Professor Tarun Kant, Present
Dean (Planning) I.I.T. Powai, Mumbai India and Ex. Head of Civil Engineering
Department who extended the facility of putting this work on web site of I.I.T. Powai,
Mumbai, India. Also I thank present Head of Civil engineering Department, I.I.T. Powai,
India. Professor G.Venkatachalam who whole heatedly extended the facility of typing
this manuscript in the departmental library by Mrs. Jyoti Bhatia and preparation of
figures in the drawing office by Mr. A.J. Jadhav and Mr. A.A. Hurzuk.
The inspiration and knowledge gained from my teachers and literature have
motivated me to condense the conventional methods of structural analysis in this work so
that a reader can get the quick insight into the essence of the subject of Structural
Mechanics.
In the end I wish to acknowledge specifically the efforts of Mrs. Jyoti Bhatia, who
faithfully, sincerely and conscientiously typed the manuscript and perfected it as far as
possible by reviewing and removing the errors.
Finally I thank Sri Anil Kumar Sahu who scanned all the figures and integrated
the same with the text and arranged the entire course material page wise.
REFERENCES
1. Andersen P. Statically Indeterminate Structures. The Ronald Press Company, New
York, U.S.A, 1953.
3. Junnarkar S.B. Mechanics of Structures, Valumes I & II. Chartor Book Stall, Anand,
India.
4. Mohan Kalani. Analysis of continuous beams and frames with bars of variable cross-
section. I. Indian Concrete Journal, March 1971.
5. Mohan Kalani. Analysis of continuous beams and frames with bars of variable cross-
section :2. Indian Concrete Journal, November 1971.
6. Norris C.H., Wilbur J.B & Utku S. Elementary Structural Analysis. McGraw – Hill
Book Company, Singapore.
8. Timoshenko S.P & Young D.H. Theory of Structures. McGraw – Hill Kogakusha
Ltd., Tokyo, Japan.
9. Vazirani V.N. & Ratwani M.M. Analysis of Structures, Khanna Publishers, Delhi,
India.
10. West H.H. Analysis of Structures. John Wiley & Sons, New York, USA.
CONTENTS
1. INTRODUCTION 1
2. CLASSIFICATION OF SKELEAL
OR FRAMED STRUCTURES 1
1 INTRODUCTION
The skeletal structure is one whose members can be represented by lines possessing
certain rigidity properties. These one dimensional members are also called bar members
because their cross sectional dimensions are small in comparison to their lengths. The
skeletal structures may be determinate or indeterminate.
1) Direct force structures such as pin jointed plane frames and ball jointed space
frames which are loaded and supported at the nodes. Only one internal force or
stress resultant that is axial force may arise. Loads can be applied directly on the
members also but they are replaced by equivalent nodal loads. In the loaded
members additional internal forces such as bending moments, axial forces and
shears are produced.
The plane truss is formed by taking basic triangle comprising of three members and three
pin joints and then adding two members and a pin node as shown in Figure 2.1. Sign
-2-
Convention for internal axial force is also shown. In Fig.2.2, a plane triangulated truss
with joint and member loading is shown. The replacement of member loading by joint
loading is shown in Fig.2.3. Internal forces developed in members are also shown.
The space truss is formed by taking basic prism comprising of six members and four ball
joints and then adding three members and a node as shown in Fig.2.4.
2) Plane frames in which all the members and applied forces lie in same plane as
shown in Fig.2.5. The joints between members are generally rigid. The stress
resultants are axial force, bending moment and corresponding shear force as shown
in Fig.2.6.
3) Plane frames in which all the members lay in the same plane and all the applied
loads act normal to the plane of frame as shown in Fig.2.7. The internal stress
resultants at a point of the structure are bending moment, corresponding shear force
and torsion moment as shown in Fig.2.8.
External forces including moments acting on a structure produce at any section along a
structural member certain internal forces including moments which are called stress
resultants because they are due to internal stresses developed in the material of member.
The maximum number of stress resultants that can occur at any section is six, the three
Orthogonal moments and three orthogonal forces. These may also be described as the
axial force F1 acting along x – axis of member, two bending moments F5 and F6 acting
about the principal y and z axes respectively of the cross section of the member, two
corresponding shear forces F3 and F2 acting along the principal z and y axes respectively
and lastly the torsion moment F4 acting about x – axis of member. The stress resultants at
any point of centroidal axis of member are shown in Fig. 3.1 and can be represented as
follows.
-3-
⎧F1 ⎫ ⎧Fx ⎫
⎪F ⎪ ⎪F ⎪
⎪ 2⎪ ⎪ y ⎪
⎪F ⎪ ⎪F ⎪
{F} = ⎪⎨ 3 ⎪⎬ OR ⎪⎨ z ⎪⎬
⎪F4 ⎪ ⎪M x ⎪
⎪F5 ⎪ ⎪M y ⎪
⎪ ⎪ ⎪ ⎪
⎪⎩F6 ⎪⎭ ⎪⎩M z ⎪⎭
Numbering system is convenient for matrix notation and use of electronic computer.
Each of these actions consists essentially of a pair of opposed actions which causes
deformation of an elemental length of a member. The pair of torsion moments cause twist
of the element, pair of bending moments cause bending of the element in corresponding
plane, the pair of axial loads cause axial deformation in longitudinal direction and the
pair of shearing forces cause shearing strains in the corresponding planes. The pairs of
biactions are shown in Fig.3.2.
In many frames some of six internal actions contribute greatly to the elastic strain energy
and hence to the distortion of elements while others contribute negligible amount. The
material is assumed linearly elastic obeying Hooke’s law. In direct force structures axial
force is primary force, shears and bending moments are secondary. Axial force structures
do not have torsional resistance. The rigid jointed plane grid under normal loading has
bending moments and torsion moments as primary actions and axial forces and shears are
treated secondary.
In case of plane frame subjected to in plane loading only bending moment is primary
action, axial force and shear force are secondary. In curved members bending moment,
torsion and thrust (axial force) are primary while shear is secondary. In these particular
cases many a times secondary effects are not considered as it is unnecessary to
complicate the analysis by adopting general method.
For the analysis of structures various loads to be considered are: dead load, live load,
snow load, rain load, wind load, impact load, vibration load, water current, centrifugal
force, longitudinal forces, lateral forces, buoyancy force, earth or soil pressure,
hydrostatic pressure, earthquake forces, thermal forces, erection forces, straining forces
etc. How to consider these loads is described in loading standards of various structures.
These loads are idealized for the purpose of analysis as follows.
-4-
Concentrated loads: They are applied over a small area and are idealized as point loads.
Line loads: They are distributed along narrow strip of structure such as the wall load or
the self weight of member. Neglecting width, load is considered as line load acting along
axis of member.
Surface loads: They are distributed over an area. Loads may be static or dynamic,
stationary or moving. Mathematically we have point loads and concentrated moments.
We have distributed forces and moments, we have straining and temperature variation
forces.
-5-
Equations of equilibrium
∑Fx = 0 ∑Mx = 0
∑Fy = 0 ∑My = 0
∑Fz = 0 ∑Mz = 0
For space frames number of equations of equilibrium is 6. Forces along three orthogonal
axes should vanish and moments about three orthogonal axes should vanish.
There are three equations of equilibrium. Forces in x and y directions should vanish and
moment about z axis should vanish.
-6-
⎧0 ⎫
⎪F ⎪
⎪ y ⎪
⎪⎪Fz ⎪⎪
Release for Axial Force (AF) Fx: ⎨ ⎬
⎪M x ⎪
⎪M y ⎪
⎪ ⎪
⎪⎩M z ⎪⎭
⎧Fx ⎫
⎪0 ⎪
⎪ ⎪
⎪⎪Fz ⎪⎪
Release for Shear Force (SF) Fy: ⎨ ⎬
⎪M x ⎪
⎪M y ⎪
⎪ ⎪
⎪⎩M z ⎪⎭
⎧Fx ⎫
⎪F ⎪
⎪ y ⎪
⎪⎪0 ⎪⎪
Release for Shear Force (SF) Fz: ⎨ ⎬
⎪M x ⎪
⎪M y ⎪
⎪ ⎪
⎪⎩M z ⎪⎭
⎧Fx ⎫
⎪F ⎪
⎪ y ⎪
⎪⎪F ⎪⎪
Release for Torsion Moment (TM) Mx: ⎨ z ⎬
⎪0 ⎪
⎪M y ⎪
⎪ ⎪
⎪⎩M z ⎪⎭
-7-
⎧Fx ⎫
⎪F ⎪
⎪ y ⎪
⎪⎪F ⎪⎪
Release for Bending Moment (BM) My: ⎨ z ⎬
⎪M x ⎪
⎪0 ⎪
⎪ ⎪
⎪⎩M z ⎪⎭
⎧Fx ⎫
⎪F ⎪
⎪ y ⎪
⎪⎪F ⎪⎪
Release for Bending Moment (BM) Mz: ⎨ z ⎬
⎪M x ⎪
⎪M y ⎪
⎪ ⎪
⎪⎩0 ⎪⎭
⎧Fx ⎫ ⎧0⎫
⎪F ⎪ ⎪0⎪
⎪ y⎪ ⎪ ⎪
⎪⎪F ⎪⎪ ⎪⎪0⎪⎪
Universal joint (Ball and socket joint) F = ⎨ z ⎬ , Cut F= ⎨ ⎬
⎪0 ⎪ ⎪0⎪
⎪0 ⎪ ⎪0⎪
⎪ ⎪ ⎪ ⎪
⎪⎩0 ⎪⎭ ⎩⎪0⎭⎪
A release does not necessarily occur at a point, but may be continuous along whole length
of member as in chain for BM. On the other hand a constraint is defined as that which
prevents any relative degree of freedom between two adjacent nodes connected by a
member or when a relative displacement of the nodes does not produce a stress resultant
in the member.
∝ s = P (M – N + 1) – r = PR – r
∝ k = P (N – 1) + r – c
-8-
∝ s + ∝ k = PM –c
R = (M – N + 1)
∝ s = M - (NDOF) N + P
NDOF = Degrees of freedom at node which is 2 for plane truss and 3 for space truss.
-9-
It may be noted that structure may be mechanism even if ∝ s > 0 if the releases are
present in such a way so as to cause collapse as mechanism. The situation of mechanism
is unacceptable.
Statical Indeterminacy
It is difference of the unknown forces (internal forces plus external reactions) and the
equations of equilibrium.
Kinematic Indeterminacy
It is the number of possible relative displacements of the nodes in the directions of stress
resultants.
Using formula:
∝ s = P (M – N + 1) – r
∝ k = P (N – 1) + r – c
P = 6, M = 7, N = 6
∝ s = P (M – N + 1) – r = 6 (7 – 6 + 1) – 3 = 9
∝ k = P (N – 1) + r – c = 6 (6 – 1) + 3 – 12 = 21
-10-
∝ s + ∝ k = PM – c = 6 x 7 – 12 = 30
Static indeterminacy can also be determined by introducing releases in the system and
rendering it a stable determinate system. The number of biactions corresponding to
releases will represent static indeterminacy. Consider a portal frame fixed at support
points as shown in Fig.6.4. The entire structure is shown in Fig.6.5 and completely stiff
structure in Fig.6.6.
∝ s = P (M – N + 1) – r
∝ k = P (N – 1) + r – c
P = 3, M = 4, N = 4, c = 3, r = 0
∝ s = 3 (4 – 4 + 1) – 0 = 3
∝ k = 3 (4 – 1) + 0 – 3 = 6
∝s + ∝k = 3 + 6 = 9
The structure can be made determinate by introducing in many ways three releases and
thus destroying its capacity to transmit internal forces X1, X2, X3 at the locations of
releases.
In figure 6.7. a cut is introduced just above clamp D that is clamp is removed. It becomes
tree or cantilever structure with clamp at A. At this cut member was transmitting three
forces X1, X2 and X3 (Two forces and one moment). Therefore ∝ s = 3. This is external
static indeterminacy.
In figure 6.8. a cut is introduced at point R on member BC. We have two trees or
cantilevers with clamps at A and D. We have three internal unknown forces X1, X2, and
X3. Thus ∝ s = 3.
In figure 6.9. three hinges are introduced. We have determinate and stable system and
there are three unknown moments X1, X2 and X3. Thus ∝ s = 3.
In figure 6.10. one roller cum hinge and one hinge is introduced. We have one unknown
force X1 and two unknown moments X2 and X3 at these releases. Thus ∝ s = 3.
The static and kinematic indeterminacies of a few structures are computed in Table 1.
-11-
Example Figure P M N R c r ∝s ∝k
No: No:
1 6.11 3 4 3 2 6 3 3 3
2 6.12 3 2 2 1 5 2 1 0
3 6.13 3 2 2 1 3 1 2 1
4 6.14 3 12 9 4 6 2 10 20
5 6.15 3 7 6 2 3 3 3 15
6 6.16 3 12 6 7 25 19 2 9
7 6.17 3 13 6 8 28 20 4 7
8 6.18 3 6 2 5 0 0 15 3
14 10 5 24 0 15 3
9 6.19 6 9 7 3 12 0 18 24
10 6.20 3 4 3 2 6 6 0 5
These are the two basic methods by which an indeterminate skeletal structure is analyzed.
In these methods flexibility and stiffness properties of members are employed. These
methods have been developed in conventional and matrix forms. Here conventional
methods are discussed.
Flexibility Method
-12-
In flexibility method since unknowns are forces at the releases the method is also called
force method. Since computation of displacement is also required at releases for
imposing conditions of compatibility the method is also called compatibility method. In
computation of displacements use is made of flexibility properties, hence, the method is
also called flexibility method.
Stiffness Method
Since nodal displacements are unknowns, the method is also called displacement method.
Since equilibrium conditions are applied at the joints the method is also called
equilibrium method. Since stiffness properties of members are used the method is also
called stiffness method.
(1) Obtain the reactions at the supports of structure applying appropriate equations of
equilibrium.
(2) Separate the members at the joints as free bodies and apply equations of equilibrium
to each member to obtain member end forces.
(3) Cut the member at a section where internal forces are required. Apply equations of
equations to any of the two segments to compute unknown forces at this section.
Example 8.1
Compute reactions for the beam AB loaded as shown in figure 8.1. Also find internal
forces at mid span section C.
-13-
Detach the beam from supports and show unknown reactions as shown in Fig.8.2
The reaction RB which is perpendicular to rolling surface is replaced with its horizontal
and vertical components RBX and BBY.
3 4
RBX = RB Sin θ = RB, RBY = RB Cos θ = RB
5 5
At A reaction in vertical direction is zero and other components are RAX and MAZ.
Resultant of triangular load W is shown acting at 8m from A and 4m from B that is
through CG of triangular loading. The free body diagram with known forces is shown in
Fig.8.3.
1
W= x 50 x 12 = 300 kN
2
∑ Fx = 0, ∑ Fy = 0 and ∑ Mz = 0
Alternatively, ∑ Fx = 0, ∑M
A
z = 0, ∑M
B
z =0
∑M
B
z = 0 gives : MAZ = 4 x 300 = 1200 kNm
2400 + 1200
∑M
A
z = 0 gives : 12 RBY = 8 x 300 + MAZ, RBY =
12
= 300 (check)
5 3
RB = x 300 = 375, RBX = x 375 = 225.
4 5
Now the beam is cut at mid span and left segment is considered as a free body.
The free body diagram of segment AC with unknown forces is shown in Fig.8.4.
1
Total triangular load = x 6 x 25 = 75 kN. It acts at 4m from A and 2m from C.
2
-14-
∑ Fy = 0 gives, RCY = 75 kN
∑M
C
z = 0 gives, MCZ = 1200 – 75 x 2 = 1050 kNm
Example 8.2
Determine the reactions for the three hinged arched frame ABC loaded as shown in Fig.
8.5. Show free body diagrams for members AB and BC and segments BD and DC.
We have three equations of equilibrium and four unknown reactions. The structure is
determinate despite four unknown reactions as the moment at hinge B is zero. The free
body diagrams of members AB and BC are shown in Fig.8.6 and Fig.8.7.
∑M
A
z = 0, 3 RBX + 4 RBY = 40 x 2.5 = 100 ……… (3)
∑M
C
z = 0, 4 RBX - 3 RBY = - 25 + 10 x 1.5 + 5 x 2 = 0 ……… (6)
-15-
1 144
From (7), RBX = [400 – 16 x 16] = = 12
12 12
9 x 16
From (8), RBX = = 12 (check)
12
The free body diagrams of members AB and BC with known forces are shown in Figures
8.8 and 8.9.
Member BDC is shown horizontally and the forces are resolved along the axis of member
(suffix H) and normal to it (suffix V) as shown in figure 8.10.
3 4
At B : RBH = 12 cos θ + 16 sin θ = 12 x + 16 x = 20
5 5
4 3
RBV = 12 sin θ - 16 cos θ = 12 x - 16 x = 0
5 5
4 3
At D : RDH = 10 sin θ - 5 cos θ = 10 x -5x =5
5 5
4 3
RDV = - 5 sin θ - 10 cos θ = - 5 x - 10 x = - 10
5 5
3 4
At C : RCH = +7 cos θ + 26 sin θ = + 7 x + 26 x = + 25
5 5
4 3
RCV = - 7 sin θ + 26 cos θ = - 7 x + 26 x = 10
5 5
It can easily be verified that equations of equilibrium are satisfied in this configuration.
By cutting the member just to left of D the free body diagrams of segments are shown in
Fig. 8.11.
-16-
The trusses are classified as determinate and indeterminate. They are also classified as
simple, compound and complex trusses. We have plane and space trusses. The joints of
the trusses are idealized for the purpose of analysis. In case of plane trusses the joints are
assumed to be hinged or pin connected. In case of space trusses ball and socket joint is
assumed which is called universal joint. If members are connected to a hinge in a
plane or universal joint in space, the system is equivalent to m members rigidly
connected at the node with hinges or socketed balls in (m-1) number of members at the
nodes as shown in figure 9.1. In other words it can be said that the members are allowed
to rotate freely at the nodes. The degree of freedom at node is 2 for plane truss (linear
displacements in x and y directions) and 3 for space truss (linear displacements in x,y and
z directions). The plane truss requires supports equivalent of three reactions and
determinate space truss requires supports equivalent of six reactions in such a manner
that supporting system is stable and should not turn into a mechanism. For this it is
essential that reactions should not be concurrent and parallel so that system will not rotate
and move. As regards loads they are assumed to act on the joints or points of concurrency
of members. If load is acting on member it is replaced with equivalent loads applied to
joints to which it is connected. Here the member discharges two functions that is function
of direct force member in truss and flexural member to transmit its load to joints. For this
member the two effects are combined to obtain final internal stress resultants in this
member.
The truss is said to be just rigid or determinate if removal of any one member destroys its
rigidity and turns it into a mechanism. It is said to be over rigid or indeterminate if
removal of member does not destroy its rigidity.
Relation between number of members and joints for just rigid truss.
Number of equivalent links or members or reactive forces to constrain the truss in space
is 6 corresponding to equations of equilibrium in space ( ∑ Fx = 0, ∑ Fy = 0, ∑ FZ = 0,
∑ Mx = 0, ∑ My = 0, ∑ MZ = 0). For ball and socket (universal) joint the minimum
number of links or force components for support or constraint of joint in space is 3
corresponding to equations of equilibrium of concurrent system of forces in space ( ∑ Fx
= 0, ∑ Fy = 0, ∑ FZ = 0). Each member is equivalent to one link or force.
Total number of links or members or forces which support j number of joints in space
truss is (m + 6). Thus total number of unknown member forces and reactions is (m + 6).
The equations of equilibrium corresponding to j number of joints is 3j. Therefore for
determinate space truss system: (m + 6) = 3j.
-17-
m = (3j – 6)
Minimum just rigid or stable space truss as shown in Fig.9.2. is a tetrahedron for which
m = 6 and j = 4. For this relation between members and joints is satisfied.
m = 3 x 4 – 6 = 6 (ok)
By adding one node and three members the truss is expanded which can be supported on
support system equivalent of six links or forces neither parallel nor concurrent. We get
determinate and stable system. As can be seen joints 5 and 6 are added to starting stable
and just rigid tetrahedron truss. Three links at each of two joints 3 and 6 corresponding to
ball and socket joint are provided.
Plane truss
The stable and just rigid or determinate smallest plane truss as shown in Fig.9.3.
comprises of a triangle with three nodes and three members. Two members and a pin
joint are added to expand the truss. Total number of non-parallel and non-concurrent
links or reactive forces required to support j number of joints is 3. Total number of
unknowns is number of member forces and reactions at the supports. Number of
available equations is 2j. Therefore for determinate plane truss system:
(m + 3) = 2j
m = (2j – 3)
Hinge support is equivalent of two reactions or links and roller support is equivalent of
one reaction or link.
Exceptions
The stable and just rigid or determinate truss is shown in figure 9.6, m = 9, j = 6, m = 2 j
– 3 = 2 x 6 – 3 = 9. The relation between members and joints will also be satisfied if
arched part is made horizontal as shown in Fig.9.7. The system has partial constraint at C
as there is nothing to balance vertical force at pin C. The two members must deflect to
support vertical load at C. In fact the rule for forming determinate simple truss is violated
as joint 1 is formed by members 1 and 2 by putting them along same line because these
are the only two members at that joint.
-18-
Compound truss
Compound plane truss is formed by joining together two simple plane trusses by three
nonparallel and nonconcurrent members or one hinge and the member. Compound truss
shown in figure 9.8 is formed by combining two simple trusses ABC and CDE by hinge
at C and member BE. It is shown supported at A and B. For purpose of analysis after
determining reactions at supports the two trusses are separated and unknown forces X1,
X2 and X3 are determined by applying equations of equilibrium to any one part. There-
after each part is analyzed as simple truss. This is shown in Fig.9.9.
Compound truss shown in figure 9.10 is formed by combining the two simple trusses by
three nonparallel and nonconcurrent members. The truss is supported by two links
corresponding to hinge support at A and one link corresponding to roller at B. By cutting
these three members the two parts are separated and the unknown forces X1, X2 and X3 in
these members are determined by equations of equilibrium and each part is analyzed as
simple truss. This is shown in Fig.9.11.
In case of compound space truss six members will be required to connect two simple
space trusses in stable manner so that connecting system does not turn into a mechanism.
Alternatively one common universal ball and socket joint and three members will be
required. The method of analysis will be same as in plane truss case.
Complex truss
A complex truss is one which satisfies the relation between number of members and
number of joints but can not be configured by rules of forming simple truss by starting
with triangle or tetrahedron and then adding two members or three members and a node
respectively for plane and space truss. A complex truss is shown in figure 9.12.
M = 9, j = 6, m = 2j – 3 = 2 x 6 – 3 = 9
There are two methods of analysis for determining axial forces in members of truss under
point loads acting at joints. The forces in members are tensile or compressive. The first
step in each method is to compute reactions. Now we have system of members connected
at nodes and subjected to external nodal forces. The member forces can be determined
by following methods.
-19-
The method of joints is used when forces in all the members are required. A particular
joint is cut out and its free body diagram is prepared by showing unknown member
forces. Now by applying equations of equilibrium the forces in the members meeting at
this joint are computed. Proceeding from this joint to next joint and thus applying
equations of equilibrium to all joints the forces in all members are computed. In case of
space truss the number of unknown member forces at a joint should not be more than
three. For plane case number of unknowns should not be more than two.
Method of sections
This method is used when internal forces in some members are required. A section is
passed to cut the truss in two parts exposing unknown forces in required members. The
unknowns are then determined using equations of equilibrium. In plane truss not more
than 3 unknowns should be exposed and in case of space truss not more than six
unknowns should be exposed.
Example 9.1
Determine forces in all the members of plane symmetric truss loaded symmetrically as
shown in figure 9.13 for all members by method of joints and in members 2,4 and 5 by
method of sections.
∑ Fx = 0 gives, R3 = 0
∑M
A
Z = 0 gives, 30 R2 = 1000 x 10 + 1000 x 20 = 30,000, R2 = 1000 kN
-20-
Method of joints
Joint A
Free body is shown in figure 9.14. Force in member 1 is assumed tensile and in member
3 compressive. Actions on pin at A are shown.
∑F
A
X = 0 : F1 – F3 cos 450 = 0,
∑F
A
Y = 0 : - F3 sin 450 + 1000 = 0, F3 = 1000 2 = 1414 kN,
1
F1 = 1000 2 x = 1000 kN
2
Since positive results are obtained the direction and nature of forces F1 and F3 assumed
are correct. At joint C there will be three unknowns, hence, we proceed to joint B where
there are only two unknowns.
Joint B
1
∑F
B
X = 0 gives, F3 cos 450 – F4 = 0, F4 = 1000 2 x
2
= 1000 kN
1
∑F
B
y = 0, gives : F6 + F3 cos 450 = 0, F6 = - 1000 2 x
2
= - 1000 kN
The negative sign indicates that direction of force assumed is wrong and it would be
opposite. It is desirable to reverse the direction of F6 here it self and then proceed to joint
C, else the value will have to be substituted in subsequent calculation with negative sign
and there are more chances of making mistakes in calculations. The corrected free body
diagram of joint B is shown in figure 9.16.
Joint C
The free body diagram for joint C is now prepared and is shown in figure 9.17.
∑F
C
Y = 0 gives, F5 cos 450 = 0, F5 = 0
∑F
C
X = 0 gives, F2 = 1000 kN
-21-
The results are shown in figure 9.18. The arrows shown at the ends of members are forces
actually acting on pin joints. The reactive forces from joints onto members will decide
whether it is tension or compression in the members. The sign convention was explained
in theory.
Method of sections
Now a section is passed cutting through members 2, 4 and 5 and left segment is
considered as a free body as shown in Fig.9.19. The unknown member forces are
assumed tensile. However, if it is possible to predict correct nature, the correct direction
should be assumed so as to obtain positive result. A critical observation of free body
indicates that F5 = 0 as its vertical component can not be balanced as remaining resultant
nodal forces in vertical direction vanish. Now equilibrium in horizontal direction
indicates that F4 = - F2. The segment is subjected to clockwise moment of 10,000 kNm,
hence, F2 and F4 should form counter clockwise couple to balance this moment. This also
indicates force F4 should have opposite direction but same magnitude. Since arm is 10 m,
F2 x 10 = 10,000, hence, F2 = 1000 kN. and F4 = - 1000 kN. By method of sections we
proceed as follows:
∑M
D
Z = 0 gives : F2 x 10 + 1000 x 10 – 1000 x 20 = 0, F2 = 1000 kN
∑M
C
Z = 0 gives : - F4 x 10 – 1000 x 10 = 0, F4 = - 1000 kN
1
∑ Fy = 0 gives : - F5 x - 1000 +1000 = 0, F5 = 0
2
1
∑ FX = 0 gives : - F5 x + F2 + F4 = 0, F5 = 0
2
θx, θy, θz = angle that axis of member AB makes with x, y and z axis respectively.
-22-
LAB = (x B − x A )2 + (y B − y A )2 + (z B − z A )2
Tension coefficient t for a member is defined as tensile force T in the member divided by
its length L.
T T
t= , tAB = AB = tension coefficient for member AB.
L L AB
(x B − x A )
TAB cos θx = TAB = tAB (xB – xA)
L AB
( yB − y A )
TAB con θy = TAB = tAB (yB – yA)
L AB
(z B − z A )
TAB cos θz = TAB = tAB (zB – zA)
L AB
TAB, TAC, TAD = Unknown tensile forces acting on members AB, AC and AD at joint A.
-23-
tAB (xB – xA) + tAC (xC – xA) + tAD (xD – xA) + QAX + PAX = 0
tAB (yB – yA) + tAC (yC – yA) + tAD (yD – yA) + QAY + PAY = 0
tAB (zB – zA) + tAC (zC – zA) + tAD (zD – zA) + QAZ + PAZ = 0
These equations can be written in compact form by identifying any member with far and
near ends.
(x B − x A )
Component of pull TAB in x-direction = TAB cos θx = TAB = tAB (xB – xA)
L AB
(y B − y A )
Component of pull TAB in y-direction = TAB cos θy = TAB = tAB (yB – yA)
L AB
LAB = (x B − x A ) 2 + ( y B − y A ) 2
Compact form of equations of equilibrium at joint A is:
-24-
Example 9.2
For the shear leg system shown in figure 9.23 determine the axial forces in legs and tie
for vertical load of 100 kN at the apex (head). Length of each leg is 5 m and spread of
legs is 4 m. The distance from foot of guy rope to center of spread is 7 m. Length of guy
rope is 10 m.
OC = 7 m, AB = 4 m, AC = BC = 2 m, OH = 10 m, AH = BH = 5 m.
CH = 52 − 22 = 21 = 4.5826 m
Cos θ =
(OH 2
+ OC 2 − CH 2
=
) ( )
10 2 + 7 2 − 21
(2OC x OH ) (2 x 7 x 10)
Cos θ = 0.9143
OD = 10 Cos θ = 9.143 m
HD = 10 Sin θ = 4.051 m
CD = 9.143 – 7 = 2.143 m
Node x y z
O 0 0 0
H 0 9.143 4.051
A -2 7 0
B 2 7 0
-25-
tHA (xA– xH) + tHB (xB – xH) + tHO (xO – xH) = 0 _____ (1)
tHA (yA – yH) + tHB (yB – yH) + tHO (yO – yH) = 0 _____ (2)
tHA (zA – zH) + tHB (zB – zH) + tHO (zO – zH) + PHY = 0 ____ (3)
From eqn (2): -2 x 2.143 tHA = 9.143 tHO, ∴ tHA = -2.1332 tHO
From Eqn (3): - 4.051 (-2.1332 – 2.1332 – 1) tHO = 100, tHO = 7.5573
10.1 Cables
Cable is a very efficient structural form as it is almost perfectly flexible. Cable has no
flexural and shear strength. It has also no resistance to thrust, hence, it carries loads by
simple tension only. Cable adjusts its shape to equilibrium link polygon of loads to which
it is subjected. A cable has a shape of catenary under its own weight. If a large point load
W compared to its own weight is applied to the cable its shape changes to two straight
segments. If W is small compared to its own weight the change in shape is insignificant
as shown in figure 10.1. From equilibrium point of view a small segment of horizontal
length dx shown in Fig.10.2 should satisfy two equations of equilibrium ∑Fx = 0 and ∑Fy
= 0. The cable maintains its equilibrium by changing its tension and slope that is shape.
One unknown cable tension T can not satisfy two equilibrium equations, hence, one
additional unknown of slope θ is required. The cables are used in suspension and cable
stayed bridges, cable car systems, radio towers and guys in derricks and chimneys. By
assuming the shape of cable as parabolic, analysis is greatly simplified.
-26-
n
∑W = ∑ W = W1 + W2 + ---- + Wi + ----- + Wn
i =1
ΣWi a i ΣWi b i
a= , b=
ΣW ΣW
∑ M = b ∑W
B
∑M
RA = Vertical reaction at A = B
- H tan ∝
L
∑M
RB = (∑W + H tan ∝ - B
)
L
Consider a point X on cable at horizontal coordinate x from A and vertical dip y from
chord.
Since cable is assumed to be perfectly flexible the bending moment at any point of cable
is zero. Considering moment equilibrium of segment of cable on left of X the relation
between H, x and y is obtained which defines general cable theorem.
-27-
H (XX1) = ∑ M - RA x
X
∑M
H (x tan ∝ - y) = ∑M - (
X
B
L
- H tan ∝) x
⎡x ⎤
Hy = ⎢ ∑ M - ∑ M ⎥
⎣L B X ⎦
b
∑M
VA = ∑ W = B
L L
∑M
MX = VAx - ∑M =X
B
L
x- ∑M
X
Thus: Hy = MX
The general cable theorem therefore states that at any point on the cable subjected to
vertical loads, Hy the product of horizontal component of tension in cable and the vertical
dip of that point from cable chord is equal to the bending moment MX at the same
horizontal coordinate in a simply supported beam of same span as cable and subjected to
same vertical loading as the cable.
RA
TA = H 2 + R 2A , θA = tan-1
H
RB
TB = H 2 + R 2B , θB = tan-1
H
ωL x2
Hy = MX = x- ω
2 2
-28-
L
At mid span x = and y = h the dip of cable.
2
ωL2
H=
8h
Hy = MX
ωL2 ωL x2
y= x-ω
8h 2 2
4h
y= x (L – x)
L2
This is the equation of cable curve with respect to cable chord. The cable thus takes the
shape of parabola under the action of udl. The same equation is valid when chord is
horizontal as shown in Fig.10.5.
⎛ ⎛ dy ⎞
2 ⎞
L
S = ∫ ds = ∫ ⎜ L
1+ ⎜ ⎟ ⎟ dx
O⎜ ⎟
O
⎝ ⎝ dx ⎠ ⎠
dy 4h
= 2 (L – 2x)
dx L
1
⎡ 16h 2 2⎤
S = ∫ ⎢1 + 2 (L - 2x ) ⎥ dx
L 2
O
⎣ L ⎦
⎡ 8 h2 32 h4 256 h 6 ⎤
S = L ⎢1 + 2
− + − ......⎥
⎣ 3 L 5 L4 7 L6 ⎦
-29-
h 1
For flat parabolic curves ≤ , only two terms are retained.
L 10
⎡ 8 h2 ⎤
S = L ⎢1 + ⎥
⎣ 3 L2 ⎦
10.5 Example
RA =
(300 x 30 + 1 x 40 x 20) = 245 N
40
RB = (300 + 40 x 1) – 245 = 95 N
H 2 + 245 2 = 1000
H = 970 N
Considering segment left of P, the clockwise moment at P is computed and set to zero
since cable is flexible.
2400
h= = 2.474 m.
970
10.6. Arches
An arch is a curved beam circular or parabolic in form supported at its ends and is
subjected to inplane loading. The internal forces developed in the arch are axial force,
shear force and bending moment. Depending upon number of hinges the arches are
-30-
classified as (1) three hinged arch (2) two hinged arch (3) single hinged arch and (4) fixed
arch as shown in figures 10.7 to 10.10. A three hinged arch is statically determinate. The
remaining three are statically indeterminate to first, second and third degree respectively.
Here only determinate three hinged arch will be considered.
Arch under vertical point loads shown in figure 10.11 is a three hinged arch subjected to
vertical loads W1 to Wn. The reactions developed at the supports are shown. It may be
noted that moment at the hinge at C in the arch is zero, hence, horizontal component of
reaction can be computed from this condition.
bΣW aΣW
RA = , RB = , ∑W = W1 + ….. + Wn
L L
L
∑ M = RA
C 2
- Hh
⎛ bΣ W ⎞
⎜ − ∑ M ⎟⎠
H= ⎝
L C
∑ M = Counter clockwise moment about C of all applied vertical loads acting left of C
C
At any point P on the arch as shown in Fig.10.12, the internal forces Fx,Fy and Mz can
easily be computed as explained previously. From Fx and Fy shear force and thrust in the
arch can be computed.
V = Shear at P
C = Thrust at P
M = Bending moment at P
M = Mz
V = Fy cos θ - Fx Sin θ
C = - Fy Sin θ - Fx cos θ
-31-
4h
y= x (L – x)
L2
ωL
RA = RB =
2
L ωL L
RA - Hh - =0
2 2 4
ωL ⎛ ωL2 ⎞ ωx 2
Mx = x - ⎜⎜ ⎟⎟ y -
2 ⎝ 8h ⎠ 2
ωL ωx 2 ωL2 ⎡ 4h ⎤
= x- - ⎢ x (L - x)⎥
2 2 8h ⎣ L2
⎦
ωL ωx 2 ωLx ωx 2
= x- - + =0
2 2 2 2
10.8. Example
A three hinged parabolic arch of 20 m span and 4 m central rise as shown in Fig.10.14
carries a point load of 40 kN at 4 m horizontally from left support. Compute BM, SF and
AF at load point. Also determine maximum positive and negative bending moments in
the arch and plot the bending moment diagram.
-32-
4h 4x4
y= 2
x (L – x) = x (20 – x)
L 400
x
y= (20 – x)
25
4 16
RB = x 40 = 8 kN, RA = x 40 = 32 kN
20 20
MC = 0, 4H = 32 x 10 – 40 x 6 = 80, H = 20 kN
0 ≤x≤ 4m
x 4 2
Mx = 32 x – 20 (20 – x) = 16 x + x
25 5
4 x 16
x = 4, Mx = 16 x 4 + = 76.8 kNm
5
4 m ≤ x ≤ 20 m
x 4 2
Mx = 32 x – 20 (20 – x) – 40 (x – 4) = 160 – 24 x + x
25 5
4
x = 15, Mx = 160 – 24 x 15 + x 225 = - 20 kNm
5
x = 20, Mx = 0 (ok)
An influence line is a graph or curve showing the variation of any function such as
reaction, bending moment, shearing farce, axial force, torsion moment, stress or stress
resultant and displacement at a given section or point of structure, as a unit load acting
parallel to a given direction crosses the structure. The influence line gives the value of the
function at only one point or section of the structure and at no other point. A separate
influence line is to be drawn for the function at any other point.
-33-
There are two methods of construction of influence lines for determinate and
indeterminate structures.
In the direct method, first response function and its sign convention are identified.
Conventional free body and equilibrium are used to obtain the value of response function
for a number of positions of unit load placed along the axis of members of structure. The
response function values are plotted as influence line curve. The response function can
also be expressed as function of coordinate x measured from a reference point for various
segments of structure and then plotted as IL.
For simply supported beam AB of span L shown in Fig.11.1, the IL diagrams for
reactions and bending moment and shear force at section X are plotted as the vertical unit
load rolls from A to B along the axis of beam.
RB =
x
, RA =
(L - x )
L L
The above equations are for straight line hence, IL will also be a straight line
x = 0, RA = 1, RB = 0
x = L, RA = 0, RB = 1
If a horizontal unit force moves along axis of member, the horizontal reaction H at the
hinge will be unity. Consequently the IL diagram will be a rectangle with ordinate unity.
x = 0 to L, H = 1
-34-
The directions identified for RA and RB are vertical upwards and direction identified for
H is horizontal to left.
The directions of internal forces Vx and Mx at section X are identified as shown in figure
11.3. Unit load is placed at coordinate x.
0 ≤x≤ a
x x bx
RB = , Vx = , Mx =
L L L
a ≤x≤ L
RA =
(L - x ) , V =
(x - L ) , M =
a( L - x )
x x
L L L
x = 0, Vx = 0, Mx = 0
a ab
x = a, Vx = , Mx = (load just to left of X)
L L
-b ab
x = a, Vx = , Mx = (load just to right of X)
L L
x = L, Vx = 0, Mx = 0
SF is positive when load is just to leave of section X and it is negative when it is just to
the right of section. The BM is positive for all positions of load.
A four panel truss of span L and height h is shown in figure 11.4. Length of each panel is
a. It is required to plot influence lines for forces in members 1,2 and 3 as a unit load
moves along bottom chord from A to B.
-35-
I L for F1
MD
F1 = (compressive)
h
MD = BM at joint D.
Since height of truss is constant the IL for F1 is obtained by drawing the IL for moment at
D and dividing its ordinates by h. The IL will be triangle with ordinate at D equal to
(2a )(2a ) = a .
hL h
IL for F2
MC
F2 = (Tensile)
h
a (3a ) 3a
IL for moment at C is a triangle with ordinate =
4ah 4h
IL for F3
The vertical equilibrium of the parts of the truss on either side of the section xx requires
that the vertical component of force F3 should balance whatever forces may be imposed
on these parts that is ∑V=0.
h
F3 sin θ + RB = 0, sin θ =
a + h2
2
a2 + h2
F3 = - RB cosec θ, cosec θ =
h
-36-
The negative sign indicates that the actual force in member 3 in compressive so long as
load is to the left of joint J. The IL in this region may therefore be drawn by drawing IL
for RB and multiplying the ordinates by cosec θ.
F3 sin θ - RA = 0, F3 = RA cosec θ
The positive sign indicates that the force in member is tensile so long as the load is right
of D. The IL in this region is drawn by plotting the IL for RA and then multiplying the
ordinates by cosec θ.
The variation is linear. In fact the IL for diagonal member is proportional to the IL for the
shear in panel.
Unit load at J
1 1
RB = , F3 = - cosec θ
4 4
Unit load at D
1 1
RA = , F3 = cosec θ
2 2
12.1. Deformations
When a structure is subjected to the action of applied loads each member undergoes
deformation due to which the axis of structure is deflected from its original position. The
deflections also occur due to temperature variations and misfit of members. The
infinitesimal element of length dx of a straight member (ds of curved member) undergoes
axial, bending, shearing and torsional deformations as shown in figure 12.1 to 12.7. It is
assumed that the material of member obeys Hooke’s law. Small displacements are
considered so that structure maintains geometry.
Axial deformation
F1x dx
∈x dx =
EA
E = Modulus of elasticity
EA = Axial rigidity
∈x = Strain in x-direction
∈x dx = α ΔT dx
Bending deformations
Bending deformations which occur about y & z axes comprise of relative rotations of the
sides of the infinitesimal element through an angle dθy and dθz respectively.
F5x dx
dθy = = kydx
EI y
-38-
F5x
ky = = Elastic curvature of axis of member in xz-plane.
EI y
F6x dx
dθz = = kzdx
EI z
F6x
kz = = Elastic curvature of axis of member in xy-plane.
EI z
If the element as shown in Fig.12.4. is subjected to linear temperature change from ΔTt at
top to ΔTb at bottom, the angle change dθz due to this effect will be
d = depth of member
Shearing deformations
The deformations dδy and dδz due to shearing forces consist of displacement dδ of one
side of element with respect to other with respect to y and z directions.
-39-
F2x dx
dδy = μy
GA
F3x dx
dδz = μz
GA
G = Shear modulus
GA = Shear rigidity
μy, μz = nondimensional factors depending solely upon the shape and size of cross-
section which accounts for the nonuniform distribution of shearing stresses. For
A
rectangular section μ = 1.2 and for circular section μ = 1.11. For I or H sections can be
μ
taken as web area or in other words μ can be taken as ratio of area of cross-section to web
⎛ A ⎞
area ⎜⎜ μ = ⎟.
⎝ A w ⎟⎠
Torsional deformation
It is given by angle of twist dθx in radians, which represents the difference in the angles
of rotation of its faces about longitudinal axis of the element.
F4x dx
dθx =
GI x
-40-
L⎛ F
2
μ y F2x2 μ z F3x2 F2 F2 F2 ⎞
U r = ∫ ⎜ 1x + + + 4x + 5x + 6x ⎟ dx
o ⎜ 2EA 2GA 2GA 2GI x 2EI y 2EI z ⎟⎠
⎝
U= ∑U r
r
= Elastic strain or potential energy of all members of structure.
The energy due to shear is neglected being very small compared to that due to other
actions.
In plane frames subjected to inplane loading primary action is bending moment only
hence energy due to axial force is neglected.
Considering only the relevant primary actions the elastic strain energy for various
structures is given by following expressions.
F1x2
∑ Ur = ∑ ∫
L
U= dx
r r
o 2EA
F6x2 dx
∑U ∑ ∫
L
r
U= =
r r
o 2EI z
-41-
⎛ F4x2 F6x2 ⎞
∑U
L
∑ ∫ ⎜⎜ + ⎟⎟ dx
r
U= =
r r
o
⎝ 2GI x 2EI z ⎠
⎛ F1x2 F4x2 F2 F2 ⎞
∑ Ur = ⎜ + 5x + 6x ⎟ dx
L
U= ∑ ∫o ⎜ 2EA
+
2GI x 2EI y 2EI z ⎟⎠
r r ⎝
There are various methods developed for computation of displacements depending upon
structural system and nature of loading but basic methods are based on energy principles
such as Castigliano’s theorem and virtual work.
The partial derivative of elastic strain energy U of the structure with respect to any
external load P is equal to displacement δ in the structure corresponding to that force. The
terms force and displacement are used in the generalized sense that is word force may
mean force or moment and word displacement may mean linear or angular displacement.
Strain energy U is function of P.
∂U
δ=
∂P
If the deflection is required at a point where there is no load, a load P is placed there and
in the expression for partial derivative of elastic energy P is set equal to zero. If the
deflection is required in the direction of a particular defined load the load is replaced with
P and finally P is set equal to prescribed value.
The deflection calculations are some what simplified if the partial derivatives are worked
out before integration.
-42-
The derivatives represent the rate of change of forces Fex with respect to P.
∂F1x ∂F ∂F ∂F
f1x = , f 4 x = 4x , f 5x = 5x , f 6x = 6x
∂P ∂P ∂P ∂P
These derivatives f are equal to values Fex as caused by a unit load (P=1) and are
represented by fex.
⎡ F1x f1x
⎢ ∫L ∫L F4x f 4x ∫L F5x f 5x ∫L F6x f 6x ⎤⎥
δ = ∑⎢ dx + dx + dx + dx ⎥
r
⎢ EA GI x EI y EI z ⎥
⎣ ⎦
The equations of dummy load method are derived from principle of virtual work, hence,
it is also called virtual work method. It is also called Maxwell-Mohr method.
The opposing pair of dummy unit loads is deployed to obtain relative displacement or
rotation of two points on the structure.
According to principle of virtual work if the second system is given a small displacement
the total work of the forces will be zero.
At a point represented by local coordinate x in a member the internal stress resultants will
be fex (e = 1,….,6) due to dummy unit action. The virtual displacements of the second
system are taken as the actual displacements of the first system. Then in accordance with
the principle of virtual work.
⎡ F f F f F f F f ⎤
1 x δ = ∑ ⎢ ∫ 1x 1x dx + ∫ 4x 4x dx + ∫ 5x 5x dx + ∫ 6x 6x dx ⎥
⎢ L EA
r ⎣ L
GI x L
EI y L
EI z ⎥⎦
-43-
It will be observed that the deflection calculations using Castigliano’s theorem are the
same as in dummy load method.
Example
Determine the horizontal and vertical deflections and the angle of rotation at free end A
of cantilever bracket shown in figure 12.8 neglecting axial and shear deformations. The
members AB and BC have flexural rigidity EI and axial rigidity EA. Determine
additional deflection at A if axial deformations are considered.
The structural system comprises of two members. The members are numbered 1 and 2.
The local axes xyz of members are shown. The common axes system for whole structure
is XYZ. The bending moment diagrams due to given loading, unit horizontal and vertical
and unit couple at A are shown in Fig.12.9. The displacements are computed in system
coordinates XYZ.
Vertical deflection at A
F6x f 6x
δY = -∑∫ dx
r L
EI
a Px x Pa 3
r = 1: ∫
o EI
dx =
3EI
h Pa a Pa 2 h
r = 2: ∫
o EI
dx =
EI
⎡ Pa 2 h Pa 3 ⎤
δY = - ⎢ +
⎣ EI 3EI ⎥⎦
Horizontal deflection at A
F6x f 6x
δX = + ∑ dx
r EI
-44-
a Px x 0
r = 1: ∫o EI
dx = 0
h Pa (h - x ) Pa h 2 Pah 2 Pah 2
r = 2: ∫o EI
dx =
EI
−
2EI
=
2EI
Pah 2
δX =
2EI
Rotation at A
F6x f 6x
θZ =− ∑ dx
r EI
a Px x 1 Pa 2
r = 1: ∫o EI
dx =
2EI
h Pa x 1 Pah
r = 2: ∫o EI
dx =
EI
⎡ Pa 2 Pah ⎤
θZ =− ⎢ + ⎥
⎣ 2EI EI ⎦
F1x f1x
δY = -∑ dx
r EA
h P x1 Ph
r = 2: ∫o EA
dx =
EA
ph
δY = -
EA
-45-
Example
RA = RB = 3 kN
3 4
Sin θ = , cos θ =
5 5
3x5
F1 sin θ = 3, F1 = = 5 kN (C)
3
3x5
F2 sin θ = 3, F2 = = 5 kN (C)
3
4
F3 = F2 cos θ = 5 x = 4 kN (T)
5
1 5
F1 sin θ = , F1 = (T)
2 6
5
F2 sin θ = 1, F2 = (T)
6
5 4 2
F3 = F1 cos θ = x = (C)
6 5 3
-46-
3
8 RB = 3 x 1, RB = ↑
8
3
8 RA = 3 x 1, RA = ↓
8
Joint A
3 5
F1 sin θ = , F1 = (T)
8 8
5 1 1
F3 = 1 - cos θ = 1 - = (T)
8 2 2
Joint B
3 5
F2 sin θ = , F2 = (C)
8 8
Joint C
5
2x cos θ = 1 (check)
8
F1 = F2 = 0
F3 = 1 (T)
-47-
1 5 5 25 25
1 (CA) 5000 150 -5 + + 0 - 0
- 6
6 8 6 8
f1XC , f1YC, f1XB = Axial force in member due to unit loads in X and Y directions at point
C and B.
3
F1f1L 2 25 1 25 1 4
δxc = ∑
r =1 AE
=2x
5
+
8
x
6
-
8
x
6
=
5
= 0.8 mm
3
F1f1L 8 2 25 1 25 1
δyc = ∑
r =1 AE
= - x
3 5
-
6
x
6
-
6
x
6
= 2.46 mm
3
F1f1L 2
δXB = ∑r =1 AE
=4x
5
= o.8 mm
Check
(ΔL)3 = ⎛⎜
FL ⎞ 4 x 8000
⎟ = 1.6 mm.
⎝ AE ⎠ 3 100 x 200
Example
-48-
a = 100 mm, h = 200 mm, φ = 100 mm, E = 200 kN/mm2, G = 70 kN/mm2, P = 200 kN
The two members are separated as shown in Fig.12.12. Their local axes x1, y1, z1 and x2,
y2, z2 are shown. The free bodies of the two members are prepared and BM and TM
diagrams are shown for P and P = 1.
⎡ F4x f 4x ⎤
⎢∫ ∫F f
⎥
2 5x 5x
δ = ∑⎢L dx + L
dx ⎥
r =1
⎢ GI x EI y ⎥
⎣ ⎦
Ix = Iy + Iz
πφ 4 π x 100 4 4 6 4
Iy = Iz = = = 4.909 x 10 6 mm , Ix = 9.818 x 10 mm
64 64
Member no I
a
Px 2 Pa 3 200 x 100 3
δA1 = ∫o EI y dx = =
3EI y 3 x 200 x 4.909 x 10 6
= 0.068 mm
Member 2
The member no.2 is subjected to primary actions of bending moment about y2-axis and
torsion moment about x2-axis.
h h
1 1 Ph 3 Pa 2 h
δA2 = ∫ Px dx + ∫ Pa dx = +
2 2
EI y o
GI x o
3EI y GI x
-49-
Check
13 NONPRISMATIC MEMBERS
A nonprismatic member has a variable section along its length. For any structural
member, the stress resultants do not remain constant throughout its length. Where there is
a significant variation of the stress resultants along the length of member, economy can
be achieved more efficiently by varying the cross-sectional area of the member, keeping
in view the extreme values of the stress resultants in the middle and end sections. Often
such variation can be adopted to add to the architectural appearance of the structure.
Generally, two types of members are used:
The axis of member with haunches or with variable depth is assumed to be the same as
for the uniform part usually in the central portion of the member. For a tapered member
the uniform part will correspond to the minimum section.
Analysis of such members shown in Fig.13.1 to 13.4 involves the determination of fixed
end reactions due to self weight and loads acting on the members and the flexibility and
stiffness properties that is force-displacement relationships. For haunched members the
fixed end moments, stiffness and flexibility properties are expressed in terms of integrals
or bar constants. The bar constants are available in handbooks.
-50-
Concentrated vertical load F and udl of intensity ω will be considered as shown in Fig.
13.5.
x (x - a )dx (x - a )dx
L L L L
xdx x 2 dx
∫o I z ∫a I z − ∫o I z ∫ Iz
(M z )A = F a
2
x 2 dx ⎛ xdx ⎞
L L L
dx
∫o I z ∫o I z − ⎜⎜⎝ ∫o I z ⎟⎟⎠
2
xdx x 3dx ⎛ x 2 dx ⎞
L L L
∫ ∫
ω o Iz o Iz
− ⎜∫ I ⎟
⎜
⎝o z ⎠
⎟
(M z )A =
2 L dx L x 2 dx ⎛ L xdx ⎞ 2
∫o I z ∫o I z − ⎜⎜⎝ ∫o I z ⎟⎟⎠
The fixed end moment at the right end B can be determined by performing the
integrations in the opposite direction. The vertical reactions at the two ends can be
determined from the equilibrium of the free body of the member. The case of bending in
xz-plane is treated similarly.
The angle of twist between any section on the uniform part of the haunched beam and the
fixed end is given by:
-51-
T ⎛ LH LU ⎞
θ= ⎜ + ⎟
G ⎜⎝ J H J U ⎟⎠
b 3 (d 2 − d1 ) d
JH = for > 1.2
⎛ d − 0.63b ⎞ b
3 log e ⎜⎜ 2 ⎟⎟
⎝ d1 − 0.63b ⎠
b 3d ⎛ b⎞ d
JU = ⎜1 - 0.63 ⎟ for > 1.2
3 ⎝ d⎠ b
The torsional stiffness factor Kt for a member with haunched and uniform parts is given
by:
G
Kt =
⎛ LU LH ⎞
⎜⎜ + ⎟⎟
⎝ JU JH ⎠
T (K t )A
TA =
(K t )A + (K t )B
T (K t )B
TB =
(K t )A + (K t )B
The fixed end reactions for the case of uniformly distributed torque can be worked out in
similar manner.
-52-
The fixed end reactions PA and PB for axial force P acting at any point on uniform part of
a haunched member as shown in Fig.13.10 are given by:
P (K F )A
PA =
(K F )A + (K F )B
P (K F )B
PA =
(K F )A + (K F )B
E
KF =
⎛ LU L ⎞
⎜⎜ + H ⎟⎟
⎝ AU AH ⎠
b(d 2 − d1 )
AH =
d
log e 2
d1
The equivalent uniform area Aeq of member of variable section is given in terms of area
Ax at coordinate x.
L
1
1 ∫A dx
= o x
A eq L
EA eq
Axial rigidity =
L
The fixed end reactions for the case of uniformly distributed axial force can be worked
out in a similar manner.
Like fixed end reactions the bar constants are also expressed in terms of integrals.
Evaluation of integrals from first principles is cumbersome hence, use is made of hand
books to save this labour.
13.2 Basic concepts and definitions of bar constants of members with variable section
-53-
The rotational stiffness factor K at one end of member, which is assumed hinged, is
defined as the moment required to produce a unit rotation at this end, while the other end
is assumed to be fixed as shown in Fig.13.11.
In hand books the stiffness factor K at an end of a haunched member is given as:
kEI min
K=
L
k = Stiffness coefficient
E = Modulus of elasticity
L = Length of member
The carry over factor C from the near end to the far end of a bar is defined as the ratio of
the moment induced at far end, which is assumed fixed, to the applied moment at the near
end which is assumed hinged. The product of the carry over factor and the rotational
stiffness factor at one end of member AB is equal to similar product at the other end.
The modified rotational stiffness factor K ′ at one end of the member is defined as the
moment required to produce unit rotation at this end which is assumed hinged when the
other end is also hinged as shown in Fig.13.12.
-54-
The lateral- stiffness factor K at one end of the member, which is assumed fixed against
rotation, is defined as the moment produced at this end when the other end also fixed
against rotation, is displaced laterally through unit distance with respect to the first end,
as shown in Fig.13.13.
K AB
K AB = (1 + CAB)
L
K BA
K BA = (1 + CBA)
L
′
The modified lateral stiffness factor K at one end of a member, which is assumed fixed
against rotation is defined as the moment produced at this end when the other end,
assumed hinged, is displaced laterally through unit distance with respect to first end as
shown in Fig.13.14.
K AB
K ′ AB = (1 - CAB CBA)
L
K BA
K ′ BA = (1 - CBA CAB)
L
Sign convention
Clockwise moments and rotations at the end of bar are treated positive. The lateral
displacement between two ends of bar resulting in positive end moments or counter
clockwise lateral displacement angle ψ is treated as positive. The shear causing clockwise
moment in the bar is treated as positive. The external loads acting vertically downwards
are treated positive as shown in Fig.14.1.
The force displacement relationship at the nodes of the elastic bar is defined by following
slope deflection equations:
Δ
MAB = MFAB + KAB θA + CBA KBA θB + KAB (1 + CAB)
L
-55-
Δ
MBA = MFBA + KBA θB + CAB KAB θA + KBA (1 + CAB)
L
Designating A as near end and B as far end of member, the two slope deflection
equations can be written as:
Δ
Mnear = MFnear + Knear θnear + Cfar Kfar θfar + Knear (1 + Cnear)
L
1
Cnear = Cfar =
2
4EI
Knear = Kfar =
L
2EI ⎡ 3Δ ⎤
Mnear = MFnear + ⎢ 2θ near + θ far +
L ⎣ L ⎥⎦
Example
Find the fixed end moments and reactions for the fixed beam AB loaded as shown in
figure 14.2 and due to downward sinking of support B of 10 mm. E = 15 kN/mm2.
1.2 240
aA = = 0.2, rA = = 0.6
6 400
1.8 160
aB = = 0.3, rB = = 0.4
6 400
For 4 kN load
-56-
For 3 kN load
250 x 400 3
Δ = 10 mm, L = 6000 mm, Imin = = 1333.333 x 106 mm4
12
Δ 10
MFAB = KAB (1 + CAB) = x 216 x 105 (1.622) = 58,392 kN mm
L 6000
Δ 10
MFBA = KBA (1 + CBA) = x 222 x 105 x 1.604 = 59348 kN mm
L 6000
Total fixed end moments due to loads and relative displacement of supports:
-57-
RA = 5.339 kN
RB = 1.661 kN
RA + RB = 7 kN (check)
Example
A 4 m span fixed tapered beam AB of section 250 x 560 mm at A and 250 x 280 mm at B
is subjected to udl of intensity 5 kN/m as shown in Fig.14.4. Using hand book of bar
constants determine fixed end moments and hence reactions.
280
aA = 1.0, aB = 0, rA = = 0.5, rB = 0
560
FEM = coeff x w x L2
1
RA = [9.728 + 5 x 4 x 2 – 4.232] = 11.374 kN
4
1
RB = [4.232 + 5 x 4 x 2 – 9.728] = 8.626 kN
4
-58-
The slope deflection method involves solution of many simultaneous equations. The
moment distribution method does not involve as many simultaneous equations. The
number of equations reduces considerately and in case of structures with no lateral
displacements of nodes of members there are no equations to be solved. Hence, the
computational work is much less in this method compared to other methods. This method
has an additional advantage as it is based on iterative technique consisting of series of
cycles of clamping and unclamping of the joints of structure, each cycle converging on
the precise final result. Therefore, the series of cycles can be terminated when the desired
degree of accuracy of results has been achieved.
The moment at each end of a member as can be seen from slope-deflection equations
comprises of four separate effects.
(1) Fixed end moment (FEM) at the end under consideration due to applied load on the
member.
In moment distribution method these effects are considered separately and results are
superimposed. Consider a joint A where four members meet rigidly as shown in figure
15.1. The far ends of members are clamped. An artificial clamp is introduced at A. The
structure is now kinematically determinate comprising of fixed ended members and the
fixed end moments can be computed for the external loads applied to members. The
degree of freedom of this structure is one that is rotation at A. Since far ends are clamps
they neither rotate nor displace laterally.
If the algebraic sum of fixed end moments at A does not vanish, the resultant moment
acting at the joint is termed unbalanced moment. The artificial clamp holds or balances
this unbalanced moment.
If now the artificial clamp is removed the joint will rotate under the action of unbalanced
moment as shown in Fig.15.2 and end moments will develop in all the members as shown
in Fig.15.3. The moments developed at A due to rotation restore the equilibrium of joint
and are called distributed moments. The moments developed at far ends are called carry-
over-moments. Same sign convention is followed as in slope deflection method.
-59-
When clamp at A is removed or relaxed as it is called the structure deforms further and
develops distributed member end moments at A and carry over moments at B1 to B4.
MAB1 = KAB1 θA
MAB2 = KAB2 θA
MAB3 = KAB3 θA
MAB4 = KAB4 θA
⎡ ⎤
θA ⎢∑ K AB ⎥ + M = 0
⎣A ⎦
-M
θA =
∑ K AB
A
- K AB1
MAB1 = M = - (DF)AB1 M
∑ K AB
A
Similarly expression for moments at end A for other members can also be written.
K AB1
(DF)AB1 =
∑ K AB
A
The distributed moment developed at the relaxed end of member under unbalanced
moment M at corresponding joint is equal to distribution factor (DF) of this member at
this joint times the unbalanced moment with the sign reversed.
The carry over moment at end B1 will be: MB1A = CAB1 MAB1
-60-
The carry-over moment to far end is equal to carry over factor times the corresponding
distributed moment and has the same sign.
If far end of a member is hinged then modified stiffness of the member is to be taken in
computation of distribution factors as shown in figure 15.4. B1 and B3 are clamps and B2
and B4 are hinges.
∑K
A
AB = K1 + K ′2 + K 3 + K ′4
K ′AB2
(DF)AB2 =
∑ K AB
A
The cantilever portions of the structure can be taken care of in two ways
(1) The effect of loads acting on cantilever can be transferred to nearest joint and
cantilever part removed as shown in figure 15.5. Nearest joint will be subjected to
externally applied vertical load and a moment. This moment is to be taken at this joint
as external moment. For the case under consideration the vertical load will be (P+Q)
and moment will be (PL+Qa) as shown in figure 15.5. Now the continuous beam
ABCD can be analyzed by moment distribution method. D is now a simple support.
Distribution factor for DC will be 1. There will be no carry over from C to D. For CD
modified stiffness is used for the moment distribution process.
(2) The continuous beam ABCDE is retained as it is and the cantilever part DE is not
knocked out. Cantilever moment is entered as FEM for DE. The DF at D will be 1 for
DC and zero for DE. Since member is now continuous over support D there will be
carry over of moments from C to D and rotational stiffness of CD will not be
modified in view of fact that D is not an end simple support as for case when
cantilever is knocked out. The moments DE and DC will be of same magnitude and
equal to cantilever moment but of opposite sign.
It may be noted that cantilever arm has no restraining effect on the rotation of joint to
which it is rigidly connected or in other words its stiffness is zero. Any unbalanced
moment is therefore carried or distributed entirely by other members meeting at this joint.
In the moment distribution process all the hinged joints should be released first for quick
convergence. As regards other joints one should start with the joint which has largest
unbalanced moment for rapid convergence. The final results are not affected by the order
of relaxation of joints.
-61-
In the structures where only rotation of joints is involved no simultaneous equations are
required to be solved. Such structures have no side sway. For structures undergoing side
sway the number of equations to be framed and solved will be equal to number of
unknown lateral displacements Δ or the side sway angles or chord angles ψ.
For problems involving side sway the solution by moment distribution is carried out in
two separate parts.
For this artificial constraints are introduced in structure so that there is no sway and
reactions are determined at these constraints.
Methodology
The methodology is explained on single bay two storey frame as shown in figure 15.6.
Degree of freedom with respect to side sway = 2.
This frame will have to be analyzed for (n+1) = (2+1) = 3 cases. The side sway is
prevented by introducing two constraints as shown in figure 15.7 and the frame is
analyzed by moment distribution for no sway case and reactions R1 and R2 determined at
the fictitious constraints. Now constraint at E is given unit displacement to right without
imparting any rotations to joints B, E, C and F as shown in Fig.15.8. The frame is
analyzed and reactions R1A and R2A are determined. Now constraint at F is given unit
-62-
R1 + R1A Δ1 + R1B Δ2 = 0
R2 + R2A Δ1 + R2B Δ2 = 0
Δ1 Δ Δ
Ψ1 = , Ψ2 = 2 , Ψ = 1
H1 H2 H
The two equations are solved for Δ1 and Δ2 and moments obtained in the two sway cases
are modified by substituting the values of Δ1 and Δ2 respectively. Now sway and non
sway moments are added to get final result.
It may be mentioned here that Kani’s method which is also based on iteration technique
does not involve solution of any simultaneous equations. Kani’s method will not be
discussed here.
Example 15.1.
22,200 33,350
(DF)BA = = 0.4, (DF)BC = = 0.6
55,550 55,550
-63-
The moment distribution for concentrated moment at joint B is carried out separately. It
is entered on top of the joint and then distributed. The moments at the ends of members
for the two cases are now added to obtain final moments.
The reactions can be computed from free body diagrams of supports shown in Fig.15.11.
JOINT A B B C
DF & COF 0.604 ← 0.4 0.6 → 0.294
-5
JOINT A B B C
DF & COF 0.604 ← 0.4 0.6 → 0.294
DC - 2.00 + 3.00
COC + 1.21 + 0.88
TOTAL + 1.21 + 2.00 + 3.00 + 0.88
-64-
Final moments
-5
JOINT A B B C
DF & COF 0.604 ← 0.4 0.6 → 0.294
Beam AB
RA =
40 x 4.8 30 x 1.8
+ +
(50.52 − 29.42) = 32 + 9 + 3.52 = 44.52 kN
6 6 6
Beam BC
24.42 + 0.2
RBC = + 5 x 2 = 6.12 + 10 = 16.12 kN
4
24.42 + 0.2
RC = 5 x 2 - = 10 – 6.12 = 3.88 kN
4
Example 15.2
The portal frame with unequal legs loaded as shown figure 15.12 is analyzed by moment
distribution method.
-65-
MFAB = MFBA = 0
4E (2I ) 3E (I )
KBC = = 2EI, KBA = = EI
4 3
∑ K = 3 EI
B
1 2
(DF)BA = , (DF)BC =
3 3
4EI 2EI
KCB = 2 EI, KCD = =
6 3
2EI 8
∑ K = 2 EI +
C 3
= EI
3
3 3 2 3 1
(DF)CB = 2 x = , (DF)CD = x =
8 4 3 8 4
Fixed-end-moments, distribution factors and carry over factors are entered as shown and
the distribution and carry over cycles of moment distribution method are carried out. The
final moments for no sway case are obtained.
-66-
JOINT A B B C C D
DF & COF 0← 1
3
2
3 →0.5 0.5 ← 3
4
1
4 →0.5
Horizontal shears are computed from free body diagrams of columns shown in Figure
15.14.
6.19
RAH = = 2.063 kN
3
30.46 - 29.05
RCH = 10 x 3 + = 30 + 0.23 = 30.23 kN
10
-67-
6EI 2
MFAB = MFBA = = EI
3x 3 3
6EI EI
MFCD = MFDC = =
6x 6 6
2 60
MFAB = MFBA = EI x = 40 kNm
3 EI
EI 60
MFCD = MFDC = x = 10 kNm
6 EI
6EIΔ
= 40, Δ = 60/EI
9
6EIΔ
= 10, Δ = 60/EI
36
60
These proportional FEMs are corresponding to Δ = .
EI
These FEMs are entered for process of moment distribution which is carried out and final
proportional moments obtained as shown in scheme of moment distribution.
-68-
JOINT A B B C C D
DF 1 2 3 1
3 3 4 4
COF 0.5← →0.5 0.5← →0.5
From the free body diagrams of columns shown in Figure 15.16 the horizontal shears are
obtained.
R AH =
13.82
= 4.61 kN, R DH =
(9.06 + 9.52) = 3.10 kN
3 6
-69-
MF = Modification factor
R + MF x R = 0
- R − 27.702
MF = = = 3.593
R − 7.71
The proportional moments for sway case are corrected by multiplying with modification
factor and added to moments obtained for nonsway case to obtain final moments for the
given problem.
MAB = 0
The continuous beams and plane frames are indeterminate structures. The following two
methods will be discussed here for analysis of such structures.
-70-
In these methods the effect of axial deformations is neglected. These are in fact
displacement methods. The unknown displacements involved are rotations and lateral
displacements at the joints. The joints of the structure are assumed to be rigid, hence, the
member end displacements for all the members meeting at the joint are the same. This is
known as the condition of compatibility. Besides compatibility, the moment equilibrium
of the joints and the lateral shear equilibrium of the storeys of frame should be satisfied.
In slope-deflection equations, the member end moments are expressed in terms of
member end rotations and relative lateral displacement and loads acting on member.
Therefore it is possible to express joint and storey equilibrium equations in terms of
unknown displacements of structure by using slope-deflection equations of all members.
This gives the system of linear algebraic equations, which is solved by conventional
methods in case of slope-deflection method and by an iterative physical relaxation
technique in case of moment distribution method. Final results are obtained in the form of
member end forces.
Slope-deflection method
Δ
Mnear = MFnear + Knear θnear + Cfar Kfar θfar + K near (1 + C near )
L
4EI
Knear = Kfar =
L
1
Cnear = Cfar =
2
The slope deflection equation for uniform member takes the following form.
2EI ⎡ 3Δ ⎤
Mnear = MFnear +
L ⎢⎣2θ near + θ far + L ⎥⎦
2EI
Mnear = MFnear + [2θnear + θfar + 3ψ]
L
-71-
The free body diagram of a rigid joint A of a plane structure is shown in figure 16.1. The
joint exerts positive clockwise moments at the ends of members meeting at this joint and
is in turn reacted by the same moments in the opposite counter clockwise directions
together with an external positive moment MextA acting on it.
An external counter clockwise moment directly acting on the joint is taken positive. The
joint moment equilibrium equation is written as follows:
∑ Mnear + Mext = 0
∑M
A
AB + MextA = 0
∑M
A
AB = The sum of the moments acting at end A of member AB for all the members
This pertains to the vertical columns of a particular storey in a multistory frame and is
derived from the shear equilibrium conditions of that storey. The columns are assumed
free of lateral loads directly acting on them. If there are lateral loads directly acting on the
columns, the same are replaced by statically equivalent horizontal loads and concentrated
moments acting on the nodes of frame and columns assumed free of external loads for
their function in the frame.
Inclined columns and columns of different heights can be handled by writing shear
equilibrium equation accordingly.
The free body diagram of columns of a particular storey is shown in figure 16.2. The
equation of equilibrium for a typical column AB shown in figure 16.3 is written as
follows:
-72-
M AB + M BA
HAB = -
h AB
Let Qr be the external storey shear for the rth storey of frame, which is taken positive
when acting to right. It is equal to sum of all the external horizontal loads acting on the
frame above this storey.
∑H
r
AB = Reactive rth storey shear which is sum of base shears of all the columns of this
storey.
Qr - ∑H
r
AB =0
(M AB + M BA )
∑r h AB
+ Qr = 0
First storey
The columns of first storey are rigidly connected to the nodes of frame at their upper
ends, their lower ends being either clamped or hinged to the foundation while the
columns of the other upper storeys are rigidly connected at both the ends. In order that
the rth storey shear equilibrium equation is also valid for the first storey, the hinged
column need to be replaced with equivalent column fixed at base. For structural similarity
of the original hinged column and the substitute fixed column the following conditions
should be satisfied as shown in figure 16.4.
MABO = MABS
θO = θS
ΔO = ΔS
-73-
3
KABS = KABO
4
3
hABS = hABO
2
KABO, KABS = Rotational stiffnesses of original hinged and substitute fixed columns.
Shown in figure 16.5 is a two span continuous beam of variable section. The continuous
beam is fixed at A and C and is roller supported at B. The degree of freedom is one that is
rotation θB at B which is unknown.
θA = θC = 0
The haunched member AB and the tapered member BC were discussed in previous
sections hence the required properties are reproduced below.
Member AB Member BC
MAB = MFAB + KAB θA + CBA KBA θB MBC = MFBC + KBC θB + CCB KCB θC
MBA = MFBA + KBA θB + CAB KAB θA MCB = MFCB + KCB θC + CBC KBC θB
MAB = - 44.46 + 13,408.8 θB MBC = - 9.37 + 33,350 θB
MBA = 39.45 + 22,200 θB MCB = 4.23 + 9804.9 θB
-74-
Moment equilibrium equation of joint B
The reactions are computed as follows using free body diagrams shown in figure 16.6.
Member AB
Member BC
RB2 =
5x4
+
(24.427 + 0.197 ) = 10 + 6.16 = 16.16 kN
2 4
RC =
5x4
+
(24.427 + 0.197) = 10 – 6.16 = 3.84 kN
2 4
Roller at B
-75-
A portal frame loaded as shown in figure 16.8 is analyzed for member end forces by
slope-deflection method.
Rotations are taken positive clockwise at the joints and lateral displacement Δ at the
beam level is taken positive to left.
Since D is clamp θD = 0
40 x 2 x 2 x 2
MFBC = - = - 20 kNm
4x4
40 x 4
MFCB = + = + 20 kNm
8
MFAB = MFBA = 0
10 x 6 x 6
MFCD = - = - 30 kNm
12
10 x 6 x 6
MFDC = + = + 30 kNm
12
2 EI ⎛ 3Δ ⎞
MAB = + ⎜ 2θ A + θ B + ⎟
3 ⎝ 3 ⎠
2 EI ⎛ 3Δ ⎞
MBA = + ⎜θ A + 2θ B + ⎟
3 ⎝ 3 ⎠
MBC = - 20 +
4 EI
(2θ B + θ C )
4
MCB = + 20 +
4 EI
(θ B + 2θ C )
4
2 EI ⎛ 3Δ ⎞
MCD = - 30 + ⎜ 2θ C + 0 + ⎟
6 ⎝ 6 ⎠
-76-
2 EI ⎛ 3Δ ⎞
MDC = + 30 + ⎜0 + θC + ⎟
6 ⎝ 6 ⎠
The unknown displacements θA, θB, θC and Δ can be obtained from three moment
equilibrium equations of joints and one shear equilibrium equation.
HA + HD – 60 = 0 ------- (4)
M BA
HA =
3
M CD + M DC + 60 x 3
HD =
6
2EI ⎛ 3Δ ⎞
⎜ 2θ A + θ B + ⎟= 0
3 ⎝ 3 ⎠
θB Δ
θA = - − ------------------ (1)
2 2
2EI
(θ A + 2θ B + Δ ) - 20 + EI (2θB + θC) = 0
3
60
2θA + 4θB + 2Δ + 6θB + 3θC =
EI
60
3θC + 10θB + 2θA + 2Δ = --- (2)
EI
EI ⎛ Δ⎞
20 + EI (θB + 2θC) – 30 + ⎜ 2θ C + ⎟ = 0
3 ⎝ 2⎠
-77-
60
6θB + 12θC + 4θC + Δ = +
EI
+ 60
6θB + 16θC + Δ = --------- (3)
EI
4EI
(θ A + 2θ B + Δ ) − 30 + EI ⎛⎜ 2θ C + Δ ⎞⎟ + 30 + EI ⎛⎜θ C + Δ ⎞⎟ = 180
3 3 ⎝ 2⎠ 3 ⎝ 2⎠
540
4θA + 8θB + 3θC + 5Δ = ----- (4)
EI
60
3θC + 9θB + Δ = --------- (5)
EI
540
6θB + 3θC + 3Δ = ----- (6)
EI
180
or 2θB + θC + Δ = ----- (6)
EI
120
7θB + 2θC = ----------- (7)
EI
120
4θB + 15θC = - ----------- (8)
EI
360 3.711
θC = - =−
97EI EI
-78-
Now member end moments are computed from slope deflection equations.
⎡ 2 x (- 16.082) 3.711⎤
MBC = - 20 + EI ⎢ − = - 55.875 kNm
⎣ EI EI ⎥⎦
⎡ - 16.082 2 x 3.711⎤
MCB = 20 + EI ⎢ − = - 3.504 kNm
⎣ EI EI ⎥⎦
EI ⎡ − 2 x 3.711 215.875 ⎤
MCD = - 30 + + = + 3.505 kNm
3 ⎢⎣ EI 2EI ⎥⎦
EI ⎡ − 3.711 215.875 ⎤
MDC = 30 + + = 64.742 kNm
3 ⎢⎣ EI 2EI ⎥⎦
Shears and axial forces at the ends of members are computed from free body diagrams of
members shown in Figure 16.9.
Member AB
55.876
HA = HB = = 18.625 kN
3
-79-
Member BC
HC = - HB = - 18.625 kN
VB =
40
+
(55.875 + 3.504) = 34.845 kN
2 4
40 (55.875 + 3.504)
VC = - = 5.155 kN
2 4
Member CD
VD = VC = 5.155 kN
10 x 6 x 3 + 64.742 + 3.505
HD = = 41.375 kN
6
Member AB
VA = VB = 34.845 kN
VA + VD = 40 kN, HA + HD = 60 kN (check)
∝E = R – r
∝ E = External indeterminacy
∝I = m – (2j – r)
∝I = Internal indeterminacy
-80-
∝S = ∝ E + ∝I = (m + r – 2j)
∝S = (m + R – 3j)
∝I = m – (3j – r)
X = Force in member which has been removed from redundant truss due to external loads
acting on indeterminate truss as shown in figure 17.3.
Fi = Total force in ith member due to external loads and the forces X corresponding to
removed member.
Fi = Fi + ki X
U=∑
m
(Fi + k i X )2 L i
i =1 2A i E
-81-
∂U
= Displacement in the direction of X
∂X
E = modulus of elasticity
XL o
ΔLo =
AoE
∂U
= - ΔLo
∂X
m
(Fi + k i X )k i L i − XL o
∑
i =1 AiE
=
AoE
∑Fk L
i =1
i i i / AiE
X=-
⎛ k i2 L i
m
L ⎞
⎜⎜ ∑ + o ⎟⎟
⎝ i=1 A i E A o E ⎠
U=∑
m
(Fi + k i X )2 L i
i =1 2A i E
∂U
=0
∂X
-82-
m
(Fi + k i X )k i Li
∑
i =1 AiE
=0
m
Fi k i L i
∑AE
X = - i=m1 2i
K i Li
∑i =1 A i E
The truss is made determinate by removing redundant members and reactions. The
determinate truss is analyzed for forces in members due to given external loads and due
to unit redundant forces X1, …. Xj, …., Xn.
n
Fi = Fi + ∑k
j=1
ij Xj
The strain energy of the truss including strain energy of redundant members is given by:
2
⎛ n ⎞
⎜ Fi + ∑ k ij X j ⎟ L i
m ⎜ ⎟ X 2j L j
U=∑⎝ ⎠
j=1
+
i =1 2A i E 2A j E
∂U ∂U ∂U ∂U
= = ….. = ……… = =0
∂X1 ∂X 2 ∂X j ∂X n
⎛ n ⎞
⎜ F + ∑ k ij X j ⎟k i L i
m ⎜ ⎟
∂U
= ∑⎝ ⎠
j=1 X jL j
+ =0
∂X j i =1 AiE A jE
-83-
For unknown forces corresponding to reactions Lj is set equal to zero. The partial
derivatives with respect to all the n-unknowns give a system of simultaneous equations,
which is solved for unknowns and the forces in indeterminate truss.
No forces are generated in members of determinate truss due to lack of fit but in case of
indeterminate truss the forces are generated in members when the members with lack of
fit are forced into position. Consider an indeterminate truss ABCD in which member BC
is shorter than its exact length by amount δ as shown in figure 17.8. The member will
come in tension when fitted.
System is made determinate by removing members having lack of fit. Unit forces X = 1
are applied as shown in figure 17.9.
Fi = kiX
U=∑
m
( k i X )2 L i
i =1 AiE
∂U
δ= = movement in the direction of force
∂X
m
k i2 XL i − XL o
∑
i =1 A i E
=
AoE
+δ
δ
X=
⎛ m k Li 2
L ⎞
⎜⎜ ∑ i
+ o ⎟⎟
⎝ i=1 A i E A o E ⎠
If member with lack of fit is longer than Lo, δ will be taken negative.
-84-
δLo = ∝Lo ΔT
U=∑
m
( k i X )2 L i
i =1 2A i E
∂U
= movement in the direction of force X
∂X
m
k i2 XL i XL o XL o
∑
i =1 A i E
= δ Lo −
AoE
= L o ∝ ΔT -
AoE
L o ∝ ΔT
X=
⎛ m k i2 L i L ⎞
⎜⎜ ∑ + o ⎟⎟
⎝ i=1 A i E A o E ⎠
The method discussed above is applicable to plane as well as space trusses and pertains to
force method.
Numerical examples
Example 17.1.
For the cantilever truss shown in figure 17.11 find the forces in members. The axial
rigidity AE of all the members is same.
-85-
There are three members meeting at joint A and two equations of equilibrium, hence,
static indeterminacy is 1. Let AC be redundant member. It is removed to obtain
determinate system as shown in figure 17.12.
Joint A
3 2 3
F1 = F3 = F3
2 2
⎛ 3 1⎞
⎜ ⎟
⎜ 2 + 2 ⎟ F3 = 1000
⎝ ⎠
2000
F3 =
( 3 +1 )
= 732.051 (comp)
F1 = 896.576 (Tension)
Joint A
k3 = 2 k1
k1 3
+ 2 k1 = 1
2 2
(
k1 1 + 3 = 2 )
-86-
2 1.4142
=
k1 =
(1 + 3 ) 2.7321
= 0.518 (comp)
2
k3 =
(
1+ 3 )
= 0.732 (comp)
Fi k i L i
∑
i =1,3 A i E F1k 1L1 + F3 k 3 L 3
X=- = -
⎛ k Li 2
L ⎞ k 12 L1 + k 32 L 3 + L o
⎜⎜ ∑ i
+ o ⎟⎟
⎝ i=1,3 A i E A o E ⎠
X=-
(- 656.798 + 618.759) = 38.039 = 19.039 kN
(0.379 + 0.619 + 1) 1.998
F2 = X = 19.039 kN
Check
x – direction
3 1
- 19.039 + 745.988 x - 886.714 x = - 19.039 + 646.045 – 627.001 = 0.004 (ok)
2 2
-87-
y-direction
1 1
886.714 x + 745.988 x - 1000 = 627.001 + 372.994 – 1000 = 0 (ok)
2 2
It can be seen that there are small computational errors as the forces are not summing up
to zero exactly. However, residuals are small and neglected.
Example 17.2.
For internally indeterminate truss shown in figure 17.15 determine forces in all members.
E = 200 kN/mm2
Member BC is removed to make the truss determinate as shown in figure 17.16. Forces
in members F are now determined.
RAV = RD = 7.5 kN
RAH = 10 kN
Joint C
F4 Cos θ = 10
5 50
F4 = 10 x = = 12.5 kN (Tension)
4 4
50 3 15
F2 = F4 Sin θ = x = = 7.5 kN (comp)
4 5 2
Joint B
F1 = F5 = 0
-88-
Joint A
50 4
F3 = 10 – F4 Cos θ = 10 - x =0
4 5
Now unit actions are applied at B and C as shown in figure 17.17 and forces are
determined in members.
Joint B
5
k5 Cos θ = 1, k5 = = 1.25 kN (comp)
4
5 3 3
k1 = k5 Sin θ = x = = 0.75 kN (Tension)
4 5 4
Joint A
3
k4 Sin θ = k1 =
4
5 3 5
k4 = x = kN (comp)
3 4 4
5 4
k3 = k4 Cos θ = x = 1 kN (Tension)
4 5
Joint C
5 3 3
K2 = k4 Sin θ = x = = 0.75 kN (Tension)
4 5 4
FkL FkL
∑ AE A
−∑
X=- =
k 2L Lo k 2L Lo
∑ AE + A E ∑ A + A
o o
-89-
L o 4000
= = 4.44444
Ao 900
− 223.437
X=- = + 4.17 kN
(49.1318 + 4.4444)
18 APPROXIMATE METHODS OF ANALYSIS OF
STATICALLY INDETERMINATE STRUCTURES
Even the so called exact analysis of structures is approximate in the broad sense because
every analysis is based on assumptions with respect to geometry and material behavior of
structure. However, the approximate analysis for statically indeterminate structures is
based on assumptions, which render the indeterminate system to a determinate system,
which can be analyzed exactly by use of equations of equilibrium. The advantage is that
analysis of determinate systems does not depend on elastic properties of members where
as analysis of statically indeterminate systems depends on elastic properties of members.
Approximate analysis is also required for assumption of initial elastic properties of
members in order to carry out design of indeterminate structure. Approximate methods
are also adopted under situation of configuration complexity of structure for which exact
method of analysis is not available. Where exact analysis requires prohibitive time and
cost it should be abandoned and approximate method adopted which generally gives
conservative design. The first step in approximate analysis of statically indeterminate
structure is to determine its statical indeterminacy and then introduce equal number of
releases to make it determinate.
-90-
For the portal shown in figure 18.2 the degree of static indeterminacy is 3. Hence,
three assumptions are required. The points of contra flexure occur at mid height of
columns. Hence, two hinges will be introduced at mid height as shown in Figure 18.3.
The third assumption will be made as equal resistance to shear by the columns. Now
analysis is possible by equations of equilibrium and is self explanatory from figure
18.4.
Here also, the hinges are introduced at mid height of columns and it is assumed that
columns carry same horizontal shear as shown in figures 18.5 and 18.6.
⎛H ⎞
P ⎜ + h⎟
Vertical reaction at hinges = ⎝
2 ⎠
L
-91-
Tower is a space truss. For approximate analysis it is necessary that each face of
tower should be a plane. Each face is separated out as a plane truss as shown in figure
18.8 and equivalent inplane loads are shown. The indeterminate plane truss is
replaced by two determinate plane trusses as shown figures 18.9 and 18.10 and the
loads are apportioned on these trusses. Final result of forces in the members is
obtained by superposing the results of individual trusses.
The analysis is carried out separately for vertical and lateral loads. For lateral loads
following methods are used.
(b) A point of inflection occurs in the girder at one-tenth point measured along span from
each end. The hinges are introduced at these points as shown in figure 18.12.
(c) The girder is analyzed for positive BM as a simply supported beam of span 0.8L and
the vertical reactions are transferred to cantilever parts of span 0.1L as shown in
Figure 18.13. The negative moments at the ends of beam are computed as the
negative moments at the cantilever clamps as shown in figure 18.14. The vertical
reactions of clamps are transferred to the columns as axial loads.
(d) At the beam column junction the clamp bending moment is distributed between the
columns above and below in proportion to their stiffness. A typical girder AB is
considered as shown in figure 18.15.
ω ( 0.8 L )2
Maximum positive BM at centre = = 0.08 ωL2.
8
Maximum reaction at hinge = 0.4ωL.
-92-
K2
MAU = MA
(K1 + K 2 )
K1
MAL = MA
(K1 + K 2 )
If at the joint their is beam on the other side also, its analysis can also be carried out
similarly and forces in column superposed.
Portal Method
(a) The point of contra flexure in each girder is located at its mid span point.
(b) The point of contra flexure in a column is located at mid height of each column.
(c) The total horizontal external shear on each storey is distributed between the column
of that storey so that each interior column carries twice as much shear as the exterior
column.
Shown in figure 18.16 is a plane frame with columns C1 to C6 and girders g1 to g4. It is
subjected to lateral loads P and Q. The hinges are assumed at mid points of all the
members and the horizontal shears are distributed in columns in accordance with the
rules as shown in figure 18.17.
-93-
BM in column
⎛ P + Q ⎞ H1 (P + Q ) H1
BM in C1 = ⎜ ⎟ =
⎝ 2 ⎠ 2 4
BM in C2 =
(PH 2 )
4
BM in girder
The girder and column moments act in apposite directions at their junction in accordance
with their elastic line or deflection curve as shown in figure 18.18.
At joint E:
Since column end moments are now known the beam end moments can be computed. At
end F of g2:
Thus proceeding across the girder of first level the girder moments are determined. The
girder moments at other level can also be determined in the same manner.
SF in girder
Considering free body of g2, as shown in figure 18.19 the shear at ends is obtained
2M g2
Vg2 =
L2
-94-
AF in columns
Axial forces in columns can be obtained by summing up from the top of the column, the
shears transferred to the column by the girders connected to this column.
AF in girder
Axial forces in girders are obtained by summing up the horizontal forces from one end of
the girder to the other end.
Cantilever method
1) As in portal method the points of contra flexure are assumed at mid points of all
beams and columns as shown in figure 18.20.
x=
(3L1 + 2L 2 + L 3 )
4
Taking a free body above the horizontal section through mid points of columns of first
storey the expression for moment about mid point of column DH is set equal to zero.
-95-
⎛ H ⎞ ⎛ H1 ⎞ F
P ⎜ H2 + 1 ⎟ + Q ⎜ ⎟ - FAE (L1 + L2 + L3) - AE (x – L1) (L2 + L3)
⎝ 2 ⎠ ⎝ 2 ⎠ x
FAE
+ (L1 + L2 - x) L3 = 0
x
This gives the value of FAE, hence, forces in columns of first storey are computed. The
forces in columns of second storey are computed in the same manner by taking moment
about mid point of column HL of the free body above the horizontal section through mid
points of columns of second storey.
Girder shears
The girder shears can be obtained from the column axial forces at the various joints as
shown in figure 18.22.
Joint I
VIJ = FIE
Joint E
Girder moments
The moment at each end of girder equals the shear in that girder multiplied by half the
length of that girder as shown in figure 18.23. For example:
L1
MEF = VEF
2
Column moments
Column moments are determined by beginning at the top of each column stack and
working progressively down to base as shown in figure 18.24.
Joint J
L1 L
MJF = VIJ + VJK 2
2 2
-96-
Since there is a point of contra flexure at the center of FJ, MFJ will be equal to MJF.
MFJ = MJF
MBF = MFB
Proceeding in this manner moments, shears and axial forces in all members are
determined.
This method is more accurate compared to portal and cantilever methods. The
assumptions in factor method are based on elastic action of structure hence results of this
method correspond to approximate slope-deflection analysis. The method is not discussed
here.
Readers are welcome to point out any errors noticed in the lecture course material
so that corrections can be incorporated. Any suggestions are welcome by the author.