NAFEMS FV32: 厚肉円板の自由振動

いい質問だね。FV32は厚肉円板の自由振動を扱うベンチマークだよ。具体的には、面外曲げの固有振動数をFEMで求めて、理論解 $f_1 = 1.4568\,\text{Hz}$ と比較する問題なんだ。

そこが核心だ。厚板ではせん断変形の影響が大きくなる。薄板の前提で使うKirchhoff理論では「断面は変形後も中立面に垂直なまま」と仮定するけど、厚板ではそれが成り立たない。Kirchhoff理論をそのまま使うと固有振動数を過大評価してしまうんだ。

FV32の条件（板厚/半径 = 0.1）だと、Kirchhoff理論で求めた1次固有振動数はMindlin理論より数%高く出る。実務で「FEMの結果が理論より高い」と思ったら、せん断変形の効果を考慮しているかを真っ先に疑うべきだね。だからこそFV32は、Mindlin板要素やソリッド要素の精度検証に広く使われているんだよ。

シンプルだけど奥が深い設定だよ。

その通り。無次元量として重要なのは $t/R = 0.1$ という比率だけだ。この値は「明確に厚板だが、極端に厚くはない」という領域で、薄板理論と厚板理論の差が顕著に出る絶妙な設定なんだよ。例えば自動車のブレーキディスクなんかもだいたい $t/R = 0.05\sim0.15$ くらいの範囲だから、実務的にも意味のある板厚比なんだ。

まずKirchhoff（古典的薄板理論）の基本仮定は「変形後も断面は中立面に垂直」、つまりせん断変形ゼロという仮定だ。これにより断面の回転角 $\theta$ が面外変位 $w$ の勾配で一義的に決まる：

まさにそのイメージだ。厚い本を曲げてみると、ページ（断面）が表紙に対して垂直を保てずに傾くだろう？ あれがせん断変形の効果だ。薄い紙1枚なら断面は常に垂直だけど、厚くなると無視できなくなる。

完璧な理解だよ。剛性が過大 → 振動数が高く出る。これはFEM実務でも非常に重要なポイントで、板・シェル要素を選ぶときは必ず「Mindlinタイプかどうか」を確認すべきなんだ。

Mindlin板の自由振動を円板座標で書くと、軸対称モード（周方向波数 $n = 0$）の場合、面外変位 $w(r,t)$ と断面回転角 $\psi(r,t)$ の連立方程式になる：

いい着眼点だね。Mindlin理論では板厚方向のせん断応力分布を一定と仮定するけど、実際は放物線分布だ。その差を補正するのが $\kappa$ で、矩形断面では厳密に $5/6$ になる。ちなみに円形断面の梁では $\kappa = 6/7$ を使うこともあるけど、板では慣例的に $5/6$ を使うよ。

NAFEMSが公表している1次固有振動数の参照解は：

Kirchhoff薄板理論による単純支持円板の1次固有振動数は：

自由振動だから外力ゼロ。変位を $\{u\} = \{\phi\} e^{i\omega t}$ と仮定すると、一般化固有値問題に帰着する：

鋭い質問だね。一般にconsistent質量マトリクスの方が精度は高い。ただしlumped質量は対角化されるため計算が軽い。FV32程度の問題規模ではconsistentを推奨する。大規模モデル（数百万DOF）でメモリが厳しい場合はlumpedでも実用的な精度が得られるよ。実際の差は0.5%未満であることが多い。

どちらでも正解が出せるけど、それぞれ注意点が違うんだ。

その通り。FV32はまさにせん断変形が主役の問題だから、ロッキングが起きると固有振動数が理論値より大幅に高くなる。完全積分の1次要素（CQUAD4完全積分、C3D8完全積分）は特に要注意だ。低減積分またはB-bar法を使う、もしくは最初から2次要素を選ぶのがベストプラクティスだよ。

Mindlin板要素（CQUAD8相当）を使った場合の典型的な収束データを見てみよう。

そうなんだ。振動解析では2次要素のコストパフォーマンスが圧倒的に高い。これはFV32が如実に示している教訓だね。実務でも「まず2次要素で少ないメッシュ → 収束確認 → 必要なら細分化」がゴールデンルートだよ。

よし、ステップバイステップで行こう。

一番多い失敗は2つ。1つ目は単位系の不整合。$E$ をGPaで入れたのに密度をkg/m³ではなくtonne/mm³で入れ忘れる、みたいなやつだ。固有振動数は $f \propto \sqrt{E/\rho}$ だから、密度が1000倍ズレると周波数が約30倍変わって「なんだこの数値は」ってなる。

2つ目は境界条件の解釈ミス。「単純支持」と言っても、板要素の場合は面外変位 $w = 0$ だけ拘束して回転は自由にする。面内変位まで拘束してしまうと「固定支持」寄りの結果になって固有振動数が高く出る。ソリッド要素の場合は外周面の全節点でz方向変位のみゼロにするのが正しい設定だよ。

その通り。軸対称モード（$(0,1)$ モード）だけを求めるなら、2次元軸対称要素（CAX8、PLANE183軸対称など）で半径方向の断面だけをモデル化すれば十分だ。DOFが劇的に減る。

主要4ソルバーでの設定を比較しよう。

Lanczos法は最もロバストで、FV32のような中小規模問題ではほぼ一択だ。大規模問題（100万DOF超）ではAMLS（Automated Multi-Level Substructuring）やShift-Invert法など、より効率的な手法が使われることもある。ただしFV32レベルの問題ではどの手法でも一瞬で終わるから、手法の差を気にする必要はないよ。

適切な要素を選べば、どのソルバーでも参照解 $f_1 = 1.4568\,\text{Hz}$ の誤差0.1%以内に収まる。これが「ベンチマーク」の存在意義だね。

FV32のように「問題設定がクリーン」な場合はソルバー間差は極めて小さい。差が出るのは、メッシュが粗かったり、せん断ロッキング対策が不十分な要素を使ったときだ。つまり「ソルバーの差」ではなく「ユーザーの設定の差」が支配的ということだね。これはCAE実務でも同じことが言えるよ。

あるよ。代表的なのはReddy（1984）の3次せん断変形理論（TSDT: Third-order Shear Deformation Theory）だ。Mindlin理論ではせん断ひずみの板厚方向分布を一定と仮定するけど、TSDTでは3次多項式で近似する。これによりせん断補正係数 $\kappa$ が不要になるんだ。

$t/R > 0.15$ くらいの極厚板や、積層複合材料（CFRPなど）の場合だ。複合材は層間のせん断剛性が各層で異なるため、Mindlinの一定せん断仮定では不十分になる。FV32自体はMindlinで十分だけど、「その先」を知っておくと実務の幅が広がるよ。

円板のような滑らかな形状にはIGAは非常に相性がいい。NURBSベースの基底関数を使うから、曲線境界の近似精度がメッシュに依存しないんだ。研究レベルでは従来FEMより少ないDOFで同等以上の精度が報告されている。ただし2026年現在、IGAを標準でサポートする商用ソルバーはまだ限定的だね。LS-DYNAやCOMSOLが一部対応している程度だ。

1.50Hz ということは約3%高い。これは典型的な「Kirchhoff要素を使ってしまった」パターンだね。使った要素は何かな？

ビンゴ。CTRIA3は薄板理論ベースでせん断変形を考慮しない。しかも1次要素だからメッシュ感度も高い。CQUAD8に変えてやり直してごらん。

そう。FEMのデバッグは「ズレ方のパターン」を知っているかどうかが全てだ。固有振動数が高い → 剛性が過大（ロッキング、拘束過多、要素理論の制約）。固有振動数が低い → 剛性不足または拘束不足。この原則はFV32に限らず、あらゆる振動解析で使えるよ。

絶対に確認すべきだ。FV32の1次モードは軸対称の面外曲げ、つまり「お椀型」に変形するモードだ。もし面内振動モードや反対称モードが1次に出てきたら、境界条件の設定ミスの可能性が高い。

その通り。「結果をポスト処理で必ず可視化する」はCAEの鉄則だ。FV32は問題が単純だからこそ、モード形状の異常にすぐ気づける。これを機に「必ずモード形状を確認する」習慣を身につけてほしい。実務の複雑なモデルではさらに重要になるからね。