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
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.
RE: Induction motor modelling.
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.
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.
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.
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.
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.
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.
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.
gokul
RE: Induction motor modelling.
RE: Induction motor modelling.
Please take a close look at the attachement:
Rotor Inertia = 87 lbs- ft2
RE: Induction motor modelling.
RE: Induction motor modelling.
gokul
RE: Induction motor modelling.
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.
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.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Induction motor modelling.
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.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Induction motor modelling.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Induction motor modelling.
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.
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.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Induction motor modelling.
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.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Induction motor modelling.
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.
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.
gokul
RE: Induction motor modelling.
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.
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.
gokul
RE: Induction motor modelling.
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.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Induction motor modelling.
gokul
RE: Induction motor modelling.
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.
=====================================
Eng-tips forums: The best place on the web for engineering discussions.
RE: Induction motor modelling.
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.
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.
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.
gokul