基于特征正交分解的桥塔风效应分析.docxVIP

基于特征正交分解的桥塔风效应分析 随着计算机技术的快速发展，基于monteithil的时间间隔抖动分析已成为计算大桥间隔振动响应的主要方法。时间间隔方法中最重要的问题是桥梁动脉风速场的随机模拟。目前此类模拟主要使用的是ARMA方法和原型谱表示法两类方法,另有小波模拟等方法。其中原型谱表示法由于具有高精度和稳定性好等特征,在工程实际中已经得到了广泛的应用。本文所描述的特征正交分解型谱表示法是近年来在对原型谱表示法进行重大改进的基础上提出的,核心是以对功率谱矩阵的特征值分解取代Cholesky分解。在随机积分的概率刻划下,它的物理机制有明确的解释。它可用来模拟一维多变量随机过程,如脉动风速场等。 当桥面和桥塔或桥墩等不处在同一高度且高差较大时,由于重力作用,风的波动从较高处向较低处的传播有一个过程,造成不同高度处的脉动风速之间存在着较大的相位不同步,对桥梁的风致振动响应存在一定的影响,称为桥塔风效应。除高差外,效应的大小还受平均风速和地貌类型的影响。当相对高差较大,一般认为超过100m时,该效应对桥梁结构风致振动的影响较为显著,而这种高差在高墩大跨桥梁飞速发展的我国并不罕见。 桥塔风效应使得两点间脉动风速交叉谱除包含常规的共谱项,还会包含余谱项,即模拟时的目标谱矩阵是一个复值矩阵,在采用原型谱表示法和特征正交分解型谱表示法时均需涉及复矩阵的分解运算,会比较耗时,而且特征正交分解本来明确的物理意义将无法表现。因此,本文基于特征正交分解型的谱表示法,推导了模拟具有桥塔风效应风场的简化公式,变复矩阵运算为实矩阵运算以提高运算效率,并以龙潭河特大桥施工阶段风场模拟为例,解释了POD振型的物理意义,并验证了该算法。 1 桥塔风效应定义 顺风向的桥梁风场可以离散化为一个一维多变量的随机向量过程V(t)=[v1(t),v2(t),…,vN(t)]T;由N个部分相干的零均值弱平稳正态随机过程vj(t)(j=1…N)组成,其功率谱密度矩阵为SVV(ω): [Svv(ω)]=[S11(ω)S12(ω)?S1Ν(ω)S21(ω)S22(ω)?S2Ν(ω)????SΝ1(ω)SΝ2(ω)?SΝΝ(ω)](1)[Svv(ω)]=??????S11(ω)S21(ω)?SN1(ω)S12(ω)S22(ω)?SN2(ω)????S1N(ω)S2N(ω)?SNN(ω)??????(1) Spq(ω)=Svpvq(ω)=√Svp(ω)?Svq(ω)?γpq(ω)(2)Spq(ω)=Svpvq(ω)=Svp(ω)?Svq(ω)????????????√?γpq(ω)(2) Svp(ω)为某点处的单点功率谱密度函数,γpq(ω)为vp和vq间的相干函数。 根据我国的桥梁抗风设计规范,单点功率谱密度函数Svv(ω)可取为Kaimal谱: Svv(ω)=122002πu2*zU(z)1[1+50ωz2πU(z)]5/3(3)Svv(ω)=122002πu2?zU(z)1[1+50ωz2πU(z)]5/3(3) 其中u*为气流的剪切摩阻速度,z为离地面高度坐标,U(z)为z高度处平均风速,均可按照规范取值。在计入桥塔风效应时,相干函数γpq(ω)可取为: γpq(ω)=exp{-ω2π[C2x(xp-xq)2+C2z(zp-zq)2]1212[Uz(zp)+Uz(zq)]}?exp[-iθpq(ω)]γpq(ω)=exp{?ω2π[C2x(xp?xq)2+C2z(zp?zq)2]1212[Uz(zp)+Uz(zq)]}?exp[?iθpq(ω)] (4) 其中xp、xq为p、q两点的侧向坐标(可参见图1),两个常数Cz、Cx分别取10和16。式(4)由两部分组成,前一项为常规的实值相干函数项,后一项为是复值项,代表桥塔风效应,即两点之间脉动风速的相位差: θpq(ω)=ω(zp-zq)vp,qapp(5) vp,qapp定义为在p、q两点间风的视在波动传播速度: vp,qapp=1Cθ[Uz(zp)+Uz(zq)]2(6) Cθ为和传播介质有关的波动速度系数,在空气中可取为Cθ=5.5。从式(5)和式(6)可以看出,桥塔风效应实际上就是在高度较低的点的脉动风速和较高的点处的风速之间存在的时间滞后,它和传播速度以及两点之间的高差有关,在同一高度之间则不存在此效应。它非常类似于地震动中的行波效应。另一种对桥塔风效应的定义方式可参见文献,其物理意义不如式(5)明确。式(5)可改写为直接由滞后时间表述: θpq(ω)=ω(Tq-Tp) (7) Τp=p-1∑i=1ω(zi-zi+1)vi,i+1app(8) Tp为p点和I点(高度基准点)之间的累计传播滞后时间。 桥塔风效应不仅会使交叉谱Spq(ω)是一个复值,它还会使p、q两点间风速的互相关函数峰值不再位于零滞后时间