3. 多自由度系统:从单点到全局的振动视角

各位工程师朋友,欢迎来到第三章。

上一章我们聊了单自由度系统,说白了就是一个质量块加一根弹簧。但现实中的多旋翼,哪有这么简单?机臂、机身、电机座、飞控板……每个部件都在振动,它们之间还会互相影响。

这就是多自由度系统要解决的问题。

3.1 多自由度系统建模方法

建模这事儿,我个人的习惯是「先简化,再精确」。你不可能把每个螺丝都建进去,那算到明年也算不完。

3.1.1 牛顿-欧拉法

这是最直接的方法。对每个质量块写牛顿第二定律,然后联立方程组。

举个例子,一个两自由度系统:

m1·x1'' + (c1+c2)·x1' - c2·x2' + (k1+k2)·x1 - k2·x2 = F1
m2·x2'' - c2·x1' + c2·x2' - k2·x1 + k2·x2 = F2

写成矩阵形式,就清爽多了:

[M]{x''} + [C]{x'} + [K]{x} = {F}

其中[M]是质量矩阵,[C]是阻尼矩阵,[K]是刚度矩阵。

我的经验: 写矩阵的时候,注意对称性。质量矩阵和刚度矩阵通常是对称的,如果不对称,八成是你写错了。

3.1.2 拉格朗日法

这个方法我更喜欢,尤其系统复杂的时候。你只需要写出系统的动能T、势能U,然后套公式:

d/dt(∂L/∂q̇) - ∂L/∂q = Q

其中L = T - U,Q是广义力。

为什么推荐这个方法?因为你不必纠结每个力的方向,能量是标量,好处理得多。

避坑指南: 我曾经在建模一个六轴无人机时,用牛顿法推导了半天,结果发现有个力的方向搞反了。换成拉格朗日法,半小时搞定。所以,复杂系统我建议用拉格朗日法。

3.2 模态分析基础

模态分析,说白了就是找系统的「固有振动模式」。每个模式都有两个关键参数:模态频率和模态振型。

3.2.1 模态振型

模态振型描述的是:在某个特定频率下,系统各部分的相对运动关系。

举个例子,一个两自由度系统:

  • 第一阶模态:两个质量块同向运动,像两个人一起荡秋千
  • 第二阶模态:两个质量块反向运动,像在玩跷跷板

你想想看,如果电机激励频率刚好接近第一阶模态频率,整个机臂就会一起晃。如果接近第二阶,机臂两端就会扭着动。

3.2.2 模态频率

模态频率就是系统「喜欢」振动的频率。求解方法很简单,解特征值问题:

([K] - ω²[M]){φ} = {0}

其中ω是模态频率(圆频率),{φ}是模态振型向量。

注意: 模态频率只取决于质量和刚度分布,跟阻尼无关。阻尼只影响振动幅值,不影响频率。

3.3 模态叠加法

这个方法的思路很巧妙:既然系统有N个模态,那任何振动都可以看成这些模态的线性组合。

数学表达就是:

{x(t)} = [Φ]{q(t)}

其中[Φ]是模态矩阵(每列是一个模态振型),{q(t)}是模态坐标。

这样做的好处是什么?

  • 把N个耦合的方程,变成N个独立的单自由度方程
  • 每个方程只对应一个模态,计算量大大降低
  • 你可以只保留前几阶模态,忽略高频模态
我的建议: 实际工程中,保留前3-5阶模态通常就够了。高阶模态贡献很小,忽略它们不会影响精度。

3.4 瑞利阻尼

阻尼这东西,说实话很难精确测量。瑞利阻尼提供了一个实用的近似方法。

它的形式很简单:

[C] = α[M] + β[K]

其中α和β是两个系数,需要根据实验数据确定。

怎么确定α和β?

  1. 测量两个不同模态的阻尼比ξ₁和ξ₂
  2. 解方程组:
ξ₁ = α/(2ω₁) + βω₁/2
ξ₂ = α/(2ω₂) + βω₂/2

解出α和β,就得到了整个系统的阻尼矩阵。

避坑指南: 我曾经在测试一款四轴时,发现仿真和实测对不上。查了半天,原来是瑞利阻尼的β系数设大了,导致高频模态被过度阻尼。记住,β主要影响高频,α主要影响低频。

3.5 实际应用中的注意事项

好了,理论讲完了,说点实际的。

问题 原因 解决方法
模态频率计算不准 边界条件简化过度 用实验模态分析修正模型
共振频率接近激励频率 结构设计不合理 调整刚度或增加阻尼
模态耦合严重 模态频率太接近 改变结构布局,拉开频率间距

最后说一句:多自由度系统分析,工具只是辅助,关键是你对物理本质的理解。我见过太多人拿着ANSYS算了一堆漂亮的模态图,结果连基本的一阶弯曲频率都解释不清楚。

嗯,这一章就到这里。下一章我们聊聊如何把这些理论用到实际的减震结构设计中去。