3 files changed, 84 insertions, 0 deletions
diff --git a/3705/CH9/EX9.1/Ex9_1.sce b/3705/CH9/EX9.1/Ex9_1.sce
new file mode 100644
index 000000000..07b55e8aa
--- /dev/null
+++ b/3705/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,33 @@
+
+clear//
+
+//Variable Declaration
+n=20 //Modular Ratio
+sigma_wd=8*10**6 //Maximum bending stress in wood in Pa
+sigma_st=120*10**6 //Maximum bending stress in steel in Pa
+
+//Cross Sectional Details
+Awd=45 //Area of wood in mm^2
+y_wd=160 //Neutral Axis of from bottom of the wooden section in mm
+Ast=15 //Area of steel in mm^2
+y_st=5 //Neutral Axis of the Steel section in mm
+//Dimensions
+ww=150 //width of wooden section in mm
+dw=300 //depth of wooden section in mm
+ws=75 //width of steel section in mm
+ds=10 //depth of steel section in mm
+
+//Calculations
+y_bar=(Awd*y_wd+Ast*y_st)*(Ast+Awd)**-1 //Location of Neutral axis from the bottom in mm
+//Moment of inertia 
+I=(ww*dw**3*12**-1)+(ww*dw*(y_wd-y_bar)**2)+(n*ws*ds**3*12**-1)+(n*ws*ds*(y_bar-y_st)**2) //mm^4
+c_top=dw+ds-y_bar //Distance from NA to top fibre in mm
+c_bot=y_bar //Distance from NA to bottom fibre in mm
+
+//The solution will be in different order 
+M1=sigma_wd*I*10**-12*c_top**-1 //Maximum Bending Moment in N.m
+M2=sigma_st*I*10**-12*c_bot**-1 //Maximum Bending Moment in N.m
+M=min(M1,M2) //Maximum allowable moment in N.m
+
+//Result
+printf("\n The Maximum Allowable moment that the beam can support is %0.1f kN.m",M)
diff --git a/3705/CH9/EX9.2/Ex9_2.sce b/3705/CH9/EX9.2/Ex9_2.sce
new file mode 100644
index 000000000..b90500567
--- /dev/null
+++ b/3705/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,28 @@
+
+clear//
+
+//Variable Declaration
+dw=8 //Depth of wooden section in inches
+da=0.4 //Depth og aluminium section in inches 
+w=2 //Width of the section in inches
+n=40*3**-1 //Modular Ratio
+Ewd=1.5*10**6 //Youngs modulus of wood in psi
+Eal=10**7 //Youngs Modulus of aluminium in psi
+V_max=4000 //Maximum shear in lb
+b=24 //Inches
+L=72 //Length in inches
+P=6000 //Load on the beam in lb
+
+//Calculations
+I=w*dw**3*12**-1+2*(n*w*da**3*12**-1+n*da*4.2**2) //Moment of Inertia in in^4
+
+//Part 1
+Q=(w*dw*0.5)*2+(n*da)*(dw*0.5+da*0.5) //First Moment in in^3
+tau_max=V_max*Q*I**-1*w**-1 //Maximum Shear Stress in psi
+
+//Part 2
+delta_mid=(P*b)*(48*Ewd*I)**-1*(3*L**2-4*b**2)
+
+//Result
+printf("\n The maximum shear stress allowable is %0.0f psi",tau_max)
+printf("\n The deflection at the mid-span is %0.4f in",delta_mid)
diff --git a/3705/CH9/EX9.4/Ex9_4.sce b/3705/CH9/EX9.4/Ex9_4.sce
new file mode 100644
index 000000000..99ca2c8ef
--- /dev/null
+++ b/3705/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,23 @@
+
+clear//
+
+//Variable Declaration
+sigma_co_w=12 //Maximum stress in compression in MPa
+sigma_st_w=140 //Maximum stress in steel in MPa
+M=90 //Moment in kN.m
+n=8 //Modular Ratio 
+
+//Calculations
+//h=0.4068d
+//bd^2=0.04266
+b=(0.04266/(1.5**2))**0.3333 //Breadth in m 
+d=1.5*b //Depth in m
+h=0.4068*d //Height in m
+
+//Area of steel
+Ast=((M*10**3)/((d-h*3**-1)*sigma_st_w*10**3))*10**3 //Area of steel in mm^2
+
+//Result
+printf("\n The dimensions of the beam are")
+printf("\n b= %0.0f mm and d= %0.0f mm",b*1000,d*1000)
+printf("\n Area of steel= %0.0f mm^2",Ast)