2 files changed, 43 insertions, 0 deletions
diff --git a/3289/CH9/EX9.1/Ex9_1.sce b/3289/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..37a1bec18
--- /dev/null
+++ b/3289/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,22 @@
+clc
+
+w=0.1 //m
+d=0.115 //m
+l=4 //m
+p=175 //kN/m
+k=14*10^6 //Pa
+E=200*10^9 //Pa
+I=(0.1*(0.15)^3)
+
+//deltav=(p/2*k)*derivative(x)*beta*exp^(betax)*(cos beta(x)+sin beta(x))
+//vA=(p/2k)*(2-exp^(betaa)*cos betaa - exp^(betab)*cos betab)
+
+beta=(k/(4*E*I/12))^(0.25)
+disp(beta,"in meter inverse is= ")
+
+vmax=(p*(2-(-0.0345)-(0.0345)))/(2*14000)
+disp(vmax,"in meter is= ")
+z=k*vmax
+disp(z,"maxi force per unit of length between beam & foundation in kN/m is= ")
+
+// Ans varies due to round of error
diff --git a/3289/CH9/EX9.4/Ex9_4.sce b/3289/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..00ffc5f91
--- /dev/null
+++ b/3289/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,21 @@
+clc
+
+a=1.5 //m
+E=206.8*10^9 //Pa
+K=10000 //N/m
+I=6*10^-6 //m^4
+P=6700 //N
+c=0.05
+
+k=K/a
+disp(k,"foundation modulus of the equivalent continuous elastic support in Pa is=")
+
+beta=(k/(4*E*I))^(1/4)
+disp(beta)
+
+//sigmamax=(M*c/I)=(P*c/4*beta*I)
+sigmamax=((P*c)/(4*beta*I))
+disp(sigmamax,"in Pa is=")
+
+vmax=(P*beta)/(2*k)
+disp(vmax,"in meter is=")