JP7324630B2 - Motor magnet temperature estimation method and magnet temperature estimation device - Google Patents

Motor magnet temperature estimation method and magnet temperature estimation device Download PDF

Info

Publication number
JP7324630B2
JP7324630B2 JP2019128757A JP2019128757A JP7324630B2 JP 7324630 B2 JP7324630 B2 JP 7324630B2 JP 2019128757 A JP2019128757 A JP 2019128757A JP 2019128757 A JP2019128757 A JP 2019128757A JP 7324630 B2 JP7324630 B2 JP 7324630B2
Authority
JP
Japan
Prior art keywords
magnet
motor
flux linkage
current
estimating
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.)
Active
Application number
JP2019128757A
Other languages
Japanese (ja)
Other versions
JP2021016226A (en
Inventor
健介 佐々木
崇 加藤
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.)
Renault SAS
Original Assignee
Renault SAS
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 Renault SAS filed Critical Renault SAS
Priority to JP2019128757A priority Critical patent/JP7324630B2/en
Publication of JP2021016226A publication Critical patent/JP2021016226A/en
Application granted granted Critical
Publication of JP7324630B2 publication Critical patent/JP7324630B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Electric Motors In General (AREA)
  • Control Of Ac Motors In General (AREA)

Description

本発明は、モータの磁石温度の推定方法、及び、磁石温度の推定装置に関するものである。 The present invention relates to a motor magnet temperature estimation method and a magnet temperature estimation device.

モータの制御精度を向上させるために、回転中のモータの磁石温度を推定し、推定した磁石温度を用いてモータを制御する技術が検討されている。例えば、特許文献1に開示されている技術によれば、コイルエンドに磁束センサを設け、ロータからの漏洩する磁石磁束を磁束センサによって検出し、予め記憶している漏洩する磁石磁束と磁石温度との関係から磁石温度が推定されている。 In order to improve the control accuracy of a motor, a technique of estimating the magnet temperature of a rotating motor and controlling the motor using the estimated magnet temperature is being studied. For example, according to the technique disclosed in Patent Document 1, a magnetic flux sensor is provided at the coil end, the magnetic flux sensor detects the magnetic flux leaking from the rotor, and the leaked magnetic flux and the magnet temperature are stored in advance. The magnet temperature is estimated from the relationship of

特開2004-222387号公報Japanese Patent Application Laid-Open No. 2004-222387

特許文献1に開示されている技術によれば、磁石磁束に加えてコイルに負荷電流が流れることによっても負荷磁束が発生しており、コイルの近傍に存在する磁束センサが負荷磁束と磁石磁束との合計磁束を検出してしまうことがある。そのため、磁石磁束の測定精度が低下してしまい、磁石温度の推定精度が低下するおそれがあるという課題があった。 According to the technique disclosed in Patent Document 1, the load magnetic flux is generated not only by the magnet magnetic flux but also by the load current flowing through the coil, and the magnetic flux sensor existing near the coil detects the load magnetic flux and the magnet magnetic flux. may detect the total magnetic flux of Therefore, there is a problem that the measurement accuracy of the magnet magnetic flux is lowered, and the estimation accuracy of the magnet temperature may be lowered.

本発明のモータの磁石温度の推定方法は、トルク指令値に応じた目標電流を流すことによりモータの回転を制御する。モータの制御方法は、目標電流を変化させ、目標電流を変化させる前の第1状態、及び、目標電流を変化させた後の第2状態のそれぞれにおいてd軸磁束鎖交数を求め、第1状態、及び、第2状態におけるd軸磁束鎖交数に基づいて、回転子が備える永久磁石の磁石磁束鎖交数を算出し、算出した磁石磁束鎖交数に応じて、永久磁石の磁石温度を推定する。 A method for estimating the magnet temperature of a motor according to the present invention controls rotation of a motor by applying a target current according to a torque command value. A method of controlling a motor includes changing a target current, obtaining a d-axis magnetic flux linkage number in each of a first state before changing the target current and a second state after changing the target current, and obtaining a first Based on the number of d-axis magnetic flux linkages in the state and the second state, the number of magnetic flux linkages of the permanent magnets provided in the rotor is calculated, and the magnet temperature of the permanent magnets is calculated according to the calculated number of magnetic flux linkages. to estimate

本発明のモータの磁石温度の推定方法によれば、磁束センサを用いることなく永久磁石の磁石磁束鎖交数を算出し、この磁石磁束鎖交数に応じて永久磁石の推定温度を求めることで、永久磁石の推定温度の精度の向上を図る。 According to the method for estimating the magnet temperature of the motor of the present invention, the magnet flux linkage number of the permanent magnet is calculated without using a magnetic flux sensor, and the estimated temperature of the permanent magnet is obtained according to the magnet flux linkage number. , to improve the accuracy of the estimated temperature of the permanent magnet.

図1は、第1実施形態のモータシステムの概略構成図である。FIG. 1 is a schematic configuration diagram of the motor system of the first embodiment. 図2は、磁石磁束鎖交数算出部の詳細な構成を示すブロック図である。FIG. 2 is a block diagram showing the detailed configuration of the magnet flux linkage number calculator. 図3は、磁束鎖交数推定部の詳細な構成を示すブロック図である。FIG. 3 is a block diagram showing the detailed configuration of the magnetic flux linkage number estimator. 図4は、制御対象へと流れる電流を示す図である。FIG. 4 is a diagram showing the current flowing to the controlled object. 図5は、第2実施形態のモータへと流れる電流を示す図である。FIG. 5 is a diagram showing the current flowing to the motor of the second embodiment. 図6は、第3実施形態のモータへと流れる電流を示す図である。FIG. 6 is a diagram showing current flowing to the motor of the third embodiment. 図7は、磁石磁束鎖交数算出部の詳細な構成を示すブロック図である。FIG. 7 is a block diagram showing the detailed configuration of the magnet flux linkage number calculator. 図8は、第4実施形態のモータシステムを構成するニューラルネットワークの構成を示す図である。FIG. 8 is a diagram showing the configuration of a neural network that constitutes the motor system of the fourth embodiment. 図9は、第5実施形態のモータシステムの一部のトルク誤差補償部に関する構成を示すブロック図である。FIG. 9 is a block diagram showing the configuration of a part of the torque error compensator of the motor system of the fifth embodiment.

本発明の実施形態における磁石温度の推定装置を有するモータシステムについて説明する。 A motor system having a magnet temperature estimation device according to an embodiment of the present invention will be described.

(第1実施形態)
図1は、第1実施形態の磁石温度の推定装置を有するモータシステムの概略構成図である。
(First embodiment)
FIG. 1 is a schematic configuration diagram of a motor system having a magnet temperature estimation device according to the first embodiment.

モータシステム10には、制御対象であるモータ1が含まれている。そのため、モータ1以外の点線で示されるA部が、モータ制御装置に相当する。なお、モータ制御装置は、モータ1のロータが備える永久磁石の磁石温度を推定する。さらに、モータ制御装置は、その推定した磁石温度を用いてモータ1の回転動作を制御する。 A motor system 10 includes a motor 1 to be controlled. Therefore, the portion A indicated by the dotted line other than the motor 1 corresponds to the motor control device. Note that the motor control device estimates the magnet temperature of the permanent magnets provided in the rotor of the motor 1 . Furthermore, the motor control device controls the rotational operation of the motor 1 using the estimated magnet temperature.

モータシステム10においては、電流指令値算出部11は、トルク指令値T*に応じてトルク電流指令値id * ,q *が算出される。なお、電流指令値算出部11には、補正部11Aが並設されており、補正部11Aは推定された磁石温度Tmagに応じた補正値を電流指令値算出部11に出力する。例えば、磁石温度Tmagが高くなると出力されるトルク値が小さくなる傾向にあるので、補正部11Aは、トルク電流指令値id * ,q *が大きくなるように正の補正値を出力する。 In the motor system 10, the current command value calculator 11 calculates torque current command values id * and iq * according to the torque command value T * . A correction unit 11A is provided in parallel with the current command value calculation unit 11, and the correction unit 11A outputs a correction value corresponding to the estimated magnet temperature Tmag to the current command value calculation unit 11. For example, when the magnet temperature T mag increases, the output torque value tends to decrease, so the correction unit 11A outputs positive correction values so that the torque current command values i d * and i q * increase. .

算出されたトルク電流指令値id * ,q *において、d軸トルク電流指令値id *に対して、微小電流重畳部12から出力される微小電流指令値Δid *が加算器13により重畳される。さらに、重畳されたトルク電流指令値id * ,q *に対して、モータ1からのフィードバック入力である電流値id,qが減算器14により負帰還される。そして、加算器13、減算器14を経たトルク電流指令値id * ,q *が負荷電流制御部15に入力される。 In the calculated torque current command values i d * and i q * , the minute current command value Δi d * output from the minute current superimposing unit 12 is added by the adder 13 to the d-axis torque current command value i d * . superimposed. Furthermore, current values id and iq , which are feedback inputs from the motor 1 , are negatively fed back by a subtractor 14 to the superimposed torque current command values id* and iq * . The torque current command values i d * and i q * that have passed through the adder 13 and the subtractor 14 are input to the load current control section 15 .

負荷電流制御部15は、一般的なPI制御器であって、モータ1における観測値である電流値id,qを用いたフィードバック制御を行う。負荷電流制御部15は、モータ1における電流値id,qがトルク電流指令値id * ,q *に追従するように電圧指令値vd * ,q *を生成して、生成した電圧指令値vd * ,q *をモータ1に出力する。 The load current control unit 15 is a general PI controller, and performs feedback control using current values id and iq , which are observed values in the motor 1 . The load current control unit 15 generates voltage command values vd * and vq* such that the current values id and iq in the motor 1 follow the torque current command values id * and iq * , and generates Then, the voltage command values v d * and v q * are output to the motor 1 .

モータ1は、インバータやセンサを含む構成であってもよい。電圧指令値vd * ,q *の入力に応じてインバータによって生成された三相交流電圧がモータ1のバスバに印加されて、モータ1のロータが回転駆動する。そして、インバータとモータ1との間の配線上に電流センサが設けられており、電流センサにより検出された三相交流電流が相変換されることで、dq軸の電流値id,qを取得することができる。 The motor 1 may be configured to include an inverter and a sensor. A three-phase AC voltage generated by an inverter in response to the input of voltage command values v d * and v q * is applied to the busbar of motor 1 to rotate the rotor of motor 1 . A current sensor is provided on the wiring between the inverter and the motor 1, and the three-phase AC current detected by the current sensor is phase-converted to obtain the current values id and iq of the dq axes. can be obtained.

微小電流重畳部12は、上位装置から磁石温度の推定要求に応じてTrigが入力されると、その期間だけ、直流の微小電流Δidが流れるように、微小電流指令値Δid *を生成する。モータシステム10においては、微小電流指令値Δid *がトルク電流指令値id * ,q *に重畳されることで、後述の磁石温度の推定が行われる。そして、微小電流指令値Δid *が重畳される期間と、重畳されない期間との2つの期間の測定値を用いて、モータ1のロータが備える永久磁石の磁石温度Tmagを推定する。磁石温度Tmagの推定方法の詳細については後に説明する。 When Trig is input from a host device in response to a magnet temperature estimation request, the minute current superimposing unit 12 generates a minute current command value Δid * so that a minute DC current Δid flows only during that period. . In the motor system 10, the magnet temperature, which will be described later, is estimated by superimposing the minute current command value Δid * on the torque current command values id * and iq * . Then, the magnet temperature T mag of the permanent magnet provided in the rotor of the motor 1 is estimated using measured values during two periods, the period during which the minute current command value Δi d * is superimposed and the period during which it is not superimposed. The details of the method of estimating the magnet temperature Tmag will be described later.

磁石磁束鎖交数算出部16には、モータ1への入力となる電圧指令値vd * ,q *、モータ1において検出される電流値id,q、及び、微小電流重畳部12に入力されるTrigに加えて、モータ1のロータの近傍に設けられるレゾルバにより検出される電気位相角θeが入力される。磁石磁束鎖交数算出部16は、これらの入力に基づいてモータ1のロータに設けられる永久磁石に起因する磁石磁束鎖交数ψaを推定する。 The magnet flux linkage calculation unit 16 includes voltage command values v d * and v q * to be input to the motor 1 , current values id and i q detected in the motor 1 , and the minute current superimposition unit 12 In addition to Trig input to , an electrical phase angle θe detected by a resolver provided near the rotor of the motor 1 is input. The magnet flux linkage number calculator 16 estimates the magnet flux linkage number ψ a resulting from the permanent magnets provided in the rotor of the motor 1 based on these inputs.

磁石温度推定部17は、予め記憶しているテーブルを用いて、磁石磁束鎖交数算出部16により推定された磁石磁束鎖交数ψaと、モータ1において観測される電流値id,qに基づいて、磁石温度Tmagを推定する。 The magnet temperature estimation unit 17 uses a table stored in advance to obtain the magnet flux linkage number ψ a estimated by the magnet flux linkage number calculation unit 16 and the current values id, i Based on q , estimate the magnet temperature Tmag .

なお、この推定された磁石温度Tmagは、補正部11Aを経て電流指令値算出部10に再帰的に入力され、磁石温度Tmagに応じてトルク電流指令値id * ,q *を変化させることで、モータ1の回転制御の精度を向上させることができる。また、第5実施形態のように、モータシステム10内にトルク誤差補正部(図9参照)を設けて、磁石温度Tmagの推定過程で求められるモータ1の磁束鎖交数を用いて、モータ1の回転制御の精度の向上を図ってもよい。 The estimated magnet temperature T mag is recursively input to the current command value calculation unit 10 via the correction unit 11A, and the torque current command values id * and i q * are changed according to the magnet temperature T mag. By doing so, the accuracy of the rotation control of the motor 1 can be improved. Further, as in the fifth embodiment, a torque error correction unit (see FIG. 9) is provided in the motor system 10, and the motor The accuracy of rotation control of 1 may be improved.

図2は、磁石磁束鎖交数算出部16の詳細な構成を示すブロック図である。 FIG. 2 is a block diagram showing a detailed configuration of the magnet flux linkage number calculator 16. As shown in FIG.

磁石磁束鎖交数算出部16は、磁束鎖交数推定部21と、磁石磁束鎖交数決定部22とを有する。 The magnet flux linkage number calculation unit 16 has a magnetic flux linkage number estimation unit 21 and a magnet flux linkage number determination unit 22 .

磁束鎖交数推定部21は、電圧指令値vd * ,q *、電流値id,q、及び、ロータの電気位相角θeの入力を受け付けると、磁束鎖交数推定値ψd ^ ,ψq ^を推定し、磁石磁束鎖交数決定部22に出力する。なお、磁束鎖交数推定部21における詳細な処理は、後に図3を用いて説明する。また、磁束鎖交数推定部21においてはdq軸磁束鎖交数推定値ψd ^ ,ψq ^が算出されるが、後段の磁石磁束鎖交数決定部22における処理には、d軸磁束鎖交数推定値ψd ^のみが用いられる。 Upon receiving inputs of the voltage command values v d * , v q * , the current values id, i q , and the rotor electrical phase angle θe, the flux linkage number estimator 21 calculates the flux linkage number estimation value ψ d ^ , ψ q ^ are estimated and output to the magnet flux linkage number determination unit 22 . Detailed processing in the magnetic flux linkage number estimating unit 21 will be described later with reference to FIG. In the magnetic flux linkage number estimation unit 21, the dq-axis magnetic flux linkage number estimated values ψ d ^ and ψ q ^ are calculated. Only the linkage number estimate ψ d ^ is used.

磁石磁束鎖交数決定部22は、磁束鎖交数推定部21により推定されたdq軸磁束鎖交数推定値ψd ^ ,ψq ^のうちのd軸磁束鎖交数推定値ψd ^、電流値id,q、及び、Trigの入力を受け付ける。図1に示された微小電流重畳部12は、Trigが入力されている期間だけ微小電流指令値Δid *を重畳するため、磁石磁束鎖交数決定部22には、Trigが入力されている期間と、入力されていない期間との2つの期間のd軸磁束鎖交数推定値ψd ^が入力される。そして、磁石磁束鎖交数決定部22は、それらの2つの期間のd軸磁束鎖交数推定値ψd ^を用いて、後述の(4)式に従って、磁石温度Tmagを推定する。 The magnet flux linkage number determination unit 22 determines the d - axis flux linkage number estimated value ψ d ^ among the dq-axis flux linkage number estimated values ψ d ^ and ψ q ^ estimated by the magnetic flux linkage number estimation unit 21 . , current values i d , i q , and Trig. Since the minute current superimposing unit 12 shown in FIG. 1 superimposes the minute current command value Δi d * only during the period when Trig is input, Trig is input to the magnet flux linkage number determining unit 22. The d-axis flux linkage number estimated values ψ d ^ for two periods, one period and the other that are not input, are input. Then, the magnet flux linkage number determination unit 22 uses the d-axis flux linkage number estimated values ψ d ^ for these two periods to estimate the magnet temperature T mag according to Equation (4) described later.

図3は、磁束鎖交数推定部21の詳細な構成を示すブロック図である。dq軸はロータと同期した回転座標系であり、αβ軸は静止座標系であるものとする。 FIG. 3 is a block diagram showing the detailed configuration of the magnetic flux linkage number estimator 21. As shown in FIG. It is assumed that the dq axes are a rotating coordinate system synchronized with the rotor, and the αβ axes are a stationary coordinate system.

磁束鎖交数推定部21は、抵抗値算出部31、dq/αβ変換部32、及び、減算器301からなる電圧モデル部と、インダクタンス算出部33、加算器302、dq/αβ変換部34、減算器303、及び、PI制御器35からなる電流モデル部とを有する。そして、電圧モデル部及び電流モデル部からの出力は加算器304において加算され、加算結果が積分器36において時間積分されることで、αβ軸磁束鎖交数推定値ψα ^ ,ψβ ^が算出される。最終的に、αβ/dq変換部37を経て、dq軸磁束鎖交数推定値ψd ^ ,ψq ^が出力される。なお、dq/αβ変換部32、dq/αβ変換部34、及び、αβ/dq変換部37における変換処理においては、レゾルバにより検出されるロータの電気位相角θeが用いられる。 The magnetic flux linkage estimation unit 21 includes a voltage model unit including a resistance value calculation unit 31, a dq/αβ conversion unit 32, and a subtractor 301, an inductance calculation unit 33, an adder 302, a dq/αβ conversion unit 34, It has a subtractor 303 and a current model section consisting of a PI controller 35 . Then, the outputs from the voltage model section and the current model section are added in the adder 304, and the addition result is time-integrated in the integrator 36, so that the αβ-axis flux linkage number estimated values ψ α ^ and ψ β ^ are Calculated. Finally, the dq-axis magnetic flux linkage number estimated values ψ d ̂ and ψ q ̂ are output through the αβ/dq conversion unit 37 . Note that the rotor electrical phase angle θe detected by the resolver is used in the conversion processing in the dq/αβ conversion unit 32, the dq/αβ conversion unit 34, and the αβ/dq conversion unit 37.

まず、電圧モデル部について説明する。 First, the voltage model section will be described.

抵抗値算出部31は、入力される電流値id,qに対して抵抗値Rd,qを乗ずることで、モータ1の巻線抵抗に生じる電圧値を算出する。この値は、後段に設けられるdq/αβ変換部32によってαβ座標の値に変換される。そして、減算器301において、電圧指令値vα * ,β *から、dq/αβ変換部32から出力される巻線抵抗の電圧値を減算することでΔvα,Δvβを算出する。 The resistance value calculator 31 multiplies the input current values id and iq by the resistance values Rd and Rq to calculate the voltage values generated in the winding resistances of the motor 1 . This value is converted into an αβ coordinate value by a dq/αβ conversion unit 32 provided at a later stage. Then, subtractor 301 calculates Δv α and Δv β by subtracting the voltage value of the winding resistance output from dq/αβ conversion unit 32 from voltage command values v α * and v β * .

そして、電圧モデル部から出力されるΔvα,Δvβは、加算器304において、後述のように電流モデル部からの出力値と加算される。積分器36は、加算器304における加算結果に対して積分を行うことで、αβ軸磁束鎖交数推定値ψα ^ ,ψβ ^を算出する。αβ/dq変換部37においては、入力されるαβ軸磁束鎖交数推定値ψα ^ ,ψβ ^に対してαβ/dq変換をすることで、dq軸磁束鎖交数推定値ψd ^ ,ψq ^を算出する。 Δv α and Δv β output from the voltage model section are added to an output value from the current model section in adder 304 as will be described later. The integrator 36 integrates the addition result of the adder 304 to calculate the αβ-axis magnetic flux linkage number estimated values ψ α ̂ and ψ β ̂ . The αβ/dq conversion unit 37 performs αβ/dq conversion on the input αβ-axis magnetic flux linkage estimated values ψ α ^ and ψ β ^ to obtain the dq-axis magnetic flux linkage estimated values ψ d ^ , ψ q ^ .

ここで、モータ1においては、磁束鎖交数の時間変化量に応じて誘起電圧が発生しており、この誘起電圧は電圧指令値と巻線抵抗の電圧値との差Δvα,Δvβに相当する。そのため、誘起電圧に相当するΔvα,Δvβを時間積分することで、モータ1における鎖交磁束数を求めることができる。そこで、Δvα,Δvβを後段の積分器36において時間積分することで、モータ1における磁束鎖交数として、αβ軸磁束鎖交数推定値ψα ^ ,ψβ ^が求められる。 Here, in the motor 1, an induced voltage is generated according to the amount of change in the magnetic flux linkage over time . Equivalent to. Therefore, the number of interlinkage magnetic fluxes in the motor 1 can be obtained by time-integrating Δv α and Δv β corresponding to the induced voltage. Therefore, by time-integrating Δv α and Δv β in the post-stage integrator 36, αβ-axis magnetic flux linkage estimated values ψ α ̂ and ψ β ̂ are obtained as the magnetic flux linkages in the motor 1 .

しかしながら、積分処理に起因して不定数である積分乗数が生じてしまう。そこで、積分乗数により生じる誤差を抑制するために、インダクタンス算出部33~PI制御器35からなる電流モデル部が構成されている。以下では、電流モデル部について説明する。 However, due to the integration process, an indefinite integral multiplier is generated. Therefore, in order to suppress the error caused by the integral multiplier, the current model section including the inductance calculator 33 to the PI controller 35 is configured. The current model section will be described below.

インダクタンス算出部33は、電流値id,qに対してモータ1のインダクタンス値Ld,qを乗ずることで、モータ1において電流の負荷に応じて発生する磁束鎖交数を算出する。この処理は、一般的な電流、インダクタンス、及び、磁束鎖交数の関係に基づいている。なお、インダクタンス値Ld,qは固定値であるものとする。 The inductance calculator 33 multiplies the current values id and iq by the inductance values Ld and Lq of the motor 1 to calculate the number of magnetic flux linkages generated in the motor 1 according to the current load. This process is based on common current, inductance and flux linkage relationships. It is assumed that the inductance values L d and L q are fixed values.

加算器302は、インダクタンス算出部33において算出された電流に起因する磁束鎖交数と、磁石磁束鎖交数固定値ψa'とを加算することで、dq軸磁束鎖交数推定値(電流モデル)ψd_i ^ ,ψq_i ^を算出する。なお、磁石磁束鎖交数固定値ψa'は、モータ1において電流が印加されない無負荷状態における永久磁石に起因する磁石磁束鎖交数であって、予め磁束鎖交数推定部21に記憶されている。このように、磁石磁束鎖交数固定値ψa'を予め加算しておくことで、初期状態がない駆動直後などの状態におる誤差を抑制できる。 The adder 302 adds the magnetic flux linkage number resulting from the current calculated in the inductance calculation unit 33 and the magnet flux linkage fixed value ψ a ' to obtain the dq-axis magnetic flux linkage estimated value (current Model) Calculate ψ d_i ^ , ψ q_i ^ . The magnet flux linkage number fixed value ψ a ' is the magnet flux linkage number caused by the permanent magnet in the no-load state in which no current is applied to the motor 1, and is stored in advance in the flux linkage number estimator 21. ing. In this way, by adding the magnet flux linkage fixed value ψ a ' in advance, it is possible to suppress an error in a state such as immediately after driving when there is no initial state.

dq/αβ変換部34においては、dq軸磁束鎖交数推定値(電流モデル)ψd_i ^ ,ψq_i ^に対してdq/αβの変換をすることで、αβ軸磁束鎖交数推定値(電流モデル)ψα_i ^ ,ψβ_i ^を算出する。 In the dq/αβ conversion unit 34 , the estimated αβ - axis magnetic flux linkage number ( current model) ψ α_i ^ , ψ β_i ^ are calculated.

減算器303においては、αβ軸磁束鎖交数推定値(電流モデル)ψα_i ^ ,ψβ_i ^から、後段の積分器36から出力されるαβ軸磁束鎖交数推定値ψα ^ ,ψβ ^を減じることで、Δψα_i,Δψβ_iを算出する。 In the subtractor 303, the αβ-axis magnetic flux linkage estimated values (current model) ψ α_i ^ , ψ β_i ^ are converted to the αβ-axis magnetic flux linkage estimated values ψ α ^ , ψ β output from the integrator 36 in the subsequent stage. Δψ α_i and Δψ β_i are calculated by subtracting ^ .

そして、PI制御器35においては、Δψα_i,Δψβ_iに対してPI制御を行うことにより、αβ軸磁束鎖交数推定値(電流モデル)ψα_i ^ ,ψβ_i ^とαβ軸磁束鎖交数推定値ψα ^ ,ψβ ^との偏差が小さくなるように制御される。 Then, the PI controller 35 performs PI control on Δψ α_i and Δψ β_i to obtain αβ-axis magnetic flux linkage estimated values (current models) ψ α_i ^ , ψ β_i ^ and αβ-axis magnetic flux linkage Control is performed so that deviations from the estimated values ψ α ^ and ψ β ^ are small.

上述のように電圧モデル部からの出力を積分器36により積分する際に、積分定数が生じて系が不安定となるおそれがある。しかしながら、電流モデル部において、PI制御器35によってフィードバック系が組まれるとともに、減算器301においては初期値として磁石磁束鎖交数固定値ψa’が入力される。このような構成となることで、初期状態を含む任意の状態において、積分定数が不定となることに起因するαβ軸磁束鎖交数推定値ψα ^ ,ψβ ^の精度悪化を抑制することができる。 When the output from the voltage model section is integrated by the integrator 36 as described above, an integration constant may be generated and the system may become unstable. However, in the current model section, a feedback system is formed by the PI controller 35, and the magnet flux linkage number fixed value ψ a ' is input to the subtractor 301 as an initial value. With such a configuration, in any state including the initial state, it is possible to suppress deterioration in the accuracy of the αβ-axis magnetic flux linkage number estimation values ψ α ^ and ψ β ^ due to the constant of integration becoming indefinite. can be done.

ここで、電圧モデル部、及び、電流モデル部それぞれの、dq軸磁束鎖交数推定値ψd ^ ,ψq ^の推定に対する寄与度は、PI制御器35において定めることができる。例えば、PI制御器35におけるゲインが低下するように設計することで、電流モデルに比べて電圧モデルがαβ軸磁束鎖交数推定値ψα ^ ,ψβ ^の推定に強く影響することになる。また、モータ1の回転数が大きい場合には、電流モデル部に設けられるPI制御器35のゲインが低下することが知られており、このような高回転時のゲイン低下が考慮されてPI制御器35が設計される。 Here, the PI controller 35 can determine the degree of contribution of each of the voltage model section and the current model section to the estimation of the dq-axis magnetic flux linkage estimated values ψ d ^ and ψ q ̂ . For example, by designing the gain in the PI controller 35 to decrease, the voltage model has a stronger influence on the estimation of the αβ-axis magnetic flux linkage number estimated values ψ α ^ and ψ β ^ than the current model. . It is also known that when the number of revolutions of the motor 1 is large, the gain of the PI controller 35 provided in the current model section decreases. A vessel 35 is designed.

図4は、モータ1へと流れる電流を示す図である。この図において、x軸は時間を示し、y軸はd軸電流値idを示す。 FIG. 4 is a diagram showing the current flowing to the motor 1. As shown in FIG. In this figure, the x-axis indicates time and the y-axis indicates the d-axis current value id .

期間1において微小電流重畳部12は何も出力しておらず、期間2においてはTrigに応じて微小電流重畳部12が微小電流指令値Δid *を出力する。そのため、期間1においてidであった電流値は、期間2においては微小電流重畳部12から微小電流指令値Δid *が重畳されることにより、d軸電流値id’となる。なお、この期間2は、Trigが微小電流重畳部12に入力される期間である。なお、期間1は第1状態の一例であり、期間2は第2状態の一例である。 During the period 1, the minute current superimposing unit 12 does not output anything, and during the period 2, the minute current superimposing unit 12 outputs the minute current command value Δi d * according to Trig. Therefore, the current value i d in the period 1 becomes the d-axis current value i d ′ in the period 2 by superimposing the minute current command value Δi d * from the minute current superimposing unit 12 . Note that this period 2 is a period during which Trig is input to the minute current superimposing unit 12 . Note that the period 1 is an example of the first state, and the period 2 is an example of the second state.

次に、磁石磁束鎖交数決定部22における磁石磁束鎖交数の算出方法について説明する。 Next, a method for calculating the magnet flux linkage number in the magnet flux linkage number determination unit 22 will be described.

上述のように、磁束鎖交数推定部21は、電流値id,qから磁束鎖交数ψd,ψqを推定する。そして、電流値id,qが流れる間において、モータ1において生じるd軸磁束鎖交数ψdは、磁石磁束鎖交数ψaと、電流に起因する磁束鎖交数との和となる。なお、電流に起因する磁束鎖交数は、d軸インダクタンスLd及びd軸電流値idの積(Ld・id)により定められる。 As described above, the magnetic flux linkage number estimator 21 estimates the magnetic flux linkage numbers ψ d and ψ q from the current values id and i q . The d-axis magnetic flux linkage number ψ d generated in the motor 1 while the current values id and i q flow is the sum of the magnet magnetic flux linkage number ψ a and the magnetic flux linkage number caused by the current. . Note that the number of magnetic flux linkages caused by the current is determined by the product (L d · id ) of the d-axis inductance L d and the d-axis current value id .

そのため、期間1におけるd軸磁束鎖交数をψd、d軸電流をidとし、期間2におけるd軸磁束鎖交数をψd’、d軸電流をid’とすると、期間1と期間2とにおいて、それぞれ次式が成立する。 Therefore, if the number of d-axis magnetic flux linkages in period 1 is ψ d and the d-axis current is id , and the number of d-axis magnetic flux linkages in period 2 is ψ d ′ and the d-axis current is id ′, period 1 and In period 2, the following equations are respectively established.

Figure 0007324630000001
Figure 0007324630000001

Figure 0007324630000002
Figure 0007324630000002

これらの式(1)、(2)のそれぞれについて、インダクタンスLが等しいと仮定すると、(1)式からインダクタンスLは以下のように示すことができる。 Assuming that the inductance L is equal for each of these equations (1) and (2), the inductance L can be expressed as follows from equation (1).

Figure 0007324630000003
Figure 0007324630000003

(3)式を(2)式に代入すれば、次の式を導出できる。 By substituting equation (3) into equation (2), the following equation can be derived.

Figure 0007324630000004
Figure 0007324630000004

磁石磁束鎖交数決定部22は、期間1、2の2つの期間における測定値を用いてd軸磁束鎖交数ψd、ψd’を算出し、(4)式に基づいて、d軸磁束鎖交数ψd、ψd’とd軸電流値id、id’とを用いて磁石磁束鎖交数ψaを推定することができる。 The magnet flux linkage number determination unit 22 calculates the d-axis magnetic flux linkage numbers ψ d and ψ d ′ using the measured values in the two periods 1 and 2, and based on the equation (4), the d-axis The magnetic flux linkage number ψ a can be estimated using the magnetic flux linkage numbers ψ d , ψ d ′ and the d-axis current values id , id ′.

磁石温度推定部17は、予め、d軸電流値id、q軸電流値iq、及び、磁石磁束鎖交数ψaと、磁石温度Tmagとを対応つけたテーブルを記憶しており、入力値とテーブルとを用いて、磁石温度Tmagを推定する。なお、磁石温度推定部17は、1つの電流動作点だけでなく多数の電流動作点に対し個々にテーブルを記憶することで、磁石温度Tmagの推定精度を向上させることができる。 The magnet temperature estimator 17 stores in advance a table in which the d-axis current value i d , the q-axis current value i q , and the magnet flux linkage number ψ a are associated with the magnet temperature T mag , Estimate the magnet temperature T mag using the input values and the table. It should be noted that the magnet temperature estimator 17 can improve the estimation accuracy of the magnet temperature T mag by individually storing tables not only for one current operating point but also for many current operating points.

なお、テーブル値は、d軸電流値id、q軸電流値iq、磁石磁束鎖交数ψa、及び、磁石温度Tmagの4つのパラメータを持つテーブルであるため、全てのパラメータを精度よく対応させたテーブルを記憶するのは難しい。そこで、入力値がテーブルに存在しない場合には、線形補完などの方法を用いてテーブルを補完する。このようにすることで、予めテーブルの作成に要する時間を短くすることができるだけでなく、磁石温度推定部17が記憶するテーブルの容量を小さくすることができる。なお、第4実施形態において、テーブルを用いずにニューラルネットワークを用いて磁石温度Tmagを推定する説明する。 Note that the table value is a table having four parameters of the d-axis current value i d , the q-axis current value i q , the magnet flux linkage number ψ a , and the magnet temperature T mag . It is difficult to remember well-matched tables. Therefore, when the input value does not exist in the table, the table is interpolated using a method such as linear interpolation. By doing so, it is possible not only to shorten the time required to prepare the table in advance, but also to reduce the capacity of the table stored in the magnet temperature estimating section 17 . In the fourth embodiment, a neural network is used to estimate the magnet temperature T mag without using a table.

なお、本実施形態においては、微小電流重畳部12を用いてd軸の微小電流を重畳させる例について説明したが、これに限らない。d軸成分を服務微小電流が重畳されてもよい。例えば、アクセル操作などに応じて所望のタイミングで電流の大きさが変化する場合には、そのタイミングに応じてd軸磁束交鎖数ψdを推定することで、磁石磁束鎖交数ψaを求めることができる。 In the present embodiment, an example of superimposing a minute current on the d-axis using the minute current superimposing unit 12 has been described, but the present invention is not limited to this. A minute current serving the d-axis component may be superimposed. For example, when the magnitude of the current changes at a desired timing according to the accelerator operation, etc., by estimating the d-axis magnetic flux linkage number ψd according to the timing, the magnet flux linkage number ψa can be calculated as can ask.

第1実施形態によれば、以下の効果を得ることができる。 According to the first embodiment, the following effects can be obtained.

第1実施形態のモータ1の磁石温度の推定方法によれば、微小電流重畳部12は、Trigの入力に応じて、目標トルクに応じたトルク電流指令値id * ,q *に対して、トルク電流指令値id * ,q *よりも小さな微小電流指令値Δid *を重畳することでトルク電流指令値id * ,q *を変化させる。そして、磁束鎖交数推定部21は、微小電流指令値Δid *が重畳される期間1(第1状態)と、重畳されない期間2(第2状態)とのそれぞれにおいてd軸磁束鎖交数推定値ψd ^を推定する。磁石磁束鎖交数決定部22は、期間1と期間2との場合におけるd軸磁束鎖交数推定値ψd ^に基づいて、回転子が備える永久磁石の磁石磁束鎖交数ψaを算出する。そして、磁石温度推定部17は、磁石磁束鎖交数ψaに基づいて永久磁石の磁石温度Tmagを推定する。 According to the method for estimating the magnet temperature of the motor 1 of the first embodiment, the minute current superimposing unit 12 generates torque current command values i d * and i q * corresponding to the target torque according to the input of Trig. , torque current command values id * and iq * , which are smaller than the torque current command values id * and iq * . Then, the magnetic flux linkage number estimator 21 calculates the d - axis magnetic flux linkage number Estimate the estimated value ψ d ^ . The magnet flux linkage number determination unit 22 calculates the magnet flux linkage number ψ a of the permanent magnets included in the rotor based on the d-axis flux linkage number estimated value ψ d ^ in the case of period 1 and period 2. do. Then, the magnet temperature estimator 17 estimates the magnet temperature T mag of the permanent magnet based on the magnet flux linkage number ψ a .

また、磁束の変化に応じて生じる誘起電圧に相当するΔvα,Δvβを求め、その誘起電圧を積分器36において積分することで、磁束鎖交数推定値ψd ^ ,ψq ^を推定する。これにより、磁束センサを用いることなく磁石磁束鎖交数ψaを求めることができる。そのため、磁束センサが磁石磁束以外の負荷磁束を測定するおそれがなくなり、磁石鎖交磁束数を精度よく求めることができ、その結果、磁石温度Tmagの推定精度を向上させることができる。その結果、高い精度で推定された磁石温度Tmagを用いて目標電流を制御することで、モータ1の回転制御の精度の向上を図ることができる。 Also, Δv α and Δv β corresponding to the induced voltages generated according to the change in the magnetic flux are obtained, and the induced voltages are integrated in the integrator 36 to estimate the flux linkage number estimated values ψ d ^ and ψ q ̂ . do. As a result, the magnet flux linkage number ψ a can be obtained without using a magnetic flux sensor. Therefore, the magnetic flux sensor is not likely to measure the load magnetic flux other than the magnet magnetic flux, and the magnet interlinkage magnetic flux number can be obtained with high accuracy. As a result, the estimation accuracy of the magnet temperature Tmag can be improved. As a result, by controlling the target current using the magnet temperature Tmag estimated with high accuracy, the accuracy of the rotation control of the motor 1 can be improved.

第1実施形態のモータ1の磁石温度の推定方法によれば、微小電流指令値Δid *の重畳の有無に応じた期間1、2において、異なるd軸磁束鎖交数推定値ψd ^が求められる。ここで、d軸磁束鎖交数推定値ψd ^の一部である磁石磁束鎖交数ψaは、微小電流指令値Δid *の重畳の有無では変化しないものと仮定する。さらに、d軸磁束鎖交数推定値ψd ^のその他の部分は、d軸電流値idとd軸インダクタンスLdとの積に相当するが、d軸インダクタンスLdは微小電流指令値Δid *の重畳の有無では変化しないものと仮定する。 According to the method for estimating the magnet temperature of the motor 1 of the first embodiment, different d - axis magnetic flux linkage number estimated values ψ d ̂ Desired. Here, it is assumed that the magnet flux linkage number ψ a , which is a part of the d-axis magnetic flux linkage number estimated value ψ d ^ , does not change depending on whether or not the minute current command value Δi d * is superimposed. Furthermore, the other part of the d-axis flux linkage estimated value ψ d ^ corresponds to the product of the d-axis current value id and the d-axis inductance L d , but the d-axis inductance L d is the minute current command value Δi It is assumed that d * does not change with or without superimposition.

そこで、まず、微小電流指令値Δid *を重畳しない期間1における場合において、d軸磁束鎖交数推定値ψd ^とd軸電流値idとを用いてd軸インダクタンスLdを算出しておく。そして、算出したd軸インダクタンスLdを微小電流指令値Δid *を重畳する期間2における場合に代入することで磁石磁束鎖交数ψaを求めることができる。詳細には、式(4)を用いることにより、磁石磁束鎖交数ψaを求めることができる。 Therefore, first, in the case of period 1 in which the minute current command value Δi d * is not superimposed, the d-axis inductance L d is calculated using the d-axis flux linkage estimated value ψ d ^ and the d-axis current value id. Keep Then, by substituting the calculated d-axis inductance L d for the case of period 2 in which the minute current command value Δi d * is superimposed, the magnet flux linkage number ψ a can be obtained. Specifically, the magnet flux linkage number ψ a can be obtained by using the equation (4).

第1実施形態のモータ1の磁石温度の推定方法によれば、磁石温度推定部17は、予め、d軸電流値id、q軸電流値iq、及び、磁石磁束鎖交数ψaと、磁石温度Tmagとの対応関係を示すテーブルを記憶しており、入力されるこれらのパラメータと、それらの対応関係を示すテーブルとを用いて、磁石温度Tmagを推定する。このようにテーブルを用いることで、磁石温度推定部17の処理負荷を高めることなく、テーブル検索のみで磁石温度Tmagを推定することができる。 According to the method for estimating the magnet temperature of the motor 1 of the first embodiment, the magnet temperature estimator 17 preliminarily calculates the d-axis current value i d , the q-axis current value i q , and the magnet flux linkage number ψ a and , magnet temperature T mag , and the magnet temperature T mag is estimated using these input parameters and the table showing their correspondence. By using the table in this way, the magnet temperature T mag can be estimated only by searching the table without increasing the processing load of the magnet temperature estimator 17 .

また、テーブルに記憶されるサンプル数が少ない場合において、入力されるパラメータがテーブル含まれない場合には、テーブルに存在しないパラメータをテーブルにおいて補完し、補完されたテーブルを用いて磁石温度Tmagを推定する。このようにすることで、磁石温度推定部17の処理負荷は若干高くなるが、テーブルを小さくすることがでるのでメモリ使用量の低減を図ることができる。 In addition, when the number of samples stored in the table is small, if the input parameter is not included in the table, the parameter that does not exist in the table is complemented in the table, and the magnet temperature T mag is calculated using the complemented table. presume. By doing so, the processing load of the magnet temperature estimating unit 17 is slightly increased, but the table can be made smaller, so that the amount of memory used can be reduced.

(第2実施形態)
第1実施形態においては、微小電流重畳部12は、期間1において微小電流を出力せず、期間2において正の微小電流を出力したがこれに限らない。第2実施形態においては、微小電流重畳部12は、微小電流の無出力、正の微小電流の重畳、及び、負の微小電流の重畳を繰り返し行ってもよい。
(Second embodiment)
In the first embodiment, the micro-current superimposing unit 12 does not output a micro-current in period 1 and outputs a positive micro-current in period 2, but it is not limited to this. In the second embodiment, the micro-current superimposing unit 12 may repeatedly perform non-output of micro-current, superimposition of positive micro-current, and superimposition of negative micro-current.

図5は、モータ1に流れる電流を示す図である。期間1において微小電流重畳部12は何も出力しておらず電流値はidである。期間2において微小電流重畳部12は正の微小電流が重畳されるように微小電流指令値Δid *を出力しているので、モータ1に流れる電流値はid2となる。期間3において微小電流重畳部12は負の微小電流が重畳されるように指令値-Δid *を出力しているので、電流値はid3となる。 FIG. 5 is a diagram showing the current flowing through the motor 1. As shown in FIG. In the period 1, the minute current superimposing unit 12 does not output anything and the current value is id . In period 2, the minute current superimposing unit 12 outputs the minute current command value Δi d * so that a positive minute current is superimposed, so the current value flowing through the motor 1 becomes id2 . In the period 3, the minute current superimposing unit 12 outputs the command value -Δi d * so that the minute negative current is superimposed, so the current value becomes i d3 .

ここで、期間1及び2と、期間1及び3とにおいては、それぞれに、式(3)の関係が成立するため、次の2つの式が成立する。 Here, in the periods 1 and 2 and in the periods 1 and 3, since the relationship of formula (3) holds respectively, the following two formulas hold.

Figure 0007324630000005
Figure 0007324630000005

Figure 0007324630000006
Figure 0007324630000006

ただし、ψa1は、期間1及び2を用いて算出した磁石磁束鎖交数、ψa2は、期間1及び3を用いて算出した磁石磁束鎖交数である。さらに、ψd1、id1は、それぞれ、期間1におけるd軸磁束鎖交数、d軸電流であり、ψd2、id2は、それぞれ、期間2におけるd軸磁束鎖交数、d軸電流であり、ψd3、id3は、それぞれ、期間3におけるd軸磁束鎖交数、d軸電流である。 However, ψ a1 is the magnet flux linkage number calculated using the periods 1 and 2, and ψ a2 is the magnet flux linkage number calculated using the periods 1 and 3. Furthermore, ψ d1 and i d1 are the d-axis flux linkage number and d-axis current in period 1, respectively, and ψ d2 and i d2 are the d-axis flux linkage number and d-axis current in period 2, respectively. , and ψ d3 and i d3 are the d-axis flux linkage number and the d-axis current in period 3, respectively.

そして、式(5)、(6)の平均をとることで、磁石磁束鎖交数ψaは次式のように求められる。 Then, by averaging the equations (5) and (6), the magnet flux linkage number ψ a can be obtained as in the following equation.

Figure 0007324630000007
Figure 0007324630000007

このような(7)式を用いることで、微小電流重畳部12の構成を大きく変えることなく、正及び負の微小電流が重畳される状態を用いて磁石磁束鎖交数ψaを算出できるので、推定精度を向上させることができる。なお、本実施形態においては、同じ大きさで符号の異なる微小電流を重畳させたが、これに限らない。異なる大きさの3以上の微小電流を重畳させ、3以上の磁石磁束鎖交数を算出して、これらの平均を求めてもよい。 By using such equation (7), the magnet flux linkage number ψa can be calculated using the state in which the positive and negative minute currents are superimposed without significantly changing the configuration of the minute current superimposing unit 12. , can improve the estimation accuracy. In this embodiment, minute currents of the same magnitude but different signs are superimposed, but the present invention is not limited to this. Three or more minute currents of different magnitudes may be superimposed, three or more magnet flux linkages may be calculated, and an average of these may be obtained.

第2実施形態によれば以下の効果を得ることができる。 According to the second embodiment, the following effects can be obtained.

第2実施形態のモータ1の磁石温度の推定方法によれば、微小電流重畳部12は、期間1においては微小電流を重畳せず、期間2においては正の微小電流指令値Δid *を重畳し、期間3において負の微小電流指令値Δid *を重畳する。そして、期間1、2のd軸磁束鎖交数ψdを用いて算出する磁石磁束鎖交数と、期間1、3のd軸磁束鎖交数ψdを用いて算出する磁石磁束鎖交数との平均を、最終的な磁石磁束鎖交数ψaとする。このようにサンプル数を増やすことにより、より高い精度で磁石磁束鎖交数ψaを推定できるので、磁石温度Tmagの推定精度の向上を図ることができる。 According to the method for estimating the magnet temperature of the motor 1 of the second embodiment, the minute current superimposing unit 12 does not superimpose a minute current in period 1, and superimposes a positive minute current command value Δi d * in period 2. , and in period 3, a negative minute current command value Δi d * is superimposed. Then, the magnet flux linkage number calculated using the d-axis flux linkage number ψ d in periods 1 and 2, and the magnet flux linkage number calculated using the d-axis flux linkage number ψ d in periods 1 and 3 and the final magnet flux linkage number ψ a . By increasing the number of samples in this way, the magnet flux linkage ψ a can be estimated with higher accuracy, so that the estimation accuracy of the magnet temperature T mag can be improved.

(第3実施形態)
第1、2実施形態においては、所定の直流の微小電流を重畳させたがこれに限らない。第3実施形態においては、交流の微小電流を重畳させる例について説明する。なお、この方法によれば、モータ1においては直流成分で制御されており、交流成分は磁石温度の推定にのみ用いられるものとして区別することができるため、磁石温度Tmagの推定精度のさらなる向上を図ることができる。
(Third embodiment)
In the first and second embodiments, a predetermined minute DC current is superimposed, but the present invention is not limited to this. In the third embodiment, an example of superimposing a very small alternating current will be described. According to this method, since the motor 1 is controlled by the DC component and the AC component can be distinguished as being used only for estimating the magnet temperature, the estimation accuracy of the magnet temperature Tmag is further improved. can be achieved.

図6は、第3実施形態のモータシステム10によってモータ1に流れる電流を示す図である。この図に示されるように、微小電流重畳部12は、Trigに応じて、期間2において周波数fHF、振幅iHFの交流電流を出力しており、交流電流がd軸電流値idに重畳される。 FIG. 6 is a diagram showing the current flowing through the motor 1 by the motor system 10 of the third embodiment. As shown in this figure, the minute current superimposing unit 12 outputs an alternating current of frequency f HF and amplitude i HF in period 2 according to Trig, and the alternating current is superimposed on the d -axis current value id. be done.

図7は、第3実施形態のモータシステム10における磁石磁束鎖交数算出部16の構成を示す図である。本実施形態の磁石磁束鎖交数算出部16は、図2に示された第1実施形態の磁石磁束鎖交数算出部16と比較すると、インダクタンス推定部71が追加されている。 FIG. 7 is a diagram showing the configuration of the magnet flux linkage number calculator 16 in the motor system 10 of the third embodiment. The magnet flux linkage number calculator 16 of the present embodiment has an additional inductance estimator 71 when compared with the magnet flux linkage number calculator 16 of the first embodiment shown in FIG.

ここで、微小電流重畳部12から交流電流が入力されている期間2においては、モータ1における電圧ベクトルには、インダクタンスLに起因して生じる電圧成分であって、重畳される電流ベクトルに対して90度の位相が進む成分vHF_imgが含まれる。このインダクタンスが寄与する成分であるvHF_imgを測定し、測定したvHF_imgを用いて、次式に基づいてd軸インダクタンスLdを算出することができる。 Here, in the period 2 in which the alternating current is input from the minute current superimposing unit 12, the voltage vector in the motor 1 has a voltage component caused by the inductance L, which is the current vector to be superimposed. A 90 degree phase lead component v HF_img is included. By measuring vHF_img , which is a component contributed by this inductance, the measured vHF_img can be used to calculate the d-axis inductance Ld based on the following equation.

Figure 0007324630000008
Figure 0007324630000008

インダクタンス推定部71には、ロックインアンプなどにより検出されたvHF_imgが入力される。さらに、インダクタンス推定部71は、予め記憶している微小電流重畳部12にて生成する交流成分の周波数fHF、振幅iHFを記憶している。インダクタンス推定部71は、Trigが入力されている区間2において、式(8)に基づいて、d軸インダクタンスLdを推定することができる。 v HF_img detected by a lock-in amplifier or the like is input to the inductance estimator 71 . Further, the inductance estimator 71 stores the frequency f HF and the amplitude i HF of the AC component generated by the minute current superimposing unit 12, which are stored in advance. The inductance estimator 71 can estimate the d-axis inductance Ld based on the equation (8) in the interval 2 in which Trig is input.

そして、磁石磁束鎖交数決定部22は、次式を用いて磁石磁束鎖交数ψaを推定する。 Then, the magnet flux linkage number determination unit 22 estimates the magnet flux linkage number ψ a using the following equation.

Figure 0007324630000009
Figure 0007324630000009

(9)式は、推定されるd軸磁束鎖交数推定値ψdは、磁石磁束鎖交数ψaと電流に起因する磁束Ld・idとの和であることに起因する。このようにしても、磁石磁束鎖交数ψaを算出することができる。なお、磁石磁束鎖交数ψaの算出は、交流信号が重畳されていない期間1において行う方が、d軸電流値idが安定しており望ましい。しかしながら、これに限らず、期間2において行ってもよい。 Equation (9) results from the fact that the estimated d-axis flux linkage number ψ d is the sum of the magnet flux linkage number ψ a and the magnetic flux L d · id caused by the current. Also in this way, the magnet flux linkage number ψ a can be calculated. It should be noted that it is preferable to calculate the magnet flux linkage number ψ a during the period 1 in which the AC signal is not superimposed, because the d-axis current value id is stable. However, it is not limited to this, and may be performed in period 2.

第3実施形態によれば以下の効果を得ることができる。 According to the third embodiment, the following effects can be obtained.

第3実施形態のモータ1の磁石温度の推定方法によれば、微小電流重畳部12において、交流の微小電流が重畳される。このような交流の電流が重畳される場合においては、モータ1において観測される電圧ベクトルにおいて、インダクタンス成分に起因して、重畳される電流ベクトルに対して90度の位相が進んだ成分vHF_imgが含まれる。そこで、この観測されたvHF_imgと、重畳された交流の微小電流とによって、d軸インダクタンスLdを算出する。このように、d軸インダクタンスLdを直接的に算出することができるため、第1、第2実施形態のようにd軸インダクタンスLdが固定であると仮定する場合よりも、磁石磁束鎖交数ψaの推定精度を向上させることができる。 According to the method for estimating the magnet temperature of the motor 1 of the third embodiment, the minute current superimposing unit 12 superimposes a minute AC current. When such an alternating current is superimposed, in the voltage vector observed in the motor 1, the component v HF_img whose phase leads the superimposed current vector by 90 degrees due to the inductance component is included. Therefore, the d-axis inductance L d is calculated from the observed v HF_img and the superimposed minute AC current. In this way, since the d-axis inductance Ld can be directly calculated, the magnet flux linkage is greater than the case where the d-axis inductance Ld is assumed to be fixed as in the first and second embodiments. The estimation accuracy of the number ψ a can be improved.

(第4実施形態)
第1乃至3実施形態においては、磁石温度推定部17において、予め、d軸電流値id、q軸電流値iq、及び、磁石磁束鎖交数ψaと、磁石温度Tmagとを対応つけたテーブルを記憶する例について説明したが、これに限らない。第4実施形態においては、磁石温度推定部17において、ニューラルネットワークを用いて磁石温度Tmagを推測する例について説明する。
(Fourth embodiment)
In the first to third embodiments, the magnet temperature estimator 17 previously associates the d-axis current value i d , the q-axis current value i q , and the magnet flux linkage number ψ a with the magnet temperature T mag . Although an example of storing the attached table has been described, the present invention is not limited to this. In the fourth embodiment, an example in which the magnet temperature estimator 17 estimates the magnet temperature Tmag using a neural network will be described.

図8は、ニューラルネットワークの構成を示す例である。ニューラルネットワークにおいては入力層(Input Layer)と出力層(Output Layer)との間に複数の隠れ層(Hidden Layer)が設けられている。本実施形態においては、入力層には磁石磁束鎖交数ψaが入力され、出力層からは磁石温度Tmagが出力される。また、隠れ層においては、ニューロン1~4が示されている。なお、この例では、ニューロンにおける発火関数は、非線形関数である1/(1+e-ax)が用いられるものとする。 FIG. 8 is an example showing the configuration of a neural network. In a neural network, a plurality of hidden layers are provided between an input layer and an output layer. In this embodiment, the magnet flux linkage number ψ a is input to the input layer, and the magnet temperature T mag is output from the output layer. Also, in the hidden layer, neurons 1-4 are shown. In this example, it is assumed that the nonlinear function 1/(1+e -ax ) is used as the firing function in the neuron.

この例においては、入力層から、それぞれ入力重み付け1~4(Win-1~Win-4)を乗じて、ニューロン1~4へと入力される。そして、ニューロン1~4においては、それぞれオフセット1~4が加算された後に、出力重み付け1~4(Wout-1~Wout-4)がなされて、出力層へと出力される。なお、出力層においても、オフセット(Offset_out)が付与される。そのため、推定される磁石温度Tmagは以下のように示すことができる。 In this example, the input layers are multiplied by input weights 1 to 4 (Win-1 to Win-4) and input to neurons 1 to 4, respectively. In neurons 1 to 4, after offsets 1 to 4 are added, output weighting 1 to 4 (Wout-1 to Wout-4) is performed and output to the output layer. Also in the output layer, an offset (Offset_out) is given. Therefore, the estimated magnet temperature T mag can be expressed as follows.

Figure 0007324630000010
Figure 0007324630000010

このように、入力に対して重みの乗算、オフセット値の加算、及び、非線形関数を用いることで、入力データから出力データを近似的に求めることができる。実際の磁石温度推定部17においては、これら重み、オフセット、非線形関数を行列として実装することでニューラルネットワークを実現できる。なおニューラルネットワークに限らずカーネル法等の他の方法を用いても良い。 Thus, by multiplying the input by the weight, adding the offset value, and using the nonlinear function, the output data can be obtained approximately from the input data. In the actual magnet temperature estimator 17, a neural network can be realized by implementing these weights, offsets, and nonlinear functions as matrices. It should be noted that other methods such as the kernel method may be used instead of the neural network.

第4実施形態によれば、このように、定数、および非線形関数を成分とする行列を用いて、ニューラルネットワーク等の方法により磁石磁束鎖交数ψaから、磁石温度Tmagを推定することで、磁石温度推定部17において磁石磁束鎖交数に関するテーブルを記憶する必要がなくなるので、磁石温度推定部17の構成を簡略化することができる。 According to the fourth embodiment, by estimating the magnet temperature T mag from the magnet flux linkage number ψ a by a method such as a neural network using a matrix having constants and nonlinear functions as components, Since the magnet temperature estimator 17 does not need to store a table relating to the magnet flux linkage number, the configuration of the magnet temperature estimator 17 can be simplified.

(第5実施形態)
第5実施形態においては、モータ1に対する負荷電流制御部15に対する補正として、推定した磁束鎖交数推定値ψd ^、ψq ^に応じた補正を行う例について説明する。なお、この補正は、q軸電流値iqを補正することにより行われる。
(Fifth embodiment)
In the fifth embodiment, as a correction to the load current control unit 15 for the motor 1, an example of performing a correction according to the estimated magnetic flux linkage number estimated values ψ d ̂ and ψ q ̂ will be described. This correction is performed by correcting the q-axis current value iq .

ここで、モータ1において生じるトルクTは、磁束鎖交数推定値ψd ^、ψq ^と電流値id、iqとの外積として推定できることが知られており、次式のように示せることが知られている。なお、次式で推定されるトルクは、1極対あたりに生じるトルクである。 Here, it is known that the torque T generated in the motor 1 can be estimated as the outer product of the estimated values of magnetic flux linkage ψ d ^ and ψ q ̂ and the current values id and i q , and can be expressed by the following equation. It is known. Note that the torque estimated by the following equation is the torque generated per pole pair.

Figure 0007324630000011
Figure 0007324630000011

図9は、負荷電流制御部15に対してトルク誤差補償を行う場合のブロック図である。 FIG. 9 is a block diagram when torque error compensation is performed on the load current control unit 15. As shown in FIG.

この構成のうち、磁束鎖交数推定部21から、乗算器901、902、加算器903、及び、乗算器904までの構成が、式(11)に相当する。すなわち乗算器901によるd軸磁束鎖交数推定値ψd ^とq軸電流値iqとの積と、乗算器902によるq軸磁束鎖交数推定値ψq ^とd軸電流値idとの積とが、加算器903によって加算されることで、1極対あたりのトルクが算出される。さらに後段の乗算器904によれば、極対数が乗ざれてトルク推定値T^が算出される。 Of this configuration, the configuration from the flux linkage number estimating section 21 to the multipliers 901 and 902, the adder 903, and the multiplier 904 corresponds to Equation (11). That is, the product of the d-axis flux linkage number estimated value ψ d ^ and the q-axis current value i q by the multiplier 901 and the q-axis flux linkage number estimate value ψ q ^ and the d-axis current value i d by the multiplier 902 is added by the adder 903 to calculate the torque per pole pair. Furthermore, according to the multiplier 904 in the latter stage, the torque estimated value T̂ is calculated by multiplying by the number of pole pairs.

そして、減算器905において、トルク指令値T*からトルク推定値T^が減じられたものがトルク誤差補償部91に入力される。トルク誤差補償部91においては、トルク推定値T^がトルク指令値T*に追従するようにq軸電流補正値が算出され、加算器906においてq軸電流指令値iq *に加算される。 Then, the subtractor 905 subtracts the estimated torque value T^ from the torque command value T * , and this result is input to the torque error compensator 91 . In the torque error compensator 91, the q-axis current correction value is calculated so that the estimated torque T^ follows the torque command value T * , and is added to the q-axis current command value iq * in the adder 906.

そして、トルク誤差補償部91の前段において、トルク指令値T*に対してトルク推定値T^を減算することでトルク誤差が算出され、トルク誤差補償部91においてはそのトルク誤差が補償されるようなq軸電流補正値が求められる。そして、トルク誤差補償部91の後段において、q軸電流指令値iq *に対してq軸電流補正値が減算補正される。q軸電流値iqはトルク制御に適しているだけでなく、磁石温度測定に用いられる微小電流はd軸電流値idに重畳されておりq軸電流値iqには電流補正値が加減されるだけの余裕がある。そのため、q軸電流値iqを制御することで、トルク誤差を補正することができる。 A torque error is calculated by subtracting the estimated torque value T^ from the torque command value T * in the preceding stage of the torque error compensating unit 91, and the torque error compensating unit 91 compensates for the torque error. q-axis current correction value is obtained. Then, in the subsequent stage of the torque error compensator 91, the q-axis current correction value is subtracted from the q-axis current command value i q * . The q-axis current value i q is suitable not only for torque control, but also because the minute current used for magnet temperature measurement is superimposed on the d-axis current value i d , the q-axis current value i q can be adjusted by the current correction value. I can afford to be. Therefore, the torque error can be corrected by controlling the q-axis current value i q .

第5実施形態によれば以下の効果を得ることができる。 According to the fifth embodiment, the following effects can be obtained.

第5実施形態のモータ1の磁石温度の推定方法によれば、推定した磁石温度Tmagを用いてモータ1において生じているトルク推定値T^を推定し、トルク推定値T^がトルク指令値T*に追従するようにq軸電流指令値iq *に対して補正を行う。磁石温度Tmagの測定は、モータ1のトルクへの影響が少ないd軸電流値idに微小電流を重畳させており、q軸電流値iqはトルク制御のための余裕があるため比較的トルク補正を行いやすい。このようにq軸電流指令値iq *に対して補正をすることにより、磁石温度Tmagが考慮されたモータ1の制御が行われるので、トルク変動を低減するとともに、静音化を図ることができる。 According to the method of estimating the magnet temperature of the motor 1 of the fifth embodiment, the estimated magnet temperature T mag is used to estimate the estimated torque value T^ occurring in the motor 1, and the estimated torque value T^ is the torque command value. The q-axis current command value i q * is corrected so as to follow T * . The magnet temperature Tmag is measured by superimposing a minute current on the d -axis current value id, which has little effect on the torque of the motor 1, and the q-axis current value iq has a margin for torque control, so it is relatively small. Torque correction is easy. By correcting the q-axis current command value i q * in this way, the motor 1 is controlled in consideration of the magnet temperature Tmag , so that torque fluctuations can be reduced and noise can be reduced. can.

以上、本発明の実施形態について説明したが、上記実施形態は本発明の適用例の一部を示したに過ぎず、本発明の技術的範囲を上記実施形態の具体的構成に限定する趣旨ではない。また、上記実施形態は、適宜組み合わせ可能である。 Although the embodiments of the present invention have been described above, the above embodiments merely show a part of application examples of the present invention, and the technical scope of the present invention is not limited to the specific configurations of the above embodiments. do not have. Moreover, the above-described embodiments can be combined as appropriate.

10 モータシステム
1 モータ
11 電流指令値算出部
11A 補正部
12 微小電流重畳部
15 負荷電流制御部
16 磁石磁束鎖交数算出部
17 磁石温度推定部
21 磁束鎖交数推定部
22 磁石磁束鎖交数決定部
REFERENCE SIGNS LIST 10 motor system 1 motor 11 current command value calculation unit 11A correction unit 12 minute current superimposition unit 15 load current control unit 16 magnet flux linkage calculation unit 17 magnet temperature estimation unit 21 magnetic flux linkage number estimation unit 22 magnet flux linkage number decision part

Claims (12)

トルク指令値に応じた目標電流を流すことによりモータの回転を制御するモータの磁石温度の推定方法であって、
前記目標電流に該目標電流よりも小さいd軸の微小電流を重畳することで該目標電流を変化させ、
前記微小電流を前記目標電流に重畳させる前の第1状態、及び、前記微小電流を前記目標電流に重畳させた後の第2状態のそれぞれにおいてd軸磁束鎖交数を求め、
前記第1状態、及び、前記第2状態における前記d軸磁束鎖交数に基づいて、回転子が備える永久磁石の磁石磁束鎖交数を算出し、
算出した前記磁石磁束鎖交数に応じて前記永久磁石の磁石温度を推定する、モータの磁石温度の推定方法。
A method for estimating a magnet temperature of a motor for controlling the rotation of the motor by applying a target current according to a torque command value,
changing the target current by superimposing a d-axis minute current smaller than the target current on the target current ;
Obtaining the d-axis magnetic flux linkage number in each of a first state before superimposing the minute current on the target current and a second state after the minute current is superimposed on the target current;
based on the d-axis magnetic flux linkage numbers in the first state and the second state, calculating the magnetic flux linkage numbers of the permanent magnets included in the rotor;
A method for estimating the magnet temperature of a motor, wherein the magnet temperature of the permanent magnet is estimated according to the calculated magnet flux linkage.
請求項1に記載のモータの磁石温度の推定方法であって、
前記モータに流れる電圧の測定値と、前記目標電流に応じた電圧指令値との差分を求め、
前記差分を積分することにより、前記d軸磁束鎖交数を含む磁束鎖交数を求める、モータの磁石温度の推定方法。
A method for estimating a magnet temperature of a motor according to claim 1,
obtaining a difference between a measured value of the voltage flowing through the motor and a voltage command value corresponding to the target current;
A method for estimating a magnet temperature of a motor, wherein a magnetic flux linkage number including the d-axis magnetic flux linkage number is obtained by integrating the difference.
請求項に記載のモータの磁石温度の推定方法であって、
前記微小電流は、直流であり、
前記第1状態の前記d軸磁束鎖交数、前記第2状態における前記d軸磁束鎖交数、及び、前記モータに流れるd軸電流に基づいて前記磁石磁束鎖交数を算出する、モータの磁石温度の推定方法。
A method for estimating the magnet temperature of a motor according to claim 2 ,
the minute current is a direct current,
calculating the magnet flux linkage number based on the d-axis magnetic flux linkage number in the first state, the d-axis magnetic flux linkage number in the second state , and the d-axis current flowing through the motor; method for estimating the magnet temperature.
請求項に記載のモータの磁石温度の推定方法であって、
前記磁石磁束鎖交数は、式(1)により算出される、モータの磁石温度の推定方法。
Figure 0007324630000012

は、前記第1状態における前記d軸電流であり、
ψは、前記第1状態における前記d軸磁束鎖交数であり、
’は、前記第2状態における前記d軸電流であり、
ψ’は、前記第2状態における前記d軸磁束鎖交数である。
A method for estimating the magnet temperature of a motor according to claim 3 ,
The method for estimating the magnet temperature of the motor, wherein the magnet flux linkage number is calculated by Equation (1).
Figure 0007324630000012

id is the d-axis current in the first state;
ψ d is the d-axis magnetic flux linkage number in the first state,
i d ' is the d-axis current in the second state;
ψ d ′ is the d-axis magnetic flux linkage number in the second state.
請求項に記載のモータの磁石温度の推定方法であって、
2以上の大きさの異なる前記微小電流を用いて、2以上の前記磁石磁束鎖交数を算出し、
算出された2以上の前記磁石磁束鎖交数を平均することで、前記磁石磁束鎖交数を算出する、モータの磁石温度の推定方法。
A method for estimating the magnet temperature of a motor according to claim 2 ,
Calculating two or more magnet flux linkages using the minute currents of two or more different magnitudes,
A method for estimating a magnet temperature of a motor, wherein the magnet flux linkage number is calculated by averaging two or more calculated magnet flux linkage numbers.
請求項に記載のモータの磁石温度の推定方法であって、
前記微小電流は、交流であり、
前記第2状態において、前記モータの電圧ベクトルにおける前記微小電流のベクトルに対する90度の進み成分と、前記微小電流の振幅及び周波数とに基づいて、前記モータのインダクタンスを求め、
求められた前記d軸磁束鎖交数から、前記インダクタンスと前記目標電流との積を減じることで、前記磁石磁束鎖交数を算出する、モータの磁石温度の推定方法。
A method for estimating the magnet temperature of a motor according to claim 2 ,
the minute current is an alternating current,
In the second state, obtaining the inductance of the motor based on the 90-degree lead component of the voltage vector of the motor with respect to the vector of the minute current and the amplitude and frequency of the minute current;
A method for estimating the magnet temperature of a motor, wherein the magnet flux linkage number is calculated by subtracting the product of the inductance and the target current from the obtained d-axis flux linkage number.
請求項に記載のモータの磁石温度の推定方法であって、
前記インダクタンスは、式(2)により求められ、前記磁石磁束鎖交数は、式(3)により算出される、モータの磁石温度の推定方法。
Figure 0007324630000013

Figure 0007324630000014

は、前記インダクタンスであり、
HF_imgは、前記第2状態において、前記モータへ印加される電圧の電流に対する90度の進み成分であり、
HFは、前記微小電流の周波数であり、
HFは、前記微小電流の振幅であり、
は、前記モータに流れるd軸電流であり、
ψは、求められた前記d軸磁束鎖交数である。
A method for estimating the magnet temperature of a motor according to claim 6 ,
A method for estimating a magnet temperature of a motor, wherein the inductance is obtained by equation (2), and the magnet flux linkage is calculated by equation (3).
Figure 0007324630000013

Figure 0007324630000014

L d is the inductance;
v HF_img is the leading component of the voltage applied to the motor by 90 degrees with respect to the current in the second state;
f HF is the frequency of said minute current;
i HF is the amplitude of said minute current;
id is the d-axis current flowing through the motor;
ψ d is the obtained d-axis magnetic flux linkage number.
請求項に記載のモータの磁石温度の推定方法であって、
あらかじめ記憶された、前記磁石磁束鎖交数、前記モータに流れるd軸電流値、及び、q軸電流値のパラメータと、前記磁石温度との対応関係を示すテーブル用いて、入力される前記パラメータに応じて、前記永久磁石の温度を推定する、モータの磁石温度の推定方法。
A method for estimating a magnet temperature of a motor according to claim 1 ,
Using a pre -stored table showing the correspondence relationship between the parameters of the magnet flux linkage, the d-axis current value and the q-axis current value flowing through the motor , and the magnet temperature, the parameters to be input are A method for estimating the magnet temperature of a motor, estimating the temperature of said permanent magnet accordingly.
請求項に記載のモータの磁石温度の推定方法であって、
入力される前記パラメータが、前記テーブルに含まれていない場合には、該含まれていない前記パラメータを前記テーブルにおいて補完する、モータの磁石温度の推定方法。
A method for estimating the magnet temperature of a motor according to claim 8 ,
A method for estimating a magnet temperature of a motor, wherein when the input parameter is not included in the table, the parameter not included is complemented in the table.
請求項に記載のモータの磁石温度の推定方法であって、
ニューラルネットワークを用いて、前記磁石磁束鎖交数に応じて前記磁石温度を推定する、モータの磁石温度の推定方法。
A method for estimating a magnet temperature of a motor according to claim 1 ,
A method for estimating the magnet temperature of a motor, wherein the magnet temperature is estimated according to the magnet flux linkage number using a neural network.
請求項に記載のモータの磁石温度の推定方法であって、
前記磁束鎖交数、及び、前記磁石温度に応じて、前記モータにおけるトルク推定値を推定し、
前記トルク指令値と、前記トルク推定値との偏差が小さくなるように、前記モータに入力されるq軸の前記目標電流を補正する、モータの磁石温度の推定方法。
A method for estimating the magnet temperature of a motor according to claim 2 ,
estimating an estimated torque value in the motor according to the magnetic flux linkage and the magnet temperature;
A method for estimating a magnet temperature of a motor, wherein the target current of the q-axis input to the motor is corrected so that a deviation between the torque command value and the torque estimation value is reduced.
トルク指令値に応じた目標電流をモータに流すことにより前記モータの回転を制御するモータの磁石温度の推定装置であって、
前記目標電流に該目標電流よりも小さいd軸の微小電流を重畳することで該目標電流を変化させる変化部と、
前記微小電流を前記目標電流に重畳させる前の第1状態、及び、前記微小電流を前記目標電流に重畳させた後の第2状態のそれぞれにおいてd軸磁束鎖交数を求める磁束鎖交数算出部と、
前記第1状態、及び、前記第2状態における前記d軸磁束鎖交数に基づいて、回転子が備える永久磁石の磁石磁束鎖交数を算出する磁石磁束鎖交数算出部と、
算出した前記磁石磁束鎖交数に応じて前記永久磁石の磁石温度を推定する磁石温度推定部と、を備えるモータの磁石温度の推定装置。
A magnet temperature estimating device for controlling the rotation of the motor by causing a target current corresponding to a torque command value to flow through the motor,
a changing unit that changes the target current by superimposing a d-axis minute current that is smaller than the target current on the target current ;
Magnetic flux linkage calculation for obtaining a d-axis magnetic flux linkage in each of a first state before superimposing the minute current on the target current and a second state after the minute current is superimposed on the target current. Department and
a magnet flux linkage number calculation unit that calculates a magnet flux linkage number of a permanent magnet included in a rotor based on the d-axis flux linkage numbers in the first state and the second state;
A magnet temperature estimating device for a motor, comprising: a magnet temperature estimating unit that estimates the magnet temperature of the permanent magnet according to the calculated magnet flux linkage.
JP2019128757A 2019-07-10 2019-07-10 Motor magnet temperature estimation method and magnet temperature estimation device Active JP7324630B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019128757A JP7324630B2 (en) 2019-07-10 2019-07-10 Motor magnet temperature estimation method and magnet temperature estimation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019128757A JP7324630B2 (en) 2019-07-10 2019-07-10 Motor magnet temperature estimation method and magnet temperature estimation device

Publications (2)

Publication Number Publication Date
JP2021016226A JP2021016226A (en) 2021-02-12
JP7324630B2 true JP7324630B2 (en) 2023-08-10

Family

ID=74531920

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019128757A Active JP7324630B2 (en) 2019-07-10 2019-07-10 Motor magnet temperature estimation method and magnet temperature estimation device

Country Status (1)

Country Link
JP (1) JP7324630B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT202100012575A1 (en) * 2021-05-17 2022-11-17 Eldor Corp Spa METHOD OF ESTIMING THE TEMPERATURE OF MAGNETS IN AN ELECTRIC MACHINE, A METHOD OF CONTROLLING THE SAID ELECTRIC MACHINE, ELECTRONIC DEVICE FOR ESTIMING THE TEMPERATURE OF MAGNETS IN AN ELECTRIC MACHINE AND CONTROL SYSTEM FOR THE SAID ELECTRIC MACHINE
JP7452511B2 (en) 2021-10-01 2024-03-19 株式会社豊田中央研究所 motor control device

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008148546A (en) 2006-11-14 2008-06-26 Mitsubishi Electric Corp Controller for synchronous machine, and method of controlling the synchronous machine
JP2014222954A (en) 2013-05-13 2014-11-27 三菱電機株式会社 Synchronous machine control device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008148546A (en) 2006-11-14 2008-06-26 Mitsubishi Electric Corp Controller for synchronous machine, and method of controlling the synchronous machine
JP2014222954A (en) 2013-05-13 2014-11-27 三菱電機株式会社 Synchronous machine control device

Also Published As

Publication number Publication date
JP2021016226A (en) 2021-02-12

Similar Documents

Publication Publication Date Title
US8115441B2 (en) On-line measurement of an induction machine's rotor time constant by small signal d-axis current injection
JP4531751B2 (en) Synchronous machine controller
US6281659B1 (en) Induction motor drive and a parameter estimation method thereof
KR102017806B1 (en) Extended Luenberger-Sliding Mode Observer Capable of Estimating Rotor flux and Rotor Resistance for Three Phase Induction Motor
JP5327700B2 (en) Induction motor control device and control method thereof
JPH09219999A (en) Variable speed drive device
WO2016125559A1 (en) Motor control device
JP7324630B2 (en) Motor magnet temperature estimation method and magnet temperature estimation device
JP3253004B2 (en) Method of estimating speed of permanent magnet type synchronous motor, method of estimating rotor misalignment angle, and method of correcting rotor position
JPWO2013141059A1 (en) Control device for three-phase AC induction motor and control method for three-phase AC induction motor
JP2002252991A (en) Motor control apparatus
Chatterjee A simple leakage inductance identification technique for three-phase induction machines under variable flux condition
JP2017216807A (en) Vector control compensation method for induction motor, and vector control device
JP6033381B2 (en) Induction motor control device
JP3920750B2 (en) DC brushless motor control device
JP5050387B2 (en) Motor control device
JP6626309B2 (en) Control device for compensating iron loss of AC motor
JP2016096666A (en) Motor control apparatus
JP2010081726A (en) Controller of ac motor and control method thereof
JP4273775B2 (en) Magnetic pole position estimation method and control device for permanent magnet type synchronous motor
JP5925058B2 (en) Induction motor control device
JP6108114B2 (en) Control device for permanent magnet type synchronous motor
JP6308894B2 (en) Control device for three-phase AC rotating machine
CN109842337B (en) Flux linkage difference processing method and motor control device
JP2013183558A (en) Motor controller

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220510

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230307

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230308

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230509

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20230704

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230731

R150 Certificate of patent or registration of utility model

Ref document number: 7324630

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150