Professional Documents
Culture Documents
Untitled
Untitled
2
Please read this disclaimer before proceeding:
This document is confidential and intended solely for the educational purpose of
RMK Group of Educational Institutions. If you have received this document
through email in error, please notify the system manager. This document
contains proprietary information and is intended only to the respective group /
learning community as intended. If you are not the addressee you should not
disseminate, distribute or copy through e-mail. Please notify the sender
immediately by e-mail if you have received this document by mistake and delete
this document from your system. If you are not the intended recipient you are
notified that disclosing, copying, distributing or taking any action in reliance on
the contents of this information is strictly prohibited.
3
ME8692 FINITE ELEMENT ANALYSIS
4
Table of Contents
2 Pre Requisites 06
3 Syllabus 07
4 Course outcomes 08
6 Lecture Plan 10
8 Lecture Notes 13
9 Assignments 45
10 Part A Q & A 47
11 Part B Qs 49
15 Assessment Schedule 56
5
I - Course Objectives
01 To introduce the concepts of Mathematical Modeling of
Engineering Problems.
02 To appreciate the use of FEM to a range of Engineering
Problems.
GE8292 Engineering
Mechanics (Semester 2)
CE8395 Strength of
Materials for Mechanical
Engineers (Semester -4)
5
III - Syllabus
UNIT I INTRODUCTION 9
Historical Background – Mathematical Modeling of field problems in Engineering
– Governing Equations – Discrete and continuous models – Boundary, Initial and
Eigen Value problems– Weighted Residual Methods – Variational Formulation of
Boundary Value Problems – Ritz Technique – Basic concepts of the Finite Element
Method.
TOTAL PERIODS : 45
7
IV - Course Outcomes
(COs)Course Outcome (CO) - Statement
CO. No. Knowledge
Level
Course Augmentation
Units I II III IV V
Assignement 09.03.2021 23.03.2021 10.04.2021 20.04.2021 07.05.2021
Quiz
8
VI - COs-POs Matrices
POs CO PO1 PO 2 PO 3 PO 4 PO 5 PO 6 PO 7 PO 8 PO 9 PO 10 PO 11 PO 12
COs
C305
.1 CO1 3 3 3
C305
.2 CO2 3 3 2
C305
.3 CO3 3 3 3
C305
.4 CO4 3 3 3
C305
.5 CO5 3 3 2
C305
.6 CO6 3 3 3
1: Slight (Low) 2: Moderate (Medium) 3: Substantial (High)
9
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.
LECTURE PLAN
Session Mode of
Topics to be covered Reference
No. Delivery*
Second Order 2D Equations involving Scalar
1 1,2
Variable Functions
2 Variational formulation 1,2
3 Finite Element formulation 1,2
4 Triangular elements Chalk and 2,3
5 Shape functions and element matrices and Talk/Online 2,3
vectors Webinars
6 Application to Field Problems 1,2
7 Thermal problems 1,2
8 Thermal problems 1,2
9 Torsion of Non circular shafts 1,2
10 Torsion of Non circular shafts 1,2
11 Quadrilateral elements 1,2
12 Higher Order Elements 1,2
*Chalk and Talk, PPT / OHP, Quiz, Group Discussion, Student, Seminar,
Case Study, Field Study, Video (NPTEL/EDUSAT/ONLINE), Crossword
Puzzle, Tutorial Session
10
ACTIVITY BASED LEARNING
13
ACTIVITY BASED LEARNING
Answer
13
UNIT IV: TWO DIMENSIONAL VECTOR VARIABLE PROBLEMS
LECTURE NOTES
Introduction—Dimensionality of a Problem
In the first case [Fig.a], the structure is very thin as compared to its other
two dimensions and is loaded in its own plane, e.g. a rotating impeller wheel. Thus the
top and bottom faces are free from any normal stress z and since the thickness is
small, we assume z is zero throughout. So also, the shearing stresses 𝜎𝑥 , 𝜎𝑦 , 𝜏𝑥 𝑦 .
Thus the only nonzero stresses are 𝜏𝑥 𝑧 𝜏𝑦 𝑧 . The stress-strain constitutive relation for a
linear, elastic, isotropic material under such plane-stress condition is given by
𝜎𝑥 1 𝑣 0 𝜎𝑥
𝐸 𝑣 1 0 ] ( { 𝜎𝑦 } − {𝜀}0 )
{ 𝜎𝑦 } = [
𝜏𝑥 𝑦 1 − 𝑉2 1−𝑣 𝛾𝑥 𝑦
0 0
2
The problem situation depicted in the second case (Fig. b) is the other extreme
of the first case, where the length of the member is so long that it is reasonable to
assume thelongitudinal strain 𝜀𝑧 to be zero, e.g. torsion of a long uniform shaft. So also
the shearing strains𝜏𝑥 𝑧 𝜏𝑦 𝑧 Thus the only nonzero strains are 𝜀𝑥 , 𝜀𝑦 , 𝜏𝑥 𝑦 . The stress-
strain constitutive relation for a linear, elastic, isotropic material under such plane-
strain condition is given by
𝜎𝑥 1−𝑣 𝑣 0 𝜀𝑥
𝐸 1−𝑣 0 ] ( { 𝜀𝑦 } − {𝜀}0 )
{ 𝜎𝑦 } = [ 𝑣
𝜏𝑥 𝑦 (1 + 𝑉)(1 − 2𝑣) 1 − 2𝑣 𝛾𝑥 𝑦
0 0
2
The problem situation depicted in the third case (Fig.c) represents an important
class of problems, commonly described as axisymmetric problems, e.g. pressure
vessels. If the geometry of the structure, material properties, support conditions and
loading are all axisymmetric, then it is reasonable to assume that the resulting
deformation will also be axisymmetric, i.e., it does not vary along the circumferential
direction. Thus, it is sufficient to determine the variation of the unknown field variable
(e.g. deflection, temperature) in any one plane as shown in the figure. In practical
pressure vessels, there may be minor geometrical features which cause deviation from
perfect axial symmetry such as helical screw threads on an end-cap. It may be
worthwhile to ignore the small helix angle and model the problem as purely 2-d
axisymmetric. However, sometimes a perfect axisymmetric structure may be subjected
to nonaxisymmetric loading such as earthquake loading. It may still be useful to
modify a two-dimensional axisymmetric analysis rather than going for a full-fledged
three- dimensional analysis. The non-axisymmetric loading, expressed as a function of
circumferential coordinate , has a periodicity 2. It can be expressed as a Fourier
series. To each harmonic of load, the resulting deflection also varies as a sine or
cosine function and it can still be analysed as a two-dimensional problem. The
Fourier harmonics of deflection are superposed to find the total response. This
technique is known as semi-analytical FEM. The nonzero stresses and strains in a
perfect axisymmetric situation are related as follows:
v v
1 0
r 1− v 1− v r
0
v
(1− v)
1 0
z = 1− v z −
(1+ V)(1− 2v)
1 0
rz
1− 2v rz
sysmmetric 2(1− v)
Each node is capable of independent motion along the x and y axes and,
typically, we have u and v as the independent nodal degree of freedom (d.o.f.) for
plane elasticity problems. However, there is only one temperature d.o.f. for 2-d heat
transfer problems. Simple lower order elements have straight sides and nodes only at
the vertices. More complex higher order elements permit curved edges (so as to model
curved geometric feature in the domain such as a circular hole) and mid-side or even
internal nodes.
It has straight sides and three nodes, one at each vertex. The nodes have
coordinates (x1, y1), (x2, y2), (x3, y3) in the global Cartesian coordinate frame OXY as
shown in the figure. Each node has just temperature d.o.f. T or two d.o.f., viz., u
and v, the translations along global X and Y azes, respectively. In what follows, we
derive the shape (interpolation) functions for the temperature field and use the same for
the displacement field.
Considering that this expression should reduce to the nodal temperature at the
nodal points, we have
𝑇1 = 𝑐0 + 𝑐1 𝑥1 + 𝑐2 𝑦1
𝑇2 = 𝑐0 + 𝑐1 𝑥2 + 𝑐2 𝑦2 }
𝑇3 = 𝑐0 + 𝑐1 𝑥3 + 𝑐2 𝑦3
𝛼1 + 𝛽1 + 𝛾1 𝑦 𝛼2 + 𝛽2 + 𝛾2 𝑦 𝛼3 + 𝛽3 + 𝛾3 𝑦
𝑇(𝑥, 𝑦) = ( ) + 𝑇1 ( ) 𝑇2 ( ) 𝑇3
2△ 2△ 2△
where
𝛼1 = 𝑥2 𝑦3 − 𝑥3 𝑦2 , 𝛽1 = 𝑦2 − 𝑦3 , 𝛾1 = 𝑥3 − 𝑥2
1 𝑥1 𝑦1
2𝛥 = [1 𝑥2 𝑦2 ] = 2(Area of
1 𝑥3 triangle123)
𝑦3
In our standard finite element notation, we write this as
1
𝑁𝑖 = {𝛼 + 𝛽𝑖 𝑥 + 𝛾𝑖 𝑦}, 𝑖 = 1,2,3
2𝛥 𝑖
If this element is used to model structural mechanics problems, each node will
have two d.o.f., viz., u and v, and we can write the displacements at any interior
point in terms of thenodal displacements using the same shape functions as
𝑢1
𝑣1
𝑢 𝑁 0 𝑁2 0 𝑁3 0 𝑢2 [𝑁]{𝛿}𝑒
{ 𝑣 } = [ 01 𝑁3 ] 𝑣2 =
𝑁1 0 𝑁2 0
𝑢3
{ 𝑣3 }
We observe that this element may only permit a linear variation of the unknown
field variable within the element and thus derivatives such as strains and heat flux will
be constant throughout the element. In structural mechanics, this element is
therefore popularly known as the constant strain triangle (CST). We now discuss a
typical rectangular element.
We observe that as we have four nodes now, we are able to take a fourth term
in our polynomial field. Considering that this expression should reduce to the nodal
temperatures at the nodal points we have, for a rectangular element of size (l × h)
as shown in above Figure,
T1 = c0, T2 = c0 + c1l,
𝑥 𝑦 𝑥𝑦 𝑥 𝑥𝑦 𝑥𝑦 𝑦 𝑥𝑦
𝑇(𝑥, 𝑦) = (1 − − + ) 𝑇1 + ( − ) 𝑇2 + ( ) 𝑇3 + ( − ) 𝑇4
ℓ ℎ ℓℎ ℓ ℓℎ ℓℎ ℎ ℓℎ
In our standard finite element notation, we write this as
𝑇1
𝑇2
𝑇(𝑥, 𝑦) = [𝑁1 𝑁2 𝑁3 𝑁4 ] { }
𝑇3
𝑇4
If this element is used to model structural mechanics problems, each node will
have two d.o.f., viz., u and v, and we can write the displacements at any interior
point in terms of thenodal displacements using the same shape functions as follows:
𝑢1
𝑣1
𝑢2
𝑢 𝑁 0 𝑁2 0 𝑁3 0 𝑁4 0 𝑣2 [𝑁]{𝛿}𝑒
{ 𝑣 } = [ 01 𝑁4 ] 𝑢3 =
𝑁1 0 𝑁2 0 𝑁3 0
𝑣3
𝑢4
{ 𝑣4 }
Since we have six nodes, we can fit a higher degree polynomial for the variation
of the unknown field variable. To be precise, following our procedure, we can determine
up to six coefficients ci(i = 0, 1, …, 5). Thus we assume that the temperature field
over the element is given by the complete quadratic polynomial
Considering that this expression should reduce to the nodal temperature at the
nodal points, we have
𝑇1 = 𝑐0 + 𝑐1 𝑥1 + 𝑐2 𝑦1 + 𝑐3 𝑥12 + 𝑐4 𝑥1 𝑦1 + 𝑐5 12
𝑦2 = 𝑐0 + 𝑐1 𝑥2 + 𝑐2 𝑦2 + 𝑐3 𝑥22 + 𝑐4 𝑥2 𝑦2 + 𝑐5 𝑦22
𝑇
𝑇3 = 𝑐0 + 𝑐1 𝑥3 + 𝑐2 𝑦3 + 𝑐3 𝑥32 + 𝑐4 𝑥3 𝑦3 + 𝑐5 𝑦32
𝑇4 = 𝑐0 + 𝑐1 𝑥4 + 𝑐2 𝑦4 + 𝑐3 𝑥42 + 𝑐4 𝑥4 𝑦4 + 𝑐5 𝑦42
𝑇5 = 𝑐0 + 𝑐1 𝑥5 + 𝑐2 𝑦5 + 𝑐3 𝑥52 + 𝑐4 𝑥5 𝑦5 + 𝑐5 𝑦52
𝑇6 = 𝑐0 + 𝑐1 𝑥6 + 𝑐2 𝑦6 + 𝑐3 𝑥62 + 𝑐4 𝑥5 𝑦5 + 𝑐5 𝑦62}
Generic Relations
We will now summarise certain generic relations which will be useful for all the
2-d finite elements for structural mechanics. For a two-dimensional structural problem,
each point on the structure may have two independent displacements, viz., u and v,
along the two Cartesian coordinates X and Y, respectively. The strain-displacement
relations are given by
𝜕𝑢
𝜀𝑥 =
𝜕𝑥
𝜕𝑣
𝜀𝑦 =
𝜕𝑦
𝜕𝑢 𝜕𝑢
𝛾𝑥 𝑦 = +
𝜕𝑦 𝜕𝑥
For axisymmetric problems, we can consider the r-z plane (analogous to x-y
plane for plane elasticity) as shown in Figure.
We observe that, even though the deformation is axisymmetric, there can be
four independent, nonzero strains (as a change in radius automatically leads to
circumferential strain 𝜀 0):
𝜕𝑢
𝜀𝑟 =
𝜕𝑟
𝑢
𝜀𝜃 =
𝑟
𝜕𝑣
𝜀𝑧 =
𝜕𝑧
𝜕𝑢 𝜕𝑢
𝛾𝑟 𝑧 = +
𝜕𝑧 𝜕𝑟
{𝜀} = [𝐵]{𝛿}𝑒
where the size of [B] and {} will be dictated by the element and [ B] will
contain the derivatives of the shape functions. The stress-strain relations for plane
stress specialised from the 3-d Hooke’s law, are given by
𝜎𝑥 1 𝑣 0 𝜀𝑥
𝐸 0 𝜀
] ( { 𝑦 } − {𝜀}0 ) + {𝜎}0
{ 𝜎𝑦 } = [𝑣 1
𝜏𝑥 𝑦 1 + 𝑉2 0 2
0 1 − 𝑣 /2 𝛾𝑥 𝑦
where {}0 and {}0 refer to the initial strain and stress, respectively.
𝜎𝑟 1 𝑣/1 − 𝑣 𝑣/1 − 𝑣 0 𝜀𝑟
𝜎𝜃 (1 − 𝑣)𝐸 1 𝑣/1 − 𝑣 0 𝜀𝜃
{ 𝜎 } = ({ 𝜀 }
𝑧 (1 + 𝑣)(1 − 2𝑣) 1 0 𝑧
𝜏𝑟 [𝑠𝑦𝑚𝑚𝑒𝑡𝑟𝑖 1 − 2𝑣/2(1 − 𝑣)] 𝛾𝑟 𝑧
− {𝜀}0 ) + {𝜎}0
1 1
𝑈= ∫ {𝜀}𝑇 {𝜎}𝑑𝑣 = ∫ {𝜀}𝑇 ([𝐷]({𝜀} − {𝜀}0 ) + {𝜎 0 }𝑑𝑣)
2 𝑣 2 𝑣
1
= ∫ ( {𝜀}𝑇 [𝐷]{𝜀} − {𝜀}𝑇 [𝐷]{𝜀} − {𝜀}0 + {𝜀}0 {𝜎 0 }) 𝑑𝑣
𝑣 2
Element level strain energy is given as
1 1
𝑈𝑒 = ∫ {𝛿}𝑒 𝑇 [𝐵]𝑇 [𝐷] [𝐵]{𝛿}𝑒 𝑑𝑣 − ∫ {𝛿}𝑒 𝑇 [𝐵]𝑇 [𝐷]
[𝐵]{𝜀}0 𝑑𝑣 2 2
𝑢
{𝛿} = { 𝑣 }
𝑞
{𝑞} = {𝑞𝑥 }
𝑣
𝑇 𝑇 𝑒𝑇
𝑉 𝑒 = − ∫ {𝛿}𝑒 [𝑁] 𝑇 {𝑞𝑣 }𝑑𝑣 − ∫ {𝛿} 𝑒 [𝑁]𝑇 {𝑄𝑠 }𝑑𝑠 − ∑ {𝛿𝑖 } [𝑁]𝑇𝑖 {𝑃𝑖 }
𝑣 𝑠
where [N]i is the value of shape functions at the location i. Thus the total potential of
an element can be written as
1 𝑇 𝑇
∏𝑒𝑝 = ∫ {𝛿}𝑒 [𝐵]𝑇 [𝐷] [𝐵]{𝛿}𝑒 𝑑𝑣 − ∫ {𝛿}𝑒 [𝐵]𝑇 [𝐷] [𝐵]{𝜀}0 𝑑𝑣
2
𝑇 𝑇 𝑒𝑇
+ ∫ {𝛿} 𝑒 𝑇 [𝐵]𝑇 [𝜎]0 𝑑𝑣 − ∫ {𝛿}𝑒 [𝑁]𝑇 {𝑞𝑣 }𝑑𝑣 − ∫ {𝛿}𝑒 [𝑁] 𝑇 {𝑄𝑠 }𝑑𝑠 − ∑ {𝛿𝑖 } [𝑁]𝑇𝑖 {𝑃𝑖 }
𝑣 𝑠
Since {}e is a vector of nodal d.o.f., it can be taken outside the integration. Thus we
get
0
1 𝑇 𝑇 𝑇
∏𝑒𝑝 = {𝛿} 𝑒 ∫ [𝐵]𝑇 [𝐷][𝐵] 𝑑𝑣){𝛿} 𝑒 − {𝛿} 𝑒 ∫ [𝐵]𝑇 [𝐷]{𝜀} 𝑑𝑣 + {𝛿}𝑒 ∫ [𝐵]𝑇 {𝜎}0 𝑑𝑣
2 𝑣 𝑣 𝑣
𝑇 𝑇 𝑒𝑇
−{𝛿}𝑒 ∫ [𝑁]𝑇 {𝑞𝑣 }𝑑𝑣 − ∫ {𝛿}𝑒 [𝑁] 𝑇 {𝑞𝑠 }𝑑𝑠 − ∑ {𝛿𝑖 } [𝑁]𝑖𝑇 {𝑃𝑖 }
𝑣 𝑠
Defining the element stiffness matrix [ k]e and the load vector {f }e as
{𝑓}𝑒 = − ∫ [𝐵]𝑇 [𝐷]{𝜀}0 𝑑𝑣 − ∫ [𝐵]𝑇 {𝜎}0 𝑑𝑣 + ∫ [𝑁]𝑇 {𝑞𝑣 }𝑑𝑣 + ∫ [𝑁]𝑇 {𝑞𝑠 }𝑑𝑠 + ∑[𝑁]𝑇𝑖 {𝑃𝑖
} 𝑣 𝑣 𝑣 𝑠
we can rewrite as
𝑒
1 𝑇 𝑇 𝑇
∏𝑒𝑝 = {𝛿} 𝑒 ∫ [𝑘]𝑒 {𝛿}𝑒 − {𝛿} 𝑒 {𝑓}
2 𝑣
Once all the element matrices are ‘assembled’ together, we obtain the total
potential of the system as
1
∏𝑝 = ∑∏𝑝𝑒 = {𝛿} 𝑇 [𝐾]{𝛿} − {𝛿}𝑇 {𝐹}
2
where the global stiffness matrix of the structure [ K] and the global load vector {F} are
given by
𝑁𝑂𝐸𝐿𝐸𝑁
[𝐾] = ∑ {𝑘}𝑒
𝑛 =1
𝑁𝑂𝐸𝐿𝐸𝑁
[𝐹] = ∑
{𝑓}𝑒
1
NOELEM refers to the number of elements and {} contains all the nodal d.o.f.
variables for the entire finite element mesh. The summations indicating assembly imply
that the individual element matrices have been appropriately placed in the global
matrices following the standard procedure of assembly.
Using the Principle of stationary total potential we set the total potential
stationary with respect to small variations in the nodal d.o.f., i.e.*
𝜕∏𝑝 = 0
𝜕{𝛿}𝑇
[𝐾]{𝛿} = {𝐹}
We will now discuss each individual 2-d element and describe how the element
stiffness matrix and nodal load vectors can be obtained.
This is one of the simplest and historically one of the first finite elements to be
formulated. A typical element is shown in Figure.
1
𝑁1 = (𝛼 + 𝛽1 𝑥 + 𝛾1 𝑦)
2𝛥 1
1
𝑁2 = (𝛼 + 𝛽2 𝑥 + 𝛾2 𝑦)
2𝛥 2
1
𝑁3 = (𝛼 + 𝛽3 𝑥 + 𝛾3 𝑦)
2𝛥 3
𝜕𝑁 1 0 𝜕𝑁2 𝜕𝑁
0 3 0
𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕𝑁 1 𝜕𝑁 2 𝜕𝑁
[𝐵] = 0 0
𝜕𝑦 𝜕𝑦 0 3
𝜕𝑦
𝜕𝑁 1 𝜕𝑁2 𝜕𝑁 𝜕𝑁 𝜕𝑁
3 3
𝜕𝑁 1 𝜕𝑥 𝜕𝑦 2
𝜕𝑥 𝜕𝑦 𝜕𝑥 }
{ 𝜕𝑦
Thus
𝛽1 0 𝛽2 0 𝛽3 0
1
[𝐵] = {0 𝛾1 0 𝛾2 0 𝛾3 }
2𝛥
𝛾1 𝛽1 𝛾2 𝛽2 𝛾3 𝛽3
We observe that the displacement field being linear, the strains are constant
(elements of [B] are all constants) within an element and thus this element is known as
a Constant Strain Triangle (CST). This may be considered the two-dimensional analogue
of the simple linear truss element, and we expect that it will not be able to accurately
model complex problems involving significant strain/stress gradients.
Element matrices
The size of the stiffness matrix is (6 × 6) and, even after exploiting the
symmetry of the stiffness matrix, we still need to evaluate 21 integrals! Since the
arguments are all constants, the integration can be readily carried out and the element
stiffness matrix obtained as
For more complex, higher order elements, the elements of [B] matrix will vary
from point to point within an element, and explicit integration of expressions in the
above equation becomes very tedious. Also, the size of the matrix itself will be larger
and we need to evaluate many integrals. Thus we look for an efficient way of
computing these integrals numerically within a computer program, which is readily
accomplished using natural coordinates.
External loads can be concentrated forces acting at the nodes along global X
and/or Y axes, and they may be distributed over the element edges as shown in fig.
Two schemes of finding equivalent nodal loads exist. In the first scheme, we
generate a lumped nodal force vector by appropriately lumping the total force. In the
second scheme, we compute the consistent nodal force vector which does the same
amount of work as the distributed force. The element nodal force vector { f }e is a
consistent nodal force vector.
Example . Illustration of the use of CST element. We wish to determine the deflection
of a thin plate subjected to extensional loads as shown in Figure. We will use two CST
elements to solve this problem. From the coordinates of the nodes, we first compute i,
i = 12, i = 0, i = −4
j = 0, j = 3, j = 0
k = 0, k = −3, k = 4
Element 2: Shape function coefficients
i = 12, i = −3, i = 0
j = 0, j = 3, j = − 4
k = 0, k = 0, k = 4
We observe that the area of both the elements is 6 cm2 and we let the thickness
to be unity for plane stress. The strain-displacement matrix [B] is now obtained as given
below.
𝟎. 𝟐𝟓 𝟎 −𝟎. 𝟐𝟓 𝟎 𝟎 𝟎
[𝑩](𝟐) = { 𝟎 𝟎. 𝟑𝟑𝟑 𝟎 𝟎 𝟎 −𝟎. 𝟑𝟑𝟑}
𝟎. 𝟑𝟑𝟑 𝟎. 𝟐𝟓 𝟎 −𝟎. 𝟐𝟓 −𝟎. 𝟑𝟑𝟑 𝟎
Assuming E = 2 × 107 N/cm2 -strain
relation matrix [D] as
The numerical values of [k]e for the two elements are obtained as given below.
1.34
1.75 𝑠𝑦𝑚𝑚𝑒𝑡𝑟𝑖
0.714
[𝑘](1) = 107 −0.824 −0.33 0.824
−0.385 −0.288 0 0.288
−0.513 −0.385 0 0.385 0.513
[−0.330 −1.47 0.33 0 0 1.47]
12 34 56
[𝐿𝑜𝑐𝑎𝑙]
12 56 34
[𝐿𝑜𝑐𝑎𝑙]
1.34
1.75 𝑠𝑦𝑚𝑚𝑒𝑡𝑟𝑖
0.714
[𝑘](2) = 107 −0.824 −0.33 0.824
−0.385 −0.288 0 0.288
−0.513 −0.385 0 0.385 0.513
[−0.330 −1.47 0.33 0 0 1.47]
12 34 56
[𝐿𝑜𝑐𝑎𝑙]
78 34 56
[𝐿𝑜𝑐𝑎𝑙]
Since the two elements are essentially identical, we observe that our stiffness
matrices are also identical.
For ease of assembly, the local and global destinations of the elements of the
stiffness matrices have also been indicated above. Following the standard procedure of
assembly and incorporating the boundary conditions that u1 = 0 = v1 = u2, we get the
final set of equations as follows:
1.75 𝑣2 0
0.714 1.34 𝑠𝑦𝑚𝑚𝑒𝑡𝑟𝑖 𝑢3 0.15
107 0 0 1.75 𝑣3 = 0
−0.385 −0.513 −0.33 1.34 𝑢4 1.5
[−0.288 −0.385 −1.47 0.714 1.75] { 𝑣4 } { 0
}
The RHS force vector is equivalent nodal forces for the distributed force on the
edge 3−4, viz., half the total force at each node in the X-direction. Solving these
equations, we get the nodal deflections as
v2 = −4.5 × 10–8 cm
u3 = 2 × 10–7 cm
v3 = 0
u4 = 2 × 10–7 cm
v4 = −4.5 × 10–8 cm
From the nodal deflections, we can compute element level forces as described
now:
Element 1: Forces
1.34
1.75 𝑠𝑦𝑚𝑚𝑒𝑡𝑟𝑖
0.714
{𝑓}(1) = [𝑘](1){𝛿}(1) = 107 −0.824 −0.33 0.824
−0.385 −0.288 0 0.288
−0.513 0.385 0 0.385 0.513
[−0.330 −1.47 0.33 0 0 1.47]
0 −1.5
0 0
2 × 10 − 7 1.5
× 0 =
0
0 0
{ −4.5 × 10−8 } { 0 }
Element 2: Forces
1.34
1.75 𝑠𝑦𝑚𝑚𝑒𝑡𝑟𝑖
0.714
{𝑓}(2) = [𝑘](2){𝛿}(2) = 107 −0.824 −0.33 0.824
−0.385 −0.288 0 0.288
−0.513 0.385 0 0.385 0.513
[−0.330 −1.47 0.33 0 0 1.47]
2 × 10 − 7 1.5
−4.5 × 10−8 0
0 −1.5
× =
−4.5 × 10−8 0
2 × 10 − 7 0
{ −4.5 × 10−8 } { 0 }
These forces are along the d.o.f. From the nodal deflections, the strains in
elements 1 and 2 can be computed as follows:
Element 1: Strains
𝜀𝑥 −0.25 0 0.25 0 0 0
[𝜀 ](1) 𝜀
={ 𝑦}= { 0 −0.333 0 0 0 0.333}
𝛾 −0.333 −0.25 0 0.25 0.333 0
𝑦
00 5 × 10−8
2 × 10 − 7 = {
× −1.5 × 10−8 }
0
0
0
{ −4.5 × 10−8 }
Element 2: Strains
𝜺𝒙
[𝜺](𝟐) = { 𝜺𝒚 }
𝜸𝒙𝒚
𝟎. 𝟐𝟓 𝟎 −𝟎. 𝟐𝟓 𝟎 𝟎 𝟎
={ 𝟎 𝟎. 𝟑𝟑𝟑 𝟎 𝟎 𝟎 −𝟎. 𝟑𝟑𝟑}
−𝟎. 𝟑𝟑𝟑 𝟎. 𝟐𝟓 𝟎 −𝟎. 𝟐𝟓 −𝟎. 𝟑𝟑𝟑 𝟎
𝟐 × 𝟏𝟎 − 𝟕
−𝟒. 𝟓 × 𝟏𝟎−𝟖
𝟓 × 𝟏𝟎−𝟖
𝟎
× −𝟒. 𝟓 × 𝟏𝟎−𝟖 = {−𝟏. 𝟓 × 𝟏𝟎−𝟖 }
𝟐 × 𝟏𝟎 − 𝟕 𝟎
{ −𝟒. 𝟓 × 𝟏𝟎 }
−𝟖
It is observed that our solution tallies perfectly with the expected solution,
viz.,
= (PL/AE) and x and y are in the ratio of v. The element stresses can be readily
obtained from these strains.
While we can get excellent results for this simple problem with the CST
element, in general, we cannot expect such a good performance because it is only a
constant strain element. In problems involving significant gradients of stress (e.g.,
stress concentration near a hole or any sudden geometry change), the CST element
will not give good results.
Four-noded Rectangular Element
It has four nodes, straight sides, and each node admits two d.o.f., viz., u and
v translations along global the X and Y axes. Our formulation of the element follows
standard procedure and is described now.
𝑥 𝑦 𝑥𝑦
𝑁1 = 1 − − +
ℓ ℎ ℓℎ
𝑥 𝑥𝑦
𝑁2 = −
ℓ ℓℎ
𝑥𝑦
𝑁3 =
ℓℎ
𝑦 𝑥𝑦
𝑁4 = −
ℎ ℓℎ
𝑢1
𝜕𝑢 𝜕𝑁 𝜕𝑁 𝜕𝑁3 𝜕𝑁 𝑣1
0 0 0
𝜕𝑥 1 2𝜕𝑥 𝜕𝑥 4𝜕𝑥 𝑢
𝜀𝑥 𝜕𝑥 𝜕𝑁4
𝜕𝑣 𝜕𝑁1 𝜕𝑁2 2 𝑣2
{𝜀} { 𝜀𝑦 } =
𝑒 = 0 0 0 0 𝜕𝑦
𝛾𝑥 𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝑢3
𝜕𝑁4 𝜕𝑁4 𝑣3
𝜕𝑢 𝜕𝑢 𝜕𝑁1 𝜕𝑁 2 𝜕𝑁 2 𝜕𝑁3
+ 𝜕𝑥
{ 𝜕𝑦 𝜕𝑥} 𝜕𝑁 1 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑦 ]3−8 𝑢4
[ 𝜕𝑦 { 𝑣4 } 8×1
Thus the strain-displacement relation matrix [ B] is given by
𝑦 1 1 𝑦 𝑦 −𝑦
( − ) 0 ( − ) 0 ( ) 0 ( ) 0
ℓℎ ℓ ℓ ℓℎ ℓℎ ℓℎ
𝑥 1 −𝑥 𝑥 1 𝑦
[𝐵] = 0 ( − 0 ( ) 0 ( ) 0 ( )
) ℓℎ ℓℎ ℎ ℓℎ
𝑥 1 𝑦 1 −𝑥 1 𝑦 𝑥 𝑦 1 𝑦 − −𝑦
[( ℓℎ − ℓ ) ( ℓℎ − ℓ ) (
ℓℎ
) ( −
ℓ ℓℎ
) ( )
ℓℎ
( )
ℓℎ
( − )
ℎ ℓℎ
( ) ]
ℓℎ
We observe that on x = constant line, the strain varies linearly with y and vice
versa. This element is known as a bilinear rectangular element since its shape
functions involve linear polynomials in both the x and y directions. The stress-strain
relation matrix [D] is given in Eq. (5.91) for plane stress and in Eq. (5.92) for plane
strain.
[k]e = ∫[B]T[D][B] dv
The size of the stiffness matrix is (8 × 8) since there are four nodes and two
d.o.f. per node and even after exploiting the symmetry of the stiffness matrix, we
still need to evaluate 36 integrals! Also, the elements of [ B] matrix now vary from
point to point within an element and explicit integration of expressions in the above
equation becomes very tedious.
External loads can be concentrated forces acting at the nodes along the
global X and/or Y axes (or) they may be distributed over the element edges as
shown in Figure.
Two schemes of finding equivalent nodal loads exist. In the first scheme, we
generate a lumped nodal force vector by appropriately lumping the total force. In
the second scheme, we compute the consistent nodal force vector which does the
same amount of work as the distributed force. The element nodal force vector { f}e is
a consistent nodal force vector.
Example:- 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.
37
βi = y j - yk
β1 = y2 - y3 = 0 − 1 = −1
β2 = y 3 - y1 = 1 + 1 = 2
β3 = y1 - y 2 = −1 − 0 = −1
γi = − (x j - xk )
γ1 = − (x2 - x3 ) = 0 − 2 = −2
γ2 = − (x3 - x1 ) = 0 − 0 = 0
γ3 = − (x1 - x2 ) = 2 − 1 = 2
. A=½xbxh=½x2x2=2
u1
v
β1 0 β1 0 β1 0 1
1
u 2
ε = 0 λ1 0 λ1 0 λ1
2Ae v2
λ1 β1 λ1 β1 λ1 β1
u3
v3
= B d
−1 0 2 0 −1 0
1
B = 0 − 2 0 0 0 2
2*2
− 2 − 1 0 2 2 − 1
1 0
D = E 2 1 0
1− 1−
0 0
2
1 0.3 0
2x105
D = 0.3 1 0
1 − 0.3
2
1 − 0.3
0 0
2
− 1 0 − 2
0 − 2 − 1 39
1 0.3
2x2x105 2 0
0 1 − 1 0 2 0 − 1 0
0 0 − 2 0 0 0 2
0
K = 0.3 1
4x0.91 0 0 2 1 − 0.3 4
− 1 0 2
0 0 − 2 − 1 0 2 2 − 1
2
0 2 − 1
6x6
Note: In order to evaluate the element stress we can use the equation
{} = [D] [B] {d}
α 2 = x3 y1 - x1 y 3 = 2*0 − 0*5 = 0
α 3 = x1 y 2 - x2 y1 = 0*0.5 − 4*0 = 0
β1 = y2 - y3 = 0.5 − 5 = −4.5
β2 = y 3 - y1 = 5 − 0 = 5
β3 = y1 - y2 = 0 − 0.5 = −0.5
γ1 = − (x 2 - x3 ) = −(4 − 2) = −2
γ2 = − (x 3 - x1 ) = −(2 − 0) = −2
γ3 = − (x 1 - x2 ) = −(0 − 4) = 4
1 x1 y1 1 0 0
2Ae = 1 x2 y 2 = 1 4 0.5 = 19cm 2
1 x3 y3 1 4 5
1 1
N 1 (x,y) = ( 1 + β1 x + γ1 y) = (19 − 4.5x − 2y)
2Ae 19
1 1
N 2 (x,y) = (α 2 + β 2 x + γ 2 y) = (5x − 2y)
2Ae 19
1 1
N 3 (x,y) = (α3 + β3 x + γ3 y) = ( − 0.5x + 4y)
2Ae 19
P (x,y) = N 1 P1 + N 2 P2 + N 3 P3
=
1
(19 − 4.5x − 2y)40 + (5x − 2y)34 + ( − 0.5x + 4y)46
19
Nodal Coordinates
42
The length of the sides are
2b = x2 – x1 = 8 – 5 = 3
2a = y4 – y1 = 5 – 3 = 2
Substituting these values in the shape
functions
x y xy
N 1 = 1 − 1 − N3 =
3 2 6
x y y x
N 2 = 1 − N 4 = 1 −
3 2 2 3
Inspection reveals that the 50°C contour line
intersects the sides 3-4 and 1-2; therefore, we
need to assume values of y and calculate
values of x. Along side 1-2, y=0 and
x x
T ( x, y ) = 1 − T1 + T2 = 50
3 3
Substituting for T1 and T2 and solving gives
x=2.0. Along side 4-2, y= 2b =3 and
x y
T ( x, y ) = T4 + 1 − T3 = 50
3 3
Substituting for T4 and T3 and solving gives x = 1.2
x
(− 42 + 54 + 56 − 46) + 1 (42 + 46) = 50
6 2
01 Calculate the element stiffness matrix and the thermal force K2 CO4
vector for the axisymmetric triangular element as shown
below. The element experiences a 150º C increase in
temperature. Take = 10 x 10-6 / ºC, E = 2 x 105 N/mm2
and µ= 0.25
02 Calculate the element stiffness matrix and the thermal force K2 CO4
vector for the axisymmetric triangular element as shown
below. The element experiences a 100º C increase in
temperature. Take = 15 x 10-6 / ºC, E = 2.1 x 105
N/mm2 and µ= 0.3
05 K2 CO4
A bilinear rectangular element has coordinates as shown in
figure and the nodal temperatures are T1 = 100ºC, T2=
60ºC, T3 = 50ºC and T4 = 90ºC. Compute the temperature
at the point whose coordinates are (2.5, 2.5). Also Compute
the 80ºC isotherm.
PART - A QUESTIONS AND ANSWERS
UNIT – II INTRODUCTION
S.N O QUESTION AND ANSWER K CO
LEVEL
1 How do you define two dimensional elements? K1 CO4
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
element.
47
PART - A QUESTIONS AND ANSWERS
UNIT – I INTRODUCTION
S.N O QUESTION AND ANSWER K CO
LEVEL
8 Write down the stress strain relationship matrix K1 CO4
for plane stress conditions.
1− 0
E
1+ 1−2 0 1 − 0
1−2
0 0
2
▪ It is a symmetric matrix
▪ The sum of elements in any column must be equal
to zero
▪ It is an unstable element. So, the determinant is
equal to zero.
18 What do you mean by the terms: C0, C1 and Cn K1 CO4
continuity?
▪ C0 – Governing differential equation is
quasiharmonic, ø has to be continuous.
▪ C1 – Governing differential equation is biharmonic,
øas well as derivative has to be continuous inside
and between the elements.
▪ Cn – Governing differential equations is polynomial.
04 Calculate the element stresses, for the plane strain element K2 CO4
shown in the figure. Take the nodal displacements as : u1=
0.005 mm, u2 = 0.002 mm, u3=0.0mm, u4 = 0.0 mm, u5 =
0.004 mm, u6 = 0.0 mm.. Take E = 200 Gpa and = 0.3.
Use unit thickness for plane strain.
40
S.No. Question K Level CO
05 Calculate the deflection at the point load and stress induced K2 CO4
for the configuration as shown in figure. Assume plane
stress condition
10 K2 CO4
Calculate the pressure at the location (7, 4) in a rectangular
plate with the data shown in Figure and also draw 50 MPa
contour line
53
REAL TIME APPLICATION IN DAY TO DAY LIFE AND TO
INDUSTRY
54
CONTENT BEYOND SYLLABUS
55
ASSESSMENT SCHEDULE
56
PRESCRIBED TEXT BOOKS & REFERENCE BOOKS
TEXT BOOKS:
T1) Reddy. J.N., “An Introduction to the Finite Element Method”, 3rd Edition, Tata
McGraw Hill, 2005
T2) Seshu, P, “Text Book of Finite Element Analysis”, Prentice-Hall of India Pvt. Ltd.,
New Delhi, 2007.
REFERENCES:
R1) Rao, S.S., “The Finite Element Method in Engineering”, 3rd Edition, Butterworth
Heinemann, 2004
R2) Logan, D.L., “A first course in Finite Element Method”, Thomson Asia Pvt. Ltd.,
2002
R3) Robert D. Cook, David S. Malkus, Michael E. Plesha, Robert J. Witt, “Concepts
and Applications of Finite Element Analysis”, 4th Edition, Wiley Student Edition,
2002.
57
MINI PROJECT SUGGESTIONS
The transfer of a model from PTC Creo to ANSYS will be demonstrated here for a
simple solid model. Model idealizations such as shells and plane will be treated to
find displacement and stress distributions.
58
Thank you
Disclaimer:
This document is confidential and intended solely for the educational purpose of RMK Group of
Educational Institutions. If you have received this document through email in error, please notify the
system manager. This document contains proprietary information and is intended only to the
respective group / learning community as intended. If you are not the addressee you should not
disseminate, distribute or copy through e-mail. Please notify the sender immediately by e-mail if you
have received this document by mistake and delete this document from your system. If you are not
the intended recipient you are notified that disclosing, copying, distributing or taking any action in
reliance on the contents of this information is strictly prohibited.