用节线法计算悬索桥缆形

刘政伟;高能祥

【摘 要】In this paper,the initial equilibrium state of suspension bridge was analyzed by nodal line method.First,primary parameters such as form of cable line,length theMATLABby

of

unstressed

section,etc.was results results

ofwere

calculated compared

as ANSYS

the

by to

using those

program,and the

MATLABprogram.Finally,the MATLABboundary

condition,the finite element model was built by using software.The

results show that nodal line method is not precise enough in analyzing the initial equilibrium state of suspension bridge.%采用节线法对悬索桥进行初始平衡状态分析，编制 MATLAB 计算程序，计算了实例的缆索线形及其无应力长度等主要参数；同时利用 MIDAS 对程序的计算结果进行了分析对比。最后以 MATLAB 计算结果为边界条件建立 ANSYS 有限元模型，进行结果分析，验证了节线法在进行悬索桥初始平衡状态分析时不很精确的结论。

【期刊名称】《交通科技》

【年(卷),期】2014(000)003

【总页数】3页(P26-28)

【关键词】悬索桥;节线法;初始平衡状态分析;MIDAS/Civil;Ansys

【作 者】刘政伟;高能祥

【作者单位】湖北省交通规划设计院 武汉 430051;华中科技大学土木工程与力学学院 武汉 430074

【正文语种】中 文

1 悬索桥的初始平衡状态分析

1.1 悬索桥初始平衡状态分析的概念

悬索桥的主缆承担桥梁自重和施工荷载，成桥后在恒载下主缆线型和吊杆内力应与设计目标相同［1］。施工过程中悬索桥的主缆和加劲梁的几何形状变化非常大，所以设计悬索桥时，要做逆施工阶段分析，进行倒拆分析需要做悬索桥在恒载下的初始平衡状态分析［2］。悬索桥在加劲梁的自重作用下产生变形后达到平衡状态，在满足设计要求的垂度和跨度条件下，计算主缆的坐标和张力的分析一般称为初始平衡状态分析［3］。这是对运营阶段进行分析的前提条件，所以应尽量使初始平衡状态分析结果与设计条件一致。

1.2 节线法进行悬索桥初始平衡状态分析

节线法采用了日本Ohtsuki博士［4］使用的计算索平衡状态方程式，是利用桥梁自重和主缆张力的平衡方程计算主缆坐标和主缆张力的方法。

一般做法，是将索分别投影到竖直面和水平面上，利用在各自平面上张力和恒荷载的平衡关系进行分析。下面分别介绍竖向和水平面内的分析过程。

1.2.1 竖向平面内的分析

假设一个跨度内的吊杆数量为N－1，则吊杆将该跨分割成N跨。

竖向上主缆形状及力的平衡图见图1。

图1 X-Z平面内力的平衡

如图，Wsi为加劲梁和吊杆恒载平均到主缆上的均布荷载；Wci为主缆自重；T i为节点i-1和i之间的主缆单元张力；I i为主缆单元长度；T x为主缆张力水平分量；T x全跨相同。

根据力的平衡条件，在X-Z平面上，在第i个节点位置的平衡方程式如下。

在Y-Z平面上的平衡方程如下。

式中：Pi为第i个吊杆的张力；h i为吊杆的长度。

由式（1）和（2）可以得到 N-1个方程。

上面公式中的未知数为z i（i＝1，2，…，N－1）和T x，共有N个未知数。作为追

加条件使用跨中的垂度f与跨中、两边吊杆的竖向坐标的关系公式：

1.2.2 水平平面内的分析

上式可得到N－1个方程，其中水平张力T x可由竖向平面内的分析获得，主缆两端的y轴坐标y 0，y N为已知值，所以共有N－1的未知数，y i（i＝1，2，3，…，N－1）可通过方程组计算。

2 利用MATLAB进行实例分析

2.1 自锚式悬索桥实例参数

跨径布置为230 m＋460 m＋230 m，主缆矢跨比为1／10，桥塔塔顶的侧向刚度k b＝k c＝8 500 k N／m，跨中及边跨最短吊杆长为2 m。中跨主缆吊杆为29根，边跨吊杆为15根，吊杆间距均为46／3 m。其他主缆及加劲梁截面材料特性及恒活载如图所注。

算例中悬索桥参数及立面布置见图2。

图2 悬索桥立面布置图（尺寸单位：m）

以下建模计算过程中，节点划分做如下说明：吊杆从左到右依次编号，左边跨为1-15号吊杆，中跨为16-44号吊杆，右边跨为45-59号吊杆；主缆段从左到右依次编号，左边跨主缆段为1-15号；中跨为16-45号；右边跨为46-60号。

2.2 MATLAB建模过程

采用节线法的MATLAB建模流程见图3。

2.3 MATLAB程序计算结果

主缆坐标及线形。MATLAB计算出的全桥缆索线形图见图4。

图3 MATLAB建模流程图

图4 MATLAB计算的主缆线形图注：以桥面为纵坐标零点，未考虑桥面纵坡；图中＊点为吊杆连接点。

将MATLAB计算出的全桥缆索线形坐标汇 总于表1。

表1 MATAB计算的主缆坐标表 m1 2 3 4 5 6 7 8 9 10纵坐标 2 2.20 2.82 3.84 5.27 7.10 9.35 12.01 15.07 18.5吊杆编号4 11 12 13 14 15 16 17 18 19 20纵坐标 22.43 26.72 31.42 36.54 42.06 48.00 42.06 36.5吊杆编号4 31.42 26.72

3 利用MIDAS对MATLAB计算结果的验证

利用MIDAS悬索桥建模助手建立实例悬索桥模型，MIDAS模型见图5。

图5 MIDAS悬索桥模型注：未考虑桥面纵坡，其中缆索和吊杆为仅受拉单元，主梁

和桥塔为梁单元。

MIDAS模型计算结果与MTALAB计算结果对比如下：（1）主缆坐标计算结果对比。见表2。

表2 MIDAS与MATLAB主缆坐标计算对比 m吊杆编号17 18 19 20 21 22 23程序结果／k N 41.9 36.396 31.303 26.618 22.343 18.475 1项目5.016 MIDAS结果／k N 42.063 36.537 31.423 26.719 22.426 18.543 15.069相对误差／%0.39 0.39 0.38 0.38 0.37 0.36 0.35

（2）主缆内拉力计算结果对比。见表3。

表3 MIDAS与MATLAB主缆拉力结果对比.44 53 327.05 53 549.08 MIDAS结果／k N 53 030.04 53 067.31 53 141.79 53 253.36 53 401.82 53 586.92 53 808.34相对误差／%1 2 3 4 5 6 7程序结果／k N 52 768.58 52 805.97 52 880.67 52 992.56 53 141项目 缆索段编号0.49 0.49 0.49 0.49 0.49 0.48 0.48

4 利用ANSYS对节线法的论证

由于MIDAS中同样采用节线法对悬索桥进行初始平衡状态分析，由表3可知，MIDAS的计算结果和编写的MATLAB程序结果吻合很好。那么节线法计算结果的精确度如何，尚待检验。

笔者利用MATLAB程序的结果参数建立了ANSYS模型（具体介绍从略）。ANSYS

计算结果显示，主梁跨中仍会发生0.134 m的纵向位移，由此说明节线法计算出的平衡条件下建立的ANSYS有限元模型并不精确平衡。

5 结语

节线法在分析悬索桥主缆的初始平衡状态缆形时，仅利用力的平衡条件建立缆索内力的平衡方程求解，概念清晰、计算量小。Midas利用节线法对悬索桥进行初始平衡状态分析，在此基础上再进行运营阶段的线性、非线性的精确分析，可以大大提高工作效率。

然而通过节线法计算的悬索桥初始参数建立ANSYS模型，计算结果发现，模型初始状态并非精确平衡，仍会产生位移和不平衡内力，所以利用节线法进行悬索桥初始平衡状态分析时，计算结果并不很精确，必须对悬索桥进行后续的精确分析来消除节线法所带来的误差。

参考文献

［1］ 沈 慧.悬索桥施工仿真计算及设计参数分析［D］长沙：.长沙理工大学，2006.

［2］ 王立峰，孙 勇，李曼曼，等.自锚式悬索桥初始平衡状态分析［J］.中外公路，2010，30（6）：97-99.

［3］ 王戒躁，钟继卫.大跨度悬索桥主缆线形主要参数的影响性分析［J］.桥梁建设，2005（3）：21-24.

［4］ 葛俊颖.桥梁工程软件midas Civil使用指南［M］.北京：人民交通出版社，2013.