積分点 — CAE用語解説
積分点
先生、積分点ってFEMで重要って聞くんですけど、節点とは何が違うんですか?
理論と物理
積分点の基本概念
「積分点」という言葉をよく聞きますが、具体的に有限要素法の中で何をしている場所なんですか?節点とは違うんですよね?
その通り、節点とは別物です。積分点は、各要素内部に設定された「サンプリング点」です。ここで、要素の変形からひずみや応力を計算し、さらに要素剛性マトリクスを求めるための数値積分を行います。例えば、2次元の4節点アイソパラメトリック要素では、
なぜわざわざ要素内部に点を設ける必要があるんですか?節点の応力を使えばいいのでは?
良い質問だ。節点の変位は連続ですが、節点の応力は要素間で不連続になることが多い。一方、積分点での応力は、要素ごとの構成則を直接評価した「一次」の情報で、最も精度が高い。例えば弾性解析では、積分点で
「ガウス積分点」とありますが、他にも種類があるんですか?位置はどう決めるんですか?
ガウス・ルジャンドル積分が最も一般的だ。これは積分精度が高いからだ。他にも、台形公式に相当する「ロバート積分点」や、要素の角や辺上に点を置く「ロバスト積分」もある。位置と重みは、要素の自然座標系(ξ, η, ζ)で定義される。1次元で2点ガウス積分なら、位置は
数値解法と実装
積分則と要素剛性
積分点の数は、どうやって決めればいいんですか?多ければ多いほど精度が上がる?
必ずしもそうではない。「完全積分」と「低減積分」という概念がある。例えば2次元4節点要素の完全積分は先ほどの2x2=4点だ。これに対し、低減積分は1x1=1点だけ使う。完全積分は剛すぎる(ロッキング)ことがあり、低減積分は柔らかすぎてゼロエネルギー変形(アワーミング)が発生するリスクがある。ソフトウェアでは「選択的低減積分」がよく使われる。
アワーミングって具体的にどういう現象で、どう防ぐんですか?
低減積分要素で、せん断変形や体積変形に寄与しない変形モード(ゼロエネルギー・モード)が発生し、要素が「ぐにゃぐにゃ」になる現象だ。防ぐ方法は主に二つ。一つは「アワーミング制御」を加えること。もう一つは、非線形解析でよく使われる「ハイブリッド要素」だ。これは体積ひずみを積分点とは別に扱い、圧力場を追加の自由度として導入する。AbaqusのCPE4HやC3D8Hがその例だ。
要素剛性マトリクスは、具体的に積分点でどう計算されるんですか?
以下の数値積分で計算される。ここがFEMの核心部分だ。
実践ガイド
メッシュと積分点の選択
実際の解析で、要素タイプを選ぶ時に積分点のことをどう考慮すればいいですか?
まず、線形要素か2次要素かで大きく変わる。Ansysの場合、線形固体要素(SOLID185)はデフォルトで完全積分(8点)だが、キーオプションで低減積分(1点)に切り替えられる。2次要素(SOLID186)は高次なので、低減積分(14点)がデフォルトだ。一般に、曲げが支配的な場合や接触問題では低減積分要素が有利だが、体積ロッキングが懸念される場合は完全積分かハイブリッド要素を選ぶ。
積分点の結果を後処理で確認する必要はありますか?何に気をつければ?
非常に重要だ。特に非線形解析(塑性、超弾性)では、材料の降伏や損傷は積分点で評価される。Abaqus/Viewerで「プローブ値」機能を使い、特定の要素の積分点応力を直接見ることができる。気をつけるのは、積分点値は要素内で変化するため、1つの要素でも複数の値がある点だ。また、極端な変形で積分点が「反転」すると、ヤコビアンの行列式が負になり、解析が停止する。これが「負の体積エラー」の原因の一つだ。
メッシュを細かくすると、積分点の数も単純に増えると考えていいですか?計算コストへの影響は?
その通り、要素数×(1要素あたりの積分点数)で総積分点数は増える。計算コストへの影響は大きい。要素剛性マトリクスの計算、材料の状態更新(非線形)、内部変数の保存は全て積分点で行われる。例えば、100万要素のモデルで要素あたり8積分点を使うと、800万個のデータポイントを毎ステップ処理する。これがメモリとCPU時間を圧迫する。だから、必要な精度を確保しつつ、適切な積分則を選ぶことが最適化の鍵になる。
ソフトウェア比較
各ソルバーでの実装の違い
Ansys、Abaqus、COMSOLで、積分点に関する設定や呼び方に違いはありますか?
大きな違いがある。まず用語。Abaqusは「積分点」、Ansys Mechanical APDLは「積分点」だがWorkbenchでは「ソリューション内点」とも呼ぶ。COMSOLは「積分点」または「ガウス点」だ。設定方法も異なり、Abaqusは要素タイプ定義時に「C3D8R」のように末尾の文字で指定(R:低減積分、なし:完全積分)。Ansys APDLは「KEYOPT」で設定する。COMSOLは「物理場」の「離散化」設定で多項式次数と積分次数を別々に選べるのが特徴だ。
非線形材料の評価は、すべてのソフトで積分点で行われるんですか?
はい、それが標準的な実装だ。例えば弾塑性材料では、各積分点で降伏関数
結果の出力で、積分点データを直接ファイルに書き出す機能はソフトごとに違いますか?
違う。Abaqusでは「*EL FILE」出力要求で、.filファイルや.odbファイルに積分点データ(S, E, PE, PEEQなど)を保存できる。Ansysでは「ETABLE」コマンドで積分点結果をテーブルに読み込むか、「NLIST」と「PRNSOL」を組み合わせる。COMSOLは「データセット」として「積分点」を選択し、そこから派生したデータを評価・エクスポートできる。第三者後処理ツールのEnSightやFieldViewも、これらの生の積分点データを読み込んで可視化できる。
トラブルシューティング
積分点に関連するエラーと対策
解析中に「積分点で塑性ひずみが発散しました」といったエラーが出ます。原因は何ですか?
主に三つの原因が考えられる。第一に、局所的なひずみ集中で、現実にはあり得ない大きな変形が生じ、材料モデルの適用限界を超えている。第二に、材料パラメータ(硬化則など)の設定ミス。第三に、大きな時間増分(Δt)による数値的不安定さだ。対策としては、メッシュを細かくする、材料パラメータを見直す、ステップの「自動増分」設定をより厳しくする(最大増分サイズを小さくする)などがある。
低減積分要素を使ったら、「ゼロエネルギー変形モードが検出されました」と警告が出ました。無視していいですか?
絶対に無視してはいけない。この警告は、モデルが物理的に意味のない変形モードで柔らかく(剛性が低く)なっている可能性を示す。放置すると、応力が過小評価されたり、収束が異常に早くなったりする。対策は、アワーミング制御(Abaqusの「Hourglass stiffness」など)を強化する、メッシュを変更する(例えば、1層だけの低減積分要素は危険)、あるいは問題の領域に完全積分要素やハイブリッド要素を使用するだ。
後処理で、応力のコンター図が要素内で「まだら」になったり、隣接要素間でギャップが生じたりします。これも積分点が関係していますか?
まさにその通りだ。これは積分点の値(不連続)を節点に外挿・平均化する過程で生じる。特に低減積分要素では積分点が少ないので、外挿の影響が大きい。対策は二つ。一つは、後処理設定で「平均化」の方法を変えること(Abaqusでは「平均化閾値」を調整)。もう一つは、解析そのものの精度を上げることだ。つまり、メッシュを細かくするか、2次要素を使う。2次要素は積分点が多いため、より滑らかな応力分布を節点に外挿できる。
過酷な接触条件で「負のヤコビアン」エラーが出ました。積分点と関係ありますか?
直接的な原因だ。要素が極度に変形し、積分点でのヤコビアン行列式
関連トピック
なった
詳しく
報告