JP5644203B2 - Torque ripple suppression control device and control method for rotating electrical machine - Google Patents

Torque ripple suppression control device and control method for rotating electrical machine Download PDF

Info

Publication number
JP5644203B2
JP5644203B2 JP2010138898A JP2010138898A JP5644203B2 JP 5644203 B2 JP5644203 B2 JP 5644203B2 JP 2010138898 A JP2010138898 A JP 2010138898A JP 2010138898 A JP2010138898 A JP 2010138898A JP 5644203 B2 JP5644203 B2 JP 5644203B2
Authority
JP
Japan
Prior art keywords
order
torque
learning control
tbn
torque pulsation
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.)
Expired - Fee Related
Application number
JP2010138898A
Other languages
Japanese (ja)
Other versions
JP2012005269A (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.)
Meidensha Corp
Original Assignee
Meidensha Corp
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 Meidensha Corp filed Critical Meidensha Corp
Priority to JP2010138898A priority Critical patent/JP5644203B2/en
Publication of JP2012005269A publication Critical patent/JP2012005269A/en
Application granted granted Critical
Publication of JP5644203B2 publication Critical patent/JP5644203B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Ac Motors In General (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Description

本発明は、モータ等の回転電気機械のトルク制御装置において、回転電気機械のトルクリプル(トルク脈動)を自動的に抑制する制御装置および制御方法に係り、特に補償電流フーリエ係数学習制御による周期性外乱抑制制御に関する。   The present invention relates to a control device and a control method for automatically suppressing torque ripple (torque pulsation) of a rotary electric machine in a torque control device of a rotary electric machine such as a motor, and more particularly to periodic disturbance by compensated current Fourier coefficient learning control. It relates to suppression control.

モータは原理的にトルクリプル(トルク脈動)を発生し、振動・騒音、乗り心地への悪影響、電気・機械共振等の種々の問題を引き起こす。特に、近年普及が進んでいる埋込磁石PMモータ(IPMSM)は、コギングトルクリプルとリラクタンストルクリプルが複合的に発生する。その対策として、トルクリプルを打ち消す補償電流をトルク指令に重畳する種々の方式が検討されている。   In principle, the motor generates torque ripple (torque pulsation), causing various problems such as vibration / noise, adverse effects on riding comfort, and electrical / mechanical resonance. In particular, an embedded magnet PM motor (IPMSM), which has become popular in recent years, generates a combination of cogging torque ripple and reluctance torque ripple. As a countermeasure, various methods for superimposing a compensation current for canceling the torque ripple on the torque command have been studied.

しかしながら、例えば数式解析モデルを用いてフィードフォワード補償を行なう方式では、解析誤差の影響が懸念される。また、定常動作点でのフィードバック学習制御結果を記憶してフィードフォワード補償する方式では動作点毎に制御パラメータを適切に調整するための時間を要し、オンライン補償が難しくなる。また、電流リプルを低減する方式は、トルクリプルの観点からは最適に抑制されているとは限らない。そのほか、トルクリプルオブザーバ補償方式も検討されているが、可変速運転時の特性やオンラインフィードバック抑制の検証が不充分なものであった。   However, for example, in the method of performing feedforward compensation using a mathematical analysis model, there is a concern about the influence of analysis errors. Further, in the method of storing the feedback learning control result at the steady operating point and performing feedforward compensation, it takes time to appropriately adjust the control parameter for each operating point, and online compensation becomes difficult. Further, the method for reducing the current ripple is not always optimally suppressed from the viewpoint of torque ripple. In addition, torque ripple observer compensation methods have been studied, but the characteristics during variable speed operation and on-line feedback suppression verification have been insufficient.

上記の課題に対し、電気・機械共振の元凶であるトルクリプルを高精度に抑制する目的で、軸トルクメータによるフィードバック抑制制御法を本願発明者等は既に提案している(例えば、非特許文献1参照)。この制御法は、トルクリプルの周期性に着目して脈動周波数成分毎に補償する制御系を構築すると共に、システム同定結果を用いて動作状態変化に即応するようパラメータを自動調整する機能を設けたものであり、この制御法の詳細を以下に説明する。   The inventors of the present application have already proposed a feedback suppression control method using an axial torque meter for the purpose of highly accurately suppressing torque ripple, which is a cause of electrical / mechanical resonance, in response to the above problems (for example, Non-Patent Document 1). reference). This control method focuses on the periodicity of torque ripple and constructs a control system that compensates for each pulsation frequency component, and also has a function that automatically adjusts parameters to respond quickly to changes in the operating state using the system identification result. The details of this control method will be described below.

(1)トルクリプル抑制制御装置の基本構成
図4は、従来のトルクリプル抑制制御装置の基本構成図である。同図は、モータによって負荷をトルク加振するシステムに適用したものである。
(1) Basic Configuration of Torque Ripple Suppression Control Device FIG. 4 is a basic configuration diagram of a conventional torque ripple suppression control device. The figure is applied to a system in which a load is torque-excited by a motor.

トルクリプルの発生源となるモータ1と、何らかの負荷装置2をシャフト3で結合し、その軸トルクをトルクメータ4で計測してトルクリプル抑制装置5に入力する。また、ロータリエンコーダ等の回転位置センサ6を用いてモータの回転子位置(位相)情報を入力する。トルクリプル抑制装置(コントローラ)5は、トルク脈動抑圧手段を搭載し、トルク指令値(あるいは速度指令値)に基づいて生成された電流指令値に、トルク脈動補償電流を上乗せした指令値をインバータ7に与える。図4の例では、インバータ7で電流ベクトル制御することを考慮して、モータの回転に同期した回転座標(直交dq軸)上のd軸、q軸電流指令値id*、iq*を与えている。 A motor 1 that is a torque ripple generation source and some load device 2 are coupled by a shaft 3, and its shaft torque is measured by a torque meter 4 and input to a torque ripple suppression device 5. Moreover, the rotor position (phase) information of a motor is input using rotational position sensors 6, such as a rotary encoder. The torque ripple suppression device (controller) 5 is equipped with torque pulsation suppression means, and a command value obtained by adding a torque pulsation compensation current to the current command value generated based on the torque command value (or speed command value) is supplied to the inverter 7. give. In the example of FIG. 4, in consideration of the current vector control by the inverter 7, the d-axis and q-axis current command values id * and iq * on the rotation coordinates (orthogonal dq axes) synchronized with the rotation of the motor are given. Yes.

トルクリプル(トルク脈動)はモータの構造上、回転子位置に応じて周期的に発生することが知られている。そこで、モータ回転に同期してトルクリプル周波数成分を抽出する手段を用い、任意次数nのトルクリプルを余弦、正弦係数TAn,TBn[Nm]に変換する。トルクリプル周波数成分の厳密な計測手段にはフーリエ変換などあがるが、演算容易性を重視すれば回転位相θ[rad]を基準とした単相の高調波回転座標系に低域通過フィルタを通すことでトルクリプル周波数成分を抽出することができる。 It is known that torque ripple (torque pulsation) periodically occurs according to the rotor position due to the structure of the motor. Therefore, means for extracting a torque ripple frequency component in synchronism with the motor rotation is used to convert the torque ripple of arbitrary order n into cosine and sine coefficients T An and T Bn [Nm]. The exact measurement means for torque ripple frequency components includes Fourier transform, but if importance is placed on the ease of operation, a low-pass filter can be passed through a single-phase harmonic rotation coordinate system based on the rotation phase θ [rad]. A torque ripple frequency component can be extracted.

トルクリプル抑制装置5では、上記の余弦、正弦係数TAn,TBnを用いてトルクリプル抑制制御を行い、任意周波数成分のトルク脈動補償電流iqc*[A]の余弦/正弦係数IAn,IBn[A]を生成する。トルク脈動補償電流iqc*への変換は、変換時と同じ回転位相θを用いて求め、このトルク脈動補償電流iqc*はq軸電流指令値に重畳して通常のベクトル制御を行う。 The torque ripple suppression device 5 performs torque ripple suppression control using the cosine and sine coefficients T An and T Bn described above, and the cosine / sine coefficients I An and I Bn [of the torque pulsation compensation current iqc * [A] of an arbitrary frequency component. A] is generated. Conversion to the torque ripple compensation current iqc * is determined using the same rotational phase θ at the time of conversion, it performs a normal vector control the torque pulsation compensation current iqc * is superimposed on the q-axis current command value.

図5は、前記の軸トルクTdetと回転位相θからトルクリップルの余弦、正弦係数TAn、TBnを求め、これらと回転速度ωから補償電流余弦係数IAnおよび正弦係数IBnを求めてトルク周期性外来抑制の補償電流を求めるトルクリプル抑制装置の制御ブロック図である。図中の記号は、以下の意味である。 In FIG. 5, the cosine and sine coefficients T An and T Bn of the torque ripple are obtained from the shaft torque T det and the rotational phase θ, and the compensation current cosine coefficient I An and the sine coefficient I Bn are obtained from these and the rotational speed ω. It is a control block diagram of the torque ripple suppression apparatus which calculates | requires the compensation current of torque periodic external suppression. The symbols in the figure have the following meanings.

*:トルク指令値、Tdet:トルク検出値、TAn:n次トルク脈動抽出成分(余弦係数)、TBn:n次トルク脈動抽出成分(正弦係数)、ω:回転数検出値、θ:回転位相検出値、iqc*:トルクリプル補償電流、id:d軸電流検出値、id*:d軸電流指令値、iq:q軸電流検出値、iq*:q軸電流指令値、iu.iv.iw:u、v、w相電流、iqo*:q軸電流指令値(補償電流重畳前)、IAn:n次補償電流余弦係数、lBn:n次補償電流正弦係数、abz:回転センサ信号。なお、添え字のnはn次成分トルクリプルである。 T * : Torque command value, T det : Torque detection value, T An : n-th order torque pulsation extraction component (cosine coefficient), T Bn : n-th order torque pulsation extraction component (sine coefficient), ω: Rotation speed detection value, θ : Rotation phase detection value, iqc * : torque ripple compensation current, id: d-axis current detection value, id * : d-axis current command value, iq: q-axis current detection value, iq * : q-axis current command value, iu. iv. iw: u, v, w-phase current, iqo * : q-axis current command value (before compensation current superposition), I An : n-order compensation current cosine coefficient, l Bn : n-order compensation current sine coefficient, abz: rotation sensor signal . Note that the subscript n is an nth-order component torque ripple.

図5において、指令値変換部11は、トルク指令値T*から、ベクトル制御における回転dq座標系のd軸およびq軸電流指令値Id*、Iqo*に変換する。一般には、最大トルク/電流制御を実現するような変換数式やテーブルなどが用いられる。電流ベクトル制御部12は、q軸電流指令値iqo*にトルク脈動補償電流iqc*を重畳したものをq軸電流指令値iq*としてトルク脈動を抑制する。図5の例ではq軸電流指令値に補償電流iqcを重畳しているが、d軸電流、あるいはd軸とq軸の両方に与えても良い。あるいは、dq軸電流の干渉が問題にならないシステムであれば、トルク指令値に対して直接的にトルク脈動補償信号を重畳しても良い。 In FIG. 5, the command value conversion unit 11 converts the torque command value T * into d-axis and q-axis current command values Id * and Iqo * in the rotation dq coordinate system in vector control. In general, conversion formulas and tables that realize maximum torque / current control are used. Current vector control unit 12 suppresses the torque ripple those superimposed torque pulsation compensation current iqc * the q-axis current command value Iqo * as q-axis current command value iq *. In the example of FIG. 5, the compensation current i qc is superimposed on the q-axis current command value, but it may be applied to the d-axis current or both the d-axis and the q-axis. Alternatively, in a system where interference of dq axis current does not become a problem, a torque pulsation compensation signal may be directly superimposed on the torque command value.

電流ベクトル制御部12は、一般的な直交回転座標系d軸q軸において、電流ベクトル制御の動作を行い、モータ(IPMSM)13をベクトル制御で駆動することで、負荷装置14を駆動する。座標変換部15は、電流センサ16で検出する3相交流電流iu、iv、iwとモータ回転位相θから、モータ回転座標に同期したdq軸直交回転座標系の電流id,iqに変換する。回転位相/速度検出部17は、エンコーダ等の回転位置センサ18の回転センサ信号abzから回転速度ωおよび回転位相θの情報に変換する。   The current vector control unit 12 drives the load device 14 by performing a current vector control operation on a general orthogonal rotating coordinate system d-axis q-axis and driving the motor (IPMSM) 13 by vector control. The coordinate conversion unit 15 converts the three-phase alternating currents iu, iv, iw and the motor rotation phase θ detected by the current sensor 16 into currents id, iq in a dq axis orthogonal rotation coordinate system synchronized with the motor rotation coordinates. The rotation phase / speed detection unit 17 converts the rotation sensor signal abz of the rotation position sensor 18 such as an encoder into information on the rotation speed ω and the rotation phase θ.

トルク脈動周波数成分抽出部19では、軸トルクメータ20で検出する軸トルク検出値Tdetと回転位相θからトルク周期性外乱を脈動周波数成分ごとに抽出する。代表的な抽出手段としてフーリエ変換がある。脈動成分の抽出方法は任意であるが、演算容易性を重視し、軸トルク検出値Tdet[Nm]に回転位相θを基準としたn次余弦波・正弦波を乗じ、それぞれに低域通過フィルタを施すことで、(1)〜(3)式に示す近似的なフーリエ変換を行う。これを、トルクリプル同期座標変換と称する。 The torque pulsation frequency component extraction unit 19 extracts torque periodic disturbance for each pulsation frequency component from the shaft torque detection value T det detected by the shaft torque meter 20 and the rotation phase θ. A typical extraction means is Fourier transform. The method of extracting the pulsation component is arbitrary, but the emphasis is placed on ease of calculation, and the shaft torque detection value T det [Nm] is multiplied by the nth-order cosine wave and sine wave based on the rotational phase θ, and each is passed through a low-frequency range. By applying the filter, the approximate Fourier transform shown in the equations (1) to (3) is performed. This is called torque ripple synchronous coordinate transformation.

Figure 0005644203
Figure 0005644203

£:ラプラス変換、GF:脈動抽出フィルタ、ωf:脈動抽出ローパスフィルタ遮断周波数[rad/s]、s:ラプラス演算子
トルクリプル抑制制御部21では、(1)〜(3)式で抽出・変換したTAn、TBnを用いてトルクリプル抑制制御を各次成分で行い、補償電流iqc *[A]のn次周波数成分iqcn *[A]の余弦係数IAn *[A]ならびに正弦係数IBn *[A]を生成する。n次補償電流iqcn *への変換は、(4)式から、トルクリプル同期座標変換時と同一の回転位相θを用いて計算する。
£: Laplace transform, G F : Pulsation extraction filter, ω f : Pulsation extraction low-pass filter cut-off frequency [rad / s], s: Laplace operator The torque ripple suppression control unit 21 uses the expressions (1) to (3) Torque ripple suppression control is performed with each order component using the converted T An and T Bn, and the cosine coefficient I An * [A] and the sine of the n-th order frequency component i qcn * [A] of the compensation current i qc * [A]. A coefficient I Bn * [A] is generated. The conversion to the n-th order compensation current i qcn * is calculated from the equation (4) using the same rotation phase θ as that at the time of torque ripple synchronous coordinate conversion.

Figure 0005644203
Figure 0005644203

トルク指令T*[Nm]は、最大トルク/電流制御を実現するdq軸電流指令値id *、iq0 *に変換し、(4)式で生成した各次補償電流の合成値iqc *をiq0 *に重畳して、ベクトル制御を行う。基本的に、トルクリプル抑制制御装置で行う演算処理内容は、軸トルク脈動成分抽出・トルクリプル抑制・補償電流信号生成であり、それ以外の処理については一般的なインバータで行っている。 The torque command T * [Nm] is converted into a dq-axis current command value i d * , i q0 * that realizes maximum torque / current control, and a combined value i qc * of each primary compensation current generated by the equation (4) . Is superposed on i q0 * to perform vector control. Basically, the contents of calculation processing performed by the torque ripple suppression control device are extraction of shaft torque pulsation component, torque ripple suppression, and compensation current signal generation, and other processing is performed by a general inverter.

上記のトルクリプル抑制制御部21の代表形態として、後述の周期性外乱オブザーバ補償法または補償電流フーリエ係数学習制御法を用いることができる。   As a representative form of the torque ripple suppression control unit 21, a periodic disturbance observer compensation method or a compensation current Fourier coefficient learning control method described later can be used.

補償電流生成部22では、前記の(4)式でトルク脈動補償電流指令値iqc*を生成し、q軸電流指令値iq0 *に重畳する。なお、図5ではq軸電流指令値iq0 *にトルク脈動補償電流指令値iqc *を重畳するが、d軸電流指令値、d軸とq軸の双方の電流指令値、トルク指令値などに置き換えることも可能である。 The compensation current generator 22 generates a torque pulsation compensation current command value iqc * by the above equation (4) and superimposes it on the q-axis current command value i q0 * . In FIG. 5, the torque pulsation compensation current command value i qc * is superimposed on the q-axis current command value i q0 *. However, the d-axis current command value, the current command values for both the d-axis and the q-axis, the torque command value, etc. It is also possible to replace with.

(2)システムの同定
図4に示すようなシステム構成は、PMモータ1,負荷装置2,卜ルクメータ4,カップリング3等の慣性モーメントにより多慣性軸ねじれ共振系となる。軸トルク検出値をフィードバックする場合は複数の共振・反共振周波数があるため、動作状態に応じて適切に抑制制御パラメータを決定しなければならない。制御パラメータの学習時間か長いと電気・機械共振現象を増大させる危険があるため、速やかな自動調整機能が必要である。
(2) System Identification The system configuration shown in FIG. 4 is a multi-inertia torsional resonance system due to the moment of inertia of the PM motor 1, the load device 2, the torque meter 4, the coupling 3, and the like. When the detected value of the shaft torque is fed back, since there are a plurality of resonance / anti-resonance frequencies, the suppression control parameter must be appropriately determined according to the operation state. If the learning time of the control parameter is long, there is a risk of increasing the electric / mechanical resonance phenomenon, so a quick automatic adjustment function is necessary.

そこで、非特許文献1では、回転速度変化に適応する可変ノミナル制御パラメータを導出するために、図4のトルクリプル抑制装置5の出力値から入力値までのシステム伝達特性、すなわち図5におけるトルク脈動補償電流指令値iqc*から軸トルク検出器4の検出値Tdetまでの周波数伝達特性を同定する。システム同定手法は一般的な技術であるため、任意手法を用いることができる。 Therefore, in Non-Patent Document 1, in order to derive a variable nominal control parameter adapted to the rotational speed change, the system transfer characteristic from the output value to the input value of the torque ripple suppression device 5 in FIG. 4, that is, the torque pulsation compensation in FIG. The frequency transfer characteristic from the current command value iqc * to the detected value T det of the shaft torque detector 4 is identified. Since the system identification method is a general technique, an arbitrary method can be used.

(3)補償電流フーリエ係数学習制御法
前記の非特許文献1に方式1として記載されるトルクリプル抑制制御法であり、トルクリプルをそれがもつ周波数成分のフーリエ係数として求め、これから前記(4)式の演算でトルク脈動補償電流iqc*を求める。この制御法では、トルクリプル周波数に同期した周波数成分のシステム伝達関数を1次元複素ベクトルで表現し、任意周波数成分のトルクリプルの実部・虚部をフーリエ変換等で抽出している。その余弦・正弦フーリエ係数を複素ベクトルの実部・虚部に当てはめて、フィードバック抑制制御系を構築する。
(3) Compensation current Fourier coefficient learning control method This is a torque ripple suppression control method described as method 1 in Non-Patent Document 1 described above. The torque ripple is obtained as a Fourier coefficient of a frequency component of the torque ripple, and from the above equation (4), The torque pulsation compensation current iqc * is obtained by calculation. In this control method, a system transfer function of a frequency component synchronized with the torque ripple frequency is expressed by a one-dimensional complex vector, and a real part and an imaginary part of the torque ripple of an arbitrary frequency component are extracted by Fourier transform or the like. By applying the cosine and sine Fourier coefficients to the real and imaginary parts of the complex vector, a feedback suppression control system is constructed.

トルク脈動補償電流係数はI−P(比例・積分)学習制御方式で求める。I−P抑制制御系の閉ループ特性は、モデルマッチング法によって任意の標準系規範モデルの極配置と一致するように比例・積分ゲインを決定する。また、それらは前記のシステム同定結果と回転速度情報を用いて自動的にパラメータを適応させるので、多慣性共振系システムヘの実装を容易にする(例えば、特許文献1参照)。   The torque pulsation compensation current coefficient is obtained by the IP (proportional / integral) learning control method. The proportional / integral gain is determined by the model matching method so that the closed loop characteristic of the IP suppression control system matches the pole arrangement of an arbitrary standard system reference model. In addition, since the parameters are automatically adapted using the system identification result and the rotational speed information, the mounting on the multi-inertia resonance system is facilitated (see, for example, Patent Document 1).

只野 他、「PMモータの周期性外乱に着目したトルクリプル抑制制御法の検討」、平成21年電気学会産業応用部門大会、I−615〜618、会期:平成21年8月31日〜9月2日、会場:三重大学Kanno et al., “Examination of torque ripple suppression control method focusing on periodic disturbance of PM motor”, 2009 IEEJ Industrial Application Division Conference, I-615-618, Date: August 31-September 2, 2009 Sun, venue: Mie University

国際公開番号WO2010/024194A1International Publication Number WO2010 / 024194A1

一般に、トルクリプルを検出し、全周波数帯域に亘るフィードバック制御でトルクリプルを抑制する装置では、高周波帯域では演算無駄時間の影響を受けてフィードバック制御応答が低下し、高周波帯域の外乱抑圧性能が低下し、所望の補償電流生成が困難となる。   In general, in a device that detects torque ripple and suppresses torque ripple by feedback control over the entire frequency band, the feedback control response decreases due to the influence of computation dead time in the high frequency band, and disturbance suppression performance in the high frequency band decreases. It becomes difficult to generate a desired compensation current.

一方、前記のI−P(比例・積分)学習制御方式などで求める周期的学習法は、周期性外乱と同一周波数の正弦波・余弦波を生成してから、その正弦・余弦係数(振幅・位相と等価)を調整対象とする。したがって、この方式は、全周波数帯域の補償電流を一括して生成する制御器よりも高周波外乱への対応が容易になり、周期外乱抑圧性能の改善が期待できる。また、対象次数毎に並列に抑制制御器を構成すれば、同時に複数の次数のトルクリプル抑制にも対応できる。   On the other hand, the periodic learning method obtained by the above-described IP (proportional / integral) learning control method generates a sine wave / cosine wave having the same frequency as that of the periodic disturbance, and then generates the sine / cosine coefficient (amplitude / cosine coefficient). (Equivalent to phase). Therefore, this method is easier to deal with high frequency disturbances than a controller that collectively generates compensation currents in all frequency bands, and can be expected to improve periodic disturbance suppression performance. Further, if a suppression controller is configured in parallel for each target order, torque ripple suppression of a plurality of orders can be handled at the same time.

また、図5ではトルクリプル周波数に同期した座標で制御系を構築するため、トルクリプル補償電流指令値iqc*から軸トルク検出値Tdetまでのシステム伝達特性を同定する。このシステム同定結果から任意の周波数伝達特性のみを抽出してフィードバック制御する。定常状態において、トルクリプル周波数に同期したシステムの振幅・位相伝達特性は1次元複素ベクトルで表現できるので、実システムの周波数伝達関数を回転速度ωに関する複素ベクトルP(jω)で以下の(5)式のとおり定義することができる。なお、(5)式中のベクトル量および以降の式中では「→」記号をPの頂部に付して示すが、明細書中では矢印を省略して記載する。 Further, in FIG. 5, in order to construct a control system with coordinates synchronized with the torque ripple frequency, the system transfer characteristic from the torque ripple compensation current command value iqc * to the shaft torque detection value T det is identified. From this system identification result, only an arbitrary frequency transfer characteristic is extracted and feedback controlled. In the steady state, the amplitude / phase transfer characteristic of the system synchronized with the torque ripple frequency can be expressed by a one-dimensional complex vector. Can be defined as follows. In addition, in the vector amount in the equation (5) and in the following equations, the symbol “→” is attached to the top of P, but in the specification, the arrow is omitted.

Figure 0005644203
Figure 0005644203

A(ω):実システムの実部、PB(ω):実システムの虚部
非特許文献1では、トルクリプル同期座標で制御系を構築するため、(5)式から任意n次成分の周波数伝達特性のみを抽出し、トルクリプル同期座標の実システムを(6)式とする。つまり、任意n次成分のシステムの振幅・位相特性は、単純な1次元複素ベクトルPnで表現できることに着目する。なお、(1)〜(4)式と同様に、回転位相θを基準として実部・虚部の軸を定義し、余弦係数は実部成分、正弦係数は虚部成分に対応させている。
P A (ω): Real part of real system, P B (ω): Imaginary part of real system In Non-Patent Document 1, in order to construct a control system with torque ripple synchronous coordinates, an arbitrary n-th order component is obtained from Equation (5). Only the frequency transfer characteristic is extracted, and the actual system of torque ripple synchronous coordinates is expressed by equation (6). That is, it is noted that the amplitude / phase characteristics of an arbitrary n-order component system can be expressed by a simple one-dimensional complex vector P n . As in the equations (1) to (4), the real part / imaginary part axes are defined with reference to the rotational phase θ, the cosine coefficient corresponds to the real part component, and the sine coefficient corresponds to the imaginary part component.

Figure 0005644203
Figure 0005644203

An:実システムn次成分実部、 PBn:実システムn次成分虚部
図5の制御系におけるシステム同定結果についても、同様に1次元複素ベクトルで(7)式のとおり定義することができる。なお、(7)式中の推定値および以降の式中では「^」記号をPの頂部に付して示すが、明細書中では「P^」のように記す。
P An : Real system n-order component real part, P Bn : Real system n-order component imaginary part The system identification result in the control system of FIG. 5 can be similarly defined as a formula (7) with a one-dimensional complex vector. it can. In addition, although the estimated value in the equation (7) and the following equations are shown with the “^” symbol added to the top of P, it is written as “P ^” in the specification.

Figure 0005644203
Figure 0005644203

P^An:同定結果のn次成分実部推定値、P^Bn:同定結果のn次成分虚部推定値
例えば、1〜l000Hzまでのシステム同定結果をIHz毎に複素ベクトルで表現した場合、1000個の1次元複素ベクトルの要素からなるテーブルを構築することができる。特許文献1ではこれを実部・虚部成分に分けて、回転速度ωに関するルックアップテーブルLUTA、LUTBとする。トルクリプルn次成分と同期した(7)式の同定システムの実部P^Anおよび虚部P^Bnは、線形補間を用いてシステム同定テーブルLUTA、LUTBから参照する。
P ^ An : n-order component real part estimated value of identification result, P ^ Bn : n-order component imaginary part estimated value of identification result For example, when the system identification result from 1 to 1000 Hz is expressed by a complex vector for each IHz, A table composed of 1000 one-dimensional complex vector elements can be constructed. In Patent Document 1, this is divided into a real part and an imaginary part component to be look-up tables LUT A and LUT B relating to the rotational speed ω. The real part P ^ An and the imaginary part P ^ Bn of the identification system of the equation (7) synchronized with the torque ripple n-order component are referred to from the system identification tables LUT A and LUT B using linear interpolation.

以上から、複雑なシステムであってもシステムモデルは常に簡素な1次元複素ベクトルとなり、モータ回転速度ωの変化(トルクリプル周波数変化)に応じて(7)式も瞬時に参照値を変更すれば可変速運転にも容易に対応できる。   From the above, even in a complex system, the system model is always a simple one-dimensional complex vector, and equation (7) can be changed instantaneously according to changes in motor rotation speed ω (torque ripple frequency change). It can easily handle variable speed operation.

ところで、特許文献1では、(3)式における脈動抽出フィルタGF(s)を1次ローパスフィルタに限定して制御系を構築していた。また、閉ルーブ応答の計算にモデルマッチング手法を使用しており、学習制御パラメータはユーザが任意に指定する方法を採用していた。 By the way, in Patent Document 1, the pulsation extraction filter G F (s) in the expression (3) is limited to a primary low-pass filter to construct a control system. In addition, a model matching method is used to calculate the closed loop response, and a method in which the user arbitrarily designates the learning control parameter has been adopted.

しかしながら、実際はローパスフィルタの次数を高くすることも考えられ、その際は閉ループ応答を再考する必要がある。例えば、特許文献1における脈動抽出フィルタとして高次ローパスフィルタを使用した際、完全なモデルマッチングができず、所望の抑制効果が得られない。   However, in practice, it is possible to increase the order of the low-pass filter, and in that case, it is necessary to reconsider the closed-loop response. For example, when a high-order low-pass filter is used as the pulsation extraction filter in Patent Document 1, perfect model matching cannot be performed, and a desired suppression effect cannot be obtained.

本発明の目的は、トルク脈動抽出フィルタのフィルタ次数を高くした場合にも、学習制御パラメータを自動調整して所望の抑制効果を得ることができるトルクリプル抑制制御装置および制御方法を提供することにある。   An object of the present invention is to provide a torque ripple suppression control device and a control method capable of automatically adjusting a learning control parameter to obtain a desired suppression effect even when the filter order of a torque pulsation extraction filter is increased. .

本発明は、前記の課題を解決するため、トルク脈動抽出フィルタとして高次ローパスフィルタを使用したトルクリプル抑制の周期的学習制御に、部分的モデルマッチング手法を用いて近似的にモデルの代表的極配置を維持するようにしたもので、以下の装置および方法を特徴とする。   In order to solve the above-mentioned problem, the present invention provides a representative pole arrangement of a model approximately by using a partial model matching method for periodic learning control of torque ripple suppression using a high-order low-pass filter as a torque pulsation extraction filter. And is characterized by the following apparatus and method.

(トルクリプル抑制制御装置の発明)
(1)駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御装置において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出する構成とし、
前記学習制御系は、前記高次複素フーリエ係数TAn、TBnを使った学習制御で求めるn次補償電流IAn*,IBn*のうち回転電気機械の応答を支配する低次の極のみをマッチングさせる部分的モデルマッチング手法で前記制御パラメータを自動調整すると共に、I−PD(比例・積分・微分)学習制御方式で構成し、その3つの制御パラメータを任意のトルク脈動抽出フィルタ(次数・フィルタ形式は問わず)に対応させる手段を備えたことを特徴とする。
(Invention of torque ripple suppression control device)
(1) A torque pulsation component included in the torque detection value Tdet of the driven rotating electrical machine is detected in the form of two n-th order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficient is detected by a learning control system. The control parameters are automatically adjusted to suppress TAn and TBn to 0, and the torque ripple of the rotating electrical machine is suppressed by using the nth-order compensation currents IAn * and IBn * of the torque pulsation obtained by the learning control system. In the torque ripple suppression control device,
The torque pulsation extraction filter is a high-order low-pass filter configured to detect two high-order complex Fourier coefficients TAn and TBn,
The learning control system is a part for matching only the low-order poles that govern the response of the rotating electrical machine among the n-order compensation currents IAn * and IBn * obtained by the learning control using the high-order complex Fourier coefficients TAn and TBn. The control parameters are automatically adjusted by a dynamic model matching method, and configured by an I-PD (proportional / integral / derivative) learning control method, and the three control parameters can be selected from any torque pulsation extraction filter (regardless of order and filter type). (1) is provided .

(2)駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御装置において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出する構成とし、
前記学習制御系は、前記高次複素フーリエ係数TAn、TBnを使った学習制御で求めるn次補償電流IAn*,IBn*のうち回転電気機械の応答を支配する低次の極のみをマッチングさせる部分的モデルマッチング手法で前記制御パラメータを自動調整すると共に、
任意の標準系モデル(規範となる参照モデル)の学習制御応答周波数ωCも制御パラメータとして、モデルマッチングを行う手段を備えたことを特徴とする。
(2) Torque pulsation components included in the detected torque value Tdet of the rotating electric machine to be driven are detected in the form of two n-th order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficients are detected by a learning control system. The control parameters are automatically adjusted to suppress TAn and TBn to 0, and the torque ripple of the rotating electrical machine is suppressed by using the nth-order compensation currents IAn * and IBn * of the torque pulsation obtained by the learning control system. In the torque ripple suppression control device,
The torque pulsation extraction filter is a high-order low-pass filter configured to detect two high-order complex Fourier coefficients TAn and TBn,
The learning control system is a part for matching only the low-order poles that govern the response of the rotating electrical machine among the n-order compensation currents IAn * and IBn * obtained by the learning control using the high-order complex Fourier coefficients TAn and TBn. Automatic adjustment of the control parameters by a dynamic model matching method,
A learning control response frequency ωC of an arbitrary standard model (a reference model serving as a reference) is also used as a control parameter, and means for performing model matching is provided.

(トルクリプル抑制制御方法の発明)
(3)駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御方法において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出し、
前記学習制御系は、I−PD(比例・積分・微分)学習制御方式で構成し、その3つの制御パラメータを任意のトルク脈動抽出フィルタ(次数・フィルタ形式は問わず)に対応させることを特徴とする。
(Invention of torque ripple suppression control method)
(3) Torque pulsation components included in the detected torque value Tdet of the rotating electric machine to be driven are detected in the form of two n-order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficients are detected by a learning control system. The control parameters are automatically adjusted to suppress TAn and TBn to 0, and the torque ripple of the rotating electrical machine is suppressed by using the nth-order compensation currents IAn * and IBn * of the torque pulsation obtained by the learning control system. In the torque ripple suppression control method,
The torque pulsation extraction filter is a high-order low-pass filter and detects in the form of two high-order complex Fourier coefficients TAn and TBn,
The learning control system is configured by an I-PD (proportional / integral / derivative) learning control system, and the three control parameters are associated with an arbitrary torque pulsation extraction filter (regardless of order and filter type). And

(4)前記学習制御系は、3つの制御パラメータKp,Ki,Kdをトルク脈動抽出フィルタに対応させるために次式1−1の演算を行うことを特徴とする。

Figure 0005644203
ただし、P^AN:同定結果のn次成分実部推定値、P^BN:同定結果のn次成分虚部推定値、ωc:学習制御応答周波数[rad/s]、ωf:フーリエ変換応答周波数、 (4) The learning control system performs a calculation of the following expression 1-1 in order to make the three control parameters Kp, Ki, Kd correspond to the torque pulsation extraction filter .
Figure 0005644203
However, P ^ AN: n-order component real part estimated value of identification result, P ^ BN: n-order component imaginary part estimated value of identification result, ωc: learning control response frequency [rad / s], ωf: Fourier transform response frequency ,

(5)駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御方法において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出し、
前記学習制御系は、前記高次複素フーリエ係数TAn、TBnを使った学習制御で求めるn次補償電流IAn*,IBn*のうち回転電気機械の応答を支配する低次の極のみをマッチングさせる部分的モデルマッチング手法で前記制御パラメータを自動調整すると共に、
任意の標準系モデル(規範となる参照モデル)の学習制御応答周波数ωCも制御パラメータとして、モデルマッチングを行うことを特徴とする。
(5) Torque pulsation components included in the detected torque value Tdet of the rotating electric machine to be driven are detected in the form of two nth-order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficients are detected by a learning control system. The control parameters are automatically adjusted to suppress TAn and TBn to 0, and the torque ripple of the rotating electrical machine is suppressed by using the nth-order compensation currents IAn * and IBn * of the torque pulsation obtained by the learning control system. In the torque ripple suppression control method,
The torque pulsation extraction filter is a high-order low-pass filter and detects in the form of two high-order complex Fourier coefficients TAn and TBn,
The learning control system is a part for matching only the low-order poles that govern the response of the rotating electrical machine among the n-order compensation currents IAn * and IBn * obtained by the learning control using the high-order complex Fourier coefficients TAn and TBn. Automatic adjustment of the control parameters by a dynamic model matching method,
The model matching is performed using the learning control response frequency ωC of an arbitrary standard system model (reference model as a reference) as a control parameter.

以上のとおり、本発明によれば、トルク脈動抽出フィルタとして高次ローパスフィルタを使用したトルクリプル抑制の周期的学習制御に、部分的モデルマッチング手法を用いて近似的にモデルの代表的極配置を維持するようにしたため、トルク脈動抽出フィルタのフィルタ次数を高くした場合にも、学習制御パラメータを自動調整して所望の抑制効果を得ることができる。   As described above, according to the present invention, the representative pole arrangement of the model is approximately maintained using the partial model matching method for the periodic learning control of torque ripple suppression using the high-order low-pass filter as the torque pulsation extraction filter. Therefore, even when the filter order of the torque pulsation extraction filter is increased, the learning control parameter can be automatically adjusted to obtain a desired suppression effect.

フーリエ係数学習制御の演算ブロック図(実施形態1)。FIG. 4 is a calculation block diagram of Fourier coefficient learning control (first embodiment). フーリエ係数学習制御の演算ブロック図(実施形態2)。FIG. 4 is a calculation block diagram of Fourier coefficient learning control (second embodiment). 本発明によるトルクリプル抑制制御の実験結果の例。The example of the experimental result of the torque ripple suppression control by this invention. 従来のトルクリプル抑制制御装置の基本構成図。The basic block diagram of the conventional torque ripple suppression control apparatus. 従来のトルクリプル抑制制御装置の制御ブロック図。The control block diagram of the conventional torque ripple suppression control apparatus.

(実施形態1)
図1は、本実施形態におけるフーリエ係数学習制御の演算ブロック図である。同図において、トルクリプル周波数成分抽出部(図5の19)34は前記の(1)〜(3)式に従った演算ブロックに構成する。トルクリプル抑制制御部(図5の21)は、トルクリプル同期座標におけるI−P制御部31ならびに非干渉化部32で構成され、トルクリプル指令値の余弦係数(実部)TAn *と正弦係数(虚部)TBn *には0を入力して抑制する。33は実システムである。
(Embodiment 1)
FIG. 1 is a calculation block diagram of Fourier coefficient learning control in the present embodiment. In the figure, a torque ripple frequency component extraction unit (19 in FIG. 5) 34 is configured as a calculation block according to the above equations (1) to (3). The torque ripple suppression control unit (21 in FIG. 5) includes an IP control unit 31 and a non-interacting unit 32 in torque ripple synchronous coordinates, and a cosine coefficient (real part) T An * and a sine coefficient (imaginary part) of the torque ripple command value. Part) 0 is input to T Bn * for suppression. Reference numeral 33 denotes an actual system.

ここで、前記の(4)式で示した、n次補償電流iqcn *を複素ベクトルで表現して前記の(6)式で表現したシステムに与えた場合、トルク検出値Tdetのn次成分Tdetnの複素フーリエ係数は、オイラーの公式を用いて(8)式のように展開できる(ただし、周期外乱成分は除く)。 Here, when the n-order compensation current i qcn * shown in the above equation (4) is expressed by a complex vector and given to the system expressed by the above equation (6), the n-th order of the torque detection value T det The complex Fourier coefficient of the component T detn can be expanded as shown in equation (8) using Euler's formula (however, the periodic disturbance component is excluded).

Figure 0005644203
Figure 0005644203

正の次数(n>0)の複素フーリエ係数について、(8)式の実部成分TAnおよび虚部成分TBnをそれぞれ係数比較すると、以下の(9)式のように表すことができる。 When the real part component T An and the imaginary part component T Bn of the expression (8) are compared with each other for the positive order (n> 0) complex Fourier coefficients, they can be expressed as the following expression (9).

Figure 0005644203
Figure 0005644203

この(9)式から分かるように、トルク検出値のn次成分Tdetnの実部・虚部のそれぞれに補償電流の実部IAn *と虚部IBn *が含まれており、干渉することが分かる。この干渉を打ち消すために、(10)式に従った非干渉化部32を構成する。 As can be seen from this equation (9), it includes a real part I An * and an imaginary part I Bn * of compensation current to each of the real and imaginary parts of the n th component T DETN of the torque detection value, interferes I understand that. In order to cancel this interference, the non-interacting unit 32 according to the equation (10) is configured.

Figure 0005644203
Figure 0005644203

同定結果が正しいとして、(10)式を(9)式に代入して非干渉化すると、(11)式のように置き換えられ、実部と虚部の補償電流は独立かつ同一に制御設計できる。   If the identification result is correct and formula (10) is substituted into formula (9) to make it non-interfering, it is replaced as formula (11), and the compensation current for the real part and the imaginary part can be controlled independently and identically. .

Figure 0005644203
Figure 0005644203

従来の特許文献1では、脈動抽出フィルタGF(s)を1次ローパスフィルタに限定して制御設計していたが、本実施形態では、トルク脈動抽出フィルタGF(s)に2次以上のローパスフィルタを用いて、比例先行形I−P抑制制御系を構成する。 In conventional Patent Document 1, control design is performed by limiting the pulsation extraction filter G F (s) to the primary low-pass filter. However, in the present embodiment, the torque pulsation extraction filter G F (s) has a second or higher order. A proportional preceding type IP suppression control system is configured using a low-pass filter.

ローパスフィルタの次数や種類は任意であるが、一事例として(12)式に示す2次の二項係数形を用いた場合について説明する。   Although the order and type of the low-pass filter are arbitrary, a case where the second-order binomial coefficient form shown in the equation (12) is used as an example will be described.

Figure 0005644203
Figure 0005644203

(12)式を(11)式に代入し、図1のトルク目標値TAn *、TBn *から検出値TAn、TBnまでの閉ループ伝達関数を求めると、以下のようになる。 Substituting the equation (12) into the equation (11) and obtaining the closed-loop transfer function from the torque target values T An * and T Bn * to the detected values T An and T Bn in FIG.

Figure 0005644203
Figure 0005644203

この目標値に対する制御対象の応答を、任意の参照モデル極配置に一致させる手法がモデルマッチング法であるが、I−P制御では制御パラメータが2つとなる。一方.(13)式から分かるように、閉ループ伝達関数の特性多項式は3次式となるため、極が3つとなる。   A model matching method is a method for matching the response of the control target with respect to the target value to an arbitrary reference model pole arrangement. In the IP control, there are two control parameters. on the other hand. As can be seen from the equation (13), the characteristic polynomial of the closed-loop transfer function is a cubic equation, so there are three poles.

この場合、代表極(すなわち制御対象の応答を支配する低次の極)のみをマッチングさせる部分的モデルマッチング手法を用いて、参照モデルの応答に近似させる。例えば、ここでは置を近(14)式に示す3次の二項係数標準系参照モデルを用いて極配置を近似する。   In this case, the response of the reference model is approximated using a partial model matching method that matches only the representative pole (that is, the low-order pole that dominates the response of the controlled object). For example, here, the pole placement is approximated using a third-order binomial coefficient standard system reference model shown in the near (14) equation.

Figure 0005644203
Figure 0005644203

ωC:学習制御応答周波数[rad/s]
(13)、(14)式に関して、それぞれ分母を分子で除した多項式で係数を比較する。
ω C : Learning control response frequency [rad / s]
Regarding the equations (13) and (14), the coefficients are compared with polynomials obtained by dividing the denominator by the numerator.

Figure 0005644203
Figure 0005644203

ただし、上述したとおりI−P制御系は制御パラメータがKp、Kiの2つであるため、sに関して係数比較する次数は1次、2次までとする。係数比較の結果、(16)式が導かれる。 However, as described above, since the IP control system has two control parameters, K p and K i , the order of coefficient comparison with respect to s is limited to the first and second orders. As a result of the coefficient comparison, equation (16) is derived.

Figure 0005644203
Figure 0005644203

(16)式より、制御パラメータKp、Kiはシステム同定テーブルから読み出すP^AnおよびP^Bnによる速度適応形の可変ゲインとなることがわかる。 From the equation (16), it can be seen that the control parameters K p and K i are variable gains of a speed adaptive type by P ^ An and P ^ Bn read from the system identification table.

したがって、部分的モデルマッチング法で設計した(16)式の可変の比例・積分ゲインを用いることで、速度変化があっても常に所望の参照モデルの極配置に近似するようにパラメータを自動調整できる。これにより、様々な動作点における卜ルクリプル抑制制御の調整を省力化できる。   Therefore, by using the variable proportional / integral gain of the equation (16) designed by the partial model matching method, the parameters can be automatically adjusted so as to always approximate the pole arrangement of the desired reference model even if there is a speed change. . Thereby, the adjustment of the torque ripple suppression control at various operating points can be saved.

なお、本実施形態ではあくまで代表事例としてフィルタ次数・種類を2次の二項係数形としたが、3次以上の次数や、バターワース形などの他のフィルタ形であっても同様に設計できる。   In the present embodiment, the filter order / type is a second-order binomial coefficient form as a representative example, but other filter forms such as a third-order order or a Butterworth form can be similarly designed.

(実施形態2)
本実施形態では、実施形態1を拡張し、フーリエ係数学習制御のI−P制御部を比例微分先行形のI−PD(比例・積分・微分)制御で構成する。この構成図を図2に示し、基本的な方式の説明は実施形態1と同様であるが、微分制御要素(Kd・s)を追加することで、応答波形の改善や、実施形態1における部分的モデルマッチング法の近似の次数(精度)を高める。
(Embodiment 2)
In this embodiment, the first embodiment is expanded, and the IP control unit of the Fourier coefficient learning control is configured by proportional differential preceding type I-PD (proportional / integral / differential) control. This configuration diagram is shown in FIG. 2 and the description of the basic method is the same as that of the first embodiment. However, by adding a differential control element (Kd · s), the response waveform can be improved, or a part in the first embodiment. The degree of approximation (accuracy) of the dynamic model matching method.

例えば、実施形態1と同様に(12)式で示した2次の二項係数形を脈動抽出フィルタGF(s)に適用した場合、目標値TAn *,TBn *から検出値TAn,TBnまでの閉ループ伝達関数を求めると、以下のようになる。 For example, when the second-order binomial coefficient form expressed by the equation (12) is applied to the pulsation extraction filter G F (s) as in the first embodiment, the detected value T An is obtained from the target values T An * and T Bn *. , T Bn , the closed-loop transfer function is obtained as follows.

Figure 0005644203
Figure 0005644203

参照モデルも、実施形態1と同様に(14)式として、(14)、(17)式に関して、それぞれ分母を分子で除した多項式で係数を比較すると、以下の(18)式になる。   Similarly to the first embodiment, when the coefficients of the reference model are compared with the polynomials obtained by dividing the denominator by the numerator with respect to the expressions (14) and (17) as the expression (14), the following expression (18) is obtained.

Figure 0005644203
Figure 0005644203

I−PD制御の場合、制御パラメータはKp、Ki、Kdの3つとなるため、(18)式のsに関する全次数の係数をマッチングできるようになる。係数比較の結果、以下の(19)式を導くことができる。   In the case of I-PD control, since there are three control parameters, Kp, Ki, and Kd, it is possible to match coefficients of all orders related to s in the equation (18). As a result of the coefficient comparison, the following equation (19) can be derived.

Figure 0005644203
Figure 0005644203

(19)式より、制御パラメータKp、Ki、Kdはシステム同定テーブルから読み出すP^AnおよびP^Bnによる速度適応形の可変ゲインとなることがわかる。 From equation (19), it can be seen that the control parameters Kp, Ki, and Kd are speed adaptive variable gains by P ^ An and P ^ Bn read from the system identification table.

したがって、部分的モデルマッチング法で設計した(16)式の可変の比例・積分・微分ゲインを用いることで、速度変化があっても常に所望の参照モデルの極配置にパラメータを自動調整できる。これにより、様々な動作点におけるトルクリプル抑制制御の調整を省力化できる。   Therefore, by using the variable proportional / integral / differential gain of the equation (16) designed by the partial model matching method, the parameter can always be automatically adjusted to the pole arrangement of the desired reference model even if there is a speed change. Thereby, adjustment of torque ripple suppression control at various operating points can be saved.

また、実施形態1は部分的なモデルマッチングであったが、本実施形態の2次フィルタについては完全モデルマッチングを実現しており、近似による応答誤差をなくすことができる。3次以上の次数については、実施形態1と同様に低次から部分的にマッチングさせることで、抑制制御系の閉ループ伝達関数を所望の応答に近似させることができる。   Further, although the first embodiment is a partial model matching, the second-order filter of the present embodiment realizes a complete model matching and can eliminate a response error due to approximation. As for the third and higher orders, the closed loop transfer function of the suppression control system can be approximated to a desired response by partially matching from the lower order as in the first embodiment.

(実施形態3)
実施形態1、2における学習制御応答周波数ωCは、ユーザが所望する抑制応答特性に従って任意に決定するパラメータであった。しかし、部分的モデルマッチング法ではパラメータが不足しているため、学習制御応答周波数ωCも閉ループ伝達関数の極配置を決定づけるパラメータとして扱い、近似精度を高めることができる。本実施形態では、学習制御応答周波数ωCを積極的に用いて極配置を決定する。
(Embodiment 3)
The learning control response frequency ω C in the first and second embodiments is a parameter that is arbitrarily determined according to the suppression response characteristic desired by the user. However, since the partial model matching method has insufficient parameters, the learning control response frequency ω C can also be treated as a parameter that determines the pole arrangement of the closed-loop transfer function, and the approximation accuracy can be improved. In the present embodiment, the pole arrangement is determined by actively using the learning control response frequency ω C.

ここでは具体例として、脈動抽出フィルタGF(s)に(20)式の4次二項係数形フィルタを用いて、比例先行形I−P制御を構成する場合の計算例を示す。 Here, as a specific example, a calculation example in the case of configuring proportional precedence type IP control using a quaternary binomial coefficient filter of equation (20) for the pulsation extraction filter G F (s) is shown.

Figure 0005644203
Figure 0005644203

実施形態1、2同様に、目標値からの閉ループ伝達関数を求めると(21)式となる。   As in the first and second embodiments, when a closed loop transfer function from the target value is obtained, equation (21) is obtained.

Figure 0005644203
Figure 0005644203

部分的モデルマッチング法より、この目標値応答を(22)式に示す5次の二項係数標準系モデルの極配置に近似する。   From the partial model matching method, this target value response is approximated to the pole arrangement of the fifth-order binomial coefficient standard system model shown in Equation (22).

Figure 0005644203
Figure 0005644203

ωC:学習制御応答周波数[rad/s]
I−P制御パラメータは2変数であるが、学習応答を決定づける学習制御応答周波数ωCも制御パラメータと捉えて3変数とし、(21)、(22)式それぞれにおいて、分母を分子で除したときの特性多項式を考えて、その低次3次までを近似する。すなわち、以下の(23)式の係数比較から(24)式の比例ゲインKp、積分ゲインKiならびに学習制御応答周波数ωCを導く。
ω C : Learning control response frequency [rad / s]
The IP control parameter has two variables, but the learning control response frequency ω C that determines the learning response is also regarded as a control parameter and is set to three variables, and the denominator is divided by the numerator in each of equations (21) and (22). Approximate up to the low order 3rd order. That is, the proportional gain Kp, integral gain Ki, and learning control response frequency ω C of equation (24) are derived from the coefficient comparison of equation (23) below.

Figure 0005644203
Figure 0005644203

Figure 0005644203
Figure 0005644203

(15)式より、ゲインKp、Kiはシステム同定テーブルから読み出すP^AnおよびP^Bnによる速度適応形の可変ゲインとなる。また、実施形態1、2では任意パラメータとして扱っていた学習制御応答周波数ωCが指定され、脈動抽出フィルタカットオフ周波数ωf以外の未定パラメータがなくなる。 From the equation (15), the gains Kp and Ki are speed adaptive variable gains by P ^ An and P ^ Bn read from the system identification table. In addition, the learning control response frequency ω C treated as an arbitrary parameter in the first and second embodiments is designated, and there are no undetermined parameters other than the pulsation extraction filter cutoff frequency ω f .

本実施形態によれば、(22)式の参照モデル応答への近似精度を高めるだけでなく、未定パラメータの設計指針を決定づけることにより、所望の応答を実現することができる。   According to the present embodiment, not only the approximation accuracy to the reference model response of Equation (22) can be improved, but also the desired response can be realized by determining the design guideline of the undetermined parameter.

なお、他の次数・フィルタ形式でも同様に制御ゲインパラメータを設計できる。本実施形態ではI−P制御形で説明したが、他のPID制御形でも同様に実現可能である。   It should be noted that the control gain parameter can be similarly designed for other orders and filter types. Although the present embodiment has been described with the IP control type, other PID control types can be similarly realized.

(24)式のパラメータを適用して実験したときの結果を図3に示す。波形は上から軸トルク検出波形、6次成分トルク検出値のフーリエ係数、6次成分補償電流のフーリエ係数、トルク脈動振幅値を示している。図3では、0秒で抑制開始すると、軸トルクの共振成分は抑制され、良好な応答で脈動成分が低減していることが分かる。なお、このときの脈動抽出フィルタは4次でカットオフ周波数1Hzとしている。   FIG. 3 shows the result when the experiment was performed by applying the parameter of the equation (24). The waveform shows from the top the shaft torque detection waveform, the Fourier coefficient of the sixth-order component torque detection value, the Fourier coefficient of the sixth-order component compensation current, and the torque pulsation amplitude value. In FIG. 3, it can be seen that when the suppression starts at 0 seconds, the resonance component of the shaft torque is suppressed and the pulsation component is reduced with a good response. Note that the pulsation extraction filter at this time is quaternary and has a cutoff frequency of 1 Hz.

11 指令値変換部
12 電流ベクトル制御部
13 モータ
14 負荷装置
15 座標変換部
16 電流センサ
17 回転位相/速度検出部
18 回転位置センサ
19 トルク脈動周波数成分抽出部
20 軸トルク検出器
21 トルクリプル抑制制御部
22 補償電流生成部
31 I−P制御部
32 非干渉化部
33 実システム
34 トルクリプル抽出部
DESCRIPTION OF SYMBOLS 11 Command value conversion part 12 Current vector control part 13 Motor 14 Load apparatus 15 Coordinate conversion part 16 Current sensor 17 Rotation phase / speed detection part 18 Rotation position sensor 19 Torque pulsation frequency component extraction part 20 Shaft torque detector 21 Torque ripple suppression control part 22 Compensation current generation unit 31 IP control unit 32 Decoupling unit 33 Real system 34 Torque ripple extraction unit

Claims (5)

駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御装置において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出する構成とし、
前記学習制御系は、前記高次複素フーリエ係数TAn、TBnを使った学習制御で求めるn次補償電流IAn*,IBn*のうち回転電気機械の応答を支配する低次の極のみをマッチングさせる部分的モデルマッチング手法で前記制御パラメータを自動調整すると共に、I−PD(比例・積分・微分)学習制御方式で構成し、その3つの制御パラメータを任意のトルク脈動抽出フィルタ(次数・フィルタ形式は問わず)に対応させる手段を備えたことを特徴とする回転電気機械のトルクリプル抑制制御装置。
Torque pulsation components included in the detected torque value Tdet of the rotating electric machine to be driven are detected in the form of two n-order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficients TAn and TBn are detected by a learning control system. Control parameter of the rotating electrical machine, which automatically adjusts the control parameter so as to suppress the torque to 0 and suppresses the torque ripple of the rotating electrical machine using the n-th order compensation currents IAn * and IBn * of the torque pulsation obtained by this learning control system In the device
The torque pulsation extraction filter is a high-order low-pass filter configured to detect two high-order complex Fourier coefficients TAn and TBn,
The learning control system is a part for matching only the low-order poles that govern the response of the rotating electrical machine among the n-order compensation currents IAn * and IBn * obtained by the learning control using the high-order complex Fourier coefficients TAn and TBn. The control parameters are automatically adjusted by a dynamic model matching method, and configured by an I-PD (proportional / integral / derivative) learning control method, and the three control parameters can be selected from any torque pulsation extraction filter (regardless of order and filter type). A torque ripple suppression control device for a rotating electrical machine, characterized in that it comprises means for handling
駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御装置において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出する構成とし、
前記学習制御系は、前記高次複素フーリエ係数TAn、TBnを使った学習制御で求めるn次補償電流IAn*,IBn*のうち回転電気機械の応答を支配する低次の極のみをマッチングさせる部分的モデルマッチング手法で前記制御パラメータを自動調整すると共に、
任意の標準系モデル(規範となる参照モデル)の学習制御応答周波数ωCも制御パラメータとして、モデルマッチングを行う手段を備えたことを特徴とする回転電気機械のトルクリプル抑制制御装置。
Torque pulsation components included in the detected torque value Tdet of the rotating electric machine to be driven are detected in the form of two n-order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficients TAn and TBn are detected by a learning control system. Control parameter of the rotating electrical machine, which automatically adjusts the control parameter so as to suppress the torque to 0 and suppresses the torque ripple of the rotating electrical machine using the n-th order compensation currents IAn * and IBn * of the torque pulsation obtained by this learning control system In the device
The torque pulsation extraction filter is a high-order low-pass filter configured to detect two high-order complex Fourier coefficients TAn and TBn,
The learning control system is a part for matching only the low-order poles that govern the response of the rotating electrical machine among the n-order compensation currents IAn * and IBn * obtained by the learning control using the high-order complex Fourier coefficients TAn and TBn. Automatic adjustment of the control parameters by a dynamic model matching method,
A torque ripple suppression control device for a rotating electrical machine, comprising means for performing model matching using a learning control response frequency ωC of an arbitrary standard model (normative reference model) as a control parameter .
駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御方法において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出し、
前記学習制御系は、I−PD(比例・積分・微分)学習制御方式で構成し、その3つの制御パラメータを任意のトルク脈動抽出フィルタ(次数・フィルタ形式は問わず)に対応させることを特徴とするトルクリプル抑制制御方法。
Torque pulsation components included in the detected torque value Tdet of the rotating electric machine to be driven are detected in the form of two n-order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficients TAn and TBn are detected by a learning control system. Control parameter of the rotating electrical machine, which automatically adjusts the control parameter so as to suppress the torque to 0 and suppresses the torque ripple of the rotating electrical machine using the n-th order compensation currents IAn * and IBn * of the torque pulsation obtained by this learning control system In the method
The torque pulsation extraction filter is a high-order low-pass filter and detects in the form of two high-order complex Fourier coefficients TAn and TBn,
The learning control system is configured by an I-PD (proportional / integral / differential) learning control system, and the three control parameters are associated with an arbitrary torque pulsation extraction filter (regardless of order and filter type). torque ripple suppression control method to be.
前記学習制御系は、3つの制御パラメータKp,Ki,Kdをトルク脈動抽出フィルタに対応させるために次式1−1の演算を行うことを特徴とする請求項3に記載のトルクリプル抑制制御方法。
Figure 0005644203
ただし、P^AN:同定結果のn次成分実部推定値、P^BN:同定結果のn次成分虚部推定値、ωc:学習制御応答周波数[rad/s]、ωf:フーリエ変換応答周波数、
4. The torque ripple suppression control method according to claim 3, wherein the learning control system performs calculation of the following expression 1-1 in order to make the three control parameters Kp, Ki, and Kd correspond to the torque pulsation extraction filter .
Figure 0005644203
However, P ^ AN: n-order component real part estimated value of identification result, P ^ BN: n-order component imaginary part estimated value of identification result, ωc: learning control response frequency [rad / s], ωf: Fourier transform response frequency ,
駆動される回転電気機械のトルク検出値Tdetに含まれるトルク脈動成分をトルク脈動抽出フィルタを通して2つのn次複素フーリエ係数TAn、TBnの形で検出し、学習制御系により前記複素フーリエ係数TAn、TBnを0に抑制するように制御パラメータを自動調整し、この学習制御系で求めるトルク脈動のn次補償電流IAn*,IBn*を使って回転電気機械のトルクリプルを抑制する回転電気機械のトルクリプル抑制制御方法において、
前記トルク脈動抽出フィルタは、高次ローパスフィルタとし、2つの高次複素フーリエ係数TAn、TBnの形で検出し、
前記学習制御系は、前記高次複素フーリエ係数TAn、TBnを使った学習制御で求めるn次補償電流IAn*,IBn*のうち回転電気機械の応答を支配する低次の極のみをマッチングさせる部分的モデルマッチング手法で前記制御パラメータを自動調整すると共に、
任意の標準系モデル(規範となる参照モデル)の学習制御応答周波数ωCも制御パラメータとして、モデルマッチングを行うことを特徴とするトルクリプル抑制制御方法。
Torque pulsation components included in the detected torque value Tdet of the rotating electric machine to be driven are detected in the form of two n-order complex Fourier coefficients TAn and TBn through a torque pulsation extraction filter, and the complex Fourier coefficients TAn and TBn are detected by a learning control system. Control parameter of the rotating electrical machine, which automatically adjusts the control parameter so as to suppress the torque to 0 and suppresses the torque ripple of the rotating electrical machine using the n-th order compensation currents IAn * and IBn * of the torque pulsation obtained by this learning control system In the method
The torque pulsation extraction filter is a high-order low-pass filter and detects in the form of two high-order complex Fourier coefficients TAn and TBn,
The learning control system is a part for matching only the low-order poles that govern the response of the rotating electrical machine among the n-order compensation currents IAn * and IBn * obtained by the learning control using the high-order complex Fourier coefficients TAn and TBn. Automatic adjustment of the control parameters by a dynamic model matching method,
A torque ripple suppression control method , wherein model matching is performed using a learning control response frequency ωC of an arbitrary standard system model (reference model as a reference) as a control parameter .
JP2010138898A 2010-06-18 2010-06-18 Torque ripple suppression control device and control method for rotating electrical machine Expired - Fee Related JP5644203B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010138898A JP5644203B2 (en) 2010-06-18 2010-06-18 Torque ripple suppression control device and control method for rotating electrical machine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2010138898A JP5644203B2 (en) 2010-06-18 2010-06-18 Torque ripple suppression control device and control method for rotating electrical machine

Publications (2)

Publication Number Publication Date
JP2012005269A JP2012005269A (en) 2012-01-05
JP5644203B2 true JP5644203B2 (en) 2014-12-24

Family

ID=45536607

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010138898A Expired - Fee Related JP5644203B2 (en) 2010-06-18 2010-06-18 Torque ripple suppression control device and control method for rotating electrical machine

Country Status (1)

Country Link
JP (1) JP5644203B2 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10216891B2 (en) * 2013-11-26 2019-02-26 Fuji Corporation Method and apparatus for supporting design improvement work of mechanical apparatus
JP6695499B1 (en) 2018-06-21 2020-05-20 三菱電機株式会社 Synchronous rotating machine control device and machine learning device
JP6641051B1 (en) 2019-04-18 2020-02-05 三菱電機株式会社 Motor control device
CN113459822B (en) * 2020-03-31 2023-04-28 安徽威灵汽车部件有限公司 Method and device for suppressing jitter of electric automobile, electric automobile and storage medium
CN113459824B (en) * 2020-03-31 2023-04-28 安徽威灵汽车部件有限公司 Method and device for suppressing jitter of electric automobile, electric automobile and storage medium
CN113985731B (en) * 2021-10-09 2024-04-26 南京理工大学 Wind power test bed time delay identification method and system improved by applying first-order digital filter

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0817585B2 (en) * 1989-02-06 1996-02-21 株式会社日立製作所 Torque control device
JP2006238663A (en) * 2005-02-28 2006-09-07 Toshiba Corp Control device of a dynamo-electric motor
WO2010024194A1 (en) * 2008-08-26 2010-03-04 株式会社明電舎 Device for suppressing pulsation of electric motor

Also Published As

Publication number Publication date
JP2012005269A (en) 2012-01-05

Similar Documents

Publication Publication Date Title
JP5446988B2 (en) Torque ripple suppression control device and control method for rotating electrical machine
JP5088414B2 (en) Disturbance suppression device and disturbance suppression method for electric motor
JP5621274B2 (en) Motor torque control device
JP5637042B2 (en) Electric motor pulsation suppressing device and electric motor pulsation suppressing method
JP4910445B2 (en) IPM motor vector control device
JP5644203B2 (en) Torque ripple suppression control device and control method for rotating electrical machine
JP5929863B2 (en) Control device
JP5800108B2 (en) Periodic disturbance automatic suppression device
JP5982901B2 (en) Electric motor control device and electric motor control method
JP5877733B2 (en) Electric motor control device
JP5488043B2 (en) Motor torque control device
JP5488044B2 (en) Torque ripple suppression control device and control method for rotating electrical machine
JP2010057217A (en) Torque pulsation suppression device of electric motor, and suppression method
JP2010057218A (en) Pulsation suppression device of electric motor
JP5936770B2 (en) Rotating machine control device
JP6183554B2 (en) Periodic disturbance automatic suppression device
JP6833071B2 (en) Drive device, fluid utilization device and air conditioner
JP6089654B2 (en) Periodic disturbance automatic suppression device
JP5998663B2 (en) AC motor drive control device
JP7520747B2 (en) Inverter control device and method for manufacturing the same

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130328

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140304

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140428

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20140428

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20141020

R150 Certificate of patent or registration of utility model

Ref document number: 5644203

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees