×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Matlab simulation of motor torque used to control damping

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.

RE: Matlab simulation of motor torque used to control damping

I imagine you are talking about a linear generator although the answer should apply to any linear damping element.
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

Whoops. Needs a correction. Let me prepare that...

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

Nope, it was right as written.   It is a perfect illustration of the electrical analogy for mechanical vibration.   The mechanical damper C acts like a electircal resistor r=1/C.  Mechanical velocity V plays the role of electrical voltage v, mechanical force F plays the role of current i.   
<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:

Quote:

You need only 6 simple rules in order to solve any sinusoidal steady state lumped linear m/k/c mechanical vibration problem using electrical circuit impedance (R/L/C) analysis techniques:

1 - mechanical force f plays the role of electrical current I

2 - mechanical velocity difference v plays the role of electr voltage difference V

3 - mechanical mobility  v/f plays the role of electrical impedance Z = V/I
(use the symbol Z throughout.   We use the term impedance to represent either electrical impedance or mechanical mobility .. avoid the term mechanical impedance which is the reciprocal of mobility).

4 - A spring (k) acts like an electrical inductor of inductance L=1/k  
The impedance is Zk = j*w/*k

5 - A damper (c) acts like an electrical resistor of resistance R = 1/c  
The impedance is Zc = 1/c

6 - A mass (m) acts like an electrical capacitor of capacitance C =m  **
The impedance is Zm = 1/(j*w*m)
** one terminal of this mass/capacitance is connected to ground.

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

My answer is based on the presumption that a generator can be used as a linear damping element characterized by it's damping constant C.    That may or may not be the case.   I tend to think it is a close approximation if your coil inductive reactance is much less than total circuit resistance (w*L << R).    

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

(OP)
Yes, Power = Force * Velocity.  Thanks, electricpete.

RE: Matlab simulation of motor torque used to control damping

(OP)
The electrical analogy to mechanical vibration is an interesting approach.

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

(OP)
electricpete said... "My answer is based on the presumption that a generator can be used as a linear damping element characterized by it's damping constant C."

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

(OP)
Here is the FBD for this project.  

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

It is a little more of a challenge to apply and communicate the analogy for this system due to a couple of wrinkes: 1 – combining linear inertia and rotary inertia; 2 – speed transformation ratio; 3 – mechanical damper has physical electrical parameters (load resistance etc)– likely to get mixed up with the analogue electrical parameters.

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

Correction

Quote:

Zi (inertia) = m + r^2*(Js + n^2*Jm)
should have been:

Quote:

Zi (inertia) = 1/(j*w*m_equivalent)
m_equivalent = m + r^2*(Js + n^2*Jm)

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

I think I have inverted the n ratio in at least one place.  

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

Maybe just one place. Should be:

Quote:

Zi (inertia) = 1/(j*w*m_equivalent)
m_equivalent = m + r^2*(Js + (1/n^2)*Jm)

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

Quote:

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 ....
If that assumption is not met, then it is a non-linear system and there is no linear circuit analogy.  I would venture to say that a numerical solution would be the best in that case.  (You don't need numerical solution if the simple linear model applies).  

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

(OP)
It is a little more of a challenge...

... 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

(OP)
The damping transfers throught the gearbox just like the inertia does.

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

(OP)
CORRECTION

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

(OP)
Which mechanical parameters am I interested in?

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

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.

If I recall it right, the resonant frequency does not depend on the damping.

RE: Matlab simulation of motor torque used to control damping

Yes, good point.

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

I would put the damping element in parallel with the mass element as I had drawn on 26 Aug 09 22:46

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):

Quote:

Derivation of the impedances and their electrical equivalencies is as follows:
Spring:
f = k x
f = k (v/jw)
Zk = v/f = j*w / k
By comparison to ZL = j*w*L,  we find L = 1/k

Damping
f = c v
Zc = v/f = 1/C
By comparison to ZR = R,  we find R = 1/c

Mass
f = m a   (where a is referenced to ground)
f = m* (jw * v)   (where v is referenced to ground)
Zm = v/f = 1/ (jw*m)
By comparison to ZC = 1/(j*w*C),  we find C = m
(where C is connected between the location of the mass and ground).

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

(OP)
electricpete said... "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? "

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

(OP)
It has been a busy couple of days for me, so I have only been able to think about this for about a minute and a quarter.  There is a lot of info here for me to digest so I have some catching up to do so we are on the same page.

Thanks again for all of your input.

RE: Matlab simulation of motor torque used to control damping

(OP)
Okay.  I understand the point that you are making about creating an extra node by drawing the damper below the mass.  I agree that the mass and damper should be in parallel to create an electrical equivalent.

RE: Matlab simulation of motor torque used to control damping

(OP)

Quote:

electricpete said "We don't exactly have a simple series or parallel RLC system (we have L in series with parallel R/C) so it will not behave exactly that way same and I do believe R does not have any influence on wd for this particular circuit."

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 can see why you drew the FBD the way you did - in the mechanical world, it is understood a lumped mass has only one position.  The position of the "input" connection to the mass is same as the position of the "output" connection to the mass.   Only in the electrical analogy when we substitute a capacitor do we find  two distinct terminals with different voltages (related to velocity and displacement) and there would be a need to draw the Cap (representing mass) in parallel with the R (represienting damping).  So the choice of how to connect that element depends on whether you are thinking electrical or mechanical.   To use the electrical analogy it needs to be drawn as I described.  If you are viewing it in the mechanical world as a mass, your FBD was fine.  So I withdraw my comment.

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

Above I had a polarity error inside the sqrt in the expression for wde (wd from electrical analogy).

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

My solution for Rmax was
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

It could have been solved much easier I think.
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

Duh, sorry, the magnitude of Z never turns negative.  I think we should just use the magnitude of thevinin impedance which is always positive:
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

To summarize:

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

One last point  I should have included in my summary:
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

.... the last point said another way:
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

(OP)

Quote:

electricpete said, "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."

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

(OP)
electricpete,

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

Thanks for the paper.

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 winky smile

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://books.google.com/books?id=Yy_8RNlD2ngC&pg=PT242&lpg=PT242&amp;dq=mechanical+impedance+analysis&source=bl&amp;ots=7G31baaUGe&amp;sig=mFGxX25oqiAJAYRB-j7c2ZlXcN0&hl=en&ei=9rqbSqu4LM7JlAefoOW9BQ&sa=X&oi=book_result&ct=result&resnum=4#v=onepage&amp;q=mechanical%20impedance%20analysis&amp;f=false
 

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Matlab simulation of motor torque used to control damping

I should have mentioned that the analogy I used is referred to in Table 4.4 of above link as "Force Current Analog".   Not to be confused with "Force Voltage Analog".
 
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

(OP)
electricpete,

Thanks for the Maple calculations.  Very nice. smile

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

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources