平面応力問題

そう。3次元の弾性体を解くと自由度が膨大になる。しかし薄い板状の構造で面内荷重だけ受ける場合、板厚方向の応力成分がゼロと見なせる。これが平面応力の仮定だ。

板厚が面内寸法に比べて十分薄い構造だ。具体的には：

板厚方向に拘束がある場合だ。例えばダムの厚い壁や長いトンネルの断面。こういう場合は平面ひずみ（板厚方向のひずみがゼロ）の仮定を使う。平面応力と平面ひずみは似て非なるもので、混同すると大きな誤差になる。

2次元の力の平衡方程式：

平面応力のフックの法則（マトリクス形式）：

いい指摘だ。ただし平面応力では横方向（$z$ 方向）の変形は自由だ。$E/(1-\nu^2)$ が現れるのは、2成分の応力が連成するからで、3次元の $E$ とは意味が少し違う。

重要なポイント：平面応力では $\sigma_{zz} = 0$ だが、$\varepsilon_{zz} \neq 0$ だ：

そう。$z$ 方向の変形が完全に拘束されるから、3次元的な拘束効果で剛性が上がる。この差は $\nu = 0.3$ で約10%。小さいように見えるが、応力の評価では無視できない差だ。

Airyの応力関数 $\phi(x,y)$ を使うと、平衡方程式と適合条件を1つの方程式にまとめられる：

古典的な弾性論の問題（穴のある無限板、くさび、半平面の接触問題）はAiry関数で解ける。Kirsch問題（無限板の円孔周りの応力集中）も Airy関数の解だ。FEMの検証問題として非常に有用だよ。

要点：

その通り。2次元要素は計算コストを劇的に削減できるが、仮定が成り立つ場面を正しく見極めることが前提だ。疑わしければ3次元で解いて2次元の結果と比較すればいい。

2次元の平面応力要素だ。各節点に2つの自由度（$u_x, u_y$）を持ち、3次元ソリッド要素の約1/3のDOFで済む。

その通り。CSTは応力集中部のような応力勾配が大きい場所では著しく精度が低い。教育用や粗いスクリーニングにしか使えない。実務では二次要素（6節点三角形 or 8節点四角形）が標準だ。

Q4（4節点四角形）の場合：

曲げ変形を受けるとき、Q4の完全積分では寄生的なせん断ひずみが発生して要素が「ロック」する。実際にはほとんど曲がらない状態になる。これがシアロッキングだ。

対策：

低減積分の積分点（1点）で応力がゼロのまま、要素が砂時計のようにジグザグに変形するモード。エネルギーがゼロなので力学的には存在しないはずの変形だが、数値的には現れてしまう。アワーグラス制御（人工的な剛性追加）で抑制する。

そう。「R = Reduced integration」。AbaqusではCPS4が完全積分、CPS4Rが低減積分。実務ではCPS4R（低減積分＋アワーグラス制御）かCPS4I（非適合モード）が推奨される。

NastranのCQUAD4は内部的に非適合モードを含んでおり、シアロッキングに対してロバスト。デフォルトのまま使って問題ないことが多い。

そう。Nastranではプロパティカード（PSHELL）で板厚を指定。Abaqusでは *SOLID SECTION で板厚を定義。Ansysでは REAL定数やSection で設定。板厚が1.0のままの解析結果は力の次元が合わないから、結果を見れば気づけるはず。

要点：

3次元解析が当たり前の現在でも、平面応力解析には使い道がある。

計算速度だ。2次元は3次元の1/100以下の計算コスト。設計パラメータを100通り変えたいとき、3次元で100ケース回すのは現実的でないが、2次元なら数分で済む。

最も典型的な適用だ。円孔のある無限板の応力集中はKirsch問題で $K_t = 3.0$（一軸引張）。これをFEMで再現するのは平面応力の基本演習だ。

実務的なポイント：

Petersonのチャートは多くの形状の $K_t$ を整理した名著だ。FEMはPetersonに載っていない形状（非標準的な切り欠き、複合形状等）の $K_t$ を個別に評価できる。Petersonの値とFEMの結果を比較するのは検証の基本だ。

3つのポイントがある：

1. 剛体移動の拘束 — 2次元なら最低3自由度の拘束で剛体移動を止める（並進2 + 回転1）

2. 対称条件の活用 — 形状と荷重が対称なら1/2, 1/4モデルにできる。計算コスト削減と同時にメッシュの精緻化が可能

3. Saint-Venantの原理 — 荷重を与える辺が荷重の細部に敏感。着目部位から十分離して荷重を与えれば、荷重の細部は影響しない

例えば引張試験片の中央部の応力を知りたいとき、端部の拘束条件が多少不正確でも、中央部の応力は正しく出る。ただし「十分離れる」距離は板幅の1〜2倍程度が目安だ。

1. 反力の確認 — 反力の合計が荷重と釣り合っているか

その通り。AbaqusならCPS（平面応力）とCPE（平面ひずみ）の1文字の差、AnsysならKEYOPTの設定の差。たった1つの設定で物理が変わるから、最初に必ず確認すること。

平面応力は全ての構造解析ソルバーの基本機能だ。差が出るのは要素の品質と使い勝手だ。

NastranのCQUAD4/CQUAD8は非常に高品質な要素として知られている。特にCQUAD4は内部的に非適合モードを含んでおり、4節点ながら曲げ精度が高い。平面応力にはPSHELLプロパティで板厚を指定し、MID3（面内剛性）を設定する。

AbaqusのCPS（Continuum Plane Stress）ファミリーは豊富：

実務推奨はCPS8R（二次四角形、低減積分）またはCPS4I（非適合モード）。

AnsysのPLANE182（線形）/ PLANE183（二次）。KEYOPTで平面応力/平面ひずみ/軸対称を切り替える。KEYOPTの設定を間違えないことが最重要。

いくつかある：

FEMの基本を学ぶには2次元平面応力が最適な題材だ。フリーソフトで孔あき板の応力集中を再現する演習は、FEMの理解に非常に効果的だ。

2次元の平面応力自体は完成した理論だが、その拡張や応用には新しい展開がある。

従来の平面応力は板厚方向に一様な応力を仮定するが、一般化平面応力（generalized plane stress）は板厚方向の応力変化を考慮する。厚板や板のエッジ近くでは $\sigma_{zz} \neq 0$ だから、通常の平面応力では不正確になる。

Asymptoticな手法で板厚方向の応力分布を展開し、平面問題として解く。板理論（Kirchhoff, Mindlin）の3次元補正としても使われる。

2次元の亀裂問題はXFEMの最も基本的な応用だ。従来のFEMでは亀裂先端にメッシュを細分化する必要があるが、XFEMでは亀裂がメッシュを貫通してもOK。

Abaqusの ENRICHMENT + CONTOUR INTEGRAL で2次元亀裂のJ積分やSIFが計算できる。メッシュ生成の手間が大幅に減る。

ペリダイナミクス（Peridynamics）は微分方程式ではなく積分方程式で連続体を記述する手法だ。2次元平面応力問題に適用すると、亀裂の発生と進展がメッシュフリーで自然に表現される。

従来のFEMでは亀裂先端に特異性があるが、ペリダイナミクスでは力の伝達を「ボンド」として記述し、ボンドが切れることで亀裂を表現する。2次元の脆性破壊シミュレーションで急速に研究が進んでいる。

もう一つ、位相場法が2次元の破壊問題で注目されている。損傷を連続的な場（0 = 健全、1 = 完全破壊）で表現し、亀裂の分岐・合流を自然に追跡できる。FEMのフレームワークで実装可能で、Abaqusのユーザーサブルーチン（UEL/UMAT）で使える。

まさにそう。新しい手法は必ず2次元平面応力で検証してから3次元に拡張する。2次元は計算が軽いから、手法の妥当性を素早く確認できる。平面応力は「枯れた理論」であると同時に、最先端の研究を支える基盤でもあるんだ。

平面応力は基本的な解析だが、初心者がつまずきやすいポイントがある。

間違いなく最多だ。結果の影響：

確認方法：

$\nu = 0.5$ では体積変化がゼロ（非圧縮）の拘束が入る。平面ひずみの完全積分要素では拘束条件が過剰になり、体積ロッキングが起きる。低減積分要素やハイブリッド要素（Abaqusの CPE4H）を使う必要がある。

原因チェック：

変わる。FEMの応力は積分点で最も正確だ。節点平均値（各要素の値を平均したもの）は滑らかになるが、応力集中のピーク値を過小評価する。非平均化（unaveraged）の応力コンターを見て、要素間の不連続がないかチェックすること。

確認項目：

2次元要素の品質指標：

要素が裏返っている（節点の順序が逆）ことを意味する。計算は走るかもしれないが結果は完全に間違い。自動メッシュ後に必ず品質チェックを行うこと。

FEMの全ての解析は基本の上に成り立っている。2次元平面応力で基本を徹底することが、複雑な3次元解析で間違いを防ぐ最良のトレーニングだ。