手算核心方法（三）：邓克利法估算多自由度系统基频

各位工程师朋友，欢迎来到第四章。

前两章我们聊了瑞利法和矩阵迭代法。说实话，瑞利法对振型假设要求高，矩阵迭代法虽然准但计算量大。那有没有一种方法，既简单，又能快速给个大概的基频？

有。就是邓克利法。

我个人习惯叫它「懒人估算法」。你不需要假设振型，也不需要反复迭代。只需要知道每个质量单独存在时的频率，就能把整个系统的基频估出来。

听起来有点神奇？我们一步步拆开看。

4.1 邓克利法的核心思想

邓克利法（Dunkerley's Method）是一种经验性的近似方法。它的核心逻辑很简单：

系统的基频平方的倒数，约等于各部件单独频率平方的倒数之和。

说白了，就是把多自由度问题，拆成多个单自由度问题。每个单自由度问题都好算，加起来就完事了。

你想想看，这比瑞利法省事多了吧？不用猜振型，不用积分，甚至不用列完整的刚度矩阵。

4.2 邓克利法的推导逻辑

虽然我们不用背推导过程，但理解它的物理意义很重要。

邓克利法的本质，是忽略了高阶模态对基频的「刚度贡献」。它假设系统基频主要由各质量的「柔度叠加」决定。

具体来说：

1. 对于一个多自由度系统，它的柔度矩阵可以写成各阶模态贡献的叠加。
2. 基频对应的柔度最大，高阶模态的贡献相对较小。
3. 邓克利法只保留各质量单独作用时的柔度项，忽略交叉耦合项和高阶项。

所以，邓克利法估算的基频，总是比真实值偏低。这是它的一个特点，也是它的一个局限。

4.3 手算步骤与示例

我们用一个三自由度弹簧-质量系统来演示。

系统参数：

参数	符号	数值
质量1	m₁	2 kg
质量2	m₂	3 kg
质量3	m₃	1 kg
弹簧刚度1	k₁	1000 N/m
弹簧刚度2	k₂	1500 N/m
弹簧刚度3	k₃	800 N/m

步骤一：计算各质量单独存在时的频率

假设只有 m₁ 存在，其他质量固定不动。此时系统的等效刚度是多少？

嗯，这里要注意：邓克利法中的「单独存在」，是指该质量在系统中独立振动，其他质量视为固定。

对于 m₁：等效刚度 k₁_eq = k₁ = 1000 N/m

ω₁₁ = √(k₁_eq / m₁) = √(1000 / 2) = √500 ≈ 22.36 rad/s

对于 m₂：等效刚度 k₂_eq = k₂ = 1500 N/m

ω₂₂ = √(1500 / 3) = √500 ≈ 22.36 rad/s

对于 m₃：等效刚度 k₃_eq = k₃ = 800 N/m

ω₃₃ = √(800 / 1) = √800 ≈ 28.28 rad/s

步骤二：代入邓克利公式

1 / ω₁² ≈ 1 / 22.36² + 1 / 22.36² + 1 / 28.28²

1 / ω₁² ≈ 1/500 + 1/500 + 1/800

1 / ω₁² ≈ 0.002 + 0.002 + 0.00125 = 0.00525

ω₁ ≈ √(1 / 0.00525) ≈ √190.48 ≈ 13.80 rad/s

步骤三：对比精确解

这个系统的精确基频（用矩阵迭代法算出来）大约是 14.52 rad/s。

邓克利法估算值 13.80 rad/s，误差约 5%。

对于手算来说，这个精度已经相当不错了。

4.4 邓克利法的适用场景与局限

我总结了一下，邓克利法最适合以下场景：

• 方案设计阶段：快速判断结构是否可能发生共振
• 质量分布较分散的系统：比如多跨梁、多盘轴系统
• 高阶模态对基频影响小：即系统刚度主要由各部件自身刚度决定

但遇到这些情况，我建议你换方法：

• 系统有强耦合：比如相邻质量之间的弹簧刚度很大
• 需要精确值：比如做疲劳分析或共振避让
• 系统有密集模态：邓克利法只给一个值，无法区分

4.5 知识结构图

下面这张图，帮你理清邓克利法的知识脉络：

4.6 邓克利法与瑞利法的对比

很多学员会问：邓克利法和瑞利法，到底该用哪个？

我的建议是：看你的需求。

对比项	邓克利法	瑞利法
估算方向	下限（偏低）	上限（偏高）
是否需要振型假设	不需要	需要
计算量	极小	中等
适用系统	弱耦合、质量分散	振型可合理假设
精度	5%~15%	1%~5%（振型假设好）

你看，这两个方法其实是互补的。邓克利法给你一个下限，瑞利法给你一个上限。如果你两个都算，那真实基频就在这两个值之间。这在工程上非常有用——你不需要精确值，只需要知道一个范围。