4 files changed, 88 insertions, 0 deletions
diff --git a/3863/CH3/EX3.13/Ex3_13.sce b/3863/CH3/EX3.13/Ex3_13.sce
new file mode 100644
index 000000000..658228413
--- /dev/null
+++ b/3863/CH3/EX3.13/Ex3_13.sce
@@ -0,0 +1,19 @@
+clear
+//
+
+//Given
+//Variable declaration
+sigma1=120       //Major tensile stress in N/sq.mm
+sigma2=-90       //Minor compressive stress in N/sq.mm
+sigma_gp=150     //Greatest principal stress in N/sq.mm
+
+//Calculation
+ //case(a):Magnitude of the shearing stresses on the two planes
+tau=(sqrt(((sigma_gp-((sigma1+sigma2)/2))**2)-(((sigma1-sigma2)/2)**2)))
+
+ //case(b):Maximum shear stress at the point
+sigmat_max=int((sqrt((sigma1-sigma2)**2+(4*tau**2)))/2)
+
+//Result
+printf("\n Shear stress on the two planes = %0.3f  N/mm^2",tau)
+printf("\n Maximum shear stress at the point = %0.3f  N/mm^2",sigmat_max)
diff --git a/3863/CH3/EX3.16/Ex3_16.sce b/3863/CH3/EX3.16/Ex3_16.sce
new file mode 100644
index 000000000..38a814b8a
--- /dev/null
+++ b/3863/CH3/EX3.16/Ex3_16.sce
@@ -0,0 +1,24 @@
+clear
+//
+
+//Given
+//Variable declaration
+sigma1=600        //Major tensile stress in N/sq.mm
+sigma2=300        //Minor tensile stress in N/sq.mm
+tau=450           //Shear stress in N/sq.mm
+theta1=45         //Angle of inclination in degrees
+theta2=135        //Angle of inclination in degrees
+
+//Calculation
+sigman1=int(((sigma1+sigma2)/2)+(((sigma1-sigma2)/2)*cos((%pi/180)*(2*theta1)))+(tau*sin((%pi/180)*(2*theta1))))
+sigman2=int(((sigma1+sigma2)/2)+(((sigma1-sigma2)/2)*cos((%pi/180)*(2*theta2)))+(tau*sin((%pi/180)*(2*theta2))))
+sigmat1=int(((sigma1-sigma2)/2*(sin((%pi/180)*(2*theta1)))-(tau*cos((%pi/180)*(2*theta1)))))
+
+sigmat2=int(((sigma1-sigma2)/2*(sin((%pi/180)*(2*theta2)))-(tau*cos((%pi/180)*(2*theta2)))))
+
+
+//Result
+printf("\n Normal stress(when theta is 45 degrees)= %0.3f  N/mm^2",sigman1)
+printf("\n Normal stress(when theta is 135 degrees)= %0.3f  N/mm^2",sigman2)
+printf("\n Tangential stress(when theta is 45 degrees)= %0.3f  N/mm^2",sigmat1)
+printf("\n Tangential stress(when theta is 135 degrees)= %0.3f  N/mm^2",sigmat2)
diff --git a/3863/CH3/EX3.8/Ex3_8.sce b/3863/CH3/EX3.8/Ex3_8.sce
new file mode 100644
index 000000000..0c4c4b647
--- /dev/null
+++ b/3863/CH3/EX3.8/Ex3_8.sce
@@ -0,0 +1,23 @@
+clear
+//
+
+//Given
+//Variable declaration
+sigma1=100           //Major principal stress in N/sq.mm
+sigma2=-60           //Minor principal stress in N/sq.mm
+theta=90-50          //Angle of inclination in degrees
+
+//Calculation
+sigman=(((sigma1+sigma2)/2)+(((sigma1-sigma2)/2)*cos((%pi/180)*(2*theta))))
+
+sigmat=((sigma1-sigma2)/2*(sin((%pi/180)*(2*theta))))
+
+sigmaR=(sqrt(sigman**2+sigmat**2))
+
+sigmat_max=int((sigma1-sigma2)/2)
+
+//Result
+printf("\n Normal stress = %0.3f  N/mm^2",sigman)
+printf("\n Shear stress = %0.3f  N/mm^2",sigmat)
+printf("\n Resultant stress = %0.3f  N/mm^2",sigmaR)
+printf("\n Maximum shear stress = %0.3f  N/mm^2",sigmat_max)
diff --git a/3863/CH3/EX3.9/Ex3_9.sce b/3863/CH3/EX3.9/Ex3_9.sce
new file mode 100644
index 000000000..8aacb78f1
--- /dev/null
+++ b/3863/CH3/EX3.9/Ex3_9.sce
@@ -0,0 +1,22 @@
+clear
+//
+
+//Given
+//Variable declaration
+sigma1=100        //Major principal stress in N/sq.mm
+sigma2=-40        //Minor principal stress in N/sq.mm
+theta=90-60       //Angle of inclination in degrees
+
+//Calculation
+sigman=((sigma1+sigma2)/2)+(((sigma1-sigma2)/2)*cos((%pi/180)*(2*theta)))
+sigmat=((sigma1-sigma2)/2*(sin((%pi/180)*(2*theta))))
+
+sigmaR=(sqrt(sigman**2+sigmat**2))
+
+sigmat_max=int((sigma1-sigma2)/2)
+phi=int((180/%pi)*(atan(sigmat/sigman)))
+
+//Result
+printf("\n Resultant stress in magnitude = %0.3f  N/mm^2",sigmaR)
+printf("\n Direction of resultant stress = %0.3f  degrees",phi)
+printf("\n Maximum shear stress = %0.3f  N/mm^2",sigmat_max)