回路図エディタでボード線図を使用してコンバータ設計の問題をトラブルシューティングする

Python 固有のライブラリ (scipy など) を使用して、電力コンバータ制御の設計に関する問題を解決する方法のデモンストレーション。

導入

注:このユーザーガイドは、C-HIL:Texas Instruments TMS320F2808カードを使用したPMSMのフィールド指向制御に関するものです。本書に記載されているすべてのテストは、TIカードを使用せずに実施されています。
注:このアプリケーション ノートで説明されているモデルは、Typhoon HIL のパッケージ マネージャー サービスにあります。このサービスには、Typhoon HIL コントロール センター ツールチェーン内のツールとして、またはパッケージ マネージャー Web サイトからアクセスできます。

電気駆動アプリケーションにおいて、整流器は直流電圧オーバーシュートの問題に悩まされる傾向があります。本稿では、図1に示すような制御されていない三相ダイオード整流器を用いて、これらの問題の発生可能性を検証します。入力インピーダンスはボード線図を用いて計算し、モデルに導入することができます。この直流リンク電圧のピークは、三相ダイオード整流器のダイオードや、オーバーシュート電圧の定格を満たしていない三相インバータのIGBTの焼損につながる可能性があります。

モデルの説明

駆動システムの整流器部分のモデルを図1に示します。

1駆動システムの整流器部分
2デフォルトパラメータでのオーバーシュート

以下のコードブロックのモデル初期化関数を使用することで、Pythonライブラリのscipymatplotlibを用いてボード線図を生成できます。scipyライブラリの最も重要な関数は、 scipy.signal.TransferFunctionscipy.signal.bodeです

モデル初期化スクリプトを以下に示します。

# Numpy module is imported as 'np'
# Scipy module is imported as 'sp'
import matplotlib.pyplot as plt
          
#parameters that chan be changeable
Rprech = 0
Ts = 100e-6
R = 0.1
L = 0.0011
C = 0.0022
RloadAnalys = 100
          
#Solution no. 1: fix with inductor
#L = C * R**2 *0.9**2
#mdl.info(L)
         
#Solution no. 2: fix with capacitor
#fsw = 100
#C = 1/(2*np.pi*0.3*fsw)**2/L
#mdl.info("C = {} F".format(C))
          
#Solution no. 3: add precharge resistor Q = 0.707
Rprech =np.sqrt(2)*(np.sqrt(L/C)-1/np.sqrt(2)*R)
mdl.info("Rprech = {} ohm".format(Rprech))
          
#bode diagram analysis
# resonant frequecny
f0 = 1/(2*np.pi*np.sqrt(L*C)) 
mdl.info("f0 = {} Hz".format(f0))
w0 = 2*np.pi*f0
          
# absolute Q factor
Qab = 1/(R+Rprech)*np.sqrt(L/C)
mdl.info("Qab = {}".format(Qab)) 
Qdb = 20 * np.log10(Qab)
# Q factor in db         
          
#bode plot using sp.signal.lti function
s1 = sp.signal.lti([1],[1/w0**2,1/(w0*Qab),1])
w, mag, phase = sp.signal.bode(s1)
         
# Bode magnitude plot
#plt.figure()
plt.semilogx(1/(2*np.pi)*w, mag)
# Bode phase plot
plt.figure(1,1)
plt.semilogx(1/(2*np.pi)*w, phase)
plt.show()        
3入力伝達パラメータのボード線図

このボード線図と計算から、 f 0 = 1 2 π L 1つの C d c = 102.3   H z スイッチング共振周波数とダイオードブリッジのスイッチング周波数(100Hz)がほぼ同じであることがわかります。これは、DC電圧に共振とオーバーシュートが発生することを意味します。ダイオード整流器を使用しているため、スイッチング周波数も変更できません。代わりに、入力インピーダンスのパラメータを変更する必要があります。モデル抵抗(R)を増やすと損失が増加するため、インダクタンス(Q)を下げてインダクタンスに着目する方が適切です。

シミュレーション

このセクションでは、DC リンクの電圧オーバーシュートを解決するための 2 つの解決策として、インダクタンスの低減と静電容量の増加について説明します。

インダクタンスを減少させる問題の解決

まず、入力伝達関数を求めてみましょう。伝達関数は図1の回路から定義されます。

G s =   1 1 + s R 1 C d c + s 2 L 1つの C d c

伝達関数の代替標準正規化形式は次のとおりです。

G ( s ) = 1 1 + s 質問 ω 0 + s ω 0 2

パラメータQは回路の品質係数と呼ばれ、システムにおける消費電力の尺度です。このパラメータは共振角周波数を表します。

入力伝達関数から、共振周波数は次のようになります。

f 0 = 1 2 π L 1つの C d c

次の式で Q 係数を簡単に見つけることができます。

質問 = 1 R 1 L 1つの C d c

この解析結果から、Q値を下げるためのインダクタの値が分かります。インダクタの値を半分に減らしてみましょう。図4を見ると、DC電圧のオーバーシュートが低減していることがわかります。

4インダクタンスが半分に減少した場合のオーバーシュート

さて、インダクタンスを3分の1に下げてみましょう。図5を見ると、DC電圧(Cdc)を目標の400V以下に下げることができていることがわかります。

5インダクタンスが3倍に減少した場合のオーバーシュート

二次的な バターワースフィルタ (つまり、通過帯域周波数応答が最も平坦な連続時間フィルタ)の減衰不足時のQ値は0.707です。したがって、DCコンデンサの電圧のオーバーシュートを低減するには、 Q係数を0.707に設定するさて、インダクタを計算してみましょう。式は次のようになります。 L 1つの c = C d c R 1 2 質問 2 計算できる場所 L 1つの c =   1.782 e - 05 .

図 6では、この問題を解決できたことがわかります。

6計算されたインダクタンスで修正された問題

図 7のボード線図では、Q 係数を下げることができることがわかります。

7 .インダクタンスが減少した入力伝達関数のボード線図

コンデンサを増やすことで問題を解決する

理論的にはインダクタンスを下げれば問題は解決しますが、インダクタはコンデンサよりもはるかに高価なため、現実的ではない場合が多くあります。さらに、モデル内のインダクタンスが電源インダクタンス(つまりグリッドインピーダンス)に起因する場合、それを変更することは不可能です。このような場合、代わりにコンデンサの値を変更することができます。これは、共振周波数をスイッチング周波数から遠ざけることで実現できます。例えば、次のように設定します。 f 0 = 0.3 f s :

C d c = 1 2 π 0.3 f s 2 L 1つの =   0.0255   F

この新しいコンデンサにより、スイッチング周波数の減衰が得られます。図8に新しいボード線図を示します。

8コンデンサの新しい値を持つボード線図

図9では、DCリンクの電圧の新しい応答が確認できます。DCコンデンサを増やすと、インダクタンスを変更したときよりも応答が少し遅くなることがわかります。

9 DCリンクの電圧応答

プリチャージ回路

DCリンクの過電圧は、プリチャージ回路によって解決できます。図10に、プリチャージ回路の実装を示します。

10 DCリンクプリチャージ回路

この場合、Q係数を設定することでプリチャージ抵抗を計算できます。

質問 = 1 2

以下は計算式です R p r e c h :

質問 = 1 R 1 + R p r e c h L 1つの C d c = 1 2

R p r e c h =   2L 1つの C d c - R 1

図11は、Q値の低下によって生じるプリチャージ抵抗の影響を示しています。図12は、 DCリンクの電圧応答を示しています。この場合、整流器のAC側への入力電流のピークが減少していることに注目することが重要です。

11プリチャージ回路を備えたボード線図
12プリチャージ回路を用いた電圧応答

このアプローチを用いることで、パワーエレクトロニクスにおける様々なコンバータやアプリケーションの設計における問題のトラブルシューティングが可能になり、特にDCリンクにおける電圧オーバーシュートを回避できます。図12では、インダクタ電流のピーク値が高いことがわかります。この電流を低減するには、プリチャージ抵抗の値を変更する必要があります。プリチャージ回路の抵抗値は、負荷容量と必要なプリチャージ時間に基づいて選択します。プリチャージサージ電流は、以下の時間後に初期値の1/eに達します。

T   =   R p r e c h     C

約 5 * T の時間が経過すると、電流は管理可能な値まで減少します。

R p r e c h = T 5 C

13新しい抵抗器を使用したプリチャージ回路のテスト。

ここで、突入電流をどのように低減できるかがわかります。

14新しい抵抗器を使用したボード線図

実際には、設計が許せば、より大きなコンデンサを配置することで過電圧を解消する選択肢を持つ方が、プリチャージ回路を追加するよりも安価になる可能性があります。これらの結果に基づいて、これを実現する方法を想像することができます。

テスト自動化

この例のテスト自動化はまだありません。ご協力いただける場合はお知らせください。アプリケーションノートへの署名を喜んで承ります。

要件の例

表1は、モデルをリアルタイムで実行するためのファイルの場所とハードウェア要件に関する詳細情報と、この最小限のハードウェア構成でモデルを実行した場合のHILデバイスのリソース使用率を示しています。この情報は、モデルの実行とカスタマイズを必要に応じて行う際に役立ちます。

1 .最小要件
ファイル
Typhoon HILファイル

TMS320F2808 サンプルモデル/ti pmsm センサー付きフォーカス*

ti pmsmセンサー付きfoc.tse

ti pmsmセンサー付きフォーカス

20140220_pmsm3_1_TMS3202808.out

設定.runx

*パッケージマネージャー経由でダウンロード可能

外部ツール ユニフラッシュ
最小ハードウェア要件
HILデバイス数 1
HILデバイスモデル HIL402
デバイス構成 1
インタフェース

TMS320F2808 (テスト対象ハードウェア)

HIL DSP180 インターフェース

HILデバイスのリソース利用
処理コア数 2
最大マトリックスメモリ使用率

2.25% (コア0)

30.81% (コア1)

最大時間枠利用率

29.38% (コア0)

52.5%(コア1)

シミュレーションステップ、電気 1マイクロ秒
実行率、信号処理 100マイクロ秒

著者

[1] シミサ・シミッチ