CN111259582B - 一种悬索桥索鞍处主缆长度快速精确计算方法 - Google Patents

Abstract

一种悬索桥索鞍处主缆长度快速精确计算方法，涉及一种主缆长度计算方法。首先，根据主缆与索鞍的几何关系推导了索鞍处主缆曲线修正算法，然后，利用牛顿‑拉菲森迭代法对所得二元非线性方程组进行求解，最后，选取常见的主索鞍与散索鞍两组算例验证了该方法的可靠性。相比于传统算法，表达形式更加明确，对参数初始值设置没有严格要求，均可达到快速收敛的效果，增强了其可操作性，大大提高了计算效率。可方便地应用于建设期间主缆曲线长度以及索鞍位置的确定，为大跨径悬索桥的施工控制提供了有力的工具，使其施工完成后更加接近设计状态。

Description

一种悬索桥索鞍处主缆长度快速精确计算方法

技术领域

本发明涉及一种主缆长度计算方法，尤其是一种悬索桥索鞍处主缆长度快速精确计算方法，属于悬索桥索鞍曲线修正计算领域。

背景技术

悬索桥由于具有受力性能合理、抗震性能较好、大跨轻盈的特点，是特大跨径桥梁的首选桥型。当桥梁跨径超过1200m时，悬索桥被认为是最具有竞争力的桥型方案，随着跨径的增大，悬索桥的优势就越明显。

主缆是悬索桥至关重要的受力构件，悬索桥施工质量的指标很大程度上要考虑桥梁施工完成后的实际主缆线形与设计主缆线形是否一致。对于特大跨径悬索桥，主缆的安全系数要求不小于2.5，这必然要对悬索桥进行精确的求解。目前的悬索桥计算分析方法已经比较成熟，但由于计算模型在假设时精度达不到要求，所以在处理一些局部受力问题时，悬索桥的计算内力和线形与实际受力与变形有一定偏差。

悬索桥的索鞍起到了转向或分散主缆的作用，并限制了主缆的变形。对于悬索桥的成桥阶段，主缆在进出鞍座处与鞍座是处于相切状态的，两者在切点间是紧密接触的，主缆在鞍座内不能有任何的相对滑移。而在施工阶段，考虑到施工的实际情况，主缆与鞍座接触的切点是不断变化位置的，主缆与鞍座间存在着接触非线性。在进行悬索桥计算时，如果不能准确考虑主缆与鞍座间复杂的接触关系，主缆和吊索的下料长度将会减小，主缆架设控制标高的精度也会相应降低。

基于有限单元法的悬索桥的施工分析计算中，大多数方法是通过利用成桥理论点来建模分析，但这不能确保主缆与鞍座处于相切状态，从而造成了主缆与鞍座出现相交或分离。为了充分考虑鞍座对主缆的约束作用，潘永仁采用杆单元来模拟，桥塔顶部采用两个杆单元模拟，通过主从约束的方式与鞍座杆单元连接，通过修改杆单元的节点坐标，来模拟鞍座的顶推作用，改变鞍座杆单元的无应力长度，从而使得主缆在脱离点的斜率等于鞍槽曲线上此点处的斜率，进而得到切点的位置。徐君兰则采用四个梁单元来模拟鞍座及其顶推作用，用一个直杆代替散索鞍，该方法与潘永仁的方法实质是相同的。上述方法在修正悬空段主缆的无应力长度时未考虑切点相对于鞍座位置的变化，使得悬索桥求解精度较低。另外，跨内主缆线形由假定的切点坐标迭代求出后，还必须验证主缆与鞍座是否相切，否则需重新迭代。因此，该方法计算效率也不高。

有形状的索鞍在其表面上包裹并支撑主缆，在施工过程中，主缆与索鞍间的切点位置在一直变化着。在悬索桥计算分析中，由理论顶点假设推导出的主缆线形而得到的考虑索鞍曲线修正后实际的主缆形状是比较繁琐的，其计算关键在于找到主缆与索鞍间的相切点。目前常用的解析算法是求解八元非线性方程组，即借助悬链线公式得到4个几何关系式、利用力平衡条件得到2个平衡方程、依据切点和索鞍圆心的距离关系得到2个公式，再加上一些工程经验性的约束条件后通过数值分析或Mathcad等软件进行计算，多位学者对此均有论述，但该方法计算效率并不高。另外，有学者采用通用有限元软件建立主缆与索鞍的接触关系进行数值仿真分析，也有部分学者考虑了主缆与鞍座之间的摩擦来提高计算精度。上述方法在精度方面虽然能满足工程需求，但计算比解析法更困难，效率更低。

目前的索鞍曲线修正计算分析理论推导繁琐、难以操作、效率较低，所以在实际工程中，索鞍曲线修正的影响基本只在悬索桥的成桥状态下考虑。但是为了满足人们建造更大跨径悬索桥的需求，仅在成桥状态进行索鞍修正无法满足施工对线形的控制要求。因此，推导一种便于计算、效率更高的索鞍曲线修正算法是十分有必要的，不仅可以实时进行索鞍曲线修正，还可为大跨径悬索桥的施工控制提供一定的参考。

传统的索鞍曲线修正计算方法理论推导复杂、不易理解、操作性不强，导致实际工程项目中基本只在成桥状态下考虑了索鞍曲线修正的影响。为此，本发明从几何关系角度对索鞍曲线修正算法进行推导，并借助牛顿-拉斐森迭代法，提供了一种可行的数值求解方法。

发明内容

针对背景技术中存在的不足，本发明提供了一种悬索桥索鞍处主缆长度快速精确计算方法。

为实现上述目的，本发明采取下述技术方案：一种悬索桥索鞍处主缆长度快速精确计算方法，包括以下步骤：

步骤一：主缆索段力学分析，包括：

将一段只受自重的缆索分割出来，H和V分别表示主缆的水平分力和竖向分力，q为缆索沿弧长均布的自重荷载，l为主缆在x轴的投影长度，c为主缆在y轴的投影长度，得到主缆在自重作用下的曲线形状为悬链线，具体表达式为：

对于原点处H、V已知，则有：

索鞍处的主缆线形x为水平向，y为竖向，假设一个虚拟的主缆交点，即理论顶点，将局部坐标系原点建立在理论顶点上，其坐标为(0,0)，左侧主缆与索鞍的切点标为切点1，其局部坐标为(x₁,y₁)，右侧主缆与索鞍的切点标为切点2，其局部坐标为(x₂,y₂)，索鞍与主缆接触面为圆弧，索鞍圆心局部坐标为(x₃,y₃)，半径为R，x₁、x₂与y₁、y₂的关系满足公式(1)，用y＝f(x)来表示，

根据左侧主缆与索鞍相切可以得到，切点1与索鞍圆心的连线斜率k₁表达式为

根据右侧主缆与索鞍相切可以得到，切点2与索鞍圆心的连线斜率k₂表达式为

切点1与切点2同在以索鞍圆心为圆心，半径为R的圆上，根据几何关系可得：

左侧主缆方程为

右侧主缆方程为

另外，可以得到

将f₁′(x₁)、f₂′(x₂)、y₁、y₂代入公式(6)后得到求解主缆与索鞍切点的二元非线性方程组

步骤二：数值计算迭代算法，

通过牛顿-拉斐森迭代法可对上述二元非线性方程组(11)进行数值求解，则有

X_i+1＝X_i-A_i ^-1G_i (14)，

取切点1、切点2的初始x坐标代入牛顿-拉斐森迭代，得到x₁、x₂的数值解，y₁、y₂代入前文提到的悬链线方程即可得到；

步骤三：根据前述理论推导，索鞍曲线修正算法简化为二元非线性方程组，将其嵌入有限元分析软件进行求解。

与现有技术相比，本发明的有益效果是：

1.根据主缆与索鞍的几何关系，推导了索鞍处主缆曲线的修正算法，相比于传统算法，减少了6个方程与6个初始输入参数，表达形式更明确；

2.利用牛顿-拉菲森迭代法，对所得二元非线性方程组进行求解，关于输入参数的初始值设置没有严格要求，具有很强的可操作性；

3.针对常见的主索鞍与散索鞍结构，相比于传统算法，迭代次数减少50％，计算时间不足10％，大大提高了计算效率，且精度可满足要求。

附图说明

图1是悬链线缆索示意图；

图2是索鞍与主缆关系示意图。

具体实施方式

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

本发明公开了一种悬索桥索鞍处主缆长度快速精确计算方法，具体为

主缆索段力学分析：

现代悬索桥主缆一般由钢丝集束形成，且其线形平顺，在转折处曲率半径很大，因此相比其轴向抗拉刚度，抗弯刚度很小，具有很强的柔性索结构特征，本发明在推导悬索桥主缆的算法时，对主缆做三项前提假设：

1、主缆是理想柔性的索，只承受轴向拉力，不承受轴向压力和弯矩；

2、在设计悬索桥时通常会给主缆拉应力设置较大的安全系数，主缆在施工和使用过程中，材料均不会达到塑性阶段，因此其符合胡克定律，是线弹性的；

3、主缆横截面积在外荷载作用下变化量微小，不予考虑。

将一段只受自重的缆索分割出来，参照图1所示的悬链线缆索示意图，H和V分别表示主缆的水平分力和竖向分力，q为缆索沿弧长均布的自重荷载，l为主缆在x轴的投影长度，c为主缆在y轴的投影长度。

显然，可以得到主缆在自重作用下的曲线形状为悬链线，具体表达式为

其中，各参数表达式为

对于原点处H、V已知，则有：

索鞍处的主缆线形参照图2所示的索鞍与主缆关系示意图，x为水平向，y为竖向，整体分析时为方便计算，假设一个虚拟的主缆交点，即图中的理论顶点，本发明将局部坐标系原点建立在理论顶点上，其坐标为(0,0)，左侧主缆与索鞍的切点标为切点1，其局部坐标为(x₁,y₁)，右侧主缆与索鞍的切点标为切点2，其局部坐标为(x₂,y₂)，索鞍与主缆接触面为圆弧，索鞍圆心局部坐标为(x₃,y₃)，半径为R。另外，x₁、x₂与y₁、y₂的关系满足悬链线公式(1)，为推导过程简洁，用y＝f(x)来表示。

根据左侧主缆与索鞍相切可以得到，切点1与索鞍圆心的连线斜率k₁表达式为

根据右侧主缆与索鞍相切可以得到，切点2与索鞍圆心的连线斜率k₂表达式为

其中，k为理论切点位置索(索鞍)的切线斜率。

切点1与切点2同在以索鞍圆心为圆心，半径为R的圆上，根据几何关系可得：

左侧主缆方程为

右侧主缆方程为

其中，公式(7)与公式(8)中的参数如下

另外，可以得到

将f₁′(x₁)、f₂′(x₂)、y₁、y₂代入公式(6)后得到求解主缆与索鞍切点的二元非线性方程组

数值计算迭代算法：

通过牛顿-拉斐森迭代法可对上述二元非线性方程组(11)进行数值求解，令

再令

其中，A_i为雅克比矩阵，矩阵各元素表达式如下：

根据牛顿-拉斐森迭代，则有

X_i+1＝X_i-A_i ^-1G_i (14)

取切点1、切点2的初始x坐标代入牛顿-拉斐森迭代，得到x₁、x₂的数值解，y₁、y₂代入前文提到的悬链线方程即可得到。值得注意的是，传统的八元非线性方程组的初始值选取原则较为繁琐，而本发明算法对初始值的选取基本没有要求。

根据前述理论推导，索鞍曲线修正算法由传统的八元非线性方程组(传统方法)简化为二元非线性方程组(本发明方法)，可较为方便地将其嵌入有二次开发功能的有限元分析软件如TDV、Ansys、Abaqus等，也可通过Matlab等编程语言方便地进行求解。

实施例1：通过Matlab选取主索鞍为算例，对其正确性与计算效率进行验证

悬索桥的其中一个主索鞍理论顶点坐标为(230.000m,131.425m)，主缆面积A₁＝A₂＝0.408973m²，主缆自重集度为q₁＝q₂＝33kN/m，弹性模量E＝198000MPa。索鞍半径R＝6m，索鞍左右两侧的主缆索力水平分量H₁＝H₂＝189500kN，V₁＝90622.7kN，V₂＝73504.1kN。

通过编程，对传统方法与本发明方法(上述理论推导)进行对比，初始取值见表1，计算过程见表2，计算结果见表3。

表1主索鞍曲线修正初始取值

表2主索鞍曲线修正计算过程

表3主索鞍位置及主缆曲线计算结果

注：表中左切点坐标(x₁,y₁)、右切点坐标(x₂,y₂)、索鞍圆心坐标(x₃,y₃)均为整体坐标系坐标。

实施例2：通过Matlab选取散索鞍为算例，对其正确性与计算效率进行验证

悬索桥的其中一个散索鞍理论顶点坐标为(0m,54m)，主缆面积A₁＝A₂＝0.408973m²，主缆自重q₁＝q₂＝33kN/m，弹性模量E＝198000MPa。索鞍半径R＝6m，索鞍左右两侧的主缆索力水平分量H₁＝H₂＝189500kN，V₁＝137557kN，V₂＝-41804.3kN。

通过编程，对传统方法与本文方法(上述理论推导)进行对比，初始取值见表4，计算过程见表5，结果见表6。

表4散索鞍曲线修正初始取值

表5散索鞍曲线修正计算过程

表6散索鞍位置及主缆曲线计算结果

从上述两个实施例的计算结果和计算效率对比可知，本发明提出的二元非线性方程组求解方法与目前普遍采用的八元非线性方程组求解方法计算结果在0.01mm的数量级上可以保持一致性，对于工程应用可以认为两种方法的计算结果是没有差别的。计算效率方面，本发明提出的方法迭代次数减少了50％，计算耗时仅为目前普遍采用方法的10％。更重要的是，在设置初始参数时，本文提出的方法要简便很多，只需要设置x₁和x₂两个参数即可，且没有特殊要求。一般情况下，为便于编写程序可以直接设置初始x₁＝x₃-R/2，x₂＝x₃+R/2，即可很快迭代收敛。然而，目前普遍采用方法，参考文献(邵旭东,程翔云,李立峰.桥梁设计与计算[M].北京:人民交通出版社,2007)中的推导过程表明，需要设置8个初始参数，这8个参数中只有x₁、x₂较易设置，其余6个参数均需有一定的经验或者经过数次试算方可获得较好的计算结果，这给该算法实现程序自动进行索鞍曲线修正计算带来很大困难。由此说明，本发明所提出方法更为方便，大大提高了计算效率，且精度可满足要求。

对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的装体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同条件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。

此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。

Claims (1)

1.一种悬索桥索鞍处主缆长度快速精确计算方法，其特征在于：所述计算方法包括以下步骤：

步骤一：主缆索段力学分析，包括：

将一段只受自重的缆索分割出来，H和V分别表示主缆的水平分力和竖向分力，q为缆索沿弧长均布的自重荷载，l为主缆在x轴的投影长度，c为主缆在y轴的投影长度，得到主缆在自重作用下的曲线形状为悬链线，具体表达式为：

对于原点处H、V已知，则有：

索鞍处的主缆线形x为水平向，y为竖向，假设一个虚拟的主缆交点，即理论顶点，将局部坐标系原点建立在理论顶点上，其坐标为(0,0)，左侧主缆与索鞍的切点标为切点1，其局部坐标为(x₁,y₁)，右侧主缆与索鞍的切点标为切点2，其局部坐标为(x₂,y₂)，索鞍与主缆接触面为圆弧，索鞍圆心局部坐标为(x₃,y₃)，半径为R，x₁、x₂与y₁、y₂的关系满足公式(1)，用y＝f(x)来表示，

根据左侧主缆与索鞍相切可以得到，切点1与索鞍圆心的连线斜率k₁表达式为

根据右侧主缆与索鞍相切可以得到，切点2与索鞍圆心的连线斜率k₂表达式为

切点1与切点2同在以索鞍圆心为圆心，半径为R的圆上，根据几何关系可得：

左侧主缆方程为

右侧主缆方程为

另外，可以得到

将f₁′(x₁)、f′₂′(x₂)、y₁、y₂代入公式(6)后得到求解主缆与索鞍切点的二元非线性方程组

步骤二：数值计算迭代算法，

通过牛顿-拉斐森迭代法可对上述二元非线性方程组(11)进行数值求解，则有

X_i+1＝X_i-A_i ^-1G_i (14)，

取切点1、切点2的初始x坐标代入牛顿-拉斐森迭代，得到x₁、x₂的数值解，y₁、y₂代入前文提到的悬链线方程即可得到；

步骤三：根据前述理论推导，索鞍曲线修正算法简化为二元非线性方程组，将其嵌入有限元分析软件进行求解。

Images