PSN-L Email List Message

Subject: Re: Shunt damping
From: Bobhelenmcclure@.......
Date: Mon, 9 Oct 2006 21:17:07 EDT


Hi Chris,
 
  For a damped pendulum, we are dealing with the equation
 
  F = m(d2x/dt2) + c(dx/dt) + kx, where F is some forcing  function.
 
 m is mass, k is stiffness, and c is drag. Feynman defines  k=m*Omega0^2, 
c=m*Gamma, and Q=Omega0/Gamma, where Omega0 is the natural  frequency.
 
 In order to calculate the effect of a shunt load on the  pickup coil for 
damping, recall that the sensor output sensitivity is given by S  = B*L, where S 
is in units of volts per meter per second, B is field strength in  Teslas, and 
L is active wire length in meters.
 
  If we pass current through the coil, the resultant motor force is f  = S*I, 
where f is in newtons, and I is amperes.
 
  If we load the moving coil with the resistance of the coil in series  with 
an external resistance, the generated voltage is E = S * v, and the current  
is I = E/R, so the resulting drag is c*v = S*I = S*E/R = S*(S*v)/R.
 
  So, c = S*S/R, where c is drag in newtons per meter per second. From  the 
Feynman definitions, c = m * Gamma = m * Omega0 / Q. Since c from shunt  
damping is S^2 / R, we have the equality S *S / R = m * Omega0 / Q.
 
  Solving for Q, we obtain:
 
Q = R * m * Omega0 / S^2
 
  For a Q of unity, 
 
R = (Ps * S^2) / (2 * Pi * m), where Ps is sensor natural period in  seconds.
 
  When I plug in my assumed values of Ps = 12 seconds, S = 85 v-s/m,  and m = 
0.1 kg, the value for R is 138 kOhms.
 
  In reality, the value I need to get a Q of 1.0 is less  than half the 
above. I have no clue as to why. I calibrated the sensor output  several ways, all 
yielding the same result. The mass was determined by weighing  the solder wire 
mass and the coil assembly. To get better agreement, I would  have to 
considerably decrease my value for S, which incidentally would also  bring my 
calibrated seismic amplitudes into closer agreement with those recorded  by station 
PAL. My present calibration gives me only 40% of the amplitude  recorded at PAL.
 
  I hope someone out there can find a flaw in the mathematics. I am  not 
happy with the apparent disagreement between theory and experiment.
 
Cheers,
 
Bob
 
 





Hi Chris,
 
  For a damped pendulum, we are dealing with the equation
 
  F =3D m(d2x/dt2) + c(dx/dt) + kx, where F is some forcing=20 function.
 
 m is mass, k is stiffness, and c is drag. Feynman defines=20 k=3Dm*Omega0^2, c=3Dm*Gamma, and Q=3DOmega0/Gamma, where Omega0 is the natur= al=20 frequency.
 
 In order to calculate the effect of a shunt load on the= =20 pickup coil for damping, recall that the sensor output sensitivity is given=20= by S=20 =3D B*L, where S is in units of volts per meter per second, B is field stren= gth in=20 Teslas, and L is active wire length in meters.
 
  If we pass current through the coil, the resultant motor force i= s f=20 =3D S*I, where f is in newtons, and I is amperes.
 
  If we load the moving coil with the resistance of the coil in se= ries=20 with an external resistance, the generated voltage is E =3D S * v, and the c= urrent=20 is I =3D E/R, so the resulting drag is c*v =3D S*I =3D S*E/R =3D S*(S*v)/R.<= /DIV>
 
  So, c =3D S*S/R, where c is drag in newtons per meter per second= .. From=20 the Feynman definitions, c =3D m * Gamma =3D m * Omega0 / Q. Since c from sh= unt=20 damping is S^2 / R, we have the equality S *S / R =3D m * Omega0 / Q.
 
  Solving for Q, we obtain:
 
Q =3D R * m * Omega0 / S^2
 
  For a Q of unity,
 
R =3D (Ps * S^2) / (2 * Pi * m), where Ps is sensor natural period in=20 seconds.
 
  When I plug in my assumed values of Ps =3D 12 seconds, S =3D 85=20= v-s/m,=20 and m =3D 0.1 kg, the value for R is 138 kOhms.
 
  In reality, the value I need to get a Q of 1.0 is less= =20 than half the above. I have no clue as to why. I calibrated the sensor outpu= t=20 several ways, all yielding the same result. The mass was determined by weigh= ing=20 the solder wire mass and the coil assembly. To get better agreement, I would= =20 have to considerably decrease my value for S, which incidentally would also=20 bring my calibrated seismic amplitudes into closer agreement with those reco= rded=20 by station PAL. My present calibration gives me only 40% of the amplitude=20 recorded at PAL.
 
  I hope someone out there can find a flaw in the mathematics. I a= m=20 not happy with the apparent disagreement between theory and experiment.
 
Cheers,
 
Bob
 
 

[ Top ] [ Back ] [ Home Page ]