平面応力問題 — 先端技術と研究動向

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

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

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

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

• 亀裂先端のエンリッチ関数 — $\sqrt{r}$ の特異場を陽に含める
• Heaviside関数による不連続 — 亀裂面でのジャンプを表現
• メッシュの再生成不要 — 亀裂進展をメッシュ非依存で追跡

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

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

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

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

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

• 一般化平面応力 — 板厚効果の精密な考慮
• XFEM — メッシュ非依存の亀裂解析。2次元が基本
• ペリダイナミクス — 積分方程式ベースの新しい破壊手法
• 位相場法 — 連続場で亀裂を追跡。FEMフレームワークで実装可能

2次元平面応力は「古い」のではなく、「基本」だからこそ最先端の手法開発に使われ続けている。