预应力模态分析
理论与物理
预应力模态分析是什么
老师，什么是“预应力模态分析”？
是存在初始应力（预应力）状态下的固有频率分析。结构处于受拉或受压状态时，固有频率会发生变化。这和拉紧琴弦音调会变高的原理相同。
就像吉他的弦。增加张力，振动频率就会上升。
正是如此。相反，如果施加压缩力，振动频率会下降。当压缩力接近临界屈曲载荷时，频率会趋近于零。频率为零 = 屈曲。
控制方程
包含预应力的特征值问题：
$[K_\sigma]$ 是几何刚度矩阵（应力刚度），与屈曲分析相同。
屈曲是 $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$，振动是 $([K_0] + [K_\sigma] - \omega^2 [M])\{\phi\} = \{0\}$。$[K_\sigma]$ 是共通的！
观察得非常完美。屈曲和振动共享相同的几何刚度矩阵。拉伸预应力使 $[K_\sigma] > 0$，增加整体刚度，从而提高频率。压缩预应力使 $[K_\sigma] < 0$，降低整体刚度，从而降低频率。
应用实例
旋转体的振动中，离心力引起的预应力很重要呢。
涡轮叶片和旋转圆盘中，离心力会产生拉伸预应力，使频率上升。这被称为旋转软化/硬化。由于每个转速下的频率都会变化，因此需要评估各速度下的固有频率。
FEM分析步骤
1. 静力分析（预加载） — 求解初始应力（压缩、拉伸、离心力等）
2. 几何刚度矩阵的构成 — 根据静力分析的应力计算 $[K_\sigma]$
3. 特征值分析 — 以 $[K_0] + [K_\sigma]$ 为刚度求解固有频率
和屈曲分析的步骤几乎一样呢。
总结
我来整理一下预应力模态分析。
要点：
- 初始应力会改变频率 — 拉伸使频率上升，压缩使频率下降
- $[K_0] + [K_\sigma] - \omega^2 [M] = 0$ — 加入了几何刚度的特征值问题
- 屈曲和振动共享相同的 $[K_\sigma]$ — 屈曲点处频率为零
- 旋转体中离心力预应力很重要 — 旋转软化/硬化
- 静力分析→特征值分析的两阶段 — 与屈曲分析步骤相同
“屈曲点处频率为零”这个关系很深刻啊。振动和屈曲被一个理论联系起来了。
这是结构力学中最美妙的关系之一。VCT（振动关联技术）就是利用这个关系，通过频率变化来无损预测屈曲载荷。
拉紧琴弦的高音与压缩弹簧的频率
对弦乐器的弦施加张力，固有频率会上升。相反，施加压缩载荷，固有频率会下降，在屈曲载荷（Pcr）处变为0。这个关系可以用 f²=f₀²(1-P/Pcr) 表示，可以通过测量的固有频率反算出压缩载荷P（难以实测的预应力量）。这个原理实际上被用于桥梁缆索的张力管理。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施加力因此加速度可忽略”的假设。冲击载荷或振动问题中绝对不能省略此项。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的性质就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形”，强度是“不易破坏”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用于整个内部的力”（体积力），轮胎压路面的力是“仅作用于表面的力”（表面力）。风压、水压、螺栓紧固力…全都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却成了“压缩”——听起来像笑话，但在3D空间中坐标系旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。弹一下吉他的弦试试。声音会一直持续吗？不，会逐渐变小。因为振动能量通过空气阻力和弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——故意吸收振动能量来提高乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置合适的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
求解器设置
请告诉我预应力模态分析的设置方法。
Nastran
```
SOL 103
CEND
SUBCASE 1
LOAD = 100 $ 预加载
METHOD = 10 $ 特征值分析
SPC = 1
```
在SOL 103中设置载荷，会自动进行预应力模态分析。
Abaqus
```
*STEP
*STATIC
*CLOAD
...
*END STEP
*STEP
*FREQUENCY
20, ,
*END STEP
```
在Static步骤后放置Frequency步骤。前一步骤的应力会自动反映到几何刚度中。
Ansys
```
/SOLU
ANTYPE, STATIC
PSTRES, ON ! 启用应力刚度
SOLVE
FINISH
/SOLU
ANTYPE, MODAL
MODOPT, LANB, 20
SOLVE
```
Ansys中PSTRES, ON是必须的对吧。和屈曲分析时一样。
忘记PSTRES, ON，应力刚度就不会被计算，会得到预应力影响为零的“普通”固有频率。和屈曲分析一样的陷阱。
预加载类型
还有温度引起的预应力啊。
结构在受约束状态下发生温度变化会产生热应力，这会影响频率。火灾时的钢梁，温度上升会产生轴向压缩力，导致频率下降→屈曲，就是这样一个过程。
总结
我来整理一下预应力模态的数值方法。
预应力模态的数值方法，整理如下。
なった
詳しく
報告