/// <summary>
/// Call this function on each step analysis to gather Bx,By data
/// </summary>
/// <param name="args"></param>
/// <param name="femm"></param>
public void GatherData(TransientStepArgs args, FEMM femm)
{
//rotate
double xRotorAngle = args.RotorAngle;
if (isRotor)
{
// normalize, make rotorAngle inside (-2alpha,2alpha)
xRotorAngle = analyser.Motor.GetNormalizedRotorAngle(args.RotorAngle);
//convert to radian
xRotorAngle *= Math.PI / 180;
}
// get flux density of all elements in stator for 1 step
for (int i = 0; i < elements.Count; i++)
{
var e = elements[i];
double x = e.center.X;
double y = e.center.Y;
//rotate
if (isRotor)
{
double xx = x * Math.Cos(xRotorAngle) - y * Math.Sin(xRotorAngle);
double yy = x * Math.Sin(xRotorAngle) + y * Math.Cos(xRotorAngle);
x = xx;
y = yy;
}
var pv = femm.mo_getpointvalues(x, y);
if (isRotor)
{
double xx = pv.B1 * Math.Cos(-xRotorAngle) - pv.B2 * Math.Sin(-xRotorAngle);
double yy = pv.B1 * Math.Sin(-xRotorAngle) + pv.B2 * Math.Cos(-xRotorAngle);
pv.B1 = xx;
pv.B2 = yy;
xx = pv.H1 * Math.Cos(-xRotorAngle) - pv.H2 * Math.Sin(-xRotorAngle);
yy = pv.H1 * Math.Sin(-xRotorAngle) + pv.H2 * Math.Cos(-xRotorAngle);
pv.H1 = xx;
pv.H2 = yy;
}
Bx[i][args.step] = pv.B1;
By[i][args.step] = pv.B2;
Hx[i][args.step] = pv.H1;
Hy[i][args.step] = pv.H2;
}
}
protected override void measureInOpenedFem(FEMM femm)
{
// Begin to measure
FEMM.LineIntegralResult lir = new FEMM.LineIntegralResult();
VPMMotor Motor = this.Motor as VPMMotor;
VPMRotor Rotor = Motor.Rotor;
Stator3Phase Stator = Motor.Stator;
GeneralParameters GeneralParams = Motor.GeneralParams;
AirgapNormal Airgap = Motor.Airgap;
PMStaticResults Results = this.Results as PMStaticResults;
double xS = Rotor.Rrotor * Math.Cos(Rotor.alpha * 0.9999);
double yS = Rotor.Rrotor * Math.Sin(Rotor.alpha * 0.9999);
// get phiD
femm.mo_addcontour(xS, yS);
//femm.mo_selectpoint(Rotor.xR, Rotor.yR);
femm.mo_addcontour(xS, -yS);
femm.mo_bendcontour(-360 / (2 * Rotor.p), 1);
lir = femm.mo_lineintegral_full();
Results.phiD = Math.Abs(lir.totalBn);
// get phiM
femm.mo_clearcontour();
femm.mo_selectpoint(Rotor.xD, Rotor.yD);
femm.mo_selectpoint(Rotor.xG, Rotor.yG);
FEMM.LineIntegralResult rr = femm.mo_lineintegral(FEMM.LineIntegralType.Bn);
Results.phiM = Math.Abs(rr.totalBn * 2);
// get phib
femm.mo_clearcontour();
femm.mo_selectpoint(Rotor.xH, Rotor.yH);
femm.mo_selectpoint(Rotor.xH, -Rotor.yH);
rr = femm.mo_lineintegral(FEMM.LineIntegralType.Bn);
Results.phib = Math.Abs(rr.totalBn);
femm.mo_clearcontour();
femm.mo_selectpoint(Rotor.xE, Rotor.yE);
femm.mo_selectpoint(Rotor.xA, Rotor.yA);
rr = femm.mo_lineintegral(FEMM.LineIntegralType.Bn);
Results.phib += Math.Abs(rr.totalBn * 2);
// get phisigmaFe
femm.mo_clearcontour();
femm.mo_addcontour(Rotor.xA, Rotor.yA);
femm.mo_addcontour(xS, yS);
rr = femm.mo_lineintegral(FEMM.LineIntegralType.Bn);
Results.phiFe = Math.Abs(rr.totalBn * 2);
// get phisigmaS
double xZ = (Stator.Rinstator + 5) * Math.Cos(2 * Math.PI / (4 * Rotor.p) - 2 * Math.PI / 180);
double yZ = (Stator.Rinstator + 5) * Math.Sin(2 * Math.PI / (4 * Rotor.p) - 2 * Math.PI / 180);
femm.mo_clearcontour();
femm.mo_selectpoint(Rotor.xS, Rotor.yS);
femm.mo_selectpoint(xZ, yZ);
rr = femm.mo_lineintegral(FEMM.LineIntegralType.Bn);
Results.phisigmaS = Math.Abs(rr.totalBn * 2);
// get FM
femm.mo_clearcontour();
femm.mo_addcontour((Rotor.xI + Rotor.xF) / 2, (Rotor.yI + Rotor.yF) / 2);
femm.mo_addcontour((Rotor.xD + Rotor.xG) / 2, (Rotor.yD + Rotor.yG) / 2);
rr = femm.mo_lineintegral(FEMM.LineIntegralType.Ht);
Results.FM = Math.Abs(rr.totalHt);
// get B_airgap
int n = 128;
Results.Bairgap = new PointD[n * 2];
double RR = (Rotor.Rrotor + Stator.Rinstator) / 2;
for (int i = 0; i < n; i++)
{
double a = -(i - n / 2.0) / n * 2 * Rotor.alpha;
double px = RR * Math.Cos(a);
double py = RR * Math.Sin(a);
FEMM.PointValues pv = femm.mo_getpointvalues(px, py);
Results.Bairgap[i].X = 2 * Rotor.alpha * RR * i / n;
Results.Bairgap[i].Y = pv.B1 * Math.Cos(a) + pv.B2 * Math.Sin(a);
if (double.IsNaN(Results.Bairgap[i].Y))
{
Results.Bairgap[i].Y = 0;
}
if (Results.Bdelta_max < Math.Abs(Results.Bairgap[i].Y))
{
Results.Bdelta_max = Math.Abs(Results.Bairgap[i].Y);
}
}
// make a mirror (odd function)
double dd = 2 * Rotor.alpha * RR;
for (int i = 0; i < n; i++)
{
Results.Bairgap[i + n].X = Results.Bairgap[i].X + dd;
Results.Bairgap[i + n].Y = -Results.Bairgap[i].Y;
}
double wd = Rotor.gammaMedeg / 180 * (Rotor.Rrotor + Airgap.delta / 2) * 2 * Math.PI / (2 * Rotor.p);
Results.Bdelta = Results.phiD / (GeneralParams.MotorLength * wd * 1e-6);
// psiM
Dictionary <String, FEMM.CircuitProperties> cps = Stator.getCircuitsPropertiesInAns(femm);
if (cps.ContainsKey("A") && cps.ContainsKey("B") && cps.ContainsKey("C"))
{
Fdq fdq = ParkTransform.abc_dq(cps["A"].fluxlinkage, cps["B"].fluxlinkage, cps["C"].fluxlinkage, 0);
Results.psiM = fdq.Magnitude;
}
else
{
Results.psiM = double.NaN;
}
femm.mo_close();
}