2 files changed, 48 insertions, 0 deletions
diff --git a/2090/CH13/EX13.1/Chapter13_example1.sce b/2090/CH13/EX13.1/Chapter13_example1.sce
new file mode 100755
index 000000000..898b4d2a8
--- /dev/null
+++ b/2090/CH13/EX13.1/Chapter13_example1.sce
@@ -0,0 +1,23 @@
+clc
+clear
+//Input data
+d=0.08;//The diameter of bore in m
+L=0.075;//The length of the stroke in m
+l=0.152;//The connecting rod length in m
+h=0.062;//Skirt length of the piston in m
+Fr=8000;//Compressive force in the connecting rod in N
+p=3000;//The pressure in the cylinder kPa
+y=0.004*10^-3;//The clearence between piston and cylinder wall in m
+U=0.006;//The dynamic viscosity of the lubricating oil in pa.s
+u=8.2;//The piston speed in m/s
+pi=3.141;//Mathematical constant of pi
+
+//Calculations
+ts=(U*u)/y;//The shear stress in N/m^2
+A=pi*d*h;//Contact area between the piston and the cylinder in m^2
+Ff=ts*A;//Friction force on the piston inN
+r=L/2;//Crank length in m
+A=atan(r/l);//The angle made by the crank in radians
+Ft=Fr*sin(A);//The side thrust in N
+//Output 
+printf(' The friction force on the piston = %3.0f N \n  The thrust force on the cylinder wall = %3.0f N',Ff,Ft)
diff --git a/2090/CH13/EX13.2/Chapter13_example2.sce b/2090/CH13/EX13.2/Chapter13_example2.sce
new file mode 100755
index 000000000..c21ba2028
--- /dev/null
+++ b/2090/CH13/EX13.2/Chapter13_example2.sce
@@ -0,0 +1,25 @@
+clc
+clear
+//Input data
+d=0.065;//The cylinder bore diameter in m
+L=6;//The length of the stroke in cm
+l=12;//The length of the connecting rod in cm
+p=50;//The pressure in the cylinder in bar
+q=90;//The crank position in power stroke of the cycle for one cylinder in degrees
+Ff=900;//Friction force in N
+pi=3.141;//Mathematical constant valu of pi
+o=0.2;//Wrist pin off set in cm
+
+//Calculations
+r=L/2;//The crank length in cm
+sineA=r/l;//The value of sine 
+cosA=(1-(sineA)^2)^(1/2);//The value of cosine
+Fr=[[(p*10^5*(pi/4)*d^2)-Ff]/cosA]/1000;//The force in the connecting rod in kN
+Ft=Fr*sineA;//The side thrust on the piston in kN
+sineA1=(r-o)/l;//The value of sine
+cosA1=(1-(sineA1)^2)^(1/2);//The value of cosine
+Fr1=[[(p*10^5*(pi/4)*d^2)-Ff]/cosA1]/1000;//The force in the connecting rod in kN
+Ft1=Fr1*sineA1;//The side thrust in kN
+
+//Output 
+printf('(a) The force in the connecting rod = %3.3f kN \n The side thrust on the piston = %3.3f kN \n (b) The side thrust on the piston = %3.3f kN ',Fr,Ft,Ft1)