CN112035928B - 大跨径悬索管道桥成桥主索线形计算方法 - Google Patents
大跨径悬索管道桥成桥主索线形计算方法 Download PDFInfo
- Publication number
- CN112035928B CN112035928B CN202010899711.2A CN202010899711A CN112035928B CN 112035928 B CN112035928 B CN 112035928B CN 202010899711 A CN202010899711 A CN 202010899711A CN 112035928 B CN112035928 B CN 112035928B
- Authority
- CN
- China
- Prior art keywords
- cable
- main
- main cable
- bridge
- sling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 42
- 239000000725 suspension Substances 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 54
- 238000013461 design Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 238000012937 correction Methods 0.000 claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000005483 Hooke's law Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bridges Or Land Bridges (AREA)
Abstract
本发明涉及一种大跨径悬索管道桥成桥主索线形计算方法,属于油气管道跨越设计领域,包括根据空间抛物线理论估算主索始端三向分力初值;利用节线法迭代计算并得到修正后的主索始端三向分力;以主索始端三向分力和设计控制点为迭代初值,利用空间分段悬链线理论计算主索线形及各索段无应力长度;利用修正影响矩阵法获得迭代变量增量,反复迭代,修正主索始端三向分力,重新计算得到主索线形及各索段无应力长度。本发明解决了以往设计不精、假设过多,计算误差太大的问题,采用本方法开展设计和计算,可以大大提高悬索管道桥成桥主索线形及无应力长度计算精度。
Description
技术领域
本发明涉及一种大跨径悬索管道桥成桥主索线形计算方法,属于油气管道跨越设计领域。
背景技术
管道悬索桥一般由主索、吊索、风索、风拉索、稳定索(共轭索)、加劲梁、索塔和锚碇等组成。由于其主要用途是输送石油、天然气或水,桥面上通常会设置支座用于架设管道,并设置简易检修道以方便管道的施工和后期维修。管道悬索桥宽跨比小,结构刚度小,抗风问题敏感,一般需设置风索体系。管道悬索桥的受力特征是:作用在桥面上的荷载和因风索产生的竖向荷载通过吊索传至主索,再传至索塔和锚碇;作用在桥面的水平荷载通过风拉索传至风索主索,再传至锚碇;传力途径明确。
管道悬索桥是以柔性索结构为主要承重结构的缆索体系。索结构受力后,反映出小应变大变形的受力特征。在荷载作用下,荷载与变形呈现出明显的非线性的关系,经典结构力学已不再适用,因为经典结构力学忽略了结构的微小变形,平衡方程建立在变形前的几何位置上。计算大变形的结构,应在变形后的位置上建立平衡方程,并且需要迭代计算。
目前,工程上设计人员主要利用平面抛物线理论对主索始端三向力估计作为迭代初值,对于空间主缆,该方法具有较大误差,可能会导致后续迭代不收敛。
发明内容
为解决上述问题,准确估计悬索跨越成桥主索始端三向力,本发明的目的在于提供一种大跨径悬索管道桥成桥主索线形计算方法。通过引入节线法对主索始端三向力进行进一步的精确估计,使得程序的收敛性大大提高。由于本方法中采用分段悬链线,用数值迭代计算空间主索线形,是一种精确的计算方法。分段悬链线方法不存在利用有限元法计算时的那种假设误差,它根据力学平衡条件和变形相容条件确定各部分索力和曲线形状,自动计入了索曲线的所有非线性,与有限元相比,计算精度大为提高。
本发明的目的是通过以下技术方案实现的:
一种大跨径悬索管道桥成桥主索线形计算方法,包括以下步骤:
步骤1,根据空间抛物线理论估算主索始端三向分力初值HFX0、HFY0和VF0;
步骤2,以主索始端三向分力初值HFX0、HFY0、VF0为迭代初值,利用节线法迭代计算并得到修正后的主索始端三向分力HFXL、HFYL、VFL;
步骤3,以主索始端三向分力HFXL、HFYL、VFL为迭代初值,由始端至末端对主索各索段进行计算,包括:
步骤3.1,根据第i索段起始端三向分力、第i索段纵桥向长度xi和主索平衡方程计算出主索第i索段无应力长度SFi及第i索段横、竖桥向长度yi、zi;
步骤3.2,对于第i根吊索,由步骤3.1中所述主索平衡方程计算得到吊索的上吊点坐标,同时根据已知的吊索下吊点坐标、下吊点设计竖向分力,计算出吊索上吊点横向分力PFDYJi和吊索的无应力长度SFDi;
步骤4,由第i索段末端节点平衡方程计算得到第i+1索段的起始端三向分力,重复步骤3.1-3.2,计算得到主索第i+1索段无应力长度SFi+1及第i+1索段横、竖桥向长度yi+1、zi+1和第i+1根吊索上吊点横向分力PFDYJi+1和第i+1根吊索的无应力长度SFDi+1,直至第n索段和第N根吊索;
步骤5,以主索末端的横、竖向坐标及跨内指定点竖向坐标为目标值,利用修正影响矩阵法修正主索始端三向分力HFXL、HFYL、VFL,获得迭代变量增量,反复迭代,直至目标变量误差小于允许值,计算得到成桥主索线形及各索段无应力长度。
进一步的,步骤1中根据空间抛物线理论估算主索始端三向分力HFX0、 HFY0、VF0初值如下:
式中,HFX0、HFY0、VF0分别为主索始端纵向、横向和竖向分力,l为主索跨径,w为主索在主索平面内沿跨长的等效均布荷载;f为主索在主索平面内垂度;yFDJi,zFDJi,yFDIi,zFDIi分别为吊索上、下吊点的横向及竖向坐标;Pyi, Pzi分别为吊索上端的横向和竖向分力,i为1~n的正整数,n为主索始端到末端索段数。
进一步的,步骤2中利用节线法迭代计算主索始端三向分力方法如下:
假设吊索数目为N,将主索分割成n段;
主索在竖向平面和水平面的平衡迭代式分别为:
其中,di为节点之间的水平长度;Ti为节点之间主索单元张力;li是节点间主索直线单元长度;Tx是主索张力水平分量;PFDi为第i根吊索的张力;hi为吊索长度;zFDJi为第i根吊索上吊点的竖向坐标;zFDIi为第i根吊索下吊点的竖向坐标;yFDIi为第i根吊索下吊点的横向坐标,yFDJi为吊索上吊点的横向坐标;WFi为主索节段自重;WFDi为将吊索力及吊索自重平均到主索上的集中荷载;
通过公式(5)、(6)求得主索节点坐标yFDIi、zFDIi,从而得到主索线形,再通过对主索系统整体受力分析,计算得到修正后的主索始端三向分力 将带入公式(5)、(6)再次迭代计算,求得未知量yFDIi+1、zFDIi+1,并计算得到修正后的主索始端三向分力反复迭代,直至相邻两次三向分力差值小于误差允许值:
进一步的,所述步骤2中的节线法存在以下假定:
1)吊索垂直于桥跨方向,且仅在横桥向有倾斜;
2)主索张力的水平分力在全跨相等;
3)主索与吊索连接点间的主索作直线考虑;
4)已知主索的两端锚点坐标、跨中垂度、桥面上的吊点位置和吊索设计竖向分力。
进一步的,计算主索第i索段横桥向长度yi和竖桥向长度zi及无应力长度 SFi方法如下:
式中,xi,yi,zi分别为主索第i索段纵桥向、横桥向和竖桥向长度;SFi,E和AF分别为主索第i索段无应力长度、弹性模量和截面面积;qF为主索索段的自重集度;HFXi-1,HFYi-1和VFi-1分别为主索第i索段起始端的纵向、横向和竖向分力; Hi-1为主索第i索段起始端的水平面合力;
其中,已知量包括主索第i索段纵桥向长度xi,主索弹性模量E和主索截面面积AF,主索索段自重集度qF,主索第i索段起始端的纵向分力HFXi-1、横向分力HFYi-1和竖向分力VFi-1。
进一步的,计算吊索上吊点横向分力PFDYJi PFDZJi和吊索的无应力长度SFDi方法如下:
主索各索段始端三向分力的递推关系如下:
HFXL、HFYL、VFL分别为主索始端纵向、横向和竖向分力;HFXi,HFYi和VFi分别为主索第i索段始端的纵向、横向和竖向分力;
对于每根吊索,已知量包括:吊索上下吊点横向距离yFDJi-yFDIi,吊索上下吊点竖向距离zFDJi-zFDIi,吊索自重集度qFD和吊索下吊点竖向分力PFDZJi。
进一步的,步骤5中利用修正影响矩阵法迭代计算得到主索线形及各索段无应力长度包括以下步骤:
步骤5.1,建立目标函数F(X),其中X=(HFXL,HFYL,VFL);
式中,HFXL、HFYL、VFL分别为主索始端纵向、横向和竖向分力;n为主索始端到末端索段数,m为主索始端到设计控制点索段数;yi,zi分别为主索第i 索段横桥向和竖桥向长度;Δy Δy和Δz分别为主索始端和末端的坐标差值,fz为设计控制点和主索始端的竖向坐标差值;ey为主索末端横向坐标误差,ez为主索末端竖向坐标误差,ef为设计控制点竖向坐标误差;
步骤5.2,求解修正影响矩阵J(X,R):
R=(ey,ez,ef),R1、R2、R3分别为R对应三个分量,E1、E2、E3分别为第 1、2、3个空间坐标向量;
步骤5.3,修正主索始端三向分力:
其中,k代表迭代次数;
反复迭代,直至误差ey,ez和ef小于允许值,得到修正后的主索始端三向分力HFXL、HFYL、VFL,重复步骤2,重新计算得到主索线形及各索段无应力长度。
本发明的有益效果为:
本发明所述用于大跨径悬索管道桥成桥主索线形及无应力长度计算方法是目前油气管道跨越行业尚无具体规范,设计尚未形成体系下悬索跨越主索系统设计的一次突破,为大跨径悬索管道桥成桥主索线形及无应力长度计算提供一套完整的设计理论和计算方法,能够结合设计控制点坐标和设计吊索力,最终设计出满足要求的成桥主索线形的计算方法。通过本发明所述计算方法能够解决了以往因主索始端三向力估值误差较大而导致算法不收敛的问题,采用了本方法开展设计和计算,可以大大提高大跨径悬索管道桥成桥主索线形算法的收敛性和适应性。
因此,通过本发明所述计算方法能够解决了以往设计不精、假设过多,计算误差太大的问题,采用本方法开展设计和计算,可以大大提高大跨径悬索管道桥成桥主索线形及无应力长度计算精度,从而进一步提高悬索管道桥主索设计合理性和降低施工难度。
本发明通过引入节线法对主索始端三向力迭代计算,使得主索始端三向力的估计精度大大提高,从而减少了影响矩阵法的迭代次数,不仅提高了计算效率,也使得算法的收敛性大大提高。
另外,传统的影响矩阵法使用时,施调量没有根据误差的大小来自动调整(如保持恒定值为1),从而使整个算法的效率很低,甚至无法收敛。传统方法只有在初值精度很高,而且三个初值相对比例在一定范围内时才可收敛。本发明引入Steffens-Newton法对影响矩阵进行修正,在计算影响矩阵时,迭代变量的改变值根据上一轮误差大小自动调整,即施调量采用误差本身R1E1,会自动根据结果的大小来不断调整,避免了整个迭代过程在计算影响矩阵时采用恒定的增量,使得迭代变量修正量的预测精度提高,从而使得整个算法的效率提高。另外,使用修正影响矩阵法使得其即使在某些弱条件下也至少是二阶收敛。同时采用差分矩阵,避免了梯度法的求偏导缺陷,也显著提高了运算速度。所以所采用修正的影响矩阵法计算空间主索线形对迭代初值的要求低,且具有较高的计算精度及收敛速度。
附图说明
图1为本发明所述大跨径悬索管道桥成桥主索线形计算方法;
图2为主索张力图解;
图3为节线法X-Z平面主索受力图;
图4为节线法Y-Z平面主索内力平衡示意图;
图5为空间索形主索计算图示;
图6空间索形吊索计算图示。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本发明对空间索形悬索桥主索缆索系统做如下假定:
(1)主索和吊索为小应变理想柔性索,其材料满足虎克定律,且忽略泊松效应。
(2)成桥状态下吊索仅沿横桥向倾斜,忽略施工过程中吊索纵桥向倾斜误差。
一种大跨径悬索管道桥成桥主索线形计算方法,包括以下步骤:
步骤1,根据空间抛物线理论估算主索始端三向分力初值HFX0、HFY0和VF0;
步骤2,以主索始端三向分力初值HFX0、HFY0、VF0为迭代初值,利用节线法迭代计算并得到修正后的主索始端三向分力HFXL、HFYL、VFL;
步骤3,以主索始端三向分力HFXL、HFYL、VFL为迭代初值,由始端至末端对主索各索段进行计算,包括:
步骤3.1,根据第i索段起始端三向分力、第i索段纵桥向长度xi和主索平衡方程计算出主索第i索段无应力长度SFi及第i索段横、竖桥向长度yi、zi;
步骤3.2,对于第i根吊索,由步骤3.1中所述主索平衡方程计算得到吊索的上吊点坐标,同时根据已知的吊索下吊点坐标、下吊点设计竖向分力,计算出吊索上吊点横向分力PFDYJi和吊索的无应力长度SFDi;
步骤4,由第i索段末端节点平衡方程计算得到第i+1索段的起始端三向分力,重复步骤3.1-3.2,计算得到主索第i+1索段无应力长度SFi+1及第i+1索段横、竖桥向长度yi+1、zi+1和第i+1根吊索上吊点横向分力PFDYJi+1和第i+1根吊索的无应力长度SFDi+1,直至第n索段和第N根吊索;
步骤5,以主索末端的横、竖向坐标及跨内指定点竖向坐标为目标值,利用修正影响矩阵法修正主索始端三向分力HFXL、HFYL、VFL,获得迭代变量增量,反复迭代,直至目标变量误差小于允许值,计算得到成桥主索线形及各索段无应力长度。
步骤1中根据空间抛物线理论估算主索始端三向分力HFX0、HFY0、VF0初值如下:
式中,HFX0、HFY0、VF0分别为主索始端纵向、横向和竖向分力,l为主索跨径,w为主索在主索平面内沿跨长的等效均布荷载;f为主索在主索平面内垂度;yFDJi,zFDJi,yFDIi,zFDIi分别为吊索上、下吊点的横向及竖向坐标;Pyi, Pzi分别为吊索上端的横向和竖向分力,i为1~n的正整数,n为主索始端到末端索段数。
步骤2中利用节线法迭代计算主索始端三向分力方法如下:
节线法采用Ohtsuki博士提出的计算索段平衡状态方程式,可以用于解决同时具有垂直和水平两个方向垂度的主索初始形状,参考图2,假定如下:
1)吊索垂直于桥跨方向,且仅在横桥向有倾斜;
2)主索张力的水平分力在全跨相等;
3)主索与吊索连接点间的主索作直线考虑;
4)已知主索的两端锚点坐标、跨中垂度、桥面上的吊点位置和设计吊索竖向分力。
一般地,把主索投影在垂直和水平面上,利用各个平面内的平衡条件分析,下面推导竖向和水平面上的分析过程。
(1)竖向平面分析
如图3所示,主索在竖向平面上的投影,假设吊索数目为N,将主索分割成n段;即N=n-1;主索第i个节点位置的平衡方程为:
在横桥向,Y-Z平面上的平衡方程如下:
主索在竖向平面的n-1平衡迭代式为:
(2)水平面内的分析
如图4所示,同竖向平面分析,水平面也可得到N-1个平衡迭代式:
其中,di为节点之间的水平长度;Ti为节点之间主索单元张力;li是节点间主索直线单元长度;Tx是主索张力水平分量;PFDi为第i根吊索的张力;hi为吊索长度;zFDJi为第i根吊索上吊点的竖向坐标;zFDIi为第i根吊索下吊点的竖向坐标;yFDIi为第i根吊索下吊点的横向坐标,yFDJi为吊索上吊点的横向坐标;WFi为主索节段自重;WFDi为将吊索力及吊索自重平均到主索上的集中荷载。
通过公式(5)、(6)及跨中垂度或设计控制点的竖向坐标ze(或者横向坐标ye):
通过公式(5)、(6)求得主索节点坐标yFDIi、zFDIi,从而得到主索线形,再通过对主索系统整体受力分析,计算得到修正后的主索始端三向分力 将带入公式(5)、(6)再次迭代计算,求得未知量yFDIi+1、zFDIi+1,并计算得到修正后的主索始端三向分力反复迭代,直至相邻两次三向分力差值小于误差允许值:
计算主索第i索段横桥向长度yi和竖桥向长度zi及无应力长度SFi方法如下,如图5所示,对空间缆索而言,由于节间只有自重作用,索段总在一个铅垂面上,只是各索段在水平面上的投影与桥轴线的夹角不同。因此,各索段在各自铅垂面内满足:
式中,xi,yi,zi分别为主索第i索段纵桥向、横桥向和竖桥向长度;SFi,E和AF分别为主索第i索段无应力长度、弹性模量和截面面积;qF为主索索段的自重集度;HFXi-1,HFYi-1和VFi-1分别为主索第i索段起始端的纵向、横向和竖向分力; Hi-1为主索第i索段起始端的水平面合力;
其中,已知量包括主索第i索段纵桥向长度xi,主索弹性模量E和主索截面面积AF,主索索段自重集度qF,主索第i索段起始端的纵向分力HFXi-1、横向分力HFYi-1和竖向分力VFi-1。
进一步的,计算吊索上吊点横向分力PFDYJi和吊索的无应力长度SFDi方法如下:
空间索形的吊索沿横桥向倾斜,只有将吊索看成弹性悬链线才能保证计算精度。如图6所示,主索各索段始端三向分力的递推关系如下:
HFXL、HFYL、VFL分别为主索始端纵向、横向和竖向分力;HFXi,HFYi和VFi分别为主索第i索段始端的纵向、横向和竖向分力;
对于每根吊索,已知量包括:吊索上下吊点横向距离yFDJi-yFDIi,吊索上下吊点竖向距离zFDJi-zFDIi,吊索自重集度qFD和吊索下吊点竖向分力PFDZJi。
公式(7)-公式(11)构成了空间索形主索的平衡方程。
步骤5中利用修正影响矩阵法迭代计算得到主索线形及各索段无应力长度包括以下步骤:
步骤5.1,建立目标函数F(X),其中X=(HFXL,HFYL,VFL);
式中,HFXL、HFYL、VFL分别为主索始端纵向、横向和竖向分力;n为主索始端到末端索段数,m为主索始端到设计控制点索段数;yi,zi分别为主索第i 索段横桥向和竖桥向长度;Δy Δy和Δz分别为主索始端和末端的坐标差值,fz为设计控制点和主索始端的竖向坐标差值;ey为主索末端横向坐标误差,ez为主索末端竖向坐标误差,ef为设计控制点竖向坐标误差;
步骤5.2,求解修正影响矩阵J(X,R):
R=(ey,ez,ef),R1、R2、R3分别为R对应三个分量,E1、E2、E3分别为第 1、2、3个空间坐标向量;
步骤5.3,修正主索始端三向分力:
其中,k代表迭代次数;
反复迭代,直至误差ey,ez和ef小于允许值,得到修正后的主索始端三向分力HFXL、HFYL、VFL,重复步骤2,重新计算得到主索线形及各索段无应力长度。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种大跨径悬索管道桥成桥主索线形计算方法,其特征在于,包括以下步骤:
步骤1,根据空间抛物线理论估算主索始端三向分力初值HFX0、HFY0和VF0;
步骤2,以主索始端三向分力初值HFX0、HFY0、VF0为迭代初值,利用节线法迭代计算并得到修正后的主索始端三向分力HFXL、HFYL、VFL:
假设吊索数目为N,将主索分割成n段;
主索在竖向平面和水平面的平衡迭代式分别为:
其中,di为节点之间的水平长度;Ti为节点之间主索单元张力;li是节点间主索直线单元长度;Tx是主索张力水平分量;PFDi为第i根吊索的张力;hi为吊索长度;zFDJi为第i根吊索上吊点的竖向坐标;zFDIi为第i根吊索下吊点的竖向坐标;yFDIi为第i根吊索下吊点的横向坐标,yFDJi为吊索上吊点的横向坐标;WFi为主索节段自重;WFDi为将吊索力及吊索自重平均到主索上的集中荷载;
通过公式(5)、(6)求得主索节点坐标yFDIi、zFDIi,从而得到主索线形,再通过对主索系统整体受力分析,计算得到修正后的主索始端三向分力将带入公式(5)、(6)再次迭代计算,求得未知量yFDIi+1、zFDIi+1,并计算得到修正后的主索始端三向分力反复迭代,直至相邻两次三向分力差值小于误差允许值:
其中,所述步骤2中的节线法存在以下假定:
1)吊索垂直于桥跨方向,且仅在横桥向有倾斜;
2)主索张力的水平分力在全跨相等;
3)主索与吊索连接点间的主索作直线考虑;
4)已知主索的两端锚点坐标、跨中垂度、桥面上的吊点位置和吊索设计竖向分力;
步骤3,以主索始端三向分力HFXL、HFYL、VFL为迭代初值,由始端至末端对主索各索段进行计算,包括:
步骤3.1,根据第i索段起始端三向分力、第i索段纵桥向长度xi和主索平衡方程计算出主索第i索段无应力长度SFi及第i索段横、竖桥向长度yi、zi;
步骤3.2,对于第i根吊索,由步骤3.1中所述主索平衡方程计算得到吊索的上吊点坐标,同时根据已知的吊索下吊点坐标、下吊点设计竖向分力,计算出吊索上吊点横向分力PFDYJi和吊索的无应力长度SFDi;
步骤4,由第i索段末端节点平衡方程计算得到第i+1索段的起始端三向分力,重复步骤3.1-3.2,计算得到主索第i+1索段无应力长度SFi+1及第i+1索段横、竖桥向长度yi+1、zi+1和第i+1根吊索上吊点横向分力PFDYJi+1和第i+1根吊索的无应力长度SFDi+1,直至第n索段和第N根吊索;
步骤5,以主索末端的横、竖向坐标及跨内指定点竖向坐标为目标值,利用修正影响矩阵法修正主索始端三向分力HFXL、HFYL、VFL,获得迭代变量增量,反复迭代,直至目标变量误差小于允许值,计算得到成桥主索线形及各索段无应力长度。
3.根据权利要求1所述的大跨径悬索管道桥成桥主索线形计算方法,其特征在于,计算主索第i索段横桥向长度yi和竖桥向长度zi及无应力长度SFi方法如下:
式中,xi,yi,zi分别为主索第i索段纵桥向、横桥向和竖桥向长度;SFi,E和AF分别为主索第i索段无应力长度、弹性模量和截面面积;qF为主索索段的自重集度;HFXi-1,HFYi-1和VFi-1分别为主索第i索段起始端的纵向、横向和竖向分力;Hi-1为主索第i索段起始端的水平面合力;
其中,已知量包括主索第i索段纵桥向长度xi,主索弹性模量E和主索截面面积AF,主索索段自重集度qF,主索第i索段起始端的纵向分力HFXi-1、横向分力HFYi-1和竖向分力VFi-1。
4.根据权利要求1所述的大跨径悬索管道桥成桥主索线形计算方法,其特征在于,步骤5中利用修正影响矩阵法迭代计算得到主索线形及各索段无应力长度包括以下步骤:
步骤5.1,建立目标函数F(X),其中X=(HFXL,HFYL,VFL);
式中,HFXL、HFYL、VFL分别为主索始端纵向、横向和竖向分力;n为主索始端到末端索段数,m为主索始端到设计控制点索段数;yi,zi分别为主索第i索段横桥向和竖桥向长度;Δy Δy和Δz分别为主索始端和末端的坐标差值,fz为设计控制点和主索始端的竖向坐标差值;ey为主索末端横向坐标误差,ez为主索末端竖向坐标误差,ef为设计控制点竖向坐标误差;
步骤5.2,求解修正影响矩阵J(X,R):
R=(ey,ez,ef),R1、R2、R3分别为R对应三个分量,E1、E2、E3分别为第1、2、3个空间坐标向量;
步骤5.3,修正主索始端三向分力:
其中,k代表迭代次数;
反复迭代,直至误差ey,ez和ef小于允许值,得到修正后的主索始端三向分力HFXL、HFYL、VFL,重复步骤2,重新计算得到主索线形及各索段无应力长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010899711.2A CN112035928B (zh) | 2020-08-31 | 2020-08-31 | 大跨径悬索管道桥成桥主索线形计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010899711.2A CN112035928B (zh) | 2020-08-31 | 2020-08-31 | 大跨径悬索管道桥成桥主索线形计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112035928A CN112035928A (zh) | 2020-12-04 |
CN112035928B true CN112035928B (zh) | 2022-04-15 |
Family
ID=73587482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010899711.2A Active CN112035928B (zh) | 2020-08-31 | 2020-08-31 | 大跨径悬索管道桥成桥主索线形计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112035928B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115130350B (zh) * | 2022-07-06 | 2023-06-20 | 广西路桥工程集团有限公司 | 一种连续索的非线性滑移计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018001147A1 (zh) * | 2016-06-29 | 2018-01-04 | 深圳市智能机器人研究院 | 一种基于优化张紧弦模型的桥索监测方法及系统 |
CN108491635A (zh) * | 2018-03-26 | 2018-09-04 | 东南大学 | 一种悬索桥吊杆力和主缆线形联合计算方法 |
CN108875283A (zh) * | 2018-08-09 | 2018-11-23 | 中交基础设施养护集团有限公司 | 一种空间索面悬索桥空缆线形的精确计算方法 |
CN110990913A (zh) * | 2019-11-13 | 2020-04-10 | 中铁大桥勘测设计院集团有限公司 | 一种采用水平面内合力控制的边跨找形方法 |
CN111523172A (zh) * | 2020-05-11 | 2020-08-11 | 重庆交通大学 | 一种空间异形索面悬索桥主缆成桥线形分析方法 |
-
2020
- 2020-08-31 CN CN202010899711.2A patent/CN112035928B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018001147A1 (zh) * | 2016-06-29 | 2018-01-04 | 深圳市智能机器人研究院 | 一种基于优化张紧弦模型的桥索监测方法及系统 |
CN108491635A (zh) * | 2018-03-26 | 2018-09-04 | 东南大学 | 一种悬索桥吊杆力和主缆线形联合计算方法 |
CN108875283A (zh) * | 2018-08-09 | 2018-11-23 | 中交基础设施养护集团有限公司 | 一种空间索面悬索桥空缆线形的精确计算方法 |
CN110990913A (zh) * | 2019-11-13 | 2020-04-10 | 中铁大桥勘测设计院集团有限公司 | 一种采用水平面内合力控制的边跨找形方法 |
CN111523172A (zh) * | 2020-05-11 | 2020-08-11 | 重庆交通大学 | 一种空间异形索面悬索桥主缆成桥线形分析方法 |
Non-Patent Citations (1)
Title |
---|
悬索桥空间缆索主缆线形的计算方法;齐东春等;《铁道建筑》;20130420(第04期);13-16 * |
Also Published As
Publication number | Publication date |
---|---|
CN112035928A (zh) | 2020-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111967085B (zh) | 大跨径悬索管道桥成桥风索线形计算方法 | |
CN112035927B (zh) | 基于修正影响矩阵的悬索管道桥成桥风索线形计算方法 | |
CN111523172B (zh) | 一种空间异形索面悬索桥主缆成桥线形分析方法 | |
CN108460229B (zh) | 连续桥面结构桥梁拉索内力的调整方法 | |
CN108875283A (zh) | 一种空间索面悬索桥空缆线形的精确计算方法 | |
CN111783201B (zh) | 一种三跨自锚式悬索桥动力特性的快速分析方法 | |
WO2022165983A1 (zh) | 一种确定三塔斜拉桥塔梁纵向约束刚度及优化基础的方法 | |
CN107142835A (zh) | 无应力状态法的肋拱桥拱肋节段分肋安装施工控制方法 | |
CN111222193B (zh) | 一种确定悬索桥主缆线形的方法 | |
CN107622174B (zh) | 基于影响矩阵修正的斜拉桥拉索张拉施工闭环控制方法 | |
CN113255029B (zh) | 一种活载作用下悬索桥结构变形及内力的确定方法 | |
CN112035928B (zh) | 大跨径悬索管道桥成桥主索线形计算方法 | |
CN111666615A (zh) | 一种基于有限元的悬索桥空间缆索找形方法 | |
CN115357965B (zh) | 自锚式悬索桥及其成桥线形确定方法 | |
CN108509710B (zh) | 一种平行双幅桥静风稳定分析方法 | |
CN113468632B (zh) | 一种偏心活载作用下的悬索桥全桥响应的确定方法 | |
CN115510535A (zh) | 一种悬索桥空缆缆形与索鞍预偏量的联合求解方法 | |
CN111931282B (zh) | 基于未知荷载系数法一次张拉斜拉扣挂索力计算方法 | |
CN116226986A (zh) | 一种拉索无应力长度计算方法 | |
CN112035929B (zh) | 一种悬索管道桥成桥风索线形计算方法 | |
CN111753435A (zh) | 基于分段悬链线和索力连续算法的缆索吊装系统计算方法 | |
CN112035926A (zh) | 基于节线法的大跨径悬索管道桥成桥风索线形计算方法 | |
CN115455547B (zh) | 一种考虑静动力效应的悬索桥施工抗风分析优化方法 | |
CN114417679A (zh) | 一种斜拉桥钢主梁无应力制造线形计算方法 | |
CN112231805B (zh) | 一种基于大数据的组合梁斜拉桥索力估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |