initial commit / add all books

8 files changed, 214 insertions, 0 deletions

diff --git a/3532/CH7/EX7.1/Ex7_1.sce b/3532/CH7/EX7.1/Ex7_1.sce
new file mode 100644
index 000000000..8ad1c9349
--- /dev/null
+++ b/3532/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,23 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.2.1\n')

+//given data

+E=1.96*10^11//youngs modulus in N/m^2

+I=4*10^-7//moment of area in m^4

+M1=100;M2=50//mass of discs 1 and 2 in Kgs

+c=0.18//distance of disc 1 from support in m

+l=0.3//distance of disc 2 from support in m

+g=9.81//aceleration due to gravity in m/sec^2

+//calculations

+a=[(c^3/(3*E*I)),(c^2/(6*E*I)*(3*l-c));(c^2/(6*E*I)*(3*l-c)),(l^3/(3*E*I))]//from SOM

+y1=g*(M1*a(1,1)+M2*a(1,2))

+y2=g*(M1*a(2,1)+M2*a(2,2))

+Wn=sqrt(g*(M1*y1+M2*y2)/(M1*y1^2+M2*y2^2))

+//now to find out lower natural frequency

+F1=M1*y1*Wn^2

+F2=M2*y2*Wn^2

+y1new=F1*a(1,1)+F2*a(1,2)

+y2new=F1*a(2,1)+F2*a(2,2)

+Wnnew=sqrt((F1*y1new+F2*y2new)/(M1*y1new^2+M2*y2new^2))//actual natural frequency in rad/sec

+//output

+mprintf(' The practical natural frequency Wn is %4.4f rad/sec,but the lower \n natural frequency Wn`is %4.4f rad/sec which is closer to the actual\n natural frequency',Wn,Wnnew)

diff --git a/3532/CH7/EX7.10/Ex7_10.sce b/3532/CH7/EX7.10/Ex7_10.sce
new file mode 100644
index 000000000..1d7ff7f30
--- /dev/null
+++ b/3532/CH7/EX7.10/Ex7_10.sce
@@ -0,0 +1,37 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.7.2\n')

+//given data

+J(1)=100//moment of inertia of first rotor in Kg-m^2

+J(2)=50//moment of inertia of second rotor in Kg-m^2

+J(3)=10//moment of inertia of third rotor in Kg-m^2

+J(4)=50//moment of inertia of fourth rotor in Kg-m^2

+Kt(1)=10^4//stiffness of shaft between 1 and 2 in N-m/rad

+Kt(2)=10^4//stiffness of shaft between 2 and 3 in N-m/rad

+Kt(3)=2*10^4//stiffness of shaft between 3 and 4 in N-m/rad

+To=10000//amplitude of applied torque in N-m

+W=5//frequency of applied torque in rad/sec

+//calculations

+b(1)=-(0.789*To)/3825//twist of shaft 1 in rad

+P(1)=J(1)*W^2

+Q(1)=P(1)*b(1)//twisting moment of shaft 1 in N-m

+R(1)=Q(1)

+S(1)=R(1)/Kt(1)//twist of shaft 1 in radians

+b(2)=b(1)-S(1)//twist of shaft 2 in rad

+P(2)=J(2)*W^2

+Q(2)=P(2)*b(2)

+R(2)=Q(1)+Q(2)+To//twisting moment of shaft 2 in N-m

+S(2)=R(2)/Kt(2)//twist of shaft 2 in radians

+b(3)=b(2)-S(2)//twist of shaft 3 in rad

+P(3)=J(3)*W^2

+Q(3)=P(3)*b(3)

+R(3)=Q(2)+Q(3)//twisting moment of shaft 3 in N-m

+S(3)=R(3)/Kt(3)//twist of shaft 3 in radians

+b(4)=b(3)-S(3)//twist of shaft 4 in rad

+P(4)=J(4)*W^2

+Q(4)=P(4)*b(4)

+R(4)=Q(3)+Q(4)//twisting moment of shaft 4 in N-m

+mprintf('The amplitudes of discs are as follows\n disc1=%4.4f rad\n disc2=%4.4f rad\n disc3=%4.4f rad\n disc4=%4.4f rad',b(1),b(2),b(3),b(4))

+mprintf('\nThe twists of shaft are as follows\nfirst shaft=%5.5f rad\nsecond shaft=%5.5f rad\nthird shaft=%5.5f rad',S(1),S(2),S(3))

+mprintf('\nThe twisting moments of shafts are as follows\nfirst shaft=%5.5f N-m\nsecond shaft=%5.5f N-m\nthird shaft=%5.5f N-m',R(1),R(2),R(3)) 

+mprintf('\nNOTE:The slight difference in values are due to the more accurate values\ncalculated by SCILAB')

diff --git a/3532/CH7/EX7.2.1/Ex7_1.sce b/3532/CH7/EX7.2.1/Ex7_1.sce
new file mode 100644
index 000000000..8ad1c9349
--- /dev/null
+++ b/3532/CH7/EX7.2.1/Ex7_1.sce
@@ -0,0 +1,23 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.2.1\n')

+//given data

+E=1.96*10^11//youngs modulus in N/m^2

+I=4*10^-7//moment of area in m^4

+M1=100;M2=50//mass of discs 1 and 2 in Kgs

+c=0.18//distance of disc 1 from support in m

+l=0.3//distance of disc 2 from support in m

+g=9.81//aceleration due to gravity in m/sec^2

+//calculations

+a=[(c^3/(3*E*I)),(c^2/(6*E*I)*(3*l-c));(c^2/(6*E*I)*(3*l-c)),(l^3/(3*E*I))]//from SOM

+y1=g*(M1*a(1,1)+M2*a(1,2))

+y2=g*(M1*a(2,1)+M2*a(2,2))

+Wn=sqrt(g*(M1*y1+M2*y2)/(M1*y1^2+M2*y2^2))

+//now to find out lower natural frequency

+F1=M1*y1*Wn^2

+F2=M2*y2*Wn^2

+y1new=F1*a(1,1)+F2*a(1,2)

+y2new=F1*a(2,1)+F2*a(2,2)

+Wnnew=sqrt((F1*y1new+F2*y2new)/(M1*y1new^2+M2*y2new^2))//actual natural frequency in rad/sec

+//output

+mprintf(' The practical natural frequency Wn is %4.4f rad/sec,but the lower \n natural frequency Wn`is %4.4f rad/sec which is closer to the actual\n natural frequency',Wn,Wnnew)

diff --git a/3532/CH7/EX7.3.1/Ex7_4.sce b/3532/CH7/EX7.3.1/Ex7_4.sce
new file mode 100644
index 000000000..2b5f13fc7
--- /dev/null
+++ b/3532/CH7/EX7.3.1/Ex7_4.sce
@@ -0,0 +1,19 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.3.1\n')

+//given data

+E=1.96*10^11//youngs modulus in N/m^2

+I=4*10^-7//moment of area in m^4

+M1=100;M2=50//mass of discs 1 and 2 in Kgs

+c=0.18//distance of disc 1 from support in m

+l=0.3//distance of disc 2 from support in m

+g=9.81//aceleration due to gravity in m/sec^2

+//calculations

+a=[(c^3/(3*E*I)),(c^2/(6*E*I)*(3*l-c));(c^2/(6*E*I)*(3*l-c)),(l^3/(3*E*I))]//from SOM

+y1=g*M1*a(1,1)//considering only M1 to be acting

+y2=g*M2*a(2,2)//considering only M2 to be acting

+W1=sqrt(g/y1)

+W2=sqrt(g/y2)

+Wn=sqrt(1/((1/W1^2)+(1/W2^2)))//applying Eqn 7.3.7,Sec7.3

+//output

+mprintf(' The natural frequency of transverse vibration obtained from \n Dunkerly method is %4.4f rad/sec which is slightly lower\n than the correct value',Wn)

diff --git a/3532/CH7/EX7.4.1/Ex7_5.sce b/3532/CH7/EX7.4.1/Ex7_5.sce
new file mode 100644
index 000000000..c14b212b0
--- /dev/null
+++ b/3532/CH7/EX7.4.1/Ex7_5.sce
@@ -0,0 +1,28 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.4.1\n')

+//given data

+E=1.96*10^11//youngs modulus in N/m^2

+I=4*10^-7//moment of area in m^4

+M1=100;M2=50//mass of discs 1 and 2 in Kgs

+c=0.18//distance of disc 1 from support in m

+l=0.3//distance of disc 2 from support in m

+g=9.81//aceleration due to gravity in m/sec^2

+//calculations

+a=[(c^3/(3*E*I)),(c^2/(6*E*I)*(3*l-c));(c^2/(6*E*I)*(3*l-c)),(l^3/(3*E*I))]//from SOM

+x1(1)=1;x2(1)=1

+for i=1:10//upto 10th iteration for more perfect answer

+F1(i)=100*x1(i)//'i' represents the dash(')

+F2(i)=50*x2(i)

+x1(i)=F1(i)*a(1,1)+F2(i)*a(1,2)

+x2(i)=F1(i)*a(2,1)+F2(i)*a(2,2)

+r=(x2(i)/x1(i))

+x2(i+1)=r

+x1(i+1)=1

+end

+x1dd=1

+W1=(x1dd/x1(10))

+W2=(r/x2(10))

+Wn=sqrt((W1+W2)/2)//natural frequency in rad/sec

+mprintf('The natural frequency of system in iilustrative example 7.2.1 obtained by\nStodala method is Wn=%f rad/sec',Wn)

+mprintf('\nNOTE:The obtained answer is more near to the perfect answer \since 10 iterations/trials\nhas been carried out.In textbook only upto 3rd iteration has been carried out')

diff --git a/3532/CH7/EX7.4/Ex7_4.sce b/3532/CH7/EX7.4/Ex7_4.sce
new file mode 100644
index 000000000..2b5f13fc7
--- /dev/null
+++ b/3532/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,19 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.3.1\n')

+//given data

+E=1.96*10^11//youngs modulus in N/m^2

+I=4*10^-7//moment of area in m^4

+M1=100;M2=50//mass of discs 1 and 2 in Kgs

+c=0.18//distance of disc 1 from support in m

+l=0.3//distance of disc 2 from support in m

+g=9.81//aceleration due to gravity in m/sec^2

+//calculations

+a=[(c^3/(3*E*I)),(c^2/(6*E*I)*(3*l-c));(c^2/(6*E*I)*(3*l-c)),(l^3/(3*E*I))]//from SOM

+y1=g*M1*a(1,1)//considering only M1 to be acting

+y2=g*M2*a(2,2)//considering only M2 to be acting

+W1=sqrt(g/y1)

+W2=sqrt(g/y2)

+Wn=sqrt(1/((1/W1^2)+(1/W2^2)))//applying Eqn 7.3.7,Sec7.3

+//output

+mprintf(' The natural frequency of transverse vibration obtained from \n Dunkerly method is %4.4f rad/sec which is slightly lower\n than the correct value',Wn)

diff --git a/3532/CH7/EX7.5/Ex7_5.sce b/3532/CH7/EX7.5/Ex7_5.sce
new file mode 100644
index 000000000..c14b212b0
--- /dev/null
+++ b/3532/CH7/EX7.5/Ex7_5.sce
@@ -0,0 +1,28 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.4.1\n')

+//given data

+E=1.96*10^11//youngs modulus in N/m^2

+I=4*10^-7//moment of area in m^4

+M1=100;M2=50//mass of discs 1 and 2 in Kgs

+c=0.18//distance of disc 1 from support in m

+l=0.3//distance of disc 2 from support in m

+g=9.81//aceleration due to gravity in m/sec^2

+//calculations

+a=[(c^3/(3*E*I)),(c^2/(6*E*I)*(3*l-c));(c^2/(6*E*I)*(3*l-c)),(l^3/(3*E*I))]//from SOM

+x1(1)=1;x2(1)=1

+for i=1:10//upto 10th iteration for more perfect answer

+F1(i)=100*x1(i)//'i' represents the dash(')

+F2(i)=50*x2(i)

+x1(i)=F1(i)*a(1,1)+F2(i)*a(1,2)

+x2(i)=F1(i)*a(2,1)+F2(i)*a(2,2)

+r=(x2(i)/x1(i))

+x2(i+1)=r

+x1(i+1)=1

+end

+x1dd=1

+W1=(x1dd/x1(10))

+W2=(r/x2(10))

+Wn=sqrt((W1+W2)/2)//natural frequency in rad/sec

+mprintf('The natural frequency of system in iilustrative example 7.2.1 obtained by\nStodala method is Wn=%f rad/sec',Wn)

+mprintf('\nNOTE:The obtained answer is more near to the perfect answer \since 10 iterations/trials\nhas been carried out.In textbook only upto 3rd iteration has been carried out')

diff --git a/3532/CH7/EX7.7.2/Ex7_10.sce b/3532/CH7/EX7.7.2/Ex7_10.sce
new file mode 100644
index 000000000..1d7ff7f30
--- /dev/null
+++ b/3532/CH7/EX7.7.2/Ex7_10.sce
@@ -0,0 +1,37 @@
+clc

+clear

+mprintf('Mechanical vibrations by G.K.Grover\n Example 7.7.2\n')

+//given data

+J(1)=100//moment of inertia of first rotor in Kg-m^2

+J(2)=50//moment of inertia of second rotor in Kg-m^2

+J(3)=10//moment of inertia of third rotor in Kg-m^2

+J(4)=50//moment of inertia of fourth rotor in Kg-m^2

+Kt(1)=10^4//stiffness of shaft between 1 and 2 in N-m/rad

+Kt(2)=10^4//stiffness of shaft between 2 and 3 in N-m/rad

+Kt(3)=2*10^4//stiffness of shaft between 3 and 4 in N-m/rad

+To=10000//amplitude of applied torque in N-m

+W=5//frequency of applied torque in rad/sec

+//calculations

+b(1)=-(0.789*To)/3825//twist of shaft 1 in rad

+P(1)=J(1)*W^2

+Q(1)=P(1)*b(1)//twisting moment of shaft 1 in N-m

+R(1)=Q(1)

+S(1)=R(1)/Kt(1)//twist of shaft 1 in radians

+b(2)=b(1)-S(1)//twist of shaft 2 in rad

+P(2)=J(2)*W^2

+Q(2)=P(2)*b(2)

+R(2)=Q(1)+Q(2)+To//twisting moment of shaft 2 in N-m

+S(2)=R(2)/Kt(2)//twist of shaft 2 in radians

+b(3)=b(2)-S(2)//twist of shaft 3 in rad

+P(3)=J(3)*W^2

+Q(3)=P(3)*b(3)

+R(3)=Q(2)+Q(3)//twisting moment of shaft 3 in N-m

+S(3)=R(3)/Kt(3)//twist of shaft 3 in radians

+b(4)=b(3)-S(3)//twist of shaft 4 in rad

+P(4)=J(4)*W^2

+Q(4)=P(4)*b(4)

+R(4)=Q(3)+Q(4)//twisting moment of shaft 4 in N-m

+mprintf('The amplitudes of discs are as follows\n disc1=%4.4f rad\n disc2=%4.4f rad\n disc3=%4.4f rad\n disc4=%4.4f rad',b(1),b(2),b(3),b(4))

+mprintf('\nThe twists of shaft are as follows\nfirst shaft=%5.5f rad\nsecond shaft=%5.5f rad\nthird shaft=%5.5f rad',S(1),S(2),S(3))

+mprintf('\nThe twisting moments of shafts are as follows\nfirst shaft=%5.5f N-m\nsecond shaft=%5.5f N-m\nthird shaft=%5.5f N-m',R(1),R(2),R(3)) 

+mprintf('\nNOTE:The slight difference in values are due to the more accurate values\ncalculated by SCILAB')