CN113468632B - 一种偏心活载作用下的悬索桥全桥响应的确定方法 - Google Patents

Abstract

本发明属于桥梁计算理论领域，特别涉及一种偏心活载作用下的悬索桥全桥响应的确定方法。该方法包括确定可全桥主要构件的内力及变形情况所需要的基本未知量；分析各跨的主缆响应；分析主梁响应；分析吊杆及桥塔响应；利用悬索桥上下游各节段主缆的两端的高度差与跨径闭合、主缆和吊杆的无应力长度守恒、主梁端部扭转角为0、主梁所受力与力矩平衡条件建立控制方程组；解得同时满足所有控制方程组的基本未知量的值；将基本未知量的值代入分析各跨的主缆响应表达式、分析主梁响应表达式、分析吊杆及桥塔响应表达式中，确定桥梁各构件的内力及变形。本方法，考虑了主梁的刚***移及吊杆倾斜，更贴近工程实际；规避了复杂计算，计算方便，便于应用。

Description

一种偏心活载作用下的悬索桥全桥响应的确定方法

技术领域

本发明属于桥梁计算理论领域，特别涉及一种偏心活载作用下的悬索桥全桥响应的确定方法。

背景技术

悬索桥凭借着跨越能力大、结构形式优美、受力性能优而逐渐称为大跨桥梁的最优选择。并且随着对悬索桥的设计理论研究的深入，通过材料特性的改善，结构尺寸的合理选用，悬索桥的技术难点，例如悬索桥的抗震抗风问题，悬索桥的结构过柔、空间缆索的设计及施工问题等正在逐一被克服。然而随着工业、交通运输业的发展，重载车辆数目增多且重量增大，重载车辆在运输过程中，一般以满载状态固定单向通行，导致桥梁上下游承受荷载不等，并且随着悬索桥向着跨径更长、桥面更宽的方向不断发展，悬索桥将承受较大的偏心活载，在这类活载作用下，悬索桥常见的扁平钢箱梁的竖向弯曲、扭转、畸变等变形问题更加突出，危及桥梁安全。近年来，由偏载引起的桥梁垮塌事故不胜枚举。对于悬索桥而言，在进行桥梁设计时一般为对称设计，桥梁上下游的主缆、吊杆的参数相同，若在长期的偏载作用下，悬索桥两侧的缆索体系受损情况不同，将缩短悬索桥在役寿命。因此，研究悬索桥在偏心活载作用下的全桥响应显得尤为必要。

近现代悬索桥的设计理论经历了弹性理论-挠度理论-有限位移理论的变革，目前，挠度理论与有限位移理论应用较多，在悬索桥静力分析的相关研究中，为充分理解活载作用下全桥响应求解方法的物理含义，现常基于挠度理论，采用解析法研究各项参数对全桥的影响。然而，一般的解析法中存在着以下几点共性问题：(1)未考虑活载的横向分布情况。活载对于悬索桥的静力分析中，一般认为荷载直接作用在桥梁纵向中轴线处。但若活载的偏载作用较大时，忽略横向分布导致的主梁扭转及全桥响应的不对称性，将大大削弱方法的精度；(2)考虑不对称活载时，忽略主梁的纵向位移。在实际工程中一般不限制主梁的纵向位移，当活载作用关于桥梁跨中位置不对称时，主梁会产生一定程度的纵向位移，进而带动吊杆发生倾斜，吊杆力不再垂直于桥面。当活载作用较大时，主梁的纵向位移较大，若不考虑将增大最终所得的响应误差；(3)将吊杆力视为均布的薄膜力，与实际不符；(4)由于考虑悬索桥的几何非线性导致最终计算过于繁琐。

发明内容

发明目的：针对上述问题，本发明提供了一种考虑活载横向分布、主梁纵向刚***移及吊杆倾斜等情况，计算简便，可以精确确定偏心活载作用下的全桥响应的确定方法。

本发明采用如下技术方案：

本发明所述的偏心活载作用下的悬索桥全桥响应的计算方法，包括如下步骤：

(1)确定可表达在偏心活载作用下全桥主要构件的内力及变形情况所需要的基本未知量；

该基本未知量包括：各跨主缆的悬链线方程参数、各节段悬链线的水平投影长度、主缆各节段水平分力、吊杆力、支反力、主梁整体水平向刚***移；

(2)分析各跨的主缆响应，用基本未知量表达：悬索桥上下游各节段主缆两端的高度差、无应力长度及主缆上的吊杆上吊点的坐标；

(3)分析主梁响应，用基本未知量表达：偏心活载作用下主梁各截面的挠度、横向扭转角及主梁上的吊杆下吊点的坐标；

(4)分析吊杆及桥塔响应，用基本未知量表达出上、下游吊杆倾角及无应力长度、桥塔各肢的偏移量；

(5)利用悬索桥上下游各节段主缆的两端的高度差与跨径闭合、主缆和吊杆的无应力长度守恒、主梁端部扭转角为0、主梁所受力与力矩平衡条件建立控制方程组；

(6)利用非线性规划方法求解步骤(5)建立的控制方程组，解得同时满足所有控制方程组的基本未知量的值；

(7)将基本未知量的值代入分析各跨的主缆响应表达式、分析主梁响应表达式、分析吊杆及桥塔响应表达式中，确定桥梁各构件的内力及变形。

本能发明所述的偏心活载作用下的悬索桥全桥响应的计算方法，

所述步骤(1)中的基本未知量表达式如下：

左边跨上下游主缆的悬链线方程参数a_u,L、a_d,L；主跨上下游左侧第一节段主缆的悬链线方程参数a_u,1、a_d,1；右边跨上下游主缆的悬链线方程参数a_u,R、a_d,R；

左边跨上下游主缆水平投影长度l_u,L、l_d,L；主跨上下游n+1个主缆节段水平投影长度l_u,1～l_u,n+1、l_d,1～l_d,n+1；右边跨上下游主缆水平投影长度l_u,R、l_d,R；

左边跨上下游主缆的水平分力H_u,L、H_d,L；主跨上下游左侧第一节段主缆的水平分力H_u,1、H_d,1；右边跨上下游主缆的水平分力H_u,R、H_d,R；

主跨上下游n根吊杆的吊杆力F_u,1～F_u,n、F_d,1～F_d,n；

主梁左端上下游支反力F_u,L、F_u,R；主梁右端上下游支反力F_d,L、F_d,R；

主梁整体水平向刚***移v；

一共有4n+23个基本未知量(n为主跨一侧索面布置的吊杆数)，下标u和d分别表示桥梁上游一侧和下游一侧的参数。

本发明所述的偏心活载作用下的悬索桥全桥响应的计算方法，

所述步骤(2)中表达上下游各节段主缆两端高度差、无应力长度以及主缆上的吊杆上吊点的坐标如下：

(2.1)每一段悬链线以其左端点为原点，建立局部坐标系，x轴水平向右，y轴竖直向下，在该坐标系下，各参数表述如下：

主跨上游一侧第i段主缆可视为悬链线，其左右端点的高度差Δy_u,i可以表示为：

式中，c_u,i＝-H_u,i/q，H_u,i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,i为主跨上游侧第i段悬链线方程参数；l_u,i为主跨上游侧第i段悬链线水平投影长度(m)；cosh为数学符号；

主跨上游一侧第i段主缆无应力长度S_u,i表示为：

式中，c_u,i＝-H_u,i/q，H_u,i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh为数学符号；a_u,i为主跨上游侧第i段悬链线方程参数；

主跨下游一侧第i段主缆左右端点的高度差Δy_d,i可以表示为：

式中，c_d,i＝-H_d,i/q，H_d,i为恒载及活载作用主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,i为主跨下游侧第i段悬链线方程参数；l_d,i为主跨下游侧第i段悬链线水平投影长度(m)；cosh为数学符号；

主跨下游一侧第i段主缆无应力长度S_d,i可以表示为：

式中，c_d,i＝-H_d,,i//q，H_d,i为恒载及活载作用主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_d,i为主跨下游侧第i段悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；a_d,i为主跨下游侧第i段悬链线方程参数；sinh为数学符号；

(2.2)在恒载和偏心活载共同作用下，上游侧第i根吊杆的上吊点坐标(x_cu,i,y_cu,i,z_cu,i)表达为：

z_cu,i＝b

式中，Δx_u,L表示左端桥塔上游侧的塔肢顶部的偏移量(m)；b为各吊杆至主梁中心的横向距离(m)；H_B表示主跨主缆左侧IP点标高(m)；H_E表示主梁左侧截面中心的标高(m)；l_u,i为主跨上游侧第i段主缆的水平投影长度(m)；Δy_u,j表示主跨上游侧第j段主缆悬链线两端高度差(m)；

下游侧第i根吊杆的上吊点坐标(x_cd,i,y_cd,i,z_cd,i)表达为：

z_cd,i＝-b

式中，Δx_d,L表示左端桥塔下游侧的塔肢顶部的偏移量(m)；b为各吊杆至主梁中心的横向距离(m)；H_B表示主跨主缆左侧IP点标高(m)；H_E表示主梁左侧截面中心的标高(m)；l_d,i为主跨下游侧第i段主缆的水平投影长度(m)；Δy_d,j表示主跨下游侧第j段主缆悬链线两端高度差(m)；

(2.3)恒载和活载共同作用下，主跨上游侧的相邻节段主缆的水平分力H和悬链线方程参数a之间的递推关系式为:

式中，F_u,i为上游侧第i根吊杆的吊杆力(kN)；

为上游侧第i根吊杆力与竖直平面的夹角(rad)；a_u,i为主跨上游侧第i段悬链线方程参数；c_u,i＝-H_u,,i/q，H_u,_i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_u,i为主跨上游侧第i段主缆的水平投影长度(m)；arcsinh为数学符号；

主跨下游侧的相邻节段主缆的水平分力H和悬链线方程参数a之间的递推关系式为：

式中，F_d,i为下游侧第i根吊杆的吊杆力(kN)；

为下游侧第i根吊杆力与竖直平面的夹角(rad)；a_d,i为主跨下游侧第i段悬链线方程参数；c_d,i＝-H_d,,i/q，H_d,_i为恒载及活载作用下主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_d,i为主跨下游侧第i段主缆的水平投影长度(m)；arcsinh为数学符号；

(2.4)上游侧左边跨的主缆段，以左锚固点为原点，纵桥向为x轴，向右为正方向，竖直向下为y轴，建立局部坐标系；

上游左边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_u,L和无应力长度S_u,L可表示为：

式中，c_u,L＝-H_u,L/q，H_u,L为恒载及活载作用左边跨上游主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,L为上游左边跨主缆悬链线方程参数；l_u,L为上游左边跨主缆水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

右边跨段的主缆以右锚固点为原点，纵桥向为x轴，向右为正，竖直向下为y轴建立局部坐标系，则上游右边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_u,R和无应力长度S_u,R可表示为:

式中，c_u,R＝-H_u,R/q，H_u,R为恒载及活载作用上游右边跨主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,R为上游右边跨悬链线方程参数；l_u,R为上游右边跨悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号；

以下游侧左边跨的主缆段的左锚固点为原点，纵桥向为x轴，向右为正方向，竖直向下为y轴，建立局部坐标系；下游左边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_d,L和无应力长度S_d,L可表示为：

式中，c_d,L＝-H_d,L/q，H_d,L为恒载及活载作用左边跨下游主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,L为下游左边跨主缆悬链线方程参数；l_d,L为下游左边跨主缆水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号；

右边跨段的主缆以右锚固点为原点，纵桥向为x轴，向右为正，竖直向下为y轴建立局部坐标系，则下游右边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_d,R和无应力长度S_d,R可表示为:

式中，c_d,R＝-H_d,R/q，H_d,R为恒载及活载作用下游右边跨主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,R为下游右边跨悬链线方程参数；l_d,R为下游右边跨悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

本发明所述偏心活载作用下的悬索桥全桥响应的计算方法，所述步骤(3)中分析主梁响应，用基本未知量表达主梁的挠度和横向扭转角以及主梁上的吊杆下吊点的坐标如下：

(3.1)以主梁左端截面形心为原点，纵桥向为x轴，竖直向下为y轴，横向为z轴建立局部坐标系；主梁受到两组吊杆力F_u,i、F_d,i(1≤i≤n)以及四个支座反力F_u,L、F_u,R、F_d,L、F_d,R的作用。

在偏心活载作用下，主梁可被n根吊杆划分为n+1个区段，每一段主梁的挠度w(x)为：

ΔF_u,L＝F_u,L-F’_u,L

ΔF_d,L＝F_d,L-F’_d,L

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力在受到活载作用前后的变化量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；F'_u,i、F'_d,i分别表示上游、下游第i根吊杆在恒载作用下的吊杆力(kN)；F'_u,L、F'_d,L分别表示上游、下游左支座在恒载作用下的支座反力(kN)；E_bI_b表示主梁抗弯刚度(kN·m²)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)；p为偏心活载集度(kN/m)；x表示所求截面与该区段主梁最左端的距离(m)；

为下游侧第i根吊杆力与竖直平面的夹角(rad)；t为积分所取微元至主梁左端的距离(m)；dt为积分变量；j为所求截面左侧距离最近的吊杆编号；

(3.2)距离主梁左端距离为x的截面Ⅰ对于左侧距该截面最近的吊杆右截面Ⅱ的相对转角为：

式中，GI_P为主梁的自由扭转刚度(kN·m/rad)；ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力在受到活载作用前后的变化量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；p为偏心活载集度(kN/m)；b为吊杆横桥向间距的一半(m)；b'为支座横桥向间距的一半(m)；e为偏心活载作用位置与桥梁纵轴线的横向距离(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)；x表示所求截面与该区段主梁最左端的距离(m)；t为积分所取微元至主梁左端的距离(m)；t为积分所取微元至主梁左端的距离(m)；dt为积分变量；j为所求截面左侧距离最近的吊杆编号；x_j表示恒载条件下第j根吊杆下吊点至原点的纵桥向距离(m)；

当x＝x_m+1时，可求得相邻两根吊杆截面的相对扭转角，根据θ(0)＝0，通过叠加即可知主梁各截面扭转角的值。

(3.3)主梁上游侧第i根吊杆的下吊点坐标(x_bu,i,y_bu,i,z_bu,i)可以表示为：

x_bu,i＝x_i+v

y_bu,i＝-w(x_i)-bθ(x_i)

z_bu,i＝b

式中，v表示主梁整体水平方向的刚***移(m)；w表示偏心活载作用下主梁发生的竖向挠度(m)；θ表示主梁发生的扭转角(rad)；b为吊杆横桥向间距的一半(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)。

主梁下游侧第i根吊杆的下吊点坐标(x_bd,i,y_bd,i,z_bd,i)可以表示为

x_bd,i＝x_i+v

y_bd,i＝-w(x_i)+bθ(x_i)

z_bd,i＝-b

式中，v表示主梁整体水平方向的刚***移(m)；w表示偏心活载作用下主梁发生的竖向挠度(m)；θ表示主梁发生的扭转角(rad)；b为吊杆横桥向间距的一半(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)。

本发明所述的偏心活载作用下的悬索桥全桥响应的计算方法，所述步骤(4)中表达上下游吊杆倾角及无应力长度、桥塔各肢的偏移量的具体步骤如下：

(4.1)根据步骤(2.2)中建立的上游侧第i根吊杆的上吊点坐标(x_cu,i,y_cu,i,z_cu,i)与步骤(3.3)中建立的上游侧第i根吊杆的下吊点坐标(x_bu,i,y_bu,i,z_bu,i)计算，在偏心活载作用下，上游侧第i根吊杆的倾角

可以表示为：

根据步骤(2.2)中建立的下游侧第i根吊杆的上吊点坐标(x_cd,i,y_cd,i,z_cd,i)与步骤(3.3)中建立的下游侧第i根吊杆的下吊点坐标(x_bd,i,y_bd,i,z_bd,i)计算，下游侧第i根吊杆的倾角

可以表示为：

(4.2)上游侧第i根吊杆的无应力长度S_su,i可以表示为：

式中，E_sA_s为吊杆抗拉刚度(kN)；S’_su,i为上游一侧第i根吊杆在恒载作用下的无应力长度(m)；F_u,i为上游侧第i根吊杆的吊杆力；x_cu,i、x_bu,i分别为上游第i根吊杆上吊点、下吊点的x坐标；y_cu,i、y_bu,i分别为上游第i根吊杆上吊点、下吊点的y坐标。

下游侧第i根吊杆的无应力长度S_sd,i可以表示为：

式中，E_sA_s为吊杆抗拉刚度(kN)；S’_sd,i为下游一侧第i根吊杆在恒载作用下的无应力长度(m)；F_d,i为下游侧第i根吊杆的吊杆力；x_cd,i、x_bd,i分别为下游第i根吊杆上吊点、下吊点的x坐标；y_cd,i、y_bd,i分别为下游第i根吊杆上吊点、下吊点的y坐标。

(4.3)左桥塔上游、下游塔肢顶部的偏移量Δx_u,L、Δx_d,L分别为：

式中，k为单肢桥塔的刚度系数(kN/m)；k_ξ为整个桥塔的抗扭刚度系数(kN/rad)；d为每肢桥塔距离桥轴线距离(m)；H_u,L、H_d,L分别为上下游左边跨主缆水平分力(kN)；H_u,1、H_d,1分别为上下游主跨左侧第一段主缆悬链线水平分力(kN)；

右桥塔上游、下游塔肢顶部的偏移量Δx_u,R、Δx_d,R分别为：

式中，k为单肢桥塔的刚度系数(kN/m)；k_ξ为整个桥塔的抗扭刚度系数(kN/rad)；d为每肢桥塔距离桥轴线距离(m)；H_u,R、H_d,R分别为上下游右边跨主缆水平分力(kN)；H_u,n+1、H_d,n+1分别为上下游主跨最右侧节段(第n+1段)主缆悬链线水平分力(kN)。

本发明所述的偏心活载作用下的悬索桥全桥响应的计算方法，所述步骤(5)中建立控制方程组的具体步骤如下：

(5.1)在恒载及活载共同作用时悬索桥上游侧左边跨、主跨及右边跨主缆的高度差

闭合的控制方程表达如下：

Δy_u,L＝H_A-H_B

Δy_u,R＝H_D-H_C

式中，H_A表示主缆左锚固点标高(m)；H_B表示主跨主缆左侧IP点标高(m)；H_C表示主跨主缆右侧IP点标高(m)；H_D表示主缆右锚固点标高(m)；

在恒载及活载共同作用时下游侧左边跨、主跨及右边跨主缆的高度差

闭合的控制方程表达如下：

Δy_d,L＝H_A-H_B

Δy_d,R＝H_D-H_C

式中，H_A表示主缆左锚固点标高(m)；H_B表示主跨主缆左侧IP点标高(m)；H_C表示主跨主缆右侧IP点标高(m)；H_D表示主缆右锚固点标高(m)。

在恒载及活载共同作用时悬索桥上游侧左边跨、主跨及右边跨主缆的跨径

闭合的控制方程表达如下：

l_u,L＝L₁+Δx_u,L

l_u,R＝L₃-Δx_u,R

式中，L₁、L₂、L₃分别表示悬索桥左边跨、主跨和右边跨的跨径(m)；Δx_u,L为左桥塔上游侧塔肢顶部纵桥向位移(m)；Δx_u,R为右桥塔上游侧塔肢顶部纵桥向位移(m)。

在恒载及活载共同作用时下游侧左边跨、主跨及右边跨主缆的跨径

闭合的控制方程表达如下：

l_d,L＝L₁+Δx_d,L

l_d,R＝L₃-Δx_d,R

式中，L₁、L₂、L₃分别表示悬索桥左边跨、主跨和右边跨的跨径(m)；Δx_d,L为左桥塔下游侧塔肢顶部纵桥向位移(m)；Δx_d,R为右桥塔下游侧塔肢顶部纵桥向位移(m)。

(5.2)在恒载及活载共同作用时上游侧左边跨、主跨各节段及右边跨主缆无应力长度(S_cu,L,S_cu,i,S_cu,R)守恒的控制方程表达如下：

S_cu,L＝S’_cu,L

S_cu,i＝S’_cu,i,(1≤i≤n+1)

S_cu,R＝S’_cu,R

式中，S'_cu,i为主跨上游侧第i节段主缆在恒载作用下的无应力长度(m)；S'_cu,L为左边跨上游侧主缆在恒载作用下的无应力长度(m)；S'_cu,R为右边跨上游侧主缆在恒载作用下的无应力长度(m)。

在恒载及活载共同作用时下游侧左边跨、主跨各节段及右边跨主缆无应力长度(S_cd,L,S_cd,i,S_cd,R)守恒的控制方程表达如下：

S_cd,L＝S’_cd,L

S_cd,i＝S’_cd,i,(1≤i≤n+1)

S_cd,R＝S’_cd,R

式中，S'_cd,i为主跨下游侧第i节段主缆在恒载作用下的无应力长度(m)；S'_cd,L为左边跨下游侧主缆在恒载作用下的无应力长度(m)；S'_cd,R为右边跨下游侧主缆在恒载作用下的无应力长度(m)。

在恒载及活载共同作用时上游侧各吊杆无应力长度(S_su,i)守恒的控制方程表达如下：

S_su,i＝S’_su,i,(1≤i≤n)

式中，S′_su,i为上游一侧第i根吊杆在恒载作用时的无应力长度(m)。

在恒载及活载共同作用时下游侧各吊杆无应力长度(S_sd,i)守恒的控制方程表达如下：

S_sd,i＝S’_sd,i,(1≤i≤n)

式中，S′_sd,i为下游一侧第i根吊杆在恒载作用时的无应力长度(m)。

(5.3)主梁端部扭转角为0控制方程表达如下：

θ(L₂)＝0

式中，L₂为悬索桥主跨跨径(m)。

(5.4)主梁y向受力(∑Y_i)平衡控制方程表达如下：

ΔF_u,R＝F_u,R-F’_u,R

ΔF_d,R＝F_d,R-F’_d,R

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；F_u,R、F_d,R分别表示上游、下游右支座在恒载和活载共同作用时的支座反力(kN)；F'_u,R、F'_d,R分别表示上游、下游右支座在恒载作用时的支座反力(kN)。

主梁x向受力(∑X_i)平衡控制方程表达如下：

式中ΔF_xu,i、ΔF_xd,i分别表示上游侧、下游侧第i根吊杆力的水平分力的改变量(kN)；F_u,_i、F_d,i分别表示上游、下游右支座在恒载和活载共同作用时的吊杆力(kN)；

为上游、下游侧吊杆力与竖直平面的夹角(rad)；

主梁绕x轴力矩(∑M_x)平衡控制方程表达如下：

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；b为吊杆横桥向间距的一半(m)；b'为支座横桥向间距的一半(m)。

主梁绕z轴力矩(∑M_z)平衡控制方程表达如下：

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；x_i表示恒载条件下，第i根吊杆下吊点至原点的纵桥向距离(m)；L₂为悬索桥主跨跨径(m)；

本发明所述偏心活载作用下的悬索桥全桥响应的计算方法，所述的步骤(3.1)中在第j-1根吊杆与第j根吊杆间的第j段主梁挠度表达式中由于不定积分产生的常数项系数与一次项系数C_j,₀、C_j,₁以及第j+1段的主梁挠度表达式中由于不定积分产生的常数项系数与一次项系数C_j+1,0、C_j+1,1的递推关系式如下：

式中，ΔF_yu,j、ΔF_yd,j分别表示上游侧、下游侧第j根吊杆力的竖向分力的改变量(kN)；x_j表示恒载条件下，第j根吊杆下吊点至原点的纵桥向距离(m)；

由主梁左端挠度w_L＝0，求得第1根吊杆左侧的主梁挠度表达式，即第1段主梁挠度中由于不定积分产生的常数项系数C_1,0，进而递推出所有段挠度表达式中的常数项系数C_i,0，包括C_n,0。

主梁右端挠度w_R＝0,可以求得第n根吊杆右侧的主梁挠度表达式，即第n+1段吊杆挠度表达式由于不定积分产生的常数项系数与一次项系数C_n+1,0、C_n+1,1的关系式，进而得到C_n+1,1，则可以反推出所有段挠度表达式中的一次项系数C_i,1。

将所求得的常数项系数及一次项系数代入相对应的主梁区段的挠度表达式，即可推得主梁任一截面的挠度。

本发明所述偏心活载作用下的悬索桥全桥响应的计算方法，步骤(6)中利用非线性规划求解控制方程组，解得同时满足所有控制方程组的基本未知量的值的具体步骤如下：

利用步骤(5)中建立的4n+23个控制方程，将步骤(2)至步骤(4)中用基本未知量表示的各项参数代入控制方程组中，通过移项使得等式右边为0，等式左边视为4n+23个函数值W_i；利用非线性GRG方法对下式进行规划求解，即可得出悬索桥在偏心活载作用下的响应；

本发明所述偏心活载作用下的悬索桥全桥响应的计算方法，将步骤(6)中所获得的基本未知量的值代入(2.1)及(2.4)中的各段主缆高度差表达式、(2.2)中吊杆上吊点坐标表达式和(2.3)中主缆水平分力的递推关系式中，确定在恒载和活载作用下主缆的内力和变形响应；

将步骤(6)中所获得的基本未知量的值代入(3.1)中的挠度表达式、(3.2)中相对转角的表达式和(3.3)中吊杆下吊点坐标表达式，确定在恒载和活载作用下主梁的内力和变形响应；

将步骤(6)中所获得的基本未知量的值代入(4.1)中的吊杆倾角表达式可确定在恒载和活载作用下吊杆的倾斜角度，代入(4.3)中各塔肢顶部桥塔偏移量的表达式可确定桥塔在恒载和活载作用下的变形响应。

有益效果

本发明提供的偏心活载作用下的悬索桥全桥响应的确定方法，将吊杆力视为离散集中力，考虑了主梁的刚***移及吊杆倾斜，所用的假设条件较少，更贴近工程实际；提出了利用非线性规划求解，可全面求解在不同位置布置偏心荷载时主梁、主缆及吊杆的响应，规避了复杂计算，计算方便，便于应用。

附图说明

图1为具体实施例中的恒载和活载共同作用下的主跨主缆示意图；

图2为具体实施例中的悬索桥全桥布置图；

图3为具体实施例中的恒载和活载共同作用下的主梁受力情况示意图；

图4为具体实施例中的吊杆下吊点的受力情况示意图；

图5为具体实施例中的恒载和活载共同作用下桥塔塔顶的示意图。

具体实施方式

为使本发明实施例的目的和技术方案更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1至图5所示：本发明所述的一种偏心活载作用下的悬索桥全桥响应的计算方法，该方法包括以下步骤：

(1)确定可表达在偏心活载作用下全桥主要构件的内力及变形情况所需要的基本未知量；

(2)分析各跨的主缆响应，用基本未知量表达出悬索桥上下游各节段主缆两端的高度差、无应力长度以及主缆上的吊杆上吊点的坐标；

(3)分析主梁响应，用基本未知量表达出偏心活载作用下主梁各截面的挠度和横向扭转角以及主梁上的吊杆下吊点的坐标；

(4)分析吊杆及桥塔响应，用基本未知量表达出上下游吊杆倾角及无应力长度、桥塔各肢的偏移量；

(5)利用悬索桥主缆的高差和跨径闭合、主缆和吊杆的无应力长度守恒、主梁端部扭转角为0、主梁所受力与力矩平衡等条件建立控制方程组。

(6)利用非线性规划方法求解控制方程组，解得同时满足所有控制方程组的基本未知量的值；

(7)将基本未知量的值代入上述各响应的表达式，确定桥梁各构件的内力及变形。

具体实施方式如下：

第一步：确定可表达在偏心活载作用下全桥主要构件的内力及变形情况所需要的基本未知量。所需要的基本未知量包括：

左边跨上下游主缆的悬链线方程参数a_u,L、a_d,L；主跨上下游左侧第一节段主缆的悬链线方程参数a_u,1、a_d,1；右边跨上下游主缆的悬链线方程参数a_u,R、a_d,R；

左边跨上下游主缆水平投影长度l_u,L、l_d,L；主跨上下游n+1个主缆节段水平投影长度l_u,1～l_u,n+1、l_d,1～l_d,n+1；右边跨上下游主缆水平投影长度l_u,R、l_d,R；

左边跨上下游主缆的水平分力H_u,L、H_d,L；主跨上下游左侧第一节段主缆的水平分力H_u,1、H_d,1；右边跨上下游主缆的水平分力H_u,R、H_d,R；

主跨上下游n根吊杆的吊杆力F_u,1～F_u,n、F_d,1～F_d,n；

主梁左端上下游支反力F_u,L、F_u,R；主梁右端上下游支反力F_d,L、F_d,R；

主梁整体水平向刚***移v；

一共有4n+23个基本未知量(n为主跨一侧索面布置的吊杆数)，下标u和d分别表示桥梁上游一侧和下游一侧的参数。

第二步：分析各跨的主缆响应，用基本未知量表达出悬索桥上下游各节段主缆两端的高度差、无应力长度以及主缆上的吊杆上吊点的坐标。

(1)分析恒载及活载共同作用下的主跨主缆(如图1所示)，主跨上游一侧第i段主缆悬链线左右端点的高差Δy_u,i可以表示为：

式中，c_u,i＝-H_u,i/q，H_u,i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,i为主跨上游侧第i段悬链线方程参数；l_u,i为主跨上游侧第i段悬链线水平投影长度(m)；cosh为数学符号。

主跨上游一侧第i段主缆无应力长度S_u,i可以表示为：

式中，c_u,i＝-H_u,i/q，H_u,i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh为数学符号；a_u,i为主跨上游侧第i段悬链线方程参数。

主跨下游一侧第i段主缆左右端点的高差Δy_d,i可以表示为：

式中，c_d,i＝-H_d,i/q，H_d,i为恒载及活载作用主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,i为主跨下游侧第i段悬链线方程参数；l_d,i为主跨下游侧第i段悬链线水平投影长度(m)；cosh为数学符号。

主跨下游一侧第i段主缆无应力长度S_d,i可以表示为：

式中，c_d,i＝-H_d,i/q，H_d,i为恒载及活载作用主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_d,i为主跨下游侧第i段悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；a_d,i为主跨下游侧第i段悬链线方程参数；sinh为数学符号。

(2)在恒载和偏心活载共同作用下，上游侧第i根吊杆的上吊点坐标(x_cu,i,y_cu,i,z_cu,i)可以表达为

z_cu,i＝b

式中，Δx_u,L表示左端桥塔上游侧的塔肢顶部的偏移量(m)；b为各吊杆至主梁中心的横向距离(m)；H_B表示主跨主缆左侧IP点标高(m)；H_E表示主梁左侧截面中心的标高(m)；l_u,i为主跨上游侧第i段主缆的水平投影长度(m)；Δy_u,j表示主跨上游侧第j段主缆悬链线两端高度差(m)。

下游侧第i根吊杆的上吊点坐标(x_cd,i,y_cd,i,z_cd,i)可以表达为

z_cd,i＝-b

式中，Δx_d,L表示左端桥塔下游侧的塔肢顶部的偏移量(m)；b为各吊杆至主梁中心的横向距离(m)；H_B表示主跨主缆左侧IP点标高(m)；H_E表示主梁左侧截面中心的标高(m)；l_d,i为主跨下游侧第i段主缆的水平投影长度(m)；Δy_d,j表示主跨下游侧第j段主缆悬链线两端高度差(m)。

(3)恒载和活载共同作用下，主跨上游侧的相邻节段主缆的水平分力H和悬链线方程参数a之间的递推关系式为

式中，F_u,i为上游侧第i根吊杆的吊杆力(kN)；

为上游侧第i根吊杆力与竖直平面的夹角(rad)；a_u,i为主跨上游侧第i段悬链线方程参数；c_u,i＝-H_u,i/q，H_u,_i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_u,i为主跨上游侧第i段主缆的水平投影长度(m)；arcsinh为数学符号。

主跨下游侧的相邻节段主缆的水平分力H和悬链线方程参数a之间的递推关系式为

式中，F_d,i为下游侧第i根吊杆的吊杆力(kN)；

为下游侧第i根吊杆力与竖直平面的夹角(rad)；a_d,i为主跨下游侧第i段悬链线方程参数；c_d,i＝-H_d,i/q，H_d,i为恒载及活载作用下主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_d,i为主跨下游侧第i段主缆的水平投影长度(m)；arcsinh为数学符号。

(4)边跨段主缆不受吊杆力作用，仅在自重作用下成悬链线形(如图2所示)。上游左边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_u,L和无应力长度S_u,L可表示为

式中，c_u,L＝-H_u,,L/q，H_u,L为恒载及活载作用左边跨上游主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,L为上游左边跨主缆悬链线方程参数；l_u,L为上游左边跨主缆水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

上游右边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_u,R和无应力长度S_u,R可表示为：

式中，c_u,R＝-H_u,R/q，H_u,R为恒载及活载作用上游右边跨主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,R为上游右边跨悬链线方程参数；l_u,R为上游右边跨悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

下游左边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_d,L和无应力长度S_d,L可表示为：

式中，c_d,L＝-H_d,L/q，H_d,L为恒载及活载作用左边跨下游主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,L为下游左边跨主缆悬链线方程参数；l_d,L为下游左边跨主缆水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

下游右边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_d,R和无应力长度S_d,R可表示为:

式中，c_d,R＝-H_d,R/q，H_d,R为恒载及活载作用下游右边跨主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,R为下游右边跨悬链线方程参数；l_d,R为下游右边跨悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

第三步：分析主梁响应，用基本未知量表达主梁的挠度和横向扭转角以及主梁上的吊杆下吊点的坐标。

(1)分析主梁受力情况(如图3所示)，主梁受到两组吊杆力F_u,i、F_d,i(1≤i≤n)以及四个支座反力F_u,L、F_u,R、F_d,L、F_d,R的作用。在仅受恒载作用下，全桥处于初始静力平衡阶段，并且主梁上各点位于设计标高，在本发明在分析主梁响应部分不考虑由于恒载以及该状态下的吊杆力等作用，所涉及到的作用与响应都是在恒载作用基础上的变化量。在偏心活载作用下，主梁可被n根吊杆划分为n+1个区段，每一段主梁的挠度表达式为

ΔF_u,L＝F_u,L-F’_u,L

ΔF_d,L＝F_d,L-F’_d,L

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力在受到活载作用前后的变化量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；F'_u,i、F'_d,i分别表示上游、下游第i根吊杆在恒载作用下的吊杆力(kN)；F'_u,L、F'_d,L分别表示上游、下游左支座在恒载作用下的支座反力(kN)；E_bI_b表示主梁抗弯刚度(kN·m²)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)；p为偏心活载集度(kN/m)；x表示所求截面与该区段主梁最左端的距离(m)；

为下游侧第i根吊杆力与竖直平面的夹角(rad)；t为积分所取微元至主梁左端的距离(m)；dt为积分变量；j为所求截面左侧距离最近的吊杆编号。

在第j-1根吊杆与第j根吊杆间的第j段主梁挠度表达式中由于不定积分产生的常数项系数与一次项系数C_j,₀、C_j,₁以及第j+1段的主梁挠度表达式中由于不定积分产生的常数项系数与一次项系数C_j+1,0、C_j+1,1的递推关系式如下：

式中，ΔF_yu,j、ΔF_yd,j分别表示上游侧、下游侧第j根吊杆力的竖向分力的改变量(kN)；x_j表示恒载条件下，第j根吊杆下吊点至原点的纵桥向距离(m)。

由主梁左端挠度w_L＝0，求得第1根吊杆左侧的主梁挠度表达式，即第1段主梁挠度中由于不定积分产生的常数项系数C_1,0，进而递推出所有段挠度表达式中的常数项系数C_i,0，包括C_n,0。此外，主梁右端挠度w_R＝0,可以求得第n根吊杆右侧的主梁挠度表达式，即第n+1段吊杆挠度表达式由于不定积分产生的常数项系数与一次项系数C_n+1,0、C_n+1,1的关系式，进而得到C_n+1,1，则可以反推出所有段挠度表达式中的一次项系数C_i,1。

将所求得的常数项系数及一次项系数代入相对应的主梁区段的挠度表达式，即可推得主梁任一截面的挠度。

(2)距离主梁左端距离为x的截面Ⅰ对于左侧距该截面最近的吊杆右截面Ⅱ的相对扭转角为

式中，GI_P为主梁的自由扭转刚度(kN·m/rad)；ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力在受到活载作用前后的变化量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；p为偏心活载集度(kN/m)；b为吊杆横桥向间距的一半(m)；b'为支座横桥向间距的一半(m)；e为偏心活载作用位置与桥梁纵轴线的横向距离(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)；x表示所求截面与该区段主梁最左端的距离(m)；t为积分所取微元至主梁左端的距离(m)；dt为积分变量；j为所求截面左侧距离最近的吊杆编号；x_j表示恒载条件下第j根吊杆下吊点至原点的纵桥向距离(m)。

当x＝x_m+1时，可求得相邻两根吊杆截面的相对扭转角，根据θ(0)＝0，通过叠加即可知主梁各截面扭转角的值。

(3)主梁上游侧第i根吊杆的下吊点坐标(x_bu,i,y_bu,i,z_bu,i)可以表示为

x_bu,i＝x_i+v

y_bu,i＝-w(x_i)-bθ(x_i)

z_bu,i＝b

式中，v表示主梁整体水平方向的刚***移(m)；w表示偏心活载作用下主梁发生的竖向挠度(m)；θ表示主梁发生的扭转角(rad)；b为吊杆横桥向间距的一半(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)。

主梁下游侧第i根吊杆的下吊点坐标(x_bd,i,y_bd,i,z_bd,i)可以表示为

x_bd,i＝x_i+v

y_bd,i＝-w(x_i)+bθ(x_i)

z_bd,i＝-b

式中，v表示主梁整体水平方向的刚***移(m)；w表示偏心活载作用下主梁发生的竖向挠度(m)；θ表示主梁发生的扭转角(rad)；b为吊杆横桥向间距的一半(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)。

第四步：分析吊杆及桥塔响应，用基本未知量表达出上下游吊杆倾角及无应力长度、桥塔各肢的偏移量。

(1)在偏心活载作用下，分析第i根吊杆下吊点(如图4所示)，根据步骤(2.2)中建立的上游侧第i根吊杆的上吊点坐标(x_cu,i,y_cu,i,z_cu,i)与步骤(3.3)中建立的上游侧第i根吊杆的下吊点坐标(x_bu,i,y_bu,i,z_bu,i)计算，在偏心活载作用下，上游侧第i根吊杆的倾角

可以表示为

根据步骤(2.2)中建立的下游侧第i根吊杆的上吊点坐标(x_cd,i,y_cd,i,z_cd,i)与步骤(3.3)中建立的下游侧第i根吊杆的下吊点坐标(x_bd,i,y_bd,i,z_bd,i)计算，下游侧第i根吊杆的倾角

可以表示为

(2)上游侧第i根吊杆的无应力长度S_su,i可以表示为

式中，E_sA_s为吊杆抗拉刚度(kN)；S’_su,i为上游一侧第i根吊杆在恒载作用下的无应力长度(m)；F_u,i为上游侧第i根吊杆的吊杆力；x_cu,i、x_bu,i分别为上游第i根吊杆上吊点、下吊点的x坐标；y_cu,i、y_bu,i分别为上游第i根吊杆上吊点、下吊点的y坐标。

下游侧第i根吊杆的无应力长度S_sd,i可以表示为

式中，E_sA_s为吊杆抗拉刚度(kN)；S’_sd,i为下游一侧第i根吊杆在恒载作用下的无应力长度(m)；F_d,i为下游侧第i根吊杆的吊杆力；x_cd,i、x_bd,i分别为下游第i根吊杆上吊点、下吊点的x坐标；y_cd,i、y_bd,i分别为下游第i根吊杆上吊点、下吊点的y坐标。

(3)左桥塔上游、下游塔肢顶部的偏移量Δx_u,L、Δx_d,L(如图5所示)分别为

式中，k为单肢桥塔的刚度系数(kN/m)；k_ξ为整个桥塔的抗扭刚度系数(kN/rad)；d为每肢桥塔距离桥轴线距离(m)；H_u,L、H_d,L分别为上下游左边跨主缆水平分力(kN)；H_u,1、H_d,1分别为上下游主跨左侧第一段主缆悬链线水平分力(kN)。

右桥塔上游、下游塔肢顶部的偏移量的偏移量Δx_u,R、Δx_d,R分别为

式中，k为单肢桥塔的刚度系数(kN/m)；k_ξ为整个桥塔的抗扭刚度系数(kN/rad)；d为每肢桥塔距离桥轴线距离(m)；H_u,R、H_d,R分别为上下游右边跨主缆水平分力(kN)；H_u,n+1、H_d,n+1分别为上下游主跨最右侧节段(第n+1段)主缆悬链线水平分力(kN)。

第五步：建立控制方程组。

(1)在恒载及活载共同作用时悬索桥上游侧左边跨、主跨及右边跨主缆的高度差

闭合的控制方程表达如下：

Δy_u,L＝H_A-H_B

Δy_u,R＝H_D-H_C

式中，H_A表示主缆左锚固点标高(m)；H_B表示主跨主缆左侧IP点标高(m)；H_C表示主跨主缆右侧IP点标高(m)；H_D表示主缆右锚固点标高(m)。

在恒载及活载共同作用时下游侧左边跨、主跨及右边跨主缆的高度差

闭合的控制方程表达如下：

Δy_d,L＝H_A-H_B

Δy_d,R＝H_D-H_C

式中，H_A表示主缆左锚固点标高(m)；H_B表示主跨主缆左侧IP点标高(m)；H_C表示主跨主缆右侧IP点标高(m)；H_D表示主缆右锚固点标高(m)。

在恒载及活载共同作用时悬索桥上游侧左边跨、主跨及右边跨主缆的跨径

闭合的控制方程表达如下：

l_u,L＝L₁+Δx_u,L

l_u,R＝L₃-Δx_u,R

式中，L₁、L₂、L₃分别表示悬索桥左边跨、主跨和右边跨的跨径(m)；Δx_u,L为左桥塔上游侧塔肢顶部纵桥向位移(m)；Δx_u,R为右桥塔上游侧塔肢顶部纵桥向位移(m)。

在恒载及活载共同作用时下游侧左边跨、主跨及右边跨主缆的跨径

闭合的控制方程表达如下：

l_d,L＝L₁+Δx_d,L

l_d,R＝L₃-Δx_d,R

式中，L₁、L₂、L₃分别表示悬索桥左边跨、主跨和右边跨的跨径(m)；Δx_d,L为左桥塔下游侧塔肢顶部纵桥向位移(m)；Δx_d,R为右桥塔下游侧塔肢顶部纵桥向位移(m)。

(2)在恒载及活载共同作用时上游侧左边跨、主跨各节段及右边跨主缆无应力长度(S_cu,L,S_cu,i,S_cu,R)守恒的控制方程表达如下：

S_cu,L＝S’_cu,L

S_cu,i＝S’_cu,i,(1≤i≤n+1)

S_cu,R＝S’_cu,R

式中，S'_cu,i为主跨上游侧第i节段主缆在恒载作用下的无应力长度(m)；S'_cu,L为左边跨上游侧主缆在恒载作用下的无应力长度(m)；S'_cu,R为右边跨上游侧主缆在恒载作用下的无应力长度(m)。

在恒载及活载共同作用时下游侧左边跨、主跨各节段及右边跨主缆无应力长度(S_cd,L,S_cd,i,S_cd,R)守恒的控制方程表达如下：

S_cd,L＝S'_cd,L

S_cd,R＝S_c'_d,R

式中，S'_cd,i为主跨下游侧第i节段主缆在恒载作用下的无应力长度(m)；S'_cd,L为左边跨下游侧主缆在恒载作用下的无应力长度(m)；S'_cd,R为右边跨下游侧主缆在恒载作用下的无应力长度(m)。

在恒载及活载共同作用时上游侧各吊杆无应力长度(S_su,i)守恒的控制方程表达如下：

S_su,i＝S’_su,i,(1≤i≤n)

式中，S′_su,i为上游一侧第i根吊杆在恒载作用时的无应力长度(m)。

在恒载及活载共同作用时下游侧各吊杆无应力长度(S_sd,i)守恒的控制方程表达如下：

S_sd,i＝S’_sd,i,(1≤i≤n)

式中，S′_sd,i为下游一侧第i根吊杆在恒载作用时的无应力长度(m)。

(3)主梁端部扭转角为0控制方程表达如下：

θ(L₂)＝0

式中，L₂为悬索桥主跨跨径(m)。

(4)主梁y向受力(∑Y_i)平衡控制方程表达如下：

ΔF_u,R＝F_u,R-F’_u,R

ΔF_d,R＝F_d,R-F’_d,R

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；F_u,R、F_d,R分别表示上游、下游右支座在恒载和活载共同作用时的支座反力(kN)；F'_u,R、F'_d,R分别表示上游、下游右支座在恒载作用时的支座反力(kN)。

主梁x向受力(∑X_i)平衡控制方程表达如下:

式中ΔF_xu,i、ΔF_xd,i分别表示上游侧、下游侧第i根吊杆力的水平分力的改变量(kN)；F_u,_i、F_d,i分别表示上游、下游右支座在恒载和活载共同作用时的吊杆力(kN)；

为上游、下游侧吊杆力与竖直平面的夹角(rad)。

主梁绕x轴力矩(∑M_x)平衡控制方程表达如下：

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；b为吊杆横桥向间距的一半(m)；b'为支座横桥向间距的一半(m)。

主梁绕z轴力矩(∑M_z)平衡控制方程表达如下：

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；x_i表示恒载条件下，第i根吊杆下吊点至原点的纵桥向距离(m)；L₂为悬索桥主跨跨径(m)。

第六步：利用非线性规划求解控制方程组，解得同时满足所有控制方程组的基本未知量的值。

将第二步至第四步中所得到的用基本未知量表示的各项参数代入第五步中建立的控制方程组中，通过移项使得等式右边为0，等式左边视为4n+23个函数值W_i。利用非线性GRG方法对下式进行规划求解，即可得出悬索桥在偏心活载作用下的响应。

第七步：将基本未知量的值代入上述各响应的表达式，确定桥梁各构件的内力及变形

将步骤(6)中所获得的基本未知量的值代入(2.1)及(2.4)中的各段主缆高度差表达式、(2.2)中吊杆上吊点坐标表达式和(2.3)中主缆水平分力的递推关系式中，可确定在恒载和活载作用下主缆的内力和变形响应；

将步骤(6)中所获得的基本未知量的值代入(3.1)中的挠度表达式、(3.2)中相对转角的表达式和(3.3)中吊杆下吊点坐标表达式，可确定在恒载和活载作用下主梁的内力和变形响应；

将步骤(6)中所获得的基本未知量的值代入(4.1)中的吊杆倾角表达式可确定在恒载和活载作用下吊杆的倾斜角度，代入(4.3)中各塔肢顶部桥塔偏移量的表达式可确定桥塔在恒载和活载作用下的变形响应。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (9)

1.一种偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：包括如下步骤：

(1)确定可表达在偏心活载作用下全桥主要构件的内力及变形情况所需要的基本未知量；

该基本未知量包括：各跨主缆的悬链线方程参数、各节段悬链线的水平投影长度、主缆各节段水平分力、吊杆力、支反力、主梁整体水平向刚***移；

(2)分析各跨的主缆响应，用基本未知量表达：悬索桥上下游各节段主缆两端的高度差、无应力长度及主缆上的吊杆上吊点的坐标；

(3)分析主梁响应，用基本未知量表达：偏心活载作用下主梁各截面的挠度、横向扭转角及主梁上的吊杆下吊点的坐标；

(4)分析吊杆及桥塔响应，用基本未知量表达出上、下游吊杆倾角及无应力长度、桥塔各肢的偏移量；

(5)利用悬索桥上下游各节段主缆的两端的高度差与跨径闭合、主缆和吊杆的无应力长度守恒、主梁端部扭转角为0、主梁所受力与力矩平衡条件建立控制方程组；

(6)利用非线性规划方法求解步骤(5)建立的控制方程组，解得同时满足所有控制方程组的基本未知量的值；

(7)将基本未知量的值代入分析各跨的主缆响应表达式、分析主梁响应表达式、分析吊杆及桥塔响应表达式中，确定桥梁各构件的内力及变形。

2.根据权利要求1所述的偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：

所述步骤(1)中的基本未知量表达式如下：

左边跨上下游主缆的悬链线方程参数a_u,L、a_d,L；主跨上下游左侧第一节段主缆的悬链线方程参数a_u,1、a_d,1；右边跨上下游主缆的悬链线方程参数a_u,R、a_d,R；

左边跨上下游主缆水平投影长度l_u,L、l_d,L；主跨上下游n+1个主缆节段水平投影长度l_u,1～l_u,n+1、l_d,1～l_d,n+1；右边跨上下游主缆水平投影长度l_u,R、l_d,R；

左边跨上下游主缆的水平分力H_u,L、H_d,L；主跨上下游左侧第一节段主缆的水平分力H_u,1、H_d,1；右边跨上下游主缆的水平分力H_u,R、H_d,R；

主跨上下游n根吊杆的吊杆力F_u,1～F_u,n、F_d,1～F_d,n；

主梁左端上下游支反力F_u,L、F_u,R；主梁右端上下游支反力F_d,L、F_d,R；

主梁整体水平向刚***移v；

一共有4n+23个基本未知量,n为主跨一侧索面布置的吊杆数，下标u和d分别表示桥梁上游一侧和下游一侧的参数。

3.根据权利要求1所述的偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：

所述步骤(2)中表达上下游各节段主缆两端高度差、无应力长度以及主缆上的吊杆上吊点的坐标如下：

(2.1)每一段悬链线以其左端点为原点，建立局部坐标系，x轴水平向右，y轴竖直向下，在该坐标系下，各参数表述如下：

主跨上游一侧第i段主缆可视为悬链线，其左右端点的高度差Δy_u,i可以表示为：

式中，c_u,i＝-H_u,i/q，H_u,i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,i为主跨上游侧第i段悬链线方程参数；l_u,i为主跨上游侧第i段悬链线水平投影长度(m)；cosh为数学符号；

主跨上游一侧第i段主缆无应力长度S_u,i表示为：

式中，c_u,i＝-H_u,i/q，H_u,i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh为数学符号；a_u,i为主跨上游侧第i段悬链线方程参数；

主跨下游一侧第i段主缆左右端点的高度差Δy_d,i可以表示为：

式中，c_d,i＝-H_d,i/q，H_d,i为恒载及活载作用主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,i为主跨下游侧第i段悬链线方程参数；l_d,i为主跨下游侧第i段悬链线水平投影长度(m)；cosh为数学符号；

主跨下游一侧第i段主缆无应力长度S_d,i可以表示为：

式中，c_d,i＝-H_d,i/q，H_d,i为恒载及活载作用主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_d,i为主跨下游侧第i段悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；a_d,i为主跨下游侧第i段悬链线方程参数；sinh为数学符号；

(2.2)在恒载和偏心活载共同作用下，上游侧第i根吊杆的上吊点坐标(x_cu,i,y_cu,i,z_cu,i)表达为：

z_cu,i＝b

式中，Δx_u,L表示左端桥塔上游侧的塔肢顶部的偏移量(m)；b为各吊杆至主梁中心的横向距离(m)；H_B表示主跨主缆左侧IP点标高(m)；H_E表示主梁左侧截面中心的标高(m)；l_u,i为主跨上游侧第i段主缆的水平投影长度(m)；Δy_u,j表示主跨上游侧第j段主缆悬链线两端高度差(m)；

下游侧第i根吊杆的上吊点坐标(x_cd,i,y_cd,i,z_cd,i)表达为：

z_cd,i＝-b

式中，Δx_d,L表示左端桥塔下游侧的塔肢顶部的偏移量(m)；b为各吊杆至主梁中心的横向距离(m)；H_B表示主跨主缆左侧IP点标高(m)；H_E表示主梁左侧截面中心的标高(m)；l_d,i为主跨下游侧第i段主缆的水平投影长度(m)；Δy_d,j表示主跨下游侧第j段主缆悬链线两端高度差(m)；

(2.3)恒载和活载共同作用下，主跨上游侧的相邻节段主缆的水平分力H和悬链线方程参数a之间的递推关系式为:

式中，F_u,i为上游侧第i根吊杆的吊杆力(kN)；

为上游侧第i根吊杆力与竖直平面的夹角(rad)；a_u,i为主跨上游侧第i段悬链线方程参数；c_u,i＝-H_u,,i/q，H_u,i为恒载及活载作用下主跨上游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_u,i为主跨上游侧第i段主缆的水平投影长度(m)；arcsinh为数学符号；

主跨下游侧的相邻节段主缆的水平分力H和悬链线方程参数a之间的递推关系式为：

式中，F_d,i为下游侧第i根吊杆的吊杆力(kN)；

为下游侧第i根吊杆力与竖直平面的夹角(rad)；a_d,i为主跨下游侧第i段悬链线方程参数；c_d,i＝-H_d,,i/q，H_d,i为恒载及活载作用下主跨下游侧第i段主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；l_d,i为主跨下游侧第i段主缆的水平投影长度(m)；arcsinh为数学符号；

(2.4)上游侧左边跨的主缆段，以左锚固点为原点，纵桥向为x轴，向右为正方向，竖直向下为y轴，建立局部坐标系；

上游左边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_u,L和无应力长度S_u,L可表示为：

式中，c_u,L＝-H_u,L/q，H_u,L为恒载及活载作用左边跨上游主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,L为上游左边跨主缆悬链线方程参数；l_u,L为上游左边跨主缆水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

右边跨段的主缆以右锚固点为原点，纵桥向为x轴，向右为正，竖直向下为y轴建立局部坐标系，则上游右边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_u,R和无应力长度S_u,R可表示为:

式中，c_u,R＝-H_u,R/q，H_u,R为恒载及活载作用上游右边跨主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_u,R为上游右边跨悬链线方程参数；l_u,R为上游右边跨悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号；

以下游侧左边跨的主缆段的左锚固点为原点，纵桥向为x轴，向右为正方向，竖直向下为y轴，建立局部坐标系；下游左边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_d,L和无应力长度S_d,L可表示为：

式中，c_d,L＝-H_d,L/q，H_d,L为恒载及活载作用左边跨下游主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,L为下游左边跨主缆悬链线方程参数；l_d,L为下游左边跨主缆水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号；

右边跨段的主缆以右锚固点为原点，纵桥向为x轴，向右为正，竖直向下为y轴建立局部坐标系，则下游右边跨主缆的端点，即锚固点与散索鞍IP点在偏心活载作用下的高度差Δy_d,R和无应力长度S_d,R可表示为:

式中，c_d,R＝-H_d,R/q，H_d,R为恒载及活载作用下游右边跨主缆的水平分力(kN)；q为主缆单位长度自重(kN/m)；a_d,R为下游右边跨悬链线方程参数；l_d,R为下游右边跨悬链线水平投影长度(m)；E_cA_c表示主缆轴向抗拉刚度(kN)；sinh、cosh为数学符号。

4.根据权利要求1所述偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：所述步骤(3)中分析主梁响应，用基本未知量表达主梁的挠度和横向扭转角以及主梁上的吊杆下吊点的坐标如下：

(3.1)以主梁左端截面形心为原点，纵桥向为x轴，竖直向下为y轴，横向为z轴建立局部坐标系；主梁受到两组吊杆力F_u,i、F_d,i(1≤i≤n)以及四个支座反力F_u,L、F_u,R、F_d,L、F_d,R的作用。

在偏心活载作用下，主梁可被n根吊杆划分为n+1个区段，每一段主梁的挠度w(x)为：

ΔF_u,L＝F_u,L-F′_u,L

ΔF_d,L＝F_d,L-F′_d,L

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力在受到活载作用前后的变化量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；F'_u,i、F'_d,i分别表示上游、下游第i根吊杆在恒载作用下的吊杆力(kN)；F'_u,L、F'_d,L分别表示上游、下游左支座在恒载作用下的支座反力(kN)；E_bI_b表示主梁抗弯刚度(kN·m²)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)；p为偏心活载集度(kN/m)；x表示所求截面与该区段主梁最左端的距离(m)；

为下游侧第i根吊杆力与竖直平面的夹角(rad)；t为积分所取微元至主梁左端的距离(m)；dt为积分变量；j为所求截面左侧距离最近的吊杆编号；

(3.2)距离主梁左端距离为x的截面Ⅰ对于左侧距该截面最近的吊杆右截面Ⅱ的相对转角为：

式中，GI_P为主梁的自由扭转刚度(kN·m/rad)；ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力在受到活载作用前后的变化量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；p为偏心活载集度(kN/m)；b为吊杆横桥向间距的一半(m)；b'为支座横桥向间距的一半(m)；e为偏心活载作用位置与桥梁纵轴线的横向距离(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)；x表示所求截面与该区段主梁最左端的距离(m)；t为积分所取微元至主梁左端的距离(m)；t为积分所取微元至主梁左端的距离(m)；dt为积分变量；j为所求截面左侧距离最近的吊杆编号；x_j表示恒载条件下第j根吊杆下吊点至原点的纵桥向距离(m)；

当x＝x_m+1时，可求得相邻两根吊杆截面的相对扭转角，根据θ(0)＝0，通过叠加即可知主梁各截面扭转角的值；

(3.3)主梁上游侧第i根吊杆的下吊点坐标(x_bu,i,y_bu,i,z_bu,i)可以表示为：

x_bu,i＝x_i+v

y_bu,i＝-w(x_i)-bθ(x_i)

z_bu,i＝b

式中，v表示主梁整体水平方向的刚***移(m)；w表示偏心活载作用下主梁发生的竖向挠度(m)；θ表示主梁发生的扭转角(rad)；b为吊杆横桥向间距的一半(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)；

主梁下游侧第i根吊杆的下吊点坐标(x_bd,i,y_bd,i,z_bd,i)可以表示为

x_bd,i＝x_i+v

y_bd,i＝-w(x_i)+bθ(x_i)

z_bd,i＝-b

式中，v表示主梁整体水平方向的刚***移(m)；w表示偏心活载作用下主梁发生的竖向挠度(m)；θ表示主梁发生的扭转角(rad)；b为吊杆横桥向间距的一半(m)；x_i表示恒载条件下第i根吊杆下吊点至原点的纵桥向距离(m)。

5.根据权利要求1所述的偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：所述步骤(4)中表达上下游吊杆倾角及无应力长度、桥塔各肢的偏移量的具体步骤如下：

(4.1)根据步骤(2.2)中建立的上游侧第i根吊杆的上吊点坐标(x_cu,i,y_cu,i,z_cu,i)与步骤(3.3)中建立的上游侧第i根吊杆的下吊点坐标(x_bu,i,y_bu,i,z_bu,i)计算，在偏心活载作用下，上游侧第i根吊杆的倾角

可以表示为：

根据步骤(2.2)中建立的下游侧第i根吊杆的上吊点坐标(x_cd,i,y_cd,i,z_cd,i)与步骤(3.3)中建立的下游侧第i根吊杆的下吊点坐标(x_bd,i,y_bd,i,z_bd,i)计算，下游侧第i根吊杆的倾角

可以表示为：

(4.2)上游侧第i根吊杆的无应力长度S_su,i可以表示为：

式中，E_sA_s为吊杆抗拉刚度(kN)；S_s’_u,i为上游一侧第i根吊杆在恒载作用下的无应力长度(m)；F_u,i为上游侧第i根吊杆的吊杆力；x_cu,i、x_bu,i分别为上游第i根吊杆上吊点、下吊点的x坐标；y_cu,i、y_bu,i分别为上游第i根吊杆上吊点、下吊点的y坐标。

下游侧第i根吊杆的无应力长度S_sd,i可以表示为：

式中，E_sA_s为吊杆抗拉刚度(kN)；S’_sd,i为下游一侧第i根吊杆在恒载作用下的无应力长度(m)；F_d,i为下游侧第i根吊杆的吊杆力；x_cd,i、x_bd,i分别为下游第i根吊杆上吊点、下吊点的x坐标；y_cd,i、y_bd,i分别为下游第i根吊杆上吊点、下吊点的y坐标；

(4.3)左桥塔上游、下游塔肢顶部的偏移量Δx_u,L、Δx_d,L分别为：

式中，k为单肢桥塔的刚度系数(kN/m)；k_ξ为整个桥塔的抗扭刚度系数(kN/rad)；d为每肢桥塔距离桥轴线距离(m)；H_u,L、H_d,L分别为上下游左边跨主缆水平分力(kN)；H_u,1、H_d,1分别为上下游主跨左侧第一段主缆悬链线水平分力(kN)；

右桥塔上游、下游塔肢顶部的偏移量Δx_u,R、Δx_d,R分别为：

式中，k为单肢桥塔的刚度系数(kN/m)；k_ξ为整个桥塔的抗扭刚度系数(kN/rad)；d为每肢桥塔距离桥轴线距离(m)；H_u,R、H_d,R分别为上下游右边跨主缆水平分力(kN)；H_u,n+1、H_d,n+1分别为上下游主跨最右侧节段(第n+1段)主缆悬链线水平分力(kN)。

6.根据权利要求1所述的偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：所述步骤(5)中建立控制方程组的具体步骤如下：

(5.1)在恒载及活载共同作用时悬索桥上游侧左边跨、主跨及右边跨主缆的高度差

闭合的控制方程表达如下：

Δy_u,L＝H_A-H_B

Δy_u,R＝H_D-H_C

式中，H_A表示主缆左锚固点标高(m)；H_B表示主跨主缆左侧IP点标高(m)；H_C表示主跨主缆右侧IP点标高(m)；H_D表示主缆右锚固点标高(m)；

在恒载及活载共同作用时下游侧左边跨、主跨及右边跨主缆的高度差Δy_d,L,

Δy_d,R闭合的控制方程表达如下：

Δy_d,L＝H_A-H_B

Δy_d,R＝H_D-H_C

式中，H_A表示主缆左锚固点标高(m)；H_B表示主跨主缆左侧IP点标高(m)；H_C表示主跨主缆右侧IP点标高(m)；H_D表示主缆右锚固点标高(m)；

在恒载及活载共同作用时悬索桥上游侧左边跨、主跨及右边跨主缆的跨径

l_u,R闭合的控制方程表达如下：

l_u,L＝L₁+Δx_u,L

l_u,R＝L₃-Δx_u,R

式中，L₁、L₂、L₃分别表示悬索桥左边跨、主跨和右边跨的跨径(m)；Δx_u,L为左桥塔上游侧塔肢顶部纵桥向位移(m)；Δx_u,R为右桥塔上游侧塔肢顶部纵桥向位移(m)；

在恒载及活载共同作用时下游侧左边跨、主跨及右边跨主缆的跨径

闭合的控制方程表达如下：

l_d,L＝L₁+Δx_d,L

l_d,R＝L₃-Δx_d,R

式中，L₁、L₂、L₃分别表示悬索桥左边跨、主跨和右边跨的跨径(m)；Δx_d,L为左桥塔下游侧塔肢顶部纵桥向位移(m)；Δx_d,R为右桥塔下游侧塔肢顶部纵桥向位移(m)；

(5.2)在恒载及活载共同作用时上游侧左边跨、主跨各节段及右边跨主缆无应力长度S_cu,L,S_cu,i,S_cu,R守恒的控制方程表达如下：

S_cu,L＝S′_cu,L

S_cu,i＝S′_cu,i,(1≤i≤n+1)

S_cu,R＝S′_cu,R

式中，S'_cu,i为主跨上游侧第i节段主缆在恒载作用下的无应力长度(m)；S'_cu,L为左边跨上游侧主缆在恒载作用下的无应力长度(m)；S'_cu,R为右边跨上游侧主缆在恒载作用下的无应力长度(m)；

在恒载及活载共同作用时下游侧左边跨、主跨各节段及右边跨主缆无应力长度(S_cd,L,S_cd,i,S_cd,R)守恒的控制方程表达如下：

S_cd,L＝S′_cd,L

S_cd,i＝S′_cd,i,(1≤i≤n+1)

S_cd,R＝S′_cd,R

式中，S'_cd,i为主跨下游侧第i节段主缆在恒载作用下的无应力长度(m)；S'_cd,L为左边跨下游侧主缆在恒载作用下的无应力长度(m)；S'_cd,R为右边跨下游侧主缆在恒载作用下的无应力长度(m)；

在恒载及活载共同作用时上游侧各吊杆无应力长度(S_su,i)守恒的控制方程表达如下：

S_su,i＝S′_su,i,(1≤i≤n)

式中，S’_su,i为上游一侧第i根吊杆在恒载作用时的无应力长度(m)；

在恒载及活载共同作用时下游侧各吊杆无应力长度(S_sd,i)守恒的控制方程表达如下：

S_sd,i＝S′_sd,i,(1≤i≤n)

式中，S’_sd,i为下游一侧第i根吊杆在恒载作用时的无应力长度(m)；

(5.3)主梁端部扭转角为0控制方程表达如下：

θ(L₂)＝0

式中，L₂为悬索桥主跨跨径(m)；

(5.4)主梁y向受力(∑Y_i)平衡控制方程表达如下：

ΔF_u,R＝F_u,R-F′_u,R

ΔF_d,R＝F_d,R-F′_d,R

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,L、ΔF_d,L分别表示上下游左支座的支座反力在受到活载作用前后的变化量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；F_u,R、F_d,R分别表示上游、下游右支座在恒载和活载共同作用时的支座反力(kN)；F'_u,R、F'_d,R分别表示上游、下游右支座在恒载作用时的支座反力(kN)；

主梁x向受力(∑X_i)平衡控制方程表达如下：

式中ΔF_xu,i、ΔF_xd,i分别表示上游侧、下游侧第i根吊杆力的水平分力的改变量(kN)；F_u,i、F_d,i分别表示上游、下游右支座在恒载和活载共同作用时的吊杆力(kN)；

为上游、下游侧吊杆力与竖直平面的夹角(rad)；

主梁绕x轴力矩(∑M_x)平衡控制方程表达如下：

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；b为吊杆横桥向间距的一半(m)；b'为支座横桥向间距的一半(m)；

主梁绕z轴力矩(∑M_z)平衡控制方程表达如下：

式中，ΔF_yu,i、ΔF_yd,i分别表示上游侧、下游侧第i根吊杆力的竖向分力的改变量(kN)；ΔF_u,R、ΔF_d,R分别表示上下游右支座的支座反力在受到活载作用前后的变化量(kN)；x_i表示恒载条件下，第i根吊杆下吊点至原点的纵桥向距离(m)；L₂为悬索桥主跨跨径(m)。

7.根据权利要求3所述偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：所述的步骤(3.1)中在第j-1根吊杆与第j根吊杆间的第j段主梁挠度表达式中由于不定积分产生的常数项系数与一次项系数C_j,0、C_j,1以及第j+1段的主梁挠度表达式中由于不定积分产生的常数项系数与一次项系数C_j+1,0、C_j+1,1的递推关系式如下：

式中，ΔF_yu,j、ΔF_yd,j分别表示上游侧、下游侧第j根吊杆力的竖向分力的改变量(kN)；x_j表示恒载条件下，第j根吊杆下吊点至原点的纵桥向距离(m)；

由主梁左端挠度w_L＝0，求得第1根吊杆左侧的主梁挠度表达式，即第1段主梁挠度中由于不定积分产生的常数项系数C_1,0，进而递推出所有段挠度表达式中的常数项系数C_i,0，包括C_n,0；

主梁右端挠度w_R＝0,可以求得第n根吊杆右侧的主梁挠度表达式，即第n+1段吊杆挠度表达式由于不定积分产生的常数项系数与一次项系数C_n+1,0、C_n+1,1的关系式，进而得到C_n+1,1，则可以反推出所有段挠度表达式中的一次项系数C_i,1；

将所求得的常数项系数及一次项系数代入相对应的主梁区段的挠度表达式，即可推得主梁任一截面的挠度。

8.根据权利要求1所述偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：步骤(6)中利用非线性规划求解控制方程组，解得同时满足所有控制方程组的基本未知量的值的具体步骤如下：

利用步骤(5)中建立的4n+23个控制方程，将步骤(2)至步骤(4)中用基本未知量表示的各项参数代入控制方程组中，通过移项使得等式右边为0，等式左边视为4n+23个函数值W_i；利用非线性GRG方法对下式进行规划求解，即可得出悬索桥在偏心活载作用下的响应；

9.根据权利要求3或4或5或7所述偏心活载作用下的悬索桥全桥响应的计算方法，其特征在于：将步骤(6)中所获得的基本未知量的值代入(2.1)及(2.4)中的各段主缆高度差表达式、(2.2)中吊杆上吊点坐标表达式和(2.3)中主缆水平分力的递推关系式中，确定在恒载和活载作用下主缆的内力和变形响应；

将步骤(6)中所获得的基本未知量的值代入(3.1)中的挠度表达式、(3.2)中相对转角的表达式和(3.3)中吊杆下吊点坐标表达式，确定在恒载和活载作用下主梁的内力和变形响应；

将步骤(6)中所获得的基本未知量的值代入(4.1)中的吊杆倾角表达式可确定在恒载和活载作用下吊杆的倾斜角度，代入(4.3)中各塔肢顶部桥塔偏移量的表达式可确定桥塔在恒载和活载作用下的变形响应。

Images