4. 圆柱壳结构模态仿真与周向波数识别
圆柱壳模态仿真：从振型到周向波数识别
两端简支圆柱壳的模态由轴向半波数 m 和周向波数 n 共同描述。仿真目标不是简单列出频率，而是从振型中识别出每一阶对应的 (m, n)，再绘制频率-波数关系图，用于判断结构在特定频段内是否存在密集模态或波型转换。
几何与材料参数（教学示例）：
- 半径 R = 0.15 m，长度 L = 0.30 m，壁厚 h = 0.002 m
- 弹性模量 E = 210 GPa，泊松比 ν = 0.30，密度 ρ = 7850 kg/m³
- 边界条件：两端简支（轴向位移约束，径向与周向自由）
- 求解频段：0–5000 Hz，提取前 30 阶模态
简支圆柱壳的周向波数 n 表现为振型沿圆周的完整波形数，轴向半波数 m 表现为沿母线的半波数。识别时先看径向位移云图：沿圆周数出波峰-波谷交替次数，即为 n；沿轴向数出半波个数，即为 m。对于 n = 0 的模态，振型呈轴对称呼吸或轴向伸缩，需单独标记。
识别流程分三步：提取模态 → 判读振型 → 归类波数。在有限元后处理中，逐阶查看径向位移云图，沿圆周路径提取位移数据并做 FFT，峰值个数即为 n；沿轴向路径提取位移数据，半波个数即为 m。对于 n = 0 的模态，圆周方向无波形，需通过轴向变形模式判断 m。
判读要点：当 n 相同而 m 不同时，频率随 m 增大而升高；当 m 相同而 n 增大时，频率先降后升，存在一个使频率最低的 n 值。频率-波数图上的曲线交叉点对应波型转换，是后续优化或避频设计的重点关注区域。
| 模态阶次 | 频率 (Hz) | 周向波数 n | 轴向半波数 m | 振型特征 |
|---|---|---|---|---|
| 1 | 842 | 2 | 1 | 周向2波，轴向1半波 |
| 2 | 1105 | 3 | 1 | 周向3波，轴向1半波 |
| 3 | 1560 | 2 | 2 | 周向2波，轴向2半波 |
| 4 | 1890 | 4 | 1 | 周向4波，轴向1半波 |
| 5 | 2150 | 3 | 2 | 周向3波，轴向2半波 |
| 6 | 2480 | 0 | 1 | 轴对称呼吸模态 |
上表为教学示例数据，用于说明归类方法。实际仿真中，频率值由求解器输出，波数由振型判读得到。将 (n, f) 数据点按 m 分组绘制，即可得到频率-波数关系图。若某阶模态的振型在圆周方向出现非整数波形或局部变形，说明网格密度不足或边界条件施加有误，需检查周向单元数量是否满足每波长至少 6–8 个单元。
常见问题：周向波数识别错误多源于网格过粗导致波形混叠，或后处理路径未沿真实圆周方向提取。建议在圆周方向至少布置 72 个单元（每 5° 一个），并在提取位移时使用柱坐标系下的径向分量，避免直角坐标分量混入轴向或周向变形。