8 files changed, 104 insertions, 0 deletions
diff --git a/3588/CH10/EX10.1/EX10_1.sav b/3588/CH10/EX10.1/EX10_1.sav
new file mode 100644
index 000000000..b44b19fee
--- /dev/null
+++ b/3588/CH10/EX10.1/EX10_1.sav
diff --git a/3588/CH10/EX10.1/EX10_1.sce b/3588/CH10/EX10.1/EX10_1.sce
new file mode 100644
index 000000000..aa12a0a1d
--- /dev/null
+++ b/3588/CH10/EX10.1/EX10_1.sce
@@ -0,0 +1,21 @@
+//Clearing console
+clc
+clear
+
+//Intializing variables
+mg = 20
+k = 25
+g = 386.4
+
+//Calculating circular frequency
+w = sqrt(k*g/mg)
+
+//Solving for constants in equation of motion
+fi = acosd(0/w)
+C = (2.3-0.8)/sind(fi)
+
+//Circular frequency in Hz
+f = w/(2*%pi)
+
+printf('\nResults\n')
+printf('\nCircular Frequency =%fHz\n Amplitude =%fin\n Phase Angle =%fdegree',f,C,fi)
diff --git a/3588/CH10/EX10.2/EX10_2.sav b/3588/CH10/EX10.2/EX10_2.sav
new file mode 100644
index 000000000..e442c23f6
--- /dev/null
+++ b/3588/CH10/EX10.2/EX10_2.sav
diff --git a/3588/CH10/EX10.2/EX10_2.sce b/3588/CH10/EX10.2/EX10_2.sce
new file mode 100644
index 000000000..98cd4ce25
--- /dev/null
+++ b/3588/CH10/EX10.2/EX10_2.sce
@@ -0,0 +1,31 @@
+//Clearing console
+clc
+clear
+
+//Intializing variables
+mg = 20
+k = 40
+g = 386.4
+
+//Calculating circular frequency
+w1 = sqrt(k*g/mg)
+w2 = sqrt(6*k*g/mg)
+
+k1 = [1 1;2 -0.5]
+f1 = [1;0.5]
+k2 = [27.8 68.1;2*27.8 0.5*68.1]
+f2 = [0;0]
+//Solving for constants in equations of motion
+u1=linsolve(k1,-f1)
+u2=linsolve(k2,-f2)
+
+fi1 = acosd(u2(1,1))
+fi2 = acosd(u2(2,1))
+
+a1 = u1(1,1)/sind(fi1)
+a2 = u1(2,1)/sind(fi2)
+
+printf('\nResults\n')
+printf('\nCircular Frequency1 =%frad/sec\n Amplitude1 =%fin\n Phase Angle1 =%fdegree',w1,a1,fi1)
+printf('\nCircular Frequency2 =%frad/sec\n Amplitude2 =%fin\n Phase Angle2 =%fdegree',w2,a2,fi2)
+
diff --git a/3588/CH10/EX10.6/EX10_6.sav b/3588/CH10/EX10.6/EX10_6.sav
new file mode 100644
index 000000000..022d3eec5
--- /dev/null
+++ b/3588/CH10/EX10.6/EX10_6.sav
diff --git a/3588/CH10/EX10.6/EX10_6.sce b/3588/CH10/EX10.6/EX10_6.sce
new file mode 100644
index 000000000..647ba56d5
--- /dev/null
+++ b/3588/CH10/EX10.6/EX10_6.sce
@@ -0,0 +1,35 @@
+//Clearing console
+clc
+clear
+
+//Intializing variables
+t = 5
+p = 7.83*10^-6
+x0 = -1
+x1 = 1
+
+//Calculating elements of mass matrix
+m11 = (150*p*t*integrate('(1-r)^2','r',x0,x1)*integrate('(1-s)^2','s',x0,x1))/16
+m12 = (150*p*t*integrate('1-r^2','r',x0,x1)*integrate('(1-s)^2','s',x0,x1))/16
+m22 = (150*p*t*integrate('(1+r)^2','r',x0,x1)*integrate('(1-s)^2','s',x0,x1))/16
+m13 = (150*p*t*integrate('1-r^2','r',x0,x1)*integrate('1-s^2','s',x0,x1))/16
+m14 = (150*p*t*integrate('(1-r)^2','r',x0,x1)*integrate('1-s^2','s',x0,x1))/16
+m23 = (150*p*t*integrate('(1+r)^2','r',x0,x1)*integrate('1-s^2','s',x0,x1))/16
+m24 = (150*p*t*integrate('1-r^2','r',x0,x1)*integrate('1-s^2','s',x0,x1))/16
+m33 = (150*p*t*integrate('(1+r)^2','r',x0,x1)*integrate('(1+s)^2','s',x0,x1))/16
+m34 = (150*p*t*integrate('1-r^2','r',x0,x1)*integrate('(1+s)^2','s',x0,x1))/16
+m44 = (150*p*t*integrate('(1-r)^2','r',x0,x1)*integrate('(1+s)^2','s',x0,x1))/16
+
+//Constructing Mass matrix
+m(1,1:8) = [m11 m12 m13 m14 0 0 0 0]
+m(2,1:8) = [m12 m22 m23 m24 0 0 0 0]
+m(3,1:8) = [m13 m23 m33 m34 0 0 0 0]
+m(4,1:8) = [m14 m24 m34 m44 0 0 0 0]
+m(5,1:8) = [0 0 0 0 m11 m12 m13 m14]
+m(6,1:8) = [0 0 0 0 m12 m22 m23 m24]
+m(7,1:8) = [0 0 0 0 m13 m23 m33 m34]
+m(8,1:8) = [0 0 0 0 m14 m24 m34 m44]
+
+printf('\nResults\n')
+printf('\nMass matrix m in (kg)')
+disp(m)
diff --git a/3588/CH10/EX10.9/EX10_9.sav b/3588/CH10/EX10.9/EX10_9.sav
new file mode 100644
index 000000000..5e3e6b0fc
--- /dev/null
+++ b/3588/CH10/EX10.9/EX10_9.sav
diff --git a/3588/CH10/EX10.9/EX10_9.sce b/3588/CH10/EX10.9/EX10_9.sce
new file mode 100644
index 000000000..747322fd4
--- /dev/null
+++ b/3588/CH10/EX10.9/EX10_9.sce
@@ -0,0 +1,17 @@
+//Clearing console
+clc
+clear
+
+//Intializing variables
+t1 = 0.03
+w1 = 5
+t2 = 0.1
+w2 = 15
+
+//Solving for Rayleigh coefficients
+k = [1/(2*w1) w1/2;1/(2*w2) w2/2]
+f = [t1;t2]
+u=linsolve(k,-f)
+
+printf('\nResults\n')
+printf('\nRayleigh coefficients\n Alpha =%f\n Beta =%f',u(1,1),u(2,1))