initial commit / add all books

4 files changed, 114 insertions, 0 deletions

diff --git a/2732/CH3/EX3.3/Ex3_3.sce b/2732/CH3/EX3.3/Ex3_3.sce
new file mode 100755
index 000000000..64fdaf34c
--- /dev/null
+++ b/2732/CH3/EX3.3/Ex3_3.sce
@@ -0,0 +1,32 @@
+clc

+// initialization of variables

+clear

+epsillon=[0.01 0 0

+     0 0.02 0.02

+     0 0.02 0.01] // dimensionless

+theta=30 // degrees

+//calculations

+theta=theta*%pi/180

+a=[cos(theta) sin(theta) 0

+  -sin(theta) cos(theta) 0

+   0 0 1]

+b=a.'

+epsillon_new=a*epsillon*b

+// calculation of strain invariants

+// for epsillon

+J1=epsillon(1,1)+epsillon(2,2)+epsillon(3,3)

+J2=epsillon(1,1)*epsillon(2,2)+epsillon(2,2)*epsillon(3,3)+epsillon(3,3)*epsillon(1,1)-2*(epsillon(1,2)^2+epsillon(2,3)^2+epsillon(3,1)^2)

+J3=epsillon(1,1)*epsillon(2,2)*epsillon(3,3)+2*epsillon(1,2)*epsillon(2,3)*epsillon(3,1)-(epsillon(1,1)*epsillon(2,3)^2+epsillon(2,2)*epsillon(3,1)^2+epsillon(3,3)*epsillon(1,2)^2)

+

+// for epsillon_new

+J11=epsillon_new(1,1)+epsillon_new(2,2)+epsillon_new(3,3)

+J22=epsillon_new(1,1)*epsillon_new(2,2)+epsillon_new(2,2)*epsillon_new(3,3)+epsillon_new(3,3)*epsillon_new(1,1)-2*(epsillon_new(1,2)^2+epsillon_new(2,3)^2+epsillon_new(3,1)^2)

+J33=epsillon_new(1,1)*epsillon_new(2,2)*epsillon_new(3,3)+2*epsillon_new(1,2)*epsillon_new(2,3)*epsillon_new(3,1)-(epsillon_new(1,1)*epsillon_new(2,3)^2+epsillon_new(2,2)*epsillon_new(3,1)^2+epsillon_new(3,3)*epsillon_new(1,2)^2)

+

+// results

+printf('The new strain tensor is');

+disp(epsillon_new);

+printf('The strain invariants of old strain tensor are J1=%0.2f J2=%.e J3=%.e \n and that of the new strain tensor are J1=%0.2f J2=%.e J3=%.e',J1,J2,J3,J11,J22,J33)

+printf('\n Hence the same strain invariants')

+

+

diff --git a/2732/CH3/EX3.4/Ex3_4.sce b/2732/CH3/EX3.4/Ex3_4.sce
new file mode 100755
index 000000000..1029d1c60
--- /dev/null
+++ b/2732/CH3/EX3.4/Ex3_4.sce
@@ -0,0 +1,36 @@
+clc

+// initialization of variables

+clear

+epsillon=[0.01 -0.02 0

+     -0.02 0.03 -0.01

+     0 -0.01 0] // dimensionless

+a_xx=0.6     

+theta=acos(a_xx) // radians

+//calculations

+// theta=theta*%pi/180

+a=[cos(theta) 0 -sin(theta)

+   0       1        0

+   sin(theta) 0 cos(theta)]

+b=a.'

+epsillon_new=a*epsillon*b

+

+// calculation of strain invariants

+// for epsillon

+J1=epsillon(1,1)+epsillon(2,2)+epsillon(3,3)

+

+J2=epsillon(1,1)*epsillon(2,2)+epsillon(2,2)*epsillon(3,3)+epsillon(3,3)*epsillon(1,1)-2*(epsillon(1,2)^2+epsillon(2,3)^2+epsillon(3,1)^2)

+

+J3=epsillon(1,1)*epsillon(2,2)*epsillon(3,3)+2*epsillon(1,2)*epsillon(2,3)*epsillon(3,1)-(epsillon(1,1)*epsillon(2,3)^2+epsillon(2,2)*epsillon(3,1)^2+epsillon(3,3)*epsillon(1,2)^2)

+

+// for epsillon_new

+J11=epsillon_new(1,1)+epsillon_new(2,2)+epsillon_new(3,3)

+

+J22=epsillon_new(1,1)*epsillon_new(2,2)+epsillon_new(2,2)*epsillon_new(3,3)+epsillon_new(3,3)*epsillon_new(1,1)-2*(epsillon_new(1,2)^2+epsillon_new(2,3)^2+epsillon_new(3,1)^2)

+

+J33=epsillon_new(1,1)*epsillon_new(2,2)*epsillon_new(3,3)+2*epsillon_new(1,2)*epsillon_new(2,3)*epsillon_new(3,1)-(epsillon_new(1,1)*epsillon_new(2,3)^2+epsillon_new(2,2)*epsillon_new(3,1)^2+epsillon_new(3,3)*epsillon_new(1,2)^2)

+

+// Results

+printf('The new strain tensor is')

+disp(epsillon_new)

+printf('The strain invariants of old stress tensor are J1=%0.2f J2=%.e J3=%.e \n and that of the new stress tensor are J1=%0.2f J2=%.e J3=%.e',J1,J2,J3,J11,J22,J33)

+

diff --git a/2732/CH3/EX3.5/Ex3_5.sce b/2732/CH3/EX3.5/Ex3_5.sce
new file mode 100755
index 000000000..7dc509127
--- /dev/null
+++ b/2732/CH3/EX3.5/Ex3_5.sce
@@ -0,0 +1,19 @@
+clc

+// initialization of variables

+clear

+epsillon_A= 700*10^-6 

+epsillon_B= 300*10^-6 

+epsillon_C= 300*10^-6 

+theta=45 // degrees

+theta=theta*%pi/180 // radians

+// calculations

+epsillon_x=epsillon_A

+epsillon_y=epsillon_C

+gamma_xy=(epsillon_B-(epsillon_x*cos(theta)^2+epsillon_y*sin(theta)^2))/(sin(theta)*cos(theta))

+epsillon_1=1/2*(epsillon_x+epsillon_y)+(1/2)*sqrt((epsillon_x-epsillon_y)^2+gamma_xy^2)

+epsillon_2=1/2*(epsillon_x+epsillon_y)-(1/2)*sqrt((epsillon_x-epsillon_y)^2+gamma_xy^2)

+phi=0.5*atan(gamma_xy/(epsillon_x-epsillon_y))

+phi=phi*180/%pi

+//results

+printf('The principal strains are %.3e, %.3e',epsillon_1,epsillon_2)

+printf('\n phi = %.2f degrees',phi)

diff --git a/2732/CH3/EX3.6/Ex3_6.sce b/2732/CH3/EX3.6/Ex3_6.sce
new file mode 100755
index 000000000..366be3d70
--- /dev/null
+++ b/2732/CH3/EX3.6/Ex3_6.sce
@@ -0,0 +1,27 @@
+clc

+// initialization of variables

+clear

+epsillon_A= 1000*10^-6 

+epsillon_B= 720*10^-6 

+epsillon_C= 600*10^-6 

+th_B=120 // degrees

+th_C=240 // degrees

+//calculations

+th_B=th_B*%pi/180

+th_C=th_C*%pi/180

+// we need to solve for epsillon_y and gamma_xy

+// Ax=B

+ep_x=epsillon_A

+A=[sin(th_B)^2 sin(th_B)*cos(th_B)

+   sin(th_C)^2 sin(th_C)*cos(th_C)]

+C=[epsillon_B-ep_x*cos(th_B)^2; epsillon_C-ep_x*cos(th_C)^2]   

+x=inv(A)*C

+ep_y=x(1,1)

+gam_xy=x(2,1)

+epsillon_x=ep_x

+epsillon_y=ep_y

+gamma_xy=gam_xy

+epsillon_1=1/2*(epsillon_x+epsillon_y)+(1/2)*sqrt((epsillon_x-epsillon_y)^2+gamma_xy^2)

+epsillon_2=1/2*(epsillon_x+epsillon_y)-(1/2)*sqrt((epsillon_x-epsillon_y)^2+gamma_xy^2)

+// Results

+printf('The principal strains are %.3e, %.3e',epsillon_1,epsillon_2)