EP3806323A1 - Method and system for estimating an initial temperature of a rotor of an electric machine - Google Patents
Method and system for estimating an initial temperature of a rotor of an electric machine Download PDFInfo
- Publication number
- EP3806323A1 EP3806323A1 EP19202515.3A EP19202515A EP3806323A1 EP 3806323 A1 EP3806323 A1 EP 3806323A1 EP 19202515 A EP19202515 A EP 19202515A EP 3806323 A1 EP3806323 A1 EP 3806323A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- rotor
- current
- temperature
- module
- voltage
- 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.)
- Withdrawn
Links
Images
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
- H02P29/00—Arrangements for regulating or controlling electric motors, appropriate for both AC and DC motors
- H02P29/60—Controlling or determining the temperature of the motor or of the drive
- H02P29/66—Controlling or determining the temperature of the rotor
-
- 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
- H02P29/00—Arrangements for regulating or controlling electric motors, appropriate for both AC and DC motors
- H02P29/60—Controlling or determining the temperature of the motor or of the drive
- H02P29/66—Controlling or determining the temperature of the rotor
- H02P29/662—Controlling or determining the temperature of the rotor the rotor having permanent magnets
Definitions
- the invention relates to a method and a corresponding system for estimating an initial rotor temperature of a rotor of an electric machine, in particular of a permanent magnet synchronous machine, before an operating phase of this electric machine, wherein the electric machine is supplied with current via current paths by means of a power inverter and wherein via a module for carrying out a so-called d/q-transformation
- IPMSM interior permanent magnet synchronous motors
- the system requirements specify torque production commonly with high "safety integrity level” - ASIL and certain accuracy for the electric drive train.
- To compute the motor torque with required accuracy in electronic control unit the influence of temperature must be taken in account. It is because permanent magnets, such as NdFeB which are used to build rotor, have their parameters, remanence and coercivity depended on temperature. Therefore, magnet flux level depends on temperature and it will influence accuracy of torque calculation.
- control unit software In control unit software, a model that estimates rotor magnet temperature is usually implemented. The model output is used for corrections of computed torque due to the temperature effect on magnets. Estimated temperature is also used in another function which prevents rotor to overheat, otherwise magnets can be irreversibly demagnetized when certain critical temperature is reached.
- Such rotor temperature models or rotor temperature estimators may be based on different principles, i.e. thermal node model - here individual motor elements/parts are represented by nodes which are interconnected by thermal resistances and capacitances and are forming a network (Foster/Cauer configuration). The inputs are electrical losses developed in motor parts and node output is the temperature.
- the node model can be of arbitrary complexity.
- flux based estimator is being estimated using motor parameters and measured quantities, currents and voltage.
- flux value is known, then using magnet characteristic as function of temperature it is possible to determine rotor temperature. If also motor emf voltage is estimated also this is related to magnet temperature.
- thermal models for rotor temperature estimation i.e. neural networks etc.
- the rotor temperature estimator If the rotor temperature estimator is used in its various forms, its output must be properly initialized at the beginning of the operating phase. There are several drawbacks associated with initialization methods in prior art. If some other measured temperature is used, it is often a case, that there is a temperature difference between temperature of interest in rotor and other temperature. Especially when cooling effects are different between two spots. Only after very long time between last and new driving cycle elapsed, the temperatures may have equalized, and error will be small.
- stator temperature measurement might be not available because it requires sensor and processing circuits, or its safety level is not sufficient.
- Last known temperature method is only suitable when time between switch off and next driving cycle is very short. The cooling which could be still operated may reduce temperature and thus temperature difference is created.
- the initialization which takes into account the last known temperature and time while the electrical machine was switched off appears as best, but to measure the time comes at a cost. There must be present some real time clock capability on device operating even when the power is off, or device is in sleep mode.
- model If the model is initialized with the temperature value having big error, then model must be capable to eliminate this error and converge to real rotor temperature over the time. During this convergence time, the torque calculation precision is impacted.
- the electrical machine preferably is an Interior Permanent Magnet Synchronous Machine (IPMSM).
- IPMSM Interior Permanent Magnet Synchronous Machine
- U dq is the machine terminal voltage
- i dq is the machine current
- ⁇ dq stator flux
- ⁇ PM permanent magnet flux
- L dq motor inductance in direct and quadrature axis respectively. If the fluxes are substituted into d voltage equation, then this results in:
- ud Rs ⁇ Id ⁇ ⁇ ⁇ Lq ⁇ Iq + d dt Ld ⁇ Id + ⁇ PM
- L d inductance term is temperature invariant, but permanent magnet flux does depend on temperature and thus temperature will affect current derivative value.
- the voltage pulse signal is applied to one voltage input while the voltage at the other voltage input is simultaneously zero.
- the one voltage input is the d input
- the one current output is the corresponding d output and the other current output is the q output.
- the voltage pulse signal has pulse components of both polarities in equal parts, but only positive part is needed for temperature estimation.
- the maximum current or current gradient (difference between maximum current and previous sample) of the current response is considered at the determination of the initial rotor temperature.
- the rotor position is determined by a sensor and taken into account in estimating the initial rotor temperature.
- the initial rotor temperature is estimated via the rotor position, the current response and corresponding table values associating the rotor position and the current response with the initial rotor temperature.
- the module for outputting switching signals is a PWM module.
- the duration of the voltage pulse signal consists of several cycles or only one cycle of switching signal, in particular PWM cycle.
- the power inverter is realized as a computer-based power inverter, in which at least the module for carrying out a so-called d/q transformation is designed as a software module.
- the voltage pulse signal is a rectangular signal, a saw tooth signal or a mixed form of rectangular and saw tooth signal.
- the system comprises a power inverter for supplying current to the electric machine via current paths, wherein this power inverter comprises a module for carrying out a so-called d/q transformation, via which module, (i) a module for outputting switching signals for the electrical machine can be controlled by means of two voltages, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs, and (ii) the current flowing along the current paths can be read out at corresponding current outputs as currents, each referred to one of the two rotor referred coordinates d and q.
- the system further comprises a temperature estimating module, which is arranged to give a voltage pulse signal to at least one of the voltage inputs before or at the beginning of the operating phase of the electrical machine and to determine the initial rotor temperature from a current response at one or more of the current outputs to this defined voltage pulse signal.
- a temperature estimating module which is arranged to give a voltage pulse signal to at least one of the voltage inputs before or at the beginning of the operating phase of the electrical machine and to determine the initial rotor temperature from a current response at one or more of the current outputs to this defined voltage pulse signal.
- the system is arranged to carry out the aforementioned method.
- Fig. 1 shows a system 10 for estimating an initial rotor temperature of a rotor 12 of an electric machine 14 before an operating phase of said machine 14.
- the electric machine 14 can be installed in a power train of a vehicle and might be designed as a permanent magnet synchronous machine 16.
- the system 10 comprises a power inverter 18 for supplying current to the electric machine via current paths 20, 22, 24, wherein this power inverter 12 comprises a module 26 for carrying out a so-called d/q transformation, via which module 26 a module 28 for outputting switching signals for the electrical machine can be controlled by means of two voltages VD, VQ, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs 30, 32.
- the corresponding current components ID, IQ of the current flowing along the current paths 20, 22, 24 can be read out at corresponding current outputs 34, 36, each referred to one of the two rotor referred coordinates d and q.
- a temperature estimating module 38 of the system 10 which is arranged to give a voltage pulse signal 40 to one of the voltage inputs 24 before or at the beginning of the operating phase of the electrical machine and to determine the initial rotor temperature from a current response 42 at one or more of the current outputs 28 to this defined voltage pulse signal.
- a well-known standard electronic power inverter 18 with appropriate software can be used to implement the proposed method for temperature identification of the initial rotor temperature.
- relevant control software blocks are described as the d/q-transformation module 26 and the switching signal module 28, which preferably is a PWM module 44 (PWM: pulse width modulator).
- PWM pulse width modulator
- the PWM module 44 converts voltage control references to switching signal for the power device electrical machine 14 and will guarantee that reference voltage is realized on the output of the power inverter 18 according to available DC voltage.
- three or two measured currents are handled by transformation blocks performing Clark-Park conversion to DQ frame.
- the angular position of the rotor 12 is determined by a sensor 46.
- inverse Park-Clark transformation blocks are processing reference voltage to the switching signal module 28 /PWM module 44.
- Rotor position is also measured by a suitable feedback device i.e. resolver, encoder. Position is needed for Park conversion and, as it will be shown later, it is essential for correct implementation of novel idea of this patent.
- a voltage pulse signal 40 is applied to VD voltage, while VQ is zero, to excite current derivative as of equation [4].
- Typical duration of voltage pulse signal is several hundreds of microseconds i.e. from 300 ⁇ s to 600 ⁇ s and covers several control cycles with active PWM. If control cycle is, say 100 ⁇ s, it lasts from 3 to 6 cycles. Amplitude of pulse depends on available DC voltage, in general it is preferred to use as much voltage as possible from DC supply.
- Shape of the pulse signal 40 may be rectangular, or as shown in Fig.2 , it may continuously rise according to current to cover drop on stator resistance. Stator resistance is usually very small in this type of traction motors and could be also neglected.
- the switching signal module 28 / PWM module 44 must ensure that commanded voltage is realized on output with most accuracy even when the supply voltage would be dropping because of increasing current and weak supply. Voltage is a sequence of positive and negative pulse. Applying pulse with both polarities has a consequence that average torque produced by the motor (considering also cross coupling with Q axis) is zero over the entire pulse duration, in other words, the rotor will not change its position.
- the current response 42 is produced in D axis which carries the information about the temperature.
- An example of current is shown in Fig.3 for voltage pulse duration 500 ⁇ s positive and 500 ⁇ s negative. It is possible to observe that current peak (at the instant of voltage inversion) depends on temperature. In addition, also current gradient, computed as a difference of new and previous current measurement sample depends on temperature. Example temperature values investigated from 20°C to 150°C are shown.
- Fig.4 is shown cross section of example IMPSM for one pole.
- the stator 48 and rotor 12 laminations are visible, with magnets 50 buried in rotor 12 and winding in stator slots 52.
- FEM flux analysis was carried out.
- the rotor flux penetration to stator 48 is not constant on stator circumference, but it depends on relative rotor-stator position and is influenced by stator slots 52. That is, flux density in stator 48 and therefore saturation is position dependent. That is true for any interior permanent magnet motor.
- Fig.4 is shown rotor 12 in a position which corresponds to minimal current peak response.
- flux lines encircle two stator slots 52 causing certain degree of saturation.
- Fig.6 shows ID current peaks detected for rotor position change corresponding to one pole.
- the slotting effect is manifested, there are 6 peaks visible in current response corresponding to 6 stator slots 52 per rotor pole.
- the shape of waveform and number of peaks is only an example for this particular rotor construction and this particular number of stator slots 52, different shapes are possible for different motors/machines 14.
- Fig. 6 shows that shape repeats regularly for this example, but in reality, there can be irregularities in current response 42 if considered entire stator circumference and all motor poles. Magnets may not be magnetized all to the same level.
- Stator 48 and rotor 12 are manufactured with some mechanical tolerances. Ideally, rotor absolute position would have to be known, to perfectly compensate for such tolerances. If absolute position over one mechanical revolution is not available, then rotor electrical position (needed for control of motor) can also be used having some averaging effect. That depends on type of feedback device employed.
- Fig.7 also shows the temperature dependency of current response 42.
- maximal and minimal current response peak of Fig.6 is plotted over the temperature variation.
- 450A current peak may correspond to either i.e. 60°C or 84°C.
- 25°C error would be committed if position is not considered.
- the rotor 12 may be in any position with respect to stator 48. If position information is considered, the temperature detection accuracy is greatly improved.
- the identified rotor temperature by short current response measurement can be used to initialize a rotor temperature model.
- the method proposed in this application uses information about rotor position of the rotor 12 to derive position depended current response 42 to identify rotor magnet temperature of IMPSM by short instantaneous measurement. Duration can be minimal, just one pulse of approximately 1ms. No need for previous rotor temperature value or information about switch off time or cooling characteristics. The rotor temperature is found with accuracy of ⁇ 5K as results from practical test measurements, which guarantees excellent torque estimation right from the start. As the rotor temperature identification is independent of other temperature measurement, i.e. stator temperature, this information can be used for sensor plausibilization purpose in "safety" software. Involved quantities, current, voltage and position are already acquired in electronic control unit of EV inverter with high safety integrity level.
- the method can be also used for instantaneous correction of temperature from thermal model, if measurement and initialization is triggered in situations when speed is zero i.e. stopping at traffic light in case of EV application.
- the method can be also used for diagnostic purpose, to detect level of demagnetization of rotor magnets if a reference temperature is known by comparing the current response 42 with one referred to some reference temperature.
- software for machine control in the ECU as described by block diagram of Fig.1 is needed, capable of producing a test pulse signal 40 as of Fig.2 .
- software must measure ID current peak and/or current gradient (difference between current peak and previous current sample).
- Electrical machine 14 under test without cooling (if present) is put into thermal chamber with constant temperature. Electrical machine 14 is heated by thermal chamber until rotor 12 reaches desired temperature. If the rotor temperature is not available simply wait sufficiently long time to equalize the temperature in all machine parts. Possibly start with highest value of temperature in the chamber, such as supported by machine parts design i.e. 150°C.
- the rotor position 0 is established. Then voltage test pulse 40 is applied by control software.
- the pulse signal 40 may be just one, but better implementation is several pulse signals, i.e. ten.
- Current response 42 is recorded for each pulse signal 40 and then is averaged over number of pulse signals applied. Block diagram of procedure is depicted on Fig.8 , the left part is active, used at calibration.
- Position step depends on particular machine design, number of stator slots per pole. Only for example purpose, a suitable value may be 3°el. Again, the test sequence is applied. The procedure is repeated for one motor revolution 360°mech., or minimally for one pole span 180°el.
- the temperature in thermal chamber is decreased and new rotor temperature is established after sufficiently long time.
- the position change procedure is repeated over from 0 position. Measurements at several temperatures are needed, 4 to 6 values. A full data set is obtained, the current response in function of position and temperatures, similar to Fig.6 .
- Fig. 8 shows a block diagram of the approach to derive the initial rotor temperature from measured quantities.
- Blocks 54 - 58 are calibration steps I and blocks 54 - 60 are initialization steps II.
- This data set can be inverted by simple mathematical operation and "surface" as of Fig.9 is obtained.
- the missing points are calculated by interpolation.
- the eventual error caused by interpolation is small because relation of current to temperature is almost linear as shown in Fig.7 .
- the shape of this two-dimensional function depends on particular rotor construction as is provided here for example Fig.4 , 5 .
- Either ID current peak or ID current derivative can be used.
- Function can be stored in the form of table and used in next step.
- Fig. 9 shows the rotor temperature as function of position and current peak (interpolation surface example).
- a second implementation is "runtime" execution and initialization procedure.
- a test pulse is applied to detect current response (peak or gradient).
- a test pulse is applied in rotor D axis but with respect to stator this can be any position. Just one pulse of approximately 1ms is sufficient, but also two or more can be used and then average value is used for further processing. The actual or averaged position is also recorded.
- This process is described in Fig. 8 on left side as a first step.
- current value and position are subsequently used in table as of Fig. 9 and by interpolation rotor temperature is obtained.
- the temperature value is very accurate depending only on current measurement accuracy for a sample machine 14.
- Fig. 10 shows a block diagram of the initialization of a thermal model for the rotor temperature T.
- Block 62 comprises condition monitoring and triggering;
- block 64 comprises initial temperature estimation or identification.
- the resulting initial temperature T i is transferred to a closed loop rotor temperature estimating system comprising a rotor thermal model (Block 66), which estimates the rotor temperature T during operation.
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
- a module (28) for outputting switching signals for the electrical machine (14) can be controlled by means of two voltages VD, VQ, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs (30, 32), and
- the current flowing along the current paths (20, 22, 24) can be read out at two corresponding current outputs (34, 36) as currents ID, IQ, each referred to one of the two rotor referred coordinates d and q.
Description
- The invention relates to a method and a corresponding system for estimating an initial rotor temperature of a rotor of an electric machine, in particular of a permanent magnet synchronous machine, before an operating phase of this electric machine, wherein the electric machine is supplied with current via current paths by means of a power inverter and wherein via a module for carrying out a so-called d/q-transformation
- a module for outputting switching signals for the electrical machine can be controlled by means of two voltages, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs, and
- the current flowing along the current paths can be read out at two corresponding current outputs as currents, each referred to one of the two rotor referred coordinates d and q.
- In electric vehicle (EV) applications nowadays interior permanent magnet synchronous motors (IPMSM) are used for their efficiency, maximal power density and suitable speed range. In such applications, the system requirements specify torque production commonly with high "safety integrity level" - ASIL and certain accuracy for the electric drive train. To compute the motor torque with required accuracy in electronic control unit, the influence of temperature must be taken in account. It is because permanent magnets, such as NdFeB which are used to build rotor, have their parameters, remanence and coercivity depended on temperature. Therefore, magnet flux level depends on temperature and it will influence accuracy of torque calculation.
- In control unit software, a model that estimates rotor magnet temperature is usually implemented. The model output is used for corrections of computed torque due to the temperature effect on magnets. Estimated temperature is also used in another function which prevents rotor to overheat, otherwise magnets can be irreversibly demagnetized when certain critical temperature is reached.
- Such rotor temperature models or rotor temperature estimators may be based on different principles, i.e. thermal node model - here individual motor elements/parts are represented by nodes which are interconnected by thermal resistances and capacitances and are forming a network (Foster/Cauer configuration). The inputs are electrical losses developed in motor parts and node output is the temperature. The node model can be of arbitrary complexity.
- Another approach is flux based estimator - here flux is being estimated using motor parameters and measured quantities, currents and voltage. When flux value is known, then using magnet characteristic as function of temperature it is possible to determine rotor temperature. If also motor emf voltage is estimated also this is related to magnet temperature.
- Yet another approach is to use a "Time constant model", here the stator temperature is measured and then, using a table based information obtained during calibration about final rotor temperature and possible time constant for heating and cooling, the rotor temperature can be estimated.
- There may be other implementations of thermal models for rotor temperature estimation - i.e. neural networks etc.
- More innovative approaches require injection of some sort of excitation signal or test voltage or current. Observing as feedback signals are modified due to machine second order properties like saturation depending on temperature it is possible to derive rotor temperature online.
- Most thermal models, require an initial value of the temperature, the initial rotor temperature, to start with, unless the algorithm identifies temperature instantaneously. Especially initial rotor temperature is essential for correct torque computation at the beginning of driving. In prior art, this initial value of temperature is obtained in several ways:
- a) some other measured value i.e. stator temperature or ambient temperature, if available;
- b) last known temperature from model (saved before switching off the unit); or
- c) calculated value from last know temperature as of point b) considering also cooling time constant and time interval during which the unit was switched off.
- If the rotor temperature estimator is used in its various forms, its output must be properly initialized at the beginning of the operating phase. There are several drawbacks associated with initialization methods in prior art. If some other measured temperature is used, it is often a case, that there is a temperature difference between temperature of interest in rotor and other temperature. Especially when cooling effects are different between two spots. Only after very long time between last and new driving cycle elapsed, the temperatures may have equalized, and error will be small.
- In some low-cost applications however, the stator temperature measurement might be not available because it requires sensor and processing circuits, or its safety level is not sufficient. Last known temperature method is only suitable when time between switch off and next driving cycle is very short. The cooling which could be still operated may reduce temperature and thus temperature difference is created.
- The initialization, which takes into account the last known temperature and time while the electrical machine was switched off appears as best, but to measure the time comes at a cost. There must be present some real time clock capability on device operating even when the power is off, or device is in sleep mode.
- If the model is initialized with the temperature value having big error, then model must be capable to eliminate this error and converge to real rotor temperature over the time. During this convergence time, the torque calculation precision is impacted.
- It is an object of the invention to provide means allowing a quick and fairly accurate estimation of the initial rotor temperature.
- This object is achieved by a method and a system for estimating an initial rotor temperature having the features according to the respective independent claims. Advantageous implementations of the invention are the subject matter of the dependent claims, of the description and of the figures.
- According to several aspects of the inventive method for estimating an initial rotor temperature of a rotor of an electric machine before an operating phase of this electric machine, in which the electric machine is supplied with current via current paths by means of a power inverter and wherein by use of a module for carrying out a so-called d/q-transformation
- (i) a module for outputting switching signals for the electrical machine can be controlled by means of two voltages, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs, and
- (ii) the current flowing along the current paths can be read out at two corresponding current outputs as currents, each referred to one of the two rotor referred coordinates d and q, it is intended before or at the beginning of the operating phase a voltage pulse signal is given to at least one of the voltage inputs and the initial rotor temperature is determined from a current response at one or more of the current outputs to this defined voltage pulse signal.
- The electrical machine preferably is an Interior Permanent Magnet Synchronous Machine (IPMSM). The basis for an explanation of the above introduced approach are the following equations of IPMSM based on flux components in direct and quadrature axis:
-
-
- Clearly the temperature effect is hidden in the derivative term on the right side. Ld inductance term is temperature invariant, but permanent magnet flux does depend on temperature and thus temperature will affect current derivative value.
- For simple understanding, the entire term in derivative is rough some equivalent inductance multiplied by a current increase L*ΔI and applied voltage Ud is constant. If the temperature is low, magnet remanence is higher, magnet flux will saturate stator core stronger, equivalent inductance will be lower, thus ΔI current term must be higher. Vice versa, if the magnet temperature is higher, its flux will be lower, saturation degree is smaller, equivalent inductance is higher and thus ΔI term must be smaller to match constant voltage term on the left side.
- In IMPSM machines for EV applications such temperature dependent saturation effect is very strong. Due to the continuous strive of designs for small motor size, the magnetic circuit is operated with highest possible flux density.
- According to a preferred embodiment of the invention the voltage pulse signal is applied to one voltage input while the voltage at the other voltage input is simultaneously zero. The one voltage input is the d input, the other voltage input is the q input (with the voltage VQ at the other voltage input being VQ = 0). The one current output is the corresponding d output and the other current output is the q output.
- Preferably the voltage pulse signal has pulse components of both polarities in equal parts, but only positive part is needed for temperature estimation.
- According to another preferred embodiment of the invention the maximum current or current gradient (difference between maximum current and previous sample) of the current response is considered at the determination of the initial rotor temperature.
- According to yet another preferred embodiment of the invention the rotor position is determined by a sensor and taken into account in estimating the initial rotor temperature.
- Preferably the initial rotor temperature is estimated via the rotor position, the current response and corresponding table values associating the rotor position and the current response with the initial rotor temperature.
- Preferably the module for outputting switching signals is a PWM module.
- According to another preferred embodiment of the invention the duration of the voltage pulse signal consists of several cycles or only one cycle of switching signal, in particular PWM cycle.
- According to yet another preferred embodiment of the invention the power inverter is realized as a computer-based power inverter, in which at least the module for carrying out a so-called d/q transformation is designed as a software module. Finally, it is advantageous that the voltage pulse signal is a rectangular signal, a saw tooth signal or a mixed form of rectangular and saw tooth signal.
- According to several aspects of the system for estimating an initial rotor temperature of a rotor of an electric machine before an operating phase of this electric machine, the system comprises a power inverter for supplying current to the electric machine via current paths, wherein this power inverter comprises a module for carrying out a so-called d/q transformation, via which module, (i) a module for outputting switching signals for the electrical machine can be controlled by means of two voltages, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs, and (ii) the current flowing along the current paths can be read out at corresponding current outputs as currents, each referred to one of the two rotor referred coordinates d and q. The system further comprises a temperature estimating module, which is arranged to give a voltage pulse signal to at least one of the voltage inputs before or at the beginning of the operating phase of the electrical machine and to determine the initial rotor temperature from a current response at one or more of the current outputs to this defined voltage pulse signal.
- According to another preferred embodiment of the system according to the invention, the system is arranged to carry out the aforementioned method.
- Further features of the invention are apparent from the claims, the figure and the description of the figure. All of the features and feature combinations mentioned above in the description as well as the features and feature combinations mentioned below in the description of the figure and/or shown in the figure alone are usable not only in the respectively specified combination, but also in other combinations or alone.
- Now, the invention is explained in more detail based on a preferred embodiment as well as with reference to the attached drawings.
- In the drawings:
-
Fig. 1 shows a control block diagram with an electric machine and a system for estimating an initial rotor temperature of a rotor of said electric machine, -
Fig. 2 shows the shape of a voltage pulse signal, -
Fig. 3 shows the dependency of an ID current peak on the magnet temperature of the rotor, -
Fig. 4 shows sector like pieces of a rotor and a stator of the electric machine in a position where ID current is minimal, -
Fig. 5 shows sector like pieces of a rotor and a stator of the electric machine in a position where ID current is maximal, -
Fig. 6 shows the ID current peak dependency on the rotor position and different rotor temperatures, -
Fig. 7 shows the ID current peak dependency on the rotor temperature and rotor position related max-min deviation, and -
Fig. 8 shows a block diagram of the approach to derive the initial rotor temperature from measured quantities. -
Fig. 9 shows the rotor temperature as function of position and current peak (interpolation surface example), and -
Fig. 10 shows a block diagram of the initialization of a thermal model for the rotor temperature. -
Fig. 1 shows asystem 10 for estimating an initial rotor temperature of arotor 12 of an electric machine 14 before an operating phase of said machine 14. The electric machine 14 can be installed in a power train of a vehicle and might be designed as a permanent magnet synchronous machine 16. Thesystem 10 comprises apower inverter 18 for supplying current to the electric machine viacurrent paths 20, 22, 24, wherein thispower inverter 12 comprises amodule 26 for carrying out a so-called d/q transformation, via which module 26 a module 28 for outputting switching signals for the electrical machine can be controlled by means of two voltages VD, VQ, each related to one of two rotor referred coordinates d and q applied to one of twocorresponding voltage inputs 30, 32. On the other hand the corresponding current components ID, IQ of the current flowing along thecurrent paths 20, 22, 24 can be read out at correspondingcurrent outputs - For estimating the initial rotor temperature Ti of the rotor, a
temperature estimating module 38 of thesystem 10, which is arranged to give avoltage pulse signal 40 to one of the voltage inputs 24 before or at the beginning of the operating phase of the electrical machine and to determine the initial rotor temperature from acurrent response 42 at one or more of the current outputs 28 to this defined voltage pulse signal. - In other words: A well-known standard
electronic power inverter 18 with appropriate software can be used to implement the proposed method for temperature identification of the initial rotor temperature. OnFig.1 , relevant control software blocks are described as the d/q-transformation module 26 and the switching signal module 28, which preferably is a PWM module 44 (PWM: pulse width modulator). The PWM module 44 converts voltage control references to switching signal for the power device electrical machine 14 and will guarantee that reference voltage is realized on the output of thepower inverter 18 according to available DC voltage. On bottom path of the d/q transformation module 26, three or two measured currents are handled by transformation blocks performing Clark-Park conversion to DQ frame. The angular position of therotor 12 is determined by asensor 46. - On top path of the d/
q transformation module 26, inverse Park-Clark transformation blocks are processing reference voltage to the switching signal module 28 /PWM module 44. Rotor position is also measured by a suitable feedback device i.e. resolver, encoder. Position is needed for Park conversion and, as it will be shown later, it is essential for correct implementation of novel idea of this patent. - A
voltage pulse signal 40 is applied to VD voltage, while VQ is zero, to excite current derivative as of equation [4]. Typical duration of voltage pulse signal is several hundreds of microseconds i.e. from 300µs to 600µs and covers several control cycles with active PWM. If control cycle is, say 100µs, it lasts from 3 to 6 cycles. Amplitude of pulse depends on available DC voltage, in general it is preferred to use as much voltage as possible from DC supply. Shape of thepulse signal 40 may be rectangular, or as shown inFig.2 , it may continuously rise according to current to cover drop on stator resistance. Stator resistance is usually very small in this type of traction motors and could be also neglected. The switching signal module 28 / PWM module 44 must ensure that commanded voltage is realized on output with most accuracy even when the supply voltage would be dropping because of increasing current and weak supply. Voltage is a sequence of positive and negative pulse. Applying pulse with both polarities has a consequence that average torque produced by the motor (considering also cross coupling with Q axis) is zero over the entire pulse duration, in other words, the rotor will not change its position. - When the
voltage pulse signal 40 is applied as described above, thecurrent response 42 is produced in D axis which carries the information about the temperature. An example of current is shown inFig.3 for voltage pulse duration 500µs positive and 500µs negative. It is possible to observe that current peak (at the instant of voltage inversion) depends on temperature. In addition, also current gradient, computed as a difference of new and previous current measurement sample depends on temperature. Example temperature values investigated from 20°C to 150°C are shown. - Current level reached during test pulse is compatible with EV types of IPMSMs and inverters. Temperature dependency effect is strong as mentioned previously and can be several hundred amps. In general, the
voltage pulse signal 40 should be kept as short as possible to limitcurrent peak response 42 to such levels where temperature effect caused current difference is well detected by processing hardware of electronic control unit. - The current level however will not be the same when rotor position of the
rotor 12 changes. InFig.4 is shown cross section of example IMPSM for one pole. Thestator 48 androtor 12 laminations are visible, withmagnets 50 buried inrotor 12 and winding instator slots 52. On such arrangement FEM flux analysis was carried out. The rotor flux penetration tostator 48 is not constant on stator circumference, but it depends on relative rotor-stator position and is influenced bystator slots 52. That is, flux density instator 48 and therefore saturation is position dependent. That is true for any interior permanent magnet motor. OnFig.4 is shownrotor 12 in a position which corresponds to minimal current peak response. Here flux lines encircle twostator slots 52 causing certain degree of saturation. - On
Fig.5 is again shownrotor 12 in different position wherecurrent response 42 is maximal, flux lines encircle just onestator slot 52 suggesting that saturation is higher as in previous case. If current response 42 (peak and gradient) depends on position, it influences heavily temperature detection accuracy. -
Fig.6 shows ID current peaks detected for rotor position change corresponding to one pole. Clearly the slotting effect is manifested, there are 6 peaks visible in current response corresponding to 6stator slots 52 per rotor pole. The shape of waveform and number of peaks is only an example for this particular rotor construction and this particular number ofstator slots 52, different shapes are possible for different motors/machines 14. - The main point here is position dependency.
Fig. 6 shows that shape repeats regularly for this example, but in reality, there can be irregularities incurrent response 42 if considered entire stator circumference and all motor poles. Magnets may not be magnetized all to the same level.Stator 48 androtor 12 are manufactured with some mechanical tolerances. Ideally, rotor absolute position would have to be known, to perfectly compensate for such tolerances. If absolute position over one mechanical revolution is not available, then rotor electrical position (needed for control of motor) can also be used having some averaging effect. That depends on type of feedback device employed. -
Fig.7 also shows the temperature dependency ofcurrent response 42. Here maximal and minimal current response peak ofFig.6 is plotted over the temperature variation. Clearly, if say 450A current peak is detected, it may correspond to either i.e. 60°C or 84°C. Such big, 25°C error would be committed if position is not considered. - At switching on of electronic control unit, the
rotor 12 may be in any position with respect tostator 48. If position information is considered, the temperature detection accuracy is greatly improved. - The identified rotor temperature by short current response measurement can be used to initialize a rotor temperature model.
- The method proposed in this application uses information about rotor position of the
rotor 12 to derive position dependedcurrent response 42 to identify rotor magnet temperature of IMPSM by short instantaneous measurement. Duration can be minimal, just one pulse of approximately 1ms. No need for previous rotor temperature value or information about switch off time or cooling characteristics. The rotor temperature is found with accuracy of ±5K as results from practical test measurements, which guarantees excellent torque estimation right from the start. As the rotor temperature identification is independent of other temperature measurement, i.e. stator temperature, this information can be used for sensor plausibilization purpose in "safety" software. Involved quantities, current, voltage and position are already acquired in electronic control unit of EV inverter with high safety integrity level. - The method can be also used for instantaneous correction of temperature from thermal model, if measurement and initialization is triggered in situations when speed is zero i.e. stopping at traffic light in case of EV application. The method can be also used for diagnostic purpose, to detect level of demagnetization of rotor magnets if a reference temperature is known by comparing the
current response 42 with one referred to some reference temperature. - Corresponding embodiments of the approach can be divided into two procedures:
- identification/measurement of position related current response relation with temperature for particular IPMSM during development, in other words, the calibration
- execution of short test measurement according to this method at ECU switch on, or in other particular instants when such measurement is permitted, using the calibration information to convert the current response value to temperature and thermal model initialization or plausibilization.
- In a first implementation, software for machine control in the ECU as described by block diagram of
Fig.1 is needed, capable of producing atest pulse signal 40 as ofFig.2 . In addition, software must measure ID current peak and/or current gradient (difference between current peak and previous current sample). - Electrical machine 14 under test, without cooling (if present) is put into thermal chamber with constant temperature. Electrical machine 14 is heated by thermal chamber until
rotor 12 reaches desired temperature. If the rotor temperature is not available simply wait sufficiently long time to equalize the temperature in all machine parts. Possibly start with highest value of temperature in the chamber, such as supported by machine parts design i.e. 150°C. - The
rotor position 0 is established. Thenvoltage test pulse 40 is applied by control software. Thepulse signal 40 may be just one, but better implementation is several pulse signals, i.e. ten.Current response 42 is recorded for eachpulse signal 40 and then is averaged over number of pulse signals applied. Block diagram of procedure is depicted onFig.8 , the left part is active, used at calibration. - A new rotor position is established. Position step depends on particular machine design, number of stator slots per pole. Only for example purpose, a suitable value may be 3°el. Again, the test sequence is applied. The procedure is repeated for one motor revolution 360°mech., or minimally for one
pole span 180°el. - With this procedure, a current response is measured for different positions at constant temperature and data as one waveform in
Fig.6 are obtained. - The temperature in thermal chamber is decreased and new rotor temperature is established after sufficiently long time. The position change procedure is repeated over from 0 position. Measurements at several temperatures are needed, 4 to 6 values. A full data set is obtained, the current response in function of position and temperatures, similar to
Fig.6 . -
Fig. 8 shows a block diagram of the approach to derive the initial rotor temperature from measured quantities. Inblock 54 voltage pulse signal(s) 40 is/are produced; in block 56 a current response (a current peak response IP) is measured; inblock 58 the angular position of therotor 12 is measured; and inblock 60 the rotor temperature is determines by use of a table + interpolation (T = f(IP, position). Blocks 54 - 58 are calibration steps I and blocks 54 - 60 are initialization steps II. - This data set can be inverted by simple mathematical operation and "surface" as of
Fig.9 is obtained. The missing points are calculated by interpolation. The eventual error caused by interpolation is small because relation of current to temperature is almost linear as shown inFig.7 . The shape of this two-dimensional function depends on particular rotor construction as is provided here for exampleFig.4 ,5 . Either ID current peak or ID current derivative can be used. Function can be stored in the form of table and used in next step.Fig. 9 shows the rotor temperature as function of position and current peak (interpolation surface example). - A second implementation is "runtime" execution and initialization procedure. At the ECU start or whenever software determines there is suitable condition i.e. zero speed, a test pulse is applied to detect current response (peak or gradient). A test pulse is applied in rotor D axis but with respect to stator this can be any position. Just one pulse of approximately 1ms is sufficient, but also two or more can be used and then average value is used for further processing. The actual or averaged position is also recorded. This process is described in
Fig. 8 on left side as a first step. In a second step, on the right side, current value and position are subsequently used in table as ofFig. 9 and by interpolation rotor temperature is obtained. The temperature value is very accurate depending only on current measurement accuracy for a sample machine 14. - The temperature value is further used to initialize rotor temperature model as shown in
Fig.10 or for checking other sensor measured value range, so called plausibilization in safety software.Fig. 10 shows a block diagram of the initialization of a thermal model for the rotortemperature T. Block 62 comprises condition monitoring and triggering; block 64 comprises initial temperature estimation or identification. The resulting initial temperature Ti is transferred to a closed loop rotor temperature estimating system comprising a rotor thermal model (Block 66), which estimates the rotor temperature T during operation.
Claims (12)
- Method for estimating an initial rotor temperature of a rotor (12) of an electric machine (14), in particular of a permanent magnet synchronous machine (16), before an operating phase of this electric machine (14), wherein the electric machine (14) is supplied with current via current paths (20, 22, 24) by means of a power inverter (18) and wherein via a module (26) for carrying out a so-called d/q-transformation- a module (28) for outputting switching signals for the electrical machine (14) can be controlled by means of two voltages VD, VQ, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs (30, 32), and- the current flowing along the current paths (20, 22, 24) can be read out at two corresponding current outputs (34, 36) as currents ID, IQ, each referred to one of the two rotor referred coordinates d and q,characterized in that before or at the beginning of the operating phase a voltage pulse signal (40) is given to at least one of the voltage inputs (30) and the initial rotor temperature is determined from a current response to this defined voltage pulse signal (40) at one or more of the current outputs (34).
- Method according to claim 1, characterized in that the voltage pulse signal (40) is applied to one voltage input (30) while the voltage VQ at the other voltage input (32) is simultaneously zero.
- Method according to claim 1 or 2, characterized in that the voltage pulse signal (40) comprises pulse components of both polarities in equal parts or only a positive pulse component.
- Method according to any one of claims 1 to 3, characterized in that a maximum current of the current response is considered at the determination of the initial rotor temperature.
- Method according to any of claims 1 to 4, characterized in that the rotor position is determined by a sensor (46) and taken into account in estimating the initial rotor temperature.
- Method according to claim 5, characterized in that the initial rotor temperature is estimated via the rotor position, the current response and corresponding table values associating the rotor position and the current response with the initial rotor temperature.
- Method according to any one of claims 1 to 6, characterized in that the module for outputting switching signals is a PWM module (44).
- Method according to any one of claims 1 to 7, characterized in that the duration of the voltage pulse signal (40) is just one switching cycle or longer than a switching cycle of the switching signal output module (28), in particular the PWM cycle.
- Method according to one of the claims 1 to 8, characterized in that the power inverter (18) is realized as a computer-based power inverter, in which at least the module for carrying out a so-called d/q transformation (26) is designed as a software module.
- Method according to one of claims 1 to 9, characterized in that the voltage pulse signal (40) is a rectangular signal, a saw tooth signal or a mixed form of rectangular and saw tooth signal.
- System (10) for estimating an initial rotor temperature of a rotor (12) of an electric machine (14), in particular a permanent magnet synchronous machine (16), before an operating phase of this electric machine (14), wherein the system (10) comprises a power inverter (18) for supplying current to the electric machine (14) via current paths (20, 22, 24), wherein this power inverter (18) comprises a module (26) for carrying out a so-called d/q transformation, via which module (26),- a module (28) for outputting switching signals for the electrical machine (14) can be controlled by means of two voltages VD, VQ, each related to one of two rotor referred coordinates d and q applied to one of two corresponding voltage inputs (30, 32), and- the current flowing along the current paths (20, 22, 24) can be read out at corresponding current outputs (34, 36) as currents ID, IQ, each referred to one of the two rotor referred coordinates d and q,characterized by a temperature estimating module (38), which is arranged to give a voltage pulse signal (40) to at least one of the voltage inputs (30) before or at the beginning of the operating phase of the electrical machine (14) and to determine the initial rotor temperature from a current response (42) to this defined voltage pulse signal (40) at one or more of the current outputs (34).
- System according to claim 11, characterized in that said system (10) is arranged to carry out the method according to one of claims 1 to 9.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP19202515.3A EP3806323A1 (en) | 2019-10-10 | 2019-10-10 | Method and system for estimating an initial temperature of a rotor of an electric machine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP19202515.3A EP3806323A1 (en) | 2019-10-10 | 2019-10-10 | Method and system for estimating an initial temperature of a rotor of an electric machine |
Publications (1)
Publication Number | Publication Date |
---|---|
EP3806323A1 true EP3806323A1 (en) | 2021-04-14 |
Family
ID=68280731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP19202515.3A Withdrawn EP3806323A1 (en) | 2019-10-10 | 2019-10-10 | Method and system for estimating an initial temperature of a rotor of an electric machine |
Country Status (1)
Country | Link |
---|---|
EP (1) | EP3806323A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102021202982A1 (en) | 2021-03-26 | 2022-09-29 | Zf Friedrichshafen Ag | Method and device for determining a rotor temperature value for an electrical machine |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130234446A1 (en) * | 2010-11-10 | 2013-09-12 | Kokusan Denki Co., Ltd | Control device for rotary electrical machine |
US20150381090A1 (en) * | 2014-06-26 | 2015-12-31 | Nidec Motor Corporation | Sensorless system and method for determining motor angle at zero or low speeds |
-
2019
- 2019-10-10 EP EP19202515.3A patent/EP3806323A1/en not_active Withdrawn
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130234446A1 (en) * | 2010-11-10 | 2013-09-12 | Kokusan Denki Co., Ltd | Control device for rotary electrical machine |
US20150381090A1 (en) * | 2014-06-26 | 2015-12-31 | Nidec Motor Corporation | Sensorless system and method for determining motor angle at zero or low speeds |
Non-Patent Citations (4)
Title |
---|
GANCHEV M ET AL: "Sensitivity and robustness aspects of sensorless rotor temperature estimation technique for Permanent Magnet Synchronous Motor", POWER ELECTRONICS AND MOTION CONTROL CONFERENCE (EPE/PEMC), 2012 15TH INTERNATIONAL, IEEE, 4 September 2012 (2012-09-04), pages LS3a - 1.2, XP032312028, ISBN: 978-1-4673-1970-6, DOI: 10.1109/EPEPEMC.2012.6397426 * |
LAZCANO U ET AL: "DC voltage signal injection for temperature monitoring in very small resistance induction machines", ELECTRICAL SYSTEMS FOR AIRCRAFT, RAILWAY AND SHIP PROPULSION (ESARS), 2012, IEEE, 16 October 2012 (2012-10-16), pages 1 - 6, XP032282798, ISBN: 978-1-4673-1370-4, DOI: 10.1109/ESARS.2012.6387467 * |
REIGOSA DAVID DIAZ ET AL: "Permanent-Magnet Temperature Estimation in PMSMs Using Pulsating High-Frequency Current Injection", IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, vol. 51, no. 4, 1 July 2015 (2015-07-01), pages 3159 - 3168, XP011663199, ISSN: 0093-9994, [retrieved on 20150715], DOI: 10.1109/TIA.2015.2404922 * |
REIGOSA DAVID ET AL: "Comparative Analysis of BEMF and Pulsating High-Frequency Current Injection Methods for PM Temperature Estimation in PMSMs", IEEE TRANSACTIONS ON POWER ELECTRONICS, INSTITUTE OF ELECTRICAL AND ELECTRONICS ENGINEERS, USA, vol. 32, no. 5, 1 May 2017 (2017-05-01), pages 3691 - 3699, XP011640391, ISSN: 0885-8993, [retrieved on 20170202], DOI: 10.1109/TPEL.2016.2592478 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102021202982A1 (en) | 2021-03-26 | 2022-09-29 | Zf Friedrichshafen Ag | Method and device for determining a rotor temperature value for an electrical machine |
US11742727B2 (en) | 2021-03-26 | 2023-08-29 | Zf Friedrichshafen Ag | Method and device for determining a rotor temperature value for an electric motor |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ofori et al. | A pulse-injection-based sensorless position estimation method for a switched reluctance machine over a wide speed range | |
Armando et al. | Experimental identification of the magnetic model of synchronous machines | |
EP0420501B1 (en) | Method and apparatus for detecting the rotor position of a brushless DC motor | |
EP1535389B1 (en) | Motor drive control | |
CN105141213B (en) | Method for determining the position of the rotor of polyphase machine | |
Antonello et al. | Benefits of direct phase voltage measurement in the rotor initial position detection for permanent-magnet motor drives | |
EP2258043B1 (en) | Sensorless control of salient-pole machines | |
JPH0951700A (en) | Controlling device of rotary electric machine | |
Wang et al. | Online multiparameter identification method for sensorless control of SPMSM | |
Armando et al. | Flux linkage maps identification of synchronous AC motors under controlled thermal conditions | |
Wu et al. | A novel and simple torque ripple minimization method of synchronous reluctance machine based on torque function method | |
JP2003180094A (en) | Method for estimating pole position of pm motor | |
Yan et al. | Initial rotor position estimation and sensorless direct torque control of surface-mounted permanent magnet synchronous motors considering saturation saliency | |
Musak et al. | Novel methods for parameters investigation of PM synchronous motors | |
EP3806323A1 (en) | Method and system for estimating an initial temperature of a rotor of an electric machine | |
Kim et al. | Torque ripple improvement for interior permanent magnet synchronous motor considering parameters with magnetic saturation | |
Kuai et al. | Sensorless control of three phase switched reluctance motor drives using an approximate inductance model | |
Yan et al. | Modeling and simulation of direct torque controlled PMSM drive system incorporating structural and saturation saliencies | |
Hall et al. | Dynamic magnetic model identification of permanent magnet synchronous machines | |
Chowdhury et al. | Parameter characterization and performance prediction of permanent magnet synchronous motor for high performance applications | |
Xu et al. | Online estimation of the stator resistances of a PMSM with consideration of magnetic saturation | |
Armando et al. | Experimental methods for synchronous machines evaluation by an accurate magnetic model identification | |
JP2009526512A (en) | Method and apparatus for determining torque of power equipment | |
JP3627426B2 (en) | Control device for rotating electrical machine | |
Feuersänger et al. | Initial rotor position identification in medium voltage synchronous machines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION HAS BEEN PUBLISHED |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20211014 |
|
RBV | Designated contracting states (corrected) |
Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
RAP3 | Party data changed (applicant data changed or rights of an application transferred) |
Owner name: VITESCO TECHNOLOGIES GMBH |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN |
|
18D | Application deemed to be withdrawn |
Effective date: 20211015 |