Equations
Equations
Equations
Equation book
2 Strength of Materials 12
3 Theory of Machines 21
4 Machine design 32
5 Production Engineering 43
6 IM & OR 75
7 Material Science 87
12 Thermodynamics 148
1
www.eggsam.com CONTENTS [2]
Important
This equation book is a compilation of equations I collected during my studies. The
aim of this book is to make revision easier. If you are studying the subjects for the first
time, then please do not use this book. This book is not for first time study.
This book is for both GATE and ESE (IES). It is useful for other Mechanical
Engineering Exams as well. I have tried to include as many topics as possible. Some of
those topics are not required for GATE. So please follow the syllabus of the exam you
are appearing for.
Projectile motion
v0 =Launching speed
θ= Launching angle w.r.t horizontal plane
3
www.eggsam.com 1.2. MOMENTUM [4]
Horizontal acceleration ax = 0
Vertical acceleration ay = −g
Horizontal velocity vx = v0 cos θ
Vertical velocity vy = v0 sin θ − gt
Horizontal displacement at t x = v0 t cos θ
1
Vertical displacement at t y = v0 t sin θ − gt2
2
2v0 sin θ
Time of flight t=
g
v02 sin2 θ
Maximum height h=
2g
v02 sin 2θ
Range R=
g
o
For maximum range θ = 45
gx2
Path y = x tan α − 2
2u cos2 α
Inclined plane
β: Inclination of plane
θ: launching angle wrt horizontal plane
2u sin(θ − β)
Time of flight t =
g cos β
2u2 sin(θ + β) cos α
Range R =
g cos2 β
Coefficient of restitution
Relative velocity after impact
e=−
Relative velocity before impact
e=0 for perfectly plastic collision
e=1 for perfectly elastic collision
Cylindrical coordinates
Position P~ = rr̂ + θθ̂
Velocity V~ = ṙr̂ + rθ̇θ̂
Acceleration ~a = (r̈ − rθ̇2 )r̂ + (rθ̈ + 2ṙθ̇)θ̂
3/2
[1 + (dy/dx)2 ] [1 + y 02 ]3/2
Radius of curvature r = =
d2 y/d2 x y 00
1.2 Momentum
Momentum P~ = mV~
-Vector
Conservation of momentum: In the absence of external forces, the total momentum
of a system is constant
www.eggsam.com 1.3. COLLISION [5]
1.3 Collision
Elastic collision
Linear momentum and kinetic energy are conserved.
m1 with velocity u1 collides with m 2 with velocity
u2
m1 − m2 2m1
Velocity after collision of m1 , v1 = u1 + u2
m1 + m2 m1 + m2
m2 − m1 2m2
Velocity after collision of m2 , v2 = u2 + u1
m1 + m2 m1 + m2
1.4 Force
Newton’s laws of motion
First law: In an inertial frame of reference, an object either remains at rest or continues
to move at a constant velocity, unless acted upon by a force.
Second law: In an inertial frame of reference, the vector sum of the forces F on an object
is equal to the mass m of that object multiplied by the acceleration a of the object.
F = ma.
Third law: When one body exerts a force on a second body, the second body simulta-
neously exerts a force equal in magnitude and opposite in direction on the first body.
d~v
Force F~ = m = m~a
dt
Torque ~τ = I~
α
mv 2 r̂
Centripetal force F~c = −
r
Force on a spring F = kx
x: extension/compression of the spring
k: spring constant or stiffness of the spring
Gm1 m2
Gravitational force F~ = − r̂
r2
G: 6.67408 × 10−11 m3 kg −1 s−2
r: Distance between the bodies
Lami’s Theorem
P Q R
= =
sin α sin β sin γ
Impulse
Integral of a force over the time interval for which it acts.
Change in momentum
www.eggsam.com 1.5. TRUSS [6]
1.5 Truss
Perfect truss, number of members m=2j-3
m>2j-3: Redundant truss
m<2j-3: Defficient truss
j: number of joints
m: number of members
1.6 Friction
For a stationary body, F~f riction = F
For a moving body, F~f riction = µmg
Angle of friction φ = tan−1 µ
Angle of Repose α = tan−1 µ
α=φ
F: horizontal force acting on the body.
Limiting friction: The maximum friction that can be generated between two static
surfaces in contact with each other. Once a force applied to the two surfaces exceeds the
limiting friction, motion will occur.
Stiction is the static friction that needs to be overcome to enable relative motion of
stationary objects in contact.
Amontons’ First Law: The force of friction is directly proportional to the applied
load.
Amontons’ Second Law: The force of friction is independent of the apparent area of
contact. [As long as there is contact]
Coulomb’s Law of Friction: Kinetic friction is independent of the sliding velocity.
1.8 Constants
π 3.14
Euler’s constant e 2.718
Plank’s constant h 6.625 × 10−34 Js
Universal gas constant 8.314 J/molK
Stefan-Boltzmann constant σ 5.67 × 10−8 Wm −2
K −4
Water
Density 1000 kg/m 3
Molecular mass 18.015 g/mol
Gas constant 461 J/kgK
Sp. heat of ice 2.108 kJ/kgK
Sp. heat of water 4.187 kJ/kgK
Sp. heat of water vapor 1.996 kJ/kgK
Latent heat of fusion 336k J/kg
Latent heat of vaporization 2264.7 kJ/kg
Air
Density 1.223 kg/m 3
Mean molecular mass 28.9 g/mol
Gas constant 287 J/KgK
Sp. heat at const volume 718 J/kgK
Sp. heat at const pressure 1005 J/kgK
1.9 Units
101 Deca (da) 10−1 deci (d)
102 Hecto (h) 10−2 centi (c)
103 kilo (k) 10−3 milli (m)
106 Mega (M) 10−6 micro (µ)
109 Giga (G) 10−9 nano (n)
1012 Tera (T) 10−12 pico (p)
1015 Peta (P) 10−15 femto (f)
1018 Exa (E) 10−18 atto (a)
1021 Zetta 10−21 zepto
1024 Yotta 10−24 yocto
7 Basic units
Quantity Name SI Unit
Mass kilo gram kg M1 L0 T0
Length meter m M0 L1 T0
Time second s M0 L0 T1
Electric current Ampere A
Temperature Kelvin K
Luminous intensity candela cd
Amount of substance mole mol
Derived units
Quantity Name Symbol In SI base units
plane angle radian rad (mm−1 )
solid angle steradian sr (m2 m −2 )
frequency hertz Hz s−1
force, weight newton N kgms−2
pressure, stress pascal Pa kgm−1 s−2
energy, work, heat joule J kgm2 s−2
power, radiant flux watt W kgm2 s−3
electric charge coulomb C As
voltage (electrical potential), emf volt V kgm2 s−3 A−1
capacitance farad F kg−1 m−2 s4 A2
resistance, impedance, reactance ohm Ω kgm2 s−3 A−2
electrical conductance siemens S kg−1 m−2 s3 A2
magnetic flux weber Wb kgm2 s−2 A−1
magnetic flux density tesla T kgs−2 A−1
inductance henry H kgm2 s−2 A−2
luminous flux lumen lm cd
illuminance lux lx m−2 cd
radioactivity (decays per unit time) becquerel Bq s−1
absorbed dose (of ionising radiation) gray Gy m2 s−2
equivalent dose (of ionising radiation) sievert Sv m2 s−2
catalytic activity katal kat mol s−1
area m2
volume m3
speed, velocity ms−1
acceleration ms−2
wavenumber m−1
density kgm−3
surface density kgm−2
specific volume m3 kg−1
current density Am−2
magnetic field strength Am−1
concentration molm−3
mass concentration kgm−3
luminance cdm−2
refractive index
relative permeability
www.eggsam.com 1.10. VECTOR [10]
1.10 Vector
P~ = p1 î + p2 ĵ + p3 k̂
Q~ = q1 î + q2 ĵ + q3 k̂
θ: Angle between the vectors
Vector ~ ~ ~
p sum R = P + Q = (p1 + q1 )î + (p2 + q2 )ĵ + (p3 + q3 )k̂
R = P 2 + Q2 + 2P Q cos θ
Dot product P~ · Q ~ = (p1 q1 ) + (p2 q2 ) + (p3 q3 ) = P Q cos θ
î ĵ k̂
Cross product P~ × Q ~ = det p1 p2 p3
q1 q2 q3
~ ~
|P × Q| = P Q sin θ
1.11 Misc
R
x̄dA
Centroid XG = R
dA
m1 m2
Reduced mass µ =
m1 + m2
Parallel axis theorem Id = ICM + md2
Perpendicular axis theorem Izz = Ixx + Iyy
Center of percussion
The point on an extended massive object attached to a pivot where a perpendicular im-
pact will produce no reactive shock at the pivot. Translational and rotational motions
cancel at the pivot when an impulsive blow is struck at the center of percussion.
The same point is called the center of oscillation for the object suspended from the pivot
as a pendulum, meaning that a simple pendulum with all its mass concentrated at that
point will have the same period of oscillation as the compound pendulum.
Varignon’s Theorem
If many concurrent forces are acting on a body, then the algebraic sum of torques of all
the forces about a point in the plane of the forces is equal to the torque of their resultant
about the same point.
Radial run-out is the result of a rotating component running off center, such as a
ball bearing with an offset center. This means that the rotating tool or shaft, instead of
being centrally aligned, will rotate about a secondary axis.
Simple pendulum
www.eggsam.com 1.11. MISC [11]
r
l
T = 2π
g
Compound
r pendulum
I
T = 2π
mgh
1 πd4 2d
Semicircle Ixx =
2 64 3π
2
Solid sphere I = mr2 Center
5
2
Hollow sphere I = mr2 Center
3
Slender rod about 1
I = ml2 Center
midpoint 12
1
Slender rod about end I = ml2 Center
3
Izz = mr2
Circular ring 1 2 Center
Ixx = Iyy = mr
2
1
Izz = mr2
Solid disc 2 Center
1 2
Ixx = Iyy = mr
4
Cylinder I = m(r12 + r22 ) Center
3R
Hemisphere
8
Chapter 2
Strength of Materials
12
www.eggsam.com 2.1. STRESS AND STRAIN [13]
Shear stress
F
τ=
A
= P/G
Volumetric stress
ev = ex + ey + ez
4V
ev =
V
v = P/K
1
= (σx + σy + σz ) (1 − 2µ)
E
Bar of uniform
h ρg istrength
Ax = A exp x
P
Tapered rod
PL
δl =
πd1 d2
E
4
Uniformly tapered rectangular bar with constant thickness t
PL a
δl = ln
Et(a − b) b
Rotating bar about one end
1
σx = ρω 2 [l2 − x2 ] (From axis)
2
1 23
δl = ρω l
3
Strain due to self weight
PL ρgL2
δl = =
2AE 2E
Thermal stress, at least one end free
σ=0
δl = Lα∆T
Impact load r
2h
Impact factor IF =1 + 1+
δstatic
σ = IF × σstatic
δ = IF × δstatic
www.eggsam.com 2.1. STRESS AND STRAIN [14]
PL
δstatic =
AE
h: Height from where the weight P is dropped.
P
Sudden load: σ = 2
A
Relation between elastic constants
E = 2G (1 + µ)
E = 3K (1 − 2µ)
9KG
E=
3K + G
Lame’s constant
Eµ
λ=
(1 + µ) (1 − 2µ)
Number of independent elastic constants
Isotropic : 2
Orthotropic : 9
Anisotropic : 21
* [0 ≤ x ≤ L/2]
Conjugate beam method
Beam Conjugate beam
Slope Shear force
Deflection Bending moment
End pin support Pin support
Internal Pin support Hinged joint
Hinged joint Internal pin support
Fixed end Free end
Free end Fixed end
Maxwell’s Reciprocal theorem
In any beam or truss, the deflection at any point D due to load W at any other point C
is the same as the deflection at C due to the same load at D.
2.5 Torsion
T τ Gθ
= =
J R L
J
Torsional section modulus Z =
R
T GJ
Torsional stiffness = =
θ L
Equivalent Moment and torsion
1 √
Me = M + M2 + T 2
2
www.eggsam.com 2.6. STRAIN ENERGY [18]
√
Te = M 2 + T 2
Pure torsion assumptions
Uniform material
Uniform twist
Shaft is of uniform circular cross section
Cross sections which are plane remains plane
All radii remain straight after twist
Compound shaft
Shafts in series Parallel shafts
θ = θ1 + θ2 θ = θ1 = θ2
T = T1 = T2 T = T1 + T2
Castigliano’s theorem
∂U
yr =
∂Pr
If the strain energy of a linearly elastic structure can be expressed as a function of gen-
eralised force Pr then the partial derivative of the strain energy with respect to Pr gives
the generalised displacement yr in the direction of Pr .
www.eggsam.com 2.7. APPLICATIONS [19]
2.7 Applications
Thin Cylinder
Pd
Hoop stress σh =
2t
Pd
Longitudinal stress σl =
4t
Pd
Longitudinal strain el = (1 − 2µ)
4tE
Pd
Hoop strain eh = (2 − µ)
4tE
Pd
Volumetric strain ev = (5 − 4µ)
4tE
Thin Sphere
Pd
Hoop stress σh =
4t
Pd
Hoop strain eh = (1 − µ)
4tE
Pd
Volumetric strain ev = 3 (1 − µ)
4tE
Thick Cylinder
B
Radial stress at radius r, σr = A − 2
r
B
Hoop stress at radius r, σh = A + 2
r
Thick Spherical shell
2B
Radial stress at radius r, σr = 3 − A
r
B
Hoop stress at radius r, σh = 3 + A
r
Springs
D
Spring index C =
d
Gd4
k=
64R3 n
WR
τmax =
πd3
16
64W R3 n
Maximum deflection δ =
Gd4
4C − 1 0.615
Wahl’s factor = +
4C − 4 C
Type le
Both ends fixed L/2√
One end fixed other end hinged L/ 2
Both ends hinged L
One end fixed, other end free 2L
Rankine’s formula
1 1 1
= +
PR PC PE
σc A
PR = 2
Le
1+a
k
σc
a = 2 : Rankine’s constant
π E
PC = Crushing load = σc A
PE = Euler’s load
PR = Rankine’s load
2.8 Misc
Axial rigidity = AE
I
Section modulus Z =
Y
J
Torsional section modulus Z =
R
GJ
Torsional stiffness=
l
le Effective length
Slenderness ratio λ = =
K Least radius of gyration
Beam
Transverse load
Bending stress on the beam
Column
Axial load
Vertical member
Compressive stress on the column
Fails by buckling if it is slender
Strut
Vertical, horizontal or inclined
Fails by buckling if it is slender
Chapter 3
Theory of Machines
A material body which is common to two or more kinematic pair is called a link
Fluid can also act as a link
Structure: DOF=0
Super structure/Preloaded structure: DOF<0
Mechanism: DOF≥1
Statically indeterminate system, DOF ≤ −1
21
www.eggsam.com 3.1. ANALYSIS OF PLANAR MECHANISMS [22]
Kinematic chain : Last link is connected to first link and no link is fixed
1 Vin Fout Torque output
Mechanical advantage = = = =
velocity ratio Vout Fin Torque input
At Toggle position : infinite mechanical advantage
4 bar mechanism
dµ ls sin θ
=
dθ pq sin µ
µ: Transmission angle
b : Coupler link
c : Output link
Grashoff ’s condition
l+s<p+q Grashoff’s I
l+s>p+q Non-Grashoff’s or Grashoff’s II
l+s=p+q Special Grashoff’s or Grashoff’s III
s: Shortest link
l: longest link
Grashoff ’s I
Shortest link fixed : crank-crank mechanism.
Adjacent link to shortest link fixed : crank-rocker mechanism.
Link opposite to shortest link fixed : rocker-rocker mechanism.
Coriolis Acceleration =2ωV
Instantaneous center
n(n − 1)
Number of instantaneous centers =
2
Where n is the number of links in relative motion.
Kennedy’s Theorem: When three rigid bodies are in relative motion, the three instan-
taneous centers shared by three bodies all lie on the same straight line.
Permanent instantaneous centers: will move, but always on a point on the link (Primary
instantaneous center)
Neither fixed nor permanent instantaneous centers: Also called secondary instantaneous
centers
Body centrode
Also called moving centrode
Locus of the instantaneous centre of the fixed body relative to the movable body.
Space centrode
Also called fixed centrode
Locus of the instantaneous centre of the moving body with respect to the fixed body.
3.1.1 Steering
w
Equation of correct steering cot φ − cot θ =
L
φ: Angle of outer wheel
θ: Angle of inner wheel
w: Pivot distance of the front stub axles
L: Wheel base
www.eggsam.com 3.1. ANALYSIS OF PLANAR MECHANISMS [24]
n = l/r
h i
p
Position x = r (1 − cos θ) + n − n2 − sin2 θ
" #
sin 2θ
Velocity v = rω sin θ + p
" 2 n2 − sin2 θ #
2
cos 2θ sin 2θ
Acceleration a = rω 2 cos θ + p − 3/2
n2 − sin2 θ 4 n2 − sin2 θ
ω cos θ
Angular velocity of connecting rod = p
n2 − sin2 θ
−ω 2 sin θ(n2 − 1)
Angular acceleration of connecting rod = 3/2
n2 − sin2 θ
When n >> 1
sin 2θ
Velocity v = rω sin θ +
2n
cos 2θ
Acceleration a = rω 2 cos θ +
n
ω
Angular velocity of connecting rod = cos θ
n 2
ω
Angular acceleration of connecting rod = − sin θ
n
Pressure force FP = Pressure × Area
FP
Crank effort FCR =
cos φ
Force on cylinder wall FW = FCR sin φ
Force on crank bearing FB = FCR cos(φ + θ)
Tangential Force on crank Ft = FCR sin(φ + θ)
Turning moment T = Ft × r
www.eggsam.com 3.2. GEARS [25]
Inversions
1st inversion Ground is fixed Slider crank mechanism
Whit-worth Quick return
mechanism
2nd inversion Crank is fixed Rotary Engine
Shaping machine
Planning machine
Crank and slotted lever
3rd inversion Rocker is fixed quick return mechanism,
Oscillating cylinder engine
Hand pump, Bull engine,
4th inversion Slider is fixed
pendulum pump
Quick return mechanism
time of cutting 180 + 2α
=
time of return 180 − 2α
Double slider crank mechanism
1st inversion Ground is fixed Elliptical trammel
2nd inversion One slider is fixed Scotch yoke mechanism, SHM
3rd inversion Connecting rod is fixed Oldham coupling
Klien’s construction
Angular acceleration of crank is zero
Green: Velocity
Red: Acceleration
3.2 Gears
Gears are Positive drive (No slip is possible)
(In negative drive, slip is possible)
T: Number of teeth on gear t: Number of teeth on pinion
D,d: Pitch diameter
Ra , ra : Addendum radius
R,r: Pitch radius
www.eggsam.com 3.2. GEARS [26]
φ: Pressure angle
Pitch circle diameter = d
Base circle diameter= d × cos φ (Constant for a gear)
d
module m =
T
πd
Circular pitch =
T
T
Diametral pitch =
d
Addendum = Ra - R
N1 ω1 T2
Gear ratio G= = =
N2 ω2 T1
input speed
Velocity ratio =
output speed
p
Gear path of approach = pRa2 − R2 cos2 φ − R sin φ
Gear path of recess = pra2 − r2 cos2 φ − r p
sin φ
Gear path of contact = Ra − R cos φ + ra2 − r2 cos2 φ − (R + r) sin φ
2 2 2
Path of contact
Arc of contact =
cos φ
Lenght of action arc of contact Angle of action
Contact ratio = = =
mπ cos φ πm Pitch angle
Working depth = Addendum + Dedendum - clearance
Working depth = Sum of addendum of both the gears
To avoid interference
2AG T
T ≥s G=
1 1
t
1+ + 2 sin2 φ − 1
G G
2AP
t≥ p
1 +pG (G + 2) sin2 φ − 1
Ra,max = (R cos φ)2 + (R sin φ + r sin φ)2
In rack and pinion,
2AR
T ≥
sin2 φ
To reduce interference
-add undercutting
-increase pressure angle
-tooth stubbing
-increase number of teeth in gear
-slightly increase centre-centre distance
-decrease gear addendum and increase pinion addendum
In an epicyclic
gear
train
T1 N2 − Na
=−
T2 N1 − Na
Helical or spiral gears
ψ: Spiral angle or helix angle
φ: Pressure angle
θ: angle between axes
www.eggsam.com 3.3. FLY WHEELS [27]
3.4 Governors
Centrifugal Governor
1
h∝ 2
N
N2 − N1
Sensitiveness=
N
Isochronous governor: Speed is constant
Effort of governor: Mean force acting on the governor
Hunting: Speed of engine fluctuates above and below the mean speed.
rω 2
Gravity controlled tan θ =
Watt Governor g
Pendulum type g 895
h= 2 = 2
ω N
Gravity controlled Frictional force f is also present at the
Watt Governor + Mass sleeve
2
2mg + (M g ± f )(1 + tan β
tan θ
)
Porter Governor (M) at the sleeve. ω =
Porter governor can not be 2mh
895 Mg ± f tan β
isochronous. h= 2 1+ (1 + tan θ )
N 2mg
Gravity controlled 895 a Mg ± f tan β
Proell Governor h= 2 1+ (1 + tan θ )
Porter governor + arms N e 2mg
Hartnell Governor Spring controlled
Spring controlled
Mg
Hartung Governor Hartnell + change in spring mrω 2 a = k(r − r0 )a + b
position 2
Hartnell governor with
Wilson Hartnell
extra spring
3 leaf springs
Pickering Governor
Used in gramophone
3.5 Balancing
Static balancing
Center of mass of the system must lie on the rotational axis
Σmr sin θ
tan θc =
pΣmr cos θ
mc rc = (Σmr sin θ)2 + (Σmr cos θ)2
Dynamic Balancing
-If there is a couple
Σmrl + mc rc lc = 0
Hammer blow: Maximum vertical unbalanced force by the mass used to balance the
reciprocating masses. = cmrω 2 p
Resultant primary unbalanced force = ((1 − c)mrω 2 cos θ)2 + (cmrω 2 sinθ)2
mrω 2 cos 2θ
Secondary unbalanced force =
n
SHM
h πω
Vmax =
2 φa 2
h πω
Acceleration fmax =
2 φa
φa : Angle of ascent
Uniform velocity
ω
V =h
φa
Constant Acceleration
ω
Vθ = 4h 2 θ
φ 2
ω
a = 4h
φ
3.7 Gyroscope
Gyroscopic torque = Iω × ωp
Active gyroscopic couple: Rotate the spin vector through 90o about precession axis to
www.eggsam.com 3.8. VIBRATION [30]
3.8 Vibration
Reileigh’s
r method
g
ωn =
deflection
Energy method
d(KE + P E)
=0
dt
Lagrange’s equation
L=KE-PE
d ∂L ∂L
− =0
dt ∂ q̇i ∂qi
Damped vibrations
mẍ +crẋ + kx = 0
k
ωn =
m r
c c 2 k
s1,2 = − ± −
2mh p 2m i m
s1,2 = ωn −ξ ± ξ 2 − 1
x = Xe−ξωn t [sin(ωd t + φ)]
√
Critical damping Cc = 2 mk = 2mωn
C
Damping ratio/factor ξ =
Cc
Degree of damping = ξ 2 p
Damped frequency ωd = ωn 1 − ξ 2
2π 1
Damped time period Td = p
ωn 1 − ξ 2
x1 2πξ
Logarithmic decrement δ = ln = ξωn Td = p
x2 1 − ξ2
Forced oscillation
(F0 /k)
Steady state Amplitude X = v(
u 2 )2 2
t 1− ω ω
u
+ 2ξ
ωn ωn
www.eggsam.com 3.8. VIBRATION [31]
ω
−2ξ
−cω ωn
tan φ = 2
= 2
k − mω ω
1−
ωn
X 1
Magnification factor = v(
Xst u 2 )2 2
t 1− ω ω
u
+ 2ξ
ωn ωn
Transmissivity ratio
s
2
ω
1 + 2ξ
FT ωn
TR = = v(
F u 2 2
) 2
t 1− ω ω
u
+ 2ξ
ωn ωn
Under damped,
ω √
TR increases if < 2
ωn
ω √
TR decreases if > 2
ωn
Force transmitted =F × T R
Vibration isolation
0 < TR√< 1
ω
> 2
ωn
Rotating unbalanced
mass
2
mr ω
M ωn
X = v(
u 2 )2 2
t 1− ω ω
u
+ 2ξ
ωn ωn
Spring,
r rolling disk, horizontal
2k
ωn =
3m
Whirling of shafts
2
ω
e
ωn
r = v(
u 2 )2 2
t 1− ω ω
u
+ 2ξ
ωn ωn
Dunkerly’s method for multiple rotor system
1 1 1
= + + ...
ωn2 ω1 ω22
2
Torsional
r vibration r
Torsional stiffness GJ
ω= =
I Il
Chapter 4
Machine design
32
www.eggsam.com 4.2. FATIGUE FAILURE [33]
2a
Maximum stress σm = σB 1+
b
1
Corrected endurance limit σe0 = σe Ka Kb Kc Kd Ke
Kf
Name Equations
Ka = 1.0 d ≤ 7.6mm
Size factor Ka = 0.85 7.6 ≤ d ≤ 50mm
Ka = 0.75 d ≥ 50mm
Kb = 1.0 For reversed bending load
Loading factor Kb = 0.8 Reversed axial load for steel
Kb = 0.577 Reversed torsional load
Surface factor Kc = 1, for polished surface
Temperature factor Kd
Ke = 1.0 50% Reliability
Reliability factor Ke = 0.89 90% Reliability
Ke = 0.7 99.99% Reliability
Fatigue strength Endurance limit for notch free speciment
Kf =
reduction factor Endurance limit for notched specimen
Kf − 1
q=
Kt − 1
max stress
Stress concentration factor Kt =
avg stress
1 √
Notch sensitivity q = , a = Neuber’s constant
1 + (a/r)2
Endurance limit σe = 0.5σut for steels
Endurance limit σe = 0.4σut for Cast Iron, Cast steels, Al alloys
Endurance limit σe = 0.3σut for Cast Al alloys
n1 n2 n3
Miners approach + + + ... = 1
L1 L2 L3
www.eggsam.com 4.3. KEYS [34]
4.3 Keys
b: width of the key
t: thickness
l: length
Usually b=D/4
Shear design
F
τ=
bl
Crushing design of key
F
σ=
lt/2
Loose fit
Feather Key
permits axial movement
Self aligning
Woodruf Key Used in tapered
curved
Square/Rectangular Normal
Taper sunk key Tapered
Barth key Edges adjusted
Kennedy Key Diamond
two keys
Tangent Keys One key can transmit
power in one dir
4.4 Rivets
√
Unwin’s formula d = 6.04 t
Bearing/crushing failure Fb = σb dt
Least Strength
Efficiency η =
Strength of plate without rivet
Rivet value R=min(Shear strength, Bearing strength, plate strength,... )
No of rivets required=F/R
www.eggsam.com 4.5. THREADED JOINTS [35]
b−d
For diamond riveting, η =
d
Eccentric loading
F A1
F1 =
ΣAi
0 F eA1 r1
F1 =
ΣAi ri2 p
Resultant force F = F12 + F102 + 2F1 F10 cos θ
F1 : Primary shear stress on rivet 1
F10 : Secondary shear stress on rivet 1
F: Force applied
ri : Distance between ith rivet and CG
Ai : Area of ith rivet
e: Distance between CG of rivets and load
4Fi
Crushing stress σc =
π(d2 − d2c )n
t: pitch
F
τmax =
0.828hl
F
For both parallel and transverse welds use the formula
0.707hl
Unsymmetrically welded and Axially loaded
F = F1 + F2
Fb
F1 =
a+b
Fa
F2 =
a+b
Eccentric loading
F
τ1 =
Throat Area
My
σbending =
I
bt3 btd2 d2
I= + ≈ bt
12 4 4
Circular weld Subjected to torsion
Mr
τ=
J
J = 2πtr3
Petroff ’s equation
r µN
2 s
Coefficient of friction, f = 2π
c P
Ns in rps
Bearing load W = Bearing pressure ×LD
Power loss = f W rω
µN
Bearing characteristic number BCN =
P
N: rpm
min(BCN)=k=Bearing modulus
Sommerfield number
-decides
r film
thickness
2 µNs
S=
c P
Ns in rps
www.eggsam.com 4.8. ROLLING CONTACT BEARING [38]
F L1/n = C
n=3 for ball bearing
n=3.33 for roller bearing
C= Dynamic load rating
Stribeck’s Equation
-gives static load capacity of the bearing
Kd2 Z
C0 =
5
Z: number of balls
www.eggsam.com 4.9. CLUTCH [39]
d: ball diameter
K: Constant
Ball bearings Light loads
Anti-friction bearings More accurate centering
Very high radial load
Deep groove bearing
Both radial and axial load
Thrust bearing Only Axial load
Oscillatory motion.
Can be mounted directly on shaft
Needle roller bearing
(piston pin bearing, rocker arms,
universal joint)
Double row roller bearing
Can carry radial and axial loads
Tapered roller bearing Should be preloaded
4.9 Clutch
Uniform pressure
3 (new clutch)
2 r1 − r23
T = µW
3 r12 − r22
W = P π (r12 − r22 )
2 r13 − r23
Friction radius =
3 r12 − r22
Uniform wear (old clutches)
Pr = C
r1 + r2
T = µW
2
r1 + r2
Friction radius =
2
W = 2πC(r1 − r2 )
Centrifugal clutch
T = µmrg rb (ω22 − ω12 )× Number of shoes
m: Mass of one shoe
Single-plate clutch : Occupies large space (Trucks)
Multiple clutch : Occupies less space (Scooter)
Centrifugal clutch : Starts after reaching a critical velocity (mopeds)
Jaw clutch : Low speeds only (Rolling mills)
It is more logical and safer to use uniform wear theory in the design of clutches.
www.eggsam.com 4.10. BRAKES [40]
4.10 Brakes
4.11 Gear
Working depth = sum of addendum of gear and pinion
d
module m=
T
πd
Circular pitch=
T
T
Diametral pitch =
d
Addendum = m
Dedendum =1.157m
Ft = bmσb Y (bombaY)
t2
Lewis form factor Y=
6hm
My Ft ht/2
σb = = 3
I bt /12
Max torque
Service factor Cs =
Rated torque
3
Velocity factor kv =
3+v
www.eggsam.com 4.12. POWER SCREW [41]
-surface
-fatigue failure
Pitting
-small cavities
-near pitch surface
Scoring/Scuffing Insufficient lubrication + metal-metal contact
Plastic flow Yielding of surface under heavy loads
Abrasion Erosion due to foreign abrasive particles
For Lifting
µ cos α + sin α
Force F = W
cos α − µ sin α
Torque T = W tan(φ + α) × D/2
For Lowering
µ cos α − sin α
Force F = W
cos α + µ sin α
Torque T = W tan(φ − α) × D/2
Self locking
tan φ ≥ tan α
WL tan α
Efficiency η = =
2πT tan(φ + α)
www.eggsam.com 4.13. MISC [42]
π φ
For maximum efficiency α = −
4 2
1 − sin φ
ηmax =
1 + sin φ
ACME threads
µ
µ0 =
cos β
4.13 Misc
To connect two joints in
Cotter Joint
tension/compression
Connects two joints and allows
Knuckle joint
angular movement at the joint
Join two rods having threads with the
Turn buckle
help of coupler nut
Riveted joints Make leak proof joints. (Boilers)
Knuckle joint
-used to transmit axial tensile force.
-unsuitable to connect two rotating shafts, which transmit torque.
flexible shafts
-low rigidity in bending.
-high rigidity in torsion.
Chapter 5
Production Engineering
[The short notes are provided for easy revision. Do not use this for first time study.]
5.1 Metrology
Active inspection or online inspection: Checking when the product is being produced
Passive inspection: Inspecting already produced products
Accuracy: Correctness
Precision: Repeatability
i = 0.45D1/3 + 0.001D µm
D in mm
1
Gauge tolerance = Work tolerance
10
1
Wear allowances = Gauge tolerance
10
Workshop gauges: never sell anything bad
Inspection gauge: Never reject anything good
ISO gauges: Go gauge is from work shop gauge and No-go gauge is from Inspection gauge
General purpose gauge
43
www.eggsam.com 5.2. CASTING [44]
Interferometry
To measure flatness
5.2 Casting
Mould box
Cope
Cheek
Drag
Pouring basin → Sprue → Splash core → Runner (Trapeziodal) → Skimbob → Ingate
→ mould cavity, Riser
Advantages of casting
-Complex shapes can be made
-Cheaper
-Both ductile and brittle can be produced
-Large sized parts can be made
Disadvantages
-Surface is not smooth
-Time consuming process
-Non-uniform properties due to non-uniform cooling
www.eggsam.com 5.2. CASTING [45]
5.2.1 Allowances
Machining/finish allowance Allowance for finishing operations
Daft/Taper allowance For easy removal of pattern without affecting the mould
Shrinkage/contraction allowance To compensate for solid shrinkage
Shake/rapping allowance For clearance between mould and pattern
Distortion/Camber allowance V or U shape castings
Highest liquid+solidification shrinkage: Aluminum
Highest total shrinkage: Steel
Shrinkage allowance of steel 1mm per 20mm (approx)
5.2.2 Parts
Patterns
Solid or single piece pattern
Split piece pattern
Gated pattern: Gate and runner are included in pattern
Loose piece pattern: For parts with internal webs
Match plate pattern
Sweep pattern
Flowbard pattern
Sprue
Best shape is parabolic tapered
Straight tapered is used
Riser
Compensates shrinkage in liquid phase and shrinkage during solidification
Vents
For easy escape of air
Cores
To produce internal cavities
Chills
To avoid hot tear
Directional solidification
Uniform solidification
Paddings
Avoid Erosion
Maximize heat transfer
Directional solidification
www.eggsam.com 5.2. CASTING [46]
Chaplets
To support core
Directional solidification
Molding sand
Refractoriness: Ability to withstand high temperature
Green strength: strength of moulding sand with moisture
Collapsibility: Ability of the mould to not to resist the shrinkage of the metal. Adding
saw dust improves collapsibility
Dry strength: Strength of the mould after the mould cavity dries when the molten metal
is poured.
Permeability: Ability to allow gases to escape (up to 8% water, permeability increases
and then decreases.)
VH
Permeability number P N =
P AT
V=2000cc
H=5.08cm
P=10g/cm2
A=20.26cm2
3007.2
PN =
T
Time in minutes
Flowability
Strength
Hardness
Adhesive property
Cohesive property
Conductivity
Thermal expansion
Collapsibility
Additives
Saw dust, wood flour Improves green strength and collapsibility
Organic binders
Starch and dextrin
Improves skin hardness
Iron oxide and
Improves hot strength
aluminum oxide
Coal dust, sea coal, improves surface finish and resistance to metal
silica flour penetration
Solidification
By skin forming or dendritic growth
Casting volume
Casting yield =
Casting volume+ gating volume
www.eggsam.com 5.2. CASTING [47]
Gating system
Ferrous castings
V
Top gate Time for filling t = √m
Ag 2gh
Turbulence and splashing
Non-ferrous castings
Time for filling
Bottom gate Am 1 √ √
t=2 √ ( h − h − hm )
Ag 2g
No splashing or turbulence
Parting gate No splashing or turbulence
Multiple ingates
Step gate Fast
No splashing or turbulence
Gating ratio = AS :AR :AG = Sprue:Runner:Ingate
Non-pressurized: 1:4:4 or 1:2:2 (For highly reactive metals like Al, Mg,...)
Pressurized: 1:2:1
Cold shut Two streams of molten metal does not fuse properly
due to cooling
Hot cracking
Hot tear residual stresses in the material causes the casting to
fail as it cools
Mold shift Due to misalignment between two halves
Core shift Core moves from its position
5.2.5 Moulds
Sand moulding
Shell moulding
Expendable moulds Investment moulding
Full moulding
CO2 moulding
Centrifugal
Permanent moulds Die casting
Slush casting
Squeeze casting
www.eggsam.com 5.2. CASTING [49]
Mass production
Fast cooling
Gravity die casting Aluminunm Piston in automobile
Fine grain
High strength and hardness
Hot Chamber die High production rate Low melting point
casting High accuracy Lead, Zinc, Tin, Magnesium
High melting point, Non-ferrous
High melting point alloys of Al,
Cu, Brass, Magnesium
Cold chamber
Carburators
Pressure die casting
Crank case and crank shaft
valve bodies
Fuel injection pump parts
Cleaning of castings
Fettling
Shot/ Sand blasting
5.2.6 Cupola
Advantage: Continuous melting
Low cost of melting
Chemical composition can be controlled
Good temperature control
www.eggsam.com 5.3. METAL CUTTING [51]
Hot hardness: Minimum temperature above which the increase in temperature causes
sudden decreases in hardness
HHS < Borazon < Ceramic < Cermet < Carbide < Diamond
HSS
W, Mo: To increase hot hardness of material
Cr: Increase strength of resistance to deformation
V: Increase wear resistance
W based HSS has more wear resistance than Mo based
Carbide tools
P type: For ferrous
K type: For non-ferrous, non-metal, cast iron
Low number: Finish machining
High number: Rough machining
Coating: Aluminum and Zirconium
Cermets
Ceramics + Metals
Surface roughness
f2
Peak to valley height Hmax =
8R
Hmax
Center line average value Ra =
4
True feed = f cos Ψ
f
Hmax =
tan Ψ + cot Ψe
Forces
Thrust/radial force Fy = Ft cos λ
Axial/ feed force Fx = Ft sin λ
Merchant’s Analysis
1 t2 cos(φ − α)
Chip thickness ratio = = ≥
r t1 sin φ
1
r cos α
tan φ =
1 − r sin α
r : Chip reduction ratio
Shear strain =cot φ + tan(φ − α)
wt1
Area of shear plane =
sin φ
wt1
Shear force Fs = τ
sin φ
t1 : Uncut chip thickness
=feed × sin λ (Turning)
=depth of cut (orthogonal machining)
Fs = Fc cos φ − FT sin φ
Ns = FT cos φ + FC sin φ
N = FC cos α − FT sin α
F = FC sin α + FT cos α
www.eggsam.com 5.3. METAL CUTTING [56]
FS cos φ − sin φ FC
=
NS sin φ cos φ FT
N cos α − sin α FC
=
F sin α cos α FT
FT
+ tan α
F
tan β = C
FT
1− tan α
FC
F
µ = tan β =
N
1
ln
If β > 45o : Classical friction theorem µ = π r
−α
2
For orthogonal cutting FR = 0
Assumptions for Merchant’s circle
Cutting edge is sharp and straight
Rigid, perfectly plastic, homogeneous material
Orthogonal cutting
Shear zone is approximated by a straight line
No BUE
Cutting Power=Fc Vc
FC
Specific cutting power =
wt1
Ernest and Merchant Theory
Minimum power consumption during machining
2φ + β − α = 90O
2φ + β − α = cot−1 K (Modified Merchant’s theory, Mohr’s theory)
K: Machining constant
Stabler Theory
2φ + α − 2β = 90o
Velocity
V Vc Vs
= =
sin(90 − (φ − α)) sin φ sin(90 − α)
VS : Shear velocity
VC : Chip velocity
www.eggsam.com 5.3. METAL CUTTING [57]
sin(φ + β − α)
Contact length lf = t1
sin φ sin β
Cutting fluids
Increase heat dissipation and decrease power consumption
Act as lubricant
Properties
-High conductivity
-Should not fume, foam
-Should not react with workpiece or tool
-Low viscosity
Cast iron No cutting fluid
Low speed : Neat oils + EP Additives
Steel Medium speed : 1:10 water emulsion + EP additives
High speed : 1:100 water emulsion
Aluminum Neat oil + EP additives at low speed. No cutting fluid at high speed
Magnesium Only neat oils
brass/Bronze Only neat oils
Tool failure
Diffusion wear
Adhesion wear (Spot welds)
Abrasion wear (Faylite pockets)
Fatigue wear
Plastic deformation: Due to the high temperature (Temperature > Hot hardness tem-
perature)
Mechanical breakage: Due to impact loads
Oxidation wear
Flank wear
Crater wear
Misc √
Maximum temperature on Rake face ∝ f eed
If cutting speed is increased, cutting force remains same. (But due to heating, it can
reduce slightly)
Broaching Super finishing operation, multipoint cutting operation, For making key ways
and internal gears
5.4 Machining
5.4.1 Tool life
Effect on tool life V > f < d
Machinability
Tool life
www.eggsam.com 5.4. MACHINING [59]
Surface finish
Cutting forces
MRR
Specific cutting energy
Shear angle (Higher shear angle-¿ Better machinability
Vt
Machinability index = 100
VS
VS : Cutting speed of standard free cutting steel for 60 minute tool life
Vt : Cutting speed of metal for 60 minute tool life
5.4.2 Grinding
• Surface grinding
• Cylindrical grinding
• Center-less grinding
• Form grinding
• Abrasive belt grinding
• Manual grinding
• Creep feed grinding
S: Structure 0 Dense
16 open
V Vitrified
S Silicate
B: Bond type B Resinoid
R Rubber
M Metal bond
Open structure: Ductile workpiece
Closed structure: Brittle workpiece
Wheel Truing: A redressing process by which the wheel is restored to its true shape
Wheel dressing: The process of making new sharp edges on grains. Required due to
grazing.
Grinding wheel wear
Grain wear, grain fracture and bond fracture
Finishing operations
Honing: To make fine surface finish for holes
Lapping: Finishing operation for flat surfaces
Polishing
Buffing
Deburring
5.4.3 Drilling
Drills are made by forging
HSS is the tool material
d
Width of chip W =
2 sin β
2β: Point angle
Rake angle: Angle formed between a plane containing the drill axis and the leading edge
of the land. Positive for right hand flute, negative for left hand flute, zero for parallel
flute.
Point angle or cutting angle. Small point angles are used for cutting ductile materials
and large point angles are used for cutting brittle materials
To make a hole
Drilling Drill is the cutting tool used
Oblique cutting process
HSS
Enlarging an existing hole+Better finish
Boring Accuracy =0.125mm
Cannot increase the length of the hole
Finishing process
Reaming Surface finish ±0.005mm
Negligible change in diameter
Reamer has multiple cutting edges
5.4.4 Milling
Peripheral milling
Slab milling
End milling
Gang milling
Straddle milling
Upmilling
No backlash
s
d d
sin θC = 2 1−
D D
s
d d
Mean chip thickness =ft 1−
D D
www.eggsam.com 5.4. MACHINING [62]
d: depth of cut
D: Diameter of cutting tool
Face milling
1 p
Compulsory approach = (D − D2 − wi2 )
2
wi = width of work + 2Offset
Hear hobbing
Fast process
Cylindrical tool with slots and gashes
Looks similar to a worm gear
Rotates continuously. A continuous process
Cannot cut internal gears
Helical, worm and spur gears can be made
Gear milling
5.4.6 Planning
More then one single point cutting tool
Work piece is reciprocating
www.eggsam.com 5.5. METAL FORMING [63]
5.4.7 Shaping
One single point cutting tool
Tool reciprocates
5.5.1 Rolling
∆h
cos α = 1 −
D
∆h: Reduction in thickness
For unaided entry µ ≥ tan α
∆hmax = µ2 R
www.eggsam.com 5.5. METAL FORMING [64]
σ0 : Flow strength
σyt ≤ σ0 ≤ σut
Maximum bite angle = tan−1 µ
A0
= 2.71
Amin
Rolling defects
Wavy edges
Spread
Crocodile crack
Alligatoring
5.5.2 Forging
Fullering or swaggering, Flattering, Finish, Cut off
Flash: extra material deposited in gutter
Drop forging
Open die forging
Cogging
Press forging (Used for making coins)
Roll forging
Precision forging
Impression forging
Forging defects
Cracks, Fold, Barreling
5.5.3 Extrusion
2
do
Extrusion ratio =
df
Ao
True strain =ln
Af
ho
Extrusion strain =
hf
Ao
Force required =KAo ln
Af
Forward or direct extrusion
Backward extrusion or indirect extrusion
Hydrostatic extrusion
Impact extrusion: Tooth paste tube
Johnson’s
equation
Ainitial
σd = σ0 a + b ln
Af inal
Defects
Pipe defect, tail pipe, Fishtailing
Surface cracking
Internal cracking
www.eggsam.com 5.5. METAL FORMING [65]
When µ = 0,
A0
σ = σy ln
A1
Tube drawing
B
1+B h1
σ = σy 1−
B h2
Defects
Center line cracks
seams
residual stresses
Punching
To punch holes
A shearing operation
Shear is given on punch
Punched out material is waste
Punch=Size of hole
Die= punch size + 2 radial clearance
Blanking
Blanked out material is the product
Shear is given on die
Punch= die size - 2 radial clearance
Die=Size of product
www.eggsam.com 5.6. WELDING [66]
Deep drawing
Press working
For making steel tumblers
Defects of drawing
Flange Wrinkle/earing: due to insufficient blank holder pressure.
Earing: due to anisotropy induced by rolling operation or due to non-uniform clearance
between tools.
Wall wrinkles
Fracture
Miss strike
Orange peal
5.6 Welding
Solid state welding Liquid state/fusion Solid/liquid
Autogeneous Homogeneous Heterogeneous
Explosive Resistance Brazing
Ultrasonic Chemical reaction Soldering
Friction -gas
Forge -thermite
Diffusion Arc
SMAW, TIG, MIG, PAW
Hand peening is a stress relieving process and it consists of hammering the weld along
the length with the peen of the hammer while joint is hot.
Weld pool
Weld bead
Reinforcement
Penetration
Root gap
Toe
Root
Throat
Deposition rate
DC Straight polarity
- Work piece is positive
- more depth of penetration
- Weld deposition rate is less
DC Reverse polarity
- Workpiece is negative
- Less depth of penetration
- Weld deposition rate is high
- Used for thin sheets
- Constant current type
- Constant voltage type
It Vt
- + =1
IShort circuit VOpen
OCV
- V=OCV- I
SCC
Arc on time
- Duty cycle =
Arc on time + idle time
Arc blow
- Deflection of electric arc due to the magnetic field formed in the material during welding
- Weld Splatter
- Provide flux coating to reduce arc blow
- Less chance in DC
Flux coating
Electrode designation
TFEPSX
T: Type of electrode manufacturing
F: Type of flux coating
E: Position of electrode
P: Polarity
S: Strength of electrode
X: Specific information about electrode
SMAW
Shielded Metal Arc Welding
Most commonly used
Electrode coating provides the shielding gas.
www.eggsam.com 5.6. WELDING [68]
TIG, GTAW
Tungsten Inert Gas
Gas Tungsten Arc Welding
Non-consumable tungsten electrode
Carried out in inert atmosphere. (He, Ar, Ne, CO2 , N2 )
Thorium and Beryllium are added to increase the thermal resistance of Tungsten
For Al and Mg alloys, AC welding is used
For all other materials, Direct current straight polarity is used
Used for thin welding
Applications: Aerospace and automobile industries
MIG, GMAW
Metal inert gas
Gas Metal Arc Welding
Consumable electrode (wire) is converted to molten drops
DC reverse polarity or AC are used for welding Al, Mg, Cu,...
Used for thick welding
Can be easily automated
For welding Stainless steels, Al, Mg, Cu, Ni alloys
Application aircraft and automobile industries
For complete combustion of 1 unit volume of acetylene, 2.5 unit volume of acetylene
is required
1 part of it is provided from the cylinder and 1.5 is obtained from atmosphere.
Gas cutting
Cutting using oxidation of Iron
Al cannot be cut using gas cutting
Spot welding
Spot welding is adopted to weld two overlapped metal pieces between two electrode
points.
Indentation is created by the force from the electrodes
The lapped pieces of metal√ are heated in a restricted area.
Diameter of nugget d = 6 t
Mainly used for lap welding thin sheets
Automobile and refrigerator bodies
www.eggsam.com 5.6. WELDING [70]
Seam welding
Electrodes are in the form of wheels
Process is similar to spot welding
Wheels roll creating series of spot welds
Creates leak proof joint
Projection welding
There are projections on one plate made by embossing
Copper plates are used instead of electrodes
Can be used to weld nuts and bolts to plates
Used to join a network of wires
Flash welding
Flash butt welding
Arcs form and soften the metal as the members come closer.
Force is applied to weld the parts by plastic deformation
Very high current 10,000A
Mild steel, Medium carbon steel, Alloy steels, Al, Ti
5.6.15 Soldering
Solder is an alloy of Lead and tin
Used in electronics industry
Solder melts at relatively low temperature
5.6.16 Brazing
Material is a alloy of Copper, Zinc and silver. It is called Spelter.
Stronger than soldering
Used to connect pipes and make leak proof joints
5.6.17 Defects
Porosity
Due to entrapment of gas bubbles
Can be reduced by proper selection of filler material, preheating the weld area, cleaning
the weld area and reducing the welding speed
Slag inclusion
Caused by materials getting trapped in the weld. (Electrode coating materials, oxides,...)
To prevent slag inclusion, clean the weld surface before next layer is deposited, provide
shield gas
www.eggsam.com 5.7. NON-TRADITIONAL MACHINING [72]
Undercut
Incomplete fusion
Overlapping
Weld spatter
Weld cracks
Weld decay
Electrolyte
-Large electric conductivity
-Good chemical stability
-Inexpensive
-Should not cause corrosion
-Low viscosity
-Non toxic
-Chloride solution in water (salt solution)
Applications
Turbine blades
Large through cavities
Blind complex cavities
www.eggsam.com 5.7. NON-TRADITIONAL MACHINING [73]
IM & OR
Output
Productivity =
Input
75
www.eggsam.com 6.1. CPM / PERT [76]
Free float
Amount of time the activity can be delayed without affecting the succeeding activity
Independent float
Time by which an activity can be adjusted without affecting the preceding or succeeding
activity
Seasonal inventory
Anticipatory inventory
Decoupling inventory
Transit or pipeline inventory
Direct inventory: Inventories that are directly a part f production and become a part
of the final product.
Ex: Raw materials, In process inventories, Purchased parts, Finished goods
r
2DCo p
r
Production model Q=
Cc − d)
(p r
√ (p − d)
TVC= 2DCo Cc
p
r r
2DCo (Cc + Cs )
With shortage Q=
Cc rCs
√ Cs
TVC= 2DCo Cc
r r Cs + Cc
2DCo Cc
OS=
Cc (Cc + Cs )
Model sensitivity
Q∗
T V C(Q) 1 Q
= +
T V C(Q∗ ) 2 Q∗ Q
Demand-profit / Static inventory
perishable items
-70-90-100
-Vital few to Trivial many
-Consumption analysis
HML
High Medium Low
-based on the unit price of the product
VED
Vital Essential Desirable
-Based on importance of the product
SDE
Scarce Difficult Easy
-Based on availability
XYZ
Based on inventory value
X: Items with high inventory value
Z: Items with low inventory value
FNSD
Fast Normal Slow Dead moving items
Based on the speed of usage of items
6.3 Forecasting
6.3.1 Qualitative methods
Educated Guess
Based on a person’s judgment based on experience and intuition
Delphi method
Panel of experts
Questionnaires are used.
Long range
For new product, technology, changes in society,...
Market research
Market trial
Exponential smoothing
Ft = Ft−1 + α(Dt−1 − Ft−1 )
Ft = Ft−1 + α(error)
Linear Regression
Σy = na + bΣx
Σxy = aΣx + bΣx2
y = a + bx
(Yc − Ȳ )2
P
2
Coefficient of determination r = P
(Y − Ȳ )2
Coefficient of correlation =r
6.3.3 Error
n
(Dt − Ft )2
P
1
Mean Square Error MSE=
n n
P
|Dt − Ft |
1
Mean Absolute Deviation MAD =
n
n
P
(Dt − Ft )
1
Bias =
n n
P
(Dt − Ft )
Cumulated deviation 1
Tracking signal = =
M AD M AD
Bias × n
Tracking signal =
M AD √
Upper limit for tracking signal =3 M SE
www.eggsam.com 6.4. BREAK EVEN ANALYSIS [82]
Z Entering variable
Basic ↓
variables
Pivot element →Leavingvariable
Feasible solution: Any values of the basic variable that obey the constraints
Unique solution: The number of zeros = number of basic variables
Multiple solution: Number of zeros > number of basic variables
Unbounded solution: All numbers in replacement ratio column is negative or infinite
No solution: Artificial variable remains in the final solution
Degenerate solution: One or more basic variable becomes zero.
6.7 Transportation
Degeneracy: Occupied cells < m+n-1
6.8 Assignment
Hungarian method
www.eggsam.com 6.9. WORK STUDY [85]
Material Science
87
www.eggsam.com 7.1. TESTS [88]
7.1 Tests
Test name Details
Tensile testing On universal testing machine
Compression
On universal testing machine
test
Measures fracture toughness
Izod Test Sample is fixed at one end and the other end is free
Non-uniform stress
Measures fracture toughness with better accuracy
Charpy Test
Sample is fixed at both the ends
Herbert
A shower of metal balls
cloudburst
To find defects
Hardness test
Spiral test Fluidity
Cupping test Formability
Dye penetrant
To find surface defects
method
To find hardenability
Jominy end
Austenite to martensite
quench test
50% - 50% Pearlite and martensite, Jominy distance
√
Gauge length =5.65 A0
www.eggsam.com 7.1. TESTS [89]
Hardness Tests
Test name Details
A qualitative test
Mho Test/ Scratch
Used as a preliminary test
Test
Used for Low-medium-high hard materials
Used for Medium hard materials
Indenter: Spherical, d=10mm
P=50-120kg
Brinell Hardness Test 2P
BHN = √
πD(D − D2 − d2 )
UTS=3.6 × BHN, for normalized plain carbon steels
UTS=3.2 × BHN, for tempered plain carbon steels
Applicable to all types
Indenter: 120o diamond cone (brail)
Rock Well Test P=1-150kg
1
Hardness ∝
t
For Medium hard materials
Indenter: 136o diamond square pyramid
Vickers Test P=50-120kg
1.854P
V HN =
davg
davg : Average diagonal indentation
For Low hard materials (Si, Ge, Ga, As,...)
Knoop Test / Micro 177o Bi-pyramidal base indenter
Hardness test P=1-1000g
KHN = 14.22P/L2
For soft plastic, thermo plastic, rubber, thin sheets ...
Shore Method
Diamond tipped indenter (hammer) in a glass tube
Shore’s Scleroscope
Height of reboundness ∝ hardness
For hard plastics, thermosets, composites,...
Barcol Method
Spring reading
P : Load at which indentation is produced
D: Diameter of indenter
d: Diameter of indentation
www.eggsam.com 7.2. PLASTICS [90]
Non-destructive testing
Visual inspection
Hang the casting in the air
Hammer test Gently strike with a
hammer and listen to the
sound
Radiography x-rays and γ rays
Defects that are open to
Liquid penetrant test the surface
Die penetrant test
Ultra-Sonic inspection
Hardness
Diamond > Silicon > Quartz > Topaz > Feldspar > Apatite > Fluorite > Calcite >
Gypsum > talcum
Ductility
Al > Cu > Zn > Mild steel
Elastic failure
Necking → Formation of small cavities → Cavities combine together → crack propaga-
tion → Fracture (Cup-cone)
7.2 Plastics
Thermoplastics Thermosetting plastics
Become soft on heating Becomes hard on heating
Recyclable Non-recyclable
Cyclic structure
Linear structure
addition polymerization
Polyvinyl Chloride (PVC) Epoxy
Polypropylene(PP) Polyester
Polyethylene(PE) Phenol formaldehyde(Bakelite)
Polystyrene(PS)
Poly Tetra Fluoro Ethylene
(PTFE, TEFLON)
Acrylic
Molecular mass
Degree of polymerization =
mer mass
www.eggsam.com 7.2. PLASTICS [91]
Condensation polymerization
Produces water or ammonia as by-product
www.eggsam.com 7.3. CERAMICS [92]
ABS Terpolymer
Polyethyene Addition polymerization
Polycarbonate Addition polymerization
Polystyrene Addition polymerization
Polyamide Natural (proteins) and artificial (Nylon)
Poly propylene Excellent fatigue strength
PTFE Low coefficient of friction
PVC Synthetic polymer, pipes, bottles,...
PMMA Poly Methyl Methacrylate
Polyether ether ketone,
PEEK bearings, piston parts, pumps, High-performance liquid
chromatography columns, compressor plate valves, and electrical
cable insulation.
Nylon Fabrics
Polyurethene Low-density flexible foam
Cyano-acrylate Adhesives
Neoprene Oil seal
Bakelite Electric switches
Araldite Adhesive
SBR Styrene buta diene rubber , Tyres
Kevlar Bullet proof vests
7.3 Ceramics
Extremely brittle
High thermal stability
High chemical stability
Corrosion resistance
High hardness
Silica (SiO2 )
Alumina
www.eggsam.com 7.4. CRYSTAL STRUCTURE AND DEFECTS [93]
Metals
-Has free electrons
Burger’s Vector
|~b| = 0: No defect
|~b| = 1: Point defect
|~b| =>: Line defect
www.eggsam.com 7.4. CRYSTAL STRUCTURE AND DEFECTS [95]
Line defects
Edge dislocation Screw dislocation
Glide and climb Glide
Dislocation lines are perpendicular to
Lies parallel to burger’s vector
Burger’s vector
Direction of movement of edge dislocation is Direction of movement is perpendicular to
in the direction of Burger’s vector Burger’s vector
Movement of edge dislocation is fast Movement of dislocation is slow
Explains plastic deformation and crystal
Explains plastic deformation
growth
Less shear force is required to make this High shear force is required to make this
defect defect
Tensile, compressive and shear stress fields
Only shear stress field is present
can be present
Termination of atomic plane in the middle
Will change surface properties significantly
of a crystal
Whisker
Movement of atomic planes is translation +
rotation
Grain boundaries restricts the motion to dislocation. It makes the material stronger.
Permanent deformation in metal or alloy is caused by movement of dislocations.
Amorphous solids
No regular arrangement of atoms
No sharp melting point
7.5 Alloys
Hume Rothary rules/Conditions
Difference in atomic radius of both atoms must be less than 15%
Valency of both the atoms should be same
Electro-negativity and electron affinity of both the atoms should be compatible
Atoms at grain boundary has more energy. So oxygen attacks there and causes cor-
rosion.
Chromium reacts with oxygen to produce Cr2 O3 , this gets in grain boundaries and pre-
vent corrosion by blocking oxygen
Weld decay: Corrosion at welded parts due to lack of Cr2 O3 , due to formation of
Chromium carbide during welding
Corsing or Miscibility gap: Due to sudden cooling, no time for diffusion, so concen-
tration gradient, causes cracks on hot working
Ni is added to stabilize austenite phase
Cr is added to stabilize ferrite phase
Carbon equivalent = %C + 31 %(Si+P)
www.eggsam.com 7.5. ALLOYS [98]
Alloying of steel
Steel Alloys
Cr: Reduces scaling
18% W or Mo
V: Abrasion resistance
HSS 4% Cr
Tungston: Hot hardness
1% V
Mo: Hardenability
0.67% C
Co: Hot hardness and wear resistance
High wear resistance
Hadfield Manganese steel ≈13% Ma High toughness
Bulldozer blades
Bronze 88% Cu
Utensils, bearings, bushes, wires,...
12% Tin
Brass Cu
Musical instruments
Zn
Chromel 90% Ni
Thermocouple
10% Cr
Constantan 55% Cu
Thermocouple
45% Ni
Duralumin 94% Al
Cooking utensils, tubes, rivets, sheets,...
4% Cu
Curie Point
-no change in crystal structure.
-Magnetic properties are changing
Ferro magnetism
Sometimes paramagnetic and sometimes diamagnetic
www.eggsam.com 7.6. PHASE DIAGRAMS [104]
Name Properties
FCC
Solid solution of Ferrite+Iron carbide
non-magnetic
Austenite (γ) in gamma iron
soft
Mn, Ni, Si are austenitic stabilizers
Not stable below 725o C
BCC
Highly magnetic
Ferrite δ, α 0.02% Carbon
soft
Ductile
BCC
δ 0.1% Carbon
1410o - 1540o
Orthorhombic
Cementite Extremely hard and brittle Fe3 C , 6.67% Carbon,
Magnetic below 200o C
α +Fe3 C
Pearlite Ferrite(87%)+Cementite Phase mixture
(13%)
Ledeburite Austenite+Cementite 4.3% carbon Mix
cooled slower than the rate required to
Hard form martensite but faster than the
Bainite
Brittle rate that would be required to form
pearlite.
Ferrite + Cementite
Martensite Hardest and brittle
Rapid cooling of HCS
Lower hardness and Formed by heat treatment of
Troosite
brittleness than martensite martensite
Lower hardness and
Sorbite Formed by heating martensite
brittleness than troosite
www.eggsam.com 7.6. PHASE DIAGRAMS [105]
Ferrimagnetic Ferromagnetic+
Anti-ferromagnetic
Reduce Hardness
Increase ductility
Annealing Heat austenite temperature and cool slowly
Improve machinability
in furnace
Relieve internal stresses
Refine grain size
Heat steel to 50o C to 70o C above the upper
critical temperature, Hold there, Slowly
Reduce
cool in furnace.
Full Annealing Hardness/brittleness
In hypoeutectoid steels, Austenite becomes
Increase
coarse pearlite and ferrite structures.
ductility/toughness
In hyper eutectoid steels, Austenite →
Pearlite+Cementite
To remove effects of cold
Process Annealing Heating below lower critical temperature,
working (relieve stress)
usually used in low carbon steels
Make soft
To increase machinability
Spheroidise annealing Heat near lower critical temperature, slow
in MC or HC steels
cooling in furnace
Increases ductility
To make uniform
Diffusion annealing 1150o C and slow cooling
composition
homogenizing
Usually done after welding
www.eggsam.com 7.7. HEAT TREATMENT [108]
Cynaniding NaCN
<3000o C Carburizing
Flame hardening Guideways of lathe 3150o C Neutral
3480o C Oxydizing
Induction hardening To harden surface
Hardness order
Nitriding > Cyaniding > Carburizing
Brine < Water +NaOH < Water < Oil < Air
Cooling rate
Air < Oil < Water < Brine
-for Al alloys
Overaging (coarsening of precipitate particles)
Artificial aging
The iron-carbon diagram is determined under equilibrium and TTT curve is determined
under nonequilibrium condition.
The main purpose of spheroidising treatment is to improve machinability of high carbon
steel.
The martensitic transformation is a process of shear, that occurs without any need for
diffusion so there will be no change in composition in this process.
7.8 Nanomaterials
At least one dimension less than 100nm
Zero size: Particles
One dimensional: only one large dimension
Two dimension: sheets
Three dimension: Cubes
Manufacturing
Top-down approach
Bulk material is converted to nano
Mechanical grinding,atomization
Bottom-top approch
Atoms combine to produce nano
Sol-gel technique
Physical/Chemical vapor deposition
7.9 Misc
Cottrell atmosphere: Due to diffusion, the interstitial Carbon gets accumulated in
dislocation sites in iron. (More energy is needed to break it and it causes upper yield
point)
Bauschinger effect:Unload the materials from the region of work hardening- reverse-
load again, Tensile yield strength increases, compressive decreases
Strain hardening
σf = Kn
n: work hardening exponent (n=0.3 for steel, n= 0.05 for Al)
K: strength coefficient
at UTS = n
Usually 0 < n < 1
www.eggsam.com 7.9. MISC [110]
Surface hardening
Shot blasting: For heavy material, steel balls
Shot peening: For small size material, manual hammering
Sand blasting: For thin material,
Hall-Patch Equation
K
σy = σ0 + √
d
σy : Yield strength
σ0 : Base strength of material
k: Constant
d: Grain diameter
Corrosion
Tin plated iron sheet: tin is anodic to iron
Galvanized: zinc is cathodic to iron
Chapter 8
8.1 Microprocessors
Microprocessor is a controlling unit of a micro-computer, fabricated on a small chip ca-
pable of performing ALU (Arithmetic Logical Unit) operations and communicating with
the other devices connected to it.
Registers
Store the data temporarily during the execution/runtime of the program
Memory
stores the information (data& instructions) in binary form.
Read Only Memory (ROM)
Stores items that the computer needs to execute when it is first turned on
Random Access Memory (RAM)
Stores user programs and datas temporarily. RAM is a volatile memory.
System Bus
1) Control bus
2) Data bus
3) Address bus
111
www.eggsam.com 8.2. MICROCONTROLLER [112]
8.1.3 Types
RISC: Reduced Instruction Set Computer
CISC: Complex Instruction Set Computer
8.2 Microcontroller
Microcontroller is integration of all microprocessor and input and memory other periph-
erals in a single chip.
Arduino is a microcontroller
Atmega: Microcontroller
8.2.1 Components
Central processing unit(CPU)
Random Access Memory)(RAM)
Read Only Memory(ROM)
Input/output ports
Timers and Counters
Interrupt Controls
Analog to digital converters
Digital analog converters
Serial interfacing ports
Oscillatory circuits
www.eggsam.com 8.2. MICROCONTROLLER [113]
Absolute encoder
Position known more accurately
360o
Resolution = N
2
Inductosyn
Position sensor
Linear and angular displacement
Based on electromagnetic induction
Piezoelectric accelerometer
Accelerometer
State space Representation
Controllability
www.eggsam.com 8.7. ACTUATORS [116]
8.7 Actuators
Hydraulic actuators
Liquid pressure energy to mechanical power
Power = Pin Q
8.8 Robotics
The notations used follows Prof. Ashitava Ghosal ’s NPTEL course.
(https://nptel.ac.in/courses/112/108/112108093/)
Students are strongly advised to visit the course at least once OR read the book
https://www.amazon.in/Robotics-Fundamental-Concepts-Ashitava-Ghosal/dp/0195673913
Cylindrical robot(R2P)
- Axes form a cylindrical coordinate system.
- Uses: assembly operations, handling at machine tools, spot welding, and handling at
die casting machines.
SCARA robot(2R1P)
-Selective Compliant Assembly Robot Arm or Selective Compliant Articulated Robot
Arm.
- It’s a robot which has two parallel rotary joints to provide compliance in a plane
- Uses: pick and place work, application of sealant, assembly operations and handling
machine tools
- This is an example of a basic SCARA robot configuration.
- Note the three rotary joints and the single vertical joint used in this horizontally artic-
ulated configuration. (3R1P)
- Can move very fast.
- Best suited to planner task
Parallel robot
Used as a mobile platform handling cockpit flight simulators
PUMA
- Industrial robot.
- Programmable Universal Machine for Assembly, or Programmable Universal Manipu-
lation Arm
- Functions like a human arm.
- A total of 6 variables are required, for specifying the position and orientation of a rigid
body in space.
- PUMA has 6 axis of rotation
Transformation
matrix
1 0 0
Rx = 0 cos θ − sin θ
0 sin θ cos θ
cos θ sin θ 0
Ry =
0 1 0
− sin θ cos θ 0
cos θ − sin θ 0
Rz = sin θ cos θ 0
0 0 1
www.eggsam.com 8.9. TRANSDUCER [118]
Cθi −Sθi 0 ai−1
i−1[T ] = Sθi Cαi−1 Cθi Cαi−1 −Sαi−1 −Sαi−1 di
i
Sθi Sαi−1 Cθi Sαi−1 Cαi−1 Cαi−1 di
0 0 0 1
Transformation based on current axis
0 0 1 2
3 [T ] =1 [T ]2 [T ]3 [T ]
Transformation based on a fixed axis
3 2 1 0
0 [T ] =3 [T ]2 [T ]1 [T ]
8.9 Transducer
The device which converts the one form of energy into another is known as the trans-
ducer.
Active transducer does not use any external power source for producing the output.
Passive transducer requires the additional energy source for working.
Piezoelectric transducer
d
V = tP
r 0
d
g=
r 0
t: Thickness of crystal
d: charge density of crystal
P: pressure on crystal
g: Voltage sensitivity of crystal
Photoelectric Transducer
Converts the light energy into electrical energy.
Made of semiconductor material. Photoemissive Cell
Photoconductive Cell
Photo-voltaic cell
Photodiode
photo-diode
Phototransistor
Chapter 9
Fluid Mechanics
9.1 Hydrostatics
1 Torr= 1mm Hg
1 bar=100kPa
1 poise=0.1Ns/ m 2
1 Stoke=10−4 m 2 / s
1kgf=9.81N
Mass
Density =
Volume
Weight
Specific weight =
Volume
1
Specific volume =
Density
Density
Specific gravity or relative density =
Density of standard fluid
Mean free path
Knudsen number Kn =
Characteristic length of flow
Kn < 0.01 : Continuous fluid
Kn > 0.01 : Continuum does not hold
119
www.eggsam.com 9.1. HYDROSTATICS [120]
4σ
Excess pressure inside a drop =
d
8σ
Excess pressure inside a bubble =
d
Capillary raise/fall
Due to Cohesion and Adhesion
4σ cos θ
h=
ρgd
Buoyancy and floatation
Buoyant force = Volume immersed × density of fluid × g
I
Meta-centric height GM = − BG
V
I
Meta-centric radius BM =
V
V: Immersed volume s
K2
Period of rolling T = 2π
gGM
Increase in GM =⇒ Decrease in time period oscillation =⇒ Increase in stability =⇒
Decrease in comfort
Decrease in GM =⇒ Increase in time period oscillation =⇒ Decrease in stability =⇒
Increase in comfort
K: Radius of gyration about axis of rolling
I: Least area moment of inertia of the body at water surface
h: Position of Center of gravity from the surface
Name A µ n Examples
S Ideal Solid - - - -
T Thixotropic τyield µ <1 Printer ink, lipstick
B Bhingam Plastic τyield µ 1 Sewage sludge, tooth paste, drilling mud
R Rheopectic fluid τyield µ >1 Gypsum
P Pseudo Plastic 0 µ <1 Blood, milk, suspension paints, paper pulp
N Newtonian 0 µ 1 Air, water
D Dilatant fluid 0 µ >1 Butter, Rice starch, Sugar sol
I Ideal fluid 0 0 - -
Thixotropic and Rheopectic fluids have time dependent viscosity.
Bhingam plastic isalso called Ideal plastic
∂u ∂v
τxy = µ +
∂y ∂x
B/T
For Liquids µ = Ae √ (Andrade’s equation)
a T
For Gases µ = (Sutherland equation)
1 + b/T
A, B, a and b are constants.
T is absolute temperature
Rheology: Study of non-Newtonian fluids
Sound r r
K ∂P
Velocity of sound in fluid = =
ρ ∂ρ
IGG sin2 θ
Center of pressure yp = h +
Ah
h : Center of gravity from surface
A : Area of the surface
IGG : Area moment of inertia about axis through CG, parallel to surface
θ : Angle of the surface with horizontal
Absolute pressure = Gauge pressure + Local atmospheric pressure
Net positive suction head P − Pv
Thoma’s cavitation number σc = 2
=
ρV /2 ρV 2 /2
Patm,abs Pvappressure
NPSH= − hsuction − hLoss −
ρg ρg
Uniform flow: Velocity at every point in the flow at a given time is same.
Non-uniform flow: Velocity may vary from one point to another in the flow.
Streak line: Locus of points of all the fluid particles that have passed continuously
through a particular spatial point in the past
www.eggsam.com 9.2. KINEMATICS OF FLUID FLOW [123]
Stream line: Curve that is instantaneously tangent to the velocity vector of the flow
Path line: Trajectory of a fluid particle
4A
Hydraulic diameter =
P
Cartesian coordinates Cylindrical coordinates
∂u ∂u ∂u ∂u
ax = +u +v +w ∂vr ∂vr vθ ∂vr ∂vr vθ2
∂t ∂x ∂y ∂z ar = + vr + + vz −
∂v ∂v ∂v ∂v ∂t ∂r r ∂θ ∂z r
ay = +u +v +w ∂vθ ∂vθ vθ ∂vθ ∂vθ vr vθ
Acceleration ∂t ∂x ∂y ∂z aθ = + vr + + vz +
∂w ∂w ∂w ∂w ∂t ∂r r ∂θ ∂z r
az = +u +v +w ∂vz ∂vz vθ ∂vz ∂vz
p∂t ∂x ∂y ∂z az = + vr + + vz
∂t ∂r r ∂θ ∂z
a = a2x + a2y + a2z
∂Ψ
=v 1 ∂Ψ
∂x ur =
∂Ψ r ∂θ
Stream = −u
∂y
function
Line of constant stream ∂Ψ
dx dy uθ = −
function = ∂r
u v
∂φ
= −u
∂x ∂φ
∂φ ur =
= −v ∂r
Potential ∂y
function ∂φ
= −w 1 ∂φ
∂z uθ =
dy u r ∂θ
Equipotential line: =−
dx v
∂u ∂v ∂vr ∂vθ
Cauchy = r =
∂x ∂y ∂r ∂θ
Reimann ∂v ∂u ∂vr ∂vθ
equations =− = −r
∂x ∂y ∂θ ∂r
Rotation of fluid particles
1
ω = ∇×V
2
1 ∂w ∂v
ωx = −
2 ∂y ∂z
1 ∂u ∂w
ωy = −
2 ∂z ∂x
www.eggsam.com 9.2. KINEMATICS OF FLUID FLOW [124]
1 ∂v ∂u
ωz = −
2 ∂x ∂y
Vorticity = 2ω
Circulation = Vorticity × area
~ ~
H
Circulation Γ = V · ds
Irrotational flow
ω=0
∇2 φ = 0
∇2 Ψ = 0
φ exists
Possible,steady,incompressible, irrotational flow if ∇2 φ = 0
Possible case of flow if Ψ exists
Vortex flow
v2
dp = ρ dr − ρgdz
r
Forced vortex Free vortex
v =ω×r vr = const
ω 2 r12 P1 ω 2 r22 P2 v 2 P1 v 2 P2
z1 − + = z2 − + z1 + 1 + = z2 + 2 +
2g ρg 2g ρg 2g ρg 2g ρg
Doesn’t follow Bernoulli equation Obeys Bernoulli equation
External force is required No external force
Bucket on rotating table Bottle, hole at bottom
Washing machine Kitchen sink
Whirlpool in a river
Tornado
ω2 2
Equation of free surface of fluid in a rotating cylinder Z = h0 − (R − 2r2 )
4g
h0 : Height of fluid at rest
ω 2 R2
Maximum height difference ∆Z =
2g
Where h0 is the initial height of fluid
Z is the height of fluid from the bottom of the cylinder
Water hammer
2L
Critical closure time Tc =
C
ρV L
Gradual closing of valve, Pressure rise =
T
Sudden closure of a rigid pipe, Pressure rise=ρvC
C: Speed of sound
Pitot tube
Stagnation pressure = static pressure + dynamic pressure
Static pressure: P
ρV 2
Stagnation pressure: P +
2
Pitot tube: Stagnation pressure - static pressure
P
Piezometric head= + Z
ρg
Piezometer tube
Direct pressure, no role for velocity
Triangular weir
8 √ θ
Q = Cd 2g tan H 5/2
15 2
Rectangular notch
2 √
Q = Cd L 2gH 3/2
3
umax
V =
2
32µV L f lV 2
Head loss hf = = [Hagen Poiseuille Formula]
ρgD2 2gD
∂τ ∂p
=
∂y ∂x
∂p r du
τ =− = −µ
∂x 2 dr
dp R
τmax = −
dx 2
P2 P1
dp = + z2 − + z1 Change in pressure due to change in piezometric head
ρg ρgp p
Friction velocity V ∗ = τ0 /ρ = V f /8
64
f=
Re
KE correction factor
KE based on actual velocity 1 R 3
α= = u dA
KE based on average velocity AU 3 A
For uniform flow, α = 1
For laminar flow inside a circular pipe α = 2
For turbulent flow inside a circular pipe α = 4/3
www.eggsam.com 9.6. TURBULENT FLOW [129]
9.7 Losses
Name Equation
f lV 2
Friction(Major loss)1 hf =
2gd
(v1 − v2 )2
Sudden expansion(Minor loss) hl =
2g
(vc − v2 )2
Sudden contraction(Minor loss) hl =
2g
v2
Entrance loss(Minor loss) hl = 0.5
2g
v2
Exit loss(Minor loss) hf =
2g
V2
Bends and other fittings hf = K
2g
1: Darcy-Weisbach formula
Coefficient of friction f 0 = 4f
Equivalent pipe
The equivalent pipe should have same head loss and flow rate of the earlier combination.
Le Li
When pipes are in series 5 = Σ 5
De Di
D
When similar pipes are parallel, d = 2/5
n
Power transmission through pipes
Power transmitted = ρg(H − hf )Q
H − hf
Efficiency =
H
H
Maximum efficiency at hf =
3
Maximum efficiency= 66.67%
Hydraulic gradient line
Shows piezometric head
www.eggsam.com 9.8. BOUNDARY LAYER FLOW [131]
Coefficient of drag
τ0
local: Cx = 1 2
2
ρu
FD
Average: CD ∗ = 1
2
ρAU 2
RL ρU 2
Drag force FD = 0 τx bdx = Cf x Bx , Bx = Ax
2
2 2
Drag force ∝ ρL V
1
Friction drag force FD = CD ρU 2 A
2
τ0
Local friction coefficient, Cf x =
ρv 2 /2
Blasius solution
5x
δ=√
Rex
0.664
Cf x = √ (Local skin friction coefficient)
Rex
1.328
CD = √ (Average drag coefficient)
ReL
τ0 : Shear stress on surface
Cubic
u 3 y 1 y 3
= −
U 2δ 2 δ
K=4.64
Linear
u y
=
U δ
K=3.46
Chapter 10
Fluid Machinery
V: Jet velocity
u: Plate velocity
Vr1 = V1 − u1 : Relative velocity of entering jet
Vr2 = V2 − u2 = KVr1 : Relative velocity of leaving jet
Vw1 = Vr1 cos β1 + u1 : Whirl velocity of entering jet
Vw2 = Vr2 cos β2 − u2 : Whirl velocity of leaving jet
V = V1 : Absolute velocity of inlet jet
V2 : Absolute velocity of leaving jet
β1 : Angle of blade at inlet
β2 : Angle of blade at outlet
α: Nozzle angle
K: Blade friction coefficient
10.2 Turbines
Hydraulic energy → mechanical energy → electrical energy
Hydraulic Power station
Reservoir → Penstock → Surge tank → Scroll casing → Guide wheel → Turbine runner
→ Daft tube
Surge tank: To reduce water hammer effect
Cavitation: Outlet of the runner or Entrance of the daft tube or suction part of pump
Wicket gate: In Francis turbine to control flow of water
Volute tube: Velocity head to pressure head
Daft tube: Convert Kinetic energy head to pressure head
Pelton Pelton2 Francis Kapplan Trend
Ns 8.5-30 30-60 60-300 300-800 Increasing
H >250 250-60 <60 Decreasing
Q . .. .... ...... Increasing
Impulse Impulse Reaction Reaction
Tangential Tangential Mixed flow Axial flow
Radial
Ns ∝ no. Adjustable
Inward
of jets Blades
flow
Pelton2 : Pelton with more than one jet
Propeller turbine: Axial √ flow Reaction turbine
N P
Specific speed Ns =
H 5/4
N is in rpm
P is in kW
H is in m
Specific speed of a turbo machine remains constant over different working conditions
Unit quantities
For a single turbine working under different conditions
N
Nu = √
H
Q
Qu = √
H
www.eggsam.com 10.2. TURBINES [136]
P
Pu = 3/2
H
For Similar turbines
Hm Hp
2 2
= 2 2
Nm Dm Np Dp
Qm Qp
3
=
Nm Dm Np Dp3
Pm Pp
3 5
= 3 5
Nm Dm Np Dp
Static pressure drop inside runner
Degree of Reaction =
Total Energy change inside runner
Vf 1
Flow ratio Ψ = √
2gH
U1
Speed ratio Φ = √
2gH
Kaplan turbine
π
Q = (Do2 − Db2 )Vf 1
4
Speed ratio for highest efficiency = 1.4 to 2
Kaplan turbine has high design efficiency and
it has constant efficiency over a wide range of
design regulation
10.3 Pumps
Foot valve + Strainer →
√ Suction pipe → impeller → Delivery pipe
N Q
Specific speed Ns =
H 3/4
Swept volume-Actual discharge
Slip S =
Swept volume
Centrifugal pump
Maximum efficiency when blades are bent backwards
www.eggsam.com 10.3. PUMPS [138]
W
Power P = (Vw2 U2 − Vw1 U1 )
g
gH H
ηm = =
Vw2 U2 He
H = Hm , manometric head
Q = πD1 B1 Vf 1 = πD2 B2 Vf 2
Vw2 U2 − Vw1 U1 Hm
Euler head He = =
g ηm
Manometric head Hm = He − loss
Vf 1
Flow ratio Ψ = √
2gH
U1
Speed ratio Φ = √
2gH
Patm,abs Pvappressure
NPSH= − hsuction − hLoss −
ρg ρg
Net positive suction head P − Pv
Thoma’s cavitation number σc = 2
=
ρV /2 ρV 2 /2
Volute casing
To collect from and deliver at constant velocity
To increase the efficiency of the pump
To reduce the loss of head in discharge
Uniform flow of fluid coming out of impeller
Forward curved blades
has obtuse angle in velocity triangle
Air vessel
-Reduces possibility of cavitation
-Smoothen the flow
-Makes it possible for the pump to run at a higher speed.
-Makes it possible to increase the suction head
-Reduces frictional loss
139
www.eggsam.com 11.1. NAMED NUMBERS [140]
Chapter 11
Heat Transfer
11.2 Conduction
Conduction happens by collision of particles and movement of electron
Diathermic: Allows heat flow.
kdiamond > kAg > kCu > kAl > kF e > kSteel
Kpure metal > kits alloy
1
kmetal ∝
T
kalloy ∝ T
1
kliquid ∝
T
knon−metal ∝ T r
T
kgas ∝
Molecular Mass
1
k ∝
M
dT
Fourier’s law of conduction: Q ∝ A
dx
dT
Q = −kA
dx
δT
Thermal resistance in conduction R =
Q
l
Plane slab R = l:Thickness of the slab
kA
ln(r2 /r1 )
Cylinder R= r2 : Outer radius
2πkl
r2 − r1
Splere R= r1 :inner radius
4πkr1 r2
Variable conductivity
If k = k0 (1 + βt)
(t2 + t1 ) (t1 − t2 )
Then Q = k0 1 + β A
2 L
1 TR2
kavg = kT dT
T2 − T1 T1
General heat conduction equations
Cartesian
coordinates
∂ ∂T ∂ ∂T ∂ ∂T ∂T
kx + ky + kz + qg = ρC
∂x ∂x ∂y ∂y ∂z ∂z ∂t
Cylindrical coordinates
∂ 2T 1 ∂T 1 ∂ 2T ∂ 2T qg ρC ∂T
2
+ + 2 2
+ 2
+ =
∂r r ∂r r ∂φ ∂z k k ∂t
Spherical coordinates
∂ 2T
1 1 ∂ ∂T 1 ∂ 2 ∂T qg ρC ∂T
2 2 2
+ 2 sinθ + 2 r + =
r sin θ ∂φ r sinθ ∂θ ∂θ r ∂r ∂r k k ∂t
1 ∂ ∂T q g ρC ∂T
r2 + =
r2 ∂r ∂r k k ∂t
www.eggsam.com 11.3. CONVECTION [142]
11.3 Convection
Newton’s law of cooling: Q ∝ A(T − T∞ )
Q = hA(T − T∞ )
δt
= P r−1/3
δ
k
Critical radius of cylinder rc =
h
2k
Critical radius of sphere rc =
h
Up to critical radius, the heat transfer rate increases
Ao − Ai
Logarithmic mean area for hollow cylinder Am =
ln(A√
o /Ai )
Logarithmic mean area, hollow sphere Am = 4πro ri = Ao Ai
Case Equations
x: distance
from T1 , Thickness
of wall = 2L
qg T2 − T1
Tx = (2L − x) + x + T1
2k 2L
Plane wall When T1 = T2 )
qg L
Twall = + T∞
h 2
qg L qL
Tmax = + + T∞
2K h
qR
Twall = + T∞
2h 2
Cylinder qg R qR
Tmax = + + T∞
qg 2 4k 2h
Tr = [R − r2 ] + Twall
4k
qR
Twall = + T∞
3h 2
Sphere qg R qR
Tmax = + + T∞
qg 2 6k 3h
Tr = [R − r2 ] + Twall
6k
www.eggsam.com 11.3. CONVECTION [143]
Forced convection
Case Equations
N ux = 0.332(Rex )1/2 (P r)1/3 (Const. temperature)
Laminar flow over flat N u = 0.664(ReL )1/2 (P r)1/3 (Const. temperature)
plates and walls N ux = 0.453(Rex )1/2 (P r)1/3 (Const. heat flux)
N u = 0.68(ReL )1/2 (P r)1/3 (Const. heat flux)
N u = 4.36 ≈ 48/11 Uniform heat flux
Laminar flow inside tubes
N u = 3.66 ≈ 48/13 Constant wall temp
Turbulent flow over flat N ux = 0.0296(Rex )0.8 P r1/3
plates N uav = 0.036(ReL )0.8 P r1/3
N uav = 0.023Re0.8 P rn
Turbulent flow in tubes n=0.3 Cooling
n=0.4 Heating
Free convection
Case Equations
General form N uav = C(GrP r)n
Vertical wall N ux = 0.378Gr1/4
Vertical plates & Laminar flow N uav = 0.59(GrP r)1/4 104 < GrP r < 109
cylinders Turbulent flow N uav = 0.10(GrP r)1/3 109 < GrP r < 1012
Laminar flow N uav = 0.54(GrP r)1/4 105 < GrP r < 2 × 107
Horizontal plates
Turbulent flow N uav = 0.14(GrP r)1/3 2 × 107 < GrP r < 3 × 1010
Laminar flow N uav = 0.53(GrP r)1/4 104 < GrP r < 109
Horizontal cylinders
Turbulent flow N uav = 0.13(GrP r)1/3 109 < GrP r < 1012
Reynold’s Analogy
Cf
St = (Assuming Pr=1)
2
Chilton and Colburn Analogy
Cf f
St.P r2/3 = =
2 8
f: Friction factor
www.eggsam.com 11.4. BOILING [144]
11.4 Boiling
Fins
Tw = Wall temperature or temperature at the base of the fin.
Ta = Ambient temperature
θ = T − Ta
θw = Tw − Ta
θ = c1 emx + c2 e−mx
θ = Arcosh mx + B sinh mx
hP
m=
kA
Actual heat transfer
Fin Effectiveness =
Heat transer without fins
Actual heat transfer Qf in
Fin efficiency η = =
Maximum heat transfer (Tw − Tα )P lh
www.eggsam.com 11.6. HEAT EXCHANGER [145]
Case Equations
h
T − Ta cosh[m(L − x)] + [sinh m(L − x)]
= km
Tw − Ta h
cosh(ml) + sinh(ml)
Finite fin km
h
√ tanh(ml) +
Q = hP kA(Tw − Ta )
km
h
1+ tanh(ml)
km
T − Ta cosh[m(L − x)]
=
Tw −√Ta cosh(ml)
Q = hP kA(Tw − Ta ) tanh(ml)
Fin with insulated tip tanh(ml)
Fin efficiency ηf in =
ml
Pl
Fin effectiveness f in = ηf in
A
ml > 2.646, tanh(ml) ≈ 1
T − Ta
= e−mx
Tw −√Ta
Infinitely long fin Q = hP kA(Tw − Ta )
1
Fin efficiency ηf in =
mlr
Pk
Fin effectiveness f in =
hA
11.7 Radiation
Irradiation (G): Rate of energy received per unit area of the surface.
Emissive power (E): Rate of energy emitted per unit area of the surface
Radiosity (J): Rate of Total energy leaving unit area of the surface
h = 6.625 × 10−34 Js
σ = 5.67 × 10−8
Grey body
α, ρ&τ same for all wavelength
α+ρ+τ =1
1 R R cosθ1 cosθ2 dA1 dA2
View factor F1−2 =
A1 A1 A2 πr2
F1−2 ⇒ From 1 to 2
Q1−2 = Q1 F1−2
F1−2 + F1−3 + F1−4 + ... = 1
A1 F1−2 = A2 F2−1
Total Emissive power E = πI = σT 4
Stefan-Boltzman law of radiation
Eb = σAT 4
E = σAT 4
E = Eb
Kirchoff ’s law
α=
www.eggsam.com 11.7. RADIATION [147]
q = E = σA(T14 − T24 )
J − Eb
Irradiation G =
1−
Radiosity J = E + ρG
σ(T14 − T24 ) Eb1 − Eb2
Q1−2,net = =
1 − 1 1 1 − 2 Rnet
+ +
A1 1 A1 F1−2 A2 2
1−
Surface resistance =
A
1
Space resistance =
A1 F12
Infinitively large parallel plates
σ(T14 − T24 )
Q1−2,net =
1 1
+ −1
1 2
Infinitely long concentric cylinders
σA1 (T14 − T24 )
Q1−2,net =
1 A1 1
+ −1
1 A2 2
Small body (1) in a large enclosure(2)
Q1−2,net = 1 σAT 4
Thermodynamics
P V γ γ=const
P T 1−γ =const
T V γ−1 =const
1 ∂V
Compressibility k = −
V ∂P T
Specific
heats
du R̄
Cv = =
dT v=c M (γ − 1)
dh
Cp =
dT p=c
Cp − nCv
Cn =
1−n
Cp
Adiabatic Index γ =
Cv
Cp − Cv = R (Meyer’s relation, ideal gases)
2
γ =1+
x
x
Internal energy u = RT [x: dof of molecule]
2
148
www.eggsam.com 12.1. BASIC PRINCIPLES [149]
Dalton’s law (also called Dalton’s law of partial pressures) states that in a mixture
of non-reacting gases, the total pressure exerted is equal to the sum of the partial pres-
sures of the individual gases.
Amagat’s law states that the volume of a gas mixture is equal to the sum of vol-
umes of the component gases, if the temperature T and the pressure p remain the same.
Gouy Stodola theorem states that rate of reversibility is proportional to the rate
of entropy generation
www.eggsam.com 12.2. WORK AND HEAT [150]
dH = mCp (T2 − T1 )
Process n W Q dS
T2
Isochoric ∞ V (P1 − P2 ) dH − W mCv ln
T1
T2
Isobaric 0 0 mCp (T2 − T1 ) mCp ln
T1
V2 V2
Isothermal 1 mRT ln W mRT ln
V1 V1
γ
Adiabatic γ (P1 V1 − P2 V2 ) 0 0
γ−1
n T2
Polytropic n (P1 V1 − P2 V2 ) dH − W mCn ln
n−1 T1
First law
• Whenever a system undergoes a cyclic change, the net work done on the system is
always equal to the amount of energy removed from the system as heat
• There exists a property of system called E such that a change in its value is equal to
the
H difference
H between the heat supplied and the work done during any change in state.
dQ = dW
• In an isolated system, the energy of the system remains constant
• Energy can neither be created nor be destroyed
• Introduces concept of internal energy
• u=q−w
• Perpetual motion machine of first kind works without power input
www.eggsam.com 12.3. LAWS OF THERMODYNAMICS [152]
Second law
• Introduces entropy
• Kelvin-Plank Statement: It is impossible to operate a cyclically operating device which
produces no effect other than the extraction of energy as heat from a single reservoir and
performs an equivalent amount of work
• Clausiu’s Statement: It is impossible to operate a cyclically operating device which
produces no effect other than the transfer of energy as heat from a low temperature body
to a high temperature body.
• PMMSK: Absorbs energy and converts the full energy to work.
Third law
• The entropy of a pure crystalline substance at absolute zero temperature is zero.
• It is impossible to achieve zero Kelvin temperature in a finite number of processes.
• The entropy measured relative to pure crystalline substance at absolute zero tempera-
ture is called absolute entropy
• Defines datum for measuring entropy
• perpectual motion machine of third kind has no friction
Carnot principle
• No heat engine operating between two given thermal reservoirs can be more efficient
than a reversible engine operating between the same two reservoirs.
• The efficiency of all reversible heat engines acting between same temperatures are same.
Turbines
Adiabatic
m(V12 − V22 )
W = H1 − H2 +
2
Diffuser
V22 − V12
H2 = H1 −
2
Throttling
• Isenthalpic process
∂T
Joule-Thomson coefficient µJT =
∂P h
Slope of constant enthalpy line in T-P diagram
µJT = 0: Inversion point
µJT > 0: Cooling
µJT < 0: Heating
Throttling of steam
Pressure decreases
Temperature increases
Dryness fraction increases
www.eggsam.com 12.4. ENTROPY [153]
Entropy increases
Enthalpy is constant
Specific volume increases
W QH − QL TH − TL
Efficiency of Carnot engine η = = =
QH QH TH
QL TL
COP of refrigerator COPr = =
QH − QL TH − TL
QH TH
COP of heat pump COPh = =
QH − QL TH − TL
12.4 Entropy
H δQ
For a system ≤ 0 (Clausius inequality)
T
H δQ
< 0 : Irreversible process
T
H δQ
= 0 : Reversible process
T
H δQ
> 0 : Impossible process
T
• Entropy is a point function and a property
• ∆Sgen (≥ 0) is not a point function
• All adiabatic processes are NOT isentropic
• All isentropic process are NOT adiabatic
• Reversible adiabatic process are isentropic
General equations
T2 V2
s2 − s1 = Cv ln + R ln
T1 V1
T2 P2
s2 − s1 = Cp ln − R ln
T1 P1
γ − n R P1
s2 − s1 = ln
γ − 1 n P2
For melting and vaporization
h1 − h2
sf g =
T
Mixing of different gases
dS = −Ru (Σni ln Xi )
Mixing same type of gas
Tf
dS = Σmi ci ln
Ti
Entropy of disorder
P2
s2 − s1 = K ln
P1
T-ds equations
T ds = du + P dv
T ds = dh − vdP
www.eggsam.com 12.5. AVAILABILITY [154]
Maxwell’s
equations
∂T ∂P
=−
∂V S ∂S V
∂P ∂S
=
∂T V ∂V
T
∂S ∂V
=−
∂P T ∂T P
∂V ∂T
=
∂S P ∂P S
12.5 Availability
• Availability is a property of system+surroundings
• High grade energy : Can be completely converted into useful work
• Low grade energy : Cannot be completely converted to work
• Exergy never increases in a process.
• Maximum possible work = Change in availability
Actual output Reversible input
Second law efficiency = or
Reversible output Actual input
η
ηII =
ηrev
T0
Available energy or exergy= 1 − δQ
T
T0
Unavailable energy or anergy= δQ = T0 ds
T
T
Exergy=mCp (T − T0 ) − mCp T0 ln
T0
Availability function for a non-flow process φ = U + P0 V − T0 S
φ2 − φ1 = (u2 − u1 ) + P0 (v2 − v1 ) − T0 (s2 − s1 )
Availability function for a flow process Ψ = H − T0 S
V 2 − V12
ψ2 − ψ1 = (h2 − h1 ) − T0 (s2 − s1 ) + 2 + g(z2 − z1 )
2
Irreversibility I = T0 ∆Su
Irreversibility I = T0 (∆Ssys + ∆Ssurroundings )
Gibbs function
G=H-TS
dG=VdP-SdT
Gibbs Helmholtz function
F=U-TS
dF=-Pdv-SdT
www.eggsam.com 12.6. PURE SUBSTANCES [155]
12.8 Psychrometry
Specific humidity
mv kg of vapour
ω= =
ma kg of dry air
Pv Pv
ω = 0.622 = 0.622
Pa Patm − Pv
Patm = Pa + Pv
Degree of saturation
ω
µ=
ωsat
Relative humidity
Pv vsat mv
φ= = =
Psat vv msat
µ
φ=
Ps
1 − (1 − µ)
Patm
Enthalpy of moist air
h = ha + ωhv
h = 1.005tdb + ω(2500 + 1.88tdb )
Cp = Cpa + ωCpv = 1.005 + 1.88ω ≈ 1.021kJ/kg-dry air-K
Tcoil − T2
Bypass factor =
Tcoil − T1
Contact factor = 1 - BPF
fg
Lewis number L =
kw cphs
For air-water mix, L=0.945
Psychrometric processes
1 : Humidification
1’: Dehumidification
2 : Sensible heating
2’: Sensible cooling
www.eggsam.com 12.9. RANKINE CYCLE [158]
Cooling tower
Approach: Temperature difference between cooled outlet water and WBT of entering air
Minimum temperature to which water can be cooled = WBT of air
Range: Temperature difference between exit cold water and incoming condenser hot water
WT − WC
Work ratio =
WT
WC
Back work ratio = =1-Work ratio
WT
Net work Wnet = WT − WC
h3 − h2
Mean temperature of heat addition Tm =
s3 − s2
T1
Thermal efficiency η = 1 −
Tm
Net work ratio is almost 1 for Rankine cycle
Rankine efficiency: Ratio of isentropic heat drop in prime mover to the amount of the
heat supplied per unit mass of steam
Rankine cycle with infinite series of regenerative feed heating has efficiency almost equal
to Carnot efficiency
Velocity compounding
Moving and fixed blades
Curtis
Pressure compounding
Moving blades and fixed nozzles
Rateau and Zoelly
V12
Nozzle efficiency ηN =
2(h0 − h1 )
Stage efficiency = ηb ηN
Reaction turbine
Q̇ = πDhVf 1
Power P = ṁ(Vw1 + Vw2 )U
2
Efficiency η = 2 −
1 + 2s cos α1 − s2
For maximum efficiency s = cos α
2 cos2 α
Maximum efficiency =
1 + cos2 α
(Vw1 + Vw2 )U
Diagram efficiency η = V 2 V 2 −V 2
1
2
+ r2 2 r1
-Adiabatic expansion
50% reaction, or identical bladings, or Parson’s turbine V1 = Vr2 V2 = Vr1 α1 =
β2 α2 = β1
P2 P3
Pressure ratio rp = =
P1 P4
V1
Compression ratio rk =
V2
V4
Expansion ratio re =
V3
www.eggsam.com 12.11. BRAYTON CYCLE, GAS TURBINES [161]
WT − WC T1 γ−1
Work ratio = = 1 − rp γ
WT T3
WC
Back work ratio = =1-Work ratio
WT
Net work Wnet = WT − WC
Wnet 1 1
Efficiency η = = 1 − γ−1 = 1 − γ−1
Qs rk
rp γ
γ
Tmax γ−1
For maximum efficiency, rpmax =
Tmin
Brayton cycle is not as efficient as Rankine cycle
Optimum work
√
rpopt = rpmax
√
T2 = T4 = √ T1 T3 √
2
Wopt = C(r T3 − T1 )
T1
η =1−
T3
Regeneration
Thermal efficiency increases
No change in Turbine work
No change in compressor work
heat supplied reduces
Mean temperature of heat addition increases
Mean temperature of heat rejection reduces
Tmin γ−1
Ideal regeneration η = 1 − rp γ
Tmax
Reheating
Turbine work increases
Thermal efficiency may or may not increase
Intercooling
Increases the net work output
Compressor work reduces
Turbine work constant
Heat supply increases
Thermal efficiency decreases
√
Perfect cooling Pi = P1 P2
12.12 Nozzle
Adiabatic, W=0
V1 <<<
p V2
V2 = v2(H1 − H2 )
u
u 2n
" n−1 #
P 2
n
V2 = t P1 v 1 1 −
n−1 P1
r
2n
V2,max = P1 v1
n−1
Turboprop engine
Has high power for take-off
Low efficiency at high altitudes and high velocity
Efficiency first increases with speed and then decreases
Flight velocity cannot exceed jet velocity
Rocket
2V0 Ve
η= 2
V0 + Ve2
Ve : Jet velocity
V0 : Vehicle velocity
Jet velocity of a rocket is independent of forward motion
Thrust=mp Ve + Ae (Pe − Pambient )
mp : mass rate of flow of propellant
Ve : Jet exit velocity
Inter-cooling/Multi-staging
" n−1 #
nN P2 nN
W = P1 (V1 − V4 ) −1
n−1 P1
Work is reduced
Weight of compressor is reduced
Flywheel weight is reduced
Volumetric efficiency is increased
In perfect inter-cooling, work is equal in both the stages
To increase efficiency
Decrease clearance ratio
Decrease delivery pressure
Multistage
Primary fuels
Wood, coal, natural gas,...
Artificial or secondary fuels
Charcoal, coal gas, coke, kerosene, diesel, petrol,...
100 8 O O
Theoretical air required= C + 8H − + S = 11.6C + 34.8(H − ) + 4.35S
23 3 8 8
kg-air/kg-fuel
11
Mass of CO2 formed = C×
3
Mass of H2 O formed = H2 × 9
www.eggsam.com 12.17. FUELS AND COMBUSTION [165]
a[O2 +3.76N2 ]
Proximate analysis
FC+VM+M+Ash=100%
Step 1: 1g at 105o C for 1 hour. Loss in weight=M
Step 2: 950o C, 7 min, covered platinum crucible. loss in weight = M+VM
Step 3: 720o C, complete burning, uncovered crucible. Remaining mass= Ash
Ultimate analysis
C+H2 +O2 +N2 +S+M+Ash=100%
Mass percentage
Orsat Apparatus
Volumetric analysis
Dry analysis
KOH (Caustic soda): CO2
Pyrogalic acid: O2
Cuprous chloride: CO
Remaining volume: N2
12.18 IC Engines
Petrol Engines(SI) Diesel Engines(CI)
Compression ratio 5-10 14-22
Speed High Relatively low
Weight Light Heavy
Peak pressure Low (To avoid self ignition) High
Thermal efficiency Low (Due to low CR) High
Cycle Otto cycle Diesel cycle
Fuel Petrol, Gasoline Diesel
Highly volatile less volatile
Ignition Spark plug is used Compression ignition
Fuel injection Air+ Fuel mixture during Fuel alone at the end of
suction stroke compression
Fuel pump and injector is
Carburetor is used
used
10% Richer flame has highest flame propagation speed and less knocking
Octane number
Iso-octane and normal heptane percentage by volume
In SI engine
BS IV: 81
BS VI: 81/85
Research Octane Number (RON) 91 91/95
Cetane number
CI engine
n-hexadecane (cetane) α methyl naphthalene
BS 4 and 6: 51
Cetane index: 46
An ignition delay in a CI engine is the time taken by the fuel to auto-ignite after
being injected into the engine cylinder.
Suction-Compression-Expansion-Exhaust
Super charging
To increase inlet air density
Uses supercharger
Best used in CI
In SI Engine, increases the chance of knocking
Ignition systems
Battery Ignition
-Battery, ignition switch, ignition coil, breaker points, condenser, distributor and spark
plugs
BP
ηmech =
IP
Brake thermal efficiency
ηrelative =
Air standard efficiency
BP
(ηth )brake =
ṁf × CV
PM EP LA N2 n
IP =
60
Actual work
Diagram factor=
Theoretical work
Vswept + Vclearance
Compression ratio =
Vclearance
Photo Chemical smog
HC and NOx
Emissions
HC follows U shaped curve (Platinum used)
NOs A shaped (Rhodium used)
CO negative exponent (Palladium used)
3 way catalytic converters
R(m−1)(n+1)p : Cm Hn Fp Clq : n + p + q = 2m + 2
R1(m−1)(n+1)p : Cm Hn Fp Clq : n + p + q = 2m
Ice plants
Reacts with copper and alloys
Requires larger displacement per TP
Ammonia Has higher compressor discharge temperature
Smell or Sulphor candle test to find leakage
Refrigerant absorber, ammonia - water
Reciprocating compressors
Window type units
Halide torch to find leakage, green flame
Freon 12
Inflammable, Non-toxic, Chemically stable
Low refrigerating effect
Low temperature cold storage
Freon 22
Halide torch to find leakage, green flame
Centrifugal systems
Freon 11
Halide torch to find leakage, green flame
Air Aircraft air conditioning
Direct contact freezing food
CO2
Reciprocating compressors
Lithium Bromide
Absorption refrigerations, solar refrigeration
Water
Azeotropes Refrigerant mix, acts like pure substance
QL TL
COP of refrigerator COPr = =
QH − QL TH − TL
Refrigerators in series acting between T1 & T and T & T2
1 1 1 1 1
= + +
COP COP1 COP2 COP1 COP2
T1 + T2
For equal work input T =
√ 2
For equal COP T = T1 T2
RSH RLH RT H
cmms,min = = =
0.0204(ti − tADP ) 50(ωi − ωADP ) 0.02(hi − hADP )
Vapor compression refrigeration
Throttle valve instead of expansion cylinder
Evaporator → Compressor → Condenser → Throttle valve
Reversed Rankine
Bell-Coleman or Joule cycle (Reverse Brayton)
Vapor absorption refrigeration system
Can use solar energy directly
Refrigerant
Low boiling point
Low freezing point
Low specific volume
Low viscosity
Low specific heat
Low positive operating pressure
High latent heat
High thermal conductivity
High density
High critical temperature
Inflammable
Non-toxic
Air conditioning
Air velocity: 6-7m/s
Air per person: 0.25 m3 /min
Sensible heat factor 0.7
Comfort conditions: 22o C and 60% RH
Air change: Air changed per hour
RSHL= 0.0204 V̇ ∆T
RLHL= 50 V̇ ∆ω
V̇ is in m3 /min
Answer is in kW
12.20 Boiler
Output of a boiler is normally stated as evaporative capacity in tonnes of steam at 100o C
that can be produced from 100o C water.
www.eggsam.com 12.20. BOILER [172]
Drum
Setting
Grate CI, above which fluid is burned
Solar Energy
Wind Energy
Biomass Energy
Geothermal (0.05 W/m2 )
Tidal energy
Ocean thermal energy conversion
Advantages
Good for environment
Unlimited supply
Cheaper
Disadvantages
Dilute form of energy
Depends on whether and location
174
www.eggsam.com 13.1. SOLAR ENERGY [175]
θ Angle of incidence Angle between incident beam and normal to the plane
θz Zenith angle Angle between the beam and normal to the horizontal plane
αa Solar altitude angle 90o − θz
γs Solar azimuth angle Angle between projection of solar beam on the horizontal and
meridian
Angular measurement of time
ω Hour angle Angle of rotation of earth since solar noon
-180o to 180o (Morning-evening)
Angle between collector plate and horizontal
β Slope of the plane 0o − 180o
(0-towards equator-90-Away from equator-180)
Vertical surface, β = 90o
cos θ = sin φ(sin δ cos β +cos δ cos γ cos ω sin β)+cos φ(cos δ cos ω cos β −sin δ cos γ sin β)+
cos δ sin γ sin ω sin β
Solar constant ISC = 1367W/m2
n: day ofthe year
360
In = ISC 1 + 0.033 cos n
365
360
Declination angle δ = 23.45o sin × (284 + n) (Cooper’s relation)
365
On March 21 ans september 21, δ ≈ 0
2hC 2
Spectral intensity Iλ (λ, T ) = Plank’s law
hC
λ5 exp −1
λkT
k: Universal Boltzman constant
C: Speed of light
h: Plank’s const
Absorber plate:Absorbs solar radiation and transfers heat to the fluid (Cu, Al, Brass,
steel, silver)
Riser tube: To absorb maximum heat from from absorber plate and transfer it to the fluid.
www.eggsam.com 13.1. SOLAR ENERGY [177]
Transparent cover: To reduce re-radiation. Uses glazed glass. High transmissivity for
for short waves and high reflectivity for long waves.
τα
(τ α)net =
1 − (1 − α)ρ
Useful energy gien to the fluid
Collector efficiency=
Total energy reaching the collector
Effective area of apperture
Concentration ratio =
Area of collector
Adding more layers of covers reduces τ α product. Reduces heat loss through convec-
tion. Reduces efficiency.
Parabolic collectors
Line concentrating
Tracing mechanism is necessary (In one axis)
Concentration ratio ≈ 100
Up to 300o C
Thermo-Chemical storage
Endothermic and exothermic reactions to store and extract energy
Thermo-Chemical Method
Biomass Gasification
Drying → Pyrolysis → Oxidation → Reduction
Produces Producer gas: CO2 , CO, N2 , H2
Partial combustion of biomass
Hydrogen-Oxygen FC Alkaline FC
Hydrogen-air FC Phosphoric acid FC
Hydrazine-oxygen FC Polymer electrolytic membrane FC
Hydrocarbon-air FC Molten carbonate FC
Synthesys gas-air FC Solid oxide FC
Ammonia-air FC
Hydrogen-Oxygen fuel cell
Fuel: Hydrogen (At negative electrode, Anode)
Oxygen (At positive electrode, Cathode)
At anode: H2 → 2H + + 2e−
1
At cathode: 2H + + 2e− + O2 → H2 O
2
Electrolyte
Conductive to ions and nonconductive to electricity
Should not get charged
Electrode
Good conductor of electricity
Should not corrode on contact with electrolyte
Stable at high temperatures
Thank you
.
I would like to thank everyone for helping me make this book a reality. I extend my
sincere thanks to all my social media supporters. The constant support of them
motivated me in making this book a reality.I especially thank the following people who
helped me by pointing out errors and mentioning the missing equations.
Anshuman Sarma
Er Yogendra Pratap Singh
Satya Panda
Hemant Sukhija
Raja Radha
Malaya ranjan
Arjun Bajpayee
Deepak Kumar
Somesh Akolkar
3 Anonymous contributors
181