JP7163641B2 - Synchronous motor controller - Google Patents

Synchronous motor controller Download PDF

Info

Publication number
JP7163641B2
JP7163641B2 JP2018128517A JP2018128517A JP7163641B2 JP 7163641 B2 JP7163641 B2 JP 7163641B2 JP 2018128517 A JP2018128517 A JP 2018128517A JP 2018128517 A JP2018128517 A JP 2018128517A JP 7163641 B2 JP7163641 B2 JP 7163641B2
Authority
JP
Japan
Prior art keywords
axis
current
synchronous motor
fundamental wave
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.)
Active
Application number
JP2018128517A
Other languages
Japanese (ja)
Other versions
JP2020010477A (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.)
Fuji Electric Co Ltd
Original Assignee
Fuji Electric Co Ltd
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 Fuji Electric Co Ltd filed Critical Fuji Electric Co Ltd
Priority to JP2018128517A priority Critical patent/JP7163641B2/en
Publication of JP2020010477A publication Critical patent/JP2020010477A/en
Application granted granted Critical
Publication of JP7163641B2 publication Critical patent/JP7163641B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Ac Motors In General (AREA)

Description

本発明は、同期電動機の制御装置が行うチューニングに関する。 The present invention relates to tuning performed by a control device for a synchronous motor.

埋め込み永久磁石型同期電動機(以下、IPMSM)において、d,q軸間の干渉を含む磁気飽和を考慮した磁束モデルを、オートチューニングにより最適化する方法が知られている。その一例として、回転子が停止した状態で、d軸とq軸とに、高周波の正弦波交番電流(以下、チューニング電流)を流すことで、IPMSMの磁気飽和特性に関するパラメータをチューニングする方法がある(例えば、特許文献1参照)。 A method of optimizing a magnetic flux model considering magnetic saturation including interference between d and q axes by auto-tuning in an embedded permanent magnet type synchronous motor (IPMSM) is known. As an example, there is a method of tuning the parameters related to the magnetic saturation characteristics of the IPMSM by applying high-frequency sinusoidal alternating currents (hereinafter referred to as tuning currents) to the d-axis and q-axis while the rotor is stopped. (See Patent Document 1, for example).

なお、磁気飽和特性とは、電流の増加に伴う電動機鉄心の磁気飽和により、d,q軸磁束とこれらに対応する各軸電流との線形性が崩れる特性をいい、d,q軸間の干渉とは、他軸電流の影響により自軸磁束が変化する特性をいう。 The magnetic saturation characteristic is the characteristic that the linearity between the d- and q-axis magnetic fluxes and the respective axis currents corresponding to them is lost due to the magnetic saturation of the motor iron core as the current increases. A characteristic is that the self-axis magnetic flux changes due to the influence of the other-axis current.

特開2015‐144502号公報JP 2015-144502 A

従来のオートチューニング方法では、チューニング電流が流れている間、同期電動機の回転子が停止していることを前提にしている。しかし、チューニング中の同期電動機にはチューニング電流に起因する高周波の振動トルクが発生することにより、磁極位置がチューニング開始時の位置からずれる場合がある。この位置ずれは、特に、慣性モーメントが小さい電動機で発生しやすい。 The conventional auto-tuning method assumes that the rotor of the synchronous motor is stopped while the tuning current is flowing. However, the magnetic pole position may deviate from the position at the start of tuning due to high-frequency vibration torque generated in the synchronous motor during tuning due to the tuning current. This positional deviation is particularly likely to occur in an electric motor having a small moment of inertia.

同期電動機が位置検出器または速度検出器を備えている場合、磁極位置ずれの検出は容易である。しかしながら、これらの検出器を使用せずに同期電動機を制御する位置・速度センサレス制御の場合、磁極位置ずれを直接検出することができないので、チューニング電流がd軸またはq軸からずれた位相に流れ、正しいチューニング結果が得られないおそれがある。 If the synchronous motor is equipped with a position detector or speed detector, it is easy to detect the magnetic pole misalignment. However, in the case of position/speed sensorless control that controls the synchronous motor without using these detectors, the magnetic pole position deviation cannot be detected directly, so the tuning current flows in a phase shifted from the d-axis or q-axis. , the correct tuning result may not be obtained.

そこで、本開示は、正しい位相にチューニング電流を流すことを可能にする、同期電動機の制御装置を提供する。 Accordingly, the present disclosure provides a controller for a synchronous motor that allows tuning currents to flow in the correct phase.

本開示の技術の一態様として、
同期電動機の回転子磁極方向に平行なd軸とこのd軸に直交するq軸とにより定められるd,q直交回転座標系とは別に定義される座標系であってγ軸とδ軸とにより定められるγ,δ直交回転座標系において、電力変換器により前記同期電動機に供給する電流及び電圧を制御する制御装置であって、
電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルに基づいて演算する演算装置を備え、
前記演算装置は、
前記同期電動機のγ軸電流の直流成分を指令値に制御し、
前記γ軸に正弦波の交番電圧を印加し、
前記d軸のインダクタンスの演算と前記磁束モデルのパラメータの演算との少なくとも一方を行い、
前記同期電動機のγ軸電圧の基本波成分のフーリエ係数を演算し、
前記同期電動機のδ軸電流の基本波成分のフーリエ係数を演算し、
前記γ軸電圧の基本波成分のフーリエ係数と、前記δ軸電流の基本波成分のフーリエ係数とを用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御することを特徴とする、同期電動機の制御装置が提供される。
As one aspect of the technology of the present disclosure,
A coordinate system defined separately from the d, q orthogonal rotating coordinate system defined by the d-axis parallel to the rotor magnetic pole direction of the synchronous motor and the q-axis orthogonal to the d-axis, and is defined by the γ-axis and the δ-axis A control device for controlling the current and voltage supplied to the synchronous motor by a power converter in a defined γ, δ orthogonal rotating coordinate system,
A computing device that performs computation based on a magnetic flux model that takes into account at least the magnetic saturation characteristics of an iron core of the motor,
The computing device is
controlling the DC component of the γ-axis current of the synchronous motor to a command value;
applying a sinusoidal alternating voltage to the γ-axis,
performing at least one of the calculation of the d-axis inductance and the calculation of the parameters of the magnetic flux model;
calculating the Fourier coefficient of the fundamental wave component of the γ-axis voltage of the synchronous motor;
calculating the Fourier coefficient of the fundamental wave component of the δ-axis current of the synchronous motor;
Using the Fourier coefficient of the fundamental wave component of the γ-axis voltage and the Fourier coefficient of the fundamental wave component of the δ-axis current, the γ, δ orthogonal rotating coordinate system is set such that the d-axis and the γ-axis are aligned. A control device for a synchronous motor is provided, characterized by controlling the phase.

また、本開示の技術の一態様として、
同期電動機の回転子磁極方向に平行なd軸とこのd軸に直交するq軸とにより定められるd,q直交回転座標系とは別に定義される座標系であってγ軸とδ軸とにより定められるγ,δ直交回転座標系において、電力変換器により前記同期電動機に供給する電流及び電圧を制御する制御装置であって、
電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルに基づいて演算する演算装置を備え、
前記演算装置は、
前記同期電動機のγ軸電流の直流成分を指令値に制御し、
前記δ軸に正弦波の交番電圧を印加し、
前記q軸のインダクタンスの演算と前記磁束モデルのパラメータの演算との少なくとも一方を行い、
前記同期電動機のδ軸電圧の基本波成分のフーリエ係数を演算し、
前記同期電動機のγ軸電流の基本波成分のフーリエ係数を演算し、
前記δ軸電圧の基本波成分のフーリエ係数と、前記γ軸電流の基本波成分のフーリエ係数とを用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御することを特徴とする、同期電動機の制御装置が提供される。
Further, as one aspect of the technology of the present disclosure,
A coordinate system defined separately from the d, q orthogonal rotating coordinate system defined by the d-axis parallel to the rotor magnetic pole direction of the synchronous motor and the q-axis orthogonal to the d-axis, and is defined by the γ-axis and the δ-axis A control device for controlling the current and voltage supplied to the synchronous motor by a power converter in a defined γ, δ orthogonal rotating coordinate system,
A computing device that performs computation based on a magnetic flux model that takes into account at least the magnetic saturation characteristics of an iron core of the motor,
The computing device is
controlling the DC component of the γ-axis current of the synchronous motor to a command value;
Applying a sinusoidal alternating voltage to the δ axis,
performing at least one of the calculation of the q-axis inductance and the calculation of the parameters of the magnetic flux model;
calculating the Fourier coefficient of the fundamental wave component of the δ-axis voltage of the synchronous motor;
calculating the Fourier coefficient of the fundamental wave component of the γ-axis current of the synchronous motor;
Using the Fourier coefficient of the fundamental wave component of the δ-axis voltage and the Fourier coefficient of the fundamental wave component of the γ-axis current, the γ, δ orthogonal rotating coordinate system is set such that the d-axis and the γ-axis are aligned. A control device for a synchronous motor is provided, characterized by controlling the phase.

また、本開示の技術の一態様として、
同期電動機の回転子磁極方向に平行なd軸とこのd軸に直交するq軸とにより定められるd,q直交回転座標系とは別に定義される座標系であってγ軸とδ軸とにより定められるγ,δ直交回転座標系において、電力変換器により同期電動機に供給する電流及び電圧を制御する演算装置を備え、
前記演算装置は、
前記γ軸と前記δ軸のうちチューニング電流を流す一方の軸に直交する他方の軸に発生する振動電流を検出し、前記振動電流の検出値に基づいて、チューニング電流を流す位相を補正する、同期電動機の制御装置が提供される。
Further, as one aspect of the technology of the present disclosure,
A coordinate system defined separately from the d, q orthogonal rotating coordinate system defined by the d-axis parallel to the rotor magnetic pole direction of the synchronous motor and the q-axis orthogonal to the d-axis, and is defined by the γ-axis and the δ-axis Equipped with a computing device that controls the current and voltage supplied to the synchronous motor by the power converter in the determined γ, δ orthogonal rotating coordinate system,
The computing device is
Detecting an oscillating current generated in the other of the γ-axis and the δ-axis perpendicular to one of the axes through which the tuning current flows, and correcting the phase of the tuning current flow based on the detected value of the oscillating current; A controller for a synchronous motor is provided.

本開示の技術によれば、正しい位相にチューニング電流を流すことが可能となる。 According to the technique of the present disclosure, it is possible to pass the tuning current in the correct phase.

本開示に係る同期電動機の制御装置の全体構成を例示するブロック図である。1 is a block diagram illustrating the overall configuration of a control device for a synchronous motor according to the present disclosure; FIG. 第1の実施形態で使用する軸ずれ補償器の構成例を示すブロック図である。FIG. 2 is a block diagram showing a configuration example of an axis deviation compensator used in the first embodiment; FIG. 第2の実施形態で使用する軸ずれ補償器の構成例を示すブロック図である。FIG. 10 is a block diagram showing a configuration example of an axis deviation compensator used in the second embodiment; 本実施形態で使用するγ‐δ軸の説明図である。FIG. 4 is an explanatory diagram of γ-δ axes used in this embodiment; 制御装置が備える演算装置のハードウェア構成を例示する図である。It is a figure which illustrates the hardware constitutions of the arithmetic unit with which a control apparatus is provided.

以下、図面を参照しながら本開示に係る同期電動機の制御装置の実施形態を説明する。同一の構成要素については同一の符号を付け、重複する説明は省略する。なお、本発明は下記の実施形態に限定されるものではなく、その要旨を変更しない範囲内で適宜変形して実施することができるものである。 Hereinafter, embodiments of a synchronous motor control device according to the present disclosure will be described with reference to the drawings. The same reference numerals are given to the same components, and redundant explanations are omitted. It should be noted that the present invention is not limited to the following embodiments, and can be modified as appropriate without changing the gist of the invention.

IPMSMやシンクロナスリラクタンスモータ(以下、SynRM)などの同期電動機のトルクを高精度に制御するためには、電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルを求め、この磁束モデルに基づいて電流制御を行うことが望ましい。磁束モデルに基づいて電流制御を行うためには、その磁束モデルの各パラメータを算出する必要がある。また、同期電動機の回転を高精度に制御しようとする場合、同期電動機のインダクタンスを算出することが必要になることがある。 In order to control the torque of a synchronous motor such as an IPMSM or a synchronous reluctance motor (hereinafter referred to as SynRM) with high accuracy, a magnetic flux model that takes into account at least the magnetic saturation characteristics of the motor core is obtained, and the current is controlled based on this magnetic flux model. It is desirable to In order to perform current control based on the magnetic flux model, it is necessary to calculate each parameter of the magnetic flux model. Moreover, when trying to control the rotation of the synchronous motor with high accuracy, it may be necessary to calculate the inductance of the synchronous motor.

本実施形態の同期電動機の制御装置は、同期電動機の回転子が静止した状態で交番電圧を印加し、その印加時に流れる電流の応答に基づいて、磁束モデルの各パラメータと同期電動機のインダクタンスとの少なくとも一方を演算するオートチューニングを行う。本実施形態の制御装置は、IPMSMやSynRMなどの同期電動機の回転子が停止した状態でチューニングを行う際に磁極位置にずれが発生した場合、チューニング電流を流す軸の対向軸に発生する高周波の振動電流を検出する。例えば、γ軸とδ軸のうち、チューニング電流を流す一方の軸がγ軸であれば、γ軸に直交する他方の軸(対向軸)はδ軸であり、チューニング電流を流す一方の軸がδ軸であれば、δ軸に直交する他方の軸(対向軸)はγ軸である。そして、本実施形態の制御装置は、その振動電流の検出値に基づき同期電動機の磁極位置のずれ量を演算し、そのずれ量の演算値に基づいてチューニング電流を流す位相を補正する。 The control device for a synchronous motor according to this embodiment applies an alternating voltage while the rotor of the synchronous motor is stationary, and based on the response of the current that flows when the alternating voltage is applied, determines the relationship between each parameter of the magnetic flux model and the inductance of the synchronous motor. Perform auto-tuning to calculate at least one of them. The control device of the present embodiment is designed to prevent high-frequency noise generated in the shaft opposite to the shaft through which the tuning current flows when the magnetic pole position is deviated when tuning is performed with the rotor of a synchronous motor such as IPMSM or SynRM stopped. Detect oscillating currents. For example, if one of the γ-axis and δ-axis through which the tuning current flows is the γ-axis, the other axis (opposite axis) orthogonal to the γ-axis is the δ-axis, and one of the axes through which the tuning current flows is the δ-axis. If it is the δ-axis, the other axis (opposite axis) perpendicular to the δ-axis is the γ-axis. Then, the control device of the present embodiment calculates the amount of deviation of the magnetic pole positions of the synchronous motor based on the detected value of the oscillating current, and corrects the phase of the tuning current based on the calculated value of the amount of deviation.

このような補正により、本実施形態の制御装置は、位置・速度検出器を持たない同期電動機の回転子が停止した状態でチューニングを行う際に磁極位置のずれが発生した場合であっても、磁極位置のずれを補償して正しい位相でチューニング電流を流すことができる。その結果、慣性モーメントが小さく磁極位置ずれが発生しやすい電動機が制御対象であっても、正しいチューニング結果を得ることが可能となる。 With such correction, the control device of the present embodiment can perform tuning even when a deviation of the magnetic pole position occurs when the rotor of a synchronous motor that does not have a position/speed detector is stopped. A tuning current can be supplied in the correct phase by compensating for the deviation of the magnetic pole position. As a result, it is possible to obtain a correct tuning result even if the control target is a motor that has a small moment of inertia and is likely to cause magnetic pole position deviation.

本実施形態の制御装置は、各パラメータのオートチューニングが完了した磁束モデルを同期電動機の電流制御に利用することで、同期電動機の回転を高精度に制御できる。例えば、制御装置は、オートチューニング後の各パラメータが反映された磁束モデルに基づいて算出されるd軸電流指令値i 及びq軸電流指令値i を用いて、周知のベクトル制御により同期電動機に供給する電流及び電圧を制御する。これにより、同期電動機の回転を高精度に制御することが可能となる。 The control device of the present embodiment can control the rotation of the synchronous motor with high precision by using the magnetic flux model for which auto-tuning of each parameter has been completed for current control of the synchronous motor. For example, the control device uses a d-axis current command value i d * and a q-axis current command value i d * calculated based on a magnetic flux model in which each parameter after auto-tuning is reflected. It controls the current and voltage supplied to the synchronous motor. This makes it possible to control the rotation of the synchronous motor with high accuracy.

次に、本実施形態の制御装置の詳細について説明する。最初に、本実施形態の制御装置が制御上で使用するγ‐δ軸について説明する。 Next, the details of the control device of this embodiment will be described. First, the γ-δ axes used for control by the control device of this embodiment will be described.

<1.γ‐δ軸>
図4は、d‐q軸とγ‐δ軸を表す図である。d‐q軸は、同期電動機の実際の磁極位置と同期するd,q直交回転座標系を定める直交軸あり、d軸を磁極位置方向と定める。電動機の回転子の磁極(N極)方向をd軸、回転子磁極方向に平行なd軸から90°進み方向をq軸と定義する。一方、γ‐δ軸は、d,q直交回転座標系とは別に定義される任意のγ,δ直交回転座標系を定める直交軸である。本実施形態の制御装置は、γ‐δ軸をd‐q軸に追従させることで位置・速度センサレス制御を実現する。
<1. γ-δ axis>
FIG. 4 is a diagram showing the dq axis and the γ-δ axis. The dq axis is an orthogonal axis that defines a d, q orthogonal rotating coordinate system synchronized with the actual magnetic pole position of the synchronous motor, and the d axis is defined as the magnetic pole position direction. The magnetic pole (N pole) direction of the rotor of the motor is defined as the d-axis, and the direction leading 90° from the d-axis parallel to the magnetic pole direction of the rotor is defined as the q-axis. On the other hand, the γ-δ axes are orthogonal axes defining an arbitrary γ, δ orthogonal rotating coordinate system defined separately from the d, q orthogonal rotating coordinate system. The control device of this embodiment realizes position/speed sensorless control by causing the γ-δ axes to follow the dq-axes.

本実施形態の制御装置は、γ‐δ軸上で同期電動機のチューニング電流を制御するが、チューニング電流の制御中に図4に示す軸ずれ量θerrが発生すると、d‐q軸とγ‐δ軸の間にずれが生じ、チューニング電流がd‐q軸とは異なる位置に制御される。次に、d軸とγ軸が一致するように、チューニング中に軸ずれ量θerrを補正する方法について説明する。 The control device of this embodiment controls the tuning current of the synchronous motor on the γ-δ axes . A shift occurs between the delta axes and the tuning current is controlled to a different position than the dq axes. Next, a method for correcting the axis deviation θ err during tuning so that the d-axis and the γ-axis are aligned will be described.

<2.第1の実施形態>
(2.1)電流及び電圧の制御方法
まず、図1は、本開示に係る制御装置を主回路と共に示したブロック図であり、以下では、永久磁石型同期電動機(以下、単に電動機又はSMともいう)の電圧及び電流の制御方法を制御装置の構成と共に説明する。なお、電力変換器により同期電動機に供給する電圧及び電流の制御演算は、d,q直交回転座標系に一致するように制御されるγ,δ直交回転座標系の上で行うこととする。
<2. First Embodiment>
(2.1) Current and voltage control method First, FIG. 1 is a block diagram showing a control device according to the present disclosure together with a main circuit. ) will be described together with the configuration of the control device. The control calculation of the voltage and current supplied to the synchronous motor by the power converter is performed on the γ, δ orthogonal rotating coordinate system controlled to match the d, q orthogonal rotating coordinate system.

第1の実施形態の制御装置100は、γ軸電流の直流成分を指令値に制御してγ軸に正弦波の交番電圧を印加しているときの電流及び電圧に基づき、電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルの各パラメータを演算するオートチューニングを行う。あるいは、制御装置100は、γ軸電流の直流成分を指令値に制御してγ軸に正弦波の交番電圧を印加しているときの電流及び電圧に基づき、同期電動機1の磁極の方向に平行な方向のd軸インダクタンスを演算(測定)するオートチューニングを行う。制御装置100は、磁束モデルの各パラメータとd軸インダクタンスとの両方を演算するオートチューニングを行ってもよい。 The controller 100 of the first embodiment controls the DC component of the γ-axis current to a command value, and controls the magnetic saturation of the motor iron core based on the current and voltage when the sinusoidal alternating voltage is applied to the γ-axis. Auto-tuning is performed to calculate each parameter of the magnetic flux model considering at least the characteristics. Alternatively, the control device 100 controls the DC component of the γ-axis current to a command value, and based on the current and voltage when a sinusoidal alternating voltage is applied to the γ-axis, the magnetic pole direction of the synchronous motor 1 is parallel to the direction of the magnetic poles of the synchronous motor 1 . Auto-tuning is performed to calculate (measure) the d-axis inductance in the direction. The control device 100 may perform auto-tuning to calculate both the parameters of the magnetic flux model and the d-axis inductance.

第1の実施形態において、例えば、γ軸電流直流成分指令値Iγ(0) を、0を含む任意の値に設定し、γ軸交番電流指令値の振幅Iγ(b1) を、0よりも大きい任意の値に設定し、δ軸交番電流指令値の振幅Iδ(b1) を、0に設定する。Iγ(0) ,Iγ(b1) ,Iδ(b1) は、それぞれ、γ軸電流の直流成分の指令値、γ軸に流す交番電流の振幅の指令値、δ軸に流す交番電流の振幅の指令値を表す。 In the first embodiment, for example, the γ-axis current DC component command value Iγ(0) * is set to an arbitrary value including 0, and the amplitude Iγ(b1) * of the γ-axis alternating current command value is set to An arbitrary value greater than 0 is set, and the amplitude I δ(b1) * of the δ-axis alternating current command value is set to 0. Iγ(0) * , Iγ(b1) * , and Iδ(b1) * are respectively the command value for the DC component of the γ-axis current, the command value for the amplitude of the alternating current flowing through the γ-axis, and the command value for the amplitude of the alternating current flowing through the δ-axis. Represents the command value for the amplitude of the alternating current.

第1の実施形態では、積分器5aは、交番電圧の角周波数ωを積分してd軸交番電流指令値の角度θを演算する。γ軸交番電流指令演算器6aは、γ軸交番電流指令値の交流成分iγh を数式1のように演算する。 In the first embodiment, the integrator 5a integrates the angular frequency ωh of the alternating voltage to calculate the angle θh of the d-axis alternating current command value. The γ-axis alternating current command calculator 6a calculates the AC component i γh * of the γ-axis alternating current command value as shown in Equation (1).

Figure 0007163641000001
Figure 0007163641000001

加算器7aは、γ軸電流直流成分指令値Iγ(0) とiγh を加算し、γ軸電流指令値iγ を演算する。 The adder 7a adds the γ-axis current DC component command value Iγ(0) * and iγh * to calculate the γ -axis current command value iγ * .

電流座標変換器8は、u相電流検出器9u及びw相電流検出器9wによりそれぞれ検出したu,w相電流検出値i,iを、軸ずれ補償器19により演算された同期電動機1の磁極位置推定値θr_estに基づいてγ,δ軸電流検出値iγ,iδに座標変換する。 The current coordinate converter 8 converts the u- and w-phase current detection values i u and i w respectively detected by the u-phase current detector 9 u and the w-phase current detector 9 w to the synchronous motor 1 based on the magnetic pole position estimation value θ r_est of , the coordinates are transformed into γ, δ axis current detection values i γ , i δ .

ローパスフィルタ10aは、γ軸電流検出値iγの高周波成分を除去してγ軸電流検出値iγfを演算する。 The low-pass filter 10a removes high-frequency components from the γ -axis current detection value to calculate the γ-axis current detection value iγf.

γ軸電流指令値iγ とγ軸電流検出値iγfとの偏差を減算器11aにて演算し、この偏差をγ軸電流調節器12aにより増幅してγ軸電圧指令値vγACR(γ軸電圧フィードバック制御値vγACR)を演算する。γ軸電流調節器12aは、γ軸電流指令値iγ とγ軸電流検出値iγfとの偏差が零になるように動作してγ軸電圧フィードバック制御値vγACRを演算する。 The deviation between the γ-axis current command value i γ * and the γ-axis current detection value i γf is calculated by the subtractor 11a, and this deviation is amplified by the γ-axis current controller 12a to obtain the γ-axis voltage command value v γACR (γ Calculate the shaft voltage feedback control value v γACR ). The γ-axis current controller 12a operates so that the deviation between the γ-axis current command value i γ * and the γ-axis current detection value i γf becomes zero to calculate the γ-axis voltage feedback control value v γACR .

磁束モデルのパラメータは、電動機鉄心の磁気飽和に起因して流れる高調波電流を利用して推定される。このため、γ軸電流調節器12aが高調波電流に作用しないようにするため、γ軸電流調節器12aの応答周波数、及び、ローパスフィルタ10aのカットオフ周波数は、交番電圧の角周波数ωの2倍よりも小さく設定される。 The parameters of the flux model are estimated using harmonic currents that flow due to magnetic saturation of the motor core. Therefore, in order to prevent the γ-axis current adjuster 12a from acting on harmonic currents, the response frequency of the γ-axis current adjuster 12a and the cutoff frequency of the low-pass filter 10a are adjusted to the angular frequency ω h of the alternating voltage. It is set smaller than twice.

直流電機子抵抗補償器13は、γ軸電機子抵抗フィードフォワード補償値vγraを数式2により演算する。 The DC armature resistance compensator 13 calculates the γ-axis armature resistance feedforward compensation value v γra by Equation (2).

Figure 0007163641000002
Figure 0007163641000002

電圧補償値演算器14は、インピーダンス推定器15で演算したd軸リアクタンス推定値Xdhest及びd軸電機子抵抗推定値Rdhestを使用し、γ軸電圧フィードフォワード補償値vγhFFを数式3により演算する。なお、d軸リアクタンス推定値Xdhest及びd軸電機子抵抗推定値Rdhestなどのインピーダンスの推定は、例えば特許文献1等に記載されている公知技術を用いてインピーダンス推定器15により行われることが可能である。 The voltage compensation value calculator 14 uses the d-axis reactance estimated value X dhest and the d-axis armature resistance estimated value R dhest calculated by the impedance estimator 15 to calculate the γ-axis voltage feedforward compensation value v γhFF according to Equation 3. do. Estimation of impedance such as the d-axis reactance estimated value X dhest and the d-axis armature resistance estimated value R dhest may be performed by the impedance estimator 15 using a known technique described in Patent Document 1, for example. It is possible.

Figure 0007163641000003
Figure 0007163641000003

また、d軸インダクタンスの演算は、公知技術を用いてインピーダンス推定器15により行われることが可能である。例えば、インピーダンス推定器15は、d軸リアクタンス推定値Xdhestを角周波数ωで除算することによって、d軸インダクタンスLdestを演算する(Ldest=Xdhest/ω)。 Also, the calculation of the d-axis inductance can be performed by the impedance estimator 15 using known techniques. For example, the impedance estimator 15 calculates the d-axis inductance L dest by dividing the d-axis reactance estimated value X dhest by the angular frequency ω h (L dest =X dhesth ).

フーリエ係数演算器16は、θ,vγ ,iγに基づき、γ軸電圧及びγ軸電流のそれぞれのフーリエ係数を演算する。フーリエ係数演算器16は、γ軸電圧の基本波成分のフーリエ余弦係数Vγ(a1)、γ軸電圧の基本波成分のフーリエ正弦係数Vγ(b1)、γ軸電流の直流成分のフーリエ係数Iγ(0)、γ軸電流の基本波成分のフーリエ余弦係数Iγ(a1)、γ軸電流の基本波成分のフーリエ正弦係数Iγ(b1)、γ軸電流の2倍高調波成分のフーリエ余弦係数Iγ(a2)、γ軸電流の2倍高調波成分のフーリエ正弦係数Iγ(b2)、γ軸電流の3倍高調波成分のフーリエ余弦係数Iγ(a3)、及び、γ軸電流の3倍高調波成分のフーリエ正弦係数Iγ(b3)を演算する。フーリエ係数の演算は、例えば特許文献1等に記載されている公知技術を用いて行われることが可能である。 The Fourier coefficient calculator 16 calculates Fourier coefficients of the γ-axis voltage and the γ -axis current based on θh, * , and . The Fourier coefficient calculator 16 calculates the Fourier cosine coefficient Vγ(a1) of the fundamental wave component of the γ-axis voltage, the Fourier sine coefficient Vγ(b1) of the fundamental wave component of the γ-axis voltage, and the Fourier coefficient of the DC component of the γ-axis current. I γ(0) , Fourier cosine coefficient I γ(a1) of the fundamental wave component of γ-axis current, Fourier sine coefficient I γ(b1) of the fundamental wave component of γ-axis current, Double harmonic component of γ-axis current Fourier cosine coefficient I γ(a2) , Fourier sine coefficient I γ(b2) of the second harmonic component of the γ-axis current, Fourier cosine coefficient I γ(a3) of the third harmonic component of the γ-axis current, and γ Calculate the Fourier sine coefficient Iγ(b3) of the triple harmonic component of the shaft current. Fourier coefficients can be calculated using a known technique described in Patent Document 1, for example.

フーリエ係数演算器16により演算された、γ軸電流の2倍高調波成分のフーリエ余弦係数Iγ(a2)、γ軸電流の2倍高調波成分のフーリエ正弦係数Iγ(b2)、γ軸電流の3倍高調波成分のフーリエ余弦係数Iγ(a3)、γ軸電流の3倍高調波成分のフーリエ正弦係数Iγ(b3)から、電機子抵抗補償器17は、数式4によりγ軸高調波電機子抵抗補償値vγhraを演算する。 The Fourier cosine coefficient Iγ(a2) of the double harmonic component of the γ-axis current, the Fourier sine coefficient Iγ(b2) of the double harmonic component of the γ-axis current, and the γ-axis calculated by the Fourier coefficient calculator 16 From the Fourier cosine coefficient I γ(a3) of the third harmonic component of the current and the Fourier sine coefficient I γ(b3) of the third harmonic component of the γ-axis current, the armature resistance compensator 17 can be calculated by Equation 4 as follows: A harmonic armature resistance compensation value vγhra is calculated.

Figure 0007163641000004
Figure 0007163641000004

磁束モデル係数推定器18には、フーリエ係数演算器16により演算された、γ軸電流直流成分のフーリエ係数Iγ(0)、γ軸電流基本波成分のフーリエ余弦係数Iγ(a1)、γ軸電流基本波成分のフーリエ正弦係数Iγ(b1)、γ軸電流の2倍高調波成分のフーリエ余弦係数Iγ(a2)、γ軸電流の2倍高調波成分のフーリエ正弦係数Iγ(b2)、γ軸電流の3倍高調波成分のフーリエ余弦係数Iγ(a3)、及び、γ軸電流の3倍高調波成分のフーリエ正弦係数Iγ(b3)が入力される。磁束モデル係数推定器18は、例えば数式5のようなd軸磁束モデルの係数(パラメータ)を演算する。 The magnetic flux model coefficient estimator 18 receives the Fourier coefficient I γ(0) of the γ-axis current DC component, the Fourier cosine coefficient I γ(a1) of the γ-axis current fundamental wave component, γ Fourier sine coefficient Iγ(b1) of the fundamental wave component of the axis current, Fourier cosine coefficient Iγ(a2) of the double harmonic component of the γ-axis current, Fourier sine coefficient Iγ(a2) of the double harmonic component of the γ -axis current b2) , the Fourier cosine coefficient Iγ(a3) of the third harmonic component of the γ-axis current, and the Fourier sine coefficient Iγ(b3) of the third harmonic component of the γ-axis current are input. The magnetic flux model coefficient estimator 18 calculates coefficients (parameters) of the d-axis magnetic flux model, such as Equation 5, for example.

Figure 0007163641000005
Figure 0007163641000005

第1の実施形態では、δ軸電流iδは0となるので、磁束モデル係数推定器18は、KLdおよびKsdを、上記のフーリエ係数を使用して、例えば代数演算や最適化演算等の公知の演算方法により演算する。例えば、磁束モデル係数推定器18は、特許文献1等に記載されている公知技術を用いて、KLdおよびKsd等の磁束モデルのパラメータを演算するオートチューニングを行う。 In the first embodiment, the δ-axis current i δ is 0, so the magnetic flux model coefficient estimator 18 calculates K Ld and K sd using the above Fourier coefficients, for example, algebraic calculations, optimization calculations, etc. is calculated by a known calculation method. For example, the magnetic flux model coefficient estimator 18 performs auto-tuning to calculate parameters of the magnetic flux model such as K Ld and K sd using the known technique described in Patent Document 1 and the like.

γ軸電圧指令値vγ は、加算器7b,7dにより、γ軸電圧フィードバック制御値vγACR、γ軸電機子抵抗フィードフォワード補償値vγra、γ軸電圧フィードフォワード補償値vγhFF、γ軸高調波電機子抵抗補償値vγhraを加算して算出される。一方、δ軸電圧指令値vδ は、零に固定する。 The γ-axis voltage command value v γ * is obtained by adders 7b and 7d into a γ-axis voltage feedback control value v γACR , a γ-axis armature resistance feedforward compensation value v γra , a γ-axis voltage feedforward compensation value v γhFF , a γ-axis It is calculated by adding the harmonic armature resistance compensation value v γhra . On the other hand, the δ-axis voltage command value v δ * is fixed at zero.

上述のように演算したγ,δ軸電圧指令値vγ ,vδ は、電圧座標変換器20によって、軸ずれ補償器19により演算された磁極位置推定値θr_estに基づいて、u,v,w相の相電圧指令値v ,v ,v に変換される。 The γ, δ axis voltage command values v γ * , v δ * calculated as described above are converted by the voltage coordinate converter 20 into u, It is converted into phase voltage command values v u * , v v * , v w * of the v and w phases.

整流回路3は、三相交流電源4からの三相交流電圧を整流して直流電圧に変換し、この直流電圧をインバータ等の電力変換器2に供給する。 The rectifier circuit 3 rectifies the three-phase AC voltage from the three-phase AC power supply 4, converts it into a DC voltage, and supplies this DC voltage to the power converter 2 such as an inverter.

PWM(Pulse Width Modulation)回路21は、相電圧指令値v ,v ,v に基づいて、電力変換器2の出力電圧を相電圧指令値v ,v ,v に制御するための複数のゲート信号を生成する。電力変換器2は、PWM回路21からの複数のゲート信号に基づいて、電力変換器2内部の複数の半導体スイッチング素子を制御することにより、IPMSMやSynRMなどの同期電動機1の端子電圧を相電圧指令値v ,v ,v に制御する。 PWM (Pulse Width Modulation) circuit 21 adjusts the output voltage of power converter 2 to phase voltage command values v u * , v v * , v based on phase voltage command values v u * , v v * , v w * . Generate a plurality of gating signals to control w * . Based on a plurality of gate signals from the PWM circuit 21, the power converter 2 controls a plurality of semiconductor switching elements inside the power converter 2 to change the terminal voltage of the synchronous motor 1 such as IPMSM or SynRM to a phase voltage. Control to command values vu * , vv * , vw * .

(2.2)軸ずれ補償器の構成
図1の軸ずれ補償器19について、第1の実施形態における構成を説明する。はじめに、γ‐δ軸とd‐q軸により、軸ずれ量θerrの演算原理を説明する。数式6は、軸ずれ発生時のd‐q軸電圧の振幅である。これにより、d‐q軸には、数式7に示す振幅を持つ基本波電流が発生する。ただし、簡単のため電機子抵抗は無視する。
(2.2) Configuration of Axial Deviation Compensator The configuration of the axial deviation compensator 19 of FIG. 1 in the first embodiment will be described. First, the calculation principle of the axis deviation amount θ err will be explained using the γ-δ axis and the dq axis. Equation 6 is the amplitude of the dq-axis voltage when the axis shift occurs. As a result, a fundamental wave current having an amplitude shown in Equation 7 is generated on the dq axes. However, for simplicity, the armature resistance is ignored.

Figure 0007163641000006
Figure 0007163641000006

Figure 0007163641000007
Figure 0007163641000007

また、γ‐δ軸電流基本波振幅は、d‐q軸電流基本波振幅と軸ずれ量θerrを用いて数式8と表現できる。数式8に数式6及び数式7を代入して整理すると、数式9及び数式10のように、γ軸電流基本波振幅およびδ軸電流基本波振幅が得られる。γ軸電圧振幅Vγを正とすると、数式9より、γ軸電流基本波振幅Iγは、軸ずれ量θerrによらず正となるが、数式10より、δ軸電流基本波振幅Iδは、d軸リアクタンスX、q軸リアクタンスXおよび軸ずれ量θerrによって極性を持つことがわかる。 Also, the γ-δ axis current fundamental wave amplitude can be expressed by Expression 8 using the d-q axis current fundamental wave amplitude and the axis deviation amount θ err . By substituting Equations 6 and 7 into Equation 8 and arranging them, the γ-axis current fundamental wave amplitude and the δ-axis current fundamental wave amplitude are obtained as shown in Equations 9 and 10. Assuming that the γ-axis voltage amplitude V γ is positive, the γ-axis current fundamental wave amplitude I γ is positive from Equation 9 regardless of the amount of axis deviation θ err . has polarities depending on the d-axis reactance X d , the q-axis reactance X q and the axis deviation θ err .

Figure 0007163641000008
Figure 0007163641000008

Figure 0007163641000009
Figure 0007163641000009

Figure 0007163641000010
Figure 0007163641000010

数式10を変形すると、数式11が得られる。ここで、2θerrが十分小さいものとして、sinθ≒θと近似すると、軸ずれ量θerrは、数式12,13で演算できる。 Transforming Equation 10, Equation 11 is obtained. Here, assuming that 2θ err is sufficiently small and approximated by sin θ≈θ, the axis deviation amount θ err can be calculated by Equations 12 and 13.

Figure 0007163641000011
Figure 0007163641000011

Figure 0007163641000012
Figure 0007163641000012

Figure 0007163641000013
Figure 0007163641000013

以上の原理に基づき、軸ずれ量θerrを補正する軸ずれ補償器19について、図2を用いて説明する。 Based on the principle described above, the axis deviation compensator 19 for correcting the axis deviation amount θ err will be described with reference to FIG. 2 .

振幅演算器22aは、数式12のδ軸電流基本波振幅Iδ(振動電流)に相当するδ軸基本波電流フーリエ正弦係数補正値Iδ(b1)'を、数式14及び数式15により演算する。これにより、前述の振幅極性を考慮した振幅を算出することが出来る。 The amplitude calculator 22a calculates a δ-axis fundamental current Fourier sine coefficient correction value I δ(b1) ' corresponding to the δ-axis current fundamental wave amplitude I δ (oscillating current) of Equation 12 using Equations 14 and 15. . As a result, the amplitude can be calculated in consideration of the aforementioned amplitude polarity.

Figure 0007163641000014
Figure 0007163641000014

Figure 0007163641000015
Figure 0007163641000015

θ演算器23aは、数式13により、δ軸電流‐軸ずれ量変換ゲインGθδを演算する。この時、d軸リアクタンスXおよびq軸リアクタンスXが未知の場合は、想定値を使用する。また、γ軸電圧振幅Vγは、数式16により演算される。 The G θ calculator 23 a calculates the δ-axis current-axis deviation amount conversion gain G θ δ by Equation (13). At this time, if the d-axis reactance Xd and the q -axis reactance Xq are unknown, assumed values are used. Also, the γ-axis voltage amplitude V γ is calculated by Equation (16).

Figure 0007163641000016
Figure 0007163641000016

上記のように演算されたIδ(b1)'とGθδとを数式12のように乗算器24aで乗算して軸ずれ補償量θerrを演算し、演算された軸ずれ補償量θerrが速度推定器25aに入力される。速度推定器25aは、例えばPI制御器で構成されており、数式17のように速度推定値ωr_estを演算する。 I δ(b1) ' and G θδ calculated as described above are multiplied by the multiplier 24a as shown in Equation 12 to calculate the compensation amount θ err for compensating for misalignment . It is input to the speed estimator 25a. The speed estimator 25a is composed of, for example, a PI controller, and calculates the speed estimated value ω r_est as shown in Equation (17).

Figure 0007163641000017
Figure 0007163641000017

積分器5bは、速度推定値ωr_estを積分し、磁極位置推定値θr_estを出力する。以上のように、軸ずれ補償器を構成することで、チューニング中に磁極位置がずれた場合であっても、チューニング電流を適切な位相に流すことが可能となる。つまり、制御装置100は、軸ずれ補償器19による補償により、d軸とγ軸が一致するようにγ,δ直交回転座標系の位相を制御できる。 The integrator 5b integrates the speed estimated value ω r_est and outputs the magnetic pole position estimated value θ r_est . By constructing the axis deviation compensator as described above, it is possible to flow the tuning current in an appropriate phase even when the magnetic pole position is deviated during tuning. In other words, the control device 100 can control the phase of the γ, δ orthogonal rotating coordinate system so that the d-axis and the γ-axis are aligned by the compensation by the axis deviation compensator 19 .

<3.第2の実施形態>
(3.1)電流及び電圧の制御方法
図1を参照して、第2の実施形態における、同期電動機1の電圧及び電流の制御方法を制御装置の構成と共に説明する。なお、電力変換器により同期電動機に供給する電圧及び電流の制御演算は、d,q直交回転座標系に一致するように制御されるγ,δ直交回転座標系の上で行うこととする。
<3. Second Embodiment>
(3.1) Current and Voltage Control Method A method for controlling the voltage and current of the synchronous motor 1 according to the second embodiment will be described together with the configuration of the control device with reference to FIG. The control calculation of the voltage and current supplied to the synchronous motor by the power converter is performed on the γ, δ orthogonal rotating coordinate system controlled to match the d, q orthogonal rotating coordinate system.

第2の実施形態の制御装置100は、γ軸電流の直流成分を指令値に制御してδ軸に正弦波の交番電圧を印加しているときの電流及び電圧に基づき、電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルの各パラメータを演算するオートチューニングを行う。あるいは、制御装置100は、γ軸電流の直流成分を指令値に制御してδ軸に正弦波の交番電圧を印加しているときの電流及び電圧に基づき、同期電動機1の磁極の方向に直交する方向のq軸インダクタンスを演算(測定)するオートチューニングを行う。制御装置100は、磁束モデルの各パラメータとq軸インダクタンスとの両方を演算するオートチューニングを行ってもよい。 The control device 100 of the second embodiment controls the DC component of the γ-axis current to a command value, and controls the magnetic saturation of the motor iron core based on the current and voltage when a sinusoidal alternating voltage is applied to the δ-axis. Auto-tuning is performed to calculate each parameter of the magnetic flux model considering at least the characteristics. Alternatively, the control device 100 controls the DC component of the γ-axis current to a command value and applies a sine-wave alternating voltage to the δ-axis, based on the current and the voltage, which is perpendicular to the direction of the magnetic poles of the synchronous motor 1 . Auto-tuning is performed to calculate (measure) the q-axis inductance in the direction of The control device 100 may perform auto-tuning to calculate both the parameters of the magnetic flux model and the q-axis inductance.

第2の実施形態において、例えば、γ軸電流直流成分指令値Iγ(0) を、0を含む任意の値に設定し、γ軸交番電流指令値の振幅Iγ(b1) を、0に設定し、δ軸交番電流指令値の振幅Iδ(b1) を、0よりも大きい任意の値に設定する。 In the second embodiment, for example, the γ-axis current DC component command value Iγ(0) * is set to an arbitrary value including 0, and the amplitude Iγ(b1) * of the γ-axis alternating current command value is set to 0, and the amplitude I δ(b1) * of the δ-axis alternating current command value is set to an arbitrary value greater than 0.

第2の実施形態では、積分器5aは、交番電圧の角周波数ωを積分してd軸交番電流指令値の角度θを演算する。δ軸交番電流指令演算器6bは、δ軸交番電流指令値の交流成分iδh を数式18のように演算する。 In the second embodiment, the integrator 5a integrates the angular frequency ωh of the alternating voltage to calculate the angle θh of the d-axis alternating current command value. The δ-axis alternating current command calculator 6b calculates the AC component i δh * of the δ-axis alternating current command value as shown in Equation (18).

Figure 0007163641000018
Figure 0007163641000018

電流座標変換器8は、u相電流検出器9u及びw相電流検出器9wによりそれぞれ検出したu,w相電流検出値i,iを、軸ずれ補償器19により演算された同期電動機1の磁極位置推定値θr_estに基づいてγ,δ軸電流検出値iγ,iδに座標変換する。 The current coordinate converter 8 converts the u- and w-phase current detection values i u and i w respectively detected by the u-phase current detector 9 u and the w-phase current detector 9 w to the synchronous motor 1 based on the magnetic pole position estimation value θ r_est of , the coordinates are transformed into γ, δ axis current detection values i γ , i δ .

ローパスフィルタ10aは、γ軸電流検出値iγの高周波成分を除去してγ軸電流検出値iγfを演算する。ローパスフィルタ10bは、δ軸電流検出値iδの高周波成分を除去してδ軸電流検出値iδfを演算する。 The low-pass filter 10a removes high-frequency components from the γ -axis current detection value to calculate the γ-axis current detection value iγf. The low-pass filter 10b removes high-frequency components from the δ -axis current detection value to calculate a δ-axis current detection value iδf.

γ軸電流指令値iγ とγ軸電流検出値iγfとの偏差を減算器11aにて演算し、この偏差をγ軸電流調節器12aにより増幅してγ軸電圧指令値vγACR(γ軸電圧フィードバック制御値vγACR)を演算する。γ軸電流調節器12aは、γ軸電流指令値iγ とγ軸電流検出値iγfとの偏差が零になるように動作してγ軸電圧フィードバック制御値vγACRを演算する。また、δ軸電流指令値iδ とδ軸電流検出値iδfとの偏差を減算器11bにて演算し、この偏差をδ軸電流調節器12bにより増幅してδ軸電圧指令値vδACR(δ軸電圧フィードバック制御値vδACR)を演算する。δ軸電流調節器12bは、δ軸電流指令値iδ とδ軸電流検出値iδfとの偏差が零になるように動作してδ軸電圧フィードバック制御値vδACRを演算する。 The deviation between the γ-axis current command value i γ * and the γ-axis current detection value i γf is calculated by the subtractor 11a, and this deviation is amplified by the γ-axis current controller 12a to obtain the γ-axis voltage command value v γACR (γ Calculate the shaft voltage feedback control value v γACR ). The γ-axis current controller 12a operates so that the deviation between the γ-axis current command value i γ * and the γ-axis current detection value i γf becomes zero to calculate the γ-axis voltage feedback control value v γACR . Further, the deviation between the δ-axis current command value i δ * and the δ-axis current detection value i δf is calculated by the subtractor 11b, and this deviation is amplified by the δ-axis current controller 12b to obtain the δ-axis voltage command value v δACR. (δ-axis voltage feedback control value v δACR ) is calculated. The δ-axis current controller 12b operates so that the deviation between the δ-axis current command value i δ * and the δ-axis current detection value i δf becomes zero to calculate the δ-axis voltage feedback control value v δACR .

磁束モデルのパラメータは、電動機鉄心の磁気飽和に起因して流れる高調波電流を利用して推定される。このため、γ軸電流調節器12a及びδ軸電流調節器12bが高調波電流に作用しないようにするため、調節器12a,12bの応答周波数、及び、ローパスフィルタ10a,10bのカットオフ周波数は、交番電圧の角周波数ωの2倍よりも小さく設定される。 The parameters of the flux model are estimated using harmonic currents that flow due to magnetic saturation of the motor core. Therefore, in order to prevent the γ-axis current regulator 12a and the δ-axis current regulator 12b from acting on harmonic currents, the response frequencies of the regulators 12a and 12b and the cutoff frequencies of the low-pass filters 10a and 10b are It is set smaller than twice the angular frequency ωh of the alternating voltage.

直流電機子抵抗補償器13は、第1の実施形態と同様に、γ軸電機子抵抗フィードフォワード補償値vγraを数式2により演算する。 The DC armature resistance compensator 13 calculates the γ-axis armature resistance feedforward compensation value v γra by Equation 2, as in the first embodiment.

電圧補償値演算器14は、インピーダンス推定器15で演算したq軸リアクタンス推定値Xqhest及びq軸電機子抵抗推定値Rqhestを使用し、δ軸電圧フィードフォワード補償値vδhFFを数式19により演算する。なお、q軸リアクタンス推定値Xqhest及びq軸電機子抵抗推定値Rqhestなどのインピーダンスの推定は、例えば特許文献1等に記載されている公知技術を用いてインピーダンス推定器15により行われることが可能である。 The voltage compensation value calculator 14 uses the q-axis reactance estimated value X qhest and the q-axis armature resistance estimated value R qhest calculated by the impedance estimator 15 to calculate the δ-axis voltage feedforward compensation value v δhFF according to Equation 19. do. Impedance estimation such as the q-axis reactance estimated value X qhest and the q-axis armature resistance estimated value R qhest may be performed by the impedance estimator 15 using the known technique described in Patent Document 1, for example. It is possible.

Figure 0007163641000019
Figure 0007163641000019

また、q軸インダクタンスの演算は、公知技術を用いてインピーダンス推定器15により行われることが可能である。例えば、インピーダンス推定器15は、q軸リアクタンス推定値Xqhestを角周波数ωで除算することによって、q軸インダクタンスLqestを演算する(Lqest=Xqhest/ω)。 Also, the calculation of the q-axis inductance can be performed by the impedance estimator 15 using known techniques. For example, the impedance estimator 15 calculates the q-axis inductance L qest by dividing the q-axis reactance estimated value X qhest by the angular frequency ω h (L qest =X qhesth ).

フーリエ係数演算器16は、θ,vδ ,iδに基づき、δ軸電圧及びγ軸電流のそれぞれのフーリエ係数を演算する。フーリエ係数演算器16は、δ軸電圧の基本波成分のフーリエ余弦係数Vδ(a1)、δ軸電圧の基本波成分のフーリエ正弦係数Vδ(b1)、δ軸電流の基本波成分のフーリエ余弦係数Iδ(a1)、δ軸電流の基本波成分のフーリエ正弦係数Iδ(b1)、δ軸電流の3倍高調波成分のフーリエ余弦係数Iδ(a3)、及び、δ軸電流の3倍高調波成分のフーリエ正弦係数Iδ(b3)を演算する。フーリエ係数の演算は、例えば特許文献1等に記載されている公知技術を用いて行われることが可能である。 The Fourier coefficient calculator 16 calculates Fourier coefficients of the δ -axis voltage and the γ-axis current based on θh, vδ * and iδ. The Fourier coefficient calculator 16 calculates the Fourier cosine coefficient V δ(a1) of the fundamental wave component of the δ-axis voltage, the Fourier sine coefficient V δ(b1) of the fundamental wave component of the δ-axis voltage, and the Fourier coefficient of the fundamental wave component of the δ-axis current. Cosine coefficient I δ(a1) , Fourier sine coefficient I δ(b1) of the fundamental wave component of the δ-axis current, Fourier cosine coefficient I δ(a3) of the triple harmonic component of the δ-axis current, and Compute the Fourier sine coefficient I δ(b3) of the triple harmonic component. Fourier coefficients can be calculated using a known technique described in Patent Document 1, for example.

フーリエ係数演算器16により演算された、γ軸電流の2倍高調波成分のフーリエ余弦係数Iγ(a2)、γ軸電流の2倍高調波成分のフーリエ正弦係数Iγ(b2)、δ軸電流の3倍高調波成分のフーリエ余弦係数Iδ(a3)、δ軸電流の3倍高調波成分のフーリエ正弦係数Iδ(b3)から、電機子抵抗補償器17は、数式20によりγ軸高調波電機子抵抗補償値vγhraおよびδ軸高調波電機子抵抗補償値vδhraを演算する。 The Fourier cosine coefficient Iγ(a2) of the double harmonic component of the γ-axis current, the Fourier sine coefficient Iγ(b2) of the double harmonic component of the γ-axis current, and the δ-axis calculated by the Fourier coefficient calculator 16 From the Fourier cosine coefficient I δ(a3) of the third harmonic component of the current and the Fourier sine coefficient I δ(b3) of the third harmonic component of the δ-axis current, the armature resistance compensator 17 calculates the γ-axis A harmonic armature resistance compensation value v γhra and a δ-axis harmonic armature resistance compensation value v δhra are calculated.

Figure 0007163641000020
Figure 0007163641000020

磁束モデル係数推定器18には、フーリエ係数演算器16により演算された、γ軸電流直流成分のフーリエ係数Iγ(0)、γ軸電流の2倍高調波成分のフーリエ余弦係数Iγ(a2)、γ軸電流の2倍高調波成分のフーリエ正弦係数Iγ(b2)、δ軸電流基本波成分のフーリエ余弦係数Iδ(a1)、δ軸電流基本波成分のフーリエ正弦係数Iδ(b1)、δ軸電流の3倍高調波成分のフーリエ余弦係数Iδ(a3)、及び、δ軸電流の3倍高調波成分のフーリエ正弦係数Iδ(b3)が入力される。磁束モデル係数推定器18は、例えば数式5のようなd軸磁束モデル及び数式12のようなq軸磁束モデルのそれぞれの係数(パラメータ)を演算する。 The magnetic flux model coefficient estimator 18 receives the Fourier coefficient I γ(0) of the γ-axis current DC component and the Fourier cosine coefficient I γ(a2 ) , Fourier sine coefficient Iγ(b2) of the double harmonic component of the γ-axis current, Fourier cosine coefficient Iδ(a1) of the fundamental wave component of the δ-axis current, Fourier sine coefficient Iδ( b1) , the Fourier cosine coefficient I δ(a3) of the triple harmonic component of the δ-axis current, and the Fourier sine coefficient I δ(b3) of the triple harmonic component of the δ-axis current are input. The magnetic flux model coefficient estimator 18 calculates respective coefficients (parameters) of the d-axis magnetic flux model such as Equation 5 and the q-axis magnetic flux model such as Equation 12, for example.

Figure 0007163641000021
Figure 0007163641000021

第2の実施形態では、磁束モデル係数推定器18は、KLd、sd、Ksdq、KLq、KsqおよびKsqdを、上記のフーリエ係数を使用して、例えば代数演算や最適化演算等の公知の演算方法により演算する。例えば、磁束モデル係数推定器18は、特許文献1等に記載されている公知技術を用いて、KLd、sd、Ksdq、KLq、KsqおよびKsqd等の磁束モデルのパラメータを演算するオートチューニングを行う。 In a second embodiment, flux model coefficient estimator 18 computes K Ld , K sd , K sdq , K Lq , K sq and K sqd using the above Fourier coefficients, e.g. Calculated by a known calculation method such as For example, the magnetic flux model coefficient estimator 18 uses the known technology described in Patent Document 1 or the like to calculate parameters of the magnetic flux model such as K Ld , K sd , K sdq , K Lq , K sq and K sqd . perform auto-tuning.

γ軸電圧指令値vγ は、加算器7b,7dにより、γ軸電圧フィードバック制御値vγACR、γ軸電機子抵抗フィードフォワード補償値vγra、γ軸高調波電機子抵抗補償値vγhraを加算して算出される。一方、δ軸電圧指令値vδ は、加算器7c,7eにより、δ軸電圧フィードバック制御値vδACR、δ軸電圧フィードフォワード補償値vδhFF、δ軸高調波電機子抵抗補償値vδhraを加算して算出される。 The γ-axis voltage command value v γ * is obtained by adding the γ-axis voltage feedback control value v γACR , the γ-axis armature resistance feedforward compensation value v γra , and the γ-axis harmonic armature resistance compensation value v γhra by adders 7b and 7d. Calculated by adding. On the other hand, the δ-axis voltage command value v δ * is obtained by adding the δ-axis voltage feedback control value v δACR , the δ-axis voltage feedforward compensation value v δhFF , and the δ-axis harmonic armature resistance compensation value v δhra by adders 7c and 7e. Calculated by adding.

上述のように演算したγ,δ軸電圧指令値vγ ,vδ は、電圧座標変換器20によって、軸ずれ補償器19により演算された磁極位置推定値θr_estに基づいて、u,v,w相の相電圧指令値v ,v ,v に変換される。 The γ, δ axis voltage command values v γ * , v δ * calculated as described above are converted by the voltage coordinate converter 20 into u, It is converted into phase voltage command values v u * , v v * , v w * of the v and w phases.

整流回路3は、三相交流電源4からの三相交流電圧を整流して直流電圧に変換し、この直流電圧をインバータ等の電力変換器2に供給する。 The rectifier circuit 3 rectifies the three-phase AC voltage from the three-phase AC power supply 4, converts it into a DC voltage, and supplies this DC voltage to the power converter 2 such as an inverter.

PWM回路21は、相電圧指令値v ,v ,v に基づいて、電力変換器2の出力電圧を相電圧指令値v ,v ,v に制御するための複数のゲート信号を生成する。電力変換器2は、PWM回路21からの複数のゲート信号に基づいて、電力変換器2内部の複数の半導体スイッチング素子を制御することにより、IPMSMやSynRMなどの同期電動機1の端子電圧を相電圧指令値v ,v ,v に制御する。 The PWM circuit 21 controls the output voltage of the power converter 2 to the phase voltage command values vu*, vv*, vw * based on the phase voltage command values vu * , vv * , vw * . Generate multiple gate signals for Based on a plurality of gate signals from the PWM circuit 21, the power converter 2 controls a plurality of semiconductor switching elements inside the power converter 2 to change the terminal voltage of the synchronous motor 1 such as IPMSM or SynRM to a phase voltage. Control to command values vu * , vv * , vw * .

(3.2)軸ずれ補償器の構成
図1の軸ずれ補償器19について、第2の実施形態における構成を説明する。はじめに、γ‐δ軸とd‐q軸により、軸ずれ量θerrの演算原理を説明する。数式22は、軸ずれ発生時のd‐q軸電圧の振幅である。これにより、d‐q軸には、数式7に示す振幅を持つ基本波電流が発生する。ただし、簡単のため電機子抵抗は無視する。
(3.2) Configuration of Shaft Offset Compensator The configuration of the shaft offset compensator 19 of FIG. 1 in the second embodiment will be described. First, the calculation principle of the axis deviation amount θ err will be explained using the γ-δ axis and the dq axis. Equation 22 is the amplitude of the dq axis voltage when the axis shift occurs. As a result, a fundamental wave current having an amplitude shown in Equation 7 is generated on the dq axes. However, for simplicity, the armature resistance is ignored.

Figure 0007163641000022
Figure 0007163641000022

また、γ‐δ軸電流基本波振幅は、d‐q軸電流基本波振幅と軸ずれ量θerrを用いて数式8と表現できる。数式8に数式22及び数式7を代入して整理すると、数式23及び数式24のように、δ軸電流基本波振幅およびγ軸電流基本波振幅が得られる。δ軸電圧振幅Vδを正とすると、数式23より、δ軸電流基本波振幅Iδは、軸ずれ量θerrによらず正となるが、数式24より、γ軸電流基本波振幅Iγは、d軸リアクタンスX、q軸リアクタンスXおよび軸ずれ量θerrによって極性を持つことがわかる。 Also, the γ-δ axis current fundamental wave amplitude can be expressed by Expression 8 using the d-q axis current fundamental wave amplitude and the axis deviation amount θ err . By substituting Equation 22 and Equation 7 into Equation 8 and arranging them, the δ-axis current fundamental wave amplitude and the γ-axis current fundamental wave amplitude are obtained as shown in Equations 23 and 24. Assuming that the δ-axis voltage amplitude V δ is positive, according to Equation 23, the δ-axis current fundamental wave amplitude I δ is positive regardless of the axis deviation amount θ err . has polarities depending on the d-axis reactance X d , the q-axis reactance X q and the axis deviation θ err .

Figure 0007163641000023
Figure 0007163641000023

Figure 0007163641000024
Figure 0007163641000024

数式24を変形すると、数式25が得られる。ここで、2θerrが十分小さいものとして、sinθ≒θと近似すると、軸ずれ量θerrは、数式26,27で演算できる。 Transformation of Equation 24 yields Equation 25. Here, assuming that 2θ err is sufficiently small and approximated by sin θ≈θ, the axis deviation amount θ err can be calculated by Equations 26 and 27.

Figure 0007163641000025
Figure 0007163641000025

Figure 0007163641000026
Figure 0007163641000026

Figure 0007163641000027
Figure 0007163641000027

以上の原理に基づき、軸ずれ量θerrを補正する軸ずれ補償器19について、図3を用いて説明する。 Based on the principle described above, the axis deviation compensator 19 for correcting the axis deviation amount θ err will be described with reference to FIG. 3 .

振幅演算器22bは、数式26のγ軸電流基本波振幅Iγ(振動電流)に相当するγ軸基本波電流フーリエ正弦係数補正値Iγ(b1)'を、数式28及び数式29により演算する。これにより、前述の振幅極性を考慮した振幅を算出することが出来る。 The amplitude calculator 22b calculates the γ-axis fundamental current Fourier sine coefficient correction value I γ(b1) ' corresponding to the γ-axis current fundamental wave amplitude I γ (oscillating current) of Equation 26 using Equations 28 and 29. . As a result, the amplitude can be calculated in consideration of the aforementioned amplitude polarity.

Figure 0007163641000028
Figure 0007163641000028

Figure 0007163641000029
Figure 0007163641000029

θ演算器23bは、数式27により、γ軸電流‐軸ずれ量変換ゲインGθγを演算する。また、δ軸電圧振幅Vδは、数式30により演算される。 The G θ calculator 23b calculates the γ-axis current-axis deviation amount conversion gain G θγ from Equation (27). Also, the δ-axis voltage amplitude V δ is calculated by Equation (30).

Figure 0007163641000030
Figure 0007163641000030

上記のように演算されたIγ(b1)'とGθγとを数式26のように乗算器24bで乗算して軸ずれ補償量θerrを演算し、演算された軸ずれ補償量θerrが速度推定器25bに入力される。速度推定器25bは、例えばPI制御器で構成されており、数式17のように速度推定値ωr_estを演算する。 I γ(b1) ' and G θγ calculated as described above are multiplied by the multiplier 24b as shown in Equation 26 to calculate the compensation amount θ err for compensating for misalignment . It is input to the speed estimator 25b. The speed estimator 25b is composed of, for example, a PI controller, and calculates the speed estimated value ω r_est as shown in Equation (17).

積分器5cは、速度推定値ωr_estを積分し、磁極位置推定値θr_estを出力する。以上のように、軸ずれ補償器を構成することで、チューニング中に磁極位置がずれた場合であっても、チューニング電流を適切な位相に流すことが可能となる。つまり、制御装置100は、軸ずれ補償器19による補償により、d軸とγ軸が一致するようにγ,δ直交回転座標系の位相を制御できる。 The integrator 5c integrates the speed estimated value ω r_est and outputs the magnetic pole position estimated value θ r_est . By constructing the axis deviation compensator as described above, it is possible to flow the tuning current in an appropriate phase even when the magnetic pole position is deviated during tuning. In other words, the control device 100 can control the phase of the γ, δ orthogonal rotating coordinate system so that the d-axis and the γ-axis are aligned by the compensation by the axis deviation compensator 19 .

図5は、制御装置が備える演算装置のハードウェア構成を例示する図である。制御装置は、電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルに基づいて演算する演算装置を備える。図5は、演算装置の一例であるマイクロコンピュータ110を示している。マイクロコンピュータ110は、メモリ121、CPU(Central Processing Unit)122、AD(Analog to Digital)変換部123、PWMモジュール124、通信部125及びタイマ126を備える。CPU122は、制御装置の制御を行うプロセッサである。通信部125は、マイクロコンピュータ110外部の上位コントローラと通信を行う。タイマ126は、タイマ値のカウントを行う。メモリ121は、プログラム等を記憶する。メモリ121内のプログラムによって、CPU122が動作する。図1の各制御ブロックの機能は、メモリ121に読み出し可能に記憶されるプログラムによってCPU122が動作することにより実現される。 FIG. 5 is a diagram illustrating a hardware configuration of an arithmetic unit included in the control device; The control device includes a computing device that performs computation based on a magnetic flux model that takes into account at least the magnetic saturation characteristics of the motor core. FIG. 5 shows a microcomputer 110, which is an example of an arithmetic device. The microcomputer 110 includes a memory 121 , a CPU (Central Processing Unit) 122 , an AD (Analog to Digital) conversion section 123 , a PWM module 124 , a communication section 125 and a timer 126 . The CPU 122 is a processor that controls the control device. The communication unit 125 communicates with a host controller outside the microcomputer 110 . The timer 126 counts timer values. The memory 121 stores programs and the like. A program in the memory 121 causes the CPU 122 to operate. The function of each control block in FIG. 1 is implemented by the CPU 122 operating according to a program stored in the memory 121 in a readable manner.

図1の各制御ブロックとは、例えば、加算器7a~7e、減算器11a,11b、交番電流指令演算器6a,6b、電流調節器12b,12c、電圧座標変換器20、電流座標変換器8、直流電機子抵抗補償器13、ローパスフィルタ10b,10c、電圧補償値演算器14、電機子抵抗補償器17、積分器5a、フーリエ係数演算器16、インピーダンス推定器15、磁束モデル係数推定器18及び軸ずれ補償器19である。 Each control block in FIG. , DC armature resistance compensator 13, low-pass filters 10b and 10c, voltage compensation value calculator 14, armature resistance compensator 17, integrator 5a, Fourier coefficient calculator 16, impedance estimator 15, magnetic flux model coefficient estimator 18 and a misalignment compensator 19 .

図1の各制御ブロックの機能は、コンピュータに各機能を実現させるプログラムによって提供可能である。また、各制御ブロックの機能は、上記のプログラムを記録したコンピュータ読み取り可能な記録媒体、又は、上記のプログラム等のコンピュータプログラムプロダクトによって提供可能である。記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、CD-ROM、磁気テープ、不揮発性のメモリカード、ROM等を用いることができる。 The function of each control block in FIG. 1 can be provided by a program that causes a computer to implement each function. Also, the function of each control block can be provided by a computer-readable recording medium recording the above program, or a computer program product such as the above program. Examples of recording media that can be used include flexible disks, hard disks, optical disks, magneto-optical disks, CD-ROMs, magnetic tapes, nonvolatile memory cards, and ROMs.

以上、同期電動機の制御装置を実施形態により説明したが、本発明は上記実施形態に限定されるものではない。他の実施形態の一部又は全部との組み合わせや置換などの種々の変形及び改良が、本発明の範囲内で可能である。 Although the control device for the synchronous motor has been described above with reference to the embodiments, the present invention is not limited to the above embodiments. Various modifications and improvements such as combination or replacement with part or all of other embodiments are possible within the scope of the present invention.

1 同期電動機
2 電力変換器
3 整流回路
4 三相交流電源
5a~5c 積分器
6a,6b 交番電流指令演算器
7a~7e 加算器
8 電流座標変換器
9u,9w 電流検出器
10a,10bローパスフィルタ
11a,11b 減算器
12a,12b 電流調節器
13 直流電機子抵抗補償器
14 電圧補償値演算器
15 インピーダンス推定器
16 フーリエ係数演算器
17 電機子抵抗補償器
18 磁束モデル係数推定器
19 軸ずれ補償器
20 電圧座標変換器
21 PWM回路
100 制御装置
110 マイクロコンピュータ
1 synchronous motor 2 power converter 3 rectifier circuit 4 three-phase AC power supply 5a-5c integrators 6a, 6b alternating current command calculators 7a-7e adder 8 current coordinate converters 9u, 9w current detectors 10a, 10b low-pass filter 11a , 11b subtractors 12a, 12b current controller 13 DC armature resistance compensator 14 voltage compensation value calculator 15 impedance estimator 16 Fourier coefficient calculator 17 armature resistance compensator 18 magnetic flux model coefficient estimator 19 axis deviation compensator 20 Voltage coordinate converter 21 PWM circuit 100 Control device 110 Microcomputer

Claims (6)

同期電動機の回転子磁極方向に平行なd軸とこのd軸に直交するq軸とにより定められるd,q直交回転座標系とは別に定義される座標系であってγ軸とδ軸とにより定められるγ,δ直交回転座標系において、電力変換器により前記同期電動機に供給する電流及び電圧を制御する制御装置であって、
電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルに基づいて演算する演算装置を備え、
前記演算装置は、
前記同期電動機のγ軸電流の直流成分を指令値に制御し、
前記γ軸に正弦波の交番電圧を印加し、
前記d軸のインダクタンスの演算と前記磁束モデルのパラメータの演算との少なくとも一方を行い、
前記同期電動機のγ軸電圧の基本波成分のフーリエ係数を演算し、
前記同期電動機のδ軸電流の基本波成分のフーリエ係数を演算し、
前記γ軸電圧の基本波成分のフーリエ係数と、前記δ軸電流の基本波成分のフーリエ係数とを用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御することを特徴とする、同期電動機の制御装置。
A coordinate system defined separately from the d, q orthogonal rotating coordinate system defined by the d-axis parallel to the rotor magnetic pole direction of the synchronous motor and the q-axis orthogonal to the d-axis, and is defined by the γ-axis and the δ-axis A control device for controlling the current and voltage supplied to the synchronous motor by a power converter in a defined γ, δ orthogonal rotating coordinate system,
A computing device that performs computation based on a magnetic flux model that takes into account at least the magnetic saturation characteristics of an iron core of the motor,
The computing device is
controlling the DC component of the γ-axis current of the synchronous motor to a command value;
applying a sinusoidal alternating voltage to the γ-axis,
performing at least one of the calculation of the d-axis inductance and the calculation of the parameters of the magnetic flux model;
calculating the Fourier coefficient of the fundamental wave component of the γ-axis voltage of the synchronous motor;
calculating the Fourier coefficient of the fundamental wave component of the δ-axis current of the synchronous motor;
Using the Fourier coefficient of the fundamental wave component of the γ-axis voltage and the Fourier coefficient of the fundamental wave component of the δ-axis current, the γ, δ orthogonal rotating coordinate system is set such that the d-axis and the γ-axis are aligned. A control device for a synchronous motor, characterized by controlling a phase.
前記演算装置は、
前記γ軸電流の基本波成分のフーリエ係数を演算し、
前記γ軸電流の基本波成分のフーリエ係数を更に用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御する、請求項1に記載の同期電動機の制御装置。
The computing device is
calculating the Fourier coefficient of the fundamental wave component of the γ-axis current;
2. The synchronous motor according to claim 1, further using a Fourier coefficient of a fundamental wave component of said γ-axis current to control the phase of said γ, δ orthogonal rotating coordinate system so that said d-axis and said γ-axis are aligned. controller.
前記演算装置は、
前記δ軸電流の基本波成分のフーリエ係数と前記γ軸電流の基本波成分のフーリエ係数とを用いて、前記δ軸に発生する振動電流を検出し、
前記γ軸電圧の基本波成分と前記振動電流の検出値とを用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御する、請求項2に記載の同期電動機の制御装置。
The computing device is
detecting an oscillating current generated in the δ-axis using the Fourier coefficient of the fundamental wave component of the δ-axis current and the Fourier coefficient of the fundamental wave component of the γ-axis current;
3. The phase of the γ, δ orthogonal rotating coordinate system is controlled so that the d-axis and the γ-axis are aligned using the fundamental wave component of the γ-axis voltage and the detected value of the oscillating current. Control device for a synchronous motor as described.
同期電動機の回転子磁極方向に平行なd軸とこのd軸に直交するq軸とにより定められるd,q直交回転座標系とは別に定義される座標系であってγ軸とδ軸とにより定められるγ,δ直交回転座標系において、電力変換器により前記同期電動機に供給する電流及び電圧を制御する制御装置であって、
電動機鉄心の磁気飽和特性を少なくとも考慮した磁束モデルに基づいて演算する演算装置を備え、
前記演算装置は、
前記同期電動機のγ軸電流の直流成分を指令値に制御し、
前記δ軸に正弦波の交番電圧を印加し、
前記q軸のインダクタンスの演算と前記磁束モデルのパラメータの演算との少なくとも一方を行い、
前記同期電動機のδ軸電圧の基本波成分のフーリエ係数を演算し、
前記同期電動機のγ軸電流の基本波成分のフーリエ係数を演算し、
前記δ軸電圧の基本波成分のフーリエ係数と、前記γ軸電流の基本波成分のフーリエ係数とを用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御することを特徴とする、同期電動機の制御装置。
A coordinate system defined separately from the d, q orthogonal rotating coordinate system defined by the d-axis parallel to the rotor magnetic pole direction of the synchronous motor and the q-axis orthogonal to the d-axis, and is defined by the γ-axis and the δ-axis A control device for controlling the current and voltage supplied to the synchronous motor by a power converter in a defined γ, δ orthogonal rotating coordinate system,
A computing device that computes based on a magnetic flux model that takes into account at least the magnetic saturation characteristics of the motor core,
The computing device is
controlling the DC component of the γ-axis current of the synchronous motor to a command value;
applying a sinusoidal alternating voltage to the δ-axis,
performing at least one of the calculation of the q-axis inductance and the calculation of the parameters of the magnetic flux model;
calculating the Fourier coefficient of the fundamental wave component of the δ-axis voltage of the synchronous motor;
calculating the Fourier coefficient of the fundamental wave component of the γ-axis current of the synchronous motor;
Using the Fourier coefficient of the fundamental wave component of the δ-axis voltage and the Fourier coefficient of the fundamental wave component of the γ-axis current, the γ, δ orthogonal rotating coordinate system is set such that the d-axis and the γ-axis are aligned. A control device for a synchronous motor, characterized by controlling a phase.
前記演算装置は、
前記同期電動機のδ軸電流の基本波成分のフーリエ係数を演算し、
前記δ軸電流の基本波成分のフーリエ係数を更に用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御する、請求項4に記載の同期電動機の制御装置。
The computing device is
calculating the Fourier coefficient of the fundamental wave component of the δ-axis current of the synchronous motor;
5. The synchronous motor according to claim 4, further using a Fourier coefficient of a fundamental wave component of said δ-axis current to control the phase of said γ, δ orthogonal rotating coordinate system so that said d-axis and said γ-axis are aligned. controller.
前記演算装置は、
前記γ軸電流の基本波成分のフーリエ係数と前記δ軸電流の基本波成分のフーリエ係数とを用いて、前記γ軸に発生する振動電流を検出し、
前記δ軸電圧の基本波成分と前記振動電流の検出値とを用いて、前記d軸と前記γ軸が一致するように前記γ,δ直交回転座標系の位相を制御する、請求項5に記載の同期電動機の制御装置。
The computing device is
detecting an oscillating current generated in the γ-axis using the Fourier coefficient of the fundamental wave component of the γ-axis current and the Fourier coefficient of the fundamental wave component of the δ-axis current;
6. The phase of the γ, δ orthogonal rotating coordinate system is controlled so that the d-axis and the γ-axis are aligned using the fundamental wave component of the δ-axis voltage and the detected value of the oscillating current. Control device for a synchronous motor as described.
JP2018128517A 2018-07-05 2018-07-05 Synchronous motor controller Active JP7163641B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018128517A JP7163641B2 (en) 2018-07-05 2018-07-05 Synchronous motor controller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018128517A JP7163641B2 (en) 2018-07-05 2018-07-05 Synchronous motor controller

Publications (2)

Publication Number Publication Date
JP2020010477A JP2020010477A (en) 2020-01-16
JP7163641B2 true JP7163641B2 (en) 2022-11-01

Family

ID=69152598

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018128517A Active JP7163641B2 (en) 2018-07-05 2018-07-05 Synchronous motor controller

Country Status (1)

Country Link
JP (1) JP7163641B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11349416B2 (en) * 2020-04-22 2022-05-31 Steering Solutions Ip Holding Corporation Feedback compensation of parameter imbalance induced current harmonics in synchronous motor drives

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009273283A (en) 2008-05-09 2009-11-19 Fuji Electric Systems Co Ltd Controller for permanent magnet type synchronous motor
JP2013066302A (en) 2011-09-16 2013-04-11 Fuji Electric Co Ltd Control device of permanent magnet type synchronous motor
JP2015144502A (en) 2014-01-31 2015-08-06 富士電機株式会社 Controller of permanent magnet synchronous motor
JP2015154620A (en) 2014-02-17 2015-08-24 富士電機株式会社 Controller of permanent magnet synchronous motor

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009273283A (en) 2008-05-09 2009-11-19 Fuji Electric Systems Co Ltd Controller for permanent magnet type synchronous motor
JP2013066302A (en) 2011-09-16 2013-04-11 Fuji Electric Co Ltd Control device of permanent magnet type synchronous motor
JP2015144502A (en) 2014-01-31 2015-08-06 富士電機株式会社 Controller of permanent magnet synchronous motor
JP2015154620A (en) 2014-02-17 2015-08-24 富士電機株式会社 Controller of permanent magnet synchronous motor

Also Published As

Publication number Publication date
JP2020010477A (en) 2020-01-16

Similar Documents

Publication Publication Date Title
JP6367332B2 (en) Inverter control device and motor drive system
KR102285041B1 (en) Inverter control unit and motor drive system
WO2016121237A1 (en) Inverter control device and motor drive system
JP2007252052A (en) Vector controller of permanent magnet motor
JPH1127999A (en) Estimating method for induced electromotive force for induction motor, speed estimating method, shaft deviation correcting method and induction motor control equipment
JP5428202B2 (en) Control device for permanent magnet type synchronous motor
JP5321792B2 (en) Control device for permanent magnet type synchronous motor
JP3951075B2 (en) Method and apparatus for controlling motor
JP2007252163A (en) Control device of voltage-type inverter
JP6248847B2 (en) Control device for permanent magnet type synchronous motor
JP2004297966A (en) Ac motor controlling device
JP2015180130A (en) Control device of permanent magnet synchronous motor
JP5104219B2 (en) Control device for permanent magnet type synchronous motor
JP5326742B2 (en) AC motor control device
JP2009273254A (en) Controller for permanent magnet type synchronous motor
JP4895120B2 (en) Control device for voltage source inverter
JP6128330B2 (en) Control device for permanent magnet type synchronous motor
JP7163641B2 (en) Synchronous motor controller
JP2000037098A (en) Power conversion apparatus using speed sensor-less vector control
JP7304891B2 (en) Rotating machine control device and electric vehicle control device
JP6971925B2 (en) Motor control device
JP2009278691A (en) Controller for permanent magnet type synchronous motor
JP2004120834A (en) Controller of dc brushless motor
JP3692085B2 (en) Motor control method and apparatus
JP6675579B2 (en) Control device for permanent magnet synchronous motor

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210614

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220406

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220412

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220517

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: 20220920

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221003

R150 Certificate of patent or registration of utility model

Ref document number: 7163641

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150