blob: 2f5513af749cdbe7296376f795c3d1891e8a54da (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
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----------------------------------------------------------------------------------
|
