1 files changed, 37 insertions, 0 deletions
diff --git a/3772/CH13/EX13.4/Ex13_4.sce b/3772/CH13/EX13.4/Ex13_4.sce
new file mode 100644
index 000000000..0cc8d1ef0
--- /dev/null
+++ b/3772/CH13/EX13.4/Ex13_4.sce
@@ -0,0 +1,37 @@
+// Problem no 13.4,Page No.303
+
+clc;clear;
+close;
+
+F=100*10**3 //N //Shear Force
+I=11340*10**-8 //m**4 //M.I
+b=20 //cm //width of Flange
+t=5 //cm //thickness of flange
+w_d=20 //cm //Depth of web
+w_t=5 //cm //thickness of web
+
+//Calculations
+
+a_1=b*t //cm**2 //Area of flange
+a_2=w_d*w_t //cm**2 //Area of web
+y_1=t*2**-1 //cm //C.G of flange
+y_2=t+w_d*2**-1
+
+//C.G of C/s
+Y=(a_1*y_1+a_2*y_2)*(a_1+a_2)**-1
+
+//Shear Stress in flange at the junction with web
+//Let tou(Shear stress)=S 
+//Change in the notifications of Shear Stress For convenience
+S_1=(F*a_1*(Y-y_1)*10**-6)*(I*b*10**-2)**-1*10**-3 
+
+//Shear Stress in web at the junction with flange 
+S_2=(F*a_1*(Y-y_1)*10**-6)*(I*w_t*10**-2)**-1*10**-3 
+
+//Max Shear Stres at N.A
+S_max=(F*(a_1*(Y-y_1)+(w_t*(Y-t))*((Y-t)*2**-1))*10**-6)*(I*w_t*10**-2)**-1*10**-3
+
+//Result
+printf("Shear Stress in flange at the junction with web %.2f",S_1);printf(" KN/m**2")
+printf("\n Shear Stress in web at the junction with flange %.2f",S_2);printf(" KN/m**2")
+printf("\n Max Shear Stress at N.A %.2f",S_max);printf(" KN/m**2")