流体压强与浮力 — 帕斯卡原理・阿基米德定律与CFD流体分析
目录
1. 流体静力学的工程挑战
潜水艇壳体设计最难的是什么？海水压强那么大，怎么保证壳体不被压垮？
最难的是抗外压屈曲。内部承受拉应力的容器，按照屈服强度设计就行；但潜水艇壳体外部受压，薄壳结构很容易在达到材料屈服强度之前就发生弹性屈曲——就像用手捏饮料瓶，在铝壳变形之前先出现凹陷。水深300米对应约3MPa外压，屈曲安全系数的FEM分析要求临界屈曲压力≥1.5倍工作压力，这对壳体厚度和肋骨间距的优化提出了极高要求。
流体压强在工程中无处不在：管道内压、水坝水压、地下工程土压、液压系统工作压力、深海装备外压……每一种都需要精确的CAE分析来保证结构安全。
2. 压强：力的面密度
压强（Pressure）定义为单位面积上的法向力：
压强单位换算：
| 单位 | 换算关系 | 工程应用场景 |
|---|---|---|
| Pa (Pascal) | 1 Pa = 1 N/m² | SI基本单位 |
| kPa | 1 kPa = 1000 Pa | 气象、低压气流 |
| MPa | 1 MPa = 10⁶ Pa = 1 N/mm² | 结构应力、液压 |
| bar | 1 bar = 10⁵ Pa ≈ 1 atm | 压力容器、管道 |
| atm | 1 atm = 101,325 Pa | 大气压基准 |
| psi | 1 psi ≈ 6895 Pa | 英制，航空液压常用 |
3. 静水压力分布
液体中深度 $h$ 处的静水压强：
其中 $p_0$ 是液面处压强（通常为大气压 101.325 kPa），$\rho$ 是液体密度，$g = 9.81$ m/s²。
海水（$\rho = 1025$ kg/m³）中压强随深度的变化：
| 深度 | 压强（表压） | 工程意义 |
|---|---|---|
| 10 m | ≈ 0.1 MPa（1 atm） | 浅水游泳/潜水 |
| 100 m | ≈ 1 MPa（10 atm） | 轻型潜水设备极限 |
| 300 m | ≈ 3 MPa（30 atm） | 常规潜水艇作战深度 |
| 1000 m | ≈ 10 MPa（100 atm） | ROV深海机器人 |
| 11,000 m（马里亚纳） | ≈ 110 MPa（1100 atm） | 全海深载人潜水器 |
4. 帕斯卡原理与液压系统
帕斯卡原理：封闭流体中某一点压强的变化，会无衰减地传递到流体的每一点。
液压缸的力放大：
汽车液压千斤顶：小活塞面积 $A_1 = 5$ cm²，大活塞面积 $A_2 = 500$ cm²，力放大100倍。以50 N的踩踏力可以举起5000 N（500 kg）的重物。
工程液压系统典型参数：
| 系统类型 | 工作压力 | 典型应用 |
|---|---|---|
| 低压液压 | <7 MPa | 翻斗车、农业机械 |
| 中压液压 | 7-21 MPa | 工程机械（挖掘机） |
| 高压液压 | 21-35 MPa | 航空飞行控制系统 |
| 超高压液压 | >35 MPa | 锻造设备、水射流切割 |
5. 阿基米德浮力定律
浸入流体中的物体受到向上的浮力，大小等于被排开流体的重力：
浮沉条件：
- $\rho_{\text{object}} > \rho_f$：沉底
- $\rho_{\text{object}} = \rho_f$：中性浮力（悬浮）
- $\rho_{\text{object}} < \rho_f$：漂浮
潜水艇通过调节压载水舱（Ballast Tank）的注排水量改变总密度，实现上浮下潜：
- 下潜：向压载水舱注入海水 → 总密度增大 → 重力 > 浮力 → 下沉
- 上浮：用高压空气将压载水舱中海水排出 → 总密度减小 → 浮力 > 重力 → 上浮
海上浮吊安装海上风机时，浮吊自身的稳定性是怎么保证的？浮力够用就不会翻？
光有浮力不够，还要看稳心高（GM）——浮心B（浮力作用点）和重心G的相对位置。G在B上方时如果倾斜，浮力会产生一个"扶正力矩"把船扶回来，这时是稳定的（GM>0）。GM越大越稳，但太大了就会摇摆剧烈，不利于精密吊装。海上风机安装要在3~4级海况（波高2m以下）才能进行，船舶的动态稳定性要用船体FEM+水动力耦合软件（如ANSYS AQWA）专门评估。
6. 伯努利方程：流体动力学入门
沿流线的伯努利方程（不可压缩、无黏、稳态流）：
三项分别代表：
- $p$：静压（压力势能）
- $\frac{1}{2}\rho v^2$：动压（动能）
- $\rho g z$：重力势能
伯努利方程在工程中的应用：
| 应用 | 原理 | 计算公式 |
|---|---|---|
| 皮托管（测速） | $p_{\text{total}} - p_{\text{static}} = \frac{1}{2}\rho v^2$ | $v = \sqrt{2\Delta p/\rho}$ |
| 文丘里管（流量计） | 截面积变化引起压差 | $Q = A_2 \sqrt{2\Delta p/(\rho(1-(A_2/A_1)^2))}$ |
| 飞机升力（简化） | 上翼面快→低压，下翼面慢→高压 | $\Delta p = \frac{1}{2}\rho(v_{\text{top}}^2 - v_{\text{bottom}}^2)$ |
注意：伯努利方程的局限性——不适用于有黏性损失、旋转流（如泵内部）、湍流或可压缩（马赫数 > 0.3）情况，这些需要CFD仿真。
7. 工程实例：潜水艇耐压壳设计
以蛟龙号载人深潜器（下潜深度7000 m）为例：
设计外压（7000 m深）：
钛合金球形耐压壳（内径2.1 m，材料Ti-6Al-4V，$\sigma_y = 900$ MPa）的壁厚估算：
薄壁球壳的环向应力（Laplace公式）：
若要求安全系数 $S = 3$：
实际还需考虑外压屈曲（弹性失稳），这才是决定壁厚的关键设计工况。使用Abaqus进行线性屈曲（特征值）分析，然后再做包含几何缺陷的非线性屈曲分析，最终确定安全壁厚。
8. 水坝静水压力结构分析
重力坝上的静水压力是典型的三角形分布载荷：
其中 $H$ 是水位高度，$z$ 是计算点高度。作用在单位宽度坝面上的合力：
合力作用点在距坝底 $H/3$ 处。
在FEM中，静水压力作为面载荷施加在上游坝面，载荷大小按深度线性变化。此外还需考虑：
- 扬压力（坝基渗流产生的向上水压力）
- 泥沙淤积压力
- 地震惯性力（Westergaard附加质量法）
- 温度梯度引起的热应力
9. CFD基础：N-S方程与压力求解
N-S（Navier-Stokes）方程是描述流体运动的控制方程：
连续性方程（质量守恒）：
动量方程（动量守恒，牛顿第二定律的流体版）：
左边是惯性项（加速度），右边依次是压力梯度、黏性力、体积力（重力）。
不可压缩流的压力求解是CFD的核心挑战之一——压力没有自己的输运方程，需要通过SIMPLE、PISO等压力-速度耦合算法迭代求解。
总结
流体压强和浮力是工程中最常遇到的外载荷形式之一。从帕斯卡原理到伯努利方程，从潜水艇外压屈曲到水坝结构分析，流体静力学和动力学的基础知识直接决定了载荷边界条件的正确设置。CFD作为流体仿真的核心工具，建立在N-S方程的数值求解之上，而N-S方程本质上是牛顿第二定律在流体中的表达。