CNC ミリング再生ビビり安定限界シミュレーター

エンドミル加工で発生する再生ビビりの安定限界を、Tlusty/Altintas の解析モデルで計算します。主軸回転数・刃数・工具固有振動数・動的剛性を変えて安定ローブ線図 (Stability Lobe Diagram) を描き、最大の軸切込みが取れる sweet spot を探す加工計画ツールです。

パラメータ設定

主軸回転数 N

rpm

現在の動作点。lobe 線図上に赤×印で表示

刃数 z

工具固有振動数 f_n

Tap test で測定する工具先端の1次共振周波数

動的剛性 k

N/m

工具先端のモード剛性。突き出しを短くすると激増

減衰比 ζ

一般工具で 0.02〜0.05、ダンパー付きで 0.08+

比切削力 K_t

N/mm²

Al≈800, 鋼≈2000, チタン≈3000, 焼入鋼≈4000

半径方向接触率 a_e/D

フル溝切りで 1.0、仕上げで 0.05〜0.10

計算結果

—

安定限界軸深 b_lim_min (mm)

—

歯通過周波数 f_TP (Hz)

—

チャター周波数 f_c (Hz)

—

最適主軸回転数 (rpm)

—

材料除去率 (cm³/min)

—

動的剛性 k (N/m)

—

加工状態可視化 — エンドミル + チップ厚変動

エンドミルが回転し、各刃が前刃の振動的に削った表面と重なります。色は現動作点が安定 (青) か不安定 (赤) かを示します。

安定ローブ線図 — b_lim vs 主軸回転数 N

工具周波数応答 — |G(f)| と Re[G(f)]

理論・主要公式

$$b_{lim,min} = \frac{-1}{2\,K_t\,G_{min}\,\rho},\qquad G_{min} = \frac{-1}{4\,k\,\zeta\,(1+\zeta)}$$

安定限界軸切込み b_lim_min は工具周波数応答 G の負の最小値で決まる (Tlusty/Altintas)。K_t：比切削力、ρ：半径方向接触率、k：動的剛性、ζ：減衰比。

$$G(j\omega) = \frac{1/k}{1 - r^{2} + 2\,i\,\zeta\,r},\qquad r = \frac{\omega}{\omega_n}$$

単一モード工具の周波数応答関数 (FRF)。r=1 で共振、Re[G]<0 領域 (r>1 弱辺) でビビりが励起される。

$$f_{TP} = \frac{N\,z}{60},\qquad N_{opt} = \frac{60\,f_n}{z\,(n+1)},\quad n=0,1,2,\dots$$

歯通過周波数と sweet spot 回転数。最も低次の n=0 lobe は f_TP=f_n で発生し、最大の安定切込みが得られる。

CNC ミリング再生ビビり安定限界 — Stability Lobe Diagram

🙋

エンドミル加工で「ビビり」って言葉、現場でよく聞きます。あれって普通の振動と何が違うんですか？回転数を上げるとむしろ消えるって聞いて、よく分かりません。

🎓

いい疑問だね。CNC で出るあのキーンという嫌な音とパキパキの加工面、あれは「再生ビビり (Regenerative Chatter)」っていう自励振動なんだ。仕組みは：前の刃が振動しながら表面に波模様を残す→次の刃が来たとき、自分の振動と前刃の波の位相関係でチップ厚が変動する→切削力が変動して工具がもっと揺れる→さらに波が深くなる…という正帰還ループ。だから普通の強制振動と違って、いったん条件が悪いと数秒で爆発的に成長して、面粗度が一気に Rz 50μm 超えとかになる。Tlusty が1965年に解明したんだ。

🙋

なるほど、前の刃が残した波と現在の刃の関係なんですね。じゃあ「回転数を上げると消える」ってのも、その位相関係で説明できる…？

🎓

そう、まさにそれが「安定ローブ線図」の核心。歯通過周波数 f_TP = N·z/60 を工具のチャター周波数 f_c とちょうど整数比にあわせると、前刃の波と現刃の振動が同位相になってチップ厚変動が打ち消されるんだ。例えば f_n=800Hz、z=4枚なら、N=12000rpm で f_TP=800Hz=f_n となって「sweet spot」になる。このとき軸切込みを 3〜5倍に取っても安定する。逆にその間の N=9000rpm 付近では位相が最悪で、ローブの「壁」になる。右のグラフでスライダーを動かすと、その谷と壁が見えるよ。

🙋

じゃあ闇雲に高速回転すればいいわけじゃなくて、その「谷」を狙う必要があるんですね。デフォルトの5000rpmって、谷から外れてるってことですか？

🎓

そう、デフォルト N=5000rpm だと f_TP=333Hz で f_n=800Hz の半分以下。一番低い n=0 lobe からも外れた「壁」付近で、b_lim_min が支配する保守的な領域だ。本ツールでも判定が warn になるはず。実務ではまず Tap test で f_n を測り、N_opt = 60·f_n/z = 12000rpm を狙う。Al 高速加工なら 20000rpm 超まで使うこともある。ただし高速化すると今度は工具寿命や主軸の動的剛性低下が新しいボトルネックになるんだ。

🙋

動的剛性 k って左のスライダーにありますね。これって普通の静剛性と違うんですか？値を上げると b_lim がぐっと伸びますけど。

🎓

大事な区別だね。静剛性は「ゆっくり押して何 N/mm 沈むか」だけど、動的剛性は工具先端を共振周波数 f_n で揺らしたときの応答振幅から決まる。一般に動的剛性 ≪ 静剛性で、長い工具・細工具では10倍以上落ちる。改善策は①突き出し量を最小化 (剛性は L³ に比例)、②ヒートシュリンクや HSK で接合剛性向上、③ダンパー付き工具で ζ を 0.03→0.08 に増やす、の3つ。チタン薄壁のブレード加工や航空機ハニカム加工では、k と ζ を稼ぐためにダンパー入り超硬工具を必ず使うよ。

🙋

最後にひとつ：比切削力 K_t を増やすと b_lim が減るのは直感的に分かるんですが、これってどう測るんですか？

🎓

K_t は「単位チップ厚あたりの切削力」で、材料と工具すくい角で決まる。実測は切削力動力計 (Kistler 等) で F_t と F_n を測り、チップ厚で割って同定する。教科書値の目安は Al6061≈800、軟鋼≈1800〜2200、SUS304≈2400、Ti6Al4V≈2800〜3200、焼入鋼SKD11≈3500〜4000 N/mm²。チタンが難加工材と言われるのは K_t が高い上に熱伝導率が低くて工具が熱劣化するから。本ツールでも材料を変えてどう sweet spot が動くか確かめてみて。

よくある質問

前の刃が削った表面の波形と、現在の刃の振動的な動きが位相差をもって重なることで、切りくず厚さが周期的に変動し、これが切削力を変動させて工具を更に振動させる自励振動です。Tlusty (1965) と Tobias が解析モデルを確立し、Altintas が周波数領域で一般化しました。摩擦由来のスティックスリップとは異なり、切削プロセス自体が増幅器として働くため、いったん発生すると数秒で爆発的に成長し、加工面が縞模様 (chatter mark) になり工具寿命を激減させます。安定ローブ線図はこの現象を主軸回転数と軸切込みの平面で表したものです。

歯通過周波数 f_TP = N·z/60 が工具のチャター周波数 f_c とちょうど整数比になる主軸回転数では、前刃が削った表面波と現刃の振動が同位相になり、切りくず厚変動が打ち消されて安定切込みが急激に増えます。これが「谷」つまり sweet spot で、その間の高い壁 (lobe boundary) では位相がずれて再生効果が最大となり不安定領域となります。最も低い lobe (n=0) は f_TP ≈ f_c のとき、すなわち N_opt = 60·f_n/z で発生し、ここでは安定切込みが最小値の2〜5倍に達することがあります。高速加工はこの sweet spot を狙う技術です。

標準的な方法は Tap test (impulse hammer test) です。主軸に工具を装着した状態で、インパクトハンマーで工具先端を XY 方向に叩き、加速度センサで応答を測定します。FFT で周波数応答関数 (FRF) を求め、共振ピークから固有振動数 f_n、半値幅から減衰比 ζ、ピーク高さから動的剛性 k = 1/(2ζ·|G(f_n)|) を抽出します。市販ソフトでは CutPro、MetalMAX、Harmonizer 等が有名で、5〜30 分で機械・工具の組合せ毎に測定できます。工具が変われば FRF も変わるため、長い工具・細い工具・テーパー違いではそれぞれ Tap test を行うのが理想です。

(1) 主軸回転数の最適化—sweet spot にあわせる (本ツールの主目的)。(2) 軸切込みを下げる—最も確実だが MRR が落ちる。(3) 不等ピッチ工具・可変リード工具を使う—周波数を分散させて再生効果を弱める。(4) 工具突き出し量を短くする—剛性は L³ に比例して向上。(5) ホルダーをヒートシュリンクや HSK に変える—接合剛性向上。(6) ダンパー付き工具 (内部質量ダンパー) を使う—航空機チタン加工で多用。(7) 切削液の量と方向を調整—摩擦由来の励振を抑える。最初に試すべきは (1) と (4) で、これだけで多くのケースが解決します。

実世界での応用

航空機構造部品の薄壁加工：翼の spar、ブレード、フレームのポケット加工では、最終肉厚が 1〜2mm の薄壁を長い工具で切削します。動的剛性が極端に低く、ビビり発生がボトルネックです。実機では Tap test で求めた FRF をもとに CAM (Mastercam、NX CAM) が自動で N と b を最適化し、ローブ図の谷を歩く NC プログラムを生成します。Al7050 で MRR 20→80 cm³/min への改善例が多数報告されています。

金型・パンチダイスの仕上げ加工：SKD11/SKH51 (HRC 60+) の高硬度材を、小径ボールエンドミル (Ø6〜Ø3) で仕上げる工程です。K_t=4000 と工具細さが重なりビビりが避けられません。本ツールでも radialEngagementRatio を 0.05〜0.10 まで下げ、主軸回転を 25000rpm 以上に上げて高次 lobe の谷を狙う設計をします。Side step の最適化は、Altintas モデルから自動算出する CutPro 等で行います。

医療デバイス・人工関節：Ti6Al4V や CoCrMo の人工骨頭・人工関節は、形状精度 ±10μm のもと、長い突き出し量で 5軸加工されます。チタンは K_t≈3000、熱伝導 6.7 W/mK と難加工で、ビビりが工具破損に直結します。Sandvik や Mitsubishi のダンパー入り超硬エンドミルで ζ を 0.08+ に高め、突き出しを最小化する戦略が標準です。

自動車量産部品 (エンジンブロック、ミッションケース)：Al-Si 合金 (ADC12) のロボットライン加工では、サイクルタイムが量産コストを支配します。本ツールが示すように、Al は K_t≈800 と低いため、剛性さえあれば 15000〜20000rpm の高速・大切込みで MRR 200 cm³/min 級が可能。MAG/MAZAK の最新加工機ではビビり検出センサと AI による主軸回転数の自動補正 (Active Chatter Suppression) も実装されています。

よくある誤解と注意点

第一に、「安定ローブ線図さえあれば全て解決する」と思い込むのは危険です。Tlusty/Altintas モデルは単一モード・単一方向・線形仮定の解析モデルで、実機では (a) 工具と部品の両方が振動モードを持つ「結合系」、(b) 高速回転で工具固有振動数が遠心力で変わる、(c) 過大切込みで非線形 (チップ分離) になる、(d) 多刃工具の各刃の力が異なる、といった現実が無視されます。本ツールの予測は「秒スケールでビビりが収束するか否かの一次近似」と捉え、必ず実機での音響・加速度センサによるリアルタイム監視と組合せてください。Sandvik CoroPlus や Bluestreak Production Studio のようなシステムが市販されています。

第二に、「Tap test を一度だけやって満足する」こと。動的剛性 k と固有振動数 f_n は、工具の突き出し量・ホルダー種類・主軸軸受の摩耗・温度で大きく変化します。新品ホルダーで測った FRF を、半年使った状態の機械にそのまま当てはめると 10〜30% ずれることはザラです。定期的 (3ヶ月 or 工具交換ごと) に Tap test を再実施し、CAM 側のローブ図を更新する運用が必要です。また同じ工具長でもチャック内突き出し量が 5mm 変わると f_n は数十 Hz ずれます。

第三に、「sweet spot に乗せれば必ず生産性が上がる」と無条件に信じること。最適主軸回転数 N_opt は機械主軸の許容上限・冷却能力・工具コーティングの熱限界を超える場合があります。例えば本ツールで f_n=2500Hz、z=2 のとき N_opt=75000rpm が出ても、主軸上限が 20000rpm ならそこで頭打ちです。また高速回転は工具すり減りが指数的に増え、工具コスト/分が悪化する場合もあります。MRR の向上分と工具コスト・電力コストの収支を、必ず別計算で評価しましょう。本ツールは「物理的安定性」を示すもので、「経済的最適」とは別物です。

CNC ミリング再生ビビり安定限界シミュレーター

CNC ミリング再生ビビり安定限界 — Stability Lobe Diagram

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点