blob: 74d9c8acf5f3f510ac7eaa586ba62d2741bf0940 (
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
|
clear
//Given
M = 50000 //ft-lb , positive bending moment applied
N = 9 // number of steel bars
n = 15 // The ratio of steel to concrete
A_s = 30 //in2 area of steel in concrete
//(10*y)*(y/2) = 30*(20-y)
//y**2 + 6*y -120
//solving quadractic equation
//
a = 1
b = 6
c = -120
// calculate the discriminant
d = (b**2) - (4*a*c)
// find two solutions
sol1 = (-b-sqrt(d))/(2*a)
sol2 = (-b+sqrt(d))/(2*a)
y = sol2 // Nuetral axis is found
l_1 = y //in- the concrete below nuetral axis is not considered
b_1 = 10 //in - width
A_1 = l_1* b_1 //in2 - area of concrete
y_1 = y/2 //in com of the concrete
y_2 = 20-y //in com of the transformed steel
I_1 = b_1*(l_1**3)/12.0 + A_1*((y_1-y)**2) //in^4 parallel axis theorm
I_2 = A_s*((y_2)**2) //in^4 first part is neglected
I_net = I_1 + I_2 //in^4 - the total moment of inertia
c_c= y //in The maximum distance in concrete
stress_concrete = M*12*c_c/I_net //psi - The maximum stress in concrete
c_s= 20- y
stress_steel =n*M*12*c_s/I_net //psi - The maximum stress in concrete
printf("\n The maximum stress in concrete %0.2f psi",stress_concrete) //
printf("\n The stress in steel %0.2f psi",stress_steel)
printf("\n answer varies due to rounding off errors")
|
