Finite Element Analysis 17ME61

Module II
One-Dimensional Elements-Analysis of Bars and Trusses, Linear interpolation polynomials in terms of
local coordinate‟s for 1D, 2D elements. Higher order interpolation functions for 1D quadratic and cubic
elements in natural coordinates, Constant strain triangle, Four-Nodded Tetrahedral Element (TET 4), Eight-
Nodded Hexahedral Element (HEXA8), 2D iso-parametric element, Lagrange interpolation functions,
Numerical integration: Gaussian quadrature one point, two point formulae, 2D integrals. Fore terms: Body
force, traction force and point loads,
Numerical Problems: Solution for displacement, stress and strain in 1D straight bars, stepped bars and
tapered bars using elimination approach and penalty approach, Analysis of trusses
08 Hours

ONE-DIMENSIONAL ELEMENTS-ANALYSIS OF BARS AND TRUSSES

Finite element formulation of 1-D Bar Element: The one-dimensional elements are used to represent axial
bar, springs, truss, beam like member for the finite element analysis purpose. finite elements formulation of
an element means derivations of elemental properties such as shape function, stiffness matrix, strain
displacement matrix, strain matrix and load vector.

Element formulation of 1-D bar element

The line elements are used to represent spring, truss, beam like members for the finite element analysis
purpose. Such elements are quite useful in analyzing truss, cable and frame structures. Such structures tend
to be well defined in terms of the number and type of elements used. For example, to represent a truss
member, a two node linear element is sufficient to get accurate results. However, three node line elements
will be more suitable in case of analysis of cable structure to capture the nonlinear effects.
Element Stiffness Matrix The stiffness matrix of a structural system can be derived by various methods like
variationalprinciple, Galerkin method etc. The derivation of an element stiffness matrix has already been
discussed in earlier lecture. The stiffness matrix is an inherent property of the structure. Element stiffness is
obtained with respect to its axes and then transformed this stiffness to structure axes.
The properties of stiffness matrix are as follows:
 Stiffness matrix issymmetric and square.
 In stiffness matrix, all diagonal elements are positive.
 Stiffness matrix is positive definite
Global Stiffness Matrix A structural system is an assemblage of number of elements. These elements are
interconnected together to form the whole structure. Therefore, the element stiffness of all the
elementsarefirst need to be calculated and then assembled together in systematic manner. It may be noted
that the stiffness at a joint is obtained by adding the stiffness of all elements meeting at that joint. To start

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
with, the degrees of freedom of the structure are numberedfirst. This numbering will start from 1 to n where
n is the total degrees of freedom. These numberings are referred to as degrees of freedom corresponding to
global degrees of freedom. The element stiffness matrix of each element should be placed in their proper
position in the overall stiffness matrix. The following steps may be performed to calculate the global
stiffness matrix of the whole structure. a. Initialize global stiffness matrix K as zero. The size of global
stiffness matrix will be equal to the total degrees of freedom of the structure. b. Compute individual element
properties and calculate local stiffness matrix k of that element. c. Add local stiffness matrixk to global
stiffness matrix K using proper locations d. Repeat the Step b. and c. till all local stiffness matrices are
placed globally.

Derivation of Shape function for 1-D linear bar element in local co-ordinate system:

Figure shows one dimensional bar element having two nodes 1 and 2, at each need one degree of freedom,
Cartesian coordinates are x1 and x2.. Let P be the point inside the element having local coordinates (L1,, L2)
and x be the Cartesian coordinates.

Since local coordinates are dimensionless numbers and whose value never exceeds unity, we define local
coordinates are
L1 + L2 =1
L1X1 + L2 X2 =x

In matrix form
0 12 3= 2 3

2 3 0 1 2 3

( )
[ ]2 3

2 3 ( )
0 1

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
The shape function and

Where ( )

Thus,
Variation of L1 and L2 is as shown in figure

Derivation of strain and [B] – Strain displacement matrix:

Figure shows one dimensional bar element having two nodes 1 and 2 whose Cartesian coordinates are x1 and
[Link] Cartesian coordinate and in natural coordinates. Position of nodes are

We know that, strain in a element is given by

Element strain in terms of partial derivatives of the displacement wrt is given by

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Strain matrix,

, -2 3

Strain matrix it can also be written as

Strain matrix, , -* +

Stress matrix

From Hooke‟s Law

, -* +

Derivation of stiffness matrix for 1-D Element:

we know that strain energy for an element is given by

S.E = ∫

For one dimensional bar element

S.E = ∫

Where , A is cross sectional area of the element,

, -* + , -* +

Thus ,substituting values of

We get
SE= ∫ , -* +

SE = {u}[K] {u}

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
Thus, Elemental stiffness matrix
[KE] = 0 1

Finite element formulation of 2-D Element:

Figure shows the triangular elements with difference numbers of nodes used for solving two dimensional
solid members.

Derivation of shape function of CST element in Local coordinates:

Figure shows the two dimensional triangular element having three nodes 1,2 and 3 .let P be the any point
inside the element having three coordinates.

Let
; and

But

Where,

A be the area of triangle 123

A2 be the area of triangle P23
A3 be the area of triangle P13
A4 be the area of triangle P12
The area of the triangles can be written using Cartesian coordinates considering x,y as coordinates of an
arbitarary point P inside or on the boundaries of the element:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

[ ]

[ ]

[ ]

[ ]

Since L1 are defined in terms of area, they are known as area coordinates. therefore they can be defined as

In matrix form

[ ]{ }= { }

{ } [ ] { }

Higher order interpolation functions for 1D quadratic and cubic elements in natural coordinates:

Rectangular elements are suitable for modeling regular geometries. sometimes, it is used along with
triangular elements to represent an arbitrary geometry .the simplest element in the rectangular family is the
four node rectangle with sides parallel to x and y axis. Figure shows rectangular elements with varying
nodes representing linear, quadratic and cubic variation of function

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

LAGRANGE INTERPOLATION FUNCTION FOR LINEAR ELEMENT:

Consider the two node bar element as shown in figue. let us consider the natural coordinate of the centre of
the element as 0, and the natural coordinate of the nodes 1 and 2 are -1 and +1 respectively. therefore the
natural coordinate at any point x can be represented by

( )

The shape function of the 1-D bar element is found by using Lagrangian polynomial in global and Cartesian
coordinates.

( ) ∏

For two noded bar element n=2.

Hence when k=1
L1(x) = N1(x) =

Hence when k=1

L2(x) = N2(x) =
In natural coordinates

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

( ) ∏

Similarly for two noded bar element n=2

Hence when k=1
L1 ( ) ( )
Shape function at node 1
N1(x) = ( )

Shape function at node 2

N1(x) = ( )

The shape function of 1-D bar element is as shown in figure

Shape function for 2-D elements-quadrilateral /rectangular element:

We can derive lagrange interpolation function for two or three dimensional element from one dimensional
element as discussed previously. those elements whose shape functions are derived from the product of 1-D
Lagrange interpolation function called Lagrange elements.

Let us consider the general quadrilateral element mapped into natural coordinates as shown in figure:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

To find shape function N1 :considering element as shown in figure

N1( ) ( ) ( )

( ) ( )=

Shape function at node 1:

N1= ( )( )

Shape function at node 2:

N1= ( )( )

Shape function at node 3:

N1= ( )( )
Shape function at node 4:

N1= ( )( )

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
Iso-parametric, Sub-parametric and Super-parametric Elements:

Iso-parametric elements: The elements for which the geometry and the displacements are represented by
the same shape functions are called iso–parametric elements.

Example: Suppose N1, N2, N3 and N4, represents the shape functions of a 4 noded quadrilateral elements.
The geometry of the element is
x= +
y= +
and the displacement is
u= +
v= +

Elements are represented by same shape functions N1, N2, N3 and N4. Hence, the 4 noded quadrilateral
element is an iso-parametric element.
It may be noted that equations (1) and (2) are of the same order and contain same number of terms. i.e., all
the nodes (terms in displacement polynomial) of the element are used to represent geometry and
displacement.

Sub-parametric elements: The elements (including higher order elements) for which the geometry is
represented by lesser terms and of lesser order when compared to the representation of displacement.

Example: For a beam element with 2 nodes, the shape functions N1 and N2, used for linear 1D element are
utilized for representing geometry.

But for representing displacements (vertical deflection and slope) some other shape functions H1, H2, H3
and H4 are used which are called Hermite cubic shape functions. These shape functions conveniently
represent deflection and slope at a node. Note that deflection is perpendicular to beam axis. It can also be
noted that H1, H2, H3 and H4 are 4 terms and of higher order (cubic). Where as N1 and N2 are 2 terms and
of lower order (linear).

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
Another example is the nine noded quadrilateral element. When curved edges are involved, the shape
functions for all nodes are used to represent geometry. When straight (linear) edges are involved, only 4
corner nodes shape functions are used to represent geometry which is sub-parametric representation.

Super-parametric elements: These are elements which are different from the iso-parametric and sub-
parametric elements. i.e., higher order functions with larger number of terms are required to represent
geometry. This means, representation of displacement is done by lesser number of terms and of lower order.
Note: All iso-parametric and sub-parametric elements satisfy the convergence requirements, where as only a few
super-parametric elements are able to satisfy convergence requirements

Types of Boundary Conditions

In structural problems we normally come across two types of boundary conditions
1. Specified displacement boundary condition
2. Multi point constraint boundary condition

If the boundary conditions are of the type Q1 = a1, Q2 = a2, Q3 = a3…………Qr = ar such are called as
specified displacement boundary condition. Here the specified displacement may be zero or non zero
values. If the specified displacement is zero it is called Homogeneous Boundary Condition. If it is non zero
it is called Non-Homogeneous Boundary Condition.

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Numerical Integration: In FEM, evaluation of certain matrices like the element stiffness matrix [K] e,
element body force vector {B}e, element traction force vector {T}e, element strain vector {E}e, etc. involves
the evaluation of integrals of matrices which are functions of two or more variables.
For example, the element stiffness matrix [K]e for a quadrilateral element is given as:

, - ∬ , -, - , - , -

Taking another example, element traction force vector {T}e for a quadrilateral element is given as:

, -= t 0∫ , - 1* +

Both the above integrals (a) & (b) may be evaluated by performing the exact integration and then applying
the limits. But this is a laborious task and not suited to be implemented as a computer program. Also, the
requirement in FEM is to evaluate the element properties at fixed point/s inside each and every discretised
element. Hence, the Numerical Integration procedure like “Gauss Quadrature” is applied.

Gauss Quadrature Approach:

Consider the problem of evaluating the one-dimensional definite integral of the form:

I=∫ ( )

By Gauss Quadrature method, the solution to this integral „I‟ can be given by a „n‟ point approximation as:

I= +

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
Theoretically, if „n‟ is set to infinity, Gauss Quadrature method will give an exact solution to the integral.
But, in actual practice, we limit „n‟ to some number which will result in an approximate solution to the
integral.
One Dimensional Integral: If the integral is with respect to only one variable, it is called One Dimensional
Integral.

Example I = ∫ ( )
I is with respect to single variable ζ

Two Dimensional Integral: If the integral is with respect to two variables, it is called Two Dimensional
Integral.

Example I = ∬ ( )
I is with respect to two variables ζ and η

One Point Formula or One Point Integration Rule:

The Gauss Quadrature method where „n‟ is set to one, i.e., n = 1, is called One Point Formula or One Point
Integration Rule.
For 1D integrals, only one weight and one sampling point is used.
For 2D integrals, also one weight and one sampling point is used and is usually
denoted as 1×1 Point Integration Rule

Two Point Formula or Two Point Integration Rule:

The Gauss Quadrature method where „n‟ is set to two, i.e., n = 2, is called Two Point
Formula or Two Point Integration Rule.
For 1D integrals, two weights and two sampling points are used.
For 2D integrals, four weights and four sampling points are used and is usually
denoted as 2×2 Point Integration Rule

Three Point Formula or Three Point Integration Rule:

The Gauss Quadrature method where „n‟ is set to three, i.e., n = 3, is called Three
Point Formula or Three Point Integration Rule.
For 1D integrals, three weights and three sampling points are used.
For 2D integrals, nine weights and nine sampling points are used and is usually
denoted as 3×3 Point Integration Rule

Condition for ‘n’ point integration rule:

Availaible at VTU HUB (Android App)
Finite Element Analysis 17ME61
„n‟ point integration rule can give exact solution to the integral only if the expression inside the integral is of
the degree (2n-1).
i.e., one-point integration rule (n = 1) can give exact answer to an integral only if the expression inside the
integral is of the degree (2n-1) = (2×1-1) = 1. (linear expression) Similarly, two-point integration rule (n = 2)
can give exact answer to an integral only if the expression inside the integral is of the degree (2n-1) = (2×2-
1) = 3. (cubic expression)

ANALYSIS OF BARS AND TRUSSES ONE DIMENSIONAL ANALYSIS

METHODS OF HANDLING BOUNDARY CONDITIONS

Elimination Method

We have two methods of handling boundary conditions namely Elimination method and penalty approach
method. Applying BC‟s is one of the vital role in FEM improper specification of boundary conditions leads
to erroneous results. Hence BC‟s need
to be accurately modelled. Elimination Method: let us consider the single boundary
conditions say Q1 = [Link] ÉÆ results in equilibrium equation.
Q = [Q1, Q2, Q3……….QN]T be the displacement vector
and
F = [F1, F2, F3…………FN] T be load vector

Say we have a global stiffness matrix [K] = [ ]

Now potential energy can written as

PHI= (Q1K11Q1 +Q1K12Q2+…..+ Q1K1NQN

+ Q2K21Q1+Q2K22Q2+………. + Q2K2NQN
…………………………………………
……………………………………… ..
- (Q1F1 + Q2F2+…………………+QNFN)

Substituting Q1 = a1 we have

PHI = . (a1K11a1 +a1K12Q2+…..+ a1K1NQN

+ Q2K21a1+Q2K22Q2+………. + Q2K2NQN

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
…………………………………………
……………………………………… ..
+QNKN1a1+QNKN2Q2+……. +QNKNNQN)
(a1F1 + Q2F2+…………………+QNFN)

Extremizing the potential energy

Where i = 2, 3...N
K22Q2+K23Q3+………. + K2NQN = F2 – K21a1 K32Q2+K33Q3+
………. + K3NQN = F3 – K31a1
………………………………………………
KN2Q2+KN3Q3+………. + KNNQN = FN – KN1a1
Now the N X N matrix reduces to N-1 x N-1 matrix as we know Q1=a1 ie first row and first column are
eliminated because of known Q1. Solving above matrix gives displacement components. Knowing the
displacement field corresponding stress can be
calculated using the relation Reaction forces at fixed end say at node1 is evaluated using the relation
R1= K11Q1+K12Q2+……………+K1NQN-F1

Penalty approach method: let us consider a system that is fixed at both the ends as shown

In penalty approach method the same system is modelled as a spring wherever there is a support and that
spring has large stiffness value as

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
Displacement at node 3. The displacement Q1 at node 1 will be approximately equal to a1, owing to the
relatively small resistance offered by the structure. Because of the spring addition at the support the strain
energy also comes into the picture of equation
.Therefore equation becomes
PHI = . QTKQ+ . C (Q1 –a1)2 - QTF
The choice of C can be done from stiffness matrix as
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 Ca1 gets added to F1 term on load vector.
That results in.

[ ][ ] =[ ]

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

QUESTION BANK With Solution

THEORY QUESTIONS:

1. Derive the shape function for a 3 node triangular element (DEC.09/JAN-10)

2. What are the different types of interpolation functions? explain briefly
3. Derive shape function for CST element and also the displacement matrix( MAY/JUNE.2010)
4. Derive shape function for a four noded quadrilateral element, in natural co-ordinates
( DEC2010)
5. Derive shape functions for a four noded quadrilateral element,in natural cordinates(DEC2010)
6. For a four noded quadrilateral element,derive an expression for jacobian matrix(DEC.2009-
JUNE.10)
7. Derive the Jacobian matrix for 2-D Triangular Element( JUNE –JULY 2005)

Numericals

1. For a triangular element shown in figure ,obtain the strain-displacement matrix’B’ and
determine the strains and nodal displacement.(DEC.08/JUNE.09)

Solution:
Strain displacement matrix for a triangular element is given by

A = Area of the element

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Strain matrix ,

1. Figure shows a one dimensional bar subject to an axial load. Taking it as single bar element,
determine: (a) Nodal displacement (b) Elemental stresses and (c) support reactions.

a. Finite elemental model

Data: Young’s Modulus, E = 200 GPa = 2x105 MPa, Area of element, A1 = 104 mm2, Length of the element, l1 = 1000
mm, load at node 2, F2 = 10kN = 10x103 N

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
b. Nodal displacement
i. Elemental matrices
Elemental stiffness matrix: , - 0 1

For element 1: , - 0 1 0 1

ii. Global stiffness matrix

Since there is only one element, global stiffness matri is same as elemental stiffness matrix.

, - , - 0 1 0 1

iii. Global nodal displacement vector

* + 2 3

iv. Global load vector

* + { } 2 3

v. Equilibrium condition
The equilibrium condition for the bar is given by,
, -* + * +
0 12 3 2 3

vi. Applying boundary conditions

Boundary conditions are: (fixed end), results in

The nodal displacement vector becomes

* + 2 3

c. Elemental stress
Stress in a 1D bar element is given by , -2 3, where i = e

, -2 3

d. Reactions at support
The reaction forces at the supports are given by: * + , -* + * +

{ } 0 12 3 2 3

{ } { }

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
2. Figure shows a one dimensional uniformly tapered bar subject to an axial load. Taking it as
two bar element, determine: (a) Nodal displacement (b) Elemental stresses and (c) support
reactions.

a. FE model
Fig

Data: at node 1 = at node 1 = .

b. Nodal displacement vector

Area at node 2,

Area of element 1 = ̅̅̅

Area of element 2 = ̅̅̅

i. Elemental stiffness matrix

Stiffness matrix of bar element is given by: , - 0 1

For element 1: , - 0 1 0 1

For element 2: , - 0 1 0 1

ii. Global stiffness matrix

, - , - , -

, - [ ]

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
iii. Global nodal displacement vector

* + { }

iv. Global load vector: * + { }

v. Equilibrium condition: , -* + * +

[ ]{ } { }

vi. Applying boundary conditions

0 12 3 2 3

The nodal displacement vector,

* + { }

c. Elemental stress
Stress in a 1D bar element is given by , -2 3, where i = e

For element 1; , -2 3

For element 2; , -{ }

d. Support reactions

…………………………………..

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

3. Figure shows a one dimensional step bar subject to an axial load. Taking it as three bar
element, determine: (a) Nodal displacement (b) Elemental stresses and (c) support reactions.

a. FE model
Fig

Data:

b. Nodal displacement vector

i. Elemental stiffness matrix
Stiffness matrix of bar element is given by: , - 0 1

For element 1: , - 0 1 0 1

For element 2: , - 0 1 0 1

For element 3: , - 0 1 0 1

ii. Global stiffness matrix

, - , - , - , -

, - [ ]

iii. Global nodal displacement vector

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

* + { }

iv. Global load vector: { } { }

v. Equilibrium condition: , -* + * +

[ ]{ } { }

vi. Applying boundary conditions

0 12 3 { }

The nodal displacement vector,

* + { }

c. Elemental stress
Stress in a 1D bar element is given by , -2 3, where i = e

Stress in the thickest section;

Element 1; , -2 3

d. Support reactions
* + , -* + * +

{ } [ ]{ } { }

{ } { }

Left support reaction,

Availaible at VTU HUB (Android App)
Finite Element Analysis 17ME61

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
Truss

4. For a two bar truss shown in the figure, Determine: (a) Nodal displacement and (b) Elemental
stresses. Take A=200 mm2 & E= 2×105 N/mm2.

a. Nodal displacement
i. Nodal coordinate data
Node No x y
1 900 0
2 500 300
3 0 300

ii. Elemental connectivity data

Element Initial node Final node Length of the Direction cosines
element (le) l m
1 1 2 500 -0.8 0.6
2 2 3 500 -1 0

Where, ; ;

√( ) ( )

iii. Elemental stiffness matrix

, - [ ]

For element 1;

, - [ ]

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

[ ]

For element 2;

, - [ ] [ ]

iv. Global stiffness matrix;

, - , - , -

, -

[ ]

v. Global force vector

* +

{ } { }

vi. Equilibrium equation

, -* + * +

[ ]{ } { }

vii. Llying boundary conditions;

The matrix reduces to; 0 12 3 2 3

Nodal displacement vector;

{ }

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
b. Stress in each element

For element 1; , -{ } , -{ }

For element 2; , -{ } , -{ }

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
5. For a three bar truss shown in the figure, Determine: (a) Nodal displacement and (b)
Elemental stresses. Take A1=1500 mm2 , A1= A3 = 2000 mm2 & E= 200 GPa

a. Nodal displacement
i. Nodal coordinate data
Node No x y

1 0 0

2 800 0

3 400 400

ii. Elemental connectivity data

Element Initial node Final node Length of the Direction cosines

element (le)
l m

1 1 2 800 1 0

2 2 3 565.68 -0.707 0.707

3 1 3 365.68 0.707 0.707

Where, ; ;

√( ) ( )

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
iii. Elemental stiffness matrix

, - [ ]

For element 1;

, - [ ]

[ ]

For element 2;

, - [ ] [ ]

For element 3;

, - [ ] [ ]

iv. Global stiffness matrix;

, - , - , - , -

, -

[ ]

v. Global force vector

* +

{ } { }

Where

vi. Equilibrium equation

, -* + * +

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

[ ]{ } { }

vii. Applying boundary conditions;

The matrix reduces to; [ ]{ } { }

Nodal displacement vector;

{ }

b. Stress in each element

For element 1; , -{ } , -{ }

For element 2;

, -{ } , -{ }

For element 3;

, -{ } , -{ }

2. Consider the Bar shown in figure,an axial load p is applied using penalty approach for
handling boundry conditions ,do the following .nodal displacement, elemental stress, reations.
(DEC.09/JAN.10,JUNE/JULY 2013)

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Solution:
(i) FE Model

The model of the stepped bar is as shown in figure

(ii) Elemental stiffness matrix

The stiffness matrix for the bar elemnt is given by

For Element 1

For Element 2

Global stiffness matrix

The global stiffness matrix for the stepped bar is given by

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Global load vector:

Equilibrium cindition:
The equilibrium condition for the stepped bar is given by

For element 1:

For element 2:

3. Figure shows a one dimensional bar subjected to an axial loading .taking it as a two bar
element, determine, Nodal displacement, Stress in Each element and reaction at the
support.(JUNE/JULY.2009)

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Solution:

Finite element FE model: Modeleing the given bar using the two bar each having two nodes and location of
nodes are as shown in the figure

Elemental matrix:
The stiffness matrix for the bar element is given by :

Element 1:

Element 2:

Global stiffness matrix for the bar given below:

Global load vector: The load vector is given by

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Substituting the corresponding values we get:

For element 1:

For element 2:

Reaction at support:

4. Figure shows a bar subjected to uniformly distributed load p,taking youngs modulus E =
70GPA , Area of cross-section A= 104mm2 .Determine the nodal displacement and stress in the
element.(MAY/JUNE.2009)

Solution:
(i) FE Model:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
The FE model of the bar is as shown in figure:

Elemental matrices:
The stiffness matrix for the bar element is given by:

For Element 1

Nodal displacement vector is given by:

Load vector due to distributed load is given by:

Equilibrium condition
Equilibrium condition or stiffness equation for continuous given by

Applying boundary conditions

For Element 1:

5. Using the direct stiffness method, determine the nodal displacement of stepped bar shown in
figure(JUNE/JULY.2009,MAY/JUNE.2009)

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Solution:
FE Model:
The FE model of the stepped bar is as shown in figure:

Element stiffness matrix:

The stiffness matrix for the bar element is given by:

For Element 1:

For Element 2:

Global stiffness matrix:

The Global stiffness matrix for the stepped bar is given by:

Global nodal displacement vector:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
The Nodal displacement vector for the stepped bar is given by:

Global load vector

The load vector for the stepped bar is given by:

Equilibrium condition:
The equilibrium condition for the stepped bar is given by:

Applying the boundary conditions:

Using elimination method of handling boundary condition,eliminating the corresponding row and column
wrt zero boundary condition

Therefore the nodal displacement vector:

6. Solve for stress in member of structure given below in figure, using penalty approach of
handling boundary condition(JUNE/JULY.2010,DEC.07/JAN.08)

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Solution:
FE Model
The FE model of the stepped bar is as shown in figure:

Element stiffness matrix

The stiffness matrix for the bar element is given by

For Element 1:

For Element 2:

Global stiffnes matrix:

The global stiffnes matrix for the stepped bar is given by

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Global load vector

The load vector for the stepped bar is given by:

Equilibrium condition
The equilibrium condition for the stepped bar is given by

Thus the equilibrium equation can be written as

Stress in each element is given by:

For element 1

For element 2.

Reaction forces:
The reaction forces at the support is given by:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

ASSIGNMENT PROBLEMS

1. Find the nodal displacement,stress in the thickest section and left support reaction for
structure shown in figure: (DEC.2007/JAN.2008)

Solution:
FE Model
The FE Model of the stepped bar is as shown in the figure:

For Element 1:

For Element 2:

For Element 3:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Global stiffness matrix

The global stiffness matrix for the stepped bar is given by

Global load vector

The load vector for the stepped bar is given by:

Equilibrium condition
The equilibrium condition for the stepped bar is given by:

Applying the boundary condition the stress in the thickest section is given by:

Support reaction

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

2. A bar is having uniform cross-sectional area of 300mm and subjected to a load p as shown in
the [Link] the displacement field,stress and support reaction in the [Link] two
element and use elimination method to handle boundary
conditions.(JUNE/[Link].2010/JAN.2010)

Solution:
FE Model
Modeling the given bar considering two bar element as shown in the figure:

Elemental stiffness matrix:

The stiffness matrix for the bar element is given by

For Element 1:

For Element 2:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Global stiffness
Global stiffness matrix for the bar element is given by

Global nodal displacement vector

The nodal displacement vector for the bar is given by

Global load vector

We can see gap between node 3 and support at the [Link] force at node 3 cannot be [Link]
load vector can be written as

Equilibrium condition the equilibrium condition for the bar is given by

Stresses in each element

The stress in each element is given by

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

For Element 1

For Element 2

Support reactions:

3. Solve the following system of simultaneous equation by Gaussian elimination method.

(DEC.2007/JAN.2008)

Solution:
The given system of equation can be expressed in matrix as

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
Step1: In this step x1 is eliminated from second row onwards using the general algorithm:

Step 2:
In this step x2 in eliminated from elimination equation3 on words using equation

Using backward substitution

4. For the two bar shown in figure,determine the nodal displacement,stresses in each element and
reaction at the support. .(JUNE/[Link].2010/JAN.2010)

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Solution:

Nodal coordinate data:

Element connectivity data:

Element stiffness matrix:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

For element 1:

For element 2:

Global stiffness matrix:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

On substitution the value of k,q and F we get

Stress in each element:

For element 1:

For element 2:

Support reaction:
The reaction forces at the supports is given by

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

5. For the two truss shown in figure,deteremine the noda; displacement and the stress in each
member also find the support reactions: (JUNE/JULY.2010,DEC.07/JAN.08)

Solution:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Nodal coordinate data:

Element connecting data:

Element stiffness matrix:

For element 1:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

For element 2:

Global stiffness matrix:

Equilibrium rquation for the entire truss:

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61

Stresses in each element:

Element 1:

Element 2:

Support reactions:
The reaction forces at the supports is given by:

Availaible at VTU HUB (Android App)

 Finite Element Analysis 17ME61

1. Evaluate the integral using 1 point and 2 point integral rule and check the solution by exact
integration. (MAY/JUNE2010’DEC2009)

I=∫ 0 1

Solution: The given integral is a 1D integral. (Integral is with respect to only one variable x)

By One-Point Integration Rule: f (x) =

One-point rule for 1D integral is: I = ∫ ( ) = W1 f (x) = 2 f (0)

By Two-Point Integration Rule:

Two point rule for 1D integral is:

I=∫ ( ) = W1 f (x) + W1 f (x) = (1) f . / + (1) f . /

√ √

I = (1) [ √ . / ] + (1) [ √ . / ]
√ √
√ √

By exact integration:

∫ 0 1 = 0 ( )1

= ( ) + ( 1+1) + (ln 3- ln1) = 8.816486

Availaible at VTU HUB (Android App)

Finite Element Analysis 17ME61
2. Evaluate the integral using 2×2 point integral rule and check the solution by
exact integration. (JAN/FEB2006)

Solution: The given integral is a 2D integral. (Integral is with respect to two variables x & y). 2×2 point
integral rule for 2D integral is

Here,

By exact integration:

Availaible at VTU HUB (Android App)