From d1e070fe2d77c8e7f6ba4b0c57b1b42e26349059 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Sat, 3 Feb 2018 10:59:42 +0530
Subject: Modified the code

---
 3776/CH11/EX11.6/Ex11_6.sce | 48 +++++++++++++++++++++++++++++++--------------
 1 file changed, 33 insertions(+), 15 deletions(-)

(limited to '3776/CH11/EX11.6/Ex11_6.sce')

diff --git a/3776/CH11/EX11.6/Ex11_6.sce b/3776/CH11/EX11.6/Ex11_6.sce
index a095bf461..a2b8da993 100644
--- a/3776/CH11/EX11.6/Ex11_6.sce
+++ b/3776/CH11/EX11.6/Ex11_6.sce
@@ -2,36 +2,54 @@ clear
 //Given
 //
 L = 15                                    //ft - The length of the each rod
-A = 46.7                                  //in2 - The length of the crossection 
+A = 46.7                                  //sq.in - The length of the crossection
 r_min = 4                                 //in - The radius of gyration
 stress_yp = 36                            //ksi - the yielding point stress
 E = 29*(10**3)                            //ksi - The youngs modulus
 C_c = ((2*(%pi**2)*E/stress_yp)**0.5) //Slenderness ratio L/R
-C_s = L*12/r_min                          // Slenderness ratio L/R of the present situation 
-//According to AISC formulas 
-if (C_s <C_c) then
-    printf ("a)The following approch is solvable")
+C_s1 = L*12/r_min                          // Slenderness ratio L/R of the present situation
+//According to AISC formulas
+printf ("a)calculated Le/r ratio is %f",C_s1)
+if (C_s1 <C_c) then
+    printf ("\n a)Since calculated Le/r ratio is less than Cc(126), we can apply the second ASD formula")
 else
-    print ("The caliculation is not possible")
+    print ("The calculation is not possible")
     end
-F_S = 5.0/3 +3*C_s/(8*C_c) -(3*C_s**3)/(8*C_c**3)        //Safety factor 
-Stress_all = (1 - (C_s**2)/(2*C_c**2))*stress_yp/F_S    //The allowable strees 
+F_S = 5.0/3 +3*C_s1/(8*C_c) -(C_s1**3)/(8*C_c**3)        //Safety factor
+Stress_all = (1 - (C_s1**2)/(2*C_c**2))*stress_yp/F_S    //The allowable stress
 printf("\n a) The allowable stress in this case is %0.2f kips",Stress_all)
+printf("\n a) The allowable pressure in this case is %0.2f kips",Stress_all*A)
+
 //Part - B
 //Given
 L = 40                                   //ft - The length of the each rod
-A = 46.7                                  //in2 - The length of the crossection 
+A = 46.7                                  //sq.in - The length of the crossection
 r_min = 4                                 //in - The radius of gyration
 stress_yp = 36                            //ksi - the yielding point stress
 E = 29*(10**3)                            //ksi - The youngs modulus
 C_c = ((2*(%pi**2)*E/stress_yp)**0.5) //Slenderness ratio L/R
-C_s = L*12/r_min                          // Slenderness ratio L/R of the present situation 
-//According to AISC formulas 
+C_s = L*12/r_min                          // Slenderness ratio L/R of the present situation
+//According to AISC formulas
+printf ("\n b)calculated Le/r ratio is %f",C_s)
 if C_s <C_c then
-    printf("b) The following approch is solvable")
+    printf("\n b)Since calculated Le/r ratio is less than Cc,  we can apply the second ASD formula")
 else
-    printf("The caliculation is not possible")
+    printf("The calculation is not possible")
 end
-F_S = 5.0/3 +3*C_s/(8*C_c) -(3*C_s**3)/(8*C_c**3)        //Safety factor 
-Stress_all = (1 - (C_s**2)/(2*C_c**2))*stress_yp/F_S    //The allowable strees 
+F_S = 5.0/3 +3*C_s/(8*C_c) -(C_s**3)/(8*C_c**3)        //Safety factor
+Stress_all = (1 - (C_s**2)/(2*C_c**2))*stress_yp/F_S    //The allowable stress
 printf("\n b) The allowable stress in this case is %0.2f kips",Stress_all)
+printf("\n b) The allowable pressure in this case is %0.2f kips",Stress_all*A)
+printf("\n Similarly, for a column fixed at one end")
+C_s = L*12/r_min*.8                          // Slenderness ratio L/R of the present situation
+//According to AISC formulas
+printf ("\n b)calculated Le/r ratio is %f",C_s)
+if C_s <C_c then
+    printf("\n b)Since calculated Le/r ratio is less than Cc,  we can apply the second ASD formula")
+else
+    printf("The calculation is not possible")
+end
+F_S = 5.0/3 +3*C_s/(8*C_c) -(C_s**3)/(8*C_c**3)        //Safety factor
+Stress_all = (1 - (C_s**2)/(2*C_c**2))*stress_yp/F_S    //The allowable stress
+printf("\n b) The allowable stress in this case is %0.2f kips",Stress_all)
+printf("\n b) The allowable pressure in this case is %0.0f kips",Stress_all*A)
-- 
cgit