生体力学・骨応力解析 戻る
バイオメカニクス

生体力学・骨応力解析シミュレーター

骨種別・体重・歩行位相・モーメントアームを変えて圧縮応力と曲げ応力をリアルタイム計算。骨断面の応力分布を可視化し、骨折リスクと疲労寿命を評価しよう。

骨・荷重設定
骨の種類
歩行位相
体重 BW (kg)
kg
曲げモーメントアーム (cm)
cm
圧縮力 F
曲げモーメント M
圧縮応力 σc
曲げ応力 σb
最大応力 σmax
安全率
疲労寿命 Nf
安全
計算結果
圧縮力 (kN)
曲げ応力 (MPa)
最大応力 (MPa)
安全率
断面
応力
理論・主要公式
$$\sigma_c = \frac{F}{A}, \quad \sigma_b = \frac{M \cdot y}{I}$$

A = π(R²−r²),I = π(R⁴−r⁴)/4
骨圧縮強度 170 MPa,引張強度 130 MPa

生体力学・骨応力解析とは

🙋
「骨応力解析」って何ですか? 骨にどれくらい力がかかるか計算するんですか?
🎓
その通り!大まかに言うと、歩いたり走ったりした時に骨のどこにどれだけ負担がかかるかを、工学の力を使って計算・予測する分野だ。例えば、このシミュレーターで「歩行位相」を「立脚中期」に、「体重」を70kgに設定してみて。関節にかかる力が体重の約3.5倍、つまり約2400ニュートンと計算されるのがわかるよ。
🙋
え、そんなに大きい力がかかるんですか!で、その力が「応力」になる?「圧縮」と「曲げ」の違いは?
🎓
いいところに気づいたね。圧縮応力は骨を真っ直ぐ押しつぶす力。曲げ応力は骨を「てこの原理」で折り曲げようとする力だ。実務では、大腿骨の頚部などはこの曲げ応力が大きくかかる場所で、骨折のリスクが高くなる。上の「曲げモーメントアーム」のスライダーを動かすと、この曲げの影響がどう変わるか、応力分布の図で直感的にわかるよ。
🙋
なるほど!で、計算された応力が赤くなってるけど、あれが危険ということ?「疲労寿命」って何ですか?
🎓
その通り、赤い部分は局所的に高い応力が集中していることを示している。骨にも金属と同じように「疲労」があるんだ。何千回、何万回も繰り返し負荷がかかると、目に見えない微細なき裂が進展して、ある日突然疲労骨折を起こす。マラソンランナーに多いあの骨折だ。このツールでは、その繰り返し負荷で何回まで耐えられるか、「疲労寿命」も推定できるんだ。

よくある質問

本シミュレーターは大腿骨長幹部を想定した簡易モデルです。実際の骨は非対称で形状が複雑ですが、応力分布の傾向を把握するには十分な精度です。より詳細な解析が必要な場合は、患者個別のCTデータを用いた有限要素法をご検討ください。
歩行位相(立脚期・遊脚期など)によって骨にかかる力の方向と大きさが変化します。例えば立脚初期では体重の約3倍の荷重がかかり、曲げモーメントも最大になります。位相を変えることで、歩行サイクル中の応力変動をリアルタイムで確認できます。
モーメントアームが長くなると、同じ力でも曲げモーメントが比例して増大し、骨表面の曲げ応力が大きくなります。その結果、特に骨外側の引張応力が高まり、骨折リスクが上昇します。例えば、股関節外転筋のモーメントアームが短いとリスクが低減します。
各歩行位相での最大応力値を基に、S-N曲線(応力-寿命曲線)を用いて疲労寿命を推定しています。骨の材料特性として、圧縮側より引張側で疲労強度が低いことを考慮しています。あくまで相対的な指標であり、絶対的な寿命予測には注意が必要です。

実世界での応用

人工関節(インプラント)の設計:大腿骨頭を人工物に置き換える手術では、残った骨と人工物の接合部にかかる応力を詳細に解析します。骨の強度(圧縮170 MPa、引張130 MPa)に対して十分な安全率を確保する設計が必須で、CAEシミュレーションはその核心ツールです。

スポーツ医学と疲労骨折の予防:長距離ランナーやバスケットボール選手など、繰り返し高負荷がかかるアスリートの骨の疲労寿命を推定します。トレーニング強度やフォームの改善により、応力集中を緩和し、疲労骨折のリスクを低減する指導に役立ちます。

骨粗鬆症治療の効果評価:骨密度が低下すると、断面二次モーメント$I$や面積$A$が減少し、同じ負荷でも応力が増大します。薬物治療により骨の幾何学的形状(中空円筒の$R$や$r$)がどう変化し、機械的強度がどう改善されるかを定量的に評価する指標となります。

リハビリテーション工程の最適化:骨折後の患者が、いつ、どの程度の負荷をかけて歩行訓練を開始すべきかの判断材料を提供します。骨癒合部にかかる応力をシミュレーションで予測し、再骨折のリスクを最小化しながら、早期の社会復帰を目指します。

よくある誤解と注意点

まず、このシミュレーターの結果を「絶対的な診断」と捉えないでください。あくまで、特定の条件下での「傾向」や「比較」を見るためのツールです。例えば、骨強度のデフォルト値は平均的な値ですが、実際の骨は個人差が非常に大きい。骨粗鬆症の方とアスリートでは、同じ応力でも骨折リスクは全く異なります。シミュレーション結果は、その人の実際の骨密度や微細構造を反映していないことを常に念頭に置きましょう。

次に、「最大応力一点」だけを見て判断しないこと。応力が局所的に高くても、その周囲の骨が十分に強ければ問題ないケースもあります。逆に、全体として応力は高くなくても、特定の方向への繰り返し負荷(サイクリックロード)が「疲労寿命」を著しく短くする場合があります。例えば、歩行時の脛骨では、圧縮応力自体は大したことなくても、前後の曲げ応力の繰り返しが「シンスプリント」や疲労骨折の原因になります。結果は分布図と寿命推定を総合的に評価してください。

最後に、境界条件の設定には細心の注意を。例えば「関節にかかる力」を入力する際、これは単に体重の倍数で済む話ではありません。動作の位相(踵接地から蹴り出しまで)によって力の方向と大きさは刻一刻と変化します。また、骨と軟骨、靭帯との相互作用を単純な「固定」や「ピン」で近似している点も限界です。実務では、これらの設定が結果に与える影響を感度解析で確認することが鉄則です。「このパラメータを10%変えたら、結果は何%変わるか?」という視点を持ちましょう。

使い方ガイド

  1. 体重スライダー(bwSlider)を50~150 kgの範囲で調整し、対象患者の実体重を入力します
  2. 腕の負荷係数スライダー(armSlider)を0~2.0の範囲で設定し、片腕懸垂などの動作負荷を指定します
  3. 「解析実行」ボタンを押すと、大腿骨近位部の圧縮力(kN)、椎体L4の曲げ応力(MPa)、最大応力およご安全率(骨強度との比)がリアルタイム表示されます

具体的な計算例

体重80 kg、腕負荷係数1.5(階段昇降時の衝撃)の場合、大腿骨頸部に約2.8 kNの圧縮力が作用し、曲げ応力は皮質骨許容値90 MPaに対して約68 MPaとなります。最大応力が発生する部位はCalcar femorale(大腿骨内側支持部)で、安全率は1.32です。健康な成人骨(Young's modulus 17 GPa)であれば骨折リスク低値ですが、骨粗鬆症患者(骨密度T値-2.5以下)では安全率が0.95まで低下し、転倒時の病的骨折リスクが高まります

実務での注意点