diff options
Diffstat (limited to '3532/CH6')
| -rw-r--r-- | 3532/CH6/EX6.15/Ex6_15.sce | 22 | ||||
| -rw-r--r-- | 3532/CH6/EX6.8.2/Ex6_15.sce | 22 |
2 files changed, 44 insertions, 0 deletions
diff --git a/3532/CH6/EX6.15/Ex6_15.sce b/3532/CH6/EX6.15/Ex6_15.sce new file mode 100644 index 000000000..0c16e4d46 --- /dev/null +++ b/3532/CH6/EX6.15/Ex6_15.sce @@ -0,0 +1,22 @@ +clc
+clear
+mprintf('Mechanical vibrations by G.K.Grover\n Example 6.8.2\n')
+//given data
+m1=250;m2=100//mass of two blocks in Kgs
+c1=80;c2=60,c=20//damping coefficients in N-sec/m
+F1=1000;F2=1500//amplitude of force acting on block 1 and 2 rsptly
+k=250000//stiffness of spring in N/m
+W=60//frequency of applied force in rad/sec
+//calculations
+M=[m1,0;0,m2];
+K=[k,-k;-k,k];
+C=[c+c1,-c;-c,c+c2];
+R=[F1;F2;0;0];
+X=K-(W^2)*M
+Y=W*C
+G=[X,-Y;Y,X]
+AB=inv(G) *R//from Eqn6.8.4 in Sec 6.8
+X1=sqrt(AB(1,1)^2 +AB(3,1)^2)
+X2=sqrt(AB(2,1)^2 +AB(4,1)^2)
+//output
+mprintf('The amplitude of vibrations are %fm for mass 1 and %fm for mass 2',X1,X2)
diff --git a/3532/CH6/EX6.8.2/Ex6_15.sce b/3532/CH6/EX6.8.2/Ex6_15.sce new file mode 100644 index 000000000..0c16e4d46 --- /dev/null +++ b/3532/CH6/EX6.8.2/Ex6_15.sce @@ -0,0 +1,22 @@ +clc
+clear
+mprintf('Mechanical vibrations by G.K.Grover\n Example 6.8.2\n')
+//given data
+m1=250;m2=100//mass of two blocks in Kgs
+c1=80;c2=60,c=20//damping coefficients in N-sec/m
+F1=1000;F2=1500//amplitude of force acting on block 1 and 2 rsptly
+k=250000//stiffness of spring in N/m
+W=60//frequency of applied force in rad/sec
+//calculations
+M=[m1,0;0,m2];
+K=[k,-k;-k,k];
+C=[c+c1,-c;-c,c+c2];
+R=[F1;F2;0;0];
+X=K-(W^2)*M
+Y=W*C
+G=[X,-Y;Y,X]
+AB=inv(G) *R//from Eqn6.8.4 in Sec 6.8
+X1=sqrt(AB(1,1)^2 +AB(3,1)^2)
+X2=sqrt(AB(2,1)^2 +AB(4,1)^2)
+//output
+mprintf('The amplitude of vibrations are %fm for mass 1 and %fm for mass 2',X1,X2)
|