summaryrefslogtreecommitdiff
path: root/3862/CH10/EX10.4/Ex10_4.sce
blob: d71fd3afc47a3d7e026fc130c47321c6ab75ff87 (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
56
57
clear
//Compare the moment carrying capacity of the section given in example 10.3 with equivalent section of the same area but (i) square section (ii) rectangular section with depth twice the width and (iii) a circular section.

//
//variable declaration

A=180.0*10.0+380.0*10.0+180.0*10.0

//If ‘a’ is the size of the equivalent square section, 

a=(sqrt(A))       //mm

I=(a*(a**3))/12  //Moment of inertia of this section** mm^4

ymax=a/2

Z=I/ymax

f=150.0 

Mcc=f*Z          //Moment carrying capacity

MccI=136985000.0

Ratio=MccI/Mcc
printf("\n (i) Moment carryingcapacity of Isection/ Moment carryingcapacityof equivalent squaresection= %0.3f  ",Ratio)


//Equivalent rectangular section of depth twice the width. Let b be the width,Depth d = 2b. Equating its area to area of I-section,we get
b=sqrt(7400/2)

ymax=b

I=b*(((2*b)**3))/12
 
M=f*I/ymax


MccI=136985000

Ratio=MccI/M
printf("\n (ii) Moment carryingcapacity of I-section/ Moment carryingcapacityof equivalent squaresection= %0.3f  ",Ratio)

//Equivalent circular section. Let diameter be d.

d=sqrt(7400*4/%pi)

I=(%pi*(d**4))/64
ymax=d/2
Z=I/ymax
fper=(150)
M=fper*Z

MccI=136985000

Ratio=MccI/M
printf("\n (i) Moment carryingcapacity of Isection/ Moment carryingcapacityof equivalent squaresection= %0.3f  ",Ratio)