1 files changed, 23 insertions, 0 deletions
diff --git a/2132/CH3/EX3.21/Example3_21.sce b/2132/CH3/EX3.21/Example3_21.sce
new file mode 100755
index 000000000..f2e65e250
--- /dev/null
+++ b/2132/CH3/EX3.21/Example3_21.sce
@@ -0,0 +1,23 @@
+//Example 3.21
+clc;
+clear;
+close;
+format('v',7);
+//Given data :
+r=2;//meter
+l=4;//meter
+A=r*l;//m^2
+xbar=2+r/2;//meter
+w=9.81;//kN/m^2
+PH=w*A*xbar;//kN
+disp(PH,"Horizontal component of resulting Pressure in kN : ");
+PV=2*r*l*w+%pi*r^2/4*l*w;//kN
+disp(PV,"Verticalal component of resulting Pressure in kN : ");
+IG=(l*r^3)/12;//in m^4
+h_bar=IG/A/xbar+xbar;//in meter
+disp(h_bar,"Position of centre of horizontal component of pressure in meter : ");
+x=(2*r+%pi*r^2/4*(4*r/3/%pi))/(2*r+%pi*r^2/4);//meter
+P=sqrt(PH^2+PV^2);//kN
+disp(P,"Resultant pressure in kN : ");
+theta=atand(PV/PH);//degree
+disp(theta,"Direction of resultant pressure in degree : ");