38 files changed, 1066 insertions, 0 deletions
diff --git a/698/CH14/EX14.1/1_derivation_of_torque_uniform_pressure.txt b/698/CH14/EX14.1/1_derivation_of_torque_uniform_pressure.txt
new file mode 100644
index 000000000..36af99bd3
--- /dev/null
+++ b/698/CH14/EX14.1/1_derivation_of_torque_uniform_pressure.txt
@@ -0,0 +1,19 @@
+Consider a differential area dA=2*pi*r dr.
+The differential normal force=dN=p dA=p(2*pi*r dR),
+The differential force dQ=f dN=f(p(2*pi*r dR)).
+The differential frictional torque=dT=r dQ=r(f*p(2*pi*r dR));
+
+Intergrating with respect to r, with p and f as constants, over r=Ri to r=Ro,
+we get the total torque as
+		T=2*pi*f*p[((Ro^3)-(Ri^3))/3]
+
+The axial force
+		F=p*pi*((Ro^2)-(Ri^2))
+from which the average pressure,
+		p=F/[pi*((Ro^2)-(Ri^2))]
+
+Substituting this value of p into
+		T=2*pi*f*p[((Ro^3)-(Ri^3))/3]
+We obtain
+
+		T=F*f*[(2/3)*((Ro^3)-(Ri^3))] = F*f*Rf
+\ No newline at end of file
diff --git a/698/CH14/EX14.1/P1_derivation_of_torque_uniform_pressure.sce b/698/CH14/EX14.1/P1_derivation_of_torque_uniform_pressure.sce
new file mode 100644
index 000000000..601fff9be
--- /dev/null
+++ b/698/CH14/EX14.1/P1_derivation_of_torque_uniform_pressure.sce
@@ -0,0 +1,22 @@
+clc
+//Example 14.1
+//Derivation of torque capacity
+
+//------------------------------------------------------------------------------
+//This example is derivation based, hence the code will comprise only of statements printed to text file
+//Printing result file to .txt
+res1=mopen(TMPDIR+'1_derivation_of_torque_uniform_pressure.txt','wt')
+mfprintf(res1,"Consider a differential area dA=2*pi*r dr.\nThe differential normal force=dN=p dA=p(2*pi*r dR),\n")
+mfprintf(res1,"The differential force dQ=f dN=f(p(2*pi*r dR)).\nThe differential frictional torque=dT=r dQ=r(f*p(2*pi*r dR));\n\n")
+mfprintf(res1,"Intergrating with respect to r, with p and f as constants, over r=Ri to r=Ro,\nwe get the total torque as\n")
+mfprintf(res1,"\t\tT=2*pi*f*p[((Ro^3)-(Ri^3))/3]\n\n")
+mfprintf(res1,"The axial force\n")
+mfprintf(res1,"\t\tF=p*pi*((Ro^2)-(Ri^2))\n")
+mfprintf(res1,"from which the average pressure,\n")
+mfprintf(res1,"\t\tp=F/[pi*((Ro^2)-(Ri^2))]\n\n")
+mfprintf(res1,"Substituting this value of p into\n")
+mfprintf(res1,"\t\tT=2*pi*f*p[((Ro^3)-(Ri^3))/3]\n")
+mfprintf(res1,"We obtain\n\n")
+mfprintf(res1,"\t\tT=F*f*[(2/3)*((Ro^3)-(Ri^3))] = F*f*Rf")
+mclose(res1)
+editor(TMPDIR+'1_derivation_of_torque_uniform_pressure.txt')
+\ No newline at end of file
diff --git a/698/CH14/EX14.10/10_determination_of_axial_thrust_and_pressure_intensity.txt b/698/CH14/EX14.10/10_determination_of_axial_thrust_and_pressure_intensity.txt
new file mode 100644
index 000000000..e2963865a
--- /dev/null
+++ b/698/CH14/EX14.10/10_determination_of_axial_thrust_and_pressure_intensity.txt
@@ -0,0 +1,7 @@
+(a)Axial thrust required to transmit the power is 848.83 N
+(b)The pressure equation is:
+	p=F/(2*pi*(Ro-Ri)*r)
+
+(c)	Maximum contact pressure occurs at inner radius, and is equal to 0.108 MPa
+   	Minimum contact pressure occurs at outer radius, and is equal to 0.072 MPa
+   	Average contact pressure is 0.086 MPa
diff --git a/698/CH14/EX14.10/P10_Determination_of_axial_thrust_and_pressure_intensity.sce b/698/CH14/EX14.10/P10_Determination_of_axial_thrust_and_pressure_intensity.sce
new file mode 100644
index 000000000..ed55773e3
--- /dev/null
+++ b/698/CH14/EX14.10/P10_Determination_of_axial_thrust_and_pressure_intensity.sce
@@ -0,0 +1,55 @@
+clc
+//Example 14.10
+//Determination of axial thrust and pressure intensity
+
+//------------------------------------------------------------------------------
+//Given Data:
+//Power to be transmitted
+P=10000 //Watt
+//Speed
+N=1000 //rpm
+//Outer and Inner diameters
+Do=0.15 //m
+Di=0.1 //m
+Ro=0.15/2 //m
+Ri=0.1/2 //m
+// number of surfaces
+n=6
+//coefficient of friction
+f=0.3
+//------------------------------------------------------------------------------
+
+// Using uniform wear theory
+// Mean radius
+Rm=(Ro+Ri)/2
+
+// Torque(T)=power/angular velocity
+T=(P*60)/(2*%pi*N)
+
+// T=F*f*Rm*n
+// Axial thrust F
+F=T/(f*Rm*n)
+
+//contact pressure at radius r:
+//p=F/(2*%pi*(Ro-Ri)*r)
+//maximum contact pressure(pmax) is at inner radius
+pmax=F/(2*%pi*(Ro-Ri)*Ri)
+//minimum contact pressure(pmin) is at outer radius
+pmin=F/(2*%pi*(Ro-Ri)*Ro)
+// average contact pressure
+pavg=F/(2*%pi*(Ro-Ri)*Rm)
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res10=mopen(TMPDIR+'10_determination_of_axial_thrust_and_pressure_intensity.txt','wt')
+mfprintf(res10,"(a)Axial thrust required to transmit the power is %0.2f N\n",F)
+mfprintf(res10,"(b)The pressure equation is:\n\tp=F/(2*pi*(Ro-Ri)*r)\n\n")
+mfprintf(res10,"(c)\tMaximum contact pressure occurs at inner radius, and is equal to %0.3f MPa\n",pmax*(10^-6))
+mfprintf(res10,"   \tMinimum contact pressure occurs at outer radius, and is equal to %0.3f MPa\n",pmin*(10^-6))
+mfprintf(res10,"   \tAverage contact pressure is %0.3f MPa\n",pavg*(10^-6))
+mclose(res10)
+if (isdef('editor') | (funptr('editor')<>0)) then
+  editor(TMPDIR+'10_determination_of_axial_thrust_and_pressure_intensity.txt')
+end
+//------------------------------------------------------------------------------
+//---------------------------------End of program-------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.11/11_outer_diameter_of_clutch.txt b/698/CH14/EX14.11/11_outer_diameter_of_clutch.txt
new file mode 100644
index 000000000..ec9dab3a8
--- /dev/null
+++ b/698/CH14/EX14.11/11_outer_diameter_of_clutch.txt
@@ -0,0 +1,3 @@
+Necessary outer diameter of disks is 70.00 mm
+Necessary axial force is 666.00 N
+Actual contact pressure is 353.32 kN/m^2
+\ No newline at end of file
diff --git a/698/CH14/EX14.11/P11_outer_diameter_of_clutch.sce b/698/CH14/EX14.11/P11_outer_diameter_of_clutch.sce
new file mode 100644
index 000000000..359e3c208
--- /dev/null
+++ b/698/CH14/EX14.11/P11_outer_diameter_of_clutch.sce
@@ -0,0 +1,53 @@
+clc
+//Example 14.11
+//Outer diameter of clutch
+
+//------------------------------------------------------------------------------
+//Given Data:
+// number of plates
+n1=5
+n2=4
+// if n is total number of surfces
+n=n1+n2-1
+
+//Total torqe transmitting capacity
+Tt=16// Nm
+//Permissible inner diameter
+Di=0.05// m
+//coefficient of friction
+f=0.1
+//average pressure
+p=350000 //N/(m^2)
+
+//------------------------------------------------------------------------------
+//Torque per pair of surfaces
+T=Tt/n
+
+//T=F*f*((Do+Di)/4)
+//T=((%pi/4)*((Do^2)-(Di^2))*p)*f*((Do+Di)/4)
+//To solve above equation for Do, it has to brought to a polynomial equation form in Do
+//((%pi*p*f)*(Do^3))+((%pi*p*f*Di)*(Do^2))-((%pi*p*f*(Di^2))*Do)-((%pi*p*f*(Di^3))+(16*T))=0
+x=poly([-((%pi*p*f*(Di^3))+(16*T)) -(%pi*p*f*(Di^2)) (%pi*p*f*Di) (%pi*p*f)],'Do','c')
+y=roots(x)
+//y will contain all roots of the polynomial, the first of which is the acceptable one
+Do=y(1)
+
+//Axial force F
+F=T/(f*((Do+Di)/4))
+
+Do=round(Do*(10^3))
+F=round(F)
+//Actual pressure
+p=F/((%pi/4)*(((Do*(10^-3))^2)-(Di^2)))
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res11=mopen(TMPDIR+'11_outer_diameter_of_clutch.txt','wt')
+mfprintf(res11,"Necessary outer diameter of disks is %0.2f mm\n",Do)
+mfprintf(res11,"Necessary axial force is %0.2f N\n",F)
+mfprintf(res11,"Actual contact pressure is %0.2f kN/m^2",p*(10^-3))
+mclose(res11)
+editor(TMPDIR+'11_outer_diameter_of_clutch.txt')
+
+//------------------------------------------------------------------------------
+//-----------------------------End of program-----------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.12/12_conditions_appropriate_for_uniform_pressure.txt b/698/CH14/EX14.12/12_conditions_appropriate_for_uniform_pressure.txt
new file mode 100644
index 000000000..fd80f2ad1
--- /dev/null
+++ b/698/CH14/EX14.12/12_conditions_appropriate_for_uniform_pressure.txt
@@ -0,0 +1,2 @@
+Uniform pressure assumption is more appropriate where the plates are
+flexible to permit deflection when wear occurs.
diff --git a/698/CH14/EX14.12/P12_Conditions_appropriate_for_uniform_pressure.sce b/698/CH14/EX14.12/P12_Conditions_appropriate_for_uniform_pressure.sce
new file mode 100644
index 000000000..fc6092a0f
--- /dev/null
+++ b/698/CH14/EX14.12/P12_Conditions_appropriate_for_uniform_pressure.sce
@@ -0,0 +1,13 @@
+clc
+//Example 14.12
+//Conditions preferring uniform pressure
+
+//------------------------------------------------------------------------------
+//This example is derivation based, hence the code will comprise only of statements printed to text file
+//Printing result file to .txt
+res12=mopen(TMPDIR+'12_conditions_appropriate_for_uniform_pressure.txt','wt')
+mfprintf(res12,"Uniform pressure assumption is more appropriate where the plates are\nflexible to permit deflection when wear occurs.\n")
+mclose(res12)
+editor(TMPDIR+'12_conditions_appropriate_for_uniform_pressure.txt')
+//------------------------------------------------------------------------------
+//-------------------------------End of program---------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.14/14_derivation_of_torque_cone_clutch.txt b/698/CH14/EX14.14/14_derivation_of_torque_cone_clutch.txt
new file mode 100644
index 000000000..7cdfdfaa3
--- /dev/null
+++ b/698/CH14/EX14.14/14_derivation_of_torque_cone_clutch.txt
@@ -0,0 +1,36 @@
+Consider a differential element bounded by circles of radii r and (r+dr).
+The area of the differential frustrum of a cone is
+		dA=2*pi*r[dr/sin (alpha)]
+
+Integrating the differential torque with respect to r over r=Ri to r=Ro,
+we get
+		T=[(2*pi*p*f)/sin (alpha)]*[((Ro^3)-(Ri^3))/3]
+
+Define the force Fn as that due to the pressure applied
+to the area as if it were stretched out into a plane:
+		F=p(2*pi*Rm*b)
+
+To relate Fn to the axial force F,consider a differential element
+with central angle d(phi).
+The differential area is
+		dA=2*pi*Rm*b*(d(phi)/(2*pi))
+		 =Rm*b*d(phi)
+The differential normal force is
+	dN=p*Rm*b*d(phi)*sin(alpha)
+The horizontal component of the differential force is dF; then
+Integrating with respect to (phi) over (phi)=0 to (phi)=2*pi
+we get
+		F=2*pi*p*Rm*b*sin(alpha)
+		 =Fn*sin(alpha)
+
+Substituting equation of p in equation of T, we get
+		T=[(Fn*f)/(Rm*b*sin (alpha))]*[((Ro^3)-(Ri^3))/3]
+
+		 =(Fn*f)*[(2/3){((Ro^3)-(Ri^3))/((Ro^2)-(Ri^2))}]
+
+		 =[(F*f)/sin (alpha)]*[(2/3){((Ro^3)-(Ri^3))/((Ro^2)-(Ri^2))}]
+
+since
+Rm=(Ro+Ri)/2
+(b*sin (alpha))=Ro-Ri
+and Fn=F/sin (alpha)
+\ No newline at end of file
diff --git a/698/CH14/EX14.14/P14_Derivation_of_torque_cone_clutch.sce b/698/CH14/EX14.14/P14_Derivation_of_torque_cone_clutch.sce
new file mode 100644
index 000000000..7db96f9f3
--- /dev/null
+++ b/698/CH14/EX14.14/P14_Derivation_of_torque_cone_clutch.sce
@@ -0,0 +1,28 @@
+clc
+//Example 14.14 
+//Derivation of torque capacity of a cone clutch
+
+//------------------------------------------------------------------------------
+//This example is derivation based, hence the code will comprise only of statements printed to text file
+//Printing result file to .txt
+res14=mopen(TMPDIR+'14_derivation_of_torque_cone_clutch.txt','wt')
+mfprintf(res14,"Consider a differential element bounded by circles of radii r and (r+dr).\nThe area of the differential frustrum of a cone is\n")
+mfprintf(res14,"\t\tdA=2*pi*r[dr/sin (alpha)]\n\n")
+mfprintf(res14,"Integrating the differential torque with respect to r over r=Ri to r=Ro,\nwe get\n")
+mfprintf(res14,"\t\tT=[(2*pi*p*f)/sin (alpha)]*[((Ro^3)-(Ri^3))/3]\n\n")
+mfprintf(res14,"Define the force Fn as that due to the pressure applied\nto the area as if it were stretched out into a plane:\n")
+mfprintf(res14,"\t\tF=p(2*pi*Rm*b)\n\n")
+mfprintf(res14,"To relate Fn to the axial force F,consider a differential element\nwith central angle d(phi).\nThe differential area is\n")
+mfprintf(res14,"\t\tdA=2*pi*Rm*b*(d(phi)/(2*pi))\n\t\t =Rm*b*d(phi)\n")
+mfprintf(res14,"The differential normal force is\n\tdN=p*Rm*b*d(phi)*sin(alpha)\nThe horizontal component of the differential force is dF; then\n")
+mfprintf(res14,"Integrating with respect to (phi) over (phi)=0 to (phi)=2*pi\nwe get\n")
+mfprintf(res14,"\t\tF=2*pi*p*Rm*b*sin(alpha)\n\t\t =Fn*sin(alpha)\n\n")
+mfprintf(res14,"Substituting equation of p in equation of T, we get\n")
+mfprintf(res14,"\t\tT=[(Fn*f)/(Rm*b*sin (alpha))]*[((Ro^3)-(Ri^3))/3]\n\n")
+mfprintf(res14,"\t\t =(Fn*f)*[(2/3){((Ro^3)-(Ri^3))/((Ro^2)-(Ri^2))}]\n\n")
+mfprintf(res14,"\t\t =[(F*f)/sin (alpha)]*[(2/3){((Ro^3)-(Ri^3))/((Ro^2)-(Ri^2))}]\n\n")
+mfprintf(res14,"since\nRm=(Ro+Ri)/2\n(b*sin (alpha))=Ro-Ri\nand Fn=F/sin (alpha)")
+mclose(res14)
+editor(TMPDIR+'14_derivation_of_torque_cone_clutch.txt')
+//------------------------------------------------------------------------------
+//-------------------------------End of program---------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.16/16_parametrs_of_cone_clutch.txt b/698/CH14/EX14.16/16_parametrs_of_cone_clutch.txt
new file mode 100644
index 000000000..e36487c37
--- /dev/null
+++ b/698/CH14/EX14.16/16_parametrs_of_cone_clutch.txt
@@ -0,0 +1,4 @@
+(a)Axial force F required to transmit the torque is 636.6 N
+(b)Axial force required to engage the clutch, enagagement taking place when clutch is not rotating, is 1799.3 N
+(c)Average normal pressure when maximum torque is being transmitted is 74.1 kN/m^2
+(d)Maximum normal pressure assuming uniform wear is 75.8 kN/m^2
diff --git a/698/CH14/EX14.16/P16_parametrs_of_cone_clutch.sce b/698/CH14/EX14.16/P16_parametrs_of_cone_clutch.sce
new file mode 100644
index 000000000..7fe18dab4
--- /dev/null
+++ b/698/CH14/EX14.16/P16_parametrs_of_cone_clutch.sce
@@ -0,0 +1,58 @@
+clc
+//Example 14.16
+//Parameterts of cone clutch
+//------------------------------------------------------------------------------
+
+//Torque
+T=200//Nm
+//Speed
+N=1250//rev/min
+//Large diameter
+Do=0.35//m
+Ro=Do/2
+//Face width
+b=0.065//m
+//coefficient of friction
+f=0.2
+//cone pitch angle
+alpha=6.25//degrees
+//converting alpha in degrees to radians
+alpha=(alpha*%pi)/180
+
+//------------------------------------------------------------------------------
+
+//mean radius
+Rm=Ro-((1/2)*b*sin (alpha))
+Rm=floor(Rm*(10^3))
+Rm=Rm*(10^-3)
+//Rm=(Ro+Ri)/2
+//smaller radius
+Ri=(2*Rm)-Ro
+Ri=floor(Ri*(10^3))
+Ri=Ri*(10^-3)
+
+//T=(F*f*Rm)/(sin alpha)
+//Axial force F required to transmit the torque
+F=(T*sin (alpha))/(f*Rm)
+
+//Axial force required to engage the clutch when clutch is not rotating
+//Fe=Fn((sin alpha)+(f*cos alpha))
+Fe=(T/(f*Rm))*((sin (alpha))+(f*cos (alpha)))
+
+//average normal pressure when maximum torque is being transmitted
+p=F/(%pi*((Ro^2)-(Ri^2)))
+
+//maximum normal pressure assuming uniform wear
+Pmax=F/(2*%pi*(Ro-Ri)*Ri)
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res16=mopen(TMPDIR+"16_parametrs_of_cone_clutch.txt",'wt')
+mfprintf(res16,"(a)Axial force F required to transmit the torque is %0.1f N\n",F)
+mfprintf(res16,"(b)Axial force required to engage the clutch, enagagement taking place when clutch is not rotating, is %0.1f N\n",Fe)
+mfprintf(res16,"(c)Average normal pressure when maximum torque is being transmitted is %0.1f kN/m^2\n",p*(10^-3))
+mfprintf(res16,"(d)Maximum normal pressure assuming uniform wear is %0.1f kN/m^2\n",Pmax*(10^-3))
+mclose(res16)
+editor(TMPDIR+"16_parametrs_of_cone_clutch.txt")
+//------------------------------------------------------------------------------
+//---------------------------------End of program-------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.17/17_cone_clutch_under_uniform_pressure.txt b/698/CH14/EX14.17/17_cone_clutch_under_uniform_pressure.txt
new file mode 100644
index 000000000..f534b086d
--- /dev/null
+++ b/698/CH14/EX14.17/17_cone_clutch_under_uniform_pressure.txt
@@ -0,0 +1,3 @@
+(a)Axial force F required to transmit the torque is 563.04 N
+(b)Axial force required to engage the clutch, enagagement taking place when clutch is not rotating, is 1591.26 N
+(c)Average normal pressure when maximum torque is being transmitted is 74.06 kN/m^2
diff --git a/698/CH14/EX14.17/P17_cone_clutch_under_uniform_pressure.sce b/698/CH14/EX14.17/P17_cone_clutch_under_uniform_pressure.sce
new file mode 100644
index 000000000..eff120d63
--- /dev/null
+++ b/698/CH14/EX14.17/P17_cone_clutch_under_uniform_pressure.sce
@@ -0,0 +1,57 @@
+clc
+//Exercise 14.17
+//same clutch as previous problem, assuming uniform pressure theory
+
+//------------------------------------------------------------------------------
+//Given Data:
+//Torque
+T=200//Nm
+//Speed
+N=1250//rev/min
+//Large diameter
+Do=0.35//m
+Ro=Do/2
+//Face width
+b=0.065//m
+//coefficient of friction
+f=0.2
+//cone pitch angle
+alpha=6.25//degrees
+//converting alpha in degrees to radians
+alpha=(alpha*%pi)/180
+
+//------------------------------------------------------------------------------
+//mean radius
+Rm=Ro-((1/2)*b*sin (alpha))
+Rm=floor(Rm*(10^3))
+Rm=Rm*(10^-3)
+//Rm=(Ro+Ri)/2
+//smaller radius
+Ri=(2*Rm)-Ro
+Ri=floor(Ri*(10^3))
+Ri=Ri*(10^-3)
+
+//Axial force reqiured to transmit the torque
+//t=F*f*((Ro^3-Ri^3)/(3*Rm*b*(sin (alpha))^2))
+F=T/(f*((Ro^3-Ri^3)/(3*Rm*b*(sin (alpha))^2)))
+
+//Axial force required to engage the clutch when clutch is not rotating
+//Fe=Fn((sin alpha)+(f*cos alpha))
+Fn=F/sin(alpha)
+Fe=Fn*((sin (alpha))+(f*cos (alpha)))
+
+//average normal pressure when maximum torque is being transmitted
+//Fn=p*2*%pi*Rm*b
+p=Fn/(2*%pi*Rm*b)
+
+//------------------------------------------------------------------------------
+//Printing file to .txt
+res17=mopen(TMPDIR+'17_cone_clutch_under_uniform_pressure.txt','wt')
+mfprintf(res17,"(a)Axial force F required to transmit the torque is %0.2f N\n",F)
+mfprintf(res17,"(b)Axial force required to engage the clutch, enagagement taking place when clutch is not rotating, is %0.2f N\n",Fe)
+mfprintf(res17,"(c)Average normal pressure when maximum torque is being transmitted is %0.2f kN/m^2\n",p*(10^-3))
+mclose(res17)
+editor(TMPDIR+'17_cone_clutch_under_uniform_pressure.txt')
+
+//------------------------------------------------------------------------------
+//------------------------------End of program----------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.18/18_design_of_cone_clutch.txt b/698/CH14/EX14.18/18_design_of_cone_clutch.txt
new file mode 100644
index 000000000..22bbfe99f
--- /dev/null
+++ b/698/CH14/EX14.18/18_design_of_cone_clutch.txt
@@ -0,0 +1,18 @@
+(a)Torque transmitting capacity of clutch is 119.366 Nm
+
+(b)To find axial thrust and dimensions, we need to solve the following equations:
+	F=pmax*2*pi*(Ri^2)*(r-1)---Eq 1
+	(F*Ri)=(2*T*sin(alpha))/(f*(1+r))---Eq 2
+
+(c)Width of face is half of mean radius:
+	b=(1/2)*Rm
+	Ro=r*Ri
+	where r=(4+sin (alpha))/(4-sin (alpha))
+
+(d)Solving the above equations, we get
+	MAIN DIMENSIONS:
+Inner radius of clutch is 125.00 mm
+Outer radius of clutch is 140.00 mm
+face width of clutch is 66.25 mm
+
+Axial force required to engage the clutch is 971.930 N
diff --git a/698/CH14/EX14.18/P18_Design_of_cone_clutch.sce b/698/CH14/EX14.18/P18_Design_of_cone_clutch.sce
new file mode 100644
index 000000000..a1f678805
--- /dev/null
+++ b/698/CH14/EX14.18/P18_Design_of_cone_clutch.sce
@@ -0,0 +1,88 @@
+clc
+//Exercise 14.18
+//Design of cone clutch
+
+//------------------------------------------------------------------------------
+//Given Data:
+//Power to be transmitted
+P=10000 //Watt
+//Speed
+N=800 //rpm
+//Limiting normal pressure
+pmax=0.09*(10^6) //Pa
+//coefficient of friction
+f=0.2
+//cone pitch angle
+alpha=24/2 //degrees
+
+//------------------------------------------------------------------------------
+//Function to standardize the dimensions
+std=[1 2 3 4 5 6 8 10 12 15 18 20 24 26 30:5:500]//array of standard dimensions
+n=length(std)
+
+funcprot(0)
+function y=stddim(x)
+    x=x*(10^3)
+    for i=1:n
+        if (x<std(i)) then
+            
+            y=std(i)
+            break
+            
+        else
+            continue
+        end
+    end 
+    y=y*(10^-3)
+endfunction
+
+
+//------------------------------------------------------------------------------
+//To find dimensions of clutch Ri,Ro,Rm,b and axial force required to engage the clutch
+//b=Rm/2
+//Rm=(Ro+Ri)/2
+//In cone clutch, (Ro-Ri)=b*sin alpha
+//Assume r=Ro/Ri
+//Solving the above equations, we obtain a relation
+//r=(4+sin alpha)/(4-sin alpha)
+r=(4+sind (alpha))/(4-sind (alpha))
+
+r=(round(r*100))*(10^-2)//rounding off to two decimal places
+
+//pmax=F/(2*%pi*(Ro-Ri)*Ri)
+
+//Torque transmitting capacity
+T=(P*60)/(2*%pi*N)
+//T=(F*f*Rm)/sin(alpha)
+//F*Rm=(T*sin(alpha))/f
+//F*((Ro+Ri)/2)=(T*sin(alpha))/f
+
+
+Ri=((2*T*sind(alpha))/(f*(1+r)*(r-1)*2*%pi*pmax))^(1/3)
+Ri=stddim(Ri)
+F=pmax*2*%pi*(Ri^2)*(r-1)
+Ro=r*Ri
+Ro=stddim(Ro)
+Rm=(Ro+Ri)/2
+b=Rm/2
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res18=mopen(TMPDIR+'18_design_of_cone_clutch.txt','wt')
+mfprintf(res18,"(a)Torque transmitting capacity of clutch is %0.3f Nm\n\n",T)
+mfprintf(res18,"(b)To find axial thrust and dimensions, we need to solve the following equations:\n")
+mfprintf(res18,"\tF=pmax*2*pi*(Ri^2)*(r-1)---Eq 1\n")
+mfprintf(res18,"\t(F*Ri)=(2*T*sin(alpha))/(f*(1+r))---Eq 2\n\n")
+mfprintf(res18,"(c)Width of face is half of mean radius:\n")
+mfprintf(res18,"\tb=(1/2)*Rm\n")
+mfprintf(res18,"\tRo=r*Ri\n\twhere r=(4+sin (alpha))/(4-sin (alpha))\n\n")
+mfprintf(res18,"(d)Solving the above equations, we get\n")
+mfprintf(res18,"\tMAIN DIMENSIONS:\n")
+mfprintf(res18,"Inner radius of clutch is %0.2f mm\n",Ri*(10^3))
+mfprintf(res18,"Outer radius of clutch is %0.2f mm\n",Ro*(10^3))
+mfprintf(res18,"face width of clutch is %0.2f mm\n\n",b*(10^3))
+mfprintf(res18,"Axial force required to engage the clutch is %0.3f N\n",F)
+mclose(res18)
+editor(TMPDIR+'18_design_of_cone_clutch.txt')
+//------------------------------------------------------------------------------
+//---------------------------------End of program-------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.19/19_comparison_between_disk_clutch_and_cone_clutch.txt b/698/CH14/EX14.19/19_comparison_between_disk_clutch_and_cone_clutch.txt
new file mode 100644
index 000000000..62dc49161
--- /dev/null
+++ b/698/CH14/EX14.19/19_comparison_between_disk_clutch_and_cone_clutch.txt
@@ -0,0 +1,5 @@
+When both clutches operate at same speed,
+subjected to same axial force
+and both have same mean diameter
+for the given dimensions,
+Power capacity of disk clutch = (1.042) times power capacity of cone clutch
+\ No newline at end of file
diff --git a/698/CH14/EX14.19/P19_comparison_between_disk_clutch_and_cone_clutch.sce b/698/CH14/EX14.19/P19_comparison_between_disk_clutch_and_cone_clutch.sce
new file mode 100644
index 000000000..2e659c69a
--- /dev/null
+++ b/698/CH14/EX14.19/P19_comparison_between_disk_clutch_and_cone_clutch.sce
@@ -0,0 +1,38 @@
+clc
+//Example 14.19
+//Comparison between disk clutch and cone clutch
+
+//------------------------------------------------------------------------------
+//Given Data:
+//Mean radius
+Rm=0.1//m
+//Axial force
+F=1000//N
+//coefficient of friction
+f=0.2
+//number of plates in disk clutch
+n1=4
+n2=3
+//total number of surfaces in contact
+n=n1+n2-1
+//angle of cone clutch
+a=20//degrees
+alpha=a/2 //pitch angle
+
+//------------------------------------------------------------------------------
+//Torque carrying capacity of a disk clutch
+Td=(n*F*f*Rm)
+//Torque carrying capacity of cone clutch
+Tc=(F*f*Rm)/sind(alpha)
+//Ratio T(disk) to T(cone)
+r=Td/Tc
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res19=mopen(TMPDIR+'19_comparison_between_disk_clutch_and_cone_clutch.txt','wt')
+mfprintf(res19,"When both clutches operate at same speed,\nsubjected to same axial force\nand both have same mean diameter\nfor the given dimensions,\n")
+mfprintf(res19,"Power capacity of disk clutch = (%0.3f) times power capacity of cone clutch",r)
+mclose(res19)
+editor(TMPDIR+'19_comparison_between_disk_clutch_and_cone_clutch.txt')
+//------------------------------------------------------------------------------
+//-----------------------------------End of program----------------------------- 
+\ No newline at end of file
diff --git a/698/CH14/EX14.2/2_derivation_of_torque_uniform_wear.txt b/698/CH14/EX14.2/2_derivation_of_torque_uniform_wear.txt
new file mode 100644
index 000000000..937d9feb2
--- /dev/null
+++ b/698/CH14/EX14.2/2_derivation_of_torque_uniform_wear.txt
@@ -0,0 +1,32 @@
+When a clutch is new, the pressure may be rather uniform.
+If the surfaces are relatively rigid, the outer portion, where velocity is high,
+will wear more than inner portion.
+After the initial wearing-in,it is reasonable to assume that the curve of the profile will maintain its shape;
+or, the wear thereafter may be considered to be uniform.
+Uniform wear can be expressed in a different way by saying that
+at any time interval, the work done per unit area is constant:
+[(frictional force)*(velocity)]/area = 
+	[(f*p*2*pi*r dr)*(rw)]/(2*pi*r dr) = constant C
+or
+	p=C'/f*r*w	 Since f and w are constant,
+	p=C/r, where C is constant.
+An alternate method of showing that pressure varies inversely as the radius is
+to consider that wear (delta) is proportional to pressure p and velocity V.
+Thus 
+	(delta)=K*p*V=K*p*(r*w),
+or
+	p=C/r		since (delta) and K are constants and w is fixed for a given clutch.
+The differential frictional torque=dT= r(f*p(2*pi*r dR));
+
+Intergrating with respect to r over r=Ri to r=Ro,
+we get the total torque as
+		T=2*pi*f*C[((Ro^2)-(Ri^2))/2]
+
+To find C, we can integrate p(2*pi*r dr) with respect to r over r=Ri to r=Ro
+We get
+		C=F/(2*pi*(Ro-Ri))
+
+Substituting this value of C into T
+We obtain
+
+		T=F*f*[(1/2)*(Ro+Ri)] = F*f*Rf
+\ No newline at end of file
diff --git a/698/CH14/EX14.2/P2_Derivation_of_torque_uniform_wear.sce b/698/CH14/EX14.2/P2_Derivation_of_torque_uniform_wear.sce
new file mode 100644
index 000000000..2d6ecd036
--- /dev/null
+++ b/698/CH14/EX14.2/P2_Derivation_of_torque_uniform_wear.sce
@@ -0,0 +1,28 @@
+clc
+//Example 14.2
+//Derivation of torque capacity for one pair of surfaces subjected to uniform wear
+
+//------------------------------------------------------------------------------
+//This example is derivation based, hence the code will comprise only of statements printed to text file
+//Printing result file to .txt
+res2=mopen(TMPDIR+'2_derivation_of_torque_uniform_wear.txt','wt')
+mfprintf(res2,"When a clutch is new, the pressure may be rather uniform.\nIf the surfaces are relatively rigid, the outer portion, where velocity is high,\nwill wear more than inner portion.\n")
+mfprintf(res2,"After the initial wearing-in,it is reasonable to assume that the curve of the profile will maintain its shape;\nor, the wear thereafter may be considered to be uniform.\n")
+mfprintf(res2,"Uniform wear can be expressed in a different way by saying that\nat any time interval, the work done per unit area is constant:\n")
+mfprintf(res2,"[(frictional force)*(velocity)]/area = \n\t[(f*p*2*pi*r dr)*(rw)]/(2*pi*r dr) = constant C\n")
+mfprintf(res2,"or\n\tp=C''/f*r*w\t Since f and w are constant,")
+mfprintf(res2,"\n\tp=C/r, where C is constant.\n")
+mfprintf(res2,"An alternate method of showing that pressure varies inversely as the radius is\nto consider that wear (delta) is proportional to pressure p and velocity V.\n")
+mfprintf(res2,"Thus \n\t(delta)=K*p*V=K*p*(r*w),\nor\n\tp=C/r\t\tsince (delta) and K are constants and w is fixed for a given clutch.\n")
+mfprintf(res2,"The differential frictional torque=dT= r(f*p(2*pi*r dR));\n\n")
+mfprintf(res2,"Intergrating with respect to r over r=Ri to r=Ro,\nwe get the total torque as\n")
+mfprintf(res2,"\t\tT=2*pi*f*C[((Ro^2)-(Ri^2))/2]\n\n")
+mfprintf(res2,"To find C, we can integrate p(2*pi*r dr) with respect to r over r=Ri to r=Ro\n")
+mfprintf(res2,"We get\n\t\tC=F/(2*pi*(Ro-Ri))\n\n")
+mfprintf(res2,"Substituting this value of C into T\n")
+mfprintf(res2,"We obtain\n\n")
+mfprintf(res2,"\t\tT=F*f*[(1/2)*(Ro+Ri)] = F*f*Rf")
+mclose(res2)
+editor(TMPDIR+'2_derivation_of_torque_uniform_wear.txt')
+//------------------------------------------------------------------------------
+//-------------------------------End of program---------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.20/20_face_width_of_cone_clutch.txt b/698/CH14/EX14.20/20_face_width_of_cone_clutch.txt
new file mode 100644
index 000000000..b05fb4128
--- /dev/null
+++ b/698/CH14/EX14.20/20_face_width_of_cone_clutch.txt
@@ -0,0 +1 @@
+The required face width is 14.0 mm
+\ No newline at end of file
diff --git a/698/CH14/EX14.20/P20_face_width_of_cone_clutch.sce b/698/CH14/EX14.20/P20_face_width_of_cone_clutch.sce
new file mode 100644
index 000000000..f683211c7
--- /dev/null
+++ b/698/CH14/EX14.20/P20_face_width_of_cone_clutch.sce
@@ -0,0 +1,68 @@
+clc
+//Example 14.20
+//Face width of cone clutch
+
+//------------------------------------------------------------------------------
+//Given Data:
+//Torque to be transmitted
+T=100//Nm
+//Available axial force
+Fe=850//N
+//Total included cone angle
+A=24//degrees
+//Pitch angle
+alpha=A/2
+//converting alpha in degrees to radians
+alpha=(alpha*%pi)/180
+//maximum pressure
+p=100*(10^3)//N/m^2
+//coefficient of friction
+f=0.20
+
+//------------------------------------------------------------------------------
+//Function to standardize the dimensions
+std=[1 2 3 4 5 6 8 10 12 14 15 16 18 20 24 26 30:5:300]
+n=length(std)
+
+funcprot(0)
+function y=stddim(x)
+    x=x*(10^3)
+    for i=1:n
+        if (x<std(i)) then
+            y=std(i)
+            break
+            
+        else
+            continue
+        end
+    end 
+    y=y*(10^-3)
+endfunction
+//------------------------------------------------------------------------------
+//Assuming uniform wear
+//To find required face width 'b'
+
+//Force required to engage the clutch
+//Fe=Fn*((sin (alpha))+(f*cos (alpha)))
+Fn=Fe/((sin (alpha))+(f*cos (alpha)))
+
+//Torque transmitting capacity
+//T=Fn*f*Rm
+//Mean radius
+Rm=T/(Fn*f)
+Rm=stddim(Rm)
+
+//Normal force
+//Fn=p*2*%pi*Rm*b
+//Face width
+b=Fn/(p*2*%pi*Rm)
+b=stddim(b)
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res20=mopen(TMPDIR+'20_face_width_of_cone_clutch.txt','wt')
+mfprintf(res20,"The required face width is %0.1f mm",b*(10^3))
+mclose(res20)
+editor(TMPDIR+'20_face_width_of_cone_clutch.txt')
+//------------------------------------------------------------------------------
+//-------------------------------End of program------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.21/21_design_of_cone_clutch_given_mean_diameter.txt b/698/CH14/EX14.21/21_design_of_cone_clutch_given_mean_diameter.txt
new file mode 100644
index 000000000..8bb6d74d0
--- /dev/null
+++ b/698/CH14/EX14.21/21_design_of_cone_clutch_given_mean_diameter.txt
@@ -0,0 +1,4 @@
+Axial force required to prevent slipping is 1176.31 N
+Face width required is 315.00 mm
+Outer radius is 110.5 mm
+Inner radius is 89.5 mm
+\ No newline at end of file
diff --git a/698/CH14/EX14.21/P21_design_of_cone_clutch_given_mean_diameter.sce b/698/CH14/EX14.21/P21_design_of_cone_clutch_given_mean_diameter.sce
new file mode 100644
index 000000000..b0871cef5
--- /dev/null
+++ b/698/CH14/EX14.21/P21_design_of_cone_clutch_given_mean_diameter.sce
@@ -0,0 +1,76 @@
+clc
+//Example 14.21
+//Design of cone clutch with given mean diameter
+
+//------------------------------------------------------------------------------
+//Given Data:
+//power to be transmitted
+P=20000//Watts
+//speed
+N=600//rpm
+//Mean diameter
+Dm=0.2//m
+//Mean radius
+Rm=Dm/2
+//Taper of cone
+a=1/15
+alpha=atand(a)
+//coefficient of friction
+f=0.18
+
+//------------------------------------------------------------------------------
+//Function to standardize the dimensions
+funcprot(0)
+function y=stddim(x)
+    x=x*(10^3)
+    std=[1 2 3 4 5 6 8 10 12 15 18 20 24 26 30:5:400]
+    n=length(std)
+    for i=1:n
+        if (x<std(i)) then
+            y=std(i)
+            break
+        else
+            continue
+        end
+    end 
+    y=y*(10^-3)
+endfunction
+
+//------------------------------------------------------------------------------
+//Torque transmitted
+T=(P*60)/(2*%pi*N)
+
+//T=(F*f*Rm)/sin(alpha)
+//Axial force transmitted
+F=(T*sind(alpha))/(f*Rm)
+
+//To design the cone clutch
+//Assuming allowable normal pressure
+p=0.09*(10^6)//Pa
+//Fn=F=p*2*%pi*Rm*b
+//face width b
+b=F/(p*2*%pi*Rm*sind(alpha))
+b=stddim(b)
+
+//Ri=inner radius 
+//Ro=outer radius 
+//(Ro+Ri)/2=Rm
+//(Ro-Ri)/b=sin (alpha)
+//solving the two equations, storing coefficients of Ri and Ro in matrix R, and constants in matrix c
+R=[1 1;1 -1]
+c=[(2*Rm);(b*sind(alpha))]
+x=R\c
+Ro=x(1)
+Ri=x(2)
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res21=mopen(TMPDIR+"21_design_of_cone_clutch_given_mean_diameter.txt",'wt')
+mfprintf(res21,"Axial force required to prevent slipping is %0.2f N\n",F)
+mfprintf(res21,"Face width required is %0.2f mm\n",b*(10^3))
+mfprintf(res21,"Outer radius is %0.1f mm\n",Ro*(10^3))
+mfprintf(res21,"Inner radius is %0.1f mm",Ri*(10^3))
+mclose(res21)
+editor(TMPDIR+"21_design_of_cone_clutch_given_mean_diameter.txt")
+//------------------------------------------------------------------------------
+//-----------------------------End of program-----------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.3/3_friction_radii_for_different_conditions.txt b/698/CH14/EX14.3/3_friction_radii_for_different_conditions.txt
new file mode 100644
index 000000000..0976bba38
--- /dev/null
+++ b/698/CH14/EX14.3/3_friction_radii_for_different_conditions.txt
@@ -0,0 +1,8 @@
+CASE I: Ro=100mm and Ri=90mm
+The friction radius based on uniform pressure theory is 95.1 mm
+The friction radius based on uniform wear theory is 95.0 mm
+
+
+CASE II: Ro=100mm ; Ri= 25mm
+The friction radius based on uniform pressure theory is 70.0 mm
+The friction radius based on uniform wear theory is 62.5 mm
+\ No newline at end of file
diff --git a/698/CH14/EX14.3/P3_Friction_radii_for_different_conditions.sce b/698/CH14/EX14.3/P3_Friction_radii_for_different_conditions.sce
new file mode 100644
index 000000000..656cb6abd
--- /dev/null
+++ b/698/CH14/EX14.3/P3_Friction_radii_for_different_conditions.sce
@@ -0,0 +1,47 @@
+clc
+//Example 14.3
+//Comparison of friction radius between clutches of different dimensions
+//Ro= outer radius; Ri=inner radius
+
+//-----------------------------------------------------------------------------
+//Case 1: Ro=100mm ; Ri= 90mm
+Ro1=100//mm
+Ri1=90//mm
+Ro31=Ro1^3
+Ri31=Ri1^3
+Ro21=Ro1^2
+Ri21=Ri1^2
+// with uniform pressure:
+Rfp1=((2/3)*(Ro31-Ri31))/(Ro21-Ri21)
+// with uniform wear
+Rfw1=(Ro1+Ri1)/2
+//-----------------------------------------------------------------------------
+
+//-----------------------------------------------------------------------------
+//Case 2: Ro=100mm ; Ri= 25mm
+Ro2=100//mm
+Ri2=25//mm
+Ro32=Ro2^3
+Ri32=Ri2^3
+Ro22=Ro2^2
+Ri22=Ri2^2
+// with uniform pressure:
+Rfp2=((2/3)*(Ro32-Ri32))/(Ro22-Ri22)
+// with uniform wear
+Rfw2=(Ro2+Ri2)/2
+//-----------------------------------------------------------------------------
+
+//-----------------------------------------------------------------------------
+//Printing result file to .txt
+res3= mopen(TMPDIR+'/3_friction_radii_for_different_conditions.txt','wt');
+mfprintf(res3,'CASE I: Ro=100mm and Ri=90mm\n')
+mfprintf(res3,'The friction radius based on uniform pressure theory is %2.1f mm\n',Rfp1);
+mfprintf(res3,'The friction radius based on uniform wear theory is %2.1f mm\n',Rfw1);
+mfprintf(res3,"\n\nCASE II: Ro=100mm ; Ri= 25mm\n")
+mfprintf(res3,"The friction radius based on uniform pressure theory is %2.1f mm\n",Rfp2)
+mfprintf(res3,"The friction radius based on uniform wear theory is %2.1f mm",Rfw2)
+mclose(res3);
+editor(TMPDIR+'/3_friction_radii_for_different_conditions.txt')
+
+//------------------------------------------------------------------------------
+//--------------------------------End of program--------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.4/4_Determination_of_pressure_intensity.txt b/698/CH14/EX14.4/4_Determination_of_pressure_intensity.txt
new file mode 100644
index 000000000..9dbdad4f6
--- /dev/null
+++ b/698/CH14/EX14.4/4_Determination_of_pressure_intensity.txt
@@ -0,0 +1,11 @@
+Pressure variation is given by 
+p=C/r =F/(2*pi*(Ro-Ri)*r)
+
+Maximum pressure intensity occurs at inner radius r=Ri
+The maximum intensity of pressure is 254.6 kN/m^2
+
+Minimum pressure intensity occurs at outer radius r=Ro
+The minimum intensity of pressure is equal to 127.3 kN/m2
+
+Average intensity of pressure
+The average intensity of pressure is equal to 169.8 kN/m2
+\ No newline at end of file
diff --git a/698/CH14/EX14.4/P4_Determination_of_pressure_intensity.sce b/698/CH14/EX14.4/P4_Determination_of_pressure_intensity.sce
new file mode 100644
index 000000000..6b8d071eb
--- /dev/null
+++ b/698/CH14/EX14.4/P4_Determination_of_pressure_intensity.sce
@@ -0,0 +1,39 @@
+clc
+//Example 14.4
+// Determination of pressure intensity
+
+//------------------------------------------------------------------------------
+//Given Data:
+//Axial force
+F=4000//N
+//Inner radius
+Ri=0.050//m
+//Outer radius
+Ro=0.100//m
+//------------------------------------------------------------------------------
+
+//Maximum pressure intensity occurs at inner radius
+Pmax=F/(2*%pi*(Ro-Ri)*Ri)
+
+//Minimum pressure intensity occurs at outer radius
+Pmin=F/(2*%pi*(Ro-Ri)*Ro)
+
+//Average intensity of pressure
+Pavg=F/(%pi*(Ro^2-Ri^2))
+//------------------------------------------------------------------------------
+
+//Printing result file to .txt
+res4= mopen(TMPDIR+'/4_Determination_of_pressure_intensity.txt','wt')
+mfprintf(res4,'Pressure variation is given by \n')
+mfprintf(res4,'p=C/r =F/(2*pi*(Ro-Ri)*r)\n\n')
+mfprintf(res4,'Maximum pressure intensity occurs at inner radius r=Ri\n')
+mfprintf(res4,'The maximum intensity of pressure is %3.1f kN/m^2\n\n',(Pmax*(10^-3)))
+mfprintf(res4,'Minimum pressure intensity occurs at outer radius r=Ro\n')
+mfprintf(res4,'The minimum intensity of pressure is equal to %3.1f kN/m2\n\n",(Pmin*(10^-3)))
+mfprintf(res4,'Average intensity of pressure\n')
+mfprintf(res4,'The average intensity of pressure is equal to %3.1f kN/m2',(Pavg*(10^-3)))
+mclose(res4);
+editor(TMPDIR+'/4_Determination_of_pressure_intensity.txt')
+
+//------------------------------------------------------------------------------
+//--------------------------------End of program--------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.5/5_relation_of_friction_radii.png b/698/CH14/EX14.5/5_relation_of_friction_radii.png
new file mode 100644
index 000000000..5afe2184b
--- /dev/null
+++ b/698/CH14/EX14.5/5_relation_of_friction_radii.png
diff --git a/698/CH14/EX14.5/P5_Relation_of_friction_radii.sce b/698/CH14/EX14.5/P5_Relation_of_friction_radii.sce
new file mode 100644
index 000000000..d735ab1bf
--- /dev/null
+++ b/698/CH14/EX14.5/P5_Relation_of_friction_radii.sce
@@ -0,0 +1,34 @@
+clc
+//Example 14.5
+//Plot ratio of friction radius to outer radius versus ratio of inside radius to outer radius [Rf/Ro vs Ri/Ro]
+//Assuming outer radius=100mm and a matrix r of multiplying constants 
+
+r=[0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9]
+Ro=100
+n=length(r)
+
+// Assuming uniform pressure
+for i=1:n
+    Ri1(i)=r(i)*Ro
+    Rf1(i)=(2/3)*(((Ro^3)-(Ri1(i)^3))/((Ro^2)-(Ri1(i)^2)))
+    xa(i)=Ri1(i)/Ro
+    ya(i)=Rf1(i)/Ro
+end
+
+//Assuming uniform wear
+for j=1:n
+    Ri2(j)=r(j)*Ro
+    Rf2(j)=(Ro+Ri2(j))/2
+    xb(j)=Ri2(j)/Ro
+    yb(j)=Rf2(j)/Ro
+end
+
+//Plotting
+plot (xa,ya,'--r')
+plot (xb,yb,'-b')
+xgrid(0.1)
+title("Plot showing relation of friction radii for friction clutches; uniform pressure compared with uniform wear",'fontsize',4)
+xlabel("Ri/Ro= Inside radius/Outside radius",'fontsize',3)
+ylabel("Rf/Ro=Friction radius/Outside radius",'fontsize',3)
+legend('By uniform pressure Rf=(2/3)*(((Ro^3)-(Ri^3))/((Ro^2)-(Ri^2)))','By uniform wear Rf=(Ro+Ri)/2')
+//-----End of program--------
+\ No newline at end of file
diff --git a/698/CH14/EX14.6/6_determination_of_torque.txt b/698/CH14/EX14.6/6_determination_of_torque.txt
new file mode 100644
index 000000000..52be13896
--- /dev/null
+++ b/698/CH14/EX14.6/6_determination_of_torque.txt
@@ -0,0 +1,9 @@
+(a) Considering uniform wear theory, the torque is 1044.91 N-m
+
+(b) Considering uniform pressure theory, the torque is 1099.35 N-m
+
+(c) From above two values of torque, it is evident that the uniform wear theory
+gives conservative result, which is better for designing of the clutch.
+As for the new clutch where uniform pressure condition prevails, the torque capacity
+will be more than designed and when the clutch becomes old,
+(where the uniform wear prevails) the clutch is still capable of transmitting required torque.
+\ No newline at end of file
diff --git a/698/CH14/EX14.6/P6_Determination_of_torque.sce b/698/CH14/EX14.6/P6_Determination_of_torque.sce
new file mode 100644
index 000000000..d9b9645a8
--- /dev/null
+++ b/698/CH14/EX14.6/P6_Determination_of_torque.sce
@@ -0,0 +1,39 @@
+clc
+//exercise 14.6
+//Determination of Torque transmitting capacity
+
+//------------------------------------------------------------------------------
+//Given Data:
+//Outside diameter
+Do=0.3//m
+Ro=Do/2//m (outside radius)
+//Inside diameter
+Di=0.13//m
+Ri=Di/2//m (inside radius)
+//coefficient of friction
+f=0.2
+//Maximum allowable pressure
+pmax=0.7*(10^6)//Pa
+//Number of surfaces
+n=2
+//Axial force
+F=2*%pi*pmax*Ri*(Ro-Ri)
+//------------------------------------------------------------------------------
+
+// Torque With Uniform Wear Theory
+Tw=(1/2)*(f*F*n)*(Ro+Ri)
+
+// Torque With Uniform Pressure Theory
+Tp=(2/3)*(f*F*n)*(((Ro^3)-(Ri^3))/((Ro^2)-(Ri^2)))
+
+//------------------------------------------------------------------------------
+//Printing the result file to .txt
+res6= mopen(TMPDIR+'/6_determination_of_torque.txt','wt');
+mfprintf(res6,"(a) Considering uniform wear theory, the torque is %0.2f N-m\n\n",Tw)
+mfprintf(res6,"(b) Considering uniform pressure theory, the torque is %0.2f N-m\n\n",Tp)
+mfprintf(res6,'(c) From above two values of torque, it is evident that the uniform wear theory\ngives conservative result, which is better for designing of the clutch.\nAs for the new clutch where uniform pressure condition prevails, the torque capacity\nwill be more than designed and when the clutch becomes old,\n(where the uniform wear prevails) the clutch is still capable of transmitting required torque.')
+mclose(res6);
+editor(TMPDIR+'/6_determination_of_torque.txt')
+
+//------------------------------------------------------------------------------
+//------------------------------End of program----------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.7/7_determination_of_power_capacity.txt b/698/CH14/EX14.7/7_determination_of_power_capacity.txt
new file mode 100644
index 000000000..7d1511062
--- /dev/null
+++ b/698/CH14/EX14.7/7_determination_of_power_capacity.txt
@@ -0,0 +1 @@
+Power capacity is 1.10 kW
+\ No newline at end of file
diff --git a/698/CH14/EX14.7/P7_determination_of_power_capacity.sce b/698/CH14/EX14.7/P7_determination_of_power_capacity.sce
new file mode 100644
index 000000000..8bd12432a
--- /dev/null
+++ b/698/CH14/EX14.7/P7_determination_of_power_capacity.sce
@@ -0,0 +1,39 @@
+clc
+//Example 14.7
+//To determine power capacity
+
+//------------------------------------------------------------------------------
+//Given Data:
+//number of steel disks
+n1=4
+//number of bronze disks
+n2=3
+//Area of contact
+a=2.5*(10^-3)//m^2
+//Mean radius
+Rf=50/1000//m
+//coefficient of friction
+f=0.25
+//Axial force
+F=350//N
+//speed
+N=400//rpm
+//------------------------------------------------------------------------------
+
+//Assuming uniform wear
+//total number of contacting surfaces
+n=n1+n2-1
+//Torque
+T=f*F*n*Rf
+//Power capacity
+P=(2*%pi*N*T)/60
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res7=mopen(TMPDIR+'7_determination_of_power_capacity.txt','wt')
+mfprintf(res7,'Power capacity is %0.2f kW',P*(10^-3))
+mclose(res7)
+editor(TMPDIR+'7_determination_of_power_capacity.txt')
+
+//------------------------------------------------------------------------------
+//------------------------------End of program----------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.8/8_determination_of_number_of_plates.txt b/698/CH14/EX14.8/8_determination_of_number_of_plates.txt
new file mode 100644
index 000000000..80fd06c12
--- /dev/null
+++ b/698/CH14/EX14.8/8_determination_of_number_of_plates.txt
@@ -0,0 +1,5 @@
+Number of contacting surfaces 11
+Number of driver plates are 6
+Number of friction plates are 6
+
+There are 11 total number of plates, with 6 pressure and 6 friction discs
+\ No newline at end of file
diff --git a/698/CH14/EX14.8/P8_Determination_of_number_of_plates.sce b/698/CH14/EX14.8/P8_Determination_of_number_of_plates.sce
new file mode 100644
index 000000000..4e630dc29
--- /dev/null
+++ b/698/CH14/EX14.8/P8_Determination_of_number_of_plates.sce
@@ -0,0 +1,52 @@
+clc
+// Example 14.8
+//Determine number of plates
+//------------------------------------------------------------------------------
+//Given Data:
+//Torque
+t=45//Nm
+//Service factor
+Ka=2
+//Outer Radius
+Ro=0.1/2//m
+//Inner radius
+Ri=0.065/2//m
+//Mean raius
+Rm=(Ro+Ri)/2//m
+//coefficient of friction
+f=0.08
+//Maximum allowable pressure
+pmax=0.7*(10^6)//Pa
+//Design torque
+T=Ka*t
+
+//------------------------------------------------------------------------------
+//Assuimg uniform wear conditions
+//Axial force
+F=2*%pi*pmax*Ri*(Ro-Ri)
+//Number of contacting surfaces
+n=T/(f*F*Rm)
+N=ceil(n)
+
+//n1= number of driver plates
+//n2=number of friction plates
+//n=n1+n2-1
+//n2=n1-1 or n2=n1
+//solving above two equations
+n1=(N/2)
+N1=ceil(n1)
+n2=n1
+N2=ceil(n2)
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res8=mopen(TMPDIR+'8_determination_of_number_of_plates.txt','wt')
+mfprintf (res8,"Number of contacting surfaces %d\n",N)
+mfprintf(res8,"Number of driver plates are %d\n",N1)
+mfprintf(res8,"Number of friction plates are %d\n\n",N2)
+mfprintf(res8,"There are %d total number of plates, with %d pressure and %d friction discs",N,N1,N2)
+mclose(res8)
+editor(TMPDIR+'8_determination_of_number_of_plates.txt')
+
+//------------------------------------------------------------------------------
+//------------------------------End of program----------------------------------
+\ No newline at end of file
diff --git a/698/CH14/EX14.9/9_parameters_of_multiplate_clutch.txt b/698/CH14/EX14.9/9_parameters_of_multiplate_clutch.txt
new file mode 100644
index 000000000..9b1a12886
--- /dev/null
+++ b/698/CH14/EX14.9/9_parameters_of_multiplate_clutch.txt
@@ -0,0 +1,6 @@
+(a) Torque capacity of one pair of surfaces, assuming uniform wear is 20.0 Nm
+(b) Total torque applied 50.9 Nm
+(c) Number of contacting surfaces 4
+(d) Actual torque per pair of surfaces is 12.7 Nm
+(e) Actual axial force required is 2315.0 N
+(f) Maximum torque occurs at inner radius,and is equal to 307.0 kN/m^2
+\ No newline at end of file
diff --git a/698/CH14/EX14.9/P9_Parameters_of_multiplate_clutch.sce b/698/CH14/EX14.9/P9_Parameters_of_multiplate_clutch.sce
new file mode 100644
index 000000000..050eed9e5
--- /dev/null
+++ b/698/CH14/EX14.9/P9_Parameters_of_multiplate_clutch.sce
@@ -0,0 +1,58 @@
+clc
+//Example 14.9
+//Parameters of a multiplate clutch
+//------------------------------------------------------------------------------
+//Given Data:
+//Power capacity
+P=4000 //Watts
+//speed
+N=750 //rev/min
+//Inner radius
+Ri=0.04 //m
+//Outer radius
+Ro=0.07 //m
+//coefficient of friction
+f=0.1
+//average allowable pressure
+pmax=350*(10^3) // N/(m^2)
+//------------------------------------------------------------------------------
+
+// Assuming uniform wear theory
+//Axial force
+F=pmax*%pi*((Ro^2)-(Ri^2))
+//Total torque applied
+T=(P*60)/(2*%pi*N)
+//Torque per pair of surfaces
+Tc=F*f*((Ro+Ri)/2)
+//number of surfaces
+n=T/Tc
+n=ceil(n)
+//This obtained value must be an even number since it is the total number of surfaces, which are in pairs
+if (n-fix(n./2).*2==0)//For checking if n is even or not. If not, adding 1 to n to make it even
+    n=n
+else 
+    n=n+1
+end
+//Actual Torque
+Ta=T/n
+//Actual axial force required:
+Fa=(2*Ta)/(f*(Ro+Ri))
+//Average actual pressure:
+Pa=Fa/(%pi*((Ro^2)-(Ri^2)))
+//Maximum actual pressure
+Pmax=Fa/(2*%pi*(Ro-Ri)*Ri)
+
+//------------------------------------------------------------------------------
+//Printing result file to .txt
+res9=mopen(TMPDIR+'9_parameters_of_multiplate_clutch.txt','wt')
+mfprintf(res9,"(a) Torque capacity of one pair of surfaces, assuming uniform wear is %0.1f Nm\n",Tc)
+mfprintf(res9,"(b) Total torque applied %0.1f Nm\n",T)
+mfprintf(res9,"(c) Number of contacting surfaces %d\n",n)
+mfprintf(res9,"(d) Actual torque per pair of surfaces is %0.1f Nm\n",Ta)
+mfprintf(res9,"(e) Actual axial force required is %0.1f N\n",Fa)
+mfprintf(res9,"(f) Maximum torque occurs at inner radius,and is equal to %0.1f kN/m^2",(Pmax*(10^-3)))
+mclose(res9)
+editor(TMPDIR+'9_parameters_of_multiplate_clutch.txt')
+
+//------------------------------------------------------------------------------
+//------------------------------------End of program----------------------------
+\ No newline at end of file