1 files changed, 55 insertions, 0 deletions
diff --git a/3864/CH6/EX6.15/Ex6_15.sce b/3864/CH6/EX6.15/Ex6_15.sce
new file mode 100644
index 000000000..a96281897
--- /dev/null
+++ b/3864/CH6/EX6.15/Ex6_15.sce
@@ -0,0 +1,55 @@
+clear
+//
+//
+
+//Initilization of Variables
+
+G=80*10**3 //N/mm**2 //Modulus of Rigidity
+d1=100 //mm //Outer diameter of hollow shft
+d2=80  //mm //Inner diameter of hollow shaft
+d=80   //mm //diameter of Solid shaft
+d3=60  //mm //diameter of Solid shaft having L=0.5m
+L1=300 //mm //Length of Hollow shaft
+L2=400 //mm //Length of solid shaft
+L3=500 //mm //LEngth of solid shaft of diameter 60mm
+T1=2*10**6 //N-mm //Torsion in Shaft AB
+T2=1*10**6 //N-mm //Torsion in shaft BC
+T3=1*10**6 //N-mm //Torsion in shaft CD
+
+//Calculations
+
+//Now Polar modulus of section AB
+J1=%pi*32**-1*(d1**4-d2**4) //mm**4 
+
+//Polar modulus of section BC
+J2=%pi*32**-1*d**4 //mm**4
+
+//Polar modulus of section CD
+J3=%pi*32**-1*d3**4 //mm**4
+
+//Now angle of twist of AB
+theta1=T1*L1*(G*J1)**-1 //radians
+
+//Angle of twist of BC
+theta2=T2*L2*(G*J2)**-1 //radians
+
+//Angle of twist of CD
+theta3=T3*L3*(G*J3)**-1 //radians
+
+//Angle of twist
+theta=theta1-theta2+theta3 //Radians
+
+//Shear stress in AB From Torsion Equation
+q_s1=T1*(d1*2**-1)*J1**-1 //N/mm**2 
+
+//Shear stress in BC
+q_s2=T2*(d*2**-1)*J2**-1 //N/mm**2 
+
+//Shear stress in CD
+q_s3=T3*(d3*2**-1)*J3**-1 //N-mm**2
+
+//As max shear stress occurs in portion CD,so consider CD
+
+//Result
+printf("\n Angle of twist at free end is %0.5f  Radian",theta)
+printf("\n Max Shear stress %0.2f  N/mm**2",q_s3)