JP2008141824A - Synchronous motor controller - Google Patents

Synchronous motor controller Download PDF

Info

Publication number
JP2008141824A
JP2008141824A JP2006323738A JP2006323738A JP2008141824A JP 2008141824 A JP2008141824 A JP 2008141824A JP 2006323738 A JP2006323738 A JP 2006323738A JP 2006323738 A JP2006323738 A JP 2006323738A JP 2008141824 A JP2008141824 A JP 2008141824A
Authority
JP
Japan
Prior art keywords
axis current
command value
value
current command
synchronous motor
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.)
Granted
Application number
JP2006323738A
Other languages
Japanese (ja)
Other versions
JP2008141824A5 (en
JP4810410B2 (en
Inventor
Yutaka Matsumoto
豊 松本
Takanori Ohashi
敬典 大橋
Masaki Sugiura
正樹 杉浦
Tomohisa Watabe
與久 渡部
Hiromichi Takano
裕理 高野
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.)
Hitachi Industrial Equipment Systems Co Ltd
Original Assignee
Hitachi Industrial Equipment Systems 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 Hitachi Industrial Equipment Systems Co Ltd filed Critical Hitachi Industrial Equipment Systems Co Ltd
Priority to JP2006323738A priority Critical patent/JP4810410B2/en
Publication of JP2008141824A publication Critical patent/JP2008141824A/en
Publication of JP2008141824A5 publication Critical patent/JP2008141824A5/ja
Application granted granted Critical
Publication of JP4810410B2 publication Critical patent/JP4810410B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Ac Motors In General (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To improve control reliability and control accuracy by making it possible to simply calculate a d-axis current command value Id<SP>*</SP>close to a theoretical value in a synchronous motor controller. <P>SOLUTION: As a current command value converter, which calculates a q-axis current command value Iq<SP>*</SP>and a d-axis current command value Id<SP>*</SP>from a torque command value T<SP>*</SP>, it is configured to calculate the q-axis current command value Iq<SP>*</SP>by multiplying the torque command value T<SP>*</SP>by a proportional constant K<SB>TI</SB>, compare this with a reference q-axis current value Iq<SB>base</SB>, calculate the d-axis current value Id<SP>*</SP>as a first linear function of the q-axis current command value Iq<SP>*</SP>, when the q-axis current command value Iq<SP>*</SP>is equal to or lower than the reference q-axis current value Iq<SB>base</SB>. When the q-axis current command value Iq<SP>*</SP>exceeds the reference q-axis current value Iq<SB>base</SB>as a result of the comparison, the converter calculates the d-axis current command value Id<SP>*</SP>as a second linear function of the q-axis current command value Iq<SP>*</SP>. <P>COPYRIGHT: (C)2008,JPO&INPIT

Description

本発明は、同期電動機制御装置に係り、特に、同期電動機の入力電流を制御する技術に関する。   The present invention relates to a synchronous motor control device, and more particularly to a technique for controlling an input current of a synchronous motor.

本発明に関連した従来技術としては、例えば、特開2000−92884号公報(特許文献1)、特開2000−358400号公報(特許文献2)及び特開平10−191700号公報(特許文献3)に記載されたものがある。特開2000−92884号公報には、サーボモータの電流位相角を0°以外の一定値に制御し、リラクタンストルクを考慮し精度の良い定トルク及び定出力駆動特性を有する制御装置を得るために、逆突極性を有するDCブラシレスモータ駆動において、トルク分電流Iqから見た位相角が一定値となるようにモータ1次電流を制御し、また、q軸インダクタンスLqの変動を補償することで、精度の良い定トルク及び定出力駆動特性を有するようにトルク分電流指令値Iq及び励磁電流指令値Idを算出するとした技術が記載され、励磁電流指令値Idをモータの基底速度ωによる2つの制御範囲に分けて定めるとしている。また、特開2000−358400号公報には、入力されたトルク指令値Tqと所定値Trefとの大小関係に応じてq軸電流指令値Iq及びd軸電流指令値Idの算出方式を切り換える、すなわちトルク指令値Tqが所定値Trefよりも小さな場合には、q軸電流指令値Iqを、Iq=Kq1×Tqにより演算し(Kq1は係数)、d軸電流指令値Idを、Id=0として演算し、トルク指令値Tqが所定値Trefよりも大きい場合には、Iq=Kq1×Trefにより演算し、d軸電流指令値Idを、Id=Kd1×(Tq−Tref)により演算する(Kd1は係数)ことで、出力トルクとトルク指令値との関係を線形化し、常に高精度のトルク制御を行えるようにするとした技術が記載され、特開平10−191700号公報には、回転子内に永久磁石を有する三相同期電動機において、定出力範囲を広くとれる特性を損なうことなく、低速回転領域においての高速応答を満たし、かつ出力トルクの低下を防ぐ制御方法を提供するために、トルク指令Tと速度フィードバック信号Nfbとに基づいてd軸電流指令値Id及びq軸電流指令Iqを設定し、運転モード切替信号が与えられたとき、q軸電流指令IqをK1Iq(K1≧1)に、d軸電流指令Idを一定値に、それぞれ切り替えるとした技術が記載されている。 As conventional techniques related to the present invention, for example, Japanese Patent Application Laid-Open No. 2000-9284 (Patent Document 1), Japanese Patent Application Laid-Open No. 2000-358400 (Patent Document 2) and Japanese Patent Application Laid-Open No. 10-191700 (Patent Document 3). There is what is described in. Japanese Patent Laid-Open No. 2000-92984 discloses a control device that controls the current phase angle of a servo motor to a constant value other than 0 °, and takes into account reluctance torque and has high accuracy constant torque and constant output drive characteristics. In the DC brushless motor drive having the reverse saliency, the motor primary current is controlled so that the phase angle viewed from the torque current Iq becomes a constant value, and the fluctuation of the q-axis inductance Lq is compensated, describes a technique to calculate a torque current command value Iq * and the excitation current command value Id * to have a good constant torque and constant output drive characteristics precision, base speed of the excitation current command value Id * motor omega b The two control ranges are determined separately. Japanese Patent Laid-Open No. 2000-358400 discloses a method of calculating the q-axis current command value Iq * and the d-axis current command value Id * according to the magnitude relationship between the input torque command value Tq * and the predetermined value Tref. When the torque command value Tq * is smaller than the predetermined value Tref, the q-axis current command value Iq * is calculated by Iq * = Kq1 × Tq * (Kq1 is a coefficient), and the d-axis current command value Id * Is calculated as Id * = 0, and when the torque command value Tq * is larger than the predetermined value Tref, Iq * = Kq1 × Tref is calculated, and the d-axis current command value Id * is calculated as Id * = Kd1 × computed by (Tq * -Tref) (Kd1 coefficient) that is to linearize the relationship between the output torque and the torque command value, always technique wherein was to allow the torque control with high precision In JP-A-10-191700, in a three-phase synchronous motor having a permanent magnet in a rotor, a high-speed response in a low-speed rotation region is satisfied without impairing characteristics that can widen a constant output range, and an output is provided. In order to provide a control method for preventing a decrease in torque, the d-axis current command value Id * and the q-axis current command Iq * are set based on the torque command T * and the speed feedback signal N fb, and the operation mode switching signal is A technique is described in which the q-axis current command Iq * is switched to K1Iq * (K1 ≧ 1) and the d-axis current command Id * is switched to a constant value when given.

特開2000−92884号公報Japanese Patent Laid-Open No. 2000-92984 特開2000−358400号公報JP 2000-358400 A 特開平10−191700号公報JP 10-191700 A

例えば、上記特開2000−92884号公報では、一定の電流位相角をβとするとき、励磁電流指令値(d軸電流指令値)Idを次の数1により演算している。
Id=tan−1β×Iq …(数1)
同一電流にて最大トルクを発生させる電流位相角βは、トルク式をβで編微分することにより、次の数2で与えられる。ただし、Фは永久磁石による鎖交磁束、Lqはq軸インダクタンス、Ldはd軸インダクタンス、Ia*は電流指令実行値である。
β=sin-1[〔-Фa+{Фa+8(Lq-Ld)Ia*2}〕/4(Lq-Ld)Ia]…(数2)
上記数2に示すように、βはIa*によって変化する。上記特開2000−92884号公報では、電流値Ia*においてモータが最大効率となる角度をβとし、励磁電流指令値Idを演算している。このように様々な近似を用いたId演算方法が考えられる。しかしながら、最近は、演算にマイコンを用いるのが一般的であり、マイコンの演算能力によって、得られる結果も異なる。このため、マイコンの能力を考慮する必要がある。
For example, in the above Japanese Patent Laid-Open No. 2000-92984, the excitation current command value (d-axis current command value) Id * is calculated by the following equation 1 when β is the constant current phase angle.
Id * = tan −1 β × Iq * (Equation 1)
The current phase angle β for generating the maximum torque with the same current is given by the following equation 2 by differentiating the torque equation with β. Here, Ф is a flux linkage by a permanent magnet, Lq is a q-axis inductance, Ld is a d-axis inductance, and Ia * is a current command execution value.
β = sin -1 [[-Фa + {Фa 2 +8 (Lq -Ld) Ia * 2} ] / 4 (Lq-Ld) Ia *] ... ( Equation 2)
As shown in Equation 2 above, β varies with Ia * . In the above Japanese Patent Laid-Open No. 2000-92984, the excitation current command value Id * is calculated with β being the angle at which the motor has the maximum efficiency at the current value Ia * . Thus, Id * calculation methods using various approximations are conceivable. However, recently, it is common to use a microcomputer for calculation, and the results obtained differ depending on the calculation capability of the microcomputer. For this reason, it is necessary to consider the capability of the microcomputer.

本発明の課題点は、上記従来技術の状況に鑑み、同期電動機制御装置において、
従来よりも簡単に、理論値に近いd軸電流指令値Idを算出できるようにし、制御の信頼性と制御精度とを改善できるようにすることである。
The problem of the present invention is that in the synchronous motor control device in view of the state of the prior art,
It is to make it possible to calculate the d-axis current command value Id * close to the theoretical value more easily than in the past, and to improve the control reliability and control accuracy.

上記課題点を解決するために、本発明の同期電動機制御装置では、トルク指令値Tに基づきq軸電流指令値Iqとd軸電流指令値Idとを演算する電流指令値変換器として、トルク指令値Tに比例定数KTIを乗算して上記q軸電流指令値Iqを演算し、これを基準のq軸電流値Iqbaseと比較し、該q軸電流指令値Iqが基準のq軸電流値Iqbase以下の場合は、上記d軸電流指令値Idを該q軸電流指令値Iqの第1の1次関数として演算し、一方、上記比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbaseを超える場合は、d軸電流指令値Idを該q軸電流指令値Iqの第2の1次関数として演算する構成を備えたものとする。上記第1の1次関数は、上記q軸電流指令値Iqに比例定数aを乗じたものであり、上記第2の1次関数は、q軸電流指令値Iqに比例定数aを乗じさらに定数bを加算したものである。かかる構成とすることで、同期電動機に対し、同一電流値において理論値に近い最大トルクを出力させる。 In order to solve the above problems, the synchronous motor control device of the present invention is a current command value converter that calculates a q-axis current command value Iq * and a d-axis current command value Id * based on a torque command value T *. , is multiplied by a proportionality constant K TI to the torque command value T * is calculated the q-axis current command value Iq *, which is compared with a reference of the q-axis current value Iq base, the q-axis current command value Iq * is If the reference q-axis current value Iq base or less, the d-axis current command value Id * is calculated as a first linear function of the q-axis current command value Iq * . If the axis current command value Iq * exceeds the q-axis current value Iq base of the reference, a configuration that calculates the d-axis current command value Id * as the second primary function of the q-axis current command value Iq * Shall be. The first primary function is one obtained by multiplying the proportional constant a 1 on the q-axis current command value Iq *, the second primary function, the q-axis current command value Iq * to the proportionality constant a 2 And a constant b is added. With this configuration, the synchronous motor is caused to output a maximum torque close to the theoretical value at the same current value.

また、電流指令値変換器として、上記トルク指令値Tに比例定数KTIを乗算して上記q軸電流指令値Iqを演算し、これを基準のq軸電流値Iqbaseと比較し、該比較の結果、q軸電流指令値Iqが基準のq軸電流値Iqbase以下の場合は、第1のd軸電流指令値Id を該q軸電流指令値Iqの第1の1次関数として演算し、一方、q軸電流指令値Iqが基準のq軸電流値Iqbaseを超える場合は、第2のd軸電流指令値Id を該q軸電流指令値Iqの第2の1次関数として演算し、さらに、交流変換される直流電圧Vdcから基準の電気角速度ωbaseを演算し、当該同期電動機の電気角速度ωと比較し、該電気角速度ωが該基準の電気角速度ωbase以上の場合に、該電気角速度ωに比例するd軸電流指令増加値Idωを演算し、上記第1のd軸電流指令値Id または上記第2のd軸電流指令値Id に上記d軸電流指令増加値Idωを加算して、該加算結果をd軸電流指令値Idとして出力する構成のものとする。かかる構成とすることで、同期電動機に対し、重負荷時、最大トルクまたはこれに近いトルクで高速回転させ、電圧飽和が生じる場合にも、安定したトルクを出力させる。 Further, as the current command value converter multiplies the proportional constant K TI to the torque command value T * is calculated the q-axis current command value Iq *, which is compared with a reference of the q-axis current value Iq base, As a result of the comparison, when the q-axis current command value Iq * is less than or equal to the reference q-axis current value Iq base , the first d-axis current command value Id * 1 is set to the first q-axis current command value Iq * . On the other hand, when the q-axis current command value Iq * exceeds the reference q-axis current value Iq base , the second d-axis current command value Id * 2 is calculated as the q-axis current command value Iq *. Is calculated from the DC voltage Vdc that is AC-converted, and a reference electrical angular velocity ω base is calculated and compared with the electrical angular velocity ω of the synchronous motor. D-axis current command increase value Id proportional to the electrical angular velocity ω when the electrical angular velocity ω base or more ω is calculated, the d-axis current command increase value Id ω is added to the first d-axis current command value Id * 1 or the second d-axis current command value Id * 2 , and the addition result is d The shaft current command value Id * is output. With this configuration, the synchronous motor is rotated at a high speed with a maximum torque or a torque close to the maximum torque under heavy load, and a stable torque is output even when voltage saturation occurs.

本発明によれば、同期電動機を安定した高トルクで運転可能な同期電動機制御装置を提供することができる。   ADVANTAGE OF THE INVENTION According to this invention, the synchronous motor control apparatus which can drive a synchronous motor with the stable high torque can be provided.

以下、本発明の実施例につき、図面を用いて説明する。
図1〜図3は、本発明の同期電動機制御装置の第1の実施例の説明図である。図1は、第1の実施例としての同期電動機制御装置の構成例図、図2は、図1の同期電動機制御装置を構成する電流指令値変換器の構成例図、図3は、図1の同期電動機制御装置におけるq軸電流指令値Iqとd軸電流指令値Idとの関係の説明図である。
Embodiments of the present invention will be described below with reference to the drawings.
1 to 3 are explanatory views of a first embodiment of the synchronous motor control device of the present invention. 1 is a configuration example diagram of a synchronous motor control device as a first embodiment, FIG. 2 is a configuration example diagram of a current command value converter constituting the synchronous motor control device of FIG. 1, and FIG. It is explanatory drawing of the relationship between the q-axis current command value Iq * and the d-axis current command value Id * in the synchronous motor control apparatus.

図1において、1は同期電動機制御装置、100は同期電動機、11は、同期電動機100のロータ(図示なし)の回転位置を検出する位置検出器としてのエンコーダ、12は、上記エンコーダ11の出力から上記ロータの回転速度を演算する速度演算器、13は、上記速度演算器12から得られる上記ロータの現在の回転速度Nと、予め設定してある速度指令値Nとから、トルク指令値Tを演算し該トルク指令値Tに対応する信号(便宜上、トルク指令値Tという)を出力する速度制御演算器、14は、上記トルク指令値Tに基づき、q軸電流指令値Iqとd軸電流指令値Idとを演算し、q軸電流指令値Iq、d軸電流指令値Idのそれぞれに対応する信号(便宜上、q軸電流指令値Iq、d軸電流指令値Idという)を形成して出力する電流指令値変換器、15は、同期電動機100の入力電流値を検出する電流検出器、16は、該検出した入力電流値をq軸電流値Iq及びd軸電流値Idに変換する軸電流値変換器としてのU−V−W/d−q変換器、17aは、上記電流指令値変換器14から出力されたq軸電流指令値Iqと、上記U−V−W/d−q変換器16から出力されたq軸電流値Iqと加え合わせ、該q軸電流指令値Iqから該q軸電流値Iqを減算する加算部、17bは、電流指令値変換器14から出力されたd軸電流指令値Idと、U−V−W/d−q変換器16から出力されたd軸電流値Idと加え合わせ、該d軸電流指令値Idから該d軸電流値Idを減算する加算部、17cは、予め設定されている同期電動機100のロータの速度指令値Nに対応する信号(便宜上、速度指令値Nという)と、上記速度演算器12から出力された同ロータの現在速度Nに対応する信号(便宜上、現在速度Nという)とを加え合わせ、該速度指令値Nから該現在速度Nを減算し速度偏差情報を出力する加算部、18は、上記q軸電流指令値Iqから上記q軸電流値Iqを減算した電流値偏差情報と、上記d軸電流指令値Idから上記d軸電流値Idを減算した電流値偏差情報とに基づき、上記同期電動機100の入力電流値を演算し出力する電流制御演算器である。該電流制御演算器18から出力される制御された該入力電流は、同期電動機100の固定子コイル(図示なし)に入力され、該同期電動機100が駆動される。上記速度制御演算器13は、上記加算部17cから出力される速度偏差情報に基づいて上記トルク指令値Tを形成して出力する。上記電流検出器15及び上記U−V−W/d−q変換器16は、上記電流制御演算器18に対してフィードバック制御系を形成する。 In FIG. 1, 1 is a synchronous motor control device, 100 is a synchronous motor, 11 is an encoder as a position detector that detects the rotational position of a rotor (not shown) of the synchronous motor 100, and 12 is an output of the encoder 11. A speed calculator 13 for calculating the rotational speed of the rotor, a torque command value T from the current rotational speed N of the rotor obtained from the speed calculator 12 and a preset speed command value N *. * the calculated signal corresponding to the torque command value T * (for convenience, referred to the torque command value T *) the speed control calculator which outputs a, 14, based on the torque command value T *, q-axis current command value Iq * And the d-axis current command value Id * are calculated and signals corresponding to the q-axis current command value Iq * and the d-axis current command value Id * (for convenience, the q-axis current command value Iq * and the d-axis current command value current command value converter d * hereinafter) formed by the output, 15, synchronous current detector for detecting the input current value of the motor 100, 16, an input current value the detected q-axis current value Iq and d A UVW / dq converter, 17a serving as an axial current value converter for converting into an axial current value Id, is a q-axis current command value Iq * output from the current command value converter 14; In addition to the q-axis current value Iq output from the UVV / dq converter 16 and adding the q-axis current value Iq to the q-axis current command value Iq *, an adder 17b The d-axis current command value Id * output from the command value converter 14 and the d-axis current value Id output from the UVW / dq converter 16 are added to the d-axis current command value Id. An adder 17c for subtracting the d-axis current value Id from * A signal corresponding to the speed command value N * of the rotor of the electric motor 100 (for convenience, the speed command value N * ) and a signal corresponding to the current speed N of the rotor output from the speed calculator 12 (for convenience, the current speed) N), and an addition unit 18 that subtracts the current speed N from the speed command value N * and outputs speed deviation information, and 18 adds the q-axis current value Iq from the q-axis current command value Iq *. Based on the subtracted current value deviation information and the current value deviation information obtained by subtracting the d-axis current value Id from the d-axis current command value Id *, a current control calculation for calculating and outputting the input current value of the synchronous motor 100 It is a vessel. The controlled input current output from the current control arithmetic unit 18 is input to a stator coil (not shown) of the synchronous motor 100, and the synchronous motor 100 is driven. The speed control calculator 13 forms and outputs the torque command value T * based on the speed deviation information output from the adder 17c. The current detector 15 and the UVW / dq converter 16 form a feedback control system for the current control calculator 18.

なお、図1の構成では、U−V−W/d−q変換器16に3相の電流検出器15から検出した電流値を入力しているが、例えば,これを2相とし、U相、W相の電流を検出して、V相電流Ivを、U相電流IuとW相電流Iwとの差とし、Iv=Iu−Iwと演算してもよい。
以下の説明中で用いる図1の構成要素には、図1の場合と同じ符号を付して用いる。
In the configuration of FIG. 1, the current value detected from the three-phase current detector 15 is input to the UVV / dq converter 16. The W-phase current may be detected, and the V-phase current Iv may be calculated as Iv = Iu−Iw with the difference between the U-phase current Iu and the W-phase current Iw.
1 used in the following description are given the same reference numerals as in FIG.

図2は、図1の同期電動機制御装置1を構成する電流指令値変換器14の構成例図である。
図2において、141は、トルク指令値Tに比例定数KTIを乗算し上記q軸電流指令値Iqを演算して出力するq軸電流指令値演算器、142は、該出力されたq軸電流指令値Iqのレベルを基準のq軸電流値Iqbaseのレベルと比較する比較部、143は、該比較部142での比較の結果、上記q軸電流指令値Iqのレベルが上記基準のq軸電流値Iqbaseのレベル以下の場合は、上記d軸電流指令値Idを該q軸電流指令値Iqの第1の1次関数として演算する第1の1次関数演算器、143aは、上記q軸電流指令値Iqに乗ずる比例定数aを設定する比例定数設定部、144は、上記比較部142での比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbaseを超える場合、上記d軸電流指令値Idを該q軸電流指令値Iqの第2の1次関数として演算する第2の1次関数演算器、144aは、上記q軸電流指令値Iqに乗ずる比例定数aを設定する比例定数設定部、144bは、該q軸電流指令値Iqと比例定数aとの積に加算される定数bを設定する定数設定部である。上記比較部142は、上記比較の結果、上記q軸電流指令値Iqのレベルが上記基準のq軸電流値Iqbaseのレベル以下の場合は、上記q軸電流指令値Iqが上記第1の1次関数演算器143に入力されるようにし、一方、上記比較の結果、上記q軸電流指令値Iqのレベルが上記基準のq軸電流値Iqbaseのレベルを超える場合は、上記q軸電流指令値Iqが上記第2の1次関数演算器144に入力されるようにする。q軸電流指令値Iqが第1の1次関数演算器143に入力された場合には、該第1の1次関数演算器143において、該q軸電流指令値Iqに比例定数aが乗算され、積がd軸電流指令値Idとして出力される。また、q軸電流指令値Iqが第2の1次関数演算器144に入力された場合には、該第2の1次関数演算器144において、該q軸電流指令値Iqに比例定数aが乗算され、さらに該積に定数bが加算されたものがd軸電流指令値Idとして出力される。
FIG. 2 is a configuration example diagram of the current command value converter 14 constituting the synchronous motor control device 1 of FIG.
2, 141, by multiplying the proportionality constant K TI to the torque command value T * q-axis current command value calculator that calculates and outputs the q-axis current command value Iq *, 142 was the output q The comparison unit 143 compares the level of the axis current command value Iq * with the level of the reference q-axis current value Iq base. As a result of the comparison by the comparison unit 142, the level of the q-axis current command value Iq * is A first linear function calculator that calculates the d-axis current command value Id * as a first linear function of the q-axis current command value Iq * when the level is equal to or lower than the reference q-axis current value Iq base . , 143a are proportional constant setting unit for setting a proportional constants a 1 multiplied to the q-axis current command value Iq *, 144 as a result of the comparison in the comparison unit 142, the q-axis current command value Iq * is the reference When the q-axis current value Iq base exceeds the d-axis current command Second primary function calculator for calculating the value Id * as the second primary function of the q-axis current command value Iq *, 144a may set the proportionality constant a 2 is multiplied to the q-axis current command value Iq * proportionality constant setting unit which, 144b is a constant setting unit for setting a constant b to be added to the product of the q-axis current command value Iq * and the proportionality constant a 2. The comparing unit 142, the result of the comparison, the q-axis current when the command value Iq * level is below level of q-axis current value Iq base of the criteria, the q-axis current command value Iq * is the first When the level of the q-axis current command value Iq * exceeds the level of the reference q-axis current value Iq base as a result of the comparison, the q The shaft current command value Iq * is input to the second linear function calculator 144. When the q-axis current command value Iq * is input to the first linear function calculator 143, the first linear function calculator 143 uses the proportional constant a 1 to the q-axis current command value Iq *. And the product is output as the d-axis current command value Id * . When the q-axis current command value Iq * is input to the second linear function calculator 144, the second linear function calculator 144 uses the proportional constant to the q-axis current command value Iq *. Multiplying a 2 and adding a constant b to the product is output as a d-axis current command value Id * .

電流指令値変換器14にて使用される上記比例定数a、a2、基準のq軸電流値Iqbaseの演算・算出方法、及び、上記第1の1次関数演算器143、上記第2の1次関数演算器144の処理内容につき、以下に説明する。
同一電流に対して最大トルクが発生される時のd軸電流指令値Idとq軸電流指令値Iqの関係は、一般に、上記数1、数2と、下記数3、数4とから導かれる数5で示される。ただし、同期電動機100のq軸インダクタンスLq、d軸インダクタンスLdは、Lq>Ldであるとする。
Iq=−Ia×sinβ …(数3)
Id=Ia×cosβ …(数4)
Id=Фa/{2(Lq-Ld)}
−〔Фa/{4(Lq-Ld)}+Iq*21/2 …(数5)
同一電流において、最大トルクを発生させるためのq軸電流指令値Iqは、厳密には上記数5から演算されるが、該演算処理に時間がかかるため、該数5を以下のように近似して、q軸電流指令値Iqを算出する。すなわち、
数5において、
Фa/{2(Lq-Ld)}=A …(数6)
とし、上記数5と数6とから、O点(Iq, Id)=(0,0)、X点(Iq, Id)=(Iqbase, Idbase)=(A,(1−21/2)A)、Y点(Iq, Id)=(Iqmax, Idmax)=(Iqmax,A−(A+Iqmax 1/2)とする。以下に示すようにIqを2つの範囲に分けて、上記O点−X点間とX点−Y点間の直線を示す1次関数から算出する。
The proportional constants a 1 and a 2 used in the current command value converter 14, the calculation and calculation method of the reference q-axis current value Iq base , the first linear function calculator 143, and the second The processing contents of the linear function calculator 144 will be described below.
In general, the relationship between the d-axis current command value Id * and the q-axis current command value Iq * when the maximum torque is generated for the same current is obtained from the above formulas 1 and 2 and the following formulas 3 and 4. It is shown by the derived number 5. However, the q-axis inductance Lq and the d-axis inductance Ld of the synchronous motor 100 are assumed to be Lq> Ld.
Iq * = − Ia * × sin β (Equation 3)
Id * = Ia * × cosβ (Expression 4)
Id * = Фa / {2 (Lq−Ld)}
-[Фa 2 / {4 (Lq−Ld) 2 } + Iq * 2 ] 1/2 (Expression 5)
Strictly speaking, the q-axis current command value Iq * for generating the maximum torque at the same current is calculated from the above equation 5. However, since the calculation process takes time, the equation 5 is approximated as follows. Then, the q-axis current command value Iq * is calculated. That is,
In Equation 5,
Фa / {2 (Lq−Ld)} = A (Expression 6)
From the above formulas 5 and 6, O point (Iq * , Id * ) = (0, 0), X point (Iq * , Id * ) = (Iq base , Id base ) = (A, (1 −2 1/2 ) A), Y point (Iq * , Id * ) = (Iq max , Id max ) = (Iq max , A− (A 2 + Iq max 2 ) 1/2 ). As shown below, Iq * is divided into two ranges and is calculated from a linear function indicating straight lines between the O point and the X point and between the X point and the Y point.

(1)0≦Iq≦Iqbaseの範囲では、
Id=(1−21/2)×Iq …(数7)
となり、該数7の直線の傾きである(1−21/2)を比例定数a1として設定する。該数7が第1の1次関数演算器143で演算される。
(2)Iqbase≦Iq≦Iqmaxの範囲では、
Id=〔{Idmax-(1-21/2)×A}/(Iqmax-A) 〕×Iq
+A×{ (1-21/2)×Iqmax-Idmax}}/( Iqmax-A) …(数8)
となり、数8の直線の傾きである〔{Idmax-(1-21/2)×A}/(Iqmax-A) 〕を比例定数a2として設定し、数8の直線のIq軸との交点であるA×{ (1-21/2)×Iqmax-Idmax)}/( Iqmax-A)を定数bとして設定する。該数8が第2の1次関数演算器144で演算される。
上記Iqmaxは、上記数2より制御装置における最大電流値制限値Imaxの時のβmaxから、上記数3を用いて求める。
(1) In the range of 0 ≦ Iq * ≦ Iq base ,
Id * = (1-2 1/2 ) × Iq * (Expression 7)
Thus, (1-2 1/2 ), which is the slope of the straight line of Equation 7, is set as the proportionality constant a 1 . Equation 7 is calculated by the first linear function calculator 143.
(2) In the range of Iq base ≦ Iq * ≦ Iq max ,
Id * = [{Id max − (1-2 1/2 ) × A} / (Iq max −A)] × Iq *
+ A × {(1-2 1/2) × Iq max -Id max}} / (Iq max -A) ... ( 8)
Then, [{Id max − (1-2 1/2 ) × A} / (Iq max −A)], which is the slope of the straight line of Formula 8, is set as the proportionality constant a 2 , and Iq * of the straight line of Formula 8 is set . A × {(1-2 1/2 ) × Iq max −Id max )} / (Iq max −A), which is the intersection with the axis, is set as a constant b. The number 8 is calculated by the second linear function calculator 144.
The above Iq max is obtained by using the above equation 3 from β max when the maximum current value limit value I max in the control device is obtained from the above equation 2.

図2において、トルク指令値T*が電流指令値変換器14に入力され、q軸電流指令値演算器141において比例定数KTIを乗じることでq軸電流指令値Iqが算出される。該算出されたq軸電流指令値Iqのレベルは、比較部142において、基準のq軸電流値Iqbaseのレベルと比較される。該比較の結果、q軸電流指令値Iqのレベルが基準のq軸電流値Iqbaseのレベル以下の場合には、第1の1次関数演算器143において、該q軸電流指令値Iqに比例定数a1が乗算され、これがd軸電流指令値Idとされて出力される。また、上記比較の結果、q軸電流指令値Iqのレベルが基準のq軸電流値Iqbaseを超える場合には、第2の1次関数演算器144において、該q軸電流指令値Iqに比例定数aが乗算され、該積にさらに定数bが加算されたものがd軸電流指令値Idとして出力される。上記比例定数KTIは上記q軸電流指令値演算器141が、上記比例定数aは上記比例定数設定部143aが、上記比例定数aは上記比例定数設定部144aが、定数bは上記定数設定部144bが、それぞれ、同期電動機制御装置100の初期化処理のときに設定する。
以下の説明中で用いる図2の構成要素には、図2の場合と同じ符号を付して用いる。
2, the torque command value T * is inputted to the current command value converter 14, the q-axis current command value Iq * is calculated by multiplying the proportional constant K TI in the q-axis current command value calculator 141. The calculated level of the q-axis current command value Iq * is compared with the level of the reference q-axis current value Iq base in the comparison unit 142. As a result of the comparison, when the level of the q-axis current command value Iq * is equal to or lower than the level of the reference q-axis current value Iq base , the first linear function calculator 143 uses the q-axis current command value Iq *. Is multiplied by a proportional constant a 1 , which is output as a d-axis current command value Id * . As a result of the comparison, when the q-axis current command value Iq * level exceeds the reference q-axis current value Iq base, in the second linear function calculator 144, the q-axis current command value Iq * proportionality constant a 2 is multiplied by and further the constant b in the laminate has been added is output as a d-axis current command value Id * to. The proportionality constant K TI is the q-axis current command value calculator 141, the proportional constant a 1 is the proportionality constant setting unit 143a is the proportionality constant a 2 is the proportional constant setting section 144a is constant b the constant The setting unit 144b sets the time when the synchronous motor control device 100 is initialized.
2 used in the following description are given the same reference numerals as in FIG.

図3は、図1の同期電動機制御装置100におけるq軸電流指令値Iqとd軸電流指令値Idとの関係の説明図である。
図3中、実線は上記数7、数8による特性線、点線は理論値となる数5の場合から算出される特性線、一点鎖線は、q軸電流指令値Iqの範囲に関係なく、一定の比例定数を該q軸電流指令値Iqに乗じてd軸電流指令値Idを算出したときの特性線である。
FIG. 3 is an explanatory diagram of the relationship between the q-axis current command value Iq * and the d-axis current command value Id * in the synchronous motor control device 100 of FIG.
In FIG. 3, the solid line is the characteristic line according to Equations 7 and 8, the dotted line is the characteristic line calculated from the case of Equation 5, which is the theoretical value, and the alternate long and short dash line is irrespective of the range of the q-axis current command value Iq * . This is a characteristic line when the d-axis current command value Id * is calculated by multiplying the q-axis current command value Iq * by a constant proportional constant.

図1の同期電動機制御装置100においては、図3において実線の特性線となり、q軸電流指令値Iqが、0≦Iq≦Iqbaseの範囲では、Id=a1×Iq=(1−21/2)×Iqとなり、Iqbase≦Iq≦Iqmaxの範囲では、Id=a×Iq+b=〔{Idmax-(1-21/2)×A}/(Iqmax-A) 〕×Iq+A×{ (1-21/2)×Iqmax-Idmax)}/( Iqmax-A)となる。 In the synchronous motor control device 100 of FIG. 1, the characteristic line is a solid line in FIG. 3, and when the q-axis current command value Iq * is in the range of 0 ≦ Iq * ≦ Iq base , Id * = a 1 × Iq * = ( 1-2 1/2 ) × Iq * , and in the range of Iq base ≦ Iq * ≦ Iq max , Id * = a 2 × Iq * + b = [{Id max − (1-2 1/2 ) × A} / (Iq max −A)] × Iq * + A × {(1-2 1/2 ) × Iq max −Id max )} / (Iq max −A).

上記第1の実施例としての同期電動機制御装置1によれば、簡単に、理論値に近いd軸電流指令値Idを算出することができ、制御の信頼性と制御精度とを改善できる。この結果、同期電動機100を安定した高トルクで運転することができる。 According to the synchronous motor control device 1 as the first embodiment, the d-axis current command value Id * close to the theoretical value can be easily calculated, and the control reliability and control accuracy can be improved. As a result, the synchronous motor 100 can be operated with a stable high torque.

図4〜図5は、本発明の同期電動機制御装置の第2の実施例の説明図である。図4は、第2の実施例としての同期電動機制御装置の構成例図、図5は、図4の同期電動機制御装置を構成する電流指令値変換器の構成例図である。
本第2の実施例の同期電動機制御装置は、同期電動機が高速運転されたときまたは重負荷運転されたときに発生する電圧飽和の抑制を可能にする構成を有する。
4-5 is explanatory drawing of the 2nd Example of the synchronous motor control apparatus of this invention. FIG. 4 is a configuration example diagram of a synchronous motor control device as a second embodiment, and FIG. 5 is a configuration example diagram of a current command value converter constituting the synchronous motor control device of FIG.
The synchronous motor control device of the second embodiment has a configuration that enables suppression of voltage saturation that occurs when the synchronous motor is operated at high speed or when it is operated under heavy load.

図4において、1'は同期電動機制御装置、100は同期電動機、11は、同期電動機100のロータ(図示なし)の回転位置を検出する位置検出器としてのエンコーダ、12'は、エンコーダ11の出力から同期電動機100のロータの回転速度及び電気角速度ωを演算する速度演算器、13は、上記速度演算器12'から得られる上記ロータの現在の回転速度Nと、予め設定してある速度指令値Nとから、トルク指令値Tを演算し出力する速度制御演算器、14'は、上記トルク指令値Tに基づき、q軸電流指令値Iqとd軸電流指令値Idとを演算し出力する電流指令値変換器、15は、同期電動機100の入力電流値を検出する電流検出器、16は、該検出した入力電流値をq軸電流値Iq及びd軸電流値Idに変換する軸電流値変換器としてのU−V−W/d−q変換器、17aは、上記電流指令値変換器14'から出力されたq軸電流指令値Iqと、上記U−V−W/d−q変換器16から出力されたq軸電流値Iqと加え合わせ、該q軸電流指令値Iqから該q軸電流値Iqを減算する加算部、17bは、電流指令値変換器14'から出力されたd軸電流指令値Idと、U−V−W/d−q変換器16から出力されたd軸電流値Idと加え合わせ、該d軸電流指令値Idから該d軸電流値Idを減算する加算部、17cは、予め設定されている同期電動機100のロータの速度指令値Nと、上記速度演算器12から出力された同ロータの現在速度Nとを加え合わせ、該速度指令値Nから該現在速度Nを減算し速度偏差情報を出力する加算部、18'は、上記q軸電流指令値Iqから上記q軸電流値Iqを減算した電流値偏差情報と、上記d軸電流指令値Idから上記d軸電流値Idを減算した電流値偏差情報とに基づき、上記同期電動機100の入力電流値を演算し出力する電流制御演算器である。該電流制御演算器18'から出力される制御された該入力電流は、同期電動機100の固定子コイル(図示なし)に入力され、該同期電動機100が駆動される。上記速度制御演算器13は、上記加算部17cから出力される速度偏差情報に基づいて上記トルク指令値Tを形成して出力する。上記電流検出器15及び上記U−V−W/d−q変換器16は、上記電流制御演算器18'に対してフィードバック制御系を形成する。 In FIG. 4, 1 ′ is a synchronous motor control device, 100 is a synchronous motor, 11 is an encoder as a position detector that detects the rotational position of a rotor (not shown) of the synchronous motor 100, and 12 ′ is an output of the encoder 11. A speed calculator 13 for calculating the rotational speed and electrical angular speed ω of the rotor of the synchronous motor 100 from the current rotational speed N of the rotor obtained from the speed calculator 12 ′ and a preset speed command value. The speed control computing unit 14 ′, which calculates and outputs the torque command value T * from N * , calculates the q-axis current command value Iq * and the d-axis current command value Id * based on the torque command value T *. A current command value converter for calculating and outputting, 15 a current detector for detecting the input current value of the synchronous motor 100, and 16 for converting the detected input current value into a q-axis current value Iq and a d-axis current value Id. Axis Current values converter as the U-V-W / d- q converter, 17a is the current command value converter 14 q-axis current command value output from 'and Iq *, the U-V-W / d -An addition unit 17b that adds the q-axis current value Iq output from the q-axis converter 16 and subtracts the q-axis current value Iq from the q-axis current command value Iq *. The output d-axis current command value Id * is added to the d-axis current value Id output from the UVV / dq converter 16, and the d-axis current command value Id * is added to the d-axis current. An adding unit 17c for subtracting the value Id adds the preset rotor speed command value N * of the synchronous motor 100 and the current speed N of the rotor output from the speed calculator 12, adder for outputting a speed deviation information by subtracting the current speed N from the speed command value N *, 1 'Includes a current value deviation information from the q-axis current command value Iq * obtained by subtracting the q-axis current value Iq, and a current value deviation information from the d-axis current command value Id * obtained by subtracting the d-axis current value Id This is a current control calculator that calculates and outputs the input current value of the synchronous motor 100 based on the above. The controlled input current output from the current control calculator 18 ′ is input to a stator coil (not shown) of the synchronous motor 100, and the synchronous motor 100 is driven. The speed control calculator 13 forms and outputs the torque command value T * based on the speed deviation information output from the adder 17c. The current detector 15 and the UVW / dq converter 16 form a feedback control system for the current control calculator 18 ′.

電流指令値変換器14'には、トルク指令値Tの他、電流制御演算器18'からの直流電圧Vdcと、速度演算器12'からの電気角速度ωが入力される。該電流制御演算器18'は、上記のように、q軸電流指令値Iqとd軸電流指令値Idとを形成して出力するが、d軸電流指令値Idについては、次のようにして形成し出力する。すなわち、q軸電流指令値Iqが基準のq軸電流値Iqbase以下の場合は、第1のd軸電流指令値Id を該q軸電流指令値Iqの第1の1次関数として演算し、一方、q軸電流指令値Iqが基準のq軸電流値Iqbaseを超える場合は、第2のd軸電流指令値Id を該q軸電流指令値Iqの第2の1次関数として演算し、さらに、上記直流電圧Vdcから基準の電気角速度ωbaseを演算し、上記速度演算器12'から入力される電気角速度ωと比較し、該電気角速度ωが該基準の電気角速度ωbase以上の場合に、該電気角速度ωに比例するd軸電流指令増加値Idωを演算し、上記第1のd軸電流指令値Id または上記第2のd軸電流指令値Id に上記d軸電流指令増加値Idωを加算して、該加算結果をd軸電流指令値Idとして出力する。 In addition to the torque command value T * , the DC voltage Vdc from the current control calculator 18 ′ and the electrical angular velocity ω from the speed calculator 12 ′ are input to the current command value converter 14 ′. As described above, the current control arithmetic unit 18 ′ forms and outputs the q-axis current command value Iq * and the d-axis current command value Id * . The d-axis current command value Id * In this way, it forms and outputs. That is, when the q-axis current command value Iq * is less than or equal to the reference q-axis current value Iq base, first linear function of the first d-axis current command value Id * 1 of the q-axis current command value Iq * calculated as while, q-axis current command value Iq * If exceeds the reference q-axis current value Iq base, the second d-axis current command value Id * 2 the q-axis current command value Iq * of the second Further, a reference electrical angular velocity ω base is calculated from the DC voltage Vdc, and compared with the electrical angular velocity ω input from the speed calculator 12 ′. When the electrical angular velocity is equal to or higher than ω base , the d-axis current command increase value Id ω proportional to the electrical angular velocity ω is calculated, and the first d-axis current command value Id * 1 or the second d-axis current command value is calculated. The d-axis current command increase value Id ω is added to Id * 2 , and the addition result is added to the d-axis current index. Output as normal value Id * .

図5は、図4の同期電動機制御装置1'を構成する電流指令値変換器14'の構成例図である。
図5において、141は、トルク指令値Tに比例定数KTIを乗算し上記q軸電流指令値Iqを演算して出力するq軸電流指令値演算器、142は、該出力されたq軸電流指令値Iqのレベルを基準のq軸電流値Iqbaseのレベルと比較する比較部、143'は、該比較部142での比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbase以下の場合、第1のd軸電流指令値Id を該q軸電流指令値Iqの第1の1次関数として演算する第1の1次関数演算器、143a'は、上記q軸電流指令値Iqに乗ずる比例定数a'を設定する比例定数設定部、144'は、上記比較部142での比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbaseを超える場合は、第2のd軸電流指令値Id を該q軸電流指令値Iqの第2の1次関数として演算する第2の1次関数演算器、144a'は、上記q軸電流指令値Iqに乗ずる比例定数a'を設定する比例定数設定部、144b'は、該q軸電流指令値Iqと比例定数a'との積に加算される定数b'を設定する定数設定部、145は、交流変換される直流電圧Vdcから基準の電気角速度ωbaseを演算し、当該同期電動機のロータの電気角速度ωと比較し、該比較の結果、該電気角速度ωが該基準の電気角速度ωbaseを超える場合は、該電気角速度ωと該基準の電気角速度ωbaseとの差に比例するd軸電流指令増加値Idωを演算して出力し、上記比較の結果、上記電気角速度ωが該基準の電気角速度ωbase以下の場合には、d軸電流指令増加値Idωをゼロレベルにする第3の1次関数演算器、145aは、電気角速度ωを絶対値|ω|に変換する絶対値変換部、145bは、直流電圧Vdcから基準の電気角速度ωbaseを演算する基準角速度演算部、145cは、上記電気角速度ωの絶対値|ω|と上記基準の電気角速度ωbaseを比較し、該絶対値|ω|が該基準の電気角速度ωbaseを超える場合は、その差の信号(便宜上、角速度差という)を出力し、一方、該絶対値|ω|が該基準の電気角速度ωbase以下の場合には、出力信号をゼロにする角速度比較部、145dは、上記角速度差(|ω|−ωbase)に乗ずる比例定数Kωを設定する比例定数設定部である。比例定数設定部145dからは、Kω×(|ω|−ωbase)のd軸電流指令増加値Idωが出力される。146は、上記第1のd軸電流指令値Id または上記第2のd軸電流指令値Id に上記d軸電流指令増加値Idωを加算し、該加算結果をd軸電流指令値Idとして出力する加算部である。
FIG. 5 is a configuration example diagram of a current command value converter 14 ′ constituting the synchronous motor control device 1 ′ of FIG.
5, 141, multiplies the proportional constant K TI to the torque command value T * q-axis current command value calculator that calculates and outputs the q-axis current command value Iq *, 142 was the output q The comparison unit 143 ′ that compares the level of the shaft current command value Iq * with the level of the reference q-axis current value Iq base indicates that the q-axis current command value Iq * is the reference value as a result of the comparison by the comparison unit 142. if: q-axis current value Iq base of the first linear function calculator that calculates a first d-axis current command value Id * 1 as a first linear function of the q-axis current command value Iq *, 143a 'is a proportionality constant a 1 is multiplied to the q-axis current command value Iq *' proportionality constant setting unit for setting a, 144 'a result of the comparison in the comparison unit 142, the q-axis current command value Iq * is If the q-axis current value Iq base is exceeded, the second d-axis current command value Id A second linear function computing unit 144a ′ that calculates * 2 as a second linear function of the q-axis current command value Iq * has a proportional constant a 2 ′ multiplied by the q-axis current command value Iq *. The proportional constant setting unit 144b ′ for setting sets the constant b ′ to be added to the product of the q-axis current command value Iq * and the proportional constant a 2 ′. A reference electrical angular speed ω base is calculated from the voltage Vdc and compared with the electrical angular speed ω of the rotor of the synchronous motor. If the electrical angular speed ω exceeds the reference electrical angular speed ω base as a result of the comparison, The d-axis current command increase value Id ω proportional to the difference between the angular velocity ω and the reference electrical angular velocity ω base is calculated and output. As a result of the comparison, the electrical angular velocity ω is less than or equal to the reference electrical angular velocity ω base . If the third of the d-axis current command increment value Id omega zero level Linear function calculator, 145a is the electrical angular velocity omega absolute value | omega | absolute value converter for converting the, 145b, the reference angular velocity calculating section for calculating the electrical angular velocity omega base reference from the DC voltage Vdc, 145c is The absolute value | ω | of the electrical angular velocity ω is compared with the reference electrical angular velocity ω base, and when the absolute value | ω | exceeds the reference electrical angular velocity ω base , the difference signal (for convenience, the angular velocity difference On the other hand, when the absolute value | ω | is equal to or less than the reference electrical angular velocity ω base , the angular velocity comparison unit 145d that makes the output signal zero outputs the angular velocity difference (| ω | -ω it is a proportionality constant setting unit for setting a proportional constant K omega to multiply the base). The proportional constant setting unit 145d outputs a d-axis current command increase value Id ω of K ω × (| ω | −ω base ). 146 adds the d-axis current command increase value Id ω to the first d-axis current command value Id * 1 or the second d-axis current command value Id * 2 , and adds the addition result to the d-axis current command. An adder that outputs the value Id * .

上記基準角速度演算部145bでは、変換係数Kを用いて直流電圧Vdcから基準の電気角速度ωbaseを下記数9により演算する。
ωbase=K×Vdc …(数9)
ここで、同期電動機100のq軸インダクタンスをLq、d軸インダクタンスをLd、永久磁石による鎖交磁束をΦ、q軸電流値の最大値をIqmax、d軸電流値の最大値をIdmax、電圧利用率をλとするとき、上記変換係数Kは、
=λ/{(Φa+Ld×Idmax+(Lq×Iqmax}1/2…(数10)
で表される値となる。該変換係数Kは、基準角速度演算部145b内に設定される。上記数9から明らかなように、基準の電気角速度ωbaseは、上記直流電圧Vdcにより変化する。
In the reference velocity calculation unit 145b, the electrical angular velocity omega base reference for calculating the following Equation 9 from the DC voltage Vdc using a conversion coefficient K f.
ω base = K f × Vdc (Equation 9)
Here, the q-axis inductance of the synchronous motor 100 is Lq, the d-axis inductance is Ld, the linkage flux by the permanent magnet is Φ, the maximum value of the q-axis current value is Iq max , and the maximum value of the d-axis current value is Id max , When the voltage utilization factor is λ, the conversion coefficient K f is
K f = λ / {(Φa + Ld × Id max ) 2 + (Lq × Iq max ) 2 } 1/2 (Equation 10)
The value represented by The conversion coefficient K f is set to the reference angular velocity calculating unit 145b. As is clear from the above formula 9, the reference electrical angular velocity ω base varies with the DC voltage Vdc.

上記角速度比較部145cでは、同期電動機100のロータが、その電気角速度ωが、基準の電気角速度ωbase以下かまたは該基準の電気角速度ωbaseを超えるかを判断する。電気角速度ωが、基準の電気角速度ωbase以下の場合は、電圧飽和状態とならないため、第3の1次関数演算器145から出力されるd軸電流指令増加値Idωがゼロレベルとなるようにする。また、電気角速度ωが、基準の電気角速度ωbaseを超える場合は、電圧飽和状態を回避すべく、該第3の1次関数演算器145からは、角速度差(|ω|−ωbase)に比例したd軸電流指令増加値Idωが出力されるようにする。すなわち、
(1)0≦|ω|≦ωbaseの場合は、電圧飽和状態となっていないため、第3の1次関数演算器145から出力されるd軸電流指令増加値Idωは、
Idω=0 …(数11)
とされる。
(2)ωbase<|ω|≦ωmaxの場合は、電圧飽和状態を回避するために、第3の1次関数演算器145から出力されるd軸電流指令増加値Idωは、
Idω=Kω×(|ω|−ωbase) …(数12)
とされる。
In the angular velocity comparator unit 145c, the rotor of the synchronous motor 100, the electrical angular velocity omega determines whether or the following criteria electrical angular velocity omega base or exceeds the electrical angular velocity omega base of the reference. When the electrical angular velocity ω is equal to or less than the reference electrical angular velocity ω base , the voltage saturation state is not reached, so that the d-axis current command increase value Id ω output from the third linear function calculator 145 becomes zero level. To. When the electrical angular velocity ω exceeds the reference electrical angular velocity ω base , the third linear function calculator 145 generates an angular velocity difference (| ω | −ω base ) to avoid a voltage saturation state. A proportional d-axis current command increase value Id ω is output. That is,
(1) In the case of 0 ≦ | ω | ≦ ω base , since the voltage is not saturated, the d-axis current command increase value Id ω output from the third linear function calculator 145 is
Id ω = 0 (Expression 11)
It is said.
(2) In the case of ω base <| ω | ≦ ω max , in order to avoid a voltage saturation state, the d-axis current command increase value Id ω output from the third linear function calculator 145 is
Id ω = K ω × (| ω | −ω base ) (Equation 12)
It is said.

上記数11、数12により演算されたd軸電流指令増加値Idωは、選択的に第3の1次関数演算器145から出力され、加算部146において、第1の1次関数演算器143'から出力された第1のd軸電流指令値Id または第2の1次関数演算器144'から出力された第2のd軸電流指令値Id に加算される。 The d-axis current command increase value Id ω calculated by the above equations 11 and 12 is selectively output from the third linear function calculator 145, and the adder 146 outputs the first linear function calculator 143. Is added to the first d-axis current command value Id * 1 output from 'or the second d-axis current command value Id * 2 output from the second linear function calculator 144'.

上記第2の実施例としての同期電動機制御装置1'によっても、上記第1の実施例の同期電動機制御装置1と同様、簡単に、理論値に近いd軸電流指令値Idを算出することができ、制御の信頼性と制御精度とを改善できる。この結果、同期電動機100を安定した高トルクで運転することができる。特に、本第2の実施例としての同期電動機制御装置1'の場合は、高速回転時にも電圧飽和状態を回避して安定した高トルク運転が可能となる。 Similarly to the synchronous motor control device 1 of the first embodiment, the synchronous motor control device 1 'as the second embodiment can easily calculate the d-axis current command value Id * close to the theoretical value. Thus, the reliability and accuracy of control can be improved. As a result, the synchronous motor 100 can be operated with a stable high torque. In particular, in the case of the synchronous motor control apparatus 1 ′ according to the second embodiment, a stable high torque operation is possible by avoiding a voltage saturation state even at high speed rotation.

第1の実施例としての同期電動機制御装置の構成例図である。It is a structural example figure of the synchronous motor control apparatus as a 1st Example. 図1の同期電動機制御装置を構成する電流指令値変換器の構成例図である。It is an example of a structure of the electric current command value converter which comprises the synchronous motor control apparatus of FIG. 図1の同期電動機制御装置におけるq軸電流指令値Iqとd軸電流指令値Idとの関係の説明図である。It is explanatory drawing of the relationship between q-axis current command value Iq * and d-axis current command value Id * in the synchronous motor control apparatus of FIG. 第2の実施例としての同期電動機制御装置の構成例図である。It is a structural example figure of the synchronous motor control apparatus as a 2nd Example. 図4の同期電動機制御装置を構成する電流指令値変換器の構成例図である。It is an example of a structure of the current command value converter which comprises the synchronous motor control apparatus of FIG.

符号の説明Explanation of symbols

1、1'…同期電動機制御装置、
100…同期電動機、
11…エンコーダ、
12、12'…速度演算器、
13…速度制御演算器、
14、14'…電流指令値変換器、
15…電流検出器、
16…U−V−W/d−q変換器、
17a、17b、17c、146…加算部、
18、18'…電流制御演算器、
141…q軸電流指令値演算器、
142…比較部、
143、143'…第1の1次関数演算器、
143a、144a、143a'、144a'、145d…比例定数設定部、
144、144'…第2の1次関数演算器、
144b、144b'…定数設定部、
145…第3の1次関数演算器、
145a…絶対値変換部、
145b…基準角速度演算部、
145c…角速度比較部。
1, 1 '... synchronous motor control device,
100 ... Synchronous motor,
11: Encoder,
12, 12 '... speed calculator,
13 ... Speed control calculator,
14, 14 '... current command value converter,
15 ... current detector,
16 ... U-V-W / dq converter,
17a, 17b, 17c, 146 ... adder,
18, 18 '... current control arithmetic unit,
141... Q-axis current command value calculator,
142 ... comparison part,
143, 143 '... the first linear function calculator,
143a, 144a, 143a ', 144a', 145d ... proportional constant setting unit,
144, 144 '... a second linear function calculator,
144b, 144b '... constant setting unit,
145 ... third linear function calculator,
145a ... absolute value converter,
145b: Reference angular velocity calculation unit,
145c ... Angular velocity comparison unit.

Claims (7)

同期電動機の入力電流を制御する同期電動機制御装置であって、
上記同期電動機のロータの回転位置を検出する位置検出器と、
上記位置検出器の出力から上記ロータの回転速度を演算する速度演算器と、
上記速度演算器から得られる上記ロータの現在の回転速度と、予め設定してある速度指令値とから、トルク指令値Tを演算する速度制御演算器と、
上記トルク指令値Tに基づき、q軸電流指令値Iqとd軸電流指令値Idとを演算する電流指令値変換器と、
上記同期電動機の入力電流値を検出する電流検出器と、
上記検出した入力電流値をq軸電流値Iq及びd軸電流値Idに変換する軸電流値変換器と、
上記q軸電流指令値Iqから上記q軸電流値Iqを減算した電流値情報と、上記d軸電流指令値Idから上記d軸電流値Idを減算した電流値情報とに基づき、上記同期電動機の入力電流値を演算し該入力電流を出力する電流制御演算器と、
を備えて成り、
上記電流指令値変換器が、
上記トルク指令値Tに比例定数KTIを乗算して上記q軸電流指令値Iqを演算するq軸電流指令値演算器と、
上記演算されたq軸電流指令値Iqを基準のq軸電流値Iqbaseと比較する比較部と、
上記比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbase以下の場合は、上記d軸電流指令値Idを該q軸電流指令値Iqの第1の1次関数として演算する第1の1次関数演算器と、
上記比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbaseを超える場合は、上記d軸電流指令値Idを該q軸電流指令値Iqの第2の1次関数として演算する第2の1次関数演算器と、
を備えた構成であることを特徴とする同期電動機制御装置。
A synchronous motor control device for controlling an input current of a synchronous motor,
A position detector for detecting the rotational position of the rotor of the synchronous motor;
A speed calculator that calculates the rotational speed of the rotor from the output of the position detector;
A speed control calculator for calculating a torque command value T * from a current rotational speed of the rotor obtained from the speed calculator and a preset speed command value;
A current command value converter for calculating a q-axis current command value Iq * and a d-axis current command value Id * based on the torque command value T * ;
A current detector for detecting an input current value of the synchronous motor;
An axial current value converter for converting the detected input current value into a q-axis current value Iq and a d-axis current value Id;
Based on the current value information obtained by subtracting the q-axis current value Iq from the q-axis current command value Iq * and the current value information obtained by subtracting the d-axis current value Id from the d-axis current command value Id *. A current control calculator that calculates the input current value of the motor and outputs the input current;
Comprising
The current command value converter is
And the q-axis current command value calculator for calculating the q-axis current command value Iq * is multiplied by a proportionality constant K TI to the torque command value T *,
A comparison unit that compares the calculated q-axis current command value Iq * with a reference q-axis current value Iq base ;
Result of the comparison, if the q-axis current command value Iq * is less than or equal to the q-axis current value Iq base of the reference, the d-axis current command value Id * and the q-axis current command value Iq * first of 1 A first linear function calculator that operates as a secondary function;
Result of the comparison, if the q-axis current command value Iq * exceeds the q-axis current value Iq base of the reference, the d-axis current command value Id * and the q-axis current command value Iq * second of 1 A second linear function calculator that operates as a secondary function;
A synchronous motor control device characterized by comprising:
上記第1の1次関数演算器は、上記q軸電流指令値Iqに比例定数aを乗じてd軸電流指令値Idを求め、上記第2の1次関数演算器は、上記q軸電流指令値Iqに比例定数aを乗じさらに定数bを加算して上記d軸電流指令値Idを求める構成である請求項1に記載の同期電動機制御装置。 Said first linear function calculator calculates a d-axis current command value Id * is multiplied by the proportional constant a 1 on the q-axis current command value Iq *, the second linear function calculator, said q 2. The synchronous motor control device according to claim 1, wherein the d-axis current command value Id * is obtained by multiplying a shaft current command value Iq * by a proportional constant a2 and further adding a constant b. 上記基準のq軸電流値Iqbase,上記比例定数aまたは上記定数bは、永久磁石による鎖交磁束Ф、q軸インダクタンスLq、d軸インダクタンスLd、q軸電流の最大値Iqmaxの関数として求められる請求項2に記載の同期電動機制御装置。 The reference q-axis current value Iq base , the proportionality constant a 2 or the constant b is a function of the linkage flux 永久 by the permanent magnet, the q-axis inductance Lq, the d-axis inductance Ld, and the maximum q-axis current value Iq max . The synchronous motor control device according to claim 2, which is required. 上記電流検出器及び上記軸電流値変換器は、上記電流制御演算器に対するフィードバック制御系を形成する請求項1、2または3に記載の同期電動機制御装置。   4. The synchronous motor control device according to claim 1, wherein the current detector and the shaft current value converter form a feedback control system for the current control arithmetic unit. 同期電動機の入力電流を制御する同期電動機制御装置であって、
上記同期電動機のロータの回転位置を検出する位置検出器と、
上記位置検出器の出力から上記ロータの回転速度を演算する速度演算器と、
上記速度演算器から得られる上記ロータの現在の回転速度と、予め設定してある速度指令値とから、トルク指令値Tを演算する速度制御演算器と、
上記トルク指令値Tに基づき、q軸電流指令値Iqとd軸電流指令値Idとを演算する電流指令値変換器と、
上記同期電動機の入力電流値を検出する電流検出器と、
上記検出した入力電流値をq軸電流値Iq及びd軸電流値Idに変換する軸電流値変換器と、
上記q軸電流指令値Iqから上記q軸電流値Iqを減算した電流値情報と、上記d軸電流指令値Idから上記d軸電流値Idを減算した電流値情報とに基づき、上記同期電動機の入力電流値を演算し該入力電流を出力する電流制御演算器と、
を備えて成り、
上記電流指令値変換器が、
上記トルク指令値Tに比例定数KTIを乗算して上記q軸電流指令値Iqを演算するq軸電流指令値演算器と、
上記演算されたq軸電流指令値Iqを基準のq軸電流値Iqbaseと比較する比較部と、
上記比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbase以下の場合は、第1のd軸電流指令値Id を該q軸電流指令値Iqの第1の1次関数として演算する第1の1次関数演算器と、
上記比較の結果、上記q軸電流指令値Iqが上記基準のq軸電流値Iqbaseを超える場合は、第2のd軸電流指令値Id を該q軸電流指令値Iqの第2の1次関数として演算する第2の1次関数演算器と、
交流変換される直流電圧Vdcから基準の電気角速度ωbaseを演算し、当該同期電動機の電気角速度ωと比較し、該比較の結果、該電気角速度ωが該基準の電気角速度ωbaseを超える場合は、該電気角速度ωと該基準の電気角速度ωbaseとの差に比例するd軸電流指令増加値Idωを演算して出力し、上記比較の結果、上記電気角速度ωが該基準の電気角速度ωbase以下の場合には、d軸電流指令増加値Idωをゼロレベルにする第3の1次関数演算器と、
上記第1のd軸電流指令値Id または上記第2のd軸電流指令値Id に上記d軸電流指令増加値Idωを加算し、該加算結果をd軸電流指令値Idとして出力する加算部と、
を備えた構成であることを特徴とする同期電動機制御装置。
A synchronous motor control device for controlling an input current of a synchronous motor,
A position detector for detecting the rotational position of the rotor of the synchronous motor;
A speed calculator that calculates the rotational speed of the rotor from the output of the position detector;
A speed control calculator for calculating a torque command value T * from a current rotational speed of the rotor obtained from the speed calculator and a preset speed command value;
A current command value converter for calculating a q-axis current command value Iq * and a d-axis current command value Id * based on the torque command value T * ;
A current detector for detecting an input current value of the synchronous motor;
An axial current value converter for converting the detected input current value into a q-axis current value Iq and a d-axis current value Id;
Based on the current value information obtained by subtracting the q-axis current value Iq from the q-axis current command value Iq * and the current value information obtained by subtracting the d-axis current value Id from the d-axis current command value Id *. A current control calculator that calculates the input current value of the motor and outputs the input current;
Comprising
The current command value converter is
And the q-axis current command value calculator for calculating the q-axis current command value Iq * is multiplied by a proportionality constant K TI to the torque command value T *,
A comparison unit that compares the calculated q-axis current command value Iq * with a reference q-axis current value Iq base ;
As a result of the comparison, if the q-axis current command value Iq * is less than or equal to the reference q-axis current value Iq base , the first d-axis current command value Id * 1 is changed to the q-axis current command value Iq * . A first linear function calculator that operates as a linear function of 1;
Result of the comparison, if the q-axis current command value Iq * exceeds the q-axis current value Iq base of the reference, the second d-axis current command value Id * 2 of the q-axis current command value Iq * A second linear function calculator that operates as a linear function of 2;
When the reference electrical angular velocity ω base is calculated from the DC voltage Vdc to be converted into AC, and compared with the electrical angular velocity ω of the synchronous motor. As a result of the comparison, when the electrical angular velocity ω exceeds the reference electrical angular velocity ω base The d-axis current command increase value Id ω proportional to the difference between the electrical angular velocity ω and the reference electrical angular velocity ω base is calculated and output. As a result of the comparison, the electrical angular velocity ω is converted into the reference electrical angular velocity ω. In the case of base or less, a third linear function calculator that sets the d-axis current command increase value Id ω to zero level;
The d-axis current command increase value Id ω is added to the first d-axis current command value Id * 1 or the second d-axis current command value Id * 2 , and the addition result is used as the d-axis current command value Id *. An adder that outputs as
A synchronous motor control device characterized by comprising:
上記第1の1次関数演算器は、上記q軸電流指令値Iqに比例定数αを乗じて上記第1のd軸電流指令値Id を求め、上記第2の1次関数演算器は、上記q軸電流指令値Iqに比例定数αを乗じ該積に定数βを加算して上記第2のd軸電流指令値Id を求める構成である請求項5に記載の同期電動機制御装置。 It said first linear function calculator calculates the first d-axis current command value Id * 1 is multiplied by a proportionality constant alpha 1 to the q-axis current command value Iq *, the second primary function operation 6. The apparatus according to claim 5, wherein the second multi-axis current command value Id * 2 is obtained by multiplying the q-axis current command value Iq * by a proportional constant α2 and adding a constant β to the product. Synchronous motor control device. 上記第3の1次関数演算器は、Kωを比例定数とするとき、上記d軸電流指令増加値Idωを、Kω(ω−ωbase)により求める構成である請求項5に記載の同期電動機制御装置。 6. The configuration according to claim 5, wherein the third linear function calculator is configured to obtain the d-axis current command increase value Id ω from K ω (ω−ω base ) when K ω is a proportional constant. Synchronous motor control device.
JP2006323738A 2006-11-30 2006-11-30 Synchronous motor controller Active JP4810410B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006323738A JP4810410B2 (en) 2006-11-30 2006-11-30 Synchronous motor controller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006323738A JP4810410B2 (en) 2006-11-30 2006-11-30 Synchronous motor controller

Publications (3)

Publication Number Publication Date
JP2008141824A true JP2008141824A (en) 2008-06-19
JP2008141824A5 JP2008141824A5 (en) 2009-06-04
JP4810410B2 JP4810410B2 (en) 2011-11-09

Family

ID=39602732

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006323738A Active JP4810410B2 (en) 2006-11-30 2006-11-30 Synchronous motor controller

Country Status (1)

Country Link
JP (1) JP4810410B2 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009058579A1 (en) 2008-12-17 2010-06-24 Aida Engineering, Ltd., Sagamihara Controller for synchronous motor, has current-preset value-connection unit for combining two d-axis-current-preset values associated to characteristic values and outputting definite d-axis-current-preset value
JP2012065464A (en) * 2010-09-16 2012-03-29 Seiko Epson Corp Motor control device
JP2012139105A (en) * 2012-04-24 2012-07-19 Seiko Epson Corp Motor control device
JP2015035923A (en) * 2013-08-09 2015-02-19 株式会社日立産機システム Synchronous motor controller
KR20150094431A (en) * 2014-02-11 2015-08-19 삼성전자주식회사 Power apparatus, motor driving apparatus therein
CN110957953A (en) * 2018-09-25 2020-04-03 欧姆龙(上海)有限公司 Control device and control method for alternating current motor

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10191700A (en) * 1996-12-20 1998-07-21 Yaskawa Electric Corp Control method for synchronous motor
JPH1141998A (en) * 1997-07-14 1999-02-12 C & S Kokusai Kenkyusho:Kk Vector control method for synchronous motor
JP2000092884A (en) * 1998-09-18 2000-03-31 Mitsubishi Electric Corp Servo motor controller
JP2000358400A (en) * 1999-06-11 2000-12-26 Okuma Corp Control device of synchronous electric motor
JP2002360000A (en) * 2001-06-01 2002-12-13 Hitachi Ltd Drive unit for synchronous motor
JP2003274699A (en) * 2002-03-12 2003-09-26 Toshiba Corp Motor controller

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10191700A (en) * 1996-12-20 1998-07-21 Yaskawa Electric Corp Control method for synchronous motor
JPH1141998A (en) * 1997-07-14 1999-02-12 C & S Kokusai Kenkyusho:Kk Vector control method for synchronous motor
JP2000092884A (en) * 1998-09-18 2000-03-31 Mitsubishi Electric Corp Servo motor controller
JP2000358400A (en) * 1999-06-11 2000-12-26 Okuma Corp Control device of synchronous electric motor
JP2002360000A (en) * 2001-06-01 2002-12-13 Hitachi Ltd Drive unit for synchronous motor
JP2003274699A (en) * 2002-03-12 2003-09-26 Toshiba Corp Motor controller

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009058579A1 (en) 2008-12-17 2010-06-24 Aida Engineering, Ltd., Sagamihara Controller for synchronous motor, has current-preset value-connection unit for combining two d-axis-current-preset values associated to characteristic values and outputting definite d-axis-current-preset value
JP2010148198A (en) * 2008-12-17 2010-07-01 Hitachi Industrial Equipment Systems Co Ltd Control device of synchronous motor
JP2012065464A (en) * 2010-09-16 2012-03-29 Seiko Epson Corp Motor control device
CN102403949A (en) * 2010-09-16 2012-04-04 精工爱普生株式会社 Electric motor control device
US8823300B2 (en) 2010-09-16 2014-09-02 Seiko Epson Corporation Electric motor control device
JP2012139105A (en) * 2012-04-24 2012-07-19 Seiko Epson Corp Motor control device
JP2015035923A (en) * 2013-08-09 2015-02-19 株式会社日立産機システム Synchronous motor controller
KR20150094431A (en) * 2014-02-11 2015-08-19 삼성전자주식회사 Power apparatus, motor driving apparatus therein
WO2015122655A1 (en) * 2014-02-11 2015-08-20 삼성전자주식회사 Power system and motor driving apparatus included in power system
CN105981289A (en) * 2014-02-11 2016-09-28 三星电子株式会社 Power system and motor driving apparatus included in power system
US9966880B2 (en) 2014-02-11 2018-05-08 Samsung Electronics Co., Ltd. Power system and motor driving apparatus included in power system
US10418920B2 (en) 2014-02-11 2019-09-17 Samsung Electronics Co., Ltd. Power system and motor driving apparatus included in power system
CN110957953A (en) * 2018-09-25 2020-04-03 欧姆龙(上海)有限公司 Control device and control method for alternating current motor
CN110957953B (en) * 2018-09-25 2022-12-20 欧姆龙(上海)有限公司 Control device and control method for alternating current motor

Also Published As

Publication number Publication date
JP4810410B2 (en) 2011-11-09

Similar Documents

Publication Publication Date Title
JP5130716B2 (en) Motor control device and electric power steering device
KR100850415B1 (en) Vector controller of induction motor
JP5130031B2 (en) Position sensorless control device for permanent magnet motor
JP4989075B2 (en) Electric motor drive control device and electric motor drive system
JP5104239B2 (en) Control device for permanent magnet type synchronous motor
US6900613B2 (en) Motor control apparatus
JP4235233B2 (en) Electric motor control device
JP4810410B2 (en) Synchronous motor controller
JP5330652B2 (en) Permanent magnet motor control device
JP5351002B2 (en) Motor control device
JP5276688B2 (en) Synchronous machine controller
JP2010200430A (en) Drive controller for motors
JP6541092B2 (en) Control device of permanent magnet synchronous motor
CN110235355B (en) Speed estimation method for induction motor and power conversion device using the same
JP2008148437A (en) Controller for permanent magnet type synchronous motor
JP4579627B2 (en) Rotating machine control device
JP5595436B2 (en) Motor control device
JP2009284598A (en) Controller for alternating-current motors
JP2005065349A (en) Synchronous motor controller
JP4543720B2 (en) Speed sensorless vector controller
JP2020036452A (en) Control device for rotary electric machine
JP5595437B2 (en) Motor control device
JP5040605B2 (en) Control device for permanent magnet type synchronous motor
JP5538657B2 (en) Motor control device
US20230198438A1 (en) Rotary machine control device

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090416

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090416

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110727

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110822

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140826

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4810410

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150