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--------------------------------