In [2]:

```
import math
# Variables:
d = 300./1000
r = d/2
l = 600./1000 #m
m = 5. #kg
N = 300. #rpm
#Solution:
#Calculating the angular speed of the disc
omega = 2*math.pi*N/60 #rad/s
#Calculating the mass moment of inertia of the disc
#about an axis through its centre of gravity and perpendicular to the plane of the disc
I = m*r**2/2 #kg-m**2
#Calculating the couple due to mass of disc
C = m*9.81*l #N-m
#Calculating the speed of precession
omegaP = C/(I*omega) #rad/s
#Results:
print " Speed of precession, omegaP = %.1f rad/s."%(omegaP)
```

In [7]:

```
import math
# Variables:
d = 150./1000
r = d/2
x = 100./1000 #m
m = 5. #kg
N = 1000.
NP = 60. #rpm
#Solution:
#Calculating the angular speed of the disc
omega = round(2*math.pi*N/60,1) #rad/s
#Calculating the speed of precession of the axle
omegaP = round(2*math.pi*NP/60,3) #rad/s
#Calculating the mass moment of inertia of the disc
# about an axis through its centre of gravity and perpendicular to the plane of disc
I = round(m*r**2/2,3) #kg-m**2
#Calculating the gyroscopic couple acting on the disc
C = round(I*omega*omegaP,1) #N-m
#Calculating the force at each bearing due to the gyroscopic couple
F = C/x #N
#Calculating the reactions at the bearings A and B
RA = m/2*9.81 #N
RB = round(RA,1) #N
#Resulmath.tant reaction at each bearing:
#Calculating the resultant reaction at the bearing A
RA1 = F+RA #N
#Calculating the resultant reaction at the bearing B
RB1 = F-RB #N
#Results:
print " Resultant reaction at the bearing A RA1 = %.1f N upwards."%(RA1)
print " Resultant reaction at the bearing B RB1 = %.1f N downwards."%( RB1)
```

In [5]:

```
import math
# Variables:
R = 50.
k = 0.3 #m
v = 200.*1000/3600 #m/s
m = 400. #kg
N = 2400. #rpm
#Solution:
#Calculating the angular speed of the engine
omega = 2*math.pi*N/60 #rad/s
#Calculating the mass moment of inertia of the engine and the propeller
I = m*k**2 #kg-m**2
#Calculating the angular velocity of precession
omegaP = v/R #rad/s
#Calculating the gyroscopic couple acting on the aircraft
C = I*omega*omegaP/1000 #kN-m
#Results:
print " Gyroscopic couple acting on the aircraft, C = %.3f kN-m."%(C)
print " The effect of the gyroscopic couple is to lift the nose upwards and tail downwards."
```

In [11]:

```
import math
# Variables:
m = 8.*1000 #kg
k = 0.6
R = 75. #m
N = 1800. #rpm
v = 100.*1000/3600 #m/s
#Solution:
#Calculating the angular speed of the rotor
omega = round(2*math.pi*N/60,1) #rad/s
#Calculating the mass moment of inertia of the rotor
I = m*k**2 #kg-m**2
#Calculating the angular velocity of precession
omegaP = round(v/R,2) #rad/s
#Calculating the gyroscopic couple
C = I*omega*omegaP/1000 #kN-m
#Results:
print " Gyroscopic couple, C = %.3f kN-m."%(C)
```

In [7]:

```
import math
# Variables:
N = 1500. #rpm
m = 750. #kg
omegaP = 1. #rad/s
k = 250./1000 #m
#Solution:
#Calculating the angular speed of the rotor
omega = 2*math.pi*N/60 #rad/s
#Calculating the mass moment of inertia of the rotor
I = m*k**2 #kg-m**2
#Calculating the gyroscopic couple transmitted to the hull
C = I*omega*omegaP/1000 #kN-m
#Results:
print " Gyroscopic couple transmitted to the hull, C = %.3f kN-m."%(C)
print " When the pitching is upward, the relative gyroscopic couple acts in the clockwise direction."
```

In [8]:

```
import math
# Variables:
m = 3500. #kg
k = 0.45 #m
N = 3000. #rpm
#Solution:
#Calculating the angular speed of the rotor
omega = 2*math.pi*N/60 #rad/s
#When the ship is steering to the left:
R = 100. #m
v = 36.*1000/3600 #m/s
#Calculating the mass moment of inertia of the rotor
I = m*k**2 #kg-m**2
#Calculating the angular velocity of precession
omegaP = v/R #rad/s
#Calculating the gyroscopic couple
C = I*omega*omegaP/1000 #kN-m
#Results:
print " Gyroscopic couple when the ship is steering to the left. C = %.2f kN-m."%(C)
print " When the rotor rotates clockwise and the ship takes a left turn.\
the effect of the reactive gyroscopic couple is to raise the bow and lower the stern."
#When the ship is pitching with the bow falling:
tp = 40. #s
#Calculating the amplitude of swing
phi = 12./2*math.pi/180 #rad
#Calculating the angular velocity of the simple harmonic motion
omega1 = 2*math.pi/tp #rad/s
#Calculating the maximum angular velocity of precession
omegaP = phi*omega1 #rad/s
#Calculating the gyroscopic couple
C = I*omega*omegaP/1000 #kN-m
#Results:
print " Gyroscopic couple when the ship is pitching with the bow falling, C = %.3f kN-m."%(C)
print " When the bow is falling, the effect of the reactive gyroscopic couple is to move\
the ship towards port side."
```

In [9]:

```
import math
# Variables:
m = 20.*1000 #kg
k = 0.6 #m
N = 2000. #rpm
phi = 6.*math.pi/180 #rad
tp = 30. #s
#Solution:
#Calculating the angular speed of the rotor
omega = 2*math.pi*N/60 #rad/s
#Maximum gyroscopic couple:
#Calculating the mass moment of inertia of the rotor
I = m*k**2 #kg-m**2
#Calculating the angular velocity of the simple harmonic motion
omega1 = 2*math.pi/tp #rad/s
#Calculating the maximum angular velocity of precession
omegaPmax = phi*omega1 #rad/s
#Calculating the maximum gyroscopic couple
Cmax = I*omega*omegaPmax/1000 #kN-m
#Calculating the maximum angular acceleration during pitching
alphamax = phi*omega1**2 #Maximum angular acceleration during pitching rad/s**2
#Results:
print " Maximum gyroscopic couple, Cmax = %.3f kN-m."%(Cmax)
print " Maximum angular acceleration during pitching = %.4f rad/s**2."%(alphamax)
print " When the rotation of the rotor is clockwise when looking from the left and when\
the bow is rising, then the reactive gyroscopic couple acts in the direction which\
tends to turn the bow towrds right."
```

In [30]:

```
import math
# Variables:
m = 5.*1000 #kg
N = 1000. #rpm
k = 0.5 #m
#Solution:
#Calculating the angular speed of the rotor
omega = 2.*math.pi*2100/60 #rad/s
#When the ship steers to the left:
v = 30.*1000/3600 #m/s
R = 60. #m
#Calculating the angular velocity of precession
omegaP = round(v/R,2) #rad/s
#Calculating the mass moment of inertia of the rotor
I = m*k**2 #kg-m**2
#Calculating the gyroscopic couple
C = I*omega*omegaP/1000 #kN-m
#Results:
print " C = %.1f k N-m"%C
#When the ship pitches with the bow descending:
phi = round(6*math.pi/180,3) #rad
tp = 20. #s
#Calculating the angular velocity of simple harmonic motion
omega1 = round(2*math.pi/tp,4) #rad/s
#Calculating the maximum velocity of precession
omegaPmax = round(phi*omega1,3) #rad/s
#Calculating the maximum gyroscopic couple
Cmax = I*omega*omegaPmax #N-m
#Results:
print " Cmax = %.f N-m"%Cmax
#When the ship rolls:
omegaP = 0.03 #rad/s
#Calculating the gyroscopic couple
C = I*omega*omegaP #N-m
#Results:
print " C = %.2f N-m"%(round(C,-1))
#Calculating the maximum angular acceleration during pitching
alphamax = phi*omega1**2 #rad/s**2
#Results:
print " Maximum angular acceleration during pitching, alphamax = %.2f rad/s**2."%(alphamax)
# rounding off error. please check.
```

In [6]:

```
import math
# Variables:
m = 2000. #kg
N = 3000. #rpm
k = 0.5
R = 100. #m
v = 16.1*1855/3600 #m/s
#Solution:
#Calculating the angular speed of the rotor
omega = 2*math.pi*N/60 #rad/s
#Gyroscopic couple:
#Calculating the mass moment of inertia of the rotor
I = m*k**2 #kg-m**2
#Calculating the angular velocity of precession
omegaP = v/R #rad/s
#Calculating the gyroscopic couple
C = I*omega*omegaP/1000 #kN-m
#Torque during pitching:
tp = 50. #s
phi = 12./2*math.pi/180 #rad
#Calculating the angular velocity of simple harmonic motion
omega1 = 2*math.pi/tp #rad/s
#Calculating the maximum angular velocity of precession
omegaPmax = phi*omega1 #rad/s
#Calculating the maximum gyroscopic couple during pitching
Cmax = I*omega*omegaPmax #N-m
#Calculating the maximum acceleration during pitching
alphamax = phi*omega1**2 #rad/s**2
#Results:
print " Torque during pitching, Cmax = %d N-m."%(Cmax)
print " Maximum acceleration during pitching, alphamax = %.5f rad/s**2."%(alphamax)
```

In [41]:

```
import math
# Variables:
m = 2500. #kg
x = 1.5
R = 30.
dW = 0.75
rW = dW/2
h = 0.9 #m
v = 24.*1000/3600 #m/s
G = 5.
IW = 18.
IE = 12. #kg-m**2
#Solution:
#Calculating the road reaction on each wheel
r = m*9.81/4 #Road reaction on each wheel N
#Calculating the angular velocity o the wheels
omegaW = round(v/rW,1) #rad/s
#Calculating the angular velocity of precession
omegaP = round(v/R,2) #rad/s
#Calculating the gyroscopic couple due to one pair of wheels and axle
CW = round(2*IW*omegaW*omegaP) #N-m
#Calculating the gyroscopic couple due to the rotating parts of the motor and gears
CE = round(2*IE*G*omegaW*omegaP) #N-m
#Calculating the net gyroscopic couple
C = CW-CE #N-m
#Calculating the reaction due to gyroscopic couple at each of the outer or inner wheels
P = round((-C)/(2*x),1) #N
#Calculating the centrifugal force
FC = round(m*v**2/R,1) #N
#Calculating the overturning couple
CO = FC*h #N-m
#Calculating the reaction due to overturning couple at each of the outer and inner wheels
Q = 2*CO/(2*x) #N
#Calculating the vertical force exerted on each outer wheel
PO = m*9.81/4-P/2+Q/2 #N
#Calculating the vertical force exerted on each inner wheel
PI = m*9.81/4+P/2-Q/2 #N
#Results:
print " Vertical force exerted on each outer wheel, PO = %.2f N."%(PO)
print " Vertical force exerted on each inner wheel, PI = %.2f N."%(PI)
# note : value of Fc is calculated wrongly. please check using calculator. so answers are different.
```

In [12]:

```
import math
# variables
R = 100 #m
IW = 2.5 # kg-m 2
dW = 0.6 #m
rW = 0.3 #m
IE = 1.2 #kg-m 2 ;
G = 3.
m = 1600. # kg
h = 0.5 #m ;
x = 1.5 #m
v = 1.
# calculations
road_reaction = m*9.81/4 #N
wW = v / rW
wP = v/R
cW = 4 * IW*wW*wP
cE = IE*G*wW*wP
C = cW + cE
Pby2 = C/(2*x)
Fc = m*v**2/R
Co = Fc * h
Qby2 = Co/(2*x)
v_2 = road_reaction/(Pby2 + Qby2)
v = math.sqrt(v_2)
ans = v*3600./1000
# result
print "V <= %.1f m/s = %.1f * 3600 /1000 = %.2f km/h"%(v,v,ans)
```

In [25]:

```
import math
# Variables:
m = 2000. #kg
mE = 75. #kg
b = 2.5 #m
x = 1.5 #m
h = 500./1000 #m
L = 1. #m
dW = 0.8 #m
rW = round(dW/2,1) #m
kE = 100./1000 #m
R = 60. #m
IW = 0.8 #kg-m**2
G = 4.
v = round(60.*1000/3600,2) #m/s
#Solution:
#Refer Fig. 14.12
#Calculating the weight on the rear wheels
W2 = (m*9.81*1)/b #N
#Calculating the weight on the front wheels
W1 = m*9.81-W2 #N
#Calculating the weight on each of the front wheels
Wf = W1/2 #Weight on each of the front wheels N
#Calculating the weight on each of the rear wheels
Wr = W2/2 #Weight on each of the rear wheels N
#Calculating the angular velocity of wheels
omegaW = v/rW #rad/s
#Calculating the angular velocity of precession
omegaP = round(v/R,3) #rad/s
#Calculating the gyroscopic couple due to four wheels
CW = round(4*IW*omegaW*omegaP,1) #N-m
#Calculating the magnitude of reaction due to gyroscopic couple due to four wheels at each of the inner or outer wheel
P = round((CW/(2*x)),2) #N
#Calculating the mass moment of inertia of rotating parts of the engine
IE = mE*(kE)**2 #kg-m**2
#Calculating the gyroscopic couple due to rotating parts of the engine
CE = round(IE*(kE)**2*G*omegaW*omegaP*100,1) #N-m
#Calculating the magnitude of reaction due to gyroscopic couple due to rotating parts of the engine at each of the inner or outer wheel
F = (CE/(2*b)) #N
#Calculating the centrifugal force
FC = round(m*v**2/R) #N
#Calculating the centrifugal couple tending to overturn the car
CO = FC*h #N-m
#Calculating the magnitude of reaction due to overturning couple at each of the inner or outer wheel
Q = round((CO/(2*x)),2) #N
#Calculating the load on front wheel 1
Fw1 = (W1/2)-(P/2)-(F/2)-(Q/2) #Load on front wheel 1 N
#Calculating the load on front wheel 2
Fw2 = W1/2+P/2-F/2+Q/2 #Load on front wheel 2 N
#Calculating the load on rear wheel 3
Rw3 = W2/2-P/2+F/2-Q/2 #Load on rear wheel 3 N
#Calculating the load on rear wheel 4
Rw4 = W2/2+P/2+F/2+Q/2 #Load on rear wheel 4 N
#Results:
print " Load on front wheel 1 = %.2f N."%(Fw1)
print " Load on front wheel 2 = %.2f N."%(Fw2)
print " Load on rear wheel 3 = %.2f N."%(Rw3)
print " Load on rear wheel 4 = %.2f N."%(Rw4)
# note : incorrect answers in the textbook
```

In [3]:

```
import math
# Variables:
m = 2000. #kg
mI = 200. #kg
x = 1.6 #m
R = 30. #m
dW = 0.7 #m
rW = dW/2 #m
k = 0.3 #m
h = 1. #m
v = 54.*1000/3600 #m/s
theta = 8. #degrees
#Solution:
#Refer Fig. 14.13
#Calculating the reactions at the wheels:
#Taking moments about B
RA = (m*9.81*math.cos(math.radians(theta))+m*v**2/R*math.sin(math.radians(theta)))*1/2+(m*9.81*math.sin(math.radians(theta)) \
-m*v**2/R*math.cos(math.radians(theta)))*h/x #N
#Resolving the forces perpendicular to the track
RB = (m*9.81*math.cos(math.radians(theta))+m*v**2/R*math.sin(math.radians(theta)))-RA #N
#Calculating the angular velocity of wheels
omegaW = v/rW #rad/s
#Calculating the angular velocity of precession
omegaP = v/R #rad/s
#Calculating the gyroscopic couple
C = mI*k**2*omegaW*math.cos(math.radians(theta))*omegaP #N-m
#Calculating the force at each pair of wheels due to the gyroscopic couple
P = C/x #N
#Calculating the pressure on the inner rail
PI = RA-P #N
#Calculating the pressure o the outer rail
PO = RB+P #N
#Results:
print " Pressure on the inner rail, PI = %.2f N."%(PI)
print " Pressure on the outer rail, PO = %.2f N."%(PO)
```

In [17]:

```
import math
# Variables:
I = 180. #kg-m**2
D = 1.8 #m
R = D/2 #m
x = 1.5 #m
v = 95.*1000/3600 #m/s
t = 0.1 #s
#Solution:
#Gyroscopic couple set up:
#Calculating the angular velocity of the locomotive
omega = v/R #rad/s
#Calculating the amplitude
A = 1./2*6 #mm
#Calculating the maximum velocity while falling
vmax = 2*math.pi/t*A/1000 #m/s
#Calculating the maximum angular velocity of tilt of the axle or angular velocity of precession
omegaPmax = vmax/x #rad/s
#Calculating the gyroscopic couple set up
C = I*omega*omegaPmax #N-m
#Calculating the reaction between the wheel and rail due to the gyroscopic couple
P = C/x #N
#Results:
print " Gyroscopic couple set up, C = %.1f N-m."%(C)
print " Reaction between the wheel and rail due to the gyroscopic couple, P = %d N."%(P)
```

In [62]:

```
import math
# Variables:
m = 250. #kg
IE = 0.3 #kg-m**2
IW = 1. #kg-m**2
G = 5.
h = 0.6 #m
rW = 300./1000 #m
R = 50. #m
v = 90.*1000/3600 #m/s
#Solution:
#Calculating the angle of inclination with respect to the vertical of a two wheeler
#Equating total overturning couple to balancing couple
tantheta = math.atan(1/(m*9.81*h))*((v**2/(R*rW)*(2*IW+G*IE))+(m*v**2/R*h)) #degrees
theta = math.degrees(math.atan(tantheta))
#Results:
print " Angle of inclination with respect to the vertical of a two wheeler, tan theta = %.2f ."%(theta)
```

In [6]:

```
import math
# Variables:
m1 = 0.5 #kg
m2 = 0.3 #kg
k = 20./1000 #m
OG = 10./1000 #m
h = OG #m
R = 50. #m
N = 3000. #rpm
v = 15. #m/s
#Solution:
#Refer Fig. 14.15 and Fig. 14.16
#Calculating the angular speed of the wheel
omega = 2*math.pi*N/60 #rad/s
#Calculating the mass moment of inertia of the gyrowheel
I = m1*k**2 #kg-m**2
#Calculating the angular velocity of precession
omegaP = v/R #rad/s
#When the vehicle moves in the direction of arrow X taking a left turn along the curve:
#Calculating the angle of inclination of the gyrowheel from the vertical
#Equating the overturning couple to the balancing couple for equilibrium condition
tantheta1 = (1/(m2*9.81*h))*(I*omega*omegaP-m2*v**2/R*h) #degrees
theta1 = math.degrees(math.atan(tantheta1))
#When the vehicle reverses at the same speed in the direction of arrow Y along the same path:
#Calculating the angle of inclination of the gyrowheel from the vertical
#Equating the overturning couple to the balancing couple for equilibrium condition
tantheta2 = round((1/(m2*9.81*h))*(I*omega*omegaP+m2*v**2/R*h),1) #degrees
theta2 = math.degrees(math.atan(tantheta2))
#Results:
print " Angle of inclination of the gyrowheel from the vertical when the vehicle moves\
in the direction of arrow X taking a left turn along the curve, theta = %.2f degrees."%(theta1)
print " Angle of inclination of the gyrowheel from the vertical when the vehicle reverses \
at the same speed in the direction of arrow Y along the same path, theta = %.f degrees."%(theta2)
# rounding error
```

In [25]:

```
import math
# Variables:
d = 0.6 #m
r = d/2 #m
m = 30. #kg
theta = 1. #degree
N = 1200. #rpm
#Solution:
#Calculating the angular speed of the shaft
omega = 2*math.pi*N/60 #rad/s
#Calculating the gyroscopic couple acting on the bearings
C = round(m/8*omega**2*r**2*math.sin(math.radians(2*theta))) #N-m
#Results:
print " Gyroscopic couple acting on the bearings, C = %d N-m."%(C)
```