St7201 Finite Element Method
St7201 Finite Element Method
St7201 Finite Element Method
YEAR/SEM: I/II
PART-A
UNIT-I
1. What is the basic concept of finite element analysis?
2. What is meant by discretization?
3. What is the requirement of displacement field to be satisfied in the use of Rayleigh Ritz
method?
4. Give examples of Eigen value problems in structural mechanics.
5. How will you classify essential and non-essential boundary condition?
6. Derive thermal load vector for 1D bar element.
7. Explain body force and surface force with examples.
8. What do you mean by convergence in finite element analysis?
9. What is the significance of weak formulation?
10. What is quadratic shape function?
11.Why polynomial shape functions are preferred?
12. What are the limitations of Galerkin formulation?
13. Write down the stiffness matrix for 2D beam element.
14. What is shape function?
15. What is displacement and shape function?
16. What is Ritz technique?
17. Write the potential energy for beam of span L simply supported at ends, subjected to a
concentrated P at midspan. Assume EI constant.
18. What are the advantages of FEA?
19. Define banded structure?
20. Define the modulus of resilience?
PART-B
1. Solve the following equations by Gauss elimination method?
i)
ii)
iii)
iv)
28r1+6r2 =1
6r1+ 24r2+6r3=0
6r2+28r3+8r4= -1
8r3+16r4 = 10
2. A simply supported beam is subjected to uniformly distributed load over entire span.
Determine the bending moment and deflection at the mid span using Rayleigh Ritz method
and compare with exact solution. Use a two term trial function y=a1sin(x/l)+a2sin(3x/l).
3. Discuss Rayleigh Ritz and Galerkin methods of formation by taking an example.
4. Show from first principle that the stiffness matrix of a general finite element can be
evaluated in the form of K= BTCBdr
V
5. Derive stiffness matrix for a one dimensional bar subjected to both distributed load and
point loads. (Note: Differential equation should be formed and then apply basic Galerkin
method on differential equation to form element stiffness matrix.)
6. Solve the following differential equation using Ritz method.
d2y/dx2 = -sin (x) boundary conditions u(0) = 0 and u(1) = 0.
7.i) Taking a differential equation , explain the process of weak formulation.
ii) Explain any two methods of weighted residuals with examples.
8. i)Derive the element stiffness matrix and element force matrices for a one dimensional line
element.
ii) A simply supported beam of span L, youngs modulus, moment of inertia I is subjected to
a uniformly distributed load of P/unit length. Determine the deflection W at the midspan. Use
Rayleigh Ritz method.
9.For the bar shown in fig, determine the nodal displacement, stress in each material and
reaction forces.
1. Aluminum
A1= 2400 mm2
E1= 70GPa
2.Steel
A2= 600mm2
E2= 200GPa
P= 200kN
10kN/m
3m
3m
A
UNIT II
PART-B
1. Derive shape functions and stiffness matrix for a 2D rectangular element.
2. Evaluate the nodal load vector due to self-weight of a four noded rectangular element with
two degrees of freedom (translations) at each node. Use Gauss quadrature method of
numerical integration.
3. Evaluate the shape functions N1, N2,N3 at the interior point P(3.85,4.8) for the triangular
element shown in fig 1.
y
3(4,7)
2(7,3.5)
1(1.5,2)
x
4. Discuss the convergence requirements of interpolation polynomials.
5. The nodal coordinates of a quadrilateral element as given as (0,7), (9,4),(7,9) and (2,8).
Evaluate the integral (x2 +y2 -3xy) dA over the area of the element of second order
Gaussian quadrature.
6. Derive shape functions for constant strain triangle element.
7. Integrate the following function using Gaussian integration. Proper Gauss points should be
specified. The x limit is varying from 0 to 2 and y limit is varying from 1 to 3 (xy) dxdy.
8. i) Derive the weights and Gauss points of Gauss one point formula and two point formula.
ii) Derive the shape functions of a four noded quadrilateral element.
11
3(40,40)
1(20,20)
2(40,20)
x
UNIT-III
PART-B
1. Determine the stiffness for the axisymmetric element shown in fig. Take E as 2.1x 105
50,0
0,0
r =-1
r=0
r=+1
x,u
X1
L/2
L/2
4. State the need for mesh refinement. Discuss the methods of mesh refinement.
5. Derive the shape functions for element shown in fig. Shape functions should be specified
in natural coordinate system..
7
555555556
6. Derive the shape functions for ID cubic element. Shape functions should specified in both
natural and global coordinate systems.
7. i)Obtain the shape functions of a nine noded quadrilateral element.
ii) Describe auto and adaptive mesh generation techniques.
8. i) Comment on discretization error with an example.
ii) Discuss p and h methods of refinement and give applications of each method.
9. Determine the shape functions of six noded triangular elements.
10. What is adoptive meshing? Explain any one algorithm for auto meshing.
UNIT-IV
2. Write a detailed note on Matrix solution techniques and Natural coordinate systems.
3. Explain the problems involved in the analysis of material non linearity and explain how a
solution procedure for search problems may be established for structures made of ductile
materials.
4. Explain how the consistent mass matrix for a pinpointed bar element is obtained.
5. Discuss the vector iteration method of Eigen value problems.
6. The cross section of a bar of length 10cm is rectangular in section of width 3cm and depth
1 cm. The bar is subjected to forced convection over its length due to flow of fluid at
temperature of 25OC. The convection coefficient is 5 W/cm2o C. Compute the thermal load
vector due to convection.
7. Determine Pcritical for a pinned pinned column of length, L. Consider 3 elements for the
analysis. And also draw mode shape.
8.i) Explain characteristic polynomial technique of Eigen value- Eigen vector evaluation.
ii) Explain any one method of handling geometric non- linearity.
9. Determine the Eigen values and Eigen vectors for the stepped bar shown in fig. Take
E= 2 x 105 N/mm2 specific weight = 78.5 kN/m3.
A1=600mm2
250mm
A2=300mm2
125mm
10. Using lumped mass approach obtain natural frequencies and shapes of flexural modes of
a fixed end beam of span 900mm.
A= 500mm2, I= 400mm2, p=7840 kg/m3 and E=200 kN/mm2
UNIT-V
1. List out the two advantages of post processing.
2. List out FEM software packages.
3. Write the equation for calculating element mass matrix in terms of shape functions.
4. State the functions of a preprocessor in FEA software package.
5. Name few software packages used for finite element analysis.
6. Specify degrees of freedom for SOLID 45 element analysis.
7. Give one dimensional heat flow equation.
1m
6. How will you solve this problem (refer fig) using finite element software? Detailed
steps should be clearly specified.
P
7.Write the step by step procedure of solving a structural problem by using any finite element
software.
8. Obtain the element matrices for the one dimensional heat conduction equation.
d/dx (KdT/dx)+Q=0 subject to boundary conditions
T/(x=0) = To q/(x=c) =n (TL-T)
9. Find the forces developed at the fixed end A due to rise of 60oC.
A= 800mm2
E= 200 kN/mm2
=10x 10-6/oC
2m
10. Find the forces in the members due to rise in temperature of 30oC.
A= 400mm2, E= 200kN/mm2 and = 70x 10-6/oC.
2m