電圧形PWMインバータ駆動永久磁石同期電動機の定
常解析
著者
篠原 勝次, 山本 吉朗, 実成 義孝, 入佐 俊幸
雑誌名
鹿児島大学工学部研究報告
巻
30
ページ
53-59
別言語のタイトル
Steady state analysis of an interior
permanent-magnet synchronous motor driven by a
PWM voltage source inverter
電圧形PWMインバータ駆動永久磁石同期電動機の定
常解析
著者
篠原 勝次, 山本 吉朗, 実成 義孝, 入佐 俊幸
雑誌名
鹿児島大学工学部研究報告
巻
30
ページ
53-59
別言語のタイトル
Steady state analysis of an interior
permanent-magnet synchronous motor driven by a
PWM voltage source inverter
電圧形PWMインバータ駆動永久磁石同期電動機の定常解析
篠原勝次・山本吉朗・実成義孝・入佐俊幸
(受理昭和63年5月31日) STEADYSTATEANALYSISOFANINTERIORPERMANENT−MAGNET SYNCHRONOUSMOTORDRIVENBYAPWMVOLTAGESOURCEINVERTER KatsujiSHINOHARA,KichiroYAMAMOTO,YoshitakaMINARI andToshiyukilRISA InteriorPermanent-Magnetmotorsfedfrominvertersarebecomingincreasinglyattractiveinawide varietyofspeedcontrolapplicationsoflowandmediumpowerrange、I
n
t
h
i
s
p
a
p
e
r
,
s
t
e
a
d
y
s
t
a
t
e
a
n
a
l
y
s
i
s
i
s
c
o
m
p
u
t
e
d
b
y
t
h
e
s
t
a
t
e
-
v
a
r
i
a
b
l
e
m
e
t
h
o
d
・
C
o
n
s
e
q
u
e
n
t
l
y
,
t
h
e
followingresultswerefound:(
1
)
A
l
t
e
r
n
a
t
i
n
g
-
c
u
r
r
e
n
t
c
o
m
p
o
n
e
n
t
s
o
f
d
a
n
d
q
s
t
a
t
o
r
a
x
e
s
c
o
r
r
e
s
p
o
n
d
t
o
a
l
t
e
r
n
a
t
i
n
g
c
o
m
p
o
n
e
n
t
s
o
f
d
andqstatorvoltageswhicharechangedbythePWMpattern.(2)Directionsofdandqdumpercurrentsareoppositethoseofdandqstatorcurrents.
1.まえがき 近年,半導体素子の発達に伴いAC可変速システ ムの高性能化が進んでおり,誘導電動機,同期電動機 共にベクトル制御を始めとする様々な高速応答制御が行われている(1)(2)。またトルクリプルや騒音の低減の
為にPWM方式を用い,より高速なスイッチングに よりその電流波形を正弦波に近づけようとする試みが なされている。このような状況にあって永久磁石同期 電動機は可変速駆動に用いた場合,誘導機と比べ一般 にかなり効率がよく,磁性材料の高性能化とあいまっ て今後その応用は一層拡大されると思われる。 筆者らは先にPWMインバータを用いて永久磁石 同期電動機のベクトル制御システムを試作し,その実 験結果とシミュレーション結果についての比較を行った(3)。解析において,電流はPWM時のリプルを含
んだ波形ではなくその基本波成分についてのみ考慮し た。今回はこの実験回路における同期電動機の定常時 の電流波形についてより厳密にその解析を行った。本 稿では,まずこの実験回路の回路動作の説明を行った 後,解析に用いた状態方程式を導出し,さらに数値計 算結果について検討を行う。 2.回路動作 図1に実験回路の構成図を示す。このシステムは, (1)式で示されるように速度指令値‘Uγ*とインクリメ ンタルエンコーダ(600パルス/回転)によって検出 される実速度①γとの偏差をP増幅してトルク指令 T*とし,さらに(2)式により電流指令j*を決定して いる。 T*=KP(のγ*−Cuγ)………・……・…・……(1)‘拳=器………(2)
ただし,KP:比例ゲイン Kr:トルク定数 このj*に回転子位置に同期した3相正弦波を乗じて 各相の電流指令ju*,jひ*,j"*を得る。正弦波はアブ ソリュートエンコーダ(8ビット)を用い,回転子の 磁極位置と対応したEP−ROMのアドレスより正弦 波データを読み出し,これをDA変換して作られる。 ju*,jひ*,j"*は回転子速度と同じ周波数を持ち,回 転子上の。軸と常に一定の位相関係を保った信号と なるため,必要とされるトルクに応じてj*つまり相 電流指令の振幅のみを制御すればよい。次段の電流制 御回路では,電流指令ju*,if,*,Z"*と絶縁アンプを表 l 同 期 電 動 機 の 定 数 54 鹿 児 島 大 学 工 学 部 研 究 報 告 第 3 0 号 ( 1 9 8 8 ) 3 . 2 P W M パ タ ー ン 及 び 電 圧 モ ー ド 図4に3相の変調波と三角波(搬送波)によって作
"
|
」
」
壷
」
割
1000 邸F '3‘」 」
A、R・E 8bit,」
絶縁アンプ 各部ゲートへ 0.01 100k.0皿F 100k鋸
12k ★吟卦 崎10k毛
虻
1
1
・10k,4.7k FETゲート回路 三相正弦波発生器 ★v 01 %変換器 0.1倍 '0sin(。‘一号燕l『× ★四 .8 lOeユ、(α‘一寺元) × 角 波 発 振 器 図 1 実 験 回 路 介して検出してきた実際の相電流との偏差をPI増幅 し,これと三角波との比較を行うことによりPWM 波形を得ている。このシステムに用いた供試機の定数 については,表1に示す。 られるPWM波形を示す(3=5.2。)。ひ狸,りひ,ひ”は 同期電動機に印加される電圧であり,りsd,ひsoはこ れら3相を図3に基いてd−q変換した2相の電圧 である。図4には電気角で60.の区間が示されている が,この次の60.の区間ではりu=−ひび,ひじ=−ひ", ひ"=一ひ泌でありりsd,ひsoは図4に示されている区間のパターンと全く同じものになる(5)。このため60.以
降は省略してある。 図5は図4の破線部を拡大したものである。E’は インバータ入力直流電圧である。図において△ひsα(t"), △ひs・(t,jはそれぞれ時刻tnでのモード変化による Usd,ひs9の瞬時の変化分を示している。 ひsd,ひs・はそのスイッチングのモードにより図6に示すような電圧ベクトルのモードをとる(6)。例えば
V6(1,0,1)はインバータのu相の上側アームのFET 電 機 子 d 軸 イ ン ダ ク タ ン ス L d ダンパ巻線インダクタンスLdr,L9『 5.3mH 相 数 、 3 16.1A 極 数 P 4 永 久 磁 石 等 価 界 磁 電 流 ル 電 機 子 抵 抗 凡 0 . 5 1 Q 電機子q軸インダクタンスL9 2.7mH ダンパ巻線抵抗Rdr,Rq『 0.2Q 電機子漏れインダクタンスLI 3.定常特性の解析 3.1永久磁石同期電動機の等価回路 回転座標系において永久磁石同期電動機の等価回路は図2のように書ける(4)。ここで用いた回転座標系
(d−q軸座標系)は,図3に示されるようにq軸が 。軸よりも90.進んでおり,反時計方向に同期速度で 回転している座標系である。d軸とu相軸との角度 をβとすると,8は同期速度⑳e,初期位相aOを用 いて次のように表せる。 β=Cuet+ao・…・……・………..(3) 図2の等価回路では,永久磁石形であるので界磁回 路を等価な定電流源乃と置き換えてある。永久磁石 の生ずる磁束蝿αgは,この〃を用いて偽αg=L,mdL と表すことができる。また使用した電動機は,かご形 巻線を有しておりこれがダンパ巻線として作用する。 等価回路ではこのダンパ巻線を等価な。,qの巻線 に分解し,それぞれ。軸,q軸の2次側の短絡され た閉回路としてこれを表している。 44.3mH 18.8mHがON,ひ相の下側アームのFETがON,”相の上
側
ア
ー
ム
の
F
E
T
が
O
N
,
つ
ま
り
'
趣
=
等
'
,
=
一
等
‘
"
=
号
の
状
態
を
示
し
て
い
る
。
二
の
電
気
角
0
.
-
6
0
.
の
期間においてはインバータはVケ,V6,V5,Vi,Vb のモードをとり,このときのりsd,ひsqは次のように なる。 V7,Vbモードのとき{
:
:
:
二
:
…
…
…
…
…
一
…
…
側
V6モードのとき 変調波 e“ 変 調 波 eu q軸篠原・山本・実成・入佐:電圧形PWMインバータ駆動永久磁石同期電動機の定常解析55
M】 リ ︵︶ d 軸 町 L‘q Rs Zsq Lkq D 0 U = C 相軸 U u U U U U U U U 〕00 RA、9○
入Sd=(LId+L耐α)jsd+Lmdjkd+LmdL 入s・=(L1.+Lm。)jsq+Lmqjkq 回転座標系における永久磁石同期電動機の等価 回路 切 相 軸 図3u,ひ,切相軸とd−q軸との関係 図2肝
搬 送 波F
R
F
H
1 1上
F
上
F
H
帝 か君 変 調 波 e,, 電気角60。h
Uul I I I I I I I I I n l
■■■■■■■■■■■■■■I
l
l
U
l
I
l
'
T
,
r
i
i
T
i
i
f
l
l
二
'
二
'
二
'
二
l
R
R
F
H
二
HHH
十
◎ハⅢU。︿叩﹀◎ハ叩U
1H捌州1H
J〃しJrTしF1〃しP
処仙畦仙畦仙
皿娠皿に皿応
ScSCSC
、″、〃n“、〃、″、〃
EEEEEE
〃﹂W一吻唱W略一Ⅲ﹄W−I↑W蛸−4W一Ⅲ↑W一一一一の一一一一の一一一一
皿”吟瞬︲湿珂吟ト︲誕諏
研皿脹︲︲︲翌、︲︲︲︲m︾一毎皿脹︲︲︲組︲︲︲︲︲、﹀一mm眠︲︲︲忽︲︲︲︲︲、︶ モ モ 腿 wW, Uu1’
U U U ・・・…・…………(5) リ銅 ……・…..…・…(6)岬Ⅲ11冊lWlm
Lg−−j t64 Z O 図4パルス幅制御電圧のd−q変換(6=5.2。) .。…………..…(7)ー − ー 56 5モード ひsd,ひs9は上式からわかるようにVケ,Vbのモード 以外は,どのモードでも次のような形をとる。
│
:
魚
鳳
これよりVケ,Vbモード以外については微分演算子を pとすると以下のような関係が成立する。{
:
:
:
:
二
雪
:
:
;
茎
…
…
…
…
…
…
…
…
…
…
…
⑧
変調波が10(Hz),三角波が630(Hz)のときには りsd,ひs9のモード変化は電気角60°の区間で63回存 在する。ここで, x=〔ひsd,ひsq,L,is。,jsq,jltd,jk。〕T、.(9) 次式で△工”を定義する。 △工"=〔Usd(t、+)−Usd(tm-),ひso(t"J−ひso(in_),0,0,0,0,0〕T ………・…(10) 、=1,2,3,……,63 リ,
/
吾
E
,
ド Usd 八 0 ド〃″
EE
信杉FFww−の峰一W”程式は次のように表せる(7)。
;1仰坤'ふ
譜則り+卜汁⑪
l
l
H
J
灘
川
嶋
轍
モ ー ド ひ6 -一一一一一一r △ひsq(t2)! 鹿 児 島 大 学 工 学 部 研 究 報 告 第 3 0 号 ( 1 9 8 8 ) l o t , 0T
−
△Usq(j3) Usq △Usq(t,)ll
L
‐
RjtdLmd/△d※
一
①
e
L
m
d
(
L
I
d
+
L
m
q
)
/
△
q
−Rl上d(Lld+Lmd)/△d cUeLmdL瓶q/△q 図5Usd,Us9の変化美
穂
卜
図 6 電 圧 ベ ク ト ル i 2 3i(13),(14)式はそれぞれ次の状態方程式を表してい る。
p工=A1工(vi,vb,v6モードのとき)…(15)
p工=A2工(M,Vbモードのとき)………(16) O O O O O O ※−Rs(Lkd+L耐d)/△.⑳e(L臆d+L、。)(L1q+L耐9)/△d※ −cue(Lld+L耐。)(Lk9+L"9)/△q−Rs(L厳9+L"9)/△q RsL耐d/△.−CUeL耐d(LI9+Lm9)/△d CueL洞9(Lld+L凧d)/△9RsL耐9/△q{
i
化
製
訓
⑫
蕊
釜
訓
汁
※ 3.4初期値の決定 (15)式および(16)式により状態推移行列を 。A,(△t")=explA,(△Zn)| oA2Wn)=eXPlA2(△t、)| と定義する。ただし△t"=tn-tn-,(、=1,2,…… 64) この状態推移行列を用いると to≦t<t,のとき工
(
t
)
=
①
A
2
(
t
)
工
(
0
+
)
=
①
A
2
工
(
0
-
)
〔
.
.
・
工
(
0
J
=
工
(
0
-
)
〕
t=t,のとき 工(t,)=△工,+①A2(△t,)工(0-) t,≦t<t2のとき 工(t)=。A,(t-t,)〔△工,+のA2(△t,)工(0−)〕 t=t2のとき工
(
t
2
)
=
△
工
2
+
の
A
l
(
△
t
2
)
〔
△
工
,
+
の
A
2
(
△
t
,
)
工
(
0
-
)
〕
t=t64-のとき工
(
t
6
4
_
)
=
の
A
2
(
△
t
6
‘
)
〔
△
工
6
3
+
①
A
l
(
△
t
6
3
)
〔
…
…
……〔△工,+①A2(△オ!)工(0-)〕〕……〕 t二二t64‐のときの式を変形する。 工(t64-)=①A2(△t64)△工63 +①A2(△t64)①A1(△t63)△工62 +①A2(△t64)のA1(△t63)……のAl(△t2)△工1+
①
A
2
(
△
t
6
4
)
①
A
1
(
△
t
‘
3
)
…
…
の
A
l
(
△
t
l
)
△
工
(
0
-
)
…
(
1
9
)
1
0 0 0 0 0 0 0 0 O Lhd+L凧d)/△d O O ※ 0(Lk9+Lmq)/△9−⑳eL,、d(L脆・+L凧9)/△・ 一 L 耐 。 / △ d O O O一Lmq/△q“eLmdL,、o/△・ (13) 0 〔 U e O O − “ e O O O O L k d + L 碗 d ) / △ d O O ※ 0(Lk9+L凧。)/△。−“eL耐d(Lj『。+L耐9)/△o −Lmd/△d O O O − L 耐 。 / △ ・ ⑳ e L m d L 洞 q / △ 9 恥恥ア々・粒.恥恥・恥 …(14) p 篠原・山本・実成・入佐:電圧形PWMインバータ駆動永久磁石同期電動機の定常解析57 恥恥r〃・地・軸・坤・恥 p O O O O O O ※RkdLmd/△.①eL耐9(Lkd+L耐d)/△d -cueL耐d(LIE9+L醜9)/△9RIwL碗9/△q −Rltd(L‘d+L"d)/△.−‘ueL,ndL嗣q/△。 〔ueL耐dL醜9/△9−R廊9(L‘’十L耐q)/△q また恥,Vbモードに対して5 58 10 この連立方程式(19)をガウスの消去法により解き 工(0-)を求める。 1.0 15 5 0 ︹く︺でせ﹃
W
I
‘
了
A
汽
十
A
I
1
W
1
3.5数値計算結果 上で求めた初期値を用い,(15),(16)式の連立微分 方程式をルンゲクッタ法により解き,数値解析を行っ た。この数値解析の結果を図7∼9に示す。これらは, 同期速度2Mrad/s)(電気角),インバータ入力直 流電圧E、=70(V)の場合の解析結果である。図7は 回転座標系における固定子電流jsd,jsqの波形であ る。図8は,ダンパ巻線電流j極,恥の波形であり, 図9は回転子座標系上にあるjsd,js9に回転→静止 変換を行い,2相3相変換によって求めたu相電流 の波形である。図4においてりsd,ひs9は3相電圧の モードパターンによって激しく変化するが,この変化 に対する各電流の変化の様子が厳密にシミュレートさ れている。ダンパ電流は,正弦波駆動の場合であれば定常運転時には零となるがPWM駆動の場合,電圧
の瞬時変化に対して過渡的に逆向きの電流が流れてい る。即ちりsd,Us9の変化による空隙磁束の変化を抑 制する方向に流れる。 0 −0.5 −1.0 1.0J
−1i
.0I
I
I
I
W
州
I
図8ダンパ巻線電流波形 8.0 50 1V 戸 、 < ー ノ 1 ; ; 4J0 10( 〔msec〕 6.0 0 鹿児島大学工学部研究報告第30号(1988) 15 図 7 固 定 子 電 流 波 形 〔mBec〕 3.0 < 〔mBec〕 01
W
‘
,
、
寂
0 10 8.0 5 −3.0 三 一 : 『 ・鳥凸 4.0 −6.0 図 9 u 相 電 流 波 形篠原・山本・実成・入佐:電圧形PWMインバータ駆動永久磁石同期電動機の定常解析59 4 . 結 論 本稿では,回転座標系における等価回路から,永久 磁石同期電動機の基本方程式,状態方程式を導出し, PWMパターンを考慮に入れて,回転座標系でその定 常特性の解析を行った。その結果,次の事が明らかに なった。 (1)固定子の。.q軸電流の交流分はPWM変調 された。.q軸電圧の交流分に応答して流れている。 (2)回転子の。.q軸ダンパ電流の変化は固定子 。.q軸電流の変化と逆方向である。 5 . 文 献 l)ThomasMJahns:“Flux-WeakeningRegimeOp‐ erationofanlnteriorPermanent-MagnetSynchro‐ nousMotorDrive”,IEEETm〃s・加。〃stγ・Ap. ”‘.,肌一〃,681(1987) 2)S・Ogasawara,M・Nishimura,H・Akagi,A、 Nabae,Y・Nakanishi:“AHighPerformanceAC ServoSystemwithPermanentMagnetSynchro‐ nousMotors'’’1EEETm”s・〃。〃stγ・Ejectγりれjcs 伽加J肋s”畑.,IE 33,87(1986) 3)篠原・山本・実成・入佐:「永久磁石同期電動機 のベクトル制御とそのシミュレーション」昭63電 気学会全国大会No.1511 4)ThomasMJahns,GeraldB・Klimann,Thomas W・Neumann:“InteriorPermanent-MagnetSyn・ chronousMotorforAdjustable-Speeddrives,,, IEEETm"s,伽"stγ、A妙"c、,Z4-22,681(1986) 5)篠原・実成・入佐・崎山:「誘導電動機のベクト ル制御とそのシミュレーション」,電気学会研究 会半導体電力変換研究会資料SPC−88−24(昭 63) 0)半導体電力変換回路,P137,電気学会 7)電気機器学,P141,電気学会