斜拉索受力状态及几何参数的快速迭代分析方法_白冰.pdf
《斜拉索受力状态及几何参数的快速迭代分析方法_白冰.pdf》由会员分享,可在线阅读,更多相关《斜拉索受力状态及几何参数的快速迭代分析方法_白冰.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 40 卷第 5 期2023 年 5 月公路交通科技Journal of Highway and Transportation esearch and DevelopmentVol.40No.5May 2023收稿日期:20220408基金项目:中央级公益性科研院所基本科研业务费专项资金项目(20219012a)作者简介:白冰(1989),男,山东菏泽人,博士,副研究员.(bbai89 )*通讯作者:杨飞(1985),男,陕西岐山人,硕士,副研究员.(yangfei202 )doi:10.3969/j.issn.10020268.2023.05.008斜拉索受力状态及几何参数的快速迭代分析方
2、法白冰1,杨飞*1,2(1.交通运输部公路科学研究院,北京100088;2.中路高科交通检测检验认证有限公司,北京100088)摘要:斜拉索是大跨度斜拉桥的关键受力构件之一,其受力状态及几何参数的准确确定是开展结构设计分析和施工控制的重要基础和前提之一。针对这一问题,基于悬索解析理论以及数值求解技术,研究了斜拉索在给定目标索力条件下受力状态及几何参数等信息的快速高精度计算方法。首先,根据悬索解析分析,对适于工程应用的斜拉索线形、倾角及无应力索长等计算控制方程进行了分析总结,重点针对控制方程参数,在分析简单迭代格式的基础上,结合弦截法数值求解技术,提出了求解方程参数的快速迭代解法,并就方法性能进
3、行了对比验证。最后,应用上述方法对实际斜拉桥的斜拉索进行了分析探讨。结果表明:所提出方法具有较高的计算精度和良好的适用性,能够应用于实际桥梁工程斜拉索的控制之中;算法计算稳定、快速,首次迭代计算值即与最终收敛结果较为接近,迭代 35 次后计算即可准确收敛;随着斜拉索索长、竖向夹角不断增大,其非线性效应不断增强。算法针对各斜拉索仍能保持良好的精度和性能,取得了较好的应用效果,这为结构的高精度分析和斜拉索的精准制造、精确安装提供了参考。关键词:桥梁工程;受力状态;数值迭代;斜拉索;悬链线理论中图分类号:U448.27文献标识码:A文章编号:10020268(2023)05005807A Metho
4、d for Fast Iteration Analysing Mechanical State and GeometricParameters of Stay CableBAI Bing1,YANG Fei*1,2(1 esearch Institute of Highway,Ministry of Transport,Beijing 100088,China;2 China-oad Transportation Verification Inspection Hi-Tech Co.,Ltd.,Beijing 100088,China)Abstract:Stay cables are one
5、of the key load-bearing components of long-span cable-stayed bridge,theaccurate determination of its stress state and geometric parameters is one of the important foundations andprerequisites for conducting structural design analysis and construction control.In review of this issue,a fastand high-pr
6、ecision calculation method for the stress state and geometric parameters of stay cables under targetforce condition is studied based on suspension cable analysis theory and numerical solution techniques.First,based on the analysis of suspension cables,the calculation control equations for the shape,
7、inclination angleand unstressed length of the stay cables suitable for engineering applications are analyzed and summarized.Focusing on the parameters of the control equations,a fast iterative solution for solving the equationparameters is proposed based on the analysis of fix point iterative format
8、 and the numerical solution techniqueof the secant method,and the performance of the proposed method is compared and verified.Then,the staycables in an actual cable-stayed bridge are analyzed and discussed by using this method.The result shows第 5 期白冰,等:斜拉索受力状态及几何参数的快速迭代分析方法that(1)The proposed method
9、 has good accuracy and applicability,which is suitable for the control of staycables in actual bridge engineering.(2)The method is stable and fast in calculation,and the first iterationvalues are relatively close to the final convergence result.After 35 iterations,the calculation can accuratelyconve
10、rge.(3)As the length and vertical angle of the stay cable continue to increase,its nonlinear effectcontinues to strengthen.The method can maintain good accuracy and performance for each stay cable,andhas achieved good application result.This provides a reference for high-precision structural analysi
11、s andprecise manufacturing and installation of stay cables.Key words:bridge engineering;stress state;numerical iteration;stay cable;catenary theory0引言斜拉索是现代斜拉桥最重要的结构组成之一1,准确厘清其受力行为特点对于大跨度斜拉桥的高质量设计、施工以及安全运营均具有十分重要的意义。对于斜拉索安装分析,其无应力长度、倾角及线形是重要的几何控制参数2,针对这方面内容,国内外学者已开展了相应研究。Jayaraman 等3 基于Irvine4 所推导的悬
12、链线索自重状态方程首先提出了弹性悬链线索单元的构建方法。张立新等5 基于idders 改进弦割法迭代技术并结合无应力索长作为迭代控制量,对连续长索的无应力索长求解进行了研究。方志等6 根据索端竖向分力与无应力索长之间关系,利用 Levenberg-Marquardt 算法,提出了一种新的非线性迭代求解技术,该方法精度较高、性能优良,但计算相对繁琐,不便于工程应用7。为解决这一问题,赵雷等8 采用梁端索夹角与索力水平分量作为迭代未知量,结合悬链线理论,对分幅斜拉桥斜拉索的无应力索长确定问题进行了研究,但该方法需要联立迭代求解非线性方程组,计算仍稍显繁琐。此外,Jung 等9 基于索几何线形方程表
13、达式,通过 Newton 迭代法对斜拉索无应力索长进行求解,取得了一定效果。上述研究对于斜拉索无应力索长计算等问题进行了深入探讨,有力推动了斜拉索构件精细化分析进展。但就深入理解其力学行为而言,无应力索长仅是表征斜拉索力学要素的一部分,其还包括索力沿索长变化规律、索线形以及索倾角方程等信息,他们互为补充,共同构成了斜拉索完整的力学要素信息。同时,对于索的几何、力学参数求解,既有方法虽然精度较高,但其一定程度上仍存在不够简便的问题,有必要对其进行进一步探讨。鉴于这一情况,首先根据悬链线理论,研究了斜拉索受力线形方程参数的快速迭代格式,建立了适于工程应用的斜拉索关键几何、力学特征的分析求解方法。在
14、此基础上,将求解结果与非线性有限元方法进行对比,验证了所提出方法的有效性。最后,以一项工程实例,对其进行了相应分析和探讨。1基于悬链线理论的斜拉索受力状态及几何参数求解方法1.1斜拉索线形方程建立对于斜拉索线形,众多学者已进行相关探索。参照其思路10 由索力几何关系出发对其进行描述,见图 1。图 1斜拉索受力及索微元几何示意图Fig.1Illustration of stresses on stay cable and geometry ofcable element鉴于索力 T 水平分量 H 沿长不变,而竖向分量V 沿长增量等于索重增量11,故有dH=d(Tcos)=cos dT Tsin
15、d=0,(1)dV=d(Tsin)=sin dT+Tcos d=dl,(2)式中,为索线重度;dl 为索长微元。将式(1)代入式(2)则Tsin2dcos+Tcos d=Tdcos=dl。(3)为简化方程,对 分析。由图1 索微元关系,可得d=d(tan1y)=1/(1+y2)dy。(4)结合微元几何关系以及 sinh1()微分形式d(sinh1(t)=1/1+t2dt,将式(4)代入式(3)并化简95公路交通科技第 40 卷T1+y2dy=Td(sinh1(y)=dl,(5)该式表征了索力 T 与其几何坐标之间的微分关系。进一步引入T 与H 关系H/T=dx/dl 12,积分可得:y=Hco
16、shxH+C()+D,(6)式中,C 和 D 为积分常数。式(6)即为斜拉索线形方程一般形式10,其在其他文献亦有介绍,但各自推导方法和表达形式不尽相同。为求解式(6)中 H,C 和 D 这 3 个未知参量,令(x1,y1)=(0,0),(x2,y2)=(xjxi,yjyi),代入式(6)4 0=Hcosh(C)+D,(7)y2=Hcoshx2H+C()+D。(8)同时,引入目标索力 Tx=0=Ti,根据 T2=V2+H2,则13 Ti=Hcosh(C)。(9)值得说明的是,上式并非唯一形式的第 3 组边界条件表达,其亦可采用索端夹角8 或索力分量合成6 等形式,但结合后续迭代格式来看,应用该
17、式可获得更简洁的表达形式,从而为后续分析提供方便。1.2参数迭代求解方法由式(7)(9)可以看出,参数 D 由于不存在耦合,因此可直接求解为:D=Ti。(10)而其余两个参数 H 和 C 由于相互耦合均需迭代求解。为此,对式(7)(9)进一步变形,可得C=cosh1Hy2+Ti()x2H,(11)H=Ti/cosh(C)。(12)此时,原超越方程组蜕变为显式迭代格式,确定初始索力水平分量 H0后,即可依次通过式(11)(12)反复迭代来获取 H 及 C 参数值。其具体做法为将 H(i)代入式(11)来获取 C(i),进而由式(12)更新第 i+1 步的迭代值 H(i+1),如此反复,直至二者收
18、敛,此即为简单迭代求解方法。该方法具有计算简便、稳定性好、迭代格式直观的优点,即使对于“松索”,也能够稳健收敛。但这一方法主要不足在于方程收敛速度相对较慢,一般需要约 10 次迭代,H 误差才能小于 1 kN。这与其他显式迭代格式7 相比,虽然计算更简单,但由于未考虑索倾角几何关系辅助迭代,因而效率相对较低。为解决这一问题,对上述迭代方法改进,将式(12)代入式(11),消去参数 C,并进一步变形可得Hcoshx2H+cosh1TiH()y2Ti=0(C 0),(13)Hcoshx2H cosh1TiH()y2Ti=0(C 0)。(14)对于上述两式,由于斜拉索构件参数 C 一般大于0(C0
19、情形多对应于垂度较大倾角为负的悬索构件),因此主要对式(13)分析。此式为一复杂超越方程,仍需迭代求解。一般而言,上述方程多采用Newton-aphson 方法8,14 进行求解,但该方法由于需要对全式求导,计算繁琐,不便于工程应用,因此采用弦截法1516 进行求解。其参数选取得当时,算法具有良好的效率和稳定性。记F(H)=Hcoshx2H+cosh1TiH()y2Ti,(15)则可建立迭代控制方程为:H(i+1)=H(i)F(H(i)F(H(i)F(H(i1)(H(i)H(i1)。(16)式(16)相应的迭代几何示意图如图 2 所示。图 2弦截法参数迭代求解示意图Fig.2Illustrat
20、ion of iterative solution of parameters bysecant method算法具体求解做法为,在第 i 步依据式(16)求得 H(i+1)后,将其与 H(i)配合而舍去 H(i1)再次代入式(16)进行第 i+1 步 H(i+2)更新,如此反复,直至相邻两步 H 结果收敛。由图 2 可以看出,弦截法实质上是通过不断截取弦长通过零点来持续更新 H(i+1),因而其收敛速度相较于之前简单迭代法大为增加。对于迭代初始值的选取,其涉及到两个迭代初始量 H(0)及 H(1)。方便起见,H(0)可依据斜拉索弦长06第 5 期白冰,等:斜拉索受力状态及几何参数的快速迭代分
21、析方法水平夹角确定17:H(0)=Tix2x22+y22。(17)而对于 H(1),由于斜拉索梁端切向倾角一般小于弦长夹角,因而 H(1)应较 H(0)偏大为宜。根据经验,其取值比 H(0)略大时,迭代收敛速度较为良好。此处建议 H(1)为:H(1)=H(0)+20 kN。(18)当计算收敛后,将 H 计算结果代入式(11)解得积分常数 C,斜拉索受力及线形方程即可完整确定,相应无应力索长、索倾角等计算方程亦可由此进一步导得。对于这一迭代方式,其具有格式简洁、无需求导的便捷特点,因而将其应用于工程实际较为方便。但值得说明的是,上述快速迭代方法对于斜拉索这种应力水平较高的主动张拉构件一般较为合适
22、,其性能表现良好。而对于部分张紧程度较小、以被动受力为主的悬索构件,该方法则存在一定的求解困难,这主要是由于索端夹角正负号以及迭代初始值的选取问题而导致,需对二者进行相应修正。但一般而言,这类构件与斜拉索受力特性具有一定差异,因而在实际斜拉桥工程中较为少见(对方法适用性影响有限),即使出现,其亦可采用之前导出的式(11)(12)简单迭代求解方法或依据式(14)形式进行相应迭代处理,该方法对于这类悬索构件计算效果良好,其与快速迭代方法互为补充,可有效解决参数的高精度求解问题。1.3斜拉索无应力长度及索力索倾角计算无应力索长 Lunstress一般由变形后索长 L 减去索弹性伸长量 L 来获得2:
23、Lunstress=LL。对于 L,将式(6)求导变形并积分,可得L=x201+y2dx=Hsinhx2H+C()Hsinh C,(19)该形式最早由靳明君等18 导得。而对于 L,有L=x20TEAdluntress=x20T/EA1+T/EAdl x20TEAdl。(20)由于索应变极其微小,故而式(20)近似不致引起较大误差11。将沿程索力和式(19)代入,则14 L=H22EA12sinh2x2H+2C()+x2HH24EAsinh(2C)。(21)同时,对式(6)求导并结合T 与H 关系,可得 13 索倾角沿程变化规律=tan1sinhxH+C(),(22)索力沿程变化规律T=Hco
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 斜拉索受力 状态 几何 参数 快速 分析 方法 白冰
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。