Added new codeHEADmaster

1 files changed, 72 insertions, 0 deletions

diff --git a/3864/CH5/EX5.18/Ex5_18.sce b/3864/CH5/EX5.18/Ex5_18.sce
new file mode 100644
index 000000000..26a9324a6
--- /dev/null
+++ b/3864/CH5/EX5.18/Ex5_18.sce
@@ -0,0 +1,72 @@
+clear
+//
+//
+
+//Initilization of Variables
+L_ED=2 //m //Length of DB & AC
+L_AC=2
+L_DB=2
+L_CD=4 //m //Length of CD
+L_CE=2 //m //Length of CE
+F_A=40 //KN //Force at C
+F_B=20 //KN //Force at A
+E=200*10**6 //KN/mm**2 //Modulus of Elasticity
+I=50*10**-6 //m**4 //M.I
+
+//Calculations
+
+//LEt V_C & V_D be the reactions at C &  D respectively
+//V_C+V_D=60
+
+//Taking Moment At D,M_D
+V_C=-(-F_A*(L_AC+L_CE+L_ED)+F_B*L_DB)*L_CD**-1
+V_D=60-V_C
+
+//Now Taking Moment at Distance x from A,
+//M_x=-40*x+50*(x-2)+10*(x-6)
+
+//EI*(d**2*y/dx**2)=-40*x+50*(x-2)+10*(x-6)
+
+//Now Integrating above Equation we get
+//EI*(dy/dx)=C1+20*x**2-25*(x-2)+5*(x-6)**2
+
+//Again Integrating above Equation we get
+//EI*y=C2+C1*x-20*3**-1*x**3+25*3**-1*(x-2)**3+5*3**-1*(x-6)**3
+
+//At
+x=0
+y=0
+//C2+2*C1=-53.33    ...............(1)
+
+//At 
+x=6
+y=0
+//C2+6*C1=906.667    ...............(2)
+
+//Subtracting Equation 1 from 2 we get
+C1=853.333*4**-1
+C2=53.333-2*C1
+x=0
+y_A=(C2+C1*x-20*3**-1*x**3+25*3**-1*(x-2)**3+5*3**-1*(x-6)**3)*(E*I)**-1
+
+//Answer For y_A is incorrect in textbook
+
+//At Mid-span
+C1=853.333*4**-1
+C2=53.333-2*C1
+x=4
+y_E=(C2+C1*x-20*3**-1*x**3+25*3**-1*(x-2)**3+5*3**-1*(x-6)**3)*(E*I)**-1
+
+//Answer For y_E is incorrect in textbook
+
+//At B
+C1=853.333*4**-1
+C2=53.333-2*C1
+x=8
+y_B=(C2+C1*x-20*3**-1*x**3+25*3**-1*(x-2)**3+5*3**-1*(x-6)**3)*(E*I)**-1
+
+
+//Result
+printf("\n Deflection relative to the level of the supports:at End A %0.4f  mm",y_A)
+printf("\n                                                 :at End B %0.4f  mm",y_B)
+printf("\n                                                 :at Centre of CD %0.4f  mm",y_E)