1 files changed, 52 insertions, 0 deletions
diff --git a/3864/CH4/EX4.4/Ex4_4.sce b/3864/CH4/EX4.4/Ex4_4.sce
new file mode 100644
index 000000000..40556ec38
--- /dev/null
+++ b/3864/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,52 @@
+clear
+//
+
+//Initilization of Variables
+
+//Flange (Top)
+b1=80 //mm //Width 
+t1=40 //mm //Thickness
+
+//Flange (Bottom)
+b2=160 //mm //width
+t2=40 //mm //Thickness
+
+//web
+d=120 //mm //Depth
+t3=20 //mm //Thickness
+
+D=200 //mm //Overall Depth
+sigma1=30 //N/mm**2 //Tensile stress
+sigma2=90 //N/mm**2 //Compressive stress
+L=6000 //mm //Span
+
+//Calculations
+
+//Distance of centroid from bottom fibre
+y_bar=(b1*t1*(D-t1*2**-1)+d*t3*(d*2**-1+t2)+b2*t2*t2*2**-1)*(b1*t1+d*t3+b2*t2)**-1 //mm
+
+//Moment of Inertia
+I=1*12**-1*b1*t1**3+b1*t1*(D-t1*2**-1-(y_bar))**2+1*12**-1*t3*d**3+t3*d*(d*2**-1+t2-(y_bar))**2+1*12**-1*b2*t2**3+b2*t2*(t2*2**-1-(y_bar))**2
+
+
+//Extreme fibre distance of top and bottom fibres are y_t and y_c respectively
+
+y_t=y_bar //mm
+y_c=D-y_bar //mm
+
+//Moment carrying capacity considering Tensile strength 
+M1=sigma1*I*y_t**-1*10**-6  //KN-m
+
+//Moment carrying capacity considering compressive strength 
+M2=sigma2*I*y_c**-1*10**-6  //KN-m
+
+//Max Bending moment in simply supported beam 6 m due to u.d.l
+//M_max=w*L*10**-3*8**-1
+//After simplifying further we get
+//M_max=4.5*w
+
+//Now Equating it to Moment carrying capacity, we get load carrying capacity
+w=M1*4.5**-1 //KN/m
+
+//Result
+printf("\n Max Uniformly Distributed Load is %0.3f  KN/m",w)