initial commit / add all books

8 files changed, 154 insertions, 0 deletions

diff --git a/3776/CH4/EX4.12/Ex4_12.sce b/3776/CH4/EX4.12/Ex4_12.sce
new file mode 100644
index 000000000..05e9b7ae8
--- /dev/null
+++ b/3776/CH4/EX4.12/Ex4_12.sce
@@ -0,0 +1,10 @@
+clear
+//Given
+dai_bc =  240   //mm- daimeter of '8'bolt circle 
+dia =  dai_bc/8 //Diameter of each bolt
+A =  0.25*(dia**2)*3.14 // Area of a bolt
+S_allow  = 40           //MPa - The maximum allowable allowable shear stress 
+P_max =  (S_allow)*A    //N - The maximum allowable force 
+D = 120.0               //mm - the distance from central axis 
+T_allow =P_max*D*8      //N-m The allowable torsion on the 8 bolt combination 
+printf("\n The allowable torsion on the 8 bolt combination %e N-m",T_allow *10**-3)
diff --git a/3776/CH4/EX4.15/Ex4_15.sce b/3776/CH4/EX4.15/Ex4_15.sce
new file mode 100644
index 000000000..0d15b319f
--- /dev/null
+++ b/3776/CH4/EX4.15/Ex4_15.sce
@@ -0,0 +1,26 @@
+clear
+//Given 
+//AISC MANUALS
+//approximated by three narrow tubes 
+//J = Bbt^3
+B = 0.33 // constant mentiones in AISC
+//three rods 
+
+//rod_1
+t_1 = 0.605 //inch - Thickness 
+b =  12.0 //inches - width 
+J_1 = B*b*(t_1**3) //in^4 - Torsion constant 
+
+//rod_2
+t_2 = 0.605 //inch - Thickness 
+b =  12 //inches - width 
+J_2 = B*b*(t_2**3) //in^4 - Torsion constant 
+
+//rod_3
+t_3 = 0.390 //inch - Thickness 
+b =  10.91 //inches - width 
+J_3 = B*b*(t_3**3) //in^4 - Torsion constant 
+
+//Equivalent
+J_eq = J_1+J_2+J_3  //in^4 - Torsion constant 
+printf("\n the Equivalent Torsion constant is  %0.2f in^4 ",J_eq)
diff --git a/3776/CH4/EX4.16/Ex4_16.sce b/3776/CH4/EX4.16/Ex4_16.sce
new file mode 100644
index 000000000..d62fbcd29
--- /dev/null
+++ b/3776/CH4/EX4.16/Ex4_16.sce
@@ -0,0 +1,14 @@
+clear
+//Given 
+dia_out = 10       //mm- outer diameter of shaft
+dia_in  = 8     //mm- inner diameter of shaft 
+c_out = dia_out/2 //mm - outer Radius of shaft  
+c_in  = dia_in/2  //mm - inner radius of shaft 
+T = 40            //N/mm -Torque in the shaft 
+//caliculations
+
+J = 3.14*((dia_out**4)- (dia_in**4))/32 //mm4
+shear_T_max = T*c_out*(10**3)/J       // The maximum torsion shear in the shaft
+shear_T_min = T*c_in*(10**3)/J        // The maximum torsion shear in the shaft
+printf("\n The maximum shear due to torsion is  %0.2f MPa",shear_T_max)
+printf("\n The minimum shear due to torsion is  %0.0f MPa",shear_T_min)
diff --git a/3776/CH4/EX4.3/Ex4_3.sce b/3776/CH4/EX4.3/Ex4_3.sce
new file mode 100644
index 000000000..7ad10bbe6
--- /dev/null
+++ b/3776/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,14 @@
+clear
+//Given 
+dia_out = 20      //mm- outer diameter of shaft
+dia_in  = 16      //mm- inner diameter of shaft 
+c_out = dia_out/2 //mm - outer Radius of shaft  
+c_in  = dia_in/2  //mm - inner radius of shaft 
+T = 40            //N/mm -Torque in the shaft 
+//caliculations
+
+J = 3.14*((dia_out**4)- (dia_in**4))/32 //mm4
+shear_T_max = T*c_out*(10**3)/J       // The maximum torsion shear in the shaft
+shear_T_min = T*c_in*(10**3)/J        // The maximum torsion shear in the shaft
+printf("\n The maximum shear due to torsion is  %0.2f MPa",shear_T_max)
+printf("\n The minimum shear due to torsion is  %0.0f MPa",shear_T_min)
diff --git a/3776/CH4/EX4.4/Ex4_4.sce b/3776/CH4/EX4.4/Ex4_4.sce
new file mode 100644
index 000000000..105784f04
--- /dev/null
+++ b/3776/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,16 @@
+clear
+//Given
+hp = 10         // horse power of motor 
+f = 30          // given 
+shear_T = 55    //MPa - The maximum shearing in the shaft 
+//caliculations
+
+T = 119*hp/f            // N.m The torsion in the shaft 
+//j/c=T/shear_T=K
+k = T*(10**3)/shear_T //mm3
+//c3=2K/3.14
+c = ((2*k/3)**0.33)   //mm - The radius of the shaft 
+diamter = 2*c           //mm - The diameter of the shaft
+printf("\n The Diameter of the shaft used is %0.2f mm",diamter)
+printf("\n For practical purposes, a 16-mm shaft would probably be selected")
+
diff --git a/3776/CH4/EX4.5/Ex4_5.sce b/3776/CH4/EX4.5/Ex4_5.sce
new file mode 100644
index 000000000..8d033dc9a
--- /dev/null
+++ b/3776/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,23 @@
+clear
+//Given 
+hp = 200          //Horse power
+stress_sh = 10000 //psi- shear stress
+rpm_1 = 20.0      // The rpm at which this shaft1 operates 
+rpm_2 = 20000.0   // The rpm at which this shaft2 operates
+T_1= hp*63000.0/rpm_1 //in-lb Torsion due to rpm1
+T_2= hp*63000/rpm_2   //in-lb Torsion due to rpm1
+//caliculations 
+
+//j/c=T/shear_T=K
+k_1= T_1/stress_sh       //mm3
+//c3=2K/3.14
+c_1= ((2*k_1/3)**0.33) //mm - The radius of the shaft 
+diamter_1 = 2*c_1        //mm - The diameter of the shaft
+printf("\n The Diameter of the shaft1 is %0.2f mm",diamter_1)
+
+//j/c=T/shear_T=K
+k_2= T_2/stress_sh       //mm3
+//c3=2K/3.14
+c_2= ((2*k_2/3)**0.33) //mm - The radius of the shaft 
+diamter_2 = 2*c_2        //mm - The diameter of the shaft
+printf("\n The Diameter of the shaft2 is %0.3f mm",diamter_2)
diff --git a/3776/CH4/EX4.7/Ex4_7.sce b/3776/CH4/EX4.7/Ex4_7.sce
new file mode 100644
index 000000000..6b437ff50
--- /dev/null
+++ b/3776/CH4/EX4.7/Ex4_7.sce
@@ -0,0 +1,21 @@
+clear
+//Given 
+T_ab = 0    //N.m - torsion in AB 
+T_bc = 150  //N.m - torsion in BC
+T_cd = 150  //N.m - torsion in CD
+T_de = 1150 //N.m - torsion in DE
+l_ab = 250 //mm - length of AB
+l_bc = 200 //mm - length of BC
+l_cd = 300 //mm - length of cd 
+l_de = 500.0 //mm - length of de
+d_1 = 25 //mm - outer diameter 
+d_2 = 50 //mm - inner diameter
+G = 80 //Gpa -shear modulus
+//Caliculations 
+
+J_ab = 3.14*(d_1**4)/32           //mm4
+J_bc = 3.14*(d_1**4)/32           //mm4
+J_cd = 3.14*(d_2**4 - d_1**4)/32  //mm4
+J_de = 3.14*(d_2**4 - d_1**4)/32  //mm4
+rad =  T_ab*l_ab/(J_ab*G)+ T_bc*l_bc/(J_bc*G)+ T_cd*l_cd/(J_cd*G)+ T_de*l_de/(J_de*G) // adding the maximum radians roteted in each module
+printf("\n The maximum angle rotated is  %0.3f radians ",rad)
diff --git a/3776/CH4/EX4.9/Ex4_9.sce b/3776/CH4/EX4.9/Ex4_9.sce
new file mode 100644
index 000000000..5260a4497
--- /dev/null
+++ b/3776/CH4/EX4.9/Ex4_9.sce
@@ -0,0 +1,30 @@
+clear
+//given 
+//its a statistally indeterminant 
+//we will take of one of the support 
+//Given 
+T_ab = 0    //N.m - torsion in AB 
+T_bc = 150  //N.m - torsion in BC
+T_cd = 150  //N.m - torsion in CD
+T_de = 1150 //N.m - torsion in DE
+l_ab = 250  //mm - length of AB
+l_bc = 200  //mm - length of BC
+l_cd = 300  //mm - length of cd 
+l_de = 500.0//mm - length of de
+d_1 = 25 //mm - outer diameter 
+d_2 = 50 //mm - inner diameter
+//Caliculations 
+
+J_ab = 3.14*(d_1**4)/32          //mm4
+J_bc = 3.14*(d_1**4)/32          //mm4
+J_cd = 3.14*(d_2**4 - d_1**4)/32 //mm4
+J_de = 3.14*(d_2**4 - d_1**4)/32 //mm4
+G = 80 //Gpa -shear modulus
+rad =  T_ab*l_ab/(J_ab*G)+ T_bc*l_bc/(J_bc*G)+ T_cd*l_cd/(J_cd*G)+ T_de*l_de/(J_de*G) 
+//now lets consider T_A then the torsion is only T_A
+// T_A*(l_ab/(J_ab*G)+ l_bc/(J_bc*G)+ l_cd/(J_cd*G)+ l_de/(J_de*G)) +rad = 0
+// since there will be no displacement 
+T_A =-rad/(l_ab/(J_ab*G)+ l_bc/(J_bc*G)+ l_cd/(J_cd*G)+ l_de/(J_de*G)) //Torsion at A
+T_B = 1150 - T_A                                                        //n-m F_X = 0 torsion at B
+printf("\n The Torsion at rigid end A is %0.2f N-m",T_A)
+printf("\n The Torsion at rigid end B is %0.2f N-m",T_B)