预应力模态分析
预应力模态分析的理论基础
预应力模态分析是什么
老师，"预应力模态分析"是什么意思？
初始应力（预应力）存在条件下的固有振动数解析。构造受到拉力或压缩力时，固有振动数会变化。弦张紧时声音变高，原理相同。
就像吉他的弦一样。拉力越大振动数越高。
完全正确。反之，压缩力作用时振动数会下降。当压缩力接近临界座屈荷重时，振动数接近零。振动数为零 = 座屈。
支配方程
含预应力的固有值问题：
$[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）处变为零。这个关系由 f²=f₀²(1-P/Pcr) 表示，从测量固有振动数可反算压缩荷重P（直接测量困难的预应力量）。这个原理在桥梁电缆张力管理中已实用化。
预应力模态分析的数值计算方法
求解器设置
请说明预应力模态分析的设置。
Nastran
```
SOL 103
CEND
SUBCASE 1
LOAD = 100 $ 预加载
METHOD = 10 $ 固有值分析
SPC = 1
```
在SOL 103中设置荷重时，自动变成预应力模态分析。
Abaqus
```
*步骤
*STATIC
*CLOAD
...
*END 步骤
*步骤
*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就不会计算应力刚度，预应力影响变为零，得到"普通的"无预应力固有振动数。这与座屈是同样的陷阱。
预加载的类型
温度产生的预应力也会影响。
结构在约束下温度变化时产生热应力，影响振动数。火灾中钢梁因温度升高产生轴向压缩力，导致振动数下降→座屈的过程。
总结
我整理预应力模态分析的数值方法。
要点：
- 静解析→固有值分析两步骤 — 与座屈相同步骤
- Nastran SOL 103只需添加荷重 — 最简单
- Abaqus: Static→Frequency — 自动反映预应力
- Ansys: 不要忘PSTRES, ON — 必须设置
- 遐心力、温度、内压都算作预加载
预应力模态分析的两步骤手法
含预应力的模态分析通过步骤1（静解析获得预应力状态）→步骤2（加入应力刚度的固有值分析）两步骤进行。ANSYS中通过PSTRES,ON设置自动将步骤1结果继承到步骤2。Abaqus中Initial Stress Stiffness选项具有同等功能。当应力刚度影响达10%以上时，必须考虑预应力。
预应力模态分析的实际应用
预应力模态分析的实际应用
预应力模态分析在实际中如何应用？
坎贝尔图（旋转机械）
旋转体（涡轮、压缩机、风扇）旋转速度变化时固有振动数会改变。坎贝尔图是旋转速度 vs. 固有振动数的图表，旋转速度的整数倍励振线（$f = n \times N_{rpm}/60$）与其交点是共振条件。
要多次改变旋转速度进行固有值分析吗？
对的。在10～20个旋转速度点求固有振动数，然后绘制图表。各点计算遐心力预应力→固有值分析。
电缆结构的振动
吊桥的电缆和PC梁的PC钢材由拉伸预应力决定振动数。通过测量电缆振动数来反推预应力（张力）的张力推定法存在。
测量振动数就能得到张力！很方便。
吊桥电缆张力管理已实用化。用加速度传感器测电缆固有振动数，用上式计算张力。
实务检查清单
预应力模态分析的检查清单请告诉我。
"拉伸上升、压缩下降"方向检查很重要。如果反了说明有问题。
这个方向检查是最简单最有效的验证。只要懂弦的物理就能直观判断。
PC桥梁固有振动数变化的实测管理
预应力混凝土桥的健全度管理利用固有振动数的经时变化。FEM分析确认PC钢材紧张力降低10%时，桥的1次固有振动数约下降3～5%。从定期计测数据能非破坏检测紧张力下降。国土交通省桥梁检查手册（2019年修订）在参考值中提及此手法。
预应力模态分析的软件对比
预应力模态分析的工具
预应力模态分析的求解器对比是？
Ansys有坎贝尔图专用工具。
Ansys Workbench旋转动力学解析模块包含Campbell Diagram Tool。自动参数化改变旋转速度，自动生成坎贝尔图。
选择指南
预应力模态分析是任何求解器都有的基本功能呢。
是的。坎贝尔图自动生成是Ansys方便，但用其他求解器手动改变旋转速度也能做。
NASTRAN SOL103 STATSUB功能
MSC NastranSOL 103加入STATSUB(SID)选项后，先行静解析（SID是子情形号）的应力场自动继承，执行模态分析。Boeing用此功能把翼的气动荷重场在模态分析中加入，用预应力考虑的固有值算出颤振判定用的固有值，用于787型式证明申请。
预应力模态分析的先端研究
预应力模态分析的先端研究
请介绍预应力模态分析的最前沿。
VCT（振动相关技术）
"振动数为零 = 座屈"关系被利用，从振动数下降非破坏预测座屈荷重的手法。逐步增加荷重同时测量振动数，$f^2$ vs. $P$ 曲线外推 $f = 0$ 的荷重。
不破坏结构就知道座屈荷重！在外壳座屈页也出现过。
VCT是预应力模态分析最重要的应用。用FEM模拟 $f^2$ vs. $P$ 关系，与实验比较验证。
非线性预应力模态
大变形预应力（如膨胀膜结构）的固有振动数。线性预应力模态不够，需在非线性静解析后线性化求固有值。Abaqus用NLGEOM=YES → *FREQUENCY实现。
数字孪生和VCT
用结构数字孪生实时监视振动数，推定座屈余裕（$P/P_{cr}$）。桥梁和储罐的结构健全度监测中应用。
总结
我总结预应力模态分析的先端研究。
旋转试验中的遐心预应力固有值
高速旋转翼（涡轮叶片、螺旋桨）中遐心力预应力很大，静止时和旋转时固有振动数差别很大。典型涡轮叶片在运转条件（3000～15000rpm）时1次固有振动数上升20～50%。IHI·Kawasaki Aerospace用旋转试验与联合仿真的固有振动数计测验证翼设计。
预应力模态分析的故障排除
预应力模态分析的故障
预应力模态分析常见的故障有哪些？
预应力效果没有显现
加了预加载但振动数没变。
确认项：
1. Ansys: PSTRES, ON设置了吗 — 忘了就应力刚度为零
2. 预加载静解析正常求解了吗 — 应力为零就 $[K_\sigma] = 0$
3. 荷重方向正确吗 — 打算拉伸反成压缩了吗
振动数为负（$\omega^2 < 0$）
$[K_0] + [K_\sigma]$ 的整体刚度变成零以下时 $\omega^2 < 0$。这表示座屈超越状态（结构不稳定）。预加载超过座屈荷重。
怎么办？
减少预加载或补强结构。$\omega^2 < 0$ 是FEM"这个荷重下结构撑不住"的警告。
旋转体结果与旋转速度无关
遐心力设置可能错了。在Nastran的RFORCE卡确认旋转速度和旋转轴正确设置了吗。
总结
我整理预应力模态分析故障处置。
预应力分析中固有振动数过高的情况
预应力考虑分析中固有振动数还是比实测高，原因可能是接触面滑动或隙缝影响。预应力可能低于假设（松弛、蠕变）。首先比较实测FRF和分析FRF的峰值形状，确认模式形状一致吗。一致则需重审边界条件·预应力量。不一致则需根本性重新建模。
大
更详细
错误