斜拉索上瞬时风速与风荷载的数值模拟.pdf

斜拉索上瞬时风速与风荷载的数值模拟 方治华 ，李 晨 (内蒙古科技大学 建筑与土木工程学院，包头 014010) 摘 要：基于线性自回归模型法，确定了随时间和高度变化的脉动风模拟方法，编制了瞬时风速和风荷载模拟程 序。将此方法应用到斜拉索结构，对苏通大桥Z34号拉索表面的风速、风荷载进行模拟。对所得脉动风时程进行谱分析 和数理统计分析。结果表明，模拟点处模拟谱与 目标谱吻合较好 ，证明了模拟结果的准确性。 关键词：斜拉索；风速时程；风压时程；功率谱密度 中图分类号：TU31 文献标识码 ：A 随着对空间结构跨越能力的要求不断提高，结构 AR法是线性滤波法的一种，它用 自回归模型来表示。 趋于柔性化，风荷载成为结构振动控制的关键因素。 本文用 自回归模型模拟多维脉动风速过程 ，设 m个点 目前，结构的抗风研究主要采用风洞试验和数值模拟 空间相关脉动风速时程 ( ，l，，Z，t)列向量的P阶模 方法 J。风洞试验需要特殊的试验设备和较高的技 型可表示为： 术，而且需要投人较多的资金，因而数值模拟方法具有 ( ，y，Z，t)= 相当的优势，得到广泛的应用。数值模拟方法主要有 P ∑ v(x，Y，Z，t—kAt)+N(t) (1) 谐波合成法和线性滤波法。线性滤波法中的AR法应 用线性 自回归模型，可根据场地、风速谱特性、建筑的 上式中：X=[ 一， ]，Y=[Y 一，Y]，Z：[ ．．， 特点等条件，使风速时程的模拟尽量接近实际值，因此 ]T，(，Y，)为空间第i点坐标(i=1，…．．，m)；p为 得到广泛的应用。例如，李元齐 模拟了大跨度空间 模型的阶数；At是模拟时间步长， 为m×m阶的自回 结构，李春祥 模拟 了超高层建筑的脉动风速时程 ， 归系数矩阵；Ⅳ(t)是方差为 1的独立随机过程向量。 YongleLil4J模拟了斜拉桥上的随机风速。本文将此 应用(1)式模拟脉动风速时程时，首先必须确定 自 ， 方法推广到模拟斜拉索表面的瞬时风速与风荷载，并 回归系数矩阵 和独立随机向量 N(t)。随机风过程 对风速功率谱、自相关函数的模拟值与相应 目标值的 的协方差矩阵R与回归系数矩阵 之间的关系可写成 吻合程度进行检验 ，验证 了用线性 自回归模型模拟索 矩阵形式 ： 类结构的准确性。 R× =[R ，0]T =[，， ，…， r P 1 自回归模型的脉动风速模拟 R =R。+∑ R(kat) (2) 风速时程模拟最常用 的方法是线性滤波法。它是 式中：J为m阶单位矩阵； 为 (p+1)m ×m阶矩阵， 将随机过程抽象为满足一定条件的白噪声，然后经过某 0。为pm×m阶零矩阵； 为(p+1)m×(P+1)m阶 自 一 假定系统进行适当变换，拟合出该过程的时域过程。 相关矩阵，其形式如下 ： R11(O) R12(At) R13(2At) R1(p+1)[pAt] R2(At) R22(0) R23(At)