Added new codeHEADmaster

1 files changed, 28 insertions, 0 deletions

diff --git a/3863/CH1/EX1.31/Ex1_31.sce b/3863/CH1/EX1.31/Ex1_31.sce
new file mode 100644
index 000000000..395352333
--- /dev/null
+++ b/3863/CH1/EX1.31/Ex1_31.sce
@@ -0,0 +1,28 @@
+clear
+//
+
+//Given
+//Variable declaration
+Dc=15             //Diameter of copper rod in mm
+Di_s=20           //Internal diameter of steel in mm
+De_s=30           //External diameter of steel in mm
+T1=10             //Initial temperature in degree celsius
+T2=200            //Raised temperature in degree celsius
+Es=2.1e5          //Youngs modulus of steel in N/sq.mm
+Ec=1e5            //Youngs modulus of copper in N/sq.mm
+alpha_s=11e-6     //Co-efficient of linear expansion of steel in per degree celsius
+alpha_c=18e-6     //Co-efficient of linear expansion of copper in per degree celsius
+
+//Calculation
+Ac=(%pi/4)*Dc**2                //Area of copper tube in sq.mm
+As=(%pi/4)*(De_s**2-Di_s**2)    //Area of steel rod in sq.mm
+T=T2-T1                             //Rise of temperature in degree celsius
+sigmas=(((alpha_c-alpha_s)*T)/(((As/Ac)/Ec)+(1/Es)))
+
+sigmac=(sigmas*(As/Ac))
+
+
+//Result
+printf("\n NOTE: The answers in the book for stresses are wrong.The correct answers are,")
+printf("\n Stress in steel = %0.3f  N/mm^2",sigmas)
+printf("\n Stress in copper = %0.3f  N/mm^2",sigmac)