ME8692 QB Finite Element Analysis
ME8692 QB Finite Element Analysis
ME8692 QB Finite Element Analysis
net
UNIT I: INTRODUCTION
et
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
.n
Technique – Basic concepts of the Finite Element Method.
PART A
1 Illustrate the methods generally associated with the finite element
pz
BT2 Understanding
analysis.
2 Write the stiffness matrix for a one dimensional 2 noded linear element. BT3 Applying
3 Point out any four advantages of finite element method BT3 Applying
4 Why polynomial type interpolation functions are mostly used in FEM. BT4 Analyzing
5 State the advantages of Rayleigh Ritz method. BT2 Understanding
ee
6 Compare the Ritz technique with the nodal approximation method. BT4 Analyzing
7 How to develop the equilibrium equation for a finite element? BT6 Creating
8 What do you mean by constitutive law? BT2 Understanding
9 List the various method of solving boundary value problems. BT1 Remembering
10 Formulate the boundary conditions of a cantilever beam AB of span L
ad
19 List out the advantages of finite element method over other numerical
BT1 Remembering
analysis method.
20 List the various weighted residual methods BT1 Remembering
w
PART B
1 i. Explain the various methods of engineering analysis with suitable
illustrations. (8) BT1 Remembering
ii. Describe the principle of stationary total potential energy. (5)
2 Using collocation method, find the solution of given governing equation BT3 Applying
www.padeepz.net
www.padeepz.net
𝑑2 Φ
+ Φ + X = 0, 0 ≤ X ≤ 1 subject to the boundary conditionsΦ(0) =
𝑑𝑋 2
Φ(1) = 0. Use X=1⁄4 and 1⁄2 as the collocation points.
3 A uniform rod subjected to a uniform axial load is illustrated in figure,
the deformation of the bar is governed by the differential equation given
below. Determine the displacement by applying Weighted Residual
Method (WRM).
et
BT3 Applying
.n
4 Describe the step by step procedure of solving FEA. BT4 Analyzing
pz
5 A tapered bar made of steel is suspended vertically with the larger end
rigidly clamped and the smaller end acted on by a pull of 10 5N. The
areas at the larger and smaller ends are 80 cm 2 and 20 cm2 respectively.
The length of the bar is 3m. The bar weighs 0.075 N/cc. Young’s BT1 Remembering
modulus of the bar material is E=2 x 107 N/cm2. Obtain an approximate
ee
expression for the deformation of the rod using Ritz technique.
Determine the maximum displacement at the tip of the bar.
6 Using any one of the weighted residue method, find the displacement of .BT3 Applying
𝑑 𝑑𝑢 2
given governing equation 𝑑𝑥 [𝑥 𝑑𝑥 ] - 𝑥2 = 0, 1 < 𝑥 < 2 at x=1 , u=2 at
𝑑𝑢 1 BT2 Understanding
ad
x=2 , x𝑑𝑥 = − 2
7 The Governing Equation for one dimensional heat transfer through a fin of
length l attached to a hot source as shown in fig is given by
.p
BT3 Applying
w
If the free end of the fin is insulated, give the boundary conditions and
determine using the Collocation technique the temperature distribution in the
fin. Report the temperature at the free end.
w
www.padeepz.net
www.padeepz.net
d 2 y dy
(ii) Find a solution of the initial value problem 2y 0 ,
dx 2 dx BT6 Creating
boundary conditions y(0) = 2, 𝑦 ′(0)=5.
et
boundary conditions as y (0) = 0 and y (1) = 0 using the Galerkin’s BT3 Applying
method with the trial function No(x) = 0; N1(x) = x (1-x2).
13 Solve the differential equation for a physical problem expressed as
.n
d2y/dx2 + 50 = 0, 0≤x≤10 with boundary conditions as y (0) = 0 and y
(10) = 0 using the trial function y = a1x (10-x) find the value of the
BT3 Applying
parameters a1 by the following methods listed below.
(i) Point collocation method (ii) Sub domain collocation method
pz
(iii) Least squares method and (iv) Galerkin method
14 Find the solution of the boundary value problem 𝑦 ′′ + 4𝑦 = 0 with
𝜋 𝜋 BT6 Creating
y( )=0, y( )=1.
8 6
BT6 Creating
.p
BT5 Evaluating
w
www.padeepz.net
www.padeepz.net
et
. UNIT 2 ONE-DIMENSIONAL PROBLEMS
One Dimensional Second Order Equations – Discretization – Element types- Linear and Higher order
Elements – Derivation of Shape functions and Stiffness matrices and force vectors- Assembly of
.n
Matrices - Solution of problems from solid mechanics and heat transfer. Longitudinal vibration
frequencies and mode shapes. Fourth Order Beam Equation –Transverse deflections and Natural
frequencies of beams
PART A
pz
1 Write the conduction, convection and thermal load matrices for 1D heat BT1
Remembering
transfer through a fin.
2 Define shape function. BT1 Remembering
3 Illustrate shape function of a two node line element BT3 Applying
4
ee
List out the stiffness matrix properties. BT1 Remembering
5 Mention the characteristics of shape functions BT2 Understanding
6 Differentiate global and local coordinate. BT2 Understanding
7 Express the element stiffness matrix of a truss element BT2 Understanding
8 Illustrate a typical truss element shown local global transformation BT3 Applying
ad
14 List out the expression of governing equation for free axial vibration of
BT1 Remembering
rod and transverse vibration of beam.
15 Determine the element mass matrix for one-dimensional dynamic
structural analysis problems. Assume the two-node, linear element. BT4 Analyzing
w
20 Illustrate the expression of longitudinal vibration of the bar element. BT3 Applying
PART B
1 Formulate the shape function for One-Dimensional Quadratic bar
BT6 Creating
element.
2 A steel bar of length 800mm is subjected to an axial load of 3kN as BT2 Understanding
www.padeepz.net
www.padeepz.net
shown in fig. Estimate the nodal displacement of the bar and load
vectors. Take E= 2 x 105 N/mm2.
3 For the bar element as shown in the figure. Calculate the nodal
et
displacements and elemental stresses. Take E= 2.1 x 105 N/mm2.
.n
BT3 Applying
pz
4 Determine the Eigen values for the stepped bar shown in figure
BT5 Evaluating
ee
5 Consider a bar as shown in figure an axial load of 200kN is applied at a
point P. Take A1=2400 mm2 , E1=70x109 N/mm2 A2=600 mm2 and E2 =
200x109 N/mm2 . Calculate the following (i) the nodal displacement at
ad
BT4 Analyzing
.p
w
model. The bar has a mass density ρ = 7800 kg/m3, the young’s
BT2 Understanding
modulus E = 2x105 MN/m2. In addition to self-weight, the bar is
subjected to a point load P= 1 kN at its Centre. Also determine the
w
www.padeepz.net
www.padeepz.net
et
.n
7 Consider a bar as shown in fig Young’s Modulus E=2 x 105 N/mm2 A1=
2 cm2;A2 = 1cm2 and force of 100N. Calculate the nodal displacement.
pz
BT3 Applying
ee
8 A metallic fin 20 mm wide and 4 mm thick is attached to a furnace
whose wall temperature is 180 °C. The length of the fin is 120 mm. if
the thermal conductivity of the material of the fin is 350 W/m °C and
BT5 Evaluating
ad
BT2 Understanding
w
w
BT4 Analyzing
dimensions are in mm E=200GPa. and ρ= 0.78 kg/cc. A=4 cm 2, length
l=500mm.
www.padeepz.net
www.padeepz.net
et
cantilever beam shown in figure by applying one 1D beam element.
BT1 Remembering
.n
pz
12 A two noded truss element is shown in figure. The nodal displacements
are u1=5mm and u2= 8mm. calculate the displacement at x=L/4, L/3 and
L/2
ee
BT3 Applying
ad
13 For the two bar truss shown in the fig, Estimate the displacements of
node 1 and the stress in element 1-3.Take E=70GPa A=200 mm2
.p
BT4 Analyzing
w
w
14 Determine the deflection in the beam loaded as shown in fig. at the mid-
span and at the length of 0.5m from the left support. Determine also the BT3 Applying
reactions at the fixed ends. E= 200 GPa. I1= 20 x 10-6 m4 I2 = 10 x 10-6 m4.
w
www.padeepz.net
www.padeepz.net
et
PART C
1 Develop the Shape function, Stiffness matrix and force vector for
BT6 Creating
one dimensional linear element.
.n
2 Consider the bar shown in figure axial force P = 30kN is applied as
shown. Determine the nodal displacement, stresses in each element and
reaction forces
pz
ee BT5 Evaluating
3 For the beam and loading as shown in figure. Calculate the slopes at
nodes 2 and 3 and the vertical deflection at the mid-point of the
ad
BT3 Applying
.p
w
4 Calculate the force in the members of the truss as shown in fig. Take
E=200 GPa
w
BT3 Applying
w
www.padeepz.net
www.padeepz.net
UNIT 3 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.
PART A
1 Show the displacement function equation for CST element. BT3 Applying
et
2 How will you modify a three-dimensional problem to a Two-
dimensional problem? BT6 Creating
3 List out the application of two-dimensional problems. BT1 Remembering
4 What is meant by steady state heat transfer? Write down its governing
.n
differential equation. BT2 Understanding
pz
7 Define QST (Quadratic strain Triangle) element. BT4 Analyzing
8 Relate path line with streamline. BT3 Applying
9 Formulate the (B) matrix for CST element. BT6 Creating
10 Express the interpolation function of a field variable for three-node
BT2 Understanding
triangular element
ee
11 List out the CST and LST elements. BT1 Remembering
12 Illustrate the shape function of a CST element. BT3 Applying
13 Distinguish between scalar and vector variable problems in 2D. BT1 Remembering
14 Express the nodal displacement equation for a two dimensional
ad
20 Write the Lagrange shape functions for a 1D, 2noded elements. BT5 Evaluating
w
PART B
1
Derive the conductance matrix for a 3 noded triangular element whose
nodal coordinates are known. The element is to be used for two BT6 Creating
w
www.padeepz.net
www.padeepz.net
et
4 Calculate the element equations for the element shown in fig. which
.n
experiences convection on the side jk and its upper face.
pz
BT4 Analyzing
ee
5 Find the temperature distribution in a square region with uniform
energy generation as shown in fig. assume that there is no temperature
ad
BT5 Evaluating
.p
w
6 Derive the expression for shape function for heat transfer in 2D element. BT2 Understanding
w
www.padeepz.net
www.padeepz.net
8 For a 4-noded rectangular element shown in fig. Infer the temperature at
the point (2.5, 2.5). The nodal values of the temperatures are T1= 100°C,
T2= 60°C and T3= 50°C and T4= 90°C. Also determine the 80°C
isotherm.
BT4 Analyzing
et
9 Calculate the element stiffness matrix and temperature force vector for
the plane stress element shown in fig. The element experiences a 20ºC
increase in temperature. Assume α= 6x10-6 C. Take E=2x105 N/mm2,
.n
v= 0.25, t= 5mm.
BT3 Applying
pz
ee
10 For a 4-noded rectangular element shown in fig. Estimate the
temperature point (7,4). The nodal values of the temperatures are T 1 =
42˚C, T2 = 54˚C and T3 = 56˚Cand T4 = 46˚C. Also determine 3 point
on the 50˚C contour line.
BT2 Understanding
ad
11 The figure below shows a shaft having rectangular cross section with 8
.p
BT5 Evaluating
w
w
www.padeepz.net
www.padeepz.net
et
coordinates of the nodes are (0,-4), (8,0) and (0,4) at nodes i, j, k. BT5 Evaluating
Assume plane stress condition E= 200 GPa, Poisson’s ratio = 0.35
.n
PART C
1 The x,y coordinates of nodes i,j and k of a triangular element are given
by (0,0) (3,0) and (1.5,4) mm respectively. Evaluate the shape functions
N1,N2 and N3 at an interior point P(2,2.5)mm of the element. Evaluate
pz
BT5 Evaluating
the strain displacement relation matrix B for the above same triangular
element and explain how stiffness matrix is obtained assuming scalar
variable problem.
2 For the triangular element shown in the figure determine the strain-
displacement matrix [B] and constitutive matrix [D]. Assume plane
ee
stress conditions. Take µ=0.3 , E=30 x 106 N/m2 and thickness t= 0.1 m.
And also calculate the element stiffness matrix for the triangular
element
BT4 Analyzing
ad
.p
BT4 Analyzing
w
w
www.padeepz.net
www.padeepz.net
element matrices are given below. Carry out the assembly and solve for
the unknown stress function values.
1 −1 0 29.1
For triangular K = 1⁄2 [−1 2 −1] r = [29.1]
0 −1 1 29.1
4 1 2 1 43.6
1 4 1 2 43.6
1
For rectangle K = ⁄6 r=
2 1 4 1 43.6
et
1 2 1 4 43.6
.n
pz
ee
UNIT 4 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
PART A
ad
www.padeepz.net
www.padeepz.net
the following structures. Explain your answer.
A wall subjected to wind load
A wrench subjected to a force in the plane of the wrench.
19 Define a plane strain with suitable example. BT1 Remembering
20 Define a plane stress problem with a suitable example. BT1 Remembering
PART B
1 Determine the stiffness matrix for the axisymmetric element shown in
figure. Take E=2.1 x 105 N/mm2, v=0.25 The coordinates are in mm.
et
BT3 Applying
.n
pz
2 Develop shape function for axisymmetric triangular elements BT6 Creating
3 Triangular element are used for the stress analysis of plate subjected to
inplane loads. The (x,y) coordinates of nodes i, j, and k of an element
are given by (2,3), (4,1), and (4,5) mm respectively. The nodal
ee
displacement are given as : BT1 Remembering
u1=2.0 mm, u2=0.5 mm, u3= 3.0 mm v1=1.0 mm, v2= 0.0 mm, v3= 0.5
mm Examine element stress. Let E=160GPa, poisson's ratio = 0.25 and
thickness of the element t=10 mm
4 The nodal coordinates for an axisymmetric triangular element are given
ad
BT5 Evaluating
.p
BT1 Remembering
w
www.padeepz.net
www.padeepz.net
7 The nodal coordinates for an axisymmetric triangular element shown in
fig are given below. Evaluate the strain-displacement matrix for that
element
BT4 Analyzing
et
8 Calculate the element strains for an axisymmetric triangular element
shown in fig the nodal displacement are. u 1= 0.001, u2 = 0.002, u3 = -
0.003, w1 = 0.002, w2 = 0.001 and w3 = 0.004 all dimensions are in mm.
.n
BT3 Applying
pz
9 Estimate the global stiffness matrix for the plate shown in fig. Taking
two triangular elements. Assume plane stress conditions
ee
BT2 Understanding
ad
10 Explain the classification of the shell elements and also brief the
BT4 Analyzing
assumptions used in Finite element Analysis of Shell element.
11 A tin plate of thickness 5mm is subjected to an axial loading as shown BT4 Analyzing
.p
BT1 Remembering
w
www.padeepz.net
www.padeepz.net
12 Evaluate the Stress-Strain relationship matrix for axisymmetric
BT5 Evaluating
triangular element
13 Develop Strain-Displacement matrix for axisymmetric triangular
BT6 Creating
element
14 Derive the Finite element equation for triangular plate bending element
BT4 Analyzing
with 9 degrees of freedom.
PART C
1 Develop the four basic sets of elasticity equation BT6 Creating
2 A long hollow cylinder of inside diameter 100mm and outside diameter
et
120mm is firmly fitted in a hole of another rigid cylinder over its full
length as shown in fig. The cylinder is then subjected to an internal
pressure of 2 MPa. By using two element on the 10mm length shown
calculate the displacements at the inner radius tame E = 210 GPa. μ =
.n
0.3
BT3 Applying
pz
3
ee
Triangular element are used for the stress analysis of plate subjected to
in plane loads. The (x,y) coordinates of nodes 1, 2, and 3 of an element
are given by (5,5), (25,5), and (15,15) mm respectively. The nodal
BT5 Evaluating
displacement are given as : u1=0.005 mm, u2=0.002 mm, u3= 0.0 mm,
u4=0.0 mm, u5= 0.005 mm, u6= 0.0 mm. Evaluate element stress. Let
ad
E= 200 GPa, poisson's ratio = 0.3 and use unit thickness of the element
4 For an axisymmetric triangular elements as shown in fig. Evaluate the
stiffness matrix. Take modulus of elasticity E = 210 GPa. Poisson’s
ratio = 0.25. the coordinates are given in millimeters
BT5 Evaluating
.p
w
Equations of elasticity – Plane stress, plane strain and axisymmetric problems – Body forces and temperature
effects – Stress calculations - Plate and shell elements.
PART A
w
www.padeepz.net
www.padeepz.net
element.
7 List out the advantages of Gauss quadrature numerical integration for
BT1 Remembering
Isoparametric element
8 Define Isoparametric element BT2 Understanding
9 Discuss about Numerical integration BT2 Understanding
10 Discuss about Gauss-quadrature method. BT4 Analyzing
11 Differentiate between implicitly and explicitly methods of numerical Analyzing
BT4
integration
12 Differentiate between geometric and material non-linearity. BT1 Remembering
13 List out the significance of Jacobian transformation BT1 Remembering
et
14 Define Isoparametric element with suitable examples. BT6 Creating
15 Develop Stress- displacement matrix for Four noded quadrilateral Creating
BT6
element using natural coordinates.
.n
16 Develop Stiffness matrix for Isoparametric quadrilateral element BT1 Remembering
17 Define Newton cotes quadrature method BT2 Understanding
18 Distinguish between trapezoidal rule and Simpson’s rule BT2 Understanding
19 Distinguish between trapezoidal rule and Gauss quadrature. BT5 Evaluating
pz
20 State the transformation for mapping x-coordinate system into a natural
coordinate system for a linear spar element and for a quadratic spar BT4 Analyzing
element ee PART B
1
For the four noded element shown in fig. determine the Jacobian
and evaluate its value at the point (0,0)
ad
BT1 Remembering
.p
BT3 Applying
w
3 Develop the shape function for 4 noded isoparametric quadrilateral element BT6 Creating
4 Develop the strain displacement matrix, stress-strain matrix and
BT6 Creating
stiffness matrix for an isoparametric quadrilateral element
5 Evaluate the Jacobian matrix at the local coordinates ε=η= 0.5 for
BT5 Evaluating
the linear quadrilateral element with its global coordinates as shown in
www.padeepz.net
www.padeepz.net
fig. Also evaluate the strain-displacement matrix
et
6 For the four noded quadrilateral element shown in fig analysis the
Jacobian and evaluate its value at the point (1/2, 1/2)
.n
BT4 Analyzing
pz
Calculate the Cartesian coordinates of the point P which has local
ee
coordinates ε = 0.8 and η = 0.6 as shown in figure
BT3 Applying
ad
x
8 1
Evaluate 4
x 2 dx by applying 3 point Gaussian quadrature. BT5 Evaluating
1
.p
9 1
Evaluate e x dx by applying 3 point Gaussian quadrature. BT5 Evaluating
1
10 1
x
Evaluate the integral, I = Cos 2 dx by applying 3 point Gaussian
w
1
BT4 Analyzing
quadrature and compare with exact solution.
11 For a four noded rectangular element shown in fig. Estimate the
w
following
a. Jacobian matrix b. Strain-Displacement matrix c. Element strain BT2 Understanding
and d. Element stress
w
www.padeepz.net
www.padeepz.net
12 For the element shown in the figure. Calculate the Jacobian matrix.
et
.n
BT3 Applying
pz
13 Consider the isoparametric quadrilateral element with nodes 1 to 4 at
(5,5), (11,7),(12,15), and (4,10) respectively. Estimate the jacobian BT2 Understanding
matrix and its determinant at the element centroid
ee
14 Tabulate the element characteristics of a four node quadrilateral element BT1 Remembering
ii) Using Gauss quadrature evaluate the following integral using 1,2 and BT6 Creating
3 point integration.
𝟏 𝐒𝐢𝐧 𝐒
∫−𝟏 𝐒 (𝟏−𝐒𝟐 ) 𝐝𝐬
3 For the four noded element shown in Fig,
w
(i) determine the Jacobian and evaluate its value at the point (1⁄3, 1⁄3)
(ii) Using energy approach derive the stiffness matrix for a 1D linear
isoparametric element.
w
BT5 Evaluating
w
www.padeepz.net
www.padeepz.net
4 For the isoparametric quadrilateral element shown in figure, the
Cartesian coordinates of point ‘P’, are (6,4). The loads 10 kN and 12 kN
are acting in x and y
direction on that point P.
Evaluate the nodal
forces.
BT4 Analyzing
et
.n
pz
ee
ad
.p
w
w
w
www.padeepz.net