《武汉工程大学学报》 2017年05期
477-481
出版日期：2017-12-19
ISSN:1674-2869
CN:42-1779/TQ
大跨径钢箱梁斜拉桥合理成桥索力的优化
斜拉桥是一种由桥塔、加劲梁、斜拉索三种基本构件组成的组合桥梁结构体系，其工作原理是利用桥塔引出的斜拉索作为主梁的弹性中间支承，借以降低主梁的截面弯矩，减轻梁的自重，提高跨越能力［1-2］. 斜拉桥属于高次超静定结构，其重要特征是索力可以调整，而索力的大小与分布决定了斜拉桥的内力和线形，因此如何通过索力调整使结构的受力处于一个相对较为合理的状态是斜拉桥设计时一个非常重要的问题. 成桥状态的恒载索力也往往成为评价斜拉桥结构优劣的重要指标，因此斜拉桥索力优化方法的研究具有重要的工程意义. 确定斜拉桥合理成桥索力的理论计算方法有很多，主要有指定状态法、无约束索力优化法、有约束索力优化法和影响矩阵法［3-6］. 这些方法各有其优缺点与相应的适用条件，同时都有成功应用的案例［7-10］，然而针对具体的桥梁，是选用其中一种方法或多种方法的组合，仍然有待进一步深入研究确定. 本文结合某大跨径钢箱梁斜拉桥，提出一种基于Midas/Civil，综合利用影响矩阵调值原理，以结构整体弯矩最小为目标、部分关键节点的位移为约束条件的斜拉桥索力优化方法. 1 成桥索力优化理论 1.1 影响矩阵法影响矩阵法［11］不仅能解决线性调值计算问题，还能通过广义影响矩阵解决非线性的调值计算问题，因此该方法是在诸多斜拉桥索力调整方法基础上发展起来的一种理论完备的索力优化方法. 影响矩阵是由施调向量、影响向量和受调向量所组成的线性方程. 其计算原理为，当结构满足线性叠加原理， 则 [CX=D]， （1）式中： [C]为影响矩阵，[X]为施调向量， [D]为受调向量. 利用有限元软件，在单位荷载或者单位变形作用下，取斜拉桥关心截面的内力、应力、位移为受调向量[D]，斜拉索索力为施调向量[X]，通过影响矩阵[C]建立受调向量与施调向量的关系，将索力优化问题转化为式（1）. 1.2 最小弯矩法最小弯矩法［12］是以结构的弯矩平方和最小为目标函数，对斜拉桥合理索力进行求解和优化，其计算原理与最小弯曲应变能法相似. 而最小弯曲应变能法［13］是以结构的弯曲应变能作为目标函数，使结构弯曲应变能最小. 对于离散的杆系结构，弯曲应变能可写成： [U=i=1mLi2EiIiM2i]， （2）式中：m为结构单元总数；[Li]、[Ei]、[Ii]分别为第i号单元的杆件长度、材料弹性模量、截面惯性矩. 将式（2）改写成矩阵形式： [U={M}TBM]， （3）式中：[{M}]是弯矩向量； [B]为对角系数矩阵， [bii=Li/2EiIii=1，2，3，?，m]. （4）令调索前弯矩向量为[M0]，施调向量为[{T}]，结合式（1）和式（3）得： [U=C0+M0TBCT+TTCTBM0+] [TTCTBCT ，] （5）式中： [C]为影响矩阵，[C0]是与[T]无关的常数. 要使索力调整后结构弯曲应变能最小，则： [ U/ Ti=0 (i=1,2,3，?，l]）. （6）将式（5）代入式（6） ，写成矩阵形式： [CTBCT+CTBM0=0]. （7）仅当EI=1和式（7）中的[B]为单位矩阵时， 求解式（7）线性方程组即可得到结构弯矩平方和最小下的成桥索力［14］. 1.3 未知荷载系数法未知荷载系数法是计算满足某种约束条件和目标函数的优化的变量值，将每个要计算的索力定义为一个荷载工况，力的大小为单位荷载；在其它荷载工况（均布荷载）和单位力作用下进行结构分析；使用约束条件建立不平衡方程，计算满足不平衡方程的解就是所求的未知荷载系数. Midas/Civil中提供的解是使目标函数最小的解，目标函数的类型有变量的绝对值之和（[i=1nXi]）、变量的平方之和（[i=1nX2i]）、变量的绝对值中最大值（[MaxX1，X2，?，Xn]）等类型；约束条件可以是主塔和主梁位移、弯矩和索力的上下限，变量值为斜拉索索力. 2 优化方法与步骤优化是在满足所有约束条件的情况下，求出一组设计变量的值，使目标函数的值达到最优，其主要包括以下几方面：设计变量、约束条件和目标函数［15］. 其中目标函数可表达为： [minfx]， （8）式中：目标函数[fx]为结构弯矩平方和，设计变量[x={x1,x2,?xn}T]为斜拉索初张拉力. 约束条件的表现形式可写成 [gux0]或[gux0]， （9）式中：u=1，2，3，···，m，m为约束条件个数，约束条件可以是主塔和主梁位移、弯矩和索力的上下限等. 具体优化步骤如下：步骤1： 根据设计参数，采用Midas/Civil建立全桥有限元模型；步骤2： 定义主梁的恒载和斜拉索单位荷载的荷载工况及边界约束条件，然后输入恒载和单位荷载，进行计算；步骤3： 将斜拉索单位拉力和恒载进行组合，生成荷载组合工况；步骤4： 利用Midas/Civil中“未知荷载系数”功能，约束条件设置为：索塔塔顶水平位移、主梁位移小于某一限值，索塔与主梁弯矩小于某一限值，求解未知荷载系数，计算初始索力；步骤5： 根据步骤4计算所得未知荷载系数，将生成的系数矩阵导入Excel中，应用Excel的规划求解功能，建立以主梁和主塔最小弯矩为目标函数，主梁和主塔位移和弯矩为约束条件的优化函数方程求解最优索力；步骤6： 根据步骤5计算索力在Midas/ Civil建立荷载组合，重新计算，在结果中查看各限制条件是否满足设定目标，如果不满足，还可以使用Midas/Civil中自带调索功能进行微调，使主梁弯矩、主梁位移、主塔位移以及索力控制在合理的范围内. 3 工程应用 3.1 工程概况某三跨双索面钢箱梁斜拉桥，顺桥向沿中跨跨中对称，横桥向对称，斜拉索每塔16对，全桥共64对，箱梁每截断长15 m，立面布置如图1所示. 桥梁基本参数：中跨跨径600 m，边跨跨径213 m，塔高107.35 m；单根斜拉索截面面积6.02×10-5 m2，弹性模量1.95×105 MPa，加劲梁截面面积3.062 6 m2，弹性模量2.05×105 MPa. 3.2 计算模型利用Midas/Civil软件建立成桥有限元模型，主塔和主梁采用空间梁单元模拟，斜拉索采用等效的桁架单元模拟. 约束条件为塔墩固结，主梁和主塔采用弹性连接中的一般连接，顺桥方向和垂直方向固定，边墩处垂直方向固定，顺桥向释放，斜拉索与主梁采用弹性连接中的刚性连接进行连接，斜拉索与主塔上的节点直接连接. 全桥计算模型由552个单元和550个节点组成. 全桥有限元模型如图2所示. 3.3 初始成桥索力采用Midas/Civil中的未知荷载系数功能，建立约束条件：加劲梁弯矩≤1×105 kN·m，主梁竖向位移≤0.1 m，主塔弯矩≤5 000 kN·m，主塔水平位移≤0.03 m；求解能满足约束条件的荷载系数. 根据未知荷载系数计算斜拉索的初拉力，如表1所示（由于结构对称，表1中仅列出左岸主塔两侧斜拉索初拉力）. 从表1可知，主塔侧R1号索力远小于R2号索力值，其结果将导致塔梁交接处主梁绝大部分重力都分配到支座上，主梁负弯矩偏大，相应地梁、塔弯矩和位移都偏大. 因此，需对计算成桥索力进行优化. 3.4 成桥索力优化将计算初始索力形成的系数矩阵导入Excel中，利用Excel中规划求解功能，设定目标函数为主梁和主塔弯矩平方和最小，限值条件：索力2 000 kN~6 500 kN，主梁位移±0.06 m、部分主梁弯矩： ±8×104 kN·m，主塔弯矩± 3 000 kN·m. 执行规划求解，不断改变约束条件，得出较为合理的索力分布，然后转入Midas/Civil环境，进行索力微调，计算结果如表2所示. 利用优化后的索力计算斜拉桥主梁、主塔弯矩和位移，如表3所示. 由表3可知，优化后，主梁弯矩减小了36%，主塔弯矩减小了70%；主梁竖向最大位移减小了51%；主塔水平最大位移减小了74%. 为进一步说明优化效果，提取主梁、主塔的弯矩、 位移值等关键参数优化前后的值，绘制分布曲线，如图3~图5所示. 从图3~图5中可以看出，优化后索力比初始索力分布更为均匀，主梁位移较之前更接近目标线型；主塔弯矩明显减小70%，主塔顶水平位移减小74%. 4 结 语 1）针对索桥结构索力可调的特点，结合某钢箱梁斜拉桥，提出了以结构整体弯矩平方和最小为目标，关键节点位移为约束条件的索力优化方法. 2）借助Midas/Civil的调索功能与Excel的优化求解工具，能够方便、快捷地实现上述索力优化方法. 3）实际应用结果表明，本文所述的索力优化方法，能有效改善斜拉桥主梁和主塔的受力和变形，使成桥状态达到设计目标.