シミュレーターのための自動車物理

タイヤのトラクション，ブレーキ，転がり抵抗，空気抵抗，遠心力を使用したシミュレーターを考える．簡単のために自動車は平面上を移動するものとし，段差を考慮しない．

概要

車体にかかる力をすべて計算して，数値計算により速度ベクトルと位置ベクトルとを計算する．

車体にかかる以下の力を計算により求める．太字表記はベクトルを表す．

トラクション F_traction[N]

空気抵抗 F_aero[N]

転がり抵抗 F_rr[N]

遠心力 F_c[N]

ブレーキ F_brake[N]

コーナリングフォース F_cf[N]

これらを足し合わせれば車体にかかっている力が求まる．上下の移動を考慮しないので，重力は計算に入っていない．F_aero と F_rr と F_brake とは F_tractionと逆向きであることに注意．

F = F_traction + F_aero + F_rr + F_c + F_brake + F_cf

力がわかれば，加速度ベクトルは運動の第２法則で計算できる．

a = F/m

車の速度ベクトル v[m/s] は加速度ベクトルを時間に関して積分することで求められる．シミュレーターで積分を計算するには，オイラー法やホイン（Heun）法，ルンゲ＝クッタ法を使う．オイラー法を使う場合は以下のようになる．

v = v + dt・a

dt：前フレームからの経過時間[s]

車の位置 pも同様に計算できる．

p = p + dt・v

加速度の省略

加速度から一度に位置を計算するには以下のようになる．

p = p + v・dt + 0.5・a・dt²

フレームレートが高い場合は加速度の項は非常に小さな値になるため，加速度の項を省くことがある．

空気抵抗

車の空気抵抗は圧力抵抗と摩擦抵抗との和で表現できる．圧力抵抗は車体前後の圧力差によって生じる抵抗で，以下の４つに分類できる．

1. 形状抵抗：車の形によってきまる抵抗．圧力抵抗のほとんどがこれである
2. 凹凸抵抗：ミラーやウイング等の突起物による抵抗
3. 吹き抜け空気抵抗：エンジンの冷却や車内の換気によっておこる抵抗．ラジエータから機関室への空気の吹き抜けの影響が大きい
4. 誘導抵抗：車体を上方へ持ち上げる揚力によって引き起こされる抵抗

摩擦抵抗は車体表面とそれに触れている空気との間の摩擦によって生じる抵抗だ．

空気抵抗におけるそれぞれの要素の影響力の割合は一般車両では以下のようになる．

空気抵抗の成分比率
『自動車工学』自動車工学編集委員会編著 山海堂発行 70頁 図4.2 より転載

これらの要素を単純化すると空気抵抗 F_aero[N] は以下の式になる．

空気抵抗の計算式

空気の密度は以下の式で計算できる．ちなみに 0℃，1 atm の状況での空気の密度 ρ は 1.293 kg/m³である．

空気密度の計算式

以下の表は空気抵抗係数 C_dの概略値である．

自動車の正面投影面積 A が明確でない場合は以下の近似式のどちらかを使う．

転がり抵抗

転がり抵抗は，タイヤの変形，路面の凸凹，タイヤが路面上を部分的に滑るためにおこる摩擦などによって引き起こされる．これを計算するのは困難なので，転がり抵抗はタイヤにかかる荷重に比例すると考え以下の式で表す．

F_rr= u・μ_rW

トラクションがかかっていなくてもタイヤが地面上を転がるなら転がり抵抗は発生する．

転がり抵抗はこのほか車速，タイヤの空気圧，タイヤ構造などによって変わる．車速に比例して転がり抵抗はわずかに上昇する．車速が 140 km/h 以上になるとスタンディングウェーブの発生により転がり抵抗が急増する．タイヤの空気圧増加率の半分くらいの割合で転がり抵抗は減少する．

低速では転がり抵抗が主な抵抗力である．高速では空気抵抗がそれにとって代わる．100 km/h ぐらいで両者が一致する[Car Physics]．

タイヤにかかる荷重

トラクション，ブレーキ，転がり抵抗等を計算するにはタイヤにかかる荷重を計算する必要がある．タイヤにかかる荷重は X 軸回り（ロール）の回転による荷重と，Y 軸周り（ピッチ）の荷重とをかけ合わせれば計算できる．

車の座標系

タイヤにかかる荷重は タイヤにかかる荷重の計算を参照．

ピッチ方向の傾きを φ [rad]，ロール方向の傾きを Θ [rad] とすると，フロントとリアの荷重配分比率は以下のようになる．実際のタイヤの荷重はそれぞれふたつの荷重式と車の重量[N]とを掛け合わせばいい．４輪の荷重をすべて合計すると車の重量と一致する．

荷重の分配比率

タイヤにかかる荷重の数式

進行方向への加速度は減速時には符号を反転させる．旋回半径は遠心力の項目を参照．

タイヤの接地判定

タイヤが接地していないとき，その荷重を０で計算する．

トラクション

エンジンがタイヤを回転させるトルク [N・m] はエンジンの回転数，エンジンの出力，現在のギア，タイヤの半径，スロットル等に影響される．エンジンが車に与える速度変化 [N] は以下の式で表現できる．

F_traction = u・T_engine・r_t・x_g・x_d・n / R_tire

トランスミッション効率は MT 95%以上，CVT 低速＝85% 高速＝75%，AT 低速＝80% 高速＝90% と言われている．

F_tractionの大きさがタイヤの μ_sW を超えるとタイヤはスリップする．この μ_s はタイヤと路面との摩擦係数で転がり抵抗係数とは別のものである．μ_s についてはタイヤのスキッド限界を参照．

ギア比やデフ比を大きくすればいくらでもトルクは大きくなる．しかしそうするとタイヤの角速度――つまり最高速度――が低下する．ギア比やデフ比を大きくするとすぐにエンジンの回転数の上限に達し，加速できなくなる．

エンジンのトルクは 1000 rpm もしくは 500 rpm 区切り（望むならもっと細かくてもいい）の配列で表現する．現在の回転数に応じてそれを線形補間してトルクを計算する．

エンジンの RPM の計算

タイヤが接地しているならば，エンジンの回転数 [rpm] は車速やタイヤ半径などの影響を受ける．現在の車速 v [m/s] ，タイヤの半径 r [m]，ギア比 x_g，デフ比 x_dとすると，エンジンの回転数 [rpm] は以下の式で計算できる．

rpm = v/(2・π・r)・60・x_g・x_d

大抵のエンジンは 1000 rpm 以下の回転数の時のトルクを定義していない．シミュレーターがエンストを許さないなら，そのような状態になったときはクラッチを切りアイドリング状態に移行する処理を追加する．

スリップ比とトラクションとの関係

スリップ比は摩擦係数だけでなくトラクションによっても変化する．現在のタイヤのトルクからスリップ比を取得し，RPM を修正する．

『Calculating Longitudinal Wheel Slip and Tire Parameters Using GPS Velocity』Shannon L. Miller ほか著 1804項 図4 より転載

このグラフを数式にすると，スリップ比 = F_traction/2000 になる．トラクションが 4000[N] のとき，スリップ比は 2％．よって RPM・1.02 が実際の RPM になる．

遠心力

ハンドルを切ると車体には遠心力 F_c[N] がかかる．これは速度の２乗に比例して大きくなる．速度が速いほどコーナーが曲がりにくくなる原因の一つが遠心力である．

旋回半径を R [m]，角速度を ω [rad/s] とすると遠心力は ω²R [m/s²]．ωR = V なのでこれは V²/R と書き換えられる．重心に働く遠心力 [N] は以下の式で計算できる．

旋回半径はアッカーマンアングルを利用して求める．これは遠心力が発生しないぐらい低速で旋回するときのモデルであるため，定常旋回時の旋回半径とは異なる．

『自動車工学』自動車工学編集委員会編著 山海堂発行 127頁 図7.10 より転載

R₁を旋回半径 R [m] として使う場合，旋回半径は以下の式で計算できる．

車体のロール/ピッチ角度の計算

遠心力を受けたり，加減速したりすると車体にモメントが発生し，回転運動を生じる．そのときの車体の傾きは，車高を下げたり，サスペンションのばね乗数を大きくしたりして調整できる．

ロール角度

サスペンションは追加のモメント ΔW [N・m]を受け，その分だけ沈む．ΔW は F_c・h + W・h・tanΘ だ．サスペンションの反発力 w_k [N]は k・x で表現できる；x はサスペンションが沈み込んだ距離[m]で，k はばね乗数[N/m]だ．サスペンションが沈み込んだ距離 x は (t/2)・tanΘで，Θ << 1 のとき tanΘ ≒ Θ[rad] で近似できて，x=tΘ/2 になる．

最終的なサスペンションの反発モメント[N・m] は t²kΘ/2 になる．F_c・h は遠心力によるモメントで，W・h・tanΘ は車体の傾きで生じる重力によるモメントだ．結局，車体のロール角度[rad] は以下の式になる．

ただしばね乗数 k は前後輪のばね乗数の合計で，前輪のばね乗数を k_f，後輪のばね乗数を k_r とすると，k = k_f + k_r である．

この t²k/2 の項は一般的にロール剛性[N/rad] と呼ばれる．

シャシー形式ごとのロールセンターの求め方はロールセンターの求め方を参照．

加減速時のピッチ角度も同様に計算できて，以下の式になる．

コーナリングフォース

横滑り角

横滑り角は実舵角と車の進行方向との間の角度である．現在フレームの車の進行方向を計算するのに横滑り角が必要なので，ひとつ前のフレームの車の進行方向と現在フレームの実舵角とから横滑り角を計算する．

ステアリングを切り，横滑り角が発生すると遠心力の働く方向にタイヤがねじれる．このときのタイヤが元の形に戻ろうとする力がコーナリングパワーだ．コーナリングパワー [N/rad] はタイヤに固有のコーナリング特性で， 単位横滑り角あたりのコーナリングフォース[N]を示している．コーナリングパワーはタイヤ荷重やタイヤ内圧，横滑り角の増加とともに増加する．ただし横滑り角がある閾値を超えると，コーナリングフォースは減少し頭打ちが発生する．引用した図（自動車の運動と制御）では，横滑り角が約４°までは直線的にコーナリングフォースが増加する．しかしそこから次第にコーナリングフォースは減少し，横滑り角が約１０°でコーナリングフォースは飽和している．コーナリングフォースはそのほかにもタイヤのプライ数，コードクラウン角度，ブレーカ，トレッドパターン，アスペクト比に影響を受ける．

『自動車工学』自動車工学編集委員会編著 山海堂発行 123頁 図7.5 より転載

『自動車工学』自動車工学編集委員会編著 山海堂発行 122頁 図7.3，7.4 より転載

『自動車の運動と制御[第２版]』安部正人著 東京電機大学出版局発行 21頁 図2.16 より転載

コーナリングフォース F_c[N] はコーナリングパワー[N/rad]と横滑り角とから計算できるが，それを単純な関数で表現するのは難しい．いくつかの直線で近似するか，テーブルにデータを保存しておき参照するかして取得する．

横滑り角からコーナリングパワーを計算する関数を f_cp(α) とするとコーナリングフォースは以下の式で計算できる．

F_cf = T・f_cp(α)

コーナリングパワーの働く方向 T は車体の進行方向と世界の上向きベクトル（たいていは（0, 1, 0））との外積をとることで計算できる．車体の進行方向は，前フレームに計算した速度ベクトルを正規化したものを使う．

ブレーキ

タイヤのスキッド限界

トラクションやコーナリングフォースがタイヤのスキッド限界を超えるとスリップする．コーナー中にトラクションを掛けたり，ブレーキしたりするとその力と遠心力とでタイヤの制駆動力を使用してしまい，コーナーを曲がるのに使うコーナリングフォースが減少する．タイヤのスキッド限界はタイヤにかかる荷重，タイヤの温度，タイヤの損耗具合，タイヤの摩擦係数に依存する．近似としては F_skid[N] は μ_sW で近似する．W はタイヤにかかる荷重．μ_s は摩擦係数で，乾いたアスファルトでは 0.7～0.8，濡れていると 0.4～0.5，雪道では 0.1～0.2 ぐらいになる．

ブレーキの発生させる制動力

ブレーキの発生させる制動力とタイヤのスキッド限界とは別物である．たとえばタイヤのスキッド限界が 4000[N] だとしよう．この車に 6000[N] の制動力をかけても無意味だ；そうしても単にタイヤがロックされてスリップするだけだからだ．もちろんスキッド限界をわずかに超える制動力を設定するのは問題ない．

制動力はタイヤの半径とタイヤの摩擦係数とタイヤの空気圧とに依存する．タイヤの半径を大きくすれば最高速が伸びる．しかしブレーキ時の地面がタイヤを引っ張る力も大きくなり，それだけ強力なブレーキが必要になる．空気圧が高ければ接地面積が狭くなり制動力はおちる，逆に空気圧が低ければ制動力は増す．

制動力 F_b[N] は以下の式で表現できる．

F_b = u・r_brake・B_max

コーナー中にトラクションを発生させたり，制動力をかけたりするとその分だけコーナリングフォースが減少する．

前輪と後輪とのブレーキの強さは普通同じにしない．減速時には前輪により多く荷重がかかるため前輪はブレーキの効きが強く，後輪はブレーキの効きが弱くなる．そのため普通は前輪のブレーキ力を強く，後輪のブレーキ力を弱く設定する．

車速が 0 付近のとき，ブレーキに特別な処理が必要になる．そうしないと車が勝手に後退してしまう．

ある減加速度 a [m/s²] が必要な時の各タイヤのブレーキ力 B_x[N] は W_xをそれぞれのタイヤの荷重とすると，以下の式で計算できる．

B_x = W_x・a/g

路面の温度が低いとタイヤが温まりにくく，グリップが悪くなる．

まとめ

車オブジェクトが持っている定数

車オブジェクトが持っている状態変数

ゲームでは１フレームの間にステアを非現実的な速度で操作できるため， それを補正するために現在のハンドルの回転度合を保存しておく．

外部から渡されるパラメータ/環境が持っている変数

ユーザーからの入力

計算過程

1. エンジンの踏み込み量とハンドルの回転角度を取得
2. トラクションを計算
3. 空気抵抗を計算
4. 転がり抵抗を計算
5. 遠心力を計算
6. コーナリングパワーを計算
7. 制動力を計算
8. 車にかかる力を計算
9. 加速度を計算
10. 速度を計算
11. 位置を計算
12. 速度からエンジンの回転数を再計算

参考文献

[自動車工学編集委員会]自動車工学

[安部正人]自動車の運動と制御

Car Physics

VehicleDynamicsForRacingGames（pdf）

Calculating Longitudinal Wheel Slip and Tire Parameters Using GPS Velocity（pdf）