1 files changed, 28 insertions, 0 deletions
diff --git a/1898/CH9/EX9.2/Ex9_2.sce b/1898/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..c34707ace
--- /dev/null
+++ b/1898/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,28 @@
+clear all; clc;
+
+
+disp("Scilab Code Ex 9.2 : ")
+
+//Given:
+phi = -30*(%pi/180);
+theta = 60*(%pi/180);
+sigma_x = -80; //MPa
+sigma_y = 50; //MPa
+tou_xy = -25; //MPa
+
+//Plane CD:
+sigma_x1 = (sigma_x+sigma_y)/2 + ((sigma_x-sigma_y)*cos(2*phi))/2 + (tou_xy*sin(2*phi)); //Eqn 9.1
+tou_xy1 = ((-(sigma_x - sigma_y)*sin(2*phi))/2) + (tou_xy*cos(2*phi)); //Eqn 9.2
+
+//Plane BC:
+sigma_x2 = (sigma_x+sigma_y)/2 + ((sigma_x-sigma_y)*cos(2*theta))/2 + (tou_xy*sin(2*theta)); //Eqn 9.1
+tou_xy2 = (-(sigma_x - sigma_y)*sin(2*theta))/2 + tou_xy*cos(2*theta); //Eqn 9.2
+
+//Display:
+
+printf('\n\nThe normal stress of plane CD inclined at 30 degrees    = %1.1f MPa',sigma_x1);
+printf('\nThe shear stress of plane CD inclined at 30 degrees     = %1.1f MPa',tou_xy1);
+printf('\nThe normal stress of plane BC inclined at 60 degrees    = %1.2f MPa',sigma_x2);
+printf('\nThe shear stress of plane BC inclined at 60 degrees     = %1.1f MPa',tou_xy2);
+
+//-------------------------------------------------------------------------END---------------------------------------------------------------------------------------