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/CH1/EX1.7/Ex1_7.sce | 45 +++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 45 insertions(+) create mode 100755 1898/CH1/EX1.7/Ex1_7.sce (limited to '1898/CH1/EX1.7') diff --git a/1898/CH1/EX1.7/Ex1_7.sce b/1898/CH1/EX1.7/Ex1_7.sce new file mode 100755 index 000000000..2f5513af7 --- /dev/null +++ b/1898/CH1/EX1.7/Ex1_7.sce @@ -0,0 +1,45 @@ +clear all; clc; + + +disp("Scilab Code Ex 1.7 :") + +//Given : +m_lamp = 80; //Mass of lamp in Kg. +d_ab = 10; // Diameter of AB in mm. +d_bc = 8; // Diameter of BC in mm. +ab_h = 60 *(%pi/180); // In degrees - Angle made by AB with the horizontal. +w = m_lamp*9.81; //N +a_bc = (%pi/4)*(d_bc^2); //m^2 Area of cross section of rod BC +a_ab = (%pi/4)*(d_ab^2); //m^2 Area of cross section of rod AB + + + +// Equations of equilibrium: Solving equilibrium equations simultaneously ,using matrices ,in the x and y directions to obtain force in BC and force in BA. + + +a = [(4/5) -(cos(ab_h)) ; (3/5) (sin(ab_h))]; +b = [0 ; w]; +f = zeros(1) + +f = a\b; +f_bc = f(1); // Force in BC in N. +f_ba = f(2); //Force in BA in N. +avg_normal_stress_a = f_ba / a_ab; //Mpa Average Normal Stress in AB +avg_normal_stress_c = f_bc/ a_bc;// Mpa Average Normal Stress in BC + + +// Displaying results: + + +printf('\n\nThe Weight of lamp = %.2f N',w); +printf('\nThe Net force in BC = %.2f N',f_bc); +printf('\nTheNet force in BA = %.2f N',f_ba); +printf('\nThe Average Normal Stress in AB when subjected to load = %.2f MPa',avg_normal_stress_a); +printf('\nThe Average Normal Stress in BC when subjected to load = %.2f MPa',avg_normal_stress_c); + +//------------------------------------------------------------------END---------------------------------------------------------------------------------- + + + + + -- cgit