Main Content

このページの翻訳は最新ではありません。ここをクリックして、英語の最新版を参照してください。

自動車のサスペンション

この例では、前後独立垂直サスペンションを含む単純化した 1/2 車両モデルのモデル化方法を示します。このモデルには、車体のピッチおよびバウンスの自由度も含まれています。この例では、乗り心地特性を調べるためにシミュレーションをどのように使用できるかを示すモデルについて説明します。このモデルをパワートレイン シミュレーションと共に使用して、スロットル設定の変更によって生じる縦方向の移動を調べることができます。

解析と物理学

1/2 車両モデルの自由体図

図は、1/2 車両のモデル化された特性を示しています。フロントおよびリア サスペンションは、バネ/ダンパー システムとしてモデル化されています。より詳細なモデルには、タイヤ モデルや、(圧縮時よりリバウンド時に減衰が大きくなる) 速度依存減衰などのダンパーの非線形性が含まれます。車体はピッチおよびバウンスの自由度をもちます。これらは、モデル内で 4 つのステート、つまり、垂直変位、垂直速度、ピッチ角変位、ピッチ角速度で表されます。6 つの自由度をもつ完全なモデルを実装するには、ベクトル代数ブロックを使用して、軸変換と、力/変位/速度の計算を行います。方程式 1 は、フロント サスペンションがバウンス (つまり、垂直の自由度) に与える影響を記述します。

$$F_{f} = 2K_f (L_f \theta - (z + h)) + 2C_f(L_f \dot{\theta} -\dot{z})$$

ここで、

$$F_{f}, F_{r} = \mbox{ upward force on body from front/rear suspension}$$

$$K_f, K_r = \mbox{ front and rear suspension spring constant}$$

$$C_f, C_r = \mbox{ front and rear suspension damping rate}$$

$$L_f, L_r = \mbox{ horizontal distance from gravity center to front/rear suspension}$$

$$\theta, \dot{\theta} = \mbox{ pitch (rotational) angle and its rate of change}$$

$$z, \dot{z} = \mbox{ bounce (vertical) distance and its rate of change}$$

$$h = \mbox{ road height }$$

方程式 2 は、サスペンションによるピッチのモーメントを記述します。

$$M_{f} = -L_{f}F_{f}$$

$$F_{r} = -2K_r (L_r\theta + (z + h)) -2C_r ( L_r \dot{\theta} + \dot{z})$$

$$M_{r} = L_r F_{r}$$

ここで、

$$M_{f}, M_{r} = \mbox{ Pitch moment due to front/rear suspension}$$

方程式 3 は、ニュートンの第 2 法則に従って、車体運動に生じる力とモーメントを解決します。

$$m_b\ddot{z} = F_{f} + F_{r} - m_b g$$

$$I_{yy} \ddot{\theta} = M_{f} + M_{r} + M_y $$

ここで、

$$ m_b = \mbox{ body mass}$$

$$ M_y = \mbox{ pitch moment induced by vehicle acceleration}$$

$$I_{yy} = \mbox{ body moment of inertia about gravity center}$$

モデル

モデルを開くには、MATLAB® コマンド ウィンドウに sldemo_suspn と入力します。

サスペンション モデルの最上位レベルのブロック線図

サスペンション モデルは 2 つの入力をもち、両方の入力ブロックがモデル ブロック線図上で青色で表示されています。最初の入力は道路高です。ここで入力されるステップは、高さが階段状に変化する路面を走行する車両に対応します。2 番目の入力は、ブレーキまたは加速操作によって生じる、車輪の中心を介して作用する水平力です。縦方向の車体運動はモデル化されていないため、この入力はピッチ軸の周りのモーメントとしてのみ表示されます。

FrontSuspension および RearSuspension サブシステムで使用されるバネ/ダンパー モデル

フロントおよびリア サスペンションをモデル化するバネ/ダンパー サブシステムを上図に示します。Front/Rear Suspension ブロックを右クリックし、[マスク][マスク内を表示] を選択してフロント/リア サスペンション サブシステムを確認します。サスペンション サブシステムは、方程式 1-3 をモデル化するために使用されます。これらの方程式は、Gain および Summation ブロックを直接使用することによって、Simulink® ブロック線図で直接実装されます。

フロントとリアの違いは次のように考慮されます。このサブシステムはマスク ブロックであるため、インスタンスごとに異なるデータセット (LK および C) を入力できます。さらに、L は、直交座標 x と見なされ、原点または重心に対して負または正になります。したがって、KfCf および -Lf がフロント サスペンション ブロックに使用されるのに対して、KrCr および Lr がリア サスペンション ブロックに使用されます。

シミュレーションの実行

このモデルを実行するには、[シミュレーション] タブで [実行] をクリックします。初期条件が sldemo_suspdat.m ファイルからモデル ワークスペースに読み込まれます。モデル ワークスペースの内容を確認するには、Simulink エディターの [モデル化] タブで、[設計] から [モデル エクスプローラー] を選択します。モデル エクスプローラーで、sldemo_suspn モデルの内容の下で [モデル ワークスペース] を選択します。モデル ワークスペースに初期条件を読み込むことで、パラメーターを誤って変更するのを防ぎ、MATLAB ワークスペースをクリーンな状態に保つことができます。

このモデルは、MATLAB ワークスペースの sldemo_suspn_output というデータ構造に関連データのログを作成することに注意してください。構造の名前を入力して、どのようなデータが含まれているかを確認します。

シミュレーション結果

シミュレーション結果が上に表示されています。結果は、sldemo_suspgraph.m ファイルによってプロットされます。既定の初期条件については、表 1 を参照してください。

表 1: 既定の初期条件

Lf = 0.9;    % front hub displacement from body gravity center (m)
Lr = 1.2;    % rear hub displacement from body gravity center (m)
Mb = 1200;   % body mass (kg)
Iyy = 2100;  % body moment of inertia about y-axis in (kg m^2)
kf = 28000;  % front suspension stiffness in (N/m)
kr = 21000;  % rear suspension stiffness in (N/m)
cf = 2500;   % front suspension damping in (N sec/m)
cr = 2000;   % rear suspension damping in (N sec/m)

モデルを閉じる

モデルを閉じ、MATLAB ワークスペースから生成されたデータを削除します。

まとめ

このモデルでは、サスペンションの減衰と剛性を変更する影響をシミュレートし、それにより、乗り心地と性能間のトレードオフを調べることができます。一般に、レーシング カーはバネ剛性が非常に高く、減衰係数が高いのに対して、乗用車のバネはより柔らかく、振動応答性がより高くなります。

関連するトピック