1 files changed, 27 insertions, 0 deletions
diff --git a/3289/CH4/EX4.5/Ex4_5.sce b/3289/CH4/EX4.5/Ex4_5.sce
new file mode 100755
index 000000000..f785f9700
--- /dev/null
+++ b/3289/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,27 @@
+clc
+
+r=0.04 //m
+t=5 //mm
+sigmae=250 //MPa
+sigmay=300 //MPa
+
+//sigmathetamax=(p*r)/t =8*p        max values of tangential stresses
+//sigmathetamin=((-p/4)*r)/t =-2*p   min values of tangential stresses
+//sigmazamax=(p*r)/2*t =4*p       axial principal stresses
+//sigmazmin=((-p/4)*r)/2*t =-p
+
+//sigmathetaa=(sigmathetamax-sigmathetamin)/2= 5p    alternating and mean stresses
+//sigmathetam=(sigmathetamax+sigmathetamin)/2= 3p
+//sigmaza=(sigmazamax-sigmazmin)/2 =2.5p
+//sigmazm=(sigmazamax+sigmazmin)/2 =1.5p
+
+//sqrt(sigmathetaa^2-sigmathetaa*sigmaza+sigmaza^2)=sigmaea
+//sqrt(sigmathetam^2-sigmathetam*sigmazm+sigmazm^2)=sigmaem
+
+//sqrt(25p^2-12.3p^2+6.25p^2)=sigmaea
+//sqrt(9p^2-4.5p^2+2.25p^2)=sigmaem           solving this equation we get
+sigmaea=4.33 //p
+sigmaem=2.60 //p
+
+p=1/((sigmaea/sigmae)+(sigmaem/sigmay))
+disp(p,"the value of p in MPa is= ")