Tutorial summary | |
---|---|
Project name | tuto_IPMSM_04 |
Machine name | machine_IPMSM_A |
Source project | tuto_IPMSM_03 |
Description | Variable speed simulation at constant current |
In the previous tutorial Fixed-speed sinusoidal current-driven partial load vibro-acoustic simulation with permeance / mmf model the Load Extrapolation Algorithm was used to quickly obtain variable speed results neglecting saturation change over speed.
It is however possible to run a specific vpeed by speed simulation while changing the saturation level as a function of speed. This can be done in a separate project tuto_IPMSM_04_test by copy/paste (or load/save with the GUI) of tuto_IPMSM_03_test.
The variable speed mode based is activated using
Input.Simu.type_varspeed = 2;
The number of speed steps, the min and max speeds are still defined
Input.Simu.N0_min=500;
Input.Simu.N0_max=7000;
Input.Simu.Nspeeds=50;
The control law to be applied at variable speed must therefore be speficied. The variation of the current magnitude and current angle with speed is obtained by specifying
Input.Simu.type_control = 8;
To follow the constant torque/speed curve specified in the reference article at 60 Nm, the following points are specified in the machine file IPMSM_48s8p:
Input.Electrical.N0spec = [0 500 1000 1500 2000 2504 3304 4104 4904 5704 7000]; % rotation speed [rpm]
Idspec=[50 40 30 15 0 -30 -57 -77 -97 -115 -150];
Iqspec=[110 100 85 75 50 39 30 26.4 23.6 23.8 20];
% RMS phase current [A]
Input.Electrical.I0spec = sqrt(Idspec.^2+Iqspec.^2);
% electrical angle [rad]
Input.Electrical.Phi0spec = mod(atan(Iqspec./Idspec),pi);
In the GUI, all these parameters are in Workflow:
You can also preview the control that you entered:
The obtained sound pressure level can still be obtained with plot_VS_ASPL_overall
One can see that the resonances close to 3000 and 6000 rpm present at no-load are still here, but the magnitude of the maximum levels are different: this is because some destructive or constructive interference occur between armature-induced force harmonics and magnet-induced force harmonics.
The valued of id and iq currents can be checked after the calculation using plot_VS_var(’Id0’) and plot_VS_var(’Iq0’):