4節点四辺形要素(QUAD4)
理論と物理
Q4要素 — 2次元FEMの主力
先生、4節点四辺形要素（Q4）って2次元FEMの基本ですよね。
形状関数
Q4の形状関数は自然座標 $(\xi, \eta)$ で双線形（bilinear）：
$$ N_i = \frac{1}{4}(1 + \xi_i \xi)(1 + \eta_i \eta) $$
Q4の形状関数は自然座標 $(\xi, \eta)$ で双線形（bilinear）：
4つの節点は $(\xi, \eta) = (\pm 1, \pm 1)$ に対応。
双線形ということは、$1, \xi, \eta, \xi\eta$ の4項。交差項 $\xi\eta$ がある。
この交差項のおかげで、Q4は3節点三角形（CST）より曲げの表現能力が高い。ただし完全な曲げ変形を表現するには $\xi^2, \eta^2$ 項が必要で、Q4にはこれがない。
シアロッキングの問題
Q4のシアロッキングについて教えてください。
完全積分（2×2 Gauss点）のQ4で純粋曲げを表現しようとすると、寄生的なせん断ひずみが発生して要素が「固まる」。これがシアロッキング。
物理的にはこうなる：
- 純粋曲げでは上面が伸び、下面が縮む
- Q4の双線形形状関数ではこの変形を表現するとき、不可避的にせん断ひずみが出る
- このせん断ひずみのエネルギーが余分に蓄えられ、要素が硬くなりすぎる
対策は？
3つの方法：
1. 低減積分（1×1 Gauss点） — せん断を積分点1つで評価。ロッキングは消えるがアワーグラスモードが発生
2. 非適合モード（Incompatible modes） — $\xi^2, \eta^2$ に対応する内部自由度を追加。NastranのCQUAD4、AbaqusのCPS4I
3. Assumed Natural Strain (ANS) — せん断ひずみを別途仮定
NastranのCQUAD4はデフォルトで非適合モードが入っているんですよね。
そう。NastranのCQUAD4はFEMの歴史で最も成功した要素の一つと言われている。MacNealとHarderが1985年に発表した改良版で、非適合モードにより曲げ精度を大幅に改善しつつ、パッチテストも満たす。
Q4の長所と短所
| 特性 | 長所 | 短所 |
|---|---|---|
| メッシュ生成 | マップドメッシュが容易 | 自由メッシュでは三角形も混ざる |
| 精度 | 非適合モードで高精度 | 歪んだ形状で精度低下 |
| 計算コスト | 低い（8 DOF） | — |
| 曲面近似 | — | 直辺のみ（曲面は折れ線近似） |
まとめ
Q4の理論を整理します。
要点：
- 双線形形状関数 — 交差項 $\xi\eta$ あり。CSTより精度高い
- シアロッキング — 完全積分で発生。低減積分 or 非適合モードで対策
- NastranのCQUAD4は非適合モード内蔵 — FEM史上最も成功した要素の一つ
- 2次元の実務標準 — 平面応力、平面ひずみ、軸対称で広く使用
- 歪んだ形状に弱い — アスペクト比、スキューネスの管理が重要
HEX8のページで学んだロッキングやアワーグラスの概念が、Q4でもそのまま当てはまるんですね。
Q4はHEX8の2次元版だから、全く同じ問題が起きる。2次元で理解しておけば、3次元に拡張するのは容易だ。
Q4要素の等パラメトリック定式
Q4要素の等パラメトリック定式は1960年代にIrons・Ergatoudisらがバーミンガム大学で開発した。自然座標（ξ, η）を用いることで任意の四辺形形状を単位正方形に写像でき、ガウス積分との組み合わせで高精度な剛性マトリクス計算を実現した。現在も商用コードの基幹要素として採用されている。
各項の物理的意味
- 慣性項（質量項）：$\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で統一 |
数値解法と実装
Q4の積分スキーム比較
Q4の各積分スキームの特徴を教えてください。
アワーグラスが3モードってことは、8自由度中3つがゼロエネルギー？
Q4の8自由度から剛体モード3つ（並進2＋回転1）を引いた5つの変形モードのうち、低減積分で評価されるのは2つだけ。残り3つがアワーグラスモードだ。
ソルバー別の要素名
| バリエーション | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 標準Q4 | CQUAD4 | CPS4（平面応力） | PLANE182 |
| 低減積分 | — | CPS4R | PLANE182(red.) |
| 非適合モード | CQUAD4（デフォルト） | CPS4I | PLANE182(EAS) |
| ハイブリッド | — | CPE4H（平面ひずみ） | — |
NastranのCQUAD4はデフォルトで非適合モードが入っているから、特に設定不要ですね。
Nastranユーザーは「CQUAD4を使えばOK」という簡潔な世界。AbaqusやAnsysではユーザーが積分スキームを選択する必要がある分、知識が求められる。
メッシュ品質の影響
Q4は形状の歪みに敏感ですか？
非常に敏感だ。理想的な正方形Q4と、歪んだ平行四辺形Q4では精度に大きな差が出る。
| 品質指標 | 理想値 | 許容範囲 | 精度への影響 |
|---|---|---|---|
| アスペクト比 | 1.0 | < 5.0 | 大きいと精度低下 |
| スキューネス | 0° | < 45° | 大きいと精度低下 |
| テーパー比 | 0 | < 0.5 | 大きいと精度低下 |
| 内角 | 90° | 45°〜135° | 範囲外で精度低下 |
内角が45°以下のつぶれた四辺形はダメということですね。
そう。つぶれた四辺形はヤコビアンが場所によって大きく変わり、積分精度が落ちる。非適合モード要素でも歪んだ形状では精度が保証されない。Q4のメッシュ品質管理は2次元FEMで最も重要な作業だ。
パッチテスト
パッチテストって何ですか？
歪んだ要素のパッチ（集合体）に一様ひずみを与えたとき、正確に一様応力が出るかを確認するテスト。FEM要素の最低限の品質保証だ。
非適合モード要素はパッチテストを厳密には満たさないことがある（MacNeal-Harderのパッチテストは通るが、一般の歪んだパッチでは微小な誤差が出る）。実用上は問題にならないが、理論的に気にする人はいる。
まとめ
Q4の数値手法、整理します。
要点：
- NastranのCQUAD4はデフォルトで非適合モード — 設定不要で高精度
- AbaqusではCPS4I（非適合モード）を選択 — CPS4R（低減積分）も可だがアワーグラスに注意
- メッシュ品質がQ4の精度を支配 — アスペクト比、スキューネスを厳しく管理
- パッチテスト — FEM要素の最低限の品質保証
- Q4 + 非適合モードが2次元FEMの最推奨 — 迷ったらこれ
Coffee Break よもやま話
Q4要素の2×2ガウス積分
Q4要素には通常2×2（4点）のガウス積分が使われる。積分点座標は±1/√3≈±0.5774で、これにより多項式次数3以下の関数を完全積分できる。1点積分（reduced integration）にすると計算速度が約4倍になる反面、砂時計モードと呼ばれる零エネルギーモードが出現するため、スタビライゼーション技術が必要となる。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
パッチテストって何ですか？
歪んだ要素のパッチ（集合体）に一様ひずみを与えたとき、正確に一様応力が出るかを確認するテスト。FEM要素の最低限の品質保証だ。
非適合モード要素はパッチテストを厳密には満たさないことがある（MacNeal-Harderのパッチテストは通るが、一般の歪んだパッチでは微小な誤差が出る）。実用上は問題にならないが、理論的に気にする人はいる。
Q4の数値手法、整理します。
要点：
Q4要素の2×2ガウス積分
Q4要素には通常2×2（4点）のガウス積分が使われる。積分点座標は±1/√3≈±0.5774で、これにより多項式次数3以下の関数を完全積分できる。1点積分（reduced integration）にすると計算速度が約4倍になる反面、砂時計モードと呼ばれる零エネルギーモードが出現するため、スタビライゼーション技術が必要となる。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
Q4の実務適用
Q4は2次元解析のどこで使われていますか？
Q4 vs. Q8 の使い分け
Q4と8節点四辺形（Q8）のどちらを使うべきですか？
| 判断基準 | Q4 | Q8 |
|---|---|---|
| 応力集中の評価 | メッシュを細かく必要 | 粗めのメッシュでも正確 |
| 曲面の近似 | 直辺（折れ線近似） | 曲辺（二次近似） |
| DOF | 8（4節点×2） | 16（8節点×2） |
| 計算コスト | 低い | Q4の2〜3倍 |
| 推奨場面 | メッシュが十分細かい場合 | 精度重視、メッシュが粗い場合 |
メッシュが十分細かいならQ4でもQ8でも結果は同じ？
収束した結果は同じだ。ただしQ4は収束が遅い（$O(h)$の応力収束）のでQ8（$O(h^2)$）の2〜5倍のメッシュが必要。トータルのDOF数ではQ8のほうが効率的なことが多い。
じゃあQ8を使えばいいのでは？
Q8が優れているのは間違いないが、Q4のメリットは自動メッシュ生成との相性だ。Q4の自動メッシュは安定しているが、Q8の自動メッシュは中間節点の配置で問題が起きることがある。また陽解法ではQ4が標準（安定時間増分が大きい）。
メッシュ生成のベストプラクティス
Q4メッシュ生成のポイントは？
O-gridって何ですか？
穴の周囲に放射状・同心円状のQ4メッシュを配置するパターン。穴の応力集中を正確に捕捉するための古典的なメッシュ技法だ。HyperMeshやAbaqus/CAEのパーティション機能で作れる。
実務チェックリスト
Q4のチェックリストをお願いします。
三角形への退化はなぜダメなんですか？
Q4を三角形に退化させると（2つの節点を同じ位置に）、ヤコビアンが退化点でゼロになり、積分精度が著しく落ちる。三角形が必要な場所では、専用の3節点三角形要素（CST/LST）を使うべきだ。
Q4要素のアスペクト比管理
NASTRANのQUAD4要素ではアスペクト比5以下を推奨する。実務上は薄板構造で面内応力を求める場合、要素辺長比3以内に収めるとWarning 3035の抑制と精度向上が両立できる。Altair HyperMeshのQuality Indexでは「Aspect Ratio」チェックを5.0の閾値で設定するのが自動車部品解析のデファクト標準となっている。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
Q4の各ソルバーでの特徴
各ソルバーのQ4要素に差はありますか？
「ドリリングDOF」って何ですか？
Q4は面内の2自由度（$u, v$）しか持たないが、フレーム解析のQ4シェルでは面外の回転自由度（drilling rotation）が必要になることがある。NastranのCQUAD4はドリリングDOFを持っており、フレーム要素（梁要素）との接続が自然にできる。
選定ガイド
Q4に関する選定ガイドは？
Q4は「FEMの基礎の基礎」。どのソルバーでも必ず使う要素ですね。
そう。Q4を正しく使えるかどうかが、FEMエンジニアの基本スキルだ。
各ソルバーのQ4要素実装比較
NastranのCQUAD4、AbaqusのS4R、AnsysのSHELL181はいずれも4節点シェルだが内部定式が異なる。S4Rは減次積分＋砂時計制御、SHELL181はReissner-Mindlinベースでせん断ロッキング対策を内蔵する。LS-DYNAのType2シェルは1点積分をデフォルトとし、高速衝突解析に特化したBelytschko-Lin-Tsay公式を採用している。
選定で最も重要な3つの問い
- 「何を解くか」：4節点四辺形要素(QUAD4)に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
Q4の先端研究
Q4はこれ以上改良の余地がありますか？
50年以上研究されてきた要素だが、まだ改良は続いている。
Smoothed FEM (S-FEM) の2D版
Edge-based S-FEM (ES-FEM)をQ4に適用すると、通常のQ4より精度が上がり、Q8に近い結果が得られる。ひずみを辺ベースで平滑化することで、交差項以上の変形モードを暗黙的に表現する。
Polygonal FEM
四辺形に限らず、5角形、6角形、N角形の2次元要素を使う多角形要素法。Voronoi分割で自然に生成されるメッシュを使えるため、メッシュ生成の制約が大幅に緩和される。
Q4もQ8も超えた一般的な多角形要素…。
TaliscaやWachspress基底関数で任意多角形の形状関数を構成する。まだ研究段階だが、将来的にはQ4/Q8の代替になる可能性がある。
物理ベース形状関数
通常のQ4は多項式の形状関数だが、支配方程式の解を形状関数に使う手法がある。例えばラプラス方程式の基本解を形状関数にした「Trefftz要素」。特定の問題（応力集中、亀裂先端）で通常のQ4より桁違いの精度が出る。
まとめ
Q4の先端研究、まとめます。
Edge-based S-FEM (ES-FEM)をQ4に適用すると、通常のQ4より精度が上がり、Q8に近い結果が得られる。ひずみを辺ベースで平滑化することで、交差項以上の変形モードを暗黙的に表現する。
四辺形に限らず、5角形、6角形、N角形の2次元要素を使う多角形要素法。Voronoi分割で自然に生成されるメッシュを使えるため、メッシュ生成の制約が大幅に緩和される。
Q4もQ8も超えた一般的な多角形要素…。
TaliscaやWachspress基底関数で任意多角形の形状関数を構成する。まだ研究段階だが、将来的にはQ4/Q8の代替になる可能性がある。
通常のQ4は多項式の形状関数だが、支配方程式の解を形状関数に使う手法がある。例えばラプラス方程式の基本解を形状関数にした「Trefftz要素」。特定の問題（応力集中、亀裂先端）で通常のQ4より桁違いの精度が出る。
Q4の先端研究、まとめます。
Q4は50年以上の歴史があるが、まだ改良の余地がある。それだけ「2次元の四辺形要素」が構造力学にとって重要だということだ。
Q4要素の非適合モード拡張
1973年にWilson・Taylorらがカリフォルニア大学バークレー校でQ4要素に非適合変位モード（Wilson-Taylor要素）を追加し、曲げ精度を劇的に改善した。追加自由度は内部で縮合（スタティック・コンデンスレーション）されるため全体剛性マトリクスのサイズは変わらない。ABAQUSのCPS4I、NastranのCQUAD4（PSHELL）の改良版にこの技術が受け継がれている。
トラブルシューティング
Q4のトラブル
Q4要素でよくあるトラブルを教えてください。
2次元FEMの基本的なトラブルがQ4に集約される。
変位が過小（シアロッキング）
曲げ問題で変位が理論値より小さいです。
完全積分のQ4でシアロッキングが起きている。対策は3次元のHEX8と同じ：
- CPS4I（非適合モード）に変更
- CPS4R（低減積分）に変更
- Q8（CPS8R）に切り替え
アワーグラスモード（変位がジグザグ）
低減積分のQ4で変位が波打ちます。
集中荷重でアワーグラスモードが励起されている。荷重を分散させるか、CPS4I（非適合モード）に切り替える。
応力コンターに「チェッカーボード」パターン
応力が市松模様になります。
2つの可能性：
1. 低減積分のアワーグラス — CPS4Iに切り替え
2. メッシュが歪んでいる — 要素品質を確認。アスペクト比 > 5 や スキューネス > 45° がないか
板厚の設定忘れ
結果の力のオーダーがおかしいです。
2次元平面応力要素では板厚の設定が必要。デフォルトの1.0のままだと、実際の板厚と異なるため力やたわみのオーダーが合わない。Abaqusの *SOLID SECTION で板厚を必ず指定すること。
境界条件での過拘束
固定端のつもりで全節点の全自由度を拘束したら、応力がおかしくなりました。
ポアソン効果で横方向に膨張する変形が拘束されてしまっている。「固定」の意味を明確にすること：
- 面外変位のみ拘束 — 最も一般的
- 面内も拘束 — ポアソン膨張も抑制。実構造に合っているか確認
まとめ
Q4のトラブル対処、整理します。
Q4のトラブルはHEX8と同じ構造を持つ。2次元で理解すれば3次元にも応用できる。
Q4要素の砂時計モード診断
1点積分Q4要素で砂時計（hourglass）モードが発生すると、隣接要素で変位が交互に正負に振れる模様がコンター図に現れる。MSC NastranではDIA=16のフラグで零エネルギーモードを検出できる。Abaqusでは人工的な砂時計制御スティフネスをデフォルト採用しており、HOURGLASS=ENHANCEDオプションが最も精度が高いとされる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——4節点四辺形要素(QUAD4)の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告