拱桥钢绞线吊杆不对称断丝后的力学行为

2019-05-09刘剑锋李元兵张启伟

同济大学学报（自然科学版） 2019年4期

刘剑锋，李元兵，张启伟

(同济大学土木工程学院,上海 200092)

拱桥短吊杆断丝的损伤力学行为控制着吊杆结构的整体力学行为.钢丝在腐蚀和疲劳交变作用下一旦发生脆断，吊杆截面内其余钢丝会在短期内进行内力重分布而进一步开裂脆断，直至吊杆脆性断裂并发生桥面垮塌事故(见图1)，如四川宜宾南门大桥(2001)、美国俄亥俄州Sliver桥(1967)、新疆孔雀河大桥(2011)、福建公馆大桥(2011)等.在这些桥梁的坍塌事故中，不仅反映出中、下承式拱桥所面临的严峻的结构安全问题，还折射出研究人员对拱、吊杆、梁组合体系中，复杂的有别于失稳或振动破坏的损伤断裂行为的认识问题.目前，桥梁工作者尤其学术界对吊杆的安全性问题越来越重视，并已开展了一些研究工作.文献[1]中通过对南门大桥短吊杆断裂事故进行研究后认为,温度、腐蚀和疲劳是导致短吊杆断裂的主要因素.徐俊[2]认为,护套破损老化、钢丝锈蚀和腐蚀疲劳开裂是斜拉索损伤退化的主要机理.Raoof等[3]发现,在钢绞线钢丝绳中，外层单根钢丝断裂后，在影响长度内随着到钢绞线断裂处距离的增加，脆断钢丝内力逐步增大.MacDougall[4]进一步指出,随着到钢绞线断裂处距离的增加，断丝内力呈指数增加，而未断钢丝则呈指数衰减.基于钢绞线吊杆自身构造的特殊性以及内部钢丝接触、挤压和滑移等因素的复杂性，刘剑锋等[5]详细研究了钢丝脆断后钢绞线吊杆对称断丝力学行为.本文主要针对不对称断丝后钢绞线吊杆损伤力学行为展开研究，着重讨论钢丝不对称脆断后钢绞线吊杆损伤力学模型及数值验证方法.

图1 短吊杆断裂事故

1 不对称断丝后的吊杆损伤力学模型

1.1 基本假定

吊杆截面内某一束钢绞线出现不对称断丝时，由于外部接触力不平衡，会产生垂直于钢绞线轴向的变形.这种变形产生以下几个影响：①增加与断丝相邻两根外层钢丝的螺旋半径；②减小其余三根外层钢丝的螺旋半径；③由于钢绞线不对称断丝，使得与断丝相邻两根外层钢丝的外部接触力增大，其余三根外层钢丝的外部接触力减小，因而不对称断丝后钢绞线内各外层钢丝内力分布不均匀.

理论推导前作如下基本假定：①所有断丝均出现在同一位置(沿钢束)；②所有未断钢丝的应变均在弹性范围内；③钢丝与钢丝间完全接触，且不考虑接触变形；④吊杆截面内摩擦因数均相等，为一定值；⑤不考虑吊杆横截面位移；⑥钢丝开裂脆断期间，吊杆截面钢丝无其余损伤.

1.2 基本平衡方程

吊杆截面内某一束钢绞线(7丝)如图2所示，断丝形式分别为1，1、4，1、4、6，1、4、7等四种.钢绞线内各钢丝半径为R，螺旋角为α，螺旋半径Rh=2R，捻距p=2πRhtanα，曲率半径ρ0=Rh/cos2α，钢丝弹性模量为E，钢丝净面积为A0，吊杆长度为L.以钢丝断裂处为参考点(坐标原点)，x轴平行于核心钢丝轴向，协同坐标s轴沿钢丝螺旋方向，两坐标轴之间的换算关系为s=x/sinα.钢丝与护套间的摩擦力较小，忽略不计.断丝后护套对断裂钢丝新的握裹力单位长度为wsb.各钢丝间摩擦因数为μ，断丝与相邻钢丝间接触力为Xc(s)，钢丝间单位长度的摩擦力f(s)=μXc(s)，断丝轴力为Nb(s)，断丝轴向应变为εb(s)，未断钢丝轴力为Nu(s).断丝前，吊杆拉力为F，核心钢丝净面积为Ac，核心钢丝轴力Nc=F/(1+6sin3α)，外层各钢丝轴力N0=Fsin2α/(1+6sin3α)[4].断n根钢丝后，吊杆拉力变为Fbn，未断裂钢丝拉力总和为Fu(x)，核心钢丝轴向应变εc(x)=Δu/Δx，其中Δu为钢绞线微段沿轴向Δx的变形.

图2 钢绞线吊杆截面内不对称断丝示意图

Fig.2 Diagram of asymmetric broken wire arrangement in cross-section of strand suspender

图3为距离断裂截面s处断丝微段Δs的静力平衡立面图，断丝微段Δs力平衡方程如下所示[4]：

图3 断丝静力平衡图

x方向上：

(1)

y方向上：

(2)

图4为断一根丝时钢绞线内核心钢丝接触力平衡图.与断丝相邻两根钢丝和核心钢丝间接触力分别为X+60°和X-60°，与断丝相对的钢丝与核心钢丝间接触力为X+180°，其余两根钢丝与核心钢丝间接触力分别为X+120°和X-120°，断丝与核心钢丝间接触力为Xb.由力平衡知：

X+60°+X-60°+Xb=X+120°+X-120°+X+180°

(3)

图4 断一根丝时核心钢丝接触力平衡图

1.3 不对称断丝损伤力学模型

因Δφ/2较小，故cos(Δφ/2)≈1，sin(Δφ/2)≈Δφ/2.联立方程(1)和(2)并求极限,可得断裂钢丝在x断面的轴力、应变,如下所示:

(4)

(5)

钢绞线轴向变形δp引起的螺旋半径变化δh分别为[6]

δh(+180°)=-δp

(6)

由钢丝螺旋半径变化δh所引起的未断裂钢丝应变分别为

(7)

钢绞线内部接触力(单位长度)可按X0=N0cos2α/Rh计算.外层钢丝与核心钢丝间的接触力分别为

(8)

各外层钢丝的轴力分别为

N±60°(x)=EA0ε±60°(x),N±120°(x)=EA0ε±120°(x),

N+180°(x)=EA0ε+180°(x)

(9)

将式(6)～(9)代入式(3)并消除δp，不对称断丝后各外层钢丝的应变分别为

(10)

由于吊杆剩余拉力

则根据式(5)、(10)可得

(11)

根据式(11)可得核心钢丝轴力,如下所示:

(12)

根据式(10)、(12)可得外层未断钢丝轴力分别为

N±60°(x)=1.25sin2αNc(x)-0.25Nb(x),

N±120°(x)=0.75sin2αNc(x)+0.25Nb(x),

N+180°(x)=0.5sin2αNc(x)+0.5Nb(x)

(13)

1.4 影响长度

影响长度La是指外层断丝与未断钢丝拉力不均匀区域(当x=±La/2时，外层各钢丝拉力相等)，可得

(14)

将式(14)代入式(4)可得不对称断丝后断丝影响长度La的表达式，如下所示：

(15)

1.5 钢铰线拉力剩余率

钢绞线拉力剩余率γ可定义为

(16)

在0≤x≤La/2时，钢绞线拉力与伸长量的关系式可表示为

(17)

由于吊杆两端外部边界约束条件不变，整根吊杆断丝前后总的伸长量不变[4]，即Δ=Δbn(Δ为吊杆断丝前总的伸长量，Δbn为吊杆断丝后总的伸长量).

当影响长度小于等于钢绞线长度，即La≤L时，钢铰线拉力剩余率

(18)

式中：

λb1=1-

当影响长度大于钢绞线长度，即La>L时，钢铰线拉力剩余率

(19)

2 数值分析与讨论

2.1 有限元模型

采用Ansys软件中单元生死功能来模拟钢丝的脆断问题[7-9]，钢丝采用solid45三维实体单元，屈服准则采用Von Mises屈服准则.根据该屈服准则，单轴拉伸状态下钢丝临界剪应力取σy/3(σy为屈服应力).完好、损伤钢丝弹性模量分别为1.996×105、1.932×105MPa，泊松比为0.27.损伤钢丝实体模型不考虑大应变属性，材料非线性选用多线性同向强化应力-应变关系，考虑钢丝间摩擦力和护套握裹力，计算参数如表1所示.有限元模型如图5所示.为降低计算成本，模拟时钢丝长度仅取350 mm.

表1 钢绞线吊杆计算参数

2.2 吊杆截面内力重分布规律

吊杆截面内力重分布的理论计算值与有限元模拟结果如图6所示.从图6可得到以下结果：

(1)断钢丝1时，在影响长度内，断裂处破断钢丝拉力为零，随着到断裂处距离的增加呈指数增加；与钢丝1相邻的外层未断钢丝拉力最大，随距离的增加呈指数衰减；与钢丝1相对的未断钢丝拉力最小，随距离的增加呈指数增加；在影响长度外，各钢丝拉力分布又变为均匀.

(2)断钢丝1、4时，与断丝相邻的钢丝2、5、6拉力增大且应变相等，并随距离的增加呈指数衰减；与断丝相对的钢丝3、7拉力减小，并随距离的增加呈指数增加.

图5 钢绞线吊杆有限元模型

(3)断钢丝1、4、6时，与断丝相邻的钢丝2、5、7拉力增大(钢丝2最显著)，并随距离的增加呈指数衰减；钢丝3拉力减小，并随距离的增加呈指数增加.

(4)断钢丝1、4、7时，非断丝拉力均增大，并随距离的增加呈指数衰减.由于有限元模拟结果相对理论推导而言有许多假定，并引入了一定的模拟误差，如钢丝间的接触变形、摩擦、滑移、挤压等均无法精确模拟，导致有限元模拟结果存在一定程度的失真，但对断丝瞬间吊杆截面内各钢丝内力重分布规律而言，两种结果近似并具有可比性，后续参数分析时重点依据理论推导公式进行分析.

a 理论值(断钢丝1,La>L)

b 理论值(断钢丝1,La≤L)

c 有限元模拟(断钢丝1)

d 有限元模拟(不对称断钢丝1、4)

e 有限元模拟(不对称断钢丝1、4、6)

f 有限元模拟(不对称断钢丝1、4、7)

图6 不对称断丝时钢绞线吊杆内力重分布规律

Fig.6 Stress redistribution of strand suspender with asymmetric broken wires

2.3 钢绞线拉力损失率

不对称断丝情况下，钢绞线拉力损失率理论计算与有限元模拟结果如表2所示.结果表明：①在影响长度范围内，钢绞线拉力损失率最低限值均高于钢绞线内钢丝断丝率，附加应力损失主要由护套应力增加或钢丝变形引起的；②钢绞线拉力损失率主要受护套握裹力、捻角、钢丝半径、摩擦因数、钢丝长度等影响；③钢丝长度越大、直径越小，摩擦因数越大、护套握裹力越大以及螺旋角越小，钢铰线拉力损失率就越小.钢丝拉力损失率参数敏感性分析如图7所示.

表2 钢绞线拉力损失率

2.4 影响长度

不对称断丝情况下，影响长度受摩擦因数、螺旋半径、螺旋角及护套握裹力等影响，如图8所示.由图8可知：影响长度与摩擦因数成反比，随摩擦因数增大，影响长度急剧降低；影响长度与螺旋半径成正比，随螺旋半径增大，影响长度基本呈线性增加；影响长度与螺旋角成正比，随螺旋角增大，影响长度基本呈线性增加；影响长度与护套握裹力成反比，随护套握裹力增大，影响长度呈指数衰减.

a 吊杆长度

b 钢丝直径

c 摩擦因数

d 护套握裹力

图7 钢铰线拉力损失率参数敏感性分析(不对称断钢丝1)

Fig.7 Parameter sensitivity analysis of strand’s remaining prestressing fraction with wire 1 broken

2.5 护套握裹力对吊杆内力重分布的影响

图9为护套握裹力分别为5、10、20、30 MPa时，不对称断钢丝1情况下各钢丝拉力重分布的演化规律.由图9可知：护套握裹力增加对钢丝1(断丝)的影响最为显著，对钢丝4和6(±60°)的影响次之，对钢丝5和7(±120°)的影响最不显著，对钢丝2(核心钢丝)和钢丝3(+180°)的影响居中.随护套握裹力的增加，钢丝1、3、5、7拉力逐步增大，而其余钢丝拉力逐步变小.这表明护套握裹力对断丝及与断丝相邻的钢丝内力重分布影响非常明显，对其余钢丝内力重分布的影响则小得多.

a 摩擦因数

b 螺旋半径

c 螺旋角

d 护套握裹力

图8 影响长度参数敏感性分析(不对称断1、4、6丝)

Fig.8 Parameter sensitivity analysis of affected length with wires 1,4 and 6 broken

a 钢丝1

b 钢丝2

c 钢丝3

d 钢丝4和6

e 钢丝5和7

图9 护套握裹力对吊杆内力重分布的影响

Fig.9 Effect of contact forces between sheathing and outer wires on stress redistribution of suspender

2.6 钢丝间摩擦力对吊杆内力重分布的影响

图10为摩擦因数分别为0.10、0.15、0.20、0.25时，不对称断钢丝1情况下各钢丝内力重分布的演化规律.由图10可知：摩擦因数增加对钢丝1(断丝)的影响最为显著，对钢丝4和6(±60°)的影响次之，对钢丝5和7(±120°)的影响最不显著，对钢丝2(核心钢丝)和钢丝3(+180°)的影响居中.随摩擦因数的增加，钢丝1、3、5、7拉力逐步增大，而其余钢丝拉力逐步变小.这表明摩擦因数对断丝及与断丝相邻钢丝的内力重分布影响非常明显，对其余钢丝内力重分布的影响则小得多.