2 files changed, 29 insertions, 0 deletions
diff --git a/1793/CH10/EX10.1/10q1.sce b/1793/CH10/EX10.1/10q1.sce
new file mode 100755
index 000000000..81acaf8c5
--- /dev/null
+++ b/1793/CH10/EX10.1/10q1.sce
@@ -0,0 +1,16 @@
+clc
+//initialisation of variables
+sx= 2000 //lb/ft^3
+sy= 2500 //lb/ft^3
+T= 800 //lb/ft^3
+t= 0.348//radians
+//calculations
+s1= (sx+sy)/2+sqrt(((sy-sx)/2)^2+T^2)
+s2= (sx+sy)/2-sqrt(((sy-sx)/2)^2+T^2)
+sn= (sx+sy)/2+(sy-sx)*cos(2*t)/2-T*sin(2*t)
+Tn= (sy-sx)*sin(2*t)/2+T*cos(2*t)
+//results
+printf ('principle stress s1 = % 2f lb/ft^3 ',s1)
+printf ('principle stress s2 = % 2f lb/ft^3 ',s2)
+printf ('normal stress = % 2f lb/ft^3 ',sn)
+printf ('shear stress = % 2f lb/ft^3 ',Tn)
diff --git a/1793/CH10/EX10.3/10q3.sce b/1793/CH10/EX10.3/10q3.sce
new file mode 100755
index 000000000..266b022a3
--- /dev/null
+++ b/1793/CH10/EX10.3/10q3.sce
@@ -0,0 +1,13 @@
+clc
+//initialisation of variables
+x= 3 //m
+y= 4 //m
+P= 5 //kN
+z= 2 //m
+//calculations
+r= sqrt(x^2+y^2)
+k= r/z
+I= 3/(2*%pi*((r/z)^2+1)^2.5)
+s= P*I/z^2
+//results
+printf ('verticle stress increase at 2m = % 4f kN/m^3 ',s)