1. 简介
对于结构风荷载的准确评估，不仅取决于所应用的荷载模型，根本上还在于如何表示进入的大气边界层 (ABL)。在实际工程设计中，风是一个高度不稳定的、湍流的和空间相关的现象。然而，大多数基于规范的方法，包括标准欧规风模型，将这种复杂性简化为等效的静态量，以确保简单性、稳健性和保守的安全裕度。
在基于计算流体力学 (CFD) 的风模拟中，入口边界条件的定义是预测表面压力、峰值荷载和动态效应准确性的决定性因素。
本文档的方法实现了一个随机的、空间相关的湍流入口生成器，位于一个垂直的入口平面（YZ 平面），旨在用于：
- 瞬态 CFD 风模拟
- 与荷载相关的压力和力的评估
- 与类似风洞的入口条件的验证
- 超出欧规保守范围的高级应用
这种方法显著超越了主要用于静态设计荷载推导的标准欧规 (EC) 风模型，而不是用于时间分辨的流体物理学。
2. 所实现方法概述
入口生成由五个耦合的组成部分组成：
2.1 空间离散化的入口平面
- 入口定义在一个 YZ 平面上，固定 𝑥=𝑥0
- 常规网格：𝑛𝑦 × 𝑛𝑧 点（例如 20 × 100 → 2,000 探针）
- 每个网格点表示一个时间相关的速度信号
- 不进行跨流向平均（对于保留波动非常重要）
这允许垂直和横向相干效应得到物理表现。
2.2 平均风剖面（幂律，物理强制）
平均风剖面使用幂律公式定义，这是一种经验工程近似，表示在大气边界层内风速随高度的增加。该模型将速度与参考风况和地形相关的指数关联。虽然比来自相似性理论的对数边界层公式更简单，但幂律方法为工程应用提供了一种实用且广泛使用的入口定义。该公式确保了一致的垂直速度分布，参考到特定高度。
2.3 湍流模型：Kaimal 频谱
在每个入口点，使用 Kaimal 频谱生成湍流，该频谱代表大气边界层行为并定义湍流能量如何在频率间分布。这种方法捕捉到风波动的频率相关特性，而不是假设恒定行为。该模型反映了湍流强度和特征长度尺度来描述速度变化的结构。结果是，所生成的湍流包含了随频率变化的物理一致的能量内容。在每个入口点，湍流使用 Kaimal 自谱生成，与大气边界层理论一致 [1]：
这产生了频率相关的能量内容，这是欧规静态方法中所缺失的。
2.4 空间相干性和交叉频谱矩阵
应用空间相干性模型以确保入口点之间的物理现实相关性，描述湍流相似性如何随距离和频率减小。该相干函数用于定义不同位置信号之间的关系，反映出空间相关的风波动。基于此相干性，构建一个完整的交叉频谱密度矩阵，将各个频谱与其空间相关性结合。所产生的矩阵既代表了能量分布，也代表了入口处速度波动的空间耦合。为了确保入口点之间的物理现实相关性，应用距离相关的相干模型：
这产生了一个完整的交叉频谱密度矩阵：
关键特征：
- 完全空间相关的湍流
- 频率相关的相关衰减
- 捕捉入口处的实际阵风结构
2.5 稳健的 Cholesky 分解（数值稳定性）
为了在交叉频谱矩阵变得病态（尤其接近地面）时保持数值稳定性，应用了稳健的 Cholesky 分解。矩阵首先对称化以确保一致性，随后通过自适应的对角线加载来提高稳定性。如果需要，使用自动重试策略并增加正则化，确保正定性。这种方法能够稳定和可靠地产生相关的时间序列数据。因为交叉频谱矩阵可以变得病态（尤其接近地面），所以使用稳健的 Cholesky 分解：
- Hermitian 强制：
- 自适应对角线加载
- 自动重试策略以及增加正则化
这保证了正定性和稳定的时间序列生成。
2.6 时间序列构造（逆 FFT）
引入随机相位以表示信号中的湍流，同时缩放和组织频谱以确保物理现实的实值结果。然后应用逆 FFT 将频域信息转换为时间相关的速度信号。最终结果是由均值流和波动成分构成的速度场。
- 每个频率应用随机相位
- 以 sqrt(2Δ𝑓) 缩放的单边频谱
- 强制 Hermitian 对称性
- 逆 FFT → 每个入口点的时间分辨率速度信号
结果：
3. 可视化和诊断
该方法能够：
- 瞬时垂直剖面 𝑈(𝑧,𝑡o)
- 多个探头线比较（无平均工件）
- 全入口平面等高线动画 𝑈(𝑦,𝑧,𝑡)
与基于 EC 的静态剖面不同，这使得流体物理直接可观察。
4. 与标准欧规风模型的比较
表 1 总结了标准欧规 (EC) 风模型和当前随机 CFD 入口生成方法之间的基本概念和方法学差异。虽然这两种方法都旨在表示结构上的风效应，但它们基于不同的建模哲学，并服务于不同的工程目的。
表 1: 与标准欧规风模型的比较
| 方面 | 欧规 (EC) | 当前 CFD 入口方法 |
|---|---|---|
| 性质 | 静态/准静态 | 完全瞬态 |
| 湍流 | 通过因子隐式表示 | 显式时间序列 |
| 空间相关性 | 未解决 | 完全解决 |
| 频率内容 | 不存在 | Kaimal 频谱 |
| 阵风结构 | 等效静态阵风 | 物理演变 |
| 压力峰值 | 经验因子 | 流动中显现 |
| 荷载路径 | 保守包络 | 基于物理学 |
| 适用性 | 规范合规 | 高级分析和验证 |
5. 解释上的关键差异
欧规方法主要旨在通过部分安全因子嵌入湍流效应来确保安全设计荷载，但并不提供详细的荷载发展时间或机制的洞见。相反，当前方法通过解决何时、何地以及为何峰值压力发生，包括阵风渗透、分离引起的压力峰值和空间相关的荷载效应，赋予更深的物理理解。这特别有价值于轻质结构、外墙和幕墙元素、膜系统、屋顶设备以及涉及动态结构响应的研究：
📘欧规
- 设计旨在产生安全设计荷载
- 湍流效应嵌入在部分安全因子中
- 没有关于荷载何时或如何发生的信息
🌬️当前方法
- 解决何时、何地以及为何形成峰值压力
- 阵风渗透
- 分离引起的压力峰值
- 大面积的相关荷载
⭐特别重要于：
- 轻质结构
- 幕墙和外墙元素
- 膜和屋顶设备
- 动态响应研究
6. 工程影响
所提议的方法并非旨在替代欧规进行监管设计；相反，应被理解为一种补充方法，提供了超出欧规程序中通常包含的假设的额外物理洞察。它特别适合于敏感性研究、研究和验证目的、优化过程以及解释可能在欧规结果与 CFD 模拟之间产生的差异。在实际工程应用中，这种方法常常有助于澄清为何 CFD 派生的压力值与欧规预测不同，特别是关于局部峰值压力和瞬态流动效应，而不一定表示任何一种方法存在错误。
7. 总结
所提出的入口生成方法引入了物理一致的湍流建模，结合空间相干性和时间分辨的风场表示，支持稳健的数值实现。与标准欧规方法相比，该方法将分析视角从保守的荷载包络转向物理驱动的流动行为理解，从而实现更深入的洞察和更先进的风工程分析。