public override double getTorqueInAns(FEMM femm)
{
femm.mo_clearblock();
femm.mo_groupselectblock(Group_BlockLabel_Magnet_Air);
femm.mo_groupselectblock(Group_BlockLabel_Steel);
double torque = femm.mo_blockintegral(FEMM.BlockIntegralType.Steady_state_weighted_stress_tensor_torque);
if (!Motor.GeneralParams.FullBuildFEMModel)
{
torque *= 2 * p;
}
return(torque);
}
private void measureData(FEMM femm)
{
Stator3Phase stator = Motor.Stator as Stator3Phase;
PM_ACAnalysisResults Results = this.Results as PM_ACAnalysisResults;
//select block
double r = (stator.xD + stator.xE) / 2;
foreach (Stator3Phase.Coil sci in stator.coils)
{
int i = stator.coils.IndexOf(sci);
//angle go clockwise from 3 o'clock (=0 degree in decarter), shift +pi/Q (to match rotor)
int nn = stator.Q / (2 * Motor.Rotor.p);
double aa = -2 * Math.PI * i / stator.Q + (nn % 2 == 0 ? Math.PI / stator.Q : 0);
double x = r * Math.Cos(aa);
double y = r * Math.Sin(aa);
if (aa > -Motor.Rotor.alpha || aa < -2 * Math.PI + Motor.Rotor.alpha)
{
femm.mo_selectblock(x, y);
}
}
double r2 = (stator.xF2 + stator.RGap) / 2;
femm.mo_selectblock(r2, 0);
// measure data
double lossMagnetic = femm.mo_blockintegral(FEMM.BlockIntegralType.Losses_Hysteresis);
double lossResistive = femm.mo_blockintegral(FEMM.BlockIntegralType.Losses_Resistive);
Results.coefficient_lossMagnetic = lossMagnetic / (Current * Current * Freq * Freq);
Results.coefficient_lossResistive = lossResistive / (Current * Current);//=3/2*R
Results.freq = Freq;
Results.current = Current;
}
