Finite Element Analysis Notes On One Dimensional Structural Analysis-Libre
Finite Element Analysis Notes On One Dimensional Structural Analysis-Libre
Now
a0=x1+x2/2
a1= x2-x1/2
U = Nq
6
U = Nq
Where N is the shape function matrix and q is displacement matrix.
Once the displacement is known its derivative gives strain and
corresponding stress can be determined as follows.
Third step in FEM is finding out stiffness matrix from the above
equation we have the value of K as
But
Fourth step is assembly and the size of the assembly matrix is given
by number of nodes X degrees of freedom, for the present example that
has four nodes and one degree of freedom at each node hence size of
the assembly matrix is 4 X 4. At first determine the stiffness matrix of
each element say k1, k2 and k3 as
The given system is modeled as three elements and four nodes we have
three stiffness matrices.
10
Given system is fixed at one end and force is applied at other end.
Since node 1 is fixed displacement at node 1 will be zero, so set q1 =0.
And node 2, node 3 and node 4 are free to move hence there will be
displacement that has to be determined. But in the load vector because
of fixed node 1 there will reaction force say R1. Now replace F1 to R1
and also at node 3 force P is applied hence replace F3 to P. Rest of the
terms are zero.
11
12
13
14
15
K=
..
+ QNKN1Q1+QNKN2Q2+. +QNKNNQN)
- (Q1F1 + Q2F2++QNFN)
Substituting Q1 = a1 we have
= (a1K11a1 +a1K12Q2+..+ a1K1NQN
+ Q2K21a1+Q2K22Q2+. + Q2K2NQN
..
+ QNKN1a1+QNKN2Q2+. +QNKNNQN)
- (a1F1 + Q2F2++QNFN)
Extremizing the potential energy
ie
d/dQi = 0 gives
Where i = 2, 3...N
K22Q2+K23Q3+. + K2NQN = F2 K21a1
K32Q2+K33Q3+. + K3NQN = F3 K31a1
Q2
Q3
.
.
QN
F2-K21a1
F3-K31a1
=
FN-KN1a1
18
We may also choose 105 &106 but 104 found more satisfactory on most
of the computers.
Because of the spring the stiffness matrix has to be modified ie the
large number c gets added to the first diagonal element of K and Ca 1
gets added to F1 term on load vector. That results in.
A reaction force at node 1 equals the force exerted by the spring on the
system which is given by
19
To solve the system again the seven steps of FEM has to be followed,
first 2 steps contain modeling and discretization. this result in
20
Similarly
21
22
Solution:
23
24
25
Solution:
26
27
28
29
We know that
Therefore
30
Therefore
32
33
Solution:
34
35
36
37
ie
38
Or
39
Graphs show the variation of shape functions within the element .The
shape function N1 is equal to 1 at node 1 and zero at rest other nodes (2
and 3). N2 equal to 1 at node 2 and zero at rest other nodes(1 and 3) and
N3 equal to 1 at node 3 and zero at rest other nodes(1 and 2)
40
By chain rule
Now
41
Therefore
Stiffness matrix:
We know the stiffness matrix equation
For an element
42
Where
and BT
43
Body force term & surface force term can be derived as same as 2
noded bar element and for quadratic element we have
Body force:
This amount of body force and surface force will be distributed at three
nodes as the element as 3 equally spaced nodes.
44
Solution:
45
46
47
48
49
ANALYSIS OF TRUSSES
A Truss is a two force members made up of bars that are
connected at the ends by joints. Every stress element is in either tension
or compression. Trusses can be classified as plane truss and space truss.
Plane truss is one where the plane of the structure remain in
plane even after the application of loads
While space truss plane will not be in a same plane
Fig shows 2d truss structure and each node has two degrees of freedom.
The only difference between bar element and truss element is that in
bars both local and global coordinate systems are same where in truss
these are different.
Consider one truss element as shown that has nodes 1 and 2 .The
coordinate system that passes along the element (x l axis) is called
local coordinate and X-Y system is called as global coordinate
system. After the loads applied let the element takes new position
say locally node 1 has displaced by an amount q 1l and node2 has
moved by an amount equal to q2l.As each node has 2 dof in
global coordinate system .let node 1 has displacements q1 and q2
along x and y axis respectively similarly q3 and q4 at node 2.
Resolving the components q1, q2, q3 and q4 along the bar we get two
equations as
51
Or
U = qTKq
For a truss element we can write
U = qlT K ql
Where ql = L q and q1T = LT qT
52
Therefore
U = qlT K ql
Taking the product of all these matrix we have stiffness matrix for truss
element which is given as
53
But strain = B ql
and ql = T q
Therefore
54
55
3
2
= 901.3 mm
Similarly calculate all the parameters for element 2 and tabulate
56
Element 1 has displacements q1, q2, q3, q4. Hence numbering scheme
for the first stiffness matrix (K1) as 1 2 3 4 similarly for K2 3 4 5 & 6 as
shown above.
Global stiffness matrix: the structure has 3 nodes at each node 3 dof
hence size of global stiffness matrix will be 3 X 2 = 6
ie 6 X 6
57
58
59
Global stiffness matrix: the structure has 4 nodes at each node 3 dof
hence size of global stiffness matrix will be 4 X 2 = 8
ie 8 X 8
60
61
Beam element
Beam is a structural member which is acted upon by a system of
external loads perpendicular to axis which causes bending that is
deformation of bar produced by perpendicular load as well as force
couples acting in a plane. Beams are the most common type of
structural component, particularly in Civil and Mechanical
Engineering. A beam is a bar-like structural member whose primary
function is to support transverse loading and carry it to the supports
62
63
But
Therefore strain energy for an element is given by
64
65
Therefore
Similarly we can derive
66
Stiffness matrix:
Once the shape functions are derived we can write the equation
of the form
But
ie
67
68
Varying load
69
Solution:
Lets model the given system as 2 elements 3 nodes finite
element model each node having 2 dof. For each element determine
stiffness matrix.
70
For element 2 that has UDL its equivalent load and moment are
represented as
ie
71
= -0.089mm
72
73
ie
74
And also we have external point load applied at node 3, it gets added to
F5 term with negative sign since it is acting downwards. Now F
becomes,
From KQ=F
75
Reference:
1) Finite elements in Engineering,
ChandrupatlaTR&A.D Belegundupears edition
2) Finite element Analysis
H.V.Laxminarayana, universities press
3) Finite element Analysis
C.S.Krishnamurthy, Tata McGraw, New Delhi
4) Finite element Analysis
P.seshu , prentice hall of India, New Delhi
5) Finite element Method
J.N.Reddy, Tata McGraw
6) Finite element Analysis
A. J. Baker & D W Pepper
76