summaryrefslogtreecommitdiff
path: root/3776/CH11/EX11.11/Ex11_11.sce
blob: 510361b4225f6754d87e189efb2a60cb809b290e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
clear
//
P = 200.0       //K The force on the beam 
L = 15          //ft - The length of the rod
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 
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")