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---------------------------------------------------------------------------------------
|