「なぜ、ここから壊れたのか」——Ktが教えてくれること
段付き軸のフィレット部から亀裂が入った。穴あき板の穴の縁から破断した。——機械設計や材料力学に関わる人なら、一度は「局所的に壊れる」現象に出くわしたことがあるはず。壊れた箇所を調べると、計算上の平均応力よりはるかに大きな力がかかっていたことがわかる。その正体が応力集中であり、その倍率を示す数値が応力集中係数Ktだ。
このツールは、段付き丸棒・穴あき平板・U溝付き棒など6形状×引張・曲げ・ねじりの3荷重モードでKtを瞬時に算出する。Petersonの線図をいちいち本で引く必要はもうない。
なぜ応力集中係数Kt計算ツールを作ったのか
開発のきっかけ
大学時代、材料力学の試験問題で「段付き軸のフィレット部の最大応力を求めよ」と出題された。教科書にはPetersonの線図が載っていたが、対数スケールの曲線から正確な値を読み取るのは至難の業。r/dやD/dの比率が線図のカーブの間にあると、内挿するしかなく、隣の人と答えが微妙にずれるという謎の体験をした。
社会人になってからも状況は変わらなかった。軸設計のたびにPilkeyの分厚い本を引っ張り出し、D/d比とr/d比のマトリクスから係数を探す作業を繰り返していた。日本語で使えるWebベースの応力集中係数計算ツールを探したが、見つからなかった。英語のツールはいくつかあったものの、形状図がなく、どのパラメータがどの寸法に対応するのか迷うことが多い。結局「自分で作るしかない」と決めた。
こだわった設計判断
- 6形状×3荷重モード対応: 段付き丸棒だけでなく、平板の穴やU溝など現場でよく出会う形状を網羅。荷重モードもねじりまでカバーした
- SVG形状図リアルタイム連動: 入力した寸法に合わせて断面図が変化し、最大応力点が赤丸で表示される。パラメータの対応関係が一目でわかる
- 有効範囲警告: Petersonの近似式はr/d比が特定範囲でしか精度が保証されない。範囲外に入ったらその旨を結果に付記するようにした
応力集中とは何か——力の流れが「渋滞」する現象
応力集中係数 とは
物体に外力がかかると、内部には応力という力の流れが発生する。均一な断面の棒であれば、応力は断面全体に均等に分布する。しかし、穴・溝・段差のような形状の急変部があると、力の流れが急に狭くなり、その部分に応力が集中する。
日常のたとえで言えば、高速道路の料金所を想像してみて。5車線の道路が2つのゲートに絞られると、ゲート手前で車が渋滞する。応力集中はまさにこれと同じ現象で、断面が急に細くなる場所では「力の渋滞」が起きて局所的に大きな応力が発生する。
応力集中係数Ktは、その渋滞の度合いを数値化したもの。定義は単純だ:
Kt = σmax / σnom
σmax: 形状急変部での最大応力
σnom: 公称応力(断面変化がないと仮定した場合の応力)
Kt=1.0なら応力集中なし。Kt=3.0なら公称応力の3倍の応力が局所的にかかっている。
穴あき板のKt=3.0——最も有名な応力集中
材料力学の教科書で最初に出会う応力集中は、無限板に小さな円孔がある場合のKt=3.0だ。Kirchhoffの弾性理論に基づく解析解で、穴の直径が板幅に比べて十分小さい場合に成り立つ。
直感的には、穴の上下では力が通れないため、左右の縁に力が迂回して集中する。ちょうど川の中に岩があると、岩の横の流速が速くなるのと同じだ。
段付き軸の応力集中——フィレットRが命
段付き軸(太い部分から細い部分に変わる軸)では、段差のフィレット半径rが応力集中を支配する。rが小さいほど力の流れの急変が激しく、Ktが大きくなる。逆にrを十分大きく取れば、力がスムーズに流れてKtは1.0に近づく。
機械設計では「フィレットRをケチるな」とよく言われるが、その理由がまさに応力集中の低減にある。ただしRを大きくすると、軸受の嵌合部や止め輪溝との干渉が起きるため、設計はトレードオフになる。
なぜKtが設計で重要か——疲労破壊の8割は応力集中部から始まる
静的強度と疲労強度の違い
静的荷重(一回きりの力)であれば、延性材料の場合は応力集中部が局所的に降伏しても、周囲が荷重を分担してくれるため、即座に破壊にはつながらない。しかし繰返し荷重がかかる場合は事情が異なる。
疲労破壊は応力が最も高い場所から微小なき裂が発生し、荷重サイクルごとに成長して最終的に破断に至る現象だ。疲労破壊事例の実に80%以上が、穴・溝・段差・溶接止端といった応力集中部から起きていると報告されている。
応力集中を無視した設計の末路
古典的な事故例としてよく引用されるのがコメット連続墜落事故(1954年)だ。初のジェット旅客機コメットの窓には角型の窓が採用されていたが、角部のKtが極めて大きく、与圧サイクルによる疲労き裂が進展して空中分解に至った。事故後、航空機の窓が丸型に変更されたのは、応力集中を低減するためだ。
Kt=3.0とKt=1.5では、疲労寿命がオーダー(10倍以上)で変わることもある。応力集中係数を正しく把握し、設計に反映することは、安全性とコストの両面で極めて重要だ。
Ktが活躍する設計シーン
軸設計 — 段付き軸のフィレットR決定
回転軸の段差部は軸受やギアの位置決めに必要だが、段差のフィレットRが疲労強度を左右する。Ktを計算して許容応力振幅を求め、疲労安全率を確認するのが標準的な設計フロー。
板金穴配置 — ボルト穴・配線穴の周辺応力
板金パネルに開けた穴の周囲はKt≈3.0の応力集中が発生する。穴径と板幅の比で正確なKtが変わるため、このツールで確認してから穴位置を決めるのが効率的。
FEA結果の妥当性検証
有限要素法(FEA)で得られた最大応力がPetersonのKtから算出される理論値と大きくずれている場合、メッシュ品質や境界条件に問題がある可能性がある。Ktによる概算値は「答え合わせ」として使える。
疲労設計 — Kf(切欠き感度係数)の出発点
疲労設計では、KtそのものではなくKf(切欠き係数)を使う。Kfは材料の切欠き感度qとKtから算出される:
Kf = 1 + q × (Kt - 1)
Ktがわからなければ始められない計算だ。疲労寿命シミュレーターと組み合わせて使うのがおすすめ。
基本の使い方
形状を選んで寸法を入れるだけの3ステップ。
Step 1: 形状と荷重を選ぶ
ドロップダウンから6形状のいずれかを選択し、荷重モード(引張/曲げ/ねじり)をタップする。形状によってはねじりが選択不可になる(平板系はねじりに対応しない)。
Step 2: 寸法を入力する
大径D・小径d・フィレット半径rをmm単位で入力する。穴あき平板の場合はD=板幅W、r=穴径d₀として入力する。入力に連動してSVG形状図がリアルタイムに更新される。
Step 3: Ktを確認する
応力集中係数Ktが即座に計算され、カード表示される。低い(青)〜高い(赤)の色分けで一目で程度がわかる。「結果をコピー」ボタンで計算条件と結果をクリップボードに保存できる。
具体的な使用例と検証データ
ケース1: 段付き丸棒 — フィレットR=5mmの軸(引張)
入力値:
- D = 50mm, d = 30mm, r = 5mm
- 荷重: 引張
計算結果:
- D/d = 1.667, r/d = 0.167
- Kt ≈ 1.67
→ 解釈: 軽度の応力集中。このフィレットRなら静的強度は十分だが、疲労設計では公称応力の1.67倍を考慮する必要がある。
ケース2: 段付き丸棒 — フィレットR=1mmの鋭い段差(曲げ)
入力値:
- D = 40mm, d = 30mm, r = 1mm
- 荷重: 曲げ
計算結果:
- D/d = 1.333, r/d = 0.033
- Kt ≈ 2.75
→ 解釈: 中程度〜高い応力集中。R=1mmは実務で見かける値だが、曲げ疲労を受ける軸としてはかなり厳しい。フィレットRの拡大(3mm以上推奨)を検討すべき。
ケース3: 穴あき平板 — 穴径10mm/板幅100mm(引張)
入力値:
- W = 100mm, d₀ = 10mm
- 荷重: 引張
計算結果:
- d₀/W = 0.1
- Kt ≈ 2.73
→ 解釈: 教科書的なKt≈3.0よりやや低い。穴径が板幅の1/10程度なら穴の影響はあるが致命的ではない。
ケース4: 穴あき平板 — 大穴 d₀=40mm/板幅80mm(引張)
入力値:
- W = 80mm, d₀ = 40mm
- 荷重: 引張
計算結果:
- d₀/W = 0.5
- Kt ≈ 2.31
→ 解釈: 穴が大きいほどKtは下がる(意外に感じるかもしれない)。ただし公称応力自体は穴で断面が減るため大幅に上がる。Ktだけでなく有効断面での応力も確認が必要。
ケース5: U溝付き丸棒 — 深いU溝(引張)
入力値:
- D = 30mm, d = 20mm, r = 3mm
- 荷重: 引張
計算結果:
- D/d = 1.5, r/d = 0.15
- Kt ≈ 1.83
→ 解釈: 軽度〜中程度。止め輪溝など軸のU溝はよくあるが、r=3mmの溝底Rがあればそこまで深刻ではない。
ケース6: 段付き丸棒 — ねじり荷重
入力値:
- D = 60mm, d = 40mm, r = 5mm
- 荷重: ねじり
計算結果:
- D/d = 1.5, r/d = 0.125
- Kt ≈ 1.44
→ 解釈: ねじり荷重のKtは同じ形状の引張・曲げより小さい傾向がある。動力伝達軸で段差がある場合、ねじり方向の応力集中は比較的穏やかだ。
仕組み・アルゴリズム——Petersonの近似多項式
3つの算出手法の比較
| 手法 | 精度 | 速度 | 入手性 |
|---|---|---|---|
| FEA(有限要素法) | ◎ 最高精度 | △ メッシュ作成に数十分 | △ 商用ソフト必要 |
| Peterson線図(目読み) | ○ 実用十分 | △ 手作業で10分 | ○ 書籍に掲載 |
| Peterson近似多項式 | ○ 線図同等 | ◎ 瞬時 | ◎ 本ツール |
本ツールはPeterson近似多項式を採用した。Pilkeyの「Peterson's Stress Concentration Factors(第4版)」に掲載されている多項式係数を使い、r/d比とD/d比からKtを算出する。
計算フロー
段付き丸棒・引張の場合:
1. 入力: D=50mm, d=30mm, r=5mm
2. 比率を計算:
D/d = 50/30 = 1.667
2r/D = 2×5/50 = 0.2
3. 多項式係数を算出:
√(D/d-1) = √0.667 = 0.8165
C1 = 0.926 + 1.157×0.8165 - 0.099×0.667 = 1.804
C2 = 0.012 - 3.036×0.8165 + 0.961×0.667 = -1.839
C3 = -0.302 + 3.977×0.8165 - 1.744×0.667 = 2.583
C4 = 0.365 - 2.098×0.8165 + 0.878×0.667 = -0.936
4. Kt = C1 + C2×t + C3×t² + C4×t³
= 1.804 + (-1.839)×0.2 + 2.583×0.04 + (-0.936)×0.008
= 1.804 - 0.368 + 0.103 - 0.007
≈ 1.53
(実際の計算はJavaScriptで浮動小数点演算するため微小な丸め誤差が生じる)
穴あき平板はHowlandの解析解
穴あき平板の引張時はPetersonの多項式ではなく、Howlandの解析解を採用している:
Kt = 3.0 - 3.13×(d₀/W) + 3.66×(d₀/W)² - 1.53×(d₀/W)³
d₀/W→0(穴が非常に小さい場合)でKt→3.0に収束し、理論解と一致する。
なぜ近似多項式を選んだか
FEAは精度最高だが、ソフトウェアとモデリングの手間が必要。線図の目読みは10分かかる上に個人差が出る。近似多項式なら瞬時に計算でき、Pilkeyの書籍で係数が公開されているため検証性も高い。「線図を引くまでもない概算」と「FEAの答え合わせ」の両方に適した手法だ。
既存ツール・手法との違い
日本語対応の専用Webツール
英語圏にはKt計算のWebツールがいくつか存在するが、日本語で使えるものはほぼ見つからない。本ツールは入力ラベル・結果・警告メッセージすべてが日本語で、JIS規格に慣れたエンジニアがそのまま使える。
形状SVG図による直感性
多くの競合ツールは数値入力と結果だけを表示する。本ツールはSVG形状図が入力に連動してリアルタイム更新され、最大応力発生点が赤丸で可視化される。「DとdとrがSVGのどの寸法に対応しているか」で迷うことがない。
6形状のワンストップ対応
段付き丸棒だけ、穴あき板だけ、というツールは多いが、6形状を1つのツールで切り替えて使えるのは珍しい。設計フェーズで複数の形状を比較検討する際に、ツールを切り替える手間がない。
応力集中にまつわる豆知識
サンブナンの原理——応力集中は局所的
サンブナンの原理によれば、局所的な荷重や形状変化による応力の乱れは、その特徴的寸法の数倍離れると消滅する。つまり、穴の直径が10mmなら、穴の縁から30〜50mm離れた場所では応力分布はほぼ均一に戻る。応力集中は「局所病」であり、全体の強度を支配するのは局所の最大応力だ。
切欠き感度——材料によってKtの効き方が違う
理論上のKtがいくら大きくても、実際に疲労強度がKt倍だけ低下するわけではない。材料の切欠き感度q(0〜1の係数)によって、実効的な応力集中係数Kfが決まる:
Kf = 1 + q(Kt - 1)
高強度鋼はq≈0.9〜1.0で切欠きに敏感。鋳鉄やアルミは q≈0.5〜0.8 で切欠きに鈍感。つまり同じKt=3.0でも、鋳鉄ならKf≈2.0〜2.6に収まることがある。切欠き感度の詳細はPetersonの教科書を参照してほしい。
Kt低減のための実践テクニック
フィレットRを最大限に取る
段付き軸のKtを下げる最もシンプルな方法は、フィレットRを大きくすること。R/d=0.1→0.2に変えるだけでKtが0.5〜1.0程度下がることも珍しくない。軸受の肩当て部では、軸受メーカーのカタログにフィレットR許容値が記載されているので確認してみて。
リリーフ溝で段差を緩和する
フィレットRを大きくできない場合(軸受の嵌合公差域が近いなど)、段差の手前にリリーフ溝を設けて力の流れを緩和する手法がある。溝自体にも応力集中は生じるが、段差直下よりKtが小さくなるケースが多い。
穴の周囲に補強リブを追加する
板金の穴周辺のKtを下げるには、穴の縁にフランジ加工(バーリング)やリブを追加する方法がある。穴径を変えずに実効的な応力集中を低減できるため、設計変更の影響が小さい。
よくある質問
Q: 動的荷重(衝撃荷重)の場合もKtをそのまま使える?
静的な応力集中係数Ktは形状のみで決まり、荷重の動的/静的は区別しない。ただし衝撃荷重では荷重値自体が静荷重の数倍になるため、Kt×衝撃荷重で最大応力を求める必要がある。衝撃係数は別途評価が必要だ。
Q: 疲労設計ではKtではなくKfを使うべき?
そのとおり。疲労設計では材料の切欠き感度を反映したKf(= 1 + q×(Kt-1))を使う。Ktは「形状だけで決まる幾何学的な係数」、Kfは「材料も考慮した実効的な係数」。高強度鋼(HRC>40)ではq≈1.0なのでKt≈Kfだが、軟鋼や鋳鉄ではKf < Ktになる。
Q: FEAの結果とKtの計算値が合わないのはなぜ?
よくある原因は3つ: (1) FEAのメッシュが粗く、応力集中部のピーク値が捉えられていない(メッシュ収束性の確認が必要)、(2) FEAモデルの形状がPetersonの想定と微妙に異なる(面取りの有無、テーパーの有無など)、(3) 非線形材料挙動(降伏後の応力再分配)をFEAが考慮している。線形弾性解析のKtとは比較可能だが、弾塑性解析の最大応力とは直接比較できない。
Q: 入力データがサーバーに送信されることはある?
すべての計算はブラウザ内のJavaScriptで完結しており、入力値は一切外部に送信されない。社内の設計データを扱う場合もプライバシーの心配なく使える。
まとめ
応力集中係数Ktは、形状急変部の最大応力を把握するための基本中の基本。6形状×3荷重モードをワンページで即算出できるこのツールを、日々の設計に役立ててほしい。
疲労寿命の評価が必要なら疲労寿命シミュレーター、軸径の選定なら軸径選定シミュレーター、ボルトまわりの強度が気になるならボルト強度・破断モード診断ツールもあわせてチェックしてみて。
不具合や要望があれば、X (@MahiroMemo)から気軽に教えて。