パラメータ設定
統計サマリー
理論式
ピストン位置：
$$x(\theta) = r\cos\theta + \sqrt{L^2 - r^2\sin^2\theta}$$
速度（近似）：
$$v \approx -r\omega\!\left(\sin\theta + \frac{r}{2L}\sin 2\theta\right)$$
加速度（近似）：
$$a \approx -r\omega^2\!\left(\cos\theta + \frac{r}{L}\cos 2\theta\right)$$
CAEへの応用：内燃機関のクランク機構では、ピストン加速度が慣性力を決め、軸受荷重・クランクシャフト疲労設計に直結します。R/L比は運動の非調和性を支配する重要なパラメータです。
位置 x(t)
速度 v(t)
加速度 a(t)
位相 θ 対 位置
クランクスライダー機構とは
🧑🎓
このシミュレーターで動いてる仕組み、自動車のエンジンみたいですね。ピストンがどうしてあんな動き方をするんですか？
🎓
そうだね、まさにエンジンの基本モデルだ。ざっくり言うと、クランクの回転運動を、コネクティングロッドを介してピストンの直線運動に変換してるんだ。上のスライダーで「クランク半径」を大きくしてみて。ピストンの動く範囲（ストローク）が大きくなるのがわかるよね？
🧑🎓
え、そうなんですか！確かに動く距離が伸びました。でも、グラフを見ると速度や加速度が綺麗なサイン波じゃないみたいです。なんで歪んでるんですか？
🎓
良いところに気づいたね。コネクティングロッドが斜めになることで生じる「非調和性」が原因だ。今度は「R/L比」を大きくしてみて。ロッドが短くなると、上死点付近でのピストンの動きが遅くなって、加速度の山が尖るのがわかる。実務ではこの非対称性が軸受の負荷に直結するんだ。
🧑🎓
なるほど！じゃあ、エンジン設計ではこの歪みを考慮するんですね。シミュレーターで「角速度」を上げると、加速度がめちゃくちゃ大きくなりますが、これが振動や騒音の原因ってことですか？
🎓
その通り！ピストン加速度 $a$ は角速度 $\omega$ の2乗に比例して大きくなる。高回転エンジンではこの慣性力が爆発的に増えて、軸受の摩耗やシャフトの疲労破壊の主原因になる。CAEではこの加速度の波形から各部品にかかる力を精密に計算して、軽量化と強度を両立させる設計をしてるんだ。
物理モデルと主要な数式
ピストンの位置は、クランクの回転角 $\theta$ を使って幾何学的に求められます。コネクティングロッドの傾きを考慮した正確な式です。
$$x(\theta) = r\cos\theta + \sqrt{L^2 - r^2\sin^2\theta}$$
$x$: ピストン位置（基準点からの変位）, $r$: クランク半径, $L$: コネクティングロッド長, $\theta = \omega t + \phi_0$: クランク角
速度と加速度は、位置を時間で微分すれば得られますが、式が複雑になるため、実用上は $r/L$ の項で展開した近似式がよく使われます。これがシミュレーターで表示される波形の元になっています。
$$v \approx -r\omega\!\left(\sin\theta + \frac{r}{2L}\sin 2\theta\right)$$
$$a \approx -r\omega^2\!\left(\cos\theta + \frac{r}{L}\cos 2\theta\right)$$
$v$: ピストン速度, $a$: ピストン加速度, $\omega$: クランク角速度。第2項 ($\sin 2\theta, \cos 2\theta$) が「非調和項」で、波形を歪ませる物理的な原因です。$r/L$ が大きいほどこの影響が強まります。
実世界での応用
自動車・船舶の内燃機関：ガソリンエンジンやディーゼルエンジンの心臓部です。ピストンの加速度から慣性力を算出し、クランクシャフトやコンロッドの強度設計、バランスウェイトの配置、軸受の寿命予測に利用されます。高回転化に伴う慣性力の増大は設計上の大きな課題です。
往復動圧縮機・ポンプ：冷蔵庫のコンプレッサーや工場の空気圧縮機などに使われます。ピストンの動きを正確に把握することで、吐出流量の脈動を低減し、効率的で騒音の少ない設計が可能になります。
蒸気機関：産業革命を支えた原動力です。現代では歴史的な再現や教育目的で解析されます。大きなストロークと低速回転が特徴で、トルク変動の平滑化が重要なテーマでした。
CAEによる機構解析：MBD（多体動力学）ソフトウェアを用いた仮想プロトタイピングの基礎モデルです。この単純な機構の解析結果を検証用ベンチマークとして使い、より複雑なリンク機構やロボットアームの動的シミュレーション精度を担保します。
よくある誤解と注意点
まず、「R/L比は小さければ小さいほど良い」と思い込むことです。確かに、ロッド長Lを長く（R/L比を小さく）すると、ピストン速度・加速度の波形は正弦波に近づき、振動面では有利です。しかし、エンジン全体の高さや重量が増大し、摩擦損失も増えます。実務では、エンジン全体のパッケージングと性能のトレードオフの中で、R/L比=1/3.5〜1/4.5あたりがよく選ばれます。例えば、クランク半径r=45mmなら、コンロッド長Lは160mm〜200mm程度が典型的です。
次に、シミュレーターの「角速度」を現実のエンジン回転数と単純に結びつけること。ツールで角速度ωを上げると加速度が急増しますが、実際のエンジン設計では「平均ピストン速度」という指標がより重要視されます。式は $V_m = 2 \times S \times N$ （S:ストローク、N:回転数）で、これが材料や潤滑の限界を決めます。例えば、市販車用ガソリンエンジンでは$V_m$が20m/sを超えないように設計されることが多いです。
最後に、静的なバランスだけを考えてしまう点。クランクシャフトの静バランス（重りで回転のつり合いを取ること）はもちろん重要ですが、往復運動するピストンとコンロッドの慣性力による偶力（シェイキングフォース）がエンジンブロックを揺さぶります。この振動を低減するためには、バンク角を付けたV型エンジンや、バランスシャフトの採用といった対策が必要になります。シミュレーターで計算される加速度は、この「揺らす力」を計算する第一歩なのです。
関連する工学分野
このシミュレーターの計算は、機構学（メカニズム・運動学）の基本ですが、その先には広大な工学分野が広がっています。まずは振動・騒音工学（NVH）。ピストンの加速度から慣性力が求められ、これがエンジンブロックを介して車体に伝わり「ブーン」という騒音や振動の根源になります。CAEでは、この力を入力としてモード解析や周波数応答解析を行い、どの回転数でどの部品が共鳴するかを予測します。
次に、材料力学・疲労強度解析。コンロッドやクランクピンには、燃焼圧力に加え、シミュレーターで見た巨大な加速度による慣性力が繰り返し加わります。この交番応力が、金属の疲労破壊を引き起こします。CAEでは、この動的荷重を時間歴として与え、部品の危険個所や寿命を疲労解析ソフトウェアで詳細に評価します。
さらに流体力学（CFD）とも深く関わります。ピストンの動き（位置、速度）は、シリンダー内のガス交換（吸気・排気）の効率を左右します。特に上死点付近でのピストン速度が遅くなる（R/L比が大きいとより顕著）と、燃焼室の乱流形成に影響し、燃焼効率やノッキング特性が変化します。エンジンシミュレーションでは、このピストン運動を境界条件として与え、内部流れを計算するのです。
発展的な学習のために
まずは、「微分」の物理的意味をグラフと一緒に体感することから始めましょう。シミュレーターで位置→速度→加速度と連続的に表示されるのは、まさに「位置を時間で微分すると速度、速度を時間で微分すると加速度」という関係を視覚化したものです。この感覚が身につけば、振動（変位の2階微分）やジャーク（加速度の微分）といったより高次の物理量も理解しやすくなります。
次に、近似式の導出を追ってみることをお勧めします。ピストン位置の正確な式 $$x(\theta) = r\cos\theta + \sqrt{L^2 - r^2\sin^2\theta}$$ の平方根部分を、二項定理を使って $r/L$ のべき級数に展開します。例えば、$\sqrt{L^2 - r^2\sin^2\theta} = L\sqrt{1 - (\frac{r}{L}\sin\theta)^2} \approx L(1 - \frac{1}{2}(\frac{r}{L}\sin\theta)^2 + \cdots)$ のように変形し、$\theta$の関数として整理すると、速度・加速度の近似式が導けます。この「非線形項を線形化・級数展開する」手法は、CAEの世界で頻繁に使われるアプローチです。
このツールに慣れたら、次のステップとして多気筒エンジンの一次・二次慣性力のバランスや、クランクシャフトの動的トルク変動について学ぶと良いでしょう。単気筒の動きを理解した上で、それらが複数気筒でどう合成され、打ち消し合うかを考えることで、実際のエンジン設計の面白さと難しさに触れることができます。これらは全て、このクランクスライダーという基本モデルの延長線上にある世界です。