Dynamic Force Analysis
Dynamic Force Analysis
Dr. A. B. Andhare
Velocity and Acceleration of Piston in
Reciprocating Engine Mechanism
A
l r l = length of connecting rod
B B1 y
β θ O r = Crank radius
IDC A1 ODC
x x = displacement of piston from inner dead center
𝑠𝑖𝑛2 𝜃
x = B1B = BO – B1O 𝑐𝑜𝑠β = 1 − 𝑠𝑖𝑛2 𝛽 𝑐𝑜𝑠𝛽 = 1−
𝑛2
= BO – (B1A1 + A1O) 𝑦2
𝑐𝑜𝑠𝛽 = 1− 2 1
= (l+r) – (l cosβ + r cosθ) 𝑙 𝑐𝑜𝑠𝛽 = 𝑛2 − 𝑠𝑖𝑛2 𝜃
𝑛
Let l/r = n, then l = n r
(𝑟𝑠𝑖𝑛𝜃)2 𝑥 = 𝑟[ 𝑛 + 1 − ( 𝑛2 − 𝑠𝑖𝑛2 𝜃 + 𝑐𝑜𝑠𝜃)]
x = (n r + r) – (n r cosβ + r cosθ) 𝑐𝑜𝑠𝛽 = 1−
𝑙2
= r [(n+1) – (n cosβ + cosθ)] 𝑥 = 𝑟[ 1 − 𝑐𝑜𝑠𝜃 + (𝑛 − 𝑛2 − 𝑠𝑖𝑛2 𝜃) ]
Velocity of piston:
𝑑𝑥 𝑑𝑥 𝑑𝜃 𝑑 𝑑𝜃
𝑣= = = 𝑟 1 − 𝑐𝑜𝑠𝜃 + 𝑛 − (𝑛2 −𝑠𝑖𝑛2 𝜃)1/2
𝑑𝑡 𝑑𝜃 𝑑𝑡 𝑑𝜃 𝑑𝑡
1
−
2 𝑠𝑖𝑛2𝜃
= 𝑟 (0 + 𝑠𝑖𝑛𝜃 + 0 − 1/2 𝑛2 − 𝑠𝑖𝑛2 𝜃 (−2𝑠𝑖𝑛θ𝑐𝑜𝑠θ)]ω 𝑣 = 𝑟𝜔 𝑠𝑖𝑛𝜃 +
2 𝑛2 − 𝑠𝑖𝑛2 𝜃
Theory of Machines - II: Dr. A. B. Andhare 3
Velocity and Acceleration of Piston in
Reciprocating Engine Mechanism
𝑠𝑖𝑛2𝜃
𝑣 = 𝑟𝜔 𝑠𝑖𝑛𝜃 +
2 𝑛2 − 𝑠𝑖𝑛2 𝜃
𝑠𝑖𝑛2𝜃
𝑣 = 𝑟𝜔 𝑠𝑖𝑛𝜃 +
2𝑛
𝑣 = 𝑟𝜔𝑠𝑖𝑛𝜃
𝑠𝑖𝑛𝜃 2
1 2 3
−2
As y = l sinβ = r sinθ 𝑠𝑖𝑛𝛽 = 𝛼𝑐 = 𝜔 −𝑐𝑜𝑠𝜃 𝑛 − 𝑠𝑖𝑛2 𝜃 −2𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜃 + 𝑛2 − 𝑠𝑖𝑛2 𝜃 −1/2
−𝑠𝑖𝑛𝜃
𝑛 2
Differentiating with time to get ωc, the 𝑐𝑜𝑠 2 𝜃 − (𝑛2 − 𝑠𝑖𝑛2 𝜃)
angular velocity of connecting rod 𝛼𝑐 = 𝜔2 𝑠𝑖𝑛𝜃
(𝑛2 −𝑠𝑖𝑛2 𝜃)3/2
𝑑𝛽 1 𝑑𝜃 𝑑𝛽 𝑐𝑜𝑠𝜃
2
𝑛2 − 1
𝑐𝑜𝑠𝛽 = 𝑐𝑜𝑠𝜃 𝑜𝑟 = 𝜔 𝛼𝑐 = −𝜔 𝑠𝑖𝑛𝜃
𝑑𝑡 𝑛 𝑑𝑡 𝑑𝑡 𝑛𝑐𝑜𝑠𝛽 𝑛2 − 𝑠𝑖𝑛2 𝜃 3/2
𝑐𝑜𝑠𝜃 𝑐𝑜𝑠𝜃 Negative sign indicates that the angular acceleration of the rod
𝜔𝑐 = 𝜔 =𝜔
1 𝑛2 − 𝑠𝑖𝑛2 θ tends to reduce the angle β.
𝑛 𝑛 𝑛2 − 𝑠𝑖𝑛2 θ
In the present case angular acceleration of the rod is clockwise.
Theory of Machines - II: Dr. A. B. Andhare 6
Engine Force Analysis
𝐴𝑛𝑔𝑢𝑙𝑎𝑟 𝑎𝑐𝑐𝑒𝑙𝑒𝑟𝑎𝑡𝑖𝑜𝑛 𝛼 𝑔𝑎
mgmg
cosθ
sinθ mg sinθ
cosθ = ≈ = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝐴𝑛𝑔𝑢𝑙𝑎𝑟 𝑑𝑖𝑠𝑝𝑙𝑎𝑐𝑒𝑚𝑒𝑛𝑡 𝜃 (𝑘 2 +𝑎2 )
1 𝑔𝑎
Motion is approximately SHM and hence frequency of oscillation, f = ( 2 2)
2𝜋 𝑘 +𝑎
1 𝑔
For simple pendulum f =
2𝜋 𝐿
𝑘 2 + 𝑎2 𝑘 2
Equivalent length of simple pendulum which gives the same frequency, 𝐿= = +𝑎
𝑎 𝑎
b a
B G D A
b d
mb md
m = Total mass of the connecting rod
l = Length of the connecting rod
Rod replaced by an equivalent massless link with two masses.
One mass is placed at the small end B and other mass is at D
mb + m d = m 𝑏+𝑑
mb * b = md * d 𝑚𝑏 =𝑚 mb = mass at B, md = mass at D
𝑑
𝑏
𝑚𝑏 + 𝑚 𝑏 =𝑚 𝑑 𝑏
𝑑 𝑚𝑏 = 𝑚 Similarly, 𝑚𝑑 = 𝑚
𝑏+𝑑 𝑏+𝑑
𝑏
𝑚𝑏 1 + =𝑚
𝑑 Theory of Machines - II: Dr. A. B. Andhare 18
Inertia of the connecting rod
l
B G
A
b a
B G D A
b d
mb md
𝑑 𝑏
𝑚𝑏 = 𝑚 𝑚𝑑 = 𝑚
𝑏+𝑑 𝑏+𝑑
𝑑 𝑏 𝑏+𝑑
Also, I = mb b2 + md d2 𝐼=𝑚 𝑏2 + 𝑚 𝑑2 𝐼 = 𝑚𝑏𝑑 = 𝑚𝑏𝑑
𝑏+𝑑 𝑏+𝑑 𝑏+𝑑
Let k = radius of gyration of the connecting rod about an axis through G perpendicular to the plane of motion
mk2 = m b d or k2 = b d
𝑘2
Equivalent length of simple pendulum is 𝐿 = +𝑏 =𝑑+𝑏
𝑏
Theory of Machines - II: Dr. A. B. Andhare 19
Inertia of the connecting rod
l
B G
A
B a
b
b B G D A
b d
G mb md
𝑑 𝑏
𝑚𝑏 = 𝑚 𝑚𝑑 = 𝑚
𝑏+𝑑 𝑏+𝑑
𝑘2
𝐼 = 𝑚𝑏𝑑 𝐿= +𝑏 =𝑑+𝑏
𝑏
b = distance of the point of suspension from the center of mass of the body
K = radius of gyration .
d+b = L is the equivalent length if the rod is suspended from the point B and D is the center of oscillation or percussion
b a
B G D A
b d
mb md
𝑑 𝑏
𝑚𝑏 = 𝑚 𝑚𝑑 = 𝑚
𝑏+𝑑 𝑏+𝑑
For connecting rod, it is convenient to assume two point
masses at the center of the two bearings end A and B .
ma + mb = m 𝐼′ = 𝑚𝑎𝑏 Assuming a > d, 𝐼′ > 𝐼
𝑏 𝑏
𝑚𝑎 = 𝑚 =𝑚 Means – By considering two masses at A & B, instead of D & B, the inertia torque is
𝑎+𝑏 𝑙
increased from the actual value (T = I αc). This error is compensated by introducing a
𝑎 𝑎 correction couple.
𝑚𝑏 = 𝑚 =𝑚
𝑎+𝑏 𝑙
Theory of Machines - II: Dr. A. B. Andhare 21
Inertia of the connecting rod
l
B G
A
b a
ma
B G D A
b d
mb md
L
𝐼′ = 𝑚𝑎𝑏 Assuming a > d, 𝐼′ > 𝐼 • Correction couple must be applied in the opposite direction to that
Correction couple of applied inertia torque.
• Direction of applied inertia torque is always opposite to the
∆𝑇 = 𝛼𝑐 𝑚𝑎𝑏 − 𝑚𝑑𝑏 = 𝑚𝑏𝛼𝑐 (𝑎 − 𝑑) direction of angular acceleration.
∆𝑇 = 𝑚𝑏𝛼𝑐 [ 𝑎 + 𝑏 − 𝑏 + 𝑑 ] • The direction of correction couple will be the same as the direction
of angular acceleration i.e. in the direction of decreasing angle β.
∆𝑇 = 𝑚𝑏𝛼𝑐 (𝑙 − 𝐿) Taking (b+ d) = L
∆𝑇
𝑇𝑐 = 𝐹𝑦 𝑟 𝑐𝑜𝑠𝜃 = × 𝑟𝑐𝑜𝑠𝜃
𝑙𝑐𝑜𝑠𝛽
∆𝑇 𝑐𝑜𝑠𝜃 ∆𝑇 𝑐𝑜𝑠𝜃
𝑇𝑐 = 𝑇𝑐 =
𝑙 𝑐𝑜𝑠𝛽 1 𝑛2 − 𝑠𝑖𝑛2 θ
𝑟 𝑛 𝑛
∆𝑇𝑐𝑜𝑠𝜃
𝑇𝑐 = Correction torque is to be deducted from the inertia torque acting on the crankshaft,
𝑛2 − 𝑠𝑖𝑛2 θ
Theory of Machines - II: Dr. A. B. Andhare 23
Inertia of the connecting rod ∆𝑇𝑐𝑜𝑠𝜃
𝑇𝑐 =
l 𝑛2 − 𝑠𝑖𝑛2 θ
L A Due to weight of the mass at A, a torque is exerted on the
ΔT D
r crankshaft, given as
B
β G θ O Ta = (ma g ) r cosθ
For vertical engines, a torque is also exerted on the
Fy Fy crankshaft due to weight of mass at B . It is given as,
l cosβ r cosθ 𝑠𝑖𝑛2𝜃
𝑇 = (𝑚𝑏 𝑔)𝑟 𝑠𝑖𝑛𝜃 +
2 𝑛2 − 𝑠𝑖𝑛2 𝜃
The net torque on the crankshaft will be the algebraic sum of,
• Turning moment due to the force of gas pressure (T)
• Inertia torque due to inertia force at the piston as a result of inertia of the reciprocating mass + the mass of part of
the connecting rod (Tb)
• Inertia torque due to weight of mass at the crank pin which the part of mass of the connecting rod taken at the crank
pin (Ta)
• Inertia torque due to the correction couple (Tc)
• Torque due to weight of the piston in case of vertical engines
La = 0.7181 m
Lb= 0.8505 m
The distance of center of mass of the connecting rod from the center of the small end bearing,
b’ = 654 – (80/2) = 614 mm
Let the second mass be at D, assume GD = d and md = mass at D
Fn
• F = net piston effort, Ft = crank effort
• Turning moment diagram: Plot of T vs. θ
• Inertia effect of connecting rod may be neglected
while plotting T-θ diagram
• Plot of Ft vs. θ (crank effort diagram) is same as T-θ
diagram
𝐴𝑟𝑒𝑎(𝑎𝑏𝑐𝑑𝑒)
𝑇𝑚𝑒𝑎𝑛 =
2𝜋
σ σ
Theory of Machines - II: Dr. A. B. Andhare 38
σ = Circumferential stress
Total centrifugal force per unit length Fc
t
r
dθ
θ
2𝜋 × 250
𝑃 = 𝑇𝜔 = 1000 × = 26180 𝑊 𝑜𝑟 26.18 𝑘𝑊
60
ii) At any instant, ΔT = T – Tmean = (1000 + 300 sin 2θ – 500 cos 2θ) -1000
= 300 sin 2θ – 500 cos 2θ
119.5 119.5
𝑒𝑚𝑎𝑥 = ∆𝑇𝑑𝜃 = (300𝑠𝑖𝑛2𝜃 − 500𝑐𝑜𝑠2𝜃)𝑑𝜃 = −150𝑐𝑜𝑠2𝜃 − 250𝑠𝑖𝑛2𝜃 119.5
29.5 29.5 = 583.1 N-m
29.5
𝑒𝑚𝑎𝑥 = 583.1
𝐾= 2 = 0.01329 𝑜𝑟 1.329 %
𝑚𝑘 2 𝜔 2 2𝜋 × 250
400 × (0.4)2 ×
60
Theory of Machines - II: Dr. A. B. Andhare 41
iii) Acceleration or retardation is produced by excess or less torque than the mean torque at any instant.
ΔT = 300 sin 2θ – 500 cos 2θ
When θ = 60°, ΔT = 259.8 – (-250) = 509.8 N-m
Or Iα = mk2 = 509.8
Or 400 x (0.4)2 x α = 509.8 ; hence α = 9.966 rad/s
30 Nm
60 Nm
T
T 45 Nm
Cylinder1 Cylinder2 Cylinder3
60 Nm T-θ diagram
For each cylinder 30 Nm
Resultant
30 Nm
30 Nm
A single cylinder four stroke petrol engine develops 20 kW at 600 rpm. The mechanical efficiency of
the engine is 80 % and the work done by gases during expansion stroke is three times the work
consumed in the compression stroke. If the total fluctuation of speed is not to exceed ±1.5 % of the
mean speed and turning moment diagram during the compression & expansion stroke is assumed
to be triangular.
Determine the mass of the flywheel if the diameter of it is 1 m. Also, determine the cross section of
the flywheel rim. The width to thickness ratio is 1.5 and mass density of the flywheel material is
7100 kg/m3.
Permissible speed variation = ±1.5 %, Therefore, Coefficient of fluctuation of speed, K = 0.03
2𝜋𝑁 2𝜋 × 600
𝜔= = = 62.83 𝑟𝑎𝑑/𝑠
60 60
𝑏𝑝 20
𝑖𝑝 = = = 25 𝑘𝑊 Indicated power
η𝑚𝑒𝑐ℎ 0.8
As θ is more than 90°, acceleration is negative or towards the left and thus inertia force is
towards the right.
Fn
Fn
Fn
The correction torque is to be deducted from the inertia torque on the crankshaft OR as the force Fy due
to ΔT (which is clockwise) is towards right on the crankshaft , the correction torque is clockwise.
Torque due to weight of the mass at A
Ta = (mg g) r cosθ
= 32.83 x 9.81 x 0.11 x cos 140°= -27.14 N-m clockwise
Hence, total inertia torque on the crankshaft
Tb - Tc + Ta
= -248 – 45.07 – 27.14 = 320.2 clockwise