Added new codeHEADmaster

8 files changed, 153 insertions, 0 deletions

diff --git a/3863/CH12/EX12.1/Ex12_1.sce b/3863/CH12/EX12.1/Ex12_1.sce
new file mode 100644
index 000000000..52e25e1a3
--- /dev/null
+++ b/3863/CH12/EX12.1/Ex12_1.sce
@@ -0,0 +1,20 @@
+clear
+//
+
+//Given
+//Variable declaration
+L=6*1000        //Length in mm
+W=50*1000       //Point load in N
+I=78e6          //Moment of Inertia in mm^4
+E=2.1e5         //Young's modulus in N/sq.mm
+
+//Calculation
+yc=((W*L**3)/(48*E*I))                      //The deflection at the centre in mm
+
+thetaB=((180/%pi)*((W*L**2)/(16*E*I)))    //The slope at the supports
+
+
+//Result
+printf("\n Deflection at the centre = %0.3f  mm",yc)
+printf("\n NOTE:The answer given for slope at the support is wrong.The correct answer is,")
+printf("\n Slope at the support = %0.3f  degree",thetaB)
diff --git a/3863/CH12/EX12.2/Ex12_2.sce b/3863/CH12/EX12.2/Ex12_2.sce
new file mode 100644
index 000000000..cdd152ddc
--- /dev/null
+++ b/3863/CH12/EX12.2/Ex12_2.sce
@@ -0,0 +1,14 @@
+clear
+//
+
+//Given
+//Variable declaration
+L=4*1000              //Length in mm
+
+//Calculation
+thetaA=((%pi/180)*(1))    //Slope at the ends in radians
+
+yc=(thetaA*(L/3))    //Deflection at the centre in mm
+
+//Result
+printf("\n Deflection at the centre = %0.3f  mm",yc)
diff --git a/3863/CH12/EX12.3/Ex12_3.sce b/3863/CH12/EX12.3/Ex12_3.sce
new file mode 100644
index 000000000..5accbe7e5
--- /dev/null
+++ b/3863/CH12/EX12.3/Ex12_3.sce
@@ -0,0 +1,14 @@
+clear
+//
+
+//Given
+//Variable declaration
+L=3*1000              //Length in mm
+
+//Calculation
+thetaA=((%pi/180)*(1))    //Slope at the ends in radians
+
+yc=(thetaA*(L/3))    //Deflection at the centre in mm
+
+//Result
+printf("\n Deflection at the centre = %0.3f  mm",yc)
diff --git a/3863/CH12/EX12.4/Ex12_4.sce b/3863/CH12/EX12.4/Ex12_4.sce
new file mode 100644
index 000000000..365ae9017
--- /dev/null
+++ b/3863/CH12/EX12.4/Ex12_4.sce
@@ -0,0 +1,25 @@
+clear
+//
+
+//Given
+//Variable declaration
+L=5*1000        //Length in mm
+W=5*1000        //Point load in N
+a=3*1000        //Distance between point load and left end in mm
+E=2e5           //Young's modulus in N/sq.mm
+I=1e8           //Moment of Inertia in mm^4
+
+//Calculation
+b=L-a           //Width in mm
+//case(i):The slope at the left support
+thetaA=-(W*a*b)/(6*E*I*L)*(a+2*b)
+//case(iii): The deflection under the load
+yc=(W*a**2*b**2)/(3*E*I*L)
+//case(iii):The maximum deflection
+y_max=((W*b)/(9*sqrt(3)*E*I*L)*(((a**2)+(2*a*b))**(3/2)))
+
+
+//Result
+printf("\n slope at the left support = %0.3f  radians",thetaA)
+printf("\n Deflection under the load = %0.3f  mm",yc)
+printf("\n Maximum deflection = %0.3f  mm",y_max)
diff --git a/3863/CH12/EX12.5/Ex12_5.sce b/3863/CH12/EX12.5/Ex12_5.sce
new file mode 100644
index 000000000..49f9db2c1
--- /dev/null
+++ b/3863/CH12/EX12.5/Ex12_5.sce
@@ -0,0 +1,27 @@
+clear
+//
+//
+
+//Given
+//Variable declaration
+b=200                //Width in mm
+d=300                //Depth in mm
+L=5*1000             //Span in mm
+L_star=5             //Span in m
+w=9*1000             //Uniformly distributed load in N/m
+E=1e4                //Youngs modulus in N/sq.mm
+
+//Calculation
+W=w*L_star                                 //Total load in N
+I=b*d**3/12                                //Moment of Inertia in mm^4
+
+//case(i):the slope at the support
+thetaA=(-(W*(L**2))/(24*E*I))
+
+
+//case(ii):maximum deflection
+yc=(W*L**3)/(E*I)*(5/384)
+
+//Result
+printf("\n Slope at the support = %0.3f  radians",-thetaA)
+printf("\n Maximum deflection = %0.3f  mm",yc)
diff --git a/3863/CH12/EX12.6/Ex12_6.sce b/3863/CH12/EX12.6/Ex12_6.sce
new file mode 100644
index 000000000..37772c289
--- /dev/null
+++ b/3863/CH12/EX12.6/Ex12_6.sce
@@ -0,0 +1,23 @@
+clear
+//Given
+//Variable declaration
+L=5*1000      //Length in mm
+L_star=5      //Length in m
+w=9           //Uniformly distributed load in kN/m
+f=7           //Bending stress in N/sq.mm
+E=1e4         //Young's modulus in N/sq.mm
+yc=10         //Central deflection in mm
+
+//Calculation
+W=w*L_star*1e3                       //Total load in N
+bd3=((W*(L**3)*12*5)/(E*yc*384))     //width X depth^3 in mm^4
+M=(W*L/8)                            //Maximum bending moment in Nmm
+bd2=(M*12/(f*2))              //width X depth^2 in mm^3
+
+d=(bd3/bd2)                   //Depth of beam in mm
+
+b=(M*12/(f*2)/d**2)          //Width of beam in mm
+
+//Result
+printf("\n Depth of beam = %0.3f  mm",d)
+printf("\n Width of beam = %0.3f  mm",b)
diff --git a/3863/CH12/EX12.7/Ex12_7.sce b/3863/CH12/EX12.7/Ex12_7.sce
new file mode 100644
index 000000000..7eb5a2292
--- /dev/null
+++ b/3863/CH12/EX12.7/Ex12_7.sce
@@ -0,0 +1,14 @@
+clear
+//Given
+//Variable declaration
+L=5*1000      //Length in mm
+f=8           //Bending stress in N/sq.mm
+yc=10         //Central deflection in mm
+E=1.2e4       //Youngs modulus in N/sq.mm
+
+//Calculation
+d=((5*L**2*(f*2*8))/(E*384*yc)*1e-1)  //Depth of beam in cm
+
+
+//Result
+printf("\n Depth of beam = %0.3f  cm",d)
diff --git a/3863/CH12/EX12.8/Ex12_8.sce b/3863/CH12/EX12.8/Ex12_8.sce
new file mode 100644
index 000000000..a5784a1e9
--- /dev/null
+++ b/3863/CH12/EX12.8/Ex12_8.sce
@@ -0,0 +1,16 @@
+clear
+//Given
+//Variable declaration
+L=6*1000       //Length in mm
+W=40*1000      //Point load in N
+a=4*1000       //Distance of point load from left support in mm
+I=7.33e7       //Moment of Inertia in mm^4
+E=2e5          //Youngs modulus in sq.mm
+
+//Calculation
+b=L-a                                 //Width of beam in mm
+yc=(-(W*a**2*b**2)/(3*E*I*L))  //Deflection under the load in mm
+
+
+//Result
+printf("\n Deflection under the load = %0.3f  mm",yc)