脉动风载下的输电线路风偏计算研究

脉动风载下的输电线路风偏计算研究

ID:34844978

大小:9.73 MB

页数:90页

时间:2019-03-12

上传者:U-56225
脉动风载下的输电线路风偏计算研究_第1页
脉动风载下的输电线路风偏计算研究_第2页
脉动风载下的输电线路风偏计算研究_第3页
脉动风载下的输电线路风偏计算研究_第4页
脉动风载下的输电线路风偏计算研究_第5页
资源描述:

《脉动风载下的输电线路风偏计算研究》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

分类号:TU973单位代码:1的35密级:学号:21212061.V、硕±学位论文'戀巧;I中文论文题目:脉动风载下的输电线■路风偏计算研究英文论文题目Studyoncaculationofconductorswingingunder:fluctuatingwindloads申请人姓名;悦杨指导教师:楼义娟教授专业名称:结构工程研究方向:结构风工程所在学院一:津筑工棉学院,,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 渐江大学巧壬学位论文AbstractAbstractInourcoxmtrythehenomenonofimbalancedowerresourcesisserious.Asthemaor,ppjowertransfermethodTransmissionsstemdevelosdramaticallbutitisseverelthreatenedp,yp乂ybflashoveraccidentswhichwerecausedbconductorswinin.Thereforetheresearchonyygg,accurateandeficientcomputationmethodsofthedisplacementofinsulatorhasgreatracticalpsncancenigi行inengineerig.Consideredtherelativevelocityofconductorandwind,过methodtocalculateconductor-swinging化atafectedbyaerodynamicdampingwasestablished.Then化etimedependingfluctuatingwi打dloadswouldact〇打tra打smissio打lineweresimulatedwiththecombinationofuas-harmonicwavesuperstitionmethodandqisteadyassumption.Anonlineardynamic-transmissionlinemodelwhichwasconsistedofeightspanelectricalconductorswasestablishedpreciselybyusingfiniteelementmethod.Furthermore,theimpactofaerodynamicdampingcausedbconductormovementsondnamicresonseofthetransmissionlinewasdiscussedinyypdetail.Itisshownthattheaerodnamicdampingcanreducethemaximumvalueofconductoryswininsinificantlwhe打calculatedislacementofinsulatorbdnamicanalsismusttakegggy,pyyytheaerodynamicdampingintoconcern.Ontheotherhand,theaerodynamicdampinggainedwiththeUStransmissionlinedesigncodeASCENo.74wasadopted化computethe姐splacementofinsulator.Theresultoftwomethods化calculatetheresponseofinsulatorunderaerodynamicdampingwerecompared.Itshows化at:化eresultsarecloseU)eacho化er,aerodynamicdampingainedwithASCENo.74canbeusedinfutureresearch.gComared-thediferencebetweenMultisanmodalandsinleendulummodalppgp,investigatetheaccuracyofsinleendulummodal.ThrouhdiscussedthediferencebetweengpgmaximalvalueandthemeanvalueofconductorswingingcalculationresuUsthisthesis,researched化e-imactofturbule打twindvelocityonthewindinduceddnamicamlificationpypefect.虹orderU)have泣beterunderstandi打gofwindloadcalculationinthecomutationofpconductorswingingthesetinofcalculationarameterswerecomaredamonfourdiferent,gppgationaesncodes-ldi:《GB50545》、《ASCE、NgNo.74》《lEC608262003》and《BS8100》.ermoreasedon-Furth,barealeightsantransmissionlinethediferencesamonthecalculationp,gresultsofdiferentcodeswereanalyzedunderthesamedesigncondition.Theresultshown化at:口sing化esinleendulummodaltocalculate化erotationanleofsusensio打insulatorstringisgpgpsafeThednamc?imanificationfactorisabout1.31.5windloaddefinedinChinesedesin;yg;gcodeisthesmallestthemai打causeliesinthele巧ofcomrehensiveconsiderationofwind,ppressureasymmetricalcoeficientanddynamicmagnificationfactor.V 浙江大学硕去学位论文Abstractfreuenc-domaAqyinanalysismethodofconductorswingingispresentedinthispaper.Geometricnonlinearityisconsideredwhen也emeanwindloadactingon化ewires,A-Mo虹si打lesantransmissionlinewastakenastheresearchmodalreoverordertoimrovegp,,pcalculationeficienctwenttenfivethreetwomodevibratio打wassuerosed.TheKsu化y,y、、、、ppshownthat:Theresponseofconductorswinginthefreuencdomainandintimedomainisqycoseoeach-ltotherthefreuencdomai打analysismethodhasrovedtobefeasible了aki打the;,qypgfirstandsecondmodelfbrS民SScombination,ahighPrecisionconductorresponsecanbegotwiththebestcalculationefficiency.ThisthesispresentsthefrequencydomaincharacterofinsulatorStri打gsswing,foundout化atwinddynamicresponsepredominantlybackground,and化ereso打ancecanbeignored.Load-Resonse-Correplatio打(LRC)methodwasadopted化calculatetheequivalentstaticwind--loadESWL)ofasinlesantransmissionlinetheresonseofthesinlesantransmission(gp,pgp-linethatactedESWLorthetimedependingfluctuatinwindloadswerecomared.Theresultgpshownthat;TheESWLofconductorthatbasedo打theLRCmethodisrecisel.pyKeywords:Conductorswinging;Aerodynamicdamping;Transmissionlinestandard;FrequencyDomainMethod;EuivalentstaticwindloadqVI 浙江大学硕±学位论文目次^I摘要IIIAbstractV0次i第1章绪论11.1研究背景11丄1我国输电线路的发展状况11丄2风偏闪络对输电线路的影响21.2研究现状31.2.1国外学者对于导线风偏研巧31.2.2国内学者采用静力方法分析风偏5123..国内学者采用动力方法分析风偏51.2.4国内外输电线路规范关于风偏计算风荷载取值的差别61.3本文主要工作71.3.1考虑气动阻尼效应的输电线路动力风偏计算71.3.2我国规范风偏计算方法及与多国规范对比71.3.3利用频域法计算导线风偏8134..导线的等效静力风荷载8第2章考虑气动阻尼效应的输电线路风偏动态计算92.1引百92.2考虑气动阻尼的动力学公式923有限元模型的建立112.3.1八跨线路的原始参数112.3.2有限元模型建立口21.3.2输电线路动力特性32.4风速的模拟142.4.1平均风速模拟152.4.2脉动贈模拟152.4.3风速模拟结果172.5有限元结果分析182乂采用ASCENo.74所提供气动阻尼202乂.1ASCENo.74中的气动阻尼公式20I 浙江大学硕壬学位论文262SCENo74..采用A.气动阻尼有限元结果212.7本章小结22第3章我国规范风偏计算方法及与多国规范对比243.1引言243225.风偏计算模型讨论3.2.1弦多边橫型253.2.2刚体直棒模型263.2.3多自由度连续多跨模型273.2.4单摆模型正确性讨论273.3动力放大系数取值283.4#国规范线路风偏风精载的对比293.4.1输电线路风前载计算么式293.4.2巧载系数/313.4.3风速的取值313.4.3体型系数C37d3.4.4巧系数对比373.4.5八跨模型按多国规范计算风荷载值403.5本章小结42第4章利用频域法计算导线风偏444.1引言444.2基于频域算法的导线风偏响应基本理论454.2.1频域法计算导线风偏响应的过程454.2.2导线使用频域算法的条件讨论484.3单跨导线算例494.3.1单跨导线的有限元建模494.3.2导线模型施加平均风荷载504.3.3频域法计算结果544.4本章小结57第5章导线的等效静力风荷载595.1引言595.2风致响应分类605.3导线风偏响应在频域上的表现615.4导线风致响应的计算方法及其对应的等效静力风荷载抗5.4.1平均响应625.4.2背景响应树 浙江大学硕壬学位论文目次564.4.3总响应565.5导线等效静力风荷载求解5-6本章小结67第6章结论与展望的6.1本文主要结论的6一.1进步王作展望71参考文献72作者简历77攻读硕壬期间科研成果77ill 浙江大学硕壬学位论文第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玄面:alumpedmassfswindpressureforceTn:tesionSs色panlength图-Tsumoto13ji所提出的导线风偏计算模型民化保证悬垂串导线与巧塔之间的最小安全间隙为设计准则,利用图解法计算了导线风偏角。一(X在风偏规范计算么式中,存在个风压不均勾系数。导线是大跨度长距离的结构,因此在一定范围内导线的平均风压并不能处处相等a就是考虑了这个因。风压不均勾系数素后对导线平均风压的折减系数。国外学者对风压不均勾系数也做了许多研究:,"C1ha/对于不同档距的单跨输电线路进行了研究,讨论了线路风荷载的巧减方法,并。指出随着风场嫌流度及档距的增加,风压不均匀系数会相应减小4[]-M-Hiratsuka对台风影响下多跨输电线路iyakoimatestline(如图14)进行了实测,j获得了lOmin的平均风速,并观察其风偏响应的结果,并利用风洞试验,研究了时距的长、短对于风压不均匀系数的影响。.辦1!薩I ̄l4Hratsuka图i所采用的实测线路4 浙江大学巧±学位论文第1章绪论*MO?〇三IIrmS11III11 ̄ ̄ ̄ ̄ ̄巧 ̄ ̄Z.IxirnIa韦?pfiMSWIO|^I.《广T一—<〇。一言^巧_>i——巧〇ghorizontal香:。'’。—’扣s皇mini![jj^旨一/?-?-'’’一…??-…?i:!i;:肚±扛±:±r》?古i韦禅年f毕I11111*-S'吾受至II!iastesaveraevaseed/s!110mgidpd{)--ratsu图15Hiratsuka实测获得的lOmin平均风速图16Hika得到的风偏角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)在有限元程序中考虑美国输电线路设计规范ASCENo.74《Guidelin故forElectricalTransmissionLineStructuralLoading》所提供的气动阻尼,获得其在时域上的风偏响应,并与本文所述考虑气动阻尼的数值摸拟结果进行对比。验证了ASCENo.74中气动阻尼计算么式的准确性。132..我国规范风偏计算方法及与多国规范对比关于规范算法研究及多图输电线路风偏风荷载计算方法对比,本文主要做了兰部分工作;(1)为考察单摆模型的准确化,首先采用静力单摆模型根据规范给出的导线风荷,载计算公式及相关计算参数的取值方法,求解了输电线路的静态风偏响应,并与采用有限元法求解得到的八跨线路的平均风偏响应进行对化。(2)将有限元模型计算得到的风偏总响应与平均响应进行比较确定了动为放大系,数的取值。7■J 浙江大学硕去学位论文第1章绪论-(3)分别介绍中国输电线路规范《GB505452010》、美国输电线路规范《ASCENo74》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软件,建立包括整个耐张段的多控距输电线路精细化有限元模型,基于非线性瞬态分析方法揭示考察输电线路自身运动引起的气动阻尼对输电线路风偏瞬态响应的作用特征,并与考虑美国输电线路设计规范ASCEN0.74脚《Guidelin放forElectricalTransmissionLineStructuralLoading》给出的气动阻尼么式计算得到的风偏值进行了对比。2.2考虑气动阻尼的动力学公式ssLotiredoSouza的研究表明:结构发生风致振动时结构与来流之间的相对运动效应,会使动为系统的阻尼增大,而新增加的阻尼即为附加于结构上的气动阻尼,该类阻尼具有随平均风速增大而增大的特征。而输电线路具有跨度长、阻尼小W及柔度大的特征,属于9 浙江大学硕击学位论文第2章考虑气■动阻尼效应的輸电线路风偏动态计算典型的对风荷载敏感的非线性结构,风偏位移大,故在求解导线风偏响应时应当将气动阻尼的影响考虑在内。乂K向//扔^扔>1图2-1多自由度导线俯视图2-示的《自由度为详细说明导线气动阻尼的产生机理和作用特征,在此引入如图1所多跨输电线路,在风荷载作用下顺风巧位移向量r可表示为:=-?2-YYYY1[??j()?为第;其中K/点的顺风向位移,相应输电线路的非线性运动方程可W写为MY=-2+C+KYYF()口)^^,j式中M和为线路结构自身的质量矩阵和阻尼矩阵;F和y分别是系统的加速度和^5有关:速度向量,与结构的位移状态;i:7{巧是考虑了几何非线性效应的切线刚度矩阵=Y-KYK+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=-FY+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!lh33!!|!!jI"L丄丄丄1丄」L!l=。*?心228L3261331LM"LM42L393^426L33l图2-2500kV某线路N3 ̄NH耐张段参数示意图?N为研究上文中所提出的气动阻尼的作用500kV实际线路N3il耐张,本文采用某段为研究对象,该耐张段全长2918m.是由两端2个耐张巧塔和中间7个直线巧塔组成的8档线路,处于B类地貌情况下,整跨巧无转角。各跨的挡距、呼高如图4所示,年平均运2506kN4XLGJ-4005-1行张力为.。导线型号为/3,其物理参数见表2。II 浙江大学硕壬学位论文第2幸考虑气■动阻尼效应的输电线路风偏动态计算2-表1LGJ400/35钢志错绞线物理参数 ̄ ̄ ̄ ̄ ̄7ZZZ2弹性系数,mm化命An?计算截面/线密度/kg/kmmm/Mpa425.2426.82650001349---mm直线塔的悬挂绝缘子串信号为FXBW500/180伞颈为190160110患棒材料为,,环氧树脂玻璃纤维引棒,结高4.7m。耐张塔与直线塔上的绝缘子串的物理参数料如下表2-2所示。.表2-2耐张段上绝缘子串物理参数长度/量/位置泊松比fmkg耐张0257—641140绝缘塔子串直f0.3534^2.3.2有限元樸型建立输电线路处于静力平衡状态下时,受到了运行水平张为W及自重的作用。在受到这两种力的作用下,导线精确的几何形态可采用悬链线方程描述。=竺-y--CO化(210acodia)[]()qL式中:y为计算点垂度,义为计算点到挂点的水平距离H是年运行水平张力;为,q'i=+导线单位长度质量;asinh[/?(C/L)/sinhy0]y(?;L为跨距,C是导线挂点高2H差。悬链线虽能准确表示导线在自然状态下状态,但计算较为复杂,,为便于分析当导线^跨中垂度跨距I之比小于化1时导线的形状可1乂近似看成拋物线与实际情况悬链线/与,fwi相差很小,采用公式:4.z王-巧/(y+三-=义口n)2LL:式中,/为跨中垂度,采用下式计算=-。/瓷(2)通过式(2--11)及(212)即可确定出W抛物线表示的导线受重为荷载悬挂情况。考虑计算效率及精确性的问题,导线的每个单元长10m。整个耐张段导线部分共采用了30012 浙江大学巧主学位论文第2章考虑^动阻尼效应的输电线路风偏动态计算个节点。这样就得到了有限元各节点的初始坐标。在导线发生风偏闪络时主要是发生顺风向的位移,输电塔风致振动响应对导线风偏的レレ影响不大,基本可乂忽略。考虑到计算分析效率乂及计算机性能条件,在建立有限元模型时不考虑输电杆塔的影响。建模时将绝缘子串和杆塔之间的连接作为固定较接。虽然在大一风情况下,多分裂导线的各子导线存在定的不同步运动或次档距振荡,但主要仍W整体运动为主一,因此在研究绝缘子串的风偏闪络问题时,将四分裂导线等效为根单导线。输电导线属于非线性特征显著的高柔度结构,因此在现实情况下仅能受拉而不能受压且无法承受弯矩,故本文采用ANSYS有限元软件提供的LinklO单元模拟;而对于风偏发生时始终受拉的绝缘子串来说,可采用Link8单元模拟,其重量W附加质量的形式融合到有限元模型中,具体实施时采用Mass21单元模拟。考虑到导线发生风偏时,绝缘子串的运动syP7整体摆动为主f],若将其简化为刚度较大的杆,结果仍然是合理的故悬垂绝缘子串,由两个节点简化而成。整个耐张段绝缘子与导线部分共由307个节点所组成通过材料将,这些节点相连,形成单元。在导线的有限元模型上施加重为荷载,并于ANSYS找形分析一后更新坐标,差看找型前后两次坐标的差值,若差值较大还需进行下次找型。导线初始张力取年平均运行张力,为25.06kN,结构阻尼取化005。通过上述过程所建立的有限元模型如下图所示:,/户■.f.2- ̄图3N3N11耐张段有限元模型当输电线路受风荷载作用时,导线轴向应变很小但位移很大,对其进行分析时需要考虑几何非线性效应。本文通过激活ANSYS有限元软件中的大变形和应为刚化选项来考虑该效应的影响。2.3.3输电线路动为特性多跨线路的前六阶振型和频率如图2-4所示L乂看出多跨线路的振型主要W平面外,可摆动为主,且对应的频率较低。同时还可发现,在直线塔绝缘子串的连接作用下,相邻跨13 浙江大学硕壬学位论文第2章考虑气■动阻尼效应的输电线路风偏动态计算。另外,其导线之间的振动模式存在显著的传递效应,因多档距输电线路平面外刚度较低2031Hz前阶频率均低于化。—二0.1338細Zb0(a阶阶.1409加Z)()_0Hd四阶0.15770HzC£阶.14662z()()e五阶0.16345H.1629姐z0z()识六阶图2>4输电线路前六阶振型2.4风速的模拟风荷载模拟准确与否是直接关系到风偏响应计算结果是否合理可靠的关键。风速时程!^可乂分为平均风速和脉动风速两部分。平均风速通常被视作不变量,采用沿高度变化的指数或对数规律来描述,基;脉动风速则采用平稳高斯随机过程模拟。在风速模拟完成之后:于准定常假设,可将风速转换为作用于结构表面的风荷载-F=DLCU+213)y非nD(:W(式中,Fy是作用在导线上的时变风荷载,戶为空气密度,D为导线等效迎风直径,I。为加载点控制长度C为阻为系数C7和fi分别为Z高度处的平均风速和脉动风速;,d,,,41 浙江大学硕±学位论文第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输电线路的设计基本风速不宜小于27m/s,故取基本风速为27m/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仍Cohe2-20()()()()yj如)()馬&jj||拓式中Co/i(co)采用Davenport推荐的风场空间相关函数:222---仍C片+'y)+Gzz如巧所*(,*)1/。。f,、,/、V=<=-2cA仍ex>口。()p52Vz其中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?Jj,,\| ̄=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广|;||||^||500100015002000时间(S)图2-5挂点4处的风速时程图2-6为脉动风速功率谱的目标值和模拟值的对比。可W发现,风功率谱的目标值和横拟值吻合较好,说明获得的风速时程能够有效反应脉动风速能量在频域内的分布特征。’10—特征点白功宰谱fJ-^-、-2*'〇!〇1〇10!0H頻率(z)图2-6特征点处模拟功率谱与目标功率谱比较17 浙江大学硕壬学位论文第2章考虑气动陌尼效店的输电线路风偏动态计算2-7-图为按式口16得出的风场浦流度模拟值和理论値的对比。不难发现,风场浦流)度的模拟值与理论值吻合较好,说明采用修正后的Davenport谱模拟得出的风速可有效反映风场瑞流度沿高度的变化规律。■016.?-0.模m辅流强度15I-.5.巧巧()1平均风削面0-.14,货--013"强奴*..襲。-难穀狂心12..产节--X01,0.'〇20354041525305■莉度(B)图2-7模拟风场浦流度与理论瑞流剖面的比较2.5有限元结果分析取风荷载与导线夹°角为〇,将谐波叠加法模拟得到的脉动风荷载时程加载至有限元模ewmark-型的各个节点上,,采用无条件稳定的N目法对非线性动力方程直接积分求解并ew-运用NtonRaphson法对毎个时间步末尾的位移进行迭代。同时为考察气动阻尼的影,响,分别采用考虑和不考虑气动阻尼两种瞬态求解方式计算多档距输电线路的风偏响应。当采用考虑气动阻尼瞬态求解时一一荷载步时导,取前荷载步作用后导线的运动速度与后-5线所受风速的相对速度,代入式(2)进行计算。止其与橫担过近闪络-。图28给风偏主要是考虑绝缘子串的位移,防,发生放电出了位于线路中间的3、4、5号挂点的风偏水平位移时程计算结果。不考皮側R1;己不々巧气动阻旧I——-6,考阮动阻尼,考由气动阻巧6-55*鸭蜘师4-寒曇予-哪想麵補義-巧3径3I''22-票罢II11000400800120016002002400K0(J2001600200024抑41时间/S时间/s挂点3挂点518 浙江大学硕壬学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算不考虑.气动阻尼—考点与动阻尼6「占'\片2-U^4008001200化0020002400时间/s挂点4图2-8互处挂点的风偏水平响拉时程-在获得绝缘子串水平位移后,根据式(224即可得到风偏角)水平位移"I、。=//IW02-26sin()绝…缘子气山串长I,度一()"动助错邑90r不考虑側KIMQO重5i—考度初as—I…考這气动阻尼巧.15抑曲二:攀卿幽幽耐0唯叫鬟30.400800120016002000400800120016002000时间(s)时间(s)挂点3挂点5不考虑气动阻尼I—考底气动阻尼73.井剛城卿戰i巧30^15-U.?■??I■I■I-Q/400800120016002000时间(s)挂点4图2-9王处挂点的风偏角度时程19 浙江大学硕壬学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算--从困28、图29不难看出,,由于平均响应是由平均风荷载所产生因此是否考虑气动阻尼对于响应的平均值无影响,而气动阻尼对动态风偏响应幅值的降低作用十分显著。2-1因此,忽略气动阻尼的影响势必高估线路的风偏角。图0给出了各绝缘子串挂点处的风偏角极值,相应计算方法为:'-=+-a>^^(227q/)■式中均值一疗为,〇为均方根,是个保证因子,称为峰值因子,为化到要求的安全保证率98.61%,选取峰值因子2.2。——一60rI不考巧气动阻尼々巧气动阻记I■40-30^*12345677123456绝缘子挂点号绝缘子挂点号(a)各挂点风偏角平巧值(b)各挂点风偏角极大值图2-10各绝缘子挂点处风偏角统计值一-分析图210可知,方面,靠近耐张塔挂点处的绝缘子串风偏角比在整个耐张段中部挂点处风偏角要小,这是由于支座约宋效应的影响,使耐张段两端附近绝缘子的风偏角一一较小,这结论在现实情况下普遍存在,在考虑气动;另外方面还可看出阻尼后,风偏角极值大幅降低,说明气动阻尼对风偏响应有着较为显著的折减作用。因此在进行脉动风载下的动态风偏计算时有必耍考虑导线与来流的相对运动所引起的气动阻尼。2.6采用ASCENo.74所提供气动阻尼21AS.6.CENo.74中的气动阻尼么式美園输电线路设计规范ASCENo.74《GuidelinesforElectricalTransmissionLine工程师协会所出版一StructuralLoading》是由美国主木。作为世界上十分重要的本输电线路规范,其中的许多公式及参数取值被他围大量借鉴与引周。在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;因此在采用,ASCEN0.74提供的气动阻尼计算导线风偏时,在非浅性运动方程的阻尼矩阵中直接考虑=+--1)气动阻尼的影响,根据3、32)。,即^是兵式(2式(2求得相应的常系数口,片通过在ANSYS中设置相应的〇,片系数巧可使得在ANSYS模型中的阻尼值为所求得的总阻尼。22ASGENo4.6.采用.7气动阻尼有限元结果根据上文所述采用ASCEN0.74中所提供的气动阻尼与结构阻尼之和作为有限元摸型的总阻尼,进行计算,将计算结果与本章所述考虑气动阻尼的方法迸行对比,如下图所示。21 浙江大学硕壬学位论文第2章考虑气动阻尼效应的输电线路风偏动态计算本文义法80rASCENo.7470-厂"'3〇^|:7苗20-UA10j■■II■I.I./I.I4008001200160020002400时间/s2-1图1挂点4风偏角时程-ASCE图211给出了挂点4位置分别采用本章所述的考虑气动阻尼的方法与考虑L.74所给出的气动阻尼所计算得到的有限元时程结果。从图中可乂看出No,采用这两种方法所得的结果十分接近,整个时程过程基本吻合。2-22-27在图1中给出按式()所计算得到的300个导线节点的水平、竖直位移的极值统计值。-^^■本文义法*3旅1--ASCENo.74?)IASCRNo.74|?I期-i哀0?0000000巧000〇501001502002日03005011523节点巧号节点巧号(a)水平位移极值(b)竖直位移极值图2-12所有节点的风偏响应统计值由这两幅图中化可W看出,本文提出的考虑相对运动速度产生的气动阻尼与ASCE-28ASCENo.74给出的么式(2)计算得到的气动阻尼对于导线的作用十分接近,说明No.74中气动阻尼计算方法可供之后风偏研究时采用。2.7本章小结本文W500kV韦塔八档输电线路为研究对象,建立了考虑气动阻尼的导线风偏非线性动为方程。通过有限元软件,采用非线性瞬态分析方法求解了输电线路的动态风偏响应,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=11'--1,\W+0.5各+-v,斗1各K(,玄32)。、、J-y合VIf^f^V+.W505++0.+J各震,*各W各*;[J^;J=S(3J)hi?.上式中:225P为空气密度,取L;4是第i个绝缘子的受风面积;V为第!个绝缘子钢片所受到的风速?。弦多边模型较好地模拟了现实中绝缘子串的变形情况,但也使设计过程增加计算量。在实际情况中,架空输电线路的绝缘子串所受风荷载及其自重相较导线而言是个小值。因此一,般来说,在设计过程中无需将绝缘子串精细化分为若干段进行计算。只有当绝缘子串的自重较大,且需要计算精确的风偏位置时,才会采用弦多边模型。32..2刚体直棒模型由于弦多边形法计算过程较为复杂,当绝缘子串较长、绝缘片数较多时所带来的计算量不容小觀,也给工程设计带来了西难。在实际工程中计算过程较为简单的刚体静力学方法应用范围更为广泛一。刚体静力学模型化称为单摆模型,是将悬挂的绝缘子串看做根受到均布荷载作用的刚体直杆。绝缘子串所受的水平及坚向为作用于刚体直杆的中部;将导线质量和所受的风荷载集中至绝缘子串的末端。并假定当风荷载与导线和绝缘子李自重达到平衡时得出的摆角即为风偏角。////■/么iIi,h-H图3-2输电线路风偏静力学单摆计算模型-20如图3:所示,绝缘子串风偏角可表示为26 浙江大学硕壬学位论文第3章巧国规范风偏计算方法及与多国规范对比-+/2|亞&=01孤()(3-4)W+G2、J式中:&和%分别为绝缘子串和导线所受的风荷载,Gv和Wv分别是绝缘子串和导线的自重。需要指出的是,计算导线风荷载时需使用绝缘子串所处巧塔的水平档距,而确定导线自重时则需采用垂直档距。对于导线所承受的风荷载,输电线路设计规程规定采用下式进行计算:;-Asin0(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,12Q1234567挂点编号图3-3多自由度连续多跨模型与单摆模型风偏角的比较由图3-3可W看出,采用单摆模型计算得到的风偏角大于多自由度连续多跨模型的风偏角,尤其是对于与终端塔相邻的挂点,单摆模型计算值偏大较多。这是由于多自由度连续多跨模型中两端的耐张塔对输电线路面外摆动产生明显的"拉扯’’作用和各按线路间的相互稱合作用,使得其响应值比单摆模型要小。这表明,就模型本身来说,采用单摆模型来求解处于整个耐张段中的绝缘子风偏角是偏于安全的,可L乂适用于工程设计。3.3动为族大系数取值风荷载包含平均风荷载与脉动风荷载。平均风荷载作用于结构时,体现静为的效应,而当导线受到脉动风荷载作用时,所发生的脉动响应体现出动力特性。我国规范在利用式-(35)计算风偏时导线上的风荷载,取风荷载调整系数Pc为1.0,即不考虑导线的脉动效应。这样的简化方法是否可行,有待研究。取第二章中所计算得到的绝缘子串挂点处风偏响应极大值与平均值做比较,就可W计28 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比-算出动力放大系数。将各挂点处风偏响应统计值列于下表31所示:表3-尼时各挂点风偏响应的统计值1考虑气动阻挂点号1234567。34..平均值/.046.255847.450.455635.24.2均方根4.24.34.64.43.5.24。.357极值/43.255.565.560.064.842.8.171角度放大系数1.271.201.211.191.161.22风荷载放大系数1391.31.391.481.441.431.451-表31中给出了所有挂点在瞬态风场下,考虑气动脏尼所得到的风偏角度统计值。在 ̄脉动风场下求得的风偏角约为静态风偏角的1.151.3倍。将得到的风偏角转换为单摆模型进行计算,即可得到对应的水平风荷载极值。可W算得,在考虑脉动后水平风荷载极值约 ̄为平均风荷载的1.31.5倍。说明在实际风偏计算过程中脉动风的动力放大作用不容忽视,规范所取的动力放大系数偏于不安全。3.4多国规范线路风偏风箭载的对比我面输电线路在匙风情况下风偏闪络灾害频频发生,说明规范中风荷载取值方法存在一一定的缺陷。如何安全合理的对线路风荷载进行计算成为个值得研究的问题。对国内外有代表性的规范进行对比有利于学习到国外规范中先进的部分,从而改进我国的输电线路设计规范中有关参数的取值。本章选取的有代表性的输电线路设计规范主要有:中华人民共和国住房和械乡建设部-发行的《n〇KV-750KV架空输电线路设计规范GB5054520》50545(10)(下文简称GB);美国主木工程协会发行的美国输电线路规范《ASCENo.74》(下文简称ASCE74);国际电P91王协会输电线路规范《圧C60826-2003》下文简称旧C60826木工程和建筑();英国上aticeower ̄ar结构标准委员会发布的《Lttsandmaststi:codeofracticeforloadin》下ppg(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加Q3-8)乂典(式中:是荷载保证因子;Cxc是导线的体型系数通常取1;Gc是考虑了导线高度、r,地貌情况与脉动影响的系数、&是考虑了导线长度的档距系数;d是导线直卷;L为导线弦长;q〇为基本风压,可表示为:(-93)上式中:是空气密度,X是空气密度修正因子,K,通过查表求得。是基本风速r为粗粮修正因子,可查表求得。BS8100中,风巧载计算时分为平均风荷载与脉动风荷载分别求解:=R(3-10)再wf巧c〇v(3-11)巧成马W;及=sn-C〇王Ci(3口)CWcc皆。=氏义戸-1〇K。()(33)nd*音式中:马W是平均风荷载;巧W是脉动风荷载;是导线平均高度处的风速;Cc是阻为系数,由雷诺数和导线类型、表面粗稚度、覆冰情况等决定;化是导线直径;k是导线的弦长;Gc是导线阵风响应因子;乂A是地面粗植系数;为风向系数;a是高度相关系数;Ag为有效高度;是荷载因子;为基本风速。W上所述的公式^均1^准定常假定作为基础具有相似的形式将上述么式归纳为,,可一个算式:^^=-FxpvCA?sinw(314)jj^30 浙江大学硕去学位论文第3章巧国规范风偏计算方法及与多国规范对比式中:是各国规范中考虑了风压不均匀性及动力放大因素的综合参数,主要受到高^度。、跨长等影响?为荷载系数,P是空气密度,V为平坦地貌下或山区地貌下的风速,受基本风速及风剖面、山区地貌等因素的影响,Q是体型系数。A是受风面积,W是导线与风向间夹角。当采用不同规范来计算同一耐张段输电线路时所取W是定值、,风荷载计,P算值的差异产生主要是由于V、Cd、n四个数值的取值不同。接下来着重讨论这些参数的取值差别。3.4.2荷载系数尸荷载系数是根据线路的重要程度而确定,用来调整线路安全等级的。国外规范中荷载系数的取值差别直接体现在设计风速的重现期上。我围规范并没有荷载系数的概念,而是在计算荷载设计值时采用结构重要性系数与风荷载分项系数的乘积来代替荷载系数。表3-2各国荷载系数;重现期/年ASCE74IEC60826BS8100GB50545i^V根据图表差蘇'100-1."值选取,其中zih-1501.1对于特另3M黑府风荷载取200-1.3重*要的线〇96 ̄400-.1.25。1.45J^500 ̄y3.4.3风速的取值风速值的大小与基本风速、所处高度、所在地貌类型等内容密切相关。其中基本风速是根据不同高度的测风仪器一、时次时距所测得的最大风速,按定的高度、时距要求统,一50545826BS换算而成的最大风速。ASCE74;GB和圧C60;8100规定基本风速取值分别按3s、lOmin、化时距,平坦开阔地貌离地10m高,重现期按50年确定。由于自然界大气边界层内的风速有脉动特性一般来说速的时距与获得的平均风速成反比,,所取风关系。时距越小,计算得到的平均风速就越大。取不同时距时,基本风速所得值的关系可=W:4::表示为圧C:出1.310.94。在大气边界层内,由于地表粗禮元的摩擦作用,风速大小随离地高度的变化而发生变化:离地高度越高。,所受到的摩擦作用减小,平均风速相应增加31 浙江大学硕壬学位论文第3章我国规范两偏计算方法及与多国规范对比风速的垂直变化不仅受到地面粗粮度的影响,而且逐会受大气稳定度的影响。因此无论利用理论进行推导或者通过实际观测结果进行分析,想获得其变化规律都是十分复杂的。目前,主要用对数律和指数律来描述。指数规律和对数规律差别并不大,所目前围内外都倾向于用计算比较简单的指数曲线来表示风速沿高度的变化规律。GB50545、ASCE74、BS8100中都是直接使用幕指函数来表这风速的垂直变化规律。GB50545中:指数型风速剖面可表示为:(ZXU二-。U(315):、寸式中,a为平均风剖面指数,[/、分别为离地高度Z、Z处的平均风速。根提,,,,-粗粮度指数、梯度风高度W及风压间的关系,由式(316)可得平坦地貌下A、B、C、D巧类地貌对应的风压高度变化系数,其中/分别取1.284、1.000、0.544、0.%2a取化口、/?;。0.150.220.30。、、(ZT=_-/i-从(316)°U〇JGB50545四类地貌下风高系数随高度的变化规律如图3-4所示。3.5—A类[B类-3—.0C类0?5.I.III.I...I..I50100150200250300350Z(ra)图3>4GB50545的风高系数取值图ASCE74中,风高系数表示为:扣。—=-ii:21(.0)(317)vZgASCE74中对输电线路建设的区域地貌分为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.....■II.tIII.t50100150200250300350Z(m)3-5图ASCE74的风高系数Kv取值图BS8100中,风高系数表示为:。=-zh0(3-心如()/118)Ac)上式中:J苗是地面粗粮系数,&为风向系数,a是高度相关系数,七为有效高度。BS8100中共将地貌分为I、II、、IV、V五种类型与我国B类地貌。ra,其中第m类地貌等同在不同地貌下:,相关参数的取值如下表所示表3-3BS8100中不同地貌下参数取值地貌类型地面粗植系数高度相关系数有效高度KarhMgI1.200.1250’II1.100.140m1.000.1650IV0.86化192V0.720.2310当不考虑风向的影响时:,可画出风高系数如下图所示33 浙江大学硕+学位论文第3章我国规范风偏计算方法及与多国规范对比25—.I类「I—"类—…类^—IV类-2.0|V类I^.11.1.1.11.1...50100150200250300350Z/m图3-6BS8100的风高系数^1^取值图圧C60826中并未直接给出响应的风压高度变化系数么式,而是给出图来进行差值选取:,Gc包含了高度风压系数与风的脉动影响两个方面的因素3Q ̄".111ITTTI11nfnrn「nnrr11111111111111iII1![「ni[二3之二;■一二二二二D--二-:=Z:::Z:Z::Z:::::::^^::::Z:Z::::Z:ZE:::^2a----::^:一-?---—__ ̄——_——.C——二—"’一—二———二,ir1^:1IIIII三三=E::_IIIIiiirI=二乙---=二--_—==:二:二二4::记_二-_一:三!::三:二-一:.一^三二一^二二二二二lit::::::;'9—'-''二IIIII二III二IIIII二二IlIIIIIIIIII!-IN]IJ]IIIIIJiIJIiJUNNIIiNNIiIIIiIIlLLL1IIUI.610IS2025303SAO45505560hleilYstx>erorxigtivgLjm图3-7lEC60826中G取值图e从图3-4至图3-7可W看出,平埋地貌下各场地类型对应的风高系数均随高度而增加;同一高度下,随着地表粗稚度的增加,风高系数逐渐碱小,所有的规范所确定的风高系数都拥有相同的趋势。若取各规范所规定的标准地貌,即类比于我国规范中B类地貌所对应的地貌,将其风高系数作图:,可表示为34 浙江大学硕±学位论义第3章我国规范风偏计算方法及与多国规范对比"*CB5054528ASCE74^■—BS81(坚J-L2.4^畫\.2^0-.8'1i.i.i.i.1.1.1.1.1,1,125507日1001251501752002252日0275Z/m图3-8标准地貌下风高系数对比由上图可见:当高度较低时,BS8100所求得的风压高度系数较大;随着高度的不断,GB50545所给出的风高系数将大于其余两国增加;相对其余两国,我国风高系数对于’高度改变更为敏感。风速值除基本风速、平坦地貌下风高系数的影响因素外还受到了山区地貌环境的影响。我国荷载规范中利用系数对风速进行修正:1当所取位置为顶部B时:)'=-—-;71+tana(\(319)b[)]l.bH式中ana一ana,t为迎风面侧的坡度,若t大于化3时,取化3;K值对山峰取2.2,14H为山对山坡取.;顶或山坡全高(m);Z是建筑物计算位置离建筑物地面的高度(m);当>=Z2.5//时,取Z2.5//。2-其他部位的修正系数,.,),可按图39所示取A、C处的修正系数为10AB间和BC间的修正系数按线性插值确定。d44.^'4—]、fI+图3-9山峰和山谷示意图ASCE中取地形影响因子也t对山区风荷载直接进行修正,也t的计算公式如下:2K=(-(1+/:/:尤)320):,|23-式中:是考虑了地形特征与最大增速的影响值方法见表34。,取=-A:l;c//L山顶的上下风向距离增加而导致风速改变的系数。,((/J)是考虑离、||35 浙江大学硕壬学位论文第3辛巧国规范风偏计算方法及与多国规范对化心是考虑在局部地形情况下,风速随高度增加而减小的系数;JC是山顶与结构之间的距离;Z为建筑所在地面W上的高度。"":凤纖:淵倘,^h!1。胃而J上&上i當所戸^皆VIj1j%9SWWW巧J风妨醒巧是。纖汾;無a悬崖情况b二维山脊或王维轴对称山丘(())3-困10山区风影响示意图表3 ̄4A参数取值表—K/H/L,(^片山丘形状地貌类别y北?CTSIBCD二1.301..维山脊.4515531.515二维悬崖0.750.850.952.51.54王维对称山.051.1541.51.5&0.951英围规范中考虑山区风的方法是根据输电线路所处山区的位置:山顶、山坡迎风面、(3-2-(-山坡背风面,分别按1)、(322)、323)计算:。-"=己也-巧邸()(321)ZrW"10=严(3-22)^吃《(皆=v^zK-y(323)‘(y10-11其中://与均是地貌指数,/与顺风坡度、地貌分类有关由图3确定。而与/3-=山的高度和地貌类别有关由图12确定。必1);A*hh36 浙江大学巧±学位论文第3牵巧国规范风偏计算方法及与多国规范对比舊曇聽藍聽012I.5t7II,l?’,HHWn,???0AiiiTjMJMI1i1.ViVl\iA,3k-3-n图//的取值图表图312的取值图表从上述介绍可W发现,美国、英画规范中考虑山区地貌影响较为详细。我国规范采用一山区地貌影响因亲较简便的方式,所考虑的,通过个么式来进行计算,参数也较为简单少。3.4.4体型系数CdASCE、旧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浦流积分尺度B4.50.01051.8C7.00.00567.1D^0.0037^ASCE74中阵风响应因子、适用于任何电压等级输电线路,并且该系数与导线的高度控距和地貌条件等因素有关,Gw包含了3s与lOmin的风速转换关系W及风的脉动影响。因此剔除风速转换参数后就可レ乂得到lOmin风速下的;系数。7在C类地貌下-(a00mZ ̄50m(b),图313)取档距S为困定值3,h取10m;中有效50m200m ̄400m高度Zh分别取固定值10m、,档距S取。从图中可W看出求得的/系数/随导线有效高度、档距长度的增加而降低。i-6r61.「-二1.55Z10mh5■.訂.和145、氏'^;IIas'、旨1.3S,..5;1.21020304050200250300350400Z/mS/m,,(a)档距不变(b)有效高度不变图3-13ASCE计算得到的阵风响应因子GwBS=在8100中电线路脉动风荷载表示为巧G中G是导线阵风响应因,输we巧W,其c子,与ASCE74不同的是,BS8100中的阵风因子只表示为脉动部分与平均响应的比值,+G因此0c)代表了系数;的含义。Gc由下式可W求得:738 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比-G广KJC(327)i式中&是导线的长度因子,用来考虑随机风速沿线路方向的空间相关性对导线风荷载的影响,,&是高度因子,考虑了导线平均高度处的奈流度对风荷载的影响这两个值分--14315别按图3、图插值选取。—_支电古占式山—一广—r巧-———————-——0.4frjII!軍?_j—…_——,—了---厂—广―—__—-了卞十下T寸丫…_----下-丄…tl十十十了I??--?-’去史山山1山山山击山Lmt}-图314长度因子?巧JM/W巧扣》"ii:iirr^1TTH"—'^- ̄ ̄-M!irxr!4I-斗灌邦糊寻Ff ̄ ̄ ̄ ̄ ̄h ̄ ̄!iiIIjIIpI\ ̄i ̄ ̄?nS^吏.tJ10去占Iw4)4o^iono300化!■>图3-15高度因子&表3-6给出了在m类地貌下有效高度分别为?400m变化时10m与50m,档距从200m,(1+&)的取值。可L乂看出1+Ge与ASCE求得的Gw有相同的变化规律,均是随着导线有效高度、档距长度的增加而降低,但其计算结果较ASCE要大。3-6(+Gc表按化计算得到的风荷载调整系数1)距^'有效高'200m250m300m350m400m!^^10m1.951.891.851.821.7950m1.811.761.721.的1.66圧C标准中G由图3-7来确定考虑了高度风压系数与地形和高度有关和风的脉动c<)()影响两个方面的因素,因此将Gc除W我国规范计算得到的就能够把风压高度系数分39 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比L;离出来,再乘乂其档距因子&,即得到/系数。表3-7旧C所取的C系数供'有效高^^200m250m300m350m400mm1.671.601.5110.751.72150m1.511.491.441.381.30我国规范中风荷载调整系数化是针对500kV和パ0kV导地线考虑脉动风荷载的放大系数,规定在风偏角计算中不考虑该系数的影响。不均勾系数a考虑了水平控距变化对整〇3-4个耐张段风荷载影响,片与a的乘积就相当于式(1)中;7系数。表3-8GB所取的cxa;/系数P档距200m250m300m350m400mxa08.74.70.67065Pc0.0.ASCElOm有症亩志2-■ASCE50.〇m宵效嵩巧1.一广BSlOro宵效商度■—^—-A?BS50m—A有效風度1?8??---二—化ClOm任巧商度▲—1EC日Om冉效商度,LbC--—?GB50545、I—心-?1■-:.4T亡窥'-^1.2-1.0^0.8^I.■?.T■.06IIIIII■2002加300350400450500档距(S)图3-61各国规范的/系数对比/将各国规范算得的系数作于图3-16上。从图中可W看出各国规范关于7系数的取巧,值趋势均是随着接距的增加而减小,除我国,其他国家规范中都随着有效高度的增加而降低。BS所取值最大,而我国规范由于未考虑风荷载调整系数,使得C系数值最小。346八跨模型按多国规范计算风荷载值..为直观地体现按多国规范所计算出的风荷载值大小差别二章2-2中所示的八,取第图跨输电线路为计算模型:,分别按照各国规范进行求解,其中的计算参数取值如下所示40 浙江大学巧去学位论义第3章我国规范风偏计算方法及与多国规范对比(1美国规范《ASCE74》:是50年重现期对应的重现期前载调整因子取1.0;基)本风速按照3s阵风风速选取,是lOmin平均风速的1.43倍,取38.61m/s;&为风压高度变化系数;地形影响因子&,,线路所在地貌对应于美国规范中的C类地貌,由于不考虑山区影响,地形影响因子取1.0;阵风响应因子Gw是根据跨长、挂点高度、持续风的幕指数表面阻为系数乂及浦流积分尺度右确定,本耐张段各挂点对应取值范围在化64-0.67之间ASCE741。;《》中导线体型系数。取值得注意的是,当取国际单位进行计算时,Q取化613,这样获得的荷载值单位为N。(2)英国规范《BS8100》在计算具体某高度处风速巧勾时,根据线路所在地貌环境特征.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)国际电气规范《圧C60826》中:输电线路体型系、数取1.0,10m高度处基本风速 ̄0取27m/s.98之间G;耐张段各越导线的摇距因子取值&在化%;c是考虑了导线高度与取值范围在2 ̄2地貌情况的风载系数对于本课题所研究的耐张段.2.38;,,其之间考虑地貌影响的系数&取1.0。(4)我国规范《GB50545》中:基本风速取27m/s,不均匀系数按摆距进行插值,取 ̄1值范围在化630.72之间的体型系数.1虑化动效应荷;导线/%取1;计算风偏时不考,风载调整系数A巧1。图3-17是取W上设计参数,并利用各国规范公式所计算得到的导线上的风荷载值。可W看出,按照多国输电线路风巧载规范计算得到的风荷载在耐张段范围内的变化趋势基同,这是由于上述国家的输电线路风荷载么式均是采用了准定常假定,且在计算设计风速时都使用了指数率的风剖面。但是由于各种参数取值的不同,使得按照不同国家输电线路一定程度的差异设计规范计算得到的风祗载值存在。41 浙江大学硕壬学位论文第3章我国规范风偏计算方法又与多国规范对化-*—SCE74IA55-^BSSIOO「■-^lEC-50-^GB50545|[.....I.I?2QItIIII1234567绝缘子挂点号图3-口规范风荷载的比较美国规范在计算风荷载时虽取3s阵风风速,但其在阵风响应因子Gw中考虑了风速的3slO转换,在计算Gw时除去阵风风速与min平均风速比值斬的平方,按ASCE计算出的风前载值相当于取lOmin平均风速,因此其计算值并非最大。中国、美国、英国和国际1关于体型系数的取值分别为1.1.>(、1.0、0.9和10,可^看出各国所取导线体型系数相差不大,不是导致风荷载差异较大的主要原因。由前文3.4.4中^系数对比可知,BS8100所取的W系数是最大的为IEC、再为ASCE最后是GB50545。这也是造成风荷载大小,其次,差距的主要因素。相比之下按照我国输电线路规范计算得到的风荷载值最小,这是导致我国输电线路风偏设计偏于不安全、闪络事故时有发生的直接原因。3.5本章小结本章着重讨论了我国规范中关于风偏的计算方法,分别从风偏计算模型、动为放大系数取值两方面进行研究,并将我国规范与若干发达国家规范进行对比,得到W下主要结论:""(1)由于多跨导线的拉扯作用,使得其在施加相同风荷载时计算出的风偏响应比单摆模型计算值要小,这表明单摆模型是偏于安全的,可W在工程计算中采用。(2)对比有限元计算模型得到的风偏响应极大值与平均值,可W发现导线受脉动风 ̄ ̄荷载作用时所得风偏角放大约1.151.3倍对应的的水平风荷载放大1.31.5倍。,,我国在风偏计算时取风荷载调整系数为1.0偏于不安全。(3)对多国规范进行了分析和参数对比,并采用八跨模型来计算风荷载,从中可W42 浙江大学硕壬学位论文第3章我国规范风偏计算方法及与多国规范对比得知:我国所取的风速、体系系数与其他国家相差不大,但关于风压不均匀性及动力放大两个因素的综合考虑与他围相比较小,使得我国风荷载计算值偏小,这是导致我国风偏闪络现象频发的主要原因。43 浙江大学硕壬学位论文第4章利用巧域法计算导线风偏第4章利用频域法计算导线风偏4.1引言目前风致响应求解的主要方法有时域法与频域法。时域法的基本思路是将计算时间段一定顺序在所求结构上进行施加巧作用在结构上的随机荷载分为多个荷载步,利用迭,按一一代巧方法,W前荷载步作用后结构响应结果作为后荷载步开始计算的初始条件,求解一个荷载步后的结构响应r-施加毎。常用的时域方法有iVewmak片法、线性叠加法、■-0法-W沁on、Runekut法及Houbolt法等时域过程计算适用于各种形式的结构g,可W考虑结构所有的非线性因素,且能获得响应在整个时程过程中的表现,表达直观,易于理解,且计算值准确,但计算过程巧杂,消耗时间长。频域法常用的方法有模态叠加法、虚揪激励法等等。当采用模态叠加法计算时,先利用广义模态主坐标将结构的基本动为方程进行转换,再根据结构不同模态振型之间的与质量矩阵和刚度矩阵的正交性对动力方程进行解賴,就能得到各阶模态主坐标的振动方程,将动力方程进行傅里叶变化将其转为频域,通过引入频响函数对目标风致响应进行求解。在频域算法时忽略各模态之间的稻合效应时(SquareRootofthesumofSquares),会使得计算量较少,但这样计算的精度无法保证。若想得到精确解则需要考虑各阶模态之间的相关性(ComleteuadraticCombination)pQ,采用这样的计算方法有较大的计算量。为了在保证求解精度的前提下提高计算效率、减小工作量林家浩亚辉提出了虚拟激励法(PseudoExcitationMethod),这种方法,,张。在计算中自动包含各振型相关性,具有理想的计算效率总体来说频域方法相对时域法具有简单、快捷的特点,在计算高层、大跨等结构的风致响应计算中已得到了广泛的应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==上式中:YfZ.Z...、Zz..分别代表当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阻尼:=CaM+0K(4.4)y式中;a、P分别是假设的与质量和刚度的比例系数。则可利用不同模态与质量矩阵和刚度矩阵的正交性质可W得到:,=兰-m+.’心c+巧0的(45)成)A的W0、义作,}{},==式中:W.+是结构阶广义质量C.4;?()7,/为/阶广义阻尼,j完巧媛巧媛j戸冷_=112=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度响应均方值时,采用CQComptC)法可W得到精确的解,表示(CQ为:=+4-芝玄4知W八(10)屋芝粉勺=^71k\J*k==口。+。s4-z(11)1么倫i艺夺啦扣而机(■y,)l=^7\1j*kj式中;n为选取的参与模态数;是广义位移相关系数;为广义加速度相关系数,表达如下:CR(/)骑鸟(/)别g々机々r-(412)片、。9八w)K(/(乃心/)皆衣lW邮CRe[(2死/T好,(/倘(/馬(/)斯]s,:*J7;(413)*、。2。。4425八批/)(/)么(/)命獻/)(/)(/)新片,防|姑在|在‘式中;//,(/)表示导线k阶机械导纳函数的共扼复数;S&(/)代表第j阶与第k阶广义风巧载互谱。若忽略模态间的相关性,可W提高计算效李,即采用SRSSSuareRootof化eSumof(qSquares)法来计算,可表示为:口=口s=-护z(414)\么\(y,)恥47 浙江大学硕壬学位论文第4章利用频域法计算导线风偏22==Z(_。41^垃知诉窗y,)■H422..导线使用频域算法的条件讨论-上述内容从式(42)开始都需要结构满足线性要求才能继续进行计算。当采用频域法求解非线性结构风致响应时,并假设其发生小位移,若结构在受荷作用下只发生小位移twi过程中结构刚度、阻尼矩阵不发生改变,即可将模型视作线性结构。导线受风荷载产生风偏过程中位移值大,因此无法直接采用频域法进行计算。风荷载由平均风荷载和脉动风荷载组成一般来说平均风荷载会比脉动风荷载大许多。导线在受,_到平均风荷载作用后产生的风偏位移占总位移的绝大部分,因此可政认为仅受到脉动风荷载影响的导线风偏只是一个小位移的过程二。取第章中所提到的八跨输电线路有限元模型-2绝缘子挂点处时程结果。对其取平均值与最大值,将结果列于图4。5.2.1二嚮滅!LI,[1/\1^/^?.I.I.1■I.I.I.I.I.I.I.I,I,I.2*612345675671234绝缘子挂点号绝缘子挂点号(a)平均风荷载引起的y向风偏值(b)总风荷载所引起的y向风偏值图4-2风偏值的对比-2-(图4(a)给出总响应的平均值,就是平均风所产生的平均响应。图42b)中给出的是各绝缘子串挂点处的风偏位移极值,是取98.61%安全保证李所得到的。将平均风荷-载所引起的风偏值占总风偏值的百分比数据列于表41所示:48 浙江大学硕壬学位论文第4章利用频域法计算导线风偏4-表1平均风引起的风偏值占总风偏值比例 ̄ ̄ ̄ ̄ ̄ ̄挂点编平均风总风偏值平均风引起的号引起风/m风偏占总风偏偏值/m值比例/%2.753.4581挂点23.574.化86挂点34.094.6090挂点43.6342686,挂点53804.3888.挂点64.084.579072.843.4284由此可见,平均风荷载所引起的响应占总响应的绝大部分。在现实情况下(即考虑气动阻尼情况)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轴端点10038.000.0038.00..1号10036790.92371420.03583号2.661.7氏34330034.巧2.巧35.59号4号40033.593.3434.895号50032.664.0434.246抓031.794.68巧.64号70030.7.99527巧.08号8号80030.265.抑拍.589号90029.596.2732.口.10号100028.9967031.72...11号11002845707313712.巧.号120027987.310713号130027.587.6630.8114号140027.247.8730.6115150026.号.9780330.45160026..16号.778143034...17号17002669820302918号176.50化.村8.2130.2819号180026.678.2130.2820190026.748..173032号21号200026.878.0730.4111亏"0027.087.9230.5523号220027.367.7230.7424230027..7174730.98号25号240028.147.1731.2726号250028.636.8131.61260029.20627号.4032.00.855.巧号270029.94324429号280030.565.4332.93302撕031.354.8633.47号300032..31号.21423340532号310033.153J634.69巧号320034.162.拍35.3834号巧0035.242.似36.1135号340036.:391.1936.90端点23巧038.000.0038.0051 浙江大学硕击学位论文第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型的频率/Hz1阶0.1973801.68552阶0.393260.巧5923阶0.巧490.巧7254阶0.4M250.440815阶0.巧3190.5065960.655710912阶.57079196067641阶..8阶0.792810.677089阶0.993树0.8485919085722…阶.000.11阶1.19651.021912阶1.19711.022313阶1.401711971.14阶1.40:371.199315阶1.61051.375516阶1.6111.巧5814917阶.82071.55118阶1.8214.5557192037517401阶..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.ntcezwou?i第3阶模态第4阶模态mNsitwi?*?*巧"*>.<SJ7uno.mur£.‘?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巧法盈吝nlu**?CQC法巧合°-‘*^1时域?1^到的掠准差ioj/\i:::/\■.IIIIII.It.I.I...Q■Q0510152025303540节点编号图4-6时域法与频域法计算得到的位移均方根对比4-时域法计算得到巧标准差是由脉动风荷载所引起17)计算得到,含义,可由下式(与由颇域法计算得封的均方根相同,两者具有可比性。2-平均响应时域法响应标准差总响应)WW-(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■II.I.t.I.I.I.I.0己10152025303540节点编号图4-7时域法与頻域法计算得到的位移极值对比SRSS巧巧合SRSS法由备—*-eCQC法巧合5.Or-?-IIl法组合hr<ii\,U..5-I。./X苗./VW■.I.i‘-■!■?I.I.I()〇..?I代?I-|??玄m.!I.I■I051015如巧3035400510巧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臣ISnIIi議^-.**。.-、*1、0*/A苗.5?????》?%/_"0I.I.I.I.II.I.1■t?0Q■Q05101日20巧30化400日1015202日30巧10节点巧号节点编号(a)均方根对比化)位移相对差值4-图9取10阶模态进行姐合--!?SSSiSPSS法殖吝__RiCi5‘〇-夺-化7"CQC法巧合rWh〔QC法组合|时试法巧到巧标准差■4.5■\r\*05101520253035400510152【)2日303540节点巧号节点巧号(a)均方根对比(b)位移相对差值图4-组合10取5阶模态进行56 浙江大学硕±学位论文第4章利用频域法计算导线风偏—SBS法坦合—SR班法品n7rII5Or,II-亡贼数■!.4;s,\义繼f:\p..".15|.2/\害芒。'I\;r一:;/i..I,I.I.?.I.1.I.I.IQKI■?-I.I.1.I.IIIQ'Q5阳520巧孤况仰051015如25303540W1节点编号节点编号(a)均方根对比(b)位移相对差值4-图11取3阶模态进行组合〇.’■[lia.巨::\E]户.lAiVy它。..1\::i::'—一 ̄ ̄ ̄ ̄ ̄ ̄ ̄—■ ̄__—?一^_—___________OL<■■,■,■■0.0t_i.I.I.I■一i.I.I.I.I.OtII广I广III05况说4005101520253035405阳巧洗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振分量的等效静为风荷载分别求解,进行组合后得到总响应的等效静力风荷载。兒gwKMtChlenAhsanKareem等学者提出将等效静为风荷载表这为平均分量、LRC法下的背,景分量和惯化力形式表达的共振分量的线性组合,并通过权重因子使等效目标响应与实际风荷载作用下响应极值完全相同。目前学者们研究关注于如何高效率、高精度地求解高层建筑大跨结构等建筑物的等效静力风荷载。作用于导线上的风荷载将传递至输电杆塔一个非常重要的水平荷,是输电塔所受到的载导线的等效静力风荷载研究,还鲜有涉及。这主要是由于导线是明显的非线性结构,,无法在基础模型上直接求解等效静为风荷载。在第四章中曾介绍:考虑导线受平均风荷载作用达到平衡位置后,再施加脉动风荷载的过程可将结构视作线性。因此可W计算出平衡位置后导线模型的等效静力风荷载再与平均风荷载组合,得到总响应的等效风荷载。本章首先介绍了结构的风致响应类别^及对应的等效静为风荷载计算方法,再对导线风偏的过程在频域上进行分析,发现风偏响应主要由背景分量所控制。然后采用LRC计算处于平衡位置的单跨导线模型等效静为风荷载。将等效目标响应与施加气动力模型的响应进行对比,两者十分接近。这表明采用LRC法来计算导线等效静为风荷载切实可行。59 浙江大学硕壬学位论文第5章导线的等效静力凤荷载5.2风致响应分类导线所受的自然风中包含着平均风与脉动风,任意t时刻结构Z高度处所受到的风速可化写为;vzt=vz+vzt-(51),,()()()其中云Z代表平均风速,其数值不随着时间的变化而改变;代表脉动风速,不();罔时刻的脉动风速值不相同:。根据准定常假定,所引起的风荷载值可表示为2=Cyzfdp(,)去=zClvzvzt-+t+vz(52)^p\{,,{)){){)^=tFz+Fz{{,))式中222,歹2是导线在高度处由平均风引起的平均风压;歹是导线在高度处的,()(〇一般情况下,脉动风速小于平均风速云Z脉动风压,p为空气密度,Cd为体型系数。,()故脉动风速的平方项将远远小于平均风速云Z的平方项W及它们的交叉项。因此在()2计算脉动风压时,通常忽略项。7;2<。导线2高度处节点的脉动风压可表示为:,/(,)^=-FztCvzvzt(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与动阻尼4035--30\2-5S20-I-1m15叫-hlil/LyulLIA.IfjiiAJlAM一一>一j0.020.040.060.080.100.120.14频率(Hz)(a)背景响应61 浙江大学硕±学位论文第5章导线的等欢静力风荷载—不考虑……考揭动r-1401-120I^1W1ICIIAiJLJikiw^Ii....L,-.,.,\Q0150)0.250.300.040i.公巧.频率(化)(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性结构。导线的背景响应可由方差分析(covarianceanalysis法来进行求解。根据结构)力学知识,当导线在脉动风巧载片作用下;C位置处的脉动响应值可表边为:,=-f^义,乂55而批()()(从(而]此处J(X;C是施加平均风荷载后的导浅模型的影响线系数,那么脉动响应的平方d0,)为:2f义=护义义如-56A2,,(。)(l而l而2如())(K(K()脉动风具有随机性:,因此考虑其统计值才有意义。响应的均方差可W表示为2f=护X护义-7兩,|兩l而l2叫口)()(K,AK而J)((()因此响应的均方根为:"=_C欠义义JX义■*义瓜把5B而,f,rl,巧.()J卢(12)(。I)(02)12({式中:为导线;C、X,两点上风荷载的协方差。从上式可W看出,背景响应i与结构的动为特性无关。Kasperski为了求解刚性低层房屋的等效静力风荷载,提出了LRC(Load-Resonse-Correlation)法LRCp,低矮刚性房屋的风致响应主要由背景响应所控制,法就是求解背景响应的等效静为风荷载。受平均风荷载影响下的导线结构可视为线性,因此其背景响应的等效静力风荷载也可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巧*(引,WW式中-2:说为背景响应的保证因子。那么(51)可w写为:^^=^gerB^.^d^例句,(0(02VVS22),()1从式5-H)不难看出FX(,rs是背景响应的等效风荷载值与发生.()就,其响应的位置及类型密切相关。5.4.3总响应忽略了共振响应导线的总响应由平均响应与背景响应组成:,可W表示为,=-rr+(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化凤拉掛I2-.4f2-12-.2/?2〇??'??记??"??/??????巧勺8\.、8-拽¥.\/16\[6/::IV/八1..巧;'.V2:::/:Vy、--00.8I1I>II■I■I■I■I.I■I..■■I■I...If05101日202日3035400日10152025303540节点编号节点编号(a)沿导线的风荷载分布化)两者相对差值65 浙江大学巧去学位论文第5章导线的等效静力风荷载0.3〇1I<二I^齡I0-.260屠言Z0-.22化-专160-.14-0.12J1■??'■'0—?l〇0.20.40.60.81.0x/L(c)导线上的分布风载图5*4等效静力风荷载与规范算得风荷载的对比.-4给图5出了等效静力风荷载与采用规范算法得到的导线风荷载,从图(a)中可W看出,,两个风荷载在沿导线跨度方向具有相同的分布,跨中位置处风荷载有明显的降低这是由于在导线建模时,跨中附近位置有若干个节点的间距较小,相应的节点控制长度较短,,导致这些节点位置处的风荷载值较小将导线上节点风荷载转换为分布风载后(图""一(C)),,就可W发现,采用规范算法算得的导线分布风载是条光滑的曲线总的等效""静为风荷载中考虑了背景响应的等效静为风荷载,因此并不光滑;当不考虑不均勾系数时,规范算得的导线风荷载与等效静为风荷载十分接近,两者最大的相对差值约为10%。将计算得到的等效静力风荷载施加到单跨导线原始模型上,得到等效目标响应,并与实际动为风荷载作用下的风偏响应统计值进行对比:,结果如下所示等效静力风荷我作用■—I实际請方风荷载作用-10I.I.I■I.I.I.I.I.I.0510152025303540节点编号图5-5各节点水平位移响应66 浙江大学硕击学位论文第5章导线的等效静力风荷載65-.■60-./55-§.\f■、'-E四、■/-■45寡./妇\-■/4、.0■和■■■若■■■■■-3.51tIIIIIII....■■.300510152025303540节点编号图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]uhsofa-uzak巧aleldlowwureconduc4Hiratsuka8MatsiYFukudaN?巧testresindresstorC.[],,,p[]Proceedna-ingsofIEEERegion10InternatiolConferenceon.IEEE,2001,2:664668.l5StClairJGClearancecalculationsofconductorsbuidi打sC.TransmissionandDistribution]g[][Conference996*IEEE1:4巧498.,,6AllenL.Calculationofhorizontaldislacementofconductorsunderwindloadintowardbuildingsand[]pgrtittuCrsPresented-thsrucres.Paat化e37化AlConference.IEEE1993Al/lAinO.o舟suppong[]pennua,:[7]TsujimotoK,YoshiokaO,OkumuraT,etal.InvestigationofconductorswingingbywindandUsapplicationfordesinofcompacttransmissionlineJ.PowerAaratusandSstemsIEEETransactionsg[]ppy,on-198211:43614369.,()8RikhVN.Conductorscingsintransmissionlinesandeffectoflonsanswithsteeslopesinhill[]pagppyterranJourna-iJ.IEIl2004855:816.[](),,()9ChayMT,AlbermaniFawesH.Windloadsontransmissionlinestructuresinsimulated[]QHdownbursts[C],FirstWorldCongressonAssetManagement.2006.立宏胡毅李景禄等2006329-。巧龙.输电线路风偏放电的影响因素巧究阴.高电压技术4:121.,,,,,()ILl.PsrwindcareadnwnuTENCONILihononetaarameterfousedovehtransmissiolinesifalt.[]g,g,gand2006.2006IEEERegion10Conference.IEEE,2006.n巧李孟春,张艳玲,杨北革等悬垂绝缘子串风偏最小间距距离计算分析[J]电测与仪表20-,口,49(555):71013王声学.500kVJ.2008[],吴广宁范建斌,等翰电线路悬垂绝缘子串风偏闪洛的研究[电网技术,,,]329-669():5.-14李瑞祥,张志远导线悬垂绝缘子串风偏角的解析计算饥.电为建设,199415(4);1921[],-15徐钟济.蒙特卡罗方法[M],上海:上海科学技术出版社1985,62357[],。句ShinozukaM.Simulationofmultivariateandmultidimensionalrandomprocess故.TheJournalofTheAcousticalSocietyofAmerica.197172 浙江大学硕击学位论文参考文献。7]口eodatisQSimulationGfergodicmultivariatesU)chasticprocess的,JournalofEngineeringMechanics,SCE28-A19%12778787.,,(),-18张相庭工程结构风荷载理论和抗风计算手册M.上海:同济大学出版社,1990,10:3637[][]王之宏11-19风荷载的模拟巧究J,99454452[],[]建筑结构学报,1,(),20lannuzziASinelliPArtificialwi打denerationandstructuralrespon化JournalofStructural[],p,g,Enneerin-gig,ASCE,1987,113(12),23822398.21LiYKareemAARMAssteminwindenineerinrobabilisticEngineeringMechanics199052[],,ygg.P,,(),50-59.-张相庭结构风压和风振计算:.M?上海:同巧大学出版社98556090[巧[],1,口3ShinozukaM,JanCM,Digitalsimulationofrandomproc的sesandi化alications,JournalofSound]pp-andVibration197225111128.,,y),PWLiYousun,KareemAhsan.Simulationofmultivariaterandomprocesses:hybridDFTand出italfilteringg-aroach.JJournalofnreccEgineeingMhanis19%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}Yan8,LinX,LuoW,巧al.NumericalStudyonDynamicSwingofSuspension虹sulatorString虹OverheadTransmissionLineunderWindLoad口.PowerDelivery,IEEETransactionson,2010,2Sy):]248-2化-—5lm货J..D99.AliltticetII.Aridi[。Ho,"巧ongwndresponseofaowerseodynamcampnganddeflecnnrucures ̄tions口.EgineerigStt,187,483488.]()[52]ASCE:Guidelin的forElectricalTransmissionLineStructuralLoa出ng.Publishedby化eAmericanSocietyofCivilEngineersNo.74.2003.53LoredoSoiiza.ThebehavioroftransmissionlinesunderhihwindD.Ontario:Theuniversitof[]g[]ywesternOntario,1996.54孙珍茂.输电线路舞动分析及防舞技术研究[D].梳州:浙江大学,201化[]74 浙江大学硕壬学位论文参考文献uX--55YANBoLiiaohui,LiuBaoan.Dicresonseanalsisofwindaeawofoverhead[],ynampygynnesnonreachan9transmissioli[J.WSEASTransactiosAppliedandTheoticlMeics,2006,1(1):0.]56GB50009-2012建筑结构荷载规范[S.[]]?57GB-505452010110750kV架空输电线路设计规范[S].[]—5J.2742832.[巧滲宗高?特高压输电线路电力建设2006:,张华,陈海波设计风速取值的探讨[],,()59]EC608262003DesinCriteriaofOverheadTransmissionLines[S].[g-werand—60BS81001:19%Latticetosmastsarti:codeofpracticeforloading[S][]p6eeBE.Theefectofturbulenceon化esurfaceressurefieldofasrerism?J.Fluid.Mech.,1975,[。Lpquap口]'69-:2的282.62]KareemA.Windexcitedmo巧onofbuildingsD.PhD(dissertation,Dept,ofCivilEngineering,Colorado[[]StateUniv.FortCollinsColo1978.,,,63.随桃振动的虚拟激励法[M].北京2004林家浩,张亚辉:科学出版社,[]-64.27223[}林家浩,钟万热关于虛撒激励法与结构隨机响应的注记机计算为学学报9%15苗1,1,()65LinJ.H.AfastCalorithmofPSDmatricesforrandomseismicressesJ.ComPutersand[]QCgpon[]S—turetuers199244:683687.,,9)W—non66LinJ.H.n.S.Willima.Feudoittilrihmfoirrandomseisi[],Zhag,sW.Psexcaonagotrstatonaymcresn-巧6ponses.EgineeringStructures199416:270机,67.隙贤川.大務度屋盖结构风致响应和等效风荷载的理论研究及应用[D]浙江大学博壬学位论文,[]2005.6.大跨度屋盖结构的抗风巧究D.浙江大学博壬学位论文2004.[巧沈国辉[],69章李刚?大型复杂建筑结构风致效应及等效静力风荷载巧究[D?浙江大学博壬学位论文2013.[),]70]鲜荣.双层J.[,廖海黎,李明水马鞍型空间网壳屋盖风振系数非线性模燕叠加法计算[]建筑科学,2008,24(9)1(M8.uma-7avenortA.ClGustloadingfactorsJ.Jolof化eStructralDivisionASCE.196793:1134.[UDp[],,M.Kasperski.Extremewindload出stributionforlinearandnonlineardesign?En.Struct"1992,17巧口]g-14l:2734()73XinzhonChenAhsanKareem.EuivalentStaticWindLoadsonbuildins:newmodelJ,J.of[]g,qg[]Srucurenneerin-ttEi200413010:14251435.gg,,()74XinzhongChen,AhsanKareem.Coupleddmicanalsisandeuivalentstaticwindloadsonbuildins[]ynayqgw-Ji2005-2ihhreedimensionalmodes.J.fSrucureEnneerin1317: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放1RJHolm货J.D.Useofeienva山说in化ecovarianceinterationmehodfordeterminationofwind[巧jggt-urna-loadefecteJ.JolofWindEnineerinand虹dustrialAerodynamics.l98313J59370.[]gg,rrofndecn-uolmesJ.DBestR.J.Anaoach化出edeteminationwiloadef化olowrisebildins.巧引H,ppg口]A-Joul.ofWindEngineeringandIndustrialerodnamics1981,7273287.rnay;76 浙江大学硕壬学位论义作者简历作者简历杨悦,男,汉族,1990年11月出生,浙江金华人2008-.092012.06浙江理工大学建筑工程学院,获得工学学去学位 ̄2012-.09今浙江大学结构工程专业,攻读硕壬学位攻读硕壬期间科巧成果1楼文娟杨悦吕中宾锋杨伦.考虑气动阻尼效应的输电线路风偏动态分析方[],,,张少,法[J].振动与冲击(录用)[2]摆文娟,杨悦,卢明,吕中宾.连续多跨输电线路动态风偏特征及计算模型研究电[化力建设)(录用77

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。
大家都在看
近期热门
关闭