AERONAUTICAL ENGINEERING – MRCET (UGC – Autonomous)

FINITE ELEMENT ANALYSIS

COURSE FILE
III B. Tech II Semester
(2017-2018)
Prepared By
Mr. S. SHAILESH BABU
Asst. Prof.

Department of Aeronautical Engineering

MALLA REDDY COLLEGE OF ENGINEERING &

TECHNOLOGY
(Autonomous Institution – UGC, Govt. of India)
Affiliated to JNTU, Hyderabad, Approved by AICTE - Accredited by NBA & NAAC – ‘A’ Grade - ISO 9001:2015
Certified)
Maisammaguda, Dhulapally (Post Via. Kompally), Secunderabad – 500100, Telangana State, India.

III– II B. Tech Finite Element Analysis By S. Shailesh Babu

 AERONAUTICAL ENGINEERING – MRCET (UGC – Autonomous)

MRCET VISION

 To become a model institution in the fields of Engineering, Technology and

Management.

 To have a perfect synchronization of the ideologies of MRCET with challenging
demands of International Pioneering Organizations.

MRCET MISSION

To establish a pedestal for the integral innovation, team spirit, originality and

competence in the students, expose them to face the global challenges and become

pioneers of Indian vision of modern society.

MRCET QUALITY POLICY.

 To pursue continual improvement of teaching learning process of Undergraduate and

Post Graduate programs in Engineering & Management vigorously.

 To provide state of art infrastructure and expertise to impart the quality education.

III– II B. Tech Finite Element Analysis By S. Shailesh Babu

 AERONAUTICAL ENGINEERING – MRCET (UGC – Autonomous)

PROGRAM OUTCOMES
(PO’s)
Engineering Graduates will be able to:
1. Engineering knowledge: Apply the knowledge of mathematics, science, engineering
fundamentals, and an engineering specialization to the solution of complex engineering
problems.
2. Problem analysis: Identify, formulate, review research literature, and analyze complex
engineering problems reaching substantiated conclusions using first principles of
mathematics, natural sciences, and engineering sciences.
3. Design / development of solutions: Design solutions for complex engineering problems
and design system components or processes that meet the specified needs with
appropriate consideration for the public health and safety, and the cultural, societal,
and environmental considerations.
4. Conduct investigations of complex problems: Use research-based knowledge and
research methods including design of experiments, analysis and interpretation of data,
and synthesis of the information to provide valid conclusions.
5. Modern tool usage: Create, select, and apply appropriate techniques, resources, and
modern engineering and IT tools including prediction and modeling to complex
engineering activities with an understanding of the limitations.
6. The engineer and society: Apply reasoning informed by the contextual knowledge to
assess societal, health, safety, legal and cultural issues and the consequent
responsibilities relevant to the professional engineering practice.
7. Environment and sustainability: Understand the impact of the professional engineering
solutions in societal and environmental contexts, and demonstrate the knowledge of,
and need for sustainable development.
8. Ethics: Apply ethical principles and commit to professional ethics and responsibilities
and norms of the engineering practice.
9. Individual and team work: Function effectively as an individual, and as a member or
leader in diverse teams, and in multidisciplinary settings.
10. Communication: Communicate effectively on complex engineering activities with the
engineering community and with society at large, such as, being able to comprehend
and write effective reports and design documentation, make effective presentations,
and give and receive clear instructions.
11. Project management and finance: Demonstrate knowledge and understanding of the
engineering and management principles and apply these to one’s own work, as a
member and leader in a team, to manage projects and in multi disciplinary
environments.
12. Life- long learning: Recognize the need for, and have the preparation and ability to
engage in independent and life-long learning in the broadest context of technological
change.

III– II B. Tech Finite Element Analysis By S. Shailesh Babu

 AERONAUTICAL ENGINEERING – MRCET (UGC – Autonomous)

DEPARTMENT OF AERONAUTICAL ENGINEERING

VISION

Department of Aeronautical Engineering aims to be indispensable source in Aeronautical

Engineering which has a zeal to provide the value driven platform for the students to
acquire knowledge and empower themselves to shoulder higher responsibility in building
a strong nation.

MISSION

The primary mission of the department is to promote engineering education and

research. To strive consistently to provide quality education, keeping in pace with time
and technology. Department passions to integrate the intellectual, spiritual, ethical and
social development of the students for shaping them into dynamic engineers.

QUALITY POLICY STATEMENT

Impart up-to-date knowledge to the students in Aeronautical area to make them quality
engineers. Make the students experience the applications on quality equipment and
tools. Provide systems, resources and training opportunities to achieve continuous
improvement. Maintain global standards in education, training and services.

III– II B. Tech Finite Element Analysis By S. Shailesh Babu

 AERONAUTICAL ENGINEERING – MRCET (UGC – Autonomous)

PROGRAM EDUCATIONAL OBJECTIVES – Aeronautical Engineering

1. PEO1 (PROFESSIONALISM & CITIZENSHIP): To create and sustain a community of
learning in which students acquire knowledge and learn to apply it professionally with
due consideration for ethical, ecological and economic issues.
2. PEO2 (TECHNICAL ACCOMPLISHMENTS): To provide knowledge based services to
satisfy the needs of society and the industry by providing hands on experience in
various technologies in core field.
3. PEO3 (INVENTION, INNOVATION AND CREATIVITY): To make the students to design,
experiment, analyze, and interpret in the core field with the help of other multi
disciplinary concepts wherever applicable.
4. PEO4 (PROFESSIONAL DEVELOPMENT): To educate the students to disseminate
research findings with good soft skills and become a successful entrepreneur.
5. PEO5 (HUMAN RESOURCE DEVELOPMENT): To graduate the students in building
national capabilities in technology, education and research

PROGRAM SPECIFIC OUTCOMES – Aeronautical Engineering

1. To mould students to become a professional with all necessary skills, personality and
sound knowledge in basic and advance technological areas.
2. To promote understanding of concepts and develop ability in design manufacture and
maintenance of aircraft, aerospace vehicles and associated equipment and develop
application capability of the concepts sciences to engineering design and processes.
3. Understanding the current scenario in the field of aeronautics and acquire ability to
apply knowledge of engineering, science and mathematics to design and conduct
experiments in the field of Aeronautical Engineering.
4. To develop leadership skills in our students necessary to shape the social, intellectual,
business and technical worlds.

III– II B. Tech Finite Element Analysis By S. Shailesh Babu

 AERONAUTICAL ENGINEERING – MRCET (UGC – Autonomous)

MALLA REDDY COLLEGE OF ENGINEERING & TECHNOLOGY

L T/P/D C
III Year B. Tech, ANE-II Sem

4 0 4

(R15A2120) FINITE ELEMENT ANALYSIS

(CORE ELECTIVE – IV)
Objectives:
 It covers the fundamental theoretical approach beginning with a review of differential equations, boundary
conditions, integral forms, interpolation, parametric geometry, numerical integration, and matrix algebra.
 Next, engineering applications to field analysis, stress analysis and vibrations are introduced. Time
dependent problems are also treated.
 Students are also introduced, by means of selected tutorials, to the commercial finite element system
SolidWorks which is similar to one they could be expected to use upon graduation. Graduate students will
also be introduced to the more powerful (and difficult to use) Ansys system.

UNIT – I
Introduction to Finite Element Method for solving field problems. Stress and Equilibrium. Strain – Displacement
relations. Stress – strain relations. One Dimensional problems : Finite element modeling coordinates and shape
functions. Potential Energy approach : Assembly of Global stiffness matrix and load vector. Finite element
equations, Treatment of boundary conditions, Quadratic shape functions.
UNIT – II
Analysis of Beams : Element stiffness matrix for two node, two degrees of freedom per node beam element. Finite
element modelling of two dimensional stress analysis with constant strain triangles and treatment of boundary
conditions.
UNIT – III
Finite element modelling of Axisymmetric solids subjected to Axisymmetric loading with triangular elements. Two
dimensional four noded isoparametric elements and numerical integration.
UNIT – IV
Steady state heat transfer analysis : one dimensional analysis of a fin and two dimensional analysis of thin plate.
Analysis of a uniform shaft subjected to torsion.
UNIT-V
Dynamic Analysis : Formulation of finite element model, element matrices, evaluation of Eigen values and Eigen
vectors for a stepped bar and a beam.

Outcomes:
 Upon completion of the course students should be able to correlate a differential equation and its
equivalent integral form.
 Understand parametric interpolation and parametric geometry enforce essential boundary conditions to a
matrix system.
TEXT BOOK:
1. Introduction to Finite Elements in Engineering / Chandraputla, Ashok and Belegundu /
Prentice – Hall.
2. The Finite Element Methods in Engineering / SS Rao / Pergamon.
3. The Finite Element Method for Engineers – Kenneth H. Huebner, Donald L. Dewhirst, Douglas E. Smith and
Ted G. Byrom / John Wiley & sons (ASIA) Pte Ltd.

REFERENCES:
1. An introduction to Finite Element Method / JN Reddy / Me Graw Hill
2. Finite Element Methods/ Alavala/TMH
3. Finite Element Analysis/ C.S.Krishna Murthy

III– II B. Tech Finite Element Analysis By S. Shailesh Babu

 UNIT I
INTRODUCTION: The finite element analysis is a numerical technique. In this method all the complexities
of the problems, like varying shape, boundary conditions and loads are maintained as they are but the
solutions obtained are approximate. The fast improvements in computer hardware technology and slashing
of cost of computers have boosted this method, since the computer is the basic need for the application of
this method. A number of popular brand of finite element analysis packages are now available commercially
Some of the popular packages are STAAD-PRO, GT-STRUDEL, NASTRAN, NISA and ANSYS. Using
these packages one can analyze several complex structures.

The finite element analysis originated as a method of stress analysis in the design of aircrafts. It started as an
extension of matrix method of structural analysis. Today this method is used not only for the analysis in
solid mechanics, but even in the analysis of fluid flow, heat transfer, electric and magnetic fields and many
others. Civil engineers use this method extensively for the analysis of beams, space frames, plates, shells,
folded plates, foundations, rock mechanics problems and seepage analysis of fluid through porous media.
Both static and dynamic problems can be handled by finite element analysis. This method is used extensively
for the analysis and design of ships, aircrafts, space crafts, electric motors and heat engines

The basic unknowns or the Field variables which are encountered in the engineering problems are
displacements in solid mechanics, velocities in fluid mechanics, electric and magnetic potentials in electrical
engineering and temperatures in heat flow problems In a continuum, these unknowns are infinite. The finite
element procedure reduces such unknowns to a finite number by dividing the solution region into small
parts called elements and by expressing the unknown field variables in terms of assumed approximating
functions (Interpolating functions/Shape functions) within each element. The approximating functions are
defined in terms of field variables of specified points called nodes or nodal points. Thus in the finite
element analysis the unknowns are the field variables of the nodal points. Once these are found the field
variables at any point can be found by using interpolation functions. After selecting elements and nodal
unknowns next step in finite element analysis is to assemble element properties for each element. For
example, in solid mechanics, we have to find the force-displacement i.e. stiffness characteristics of each
individual element. Mathematically this relationship is of the form
[k]e {δ}e = {F}e
where [k]e is element stiffness matrix, {δ }e is nodal displacement vector of the element and {F}e is nodal
force vector. The element of stiffness matrix kij represent the force in coordinate direction ‘i’ due to a unit
displacement in coordinate direction ‘j’. Four methods are available for formulating these element properties
viz. direct approach, variational approach, weighted residual approach and energy balance approach. Any
one of these methods can be used for assembling element properties. In solid mechanics variational
approach is commonly employed to assemble stiffness matrix and nodal force vector (consistant loads).
Element properties are used to assemble global properties/structure properties to get system equations [k]
{} = {F}. Then the boundary conditions are imposed. The solution of these simultaneous equations give
the nodal unknowns. Using these nodal values additional calculations are made to get the required values
e.g. stresses, strains, moments, etc. in solid mechanics problems.

Thus the various steps involved in the finite element analysis are:
(i) Select suitable field variables and the elements.
(ii) Discritise the continua.
(iii) Select interpolation functions.
(iv) Find the element properties.
(v) Assemble element properties to get global properties.
(vi) Impose the boundary conditions.
(vii) Solve the system equations to get the nodal unknowns.
(viii) Make the additional calculations to get the required values.

Methods of Engineering Analysis

There are three methods are adopted for analyzing the product
1.Experimental methods
2.Analytical methods
Numerical methods
Experimental methods
In these methods the actual products or their proto type models or atleast their material specimen are tested
by using some equipments
Ex: UTM, Rockwell hardness tester
Analytical methods
These methods are theoretically analyzing methods. Only simple and regular shaped products like beams,
shafts, plates can be analyzed by these methods
Numerical methods
For the products of complicated sizes and shapes with complicated material properties and boundary
conditions getting solution using analytical methods is highly difficult. In such situation the numerical
method can be employed
There are three numerical methods
i)Functional approximating methods
ii) Finite element method
iii)Finite difference method

Application of FEM

S.No Area of Study Analysing problem

1 Civil Engineering structures Analysis of trusses, folded plates, shell roofs, bridges and
prestressed concrete structures
2 Aircraft structures Analysis of aircraft wings, fins, rockets, space craft and
missile structures
3 Mechanical Design Stress analysis of pressure vessels, pistons, composite
materials,Linkages and gears
4 Heat Conduction Temperature distribution in solida and fluids
5 Hydraulic and water resources Analysis of potential flows,free suface fkows,viscous
engineering flows,analysis of hydraulic structures and dams
6 Electrical Machines and Analysis of synchronous and induction machines eddy current
electromagnetics and core losses in electric machines
7 Nuclear Engineering Analysis of nuclear pressure vessels and containment
structures
8 Geomechanics Stress analysis in soils,dams,layered piles and machine
foundations

Advantages and disadvantages of FEM

Advantages
Using FEM we are able to
1.model irregular shaped bodies quite easily
2.handle general load conditions without difficulty
3.model bodies composed of several different materials because the element equations are evaluated
individually
4. handle unlimited numbers and kinds of boundary conditions
5. vary the size of the element to make it possible to use small elements
6. alter the finite element model easily and cheaply
7. include dynamic effects
Disadvantages
1.The finite element method is time consuming process
2. FEM cannot produce exact results as those of analytical methods

Equations of Equilibrium for 3D Body

Typical three dimensional element of size dx × dy × dz. Face abcd may be called as negative face of x and
the face efgh as the positive face of x since the x value for face abcd is less than that for the face efgh.
Similarly the face aehd is negative face of y and bfgc is positive face of y. Negative and
positive faces of z are dhgc and aefb. The direct stresses σ and shearing stresses τ acting on the negative
faces are shown in the Fig. with suitable subscript. It may be noted that the first subscript of shearing stress
is the plane and the second subscript is the direction. Thus the τ xy means shearing stress on the plane where
x value is constant and y is the direction.
Let the intensity of body forces acting on the element in x, y, z directions be X, Y and Z respectively as
shown in Fig The intensity of body forces are uniform over entire body. Hence the total body force in x,
y, z direction on the element shown are given by
(i) X dx dy dz in x – direction
(ii) Y dx dy dz in y – direction and
(iii) Z dx dy dz in z – direction

Equations of Equilibrium
Considering all forces are acting we can write the equilibrium equations for the element

Simplifying and dividing throughout by dx dy dz

 Similarly Fy=0 and Fz=0 Equilibrium conditions give

Mx=0

Neglecting small quantity then

zy= yz

My=0 then we will get

xz= zx

MZ=0 then we will get

xy= yx
[] =[x y z  xy yz  xz]
T

and the equilibrium equations are

Strain Displacement equations

Taking displacement components in x, y ,z directions as u, v, and w respectively, the relations among

components of strain and the components of displacement are

strains are expressed up to the accuracy of second order (quadratic) changes in displacements. These
equations may be simplified to the first (linear) order accuracy only by dropping the second order changes
terms. Then linear strain – displacement relation is given by:
LINEAR CONSTITUTIVE EQUATIONS

The constitutive law expresses the relationship among stresses and strains. In theory of elasticity, usually it
is considered as linear. In one dimensional stress analysis, the linear constitutive law is stress is proportional
to strain and the constant of proportionality is called Young’s modulus. It is very well known as Hooke’s
law.
The similar relation is expressed among the six components of stresses and strains and is called
‘Generalized Hookes Law”. This may be stated as:

where D is 6 × 6 matrix of constants of elasticity to be determined by experimental investigations for each

material. As D is symmetric matrix [Dij = Dji], there are 21 material properties for linear elastic
Anisotropic Materials. Certain materials exhibit symmetry with respect to planes within the body. Such
materials are called Ortho tropic materials. Hence for orthotropic materials, the number of material
constants reduce to 9 as shown below:
Note that there are 12 material properties in above equations However only nine of these are independent
because the following relations exist

For Isotropic Materials the above set of equations are further simplified. An isotropic material is the one
that has same material property in all directions. In other word for isotropic materials,

Hence for a three dimensional problem, the strain stress relation for isotropic material is,
PLANE STRESS PROBLEM
The thin plates subject to forces in their plane only, fall under this category of the problems. Fig. shows a
typical plane stress problem. In this, there is
PLANE STRAIN PROBLEM

A long body subject to significant lateral forces but very little longitudinal forces falls under this category of
problems. Examples of such problems are pipes, long strip footings, retaining walls, gravity dams, tunnels,
etc. In these problems, except for a small distance at the ends, state of stress is represented by
any small longitudinal strip. The displacement in longitudinal direction (z-direction) is zero in typical strip.
Hence the strain components,
Functional Approximation Methods

The nature of the problems for which the solutions to be found out are

i) Equilibrium problems

ii)Eigen value problems

iii)probagation problems

The functional approximation methods for solving the above types of problems are classified in to major
types

i)Variational methods

ii)Weighted residual methods

Rayleigh-Ritz method is good example for variational method

Weighted residual method

 Point collocation method

 sub domain collocation method
 Least square method
 Galerkin's method

Rayleigh-Ritz Method

Rayleigh -Ritz method is a typical variational method in which principle of integral approach is adopted for
solving the complex structural problems

i) Minimum potential energy method

ii)Integral approach method

Minimum potential energy method

In this method the total potential energy '' is considered as the function of generalized coordinated which
are exactly equal to the number of degrees of freedom

= U-W

U=Internal energy

W=work done by the external force

Polynomial series

y(x)=a1+a2 x+a3 x2+ ---------------------

a1, a2, a3--------------- are Ritz parameters

Integral approach method

Differential equation is
𝑑2𝑦
D 𝑑𝑥2 +Q=0

𝑙
I∫0 [𝐷/2(𝑑𝑦/𝑑𝑥)𝟐-Qy ] dx
ONE DIMENSIONAL PROBLEMS

Bar and beam elements are considered as One Dimensional elements. These elements are often used to
model trusses and frame structures

Types of Loading
i) Body force (f)

It is a distributed force acting on every elemental volume of the body. Unit is Force / Unit volume. Ex: Self
weight due to gravity.

ii) Traction (T)

It is a distributed force acting on the surface of the body. Unit is Force / Unit area. But for one dimensional
problem, unit is Force / Unit length. Ex: Frictional resistance, viscous drag and Surface shear.

iii) Point load (P)

It is a force acting at a particular point which causes displacement.

Finite Element Modeling

It has two processes. (1) Discretization of structure (2) Numbering of nodes.

CO – ORDINATES
(A) Global co – ordinates, (B) Local co – ordinates and (C) Natural co – ordinates.

 Equation of Stiffness Matrix for One dimensional bar element

For a stepped bar loaded as shown in figure.Determine a) Nodal displacements

b) support Reactions c)Element Sress

Solution

Element 1 Element 2 Element 3

A1= 15 cm2 A2= 15 cm2 A3= 24 cm2
E1= 20 X 106 N/cm2 E2= 20 X 106 N/cm2 E3= 20 X 106 N/cm2
L1= 75 cm L2= 75 cm L3= 60 cm
α 1= 11 X 10 -6 cm/cm0 C α 2= 11 X 10 -6 cm/cm0 C α 3= 11 X 10 -6 cm/cm0 C
∆T = 10 0C ∆T = 10 0C ∆T = 10 0C

F0(1)= A1 E1α 1∆T=33000 N

F0(2)= A2 E2α 2∆T=33000 N

F0(3)= A3 E3α 3∆T=52800 N

The Nodal Forces are

F1=R1+P- F0(1)= R1-33000

F2=P2 + F0(1) - F0(2) = 10000

F3=P3 + F0(2) - F0(3)= - 39800

F4= R4 + P3 + F0(2) - F0(3) = R4 + 52800

The stiffness values are

k1=A1E1/L1=4 X 106 N/cm

k2=A2E2/L2=4 X 106 N/cm

k3=A3E3/L3=8X 106 N/cm

the nodal conditions are u1=0 and u4=0

solve the above matrix then you will get the values of u 2 and u3 as - 3.48 X 10-3 cm and as - 0.49 X 10-1 cm

R1 = 34960 N

R4 = - 24960 N

r(1) =  (1) - 0(1) = -2330.7 N/cm

2

r(2) =  (2) - 0(2) = -2997.3 N/cm2

r(3) =  (3) - 0(3) = -1010 N/cm2

 CHAPTER – II
One could obtain the global stiffness matrix of a continuous beam from assembling member stiffness matrix of
individual beam elements. Towards this end, we break the given beam into a number of beam elements. The
stiffness matrix of each individual beam element can be written very easily. For example, consider a
continuous beam ABCD as shown in Fig. 1a. The given continuous beam is divided into three beam elements
as shown in Fig. 1b. It is noticed that, in this case, nodes are located at the supports. Thus each span is treated
as an individual beam. However sometimes it is required to consider a node between support points. This is
done whenever the cross sectional area changes suddenly or if it is required to calculate vertical or rotational
displacements at an intermediate point. Such a division is shown in Fig. 1c. If the axial deformations are
neglected then each node of the beam will have two degrees of freedom: a vertical displacement
(corresponding to shear) and a rotation (corresponding to bending moment). In Fig. 1b, numbers enclosed in
a circle represents beam numbers. The beam ABCD is divided into three beam members. Hence, there are four
nodes and eight degrees of freedom. The possible displacement degrees of freedom of the beam are also
shown in the figure. Let us use lower numbers to denote unknown degrees of freedom (unconstrained degrees
of freedom) and higher numbers to denote known (constrained) degrees of freedom. Such a method of
identification is adopted in this course for the ease of imposing boundary conditions directly on the structure
stiffness matrix. However, one could number sequentially as shown in Fig. 1d. This is preferred while solving
the problem on a computer.
In the above figures, single headed arrows are used to indicate translational and double headed arrows are used
to indicate rotational degrees of freedom.

Beam Stiffness Matrix:

Fig. 2 shows a prismatic beam of a constant cross section that is fully restrained at ends in local orthogonal co-
ordinate system x' y' z'. The beam ends are denoted by nodes j and k. The x' axis coincides with the centroidal
axis of the member with the positive sense being defined from j to k. Let L be the length of the member, A area
of cross section of the member and I zz is the moment of inertia about z' axis.
Two degrees of freedom (one translation and one rotation) are considered at each end of the member. Hence,
there are four possible degrees of freedom for this member and hence the resulting stiffness matrix is of the
order 4 X 4. In this method counterclockwise moments and counterclockwise rotations are taken as positive.
The positive sense of the translation and rotation are also shown in the figure. Displacements are considered as
positive in the direction of the co- ordinate axis. The elements of the stiffness matrix indicate the forces
exerted on the member by the restraints at the ends of the member when unit displacements are imposed at
each end of the member. Let us calculate the forces developed in the above beam member when unit
displacement is imposed along each degree of freedom holding all other displacements to zero. Now impose a
unit displacement along y' axis at j end of the member while holding all other displacements to zero as shown
in Fig.a. This displacement causes both shear and moment in the beam. The restraint actions are also shown in
the figure. By definition they are elements of the member stiffness matrix. In particular they form the first
column of element stiffness matrix.
In Fig.b, the unit rotation in the positive sense is imposed at j end of the beam while holding all other
displacements to zero. The restraint actions are shown in the figure. The restraint actions at ends are
calculated referring to tables given in lesson …
In Fig. 3c, unit displacement along y' axis at end k is imposed and corresponding restraint actions are
calculated. Similarly in Fig.d, unit rotation about z' axis at end k is imposed and corresponding stiffness
coefficients are calculated. Hence the member stiffness matrix for the beam member is

The stiffness matrix is symmetrical. The stiffness matrix is partitioned to separate the actions associated with
two ends of the member. For continuous beam problem, if the supports are unyielding, then only rotational
degree of freedom shown in Fig. is possible. In such a case the first and the third rows and columns will be
deleted. The reduced stiffness matrix will be,


Instead of imposing unit displacement along y' at j end of the member in Fig.a, apply a displacement u'1 along y'
at j end of the member as shown in Fig. a, holding all other displacements to zero. Let the restraining forces
developed be denoted by q11, q21, q31 and q41.

The forces are equal to,

q11 = k11u'1; q21 = k21u'1; q31 = k31u'1; q41 = k41u'1

Now, give displacements u'1, u'2, u'3 and u'4 simultaneously along displacement degrees of freedom 1, 2, 3 and 4
respectively. Let the restraining forces developed at member ends be q1, q2 , q3 and q4 respectively as shown
in Fig. b along respective degrees of freedom. Then by the principle of superposition, the force displacement
relationship can be written as,
Beam (global) Stiffness Matrix:
The formation of structure (beam) stiffness matrix from its member stiffness matrices is explained with help of
two span continuous beams shown in Fig. a. Note that no loading is shown on the beam. The orthogonal co-

ordinate system xyz denotes the global co-ordinate system.

For the case of continuous beam, the x - and x’ - axes are collinear and other axes ( y and y' , z and z' ) are
parallel to each other. Hence it is not required to transform member stiffness matrix from local co-ordinate
system to global coordinate system as done in the case of trusses. For obtaining the global stiffness matrix, first
assume that all joints are restrained. The node and member numbering for the possible degrees of freedom are
shown in Fig b. The continuous beam is divided into two beam members. For this member there are six
possible degrees of freedom. Also in the figure, each beam member with its displacement degrees of freedom
(in local co ordinate system) is also shown. Since the continuous beam has the same moment of inertia and span,
the member stiffness matrix of element 1 and 2 are the same. They are,

The local and the global degrees of freedom are also indicated on the top and side of the element stiffness
matrix. This will help us to place the elements of the element stiffness matrix at the appropriate locations of the
global stiffness matrix. The continuous beam has six degrees of freedom and hence the stiffness matrix is of the
order6. Let [K] denotes the continuous beam stiffness matrix of order 6X6. From Fig., [K] may be written as,
The 4X4 upper left hand section receives contribution from member 1 (AB) and 4X4 lower right hand section of
global stiffness matrix receives contribution from member 2. The element of the global stiffness matrix
corresponding to global degrees of freedom 3 and 4 receives element from both members 1 and 2.

FORMATION OF LOAD VECTOR:

Consider a continuous beam ABC as shown in Fig.

We have two types of load: member loads and joint loads. Joint loads could be handled very easily as done in
case of trusses. Note that stiffness matrix of each member was developed for end loading only. Thus it is
required to replace the member loads by equivalent joint loads. The equivalent joint loads must be evaluated
such that the displacements produced by them in the beam should be the same as the displacements produced
by the actual loading on the beam. This is evaluated by invoking the method of superposition.
The loading on the beam shown in Fig. (a), is equal to the sum of Fig. (b) and Fig. (c). In Fig. (c), the joints are
restrained against displacements and fixed end forces are calculated. In Fig. (c) these fixed end actions are
shown in reverse direction on the actual beam without any load. Since the beam in Fig. (b) is restrained (fixed)
against any displacement, the displacements produced by the joint loads in Fig. (c) must be equal to the
displacement produced by the actual beam in Fig. (a). Thus the loads shown in Fig. (c) are the equivalent joint
loads .Let, p1, p2 , p3 , p4 , p5 and p6 be the equivalent joint loads acting on the continuous beam along
displacement degrees of freedom 1,2,3,4,5 and 6 respectively as shown in Fig. (b). Thus the global load vector is,
 SOLUTION OF EQUILIBRIUM EQUATIONS:
After establishing the global stiffness matrix and load vector of the beam, the load displacement relationship for
the beam can be written as

Where is the global load vector, { P } { u } is displacement vector and is the global stiffness matrix. In the
above equation some joint displacements are known from support conditions. The above equation may be
written as

Where {pk} and { uk} denote respectively vector of known forces and known displacements. And {pk} and { uk}
denote respectively vector of unknown forces and unknown displacements respectively. Now expanding
equation

{ p k }  k11 {u u }  k12 {u k }

{ p u }  k 21 {u u }  k 22 {u k }
Since u k is known, the unknown joint displacements can be evaluated. And support reactions are

evaluated from equation, after evaluating unknown displacement vector.

Let R1, R3 and R5 be the reactions along the constrained degrees of freedom. Since equivalent joint loads
are directly applied at the supports, they also need to be considered while calculating the actual
reactions. Thus,



The reactions may be calculated as follows. The reactions of the beam shown in Fig. a are equal to
the sum of reactions shown in Fig. b, Fig. c and Fig. d.
Member end actions q1, q2, q3 and q4 are calculated as follows. For example consider the first element 1
Finite Elements for 2‐D Problems
General Formula for the Stiffness Matrix

Displacements (u, v) in a plane element are interpolated from nodal displacements (ui, vi) using shape functions
Ni as follows,

where N is the shape function matrix, u the displacement vector and d the nodal displacement vector. Here we
have assumed that u depends on the nodal values of u only, and v on nodal values of v only. Most commonly
employed 2‐D elements are linear or quadratic triangles and quadrilaterals.

Constant Strain Triangle (CST or T3)

This is the simplest 2‐D element, which is also called linear triangular element.

For this element, we have three nodes at the vertices of the triangle, which are numbered around the element in
the counter clockwise direction. Each node has two degrees of freedom (can move in the x and y directions). The
displacements u and v are assumed to be linear functions within the element, that is,

where bi (i = 1, 2, ..., 6) are constants. From these, the strains are found to be,

which are constant throughout the element.

Notes:
 ∙ Q4 and T3 are usually used together in a mesh with linear elements.
 ∙ Q8 and T6 are usually applied in a mesh composed of quadratic elements.
 ∙ Quadratic elements are preferred for stress analysis, because of their high accuracy and the flexibility in
modelling complex geometry, such as curved boundaries.
 Unit – III
Axisymmetric solids subjected to Axisymmetric loading with triangular elements. Two dimensional four noded
isoparametric elements and numerical integration.

Elasticity Equations
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. The types of elasticity equations are
1. Strian – Displacement relationship equations

2. Sterss – Strain relationship equation

1
σ – Stress, τ – Shear Stress, E – Young’s Modulus, v – Poisson’s Ratio, e – Strain, γ- Shear Strain.
3. Equilibrium equations

σ – Stress, τ – Shear Stress, B x - Body force at X direction,

B y - Body force at Y direction, B z - Body force at Z direction.

4. Compatibility equations

There are six independent compatibility equations, one of which is

The other five equations are similarly second order relations.

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 xisymmetric element.

2











Axisymmetric Formulation

The displacement vector u is given by

The stress σ is given by

3
The strain e is given by

Equation of shape function for Axisymmetric element

Shape function,

 Equation of Strain – Displacement Matrix [B] for Axisymmetric element

4
 Equation of Stress – Strain Matrix [D] for Axisymmetric element



 Equation of Stiffness Matrix [K] for Axisymmetric element

 Temperature Effects

The thermal force vector is given by

5
 Problem (I set) 1. For the given element, determine the stiffness matrix. Take E=200GPa and υ= 0.25. 2.




2. For the figure, determine the element stresses. Take E=2.1x105N/mm2 and υ= 0.25. The co – ordinates are in mm. The nodal
displacements are u1=0.05mm, w1=0.03mm, u2=0.02mm, w2=0.02mm, u3=0.0mm, w3=0.0mm.

3. A long hollow cylinder of inside diameter 100mm and outside diameter 140mm is subjected to an internal pressure of 4N/mm2.
By using two elements on the 15mm length, calculate the displacements at the inner radius.

6


 Isoparametric element

Generally it is very difficult to represent the curved boundaries by straight edge elements. A large number of elements may be used to
obtain reasonable resemblance between original body and the assemblage. In order to overcome this drawback, isoparametric elements
are used.

If the number of nodes used for defining the geometry is same as number of nodes used defining the displacements, then it is known
as isoparametric element.

Superparametric element
If the number of nodes used for defining the geometry is more than number of nodes used for defining the displacements, then it is
known as superparametric element.

7
 Subparametric element

If the number of nodes used for defining the geometry is less than number of nodes used for defining the displacements, then it is
known as subparametric element.

 Equation of Shape function for 4 noded rectangular parent element

8
 Equation of Stiffness Matrix for 4 noded isoparametric quadrilateral element

9
Equation of element force vector

N – Shape function, Fx – load or force along x direction, Fy – load or force along y direction.

Numerical Integration (Gaussian Quadrature)

The Gauss quadrature is one of the numerical integration methods to calculate the definite integrals. In FEA, this Gauss quadrature
method is mostly preferred. In this method the numerical integration is achieved by the following expression ,

10
Table gives gauss points for integration from -1 to 1.

 Problem (I set) 1. Evaluate, by applying 3 point Gaussian quadrature and compare with exact solution.

11
2. Evaluate, using one point and two point Gaussian quadrature. Compare with exact solution.

3. For the isoparametric quadrilateral element shown in figure, determine the local co –ordinates of the point P which has Cartesian co-
ordinates (7, 4).

4. A four noded rectangular element is in figure. Determine (i) Jacobian matrix, (ii) Strain – Displacement matrix and (iii) Element
Stresses. Take E=2x105N/mm2,υ= 0.25, u=[0,0,0.003,0.004,0.006, 0.004,0,0] T, Ɛ= 0, ɳ=0. Assume plane stress condition.

12
 MALLAREDDY COLLEGE OF ENGINEERING AND
TECHNOLOGY
MECHANICAL DEPARTMENT SUB: FEM

MODEL PAPER 1
PART A (25 MARKS)
1. a. What is meant by Engineering analysis and specify its Types (2M)
b. What is Hermite shape function (3M)
c. Write the equilibrium equations for 3D body (2M)
d. What is coordinate system. specify the types and explain (3M)
e. What is meant by axi-symmetric problems (2M)
f. Derive the shape functions for 2D truss element (3M)
g. What is the degree of freedom for the thermal problems (2M)
h. Distinguish between CST and LST (3M)
i. Write the dynamic equation of motion for the undamped free vibrations (2M)
j. Determine the Area of the triangle A(2,2),B(7,4) ,C(3,6) (3M)

PART B (10X5=50 MARKS)

2. a) Derive the equations of equilibrium in case of a three dimensional stress system.

b) Discuss the advantages and disadvantages of FEM over

(i) Classical method
(ii) Finite difference method.

OR

3.a) Solve the differential eqationfor the physical problem expressed as d2y/dx2+100=0when 0=x=10
with boundary
condition as y(0)=0 and y(10)=0 using i) point collocation ii) sub-domain collocation iii) least
square method
iv) galarkin method

b) Write the Strain displacement equations for three dimensional system

4. a) Determine the nodal displacement, Element stresses for axially loaded bar as shown in the fig.
below
 b) Derive the strain displacement matrices for triangular element of revolving body.

OR

5 a) For the beam shown in Figure below, determine the following: a) Slopes at nodes 2 and 3

b) Vertical deflection at the mid-point of the distributed load. Consider all the elements
have E=200GPa, I=5X106 mm4

6. a) For the element shown in the figure, assemble Jacobian matrix and strain displacement matrix

b) Determine the shape functions for a 8 node quadratic quadrilateral Evaluation element(boundary
noded).

OR

7. a) Establish the shape functions for a 3 – noded triangular element

 b) Find the deformed configuration, and the maximum stress and minimum stress locations for the
rectangular plate loaded as shown in the fig. Solve the problem using 2 triangular elements. Assume
thickness = 10cm; E = 70 Gpa, and υ = 0.33

8, The composite wall consists of three materials shown in figure. The inside wall temperature is at
2000C and the outside air temperature is 500C with a convection coefficient of 10 W/m2 0C.Determine
the temperature along the composite wall

9.a) Derive one dimensional steady state heat conduction equation.

b) An axisymmetric triangular element is subjected to the loading as shown in fig. the load is
distributed throughout the circumference and normal to the boundary. Derive all the necessary
equations and derive the nodal point loads.
10. a)Determine the strain displacement matrix for the TETRAHEDRAL element as shown in fig

b) Explain the concept of numerical integration and its utility in generating Isoperimetric finite
element matrices

OR

11. a) What are the necessary requirements for convergence and explain about h- and p

requirements

b) Derive the stiffness matrix for truss element in case of linear and quadratic shape functions At
20 C an axial load P = 300 x 103 is applied to the rod as shown in Fig. The temperature is then raised
0

to 600C. Assemble the element stiffness matrix and the element temperature force matrix (F).Again
determine the nodal displacements and element stresses. Also find element strains. Assume E1= 70 x
109 N/m2, 1= 900 mm2, α1 = 23 x 10-6/ 0C, E2 = 200 X 109 N/m2 , A2 = 1200

mm2 ,α 2 = 11.7 x 10-6 / 0C.

 MALLAREDDY COLLEGE OF ENGINEERING AND
TECHNOLOGY
MECHANICAL DEPARTMENT SUB: FEM

MODEL PAPER 2
PART A (25 MARKS)
1. a. What is FEM (2M)
b. Write the advantages of FEM (3M)
c. What is CST (2M)
d. Write the strain relations of three dimensional system (3M)
e. What is local coordinate System (2M)
f. Write the Eigen values and Eigen vectors for a stepped bar (3M)
g. What is the degree of freedom for the thermal problems (2M)
h. Define principle of virtual work. Describe the FEM formulation for 1D bar element (3M)
i What is dynamic analysis (2M)
j. Discuss Mess generation(3M)

PART B (10X5=50 MARKS)

2.. a) Differentiate among Bar element, Truss element and Beam element indicating D.O.F and
geometry
characteristics.
b) An axial load P =300X103N is applied at 200 C to the rod as shown in Figure below. The
temperature is the raised to 600 C a) Assemble the K and F matrices
b) Determine the nodal displacements and stresses

200 mm 300 mm

E1=70X109N/mm2 E2=200X109N/mm2
A1=900mm2 A2=1200mm2

α1=23X10-6 α2=11.7X10-6

3. a) Discuss in detail about the concepts of FEM formulation How is that emerged as powerful
tool. Discuss in detail about applications of finite element method
b) Derive an equation for finding out the potential energy by Rayleigh –Ritz method Using
Rayleigh – Ritz method, find the displacement of the midpoint of the rod shown in Fig.
Assume E = 1,A=1,ρg=1 by using linear and quadratic shape function concept
 4 a) Discuss in detail about Linear and Quadratic shape functions with examples
b) Consider axial vibration of the Aluminum bar shown in Fig., (i) develop the global stiffness and
(ii) determine the nodal displacements and stresses using elimination approach and with help of linear
and quadratic shape function concept. Assume Young’s Modulus E = 70Gpa

OR
5. a)Describe Rayleigh-Ritz method
b) Abeam is fixed at one end and supported by roller at the other end has 20KN load
applied at the center of the span of 10m.Calculate deflection and slope and also
construct shear fore and bending moment diagrams
6. a) State the properties and applications of CST

b)The nodal coordinates of the triangular element shown in figure at the interior point P.the x
coordinate is 3.3 and the shape function at node 1 is N1is 0.3.determine the shape functions at nodes 2
and 3 also find the 'y' coordinate of P

OR

7. a) Determine the stiffness and Jacobian matrix for the iso parametric quadrilateral element starting
from fundamentals.
b) Differentiate between axi- symmetric boundary condition and polar symmetric
boundary condition.
c) Derive the load vector for the axi-symmetric triangular element with the variable surface load on
the surface.
8. a) derive one dimensional steady state heat conduction equation
b)An axi symmetric element subjected to loading as shown in figure .The load is distributed
througout the circumference and normal to the boundary. Derive all necessary equations and derive
nodal point loads

OR

9. Calculate the conductance matrix [K(e)] and load vector fF(e)g for the triangle
element shown in figure 8. The thermal conductivities are kx = ky = 4 W=cm�0C
and h = 0.3 W/cm2 0C. Thickness of the element is 1cm. All coordinates are given
in cms. Convection occurs on the side joining modes i and j.

10. For the stepped bar shown in figure develope the global stiffness matrix and mass matrices
and determine the natural frequencies and mode shapes Assume E=200GPa and mass density
is 7850 Kg/m3 L1=L2=0.3 m A1=350 mm2and A2 =600 mm2

OR
11. a)Derive the shape functions for the four noded tetrahedron element from the first principles
b)discuss the importance of semi automatic meshing and practical applications
 MALLAREDDY COLLEGE OF ENGINEERING
AND TECHNOLOGY
MECHANICAL DEPARTMENT SUB: FEM

MODEL PAPER 3
PART A (25 MARKS)
1. a. List the various weighted residual methods (2M)
b. Write the properties of shape function (3M)
c. What the advantages of natural coordinate system (2M)
d. Write analogies between structural, heat transfer and fluid mechanics (3M)
e. Name few FEA packages (2M)
f. Derive the mass matrix for a 1D linear bar element (3M)
g. What are the properties of stiffness matrix (2M)
h. Explain about plain stress and plain strain conditions (3M)
i. Write down the conduction matrix for a three noded triangular element (2M)
j. Distinguish between Error in solution and Residual (3M)

PART B (10X5=50)

2. a)Determint the nodal displacement ,stress and strain for the bar shown in fig

b) Using potential energy approach, describe FE formulation for plane truss Element

OR

3. a) Solve the differential equation for the physical problem expressed as d2y/dx2+100=0when
0=x=10 with boundary condition as y(0)=0 and y(10)=0 using i) point collocation ii) sub-domain
collocation iii) least square method
 iv) galarkin method.

b) Explain the concept of FEM briefly .outline the steps involved in FEM along with applications

4. a)For a beam and loading shown in fig., determine the slopes at 2 and 3 and the vertical deflection
at the midpoint of the distributed load

b) Establish the shape functions for a 3 – noded triangular element.

OR

5. Calculate the nodal displacement, stresses and support reactions for the truss shown in figure

6.a) Evaluate the element stiffness matrix for the triangular element shown under plane strain
condition. Assume the following values E=200 GPa, µ=0.25, t=1 mm
b) For the element shown in the figure, assemble Jacobian matrix and strain displacement matrix

OR

7. a) Derive the a)shape function and b) strain displacement matrices for triangular element of
revolving body

b) for the Isoparametric quadrilateral element shown in fig , determine the local co-ordinates of the
point P whose Cartesian co=ordinates as(6,4)

8 a) Determine the temperature at the nodal interfaces for the two layered wall shown in fig.the left
face is supplied with

heat flux of Q11=5 W/cm2 and the right face is maintained at 200C

Qu2=5w/c
m
K1=0.2W/cm0C0
K2=0.06W/cm C
A=1cm2.
b) Derive the Strain displacement Matrix for 2D-Thin plate. Consider the temperature field with in the
triangular element is given by T= N1T1 + N2T2 + N3T3

OR

9. Determine the temperature distribution through the composite wall shown in figure, when
convection heat loss occurs on the left surface. Assume unit area Assume wall thickness t1 = 4cm, t2 =
2cm, k1 = 0.5w/cm0c, k2 = 0.05w/cm0c

h= 0.1w/cm2 0c and Tα= -50c

10. a) Determine the eigen values and the associated Eigen vectors of the matrix [A] given by

b) State the properties of Eigen Values.

OR

11.a) Explain difference between Lumped Mass and Consistent Mass

b) Determine the Natural frequency of the beam shown in the figure

 MALLAREDDY COLLEGE OF ENGINEERING
AND TECHNOLOGY
MECHANICAL DEPARTMENT SUB: FEM

MODEL PAPER 4
PART A (25 MARKS)
1. a. Give the limitations of FEM. (2M)
b. Briefly explain the role of shape function FEM analysis (3M)
c. Write a short note on numerical integration. (2M)
d. What do you understand by discritisation of the domain (3M)
e. What is Jocobian matrix (2M)
f. Determine the matrix relating strains and nodal displacements for an axisymmetric
triangular element. (3M)
g. What is difference between CST and LST (2M)
h. What are the ways by which a 3D problem can be reduced to a 2D problem (3M)
i. Name Few FEA packages (2M)
j. Derive the convection matrix for a 1D linear bar element (3M)

PART B (10X5=50)

2. a) Write down expressions for the element stiffness matrices and element body force vectors

b) Evaluate the stresses in each element Determine the reaction force at the support. consider
1in=1cm for SI UNITS

OR

3.a) Determine the nodal displacement, Element stresses f r axially loaded bar as shown in the fig.
below
b) Explain the elimination method and penalty method for imposing specified displacement
boundary conditions

4. a) Obtain the forces in the plane Truss shown in Figure below and determine the support reactions
also. Take E=200GPa and A= 2000mm2

b) Derive the Hermite shape functions for beam element.

OR

5.a) Analyze the beam shown in Figure method and determine the end reactions. Also determine the
below by finite element deflections at mid spans given E=2X105N/mm2, and I=5X106 mm4

b) What are the general features of a bar Element?

6. a) Formulate the finite element equations for Constant strain triangle shown in fig. Assume plane
stress E=200Gpa,ν=0.25,thickness=5mm,nodal co-ordinates. Pressure on 1-2 edge is 5N/mm2

X1=1 X2=5 X3=3

Y1=2 Y2=4 Y3=6
b) Write the Advantages of iso-parametric elements

OR
7.a) For the configuration shown in figure, determine the deflection at the point load application using

a one element model. T = 10 mm , E = 70 G Pa ,ν = .3

b) Derive the strain displacement matrix for triangular element.

8. a) The plane wall shown in fig. The thermal conductivity K = 25 /m0C and there is a uniform
generation of heat in the wall Q = 400 W /m3. Determine the temperature distribution at five nodes
(include two sides of the walls) in equal distances through the wall thickness

b) Derive Approximate the first two natural frequencies of a cantilever beam using one element model.
EI=Flexural rigidity

9. a) A metallic fin with thermal conductivity K=360W/m0c, 1mm thick and100mm long extends
from a plane wall whose temperature is 2350c. Determine the distribution and amount of heat
transferred from the fin to air at 200c with h= 9W/m20c take width of the fin is 1000 mm.
Assume tip is insulted

b) Explain the concept of numerical integration and its utility in generating Isoperimetric finite
element matrices
10.a) Determine the strain displacement matrix for the TETRAHEDRAL element as shown in fig

.
b) Determine the approximate first two natural frequeneves of a simply supported beam using on a element. Flexural
Rigidity =EI; Density = Ƥ Cross-sectional area=A

OR

11.a) State the method used for obtaining natural frequencies and corresponding eigen vectors.

b) Evaluate natural frequencies for the CANTI LEVER beam shown in fig USING ONE ELEMENT

L
 MALLAREDDY COLLEGE OF ENGINEERING
AND TECHNOLOGY
MECHANICAL DEPARTMENT SUB: FEM

MODEL PAPER 5
PART A (25 MARKS)
1. a. What is the principle of FEM (2M)
b. Write the stress strain relations for 2D plane stress and pane strain conditions (3M)
c. Differentiate between truss and beam element based on degree of freedom. (2M)
d. What is Hermite shape function (3M)
e. Write the formula for the load vector of triangular element subjected to body force(2M)
f. What is the size of stiffness matrix for axisymmetric triangular element (3M)
g. What is the degree of freedom for the thermal problems (2M)
h. Where do you apply (3M)
i. Name Few FEA packages (2M)
j. Explain the importance of lumped mass matrix (3M)

PART B 10X5=50

2. a)Why polynomial type of interpolation function is preferred over trigonometric functions? Explain

b)Draw the Pascal's triangle and Pascal's tetrahedron for understanding the interpolations functions.
Explain the salient features

OR

3. a)Explain the steps involved in obtaining an appropriate solution using AEM

b) Explain the equilibrium state of the system, when the system is subjected to different types of loads
and explain the stress and equilibrium relations

4.For a two dimensional structure as shown in figure. determine displacement of the nodes and normal
stresses developed in the members using FE. Use E = 39x 106 N/cm2 and the diameter of the cross-
section of 0.25 cm.
 OR

5. A beam is fixed at one end and supported by a roller at the other end, has a 20 KN concentrated
load applied at the center of the span of 10 m. Calculate the deflection and slope and also construct
shear force and bending moment diagrams and take I=2500 cm4

6.a) Evaluate the axisymmetric stiffness matrix K of the triangular element shown in the figure
.Consider the coordinates of the nodes (2,1), (4,00, and (3,2).also assume E= 2.6 GPa and v=0.2

r3 3

r1 1

r2 2

b) Difference between CST and LST with respect to the triangular element.

OR

7.Derive the stiffness matrix for the four noded quadrilateral element in terms of natural coordinate
systems

8.consider a brick wall of thickness 0.3 m ,k=0.7 W/m/K. The inner surface is at 280C and the outer
surface is exposed to cold air at -150C.The heat transfer coefficient associated with the outside surface
is 40W/m2K.Determine the steady state temperature distribution within in the wall and also the heat
flux through the wall. Use two elements and obtain the solution

OR

9.Derive the conductivity matrix for two dimensional triangular element subjected to convection on
one face of the element

10. For the stepped bar shown in figure. Develop the global stiffness and mass matrices and also
determine the natural frequencies and mode shapes. Assume E=200 GPa and mass desity =7850
Kg/m3 L1=L2=0.3 m A1=350mm2 A2=600 mm2

L1 L2

OR

11.a)Derive the shape functions for the four noded tetrahedron element from the first principle
b)discuss the importance of semi automatic meshing and auto mesh along with the practical
applications