流体压力与浮力 — 帕斯卡原理·伯努利到CFD分析
1. 流体压力与CAE — 工程师为何关注压力
水、空气、油等流体涉及所有工程系统。压力是流体向构造物施加力的根源，决定了设计的可行性。管道破裂、大坝溃堤、飞机增压舱破坏、深海探测器座屈——这些都是压力估算错误或意外压力载荷引起的事故。
对于CAE工程师来说，"压力"不仅仅是一个数值，而是直接输入构造物边界条件的载荷。在CFD（计算流体力学）中，求解压力场本身是分析的目的，在流体-结构耦合分析（FSI）中，将流体侧求得的压力转移到结构侧进行FEM分析。深刻理解压力物理是精确设置分析模型边界条件的基础能力。
2. 压力的定义和单位
压力定义为单位面积上作用的力。
$$p = \frac{F}{A} \quad \text{[Pa = N/m}^2\text{]}$$压力是标量量，无方向。但压力作用在面上时，力沿着面的法线方向（向外）作用。
绝对压力·表压·真空压
压力有3种表示方法，在CAE边界条件设置中容易混淆。
- 绝对压力（Absolute Pressure）：以真空（$p = 0$）为基准的压力。热力学计算中使用。
- 表压（Gauge Pressure）：以大气压为基准（$p_{gauge} = 0$）的压力。压力表显示的值。
$p_{abs} = p_{gauge} + p_{atm}$ - 真空压（Vacuum Pressure）：低于大气压时的差分。$p_{vacuum} = p_{atm} - p_{abs}$
标准大气压为 $p_{atm} = 101.325\,\text{kPa} = 1\,\text{atm} = 1013.25\,\text{hPa}$
常见压力量级
| 对象 | 压力（绝对压） | 说明 |
|---|---|---|
| 标准大气压 | 101.3 kPa | 海平面大气的重量 |
| 汽车轮胎（表压） | ~230 kPa | 约2.3个大气压（表压） |
| 血压（收缩期） | ~16 kPa above atm | 120 mmHg |
| 水深100 m | ~1.1 MPa | 大气压 + $\rho g h$ |
| 深海10,000 m（马里亚纳海沟） | ~100 MPa | 约1000个大气压 |
| 核反应堆冷却水 | ~15 MPa | 300℃下水仍不沸腾 |
| 燃气轮机燃烧室 | 1〜4 MPa | 高压下燃烧效率提升 |
3. 流体静力学（静水压）
静止流体中的压力分布用流体静力学基本式描述。
$$p = p_0 + \rho g h$$其中 $p_0$ 是液面压力（通常为大气压），$\rho$ 是流体密度 [kg/m³]，$g$ 是重力加速度 [m/s²]，$h$ 是液面以下的深度 [m]。
重要性质：压力在水平面内均匀分布，仅随深度变化。容器形状（宽或窄）无关。这有时称为"帕斯卡悖论"的基础。
大坝的水压分布和合力
高度为 $H$ 的竖直面（大坝墙面等）上的静水压合力 $F$ 通过对压力分布积分求得。
$$F = \int_0^H \rho g h \cdot w \, dh = \frac{1}{2}\rho g H^2 w$$其中 $w$ 是大坝宽度。合力的作用点（压力中心）位于底面以上 $H/3$ 处（三角形分布的重心）。
*DLOAD, TYPE=HYDRO 指定水面位置和流体密度，它会自动计算随深度变化的压力，分配到各个面单元上。
压力水头·速度水头·位置水头
静水压式也可表示为压力水头（m）的形式。
$$\frac{p}{\rho g} + z = \text{const（静止流体）}$$这是伯努利定理（流速为零的情况）的特殊形式。
4. 帕斯卡原理和液压系统
帕斯卡原理：密闭的静止流体中某一点的压力变化会等量地传递到流体各部分。
这是液压系统的基础。通过"液压增力"能用小力产生大力，原理如下。
$$\frac{F_2}{F_1} = \frac{A_2}{A_1}$$面积 $A_1 = 1\,\text{cm}^2$ 的液压缸施加 $F_1 = 100\,\text{N}$ 的力，面积 $A_2 = 100\,\text{cm}^2$ 的液压缸可产生 $F_2 = 10{,}000\,\text{N} = 10\,\text{kN}$ 的力。
飞机的液压执行器
民用飞机的飞行控制（副翼、升降舵、方向舵）通常采用工作压力为 20～34 MPa 的液压系统。飞行员轻轻推动操纵杆就能产生数kN的舵面载荷。现代飞机将该液压系统冗余设计为3套以上，保证单一故障时安全。
与CAE流体-结构耦合分析的连接
液压系统分析需要用CFD计算配管内的非定常压力波动（水击现象），再将该压力荷载输入配管FEM模型进行流体-结构耦合（FSI）分析。快速关闭阀门会产生 $\Delta p = \rho c \Delta v$（$c$：音速）的压力尖峰，可能损坏配管。
5. 阿基米德原理和浮力
浸入流体中的物体受到等于该物体排除流体重量的向上浮力。
$$F_B = \rho_f g V_{sub}$$其中 $\rho_f$ 是流体密度，$V_{sub}$ 是物体的水没体积。
浮起条件： $\rho_{obj} < \rho_f$（平均密度小于流体）
船舶的吃水线和排水量
船舶沉到自重等于浮力的深度（吃水）处。
$$W_{ship} = \rho_{seawater} g V_{displaced}$$大型油轮（排水量30万吨）的排除水量为 $V_{displaced} = 300{,}000 / 1.025 \approx 293{,}000\,\text{m}^3$。装货时吃水增加——这个简单关系是船舶工程的出发点。
海洋平台的FEM浮力荷载设置
海上风力发电浮体式平台或石油钻探平台的浮力直接影响结构。
- 水没部材受静水压作用（压缩载荷）。
- 向上浮力作为等效竖向压力荷载输入FEM。
- 波浪导致的动态浮力变化通过考虑波高、波周期的非定常CFD分析或莫里森方程（后述）评估。
6. 伯努利定理
对非粘性·不可压缩·稳定流沿流线应用能量守恒，得伯努利定理。
$$p + \frac{1}{2}\rho v^2 + \rho g z = \text{const}$$三项分别对应压力能·动能·位能（单位体积）。流速 $v$ 大的地方压力 $p$ 小——虽然反直觉，但可理解为能量的相互转换。
forceCoeffs 函数对象自动计算升力系数 $C_L$ 和阻力系数 $C_D$。湍流模型（一般用 $k$-$\omega$ SST）的选择对精度影响很大，需注意。
伯努利定理的应用例
| 应用 | 原理 | 工程意义 |
|---|---|---|
| 文丘里管（流量计） | 截面缩小→流速增→压力降 | 由 $\Delta p$ 计算流量 |
| 翼的升力（粗略理解） | 上面高流速→低压 | 飞机·风力发电叶片设计 |
| 卡门涡（流速测量） | 后流压力变动频率 | 超声波流量计的一种原理 |
| 皮托管（飞行速度计） | 滞止点压 = $p_0 + \frac{1}{2}\rho v^2$ | $v = \sqrt{2(p_0-p)/\rho}$ |
7. 表面张力·毛细现象
液体表面分子间力不对称，产生表面张力 $\gamma$ [N/m]。20℃的水表面张力约为 0.072 N/m。
$$\Delta p = \frac{2\gamma}{r} \quad \text{（球形液滴的拉普拉斯压）}$$半径 $r$ 越小，气泡或液滴内外的压力差（拉普拉斯压）越大。半径1 μm的气泡，$\Delta p = 2 \times 0.072 / 10^{-6} = 144\,\text{kPa}$——远超大气压。
毛细现象和接触角
细管内液面上升高度为:
$$h = \frac{2\gamma \cos\theta}{\rho g r}$$$\theta$ 是接触角（表示润湿性）。$\theta < 90°$（亲水面）时液面上升，$\theta > 90°$（疏水面）时下降。
微流体芯片的CAE分析
医疗诊断芯片（LOC：Lab-on-a-Chip）流道宽度为10～100 μm。这种微观尺度上表面张力压倒重力和粘性力，控制流动行为。用COMSOL Multiphysics或OpenFOAM的VOF（体积分数法）追踪气液界面，模拟芯片内液体输运。
8. 粘性流体与流动CAE模型
真实流体具有粘性。用动力粘性系数 $\mu$ [Pa·s]和运动粘性系数 $\nu = \mu/\rho$ [m²/s]表示流体的"粘性"程度。
Navier-Stokes方程的压力项
不可压缩粘性流体的连续方程和运动方程（Navier-Stokes方程）：
$$\nabla \cdot \vec{v} = 0 \quad \text{（连续方程，不可压缩）}$$ $$\rho \left(\frac{\partial \vec{v}}{\partial t} + (\vec{v}\cdot\nabla)\vec{v}\right) = -\nabla p + \mu \nabla^2 \vec{v} + \rho \vec{g}$$右侧第一项 $-\nabla p$ 是压力梯度力。流体沿压力从高向低的方向加速。CFD将这个压力场和速度场联立求解——OpenFOAM的SIMPLE算法是典型手法。
OpenFOAM中的压力边界条件
| 边界类型 | 压力边界条件 | 速度边界条件 |
|---|---|---|
| 流入口（inlet） | zeroGradient | 固定值（$v = v_{in}$） |
| 流出口（outlet） | 固定值（$p = 0$） | zeroGradient |
| 壁面（wall） | zeroGradient | 无滑移（$v = 0$） |
| 对称面（symmetry） | symmetryPlane | symmetryPlane |
压力-速度耦合的注意事项
在不可压缩流中，压力方程（泊松方程）由速度发散为零的条件决定压力。OpenFOAM的 p_rgh（表压减去静水头的修正压力）可简化有重力影响的流动（自然对流等）的设置。
9. 实践：海洋风力发电单桩基础的波浪荷载
海上风力发电最常见的基础是单桩（Monopile）方式。直径6～10 m、长40～80 m的钢管打入海底，支撑风机塔架。波浪荷载计算广泛采用莫里森方程。
莫里森方程
$$F = C_M \rho \frac{\pi D^2}{4}\dot{u} + C_D \frac{1}{2}\rho D |u|u$$第一项是惯性力（$C_M$：惯性系数，$\approx 2.0$），第二项是阻力（$C_D$：阻力系数，$\approx 1.0$）。$u$ 是水质点水平速度，$\dot{u}$ 是其时间导数（加速度），$D$ 是构件直径。
用波浪理论（线性波理论）计算水质点速度，代入莫里森式得沿桩长的分布荷载，输入FEM梁模型计算根部弯矩和疲劳荷载。
避免共振的设计
单桩的固有振动频率应远离波的主周期（5～15秒）和叶轮转速频率（1P：~0.2 Hz）。若固有频率与波或旋转励振重合，会产生大的动态放大，疲劳寿命急剧下降。
- 软-刚设计：1P与3P（三叶片）之间配置固有频率（约0.25～0.35 Hz）
- 用FEM模态分析确认，通过优化桩径、壁厚、埋深调整