《脉动风载下的输电线路风偏计算研究》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
分类号:TU973单位代码:1的35密级:学号:21212061.V、硕±学位论文'戀巧;I中文论文题目:脉动风载下的输电线■路风偏计算研究英文论文题目Studyoncaculationofconductorswingingunder:fluctuatingwindloads申请人姓名;悦杨指导教师:楼义娟教授专业名称:结构工程研究方向:结构风工程所在学院一:津筑工棉学院,,y论文提交日期2015年1月-节、,爾VV‘.成戶V:''-■*^。-.人 ̄'-’*'*々...^声.;.巧 化动风载下的输电繞路风偏计义研究论文作者签名:如知指导教师签名:為论文评阅人1;孙炳楠/教授/浙江大学评阅人2;叶尹/教授/浙江大学评阅人3;隐名评审评阅人4=5评阅人:答辩委员会主席:孙炳楠/教授/浙江大学委员1;楼文娟/教授/浙江大学委员2:叶尹/教高/浙江省电力设计研究院委员3;沈国辉/副教授/浙江大学委员4;黄铭枫/副教授/浙江大学委员5:答辩日期:2015年3月9 化动风载下的输电钱骆风偏计算巧究论文作者签名;指导教师签名:论文评阅人1;孙炳楠/教授/浙江大学评阅人2:叶尹/教授/浙江大学评阅人3;隐名评审评阅人4^5评阅人;答辩委员会主席:孙炳楠/教授/浙江大学委员1:楼文娟/教授/浙江大学委员2:叶尹/教高/浙江省电力设计研究院委员3:沈国辉/副教授/浙江大学委员4:黄铭枫/副教授/浙江大学委员5;答辩日期:2015年3月9 浙江大学硕±学位论文独创性声明浙江大学研究生学位论文独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研巧工作及取得的研巧成果。除了文中特别加W标注和致谢的地方外,论文中不包含其他人己经发表或撰写过的研究成果一,也不包含为获得浙江大学或其他教育机构的学位或证书而使用过的材料。与我同王作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。L学位论文作者签名:签字曰期:年月fe^学位论文版权使用授权书本学位论文作者完全了解浙江大学有权保留并向国家有关部口或机构送交本论文的复印件和磁盘,允许论文被查阅和借阅。本人授权浙江大学可W将学位论文的全部或部分内容编入有关数据库进行检索和传播,可W采用影印、缩印或扫描等复制手段保存、汇编学位论文。(保密的学位论文在解密后适用本授权书)学位论文作者签名:导师签名:心f户年/签字日期:年5月日签字日期:於3月日 浙江大学硕壬学位论文致谢王年前。,我有幸来到浙大,来到美丽的浙大紫金港校区,开始我的研究生生涯我知一一道这是一白毎天都需要努为充实自己段美妙的旅程,才能够成为个合格的,我也明,浙大毕业生。时光崔茸一,白骑过隙。下子兰年的时光就巴经走到了尾声。回想起当初刚进入浙大时对自己所提出的要求,我认为这互年过的并没有遗憾。无论是上课、翻阅文献、做风洞试验一一、写小论文亦或是在球场挥洒汗水,在桌前与同学们玩桌游斗智斗勇,幕幕都近一在暇前,毎每想起都会瓜笑。一个幸运的人我想我是,因为我能够在研究生如此重要的阶段遇到了楼文娟老师。楼老师学术能为卓越,总能对问题提出独到的见解,并且为对学生提供新的研究思路。当遇到了学术瓶颈时与楼老师的讨论总能让我有柳暗花明又一村的感觉。当我们因为年轻而,一颗包容的瓜来指导我们产生各种各样的错误时,化我们能够从根本上认,楼老师总是用识到错误,并改正。楼老师不仅是学未上的涵楚,而且也是生活的达人。劳逸结合的生活理念,平和乐观的生活态度使得楼老师的脸上总能法满笑容。在即将毕业之际,我由衷的对老师表示感谢,也祝福楼老师能够身体健康,工作顺利。同时也希望将来在我的人生路途中与老师还能够经常来狂!,继续向老师学习在研究生期间我遇到了一群志同道合的兄弟姐妹与他们一起工作、生活让我觉得人,生是彩龟的。其中我最想感谢的是杨伦博去,,杨伦师兄王作认真,性格开朗学术能为卓一起工作的过程中学习到了许#专业的知识越,并且非常乐于帮助同学。我在与他。此外还有章李刚博±泛在一起总是充满欢声笑语。,章李刚师兄爱巧广,科研能为出众,和他还有徐海巍、姜雄、夏亮、王东、罗置、涂志斌、潘小涛、吕江、段志勇、余江、李天吴、李正吴、梁洪超、刘萌萌等诸位师兄弟,在我的学习上、生活上都给予过我许多帮助,在此一并感谢。我还想在此感谢我的父母,感谢你们二千多年来对于我的养育与照顾,无条件的支持我,在,提供物质条件让我可W专做自己想做的事。在我失落时给我加油鼓劲巧得意时默默在我身后替我感到高兴。你们才是最伟大的人!即将从单纯美好的校园离开,成理有千万不舍。但人总是要成长,即将进入社会的我一将会面对更多的挑战。愿我能够直保持在校园中的这份热情与纯真,W乐观的々态来处I 浙江大学硕壬学位论文致巧一—理今后每个难题,;W颗年轻的来感受这个社会创造出属于我自己的未来。最后!,感谢评阅和出席论文答辩的各位专家老师给予的指导杨悦2014年12月于浙大紫金港II 浙江大学硕壬学位论文^摘要我国疆上辽阔,东西部电为资源分布和经济发展不相匹配。架空输电线路作为主要的电能传输方式其发展日新月异,而导线风偏闪络问题严重影响输电线路的安全运行,因此准确、高效地计算导线的风偏响应具有实际的工程意义。本文在导线非线性动力方程中引入了结构运动速度与风速的相对速度,建立了考虑今动阻尼的风偏计算方法。运用谐波叠加法模拟瞬态风场,W此构建作用于导线上的脉动风荷载。针对某实际八跨线路,建立精细化的多跨导线有限元樸型,基于非线性購态分析方法考察了输电线路自身运动引起的气动阻尼对输电线路风偏瞬态响应的作用特征,结果表明:气动阻尼作用明显,在风偏动态分析中必须考虑气动阻尼的影响。将本文考虑气动阻尼的风偏计算结果与考虑美国输电线路设计规范ASCEN0.74给出的气动阻尼公式计算得到的风偏值进行对比一。结果表明这两种考虑气动阻尼的方法所得的风偏计算结果基本致,ASCEN0.74中所给出的气动阻尼公式可供日后风偏计算时采用。论文研究了目前常用的风偏计算单摆模型的准确性,对比单摆模型与多自由度连续多跨模型风偏响应结果;考察了动态风荷载作用下多自由度连续多跨模型的风偏总响应与平均响应值,确定了动为放大系数取值;将我国规范中导线风偏的风荷载计算方法与美国、英国及国际通用规范进行对比,并吹实际八跨线路为例,计算了各稻导线的风荷载值。结果表明:单摆模型相对多自由度连续多跨模型偏于安全,可在实际设计过程中采用;脉动-风荷载对于导线风偏影响不容忽视..5,其荷载的动为放大系数为131;我国在计算风荷载时所取的风速、体系系数与其他国家相差不大,但关于风压不均匀性及动力放大两个因素的综合考虑与他国相比较小,使得我国风荷载计算值偏小。本文提出导线风偏响应的频域分析方法。考虑导线在平均风荷载作用下的几何非线性,W平均风载作用下的位置座位初始状态,进行脉动风载下的频域分析,W某单跨导线为例进行计算;分别取前20阶、10阶、5阶、3阶、2阶模态进行组合,考察其计算精度。结果表明:频域计算结果与时域结果相近,证明本文所述采用频域法计算导线的风偏响应是可行的;当取前2阶模态作为参与模态进行SRSS法组合,其所得的响应值与时程法计算结果的相对差值约5%,精度较好且计算效率最高。研究导线风偏响应的频域特征,发现其响应主要是背景响应,共振响应很小。采用LRC计算处于平衡位置的单跨导线模型等效静为风荷载。将等效目标响应与施加气动力模III 浙江大学硕壬学位论文楠要型的响应进行对比,结果表明:采用LRC法计算得到的导线等效静为风荷载精确度较高。关键词:风偏;气动阻尼;输电线路规范;等效;频域法静为风荷载IV 渐江大学巧壬学位论文AbstractAbstractInourcoxmtrythehenomenonofimbalancedowerresourcesisserious.Asthemaor,ppjowertransfermethodTransmissionsstemdevelosdramaticallbutitisseverelthreatenedp,yp乂ybflashoveraccidentswhichwerecausedbconductorswinin.Thereforetheresearchonyygg,accurateandeficientcomputationmethodsofthedisplacementofinsulatorhasgreatracticalpsncancenigi行inengineerig.Consideredtherelativevelocityofconductorandwind,过methodtocalculateconductor-swinging化atafectedbyaerodynamicdampingwasestablished.Then化etimedependingfluctuatingwi打dloadswouldact〇打tra打smissio打lineweresimulatedwiththecombinationofuas-harmonicwavesuperstitionmethodandqisteadyassumption.Anonlineardynamic-transmissionlinemodelwhichwasconsistedofeightspanelectricalconductorswasestablishedpreciselybyusingfiniteelementmethod.Furthermore,theimpactofaerodynamicdampingcausedbconductormovementsondnamicresonseofthetransmissionlinewasdiscussedinyypdetail.Itisshownthattheaerodnamicdampingcanreducethemaximumvalueofconductoryswininsinificantlwhe打calculatedislacementofinsulatorbdnamicanalsismusttakegggy,pyyytheaerodynamicdampingintoconcern.Ontheotherhand,theaerodynamicdampinggainedwiththeUStransmissionlinedesigncodeASCENo.74wasadopted化computethe姐splacementofinsulator.Theresultoftwomethods化calculatetheresponseofinsulatorunderaerodynamicdampingwerecompared.Itshows化at:化eresultsarecloseU)eacho化er,aerodynamicdampingainedwithASCENo.74canbeusedinfutureresearch.gComared-thediferencebetweenMultisanmodalandsinleendulummodalppgp,investigatetheaccuracyofsinleendulummodal.ThrouhdiscussedthediferencebetweengpgmaximalvalueandthemeanvalueofconductorswingingcalculationresuUsthisthesis,researched化e-imactofturbule打twindvelocityonthewindinduceddnamicamlificationpypefect.虹orderU)have泣beterunderstandi打gofwindloadcalculationinthecomutationofpconductorswingingthesetinofcalculationarameterswerecomaredamonfourdiferent,gppgationaesncodes-ldi:《GB50545》、《ASCE、NgNo.74》《lEC608262003》and《BS8100》.ermoreasedon-Furth,barealeightsantransmissionlinethediferencesamonthecalculationp,gresultsofdiferentcodeswereanalyzedunderthesamedesigncondition.Theresultshown化at:口sing化esinleendulummodaltocalculate化erotationanleofsusensio打insulatorstringisgpgpsafeThednamc?imanificationfactorisabout1.31.5windloaddefinedinChinesedesin;yg;gcodeisthesmallestthemai打causeliesinthele巧ofcomrehensiveconsiderationofwind,ppressureasymmetricalcoeficientanddynamicmagnificationfactor.V 浙江大学硕去学位论文Abstractfreuenc-domaAqyinanalysismethodofconductorswingingispresentedinthispaper.Geometricnonlinearityisconsideredwhen也emeanwindloadactingon化ewires,A-Mo虹si打lesantransmissionlinewastakenastheresearchmodalreoverordertoimrovegp,,pcalculationeficienctwenttenfivethreetwomodevibratio打wassuerosed.TheKsu化y,y、、、、ppshownthat:Theresponseofconductorswinginthefreuencdomainandintimedomainisqycoseoeach-ltotherthefreuencdomai打analysismethodhasrovedtobefeasible了aki打the;,qypgfirstandsecondmodelfbrS民SScombination,ahighPrecisionconductorresponsecanbegotwiththebestcalculationefficiency.ThisthesispresentsthefrequencydomaincharacterofinsulatorStri打gsswing,foundout化atwinddynamicresponsepredominantlybackground,and化ereso打ancecanbeignored.Load-Resonse-Correplatio打(LRC)methodwasadopted化calculatetheequivalentstaticwind--loadESWL)ofasinlesantransmissionlinetheresonseofthesinlesantransmission(gp,pgp-linethatactedESWLorthetimedependingfluctuatinwindloadswerecomared.Theresultgpshownthat;TheESWLofconductorthatbasedo打theLRCmethodisrecisel.pyKeywords:Conductorswinging;Aerodynamicdamping;Transmissionlinestandard;FrequencyDomainMethod;EuivalentstaticwindloadqVI 浙江大学硕±学位论文目次^I摘要IIIAbstractV0次i第1章绪论11.1研究背景11丄1我国输电线路的发展状况11丄2风偏闪络对输电线路的影响21.2研究现状31.2.1国外学者对于导线风偏研巧31.2.2国内学者采用静力方法分析风偏5123..国内学者采用动力方法分析风偏51.2.4国内外输电线路规范关于风偏计算风荷载取值的差别61.3本文主要工作71.3.1考虑气动阻尼效应的输电线路动力风偏计算71.3.2我国规范风偏计算方法及与多国规范对比71.3.3利用频域法计算导线风偏8134..导线的等效静力风荷载8第2章考虑气动阻尼效应的输电线路风偏动态计算92.1引百92.2考虑气动阻尼的动力学公式923有限元模型的建立112.3.1八跨线路的原始参数112.3.2有限元模型建立口21.3.2输电线路动力特性32.4风速的模拟142.4.1平均风速模拟152.4.2脉动贈模拟152.4.3风速模拟结果172.5有限元结果分析182乂采用ASCENo.74所提供气动阻尼202乂.1ASCENo.74中的气动阻尼公式20I 浙江大学硕壬学位论文262SCENo74..采用A.气动阻尼有限元结果212.7本章小结22第3章我国规范风偏计算方法及与多国规范对比243.1引言243225.风偏计算模型讨论3.2.1弦多边橫型253.2.2刚体直棒模型263.2.3多自由度连续多跨模型273.2.4单摆模型正确性讨论273.3动力放大系数取值283.4#国规范线路风偏风精载的对比293.4.1输电线路风前载计算么式293.4.2巧载系数/313.4.3风速的取值313.4.3体型系数C37d3.4.4巧系数对比373.4.5八跨模型按多国规范计算风荷载值403.5本章小结42第4章利用频域法计算导线风偏444.1引言444.2基于频域算法的导线风偏响应基本理论454.2.1频域法计算导线风偏响应的过程454.2.2导线使用频域算法的条件讨论484.3单跨导线算例494.3.1单跨导线的有限元建模494.3.2导线模型施加平均风荷载504.3.3频域法计算结果544.4本章小结57第5章导线的等效静力风荷载595.1引言595.2风致响应分类605.3导线风偏响应在频域上的表现615.4导线风致响应的计算方法及其对应的等效静力风荷载抗5.4.1平均响应625.4.2背景响应树 浙江大学硕壬学位论文目次564.4.3总响应565.5导线等效静力风荷载求解5-6本章小结67第6章结论与展望的6.1本文主要结论的6一.1进步王作展望71参考文献72作者简历77攻读硕壬期间科研成果77ill 浙江大学硕壬学位论文第1章绪论第1章绪论1.1研究背景1丄1我国输电线路的发展状况我国疆上幅员迂阔,东西跨度约5200么里。东部地区交通便利、经济发达,发展过程中所消耗的能源约占全国总能源消耗量的兰分之二,而我国大约四分之吉的煤炭资源与五分之四的水能资源集中在中西部地区。利用西部资源,建立水电和火电基地通过清洁能源""W电能,实行西电东送来解决能源基地与电为需求逆向分布的问题,不仅可W满足东部地区发展的电为需求一,还能带动西部的经济发展,是国家的项长期战略。为了能够切"一-实贯彻西电东送如图11)这战略,自改革开放W来我国输电线路建设速度惊人,从1949年直交流输电总线路长度6500km200785、,发展至年仅直流输电线路就有70km,一一跃成为直流输电线路总长度世界第。再到2013年底220kv,及W上输电线路总长度这53%000km:其中包括交流线路51810(X)km,直流线路217000km。至此我国电网规模已稳居世界首位。由于大规模的电网建设,我国输电线路大多处于长距离大容量的送电环境,且所经过的各地区地貌、气象等自然环境差别较大,随之而来的恶劣天气、人为行动、鸟害等等都对架空输电线路的安全运行造成了巨大的挑战。其中输电线路风偏引起的闪络问题发生频繁且造成的危害较大,越来越受到人们的重视。?-、公:::::繫^^5細/二;*葉&^'节I'U媒,—…图1-1西电东送示意图I 浙江大学硕±学位论文第1章绪论112..风偏闪络对输电线路的影响输电线路是电网的重要部分,在田野、山谷、丘陵、河流等地带均有分布,具有量大面广的特点。架空输电线路是由输电塔、输电导线、避雷线、绝缘子和线路金具等组成。架空输电线路在大风降雨等极端恶劣气象条件作用下,悬垂绝缘子串与导线将产生出平面的摇摆。当摇摆幅度过大,使不同相的导线或导线与输电杆塔之间的距离小于最小安全电(-2气间隙时,,引起,就会发生放电闪络事故导致构件表面烧蚀如图1所示)线路跳间、电网断电等事故,造成巨额的经济损失和重大的社会影响。S ̄f3IIJ/rl?遍、.'51(a)导线上闪络痕迹(b)塔材闪络痕迹瑪.^铜- ̄!户懸(C)塔材闪络痕迹(d)闪络位置远景图片图-12输电线路风偏闪络事故后照片输电线路的风偏放电是影响线路安全运行的主要问题。风偏闪络事故若发生于主干输999一2003电线路上将造成大面积停电,500kV,严重影响电力系统供电。在1年年间输W电线路风偏仅发生31次。随着电网规模的扩大,长距离、高电压的输电线路不断落成,2005一2011年间也使得风偏导致的断电现象频频发生。据不完全统计,年全国11066kV,()一1^751376处500及乂上输电线路发生风偏跳闽运条次,其kV,故障停运这条次。仅河南2 浙江大学硕壬学位论文第1章绪论线路在近10年来发生风偏跳阔20余次。当架空输电线路发生的故障清除后,在短时间内通过闭合断路器进行重合间,可使输电线路在短暂故障后继续工作,这是线路在运行过程中常采用的自恢复供电方法,而输电线路在发生风偏闪络后重合巧成功率很低。这是由于风偏跳间大多是在强风天气或中小尺度的局部地区产生贿线风的条件下发生,输电线路受此类风的影响持续时间较长,使得整个重合间动作过程都处于风偏状态下,放电间隙较小;同时,输电线路在进行重合间动作一时,出现定幅值的操作过电压,这将导致再次放电。因此,线路若发生风偏跳间,重合W闽成功率极低,很大程度下影响了供电的可靠性。由上述分析可知,风偏事故己是影响电网正常运行的重大安全隐患。同时受风偏影响的输电线路难W合阀成功路一旦发生风偏闪络现象即会造成可观的经济损失。悬垂绝,线缘子风偏角度过大是导致风偏事故发生的直接原因,通过计算得到的风偏值大小直接关系到电气间隙设定和输电线路建设成本,因此深入研究输电线路风偏闪络的原因并提出简单可靠的风偏计算方法W及合理的相关参数取值,有着非常重要的工程实际意义和技术经济效盜。针对近年来发生的风偏跳间事故内相关领域的专家进行了研究与分析,国,认为主要是由于外因和内因两方面因素造成的。外因是自然界发生的强风和暴雨天气;内因则为输电线路抵御强风的能力不足。找出影响风偏跳阔的关键因素,采取有针对性的方法和措施,便可提高线路的正常运行水平。1.2研究现状为避免频繁地发生风偏现象一,国内外相关学者都对此做出了努为并开展了系列的,研究王作,大多数研究的最终目的就是希望能够确定如何正确计算风偏值。现如今风偏的计算方法主要有静为计算法W及动力计算法两种。1.2.1国外学者对于导线风偏研究国外对输电线路风偏计算值的研究起步较早,其中日本和部分北美国家针对输电线路风偏W’yW问题率先展开了实测研究和理论分析,为日后风偏研究奠定了基础。7tlTsumo〇-jit基于真型输电线路的风偏观测W及理论分析,提出了W图13所示离散自由度导线模型来计算导线风偏值,并使用计算机有限元软件模拟了某单跨导线,发现与实测结果接近。3 浙江大学硕去学位论文第1章绪论^S>/化1V玄面:alumpedmassfswindpressureforceTn:tesionSs色panlength图-Tsumoto13ji所提出的导线风偏计算模型民化保证悬垂串导线与巧塔之间的最小安全间隙为设计准则,利用图解法计算了导线风偏角。一(X在风偏规范计算么式中,存在个风压不均勾系数。导线是大跨度长距离的结构,因此在一定范围内导线的平均风压并不能处处相等a就是考虑了这个因。风压不均勾系数素后对导线平均风压的折减系数。国外学者对风压不均勾系数也做了许多研究:,"C1ha/对于不同档距的单跨输电线路进行了研究,讨论了线路风荷载的巧减方法,并。指出随着风场嫌流度及档距的增加,风压不均匀系数会相应减小4[]-M-Hiratsuka对台风影响下多跨输电线路iyakoimatestline(如图14)进行了实测,j获得了lOmin的平均风速,并观察其风偏响应的结果,并利用风洞试验,研究了时距的长、短对于风压不均匀系数的影响。.辦1!薩I ̄l4Hratsuka图i所采用的实测线路4 浙江大学巧±学位论文第1章绪论*MO?〇三IIrmS11III11 ̄ ̄ ̄ ̄ ̄巧 ̄ ̄Z.IxirnIa韦?pfiMSWIO|^I.《广T一—<〇。一言^巧_>i——巧〇ghorizontal香:。'’。—’扣s皇mini![jj^旨一/?-?-'’’一…??-…?i:!i;:肚±扛±:±r》?古i韦禅年f毕I11111*-S'吾受至II!iastesaveraevaseed/s!110mgidpd{)--ratsu图15Hiratsuka实测获得的lOmin平均风速图16Hika得到的风偏角122..国内学者采用静力方法分析风偏国内的学者化对风偏计算研究做出了许多贡献。相对于采用动力学计算方法,此早期的风偏问题静为学方法拥有简单易懂,因,适用于实际工程计算的独特优点研究多采用静力模型分析的方法:龙立宏、胡毅等从实际问题出发,分析了发,生风偏闪络事故的主要原因,并应用静为学方法计算在不同工况下的风偏值大小得出风速、风压高度变化系数等四个参数对于、风压不均勾系数、风向与导线夹角12[1分析了已有的绝缘子串风偏角静力计算模型风偏值的影响,并考虑实际运行;李孟春中的输电线路绝缘子串受为状况,提出了新的绝缘子串风偏角静力计算模型,并进行对比一分析13中提出种最大风偏角的修正方,从而减少模型误差。王声学,吴广宁等在义献[]法并且骄弃了计算杆塔与导线之间最小空间间隙复杂的传统几何方法,取而代之的是在,二,而且还抒咨的模型之中引入笛卡尔维坐标,这样不仅能够使得计算更加简洁方便能够14[]提高计算的精度、张志远认为间隙圆图解法这种传统的计算风偏角的方法。此外李瑞祥工作量大一种新的图解法来计算输电线路在运行电,精度差且不利于编程求解,而提出了压、巧过电压、外过电压、带电作业等工况下的风偏角度。1.2.3围内学者采用动力方法分析风偏虽然静为法拥有独特的计算优势,但在实际情况中输电线路风偏运动却表现出明显的动为特性。自然风包含了两部分:。这是因为导线是受到了自然风的影响而产生风偏运动一平均风与脉动风。平均风可简单的认为是在段时间内速度不改变的风,其对导线的作用可近似为静态;脉动风是指风速随时间不断改变的风,其对导线的动态影响就不可忽略。一,此外输电线路是种明显非线性连续结构,无法直接采用频域法来计算风偏值因此需要利用时域法来考虑其动为特性。采用时域法时首先遇到的问题就是如何通过数值模拟来获5 浙江大学硕去学位论文第1章绪论—15得包含空间相关性的时变脉动风荷载[]。目前国内外采用MonteCarlo思想来模拟脉动风-i623t]一速的主要有:谐波叠加法与线性滤波法等两种方法。谐波叠加法是通过叠加组具有随机相位角的王角级数来模拟随机变量的过程,这种方法数学逻揖严密,算法直观,当采PW用FFT技术之后可在很大程度上提高谐波叠加法的计算效率。线性滤波法又称为白噪声滤波法一噪声进行适当。其过程是把随机过程抽象成组满足条件的白噪声,然后将白PS1变换而拟合出时程。拥有了模拟风荷载的理论知识,且随着计算机水平W及有限元软件的不断升级,越来越多的国内学者开始基于ANSYS、ABAQUS等有限元软件建立线路模型进行导线风偏值的计算。针对传统静力计算方法的不足,文献口巧在ABAQUS中建立了输电线路的有限元模型,在模型上分别施加典型阵风与脉动风,并将计算出的风偏值与输由线路设计规范方法确定的风偏值进行比较,证明了在脉动风的作用下有限元动为模型计算得到的风偏值一将大于规范值,推断这可能是风偏闪络事故频发的原因之。这是国内较早进行动态风偏研究的成果。在此基础上,林雪松模拟脉动风时采用Kaimal谱和Davenport相干亟数,弥补了义献口巧中未考虑风速随高度变化W及风的相干性的问题,再将模拟得到的脉动风加载到有限元模型上,与规范算法进行对比,讨论了现行架空输电线路规范中风偏角一计算的不足之处。此后还有:孙宝强建立了个兰自由度的多挡距动志模型并在模型,上施加不同风速的风荷载,从而得到输电线路风偏的动态特性,并成功地找到了最小间隔距离与风速的关系PW;邵瑰讳巧约借助于有限元模型分析不同愚挂方式的绝缘子串对于,PUPSs导线挂点处风荷载的影响艳IPiPWPy;郑佳、刘小会、罗先国、李黎等人也都使用了动为分析方法来计算架空输电线路的风偏问题,但是在计算过程中都未能考虑由导线运动所产生的气动阻尼影响。楼文娟、杨悦首次提出了考虑导线与风速的相对运动所产生的气动阻尼,发现气动阻尼对风偏响应的影响不容忽视,在采用动态方法计算风偏响应时必须考虑气动阻尼的影响。1.2.4围内外输电线路规范关于风偏计算风荷载取值的差别我国风偏事故频发:,但有调研表示在美、日等国家输电线路发生风偏闪络的现象W1较少,而我国规范所采用的计算模型与国外规范提供的计算模型相同,因此学者们开始研究国内外规范中:设计基本风速、风压高度变化系数、风压不均匀系数、微地形影响等计算参数的差别,试图从中找出我国风偏事故多发的原因。闽鉤,邵瑰讳等通过对出中、美、日规范中关于基础风速的设定差别,得出我国风速换算时距长,设计基准风速值较小;6 浙江大学硕去学位论文第1章绪论罗列出几个规范中不均匀系数的取值么式,指出我国关于不均匀系数的取值较为合理;对比微地形参数,指明我国规范设计的安全度相对国外规范要小等结论。还有《位学者化进行了围内外风荷载规范的对比,但其着重点在于对比适用于输电塔这类构筑物或高层建筑物的风荷载规范的差别。关于导线风偏荷载计算的规范对比,相关学者研究较少且不够全面。13本文主工作.要1.31.考虑气动阻尼效应的输电线路动为风偏计算在研究考虑气动阻尼效应的输电线路风偏动患计算时,本文主要做了W下工作:(1)针对连续多跨输电线路在瞬态风场作用下的风偏问题对导线非线性动为方程进行,研究,考虑结构运动速度与风速的相对速度,提出了考虑气动阻尼效应的输电线路风偏动态分析方法,为开展导线风偏的数值模拟提供了依据。kV八跨线路一(2)W某实际500为原型建立了个精细化的非线性动力有限元模型采用,谐波叠加法模拟了导线各点的脉动风速,在时域内计算考虑气动阻尼与不考虑气动阻尼的顺风向风偏响应,并对结果进行分析。(3)在有限元程序中考虑美国输电线路设计规范ASCENo.74《Guidelin故forElectricalTransmissionLineStructuralLoading》所提供的气动阻尼,获得其在时域上的风偏响应,并与本文所述考虑气动阻尼的数值摸拟结果进行对比。验证了ASCENo.74中气动阻尼计算么式的准确性。132..我国规范风偏计算方法及与多国规范对比关于规范算法研究及多图输电线路风偏风荷载计算方法对比,本文主要做了兰部分工作;(1)为考察单摆模型的准确化,首先采用静力单摆模型根据规范给出的导线风荷,载计算公式及相关计算参数的取值方法,求解了输电线路的静态风偏响应,并与采用有限元法求解得到的八跨线路的平均风偏响应进行对化。(2)将有限元模型计算得到的风偏总响应与平均响应进行比较确定了动为放大系,数的取值。7■J 浙江大学硕去学位论文第1章绪论-(3)分别介绍中国输电线路规范《GB505452010》、美国输电线路规范《ASCENo74》60826 ̄2003?.电工协会输电线路规范《圧C英国输电线路规范《BS8100》;国际;中关于导线风偏时风荷载的计算么式一个算式进行参数的对比分,并将这些公式归纳为,析。最后W八跨线路为例,分别根据各国规范计算各档导线的风荷载值,直观地表现出我国导线风荷载计算值与他国的差别。1.3.3利用频域法计算导线风偏在利用频域法来计算导线风偏响应值时,本文主要工作为:(1)提出可用于频域法计算风偏响应的导线摸型并通过某单跨导线模型进行检验。.,(2)在频域法计算中分别取20阶、10阶、5阶、3阶、2阶模态作为参与模态进行组合,得到在保证精确度的前提下最具计算效率的姐合模态阱数。1.3.4导线的等效静力风荷载本文求解了导线的等效静为风荷载,研究主要内容如下:(1)考察导线风偏响应的过程在频域上的表现,研究背景分量与共振分量对于总响应的贡献。(2)采用LRC计算处于平衡位置的单跨导线模型等效静力风荷载。并将等效目标响应与施加气动为模型的响应进行对比,验证采用LRC法来计算导线等效静力风荷载的有效性。8 浙江大学硕±学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算第2章考虑气动脏尼效应的输电线路风偏动态计算2.1引言设计风偏角的大小直接关系到电气间隙和输电线酪建设成本,因此风偏计算方法的精确性和相关参数的取值合理与否,对输电线路的安全运行和经济合理至关重要。国内外均对输电线路风偏做了许多研究,但研究过程主要采用静为模型进行风偏计算。导线是明显的非线性结构,当其发生风偏时表现出巧显的非线性特征。现如今计算机性能飞速发展、有限元软件的效率也越来越高,通过有限元软件来模拟导线发生风偏的过程,能够较高程度的拟合出现实中的情况。因此近年来国巧多位学者提出使用有限元软件,且考虑脉动风速对输电线路的动力效应,进行了风偏动奈响应分析但其在棋拟过程中,忽略了气动阻尼的影响,从而高估了风偏的动态响应。PilHomles曾通过数学方法给出与塔顶位移有关的格构式输电塔气动阻尼计算么式,并指出在研究风致输电塔气动效应时,气动阻尼作用显著。导线的风偏位移相对格构式输电塔要大,因此气动阻尼效应更加明显,对风偏的影响不可忽略。本章研究在输电线路非线性动为方程中,引入结构与来流之间的相对运动,建立了考虑气动阻尼的风偏计算方法。运用多变量谐波叠加法,考虑风场的时间和空间相关性特征,模拟强风气象条件下的畴恣风场,并W此构建作用于输电线路上的脉动风荷载分布模式;同时W风偏闪络事故多发的某500kV输电线路为研究对象,考虑支座高差、绝缘子串长度和几何非线性效应等因素的影响,采用ANSYS软件,建立包括整个耐张段的多控距输电线路精细化有限元模型,基于非线性瞬态分析方法揭示考察输电线路自身运动引起的气动阻尼对输电线路风偏瞬态响应的作用特征,并与考虑美国输电线路设计规范ASCEN0.74脚《Guidelin放forElectricalTransmissionLineStructuralLoading》给出的气动阻尼么式计算得到的风偏值进行了对比。2.2考虑气动阻尼的动力学公式ssLotiredoSouza的研究表明:结构发生风致振动时结构与来流之间的相对运动效应,会使动为系统的阻尼增大,而新增加的阻尼即为附加于结构上的气动阻尼,该类阻尼具有随平均风速增大而增大的特征。而输电线路具有跨度长、阻尼小W及柔度大的特征,属于9 浙江大学硕击学位论文第2章考虑气■动阻尼效应的輸电线路风偏动态计算典型的对风荷载敏感的非线性结构,风偏位移大,故在求解导线风偏响应时应当将气动阻尼的影响考虑在内。乂K向//扔^扔>1图2-1多自由度导线俯视图2-示的《自由度为详细说明导线气动阻尼的产生机理和作用特征,在此引入如图1所多跨输电线路,在风荷载作用下顺风巧位移向量r可表示为:=-?2-YYYY1[??j()?为第;其中K/点的顺风向位移,相应输电线路的非线性运动方程可W写为MY=-2+C+KYYF()口)^^,j式中M和为线路结构自身的质量矩阵和阻尼矩阵;F和y分别是系统的加速度和^5有关:速度向量,与结构的位移状态;i:7{巧是考虑了几何非线性效应的切线刚度矩阵=Y-KYK+K+KY3,(),,()^i)口)。f其中&、尤:线性刚度矩阵,非线性刚度矩阵和应为刚化矩阵;v{巧和估(0分别是为作用于输电线路结构上的风荷载时程矩阵。当不考虑结构与来流的相对运动状态时,可W表示为:222T-=...4F^C+,口+6,口+/〇M)22),兩)]口),吹,却心:.心;,。^?式中口.、,分别、5和4分别是第/点的平均风速脉动风速和有效迎风面积庐和。W,2-4不难发现是空气密度和导线阻为系数。由式(,不考虑气动阻尼时,作用于导线上的风)荷载仅与风速和截面的气动参数有关。2-:当考虑导线与来流的相对运动效应时,式(巧右端的风荷载F为F=+-口+G-[,成。詳卽,4欣1却心;22备..。。2ZO(J—-???八朽+"卽^,。。10 浙江大学硕±学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算将式口-5展开-5),得到与导线自身运动速度}相关的项即为气动阻尼。将式(2)代入式-P巧,整理后可得:MY=-FY+C+C)Y+KY+F6(,,()?_,,口)上式右边jFi和&分别为:22=+,+^^,巧成。[4巧|却)年巧22),,臺_7口)…+5,灯城"。2:=_F成…2,。(如,如:如許(別^由于是导线运动速度的平方项?,为小量可将其忽略。则左边C。。由导线自身,。即为运动状态所致的气动阻尼矩阵:-C=航础口+"口+()()。"。成。4,,,,4,2馬'…+,,4W]■吹,.由上式不难发现,气动阻尼的大小与输电线路周边的风速密切相关。另外,由于线路各点的平均风速均巧大于脉动风速1^的绝对值,因此附加于结构上的气动阻尼始终为正。,,23有限元模型的建立.2.3.1八跨线路的原始参数挂点7挂巧挂巧2^^=I====i48h42h36I33h33hM2-i!lh33!!|!!jI"L丄丄丄1丄」L!l=。*?心228L3261331LM"LM42L393^426L33l图2-2500kV某线路N3 ̄NH耐张段参数示意图?N为研究上文中所提出的气动阻尼的作用500kV实际线路N3il耐张,本文采用某段为研究对象,该耐张段全长2918m.是由两端2个耐张巧塔和中间7个直线巧塔组成的8档线路,处于B类地貌情况下,整跨巧无转角。各跨的挡距、呼高如图4所示,年平均运2506kN4XLGJ-4005-1行张力为.。导线型号为/3,其物理参数见表2。II 浙江大学硕壬学位论文第2幸考虑气■动阻尼效应的输电线路风偏动态计算2-表1LGJ400/35钢志错绞线物理参数 ̄ ̄ ̄ ̄ ̄7ZZZ2弹性系数,mm化命An?计算截面/线密度/kg/kmmm/Mpa425.2426.82650001349---mm直线塔的悬挂绝缘子串信号为FXBW500/180伞颈为190160110患棒材料为,,环氧树脂玻璃纤维引棒,结高4.7m。耐张塔与直线塔上的绝缘子串的物理参数料如下表2-2所示。.表2-2耐张段上绝缘子串物理参数长度/量/位置泊松比fmkg耐张0257—641140绝缘塔子串直f0.3534^2.3.2有限元樸型建立输电线路处于静力平衡状态下时,受到了运行水平张为W及自重的作用。在受到这两种力的作用下,导线精确的几何形态可采用悬链线方程描述。=竺-y--CO化(210acodia)[]()qL式中:y为计算点垂度,义为计算点到挂点的水平距离H是年运行水平张力;为,q'i=+导线单位长度质量;asinh[/?(C/L)/sinhy0]y(?;L为跨距,C是导线挂点高2H差。悬链线虽能准确表示导线在自然状态下状态,但计算较为复杂,,为便于分析当导线^跨中垂度跨距I之比小于化1时导线的形状可1乂近似看成拋物线与实际情况悬链线/与,fwi相差很小,采用公式:4.z王-巧/(y+三-=义口n)2LL:式中,/为跨中垂度,采用下式计算=-。/瓷(2)通过式(2--11)及(212)即可确定出W抛物线表示的导线受重为荷载悬挂情况。考虑计算效率及精确性的问题,导线的每个单元长10m。整个耐张段导线部分共采用了30012 浙江大学巧主学位论文第2章考虑^动阻尼效应的输电线路风偏动态计算个节点。这样就得到了有限元各节点的初始坐标。在导线发生风偏闪络时主要是发生顺风向的位移,输电塔风致振动响应对导线风偏的レレ影响不大,基本可乂忽略。考虑到计算分析效率乂及计算机性能条件,在建立有限元模型时不考虑输电杆塔的影响。建模时将绝缘子串和杆塔之间的连接作为固定较接。虽然在大一风情况下,多分裂导线的各子导线存在定的不同步运动或次档距振荡,但主要仍W整体运动为主一,因此在研究绝缘子串的风偏闪络问题时,将四分裂导线等效为根单导线。输电导线属于非线性特征显著的高柔度结构,因此在现实情况下仅能受拉而不能受压且无法承受弯矩,故本文采用ANSYS有限元软件提供的LinklO单元模拟;而对于风偏发生时始终受拉的绝缘子串来说,可采用Link8单元模拟,其重量W附加质量的形式融合到有限元模型中,具体实施时采用Mass21单元模拟。考虑到导线发生风偏时,绝缘子串的运动syP7整体摆动为主f],若将其简化为刚度较大的杆,结果仍然是合理的故悬垂绝缘子串,由两个节点简化而成。整个耐张段绝缘子与导线部分共由307个节点所组成通过材料将,这些节点相连,形成单元。在导线的有限元模型上施加重为荷载,并于ANSYS找形分析一后更新坐标,差看找型前后两次坐标的差值,若差值较大还需进行下次找型。导线初始张力取年平均运行张力,为25.06kN,结构阻尼取化005。通过上述过程所建立的有限元模型如下图所示:,/户■.f.2- ̄图3N3N11耐张段有限元模型当输电线路受风荷载作用时,导线轴向应变很小但位移很大,对其进行分析时需要考虑几何非线性效应。本文通过激活ANSYS有限元软件中的大变形和应为刚化选项来考虑该效应的影响。2.3.3输电线路动为特性多跨线路的前六阶振型和频率如图2-4所示L乂看出多跨线路的振型主要W平面外,可摆动为主,且对应的频率较低。同时还可发现,在直线塔绝缘子串的连接作用下,相邻跨13 浙江大学硕壬学位论文第2章考虑气■动阻尼效应的输电线路风偏动态计算。另外,其导线之间的振动模式存在显著的传递效应,因多档距输电线路平面外刚度较低2031Hz前阶频率均低于化。—二0.1338細Zb0(a阶阶.1409加Z)()_0Hd四阶0.15770HzC£阶.14662z()()e五阶0.16345H.1629姐z0z()识六阶图2>4输电线路前六阶振型2.4风速的模拟风荷载模拟准确与否是直接关系到风偏响应计算结果是否合理可靠的关键。风速时程!^可乂分为平均风速和脉动风速两部分。平均风速通常被视作不变量,采用沿高度变化的指数或对数规律来描述,基;脉动风速则采用平稳高斯随机过程模拟。在风速模拟完成之后:于准定常假设,可将风速转换为作用于结构表面的风荷载-F=DLCU+213)y非nD(:W(式中,Fy是作用在导线上的时变风荷载,戶为空气密度,D为导线等效迎风直径,I。为加载点控制长度C为阻为系数C7和fi分别为Z高度处的平均风速和脉动风速;,d,,,41 浙江大学硕±学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算2.4.1平均风速模拟56[]计算平均风速时根据《建筑结构荷载规范-,》按式(214)计算。U=x-(214:K〇亩()^上式中,2为各点高度,^1〇为1〇111高度处风速,a是浦流度,根据规范规定B类地2'貌情况下选取0.15。根据此耐张段所在地区,选取基本风压化45kN/m,计算得到10m高57;-[]度处的平均风速,为%.83m/s结合110kV750kV架空输》。00kV,《电线路设计规范5输电线路的设计基本风速不宜小于27m/s,故取基本风速为27m/s。2.4.2脉动风速模拟-式(213)中的脉动风速可基于Davenport脉动风功率谱,采用谐波叠加法模拟:'X-(215)24kvi〇l+x(卢=W为式中,x,圆频率,n为频率,A为与地貌有关的常数。V|〇需要指出的是,虽然Davenort谱具有典型的区域代表性p,但是该风谱所代表的能量在不同高度处是完全相同的,因此不能反映风场浦流度沿高度的变化规律。另外,虽然相-5关资料已给出式(21)中地貌系数A的建议取值但是脉动风速模拟值的均方根等统计特,征并不一定能够与理论澈流剖面相吻合上所述。因此,综,需要根据结构所在地域大气瑞流沿高度的变化特征,对Davenport谱进行修正。众所周知,端流度/与风谱存在W下关,系:=。-1216王"化()仍y2-从17式(()--结合式口1巧和(217)虑沿高度变化的端流度和风速剖面,考,可定义风谱修正系数狂:=逊-2-18&苦()S(〇d(〇£()式中:/为规范给出的滿流剖面理论值。修正后的Daven2;rort谱可表示为;,p15 浙江大学硕壬学位论文第2章考虎气动阻尼效应的输电线路风偏动态计算S=-za>S(〇219{,)Q{)()^由上式可W看出enort谱可有效表征脉动风速能量沿高度的变化规律。,修正后的Davp一脉动风速互功率谱密度函数表示空间两点脉动风速之前相关程度,般为由复数组成的矩阵:'W。、W=weww=o>S仍Cohe2-20()()()()yj如)()馬&jj||拓式中Co/i(co)采用Davenport推荐的风场空间相关函数:222---仍C片+'y)+Gzz如巧所*(,*)1/。。f,、,/、V=<=-2cA仍ex>口。()p52Vz其中c,、c、G分别为沿导线方向;水平垂直于导线方向;竖直方向的衰减系数。,一===这竖值都需要通过试验确定。般情况下推荐可取&16,C6,G10。,根据,'一一=...考虑组《个维随机过程"12》其互谱密度矩阵:(,/〇_/,,,)…、|(的&2(的50Sa???5?2,()^i)2?()=-22A曲口))■■■S仍…S仍(的,2()。。()=式中:0为圓频率,似人义...,/〇为互功率谱函数。将S(份)进行Cholesk分解;()(/y*rS似=//份好仍2-23()()()()其中:■00???0巧1()H...0物。(的H=2-24()抑H仿曲…H仍J)()m()_为下王角矩阵a一。根据Shinozuk理论维多变量脉动风速场可采用下式进行,模拟:=--Vzt2(Aco)H〇})cos(〇)0co)+25?)口){y,,)(?,(?^?jjj^X;^j?Jj,,\| ̄=ml/1:2.3.2节中输电线路的动力特性分析线风偏式中,考虑到导A为截止圆频率,根据=主要W面外摆动为主,同时前20阶振型对应的频率均小于0.31Hz,取嘴27crad/s;是16 浙江大学巧±学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算频率划分段数jV越大模拟结果越合理,为保证脉动风速时程中低频部分不被过滤掉同时,/^Aw=V=428x考虑计算效率取W2048Aw为频率间隔w/?.881257rl(rHz,;,;广?/扭是均匀分布于0w0=[0,2对的随即相位角;气?(是巧?(似)的复角,i/()是实数矩阵,故_?()〇。)巧根据采样定理,风速模拟模拟时间长度r欲1/4户2048s;为防止风速模拟过程中发生混叠和失真,时间划分段数M应满足iV户2A,结合导线有限元时域分析的收敛性条件,在风^!=速模拟时取M32i65536r.125yf;相应的时间间隔J《003。2.4.3风速模拟结果留2-5给出了输电线路特征点处耐张段中间位置处的悬垂绝缘子挂点2-2(,对应图中的挂点4一)考虑脉动后的风速时程。值得提的是,为消除导线模型上突加荷载冲击放大效应的影响,在风速时程的前10化加入风速由0增至风速平均值的线性增长过程。50广|;||||^||500100015002000时间(S)图2-5挂点4处的风速时程图2-6为脉动风速功率谱的目标值和模拟值的对比。可W发现,风功率谱的目标值和横拟值吻合较好,说明获得的风速时程能够有效反应脉动风速能量在频域内的分布特征。’10—特征点白功宰谱fJ-^-、-2*'〇!〇1〇10!0H頻率(z)图2-6特征点处模拟功率谱与目标功率谱比较17 浙江大学硕壬学位论文第2章考虑气动陌尼效店的输电线路风偏动态计算2-7-图为按式口16得出的风场浦流度模拟值和理论値的对比。不难发现,风场浦流)度的模拟值与理论值吻合较好,说明采用修正后的Davenport谱模拟得出的风速可有效反映风场瑞流度沿高度的变化规律。■016.?-0.模m辅流强度15I-.5.巧巧()1平均风削面0-.14,货--013"强奴*..襲。-难穀狂心12..产节--X01,0.'〇20354041525305■莉度(B)图2-7模拟风场浦流度与理论瑞流剖面的比较2.5有限元结果分析取风荷载与导线夹°角为〇,将谐波叠加法模拟得到的脉动风荷载时程加载至有限元模ewmark-型的各个节点上,,采用无条件稳定的N目法对非线性动力方程直接积分求解并ew-运用NtonRaphson法对毎个时间步末尾的位移进行迭代。同时为考察气动阻尼的影,响,分别采用考虑和不考虑气动阻尼两种瞬态求解方式计算多档距输电线路的风偏响应。当采用考虑气动阻尼瞬态求解时一一荷载步时导,取前荷载步作用后导线的运动速度与后-5线所受风速的相对速度,代入式(2)进行计算。止其与橫担过近闪络-。图28给风偏主要是考虑绝缘子串的位移,防,发生放电出了位于线路中间的3、4、5号挂点的风偏水平位移时程计算结果。不考皮側R1;己不々巧气动阻旧I——-6,考阮动阻尼,考由气动阻巧6-55*鸭蜘师4-寒曇予-哪想麵補義-巧3径3I''22-票罢II11000400800120016002002400K0(J2001600200024抑41时间/S时间/s挂点3挂点518 浙江大学硕壬学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算不考虑.气动阻尼—考点与动阻尼6「占'\片2-U^4008001200化0020002400时间/s挂点4图2-8互处挂点的风偏水平响拉时程-在获得绝缘子串水平位移后,根据式(224即可得到风偏角)水平位移"I、。=//IW02-26sin()绝…缘子气山串长I,度一()"动助错邑90r不考虑側KIMQO重5i—考度初as—I…考這气动阻尼巧.15抑曲二:攀卿幽幽耐0唯叫鬟30.400800120016002000400800120016002000时间(s)时间(s)挂点3挂点5不考虑气动阻尼I—考底气动阻尼73.井剛城卿戰i巧30^15-U.?■??I■I■I-Q/400800120016002000时间(s)挂点4图2-9王处挂点的风偏角度时程19 浙江大学硕壬学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算--从困28、图29不难看出,,由于平均响应是由平均风荷载所产生因此是否考虑气动阻尼对于响应的平均值无影响,而气动阻尼对动态风偏响应幅值的降低作用十分显著。2-1因此,忽略气动阻尼的影响势必高估线路的风偏角。图0给出了各绝缘子串挂点处的风偏角极值,相应计算方法为:'-=+-a>^^(227q/)■式中均值一疗为,〇为均方根,是个保证因子,称为峰值因子,为化到要求的安全保证率98.61%,选取峰值因子2.2。——一60rI不考巧气动阻尼々巧气动阻记I■40-30^*12345677123456绝缘子挂点号绝缘子挂点号(a)各挂点风偏角平巧值(b)各挂点风偏角极大值图2-10各绝缘子挂点处风偏角统计值一-分析图210可知,方面,靠近耐张塔挂点处的绝缘子串风偏角比在整个耐张段中部挂点处风偏角要小,这是由于支座约宋效应的影响,使耐张段两端附近绝缘子的风偏角一一较小,这结论在现实情况下普遍存在,在考虑气动;另外方面还可看出阻尼后,风偏角极值大幅降低,说明气动阻尼对风偏响应有着较为显著的折减作用。因此在进行脉动风载下的动态风偏计算时有必耍考虑导线与来流的相对运动所引起的气动阻尼。2.6采用ASCENo.74所提供气动阻尼21AS.6.CENo.74中的气动阻尼么式美園输电线路设计规范ASCENo.74《GuidelinesforElectricalTransmissionLine工程师协会所出版一StructuralLoading》是由美国主木。作为世界上十分重要的本输电线路规范,其中的许多公式及参数取值被他围大量借鉴与引周。在ASCENo.74中根据Davenport的研究成果,给出了导线气动阻尼的计算方法:20 浙江大学硕壬学位论文第2章考虑气■动阻尼效应的输电线路风偏动态计算=-0.000048r/X/12X矣C22(〇九W))(巧/式中;F为导线有效高度处lOmin的平均风速,单位为feet/s处特别指出各参数〇,此的单位,是因为美国常用单位与我国不同,若有疏忽容易产生错误。d是导线的直径,单位为in沈,。为导线的体型系数,/w为平面外摆动的基础频率,单位Hz,采用下式计算得到;.5=5。-1/分口2九(种—式中:W为中跨的弧垂长度feet。各,单位为此么式基于输电导线风偏W阶平面外摆动为主的基本假设进行推导得到的。对于本章所研究的八跨输电线路,在设计风速=27m/s2-巧时,由公式(2%)可得此耐张段的气动阻尼为化08987。〇在对导线的非线性运动方程进行求解时,为便于分析常将其阻尼矩阵C视作质量矩阵M与刚度矩阵if的线性组合,表达为:=CaM+0K(2-30)y口式中,常系数,片是由结构总阻尼比《、前两阶基本频率份2所决定的。表示为:2-3(1)的+田12=^—-32)片《(2巧+"2在非线性运动方程-(26)的求解时,输电线自身结构阻尼《取化005;因此在采用,ASCEN0.74提供的气动阻尼计算导线风偏时,在非浅性运动方程的阻尼矩阵中直接考虑=+--1)气动阻尼的影响,根据3、32)。,即^是兵式(2式(2求得相应的常系数口,片通过在ANSYS中设置相应的〇,片系数巧可使得在ANSYS模型中的阻尼值为所求得的总阻尼。22ASGENo4.6.采用.7气动阻尼有限元结果根据上文所述采用ASCEN0.74中所提供的气动阻尼与结构阻尼之和作为有限元摸型的总阻尼,进行计算,将计算结果与本章所述考虑气动阻尼的方法迸行对比,如下图所示。21 浙江大学硕壬学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算本文义法80rASCENo.7470-厂"'3〇^|:7苗20-UA10j■■II■I.I./I.I4008001200160020002400时间/s2-1图1挂点4风偏角时程-ASCE图211给出了挂点4位置分别采用本章所述的考虑气动阻尼的方法与考虑L.74所给出的气动阻尼所计算得到的有限元时程结果。从图中可乂看出No,采用这两种方法所得的结果十分接近,整个时程过程基本吻合。2-22-27在图1中给出按式()所计算得到的300个导线节点的水平、竖直位移的极值统计值。-^^■本文义法*3旅1--ASCENo.74?)IASCRNo.74|?I期-i哀0?0000000巧000〇501001502002日03005011523节点巧号节点巧号(a)水平位移极值(b)竖直位移极值图2-12所有节点的风偏响应统计值由这两幅图中化可W看出,本文提出的考虑相对运动速度产生的气动阻尼与ASCE-28ASCENo.74给出的么式(2)计算得到的气动阻尼对于导线的作用十分接近,说明No.74中气动阻尼计算方法可供之后风偏研究时采用。2.7本章小结本文W500kV韦塔八档输电线路为研究对象,建立了考虑气动阻尼的导线风偏非线性动为方程。通过有限元软件,采用非线性瞬态分析方法求解了输电线路的动态风偏响应,22 浙江大学硕±学位论文第2章考虑气动阻尼效应巧翰电线路风偏动态计算并与考虑美国规范ASCEN0.74所提出的气动阻尼有限元结果进行对比。研究结果表明:(1)输电线路在受风荷载作用情况下时,导线与气流的相对运动会产生显著的气动一阻尼,在高风速下该阻尼般大于导线本身的结构阻尼当考虑此因素时能够显著降低,,导线风偏响应的脉动值。因此在进行风偏动态计算时,有必要考虑导线与来流的相对运动所引起的气动阻尼。(2)通过有限元软件模拟结果可知:在现实情况下,靠近耐张塔挂点处的绝缘子串风偏角会比耐张段中部挂点处风偏角要小,这是由于支座约来效应的影响,使耐张段两端附近绝缘子的风偏角较小。(3)ASCENo.74中提供了导线气动阻尼的计算方法,通过此算法算得的线动阻尼值,并在有限元软件中运用此结果,将此方法计算得到的风偏响应时程值、统计值与本章所提出的考虑气动阻尼的方法有限元结果进行对比,发现两种算法算得的风偏响应十分相近。巧明了ASCEN0.74所提供的简化计算方法结果较为准确,可供之后风偏研究时采用。23 浙江大学硕+学位论文第3幸我国规范风偏计算方法及与多国规范对比第3章我国规范风偏计算方法及与多国规范对比3.1引言风偏闪络的发生将会带来巨大的经济损失W及严重的社会影响。在输电线路设计过程中,从节约计算成本的翔度出发,工程师大多根据规范中的相关条款及规定,采用简化算法对输电线路的风致响应进行分析,并开展相关的设计工作。规范的准确性及可靠性对于架空输电线路能否安全运行起着举足轻重的作用。因此,深入研究我国规范很有必要。一?目前我国输电线路规范中存在两个较为明显的疑问:首先架空输电线路是个连续的多自由度非线性结构,在整个耐张段内各档导线之间存在着精合作用,而规范中风偏计算一个绝缘子挂点进行受力分析所采用的单摆樸型只孤立地对其中,无法考虑各档导线相互之间的影响。采用单摆模型计算得到的风偏角与实际的多跨模型计算结果的差别未有文献提及单摆模型的准确性仍有待研究一。另,方面,输电线路受到的风荷载包含明显的脉动效应,在进行输电塔强度设计时规范采用风荷载调整系数来考虑导线所受风荷载的动为效应影响.;在计算导线风偏时,规范规定风荷载调整系数取10大的因素。,即不考虑动力放w3ww58国内多位学者]从与国外规范比较的角度进行阐述,表明需要考虑动为效应的影响,却由于缺乏试验和数值模拟分析,未能给出相应的动力放大系数值。相对于其他一些发这国家。,我国属于风偏多发地区将我国与欧美国家的输电线路设计规程进行对比,巧究不同规范中具有相同含义的参数取值差别,有助于了解我国参数取值不足之处,确定风偏事故频发的原因。本章首先基于静力单摆模型,采用规范给出的导线风荷载计算么式及相关计算参数的取值方法,求解了输电线路的静态风偏响应,并与采用有限元分析法求解得到的八跨线路的平均风偏响应进行对比,考察了单摆模型的准确性。随后,将有限元软件计算得到的风偏总响应与平均响应进行比较,确定了动为放大系数的取值。然后将我国规范中导线风偏的风荷载计算方法与美国一、英国及国际通用规范进行对比,并将这些风荷载么式归纳为个算式。L,进行参数的对比分析最后X第二章中实际超高压八跨线路为例分别根据各国,规范计算各档导线的风荷载值,直观地表现出我国导线风荷载计算值与他国的差别。24 浙江大学硕±学位论文第3章我国规范风偏计算方法及与多国规范对比3.2风偏计算模型讨论31.2.弦多边模型在实际工程中,因风偏静力方法具有计算简便、原理易懂的优点而被广泛采用。常用的风偏静为计算模型主要有弦多边形法和刚体静力学方法。弦多边形法假定各个绝缘片是相互较接的刚体-,具体受力模式见图31.gv〇(\^图3-1弦多边悬垂绝缘子串风偏受为简图3-图1中悬垂绝缘子串末端是导线的挂点位置W,受到了导线水平向风荷载及竖向a?重力作用Wk集中作用。假设绝缘子串共分为n个刚体第?1个刚体长度为;I其竖直方,,,向所受的荷载为?gw,所受风荷载大小为斯,。这些荷载作用的位章均为各分段刚体部件的-中央-。根据绝缘子串受风荷载影响产生风偏后的平衡条件31),利用式(和(32)即可得到水平、垂直风偏响应4和Av。=^Asin,Z,0=1,乂狀+0.5+1各(3-,*扔,*1)。|^V■)-V眠+0.5各'+《+狀+0.5妍'+各""完"**,完^))25 浙;工大学硕+学位论文第3章我国规范风偏计算方法及与多国规范对比n义=乂c〇vZ,W=11'--1,\W+0.5各+-v,斗1各K(,玄32)。、、J-y合VIf^f^V+.W505++0.+J各震,*各W各*;[J^;J=S(3J)hi?.上式中:225P为空气密度,取L;4是第i个绝缘子的受风面积;V为第!个绝缘子钢片所受到的风速?。弦多边模型较好地模拟了现实中绝缘子串的变形情况,但也使设计过程增加计算量。在实际情况中,架空输电线路的绝缘子串所受风荷载及其自重相较导线而言是个小值。因此一,般来说,在设计过程中无需将绝缘子串精细化分为若干段进行计算。只有当绝缘子串的自重较大,且需要计算精确的风偏位置时,才会采用弦多边模型。32..2刚体直棒模型由于弦多边形法计算过程较为复杂,当绝缘子串较长、绝缘片数较多时所带来的计算量不容小觀,也给工程设计带来了西难。在实际工程中计算过程较为简单的刚体静力学方法应用范围更为广泛一。刚体静力学模型化称为单摆模型,是将悬挂的绝缘子串看做根受到均布荷载作用的刚体直杆。绝缘子串所受的水平及坚向为作用于刚体直杆的中部;将导线质量和所受的风荷载集中至绝缘子串的末端。并假定当风荷载与导线和绝缘子李自重达到平衡时得出的摆角即为风偏角。////■/么iIi,h-H图3-2输电线路风偏静力学单摆计算模型-20如图3:所示,绝缘子串风偏角可表示为26 浙江大学硕壬学位论文第3章巧国规范风偏计算方法及与多国规范对比-+/2|亞&=01孤()(3-4)W+G2、J式中:&和%分别为绝缘子串和导线所受的风荷载,Gv和Wv分别是绝缘子串和导线的自重。需要指出的是,计算导线风荷载时需使用绝缘子串所处巧塔的水平档距,而确定导线自重时则需采用垂直档距。对于导线所承受的风荷载,输电线路设计规程规定采用下式进行计算:;-Asin0(35)AeA化P:a为风压不均匀系数式中,根据表格插值得到给出了两种取值方法,分别用,我国于设计和验算。//sc为导线的体型系数当导线的线爸小于17mm或覆冰时不论线径大小().2线径大于或等于17mm时取1取1.1。。,Pc为风巧载调整系数,在风偏计算过程中不考虑d为导线的外径,L是直线塔绝缘子串所处输电塔的水平挡距,0为风向与线路走线之间p的夹角,w〇为基本风压。绝缘子串上的风荷载采用下式计算: ̄-(36)式中,今是绝缘子串受风面积。由上述介绍可W看出,刚体直棒模型受力情况明确,计算过程简单。因此在架空输电线路设计过程中受到广泛采用。3.2.3多自由度连续多跨摸型第二章中所建立的八跨导线模型是一个典型的多自由度连续多跨模型。此类模型通过有限元软件,根据现实中导线的情况,采用单元对导线进行精细化的模拟建模所获得。多白由度连续《跨模型是对现实导线的一模拟,棋型情况与实际情况的导线基本致,当其受力时,各跨导线模型之间存在着精合作用。在计算导线风偏过程中采用多自由度连续模型可较好地模拟实际情况下多跨输电线路,计算得到准确的风偏响应值。3.2.4单摆模型正确性讨论实际情况中输电线路多为连续多档,各相邻档之间存在相互作用和张为传递效应,按照多自由度连续多跨模型进行计算可W获得合理结果。传统的静力学单摆摸型如上文所介绍的刚体直棒模型法考虑整个耐张段"",无内各摆导线之间的拉扯作用,因此有义要通过27 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与#国规范对比有限元模型与单摆模型计算得到的风偏角对比,进行单摆模型的可靠性研究。为方便比较,确定选取不同计算模型对于计算结果的影响,在单摆模型与多自由度连续多跨模型上施加相同的静为荷载:,即-(1)按单摆模型计算时,根据规范所提供的风载公式(35),计算出导线风荷载,在计算过程中取不均匀系数a为1.0。一(2)在采用多自由度连续多跨模型计算时致的平均风荷载。,只施加与单摆模型相取10m高度处平均风速为27m/s,考虑风高系数。基于两种计算模型求得的各绝缘子串挂点处的风偏角如图3-3所示。-?一■多自由'度连缕検型一I单巧垣瑣巧-?/\k:/\,,,,,,I11t11,12Q1234567挂点编号图3-3多自由度连续多跨模型与单摆模型风偏角的比较由图3-3可W看出,采用单摆模型计算得到的风偏角大于多自由度连续多跨模型的风偏角,尤其是对于与终端塔相邻的挂点,单摆模型计算值偏大较多。这是由于多自由度连续多跨模型中两端的耐张塔对输电线路面外摆动产生明显的"拉扯’’作用和各按线路间的相互稱合作用,使得其响应值比单摆模型要小。这表明,就模型本身来说,采用单摆模型来求解处于整个耐张段中的绝缘子风偏角是偏于安全的,可L乂适用于工程设计。3.3动为族大系数取值风荷载包含平均风荷载与脉动风荷载。平均风荷载作用于结构时,体现静为的效应,而当导线受到脉动风荷载作用时,所发生的脉动响应体现出动力特性。我国规范在利用式-(35)计算风偏时导线上的风荷载,取风荷载调整系数Pc为1.0,即不考虑导线的脉动效应。这样的简化方法是否可行,有待研究。取第二章中所计算得到的绝缘子串挂点处风偏响应极大值与平均值做比较,就可W计28 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比-算出动力放大系数。将各挂点处风偏响应统计值列于下表31所示:表3-尼时各挂点风偏响应的统计值1考虑气动阻挂点号1234567。34..平均值/.046.255847.450.455635.24.2均方根4.24.34.64.43.5.24。.357极值/43.255.565.560.064.842.8.171角度放大系数1.271.201.211.191.161.22风荷载放大系数1391.31.391.481.441.431.451-表31中给出了所有挂点在瞬态风场下,考虑气动脏尼所得到的风偏角度统计值。在 ̄脉动风场下求得的风偏角约为静态风偏角的1.151.3倍。将得到的风偏角转换为单摆模型进行计算,即可得到对应的水平风荷载极值。可W算得,在考虑脉动后水平风荷载极值约 ̄为平均风荷载的1.31.5倍。说明在实际风偏计算过程中脉动风的动力放大作用不容忽视,规范所取的动力放大系数偏于不安全。3.4多国规范线路风偏风箭载的对比我面输电线路在匙风情况下风偏闪络灾害频频发生,说明规范中风荷载取值方法存在一一定的缺陷。如何安全合理的对线路风荷载进行计算成为个值得研究的问题。对国内外有代表性的规范进行对比有利于学习到国外规范中先进的部分,从而改进我国的输电线路设计规范中有关参数的取值。本章选取的有代表性的输电线路设计规范主要有:中华人民共和国住房和械乡建设部-发行的《n〇KV-750KV架空输电线路设计规范GB5054520》50545(10)(下文简称GB);美国主木工程协会发行的美国输电线路规范《ASCENo.74》(下文简称ASCE74);国际电P91王协会输电线路规范《圧C60826-2003》下文简称旧C60826木工程和建筑();英国上aticeower ̄ar结构标准委员会发布的《Lttsandmaststi:codeofracticeforloadin》下ppg(w文简称t]BS8100。)3.4.1输电线路风荷载计算公式首先介绍各国的风荷载计算么式-。我国规范中风荷载计算按前文所述式(35)进行计算,在此不再赛述。29 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比ASCE74中,风荷载取值可由下式求得;iF=cos-VGCA(37)rwQK:K:,?Wf式中:F为风方向的风荷载;Yw为重现期巧载调整因子;Q为空气密度常数,按照不同的温度与压强情况下进行选取;Kz为风压高度变化系数;Ka为地形影响因子;V为50年重现期3s阵风风速;Gw为阵风响应因子;V为风向与垂直导地线的法线之间的夹角;C一=A为水平向投影面积f为风力系数体型系数般情况均取Cl。(),f;圧C60826中,风荷载取值可由下式求得:】G加Q3-8)乂典(式中:是荷载保证因子;Cxc是导线的体型系数通常取1;Gc是考虑了导线高度、r,地貌情况与脉动影响的系数、&是考虑了导线长度的档距系数;d是导线直卷;L为导线弦长;q〇为基本风压,可表示为:(-93)上式中:是空气密度,X是空气密度修正因子,K,通过查表求得。是基本风速r为粗粮修正因子,可查表求得。BS8100中,风巧载计算时分为平均风荷载与脉动风荷载分别求解:=R(3-10)再wf巧c〇v(3-11)巧成马W;及=sn-C〇王Ci(3口)CWcc皆。=氏义戸-1〇K。()(33)nd*音式中:马W是平均风荷载;巧W是脉动风荷载;是导线平均高度处的风速;Cc是阻为系数,由雷诺数和导线类型、表面粗稚度、覆冰情况等决定;化是导线直径;k是导线的弦长;Gc是导线阵风响应因子;乂A是地面粗植系数;为风向系数;a是高度相关系数;Ag为有效高度;是荷载因子;为基本风速。W上所述的公式^均1^准定常假定作为基础具有相似的形式将上述么式归纳为,,可一个算式:^^=-FxpvCA?sinw(314)jj^30 浙江大学硕去学位论文第3章巧国规范风偏计算方法及与多国规范对比式中:是各国规范中考虑了风压不均匀性及动力放大因素的综合参数,主要受到高^度。、跨长等影响?为荷载系数,P是空气密度,V为平坦地貌下或山区地貌下的风速,受基本风速及风剖面、山区地貌等因素的影响,Q是体型系数。A是受风面积,W是导线与风向间夹角。当采用不同规范来计算同一耐张段输电线路时所取W是定值、,风荷载计,P算值的差异产生主要是由于V、Cd、n四个数值的取值不同。接下来着重讨论这些参数的取值差别。3.4.2荷载系数尸荷载系数是根据线路的重要程度而确定,用来调整线路安全等级的。国外规范中荷载系数的取值差别直接体现在设计风速的重现期上。我围规范并没有荷载系数的概念,而是在计算荷载设计值时采用结构重要性系数与风荷载分项系数的乘积来代替荷载系数。表3-2各国荷载系数;重现期/年ASCE74IEC60826BS8100GB50545i^V根据图表差蘇'100-1."值选取,其中zih-1501.1对于特另3M黑府风荷载取200-1.3重*要的线〇96 ̄400-.1.25。1.45J^500 ̄y3.4.3风速的取值风速值的大小与基本风速、所处高度、所在地貌类型等内容密切相关。其中基本风速是根据不同高度的测风仪器一、时次时距所测得的最大风速,按定的高度、时距要求统,一50545826BS换算而成的最大风速。ASCE74;GB和圧C60;8100规定基本风速取值分别按3s、lOmin、化时距,平坦开阔地貌离地10m高,重现期按50年确定。由于自然界大气边界层内的风速有脉动特性一般来说速的时距与获得的平均风速成反比,,所取风关系。时距越小,计算得到的平均风速就越大。取不同时距时,基本风速所得值的关系可=W:4::表示为圧C:出1.310.94。在大气边界层内,由于地表粗禮元的摩擦作用,风速大小随离地高度的变化而发生变化:离地高度越高。,所受到的摩擦作用减小,平均风速相应增加31 浙江大学硕壬学位论文第3章我国规范两偏计算方法及与多国规范对比风速的垂直变化不仅受到地面粗粮度的影响,而且逐会受大气稳定度的影响。因此无论利用理论进行推导或者通过实际观测结果进行分析,想获得其变化规律都是十分复杂的。目前,主要用对数律和指数律来描述。指数规律和对数规律差别并不大,所目前围内外都倾向于用计算比较简单的指数曲线来表示风速沿高度的变化规律。GB50545、ASCE74、BS8100中都是直接使用幕指函数来表这风速的垂直变化规律。GB50545中:指数型风速剖面可表示为:(ZXU二-。U(315):、寸式中,a为平均风剖面指数,[/、分别为离地高度Z、Z处的平均风速。根提,,,,-粗粮度指数、梯度风高度W及风压间的关系,由式(316)可得平坦地貌下A、B、C、D巧类地貌对应的风压高度变化系数,其中/分别取1.284、1.000、0.544、0.%2a取化口、/?;。0.150.220.30。、、(ZT=_-/i-从(316)°U〇JGB50545四类地貌下风高系数随高度的变化规律如图3-4所示。3.5—A类[B类-3—.0C类0?5.I.III.I...I..I50100150200250300350Z(ra)图3>4GB50545的风高系数取值图ASCE74中,风高系数表示为:扣。—=-ii:21(.0)(317)vZgASCE74中对输电线路建设的区域地貌分为B、C、D王类。其中C类地貌代表着开阔地貌.,与我国B类地貌相类似。对应于B、C、DH类地貌:a分别取7、9.5、115梯;度高度Zg分别取365.76m、274.32m、213.36m。Z的取值范围要求:10含zSZ。g32 浙江大学硕去学位论文第3章我围规范风偏计算方法及与多国规范对比3.0「■—B类。2.二.51達^1.0.....■II.tIII.t50100150200250300350Z(m)3-5图ASCE74的风高系数Kv取值图BS8100中,风高系数表示为:。=-zh0(3-心如()/118)Ac)上式中:J苗是地面粗粮系数,&为风向系数,a是高度相关系数,七为有效高度。BS8100中共将地貌分为I、II、、IV、V五种类型与我国B类地貌。ra,其中第m类地貌等同在不同地貌下:,相关参数的取值如下表所示表3-3BS8100中不同地貌下参数取值地貌类型地面粗植系数高度相关系数有效高度KarhMgI1.200.1250’II1.100.140m1.000.1650IV0.86化192V0.720.2310当不考虑风向的影响时:,可画出风高系数如下图所示33 浙江大学硕+学位论文第3章我国规范风偏计算方法及与多国规范对比25—.I类「I—"类—…类^—IV类-2.0|V类I^.11.1.1.11.1...50100150200250300350Z/m图3-6BS8100的风高系数^1^取值图圧C60826中并未直接给出响应的风压高度变化系数么式,而是给出图来进行差值选取:,Gc包含了高度风压系数与风的脉动影响两个方面的因素3Q ̄".111ITTTI11nfnrn「nnrr11111111111111iII1![「ni[二3之二;■一二二二二D--二-:=Z:::Z:Z::Z:::::::^^::::Z:Z::::Z:ZE:::^2a----::^:一-?---—__ ̄——_——.C——二—"’一—二———二,ir1^:1IIIII三三=E::_IIIIiiirI=二乙---=二--_—==:二:二二4::记_二-_一:三!::三:二-一:.一^三二一^二二二二二lit::::::;'9—'-''二IIIII二III二IIIII二二IlIIIIIIIIII!-IN]IJ]IIIIIJiIJIiJUNNIIiNNIiIIIiIIlLLL1IIUI.610IS2025303SAO45505560hleilYstx>erorxigtivgLjm图3-7lEC60826中G取值图e从图3-4至图3-7可W看出,平埋地貌下各场地类型对应的风高系数均随高度而增加;同一高度下,随着地表粗稚度的增加,风高系数逐渐碱小,所有的规范所确定的风高系数都拥有相同的趋势。若取各规范所规定的标准地貌,即类比于我国规范中B类地貌所对应的地貌,将其风高系数作图:,可表示为34 浙江大学硕±学位论义第3章我国规范风偏计算方法及与多国规范对比"*CB5054528ASCE74^■—BS81(坚J-L2.4^畫\.2^0-.8'1i.i.i.i.1.1.1.1.1,1,125507日1001251501752002252日0275Z/m图3-8标准地貌下风高系数对比由上图可见:当高度较低时,BS8100所求得的风压高度系数较大;随着高度的不断,GB50545所给出的风高系数将大于其余两国增加;相对其余两国,我国风高系数对于’高度改变更为敏感。风速值除基本风速、平坦地貌下风高系数的影响因素外还受到了山区地貌环境的影响。我国荷载规范中利用系数对风速进行修正:1当所取位置为顶部B时:)'=-—-;71+tana(\(319)b[)]l.bH式中ana一ana,t为迎风面侧的坡度,若t大于化3时,取化3;K值对山峰取2.2,14H为山对山坡取.;顶或山坡全高(m);Z是建筑物计算位置离建筑物地面的高度(m);当>=Z2.5//时,取Z2.5//。2-其他部位的修正系数,.,),可按图39所示取A、C处的修正系数为10AB间和BC间的修正系数按线性插值确定。d44.^'4—]、fI+图3-9山峰和山谷示意图ASCE中取地形影响因子也t对山区风荷载直接进行修正,也t的计算公式如下:2K=(-(1+/:/:尤)320):,|23-式中:是考虑了地形特征与最大增速的影响值方法见表34。,取=-A:l;c//L山顶的上下风向距离增加而导致风速改变的系数。,((/J)是考虑离、||35 浙江大学硕壬学位论文第3辛巧国规范风偏计算方法及与多国规范对化心是考虑在局部地形情况下,风速随高度增加而减小的系数;JC是山顶与结构之间的距离;Z为建筑所在地面W上的高度。"":凤纖:淵倘,^h!1。胃而J上&上i當所戸^皆VIj1j%9SWWW巧J风妨醒巧是。纖汾;無a悬崖情况b二维山脊或王维轴对称山丘(())3-困10山区风影响示意图表3 ̄4A参数取值表—K/H/L,(^片山丘形状地貌类别y北?CTSIBCD二1.301..维山脊.4515531.515二维悬崖0.750.850.952.51.54王维对称山.051.1541.51.5&0.951英围规范中考虑山区风的方法是根据输电线路所处山区的位置:山顶、山坡迎风面、(3-2-(-山坡背风面,分别按1)、(322)、323)计算:。-"=己也-巧邸()(321)ZrW"10=严(3-22)^吃《(皆=v^zK-y(323)‘(y10-11其中://与均是地貌指数,/与顺风坡度、地貌分类有关由图3确定。而与/3-=山的高度和地貌类别有关由图12确定。必1);A*hh36 浙江大学巧±学位论文第3牵巧国规范风偏计算方法及与多国规范对比舊曇聽藍聽012I.5t7II,l?’,HHWn,???0AiiiTjMJMI1i1.ViVl\iA,3k-3-n图//的取值图表图312的取值图表从上述介绍可W发现,美国、英画规范中考虑山区地貌影响较为详细。我国规范采用一山区地貌影响因亲较简便的方式,所考虑的,通过个么式来进行计算,参数也较为简单少。3.4.4体型系数CdASCE、旧C中体型系数Cd直接可取1.0。BS中Cd取值较为详细,考虑了雷诺数和导一化6?1.3。线类型、表面粗植度、覆冰情况等因素的影响。常用的导线其值般取我国导线体型系数根据导线直径所确定:当导线的线径小于17mm或覆冰时不论线径大小)取(1.2,线径大于或等于17mm时取1.1。可见除英国规范外其余国家关于体型系数的取值差别并不大。3..45。系数对比实际风荷载对导线作用过程中具有两个特点:上的随机性导线不同位置所受到的风荷载值一1.风对于导线的影响具有空间和时间,般不同;一2.风包含着脉动的特性,对结构的影响般需根据复杂的动力学分析确定,而在实际工程中常采用拟静为的方法来考虑脉动风的动为放大作用。-式(314)中系数。;7即是考虑了这两个因素的综合参数美国规范是采用阵风响应因子G:w来代表系数的含义+2.1EJ-馬C324)37 浙江大学硕壬学位论文第3章巧国规范风偏计算方法及与《国规范对比丄。?=—-左4.9V^(325)()及二。-W’"^7—1+4:5为导线和地线的水平档距式中;&为电线有效高度单位为ft为所有导线挂,,即点的平均高度;为10m高3s阵风风速与lOmin平均风速比值,其值为1.43;OI/W为持续风的幕指数A为表面阻力系数L为浦流积分尺度这兰个参数的取值见表3-5。;;,,表3-5阵风响应因子参数取值地貌类劝持续风的幕指数aw表面阻力系数k浦流积分尺度B4.50.01051.8C7.00.00567.1D^0.0037^ASCE74中阵风响应因子、适用于任何电压等级输电线路,并且该系数与导线的高度控距和地貌条件等因素有关,Gw包含了3s与lOmin的风速转换关系W及风的脉动影响。因此剔除风速转换参数后就可レ乂得到lOmin风速下的;系数。7在C类地貌下-(a00mZ ̄50m(b),图313)取档距S为困定值3,h取10m;中有效50m200m ̄400m高度Zh分别取固定值10m、,档距S取。从图中可W看出求得的/系数/随导线有效高度、档距长度的增加而降低。i-6r61.「-二1.55Z10mh5■.訂.和145、氏'^;IIas'、旨1.3S,..5;1.21020304050200250300350400Z/mS/m,,(a)档距不变(b)有效高度不变图3-13ASCE计算得到的阵风响应因子GwBS=在8100中电线路脉动风荷载表示为巧G中G是导线阵风响应因,输we巧W,其c子,与ASCE74不同的是,BS8100中的阵风因子只表示为脉动部分与平均响应的比值,+G因此0c)代表了系数;的含义。Gc由下式可W求得:738 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比-G广KJC(327)i式中&是导线的长度因子,用来考虑随机风速沿线路方向的空间相关性对导线风荷载的影响,,&是高度因子,考虑了导线平均高度处的奈流度对风荷载的影响这两个值分--14315别按图3、图插值选取。—_支电古占式山—一广—r巧-———————-——0.4frjII!軍?_j—…_——,—了---厂—广―—__—-了卞十下T寸丫…_----下-丄…tl十十十了I??--?-’去史山山1山山山击山Lmt}-图314长度因子?巧JM/W巧扣》"ii:iirr^1TTH"—'^- ̄ ̄-M!irxr!4I-斗灌邦糊寻Ff ̄ ̄ ̄ ̄ ̄h ̄ ̄!iiIIjIIpI\ ̄i ̄ ̄?nS^吏.tJ10去占Iw4)4o^iono300化!■>图3-15高度因子&表3-6给出了在m类地貌下有效高度分别为?400m变化时10m与50m,档距从200m,(1+&)的取值。可L乂看出1+Ge与ASCE求得的Gw有相同的变化规律,均是随着导线有效高度、档距长度的增加而降低,但其计算结果较ASCE要大。3-6(+Gc表按化计算得到的风荷载调整系数1)距^'有效高'200m250m300m350m400m!^^10m1.951.891.851.821.7950m1.811.761.721.的1.66圧C标准中G由图3-7来确定考虑了高度风压系数与地形和高度有关和风的脉动c<)()影响两个方面的因素,因此将Gc除W我国规范计算得到的就能够把风压高度系数分39 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比L;离出来,再乘乂其档距因子&,即得到/系数。表3-7旧C所取的C系数供'有效高^^200m250m300m350m400mm1.671.601.5110.751.72150m1.511.491.441.381.30我国规范中风荷载调整系数化是针对500kV和パ0kV导地线考虑脉动风荷载的放大系数,规定在风偏角计算中不考虑该系数的影响。不均勾系数a考虑了水平控距变化对整〇3-4个耐张段风荷载影响,片与a的乘积就相当于式(1)中;7系数。表3-8GB所取的cxa;/系数P档距200m250m300m350m400mxa08.74.70.67065Pc0.0.ASCElOm有症亩志2-■ASCE50.〇m宵效嵩巧1.一广BSlOro宵效商度■—^—-A?BS50m—A有效風度1?8??---二—化ClOm任巧商度▲—1EC日Om冉效商度,LbC--—?GB50545、I—心-?1■-:.4T亡窥'-^1.2-1.0^0.8^I.■?.T■.06IIIIII■2002加300350400450500档距(S)图3-61各国规范的/系数对比/将各国规范算得的系数作于图3-16上。从图中可W看出各国规范关于7系数的取巧,值趋势均是随着接距的增加而减小,除我国,其他国家规范中都随着有效高度的增加而降低。BS所取值最大,而我国规范由于未考虑风荷载调整系数,使得C系数值最小。346八跨模型按多国规范计算风荷载值..为直观地体现按多国规范所计算出的风荷载值大小差别二章2-2中所示的八,取第图跨输电线路为计算模型:,分别按照各国规范进行求解,其中的计算参数取值如下所示40 浙江大学巧去学位论义第3章我国规范风偏计算方法及与多国规范对比(1美国规范《ASCE74》:是50年重现期对应的重现期前载调整因子取1.0;基)本风速按照3s阵风风速选取,是lOmin平均风速的1.43倍,取38.61m/s;&为风压高度变化系数;地形影响因子&,,线路所在地貌对应于美国规范中的C类地貌,由于不考虑山区影响,地形影响因子取1.0;阵风响应因子Gw是根据跨长、挂点高度、持续风的幕指数表面阻为系数乂及浦流积分尺度右确定,本耐张段各挂点对应取值范围在化64-0.67之间ASCE741。;《》中导线体型系数。取值得注意的是,当取国际单位进行计算时,Q取化613,这样获得的荷载值单位为N。(2)英国规范《BS8100》在计算具体某高度处风速巧勾时,根据线路所在地貌环境特征.1.21,风速安全保证因子取1,相当于将基本风压放大了1倍段各档导线的档;耐张-距因子取值&在0.670.77之间地貌特点,地表粗複度影响系数取1.0;根据线路所在;英国基本风速是选取化时距,相应平担开阔地貌离10m高处、重现期为50年的风速是lOmin时距风速的化94倍,故在计算中基本风速取25.38m/s;《BS8100》中规定导线体型J-400取C系数应根据雷诺数来进行计算,根据LG/35的物理特征,c为化9。(3)国际电气规范《圧C60826》中:输电线路体型系、数取1.0,10m高度处基本风速 ̄0取27m/s.98之间G;耐张段各越导线的摇距因子取值&在化%;c是考虑了导线高度与取值范围在2 ̄2地貌情况的风载系数对于本课题所研究的耐张段.2.38;,,其之间考虑地貌影响的系数&取1.0。(4)我国规范《GB50545》中:基本风速取27m/s,不均匀系数按摆距进行插值,取 ̄1值范围在化630.72之间的体型系数.1虑化动效应荷;导线/%取1;计算风偏时不考,风载调整系数A巧1。图3-17是取W上设计参数,并利用各国规范公式所计算得到的导线上的风荷载值。可W看出,按照多国输电线路风巧载规范计算得到的风荷载在耐张段范围内的变化趋势基同,这是由于上述国家的输电线路风荷载么式均是采用了准定常假定,且在计算设计风速时都使用了指数率的风剖面。但是由于各种参数取值的不同,使得按照不同国家输电线路一定程度的差异设计规范计算得到的风祗载值存在。41 浙江大学硕壬学位论文第3章我国规范风偏计算方法又与多国规范对化-*—SCE74IA55-^BSSIOO「■-^lEC-50-^GB50545|[.....I.I?2QItIIII1234567绝缘子挂点号图3-口规范风荷载的比较美国规范在计算风荷载时虽取3s阵风风速,但其在阵风响应因子Gw中考虑了风速的3slO转换,在计算Gw时除去阵风风速与min平均风速比值斬的平方,按ASCE计算出的风前载值相当于取lOmin平均风速,因此其计算值并非最大。中国、美国、英国和国际1关于体型系数的取值分别为1.1.>(、1.0、0.9和10,可^看出各国所取导线体型系数相差不大,不是导致风荷载差异较大的主要原因。由前文3.4.4中^系数对比可知,BS8100所取的W系数是最大的为IEC、再为ASCE最后是GB50545。这也是造成风荷载大小,其次,差距的主要因素。相比之下按照我国输电线路规范计算得到的风荷载值最小,这是导致我国输电线路风偏设计偏于不安全、闪络事故时有发生的直接原因。3.5本章小结本章着重讨论了我国规范中关于风偏的计算方法,分别从风偏计算模型、动为放大系数取值两方面进行研究,并将我国规范与若干发达国家规范进行对比,得到W下主要结论:""(1)由于多跨导线的拉扯作用,使得其在施加相同风荷载时计算出的风偏响应比单摆模型计算值要小,这表明单摆模型是偏于安全的,可W在工程计算中采用。(2)对比有限元计算模型得到的风偏响应极大值与平均值,可W发现导线受脉动风 ̄ ̄荷载作用时所得风偏角放大约1.151.3倍对应的的水平风荷载放大1.31.5倍。,,我国在风偏计算时取风荷载调整系数为1.0偏于不安全。(3)对多国规范进行了分析和参数对比,并采用八跨模型来计算风荷载,从中可W42 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比得知:我国所取的风速、体系系数与其他国家相差不大,但关于风压不均匀性及动力放大两个因素的综合考虑与他围相比较小,使得我国风荷载计算值偏小,这是导致我国风偏闪络现象频发的主要原因。43 浙江大学硕壬学位论文第4章利用巧域法计算导线风偏第4章利用频域法计算导线风偏4.1引言目前风致响应求解的主要方法有时域法与频域法。时域法的基本思路是将计算时间段一定顺序在所求结构上进行施加巧作用在结构上的随机荷载分为多个荷载步,利用迭,按一一代巧方法,W前荷载步作用后结构响应结果作为后荷载步开始计算的初始条件,求解一个荷载步后的结构响应r-施加毎。常用的时域方法有iVewmak片法、线性叠加法、■-0法-W沁on、Runekut法及Houbolt法等时域过程计算适用于各种形式的结构g,可W考虑结构所有的非线性因素,且能获得响应在整个时程过程中的表现,表达直观,易于理解,且计算值准确,但计算过程巧杂,消耗时间长。频域法常用的方法有模态叠加法、虚揪激励法等等。当采用模态叠加法计算时,先利用广义模态主坐标将结构的基本动为方程进行转换,再根据结构不同模态振型之间的与质量矩阵和刚度矩阵的正交性对动力方程进行解賴,就能得到各阶模态主坐标的振动方程,将动力方程进行傅里叶变化将其转为频域,通过引入频响函数对目标风致响应进行求解。在频域算法时忽略各模态之间的稻合效应时(SquareRootofthesumofSquares),会使得计算量较少,但这样计算的精度无法保证。若想得到精确解则需要考虑各阶模态之间的相关性(ComleteuadraticCombination)pQ,采用这样的计算方法有较大的计算量。为了在保证求解精度的前提下提高计算效率、减小工作量林家浩亚辉提出了虚拟激励法(PseudoExcitationMethod),这种方法,,张。在计算中自动包含各振型相关性,具有理想的计算效率总体来说频域方法相对时域法具有简单、快捷的特点,在计算高层、大跨等结构的风致响应计算中已得到了广泛的应7说的][]【]肿。在频域范围内采用模态叠加法计算响应只适用于线性结构。导线受风荷载影响产生风偏这个过程表现出明显的非线性特征,若直接采用频域法进行风偏值计算将出现问题,因此国内外采用频域法来计算导线风偏响应的研究还鲜有涉及。本章首先阐述频域法的基本原理,分析非线性结构在采用频域法计算响应时可视作线性的基本条件及假设,考察了风偏响应中平均响应占总响应的比率,并W此提出可W利用频域法来计算导线风偏的方法。将频域计算结果与时域结果进行对比,发现采用频域法的计算结果与时域法结果相近,证明利用本文方法来进行导线风偏的频域计算是可行的。同时为提高计算效率取10阶、5阶、3阶、2阶模态作为参与模态进行组合结果表明:,再,44 浙江大学巧壬学位论文第4章利用频域法计算导线风偏当取2阶模患采用SRSS法进行组合时具有较高的计算精度,且计算效率最高。4.2基于频域算法的导线风偏响应基本理论频域算法具有计算效率高,计算工作量小的优点。利用有限元软件对结构进行模态分析得到自振频率、,、质量各模态下质点位移等模态信息后可W直接通过模态叠加法对结构的动力方程进行解稍,从而获得模态主坐标下的动为方程。4.2.1频域法计算导线风偏响应的过程导线在受到风荷载影响产生风偏响应时,,主要产生出导线平面方向(y轴方向)位移W及竖直方向(Z轴方向)的位移,其沿导线方向(X轴方向)及自身扭转(6角方向)位一个连续均布质量的非线性结构移较小,在进行风偏响,故在风偏计算时不考虑。导线是应分析时可将其简化:将导线的均布质量及其所受荷载集中于这段导线的中点,从而形成一系列质点。风偏过程中输电导线的简化模型W及坐标轴设定如下图所示:A%图4.1风偏示意图及坐标轴设定基于随机振动理论:,导线的简化模型在发生风偏响应时的动为方程可W表示为M。CCKU"。+=_(4。^[化]巧[叫巧z巧〇mzCzKKfLJ[J,J[J[y,,J[J[J==上式中:YfZ.Z...、Zz..分别代表当t时刻导线各[乂佩扔的,,[(),(〇,,w(〇r抑(of,2质点在Y轴与Z轴方向位移的时程向量;Y、么代表导线在两轴方向的速度;乂、么代=d?ia表加速度;W是离散化的导线质点个数;Mg[/?代表导线的质量矩阵为第I段导;],=F=FF线跨中集中质量..f;F的F(〇...F(〇f[的.分别为各质点y,(0,,w(〇f巧,,y2,,,w,z,z2、沿y,z方向的水平风荷载和竖向风荷载;K,K,代表结构的刚度矩阵,zzu45 浙江大学硕击学位论义第4章利用频域法计算导线风巧均为AxTTA矩阵。导线在仅受到水平方向(沿Y轴方向)荷载时,不仅会发生沿荷载方向的位移,而且还会发生竖直方向(沿Z轴方向)的位移,表现出明显的賴合性质。因此在刚度矩阵中Ku、表示考虑了賴合作用。Cw,C^,Cu,是结构的阻尼矩阵,也均是考虑了精合作用的AWV矩阵。若所分析的结构可W视为线性:,则可将结构的变形表示为所有模态的线性组合Y=++...+t+...=-2t如t+tt(4)()li()知12(〇托q。()Mw()知()2]=nl针对第i个质点而言其运动可W表示为:,=...-++++...+Ytttt!>t(43)i()^qi()^q2()^qn()<wqw()上式中:i个质点在第n个模态上的位置量知代表第;gO为导线阶广义模态主坐/_/标时程t.34.1进行,Y是第i个质点在t时刻位移量。将式4代入式解賴,若阻尼矩阵i()为民ayleigh阻尼:=CaM+0K(4.4)y式中;a、P分别是假设的与质量和刚度的比例系数。则可利用不同模态与质量矩阵和刚度矩阵的正交性质可W得到:,=兰-m+.’心c+巧0的(45)成)A的W0、义作,}{},==式中:W.+是结构阶广义质量C.4;?()7,/为/阶广义阻尼,j完巧媛巧媛j戸冷_=112=it2;rw是结构的第阶广义刚度(/),和分别是结构的第阶自振频率和阻尼比;,}_/《y,,是第阶结构模态向量一,是个八W炬阵,矩阵中毎个元素代表在j阶模态中质点i_/一^x.=的y、z向位移;个Wl矩阵;的^的}为第{的}是&巧_/阶广义荷载时程,可{}{W由导线上的风荷载时程向模态向量投影得到。对时程运动方程进行傅里叶变换;,就可W将其转换成频域范围内频谱之间的关系=(4-6))K(/)(/)Wf\看式中;与是分别是由傅里叶变换得到的导线第阶模态的位移响应功率&//)y^!乂及广义风荷载的功率谱密度谱密度;好是导线结构第j阶机械导纳函数^/,按照下式)计算:46 浙江大学巧击学位论文第4章利用频域法计算导线风偏=//(4-7)(/)J,^2-1/+//'(/乂)/%,由于激励谱是双边谱,结构第阶广义位移和加速度可由功率谱在频域范围内_/积分得到:=WK(4-8)<而(/W!而軒4=巧-乃乃斯(49)而加而和防钟V)当考虑各阶模态之间的相关性来计算第/个导线离散自由度的位移响应均方值和加速leteuadraicombination度响应均方值时,采用CQComptC)法可W得到精确的解,表示(CQ为:=+4-芝玄4知W八(10)屋芝粉勺=^71k\J*k==口。+。s4-z(11)1么倫i艺夺啦扣而机(■y,)l=^7\1j*kj式中;n为选取的参与模态数;是广义位移相关系数;为广义加速度相关系数,表达如下:CR(/)骑鸟(/)别g々机々r-(412)片、。9八w)K(/(乃心/)皆衣lW邮CRe[(2死/T好,(/倘(/馬(/)斯]s,:*J7;(413)*、。2。。4425八批/)(/)么(/)命獻/)(/)(/)新片,防|姑在|在‘式中;//,(/)表示导线k阶机械导纳函数的共扼复数;S&(/)代表第j阶与第k阶广义风巧载互谱。若忽略模态间的相关性,可W提高计算效李,即采用SRSSSuareRootof化eSumof(qSquares)法来计算,可表示为:口=口s=-护z(414)\么\(y,)恥47 浙江大学硕壬学位论文第4章利用频域法计算导线风偏22==Z(_。41^垃知诉窗y,)■H422..导线使用频域算法的条件讨论-上述内容从式(42)开始都需要结构满足线性要求才能继续进行计算。当采用频域法求解非线性结构风致响应时,并假设其发生小位移,若结构在受荷作用下只发生小位移twi过程中结构刚度、阻尼矩阵不发生改变,即可将模型视作线性结构。导线受风荷载产生风偏过程中位移值大,因此无法直接采用频域法进行计算。风荷载由平均风荷载和脉动风荷载组成一般来说平均风荷载会比脉动风荷载大许多。导线在受,_到平均风荷载作用后产生的风偏位移占总位移的绝大部分,因此可政认为仅受到脉动风荷载影响的导线风偏只是一个小位移的过程二。取第章中所提到的八跨输电线路有限元模型-2绝缘子挂点处时程结果。对其取平均值与最大值,将结果列于图4。5.2.1二嚮滅!LI,[1/\1^/^?.I.I.1■I.I.I.I.I.I.I.I,I,I.2*612345675671234绝缘子挂点号绝缘子挂点号(a)平均风荷载引起的y向风偏值(b)总风荷载所引起的y向风偏值图4-2风偏值的对比-2-(图4(a)给出总响应的平均值,就是平均风所产生的平均响应。图42b)中给出的是各绝缘子串挂点处的风偏位移极值,是取98.61%安全保证李所得到的。将平均风荷-载所引起的风偏值占总风偏值的百分比数据列于表41所示:48 浙江大学硕壬学位论文第4章利用频域法计算导线风偏4-表1平均风引起的风偏值占总风偏值比例 ̄ ̄ ̄ ̄ ̄ ̄挂点编平均风总风偏值平均风引起的号引起风/m风偏占总风偏偏值/m值比例/%2.753.4581挂点23.574.化86挂点34.094.6090挂点43.6342686,挂点53804.3888.挂点64.084.579072.843.4284由此可见,平均风荷载所引起的响应占总响应的绝大部分。在现实情况下(即考虑气动阻尼情况)5--,平均响应约占8%90%,相应的脉动响应只占10%15%。若将平均响应和脉动响应分开考虑,可认为脉动响应时导线发生的是小位移,化动响应部分可W采用频域法进行计算。即可先在导线上施加平均风荷载,计算得到平均响应,再W施加平均风荷载后的导线模型作为新的计算模型,使用频域法计算出脉动响应,两者相加可得到总响应值。4.3单跨导线算例41.3.单跨导线的有限元建模?|備自v曲I*巧巧化抽Saty<S3m^p^麵///////////////////////////////图4-3单跨导线的形式49 浙江大学硕壬学位论文第4章利用频域法计算导线风偏-本章W图43所示的单跨导线作为研究对象。此导线跨长为353m,挂点高度38m。222549N8X0D.ASCR面积480.8mm。,年水平运行张力,导线型号为,单根截为便于分析八分裂导线等效为一建模时根据导线结构参数的特征。,将根考虑到输电线路仅能承受拉力而不能承受压为的特点,采用LinklO单元离散导线。通过在ANSYS中设置LinklO单元的特征参数,使其满足的导线物理特性。根据第二章所介绍的导线悬链线找型方法先确定导线初始形状一,在建模时每隔10m设个节点,跨长不足10m或遇特殊节点时取相应特一殊长度,共建立35个节点及左右两端各个固定点,用LinklO单元将其连接。导线结构阻尼比取0.005。4.3.2导线模型施加平均风荷载导线处于B类地貌下,取10m高度处的风速为27m/s。按施加重为加速度后所得的模-型高度,W么式(415)计算得各离散自由度处的平均风速。。=—氏x(4-15)K()’〇10上式中,Z为各点高度,Fio为10m高度处风速,a是澈流度,根据我国荷载规范,针对B类地貌,a取化15。巧运用准定常假设,按下式计算得到各节点上平均风荷载2=C-F心扔况(416)d去上式中=x:Q是导线的体型系数,P为空气密度面积。Ad^d是导线,A是受风的的直径,L为节点两端各1/2单元长度之和。0为风荷载方向与导线夹角,当风荷载W垂°直导线方向施加时餐90。在各节点位置上施加所得到的平均风荷载一,再对此模型在ANSYS软件中进行次找型,代替原模型(下文中未施加平均风荷载的模型称为原模型,施加平均风荷载后的模型称为新模型)。新模型就是导线施加了平均风荷载后的情况,若在新模型上施加脉动风荷载则可W认为只发生小位移变化。表4-2给出了施加平均风荷载前后模型各节点坐振值风巧载W垂直于导线的方,平均向(沿Y轴方向)施加,因此沿导线方向(沿X轴方向)不会改变。50 浙江大学硕击学位论文巧4章利用频域法计算导线风偏表4-2施加平均风荷载前后模型的坐标离散自施加平均风荷载前施加平均风荷载后由度编号X轴Y轴Z轴Y轴Z轴端点10038.000.0038.00..1号10036790.92371420.03583号2.661.7氏34330034.巧2.巧35.59号4号40033.593.3434.895号50032.664.0434.246抓031.794.68巧.64号70030.7.99527巧.08号8号80030.265.抑拍.589号90029.596.2732.口.10号100028.9967031.72...11号11002845707313712.巧.号120027987.310713号130027.587.6630.8114号140027.247.8730.6115150026.号.9780330.45160026..16号.778143034...17号17002669820302918号176.50化.村8.2130.2819号180026.678.2130.2820190026.748..173032号21号200026.878.0730.4111亏"0027.087.9230.5523号220027.367.7230.7424230027..7174730.98号25号240028.147.1731.2726号250028.636.8131.61260029.20627号.4032.00.855.巧号270029.94324429号280030.565.4332.93302撕031.354.8633.47号300032..31号.21423340532号310033.153J634.69巧号320034.162.拍35.3834号巧0035.242.似36.1135号340036.:391.1936.90端点23巧038.000.0038.0051 浙江大学硕击学位论文第4章利用巧域法计算导线风偏-从表42可1^Y轴方乂不难看出,原模型在施加了平均风荷载后除端点外,其余各节点向与Z轴方向的坐标都产生了改变总体而言:由于在跨中的节点约束较端点要小Y轴,,与Z轴坐标改变值较端点处的改变值要大。当导线受到垂直于导线方向的平均风荷载作用后导线节点的竖向高度会发生改变竖向高度改变值相对于节点高度是一个小值因此,,,而产生的平均风速改变可W忽略不计。上表中18号、19号节点竖向位移改变最大,为3.61m这两个点在位置改变前后的平均风速分别为:31/s.28m/s与31.88m改变值很小。,,导线结构的自振模态分布密集,高阶模态对于风致响应也有贡献,因此在运用频域法计算导线风偏时考虑多阶棋态的参与,W及不同模态之间的精合效应。作为尝试,选取施加节点f均风荷载后的新模型前20阶模态作为参与模态进行频域法计算。前20阶模态的-频率如表43所示:表4-3施加平均风荷载前后模型的自振频李对比棋态阶数施加平均风后导线樸未施加平均风导线棋型的频率/Hz型的频率/Hz1阶0.1973801.68552阶0.393260.巧5923阶0.巧490.巧7254阶0.4M250.440815阶0.巧3190.5065960.655710912阶.57079196067641阶..8阶0.792810.677089阶0.993树0.8485919085722…阶.000.11阶1.19651.021912阶1.19711.022313阶1.401711971.14阶1.40:371.199315阶1.61051.375516阶1.6111.巧5814917阶.82071.55118阶1.8214.5557192037517401阶..20拚2.03791.7404由上表可^乂看出施加平均风后导线模型各晚频率值均要大于未施加平均风前的频率,这是""由于施加平均风后,使得导线结构处于納紧状态,提高了导线的刚度,使得导线的-5白振周期变短;,从而提高了频李值。前6阶模态形式如图4所示52 浙江大学硕壬学位论文第4章利用频域法计算导线风偏呵'n;^^mw.>巧1:^誌:;腾品?^w.ecs?.mom*:e?OMts第1阶模态第2阶模态":;:;;:;^^^new仇hi?>TX???}*.4巧*CK"S3"ms?9C*、?u.?r.<7.ntcezwou?i第3阶模态第4阶模态mNsitwi?*?*巧"*>.<SJ7uno.mur£.‘?K.?>??wK.mm第5阶模态第6阶模态图4-5采用频域法计算所取的前6阶模态一-图45给出的是此跨导线在施加平均风荷载后前6阶模态的振型图,从图中可W看一出导线的的模态振型形式:第种如第1阶、第3阶、第5阶等等WY轴,主要分为两种位移为主的模态;第二种如第2阶、第4阶、第6阶等等L)(Z轴位移为主的模态。在现实情况下,导线受到的风荷载W水平向为主(即平行于Y轴方向),那么导线的风偏响应将主要受到水平振型方向为主的模态所控制。53 浙江大学硕击学位论文第4章利用频域法计算导线凤偏4.3.3颊域法计算结果将上述施加平均风荷载后的有限元模型参数,取前20阶模态作为参与模态,运用模一二态叠加法atbb编程软件中进行编程计算。值得提的是,在第章中己证明:在,在MASCENo74时域计算过程中中采用美国规范.所提供的气动阻尼值与本文所提出的考虑气动阻尼方法所得的计算结果十分接近。因此在频域法计算过程中,取结构阻尼与ASCE。N0.74么式计算得到的气动阻尼的和作为结构的总阻尼,来考虑气动阻尼的影响再按第二章所述风荷载模拟方法,获得此单跨导线的时变风荷载,在ANSYS中W考虑气动阻尼的计算方法求得各节点的时程响应,并对时域计算结果与频域计算结果进行对比。下图所列出的就是时域法算得的位移响应标准差与频域法算得的位移响应均方根对比。-?-SR巧法盈吝nlu**?CQC法巧合°-‘*^1时域?1^到的掠准差ioj/\i:::/\■.IIIIII.It.I.I...Q■Q0510152025303540节点编号图4-6时域法与频域法计算得到的位移均方根对比4-时域法计算得到巧标准差是由脉动风荷载所引起17)计算得到,含义,可由下式(与由颇域法计算得封的均方根相同,两者具有可比性。2-平均响应时域法响应标准差总响应)WW-(417):-417)中N代表各节点处所受的总荷载步数,式(。为了防止导线风偏闪络现象的发生一L-1在架空输电线路设计时般乂导线风偏的最大值作为控制条件,式(48)可W计算出响54 浙江大学巧击学位论文第4章利用频域法计算导线风偏应极大值;4-1(8)-式(418)中,r是响应平均值,在计算频域结果时取施加平均风荷载后的新导线模型Y轴坐标值。CTf是脉动响应的均方值,g为保证因子,当安全保证率为98.61%时,g。取2.2-■一SRSS症姐爸^* ̄CQC法姐合r-^1时祐洁结染%!\I■II.I.t.I.I.I.I.0己10152025303540节点编号图4-7时域法与頻域法计算得到的位移极值对比SRSS巧巧合SRSS法由备—*-eCQC法巧合5.Or-?-IIl法组合hr<ii\,U..5-I。./X苗./VW■.I.i‘-■!■?I.I.I()〇..?I代?I-|??玄m.!I.I■I051015如巧3035400510巧20巧泌巧40节点编号节点巧号(a)绝对差值化)相对差值图4-8时域法与频域法位移响应的绝对差值与相对差值图4-7给出了采用频域法与时域法所求得的位移响应极值对比。为确定频域法与时域-CQC法与SRSS法两种不同组合方法计算法所得响应的具体差值,图48给出各节点上按得到的响应极值与时域法计算值的绝对差值和相对差值。从图可见当选用CQC法组合计算的响应值较SRSS法更加精确,最大差值绝对值约为8mm。端点处由于约宋较好,总响应值小,相对差值最大,约为2%。选用SRSS法组合时精度也令人满意,最大位移差值约22mm,最大差值比率5%。这说明按照本文所提出的频域法计算响应值可W满足计算55 浙江大学硕击学位论文第4章利用频域法计算导线风偏一的要求,仅受脉动风荷载影响的导线结构可视作线性结构这假设成立。采用20阶模态进行组合时,,精确度较高,但计算过程也较为复杂。为确定在保化精确度的前提下选取1、2几阶模态进行计算有较高的计算效率,分别再选取0阶、5阶、3阱阶模态作为参与模态进行组合。前文中己介绍,此单跨模型前2阶模态分别是WY轴、Z轴位移为主的模态,导线主要考虑这两个方向的运动,因此至少取2阶进行组合。将获得的位移均方根与时域法进行对比,并给出时域、频域计算响应结果的相对差值,如下图所示:^^IV臣ISnIIi議^-.**。.-、*1、0*/A苗.5?????》?%/_"0I.I.I.I.II.I.1■t?0Q■Q05101日20巧30化400日1015202日30巧10节点巧号节点编号(a)均方根对比化)位移相对差值4-图9取10阶模态进行姐合--!?SSSiSPSS法殖吝__RiCi5‘〇-夺-化7"CQC法巧合rWh〔QC法组合|时试法巧到巧标准差■4.5■\r\*05101520253035400510152【)2日303540节点巧号节点巧号(a)均方根对比(b)位移相对差值图4-组合10取5阶模态进行56 浙江大学硕±学位论文第4章利用频域法计算导线风偏—SBS法坦合—SR班法品n7rII5Or,II-亡贼数■!.4;s,\义繼f:\p..".15|.2/\害芒。'I\;r一:;/i..I,I.I.?.I.1.I.I.IQKI■?-I.I.1.I.IIIQ'Q5阳520巧孤况仰051015如25303540W1节点编号节点编号(a)均方根对比(b)位移相对差值4-图11取3阶模态进行组合〇.’■[lia.巨::\E]户.lAiVy它。..1\::i::'—一 ̄ ̄ ̄ ̄ ̄ ̄ ̄—■ ̄__—?一^_—___________OL<■■,■,■■0.0t_i.I.I.I■一i.I.I.I.I.OtII广I广III05况说4005101520253035405阳巧洗2节点巧号节点巧号(a)均方根对比(b)位移相对差值图4-12取2阶模态进行组合-9至图4-从图412可W看出:频域法进行计算时,当选取的参与模态减少,计算得到的位移响应均方根与时域计算结果差别增大,采用CQC法与SRSS法计算结果差别变小,当选取2阶模态作为参与模态时,CQC法与SRSS法结果几乎完全重合。这是由于选取的模态越少,CQC法所能考虑的模态间稍合作用越小,故与SRSS法计算结果接近。取20阶模态进行組合的SRSS法计算值与取2阶模态SRSS法结果十分接近,说明导线风偏过程主要由前两阶模态所控制。当取前2阶模态进行组合时,频域与时域计算出的位移响应相对差值均小于5%。说明只考虑2阶模态的姐合即可获得较高的精确度,这样的计算方法能大大减少计算量。4.4本章小结V本章分析了非线性结构在使用频域法计算风致响应时的条件及假设,提出了导线采用L一频域法计算风偏响应的方法。乂跨导线结构作为算例,按本章方法计算风偏响应,并与57 浙江大学硕击学位论文第4章利用频域法计算导线风偏时域法结果进行对比。得到了W下主要结论:(1)平均风荷载所引起的响应占总响应的绝大部分,受脉动风荷载影响时导线只产生小位移。因此可将只受到脉动风荷载作用的导线模型视作线性结构。在考虑了平均风荷载作用的导线模型上采用频域法计算风偏的脉动响应,再与平均响应相加即可得到总响应值。(2)本章所提出的导线频域计算方法所计算得到的风偏响应结果与时域法结果相近,说明本章所述方法可L乂在实际工程中采用。用频域法来计算导线的风偏能够很大程度的提高计算效率。(3)导线各阶模态间存在着相关性,因此采用考虑模态间相关性的CQC法计算精度较不考虑相关性的SRSS法精度高,1:^选用SRSS法组合时精度也令人满意。CQC法计算需考虑模态相关项,带来许多计算量,因此若在计算时需要较高的计算效率,可采用SRSS法进行组合‘。(4)取前20阶模态进行组合的SRSS法计算值与取前2阶模态的SRSS结果十分接近,这说明导线风偏过程主要由前两阶模态所控制。(5)取前2阶模态进行组合时5%。,频域与时域计算出的导线位移响应相对差值均小于在实际工程计算过程中,为提高计算效率,可只考虑前2阶模态的影响,采用SRSS法进行組合。这样计算的结果也拥有较高的精确度。58 浙江大学硕壬学位论文第5章导线的等效静力风荷载第5章导线的等效静为风婿载5.1引言一自然风具有脉动特性,对结构的作用般需要进行复杂的动为学分析。在实际工程中,为方便计算,设计人员常将风荷载与其他荷载进行组合后采用静为分析。为使随机动力风荷载数据能够直接应用于实际工程中,学者提出等效静力风荷载的概念:将等效静力风荷载作用于结构上时结构所产生的如位移一、应力、应变等效应与外加动为荷载时定保证,一率下的风致峰值响应基本致。tW阵风荷载因子(GLF)法是由Davenport所提出的最早求解等效静为的方法,其基本思路是通过阵风因子来表这结构对脉动风的动为放大作用。但当平均风荷载或平均响应72[Kl接近于零时,阵风荷载因子法就不能够计算出合理的等效静为风荷载值。asperski通过--对于低矮房屋的等效静力风荷载研究,,提出位移响应相关(LRC)法其实质是求解了背景响应的等效静力风巧载。LRC法提出后,学者们意识到应对风致响应中背景分量与共nzhon振分量的等效静为风荷载分别求解,进行组合后得到总响应的等效静力风荷载。兒gwKMtChlenAhsanKareem等学者提出将等效静为风荷载表这为平均分量、LRC法下的背,景分量和惯化力形式表达的共振分量的线性组合,并通过权重因子使等效目标响应与实际风荷载作用下响应极值完全相同。目前学者们研究关注于如何高效率、高精度地求解高层建筑大跨结构等建筑物的等效静力风荷载。作用于导线上的风荷载将传递至输电杆塔一个非常重要的水平荷,是输电塔所受到的载导线的等效静力风荷载研究,还鲜有涉及。这主要是由于导线是明显的非线性结构,,无法在基础模型上直接求解等效静为风荷载。在第四章中曾介绍:考虑导线受平均风荷载作用达到平衡位置后,再施加脉动风荷载的过程可将结构视作线性。因此可W计算出平衡位置后导线模型的等效静力风荷载再与平均风荷载组合,得到总响应的等效风荷载。本章首先介绍了结构的风致响应类别^及对应的等效静为风荷载计算方法,再对导线风偏的过程在频域上进行分析,发现风偏响应主要由背景分量所控制。然后采用LRC计算处于平衡位置的单跨导线模型等效静为风荷载。将等效目标响应与施加气动力模型的响应进行对比,两者十分接近。这表明采用LRC法来计算导线等效静为风荷载切实可行。59 浙江大学硕壬学位论文第5章导线的等效静力凤荷载5.2风致响应分类导线所受的自然风中包含着平均风与脉动风,任意t时刻结构Z高度处所受到的风速可化写为;vzt=vz+vzt-(51),,()()()其中云Z代表平均风速,其数值不随着时间的变化而改变;代表脉动风速,不();罔时刻的脉动风速值不相同:。根据准定常假定,所引起的风荷载值可表示为2=Cyzfdp(,)去=zClvzvzt-+t+vz(52)^p\{,,{)){){)^=tFz+Fz{{,))式中222,歹2是导线在高度处由平均风引起的平均风压;歹是导线在高度处的,()(〇一般情况下,脉动风速小于平均风速云Z脉动风压,p为空气密度,Cd为体型系数。,()故脉动风速的平方项将远远小于平均风速云Z的平方项W及它们的交叉项。因此在()2计算脉动风压时,通常忽略项。7;2<。导线2高度处节点的脉动风压可表示为:,/(,)^=-FztCvzvzt(53),^p,{){)[)当平均风压Fz与脉动风压作巧于导线上时,将分别产生平均响应与脉动响(;)应。脉动风压对于导线结构而言是动力荷载,导线受到其作用时将在自振频率处发生共振,这部分响应称作共振响应,共振响应与结构的模态、阻尼、质量、刚度等动力特性密切相关,这;在非共振频率处,脉动风压也会产生响应部分响应成为背景响应,背景响应与结构动为特性无关主要受到风荷载频率分布特征的影响一,是种准静态的响应。因此风致,响应可W表示成如下形式:'平均响应f共振响应风致响应脉动响应<背景响应60 浙江大学硕壬学位论文第5章导线的等效静力风荷载5-:若将风致响应体现在功率谱上,1所,W谱密度的形式表达如图示5-图1背景响应及共振响应5-图1中JT代表背景响应;乂代表共振响应,其尖端位置所对应的频率是结构的自。g振频率。5.3导线风偏响应在频域上的表现为确定导线风偏响应在频域上的分布2-24位置处绝缘子,将图所示的八跨导线法点5-2(a)串风偏角时程结果进行傅里叶变换转换为功率谱。在图、(b)中分区段给出。—不考虑气动阻尼…考i与动阻尼4035--30\2-5S20-I-1m15叫-hlil/LyulLIA.IfjiiAJlAM一一>一j0.020.040.060.080.100.120.14频率(Hz)(a)背景响应61 浙江大学硕±学位论文第5章导线的等欢静力风荷载—不考虑……考揭动r-1401-120I^1W1ICIIAiJLJikiw^Ii....L,-.,.,\Q0150)0.250.300.040i.公巧.频率(化)(b)共振响应■图5-2挂点4的风偏角响应在频域上表现-2化-图5)中的0.160.1細Z区段内功率谱有明显的峰值,结合导线的自振频率可知,5-2(a)中也存在着响应这些尖峰代表动态响应的共振分量。在图,这部分主要为背景分量。5-2从图可W看出:无论是否考虑气动阻尼,在风偏的背景分量频域范围内,两者的功率谱曲线均集中在风致振动响应的低频部分且几乎重合;在共振分量频域范围内,不考虑气动阻尼所得到的风偏响应共振分量有若干个十分明显的峰值,引入气动阻尼后,共振分量则被大幅削减。现实情况下,导线受风荷载作用产生风偏响应的过程中总是受到气动阻尼影响,因此在实际工程应用中可W忽略共振分量对风偏的贡献。另外需要说明的是,输电线路属于非线性特征显著的高柔性结构,在风荷载作用下产生的应为刚化效应会使得结构共振响应频率略高于静为平衡状态下的自振频率。结合图-52b-二()不难看出,发生风偏时共振分量对应的频率范围为化160.18HZ,而不在第章 ̄中所给出的导线初始自振频率化130.16Hz之间。5.4导线风致响应的计算方法及其对应的等效静为风荷载5.3节中己介绍:实际情况下,输电线路的风偏响应主要为背景响应,共振响应很小可W忽略。因此,在计算时主要考虑导线的平均响应与背景响应及其等效静为风荷载。5.4.1平均响应首先讨论导线平均响应的等效静力风荷载。根据静力学导线在平均风荷载只X的,();62 浙江大学硕击学位论文第5章导线的等效轮力风巧载x?作用下〇位置处的/响店为:,=^护义/义化5-4而,)()()()心〇_[式中:r代表了导线响应的类型,输电线路在建设过程中主要考虑避免导线间距过近发生风偏闪络一,因此般讨论的导线响应是其位移;足杉〇,,句为初始导线模型的影响线系数其物理含义是在任意X位置施加与歹(句同向的单位力时,在;C0位置产生的r响应值的大小。平均响应和平均风荷载都是静为值,平均风荷载就是平均响应的等效风荷载。5.4.2背景响应第四章中邑证巧:考虑了平均风荷载影响的导线模型在脉动风荷载影响下,可视作线8283tKl性结构。导线的背景响应可由方差分析(covarianceanalysis法来进行求解。根据结构)力学知识,当导线在脉动风巧载片作用下;C位置处的脉动响应值可表边为:,=-f^义,乂55而批()()(从(而]此处J(X;C是施加平均风荷载后的导浅模型的影响线系数,那么脉动响应的平方d0,)为:2f义=护义义如-56A2,,(。)(l而l而2如())(K(K()脉动风具有随机性:,因此考虑其统计值才有意义。响应的均方差可W表示为2f=护X护义-7兩,|兩l而l2叫口)()(K,AK而J)((()因此响应的均方根为:"=_C欠义义JX义■*义瓜把5B而,f,rl,巧.()J卢(12)(。I)(02)12({式中:为导线;C、X,两点上风荷载的协方差。从上式可W看出,背景响应i与结构的动为特性无关。Kasperski为了求解刚性低层房屋的等效静力风荷载,提出了LRC(Load-Resonse-Correlation)法LRCp,低矮刚性房屋的风致响应主要由背景响应所控制,法就是求解背景响应的等效静为风荷载。受平均风荷载影响下的导线结构可视为线性,因此其背景响应的等效静力风荷载也可W采用这种方法进行求解,主要计算过程如下。5-8)背景响应的均方根式(可写为:63 浙江大学硕击学位论文第5章导线的等效静力风荷载=X义如_9U而,而,,口),21兩12)片()[侣(化()於5-9式()中括号内的积:分值就是响应与荷载的协方差,即。=KXXZXdx(,,:ffi2)l,(〇,),,{5-10()=义义。而。义^(。,),。2.(.片)卢()式中:是背景响应与荷载的互相关系数W写为:,可X/jfX(kI馬(,,兩)"(〇,,)i=- ̄ ̄—-P口11F)r’"〇Xr.aU)>(2)5-105-将式()代入到(9)后便有:^=xX^^^5-^12rBh,)\{^P22(.M〇)rFF(),j若定义:=5-1巧*(引,WW式中-2:说为背景响应的保证因子。那么(51)可w写为:^^=^gerB^.^d^例句,(0(02VVS22),()1从式5-H)不难看出FX(,rs是背景响应的等效风荷载值与发生.()就,其响应的位置及类型密切相关。5.4.3总响应忽略了共振响应导线的总响应由平均响应与背景响应组成:,可W表示为,=-rr+(X5"gB,e()上式可写为:r=x5-Fjc/A:xdx+Fx/jcdjc16()^(〇,)()^,(?,)()jj^在初始模型上分别施加F切和巧,的得到的总响应可W认为是准确值,由于导线是非线化结构一,在初始模型上直接施加歹(JC)与F(;c)的和所得的响应值与准确值存在定fj>差距,但影响很小,因此总响应的等效静力风荷载可写为:fz=护z5-z+17r()()巧g()().64 浙江大学硕击学位论文第5章导线的等效静力风荷载5.5导线等效静力风荷载求解W第四章中介绍的平均风荷载作用后的单跨导线新模型作为算例。取出其水平方向位19移影响线矩阵,并提取出所施加的脉动风荷载时程,W水平位移响应最大值的跨中号节点作为位移等效点,为方便比,采用LRC法计算导线风偏背景响应的等效静力风荷载较,在利用LRC法计算的过程中取背景分量保证因子g与计算风偏响应时程统计值所取一=225-17)的峰值因子相致,.。忽略导线的共振响应再通过式(进行姐合,即可得到g,此一跨模型的导线风偏总响应的等效静为风荷载。导线的等效静为风荷载与输电杆塔上其他荷载狙合后供设计抒塔时采用。规范中化给出了设计巧塔时导线风荷载的计算方法,,可可将两者进行对比。,为方便比较本章在计算导线等效静为风荷载时不考虑在跨度范围内平均风的折减,《=当采用规范方法计算此单跨导线上各节点的风荷载时,取不均匀系数1.〇。此外,还有两点值得注意:(1)不同于计算导线风偏时A取1.0,此时风荷载调整系数A应按照作用于杆塔上的导线风荷载取值,对应于lOm高度处基本风速27m/s时,A取1.2。C=(ll2)当采用规范么式计算导线风荷载时按照规范要求取导线的体型系数d.,而一根据些风洞试验结果.。,导线的实际体系约为10因此,在计算等效静力风荷载时体型=系数取Cdl.〇。-?一等效廚力风荷扱?,曲巧a化凤拉掛I2-.4f2-12-.2/?2〇??'??记??"??/??????巧勺8\.、8-拽¥.\/16\[6/::IV/八1..巧;'.V2:::/:Vy、--00.8I1I>II■I■I■I■I.I■I..■■I■I...If05101日202日3035400日10152025303540节点编号节点编号(a)沿导线的风荷载分布化)两者相对差值65 浙江大学巧去学位论文第5章导线的等效静力风荷载0.3〇1I<二I^齡I0-.260屠言Z0-.22化-专160-.14-0.12J1■??'■'0—?l〇0.20.40.60.81.0x/L(c)导线上的分布风载图5*4等效静力风荷载与规范算得风荷载的对比.-4给图5出了等效静力风荷载与采用规范算法得到的导线风荷载,从图(a)中可W看出,,两个风荷载在沿导线跨度方向具有相同的分布,跨中位置处风荷载有明显的降低这是由于在导线建模时,跨中附近位置有若干个节点的间距较小,相应的节点控制长度较短,,导致这些节点位置处的风荷载值较小将导线上节点风荷载转换为分布风载后(图""一(C)),,就可W发现,采用规范算法算得的导线分布风载是条光滑的曲线总的等效""静为风荷载中考虑了背景响应的等效静为风荷载,因此并不光滑;当不考虑不均勾系数时,规范算得的导线风荷载与等效静为风荷载十分接近,两者最大的相对差值约为10%。将计算得到的等效静力风荷载施加到单跨导线原始模型上,得到等效目标响应,并与实际动为风荷载作用下的风偏响应统计值进行对比:,结果如下所示等效静力风荷我作用■—I实际請方风荷载作用-10I.I.I■I.I.I.I.I.I.0510152025303540节点编号图5-5各节点水平位移响应66 浙江大学硕击学位论文第5章导线的等效静力风荷載65-.■60-./55-§.\f■、'-E四、■/-■45寡./妇\-■/4、.0■和■■■若■■■■■-3.51tIIIIIII....■■.300510152025303540节点编号图5-6水平位移相对差值5-5从图可看出:受等效风荷载作用下的水平位移响应与实际动为风荷载作用下的响应基本一致,证明了在施加平均风荷载后的导线模型上采用LRC计算等效静力风荷载的可行性19,但等19;号挂点虽是目标位移等效点效静力风荷载作用下的号挂点并未与实际情况相重合,这是由两部分原因所造成的:首先,导线风偏过程中共振响应虽然很小,但对于总响应仍有影响。LRC法只计算了背景响应部分的等效静为风荷载,只受此等效静力风荷载作用的风偏响应值将小于实际风偏响应,无法完全重合;其次,本文采用LRC法计算得到的背景响应等效静为风荷载是基于受平均风荷载作用下的导线模型,因此在初始模型上依次施加平均风荷载与背景响应等效静为风荷载可W获得不考虑共振分量的总响应准确值一。由于初始导线模型是非线性结构,在此模型上次性施加平均风荷载与背景响应等效静为风荷载的和所得到的响应值与准确值存在一定差距。从图5-6可レ乂看出:吕标位移重合点:19号挂点,其等效目标响应与施加气动为模型的响应相对差值最小,约为3.5%,最大值小于6.5;相对差值从19号挂点向两端不断增加%,这说明采用L民C法来计算导线等效静力风荷载精确度较高。5.6本章小结本章考察了导线风偏响应在频域上的表现,再W单跨导线模型为例,采用LRC法求解其等效静为风荷载得到L乂下主要结论:,(1)现实情况下,导线受风荷载作用产生风偏响应的过程中,主要由背景响应所控制,共振分量受到气动阻尼的影响而被削减,可W忽咯。导线受平均风荷载影响后所形成的新模型,其施加脉动风荷载产生的脉动响应的等效静力风荷载可由LRC法计算。67 浙江大学硕±学位论文第5章导线的等效静力风荷载(2)当不考虑不均勾系数时规范算得的导线风荷载与等效静为风荷载十分接近,,两者最大的相对差值约为10%。(3)由于LRC法只计算了背景响应部分的等效静力风荷载,只受此等效静为风荷载作用的风偏响应值将小于实际风偏响应;其次,本文采用LRC法计算得到的背景响应等效静为风荷载是基于受平均风荷载作用下的导线模型,因此在初始模型上依次施加平均风荷载与背景响应等效静为风荷载可W获得不考虑共振分量的总响应准确值。由于初始导线模型是非线性结构一,在此模型上次性施加平均风荷载与背景响应等效静为风荷载的和所得到的响应值与准确值存在一定差距。这两个原因导致目标位移等效点的响应无法与施加气动为模型的响应完全重合,但相对差值较小,可W满足实际工程计算需要。68 浙江大学硕壬学位论文第6章结论与展望第6章结论与展望6.1本文主要结论在考虑气动阻尼效应的输电线路风偏动态计算方面:,本文主要结论包括1)(输电线路在受风荷载作用下时,导线产生风偏响应。在产生响应的过程中导线与气流的相对运动会产生显著的气动正阻尼一,在高风速下该阻尼般大于导线本身的结构阻尼。,不可忽略。当考虑此因素时,能够显著降低导线风偏响应的脉动值.(2)根据有限元软件计算得到的绝缘子挂点处风偏角响应可知:靠近耐张塔挂点处的绝缘子串风偏角会比耐张段中部挂点处风偏角要小,这是由于支座约宋效应的影响,对""作用一耐张段两端附近绝缘子存在着拉扯。实际线路中也普遍存在这现象,在日常设计一过程中若考虑这因素的影响。,可W降低输电线路的建造成本(3)美国输电线路规范《ASCEN0.74》中提供了导线气动阻尼值的计算么式,通过此么式算得八跨导线的气动阻尼值一部分阻尼,在有限元程序中通过瑞利阻尼么式考虑这值,并计算出风偏响应时程值、统计值与本章所提出的考虑气动阻尼的方法有限元结果进行对比,发现两种算法算得的风偏响应十分相近。证明了ASCENo.74所提供的简化计算方法结果较为准确,可供日后风偏研究时采用。关于我国规范风偏计算方法及与多国规范对比:,本文主要结论包括(1)在计算导线风偏过程中,规范推荐采用的单摆模型具备原理明确、计算简单的特点,而实际情况下导线多为连续多跨,各导线之间存在着賴合作用。通过本文计算发现,"拉扯"效应即导线的风偏响应受到相邻导线的抑制这种稿合作用体现为,,多自由度连续模型算得的风偏角比单摆模型算得的风偏角要小。从而证明在实际生产设计中,采用单摆。模型计算导线风偏响应偏于安全,可在工程计算中采用(2)有限元软件算得的风偏响应极大值与平均值的比值,即是考虑脉动对于风偏响应的影响系数。通过文章中的对比可W发现,导线受脉动风荷载作用时,所得风偏角放大约 ̄ ̄1.151.3倍,对应的的水平风荷载放大1.31.5倍。我国在风偏计算时取风荷载调整系数为1.0偏于不安全。(3)各国规范中关于导线风偏的风荷载均基于准定常假设,因此可将各国风前载规范表达为同一个算式。对比各国规范关于算式中参数取值的不同,可W发现,我国所取的69 浙江大学硕击学位论文第6章结论与展望风速、体系系数与其他国家相差不大,但关于风压不均匀性及动力放大两个因素的综合考虑与他国相比较小,使得我国风荷载计算值偏小,这是导致我国风偏闪络现象频发的主要原因。关于利用频域法计算导线风偏:,本文主要结论包括(1)导线风偏过程中平均风荷载作用而产生的平均响应占总响应的绝大部分,,受受脉动风荷载影响时导线做小位移运动。因此可将只受到脉动风荷载作用的导线模型视作。线性结构,在频域范围内计算脉动响应值,再与平均响应相加即可得到导线的总响应值(2)对比频域法与时域法计算结果,可W发现两者所得的风偏响应值接近,说明可W采用频域方法*计算导线的风致响应值,这样的计算方法能够提高计算效率。(3)导线各阶模态间存在着相关性,因此采用考虑模态间相关性的CQC法计算精度较不考虑相关性的SRSS法精度高,采,但两者相差不大用SRSS法组合时精度也令人满意。(4)取前20阶模态进行组合的SRSS法计算值与取前2阶模态的SRSS结果十分接近,这说明导线风偏过程主要由前两阶樸态所控制。当取前2阶樸态进行组合时,频域与时域计算出的导线位移响应相对差值均小于5%。因此在实际工程中,可只考虑前2阶模态对于导线风偏响应的影响,采用SRSS法进行组合。这样的计算方法能够在保证精度的前提下拥有极高的计算效率。关于导续的等效静为风荷载:,本文主耍结论包括(1)通过研:究导线风偏响应在频域范围内的表现,可W发现现实情况下导线受,风荷载作用产生风偏响应的过程中,主耍由背景响应所控制,共振分量受到气动阻尼的影响而被削减,可W忽略。因此,可采用LRC法计算导线受平均风荷载影响后所形成的新模型的等效静力风荷载。(2)由于LRC法只计算了背景响应部分的等效静力风荷载,只受此等效静为风荷载作用的风偏响应值将小于实际风偏响应;其次,本文采用LRC法计算得到的背景响应等效静力风荷载是基于受平均风荷载作用下的导线模型,因此在初始模型上依次施加平均风荷载与背景响应等效静力风荷载可W获得不考虑共振分量的总响应准确值。由于初始导线一模型是非线性结构,在此模型上次性施加平均风荷载与背景响应等效静为风荷载的和所得到的响应值与准确值存在一定差距。这两个原因导致目标位移等效点的响应无法与施加气动力模型的响应完全重合。,但相对差值较小,可W满足实际工程计算需要70 浙江大学硕+学位论文第6章结论与展望6一.2进步工作展望(1)实际情况下,导线受风荷载产生风偏过程中,导线的运动使得风荷载各时刻的。,风攻角发生改变,从而体型系数会产生变化本文在计算过程中体型系数考虑为常数因一一定出入而与实际情况存在。同时在现实输电线路中导线处于个王维风场环境中,本文只考虑了其主耍风向(即垂直于导线的方向)的影响,更精细化的有限元计算值得深入分析。一(2)输电线路规范中些风巧载计算参数包含了对于多个影响因素的考虑不同,而国家对于相同影响因素考虑可能存在于不同的参数中,且并未给出考虑单独影响因素的计算方法,这给对比工作带来了困难。若可W将规范中各影响因素分离成单独的参数进行对比,则会使得结果更为直观,有利于我围规范的改进。(3)本文在频域法计算导线风偏时采用了单跨模型,这样的研究过程具有基础性。一般而言,实际输电线路为^跨线路,后续可继续研究频域法计算多跨线路的精确度W及计算效率问题。71 浙江大学硕壬学位论文参考文献参考文献-I.张龙伟.特高压共走廊输电线路电滋稱合分析及改善措施研究[D].成都:西南交通大学2014:12,[]-.我国500kV输电线路风偏闪络分析J.2005。9(76573[电网技术).巧张禹芳]-500kJ.200923(59902.口]肖东坡.V输电线路风偏故障分析及对策[电网技术:)1]uhsofa-uzak巧aleldlowwureconduc4Hiratsuka8MatsiYFukudaN?巧testresindresstorC.[],,,p[]Proceedna-ingsofIEEERegion10InternatiolConferenceon.IEEE,2001,2:664668.l5StClairJGClearancecalculationsofconductorsbuidi打sC.TransmissionandDistribution]g[][Conference996*IEEE1:4巧498.,,6AllenL.Calculationofhorizontaldislacementofconductorsunderwindloadintowardbuildingsand[]pgrtittuCrsPresented-thsrucres.Paat化e37化AlConference.IEEE1993Al/lAinO.o舟suppong[]pennua,:[7]TsujimotoK,YoshiokaO,OkumuraT,etal.InvestigationofconductorswingingbywindandUsapplicationfordesinofcompacttransmissionlineJ.PowerAaratusandSstemsIEEETransactionsg[]ppy,on-198211:43614369.,()8RikhVN.Conductorscingsintransmissionlinesandeffectoflonsanswithsteeslopesinhill[]pagppyterranJourna-iJ.IEIl2004855:816.[](),,()9ChayMT,AlbermaniFawesH.Windloadsontransmissionlinestructuresinsimulated[]QHdownbursts[C],FirstWorldCongressonAssetManagement.2006.立宏胡毅李景禄等2006329-。巧龙.输电线路风偏放电的影响因素巧究阴.高电压技术4:121.,,,,,()ILl.PsrwindcareadnwnuTENCONILihononetaarameterfousedovehtransmissiolinesifalt.[]g,g,gand2006.2006IEEERegion10Conference.IEEE,2006.n巧李孟春,张艳玲,杨北革等悬垂绝缘子串风偏最小间距距离计算分析[J]电测与仪表20-,口,49(555):71013王声学.500kVJ.2008[],吴广宁范建斌,等翰电线路悬垂绝缘子串风偏闪洛的研究[电网技术,,,]329-669():5.-14李瑞祥,张志远导线悬垂绝缘子串风偏角的解析计算饥.电为建设,199415(4);1921[],-15徐钟济.蒙特卡罗方法[M],上海:上海科学技术出版社1985,62357[],。句ShinozukaM.Simulationofmultivariateandmultidimensionalrandomprocess故.TheJournalofTheAcousticalSocietyofAmerica.197172 浙江大学硕击学位论文参考文献。7]口eodatisQSimulationGfergodicmultivariatesU)chasticprocess的,JournalofEngineeringMechanics,SCE28-A19%12778787.,,(),-18张相庭工程结构风荷载理论和抗风计算手册M.上海:同济大学出版社,1990,10:3637[][]王之宏11-19风荷载的模拟巧究J,99454452[],[]建筑结构学报,1,(),20lannuzziASinelliPArtificialwi打denerationandstructuralrespon化JournalofStructural[],p,g,Enneerin-gig,ASCE,1987,113(12),23822398.21LiYKareemAARMAssteminwindenineerinrobabilisticEngineeringMechanics199052[],,ygg.P,,(),50-59.-张相庭结构风压和风振计算:.M?上海:同巧大学出版社98556090[巧[],1,口3ShinozukaM,JanCM,Digitalsimulationofrandomproc的sesandi化alications,JournalofSound]pp-andVibration197225111128.,,y),PWLiYousun,KareemAhsan.Simulationofmultivariaterandomprocesses:hybridDFTand出italfilteringg-aroach.JJournalofnreccEgineeingMhanis19%119510781098.pp[],()25袁波,应惠清徐佳讳.基于线性滤波法的脉动风速模拟及其MATLAB程序的实现阴.结构工程[],师2007-:5561.,04,口句謀佳艳,严波,刘小会,等,悬垂绝缘子串动态风偏响应的数值模拟阴.重庆大学学报(自然科学20062—版):100103,,%口)27林雪松严波刘仲全.220kV高压输电线路风偏有限元模拟研究[J.应用力学学报2009261[],,],,(),0-24121.2牵风偏角计算公式中风荷载调整系数研究D重庆:重庆大学2008.[巧林雪松悬垂绝缘子[],29-孙保强.201361128082813.[],侯鍾,孟晓波等不同风速下导线风偏动力响应分析P]高电压蚊术,化(),瑰讳,巧约董彦武等.输电导线悬挂点处等效风荷载仿真分析[J.高电201238002口巧邵,,]压技术,(),,476-482.[3U郑佳艳.动态风作用下悬垂绝缘子牵风偏计算研究.[D].重庆:重庆大学,2006.松等-3.500kV超高压输电线路风偏数值模拟研究.J工程力学2009244249.,严奴林雪,,[q刘小会[]y)33罗先国.特高压绝缘子串风偏研究.D,武汉.华中科技大学,2012[][]34李黎泻林海罗先国等.特高压绝缘子串的风偏计算方法饥.高电压技术2013392],,,,,(1,[)2924-2932.-35孔德怡..200829959[],李黎,龙晓鸿,等悬垂绝缘子串动态风偏角有限元分析机电为建设,();36楼文娟.[],杨悦,吕中宾,张少锋,杨伦考虑气动阻尼效应的输电线路风偏动态分析方法饥.振动与73 浙江大学硕壬学位论文参考文献冲击(录用)7巧询邵瑰讳文志巧等.国内外悬垂绝缘子串风偏设计参数对比与分析口],,[J],电力建设20-133441926,();3..昆明.昆明理工大学2013,[巧陈博氣中越风荷载与地震荷载对比研究[D]39红河J2003-袁勇.中越风荷载规范对比研究[.结构工程师1604:1720.[]阮,],,()40J.夏瑞光.[],范存新中澳欧风荷载规范的对比研究[]苏州科技学院学报(工程技术20-版),12,49(03):1318.-"屈论昭.国内外输电塔风巧载技术标准比较分析[J,电力建设2〇n05:22:29.,,[]]42张盈哲廖宗高谢麗-拍.输电线路设计规范中风荷载计算方法的比较J,电力设201307:57.,,,,[][]^[43]徐彬,冯衡,曾德森.国内外输电线路设计规范风荷载比较阴.华中电力,2012,86(02):7氏81.445034--薛芙.采用欧洲标准EN11计算输电线路风巧载的基本算法机.新疆电为技术,2010,04:9597.[]45]安旭文,王吴深,侯建国.国内外输电线路设计规范风荷载取值标准的比较机.武汉大学学报(工学[-版2011764〇6:740743),,()4汇国巧外输电铁堪设计标准对比分析化特种结构20920-,朱,化:102107.[巧王永华(句47潘峰,陈稼苗,轰建波,郭勇应建国.国内外规范输电线路铁堪风荷载特性对比阴.中国电[],力2饥"304:3742.,,-4.J.[巧袁玲,李庆祥,许伟,张春梅中日巧风荷载规范中风荷载与风振响应的比较[]广州大学学报(自然-%科学版)方):5014,01.4,邸玉贤,李新民,等.安装巧间间隔棒的输电线防.高电压技术[巧朱宽军风偏设计有限元分析m,364-:10381043.()口0}Yan8,LinX,LuoW,巧al.NumericalStudyonDynamicSwingofSuspension虹sulatorString虹OverheadTransmissionLineunderWindLoad口.PowerDelivery,IEEETransactionson,2010,2Sy):]248-2化-—5lm货J..D99.AliltticetII.Aridi[。Ho,"巧ongwndresponseofaowerseodynamcampnganddeflecnnrucures ̄tions口.EgineerigStt,187,483488.]()[52]ASCE:Guidelin的forElectricalTransmissionLineStructuralLoa出ng.Publishedby化eAmericanSocietyofCivilEngineersNo.74.2003.53LoredoSoiiza.ThebehavioroftransmissionlinesunderhihwindD.Ontario:Theuniversitof[]g[]ywesternOntario,1996.54孙珍茂.输电线路舞动分析及防舞技术研究[D].梳州:浙江大学,201化[]74 浙江大学硕壬学位论文参考文献uX--55YANBoLiiaohui,LiuBaoan.Dicresonseanalsisofwindaeawofoverhead[],ynampygynnesnonreachan9transmissioli[J.WSEASTransactiosAppliedandTheoticlMeics,2006,1(1):0.]56GB50009-2012建筑结构荷载规范[S.[]]?57GB-505452010110750kV架空输电线路设计规范[S].[]—5J.2742832.[巧滲宗高?特高压输电线路电力建设2006:,张华,陈海波设计风速取值的探讨[],,()59]EC608262003DesinCriteriaofOverheadTransmissionLines[S].[g-werand—60BS81001:19%Latticetosmastsarti:codeofpracticeforloading[S][]p6eeBE.Theefectofturbulenceon化esurfaceressurefieldofasrerism?J.Fluid.Mech.,1975,[。Lpquap口]'69-:2的282.62]KareemA.Windexcitedmo巧onofbuildingsD.PhD(dissertation,Dept,ofCivilEngineering,Colorado[[]StateUniv.FortCollinsColo1978.,,,63.随桃振动的虚拟激励法[M].北京2004林家浩,张亚辉:科学出版社,[]-64.27223[}林家浩,钟万热关于虛撒激励法与结构隨机响应的注记机计算为学学报9%15苗1,1,()65LinJ.H.AfastCalorithmofPSDmatricesforrandomseismicressesJ.ComPutersand[]QCgpon[]S—turetuers199244:683687.,,9)W—non66LinJ.H.n.S.Willima.Feudoittilrihmfoirrandomseisi[],Zhag,sW.Psexcaonagotrstatonaymcresn-巧6ponses.EgineeringStructures199416:270机,67.隙贤川.大務度屋盖结构风致响应和等效风荷载的理论研究及应用[D]浙江大学博壬学位论文,[]2005.6.大跨度屋盖结构的抗风巧究D.浙江大学博壬学位论文2004.[巧沈国辉[],69章李刚?大型复杂建筑结构风致效应及等效静力风荷载巧究[D?浙江大学博壬学位论文2013.[),]70]鲜荣.双层J.[,廖海黎,李明水马鞍型空间网壳屋盖风振系数非线性模燕叠加法计算[]建筑科学,2008,24(9)1(M8.uma-7avenortA.ClGustloadingfactorsJ.Jolof化eStructralDivisionASCE.196793:1134.[UDp[],,M.Kasperski.Extremewindload出stributionforlinearandnonlineardesign?En.Struct"1992,17巧口]g-14l:2734()73XinzhonChenAhsanKareem.EuivalentStaticWindLoadsonbuildins:newmodelJ,J.of[]g,qg[]Srucurenneerin-ttEi200413010:14251435.gg,,()74XinzhongChen,AhsanKareem.Coupleddmicanalsisandeuivalentstaticwindloadsonbuildins[]ynayqgw-Ji2005-2ihhreedimensionalmodes.J.fSrucureEnneerin1317:1071108.tt[]ottgg,,()75 浙江大学硕击学位论文参考义献75鞠红梅.高层建筑等效静力风荷载分析D.北京交通大学2010.[][],7..超高层建筑风重賴合效应及等效静力风荷载研究2014.[D]浙江大学,[词钟振宇77邹良浩..[],梁枢果,汪大海,熊铁华基于风洞试验的对称截面高层建筑王维等效静力风荷载研究m建20-筑结构学报,1211:2735.,78李玉学.大跨屋盖结构风振响应和等效静力风荷载关键性问题研究D.2010.]北京交通大学,[[]79罗楠.大跨屋盖结构风致振动分析及等效静为风荷载研究[D].西南交通大学,2012.[J80梁枢果,吴海洋郭必武,周向阳.大跨度屋盖结构等效静力风荷载数值计算方法机.华中科技大学学[],报200804-白然科学版:110114.(),,-眶毅..200701:1251巧.巧。顾明,周大跨度屋盖结构等效静力风荷载方法及应用m建筑结构学报,,86放1RJHolm货J.D.Useofeienva山说in化ecovarianceinterationmehodfordeterminationofwind[巧jggt-urna-loadefecteJ.JolofWindEnineerinand虹dustrialAerodynamics.l98313J59370.[]gg,rrofndecn-uolmesJ.DBestR.J.Anaoach化出edeteminationwiloadef化olowrisebildins.巧引H,ppg口]A-Joul.ofWindEngineeringandIndustrialerodnamics1981,7273287.rnay;76 浙江大学硕壬学位论义作者简历作者简历杨悦,男,汉族,1990年11月出生,浙江金华人2008-.092012.06浙江理工大学建筑工程学院,获得工学学去学位 ̄2012-.09今浙江大学结构工程专业,攻读硕壬学位攻读硕壬期间科巧成果1楼文娟杨悦吕中宾锋杨伦.考虑气动阻尼效应的输电线路风偏动态分析方[],,,张少,法[J].振动与冲击(录用)[2]摆文娟,杨悦,卢明,吕中宾.连续多跨输电线路动态风偏特征及计算模型研究电[化力建设)(录用77
此文档下载收益归作者所有