1 files changed, 39 insertions, 0 deletions
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