1 files changed, 39 insertions, 0 deletions
diff --git a/3864/CH7/EX7.22/Ex7_22.sce b/3864/CH7/EX7.22/Ex7_22.sce
new file mode 100644
index 000000000..99d992598
--- /dev/null
+++ b/3864/CH7/EX7.22/Ex7_22.sce
@@ -0,0 +1,39 @@
+clear
+//
+//
+
+//Initilization of Variables
+
+//strains
+e_A=500 //microns
+e_B=250 //microns
+e_C=-150 //microns
+E=2*10**5 //N/mm**2 //Modulus of Elasticity
+mu=0.3 //Poissons ratio
+theta=45 //Degrees
+
+//Calculations
+e_x=500
+e_A=500
+e_45=250
+e_B=250
+e_y=-150 
+e_C=-150 
+
+//e_45=(e_x+e_y)*2**-1+(e_x-e_y)*2**-1*cos(2*theta)+rho_x_y*2**-1*sin(2*theta)
+//After sub values and further simplifying we get
+rho_x_y=(e_45-(e_x+e_y)*2**-1-(e_x-e_y)*2**-1*cos(2*theta*%pi*180**-1))*(sin(2*theta*%pi*180**-1))**-1*2
+
+//Principal strains are given by
+e1=(e_x+e_y)*2**-1+(((e_x-e_y)*2**-1)**2+(rho_x_y*2**-1)**2)**0.5 //microns
+e2=(e_x+e_y)*2**-1-(((e_x-e_y)*2**-1)**2+(rho_x_y*2**-1)**2)**0.5 //microns
+
+//Principal Stresses
+sigma1=E*(e1+mu*e2)*(1-mu**2)**-1*10**-6 //N/mm**2
+sigma2=E*(e2+mu*e1)*(1-mu**2)**-1*10**-6 //N/mm**2
+
+//Result
+printf("\n Principal Strains are:e1 %0.2f  N/mm**2",e1)
+printf("\n                      :e2 %0.2f  N/mm**2",e2)
+printf("\n Principal Stresses are:sigma1 %0.2f  N/mm**2",sigma1)
+printf("\n                       :sigma2 %0.2f  N/mm**2",sigma2)