桥梁结构分析理论的基本问题是桥梁体系、构件中力与变形的求解。传统理论主要涵盖了桥梁结构的空间分析理论、各类梁桥(包括钢桥、混凝土梁桥及组合梁桥)、拱桥及大跨度桥梁(缆索承重桥梁为主)的静力计算理论,桥梁施工过程的时变力学问题也是桥梁结构分析理论研究的一个方面[1]。近年来,随非线性理论的成熟及通用有限元软件的应用推广,桥梁精细化分析理论发展迅速,基于三维实体有限元模型的桥梁抗火、船撞、车撞、落石冲击、波浪冲击及桥梁侧翻等方面的研究丰富了桥梁结构理论的范畴[2]。
论文综述了传统桥梁结构分析理论的研究进展,包括箱梁结构的空间分析理论及缆索承重桥梁计算理论两个分支,只包含静力分析。每部分对热点研究问题进行重点阐述,包括波形钢腹板箱梁理论、传统箱梁空间分析理论的拓展、多塔悬索桥中塔效应、缆索承重桥梁的极限承载力及钢-UHPC桥面板,分析存在的不足之处,总结相关理论的工程应用情况,在此基础上提出后续研究的建议。
高性能装配式组合桥梁结构在桥梁工程学科工程开发前沿占据重要一席,波形钢腹板的使用对高性能装配式组合桥梁结构的发展创新具有重要的作用[3]。波形钢腹板箱梁是波形钢腹板应用的主要形式之一,具有材料组合方式、板件连接方式灵活及体外预应力技术协作方便等优点,目前已在滁河大桥、鄄城黄河大桥、桃花峪大桥引桥、头道河大桥等工程中得到应用。近年来,相关研究重点关注了几种新型箱梁(图 1)的基本弯、剪、扭及横向受力的分析方法及力学行为。
张紫辰等[4]与王根会等[5]对波形钢腹板组合箱梁(底板为钢板)的剪力滞进行了研究,采用能量变分法推导了该结构形式的剪力滞后控制方程,并采用试验模型与有限元模型进行分析对比。Chen[6]等则通过模型试验、数值模拟及理论推导研究了如图 1(a)所示的波形钢腹板箱梁的剪力滞行为,这种箱梁形式充分利用了钢管混凝土的高承载力优势和桁架的减重功能,是一种轻型波形钢腹板箱梁。李运生等[7]考虑了曲率影响、截面剪力滞效应、波形腹板剪切效应、扭转和畸变效应,采用能量变分法推导了波形钢腹板简支曲线结合梁在弯扭效应下的解析解,随后通过固定跨度和截面模型研究了曲线半径和圆心角对弯扭效应的影响,得到了跨中挠度、扭转角等畸变角随圆心角的变化规律。姜瑞娟等[8]提出了基于3个广义位移函数的波形钢腹板组合箱梁剪力滞效应微分方程和基于有限差分法的半解析半数值解法。李丽园等[9]基于能量变分原理推导了简支箱梁的剪切附加挠度计算方法,结果表明,该方法与ANSYS空间有限元模型计算结果及实测结果吻合较好。Zhou等[10]提出了一种用于评价波形钢腹板预应力组合箱梁纯扭性能的改进柔性扭转膜模型,该模型可以有效地考虑混凝土受拉强度、受拉软化行为及混凝土的抗剪行为,经过数值模型及试验数据的对比,它的合理性与准确性得到了相应的验证,如图 2所示。Shen等[11]同样建立了用于波纹钢腹板预应力组合箱梁扭转全过程分析的软化膜模型。赵品等[12]提出了一种预测单箱双室波形钢腹板箱梁横向内力及变形的刚架模型,对比实验结果、传统箱梁框架模型和公路桥规中的简支板与连续板模型的横向内力计算结果,发现刚架模型计算结果较规范结果偏小。王兆南等[13]等基于框架分析法建立了波形钢腹板箱梁横向内力解析计算公式,并分析了波形钢腹板厚度、波折角、波长及箱梁宽高比等参数对横向内力的影响规律。
对经典箱梁理论的延伸研究也有一些报道,可以归纳成两条研究主线:1)基于解析法的传统箱梁理论的延伸。如郭增伟等[14]通过比拟杆法建立了变截面悬臂箱梁的剪力滞求解方法;Li等[15]采用能量法建立了适用于任意宽度翼缘的箱梁剪力滞解析法;张玉元等[16]应用能量变分法建立了以附加挠度为未知量的控制微分方程及边界条件,并导出了简支箱梁和两跨连续箱梁剪力滞附加挠度和翘曲应力的解析解;2)将箱梁空间分析理论与有限元理论结合起来发展箱梁的高等梁单元模型。例如,夏桂云等[17]利用初参数法和传递矩阵,建立了薄壁箱梁约束扭转分析的有限元列式;Lezgy等[18]采用能量原理推导了可考虑剪力滞后及钢混连接界面滑移行为的25自由度模型梁单元模型;David等[19]提出了一种可以准确分析混凝土徐变、截面翘曲及剪力滞后的通用梁理论(GBT);朱力等[20]提出了曲线钢-混凝土组合箱形梁考虑翘曲、畸变及双向滑移的22自由度有限梁单元模型。
对波形钢腹板箱梁及普通箱梁理论的研究主要采用解析法(能量法、比拟杆法等)对材料处于弹性范围内的情况而开展,对影响箱梁空间受力行为的许多因素均有不同程度的考虑,但仍存在方法局限性大、不易应用推广等缺点。文献中建立的有限元模型以板壳/实体有限元模型为主,对箱梁空间模型的探讨及应用不足,如单梁、多梁、梁格模型等空间杆系模型的研究缺乏。这些模型虽在计算精度上存在一定的缺失,但从文献[17-20]来看,采用箱梁的高等梁理论(通过引入扭转、畸变、翘曲、剪力滞等自由度形成的多自由度梁模型)来研究箱梁的空间效应,箱梁空间杆系模型的计算精度可以得到有效的提高,在兼有高效计算效率的优势下,形成箱梁结构高效精细的计算方法具有较大潜力。但是,还需考虑复杂的全过程受力行为,使箱梁理论进一步对接桥梁结构体系弹塑性分析理论,拓宽分析理论的适用范围。
中国海洋桥梁工程的调研报告指出[21],多塔长联缆索承重桥梁是一种具有显著技术经济环保优势的桥型。多塔悬索桥在该类型桥梁中占有举足轻重的地位。中国已建成3座多塔悬索桥,包括马鞍山长江大桥(主跨1 080 m)、江苏泰州长江大桥(主跨1 080 m)、武汉鹦鹉洲长江大桥(主跨850 m)。与传统单主跨悬索桥相比,该桥型的两个主跨在受活载偏一侧的作用下,具有较为明显的“中塔效应”(图 3),具体体现为:采用柔性中间桥塔使桥梁的整体刚度降低,而采用刚性中间桥塔则存在主缆与鞍槽的滑移风险。多塔悬索桥中间桥塔的设计涉及桥塔稳定性、主缆滑移及中塔适宜纵向刚度取值等关键问题。
沈锐利、王秀兰及Zhang等学者对多塔悬索桥的中塔效应进行了研究。沈锐利等[22]构建了主跨500~1 500 m范围内的6座三塔悬索桥计算模型,通过改变中塔纵向刚度,对加劲梁竖向挠度及主缆抗滑安全系数进行计算分析,建议了不同竖向刚度要求下中塔纵向刚度的取值范围,如图 4所示。从图中可以看出,中塔刚度Kt的对数值与结构刚度基本呈线性关系。王秀兰等[23]以挠跨比和主缆抗滑安全系数为控制指标,提出了多塔悬索桥中塔刚度上、下限值的解析算法,给出了跨径在1 000~2 000 m范围内的中塔刚度取值范围。Zhang等[24]基于分段悬链线理论推导了不等跨多塔悬索桥中间桥塔刚度的解析算法,并对一座跨度组合为248 m+1 060 m+ 1 360 m+380 m的背景桥梁进行了分析,根据计算分析得到中间桥塔的刚度取值范围为6.615~16.422 MN/m。
由于桥塔形式多样、分析方法不统一等,现阶段对几座已建多塔悬索桥中塔刚度取值范围并没有完全达成共识,对“中塔效应”的认识存在一桥一议的事实。但自多塔悬索桥建设以来,随相关科学研究工作的不断深入,逐步认识到多塔悬索桥中间桥塔稳定性特征[25]、侧向力对主缆抗滑不可忽略的贡献[26]等,多塔悬索桥中间桥塔(包括主缆抗滑设计在内)的设计正逐步从粗犷、经验、保守的方式转变为精细、科学、合理的方式,相关规范条文也正在积极补充。
实际桥梁结构的稳定性问题以第二类稳定问题为主,需全面考虑若干非线性因素的影响,特别是几何与材料非线性的影响,属于极限承载力的范畴。桥梁结构应按弹塑性理论分析获得典型的荷载变形曲线及加载过程中桥梁结构的应力、变形等状态,从而充分把握其极限承载模式,明确桥梁设计的安全储备,以保障全桥的静力安全性。关于桥梁极限承载力的研究主要以在建桥梁为背景而展开,笔者主要论述大跨径缆索承重桥梁极限承载力的研究进展。对拱桥的相关研究情况可参考文献[27-30]。
吕梁等[31]以南京长江五桥(三塔斜拉桥)为背景,采用ANSYS软件建立全桥杆系有限元模型,分析了桥梁施工阶段稳定系数的变化规律和斜拉索、桥塔、主梁的极限承载力及组合桥塔外包钢壳对极限承载力的影响。Kim等[32]采用“成桥状态确定-倒拆形成施工阶段-极限承载力分析”的三步法分析了斜拉桥施工阶段的稳定性。Yu等[33]对世界上最长的钢桁架斜拉桥——鸭池河桥面板施工过程稳定性能进行了研究,研究了几何非线性、材料非线性、温度变化和静态横向风对结构稳定性能的影响。王春江等[34]采用ABAQUS软件建立全桥多尺度有限元模型(跨中段钢箱梁建立局部板壳精细模型)对重庆鹅公岩轨道专用桥(主跨600 m的自锚式钢箱梁悬索桥)的极限承载力进行了分析,研究了该桥的整体稳定和局部稳定问题及其相互影响规律,结果发现,箱梁中薄弱的加劲钢腹板会在全桥破坏前屈曲。为进一步完善自锚式悬索桥的稳定理论,白伦华等[35]用结构稳定性概念定性分析了自锚式悬索桥不会发生面内弹性稳定性的可能(图 5),进而通过在自锚式悬索桥挠度方程中引入位移干扰量,以幂级数干扰位移形式反证其不会出现面内弹性分岔失稳。同时,以重庆鹅公岩轨道专用桥为例,通过数值模型按弹性及弹塑性稳定理论计算分析了桥梁的荷载系数、破坏模式等。
从多塔斜拉桥组合桥塔外包钢壳效应[31]、基于施工过程分析的斜拉桥稳定性计算方法[32]、自锚式钢箱梁悬索桥柔弱腹板局部屈曲效应[34]等方面的研究来看,大跨径桥梁极限承载力的分析理论正在向精细化方向逐步深化,对影响全桥极限承载力的诸多复杂因素考虑得更加全面。然而,桥梁局部-整体破坏精细化分析方法仍有诸多问题亟待解决,如主梁采用钢箱梁形式的桥梁局部-整体相关屈曲问题不够精确,组合结构桥梁中构件不同材料对极限承载力的贡献比不清楚,极限承载力分析所采用的材料本构模型不统一等。下面对主梁采用钢箱梁形式的桥梁局部-整体相关屈曲问题的研究作简单分析。
对于主梁采用钢箱梁形式的大跨径缆索承重桥梁,上述研究按通用有限元程序的梁柱单元理论,没有考虑几何缺陷、残余应力及局部屈曲的影响,这也是当前钢结构桥梁极限承载力分析普遍存在的问题。如何全面计入影响钢桥稳定性的非线性因素,使钢桥的极限承载力分析更加精确是值得探索与研究的问题。文献[34]虽然采用节段板壳有限元来模拟桥梁中容易破坏,但并没有计入几何缺陷、残余应力的影响,也没有相应的试验结果验证,所建立的多尺度模型在自锚式钢箱梁悬索桥的极限承载力分析中仍需进一步校验其准确性及有效性。进一步回顾采用板壳有限元法研究全桥极限承载力的几篇文献报道[36-38],由于桥梁结构规模庞大,对几何缺陷及残余应力均没有完全模拟,而是采用了不同的简化方式进行处理(图 6),在几何缺陷的考虑上,主要采用了模态引入法,对残余应力的模拟则是通过在单元积分点处添加初应力来实现。但从几种处理残余应力的方式来看,基于板壳有限元模型的钢桥极限承载力分析方法并不成熟,方法之间的差异也缺乏参照对比。对钢桥的局部-整体相关屈曲问题并没有得到很好的解决。对于宽幅箱梁结构,剪力滞后及桥面板横向受力对极限承载力的影响也没有很好地考虑。除此之外,桥塔或主梁采用混凝土的区域,对混凝土的材料行为,如塑性损伤、开裂效应等也不在考虑之中。桥梁局部-整体破坏关系的分析方法可向全面计入钢板结构的局部屈曲、箱梁的空间效应及材料的细观力学行为的桥梁高等二阶分析方法方面发展。
正交异性钢桥面板在长期服役过程中暴露出了疲劳开裂和铺装层损坏两个难以解决的问题,严重影响了缆索承重桥梁结构的正常使用。近年来,学者们以新材料新型桥面结构体系为主要方法从根本上克服钢桥面系统的顽疾。钢-超高性能混凝土(UHPC)桥面板经过方案提出、试验验证、理论研究及综合工程经济效益对比的系列研究,在缆索承重桥梁结构中已取得了初步的应用[39]。以湖南大学邵旭东教授团队为代表对钢-UHPC桥面板的研究较为系统,在静力性能方面,对栓钉、螺栓、短钢筋等连接件的抗剪性能、UHPC板的抗弯性能等进行了深入的研究[40-45]。
目前,对钢-UHPC桥面板的研究主要集中在构造、连接、构件方面,基本力学性能的探究较多,对相关的体系模型计算研究还比较有限。如何在体系模型中充分考虑UHPC铺装方式、徐变效应及钢-UHPC桥面板剪力滞后等的影响,并探讨缆索承重桥梁结构性能,都是值得进一步深入研究的问题[46-47]。
箱梁与缆索承重桥梁分析理论是与工程实际联系十分紧密的课题,它们关注桥梁建设中最基本的施工与设计,与桥梁工程建设相生相伴。
1) 箱梁理论的应用:箱梁空间分析理论中最为突出的问题是顶底板的剪力滞后效应,通过能量法或比拟杆法重新建立波形钢腹板箱梁的剪力滞计算方法,考察了不同形式中剪力滞效应受构造参数影响的情况,如桁式底板波形钢腹板箱梁中底部桁架管件中混凝土的影响、波形钢腹板的褶皱效应等,这些研究结论为该类结构的优化设计提供了指导思路。
2) 多塔悬索桥分析理论的应用:目前,有关中塔适宜刚度理论、主缆与中塔塔顶鞍槽滑移理论及抗滑技术及人字形钢桥塔的稳定理论等方面的丰硕成果为中国几座已建及在建的多塔悬索桥提供了理论支撑与技术保障。
3) 桥梁极限承载力分析理论的应用:以非线性有限元理论为基础的桥梁极限承载力计算方法具有较多的实现形式,如实体有限元模型、多尺度有限元模型及全杆系有限元模型,以这些模型为载体对斜拉桥、自锚式悬索桥及拱桥进行分析,丰富了各类桥型极限承载力理论的内涵,对一些特殊问题的把握,如多塔斜拉桥中间组合桥塔的承载模式、自锚式悬索桥的面内弹性稳定问题等,不仅有效地消除了设计隐患,也有利于更加客观地认识桥梁结构的非线性物理特征。该理论的发展得益于新的桥梁结构形式的出现,也很好地直接指导了桥梁的建设(图 7)。
4) 钢-UHPC桥面系技术理论的应用:钢-UHPC技术已在中国17座桥梁得到了应用[39],涵盖多种桥梁结构形式,涉及自锚式悬索桥、地锚式钢桁梁悬索桥、拱桥及斜拉桥等缆索承重桥梁体系,为新型大跨径组合结构桥梁体系发展奠定了良好的实践基础。
新型桥梁结构体系,如波形钢腹板箱梁桥、多塔悬索桥、多塔斜拉桥、腹板易屈曲钢箱梁自锚式悬索桥、超高性能混凝土组合桥梁等,催化了桥梁结构分析理论——箱梁与大跨径桥梁理论的发展,波形钢腹板箱梁理论得到了进一步的完善,多塔悬索桥中塔效应问题在更广的参数范围内被讨论,经典桥梁极限承载力理论应用到更多的实际工程中。结合先前关于箱梁与缆索承重桥梁分析理论的讨论,对后续研究提出以下建议:
1) 针对箱梁,构件的竖向弯曲、扭转及弯扭耦合效应的研究还需进一步积极开展与完善,需要考虑材料非线性的影响,这些理论属于薄壁构件理论,可以与传统梁理论进行接轨,形成高等梁柱有限元单元,建立高效、精细的箱梁模型,用于构件压弯剪扭复杂力学行为的全过程分析。
2) 将箱梁的高等梁理论进一步融入到缆索承重桥梁的非线性分析理论中,同时,应充分重视如混凝土的开裂、钢混组合结构的界面行为及钢板的局部屈曲等力学行为,使桥梁体系杆系模型更加精细。
3) 缆索承重桥梁的精细化分析理论还需囊括局部关键区域的分析方法,如悬索桥锚跨区域索股的应力分析、自锚式悬索桥锚固段受力的精细化模拟方法、多塔悬索桥中塔主缆与索鞍相对滑移的精细化计算方法等。
4) 随着海洋桥梁工程事业的推进,非等跨三塔悬索桥、超三塔悬索桥等不同类型多塔悬索桥的理论研究需迫切开展。
5) 针对UHPC铺装方式、徐变效应及钢-UHPC桥面板剪力滞后等,建立并完善钢-UHPC桥面缆索承重桥梁的体系分析计算方法。