diff options
Diffstat (limited to '3776/CH11/EX11.11/Ex11_11.sce')
| -rw-r--r-- | 3776/CH11/EX11.11/Ex11_11.sce | 91 |
1 files changed, 40 insertions, 51 deletions
diff --git a/3776/CH11/EX11.11/Ex11_11.sce b/3776/CH11/EX11.11/Ex11_11.sce index 510361b42..ba17c5e84 100644 --- a/3776/CH11/EX11.11/Ex11_11.sce +++ b/3776/CH11/EX11.11/Ex11_11.sce @@ -1,59 +1,48 @@ clear // -P = 200.0 //K The force on the beam +P = 200.0 //K The force on the beam L = 15 //ft - The length of the rod -F_y = 50.0 //ksi +F_y = 50.0 //ksi F_a = F_y/(5.0/3) //ksi -AISC MANUAL ,allowable axial stress if axial force is alone F_b = F_a //Allowable compressive bending stress M_1 = 600.0 //k-in - The moment acting on the ends of the rod M_2 = 800.0 //k-in - the moment acting on the other end of teh rod -B_x = 0.264 //in - Extracted from AISC manual -E = 29*(10**3) -A = P/F_a + M_2*B_x/F_b //in2- The minimum area +B_x = 0.264 //in - Extracted from AISC manual +E = 29*(10**3) +A = P/F_a + M_2*B_x/F_b //sq.in- The minimum area printf("\n \n The minimum area is %0.2f in^2",A) -//we will select W10x49 section -A_s = 14.4 //in2 - The area of the section -r_min = 2.54 //in The minimum radius -r_x = 4.35 //in -f_a = P/A_s //Ksi- The computed axial stress -f_b = M_2*B_x/A_s //Computed bending stess -C_c = ((2*(%pi**2)*E/F_y)**0.5) //Slenderness ratio L/R -C_s = L*12/r_min // Slenderness ratio L/R of the present situation -if C_s <C_c then - printf("\n The following approch is solvable") -else - printf("\n The caliculation is not possible") - end -F_a_1 = 19.3 //Ksi - AISC lets try this -c_m = 0.6 - 0.4*(-M_1/M_2) -F_e = (12*(%pi**2)*E)/(23*(L*12/r_x)**2) -k = f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b //Condition mentioned in AISC -if k>1 then - printf("\n The following W10x49 section is not satisfying our constraints since f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b %0.3f >1",k) -else - printf("\n The following W10x49 section is satisfying our constraints since f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b %0.3f <1",k) - end -//trail - 2 -//Lets take W10 x 60 -A_s = 17.6 //in2 - The area of the section -r_min = 2.57 //in The minimum radius -r_x = 4.39 //in -f_a = P/A_s //Ksi- The computed axial stress -f_b = M_2*B_x/A_s //Computed bending stess -C_c = ((2*(%pi**2)*E/F_y)**0.5) //Slenderness ratio L/R -C_s = L*12/r_min // Slenderness ratio L/R of the present situation -if C_s <C_c then - printf("\n The following approch is solvable") -else - printf("\n The caliculation is not possible") - end -F_a_1 = 19.3 //Ksi - AISC lets try this -c_m = 0.6 - 0.4*(-M_1/M_2) -F_e = (12*(%pi**2)*E)/(23*(L*12/r_x)**2) -k = f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b //Condition mentioned in AISC -if k>1 then - printf("\n The following W10x49 section is not satisfying our constraints since f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b %0.3f >1",k) -else - printf("\n The following W10x49 section is satisfying our constraints since f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b %0.3f <1",k) - end -printf("\n small variation due to rounding off errors")
\ No newline at end of file + + + + +for i=1:2 + + st =['W10x49', 'W10x60'] + printf("\n we will select %s section",st(i)) + + A_s = [14.4, 17.6 ] //sq.in - The area of the section + r_min = [2.54 , 2.57 ] //in The minimum radius + r_x = [4.35 ,4.39] //in + f_a = P/A_s(i) //Ksi- The computed axial stress + f_b = M_2*B_x/A_s(i) //Computed bending stess + C_c = ((2*(%pi**2)*E/F_y)**0.5) //Slenderness ratio L/R + C_s = L*12/r_min(i) // Slenderness ratio L/R of the present situation + if C_s <C_c then + printf("\n Since calculated Le/r ratio is less than Cc, we can apply the second ASD formula") + else + printf("\n The calculation is not possible") + end + F_a_1 = 19.3 //Ksi - AISC lets try this + c_m = 0.6 - 0.4*(-M_1/M_2) + F_e = (12*(%pi**2)*E)/(23*(L*12/r_x(i))**2) + k = f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b //Condition mentioned in AISC + if k>1 then + printf("\n The following %s section is not satisfying our constraints since f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b %0.3f >1",st(i),k) + else + printf("\n The following %s section is satisfying our constraints since f_a/F_a_1 + c_m*f_b*(1-(f_a/F_e))/F_b %0.3f <1",st(i),k) + end + +end + + +printf("\n small variation due to rounding off errors") |