非線形なシステムへの対応

$$F_{t r}=F_{r r}+K_{a e r o} V^{2}+M \frac{d V}{d t}$$

空気抵抗は車速の二乗に比例する非線形な特性です。
次にこの車両の駆動力と走行抵抗の物理式から、$V$と$V^{2}$を状態変数に選び、状態方程式を生成します。

$$\frac{d V}{d t}=-\frac{K_{\text {aero }}}{M} V^{2}+\frac{1}{M}\left(F_{t r}-F_{r r}\right)$$

$$\frac{d V^{2}}{d t}=\frac{d V^{2}}{d V} \frac{d V}{d t}=2 V \frac{d V}{d t}=-\frac{2 V K_{\text {aero }}}{M} V^{2}+\frac{2 V}{M}\left(F_{t r}-F_{r r}\right)$$

状態方程式にまとめると

$$\frac{d}{d t}\left[\begin{array}{l}V \\ V^{2}\end{array}\right]=\left[\begin{array}{cc}0 & -\frac{K_{\text {aero }}}{M} \\ 0 & -\frac{2 V K_{\text {aero }}}{M}\end{array}\right]\left[\begin{array}{l}V \\ V^{2}\end{array}\right]+\left[\begin{array}{c}\frac{1}{M} \\ \frac{2 V}{M}\end{array}\right]\left(F_{t r}-F_{r r}\right)$$

平衡動作点の近傍での作動であれば、平衡動作点で非線形性を線形近似して状態方程式に取り込めばよいです。

例えば、速度30[m/s]で走行すると、駆動力750[N]、転がり抵抗 𝐹𝑟𝑟= 300[N] 、空気抵抗 𝐹𝑎𝑒𝑟𝑜= 450[N] で平衡状態になります。速度30[m/s]の平衡動作点での挙動を確認したい場合は、速度30[m/s]での状態方程式のAとBを離散化し、得られた逐次式のPとQで時間応答を確認できます。

Excelファイルのダウンロード

非線形モデルのEXCELファイルはこちらからダウンロードできます。

EXCELの状態方程式のシートでは、状態方程式を離散時間モデルの逐次式に変換しています。
シミュレーション (動作平衡点で近似)のシートは、特定の平衡動作点で線形近似した場合のシミュレーション結果を示しています。
シミュレーション (離散時間ごとにPQ更新)のシートは、離散化したすべての動作点でPとQを再計算して線形近似した場合のシミュレーション結果を示しています。

誤差グラフのシートでは、上記の計算方法での空気抵抗の計算結果の精度比較を示しています。

平衡動作点の近傍で時間応答を算出

平衡動作点が特定できないケースの時間応答算出方法

平衡動作点が特定できない場合は、動作点ごとに非線形特性を線形近似した値を状態方程式のAとBのパラメータに設定する必要があります。時間応答の確認のためには、AとBのパラメータの更新に合わせ、離散化した逐次式のPとQのパラメータも動作点ごとに更新する必要があります。

この場合、時間応答の計算は、EXCELシート上での計算だけでは無理で、ちょっとVBAも使って離散時間（動作点）ごとにPとQを計算しています。

シミュレーション（離散時間ごとにPQ更新）のシートは、速度30[m/s]で走行している状態で駆動力$F_{tr}$を750[N]から1200[N]に上げた場合について、PとQのパラメータを動作点ごとに更新した場合のシミュレーション結果を示しています。

システムの非線形性に対しては、どのように対応するかはケースバイケースで考える必要があります。

VBAプログラムの説明

VBAプログラムでPとQのパラメータを動作点ごとに更新しています。
※EXCELの「開発」タグを選択し、「Visual Basic」からプログラムを確認できます。

シミュレーション結果

誤差グラフのシートでは、特定の平衡動作点で線形近似した場合と、動作点ごとに線形近似した場合について、算出される空気抵抗の誤差を比較しています。

※EXCELシート上部のスライドボリュームで駆動力 Ftr を変更できますが、変更した場合は、「非線形モデル時間応答」のボタンを押して再計算する必要があります。

動作点ごとに線形近似した場合は、動作点が動いても誤差は小さいことが分かります。
車両の駆動力と走行抵抗のモデルでは、平衡動作点が特定できないため、動作点ごとに離散化した逐次式のP,Qのパラメータを算出した方がより適切なシミュレーション結果になります。

非線形モデルをどのように状態方程式に取り込むかはシステムの特性や作動条件に応じて検討する必要があります。

最後までお読みいただき、ありがとうございました。

次は現代制御のオブザーバ（状態観測器）になります。