4節点四辺形要素(QUAD4)
4節点四辺形要素(QUAD4)の理論基礎
Q4要素 — 2次元FEMの主力
先生、4節点四辺形要素（Q4）って2次元FEMの基本ですよね。
形状関数
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らがバーミンガム大学で開発した。自然座標（ξ, η）を用いることで任意の四辺形形状を単位正方形に写像でき、ガウス積分との組み合わせで高精度な剛性マトリクス計算を実現した。現在も商用コードの基幹要素として採用されている。
4節点四辺形要素(QUAD4)の数値計算手法
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の最推奨 — 迷ったらこれ
Q4要素の2×2ガウス積分
Q4要素には通常2×2（4点）のガウス積分が使われる。積分点座標は±1/√3≈±0.5774で、これにより多項式次数3以下の関数を完全積分できる。1点積分（reduced integration）にすると計算速度が約4倍になる反面、砂時計モードと呼ばれる零エネルギーモードが出現するため、スタビライゼーション技術が必要となる。
4節点四辺形要素(QUAD4)の実務適用
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の閾値で設定するのが自動車部品解析のデファクト標準となっている。
4節点四辺形要素(QUAD4)のソフトウェア比較
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公式を採用している。
4節点四辺形要素(QUAD4)の先端研究
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の先端研究、まとめます。
Q4は50年以上の歴史があるが、まだ改良の余地がある。それだけ「2次元の四辺形要素」が構造力学にとって重要だということだ。
Q4要素の非適合モード拡張
1973年にWilson・Taylorらがカリフォルニア大学バークレー校でQ4要素に非適合変位モード（Wilson-Taylor要素）を追加し、曲げ精度を劇的に改善した。追加自由度は内部で縮合（スタティック・コンデンスレーション）されるため全体剛性マトリクスのサイズは変わらない。ABAQUSのCPS4I、NastranのCQUAD4（PSHELL）の改良版にこの技術が受け継がれている。
4節点四辺形要素(QUAD4)のトラブル対応
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オプションが最も精度が高いとされる。
関連トピック
なった
詳しく
報告