リナンバリング — CAE用語解説

それは、有限要素法で解く連立一次方程式の係数マトリックス、つまり剛性マトリックスの「帯域幅」に直接影響するからです。節点番号がランダムだと、非ゼロ要素がマトリックス内で散らばり、帯域幅が広がります。例えば、100万自由度の問題で帯域幅が10000から500に縮小されると、必要なメモリ量は単純計算で20分の1近くになり、ソルバーの計算時間も劇的に短縮されます。

その通りです。直接法ソルバー（例えばスカイライン法やLDL^T分解）では、帯域内のすべての要素をメモリに保持し、演算します。帯域幅

基本的な考え方は、「隣接する節点同士の番号差をできるだけ小さくする」ことです。例えば、細長い梁を1次元要素でメッシュ分割する場合、端から順番に1, 2, 3...と番号を振るのが最適です。これに対して、ランダムに番号を振ると、隣接する節点1と1000が剛性マトリックス上では遠く離れた位置に配置され、帯域幅が最大節点差に近い値になってしまいます。

いくつかの古典的かつ強力なアルゴリズムがあります。代表的なのは「Cuthill-McKee (CM) アルゴリズム」とその逆順版「Reverse Cuthill-McKee (RCM)」です。これは、グラフ理論を用いて、次数（接続する辺の数）の小さい節点から探索を始め、隣接節点を番号付けしていく方法です。多くの商用ソフトウェアのデフォルトはこのRCMです。

はい。「Gibbs-Poole-Stockmeyer (GPS) アルゴリズム」や「Sloanのアルゴリズム」などがあります。これらはRCMよりも計算コストは高いものの、一般により狭い帯域幅を生成することが知られています。例えば、ある航空機エンジンブラケットのモデルで、RCMが帯域幅を1200から350に縮小したのに対し、Sloanのアルゴリズムは280まで縮小できた、という事例があります。

上位のソルバーや前処理モジュールでは可能です。例えば、MSC Nastranでは`PARAM,GRDPNT`や`PARAM,PRTMAXIM`などのカードでリナンバリングの出力レベルを制御し、内部で適切なアルゴリズムを選択します。オープンソースのCalculiXでは、`ccx`の実行時に`-i`オプションで入力ファイルを渡すと自動でRCMが適用されますが、ユーザーがアルゴリズムを直接指定するインターフェースはあまり一般的ではありません。

ほとんどのプリポストソフト（Ansys Mechanical, Abaqus/CAE）では、ソルバーが入力ファイル（.dat, .inp）を読み込んだ直後に自動で実行します。ユーザーが確認すべきは、ソルバーのログファイル(.log, .msg)です。例えばAbaqus/Standardの.msgファイルには「BEFORE RENUMBERING」と「AFTER RENUMBERING」のセクションがあり、行列のプロファイル（半帯域幅）や項数が記載されています。項数が大幅に減っていれば、リナンバリングが効いています。

よくある問題です。例えば、CATIAで作成したメッシュをNASTRAN形式でエクスポートすると、CATIA内部の幾何学的な作成順序に依存した変な番号付けになることがあります。これを「as-is」でソルバーに渡すと効率が悪い。対策は二つ：1) プリプロセッサ（FemapやAnsys Meshing）でメッシュを再生成する。2) ソルバーの入力デッキにリナンバリングを強制するコマンドを記述する（NASTRANなら`PARAM,OLDNA,NO`など）。

一点、注意点があります。それは「結果のマッピング」です。リナンバリング後の節点番号で結果が出力されるため、リナンバリング前のオリジナルモデル（CAD連携用の番号）と結果を対応づけるには、ソルバーが内部で保持しているマッピングテーブルが必要です。この対応が失われると、後処理で特定の点の結果を追えなくなります。ただし、通常のプリポスト連携環境では自動で処理されるので、ユーザーが意識することは稀です。

大きな違いは、自動化の度合いとユーザー制御の可否です。Ansys Mechanical はワークベンチ環境下では完全にブラックボックス化されており、ユーザーがアルゴリズムを選択することはほぼできません。一方、MSC Nastran は歴史的にバルクデータデッキを直接触る文化があるため、`PARAM`カードを用いた詳細な制御（`PARAM,PRGPST,2` でGPSアルゴリズム選択など）が可能です。

Abaqus/Standard は非常に強力な自動リナンバリングを内部で行います。ユーザーインターフェースから制御するオプションはほとんどありませんが、その代わり、マルチポイント拘束(MPC)や接触が定義された極めて複雑な接続性を持つモデルに対しても、頑健に帯域幅を最小化します。COMSOL Multiphysics は、選択した物理場と離散化方法に応じて最適な番号付けを自動で決定します。有限要素法と有限体積法が混在するような複合問題では、各物理場ごとに独立した番号付けを行う「ブロック編成」を採用し、全体の行列をブロック構造化することで効率を上げています。

CalculiX (ccx) はデフォルトでRCMを使用します。ソースコードを読むと、`reorder.c`というファイルにアルゴリズムが実装されています。Code_Aster も同様に自動リナンバリングを行いますが、コマンド（`AFFE_MODELE`）の中でメッシュを再番号付けする`RENUM`オプションを明示的に指定することもできます。これらは商用ソルバーに引けを取らない高度なアルゴリズムを実装しています。

直接的なエラーより、「気付かない非効率」が問題です。例えば、大規模モデルでメモリ不足エラーが出た場合、リナンバリングが適切に行われていない可能性があります。Ansysの`SOLVE`コマンド前の`/STATUS, SOLU`出力や、NastranのF06ファイルの「MATRIX STATISTICS」セクションで、行列の項数(Terms)や最大帯域幅を確認するのが第一歩です。項数が想定より一桁多い場合は、リナンバリングが機能していないサインです。

稀ですが、あります。主に2ケース：1) デバッグ時。入力ファイルの特定の節点/要素に問題がある場合、ソルバー出力の番号が書き換えられると原因箇所の追跡が困難になります。Nastranでは`PARAM,OLDNA,YES`で無効化できます。2) 外部ツールとの連携時。自社開発の後処理ツールが固定の節点番号を前提としている場合などです。ただし、この場合は連携ツール側を修正するのが本来の解決策です。

大きく変わります。並列計算では、「帯域幅の最小化」よりも「ドメイン分割」が重要になります。目標は、各プロセッサに割り当てられる自由度の数を均等にしつつ、プロセッサ間のインターフェース（通信が必要な節点）の数を最小化することです。このためのアルゴリズムは「グラフ分割」と呼ばれ、METISやParMETISといったライブラリが事実上の標準です。AnsysやAbaqusの並列ソルバーは、内部でこれらのライブラリを呼び出して、効率的な領域分割を行った上で、各領域内で局所的なリナンバリングを行います。

良い指摘です。メッシュを細かくすると、問題の規模（自由度N）が増えるだけでなく、行列の条件数