summaryrefslogtreecommitdiff
path: root/1898/CH6/EX6.23/Ex6_23.sce
blob: 0045d814cc2952ba37e35ac76bd441143a24d2d5 (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
clear all; clc;

disp("Scilab Code Ex 6.23 : ")

//Given:
M = 60; //kNm
Est = 200; //GPa
Econc = 25; //GPa
d = 25;//mm
r = d/2;
w = 300;//mm
ht =400; //mm

//Section Properties:
n = Est/Econc;
Ast = 2*%pi*r^2;
A = n*Ast;

p = [1 52.37 -20949.33]
h = roots(p)
h = h(2);

I = (1/12)*(w*h^3) +w*h*(h/2)^2 + A*(ht - h)^2;

//Normal Stress:
sigma_conc_max = (M*1000*h*1000)/(I);
sigma_conc = (M*1000*(ht-h)*1000)/(I);
sigma_st = n*sigma_conc;

//Display:

    
 printf("\n\nThe normal stress in each steel reinforcing rod = %1.2f MPa',sigma_st);
 printf("\nThe maximum normal stress in the concrete = %1.2f MPa',sigma_conc_max);
 //------------------------------------------------------------------------END---------------------------------------------------------------------------------------