Handbook - Mechanical
Handbook - Mechanical
Handbook - Mechanical
Mechanical Engineering
HAND BOOK
2020
Department of
Mechanical & Manufacturing Engineering
Table of Contents
2) Let ∅ = ∅(𝑥, 𝑦, 𝑧) is the scalar field, the gradient of ∅ at the point (𝑥, 𝑦, 𝑧) is
𝜕∅ 𝜕∅ 𝜕∅ 𝜕∅
𝑔𝑟𝑎𝑑∅ = ∇∅ = 𝜕𝑥 𝒊 + 𝜕𝑦 𝒋 + 𝜕𝑧 𝒌 = ∑ 𝜕𝑥 𝒊
3) Let 𝑓⃗ = 𝑓1 𝒊 + 𝑓2 𝒋 + 𝑓3 𝒌 is the vector field, the divergence of 𝑓⃗ at the point (𝑥, 𝑦, 𝑧) is,
𝜕𝒇 𝜕𝒇 𝜕𝒇 𝜕𝒇
𝑑𝑖𝑣𝑓⃗ = ∇. ⃗⃗
𝒇 = 𝜕𝑥𝟏 + 𝜕𝑦𝟐 + 𝜕𝑧𝟑 = ∑ 𝜕𝑥𝟏
4) If 𝜙1 (𝑥, 𝑦, 𝑧) and 𝜙2 (𝑥, 𝑦, 𝑧) are the two surfaces, then the angle between their surfaces at
(𝑥1 , 𝑦1 , 𝑧1 ) is
∇𝜙 .∇𝜙
cos 𝜃 = |∇𝜙 1||∇𝜙2 |.
1 2
5) If 𝜙(𝑥, 𝑦, 𝑧) is a scalar function and 𝑑⃗ is the given direction vector, then the directional derivative
𝑑⃗
of 𝜙 along 𝑛̂ is ∇𝜙. 𝑛̂, where 𝑛
̂ = |𝑑⃗|.
7) The directional derivative of a scalar function 𝜙 at any point is maximum along ∇𝜙.
𝜕2 𝜕2 𝜕2
8) The Laplacian operator ∇2 is defined as ∇2 = 𝜕𝑥 2 + 𝜕𝑦 2 + 𝜕𝑧 2
𝜕2 ∅ 𝜕2 ∅ 𝜕2 ∅
9) The Laplacian of a scalar function ∅ as, ∇2 ∅ = 𝜕𝑥 2 + 𝜕𝑦 2 + 𝜕𝑧 2
10) Let 𝑓⃗ = 𝑓1 𝒊 + 𝑓2 𝒋 + 𝑓3 𝒌 is the vector field, the curl of 𝑓⃗ at the point (𝑥, 𝑦, 𝑧) is
𝑖 𝑗 𝑘
𝜕 𝜕 𝜕
𝑐𝑢𝑟𝑙 𝑓⃗ = ∇ × 𝑓⃗ = | |
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝑓1 𝑓2 𝑓3
11) If ∇ × 𝑓⃗ = 𝟎, then the vector field 𝑓⃗ is irrotational
12) If ∇. 𝑓⃗ = 𝟎, then the vector field 𝑓⃗ is solenoidal
Vector Integration
1) Suppose r(u) = x(u)𝐢 + y(u)𝐣 + z(u)𝐤 is the position vector of points P(x,y,z) and suppose
r(u) defines a curve C joining points 𝑃1 and 𝑃2 where 𝑢 = 𝑢1 and 𝑢 = 𝑢2 respectively. Let
𝐀(x, y, z) = A1 𝐢 + A2 𝐣 + A3 𝐤 be a vector function of position defined and continuous along C.
P2
2) If A is the force F on a particle moving along C, then the line integral represents the work done
by the force Work done= F .dr
C
3) If C is a closed curve, the integral around C is often denoted by
A dr A dx A dy A dz
1 2 3
4) In Aerodynamics and fluid dynamics, this line integral is called the circulation of A about C, where
A represents the velocity of a fluid.
5) Let S be a two- sided surface. A unit normal n to any point of the positive side of S is called a
positive or outward drawn unit normal. Associate with the differential of surface area dS a vector
dS and whose direction is that of n. Then dS = n dS. The integral A. dS A ndS is
S S
called surface integral or the flux of A over S.
dxdy
6) Suppose that the surface S has projection R on the xy-plane, then A dS A n n k
S R
A dr ( A) ndS ( A) dS
C S S
where C is traversed in the positive direction.
10) Green’s Theorem in the plane: Suppose R is a closed region in the xy-plane bounded by a
simple closed curve C, and suppose M and N are continuous functions of x and y having
N M
continuous derivatives in R. Then Mdx Ndy x
C R
y
dxdy where C is traversed in
the positive (counter- clockwise) direction.
FOURIER SERIES
A function 𝑓(𝑡) is said to be periodic if there exists a real number T such that 𝑓(𝑡 + 𝑇) = 𝑓(𝑡) for all t.
Smallest such number is called period of the function.
Note: If 𝑓(𝑡) is periodic with period 𝑇, then 𝑓(𝑡 + 𝑛𝑇) = 𝑓(𝑡) for all integers n.
If 𝑓(𝑥) is a periodic function with period 2l and is known in the interval 𝛼 < 𝑥 < 𝛼 + 2𝑙 , then the Fourier
series expansion of f(x) in terms of an infinite sum of sines and cosines is given by
𝑛𝜋𝑥 𝑛𝜋𝑥
𝑓(𝑥) = 𝑎0 + ∑∞
𝑛=1 𝑎𝑛 cos + ∑∞
𝑛=1 𝑏𝑛 sin
𝑙 𝑙
where
1 𝛼+2𝑙 1 𝛼+2𝑙 𝑛𝜋𝑥 1 𝛼+2𝑙 𝑛𝜋𝑥
𝑎0 = 2𝑙 ∫𝛼 𝑓(𝑥) 𝑑𝑥, 𝑎𝑛 = 𝑙 ∫𝛼 𝑓(𝑥) cos 𝑑𝑥, 𝑏𝑛 = 𝑙 ∫𝛼 𝑓(𝑥) sin 𝑑𝑥
𝑙 𝑙
Note: The graph of an even function is symmetric about 𝑦 −axis, whereas the graph of an odd function is
symmetric about the origin.
Consider a piecewise continuous function 𝑓(𝑥), defined in a finite interval (0, 𝑙). Then it is possible to extend
𝑓(𝑥) to a periodic function, which is even or odd.
𝑓(𝑥), 0<𝑥<𝑙
𝑔(𝑥) = { ; 𝑔(𝑥 + 2𝑙) = 𝑔(𝑥).
𝑓(−𝑥), −𝑙 < 𝑥 < 0
Then 𝑔(𝑥) is called an even periodic extension of 𝑓(𝑥).
FOURIER TRANSFORMS
Consider the Fourier integral representation of the function 𝑓(𝑥) given by
1 ∞ ∞ 1 ∞ 1 ∞
𝑓(𝑥) = ∫ ∫ 𝑓(𝑡)𝑒 𝑖𝑠(𝑥−𝑡) 𝑑𝑡
√2𝜋 −∞ −∞
𝑑𝑠 = ∫ 𝑒 𝑖𝑠𝑥 [√2𝜋 ∫−∞ 𝑓(𝑡)𝑒 −𝑖𝑠𝑡 𝑑𝑡] 𝑑𝑠
√2𝜋 −∞
1 ∞
Let 𝐹(𝑠) = ∫ 𝑓(𝑡)𝑒 −𝑖𝑠𝑡 𝑑𝑡
√2𝜋 −∞
…………………………………………………………………(1)
1 ∞
Then 𝑓(𝑥) = ∫ 𝐹(𝑠)𝑒 −𝑖𝑠𝑥 𝑑𝑠…………………………………………………………………(2)
√2𝜋 −∞
The integral defined by (1) is called the Fourier transform of the function 𝑓(𝑥) and is denoted by 𝐹{𝑓(𝑥)} .
Given 𝐹(𝑠) = 𝐹{𝑓(𝑥)} , the formula (2) defined 𝑓(𝑥), which is called the inverse Fourier transform of
𝐹(𝑠)and is denoted by 𝐹 −1 {(𝐹(𝑠)}.
Note: A function 𝑓(𝑥) is said to be self-reciprocal under Fourier transforms if 𝐹{𝑓(𝑥)} = 𝑓(𝑠).
Properties of the Fourier transforms:
If 𝐹{𝑓(𝑥)} = 𝐹(𝑠), then
1. 𝐹{𝑒 𝑖𝑎𝑥 𝑓(𝑥)} = 𝐹(𝑠 − 𝑎).
2. 𝐹{𝑓(𝑥 − 𝑎)} = 𝑒 𝑖𝑠𝑎 𝐹(𝑠).
1 𝑠
3. 𝐹{𝑓(𝑎𝑥)} = 𝑎 𝐹 (𝑎) , 𝑎 > 0.
Convolution:
For functions 𝑓(𝑥) & 𝑔(𝑥) , we define the convolution of 𝑓(𝑥) & 𝑔(𝑥) denoted
1 ∞
by (𝑓 ∗ 𝑔)(𝑥) as (𝑓 ∗ 𝑔)(𝑥) = ∫−∞ 𝑓(𝑡)𝑔(𝑥 − 𝑡)𝑑𝑡 ;provided the integral exists.
√2𝜋
Note that 𝑓 ∗ 𝑔 = 𝑔 ∗ 𝑓.
Convolution Theorem:
Parseval’s Identity:
∞ ∞
If 𝐹{𝑓(𝑥)} = 𝐹(𝑠), then ∫−∞| 𝑓(𝑥)|2 𝑑𝑥 = ∫−∞|𝐹(𝑠)|2 𝑑𝑠 .
Fourier Cosine and Sine transforms
2 2 ∞
= √𝜋 [√𝜋 ∫0 𝑓(𝑡) cos 𝑠𝑡 𝑑𝑡] cos 𝑠𝑥 𝑑𝑠
2 ∞
Let 𝐹𝑐 { 𝑓(𝑥)} = √𝜋 ∫0 𝑓(𝑡) cos 𝑠𝑡 𝑑𝑡 = 𝐹𝑐 (𝑠)……………………………………………..(A)
2 ∞
Then 𝑓(𝑥) = √𝜋 ∫0 𝐹𝑐 { 𝑓(𝑥)} cos 𝑠𝑥 𝑑𝑠 ……………………………………………………(B)
The transform 𝐹𝑐 {𝑓(𝑥)} defined by (A) is called the Fourier cosine transform of f(x).The formula (B) is
called the inverse Fourier cosine transform of 𝐹𝑐 {𝑓(𝑥)} = 𝐹𝑐 (𝑠) and is denoted by 𝑓(𝑥) = 𝐹𝑐−1 {𝐹𝑐 (𝑠)}.
2 ∞
𝑓(𝑥) = 𝐹𝑠−1 {𝐹𝑠 (𝑠)} = √𝜋 ∫0 𝐹𝑠 (𝑠) sin 𝑠𝑥 𝑑𝑠.
Definition: A function f(x) is said to be self-reciprocal under Fourier cosine (sine) transform if
2
𝐹𝑐 {𝑓 ′ (𝑥)} = −√𝜋 𝑓(0) + 𝑠𝐹𝑠 (𝑠) and 𝐹𝑠 {𝑓 ′ (𝑥)} = −𝑠𝐹𝑐 (𝑠)
2 2
(6) 𝐹𝑐 {𝑓 ′′ (𝑥)} = −√𝜋 𝑓 ′ (0) − 𝑠 2 𝐹𝑐 (𝑠) and 𝐹𝑠 {𝑓 ′′ (𝑥)} = √𝜋 𝑠𝑓 (0) − 𝑠 2 𝐹𝑠 (𝑠) provided 𝑓(𝑥)
and 𝑓′(𝑥) → 0 as 𝑥 → ∞.
𝑑𝐹𝑠 𝑑𝐹𝑐
(7) 𝐹𝑐 {𝑥𝑓(𝑥)} = and 𝐹𝑠 {𝑥𝑓(𝑥)} =− .
𝑑𝑠 𝑑𝑠
(8) If 𝐹𝑐 {𝑓(𝑥)} = 𝐹𝑐 (𝑠), 𝐹𝑐 {𝑔(𝑥)} = 𝐺𝑐 (𝑠), 𝐹𝑠 {𝑓(𝑥) = 𝐹𝑠 (𝑠) and 𝐹𝑠 {𝑔(𝑥)} = 𝐺𝑠 (𝑠) exist, then
∞ ∞ ∞
∫0 𝐹𝑐 (𝑠)𝐺𝑐 (𝑠)𝑑𝑠 = ∫0 𝐹𝑠 (𝑠)𝐺𝑠 (𝑠)𝑑𝑠 = ∫0 𝑓(𝑥)𝑔(𝑥)𝑑𝑥 and
∞ ∞ ∞
∫0 |𝐹𝑐 (𝑠)|2 𝑑𝑠 = ∫0 |𝐹𝑠 (𝑠)|2 𝑑𝑠 = ∫0 |𝑓(𝑥)|2 𝑑𝑥 which is called Perceval’s Identity.
Partial differential equation is an equation involving partial derivatives of an unknown functions of two or
more independent variables.
The degree of the Partial differential equation is the degree of highest derivative after clearing the
fractional power.
Partial differential equation is linear if it is of the first degree in the dependent variable and its partial
derivatives.
2u 2 u
2 One dimensional Homogeneous wave equation
c
t 2 x 2
u 2u One dimensional Homogeneous heat equation
c2 2
t x
2u 2u Two dimensional Homogeneous Laplace equation
0
x 2 y 2
2u 2u Two dimensional non Homogeneous with (f≠ 0) Poisson equation
f ( x, y )
x 2 y 2
2u 2u 2u Three dimensional Laplace equation
0
x 2 y 2 z 2
Then, the solution for the p.d.e in is in the form 𝑢(𝑥, 𝑡) = 𝑋(𝑥) . 𝑇(𝑡) where X and T are the function of x
and t respectively.
Substitute the partial derivatives involved in the equation and separate the variables, and equate to a
common constant 𝑐 2 .
2 y 2 y T
1. One-dimensional Wave equation: c2 where c2
t 2 x 2 m
T-tension of the string, m-mass per unit length of the string.
u 2u k
3. One-dimensional heat-flow equation: c2 2 where c
2
.
t x s
-density, s-the specific heat, k-thermal conductivity.
4. Solution of one-dimensional heat-flow equation:
Numerical Methods
Numerical Solution of Boundary Value Problems:
Boundary value problems are of great importance in science in engineering. In this section, we have
elaborated Numerical methods based on finite difference scheme for the solution of following problems:
1. Boundary value problems in second order ordinary differential equations
2. Boundary value problems governed by linear second order partial differential equations:
Laplace equation and Poisson equation.
3. Initial boundary value problems governed by linear second order partial differential
equations: One dimensional Heat and Wave equation
For our discussion, we shall consider only the linear second order ordinary differential equations
y p(x) y(x) q(x) y r(x) , x [a, b] (1)
Since the ordinary differential equation is second order, we need to prescribe two conditions to obtain a
unique solution of the problem. If the conditions are prescribed at the end points x = a and x = b, then it is
called a two point boundary value problem. The two conditions required to solve (1), can be prescribed in
one of the following three boundary conditions:
Approximations to y(x i ) at x x i :
y(x i )
1
y(x i1 ) - y(x i ) , or yi
1
yi1 - yi (5)
h h
(ii) Backward difference approximation of first order or O(h) approximation:
y(x i )
1
y(x i ) - y(x i-1 ) , or yi
1
yi - yi-1 (6)
h h
(iii) Central difference approximation of second order or O(h2) approximation:
y(x i )
1
y(x i1 ) - y(x i-1 ) , or yi
1
yi1 - yi-1 (7)
2h 2h
Approximations to y(x i ) at x x i :
The finite difference solution of a boundary value problem (1) is obtained by replacing the differential
equation at each nodal point by a difference equation ((5) – (8)) along with given boundary conditions ((2)
– (4))
Boundary value problems governed by linear second order partial differential equations:
Over a two dimensional Cartesian domain, let 𝑢 be the dependent variable. Then a general second order
partial differential equation may be written as,
𝜕2 𝑢 𝜕2 𝑢 𝜕2 𝑢
𝐴 𝜕𝑥 2 + 𝐵 𝜕𝑥𝜕𝑦 + 𝐶 𝜕𝑦 2 + 𝐹(𝑥, 𝑦, 𝑢, 𝑢𝑥 , 𝑢𝑦 ) + 𝐺 = 0 ----------------------(1)
where 𝐴, 𝐵, 𝐶 are functions of (𝑥, 𝑦) and 𝐹 may be non-linear function, then equation (1) is called a quasi-
linear partial differential equation. If 𝐹 is also a linear function then equation (1) is called a linear partial
differential equation. If 𝐺 = 0 then equation (1) is homogeneous otherwise non-homogeneous.
Equation (1) is said to be elliptic, parabolic and hyperbolic depending on 𝐵2 − 4𝐴𝐶 < 0, 𝐵2 − 4𝐴𝐶 = 0
and 𝐵2 − 4𝐴𝐶 >, 0 at a point or in a domain.
𝜕𝑢 1
Backward: ( ) = (𝑢𝑖,𝑗 − 𝑢𝑖−1,𝑗 ) + 𝑂(ℎ)
𝜕𝑥 𝑖,𝑗 ℎ
𝜕𝑢 1
Central: ( ) = (𝑢𝑖+1,𝑗 − 𝑢𝑖−1,𝑗 ) + 𝑂(ℎ2 )
𝜕𝑥 𝑖,𝑗 2ℎ
Similarly,
𝜕𝑢 1
Forward: ( ) = 𝑘 (𝑢𝑖,𝑗+1 − 𝑢𝑖,𝑗 ) + 𝑂(𝑘)
𝜕𝑦 𝑖,𝑗
𝜕𝑢 1
Backward: ( ) = (𝑢𝑖,𝑗 − 𝑢𝑖,𝑗−1 ) + 𝑂(𝑘)
𝜕𝑦 𝑖,𝑗 𝑘
𝜕𝑢 1
Central: ( ) = (𝑢𝑖,𝑗+1 − 𝑢𝑖,𝑗−1 ) + 𝑂(𝑘 2 )
𝜕𝑦 𝑖,𝑗 2𝑘
The finite difference approximations for the second order partial derivatives are as follows.
𝜕 2𝑢 1
( 2 ) = 2 (𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗 ) + 𝑂(ℎ2 )
𝜕𝑥 𝑖,𝑗 ℎ
𝜕2 𝑢 𝜕2 𝑢
Similar expressions can be written for ( ) , (𝜕𝑥𝜕𝑦).
𝜕𝑦 2
𝜕2 𝑢 𝜕2 𝑢
+ 𝜕𝑦 2 = 𝑓(𝑥, 𝑦) or ∇2 𝑢 = 𝑓(𝑥, 𝑦), 𝑎 ≤ 𝑥 ≤ 𝑏, 𝑐 ≤ 𝑦 ≤ 𝑑
𝜕𝑥 2
1
𝑢𝑖,𝑗 = (𝑢 + 𝑢𝑖+1,𝑗 + 𝑢𝑖,𝑗+1 + 𝑢𝑖,𝑗−1 )
4 𝑖−1,𝑗
Diagonal 5-point formula:
1
𝑢𝑖,𝑗 = (𝑢 + 𝑢𝑖+1,𝑗−1 + 𝑢𝑖+1,𝑗+1 + 𝑢𝑖−1,𝑗−1 )
4 𝑖−1,𝑗+1
Solution to Poisson equation:
where 𝑐 2 is the diffusivity of the substance, 𝑢(𝑥, 𝑡) is a temperature function which is defined for values of
𝑥 from 0 to 𝑙(length of the bar)and for values of time 𝑡 from 0 to ∞.
Schmidt method:
kc2 1
ui,j+1 = λui−1,j + (1 − 2λ)ui,j + λui+1,j where λ = ;0<λ≤2
h2
Where λ is called the mesh ratio parameter, h mesh length along x-axis and k mesh length
along t-axis
1
ui,j+1 = (u + ui+1,j )
2 i−1,j
1
Note : For λ = 6, the error in Schmidt’s formula is least
Consider the problem of vibrations of an elastic string governed by the one dimensional wave equation:
𝜕 2𝑢 𝜕 2𝑢
= 𝑐 2
, 0xl , t0
𝜕𝑡 2 𝜕𝑥 2
Subject to the initial and boundary conditions:
(we consider the case when the ends of the string are fixed)
Where c2 is a constant and depends on the material properties of the elastic string.
𝑐 2𝑘2
𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1 = [𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗 ]
ℎ2
The above method is stable when r ≤ 1, Higher order method uses the value of r = 1
Therefore, we have 𝑢𝑖,𝑗+1 = [𝑢𝑖+1,𝑗 + 𝑢𝑖−1,𝑗 − 𝑢𝑖,𝑗−1 ], is the explicit finite difference
method for wave equation at higher level.
1
To get the value at first level, we have 𝑢𝑖,1 = [𝑢𝑖−1,0 + 𝑢𝑖+1,0 ] + k𝑔𝑖 .
2
𝜕𝑢
If the initial condition is prescribed as (𝑥, 0) = 0, that is g(x) = 0, then above formula simplifies
𝜕𝑡
to
1
𝑢𝑖,1 = [𝑢𝑖−1,0 + 𝑢𝑖+1,0 ] , gives the value of ‘u’ at first level.
2
1. Basic Concepts
1.1 Relation between the number of links and the number of joints for a kinematic chain having
constrained motion
2. Velocity Analysis
2.1. Addition and Subtraction of Vectors
3. Acceleration Analysis
3.1. Tangential Acceleration, 𝑓 𝑡 (𝑚/𝑠 2 )
𝑓 𝑡 = 𝛼. 𝑟 𝛼 = Angular Acceleration of link (rad/s2)
r = length of link (m)
4.2 Length of path of contact (LPOC), length of path of approach (LPOA) and length of path of
recess (LPOR)
2
LPOA = √𝑅𝑎2 − (𝑅2 cos ∅)2 − (𝑅2 sin ∅) 𝑅𝑎1 and 𝑅𝑎2 = radii of addendum circle on
2 gear 1 and gear 2 respectively (mm)
LPOR = √𝑅𝑎1 − (𝑅1 cos ∅)2 − (𝑅1 sin ∅)
𝑅1 and 𝑅2 = radii of pitch circle on gear 1 and
LPOC = LPOA + LPOR gear 2 (mm)
∅ = pressure angle (degrees)
4.3 Length of arc of contact (LAOC), length of arc of approach (LAOA) and length of arc of recess
(LAOR)
4.6 Minimum number of teeth on the pinion to avoid interference when meshing with a rack
5. Gear Trains
5.1. Train Value
1
𝑆𝑝𝑒𝑒𝑑 𝑟𝑎𝑡𝑖𝑜 =
𝑡𝑟𝑎𝑖𝑛 𝑣𝑎𝑙𝑢𝑒
6.12 Maximum velocity in the belt (𝑣𝑚𝑎𝑥 ) for maximum power transmission
7. Friction
7.1 Frictional torque on a flat pivot bearing
(𝑟𝑜 + 𝑟𝑖 )
𝑇 = 𝜇𝑊
2
1.1. Grain Fineness Number (GFN) where, Mi = multiplying factor for the ith sieve
fi = percentage sand retained on the ith sieve
𝐺𝐹𝑁 = ∑ 𝑀𝑖 𝑓𝑖 and the pan
2. Lathe operations
Figure 2.1: Schematic of calculating the swivel angle for taper turning
2.3.2. Time for a single pass where t = Time required for single pass (min)
L = Length of the job (mm)
𝐿 + 𝐿𝑂 LO = Over travel of the tool beyond the length
𝑡= of the job (mm)
𝑓𝑁
f = Feed rate (mm/rev)
N = Rotational speed of the workpiece (rpm)
2.3.3 Number of roughing passes where Pr = Number of passes for roughing operation
A = Total machining allowance (mm)
𝐴 − 𝐴𝑓 Af = Finish machining allowance (mm)
𝑃𝑟 =
𝑑𝑟 dr = Depth of cut in roughing (mm)
2.3.4 Number of finishing passes where Pf = Number of passes for finishing operation
Af = Finish machining allowance (mm)
𝐴𝑓
𝑃𝑓 = df = Depth of cut in finishing (mm)
𝑑𝑓
3. Milling operations
𝑁×|𝑁1 −𝑁2|
40×𝑁1 ×𝑁2
= 𝑛1 where, N1 and N2 are the hole circles in
the same plate
40 𝑛 𝑛
= ±
𝑁 𝑁1 𝑁2
Table 4.1: Mean cutting process parameter for drilling using HSS Tool
4.1.3 Total length of tool travel where, L = Total length of tool travel (mm)
𝐿 = 𝐴+𝑙+2 l = length of the hole (mm)
Crystal Structure
Phases in solids
Lever Rule:
S = Solidus line
L = Liquidus line
F = Fulcrum
Ws = Weight percent of solid
Wl = Weight percent of liquid
𝐹𝐿
Ws fomed = 𝑋 100 (Proeutectic Formed)
𝑆𝐿
𝐹𝑆
Wl remaining = 𝑋 100 (Wt. % which will convert into eutectic mixture)
𝑆𝐿
Weight % of eutectic mixture formed: Wt. percent of liquid present at eutectic temperature before
the eutectic reaction begins = 100 - Ws
𝐹𝐷
Weight % of A or α formed = 𝑋 100
𝐶𝐷
𝐶𝐹
Weight % of B or β formed = 𝑋 100
𝐶𝐷
𝐶𝐹 𝐹𝐷
Weight ratio of two solids after the reaction: A/B= or B/A =
𝐹𝐷 𝐶𝐹
Equilibrium Condition:
Load (N)
P P1 P2 l1 , l2 = Elongation of respective
material (mm)
Elongation (mm)
l l1 l1 P1 , P2 = Load of respective material (N)
1.15.2 Thermal stress developed when the bar is t = Temperature rise (OC)
fully constrained (MPa): = Linear coefficient of expansion (OC)
th Et E = Young’s Modulus (MPa)
1.15.3 Thermal stress developed when the bar is E = Young’s Modulus (MPa)
partially constrained (MPa): L = Length of the member (mm)
E t = Temperature rise (OC)
th (tL ) = Linear coefficient of expansion (OC)
L
= Restricted elongated length of the
member (mm)
2E
6G
L2
M w
2
4. Deflection of beams
4.1 Differential equation of the beam: E = Young’s Modulus (MPa)
d2y I = Moment of inertia of the
EI M cross-section (mm4)
dx 2
M = Bending moment (Nmm)
4.2 Cantilever beam with moment at free end M O = Externally applied moment
M 0L (Nmm)
Maximum slope (in radians) = At the free end
EI L = Length of the beam (mm)
2
MoL W = Applied Load (N)
Maximum deflection (in mm) = At the free end
2 EI w = Distributed load (N/mm)
WL2
Maximum slope (in radians) = At the free end
2 EI
WL3
Maximum deflection (in mm) = At the free end
3EI
wL3
Maximum slope (in radians) = At the free end
6 EI
wL4
Maximum deflection (in mm) = At the free end
8EI
wL3
Maximum slope (in radians) = At the free end
24 EI
wL4
Maximum deflection (in mm)= At the free end
30 EI
WL2
Maximum slope (in radians) = At the support
16 EI
WL3
Maximum deflection = At the centre
48EI
wL3
Maximum slope (in radians) = At the support
24 EI
5wL4
Maximum deflection (in mm) = At the centre
384 EI
wL4
Maximum deflection (in mm) = 0.006523
EI
At a distance of 0.5193 L from the end with load
intensity zero.
5. Torsion of shafts
5.1 Torsion equation: T = Torque (Nmm)
T G J = Polar Moment of inertia (mm4)
J r L r = Radius of the shaft (mm)
= Shear stress corresponding to radius (MPa)
G = Modulus of rigidity (MPa)
= Angular twist (radians)
L = Length of the shaft (mm)
Crushing Rankine’s
Material
stress (MPa) constant
Mild steel 320 1/7500
Cast iron 550 1/1600
Wrought iron 250 1/9000
Timber 50 1/750
c
pd d = Change in diameter (mm)
2t d = Initial diameter (mm)
7.1.2 Longitudinal stress (MPa) l = Change in length (mm)
pd
l l = Initial length (mm)
4t
V = Change in volume (mm3)
7.1.3 Circumferential strain
V = Original Volume (mm3)
d pd
c (2 ) t = Thickness of the cylinder/sphere (mm)
d 4tE
E = Young’s Modulus (MPa)
7.1.4 Longitudinal strain = Poisson’s ratio
l pd
l (1 2 )
l 4tE
7.1.5 Volumetric strain
V pd
V (5 4 )
V 4tE
8. Compound stresses
8.1 Normal Stress on plane AB (MPa) x = Direct stress along X direction (MPa)
x y x y y = Direct stress along Y direction (MPa)
n cos 2 sin 2
2 2
= Shear stress (MPa)
8.2 Tangential Stress on plane AB (MPa)
= angle of reference plane
x y
t sin 2 cos 2
2
8.3 Major & Minor Principal Stresses (MPa)
y x y
2
1, 2 x
2
2 2
R n2 t2
8.5 Resultant stress inclination
n
tan
t
8.6 Inclination of major principal plane
2
tan 21
x y
8.7 Inclination of minor principal plane
2 1 90
8.8 Maximum tangential Stress (MPa)
1 2
max
2
8.9 Inclination of plane of maximum tangential stress
3 1 45
2
𝑣2
𝑤1−2 = ∫ 𝑃𝑑𝑣 = 𝑝1 𝑣1 ln
1 𝑣1
2 n = Polytropic index
𝑃1 𝑣1 − 𝑃2 𝑣2
𝑤1−2 = ∫ 𝑃𝑑𝑣 =
1 𝑛−1
2 𝛾 = 𝐴𝑑𝑖𝑎𝑏𝑎𝑡𝑖𝑐 𝑖𝑛𝑑𝑒𝑥
𝑃1 𝑣1 − 𝑃2 𝑣2
𝑤1−2 = ∫ 𝑃𝑑𝑣 =
1 𝛾−1
4 Entropy
.
4.1 General formula for change in entropy
2 𝑑𝑞𝑅 = Heat transfer during process
𝑑𝑞𝑅
𝑠2 − 𝑠1 = ∫ T= Temperature (K)
1 𝑇 1 & 2 represents initial and final state
4.2 When two fluids are mixed together change in entropy of the universe
𝑇𝑓 𝑇𝑓 (∆𝑆)𝑈𝑛𝑖 = 𝐶ℎ𝑎𝑛𝑔𝑒 𝑖𝑛 𝑒𝑛𝑡𝑟𝑜𝑝𝑦 𝑜𝑓 𝑢𝑛𝑖𝑣𝑒𝑟𝑠𝑒 (𝐽/𝑘𝑔𝐾)
(∆𝑆)𝑈𝑛𝑖 = 𝑚1 𝑐1 ln + 𝑚2 𝑐2 ln 𝑚1 = 𝑚𝑎𝑠𝑠 𝑜𝑓 𝑡ℎ𝑒 𝑓𝑖𝑟𝑠𝑡 𝑓𝑙𝑢𝑖𝑑 (𝑘𝑔)
𝑇1 𝑇2 𝑐1 = 𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐 ℎ𝑒𝑎𝑡 𝑜𝑓 𝑡ℎ𝑒 𝑓𝑖𝑟𝑠𝑡 𝑓𝑙𝑢𝑖𝑑 (𝐽/𝑘𝑔𝐾)
𝑚2 = 𝑚𝑎𝑠𝑠 𝑜𝑓 𝑡ℎ𝑒 𝑠𝑒𝑐𝑜𝑛𝑑 𝑓𝑙𝑢𝑖𝑑 (𝑘𝑔)
𝑐2 = 𝑠𝑝𝑒𝑐𝑖𝑓𝑖𝑐 ℎ𝑒𝑎𝑡 𝑜𝑓 𝑡ℎ𝑒 𝑠𝑒𝑐𝑜𝑛𝑑 𝑓𝑙𝑢𝑖𝑑 (𝐽/𝑘𝑔𝐾)
𝑇1 = 𝐼𝑛𝑖𝑡𝑖𝑎𝑙 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑒 𝑜𝑓 𝑓𝑖𝑟𝑠𝑡 𝑓𝑙𝑢𝑖𝑑 (𝐾)
𝑇2 = 𝐼𝑛𝑖𝑡𝑖𝑎𝑙 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑒 𝑜𝑓 𝑠𝑒𝑐𝑜𝑛𝑑 𝑓𝑙𝑢𝑖𝑑 (𝐾)
𝑇𝑓 = 𝐸𝑞𝑢𝑖𝑙𝑖𝑏𝑟𝑖𝑢𝑚 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑒 𝑜𝑓 𝑚𝑖𝑥𝑡𝑢𝑟𝑒 (𝐾)
4.3.1 When two finite bodies connected to a heat engine at temperature T1 and T2 with T1>T2
5 Pure substance
5.1 General property relation in the wet region
5.2.3 Initial dryness fraction between states 1 and 2 using Separating and throttling calorimeter
𝐶𝑝
𝛾=
𝐶𝑣
∆ℎ = ℎ2 − ℎ1 = 𝑐𝑝 (𝑇2 − 𝑇1 ) 𝐽/𝑘𝑔
6.1.5 Entropy change for any ideal gas and for any process
(𝛾 − 𝑛) 𝑛 = 𝑝𝑜𝑙𝑦𝑡𝑟𝑜𝑝𝑖𝑐 𝑖𝑛𝑑𝑒𝑥
𝑞𝑅 = [ 𝑅(𝑇1 − 𝑇2 ) ] 𝐽/𝑘𝑔
(𝛾 − 1) (𝑛 − 1)
𝑛𝑖
=
𝑛𝑚
6.3.2 Molecular weight of the mixture (Mm)
1.1.1 𝑖=𝑘
𝑀𝑚 = ∑ 𝑦𝑖 ∗ 𝑀𝑖
𝑖=1
1.1.3 𝑚1 𝑢1 + 𝑚2 𝑢2 + 𝑚3 𝑢3
𝑢𝑚 = [ ] 𝐽/𝑘𝑔
𝑚1 + 𝑚2 + 𝑚3
𝑚1 ℎ1 + 𝑚2 ℎ2 + 𝑚3 ℎ3
ℎ𝑚 = [ ] 𝐽/𝑘𝑔
𝑚1 + 𝑚2 + 𝑚3
Probability
Addition rule
The conditional probability of an event B, given that the event A already taken place is
P( A B)
P( B / A) , P( A) 0.
P( A)
Baye’s Theorem
Let B1 , B2 ,...Bk are the partitions of S with P( Bi ) 0, i 1, 2,...k and A be any event of S, then
P( A / Bi ) P( Bi ) .
P( Bi / A) k
P( A / B ) P( B )
i 1
i i
𝑃(𝐴)𝑃(𝐵|𝐴), 𝑖𝑓 𝑃(𝐴) ≠ 0
The multiplicative rule of probability : 𝑃(𝐴 ∩ 𝐵) = {
𝑃(𝐵)𝑃(𝐴|𝐵), 𝑖𝑓 𝑃(𝐵) ≠ 0
If 𝑃(𝐴 ∩ 𝐵) = 𝑃(𝐴)𝑃(𝐵), then A and B are independent.
Random Variable: Let S be the sample of space of a random experiment. Suppose with each
element s of S, a unique real number X is associated according to some rule then X is called random
variable. There are two types of random variable, i) Discrete and ii) Continuous.
Discrete Random Variable: A random variable X is said to be discrete, if the number of possible
values of X is finite or countably infinite. The probability distribution function (pdf) is named as
probability mass function (PMF). The Probability mass function is defined as, let X be a random
variable, hence the range space 𝑅𝑋 of consists of atmost a countably infinite number of values. The
probability mass function is defined as
ii) ∑𝑘
𝑖=1 𝑝(𝑥𝑖 ) = 1.
Continuous Random Variable: A random variable X is said to be continuous if it can take all
possible values between certain limits, here the range space of X is infinite. Therefore the probability
distribution function named for such random variable is Probability density function (PDF), which is
defined as the pdf of X is a function 𝑓(𝑥) satisfying the following properties i) 𝑓(𝑥) ≥ 0
∞
ii) ∫
−∞
𝑓(𝑥)𝑑𝑥 =1
𝑏
iii) Pr{𝑎 ≤ 𝑋 ≤ 𝑏} = ∫𝑎 𝑓(𝑥)𝑑𝑥 for any a, b such that −∞ < 𝑎 < 𝑏 < ∞.
Note:
Discrete 2D: If the possible values of (X, Y) are finite or countably infinite then (X, Y) is called discrete
and it is defined as 𝑃(𝑥𝑖 , 𝑦𝑗 ) satisfying the following condition,
i) 𝑃(𝑥𝑖 , 𝑦𝑗 ) ≥ 0 and
ii) ∑∞ ∞
𝑗=1 ∑𝑖=1 𝑃(𝑥𝑖 , 𝑦𝑗 ) = 1. The function 𝑃(𝑥𝑖 , 𝑦𝑗 ) defined is called as Joint probability
Continuous 2D: If (X, Y) is a continuous random variable assuming all values in some region R of
the Euclidean plane, then the Joint probability density function 𝑓(𝑥, 𝑦) is a function satisfying the
following conditions
Case ii) In the continuous (X, Y), it is defined as the marginal probability function of X is defined as
∞ ∞
𝑔(𝑥) = ∫−∞ 𝑓(𝑥, 𝑦)𝑑𝑦 and the marginal probability function of Y is defined as ℎ(𝑦) = ∫−∞ 𝑓(𝑥, 𝑦)𝑑𝑥.
To calculate the conditional probability:
𝑃(𝑥𝑖 ,𝑦𝑗 )
Case i) Discrete: Probability of 𝑥𝑖 given 𝑦𝑗 is defined as = , 𝑞(𝑦𝑗 ) > 0
𝑞(𝑦𝑗 )
𝑃(𝑥𝑖 ,𝑦𝑗 )
Probability of 𝑦𝑗 given 𝑥𝑖 is defined as = , 𝑝(𝑥𝑖 ) > 0
𝑝(𝑥𝑖 )
𝑓(𝑥,𝑦)
Case ii) Continuous: The pdf of X for given Y=y is = , ℎ(𝑦) > 0
ℎ(𝑦)
𝑓(𝑥,𝑦)
The pdf off Y for given X=x is = , 𝑔(𝑥) > 0.
𝑔(𝑥)
Independent Random variable: If X and Y are independent random variable then two dimensional
random variable in case of discrete is defined as 𝑃(𝑥𝑖 , 𝑦𝑗 ) = 𝑝(𝑥𝑖 ). 𝑞(𝑦𝑗 ) for all the values of i and
j. In case of Continuous it is defined as 𝑓(𝑥, 𝑦) = 𝑔(𝑥). ℎ(𝑦).
Mathematical Expectation: If X is a discrete random variable with pmf p(x), then the expectation of X
is given by 𝐸(𝑋) = ∑𝑥 𝑥𝑝(𝑥), provided the series is absolutely convergent.
DISTRIBUTIONS
Distribution PMF/PDF Mean Variance
Binomial P( x) Ck p (1 p)
n k n k
, k 0,1, 2,..., n E ( x) np V ( x) np(1 p)
distribution
X ~ B(n, p)
Poisson’s e k E ( x) np V ( x) np
Distribution P( x) , k 0,1, 2,..., 0
X ~ P( ) k!
Uniform 1 ba (b a)2
, a xb E ( x) V ( x)
f ( x) b a
Distribution
2 12
X ~ U (a, b)
0, otherwise
Normal
1
1 ( x )2 E ( x) V ( x) 2
Distribution f ( x) e2 2
, x , , 0
X ~ N ( , 2 ) 2
1
Exponential e x , x0 E ( x)
1 𝑉(𝑋) =
Distribution f ( x) 𝜆2
X ~ E ( ) 0, otherwise
Gamma x r 1e x r E ( x)
r
V ( x)
r
, x 0, , r 0
f ( x) (r )
Distribution
X ~ G (r , ) 2
0, elsewhere
Chi-square n
1
x E ( x) n V ( x) 2n
x e
Distribution 2 2
, x0
X ~ 2 ( n) f ( x) n
(n / 2)2 2
0,
elsewhere
Uniform distribution on a two dimensional set: If R is a set in the two-dimensional plane, and R has a
finite area, then we may consider the density function equal to the reciprocal of the area of R inside R,
and equal to 0 otherwise:
1
; 𝑖𝑓 (𝑥, 𝑦) ∈ 𝑅
𝑓(𝑥, 𝑦) = {𝑎𝑟𝑒𝑎 𝑅 .
0 𝑂𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒
Chebyshev’s inequality:
Let x be random variable with mean and variance 2 then for any positive real number k(k>0)
2
p X k (Upper bound)
k2
2
p X k 1 (Lower bound)
k2
p X k p x k 1 (Upper bound)
1
1. and
k2
p X E ( x c) 2 and p x 1
1
2.
2
Covariance:
Cov( X , Y ) E ( XY ) E ( X ) E (Y )
Correlation coefficient:
E ( XY ) E ( X ) E (Y )
xy
V ( X )V (Y )
Properties:
1. 𝐸(𝑐) = 𝑐, where c is a constant.
2. 𝑉(𝑐) = 0, where c is a constant.
3. If E ( XY ) 0 then X and Y are orthogonal.
6. V ( AX BY ) A2V ( X ) B 2V (Y ) 2 ABCOV ( X , Y )
Regression:
If the straight line chosen such that the sum of the squares of errors parallel to y axis is minimum then
it is called line of regression of y on x and it gives the best estimate of y for any given value of x. It’s
given by
y y
y y ( x x) , b yx is the regression coefficient of Y on X
x x
Therefore y y byx ( x x)
x x
xx ( y y) , bxy is the regression coefficient of X on Y
y y
Therefore x x bxy ( y y)
Note:
x x
2 2
2
E ( x 2 ) ( E ( x)) 2 is the variance of X
x
n n
y y
2 2
2
E ( y 2 ) ( E ( y )) 2 is the variance of Y
y
n n
E ( xy ) E ( x) E ( y )
V ( x)V ( y )
Also, bxy b yx . The correlation coefficient is positive if the two regression coefficients are
positive, and negative if the two regression coefficients are negative.
Let 𝑋 be a random variable defined on 𝑆. Let 𝑦 = 𝐻(𝑥) is a real valued function of 𝑥. Then 𝑌 = 𝐻(𝑋) is
a random variable on 𝑆. i.e., for each element s S, there is a real number associated, say 𝑦 = 𝐻(𝑋(𝑠)).
Here 𝑌 is called a function of the random variable 𝑋.
Notations:
1. 𝑅𝑋 – the set of all possible values of the function 𝑋, called the range space of the random
variable 𝑋.
2. 𝑅𝑌 – the set of all possible values of the function 𝑌 = 𝐻(𝑋), called the range space of the
random variable 𝑌.
Equivalent Events: Let C be an event associated with the range space 𝑅𝑌 . Let 𝐵 ⊂ 𝑅𝑋 defined by 𝐵 =
{𝑥 ∈ 𝑅𝑋 ; 𝐻(𝑥) ∈ 𝐶}, then 𝐵 and 𝐶 are called equivalent events.
Case 1: Let 𝑋 be a discrete random variable with p.m.f. 𝑝(𝑥𝑖 ) = 𝑃(𝑋 = 𝑥𝑖 ) for 𝑖 = 1,2,3, … Let 𝑌 =
𝐻(𝑋) then 𝑌 is also a discrete random variable. If 𝑌 = 𝐻(𝑋) is a one to one function then the probability
distribution of 𝑌 is as follows:
For the possible values of 𝑦𝑖 = 𝐻(𝑥𝑖 ) for 𝑖 = 1,2,3, …. The p.m.f. of 𝑌 = 𝐻(𝑋) is 𝑞(𝑦𝑖 ) = 𝑃(𝑌 = 𝑦𝑖 ) =
𝑃(𝑋 = 𝑥𝑖 ) = 𝑝(𝑥𝑖 ) for 𝑖 = 1,2,3, ….
Case 2: Let 𝑋 be a discrete random variable with p.m.f. 𝑝(𝑥𝑖 ) = 𝑃(𝑋 = 𝑥𝑖 ) for 𝑖 = 1,2,3, … Let 𝑌 = 𝐻(𝑋)
then 𝑌 is also a discrete random variable. Suppose tat for one value of 𝑌 = 𝑦𝑖 there corresponds several
values of 𝑋 say 𝑥𝑖1 , 𝑥𝑖2 , … , 𝑥𝑖𝑗 , … then the p.m.f. of of 𝑌 = 𝐻(𝑋) is
Case 3: Let 𝑋 be a continuous random variable with p.d.f. 𝑓(𝑥). Let 𝑌 = 𝐻(𝑋) be a discrete random
variable. Then if the set {𝑌 = 𝑦𝑖 } is equivalent to an event 𝐵 ⊆ 𝑅𝑋 then the p.m.f. of 𝑌 is
Case 4: Let 𝑋 be a continuous random variable with p.d.f. 𝑓(𝑥). Let 𝑌 = 𝐻(𝑋) be a continuous random
variable. Then the p.d.f. of 𝑌, say 𝑔 is obtained by the following procedure:
Theorem: Let 𝑋 be a continuous random variable with p.d.f. 𝑓(𝑥) where 𝑓(𝑥) > 0 for 𝑎 < 𝑥 < 𝑏.
Suppose that 𝑌 = 𝐻(𝑋) is strictly monotonic function on [𝑎, 𝑏]. Then the p.d.f. of the random variable
𝑌 = 𝐻(𝑋) is given by
𝑑𝑥
𝑔(𝑦) = 𝑓(𝑥) | |
𝑑𝑦
If 𝑌 = 𝐻(𝑋) is strictly increasing then 𝑔(𝑦) > 0 for 𝐻(𝑎) < 𝑦 < 𝐻(𝑏).
If 𝑌 = 𝐻(𝑋) is strictly decreasing then 𝑔(𝑦) > 0 for 𝐻(𝑏) < 𝑦 < 𝐻(𝑎).
Theorem: Let 𝑋 be a continuous random variable with p.d.f. 𝑓(𝑥). Let 𝑌 = 𝑋 2 then the p.d.f. of 𝑌 is
1
𝑔(𝑦) = [𝑓(√𝑦) + 𝑓(−√𝑦) ]
2√𝑦
𝑊 = 𝐺(𝑋, 𝑌) and obtain the joint p.d.f. of the two dimensional random variable (𝑍, 𝑊), say 𝑘(𝑧, 𝑤). From
this distribution, the p.d.f. of 𝑍 can be obtained by integrating 𝑘 with respect to 𝑤.
Theorem: Suppose (𝑋, 𝑌) is a two dimensional continuous random variable with joint p.d.f. 𝑓(𝑥, 𝑦)
defined on a region 𝑅 of the XY-plane. Let 𝑍 = 𝐻1 (𝑋, 𝑌) and 𝑊 = 𝐻2 (𝑋, 𝑌). Suppose that 𝐻1 and 𝐻2
satisfies the following conditions;
(i) 𝑧 = 𝐻1 (𝑥, 𝑦) and 𝑤 = 𝐻2 (𝑥, 𝑦) may be uniquely solved for 𝑥, 𝑦 in terms of 𝑧 & 𝑤 say,
𝑥 = 𝐺1 (𝑧, 𝑤) and 𝑦 = 𝐺2 (𝑧, 𝑤).
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
(ii) The partial derivatives , , 𝜕𝑧 , 𝜕𝑤 exist and are continuous
𝜕𝑧 𝜕𝑤
(𝑥, 𝑦) ↦ (𝑧, 𝑤). Also, 𝑘(𝑧, 𝑤) > 0 for those values of (𝑧, 𝑤) corresponding to the values of (𝑥, 𝑦) for
which 𝑓(𝑥, 𝑦) > 0.
Properties of m.g.f.: Let 𝑋 be any one dimensional random variable and 𝑀𝑋 (𝑡) be the m.g.f. of 𝑋 then
1. 𝑀𝑋𝑛 (0) = 𝐸(𝑋 𝑛 ) where 𝑀𝑋𝑛 (0) is the nth derivative of 𝑀𝑋 (𝑡) at 𝑡 = 0.
i.e.; 𝑀𝑋′ (0) = 𝐸(𝑋)
𝑀𝑋′′ (0) = 𝐸(𝑋 2 )
2
2. 𝑉(𝑋) = 𝑀𝑋′′ (0) − (𝑀𝑋′ (0))
3. Let 𝑋 be any one dimensional random variable and 𝑀𝑋 (𝑡) be the m.g.f. of 𝑋. Let 𝑌 = 𝛼𝑋 + 𝛽.
Then the m.g.f. of 𝑌 is
𝑀𝑌 (𝑡) = 𝑒 𝛽𝑡 𝑀𝑋 (𝛼𝑡).
4. Suppose that 𝑋 and 𝑌 are independent random variables. Let 𝑍 = 𝑋 + 𝑌. Let 𝑀𝑋 (𝑡), 𝑀𝑌 (𝑡)
and 𝑀𝑍 (𝑡) be the m.g.f.’s of the random variables 𝑋, 𝑌 and 𝑍 respectively. Then
𝑀𝑍 (𝑡) = 𝑀𝑋 (𝑡)𝑀𝑌 (𝑡)
5. Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be 𝑛 independent random variables which follows a normal distribution
𝑁(𝜇𝑖 , 𝜎𝑖2 ) for 𝑖 = 1,2,3, . . , 𝑛. Let 𝑍 = 𝑋1 + 𝑋2 + ⋯ + 𝑋𝑛 then
𝑍 → 𝑁(𝜇1 + 𝜇2 + ⋯ + 𝜇𝑛 , 𝜎12 + 𝜎22 + ⋯ + 𝜎𝑛2 ).
6. Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be 𝑛 independent random variables which follows a Poisson distribution with
parameter 𝛼𝑖 for 𝑖 = 1,2,3, . . , 𝑛. Let 𝑍 = 𝑋1 + 𝑋2 + ⋯ + 𝑋𝑛 then 𝑍 has a Poisson distribution with
parameter
𝛼 = 𝛼1 + 𝛼2 + ⋯ + 𝛼𝑛 .
7. Let 𝑋1 , 𝑋2 , … , 𝑋𝑘 be 𝑘 independent random variables which follows a Chi-square distribution with
degrees of freedom 𝑛𝑖 for 𝑖 = 1,2,3, . . , 𝑘. Let 𝑍 = 𝑋1 + 𝑋2 + ⋯ + 𝑋𝑘 then 𝑍 has a Chi-square
distribution with degrees of freedom
𝑛 = 𝑛1 + 𝑛2 + ⋯ + 𝑛𝑘 .
8. Let 𝑋1 , 𝑋2 , … , 𝑋𝑘 be 𝑘 independent random variables, each having distribution 𝑁(0,1). Then 𝑆 =
𝑋12 + 𝑋22 + ⋯ + 𝑋𝑘2 has a Chi-square distribution with degrees of freedom 𝑘.
9. Let 𝑋1 , 𝑋2 , … , 𝑋𝑟 be 𝑟 independent random variables, each having exponential distribution with
the same parameter 𝛼. Let 𝑍 = 𝑋1 + 𝑋2 + ⋯ + 𝑋𝑟 then 𝑍 has a Gamma distribution with
parameters 𝛼 and 𝑟.
∝
4. Exponential Distributions: 𝑀𝑋 (𝑡) = ∝−𝑡
∝𝑟
5. Gamma Distributions: 𝑀𝑋 (𝑡) = (∝−𝑡)𝑟
Sampling Theory
In statistical investigation, the characteristics of a large group of individuals (called population) is
studied. Sampling is a study of the relationship between a population and samples drawn from it.
The population mean and the population variance are denoted by 𝜇 and 𝜎 2 respectively.
Sample mean and sample variance: Let 𝑋 be the random variable which denotes the population
with mean 𝜇 and variance 𝜎 2 . Let (𝑋1 , 𝑋2 , … , 𝑋𝑛 ) be a random sample of size 𝑛 from 𝑋. Then,
∑𝑛
𝑖=1 𝑋𝑖
̅
Sample mean, 𝑋 = and
𝑛
∑𝑛 ̅ 2
𝑖=1(𝑋𝑖 −𝑋)
Sample variance, 𝑠 2 = 𝑛
Central Limit Theorem: Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be 𝑛 independent random variables all of which have the
same distribution. Let 𝜇 = 𝐸(𝑋𝑖 ) and 𝜎 2 = 𝑉(𝑋𝑖 ) be the common expectation and variance. Let 𝑆 =
𝑆−𝐸(𝑆)
∑𝑛𝑖=1 𝑋𝑖 then 𝐸(𝑆) = 𝑛𝜇 and 𝑉(𝑆) = 𝑛𝜎 2 then for large values of 𝑛, the random variable 𝑇𝑛 = has
√𝑉(𝑆)
approximately the distribution 𝑁(0,1).
Curve fitting
Let y = f(x) be the curve to be fit with n set of points and Y = F(x) is actual values, then the estimated
error is 𝑒𝑖 = 𝑦𝑖 − 𝑌𝑖 .
i) Fitting a straight line: For a straight line 𝑦 = 𝑎𝑥 + 𝑏, where a, b are the parameters, the normal
equations are to find a and b is given by
ii) Fitting a Parabola: For a parabola 𝑦 = 𝑎𝑥 2 + 𝑏𝑥 + 𝑐, the normal equations are given by
∑ 𝑌 = 𝑎 ∑ 𝑋 2 + 𝑏 ∑ 𝑋 + 𝑛𝑐
∑ 𝑋𝑌 = 𝑎 ∑ 𝑋 3 + 𝑏 ∑ 𝑋 2 + 𝑐 ∑ 𝑋
∑ 𝑋 2𝑌 = 𝑎 ∑ 𝑋 4 + 𝑏 ∑ 𝑋 3 + 𝑐 ∑ 𝑋 2
iii) Fitting a curve of the form 𝑦 = 𝑎𝑒 𝑏𝑥 , the normal equations are given by taking log on both
sides, 𝑙𝑜𝑔𝑒 𝑦 = 𝑙𝑜𝑔𝑒 (𝑎𝑒 𝑏𝑥 )
𝑙𝑜𝑔𝑒 𝑦 = 𝑙𝑜𝑔𝑒 𝑎 + 𝑏𝑥, Now consider 𝑌 = 𝑙𝑜𝑔𝑒 𝑦, 𝐴 = 𝑙𝑜𝑔𝑒 𝑎, 𝑏 = 𝐵, 𝑋 = 𝑥.
Therefore equation is 𝑌 = 𝐴𝑋 + 𝐵.
Special functions
𝑥 1
1. Bessel function 𝐽𝑛 (𝑥) = ∑∞ 𝑟
𝑟=0(−1) (2)
𝑛+2𝑟
𝑟! 𝛾(𝑛+𝑟+1)
0,
1
d2y dy
5. The differential equation of the form (1 x ) 2 2 x n(n 1) y 0
2
is the
dx dx
Legendre’s equation whose general solution when n is a positive integer, is polynomial solution
and an infinite series solution. The polynomial solution is called Legendre’s polynomial of order
n and infinite series solution is called Legendre’s function of the second kind.
1. Static Forces
1.1 Static Equilibrium
2. Inertia Forces
2.1 Gas Force, Fg (N)
𝑠𝑖𝑛2𝜃
T = 𝑃𝑟 [𝑠𝑖𝑛𝜃 + ] P = Net piston effort in N
2√𝑛2 − 𝑠𝑖𝑛2 𝜃
r = Radius of the crank (m)
θ = Angular position of the crank from horizontal
reference (degrees)
n = Ratio of the length of connecting rod to crank
radius
E = Tmean * θ
Tmean = Mean Torque (Nm)
θ = Angular position of the crank (radians)
OR
𝜔1 −𝜔2
𝜔1 = 𝜔𝑚𝑎𝑥 = Max. angular speed (rad/s)
𝐾𝑠 = 𝜔2 = 𝜔𝑚𝑖𝑛 = Min. angular speed (rad/s)
𝜔0
𝜔0 = 𝜔𝑚𝑒𝑎𝑛 = Mean angular speed (rad/s)
1
𝐾𝐸𝑚𝑒𝑎𝑛 = 2
𝐼𝜔𝑚𝑒𝑎𝑛 I = Moment of Inertia (N-m-s2)
2
𝜔𝑚𝑒𝑎𝑛 = Mean angular speed (rad/s)
OR
𝐾𝑒 𝐸
𝐾𝐸𝑚𝑒𝑎𝑛 = 2 𝐾𝑠
E = Work done per cycle (Nm)
K s = Fluctuation of speed (Unit less)
Ke = Coefficient of fluctuation of energy (Unit less)
PUNCHING PRESS:
Wp = 𝜋 ∗ 𝑑 ∗ 𝑡 ∗ 𝑋
d = Diameter of the punch (cm)
1
𝑊𝑝 = ( 𝜋 ∗ 𝑑 ∗ 𝑡 ∗ 𝜏𝑢 ) ∗ 𝑡 d = Diameter of the punch (cm)
2
t = Thickness of the plate (cm)
𝜏𝑢 = Ultimate shear strength of the plate (N/cm 2)
OR
𝐸𝑓 = 𝑊𝑝 [1 − { 2
𝜃 − 𝜃1
}] 𝑊𝑝 = Work required per punching (Nm)
2𝜋 𝜃1 = Initial angular position of the crank (radians)
𝜃2 = Final angular position of the crank (radians)
OR
3. Governors
3.1 Height of the Watt Governor, ℎ (m)
895
ℎ= h = Height of the governor (m)
𝑁2
N = Speed of fly ball in revolutions per minute
895 2𝑚𝑔+(𝑀𝑔±𝑓)(1+𝑘)
𝑁2 = [ ]
ℎ 2𝑚𝑔
h = Height of the governor (m)
M = Mass of the sleeve (kg)
m = Mass of each ball (kg)
f = Force of friction at the sleeve (N)
g = Acceleration due to gravity (m/s2)
k = Ratio of tan β / tan θ
β = Angle of inclination of the lower link to vertical
(degrees)
θ = Angle of inclination of the upper link to vertical
(degrees)
895 𝑎 2𝑚𝑔+(𝑀𝑔±𝑓)(1+𝑘)
𝑁2 = [ ]
ℎ 𝑒 2𝑚𝑔
𝐹𝑠2 −𝐹𝑠1
𝑠= hL = lift of the sleeve (mm)
ℎ𝐿
Fs2 = spring force on sleeve at maximum radius (N)
OR Fs1= spring force on sleeve at minimum radius (N)
𝑎2 𝐹 −𝐹
𝑠 = 2 [( 2 ) ∗ ( 2 1)]
𝑏 𝑟 −𝑟 2 1
𝐸
3.5 Effort of a Porter Governor, (N)
2
𝐸 𝑐𝑔
2
= (1+𝑘) [2𝑚 + 𝑀(1 + 𝑘)]
c = Fractional increase of speed (Unit less)
m = Mass of each fly ball (kg)
M = Mass of the sleeve (kg)
g = Acceleration due to gravity (m/s2)
k = Ratio of tan β / tan θ
β = Angle of inclination of the lower link to vertical
(degrees)
θ = Angle of inclination of the upper link to vertical
(degrees)
𝐸
3.6 Effort of a Hartnell Governor, (N)
2
𝐸
= 𝑐(𝑀 ∗ 𝑔 + 𝐹𝑠 ) c = Fractional increase of speed (Unit less)
2
M = Mass of the sleeve (kg)
Fs = Spring force exerted on sleeve (N)
g = Acceleration due to gravity (m/s2)
𝑀 4𝑐2
If 𝑘 ≠ 1; 𝑃 = [𝑚 + 2
(1 + 𝑘)]𝑔ℎ [1 +2𝑐
]
4. Gyroscope
5. Balancing
FOR ROTATING MASSES:
5.1 Balancing of a single rotating mass by a single mass rotating in the same plane.
2 + 𝐹2
𝐹𝑝𝑣 = Vertical component of primary disturbing
𝐹𝑝𝑜 = √𝐹𝑝𝑣 𝑝ℎ force (N)
𝐹𝑝ℎ = Horizontal component of primary disturbing
force (N)
5.14 Primary force due to each direct (or reverse) crank, 𝐹𝑝−𝑑𝑖𝑟𝑒𝑐𝑡 ( 𝑜𝑟 𝐹𝑝−𝑟𝑒𝑣𝑒𝑟𝑠𝑒 ) (N)
𝐹𝑝−𝑚𝑎𝑥 = 𝐹𝑝−𝑑𝑖𝑟𝑒𝑐𝑡 + 𝐹𝑝−𝑟𝑒𝑣𝑒𝑟𝑠𝑒 𝐹𝑝−𝑑𝑖𝑟𝑒𝑐𝑡 = Primary force due to each direct crank (N)
𝐹𝑝−𝑟𝑒𝑣𝑒𝑟𝑠𝑒 = Primary force due to each reverse crank (N)
Also 𝐹𝑠−𝑑𝑖𝑟𝑒𝑐𝑡 = Secondary force due to each direct crank (N)
𝐹𝑠−𝑟𝑒𝑣𝑒𝑟𝑠𝑒 = Secondary force due to each reverse crank (N)
𝐹𝑠−𝑚𝑎𝑥 = 𝐹𝑠−𝑑𝑖𝑟𝑒𝑐𝑡 + 𝐹𝑠−𝑟𝑒𝑣𝑒𝑟𝑠𝑒
𝐹𝑝−𝑚𝑖𝑛 = 𝐹𝑝−𝑑𝑖𝑟𝑒𝑐𝑡 − 𝐹𝑝−𝑟𝑒𝑣𝑒𝑟𝑠𝑒 𝐹𝑝−𝑑𝑖𝑟𝑒𝑐𝑡 = Primary force due to each direct crank (N)
𝐹𝑝−𝑟𝑒𝑣𝑒𝑟𝑠𝑒 = Primary force due to each reverse crank (N)
Also 𝐹𝑠−𝑑𝑖𝑟𝑒𝑐𝑡 = Secondary force due to each direct crank (N)
𝐹𝑠−𝑟𝑒𝑣𝑒𝑟𝑠𝑒 = Secondary force due to each reverse crank (N)
𝐹𝑠−𝑚𝑖𝑛 = 𝐹𝑠−𝑑𝑖𝑟𝑒𝑐𝑡 − 𝐹𝑠−𝑟𝑒𝑣𝑒𝑟𝑠𝑒
1. Properties of Fluids
1.1 Mass density or Specific mass (kg/m 3)
𝑀 M = Mass (kg)
𝜌= 𝑉 V = Volume (m3)
For water at 4oC, ρ = 1000 kg/m 3
For air at 15oC at sea level, ρ = 1.22 kg/m3
1.2 3
Specific weight or Weight density (N/m )
𝑊 W = Weight (N)
𝛾= = 𝜌𝑔
𝑉 𝑉 = Volume (m3)
= Density (kg/m3)
g = Acceleration due to gravity = 9.81 m/s2
1.3 Specific gravity or Relative density
𝜌𝑓𝑙𝑢𝑖𝑑 ρstandard fluid = ρwater = 1000 kg/m3
𝑆=𝜌
𝑠𝑡𝑎𝑛𝑑𝑎𝑟𝑑 𝑓𝑙𝑢𝑖𝑑 S = 1 (for water at 4oC)
S = 13.6 (for mercury)
1.4 Specific volume (m 3/ kg)
𝑉 M = Mass (kg)
𝑣=𝑀
𝑉 = Volume (m3)
5. Fluid Kinematics
5.1 Velocity field
⃗ × 𝑑𝑠 = 0
𝑉 ⃗ = Velocity vector (m/s)
𝑉
For 2D flow: 𝑑𝑠 = vector element along stream line (m)
𝑑𝑦 𝑣 𝑢, 𝑣 = Scalar components of velocity along x and y
= directions respectively (m/s)
𝑑𝑥 𝑢
𝜕𝑢 𝜖𝑥̇ , 𝜖𝑦̇ and 𝜖𝑧̇ = Linear strain rates of fluid element along
𝜖𝑥̇ =
𝜕𝑥 𝑥, 𝑦 and 𝑧 directions respectively
𝜕𝑣 𝑢, 𝑣 and 𝑤 = Scalar components of velocity in x, y and
𝜖𝑦̇ =
𝜕𝑦 z directions respectively (m/s)
𝜕𝑤
𝜖𝑧̇ =
𝜕𝑧
𝜕 𝜕
𝑢 = − 𝜕𝑦 and 𝑣 = 𝜕𝑥
𝜕𝜙 𝜕𝜙 𝜕𝜙
𝑢= 𝜕𝑥
, 𝑣 = 𝜕𝑦 , 𝑤 = 𝜕𝑧
6. Fluid Dynamics
6.1 Euler’s equation of motion
𝐼𝑛𝑒𝑟𝑡𝑖𝑎 𝑓𝑜𝑟𝑐𝑒 𝑉
𝐹𝑒 = √ =
𝐺𝑟𝑎𝑣𝑖𝑡𝑦 𝑓𝑜𝑟𝑐𝑒 √𝐿𝑔
𝐼𝑛𝑒𝑟𝑡𝑖𝑎 𝑓𝑜𝑟𝑐𝑒 𝑉
𝐸𝑢 = √ =
𝑃𝑟𝑒𝑠𝑠𝑟𝑒 𝑓𝑜𝑟𝑐𝑒
√𝑝⁄𝜌
𝐼𝑛𝑒𝑟𝑡𝑖𝑎 𝑓𝑜𝑟𝑐𝑒 𝑉
𝑊𝑒 = √ =
𝑆𝑢𝑟𝑓𝑎𝑐𝑒 𝑡𝑒𝑛𝑠𝑖𝑜𝑛 𝑓𝑜𝑟𝑐𝑒
√𝜎⁄𝜌𝐿
8.5 Mach number
𝐼𝑛𝑒𝑟𝑡𝑖𝑎 𝑓𝑜𝑟𝑐𝑒 𝑉
𝑀=√ =
𝐸𝑙𝑎𝑠𝑡𝑖𝑐 𝑓𝑜𝑟𝑐𝑒 𝐶
9. Viscous Flow
9.1 Flow of viscous fluid through a circular pipe:
𝜋 𝜕𝑝 p
= Pressure gradient (N/m 3)
𝑄= (− ) 𝑅 4 x
8𝜇 𝜕𝑥
= Dynamic viscosity (Ns/m2)
9.1.4 Drop of pressure over a given length of pipe hf = Drop of pressure head or head lost
(Hagen-Poiseuille formula) due to friction (m)
ρ = Density of fluid (kg/m 3)
32uL 128QL g = Acceleration due to gravity = 9.81
p p2 p1 m/s2
D2 D4
u = Velocity at radial distance r (m/s)
p 32 uL 128QL
hf u = Average velocity (m/s)
g gD 2 g D 4 FD = Total drag force (N)
2 x (0 ≤ y ≤ t)
p
9.2.3 Volume flow rate per unit width = Pressure gradient (N/m 3)
x
1 𝜕𝑝 3 u = Velocity at distance y (m/s)
𝑄=− ( )𝑡
12𝜇 𝜕𝑥 u = Average velocity (m/s)
9.2.4 Drop in pressure for a length L of plates (Couette p = Pressure drop (N/m 2)
flow) p1 and p2 = Pressure intensities at sections
1 and 2
12uL
p p1 p2
L = Length of the plates (m)
t2 = Dynamic viscosity (Ns/m 2)
hf = Drop of pressure head or head lost
p 12uL 12 LQ
hf due to friction (m)
g gt 2 gbt 3 b = Width of the plates (m)
ρ = Density of fluid (kg/m 3)
9.2.5 The total drag force on the plates of length L per g = Acceleration due to gravity = 9.81
unit width m/s2
FD max 2L
9.3 Power required to maintain a laminar flow between two sections 1 and 2
P p2 p1 Q P = Power (W)
p1 and p2 = Pressure intensities at sections
1 and 2 (N/m2)
Q = Volume flow rate (m 3/s)
10.2.1 Fanning friction factor or Skin friction Cf = Fanning friction factor or Skin friction
coefficient coefficient
1⁄ (𝐷ℎ ) ∆𝑝∗ hf = Head lost due to friction (m)
𝜏𝑤 4 𝐿
𝐶𝑓 = = w = Wall shear stress resistance to flow (N/m 2)
1 2 1 2
2 𝜌𝑉 2 𝜌𝑉 Dh =
4𝐴
= Hydraulic diameter (m)
1 ℎ𝑓 𝑆
=( ) where, A = Area of cross section of flow (m 2)
4 𝐿 𝑉2
( )( ) S = Wetted perimeter (m)
𝐷ℎ 2𝑔
For circular pipe, Dh = D = pipe diameter (m)
MME 2252: Fluid Mechanics 77
10.2.2 Head lost due to friction = Density of fluid (kg/m3)
V = Average velocity of flow (m/s)
∆𝑝∗ 4𝐶𝑓 𝐿𝑉 2 L = Length of pipe (m)
ℎ𝑓 = =
𝜌𝑔 2𝑔𝐷ℎ p* = Piezo metric pressure drop over length L
(N/m2)
10.2.3 Darcy-Weisbach equation (if pipe is horizontal, p* = static pressure drop,
p)
𝑓𝐿𝑉 2 g = Acceleration due to gravity = 9.81 m/s2
ℎ𝑓 =
2𝑔𝐷ℎ f = Friction factor
ℎ𝑓
𝑓 = 4𝐶𝑓 =
𝐿 𝑉2
( )( )
𝐷ℎ 2𝑔
𝑓 = 𝜑(𝑅𝑒, 𝜀)
10.4 Fig. 10.1 gives the Moody’s chart which is a plot of curves between friction factor f and relative
roughness / D. It is used for finding the friction factor for circular pipes with smooth and rough
walls.
10.4.1 Head lost due to sudden enlargement V = Average velocity of flow (m/s)
(V1 V2 ) 2 V1 and V2 = Velocity of flow immediately before
he and after the enlargement/contraction of pipe
2g (m/s)
Vo = Exit velocity (m/s)
10.4.2 Head lost due to sudden contraction
kb = Coefficient of bend
kV2 2
ho kob = Coefficient of obstruction
2g kfit = Coefficient of pipe fittings
1
2 g = Acceleration due to gravity = 9.81 m/s2
where, k 1
Cc
Area at vena contracta
Cc
Area after contraction
If Cc is not known, take k = 0.5
c
C ( A a ) ( A a)
10.4.7 Head lost due to an obstruction in the A = Cross sectional area of pipe (m ) 2
pipe a = Maximum area of obstruction (m2)
kobV 2 ac = Cross sectional area at vena contracta (m2)
hob
2g
Fig 10.1 The Moody chart for pipe friction with smooth and rough walls
(Source: Bruce R. Munson, Donald F. Young and Teodore H. Okiishi, Fundamentals of Fluid Mechanics, Wiley)
Energy thickness
u
u 2
** 1 dy
0U U
11.4 Von Karman momentum integral equation
0.455
CDL
log10 ReL
2.58
FD p cos dA o sin dA FD = Total drag force acting on the solid body (N)
FL = Total lift force acting on the solid body (N)
U 2 FR = Resultant force exerted by fluid on body (N)
FD CD A
2 p = Pressure exerted by fluid on the solid body, normal
to its surface (N/m2)
o = Shear stress exerted by fluid on the solid body,
Total lift force tangential to its surface (N/m 2)
FL p sin dA o cos dA = Inclination of pressure force on a differential area
dA of solid body surface (degree)
U 2
FL CL A CD = Coefficient of drag
2
CL = Coefficient of lift
= Density of the fluid (kg/m 3)
Resultant force U = Free stream velocity of fluid (m/s)
A = Projected area of the solid body perpendicular to
FR FD2 FL2
the fluid flow OR largest projected area of immersed
body (m2)
1.2 Sensitivity:
Change of output signal
Static Sensitivity, K =
𝐶ℎ𝑎𝑛𝑔𝑒 𝑖𝑛 𝑖𝑛𝑝𝑢𝑡 𝑠𝑖𝑔𝑛𝑎𝑙
1.3 Span:
1.4 Linearity:
1.5 Error:
2. Pressure Measurement
2.1 The following are the units and conversion factors that are normally used:
1 Pa = 1 N/m2
1 atmosphere = 760 mm of Hg
1 mm of Hg = 1 Torr
1 Torr = 1.316 x 10-3 atmosphere = 133.3 Pa
1 bar = 105 Pa
X = 0.453 ∗ P ∗ b ∗ n ∗ D2 √(1 − μ2 ) * Et 3
Where P = Pressure (N/m2)
b = radius of each corrugation (m)
n = number of semi-circular corrugations
t = thickness of wall (m)
D= mean diameter (m)
E = modulus of elasticity (N/m2)
µ = Poisson`s ratio
3 P
Ymax = ∗ 3 ∗ R4 (1 − μ2 )
16 Et
a∗h2
P1 =
V1
Where, a = Area of cross section of the measuring capillary
h = Height of compressed gas in the measuring capillary
V1 = Bulb volume including the volume of measuring capillary
= V + (a * hc)
hc = Height of measuring capillary tube
dR
R 2
2.6.1 Sensitivity = = =α
𝑃 𝐸
α = Pressure coefficient of resistance of wire material
Bridgeman gauge using manganin element has a pressure coefficient of resistance of 2.3 x 10 -11 m2/N
3. Temperature Measurement
3.1 For an Electric resistance thermometer:
Over a limited temperature range around 0 0 C (273 K), the following linear relationship can be applied:
Rt = Ro (1 + α t)
Where, α = the temperature coefficient of resistance of material in (Ω/Ω)/ 0C,
Ro = resistance at 0 0C,
t = temperature relative to 0 0C.
Rt = resistance at t 0C,
4. Strain Measurement
Where, L = Length of the electrical conductor, A = Cross sectional area of the conductor
ΔD
D
Poisson’s ratio = μ = Lateral strain/Longitudinal strain = - ΔL
𝐿
Ei
The value of unbalanced current IG can be calculated for any change in strain gauge resistance R 1.
Taking a special case when R1 = R2 = R3 = R4 and if R1 changes to R1 + ∆R1,
IG = (- Ei Gf Є1)/ [4 (R1 + RG)] where Є1 is the strain which causes ∆R1, and Є1 = ∆R1/ R1 Gf
If all the four arms have strain gauge whose resistances change due to strains, it can be shown that
IG = {(- Ei)/ [4 (R1 + R4)]} [∆R1/R1 - ∆R2/R2 + ∆R3/R3 - ∆R4/R4]
4.3 The measuring instrument has an infinitely high internal resistance (for an infinite impedance at
output), the change in output-voltage due to applied strain is given by,
A tensile compressive load cell uses four strain gauges each mounted at 900 to each other on a steel
cylinder.
The change in resistance of the strain gauges are measured by the output voltage of the wheat stone
bridge is given by,
[2(1 + μ) ∗ Gauge factor ∗ Є ∗ E]
Eo =
4
2πNT
Power =
60
3
𝑓𝑠 π𝑑
Also, Torque, T= where fs = Shear stress induced in the shaft
16
If F is the force at the support, the torque T transmitted by the dynamometer is given by T= F * L.
Thus by measuring the force F, the transmitted torque is measured.
The power transmitted can be calculated from the torque, using the equation
P=ω*T
Where P is the power (W), T the torque (N-m) and ω the angular speed (rad /s)
One end of the rope is connected to a mass while the other end is connected to a spring balance
is wound on a pulley with radius r.
2πN(W – S) ∗ r
The power is, P= where N is the speed in rpm
60
𝐴2 2(𝑃1 −𝑃2 )
𝑄 = 𝐶𝑑 × ×√
2 𝜌
√1−(𝐴2 )
𝐴1
2
𝑣 = √𝜌 × (𝑃2 − 𝑃1 )
6.3 Rotameter:
(𝐴𝑡 − 𝐴𝑓 ) 𝑉𝑓 𝜌𝑓 − 𝜌𝑓𝑓
𝑄 = 𝐶𝑑 ∗ √2𝑔 ∗ √ ∗√
2 𝐴𝑓 𝜌𝑓𝑓
√1 − (𝐴𝑡 − 𝐴𝑓 )
[ (𝐴𝑡 )2 ]
3
𝑖 = 0.45 √𝐷 + 0.001𝐷
where, D – Diameter in mm: D = √𝐷𝑖 × 𝐷𝑓
Di – Initial value of Standard Diameter range in mm
Df – Final value of Standard Diameter range in mm
The various standard diameter main ranges specified by I.S.I, are: 1—3, 3—6, 6—10, 10—18, 18—30,
30—50, 50—80, 80—120, 120—180, 180—250, 250—315, 315-400, 400—500, 500—630, 630—800,
800—1000, 1000—1250, 1250—1600, 1600—2000, 2000—2500, 2500—3150 mm
Tolerance
IT5 IT6 IT7 IT8 IT9 IT10 IT11 IT12 IT13 IT14 IT15 IT16
Grade
Magnitude 7i 10i 16i 25i 40i 64i 100i 160i 250i 400i 640i 1000i
Interference fit:
Maximum interference = LL of hole – UL of shaft
Minimum interference = UL of hole – LL of shaft
Transition fit:
Maximum clearance = UL of hole – LL of shaft
Maximum interference = LL of hole – UL of shaft
Here, UL is the Upper limit size & LL is the Lower limit size
8. Gauges
Gauge Makers Tolerance (GMT):
1
GMT = X Work Tolerence
10
1
Wear Allowance = X GMT
10
Table 8.1: Slip gauges – M112 Set
Range (mm) Step (mm) No. of Pieces
1.001 to 1.009 0.001 09
1.010 to 1.490 0.010 49
0.50 to 24.50 0.500 49
25, 50, 55, 100 25 04
1.0005 - 01
Total 112
2π
1 minute of arc = radians
360×60
T = M - 2d
d = Diameter of the best-size wire (mm)
M = Dimension over the wires (mm)
p θ θ
P = [ Cot ] − d[Cosec − 1]
2 2 2
p = Pitch of thread (mm)
θ = Angle of thread
θ p θ
E = M − d [1 + Cosec ] + [ Cot ]
2 2 2
M = Dimension over the wires (mm)
p = Pitch of thread (mm)
θ = Angle of thread
p
db = 2 [ ] Sec θ
4
1 𝑟𝑣 = compression ratio
𝑃𝑚 = (𝛾−1)(𝑟 −1) × 𝑝1 × 𝑟𝑣 (𝛽 − 1)(𝑟𝑣 𝛾−1 − 1)
𝑣 𝛽 = explosion ratio
𝑝1 = Initial Pressure (Pa)
1.2.2 Diesel cycle
𝑝 ×𝑟 𝛾 𝑟𝑣 = compression ratio
1 𝑣
𝑃𝑚 = (𝛾−1)(𝑟 [𝑟 1−𝛾 (1 − 𝛼 𝛾 ) + 𝛾(𝛼 − 1)]
−1) 𝑣
𝑣 𝛼 = cut-off ratio
𝑝1 = Initial Pressure (Pa)
1.2.3 Dual cycle
𝑝 ×𝑟 𝛾 𝑟𝑣 = compression ratio
1 𝑣
𝑃𝑚 = (𝛾−1)(𝑟 [(𝛽 − 1) + 𝛾𝛽(𝛼 − 1) − 𝑟𝑣 1−𝛾 (𝛽𝛼 𝛾 −
−1)
𝑣 𝛼 = cut-off ratio
1)] 𝛽= explosion ratio
𝑝1 = Initial Pressure (Pa)
1 𝑟𝑝 = Pressure ratio
𝜂 = 1− 𝛾−1
(𝑟𝑝 ) 𝛾
ℎ1 −ℎ2𝑠 𝑇1 −𝑇2𝑠
𝜂𝑐 = =
ℎ1 −ℎ2𝑎 𝑇1 −𝑇2𝑎
𝑛−1
1.1.2 𝑛 𝑃2 𝑛
𝑉𝑒𝑠 =Effective Swept volume in (m3/kg)
𝑛
𝑊 = 𝑛−1 𝑃1 𝑉𝑒𝑠 [(𝑃 ) − 1] = 𝑛−1 𝑅[𝑇2 − 𝑇1 ] J
1
1.1.3
3.2 Two-stage Compression
𝑛−1 𝑛−1
𝑛 𝑃 𝑛 𝑛 𝑃 𝑛
𝑊 = 𝑛−1 𝑃1 𝑉1 [(𝑃2 ) − 1] + 𝑛−1 𝑃2 𝑉2′ [(𝑃3 ) − 1] J P3=Discharge pressure (Pa)
1 2
P2=Intermediate pressure (Pa)
V1= Suction volume of first stage (m3)
V2′ = Suction volume of second stage (m3)
𝑛−1 𝑛−1
𝑛 𝑃2 𝑛 𝑃3 𝑛
𝑊 = 𝑛−1 𝑃1 𝑉1 [(𝑃 ) + (𝑃 ) − 2] J
1 2
𝑃2 𝑃3
= = 𝑍 =Pressure ratio per stage
𝑃1 𝑃2
Or, 𝑃2 = √𝑃1 𝑃3
Also, for ‘x’ stages of compressor, the pressure ratio per stage
is given by,
1
𝑃𝑥+1 𝑥
𝑍=( )
𝑃1
𝑛−1
2𝑛 𝑃2 𝑛
𝑊 = 𝑛−1 𝑃1 𝑉1 [(𝑃 ) − 1] J, for 2-stages. P3=Discharge pressure (Pa)
1
P1=Suction pressure (Pa)
Or P2=Intermediate pressure (Pa)
V1=Suction volume (m3)
𝑛−1
2𝑛 𝑃 2𝑛
𝑊 = 𝑛−1 𝑃1 𝑉1 [(𝑃3 ) − 1] J, for 2-stages.
1
𝑛−1
𝑥𝑛 𝑃𝑥+1 𝑥𝑛
𝑊 = 𝑛−1 𝑃1 𝑉1 [( ) − 1] J
𝑃1
𝐵𝑃
𝜂𝑚 = ( 𝐼𝑃 )100 %
5. Refrigeration
Bell-Coleman cycle
𝛾 = Specific heat ratio.
1 P2 = Pressure after compression (Pa)
𝐶𝑂𝑃 = 𝛾−1 P1 = Pressure before compression (Pa)
𝑃2 𝛾
[( ) ]−1
𝑃1
𝑇𝑐 2
= Tc = Critical temperature (K)
𝑇𝑜 𝛾 + 1
ρc = Critical density (kg/m3)
1
Pc = Critical pressure (Pa)
𝜌𝑐 2 𝛾−1
=[ ]
𝜌𝑜 𝛾+1
𝛾
𝑃𝑐 2 𝛾−1
=[ ]
𝑃𝑜 𝛾+1
2.2 Circle:
2.2.1 Parametric equation of a circle
Xn = Xc + R * Cos u where 0 ≤ u ≤ 2 and Xn, Yn, and Zn are
the coordinates of any point Pn on the
Yn = Yc + R * Sin u
circumference of circle with radius = R;
Zn = Zc
Xn+1 = Xc + (Xn – Xc) Cos u – (Yn – Yc) Sin u u = incremental angle subtended at the
center of the circle;
Yn+1 = Yc + (Yn – Yc) Cos u + (Xn – Xc) Sin u
Xn, Yn, and Zn = the coordinates of any
Zn+1 = Zc,
point Pn on the circumference of circle
with radius = R;
Xc,Yc and Zc = centre point coordinates
of the circle.
2.3 Ellipse:
2.3.1 Parametric equation of an Ellipse
Xn = Xc + A * Cos u where 0 ≤ u ≤ 2
Yn = Yc + B * Sin u Xn, Yn, and Zn = the coordinates of any
point Pn on the circumference of the
Zn = Zc
ellipse.
A = half of the length of the major axis
and
B = half of the length of the minor axis;
2.3.2 Recursive equation of an Ellipse
Xn+1 = Xc + (Xn – Xc) Cos u – A/B (Yn – Yc) Sin u where 0 ≤ u ≤ 2and u= incremental
angle
Yn+1 = Yc + (Yn – Yc) Cos u + B/A (Xn – Xc) Sin u
Xn, Yn, and Zn = the coordinates of any
Zn+1 = Zc,
point Pn on the circumference of ellipse.
A = half of the length of the major axis
and
B = half of the length of the minor axis;
u = the angle made by any point P on
the circumference with reference to x-
axis
Xc,Yc and Zc = the center point
coordinates of the ellipse;
2.3.3 If the major axis of the ellipse is inclined at an angle α relative to the X axis, then the
parametric relationship changes to
Yn = Yc + A Cos u Sin + B Sin u Cos A = half of the length of the major axis
and
Zn = Zc B = half of the length of the minor axis;
u = the angle made by any point P on
the circumference with reference to x-
axis
= angle made by the ellipse with
reference to X axis
Xc,Yc and Zc = the center point
coordinates of the ellipse
Yn = Yc + A Cos (un + u) Sin + B Sin (un + u) A = half of the length of the major axis
Cos and
UH = √𝑋𝐻𝑊 /𝐴 UL = √𝑋𝐿𝑊 /𝐴
P(u) = (2u3 – 3u2 +1) P0 + (-2u3 + 3u2) P1 + (u3 – 2u2 + u) P0& P1 = position vector at the
P0 + (u3-u2) P1 end points of the curve
P0’ & P1’ = tangent vector at the
end points of the curve
3.2.2 The tangent vectors at the first and last endpoints are: P1& P0 = position vector at the
first two points on the curve
P0 = n(P1 – P0)
Pn – Pn-1 = position vector at the
P1 = n(Pn – Pn-1)
last two points on the curve
n (degree of the curve) =
number of points defining the
curve – 1
3.3 B-spline curve:
3.3.1 Position vector of any point on the B- Spline curve is given Pi = position vector at ith point
by on the curve
P(u) = ∑𝑛
𝑖=0 Pi Ni,k (u) Ni,k (u) = Bspline coefficient
0≤u≤1
4.1.1 The position vector of any point P on the plane surface P0, P1 and P2 = points on the
surface that are non collinear
P(u,v) = P0 + u(P1 – P0) + v(P2 – P0)
0 ≤ u ≤ 1, 0≤v≤1
4.1.2 The tangent vectors at point P on the plane surface P0, P1 and P2 = points on the
surface that are non collinear
Pu(u,v) = (P1 – P0)
0 ≤ u ≤ 1, 0≤v≤1
Pv(u,v) = (P2 – P1)
6.2 List of ‘M’ codes and their function in CNC turning operation
M CODES FUNCTION M CODES FUNCTION
M00 Optional program stop automatic M01 Optional program stop request
M02 Program end M03 Spindle ON clock wise (CW)
M04 Spindle ON counter clock wise M05 Spindle stop
(CCW)
M06 Tool change M07 Mist coolant ON (coolant 1 ON)
M08 Flood coolant ON (coolant 2 ON) M09 Coolant OFF
M19 Spindle orientation M30 End of program, Reset to start
M98 Sub program call M99 Sub program end
7.2 Step turning operation using the box turning cycle (G90 cycle)
The syntax is: X is the diameter to which the tool
movement is being made
G90 X… Z… F…
Z is the Z axis coordinate to which the tool
movement is being made
F is the feed rate
7.3 Taper turning cycle (G90 cycle)
The syntax is: X is the diameter to which the tool
movement is being made.
G90 X… Z… R… F…
Z is the Z axis coordinate to which the tool
movement is being made.
F is the feed rate being used.
R is the amount of taper =
(Tool diameter at the entry – Tool diameter
at the exit)/ 2
{𝑓} = [𝑘]{𝑞}
f q
f 1 &q 1
f2 q2
1.3 Displacement q at any arbitrary point for a horizontal 1D bar Where q1 and q2 are displacements
element is given by at node 1 and 2 for a two noded bar
element (m)
𝑞 = 𝑁1 𝑞1 + 𝑁2 𝑞2 N1 and N2 are Shape functions for
bar element.
1 1 is the natural coordinate
N1 & N2
2 2
1.4 Stiffness matrix for 1D bar element: A = Cross sectional area of element
(m2)
𝐴𝐸 1 −1 E = Young’s Modulus (Pa or N/m2)
[𝑘] = [ ]
𝑙𝑒 −1 1 le = Length of element (m)
{𝑓} = [ 𝑘] {𝑞}
f q
f 1 &q 1
f2 q2
𝜎 = 𝐸𝜖
f2 x
f2 y
And element displacement vector is given by:
q1x
q
q
1y
q2 x
q2 y
E
C S C Sq
L
1.5 1D steady state heat transfer due to conduction: A = area of cross section of
Element Conduction stiffness matrix is given by: element (m2)
k = Thermal conductivity of material
Ak 1 1 (W/m.K)
[ kc ] 1 1 L = Length of element (m)
L
{Qc} = Thermal load vector
Element equations can be written as {T} = Temperature vector (K)
Qc [kc ]T
Where
Qc1
{Qc } Represents the heat transfer by conduction
Qc 2
at node 1 and node 2 of element respectively.
and
T
{T } 1
T2
Represents temperature at node 1 and node 2 of element
respectively.
1.6 Steady state heat transfers due to conduction and {fQ} =Heat source force (W/K)
convection {fq} = Heat flux force (W/K)
{fh} = Convection force (W/K)
Element stiffness equations can be written as {fh}end = End convection force (W/K)
f k T
e e e
A = area of cross section of
element (m2)
Thermal load vector due to conduction and convection is
L = Length of element (m)
given as:
Q = Heat source (W/m3)
q = Heat flux (W/m2)
f fQ f q f h f h end
e
P = Perimeter of element (m)
h = Convective heat transfer
coefficient (W/m2.K)
T = Ambient temperature (K)
k e kc kh kh end
1 1
Ak
[ kc ] 1 1
L
hPL 2 1
[kh ]
6 1 2
1 0
[kh ]end1 hA
0 0
0 0
[kh ]end 2 hA
0 1
𝑣1
𝑞1
𝑞2 (𝑑𝑣⁄𝑑𝑥) q1 & q2 are DOF at node 1 and q3 & q4
1 are DOF at node 2 representing linear
{𝑞} = {𝑞 } = 𝑣2
3 and rotary displacements (or slope)
𝑞4 𝑑𝑣
{( ⁄𝑑𝑥 )2 }
respectively
1 1
𝐻1 = (2 − 3 + 3 ) = (1 − )2(2 + )
4 4
1 1 H1 & H2 are shape functions at node 1
𝐻2 = (1 − − 2 + 3 ) = (1 − )2 (1 + ) and H3 & H4 are shape functions at
4 4
1 1 node 2, representing deflection and
𝐻3 = (2 + 3 − 3 ) = (1 + )2 (2 − ) slope respectively
4 4
1 1
𝐻4 = (−1 − + 2 + 3 ) = (1 + )2 ( − 1)
4 4
Consideration of distributed load ‘p’ when constant, i.e., Uniformly distributed load (UDL):
T
pl pl
2
pl
2
ple
f e e e
e
2 12 2 12
1 𝜐 0
𝐸 𝜐 1 0 E = Young’s Modulus (N/m2)
[𝐷] = ⌊ 1 − 𝜐⌋
2
(1 − 𝜐 ) = Poisson's Ratio
0 0
2
Constitutive matrix for Plane Strain condition:
1−𝜐 𝜐 0
𝐸 𝜐 1−𝜐 0 E = Young’s Modulus (N/m2)
[𝐷] = ⌊ 1 − 2𝜐⌋
(1 + 𝜐)(1 − 2𝜐) = Poisson's Ratio
0 0
2
Stress-Strain relationship:
{} = Stress vector
𝜎𝑥𝑥 {} = Strain vector
{𝜎} = {𝜎𝑦𝑦 } = [𝐷]{𝜀} = [𝐷][𝐵]{𝑞} {q} = Displacement vector
𝜏𝑥𝑦 [B] = Strain-Displacement Interaction
Matrix
𝑦 = 𝑁1 𝑦1 + 𝑁2 𝑦2 + 𝑁3 𝑦3
Displacement components:
‘u’ and ‘v’ are displacement
𝑢 = 𝑁1 𝑞1 + 𝑁2 𝑞3 + 𝑁3 𝑞5 components of any arbitrary point
within the triangular element in x and
𝑣 = 𝑁1 𝑞2 + 𝑁2 𝑞4 + 𝑁3 𝑞6 y directions respectively
e
b1 bx
b b
2 y
b3 Aete bx
{b}
e
b
4 3 by
b5 bx
b6 by
Displacement vector:
Shape Functions:
𝑦 = 𝑁1 𝑦1 + 𝑁2 𝑦2 + 𝑁3 𝑦3 + 𝑁4 𝑦4
Displacement components:
‘u’ and ‘v’ are displacement
𝑢 = 𝑁1 𝑞1 + 𝑁2 𝑞3 + 𝑁3 𝑞5 + 𝑁4 𝑞7 components of any arbitrary point
within the quadrilateral element in x
𝑣 = 𝑁1 𝑞2 + 𝑁2 𝑞4 + 𝑁3 𝑞6 + 𝑁4 𝑞8 and y directions respectively
[ B] B1 B2 B3 B4
Where, B1, B2, B3 & B4 are sub-matrices of [B]
𝜕𝑦 𝜕𝑁𝑖 𝜕𝑦 𝜕𝑁𝑖
( − ) 0
𝜕𝜂 𝜕𝜁 𝜕𝜁 𝜕𝜂
1 𝜕𝑥 𝜕𝑁𝑖 𝜕𝑥 𝜕𝑁𝑖
[𝐵𝑖 ] = 0 ( − )
𝐷𝑒𝑡[𝐽] 𝜕𝜁 𝜕𝜂 𝜕𝜂 𝜕𝜁
𝜕𝑥 𝜕𝑁𝑖 𝜕𝑥 𝜕𝑁𝑖 𝜕𝑦 𝜕𝑁𝑖 𝜕𝑦 𝜕𝑁𝑖
( − ) ( − )
[ 𝜕𝜁 𝜕𝜂 𝜕𝜂 𝜕𝜁 𝜕𝜂 𝜕𝜁 𝜕𝜁 𝜕𝜂 ]
Where i = 1, 2, 3, 4 and,
𝜕𝑥 𝜕𝑦
𝜕𝜁 𝜕𝜁
[𝐽] =
𝜕𝑥 𝜕𝑦
[𝜕𝜂 𝜕𝜂]
Constitutive matrix:
1−𝜗 𝜗 𝜗 0 0 0
𝜗 1−𝜗 𝜗 0 0 0
𝜗 𝜗 1−𝜗 0 0 0
𝐸 (1 − 2𝜗)
[𝐷] = 0 0 0 0 0
(1 + 𝜗)(1 − 2𝜗) 2
(1 − 2𝜗)
0 0 0 0 0
2
(1 − 2𝜗)
[ 0 0 0 0 0
2 ]
Displacement vector:
{𝑞} = {𝑞1 𝑞2 𝑞3 𝑞4 𝑞5 𝑞6 𝑞7 𝑞8 𝑞9 𝑞10 𝑞11 𝑞12 }𝑇 q3i-2, q3i-1 and q3i are DOF’s at each
node (i = 1, 2, 3, 4)
Shape Functions
𝑥 = 𝑁1 𝑥1 + 𝑁2 𝑥2 + 𝑁3 𝑥3 + 𝑁4 𝑥4
𝑦 = 𝑁1 𝑦1 + 𝑁2 𝑦2 + 𝑁3 𝑦3 + 𝑁4 𝑦4
𝑧 = 𝑁1 𝑧1 + 𝑁2 𝑧2 + 𝑁3 𝑧3 + 𝑁4 𝑧4
Displacement components:
𝑢 = 𝑁1 𝑞1 + 𝑁2 𝑞4 + 𝑁3 𝑞7 + 𝑁4 𝑞10
Where, ‘u’, ‘v’ and ‘w’ are
𝑣 = 𝑁1 𝑞2 + 𝑁2 𝑞5 + 𝑁3 𝑞8 + 𝑁4 𝑞11 displacement components of any
arbitrary point within the element in x
𝑤 = 𝑁1 𝑞3 + 𝑁2 𝑞6 + 𝑁3 𝑞9 + 𝑁4 𝑞12 and y directions respectively
where 𝐴1 = 𝐴11 + 𝐴12 + 𝐴13 , 𝐴2 = 𝐴21 + 𝐴22 + 𝐴23 , 𝐴3 = 𝐴31 + 𝐴32 + 𝐴33
𝐴11 𝐴12 𝐴13 𝑦24 𝑧34 − 𝑦34 𝑧24 𝑦34 𝑧14 − 𝑦14 𝑧34 𝑦14 𝑧24 − 𝑦24 𝑧14
1
𝐴
[ 21 𝐴22 𝐴23 ] = [ 𝑧24 𝑥34 − 𝑧34 𝑥24 𝑧34 𝑥14 − 𝑧14 𝑥34 𝑧14 𝑥24 − 𝑧24 𝑥14 ]
𝐴31 𝐴32 𝐴33 det[𝐽] 𝑥24 𝑦34 − 𝑥34 𝑦24 𝑥34 𝑦14 − 𝑥14 𝑦34 𝑥14 𝑦24 − 𝑥24 𝑦14
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝜉 𝜕𝜉 𝜕𝜉 𝑥14 𝑦14 𝑧14
𝜕𝑥 𝜕𝑦 𝜕𝑧
Jacobian [𝐽] = 𝜕𝜂 𝜕𝜂 𝜕𝜂
= [𝑥24 𝑦24 𝑧24 ]
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝑥34 𝑦34 𝑧34
[𝜕𝜁 𝜕𝜁 𝜕𝜁 ]
𝑏𝑥
𝑏𝑦 bx, by, bz are body force component in
x, y & z direction respectively at each
𝑏𝑧
𝑉 node
𝑒 𝑏 𝑥
{𝑏}𝑒 =
4 𝑏𝑦
𝑏𝑧
⋮
{ 𝑏𝑧 }12×1
𝐹 = ∑({𝑏}𝑒 + {𝑇}𝑒 + 𝑃)
𝑒
Shape Functions:
1
𝑁𝑖 = (1 + 𝜉𝜉𝑖 )(1 + 𝜂𝜂𝑖 )(1 + 𝜁𝜁𝑖 )
8
1 where i = 1 to 8 and (ξi, ηi, ζi)
𝑁1 = (1 − 𝜉)(1 − 𝜂)(1 − 𝜁)
8 represents the coordinates of node i
1 of the element in the (ξ, η, ζ) system
𝑁2 = (1 + 𝜉)(1 − 𝜂)(1 − 𝜁)
8
1
𝑁3 = (1 + 𝜉)(1 + 𝜂)(1 − 𝜁)
8
1
𝑁4 = (1 − 𝜉)(1 + 𝜂)(1 − 𝜁)
8
1
𝑁5 = (1 − 𝜉)(1 − 𝜂)(1 + 𝜁)
8
1
𝑁6 = (1 + 𝜉)(1 − 𝜂)(1 + 𝜁)
8
1
𝑁7 = (1 + 𝜉)(1 + 𝜂)(1 + 𝜁)
8
1
𝑁8 = (1 − 𝜉)(1 + 𝜂)(1 + 𝜁)
8
𝑥 = 𝑁1 𝑥1 + 𝑁2 𝑥2 + 𝑁3 𝑥3 + ⋯ + 𝑁8 𝑥8
𝑦 = 𝑁1 𝑦1 + 𝑁2 𝑦2 + 𝑁3 𝑦3 + ⋯ + 𝑁8 𝑦8
𝑧 = 𝑁1 𝑧1 + 𝑁2 𝑧2 + 𝑁3 𝑧3 + ⋯ + 𝑁8 𝑧8
Displacement components:
𝑢 = 𝑁1 𝑞1 + 𝑁2 𝑞4 + 𝑁3 𝑞7 + ⋯ + 𝑁8 𝑞22
𝑣 = 𝑁1 𝑞2 + 𝑁2 𝑞5 + 𝑁3 𝑞8 + ⋯ + 𝑁8 𝑞23
𝑤 = 𝑁1 𝑞3 + 𝑁2 𝑞6 + 𝑁3 𝑞9 + ⋯ + 𝑁8 𝑞24
[𝐵] = [𝐵1 𝐵2 𝐵3 𝐵4 𝐵5 𝐵6 𝐵7 𝐵8 ]
and i = 1 to 8;
Strain Vector:
u
xx
x
v
yy
y
w
zz z
{ }
u v
xy y x
v w
yz z y
w u
zx
x z
Jacobian Matrix:
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝜉 𝜕𝜉 𝜕𝜉
𝜕𝑥 𝜕𝑦 𝜕𝑧
[𝐽] =
𝜕𝜂 𝜕𝜂 𝜕𝜂
𝜕𝑥 𝜕𝑦 𝜕𝑧
[𝜕𝜁 𝜕𝜁 𝜕𝜁]
−1 −1
1D Problem:
𝜌 𝐴 𝑙𝑒 2 1
[𝑚]𝑒 = ⌊ ⌋
6 1 2
Truss Element:
2 0 1 0
𝜌 𝐴 𝑙𝑒 0 2 0 1
[𝑚]𝑒 = [ ]
6 1 0 2 0
0 1 0 2
Beam Element:
1D Problem:
𝜌 𝐴 𝑙𝑒 1 0
[𝑚]𝑒 = ⌊ ⌋
2 0 1
Truss Element:
1 0 0 0
𝜌 𝐴 𝑙𝑒 0 1 0 0
[𝑚]𝑒 = [ ]
2 0 0 1 0
0 0 0 1
Beam Element:
1 0 0 0
𝜌 𝐴 𝑙𝑒 0 0 0 0
[𝑚]𝑒 = [ ]
2 0 0 1 0
0 0 0 0
P P = Power (W)
Power coefficient = Density (kg/m3)
N 3 D5
VD V = Flow velocity (m/s)
Reynolds number = Absolute viscosity of the fluid (kg/ms)
Q Qu = Unit discharge
Qu Q = Volume flow rate (m3/s)
h h = Head (m)
N Nu = Unit speed
Nu N = Rotational speed (RPM)
h
P Pu = Unit power
Pu 3/ 2 P = Power (kW)
h
1-1 D2
n
1 = Efficiency of model
For turbines: 2
1-2 D1
= Efficiency of prototype
D1 = Diameter of the model (m)
D2 = Diameter of the prototype (m)
n = Empirical constant
n Re1 = Reynolds number of flow for the model
1-1 Re2
For pumps: Re2 = Reynolds number of flow for the prototype
1- 2 Re1
R
u 2
1 u22 w22 w12 R = Degree of reaction
C 2
1 C22 u12 u22 w22 w12
tanβ1' +tan β 2' c 1' = Blade angle at inlet with respect to axial
c
R a a tanβ mean
'
direction (deg)
u 2 u
2' = Blade angle at exit with respect to axial
direction (deg)
C12 C22
C12 RC22
ln
T2 p = Polytropic efficiency
T1
p 1
T1 = Static temperature at the inlet of turbine (K)
T2 = Static temperature at the exit of turbine (K)
ln
p2 p1 = Static pressure at the inlet of turbine (Pa)
p1 p2 = Static pressure at the exit of turbine (Pa)
p2
rp = Stage pressure ratio =
1
p
p1
T 1 = Ratio of specific heats
1 2 1
T1 rp s = Stage efficiency
s
T 1
rpov = Overall pressure ratio = (rp)N
1 2s 1 N = Number of stages
T1 1 ov = Overall efficiency
rp
n = Polytropic index
n 1
p
n 1
1
N
1 r
ov p
N
1 1
1 1 1 rp
s
n 1
p
n 1
Cm Flow coefficient
Cm = Meridional component of absolute velocity
2 ghav
(m/s)
p2 = Static pressure at turbine exit (Pa)
u1C1u u2C2u
hyd C3 = Velocity of water at the exit of draft tube
ghav (m/s)
hfdr = Frictional head loss in the draft tube (m)
p2 p C 2 C32 hs = Distance between turbine exit and water
atm 2 h fdr hs level in the tail race (m)
g g 2g C2 = Absolute velocity of fluid at runner exit
(m/s)
C22 C32 C3 = Absolute velocity of fluid at draft tube exit
h fdr (m/s)
2g
dr dr Draft tube efficiency
C22
2g
Ca Flow coefficient
2 ghav Blade speed ratio
D = Diameter of the impeller (m)
u
dh = Hub diameter (m)
2 ghav
Q
4
D 2
dh2 2 ghav
u C1u C2u
hyd
ghav
7. Steam Turbines
7.1 Single stage impulse turbine
C12 u
Blade speed ratio
2 C1
1 = Flow angle at inlet (deg)
Depreciation
i. Straight Line Depreciation Dn = Annual depreciation amount
𝑃𝑢𝑟𝑐ℎ𝑎𝑠𝑒 𝑃𝑟𝑖𝑐𝑒−𝑆𝑎𝑙𝑣𝑎𝑔𝑒 𝑉𝑎𝑙𝑢𝑒 I = Purchase price of the equipment
Annual Depreciation = 𝐷𝑛 = n = Number of Years
𝑌𝑒𝑎𝑟𝑠 𝑜𝑓 𝑢𝑠𝑒𝑓𝑢𝑙 𝑙𝑖𝑓𝑒
Book Value = 𝐼 − (𝑛 ∗ 𝐷𝑛 )
Ratio Analysis
Types of Financial Ratios
i. Liquidity Ratios:
Current assets
𝐶𝑢𝑟𝑟𝑒𝑛𝑡 𝑅𝑎𝑡𝑖𝑜 =
Current liabilities
Structural Ratio
Total debt or long term debt
𝐷𝑒𝑏𝑡 𝑡𝑜 𝑒𝑞𝑢𝑖𝑡𝑦 𝑟𝑎𝑡𝑖𝑜𝑠 =
Shareholder’s equity
1. Inventory turnover
Cost of goods sold
𝐼𝑛𝑣𝑒𝑛𝑡𝑜𝑟𝑦 𝑡𝑢𝑟𝑛𝑜𝑣𝑒𝑟 =
Average inventory
OR
Sales
In the absence of data, Inventory turnover =
Closing Inventory
Total sales
In the absence of data, Debtors turnover =
(Debtors + Bills receivable)
3. Asset Turnover
Costs of goods sold
Fixed Asset Turnover =
Avg. fixed assets
Gross profit
1. Gross profit margin = ∗ 100
Sales
Return on Investment
Profits of firm to its investment
Return on Assets
Return on equity
1. Introduction
1.1 Heat flow rate (W)
Tw = Wall temperature
T∞ = Fluid temperature
h = Convective heat transfer coefficient (W/m2K)
1
𝑅=
ℎ𝐴
1.1.3 Radiation
1
𝑅=
𝜎𝐴𝐹(𝑇1 + 𝑇2 2 )(𝑇1 + 𝑇2 )
2
1 1
1 1 1 𝑇𝑖 −𝑇𝑟 ( − )
𝑟1 𝑟
𝑅= ( − 𝑟 ), = 1 1
4𝜋𝑘 𝑟1 2 𝑇𝑖 −𝑇𝑜 (𝑟 −𝑟 )
1 2
1 1 𝐿1 𝐿2 𝐿3 1
𝑅= [ + + + + ]
𝐴 ℎ𝑎 𝑘1 𝑘2 𝑘3 ℎ𝑏
1 1 1 𝑟 1 𝑟 1 𝑟 1
𝑅= [ + 𝑙𝑛 ( 2) + 𝑙𝑛 ( 3) + 𝑙𝑛 ( 4) + ]
2𝜋𝐿 ℎ𝑎 𝑟1 𝑘1 𝑟 𝑘 1 𝑟 𝑘2 𝑟 ℎ 2 3 3 𝑏 𝑟4
1 1 1 1 1 1 1 1 1
𝑅= [ 2
+ (𝑟 − 𝑟 ) + 𝑘 (𝑟 − 𝑟 ) + ℎ 2
]
4𝜋 ℎ𝑖 𝑟𝑖 𝑘1 1 2 2 2 3 𝑜 𝑟3
𝑅 = 𝑅1 + 𝑅2 + 𝑅3 + 𝑅4
𝑅 𝑅 𝑅𝐷 𝑅𝐸 𝑅𝐹
𝑅1 = 𝑅𝐴 , 𝑅2 = 𝑅 𝐵+𝑅𝐶 , 𝑅3 = 𝑅 , 𝑅4 = 𝑅𝐺
𝐵 𝐶 𝐷 𝑅𝐸 +𝑅𝐸 𝑅𝐹 +𝑅𝐹 𝑅𝐷
𝐿𝐴 𝐿𝐺
𝑅𝐴 = 𝑘 ,… 𝑅𝐺 = 𝑘 , 𝐴 = 𝐴𝐴 = 𝐴𝐵 + 𝐴𝐶 = 𝐴𝐷 + 𝐴𝐸 + 𝐴𝐹 = 𝐴𝐺
𝐴 𝐴𝐴 𝐺 𝐴𝐺
𝑞̇ 𝑇𝑥 −𝑇𝑜 𝑥 2 𝑞̇
𝑇𝑜 = 𝑇𝑤 + 𝐿2 , =( ) , 𝑇𝑥 = 𝑇𝑜 − 𝑥 2 , 𝑞𝑥 = 𝑞̇ 𝑥
2𝑘 𝑇𝑤 −𝑇𝑜 𝐿 2𝑘
𝑞̇ 𝐿
𝑇𝑤 = 𝑇∞ +
ℎ
𝑞̇ 𝑥 1
𝑇𝑥 = 2𝑘 (𝐿2 − 𝑥 2 ) + 2𝐿 (𝑇𝑤 2 − 𝑇𝑤 1 ) + 2 (𝑇𝑤 1 + 𝑇𝑤 2 )
𝑘 𝑞̇ 𝐿2 𝑘 2 1
𝑥𝑚𝑎𝑥 = 2𝑞̇ 𝐿 (𝑇𝑤 2 − 𝑇𝑤 1 ) , 𝑇𝑚𝑎𝑥 = + 8𝑞̇ 𝐿2 (𝑇𝑤 1 − 𝑇𝑤 2 ) + 2 (𝑇𝑤 1 + 𝑇𝑤 2 )
2𝑘
𝑞̇ 𝑅 2 𝑇𝑟 −𝑇𝑤 𝑟 2 𝑞̇
𝑇𝑜 = 𝑇𝑤 + , = 1 − (𝑅) , 𝑇𝑟 = 𝑇𝑤 + 4𝑘 (𝑅 2 − 𝑟 2 ) , 𝑞𝑟 = 𝑞̇ 𝜋𝑟 2 𝐿
4𝑘 𝑇𝑜 −𝑇𝑤
𝑘
𝑟𝑐 =
ℎ
𝑞̇ 𝑅
𝑇𝑤 = 𝑇∞ +
2ℎ
𝑇𝑖 −𝑇𝑟 𝑟 2 𝑞̇ 4 𝑞̇
= (𝑅 ) , 𝑇𝑖 = 𝑇𝑜 + 6𝑘 𝑅𝑜 2 , 𝑞𝑟 = 𝑞̇ 3 𝜋𝑟 3 , 𝑇𝑟 − 𝑇𝑜 = 6𝑘 (𝑅𝑜 2 − 𝑟 2 )
𝑇𝑖 −𝑇𝑜 𝑜
𝑞̇ 2𝑘
𝑇𝑟 = 𝑇𝑖 − 6𝑘 𝑟 2 , 𝑟𝑐 = ℎ
𝑞̇ 𝑅𝑜
𝑇𝑜 = 𝑇∞ +
3ℎ
ε = Fin effectiveness
𝑚 = √(ℎ𝑃 ⁄𝑘𝐴) θ1 = Temperature difference (T1-Ta)
θ2 = Temperature difference (T2-Ta)
2.3.1 Long fin (TL = T∞)
𝑇−𝑇∞ 1
= 𝑒 −𝑚𝑥 , 𝑄 = (𝑇𝑏 − 𝑇∞ )(ℎ𝑃𝑘𝐴)0.5 , 𝜂𝑓𝑖𝑛 = 𝑚𝐿 , 𝜀 = (𝑃𝑘⁄ℎ𝐴)0.5
𝑇𝑏 −𝑇∞
𝑇−𝑇∞ 𝑐𝑜𝑠ℎ(𝑚[𝐿−𝑥])
= , 𝑄 = (𝑇𝑏 − 𝑇∞ )(ℎ𝑃𝑘𝐴)0.5 𝑡𝑎𝑛ℎ(𝑚𝐿)
𝑇𝑏 −𝑇∞ 𝑐𝑜𝑠ℎ(𝑚𝐿)
𝑡𝑎𝑛ℎ(𝑚𝐿)
𝜂𝑓𝑖𝑛 = , 𝜀 = 𝑡𝑎𝑛ℎ(𝑚𝐿)(𝑃𝑘⁄ℎ𝐴)0.5
𝑚𝐿
2.3.3 Short fin (convective loss from end)
𝜃1 𝑠𝑖𝑛ℎ(𝑚[𝐿−𝑥])+𝜃2 𝑠𝑖𝑛ℎ(𝑚𝑥)
𝑇 − 𝑇𝑎 = 𝑠𝑖𝑛ℎ(𝑚𝐿)
𝑐𝑜𝑠ℎ(𝑚𝐿)−1
𝑄 = (𝜃1 + 𝜃2 ) (ℎ𝑃𝑘𝐴)0.5
𝑠𝑖𝑛ℎ(𝑚𝐿)
𝛿ℎ𝑥 = 𝛿𝑇𝑥
𝑁𝑢𝐿 = 1.25𝑁𝑢𝑥
L = Lc
0.25
Laminar: Gr Pr<109 𝑁𝑢 = 0.59(𝐺𝑟𝑃𝑟)
𝑮𝒓𝑫 𝑷𝒓 C m
102 to 104 0.85 0.188
104 to 107 0.48 0.250
107 to 1012 0.125 0.333
𝐴𝑠⁄
𝐿𝑐 = 𝑃
𝐿𝑐 = 𝐷
𝐿𝑐 = 𝐷
Plate:
Upper surface heated: 𝑁𝑢 = 0.54(𝐺𝑟𝑃𝑟)0.25
Lower surface heated: 𝑁𝑢 = 0.27(𝐺𝑟𝑃𝑟)0.25
Sphere: 𝑁𝑢 = 0.6(𝐺𝑟𝑃𝑟)0.25
−𝑙𝑛[1−(1+𝐶)𝜀]
𝐶𝑚𝑖𝑛 = Minimum fluid capacity rate (W/K) = 𝐶𝑝 𝑚̇
𝑁= (1+𝐶) 𝐶𝑚𝑎𝑥 = Maximum fluid capacity rate (W/K) = 𝐶𝑝 𝑚̇
4.3 Overall heat transfer coefficient 𝑅𝑓𝑖 = Inner surface fouling resistance (m 2K/W)
5.1 Introduction
𝑐 𝜆 = Wavelength (m)
5.1.1 Wave length, 𝜆=𝜈 = Speed of light = 3*108 m/s
𝑐
5.1.2 Radiation property, 𝜌 + 𝛼 + 𝜏 =1 𝜈 = Frequency
𝜀𝜆 = Monochromatic emissivity
5.1.5 Stefan-Boltzmann law for black bodies,
𝛼𝜆 = Monochromatic absorptivity
𝐸𝑏 = 𝜎𝑇 4 𝜀 = Emissivity = 𝐸 ⁄𝐸𝑏
𝐶 𝜆−5 𝐸 = Emissive power (W/m 2)
5.1.6 Plank’s law, 𝐸𝑏𝜆 = 𝑒𝑥𝑝(𝐶1 ⁄𝜆𝑇)−1
2
𝐸𝑏 = Emissive power of black body (W/m2)
𝐴 = Radiated area (m 2)
5.2.3 Net heat exchange between two black
𝐽 = Radiosity (W/m 2)
bodies,
𝐺 = Irradiation (W/m 2)
𝑄12 = 𝐴1 𝐹1−2 𝜎(𝑇1 4 − 𝑇2 4 )
𝐹1−2 = Shape factor
𝐹2−1 = Shape factor