脉动风速谱matlab程序,1.脉动风速谱.PDF

1. 脉动风速谱

根据统计分析可知，任意一点的风速 V ( z, t) 为平稳Gauss随机过程，可以分成平均风速

[5]

v (z) 与零均值脉动风速v(z, t) 。平均风速的高度变化律可以采用对数或指数形式 。

V ( z, t) v (z ) +v(z , t) (1)a

ln(z / z )

v (z) v 0 (1)b

s ln(z / z )

s 0

式中：v 为标准高度处的平均风速，z 为标准高度 (我国取10m)，z 为风速为 0 的高

s s 0

度(与地貌有关)。

脉动风的风速谱密度模型主要有 Davenport 谱、Kaimal 谱、Vickery 谱等，本文选择

Davenport 谱，其表达形式为：

2

2 r ωL

Sv (ω) 4Kv10 2 4 / 3 ，r (2)

r

ω(1+ ) v10

式中：K 为地面粗糙度系数，L 为湍流积分尺度，此处取为常数 1200m，v10 为高度 10m

处的平均风速。

由于 Davenport 谱是与高度无关的，所以引入反映与高度相关的随机函数，其均值为 0，

某二点的互相关函数取为：

⎧⎪ (z2 −z1 )2 ⎫⎪

Rz (z1, z 2 ) exp ⎨=− ⎬ (3)

⎩⎪ 3600 ⎭⎪

其中z1 、z2 为两点的高度。

2. 随机振动求解

[6]

由连续体系微元的动力平衡关系，可以推导出连续体系动力学方程 ：

2 2 ⎡ 2 ⎤

∂ d (z, t) ∂d (z ,t) ∂ ∂ d (z, t)

m(z) 2 +c(z ) + 2 ⎢EI (z )