運動方程式 (2DOF)
$\dot{\beta}= -r + \frac{C_f+C_r}{mV}\beta + \frac{l_f C_f - l_r C_r}{mV^2}\delta$
$\dot{r}= \frac{l_f C_f - l_r C_r}{I_z}\beta - \frac{l_f^2 C_f + l_r^2 C_r}{I_z V}r + \frac{l_f C_f}{I_z}\delta$
線形自転車モデル(2DOF)でヨーレートと横すべり角の過渡応答を計算。車両パラメータを変えてアンダーステア・オーバーステアの特性をリアルタイムに体験しよう。
$\dot{\beta}= -r + \frac{C_f+C_r}{mV}\beta + \frac{l_f C_f - l_r C_r}{mV^2}\delta$
$\dot{r}= \frac{l_f C_f - l_r C_r}{I_z}\beta - \frac{l_f^2 C_f + l_r^2 C_r}{I_z V}r + \frac{l_f C_f}{I_z}\delta$
このモデルの核心は、横すべり角 $\beta$ とヨーレート $r$ の時間変化を記述する2つの連立微分方程式です。車両の横方向の力のつり合いと、ヨー方向のモーメントのつり合いから導出されます。
$$ \begin{aligned}\dot{\beta}&= -r + \frac{C_f + C_r}{m V}\beta + \frac{l_f C_f - l_r C_r}{m V^2}\delta \\ \dot{r}&= \frac{l_f C_f - l_r C_r}{I_z}\beta - \frac{l_f^2 C_f + l_r^2 C_r}{I_z V}r + \frac{l_f C_f}{I_z}\delta \end{aligned}$$$\beta$: 車体横すべり角 [rad]、 $r$: ヨーレート [rad/s]、 $m$: 車両質量 [kg]、 $I_z$: ヨー慣性モーメント [kg·m²]、 $V$: 車速 [m/s]、 $\delta$: 前輪舵角 [rad]、 $C_f, C_r$: 前後コーナリングパワー [N/rad]、 $l_f, l_r$: 前後軸から重心までの距離 [m]
操縦安定性の指標となる「アンダーステア勾配 $K$」は、以下の式で定義されます。この値の正負でアンダーステア(US)かオーバーステア(OS)かが決まります。
$$ K = \frac{m}{L}\left( \frac{l_r}{C_f}- \frac{l_f}{C_r} \right) \quad (L = l_f + l_r) $$$K > 0$ のときアンダーステア、$K < 0$ のときオーバーステアです。市販車は安全のため $K > 0$ (アンダーステア)にチューニングされることがほとんどです。この式から、前輪のグリップ($C_f$)が弱まると $K$ が大きくなり、アンダーステア傾向が強まることがわかります。
車両開発・チューニング:新車の基本設計段階で、目標とする操縦安定性(アンダーステア量)を達成するために、サスペンション特性や重量配分を決定する際の基礎データとして広く用いられます。シミュレータでパラメータを変える作業は、開発現場そのものです。
ドライビングシミュレータ:ゲームやプロドライバー訓練用のシミュレータの物理エンジンでは、この線形モデルを発展させた非線形モデルが使われています。リアルな挙動再現の出発点となる重要なモデルです。
安全性解析:急ハンドルや緊急回避時の車両の挙動を予測し、ESC(横滑り防止装置)の制御ロジックを設計するための基礎理論として応用されています。横すべり角の応答は特に重要です。
学生教育・研究:自動車工学や機械力学の授業で、操縦安定性の概念を理解するための最も標準的な教材です。数式だけではなく、このようなシミュレータで直感的に学ぶことができます。
このシミュレーターを使い始めるとき、いくつか気をつけてほしいポイントがあるよ。まず、これは線形モデルだってことを忘れないで。つまり、タイヤのグリップ力(コーナリングパワー)は横すべり角に比例すると仮定しているんだ。でも実際のタイヤは、ある程度以上すべるとグリップ力が頭打ちになったり減ったりする(非線形性)。だから、このツールで大きな舵角を入力すると、計算上はヨーレートが際限なく大きくなることがあるけど、それは物理的にはありえないんだ。実務では、このモデルは小さいステア入力での基本的な挙動傾向を評価するための第一歩と心得よう。
次に、パラメータ設定の落とし穴。例えば「前コーナリングパワー Cf」の値をいじると、アンダーステア傾向が変わるのはわかるよね。でも、実車ではタイヤの空気圧を変えたり、サスペンションのキャンバー角を変えたりすることで、この値が変化するんだ。シミュレーション上で「Cfを小さくする」のは、「前輪のグリップを意図的に悪くしている」のと同じことだとイメージすると、実車チューニングとの結びつきが理解しやすいよ。
最後に、速度Vの重要性。このモデルは一定速度での応答を計算している。加速や減速を伴う実際のコーナリングとは状況が違う。特に「特性速度」の概念は、あくまでこの線形定常円旋回モデルから導かれる理論値だ。実際のスポーツカー開発では、この考え方を発展させて、様々な速度域でのバランスを総合的に評価するんだ。
この「線形自転車モデル」で学ぶ力学的な考え方は、CAEの世界のあちこちで顔を出すんだ。まず真っ先に挙がるのはサスペンション設計だね。モデルに出てくるコーナリングパワーは、タイヤ単体の値だけど、サスペンションのキングピン傾角やキャスター角などが、接地部の特性にどう影響するかを考える上での基礎になる。
もっと直接的に発展するのが、ESC(横滑り防止装置)やEPS(電動パワーステアリング)の制御ロジック開発だ。これらの先進運転支援システムは、車両の目標とする挙動(このモデルで計算されるような理想的応答)と、実際の挙動(センサーで検知した値)の差を検出して、ブレーキやステアリングトルクで修正をかける。その「目標モデル」として、この2DOFモデルがシンプルで有用なんだ。
意外なところではロボットの自律移動にも応用される。二輪のモデルだけど、四輪ロボットの運動計画でも、車体の向きと進路のズレ(横すべり角)を考慮した制御を行うことがある。また、トラクターとトレーラーのような連結車両の運動安定性解析も、このモデルを拡張して後続車両を追加することで考えられるんだ。
このツールに慣れてきたら、次のステップに進んでみよう。まず数学的には、シミュレーターの背後にある連立微分方程式をラプラス変換して伝達関数を求めることをおすすめする。例えば、ステア入力δに対するヨーレートrの伝達関数 $G_{r\delta}(s) = r(s)/\delta(s)$ を導出してみるんだ。そうすると、応答の「速さ」や「落ち着き具合」を決める固有振動数や減衰係数といった指標が得られて、車両のダイナミクスをより深く評価できるようになるよ。
モデルの拡張としては、「3DOFモデル」を学ぶのが定番コースだ。これは横方向・ヨー方向に加えて、前後方向(駆動力・制動力)の自由度を追加したモデル。これにより、加速やブレーキング時の荷重移動がコーナリングパワーに与える影響(例えば、アクセルを踏み込むとオーバーステア傾向が強まるなど)を考慮できるようになる。これが理解できれば、実車の挙動にまた一歩近づける。
最終的には、この線形領域の知見を土台にして、非線形領域のシミュレーションに挑戦しよう。具体的には、「マジックフォーミュラ」と呼ばれる非線形タイヤモデルを使ったシミュレーションや、閉ループドライバーモデルを組み合わせたコース走行シミュレーションだ。この基礎ツールで感じたパラメータ変化の感覚は、その先のより複雑な世界でも必ず役に立つ、エンジニアの財産になるからね。