塔架涡激振动的计算方法——基于EuroCode和GB

1. 涡激振动简介

涡激振动是结构在一定条件下，由于旋涡周期性脱落带来的横风向周期性激励与结构自振频率一致导致的共振现象。

发生涡激振动需要满足两个条件：
1、一个是在结构高度范围内存在共振区，即结构顶点风速大于共振风速（也称为临界风速）；
2、发生强风共振，即结构处于跨临界区（临界区的定义可以参考任意一本风工程教材）。

本文主要针对混凝土塔架在混凝土段吊装完毕以后，未张拉预应力钢绞线前可能会发生的涡激振动进行计算。

2. EuroCode（欧标）中涡激振动的计算流程

EuroCode中计算涡激振动方法可以参考EuroCode1-4 附录E.1。

2.1 发生涡激共振的判定（检验是否需要进一步计算）

按照EC1-4，满足下式可以不再考虑涡激振动：
$$v_{crit,i}>1.25\cdot v_m\text{\hspace{1em}(1)}$$

其中，$v_m$为验算高度处（一般取结构顶部）10min平均风速，风电行业的取值一般由风资源直接提供，不考虑其他修正；$v_{crit,j}$为临界风速，它是导致漩涡脱落频率与结构自振频率相等的风速，表达式如下：
$$v_{crit,i}=\frac{d\cdot n_{i,y}}{S_t}\text{\hspace{1em}(2)}$$
$S_t$——斯托拉哈数，取0.18；
$d$——取3/4高度处塔筒外直径；
$n_{i,y}$——结构$i$阶固有频率。

2.2 惯性力的计算

本节先看两个无量纲常数——斯科拉顿数和雷诺数，然后给出惯性力及其依赖参数的计算方法。

2.2.1 斯科拉顿数

$$S_c=\frac{2\cdot {\delta }_s\cdot m_{i,e}}{\rho \cdot d^2}\text{\hspace{1em}(3)}$$
${\delta }_s$——结构阻尼，取值0.03（按照EuroCode1-4 Table F.2取值）；
$m_{i,e}$——模态$i$的单位长度等效质量；
$\rho$——空气密度；

2.2.2 雷诺数

$$Re(v_{crit,i})=\frac{d\cdot v_{crit,i}}{\nu }\text{\hspace{1em}(4)}$$
$\nu$——空气运动粘滞系数($\nu \approx 1.5\times 10^{-6}{m}^2/s$)

2.2.3 单位长度上的惯性力

$$F_w(s)=m(s)\cdot (2\cdot \pi \cdot n_{i,y}{)}^2\cdot {\varphi }_{i,y}(s)\cdot y_{F,max}\text{\hspace{1em}(5)}$$

2.2.4 最大位移的计算

2.2.4.1 方法1

$$\frac{y_{F,max}}{d}=\frac{1}{S_t^2}\cdot \frac{1}{S_c}\cdot K\cdot K_w\cdot c_{lat}\text{\hspace{1em}(6)}$$

EuroCode提供了两种计算最大位移的方式，这里给出的是普遍性的计算方法，另一方法在本节末尾2.2.5小节中给出。

其中，
$$c_{lat}=\{\begin{array}{c}{c}_{lat,0},v_{crit,i}/v_{m,Lj}⩽0.83\\ (3-2.4\cdot \frac{v_{crit,i}}{v_{m,Lj}})\cdot c_{lat,0},0.83<v_{crit,i}/v_{m,Lj}<1.25\\ 0,\end{array}$$