风车桥理论简介.docVIP

第二节 考虑脉动风作用的车桥系统振动响应分析方法 列车过桥时列车一桥梁系统是一个时变系统，而风对结构的作用以及风荷载的计算又是非常复杂的空气动力学问题，尤其对于大跨度桥梁许多国家的规范多主张对其应另行处理，许多问题都需要进一步总结研究[37]，所以对脉动风荷载作用下大跨度铁路桥梁与列业系统振动响应分析难度是相当大的。学生认为采用时域分析法用计算机模拟求解是个比较好的方法，首先根据风速谱函数，针对大跨度桥梁的特点，通过数字模拟的方法获得考虑了空间相关特性的风速时程曲线，以些再根据结构的空气力系数模拟产生脉动风荷载，在前述列车一桥梁时变系统分析模型的基础上，将它们加到车桥系统动力学方程中求解。采用这种计算机仿真模拟方法，基本能反映车桥系统振动的主要特征。 （一）大跨度桥梁随机风场的模拟 大跨度桥主梁上的风速不仅包含了在x, y, z三个方向上的脉动，而且还沿着主梁长度方向变化。本文仅考虑水平方向上的脉动分量，将风场近似看作沿梁长度方向上若干点处的随机风皮的总和，则该随机风场可以看作是一个多维多变量的平稳高斯随机过程。 模拟随机过程的一类方法是基于三角级数叠加的谐波合成法。考虑一个一维，n变量，零均值的高斯随机过程。它包含个变量。 其互谱密度矩阵为：（3.1） 根据Shinozuka的理论，随机过程的样本可以由下式来模拟[31]： （3.2） 式中各变量意义如下： N为一充分大的正整数， ，为频率增量， 为载止圆频率，即当时， 为均匀分布于（0，2）区间的随机相位， 是矩阵是的Cholesky分解，即： （3.3） （3.4） 由于通常情况下为一复数矩阵，且不一定正定，因此通常也得矩阵，其对角元素为实数，非对角元素为复数。是其共轭转置矩阵。为的复角 根据Shinozuka等人的研究，为增大模拟样本的周期，可按如下取值[32]： （3.5） 可以证明，当时，式（3.2）模拟的随机过程满足式（3.1）的目标谱。 避免式（3.2）的模拟结果失真，的增量须满足以下条件： （3.6） 式（3.2）模拟的随机过程的周期为： （3.7） 由以上可知，只要已知了，恰当地选择了N，和，就可以获得好的随机过程的样本。但要对这样的复数矩阵进行Cholesky分解是比较繁琐的，通常要借助于递推公式来计算。由于是的函数，由式（3.2）可知对每个频率进行一次余弦项叠加时，都要进行一次的Cholesky分离，其计算量是可观的。但对于本文研究的大跨度桥梁，该计算可以大大简化[37]。 在工程实验中，空间不同点之间的互谱密度通常用实数函数来表达[36]，因此式（3.1）可以简化为实数矩阵。假设我们要模拟桥主梁上若干点的水平风速，则可以将它们看作是一个一维多变量的高斯随机过程。一般情况下（即地形起伏不大，离在一定高度），我们可以近似认为水平风谱沿水平方向是不变的[36]。因此主梁上各点的水平风谱是相同的，即有下式成立： （3.8） （3.9） 式中为j点和m点之间的水平距离，为相干函数。 根据Davenport的研究可得[36]： （3.10） 式中为桥面处的平均风速，C上的函数。 （3.11） 将式（3.8）~（3.11）代入式(3.1)中可得： （3.12） 式（3.12）的Cholesky分解可以显式地写为。 式中： （3.14） 写成解析式则为： （3.15） 由于，故为实矩阵，因此有 （3.16） 本文学生运用以上方法，针对我国京沪高速铁路南京长江大桥三塔斜拉桥方案（L = 84+160+2*488+160+84=1464米），模拟了桥面上沿跨度方分布，间距为50m的30个点的水平脉动风速。模拟计算的主要参数和结果如下：长度：L = 1464m；主梁离地高度：z = 41.0m；地面精糙度：z0 = 0.01m；主梁处平均风速：U(z)=25.0m/s；模拟点数：n = 30；模拟点间距：；截止频率：；频率等分数：N = 256；模拟采样时距：；目标谱：Kaimal谱，Kaimal水平脉动风速谱的计算公式为[36]： （3.17） 其中 图3-1模拟得到的第1、2、16点的风速时程的片段。由图3-1可以看出，点1、2之