blob: ed5476cb87f5f2facbb7b8cbfad18aedce03a1eb (
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
|
clc
//
// variable declaration
E=(29*(10**6)) // Modulus of elastoplasticity(psi)
sY=50 // Stress(ksi)
// Calculation
//Case(a) Onset Of Yield
I=((1/12.0)*(12)*((16**3))-(1/12.0)*(12-0.75)*((14**3))**0) // Centroidal moment of inertia(in**4)
//Bending Moment
sMAX=sY // Stress(ksi)
c=8.0 // Distance(in)
My=(sY*I)/c // Bending moment(kip.in)
//Radius of Curvature
Ey=sY/(E*(1.0)) // Strain
pY=(c/Ey)/(1000.0) // Radius of curvature(in)
//Case(b) Flanges Fully Plastic
R1=50*12*1 // Compressive forces on top(kips)
R4=R1 // Compressive forces on top(kips)
R2=((1/2.0)*(50)*(7)*(0.75)+0.05) // Compressive forces on top half(kips)
R3=R2 // Compressive forces on top half(kips)
//Bending Moment
M=2*((R1*7.5)+(R2*4.67)) // Bending moment(kip.in)
//Radius of Curvature
p=(((7/0.001724)*0.0833)) // Radius of curvature(ft)
// Result
printf("\n Case(a) Bending moment = %0.3f kip.in' ,My)
printf("\n Case(a) Radius of curvature = %0.3f in' ,pY)
printf("\n Case(b) Bending moment = %0.3f kip.in' ,M)
printf("\n Case(b) Radius of curvature = %0.3f ft' ,p)
|
