blob: 3121f5351d89d4f81fdbff260ba3dde761e4ad05 (
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
46
47
48
49
50
51
52
53
54
55
|
clear all; clc;
disp("Scilab Code Ex 4.15 : ")
//Given
weight = 15; //kN
l_ab = 5; //m
l_ac= 5.0075; //m
area = 30; //mm^2
//calculations:
strain_ab = (l_ac-l_ab)/l_ab;
max_strain = 0.0017;
stress_ab = (350*strain_ab)/max_strain;
F_ab = stress_ab*area; // F= stress*area
E_st = 350/max_strain; //Modulus ofelasticity
del1 = l_ab/(area*10^-6*E_st*10^3); //del = PL/AE
del2 = l_ac/(area*10^-6*E_st*10^3); //del = PL/AE
//Eqn1 = T_ab + T_ac = weight
//Eqn2 = del1*T_ab - del2*T_ac = (l_ac-l_ab)
//Solving using matrices:
A = [1 1;del1 -del2];
b = [weight; (l_ac-l_ab)];
T = A\b;
T_ab = T(1);
T_ac = T(2);
stress_in_ab = (T_ab*10^3)/area;
if(stress_in_ab>350)
T_ab = (350*area)/1000;
end
T_ac = 15-T_ab;
stress = (T_ac*10^3)/area;
strain_ac = (stress*max_strain)/350;
elong_ac = strain_ac*l_ac; //m
elong_ab = (l_ac-l_ab)+elong_ac; //m
//Display:
printf('\n\nThe force experienced by wire AB = %1.1f kN',T_ab);
printf('\nThe force experienced by wire AC = %1.1f kN',T_ac);
printf('\nThe elongation in wire AB = %1.5f m',elong_ab);
printf('\nThe elongation in wire AC = %1.5f m',elong_ac);
//---------------------------------------------------------------------------END------------------------------------------------------------------------
|
