片持ち梁の曲げ（集中荷重） — 実践ガイドとベストプラクティス

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

1. 問題定義: 形状（$L=1$ m, $b=0.1$ m, $h=0.05$ m）、材料（$E=200$ GPa, $\nu=0.3$）、荷重（$P=1000$ N）を明示的に文書化する

2. 理論解の算出: $\delta_{tip}$、$\sigma_{max}$、反力を手計算で求め、参照値として記録

3. メッシュ収束スタディ: 最低3水準のメッシュ（要素サイズ比 $r = 2$ が推奨）で系統的に計算

4. GCIの算出: 観測収束次数と離散化誤差の95%信頼区間を報告

5. 結果の文書化: 入力ファイル、メッシュ、結果をバージョン管理下に置く

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

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

• 全節点固定: 固定端面の全節点で3方向変位をゼロに拘束。最も一般的だが、Poisson効果による横方向の拘束が付随的に生じ、固定端近傍で応力が理論値からずれる
• RBE2/MPC: 端面の節点をマスター節点に剛体結合し、マスター節点を拘束。回転の拘束方法が明確になる
• 分布拘束: 端面に一様変位拘束を適用。Abaqusの COUPLING + KINEMATIC がこれに相当

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

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

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

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

• 問題の明確な定義（形状、材料、荷重、境界条件の図示）
• 使用したソルバーとバージョン番号
• 要素タイプ、メッシュ密度、積分スキームの仕様
• 理論解の導出過程
• メッシュ収束データ（表とグラフ）
• GCIの計算過程と結果
• 入力ファイルの全文または参照先（再現性の保証）

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