summaryrefslogtreecommitdiff
path: root/1898/CH9/EX9.9/Ex9_9.sce
diff options
context:
space:
mode:
Diffstat (limited to '1898/CH9/EX9.9/Ex9_9.sce')
-rwxr-xr-x1898/CH9/EX9.9/Ex9_9.sce32
1 files changed, 32 insertions, 0 deletions
diff --git a/1898/CH9/EX9.9/Ex9_9.sce b/1898/CH9/EX9.9/Ex9_9.sce
new file mode 100755
index 000000000..317baf009
--- /dev/null
+++ b/1898/CH9/EX9.9/Ex9_9.sce
@@ -0,0 +1,32 @@
+clear all; clc;
+
+disp("Scilab Code Ex 9.9 : ")
+
+//Given:
+sigma_x = -12; //MPa
+sigma_y = 0;
+tou_xy = -6; //MPa
+
+//Construction of the circle:
+sigma_avg = (sigma_x+sigma_y)/2;
+R = sqrt((-sigma_x+sigma_avg)^2 + (tou_xy)^2);
+
+//Principal Stresses:
+sigma2 = -R+sigma_avg;//From the Mohr's circle
+sigma1 = R+sigma_avg;
+
+theta_p2 = atan((-tou_xy)/(-sigma_x+sigma_avg));
+theta_p2 = theta_p2/2*(180/%pi);
+
+//Display:
+
+printf('\n\nThe first principal stress is = %1.2f MPa',sigma1);
+printf('\nThe second principal stress is = %1.1f MPa',sigma2);
+printf('\nThe direction of the principal plane is = %1.1f degrees',theta_p2);
+
+
+//----------------------------------------------------------------------END--------------------------------------------------------------------------------
+
+
+
+