Matlab simulation of motor torque used to control damping
Matlab simulation of motor torque used to control damping
(OP)
I am trying to simulate a spring-mass-damper system which uses an electric generator connected to a resistive load as a variable damper. The damper is varied by controlling the current to the resistive load.
I seem to have a mental deficiency pertaining to the correlation between the damping and power.
The system is forced to oscillate by a sine wave with amplitude = h and period = T. The force from damping is usually reperesented as Force = damping * velocity or F=c*v.
How can I express the power consumed by the damping action in terms of h, T, c and v?
Please assume everything is in phase just to keep this discussion simple.
I seem to have a mental deficiency pertaining to the correlation between the damping and power.
The system is forced to oscillate by a sine wave with amplitude = h and period = T. The force from damping is usually reperesented as Force = damping * velocity or F=c*v.
How can I express the power consumed by the damping action in terms of h, T, c and v?
Please assume everything is in phase just to keep this discussion simple.





RE: Matlab simulation of motor torque used to control damping
Assume the velocity V(t) appears across the damper
V(t) = V0*sin(w*t)
The force is
F = C*V = C*V0*sin(w*t)
Power = F(t)*V(t) = C*V0*sin(w*t) * V0*sin(w*t)
Power(t) = C*V0^2*sin^2(w*t)
<Power> =0.5* C*V0^2
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
<Power> = 0.5*C*V^2 (mechanical) = 0.5*(1/r)*v^2 = vrms^2/r (electrical)
We could also express it as
<Power> = 0.5*F^2/C (mechanical) = 0.5*r*i^2 = irms^2*r (electrical)
Here is the electrical analogy for mechanical vibration, according to electricpete:
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
RE: Matlab simulation of motor torque used to control damping
The next step in this problem is to choose m, k and c to maintain resonance over a range of wave periods by varying the power that is extracted power from the damper. The electrical analogy might help.
I will try to upload a free body diagram for discussion.
RE: Matlab simulation of motor torque used to control damping
This is a valid presumption because the time constant of the generator and circuit is much, much smaller that the wave period. In this case, tenths of second versus 3 to 5 seconds.
RE: Matlab simulation of motor torque used to control damping
The equivalent inertia of the motor, gearbox and sheave is lumped together and called Js.
I will attempt to convert this mechanical system into an electrical equivalent.
Please post an electrical equivalent if anyone is comfortable with this method of analysis.
Thank you in advance.
RE: Matlab simulation of motor torque used to control damping
I am almost inclined to say we should avoid the analogy in the face of these difficulties due to increased potential for confusing what the variables represent. Nevertheless, it can be done.
1 – convert all inertia's to equivalent linear inertia
2 – convert everything to left side of the gears
Here is the electrical analogue circuit I would use:
Ground===Vapplied ===Zk =====Zi === Ground
|
===Zd=====Ground
Where
Vapplied = d/dt(x) = w*H*cos(w*t) / 2
Zk = j*w / k
Zi (inertia) = m + r^2*(Js + n^2*Jm)
Zd (damping) = 1/Ceq
r = radius of shaft upon which you've wound that string
Ceq = is equivalent damping factor seen on left side of gears
I don't know what your bm is supposed to mean but we need a model of your generator.
Let's say your generator voltage is Eg = K' * w' where w' is speed on left side of the gears (w' = w*n if w is speed of generator on right side of gears)
Expressing in terms of velocity on left side v' = w'/r...
Eg = K' * v' * r
Let's say this voltage is applied to resistance R
P = Eg^2 / R = K'^2 * v'^2 * r^2
Ceq = P/v'^2 = K'^2 * r^2
Once you solve the system, the current flowing from Vapplied is the spring force. All other mechanical parameters can also \be solved from the elctrical solution as well although not as easily... not sure which mechanical parameters you are interested in.
I have assumed the "string" that you showed does not stretch. I also assumed the angular velocity of the shaft is low enough that the string always remains tight between shaft and mass (i.e. the downward acceleration of the mass below always remains less than 1 g). Under this assumption we have ideally a linear system with two forcing functions (displacement on top and gravity below) and we divide the solutions by superposition and discard the uninteresting static gravity response.
By the way, what do you mean by "wave period"....what wave?
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
should have been:
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
... but after all of the variables of the mechanical parameters are transfered to the left of the gearbox and combined then the model becomes more straight forward.
The variable 'bm' was supposed to be the motor damping. I should have wrote 'cm'
The string does not stretch and the string always remains tight.
The gear ratio of 1/n is correct. The standard seems to be that a speed reducer has a gear ratio greater that 1 and a speed increaser has a gear ratio less than one.
The wave period that I mentioned refers to the sine wave displacement applied to the free end of the spring.
RE: Matlab simulation of motor torque used to control damping
Ce = c + r^2*(cs + n^2*cm)
In this case I choose to ignore the damping in the rope, rope sheave, and gear box so c = 0 and cs = 0 and then,
Ce = r^2*n^2*cm
RE: Matlab simulation of motor torque used to control damping
The damping transfers throught the gearbox just like the inertia does.
Ce = c + r^2*(cs + 1/n^2*cm)
In this case I choose to ignore the damping in the rope, rope sheave, and gear box so c = 0 and cs = 0 and then,
Ce = r^2*1/n^2*cm
RE: Matlab simulation of motor torque used to control damping
The original goal of this experiment was to maintain resonance by varying Ce. In other words, vary Ce so that the resonant frequency of the mechanical system is the same as the frequency of the displacement that is driving the free end of the spring.
Ideally, the power absorbed by Ce would be 1000W, but a range from 500W to 3000W is acceptable. There are other physical constraints that limit k to a maximum of 1000 kg/m. The equivalent mass, Me, can be any value.
RE: Matlab simulation of motor torque used to control damping
If I recall it right, the resonant frequency does not depend on the damping.
RE: Matlab simulation of motor torque used to control damping
As we know, in a simple series or parallel RLC system, while the undamped natural frequency w0=1/sqrt(LC) does not depend on damping R, the actual damped natural frequency wd does depend on R. You can decrease wd from w0 all the way down to 0 (critical damping) by increasing damping (resistance) for the simple system.
We don't exactly have a simple series or parallel RLC system (we have L in series with paralle R/C) so it will not behave exactly that way same and I do believe R does not hae any influence on wd for this particular circuit. But it can certainly have a big impact on the magnitude of the impedance at a particular frequency.
It leads to an important question - what are we really trying to achieve....what is the purpose of varying resonant frequency? Are we trying to minimize movement of the mass inertia elements?
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
and
Equivalent Free Body Diagram
RE: Matlab simulation of motor torque used to control damping
Outside of ground and the applied displacement / velocity, there is only one nodal position (displacement or velocity) that matters - it is the node between the spring and where the damper / mass connect. The extra node you created between mass and damper has no meaning and would cause an error. The force on the mass is 2nd derivative of position - no reference other than ground and position of the mass is required. The force on the damper is 1st derivative of position (or it's rotational analogue = angle). I imagine you have a rotating magnet and stationary coil which feeds resistor. The stationary coil forms the reference on which the motion is based. Therefore once again the only reference points for damping are the point of connection to the system and ground (ground = stationary reference frame).
Here is the derivation associated with the 6 rules of the analogy above (I should have posted this before):
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
I would like to extract the maximum power out of the system by way of the generator/damper for a displacement wave input on the spring in the range of 3 to 5 second periods.
I think this can be done by keeping the force vector and the velocity vector in the same direction. In other word, make the angle between force and velocity close to zero. In other words, make the damped resonant frequency equal to the frequency of the input dispacement wave.
It is not hard to imagine that if the mass was not oscillating at the same frequency as the displacement, that there would be moments in time where the movement of the mass and damper would diminish to zero and the power extraction would also diminish to zero.
At the damper... Power = Force * Velocity
RE: Matlab simulation of motor torque used to control damping
Thanks again for all of your input.
RE: Matlab simulation of motor torque used to control damping
RE: Matlab simulation of motor torque used to control damping
The bode plot of the transfer function of the spring/mass/damper definitely shows how changing the damping effects the resonant peak of the system. I don't know how to explain this in terms of the electrical equivalent circuit.
Here is a m-file that varies Ce from zero damping to over damping.
RE: Matlab simulation of motor torque used to control damping
I agree also for this system changing damping changes the damped resonant frequencies. I think I was previously working with a slightly different transfer function before I completely understood what are the inputs and ouptuts of your physical problem.
Attached I have done an analysis of the systems using mechanical relationships (part 2) and using the electrical analogy (part 3). I was glad to see the results (transfer function H) are the same with both methods. In part 4 I used the electrical analogy system to determine the "Rmax" which is the value of electrica-analogue parameter R (related to mechanical c) that will maximize power transfer to the generator. In part 5, I show that setting R to maximize power transfer (R=Rmax) does not appear to be equivalent to adjusting R toward resonance. You are welcome to double-check my logic... maybe I missed something.
Also I wonder if the wave excitation is truly independent of the energy extracted... or does the system interact with the waves to change them.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
I have corrected this in the attached. The solution Rmax makes sense only for w below the undamped natural freq of the system. I am not sure what value to choose if w is above that undamped natural freq... would require some more thought and possibly analysis. Probably not today.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
Rmax = w L/ (1 − w2 L C)
and it applies far below w=1/sqrt(LC), not sure about above yet.
One please reality check is to set w very small w<1/sqrt(LC). The the denominator is approx 1 and the numerator approx w*L. For this condition the parallel capacitive impedance is aprox infinite and makes no affect on the circuit. The optimum R is that which matches the supply impedance w*L. That's a good sign that the solution is correct.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
Find the Thevinin equivalent of the circuit supplying R. That is the impedance looking at the rest of the circuit with the voltage source shorted - it is ZC in parallel with ZL.
Zth = ZL*ZC/(ZL+ZC) = j*w*L/(1-w^2*L*C)
The R to maxmimize power transfer has a value equal to the magnitude of Zth
Rmax = w*L/(1-w^2*L*C)
This confirms the analysis for w below undamped resonant frequency which gives postive R. I am not sure what to make of the above-undamped-resonant-frequency case where the denominator turns negative resulting in negative R.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
Rmax = |w*L/(1-w^2*L*C)|
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
We started with mechanical system:
Ground===dapplied==k===meq====ceq====Ground
where dapplied=applied displacement, k=spring, meq=mass, ceq=damping
We converted the input by known relationship d->v:
Ground===vapplied==k===meq====ceq====Ground
where vapplied= applied velocity
We converted to electrical analogy using 6 "rules" above:
Ground===Vapplied ===L ======R === Ground
|
===C=====Ground
where L = 1/k, R = 1/c, C = m, Vapplied = applied voltage
We converted to Thevinin equivalent:
Ground===Vth ===Zth ======R === Ground
Where Zth = j*w*L/(1-w^2*L*C) is impedance from the middle node to ground with Vapplied shorted.
Vth – can be determined by voltage divider to find voltage betwen L and C with R removed, but doesn't matter for purposes of finding the optimum R.
Maximim power transfer occurs when R="Rmax" = |Zth|= |w*L/(1-w^2*L*C)|
This corresponds to c="cmax" = |w/(k-w^2*m)|
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
adjusting damping (R or c) to create max power transfer for a given applied frequency does not lead to the same solution (of R or c) as adjusting damping (R or c) to change the system natural frequency to be equal to the applied frequency.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
Choosing R to maximize I (=current thru resistor) or to maximize V (voltage accross the resistor)
is not the same as
Choosing R to maximize I^2*R
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
If the energy in a wave is to be completely extracted, then the wave is completely destroyed. To destroy a wave is to make a wave which is equal and opposite.
Falnes does a good job explaining the idea in the attached paper. See Section 3 on page 190.
RE: Matlab simulation of motor torque used to control damping
This electrical equivalent stuff seems to be quite useful.
What is your source of information for this technique?
Can you suggest a good book or a few good papers?
Thanks
Paul
RE: Matlab simulation of motor torque used to control damping
My original reference for the vibration / electrical analogy was Harris' Shock and Vib Handbook chapter on "Mechanical Impedance Analysis". He works several problems. The bizarre thing is that he never once mentions anything electrical (not voltage, current, resistance, inductance, or capacitance), even though he uses terms like Thevinin equivalent, Kirchoff's laws etc. Must have been written by a mechanical type
From studying Harris, I wrote down the analogy in terms of the 6 rules (posted 25 Aug 09 11:32) and the derivation (posted 28 Aug 09 17:48). I keep those in a word file for easy reference – for me that captures everything I need to know about it.
I googled "Mechanical Impedance Analysis" and came up with a better reference than Harris (see particularly Table 4.4 which summarizes the analogy):
http:
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
For me the Force Current Analog has good intuitive advantage over the Force Voltage Analogue. In the force current analog, forces flow through branches like current... and relative velocity is a difference between nodes like voltage. The only downside is that the thing called "mechanical impedance" transforms to the inverse of electrical impedance under the Force Current Analog. I just tend to use Z to represent electrical impedance and ignore the conflict with the term "mechanical impedance" (I don't' use that term).
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Matlab simulation of motor torque used to control damping
Thanks for the Maple calculations. Very nice.
I almost finished going throught the calcuations, when my lovely Jeannabelle's water broke and we had to leave for the birthing center. I may have one question about the solution using the equivalent circuit. I will re-examine the calculations after the three of us are back home.
Thanks again for walking me through the electrical equivalent.
Paul