剛性マトリクスとは？有限要素法での役割をわかりやすく解説！（stiffness matrix：節点：変位：荷重：構造解析など）
橋・ビル・機械部品・航空機など現代の構造物の設計・解析において欠かせないコンピューター計算手法が「有限要素法（FEM：Finite Element Method）」です。
その有限要素法の中核をなす数学的構造が「剛性マトリクス（Stiffness Matrix）」です。
剛性マトリクスとは何か、どのように構築され、構造解析においてどのような役割を果たすのかを理解することは、構造力学・CAE（コンピューター支援エンジニアリング）を使いこなす上で不可欠な基礎知識です。
本記事では、剛性マトリクスの定義・物理的意味・有限要素法での構築方法・荷重と変位の関係・全体剛性マトリクスの組み立てまで、わかりやすく体系的に解説していきます。
剛性マトリクスとは？その本質をひとことで言えば
それではまず、剛性マトリクスの定義とその物理的な意味について解説していきます。
剛性マトリクス（Stiffness Matrix）とは、構造物の節点（nodes）に作用する力（荷重ベクトル）と節点変位（displacement vector）の関係を記述する正方行列です。
英語では「stiffness matrix（スティフネスマトリクス）」または「K matrix（Kマトリクス）」とも呼ばれます。
剛性マトリクスの基本方程式
剛性マトリクスの基本方程式
{F} = [K]{u}
{F}：節点荷重ベクトル（各節点に作用する外力・モーメントを並べたベクトル）
[K]：剛性マトリクス（n×n の正方行列、nは節点の自由度の総数）
{u}：節点変位ベクトル（各節点の変位・回転角を並べたベクトル）
既知の荷重{F}から変位{u}を求める場合：
{u} = [K]⁻¹{F}（剛性マトリクスの逆行列を計算）
この方程式は「スプリング（ばね）のフックの法則 F = k・u（F：力、k：ばね定数、u：変形量）」の多自由度版として理解できます。
剛性マトリクス[K]の各要素K_ijは「j番目の節点自由度に単位変位を与えたとき、i番目の節点自由度に生じる反力」を意味します。
剛性マトリクスの物理的意味
剛性マトリクスの各要素K_ijの物理的意味をより具体的に説明します。
K_ijとは「節点jの自由度iに単位変位を与え、他の節点は変位ゼロに固定したとき、節点iの自由度iに現れる反力の大きさ」です。
この定義から、剛性マトリクスの対角成分K_iiは正の値（自分自身への剛性）、非対角成分K_ijは結合剛性（他の節点との連成の強さ）を表すことがわかります。
剛性マトリクスの数学的性質
剛性マトリクスには以下の重要な数学的性質があります。
| 性質 | 内容 | 意味 |
|---|---|---|
| 対称性 | [K] = [K]^T（転置行列と等しい） | 相反定理（Maxwell-Betti）から導かれる |
| 正値半正定値性 | 任意ベクトル{u}に対して {u}^T[K]{u} ≥ 0 | 変形エネルギーが非負であることを保証 |
| 疎行列（スパース） | 多くの要素がゼロ | 直接接続されない節点間はK_ij = 0 |
| 境界条件適用前の特異性 | 拘束なしでは行列式 = 0 | 剛体移動が可能な状態（拘束なし）を示す |
有限要素法における剛性マトリクスの構築
続いては、有限要素法（FEM）において剛性マトリクスがどのように構築されるかについて確認していきます。
FEMの本質はこの剛性マトリクスの構築と求解プロセスにあります。
有限要素法の基本的な考え方
有限要素法とは、複雑な形状の連続体（梁・板・三次元固体など）を「有限な数の小さな要素（Elements）」に分割し、各要素の挙動を記述する剛性マトリクスを組み合わせることで全体の挙動を近似する数値解析手法です。
この分割操作を「メッシュ（Mesh）分割」と呼び、FEMの前処理（プリプロセス）における重要な作業です。
メッシュを細かくするほど解析精度は向上しますが、計算コスト（時間・メモリ）も増大するというトレードオフが存在します。
要素剛性マトリクスの導出
有限要素法では、まず各要素（Element）について「要素剛性マトリクス（Element Stiffness Matrix）」を導出します。
1次元の棒要素（軸力のみを伝達する要素）の要素剛性マトリクスを例に示します。
棒要素の要素剛性マトリクス（最も単純な例）
断面積A・ヤング率E・長さLの棒要素（節点1・節点2を持つ）の場合：
[k_e] = (AE/L) × [ 1 −1 ]
[−1 1 ]
節点変位ベクトル：{u_e} = {u₁, u₂}^T
節点荷重ベクトル：{f_e} = {f₁, f₂}^T
{f_e} = [k_e]{u_e}
AE/Lは棒のばね定数（剛性）に相当する。
より複雑な要素（梁要素・三角形板要素・四面体要素・六面体要素など）では、同様の原理でより大きな要素剛性マトリクスが導出されます。
全体剛性マトリクスのアセンブリ
各要素の要素剛性マトリクスを組み合わせて、構造全体の「全体剛性マトリクス（Global Stiffness Matrix）」を構築する操作を「アセンブリ（Assembly）」と呼びます。
アセンブリの基本原則は「共有節点での変位の適合条件（各要素が共有節点で同じ変位を持つ）と力の平衡条件（節点に作用する全要素の力の和が外力と等しい）」です。
全体剛性マトリクスは通常、対角成分付近にしか非ゼロ要素を持たない疎行列（スパース行列）となり、この特性を利用した効率的な連立方程式の求解アルゴリズムが有限要素ソフトウェアの性能を決定します。
境界条件の適用と連立方程式の求解
続いては、剛性マトリクス方程式の求解に必要な境界条件の適用と連立方程式の解法について確認していきます。
正しい境界条件の設定が有限要素解析の精度を左右します。
境界条件の種類と適用方法
有限要素法での境界条件には主に2種類あります。
強制変位境界条件（Dirichlet境界条件）は、特定の節点の変位・回転角をゼロまたは指定値に固定する条件です（例：固定端で変位＝0・回転角＝0）。
荷重境界条件（Neumann境界条件）は、特定の節点に作用する外力・外モーメントを指定する条件です（例：自由端での集中荷重）。
境界条件を適用することで、全体剛性マトリクスから拘束自由度に対応する行・列を除去（または修正）し、正則な連立方程式として解けるようになります。
連立方程式の求解手法
境界条件適用後の剛性マトリクス方程式{F} = [K]{u}は、線形代数の連立方程式として求解されます。
主要な求解手法として、直接法（ガウスの消去法・コレスキー分解）と反復法（共役勾配法・前処理付き共役勾配法）があります。
大規模なFEMモデル（節点数が数百万以上）では反復法が効率的であり、疎行列の性質を活用した高速アルゴリズムが商用FEMソフトウェアに実装されています。
剛性マトリクスの応用範囲と実務での利用
剛性マトリクスの概念は構造解析だけでなく、熱伝導解析・流体解析・電磁場解析など多くの物理現象のFEM解析に拡張して適用されています。
ANSYS・NASTRAN・ABAQUS・COMSOL・Autodesk Simulationなどの商用CAEソフトウェアは、すべて剛性マトリクス方程式を核心的な計算エンジンとして動作しています。
剛性マトリクスが活用される主要な解析分野
・線形静的構造解析：変位・応力・ひずみの計算（最も基本的な用途）
・固有値解析（振動解析）：固有振動数・振動モード形の計算
・座屈解析：圧縮荷重による座屈荷重・座屈モードの計算
・熱応力解析：温度変化による熱応力・熱変形の計算
・動的解析：時刻歴応答・周波数応答の計算（質量マトリクスと減衰マトリクスと組み合わせ）
まとめ
本記事では、剛性マトリクスの定義・物理的意味・数学的性質・有限要素法における要素剛性マトリクスの導出・全体剛性マトリクスのアセンブリ・境界条件の適用・求解手法・工学応用まで幅広く解説してきました。
剛性マトリクスとは節点荷重と節点変位の関係{F} = [K]{u}を表す正方行列であり、有限要素法という現代工学の最重要数値解析手法の数学的核心です。
対称性・スパース性・正値半正定値性という数学的性質が、FEMの効率的な計算アルゴリズムの基盤となっています。
構造設計・機械設計・建築・土木・航空宇宙など現代のあらゆる工学分野において、剛性マトリクスの概念はCAE解析の根底として機能し続けており、工学を学ぶ上で理解すべき最重要基礎知識のひとつといえるでしょう。
