有限元程序对拉索的模拟计算方法探讨.doc

有限元程序对拉索的模拟计算方法探讨有限元程序对索单元的模拟计算方法探讨邓 科，李 翠，郝海龙，王灿东，郑徳江(中交公路规划设计院有限公司，北京100088)摘要：本文讨论了目前冇限元软件对索单元模拟计算的考虑垂度的Ernst公式刚度折减和考虑初应力效应的儿何刚度方法，并给出了精确描述 索单元力学行为的解析方法根据目前桥梁工程中主缆和斜拉索常用的力 学和尺寸参数进行计算，分析两种模拟计算方法与解析解之间的误差，最 后总结了悬索桥和斜拉桥计算时对主缆和斜拉索应采用的合理的模拟计 算方法关键词：索单元，垂度，儿何刚度Discussion on FEM simulating method of cableDENG Ke, LI Cui, HAO Hai-Iong, WANG Can・dong, ZHENG De-jiang(CCCC Highway Consultants CO., Ltd.)Abstract: This paper discussed the simulation of cable elements andanalyzed with finite element software. The simulation of cable elements concluded Ernst stiffness reduction method considering sag of cable, geometric stiffness method considering the effect of initial stress and an accurate method of description for the cable mechanical behavior was proposed. According to the calculation with mechanics and dimensional parameters commonly used in main cable and stayed-cable, errors between the simulati on method and the an alytic soluti on were analyzed. Fin ally the reas on able simulati on method for the main cable and stayed-cable whileanalyzing suspension bridge and cable-stayed bridge was proposed.Keywords: cable, sag, geometric stiffness重力作用下存在垂度；另一方面，结构建模1概述索是指仅有轴向拉压刚度而没有抗弯刚度的受力构件，在桥梁工程中大量应用这种构件，如悬索桥中的主缆、吊索，斜拉桥中的斜拉索。

在桥梁设计时，目前的冇限元程序都是将主缆、吊杆和斜拉索采用桁架单元模拟一方面，索单元由于其不存在抗弯刚度，在都是基于成桥状态模拟结 构构形，此吋索单元内部存在初应力，这使得索单元的非线性行为非常复 杂在斜拉桥计算时，常釆用Ernst公式对斜拉索刚度进行折减；而悬索 桥计算吋又采用冇限位移理论，考虑主缆和吊杆的初应力來提高缆索系统 的刚度同样是索单元，何时考虑刚度折减何时考虑刚度提高，这是一个 让桥梁设计工程师倍感困惑的问题木文通过理论分析和数值验证，给出这一问题的解答2索单元的模拟计算方法2.1考虑垂度的刚度折减索单元在自重影响下会产生一定的垂度，这一垂度的大小与索力冇关,而且垂度与索力呈非线性关系，推导如下图1考虑垂度的索单元假设索的重力集度为q,弹性模量为E,索原长SO,变形后端点之间 直线长度为S,索与水平面的夹角为ci,索内的轴力为T,索中部m的垂 度为f,取索的一半为隔离体，对索端部取矩，不难得到：Tf=l8qS2cosa(1)当垂度较小时，可认为索的形状为抛物线，那么索的形状长度为：2232P 二 S+8fq3S=S+lcosa24T2(2)又根据索的轴力为T,那么其形状长度为：P???T?1??EA??SO(3)由(2)和⑶,可以得到：SST q2S3cos2aO=AS=EAS 24T2⑷根据弹性模量的定义：EdodTSSOeff=de=0dASA=A dT=E⑸+qSAcoslal2SEOT近似认为S~SO,其引起的误差可忽略不计，将⑸整理后，得到EEeff=q2S2Acos21+aE⑹12T这就是考虑垂度效应后对拉索轴向刚度进行修正的的Ernst公式，从这个公式不难发现，无论索两端轴力为何值，其等效刚度Eeff总小于弹性刚度E,而且随着拉力的减小，其刚度减小更多。

2.2考虑初应力的刚度提高众所周知，构件受到轴向压力后，其侧向刚度减小，当轴向压力达到构件的压曲荷载时，英侧向刚度为零, 此时构件发生失稳而构件受到轴向拉力时，其侧向刚度提高，抵抗侧向 荷载的能力加强这个特性是由于结构的几何刚度改变所引起索的几何刚度矩阵可以简单推导如下TT■J图2考虑初拉力的索单元如图2所示，拉索两端受到拉力为T,若拉索两端产生横向位移A i 和Aj,在小位移情形，假设拉索轴力保持T不变，则对i点取矩可得：V?j??ij?IT(7)其中为I拉索单元长度又根据竖向受力平衡，不难知道Vi二Vj,由此可以得到拉索抵抗侧向荷 载的平衡方程为：??Vi???l?l???i??V??Tj?l其中T??l?l?I?□?便是拉索的几何刚度??矩阵不难看出，拉索初拉力的存在使得其产生横向位移需耍施加侧 向力，而且随着拉索拉力T的增加，使拉索两端产生横向位移所需的侧向 力相应加大，意味着拉索横向刚度的提高2.3拉索力学行为的解析解事实上，索的受力行为具冇高度的非线性，对其力学行为的真实描述 不能仅考虑垂度和初拉力，应结合索的基本属性，通过受力分析加以确定对索单元，引入以下假定：⑴忽略索的抗弯刚度；(2)材料满足胡克定律;(3)索的抗拉刚度变形前后不发生变化。

设索弹性模量为E,横截面积为A,延米自重为q,水平长度为I,竖向高度为h,无应力长度为SO,变形后形 状长度为P0,左端的水平力和竖向力分别为H、Vi,右端的水平力和竖向 力分别为H、Vj,其中由于自重作用方向竖直向下，因此索段处处水平力 分量相同，均为H■JHVi图3索单元解析解设索段上一点M,它与左端点0之间的索段的无应力长度为s,变形 后形状长度为P，在笛卡尔坐标系中的坐标为x、y,该处的轴力，水平分 力和竖向分力分别为T、H、Vo那么根据索段长度定义，可知?22??dx??dy??dp???????dp????空9)根据受力平衡：dxdp?HT 和 dyVVi?qsdp?T?T(10)根据胡克定律，索段的伸长dp?T?ds???EA?l??(11)将(9)和(10)联立，不难得到T??H?2??Vi?qs?2(12)(10)、(11)利12)联立,可知 dxdxds?dpdpds (13)?HHEA?H2?Vi?qs2分离变量并积分，可得l?HSOEA■H?q??ln?22???Vi?qS0?H??Vi?qS0???(14)?ln???Vi?H2?V2i??????类似的dyds?dydpdpds?Vi?psVi?ps(15)EA■?H?2??V2i?qs?分离变量并积分，可得 h? qS20?2ViS02EA (16) ?l?H2?V22q??i?qS0?H2?Vi???在1和h确定的前提下，H、Vi和SO三个变量只有一个是独立的, 当其中一个确定后，剩下两个变量可以通过(14)和(16)联立求得。

dl, dh, dH和dVi的相互关系推导如下dl??l?HdH??l?VdVii???IH2?H■• •??• •1??Vi?qS?H2??Vi?qSO?20?1H2??V2i?qS0???(17)?11??VH2?V2?H2H??l???H2?Vi?qS021• •?dVH2?V2ii???BldH?B2dVidh??h??HdH?h?VdVii■H?• • • •H2?V?qS2iO?l??2?dHH?V2?????SO?l?Vi?qS0??EAq??H2??Vi?qS0?2?Vi????H2?V2?dVi?B3dH?Bi?4dVi???2.4各种模拟方法的对比研究为了对索的三种模拟计算方法进行对比研究，我们取用工程实践中常采用的平行钢丝拉索或索股的力学参数进行计算研究索的弹性模量取200GPa；索单元水平投影长度I取10-500m；索单元与水平面的夹角a取10c ~80拉索上端固定，下端水平约束，竖向自由，计算其竖向 发生单位位移时所需力，即为其下端的竖向刚度，如图4所示HH图4索单元受力模式若采用桁架单元不考虑垂度效应，也不考虑初拉力提高刚度，那么拉 索下端的竖向刚度为：KEAcos?sin2?0?1(19)对于考虑垂度影响的Ernst公式计算拉索下端的竖向刚度：K?sin2?1?EeffAcosI(20)对于考虑初应力刚度提高的拉索模拟方法，下端的竖向刚度：23KEAcos?sin?2?l?Tcos?(21)对于索单元采用解析方法，根据(17) (18),可知：K3?dVidh?BlB(22)1B47B2B3图5~图8给出了当索内初应力为600MPa、400MPa> 200MPa,不同 水平投影长度及倾角的KO, KI, K2与K3的比值。

3K/iKI (m)图 5 竖向刚度比(o =600MPa/a =10> 20° )3K/iKI (m)图 6 竖向刚度比(o =600MPa,a =30> 40° )3K/iKI (m)图 7 竖向刚度比(o=600MPa,a=50> 60° )3K/iKI (m)图 8 竖向刚度比(a=600MPa,a=70> 80° )3K/iKI (m)图 9 竖向刚度比(o =400MPa/a =10> 20° ) 3K/iKI (m)图］0 竖向刚度比(o =400MPa,a =30> 40° )3K/iKI (m)S 11 竖向刚度比(o =400MPaza =50> 60° )3K/iKI (m)图］2 竖向刚度比(a =400MPa,a =70> 80° )3K/iKI (m)图 13 竖向刚度比(o =200MPa,a =10> 20° )3K/iKI (m)图 14 竖向刚度比(o =200MPa,a =30> 40° )3K/iKI (m)图］5 竖向刚度比(o =200MPa,a =50> 60° )3K/iKI (m)图 16 竖向刚度比(o =200MPaza =70> 80° )从计算结果不难看出，当索的水平投影长度I较短，与水平面的夹角 Q较小，索体内初应力较大时，若不考虑初应力对刚度的提高，而单纯采用桁。