1 files changed, 21 insertions, 0 deletions
diff --git a/3532/CH4/EX4.1/Ex4_1.sce b/3532/CH4/EX4.1/Ex4_1.sce
new file mode 100644
index 000000000..8dfd86913
--- /dev/null
+++ b/3532/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,21 @@
+clc
+clear
+mprintf('Mechanical vibrations by G.K.Grover\n Example 4.2.1\n')
+//given data
+//T=To*sin(W*t)
+To=0.588 //maximum value of periodic torque in N-m 
+W=4// freqency of applied force in rad/sec
+J=0.12//moment of inertia of wheel in kg-m^2
+Kt=1.176//stiffness of wire in N-m/rad
+Ct=0.392/1 //damping coefficient in N-m_sec/rad
+//calculations
+theta=To/sqrt((Kt-J*W^2)^2+(Ct*W)^2)//Equation for torsional vibration amplitude from Fig (4.2.2) and Eqn (4.2.5)
+MaxDcoup=Ct*W*theta//maximum damping couple in N-m
+if atan((Ct*W)/(Kt-J*W^2))>0  then
+    phiD=(180/%pi)*atan((Ct*W)/(Kt-J*W^2));//from eqn 4.2.6(in degrees)
+else
+    phiD=180+(180/%pi)*atan((Ct*W)/(Kt-J*W^2));
+    
+end
+//output
+mprintf(' a)The maximum angular displacement from rest position is %4.4f radians\n b)The maximum couple applied to dashpot is %4.4f N-m\n c)angle by which the angular displacement lags the torque is %4.4f degrees',theta,MaxDcoup,phiD)