UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |1

V.R.S. College of Engineering & Technology, Arasur-607107

(Accredited by NAAC and an ISO 9001:2008 Recertified Institution)

Subject Name : FINITE ELEMENT ANALYSIS

Subject code : ME6603

Year : III

Semester : VI

UNIT III
TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS
Second Order 2D Equations involving Scalar Variable Functions – Variational formulation –Finite
Element formulation – Triangular elements – Shape functions and element matrices and vectors.
Application to Field Problems - Thermal problems – Torsion of Non circular shafts –Quadrilateral
elements – Higher Order Elements.

1. Write short notes on second order 2D equations involving scalar variable functions
SECOND ORDER 2D EQUATIONS INVOLVING SCALAR VARIABLE FUNCTIONS
Two dimensional elements Two dimensional elements are defined by three or more nodes in a two
dimensional plane (i.e., x, y plane). The basic element useful for two dimensional analysis is the
triangular element.

Many engineering structures and mechanical components are subjected to loading in two directions.
Shafts, gears, couplings, mechanical joints, plates, bearings, are few examples. Analysis of many three
dimensional systems reduces to two dimensional, based on whether the loading is plane stress or plane
strain type.
Triangular elements or Quadrilateral elements are used in the analysis of such components and systems.
The various load vectors, displacement vectors, stress vectors and strain vectors used in the analysis are
as written below, the displacement vector u = [u, v]T ,
u is the displacement along x direction, v is the displacement along y direction, the body force vector f =
[ fx , fy]T

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |2

fx , is the component of body force along x direction, fy is the component of body force along y
direction.
the traction force vector T = [ Tx , Ty]T
Tx , is the component of body force along x direction, Ty is the component of body force along y
direction.
The element having mid side nodes along with corner nodes is a higher order element.
Element having a curved side is also a higher order element. A simple quadrilateral element has straight
edges and corner nodes. This is also a linear element. It can have constant thickness or variable
thickness. The quadrilateral having mid side nodes along with corner nodes is a higher order element.
Element having a curved side is also a higher order element.
The given two dimensional components is divided in to number of triangular elements or quadrilateral
elements. If the component has curved boundaries certain small region at the boundary is left uncovered
by the elements. This leads to some error in the solution.

Two dimensional stress strain equations

From theory of elasticity for a two dimensional body subjected to general
loading the equations of equilibrium are given by

The strain displacement relations are given by

The stress strain relationship for plane stress and plane strain conditions are given by the matrices shown
in the next page. are usual stress strain components, v is the poisons ratio. E is
young’s modulus. Please note the differences in [ D] matrix .
Plane Stress and Plane Strain
The 2D element is extremely important for the Plane Stress analysis and Plane Strain analysis.
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |3

Plane Stress Analysis:

It is defined to be a state of stress in which the normal stress (ċ) and shear stress (Ċ) directed
perpendicular to the plane are assumed to be zero.

Plane Strain Analysis:

It is defined to be a state of strain in which the normal to the xy plane and the shear strain are assumed to
be zero.

VARIATIONAL FORMULATION
Variational formulation refers to the construction of a functional or a variational principle that is
equivalent to the governing equations of the problem. It is nothing but the formation I which the
governing equations are translated into equivalent weighted integral statements that are not necessarily
equivalent to a variational principle.
It is common, especially in structural mechanics to express finite element formulation in vector notation
(i.e., in terms of matrices).
We know that,
B(w,u) = l(w)
Let,

 a11 a12 0  / x 
C = a 21 a 22 0  D =  / y 
 0 0 a00   1 

B and l are expressed as,

w / x   a11 a12 0 u / x 

B (w,u) = ∫ w / y  T a
 21 a 22 0  u / y  dx dy
 
 w   0 0 a00   u 

It is, B (w,u) = ∫(𝐷𝑊)T CD u dx dy

l(w) = ∫(𝑊)T f dx dy + ∫ 𝑊 T qds is the variational formulation of a two dimensional problem.

2. Explain the finite element formulation for 2D in detail

FINITE ELEMENT FORMULATION
The major steps are as follows
1. Discretization of domain into finite elements.
2. Weak formulation of governing differential equation.
3. Derivation of finite element interpolation functions.
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |4

4. Develop finite element model using weak form.

5. Assembly of finite elements to obtain the global system of algebraic equations.
6. Imposition of boundary conditions.
7. Solution of equations.
8. Post computation of solution.
Here, steps 6 and 7 remain unchanged from one-dimensional finite element analysis.
The basic element useful for two dimensional analysis is the triangular element. The simple two
dimensional elements have corner nodes.
a) Triangular element
b) Rectangular element
c) Quadratic element
d) Parallelogram element
Triangular elements: Nodal displacements of the element represent the displacements at points inside
an element. As discussed earlier, the finite element method uses the concept of shape function in
systematically developing these interpolations.
Generally triangular elements are classified into two types, they are:
(i) Constant Strain Triangular element
(ii) Linear strain triangular element

Constant Strain Triangular (CST) Element

A three noded triangular element is known as constant strain triangular (CST) element. It has six
unknown displacement degrees of freedom (u1,v1, u2,v2, u3,v3).

Displacement function for the CST element

Linear Strain Triangular (LST) Element:

Linear strain triangular element is those, which has six nodes. The degrees of freedom of an LST
element are twelve, because it has twelve unknown displacement.

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |5

The development of stiffness matrix in LST element is same as that of the CST element. The difference
is that more number of equations and it is a tedious process to solve those equations. Hence it is solved
by computer by using mathematical equations.
For plane stress applications, LST elements are preferred that the CST element. When large numbers of
nodes are used, LST element is not preferred, since cost of formation of the element stiffness equation
bandwidth is high compared to CST element.

3. Derive the section with the development of the shape function for CST element for strain
displacement matrix. (AU-NOV/DEC-2010)

Consider a typical CST element with nodes 1,2 and3 as shown in fig. let the nodal displacement,

 u1 
u 
 2
u 
displacement ,{u}   3 
 v1 
 v2 
 
 v3 

Since CST element has got two degrees of freedom at each node (u,v) the degree of freedom is 6.
Hence it has 6 generalized coordinates.

Let

u  a1  a2 x  a3 y
v  a4  a5 x  a6 y

Where a1, a2, a3, a4, a5 and a6 are globalized coordinates

u1  a1  a2 x1  a3 y1
u2  a1  a2 x2  a3 y2
u3  a1  a2 x3  a3 y3

Write above equation in matrix form

 u1  1 x 1 y1   a1 
    
u2   1 x2 y2  a2 
u  1 x y3   a3 
 3  3
1
 a1  1 x 1 y1   u1 
    
a2   1 x2 y2  u2 
 a  1 x y3  u 
 3  3  3
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |6

1 x 1 y1 
D  1 x2 y2 
1 x3 y3 
CT
D 1 
Let D

The co-factors of matrix D is

( x2 y3  x3 y2 ) ( y2  y3 ) ( x3  x2 ) ( x2 y3  x3 y2 ) ( x3 y1  x1 y3 ) ( x1 y2  x2 y1 )
C  ( x3 y1  x1 y3 ) ( y3  y1 ) ( x1  x3 ) C T  ( y2  y3 ) ( y3  y1 ) ( y1  y2 )
( x1 y2  x2 y1 ) ( y1  y2 ) ( x2  x1 ) ( x3  x2 ) ( x1  x3 ) ( x2  x1 )

We know that

1 x 1 y1 
D  1 x2 y2 
1 x3 y3 
D  1( x2 y3  x3 y2 )  x1 ( y3  y2 )  y1 ( x3  x2 )

Substitute CT and D values in equation

( x2 y3  x3 y2 ) ( x3 y1  x1 y3 ) ( x1 y2  x2 y1 ) 
D 1 
1  (y  y ) ( y3  y1 ) ( y1  y2 ) 
1( x2 y3  x3 y2 )  x1 ( y3  y2 )  y1 ( x3  x2 ) 
2 3

 ( x3  x2 ) ( x1  x3 ) ( x2  x1 ) 
substitute D-1 value in the above equation
1
 a1  1 x 1 y1   u1 
    
a2   1 x2 y2  u2 
 a  1 x y3  u 
 3  3  3
 a1  ( x2 y3  x3 y2 ) ( x3 y1  x1 y3 ) ( x1 y2  x2 y1 )  u1 
   (y  y )  
( y1  y2 )  u2 
1
a2    ( y3  y1 )
 a  1( x2 y3  x3 y2 )  x1 ( y3  y2 )  y1 ( x3  x2 )  ( x  x )
2 3

 3  3 2 ( x1  x3 ) ( x2  x1 )  u3 

The area of the triangle can be expressed as function of the x, y coordinates of the nodes 1, 2 and 3
1 x 1 y1  1
1
A  1 x2 y2  A  [1( x2 y3  x3 y2 )  x1 ( y3  y2 )  y1 ( x3  x2 )]
2
2
1 x3 y3  2 A  ( x2 y3  x3 y2 )  x1 ( y3  y2 )  y1 ( x3  x2 )

Substitute 2A value in equation

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |7

 a1  ( x2 y3  x3 y2 ) ( x3 y1  x1 y3 ) ( x1 y2  x2 y1 )  u1 
  1   
a2    ( y2  y3 ) ( y3  y1 ) ( y1  y2 )  u2 
a  2 A  ( x  x ) ( x1  x3 ) ( x2  x1 )   
 3  3 2 u3 
 a1   p1 p2 p3   u1 
  1   
 a2    q1 q2 q3  u2 
a  2 A  r u 
r3  
 3  1 r2 3

p1  ( x2 y3  x3 y2 )
p2  ( x3 y1  x1 y3 )
p3  ( x1 y2  x2 y1 )
q1  ( y2  y3 )
q2  ( y3  y1 )
q3  ( y1  y2 )
r1  ( x3  x2 )
r2  ( x1  x3 )
r3  ( x2  x1 )

From the equation we know that

u  a1  a2 x  a3 y
We can write this equation in matrix form
 a1 
 
u  1 x y  a2 
a 
 3
 a1 
 
Substitute  a2  value in the above equation
a 
 3
 p1 p2 p3   u1 
1   
u  1 x y q1 q2 q3  u2 
2A 
 r1 r2 r3  u3 
 p1 p2 p3   u1 
 
y   q1 q3  u2 
1
u 1 x q2
2A
 r1 r2 r3  u3 
 u1 
1  
u  p1  q1 x  r1 y p2  q2 x  r2 y p3  q3 x  r3 y  u2 
2A u 
 3

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |8

 u1 
 p  q x  r1 y p2  q2 x  r2 y p3  q3 x  r3 y   
u 1 1  u2 
 2A 2A 2A  
u3 
The equation in the form of
 u1 
 
u   N1 N 2 N3  u2 
u 
 3
Similarly
 v1 
 
v   N1 N2 N3  v2 
v 
 3
Shape function
p  q x  r1 y
N1  1 1
2A
p  q x  r2 y
N2  2 2
2A
p  q x  r3 y
N3  3 3
2A

Assembling the equation in matrix form

Displacement function

 u1 
v 
 1
u ( x, y )   N1 0 N2 0 N3 0  u2 
u   
 v ( x, y )   0 N1 0 N2 0 N3   v2 
u3 
 
 v3 

4. Derive the Strain – Displacement matrix [B] for CST element.

Strain – Displacement matrix [B] for CST element

Displacement function

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS Page |9

 u1 
v 
 1
u ( x, y )   N1 0 N2 0 N3 0  u2 
u   
 v ( x, y )   0 N1 0 N2 0 N3   v2 
u3 
 
 v3 
u=N1u1=N2u2+N3u3
v=N1v1=N2v2+N3v3
The strain components for CST element are;
 ex 
 
 ey 
 xy 
 

Where,
e x and e y --- Normal strain
 xy --- Shear Strain
e x  u / x
Normal Strain,
u / x   / x( N1u1  N2u2 + N3u3)
e y  v / y
Normal Strain ,
v / y   / y ( N1v1  N2v2 + N3v3)
Shear strain,  xy = u / y + v / x
 xy = u1 (N1 /y)  u 2 (N 2 /y)  u3 (N 3 /y) + v1 (N1 /x)  v2 (N 2 /x)  v3 (N 3 /x)
On rearranging

 u1 
v 
 ex 
N1 /x 0 N 2 /x 0 N 3 /x 0   1
 0 u 2 
 
 N1 /y 0 N 2 /y 0 N 3 /y  v 
ey  
 xy  N1 /y N1 /x N 2 /y N 2 /x N 3 /y N 3 /x   2
 
u 3 
 v 3 

The above equation is of the form {e} = [B] (u)

We know that,

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 10

p1  q1 x  r1 y
N1 
2A
p2  q2 x  r2 y
N2 
2A
p3  q3 x  r3 y
N3 
2A On substituting the shape functions we get,

[B] is the strain displacement matrix for the CST element. This equation is the element strain equation.
5. Derive the stress – strain relationship matrix for two dimensional elements.

Hooke’s law states that when a material is loaded within its elastic limit, the stress is directly
proportional to the strain.

Stress is directly proportional to strain

 e
 =Ee

e=
E
where
e  strain
  stress
E  young ' s mod ulus

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 11

x
e'x =
Strain E

The above fig shows the stress in y direction produces a negative strain in the x direction as a result of
poisson’s effect which is given by

v x
-e''x =
E
v
e''x =- x
E
v  poisson ' s ratio

Similarly the stress in the z direction produces a negative strain in the x direction as shown in fig

v z
-e'''x =
E
v
e'''x   z
E

By applying superposition principle to the equation we get

x y z
ex  v v
E E E
The above equation is strain equation in x direction
Similarly the strain in y and z direction are
x y z
e y  v  v
E E E
x y z
ez  v v 
E E E
Solving the above equation for normal stresses (σx, σy and σz).

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 12

E
x  ex (1  v)  ve y  vez 
(1  v)(1  2v) 
E
y  vex  (1  v)e y  vez 
(1  v)(1  2v) 
E
z  vex  vey  (1  v)ez 
(1  v)(1  2v) 

The shear stress and shear strain relationship is given by,

  G

The expressions for the three different sets of shear stresses are

E
 xy   xy
2(1  v)
E  1  2v 
 xy     xy
(1  v)(1  2v)  2 
E
 yz   yz
2(1  v)
E  1  2v 
 xy  G xy  yz     yz
(1  v)(1  2v)  2 
 yz  G yz
E
 zx  G zx  zx   zx
2(1  v)
where
E  1  2v 
G  mod ulus of rigidity =
E  zx     zx
2(1  v ) (1  v)(1  2v)  2 

Assembling the above equation

1  v v v 0 0 0 
 v 1 v 0   e 
 x   v 0 0
x
   v v 1 v 0 0 0 e 
 y   y 
 z  
0   ez 
E  0 1 2v
0 0 0
     
 xy  (1  v)(1  2v) 
2
1  2v   xy 
 yz   0 0 0 0 0   yz 
   2  
 zx   
1  2v   zx 
 0 0 0 0 0 
 2 
The above equation is in the form of {σ}=[D]{e}

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 13

6. Write down the Derivation for plane stress and plane strain in two dimensional planes.
(AU-NOV/DEC-2011)

PLANE STRESS

x y z
ex  v v
E E E
x y z
e y  v  v
E E E

For two dimensional plane stress problems, the normal stress, σz and the shear stress τxz, τyz are zero.

τxz=τyz=σz=0

The shear strain γxz, γyz are zero but ez≠0

γxz= γyz=0
substitute σz=0 in ex equation
x y
ex  v -----------1
E E
substitute σz=0 in ey equation
x y
e y  v  ------------2
E E
Solving equation 1and 2
E
x  (ex  vey )
(1  v 2 )
E
y  (vex  ey )
(1  v 2 )
We know that,
Shear stress, τxy = G γxy
G-modulus of rigidity = E/2(1+v)
γxy-shear strain
v-poisson’s ratio

E (1  v)
 xy   xy
(1  v)(1  v) 2
E E (1  v)
 xy   xy  xy   xy
2(1  v) (1  v 2 ) 2

Arranging the equations in matrix form

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 14

 
 x  1 v 0   ex 
  E   
 y   v 1 0   ey 
  1  v2
 xy 
 1 v   
0 0   xy 
 2 

The above equation is in the form of {σ}=[D]{e}

The equation gives the two dimensional stress strain relationship for plane stress problem.

Where [D]=stress-strain relationship matrix

 
1 v 0 
E  
 D 
1  v2 
v 1 0 
 1 v 
0 0 
 2 

PLANE STRAIN

For plane strain we assume the following strain to be zero ez=γxz=γyz=0

the shear stresses τxz=τyz=0, but σz≠0

from equation we know that

1  v v v 0 0 0 
 v  
 x   1 v v 0 0 0   ex 
   v v 1 v 0 0 0 e 
 y   y 
 z 
 1  2v
 E  0 0 0 0 0   ez 
     
 xy  (1  v)(1  2v) 
2
1  2v   xy 
 yz   0 0 0 0 0   yz 
   2  

 zx   zx 
  1  2v   
 0 0 0 0 0 
 2 

Simplifying the above matrix we get

 
 x  1  v v 0   ex 
  E   
 y    v 1 v 0   ey 
  (1  v)(1  2v)  1  2v   xy 
 xy   0 0 
 2 

The above equation is in the form of {σ}=[D]{e}.

The equation gives the two dimensional stress strain relationship for plane strain problem

7. Derive the stiffness matrix equation for two dimensional element (CST element)
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 15

We know that,
Stiffness matrix, k   B  DB dv
T

1 x 1 y1 
1
Where, A- area of the triangular element= A  1 x2 y2 
2
1 x3 y3 
t-thickness of the element
[B]-strain-displacement matrix
q1 0 q 2 0 q3 0 
1 
[B]= 0 r1 0 r2 0 r3 
2A 
 r1 q1 r2 q 2 r3 q3 

q1  ( y2  y3 )  0  4  4
q2  ( y3  y1 )  4  0  4
q3  ( y1  y2 )  0  0  0
r1  ( x3  x2 )  1.5  3  1.5
Where, r2  ( x1  x3 )  0  1.5  1.5
r3  ( x2  x1 )  3  0  3

[D]-stress-strain relationship matrix.

For plane stress problems,
 
1 v 0 
D  E 2 v 1 0 
1 v  1 v 
0 0 
 2 
Where, E-young’s or modulus of elasticity,
v-Poisson’s ratio.
8. Determine the stiffness matrix for the CST element shown in fig. the coordinates are given
in units of millimeters. Assume plane stress conditions. Take E=210 GPa, v = 0.25 and t =
10 mm. (April/May 2012)

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 16

Given
x1=20mm; y1=30mm
x2=80mm; y2=30mm
x3=50mm; y3=120mm
young’s modulus, E = 210 GPa = 210*109Pa = 210*103N/mm2

poisson’s ratio, v = 0.25

thickness t = 10mm

To find: Stiffness matrix [k]

Solution: we know that stiffness matrix, [k] = [B]T [D] [B] At

1 x 1 y1  1 20 30
1
y2 
1
A  1 x2 A  1 80 30
2 2
area 1 x3 y3  1 50 120

1
A 1(80*120  50*30)  20(120  30)  30(50  80)
2
A  2700mm2

We know that strain-displacement

q1  ( y2  y3 )  30  120  90
q2  ( y3  y1 )  120  30  90
 q1 0 q2 0 q3 0
1  q3  ( y1  y2 )  30  30  0
[ B]  0 r1 0 r2 0 r3 
2A  r1  ( x3  x2 )  50  80  30
 r1 q1 r2 q2 r3 q3 
r2  ( x1  x3 )  20  50  30
r3  ( x2  x1 )  80  20  60

Substitute the above values in the matrix equation

 90 0 90 0 0 0
1 
[ B]   0 30 0 30 0 60 
2A
 30 90 30 90 60 0 

 90 0 90 0 0 0
[ B] 
1  0 30 0 30 0 60 
2*2700 
 30 90 30 90 60 0 
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 17

 3 0 3 0 0 0   3 0 3 0 0 0 
30 
[ B]   0 1 0 1 0 2  [ B]  5.55*10  0 1 0 1 0 2 
3

2*2700
 1 3 1 3 2 0   1 3 1 3 2 0 

We know that, Stress – strain relationshipmatrix for plane stress is

   
1 v 0   1 0.25 0 
E   2.1*105  
 D  2 v 1 0   D  2 
0.25 1 0 
1 v  1  (0.25) 
1 v  1  0.25 
0 0   0 0 
 2   2 

4 1 0 
 D  56*10 1 4 0 
3

0 0 1.5

4 1 0   3 0 3 0 0 0 

 D[ B]  56*10 1 4 0  *5.55*10  0 1 0 1 0 2
3  3

0 0 1.5  1 3 1 3 2 0 

 4 1 0   3 0 3 0 0 0
 D[ B]  311.08 1 4 0   0 1 0 1 0 2 
0 0 1.5  1 3 1 3 2 0 
 12 1 12 1 0 2
 0 8 
 D[ B]  311.08  3 4 3 4
 1.5 4.5 1.5 4.5 3 0 

We know that

 3 0 1
 0 1 3
 3 0 3 0 0 0   
  
[ B]  5.55*10  0 1 0 1 0 2 
3 0 1
3
[ B]T  5.55*103  
 1 3 1 3 2 0   0 1 3 
0 0 2
 
0 2 0

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 18

 3 0 1
 0 1 3
   12 1 12 1 0 2 
 3 0 1 
[ B]  D [ B]  5.55*10 
T 3
 311.08  3 4 3 4 0 8 
 0 1 3   1.5 4.5 1.5 4.5 3 0 
0 0 2
 
0 2 0

 37.5 7.5 34.5 1.5 3 6 

 7.5 17.5 1.5 9.5 9 8

 34.5 1.5 37.5 7.5 3 6
[ B]  D [ B]  1.728 
T

 1.5 9.5 7.5 17.5 9 8
 3 9 3 9 6 0
 
 6 8 6 8 0 16 

Substitute [B]T[D][B] and A, t values in equation

 37.5 7.5 34.5 1.5 3 6 

 7.5 17.5 1.5 9.5 9 8

 34.5 1.5 37.5 7.5 3 6
[k ]  1.728   * 2700*10 N / mm
 1.5 9.5 7.5 17.5 9 8
 3 9 3 9 6 0
 
 6 8 6 8 0 16 

 37.5 7.5 34.5 1.5 3 6 

 7.5 17.5 1.5 9.5 9 8 

 34.5 1.5 37.5 7.5 3 6
[k ]  46.656*103   N / mm
 1.5 9.5 7.5 17.5 9 8 
 3 9 3 9 6 0
 
 6 8 6 8 0 16 

9. Evaluate the stiffness matrix for the CST element shown in fig. the coordinates are given
in units of millimeters. Assume plane stress conditions. Take E=2*10 5 N/mm2, v = 3 and t =
10 mm. (AU-Nov/Dec-2013, Jan 2006)

Given

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 19

x1=0mm; y1=0mm
x2=3mm; y2=0mm
x3=1.5mm; y3=4mm
young’s modulus, E = 2*105N/mm2, poisson’s ratio, v = 0.3, thickness t = 10mm

To find: stiffness matrix [k]

Solution: we know that stiffness matrix, [k] = [B]T [D] [B] At

1 0 0 1
1
A 1 3 0
A 1(12  0)  0  0
2
2
1 1.5 4 A  6mm2

We know that strain-displacement

q1  ( y2  y3 )  0  4  4
q2  ( y3  y1 )  4  0  4
 q1 0 q2 0 q3 0
1  q3  ( y1  y2 )  0  0  0
[ B]  0 r1 0 r2 0 r3 
2A  r1  ( x3  x2 )  1.5  3  1.5
 r1 q1 r2 q2 r3 q3 
r2  ( x1  x3 )  0  1.5  1.5
r3  ( x2  x1 )  3  0  3

Substitute the above values in the matrix equation

 4 0 4 0 0 0  4 0 4 0 0 0
1 
1.5 0 3
1 
[ B]  0 1.5 0 [ B]  0 1.5 0 1.5 0 3
2A  2*6 
 1.5 4 1.5 4 3 0  1.5 4 1.5 4 3 0 

 4 0 4 0 0 0  4 0 4 0 0 0
1 
1.5 0 3
1 
[ B]   0 1.5 0 [ B]   0 1.5 0 1.5 0 3
2*6 12
 1.5 4 1.5 4 3 0   1.5 4 1.5 4 3 0 

We know that Stress – strain relationshipmatrix for plane stress is

   
1 v 0   1 0.3 0 
E   2*105  
 D  2 v 1 0   D  2 
0.3 1 0 
1 v  1  (0.3) 
1 v  1  0.3 
0 0  0 0 
 2   2 
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 20

 1 0.3 0 
 D  219.78*10 0.3 1 0 
3

 0 0 0.3

 1 0.3 0   4 0 4 0 0 0
 
 D[ B]  219.78*10 0.3 1 0  *  0 1.5 0 1.5 0 3
3 1 
12
 0 0 0.3  1.5 4 1.5 4 3 0 

 1 0.3 0   4 0 4 0 0 0
3  
 D[ B]  18.32*10 0.3 1 0   0 1.5 0 1.5 0 3
 0 0 0.35  1.5 4 1.5 4 3 0

 4 0.45 4 0.45 0 0.9 

 D[ B]  18.32*103  1.2 1.5 1.2 1.5 0 3 
 0.525 1.4 0.525 1.4 1.05 0 

We know that

 4 0 1.5
 0 1.5 4 
 4 0 4 0 0 0  
1    
1.5 0 3
1 4 0 1.5
[ B]   0 1.5 0 [ B]T   
12 12  0 1.5 4 
 1.5 4 1.5 4 3 0 
0 0 3 
 
0 3 0 

 4 0 1.5
 0 1.5 4 
   4 0.45 4 0.45 0 0.9 
 1.5 3 
3 
1 4 0
[ B ]  D [ B ]  
T
 18.32*10  1.2 1.5 1.2 1.5 0
12  0 1.5 4 
 0.525 1.4 0.525 1.4 1.05 0 
0 0 3 
 
0 3 0 

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 21

 16.78 3.9 15.21 0.3 1.57 3.6 

 3.9 7.85 0.3 3.35 4.2 4.5 

 15.21 0.3 16.78 3.9 1.578 3.6 
[ B]T  D [ B]  1.526*103  
 0.3 3.35 3.9 7.85 4.2 4.5
 1.57 4.2 1.57 4.2 3.15 0 
 
 3.6 4.5 3.6 4.5 0 9 

Substitute [B]T[D][B] and A, t values in equation

 16.78 3.9 15.21 0.3 1.57 3.6 

 3.9 7.85 0.3 3.35 4.2 4.5

 15.21 0.3 16.78 3.9 1.578 3.6 
[k ]  1.526*103   *10*6 N / mm
 0.3 3.35 3.9 7.85 4.2 4.5
 1.57 4.2 1.57 4.2 3.15 0 
 
 3.6 4.5 3.6 4.5 0 9 

 16.78 3.9 15.21 0.3 1.57 3.6

 3.9 7.85 0.3 3.35 4.2 4.5

 15.21 0.3 16.78 3.9 1.578 3.6 
[k ]  91.6*103   N / mm
 0.3 3.35 3.9 7.85 4.2 4.5
 1.57 4.2 1.57 4.2 3.15 0 
 
 3.6 4.5 3.6 4.5 0 9 

10.For the plane strain element shown in fig. the nodal displacement are:

u1=0.005mm;v1=0.002mm
u2=0.0mm; v2=0.0mm
10.
u3=0.005mm; v3=0.0mm. determine the element stresses  x ,  y ,  xy  1 2 andp
theprincip leangle , E=70Gpa, and v=0.3 and use unit thickness for plane strain. All co-
ordinates are in mm. (AU-APR/MAY-2010)

Given: nodal displacements:

u1=0.005mm ; v1=0.002mm
u2=0.0mm ; v2=0.0mm
u3=0.005mm ; v3=0.0mm
x1=5mm ; y1=15mm.
x2=15mm ; y2=5mm.

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 22

x3=25mm ; y3=15mm.
young’s modulus, E=70Gpa=70X109 pa.=70X103 N/mm2
possion’s ratio, v=0.3
unit thickness i=1mm
to find: (1) element stresses
a. Normal stresses , σx
b. Normal stresses, σy
c. Shear stress, τxy
d. Maximum mormal stress, σ1
e. Minimum normal stress, σ2
f. Principal angle, øP
Solution:

We know that,
1 x1 y1 1 15 15
1 1
Area of the element, A= 1 x2 y2 = 1 15 5
2A 2A
1 x3 y3 1 25 15

X 1X 15 X 15  5 X 25  51X 15 X  5 X 1  151X 25  15 X 1

1
= A=100mm2
2
q1 0 q 2 0 q3 0 
1 
Strain-displacement, [B]=  0 r1 0 r2 0 r3 
2A
 r1 q1 r2 q 2 r3 q3 
Where, q1=y2-y3=5-15=-10
q2=y3-y1=15-15=0
q3=y1-y2=15-5=10
r1=x3-x2=25-15=10
r2=x1-x3=5-25=-20
r3=x2-x1=15-5=10
substitute the above values in equation

 10 0 q2 0 10 0
1 
[B]= 0 10 0  20 0 10
2A 
 10  10  20 0 10 10
Substitute area, a value
 10 0 q2 0 10 0 
1 
[B]= 0 10 0  20 0 10
2 X 100 
 10  10  20 0 10 10
We know that,
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 23

Stress-strain relationship matrix [D] for plane strain problem is,

 
1 v v 0 
D  E
 v 1 v 0 
1  v 1  2v   1  2v 
 0 0 
 2 
 
3 1  0.3 0.3 0 
D  70 X 10
 0..3 1  0.3 0 
1  0.31  2 X 0.3  1  2 X 0.3 
 0 0 
 2 
 0.7 0.3 0  3.5 1.5 0
D  134.615 X 10 0..3 0.7 0 
3
D  134.615 X 10 X 0.21.5 3.5 0
3

 0 0 0.2  0 0 1
3.5 1.5 0
D  26.923 X 10 1.5 3.5 0
3

 0 0 1
3.5 1.5 0  10 0 q2 0 10 0 
 
DB  26.923 X 10 1.5 3.5 0
3 1 
0 10 0  20 0 10
2 X 100 
 0 0 1  10  10  20 0 10 10
 35 15 0  30 35 15 

DB  134.615  15 35 0  70 15 35
 10  10  20 0 10 10 
We know that,
Stress, {σ}=[D] [B] {u}
 u1  0.005
v  0.002
 1  35 15 0  30 35 15   
u 2  
 
 134.615  15 35  70 15 35 X 
0.0
=[D] [B]   0 
v 2   10  10  20 0 10 10   0.0 
u 3  0.005
   
 v 3   0.0 
 35 X 0.005  15 X 0.002  35 X 0.005 0.003
   
 134.615  15 X 0.005  35 X 0.002  15 X 0.005  134.615 0.07 
 10 X 0.005  10 X 0.002  10 X 0.005   0.08 
   

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 24

 4.038   x   4.038 
     
  9.4223  y    9.423 
10.769    10.769
   xy   
a. Normal stresses , σx=1.038 N/mm2
b. Normal stresses, σy=9.423 N/mm2
c. Shear stress, τxy=10.769 N/mm2
We know that,
 x  y  x  y 
Maximum normal stress, σmax=σ1=      2 xy
2  2 
4.038  9.423  4.038  9.423 
    10.769
2
=
2  2 
σ1=17.83 N/mm 2

 x  y  x  y 
Minimum normal stress, σmax=σ2=      2 xy
2  2 
4.038  9.423  4.038  9.423 
    10.769
2
=
2  2 
σ2=-4.369 N/mm 2

we know that,

2 xy  2 X 10.769 
priciple angle, tan 2 øp=  tan 1  
 x  y  10.038  9.423 

2 øp=-75.96º

øp=-37.98º

Result:

a. Normal stresses , σx=1.038 N/mm2

b. Normal stresses, σy=9.423 N/mm2
c. Shear stress, τxy=10.769 N/mm2
d. Maximum mormal stress, σ1=17.83 N/mm2
e. Minimum normal stress, σ2=-4.369 N/mm2
f. Principal angle, øP=37.98º

11. The two dimensional propped beam shwon in fig. is divided in two CST element.
Determine the nodal displacement and element stresses using plate stress costions. Body
force is negleced in compression with the external forces.Take: thickness,t=10mm, young’s
moduls , E=2X105 N/mm2, poisson’s ratio,v=0.25
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 25

Given:

Thickness, t=10mm

Young’s modulus, E=2X105 N/mm2

Poisson’s ratio, v=0.25

To find: (i) nodal displacement u1, v1, u2,v2,u3,v3 and u4,v4

(ii) element stress, σ1,σ2

Solution:

Consider (1): nodal displacements u1, v1, u2,v2,u3,v3 and u4,v4

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 26

Take node 1 as origin.

For node 1 : 0,0

For node 2: 1000, 750

For node 3: 0,750

Stiffness matriux, [K]= [B]T [D] [B] A t

Where, A= area of the triangular element

1 x1 y1 1 0 0
1 1
= 1 x2 y2 = 1 1000 750
2A 2A
1 x3 y3 1 0 750

1 1000 X 750
 X 1X 1000 X 750  0  
2 2
A  375 X 10 mm
3 2

q1 0 q2 0 q3 0
1 
Strain- displacement matrix, [B] = 0 r1 0 r2 0 r3 
2A 
 r1 q1 r2 q2 r3 q3 

Where, q1=y2-y3=750-750=0
q2=y3-y1=750-0=750
q3=y1-y2=0-750=-750
r1=x3-x2=0-1000=-1000
r2=x1-x3=0-0=0
r3=x2-x1=1000-0=1000
substitute the above values in equation

 0 0 750 0 750 0 
1 
[B]= 0 1000 0 0 0 1000 
2A 
 1000 0 0 750 1000 750 
Substitute area, a value
 0 0 750 0 750 0 
[B]=
1  0 1000 0 0 0 1000 
3 
2 X 375 X 10
 1000 0 0 750 1000 750

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 27

 0 0 3 0 3 0 
[B]=
250  0 4 0 0 0 4 
3  
2 X 375 X 10
 4 0 0 3 4 3
Stress=strain relationship matrix [D] for plane stress problem is,

   
1 v 0   1 0.25 0 
E   2 X 105  
 D   D 
2 
1   0.25 
v 1 0  0.25 1 0 
1  v   1 v 
2
 1  0.25 
0 0   0 0 
 2   2 

 1 0.25 0  4 1 0 
2 X 105  2 X 105 X 0.25 
 D  0.25 1 0   D  1 4 0 
0.9375  0.9375 
 0 0 0.375 0 0 1.5

4 1 0 
 D  2 X 10 X 0.2667 1 4 0 
5

0 0 1.5

4 1 0   0 0 3 0 3 0 
 B D  2 X 10 X 0.2667 1 4 0  X  0 4 0 0 0 4 
5 250
2 X 375 X 103  
0 0 1.5  4 0 0 3 4 3

4 1 0   0 0 3 0 3 0 
250 X 2 X 105 X 0.2667 
 B D  3  1 4 0  X  0 4 0 0 0 4 
2 X 375 X 10
0 0 1.5  4 0 0 3 4 3
 0  0  0 0  4  0 12  0  0 0  0  0 12  0  0 0  4  0 
 B D  17.78  0  0  0 0  16  0 3  0  0 0  0  0 3  0  0 0  16  0 
0  0  16 0  0  0 0  0  0 0  0  4.5 0  0  6 0  0  4.5
 0 4 12 0 12 4 

 B D  17.78  0 16 3 0 3 16 
 16 0 0 4.5 6 4.5

We know that,

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 28

 0 0 4 
 0 4 0 
 0 0 3 0 3 0   
250  0 4 0 0 0 4  250  3 0 0 
3  
[B]= [B]T=
2 X 375 X 103   2 X 375 X 10  0 0 3 
 4 0 0 3 4 3
 3 0 4 
 
 0 4 3

[B]T[D]

 0 0 4   0 0 4 
 0 4 0   0 4 0 
     0 4 12 0 12 4 
  3 0 0  3 16 
250 3 0 0 250
 X   X 17.78  0 16 3 0
3
2 X 375 X 10  0 0 3  2 X 375 X 103 0 0 3  16 0 0 4.5 6 4.5
 3 0 4   3 0 4 
   
 0 4 3  0 4 3

 24 0 0 18 24 18 
 0 64 12 0 12 64 

 0 12 36 0 36 12 
[K]= 5.927 X103  
 18 0 0 13.5 18 13.5
 24 12 36 18 60 30 
 
 18 64 12 13.5 30 77.5 

 53.28 0 0 39.96 53.28 39.96 

 142.08
 0 142.08 26.64 0 26.64
 
4  
stiffnessmatri,[k ]  1X 10  0 26.64 79.92 0 79.92 26.64 
 39.96 0 0 29.97 39.96 29.97 
 
 53.28 26.64 79.92 39.96 133.2 66.6 
 39.96 142.08 26.64 29.97 66.6 172.05 
 

Consider (2): nodal displacements u1, v1, u2,v2,u3,v3 and u4,v4

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 29

Take node 1 as origin.

For node 1 : 0,0

For node 2: 1000, 750

For node 3: 0,750

Stiffness matriux, [K2]= [B]T [D] [B] A t

Where, A= area of the triangular element

1 x1 y1 1 0 0
1 1
= 1 x2 y2 1 1000 0
2A 2A
1 x3 y3 1 1000 750
1 1000 X 750
 X 1X 1000 X 750  0  
2 2
A  375 X 10 mm
3 2

q1 0 q2 0 q3 0
1 
Strain- displacement matrix, [B] = 0 r1 0 r2 0 r3 
2A 
 r1 q1 r2 q2 r3 q3 

Where, q1=y2-y3=0-750=-750
q2=y3-y1=750-0=750
q3=y1-y2=0-0=0
r1=x3-x2=1000-1000=0
r2=x1-x3=0-1000=-1000
r3=x2-x1=1000-0=1000
substitute the above values in equation

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 30

 750 0 750 0 0 0 
1 
[B]=
 0 0 0 1000 0 1000
2A
 0 750 1000 750 1000 0 
Substitute area, a value
750 0 750 0 0 0
1
[B]= 0 0 0 1000 0 1000
2 X 375 X 103
0 750 1000 750 1000 0
 3 0 3 0 0 0 
[B]=
250  0 0 0 4 0 4 
3  
2 X 375 X 10
 0 3 4 3 4 0 
Stress=strain relationship matrix [D] for plane stress problem is,

   
1 v 0   1 0.25 0 
E   2 X 105  
 D   D 
2 
1   0.25 
v 1 0  0.25 1 0 
1  v   1 v 
2
 1  0.25 
0 0   0 0 
 2   2 

 1 0.25 0  4 1 0 
2 X 105  2 X 105 X 0.25 
 D   0.25 1 0   D   1 4 0 
0.9375 0.9375
 0 0 0.375 0 0 1.5

4 1 0 
 D  2 X 10 X 0.2667 1 4 0 
5

0 0 1.5

4 1 0   3 0 3 0 0 0 
 B D  2 X 10 X 0.2667 1 4 0  X  0 0 0 4 0 4 
5 250
2 X 375 X 103  
0 0 1.5  0 3 4 3 4 0 

 12 0 12 4 0 4 

 B D  17.78  3 0 3 16 0 16
 0 4.5 6 4.5 6 0 

We know that,

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 31

 3 0 0 
 0 0 3
 3 0 3 0 0 0   
250  0 0 0 4 0 4  250  3 0 4 
3  
[B]= [B]T=
2 X 375 X 103   2 X 375 X 10  0 4 3 
 0 3 4 3 4 0 
0 0 4
 
0 4 0

 3 0 0 
 0 0 3
   12 0 12 4 0 4 
    3 16 0 16 
250 3 0 4
[B]T[D] [B]=   X 17.78  3 0
2 X 375 X 103  0 4 3 
 0 4.5 6 4.5 6 0 
0 0 4
 
0 4 0

 36 0 36 12 0 12 
 0 13.5 18 13.5 18 0 

 36 18 60 30 24 12 
[K]= 5.927 X10 3  
 12 13.5 30 77.5 18 64
 0 18 24 18 24 0 
 
 12 0 12 64 0 64 

 79.92 0 79.92 26.64 0 26.64 

 0 
 0 29.97 39.96 29.97 39.96
 
4  
stiffnessmatri,[k ]  1X 10  79.92 39.96 133.2 66.6 53.28 26.64 
 26.64 29.97 66.6 172.05 39.96 142.08
 
 0 39.96 53.28 39.96 53.28 0 
 26.64 26.64 142.08 142.08 
 0 0

Global stiffness matrix [K]

Assemble the stiffness matrix equation

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 32

Global stiffness matrix, [k]=1x104

U1 V1 U2 V2 U3 V3 U4 V4

53.28+ 0+0 -79.92 26.64 0+0 -39.96 -53.28 39.96 U1

79.92 -26.64

0+0 +142.08 39.96 -29.97 -26.64+ 0+0 26.64 -142.28 V1

+29.97 -39.96

-79.92 39.96 133.2 -66.6 -53.26 26.64 0 0 U2

26.64 -29.97 -66.6 172.06 39.96 -142.08 0 0 V2

0+0 -26.64+ -53.28 39.96 79.92 0+0 -79.92 26.64 U3

-39.96 +53.28

- 0+0 26.64 -142.28 0+0 29.97+1 39.96 -29.96 V3

39.96+ 42.08

-26.64

-53.28 26.64 0 0 -79.92 39.96 133.2 -66.6 U4

39.96 -142.08 0 0 26.64 -29.97 -66.6 172.05 V4

Global stiffness matrix [k]=1X104

 u1 v1 u2 v2 u3 v3 u4 v4 
 133.2 0 79.92 26.64 0 66.6  53.28 39.96  u1
 v
 0 172.05 39.96 29.97 66.6 0 26.64 142.08
1

  u2
 79.922 39.96 133.2 66.6 53.28 26.64 0 0 
v2
 26.64 29.97 66.6 172.06 39.96 142.08 0 0 
  u3
 0 66.6 53.28 39.96 133.2 0 79.92 26.64 
 66.6 v3
0 26.64 142.08 0 172.08 39.96 29.97 
  u4
 53.28 26.64 0 0 79.92 39.96 133.2 66.6 
 39.96 142.08 v4
 0 0 26.64 29.97 66.6 172.05 

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 33

we know that, general force equation is {F}=[k] {u}

 
 F1 X   133.2  u1 
F   0 79.92 26.64 0 66.6 53.28 39.96   
 1 y  v1
 0 172.05 39.96 29.97 66.6 0 26.64 142.08  
 F2 x    u2 
   79.922 39.96 133.2 66.6 53.28 26.64 0 0  
 F2 y   v2 
  1X 10  26.64 29.97 66.6 172.06 39.96 142.08 0  
4
 0
 F3 x    u3
 0 66.6 53.28 39.96 133.2 0 79.92 26.64   
 F3 y   v3 
   66.6 0 26.64 142.08 0 172.08 39.96 29.97   
 F4 x    u4 
 F4 y    53.28 26.64 0 0 79.92 39.96 133.2 66.6  
   39.96 142.08  v4 
 0 0 26.64 29.97 66.6 172.05 

Applying boundary conditions

1. Node 1, and node 4 are fixed. So, u1, v1, and u4,v4 are =0

2. Node 2, is moving in x direction, so, u2=0 but, v2=0

3. At node 3, a point load of 75000N is acting in x direction, so, F3-x=75,000N

4. Body force is neglected. So, the remaining forces are zero.

F1 X
F1 y
F2  x
F2  y
F3 x
5. F3 y =0
F4  x
F4  y

Substitute the above values in equation

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 34

 
 0   133.2 0 
 0   0 79.92 26.64 0 66.6 53.28 39.96 0 
   0 172.05 39.96 29.97 66.6 0 26.64 142.08  
 0    u 
   79.922 39.96 133.2 66.6 53.28 26.64 0 0  2
 0  0
 3
 1X 104  26.64 29.97 66.6 172.06 39.96 142.08 0 0  
75 X 10    u
 0 66.6 53.28 39.96 133.2 0 79.92 26.64   3 
 0  v 
   66.6 0 26.64 142.08 0 172.08 39.96 29.97   3 
 0   0 
 0   53.28 26.64 0 0 79.92 39.96 133.2 66.6   
   39.96 142.08 0
 0 0 26.64 29.97 66.6 172.05 

In the above equation u1, v1, u2,v2,u3,v3 and u4,v4=0 so, delete the corresponding row and column of [K]

matrix. Hence the equation reduces to

 0   133.2 53.28 26.64  u2 
 3 4   
75 X 10   1X 10  53.28 133.2 0  u3 
 0  172.05  v3 
   26.64 0

 0   133.2 53.28 26.64  u2   0   133.2 53.28 26.64  u2 

         
7.5   53.28 133.2 0  u3  18.75   53.28 133.2 0  u3 
 0   26.64 172.05  v3   0   26.64 172.05  v3 
   0    0

 0  133.2 53.28 26.64  u2 

    
18.75   0 279.72 26.64  u3 
18.75  0 4349.81  v3 
   0

-4349.81 v3=18.75

V3=-0.00431mm

279.72u3+26.64v3=18.75

279.72u3+26.64X(-0.00431)=18.75

u3=0.067mm

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 35

u2=0.02766mm

nodal displacements:

u1=0mm; v1=0mm

u2=0.02766mm; v2 =0mm

u3 =0.067mm ,v3 =-0.00431

u4 =0mm; v4=0mm

stress in each element:

we kow that. Stress, {σ}=[D] [B] {u}

for element (1) nodal displacement u1, v1, u3,v3 and u4,v4

 u1 
v 
 1
 
 0 4 12 0 12 4     x  14.295 
 0 16 3 0      
  u   y    3.574  N / mm
2
 3 16
 6 0 0 4.5 6 4.5  3    0.345
 xy   
stress , {σ}=17.78  v3 
 
u4 
v 
 4

For element (2)

 u1 
v 
 1
u 
 12 0 12 4 0 4   2   x   5.595 
 3   v2     
  y   0.2492  N / mm
2
 0 3 16 0 16  u 
 0 4.5 6 4.5 6 0   3     4.196 
stress , {σ}=17.78   v3   xy   
 
 
 
 
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 36

Result:

1. Nodal displacement: u1=0, v1=0, u2 =0.02766,v2 =0,u3 =0.067,v3 =-0.00431 and u4 =0,v4 =0

2. Element stresses:
 x  14.295 
   
 y    3.574  N / mm
2

  0.345
For elements (1)  xy   
 x   5.595 
   
 y   0.2492  N / mm
2

   4.196 
For elements (2)  xy   

12. A thin plate is subjected to surface traction as shown in fig. calculate the global stiffness
matrix. (Nov/Dec 2011, April/May 2010)

Take: t=25mm, young’s moduls , E=2X105 N/mm2, poisson’s ratio,v=0.30.

assume plane stress condition.

Given :

Thickness, t:25mm

Young’s modulus, E=2X105 N/mm2

Poisson’ ratio, v=0.30

Breadth=250mm

Length, l=500mm

Tensile surface traction, T=0.4 N/mm2

The tensile surface traction is converted into nodl force.

F=1/2 TA=1/2 X T X (bxt)

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 37

=1/2X0.4X250X25

Nodal force, F=1250N

To find:

Stiffness matrix [k]

Solution:

Consider (1): nodal displacements u1, v1, u2,v2,u3,v3 and u4,v4

Take node 1 as origin.

For node 1 : (0,0)

For node 2: (500, 250)

For node 3: (0,250)

Stiffness matriux, [K]= [B]T [D] [B] A t

Where, A= area of the triangular element

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 38

1 x1 y1 1 0 0 1
1 1  X 1X  500 X 250  0 
= 1 x2 y2 1 500 250 2
2A 2A
1 x3 y3 1 0 250 A  62.5 X 103 mm2

q1 0 q2 0 q3 0
1 
Strain- displacement matrix, [B] = 0 r1 0 r2 0 r3 
2A 
 r1 q1 r2 q2 r3 q3 

Where, q1=y2-y3=250-250=0
q2=y3-y1=250-0=250
q3=y1-y2=0-250=-250
r1=x3-x2=0-500=-500
r2=x1-x3=0-0=0
r3=x2-x1=500-0=500
Substitute the above values in equation

 0 0 500 0 250 0 
1 
[B]=
 0 500 0 0 0 500 
2A
 500 0 0 250 500 250 
Substitute area, a value
 0 0 500 0 250 0 
[B]=
1  0 500 0 0 0 500 
3 
2 X 62.5 X 10
 500 0 0 250 500 250
 0 0 1 0 1 0 
[B]=
250  0 2 0 0 0 2 
3  
2 X 62.5 X 10
 2 0 0 1 2 1
Stress=strain relationship matrix [D] for plane stress problem is,

   
1 v 0   1 0.3 0 
E   2 X 105  
 D   D 
2 
1   0.3 
v 1 0  0.3 1 0 
1  v   1 v 
2
 1  0.3 
0 0  0 0 
 2   2 

 1 0.3 0 
2 X 105 
 D  0.3 1 0 
0.91 
 0 0 0.35

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 39

 1 0.3 0   0 0 1 0 1 0 
2 X 105   250  0 2 0 0 0 2 
 B D   0.3 1 0 X 3  
0.91 2 X 62.5 X 10
 0 0 0.35  2 0 0 1 2 1

 0 0.6 1 0 1 0.6 
 2 
 B D  439.56  0 2 0 0 0
 0.7 0 0 0.35 0.7 0.35

We know that,

 0 0 2 
 0 2 0 
 0 0 1 0 1 0   
250  0 2 0 0 0 2  250  1 0 0 
[B]= 3   [B]T=  
2 X 62.5 X 10 2 X 62.5 X 103  0 0 1 
 2 0 0 1 2 1
 1 0 2 
 
 0 2 1
[B]T[D][B]=
 0 0 2 
 0 2 0   
   1 0.3 0   0 0 1 0 1 0 
250 1 0 0 2 X 105   250  0 2 0 0 0 2 
 
  0.3 1
X 0 X  
2 X 62.5 X 103  0 0 1  1   0.32   2 X 62.5 X 103
1  0.3  2 0 0 1 2 1
 1 0 2  0 0 
   2 
 0 2 1

 1.4 0 0 0.7 1.4 0.7 

 0 4 0.6 0 0.6 4 

 0 0.6 1 0 1 0.6 
[K]= 1373.59 X 103  
 0.7 0 0 0.35 0.7 0.35
 1.4 0.6 1 0.7 2.4 1.3 
 
 0.7 4 0.6 0.35 1.3 4.35 

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 40

 1923.08 0 0 961.513 1923.026 961.513 

 5494.36 
 0 5494.36 824.154 0 824.154
 
 
stiffnessmatri,[k ]  1X 103  0 824.24 1373.59 0 1373.59 824.24 
 961.513 0 0 480.7565 961.513 480.7565 
 
 1923.026 84.24 1373.59 961.513 3296.616 1785.6697 
 961.513 5494.36 824.154 480.756 1785.67 5975.1165 
 

Consider (2): nodal displacements u1, v1, u2,v2,u3,v3 and u4,v4

Take node 1 as origin.

For node 1 : (0,0)

For node 2: (500, 0)

For node 3: (500,250)

Stiffness matriux, [K2]= [B]T [D] [B] A t

Where, A= area of the triangular element

1 x1 y1 1 0 0 1 500 X 250
1 1  X 1X  500 X 250  0  
= 1 x2 y2 1 500 0 2 2
2A 2A
1 x3 y3 1 500 250 A  62.55 X 10 mm
3 2

q1 0 q2 0 q3 0
1 
Strain- displacement matrix, [B] = 0 r1 0 r2 0 r3 
2A 
 r1 q1 r2 q2 r3 q3 

Where, q1=y2-y3=0-250=-250 q2=y3-y1=250-0=250 q3=y1-y2=0-0=0

r1=x3-x2=500-500=0 r2=x1-x3=0-500=-500 r3=x2-x1=500-0=500
Substitute the above values in equation
 250 0 250 0 0 0 
1 
[B]=  0 0 0 500 0 500 
2A
 0 250 500 250 1000 0 
 1 0 1 0 0 0 
Substitute area, a value [B]=
1  0 0 0 2 0 2 
2 X 62.5 X 103  
 0 1 2 1 2 0 
Stress=strain relationship matrix [D] for plane stress problem is,

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 41

   
1 v 0   1 0.3 0   1 0.3 0 
E   2 X 105   2 X 105 
 D   D   D  0 
2 
1   0.3 
v 1 0  0.3 1 0  0.3 1
1  v   1 v 
2
 1  0.3 
0.91 
 0 0 0.35
0 0  0 0 
 2   2 

 1 0.3 0   1 0 1 0 0 0 
2 X 105 
0  X  0 0 0 2 0 2
250
 B D   0.3 1  
0.91 2 X 62.5 X 103
 0 
0 0.35  0 1 2 1 2 0 

 1 0 1 0.6 0 0.6
 2 
 B D  439.56 3 0 0.3 2 0
 0 0.35 0.7 0.35 0.7 0 

 1 0 0 
 0 0 1
 1 0 1 0 0 0   
250  0 0 0 2 0 2   1 0 2 
[B]= 3 
[B]T= 2 X 103  

2 X 62.5 X 10
 0 1 2 1 2 0   0 2 1 
0 0 2
We know that,  
0 2 0

 1 0 0 
 0 0 1
   1 0 1 0 0 0 
 1 0 2   
[B]T[D] [B]= 2 X 103 X 439.96   X  0 0 0 2 0 2 
 0 2 1   0 1 2 1 2 0 
0 0 2  
 
0 2 0

 1 0 1 0.6 0 0.6 
 0 0.35 0.7 0.35 0.7 0 

 1 0.7 2.4 1.3 1.4 0.6 
[K]= 0.8791  
 0.6 0.35 1.3 4.35 0.7 4 
 0 0.7 1.4 0.7 1.4 0 
 
 0.6 0 0.6 4 0 4 

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 42

 1373.59 0 1373.59 824.154 0 824.154 

 0 
 0 480.7565 961.513 480.7565 961.13 
 
 
stiffnessmatri,[k2 ]   1373.59 961.513 3296.616 1785.667 1923.026 824.154 
 824.154 480.7565 1785.667 5975.165 961.513 5494.36 
 
 0 961.513 1923.026 961.513 1923.026 0 
 824.154 5494.36 5494.36 
 0 824.154 0

Global stiffness matrix [K] Assemble the stiffness matrix equation

Global stiffness matrix, [k]=1x103

U1 V1 U2 V2 U3 V3 U4 V4

1923.026+ 0+0 -1373.59 824.154 0+0 - -1923.026 961.513 U1

1373.59 961.513+

-824.154

0+0 5494.36+ 961.513 - - 0+0 824.154 -5494.36 V1

480.7565 480.7565 824.154+-

961.513

-1373.59 961.513 3296.61 - -1923.026 824.154 0 0 U2

6 1785.667

824.24 - - 5975.116 961.513 -5494.36 0 0 V2

480.7565 1785.66 5

0+0 - 0+- 0+961.51 1373.59+1 0+0 -1373.59 824.154 U3

824.154+ 1923.02 3 923.026

-961.513 6

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 43

-961.513+- 0+0 0+824.1 0+- 0+0 961.513 3296.616 -1785.667 V3

824.154 54 5494.36

-1923.026 824.154 0 0 -1373.59 961.513 3296.616 -1785.667 U4

961.513 -5494.36 0 0 824.154 - -1785.667 -5975.1165 V4

480.7565

Global stiffness matrix [k]=1X103

 u1 v1 u2 v2 u3 v3 u4 v4 
 3296.616 0 1373.59 824.154 0 1785.667 1923.026 961.513 
u1
 v1
 0 5975.1165 961.513 480.7565 1785.667 0 824.154 5496.36 
  u2
 1373.59 961.513 3296.616 1785.667 1923.026 824.154 0 0  v2
 824.154 480.7565 1785.667 5975.1165 961.513 5494.36 0 0 
  u3
 0 1785.667 1923.026 961.513 3296.616 0 1373.59 824.154 
 1785.667 v3
0 824.154 5494.36 0 5975.165 961.513 480.7565 
  u4
 1923.026 824.24 0 0 1373.59 961.513 3296.616 1785.667 
 961.513 v4
 5494.36 0 0 824.154 480.7565 1785.667 5975.1165 

 k   1X 103 X
 
 3296.616 0 1373.59 824.154 0 1785.667 1923.026 961.513 

 0 5975.1165 961.513 480.7565 1785.667 0 824.154 5496.36 
 
 1373.59 961.513 3296.616 1785.667 1923.026 824.154 0 0 
 824.154 480.7565 1785.667 5975.1165 961.513 5494.36 0 0 
 
 0 1785.667 1923.026 961.513 3296.616 0 1373.59 824.154 
 1785.667 0 824.154 5494.36 0 5975.165 961.513 480.7565 
 
 1923.026 824.24 0 0 1373.59 961.513 3296.616 1785.667 
 961.513 5494.36 480.7565 1785.667 5975.1165 
 0 0 824.154

THERMAL PROBLEM
Temperature effect:

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 44

Distribution of the change in temperature (ΔT) is known as strain. Due to the change in temperature can
be considered as an initial strain e0.

T 
{ e0} = T 
 
 0  for general and plain stress problem

Where, ΔT --- change in temperature

α --- co-efficient of thermal expansion

For Plain strain problem

T 
{ e0} = (1+µ) T 
 
 0 

Where,e0 ----initial strain,

µ ----- Poissons ratio

ΔT --- change in temperature

α --- co-efficient of thermal expansion

the stress and strain are releated by the following relation,
e = D (e-e0)
σ = D (Bu - e0)
where, e0--- initial strain
σ---- stress
D----stress strain relationship matrix
B----strain displacement relationship matrix
u---displacement
Heat transfer in two dimensional elements: the shape function and formulae are similar to that of the
CST element.

Governing Equation for 2D Heat transfer by conduction and convection

  2T  2T 
k  2  2   h(T  T )  0
 x y 
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 45

Weak form of the equation

T w T w
 x x
dxdy  
y y
dxdy   hTw( x, y )dxdy

  hT w( x, y )dxdy

 N N j N N j 
K ij conv  k  i dxdy   i dxdy
 x x y y 

1
N i (x,y) = (αi + βi x + γi y)
2 Ae

 12   12 1  2   1 2 1  3   1 3 
k  
K  conv
e
 1  2   1 2  22   22  2  3   2 3 
4A 
 1  3   1 3  2  3   2 3  32   32 

 q1  0
  hlT  
q 2   1
q  2  
 3 1

hpl 2 1
k conv 
6 1 2
 p 1

13. Problem: determine the temperature distribution on the plate

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

 UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 46

  12   12 1  2   1 2 1  3   1 3 
k  
[k ]e   2  1   1 2  22   22  2  3   2 3 
4A 
  1  3   1 3  2  3   2 3  32   32 


Element 1 and 2 Element 3

β1 = -0.15 , 1 = 0, β1 = 0.15 , 1 = -0.4

β2 = 0.15 , 2 = -0.4 β2 = 0.3 , 2 = 0

β3 = 0 , 3 = 0.4 β3 = -0.15 , 3 = 0.4

 0.0225  0.0225 0 
 [k ] 
1.5  
 0.0225 0.1825  0.16
2 1
[k ] cond cond
4
 0.4  0.15  0  0.16 0.16 
2

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

 UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 47

 0.028125  0.028125 0 

 10  0.028125 0.228125  0.2
 0  0.2 0.2 

 0.28125  0.28125 0 
  0.28125 2.28125  2
 0 2 2 

 0.1825  0.045  0.1825

[k ] 3

1.5   0.045 0.09  0.045 
cond
1 
4  2   0.4  0.15  0.1825  0.045 0.1825 
2
1.14  0.28125  0.86 

 0.5625  0.28125
 1.14 

k conv  hpl 
2 1
6 1 2
 p 1

Element i j k

1 1 2 3

2 5 4 3

3 1 3 5

0 0 0 0 0 0 
 0 2 1
hl 
[k ]  [k ]
2
conv
1
conv

 0 2.5 1.25
6
0 1 2 0 1.25 2.5 
0  0 
hlT    
Q 1  93.75
2    
1 93.75
 [k ]condn  [k ]conv
e
[k ] Thermal

 0.28125  0.28125 0 
1 
[k ]  [k ]   0.28125
2
4.78  0.75
th th

 0  0.75 4.5 
 1.14  0.28125  0.86 

[k ]th   0.28125
3 ME6603 FEA 0.5625  0.28125 Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
  0.86  0.28125 1.14 
 UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 48

T1   0 
T   
 2   93.75 
[k ]th T3   [Q]G  [Q ] G
  93.75  93.75 
T   93.75 
 4  
T5   0 

 1.42125  0.28125  0.28125 0  0.86 

 0.28125 4.78  0.75 0 0 
 
[k ]G   0.28125  0.75 9.5625  0.75  0.28125
 
 0 0  0.75 4.78  0.28125
  0.86 0  0.28125  0.28125 1.42125 
 1.42125  0.28125  0.28125 0  0.86  T1  0
 0.28125  0.75    1
 4.78 0 0  T2   
 0.28125  0.75 9.5625  0.75  0.28125 T3   93.752
 
 0 0  0.75 4.78  0.28125 T4  1
   
  0.86 0  0.28125  0.28125 1.42125  T5  0
Substitute forT4 & T5 as 180º and evaluate T1, T2 ,T3

STIFFNESS MATRIX FOR BI LINEAR RECTANGULAR ELEMENT

 x  y   x  y 
N1  1  1   N3    
 2a  2b   2a  2b 

 x  y   x  y 
N 2   1   N 4  1   
 2a  2b   2a  2b 
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 49

 2(a 2  b 2 ) a 2  2b 2 ( a 2  b 2 ) (b 2  2a 2 ) 
 2 (a 2  b 2 ) 
k  a  2b 2(a 2  b 2 ) (b 2  2a 2 )
2


6ab  (a 2  b 2 ) (b 2  2a 2 ) 2(a 2  b 2 ) (2b 2  a 2 ) 
 2 
 (b  2a ) ( a 2  b 2 ) a 2  2b 2 2(a 2  b 2 ) 
2

14. Explain about the torsion of non-circular shaft

TORSION OF NON-CIRCULAR SHAFT
The first application area is the torsion of Non-Circular sections. The governing differential equation is
1  2 1  2
+ + 2 = 0
G x 2 G y 2
where G - shear modulus of the material

 - is the angle of twist.

The governing equation for the torsion problem is given by

1  2 1  2
+ + 2 = 0
G x 2 G y 2

 2  2
+ = - 2Gθ
x 2 y 2

 
τ zx = τ zy = -
y x
On the free boundary  = 0.

Here  - is a stress function

The shear stresses within the shaft are related to the derivatives of  with respect to x and y.

On the free boundary  = 0. This is the case of a Poisson’s Equation

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 50

To derive the weak form multiply the governing equation with a weighting function w(x,y)
 2  2
+  2Gθ  0
x 2 y 2

 2  2
 ( x 2
+
y 2
 2Gθ ) w( x, y )dxdy  0

On simplifying, the shape function can be expressed as

1
N i (x,y) = (αi + βi x + γi y)
2 Ae

Finite element equation for torsional triangular element

 f1  1
  A 
 f 2   2Gθ 1
f  3
 3 1

Where, θ angle o twist, f force vector

QUADRILATERAL ELEMENTS:
Each node is allowed to move in only one direction ‘+x’ in one dimensional problems, but in two
dimensional problems, each node is permitted to move in the two directions (i.e) x and y. Each node has
two degrees of freedom (Nodal Displacements)

Derivation of shape function for the quadrilateral elements:

Consider a four node rectangular element, the coordinates of the nodes are
Node 1: x1,y1 Node 2: x2,y2 Node 3: x3,y3 Node 4: x4,y4,
It has eight unknown displacement degrees of freedom.
The shape functions are
N1 =(1-x/l-y/h+xy/lh), N2 =(x/l-xy/lh) N3 =(xy/lh) N4 =(y/h-xy/lh)
Assembling the equations in matrix form
Displacement form,

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 51

u1 
v 
 1
u 2 
 
u ( x, y )  N 1 0 N2 0 N3 0 N4 0  v 2 
u   
 v ( x, y )   0 N1 0 N2 0 N3 0 N 4  u 3 
v3 
 
u 4 
v 
 4

15. Derive shape functions for 4 noded rectangular parent element by using natural co-

ordinate system and coordinate transformation.

Consider a four noded rectangular element as shown in fig. the parent element is defined in ε and η

co-ordinates natural co-ordinates ε is varying from -1 to 1 and η is also varying -1 to 1

We know that,

Shape function value is unity at its own node and its value is zero at other nodes.

At node 1:

(co-ordinates ε=-1, η=-1)

Shape function N1=1 at node 1

N1=0 at node 2, 3, and 4

N1 has to be in the form oof N1=C(1-ε) (1-η)

Where C is cnsatant

Substitute ε=-1 and η=-1 in equation N1=C(1+1) (1+1)

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 52

N1=4C 1=4C C=1/4

1
N1  1   1   
4

At node 2:

(co-ordinates ε=1, η=-1)

Shape function N2=1 at node 2

N2=0 at node 1, 3, and 4

N2 has to be in the form oof N2=C(1+ε) (1-η)

Where C is cnsatant

Substitute ε=1 and η=-1 in equation N2=C(1+1) (1+1)

N2=4C 1=4C C=1/4

1
N 2  1   1   
4

At node 3:

(co-ordinates ε=1, η=1)

Shape function N3=1 at node 3 N3=0 at node 1, 2, and 4

N3 has to be in the form of N3=C(1+ε) (1+η)

Where C is cnsatant

Substitute ε=1 and η=1 in equation N3=C(1+1) (1+1)

N3=4C 1=4C C=1/4

1
N3  1   1   
4

At node 4:

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 53

(co-ordinates ε=-1, η=1)

Shape function N4=1 at node4 N4=0 at node 1, 2, and 3

N4 has to be in the form of N4=C(1-ε) (1+η)

Where C is cnsatant

Substitute ε=1 and η=1 in equation N4=C(1+1) (1+1)

N3=4C 1=4C C=1/4

1
N 4  1   1   
4
u 
 
v 
Consider a point P with coordinate (ε,η) if displacement function u= represents to be

displacement components of a point located at (ε,η) then,

u=N1u1+N2 u2=N3 u3+N4 u4 and v=N1v1+N2 v2=N3v+N4 v4

it can be in matrix form as,

 u1 
v 
 1
u2 
 
u   N 0 N2 0 N3 0   v2 
u  1  
v   0 N1 0 N2 0 N3  u3 
v3 
 
u4 
v 
 4

In the isoparametric formulation. For global system, the co0ordinates of the nodal points are (x 1,y1),

(x2,y2) , ), (x3,y3) and ), (x4,y4) in order to get mapping the co-ordinate point p is defined as

x=N1x1+N2 x2=N3 x3+N4 x4 and y=N1y1+N2 y2=N3y+N4 y4

The above equations can be written in matrix form as

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 54

 x1 
y 
 1
 x2 
 
x N 0 N2 0 N3 0   y2 
u  1   
 y  0 N1 0 N2 0 N 3   x3 
 y3 
 
 x4 
y 
 4
HIGHER ORDER ELEMENTS:

The higher order element can be either a complex or a multiple element.

Higher order elements are nothing but if the interpolation polynomial is the order of two or more
elements.
In the higher order elements, some secondary node are introduced in addition to the primary nodes in
order to match the number of nodal degrees of freedom with the number of constants in the interpolation
polynomial.

16. Derive the shape function for eight noded triangular element.
Consider a eight noded rectangular element. It belongs to the serendipity family of elements. It
consists of eight nodes, which are located on the boundary.

1 1
N1   (1   )(1   )(1     ) N 2   (1   )(1   )(1     )
4 4

1 1
N3   (1   )(1   )(1     ) N 4   (1   )(1   )(1     )
4 4

1 1
N5  (1   )(1   )(1   ) N6  (1   )(1   )(1   )
2 2

1 1
N7  (1   )(1   )(1   ) N8  (1   )(1   )(1   )
2 2

17. Derive the shape function for six noded triangular element.

Consider a six noded triangular element. It belongs to the serendipity family of elements. It consists of
six nodes, which are located on the boundary.

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 55

N1  L1  2L1 1 N2  L2  2L2  1 N3  L3  2L3 1

N 4  4 L1L2 N5  4 L2 L3 N6  4 L1L3

Two Marks Question and Answers.

UNIT-III
1. What meant by plane stress analysis?

Plane stress is defined to be a state of stress in which the normal stress and shear stress directed
perpendicular to the plane are assumed to be zero.

2. Define plane strain analysis. (Nov 2011)

Plane strain is defined to be state of strain normal to the xy plane and the shear strains are assumed to be zero.

3. Write down the stiffness matrix equation for two dimensional CST elements.


Stiffness matrix     t

T
  -Strain displacement    -Stress strain matrix  -Strain displacement matrix

4. How do you define two dimensional elements?( AU-NOV/DEC-2010)

Two dimensional elements are define by three or more nodes in a two dimensional plane. The
basic element useful for two dimensional analysis is the triangular elements. The elements are
analysed in two different axis as x and y axis. the two dimensional elements are used to analyse
bar element, beam elemnt and the truss element. The structed analysis are carried out by means
of governing equation and the displacement function with its boundary conditions.
5. What is CST element?
Three noded triangular element is known constant strain triangle (CST) which is shown in fig. it
has six unknown displacement degree of freedom(u1 v1,u2 v2,u3 v3).the element is called CST
because it has a constant strain throughout it.
Merits:
Calculation of stiffness matrix is easier.
Demerits:
The strain variation within the element is considered as constant. So, the results will be
poor.
6. What is LST element?( AU-NOV/DEC-2011)
Linear strain triangular element:
Six noded triangular element is known as linear strain triangular (LST). it has twelve unknown
displacement degree of freedom. The displacement function for the element are quadratic instead

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 56

of linear as in the CST.the nodes are arise in between the nodes in the CST(constant strain
triangular elemnt. the elements are analysed and are a time consuming one.
7. What is QST element?
Ten noded triangular element is known as quadratic strain triangle (QST) it is also called as
cubic displacement triangle.they are used to analyse the large structure to find the closest value
and to rectify the errors occur due to analysis. increasing the nodes will always gives us the
closest or nearest value to the original solution and we may get small differences like tolerance
value while we Are diving it as many element.
8. Write a displacement function equation for CST element.
 u1 
v 
 1
u  x , y  
 N1 0 N 2 0 N 3 0 u 2  
Displacement function u=    
vx, y   0 N 1 0 N 2 0 N 3   v 2 
u 3 
 
 v3 
Where N1, N2, N3 are shape functions.
9. Write a stain-displacement matrix for CST element. (AU-NOV/DEC-2013)
Strain –displacement matrix for CST element is,
q1 0 q 2 0 q 3 0 
B   0 r1 0 r2 0 r3 
1 
2A
 r1 q1 r2 q 2 r3 q 3 
Where, A= area of the element
q1=Y2-Y3; q2=y3-y1; q3=y1=y2
r1=x3-x2; r2=x1-x3; r3=x2-x1.
10. Write down the stiffness equation for two dimensional CST element.
Stiffness matrix,
k   BT DAt
Where, [B]-strain –displacement matrix
[D]-Stress-strain matrix
A-area of the element
t-thickness of the element.
11. Write down the expression for the shape functions for a constant strain triangular element.
(AU-APR/MAY-2010)
For CST element,
p q x  r1 y
Shape function, N1= 1 1
2A
p q x  r2 y
N2= 2 2
2A
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 57

p 3 q3 x  r3 y
N3=
2A
Where, p1=x2y3-x3y2
P2=x3y1-x1y3
P3=x1y2-x2y1
q1=Y2-Y3; q2=y3-y1; q3=y1=y2
r1=x3-x2; r2=x1-x3; r3=x2-x1.
12. Differentiate CST and LST element.
Constant strain triangular element (CST):
CST element is a two dimensional linear element (simplex element). It has only three
primary nodes at the corners. The strain is constant throughout the element. The polynomial
function for CST element is u=a1+a2x+a3y; v=a4+a5x+a6y.
Linear strain triangular element (LST):
LST is a two dimensional non-linear element (complex element). It has three primary
nodes at the corners and three secondary nodes at the midsides. The strain is varying linearly
inside the element. The polynomial function for the LST element are
u=a1+a2x+a3y+a4x2+a5y2+a6xy; v=a7+a8x+a9y+a10x2+a11y2+a12xy.
13. What is meant by two dimensional scalar variable problem?
Two dimensional scalar variable problem:
If the geometry and material properties of any element are described by two spatial
coordinates then that element is referred as two dimensional finite element and in a problem
containing that element if the measured parameter is having only one quantity (magnitude only)
and not having direction of application then it is referred to as two – dimensional scalar variable
problem.
14. Specify the applications of two –dimensional problems.
Applications of two dimensional problems:
1. The plates under bi-axial loading, to find the load, stress, strain and displacement of
the plates.
2. The bending of plates to find the load and moment acting on the beam and the
displacement function.
3. The temperature distribution on the surface due to heat transfer to find the
temperature distribution in the element.
15. Write short note on finite element modeling of two dimensional element.
Finite element modeling of two dimensional element:
Finite element modeling is the discretization of bigger element of irregular shape into
many number of calculatable regular shapes of small sized elements. This process is followed in
one dimensional object or element to reform the tapered rod into cylindrical rod. In two
dimensional element they are discretised as triangular element to get the regular interpretation of
the two dimensional element.
16. Differentiate simplex and complex elements?

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 58

Simplex element:
The non-structured problems are elements are discretized in different parameters. If it is
discretisized as a triangular element then the element is known to be simplex element.
Complex element:
The complicated structure which are in need or the easy to solve or find the solution some
of the element are discretisized as rectangular, quadrilateral, parallelogram such kind of
elements are known to be complex element.
17. What are structural and non-structural problems?

Structural problems: In structural problems, displacement at each nodal point is obtained. By

using these displacement solutions, stress and strain in each element can be calculated.
Non-structural problems: In non-structural problems, temperatures or fluid pressure at each nodal
point is obtained. By using these values, properties such as heat flow, fluid flow, etc., for each
element can be calculated.

18. Define shape function. (AU MAY 2008)

In finite element method, field variables with in an element are generally expressed by the
following approximate relation:
Ф(x, y) = N1(x, y) ф1 + N2(x, y) ф2 + N3(x,y) ф3
N1 N2&N3 are also called shape functions because they are used to express the geometry or shape of
the element.
ф1 ф2 & ф3 are the values of the field variable at the nodes and N1 N2&N3 are the interpolation
functions
19. If a displacement field in x direction is given by u: 2 x2 +4 y2 + 6xy. Determine the Strain
in x direction.
U: 2x2+4y2+6xy
Strain, e ∂u/∂x = 4x+6y
20. What are the ways in which a three dimensional problem can be reduced to a two
dimensional approach.
Plane stress: one dimensional is too small when compared to other two dimensions.
Example: Gear thickness is small

Plane strain: one dimensional is too large when compared to other two dimensions.
Example: long pipe [length is long compared to diameter] Axisymmetric: geometry
is symmetry about the axis.
Example: Cooling tower.
21. What is the purpose of isoparametric elements? (AU DEC 2007)
It is difficult to represent the curved boundaries by straight edge finite elements. A large
number of finite elements may be used to obtain reasonable resemblance between original body
and the assemblage. In order to overcome this drawback, isoparametric elements are used. i.e., for
problem involving curved boundaries, a family to elements is known as “isoparametric elements are
used.

22. Give examples for essential (forced or geometric) and non-essential boundary

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 59

conditions. (AU DEC 2010)

The geometric boundary conditions are displacement, slope, etc. the natural boundary
conditions are bending moment, shear force, etc.
23. What are the types of non-linearity? (AU MAY 2010)
i. Non – linearity in material behavior from point to point.
ii. Non – linearity in loading- deformation relation.
iii.Geometric Non – linearity
iv.Change in boundary condition for different loading

24. Define frequency of vibration.

It is the number of cycles described in one second. Unit is Hz.

25. Define Damping ratio.
It is defined as the ratio of actual damping coefficient (C) to the critical damping
coefficient (Cc).
Damping ration ε = C / Cc = C / 2mωn
26. What is meant by longitudinal vibrations?
When the particles of the shaft or disc moves parallel to the axis of the shaft, then the
vibrations are known as longitudinal vibrations.

27. What is meant by transverse vibrations?

When the particles of the shaft or disc move approximately perpendicular to the axis
of the shaft, then the vibrations are known as transverse vibrations.

28. Define magnification factor.

The ratio of the maximum displacement of the forced vibration (xmax) to eh static deflection under the
static force (Xo) is known as magnification factor.
29. Write down the expression of longitudinal vibration of bar element.
Free vibration equation for axial vibration of bar element is,
[K] {u} = ω2[m]{u}
Where, u – displacement
[K] – stiffness matrix
ω – Natural frequency
[m] – Mass Matrix
30. Write down the expression of governing equation for free axial vibration of rod.
The governing equation for free axial vibration of a rod is given by,
Where, E – young’s modulus,
A – Cross-sectional area
ρ - Density
31. Write down the expression of governing equation for transverse vibration of beam.
The governing equation for free transverse vibration of a beam

Where, E – young’s modulus

I – moment of inertia
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 60

A – Cross-sectional area
ρ – Density
32. Write down the expression of transverse vibration of beam element.
Free vibration equation for transverse vibration of beam element is,
[K] {u} = ω2[m]{u}
Where, [K] = stiffness matrix for beam element
33. What are the types of Eigen value problems?
There are essentially three groups of method of solution,
1. Determinant based methods.
2. Transformation based methods.
3. Vector iteration methods.
34. State the principle of superposition.
It states that for linear systems, the individual responses to several disturbances or
driving functions can be superposed on each other to obtain the total response of the system.
35. Define resonance.
When the frequency of external force is equal to the natural frequency of a vibrating
body, the amplitude of vibration becomes excessively large. This phenomenon is known as
resonance.
36. Define Dynamic Analysis.
When the inertia effect due to the mass of the components is also considered in
addition to the externally applied load, then the analysis is called dynamic analysis.
37. What are the methods used for solving transient vibration problems?
There are two methods for solving transient vibration problems. They are:
a) Mode superposition method
b) Direct integration method
38. State the two difference between direct and iterative methods for solving system of
equations.
Direct Method
i) It gives exact value.
ii) Simple, take less time.
iii) Determine all the roots at the same time.
Iterative Method
i) It gives only approximate solution.
ii) Time consuming and labourious.
iii) Determine only one root at the time.
39. Define linear dependence and independence of vectors.
Linear dependence : The vectors X1, X2,…Xn are said to be linearly dependent if there exist scalars
λ1, λ2, . . . λn (not all zero) such that,
λ1X1 + λ2X2 + . . . + λnXn = 0
Independence: The vectors X1, X2,…Xn are said to be linearly independent if,
λ1X1 + λ2X2 + . . . + λnXn is equal to zero such that
λ1= 0, λ2= . . . . . . = λ2

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 61

40. Define Heat transfer.

Heat transfer can be defined as the transmission of energy from one region to another region
due to temperature difference.
41. Write down the stiffness matrix equation for one dimensional heat conduction element.
Stiffness matrix, [K] =
Where,
A = area of the element, m2
K = thermal conductivity, W/mK
l = length of the element, m
42. Write down the expression of shape function, N and temperature function, T for one
dimensional heat conduction element.
For one dimensional heat conduction element,
Temperature function, T = N1T1 + N2T2

43. Write down the finite element equation for one dimensional heat conduction with free end
convection.
Finite element equation for one dimensional element heat conduction with free end
convection is given by, T = N1T1 + N2T2+ N3T3
44. Define path line.
A path line is defined as locus of points through which a fluid particle of fixed identity passes
as it moves in space.
45. Define streamline.
A streamline is an imaginary line that connects a series of points in space at a given instant in
such a manner that all particles falling on the line at that instant have velocities whose vectors are
tangent to the line.
46. Define Inviscid flow.
A inviscid flow is a frictionless flow characterized by zero viscosity. A viscous flow is one in
which the fluid is assumed to have non-zero viscosity.

47. Define two dimensional element.

Two dimensional elements
Two dimensional elements are defined by three or more nodes in a two dimensional plane
(i.e., x, y plane). The basic element useful for two dimensional analysis is the triangular
element.

48. Define Plane Stress and Plane Strain

The 2d element is extremely important for the Plane Stress analysis and Plane
Strain analysis.
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 62

Plane Stress Analysis:

It is defined to be a state of stress in which the normal stress () and shear stress
() directed perpendicular to the plane are assumed to be zero.
Plane Strain Analysis:
It is defined to be a state of strain in which the normal to the xy plane and the shear
strain are assumed to be zero
49. Write the expression for Temperature Effects
Distribution of the change in temperature (ΔT) is known as strain. Due to the
change in temperature can be considered as an initial strain e0.
σ = D (Bu - e0)
50. Define variational formulation
Variational formulation refers to the construction of a functional or a variational principle that is
equivalent to the governing equations of the problem. It is nothing but the formation I which the
governing equations are translated into equivalent weighted integral statements that are not
necessarily equivalent to a variational principle.
51. Write the Stress-strain relationship matrix [D] for plane strain problem.

Stress-strain relationship matrix [D] for plane strain problem is,

 
1 v v 0 
D  E
 v 1 v 0 
1  v 1  2v   1  2v 
 0 0 
 2 
52. Write the Stress-strain relationship matrix [D] for plane stress problem.

 
1 v 0 
E  
 D  2 
v 1 0 
1  v   1 v 
0 0 
 2 

53. Write the Strain displacement relationship matrix [B] for CST element.
[B]-strain-displacement matrix
q1 0 q2 0 q3 0
1 
[B]= 0 r1 0 r2 0 r3 
2A 
 r1 q1 r2 q2 r3 q3 

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET

UNIT-III / TWO DIMENSIONAL SCALAR VARIABLE PROBLEMS P a g e | 63

ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET