diff options
Diffstat (limited to '3750/CH1/EX1.11/Ex1_11.sce')
| -rw-r--r-- | 3750/CH1/EX1.11/Ex1_11.sce | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/3750/CH1/EX1.11/Ex1_11.sce b/3750/CH1/EX1.11/Ex1_11.sce new file mode 100644 index 000000000..f2d0a5dc6 --- /dev/null +++ b/3750/CH1/EX1.11/Ex1_11.sce @@ -0,0 +1,28 @@ +//Strength Of Material By G.H.Ryder
+//Chapter 1
+//Example 11
+// To Calculate thermal Stress in rod & tube
+SteelOD=2.4; //External diameter of steel tube, Unit in cm
+SteelID=1.8; //Internal diameter of steel tube, Unit in cm
+CopperDia=1.5; //Diameter of copper rod, unit in mm
+Es=210,000; //Young's Modulus for steel , Unit in N/mm^2
+Ec=100,00; //Young's Modulus for copper , Unit in N/mm^2
+alphaS=11e-6; //co-efficient of linear expansion for steel, Unit in perdegreeC
+alphaC=18e-6; //co-efficient of linear expansion for copper, Unit in perdegreeC
+AreaSteel=%pi*(SteelOD^2-SteelID^2)/4; //cross section Area of steel tube, Unit in cm^2
+AreaCopper=%pi*CopperDia^2/4; //cross section Area of copper bar, Unit in cm^2
+
+//Equillibrium Equation : SigmaC*AreaCopper=SigmaS*AreaSteel
+Ti=10; //Initial Temperature, Unit in perdegreeC
+Tf=200;//Final Temperature, Unit in perdegreeC
+//Compatibility Equation:
+ //alphaC*(Tf-Ti)-SigmaC/Ec=alphaS*(Tf-Ti)-SigmaS/Es
+SigmaS=(alphaS+alphaC)*(Tf-Ti)/((1/Es)*AreaSteel/(AreaCopper*Ec));
+//Stress in steel, Expression from compatability & Equillibrium Equation, Unit in N/mm^2
+
+SigmaC=AreaSteel*SigmaS/AreaCopper;
+printf("Stress in Copper rod= %f N/mm^2\n",SigmaC)
+printf("Stress in Steel Tube= %f N/mm^2\n",SigmaS)
+
+
+
|