blob: f097d8ab808ddf9bc26e127af324fb9a2ff8ffad (
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
//
//To determine reactions, consider equilibrium equations
//variable declaration
//all Vertical loading are in KN
PL1=200.0
PL2=200.0
PL3=150.0
PL4=100.0
PL5=100.0
//length in m
UL1=6.0
UL2=8.0
UL3=9.0
UL4=UL2
UL5=UL1
L1=6.0
L2=6.0
L3=6.0
L4=6.0
L5=6.0
L6=6.0
//moment at point LO
R2=(PL1*L1+PL2*(L1+L2)+PL3*(L1+L2+L3)+PL4*(L1+L2+L3+L4)+PL5*(L1+L2+L3+L4+L5))/(L1+L2+L3+L4+L5+L6)
R1=PL1+PL2+PL3+PL4+PL5-R2
//Take the section (1)–(1) and consider the right hand side part.
U3U4=sqrt((1**2)+(UL1**2))
theta1=asin(1/U3U4)
L3U4=sqrt((UL1**2)+(UL2**2))
theta2=asin(6/L3U4)
//moment at U4
FL3L4=(R2*(L5+L6)-PL4*L4)/UL4
printf("\n FL3L4= %0.1f KN (Tension)",FL3L4)
//moment at L3
FU4U3=(-PL4*L4-PL5*(L4+L5)+R2*(L4+L5+L6))/(cos(theta1)*UL3)
printf("\n FU4U3= %0.1f KN (Comp.)",FU4U3)
//sum of horizontal forces
FL4L3=FL3L4
FU4L3=(-FL4L3+FU4U3*cos(theta1))/sin(theta2)
printf("\n FU4L3= %0.1f KN (Tension)",FU4L3)
|
