锚拉板索梁锚固系统几何参数精确求解方法
2021-09-22李渊詹璐张铭
李渊,詹璐,张铭
(湖北省交通规划设计院股份有限公司,湖北 武汉 430051)
1 概述
在大跨径斜拉桥中,斜拉索作为联系桥塔和主梁的重要受力构件,其锚点坐标、无应力索长、安装角度等参数对斜拉桥的设计和施工至关重要。关于拉索无应力索长、倾角等参数的近似求解法、高精度近似求解法、精确求解法,许多文献已有充分研究,但都是基于拉索梁端及塔端锚点坐标已知的前提。对于如图1所示采用锚拉板索梁锚固系统的斜拉索,为使结构受力合理,需要保证锚拉板轴线方向与斜拉索梁端锚点受力切线方向一致,然而此二者相互影响,事先却仅能确定拉索塔端锚点、锚拉板轴线与主梁交点、拉索锚点距主梁高度d等参数,拉索梁端锚点坐标无从得知,更无法得到拉索无应力索长等参数。
图1 锚拉板索梁锚固系统
针对采用锚拉板索梁锚固系统且斜拉索塔端张拉力已知的斜拉桥,该文基于精确悬链线索形理论,采用改进的牛顿下山法,通过建立内外双重循环迭代模型,在确保锚拉板轴线方向与斜拉索梁端受力切线方向一致(Δ=0)的前提下,求解得到斜拉索梁端锚点坐标及无应力索长、塔端倾角、梁端倾角等参数的精确解,其结果有助于判断近似解的合理性及精度。
2 求解悬链线索形方程
假定拉索为理想柔性,且横截面保持不变,对于如图2所示的拉索,由微元体平衡条件可得:
(1)
式中:q为拉索的重力集度;H为索端的水平分力;y为拉索各点的竖向坐标值;x为拉索水平方向上各点的坐标。由式(1),可求得悬链线索形方程:
(2)
图2 拉索计算简图
拉索各点切线方向与水平方向的夹角可通过对悬链线索形方程求导得到:
(3)
则各点处的张力:
(4)
其中,坐标原点处即桥塔处的拉索张力:
(5)
该公式可进一步改写为:
(6)
针对塔端(即坐标原点处)拉索张拉力Tt已知的情况,对方程:
(7)
进行求解,得到满足f(b)=0时b的值,将该值代入式(2)即可求得拉索的悬链线方程。
式(7)是关于b的超越方程,其值一般需要通过迭代方法求解。通过塔端张拉力在拉索弦线水平方向的分力H确定初值b0,则:
(8)
对式(7)求导得到:
(9)
通过改进的牛顿下山法迭代求解方程f(b),具体流程如下:
(1)选取b0为初始值。
(2)取下山因子λ0=1。
(4)计算f(bn+1),并比较|f(bn+1)|与|f(bn)|的大小,分以下两种情况:
① 若|f(bn+1)|<|f(bn)|,则当|bn+1-bn|<μ2时(μ2为适当小的正数),取b≈bn+1,计算过程结束;当|bn+1-bn|≥μ2时,把bn+1作为新的bn值,并重复回到(3)。
3 求解拉索各关键参数
将求解所得的b值代入式(3)可分别得到拉索在塔、梁锚固点处的切线角,其中梁端拉索切线方向与水平方向的夹角为:
θl=y′l=sh(a-2b)
(10)
塔端拉索切线方向与水平方向的夹角为:
θ0=y′0=sha
(11)
拉索的悬链线索长、弹性伸长量及无应力索长可由下列公式求得:
悬链线索长:
(12)
拉索弹性伸长:
(13)
从而拉索无应力索长:
S0=S-ΔS
(14)
对于图3所示的拉索空间角度,可通过C、B两点(图1)坐标求得αH,根据空间角度之间的关系可得:
(15)
cosα=cosαH·cosγ0=cosβ0·cosαV
(16)
cosβ=cosβH·cosγ0
(17)
cosγ=cosα0·cosγW
(18)
则:
(19)
(20)
其中γ0=θl。对于拉索在塔端的锚固点,上述公式同样适用,只需注意塔端锚固处γ0=θ0。
图3 拉索空间角度示意图
4 求解锚拉板安装角度
对于主梁采用锚拉板索梁锚固形式的情况,一般情况下,锚拉板上拉索锚固点相对主梁的高度d是固定值,但拉索所在平面内,锚拉板与水平方向的夹角并不明确,要使锚拉板安装角度与拉索在主梁锚点切线角度一致,可通过迭代求解,具体做法如下:
(1)取锚拉板安装角度初值η0为锚拉板轴线与主梁交点C和拉索在桥塔锚点B连线与水平面的夹角,该值也是合理范围内锚拉板安装角度所能取得的最大值。
(2)通过ηn、αH、d,即可求得拉索在锚拉板上的锚点(xn,yn,zn),结合拉索在塔端锚点(xt,yt,zt),可求得拉索在竖直和水平方向投影的初值hn和ln,将该值代入式(7)~(9),采用改进的牛顿下山法迭代求解拉索的悬链线索形方程,进而求得拉索在锚拉板处锚点的切线角度θln,此时求得的θln即是合理范围内拉索梁端最小的切线角。
(4)用步骤(3)中求得的η求解拉索在锚拉板上的锚点(xl,yl,zl),并基于该锚点求解拉索悬链线索形方程,进而求解拉索各相关参数。
计算流程如图4所示。
图4 计算流程图
5 工程算例
以湖北香溪河大桥为例,运用该文提出的方法对该斜拉桥中跨组合梁斜拉索EC01和EC17进行计算分析。香溪河大桥为主跨470 m的双塔组合梁斜拉桥,跨径布置为(48+48+78+470+78+48+48)m,其中,边跨100.75 m范围内为Π形混凝土边主梁,其余为Π形工字钢组合梁,索塔横桥向为钻石形,全桥采用17对空间斜拉索,拉索和组合梁之间采用锚拉板式索梁锚固形式,斜拉索与主塔之间采用钢锚梁锚固形式,桥型布置如图5所示。
坐标系统以顺桥向中心线为x轴(里程增加方向为正向),以竖直方向为z轴(向上为正向),桥塔横桥向中心线为y轴(按右手螺旋法则确定正向)。拉索计算所用参数如表1所示,其中拉索弹性模量取1.95×105MPa,采用该文计算方法,得到斜拉索各参数如表2、3所示。求解主梁-锚拉板交点坐标和拉索梁端锚点坐标所连直线的空间角度,可知该锚拉板轴线倾角与拉索的梁端倾角是一致的,这一点在实桥中也得以印证。
工程实践并不严格要求采用精确解,但如果不知道精确解,近似解的合理性及精度也就无从得知。实际应用中存在直接连接拉索的塔端锚点B、锚拉板轴线与主梁的交点C(图1)得到直线,并基于该直线与xy平面的夹角计算锚拉板轴线角度,进而求解拉索梁端锚点坐标和其他参数的做法。采用该近似法求解得拉索梁端相关参数如表4所示。对比两种计算方法,得到拉索梁端关键参数计算差值如表5所示。
图5 香溪河大桥桥型布置图 (单位:cm)
表1 拉索参数
表2 拉索空间角度
表3 拉索无应力索长及锚点坐标
表4 近似法计算结果
表5 计算差值对比
由表5可知:对于短索,估算方法与精确解相对接近,但对于长索,估算方法误差较大,难以满足工程要求。
6 结论
(1)采用改进的牛顿下山法,建立内外双重循环迭代模型求解悬链线索形方程,得到锚拉板索梁锚固系统下斜拉索梁端锚点坐标及其他几何参数的精确解。
(2)通过对比发现,拉索较长的情况下,采用直线假定得到的近似解与精确解的差距较大。因此对于采用锚拉板索梁锚固系统的大跨径斜拉桥,采用精确求解方法是非常必要的。