JP7450150B2 - A mathematical model of a synchronous motor that generates reluctance torque, and a simulation, characteristic analysis, and control device based on the same model. - Google Patents

A mathematical model of a synchronous motor that generates reluctance torque, and a simulation, characteristic analysis, and control device based on the same model. Download PDF

Info

Publication number
JP7450150B2
JP7450150B2 JP2020025695A JP2020025695A JP7450150B2 JP 7450150 B2 JP7450150 B2 JP 7450150B2 JP 2020025695 A JP2020025695 A JP 2020025695A JP 2020025695 A JP2020025695 A JP 2020025695A JP 7450150 B2 JP7450150 B2 JP 7450150B2
Authority
JP
Japan
Prior art keywords
reluctance torque
expressed
axis
phase
equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2020025695A
Other languages
Japanese (ja)
Other versions
JP2021121163A (en
Inventor
新二 新中
Original Assignee
有限会社シー・アンド・エス国際研究所
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 有限会社シー・アンド・エス国際研究所 filed Critical 有限会社シー・アンド・エス国際研究所
Priority to JP2020025695A priority Critical patent/JP7450150B2/en
Publication of JP2021121163A publication Critical patent/JP2021121163A/en
Application granted granted Critical
Publication of JP7450150B2 publication Critical patent/JP7450150B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Ac Motors In General (AREA)

Description

本発明は、回転子が固定子側の電力周波数と平均的に同じ電気速度でかつリラクタンストルクを発生しつつ回転する同期電動機の数学モデル、および同モデルに立脚した模擬装置、特性解析装置、制御装置に関する。この種の同期電動機としては、同期リラクタンス電動機、永久磁石同期電動機、巻線形同期電動機、回転子に永久磁石と励磁巻線の両者を備えこれらにより界磁を生成するハイブリッド界磁同期電動機などが知られている。The present invention relates to a mathematical model of a synchronous motor in which the rotor rotates at the same electric speed as the stator power frequency on average and while generating reluctance torque, as well as a simulation device, a characteristic analysis device, and a control device based on the model. Regarding equipment. Known examples of this type of synchronous motor include synchronous reluctance motors, permanent magnet synchronous motors, wound type synchronous motors, and hybrid field synchronous motors that have both permanent magnets and excitation windings in their rotor and generate a field using them. It is being

本発明では、同期電動機において三相巻線が施された部分を「固定子」と呼称する。本発明における「固定子」は、「電機子」と同義である。In the present invention, a portion of a synchronous motor provided with three-phase windings is referred to as a "stator." "Stator" in the present invention has the same meaning as "armature".

固定子電流により発生する磁束の名称として、「電機子反作用磁束」、「固定子反作用磁束」などの用語が、広く同義で使用されている。本発明では、この磁束を、原則として、固定子反作用磁束と呼称する。Terms such as "armature reaction magnetic flux" and "stator reaction magnetic flux" are widely used interchangeably as names for the magnetic flux generated by stator current. In the present invention, this magnetic flux is, in principle, referred to as stator reaction magnetic flux.

本発明では、固定子インダクタンスを簡単にインダクタンスを呼称する。d軸固定子インダクタンス、q軸固定子インダクタンスを、簡単のため、各々、d軸インダクタンス、q軸インダクタンスと呼称する。In the present invention, stator inductance is simply referred to as inductance. For simplicity, the d-axis stator inductance and the q-axis stator inductance will be referred to as d-axis inductance and q-axis inductance, respectively.

本発明では、2次元空間(平面)を極座標的に捉え、角度、空間的位置、空間的位相の3用語を同義で使用する。これらの単位は「ラジアン(rad)」または「度(degree)」である。本発明における角度、空間的位置、空間的位相の正方向は、左周り(反時計周り)、右周り(時計周り)のいずれに定義してもよい。ただし、本明細書では、説明の簡明性を維持すべく、角度、空間的位置、空間的位相の正方向は左周り(反時計周り)と定義し、本発明を説明する。これにより、本発明の一般性を失うことはない。In the present invention, a two-dimensional space (plane) is viewed in terms of polar coordinates, and the three terms angle, spatial position, and spatial phase are used synonymously. These units are "rad" or "degree." The positive direction of the angle, spatial position, and spatial phase in the present invention may be defined as either left-handed (counterclockwise) or right-handed (clockwise). However, in this specification, in order to maintain the simplicity of the explanation, the positive direction of the angle, spatial position, and spatial phase is defined as counterclockwise rotation, and the present invention will be described. This does not diminish the generality of the invention.

本発明における「空間依存して変化する固定子インダクタンス」とは、「d軸インダクタンス、q軸インダクタンス、d軸とq軸の両インダクタンスから決められるインダクタンスの少なくともいずれか1つのインダクタンスの値が、回転子の空間位相に応じて変化するインダクタンス」を意味する。In the present invention, the term "stator inductance that varies depending on space" means that the value of at least one of the inductances determined from the d-axis inductance, the q-axis inductance, and the inductance of both the d-axis and the q-axis depends on the rotation. Inductance that changes depending on the spatial phase of the child.

空間依存して変化するd軸インダクタンスを「Ld(θα)」で表現し、空間依存して変化するq軸インダクタンスを「Lq(θα)」で表現し、空間依存して変化するd軸インダクタンスの空間微分値を「Ld(θα)」で表現し、空間依存して変化するq軸インダクタンスの空間微分値を「Lq(θα)」で表現するとき、d軸、q軸インダクタンスの空間微分値は、数式による場合、以下のように定義される(回転子位相を意味するθαは、図1を用いて後に詳細説明)。

Figure 0007450150000001
他のインダクタンスに対する「空間微分値」の意味、定義も(1)式と同様である。The d-axis inductance, which changes depending on space, is expressed as "Ld (θα)," and the q-axis inductance, which changes depending on space, is expressed as "Lq (θα)." When the spatial differential value is expressed as “L ~ d (θα)” and the spatial differential value of the q-axis inductance, which changes depending on space, is expressed as “L ~ q (θα),” the d-axis and q-axis inductances are The spatial differential value is defined as follows when using a mathematical formula (θα, which means the rotor phase, will be explained in detail later using FIG. 1).
Figure 0007450150000001
The meaning and definition of "spatial differential value" with respect to other inductances are also the same as in equation (1).

本発明では、原則として、d軸、q軸の直交2軸からなるdq同期座標系(図1を用いて後に詳細説明)において、各軸と関連したパラメータ、物理量等には、脚符d、qを付して各軸との関係を明示する。In the present invention, in principle, in a dq synchronous coordinate system consisting of two orthogonal axes, the d-axis and the q-axis (details will be explained later with reference to FIG. 1), parameters, physical quantities, etc. related to each axis are represented by foot marks d, Add q to clearly indicate the relationship with each axis.

本発明の対象とする同期電動機は、リラクタンストルクを発生する同期電動機である。この種の同期電動機としては、同期リラクタンス電動機、永久磁石同期電動機、巻線形同期電動機、ハイブリッド界磁同期電動機などが存在する。これら同期電動機の中で、通電により、リラクタンストルクのみを発生する電動機は、同期リラクタンス電動機である。他の同期電動機は、通電により、リラクタンストルクに加えて、マグネットトルクを発生する。本発明はリラクタンストルクを発生するすべての同期電動機に適用される。しかしながら、本発明の核心部分の簡明な説明には、リラクタンストルクのみを発生する同期リラクタンスが好適でる。この点を考慮し、以降の説明では、同期電動機として、主として同期リラクタンス電動機を用い、本発明を説明する。A synchronous motor to which the present invention is applied is a synchronous motor that generates reluctance torque. Examples of this type of synchronous motor include a synchronous reluctance motor, a permanent magnet synchronous motor, a wound synchronous motor, and a hybrid field synchronous motor. Among these synchronous motors, a motor that generates only reluctance torque when energized is a synchronous reluctance motor. Other synchronous motors generate magnet torque in addition to reluctance torque when energized. The present invention applies to all synchronous motors that generate reluctance torque. However, for a simple explanation of the core of the invention, a synchronous reluctance system that generates only reluctance torque is preferable. In consideration of this point, in the following description, the present invention will be explained using mainly a synchronous reluctance motor as the synchronous motor.

本発明の対象とする同期電動機の数学モデルに関する先行発明、同数学モデルに立脚した模擬装置、特性解析装置、制御装置としては、例えば、非特許文献1がある。非特許文献1などに用いられた数学モデルの説明に先立って、数学モデルが定義される座標系を説明する。図1を考える。同期リラクタンス電動機の回転子の位相θα、θγとして、同図(a)は回転子の順突極(正突極)位相を採用し、同図(b)は回転子の逆突極(負突極)位相を採用している。両図には、同期リラクタンス電動機のための2軸直交座標系として、電気速度ω2nで回転する回転子に位相差なく同期したdq同期座標系、αβ固定座標系、任意速度ωγで回転するγδ一般座標系の3座標系を描画している。両図とも、基軸から副軸の方向を正方向としている。したがって、副軸は主軸に対してπ/2[rad]位相進みの位置にある。回転子位相θα、θγは、各々α軸、γ軸を基準にしている。γδ一般座標系は、αβ固定座標系、dq同期座標系を特別の場合として包含する、最も一般性の高い座標系である。なお、位相算定の基準であるαβ固定座標系の基軸・α軸の位置は、一般的には、uvw三相巻線のu相巻線の中心位置に選定される。上記説明より既に明白なように、本発明においては、回転子位相と回転子の突極位相とは同義である。Prior inventions related to mathematical models of synchronous motors, which are the object of the present invention, as well as simulation devices, characteristic analysis devices, and control devices based on the mathematical models, include, for example, Non-Patent Document 1. Prior to explaining the mathematical model used in Non-Patent Document 1, etc., the coordinate system in which the mathematical model is defined will be explained. Consider Figure 1. As the phases θα and θγ of the rotor of the synchronous reluctance motor, the figure (a) adopts the forward salient pole (positive salient pole) phase of the rotor, and the figure (b) adopts the rotor's reverse salient pole (negative salient pole) phase. Polar) phase is adopted. In both figures, two axes orthogonal coordinate systems for a synchronous reluctance motor include a dq synchronous coordinate system that is synchronized with the rotor that rotates at an electrical speed of ω2n without a phase difference, an αβ fixed coordinate system, and a general γδ that rotates at an arbitrary speed ωγ. A three-coordinate system is drawn. In both figures, the direction from the base axis to the minor axis is the positive direction. Therefore, the secondary axis is at a phase lead of π/2 [rad] with respect to the main axis. The rotor phases θα and θγ are based on the α axis and the γ axis, respectively. The γδ general coordinate system is the most general coordinate system that includes the αβ fixed coordinate system and the dq synchronous coordinate system as special cases. Note that the position of the base axis/α axis of the αβ fixed coordinate system, which is the reference for phase calculation, is generally selected to be the center position of the u-phase winding of the uvw three-phase winding. As is already clear from the above description, in the present invention, the rotor phase and the salient pole phase of the rotor have the same meaning.

同期電動機は、電気回路であり、トルク発生機であり、電気エネルギーを機械エネルギーへ変換するエネルギー変換機でもある。同期電動機のこの本質に起因して、同期電動機の数学モデルは、厳密には、電気回路としての動的特性を記述した回路方程式(第1基本式)、トルク発生機としてのトルク発生関係を記述したトルク発生式(第2基本式)、エネルギー変換機としての動的関係を記述したエネルギー伝達式(第3基本式)の3基本式より構成される。また、数学モデルが合理的モデルとして工学的意味をもつためには、単一の同期電動機を異なる観点から記述したこれら3基本式は互いに数学的に整合するものでなくてはならない。本整合性は、「自己整合性」と呼ばれる。非特許文献1によれば、同期リラクタンス電動機のγδ一般座標系上の数学モデルとして、次のものが知られている。A synchronous motor is an electric circuit, a torque generator, and an energy converter that converts electrical energy into mechanical energy. Due to this essence of a synchronous motor, the mathematical model of a synchronous motor is strictly a circuit equation (first basic equation) that describes the dynamic characteristics as an electric circuit, and a torque generation relationship as a torque generator. It is composed of three basic equations: a torque generation equation (second basic equation), and an energy transfer equation (third basic equation) that describes the dynamic relationship as an energy converter. Furthermore, in order for the mathematical model to have engineering meaning as a rational model, these three basic equations that describe a single synchronous motor from different viewpoints must be mathematically consistent with each other. This consistency is called "self-consistency." According to Non-Patent Document 1, the following is known as a mathematical model of a synchronous reluctance motor on the γδ general coordinate system.

回路方程式(第1基本式)

Figure 0007450150000002
Figure 0007450150000003
トルク発生式(第2基本式)
Figure 0007450150000004
エネルギー伝達式(第3基本式)
Figure 0007450150000005
Circuit equation (first basic equation)
Figure 0007450150000002
Figure 0007450150000003
Torque generation formula (second basic formula)
Figure 0007450150000004
Energy transfer formula (third basic formula)
Figure 0007450150000005

上記数学モデルにおける2×1ベクトルv1、i1、φiは、γδ一般座標系上で定義された固定子電圧、固定子電流、固定子反作用磁束である。τrは発生トルク(リラクタンストルク)であり、ω2n、ω2mは回転子の電気速度および機械速度である。また、Np、R1、Li、Lmは、極対数、固定子の抵抗、同相インダクタンス、鏡相インダクタンスである。なお、同相、鏡相インダクタンスは、d軸、q軸インダクタンスと次の関係を有する。

Figure 0007450150000006
従前数学モデルにおいては、(5)式等におけるd軸、q軸インダクタンス、同相、鏡相インダクタンスは「一定」であり、ひいては「空間依存性を有しない」ことを前提としている。なお、鏡相インダクタンスLmは、(5a)式より明白なように、回転子位相を順突極位相に選定する場合(図1(a))には正となり、逆突極位相に選定する場合(図1(b))には負となる。The 2×1 vectors v1, i1, and φi in the above mathematical model are the stator voltage, stator current, and stator reaction magnetic flux defined on the γδ general coordinate system. τr is the generated torque (reluctance torque), and ω2n and ω2m are the electrical speed and mechanical speed of the rotor. Further, Np, R1, Li, and Lm are the number of pole pairs, stator resistance, common-mode inductance, and mirror-phase inductance. Note that the in-phase and mirror-phase inductances have the following relationship with the d-axis and q-axis inductances.
Figure 0007450150000006
In the conventional mathematical model, it is assumed that the d-axis, q-axis inductance, in-phase, and mirror-phase inductance in equation (5) etc. are "constant" and, furthermore, "have no spatial dependence." Note that, as is clear from equation (5a), the mirror phase inductance Lm is positive when the rotor phase is selected as the forward salient pole phase (Fig. 1(a)), and when the rotor phase is selected as the reverse salient pole phase. (FIG. 1(b)) is negative.

回転子の位相と速度に関しては、次の関係が成立している。

Figure 0007450150000007
また、記号「s」は微分演算子「d/dt」を意味し、種々の2×2行列は、以下のように定義されている。
Figure 0007450150000008
Figure 0007450150000009
Regarding the phase and speed of the rotor, the following relationship holds true.
Figure 0007450150000007
Further, the symbol "s" means a differential operator "d/dt", and various 2×2 matrices are defined as follows.
Figure 0007450150000008
Figure 0007450150000009

γδ一般座標系上の数学モデルに、αβ固定座標系の条件(θγ=θα、ωγ=0)を付与すると、これはαβ固定座標系上の数学モデルとなる。また、γδ一般座標系上の数学モデルに、dq同期座標系の条件(θγ=0、ωγ=ω2n)を付与すると、これはdq同期座標系上の数学モデルとなる。例えば、(3)式のトルク発生式は、dq同期座標系の条件(θγ=0、ωγ=ω2n)を付与すると、次式となる。

Figure 0007450150000010
ここに、id、iqは各々d軸電流、q軸電流を意味する。When the conditions of the αβ fixed coordinate system (θγ=θα, ωγ=0) are given to a mathematical model on the γδ general coordinate system, this becomes a mathematical model on the αβ fixed coordinate system. Furthermore, when the conditions of the dq synchronous coordinate system (θγ=0, ωγ=ω2n) are given to the mathematical model on the γδ general coordinate system, this becomes a mathematical model on the dq synchronous coordinate system. For example, the torque generation equation of equation (3) becomes the following equation when the conditions of the dq synchronous coordinate system (θγ=0, ωγ=ω2n) are given.
Figure 0007450150000010
Here, id and iq mean a d-axis current and a q-axis current, respectively.

「インダクタンスは一定」とする(8)式は、「d軸電流、q軸電流が一定の場合には、発生のリラクタンストルクは一定である」ことを示している。ところが、「一定のd軸電流、q軸電流に対して、実際の同期リラクタン電動機は空間依存のトルクリップルを発生する」ことが、実験的に知られている。本特性は、リラクタンストルクを発生する他の同期電動機も同様である。これらの事実は、以下を意味してている。(a)従前数学モデルに立脚した模擬装置では、空間依存のリラクタンストルクリプルが模擬できない。(b)従前数学モデルに立脚した特性解析装置では、空間依存のリラクタンストルクリプルを解析できない。(c)従前数学モデルに立脚した制御装置では、空間依存のリラクタンストルクリプルを適切に制御できない。Equation (8), which assumes that "inductance is constant," indicates that "when the d-axis current and the q-axis current are constant, the generated reluctance torque is constant." However, it has been experimentally known that ``actual synchronous reluctance motors generate space-dependent torque ripples for constant d-axis and q-axis currents.'' This characteristic is similar to other synchronous motors that generate reluctance torque. These facts mean the following. (a) Space-dependent reluctance torque ripple cannot be simulated with a simulation device based on conventional mathematical models. (b) Characteristic analysis devices based on conventional mathematical models cannot analyze space-dependent reluctance torque ripple. (c) A control device based on a conventional mathematical model cannot appropriately control space-dependent reluctance torque ripple.

新中新二:「同期モータのベクトル制御技術」、東京電機大学出版局(2119)Shinji Shinnaka: “Vector control technology for synchronous motors”, Tokyo Denki University Press (2119)

本発明は上記背景の下になされたものであり、その目的は、リラクタンストルクを発生する同期電動機のための数学モデルであって、同期電動機が空間依存のリラクタンストルクを発生する場合にも、模擬、特性解析、制御への利用可能な数学モデルを提供することにある。さらには、同モデルに立脚した模擬装置、あるいは特性解析装置、あるいは制御装置を提供することにある。The present invention has been made against the above background, and its purpose is to provide a mathematical model for a synchronous motor that generates reluctance torque, and to simulate the case where the synchronous motor generates a spatially dependent reluctance torque. The objective is to provide a mathematical model that can be used for characterization, analysis, and control. Furthermore, it is an object of the present invention to provide a simulation device, a characteristic analysis device, or a control device based on the same model.

上記目的を達成するために、請求項1の発明は、回転子が固定子側の電力周波数と平均的に同じ電気速度でかつリラクタンストルクを発生しつつ回転する同期電動機の模擬、あるいは特性解析、あるいは制御のための数学モデルであって、該数学モデルが、該リラクタンストルクの発生モデルとして、空間依存して変化するインダクタンスに比例したリラクタンストルクの発生を示す比例リラクタンストルクの発生モデルと、空間依存して変化するインダクタンスの空間微分値に比例したリラクタンストルクの発生を示す微分リラクタンストルクの発生モデルとの両発生モデルを備えることを特徴とする。In order to achieve the above object, the invention of claim 1 provides simulation or characteristic analysis of a synchronous motor in which the rotor rotates at an electrical speed that is on average the same as the power frequency of the stator and while generating reluctance torque. Alternatively, the mathematical model is a mathematical model for control, and the mathematical model includes, as the reluctance torque generation model, a proportional reluctance torque generation model that shows the generation of reluctance torque proportional to the inductance that changes in a spatially dependent manner, and a spatially dependent The present invention is characterized by having both a differential reluctance torque generation model and a differential reluctance torque generation model that shows the generation of reluctance torque proportional to the spatial differential value of inductance that changes.

請求項2の発明は、請求項1記載の同期電動機の模擬、あるいは特性解析、あるいは制御のための数学モデルであって、該固定子の特定箇所から評価した該回転子の突極の位相を「θα」で表現し、突極の位相をd軸の位相とし、d軸に対してπ/2[rad]の位相進みにq軸をもつ2軸直交座標系をdq同期座標系とし、空間依存して変化するd軸インダクタンスを「Ld(θα)」で表現し、空間依存して変化するq軸インダクタンスを「Lq(θα)」で表現し、空間依存して変化するd軸インダクタンスの空間微分値を「Ld(θα)」で表現し、空間依存して変化するq軸インダクタンスの空間微分値を「Lq(θα)」で表現し、d軸電流を「id」で表現し、q軸電流を「iq」で表現し、該同期電動機の極対数を「Np」で表現し、さらには、リラクタンストルクを「τr」で表現し、比例リラクタンストルクを「τrp」で表現し、微分リラクタンストルクを「τrd」で表現するとき、該リラクタンストルクの発生モデルを、dq同期座標系上での評価において、次式

Figure 0007450150000011
とすることを特徴とする。The invention according to claim 2 is a mathematical model for simulating, analyzing characteristics, or controlling the synchronous motor according to claim 1, which calculates the phase of the salient poles of the rotor evaluated from a specific location on the stator. Expressed as "θα", the phase of the salient pole is the phase of the d-axis, and a two-axis orthogonal coordinate system with the q-axis leading the d-axis by π/2 [rad] is a dq-synchronous coordinate system, and the space The d-axis inductance, which changes depending on space, is expressed as "Ld (θα)," and the q-axis inductance, which changes depending on space, is expressed as "Lq (θα)." The differential value is expressed as “L ~ d (θα)”, the spatial differential value of the q-axis inductance that changes depending on space is expressed as “L ~ q (θα)”, and the d-axis current is expressed as “id”. Then, the q-axis current is expressed as "iq", the number of pole pairs of the synchronous motor is expressed as "Np", the reluctance torque is expressed as "τr", and the proportional reluctance torque is expressed as "τrp". , when the differential reluctance torque is expressed as "τrd", the generation model of the reluctance torque is evaluated on the dq synchronous coordinate system as follows:
Figure 0007450150000011
It is characterized by:

請求項3の発明は、請求項1、あるいは請求項2記載の同期電動機の模擬、あるいは特性解析、あるいは制御のための数学モデルであって、該固定子の特定箇所から評価した該回転子の突極の位相を「θα」で表現するとき、空間依存して変化する固定子インダクタンスを、突極位相θαを変数とする三角関数を用いて近似表現したことを特徴とする。The invention according to claim 3 is a mathematical model for simulating, characterizing, or controlling the synchronous motor according to claim 1 or 2, wherein the rotor is evaluated from a specific location on the stator. When the phase of the salient pole is expressed as "θα," the stator inductance, which changes spatially, is approximately expressed using a trigonometric function with the salient pole phase θα as a variable.

請求項4の発明は、回転子が固定子側の電力周波数と平均的に同じ電気速度でかつリラクタンストルクを発生しつつ回転する同期電動機のための、数学モデル立脚形の模擬装置、あるいは特性解析装置、あるいは制御装置であって、立脚の該数学モデルに、空間依存して変化する固定子インダクタンスに比例したリラクタンストルクの発生を示す比例リラクタンストルクの発生モデルと、空間依存して変化する固定子インダクタンスの空間微分値に比例したリラクタンストルクの発生を示す微分リラクタンストルクの発生モデルとの両発生モデルを備えさせたことを特徴とする。The invention of claim 4 is a mathematical model pedestal type simulating device or characteristic analysis for a synchronous motor in which the rotor rotates at the same electrical speed as the stator side power frequency on average and while generating reluctance torque. A device or a control device, which includes a proportional reluctance torque generation model that shows the generation of reluctance torque proportional to a stator inductance that changes depending on space, and a stator that changes depending on space in the mathematical model of the stance leg. The present invention is characterized by having both a differential reluctance torque generation model and a differential reluctance torque generation model showing the generation of reluctance torque proportional to the spatial differential value of inductance.

本発明の効果を説明する。後掲の実施例で詳しく説明するように、請求項1の本発明によれば、空間依存のリラクタンストルクリプルを発生する同期電動機に対して、この模擬、特性解析、制御を可能とする数学モデルを、最も一般性の高いγδ一般座標系上で、さらには自己整合性を有する形で構成できるようになると言う効果が得られるようになる。すなわち、電気回路としての動的特性を記述した回路方程式(第1基本式)、トルク発生機としてのトルク発生関係を記述したトルク発生式(第2基本式)、電気エネルギーを機械エネルギーへ変換するエネルギー変換機としての動的関係を記述したエネルギー伝達式(第3基本式)の3基本式よりなり、さらにはこれら3基本式が互いに数学的に整合した数学モデルが構成できるようになるという効果が得られる。ひいては、空間依存のリラクタンストルクリプルを発生する同期電動機に対して、この模擬、特性解析、制御を合理的に遂行できるようになるという効果が得られる。The effects of the present invention will be explained. As will be explained in detail in the examples below, according to the present invention of claim 1, a mathematical model that enables simulation, characteristic analysis, and control of a synchronous motor that generates space-dependent reluctance torque ripple is provided. , the effect of being able to be constructed on the most general γδ general coordinate system and in a self-consistent manner. In other words, the circuit equation (first basic equation) describes the dynamic characteristics of an electric circuit, the torque generation equation (second basic equation) describes the torque generation relationship as a torque generator, and the converter converts electrical energy into mechanical energy. The effect is that it consists of three basic equations, the energy transfer equation (third basic equation), which describes the dynamic relationship as an energy converter, and furthermore, it becomes possible to construct a mathematical model in which these three basic equations are mathematically consistent with each other. is obtained. As a result, it is possible to rationally perform simulation, characteristic analysis, and control of a synchronous motor that generates space-dependent reluctance torque ripple.

つづいて、請求項2の発明の効果を説明する。請求項1の発明による数学モデルは、γδ一般座標系、αβ固定座標系、dq同期座標系のいずれの座標系の上でも構築可能である。しかしながら、構築の難易度は同一ではない。3座標系の中ではdq同期座標系上での構築がもっとも平易である。ひいては、請求項2の発明によれば、同期電動機の模擬、あるいは特性解析、あるいは制御に供しうる数学モデルを最も簡単に構築できると言う効果が得られる。この結果、請求項2の発明によれば、請求項1の発明の効果を高めることができると言う効果も得られる。Next, the effect of the invention of claim 2 will be explained. The mathematical model according to the invention of claim 1 can be constructed on any of the γδ general coordinate system, αβ fixed coordinate system, and dq synchronous coordinate system. However, the difficulty of construction is not the same. Among the three coordinate systems, construction on the dq synchronous coordinate system is the simplest. Furthermore, according to the second aspect of the invention, it is possible to most easily construct a mathematical model that can be used for simulation, characteristic analysis, or control of a synchronous motor. As a result, according to the invention of claim 2, the effect of increasing the effect of the invention of claim 1 can also be obtained.

つづいて、請求項3の発明の効果を説明する。大半の同期電動機において、リラクタンストルクの空間依存リプル主成分は、回転子電気速度の6次、12次の高調波成分であることが実験的に知られている。主成分を中心とした同期電動機の模擬、特性解析、さらには制御には、空間依存して変化するインダクタンスを三角フーリエ級数展開し、0次、6次、12次の成分で近似表現するのが合理的である。請求項3の発明によれば、突極位相θαを変数とする三角関数を用いてインダクタンスを近似表現できるので、三角フーリエ級数展開を高い数学的合理性をもって適用できるようになる。この結果、請求項3の発明によれば、空間依存リラクタンストルクリプルの主成分を中心とした同期電動機の模擬、特性解析、さらには制御を可能とするコンパクトな数学モデルが構築できると言う効果が得られる。Next, the effects of the invention according to claim 3 will be explained. It is experimentally known that in most synchronous motors, the main spatially dependent ripple components of reluctance torque are the 6th and 12th harmonic components of the rotor electrical speed. For simulation, characteristic analysis, and even control of a synchronous motor centered on its principal components, it is recommended to expand the space-dependent inductance into a triangular Fourier series and approximately express it using 0th, 6th, and 12th order components. Reasonable. According to the third aspect of the invention, since inductance can be approximately expressed using a trigonometric function with the salient pole phase θα as a variable, triangular Fourier series expansion can be applied with high mathematical rationality. As a result, according to the invention of claim 3, it is possible to construct a compact mathematical model that enables simulation, characteristic analysis, and control of a synchronous motor centered on the main component of space-dependent reluctance torque ripple. It will be done.

つづいて、請求項4の発明の効果を説明する。請求項4の発明は、請求項1の発明による数学モデルに立脚した模擬装置、あるいは特性解析装置、あるいは制御装置の構築を可能とするものである。この結果、請求項4の発明によれば、請求項1の発明の効果を、当該の模擬装置、あるいは特性解析装置、あるいは制御装置に反映さすことができると言う効果を得ることができる。すなわち、請求項4の発明によれば、リラクタンストルクリプルを発生する同期電動機のための模擬装置、あるいは特性解析装置、あるいは制御装置を構築できると言う効果を得ることができる。Next, the effect of the invention of claim 4 will be explained. The invention of claim 4 makes it possible to construct a simulation device, a characteristic analysis device, or a control device based on the mathematical model according to the invention of claim 1. As a result, according to the invention of claim 4, it is possible to obtain the effect that the effect of the invention of claim 1 can be reflected in the simulation device, characteristic analysis device, or control device. That is, according to the fourth aspect of the present invention, it is possible to construct a simulation device, a characteristic analysis device, or a control device for a synchronous motor that generates reluctance torque ripple.

「3種の2軸直交座標系の関係を示す図」 "Diagram showing the relationship between three types of two-axis orthogonal coordinate systems" 「本発明によるγδ一般座標系上の模擬装置あるいは特性解析装置の構成例を示す図」 "A diagram showing an example of the configuration of a simulating device or a characteristic analysis device on the γδ general coordinate system according to the present invention" 「本発明による制御装置を用いた駆動システムの構成例を示す図」 "Diagram showing an example of the configuration of a drive system using a control device according to the present invention"

以下、図面を用いて、本発明の好適な実施態様を具体的に説明する。Hereinafter, preferred embodiments of the present invention will be specifically described using the drawings.

同期リラクタンス電動機に対して、請求項1の発明に基づくγδ一般座標系上の数学モデルを以下に示す。A mathematical model on a γδ general coordinate system based on the invention of claim 1 for a synchronous reluctance electric motor is shown below.

回路方程式(第1基本式)

Figure 0007450150000012
Figure 0007450150000013
トルク発生式(第2基本式)
Figure 0007450150000014
エネルギー伝達式(第3基本式)
Figure 0007450150000015
Circuit equation (first basic equation)
Figure 0007450150000012
Figure 0007450150000013
Torque generation formula (second basic formula)
Figure 0007450150000014
Energy transfer formula (third basic formula)
Figure 0007450150000015

(10)式の回路方程式、(11)式のトルク発生式においては、同相インダクタンス、鏡相インダクタンスが空間依存性をもつものとして、回転子位相θαの関数としてこれらを表現している。空間依存性をもつ同相、鏡相インダクタンスは、空間依存性をもつd軸、q軸インダクタンスと次の関係をもつ。

Figure 0007450150000016
「Li(θα)」、「Lm(θα)」は、空間依存性をもつ同相、鏡相インダクタンスの空間微分値を意味し、数式では、以下のように定義される((1)式参照)。
Figure 0007450150000017
(14)式の第2式は、空間依存性をもつd軸、q軸インダクタンスの空間微分値との関係を示している。なお、(10)~(12)式の数学モデルにおいても、回転子の位相と速度に関する(6)式、2×2行列の定義に関する(7)式は適用される。In the circuit equation of equation (10) and the torque generation equation of equation (11), the in-phase inductance and mirror-phase inductance are assumed to have spatial dependence, and are expressed as a function of the rotor phase θα. The spatially dependent in-phase and mirror-phase inductances have the following relationship with the spatially dependent d-axis and q-axis inductances.
Figure 0007450150000016
“L ~ i (θα)” and “L ~ m (θα)” mean the spatial differential values of the in-phase and mirror-phase inductances that have spatial dependence, and are defined in the formula as follows ((1 ) formula).
Figure 0007450150000017
The second equation of equation (14) shows the relationship with the spatial differential values of the d-axis and q-axis inductances, which have spatial dependence. Note that equation (6) regarding the phase and speed of the rotor and equation (7) regarding the definition of the 2×2 matrix are also applied to the mathematical models of equations (10) to (12).

(11)式のトルク発生式においては、(11b)式が、空間依存して変化するインダクタンスに比例したリラクタンストルクの発生を示す比例リラクタンストルクの発生モデルである。(11b)式の第2式が、インダクタンス(特に、鏡相インダクタンス)に比例した比例リラクタンストルクの発生を明瞭に示している。また、(11c)式が、空間依存して変化するインダクタンスの空間微分値に比例したリラクタンストルクの発生を示す微分リラクタンストルクの発生モデルである。(11c)式の第2式が、インダクタンス(同相インダクタンス、鏡相インダクタンス)の空間微分値に比例したリラクタンストルクの発生を示す微分リラクタンストルクの発生を明瞭に示している。(11a)式は、「全リラクタンストルクは、比例リラクタンストルクと微分リラクタンストルクの総和である」ことを示している。これは、請求項1の本発明に正確に従っている。In the torque generation equation (11), equation (11b) is a proportional reluctance torque generation model that indicates the generation of reluctance torque proportional to the inductance that changes depending on space. The second equation (11b) clearly shows the generation of proportional reluctance torque proportional to inductance (especially mirror inductance). Further, equation (11c) is a differential reluctance torque generation model that indicates the generation of reluctance torque that is proportional to the spatial differential value of inductance that changes depending on space. The second equation (11c) clearly shows the generation of differential reluctance torque that is proportional to the spatial differential value of inductance (in-phase inductance, mirror-phase inductance). Equation (11a) indicates that "the total reluctance torque is the sum of the proportional reluctance torque and the differential reluctance torque." This is exactly in accordance with the invention of claim 1.

(10a)式の回路方程式の物理的意味は、従前の(2a)式と同様であり、当業者には周知であるので、この説明は省略する。(12)式のエネルギー伝達式においては、同式左辺は同期電動機への瞬時印加電力を意味し、同式右辺は電動機内部の瞬時電力の様子を意味している。すなわち、右辺第1項は瞬時銅損を、第2項はインダクタンスに蓄積された磁気エネルギーの瞬時変化を、第3項は回転子から出力される瞬時機械的電力を意味している。The physical meaning of the circuit equation (10a) is the same as that of the previous equation (2a) and is well known to those skilled in the art, so a description thereof will be omitted. In the energy transfer equation (12), the left side of the equation means the instantaneous power applied to the synchronous motor, and the right side of the equation means the state of the instantaneous power inside the motor. That is, the first term on the right side means instantaneous copper loss, the second term means instantaneous change in magnetic energy stored in the inductance, and the third term means instantaneous mechanical power output from the rotor.

(11)式のトルク発生式に用いた空間依存の同相インダクタンス、鏡相インダクタンスは、(10)式の回路方程式に用いたものと同一である。さらには、(11)式のトルク発生式に用いた固定子反作用磁束は、(10)式の回路方程式に用いたもの、(12)式のエネルギー伝達式に用いたものと同一である。また、(12)式のエネルギー伝達式の右辺第3項におけるリラクタンストルクτrは、(11)式のトルク発生式に規定した全リラクタンストルクτrである。本事実は、請求項1の発明に基づく(11)式のトルク発生式は、(10)式の回路方程式と(12)式のエネルギー伝達式と数学的に整合していること、3基本式からなる(10)~(12)式の数学モデルは自己整合性を有することを意味する。以上より明白なように、請求項1の本発明によれば、γδ一般座標系上で、自己整合性を備えた数学モデル(特に、リラクタンストルクリプルを模擬、解析等ができる数学モデル)を構築することができる。3基本式の自己整合性を備えた数学モデルを得るには、請求項1の発明による比例リラクタンストルクと微分リラクタンストルクとが同時に必要である。The space-dependent common-mode inductance and mirror-phase inductance used in the torque generation equation of equation (11) are the same as those used for the circuit equation of equation (10). Furthermore, the stator reaction magnetic flux used in the torque generation equation (11) is the same as that used in the circuit equation (10) and the energy transfer equation (12). Further, the reluctance torque τr in the third term on the right side of the energy transfer equation (12) is the total reluctance torque τr defined in the torque generation equation (11). This fact indicates that the torque generation equation of equation (11) based on the invention of claim 1 is mathematically consistent with the circuit equation of equation (10) and the energy transfer equation of equation (12), and that the three basic equations This means that the mathematical model of equations (10) to (12) consisting of is self-consistent. As is clear from the above, according to the present invention of claim 1, a self-consistent mathematical model (particularly a mathematical model that can simulate and analyze reluctance torque ripple) is constructed on the γδ general coordinate system. be able to. In order to obtain a mathematical model with self-consistency of the three basic equations, the proportional reluctance torque and the differential reluctance torque according to the invention of claim 1 are required at the same time.

(10)~(12)式のγδ一般座標系上の数学モデルに、αβ固定座標系の条件(θγ=θα、ωγ=0)を付与すると、以下に示すαβ固定座標系上の数学モデルを得ることができる。When the conditions of the αβ fixed coordinate system (θγ=θα, ωγ=0) are given to the mathematical model on the γδ general coordinate system of equations (10) to (12), the following mathematical model on the αβ fixed coordinate system can be obtained. Obtainable.

回路方程式(第1基本式)

Figure 0007450150000018
トルク発生式(第2基本式)
Figure 0007450150000019
エネルギー伝達式(第3基本式)
Figure 0007450150000020
Circuit equation (first basic equation)
Figure 0007450150000018
Torque generation formula (second basic formula)
Figure 0007450150000019
Energy transfer formula (third basic formula)
Figure 0007450150000020

請求項1の発明に基づくリラクタンストルクの発生式である(16)式が、他の基本式である回路方程式、エネルギー伝達式と整合していることは、γδ一般座標系上の数学モデルと同様である。特に、αβ固定座標系上では、エネルギー伝達式における磁気エネルギーとトルク発生式における全リラクタンストルクの間には、θmを回転子の機械位相とするとき、磁気エネルギーの空間微分が全リラクタンストルクに等しいとする次の関係が成立している。

Figure 0007450150000021
以上より明白なように、請求項1の本発明によれば、αβ固定座標系上で、自己整合性を備えた数学モデル(特に、リラクタンストルクリプルを模擬、解析等ができる数学モデル)を構築することができる。αβ固定座標系上で、3基本式の自己整合性を備えた数学モデルを得るには、請求項1の発明による比例リラクタンストルクと微分リラクタンストルクとが同時に必要である。The fact that equation (16), which is the reluctance torque generation equation based on the invention of claim 1, is consistent with the other basic equations, the circuit equation and the energy transfer equation, is similar to the mathematical model on the γδ general coordinate system. It is. In particular, on the αβ fixed coordinate system, between the magnetic energy in the energy transfer equation and the total reluctance torque in the torque generation equation, the spatial differential of the magnetic energy is equal to the total reluctance torque, where θm is the mechanical phase of the rotor. The following relationship holds true.
Figure 0007450150000021
As is clear from the above, according to the present invention of claim 1, a self-consistent mathematical model (particularly a mathematical model that can simulate and analyze reluctance torque ripple) is constructed on an αβ fixed coordinate system. be able to. In order to obtain a mathematical model with self-consistency of the three basic equations on the αβ fixed coordinate system, the proportional reluctance torque and the differential reluctance torque according to the invention of claim 1 are required at the same time.

(10)~(12)式のγδ一般座標系上の数学モデルに、dq同期座標系の条件(θγ=0、ωγ=ω2n)を付与すると、以下に示すdq同期座標系上の数学モデルを得ることができる。これは、請求項1と請求項2の発明に基づく数学モデルでもある。When the conditions of the dq synchronous coordinate system (θγ = 0, ωγ = ω2n) are given to the mathematical model on the γδ general coordinate system of equations (10) to (12), the following mathematical model on the dq synchronous coordinate system can be obtained. Obtainable. This is also a mathematical model based on the invention of claims 1 and 2.

回路方程式(第1基本式)

Figure 0007450150000022
トルク発生式(第2基本式)
Figure 0007450150000023
Figure 0007450150000024
エネルギー伝達式(第3基本式)
Figure 0007450150000025
Circuit equation (first basic equation)
Figure 0007450150000022
Torque generation formula (second basic formula)
Figure 0007450150000023
Figure 0007450150000024
Energy transfer formula (third basic formula)
Figure 0007450150000025

請求項1、請求項2の発明に基づくリラクタンストルクの発生式である(20)式は、請求項2の発明によるリラクタンストルクそのもの、すなわち(9)式と同一である。(20)式が、他の基本式である回路方程式、エネルギー伝達式と整合していることは、γδ一般座標系上の数学モデルと同様である。以上より明白なように、請求項1、請求項2の発明によれば、dq同期座標系上で、自己整合性を備えた数学モデル(特に、リラクタンストルクリプルを模擬、解析等ができる数学モデル)を構築することができる。dq同期座標系上で、3基本式の自己整合性を備えた数学モデルを得るには、請求項1、請求項2の発明による比例リラクタンストルクと微分リラクタンストルクとが同時に必要とされる。Equation (20), which is the reluctance torque generation equation based on the invention of claims 1 and 2, is the same as the reluctance torque itself, ie, equation (9), according to the invention of claim 2. The fact that equation (20) is consistent with other basic equations, such as the circuit equation and the energy transfer equation, is similar to the mathematical model on the γδ general coordinate system. As is clear from the above, according to the inventions of claims 1 and 2, a mathematical model (particularly a mathematical model that can simulate and analyze reluctance torque ripple) that has self-consistency on a dq synchronous coordinate system. can be constructed. In order to obtain a mathematical model with self-consistency of the three basic equations on the dq synchronous coordinate system, the proportional reluctance torque and the differential reluctance torque according to the inventions of claims 1 and 2 are required at the same time.

当業者は周知のように、リラクタンストルクを発生する同期電動機においては、空間依存リラクタンストルクリプルの主成分は6次、12次成分である。トルクリプルの主要成分を容易に再現できるようにするには、請求項3の発明を適用するばよい。すなわち、請求項3の発明に基づき、空間依存して変化するインダクタンスを、突極位相θαを変数とする三角関数を用いて、例えば次式のように、近似表現すればよい。

Figure 0007450150000026
上式におけるLdfはd軸インダクタンスの空間平均値を、また、wd6c、wd6sはd軸インダクタンスの6次高調波成分(余弦成分、正弦成分)振幅の平均値Ldfに対する相対比を意味し、wd12c、wd12sはd軸インダクタンスの12次高調波成分(余弦成分、正弦成分)振幅の平均値Ldfに対する相対比を意味する。q軸インダクタンスに関しても、同様である。As is well known to those skilled in the art, in a synchronous motor that generates reluctance torque, the main components of the space-dependent reluctance torque ripple are the 6th and 12th order components. In order to easily reproduce the main components of torque ripple, the invention of claim 3 may be applied. That is, based on the invention of claim 3, the inductance that changes depending on space may be approximately expressed using a trigonometric function using the salient pole phase θα as a variable, for example, as in the following equation.
Figure 0007450150000026
In the above formula, Ldf is the spatial average value of the d-axis inductance, and wd6c and wd6s are the relative ratios of the amplitudes of the 6th harmonic components (cosine component and sine component) of the d-axis inductance to the average value Ldf, and wd12c, wd12s means the relative ratio of the amplitude of the 12th harmonic component (cosine component, sine component) of the d-axis inductance to the average value Ldf. The same applies to the q-axis inductance.

(22)式に対応したインダクタンス空間微分値は,次式となる。

Figure 0007450150000027
Figure 0007450150000028
The inductance spatial differential value corresponding to equation (22) is as follows.
Figure 0007450150000027
Figure 0007450150000028

(22)式、(23)式では、空間依存して変化するインダクタンスとして、d軸、q軸インダクタンスを選定し、これを、突極位相θαを変数とする三角関数を用いて近似表現した。これに代わって、空間依存して変化するインダクタンスとして、同相、鏡相インダクタンスを選定し、これを、突極位相θαを変数とする三角関数を用いて近似表現することもできる。(22)式を(13a)式に、(23)式を(14)式第2式に用いれば、所期の三角関数を用いた近似表現を得ることができる。In equations (22) and (23), the d-axis and q-axis inductances are selected as the inductances that change depending on space, and these are approximately expressed using trigonometric functions with the salient pole phase θα as a variable. Alternatively, in-phase or mirror-phase inductance may be selected as the inductance that changes depending on space, and this may be approximately expressed using a trigonometric function using the salient pole phase θα as a variable. By using equation (22) as equation (13a) and using equation (23) as equation (14), it is possible to obtain an approximate expression using the desired trigonometric functions.

請求項4の発明に(11)式に加え(10)式を適用すると、同期リラクタンス電動機のためのγδ一般座標系上の模擬装置あるいは特性解析装置として、図2を得る。同図(a)は、装置の全体構成を、同図(b)は、(7a)式に定義した2×2D行列の逆行列実現の1例を示している。図2では、図中における太い信号線は2×1のベクトル信号を意味する。また,「1/s」ブ

Figure 0007450150000029
の場合にはスカラ信号によるベクトルの各要素との乗算を実行するベクトル乗算器を,入力が2個のベクトル信号の場合には内積演算を遂行し結果をスカラ信号として出力する内積器を意味するものとしている。なお,図2では,図の輻輳を避けるため,加算器への入力信号の極性は,正の場合は極性記述を省略し,負の場合のみ極性反転記号「-」を付している。また,同様の理由により,関連ブロックのγδ一般座標系の速度ωγ,同じく回転子位相θγ、θαへの依存性は,貫徹矢印で表現している。図2の模擬装置あるいは特性解析装置においては、比例リラクタンストルク、微分リラクタンストルク、全リラクタンストルクを、一点鎖線の矢印で示した。同図では、比例リラクタンストルクと微分リラクタンストルクに用いられる極対数Npを共有するコンパクトな形式を採用している。なお、全リラクタンストルク(比例リラクタンストルクと微分リラクタンストルクの和)が入力印加されているブロックは、機械負荷系を示している。本発明では、機械負荷系は主眼ではないので、慣性モーメントJmと粘性摩擦係数Dmを用いた簡易な形式で示している。When formula (10) is applied in addition to formula (11) to the invention of claim 4, FIG. 2 is obtained as a simulation device or characteristic analysis device on the γδ general coordinate system for a synchronous reluctance motor. FIG. 5(a) shows the overall configuration of the device, and FIG. 1(b) shows an example of realizing an inverse matrix of the 2×2D matrix defined by equation (7a). In FIG. 2, thick signal lines in the figure mean 2×1 vector signals. Also, the “1/s” block
Figure 0007450150000029
In this case, it means a vector multiplier that performs multiplication with each element of a vector by a scalar signal, and when the input is two vector signals, it means an inner product that performs an inner product operation and outputs the result as a scalar signal. I take it as a thing. In FIG. 2, in order to avoid congestion in the diagram, the polarity description of the input signal to the adder is omitted when it is positive, and a polarity inversion symbol "-" is added only when it is negative. Furthermore, for the same reason, the dependence of the related blocks on the speed ωγ of the γδ general coordinate system and the rotor phases θγ and θα is expressed by penetrating arrows. In the simulation device or characteristic analysis device shown in FIG. 2, the proportional reluctance torque, differential reluctance torque, and total reluctance torque are indicated by dashed-dotted arrows. In the figure, a compact format is adopted in which the number Np of pole pairs used for proportional reluctance torque and differential reluctance torque is shared. Note that the block to which the total reluctance torque (the sum of the proportional reluctance torque and the differential reluctance torque) is inputted indicates a mechanical load system. In the present invention, since the mechanical load system is not the main focus, it is shown in a simple format using the moment of inertia Jm and the coefficient of viscous friction Dm.

図2の模擬装置あるいは特性解析装置によれば、印加された固定子電圧v1に対し、固定子電流i1、固定子電流と関係した固定子反作用磁束φiの瞬時値を模擬・解析できる。また、発生リラクタンストルク(比例リラクタンストルク、微分リラクタンストルク、全リラクタンストルク)の瞬時値も模擬・解析できる。さらには、模擬・解析されたこれの物理量を(12)式に適用すると、電動機印加の入力電力、電動機内部の銅損、磁気エネルギー、磁気エネルギーの時間微分、機械的電力の瞬時値を模擬・解析できる。According to the simulation device or characteristic analysis device of FIG. 2, it is possible to simulate and analyze the stator current i1 and the instantaneous value of the stator reaction magnetic flux φi related to the stator current with respect to the applied stator voltage v1. It is also possible to simulate and analyze the instantaneous value of the generated reluctance torque (proportional reluctance torque, differential reluctance torque, total reluctance torque). Furthermore, by applying the simulated and analyzed physical quantities to equation (12), the input power applied to the motor, the copper loss inside the motor, the magnetic energy, the time derivative of the magnetic energy, and the instantaneous value of the mechanical power can be simulated and analyzed. Can be analyzed.

図2の模擬装置あるいは特性解析装置に、αβ固定座標系の条件(θγ=θα、ωγ=0)を付与すると、これはαβ固定座標系上で定義された電圧、電流、磁束を用いた模擬装置あるいは特性解析装置となる。同様に、図2の模擬装置あるいは特性解析装置に、dq同期座標系の条件(θγ=0、ωγ=ω2n)を付与すると、これはdq同期座標系上の電圧、電流、磁束を用いた模擬装置あるいは特性解析装置となる。When the conditions of the αβ fixed coordinate system (θγ = θα, ωγ = 0) are given to the simulation device or characteristic analysis device in Figure 2, this is a simulation using the voltage, current, and magnetic flux defined on the αβ fixed coordinate system. It becomes a device or a characteristic analysis device. Similarly, if the conditions of the dq synchronous coordinate system (θγ = 0, ωγ = ω2n) are given to the simulation device or characteristic analysis device in Fig. 2, this will result in a simulation using the voltage, current, and magnetic flux on the dq synchronous coordinate system. It becomes a device or a characteristic analysis device.

請求項4の発明に基づく制御装置の実施形態の1例を示す。請求項4の発明に基づく制御装置を利用した駆動システムの1例を図3に示した。図3は、制御装置を含む駆動システムの全体を描画している。駆動システムは、大きくは、同期電動機1、電力変換装置2(破線ブロック表示)、制御装置3(破線ブロック表示)から構成されている。電力変換装置は、電力変換器21、電流検出器22から構成されている。制御装置3は、大きくは、指令変換器31と、電流制御部32(破線ブロック表示)から構成されている。電流制御部32では、電流制御が遂行できるように、3相2相変換器321a、2相3相変換器321b、ベクトル回転器322a、322b、電流制御器323を所有している。電流制御部は、さらに、回転子位相を検出のための位相検出器(レゾルバ、エンコーダ等)324、位相検出器で得た回転子位相を近似微分処理して速度を検出する速度検出器325、位相検出器で得た回転子位相を処理してこの余弦信号と正弦信号を発生する余弦正弦信号発生器326も所有している。なお、同図では、簡明のため、複数のスカラ信号を1つのベクトル信号として捉え、複数のスカラ信号線を1本の太い信号線で表現している。三相二相変換器、二相三相変換器から右側に存在する三相信号(すなわち、3×1ベクトル信号)、同じく、左側に存在する二相信号(すなわち、2×1ベクトル信号)は、1本の太い信号線で表現している。なお、ベクトル信号の脚符t、r、sは、各々、uvw座標系上の信号、αβ固定座標系上の信号、dq同期座標系上の信号であることを示している。本実施例の核心は、指令変換器31にある。本器を除く諸機器は基本的に先行発明のものと同様であるので、これら諸機器の説明は省略する。次に、本実施例の核心である指令変換器について説明する。An example of an embodiment of a control device based on the invention of claim 4 is shown. An example of a drive system using a control device based on the invention of claim 4 is shown in FIG. FIG. 3 depicts the entire drive system including the control device. The drive system is broadly composed of a synchronous motor 1, a power converter 2 (represented by a broken line block), and a control device 3 (represented by a broken line block). The power conversion device includes a power converter 21 and a current detector 22. The control device 3 is mainly composed of a command converter 31 and a current control section 32 (indicated by a broken line block). The current control unit 32 includes a three-phase two-phase converter 321a, a two-phase three-phase converter 321b, vector rotators 322a and 322b, and a current controller 323 to perform current control. The current control unit further includes a phase detector (resolver, encoder, etc.) 324 for detecting the rotor phase, a speed detector 325 that performs approximate differential processing on the rotor phase obtained by the phase detector, and detects the speed. It also has a cosine and sine signal generator 326 that processes the rotor phase obtained by the phase detector and generates the cosine and sine signals. Note that in the figure, for simplicity, a plurality of scalar signals are treated as one vector signal, and a plurality of scalar signal lines are expressed as one thick signal line. From a three-phase two-phase converter, the three-phase signal (i.e., 3 × 1 vector signal) present on the right side from the two-phase three-phase converter, and the two-phase signal (i.e., 2 × 1 vector signal) present on the left side are , is expressed by one thick signal line. Note that the foot marks t, r, and s of the vector signal indicate that the vector signal is a signal on the uvw coordinate system, a signal on the αβ fixed coordinate system, and a signal on the dq synchronous coordinate system, respectively. The core of this embodiment lies in the command converter 31. Since the various devices other than this device are basically the same as those of the prior invention, explanations of these devices will be omitted. Next, the command converter, which is the core of this embodiment, will be explained.

指令変換器31は、トルク指令値τr*からdq同期座標系上の電流指令値i1r*(2×1ベクトル)を生成する役割を担っている。例えば、請求項2の発明による(9)式、(20)式のトルク発生式が示すように、d軸、q軸電流が一定の場合、これに対応したリラクタンストルクは、空間依存(回転子位相θαに応じた)のトルクリプルを発生する。請求項2の発明による(9)式、(20)式のトルク発生式は、同時に、「空間依存した形でd軸、q軸電流を通電する場合には、発生リラクタンストルクの平均値はトルク指令値に合致させた上で、発生リラクタンストルクにおけるリラクタンストルクリプルを低減あるいは消滅できる」ことを示している。図3の指令変換器はこの例を示したものである。これを裏付けるように、指令変換器には、トルク指令値τr*とともに、本目的達成に必須の回転子位相θαが入力されている。本実施例における指令変換器の具体的な構成に必要なモデルが、(9)式、(20)式のトルク発生式である。The command converter 31 has the role of generating a current command value i1r* (2×1 vector) on the dq synchronous coordinate system from the torque command value τr*. For example, as shown in the torque generation equations (9) and (20) according to the invention of claim 2, when the d-axis and q-axis currents are constant, the corresponding reluctance torque is spatially dependent (rotor A torque ripple corresponding to the phase θα is generated. The torque generation equations of equations (9) and (20) according to the invention of claim 2 simultaneously state that ``When d-axis and q-axis currents are applied in a space-dependent manner, the average value of the generated reluctance torque is This shows that the reluctance torque ripple in the generated reluctance torque can be reduced or eliminated after matching the command value. The command converter in FIG. 3 illustrates this example. To support this, the rotor phase θα, which is essential for achieving this objective, is input to the command converter together with the torque command value τr*. The models required for the specific configuration of the command converter in this embodiment are the torque generation equations of equations (9) and (20).

請求項2の発明による(9)式、(20)式のトルク発生式は、空間依存のリラクタンストルクリプル低減に加え、「エネルギー伝達式に含まれる銅損を最小化できる、高い自由度の固定子電流が存在する」ことも示している。(9)式、(20)式のトルク発生式には、(21)式のエネルギー伝達式が対応する。指令変換器31を、数学モデルを構成するトルク発生式(第2基本式)を中核に、銅損最小化をも考慮して、トルク指令値から電流指令値を生成するように構成すれば、空間依存のトルクリプルの消滅・低減に加え、銅損の最小化・低減化をも図ることができる。The torque generation equations of equations (9) and (20) according to the invention of claim 2 not only reduce space-dependent reluctance torque ripple, but also provide a stator with a high degree of freedom that can minimize the copper loss included in the energy transfer equation. It also shows that there is an electric current. The torque generation equations of equations (9) and (20) correspond to the energy transfer equation of equation (21). If the command converter 31 is configured to generate the current command value from the torque command value, with the torque generation formula (second basic formula) that constitutes the mathematical model as the core, and also taking copper loss minimization into consideration, In addition to eliminating and reducing space-dependent torque ripple, it is also possible to minimize and reduce copper loss.

実施例1~8においては、リラクタンストルクを発生する同期電動機として同期リラクタン電動機を用いた。請求項1~4の本発明は、同期リラクタンス電動機に限定されるものではなく、リラクタンストルクを発生する同期電動機(永久磁石同期電動機、巻線形同期電動機、ハイブリッド界磁同期電動機など)には、同様に適用されることを指摘しておく。適用可能性を裏付ける証として、「磁気エネルギーの空間微分が全リラクタンストルクに等しい」とする(18)式の関係が、これら同期電動機においても同様に成立することを指摘しておく。In Examples 1 to 8, a synchronous reluctance motor was used as the synchronous motor that generates reluctance torque. The present invention according to claims 1 to 4 is not limited to synchronous reluctance motors, but can similarly apply to synchronous motors that generate reluctance torque (permanent magnet synchronous motors, wound-rotor synchronous motors, hybrid field synchronous motors, etc.). It should be pointed out that this applies to As proof of its applicability, it should be pointed out that the relationship in equation (18), which states that "the spatial differential of magnetic energy is equal to the total reluctance torque," holds true in these synchronous motors as well.

本発明は、リラクタンストルクを活用する同期電動機の模擬、特性解析、制御に、さらにはこのための模擬装置、特性解析装置、制御装置に好適である。INDUSTRIAL APPLICABILITY The present invention is suitable for simulation, characteristic analysis, and control of a synchronous motor that utilizes reluctance torque, and furthermore for a simulation device, a characteristic analysis device, and a control device for this purpose.

1 同期電動機
2 電力変換装置
21 電力変換器
22 電流検出器
3 制御装置
31 指令変換器
32 電流制御部
321a 三相二相変換器
321b 二相三相変換器
322a ベクトル回転器
322b ベクトル回転器
323 電流制御器
324 位相検出器
325 速度検出器
326 余弦正弦信号発生器
1 Synchronous motor 2 Power converter 21 Power converter 22 Current detector 3 Control device 31 Command converter 32 Current control section 321a Three-phase two-phase converter 321b Two-phase three-phase converter 322a Vector rotator 322b Vector rotator 323 Current Controller 324 Phase detector 325 Speed detector 326 Cosine sine signal generator

Claims (2)

回転子が固定子側の電力周波数と平均的に同じ電気速度でかつリラクタンストルクを発生しつつ回転する同期電動機の模擬、あるいは特性解析、あるいは制御のための数学モデルであって、
該数学モデルが、該リラクタンストルクの発生モデルとして、
空間依存して変化する固定子インダクタンスに比例したリラクタンストルクの発生を示す比例リラクタンストルクの発生モデルと、
空間依存して変化する固定子インダクタンスの空間微分値に比例したリラクタンストルクの発生を示す微分リラクタンストルクの発生モデルとの
両発生モデルを備え、
このうえで、
該固定子の特定箇所から評価した該回転子の突極の位相を「θα」で表現し、
突極の位相をd軸の位相とし、d軸に対してπ/2[rad]の位相進みにq軸をもつ2軸直交座標系をdq同期座標系とし、
空間依存して変化するd軸固定子インダクタンスを「Ld(θα)」で表現し、空間依存して変化するq軸固定子インダクタンスを「Lq(θα)」で表現し、
空間依存して変化するd軸固定子インダクタンスの空間微分値を「Ld(θα)」で表現し、
空間依存して変化するq軸固定子インダクタンスの空間微分値を「Lq(θα)」で表現し、
d軸固定子電流を「id」で表現し、q軸固定子電流を「iq」で表現し、
該同期電動機の極対数を「Np」で表現し、
さらには、リラクタンストルクを「τr」で表現し、比例リラクタンストルクを「τrp」で表現し、微分リラクタンストルクを「τrd」で表現するとき、
該リラクタンストルクの発生モデルを、dq同期座標系上での評価において、次式
Figure 0007450150000030
とすることを特徴とする同期電動機の模擬、あるいは特性解析、あるいは制御のための数学モデル。
A mathematical model for simulation, characteristic analysis, or control of a synchronous motor in which the rotor rotates at the same electrical speed as the stator side power frequency on average and while generating reluctance torque,
The mathematical model is a generation model of the reluctance torque,
A proportional reluctance torque generation model that shows the generation of reluctance torque proportional to the stator inductance that changes depending on space;
Equipped with both a differential reluctance torque generation model that shows the generation of reluctance torque proportional to the spatial differential value of stator inductance that changes depending on space,
After this,
The phase of the salient poles of the rotor evaluated from a specific location on the stator is expressed as "θα",
Let the phase of the salient pole be the phase of the d-axis, and let the dq-synchronous coordinate system be a two-axis orthogonal coordinate system with the q-axis having a phase lead of π/2 [rad] with respect to the d-axis,
The d-axis stator inductance, which varies depending on space, is expressed as "Ld (θα)," and the q-axis stator inductance, which varies depending on space, is expressed as "Lq (θα)."
The spatial differential value of the d-axis stator inductance, which changes depending on space, is expressed as “L ~ d(θα)”,
The spatial differential value of the q-axis stator inductance, which changes depending on space, is expressed as "L ~ q (θα)",
The d-axis stator current is expressed as "id", the q-axis stator current is expressed as "iq",
The number of pole pairs of the synchronous motor is expressed as "Np",
Furthermore, when the reluctance torque is expressed by "τr", the proportional reluctance torque is expressed by "τrp", and the differential reluctance torque is expressed by "τrd",
In evaluating the reluctance torque generation model on the dq synchronous coordinate system, the following formula is used:
Figure 0007450150000030
A mathematical model for simulation, characteristic analysis, or control of a synchronous motor, characterized by:
該固定子の特定箇所から評価した該回転子の突極の位相を「θα」で表現するとき、
空間依存して変化する固定子インダクタンスを、突極位相θαを変数とする三角関数を用いて近似表現したことを特徴とする請求項記載の同期電動機の模擬、あるいは特性解析、あるいは制御のための数学モデル。
When the phase of the salient poles of the rotor evaluated from a specific location on the stator is expressed as "θα",
2. A method for simulating, analyzing characteristics, or controlling a synchronous motor according to claim 1 , wherein the stator inductance that changes depending on space is approximately expressed using a trigonometric function with a salient pole phase θα as a variable. mathematical model.
JP2020025695A 2020-01-30 2020-01-30 A mathematical model of a synchronous motor that generates reluctance torque, and a simulation, characteristic analysis, and control device based on the same model. Active JP7450150B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020025695A JP7450150B2 (en) 2020-01-30 2020-01-30 A mathematical model of a synchronous motor that generates reluctance torque, and a simulation, characteristic analysis, and control device based on the same model.

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020025695A JP7450150B2 (en) 2020-01-30 2020-01-30 A mathematical model of a synchronous motor that generates reluctance torque, and a simulation, characteristic analysis, and control device based on the same model.

Publications (2)

Publication Number Publication Date
JP2021121163A JP2021121163A (en) 2021-08-19
JP7450150B2 true JP7450150B2 (en) 2024-03-15

Family

ID=77270021

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020025695A Active JP7450150B2 (en) 2020-01-30 2020-01-30 A mathematical model of a synchronous motor that generates reluctance torque, and a simulation, characteristic analysis, and control device based on the same model.

Country Status (1)

Country Link
JP (1) JP7450150B2 (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001359262A (en) 2000-06-13 2001-12-26 Hitachi Ltd Reluctance motor, rotor thereof, and member thereof
JP2008295200A (en) 2007-05-24 2008-12-04 Aisin Seiki Co Ltd Synchronous motor control device and method for optimizing synchronous motor control
JP2013039020A (en) 2011-07-13 2013-02-21 Suri-Ai:Kk Transverse flux machine

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001359262A (en) 2000-06-13 2001-12-26 Hitachi Ltd Reluctance motor, rotor thereof, and member thereof
JP2008295200A (en) 2007-05-24 2008-12-04 Aisin Seiki Co Ltd Synchronous motor control device and method for optimizing synchronous motor control
JP2013039020A (en) 2011-07-13 2013-02-21 Suri-Ai:Kk Transverse flux machine

Also Published As

Publication number Publication date
JP2021121163A (en) 2021-08-19

Similar Documents

Publication Publication Date Title
JP3328636B2 (en) Vector control method for synchronous reluctance motor
JP3358390B2 (en) Synchronous motor rotation angular velocity detection device, rotation angle detection device, control device and control method for synchronous motor
JP2016171741A (en) Method for detecting initial magnetic-pole position of rotor in permanent magnet synchronous motor
Mevey Sensorless field oriented control of brushless permanent magnet synchronous motors
JP6806983B2 (en) Double three-phase winding permanent magnet synchronous motor drive system
Singh et al. Performance investigation of permanent magnet synchronous motor drive using vector controlled technique
KR20190143630A (en) Sensorless control system and method for permanent magnet synchronous motor
Awaze Four quadrant operation of BLDC motor in MATLAB/SIMULINK
JP7450150B2 (en) A mathematical model of a synchronous motor that generates reluctance torque, and a simulation, characteristic analysis, and control device based on the same model.
Veena et al. Vector control of three phase PMSM drive using power transformations for future spacecraft application
Elsayed et al. Simulation study of a new approach for field weakening control of PMSM
JP6893293B2 (en) Double three-phase winding permanent magnet synchronous motor drive system
JP6971439B2 (en) Dual three-phase winding permanent magnet Simulated / characteristic analysis / control device based on the same model as the mathematical model of synchronous motors
Heidari et al. A comparison of the vector control of synchronous reluctance motor and permanent magnet-assisted synchronous reluctance motor
JP2013039015A (en) Method for controlling driving of permanent-magnet synchronous motor
JP7302140B2 (en) Drive for synchronous reluctance motor
Shinnaka New dynamic mathematical model and new dynamic vector simulators of hybrid-field synchronous motor
JP7190098B2 (en) Phase estimator for synchronous reluctance motor
Ramirez-Leyva et al. Direct Torque Control of a Permanent-Magnet Synchronous Motor with Neural Networks
Mirzaev et al. Mathematical Model of an Asynchronous Motor in Full-Phase Operation
Beccue et al. A coupled piezoelectric/single-hall-sensor position observer for permanent magnet synchronous machines
Ghosh et al. Comparative analysis of sensorless DTC schemes for a synchronous reluctance motor drive
Yang et al. Design and implementation of a FPGA-based AC servo system
JP2018201319A (en) Drive system of double three-phase wiring permanent magnet type synchronous motor
JP4448351B2 (en) Control device for power converter

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221227

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231025

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231107

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231216

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240210

R150 Certificate of patent or registration of utility model

Ref document number: 7450150

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150