×
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

Induction motor modelling.

Induction motor modelling.

Induction motor modelling.

(OP)
Hi everyone.

I need some help with modelling an induction motor. All i want to do is to get a model of an induction motor which can simulate the starting time(Time it takes to get to the full speed) with no load attached to it. I tried to use matlab simulink but there is a problem. I have the motor parameters like R and X from the manufacturer. They have given me 2 sets of data(Data Attached). I dont know why they have given me 2 sets for full load and Locked rotor. But when i try these values i get two different starting times. I believe its giving the deep bar effect of the motor. But i want to have an idea from you guys. please help.

gokul

RE: Induction motor modelling.

I only see 1 set of data.  It covers full load and starting.

RE: Induction motor modelling.

(OP)
ya i was saying two sets for Full load and starting only.

My question is R and X should be one value right.

I dont know if i am wrong. But i think that how can the same winding can have different values while starting and at full load.

  

RE: Induction motor modelling.

I would say mabye R1 changed due to temperature (assume higher temperature running than starting).

R2/s of course changed due to change in slip.

Varying levels of sSaturation of certain areas may affect X1, X2, Xm.

Deep bar effect may play some role as well for R2 and maybe X2.

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

RE: Induction motor modelling.

(OP)
Ya I agree that slip changes. But when we calculate with the corresponding slip value too R2 is differing.

Major doubts i have is that i dont know if this motor has a rotor which accounts for deep bar effect.

And suppose it provides deep bar effect then it should have a higher resistance and reactance at the starting right. But for my agony :) it is in the other way around. Please give me some reason for that.

thanks
gokul

RE: Induction motor modelling.

This does in fact show higher rotor resistance at starting as expected by deep bar effect if you account for slip.

Locked rotor:
R2/s = 0.0769
s=1
R2 = R2/s * s = 0.0769

Running
R2/s = 2.8292
s = (1800-1785)/1800 = 0.01
R2 = s* R2/s = 0.0283

R2 higher during locked rotor due to deep bar effect.
 

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

RE: Induction motor modelling.

(OP)
ok. Thank you for that clarification. Lets take that that motor accounts for deep bar effect.

Now my question is why is then the reactances x1 adn x2 are greater when they are at full load? They are also kind of resistances and hence should be high at starting and less at running correct.

This is where this issue confuses me

RE: Induction motor modelling.

As I said the reactances are a little trickier.

Here is my take fwiw.  

The leakage reactance paths saturate during start due to high current.  For example in closed rotor slots that would be the bridge section.  In open or semi-open slots, the tooth and tooth ear corners.    This is the iron part where flux the flux that encircles the conductor (leakage flux) IS constricted (high flux density).  Saturation of leakage reactance during start means the effective leakage inductances are lower during start.

On the other hand, the airgap flux is higher during run than start and the toothtops tend to move more into saturation.  That will make running magnetizing reactance slightly higher.
 

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

RE: Induction motor modelling.

Two corrections in bold below.

Quote:

The leakage reactance paths saturate during start due to high current.  For example in closed rotor slots that would be the bridge section.  In open or semi-open slots, the tooth ear corners.    This is the iron part where flux the flux that encircles the conductor (leakage flux) IS constricted (high flux density).  Saturation of leakage reactance during start means the effective leakage inductances are lower during start.

On the other hand, the airgap flux is higher during run than start and the toothtops tend to move more into saturation.  That will make running magnetizing reactance slightly lower.
The two corrections account for the following:

1 - The area that saturates for semi-open or open slots is tooth ear corners (not teeth in general).
2 - Whichever one saturates gets lower effective reactance.  So the running magnetizing reactance is lower (since it is more saturated.
 

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

RE: Induction motor modelling.

(OP)
Thank you so much for that heads up. Do you have some specific literature explaining these things online so that i can go through them and get known. Thanks.

gokul

RE: Induction motor modelling.

Gokul..Typically to calculate the starting time, the driving torque & driven torque are needed. Do you have the Moment of Inertia of the motor only?

RE: Induction motor modelling.

cherry2000,
Please take a close look at the attachement:
Rotor Inertia = 87 lbs- ft2

RE: Induction motor modelling.

One more data required for calculating the starting time is the Motor speed torque characteristic. The no load curve also will be required. If both can be plotted together,the starting time can be calculated by deriving the plotted values on an excel sheet and using some fundamental formulae

RE: Induction motor modelling.

the  one set of data in locked rotor when  s= 1
r2 vary with s in locked rotor or at start r2 is minimum
the secend set is in full load  when s=sr=ns-nr/ns=min
r2 =max  

RE: Induction motor modelling.

You want time to start an unloaded motor... all it requires is torque speed curve and inertia.

Start from first principles:
dw/dt = (Telec-Tmech)/J     
dt = [J/(Telec-Tmech)] dw    
time =  [J/(Telec-Tmech)] dw, w=0..wfinal    

For N in rpm, J in lbm-ft^2, T in ft-lbf, time in seconds, we can express the above equation as follows :    
time =  [J/<(Telec-Tmech)*307>] dN, N=0..Nfinal    
The constant 307 = 60*32.2/(2Pi)
Of course in your case Tmech=0.

The integration can be carried out numerically.  A long time ago I created a a spreadsheet that does the job.  Requires reading off of the torque speed curve.  I have plugged in coarsely-spaced data from your torque speed curve along with inertia = 87 lbm=ft^2 mentioned above.   I get 0.54 seconds acceleration time at full voltage and 0.846 seconds at 80% voltage.

You will certainly want to check the calculation carefully for yourself to make sure I have correctly represented your problem and make sure you agree with the calcualtion approach and make sure I didn't make any mistakes.

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

RE: Induction motor modelling.

I left out the integral sign in expression for time - I think you can figure out where it goes.

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

RE: Induction motor modelling.

The endpoint of the integration deserves a little more discussion.    You actually cannot integrate all the way to the point of 0 accelerating torque  since it would require smaller and smaller step size and the time to reach 0 accelerating torque is theoretically infinite (integral 1/x, x=0..x0  is infinite) even though for all measurable and practical purposes the speed becomes steady.  

To get around this, we need to select an endpoint which is somewhere short of the 0-accelerating torque point (for example within some tolerance).   In the case of no-load motor full load start, assuming we are looking to size relays,  a convenient place to stop the integration might be the speed corresponding to full-load torque, since current will drop to FLA at this point.   So to be more precise I should have divided the interval beyond breakdown torque into more slices and stopped at the speed corresponding to full load and stated it was the time to accelerate to full load speed.  If you perform that excercise, I don't think you will see much change in the result.
 

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

RE: Induction motor modelling.

(OP)
Thank you so much for your help. Sorry for the late reply. Somehow the notification went to the spam. I went through your spreadsheet. I have two questions for you on the spreadsheet.

Why do you find the average of the torque from one speed to the next?

Then, you do you multiply torque with voltage square?

Once again thank you so much for the help. It was very helpful.

gokul

RE: Induction motor modelling.

Quote:

Why do you find the average of the torque from one speed to the next?
It is the trapezoidal rule for integration.   To integrate we need to find the area under the curve.   Column B is snapshots of the curve at specific speeds. The trapezoidal rule is equivalent to drawing a straight line in between each data point and finding the area under that line.  Then the average value over the interval is simply the average of the endpoints.

Quote:

Then, you do you multiply torque with voltage square?
Torque is proportional to voltage squared.  We can see this from the equivalent circuit.   The power output is given by the power dissipated in the element with resistance
R2*(1-s)/s.   Pick a given speed and keep it fixed (s is fixed also). If we multiply the input voltage by 0.9, then every voltage in the system changes by 0.9.  The voltage accross the resistance R2*(1-s)/s will change by 0.9.  The power consumed by that element (delivered to the load) changes according to V^2/R (where V is voltage accross that element and R = R2*(1-s)/s).  So the output power will go to 0.9^2 = 0.81.  Since the speed is fixed, the torque (Torque = P/w) will also increase by 0.l9^2 =0.81.  So it is seen from the equivalent circuit that torque varies with voltage squared.

I seem to recall for some motors the real exponent is closer to 1.8, but I'm not sure about that last part.  Does anyone use 1.8 instead of 2.0?

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

RE: Induction motor modelling.

typographical correction:

Quote (electricpete):

Since the speed is fixed, the torque (Torque = P/w) will also decrease by 0.9^2 =0.81

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

RE: Induction motor modelling.

(OP)
Thank you so much for the information.

Actually i am trying to theoretically prove a solution to a problem which i found out using trial and error basis.

It is for this project the induction motor time is going to be used.

I would like to discuss the scenario with you.

We have a conveyor which is run by a single drive of 250HP.

This drive is coupled to the load by a Fluid coupling.

Also, This drive is started using a soft starter.

Problem this conveyor was giving to us was that it didnt start when it was started with load on it. But it would start if there is no load on it and it would run with load.

So we figured that it was a starting torque issue.

Now after discussing this issue long back with this forum and having got a lot of inputs from here we figured that we need to have the fluid coupling couple the load to the motor exactly at the maximum torque level.

So we increased the starting voltage applied to the motor and tried to give more time delay to the fluid coupling.

Before this problem the soft starter was set to 4xFLA.

We went to 5xFLA and to the smallest plug possible to give time delay of 5 seconds.

It then started the conveyor with no prob.

Now i want to prove this theoritically.

That is why i am trying to simulate this motor time, to show that at 4xFLA, the time to start the motor is more and the fluid coupling starts coupling the load soon.

So when we increase the voltage the more gets to full speed soon so that it will give the max torque when the load is applied.

i will let you know the outcome after i come with the results.

Let me know if you find any issues with my proof.

thanks


gokul



 

RE: Induction motor modelling.

Beats me.  We don't have any fluid couplings or soft-starters at my plant.  Maybe someone else can comment on your scenario.

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

RE: Induction motor modelling.

(OP)
Hi electricpete,

I have another question for you. Is there a way i can incorporate the effect of mechanical load added to the motor. I want to increase the Tm(to a value equal to the rated torque of the motor) linearly as the motor start  and calculate the time taken by the motor to reach the full speed. Can you help me on this please.

gokul

RE: Induction motor modelling.

Attached is an analysis of your motor (based on data provided above) accelerating a load whose torque varies linearly from 0 to 735.6 ft-lbf (rated motor torque) over a speed range from 0 to 1785 (rated motor speed).

The summary is shown on the graphic tab (speed vs time).

The important change made in the calculations tab was simply to use Taccelerating = Tmotor-Tload for computing acceleration, in place of Tmotor that was used in the previous unloaded start case.

There were some housekeeping changes required as well. To develop a good plot I changed the averaging interval from "this-speed-to-the-next" to "last-speed-to-this".  That change was a prerequisite to providing a traditional plot since the logical plotting point involves the endpoint of the speed interval rather than the beginning.

As you may notice the speed approaches 1785 rpm asymptotically.  We reached 1784 rpm in 0.95 seconds.  We cannot continue the numerical solution to 1785 because we always need smaller and smaller intervals since the rate of increase is slowing.   

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

RE: Induction motor modelling.

(OP)
Thank you electricpete. I guess i am closing in to a close good solution. I am talking to the fluid coupling guys to get their input. Possibly will figure out the solution. Will let you know the result once i get all stuffs covered. Thank you once again for the sheet.

gokul

RE: Induction motor modelling.

One thing to mention is that I used the motor inertia for the calculation.  I don't think you told me any load inertia... I assumed it was small in comparision to the motor inertia (but we all know what they say about assmptions)

If not, the times reported above should be adjusted proportionally to the ratio of total inertia/motor inertia. i.e.
TimeToReach1784rpm=0.95sec*(MotorInertia+LoadInerta)/MotorInertia

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

RE: Induction motor modelling.

My integration method should converge to the right solution if step size is small enough.  But it is a little clumsy since we are assuming:
Integral (1/Taccel) = 1/Tavg where Tavg is the average value of the torque over the interval.  It works well when the change in torque over the interval is a small fraction of the torque.  But not so well when the torque is changing.

I think I will try a new approach.  On the interval N1 to N2, torque is given by a linear curve:
T = T1 + (N-N1)/(N2-N1) * (T2-T1)

1/T = 1/[T1 + (N-N1)/(N2-N1) * (T2-T1)]

We can integrate the above expression with respect to N over the interval N1 to N2 and should get a result in terms of ln function (since the variable of integration N appears in the denominator).  I will do that and post a new spreadsheet shortly.

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

RE: Induction motor modelling.

(OP)
sounds great. Thanks for that.

gokul

RE: Induction motor modelling.

Attached is a revised spreadsheet as discussed above.  Now the time to get to 1784rpm is 0.89 seconds.

This approach I'm sure is more accurate than my previous approach and also does not require the user to enter as many points in areas where the torque curve changes by a significant fraction  (i.e where it approaches 0).  This numerical integration is accurate to the extent that the torque changes linearly over the speed range of interest.  Thus it does very well in the previously-problematic speed range adjacent to 1785rpm.... since torque changes linearly in this range.

Derivation is shown in the derivation tab.  The resulting formula for time dT12 to accelerate thru a given  speed range N1..N2 with torque T1 and T2 on the endpoints is:
dT12 = J*delta/(307*(T2-T1)) * [ln(T2/T1)]
where delta = N2-N1

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

RE: Induction motor modelling.

One more comment - if your motor is driving a centrifugal pump or fan, the load torque is proportional to speed^2.  It would be very easy to modify the spreadsheet to cover that.

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

RE: Induction motor modelling.

(OP)
I have a comment to make. Now when i correct the voltage for reduced voltage ( say .8V) then at some point the load torque becomes more than the motor torque and the resultant goes negative. What should i do then, jus make another column which gives the abs value of the Tacc and proceed. Is that a correct way to do?

gokul

RE: Induction motor modelling.

No.   The point where the accelerating torque crosses 0 (in direction from positive to negative) is the intersection of the torque speed curve of the motor and load - this is the operating point. Speed will not go higher than this and only the time leading up to this crossover has any physical significance.  In the case of parameters you provided where load torque was equal to motor nameplate torque at motor nameplate speed, if we reduce voltage below nameplate then the motor will not make it to nameplate speed, it will be somewhere below nameplate and the load will be somewhere above nameplate (overloaded).

I re-did the simulation using V=0.8 for the data.  It is easy to find the point where the accelerating torque turns negative because we get #num (errors) due negative number input to the ln function.  I filled in a higher point density just before it turns negative to plot the curve leading up to the operating point.  The results are shown in attached.  The operating speed is 1777rpm and it takes 3 seconds to get there.  The motor is loaded to roughly =(1800-1777)/(1800-1785) = 153% and would probably trip on overload not too long after start.

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

RE: Induction motor modelling.

Sorry, forget the part about being overloaded, that was incorrect..  The HP drawn would be less than nameplate if the speed is less than nameplate.  And also if we really wanted to determine the operating point we need the operating voltage instead of starting voltage... operating voltage would be higher than starting resulting in higher operating speed.  My spreadsheet is not equipped to predict changes in voltage (depends on the system impedance and current drawn by the motor).

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

RE: Induction motor modelling.

I agree with electricpete of course.
The equivalent circuit is the extension of the Steinmetz model for transformer.
Using this model one will get about 10 to 25% errors.
The maximum[break down ,pull out] torque does not depends on rotor resistance but the rotor slip corresponding to this maximum torque does ,then one may increase the rotor resistance in order to increase the starting torque up to maximum.
As in deep slot bar the resistance is elevated at start the deep slot is used -in most of induction motors- in order to improve the starting conditions: limiting the starting [lock out] current and increase the starting torque.
In this case the resistance and the reactance of the rotor are related to the slip in a nonlinear dependence. The conventional model is not applicable entirely.
Some people uses two series of parameters one for start and other for rated speed.
Other problem is the excitation current which depends on supply voltage and emf produced in stator windings and it is strong influenced by core saturation. Emf contains then saturation harmonics but also space harmonics and slot harmonics.
 

RE: Induction motor modelling.

(OP)
Hi electripete,

I have another question. When i was going through the calculation i found you were using 32.2 constant to arrive at 307. what is this factor 32.2 account for.

gokul

RE: Induction motor modelling.

Hi gokulkrish2!

The short answer:  32.2 is a unit conversion associated with British/English/American units.  Some people carry around gc but I don't.  The way I look at it  is simply 1 = [lbf  sec^2 / (32.2  lbm ft)]  (easy to remember since F=MA  =>  1=F/MA and since we know that  a force on one lbf occurs when a mass of 1 lbm is subjected to acceleration of gravity which is 32.2 ft/sec^2).   

The long answer is to repeat the full derivation from the beginning, this time accounting explicitly for the units.     Here it goes:
Let T(w) = accelerating torque = Telec – Tmech
The rotational equivalent of F = m*a is:
T(w) = J * dw/dt

Rearrange:
dw/dt = T(w)/J
dt = [J/T(w)] dw
t =  Int [J/(T(w)] dw    [equation 1]

Define the right hand side of equation 1 as "Integral 1"
Integral 1 = Int [J/(T(w)] dw

We know speed N = w/(2*pi)
dN = dw/(2*pi)
Substitute dw = (2*pi)*dN
Integral 1 = Int [J/(T(w)] dw = 2 * Pi * Int [J/(T(w)] dN

Multiply by 1 = [lbf  *sec^2 / (32.2  lbm *ft)]  
Integral 1 = 2 * Pi * Int [J/(T(w)] dN    * lbf  *sec^2 / (32.2  lbm *ft)

Multiply by ft/ft
Integral 1 = 2 * Pi * Int [J/(T(w)] dN    * ft-lbf  *sec^2 / (32.2  lbm *ft^2)

Rearrange:
Integral 1 = 2 * Pi * Int [{/J/<lbm-ft^2>} (T/<ft-lbf>) )] dN   * sec^2 / 32.2

multiply by minute/<60*sec>
Integral 1 = 2 * Pi * Int [{/J/lbm-ft^2} (T/ft-lbf) )] dN * minute   * sec / (32.2 *60)

Substitute minute = 1/min^-1
Integral 1 = 2 * Pi * Int [{/J/lbm-ft^2} (T/ft-lbf) )] d{N/min^-1} * sec / (32.2 *60)

Put Integral 1 back into equation 1:
t= 2 * Pi * Int [{/J/lbm-ft^2} (T/ft-lbf) )] d{N/min^-1} * sec / (32.2 *60)

Divide both sides by sec
{t/sec}= 2 * Pi * Int [{/J/lbm-ft^2} (T/ft-lbf) )] d{N/min^-1}  / (32.2 *60)

Collect Constants at the end
{t/sec}=  Int [{/J/lbm-ft^2} (T/ft-lbf) )] d{N/min^-1}  *[2*Pi]/ (32.2 * 60)

Now we only need recognize that t/sec is a unitless quantity we would call "time in seconds".   J/lbm-ft^2 is a unitless quantity we would call "rotating inertia in lbm-ft^2".  T/ft-lbf is a unitless quantity we can call "Torque in ft-lbf".   N/min^-1 is a unitless quantity we can call speed in rpm.

Based on the above we can rewrite in a more compact form:
t =  Int [J/T)] dN * (2*pi) / (32.2 * 60)
t =  Int [J/T)] dN *  / 307     [equation 2]
where t in seconds, J in lbm-ft^2,  T in lbf-ft,  N in rpm.

Now look at some small interval from speed N1 to N2 over which we will perform numerical integration.  From here to the end we will ASSUME the torque changes linearly with speed over each interval as given by:    
    T = T1 + (N-N1)/(N2-N1) * (T2-T1)
Let delta = N2-N1 be the width of the interval
    T = T1 + (N-N1)/ delta * (T2-T1)
Call this linear function f(N) = T1 + (N-N1)/ delta * (T2-T1)    

Let dt12 be the time to get from N1 to N2            
dt12 = J / 307 * Int ( 1/[f(N) ), N=N1..N2  from equation 2.
    Evaluate integral by chain rule:        
dt12 = J/307 * [ln(f(N) / df/dN]  evaluated at N2 minus the same thing evaluated at N1 (provided that df/dN is not itself a function of N)            
    df/dN = (T2-T1)/delta    (we now confimred that df/dN is not a function of N)        
dt12 = J*delta/(307*(T2-T1)) * [ln(f(N)]  evaluated at N2 minus the same thing evaluated at N1            
    But we know f(N) evaluated at N2 and N1... it is T2 and T1!        
dt12 = J*delta/(307*(T2-T1)) * [ln(T2) -ln(T1)]             
    
    Simplify by noting ln(T2)-ln(T1) is ln(T2/T1)        
dt12 = J*delta/(307*(T2-T1)) * [ln(T2/T1)]       [equation 3]        
where
dt12 - time to accelerate through speed range (in seconds)
J  - rotating inertia (in lbm-ft^2)
T1, T2  = Torque at beginning and end of speed range (in lbf-ft)
delta = width of the speed range (in rpm)

Should we prefer to avoid British units, instead of having unit factor the factor 307 =(32.2 * 60)/(2*pi) in the denominator, we would not need the 32.2 and instead have a factor (60)/(2*pi) = 9.55 in the denominator.

dt12 = J*delta/(9.55*(T2-T1)) * [ln(T2/T1)]             
where
dt12 - time to accelerate through speed range (in seconds)
J  - rotating inertia (in kg-m^2)
T1, T2  = Torque at beginning and end of speed range (in N-m)
delta = width of the speed range (in rpm)
 

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

RE: Induction motor modelling.

(OP)
Thank you so much. I got that.

gokul

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