Indian Institute of Technology Bhubaneswar

School of Mechanical Sciences

End Semester (Spring) Tixamination - 2024

Subject Name: Finite Elenments Methods in Engineering Subject Code: ME6LO52

Date: 26.04.2024
Duration:3 ILours Full Marks: 100
Answer all questions.
Answer a new question on a fresh page (strictly) and all parts of a question together.

1. (a) Derive the strain measures in the Eulerian and Lagrangian presentations considering
the displacement field of the body. [8]
(b) Derive thie relationshiy among different quantities mentioned below in the present and
reference coordinate systems: [8]
) Volume
(ü) Area
(im) Cauchy Stresses, 1t Piola Krichoff Stress, 2nd Piola Krichof Stress.
() Why linear strain tensor is not an accurate strain measure as compared to Green
Lagrange sirain tensor? Explain with an example.

2. (a) Why Cauchy Stress tensor and 2nd Piola Krichoff Stress tensors are symmetric while
the l PiolaKrichoff Stress tensor is not symmetric. 5]
(b) How rany independent material properties are essential for linear isotropic,
orthottopic, and anisotropic materials? Explain why? 5]
(c) What is he difference between plane stress and plane strain condition? Explain with
examples. Also write down the equilibrium equation, compatibility equations and
constitutive reiation of a linear isotropic body in plane stress and plane strain condition.
Derive the Beitrami Michell equations of a body in plane stress and plane strain condition.
[16]
d Write down the strain compatibility equations for a linear isotropic body and explain
its significance. [4]
e (4) What is Airy's stress function? How does it help to solve problems insolid mechanics?
How to use Airy's stress function to obtain stresses in Cartesian and Polar coordinate
systems.

PT.O.
t (eHow to decide the order of Airy's stress function in solving solid mechanics problems.
Explain this with the example of asimpiy supported beam subjected to alinearly varying
[6]
transverse load.
3. (a) Athin plate of unit thickness subjected to different loads as shown in Fig.1. Find out
the distribution of the components of stress, strain, and displacements within the plate.
Make suitable (20]
assumptionsif any.

F1

2c X1

2a

Fig.1: Plate subjected to constant loads along its edges

(b) A thin plate with a hole of radius "r subjected to different loads as shown in Fig.2.
Also, FE a>> r and 2c >>r. Find out the distribution of the components of the stress
within the plate. The angle between the coordinate axis x) and the edge of the plate is 30°.
Make suitable assumptions if any. [20]

20
2a

Fig.2. Thin plaie with a central hole of diameter 'r subjected to loadings at its edges
 Indian Institute of Technology Bhubaneswar
School of Mechanical Sciences
Mechanical Systems Simulation Lab II (MEGPO03)
Mid-Semester Examination

MSSL-II Spring 2023-2024 Max mark: 30

1. Asimply supported beam is loaded as shown in Fig. . Using singularity functions. the
displacement along the beam can be expressed by the equation
15
uy (z) = -0)-(-5)]+ 6
(-8) +75 (z -7) +-238.25z
By definition, the singularity function can be expressed as,

(r-a)" =
(a-a)" if r>
if z<0

20

150

22

N/m; moment is 150 Nm; point

Figure l: Simply supported beam. Distributed force is 20
load 15 N

displacement versus distance along the beam,

Develop an M-file that creates a plot of
the beam.
T. Note that z =0at the left end of
jumper as a function of time is
2. The velocity of the free falling bungee
gm
tanh
(t)=ca
jumper after time t is given by,
the vertical distance travelled by the
 parameters in the equations are g = 9.81 m/s Tm 65 ko CA
Assume the any naner
kg/m. Find the distance travelled by the jumper after 3 8 using
0.25
method. writ» your own
function for the method and compare with result of the
cal
in-build Matlab function.
this problem is z() = then fcosh ( / t ) }. should Plot the numerica
Exact solutionfor top of exact solution. plot be labelled,
solution at t =1, 2, 3, 4, 5, 6, 7, 8 s on
appropriately.
titled etc. give the legends
(along x-direction) and 6 m
Temperature of the heaated rectangular plale of 8 mlong
3. y-direction) is
wide (along
T(, y) = 2ry +2r--2y+72.
temperature.
compute the average
structure under
spring-mass model is used to gain insight into the dynamics of the
4. A the force balance equations are,
earthauake. For the simplest model shown in Fig. 2
k
dt2
k3
-(T1 T2) +m2
dt2
k3
dt? -(2
dynamics of the structure from t=0 to 20 s, given the initial velocity
Simulate the and
floor is = l m/s and all other initial values of the displacement
of the ground (b)
time series plots of (a) displacement
dt
Present your result in two
velocities are zero.
velocities.

m, = 8,000 kg

|k, = 1800 KN/m

m, = 10,000 kg
k, = 2400 KN/m

m, = 12,000 kg
k, = 3000KN/m

structure under earthquake

Figure 2: Spring-mass model for

Page 2
 Indian Institute of Technology Bhubaneswar
School of Mechanical Sciences

Mid-Spring Semester Examination - 2024

Subject Name: Experimental Modal Analysis Code: ME6LO58

Date: February 19, 2024 (Monday) Semester: gh Semester(DD) +2 Semester M.Tech

Duration: 2 hours (10:00 AM - 12:00 Noon) Full Marks: 30+3 (bonus)

*Instructions:*
This question paper contains six questions in twO pages. Figures in the right hand margin indicate
marks. Read the question carefully before starts answering it.
All questions are to be answered. Try to answer a new question on a fresh page (strictly) and
all parts of a question together. Draw all the figures as clearly as possible.
Write all the steps, units clearly and give explanation for complete credit.
Assumptions made if any is valid till taking exam.

Problem 1

The modal model for a three DoF system is as given below. The mode shape is mass normalized
[21.0 [0.749 0.638 0.208
96.6 =0.478 -0.384 -0.534
0 2124! Lo. 223 -0.431 0.514

Two harmonic forces of 250 N and 300 Namplitude and of frequency 3.0 and 5.0 rad/sec
respectively are acting at 1* and 2d DoF respectively. Find the acceleration response at 3r
DoF. [7.5]

Problem 2

Two DoF system shown below. Determine a1 or az2 by direct as well as indirect method.
f() f2(t)
k k, K3

m W
x07
Given. m, = 5 kg m, = 10 kg kË = kz = 2 N/m kz = 4 N/m [7.5]
1
Problem3
analysis tocalculate the solutionof
Use modal

[6]
(0 1 mm and X(0) = (0 0)
for the initial conditions X(0) =

Problen 4
viscous damped MDOF system. (4]
Derive the expression for FRF matrix for a proportionally

Problen 5
answer. [3]
dynamic stiffness matrix? Can it be measured? Justify your (3]
a) What is modal analysis? How do they differ?
experimental
b) What is analytical and
Problem 6 (Bonus)
given in the figure below in log-log scale. Using
Receptance function of a SDOF system is stiffness and loss factor of
the
figure, determine the mass,
the information given in (3]
system.

10

10
[m/N]
Receptance

105

102 109
106
101 100 10
Frequency (rad/s)

2
 Indian Instituteof Technology Bhubaneswar
School of Mechanical Sciences
Mid Semester (spring) Examination -2024
Subject Name: Applied Elasticity Subject Code: ME6LO52
Date: 21.02.2024
Duration: 2 Hours Full Marks: 30

Answer all questions

Answer a new question on a fresh page (strictly) and all parts of a question together

1, (a) What is stress and how it is different from pressure? (10)

(b) What is the difference between stress vector and stress tensor?
(c) What are strain measure and strain displacement relations?
(d) Comments on similarities among the Eigen values and Eigen vector of the deformation
gradient tensor (E), Green Lagrange strain tensor (E), and Right Cauchy Green tensor C)
considering all the components are real.
(e) Write down all the field equations of linear elasticity.
2 An elastic body under the action of load undergoes the following motion: (8)
X1==5X-0.2X,-0.2X3, X;= -0.2X;+5X,-0.2X3, X3 =-0.2X1-0.2X_+5X3
reference
where x and X (i =1,2,3) denote the position of apoint in the present and
configuration, respectively. The units are in meters and seconds.
Lagrange strain
(a) Evaluate the deformation gradient (F), displacement gradient (D), Green
(E) and Linear strain tensor (E) at that point.
(dX), in that elastic body when
(b) Find out the deformed shape (dz) to the line segments
(i) dX =1e, (ii) dX =1e1+1e2, (ii) dX =1e,+1e,.les.
= (12)
3 Considera prismatic bar of square cross section (4m) and length L
this load, the
10m subjected to a load of magnitude f. Because of
Considering
change in the length of the bar at its free end is a =0.01m. Im

young's
the material is made of isotropic homogeneous material with L
0.3. Find out the
modulus E =200 GPA and Poisson's ratio V=
(E) and
deformation gradient tensor (E), Green Lagrange strain tensor
(placed at the centre of the
the stretch of the diameter of a circle
and e; direction.
surface of the bar) as shown in the figure along the e1
 Indian Institute of Technology Bhubaneswar
School of Mechanical Sciences
Mechanical Systems Simulation Lab II(ME6P003)
Semester Examination

MSSL-II Spring 2023-2024 Max mark: 30

1. Consider the spring system composed of three springs as shown in Fig. K, = 40 N/mm, K, =
60 N/mm, K, = 80 N/mm K, = 50 N/mm,the applied loads are F = 100 N,F, =
60 N. Determine the following using a MATLAB script,

" the global stifnesS matrix for the system.

" the displacements at nodes.
" the reactions at nodes.
" the deflection in each spring.
2. Two copper wires (1 and 3) and a steel wire (2) are connected by a body as shown in
the figure. The three wires are equal length of 50 cm and in the same cross-sectional
area of A = 65 mm'. When a load F = 15000 N is applied, find using a MATLAB
script the stresses in the copper and steel wires using three bar elements. For the
copper wire, use the elastic modulus E. = 80 Gpa. For the steel wire, use the elastic
modulus E, = 200GPa.

2100Pa
WaS a

Raen
5m

(a) Question 2. (b) Question 3

3. Determine the location and magnitude of Maximum Von-Mises stress in a plate as

shown in figure. Also find the stress concentration factor.
 Indian Institute of Technology Bhubaneswar
School of Mechanical Seiences
M.Tech. End Semester (Spring) Examination - 2024
Subject Name: Soft Computing and Applications Code: ME 6L060
Date: 29-04-2024
Sermester: Two
Duration: 3 hour Full Marks: 100
*Instructions:
All questions are to be answered. Try to answer a new question on a fresh page (strictly) and all
parts of a
question together. Draw all the figures as clearly as possible. Assune suitable data, if not given.

1. () Define the following terms related to optimization by taking a suitable example: (a) Decision variable (b)
Objective function (c) functional constraints (d) geometric constraints. [8]
(i) Minimize f (xl,x2) = 4x + 3x -6x] X2 - 4x1 [71
Subject to -10.0 < X1, X2 S 10.0
Using random walk method. Assume à = 1.2, e - 0.075, N = 100.

2. () Wiite nTtes on various crossover schemes used in real-coded genetic algorithm.

(ü) Write short notes on parametric study in GA. Howthe optimal parameters of GA is arrived. [5]

3. () Write the steps involved in the operation of NSGA.

(ii) Perforr: one iteration for the muiti-objective optimization problem using NSGA-II for the problem formulated
as below: Min: f; X1, f; =(1+ X)X, for the narent (Pt) and chile (Ot) nonulation given beiow. [10]
P
Qt
Soln XË X2 X1 X2
Soln
0.31 0.09
1 0.21 0.24
0.41 1.92
2 0.79 2.14
0.22 0.56 2.32
3 0.51
0.59 3.63 0.27 0.87
4
d
0.66 1.41 0.58 1.62
5
0.03 2.51 1.05
6 0.24

4 Identify the optimal parameters of the manufacturing process responsible for optimizing the hardness and surface
roughness of the component using Grey Relational Analysis. [10]

Trail Speed (RPM) Feed (mm/rev) DOC(mm) Hardness Roughness

288N 1.0 87.33 6.000
598 0.30 86 4.066
938 0.45 95.33 3.066
4 238 0.15 1.5 97.66 4.466
598 0.30 2.0 86.66 4100
6 938- 045 10 83 2.833
288 0.15 2.0 86.66 5.166
598 J 0.30 10 96.33 3.400
938 0.45 15 8L.66 3.366
 Indian Institute of Technology Bhubaneswar
Mechanical Sciences
School of
Spring End-Semester Examination-2024
Subject Code: ME6LAS1
Control of Mechanical Systems
Subject Name: Dynamies and

Date: 2nd May, 2024 FullMarks: 100

Duration: Three Hours
and all parts ofa question together.
Answer a new question on a fresh page (strictly)
Answer all questions.
feedback form and find Y/X. (10)
Figure 1 to unity
1. Reduce the block diagram shown in
y)

Figure 1 Figure 2

Figure 2 shows a simple vibration absorber model. Assuming the harmonic force F(t) = Asin(ot) is the
2.
equations for the system. Clearly mention the
disturbance applied to the mass M, derive tie state space
(12)
choice of state variables.

3. The differential equation of a linear system is given as

d'y(0) d'y(t) dy (0)
+5 dt2 .+6 + 10y(t) =r(t)
dt3 dt

where y(t) is the output, and r(t) is the input.

function, i.e., r(t) = u(t).
(a) Find the output y(t) for t > 0, when the input is unit step
(14)
(b) Obtain the final value of y (t) by using the final-value theorem.
 A The hlock diagram of acontrol system is
shown in Figure3. The error signalis defined to be e(t). and
100
G,()=
(1+0.1s) (1 +0.5s)
Find the step (K,), ramp (K,), and
parabolic (K) error constants in terms of Kand K. What
must be made, if any,on the values of Kand K,, so constraints
that the answers are valid (from the system
of view)? stability point
(14)
R(s) E)
Ys)
20:

K,

Figure 3
5. For the unity feedback system with
forward path transfer function G(S) = K(s+2)(s+1) (a) sketch the
(s2-2s+2)
complete root locus, (b) find the range of Kto keep the system stable, (c)
obtain the value of Kthat yíelds a
stable system with critically damped second-order poles (using root
locus diagram). (18)
6. Sketch the asymptotic Bode plots for the open loop transfer
function G(s)H (s) =TO in the semi
s(s+2) (s+4)
log graph papers. Deteriaine the values ef gain-margin and phase
margin from these plots. (18)
Draw the Nyquist plot fe thesystem shown in Figure 4. Using the
Nyquist stability criterion, find the range
of K for stability of this system.
(14)
R(S) + K C(s)
(s+2)

(s + 4)(s+6)

Figure 4
 Class Test 2
Dynamics and Control of Mechanical Sysens (ME6LO51)
1. Use the Routh-Hurwitz method to find the range of pronortijonalcontroller gain Kfor
which the following system willbe stable. (8)
P controller ciant
R(s}
K

2. Plot the root locus for the following open loop transfer function
K(s + 2s + 2)
G(S)H(s)=
s(s +1)(s +3)
Also comment on the stability of the system for range of values of K, based on root
locus. (12)
 Class Test 1
Dynamics and Control of Mechanical Systems (ME6LOS1 )
1. A solid sphere of radius r and mass m is projected along a rough horizontal surface
with the initial velocities indicated. If the final velocity of the sphere is to be zero.
express (in terms of vo, r, and 4z), (a) required magnitude of a,. (b) time t, required
for the sphere to come to est. The radius of gyrationfor solid sphere is r. (12)

Vo

2. Two uniform rods AB and CE, cach of mass 30 kg and length 24 m, are welded to each
other at their midpoints. Knowing that this assembly has an angular velocity of constant
magnitude w = 12 rad/s, determine (a) the angular momentum Lp of the assembly
about D, (b) bearing reaction forces due to rotational unbalance. (18)

D
B

X
9m
3m F
9m
3m
 Indian Institute of Technology Bhubaneswar
Sehool of Mechanical Sciences
Spring Mid-Semester Examination-2024
Subject Name: Dynamics andControl of Mechanical Systems
Subject Code: ME6L051

Date: 27h February, 2024

Duration: Two Hours FullMarks: 60
Answer a new question on a fresh page (strictly) and all parts of a question together.
Answer all questions.
1. The 200 kg satellite, shown in Figure 1, has its centre of mass at point G. Its radii of gyration about Z, X.
y' axes are k,, = 300 mm, kyy = kyy = 500 mm, respectively. At the instant shown, the satellite rotates

about x, y, and z' axes with the angular velocity shown, and its centre of mass Ghas avelocity of v, =
(-250 i + 200j + 120 k) m/s (in xyz frame). Determine the angular momentum of the satellite about
point A at this instant. Also determine the kinetic energy of the satellite at this instant. (15)
o-2Sok t2oft12ay
w = 1250rad/s
1.30-m

vc= 100 km/h

o= 600 rad/s 80 m
300 rad/s
800 mm

Figure 1 Figure 2
2. The car travels at a constant speed of v = 100 km/h around the horizontal curve having a radius of 80 m,
as shown in Figure 2. If each wheel has a mass of 16 kg, a radius of gyrationk=300 mm about its spinning
axis, anda radiusof 400 mm, determine the difference between the normal forces of the rear wheels, caused
by the gyroscopic effect. The distance between the wheels is 1.30 m. (8)
3. Describe the Holonomic and Non-holonomic constraints with suitable examples.
(2+2)
4. Using Lagrange's's equation, derive the cquations of motion of the system given in Figure 3. The mass of the
sliding cart is M. The dashpot is connected at the mid-point (G) of the oscillating bar of mass mand length
b. Assume that the spring (of stiffiness k) and the dashpot (of damping constant C) deflect only horizontaly
and that the force F is always applied horizontally at free end offthe bar. Clearly mention the Generalised
coordinates chosen for solving the problem. (15)

k M

Figure 3

boundary conditions for the transverse

5. Using Hamilton's principle, derive the equation of motion and
discrete mass
vibration of a uniform Euler-Bernoulli beam shown in Figure 4. The beam is connected to a
dashpot of damping
m and a spring of stiffness k at x = 0, and is attached to another discrete mass m, and a
constant cat x = l. Assume the necessary properties of the bar as shown inFigure. (18)

Beam, EI

I=0

Figure 4