6 files changed, 70 insertions, 0 deletions
diff --git a/564/CH10/EX10.1/10_1.sce b/564/CH10/EX10.1/10_1.sce
new file mode 100755
index 000000000..e97916700
--- /dev/null
+++ b/564/CH10/EX10.1/10_1.sce
@@ -0,0 +1,13 @@
+pathname=get_absolute_file_path('10_1.sce')
+filename=pathname+filesep()+'10_1data.sci'
+exec(filename)
+lambda=(m*l^3)/(3*48*EI);//λ
+x=poly(0,'x');
+y=21*x^2 -22*x +1;
+m=roots(y);
+omega1=(m(1,:)/lambda)^0.5;//ω1
+omega2=(m(1,:)/lambda)^0.5;//ω2
+f1=omega1/(2*%pi);
+f2=omega2/(2*%pi);
+printf("\nf1:%f",f1);
+printf("\nf2:%f",f2);
+\ No newline at end of file
diff --git a/564/CH10/EX10.2/10_2.sce b/564/CH10/EX10.2/10_2.sce
new file mode 100755
index 000000000..13348740a
--- /dev/null
+++ b/564/CH10/EX10.2/10_2.sce
@@ -0,0 +1,12 @@
+pathname=get_absolute_file_path('10_2.sce')
+filename=pathname+filesep()+'10_2data.sci'
+exec(filename)
+lambda=(m*l^3)/(6*EI);//λ
+a=188,b=-44,c=1
+m1=(-b+(b*b -4*a*c)^0.5)/(2*a);
+m2=(-b-(b*b -4*a*c)^0.5)/(2*a);
+omega1=(m1/lambda)^0.5;//ω1
+omega2=(m2/lambda)^0.5;//ω2
+f1=[omega1;omega2];
+f=min(f1)/(2*%pi);
+printf("\nlowest natural frequency is: %f",f);
+\ No newline at end of file
diff --git a/564/CH10/EX10.4/10_4.sce b/564/CH10/EX10.4/10_4.sce
new file mode 100755
index 000000000..667240355
--- /dev/null
+++ b/564/CH10/EX10.4/10_4.sce
@@ -0,0 +1,9 @@
+pathname=get_absolute_file_path('10_4.sce')
+filename=pathname+filesep()+'10_4data.sci'
+exec(filename)
+for i=1:3
+    f(i)=(0.5*%pi*i/L)*(EI/(rho*A))^0.5;
+end
+printf("\nf1: %f",f(1));
+printf("\nf2: %f",f(2));
+printf("\nf3: %f",f(3));
+\ No newline at end of file
diff --git a/564/CH10/EX10.5/10_5.JPG b/564/CH10/EX10.5/10_5.JPG
new file mode 100755
index 000000000..024539a6d
--- /dev/null
+++ b/564/CH10/EX10.5/10_5.JPG
diff --git a/564/CH10/EX10.5/10_5.sce b/564/CH10/EX10.5/10_5.sce
new file mode 100755
index 000000000..d3e33d0e3
--- /dev/null
+++ b/564/CH10/EX10.5/10_5.sce
@@ -0,0 +1,31 @@
+pathname=get_absolute_file_path('10_5.sce')
+filename=pathname+filesep()+'10_5data.sci'
+exec(filename)
+for i=1:3
+    lambda(i)=(i-0.5)*%pi/L;
+    omega(i)=(lambda(i)^2)*(EI/(rho*A))^0.5;
+    fre(i)=omega(i)/(2*%pi);
+    K(i)=(cos(L*lambda(i))+cosh(L*lambda(i)))/(sin(L*lambda(i))+sinh(L*lambda(i)));
+end
+function[V1]=f(z)
+    V1=K(1)*(cosh(z*lambda(1)) -cos(z*lambda(1)) -K(1)*(sinh(z*lambda(1)) -sin(z*lambda(1))));
+endfunction
+function[V2]=f1(z)
+    V2=K(2)*(cosh(z*lambda(2)) -cos(z*lambda(2)) -K(2)*(sinh(z*lambda(2)) -sin(z*lambda(2))));
+endfunction
+function[V3]=f2(z)
+    V3=K(3)*(cosh(z*lambda(3)) -cos(z*lambda(3)) -K(3)*(sinh(z*lambda(3)) -sin(z*lambda(3))));
+endfunction
+z=linspace(0,L,10*L);
+subplot(3,1,1),xgrid(3),xtitle('1st Normal Mode Of Vibration','z','V')
+V1=feval(z,f);
+plot2d(z,V1);
+subplot(3,1,2),xgrid(3),xtitle('2nd Normal Mode Of Vibration','z','V')
+V2=feval(z,f1);
+plot2d(z,V2);
+subplot(3,1,3),xgrid(3),xtitle('3rd Normal Mode Of Vibration','z','V')
+V3=feval(z,f2);
+plot2d(z,V3);
+printf("\nf1: %f ",fre(1));
+printf("\nf2: %f ",fre(2));
+printf("\nf3: %f ",fre(3));
+\ No newline at end of file
diff --git a/564/CH10/EX10.6/10_6.sce b/564/CH10/EX10.6/10_6.sce
new file mode 100755
index 000000000..313569b49
--- /dev/null
+++ b/564/CH10/EX10.6/10_6.sce
@@ -0,0 +1,5 @@
+pathname=get_absolute_file_path('10_6.sce')
+filename=pathname+filesep()+'10_6data.sci'
+exec(filename)
+omega=1.1584*(EI/(m*l^3));
+printf("\nω1: %f",omega)
+\ No newline at end of file