model SM_ElectricalExcited "Electrical excited synchronous machine with optional damper cage"
extends Modelica.Magnetic.FundamentalWave.Interfaces.PartialBasicInductionMachine(is(start = zeros(m)), Rs(start = 0.03), Lssigma(start = 0.1 / (2 * pi * fsNominal)), final L0(d = 2 * Lmd / m / effectiveStatorTurns ^ 2, q = 2 * Lmq / m / effectiveStatorTurns ^ 2), redeclare final Modelica.Electrical.Machines.Thermal.SynchronousInductionMachines.ThermalAmbientSMEE thermalAmbient(final useDamperCage = useDamperCage, final Te = TeOperational, final Tr = TrOperational), redeclare final Modelica.Electrical.Machines.Interfaces.InductionMachines.ThermalPortSMEE thermalPort(final useDamperCage = useDamperCage), redeclare final Modelica.Electrical.Machines.Interfaces.InductionMachines.ThermalPortSMEE internalThermalPort(final useDamperCage = useDamperCage), redeclare final Modelica.Electrical.Machines.Interfaces.InductionMachines.PowerBalanceSMEE powerBalance(final lossPowerRotorWinding = damperCageLossPower, final powerExcitation = ve * ie, final lossPowerExcitation = excitation.resistor.LossPower, final lossPowerBrush = brush.lossPower, final lossPowerRotorCore = 0));
parameter Modelica.SIunits.Inductance Lmd(start = 1.5 / (2 * pi * fsNominal)) "Stator main field inductance, d-axis"
annotation (Dialog(
tab = "Nominal resistances and inductances",
groupImage = "modelica://Modelica/Resources/Images/Electrical/Machines/SMEE.png"));
parameter Modelica.SIunits.Inductance Lmq(start = 1.5 / (2 * pi * fsNominal)) "Stator main field inductance, q-axis"
annotation (Dialog(tab = "Nominal resistances and inductances"));
parameter Boolean useDamperCage(start = true) "Enable/disable damper cage"
annotation (Dialog(
tab = "Nominal resistances and inductances",
group = "DamperCage"));
parameter Modelica.SIunits.Inductance Lrsigmad(start = 0.05 / (2 * pi * fsNominal)) "Rotor leakage inductance, d-axis, w.r.t. stator side"
annotation (Dialog(
tab = "Nominal resistances and inductances",
group = "DamperCage",
enable = useDamperCage));
parameter Modelica.SIunits.Inductance Lrsigmaq = Lrsigmad "Rotor leakage inductance, q-axis, w.r.t. stator side"
annotation (Dialog(
tab = "Nominal resistances and inductances",
group = "DamperCage",
enable = useDamperCage));
parameter Modelica.SIunits.Resistance Rrd(start = 0.04) "Rotor resistance, d-axis, w.r.t. stator side"
annotation (Dialog(
tab = "Nominal resistances and inductances",
group = "DamperCage",
enable = useDamperCage));
parameter Modelica.SIunits.Resistance Rrq = Rrd "Rotor resistance , q-axis, w.r.t. stator side"
annotation (Dialog(
tab = "Nominal resistances and inductances",
group = "DamperCage",
enable = useDamperCage));
parameter Modelica.SIunits.Temperature TrRef(start = 293.15) "Reference temperature of damper resistances in d- and q-axis"
annotation (Dialog(
tab = "Nominal resistances and inductances",
group = "DamperCage",
enable = useDamperCage));
parameter Modelica.Electrical.Machines.Thermal.LinearTemperatureCoefficient20 alpha20r(start = 0) "Temperature coefficient of damper resistances in d- and q-axis"
annotation (Dialog(
tab = "Nominal resistances and inductances",
group = "DamperCage",
enable = useDamperCage));
parameter Modelica.SIunits.Temperature TrOperational(start = 293.15) "Operational temperature of (optional) damper cage"
annotation (Dialog(
group = "Operational temperatures",
enable = not useThermalPort and useDamperCage));
parameter Modelica.SIunits.Temperature TeOperational(start = 293.15) "Operational excitation temperature"
annotation (Dialog(
group = "Operational temperatures",
enable = not useThermalPort));
parameter Modelica.SIunits.Voltage VsNominal(start = 100) "Nominal stator voltage"
annotation (Dialog(tab = "Excitation"));
parameter Modelica.SIunits.Current IeOpenCircuit(start = 10) "Open circuit excitation current @ nominal voltage and frequency"
annotation (Dialog(tab = "Excitation"));
parameter Modelica.SIunits.Resistance Re(start = 2.5) "Warm excitation resistance"
annotation (Dialog(tab = "Excitation"));
parameter Modelica.SIunits.Temperature TeRef(start = 293.15) "Reference temperature of excitation resistance"
annotation (Dialog(tab = "Excitation"));
parameter Modelica.Electrical.Machines.Thermal.LinearTemperatureCoefficient20 alpha20e(start = 0) "Temperature coefficient of excitation resistance"
annotation (Dialog(tab = "Excitation"));
parameter Real sigmae(min = 0, max = 1, start = 0.025) "Stray fraction of total excitation inductance"
annotation (Dialog(tab = "Excitation"));
parameter Modelica.Electrical.Machines.Losses.BrushParameters brushParameters "Brush loss parameter record"
annotation (Dialog(tab = "Losses"));
output Modelica.SIunits.Voltage ve = pin_ep.v - pin_en.v "Excitation voltage";
output Modelica.SIunits.Current ie = pin_ep.i "Excitation current";
Modelica.Blocks.Interfaces.RealOutput ir[2](start = zeros(2), each final quantity = "ElectricCurrent", each final unit = "A") if useDamperCage "Damper cage currents"
annotation (
Placement(visible = false),
Dialog(showStartAttribute = true));
Modelica.Blocks.Interfaces.RealOutput irRMS(final quantity = "ElectricCurrent", final unit = "A") if useDamperCage "Damper cage RMS current"
annotation (Placement(visible = false));
FundamentalWave.Components.Short short if not useDamperCage "Magnetic connection in case the damper cage is not present"
annotation (Placement(transformation(
origin = {10, -40},
extent = {
{10, 10},
{-10, -10}},
rotation = 270)));
Components.SaliencyCageWinding rotorCage(final Lsigma(d = Lrsigmad, q = Lrsigmaq), final effectiveTurns = sqrt(1.5) * effectiveStatorTurns, final useHeatPort = true, final TRef = TrRef, final TOperational = TrOperational, final RRef(d = Rrd, q = Rrq), final alpha20 = alpha20r) if useDamperCage "Symmetric rotor cage winding including resistances and stray inductances"
annotation (Placement(transformation(
extent = {
{-10, -10},
{10, 10}},
rotation = 90,
origin = {30, -40})));
Components.SinglePhaseWinding excitation(final orientation = 0, final RRef = Re, final TRef = TeRef, final Lsigma = Lesigma, final effectiveTurns = 0.5 * (effectiveStatorTurns * turnsRatio * m), final useHeatPort = true, final TOperational = TeOperational, final alpha20 = alpha20e) "Excitation winding including resistance and stray inductance"
annotation (Placement(transformation(extent = {
{-30, -50},
{-10, -30}})));
Modelica.Electrical.Analog.Interfaces.PositivePin pin_ep "Positive pin of excitation"
annotation (Placement(transformation(extent = {
{-110, 70},
{-90, 50}})));
Modelica.Electrical.Analog.Interfaces.NegativePin pin_en "Negative pin of excitation"
annotation (Placement(transformation(extent = {
{-90, -50},
{-110, -70}})));
Electrical.Machines.Losses.DCMachines.Brush brush(final brushParameters = brushParameters, final useHeatPort = true) annotation (Placement(transformation(
extent = {
{10, -10},
{-10, 10}},
rotation = 90,
origin = {-80, 40})));
protected
final parameter Real turnsRatio = sqrt(2) * VsNominal / (2 * pi * fsNominal * Lmd * IeOpenCircuit) "Stator current / excitation current";
final parameter Modelica.SIunits.Inductance Lesigma = 0.5 * (Lmd * turnsRatio ^ 2 * 3) * sigmae / (1 - sigmae) "Leakage inductance of the excitation winding";
Modelica.Blocks.Interfaces.RealOutput damperCageLossPower(final quantity = "Power", final unit = "W") "Damper losses";
equation
connect(ir,rotorCage.i);
connect(irRMS,rotorCage.iRMS);
connect(damperCageLossPower,rotorCage.lossPower);
if not useDamperCage then
damperCageLossPower = 0;
end if;
connect(pin_en,excitation.pin_n) annotation (Line(
points = {
{-100, -60},
{-100, -60},
{-100, -50},
{-30, -50}},
color = {0, 0, 255}));
connect(pin_ep,brush.p) annotation (Line(
points = {
{-100, 60},
{-80, 60},
{-80, 50}},
color = {0, 0, 255}));
connect(airGap.port_rn,excitation.port_p) annotation (Line(
points = {
{-10, -10},
{-10, -20},
{-10, -20},
{-10, -30}},
color = {255, 128, 0}));
connect(airGap.port_rp,rotorCage.port_n) annotation (Line(
points = {
{10, -10},
{10, -30},
{30, -30}},
color = {255, 128, 0}));
connect(brush.n,excitation.pin_p) annotation (Line(
points = {
{-80, 30},
{-80, -30},
{-30, -30}},
color = {0, 0, 255}));
connect(brush.heatPort,internalThermalPort.heatPortBrush) annotation (Line(
points = {
{-70, 50},
{-40, 50},
{-40, -90}},
color = {191, 0, 0}));
connect(rotorCage.port_p,excitation.port_n) annotation (Line(
points = {
{30, -50},
{-10, -50}},
color = {255, 128, 0}));
connect(rotorCage.heatPortWinding,internalThermalPort.heatPortRotorWinding) annotation (Line(
points = {
{40, -40},
{40, -80},
{-40, -80},
{-40, -90}},
color = {191, 0, 0}));
connect(short.port_n,airGap.port_rp) annotation (Line(
points = {
{10, -30},
{10, -10}},
color = {255, 128, 0}));
connect(short.port_p,excitation.port_n) annotation (Line(
points = {
{10, -50},
{-10, -50}},
color = {255, 128, 0}));
connect(excitation.heatPortWinding,internalThermalPort.heatPortExcitation) annotation (Line(
points = {
{-20, -50},
{-20, -80},
{-40, -80},
{-40, -90}},
color = {191, 0, 0}));
annotation (
defaultComponentName = "smee",
Icon(graphics = {
Ellipse(
extent = {
{-134, 34},
{-66, -34}},
lineColor = {0, 0, 255}),
Line(
points = {
{-100, 50},
{-100, 20},
{-130, 20},
{-130, -4}},
color = {0, 0, 255}),
Line(
points = {
{-130, -4},
{-129, 1},
{-125, 5},
{-120, 6},
{-115, 5},
{-111, 1},
{-110, -4}},
color = {0, 0, 255}),
Line(
points = {
{-110, -4},
{-109, 1},
{-105, 5},
{-100, 6},
{-95, 5},
{-91, 1},
{-90, -4}},
color = {0, 0, 255}),
Line(
points = {
{-90, -4},
{-89, 1},
{-85, 5},
{-80, 6},
{-75, 5},
{-71, 1},
{-70, -4}},
color = {0, 0, 255}),
Line(
points = {
{-100, -50},
{-100, -20},
{-70, -20},
{-70, -2}},
color = {0, 0, 255})}),
Documentation(info = "<html>\n<p>\nThe symmetry of the stator is assumed. For rotor asymmetries can be taken into account by different resistances and stray inductances in the d- and q-axis. The machine models take the following loss effects into account:\n</p>\n\n<ul>\n<li>heat losses in the temperature dependent stator winding resistances</li>\n<li>heat losses in the temperature dependent excitation winding resistance</li>\n<li>optional, when enabled: heat losses in the temperature dependent damper cage resistances</li>\n<li>brush losses in the excitation circuit</li>\n<li>friction losses</li>\n<li>core losses (only eddy current losses, no hysteresis losses)</li>\n<li>stray load losses</li>\n</ul>\n\n<h4>See also</h4>\n<p>\n<a href=\"modelica://Modelica.Magnetic.FundamentalWave.BasicMachines.SynchronousInductionMachines.SM_PermanentMagnet\">SM_PermanentMagnet</a>,\n<a href=\"modelica://Modelica.Magnetic.FundamentalWave.BasicMachines.SynchronousInductionMachines.SM_ReluctanceRotor\">SM_ReluctanceRotor</a>,\n</p>\n</html>"));
end SM_ElectricalExcited;