UNIT-IV / TWO DIMENSIONAL VECTOR 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 IV
TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS
Equations of elasticity – Plane stress, plane strain and axisymmetric problems – Body forces and
temperature effects – Stress calculations - Plate and shell elements.

1. Explain the equation of elasticity in detail.

Elasticity
Elasticity is the property of a deformable body due to which the body recovers its original shape upon
the removal of forces causing deformation.
Assumptions made on elasticity:
 Perfectly elastic
 Homogeneous
 Isotropic
Elasticity equations are used for solving structural mechanical problems. There are four basics sets of
elasticity equation. They are,
i. Equilibrium equation
ii. Compatibility equation
iii. Strain-displacement relationship equation
iv. Stress-Strain relationship equation
Equilibrium equation

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |2

Compatibility equation

Strain-displacement relationship equation

Stress-Strain relationship equation

2. Differentiate plane stress and plane strain.

PLANE STRESS: - A 3D problem can be reduced to a plane stress condition if it is characterized by
very small dimensions in one of the normal directions. A thin plate with a cut out subjected to in-plane
loading. Thin plate subjected to in-plane loading

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |3

In these cases the stress components z, xz, & yz are zero and it is assumed that no stress component
varies across the thickness. The state of stress is then specified by x, y and xy only, (functions of x &
y) and is called plane stress. The stress strain relations are given by

PLANE STRAIN:-
There exist problems involving very long bodies i.e. a body whose geometry and loading do not vary
significantly in the longitudinal direction. Such problems are referred to as plane strain problems.
Some typical examples include a long cylinder such as a tunnel, culvert or buried pipe, a laterally loaded
retaining wall, a long earth dam, and a loaded semi-infinite half space such as a strip footing on a soil
mass. In all these problems, the dependant variable can be assumed to be functions of only x & y co-
ordinates provided that we consider a cross-section some distance away from the two ends.
If we further assume that ‘w’ the displacement component in the ‘z’ direction is zero at every cross-
section, then the non-zero strain components will be
x = u/x ; y = v/y ; xy = u/y + v/x
and the strain components z, xz, yz will vanish. The dependant stress variables are x, y & xy and the
constitutive relation for an elastic isotropic material is given by

It is important to note here that only εz = 0 but z  0.

εz = z/E - /E ; x  /E; y  

z  - x  y

Axisymmetric Elements Most of the three dimensional problems are symmetry about an axis of
rotation. Those types of problems are solved by a special two dimensional element called as
axisymmetric element.

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |4

Many engineering problems involve solids of revolution (axisymmetric solids) subject to axially
symmetric loading.
Examples are a circular cylinder loaded by uniform internal or external pressure or other axially
symmetric loading as shown in

Because of symmetry the stress components are independent of the angular co-ordinate ‘’ and hence all
the derivatives with respect to ‘ ‘ vanish and the components x , r are zero. The strain displacement
relation are given by
εr = u/ x ; ε = u/r; εz = w/z ; rz = u/z + w/r

Axisymmetric Formulation
The displacement vector u is given by

The stress σ is given by

The strain e is given by

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |5

Temperature Effects
The thermal force vector is given by

3. Problem :- Assuming plane stress conditions evaluate the stiffness matrix for the element shown
in Fig. Assume E= 2 x 105 N/cm2 and =0.3. u1=0.000, v1=0.0025, u2=0.0012, v2=0.000, u3
=0.0000 & v3= 0.0025.

1 = y2 – y3 = 0 – 1 = -1 1 = -(x2 – x3) = 0 – 2 = -2

2 = y3 – y1 =1+1= 2 2 = -(x3 – x1) = 0 – 0 = 0

3 = y1 – y2 = -1– 0 = -1 3 = -(x1 – x2) = 2 – 0 = 2

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |6

4. Derivation of shape function for axisymmetric element. (triangular element)

Consider an axisymmetric triangular element with nodes 1, 2 and 3 shown in fig.
Let the nodal displacement be u1 w1, u2 w2,and u3 w3.
 u1 
w 
 1
 u 
Displacement, {u}=  2 
w2 
u3 
 
 w3 
Since the triangular element has two degree of freedom at each node, it has 6 generalized co-
ordinates.
Displacement functions, u=a1+a2r+a3z

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |7

W=a4+a5r+a6z
Where, a1, a2, a3, a4, a5, a6 are global or generalized co-ordinates.
Let u1= a1+a2r1+a3z1
u2= a1+a2r2+a3z2
u3= a1+a2r3+a3z3
write the above equations in matrix form,
u1  1 r1 z1   a1 
    
u 2   1 r2 z 2  a 2 
u  1 r z  a 
 3  3 3  3 
1
 a1  1 r1 z1   u1 
    
a 2   1 r2 z 2  u 2 
a  1 r z 3  u 
 3  3  3
1 r1 z1 
Let D= 1 r2 z 2 
1 r3 z3 

-1 CT
D =
D

Find the co-factors of matrix D.

r2 z2 1 z2
C11=+ =(r2z3-r3z2) C12 =- =-(z3-z2)=(z2-z3)
r3 z3 1 z3

1 r2 r1 z1
C13=+ =+(r3-r2) C21= =-(r1z3-r3z1)=r3z1-r1z3
1 r3 r3 z3

1 z1 1 r1
C22=+ =(z3-z1) C23=- =-(r3-r1)=(r1-r3)
1 z3 1 r2

r1 z1 1 z1
C31=+ =r1z2-r2z1 C32=- =-(z2-z1)=(z1-z2)
r2 z2 1 z2

1 r1
C33=+ =(r2-r1)
1 r2

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |8

r2 z 3  r3 z 2  z 2  z 3  r3  z 2 
C= r3 z1  r1 z 3  z 3  z1  r1  r3 
r1 z 2  r2 z1  z1  z 2  r2  r1 
r2 z 3  r3 z 2 z 2  z3 r3  z 2
C = r3 z1  r1 z 3
T
z 3  z1 r1  r3
r1 z 2  r2 z1 z1  z 2 r2  r1

We know that,

1 r1 z1  1 r1 z1 
D= 1 r2 z 2  D = 1 r2 z 2 
1 r3 z3  1 r3 z3 

D =1 r2 z 3  r3 z 2  -r1 z 3  z 2  +z1 r3  r2 

Substitute CT and D values in equation

r2 z 3  r3 z 2 z 2  z3 r3  z 2
1
D = -1
x r3 z1  r1 z 3 z 3  z1 r1  r3
r2 z 3  r3 z 2  - r1 z 3  z 2  + z1 r3  r2 
r1 z 2  r2 z1 z1  z 2 r2  r1

1
 a1  1 r1 z1   u1 
    
a 2   1 r2 z 2  u 2 
a  1 r z 3  u 
 3  3  3

r2 z 3  r3 z 2 z 2  z3 r3  z 2  u1 
1  
= x r3 z1  r1 z 3 z 3  z1 r1  r3 x u 2 
r2 z 3  r3 z 2  - r1 z 3  z 2  + z1 r3  r2 
r1 z 2  r2 z1 z1  z 2 r2  r1 u3 

Then area of the triangle can be expressed as a function of the r,z co-ordinates of the nodes 1,2 and 3.

1 r1 z1 
1 
A= 1 r2 z 2 
2
1 r3 z3 

[1 r2 z 3  r3 z 2  -r1 z 3  z 2  +z1 r3  r2  ]

1
A=
2

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS Page |9

2A= r2 z 3  r3 z 2  -r1 z 3  z 2  +z1 r3  r2 

Substituting equation,

 a1  r2 z 3  r3 z 2 z 2  z3 r3  z 2  u1 
  1  
a 2   r3 z1  r1 z 3 z 3  z1 r1  r3 x u 2 
a  2 A r z  r z z1  z 2 r2  r1 u3 
 3 1 2 2 1

 a1   1  2  3  u1 
  1   
a 2   1  2  3  x u2 
a  2 A    3  u3 
 3  1 2

Where, 1 = r2 z3  r3 z2 ;  2 = r3 z1  r1 z 3 ;  3 = r1 z 2  r2 z1

1  z 2  z1 ;  2  z 3  z1 ;  3  z1  z 2 ;

 1  r3  r2 ;  2  r1  r3 ;  3  r2  r1 ;

u=a1+a2r+a3z

we can write this equation in matrix form,

 a1 
 
u= 1 r z a 2 
a 
 3

 1  2  3   u1 
1   
= 1 r z 1  2  3  x u 2 
2A 
 1  2  3  u 3 

 u1 
1  
= [ 1   1 r   1 z  2   2 r   2 z  3   3 r   3 z ] x u 2 
2A u 
 3

 u1 
 1  1 r   1 z  2   2 r   2 z  3   3 r   3 z  
U= x u 2 
2A 2A 2A u 
 3

The above equation is in the form of

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 10

 u1   w1 
   
U= N1 N2 N 3  u 2  W= N1 N2 N 3   w2 
u  w 
 3  3

 1  1 r   1 z 2  2r   2 z 3  3r   3 z
Where, shape function,N1= N2= N3 =
2A 2A 2A

We can write equation as follow

U=N1u1+N2u2+N3u3

W=N1w1+N2w2+N3w3

Assembling the equation in matrix form,

 u1 
w 
 1
 u r , z   N 1 0 N2 0 N3 0  u 2 

u(r,z)=     
wr , z   0 N1 0 N2 0 N 3   w2 
u3 
 

 w3 


5. Derive the strain-displacement matrix [B] for axisymmentric triangular element

Displacement function for axisymmetric triangular element is given by,

 u1 
w 
 1
 u r , z   N 1 0 N2 0 N3 0  u 2 

Displacement function, u(r,z)=     
wr , z   0 N1 0 N2 0 N 3   w2 
u3 
 

 w3 


We can write, u =N1u1+N2u2+N3u3

W=N1w1+N2w2+N3w3

The strain components are,

u N1 N 2 N 3
Radial strain, er=  u1  u2  u3
r r r r

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 11

N1 N 2 N 3
er= u1  u2  u3
r r r

N1 N N
circumferential strain, eϴ= u1  2 u 2  3 u 3
r r r

w
longitudinal strain, ez=
z

N1 N 2 N 3
ez= w1  w2  w3
r r r

u w
shear strain,  rz  
z r

N1 N 2 N 3 N1 N 2 N 3
 rz  u1  u2  u3  w1  w2  w3
r r r r r r

Arranging equation in matrix form

 N 1 N 2 N 3 
 r 0 0 0   u1 
r r  
 er   N N2 N3   w1 
e   1 0 0 0  
   r r r u 2 
  N 1 N 2 N 3  w 
 ez   0 0 0  2
 rz   z z z   u 
N N 1 N 2 N 2 N 3 N 3   3 
 1 w 
 z r z r z r   3 

From equation, we know that

 1  1 r   1 z 2  2r   2 z 3  3r   3 z
N1= N2= N3=
2A 2A 2A

N 1  1 N 2  2 N 3  3
Partial differentiation-   
r 2A r 2A r 2A

N1 1  1  z N2 1 2  z
   1  1     2  2 
r 2A  r r  r 2A  r r 

N3 1 3  z
   3  3 
r 2A  r r 

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 12

N 1  1 N 2  2 N 3  3
  
r 2A r 2A r 2A

Substitute the values.

 
   u1 
 1 0 2 0 3 0   w1 
  
3  3z
B  1  1   1   1 z 0
2
 2 
 2z
0  3 
u 
0  2 
 w2 
2A r r r r r r
 0 1 0 2 0  3  
  u3
 1 1 2 2 3  3  
 w3 

The above equation n the form of

{e}=[B] {u}

Where, [B]=strain displacement metrix,

 
 1 0 2 0 3 0
 
1  1  z 2  2z 3  3z
B   1  1 0  2  0  3  0
2A  r r r r r r 
 0 1 0 2 0 3 
 
 1 1 2 2 3 3 

 1  r2 z 3  r3 z 2  2  r3 z1  r1 z 3  3  r1 z 2  r2 z1

1  z 2  z 3  2  z 3  z1  3  z1  z 2

 1  r3  r2  2  r1  r3  3  r2  r1

6. Derive the stress - strain matrix [D] for axisymmentric triangular element.

By using Hooke’s law, we derived the following stresses equations,

x 
E
1  v 1  2v  x

e 1  v   vey  vez 
y 
E
1  v 1  2v  x

ve 1  v   vey  vez 
z 
E
1  v 1  2v  x

ve 1  v   vey  vez 

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 13

E  1  2v 
 xz  X  X
1  v 1  2v   2  xz
Substitute x=r and y=ϴ in the above equations,
 
E
e 1  v   ve  vez 
Radial stress, r
1  v 1  2v  r
 
E
ve 1  v   ve  vez 
Circumferential stress, 
1  v 1  2v  r
Longitudinal stress,  z  1  v 1  2v  ver  ve  1  v vez 
E

E  1  2v 
  X  X
Shear stress, xz 1  v 1  2v   2  xz
Arranging the above equations, in matrix form
 r  1  v v v 0  e 
   v 1  v v 0  r 
  E    e 
     e 
 z  1  v 1  2v  
v v 1 v 0
1  2v   z 
 rz    
2   rz 
0 0 0

The above equation is in the form of,

{σ}=[D] {e}
Where, [D]=stress –strain relationship matrix,
1  v v v 0 
 v 1 v v 0 
E 
= 1  v 1  2v   v v 1 v 0 
 1  2v 
 0 0 0
2 

7. Derive element stiffness matrix for axi-symmetric triangular element.

Element stiffness matrix [K] for Axi-symmetric triangular element We know that,

Stiffness matrix, [k] =  [ B]T [ D][ B]dv

v

Where, v---volume=2ΠrA Coordinates,r = (r1+r2+r3)/3

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 14

 
 1 0 2 0 3 0
 
1  1  z 2  2z 3  3z
B   1  1 0  2  0  3  0
2A  r r r r r r 
 0 1 0 2 0 3 
 
 1 1 2 2 3 3 

 1  r2 z 3  r3 z 2  2  r3 z1  r1 z 3  3  r1 z 2  r2 z1

1  z 2  z 3  2  z 3  z1  3  z1  z 2

 1  r3  r2  2  r1  r3  3  r2  r1

BODY FORCE
A body force is distributed force acting on every elemental volume of the body
Unit: Force per unit volume.
Example: Self weight due to gravity
 u1 
 w1 
u 2 
 
2  u T f r dA   w2 f
 
u3 
 w3 

TEMPERATURE EFFECTS:

When the free expansion is prevented in the axisymmetric element, the change in temperature causes
stresses in the element.

Let ΔT be the raise in temperature and α be the co-efficient of thermal expansion. The thermal force
vector due to raise in temperature is given by

{F}t = [B]T [D] {e}t X 2ΠrA

For axi-symmetric triangular element

 F1u 
 F1w   
   T 
F t   F2u  and strain et  T 
F

 2w  
0

 3u 
F T 
F 
 3w 

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 15

STRESS CALCULATION

Elemental Stresses In Any Axi-Symmetric Problems Are As Follows

 Radial strain
 Circumferential strain
 Longitudinal strain
 Shear strain

u N1 N 2 N 3
Radial strain, er=  u1  u2  u3
r r r r

N1 N 2 N 3
er= u1  u2  u3
r r r

N1 N N
circumferential strain, eϴ= u1  2 u 2  3 u 3
r r r

w N1 N 2 N 3
longitudinal strain, ez= ez= w1  w2  w3
z r r r

u w
shear strain,  rz  
z r

8. Problem: The nodal co-ordinates for an axisymmetric triangular element are given below:

r1 = 10 mm z1 = 10 mm r2 = 30 mm z2 = 10 mm

r3 = 30 mm z3 = 40 mm Evaluate [B] matrix for the element

Given : Co-ordinates: r1 = 10 mm, z1 = 10 mm; r2 = 30 mm, z2 = 10 mm, r3 = 30 mm, z3 = 40 mm

To find: Displacement matrix [B].

Because of symmetry the stress components are independent of the angular co-ordinate ‘’ and hence all
the derivatives with respect to ‘ ‘ vanish and the components x , r are zero. The strain displacement
relation are given by

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 16

SOLUTION

Strain Displacement matrix

1
 [( r2 z 3  r3 z 2 )  r1 ( z 3  z 2 )  z1 (r3  r2 )]
2
1
 [(30  40)  (30  10)  10(40  10)  (30  30)]
2
1
  (600)
2
 A  300 mm 2

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 17

r1  r2  r3 10  30  30
Co-ordinates: r 
3 3
r  23.3334
z  z 2  z 3 10  10  40
z 1 
3 3
z  20

1  1z
900 0  20
  1  
 (30)   8.571
r r 23.334 23.334
  z  100 (20  20)
 2  2  2   (30)   8.571
r r 23.334 23.334
  z  200 20  20
 3  3  3   (0)   8.571
r r 23.334 23.334

Substituting the values in strain displacement matrix

[D]

9. Differentiate plate and shell elements used in FEA

PLATE AND SHELL ELEMENTS

Plate is a flat surface having considerably large dimensions as compared to its thickness. Slabs in civil

engineering structures, bearing plates under columns, many parts of mechanical components are the
common examples of plates. In this chapter, we are considering bending of such plates under lateral
loads. The bending properties of a plate depend greatly on its thickness. Hence in classical theory we
have the following groups
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 18

(i) thin plates with small deflections

(ii) thin plates with large deflections, and

(iii) thick plates

In thin plates with small deflections theory, the following assumption are made

(a) There is no deformation in the middle plane of the plate. This plane remains neutral during bending.

(b) Points of the plate lying initially on a normal to the middle surface of the plate remain on the normal
to the same surface even after bending.

(c) The normal stresses in the direction transverse to the plate are negligible.

This theory is satisfactory for plates with ratio of thickness to span exceeding 1/10 and the ratio of
maximum deflection to thickness less then 1/5

. Many engineering problems lie in the above category and satisfactory results are obtained by classical
theories of thin plates.

Stresses in the middle plane are negligible, if the deflections are small in comparison with thickness. If
the deflections are large, the in plane stresses developed in the so called neutral surfaces are to be
considered. This gives rise to theory of thin plates with large deflections, in which geometric non-
linearity is incorporated.

Displacement models for plate analysis

Category I: C2-Continuity element i.e. second order continuity elements in which second derivates of
‘w’ are also nodal unknowns.
Category II: C1-Continuity elements i.e. first order continuity elements in which highest order of
derivatives of ‘w’ is one only.
Category III: C0-Continuity element i.e. the elements in which only continuity of nodal variables are to
be ensured.
Analysis of shell element:
A shell is a curved surface. Due to their shape they transfer most of the load applied on their surface as
in plane forces (membrane forces) rather than by flexure. Hence the shells are examples of strength
through form rather than mass.
Civil engineers use them as roofs to get large column free areas covered. Cylindrical shells, domes
hyperbolic parabolic shells etc. are common examples of shell roofs.
Cooling towers, conical shells are also commonly used shells.
Mechanical and chemical engineers use shells as pressure vessels and as components of many
machines.
Shells may be classified as singly curved or doubly curved. Classification of shell surfaces is attempted
on the basis of Gauss curvature (product of principle curvature in two perpendicular directions).

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 19

If the Gauss Curvature is positive, zero, negative the surface will be classified as synclastic,
developable, anticlastic respectively. Further classification is possible depending upon whether a shell is
translational surface, a ruled surface or a surface of revolution.

Forces on Shell element

A typical shell element and various stress resultants acting on it. It may be noted that the sign
convention is:

(i) Coordinate direction are as per right hand thumb rule

(ii) A force acting on +ve face in +ve direction or –ve face –ve direction is +ve
(iii) A +ve force acting on +ve z-direction produces +ve moment, about mid surface

The four different approaches used to generate the shell are

1. Flat Elements
2. Curved Elements
3. Solid Elements
4. Degenerated Solid Elements.

Flat quadrilateral shell elements have use limitations even in linear analysis, since a mesh that consists
of strictly flat elements may be impossibly to construct over a doubly-curved shell surface. For large
deflection nonlinear analysis this deficiency becomes more pronounced. Even if the initial mesh satisfies
the flat element restriction, the deformations can become so large that warping of the elements can be
significant. Finding ways of handling warped element geometries is thus of fundamental importance for
quadrilateral shell elements.

The current approach to deriving the quadrilateral plate bending element utilizes reference lines.
Hrennikoff [00] first used this concept for plate modeling where the goal was to come up with a beam
framework useful as a model for bending of flat plates.

Park and Stanley [ 00, 00] used the reference line concept in their development of several plate and shell
elements based on the ANS formulation. The reference lines were used to find beam-like curvatures;
these curvatures were then used to find the plate curvatures through various Assumed Natural Strain
distributions. These plate and shell elements were of Mindlin-Reissner type, and the reference lines were
treated as Timoshenko beams. The present element is a Kirchhoff type plate and the reference lines are
thus treated like Euler-Bernoulli (or Hermitian) beams.

10. A hollow cylinder of inside diameter 100mm and outside diameter 140mm is subjected to
an internal pressure of 4 N/mm2. As shown in fig. by using two elements on the 15mm
length shown in fig. calculate the displacements at the inner radius.
Take E=2x105 N/mm2 and v=0.3.

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 20

(16)
Given:

Inner diameter, de=100mm

Inner radius, re=50mm
Outer diameter, De=140mm
Outer radius, Re=70mm
Internal pressure, p=4 N/mm2
Young’s modulus, E=2x105 N/mm2
Poisson’s ratio=0.3
To find: nodal displacement: u1,w1,u2,w2,u3,w3,u4,w4.
Solution:

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 21

For element (1) nodal displacements, u1,w1,u2,w2,u3,w3,u4,w4.

Co-ordinates:
At node 1:
r1=50mm
z1=15mm
At node 2:
r2=50mm
z2=0mm
At node 3:
r3=70mm
z3=15mm
r  r  r3 50  50  70
we know that, r= 1 2 =
3 3
r=56.66667mm
z  z 2  z 3 15  0  15
z= 1 =
3 3
z=10mm
area of the triangle element=1/2xbreath x height=1/2x(70-50)x15
A=150mm
We know that matrix for axisymmetric triangular element (1)
[K]1=2πrA [B]T [D] [B]
1  v v v 0 
 v 1 v v 0 
E 
Stress-strain relationship matrix, [D]=
1  v 1  2v   v v 1 v 0 
1  2v 
 0 0 0
2 

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 22

1  0.3 0.3 0. 0 
 
2 x10 5  0.3 1  0.3 0.3 0 
[D]=  0.3 0.3 1  0.3 0 
0.5
 1  2 x0.3 
 0 0 0
2 

1  0.3 0.3 0. 0 
 
2 x10 5  0.3 1  0.3 0.3 0 
  0.3 0.3 1  0 .3 0 
0.5
 1  2 x0.3 
 0 0 0
2 

We know that,

Strain-displacement matrix,
 
 1 0 2 0 3 0
 
1  1  1z 2  2z 3  3z
B   1  0  2  0  3  0
2A  r r r r r r 
 0 1 0 2 0 3 
 
 1 1 2 2 3 3 
 1  r2 z 3  r3 z 2  50 X 15  70 X 0
 1  750mm 2
 2  r3 z1  r1 z 3  70 X 1550 X 15
 2  300mm 2
 3  r1 z 2  r2 z1  50 X 0  50 X 15
 3  750mm 2
 1  z 2  z 3  0  15
 1  15mm
 2  z 3  z1  15  15
 2  0mm
 3  z1  z 2  15  0
 3  15mm
 1  r3  r2  70  50
 1  20mm

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 23

 2  r1  r3  50  70
 2  20mm
 3  r2  r1  50  50
 3  0mm
1 1z 750 20 X 10
 1    (15)   1.7647mm
r r 56.667 56.6667
2  z 300  20 X 10
 2  2  0  1.7647mm
r r 56.667 56.6667
3  z  750
 3  3   15  0  1.7647mm
r r 56.667
Substitute these values.

 
  15 0 0 0 15 0
B  1 1.7647 0 1.7647 0 1.7647 0 
2 X 150 0 20 0  20 0 0
 
 20  15  20 0 0 15
 
  15 0 0 0 15 0
B  3.333 X 10 3 1.7647 0 1.7647 0 1.7647 0 
0 20 0  20 0 0
 
 20  15  20 0 0 15
0.7 0.3 0. 0  
0.3 0.7 0.3 0    15 0 0 0 15 0
 
BD  384.6153 X 10 3 0.3 0.3 0.7 0  X 3.333 X 10 3 1.7647 0 1.7647 0 1.7647 0 
 0  0 20 0  20 0 0
 
 0 0 0     
 0.2  20 15 20 0 0 15
 
 9.9706 6 .5294  6 11.0294 0
BD  1.282 X 10 3  3.2647 6 1.2353  6 5.7353 0
 3.9706 14 0.5294  14 5.0294 0
 
4 3 4 0 0 3

We know that,

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 24

 
  15 0 0 0 15 0
B  3.333 X 10 3 1.7647 0 1.7647 0 1.7647 0 
0 20 0  20 0 0
 
 20  15  20 0 0 15
 15 1.7647 0 20 
 0 0 20  15 

 0 1.7467 0  20
BT  3.333 X 10 3  
 0 0  20 0 
 15 1.7467 0 0 
 
 0 0 0 15 
 15 1.7647 0 20 
   20  15 
 9.9706 6 .5294  6 11.0294 0  0 0
3
 0 1.7467 0  20
B BD  1.282 X 10  3.2647 6 1.2353  6 5.7353 03.333 X 10 3 
T

 3.9706 14 0.5294  14 5.0294 0  0 0  20 0 
 
4 3 4 0 0 3  15 1.7467 0 0 
 
 0 0 0 15 
 223.7978  139.4118  85.7611 79.4118  155.32 60 
  139.412 325 70.588  280 100.558  45

  85.7612 70.588 82.18  10.588 10.211  60
B BD  4.27333
T

 79.412  280  10.58 280  100.588 0 
 155.3202 100.5882 10.1210  100.588 175.5621 0 
 
 60  45  60 0 0 45 

Substitute these values

 223.7978  139.4118  85.7611 79.4118  155.32 60 
  139.412 325 70.588  280 100.558  45

  85.7612 70.588 82.18  10.588 10.211  60
K 1  2 XX 56.667 X 150 X 4.27333 
 79.412  280  10.58 280  100.588 0 
 155.3202 100.5882 10.1210  100.588 175.5621 0 
 
 60  45  60 0 0 45 

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 25

 223.7978  139.4118  85.7611 79.4118  155.32 60 

  139.412 325 70.588  280 100.558  45

  85.7612 70.588 82.18  10.588 10.211  60
K 1  228.2246 X 103 X  
 79.412  280  10.58 280  100.588 0 
 155.3202 100.5882 10.1210  100.588 175.5621 0 
 
 60  45  60 0 0 45 
 51.076  31.817  19.573 18.124  35.448 13.693 
 31.817 74.173 16.10  63.903 22.957  10.270

  19.573 16.110 18.755  2.46 2.310  13.693
K 1  106 X  
 18.124  63.903  2.416 63.903  22.957 0 
 35.448 22.957 2.310  22.957 40.068 0 
 
 13.693  10.270  13.693 0 0 10.270 

For element:2

Co-ordinates:
At node 1:
r1=50mm
z1=0mm
At node 2:
r2=70mm
z2=0mm
At node 3:
r3=70mm
z3=15mm
r  r  r3 50  70  70
we know that, r= 1 2 =
3 3
r=63.3333mm

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 26

z1  z 2  z 3 0  0  15
z= =
3 3
z=5mm
area of the triangle element=1/2xbreath x height=1/2x(70-50)x15
A=150mm
We know that matrix for axisymmetric triangular element (1)
[K]1=2πrA [B]T [D] [B]
1  v v v 0 
 v 1 v v 0 
E 
Stress-strain relationship matrix, [D]=
1  v 1  2v   v v 1 v 0 
1  2v 
 0 0 0
2 
1  0.3 0.3 0. 0 
 
2 x10 5  0.3 1  0.3 0.3 0 
[D]=  0.3 0 .3 1  0.3 0 
0.52
 1  2 x0.3 
 0 0 0
2 

1  0.3 0.3 0. 0 
 
2 x10 5  0.3 1  0.3 0.3 0 
  0.3 0 .3 1  0 .3 0 
0.52
 1  2 x0.3 
 0 0 0
2 

We know that,

Strain-displacement matrix,
 
 1 0 2 0 3 0
 
3  3z
B  1  1   1   1 z 0 2
 2 
 2z
0  3  0
2A  r r r r r r 
 0 1 0 2 0 3 
 
 1 1 2 2 3 3 
 1  r2 z 3  r3 z 2  70 X 15  70 X 0
 1  1050mm 2
 2  r3 z1  r1 z 3  70 X 0  50 X 15
 2  750mm 2

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 27

 3  r1 z 2  r2 z1  50 X 0  70 X 0
 3  0mm 2
 1  z 2  z 3  0  15
 1  15mm
 2  z 3  z1  15  0
 2  15mm
 3  z1  z 2  0  0
 3  0mm
 1  r3  r2  70  70
 1  0mm
 2  r1  r3  50  70
 2  20mm
 3  r2  r1  70  50
 3  20mm
1 1z 1050
 1    (15)  0  1.5647mm
r r 63.333
2  z  750  20 X 5
 2  2   15   1.5647mm
r r 63.367 63.6336
3  z 20 X 5
 3  3  0  0   1.579mm
r r 63.33
Substitute these values.

 
  15 0 0 0 15 0
B  1 1.7647 0 1.7647 0 1.7647 0 
2 X 150 0 20 0  20 0 0
 
 20  15  20 0 0 15
 
  15 0 0 0 15 0
B  3.333 X 10 3 1.7647 0 1.7647 0 1.7647 0 
0 20 0  20 0 0
 
 20  15  20 0 0 15

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 28

0.7 0.3 0. 0  
0.3 0.7 0.3 0    15 0 0 0 15 0
 
BD  384.6153 X 10 3 0.3 0.3 0.7 0  X 3.333 X 10 3 1.7647 0 1.7647 0 1.7647 0 
 0  0 20 0  20 0 0
 
0 0 0     
 0.2  20 15 20 0 0 15
 
 9.9706 6 .5294  6 11.0294 0
BD  1.282 X 10 3  3.2647 6 1.2353  6 5.7353 0
 3.9706 14 0.5294  14 5.0294 0
 
4 3 4 0 0 3

We know that,
 
  15 0 0 0 15 0
B  3.333 X 10 3 1.7647 0 1.7647 0 1.7647 0 
0 20 0  20 0 0
 
 20  15  20 0 0 15
 15 1.7647 0 20 
 0 0 20  15 

 0 1.7467 0  20
BT  3.333 X 10 3  
 0 0  20 0 
 15 1.7467 0 0 
 
 0 0 0 15 
 15 1.7647 0 20 
   0 0 20  15 
 9.9706 6 .5294  6 11.0294 0 
3
 0 1.7467 0  20
B BD  1.282 X 10  3.2647 6 1.2353  6 5.7353 03.333 X 10 3 
T

 3.9706 14 0.5294  14 5.0294 0  0 0  20 0 
 
4 3 4 0 0 3  15 1.7467 0 0 
 
 0 0 0 15 
 223.7978  139.4118  85.7611 79.4118  155.32 60 
  139.412 325 70.588  280 100.558  45

  85.7612 70.588 82.18  10.588 10.211  60
B BD  4.27333
T

 79.412  280  10.58 280  100.588 0 
 155.3202 100.5882 10.1210  100.588 175.5621 0 
 
 60  45  60 0 0 45 

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 29

Substitute these values

 223.7978  139.4118  85.7611 79.4118  155.32 60 
  139.412 325 70.588  280 100.558  45

  85.7612 70.588 82.18  10.588 10.211  60
K 1  2 XX 56.667 X 150 X 4.27333 
 79.412  280  10.58 280  100.588 0 
 155.3202 100.5882 10.1210  100.588 175.5621 0 
 
 60  45  60 0 0 45 
 223.7978  139.4118  85.7611 79.4118  155.32 60 
  139.412 325 70.588  280 100.558  45

  85.7612 70.588 82.18  10.588 10.211  60
K 1  228.2246 X 103 X  
 79.412  280  10.58 280  100.588 0 
 155.3202 100.5882 10.1210  100.588 175.5621 0 
 
 60  45  60 0 0 45 
 51.076  31.817  19.573 18.124  35.448 13.693 
 31.817 74.173 16.10  63.903 22.957  10.270

  19.573 16.110 18.755  2.46 2.310  13.693
K 1  106 X  
 18.124  63.903  2.416 63.903  22.957 0 
 35.448 22.957 2.310  22.957 40.068 0 
 
 13.693  10.270  13.693 0 0 10.270 

Two Marks Question and Answers.

UNIT-III
1. What is axisymmetric element?
Many three dimensional problem in engineering exhibit symmetry about an axis of rotation such type of
problem are solved by special two dimensional element called the axisymmetric element
2. What are the conditions for a problem to be axisymmetric? (Apr 2010)
The condition to be axi-symmetrix is as follows
 The problem domain must be symmetric about the axis of revolution
 All boundary condition must be symmetric about the axis of revolution
 All loading condition must be symmetric about the axis of revolution
3. Give the stiffness matrix equation for an axisymmetrc triangular element.
Stiffness matrix  K    B   D  B  2 rA
T

A-Area
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 30

r-Radius

4. What is the purpose of Isoparametric element?

It is difficult to represent the curved boundaries by straight edges finite elements. A large number of finite
elements may be used to obtain reasonable resemblance between original body and the assemblage.
5. Define super parametric element?

If the number of nodes used for defining the geometry is more than of nodes used for defining the
displacement is known as super parametric element
6. What is meant by Isoparametric element?
If the number of nodes used for defining the geometry is same as number of nodes used for defining the
displacement is known as Isoparametric element
7. Define sub parametric element?
If the number of nodes used for defining the geometry is less than of nodes used for defining the displacement
is known as sub parametric element
8. Is beam element an Isoparametric element?
Beam element is not an Isoparametric element since the geometry and displacement are defined by different
order interpretation functions.
9. What is meant by degrees of freedom?
When the force or reaction act at nodal point node is subjected to deformation. The deformation includes
displacement rotation, and or strains. These are collectively known as degrees of freedom.

10. State the principles of virtual energy?

A body is in equilibrium if the internal virtual work equals the external virtual work for the every
kinematically admissible displacement field.
11. What is homogeneous form?
When the specified values of dependent variables is zero, the boundary condition are said to be
homogeneous.
12. What is non-homogeneous form?
When the specified values of dependent variables are non-zero, the boundary conditi said to be non-
homogeneous.
13. Write down the shape functions for 4 noded rectangular elements using natural
coordinate system?

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 31

14. Write down the stress strain relationship matrix for plane stress conditions?

E-youngs modulus
v-poisson ratio
15. Write down Jacobian matrix for 4 noded quadrilateral elements?

16. What is Aspect ratio?

It is defined as the ratio of the largest dimension of the element to the smallest dimension. In many
cases, as the aspect ratio increases the in accuracy of the solution increases. The conclusion of many
researches is that the aspect ratio should be close as unity as possible.
17. What is higher order element ? (Nov 2011)
When the finite element approximation is gradually refining the mesh in comparing the solution then it is called
higher order elements
18. What are the four basic sets of elasticity equations? (May 2012)
Elasticity equations are used for solving structural mechanics problems. These equations must be satisfied if an
exact solution to a structural mechanics problem is to be obtained. There are four basic sets of elasticity equations.
They are:
 Strain- Displacement relationship equations
 Stress-Strain relationship equations
 Equilibrium equations
 Compatibility equations
19. What are the types of Non-Linearity ? (May 2012)

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 32

 Non-Linearity in material behavior from point to point

 Non-Linearity in loading deformation relation
 Geometric non-Linearity
 Change in boundary condition for different loading
20. What are higher order elements and why are they preferred? (April 2011)
 For any element, if the interpolation polynomial is the order of two or more, that element is known as
higher order elements
 It is used to represent the curved boundaries
 The number of elements are reduced when compared with straight edge elements to model geometry
21. Give four applications where axi-symmetric elements can be used? (April 2011)
 Pressure Vessels
 Rocket Castings
 Cooling Towers
 Sub-Marine hulls
 Belleville springs

22. What is 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.plane stress if the stress vector is zero
across a particular surface. When that situation occurs over an entire element of a structure, as is
often the case for thin plates, the stress analysis is considerably simplified, as the stress state can
be represented by a tensor of dimension 2
23. Define plane strain analysis.( AU-NOV/DEC-2012)
Plane strain is defined to be a state of strain in which the strain normal to the xy plane and the
shear strains are assumed to be zero.Plane strain is applicable to rolling, drawingand forging
where flow in a particular direction is constrained by the geometry of the machinery, e.g. a well-
lubricated die wall.
24. Write down the stress-strain relationship matrix for plane stress conditions.
For plane stain problem, stress strain relationship matrix is,
 
1 v 0 
D  E v 1 0 
1  v2  1 v
0 v 
 2 
Where E=young’s modulus
V=Poisson’s ratio.
25. Write down the stress-strain relationship matrix for plane strain condition.

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 33

For plane strain problems, stress –strain relationship matrix is,

 
1  v  v 0 
 
D  E
 v 1  v  0 
1  v 1  2v    1  2v 
 0 0  
  2 
3 noded triangular elements – four noded rectangular elements – higher order elements
26. What is axisymmetric element?(AU-APR/MAY-2014)
Many three dimensional problems in engineering exhibit symmetry about an axis of
rotation. Such types of problems are solved by a special two dimensional element called as
axisymmetric element. The axisymmetric problemor elements which is formulated from the three
dimensional element into two dimensional element because of the symmetric axis. The
symmetric axis which explains the side a is equal to the side a’ which refers us to considerthe
three dimensional problem into two dimensional one.
27. What are the conditions for a problem to be axisymmetric?
Conditions for a problem to be axisymmetric:
1.The problem domain must be symmetric about the axis of revolution.
2. All boundary conditions must be symmetric about the axis of revolution.determination of
boundary condition is important in the axisymmetric probles.
3. All loading conditions must be symmetric about the axis of revolution.
28. Write down the displacement equation for an axisymmetric triangular element.

 ur , z   N1 0 N2 0 N3 0
Displacement function, u(r,z)=j  
wr , z   0 N1 0 N2 0 N 3 

r,z represents the direction

u – represents displacement function in x axis

w-represents displacement function in y direction.

29. Write down the shape function for an axisymmetric triangular element. (AU-APR/MAY-
2013)
  1 r   1 z
Shape function, N1= 1
2A
  2r   2 z
N2= 2
2A
  3r   3 z
N3= 3
2A
Where,  1  r2 z 3  r3 z 2

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

UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 34

 2  r3 z1  r1 z 3
 3  r1 z 2  r2 z1
1  z 2  z 3
 2  z 3  z1
 3  z1  z 2
 1  r3  r2
 2  r1  r3
 3  r2  r1
30. Give the strain-displacement matrix equation for an axisymmetric triangular element.
 
 1 0 2 0 3 0
 
B  1  1   1   1 z 0  2   2   2 z 0  3   3   3 z 0 
2A  r r r r r r 
 0 1 0 2 0  3

 
 1 1 2 2 3 3 
r  r  r3
Where, co-ordinate, r= 1 2
3
z  z 2  z3
Co-ordinate, z= 1
3
Where, 1 ,  2  3 ,  1 ,  2  3 , 1 ,  2 ,  3 -co-ordinates
31. What are the types of non-linearity?
Types of non-linearity:
(i) Non linearity in material behavior from point to point.
(ii) Non-linearity loading-deformation reaction.
(iii) Geometric non-linearity.
(iv) Change in boundary condition for different loading.
natural coordinates and coordinate transformations – triangular and quadrilateral
elements
32. Write down the stress-strain relationship matrix for an axisymmetric for an axisymmetric
triangular element. (AU-APR/MAY-2012)
1  v v v 0 
 v 1 v v 0 

Stress-strain relationship matrix, D  
E
1  v 1  2v   v v 1 v 0 
1  2v 
 0 0 0
2 
Where, E,-Young;s modulus
v-poisson’s ratio.
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 35

33. Give stiffness matrix equation for an axisymmetric triangular element

Stiffness matrix,
K   2rABT DB
Where, co-ordinate
r  r  r3
r= 1 2
3
A-area of the triangular element matrix.
34. What are the ways in which a three dimensional problem can be reduced to a two
dimension approach.
Reducing three dimensional element into two dimensional element:
(i) Plane stress; one dimensional is too small when compared to other two ic
(ii) dimensions.
(iii) Plane strain: one dimensional is too large when compared to other two dimensions
(iv) Axisymmetric: geometry about the axis.
35. What are natural co-ordinates? (AU-NOV/DEC-2013)
Natural coordinates:
A simple natural coordinate system is a kind of local coordinate system that permits the
specification of a point within the element by a set of dimensionless numbers whose magnitude
varies from -1 to +1. The actual finite element say two dimensional element which may be a
rectangle or quadrilateral or any curved sided and specified any global coordinates system is
represented in natural coordinate system.
36. What are the advantages of natural co-ordinates?
Advantages of natural coordinates:
i. Natural coordinate system enable us to formulate shape function easily.
ii. The shape function for the complicated structure also very easy in the natural
coordinate system.
iii. The natural coordinate system helps us to solve the rectangle, quadrilateral, curve
sided element in a simple way.
37. What is meant by axisymmetric solid? (AU-NOV/DEC-2012)
Axisymmetric solid:
In some three dimensional solid like cylinder, flywheel, turbine disc the material content
is symmetric with respect to their axes. The material content is equal in equal distance of
opposite sides. Hence the stress developed, displacement produced are considered as symmetric.
Such solids are known as axisymmetric solids.
Due to this axisymmetric condition these three dimensional solids can be treated as two
dimensional solids for analysis.
38. Specify the machine component related with axisymmetric concept. (AU-NOV/DEC-2011)
Machine component related with axisymmetric concept:
Pressure vessels- the stress deformation due to pressure are analysed here.
Pressure cylinders- the cylinders which carry internal or external loads.
ME6603 FEA Prepared By Er.M.GOPALA KRISHNAN/AP/Mech/VRSCET
UNIT-IV / TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS P a g e | 36

Flywheel – energy storing element which may deformed during torque transmission.
Turbine discs- the rotary disc either pressure due to angle or force.
Soil mass – structural problem defines base plan of constructions.
39. Write short note on axisymmetric formulation.
Axisymmetric formulation:
 u 
 r 
e r   u 
e   
   u u w  u w  
T
  r
{e}         
e
  
z w   r r z  z r  
 rz   z 
 u w 
  
 z r 
40. Briefly discuss about finite element modeling for axisymmetric formulation.
Finite element modeling for axisymmetric formulation:
For an axisymmetric element the two dimensional region defined by the revolving area
can be divided into triangular elements or quadrilateral elements.
{F} = [K] {u}
[k ]   [ B]T [ D][ B]dv
v

[ D]  stress  strain matrix

[B] - strain displacement matrix
[k] - stiffness matrix
{F} – Force
{u} - displacement

41. What are the non zero strain and stress components of axisymmetric element?
Non zero strain and stress components of axisymmetric element:

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