[发明专利]一种平行双幅桥静风稳定分析方法

权利要求书

1.一种平行双幅桥静风稳定分析方法，其特征在于，包括如下步骤：

1)、构建数值风洞或物理风洞模型，然后基于该模型获得上游主梁和下游主梁的三分力系数与四个参数的对应关系；所述三分力系数包括阻力系数、升力系数和升力矩系数，所述四个参数包括上游主梁的有效风攻角α、下游主梁的有效风攻角α'、两个主梁的扭转中心间距d和上下高差h；所述获得上游主梁和下游主梁的三分力系数与其四个参数的对应关系的具体步骤如下：将上游主梁断面和下游主梁断面同时纳入数值风洞或物理风洞模型，并变换上述四个参数，测得不同组合下的三分力系数；将四个参数不同的组合及其对应的上、下游主梁断面三分力系数以可调用数组方式进行存储并获得上游主梁和下游主梁的三分力系数与其四个参数的对应关系；

2)、建立空间有限元模型，同时包含上游桥梁和下游桥梁的信息，上游桥梁和下游桥梁的相对位置关系与实桥一致，然后进行自重作用下几何非线性求解；

3)、设定初始风速V₀和风速步长ΔV，当前风速V_i＝V₀，并设定迭代次数上限N_max；

4)、从空间有限元模型中提取上游主梁各单元的扭转角θ_i、下游主梁各单元的扭转角θ_i'、上游主梁各单元与下游主梁对应单元的中心间距d_i和上下高差h_i；

5)、根据步骤4)中获得的上游主梁各单元的扭转角θ_i、下游主梁各单元的扭转角θ_i'以及初始的风攻角α₀，计算上游主梁各单元的有效风攻角α_i(＝α₀+θ_i)和下游主梁各单元的有效风攻角α_i'(＝α₀+θ_i')；然后根据四参数组合(α_i,α_i',d_i,h_i)利用五维空间内插法分别计算上游主梁和下游主梁各单元的三分力系数；

6)、在当前风速V_i下，计算作用在上游主梁各单元的横向风荷载P_Hi、竖向风荷载P_Vi和扭转力矩P_Mi，及下游主梁各单元的横向风荷载P′_Hi、竖向风荷载P′_Vi和扭转力矩P′_Mi；

7)、在上游主梁、下游主梁的各单元上分别施加横向风荷载P_Hi和P′_Hi、竖向风荷载P_Vi和P′_Vi、扭转力矩P_Mi和P′_Mi，进行桥梁结构几何非线性求解，获得上游主梁各单元的扭转角θ_i、下游主梁各单元的扭转角θ_i'、上游主梁各单元与下游主梁对应单元的扭转中心间距d_i和上下高差h_i，根据下式判断这四个参数的欧几里得范数是否分别小于等于允许值ε₁、ε₂、ε₃和ε₄：

式中，N为主梁的单元总数；k为当前荷载步编号；i为梁单元序号；

8)、如果上述四式中的任一式不成立，则重复步骤5)-7)；如果迭代次数达到迭代次数上限N_max，则当前风速难以收敛，此时令当前风速V_i+1＝V_i-ΔV，然后缩短风速步长，返回步骤5)，重复步骤5)-7)；如果风速步长小于预定值，则计算结束；如果上述四式全部成立，则当前风速计算结果收敛，输出计算结果，其中，计算结果中包括桥梁结构变形参数，此时令当前风速V_i+1＝V_i+ΔV，重复步骤5)-7)；

9)、根据步骤8)得到的计算结果，得到上游桥梁和下游桥梁的静风失稳临界风速和失稳形态。