模态阻尼识别与设置
模态阻尼识别与设置的理论基础
什么是阻尼
教授,在结构力学中如何处理"阻尼"?
阻尼是振动能量散逸的效应。如果没有阻尼,结构将永远振动不止。实际结构总是存在阻尼,它使得共振时的振幅保持有限。
运动方程
带阻尼的运动方程:
$[C]$ 是阻尼矩阵。$[M]$ 和 $[K]$ 在物理上很清楚,但$[C]$ 通常有很大的不确定性。
确定阻尼矩阵是问题所在吧。
没错。阻尼的建模是FEM动力学分析中最不确定的参数。
阻尼模型
主要的阻尼模型:
1. 模态阻尼(Modal Damping)
为每个模态直接指定阻尼比 $\zeta_i$。这是最通用的方法。
$q_i$ 是模态坐标,$\omega_i$ 是固有圆频率。
2. Rayleigh阻尼(Rayleigh Damping)
$[C] = \alpha [M] + \beta [K]$。$\alpha$ 和 $\beta$ 通过在两个振动频率处匹配 $\zeta$ 来确定。
3. 结构阻尼(Structural Damping)
也称为滞后阻尼。与频率无关的阻尼,用复刚度 $[K^*] = [K](1 + ig)$ 表示。$g$ 是结构阻尼系数。
三种模型怎么区分使用?
阻尼比的典型值
| 结构 | 阻尼比 $\zeta$ |
|---|---|
| 钢结构(焊接) | 0.5~1% |
| 钢结构(螺栓连接) | 1~2% |
| 混凝土结构 | 3~5% |
| 隔震结构 | 10~30% |
| 机械结构 | 1~3% |
| 复合材料结构 | 0.5~2% |
钢结构中0.5~1%...非常小啊。
钢的内部阻尼很小。因此钢结构在共振时容易产生较大振幅,阻尼的设置对结果的影响很大。
总结
模态阻尼的理论要点整理一下。
要点:
- 阻尼是动力学分析中最不确定的参数 — 需要进行敏感性分析
- 三种模型 — 模态阻尼、Rayleigh阻尼、结构阻尼
- 模态阻尼最为通用 — 为每个模态指定 $\zeta_i$
- Rayleigh阻尼用于时间域 — 用 $\alpha, \beta$ 在两个频率处匹配
- 阻尼比的典型值 — 钢:1%、混凝土:5%。根据用途而异
阻尼设置会导致结果相差数倍。共振时的振幅与 $1/(2\zeta)$ 成正比,所以 $\zeta = 1\%$ 和 $\zeta = 2\%$ 的振幅相差2倍。
因此阻尼是"影响最大且最不确定的参数"。没有对阻尼进行敏感性分析就不应该相信动力学分析的结果。
"2%阻尼"这个"魔法数字"
结构阻尼比ζ=2%作为设计约定俗成的标准广泛应用,但实际钢结构物阻尼在0.5~5%之间差异很大。这个值是Lankford在1954年从建筑物实测数据中统计提出的中位数。从1970年代被美国建筑基准法(UBC)采用后便作为"标准值"沿用至今,但焊接结构不足1%,螺栓连接则为3~5%,完全不同,不能忽视。
模态阻尼识别与设置的数值计算方法
阻尼的设置方法
在FEM中怎样设置阻尼?
Nastran中的模态阻尼
```
TABDMP1, 1, CRIT
, 0., 0.02, 100., 0.02, ENDT
```
为所有模态设置 $\zeta = 2\%$。也可为不同频率范围设置不同的 $\zeta$。
Abaqus中的模态阻尼
```
*MODAL DAMPING
1, 50, 0.02
```
为模态1~50一次设置 $\zeta = 2\%$。
Rayleigh阻尼的设置
确定 $\alpha$ 和 $\beta$ 的方法:
当两个振动频率 $f_1, f_2$ 处 $\zeta_1 = \zeta_2 = \zeta$ 时:
$f_1$ 和 $f_2$ 怎样选择?
选择关注的振动频率范围的下限和上限。例如地震响应分析针对1~10 Hz,则选 $f_1 = 1$ Hz, $f_2 = 10$ Hz。超出这个范围的阻尼比会有偏差,需要注意。
Abaqus中的Rayleigh阻尼
```
*DAMPING, ALPHA=0.5, BETA=0.001
```
阻尼的识别
实际结构的阻尼比怎样测量?
通过实验模态分析来测量:
1. 锤击激励法 — 用冲击锤激励,采集加速度
2. 激励器法 — 用正弦波/随机激励,获取频率响应函数(FRF)
3. 阻尼识别 — 从FRF的半幂带宽法,或曲线拟合
半幂带宽法:从共振峰的振幅为 $1/\sqrt{2}$ 的两个频率 $f_1, f_2$ 计算:
可以简单地测量啊。
原理很简单,但试验精度(激励点、采集点、噪声处理)会影响结果。用多种方法进行交叉验证比较安全。
总结
阻尼设置的数值方法整理一下。
要点:
- TABDMP1(Nastran), *MODAL DAMPING(Abaqus) — 模态阻尼的设置
- $\alpha, \beta$ 的确定 — 在两个频率处匹配
- 通过实验模态分析测量阻尼比 — 半幂带宽法是基础
- 注意Rayleigh阻尼的频率范围 — 超出范围会偏差
- 没有实验数据时,用文献值并进行敏感性分析
半幂带宽法和对数衰减率的实测步骤
从实验获得阻尼比的半幂带宽法中,对频率响应函数(FRF)的共振峰,从峰值频率fr处±3dB点的间隔Δf计算ζ=Δf/(2fr),精度取决于频率分辨率,需要不小于Δf的五分之一的分辨率。冲击锤试验中,从自由衰减波形的包络线对数求对数衰减率来直接求解阻尼比在某些情况下精度更高。
模态阻尼识别与设置的实务应用
阻尼的实务应用
在实务中阻尼怎样设置?
设计规范中的阻尼值
在某些设计规范中规定了阻尼比:
| 规范 | 结构 | 阻尼比 $\zeta$ |
|---|---|---|
| 建筑基准法(日本) | 混凝土结构 | 5% |
| 建筑基准法(日本) | 钢结构 | 2% |
| 欧洲规范8 | 混凝土结构 | 5% |
| ASCE 7 | 全结构 | 5%(默认) |
| API 617 | 旋转机械 | 实测或1% |
地震响应中混凝土结构的5%是默认值啊。
混凝土的5%代表"包含裂纹的混凝土结构的等效阻尼"。无裂纹的混凝土结构阻尼仅1~2%,但地震中出现微小裂纹后能量散逸增加,等效阻尼达到5%。
没有实验数据时
实务检查清单
请给出阻尼设置的检查清单。
如果把阻尼比改为 $\zeta/2$,共振振幅就变成2倍...敏感性分析是必须的啊。
阻尼是FEM动力学分析的最大不确定因素。忽视敏感性分析,结果就失去意义。
隔震建筑的阻尼器优化
超高层建筑的减振设计中,粘性阻尼器或摩擦阻尼器转换为模态阻尼比形式,来实现目标阻尼比(通常3~5%)的阻尼器配置最优化。六本木Hills森塔(238m)在安装粘性减振阻尼器后,低阶模态的阻尼比从1%上升到4%,地震应答加速度降低40%。
模态阻尼识别与设置的软件比较
阻尼设置工具
在求解器间的阻尼设置有什么差异?
可以为不同材料设置不同的阻尼吗?
Nastran的MAT1卡片的GE字段可以为每种材料设置结构阻尼系数。Abaqus和Ansys也有类似功能。这对钢和橡胶混合结构中设置不同阻尼很有用。
实验模态分析工具
| 工具 | 特点 |
|---|---|
| LMS Test.Lab | 从FRF识别阻尼。与FEM的相关性 |
| ME'scope | 模态形状与阻尼的可视化 |
| ModalVIEW | Excel基础。简便的模态分析 |
选型指南
根据分析方法能用的阻尼模型不同啊。
模态阻尼用于模态重叠法,Rayleigh阻尼用于直接积分,结构阻尼用于频率响应。用错了就会得到物理上没有意义的结果。
MSC Nastran模态阻尼输入形式
Nastran SOL 108(频率响应)的TABDMP1卡片可以定义模态阻尼比的频率依赖性。用CRITDAMP=MODE选项为每个模态单独设置阻尼比,可以直接输入实测值。LMS(现Siemens Simcenter)公司利用这一功能建立了试验-分析相关流程,试验后1小时内就能更新FEM的阻尼模型,为汽车行业提供此系统。
模态阻尼识别与设置的先进研究
阻尼的先进研究
请讲讲阻尼模型的最新进展。
非粘性阻尼(Non-viscous Damping)
Rayleigh阻尼假设粘性阻尼($[C]\{\dot{u}\}$,与速度成正比),但实际结构的阻尼通常是非粘性的。摩擦阻尼(接合部微小滑动)、滞后阻尼(材料非线性历程)无法用粘性模型准确表示。
广义阻尼模型:
$[G(t)]$ 是阻尼核函数。粘性阻尼是 $[G(t)] = [C]\delta(t)$ 的特殊情形。
接合部阻尼
螺栓接合部的微小滑动(微动)是结构整体阻尼的主要来源。接合部阻尼模型(Iwan单元、Jenkins单元)正在被研究。
接合部微小滑动是阻尼的主要原因?
焊接的钢结构中,材料本身的内部阻尼非常小($\zeta < 0.1\%$)。实测的1~2%的阻尼主要来自接合部的微动和非结构构件(涂层、电缆等)的摩擦。
能动阻尼(Active Damping)
用压电素子或MR阻尼器能动控制阻尼。在FEM中可以耦合传感器-致动器的反馈控制与结构分析。
总结
阻尼先进研究要点整理。
阻尼是结构力学的"最柔和话题",却是振动控制实效性的最关键参数。
材料损失因子与模态阻尼的关系
整体结构的模态阻尼比ζn由各部件的材料损失因子η与部件在该模态下储存的应变能的比率决定(模态应变能法)。在高阻尼橡胶(η=0.3~1.0)中少量夹入就能大幅提高整体结构阻尼比,原理就是这样。制振钢板(η=0.1~0.3)用在汽车地板上可以降低室内噪声3~6dB。
模态阻尼识别与设置的故障排查
阻尼设置故障
阻尼设置中常见的问题是?
共振振幅不切实际
共振处的响应非常大。
阻尼太小了。确认阻尼值。$\zeta = 0$(无阻尼)时共振振幅无限大。
检查:
- 是否设置了阻尼(最多见的是忘记设置)
- 阻尼比值是否合适(与文献值、规范值对比)
- $\zeta = 0.01$ 时振幅为 $1/(2 \times 0.01) = 50$ 倍。这是否合理
Rayleigh阻尼中低频/高频响应异常
Rayleigh阻尼中 $\alpha$ 项在低频占主导(阻尼大),$\beta$ 项在高频占主导(阻尼大)。指定的两个频率间 $\zeta$ 合适,但范围外阻尼会过大。
怎样处理?
在时间域中使用了结构阻尼
在时间域应该用模态阻尼或Rayleigh阻尼吧。
对。频率域 = 结构阻尼OK。时间域 = Rayleigh或模态阻尼。用错了会因果律破坏,导致非物理的响应。
总结
阻尼设置故障排查要点整理。
模态阻尼比被低估的情况
FEM分析的模态阻尼比通常比实测偏低。这是因为FEM不能表现接合部(螺栓、铆钉、焊缝)处摩擦和微动造成的能量散逸。接合部阻尼可以用集中阻尼单元(C单元)来补充,或者直接把实测阻尼比输入模态分析,这是实务的实用做法。
更详细
错误