diff options
Diffstat (limited to '1898/CH12/EX12.10/Ex12_10.sce')
-rwxr-xr-x | 1898/CH12/EX12.10/Ex12_10.sce | 38 |
1 files changed, 38 insertions, 0 deletions
diff --git a/1898/CH12/EX12.10/Ex12_10.sce b/1898/CH12/EX12.10/Ex12_10.sce new file mode 100755 index 000000000..7446536d2 --- /dev/null +++ b/1898/CH12/EX12.10/Ex12_10.sce @@ -0,0 +1,38 @@ +clear all; clc;
+
+disp("Scilab Code Ex 12.10 : ")
+
+//Given:
+E = 200*10^6; //kN/m^2
+I = 17*10^-6;//mm^4
+l_ac = 2; //m
+l_cF = 4; //m
+l_Fb = 2; //m
+l_cb = 6; //m
+l_aF = 6; //m
+l_ab = 8; //m
+F = 16; //kN
+R_b = (F*l_cb)/l_ab;
+R_a = F - R_b;
+
+mc = R_a*l_ac;
+mf = R_b*l_Fb;
+theta_ca = (0.5*l_ac*mc)/(E*I);
+
+A1 = 0.5*l_aF*mf;
+t1_ba = (l_Fb + l_aF/3)*(A1);
+
+A2 = 0.5*l_Fb*mf;
+t2_ba = (l_Fb*2*A2)/3;
+
+t_ba = (t1_ba+t2_ba)/(E*I);
+
+theta_c = (t_ba/l_ab)-(theta_ca);
+
+//Display:
+
+printf("\n\nThe slope at point C of the steel beam = %1.5f rad',theta_c);
+
+
+//----------------------------------------------------------------------END--------------------------------------------------------------------------------
+
|