集中質量要素
理論と物理
集中質量要素とは
先生、「集中質量要素」って何ですか？ 静解析で質量が必要になることはあるんですか？
集中質量要素は1つの節点に質量を付加する要素だ。静解析では自重（重力荷重）の計算に使うが、主な用途は動的解析（固有振動、時刻歴応答、衝撃）だ。
質量マトリクス
集中質量要素の質量マトリクスは対角行列：
$$ [M] = \begin{bmatrix} m & 0 & 0 & 0 & 0 & 0 \\ 0 & m & 0 & 0 & 0 & 0 \\ 0 & 0 & m & 0 & 0 & 0 \\ 0 & 0 & 0 & I_{xx} & 0 & 0 \\ 0 & 0 & 0 & 0 & I_{yy} & 0 \\ 0 & 0 & 0 & 0 & 0 & I_{zz} \end{bmatrix} $$
集中質量要素の質量マトリクスは対角行列：
ここで $m$ は並進質量、$I_{xx}, I_{yy}, I_{zz}$ は回転慣性モーメント。
回転慣性モーメントも設定できるんですね。
モーターや大型機器のような回転部品を質量要素で代替する場合、回転慣性がないと固有振動モードが不正確になる。並進質量だけでなく回転慣性も正しく設定することが重要だ。
用途
集中質量要素の主な用途：
| 用途 | 説明 |
|---|---|
| 非構造質量の付加 | 配管内の流体、装置の重量 |
| 質量のバランス調整 | FEMモデルの全体質量を設計値に合わせる |
| 装置のモデル化 | モーター、バルブ等を質量点で代替 |
| 重心の調整 | 重心位置を実構造に合わせる |
| 質量スケーリング（陽解法） | 安定時間増分を大きくする |
配管内の流体の質量はどう入れるんですか？
配管要素の各節点に流体の単位長さあたりの質量を集中質量として追加する。AbaqusではNONSTRUCTURAL MASSが便利で、要素に面密度（kg/m²）や線密度（kg/m）を追加できる。
ソルバー別の要素名
| ソルバー | 要素名 | 備考 |
|---|---|---|
| Nastran | CONM2 | 6DOFの集中質量。オフセット可能 |
| Abaqus | MASS / ROTARY INERTIA | 並進と回転を別々に定義 |
| Ansys | MASS21 | KEYOPT(3)で回転慣性の有無を選択 |
NastranのCONM2は「オフセット可能」？
CONM2は質量の重心を節点からオフセットできる。これにより重心が節点の位置になくても正しい慣性効果を表現できる。大型機器の取付点と重心がずれている場合に有用だ。
まとめ
集中質量要素を整理します。
要点：
- 1節点に質量を付加する要素 — 剛性は持たない
- 動的解析で主に使用 — 固有振動、時刻歴応答
- 並進質量と回転慣性 — 両方設定すべき（回転慣性を忘れない）
- 非構造質量の付加 — 流体、装置、付帯物の質量
- CONM2（Nastran）のオフセット — 重心位置の精密なモデル化
静解析では自重にしか関係しないけど、動的解析では固有振動数を支配するから超重要ですね。
その通り。構造の質量分布が間違っていたら固有振動数も間違う。質量要素の設定は動的解析の基本中の基本だ。
集中質量の理論的背景
集中質量行列（Lumped Mass Matrix）は一貫性質量行列（Consistent Mass Matrix）に対し、要素の質量を節点に集中配分する近似法だ。1968年にHintonらは集中質量行列の方が数値減衰が少なく陽的時間積分に適することを示した。Abaqus Explicitは全要素でこの方式を採用しており、現代の爆発・衝突解析の計算効率を支えている。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
集中質量と分布質量
FEMでは「集中質量」と「分布質量（一貫質量）」がありますよね。何が違いますか？
通常のFEM要素（梁、シェル、ソリッド）は材料密度 $\rho$ から質量マトリクスを自動生成する。これが一貫質量マトリクス（consistent mass matrix）。対して集中質量要素は追加の質量を特定の節点に付与する。
さらに、一貫質量マトリクスを対角化した集中質量マトリクス（lumped mass matrix）がある。
| 質量マトリクス | 特徴 | 用途 |
|---|---|---|
| 一貫質量 | 非対角成分あり。精度高い | 陰解法の動的解析 |
| 対角集中質量 | 対角のみ。計算効率高い | 陽解法（Explicit） |
| 追加集中質量（CONM2等） | ユーザーが追加する質量 | 非構造質量の付加 |
陽解法では対角集中質量が必須なんですか？
そう。陽解法の時間積分（中央差分法）は質量マトリクスの逆行列を使うから、対角でないと効率的に計算できない。LS-DYNAやAbaqus/Explicitでは質量マトリクスを自動的に対角化する。
質量の確認方法
モデルの全体質量が正しいかどうかはどう確認しますか？
ソルバーの質量サマリー出力を確認する：
確認項目：
- 全体質量が設計値（図面値、実測値）と一致するか
- 重心位置（$X_{CG}, Y_{CG}, Z_{CG}$）が妥当か
- 各方向の慣性モーメントが妥当か
重心位置の確認は見落としがちですね。
全体質量が合っていても重心がずれていたら、動的応答（特に転倒モーメント、偏心効果）が間違う。質量サマリーは解析の最初に必ず確認すべき。
質量スケーリング
「質量スケーリング」って何ですか？
陽解法の安定時間増分は $\Delta t \propto L / c$（$L$: 最小要素サイズ、$c$: 音速）で決まる。小さい要素があると $\Delta t$ が極端に小さくなる。質量を人為的に増やして $c$ を下げ、$\Delta t$ を大きくするのが質量スケーリング。
質量を増やすと慣性効果が変わりませんか？
変わる。だから質量スケーリングは準静的問題にのみ適用すべきで、追加された質量（kinetic energy）が全エネルギーの5%以下であることを確認する必要がある。
まとめ
集中質量の数値手法、整理します。
要点：
- 一貫質量 vs. 集中質量 — 陰解法は一貫質量、陽解法は集中質量
- 追加集中質量 — CONM2, *MASS で非構造質量を付加
- 質量サマリーの確認 — 全体質量、重心、慣性モーメント
- 質量スケーリング — 陽解法の時間増分を大きくするが慣性効果に注意
RBAR/RBEと集中質量の組合せ
試験で計測困難な付属物の質量をFEMに取り込む場合、剛体要素（RBE2/RBE3）と集中質量要素（CONM2）の組合せが実務標準だ。NastranのCONM2要素は質量と慣性モーメント6成分を指定でき、1975年の初期MSC Nastranから存在する歴史的要素だ。自動車ドアの開閉ハンドル・ヒンジの質量モデリングなどに現役で活躍する。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
集中質量の実務適用
集中質量要素の実務での使い方を教えてください。
最も一般的な適用：
機器のモデル化
工場やプラントの構造解析で、配管上の大型バルブやポンプを詳細にモデル化するのは非現実的。集中質量要素で質量と重心位置だけ表現するのが実務的だ。
手順：
1. 機器の質量 $m$、重心位置 $(x_g, y_g, z_g)$、慣性モーメント $(I_{xx}, I_{yy}, I_{zz})$ を図面から取得
2. 取付点にCONM2/MASSを配置
3. 必要に応じてRBE2/RBE3で取付点を分配
自動車のモデル化
NVH解析（振動騒音）では、エンジン、バッテリー、乗員をそれぞれ集中質量で表現することが多い。ボディのシェルモデルに質量要素を追加して全体の固有振動を評価する。
建築の付帯物
建築の耐震解析では、設備機器、家具、仕上げ材の質量を各階の床に分配する。床スラブの節点に集中質量を追加するか、面密度（NONSTRUCTURAL MASS）を設定する。
質量分布の妥当性確認
質量分布が正しいかどうか、どう確認しますか？
固有振動解析（モード解析）で確認するのが最良だ。1次モードの固有振動数が予想と大きくずれていたら、質量分布に問題がある。
簡易チェック：
- 片持ち梁 — $f_1 = 3.516/(2\pi) \sqrt{EI/(\rho A L^4)}$
- 単純梁 — $f_1 = \pi^2/(2\pi) \sqrt{EI/(\rho A L^4)}$
- 先端に集中質量 $M$ の片持ち梁 — $f_1 \approx 1/(2\pi) \sqrt{3EI/(M L^3)}$
実務チェックリスト
集中質量要素のチェックリストをお願いします。
回転慣性の設定忘れが一番多いミスですか？
そう。並進質量だけ設定して回転慣性を忘れると、ねじりモードや回転モードの振動数が過大になる。大型機器のモデル化では回転慣性が振動特性に大きく影響する。
人工衛星の質量モデリング
JAXA「はやぶさ2」（2014年打上）のFEMモデルでは、イオンエンジン・太陽電池パドル展開機構など個別質量をCONM2要素で表現し、総質量600kgのモデルを約15万自由度に抑えた。これにより振動試験との固有値一致精度±3%以内を実現し、打上荷重解析の信頼性を確保した。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
集中質量のツール比較
各ソルバーの集中質量要素の違いは？
NastranのCONM2が最も高機能ですか？
CONM2は1枚のカードで並進質量、回転慣性、オフセットを全て定義できる。構造解析の長い歴史で磨かれた要素だ。AbaqusではMASSとROTARY INERTIAを別々に定義する必要があり、やや手間がかかる。
選定ガイド
質量要素は基本的な要素だから、どのソルバーでも十分ですね。
そう。質量要素で差がつくのはソルバーの機能ではなく、エンジニアが質量分布を正しく把握しているかどうかだ。
各ソルバーの質量要素仕様比較
NastranのCONM2は6×6慣性マトリクス全成分を指定可能。AbaqusのPOINT MASS/ROTARYINERTIAは別々に定義する仕様で、慣性テンソルの指定方法が異なる。ANSYSのMASSE21は21番の要素タイプとして1980年代から存在し、KEY OPTION 1でRotary Inertiaのon/offを切り替える。ソルバー間変換時にこの差異でモード形状が変わる事例が知られる。
選定で最も重要な3つの問い
- 「何を解くか」：集中質量要素に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
質量要素の先端トピック
集中質量要素に先端研究はありますか？
質量要素自体はシンプルだが、質量の最適配置やモデル更新が活発だ。
モデル更新（Model Updating）
実験モード解析（振動試験）の結果とFEMの固有振動数・モード形状を一致させるために、集中質量や剛性を調整する手法。
手順：
実験モード解析（振動試験）の結果とFEMの固有振動数・モード形状を一致させるために、集中質量や剛性を調整する手法。
手順：
1. FEMで固有振動解析
2. 実験で固有振動数とモード形状を測定
3. FEMと実験の差が最小になるよう質量・剛性を最適化
4. 更新されたFEMで予測解析
振動試験の結果をFEMにフィードバックするんですね。
自動車のNVH開発では標準的なワークフローだ。試作車の振動試験結果でFEMを更新し、更新されたFEMで次の設計変更の効果を予測する。
トポロジー最適化と質量
構造のトポロジー最適化では質量制約が一般的：「全体質量が指定値以下で、剛性を最大化する」。集中質量要素は最適化の質量計算に含まれるため、正しく設定されていないと最適化結果が狂う。
デジタルツインと質量モデル
デジタルツインではリアルタイムのセンサーデータで質量分布を更新する。例えば橋梁のデジタルツインで、通行車両の質量と位置をリアルタイムで更新し、構造の応答を予測する。集中質量要素の動的な更新がこの基盤技術になる。
まとめ
集中質量の先端トピック、まとめます。
質量要素は「シンプルだが正確な設定が求められる」要素だ。設定の正確さが動的解析の全てを支配する。
非対角質量行列と結合慣性
集中質量要素に慣性テンソルの非対角成分（積慣性）を設定すると、FEM質量行列の非対角項を通じて振動モードが結合する。1990年代の工作機械主軸の解析では、この積慣性の無視が固有振動数を最大8%誤らせることをOkumaらが実証。以後、精密機械解析ではCONM2の慣性テンソル完全指定が標準化された。
トラブルシューティング
集中質量のトラブル
集中質量要素でよくあるトラブルを教えてください。
質量要素のトラブルは設定忘れか値の間違いがほとんどだ。
固有振動数が予想と大きく違う
固有振動数がFEMと手計算で全然合いません。
確認項目：
1. 質量が正しく設定されているか — 質量サマリーで全体質量を確認
2. 回転慣性を忘れていないか — ねじり/回転モードの振動数に影響
3. 質量の位置が正しいか — CONM2のオフセット、またはRBE2/RBE3での分配
4. 単位系 — kg? ton? 質量と力の単位が整合しているか（mm系ではkg→ton）
5. 材料密度 — 自動計算される質量と追加質量の重複がないか
単位系の問題は質量で特に多いですか？
mm-N-MPa系では質量の単位がton（= 1000 kg）になる。kg で入力すると1000倍の質量になり、固有振動数が $\sqrt{1000} \approx 32$ 倍ずれる。これは非常によくある間違いだ。
静解析で質量要素の影響がない
集中質量を追加したのに静解析の結果が変わりません。
静解析では質量は重力荷重（自重）にしか影響しない。重力荷重を設定していなければ質量要素は結果に影響しない。
確認：
- 重力加速度を荷重として入力しているか（*DLOAD, GRAV / GRAV カード / ACEL）
- 重力の方向と大きさが正しいか
質量の二重計上
質量が二重にカウントされることはありますか？
ある。材料密度から自動計算される質量と、追加の集中質量の両方がカウントされる。例えば配管の壁面はシェル要素で密度が定義されており、さらに配管内の流体質量をCONM2で追加する場合は正しい。しかしモーター全体をソリッドでモデル化しつつ、同じ質量をCONM2でも追加すると二重計上になる。
質量サマリーで全体質量を確認し、設計値と一致するか必ずチェックすること。
まとめ
集中質量のトラブル対処、整理します。
「質量サマリーを確認する」がデバッグの第一歩ですね。
質量サマリーはFEMモデルの「健康診断」だ。解析を回す前に必ず確認する習慣をつけよう。
質量要素の誤配置による振動異常
集中質量要素を構造の慣性重心から離れた節点にオフセットなしで配置すると、慣性力による偽モーメントが生じて固有振動数が実験値と乖離する。航空機の機内装備品モデリングで2008年にBoeing 787の設計解析チームがこの問題を特定し、RBE3+CONM2のオフセット付きモデリング手順を社内ガイドラインに追加した。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——集中質量要素の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告