NAFEMS FV32: 厚肉圆板的自由振动

很好的问题。FV32是**厚肉圆板的自由振动**基准问题。具体来说，就是通过FEM分析求圆板的面外弯曲固有频率，然后与理论解 $f_1 = 1.4568\,\text{Hz}$ 进行比较的问题。

这正是核心所在。厚板中**剪切变形的影响**变得很大。薄板理论用Kirchhoff理论，假设"变形后截面仍然垂直于中性层"，即剪切变形为零。但对于厚板，这个假设不成立。用Kirchhoff理论求解会**高估固有频率**。

在FV32的条件下（厚度/半径 = 0.1），用Kirchhoff理论得到的一次固有频率比Mindlin理论高出数个百分点。实际工作中，如果"FEM结果比理论值高"，第一时间要怀疑是否考虑了剪切变形。这正是FV32被广泛用来**验证Mindlin板单元和实体单元精度**的原因所在。

虽然看起来很简单，但设置相当巧妙。

没错。真正重要的是无量纲比例 $t/R = 0.1$。这个值处于"明显是厚板，但又不至于极厚"的区域，正好能突出薄板理论与厚板理论之间的差异。比如汽车制动盘的厚度比也大约在 $t/R = 0.05\sim0.15$ 范围内，所以从工程应用角度讲也很有代表性。

好。首先Kirchhoff（古典薄板理论）的基本假设是**"变形后截面仍垂直于中性层"**，也就是说剪切变形为零。这样断面的转角 $\theta$ 就由面外位移 $w$ 的梯度唯一确定：

完全正确。你拿一本厚书弯曲时，页面（截面）就不再垂直于书脊了，对吧？这就是剪切变形的效果。薄纸张时截面始终垂直，但厚度增加后就管不了。

正是如此。刚度过大 → 振动频率升高。这在FEM实务中非常关键，选板或壳单元时**必须确认是Mindlin类型**。

Mindlin板在圆形坐标系下的自由振动方程，对轴对称模态（周向波数 $n = 0$），面外位移 $w(r,t)$ 和断面转角 $\psi(r,t)$ 满足联立方程：

好眼光。Mindlin理论假设厚度方向剪切应力分布均匀，但实际是抛物线分布。$\kappa$ 就是用来修正这个差异的。对矩形截面，理论值恰好是 $5/6$。圆形梁有时用 $\kappa = 6/7$，但板的惯例是 $5/6$。

NAFEMS公布的一次固有频率参考值是：

Kirchhoff薄板理论的简支圆板一次固有频率为：

自由振动，外力为零。假设位移为 $\{u\} = \{\phi\} e^{i\omega t}$，就归结为一般化特征值问题：

一致质量一般精度更高。集中质量虽然计算轻，但对FV32这样的规模，推荐用一致质量。两者的差别通常不超过0.5%，大规模模型（数百万DOF）内存吃紧才用集中质量。

两种都能得到正确答案，但要点各不相同。

绝对是。FV32的核心就是剪切变形，所以锁定会致命。完全积分1阶单元（CQUAD4、C3D8完全积分）会让固有频率远超理论值。用低减积分或直接选2阶单元是标准做法。

看用Mindlin板单元（CQUAD8等级）的典型收敛数据：

就是这样。振动解析中2阶单元的性价比压倒性优势。实务的黄金法则："先用2阶单元和粗网格 → 验证收敛性 → 必要时细化"。

好，逐步来。

最常见的失误两个。第一是**单位制混用**。比如 $E$ 用GPa，密度却用 $\text{tonne/mm}^3$ 而不是 $\text{kg/m}^3$，会导致频率差30倍。频率正比 $\sqrt{E/\rho}$，密度错一千倍，频率就错30倍。

第二是**边界条件的误解**。"简支"在板单元上是只约束面外位移 $w = 0$，转动自由。但容易把面内位移也锁住，变成"固定支持"，这样固有频率就会偏高。实体单元也是，只锁z方向位移，不要管r和θ方向。

完全同意。要只求轴对称模态（$(0,1)$ 模式），用**2维轴对称单元**（CAX8、PLANE183轴对称等）建半径截面就够了。DOF会剧减。

四大求解器的设置对比：

Lanczos是FV32这样中小规模问题的首选，鲁棒性最好。超大规模（百万DOF以上）才考虑AMLS或Shift-Invert这样的高级方法。但对FV32，再好的算法也要不了几秒。

用对了单元的话，各求解器都能到参考值 $f_1 = 1.4568\,\text{Hz}$ 的0.1%以内误差。这正是"基准问题"存在的意义。

问题设置"干净"就是这样。真正的差异在你的网格质量、单元选择、边界条件设置上。"求解器间差异"还不如"使用者设置差异"大。这个道理在所有CAE工作中都适用。

有的。Reddy（1984）的**三阶剪切变形理论（TSDT）**就是。Mindlin理论假设厚度方向剪切应变均匀，TSDT用三次多项式来近似，这样就不需要剪切修正系数 $\kappa$ 了。

$t/R > 0.15$ 左右的超厚板，或者复合材料（CFRP等）层间剪切刚度差异大的情况。复合材料里Mindlin的均匀剪切假设就不适用了。但FV32本身Mindlin足够，这是对"更深层"的认知。

圆板这样的光滑形状用IGA特别合适。基于NURBS的基函数不受网格近似的制约，边界拟合精度天生高。研究级报告显示，IGA在同等DOF下精度比传统FEM更高。但商用软件中IGA支持还有限，只有LS-DYNA、COMSOL等少数几个。

1.50Hz比参考值高3%左右。典型的"用了Kirchhoff单元"症状。什么单元？

就是了。CTRIA3是薄板理论，不考虑剪切变形，而且1阶单元网格敏感性也高。换成CQUAD8重做。

FEM调试就是这样。知道"怎样变形代表什么问题"就能快速排查。高频 → 刚度过大；低频 → 刚度不足或约束不足。FV32这样简洁的问题最适合掌握这个本领。实务的复杂模型才能派上大用场。

必须的。FV32的一次模态是**轴对称面外弯曲**，像碗的形状。如果出来是面内振动或非对称的，那边界条件肯定有问题。

完全正确。"结果一定要后处理可视化"是CAE铁律。FV32正因为简单，模态形状的异常才显而易见。这个习惯在复杂模型里才最关键。