blob: 14e5cb04eb44e87270cb3cf5edd63ced9ebe8bcf (
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
|
clear all; clc;
disp("Scilab Code Ex 4.3 : ")
//Given:
d_ac = 20; //mm
r_ac = d_ac/(2*1000); //radius in m
d_bd =40; //mm
r_bd = d_bd/(2*1000); //radius in m
P = 90; //kN
E_st = 200*(10^9); //Pa
E_al = 70*(10^9); //Pa
l_af = 200; //mm
l_fb = 400; //mm
l_bd = 300; //mm
l_ac = l_bd;
//Calculations:
//Internal Force:
P_ac = 60; //kN
P_bd = 30; //kN
//Displacement:
//Post AC: delta = PL/AE
num1 = -(P_ac*10^3*(l_ac/1000));
denom1 = %pi* r_ac^2*E_st;
delta_a = -num1/denom1; //downwards
delta_a = delta_a*1000; //in m
//Post BD: delta = PL/AE
num2 = -(P_bd*10^3*(l_bd/1000));
denom2 = %pi* r_bd^2*E_al;
delta_b = -num2/denom2; //downwards
delta_b = delta_b*1000; //in m
delta_f = delta_b + (0.184)*(l_fb/(l_af+l_fb)); //By similar triangles from the figure.
//Display:
printf('\n\nThe displacement of Post AC = +%1.3f mm downwards',delta_a);
printf('\nThe displacement of Post BD = +%1.3f mm downwards',delta_b);
printf('\nnThe displacement of point F = +%1.3f mm downwards',delta_f);
//------------------------------------------------------------------------------END-----------------------------------------------------------------------------------------
|
