From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 3289/CH2/EX2.1/Ex2_1.sce | 30 ++++++++++++++++++++++++++++++
 3289/CH2/EX2.2/Ex2_2.sce | 37 +++++++++++++++++++++++++++++++++++++
 2 files changed, 67 insertions(+)
 create mode 100755 3289/CH2/EX2.1/Ex2_1.sce
 create mode 100755 3289/CH2/EX2.2/Ex2_2.sce

(limited to '3289/CH2')

diff --git a/3289/CH2/EX2.1/Ex2_1.sce b/3289/CH2/EX2.1/Ex2_1.sce
new file mode 100755
index 000000000..1c53af2b5
--- /dev/null
+++ b/3289/CH2/EX2.1/Ex2_1.sce
@@ -0,0 +1,30 @@
+clc;
+
+radius=((sqrt(195^2+130^2))*10^(-6));
+disp(radius,"radius of the circle in degree = ")
+theta1=(atand(130/195))
+disp(theta1,"pricipal stresses in degree");// displaying result
+epsilonx=510*10^(-6)
+epsilony=120*10^(-6)
+epsilon=(epsilonx+epsilony)/2
+disp(epsilon,"distance between O and c=")
+
+//solution a
+angle=60- theta1
+disp(angle,"angle of ACA1 in degree = ")// displaying result
+epsilonx1=epsilon+radius*cosd(26.3)
+disp(epsilonx1,"strains in x axis= ")// displaying result
+epsilony1=epsilon-radius*cosd(26.3)
+disp(epsilony1,"strains in y axis= ")// displaying result
+gammaxy=-2*(radius*sind(26.3))
+disp(gammaxy,"shear strain")// displaying result
+
+//solution b
+epsilon1=epsilon+radius
+disp(epsilon1,"strains in x axis= ")// displaying result
+epsilon2=epsilon-radius
+disp(epsilon2,"strains in x axis= ")// displaying result
+
+//solution c
+gammamax=-+468*10^(-6)
+disp(gammamax,"maxi shear stress= ")
diff --git a/3289/CH2/EX2.2/Ex2_2.sce b/3289/CH2/EX2.2/Ex2_2.sce
new file mode 100755
index 000000000..8a5b21d21
--- /dev/null
+++ b/3289/CH2/EX2.2/Ex2_2.sce
@@ -0,0 +1,37 @@
+clc;
+
+epsilon0=190*10^(-6)
+epsilon60=200*10^(-6)
+epsilon120=-300*10^(-6)
+E=200// GPa
+v=0.3
+epsilonx=epsilon0
+disp(epsilonx,"value of epsilonx is= ")
+
+// epsilon60=((epsilonx+epsilony)/2)-((epsilonx-epsilony)/4)+(gammaxy*sqrt(3))/4 eqn 1
+// epsilon120=((epsilonx+epsilony)/2)-((epsilonx-epsilony)/4)-(gammaxy*sqrt(3))/4 eqn 2
+
+epsilony=(2*(epsilon60+epsilon120)-epsilon0)/3
+disp(epsilony,"value of epsilony is= ")
+gammaxy=(2/sqrt(3))*(epsilon60-epsilon120)// from eqn 1 and eqn 2
+disp(gammaxy,"value of gammaxy is= ")
+epsilon1=((epsilonx+epsilony)/2)+sqrt(((epsilonx-epsilony)/2)^2+(gammaxy/2)^2)// epsilony value is in negative so the sign changes in the eqn
+disp(epsilon1,"value of epsilon1 is= ")
+epsilon2=((epsilonx+epsilony)/2)-sqrt(((epsilonx-epsilony)/2)^2+(gammaxy/2)^2)//epsilony value is in negative so the sign changes in the eqn
+disp(epsilon2,"value of epsilon2 is= ")
+
+gammamax=(2*10^-6)*sqrt(((epsilonx-epsilony)/2)^2+(gammaxy/2)^2)
+disp(gammamax,"max shear strain is= ")
+thetap=atand(577/320)/2
+disp(thetap,"orientations of principal axes is= ") // or
+thetap1=atand(577/320)*2
+disp(thetap1,"orientations of principal axes is= ")
+sigma1=(200*10^9/(1-0.09))*(epsilon1+0.3*epsilon2)
+disp(sigma1,"plane stresss is Pa= ")
+sigma2=(200*10^3/(1-0.09))*(epsilon2+0.3*epsilon1)
+disp(sigma2,"plane stresss is MPa= ")
+
+taumax=(200*10^9/(2*(1+0.3)))*gammamax
+disp(taumax,"plane stresss is MPa= ")
+
+
-- 
cgit