diff options
Diffstat (limited to '1898/CH7')
-rwxr-xr-x | 1898/CH7/EX7.1/Ex7_1.sce | 40 | ||||
-rwxr-xr-x | 1898/CH7/EX7.2/Ex7_2.sce | 50 | ||||
-rwxr-xr-x | 1898/CH7/EX7.3/Ex7_3.sce | 39 | ||||
-rwxr-xr-x | 1898/CH7/EX7.4/Ex7_4.sce | 48 | ||||
-rwxr-xr-x | 1898/CH7/EX7.5/Ex7_5.sce | 32 | ||||
-rwxr-xr-x | 1898/CH7/EX7.6/Ex7_6.sce | 33 | ||||
-rwxr-xr-x | 1898/CH7/EX7.7/Ex7_7.sce | 33 |
7 files changed, 275 insertions, 0 deletions
diff --git a/1898/CH7/EX7.1/Ex7_1.sce b/1898/CH7/EX7.1/Ex7_1.sce new file mode 100755 index 000000000..921a97e4d --- /dev/null +++ b/1898/CH7/EX7.1/Ex7_1.sce @@ -0,0 +1,40 @@ +clear all; clc;
+
+disp("Scilab Code Ex 7.1 : ")
+
+//Given:
+V = 3; //kN
+h = 125; //mm
+b = 100; //mm
+y_top = 50; //mm
+x_right = 37.5; //mm
+
+//Part (a):
+
+//Section Properties:
+I = (b*h^3)/12;
+y_dash_1 = ((h-y_top)-(h/2));
+A = y_top*b;
+Q = (y_dash_1+(y_top/2))*A;
+
+//Shear Stress:
+tou_p = (V*Q)/(I*b); //tou = VQ/It
+tou_p = tou_p*10^3;
+
+//Part (b):
+
+//Section Properties:
+y_dash_2 = (y_dash_1+(y_top));
+a_dash= b*y_dash_2;
+Q_dash =(y_dash_2*a_dash)/2;
+
+//Shear Stress:
+tou_max = (V*Q_dash)/(I*b);
+tou_max = tou_max*10^3;
+
+//Display:
+
+printf("\n\nThe shear stress in the beam at point P = %1.3f MPa',tou_p);
+printf('\nThe maximum shear stress in the beam = %1.3f MPa',tou_max);
+
+//----------------------------------------------------------------------END--------------------------------------------------------------------------------
diff --git a/1898/CH7/EX7.2/Ex7_2.sce b/1898/CH7/EX7.2/Ex7_2.sce new file mode 100755 index 000000000..8b624ffb1 --- /dev/null +++ b/1898/CH7/EX7.2/Ex7_2.sce @@ -0,0 +1,50 @@ +clear all; clc;
+
+disp("Scilab Code Ex 7.2 : ")
+
+//Given:
+V = 80; //kN
+thick_1 = 20/1000; //m
+thick_2 = 15/1000; //m
+l = 300/1000; //m
+y = 100/1000; //m
+h = 2*y;
+y_dash = y +thick_1/2;
+
+//Part(a):
+
+I1 = (thick_2*(h^3))/12;
+I2 = (l*(thick_1^3))/12;
+I3 = (l*thick_1*(y_dash)^2);
+I = I1+2*(I2+I3); //Moment of inertia
+
+Q_b = y_dash*l*thick_1;
+//At B'
+tou_b_dash = (V*Q_b)/(I*l*1000);
+//At B
+tou_b = (V*Q_b)/(I*thick_2*1000);
+
+//At C:
+Q_c = (y_dash*l*thick_1)+(y*thick_2*y/2);
+tou_c = (V*Q_c)/(I*thick_2*1000);
+
+//Part(b)
+
+
+y0 = -0.1;
+y1 = 0.1;
+
+function Q =f(y),Q = ((0.735 - (7.5*y*y))*10^-3),
+endfunction
+Int =intg(y0,y1,f)
+
+V_w = (V*Int*thick_2)/(I*thick_2);
+
+//Display:
+
+printf("\n\nThe shear stress at B dash = %1.2f MPa',tou_b_dash);
+printf("\nThe shear stress at B = %1.1f MPa',tou_b);
+printf("\nThe shear stress at C = %1.1f MPa',tou_c);
+printf("\nThe shear force resisted by the web = %1.1f kN',V_w);
+
+//------------------------------------------------------------------------END----------------------------------------------------------------------------------
diff --git a/1898/CH7/EX7.3/Ex7_3.sce b/1898/CH7/EX7.3/Ex7_3.sce new file mode 100755 index 000000000..f236accae --- /dev/null +++ b/1898/CH7/EX7.3/Ex7_3.sce @@ -0,0 +1,39 @@ +clear all; clc;
+
+disp("Scilab Code Ex 7.3 : ")
+
+//Given:
+udl = 6.5; //kN
+l_bc = 8; //m
+l = 150/1000;//m
+t = 30/1000;//m
+
+//Internal Shear:
+w = udl*l_bc/2;
+l_wc = l_bc/4;
+l_bw = l_bc - l_wc;
+V = (w*l_bw)/l_bc;
+R_b = w - V;
+
+//Section Properties:
+y1= l/2;
+A = (l*t);
+y2= l+(t/2);
+y_dash = (y1*A + y2*A)/(2*A);
+I1 = (t*l^3)/12;
+I2 = (A*(y_dash-y1)^2);
+I3 = (l*t^3)/12;
+I4 = (A*(y2 - y_dash)^2);
+I = I1+I2+I3+I4;
+
+Q = ((l+t)-(t/2)-y_dash)*A;
+
+//Shear Stress:
+tou_max = (V*Q)/(I*t*1000);
+
+//Display:
+
+printf("\n\nThe maximum shear stress in the glue necessary to hold the boards together = %1.2f MPa',tou_max);
+
+
+//----------------------------------------------------------------------END--------------------------------------------------------------------------------
diff --git a/1898/CH7/EX7.4/Ex7_4.sce b/1898/CH7/EX7.4/Ex7_4.sce new file mode 100755 index 000000000..73c377ca9 --- /dev/null +++ b/1898/CH7/EX7.4/Ex7_4.sce @@ -0,0 +1,48 @@ +clear all; clc;
+
+disp("Scilab Code Ex 7.4 : ")
+
+//Given:
+
+V = 850; //kN
+l1 =250/1000; //m
+l2 = 300/1000; //m
+l3 = 125/1000;//m
+t = 10/1000; //m
+h = 200/1000; //m
+
+A1 = l1*t;
+A2 = l2*t;
+A3 = l3*t;
+
+y1 = l2+(t/2);
+y2 = l2/2;
+y3 = h+(t/2);
+
+y_dash = (2*y2*A2 + A1*y1 + A3*y3)/(2*A2 + A1 + A3);
+
+I1 = ((l1*t^3)/12) +(A1 * (l2+(t/2)-y_dash)^2);
+I2 = ((t*l2^3)/12) +(A2 * (y_dash - (l2/2))^2);
+I3 = ((l3*t^3)/12) +(A1 * (h+(t/2)-y_dash)^2);
+I = 2*I2 + I1 + I3;
+
+Q_b = (l2+(t/2) - y_dash)*A1; //Q = y'A'
+Q_c = (h+(t/2) - y_dash)*A3; //Q = y'A'
+
+//Shear Flow:
+
+q_b = (V*Q_b)/I;
+q_c = (V*Q_c)/I;
+
+q_b = q_b/(2*1000);
+q_c = q_c/(2*1000);
+
+//Display:
+
+printf("\n\nThe shear flow at B, resisted by the glue is = %1.2f MN/m',q_b);
+printf("\nThe shear flow at C, resisted by the glue is = %1.4f MN/m',q_c);
+
+
+
+
+//----------------------------------------------------------------------END--------------------------------------------------------------------------------
diff --git a/1898/CH7/EX7.5/Ex7_5.sce b/1898/CH7/EX7.5/Ex7_5.sce new file mode 100755 index 000000000..de724633b --- /dev/null +++ b/1898/CH7/EX7.5/Ex7_5.sce @@ -0,0 +1,32 @@ +clear all; clc;
+
+disp("Scilab Code Ex 7.5 : ")
+
+//Given:
+V = 80; //N
+t = 1.5; //cm
+a = 7.5; //cm
+b = a-2*t; //cm
+F_nail= 30; //N
+
+//Section Properties:
+I = (a*a^3 - b*b^3 )/12;
+Q_b = (((a-2*t)/2)+(t/2))*a*t; //Q = y'A'
+Q_c = (((a-2*t)/2)+(t/2))*(a-2*t)*t; //Q = y'A'
+
+//Shear Flow:
+q_b = (V*Q_b)/I;
+q_c = (V*Q_c)/I;
+
+s_b = F_nail/(q_b/2);
+s_c = F_nail/(q_c/2);
+
+//Display:
+
+
+printf("\n\nThe maximum spacing of nails required at B is = %1.0f cm',s_b);
+printf("\nThe maximum spacing of nails required at C is = %1.1f cm',s_c);
+
+
+
+//----------------------------------------------------------------------END--------------------------------------------------------------------------------
diff --git a/1898/CH7/EX7.6/Ex7_6.sce b/1898/CH7/EX7.6/Ex7_6.sce new file mode 100755 index 000000000..db891753f --- /dev/null +++ b/1898/CH7/EX7.6/Ex7_6.sce @@ -0,0 +1,33 @@ +clear all; clc;
+
+disp("Scilab Code Ex 7.6 : ")
+
+//Given:
+F = 40; //N
+s = 9; //cm
+h = 5; //cm
+t = 0.5; //cm
+w = 3; //cm
+w_3 = w/3; //cm
+
+//Calculations:
+
+I = (w*h^3)/12 - (2*w_3*(h - 2*t)^3)/12;
+
+//Case 1:
+
+Q1 = ((h-t)/2)*(w*t);
+V1 =((F/s)*I)/Q1 ; //q = VQ/I
+
+//Case2:
+
+Q2 = ((h-t)/2)*(w_3*t);
+V2 =((F/s)*I)/Q2 ; //q = VQ/I
+
+//Display:
+
+
+printf("\n\nThe largest vertical shear that can be supported in Case 1 = %1.1f N',V1);
+printf("\nThe largest vertical shear that can be supported in Case 2 = %1.1f N',V2);
+
+//-------------------------------------------------------------------------END---------------------------------------------------------------------------------------
diff --git a/1898/CH7/EX7.7/Ex7_7.sce b/1898/CH7/EX7.7/Ex7_7.sce new file mode 100755 index 000000000..6e752ac58 --- /dev/null +++ b/1898/CH7/EX7.7/Ex7_7.sce @@ -0,0 +1,33 @@ +clear all; clc;
+
+disp("Scilab Code Ex 7.7 : ")
+
+//Given:
+V = 10; //kN
+b1 = 6; //cm
+h1 = 8; //cm
+t = 1; //cm
+b2 = b1-2*t;
+h2 = h1-2*t; //cm
+b3 = 4; //cm
+
+//Calculations:
+I = ((b1*h1^3)/12) - ((b2*h2^3)/12);
+
+q_b = 0;
+
+Q_c = ((b1/2)+(t/2))*(b3+(t))*t;
+q_c = (V*Q_c*100)/(I); //Q = VQ/I
+
+Q_d = (2*h1/4*t*b3) + ((b1/2)+(t/2))*b3*t;
+q_d = (V*Q_d*100)/(I); //Q = VQ/I
+
+//Display:
+
+
+printf("\n\nVariation of shear flow at B = %1.1f N/mm',q_b);
+printf("\nVariation of shear flow at C = %1.1f N/mm',q_c);
+printf('\nVariation of shear flow at D = %1.1f N/mm',q_d);
+
+//-------------------------------------------------------------------------END---------------------------------------------------------------------------------------
+
|