diff options
Diffstat (limited to '3776/CH7')
| -rw-r--r-- | 3776/CH7/EX7.1/Ex7_1.sce | 23 | ||||
| -rw-r--r-- | 3776/CH7/EX7.2/Ex7_2.sce | 23 | ||||
| -rw-r--r-- | 3776/CH7/EX7.6/Ex7_6.sce | 15 | ||||
| -rw-r--r-- | 3776/CH7/EX7.8/Ex7_8.sce | 18 | ||||
| -rw-r--r-- | 3776/CH7/EX7.9/Ex7_9.sce | 23 |
5 files changed, 102 insertions, 0 deletions
diff --git a/3776/CH7/EX7.1/Ex7_1.sce b/3776/CH7/EX7.1/Ex7_1.sce new file mode 100644 index 000000000..22b3a758c --- /dev/null +++ b/3776/CH7/EX7.1/Ex7_1.sce @@ -0,0 +1,23 @@ +clear +//Given +shear_v = 3000 //N - Transmitted vetical shear +shear_al = 700 //N - The maximum allowable +//We will divide this into two parts +l_1 = 50.0 //mm +l_2 = 200.0 //mm +b_1 = 200.0 //mm +b_2 = 50.0 //mm +A_1 = l_1* b_1 //mm2 - area of part_1 +y_1 = 25.0 //mm com distance +A_2 =l_2*b_2 //mm2 - area of part_1 +y_2 = 150.0 //in com distance +y_net = (A_1*y_1 +A_2*y_2)/(A_1+A_2) //mm - The com of the whole system +c_max = (4-y_net) //mm - The maximum distace from com to end +c_min = y_net //mm - the minimum distance from com to end +I_1 = b_1*(l_1**3)/12 + A_1*((y_1-y_net)**2) //Parallel axis theorm +I_2 = b_2*(l_2**3)/12 + A_2*((y_2-y_net)**2) +I_net = I_1 + I_2 //mm4 - the total moment of inertia +Q = A_1*(-y_1+y_net) //mm3 +q = shear_v*Q/I_net //N/mm - Shear flow +d = shear_al/q // The space between the nails +printf("\n The minimal space between the nails %0.1f mm",d) diff --git a/3776/CH7/EX7.2/Ex7_2.sce b/3776/CH7/EX7.2/Ex7_2.sce new file mode 100644 index 000000000..71c8e8755 --- /dev/null +++ b/3776/CH7/EX7.2/Ex7_2.sce @@ -0,0 +1,23 @@ +clear +//Given +l = 6 //m -length of the beam +p = 3 //KN-m _ the load applied +R_a = l*p/2 //KN -The reaction at a, Since the system is symmetry +R_b = l*p/2 //KN -The reaction at b +l_s = 10 //mm - The length of the screw +shear_al = 2 //KN - The maximum load the screw can take +I = 2.36*(10**9) //mm2 The moment of inertia of the whole system +//We will divide this into two parts +l_1 = 50.0 //mm +l_2 = 50.0 //mm +b_1 = 100.0 //mm +b_2 = 200.0 //mm +A_1 = l_1* b_1 //in2 - area of part_1 +y_1 = 200.0 //mm com distance +A_2 =l_2*b_2 //mm2 - area of part_1 +y_2 = 225.0 //in com distance +Q = 2*A_1*y_1 + A_2*y_2 // mm**3 For the whole system +q = R_a*Q*(10**3)/I //N/mm The shear flow +d = shear_al*(10**3)/q //mm The space between the nails +printf("\n The minimal space between the nails %0.0f mm",d) +printf("\n Similar calculation for 4.5 KN gives spacing of 246mm") diff --git a/3776/CH7/EX7.6/Ex7_6.sce b/3776/CH7/EX7.6/Ex7_6.sce new file mode 100644 index 000000000..89fca203f --- /dev/null +++ b/3776/CH7/EX7.6/Ex7_6.sce @@ -0,0 +1,15 @@ +clear +//Given +//we will divide this into two equal parts and other part +l = 10.0 // in - The height +t = 0.1 // in - The width +b = 5.0 //mm- The width of the above part +A = t* b //in2 - area of part +y_net = l/2 // The com of the system +y_1 = l // The position of teh com of part_2 +I_1 = t*(l**3)/12 //in^4 The moment of inertia of part 1 +I_2 = 2*A*((y_1-y_net)**2) //in^4 The moment of inertia of part 2 +I = I_1 + I_2 //in^4 The moment of inertia +e = (b**2)*(l**2)*t/(4*I) //in the formula of channels +l_sc = e - t/2 //in- The shear centre +printf("\n The shear centre from outside vertical face is %0.3f in",l_sc ) diff --git a/3776/CH7/EX7.8/Ex7_8.sce b/3776/CH7/EX7.8/Ex7_8.sce new file mode 100644 index 000000000..67646ba58 --- /dev/null +++ b/3776/CH7/EX7.8/Ex7_8.sce @@ -0,0 +1,18 @@ +clear +//Given +dia = 10.0 //mm - The diameter of the cylinder +c = dia/2 //mm - the radius of the cylinder +A = 3.14*(c**2) //mm2 The area of the crossection +y = 4*c/(3*3.14) //mm The com of cylinder +I = 3.14*(c**4)/4 //mm4 - The moment of inertia of the cylinder +j = 3.14*(dia**4)/32 //mm4 +T = 20.0 //N.m - The torque +V = 250.0 //N - The shear +M = 25.0 //N-m The bending moment +Q = A*y/2 //mm +stress_dmax = 4*V/(3*A) //V*Q/(I*d) //MPa The direct maximum stress +stress_tmax = T*c*(10**3)/j //-MPa The torsion maximum stress +stress_total = stress_dmax + stress_tmax //MPa The total stress +printf("\n The direct maximum stress %0.2f MPa",stress_dmax) +printf("\n The torsion maximum stress %0.2f MPa",stress_tmax) +printf("\n The total stress %0.2f MPa",stress_total) diff --git a/3776/CH7/EX7.9/Ex7_9.sce b/3776/CH7/EX7.9/Ex7_9.sce new file mode 100644 index 000000000..4822c44e0 --- /dev/null +++ b/3776/CH7/EX7.9/Ex7_9.sce @@ -0,0 +1,23 @@ +clear +//Given +dia = 15 //mm - The diameter of the rod +h = 0.5 //mt - The freely falling height +A = 3.14*(dia**2)/4 //mm2 The area of the crossection +E = 200 //Gpa -Youngs modulus +L = 750 //mm - The total length of the rod +G = 80 //gpa - Shear modulus +N = 10 //number of live coils +d = 5 //mm the diameter of live coil +m = 3 // the mass of freely falling body +H = 500 //mm -from mass to spring +F= m*9.81 //Kg the force due to that mass +p = 3 //KN-m _ the load applied +//e = e_rod + e_spr +//e_rod +e_rod = p*L*(10**-3)/(A*E) //mm The elongation due to freely falling body +//e_spr +e_spr = 64*F*(dia**3)*N*(10**-3)/(G*(d**4)) //mm The elongation due to spring +e = e_rod + e_spr //mm The total elongation +p_dyn =F*(1+((1+(2*H/e))**0.5)) +Stress_max = p_dyn/A //MPa - The maximum stress in the system +printf("\n The maximum stress in the system %0.2f MPa",Stress_max) |