NAFEMS FV32：厚圆盘自由振动

问得好。FV32是一个处理厚圆盘自由振动的基准测试。具体来说，是用有限元法求解面外弯曲的固有频率，并与理论解 $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$ 的梯度唯一确定：

另一方面，Mindlin（Reissner-Mindlin理论）将截面转角作为独立变量处理，允许剪切应变 $\gamma$ 存在：

正是这个意思。试着弯曲一本厚书，书页（截面）会无法保持与封面垂直而倾斜，对吧？那就是剪切变形的效果。如果是薄薄的一张纸，截面总是垂直的，但变厚后就无法忽略了。

理解得非常完美。刚度被高估 → 频率偏高。这也是FEM实际应用中非常重要的一个点，选择板/壳单元时，一定要确认“是否是Mindlin类型”。

将Mindlin板的自由振动用圆盘坐标表示，对于轴对称模态（周向波数 $n = 0$），会得到面外位移 $w(r,t)$ 和截面转角 $\psi(r,t)$ 的联立方程：

其中：

• $D = \dfrac{Et^3}{12(1-\nu^2)}$：板的弯曲刚度
• $G = \dfrac{E}{2(1+\nu)}$：剪切弹性模量
• $\kappa = 5/6$：Mindlin剪切修正系数

问得好。Mindlin理论假设板厚方向的剪切应力分布是均匀的，但实际上是抛物线分布。$\kappa$ 就是用来修正这个差异的，对于矩形截面，严格来说是 $5/6$。顺便提一下，对于圆形截面的梁，有时也使用 $\kappa = 6/7$，但对于板，习惯上使用 $5/6$。

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

这对应于面外弯曲的轴对称一次模态（$(0,1)$ 模态，无节线）。它是基于Mindlin厚板理论的解析解推导出来的。

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

其中 $\lambda_{01} \approx 4.935$（简支的一次根）。用这个公式计算，大约在 $1.50\,\text{Hz}$ 前后，比Mindlin理论的 $1.4568\,\text{Hz}$ 高出约3%。可能觉得只有3%，但在实际工作中，无法区分“是3%的误差，还是3%的理论差异”是致命的。

因为是自由振动，所以外力为零。假设位移为 $\{u\} = \{\phi\} e^{i\omega t}$，则归结为广义特征值问题：

其中 $[K]$ 是整体刚度矩阵，$[M]$ 是整体质量矩阵，$\omega = 2\pi f$ 是角频率，$\{\phi\}$ 是模态向量。

问得很敏锐。一般来说，一致质量矩阵精度更高。但集中质量矩阵是对角化的，计算量小。对于FV32这种规模的问题，推荐使用一致质量。如果是大规模模型（数百万自由度）且内存紧张，使用集中质量也能获得实用的精度。实际差异通常在0.5%以内。

两种都能得到正确答案，但各自的注意事项不同。

没错。FV32正是剪切变形唱主角的问题，如果发生自锁，固有频率会比理论值高出很多。完全积分的一次单元（CQUAD4完全积分、C3D8完全积分）尤其要注意。使用减缩积分或B-bar法，或者一开始就选择二次单元是最佳实践。

我们来看看使用Mindlin板单元（相当于CQUAD8）时的典型收敛数据。

是的。在振动分析中，二次单元的性价比极高。这是FV32明确展示的教训。在实际工作中，“先用二次单元进行少量网格划分 → 确认收敛 → 必要时细化”也是黄金法则。