WO2024259419A2 - Commande de moteur à l'aide d'une génération de référence d'efficacité optimale - Google Patents

Commande de moteur à l'aide d'une génération de référence d'efficacité optimale Download PDF

Info

Publication number
WO2024259419A2
WO2024259419A2 PCT/US2024/034328 US2024034328W WO2024259419A2 WO 2024259419 A2 WO2024259419 A2 WO 2024259419A2 US 2024034328 W US2024034328 W US 2024034328W WO 2024259419 A2 WO2024259419 A2 WO 2024259419A2
Authority
WO
WIPO (PCT)
Prior art keywords
motor
control parameter
reference frame
function
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/US2024/034328
Other languages
English (en)
Other versions
WO2024259419A3 (fr
Inventor
Matthias PREINDL
Bernard William STEYAERT
Ethan Bagget SWINT
Walter Wesley PENNINGTON III
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Columbia University in the City of New York
Tau Motors Inc
Original Assignee
Columbia University in the City of New York
Tau Motors Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Columbia University in the City of New York, Tau Motors Inc filed Critical Columbia University in the City of New York
Priority to KR1020267001514A priority Critical patent/KR20260046084A/ko
Priority to EP24824347.9A priority patent/EP4728634A2/fr
Publication of WO2024259419A2 publication Critical patent/WO2024259419A2/fr
Publication of WO2024259419A3 publication Critical patent/WO2024259419A3/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/0085—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for high speeds, e.g. above nominal speed
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/022—Synchronous motors

Definitions

  • a synchronous motor is an alternating current (AC) motor having a stator that is driven by AC supply signals (e.g., one signal for each phase of the stator) to cause rotation of a rotor.
  • AC supply signals e.g., one signal for each phase of the stator
  • the AC supply signals in stator windings of the stator generate magnetic fields that interact with a magnetic field or fields of the rotor to cause rotation of the rotor.
  • the rotation of the rotor is generally synchronous with the frequency of the AC supply current.
  • the rotor may be a permanent magnet rotor, a wound field rotor, or a hybrid rotor including both wound fields and permanent magnets. In the case of a permanent magnet rotor, one or more permanent magnets of the rotor generate the magnetic field or fields of the rotor.
  • a motor controller may control an inverter to provide an AC signal to each phase of the motor based on current rotor position and other -1- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 characteristics of the motor.
  • the physics of the magnetic fields of each stator winding interacting with the rotating rotor can lead to complex mathematics problems that are challenging to create and solve to address factors that lead to efficient driving of the motor, and these challenges can be exacerbated in the case of a wound field synchronous (WFS) motor because of the added wound field rotor.
  • WFS wound field synchronous
  • a WFS motor may also be referred to as WFS machine, a wound rotor synchronous machine (WRSM), a wound field synchronous machine (WFSM), a wound rotor synchronous generator (WRSG), a wound field synchronous generator (WFSG), as well as several other names.
  • WFS machine as a power-dense, permanent magnet free, synchronous machine, has gained significant interest in recent years in the field of transportation electrification.
  • a motor controller for a WFS machine may receive a control input, e.g., a reference current or flux, and control the motor in an attempt to achieve an actual motor current or flux that matches the reference current or flux.
  • the control input may be generated by a reference generation map (or reference map).
  • the reference map may itself receive a control input (e.g., a reference torque (T*) or reference motor speed ( ⁇ *), for example, from a user input (e.g., accelerator pedal or other throttle or torque control) or memory. Based on this control input (e.g., T* or ⁇ *), the reference map may generate as output the control input for the motor controller (or an intermediate value that is further translated to the control input).
  • a control input e.g., a reference torque (T*) or reference motor speed ( ⁇ *)
  • T* or ⁇ * reference torque
  • the reference map may generate as output the control input for the motor controller (or an intermediate value that is further translated to the control input).
  • reference generation maps for electric machines may take some combination of torque and/or speed and output a set of currents that attempt to minimize the electrical losses of the machine.
  • the ability to control the machine at high efficiency by using as a reference map either a static map (like a maximum torque per ampere (MTPA) map) or a dynamic optimization problem (like direct torque model predictive control (MPC)) use both loss models of the machine and a mechanism to operate the efficiency map in real time.
  • the dominant losses are copper loss and core loss, which are generally proportional to torque and speed, respectively.
  • the integration of core losses into a reference map is generally neglected in literature because of computational complexity and domination of copper losses over core losses.
  • reference maps may output reference values (e.g., current or flux values) that do not minimize losses, particularly at higher motor speeds where core losses can increase.
  • OERG optimal efficiency reference generation
  • the OERG-based control is based on an optimization problem (or cost function) that uses a convex loss function. Coefficients of the loss function may be determined using finite element analysis (FEA) data, and may be solved over a wide range of inputs (e.g., torques and speeds), showing different output trajectories (e.g., current trajectories).
  • FFA finite element analysis
  • Machine design engineers may design machines to minimize their core loss by analyzing the effects of eddy currents, hysteresis, and armature reaction effects with different geometries, materials, and laminations. For WFS machines, this can be important as its primary application until recently has been larger, megavolt-ampere (MVA)-sized machines for power generation.
  • MVA megavolt-ampere
  • WFS machines have seen a recent increase in popularity in automotive applications, as a WFS machine is a compromise between two popular machines types in the space: a high power density permanent magnet synchronous machine (PMSM) that may be efficient but expensive, and low power density induction machine (IM) that may be cheap but inefficient.
  • PMSM permanent magnet synchronous machine
  • IM low power density induction machine
  • Some WFS machines use hairpin windings to increase slot fill factor; but, this approach increases core losses and introduces additional manufacturing complexity.
  • -3- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [0011] It may be desirable, particularly in automotive applications, for machines to operate efficiently for a wide range of speed and torque.
  • the proposed models have 12% and 53% average error when compared to FEA.
  • the two models use just 15 and 12 floating point operations each, and use 9 or 729 coefficients each.
  • Example use-cases of the two models are maximum efficiency point selection, real-time control, and FEA outlier detection.
  • some embodiments provided herein are directed to OERG-based motor control.
  • OERG-based motor control that use one of the core loss models described herein.
  • OERG-based motor control is also applicable to other motor types, including other permanent magnet motors, brushless motors with permanent magnet rotors, induction motors, universal motors, reluctance motors (synchronous and switched), and the like.
  • an electric machine serving as an electric motor that outputs mechanical power from input electric power may also operate in reverse and serve as an electric generator that outputs electric power from input mechanical power.
  • a motor system includes a power switching network configured to be coupled to a power supply and to a motor; and an electronic controller.
  • the electronic controller is configured to: determine current values for the motor in a rotational reference frame, each current value associated with a dimension of a set of dimensions -4- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 of the rotational reference frame; determine, based on a desired control parameter, a target motor control parameter value for each dimension of the set of dimensions of the rotational reference frame using an optimization cost function that considers motor speed, copper loss, and core loss; and control the power switching network based on the current values and the target motor control parameter values.
  • a method of controlling a motor is provided.
  • the method includes: determining, by an electronic controller, current values for a motor in a rotational reference frame, each current value associated with a dimension of a set of dimensions of the rotational reference frame; determining, by the electronic controller and based on a desired control parameter, a target motor control parameter value for each dimension of the set of dimensions of the rotational reference frame using an optimization cost function that considers motor speed, copper loss, and core loss; and controlling, by the electronic controller, a power switching network based on the current values and the target motor control parameter values.
  • a non-transitory computer-readable medium storing computer- executable instructions, where the instructions are for causing a processor to: determine current values for a motor in a rotational reference frame, each current value associated with a dimension of a set of dimensions of the rotational reference frame; determine, based on a desired control parameter, a target motor control parameter value for each dimension of the set of dimensions of the rotational reference frame using an optimization cost function that considers motor speed, copper loss, and core loss; and control a power switching network coupled to the motor based on the current values and the target motor control parameter values.
  • FIG.1 illustrates a motor system according to some embodiments. -5- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [0020]
  • FIG. 2 illustrates a motor control system implementing optimal efficiency reference generation (OERG) according to some embodiments.
  • OERG optimal efficiency reference generation
  • FIG. 3 illustrates a process for implementing OERG-based motor control according to some embodiments.
  • FIGS.4A and 4B illustrate electrical losses and efficiencies of a wound field synchronous (WFS) motor for raw finite element analysis (FEA) data compared to analytical loss models.
  • FIG.5 illustrates a solution set of current trajectories for an OERG optimization problem.
  • FIG. 6 illustrates a flux map for a WFS motor showing cross coupling modeled by a continuous linear function.
  • FIG.7A illustrates a current-speed domain (left) and a torque speed-domain (right) for a WFS motor, according to some examples.
  • FIG.6 illustrates a current-speed domain (left) and a torque speed-domain (right) for a WFS motor, according to some examples.
  • FIG. 7B illustrates a function relating currents to fluxes for a WFS motor showing saturation and cross saturation, according to some examples.
  • FIG. 7C illustrates power efficient current trajectories from solving an optimization problem with approximation using ⁇ ⁇ ⁇ ⁇ , according to some examples.
  • FIG. 7D illustrates power efficient current trajectories from solving an optimization problem with approximation using ⁇ ⁇ ⁇ ⁇ , according to some examples.
  • FIG. 8B-C illustrate matrix coefficients of matrix G for global core loss model and binned core loss model multiplied by ⁇ 2 .
  • FIG.9A illustrates a trend of a global core loss model against torque, speed, and flux.
  • FIG. 9B illustrates core losses in the flux domain based on FEA (first row), the global core loss model (second row), and the binned core loss model (third row), and error between the FEA and the global core loss model (fourth row) and error between the FEA and the binned core loss model (fifth row).
  • -6- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [0033]
  • FIG.9A illustrates a trend of a global core loss model against torque, speed, and flux.
  • FIG. 9B illustrates core losses in the flux domain based on FEA (first row), the global core loss model (second row), and the binned core loss model (third row), and
  • FIG. 10 illustrates core losses in the torque-speed domain for the global core loss model (top) and the binned core loss model (bottom).
  • FIG.11A illustrates boxplots showing error of the global and binned core loss models.
  • FIG.11B illustrates average core loss error between FEA and global and binned analytics models according to speed.
  • FIG.12A illustrates a plot of an example piecewise affine (PWA) function.
  • FIG.12B illustrates a plot of an example piecewise quadratic (PWQ) function.
  • FIG. 13 illustrates constraints used to construct a piecewise map according to some examples.
  • FIG.14A illustrates FEA datapoints ⁇ sliced on speeds, according to some examples.
  • FIG. 14B illustrates pareto-optimal datapoints ⁇ ⁇ with iso-power curves, according to some examples.
  • FIG.14C illustrates a pareto-optimal surface ⁇ ⁇ ⁇ , according to some examples.
  • FIG.14D illustrates electrical losses from experimental testing with a machine controlled according to an example simplical complex formed using surface reconstruction.
  • FIG. 15A-15D illustrate mesh reduction applied to pareto-optimal surfaces, according to some examples. DESCRIPTION [0044] One or more embodiments are described and illustrated in the following description and accompanying drawings.
  • non-transitory computer-readable medium comprises all computer-readable media but does not consist of a transitory, propagating signal. Accordingly, non-transitory computer-readable medium may include, for example, a hard disk, a CD-ROM, an optical storage device, a magnetic storage device, a ROM (Read Only Memory), a RAM (Random Access Memory), register memory, a processor cache, or any combination thereof.
  • ROM Read Only Memory
  • RAM Random Access Memory
  • register memory a processor cache
  • connection and “coupled” are used broadly and encompass both direct and indirect connecting and coupling, and may refer to physical or electrical connections or couplings.
  • phase "and/or” used with two or more items is intended to cover the items individually and together.
  • a and/or b is intended to cover: a; b; and a and b.
  • Flux linkage may be described as the change in magnetic field that can be detected as a voltage between two ends of a conductive element.
  • the term “flux” is used herein as abbreviated or shorthand notation for "flux linkage” when discussing the relationship between the magnetic field and electrical circuit within an electromagnetic mechanical machine.
  • Inductance is a quantity derived from the relationship between the flux linkage across an electrical element and the current through that electrical element. Being a non- linear relationship, such inductance may be described as the instantaneous change in flux linkage with respect to current (also referred to as “incremental inductance”); as relative to the total flux linkage ( ⁇ ) at some current (i), where ⁇ / i is the "apparent inductance”; or as relative to the total field energy at some current (i), which is determined by ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ (also referred to as "energy-equivalent inductance").
  • Embodiments described herein provided an reference generation (OERG)-based motor control that integrates core losses into the generation of reference values.
  • OERG reference generation
  • the OERG-based motor control is based on an optimization problem (or cost function) that uses a convex loss function.
  • the core loss models may use the least squares method for determining a quadratic core loss function.
  • motor controllers operate using a rotating reference frame to simplify the motor control.
  • motor characteristics in a stationary reference frame may be measured and transformed into a direct-quadrature-Null (DQN) space, or DQN + rotor (R) space or reference frame (also referred to as the DQNR, RDQNull, and RDQ ⁇ reference frame), using a transform based on the Clarke and Park transforms.
  • the motor characteristics e.g., stator currents, rotor currents, and rotor position
  • FIG. 1 illustrates a motor system 100, according to some embodiments.
  • the motor system 100 includes a power supply 105, a motor drive circuit 110, an electric machine 115 (also referred to as an electric motor or motor 115), and a motor controller 120.
  • the power supply 105 provides direct current (DC) power to the motor drive circuit 110.
  • the motor controller 120 is configured to control the motor drive circuit 110 to apply power from the power supply 105 to the motor 115 to drive rotation of the motor 115.
  • the motor controller 120 is configured to control the motor drive circuit 110 to apply electric power from the motor 115 to the power supply 105.
  • the power supply 105 includes a DC power source that provides the DC power to the motor drive circuit 110.
  • the DC power source may be, for example, one or more batteries, photovoltaic cells, or the like.
  • the power supply 105 -9- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 includes an AC/DC rectifier that receives alternative current (AC) power from an AC power source, which may be a utility grid or external generator. In these embodiments, the AC/DC rectifier outputs the DC power to the motor drive circuit 110.
  • the AC power source is part of the power supply 105 (e.g., in the case of an on-site wind turbine or generator).
  • the power supply 105 includes both the DC power source and the AC/DC rectifier, and the DC power from the power supply 105 to the motor drive circuit 110 is provided from one or both sources.
  • the motor controller 120 includes an electronic processor 125 and a memory 130 (collectively, processing circuitry). Generally, the motor controller 120 monitors characteristics of the motor 115 based on signals received from one or more motor sensors and, based on these characteristics, provides control signals to the motor drive circuit 110.
  • the memory 130 includes one or more of a read only memory (ROM), random access memory (RAM), or other non- transitory computer-readable media.
  • the electronic processor 125 is configured to, among other things, receive instructions and data from the memory 130 and execute the instructions to, for example, carry out the functionality of the motor controller 120 described herein.
  • the memory 130 includes control software defining, among other things, control techniques for the motor 115.
  • the electronic processor 125 may be configured to execute the control software to monitor characteristics of the motor 115, receive operational parameters (e.g., motor commands from an input device (not shown)), and to drive the motor drive circuit 110 in accordance with the operational parameters and monitored characteristics.
  • the input device may be or include, for example, an accelerator pedal of an electric vehicle, a trigger, a dial, a keypad, laptop, smartphone, or the like that outputs one or more operational parameters to the motor controller 120 (e.g., encoded in an analog or digital signal).
  • Example operational parameters that may be input and received by the motor controller 120 include torque commands and/or speed commands.
  • the motor controller 120, the electronic processor 125, and the memory 130 are each illustrated as a respective, single unit, in some embodiments, one or more of these components is a distributed component.
  • the electronic processor 125 includes one or more microprocessors and/or hardware circuit elements
  • the memory 130 includes one or more memories
  • the motor controller 120 includes one or more motor controllers (e.g., each with respective processors and memories).
  • the motor 115 includes a stator assembly and a rotor assembly.
  • the motor 115 may be synchronous motor, for example, a wound field synchronous (WFS) motor, a permanent magnet synchronous (PMS) motor, or a hybrid synchronous motor with a rotor having both wound field(s) and permanent magnet(s).
  • WFS wound field synchronous
  • PMS permanent magnet synchronous
  • hybrid synchronous motor with a rotor having both wound field(s) and permanent magnet(s).
  • the stator assembly includes a stator core and a plurality of stator windings on the stator core that are selectively driven with current to induce magnetic fields that rotate the rotor assembly.
  • the stator core may be, for example, a lamination stack formed by a plurality of laminations.
  • the lamination stack may include a generally annular profile with teeth extending radially inward (in the case of an outer stator) or radially outward (in the case of an inner stator).
  • the stator windings may be wrapped around the teeth or may include conductors that otherwise fill the slots between teeth (i.e., the windings may, in some examples, not actually be wound around another object).
  • the rotor assembly includes a rotor core and one or more field windings that are selectively driven with current to induce magnetic fields that interact with the magnetic fields of the stator assembly to rotate the rotor assembly.
  • the rotor core may be, for example, a lamination stack formed by a plurality of laminations.
  • the lamination stack may include a generally annular profile with teeth extending radially inward (in the case of an outer rotor) or radially outward (in the case of an inner rotor).
  • the rotor windings may be wrapped around the teeth or may include conductors that otherwise fill the slots between teeth.
  • the rotor assembly includes a combination of a permanent magnets and field windings. In embodiments in which the motor 115 is a PMS motor, the rotor assembly includes one or more permanent magnets and is without rotor field windings.
  • the motor 115 is described herein primarily as a synchronous motor, in some examples, the motor 115 is of another type, such as an induction motor, a universal motor, a switched reluctance motors, or another type. Although this example of the motor 115 is described as including teeth, in some examples, the motor 115 does not include teeth, for example, when implemented as a slotless motor.
  • the motor 115 (or, electromagnetic mechanical machine), utilizes one or more controllable magnetic fields that are constructed, or energized, in such a manner as to provide a force or torque between two or more components.
  • the force or torque may arise from the interaction of two or more magnetic fields (at least one of which is controllable) such that the relative motion of one component results in a -11- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 lower energy state due to reduced interference between the fields.
  • the force or torque may also arise from a circuit that a given magnetic field, or combination of fields, must take through the materials of two or more components, such that the relative motion of one or more components results in a lower energy state due to lower reluctance of the magnetic circuit, where reluctance is the ratio of magnetomotive force to the magnetic field strength.
  • the motor drive circuit 110 includes a stator drive circuit coupled to one or more stator windings of the motor 115 and a rotor drive circuit coupled to one or more rotor windings of the motor 115.
  • the motor drive circuit 110 includes a stator drive circuit coupled to one or more stator windings of the motor 115, but does not include a rotor drive circuit.
  • the stator drive circuit includes, for example, a plurality of power switching elements connected in a bridge configuration.
  • the power switching elements are semiconductor switching devices such as, for example, a field effect transistor (FET) (e.g., a metal-oxide-semiconductor field effect transistors (MOSFETs)), a bipolar junction transistor (BJT), or insulated gate bipolar transistor (IGBT).
  • FET field effect transistor
  • MOSFETs metal-oxide-semiconductor field effect transistors
  • BJT bipolar junction transistor
  • IGBT insulated gate bipolar transistor
  • the stator drive circuit may include an output terminal for each phase of the stator assembly of the motor 115.
  • the stator drive circuit may include three output terminals, each connected to a terminal of a respective phase of the stator assembly.
  • the stator drive circuit receives DC power from the DC power supply 105 and control signals from the motor controller 120.
  • the control signals which may be pulse-width modulated control signals having respective duty cycles, control the power switching elements to turn on and off in a coordinated manner to drive the stator windings of the motor 115.
  • the motor controller 120 via the control signals, may control the stator drive circuit to generate a sinusoidal drive signal at each output terminal to drive each phase of the stator assembly of the motor 115 with a respective sinusoidal drive signal.
  • the stator drive circuit may also be referred to as a DC-to-AC inverter.
  • Each phase of the stator assembly of the motor 115 may be associated with one or more stator windings.
  • -12- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [0059]
  • the rotor drive circuit when present, includes, for example, a further one or more power switching elements.
  • the power switching elements of the rotor drive circuit may also be connected in a bridge configuration.
  • the rotor drive circuit may include an output terminal pair coupled across each controllable rotor winding of the rotor assembly of the motor 115.
  • the rotor drive circuit receives DC power from the DC power supply 105 and control signals from the motor controller 120.
  • the control signals which may be pulse-width modulated control signals having respective duty cycles, control the power switching elements of the rotor drive circuit to turn on and off in a coordinated manner to drive the rotor windings of the motor 115.
  • the motor controller 120 via the control signals, may control the rotor drive circuit to generate a DC voltage across each rotor winding.
  • the rotor drive circuit includes a single power switching element, single passive element (e.g., a diode), or a plurality of passive elements (e.g., a plurality of diodes) arranged to control the current through the rotor winding(s).
  • the rotor drive circuit provides a power coupling between the power supply 105, which is stationary (i.e., non-rotating), and the one or more windings of the rotor assembly, which rotates.
  • the rotor drive circuit may include a stationary portion and a rotary portion.
  • the rotor drive circuit may include a slip ring and brushes that provides a conductive connection between the stationary portion and the rotary portion.
  • the rotor drive circuit includes another power coupling type.
  • the rotor drive circuit is or includes a DC-to-DC converter that steps down or steps up DC voltage received from the DC power supply 105 to a desired voltage level for the rotor winding(s).
  • Optimal Efficiency Reference Generation (OREG)-based Motor Control [0062]
  • FIG. 2 illustrates a particular example of the motor system 100, identified as motor system 200, implementing such an OERG control scheme, according to some embodiments.
  • the description of components for FIG. 1 above similarly applies to the components in FIG. 2 sharing the same element numbers or names, except as otherwise provided herein.
  • the motor controller 120 is illustrated as a collection of functional blocks with respective inputs and outputs.
  • Each of the functional blocks may be implemented by a dedicated hardware circuit of the electronic processor 125 of the controller 120, by a block of software or instructions stored -13- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 in the memory 130 and executed by the electronic processor 125, or a combination thereof.
  • the motor drive circuit 110 is further illustrated as including a stator drive circuit 205 and a rotor drive circuit 210.
  • the motor drive circuit 110, stator drive circuit 205, and the rotor drive circuit 210 may each be referred to individually or collectively as a power switching network.
  • a power switching network such as these circuits, may be configured to switch voltage, current, and/or power.
  • the motor 115 is illustrated as a WFS motor having a three-phase stator with three phases (A, B, C) and a rotor field winding (R).
  • the motor 115 is a permanent magnet synchronous motor, a hybrid synchronous motor, or another motor type.
  • the motor controller 120 may not sense or control current through the rotor field winding (R), and the control blocks of the motor controller 120 may not receive, process, or generate rotor field components.
  • the MTPA problem becomes more complex with the addition of the strong saturation in magnetic flux of a WFS motor that operates in the linear and non-linear magnetic regimes to prevent high flux error during saturation and cross-saturation, where the problem compounds across varying speeds.
  • Existing online MTPA methods may be computationally expensive on a controller, while offline MTPA methods include adding a cross-coupling torque term and additional variables that decrease the inductance in saturation to approximate saturation effects, which produces a difficult-to-optimize equation and large lookup table.
  • conventional solutions are limited to well-behaved loss mechanisms such as copper losses – whereas non-linear or higher dimensional loss mechanisms such as core losses, windage losses, bearing losses, etc. are not included.
  • the motor system 100 implements an optimal efficiency reference generation (OERG) control scheme.
  • the OERG control scheme can reduce or minimize electrical losses of the motor 115 given a reference torque (e.g., an input torque command indicating a desired output torque of the motor 115) and a motor speed ( ⁇ ) of the motor 115.
  • a reference torque e.g., an input torque command indicating a desired output torque of the motor 115
  • ⁇ motor speed
  • an OERG optimization problem (see, e.g., equation (7) below) is solved, for example, in real time by the motor controller 120 to generate reference values (e.g., current or flux values) for the motor controller 120 that minimize core losses.
  • the OERG control scheme for reference generation may be online (e.g., embedded in a function and solved real time). In other examples, however, the OERG control scheme for reference generation is offline (e.g., encoded in a map, or lookup table that is referenced during operation of the motor in real time).
  • the term “optimal,” as used herein with respect to efficiency reference generation, may refer to a reference value that is calculated or determined according to one of the techniques described herein, which, as also described herein, can be used in a motor control scheme to provide for a more efficient or optimized motor operation.
  • Optimal efficiency reference generation may also be referred to as accurate efficiency reference generation and/or computationally accurate efficiency reference generation.
  • the OERG techniques may also be referred to as use of a computational twin for efficiency reference generation.
  • this description focuses on WFS motors, similar concepts are applicable to other motor types, including PMS motors, hybrid synchronous motors, universal motors, induction motors, and reluctance motors (synchronous and switched). [0067] In FIG.
  • the functional blocks of the motor controller 120 include a Clarke-Park current transform block 212, current-to-flux linkage map 214, a reference generation function block 215 (also referred to as an OERG function block 215), current-to-flux linkage map 220, a difference calculation block 235, a flux controller 240, an inverse Clarke-Park voltage transform block 245, and a pulse width modulation (PWM) generation block 250.
  • one or more of the functional blocks are combined together or distributed into sub-blocks.
  • An example -15- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 of operation of the motor system 100 and the motor controller 120 of FIG. 2 is provided below with respect to FIG.3.
  • FIG. 3 illustrates a process 300 for implementing OERG motor control.
  • the process 300 is described as being carried out by the motor system 200 of FIG. 2. However, in some embodiments, the process 300 may be implemented by another motor system, for example, another example of the motor system 100. Additionally, although the blocks of the process 300 are illustrated in a particular order, in some embodiments, one or more of the blocks may be executed partially or entirely in parallel, may be executed in a different order than illustrated in FIG.3, or may be bypassed.
  • the motor controller 120 determines current values for the motor 115 in a rotational reference frame, such as the RDQN reference frame.
  • Each current value is associated with a dimension (or axis) of a set of dimensions of the RDQN reference frame.
  • the set of dimensions includes the R (or field (f)), D, and Q dimensions (e.g., i f , i d , i q , also referenced as if,dq).
  • the variables F, f, R, and r are used interchangeably to refer to rotor field characteristics.
  • rotor field current may be expressed as ir or as if
  • rotor field flux linkage may be expressed as ⁇ r or as ⁇ f .
  • the motor controller 120 may determine electrical operational characteristics of the motor 115 in a stationary reference frame; determine a rotational position of the motor 115 (e.g., of the rotor of the motor 115); and transform the electrical operational characteristics and the rotational position to the current values for the motor 115 in the rotational reference frame.
  • the controller 120 may receive current measurements from current sensors 255 configured to sense current of each phase of the stator windings (e.g., ia, ib, ic, also collectively referred to as iabc) and current of the rotor winding(s) (e.g., i f , sometimes referred to as i r ) of the motor 115.
  • the controller 120 e.g., at Clarke-Park transform block 212 may also receive rotational position measurements ( ⁇ ) from a rotational position sensor 260 configured to measure the rotational position of the rotor.
  • the controller 120 may determine current and rotational motor position using other techniques. For example, the controller 120 may use a "sensorless" design to determine the rotor position, for -16- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 example, by inferring rotor position by detecting zero-crossings, peaks, and/or valleys of back electromotive force (emf) signals on the stator windings. Further, the controller 120 may calculate the current values from voltage measurements on the stator windings and/or rotor winding(s) provided by voltage sensors.
  • emf back electromotive force
  • the motor controller 120 may perform, via Clarke-Park transform block 212, a Clarke-Park transform on the determined current i abc and rotational position ( ⁇ ) of the motor 115.
  • the motor controller 120 via current-to-flux linkage transform block 214 (also referred to as current-to-flux linkage map), may determine flux linkage values of the motor 115 based on the current values output by the Clarke-Park transform block 212.
  • the transform block 214 may map input current values to corresponding flux linkage values.
  • the current-flux map of transform block 214 may be obtained with finite element analysis (FEA) or experimental measurements.
  • mapping function may include a lookup table (mapping input current to flux linkage) or a function may be fitted to the resulting data points to provide, where the function receives current as input and provides an approximate flux linkage as an output.
  • the function in some examples, may be a piecewise function (e.g., a piecewise affine function).
  • the motor controller 120 determines, based on a desired control parameter, a target motor control parameter value for each dimension of the set of dimensions of the rotational reference frame using an optimization cost function that considers motor speed, copper loss, and core loss (e.g., using the OERG block 215).
  • the motor controller 120 e.g., at OERG block 215) may receive the desired control parameter in the form of an input command or reference value, which may indicate a desired motor torque value (T*) and/or speed value ( ⁇ *).
  • the desired control parameter may be retrieved from a memory (e.g., the memory 130) or received via an input/output device of the motor controller 120 (e.g., from a user operating a keyboard, pushbutton, level, dial, etc.).
  • the motor controller 120 e.g., at OERG block 215) may further receive a motor speed ( ⁇ ) or torque (T) of the motor 115.
  • the motor torque (T) may be sensed or -17- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 inferred (e.g., from a motor current signal).
  • the motor controller 120 e.g., at OERG block 215) may further receive an indication of the DC voltage (VDC) for the drive circuit 110, which may be, for example, retrieved from a memory or sensed by a voltage sensor. [0073]
  • the motor controller 120 may then apply the OERG function of OERG block 215 to the desired control parameter (T* and/or ⁇ *), the motor speed ( ⁇ ) or torque (T) if either is not a desired control parameter, and DC voltage.
  • the OERG block 215 may solve a real time optimization problem (see equation (7)) based on the desired control parameter (T*), motor speed ( ⁇ ), and DC voltage to generate target current values i*r,dq, as described in further detail below.
  • the OERG block 215 may solve a real time optimization problem (see equation (7) or (45)) based on the desired control parameter (T*), motor speed ( ⁇ ), and DC voltage to generate target flux linkage values ⁇ *r,dq.
  • the optimization cost function considers motor speed, copper loss, and core loss because, for example, the optimization cost function includes these elements as parameters.
  • the output current value of the function block 215 is an intermediate target motor control parameter value that is then further translated to a (final) target motor control parameter value by the current-to-flux linkage map 220.
  • the current-to-flux linkage map 220 may be similar in construction and operation as the current-to-flux linkage map 214. In other examples, such as where the controller 240 of FIG.
  • the motor controller 120 controls the power switching network based on the current values and the target motor control parameter values. For example, the motor controller 120 may generate control signals in the stationary reference frame to drive the motor 115 based on a difference between the target motor control parameter value (e.g., flux linkage value output by block 220) and the flux linkage value for each dimension (e.g., output by the block 212).
  • the target motor control parameter value e.g., flux linkage value output by block 220
  • the flux linkage value for each dimension e.g., output by the block 212).
  • the flux controller 240 may be, for example, a proportional integral derivative (PID) controller, a proportional integral (PI) controller, a lookup table, a model-based controller (e.g., implementing model predictive control (MPC), as described in further detail below), or another regulating control device.
  • the motor controller 120 e.g., via the flux controller 240
  • the flux controller 240 may output voltage commands V f , V d , and V q (also referred to collectively as V f,dq ).
  • the flux controller 240 may determine the output voltage commands Vf,dq so that the difference between each flux linkage value and target flux linkage value is minimized.
  • the motor controller 120 may then transform the voltage commands from the rotational reference frame to the stationary reference frame.
  • the motor controller 120 using the inverse Clarke-Park transform block 245, may perform an inverse Clarke-Park transform on the voltage commands V f,dq to generate voltage commands V f , V a , V b , and V c (also referred to collectively as V f,abc ) in the stationary reference frame.
  • the motor controller 120 may then, using PWM generation block 250, generate a pulse width modulated control signal for each dimension of the stationary reference frame to control the power switching network to drive a stator of the motor.
  • the PWM generation block 250 may implement, and the motor controller 120 may access, respective lookup tables for each of the stator phases and the rotor field winding, where the motor controller 120 provides a voltage command to the respective lookup tables of the PWM generation block 250 (e.g., V a to the lookup table for stator phase A, Vb to the lookup table for stator phase B, Vc to the lookup table for stator phase C, and Vf to the lookup table for the rotor field winding).
  • the PWM generation block 250 may return control signal parameters (e.g., a duty cycle for each PWM signal) for each stator phase and the rotor field winding.
  • the motor controller 120 may provide control signals to the drive circuit 110 including stator drive control signals D a , D b , D c (also collectively referred to as D abc ) and rotor drive control signals D f (sometimes referred to as D r ) in accordance -19- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 with the control signal parameters (e.g., at the particular duty cycles indicated by the voltage commands).
  • the motor drive circuit 110 when in a generator operational mode, based on the control signals, the motor drive circuit 110 is controlled to apply electric power from the motor 115 to the power supply 105 (e.g., to charge the power supply) and/or to another electrical load.
  • the flux controller 240 within the controller 120 may implement current-based motor control (i.e., as a current controller), rather than flux linkage-based control as illustrated in FIG.2.
  • the current-to-flux linkage maps 214 and 220 may not be present in the controller 120, and the target current values i* r,dq and measured current values i r,dq may be provided to the difference calculation block 235 of the controller 120.
  • the difference calculation block 235 of the controller 120 may then provide difference values indicating the differences between the target current values i* r,dq and measured current values i r,dq to the current-based controller block that is provided in place of the flux controller block shown in FIG.2.
  • the current-based controller block may be, for example a proportional integral derivative (PID) controller, a PI controller, a lookup table, or another control device.
  • PID proportional integral derivative
  • the current-based controller block (and, thus, the motor controller 120) may then generate a voltage command for each dimension of the set of dimensions of the rotational reference frame based on the received difference values and the rotational position of the motor ( ⁇ ).
  • the current-based controller may output voltage commands Vf, Vd, and Vq (also referred to collectively as Vf,dq or Vr,dq).
  • the voltage commands may then be used to control the motor similar to as described above with respect to FIG.2 (e.g., via inverse Clarke-Park transform block 245, PWM generation block 250, and motor drive circuit 110).
  • Control of a wound rotor synchronous (WRS) motor may be based on a torque function, ⁇ ⁇ ⁇ ⁇ , ⁇ :R ⁇ ⁇ R, which is a function of flux ⁇ and current ⁇ .
  • ⁇ is the stator cross product matrix 0 0 0 0 ⁇ ⁇ ⁇ 0 0 ⁇ 1 ⁇ (2) 0 1 0 and ⁇ is the number of pole pairs of the machine.
  • Flux is limited by the nonlinear function ⁇ ⁇ ⁇ ⁇ and the current range to a set ⁇ .
  • the map ⁇ ⁇ can include saturation in the form of piecewise affine maps and exhibits saturation as well as a strong cross saturation between the rotor and stator d-axis.
  • the machine’s electrical speed is denoted ⁇ ⁇ R and the machine’s DC voltage (VDC) is denoted ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ R ⁇ .
  • the discrete time state equation with flux as the state variable is ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ (3)
  • ⁇ ⁇ ⁇ R is the sampling time
  • ⁇ ⁇ R ⁇ is the identity matrix.
  • the dq voltage is limited by the DC bus voltage (VDC) of the inverter and the modulation strategy to some ⁇ ⁇ ⁇ .
  • VDC DC bus voltage
  • the rotational position of the motor ( ⁇ ) is also a variable in equation (3) and considered in the state equation.
  • the OREG control described herein considers the machine’s electrical losses including copper losses ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ :R ⁇ ⁇ R ⁇ and iron losses ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ :R ⁇ ⁇ R ⁇ . These losses are -21- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 generally a function of current, flux, and speed of the machine.
  • the quadratic term included is typically the most dominant term. These loss models are especially useful for optimization problems because they are convex due to ⁇ ⁇ 0 and ⁇ ⁇ 0.
  • the losses can be added to make a generalized loss function ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ (6) [0085] As noted, the OERG (FIG. 2) generates a target motor control parameter value that targets minimizing losses ⁇ ⁇ ⁇ ⁇ , ⁇ , which may include the sum of winding (copper) losses ⁇ ⁇ ⁇ ⁇ , ⁇ and core losses ⁇ ⁇ ⁇ ⁇ , ⁇ , for a given reference torque T* and motor speed ( ⁇ ).
  • the output of the OERG block 215 may be a reference current i* (as shown in FIG.2) and/or a reference flux ⁇ *.
  • the winding losses may be defined as ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
  • the core losses may be defined as ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , where i is ⁇ is flux, T is torque, R is a matrix defining winding DC and (skin effect and proximity effect) AC winding losses, and the core conductance that approximates (eddy current and hysteresis effect) core losses.
  • the function (optimization problem) solved by the OERG block 215 may be stated as follows: for a torque reference T* and motor speed ( ⁇ ), the OERG reference current i* and reference flux ⁇ * are: ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ I, ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ (7) ⁇ ⁇ ⁇ ⁇ . ⁇ ⁇ ⁇ (8) (9) -22- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 ⁇ ⁇ ⁇ ⁇ .
  • An additional parameter ⁇ can be added to (10) which is minimized, or ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , and both ⁇ ⁇ and ⁇ can be minimized.
  • the current constraint I can be modified to have a strictly positive rotor current, i.e., ⁇ ⁇ ⁇ 0, in this way the solver will avoid symmetric solutions.
  • an initial guess ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ which can be chosen based on predicted efficient points, or based on previous optimization iterations can be loaded into the solver.
  • equation (9) may define a current-flux relationship for the motor, such as defined by, for example, a current-flux map as implemented by map blocks 214 and 220 (FIG. 2).
  • the current-flux map may be obtained with finite element analysis (FEA) or experimental measurements.
  • FEA finite element analysis
  • An equation that may be used for relating current to flux linkage in an WFS motor without saturation has the form ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , where L is the inductance matrix and ⁇ is the flux-offset vector: ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ .
  • the rotor (r) variables are not used, for example, for motors without a rotor or field winding (e.g., permanent magnet synchronous motors). Additionally, in some examples, multiple of these matrices may be stitched together into a piecewise flux map.
  • the rotational position of the motor ( ⁇ ) is also a variable in equation (4), (5), and/or (7), and considered as part of determining the losses and/or speed or torque reference.
  • the motor controller 120 e.g., the processor 125
  • the motor controller 120 is operable to solve the optimization problem (7) for OERG-based control in real time.
  • the motor controller 120 may implement an online, real-time solver.
  • the real-time solver may be a constrained gradient solver, primal dual interior point solver, or a numerical solver, or the like.
  • the motor controller 120 solves the optimization problem (7) in real time by accessing a map or lookup table generated in advance offline and stored in a memory (e.g., the memory 130).
  • the optimization problem (7) is solved for a range of operating points for the motor offline to generate a set of data points, which are then mapped to a piecewise function (e.g., piecewise affine, piecewise quadratic, piecewise cubic) with domains divided by, for example, motor speed ( ⁇ ), to approximate the optimization problem (7).
  • a piecewise function e.g., piecewise affine, piecewise quadratic, piecewise cubic
  • the piecewise function is stored in the motor controller 120 and executed (i.e., solved) in real-time (online) based on input parameters (e.g., torque reference (T*) and motor speed ( ⁇ )). Additional discussion for generating such piecewise functions, including examples using surface reconstruction techniques and/or mesh reduction techniques, is provided below.
  • input parameters e.g., torque reference (T*) and motor speed ( ⁇ )
  • T* torque reference
  • ⁇ motor speed
  • OERG -24- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 control scheme uses constraints, limiting the function to the explicit parameters and limiting the problem to their feasible sets, which implies that voltages, currents, fluxes, and the relationships between them are well behaved, or well-formed.
  • Experimental Results for OERG [0093] Experimental results for the OERG-based reference generation technique described above with respect to an example of the OERG block 215 (FIG. 2) and block 310 (FIG. 3) are provided below. FEA data for a 65kW WRSM with parameters shown in Table 1 are used to calculate the loss coefficients.
  • Table 1 WRSM Parameters Parameter Value Pole pairs ⁇ 2 Nameplate r-axis inductance ⁇ ⁇ 1.956 mH Nameplate d-axis inductance ⁇ ⁇ 2.420 mH N ameplate q-axis inductance ⁇ ⁇ 0.789 mH Base speed 30001/min Max speed 120001/min DC-link voltage 325 V Maximum power 65 kW Maximum torque 220 Nm [0094]
  • ⁇ and ⁇ were computed by using the least squares approach, e.g., for ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
  • similarly for ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
  • the matrices are -25- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 5 .6 0.0 0.0 0.0 0.0 0.0 0.0 ⁇ ⁇ ⁇ 0.0 0.045 0.0 ⁇ , ⁇ ⁇ ⁇ 0.0 0.0033 0.0 0 .0 0.0 0.045, 0.0 0.0 [0095]
  • the optimization problem (7) can be solved using a solver such as, for example, Matlab’s fmincon over the full operating of range of torques ⁇ ⁇ ⁇ that are within bounds of ( ⁇ ⁇ I) and ( ⁇ ⁇ ⁇ ) via (1) and the field weakening is enforced by the equation (8).
  • the static outputs are local or globally optimal operating points of the machine.
  • Optimal Efficiency Reference Generation Function with Multiple Affine Models
  • producing optimal reference currents that minimize copper loss and core loss for a combination of torque and speed is generally a difficult problem to solve analytically given the many non-linearities in a machine, but can be necessary for efficient operation of a machine.
  • traditional methods of mapping a motor, or non-linear power converter system, in a computationally efficient manner is limited – particularly in a highly dimensional space. For example, the magnetic behavior of a wound rotor synchronous (WRS) machine changes between zero torque and rated torque.
  • RFS wound rotor synchronous
  • a machine at zero torque, a machine may have a saliency ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , and at rated torque, the WRS machine may have a saliency ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ .
  • This variation in behavior may result from the phenomena of magnetic saturation in the WRS machine.
  • an affine magnetics model is created at each of these and an optimization problem for reference generation may be solved at each of these two points. The solution sets may ultimately be used to control the WRS machine.
  • an affine magnetics model is created at more than two points, the optimization problem for reference -26- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 generation is solved at each of these points, and resulting solution sets may ultimately be used to control the WRS machine.
  • the dq-axis stator current of the WRS machine (using the power-invariant Clarke-Park transform) may be limited by a stator rated current ⁇ ⁇ , ⁇ , while the rotor axis current is limited by a rated rotor current ⁇ ⁇ , ⁇ . These limits may be set by thermal constraints.
  • the current set I is thus constrained by a cylindrical shape (shown in FIG.7A).
  • inductances may be represented as a matrix ⁇ ⁇ R ⁇ and flux offsets may be expressed as ⁇ ⁇ R as in (9*).
  • V ⁇ ⁇ ariables ⁇ ⁇ and ⁇ ⁇ refer inductances, which are the diagonal terms of these inductance matrices.
  • the torque set ⁇ ⁇ is ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ R
  • the negative with the positive can achieve ⁇ ⁇ ⁇ , ⁇ ).
  • the magnitude of the stator voltage ( ⁇ ⁇ ⁇ ⁇ ) in the ⁇ ⁇ or dq reference frame for synchronous machines is bounded by ⁇ ⁇ ⁇ ⁇ ⁇ /2 or ⁇ ⁇ ⁇ ⁇ ⁇ / ⁇ 3 (if third harmonic injection is used) in a hexagonal shape.
  • the base speed of the machine ⁇ ⁇ is defined as when the product of the stator flux and electrical speed reaches the max stator voltage, or -28- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [00107] Beyond base speed ⁇ ⁇ the stator flux is decreased (flux weakening) below ⁇ ⁇ , ⁇ while at maximum voltage, allowing for higher speeds, at the cost of decreased torque by (1) of the previous section.
  • the maximum speed of the WRS motor, ⁇ ⁇ may be set by mechanical limits.
  • FIG. 7A shows the torque speed-domain for one quadrant, including flux weakening.
  • a torque (except zero torque and maximum torque) produced by (1) has a non-unique set of currents and fluxes that can produce it. Given a reference (or feedback) torque and reference (or feedback) speed (depending on if using a torque or speed controller), producing the set of currents and fluxes that minimize the electrical losses in the machine can provide for optimal efficiency control.
  • the optimization problem (7) of the previous section may be used. However, in some examples, the optimization problem is modified to not include the constraint (8), but may still use the other constraints (9)-(11). In either case, the optimization problem (7) is solved two separate times for each of the two affine current to flux approximations in (7*) and (8*), once where constraint (9) is ⁇ ⁇ ⁇ ⁇ and once where constraint (9) is ⁇ ⁇ ⁇ ⁇ .
  • the torque equation (1) in constraint (10) is fixed to a reference torque, and constraint (11) fixes the speed to a constant.
  • the optimization problem (cost function) (7) is quadratic and the constraints are all affine except the torque constraint (10), which is quadratic and not convex. For this reason, additional considerations can be included when solving to help the numerical solver to reach a feasible solution.
  • An additional parameter ⁇ can be added to (11), which is minimized, or ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , and both ⁇ ⁇ and ⁇ can be minimized.
  • the current constraint I can be modified to have a strictly positive rotor current, i.e., ⁇ ⁇ ⁇ 0, in this way the solver will avoid symmetric guess ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ which can be chosen based on predicted efficient points, or based on previous optimization iterations can be loaded into the solver.
  • the motor controller 120 e.g., the processor 125
  • the motor controller 120 is operable to solve the optimization problem (7) for control in real time, like described in the previous OERG section.
  • the motor controller 120 may determine which of the affine models to utilize (e.g., at each operation point when a reference is to be generated). For example, the motor controller 120 may detect a motor characteristic of the motor during operation (e.g., motor current or motor torque) and then select the affine model to use to generate the reference based on the motor characteristic.
  • the motor characteristic may correlate to magnetic saturation of the motor.
  • the motor controller 120 may solve the optimization problem (7) where the constraint (9) is ⁇ ⁇ ⁇ ⁇ (corresponding to the first affine model).
  • the motor controller 120 may solve the optimization problem (7) where the constraint (9) is ⁇ ⁇ ⁇ ⁇ (corresponding ot the second affine model).
  • the motor controller 120 may implement an online, real-time solver.
  • the real-time solver may be a constrained gradient solver, primal dual interior point solver, or a numerical solver, or the like.
  • the motor controller 120 e.g., via block 215) -30- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 solves the optimization problem (7) with the selected constraint in real time by accessing a map or lookup table corresponding to the optimization problem (7) with the selected constraint generated in advance offline and stored in a memory (e.g., the memory 130).
  • the optimization problem (7) is solved for a range of operating points for the motor offline to generate a set of data points, which are then mapped to a respective piecewise function (e.g., piecewise affine, piecewise quadratic, piecewise cubic) with domains divided by, for example, motor speed ( ⁇ ), to approximate the optimization problem (7).
  • a respective piecewise function e.g., piecewise affine, piecewise quadratic, piecewise cubic
  • ⁇ motor speed
  • each piecewise function corresponds to one of the affine models (e.g., a first piecewise function corresponding to the first affine model and a second piecewise function corresponds to the second affine model).
  • the piecewise functions are stored in the motor controller 120 and, during operation of the motor, one is selected (e.g., based on saturation as indicated by current or torque relative to a threshold.
  • the motor controller 120 may execute (i.e., solve) in real-time (online) the selected piecewise function based on input parameters (e.g., torque reference (T*) and motor speed ( ⁇ )). Additional discussion for generating such piecewise functions, including examples using surface reconstruction techniques and/or mesh reduction techniques, is provided below. [00114] In some examples, the optimization problem (7) is solved in real time for each constraint (e.g., solved with ⁇ ⁇ ⁇ ⁇ and also solved with ⁇ ⁇ ⁇ ⁇ ), and the motor controller 120 selected to affine model to use in reference by selecting the particular solution corresponding to the selected affine model to use for the reference generation.
  • T* torque reference
  • ⁇ motor speed
  • the motor controller 120 may select the solution to use based on the saturation of the motor, which may be indicated by motor current or torque (e.g., being above or below a threshold), as described above.
  • This section describes use of two affine models, each corresponding to a motor saturation level as indicated by motor current or motor torque, that the motor controller 120 selects between to generate a reference.
  • more than two affine models are used, where each affine model corresponds to a magnetic saturation level (e.g., as indicated and defined by a motor current or motor torque range).
  • the motor controller 120 may select a first affine model when motor current (or torque) is between 0 and a first threshold, may select a second affine model when motor current (or torque) is between the first threshold and a second (higher) threshold, and may select a third affine model when the motor current (or torque) is above the second threshold.
  • a -31- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 corresponding threshold may be included such that each affine model corresponds to a magnetic saturation level range (e.g., as defined by a range of current or torque values).
  • the variables ⁇ and ⁇ are computed by using a least squares approach, e.g., where, for ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
  • An FEA dataset sweeps the parameters ⁇ , ⁇ , ⁇ , ⁇ , and outputs ⁇ ⁇ , ⁇ ⁇ , ⁇ ⁇ .
  • the matrices are 0 .004 0.0 0.0 ⁇ ⁇ ⁇ , ⁇ .
  • Table 2 WRS Motor Drive Parameters Parameter Value Turns ratio ⁇ ⁇ / ⁇ ⁇ 39 Pole pairs ⁇ 2 Stator resistance ⁇ ⁇ 11.732 m ⁇ Rotor resistance (stator referred) ⁇ ⁇ 5.461 m ⁇ Shaft inertia 22.76E-3 kg m ⁇ Switching frequency 10 kHz Sampling frequency 20 kHz Nameplate r-axis inductance ⁇ ⁇ 1.956 mH -32- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 N ameplate d-axis inductance ⁇ 2.420 mH Nameplate q-axis inductance ⁇ ⁇ 0.789 mH Base speed 30001/min Max speed 120001/min DC-link voltage 325 V Maximum power 65 kW Maximum torque 220 Nm [00118] Additionally, some operating points of interest for the WRS motor at zero torque and peak torque are provided below in Table 3.
  • the inductance matrices may be 2 .07 2.12 0.0 0.19 0.19 0.0 ⁇ ⁇ ⁇ 2.07 2.42 0.0 ⁇ , ⁇ ⁇ ⁇ ⁇ 0.18 0.28 0.0 ⁇ . 0 .0 0.0 0.79 0.0 0.0 0.31 -33- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [00120]
  • the error is 278 Vs or 94% for d-axis flux and 101 Vs or 80% for q-axis flux.
  • the optimization problem (7) is solved twice using Matlab’s fmincon over the full operating of range of torques ⁇ ⁇ ⁇ that are within bounds of ( ⁇ ⁇ I) via (1) and the field weakening is enforced by (9).
  • the static outputs are local or globally optimal operating points of the machine.
  • the solution set using the zero torque approximation ⁇ ⁇ ⁇ ⁇ function is shown in FIG. 7C. At low speed, the trajectories generally follow a straight path of positive ⁇ ⁇ , ⁇ ⁇ , ⁇ ⁇ ; then, at higher speeds, the machine field weakens and d-axis current decreases.
  • OERG-based control as described herein, generates reference values that provide efficient motor operation, with more accurate estimations of non- linear losses (e.g., core losses) or general machine behavior, while considering saturation, with less data and computations, and that the OERG-based control is able to be implemented by a motor controller (e.g., microcontroller) in real time.
  • a motor controller e.g., microcontroller
  • These core loss estimation techniques may be used as the core loss term ⁇ ⁇ ⁇ ⁇ , ⁇ in the OERG optimization problem (7) to generate the target motor control parameter value.
  • the two models use just 15 and 12 floating point operations each, and use 9 or 729 coefficients each.
  • the analytical models have been validated through FEA simulation for a wound rotor synchronous machine (WRSM).
  • the core loss estimation methods have 12% and 53% average error over all operating points of the machine. Additionally, these methods are extremely light computationally and use very few coefficients, making them well-suited for real-time controllers for various applications, including in the OERG block 215 of the motor controller 120 (FIG. 2).
  • Example use-cases of the two models include maximum efficiency point selection, use in real- time control, and FEA outlier detection.
  • the current in a three-phase WRSM has two parts, the AC stator current ⁇ ⁇ which utilizes the dq-axis from the power-invariant Clarke-Park transform, and DC rotor (sometimes called field) current ⁇ ⁇ . The rotor is aligned to the stator d-axis.
  • the torque per pole pair of the machine is defined ⁇ :R ⁇ ⁇ R ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ where ⁇ is the stator cross product matrix 0 0 0 and ⁇ is the number of pole pairs of [00128]
  • the maximum torque ⁇ ⁇ and speed ⁇ ⁇ are generally limited by mechanical constraints.
  • the rotor and stator can be “flux weakened" in the sense of a PMSM such that electrically there is one maximum torque (at ⁇ ⁇ , ⁇ and ⁇ ⁇ , ⁇ ) and a theoretically unlimited electrical speed.
  • Core loss also referred to as iron loss, or ⁇ ⁇ [W]
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] Core loss
  • ⁇ ⁇ [W] may be modelled using the Steinmetz Equation, which in its simplest form is ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ (19) where ⁇ is a coefficient, ⁇ ⁇ is the switching frequency, and ⁇ ⁇ is the peak value of the magnetic flux density.
  • the machine speed is the rate at which the magnetic flux of the core material changes, so ⁇ replaces ⁇ ⁇ .
  • Flux density ⁇ ⁇ is proportional to the more commonly used machine flux ⁇ , which leads to the equation ⁇ ⁇ ⁇ ⁇ ⁇
  • equation (20) [00130] While this equation may be too general to apply to a real-world system, a choice of exponents ⁇ and ⁇ can be chosen using some estimations to the underlying physics of the machine. There are many variations of equation (20)) used in motor loss modelling.
  • One example, called the Bertoti iron loss formula uses terms representing hysteric loss, lamination thickness, and excess loss with coefficients ⁇ ⁇ , ⁇ of ⁇ 2,1 ⁇ , ⁇ 2,2 ⁇ , ⁇ 1.5,1.5 ⁇ as ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ . ⁇ ⁇ . ⁇ . ⁇ . ⁇ .
  • the first core loss model is ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , (22) and, when considering ⁇ is a ⁇ 3 ⁇ 1 ⁇ matrix from (16), becomes ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ (23) Where the coefficient ⁇ is distributed ⁇ term of flux ⁇ is simple to compute. A linear term ⁇ ⁇ 1 could potentially be added.
  • This model can be called the global model, or ⁇ ⁇ , ⁇ . Torque and speed both contribute to this ⁇ ⁇ equation, with speed proportional to the ⁇ ⁇ term and torque as part of the flux term ⁇ in equation (17).
  • FIG.9A illustrates a trend of the global model against torque, speed, and flux.
  • the second core loss model is ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ for a set of discrete there is a separate ⁇ ⁇ matrix per discrete speed.
  • This core loss model is binned by speed, and thus denoted ⁇ ⁇ , ⁇ .
  • the second core loss model may also be represented as: ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ cover the -37- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 machine speeds ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ... ⁇ ⁇ ⁇ .
  • the model has n piecewise quadratic equations which each have three matrices of coefficients corresponding to the quadratic dependence on speed ( ⁇ ⁇ ), linear dependence on speed ( ⁇ ⁇ ), and no dependence on speed ( ⁇ ⁇ ), and can be formulated to be continuous.
  • This core loss model is also binned by speed.
  • the matrix ⁇ may be obtained by first having a set of available loss datapoints ⁇ ⁇ ⁇ , ⁇ , ⁇ and solving the following convex optimization problem using all points ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
  • the two core loss models were run on a 65 kW WRSM with parameters shown in TABLE 4, where any such parameter could also be considered a feasible dimension in certain embodiments.
  • the FEA dataset used had 498,606 FEA datapoints ⁇ ⁇ ⁇ , ⁇ , ⁇ corresponding to the full current range of the machine ⁇ ⁇ ⁇ I and 81 specific speeds ⁇ ⁇ .
  • the values for ⁇ were negligible for all terms except ⁇ ⁇ , ⁇ and ⁇ ⁇ , ⁇ which are the self-induced stator core losses of the d-axis and q-axis respectively. This is because for this specific WRSM, the rotor is excited by DC current.
  • the values are shown in TABLE 5 and FIGS.8B-C.
  • core loss error for the global model is shown in the upper plots of FIG.10, and core loss error for the binned model is shown in the lower plots of FIG.10.
  • Boxplots showing the error of both core loss models are shown in FIG. 11A, and average core loss error between FEA and the global and binned analytics models are shown according to speed in FIG. 11B.
  • the error tends to decrease dramatically for both models as speed is increased.
  • the average error for ⁇ ⁇ , ⁇ is 12% compared to an average error of 53% for ⁇ ⁇ , ⁇ .
  • ⁇ ⁇ , ⁇ is much more accurate, but considerably slower than ⁇ ⁇ , ⁇ . It is likely that not all 81 speeds in the piecewise function are necessary to have a reasonably accurate core loss model. Accordingly, in some examples, the piecewise function has fewer than 81 speeds (i.e., fewer bins) and, thus, fewer than 81x more coefficients.
  • Example benefits of these two core loss models versus more complex models are 1) increased speed of computation 2) relatively low error 3) no need to know complex machine geometry, and 4) includes core losses from coupled flux. These benefits allow for a wide range of potential applications including fast maximum efficiency point selection, use a cost in a real- time controller when moving between reference speed-torques, and FEA outlier detection.
  • the state space model of the system uses the flux ⁇ ⁇ ⁇ R ⁇ in the stator dq-axis (using the magnitude-invariant Clarke-Park transform) and rotor axis (aligned to the d-axis) as the state variable ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , (30) ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , (31) ⁇ ⁇ ⁇ ⁇ [00147]
  • the inputs are a the stator and rotor voltages ⁇ ⁇ .
  • ⁇ ⁇ R is the mechanical speed (1/min) multiplied by ⁇ ⁇ ⁇ ⁇ , where ⁇ is the number of pole pairs of the machine.
  • -42- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [00148]
  • the relationship between the current ⁇ ⁇ ⁇ R ⁇ and flux ⁇ ⁇ is nonlinear, and has saturation and cross saturation effects. This can be modelled by the continuous nonlinear function ⁇ ⁇ , , and the inverse can be modelled by ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ .
  • the maximum value of the torque function is dependent on speed by m ax ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ / ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ . (43) Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 [00154]
  • the domain (inputs) of optimal generation map are speeds and torques bounded by (43), the range (output) is a set of currents ⁇ ⁇ that may be requested (e.g., of the controller 240).
  • the largest sources of electrical losses in a machine are copper loss ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ R ⁇ ⁇ R ⁇ and iron loss ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ R ⁇ ⁇ R ⁇ .
  • Copper loss has a typically square dependence on current and linear by resistance.
  • the resistance of the stator and rotor vary non- linearly with machine temperature and speed, and are also frequency dependent.
  • the core loss is even more difficult to model with an analytical function, it typically has a square dependence on flux ⁇ (which is nonlinearly dependent on current) and speed ⁇ .
  • the two largest sources of core loss are eddy currents and hysteresis, which are both difficult to model.
  • problem (45) can be solved for all combinations of (43).
  • the solution set for all torque and speed combinations can be assigned to a continuous function ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ R ⁇ ⁇ R ⁇ . Because there is no equation for the objective function, there is also no analytical solution.
  • a pareto frontier can be constructed using the points in ⁇ ⁇ ⁇ which minimizes electrical loss ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ .
  • Pareto frontiers are a collection of pareto optimal points from a set of discrete datapoints that minimize one dimension of the objective function.
  • the pareto optimal FEA datapoints are denoted ⁇ ⁇ , and ⁇ ⁇ ⁇ ⁇ .
  • ⁇ ⁇ is the discretized FEA solution to (45).
  • the discrete pareto-optimal points ⁇ ⁇ may be used to create a continuous pareto- optimal surface ⁇ ⁇ ⁇ (simplical mesh) that best approximates the ideal surface ⁇ ⁇ .
  • a surface reconstruction technique may be used to translate the discrete pareto-optimal points ⁇ ⁇ to a continuous pareto-optimal surface ⁇ ⁇ ⁇ .
  • Such a surface reconstruction technique may depend on the original surface, ⁇ ⁇ , being a surface (two-dimensional manifold) that is compact, connected, and orientable.
  • a surface reconstruction technique such as described in “Surface Reconstruction from Unorganized Points” (Hoppe et al. 1992), takes as input an organized set of points on or near an unknown manifold M and produces as output a simplicial surface that approximates M.
  • the surface reconstruction technique that is employed includes a first stage to define a function f that estimates the signed geometric distance to the unknown surface M, and a -45- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 second stage that uses a contouring algorithm to approximate Z(f) by a simplical surface, where zero set Z(f) is an estimate for M.
  • an oriented plane may be associated with each of the data points.
  • Each plane (referred to as tangent planes) may serve as a local linear approximation to the surface.
  • the tangent planes may not directly define the surface because their union may have a complicated non-manifold structure.
  • the tangent planes may define the signed distance function to the surface.
  • other surface reconstruction techniques may be employed.
  • the FEA method will produce a set of discrete points ⁇ ⁇ ⁇ ⁇ which will have some sampling density ⁇ and noise factor ⁇ .
  • a sampled space is said to be ⁇ ⁇ dense if for any sphere with radius ⁇ ⁇ R ⁇ there is at least one sample point ⁇ . If the original sample ⁇ has some ⁇ ⁇ density, then the pareto points ⁇ ⁇ will have a ⁇ ⁇ density less than or equal to the original ⁇ ⁇ density, as ⁇ ⁇ ⁇ .
  • ⁇ ⁇ ⁇ density in ⁇ it may be necessary to have a higher ⁇ ⁇ density than in ⁇ , requiring denser sampling in the FEA.
  • Any pareto point ⁇ ⁇ , ⁇ will be equal to the ideal surface with some added error, or ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ .
  • a sampled space is called ⁇ ⁇ noisy if
  • the ⁇ value (or maximum error) for an FEA simulation is generally known, and decreases with the size of the simplical mesh.
  • the discrete points produced by the FEA method will serve as an input to the surface reconstruction technique to generate the be used to translate the discrete pareto-optimal points ⁇ ⁇ to a continuous pareto-optimal surface ⁇ ⁇ ⁇ , which is a simplical mesh.
  • an optimization problem is employed (e.g., in OERG block 215 of FIG.2 and/or in block 310 of FIG.3) that does not consider core losses and, rather, focuses on minimizing copper losses. Such an optimization problem may provide a less complex function while still providing efficient motor operation.
  • the process 300 of FIG.3 and operation of the system 200 may otherwise proceed similarly to the other examples discussed herein.
  • the following optimization problem may minimize copper losses given any torque and any speed that satisfies equation (43): ⁇ ⁇ ⁇ ⁇ min ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ (51) ⁇ ⁇ ⁇ ⁇ .
  • FEA method may be employed to approximate loss using this copper loss-focused optimization problem.
  • the FEA method locally linearizes nonlinear thermal and magnetic equations using simplical meshes given a set of inputs.
  • An example may be inputs of fixed current ⁇ ⁇ I and speed ⁇ ⁇ ⁇ , and outputs may be copper loss ⁇ ⁇ , resistance R, and flux ⁇ .
  • Pareto frontiers are a collection of pareto optimal points from a set of discrete datapoints that minimize one dimension of the objective function[42].
  • the pareto optimal FEA datapoints are denoted ⁇ ⁇ , and ⁇ ⁇ ⁇ ⁇ .
  • ⁇ ⁇ is the discretized FEA solution to (51).
  • a simplical complex in R ⁇ can be constructed using the points ⁇ ⁇ and a triangulation algorithm such as the Delaunay triangulation.
  • the simplical complex is denoted ⁇ ⁇ ⁇ . [00171] Regardless of whether using the optimization problem of (7), (45) or (51), the FEA method may be applied to provide discrete pareto-optimal points ⁇ ⁇ that may be translated into a simplical complex ⁇ ⁇ ⁇ .
  • a surface reconstruction technique may be employed to generate the simplical complex denoted ⁇ ⁇ ⁇
  • Delaunay triangulation may be employed to generate the simplical complex denoted ⁇ ⁇ ⁇
  • another decomposition technique may be employed to generate the simplical complex denoted ⁇ ⁇ ⁇ .
  • Delaunay triangulation is a known mathematical meshing algorithm or technique, and quad tree, box tree, KD-tree, and alpha shape are also known domain decomposition algorithms or techniques.
  • a Voronoi diagram of the current points may be constructed.
  • the Voronoi -47- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 diagram splits the current space into N Voronoi cells, where all points in a Voronoi cell are closer to a single point of the original set of current points than any other point.
  • the Delaunay triangulation one may find the dual of the Voronoi diagram.
  • the Delaunay triangulation maximizes the minimum angle in the simplices it creates, thus reducing “skinny” simplices.
  • “skinny” simplices may be described by their aspect ratio. In other words, a simplex having an aspect ratio above a threshold amount may be considered “skinny,” while a simplex having an aspect ratio below the threshold amount may be considered “not skinny.” [00173] For the resulting current simplices generated, each simplex is connected at a shared boundary to another simplex so that all simplices are connected within the domain, and no simplex overlaps another simplex within the domain. Additionally, the simplices may be defined such that they are closed domains on one side and open domains on the other such that any arbitrary point in the domain will belong to one and only one simplex within the Delaunay construction (even if the point is on a boundary).
  • the resulting simplical complex ⁇ ⁇ ⁇ may be a collection or mesh of simplices (e.g., of two-dimensional simplices in three-dimensional space).
  • a piecewise map or function may be fitted to the resulting simplical complex ⁇ ⁇ ⁇ , where the function may then be used to approximate the solution set.
  • the surface reconstructed as described above may be a collection of simplices in R ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ 1/ ⁇ ⁇ ⁇ space.
  • the optimal output current set ⁇ ⁇ is generated by using a piecewise affine function defined by the vertices of ⁇ ⁇ ⁇ .
  • Each simplex will have an affine equation assigned to it such that the overall function will be closed and continuous.
  • Each simplex (plane) in ⁇ ⁇ ⁇ is defined as the convex hull of three of three points ⁇ , ⁇ ⁇ H ⁇ ⁇ ⁇ ⁇ ⁇ , ⁇ , ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ .
  • FIG. 14B shows pareto-optimal datapoints ⁇ ⁇ with iso- power curves
  • FIG.14C shows pareto-optimal surface ⁇ ⁇ ⁇ .
  • FIG.14D illustrates electrical losses from experimental testing with a machine controlled according an example simplical complex formed using surface reconstruction. In this testing, the machine is controlled with a 120 second drive cycle with positive and negative torques, where the machine includes -49- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 parameters as shown in Table 4 above. The current from h ⁇ ⁇ is mapped back to torque to show the difference in torque to the original requested torque.
  • Piecewise Maps [00181] As noted above, in some examples, the OERG-based motor control described herein uses a piecewise map (also referred to as a piecewise function).
  • a piecewise map may be a function that is fitted to a solution set of data points, where the function may then be used to approximate the solution set.
  • Piecewise maps divide a nonlinear map into M domains, where each of the M domains is made up of a (sub) function. In other words, sub-functions are pieces of the piecewise map and, collectively, form the piecewise map.
  • Piecewise maps may be classified as a piecewise constant map, piecewise affine map, piecewise quadratic map, piecewise cubic map, or a piecewise map with functions having an order higher than three.
  • the nonlinear map is divided into M domains over which the function may be constant values.
  • the nonlinear map is divided into M domains over which the function may be linearized.
  • the nonlinear map is divided into M domains over which the function may be quadratic.
  • the nonlinear map is divided into M domains over which the function may be cubic.
  • Piecewise maps of a higher order are similarly divided into M domains over which the function may be of the higher order. [00182] Piecewise maps divide the original domain into M domains or sets.
  • FIGS.12A-12B illustrate two piecewise maps. More particularly, FIG.12A illustrates a piecewise affine map (PWA map) and FIG. 12B illustrates a piecewise quadratic map (PWQ map).
  • Piecewise maps may be constructed with regularly or irregularly sampled points for data in any dimension, although there are some limitations for higher order polynomials.
  • to generate a piecewise map may include letting ⁇ ⁇ R ⁇ and ⁇ :R ⁇ ⁇ R be a potentially unknown C ⁇ function that is irregularly sampled m times subject to ⁇ ⁇ ⁇ -50- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 ... ⁇ ⁇ ⁇ ⁇ , with ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ... ⁇ and ⁇ ⁇ ⁇ ⁇ ⁇ ... ⁇ . Attention may be restricted to the convex hull ⁇ ⁇ ⁇ ⁇ hull ⁇ ⁇ ⁇ ⁇ R ⁇ .
  • the samples ⁇ may be triangulated into l simplices, for example, using the n-dimensional Delaunay method.
  • ⁇ ⁇ :l ⁇ ⁇ ⁇ ⁇ is connected, non-overlapping, and convex.
  • Each ⁇ is defined by ⁇ ⁇ 1 samples and ⁇ ⁇ hull ⁇ ⁇ ) with ⁇ ⁇ ⁇ ⁇ ⁇ , ... ⁇ ⁇ ⁇ .
  • piecewise multivariate polynomials ⁇ : ⁇ ⁇ R may be defined to approximate ⁇ .
  • piecewise multivariate polynomisals may be defined to approximate ⁇ as: ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ [00185]
  • Example polynomials are form, in Table 7. Table 6 – Example Polynomials ⁇ of parameter Order Type Function Gradient Hessian ⁇ ⁇ 0 constant ⁇ ⁇ 0 0 1 1 affine ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ 0 ⁇ ⁇ 1 2 quadratic ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ 1 cubic 3 ...
  • l ⁇ ⁇ ⁇ ⁇ will not leave any degrees of freedom to a fitting function. Therefore, it may be required that l ⁇ ⁇ . This can be useful for fitting to a low number of points and/or simplices.
  • An example relevant special case includes piecewise cubic functions (PWC).
  • piecewise fitting may include, where ⁇ are coordinates in ⁇ - dimensional space, ⁇ are the number of coordinates and values ⁇ ⁇ ⁇ , ⁇ ⁇ ⁇ , and l are the number of (Delaunay) simplices that triangulate the space, the following: ⁇ ⁇ vertices per simplex: ⁇ ⁇ ⁇ ⁇ ⁇ 1 ⁇ ⁇ parameters ⁇ per simplex of ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ "Full ⁇ ⁇ ": ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ 1 ⁇ ⁇ ⁇ ⁇ 1 ⁇ ⁇ 1 "sym ⁇ ⁇ ": ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ 1 ⁇ ⁇ ⁇ values: ⁇ ⁇ ⁇ ⁇ ⁇ ⁇
  • Function value constraint ⁇ ⁇ ⁇ ⁇ ⁇
  • a data set of input operational points and corresponding output operational points is generated.
  • the data set may be generated through simulation (e.g., using finite element analysis (FEA)), through experimentation, or through a combination of simulation and experimentation.
  • FEA finite element analysis
  • a domain decomposition algorithm is applied to the data set to generate simplices.
  • Various domain decomposition algorithms or techniques also referred to as domain subdivision algorithms, may be applied to generate the simplices.
  • the domain decomposition algorithm or technique may be Delaunay triangulation or may be a surface reconstruction technique as described above.
  • the domain decomposition algorithm or technique may be an irregularly sampled, but rectangular, decomposition, for example, a quad tree algorithm, a box tree algorithm (also referred to as oct tree), or KD-tree algorithm (depending upon the number of independent dimensions).
  • the domain decomposition algorithm is an alpha shape algorithm or technique. For the resulting simplices generated, each simplex is connected at a shared boundary to another simplex so that all simplices are connected within the domain, and no simplex overlaps another simplex within the domain.
  • the simplices may be defined such that they are closed domains on one side and open domains on the other such that any arbitrary point in the domain will belong to one and only one simplex within the Delaunay construction (even if it is on a boundary). Examples of such connected simplices are shown in the PWA map of FIG.12A and the PWQ map of FIG.12B.
  • the binned core loss model (24) is implemented using a piecewise map.
  • the piecewise map is a piecewise quadratic (PWQ) map.
  • each of the M domains of the PWQ map corresponds to a motor speed range (e.g., ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ).
  • the motor controller 120 may select the M domain of the PWQ map (and, thus, the sub-function of the PWQ map) based on the motor speed of the motor 115. For example, when motor speed ( ⁇ ) is between ⁇ ⁇ and ⁇ ⁇ , the controller 120 will select and solve ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ to determine the core loss. This core loss may then be used in the reference generation by the OERG block 215 (e.g., when solving the optimization problem (7)).
  • the domains may be according to reference torque, motor speed, and additional parameters, such as, for example, current and/or voltage limit.
  • This technique in some examples, can be viewed as taking a system that has more than two degrees of freedom beyond torque and speed, and running an optimization to describe how best to collapse those additional degrees of freedom into torque and speed, so that the system operates withing its constraints.
  • the current-flux map (equation (9)) is implemented using a piecewise map, which may be a PWA map, PWC map, PWQ map, etc.
  • Using decomposition techniques, like Delaunay triangulation, to create a mesh over a multidimensional space e.g., dq0 or rdq0, for instance
  • a multidimensional space e.g., dq0 or rdq0, for instance
  • dq0 or rdq0 and speed, or rdq0, speed, and core losses, etc. the meshing becomes more challenging.
  • the mesh when relying strictly on the data to create the mesh, the mesh can become ill-formed or "not smooth.” Strictly relying on the data may refer to using FEA data to construct a motor model (e.g., describing the machine, which is used for the purposes of control). When operating a machine across a trajectory that crosses an ill-formed mesh, a proper response is not formulated. As a result, the ill-formed mesh affects control and motor dynamics because current may reverse, may jump from "peaks" to "valleys,” or the like, across this ill-formed (noisy) mesh.
  • constraints on the system that generate the meshes e.g., optimization problem (7)
  • constraints on the system that generate the meshes can create a much smoother surface, particularly when compared to a system relying only on the data (without such constraints) and/or without considering core losses to construct the mesh.
  • a trajectory -58- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 across the mesh (e.g., with piecewise functions) provides a more effective traversal.
  • a trajectory for example, identifies the most advantageous set point or operating points for the motor to produce certain torque and speed across a drive cycle or operating condition.
  • the motor controller 120 can consider multiple speeds and torques, and create a trajectory across flux and speed that operates the machine to accomplish that output at a high efficiency.
  • the PWA function may be generated using a separate computing device.
  • a computing device e.g., server, desktop, laptop, etc. having a memory and a processor, where the processor executes instructions retrieved from the memory to perform the various processing steps, algorithms, and techniques described above (e.g., FEA analysis, surface reconstruction, mesh reduction (described below), and the like) to generate the piecewise function.
  • the piecewise function may then be transmitted by the computing device (or another intermediary device) to the motor controller 120 for storage on the memory 130.
  • Mesh Reduction [00209]
  • a mesh reduction algorithm may be applied to the simplical complex.
  • each simplex corresponds to a domain or function of the piecewise function
  • the complexity and size of the ultimate piecewise function may be reduced.
  • less memory space may be used to store the piecewise function and a controller may execute the piecewise function (e.g., determine a target motor control parameter value based on a desired control parameter) more quickly.
  • the mesh reduction algorithm simplification is configured to maintain sufficient accuracy in its approximation of the optimization problem to remain effective and provide efficient reference generation.
  • the mesh reduction algorithm may be, for example, an edge contraction algorithm, a vertex contraction -59- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 algorithm, and/or a vertex decimation algorithm.
  • a vertex contraction algorithm may be based on an iterative contraction of vertex pairs where, to contract a vertex pair, the vertices of the pair are moved to a new position, the pair’s incident edges are connected to one vertex of the pair, the other vertex of the pair is deleted, and, subsequently, edges or faces that have become degenerate are removed.
  • the simplical complex (also called simplical mesh) ⁇ ⁇ ⁇ has a connected domain; that is, in the domain ( ⁇ ⁇ , ⁇ ), there are no gaps in the triangles.
  • a connected domain is generally guaranteed for the Delaunay triangulation method, but, not for all surface reconstruction methods.
  • the domain ( ⁇ ⁇ , ⁇ ) is bounded by (43), so the surface is open.
  • ⁇ ⁇ ⁇ is a connected, open two-dinemsional manifold.
  • Mesh reduction algorithms aim to reduce the number of simplices in a simplical complex while preserving the general shape. They are either topology preserving or non- topology preserving. Non-topology preserving algorithms may change the topological properties of the surface.
  • non-topology preserving algorithms include vertex contraction and vertex clustering.
  • vertex contraction and vertex clustering For the PWA map, it may be detrimental for the map to go from connected to unconnected, as the output reference currents would be undefined.
  • Two iterative approaches that mesh reduction algorithms may follow include: 1) setting the maximum number of simplices or 2) setting the maximum allowable error. Some algorithms work with both such as edge contraction, vertex contraction, and vertex decimation. The memory and time constraints may be directly dependent on the number of simplices ⁇ . Accordingly, in some examlpes, the first option is used to specify the maximum number of simplices. Some algorithms, for example, simplification envelopes, set a maximum Euclidean distance ⁇ between the original and reduced meshes and reduce until that distance is met.
  • Each simplex in the simplical complex ⁇ ⁇ ⁇ may have a set of three-dimensional affine coefficients (slope ⁇ ⁇ and intercept ⁇ ⁇ ) as per (60).
  • the input dimension ⁇ (dimension of ⁇ ) and output dimension ⁇ (dimension of ⁇ ) will produce a PWA function with a slope matrix sized ⁇ ⁇ ⁇ ⁇ and offset vector ⁇ ⁇ ⁇ 1 ⁇ .
  • the boundaries of each simplex are defined by the ⁇ ⁇ 1 affine equations defining the H-notation of the simplex ⁇ ⁇ ⁇ ⁇ .
  • the overall number of coefficients may be ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ 2 ⁇ ⁇ ⁇ 1 ⁇ ⁇ ⁇ ⁇ ⁇ 2 ⁇ ⁇ ⁇ ⁇ 2 ⁇ , or the number of -60- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 coefficients increases linearly with the number of simplices ⁇ ⁇ .
  • the simplical complex may be stored in a tree structure whereby the search time complexity is ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ without warm starting and ⁇ 1 ⁇ with warm starting. [00214] Given an allotted amount of memory on a microcontroller, the number of simplices ⁇ ⁇ that can fit in memory can be predicted.
  • the number of simplices ⁇ ⁇ that can be used without running out of time can be predicted.
  • the smaller of ⁇ ⁇ and ⁇ ⁇ may be selected and used as the target number of simplices.
  • certain mesh reduction algorithms are used for online level-of-detail (LOD) modelling, and are designed to perform quickly, but may sacrifice some accuracy.
  • the MTPA PWA map is computed offline and the static map is loaded onto a controller (e.g., the motor controller 120). In such examples, fast computation may be less of a priority.
  • the mesh reduction algorithm used is vertex contraction, which can change the topological properties of the mesh, but joins surfaces, rather than separate surfaces.
  • the resulting piecewise map corresponding to a simplical complex output or provided by application of the mesh reduction algorithm may be employed by the motor controller 120 to use or solve the optimization problem (e.g., (7), (45), or (51)) to determine the target motor control parameter values (e.g., as described with respect to block 310 of FIG. 3).
  • the piecewise function may be stored in the memory 130 of the motor controller 120 and, to execute OERG block 215 of FIG. 2 and/or block 310 of FIG.
  • the motor controller 120 solves the piecewise function with the desired control parameter (torque and/or speed) as an input to the piecewise function.
  • the output or solution of the piecewise function may be, for example, the reference current ir* or, when OERG block 215 is integrated with the flux linkage map block 220, the reference flux ⁇ r * (see, e.g., FIG.2).
  • FEA analysis on a 65kW WRSM was performed.
  • the copper loss datapoints and pareto-optimal points in ( ⁇ ⁇ , ⁇ , ⁇ ⁇ ) space are shown in FIGS. 14A and 14B, respectively.
  • An initial simplical complex was generated that includes 44,640 simplices (see FIG.
  • a vertex contraction mesh reduction algorithm was applied to the initial simplical complex to reduce the number simplices to 4,463 in a first iteration (see FIG. 15B), to 446 -61- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 simplices in a second iteration (see FIG.15C), and to 45 simplices in a third iteration (see FIG. 15D). Accordingly, after three mesh simplification steps (e.g., three iterations of mesh reduction via vertex contraction), the number of simplices was reduced by 1000x to generate a simplical complex of 45 faces (see FIG.15D).
  • the motor controller 120 in any of its various configurations described herein (see, e.g., FIG. 2) is implemented as a set of instructions stored on a nontransitory computer readable medium, where the instructions are for execution by a processor.
  • the processor may be configured to (or may be connected to another device configured to) simulate a motor, power supply, and power switching network (simulating an arrangement similar to, for example, the arrangement in FIG. 2).
  • the processor through execution of the set of instructions, may be configured to monitor and control a motor where the motor is a simulated motor coupled to a simulated power supply via a simulated power switching network.
  • Example 1 A method, apparatus, and non-transitory computer-readable medium for motor control comprises: a power switching network configured to be coupled to a power supply and to a motor; and an electronic controller configured to: determine current values for the motor in a rotational reference frame, each current value associated with a dimension of a set of dimensions of the rotational reference frame; determine, based on a desired control parameter, a target motor control parameter value for each dimension of the set of dimensions of the rotational reference frame using an optimization cost function that considers motor speed, copper -62- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 loss, and core loss; and control the power switching network based on the current values and the target motor control parameter values.
  • Example 2 The method, apparatus, and non-transitory computer-readable medium according to Example 1, wherein the optimization cost function considers core loss by using a global core loss model with a matrix G of coefficients applicable regardless of motor speed of the motor.
  • Example 3 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 to 2, wherein the optimization cost function considers core loss by using a binned core loss model with a matrix of coefficients that depends on a speed of the motor.
  • Example 4 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 to 3, wherein the optimization cost function considers core loss by using a binned core loss model, wherein the binned core loss model is implemented as a piecewise function with M domains defined by motor speed, each of the M domains corresponding to a motor speed range and a matrix G of coefficients.
  • Example 5 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 to 4, wherein, a solution set of the optimization cost function is defined as a piecewise function with M domains, each of the M domains corresponding to a motor speed range and a motor torque range.
  • Example 6 The method, apparatus, and non-transitory computer-readable medium according to Example 5, wherein each of the M domains of the piecewise function corresponds to a simplex of a surface reconstructed from a set of pareto optimal points of datapoints derived from a set of inputs applied to the optimization cost function.
  • Example 7 The method, apparatus, and non-transitory computer-readable medium according to Example 5, wherein each of the M domains of the piecewise function corresponds to a simplex of a reduced simplical complex of a simplical complex, where simplical complex was formed from a set of pareto optimal points of datapoints derived from a set of inputs applied to the optimization cost function.
  • Example 8 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 6 or 7, wherein the simplical complex was formed from the set of pareto optimal points using at least one selected from a group of a surface reconstruction technique and a triangulation technique.
  • Example 9 The method, apparatus, and non-transitory computer-readable medium according to Example 7, wherein the reduced simplical complex was formed using a mesh reduction technique to reduce the simplical complex to a target number of simplices.
  • Example 10 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 5 to 9, wherein the piecewise function is stored in a memory of the electronic controller and, to determine the target motor control parameter value, the electronic controller solves the piecewise function with the desired control parameter as an input to the piecewise function.
  • Example 11 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 10, wherein the optimization cost function is associated with a first model corresponding to a first magnetic saturation level of the motor and with a second model corresponding to a second magnetic saturation level of the motor, wherein a first solution set of the optimization cost function for the first model is defined as a first piecewise function with domains, each of the domains corresponding to a respective motor speed range and a respective motor torque range, wherein a second solution set of the optimization cost function for the second model is defined as a second piecewise function with further domains, each of the further domains corresponding to a respective motor speed range and a respective motor torque range, and wherein, to determine the target motor control parameter value using the optimization cost function, the electronic controller is configured to: select a piecewise function from the first piecewise function or the second piecewise function to use based on a motor characteristic of the motor during operation; and solve the piecewise function with the desired control parameter as an input to the piecewise function
  • Example 12 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 11, wherein the optimization cost function is associated with a first model corresponding to a first magnetic saturation level of the motor and with a second model corresponding to a second magnetic saturation level of the motor, and wherein, to -64- Q B ⁇ 175073.00216 ⁇ 90201874.3 Attorney Docket No.: 175073.00216 determine the target motor control parameter value using the optimization cost function, the electronic controller is configured to: select a model from the first model or the second model to use based on a motor characteristic of the motor during operation; and use a solution of the model based on the desired control parameter as an input to the model.
  • Example 13 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 12, wherein, to determine current values for the motor in a rotational reference frame, the electronic controller is configured to: determine electrical operational characteristics of the motor in a stationary reference frame; determine a rotational position of the motor; and transform the electrical operational characteristics and the rotational position to the current values for the motor in the rotational reference frame.
  • Example 14 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 13, the electronic controller further configured to: determine, based on the current values, a flux linkage value for each dimension of the set of dimensions of the rotational reference frame, and wherein, to control the power switching network based on the current values, the electronic controller is configured to control the power switching network based on the flux linkage values determined from the current values.
  • Example 15 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 14, wherein the desired control parameter is a target torque value for the motor.
  • Example 16 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 15, wherein, to control the power switching network based on the current values and the target motor control parameter values, the electronic controller is configured to: generate a voltage command for each dimension of the set of dimensions of the rotational reference frame based on a difference between the target motor control parameter value and a motor parameter indicated by the current value for the dimension; transform the voltage commands in the rotational reference frame to the stationary reference frame; generate a pulse width modulated control signal for each dimension of the stationary reference frame to control the power switching network to drive a stator of the motor; and generate a rotor control signal to control driving of a rotor field winding.
  • Example 17 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 16, wherein, to control the power switching network based on the current values and the target motor control parameter values, the electronic controller is configured to: generate control signals in the stationary reference frame to drive the motor based on a difference between the target motor control parameter value and a motor parameter indicated by the current value for the dimension.
  • Example 18 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 17, wherein the motor is a wound field synchronous motor comprising at least three stator phases and at least one rotor field winding.
  • Example 19 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 18, wherein the power switching network includes an inverter switch bridge including a plurality of power switching elements, the inverter switch bridge configured to receive DC power and output AC power to windings of the stator based on pulse width modulated control signals from the electronic controller.
  • Example 20 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 19, further comprising a DC/DC converter configured to receive input DC power and to provide output DC power to at least one rotor field winding in accordance with a pulse width modulated rotor control signal from the electronic controller.
  • Example 21 The method, apparatus, and non-transitory computer-readable medium according to any of Examples 1 or 20, wherein the motor is at least one selected from the group of a wound field synchronous motor, a hybrid synchronous motor, a permanent magnet synchronous motor, an induction motor, a universal motor, or a reluctance motor.
  • the motor is at least one selected from the group of a wound field synchronous motor, a hybrid synchronous motor, a permanent magnet synchronous motor, an induction motor, a universal motor, or a reluctance motor.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

L'invention concerne des systèmes et des procédés de commande de moteur utilisant une génération de référence d'efficacité optimale. Un dispositif de commande électronique peut déterminer des valeurs de courant pour un moteur dans un cadre de référence de rotation. Chaque valeur de courant peut être associée à une dimension d'un ensemble de dimensions du cadre de référence de rotation. Le dispositif de commande électronique peut en outre déterminer, sur la base des valeurs de courant, une valeur de liaison de flux pour chaque dimension de l'ensemble de dimensions de la trame de référence de rotation à l'aide d'une fonction de coût d'optimisation qui considère la vitesse du moteur, la perte de cuivre et la perte de noyau. Le dispositif de commande électronique peut en outre déterminer une valeur de liaison de flux cible pour chaque dimension de l'ensemble de dimensions du cadre de référence de rotation. Le dispositif de commande électronique peut ensuite commander un réseau de commutation de puissance couplé entre une alimentation électrique et le moteur sur la base des valeurs de liaison de flux et des valeurs de liaison de flux cibles.
PCT/US2024/034328 2023-06-15 2024-06-17 Commande de moteur à l'aide d'une génération de référence d'efficacité optimale Ceased WO2024259419A2 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020267001514A KR20260046084A (ko) 2023-06-15 2024-06-17 최적 효율 레퍼런스 생성을 이용한 모터 제어
EP24824347.9A EP4728634A2 (fr) 2023-06-15 2024-06-17 Commande de moteur à l'aide d'une génération de référence d'efficacité optimale

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US202363521261P 2023-06-15 2023-06-15
US63/521,261 2023-06-15

Publications (2)

Publication Number Publication Date
WO2024259419A2 true WO2024259419A2 (fr) 2024-12-19
WO2024259419A3 WO2024259419A3 (fr) 2025-04-17

Family

ID=93852719

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2024/034328 Ceased WO2024259419A2 (fr) 2023-06-15 2024-06-17 Commande de moteur à l'aide d'une génération de référence d'efficacité optimale

Country Status (3)

Country Link
EP (1) EP4728634A2 (fr)
KR (1) KR20260046084A (fr)
WO (1) WO2024259419A2 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN120255356A (zh) * 2025-04-02 2025-07-04 临沂大学 一种轴向磁通电机的控制器优化方法、设备及介质
CN120993749A (zh) * 2025-09-25 2025-11-21 上海交通大学 一种降低损耗的高动态电机控制参数实时优化方法及系统

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7023168B1 (en) * 2004-09-13 2006-04-04 General Motors Corporation Field weakening motor control system and method
CA2659088C (fr) * 2006-07-24 2013-07-09 Kabushiki Kaisha Toshiba Systeme d'entrainement de moteur a flux variable
JP5167631B2 (ja) * 2006-11-30 2013-03-21 株式会社デンソー モータの制御方法及びそれを利用するモータ制御装置
US11159112B2 (en) * 2018-11-30 2021-10-26 The Trustees Of Columbia University In The City Of New York Systems and methods for high performance filtering techniques for sensorless direct position and speed estimation
US11233473B2 (en) * 2019-03-01 2022-01-25 Deere & Company Method and system for controlling a permanent magnet machine without a mechanical position sensor
EP4029576B1 (fr) * 2021-01-13 2024-12-18 Hydropool Inc. Piscine/spa pour nager sur place

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN120255356A (zh) * 2025-04-02 2025-07-04 临沂大学 一种轴向磁通电机的控制器优化方法、设备及介质
CN120993749A (zh) * 2025-09-25 2025-11-21 上海交通大学 一种降低损耗的高动态电机控制参数实时优化方法及系统

Also Published As

Publication number Publication date
WO2024259419A3 (fr) 2025-04-17
KR20260046084A (ko) 2026-04-06
EP4728634A2 (fr) 2026-04-22

Similar Documents

Publication Publication Date Title
Qu et al. Loss-minimizing flux level control of induction motor drives
Habetler et al. Stator resistance tuning in a stator-flux field-oriented drive using an instantaneous hybrid flux estimator
US12348161B2 (en) Motor control using piecewise affine model
WO2024259419A2 (fr) Commande de moteur à l'aide d'une génération de référence d'efficacité optimale
Perera Sensorless control of permanent-magnet synchronous motor drives
Nie et al. Deadbeat-direct torque and flux control for wound field synchronous machines
De Kock et al. Optimal torque control of synchronous machines based on finite-element analysis
Hua et al. Improved model‐predictive‐flux‐control strategy for three‐phase four‐switch inverter‐fed flux‐reversal permanent magnet machine drives
Merlyn et al. Review of control topologies for flux switching motor
Zhao Position/speed sensorless control for permanent-magnet synchronous machines
Dwivedi et al. Review on control strategies of permanent magnet-assisted synchronous reluctance motor drive
US20250007434A1 (en) System and method for controlling a motor
Qu et al. Minimizing losses of a synchronous reluctance motor drive taking into account core losses and magnetic saturation
Kumar et al. Modified direct torque control of three-phase induction motor drives with low ripple in flux and torque
Li et al. Non‐linear deadbeat direct torque and flux control for switched reluctance motor
Mohith et al. Comparative analysis of different control techniques for six-phase PMSM as an application to HEV
Vanamala et al. Speed control of PMSM drive using model predictive control based field oriented control
Casadei et al. Unified model and field oriented control algorithm for three-phase AC machines
CN118120145A (zh) 用于控制三相电机的方法、控制单元及其系统
Bramerdorfer et al. Multi-harmonic design and optimization of PMSMs
Tahim et al. Efficient implementation of high-fidelity models of IPMSM drive systems in offline and real-time simulators
İnan et al. Speed-Sensorless DTC of BLDC Motor with EKF-based Estimator Capable of Load Torque Estimation for Electric Vehicle
Correa et al. Sensorless control strategies for single-phase induction motor drive system
De Kock et al. Optimal torque control of interior permanent magnet synchronous machines in the full speed range
Iyer Six-Step Inverter-Fed Permanent Magnet Synchronous Motor and Brushless DC Motor Drives

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 24824347

Country of ref document: EP

Kind code of ref document: A2

WWE Wipo information: entry into national phase

Ref document number: 202547134691

Country of ref document: IN

WWP Wipo information: published in national office

Ref document number: 202547134691

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: 2024824347

Country of ref document: EP

Ref document number: 1020267001514

Country of ref document: KR

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2024824347

Country of ref document: EP

Effective date: 20260115

ENP Entry into the national phase

Ref document number: 2024824347

Country of ref document: EP

Effective date: 20260115

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 24824347

Country of ref document: EP

Kind code of ref document: A2

ENP Entry into the national phase

Ref document number: 2024824347

Country of ref document: EP

Effective date: 20260115

ENP Entry into the national phase

Ref document number: 2024824347

Country of ref document: EP

Effective date: 20260115

WWP Wipo information: published in national office

Ref document number: 2024824347

Country of ref document: EP