diff options
Diffstat (limited to '3764/CH4/EX4.5/Ex4_5.sce')
| -rw-r--r-- | 3764/CH4/EX4.5/Ex4_5.sce | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/3764/CH4/EX4.5/Ex4_5.sce b/3764/CH4/EX4.5/Ex4_5.sce new file mode 100644 index 000000000..ed5476cb8 --- /dev/null +++ b/3764/CH4/EX4.5/Ex4_5.sce @@ -0,0 +1,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) |