pmsm drive characteristics and constraint curves -凯发k8网页登录

this example uses motor control blockset™ to show how to use the pmsm characteristic plotting and pmsm milestone speed identification functions to obtain a control trajectory.

required mathworks® products

motor control blockset
simulink®

pmsm drive characteristics and constraint curves

the permanent magnet synchronous motors (pmsms) come in different configurations.

when the permanent magnets (pms) are mounted on the surface of the rotor, the motor is called a surface-mounted pmsm (spmsm).
when the pms are embedded inside the rotor below the surface, the motor is called an interior pmsm (ipmsm).

various parameters affect whether you select an spmsm or ipmsm for a given application.

in the $d$ - $q$ rotor reference frame for pmsms, the $d$ -axis is the axis parallel to the rotor's magnetic orientation and the $q$ -axis is the perpendicular axis, leading the $d$ -axis. an spmsm has equal $d$ -axis and $q$ -axis inductances. an ipmsm has a higher $q$ -axis inductance than its $d$ -axis inductance. this helps the ipmsm to utilize the reluctance torque in addition to the magnetic torque.

the corner speed of a pmsm is the speed at which the torque-vs-speed (drive characteristics) curve of the pmsm changes shape (has a corner) for a given current limit. when you plot the drive characteristics for rated torque (at rated current), the corner speed is also known as the rated speed or base speed.

you draw the constraint curves for a pmsm to understand the possible operational area and boundaries.

maximum torque per ampere (mtpa) trajectory of an ipmsm provides the highest torque for the allowed stator current.
maximum torque per voltage (mtpv) trajectory of an ipmsm provides the highest torque while satisfying the stator voltage constraint.
voltage limit curve (ellipse for an ipmsm or circle for spmsm) is centered around a point that is called the characteristic current. this curve shrinks with increasing rotor speed.

these figures show three segments corresponding to the motor operation. the motor operates in each segment as follows.

segment i

to begin the operation, the motor increases the torque from zero to the rated torque.
the operating point shifts from point o to a along the mtpa trajectory (line segment oa). this corresponds to the $i_{d}$ and $i_{q}$ values at the tip of mtpa.
the motor speed $ω_{1}$ , at which the voltage limit ellipse touches the mtpa curve and the current limit circle, is the rated speed. this is the first milestone speed of the characteristic curve.

segment ii

at the operating point a, the control reaches the maximum possible current, and the field weakening strategy begins.
the current trajectory moves along the current limit circle while maintaining the current constraint from point a to b (line segment ab).
the operating point b is the speed at which the voltage ellipse touches the current limit circle and the mtpv curve. this is the second milestone speed of the characteristic curve.

segment iii

beyond the point b, the voltage ellipse shrinks further and the mtpv trajectory provides the optimal torque until the maximum speed is reached at the operating point c (line segment bc).
the theoretical maximum speed for a friction-free motor with mtpv is infinite. practically, the frictional components in the motor limit the maximum speed to a finite value. in motor control blockset™ software, the maximum speed, including the friction, is when the frictional torque balances the motor torque produced at that speed. the practical speed limit of the motor is usually lower than this maximum speed value. several factors regulate this limit, including mechanical stability, lower efficiency at higher speeds, and so on.
the calculated maximum speed is the third milestone speed. in some motor configurations, mtpv is not possible. in such cases, the frictional torque limits the operating point b, and there is no operating point c. such motors have two milestone speeds.

additionally, note that in these figures:

the drive characteristic curves show the operating points o, a, and b. these curves do not show the operating point c because it is far away from the remaining points.
the operating point a is also at the rated speed $ω_{1}$ . the second speed $ω_{2}$ for the curve is higher than rated speed, but lower than the speed at operating point b. the speed at the operating point b is also the speed at which the current trajectory changes direction from following the current limit circle to following the mtpv trajectory.
the drive characteristics do not show segment i because the motor maintains torque at maximum value, and the $i_{d}$ and $i_{q}$ values remain at values corresponding to the rated torque.

this table shows the functions used in this example.

the pmsm and inverter arguments are structures. this table shows the fields in each structure.

maximum and milestone speeds of pmsm

to get all the milestone speeds of the pmsm, use the mcbpmsmspeeds function. the operable speeds are different for different field weakening control (fwc) strategies. the default fwc method is set to voltage current limited maximum torque (vclmt) control, also called optimal current vector control. the function returns the speed values in rpm.

set the fields for the pmsm and inverter structures and calculate the milestone speeds.


%inverter = mcb_setinverterparameters( 'boostxl-drv8305'); % use an inverter structure from a template
inverter.v_dc=24;
pmsm     = mcb_setpmsmmotorparameters('bly171d');
milestone_speeds=mcbpmsmspeeds(pmsm,inverter);
disp(milestone_speeds);

        5393
        9373

pmsm.i_rated=8;% assuming a higher i_rated to be able to get mtpv trajectory
[milestone_speeds]=mcbpmsmspeeds(pmsm,inverter);
disp(milestone_speeds);

        2342
        2958
       24369

speed milestones of pmsm with different fw control methods

calculate speed milestones for different field weakening methods. when you set the verbose option to 1, each command displays the milestone speeds with the description.

set the fields for the pmsm and inverter structures and calculate the milestone speeds.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
verbose=0; % setting the verbose=1 will also print the messages

milestone_speeds=mcbpmsmspeeds(pmsm,inverter,'verbose',verbose,'fwcmethod','vclmt');
disp(milestone_speeds);

        5393
        9373

milestone_speeds=mcbpmsmspeeds(pmsm,inverter,'verbose',verbose,'fwcmethod','cvcp');
disp(milestone_speeds);

        5393
        7008
        8361

milestone_speeds=mcbpmsmspeeds(pmsm,inverter,'verbose',verbose,'fwcmethod','cccp');
disp(milestone_speeds);

        5393
        7008

milestone_speeds=mcbpmsmspeeds(pmsm,inverter,'verbose',verbose,'fwcmethod','none'); % no id control, id is always set to 0. 
disp(milestone_speeds);

        5393
        6226

$i_{d}$ and $i_{q}$ at different milestone speeds for different fw control methods

calculate $i_{d}$ and $i_{q}$ values at the speed milestones for different field weakening methods. each command outputs the array with columns for each milestone speed, and each row of the column provides $i_{d}$ , $i_{q}$ , and speed (in electrical rad/s).

set the fields for the pmsm and inverter structures and calculate currents and milestone speeds.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');

[milestone_speeds]=mcbpmsmspeeds(pmsm,inverter,'fwcmethod','vclmt','outputall',0);
[outputarray]=mcbpmsmspeeds(pmsm,inverter,'fwcmethod','vclmt','outputall',1);
disp(outputarray);

   1.0e 03 *
         0   -0.0018
    0.0018    0.0004
    2.2591    3.9262

disp(round([milestone_speeds(1) outputarray(3,1)*(60/(2*pi))/pmsm.p;    milestone_speeds(2) outputarray(3,2)*(60/(2*pi))/pmsm.p])); % comparing the corner rpm & max rpm, and the outputarray's corresponding field.

        5393        5393
        9373        9373

outputarray=mcbpmsmspeeds(pmsm,inverter,'fwcmethod','cvcp','outputall',1);
disp(outputarray);

   1.0e 03 *
         0   -0.0011   -0.0018
    0.0018    0.0014    0.0003
    2.2591    2.9355    3.5023

outputarray=mcbpmsmspeeds(pmsm,inverter,'fwcmethod','cccp','outputall',1);
disp(outputarray);

   1.0e 03 *
         0   -0.0011
    0.0018    0.0014
    2.2591    2.9355

outputarray=mcbpmsmspeeds(pmsm,inverter,'fwcmethod','none','outputall',1);
disp(outputarray);

   1.0e 03 *
         0         0
    0.0018    0.0002
    2.2591    2.6081

plot drive characteristics

the drive characteristics display the torque-vs-speed, power-vs-speed, and current-vs-speed characteristic plots of a motor under the given operating constraints. use the mcbpmsmcharacteristics function to plot the drive characteristics.

plot the drive characteristics of the pmsm under different field weakening control methods. you can also plot the constraint curves simultaneously, in which case, you also plot the current trajectory (to get maximum torque for a given speed) for the chosen drive characteristics. when you set the drivecharacteristics option to 2, the function plots an additional figure in the $v_{d}$ - $v_{q}$ space with the voltage constraint curve and the voltage trajectory. the default field weakening control is set to vclmt.

set the fields for the pmsm and inverter structures and plot the characteristics.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
mcbpmsmcharacteristics(pmsm,inverter,'drivecharacteristics',2,'constraintcurves',1,'fwcmethod','vclmt')

mcbpmsmcharacteristics(pmsm,inverter,'drivecharacteristics',1,'constraintcurves',1,'fwcmethod','cvcp')

mcbpmsmcharacteristics(pmsm,inverter,'drivecharacteristics',1,'constraintcurves',1,'fwcmethod','cccp')

mcbpmsmcharacteristics(pmsm,inverter,'drivecharacteristics',1,'constraintcurves',1,'fwcmethod','none')

plot drive characteristics for $i_{\max}$

when the motor operates continuously supporting a constant load torque, that condition is the rated operating condition, and the operating current is the rated current. the short-term rating of the motor is higher than the continuous-current rating.

for this example, pmsm.i_rated is set to 1.8 a and pmsm.i_max is set to 8 a. you can see the comparison between the drive characteristics and use this to determine the short-term operating characteristics.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
mcbpmsmcharacteristics(pmsm,inverter,"fwcmethod","vclmt","drivecharacteristics",1,"imax",8,"constraintcurves",0)

next, assume a high continuous rating to let the motor run with mtpv under rated conditions. the maximum current $i_{\max}$ also operates the motor under the mtpv condition, resulting in no change in the maximum speed.

pmsm.i_rated=6;
mcbpmsmcharacteristics(pmsm,inverter,"fwcmethod","vclmt","drivecharacteristics",1,"imax",10,"constraintcurves",0)

when you specify the field weakening method as none, you might not see a change in maximum speed for the original motor (unchanged rated current) even with higher current limit, but only the peak torque and power changes.

pmsm.i_rated=1.8; % setting the rated current to the normal (low) value)
mcbpmsmcharacteristics(pmsm,inverter,"fwcmethod","cvcp","drivecharacteristics",1,"imax",8,"constraintcurves",0)

mcbpmsmcharacteristics(pmsm,inverter,"fwcmethod","cccp","drivecharacteristics",1,"imax",8,"constraintcurves",0)

mcbpmsmcharacteristics(pmsm,inverter,"fwcmethod","none","drivecharacteristics",1,"imax",8,"constraintcurves",0)

customize drive characteristics

customize the plots using the optional name-value arguments for the characteristics function.

the pmsm structure expects values for the p, rs, ld, lq, fluxpm, b and i_rated fields. the inverter structure expects a value for the v_dc field. set the field values and plot the characteristics.


inverter.v_dc=24;
pmsm.p=4;
pmsm.rs=0.45;
pmsm.ld=1e-3;
pmsm.lq=pmsm.ld*1.4;
pmsm.fluxpm=5.2e-3;
pmsm.b=1.16e-5;
pmsm.i_rated=2;
w_rpm=7000;
t_load=0.01;
fwcmethod="vclmt";
mcbpmsmcharacteristics(pmsm,inverter,'speed',w_rpm,'torque',t_load,'fwcmethod',fwcmethod, 'constraintcurves',0,'drivecharacteristics',1);

plot constraint curves in $i_{d}$ - $i_{q}$ space

you can plot the constraint curves for pmsm in the $i_{d}$ - $i_{q}$ space using the mcbpmsmcharacteristics function.

set the fields for the pmsm and inverter structures and plot the characteristics.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
pmsm.rs=0.1;
pmsm.lq=pmsm.ld*1.8;
w_rpm=5000;
t_load=0;
mcbpmsmcharacteristics(pmsm,inverter,'speed',w_rpm,'torque',t_load);
legend("position", [-0.044093,0.16512,0.60536,0.2])
xlim([-14.7 2.1])
ylim([-12.5 7.6])

plot constraint curves at maximum speed

the constraint curves provide a secondary check for whether the calculated speeds are correct. when you plot the constraint curves at maximum speed, the plot shows three intersecting curves: current limit circle, voltage limit curve, and constant torque curve.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
pmsm.rs=0.1;
pmsm.lq=pmsm.ld*1.8;
milestone_speeds=mcbpmsmspeeds(pmsm,inverter);
mcbpmsmcharacteristics(pmsm,inverter,'speed',milestone_speeds(end),torque=0)
xlim([-9.1 2.2])
ylim([-5.9 6.0])
legend("position",[0.57813,0.64246,0.59643,0.2869])

plot constraint curves at milestone speeds

instead of plotting the constraint curves for all the speed milestones in different figures, you can plot them together in the same characteristic plot. in addition to the intersecting curves at maximum speed, notice that at the corner speed, the mtpa curve, the current limit circle, and the voltage constraint curve intersect.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
pmsm.lq=pmsm.ld*1.8;
pmsm.rs=0.1;
[milestone_speeds]=mcbpmsmspeeds(pmsm,inverter,'constraintcurves',1);
legend("position",[0.61205,0.66151,0.375,0.22976])
ax = gca;
chart = ax.children(3);
datatip(chart,-3.099,-3.697,"location","southwest");
chart = ax.children(9);
datatip(chart,-2.764,3.847,"location","northwest");
xlim([-15.0 1.8])
ylim([-12.2 12.3])

disp(milestone_speeds);

        5645
        9641

plot constraint curves with different speeds in same figure

plot the constraint curves of the motor at different speeds in the same figure.

use the opacityoption to set the opacity of the presently plotted constraint curve.

when you set opacity to 1 (default), the function creates a new figure window for plotting.
when you set opacity to a value less than 1 (but more than 0), the function plots the voltage constraint, mtpa, and mtpv curves with changing opacity.

the opacity option is useful when you want to identify changes and track trends in these curves.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
pmsm.rs=0.1;
pmsm.lq=pmsm.ld*1.8;
mcbpmsmcharacteristics(pmsm,inverter,'speed',4000,'torque',0,'opacity',1)
mcbpmsmcharacteristics(pmsm,inverter,'speed',5000,'opacity',0.8)
mcbpmsmcharacteristics(pmsm,inverter,'speed',6000,'opacity',0.7)
mcbpmsmcharacteristics(pmsm,inverter,'speed',7000,'opacity',0.6)
mcbpmsmcharacteristics(pmsm,inverter,'speed',8000,'opacity',0.5)
mcbpmsmcharacteristics(pmsm,inverter,'speed',9000,'opacity',0.4)
mcbpmsmcharacteristics(pmsm,inverter,'speed',10000,'opacity',0.3)
mcbpmsmcharacteristics(pmsm,inverter,'speed',11000,'opacity',0.2)
mcbpmsmcharacteristics(pmsm,inverter,'speed',12000,'opacity',0.1)
xlim([-7.61 1.23])
ylim([-1.55 9.05])
ax2 = gca;
chart2 = ax2.children(3);
datatip(chart2,-7.541,-0.4029,"location","southwest");
chart2 = ax2.children(51);
datatip(chart2,-3.893,5.438,"location","northwest");
xlim([-10.4 2.9])
ylim([-6.0 9.9])

actual vs approximate voltage equations for maximum speed

usually, the pmsm characteristic equations for $v_{d}$ and $v_{q}$ have an $i \times r$ term, which the computations ignore because this value becomes negligible at higher speeds. in this example, you can include or exclude the $i \times r$ term from the $v_{d}$ and $v_{q}$ equations while calculating the maximum speeds. the actual form of the equation is as follows.

$v_{d} = i_{d} r - ω l_{q} i_{q}$

$v_{q} = i_{q} r ω l_{d} i_{d} ω ψ_{m}$

this figure shows the difference in the drive characteristics with actual and approximate equations.

set the fields for the pmsm and inverter structures and plot the characteristics with both the actual and approximate equations.


inverter = mcb_setinverterparameters( 'boostxl-drv8305');
pmsm     = mcb_setpmsmmotorparameters('bly171d');
max_speed_actual_eqn=mcbpmsmspeeds(pmsm,inverter,'verbose',verbose,'voltageequation','actual');
max_speed_approx_eqn=mcbpmsmspeeds(pmsm,inverter,'verbose',verbose,'voltageequation','approximate');
disp([max_speed_actual_eqn max_speed_approx_eqn])

        5393        5359
        9373        8775

temporary_speed=6000;
mcbpmsmcharacteristics(pmsm,inverter,'speed',temporary_speed,'voltageequation','actual',drivecharacteristics=1)
mcbpmsmcharacteristics(pmsm,inverter,'speed',temporary_speed,'voltageequation','approximate','opacity',0.5,drivecharacteristics=1)

xlim([-9.7 2.1])
ylim([-7.1 7.1])
legend("position", [0.65054,0.72607,0.25,0.19286])
ax3 = gca;
chart3 = ax3.children(13);
datatip(chart3,-7.744,2.653,"location","southwest");
chart3 = ax3.children(5);
datatip(chart3,-8.518,3.158,"location","northwest");

customize constraint curves

change the parameters to plot the constraint curves.


inverter.v_dc=24;
pmsm.p=4;
pmsm.rs=0.15;
pmsm.ld=1e-3;
pmsm.lq=pmsm.ld*2.3;
pmsm.b=1.16e-5;
pmsm.fluxpm=5.2e-3;
w_rpm=3500;
t_load=0.01;
pmsm.i_rated=2.5;
mcbpmsmcharacteristics(pmsm,inverter,'speed',w_rpm,'torque',t_load);
xlim([-10.9 6.0])
ylim([-12.0 7.7])
legend("position",[0.32518,0.11984,0.61786,0.2869])

case studies

this section provides case studies that show how to find a possible control trajectory for an ipmsm and spmsm. the case studies use the following workflow:

plot the drive characteristics (including mtpv) and the constraint curves using the mcbpmsmcharacteristics function. for this plot, use the corner speed.
use the mcbpmsmcharacteristics function with drivecharacteristics set to 0 and opacity set to 0.5 to plot the constraint curves over the figure generated in the previous step. for this plot, use the speed at which mtpv starts.
use the mcbpmsmcharacteristics function with drivecharacteristics set to 0 and opacity set to 0.3 to plot the constraint curves over the figure generated in the previous step. for this plot, use the maximum speed.

doing so plots the three voltage limiting curves overlapping the current limit circle at appropriate points in the final figures for both motors. these figures show the current trajectory for field weakening along the highlighted cyan line. this line follows the mtpv path. when the motor resistance is high enough, the mtpv curve shifts dynamically along with the speed, resulting in the current trajectory not following any single mtpv curve. the resulting cyan curve is the achievable mtpv curve because the example uses the actual equations to dynamically solve for operating points and plot the curves.

to observe the shift in the mtpv curve at different resistance values, use the drop-down menu in the following sections to pick a division factor for the resistance. with a higher resistance, the mtpv curves are farther apart, and with a lower resistance, the mtpv curves are closer to each other.

plot these characteristics for each type of motor.

ipmsm case study

set the field values and plot the characteristics.


inverter.v_dc=24;
pmsm.p=4;
pmsm.rs=0.45;
pmsm.rs=pmsm.rs/1;
pmsm.ld=1e-3;
pmsm.lq=1.25*pmsm.ld;
pmsm.fluxpm=5.2e-3;
pmsm.i_rated=10;
pmsm.b=1.16e-5;
[milestone_speeds]=mcbpmsmspeeds(pmsm,inverter);
mcbpmsmcharacteristics(pmsm,inverter,'speed',milestone_speeds(1),'drivecharacteristics',1,'constraintcurves',1)

mcbpmsmcharacteristics(pmsm,inverter,'speed',milestone_speeds(2),'drivecharacteristics',0,'constraintcurves',1,'opacity',0.5)
mcbpmsmcharacteristics(pmsm,inverter,'speed',milestone_speeds(end),'drivecharacteristics',0,'constraintcurves',1,'opacity',0.3)

spmsm case study

set the field values and plot the characteristics.


inverter.v_dc=24;
pmsm.p=4;
pmsm.rs=0.45;
pmsm.rs=pmsm.rs/1;
pmsm.ld=1e-3;
pmsm.lq=pmsm.ld;
pmsm.fluxpm=5.2e-3;
pmsm.i_rated=10;
pmsm.b=1.16e-5;
[milestone_speeds]=mcbpmsmspeeds(pmsm,inverter);
mcbpmsmcharacteristics(pmsm,inverter,'speed',milestone_speeds(1),'drivecharacteristics',1,'constraintcurves',1)

mcbpmsmcharacteristics(pmsm,inverter,'speed',milestone_speeds(2),'drivecharacteristics',0,'constraintcurves',1,'opacity',0.5)
mcbpmsmcharacteristics(pmsm,inverter,'speed',milestone_speeds(end),'drivecharacteristics',0,'constraintcurves',1,'opacity',0.3)

note: you can use the tab button to get a list of supported optional arguments.

in the command line, in the live-script environment, or in the editor, the function provides suggestions for the name-value arguments.

references

[1] mihailovic, zoran. “modeling and control design of vsi-fed pmsm drive systems with active load.” thesis, virginia tech, 1998. .

[2] li, muyang. "flux-weakening control for permanent-magnet synchronous motors based on z-source inverters." thesis, 2014. http://epublications.marquette.edu/theses_open/284.

pmsm drive characteristics and constraint curves -凯发k8网页登录

required mathworks® products