7 files changed, 147 insertions, 0 deletions
diff --git a/2732/CH2/EX2.10/Ex2_10.sce b/2732/CH2/EX2.10/Ex2_10.sce
new file mode 100755
index 000000000..cd4721315
--- /dev/null
+++ b/2732/CH2/EX2.10/Ex2_10.sce
@@ -0,0 +1,15 @@
+clc
+// initialization of variables
+clear
+tau_xx= -1 // kgf/cm^2
+tau_yy= 0 // kgf/cm^2
+tau_xy= 7 // kgf/cm^2
+// calculations 
+sigma_1=(tau_xx+tau_yy)/2+sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+sigma_2=(tau_xx+tau_yy)/2-sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+x=sigma_1 // positive one is tension
+if(sigma_2>sigma_1)
+    x=sigma_2
+end
+// Results
+printf('The diagonal tension is %.3f kgf/cm^2',x)
diff --git a/2732/CH2/EX2.11/Ex2_11.sce b/2732/CH2/EX2.11/Ex2_11.sce
new file mode 100755
index 000000000..56fd8a754
--- /dev/null
+++ b/2732/CH2/EX2.11/Ex2_11.sce
@@ -0,0 +1,17 @@
+clc
+// initialization of variables
+clear
+d=2 // m
+l=10 // m
+t=1 // cm
+p=15 // kgf/cm^2
+pitch= 2*%pi //m
+//calculations
+w=2*%pi*d/2 // m
+theta=atan(w/(2*%pi))
+sigma_z=p*d*100/(4*t)
+sigma_th=p*d*100/(2*t)
+sigma_th_new=(sigma_th+sigma_z)/2+(sigma_th-sigma_z)/2*cos(2*theta)
+tau_thz=(sigma_z-sigma_th)*sin(2*theta)/2
+// results
+printf('At the junction, the normal and shear stresses are %d and %d kgf/cm^2 \n respectively, and the rivets must be designed for this',sigma_th_new,-tau_thz)
diff --git a/2732/CH2/EX2.4/Ex2_4.sce b/2732/CH2/EX2.4/Ex2_4.sce
new file mode 100755
index 000000000..0d01f3611
--- /dev/null
+++ b/2732/CH2/EX2.4/Ex2_4.sce
@@ -0,0 +1,20 @@
+clc
+// initialization of variables
+clear
+tau=[200 100 0
+     100 0 0
+     0 0 500] // some units
+theta=60 // degrees
+//calculations
+theta=theta*%pi/180
+a=[cos(theta) sin(theta) 0
+  -sin(theta) cos(theta) 0
+   0 0 1]
+b=a.'
+tau_new=a*tau*b
+// Results
+printf('The new stress tensor is')
+disp(tau_new)
+
+   
+
diff --git a/2732/CH2/EX2.5/Ex2_5.sce b/2732/CH2/EX2.5/Ex2_5.sce
new file mode 100755
index 000000000..cef88f7be
--- /dev/null
+++ b/2732/CH2/EX2.5/Ex2_5.sce
@@ -0,0 +1,10 @@
+clc
+// initialization of variables
+clear
+sigma_1=100 //kg*f/cm^2
+sigma_2=100 //kg*f/cm^2
+sigma_3=-200 //kg*f/cm^2
+// calculations
+tau_oct=1/3*sqrt((sigma_1-sigma_2)^2+(sigma_2-sigma_3)^2+(sigma_3-sigma_1)^2)
+// Results
+printf('Octahedra shear stress at the point is=%.1f kgf/cm^2',tau_oct)
diff --git a/2732/CH2/EX2.7/Ex2_7.sce b/2732/CH2/EX2.7/Ex2_7.sce
new file mode 100755
index 000000000..9bece8d3a
--- /dev/null
+++ b/2732/CH2/EX2.7/Ex2_7.sce
@@ -0,0 +1,31 @@
+clc
+// initialization of variables
+clear
+tau=[200 100 0
+     100 0 0
+     0 0 500] // some units
+theta=60 // degrees
+//calculations
+theta=theta*%pi/180
+a=[cos(theta) sin(theta) 0
+  -sin(theta) cos(theta) 0
+   0 0 1]
+b=a.'
+tau_new=a*tau*b
+
+// stress invariants :old 
+I1=tau(1,1)+tau(2,2)+tau(3,3)
+I2=tau(1,1)*tau(2,2)+tau(2,2)*tau(3,3)+tau(3,3)*tau(1,1)-(tau(1,2)^2+tau(2,3)^2+tau(3,1)^2)
+I3=tau(1,1)*tau(2,2)*tau(3,3)+2*tau(1,2)*tau(2,3)*tau(3,1)-(tau(1,1)*tau(2,3)^2+tau(2,2)*tau(3,1)^2+tau(3,3)*tau(1,2)^2)
+
+// stress invariants :new
+I11=tau_new(1,1)+tau_new(2,2)+tau_new(3,3)
+I22=tau_new(1,1)*tau_new(2,2)+tau_new(2,2)*tau_new(3,3)+tau_new(3,3)*tau_new(1,1)-(tau_new(1,2)^2+tau_new(2,3)^2+tau_new(3,1)^2)
+I33=tau_new(1,1)*tau_new(2,2)*tau_new(3,3)+2*tau_new(1,2)*tau_new(2,3)*tau_new(3,1)-(tau_new(1,1)*tau_new(2,3)^2+tau_new(2,2)*tau_new(3,1)^2+tau_new(3,3)*tau_new(1,2)^2)
+
+// Results
+printf('The invariants of old stress tensor are I1=%0.2f I2=%.e I3=%.e \n and that of the new stress tensor are I1=%0.2f I2=%.e I3=%.e',I1,I2,I3,I11,I22,I33)
+printf('\n Hence the same stress tensor invariants')
+
+
+
diff --git a/2732/CH2/EX2.8/Ex2_8.sce b/2732/CH2/EX2.8/Ex2_8.sce
new file mode 100755
index 000000000..717fd8045
--- /dev/null
+++ b/2732/CH2/EX2.8/Ex2_8.sce
@@ -0,0 +1,13 @@
+clc
+// initialization of variables
+clear
+sigma_3=0 // kgf/cm^2
+tau_oct=1500 // kgf/cm^2
+n=2 // given that sigma_1=n*sigma_2
+// calculations
+sigma_2=1500*3/(sqrt(2*n^2-2*n+2)) // // kgf/cm^2
+sigma_1=n*sigma_2 // kgf/cm^2 
+sigma_0=4500/sqrt(2) // kgf/cm^2
+// Results
+printf('The necessary stresses sigma_1, sigma_2 for biaxial yielding are \n %d kgf/cm^2, %d kgf/cm^2 and for uniaxial yielding sigma_0 %0.2f kgf/cm^2.',sigma_1,sigma_2,sigma_0);
+
diff --git a/2732/CH2/EX2.9/Ex2_9.sce b/2732/CH2/EX2.9/Ex2_9.sce
new file mode 100755
index 000000000..f5db1a9fd
--- /dev/null
+++ b/2732/CH2/EX2.9/Ex2_9.sce
@@ -0,0 +1,41 @@
+clc
+// initialization of variables
+clear
+// part (a)
+tau_xx=300 // kgf/cm^2
+tau_yy=0 // kgf/cm^2
+tau_xy=600 // kgf/cm^2
+// calculations
+sigma_1=(tau_xx+tau_yy)/2+sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+sigma_2=(tau_xx+tau_yy)/2-sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+Beta=atan(2*tau_xy/(tau_xx-tau_yy))
+Beta=Beta*180/%pi
+// Results
+printf('Part (a) \n The magnitude of principal stresses are %d,%d kgf/cm^2 and \n the direction is given by 2*beta=%.2f degrees',sigma_1,sigma_2,Beta)
+
+//part (b)
+tau_xx=1000 // kgf/cm^2
+tau_yy=150 // kgf/cm^2
+tau_xy=450 // kgf/cm^2
+// calculations
+sigma_1=(tau_xx+tau_yy)/2+sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+sigma_2=(tau_xx+tau_yy)/2-sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+Beta=atan(2*tau_xy/(tau_xx-tau_yy))
+Beta=Beta*180/%pi
+// Results
+printf('\n Part (b) \n The magnitude of principal stresses are %d,%d kgf/cm^2 and \n the direction is given by 2*beta=%.2f degrees',sigma_1,sigma_2,Beta)
+
+// part (c)
+tau_xx=-850 // kgf/cm^2
+tau_yy=350 // kgf/cm^2
+tau_xy=700 // kgf/cm^2
+// calculations
+sigma_1=(tau_xx+tau_yy)/2+sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+sigma_2=(tau_xx+tau_yy)/2-sqrt((1/2*(tau_xx-tau_yy))^2+tau_xy^2)
+Beta=atan(2*tau_xy/(tau_xx-tau_yy))
+Beta=Beta*180/%pi
+// Results
+printf('\n Part (c) \n The magnitude of principal stresses are %d,%d kgf/cm^2 and \n the direction is given by 2*beta=%0.2f',sigma_1,sigma_2,-Beta)
+
+// wrong answers were given in textbook for part (b)
+