跳到主要内容
Sandbox Physics

M067 · 振动设计

调谐质量阻尼器设计

给主结构添加滑动质量、弹簧和阻尼器。直接拖动结构释放,同步检查两块质量的运动,将解析力响应曲线与三个周期的采集测量比较,并在相对行程约束下搜索调谐与阻尼。

交互模型调谐质量阻尼器设计
已采集时间0 s0\,\mathrm{s}
主质量位移0 m0\,\mathrm{m}
主质量速度0 m s−10\,\mathrm{m\,s^{-1}}
机械能0 J0\,\mathrm{J}
积分有向执行器做功0 J0\,\mathrm{J}
独立黏性耗散0 J0\,\mathrm{J}
独立能量收支误差0 J0\,\mathrm{J}
步长减半位移差 · 最初四个周期0 m0\,\mathrm{m}
实际积分步长0 ms0\,\mathrm{ms}
模型范围与记录状态Ready\text{Ready}
吸振质量绝对位移00
有向相对行程00
解析稳态主振幅00
解析稳态相对行程振幅00
裸结构稳态振幅00
三个周期采集的主振幅00
采集谐波残差均方根00
频段采样最大主振幅00
频段采样最大相对行程00
约束搜索结果00

物理教程

保护结构,同时预算行程

背景次级质量通过弹簧和阻尼器相对主结构运动,其反作用改变主响应并形成两个共振。

为什么值得理解MIT Engineering Dynamics 第 26 讲推导双自由度谐波响应并演示动态吸振器。本工作台加入有限采集窗与行程约束下的有界网格搜索。

先抓重点

核心问题
压低一个共振,是否会产生更严重的邻近峰值或过大行程?
一句话直觉
同时比较整个显示频段与相对行程,不能只看一个频率下的主振幅。

核心数学模型

两个耦合受力方程

Mx¨+Cx˙+Kx+c(x˙−y˙)+k(x−y)=F0cos⁡ωt,my¨+c(y˙−x˙)+k(y−x)=0M\ddot x+C\dot x+Kx+c(\dot x-\dot y)+k(x-y)=F_0\cos\omega t,\quad m\ddot y+c(\dot y-\dot x)+k(y-x)=0

两个位置均为绝对坐标;连接弹簧与阻尼器使用相对运动。

调谐与质量比

ω0=K/M,ωa=k/m,μ=m/M\omega_0=\sqrt{K/M},\quad \omega_a=\sqrt{k/m},\quad \mu=m/M

辅助频率属于独立吸振器,耦合后实际模态会改变。

复动态刚度

[K+k−Mω2+iω(C+c)−k−iωc−k−iωck−mω2+iωc][XY]=[F00]\begin{bmatrix}K+k-M\omega^2+i\omega(C+c)&-k-i\omega c\\-k-i\omega c&k-m\omega^2+i\omega c\end{bmatrix}\begin{bmatrix}X\\Y\end{bmatrix}=\begin{bmatrix}F_0\\0\end{bmatrix}

求解矩阵得到独立的无限时间谐波参考。

独立能量账本

E=12Mx˙2+12my˙2+12Kx2+12k(y−x)2,E−E0=W−QE=\tfrac12M\dot x^2+\tfrac12m\dot y^2+\tfrac12Kx^2+\tfrac12k(y-x)^2,\quad E-E_0=W-Q

耗散分别积分两个阻尼器的功率,不由能量差定义。

采样约束设计

min⁡ra,ζamax⁡ω∈G∣X(ω)∣,max⁡ω∈G∣Y(ω)−X(ω)∣≤smax⁡\min_{r_a,\zeta_a}\max_{\omega\in\mathcal G}|X(\omega)|,\quad \max_{\omega\in\mathcal G}|Y(\omega)-X(\omega)|\le s_{\max}

有限频率与参数网格定义搜索,网格间的极值可能更大。

常见难点

调谐不会消除两个模态

容易误解完美凹口能保证所有频率安全。

正确理解凹口位于分裂共振之间,需要检查两个峰和行程曲线。

坐标与行程不同

容易误解吸振质量的位置就是机械行程。

正确理解行程相对于运动的主结构。

网格设计有边界

容易误解搜索找到全局最优阻尼器。

正确理解它检验固定有限频段、网格以及谐波载荷约定。

动手实验

  1. 01

    比较调谐

    采集调谐与失谐记录。

    应该观察到: 两块质量的运动和共振峰一同变化。
  2. 02

    定位反共振

    选择理想反共振并采集完整记录。

    应该观察到: 解析主振幅凹口可与衰减的有限时间瞬态同时存在。
  3. 03

    约束吸振器

    降低允许行程并运行搜索。

    应该观察到: 可行网格解可能改变或消失。
  4. 04

    核验计算

    使用保守能量交换,再减半积分步长。

    应该观察到: 能量交换仍然存在,账本误差随细化减小。