片持ち梁の曲げ（集中荷重）

片持ち梁の曲げは、FEA検証の出発点として業界で広く利用されているベンチマーク問題だ。先端たわみ $\delta = PL^3/(3EI)$、固定端最大応力 $\sigma_{max} = PLc/I$ という厳密解が存在するから、数値解法の実装精度を定量的にチェックできる。NAFEMSの入門ベンチマーク集にもこの問題が収録されている。

まさにその通り。ASME V&V 10-2006では、解析コードの数学的正確性を確認する段階としてCode Verificationを位置づけている。厳密解を持つ問題で数値解との一致を示すことが第一歩だ。片持ち梁はその入り口として最適で、Euler-Bernoulli梁理論の仮定が成り立つ範囲なら1次元梁要素で厳密一致する。

Euler-Bernoulli梁理論に基づくたわみ曲線は次の通りだ。

曲げモーメントは固定端で最大 $M_{max} = PL$ だから、最大曲げ応力は固定端の最外縁で生じる。

Timoshenko梁ではせん断変形を考慮する。せん断変形によるたわみ増分 $\delta_s = \kappa PL/(GA)$ が加わる。$\kappa$ はせん断補正係数（矩形断面で $5/6$ ）だ。スパン/せい比 $L/h$ が10以上ならせん断変形の寄与は1%未満になるから、Euler-Bernoulli理論で十分だ。$L/h < 5$ の深梁ではTimoshenkoか3Dソリッド要素が必須になる。

標準設定は $L = 1$ m、$b = 0.1$ m、$h = 0.05$ m、$P = 1000$ N、$E = 200$ GPa、$\nu = 0.3$ だ。このとき理論値は $\delta_{tip} = 0.160$ mm、$\sigma_{max} = 240$ MPa になる。

HEX8は完全積分だとせん断ロッキングが発生し、梁が実際より硬く振る舞う。曲げ支配の問題では低次六面体は本質的に不利で、低減積分やB-bar法を使わない限り収束が遅い。二次要素は中間節点が曲げの変形モードを正確に表現できるから、粗いメッシュでも高精度だ。

ある。FEMの誤差評価定理（Céaの定理から導かれるアプリオリ誤差評価）によれば、$p$ 次要素でエネルギーノルム誤差は $O(h^p)$ で減少する。つまり線形要素なら要素サイズを半分にすると誤差は約半分、二次要素なら約1/4になる。この理論的収束率が実際のメッシュ収束で再現できるかどうかがCode Verificationの本質だ。

典型的なのは応力特異点がある場合だ。片持ち梁の固定端は実際には応力特異点ではないが、固定拘束の実装方法によっては局所的に応力集中が生じて収束率が落ちることがある。対処法としてはSaint-Venantの原理を意識して、拘束端から十分離れた位置で評価するか、分布拘束を使うことだ。

目的によって使い分ける。梁要素ならEuler-Bernoulli理論を直接離散化するから最小の自由度で厳密解に到達する。検証目的でシェルやソリッドを使う場合は、Galerkin法による弱形式の離散化が基本だ。

その通り。片持ち梁程度の規模なら直接法（Cholesky分解）で問題ない。DOFが数万を超えたら前処理付きCG法の方がメモリ効率が良い。この問題の本質は解法ではなく要素定式化の精度確認だから、ソルバー選択よりも要素タイプとメッシュ密度に集中すべきだ。

NastranならCBEAM要素で $w_{tip}$ が厳密一致する。Abaqusだと B31（Timoshenko梁）で同様。ソリッド検証をする場合、AbaqusのC3D20R（二次六面体・低減積分）が鉄板だ。NastranならCHEXA（20節点）。

完全積分のHEX8（2×2×2 Gauss点）は曲げ問題でロッキングする。低減積分（1×1×1）にするとロッキングは解消するが、ゼロエネルギーモード（アワーグラスモード）のリスクがある。B-bar法やEAS（Enhanced Assumed Strain）法はロッキングを回避しつつアワーグラスも抑制する。Abaqusの C3D8I（非適合モード）もこの問題に有効だ。

3水準以上のメッシュで計算し、Richardson外挿で漸近解を推定する。メッシュ比 $r$、2つの解 $f_h$ と $f_{rh}$ から

Grid Convergence Index は Roache が提案した指標で、ASME V&V 20でも採用されている。

Solution Verificationの一環だ。単一ソルバーの結果だけでは実装バグに気づけないことがある。同一メッシュ（UNVやMED形式で共有）を使って Nastran、Abaqus、CalculiX で解き、全員が理論解に収束するか確認する。1つだけ外れていたらそのソルバーの要素定式化か境界条件の適用に問題がある。

節点番号の順序がソルバーごとに異なる。Abaqusは反時計回り、Nastranは時計回りなど。変換時に法線が反転して荷重方向が逆になる事故が起きうる。Gmshで出力形式を切り替えるのが安全だ。また、要素タイプの厳密な対応関係（例: AbaqusのC3D20とNastranのCHEXA-20は節点順序が違う）を確認しておく必要がある。

ASME V&V 10の枠組みに沿って進める。

片持ち梁なら構造格子を使うべきだ。六面体マッピングメッシュで要素サイズを均一にすることで、Richardson外挿の前提条件（均一なメッシュ比）を満たせる。自動テトラメッシュだとメッシュ比が局所的にばらつき、GCIの計算が不安定になる。

本当だ。梁要素なら全自由度拘束で問題ないが、ソリッド要素では実装方法が結果に影響する。

正解はないが、理論解との比較が目的なら「固定端から十分離れた位置で結果を評価する」のが正攻法だ。Saint-Venantの原理により、固定端から梁せい $h$ の2〜3倍離れれば拘束方法の影響はほぼ消える。評価位置を明記して比較することが重要だ。

以下を必ずチェックする。

反力が一致しない場合、境界条件か荷重の設定に根本的な間違いがある。これは数秒で確認できて致命的なバグを検出できるから、最優先で実施すべきだ。新人がよくやるミスとして、Nastranで FORCE カードの座標系を間違えて荷重方向が意図と違うケースがある。反力チェックで即座に発覚する。

NAFEMS QSS（Quality System Supplement）に沿ったレポートには以下を含める。

航空宇宙や原子力の規制対応では必須だ。一般産業でも、5年後に同じ結果を再現できることがV&Vの核心だから、入力ファイルのアーカイブは省略すべきでない。Gitリポジトリへのハッシュ付きリンクでも可だ。

SOL 101（線形静解析）を使う。CBEAM要素なら PBEAM で断面を定義し、SPC1 で固定端を拘束、FORCE で先端荷重を印加する。結果はたわみが DISPLACEMENT、応力が STRESS で .f06 に出力される。

CHEXA（20節点）を使い、荷重をRBE3で分配する。先端面の全節点をRBE3の従属節点にし、独立節点に荷重を印加する。RBE2にすると端面が剛体化して局所応力に影響するから注意。PSOLID カードで積分点数を制御でき、デフォルトの2×2×2（完全積分）が推奨だ。

STEP, NAME=STATIC, PERTURBATION で線形摂動ステップを使う。B31要素なら BEAM SECTION で断面定義、C3D20R なら SOLID SECTION で材料参照を設定。固定端は BOUNDARY で ENCASTRE（全自由度拘束）、荷重は *CLOAD で節点力を印加する。

応力の出力位置だ。Nastranは要素中心の応力をデフォルトで出力するが、Abaqusは積分点の値を出力する。節点外挿のアルゴリズムも異なる。正確な比較をするには、同一座標点での応力を両者から抽出して比較する必要がある。パス上の応力プロットを取ると差異が可視化しやすい。

WorkbenchのStatic Structuralモジュールで設定するのが一般的だが、検証目的ならAPDLコマンドで直接制御する方が透明性がある。BEAM188（Timoshenko梁）かSOLID186（二次六面体）を使う。APDLの ET,1,SOLID186 で要素タイプを指定し、KEYOPT(1,2) で積分オプションを制御する。

その通り。WorkbenchはGUIで簡単に設定できる反面、デフォルト設定が何なのか把握しづらい。V&Vでは全設定を明示する必要があるから、APDL Command Snippetを挿入して要素オプションを明示的に指定するか、最初からAPDLで書く方が確実だ。

CalculiXはAbaqus互換の入力形式を採用しているから、.inp ファイルの書式はほぼ同じだ。C3D20要素を使い、BOUNDARY で固定拘束、CLOAD で荷重印加する。

接触や大変形では定式化の差異が出るが、線形静解析の片持ち梁では両者ほぼ一致する。ただし CalculiX の C3D20 はAbaqusの C3D20 と節点順序が微妙に異なることがあるから、Gmshで出力するときは -format inp オプションで直接CalculiX向けに出すのが安全だ。

同一メッシュ（HEX20、20×4×2）での比較結果はこうなる。

Method of Manufactured Solutions（MMS）だ。片持ち梁の厳密解は特定の荷重・境界条件でしか使えないが、MMSでは任意の解を仮定してそこから逆算した外力項を体積力として与える。こうすれば複雑な解析コード内の離散化の正しさを網羅的に検証できる。

例えば、単純な多項式ではなく三角関数を含む変位場 $u(x,y) = A\sin(\pi x/L)\cos(\pi y/h)$ を仮定し、弾性体の平衡方程式に代入して得られるソース項を体積力として与える。この製造解と数値解の差を評価すれば、要素定式化のバグや積分精度の問題を発見できる。梁理論では検出できないソリッド要素固有のバグ（例えばヤコビアン計算の誤り）を炙り出せる。

段階的に複雑化するのが定石だ。

大たわみでは荷重-変位曲線の全体形状（解析解との比較）、弾塑性では降伏開始荷重と全塑性モーメント $M_p = \sigma_y bh^2/4$、接触では反力分布と境界条件の遷移だ。各段階で理論値と比較できる量が存在することがポイントで、比較対象がない非線形問題にいきなり飛ぶのは危険だ。

CI/CDパイプラインに回帰テストとして組み込むのが最善だ。ソルバーのバージョンアップ時やカスタム要素の開発時に、片持ち梁を含む標準ベンチマーク群を自動実行し、結果が許容範囲内にあることを確認する。

変位の相対誤差 0.1% 以内、応力の相対誤差 1% 以内を目安にする。メッシュは十分に収束した水準で固定する。基準を厳しくしすぎると浮動小数点の丸め誤差やプラットフォーム差異で偽陽性が出るから、有効数字4〜5桁の一致で十分だ。

以下の順で切り分ける。

入力ファイルの冒頭にコメントで使用単位系を明記する慣習をつけることだ。例えば Nastran なら $ UNITS: mm, N, MPa, tonne, s と書いておく。Abaqusには単位系の概念がないから、ユーザーが一貫性を保証する責任がある。最近はCOMSOLのように単位系をGUIで明示管理するソルバーも増えている。

ソリッド要素モデルで固定端面の角部に応力集中が出ている可能性が高い。梁理論は平面保持の仮定に基づくから、端面での局所的な3D効果は含まない。

片持ち梁の固定端は幾何学的には特異点ではないが、境界条件の実装によっては人工的な特異点を作ってしまうことがある。典型的なのは、固定端面の角部の1節点だけに荷重が集中するケースだ。これは節点力の特異点になる。評価位置を角部から離すか、RBE3で荷重を分配すれば解消する。

直接法なら「収束しない」ではなく「特異行列」エラーになる。これは拘束条件の不足で剛体移動が可能な状態を意味する。片持ち梁で起きる典型パターンは:

1. .f06 の AUTOSPC TABLE でどの節点のどの自由度が追加拘束されたか確認

Q: シェル要素でオフセットを使うべきか？