CN108491635A - 一种悬索桥吊杆力和主缆线形联合计算方法 - Google Patents

一种悬索桥吊杆力和主缆线形联合计算方法 Download PDF

Info

Publication number
CN108491635A
CN108491635A CN201810250670.7A CN201810250670A CN108491635A CN 108491635 A CN108491635 A CN 108491635A CN 201810250670 A CN201810250670 A CN 201810250670A CN 108491635 A CN108491635 A CN 108491635A
Authority
CN
China
Prior art keywords
towing rope
main
main push
bridge
cable saddle
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.)
Granted
Application number
CN201810250670.7A
Other languages
English (en)
Other versions
CN108491635B (zh
Inventor
张文明
刘钊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201810250670.7A priority Critical patent/CN108491635B/zh
Publication of CN108491635A publication Critical patent/CN108491635A/zh
Application granted granted Critical
Publication of CN108491635B publication Critical patent/CN108491635B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

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

本发明公开了一种悬索桥吊杆力和主缆线形联合计算方法,该方法包括以下步骤:(1)以悬索桥中跨左切点和各个吊点为坐标原点建立坐标系,表示各段主缆的悬链线;(2)假定一组吊杆力,即为每个所述吊杆力赋予一个初值;(3)利用解析法计算中跨主缆的成桥线形和内力;(4)根据步骤(3)计算的数值建立所述悬索桥有限元模型,计算吊杆力,用欧几里得范数检验与所述吊杆力假定值的误差,若所述误差超过限值,带着新算的吊杆力返回步骤(2),否则得到所述悬索桥的中跨主缆线形和精确的吊杆力;(5)利用所述中跨主缆内力的水平分量作为已知条件计算边跨主缆线形。该发明可以获得悬索桥精确的吊杆力及相应的主缆成桥线形,其中主缆成桥线形中包含了切点的精确位置。

Description

一种悬索桥吊杆力和主缆线形联合计算方法
技术领域
本发明涉及一种悬索桥设计和施工控制计算方法,具体涉及一种悬索桥吊杆力和主缆成桥线形联合计算方法。
背景技术
悬索桥的跨越能力在常见桥型中独占鳌头。随着桥梁设计和施工水平的不断提高,悬索桥跨度记录不断被刷新,如已建成的日本明石海峡大桥(主跨1991m)和中国西堠门大桥(主跨1665m),以及正在施工的意大利墨西拿海峡大桥(主跨3300m)和中国的杨泗港长江大桥(主跨1700m)。
悬索桥与斜拉桥不同,不能通过施工阶段的跟踪调整来实现设计线形,它的主缆线形必须通过可靠的计算参数,采用精确的计算方法计算出设计线形,同时通过严格的施工控制保证最终的成桥线形与设计线形相符。成桥状态主缆的设计线形和吊杆力是计算施工控制参数的先决条件。
悬索桥的吊杆力在成桥状态是不均匀的,尤其是桥塔附近的几根吊杆,不能简单地用恒载集度按照吊杆间距直接分配给各吊杆,也不能用多点刚性支撑连续梁的支反力代替吊杆力。吊杆力的求解精度直接影响到主缆成桥线形、空缆线形、主缆无应力长度和吊杆无应力长度。如果把加劲梁比拟成一个多点弹性支撑连续梁,支撑反力(即吊杆力)取决于多点支撑刚度的比例关系。每个点的支撑刚度取决于吊杆的刚度和主缆对该点竖向支撑刚度的贡献。对于悬索桥成桥状态来说,一般在设计时只知道全桥各跨跨径、主缆矢高、主缆跨中标高、各索鞍理论交点坐标、吊杆间距以及加劲梁线形等。而主缆各吊点的标高是未知的,因此无法确定吊杆长度。在主缆线形和吊杆长度都未确定的情况下,加劲梁的多点弹性支撑刚度是未知的。因此,悬索桥吊杆力和主缆线形之间需要迭代求解。
有限元法和解析法是目前悬索桥线形和内力计算的常用方法,两者各有优缺点。利用有限元法计算吊杆力是方便的,但是鞍座的模拟是困难的,无论采用多个刚性的杆件还是用将鞍座和索融为一体的鞍座-索单元都不能让主缆与鞍座的切点随着主缆的受力情况而不断改变。解析法计算过程明了,可以方便地处理主缆与索鞍的切点问题,但是不便于计算吊杆力。
发明内容
发明目的:为了克服现有技术的不足,本发明提供一种悬索桥吊杆力和主缆成桥线形联合计算方法,该方法可以精确计算悬索桥吊杆力和主缆成桥线形。
技术方案:本发明所述的悬索桥吊杆力和主缆成桥线形联合计算方法,该方法包括以下步骤:
(1)以悬索桥中跨左切点和各个吊点为坐标原点建立坐标系,表示各段主缆的悬链线;
(2)假定一组吊杆力,即为每个所述吊杆力赋予一个初值;
(3)利用解析法计算中跨主缆的成桥线形和内力;
(4)根据步骤(3)计算的数值建立所述悬索桥有限元模型,计算吊杆力,用欧几里得范数检验与所述吊杆力假定值的误差,若所述误差超过限值,将新计算的吊杆力返回步骤(2)迭代计算,否则,计算结束,得到所述悬索桥的中跨主缆线形和精确的吊杆力;
(5)利用所述中跨主缆内力的水平分量作为已知条件计算边跨主缆线形。
优选的,步骤(1)中,所述坐标系以x轴水平向右,y轴竖直向下,任一段所述主缆的悬链线方程可表达为:
其中,c=-H/q,H为成桥状态主缆水平力,q为成桥状态主缆自重荷载集度,ai和bi是第i段悬链线方程的参数。
优选的,步骤(3)中,所述中跨主缆的成桥线形和内力的计算方法包括:
根据主缆跨中高程、切点高差的闭合和最右侧一段悬链线主缆的水平投影长度ln+1满足设计要求等三个条件分别建立三个方程:
其中,m是左切点O0到跨中点Om的主缆段数;n是吊点数量;Δhi是任一段主缆两端点的高差;是左切点O0与跨中点Om之间的高差;是左索鞍上切点O0与右索鞍上切点On+1之间的高差;是最右侧一根吊杆与右索鞍实际顶点D'之间的设计水平距离;Δ1是右侧索鞍上切点On+1与实际顶点D'之间的水平距离;
上述三个方程可表达成以下三个函数形式:
f1(H,a1,ln+1)=0
f2(H,a1,ln+1)=0
f3(H,a1,ln+1)=0
在利用非线性规划求解法求解方程时,以H、a1和ln+1为变量,目标函数为:
f1 2+f2 2+f3 2=0
解所述函数可得到,成桥状态主缆水平力H、第一段悬链线主缆线形方程中系数a1、最右侧一段悬链线主缆的水平投影长度ln+1;进而得到中跨主缆的成桥线形和内力。
优选的,所述成桥线形和内力包括所述中跨主缆的各切点位置、各吊点高程和各索段的应力和应变。
优选的,所述任一段主缆两端点的高差Δhi可表达为:
其中,li为第i段主缆的水平投影长度,ai是第i段悬链线方程的系数;
所述左索鞍上切点O0与右索鞍上切点On+1之间的高差表达为:
其中,hC'是圆心C'在大地坐标系中的设计高程;R'是右索鞍在主缆截面中心处的圆弧半径;
所述右侧索鞍上切点On+1与实际顶点D'之间的水平距离Δ1表达为:
其中,Δ2是右索鞍圆心C'与实际顶点D'的水平距离。
优选的,所述第i段悬链线方程的系数ai和第i+1段悬链线方程的系数ai+1的关系为:
其中,1≤i≤n+1,li为第i段主缆的水平投影长度,pi是第i段的吊杆力,H为成桥状态主缆水平力。
优选的,步骤(4)中,所述有限元模型包括模拟中跨的主缆、吊杆和主梁,主缆的端点分别是主缆在左、右两个主索鞍上的切点,用铰接约束,所述主缆用索单元模拟,并输入各索段的初应变,所述主梁用梁单元模拟,其节点坐标依据桥面设计线形确定。
优选的,步骤(4)中,所述欧几里得范数检验与吊杆力假定值的误差公式表示为:
其中,Pi,k是用有限元法计算出的吊杆力;Pi,k-1是吊杆力假定值;ε是误差限值。
优选的,步骤(5)中,所述边跨主缆线形的计算方法为:
所述边跨主缆与主索鞍的切点为Q,与散索鞍的切点为M,设Q、M两点之间的水平距离为s,以Q点为原点建立坐标系,x轴水平向右,y轴竖直向下,方程表达为:
其中,c=-H/q,H为成桥状态主缆水平力,q为成桥状态主缆自重荷载集度;
根据主索鞍IP点与散索鞍IP点之间的水平距离和高差得出:
Δ23+s+Δ4=L
h1+h2+h3=h
其中,Δ2是主索鞍圆心C'与桥塔中心线之间的水平距离,为已知量;Δ3是主索鞍圆心C'与边跨主缆在主索鞍上的切点Q之间的水平距离;Δ4是主缆在散索鞍上的切点M与散索鞍IP点之间的水平距离;L是悬索桥边跨跨径,h1是主索鞍IP点与边跨主缆在主索鞍上的切点Q之间高差;h2是边跨主缆在主索鞍上的切点Q与在散索鞍上的切点M之间的高差;h3是切点M与散索鞍IP点之间的高差;h是主索鞍IP点与散索鞍IP点之间的高差。
优选的,所述主索鞍表示为:
Δ3=R'sinβ
其中,R'是所述主索鞍半径;β是所述直线C'Q与竖直线的夹角;
所述主索鞍IP点与边跨主缆在主索鞍上的切点Q之间高差h1表示为:
h1=(Δ23)tanβ
所述边跨主缆,两端点高差h2表示为:
所述散索鞍M、N两点之间的水平距离Δ4和高差h3分别表达为:
其中,r是所述散索鞍半径,ω是散索鞍圆心P点与IP点连线与竖直线的夹角,是所述散索鞍圆心P点与切点M连线与竖直线的夹角,从而得到边跨主缆的线形和切点位置。
有益效果:本发明与现有技术相比,其显著优点是通过悬索桥吊杆力和主缆线形联合计算方法可以获得悬索桥精确的吊杆力及相应的主缆成桥线形,其中主缆成桥线形中包含了切点的精确位置。
附图说明
图1为本发明所述的成桥状态中跨的主缆线形示意图;
图2为本发明所述的成桥状态左索鞍布置图;
图3为本发明所述的成桥状态右索鞍布置图;
图4为本发明所述的吊点处力的平衡关系图;
图5为本发明所述的有限元模型示意图;
图6为本发明所述的边跨主缆线形示意图。
具体实施方式
本发明所述的方法利用解析法计算主缆线形,利用有限元法计算吊杆力,两者互为对方的输入条件,迭代计算。
由于吊杆力的作用,成桥状态的主缆线形是带折角的多段悬链线,即相邻吊点之间的线形是悬链线。计算成桥主缆线形一般先计算中跨。如图1所示,分别以左切点和各吊点为坐标原点建立坐标系,x轴水平向右,y轴竖直向下,那么任一段主缆的悬链线方程可表达为:
式中,c=-H/q,H为成桥状态主缆水平力(kN),q为成桥状态主缆自重荷载集度(kN/m),ai和bi是第i段悬链线方程的参数。
根据边界条件y(0)=0,可得bi=-ccoshai。代入上式,任一段主缆的悬链线方程可改写为
如图2和图3所示,D和D'所在的竖直线为桥塔中心线,该线图2中左侧区域为边跨,右侧区域为中跨,图3中左侧区域表示中跨,右侧区域表示边跨,可根据主缆跨中高程、切点高差的闭合和最右侧一段悬链线主缆的水平投影长度ln+1满足设计要求等三个条件分别建立三个方程:
式中,m是左切点O0到跨中点Om的主缆段数;n是吊点数量;Δhi是任一段主缆两端点的高差;是左切点O0与跨中点Om之间的高差;是左索鞍上切点O0与右索鞍上切点On+1之间的高差;是最右侧一根吊杆与右索鞍实际顶点D'之间的设计水平距离;Δ1是右侧索鞍上切点On+1与实际顶点D'之间的水平距离。
这三个方程含了三个未知量:成桥状态主缆水平力H、第一段悬链线主缆线形方程中系数a1、最右侧一段悬链线主缆的水平投影长度ln+1。接下来,可将方程中的参数表达成这三个未知量的函数。
任一段主缆两端点的高差可表达为
式中,li为任一段主缆的水平投影长度。
对于左索鞍,已知圆心C在大地坐标系中的设计高程hC,那么切点O0在大地坐标系中的高程可表达为
式中,R是左索鞍在主缆截面中心处的圆弧半径;θ是切点O0和圆心C的连线与竖直线BC的夹角。
由于
所以有
cosθ=secha1 (7)
将式(7)代入式(5),得
已知主缆跨中点Om在大地坐标系中的设计高程那么左切点O0与跨中点Om之间的高差可表达为
对于右索鞍,已知圆心C'在大地坐标系中的设计高程hC',那么切点On+1在大地坐标系中的高程可表达为
式中,R'是右索鞍在主缆截面中心处的圆弧半径;θ'是切点On+1和圆心C'的连线与竖直线B'C'的夹角。
由于
所以有
将式(12)代入式(10),得
那么,左索鞍上切点O0与右索鞍上切点On+1之间的高差为:
右侧索鞍上切点On+1与实际顶点D'之间的水平距离Δ1可表达为
式中,Δ2是右索鞍圆心C'与实际顶点D'的水平距离。
在任一吊点处,将主缆轴向拉力分解为水平分力和竖向分力,如图3所示。利用竖向力的平衡可得:
Htanδ=Htanφ+Pi (17)
其中,Pi是吊杆力;δ和φ分别为吊点左右两侧主缆索段的倾角。
将tanδ=sinh(li/c+ai)和tanφ=sinhai+1代入上式可得
Hsinh(li/c+ai)=Hsinhai+1+Pi (18)
从而可得
将式(4)、(9)、(15)、(16)和(19)代入方程组(3),可得到三个耦合的方程,然后利用非线性规划求解法求解。方程组(3)的三个方程可表达成以下三个函数形式:
f1(H,a1,ln+1)=0 (20.1)
f2(H,a1,ln+1)=0 (20.2)
f3(H,a1,ln+1)=0 (20.3)
在利用非线性规划求解法求解方程时,以H、a1和ln+1为变量,以下式为目标函数
f1 2+f2 2+f3 2=0 (21)
可解得三个未知量:成桥状态主缆水平力H、第一段悬链线主缆线形方程中系数a1、最右侧一段悬链线主缆的水平投影长度ln+1。进而可以获得中跨主缆的成桥线形和内力,具体包括:切点位置、各吊点高程、各索段的应力和应变。
边跨主缆成桥线形的计算方法与中跨基本相同,只是已知条件略有不同。计算中跨时已知矢跨比和跨中主缆高程,而计算边跨时未知。但边跨主缆的水平力H可根据塔顶主鞍座的平衡条件求得。一般假定桥塔不承受主缆水平力,所以计算边跨时,认为其主缆水平力与中跨相等。因此边跨主缆线形的计算更加方便。
在上述用解析法计算主缆线形的过程中,需用到吊杆力。但是吊杆力是未知的,可先假定一组吊杆力。至于假定值是否准确,还需用有限元法进行检验。
如图4所示,在建立有限元模型时,只需模拟中跨的主缆、吊杆和主梁,不用模拟桥塔,这样可以避免桥塔弹性压缩带来的麻烦。主缆的两端点分别是主缆在左、右两个主索鞍上的切点,用铰接约束。主缆上各节点坐标就是各吊点的坐标,来自于上述用解析法计算主缆线形的结果。用索单元模拟主缆,并输入各索段的初应变。用梁单元模拟加劲梁,其节点坐标依据桥面设计线形确定。依据梁端支座布置情况,在有限元模型中约束梁端相应的自由度。在主缆吊点和加劲梁上吊点之间用索单元模拟吊杆,吊杆中可以输入初应变,这样可以避免加劲梁在自重作用下产生过大的挠度变形。初应变可以根据假定的吊杆力求出。经静力计算,可获得悬索桥在自重作用下的吊杆力。新算出的吊杆力可能与先前的假定值有差异,于是需要迭代计算。迭代步骤如下:
(1)假定一组吊杆力;
(2)利用解析法求解中跨主缆的成桥线形和内力,具体包括:切点位置、各吊点高程、各索段的应力和应变。
(3)利用有限元法计算吊杆力,然后用欧几里得范数检验与吊杆力假定值的误差:
式中,Pi,k是用有限元法计算出的吊杆力;Pi,k-1是吊杆力假定值;ε是误差限值。
如果吊杆力的欧几里得范数不满足要求,把Pi,k作为吊杆力假定值,返回步骤(2);否则,计算结束,求出了悬索桥的主缆线形和吊杆力。
在中跨主缆线形求解完成之后进行。假定桥塔不承受主缆水平力,所以计算边跨时,认为其主缆水平力与中跨相等,即中跨求得的主缆水平力可以作为边跨主缆线形求解的已知条件。悬索桥边跨一般没有吊杆。本实施例以右边跨为例。
如图6所示,D'所在的竖直线为桥塔中心线,该线左侧区域表示中跨,右侧区域表示边跨,边跨主缆与主索鞍的切点为Q,与散索鞍的切点为M,假设Q、M两点之间的水平距离为s。以Q点为原点建立坐标系,x轴水平向右,y轴竖直向下,那么边跨主缆(Q、M两点之间)的悬链线方程可表达为:
式中,c=-H/q,H为成桥状态主缆水平力(kN),q为成桥状态主缆自重荷载集度(kN/m)。由于H和q均已知,因此c是已知量。
根据边界条件y(0)=0,可得b=-ccosha。代入上式,边跨主缆的悬链线方程可改写为
可根据主索鞍IP点(T)与散索鞍IP点(N)之间的水平距离和高差列出以下方程:
Δ23+s+Δ4=L (25.1)
h1+h2+h3=h (25.2)
式中,Δ2是主索鞍圆心C'与桥塔中心线之间的水平距离,为已知量;Δ3是主索鞍圆心C'与边跨主缆在主索鞍上的切点Q之间的水平距离;Δ4是主缆在散索鞍上的切点M与散索鞍IP点(N)之间的水平距离;L是悬索桥边跨跨径,即主索鞍IP点(T)与散索鞍IP点(N)之间的水平距离,为已知量。h1是主索鞍IP点(T)与边跨主缆在主索鞍上的切点Q之间高差;h2是边跨主缆在主索鞍上的切点Q与在散索鞍上的切点M之间的高差;h3是切点M与散索鞍IP点(N)之间的高差;h是主索鞍IP点(T)与散索鞍IP点(N)之间的高差。
在该方程组中,暗含两个未知量:a和s。为了求解方程组,需先将方程组中的各参数表达成未知量a和s的函数。
在主索鞍上,
Δ3=R'sinβ (26)
式中,R'是主索鞍半径;β是直线C'Q与竖直线的夹角。
由于
所以有
sinβ=tanha (28)
将式(28)代入式(26),可得
Δ3=R'tanha (29)
主索鞍IP点(T)与边跨主缆在主索鞍上的切点Q之间高差h1可表述为
h1=(Δ23)tanβ=(Δ2+R'tanha)sinha (30)
对于边跨主缆(Q、M两点之间),两端点高差h2可表示为
对于散索鞍,M、N两点之间的水平距离Δ4和高差h3可分别表达为
式中,r是散索鞍半径;ω是散索鞍圆心P点与IP点(N)连线与竖直线的夹角,是已知量;是散索鞍圆心P点与切点M连线与竖直线的夹角。
由于
所以有
分别代入式(32)和(33),可得
将上述式(29)和(38)代入式(25.1),将式(30)、(31)和(39)代入式(25.2),然后求解方程组,可解得未知量a和s。因此,获得边跨主缆的线形和切点位置。

Claims (10)

1.一种悬索桥吊杆力和主缆线形联合计算方法,其特征在于,该方法包括以下步骤:
(1)以悬索桥中跨左切点和各个吊点为坐标原点建立坐标系,表示各段主缆的悬链线;
(2)假定一组吊杆力,即为每个所述吊杆力赋予一个初值;
(3)利用解析法计算中跨主缆的成桥线形和内力;
(4)根据步骤(3)计算的数值建立所述悬索桥有限元模型,计算吊杆力,用欧几里得范数检验与所述吊杆力假定值的误差,若所述误差超过限值,将新计算的吊杆力返回步骤(2)迭代计算,否则,计算结束,得到所述悬索桥的中跨主缆线形和精确的吊杆力;
(5)利用所述中跨主缆内力的水平分量作为已知条件计算边跨主缆线形。
2.根据权利要求1所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,步骤(1)中,所述坐标系以x轴水平向右,y轴竖直向下,任一段所述主缆的悬链线方程可表达为:
其中,c=-H/q,H为成桥状态主缆水平力,q为成桥状态主缆自重荷载集度,ai和bi是第i段悬链线方程的参数。
3.根据权利要求1所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,步骤(3)中,所述中跨主缆的成桥线形和内力的计算方法包括:
根据主缆跨中高程、切点高差的闭合和最右侧一段悬链线主缆的水平投影长度ln+1满足设计要求等三个条件分别建立三个方程:
其中,m是左切点O0到跨中点Om的主缆段数;n是吊点数量;Δhi是任一段主缆两端点的高差;是左切点O0与跨中点Om之间的高差;是左索鞍上切点O0与右索鞍上切点On+1之间的高差;是最右侧一根吊杆与右索鞍实际顶点D'之间的设计水平距离;Δ1是右侧索鞍上切点On+1与实际顶点D'之间的水平距离;
上述三个方程可表达成以下三个函数形式:
f1(H,a1,ln+1)=0
f2(H,a1,ln+1)=0
f3(H,a1,ln+1)=0
在利用非线性规划求解法求解方程时,以H、a1和ln+1为变量,目标函数为:
f1 2+f2 2+f3 2=0
解所述函数可得到,成桥状态主缆水平力H、第一段悬链线主缆线形方程中系数a1、最右侧一段悬链线主缆的水平投影长度ln+1;进而得到中跨主缆的成桥线形和内力。
4.根据权利要求3所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,所述成桥线形和内力包括所述中跨主缆的各切点位置、各吊点高程和各索段的应力和应变。
5.根据权利要求3所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,所述任一段主缆两端点的高差Δhi表示为:
其中,li为第i段主缆的水平投影长度,ai是第i段悬链线方程的系数;
所述左索鞍上切点O0与右索鞍上切点On+1之间的高差表达为:
其中,hC'是圆心C'在大地坐标系中的设计高程;R'是右索鞍在主缆截面中心处的圆弧半径;
所述右侧索鞍上切点On+1与实际顶点D'之间的水平距离Δ1表达为:
其中,Δ2是右索鞍圆心C'与实际顶点D'的水平距离。
6.根据权利要求5所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,所述第i段悬链线方程的系数ai和第i+1段悬链线方程的系数ai+1的关系为:
其中,1≤i≤n+1,li为第i段主缆的水平投影长度,pi是第i段的吊杆力,H为成桥状态主缆水平力。
7.根据权利要求1所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,步骤(4)中,所述有限元模型包括模拟中跨的主缆、吊杆和主梁,主缆的端点分别是主缆在左、右两个主索鞍上的切点,用铰接约束,所述主缆用索单元模拟,并输入各索段的初应变,所述主梁用梁单元模拟,其节点坐标依据桥面设计线形确定。
8.根据权利要求1所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,步骤(4)中,所述欧几里得范数检验与吊杆力假定值的误差公式表示为:
其中,Pi,k是用有限元法计算出的吊杆力;Pi,k-1是吊杆力假定值;ε是误差限值。
9.根据权利要求1所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,步骤(5)中,所述边跨主缆线形的计算方法为:
所述边跨主缆与主索鞍的切点为Q,与散索鞍的切点为M,设Q、M两点之间的水平距离为s,以Q点为原点建立坐标系,x轴水平向右,y轴竖直向下,方程表达为:
其中,c=-H/q,H为成桥状态主缆水平力,q为成桥状态主缆自重荷载集度;
根据主索鞍IP点与散索鞍IP点之间的水平距离和高差得出:
Δ23+s+Δ4=L
h1+h2+h3=h
其中,Δ2是主索鞍圆心C'与桥塔中心线之间的水平距离,为已知量;Δ3是主索鞍圆心C'与边跨主缆在主索鞍上的切点Q之间的水平距离;Δ4是主缆在散索鞍上的切点M与散索鞍IP点之间的水平距离;L是悬索桥边跨跨径,h1是主索鞍IP点与边跨主缆在主索鞍上的切点Q之间高差;h2是边跨主缆在主索鞍上的切点Q与在散索鞍上的切点M之间的高差;h3是切点M与散索鞍IP点之间的高差;h是主索鞍IP点与散索鞍IP点之间的高差。
10.根据权利要求9所述的悬索桥吊杆力和主缆线形联合计算方法,其特征在于,所述主索鞍表示为:
Δ3=R'sinβ
其中,R'是所述主索鞍半径;β是所述直线C'Q与竖直线的夹角;
所述主索鞍IP点与边跨主缆在主索鞍上的切点Q之间高差h1表示为:
h1=(Δ23)tanβ
所述边跨主缆,两端点高差h2表示为:
所述散索鞍M、N两点之间的水平距离Δ4和高差h3分别表达为:
其中,r是所述散索鞍半径,ω是散索鞍圆心P点与IP点连线与竖直线的夹角,是所述散索鞍圆心P点与切点M连线与竖直线的夹角,从而得到边跨主缆的线形和切点位置。
CN201810250670.7A 2018-03-26 2018-03-26 一种悬索桥吊杆力和主缆线形联合计算方法 Active CN108491635B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810250670.7A CN108491635B (zh) 2018-03-26 2018-03-26 一种悬索桥吊杆力和主缆线形联合计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810250670.7A CN108491635B (zh) 2018-03-26 2018-03-26 一种悬索桥吊杆力和主缆线形联合计算方法

Publications (2)

Publication Number Publication Date
CN108491635A true CN108491635A (zh) 2018-09-04
CN108491635B CN108491635B (zh) 2022-04-12

Family

ID=63337627

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810250670.7A Active CN108491635B (zh) 2018-03-26 2018-03-26 一种悬索桥吊杆力和主缆线形联合计算方法

Country Status (1)

Country Link
CN (1) CN108491635B (zh)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271733A (zh) * 2018-10-09 2019-01-25 东南大学 一种悬索桥空缆安装线形和索鞍初始安装位置的确定方法
CN109583120A (zh) * 2018-12-11 2019-04-05 佛山科学技术学院 一种用于大跨度悬索桥主缆内力识别的方法
CN109799011A (zh) * 2019-03-27 2019-05-24 东南大学 一种悬索桥吊杆力测定装置
CN109815588A (zh) * 2019-01-24 2019-05-28 东南大学 一种不等主跨的三塔悬索桥的主缆线形构建方法
CN110345890A (zh) * 2019-06-24 2019-10-18 中铁大桥勘测设计院集团有限公司 一种成桥状态散索鞍位置和锚跨线形的联合确定方法
CN110939066A (zh) * 2019-12-06 2020-03-31 中铁大桥勘测设计院集团有限公司 一种锚跨丝股张拉力的确定方法
CN111241616A (zh) * 2020-01-08 2020-06-05 东南大学 一种三主缆悬索桥恒载横桥向均匀分配方法
CN111259582A (zh) * 2020-01-12 2020-06-09 哈尔滨工业大学 一种悬索桥索鞍处主缆长度快速精确计算方法
CN111523172A (zh) * 2020-05-11 2020-08-11 重庆交通大学 一种空间异形索面悬索桥主缆成桥线形分析方法
CN111967085A (zh) * 2020-08-31 2020-11-20 中国石油天然气集团有限公司 大跨径悬索管道桥成桥风索线形计算方法
CN112012110A (zh) * 2020-08-31 2020-12-01 东南大学 一种使三主缆悬索桥恒载横桥向均匀分配的装置及方法
CN112035928A (zh) * 2020-08-31 2020-12-04 中国石油天然气集团有限公司 大跨径悬索管道桥成桥主索线形计算方法
CN112815840A (zh) * 2020-12-31 2021-05-18 中交第二航务工程局有限公司 一种基于机器视觉的悬索桥一般索股线形测量方法
CN112907505A (zh) * 2021-01-11 2021-06-04 南京工程学院 一种水下机器人电缆三维形状估计方法
CN113255029A (zh) * 2021-04-19 2021-08-13 东南大学 一种活载作用下悬索桥结构变形及内力的确定方法
CN113468632A (zh) * 2021-05-26 2021-10-01 东南大学 一种偏心活载作用下的悬索桥全桥响应的确定方法
CN114910031A (zh) * 2022-06-27 2022-08-16 重庆交通大学 一种悬索桥健康监测方法、系统、计算机设备和存储介质
CN115357965A (zh) * 2022-06-27 2022-11-18 中建三局第一建设工程有限责任公司 自锚式悬索桥及其成桥线形确定方法
CN118246129A (zh) * 2024-05-28 2024-06-25 四川省交通勘察设计研究院有限公司 大跨度悬索桥主缆成桥线形分析方法、系统、终端及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4612680A (en) * 1984-03-26 1986-09-23 Harumoto Iron Works Co., Ltd. Cover joint and armor for bridge cable
JP2008007998A (ja) * 2006-06-28 2008-01-17 Se Corp 斜張ケーブルの架設方法
CN104594179A (zh) * 2015-02-06 2015-05-06 中建三局集团有限公司 一种纵向双吊杆自锚式叠合梁悬索桥吊杆张力的控制方法
CN106777778A (zh) * 2017-01-05 2017-05-31 南京林业大学 一种单侧悬吊的曲梁悬索桥建模方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4612680A (en) * 1984-03-26 1986-09-23 Harumoto Iron Works Co., Ltd. Cover joint and armor for bridge cable
JP2008007998A (ja) * 2006-06-28 2008-01-17 Se Corp 斜張ケーブルの架設方法
CN104594179A (zh) * 2015-02-06 2015-05-06 中建三局集团有限公司 一种纵向双吊杆自锚式叠合梁悬索桥吊杆张力的控制方法
CN106777778A (zh) * 2017-01-05 2017-05-31 南京林业大学 一种单侧悬吊的曲梁悬索桥建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HONGYOU CAO 等: "FEForm-finding analysis of suspension bridges using an explicit Iterative approach", 《STRUCTURAL ENGINEERING AND MECHANICS》 *
吴章旭: "基于ANSYS的悬索桥主缆优化找形计算方法", 《桥隧工程》 *

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271733B (zh) * 2018-10-09 2023-05-16 东南大学 一种悬索桥空缆安装线形和索鞍初始安装位置的确定方法
CN109271733A (zh) * 2018-10-09 2019-01-25 东南大学 一种悬索桥空缆安装线形和索鞍初始安装位置的确定方法
CN109583120A (zh) * 2018-12-11 2019-04-05 佛山科学技术学院 一种用于大跨度悬索桥主缆内力识别的方法
CN109815588A (zh) * 2019-01-24 2019-05-28 东南大学 一种不等主跨的三塔悬索桥的主缆线形构建方法
CN109815588B (zh) * 2019-01-24 2022-02-11 东南大学 一种不等主跨的三塔悬索桥的主缆线形构建方法
CN109799011A (zh) * 2019-03-27 2019-05-24 东南大学 一种悬索桥吊杆力测定装置
CN110345890A (zh) * 2019-06-24 2019-10-18 中铁大桥勘测设计院集团有限公司 一种成桥状态散索鞍位置和锚跨线形的联合确定方法
CN110939066A (zh) * 2019-12-06 2020-03-31 中铁大桥勘测设计院集团有限公司 一种锚跨丝股张拉力的确定方法
CN111241616A (zh) * 2020-01-08 2020-06-05 东南大学 一种三主缆悬索桥恒载横桥向均匀分配方法
CN111241616B (zh) * 2020-01-08 2023-08-01 东南大学 一种三主缆悬索桥恒载横桥向均匀分配方法
CN111259582B (zh) * 2020-01-12 2020-08-28 哈尔滨工业大学 一种悬索桥索鞍处主缆长度快速精确计算方法
CN111259582A (zh) * 2020-01-12 2020-06-09 哈尔滨工业大学 一种悬索桥索鞍处主缆长度快速精确计算方法
CN111523172B (zh) * 2020-05-11 2022-10-04 重庆交通大学 一种空间异形索面悬索桥主缆成桥线形分析方法
CN111523172A (zh) * 2020-05-11 2020-08-11 重庆交通大学 一种空间异形索面悬索桥主缆成桥线形分析方法
CN112035928A (zh) * 2020-08-31 2020-12-04 中国石油天然气集团有限公司 大跨径悬索管道桥成桥主索线形计算方法
CN111967085B (zh) * 2020-08-31 2022-03-01 中国石油天然气集团有限公司 大跨径悬索管道桥成桥风索线形计算方法
CN112035928B (zh) * 2020-08-31 2022-04-15 中国石油天然气集团有限公司 大跨径悬索管道桥成桥主索线形计算方法
CN112012110A (zh) * 2020-08-31 2020-12-01 东南大学 一种使三主缆悬索桥恒载横桥向均匀分配的装置及方法
CN111967085A (zh) * 2020-08-31 2020-11-20 中国石油天然气集团有限公司 大跨径悬索管道桥成桥风索线形计算方法
CN112815840A (zh) * 2020-12-31 2021-05-18 中交第二航务工程局有限公司 一种基于机器视觉的悬索桥一般索股线形测量方法
CN112815840B (zh) * 2020-12-31 2022-03-29 中交第二航务工程局有限公司 一种基于机器视觉的悬索桥一般索股线形测量方法
CN112907505A (zh) * 2021-01-11 2021-06-04 南京工程学院 一种水下机器人电缆三维形状估计方法
CN112907505B (zh) * 2021-01-11 2024-03-26 南京工程学院 一种水下机器人电缆三维形状估计方法
CN113255029A (zh) * 2021-04-19 2021-08-13 东南大学 一种活载作用下悬索桥结构变形及内力的确定方法
CN113255029B (zh) * 2021-04-19 2023-12-12 东南大学 一种活载作用下悬索桥结构变形及内力的确定方法
CN113468632A (zh) * 2021-05-26 2021-10-01 东南大学 一种偏心活载作用下的悬索桥全桥响应的确定方法
CN113468632B (zh) * 2021-05-26 2022-12-09 东南大学 一种偏心活载作用下的悬索桥全桥响应的确定方法
CN114910031B (zh) * 2022-06-27 2023-07-04 重庆交通大学 一种悬索桥健康监测方法、系统、计算机设备和存储介质
CN115357965B (zh) * 2022-06-27 2023-10-27 中建三局第一建设工程有限责任公司 自锚式悬索桥及其成桥线形确定方法
CN115357965A (zh) * 2022-06-27 2022-11-18 中建三局第一建设工程有限责任公司 自锚式悬索桥及其成桥线形确定方法
CN114910031A (zh) * 2022-06-27 2022-08-16 重庆交通大学 一种悬索桥健康监测方法、系统、计算机设备和存储介质
CN118246129A (zh) * 2024-05-28 2024-06-25 四川省交通勘察设计研究院有限公司 大跨度悬索桥主缆成桥线形分析方法、系统、终端及介质

Also Published As

Publication number Publication date
CN108491635B (zh) 2022-04-12

Similar Documents

Publication Publication Date Title
CN108491635A (zh) 一种悬索桥吊杆力和主缆线形联合计算方法
CN107587429B (zh) 悬索桥索夹安装位置修正方法
Mesnil et al. Linear buckling of quadrangular and kagome gridshells: a comparative assessment
CN107988908A (zh) 钢管拱总拼装方法
Müller et al. Low Reynolds number drag and particle collision efficiency of a cylindrical fiber within a parallel array
Beghini et al. Rankine’s theorem for the design of cable structures
Kong et al. Elastic instability analysis for slender lattice-boom structures of crawler cranes
CN116698471A (zh) 一种飞行器舵面静强度试验方法
Shea et al. Ascent Aerodynamic Force and Moment Database Development for the Space Launch System
CN102779356A (zh) 曲面网格化系统及方法
CN110487576B (zh) 损伤状态倾角对称斜率的等截面梁损伤识别方法
Xiao et al. Determination of the reasonable state of suspension bridges with spatial cables
CN109977432A (zh) 塔架门洞疲劳安全系数计算方法和计算装置
CN113139218B (zh) 盾构管片外径边线绘制方法、装置和计算机设备
Lamb et al. Effects of packed bed properties on local concentration and temperature patterns
Oliveira et al. A local mesh free method with the singularity subtraction technique
Rozvany et al. A computer algorithm for deriving analytically and plotting optimal structural layout
Wilcock et al. Methodology for stability assessment of discretised shell structures during robotic assembly
Wu et al. An alternative approach to the structural motion analysis of wedge-beam offshore structures supporting a load
Rivero-Rodriguez et al. Pdes on deformable domains: Boundary arbitrary lagrangian-eulerian (bale) and deformable boundary perturbation (dbp) methods
Chen et al. Modified approaches for calculation of effective length factor of frames
Janus et al. Numerical modeling of flow phenomena in mine drift using laser scan results
Kounadis et al. Dynamic buckling of a 2-DOF imperfect system with symmetric imperfections
CN108388698A (zh) 一种改进的力密度平衡形态分析方法
WO2019233091A1 (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