4 files changed, 127 insertions, 0 deletions
diff --git a/3705/CH10/EX10.1/Ex10_1.sce b/3705/CH10/EX10.1/Ex10_1.sce
new file mode 100644
index 000000000..d2637d4a7
--- /dev/null
+++ b/3705/CH10/EX10.1/Ex10_1.sce
@@ -0,0 +1,28 @@
+
+clear//
+
+//Variable Declaration
+Le=7 //Effective length in m
+P=450 //Applied axial Load in kN
+FOS=3 //Factor of safety 
+sigma_pl=200*10**6 //Stress allowable in Pa
+E=200*10**9 //Youngs Modulus in Pa
+end_cond=0.7 //End Condition factor to be multiplied
+
+//Calculations
+Pcr=P*FOS //Critical Load in kN
+A=Pcr*sigma_pl**-1*10**9 //Area in mm^2
+
+//Part 1
+I1=10**15*(Pcr*Le**2)*(%pi**2*E)**-1 //Moment of Inertia Required in mm^4
+//From table selecting appropriate Section W250x73
+
+//Part 2
+I2=10**15*(Pcr*end_cond**2*Le**2)*(%pi**2*E)**-1 //Moment of Inertia Required in mm^4
+//From table selecting appropriate Section W200x52
+
+//Lightest Section that meets these criterion is W250x58 section
+
+
+//Result
+printf("\n From the above computation we select W250x58 section")
diff --git a/3705/CH10/EX10.2/Ex10_2.sce b/3705/CH10/EX10.2/Ex10_2.sce
new file mode 100644
index 000000000..a320e3bbe
--- /dev/null
+++ b/3705/CH10/EX10.2/Ex10_2.sce
@@ -0,0 +1,39 @@
+
+clear//
+
+//Variable Declaration
+E=200*10**9 //Youngs Modulus in Pa
+sigma_yp=380*10**6 //Stress allowable in Pa
+Le=10 //Length in m
+end_cond=0.5 //Support condition factor to be ,ultiplied to length
+A=15.5*10**-3 //Area in m^2
+
+//Calculations
+Cc=sqrt((2*%pi**2*E)*sigma_yp**-1) //Slenderness Ratio
+
+//Part 1
+S_R1=142.9 //Slenderness ratio 
+sigma_w=(12*%pi**2*E)/(23*S_R1**2) //Allowable Working Stress in Pa
+P=sigma_w*A //Maximum Allowable Load in kN
+
+//Part 2
+S_R2=79.37 //Slenderness ratio 
+N=5*3**-1+((3*S_R2)/(8*Cc))-(S_R2**3*(8*Cc**3)**-1) //Factor Of Safety
+
+sigma_w2=(1-(S_R2**2*0.5*Cc**-2))*(sigma_yp*N**-1) //Allowable working Stress in Pa
+P2=sigma_w2*A //Allowable Load in kN
+
+//Part 3
+S_R3=55.56 //Slenderness Ratio
+N3=5*3**-1+((3*S_R3)/(8*Cc))-(S_R3**3*(8*Cc**3)**-1) //Factor Of Safety
+
+sigma_w3=(1-(S_R3**2*0.5*Cc**-2))*(sigma_yp*N3**-1) //Allowable working Stress in Pa
+P3=sigma_w3*A //Allowable Load in kN
+
+//Result
+printf("\n The results for Part 1 are")
+printf("\n Maximum Allowable Load P= %0.0f kN",P*10**-3)
+printf("\n Part 2")
+printf("\n Maximum Allowable Load P= %0.0f kN",P2*10**-3)
+printf("\n Part 3")
+printf("\n Maximum Allowable Load P= %0.0f kN",P3*10**-3)
diff --git a/3705/CH10/EX10.3/Ex10_3.sce b/3705/CH10/EX10.3/Ex10_3.sce
new file mode 100644
index 000000000..2f4065870
--- /dev/null
+++ b/3705/CH10/EX10.3/Ex10_3.sce
@@ -0,0 +1,31 @@
+
+clear//
+//
+
+//Variable Declaration
+E=29*10**6 //Youngs Modulus in psi
+sigma_yp=36*10**3 //Stress in psi
+L=25 //Length in ft
+A=17.9 //Area in in^2
+Iz=640 //Moment of inertia in in^4
+Sz=92.2 //Sectional Modulus in in^3
+P=150*10**3 //Load in lb
+e=4 //Eccentricity in inches
+
+//Calculations
+
+//Part 1
+a=1.09836
+sigma_max=P*A**-1+(P*e*Sz**-1)*a //Maximum Stress in psi
+
+//Part 2
+//After simplification we get the equation to compute N
+N=2.19 //Trial and Error yields
+
+//Part 3
+v_max=e*((cos(sqrt((P*L**2*12**2)*(4*E*Iz)**-1)))**-1-1)
+
+//Result
+printf("\n The maximum compressive stress in the Column is %0.2f psi",sigma_max)
+printf("\n The factor of safety is %0.3f ",N)
+printf("\n The maximum lateral dfelection is %0.3f in",v_max)
diff --git a/3705/CH10/EX10.4/Ex10_4.sce b/3705/CH10/EX10.4/Ex10_4.sce
new file mode 100644
index 000000000..9a1ca8516
--- /dev/null
+++ b/3705/CH10/EX10.4/Ex10_4.sce
@@ -0,0 +1,29 @@
+
+clear//
+
+//Variable Declaration
+Le=7 //Effective Length in m
+N=2 //Factor of Safety
+h_max=400 //Maximum depth in mm
+E=200 //Youngs Modulus in GPa
+sigma_yp=250 //Maximum stress in yielding in MPa
+P1=400 //Load 1 in kN
+P2=900 //Load 2 in kN
+x1=75 //Distance in mm
+x2=125 //Distance in mm
+
+//Calculations
+e=(P2*x2-P1*x1)*(P1+P2)**-1 //Eccentricity in mm
+P=N*(P1+P2) //Applied load after factor of safety is considered in kN
+
+//Part 1 is not computable
+I=415*10**-6 //Moment of inertia from the table in mm^4
+
+//Part 2
+P_cr=%pi**2*E*10**9*I*Le**-2 //Critical load for buckling in kN
+FOS=P_cr*10**-3/(P1+P2) //Factor of safety against buckling in y-axis
+
+
+//Result
+printf("\n The critical load for buckling is %0.0f kN",P_cr*10**-3)
+printf("\n The factor of safety is %0.1f ",FOS)