CN111611644A - 施工过程索力最优的斜拉桥斜拉索张拉方法 - Google Patents
施工过程索力最优的斜拉桥斜拉索张拉方法 Download PDFInfo
- Publication number
- CN111611644A CN111611644A CN202010466510.3A CN202010466510A CN111611644A CN 111611644 A CN111611644 A CN 111611644A CN 202010466510 A CN202010466510 A CN 202010466510A CN 111611644 A CN111611644 A CN 111611644A
- Authority
- CN
- China
- Prior art keywords
- cable
- bridge
- force
- stay
- stay cable
- 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
Links
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Bridges Or Land Bridges (AREA)
Abstract
本发明公开了一种施工过程索力最优的斜拉桥斜拉索张拉方法,该方法只需要在施工节段一次张拉斜拉索,可实现过程中索力变化小,成桥后索力的均匀;同时提出的方法可以避免重复的迭代分析,由于考虑了过程中的收缩徐变、材料非线性等影响,优化的索力结果会更可靠,对于同类桥梁斜拉索的张拉和指导现场施工有重要价值和意义;本发明提出了新的斜拉桥索力优化方法,以斜拉桥为工程对象进行分析,控制了张拉过程中索力的变化量最小,同时使成桥后的斜拉索的索力值均匀;该方法适用于所有悬拼的斜拉桥桥面系施工,适用范围非常广泛,为张拉斜拉索的优化提供依据。
Description
技术领域
本发明涉及交通运输桥涵工程技术领域,具体涉及一种施工过程索力最优的斜拉桥斜拉 索张拉方法。
背景技术
斜拉桥是高次超静定结构,其主要承重结构为斜拉索和主塔,主塔的受力又是斜拉索拉 力决定的,因此斜拉索的受力好坏决定了斜拉桥最终的使用寿命和年限。历史上由于斜拉索 断裂而导致的事故屡见不鲜,造成的经济损失不可估量。研究者们一般从两个方面来分析和 解决这个难题,其一是从材料上,可以使用强度更高,抗疲劳更好的材料;其二是从结构受 力方面,优化斜拉索拉力,在保证桥面平顺的情况下,使得索力均匀。由于材料方面受到当 前技术水平的限制,常用的也就是1860MPa,1960MPa的高强度钢绞线。在斜拉索的索力优化 方面有很多方法,譬如最小弯曲能法、矩阵法、刚性支撑连续梁法等,这些方法各有优缺点, 但有个共同点:需要成桥后二次调索。众所周知,斜拉桥成桥后二次调索的工作量非常大, 其主要原因是成桥后结构形式已经固定,不再像悬拼节段那么好调整,由于调索效果相对较 差,很多时候需要反复的调整,浪费人力、物力和工期,由于不断的张拉和放张,锚头处容 易磨损,有潜在的安全风险。因此急需一种结果可靠、效果好的斜拉索索力优化方法。
发明内容
本发明的目的在于针对现有技术的不足,提出一种施工过程索力最优的斜拉桥斜拉索张 拉方法。该方法是通过对实际桥梁的设计复核后,建立正确的有限元模型,构造出相应的优 化体系,确定多目标优化方程和约束函数,最终优化出“过程索力变化最小,结果可靠”的 施工斜拉索索力,指导现场施工和拼装工作;本方法最大的优点就是过程中索力变化小,疲 劳损伤小,结果可靠,不需要二次调索,带来很大的经济价值。
为了实现上述目的,本发明采用了以下技术方案:
一种施工过程索力最优的斜拉桥斜拉索张拉方法,包括以下步骤:
步骤一,建立一个有限元模型,并且有限元模型中的材料特性、几何特性、边界条件与 外部荷载信息应和设计图纸的相应数据一致;
步骤二,确定一组斜拉索初拉力,并将确定的这一组斜拉索初拉力代入步骤一中的有限 元模型进行施工阶段正装分析,计算得到初始值向量u0、M0;
步骤三,建立一个优化体系,所述的优化体系要同时具备约束函数、优化方程和参数变 量,其优化公式为:
式中,u0表示斜拉索初拉力作用下成桥后的桥面控制测点位移值;Ti end表示成桥后第i 根斜拉索的索力值;T表示优化方程中待求解的索力增量值,为一个向量;Ti表示T中第i 根待求解的斜拉索的索力增量值;ui表示悬拼梁控制点成桥后的位移(i=1,2,3…n);M0表 示悬拼梁的位移增量矩阵;Δ1表示位移收敛允许值,Δ2表示斜拉索索力的不均匀度收敛允许 值,ε1和ε2是对应优化体系得出的最终位移收敛值和索力不均匀度收敛值;σlim表示斜拉索 的破断强度,单位MPa;dest表示控制点目标位移向量;Ai表示第i根斜拉索计算公称面积, 单位mm2;uij表示第j根斜拉索张拉单位力对第i节段悬拼梁位移增量值;
步骤四,假定一组斜拉索初拉力T0,在T0的作用下得到悬拼梁控制点的成桥后的位移ui (i=1,2,3…n);
步骤五,通过不断的变换索力T0+ΔT0得到相应的索力对成桥节段悬拼梁的位移增量矩阵 M0;
步骤六,引入收敛允许值Δ1和Δ2,均为约束条件,;
步骤七,设定一个控制点目标位移向量dest;
步骤八,将已经确定的各个参数代入优化公式计算得到最优的斜拉索拉力,进而应用于 斜拉桥斜拉索施工作业中,该优化后的索力是斜拉桥张拉控制指令的重要组成部分,是指导 施工的关键。
本发明提出在斜拉桥悬拼梁的施工过程中优化索力,把过程中最小的索力增量作为目标 函数,把控制成桥后索力均匀且悬拼梁的位移满足设计要求,作为约束条件进行分析。所述 的基于施工过程索力最优的斜拉桥斜拉索张拉方法适用于不同施工顺序和跨径的悬拼施工的 斜拉桥,适用范围非常广泛,优化结果显示成桥后不需要二次调索。具有非常大的实际推广 价值及广阔的工程运用前景,目前在多座斜拉桥施工控制中应用。
本发明进一步说明,所述的斜拉索初拉力,采用当前吊装重量的0.3~0.7倍作为初值。
本发明进一步说明,所述的控制点目标位移向量dest,在(-5,5)区间取值,单位为mm。
本发明进一步说明,各个矩阵维数要一致。
本发明还提供施工过程索力最优的斜拉桥斜拉索张拉优化体系,所述的优化体系要同时 具备约束函数、优化方程和参数变量,其优化公式为:
式中,u0表示斜拉索初拉力作用下成桥后的桥面控制测点位移值;Ti end表示成桥后第i 根斜拉索的索力值;T表示优化方程中待求解的索力增量值,为一个向量;Ti表示T中第i 根待求解的斜拉索的索力增量值;ui表示悬拼梁控制点成桥后的位移(i=1,2,3…n);M0表 示悬拼梁的位移增量矩阵;Δ1表示位移收敛允许值,Δ2表示斜拉索索力的不均匀度收敛允许 值,ε1和ε2是对应优化体系得出的最终位移收敛值和索力不均匀度收敛值;σlim表示斜拉索 的破断强度,单位MPa;dest表示控制点目标位移向量;Ai表示第i根斜拉索计算公称面积, 单位mm2;uij表示第j根斜拉索张拉单位力对第i节段悬拼梁位移增量值。
与现有技术相比较,本发明具备的有益效果:
1.本发明采用的方法优化后使斜拉索的索力值均匀,同时保证悬拼桥面梁的施工中斜拉 索的索力值变化最小,实现了过程中斜拉索的索力最优。
2.本发明考虑了整个现场的施工过程,避免了传统方法无法考虑桥梁的混凝土收缩、徐 变和结构几何非线性问题,在具体实施中简单实用;克服了悬拼过程中一次、二次张拉调索 和成桥后由于索力不均匀再次调索的繁琐施工,大大减少对以成型的桥梁结构的扰动;相对 于目前大部分学者研究的仅在成桥节段进行索力优化,更能指导施工和改善桥梁的受力,具 有重要意义。
3.本发明采用的优化方法和计算过程简单、易行,无论跨径、施工顺序如何改变都可以 进行优化分析,具有非常好的实际工程应用价值。
4.克服了成桥斜拉索索力不均匀,安全系数差异较大,部分斜拉索应力水平过高、降低 斜拉索疲劳寿命,解决了斜拉索使用寿命过短等工程难题。
附图说明
图1为某桥桥型设计立面图。
图2为某桥横梁构造图。
图3为某桥分析模型图。
图4为某桥成桥后斜拉索索力优化值和实测值对比图。
图5为本发明的斜拉桥斜拉索张拉方法流程图。
具体实施方式
下面结合附图对本发明进一步说明。
实施例:
某斜拉桥全长548m,桥跨组成为:46m+88m+280m+88m+46m。某斜拉桥设计为双塔、双索面预应力混凝土梁结构,主梁采用边主梁结构,斜拉索采用现场组装的镀锌钢绞线体系, 主塔采用弧线形门式框架钢筋混凝土结构,见附图1;“工”字形钢横梁根据横桥向宽度分为 I、II、III、IV共4种类型,钢横梁中心线处梁高2315mm,顶板600mm×24mm、底板600mm×32mm、腹板厚12mm、加劲肋厚10mm,材质为Q345qD,每片标准钢横梁重约13.3t,横 向布置图见附图2;某斜拉桥的斜拉索为空间双索面扇形布置,采用Φs15.2mm镀锌钢绞线 索,标准强度1860MPa,钢绞线采用PE护套单根防护后,再用外表面带双螺旋线的HDPE索 套管整体防护。全桥共有斜拉索136根,17种规格:Φs15.2-30、Φs15.2-31、Φs15.2-33、 Φs15.2-36、Φs15.2-40、Φs15.2-43、Φs15.2-46、Φs15.2-49、Φs15.2-52、Φs15.2-55、Φs15.2-57、 Φs15.2-59、Φs15.2-63、Φs15.2-66、Φs15.2-70、Φs15.2-76、Φs15.2-81。该工程采用本发 明提出的新方法进行优化分析。计算模型见图2所示。
计算跨径L=280000mm,Δ1=2mm,dest=0mm,Δ2=100kN,
T0=[150,300,600,900,1200,1500,1700,1950,2150,2500,2700,2850,3400,3600,4100,4300,4600]TkN,
u0=[5,13,21,30,37,45,56,66,73,86,93,103,112,125,136,144,151]TkN,
uii=[0.06,0.11,0.21,0.33,0.43,0.56,0.65,0.77,0.86,0.66,2.10,2.27,2.37,2.49,2.56,2.56,2.56]Tmm
通过优化分析得到最终斜拉索索力Ti end为:
[2662 2786 3077 3370 3660 3931 4195 4430 4653 4979 5176 5360 59006120 6572 6827,7066]T
成桥索力和常规方法对比曲线见图4所示。
可见本方法得到的成桥斜拉索索力变化均匀,钢绞线的应力值不均匀度明显小于常规方 法优化的斜拉桥索力,由于索力无突变,抗车辆冲击荷和风荷载能力更强,从而大大增加斜 拉索疲劳寿命,延长了斜拉索使用寿命。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技 术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不 会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (5)
1.一种施工过程索力最优的斜拉桥斜拉索张拉方法,其特征在于,包括以下步骤:
步骤一,建立一个有限元模型,并且有限元模型中的材料特性、几何特性、边界条件与外部荷载信息应和设计图纸相应的数据一致;
步骤二,确定一组斜拉索初拉力,并将这一组斜拉索初拉力代入步骤一中的有限元模型进行施工阶段正装分析,计算得到初始值向量u0、M0;
步骤三,建立一个优化体系,所述的优化体系要同时具备约束函数、优化方程和参数变量,其优化公式为:
式中,u0表示斜拉索初拉力作用下成桥后的桥面控制测点位移值;Ti end表示成桥后第i根斜拉索的索力值;T表示优化方程中待求解的索力增量值,为一个向量;Ti表示T中第i根待求解的斜拉索的索力增量值;ui表示悬拼梁控制点成桥后的位移(i=1,2,3…n);M0表示悬拼梁的位移增量矩阵;Δ1表示位移收敛允许值,Δ2表示斜拉索索力的不均匀度收敛允许值,ε1和ε2是对应优化体系得出的最终位移收敛值和索力不均匀度收敛值;σlim表示斜拉索的破断强度,单位MPa;dest表示控制点目标位移向量;Ai表示第i根斜拉索计算公称面积,单位mm2;uij表示第j根斜拉索张拉单位力对第i节段悬拼梁位移增量值;
步骤四,假定一组斜拉索初拉力T0,在T0的作用下得到悬拼梁控制点的成桥后的位移ui(i=1,2,3…n);
步骤五,通过不断的变换索力T0+ΔT0得到相应的索力对成桥节段悬拼梁的位移增量矩阵M0;
步骤六,引入收敛允许值Δ1和Δ2,均为约束条件;
步骤七,设定一个控制点目标位移向量dest;
步骤八,将已经确定的各个参数代入优化公式计算得到最优的斜拉索拉力,进而应用于斜拉桥斜拉索施工作业中,该优化后的索力是斜拉桥张拉控制指令的重要组成部分,是指导施工的关键。
2.根据权利要求1所述的施工过程索力最优的斜拉桥斜拉索张拉方法,其特征在于:所述的斜拉索初拉力,采用当前吊装重量的0.3~0.7倍作为初值。
3.根据权利要求1所述的施工过程索力最优的斜拉桥斜拉索张拉方法,其特征在于:所述的控制点目标位移向量dest,在(-5,5)区间取值,单位为mm。
5.根据权利要求1所述的施工过程索力最优的斜拉桥斜拉索张拉方法,其特征在于:各个矩阵维数要一致。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010466510.3A CN111611644B (zh) | 2020-05-28 | 2020-05-28 | 施工过程索力最优的斜拉桥斜拉索张拉方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010466510.3A CN111611644B (zh) | 2020-05-28 | 2020-05-28 | 施工过程索力最优的斜拉桥斜拉索张拉方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111611644A true CN111611644A (zh) | 2020-09-01 |
CN111611644B CN111611644B (zh) | 2023-03-21 |
Family
ID=72205590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010466510.3A Active CN111611644B (zh) | 2020-05-28 | 2020-05-28 | 施工过程索力最优的斜拉桥斜拉索张拉方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111611644B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112211114A (zh) * | 2020-10-20 | 2021-01-12 | 中国矿业大学 | 一种塔梁固结体系斜拉桥的悬臂拼装施工方法 |
CN112231805A (zh) * | 2020-09-28 | 2021-01-15 | 天津市市政工程设计研究院 | 一种基于大数据的组合梁斜拉桥索力估计方法 |
CN114741938A (zh) * | 2022-06-13 | 2022-07-12 | 中铁大桥科学研究院有限公司 | 一种斜拉桥调索非弹性收缩量智能监测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105568864A (zh) * | 2016-03-02 | 2016-05-11 | 东北林业大学 | 确定斜拉桥合理施工索力的综合算法 |
CN108446438A (zh) * | 2018-02-09 | 2018-08-24 | 广西交通科学研究院有限公司 | 刚构-拱组合桥成桥最优索力确定及快速实现方法 |
CN108460229A (zh) * | 2018-03-21 | 2018-08-28 | 武汉理工大学 | 连续桥面结构桥梁拉索内力的调整方法 |
CN110472306A (zh) * | 2019-07-26 | 2019-11-19 | 武汉工程大学 | 一种斜拉桥索力优化方法、装置、设备和可读存储介质 |
CN110700072A (zh) * | 2019-11-08 | 2020-01-17 | 中国矿业大学 | 一种斜拉桥缩尺试验模型斜拉索的安装方法 |
CN110765534A (zh) * | 2019-10-28 | 2020-02-07 | 江苏科技大学 | 一种斜拉桥成桥索力优化方法 |
CN110807221A (zh) * | 2019-11-04 | 2020-02-18 | 珠海交通工程技术有限公司 | 基于等效力位移法的索力计算方法 |
CN111008500A (zh) * | 2019-12-16 | 2020-04-14 | 郑州大学 | 一种斜拉桥的斜拉索初张拉力的计算方法 |
-
2020
- 2020-05-28 CN CN202010466510.3A patent/CN111611644B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105568864A (zh) * | 2016-03-02 | 2016-05-11 | 东北林业大学 | 确定斜拉桥合理施工索力的综合算法 |
CN108446438A (zh) * | 2018-02-09 | 2018-08-24 | 广西交通科学研究院有限公司 | 刚构-拱组合桥成桥最优索力确定及快速实现方法 |
CN108460229A (zh) * | 2018-03-21 | 2018-08-28 | 武汉理工大学 | 连续桥面结构桥梁拉索内力的调整方法 |
CN110472306A (zh) * | 2019-07-26 | 2019-11-19 | 武汉工程大学 | 一种斜拉桥索力优化方法、装置、设备和可读存储介质 |
CN110765534A (zh) * | 2019-10-28 | 2020-02-07 | 江苏科技大学 | 一种斜拉桥成桥索力优化方法 |
CN110807221A (zh) * | 2019-11-04 | 2020-02-18 | 珠海交通工程技术有限公司 | 基于等效力位移法的索力计算方法 |
CN110700072A (zh) * | 2019-11-08 | 2020-01-17 | 中国矿业大学 | 一种斜拉桥缩尺试验模型斜拉索的安装方法 |
CN111008500A (zh) * | 2019-12-16 | 2020-04-14 | 郑州大学 | 一种斜拉桥的斜拉索初张拉力的计算方法 |
Non-Patent Citations (4)
Title |
---|
戈士寅: "斜拉桥成桥和施工阶段索力优化应用研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
戴杰等: "斜拉桥成桥索力优化方法研究综述", 《中国公路学报》 * |
谢小辉: "多塔斜拉桥成桥索力优化及合理施工状态的确定", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
赵煜炜: "支架施工斜拉桥的施工控制与索力优化研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112231805A (zh) * | 2020-09-28 | 2021-01-15 | 天津市市政工程设计研究院 | 一种基于大数据的组合梁斜拉桥索力估计方法 |
CN112231805B (zh) * | 2020-09-28 | 2022-04-19 | 天津市市政工程设计研究院 | 一种基于大数据的组合梁斜拉桥索力估计方法 |
CN112211114A (zh) * | 2020-10-20 | 2021-01-12 | 中国矿业大学 | 一种塔梁固结体系斜拉桥的悬臂拼装施工方法 |
CN114741938A (zh) * | 2022-06-13 | 2022-07-12 | 中铁大桥科学研究院有限公司 | 一种斜拉桥调索非弹性收缩量智能监测方法 |
CN114741938B (zh) * | 2022-06-13 | 2022-08-19 | 中铁大桥科学研究院有限公司 | 一种斜拉桥调索非弹性收缩量智能监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111611644B (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111611644B (zh) | 施工过程索力最优的斜拉桥斜拉索张拉方法 | |
CN111581714B (zh) | 施工过程纠偏的拱桥吊杆拉力优化方法 | |
CN102864742B (zh) | 一种斜拉桥先滑后固型组合钢锚箱锚固装置及其施工方法 | |
CN213038218U (zh) | 一种预见性的可补张拉预应力的连续刚构桥 | |
CN103643637A (zh) | 大跨度钢管混凝土拱桥先竖转后原位整体提升施工方法 | |
CN111254803A (zh) | 一种三跨连续钢箱梁悬索桥的无索区线形控制方法 | |
CN111321846B (zh) | 一种装配式h形预应力撑杆钢柱 | |
CN114541776B (zh) | 一种拉索提升张拉转换装置及其施工方法 | |
CN109797716A (zh) | 一种防渗漏的渡槽及其施工方法 | |
CN109711041B (zh) | 温度自适应塔梁顺桥向约束方法及系统 | |
KR102066355B1 (ko) | 측면 좌굴을 최소화한 강합성 거더 제조방법 | |
CN110528401A (zh) | 悬浇施工挂篮系统 | |
CN214459551U (zh) | 一种节段间栓焊式组合连续梁 | |
CN111625894B (zh) | 位移最优的斜拉桥悬拼施工控制方法 | |
CN202830763U (zh) | 一种斜拉桥先滑后固型组合钢锚箱锚固装置 | |
CN111576230B (zh) | 施工扰动小的拱桥格子梁吊装位移控制方法 | |
CN112411353A (zh) | 一种纵向分段栓接的预制装配式钢混组合梁及施工方法 | |
CN112069584A (zh) | 一种铁路用延性构造墩的设计方法 | |
CN214033422U (zh) | 一种适用于铰缝改造的加固构造 | |
LU501216B1 (en) | Simply supported-continuous steel-concrete composite box beam partly shear-resistant and non-pull-resistant connecting piece structure | |
CN114673066B (zh) | 一种地锚式斜腿刚构桥及其施工方法 | |
CN219195625U (zh) | 一种斜拉索与主梁的锚固细部结构 | |
CN212270680U (zh) | 一种快速施工的钢-混组合梁混凝土桥面板及其吊装系统 | |
CN214459550U (zh) | 一种纵向分段栓接的预制装配式钢混组合梁 | |
CN110258287B (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 |