From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 1898/CH7/EX7.2/Ex7_2.sce | 50 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 50 insertions(+) create mode 100755 1898/CH7/EX7.2/Ex7_2.sce (limited to '1898/CH7/EX7.2') 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---------------------------------------------------------------------------------- -- cgit