initial commit / add all books

20 files changed, 708 insertions, 0 deletions

diff --git a/1898/CH6/EX6.11/Ex6_11.sce b/1898/CH6/EX6.11/Ex6_11.sce
new file mode 100755
index 000000000..195fd6560
--- /dev/null
+++ b/1898/CH6/EX6.11/Ex6_11.sce
@@ -0,0 +1,25 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.11 : ")

+

+//Given:

+l = 4.5; //m

+R1 = 1.5; //kN

+R2 = 3; //kN

+uvl = 2; //kN/m

+

+//Shear diagram:

+x = sqrt((2*R1*l)/(uvl));

+M = (R1*x) - (0.5*uvl*x^3)/(3*l);

+

+//Display:

+

+    

+ printf('\n\nV becomes zero at x = %1.1fm',x);

+ printf('\nThe magnitude of the maximum moment = %1.1f kNm',M);

+     

+//-----------------------------------------------------------------END--------------------------------------------------------------------------

+

+

+

+

diff --git a/1898/CH6/EX6.13/Ex6_13.sce b/1898/CH6/EX6.13/Ex6_13.sce
new file mode 100755
index 000000000..421d8ad0d
--- /dev/null
+++ b/1898/CH6/EX6.13/Ex6_13.sce
@@ -0,0 +1,27 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.13 : ")

+

+//Given:

+l_ab = 4; //m

+l_cd = 4; //m

+l_bc = 6; //m

+Rb = 8; //kN

+uvl = 2; //kN/m

+

+//Moment diagram:

+p = [-1/18 0 -3.6 17.6]

+x = roots(p)

+y = x(3);

+

+//Display:

+    

+ printf('\n\nV becomes zero at x = %1.2f m',y);

+

+     

+//-----------------------------------------------------------------END--------------------------------------------------------------------------

+

+

+

+

+

diff --git a/1898/CH6/EX6.14/Ex6_14.sce b/1898/CH6/EX6.14/Ex6_14.sce
new file mode 100755
index 000000000..40a46c0a7
--- /dev/null
+++ b/1898/CH6/EX6.14/Ex6_14.sce
@@ -0,0 +1,35 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.14 : ")

+

+//Given:

+b = 60; //mm

+h = 120; //mm

+sigma_max = 20; //N/mm^2

+c = b;

+

+//Part (a):

+I = (1/12)*b*h^3;

+M1 = (sigma_max*I)/(c); //sigma_max = Mc/I Flexure Formula

+M1 = M1*10^-6; //in kN/m

+

+//Part (b):

+y0=60;

+y1=-60

+

+M2 = integrate('-(20*y^2)','y',y0,y1);

+M2 = M2*10^-6;

+

+F = (0.5*sigma_max*b*b);

+c = 2*(60 -(0.5*b)); //distance between centroids of both the volumes.

+M = F*c/1000;

+

+//Display:

+    

+ printf("\n\nThe internal moment M calculated using : ");

+ printf('\na)The flexure formula = %1.2f kNm',M1);

+ printf('\nb)The resultant of the stress distribution using the basic principles = %1.2f kNm',M2);

+

+     

+//-----------------------------------------------------------------END--------------------------------------------------------------------------

+

diff --git a/1898/CH6/EX6.15/Ex6_15.sce b/1898/CH6/EX6.15/Ex6_15.sce
new file mode 100755
index 000000000..19cc03020
--- /dev/null
+++ b/1898/CH6/EX6.15/Ex6_15.sce
@@ -0,0 +1,32 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.15 : ")

+

+//Given:

+udl = 5; //kN/m

+l1 = 3;//m

+l2 = 6; //m

+t = 20/1000; //mm

+yb = 0.15;//m

+

+//Section Property:

+I_bar1 = (1/12)*(0.25)*(0.02^3);

+Ad2 = (0.25)*(0.02)*(yb+(t/2))^2;

+I_bar2 = (1/12)*(0.02)*(0.3^3);

+I = 2*(I_bar1 + Ad2) + I_bar2;

+

+//Bending stress:

+c = 0.15 + t;

+M= 22.5; //kNm

+

+sigma_max = (M*c)/(I*1000);

+

+sigma_B = (M*yb)/(I*1000);

+

+//Display:

+

+ printf('\n\nThe absolute maximum bending stress is = %1.1f MPa',sigma_max);

+

+     

+//-----------------------------------------------------------------END--------------------------------------------------------------------------

+

diff --git a/1898/CH6/EX6.16/Ex6_16.sce b/1898/CH6/EX6.16/Ex6_16.sce
new file mode 100755
index 000000000..22a9b104b
--- /dev/null
+++ b/1898/CH6/EX6.16/Ex6_16.sce
@@ -0,0 +1,38 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.16 : ")

+

+//Given:

+t1 = 15/1000; //m

+t2 = 20/1000; //m

+l = 250/1000; //m

+b = 200/1000; //m

+P = 2.4; //kN

+l_a = 2; //m

+l_b = 1; //m

+

+//Internal Moment:

+y1 = b/2;

+y2 = t2/2;

+A = (2*t1*b)+(t2*l);

+y_bar = ((2*y1*t1*b)+(y2*t2*l))/A;

+

+M = (P*l_a)+(1*y_bar);

+

+//Section Property:

+I1 = (1/12)*(l*t2^3) + (l*t2*(y_bar - y2)^2);

+I2 = (1/12)*(t1*b^3) + (t1*b*(y1 - y_bar)^2);

+I =I1+ 2*I2;

+

+//Maximum Bending Stress:

+c = b - y_bar;

+sigma_max = (M*c)/(I*1000);

+

+//Display:

+    

+  printf('\n\nThe maximum bending stress at section a-a = %1.1f MPa',sigma_max);

+

+     

+//-----------------------------------------------------------------END--------------------------------------------------------------------------

+

+

diff --git a/1898/CH6/EX6.17/Ex6_17.sce b/1898/CH6/EX6.17/Ex6_17.sce
new file mode 100755
index 000000000..0a3f795a6
--- /dev/null
+++ b/1898/CH6/EX6.17/Ex6_17.sce
@@ -0,0 +1,42 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.17 : ")

+

+//Given:

+b = 60/1000; //m

+h = 30/1000; //m

+M = 40; //Nm

+c1= h/2;

+rib_t = 5/1000; //m

+rib_w = 10/1000;//m

+

+//Without Ribs:

+I1 = (1/12)*(b*h^3);

+sigma_max1 = (M*c1)/(I1*10^6);

+

+//With Ribs:

+y1 = c1;

+y2 = h+(rib_t/2);

+A1 = h*b;

+A2 = rib_t*rib_w;

+y_bar = ((y1*A1)+2*(y2*A2))/(A1 + 2*A2);

+

+c2 = h+rib_t - y_bar;

+I2 = I1 + (b*h*(y_bar - y1)^2);

+I3 = (1/12)*rib_w*rib_t^3 + (rib_w*rib_t*(y2 - y_bar)^2);

+I = I2 + 2*I3;

+

+sigma_max2 = (M*c2)/(I*10^6);

+

+if(sigma_max2>sigma_max1)

+    

+  printf("\n\nThe maximum normal stress in the member without ribs = %1.2f MPa',sigma_max1);

+ printf("\nThe maximum normal stress in the member with ribs = %1.2f MPa',sigma_max2);

+ printf("\nThe ribs should be omitted.");

+ 

+ end

+

+     

+//-----------------------------------------------------------------END--------------------------------------------------------------------------

+

+

diff --git a/1898/CH6/EX6.18/Ex6_18.sce b/1898/CH6/EX6.18/Ex6_18.sce
new file mode 100755
index 000000000..c5938e163
--- /dev/null
+++ b/1898/CH6/EX6.18/Ex6_18.sce
@@ -0,0 +1,46 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.18 : ")

+

+//Given:

+M = 12; //kNm

+l_bc = 0.2; //m

+l_be = 0.4; //m

+

+//Internal Moment Components:

+My = (-4/5)*M;

+Mz = (3/5)*M;

+

+Iy = (1/12)*(l_be*l_bc^3);

+Iz = (1/12)*(l_bc*l_be^3); 

+

+//Bending Stress:

+sigma_B = (-Mz*1000*(l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;

+sigma_B = sigma_B/10^6;

+sigma_C = (-Mz*1000*(l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;

+sigma_C = sigma_C/10^6;

+sigma_D = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;

+sigma_D = sigma_D/10^6;

+sigma_E = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;

+sigma_E = sigma_E/10^6;

+

+//Orientation of Nuetral Axis:

+z = (0.45)/(sigma_E + sigma_B);

+

+//theta = -atan(4/3);

+tanA = (Iz/Iy)*(-4/3);

+alpha = atan(tanA);

+alpha = alpha*(180/%pi);

+

+

+//Display:

+

+    

+ printf("\n\nThe normal stress at B = %1.2f MPa',sigma_B);

+ printf("\nThe normal stress at C = %1.2f MPa',sigma_C);

+ printf("\nThe normal stress at D = %1.2f MPa',sigma_D);

+ printf("\nThe normal stress at E = %1.2f MPa',sigma_E);

+ printf("\nThe orientation of the nuetral axis = %1.1f degrees',alpha);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

diff --git a/1898/CH6/EX6.19/Ex6_19.sce b/1898/CH6/EX6.19/Ex6_19.sce
new file mode 100755
index 000000000..c5d5e7907
--- /dev/null
+++ b/1898/CH6/EX6.19/Ex6_19.sce
@@ -0,0 +1,48 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.19 : ")

+

+//Given:

+theta = 30*(%pi/180);

+M = 15; //kNm

+My = M*cos(theta); 

+Mz = M*sin(theta); 

+b = 0.1; //m

+t1 = 0.04;//m

+t2 = 0.03;//m

+

+

+//Section Properties:

+y1 = b/2;

+y2 = b + t2/2;

+A1 = (b*t1);

+A2 = (b*2*t2);

+z_bar = (y1*A1 + y2*A2)/(A1+A2);

+

+Iz = (1/12)*(b*t1^3) + (1/12)*(t2*(2*b)^3);

+Iy = (1/12)*(t1*b^3) + b*t1*(z_bar - y1)^2 + (1/12)*(2*b*t2^3) + 2*b*t2*(y2 - z_bar)^2;

+

+//Maximum Bending Stress:

+l_b = b+t2 - z_bar;

+sigma_B = (-Mz*1000*(-b))/Iz + (My*1000*(l_b))/Iy;

+sigma_B = sigma_B/10^6;

+sigma_C = (-Mz*1000*(t1/2))/Iz + (My*1000*(-z_bar))/Iy;

+sigma_C = sigma_C/10^6;

+

+sigma = max(abs(sigma_B),abs(sigma_C));

+

+//Orientation of the nuetral axis:

+theta1 = 60*(%pi/180);

+alpha = atan((Iz/Iy)*tan(theta1));

+alpha = alpha*(180/%pi);

+

+//Display:

+

+    

+ printf("\n\nThe maximum normal stress in the beam = %1.2f MPa',sigma);

+ printf("\n The orientation of the nuetral axis = %1.1f degrees',alpha);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

+

+

diff --git a/1898/CH6/EX6.20/Ex6_20.sce b/1898/CH6/EX6.20/Ex6_20.sce
new file mode 100755
index 000000000..56f0f586e
--- /dev/null
+++ b/1898/CH6/EX6.20/Ex6_20.sce
@@ -0,0 +1,39 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.20 : ")

+

+//Given:

+M =20; //kN

+Iy = 0.96*10^-3; //m^4

+Iz = 7.54*10^-3; //m^4

+theta = 57.1*(%pi/180);

+

+

+//Internal moment Components:

+My = M*sin(theta); 

+Mz = M*cos(theta); 

+

+//Bending Stress:

+y_p = -0.2; //y Coordinate of P

+z_p = 0.35; //z Coordinate of P

+

+theta1 = (%pi/2)-(theta);

+yp = -z_p*sin(theta1)+ y_p*cos(theta1);

+zp = z_p*cos(theta1) + y_p*sin(theta1);

+

+//Eq 6-17

+

+sigma_p = ((Mz*-yp)/Iz) + ((My*zp)/Iy) ;

+sigma_p = sigma_p/10^3;

+

+//Orientation of the Nuetral Axis:

+alpha = atan((Iz/Iy)*tan(theta));

+alpha = alpha*(180/%pi);

+

+//Display:

+

+    

+ printf("\n\nThe maximum normal stress at point P = %1.2f MPa',sigma_p);

+ printf("\nThe orientation of the nuetral axis = %1.1f degrees',alpha);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

diff --git a/1898/CH6/EX6.21/Ex6_21.sce b/1898/CH6/EX6.21/Ex6_21.sce
new file mode 100755
index 000000000..3aec6bca9
--- /dev/null
+++ b/1898/CH6/EX6.21/Ex6_21.sce
@@ -0,0 +1,43 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.21 : ")

+

+//Given:

+M = 2; //kNm

+Ew = 12; //GPa

+Est = 200; //GPa

+bw = 150/1000; //m

+t = 20/1000; //m

+rib = 9/1000;//m

+

+//Section Properties:

+n = (Ew/Est);

+bst = n*bw;

+

+y1 = t/2;

+A1 = t*bw;

+y2 = bw/2 + t;

+A2 = rib*bw;

+

+y_bar = (y1*A1 +y2*A2)/(A1+A2);

+

+I1 = (1/12)*(bw)*(t^3) + A1*(y_bar - y1)^2;

+I2 = (1/12)*(rib)*(bw^3) + A2*(y2-y_bar)^2;

+Ina = I1+I2;

+

+//Normal Stress:

+sigma_B = (M*(bw+t-y_bar))/(Ina*1000);

+sigma_C = (M*(y_bar))/(Ina*1000);

+

+//Normal Stress in the wood:

+sigmaB = n*sigma_B;

+

+//Display:

+

+    

+ printf("\n\nThe normal stress at point B = %1.1f MPa',sigma_B);

+ printf("\nThe normal stress at point C = %1.2f MPa',sigma_C);

+ printf("\nThe normal stress at point B in the wood = %1.2f MPa',sigmaB);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

diff --git a/1898/CH6/EX6.22/Ex6_22.sce b/1898/CH6/EX6.22/Ex6_22.sce
new file mode 100755
index 000000000..51deaa8c9
--- /dev/null
+++ b/1898/CH6/EX6.22/Ex6_22.sce
@@ -0,0 +1,50 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.22 : ")

+

+//Given:

+sigma_allow_st = 168; //MPa

+sigma_allow_w = 21; //MPa

+Est = 200; //GPa

+Ew = 12; //GPa

+Iz = 7.93*10^6; //mm^4

+A1 = 5493.75; //mm^2

+t = 5; //mm

+h = 100; //mm

+

+//Without Board:

+c = h+t;

+M1 = (sigma_allow_st*Iz)/(c*10^6);

+

+//With Board:

+bw = 300;//mm

+n = (Ew/Est); 

+bst = n*bw;

+

+//For the transformed section:

+y1 = 0;

+y2 = 55;

+A2 = bst*h;

+

+y_bar = (y1*A1 + y2*A2)/(A1+A2);

+

+I1 = Iz + A1*y_bar^2;

+I2 = (1/12)*(bst*h^3) + (A2*(y2-y_bar)^2);

+I = I1+I2;

+

+c = c+y_bar;

+M2 = (sigma_allow_st*I)/(c*10^6);

+

+cw = c - y_bar;

+Mw = (sigma_allow_w*I)/(n*cw*10^6);

+

+M = min(Mw,M2);

+

+//Display:

+

+  printf("\n\nThe maximum bending moment without re-inforcement = %1.3f kNm',M1);

+ printf("\nThe maximum bending moment with re-inforcement     = %1.2f kNm',M);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

+

diff --git a/1898/CH6/EX6.23/Ex6_23.sce b/1898/CH6/EX6.23/Ex6_23.sce
new file mode 100755
index 000000000..0045d814c
--- /dev/null
+++ b/1898/CH6/EX6.23/Ex6_23.sce
@@ -0,0 +1,37 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.23 : ")

+

+//Given:

+M = 60; //kNm

+Est = 200; //GPa

+Econc = 25; //GPa

+d = 25;//mm

+r = d/2;

+w = 300;//mm

+ht =400; //mm

+

+//Section Properties:

+n = Est/Econc;

+Ast = 2*%pi*r^2;

+A = n*Ast;

+

+p = [1 52.37 -20949.33]

+h = roots(p)

+h = h(2);

+

+I = (1/12)*(w*h^3) +w*h*(h/2)^2 + A*(ht - h)^2;

+

+//Normal Stress:

+sigma_conc_max = (M*1000*h*1000)/(I);

+sigma_conc = (M*1000*(ht-h)*1000)/(I);

+sigma_st = n*sigma_conc;

+

+//Display:

+

+    

+ printf("\n\nThe normal stress in each steel reinforcing rod = %1.2f MPa',sigma_st);

+ printf("\nThe maximum normal stress in the concrete = %1.2f MPa',sigma_conc_max);

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

+

diff --git a/1898/CH6/EX6.24/Ex6_24.sce b/1898/CH6/EX6.24/Ex6_24.sce
new file mode 100755
index 000000000..52ba2562c
--- /dev/null
+++ b/1898/CH6/EX6.24/Ex6_24.sce
@@ -0,0 +1,41 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.24 : ")

+

+//Given:

+sigma = 140; //Mpa

+ri = 90; //mm

+ro = 110; //mm

+a = 20; //mm

+

+//Section Properties:

+ 

+y = integrate('a*(1/r)','r',ri,ro)

+R = (a*a)/y;

+

+r_avg = (ri+ro)/2;

+M1 = (-sigma*a*a*ro*(r_avg - R))/(R-ro);

+M1 = M1*10^-6;

+

+M2 = (sigma*a*a*ri*(r_avg - R))/(R-ri);

+M2 = M2*10^-6;

+

+M = min(M1,M2);

+

+sigma1 = (M*(R - ro))/(a*a*ro*(r_avg - R));

+

+//For a straight Bar:

+I = (1/12)*(a*a^3);

+c = 10; //mm

+M_strt= (sigma*I)/c;

+M_strt = M_strt*10^-6;

+

+//Display:

+    

+ printf("\n\nThe maximum bending moment that can be applied to the bar = %1.3f kNm',M);

+ printf("\nThe maximum bending moment that can be applied to a straight bar = %1.3f kNm',M_strt);

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

+

+

+

diff --git a/1898/CH6/EX6.25/Ex6_25.sce b/1898/CH6/EX6.25/Ex6_25.sce
new file mode 100755
index 000000000..89531fecd
--- /dev/null
+++ b/1898/CH6/EX6.25/Ex6_25.sce
@@ -0,0 +1,39 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.25 : ")

+

+//Given:

+ri = 200/1000; //m

+r1 = 250/1000; //m

+ro = 280/1000; //m

+M = 4; //kNm

+a = 0.05; //m

+h = 0.03; //m

+

+//Section Properties:

+A1 = a^2 ;

+A2 =  (0.5*a*h);

+A = A1+A2;

+r_avg1 = (r1+ri)/2;

+r_avg2 = r1+(h/3);

+r_bar =((r_avg1*A1)+(r_avg2*A2))/A;

+

+int_dA_r1 = a*log(r1/ri);

+int_dA_r2 = (a*ro*log(ro/r1))/(ro-r1) - a;

+R = (A)/(int_dA_r1+ int_dA_r2);

+k= r_bar - R;

+

+//Normal Stress:

+sigma_B = (-M*(R-ri))/(A*ri*k*1000);

+sigma_A = (-M*(R-ro))/(A*ro*k*1000);

+

+sigma = max(abs(sigma_B),abs(sigma_A))

+

+

+//Display:

+    

+ printf("\n\nThe maximum normal stress in the bar = %1.0f MPa',sigma);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

+

diff --git a/1898/CH6/EX6.26/Ex6_26.sce b/1898/CH6/EX6.26/Ex6_26.sce
new file mode 100755
index 000000000..b7b754a12
--- /dev/null
+++ b/1898/CH6/EX6.26/Ex6_26.sce
@@ -0,0 +1,27 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.26 : ")

+

+//Given:

+M = 5; //kNm

+sigma_y = 500; //MPa

+r = 16; //mm

+h = 80; //mm

+w = 120; //mm

+r_h = r/h;

+w_h = w/h;

+k = 1.45; 

+c = h/(2000);

+t = 20/1000; //m

+

+//Calculations:

+I = (1/12)*(t)*(h/1000)^3

+sigma_max = (k*M*c)/(I*1000);

+

+//Display:

+    

+ printf("\n\nThe maximum normal stress in the steel = %1.0f MPa',sigma_max);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

+

diff --git a/1898/CH6/EX6.27/Ex6_27.sce b/1898/CH6/EX6.27/Ex6_27.sce
new file mode 100755
index 000000000..b73c9c1b0
--- /dev/null
+++ b/1898/CH6/EX6.27/Ex6_27.sce
@@ -0,0 +1,33 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.27 : ")

+

+//Given:

+sigma_y = 250; //MPa

+t = 12.5; //mm

+w = 200; //mm

+h = 225; //mm

+

+//Maximum Elastic Moment:

+yy = (h+t)/2;

+I1 = (1/12)*(w*t^3) + (w*t*yy^2);

+I = (1/12)*(t*h^3) + 2*(I1);

+c = 125; //mm

+

+My = (sigma_y*I)/(c); //Flexure Formula

+

+//Plastic Moment:

+C1= sigma_y*t*(h/2);

+C2= sigma_y*t*(w);

+Mp = (2*56.25*C1) + (2*yy*C2);

+

+//Shape Factor:

+k = Mp/My;

+

+//Display:

+

+    

+ printf("\n\nThe shape factor for the beam = %1.2f ',k);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

diff --git a/1898/CH6/EX6.28/Ex6_28.sce b/1898/CH6/EX6.28/Ex6_28.sce
new file mode 100755
index 000000000..46c556012
--- /dev/null
+++ b/1898/CH6/EX6.28/Ex6_28.sce
@@ -0,0 +1,27 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.28 : ")

+

+//Given:

+sigma_y = 250; //MPa

+t = 15/1000; //m

+w = 100/1000; //m

+h = 120/1000; //m

+c = 10/1000; //m

+

+//Calculations:

+d = ((sigma_y*t*w)+(sigma_y*t*h))/(sigma_y*t*2);

+

+T = sigma_y*t*d*10^3;

+C1 = sigma_y*t*c*10^3;

+C2 = sigma_y*t*w*10^3;

+

+Mp = (T*d/2)+(C1*c/2)+(C2*(c+t/2));

+

+//Display:

+

+    

+ printf("\n\nThe plastic moment that can be resisted by the beam = %1.1f kNm',Mp);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

diff --git a/1898/CH6/EX6.29/Ex6_29.sce b/1898/CH6/EX6.29/Ex6_29.sce
new file mode 100755
index 000000000..a40ed93f7
--- /dev/null
+++ b/1898/CH6/EX6.29/Ex6_29.sce
@@ -0,0 +1,34 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.29 : ")

+

+//Given:

+ep1 = 0.01;

+ep2 = 0.05;

+sig1 = 1050;//N/mm^2

+sig2 = 1330;//N/mm^2

+sig3 = 280; //N/mm^2

+y = 0.3; //cm

+h = 3; //cm

+w = 2; //cm

+

+//Calculations:

+yy = (h/2)-y

+T1 = (1/2)*(sig3*yy*w);

+y1 = y +(2/3)*(yy);

+T2 = yy*sig1*w;

+y2 = y+(0.5*yy);

+T3 = (0.5*y*sig1*w);

+y3 = (2/3)*(y);

+

+M = 2*(T1*y1 + T2*y2 + T3*y3);

+M = M/1000;

+

+//Display:

+

+    

+ printf("\n\nThe bending moment applied that will cause a strain of 0.05mm/mm = %1.2f kNm',M);

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

+

+

diff --git a/1898/CH6/EX6.30/Ex6_30.sce b/1898/CH6/EX6.30/Ex6_30.sce
new file mode 100755
index 000000000..6e72712b4
--- /dev/null
+++ b/1898/CH6/EX6.30/Ex6_30.sce
@@ -0,0 +1,23 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.30 : ")

+

+//Given:

+sigma_y = 250; //MPa

+t = 12.5; //mm

+w = 200; //mm

+h = 225; //mm

+c = (h/2)+t;

+I = 82.44*10^6;//mm^4

+Mp = 188; //kN

+

+//Calculations:

+sigma_allow = (Mp*10^6*c)/(I);

+y = (sigma_y*c)/(sigma_allow);

+

+//Display:

+    

+ printf("\n\nThe point of zero normal stress = %1.2f mm',y);

+ printf("\nThe Residual Stress distribution is shown in the text book.");

+ 

+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------

diff --git a/1898/CH6/EX6.5/Ex6_5.sce b/1898/CH6/EX6.5/Ex6_5.sce
new file mode 100755
index 000000000..183e0d83d
--- /dev/null
+++ b/1898/CH6/EX6.5/Ex6_5.sce
@@ -0,0 +1,22 @@
+clear all; clc;

+

+disp("Scilab Code Ex 6.5 : ")

+

+//Shear and Moment Diagrams:

+p = [-1/9 -2 30]

+x = roots(p)

+y = (x(2));

+

+   

+   M = (30*y) - (y^2) - (y^3)/27;

+

+

+

+//Display:

+    

+printf("\n\nThe magnitude of the maximum moment is = %1.0f kNm', M);

+printf('\nRefer to the shear and moment diagrams in the book.');

+

+

+//---------------------------------------------------------------------------END-----------------------------------------------------------------------------

+