構造解析の分野で欠かせない手法である有限要素法（FEM：Finite Element Method）において、剛性マトリクスは解析の核心をなす重要な概念です。
剛性マトリクスとは何か、どのように組み立てられ、構造解析においてどのような意味を持つのかを理解することは、FEMを使いこなすための基礎となります。
本記事では、剛性マトリクスの基本的な定義と物理的な意味から、有限要素法における組み立て方・節点変位と節点力の関係・計算方法・具体的な計算例まで、わかりやすく体系的に解説していきます。
FEMを学んでいる学生の方・構造解析に取り組むエンジニアの方・剛性マトリクスの理解を深めたい方にとって、役立つ内容となっているでしょう。
目次
剛性マトリクスとは？物理的な意味と定義
それではまず、剛性マトリクスの定義と物理的な意味を解説していきます。
剛性マトリクス（剛性マトリックス・stiffness matrix）とは、構造物の各節点における力と変位の関係を行列形式で表したものです。
簡単に言えば、「どれだけの変位を与えると、どれだけの力が発生するか」の関係を行列の形で表現したものです。
剛性マトリクスを用いることで、複雑な構造物における多数の節点の力と変位の関係を一つの行列方程式として整理でき、コンピュータによる数値解析を効率的に行うことが可能になります。
剛性マトリクスの基本方程式
[K]{u} = {f}
[K]：剛性マトリクス（n×n の正方行列、n は自由度数）
{u}：節点変位ベクトル（各節点の変位・回転角を並べたベクトル）
{f}：節点力ベクトル（各節点に作用する力・モーメントを並べたベクトル）
この方程式を解くことで、外力が既知の場合に各節点の変位を求めることができます。
剛性マトリクスの物理的意味
剛性マトリクスの各成分Kijには明確な物理的意味があります。
Kijは、第j自由度に単位変位を与え他のすべての自由度を固定した場合に、第i自由度に発生する反力を意味します。
これは剛性の逆数である柔性マトリクス（コンプライアンスマトリクス）と対比して理解すると分かりやすいでしょう。
剛性マトリクスは対称行列（Kij＝Kji）であり、これはマクスウェルの相反定理によって保証された性質です。
また、剛性マトリクスは正定値行列（すべての固有値が正）であり、これは物理的に意味のある変形（エネルギーが蓄積される変形）が必ず正の仕事を要することを反映しています。
1次元ばね要素の剛性マトリクス
剛性マトリクスの概念を最も簡単な形で理解するために、1次元ばね要素から始めましょう。
1次元ばね要素の剛性マトリクス
ばね定数kを持つばね要素（節点1・節点2）の場合
[k_e] = k × [ 1, −1; −1, 1 ]
節点力と節点変位の関係：
{ f₁ } k [ 1 −1 ] { u₁ }
{ f₂ } = k [−1 1 ] { u₂ }
k：ばね定数 [N/m]
u₁、u₂：節点1、節点2の変位
f₁、f₂：節点1、節点2の節点力
このばね要素の剛性マトリクスを理解することで、有限要素法における要素剛性マトリクスの概念の基礎が身につきます。
節点1を固定（u₁=0）して節点2に変位u₂を与えると、f₂=k×u₂ という直感的に理解できる結果が得られます。
全体剛性マトリクスの組み立て
有限要素法では、各要素の剛性マトリクス（要素剛性マトリクス）を組み合わせて全体剛性マトリクスを組み立てるアセンブリプロセスが重要です。
各要素の剛性マトリクスは要素の局所座標系で定義されており、必要に応じて全体座標系に変換してから全体剛性マトリクスの対応する位置に加算されます。
このプロセスは「直接剛性法」と呼ばれ、現代の有限要素法プログラムの基本的な計算手順です。
全体剛性マトリクスは各要素の剛性の重ね合わせとして組み立てられるため、疎行列（ほとんどの成分がゼロの行列）になります。
この疎行列性を活用した効率的な連立方程式の解法（スカイラインソルバー・スパースソルバーなど）により、大規模構造物の解析が可能になっています。
有限要素法における要素剛性マトリクス
続いては、有限要素法における各種要素の剛性マトリクスを確認していきます。
有限要素法では構造物をいくつかの要素（elements）に分割し、各要素の剛性マトリクスを組み合わせて全体の解析を行います。
トラス要素の剛性マトリクス
トラス（軸力のみを伝達する部材）の剛性マトリクスを求める方法を確認します。
2次元トラス要素では、各節点に水平方向（x）と鉛直方向（y）の2自由度があり、1要素に2節点×2自由度＝4×4の剛性マトリクスになります。
2次元トラス要素の剛性マトリクス
断面積A、ヤング率E、長さL、部材角度θのトラス要素の剛性マトリクス（全体座標系）
k_e = (AE/L) × [c², cs, −c², −cs; cs, s², −cs, −s²; −c², −cs, c², cs; −cs, −s², cs, s²]
c = cosθ（水平方向コサイン）
s = sinθ（鉛直方向コサイン）
自由度の順：[u₁, v₁, u₂, v₂]（節点1のx・y変位、節点2のx・y変位）
部材角度θ＝0（水平部材）の場合、非零成分は端点のx方向変位にのみ現れ、y方向の剛性項はゼロになります。
これはトラス部材が軸方向にのみ力を伝達するという物理的事実を正確に反映しています。
はり要素の剛性マトリクス
曲げモーメントと剪断力を伝達するはり要素の剛性マトリクスは、トラス要素より複雑になります。
オイラー-ベルヌーイのはり要素では、各節点に鉛直変位vと回転角θの2自由度があるため、2節点要素の剛性マトリクスは4×4行列です。
オイラー-ベルヌーイはり要素の剛性マトリクス
[k_e] = (EI/L³) × [12, 6L, −12, 6L; 6L, 4L², −6L, 2L²; −12, −6L, 12, −6L; 6L, 2L², −6L, 4L²]
EI：曲げ剛性
L：要素長さ
自由度の順：[v₁, θ₁, v₂, θ₂]（節点1のたわみ・回転角、節点2のたわみ・回転角）
このはり要素剛性マトリクスの各成分は、前述のたわみ公式（片持ちはり・単純支持はりのたわみ公式）と整合しており、曲げ剛性EIとスパン長Lが剛性の大きさを決定することが行列の形からも確認できます。
2次元・3次元ソリッド要素の剛性マトリクス
実際の構造解析では、より複雑な2次元・3次元ソリッド要素が多く使われます。
2次元平面応力・平面ひずみ問題に使われる三角形定ひずみ要素（CST要素）の場合、3節点×2自由度（x・y変位）＝6×6の剛性マトリクスになります。
要素剛性マトリクスは一般的に次の形式で計算されます。
ソリッド要素剛性マトリクスの一般式
[k_e] = ∫∫∫ [B]ᵀ [D] [B] dV
[B]：ひずみ-変位行列（形状関数の微分から計算）
[D]：弾性定数マトリクス（材料の応力-ひずみ関係を表す）
この積分は数値積分（ガウス求積法）によって計算されます。
[B]行列は要素内の変位場を記述する形状関数（補間関数）の微分から導出され、[D]行列はフックの法則に基づく弾性係数テンソルを行列形式にしたものです。
ガウス積分点の数と配置が解析精度に大きく影響するため、要素タイプと次数に応じた適切な積分則の選択が必要です。
全体剛性マトリクスの組み立てと境界条件の処理
続いては、要素剛性マトリクスから全体剛性マトリクスを組み立てる手順と、境界条件の処理方法を確認していきます。
アセンブリプロセス（直接剛性法）
全体剛性マトリクスの組み立ては、各要素の剛性マトリクスを全体の自由度番号に対応させて積み上げるアセンブリ（組み立て）プロセスによって行われます。
アセンブリプロセスの手順
1. 全体剛性マトリクス[K]を零行列で初期化する（サイズ：全自由度数×全自由度数）
2. 各要素について要素剛性マトリクス[k_e]を計算する
3. 要素の局所自由度と全体自由度の対応関係（自由度マッピング）を定義する
4. [k_e]の各成分を対応する全体剛性マトリクスの位置に加算する（重ね合わせ）
5. すべての要素について2〜4を繰り返す
このプロセスによって、各節点で隣接する要素の剛性が自動的に積み上げられた全体剛性マトリクスが完成します。
共有節点を持つ要素同士の剛性は全体マトリクス上で加算されるため、接続部分の連続性が自然に表現されます。
境界条件の処理方法
アセンブリ後の全体剛性マトリクスは、変位拘束（境界条件）が与えられていない状態では正則でなく（特異行列）、そのままでは連立方程式を解くことができません。
これは剛体変位（位置・回転が全体的にずれる運動）が拘束されていないためであり、境界条件の導入によって剛体変位自由度を拘束する必要があります。
境界条件の処理方法としては、消去法（拘束自由度の行・列を削除する方法）とペナルティ法（対角成分に大きな値を加算する方法）が一般的です。
消去法は理論的に明確ですが、自由度番号の再割り当てが必要です。
ペナルティ法は元の行列サイズを変えずに境界条件を処理できるため、プログラム実装が容易ですが、ペナルティ係数の大きさが数値精度に影響します。
連立方程式の求解
境界条件処理後の剛性方程式 [K]{u}={f} を解くことで節点変位ベクトル{u}が求まります。
大規模FEM解析では連立方程式のサイズが数十万〜数億自由度に達することもあり、効率的な求解アルゴリズムの選択が解析時間を大きく左右します。
直接法（LU分解・コレスキー分解など）は解の精度が高く信頼性があり、比較的小規模な問題に適しています。
反復法（共役勾配法・GMRES法など）は大規模問題で直接法より効率的になる場合が多く、前処理（プリコンディショニング）との組み合わせで収束性を改善することが一般的です。
剛性マトリクスが疎行列であることを活用した疎行列ソルバー（PARDISO・MUMPS・SuperLUなど）が広く利用されています。
剛性マトリクスの特性と数値解析における注意点
続いては、剛性マトリクスの数学的特性と数値解析における重要な注意点を確認していきます。
剛性マトリクスの対称性・正定値性
剛性マトリクスはいくつかの重要な数学的特性を持っており、これらを理解することがFEM解析の品質管理に役立ちます。
まず対称性（Kij＝Kji）については、線形弾性体でエネルギー保存則が成立する場合、レシプロカリティ（相反定理）から剛性マトリクスは必ず対称行列となります。
この性質により、行列の上三角部分のみ計算・格納すればよいという計算・メモリ効率上の利点があります。
正定値性（すべての固有値が正）は、物理的に意味のある弾性変形に対してエネルギーが正であることを保証します。
剛体変位が残っている場合や材料特性の設定に問題がある場合には正定値性が失われ、解析が正常に収束しないことがあります。
要素の選択とロッキング問題
FEM解析の精度は使用する要素の種類・形状・数に大きく依存します。
特に注意すべき問題としてロッキング（locking）現象があります。
体積ロッキングは非圧縮性に近い材料（ゴム・塑性変形状態の金属など、ポアソン比が0.5に近い材料）を低次要素で解析する際に発生し、変位が実際より小さく計算される誤差です。
剪断ロッキングは薄い板・シェル・はりを解析する際に発生し、不必要な剪断変形エネルギーが蓄積されることで曲げ変形が過小評価される現象です。
これらのロッキング問題を回避するために、低減積分（Reduced Integration）・非適合要素・拡張有限要素法（XFEM）・高次要素の使用などの手法が開発されています。
メッシュ依存性と解析精度の検証
有限要素法の解析結果はメッシュの分割数（メッシュ密度）に依存するため、適切なメッシュ設定が解析精度を確保するうえで重要です。
一般的に要素を細かく分割するほど厳密解に近い結果が得られますが、計算コストも増加します。
解析精度の検証方法として「メッシュ収束性確認」があり、メッシュを順次細かくして解析結果が収束しているかを確認することで、メッシュ設定の適切さを評価します。
応力集中部・き裂先端・接触境界などの応力勾配が大きい領域では局所的にメッシュを細かくする適応的なメッシュ分割（アダプティブメッシュリファインメント）が有効です。
また、解析結果の妥当性は手計算による概算値との比較・対称性条件の確認・反力の釣り合い確認・実験値との比較などで検証することが重要です。
剛性マトリクスの応用：動的解析と非線形解析
続いては、剛性マトリクスを応用した動的解析と非線形解析について確認していきます。
動的解析における質量マトリクスと剛性マトリクス
動的問題（振動・衝撃・地震応答など）の解析では、剛性マトリクスに加えて質量マトリクス[M]と減衰マトリクス[C]が必要になります。
動的解析の基本方程式（運動方程式）
[M]{ü} + [C]{u̇} + [K]{u} = {f(t)}
[M]：質量マトリクス
[C]：減衰マトリクス
[K]：剛性マトリクス
{ü}：節点加速度ベクトル
{u̇}：節点速度ベクトル
{u}：節点変位ベクトル
{f(t)}：時刻tにおける外力ベクトル
自由振動解析（固有値解析）では外力を零として [K]{u} = ω²[M]{u} の固有値問題を解くことで固有振動数ωと振動モードを求めます。
固有振動数はEIが大きいほど（剛性が高いほど）高くなり、質量が大きいほど低くなるという関係があります。
地震応答解析・フラッター解析・回転機械の振動解析などで固有値解析は広く活用されています。
非線形解析における接線剛性マトリクス
大変形・材料非線形（塑性・クリープ・超弾性）・接触問題などの非線形解析では、剛性マトリクスが変形状態によって変化する接線剛性マトリクスを用いた反復計算が必要です。
ニュートン-ラフソン法は非線形FEM解析で最もよく使われる反復解法であり、各反復ステップで接線剛性マトリクスを更新しながら解に収束させます。
幾何学的非線形（大変形問題）では変形後の形状に基づく更新ラグランジュ法・全ラグランジュ法などの定式化が用いられ、接線剛性マトリクスには弾性剛性成分に加えて初期応力による幾何剛性マトリクス（応力剛性マトリクス）が加わります。
幾何剛性マトリクスは、引張応力によって剛性が増加し（張力場効果）、圧縮応力によって剛性が低下する（座屈）という現象を定量的に表現します。
FEMソフトウェアと剛性マトリクスの実装
現代の構造解析では、ANSYS・ABAQUS・LS-DYNA・MSC Nastran・Calculix・OpenFOAMなどの商用・オープンソースFEMソフトウェアが広く利用されています。
これらのソフトウェアでは剛性マトリクスの組み立て・境界条件処理・連立方程式の求解が自動的に行われるため、ユーザーはモデル定義・解析条件設定・結果評価に集中できます。
ただし、ソフトウェアが自動的に計算しているプロセスの背後にある剛性マトリクスの理論を理解していることが、適切な解析モデルの構築・エラー診断・結果の妥当性評価において非常に重要です。
近年は機械学習とFEMを組み合わせたデータ駆動型の構造解析も研究されており、大規模問題の高速化・精度向上に向けた新たなアプローチとして注目されています。
まとめ
本記事では、剛性マトリクスとは何か、有限要素法における意味・計算方法・組み立て手順・境界条件の処理・動的解析や非線形解析への応用まで幅広く解説してきました。
[K]は[K]{u}={f}という基本方程式の核心であり、各要素の要素剛性マトリクスをアセンブリすることで全体構造の力と変位の関係を行列形式で表現します。
1次元ばね要素・トラス要素・はり要素・ソリッド要素など、要素の種類によって剛性マトリクスの形式は異なりますが、[B]ᵀ[D][B]の積分形式という共通の枠組みで導出されます。
メッシュ設定・要素選択・境界条件の処理・求解アルゴリズムなど、実際のFEM解析では理論的な理解と実務的な知識の両方が求められます。
剛性マトリクスの理論を深く理解することが、高品質なFEM解析・適切な結果評価・創造的な構造設計につながるでしょう。
本記事が構造解析・材料力学・有限要素法の学習・実務に役立てば幸いです。
