TyphoonSim ソルバーエンジン

TyphoonSim ソルバー エンジンのアルゴリズムと機能、およびモデリングの原則とシミュレーション ワークフローの概要。

TyphoonSim のみ:このドキュメントは TyphoonSim シミュレーションにのみ有効です。

TyphoonSimは、ハイブリッド微分方程式で記述される連続および離散ダイナミクスの両方を示す複雑なモデルのシミュレーションを可能にする最先端のソルバーを搭載しています。ソルバー設定の詳細については、「モデル設定」をご覧ください。電気および信号処理領域のモデリングに適用される基本原理については、 「モデリング原理」をご覧ください。

ソルバーアルゴリズム

TyphoonSimソルバーエンジンは現在、 DAEソルバーODEソルバーの2つのソルバーオプションを提供しています。デフォルトオプションであるDAEソルバーは、微分代数方程式(DAE)を形成し、指数削減法を用いて指数1のDAEに変換します。一方、ODEソルバーは、システム方程式を常微分方程式(ODE)に変換します。どちらの場合も、微分方程式はユーザーが指定した可変ステップ積分法によって積分されます。

TyphoonSimは現在、 BDFBDF-RobustTrapezoidalTrapezoidal-Robust の4つの積分法を提供しています。 BDF 法と BDF-Robust 法は、硬い微分方程式を解くために一般的に使用される後退微分公式(BDF)に基づく暗黙的な可変順序可変ステップ法です[ 1 ]。 Trapezoidal 法と Trapezoidal-Robust 法も、標準的な台形積分法に基づく暗黙的な可変ステップ法です。 BDF-Robust や Trapezoidal-Robust などのロバスト積分法は、微分方程式に追加のスケーリング手法を適用し、積分中およびハイブリッドイベントの処理中に、より厳しい収束許容誤差を適用します。通常のオプションが収束しない場合は、ロバストオプションを使用する必要があります。 DAE ソルバーは、大規模な回路のシミュレーションに推奨されます。たとえば、インバーター インターフェイスの分散型発電が広く普及している電力網モデルなど、大きな代数ループを持つ回路は、DAE ソルバーによってより効率的に処理され、大幅に高速に解決されます。

モデリングの原則

修正節点解析

TyphoonSimは、修正節点解析(MNA)を利用して、システム方程式をDAE形式[ 2 ]で構築します。これらの方程式は、回路に存在する可能性のある様々な位相的な矛盾のため、通常は高指数DAEとなります。高指数DAEは数値的に解くのがより困難であるため[ 1 ]、指数1のDAEまたは常微分方程式に変換されます。

電力コンバータ

TyphoonSim には、リアルタイム シミュレーション向けの Typhoon HIL コア ライブラリのプログラム済みおよびパッケージ済みの電力コンバータ コンポーネント ( NPC レッグなど) を使用する代わりに、ダイオードサイリスタIGBTなどの基本的なスイッチング エレメントを使用して電力コンバータの設計とシミュレーションを実行できる高度なディスクリート イベント ハンドラーが搭載されています。ただし、TyphoonSim とリアルタイム/VHIL シミュレーション間のモデルの連続性を確保するために、パッケージ済みのコンポーネントも TyphoonSim を使用してシミュレーションできます。内部的には、それらは個々のスイッチ レベルに分解され、単一のスイッチを使用して構築された場合と同じようにシミュレーションされます。ディスクリート イベント ハンドラーは、最先端の推定手法を使用して、シミュレーション実行中、特にダイオードなどのライン整流スイッチなどのスイッチング デバイスの状態を正確に推定します。

ODE solver models basic switching devices such as IGBTs as ideal switches (Ron = 0, Roff = inf). On the other hand, DAE solver models a switching device as a variable resistor (Ron > 0, Roff < inf), where Ron takes a low resistance value (the default value is 1 mΩ), when a switch is on, and a high resistance value (the default value is 1 MΩ), when a switch is off. Users can also set the turn-on voltage for line-commutated switches, namely, diodes and thyristors.

機械

機械方程式は、明示的な積分法を使用する専用の機械ソルバーによって解かれます。明示的な積分法の実行は、離散イベント ハンドラーによって 1 µs ごとに定期的に、または機械端末に関連付けられた電気信号に十分な変化があったときにトリガーされます。

シミュレーションワークフロー

図 1 は、 TyphoonSim でシミュレーションを実行するために活用される基本的なワークフローを示しています。

1回路シミュレーションフローチャート

回路シミュレーション ワークフローは、初期化フェーズと、以下で説明するいくつかの機能ユニットに分割できるメイン統合ループで構成されています。

  • 状態変数と導関数の初期化:このユニットは、状態変数と微分変数の時間導関数の一貫した初期値を取得します。
  • 状態変数と微分値の計算:可変ステップ積分法は、組み込みの誤差制御スキームを用いて、状態変数と微分変数の時間微分値を計算します。ステップサイズは、状態ベクトルが許容誤差範囲内に収まるように調整されます。
  • 出力の評価:測定コンポーネントに対応する出力信号は、ユーザーが指定したレートで評価されます。
  • 入力の評価:このユニットは、IGBT のゲート信号やマシン回路インターフェースなどの制御可能なソースおよびスイッチング デバイスの入力信号を評価します。
  • 離散イベント検出および処理:このコンポーネントは、システム挙動に不連続性をもたらす離散スイッチングイベントを検出し、処理します。ルート探索アルゴリズムを用いて、ダイオードなどの整流スイッチングデバイスのゼロクロス点を特定し、この情報をモデリングモジュールとソルバーモジュールに提供してモデルを更新し、積分ステップを完了します。また、初期化フェーズなどにおいて、整流スイッチングデバイスの状態を推定します。

参考文献

  1. LR Petzold、 「微分代数方程式の数値解析法:現状と将来の方向性」、 1989 年、 https://www.osti.gov/servlets/purl/5818293。
  2. FN Najm、 「回路シミュレーション」、 John Wiley & Sons、Ltd、2010 年。