热门搜索 :
考研考公
您的当前位置:首页正文

频激励下干摩擦阻尼叶片的共振响应

来源:伴沃教育
天津大学硕士学位论文

多频激励下干摩擦阻尼叶片的共振响应

姓名:孙艳红申请学位级别:硕士专业:一般力学与力学基础

指导教师:丁千20040101

中文摘要叶片事故的统计分析表明,叶片损坏大多数是由于叶片振动产生的动应力过大所致。为了降低叶片的振动应力,汽轮机厂家越来越多地采用带有阻尼围带或凸肩的叶片结构,即阻尼叶片。阻尼叶片通过围带或凸肩间的相互干摩擦来消耗叶片能量,起到减振的作用。本文对干摩擦阻尼叶片多谐波激振Ⅱ向应进行分析研究,主要工作分为以下几个方面:第一部分:叶片阻尼机理和激振力的介绍分析介绍分析了各种阻尼产生的机理,以及相应的模型,特别是对界面摩擦阻尼的机理及其数学模型进行了较为详细的探讨。对摩擦面间正压力变化时摩擦力进行分析。激振力研究是阻尼叶片动力特性研究的一个重要组成部分。在激振力的各阶谐波分量中,影响叶片安全运行的激振力仅是接近叶片自振频率的分量,本文考虑第一阶、第二阶和第Z阶谐波分量对叶片动力特性的影响。第二部分:叶片单自由度和两自由度模型分析首先,建立叶片单自由度模型,利用平均法求解叶片的运动方程,研究一阶和二阶主共振情况下幅频响应曲线,分析各个主要系统参数对频响曲线的影响规律。然后,考虑弯曲和扭转耦合情况下的叶片模型,建立两自由度的叶片运动方程,应用谐波平衡法得出运动方程的周期解,分析一、二阶主共振情况下系统参数对频响曲线的影响,并与数值解进行比、较验证解析解的准确性。最后对全文进行总结。本文研究得出,增加叶片的界面摩擦阻尼是一种非常有效的叶片减振方法。由于界面干摩擦的介入,使得阻尼叶片具有非线性的动力特性。与线性动力系统不同,阻尼叶片的频率不再是它的固有特性,其频率受外界激振力的影响。本文对阻尼器接触端正压力、汽流激振力、摩擦晃面的剪切刚度、摩擦系数等参数对叶片动力特性的影响进行了理论、计算和设计应用的研究。关键词:汽轮机阻尼叶片多频激励共振弯扭耦合ABSTRACTOneofthecommonstructu删failureofbladesisduetoexcessivestressescausedbylargeamplitudevibrations.Toreducethevibratorystressofblades,thedampedbladeswithintegralwireorshroudarewidelyusedbyturbineworks.Asthevibratoryenergyisdissipatedbythedryfrictionbetweenthedampers.thevibrationofbladesisreduced.Thedynamiccharacteristicsofthedryfrictiondampedbladesinfluencedbythemulti-harmonicexcitationswillbediscussedinthispaper.Partone,introductionandanalysisofthedampedmechanismandtheexcitationforceofturbineblades.Themechanismandcorrespondingmodels,especiallytheinterfacefriction,ofmanykindsofdampersarcintroducedandanalyzed.Theanalyticalfrictionforcemodelwitllvariationofnormalpressurebetweenfrictionsurfaceswillbeconsideredinthefollowinganalysis.Theexcitationforceisoneofthemostimportantelementsforresearchofdynamiccharacteristicsofdampedblades.Amongvarious-orderharmonics,theharmoniccomponent,whichapproachesthenatural-frequencyofblade,canhavethemaininfluenceonthesafetyoftheblades.Thefirst.secondandZ-orderharmonicsareconsideredinthispapertoanalyzethedynamicsofblades.Parttwo,dynamicsoftheone—andtwo-dofmodels.Firstly,theone—dofblademodelisstudiedbytheaveragemethod.Theaverageequationsaredetermined,andtheamplitude—frequencycurvesinthefirstandsecond—orderresonancesarepresented.111einfluencesofthesystemparametersontheamplitude-frequencycurvesarediscussed.Secondly,takingintoaccountthecouplingofflexion.torsion.thetwo.dofblademodeliSstudiedbytheharmonicbalancemethodtoobtaintheresponseequations.Theinfluencesofthesystemparametersontheamplitude—frequencycurvesofthefirstandsecond—orderresonancesarcdiscussed.。Numericalsolutionsareusedtoprovetheaccuracyoftheanalyticalanalyses.111eresultsinthispapershowthatthedryfrictioninbladesplaysaveryimportantrolereducingvibrationamplitudeoftheblades.ThenormalforcebetweentherubbingCaninfluencethedampingeffectofbladej.Throughthestudywefoundthatthereanoptimumvalueforthenormalforce,underwhichthevibrationofbladesearlbeeffectively.Thedampedbladesisanonlineardynamicsystem.SOitsnaturalCanbeinfluencedbyexternalforce.Besides,theinfluenceofshearstiffnessofsurfaces,theplaceofthedamperontheblades,andthefrictioncoefficientoncharacteristicsofbladesarealsodiscussed.words:Turbine,Dampedblade,Multi—harmonicexcitation,Resonance,flexion二torsioncouplinginsurfacesisdampedfrequencyrubbingdynamicKey独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作和取得的研究成果,除了文中特别加以标注和致谢之处外,论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得基鲞盘鲎或其他教育机构的学位或证书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示了谢意。学位论文作者签名:劫搏静_L签字日期:2孵年f月f2日学位论文版权使用授权书本学位论文作者完全了解鑫洼盘茔有关保留、使用学位论文的规定。特授权盘注盘堂可以将学位论文的全部或部分内容编入有关数据库进行检索,并采用影印、缩印或扫描等复制手段保存、汇编以供查阅和借阅。同意学校向国家有关部门或机构送交论文的复印件和磁盘。(保密的学位论文在解密后适用本授权说明)学位论文作者签名:孙彩L导师签名/TJ4l签字日期:2a叶年f月}2日签字日期:妒哗年1月J2日第一章绪论绪论1.1课题背景及意义叶片是透平机械的心脏部件,因承受气流周期性扰动力作用而振动,加之复杂恶劣的工作环境,常常造成口f片失效,如由叶片异常振动或承受共振应力带来的高周或低周疲劳引起的断裂。据统计,叶片事故占汽轮机事故的40%。叶片是汽轮机的心脏,也是事故最多阿关键部件,它的安全可靠直接关系到汽轮机和整个电站的安全。由于电力是国民经济的命脉,而全世界约有75%的电力来自火电站中的锅炉一汽轮机一发电机以及核电站中的核反应堆一汽轮机一发电机,因此提高汽轮机的叶片安全可靠性对于满足不断高涨的电力需求,适应国民经济的发展,有着不容忽视的作用。在石油、化工、轻工等重要经济部门中,采用工业汽轮机作为关键动力设备的现代大型企业日益增多,而由于汽轮机叶片损坏造成工厂全面停产、引起严重后果的实例,在国内外都不少见…。自1884年汽轮机诞生至今已经有一百多年,在这期间,汽轮机经历了巨大的变化。最为突出的是单机功率增长了10万倍以上一一由不到10kW发展到136万kw一~平均每6年就会翻一番。同时,蒸汽参数也由初压1MPa、初温2000C左右,分别提高到最高达34MPa、6000C以上的水平,使热效率提高了3倍。这些成就的取得在很大程度上依赖于叶片的开发和研究。首先,单机功率的飞速增长,要求汽轮机通流能力相应提高。而流通能力又取决于叶片高度所构成的通流环形面积,因此,增加叶片高度是扩大单机功率的必由之路,成为各国汽轮机技术发展的共同方向,末级长叶片也就成为衡量汽轮机技术水平的主要标志。由于叶片高度和蒸汽参数的提高,叶片的工作条件也愈来愈严酷:在进气端的调节级叶片,要承受最高近6000C的高温和喷嘴弧段的巨大冲击力:在排气端,则是巨大的离心力和接近两倍音速的湿蒸汽流的冲刷;而且所有动叶片都承受着多种形式的周期性或随机性激振力作用而处于强迫振动之中。这些不利条件使叶片成为汽轮机中影响安全的最主要的因素。据1970年前后英国、法国、原苏联等国的统计,汽轮机事故约有45~72%是叶片造成的…。美国和联邦德国的统计表明,叶片事故在汽轮机各部件中居首位,美国电力研究所EPRI.(ElectricPowerResearchInstitute)指出,美国电站汽轮机强迫停运率的70%与叶片损坏有关。各国统计还一致反映,叶片事故引起的损失往往占全部损失的一半左右【21。:筮二童鳖堡.中国也不例外,根据研究机构对国产6MW和进口10MW以上电站汽轮机的统计,从1973年到1977年这段时期,每年发生叶片事故的机组不下100台次,电站汽轮机的叶片事故率最高可达11%,其中进口机组的事故率还高于国产机组。对照国际上的最先进水平,这样的事故率是偏高的。从1978年开始,事故台数、事故台次、损坏的总级次以及事故率均有很明显的持续下降趋势。提高发电设备的安全性和可靠性有着巨大的直接经济效益和社会效益。特别在普遍供电紧张的情况下,备用容量小,叶片事故往往引起停机以致大面积停电,社会影响极大。EPRI的调查还表明,美国每年平均发生的40余起重大叶片事故,平均每起的修理费用为20万美元,因每次停机引起的替代能源耗费为250万美元,损失惨重。叶片损坏的原因是多方面的。例如、,从发生的机理分,60~80%的损坏原因是振动…;从责任范围分,可以归纳为设计制造、电站运行、外来原因、老化原因等。其中设计原因一般占80%。目前世界上的汽轮机向大容量化、高速化发展,汽轮机的大容量化使得叶片需要比以往承受更大的汽流力、离心力和工作温度。受汽流力的激振,叶片容易因动应力过大而导致损坏。安全和高效是设计人员的目标,而叶片的安全性则被考虑放在首位。叶片设计者采用各种手段来降低叶片动应力,其主要措施是增强叶片的阻尼环节。.对产生阻尼的机理进行分析表明,叶片凸肩及平台处的干摩擦阻尼的减振效果最为突出。由于阻尼叶片良好的减振性,阻尼结构越来越多地应用于燃气及蒸汽轮机之中。阻尼叶片可以很有效地克服叶片动应力过大的情况,并可以起到调频作用。目前使用最多的阻尼结构是预扭围带结构和叶片凸肩结构。国内外各大汽轮机公司出现的具体结构比较多。如美国西屋公司的Z—lock型,文献¨’阐述了此种阻尼结构对扭转恢复的抑制作用,并用计算软件NASTRAN对带此种阻尼件的叶片迸行了计算。计算结果表明“Z-Lock”型阻尼减振器对叶片的扭转恢复有较好的抑制作用和良好的成组特性,并影响叶片的振动频率。此机构的特点是阻尼减振器的接触面基本是与轴向垂直的面。GEC/ALSTHOM公司的翅翼型与z~lock型有所差别,此种阻尼连接件的相互接触面与轴向成一角度,接触面也较大,可以提供更为稳定的摩擦力n理论和实验均表明依靠接触面间相互滑动产生摩擦来增加阻尼是一个极为有效的方法【41。GE公司采用了一种套筒凸台连接形式的阻尼结构,这种结构形式通蔓二童缝堡。过叶片上的凸台和套筒连成整圈,安装时使叶片先预扭,正常运行时,就处于松连接的状态,用这种连接形式可以确保动叶在振动中有足够的阻尼,以利于改善动态特性"3ALSTHOM公司采用鳍片式阻尼结构形式,叶片在连接成组时凸面嵌入凹面内,连接件之间有两个交叉的斜接触面。振动时,阻尼减振器的接触面将产生切向与轴向的位移耦合,无论哪一个面接触,切向滑动都会导致轴向滑动,轴向和切向相对位移满足几何比例关系,径向的相对滑动相对独立,鳍片式阻尼结构对扭转恢复的抑制及减振效果都很好[61。目立公司在其新设计的26英寸末级叶片上采用了s型复环成圈阻尼结构,这种结构是在末级叶片上整体铣出S型复环,叶片安装后,阻尼器的凸凹面相互嵌入,增加叶片间接触面积,并且采用整圈连接复环结构,提高了叶片的安全可靠性17】。三菱重工在七十年代其开发的末级叶片中使用了Z型阻尼连接技术,该结构与GEC/ALSTHOM公司的翅翼型结构相似,但是三菱公司的阻尼结构翅翼较瘦小,相互的接触面也较小。【6]介绍了此种结构,其形式是数片叶片焊接成组,然后叶片组之间通过有间隙的阻尼结构相连接。阻尼端在叶片静止时是相互不接触,运行时由于叶片扭转恢复的作用而相互接触,从而通过接触端之间的摩擦作用达到减振的目的,三菱公司做的实验表明,采用这种阻尼结构的减振器可以收到极为明显的减振效果。此外,三菱公司还用到另外一种整体围带式阻尼连接结构,此种结构在18j中有介绍,该结构是使叶片与叶片之间全部都用阻尼结构相连接,当叶片旋转时,在扭转恢复的作用下整圈的叶片产生相互的耦合。与传统的成组叶片相比,它提供了更多的阻尼环节,因而动应力大幅度下降,结构的可靠性得以显著提高。在一个实际尺寸的试验透平上迸行了实际负荷的实验,结果表明这种阻尼减振结构叶片的对数衰减率比一般成组叶片高5—10倍,即使是在共振状态下其动应力水平也很低19】。国内,哈尔滨汽轮机厂设计了自带围带预扭叶片,预扭原理与ALSTHOM公司的预扭叶片原理相似。预扭围带采用平行四边形自带围带,因其具有良好的成组及阻尼作用,可使共振危险性大为减少,与铆接围带相比可使围带应力减少到原来的1/31”}。研究表明,此种围带有以下几个优点:1.能满足两种叶根(T型叶根及枞树型叶根形式)两个装配方向(切向和轴向)的要求;2.轴向连接刚性大,且随着节径数增加,频率升值快,对避开低节整二童堑鲨。——一径三重点共振有利(ALSTHOM公司元宝山600MW机组资料证明);3.围带工作接触面大,对一定的挤紧力,摩损小,表面不必硬化;4.允许围带接触面存在局部间隙,便于拂配及预扭安装。这些都是叶片阻尼结构,当叶片旋转起来,离心力使叶片发生扭转恢复而使得相邻阻尼接触端相互贴合,形成摩擦面,起到减振的作用。近年来,在阻尼叶片方面的研究有很多进展,但大多只是针对某些方面,没有系统化,而且定量研究很少,计算也限于简单结构。文献…1用NASA.C.R软件计算过单叶片和类似叶片的结构,结果与实验研究的结果相符,褥出关于阻尼件的几个关键参数:接触压力,接触表面粗糙度,外力水平,阻尼器位置,阻尼器刚度。1.2国内外阻尼叶片的研究现状阻尼叶片(在叶片顶部加工出整体围带、在叶型上加工出凸肩或装加松拉金)越来越多地被应用于燃气或蒸汽轮机。理论和实验研究证明。阻尼叶片可以有效地降低动应力及对叶片起调频作用。因此,深入研究阻尼叶片,探讨其各个关键参数对叶片动力特性的影响作用,以指导实际设计工作势在必行。国内外学者相继开始了对阻尼叶片的研究探索。1.2.1界面摩擦模型在摩擦模型的研究上,以Oden…11为代表认为,两摩擦面间的摩擦过程是先经历一个弹性过渡区,然后再发生滑动,即摩擦面间的摩擦过程可以分成静止、弹性变形和滑动三个阶段,而不是一般库仑摩擦行为所定义的静止和滑动两状态。A,V.SrinivasanIl3】发展了Oden模型,将滑动速度对摩擦力的影响也考虑进去。C.H.Menq[141在前人的基础上提出了微滑移模型,微滑移与宏滑移不同的是,它考虑了摩擦面各点处的摩擦行为的不同,比如在一定的受力情况下,可能是某些区域发生滑动,而另一部分静止。微滑移模型可以较好地模拟这样地情况。C.H.Idenq【151提出了二维滑动的模型,对摩擦滑移的描述上又更进~步。1.2.2叶片建模从80年左右,以J.H.Griffin【I6]为代表的一些学者开始了对阻尼叶片动力特性的研究。他发展了摩擦力一位移迟滞模型,提出了弹簧、质量和阻尼器的阻尼叶片描述系统,得出了一些阻尼叶片的动力特性,讨论了阻尼端正压力对系统频率响应的影响,并不断提出新的研究成果。篁二童缝鲨一一——一VasutomoKANEKO(1993)用多自由度质量、弹簧阻尼系统研究了阻尼叶片。国外学者大都采用单自由度或多自由度弹簧一质量阻尼系统。虽然这些解的物理概念比较清楚,但对解决工程实际问题略显欠缺。实际叶片几何形状复杂,目前较有效的工具方法是有限单元法。Srinivasan【l71曾用有限元法计算过阻尼叶片,并与实验值进行了比较,得出了有益结论,但是结论的轮廓欠清晰,也没有说明具体的计算方法。国内任勇生118】研究阻尼叶片较早,提出了有限元方法研究思路,并认为阻尼叶片的非线性特性是一个局部非线性动力问题。1.2.3阻尼叶片的计算分析在计算分析上,NAYFEH[”】提出了谐波平衡法来分析非线性动力问题。谐波平衡法是比较理想的求解非线性动力方程的近似解法,在满足工程精度的前提下,计算速度快。为了提高计算精度Pierre【20】提出了增量多项谐波平衡法,Wang[21】提出了多项谐波平衡法。Cameron[22J提出了用时域一频域交替算法,利用了在时域和频域解问题的各自优势,并给出了相应的算例(限于单自由度)。1.2.4阻尼叶片的设计‘在阻尼叶片设计方面,以Griffin为代表的学者较系统地研究了阻尼叶片系统地一些非线性动力特性,得出了一些有益地结论。在此基础上,Kielbl23】提出了阻尼叶片的一些设计思路,向工程实际迈进了一步。国内外对叶片阻尼器的研究有了一些结论,但是都处于比较抽象的理论探讨之中,真正实用性的研究还是很少,实验研究就更少一些。阻尼叶片已经在实际机组上运行,如何设计阻尼叶片,使其最有效地降低叶片地动应力,以及由于阻尼器的应用,它会给叶片的频率及阻尼特性带来什么样的影响,将是迫切需要解决的问题。要全面研究阻尼叶片,叶片汽流激振力的研究也是~个必不可少的研究环节。阻尼叶片的非线性动力特性与外部激振力有密切的关系,而叶片激振力的研究是一个很棘手的研究课题。一个新的实际工程问题的出现总是需要理论和实验研究的进步,也是促进理论和实验研究的进步,对于阻尼叶片的研究也是如此。因为理论研究在建立模型时,总是将一些因素简化或是忽略,这必然给理论分析带来误差,而且在对问题还不是很清楚之前,建立的模型可能忽略了~些重要因素,而使得理论分析失败。实验研究可以避免理论分析的缺陷,将所有实际因素都再现在实验结果之中,并且实验研蔓二茎究结果可以验证和指导理论分析的进展。阻尼叶片几何形状复杂,在汽流力的激振下,其振动特性非常复杂,非线性振动理论和解法也不成熟,有很多的问题需要解决。因此,阻尼叶片的动力特性研究是一个非常崭新的挑战性前沿课题。为了减少流动损失和增加叶片刚性,在涡轮机中,往往采用围带和拉筋将几个叶片连接成叶片组。几个叶片与围带和拉筋组成的振动系统塑丝——一一具有多样的振型和复杂的频谱。对于等截面叶片组振动计算,普遍地采用系数图线法。对于变截面叶片组振动计算,早在五十年代,Prohl就曾经提出过围带叶片组的振动计算方法【24。,它忽略了切向与轴向弯曲振动的耦合,只将轴向弯曲振动与扭转振动耦合在一起,这就局限它只适用于截面的最小主形心轴与轴向重合的叶片。同时它还忽略了剪切变形、转动惯量和离心力的影响,因此无论将它用于短叶片或长叶片,计算值都不够精确。Deak等提供了曾在美国Allis—Chalmers公司多年应用的汽轮机末级叶片组的振动计算方法[251,它将切向振动、轴向振动和扭转三者耦合在一起,并考虑了边界叶片和中间叶片的差异以及离心力的影响。它既可用于围带叶片组,又可用于几排拉筋连接成的叶片组。但是没有考虑剪切变形和转动惯量的影响,用于高压级叶片受到限制。Thomas等应用有限元法计算叶片组的自由振动频率【2“,用Sturm序列和反迭代法求矩阵的特征值和特征向量。可惜它只是考虑一个平面上的弯曲振动,也忽略了剪切变形和转动惯量的影响,几排拉筋连接成的叶片组的情况也没有考虑。林俊灿12。7】采用传递矩阵法,它将叶片作为变截面梁的非对称弯曲,将扭转振动与两个平面内的弯曲振动通过拉筋耦合在一起。并作了以下改进:考虑了叶片和拉筋的剪切变形和转动惯量的影响;合理地考虑了离心力的影响;在动态下计算叶片与拉筋之间的作用力和力矩,这样既提高了它的计算精度,又扩大了计算机程序的通用性。1.3本文的工作内容和安排摩擦阻尼叶片动力学的理论研究模型有多种形式,最常用的是平台阻尼器模型,且多是研究单摸态动力学模型,即单自由度质量一线性弹簧一阻尼系统在其共振频率附近的动力学行为。Wang[2副发现,共振频率漂移和最大振幅与速度相关摩擦力的系数有很大关系。Ferril2引用奇异摄动理论研究简化的0质量(无限刚度)和小质量系统的周期解,并将在一定参数范围内的结果与原系统的数值结果进行比较。BerthillierI30】茎二童笪鲨一一用多谐波频域分析结合C-B的局部模态综合方法,研究发现高阶谐波振动成分对叶片响应有很大影响。他们‘31l还用Newmark方法研究了三维结构的弯扭耦合振动。其中,将干摩擦连接件看成几个平行的离散单元。Wang【321用FET研究了分别作用单、多谐波激励的0质量系统黏滑振动响应。在一定激励参数范围内,干摩擦连接件的减振效果颇为明显。季保华[331实验研究了摩擦面压力、接触状态和激振力对叶片动态特性的影响,结果表明摩擦对减振是有效的。近期还有多篇文献研究轴系扭振与叶片振动的耦合问题‘34、351,但连接处简化为固连,或者用摆簧弹性联结,或粘性阻尼连接等等,还未涉及到干摩擦粘滑振动对叶片振动的影响。本文的工作内容和安排如下:l阻尼叶片机理和数学模型的建立叶片阻尼是叶片动力特性的重要参数,主要分为四种类型。本文介绍了各种阻尼,特别是界面摩擦阻尼产生的机理及相应的数学模型,为阻尼叶片运动方程的建立奠定了基础。本文还逐一分析了叶片中多处的干摩擦阻尼环节,比如叶根、叶片平台、围带(凸肩)、拉筋等地方。2叶片激振力的研究对激振力的研究是阻尼叶片动力特性研究的~个重要组成部分。阻尼叶片的非线性动力特性不仅与叶片自身的材料、约束条件有关,还与激振力有关。在激振力的各阶谐波分量中,影响叶片安全运行的激振力仅是接近叶片自振频率的分量,其它分量虽使叶片做强迫振动,但是动力响应与前者相比很小。因此,降低接近叶片自振频率的激振力谐波分量是主要研究的问题。3叶片单自由度方程的求解建立单自由度的千摩擦阻尼叶片振动模型,应用平均法分析了一阶和二阶主共振的情况,通过幅频响应曲线得出主要的系统参数对频响曲线的影响规律,得出对实际应用具有指导意.义的结论。4考虑叶片弯扭耦合的模型在阻尼叶片的模型中采用两个并联的阻尼器,并且考虑弯曲和扭转的耦合,建立了两个自由度的运动方程。应用谐波平衡法进行求解,分析了频响曲线的形状,得出主要系统参数对频响曲线形状的影响,并通过与数值解的比较验证结论的准确性。5对全文进行工作总结,指出可能存在的问题和今后研究的方向。第二章叶片阻尼机理及激振力分析第二章2,1概述叶片阻尼机理及激振力分析阻尼是叶片动力系统的一个关键参数,它对估算叶片的共振动应力起到非常重要的作用。目前的研究表明,叶片的阻尼环节大致可以分为以下四种:材料阻尼:由叶片材料及质量分布所决定的固有的阻尼特性;摩擦阻尼:叶片在接触面(叶根、叶根平台、围带或凸肩、拉筋等部位)的摩擦而产生的对叶片能量的消耗;气动阻尼:流体对动叶作用而产生的阻尼,包括粘性阻尼及气动弹性阻尼;冲击阻尼:叶片与叶片或叶片连接件之间在叶片振动过程中的相互碰撞所造成的阻尼。上面四种叶片阻尼中,叶片材料阻尼和摩擦阻尼是主要的,对叶片阻尼的贡献最大。叶片粘性阻尼约占叶片材料阻尼的10%,气动弹性阻尼所占的比例很小,但是气动弹性阻尼有时候会引起叶片的颤振。干摩擦阻尼是叶片阻尼中最有效的阻尼方式。要准确地估计干摩擦对阻尼叶片的减振贡献,必须建立起对它的合理数学描述。干摩擦阻尼作为一种摩擦物理现象,它不仅表现摩擦面之间的力的关系,还表现在对能量的消耗上。文献[12-14]对摩擦模型的建立做了有益的工作,但是他们的工作都是建立在单自由度系统的基础上,不容易推广使用,因此还有很多的工作要做。。2.2材料阻尼及其数学表达很多研究发现,叶片的材料阻尼是叶片的固有特性,它分布在整个叶片结构之中”“。叶片的材料应具有高阻尼特性,保证叶片有足够的抗振强度,避免损坏事件的发生。研究材料的力学特性时,往往把材料看成是纯弹性性质的,即完全遵循胡克定律,特别在低应力情况下是这样的。但事实上,即使在低压力条件下,材料的应力和应变也不完全遵循胡克定律。宏观上连续的金属材料会在微观上因应力或交变应力的作用而产生分子或晶界之间的位错运动、塑性滑移或其它原因耗损能量产生阻尼。在低应力状况下由金第二章叶片阻尼机理及激振力分析属的微观运动产生的阻尼耗能,称为金属的滞弹性,可以由图2.1加以说明。当金属材料在周期性的应力作用下,加载线OPA因上述原因形成略有上凸的曲线而不再是直线,而卸载线AB将低于加载线OPA。于是在一次周期的应力循环中,构成了应力一应变的封闭回线ABCDA。对于阻尼等于零的完全弹性材料,封闭回线将退化为面积为零的直线OAOCO。产生阻尼耗能的力学现象是应变滞后于应力,才构成这一回线,所以称之为滞迟回线,阻尼消耗能量的值正比例于封闭回线的面积,这一现象是由金属的滞弹性所造成的。严格地说,金属地滞弹性与频率和温度是有关的,在特定的温度和频率下,金属内部分子易于产生位错现象面耗能;这些耗能现象当然都是指机械能转换成热能,然后向周围环境耗损。在特定的温度及频率下,金属的微观结构也会将振动直接转变为热流而耗损能量,所以,金属的这种滞弹性虽然和应力幅值及疲劳周数无关,却和频率与温度有关。l盯、,,研c么。佰i/图2.1应力应变迟滞回线当金属材料承受的应力逐步上升,达到疲劳极限80%左右的高应力状况时,晶体内部应力的不均匀性及晶界之间的应力集中,会使某些部分的应力达到足以引起内部塑性应变的水平。即使宏观的或平均的应力很低时,这种局部的塑性应变也已经发生,而在达到高应力状况以至接近疲劳极限时,这种局部塑性应变的体积占总体积的比重便急剧增加,由此而使机械能的耗损大幅度增加,成为产生阻尼的又一个原因,或又一种机理。所以,金属在低应力状况下,是其滞弹性产生阻尼,而在应力增大时,局部塑性应变渐渐变得重要起来,以至高应力状况时成为产生阻尼的原因,其间没有明显的分界。而且,因为这种机理在应力增长过程中都在起作用而且发生变化,所以,金属材料的阻尼在应力变化过程中不为常值,而在离应力或大振幅时呈现出较大的阻尼。9第二章叶片阻尼机理及激振力分析B.J.Lanza等人对叶片的材料阻尼进行了大量的研究,发现材料阻尼与金属的成分、热处理方式、动静应力和温度有关。研究发现叶片材料阻尼的对数衰减率与叶片的振动幅值有定量关系,材料阻尼公式为:d=Jcr”(2—1)上式中d是单位体积的材料耗能,盯是局部应力,H是指数,L,是比例常数,整个叶片系统的材料耗能为:D=Iddv=肛”dv(2—2)叶片系统的应变能为:u=i1肛。(2—3)那么,叶片材料的对数衰减率为:一:一D5=——2U(一)(2—4)、7对于AISl403不锈钢的均匀粱,Lazan取":2.9,J:2.0×10—13。E=28.4x106lb/m2,如果根部应力为5000psi,计算万=0.014,这种阻尼值的估计方法一般与试验值吻合较好。文献‘”1给出了叶片材料阻尼的计算方法。2.3气动阻尼气动阻尼是动叶与流体(蒸汽和燃气)相互作用而产生的对动叶的阻尼。气动阻尼进一步可以分为粘性阻尼和气弹阻尼,它们的产生机理有所不同。2.3.1粘・眭阻尼‘叶片周围的工作介质对叶片产生的流体粘滞阻尼,称为粘性阻尼[38J。由于流体的粘性,在叶片壁面上产生固液界面的摩擦,并且叶片在振动同时,要带动周围介质一起运动,即一部分能量将传给气体介质,叶片的振动被衰减。2.3.2气弹阻尼除了粘性阻尼贡献外,介质对叶片阻尼还有一部分阻尼贡献,即气O第二章叶片阻尼机理及激振力分析弹阻尼。气弹阻尼是指叶片振动与气体流动相互作用所产生气动弹性作用对叶片振动的能量衰减。由于叶片的振动,引起叶片周围介质流动参数的变化,从而引起叶栅升力(推动力)的变化,表现为一种阻尼形式,称为气弹阻尼。一般情况下,气弹是阻止叶片振动的,但是在某些情况下,它会促使叶片振动而产生叶片的自激振动。在汽轮机或燃气轮机叶片中很少会出现自激振动现象13…。实际运行经验和实验研究证明气动阻尼在叶片阻尼中占的比例比较小。这部分阻尼可以用实验的方法测得,方法是分别在介质和真空中测得叶片的阻尼值,然后比较测试结果,便能得出叶片各阶振动的阻尼值[39J。文献[401通过实验测得了一个模型叶片组的气动阻尼比,测试结果表明在一定的叶片振动幅值情况下,叶片的气动阻尼与流体马赫数有关。在测量范围内叶片的气动阻尼比有正有负,与马赫数有关。文献【41】介绍了跨音速压缩机叶轮叶片阻尼的测试方法。2.4冲击阻尼在实际运行中,叶片振动时还存在着一种撞击阻尼,通常称之为冲击阻尼。这种形式阻尼在叶根处、拉筋或叶片凸肩与围带之间。当叶片振动幅值比较大时,如机组升速或降速过l|缶界转速时,叶片与拉筋之间或叶片连接件之间有脱离发生,发生冲击碰撞,这种碰撞是一种非完全弹性碰撞从而消耗能量。碰撞时发出的声音也是一种能量形式被释放出来。目前对这种阻尼的研究还比较少,定量研究就更少一些,需要深入的研究。但一般认为,这种阻尼在机组正常运行时对叶片阻尼贡1献较少。2.5界面摩擦阻尼两个相互贴紧的接触面承受一个附加的激振力t,当疋逐渐增大时,假设阻尼结构不发生变形,但是在结合面中仍将产生相对变形△或结合面内的接触应力和应变。图2.2表示最初的相对位移和变形与力呈线性关系,是弹性变形OP;然后,进入非线性变形区PB,随着F口的进一步增大,达到了结合面的正向压力所允许的某一极限值,在结合面中产生了相对的宏观滑移,这时两结合面之间在接触部分发生了干摩擦或库仑摩擦。到达终极位置C以后外力E反向,力一位移线沿OP平行方向下降,至p’点后有一段非线性卸载线P、B、,到达反向极限以后,再次发生相对宏观滑移13’C‘,力的循环形成力~位移的一次周期及回线C‘BCB’C‘。对第二章叶片阻尼机理及激振力分析应于图2.2中,如果在B点以前名即开始反向,便形成如图2.3所示的力一位移周期回线ADA’D、A。交变力e与咬变弹性位移△在结合面中所损耗的能量正比例于回线所包围的面积,产生结合面滞弹接触阻尼。相反,由宏观滑移BC及B、C、与交变力所损耗的能量,正比例于连接BCB、C’所构成的菱形面积,产生库仑摩擦阻尼。|『^clfnfn砖//lO・/.一c:。荔.#一一67晌。O/DAA,厶l/图2.2摩擦力与相对位移迟滞回线图2.3摩擦力与相对弹塑性变形迟滞回线当结合面之间正向应力增加时,首先相互接触的峰顶部分产生弹性变形,继续增加正向压力,峰顶部分到达屈服强度而产生塑性变形,但是周围区域仍然处于弹性变形阶段,因此,接触区便形成弹性变形和塑性变形共存的状态。任何一个接触面的表面都是高低不平的,有很多的接触点,这些很多的接触点表现出来的力学特性自然就是初期为线性,然后呈非线性。在正压力作用下,结合面的各不平度峰顶接触部分相互嵌入,在承受切向交变力t时,结合面之间就产生了交变的和微观的变形与位移;显然,这一变形或位移应该呈现出非线性性质,并形成可以耗能的滞迟回线,这~种变形亦可以称之为周期性的滞迟变形。其本质就是接触区的局部塑性应变产生的能量耗损。从表面现象上看,库仑摩擦阻尼和结合面接触阻尼有相似之处,它们都来源于结合面之间的相对运动,但是从阻尼的产生机理上比较,两者之间存在差别。如果说结合面接触阻尼是由微观的变形产生的,那么,库仑摩擦阻尼则由结合面之间相对宏观运动的干摩擦耗能产生的,它所耗能量可以通过分析计算,由摩擦力一位移滞迟回线所包围的面积得到。图2.2滞迟回线所包围的面积,实际包含了结合面接触阻尼和库仑摩擦阻尼两部分耗能。在各类库仑摩擦阻尼器中,都设计了很大面积的摩擦塑三童堕苎堕星垫墨垦塑塑杰坌近一耗能面,相对位移的大小远远大于结合面之间的微观变形。在叶片系统中有结合面摩擦的主要部位有:叶根处,叶根平台,拉筋,叶片围带或凸肩处(专门作成的阻尼结构)。下面分别对其机理进行分析。关于界面摩擦阻尼的微观和宏观机理可以参考文献¨2171比较详细的论述。2.5.1叶片摩擦阻尼的主要分类l叶根处的摩擦阻尼当叶片装入叶轮轮槽之中,叶根与轮缘的接触部分就提供了相互滑动的摩擦面。研究表明,叶根处的摩擦在机组转速比较低的时候对叶片的阻尼贡献比较大,而在机组转速比较高的时候,叶根处的摩擦对叶片的阻尼贡献就比较小了。当机组转速比较低的时候,叶片的离心力比较小,叶根与轮缘之间的接触压力也就比较小。叶片在振动时,叶根与轮缘之间比较容易产生相互摩擦滑动,对叶片振动起到比较大的阻尼作用。而当机组转速逐渐上升时,叶片的离心力也逐渐变大,叶根与轮缘之间的接触力也变大,叶根与轮缘之间的滑动就比较困难。这时,叶根处的摩擦就变的比较薄弱。J.S.RAOH8J实验研究了自由叶片的阻尼比与汽轮机转速之间的定量关系。实验时将叶片转子置于真空中,通过用电磁激振器来激励叶片,测量叶片根部阻尼并给出了实验曲线,。如图2.4所示。从图上可以看出,各阶模态曲线的趋势是随着机组转速的上升,叶片的阻尼比里减小的趋势,特别在转速为400RPM左右,叶片的阻尼比有一个突降,文中称其为阀值。这个阀值就是叶根与轮缘之间的锁定值,即由于较大的离心力,叶根与轮缘之间的滑动摩擦被“制动”的速度。剩余的阻尼主要是叶片的材料阻尼,而在阀值之前,叶片的阻尼主要是叶根与轮缘之间的摩擦阻尼。・2叶片平台处的摩擦阻尼叶片平台阻尼装置常见的形式有以下三种。第一种是用专门设计的阻尼结构,使阻尼结构与叶片安装平台下部发生摩擦,这种阻尼结构是靠离心力作用使阻尼器与叶根平台下部发生紧密接触的[391。第二种形式是叶根平台与叶轮某固定处发生摩擦,许多文献称之为叶片一基础摩擦(BLADE—GROUND)[t3J。第三种形式就是叶根平台与叶根平台之间的第二章叶片阻尼机理及激振力分析摩擦(平台之间装有摩擦垫片),特别是同向装配的叶片。对于叶根平台的摩擦运动,叶根平台与相应摩擦面之间的正压力起着非常重要的作用,这影响到阻尼结构的耗能值。理论及实验表明叶片平台与相应摩擦面之间存在着最优正压力,在最优正压力下阻尼器对叶片的耗能最大。Y~~~厂:■可1gl口善e暑&一定.二=,.、.Zi枷R露面(删!图2.4叶根阻尼示意图图2,5叶片围带(凸肩)示意图3阻尼拉筋与叶片拉筋孔之间的摩擦阻尼在机组起动后,由于离心力作用,拉筋与叶片拉筋孔之间相互贴合面之间有一接触压力。当转速较小时,接触压力也较小,拉筋与叶片之间可以相互摩擦滑动。当转速较高时,由于较大的离心力而使拉筋与叶片之间压紧力变得很大,拉筋孔与叶片之间的摩擦滑动变得困难,摩擦对阻尼的贡献便变得很小,这时拉筋对叶片(叶片组)的频率及模态影响比较大。因此,很多的研究学者认为阻尼拉筋在机组起动或停机时对叶片的减振效果比较好,可以使叶片避免过大的瞬态动应力。4叶片凸肩、围带处的摩擦阻尼前面所述三种结构摩擦环节是前人研究较早的结构,它们之间有一个很相似的特点,就是摩擦面之间正压力与转速有关,并且这个正压力不容易实现调节。研究表明摩擦阻尼对应着一个最优的正压力,在最优正压力的情况下,阻尼效果最好。因此,要获得最优的阻尼效果,就要实现阻尼摩擦面正压力的方便调节。叶片凸肩、围带阻尼结构就可以实现这种调节。如图2.5所示,在机组静止时,叶片围带(凸肩)之间有一间隙,而当机组旋转耐,由于叶片的扭转恢复作用,相邻叶片围带或4第二章叶片阻尼机理及激振力分析凸肩相互接触贴合,当叶片振动时,围带或凸肩之间便发生相互摩擦滑动。理论及实验研究表明,叶片凸肩或围带之间存在最优的正压力,因此,如果想要调整阻尼接触端的正压力值,使叶片的振动应力减至最低,只需方便地调整叶片凸肩或围带之间地间隙,就可以达到。2.5.2摩擦面力学行为的描述在切向力的作用下,如果两个接触面之间存在正压力,那么接触面上就会产生摩擦力。在切向力小于两接触面间最大静摩擦力时,摩擦面之间的静摩擦力与切向力相等。此时,两接触面间虽然没有宏观滑移,但存在弹性变形。如果切向力超过了两接触面之间的最大静摩擦力,则两接触面会发生宏观的相对滑移。前面的分析中,没有考虑摩擦面上局部细微的变化。目前主要存在着两种摩擦模型,一种是整体滑移模型,这种模型不考虑摩擦面在摩擦过程中的局部变化,而是将整个摩擦面同一对待;另一种是局部滑移模型,这种模型认为由于摩擦面间的正压力分布不均匀或者摩擦面各处的材料力学性质分布不均匀,造成摩擦面在相互摩擦时各处的变形和状态不一致。1整体滑移模型整体滑移模型认为摩擦面相互摩擦时,摩擦面上的应力、应变场分布是均匀的,在外力的作用下,整个接触面上的各点的变形是相同的,如图2.6(a)所示。因此,在两接触面发生宏观相对位移的临界点,接触面上的各点或者同时保持弹性变形状态,或者同时发生宏观的相对滑移;这样,可以将整个接触摩擦面作为一个点来对待,用这个点来代表整个摩擦面。如图2.6(b)中A点所示。整个滑移状态有两个,一是不滑动,二是滑动,摩擦面之间的摩擦力可以表示为‘a’‘b)图2.6整体滑移模型筇二章叶片阻尼机理及激振力分析。f:渺发生宏观滑移时IK^Au没有宏观滑移时(2—5)上式中“是摩擦面间的摩擦系数,_v是摩擦面间的法向压力,即正压力,K。是摩擦面间的剪切弹性系数,△“是摩擦面间的相对位移。2局部滑移模型实验证明整体滑移模型在接触面正压力比较小的时候,是比较符合实际情况的。但是对于接触面正压力比较大且分布不均匀的时候,问题变得比较复杂,整体滑移模型就不能准确反映实际情况,需要发展一种新的模型,即局部滑移模型。与整体滑移模型不同的是,局部滑移模型将摩擦面间的滑移分成三个不同的阶段:(1)弹性变形阶段在切向力比较小的情况下,整个接触面还没有发生宏观滑移,但是接触面之间已经有弹性变形了,如图2.7‘(a)所示。图中q是压力分布密度。摩擦面间的弹性变形是当地材料和应力的函数,用占,表示,两摩擦面间的摩擦力是摩擦面各处的切向弹性应力在整个面S上的积分:,=体(,)・日(r)tds(2—6)S上式中的墨(r)是摩擦面间的局部剪切弹性系数。q(a)‘h)图2.7局部滑移模型的弹性变形阶段(2)部分滑移阶段随着切向力增大,摩擦面之间的变形便由弹性变形阶段过渡到了部分滑移阶段,即虽然接触面间没有发生整体的滑移,但是有一部分接触面之间已发生了相互滑移,如图2.8(O-)所示。在弹性变形部分,摩擦力是切向弹性应力在S上的积分:在滑移部分,摩擦力是摩擦切向应力在区域足上的积分。整个摩擦面上的摩擦力是弹性变形部分和滑移部分摩擦力的和,如图2.8(b)所示。摩擦力可以用下式表示:第二章叶片阻尼机理及激振力分析,=忙(,)・。(r)・出+∥IO"I(r)・ds(2—7)马岛式中q(r)是局部法向压应力。芒生④m讯{胁k—l!JH-二(。。I%)\尸——}_—,・:Fth'图2.8局部滑移模型的部分滑移阶段图2.9局部滑移模型的全部滑移阶段(3)全部滑移阶段当切向力增大到一定程度后,在部分滑移阶段的弹性变形区域也发生了相对的滑移,此时整个接触面开始全部的滑移,如图2.9所示。摩擦面间的摩擦力为摩擦面间的切向应力在整个摩擦面上的积分:,=∥J乃p)-ds(2—8)S3干犀擦阻尼力的非线性性质根据整体滑移模型和局部滑移模型对接触面间的摩擦力的表达式,可以方便的推导出摩擦力在摩擦过程中的变化过程,如图2.10所示。图中oa段表示接触面还没有发生整体滑移,只是弹性变形或者是部分滑移,摩擦面问的作用力还没有达到滑动摩擦力州。当曲线超过了a点时,摩’擦面之间便会发生宏观的滑移,摩擦力为f=∥(2--9)摩擦力的变化过程是非线性的变化过程,.在oa段上摩擦力是摩擦面间相对位移的函数,显示了干摩擦阻尼力的非线性的性质m,。.rF0Au图2.10摩擦力的非线-眭特性第二章叶片阻尼机理及激振力分析4滑移面间隙的描述1在叶片振动的时候,两个相互摩擦的接触面有时可能也会发生分离的现象。叶片振动时振幅过大而造成阻尼器的摩擦面相互分离,如图2.1l所示。阻尼叶片在安装时,叶片阻尼器间隙值较大也会造成叶片振动时摩擦面分离(例如某机组的整体围带阻尼叶片口们)。摩擦面间隙的存在将会对阻尼叶片的频率和响应产生较大的影响¨…,进而影响到叶片的内部动应力,阻尼叶片的间隙可以用图2.12的质量一弹簧系统表示。摩擦面间的摩擦力表示为:图2.11阻尼叶片振动过程中的分离现象图2.12阻尼器间隙模型,=k,‰棚妞<Ds(2一lo)‘4Au≥Ds2.5.3摩擦力的数学表述前面介绍了摩擦面间的滑移力学模型,摩擦面问的滑移有两种,即整体滑移模型和局部滑移模型,这两种模型各有优缺点和其使用场合。整体滑移模型适用于两摩擦面间正压力比较小,且压力在整个区域内分布比较均匀的情况。整体滑移模型比较简单,计算量小,对于非线性动力问题的迭代求解来说,比较适用。而且。由于整体滑移模型将整个滑移面作为一个点来处理,月梁单元、扭曲梁单元作有限元计算更是显出其优越性。局部滑移模型适用于接触面正压力较大,分布不均匀的情况,由于它能考虑摩擦接触面的局部细节,因此模型更精确。本文在全面分析阻尼叶片摩擦面所受正压力的情况及权衡了计算量和精确度的基础上,选择了整体滑移模型。1摩擦面正压力不变时的摩擦力表达式接触面间的摩擦力的产生原因是摩擦接触面受到了与接触面平行的切向力的作用,如果此切向力是一个周期作用力,那么两接触面便会发第二章叶片阻尼机理及激振力分析生往复的周期摩擦运动。由于本文研究的是振动问题,因此就以系统受到周期作用力作为研究对象。不失一般性,设系统受到简谐激励。在周期性作用力下,两接触面间的力与位移的关系是一条迟滞回线,如图2.13(a)所示。图2.13(a)与图1不同,在图2.13(a)中,加载线和卸载线是对称的。为了方便计算,可以将图2.13(a)化为图2.13(b)。图2.13(b)中,A点表示口=0,B点表示0=口+,C点表示口=g,D点表示口=0‘+疗。设两接触面间的相对运动为简谐振动,lf/刃一_£——_/i/‘a)‘b)图2.13正压力不变时的摩擦力与位移迟滞回线即表示为Au=Bcos8(2—11)其中口=耐一伊,式中曰为摩擦面间稳态相对振动幅值,埘为激励频率,p为相位差。则摩擦力的表达式为:Ka(6u—B)+州0≤口<口’f=一心日‘≤口<石一Kd(Au一占)一aN厅≤目<目・+石(2—12)幽目++石≤目≤2re上式中局为接触面间的剪切刚度,口‘可以表示为:0‘=COS。(1—2∥I(Kd・B))(2一l3)(2—12)式是一个分段函数,很难应用到系统有限元方程中。由于系统所受激励为周期作用力,可以应用谐波平衡法对摩擦力进行处理151】。谐波平衡法是一种基于傅立叶级数展开的方法,根据此方法将摩擦力进行傅立叶级数展开,并取其合适的前几阶谐波。在保证精度的情况下取一阶谐波,即前两项为Z。=Fo(B)cos0+E(B)sin819篁三皇吐堕堕星垫堡垦塑塑垄坌盟——————其中E(B)--,2。fe.c。s伽=X疗d__BB(8'-sinS"cos8")(2--15)胛)=昙r埘n础护=弓舢一茜)2摩擦面正压力变化时的摩擦力表达式fcz叫)D/2,uI,I.。|j|秤|一rB2n冉一c纱Au图2.14正压力变化时的摩擦力与位移迟滞回线对于摩擦接触面为一斜面的阻尼结构,其摩擦接触面与叶片振动方向成口夹角(0<口<90)。当叶片阻尼器相对运动时,其接触面压力是变化的【521,变化的正压力作用在叶片上,会对叶片频率和响应产生进一步的影响。设阻尼结构间的正压力为相对位移的线性函数,相对位移△“仍然为时间的简谐函数,则接触面间的正压力表示为:N=No+^。Au(2—17)上式中No为初始正压力,丑。为正压力增长系数。阻尼端的摩擦力与位移的迟滞回线可以用图2.14表示,图中A点表示目=0,B点表示0=0+,C点表示0=石,D点表示日=0’+丌,则摩擦力的表达式为:Kd(Au—B)+/.t(No+A。Boos0)OS口<口‘口+≤口<石f=一∥(Ⅳo+A。Bcos0)一Kd(Au一占)一∥(No+旯。Beos0)(2—18)石≤0<0++厅秽++厅兰拶≤2x拱@n七A。BcosO)p・≠cos一1f1一—2,u(No—+AeB).)、(Ka+2∞。)B。同样用谐波平衡法将上式线性化:(2—19)第二章叶片阻尼机理及檄振力分析^。=F。(B)cose+F,(B)sin0(2—20)其中疋(驴昙r厶c。sOdO=姜【三(置。B+2肼pB)O++(2/.aVo一足。召)sin口’+I(KaB+2弘A。B)sin20"lp.).pBl只(耻昙jr/:rsin舭=i2‘1(世。口+2,%B)一Ka曰十{(畅曰+2,以,国cos2口+一(2∥o一亿B)cosa+]另外,叶片阻尼接触面除承受变化的摩擦力外,还承受着变化的正压力,这个变化的正压力对叶片来说,也是一种外部激振力。对变化的接触面正压力进行傅立叶分解得:N=只cos8+只sinO只=委IF(Ⅳo+)tBcos0)cos甜目=2vB只=昙lr(Ⅳo+2pBcos州nodO=。从上面的推导结果可以看出,正压力得变化只是影响叶片系统的刚度对阻尼矩阵没有影响。3摩擦过程中摩擦面发生脱离时的摩擦力表达式前面已经叙述了滑移面间隙模型,式(2一l0)表示在接触面相互摩擦的过程中,有一段时间两接触面是分离的,它们之间没有摩擦力。摩擦面之间的力与位移的迟滞回线与图2.12有很大的不同,如图2.15所示。设两接触面相对振动初始有一间隙D。,振动从A点1(0=0)开始描述,沿卸载线AB到B点(0=0+)(两接触面由最大相对正位移开始反向弹性移动),由B点到C点(0=鼠)两接触面开始反向相对滑移,到C点后两接触面脱开,摩擦力降为零,即为D点,沿DE线两接触面继续反向运动到E点,到最大相对负位移后又开始相向运动,沿ED、线到D’点(0=02)。到D、点后两接触面又开始接触并相互摩擦(加载)到F点(0=03),到F点后两接触面又开始正向相对滑移。它的力和位移的迟2l第二章叶片阻尼机理及激振力分析滞回线是比较特殊的,因此其摩擦力的表达式为:E,-历。缓盈.05图2.15阻尼器间有间隙时的摩擦力与位移迟滞回线Kdr(Au—B、+/2,v0兰0<0’一/2N0+≤0<qf=O01≤0<02(2—26)局(△“一喀)02≤0<岛懋03≤0≤2zc其中㈣os一1(务8・=2rr-岛,03=cos-'(--警-+茜)(2~27,28,29)0’与前面的相10。对上式用谐波平衡法线性化为:'xn=Fc(B)eos0七Fs(B)sin0(2—30)式中Fo(B)=丢f5工c。s8dO=击【也占(sin20++sjn2岛一sjn2岛)】+毒三二【—kB(岛一岛+目+)】+土(,W—KdB)sinOY+三,W(sin0‘一sinB)(2-31)z汀。耳行+i—lKdDs(sin02一sin03)1一』∥矿sin岛万。疗FA曰)=妻r”^sin鲥目=-击KaB(c。s2岛-cos203一cos口+)一当(州一KaB)c。s目’上∥(c。sOi—c。s目‘)+兰彪dDs(c。s岛一c。s02)(2—32)石lZ"。疗”+当弘Ⅳc。s岛+a-占KaB+三(#Ⅳ一KaB)一土弘Ⅳ死…耳冗第二章叶片阻尼机理及激振力分析2.6激振力分析在机组运行过程中,汽轮机叶片将不可避免地受到激振力的作用。激振力是叶片动力响应分析的前提,也直接关系到叶片的安全运行。但在激振力的各阶谐波分量中,影响叶片安全运行的激振力仅是接近叶片自振频率的分量,其它分量虽然使叶片做强迫振动,但是动力响应较前者小很多。因此,降低接近叶片自振频率的激振力分量是关键问题。2.6.1激振源分析作用于汽轮机动叶片的激振力主要由下列因素产生。l结构因素进排汽管;级间抽汽;喷嘴配汽;尾缘厚度;隔板加强筋。2制造安装因素喷嘴节距和喉宽沿圆周分布不均;隔板中分面位错。3气动特性因素叶型表面产生附面层甚至脱流;汽道两端二次流动;激波扰动;动叶片自激振动。在小负荷、高背压运行工况下,动叶片在远离设计工况的大负冲角下运行时,动叶片容易在振动过程中不断地从汽流中吸收能量使振动加剧而发生颤振。4运行因素动叶顶部径向汽封磨损而产生非均匀周向激励;喷嘴出汽边冲蚀削落、动叶顶部围带开裂、叶顶汽封环缺损、喷嘴遭异物打伤等均将产生非稳定激振力。在以上各种因素中,主要产生高频激励的因素是喷嘴尾迹。其余各因素主要产生低频分量。2.6.2影晌激振因子的因素激振因子表示叶片所处汽流场的不均匀性,即汽流激振力的大小。1喷嘴尾缘厚度的影响静叶尾缘厚度与~阶高频激振因子近似成正比关系。随着尾缘厚度的逐渐增大,激振因子也随之增大。因此,合理设计喷嘴尾缘厚度来改善气动状况,对减小激振因子非常重要。2动静轴向间隙的影响第二章叶片阻尼机理及激振力分析动静叶栅之间相对轴向间隙的变化与激振力的关系而言,一般的结论是间隙越大则激振力越小。因为尾迹将随着间隙增大而逐渐拉平。但是(53】的计算结果并不是如此。从图2.16可以看出,先将动静间隙从设计值10mm处增大,则可看出一阶高频激振因子随之减小,与一般结论相符;但当将动静间隙从设计值10mm处减小时,最初激振因子一直增大,当动静间隙减小到一定程度时激振因子开始减小。这说明在设计参数合理选取动静叶栅之间的轴向间隙是非常重要的,因为在不同的范围取值其变化趋势是不同的。g喜囊墨l图2.16激振因子与动静轴向间隙的关系3喷嘴喉宽节距制造偏差的影响在喷嘴的制造安装过程中,不可避免地存在误差,但是误差以何种形式出现也将决定着激振因子的分布。当喷嘴喉宽节距无制造安装偏差时,低频激振因子很低,只有在高频激振因子处出现峰值。当只有若干个喷嘴喉宽存在误差时,其激振因子的分布形式发生变化。人为的将第一个喷嘴喉宽减小,比其余的都小,则低频激振因子增大,且从2~6阶逐渐减小,到高频激振因子处出现峰值。为了更深入的研究喷嘴喉宽在随机状态下的不同形式对激振因子的影响,人为的改变其大小,令其具有一阶阶跃分布的特征,即喷嘴总数的前半段与后半段分别具有不同的值,如图2.17所示。‘此时激振因子的分布如图2.18所示。从图中可以看出,在奇数阶的激振因子及高频激振因子处出现峰值,而其余偶数阶的激振因子很小,几乎为零。用Fourier验证说明结果正确。24:F,第二章叶片阻尼机理及激振力分析图2.17喷嘴喉宽沿圆周图2.18喉宽一阶阶跃分布时的激振因子2.6.3激振力谐波分析动叶片在旋转过程中,将受到以转速为基频的周期作用力的激励,基频激振力的作用周期为丁2哆么。。对沿圆周分布的激励载荷可用付氏级数展开为转速倍率的各阶谐波分量。厂(f)=饨+nT)0=1,2…)巾)=孚+言o。c。sGw)+乩sinG‰r))记4=粤,-。=厢,吼:喀一-生ZOn则丸)=A。+∑以cos(ncaot+q’。)H=I为采用快速付里叶变换(FFT),应用欧拉公式将式(2—36)作变换,(2—36)式可以表示为:厂O)=∑C。・e…‘式卢G||%一一m矽w二:吨一!卜斑显然=第盯4次营q波的振幅为该阶复数模㈨两倍:笙三童生苎堕星垫堡丛塑堑查坌堑A。=2・蚓=√口:+醒采用FFT算法,可以快捷地求得C。的实部和虚部,进而求出a。和b。鲁=Real乜),二箬=Imag(C。)1…71…(2—40)则n次谐波的振幅和相位角分别为:』。:厢=舻辔。鲁a叫铡Keal]2.4—(Real(C,)2+—(Imag(C.)2。\卜。,/将,40,A。,及吼代入式(2—36)即可获得激振力的各阶谐波分量.‘,(f)=focos(coT)・・・吒=。,p。=。,一。=2‘l亡。12;尚c。s(n万)sinp石)(2--43)-・巴2斋与c。s㈦s;n㈤,,2iCO考虑对汽轮机叶片振动影响大的激振力第一阶、第二阶和第z阶谐波分量(n2=29)l,Qz=zQl):focos白r)=Alcoscolr)+爿2COS((D2r)+A:cos∞:T)(2—44)本文首先介绍分析了汽轮机叶片阻尼的分类,以及叶片各种阻尼的机理。研究表明,在叶片阻尼中占主要地位的是叶片的材料阻尼和界面摩擦阻尼。界面干摩擦力是阻尼叶片的一个重要外部作用力,针对实际应用中对摩擦面的受力情况和所考虑的侧重点的不同,对其可以建立各种类型的滑移模型。在此基础上,可以推导出于摩擦力的表达式,并应用谐波平衡法对其进行线性化处理。然后,本文对产生汽流激振力的激振源进行了介绍分析,讨论了影响激振因子的主要因素,并对激振力进行谐波分析。在以后各章中,将取含第一阶、第二阶和第z阶谐波的激振力进行叶片动力学研究第三章单自由度模型分析第三章单自由度模型分析3.1概述叶片是汽轮机的重要零部件,因承受气流周期性扰动力作用而振动,加之复杂恶劣的工作环境,常常会造成叶片失效。振动是影响叶片安全的重要因素,在很多’隋况下叶片失效都是因为叶片异常振动或承受共振应力带来的高周或低周疲劳而引起的。目前,世界上的汽轮机向大容量化、高速化发展。汽轮机的大容量化使得叶片需要比以往承受更大的气流力、离心力和更高的工作温度,因此对叶片的安全性提出了更高的要求。摩擦阻尼叶片具有良好的减振效果H“,对其动力学行为的研究工作已经有很多。Berthillier“”用多谐波频域分析结合C-B的局部模态综合方法,研究发现强迫力中的高阶成分对叶片的振动响应有很大的影响。季保华¨“实验研究了摩擦面压力、接触状态和激振力对叶片动态特性的影响,结果表明摩擦对减振是有效的。丁千¨”研究了干摩擦力对自激振动的影响,用平均法结合奇异性理论发现了自治系统的幅值跳跃和滞后。本文研究汽轮机干摩擦阻尼叶片在多谐波激励作用下的振动。用平均法””求出系统第一阶谐波和第二阶‘谐波的主共振稳态响应方程,分析了系统各参数与响应之间的关系,特别是考查干摩擦阻尼器参数(主要是阻尼端初压力值)抑制叶片振动的效果。所缛结论对工程应用有一定的指导意义。3.2阻尼叶片模型分析可以用质量弹簧系统来解释阻尼叶片的动力特性。叶片用质量块///来表示,阻尼器用无质量的弹簧和摩擦元件表示,如图3.1所示¨”,当叶片受稳态激励focos(coT)时,叶片的运动方程为:・肌膏+cic+kx=foCOS(mT)一无(3—1)其中正为摩擦力。在正压力作用下,两结合面相互压紧和贴合,结合面凸凹处相互嵌入,在较小的切向拉力Z作用下,两结合面之间先有一弹性至塑性的微变形,在这个阶段,叶片阻尼结构接触端之间没有发生滑动。当^增大至一定程度,两接触面之间“接触层”被剪力拉断,叶片阻尼结构间便开始相互滑动。如果将阻尼端正压力Ⅳ在叶第三章单自由度模型分析片振动过程中的变化考虑进去(设为线性变化),则摩擦力^与相对位秽x形成如图2.14所示的迟滞回线l54】,其关系如式(2—18)。在机组运行过程中,汽轮机叶片将不可避免的受到激振力的作用。激振力是叶片动力响应分析的前提,也直接关系到叶片的安全运行。但在激振力的各阶谐波分量中,影响叶片安全运行的激图3。1阻尼叶片物理模型振力仅是接近叶片自振频率的分量,其他分量虽使叶片做强迫振动,但动力响应远比前者小1591。因此,我们只需降低接近叶片自振频率的激振力谐波分量。考虑对汽轮机叶片振动影响大的激振力第一,第二和第Z阶谐波分量(Q2=2f2l,Qz=Zf21)则:focos(roT)=AlCOS(0)lr)+彳2cos(c02丁)+4cos如乞T)(3—2)3.3平均法解方程叶片的运动方程为:m2+西+舡=f0cos(砑)一/:rG)(3—3)引入无量纲化变量:卜co.t,X=言芦志一2詈Pnm西。x“”‘删乒赤24mk一悟E=鲁foVm则原方程变为:童+2劈+J=Vocos(Q丁)一只∽)(3—4’其中Focos(f】T)=Flcos(1r)+Ecos(sq2f)+疋cos(f2zf)F=丢,只伍)=砉为应用平均法求解一阶或二阶谐波主共振响应,即Q,*l或Q:a1的情第三章单自由度模型分析况,令Q,=1—60"(i=1,2),占为小参数,盯为调谐参数。在相对小项前标以小参数占,将方程(3—4)写成:戈+X=dEcos(f!。t)+Ecos(n:t)+Fzcos(f2:t)一只(X)一2斟】(3—5)取a,0作为新变量(分别代表振幅与相位),取以下变换:X=acos彬,工=一asin妒,妒=Q,,+曰(3—6)将(3-5)化为标准形式:=一S恤一出一讲卯一卉=1一卟y三d,∞矿协当:黜y凋:宠凇Y曾l臼='9+占K(f,,19)+占2%(r,,圆c,吲‘’。4j瓦dy2蹦(y湖¨2轴搠托誓w艄占’I警喝(y,回+9222("m狮们占);3-9)““其中五,E,Z。,z:不显含t,U1,U2,K,心,r+,Z+为t9的以2z为周期的周期函数践+占2E+93y,+eO西Uls百aUt(碱+s2E+93y*)+占等(髓I+ce222+c3Z*)+占:孕+占z_OU2(碣+占:匕+占,r・)Ot钟膏等(函+,CZ2+e3Z+)=砒∥¨应,∥z:∥z+s婴占警(硝膏E+83y*)+占等(理。∥z2+E3Z"Ot)OVd拶+若:掣+占z婴(蹦+s:E+占,J,‘)Ot洲¨2斋(露1+占222+e3Z*)=硼托2衍^..————————————————箜三兰一皇旦自廑蕉型坌堑(3—10)其中丸,簖,蛾,硝为(3—7)式右端对s展开的泰勒级数的系数。令等式两端占一次方的系数相等,得到:z+等=唬互+鲁=簖为了使K,Z,满足不显含t的条件,则K2芴1舷d∥ZI=去肼d∥U1=f(丸一r,)dt+Uj。K:j(硝一Z,)dt+巧。・墨+等“等z。+警=西c缈d.V’疗,¨z:+等巧+嘉五+鲁=衍砂d拶‘研~E。去M~等一一等动d妒z:2去胁一等z一嚣z.脚%=去f(一a砂u,.~OaU拶,Z。一‘矽也。K2去f(卜a砂v,r.,一嚣Z】吲西叱X一J∽叫m}8、一=粥@口)砂一以棚一出●^=髓@口)(3~12)(3—13)(3~14)(3—15)(3~16)(3—17)从(3一11)可以求出Ul,K:再令等式两端占2的系数相等:为了使匕,Z2满足不显含,的条件,则:一次近似解是:二次近似解是第三章单自由度模型分析Ix=acos(Qr+0){Ⅱ=Y+eUl(y,t,拶)I曰=口+cVl(y,t,口)而李:占I+占:Edt—d—O:钇.+占2Z。击1(3—19)。3.4第一阶谐波的主共振(n.。1)令fll=l一盯,盯为调谐参数。且Q2=2Ql,Q:=ZO.。由平均法解方程的步骤nJ知,丸,靠,缟,行为(3—7)式石端对占展开的黍勒级数的系数丸=一丢_(sin(2Q。卜P19)+sin(卵)一三厶(sin(3Q.f+19)一sin(Q。f—t9))一昙厶(sin((Z+1)Q。f+口)一sin((z一1)Q,f一.9))一2印。y(三一三cos(2f毛t+2跏切∞簖:仃一去z(c。s(2Qlf+占)+COS(19))一丢^(c。s(30】r+19)zYZv+c0S(f2lt-,9))一专^(c。s((z十1)f21t埘+cos((z一1)qt-3”一髀1sin(2Qlr+23))-t-fB(t)办=【三(2卢+口)sin(2Qit+20)+(∥一a)sin(f1.f+印一笋1。+击coscf,jt-6I沪圭聃掣舞竽+髀・cos(2Q!r+2印】[乏}_c一三正(一击c。s(3fld+3):塑三主璺自由匡垡型坌堑+—c—os—(—(-—1—+—Z:;)ik'1;:il—t=旦)一.写:1:_ycos(Qlf+19)(2∥+口)c。s(Q1f+t9)(-17+z)Ql’2Ql一…。7、1十2(∥一口))+:1争sin(2nlf+20)+.2a_A(卢一口)(y-a)】十2(∥一口))+j争8i“(2n1‘++’v疗(卢一口)】+[一吉^(c。s(2n,f+拶)+c。s(t9))一丢^(c。s(3Q。r十印+c。s(啡棚)一言五(c。《z+1)Qit+占)+c。s“z咖啡一劝+y(2fl+a)cos(2Qlt+2国+y(∥一a)cos(Qlt+t9)一2《茧m1sin(2f2tt+2,9)][一击/iSin(2Q・川)一万1拟面1sm(3fkt+,9)+等sin(f毛t+8)+圭如(2qt+25){p一石1邶叫+高sin((_1+zH卜劝+笔等(sin(2哪Ⅵ趴Qpl9)+击s-n(叩卅))一fZ(面击sin((z+1)Q”19)(∥2fl—Tir-2≯'3+sin(2%))]衍吖专他os(2脚旧+cos(劝+可1,2(c0S(3啡删+cos(n21t-圆)+。专^(c。s((z+1)Q—f+.9)+cos((z—1)Q,卜-.9))】一三厶(一赤c。s((z十1)f≈t+8)+高c。s((一l十z)f毛t-8”【击舾s(2Qp即三斛击c。s(3Q,H印+击c。s(f毛t-8))~击旭∥刊cos(f2,t+8)2一击邶训cos(f)tt+ou)+三㈣n(2f),t+28)32第三章单自由度模型分析+嚣(卢刊(y-口)】+[专舶in(2QIf+即sin(跏~万1^(_sin(3f2,t+oa)1+8in(Ql卜8))一寺^(一8in((z+1)Q1HJ)+8in((z一1)f2,t-t9))一(2fl+a)sin(2f2It+2,9)一(卢一a)sin(f!lf+口)一2孝Qlcos(2f21r+20)】[一击^血(2f2tt+oa)一专^(击sin(3f2tt+8)+击s血(f21t-t9))一专厶(赤sin((z+1)Q,r+回+高sin((-I+Z)f-t-8))+击(2∥托)哇sin(2f2一t+2,9)+f2tt+,9)+壶(∥训sin(QlⅢ)+丢细s(20tt+2¥一j1辟一石I,够刊毪一2%+sin(2㈨)J由式f3—121得军II._=一lflsin(19)一舌QtY+A(y)z.:仃一去工c。s(19)+塑‘y,其中Y和口分别是U和0的一次平均值f陟cos劬)【O,石]y∈(0,口)E={aycos(妒)+a(fl一口)【0,仍】Y∈(口,a+6)【彦cos(伊)一(Or"一f1)(y一口)【仍,石】彤)2去r‘f.(ycosg)sin‘丫)dy=一去f2f.(ycos妒)sin(州PB(y)2去r4^(yeos妒')c。s缈)d妒=击r“^(yeosq,)。?8(们咖其中仍=arccos(等一1】,辔“口是图2.14里Bc与x轴的夹角,tg-I卢表1F(rCOS们、爿fvl和肼v、积分结果令y=2x—p,有:伊为自工轴起逆时针转过的角度。考虑到滞后弹性力F(砷的对称性,积分在区间[o,石】进行,积分结果见表1。舳2.14里AB与x轴的夹角一紫一争伊华。。第三章单自由度模型分析yF(ycos伊)flycos(妒)爿(_y)L曰(,)(O,Ⅱ)O,譬y【口,o。)aycos(9)+口(口~口)一2a(fl一口)(y-a)寺(脚)(等却,#ycos(力一(口一f1)(Y一口)ny+sm(2p3)j田虱L3一l3j口J以r付出:u一三一(一击c。s(2哪岫“n(黝)一言兀(一击cos(3Qit+,9)+击cos(n,t-8))一三^(一志cos((z+1)flit+国+赤cos((弘1)f21t-00))一焘(2∥怕)cos(哪瑚)2一言(∥叫)cos(Q”回_2扣。哇f一击sin(2qt+2,9))+三石sin(O)r+亩n.,+2at(∥一口)(y-a)ZV石K=一参击s叫2Qt¨回+c。s(印)一万fz酉/1sin(QIf-跏一岛(盎sin((z+1)啡埘+赤sin(3叩十∽+击8in((z-1)Ql’一印)+j毒(口+2∥)(主8in(2Q1‘+2_9)+Q1,+固+等sin(Q1f+趴丢细s(2f毛t+2,9)+万fltcos(垆丢厣一去(p训(差却3+Si卿3))由式(3—15)得到艺,Z2,其表达式非常复杂,见附录。我们发现,在多频激励系统发生主共振时,平均法的一次近似解中只体现一阶激励项的影响,但二次近似解中包含了各个激励项的影响。给定计算中系统参数为:∥=O.01,五,/k=10,kd/k=1,Ⅳo/A=80,Z=44,占=o.I,孝=0.5,Z=5.5,厶=0.01,厶=O.1。代入二次近似方程的定常解:第三章单lit由度模型分析a_y,=era(y,'9)+占2K(y,口)=0at1芋=厦l(y,回+9222(y,_9)=0al(3—21)计算得到如图3.2(a)所示的幅频响应曲线(实线),点是数值精确解。可以看出,用平均法得到得的共振区内的周期解与数值解是相当接近的。但远离共振区,则有较大差别。其他参数固定时,振幅随着一阶谐波幅值^增大而增大,但随着阻尼系数孝的增大而减小。本文的主要目的是研究千摩擦阻尼器抑制叶片振动的效果。在设计和安装时,阻尼端初压力值是较易调节控制的参数。为此,考察初压力值Ⅳ0对共振响应幅值的影响(ka/k=O.5和1),见图3.2(b)、(c)。随着k。/k增大,不同正压力下的频响曲线移开较明显,且从图上可以看出k。/k较大时频响曲线峰值有所下降。在Ⅳ。/.厶=350时,k。/k取不同值时曲线峰值波动较大。图3.2(b)中,当Ⅳ0/.矗从O逐步增加至1OO的过程中,幅频响应曲线的峰值逐渐下降,表明干摩擦阻尼器的阻尼机制起到了消耗系统能量的作用。而当No继续提高后(从帆/兀=120之后),阻尼器接触面的相对滑动能力会随之减弱,甚至完全不能滑动,因此对系统能量的消耗反而减小,响应幅值又会回升。至No/‘=350时,因为a=7已经大于此时的响应峰值,故阻尼器接触面已完全不能滑动。而且在眠/工逐渐增大的同时其峰值向右偏移,说明同时也使得系统的固有频率有所提高,共振时间后移。从这些现象可以得出以下结论:在给定的参数范围内,阻尼端正压力对响应的影响很明显。阻尼端正压力影响阻尼叶片的频响曲线,即影响叶片的频率、响应,并且阻尼端正压力存在一个最佳的值,在此值下阻尼叶片的阻尼最大,减振效果最好(与文献【60】结论相同)。kd/k=O.5时,最佳压力值为“/厶=120;ka/k=l时,最佳压力值则为Ⅳ。/厶=220。通过图3.2(b)、(C)的幅频曲线比较还可以看出,随着b/七的增大,频响曲线峰值有所下降且偏移量增大。图3.2(d)是考虑没有阻尼器的情况,振幅随着阻尼系数掌的增大而减小,即随着孝的增大幅频响应曲线趋于平滑,随着,:的增大而增大。与有阻尼器情况下幅频曲线的比较可以看出,系统的共振频率位于1附近,共振时间没有发生漂移。由图3.2(a)、(d)可以看出阻尼系数孝,一阶谐波幅值工只是影响峰值的大小,而k。/k和Ⅳ0/矗还会对频响曲线的偏移量有所影响。第三章单自由度模型分析(a2)她5风夕≮~……~A一~,I(b1)一酝-I-NO=O2-NO=IO图3.2主共振幅频响应曲线及各参数的的影响笙三童璺自由崖蕉型坌堑3.5第二阶谐波的主共振(Q:“1)在工程实际中,有很多汽轮机某级叶片固有频率与第二阶激振力的频率接近的情况,故以下考虑Q:*1时的主共振情形(Q,=丢。:,Q:=丢盈:,。。令Q:=1一盯,盯为调谐参数,推导出丸,硝,破,群为(3—7)式右端对s展开的泰勒级数的系数:丸=一iizsin(3f)2t+∽一ilzsin(吉Q:f+t9)一i1^sin(2Q:f+圆一圭^stn(。)一丢艿s;nn=t+a)+圭,zs;哔Et2t--,9)+sin(f12t+圆力O)一争Q2+争n2cos(2f)2r+28)蝣=口一六zc。s专Q:H・回一六^c。s(兰Q2t--j+回一六五c。s(.9)一专^cos(:Q:H印一专厶cos(丝笋Qz卜印一专厶cos(堡笋Q:H圆+co“叫枷)等母:sin(2f如t+2,9)警一迦(-等1警奇加,,!兰!!遂竺竺+苎::蔓!!!竺+—cos(2—fbt+O)+(1+音z)Q2+{z)Q2…2yc。s(Q:f+口)2一去J,(p一口)c。s(Q:H-圆+i1争sin(2fkt+2圆+2。a一石t(,8一口)(y一日)]+【一圭一c。s(三Q:f+19)一三zV石ZZZc。s(丢Q:r+19)Z一三五c。s(2Q:,+t9)一i1五c。s(19)一j1以c。s((1+lz)n:f+回一委六c。s((一1+z)n:f一回+(La+口)yc。s(2Q:r+2。9)+y(fl—cr)第三章单自由度模型分析4-。。8(Q2H口)一2姻z8in(2fa2t+2,9)]【一—s—in(寺矗‘2_t一:19)1工对n呸f22t+8)一三厶!!呈!磐21±旦1fzsin((一l+三z)Q扩回一专竺4归22y(-l+1z)Q22y型+警c言sin(2f,2t+,9)+扣埘,(1+Z)fl。f22、42。。‘2‘。。’’2+等Sin(f22t+,9)+争s(2Qzt+20a)一互1卢一石1r(fl刊(∥2fl一#i一2%+sin(2圳】-群=[专/i(cos哇Q:r+'9)+c。s哇n:r+.9))+互歹1^(c。s(t9)。+cos(2Q:f+圆)+云_厶(c。s((一1+吉z)C2zt-,9)+c。s((1+i1z)Qzr+回)】lf1cos(3fa2‘+J)ft3——^——cos(12f+—9)Q.^—14——f2cos(2—Ca2t+,9)Q2l,zc。s((一l+互1z)Q2卜口)2(=l+圭z飓一去(2∥+咖cos(k'12l+8)2--击(fl一咖c。s(f22t+¥)+吉鲫n(2f22t+2印+丝yz(∥训(y叫】+【圭删n哇Q:,+圆ZZZ+丢埘n(3f2zt+O)+i1^sin(2Qzt+)+丢龙sin(J)一三厶sin((小2Z)f22t-t9)+j1厶sin((1+吉z)Q:f埘一y(2fl+a)sin(2f22,+2t9)一_y(∥-a)sin(fj2r+.9)一2c蛐zcos(2f12t+2州一去sin(圭呦呐VIZ,Z一击s;n毫Q:M,一彘sincz叼瑚,一可靠第三章单自由度模型分析sin(一l十三z)Q:,一.9)一互歹靠sin(1+互1z)Qz,+19)+考茅c扣2啄啪+扣蜊+警sinc印+¥+专细s(2Qzm旷专肛一砺1舻刊(筹一2%“n(2㈨)】由式(3—12)解得:_:一昙^sin(.9)一尊o:+4(力(3—22)Zz。:盯一:一厶c。s(19)+』堕堕ZVV其中Y和.9分别是“和0的一次平均值。同样由式(3一13)得到:u=去一c。s(3[.22t旧+cos(1D2t+oa)+去,2cos(29)2t+)+击蠹fzcos((一Z?蚴+圆一丢矗厶cosc罕邬讷一盍(2卢圳cos(Q:H∞2一毒(∥叫cos(Q:f呐+三争sin(2n:f+zo)+丝(卢一口)(y-a)ZV厅忙专岳s抽扣旧岳咖扣例一击删2Q2t+¥一立2y(丢忘sin(坚塑2叩+圆+杀Z啬2)n、’sin(堡望2叩咖)+、、(Z+2)Q2‘。(一‘“2去(口+2∥)哇sin(2f22t+28)峨r+圆+譬sin(毗印+三孝cos(2Q2t+2¥一吉芦一去(∥一a)(云笔一2%+sin(2%))由式(3—15)得到E,z:非常复杂的表达式,见附录。取系统参数为:∥=0.01,以/k=10,b/k=1,No,厶=80,z=44,占=O.1,孝=0.5,^=5.5,/i=0.01,厶=o.1。代入二次近似方程的定常解害=磷(儿卵+占2rAy,O)=。警=函(弘圆+s222(弘¥=o(3--23)计算得到如图3.3(a)所示的幅频响应曲线.(实线),点是数值法得到的精确解。可以看出,用平均洼得到得的其据阿内的周期锯与斯信艇旱相当第三章单自由度模型分析接近的。但远离共振区,则有较大差别。其他参数固定.时,振幅随着一阶谐波幅值厶增大而增大,但随着阻尼系数善的增大而减小。取k。/k=O.5和1.0,计算得到不同初压力值的频响曲线,如图3.3(b)、(c)所示。可以看出,随着k。/k增大,不同正压力下的频响曲线移开较明显,且k。/k较大时频响曲线峰值有所降低。与第二阶谐波主共振的情况类似,在给定的参数范围内,阻尼端正压力对幅频响应曲线的影响很明显。如果阻尼端初始的正压力较小,那么系统的共振频率也较低,系统响应接近无干摩擦系统的响应特性。随着正压力的增大,阻尼器接触面之间阻碍相对滑动的能力也随之增强,使得系统固有频率提高,阻尼对系统能量的消耗增加。正压力增大到一定程度后,阻尼对系统能量的消耗反而降低,系统的响应接近于接触面间无相对滑动的线性系统响应,由此可以推断,一个最优的摩擦接触状态一定存在,在此状态下,系统能量消耗最大,阻尼减振效果最好。图3.3(d)是没有阻尼器的情况下频响曲线的形状,振幅随着阻尼系数孝的增大而减小,即随着孝的增大幅频响应曲线趋于平滑。与有阻尼器情况下频响曲线的比较可以看出,系统的共振频率在1附近,没有发生漂移。由以上的计算分析可以看出,各系统参数对第一阶谐波和第二阶谐波的主共振幅频响应曲线的影响是一致的。3.6结论本文研究汽轮机千摩擦阻尼叶片在第一,二和Z阶谐波激励下的振动,考查干摩擦阻尼器参数(主要是阻尼端初压力值)抑制叶片振动的效果。首先用平均法求出系统的第一阶和第二阶谐波的主共振稳态响应方程,然后分析主共振响应与阻尼系数,阻尼端正压力初始值,摩擦系数,阻尼器刚度及各阶激励的幅值的关系。结果表明,阻尼器刚度增大会使阻尼叶片的频率升高,并降低阻尼叶片的响应。阻尼器接触面正压力的变化会显著影响阻尼叶片的频率和响应幅值,即频率会随着正压力升高,而对于幅值,则存在一个最优的接触压力,使得幅值取得最小值。通过比较知道,各系统参数对第一阶谐波和第二阶谐波的主共振幅频响应曲线的影响是一致的。本章结果说明,要获得较好的减振效果,在实际设计中,阻尼器的刚度、接触端的正压力、粗糙度及阻尼叶片的外部激振力都是必须要考虑的因素。第三章单自由度模型分析1-f2=5.5(81)一.A(a2)<b1)(b2)(C1)(d1)(d2)图3.3主共振幅频响应曲线及各参数的的影响4l第四章弯扭耦合模型分析第四章弯扭耦合模型分析4.1前言在涡轮机设计、制造和运行中,对于动叶片都给予特别的重视,并为此开展了大量的研究工作。尽管如此,叶片毁坏所造成的停机事故还是不断发生。据统计,叶片毁坏所造成的停机事故占国内外全部停机事故的30%以上。如果将叶片毁坏原因分成制造缺陷、加工装配缺陷和材料缺陷,那么设计原因是主要的。对于固定式发电机组来说,随着单机容量的直线上升,停机事故造成的损失也就急剧增加。这就给叶片设计工作者提出了更高的要求:在设计阶段中更精确地确定叶片振动的特性。为了减少流动损失和增加叶片刚性,在涡轮机中,往往采用围带和拉筋将几个叶片连接成叶片组。几个叶片与围带和拉筋组成的振动系统具有多样的振型和复杂的频谱。对于等截面叶片组振动计算,普遍地采用系数图线法。对于变截面叶片组振动计算,早在五十年代,Prohl就曾经提出过围带叶片组的振动计算方法呻“,它忽略了切向与轴向弯曲振动的耦合,只将轴向弯曲振动与扭转振动耦合在一起,这就局限它只适用于截面的最小主形心轴与轴向重合的叶片。同时它还忽略了剪切变形、转动惯量和离心力的影响,因此无论将它用于短叶片或长叶片,计算值都不够精确。Deak等提供了曾在美国AlliS—Chalmers公司多年应用的汽轮机末级叶片组的振动计算方法哺”,它将切向振动、轴向振动和扭转三者耦合在一起,并考虑了边界叶片和中间叶片的差异以及离心力的影响。它既可用于围带叶片组,又可用于几排拉筋连接成的叶片组。由于没有考虑剪切变形和转动惯量的影响,因此用于高压级叶片时受到限制。Thomas等应用有限元法计算叶片组的自由振动频率¨”,。用Sturm序列和反迭代法求矩阵的特征值和特征向量。可惜它只是考虑一个平面上的弯曲振动,也忽略了剪切变形和转动惯量的影响,几排拉筋连接成的叶片组的情况也没有考虑。林俊灿”"采用传递矩阵法,它将叶片作为变截面梁的非对称弯曲,将扭转振动与两个平面内的弯曲振动通过拉筋耦合在一起,并作以下改进:考虑了叶片和拉筋的剪切变形和转动惯量的影响:合理地考虑了离心力的影响;在动态下计算叶片与拉筋之间的作用力和力矩,这样既提高了它的计算精度,又扩大了计算机程序的通用性。谐波平衡法是求解非线性振动问凰常用的一种近似解解析法。与其它各种摄动法不同,谐波平衡法的应用不限于弱非线性系统,此外,它第四章弯扭耦合模型分析的求解过程归结为代数方程组的求解,而其它方法则需求解微分方程组或积分一微分方程组。本文考虑了汽轮机阻尼叶片弯曲和扭转的耦合,用谐波平衡法对建立的方程进行了分析,同时考虑Q。=q和Q:=03。的情况下模型的动力特性,特别是分析了参数对幅频曲线的影响。所得结论可以指导工程实际。4.2弯扭模型的建立考虑如图4.1所不的弯曲和扭转耦合模型,叶片用质量块历来表不,阻尼器用无质量的弹簧和摩擦元件表示,图中采用了两个并联的阻尼器[651,当叶片受稳态激励focos(COT)和力矩m。cos(coT)时,叶片的运动方程为:Im茧+c-2+缸+正1+五2=f0cos(cot)1L≯+q痧+红妒+六1il一厶2il=m0引入无量纲化变量:xcos(cot)(4_1),砜r㈡2言,拈乓mco;l,n,=罟nl,萏2赤mk∽,2佶m,弘詈,国√Vp万2景焉2焘一:=居吒=鲁一盟mo—i(Onl则方程(4.1)变为:f置+点力1+xI+‰+E2=Focos(QT)1攒:+劈:鼻:城岷寺一c:妻:M。。。。(叩在机组运行过程中,汽轮机叶片将不可避免地受到气流激振力作用。激振力是叶片动力响应分析的前提,也直接关系到叶片的安全运行。在激振力的各阶谐波分量中,影响叶片安全运行的激H_2’振力一般是接近叶片自振频率的分量,其他分量虽使叶片做强迫振动,但动力响应远比前者小[591。因此,我们只需考图4.1弯扭耦合模型虑接近叶片自振频率的激振力谐波分量。塑婴童变塑塑盒堡型坌堑——————考虑对汽轮机叶片振动影响大的激振力第一,第二和第Z阶谐波分量(Q2=2Ql,Qz=ZQl):^cos(mr)=/ieosh,)+五cos如k丁)+五cos(cD:T)m。cos白r)=mtcos(mtT)+m2cos(c02r)+m:cos(co:T)(4—3)所以有F0cos(flT)=Ecos(QJT)+疋cos(f22T)+疋cos(f2zT)Mocos(QT)=Mlcos(f2lT)+M2cos(f22T)+Mzcos(f2zT)(4—4)F:五,材.:旦‘‘/omo^。,无:为接触面摩擦力,如果将阻尼端正压力Ⅳ在叶片振动过程中的变化考虑进去(设为线性变化),则摩擦力厶与相对位移z形成如图2.14所示的迟滞回线[5”,其关系如式(2—18)。用谐波平衡法对摩擦力进行处理:在保证精度的情况下取一阶谐波,摩擦力经傅立叶级数展开取前两项得到:{,={cCOS3七{ssin8一fs。Scos.a+击ms哟=扣去主,“叫’上式中正=昙r^c。s-射口=昙[丢k。s+z肛,s扣++‘(2∥。一定一s)stn占++I(KaS+2以s)sinz.9+一吉删捌正=昙r^sin黝拶=昙哇敝。s+2越,s)一亿s+;kas+2肼,s)c。s2占+(4—6)一(2掣Vo—KdS)cos,9’】故加鲁(”≯1)+去(毫宁lz),工z=每“一弘1)+去(南%2—ff.1(X1+妻五)2去瞬辱五畦鲁≯+%墙墨+i1百Lmo丽mo删‘c。吲’=QIXI+Gx2+Q3墨+奶名244第四章弯扭耦合模型分析%=—●々/:,2(z。一L.X2)2去唼卑午三等詈¨+面f*z。.丽fo—t.i1百Lmo丽mo鼢】(4--8)=R1xl—R2X2+月,墨一R4膏2其中:92惫,Q=三去鲁等詈,g=习赢fsl两,幺=去圭等最耻惫如=三去鲁鲁詈,岛2习萧is2画心2去三争焘方程组(2)可以写为:Xl+晶工l+X1+Q1xl+Q2x2+Q32l+Q4力2+尺lXI—R2X2+R3墨一R4只2=Ecos(nlT)+Ecos(f22r)+rzcos(f2z丁)占2碧:+鸳:j:+卫:+(Qlx.+Q:X:+Q3X。+Q.2:)iL一(月。X1一Rz石:+皿j.一R。膏:)姜=M1cos(f2lT)+Ⅳ2cos(f22T)+Mzcos(Qzr)把方程(4—7)简化后得:墨+NIXl+Ⅳ2五2+Ⅳ3置+Ⅳ4肖2=Ecos(f2Ir)+最cos(Q2T)+Fzcos(Qzr)d222+只x2+尸2xI+P3X2+只zl=Mlcos(Olr)+^彳2cos(f22丁)+^如cos(f2zr)‘(4—11)其中:Ⅳl2最+Q3+R3,N2=Q4一R4,N3=1+鸟+R1,N4=Q2一R2鼻2鹄+Q4詈屿专,只=Q3砉一玛孝,B=1+92iL蝎喜。只=QI詈-R1专由谐波平衡法可以设解为:XI(T)=Alcos(92l丁)+日sin(Ol丁)+A2cos(922r)+B2sin(f22,)。+Azcos(f2zT)+Bzsin(QzT)XAT)=CIcos(f11r)+DIsin(nlr)+Czcos(f】2丁)+D2sin(Q2r)+Czcos(QzT)+Dzsin(f/zT)把启。(丁),置(r),岩:(丁),j:(r)求出代入方程组(4—7),使等式两端同次谐波的系数相等,整理得出各阶谐波分量的方程组。蔓婴童变垫塑盒堕型坌塑——4.3考虑一阶谐波Q,=1(其中有09n,=%,)此时有Q:=2Q,,Q:=Zn,,带入式(4—12)用谐波平衡法进行分析。4.3.1用谐波平衡法解方程由一阶谐波系数相等得到:~AtQ;+3Il蜀QI+N2DjQI+N3AI+Ⅳ4cI=E~BlQ;一NlAlQl—N2c1QI+N381+Ⅳ4D1=0一占2clQ:+只DIQJ+尸2BIQl+Bcl+只4=M一万2D,Q;一只c1Q,一只01Q。+忍Dl+只Bl=0由二阶谐波系数相等得到:(4—13)-4.42Q;+2NIB2Ql+2N2D2QI+Ⅳ3—2+Ⅳ4c2=E一482Q;-2NlA2nl-2N2C2Ql+N382+N4D2=0—462C2Q;+2只D2Ql+2只呸Ql+只c2+只爿2=M2—452D2Q;~2Ec2Q1—2马A2Ql+只D2+P482=0由Z阶谐波系数相等得到:(4一14)一z2AzQ;+zⅣI岛QI+zⅣ2D2QI+Ⅳ34z+^r4cz=兄一Z2BzQ;一zⅣIAzQI—ZN2CzQI+Ⅳ3易+Ⅳ4Dz=0—6222CzQ;+z只DzQl+zlP2BzQl+马cz+只一z=Mz—d222DzQ;一zBczQl一zP2彳zQl+BD2+只丑2=0各个系数的表达式很复杂,见附录。(4一l5)由方程(4—9,10,11),可分别解得:4l,B1,C1,Dl;A2,B2,C2,D2和4z,丑2,cz,Dz。故可以得到X。(,),五(,):XI(丁)=Atcos(QI丁)+Bjsin(nl丁)+A2cos(29/1丁)+B2sin(2QIT)+彳zeos(ZQlr)+Bzsin(ZQI丁)(4—16)工2(r)=C1cos(Q1T)+D1sin(Q1r)+C2C03(2Qlr)+D2sin(2fjlT)+Czcos(ZQlT)+Dzsin(ZQlr)给定系统参数为:m=1,k=1,也=1,J;=1,兀=1,珊o=1,占=1,∥=O.02,≈d=0.5,旯。=10,No=6,Z=44,孝=O.5,E=2,F2=O.01,疋=O.1,M1=2,村,=0.OI,M,=O.I代入方程(4一16)计算得到激励力幅值和阻尼器参数对叶片幅频曲线的影响,见图4.2第四章弯扭耦合模型分析(a1)(b1)(c1)(a2)(b2)(c2)47第四章弯扭耦合模型分析(e1)(e2)(f1)(f2)图4.2各参数影响下的幅频响应曲线第四章弯扭耦合模型分析由图4.2可以看出各个主要参数对幅频曲线的影响是很明显的。当其他参数固定时,振幅随着一阶谐波幅值鼻、M。增大,随着阻尼系数茧,岛的增大而减小。但在给定参数下,丘的幅值总是比X.的值大。为了验证剪切刚度k。对响应并,,z:的影响,给定kd两个不同的值,k。=0.5,“=O,8。在不同k情况下取了3个不同的初压力值,结果如图4.2(e)、(d)所示,如果k。,增大.不同正压力下的频响监线移开较明显,且从图上可以看出k。较大时频响曲线峰值反而下降。在图4,2fd)中,当No=6时频响曲线的,峰值较高。当No=30时频响曲线的峰值降低,但是当N。=60时峰值回升很多,同时其峰值向右偏移,表明共振时间滞后。从这些现象可以得出以下结论:在给定的参数范围内,阻尼端正压力对响应的影响很明显。阻尼端正压力影响阻尼叶片的频响曲线,即影响叶片的频率、响应,并且阻尼端正压力存在一个最佳的值,在此值下阻尼叶片的阻尼最大,减振效果最好。图4.2(e)、(f)同样显示随着k。的增大频响曲线峰值降低,且幅频曲线向右偏移。且E,M,的值增大使得幅频曲线中出现了两个峰值。4.3.2通过数值解进行比较图4.3~4.7中的点为数值解,线为解析解,显示出解析解的准确性。图4.3第四章弯扭耦合模型分析图4.4图4.5图4.6图4.7第四章弯扭耦合模型分析4.4考虑二阶谐波Q:=1此时有Q,=三Q:,9。=吾Q:,带入式(4—12)用谐波平衡法进行求解。4.4.1谐波平衡法解方程由一阶谐波系数相等得到:叫譬+NlBl-鲁-+Ⅳ2DI-譬--+N3-41+Ⅳ4CI=鼻喝等一MD22N:cl导+Na马+Ⅳ4Dl=。∥c1垡4+Bq孚+P2BI-譬-郴心mI=M一占2DI譬一Ecl导一只4粤+B日+只骂=。由二阶谐波系数相等得到:一A2D;+Ⅳ1占2Q2+Ⅳ2D2Q2+Ⅳ3爿2+Ⅳ4c2=E—B2Q:一NIA2Q2一N2C2Q2+州B2+_Ⅳ4D2=0一万2qQ;+PID2f22+尸282Q2+Bc2+只一2=M2一J2见Q;一p,c2Q2一只彳2Q2+602+只B2=0’一;也Q;+ZNlBzC22+ZN:见Q:+M如+N4Cz=兄一百Z2BzQ22Z。NIAzQ2-ZzN2CzD2+N3Bz+N4D1=o搿等c捌+导鼻见吣詈B瞩+BCz+驰也一占2百Z2见Q22一Zpicz呸Zp,彳zQ2+只。z+只Bz=0由方程(4—17,18,19)可分别解得:A1,马,C1,Di;A2,B2,C2,D2和Az,Bz,Cz,Dz。故可以得到z.m.丑∽:5l(4—17)(4—18),(4一19)第四章弯扭耦合模型分析x,(D=4co《譬r)+Blsin(号-r)+4:cos(Q:r)+B2sin(I'l:,)+A:cos(Zfl:即+Bzsin(手叩)引耻clcos(孚m。,sin(-警T)+C2cos(tq:丁)(4—20)+D2sin(f2:71)+c:c。s(吾n:n+。:sin(詈n:r)给定系统参数为:m=1,t=1,kc=l,』。=1,fo=1,mo=1,占=1,∥=0.02,kd=O.5,^,=10,No=6,z=44,孝=o.5,E=2,曩=0.01,疋=0.1M2=2,M。=0.01,M:=0.1代入方程(4—20)计算得到主要的系统参数在其它参数固定的时对幅频曲线的影响。(a1)(a2)(b2)第四章弯扭耦合模型分析(c1)(e1)(c2)(e2)53第四章弯扭耦合模型分析(f】)图4.8各参数对幅频曲线的影响由图4.8可以看出各个主要参数对幅频曲线的显著影响。其他参数固定时,振幅随着二阶谐波幅值E、M,增大而增大,随着阻尼系数点,岛的增大而减小,但是在参数一定的情况下,x,的幅值总是比J。的值大。为了验证剪切刚度k。对响应X,的影响,给定k。两个不同的值,t。=o.5,ka=0.8。在不同k。情况下取了3个不同的初压力值,结果如图4.8(c)、(d)所示,如果屯增大,不同正压力下的频响曲线移开较明显,且从图上可以看出亿较大时频响曲线峰值下降。图4.8(c)中,当N。=6时频响曲线的峰值较高。当Ⅳ。=30时频响曲线的峰值反而降低,但是当No=60对峰值回升很多,同时其峰值向右偏移,说明共振时间后移。从这些现象可以得出以下结论:在给定的参数范围内,阻尼端正压力对响应的影响很明显。阻尼端正压力影响阻尼叶片的频响曲线,即影响叶片的频率、响应,并且阻尼端正压力存在一个最佳的值,在此值下阻尼叶片的阻尼最大,减振效果最好。4.4,2通过数值解的比较图4.9Xl和X2的时间历程(Q2=0.6,kd=0.5,No=6,R=M2=2)第四章弯扭耦合模型分析图4.10X1和x2的时间历程(Q2=0.9,.kd=0.5,No=6,E=M2=2)图4.11X1和X2的时间历程(Q2=1.6,ka=0.5,No=6,E=M2=2)X1和X2的时间历程(Q2=1.5,ka=0.5,No=30,E=M2=2)图4.14一趣第四章弯扭耦合模型分析=;蹶‘=;㈣j黼三嗣第四章弯扭耦合模型分析图4.19x1和X2的时间历程(n:=1.05,%=0.5,No=30,E=M2=2,E=M=1)图4・23xl和X2的时间历程(Q2=1.65,kd=0.8,No=80,只=M2=2,曩=Ml=1)通过分析各图,验证了理论解的准确性。具体如下:图4.9~4.11中X1,X2时间历程的幅值与图4.8(a)曲线1中相应的幅值很接近(Ql=0.6,Ql=0.9,Q.s1.6);图4.12~4.14中X1,X2的时间历程的幅值与图4.8(c)曲线2中相应57第四章弯扭耦合模型分析各点的幅值很接近(Q。=0.6,QI=1,Ql=1.5):图4.15~4.17中X1,X2的时间历程的幅值与图4.8(d)曲线3中相应各点的幅值很接近(Q。=0’8,Q。=1.15,Q.=1.65)。图4.18~4.20中X1,X2的时间历程的幅值与图4.8(e)曲线2中相应各点的幅值很接近(n,=o.5,Q。=1.05,Q。=l,65)。图4.21~4.23中x1,x2的时间历程的幅值与图4.8(f)曲线3中相应各点的幅值很接近(Q,=o.65,Q,=1.15,Q.=1.65)。4.5结论本章研究弯曲扭转耦合情况下的汽轮机叶片模型,建立了弯扭耦合方程,并分析两个阻尼器并联时模型的动力特性,摩擦力将阻尼端正压力Ⅳ在叶片振动过程中的变化考虑进去(设为线性变化)。运用谐波平衡法对运动方程进行分析,分别求出一阶谐波Q,=l和二阶谐波Q,=l的稳态解,然后分析阻尼系数,阻尼端正压力初始值,摩擦系数,阻尼器刚度及各阶激励的幅值与频响曲线的关系。结果表明,阻尼器刚度增大会使阻尼叶片的频率升高,并降低阻尼叶片的响应。阻尼器接触面正压力的变化会显著影响阻尼叶片的频率和响应幅值,即频率会随着正压力升高,而对于幅值,则存在一个最优的接触压力,使得幅值取得最小值。通过与数值解的比较知道解析解的准确性。第五章全文总结第五章全文总结本文从叶片阻尼机理和叶片激振力的分析建模,单自由度和弯曲扭转耦合阻尼叶片运动方程的建立,分析等方面,对阻尼叶片非线性动力特性进行了计算分析,得出了很有意义的成果。1.研究了阻尼叶片的阻尼机理,给出了叶片各种阻尼的产生机理及数学模型。着重研究了干摩擦阻尼产生的机理,及其对叶片振动能量消耗的机制。推导出了各种摩擦模型下的摩擦力的表达式,并利用谐波平衡法对摩擦力进行了合理的简化以有利于理论分析。2.多频激振力下单自由度阻尼叶片的动力特性研究。应用平均法求出了叶片运动方程的解析解,并分析了系统参数,主要是阻尼器参数如阻尼端正压力、摩擦面剪切刚度、摩擦系数,以及激振力的大小等,对阻尼叶片的动力特性的影响规律。并通过数值解的比较验证由解析解得出的结论。3.弯曲扭转耦合模型的建立和动力特性分析。考虑弯曲扭转耦合情况下阻尼叶片的力学模型,建立阻尼叶片的运动方程,采用谐波平衡法对其进行分析,得出了阻尼器和激振力参数对叶片动力特性的影响规律。对多频激励下阻尼叶片的动力特性的理论研究,目前做的工作不是很多。本文的研究虽然还只是初步的,但得出了一些对实际应用有意义的结论。阻尼叶片的动力学特性有着很丰富的内容,其中蕴涵着大量的现象需要研究,还有很多领域的工作未能展开,需要做更进一步的研究。1.用三维有限元实体模型模拟阻尼叶片,对叶片和阻尼结构进行更加精确的描述。随着计算机技术的迅速发展和计算机计算速度的提高,这项工作会在不久的将来得以实现。2.对阻尼叶片进行整圈计算,即对整圈叶片进行网络剖分,建立有限元方程,然后求解方程。现在的模型认为叶片的边界条件为周期性边界条件,这样假设不可避免的漏掉某些模态,给计算带来了误差,如果整圈计算就可以避免此问题。整圈计算的难题主要是在计算量及判断解是否收敛上。3.在实验研究上,要用真实的叶片进行研究,最好在动平衡机上做,可以获得叶片的重要参数,得到的结果将对理论分析和叶片设计具有极为重要的价值。参考文献参考文献1.杨海光,汽轮机叶片安全准则的研究,扔微乙穆溯,1983,19(4)2.BatesR.,efa1.SteamTurbineBlades:ConsiderationinDesignandaSurveyofBladeFailures,EPRIResearchReport,19813.Wu—YangTseng.Untwistingofthe31Mod7RetrofitBladeandtheEffectofC0nstraininguntwist,EngineeringReportEM一1314,盼stingHouseE10etricCorporation,19734.Beards.C.F&Williams.J.L,ThedampingofStructuralVibrationbyRotationalSlipinJoints.JournalofSoundandHbration,1977,53(3).333N3405.黄永东译,新26英寸末级叶片的可靠性和经济性,衣方汽轮扔,1994,NO.3~4,3l~436.李小华,汽轮机长叶片振动特性及响应的研究,西安霓透7t.学砑±趁史j19947.村田肇等,日本汽轮机剩余寿命诊断和延长使用寿命的措施,燕力叠萤砬≠蒡俎,1992,NO.6,69~838.李辛毅,季葆华,盂庆集,采用整体围带或翅翼型凸肩提供摩擦阻尼的叶片典型结构,垮花沈茹专褊1996,38(4),217~22l9.E.Watanabe等,张云燕译,高效安全的新型低压末级自带围带叶片,务方圪轮扣,1994,3~4,59~6910.周宏利等,300MW优化机组自等围带预扭叶片设计,汽轮扣友刁e’1992,34(3),1~411.李小华,汽轮机长叶片振动特性及响应的研究,西安交通大学硕士论文,199412.J.T.Oden,E.B.Pires。NolocalandNonlinearFrictionLawsandVariationalPrinciplesforContactProblems;inElasticity,ASME.JournalofAppliedMechanics1983,50(1),67~7613.A.V.Srinivasan,B.N.Cassenti,ANonlinearTheoryofDynjmicSystemswithDryFrietiOnForces.ASME.JournaiofEngineeringforGas7urbinesandPower,1986,108(3),525~53014.C.H.Menq,J.BielakandJ.H.Griffin,TheInfluenceofMicrosliponVibratoryResponse,Partl:ANewMicroslipModel,JournalofSoundandVdbratio&1986,107(2),279~29315.C.H.Menq,P.ChidamparamandJ.H.Griffin,FrictionDampingofTwoDimensionalMotionanditsApplicationinVibrationControl,JournalofSoundandVibration,1991,144(3)。427~44716.C.H.Menq,J.H.Griffin,AComparsionofTransientandSteadyStateFiniteElementAnalysiSoftheForcedResponseofaFrictionallyDampedBeam,ASME.JournalofV2"bratfort,Acoustics,Stres&andReliabilityinDesign,1985,V01.107(1),19~2517.SrinivasanA.V.,CuttsD.G.,DryFrictionDampingMechanismsinEnginesBlade,ASME.JournalofEngineeringforPowel%1983,105(2),332~341参考文献18,任勇生等,带凸肩叶片的稳态响应计算及其阻尼减振分析,振动工程学报,1994,7(2),103~11119.Nayfen,A.H,andMook,D.T.,NonlinearOscillations。JohnWileyandSons,NewYork,197920.C.Pierre,A.A.Ferri,E.H.Dowell,Multi—HarmonicAnalysiSofDryFrictionDampedSystemsUsinganIncrementalHarmonicBalanceMethod,ASME.JournalofAppliedMechanics,1985,52(4),958~96.421.J.H.Wang,W.K.Chen,InvestigationoftheVibrationofaBladewithFrictiondamperbyHBM,ASME.JournalofEngineerJngforGasTurbinesandPower,1993,115(2):296~29922.T.M.Cameron,J.H.Griffin,AnA1ternatingFrequency/TimeDomainMethodforCalculatingtheSteady—stateResponseofNonlinearDynamicsSystems,ASME.JbUrnaIofNonlinearDynamicsSystems,ASME.JournalofAppliedMechanics1989,56(1),149~15423.R.Kielb,J.H.Griffin,C.H.Menq,EvaluationofaTurbineofaTurbineBladeDamperUsinganIntegratedApproach.AIAA88—2400,AIAA/ASME/ASCE/AHS.29thStructures,Structrual‘DynamicsandMaterialsConference,Williamsburg,VA/Aprial18~20,t98824.M.A.Prom。AMethodforCalculatingVibrationFrequencyandStressofaBandedGroupofTurbineBuckets,Trans爿.蜀自EV.80,NO.1,1958,169~18025.A.L.DeakandR.D.Baird.AProcedureforCalculatingthePacketFrequenciesofSteamTurbineExhaustBlades.TransASMESeriesA:JournalofEngineeringforPower,1963,324~33026.J.ThomasandH.J.Belak,FreeVibrationofBladePackets,,舱矾Eng.&五19(1),1977.1327.林俊灿,叶片组振动计算,汽轮机,198028.Wang,J.H.&Shieh,W.L.,1991,Theinfluenceofavariablefrictioncoefficientonthedynamicbehaviorofabladewithafrictiondamper..厂Sound妇反149(1):137~14529.Ferri,A.A.&Heck,B.S,1998.VibrationanalysisofdryfrictiondamperturbinebladesUSingsingularperturbationtheory,TransactionsoftheASME,。,FibAcoustics,120:588-59530.Berthi儿ier,M.,Dupont,C.,Mondal,R.&Barrau,J.J.,1998,BladesforcedresponseanalysiSwithfrictiondamper,TransactionsoftheASME,,gibdcousticS120:468~47431.Toufine,A,Barrau,J.J.&Berthillier,M.,1999,Dynamicstudyofastructurewithflexion—torsioncouplinginthepresenceofdryfriction,Nonlinear。DynamicS18:321~33732.Wang,Y.,1997.Stick—slipmotionoffrictionallYdampedturbineairfoils:afiniteelementintime(FEI)approach,TransactionsoftheASME,,盱66参考文献19:236—24233.季保华等,1997。阻尼连接叶片的试验研究。隋拗Z岌习e,39(4):238—243Acoustics,134.晏水平,黄树红等,2001,汽轮发电机组轴系扭振引起的叶片响应计算,动力r窟21(4):1288~129135.鲍文等,2000,300MW汽轮机低压缸末级长叶片与轴系扭振的耦合振动,汽轮扔貅,42(2).66~6836.N.F.RiegerandR.Subbiah,TechnologyManualforSteamTurbineBladeAnalysiS,Stress殆chnologyIncorporated,Roxhester,NY,March21,198637.RajagopalSubbiah,BladamComputerCodefortheEvaluationofDampinginTurbineBladeComponents。Report‘i.STI.198738.戴沛德,阻尼技术的工程应用,北京清华大学出版社,199139.J.S.Rao,TurbineBladeVibration,JohnWilon&SUN,199140.Z.Kovats.DynamicCascadeFacilityandMethodsforInvestigatingF10W—ExcitedVibrationandAerodynamicDampingofModelLowPressureBladeGroups,ASMEpaper,79一wA/GT一641.E.F.Crawley,AerodynamicDampingMeasureJl】entsinaTransonicCompressor,ASME.JournalofEngineeringforPower,1983,108(3),575~58442.邵荷生,曲敬信,许小棣等,摩擦与磨损,北京煤炭工业出版社,199243.孙家枢,金属的磨损,北京冶金工业出版社,199244.沈心敏,摩擦学基础,北京航空航天大学出版社,199145.陆大雄,摩擦学导论,北京出版社,1990.46.茹铮,塑性加工摩擦学,北京科学出版社,199247.金永昕,工程摩擦学,杭州浙江大学出版社,1994掘J.S.Rao,N.S.Vyas,DeterminationofBladeStressesUnderConstantSpeedandTransientConditionswithNonlinearDamping,ASME.JournalofEngineeringforgasTurbinesandPowe]-,1996,l18(2),424~43349.J.T.Oden,E.B.Pires,NolocalandNonlinearFrictionLawsandVariationalPrinciplesforContactProblemsinElasticity,ASME.JournalofAppliedMechanics1983,50(1),67~7650.LF.Wagner,J.H.Griffin,B1adeVibrationwithNonlinearTipConstraint:ModelDevelopment,ASME.JournalofTurbomachinery,1990.儿2(4),778~78551.朱因远,周纪卿,非线性振动和运动稳定性,西安交通大学出版社,1992,32~3452.C.一H.Menq,J.H.Griffin,J.Bielak,TheInfluenceofaVariableNormalLoadontheForcedVibrationofaFrictionallyDampedStructure,ASMEJournalofEngineeringforGasTurbinesand总weT,1986,Vol108(2)。300’30553.李锋,季葆华,谢浩,孟庆集,汽轮机叶片激振因子影响因素的分析研究,汽努祝疫蔫1999.10,V01.41No.554.季葆华等,汽轮机叶片干摩擦阻尼机理及其数学描述汽轮机技术1997volI39参考文献no.6345—34955.Berthi儿ier,M.,Dupont,C.,Mondal,R.&Barrau,J.J.,1998.Bladswithforcedresponseanalysisfrictiondamper.TransactionsOftheASME,dVdbdcoustic&120:468~47456.季葆华等,阻尼连接叶片的实验研究汽轮机技术,199757.Ding,Q.,Leung,A.Y.T.&Cooper,J.E.,2001.Dynamicexcitedhystereticsystem.,Soundv01.39nO.4238—243ofaanalysiSself—rib,245(1):151—164V01.18no.158.陈予恕.非线性振动.天津:天津科技出版社,198359.刘东远,孟庆集,汽轮机叶片激振力分析及优化动力工程199812—1860.季葆华等,现代汽轮机叶片阻尼结构形式设计方法的研究西安交通大学学报1999v01.33no.840一4361.M.A.Prohl,AMethodforCalculatingVibrationFrequencyandStressofaBandedGroupofTurbineBuckets,TransASME,V.80,NO.1,1958,169~18062.A.L.DeakandR.D.Baird.AProcedureforCalculatingthePacketFrequenciesofSteamTurbineExhaustBlades,TransASMESeriesd.'JournalofEngineeringforPowe91963,324~33063.J.Thomas13andH.J.Belak,FreeVibrationofBladePackets,d.Mech.Eng.Sci,19(1),1977,64.林俊灿,叶片组振动计算,苊话祝,198065.ALAINTOUFINE,JEANStudyofaJACQUESBARRAU,MARCBERTHILLIER,DynamicPresenceofDryStructurewithFlexion—TorsionCouplingintheFriction,NonlinearDynamics.18:321~337,1999作者读硕士期间发表的论文及科研・情况发表的论文1.孙艳红,丁干,叶敏,两自由度复合材料层合板的非线性参激振动,工程力学增刊,V01.20,20032.SunYanhong,DingQian,YeMin,ApproximateSolutionofaTwodegreeoffreedomNonlinearsystemwithParametricExcitation,Thel】thWorldCongressInMechanismandMachingScience,2004.43.丁千,孙艳红,汽轮机干摩擦阻尼叶片多谐波激振响应分析,振动与冲击(未发表)科研情况干摩擦阻尼器透平叶片系统的黏滑振动动力学研究,国家自然科学基金项目(10272078);堕茎——附录I单自由度模型分析1.第一阶谐波的主共振Q.。1E=-;(28v2a2y2石一4a2(fl—a)一sin(,9)yz,9-4a2颐sin(,9)yz2十16ny3掌Q1妒口+38aflay3石一lOaa2Y3石一222y3f2flsin(2'9)一2sap2y37/"+2y4fleal7/"2+y4c嘻Ql疗2—38afla2J,2万+lOa2口2y2石+32a4flair+16a4口219+16a3口2ny+16a3∥2y(z一.9)+16a4卢219+24a2cr‘;nlY2zO一16a3口20y一24a2/臼芎Q1_y2z,9+24a2脾1Y2牙2—16a4∥27/"一32a4膨J+32a3flay(,9一/t")一24a2(姻ly2厅2—16a4at'2万一16azct,9fy3Q1—16z2哟301+1622Y3{nl口口+71"2Y3甄Qlcos(,9)+222y3f2asin(2,9)+4a2觑ylr2sin(,9))/(z2QlY3)Z2=(-24Ay2V3a,9sin(,9)+16f2y2%膨2cos(2'9)一6AY2a,gZ2,cos(2y3一'9)+24fly2V3azsin(,9)一48fry2zfl,gsin(,9)一16f,y2∥3aZ2cos(2,9)+4f2y2aZ2sin(2缈3一19)一6AY2flr,z2cos(2v3一。9)+6frY2aZ2zcos(2v3一口)一24fly2石2asin(,9)十48Ay2妒3(fl—a)cos(Za)+6ZY2a'gcos(2V3一19)-3fly2(口一声)Z2sin(2p,3一19)一6flY2fl,gcos(2妒3一19)+6Zy2∥理2cos(2妒3+t9)+72y3∥2sin(2v))一16f2y2矿3(fl—a)cos(2t9)一144y3y3卢2+192y3∥2厅+84y3球2万一3』y2aZ2sin(2v3+口)一6Zy2zaZ2cos(2弘/3十19)-6fKy2,gaeos(2v3+∞+12Ly2妒3fleas(,9)+3工y2胆2sin(27t3+回+6AY2fl,gcos(2V3+|9)+6frY2∥。舅z2cos(2V3一圆+24fly2c趼2z2sin(t9)+72flay3Z2sin(2y3)一6AY2∥韶2cos(2v3+回一12/;y2膨2妒3cos(,9)+6flY2a"922cos(2{u3+口)+48fzy2azcos(Z'9)+4f2ay222sin(2v3+2圆+24Ay2(口一f1)az2y3sin(O)+24fly2(fl—a)Z2y3sin(,9)z十48^y2flz矗Z2%sin(a)一120f:ayasin(2圆+128f2a2z2(fl—a)sin(2'9)+96fla2a,9cos(O)+96fla2,&rcos(O)+192fzZa2(卢一口)sin(Z'9)一40f2y2xacos(2,9)一96fla2卢19cos(,9)+48fta(∥一口)ysin(圆一48flz2a2(a--f1)sin(,9)一12fly2axcos(8)+24f2y。8冗cos(2,9)一96ftazazeos(a)+128f2附录口2(口一f1)sin(23)一48fIa2(fl—a)sin(3)+40f2y2azZ2cos(2a)+48fly2卢疵2cos(,9)+48flZ2aaysin(3)十192fzZaaysin(Z19)+24甑y2刀QIsin(G)+96flayazcos(8)+96fla222azcos(3)+12fjy222azcos(a)+12ffy222ay3cos(—9)~4f2Y222∥sin(2矽"3—23)+6ZY2fln'cos(2妒3一圆一6fly2p∥cos(2≯'3十19)一12fry2y3口cos(3)+24fly2y3flasin(oa)一96f1..vZ2a(a一卢)口cos(8)+96flZ2a2pScos(3)一96^aayacos(8)+128f2ayflsin(2,9)一24fly2y3,8xsin(O)一4厶y2z2,8sin(2y3+23)+144y3妒3c妒一96y3xa,8+48f22y3手2疗一84y322口2丌一192y3z2∥2玎+144y3y,z2卢2+48y’口研+96y3I珈一6ZY2azeos(2y3一口)+6Zy2azcos(2矿3+t9)一144y’yaaflZ2—48y3,zQ;毒222+96y322,aax+4Ay2flsin(2¥q+2国+4f2y2∥sin(2P3—219)一72y3,aasin(2矿3)一125Y2(fl—a)sin(Z8+2y3)一3f,y2(卢一a)sin(3+2妒3)一3石y2(fl一口)sin(2p"3一|9)一12厶Y2(∥一a)sin(-Z3+2y3)~72y3∥222sin(2妒3)一48y3ao-彪72—9@3Zb咒z2一幔y2∞in(2矿3—23)--4f2y2asin(2_}uj+2卵一96ftZ2ay(fl—a)3cos(3)一96flay(z一圆Boos(a)一96fZ2口2aacos(3)一24舅ty。Z2刀Qlsin(3)一128Lz2ay(fl—a)sin(23)一192f2,Zayflsin(Z3)一48fjZ2aflvsin(,9)一24f2y222巧口cos(2¥一96fla222口厅cos(,9)一48fly2pncos(,9))/(96ny3Q。(Z2—1))2.第二阶谱波的主共振Q:“1巧2雨丽1【10080aayfla,9+7560J(fl一口)却zy2矿+5040as口2归.-10080a3aflyz+5040a4∥2口+10080a4fla(,r—19)一504033口2ya+5040口3fl2y(z一.9)一7560a2(fl一口)国2y2棚+5040a4口2(19一万)~5040a4f12刀+1260(fl—a)a2f2yz2sin(a)+1260(a—f1)a2f2yz3sin(a)一700flY2艺石(cos咕回+c。s毫.9))+2800‘y3a筇c。s哇t9)--1400‘y3(fl一口)妒3cos(1口)+7792‘y3肛os(吉19)-1536fly3azcos(I,.0)-8820掣3f12z-504‰2y3如2(fl刊+4096ny3‘抑2sin哇趴31522y3璺")2,2c。s(即3150口2y2a2万一350flY3(仆f1)sin(2y3+吾.9)一350fly3(a-f1)sin(2y3一吾.9)+11970伽3石附录+630争4脾2石2—11970a2肋2石一3150口口2y3z_5040y3.z锄2a(a一卢).9+8820a2y2f12z+315y4勰2石2】z2=(11520f1222y3+4320(fl一口)∥纠3sin(2妒3)一360f2y2肛2c。s(2y3+19)+360f2y2(fl一口)z2cos(2v/3—8)4-960fly2(a|9+肋)sin(2缈3一丢占)+960fly2fl(石一跏咄2y3+吉趴360f2y2彻2cos(293+uo)+2880fZy2嘶2c。s‘i1zo)-2880f12y32222+10240fltrny2sin(尹1—1080f12y3尼2sin(2矿3)-1440f2y2簖3sm(#)+lSOf2y2筋sin(293+8)+720fzy2鲫sin(2∥3一丢z圆+96。‘y2a,gsin(2%+丢t9)一3840‘y2册r2cos(161)+720,Z),2钎sin(2%+吾z回+180f2(口一f1)y2zsin(2|;c/3-.9)一180如fly2zsin(2%+占)一720fzY2肛sin(2妒3+昙zt9)一720fzy2flzsin(2y3一i11z回一960‘y2嬲siil(2吩一吉占)一?6。五y2嬲sin(2%+吉t9)I1一一96%y‘flOsin(2≯'3一吉J)+720f2aYaZZ2sin(O)一5760f2a2∞一卢)石2c。s(回一7680f1(口19+励一flS)a2sin0口)一720f2Y2伽2cos(8)一2880f2y2fln:2cos(8)一2880f2a2fl石sin(O)一480‘易y心n(寻t9)+1920石正y石s协((i+z)主圆一144。‘易y石sin哇19)一3584石y2卵sin(三0—1920ffz_y万sin((z一1)妻.9)+90f2y222a19万cos(2弘,3+|9)一960zy222芦|9妒3cosel9)+1920fy222fl&rcos(寺8)+1440f2y2石(尸一口)口y3sin(印+360f2yrZ2(fl-a)y3石2sin(#)+90≯2y222(∥一a),gzcos(2v/3一口)+180Ly222a93zcos(,9)-960fy222(cm7一肛一口t9)y3石cos(三19)+720Ly222fl,gz2sin(a)一90f2y222flSzcos(2妙3+t9)一180厶多222%fin"cos(8)+1440f2yZ2za6lacos(#)+1440f222a27/-2acos(3)一7680fay(口一∥)厅siil(÷t9)+5760f2aypz8cos69)一1440f222口2az8cos(8)一360五附录y222万2{n2a2ax'9cos(,9)-?2880f2ayan-sin(回+5760f2ay钎2cos(,9)+360ff2yzrrZ2sin(O)+720Ly。z2巧2flcos(8)一720正Ⅱ222嬲sin(O)+5760Asin(吉'9)+2880f2a2伽sin(,9)+1920flZ2aya,rsin(去t9)一1440f222ayfl,9石cos(,9)+5760fzZa2/勋sill(寺z19)+1920zz2口(yfl*9一口钎)sinG圆+7680石a2a//-sin(寺口)一5760f2a2flsr*9cos(占)+1920/;Z2口2Ja厅sm(去回一7680fay,9(fl一口)siI《妻口)一1920fZ2a(ya,9一aflTr)sin(毒,9)+180f222y2蝴2cos(,9)+2240fy222孝衲2∞sG圆一1440f2a222厅2∥cos(¥一576%缈溉9cos(|9)一1920fZ2ayfllrsin(寺回~5760fzZa2凹sin(去zt9)+13184fy2嬲sin(寺印一1920fZ2a2(fl~口)|9sin(砉口)十1440f2ayZ2仃2(fl—d)cos('9)+144‘0f2a222巧日19cos(19)一5760f2锄2cos(,9)+1440fzy2、石2考Q2sin(回+896zy222,f奶sin(去口)+2880Layfl万sin(t9)一8960fy2勿Q2cos(寺19)+720f2a222励’sin(拶)+120一^y尼2sin∈固一5760fzZay(fl—a)石sin(寺zt卵~240fly222甜(19—7r)sin(2妒3+去口)删协塞枷p矛铲H船回回。1萨够谢所勋㈨砒弘抱町广弘,3』玩cos(寺,9)+720f2y2If,3触cos(|9)一240fy222∥(疗一口)sin(2驴,3+去J)一90Ly2zr(p-a)疗2eos(2y,一,9)+360f2y222觚3sin('9)+1080y3z‰flsin(2∥3)一2560Ay2221fo-sin({圆+240fy222(fl一口)口sin(2妒3一妻¥一45fzy222(口一∥)万sin(2弘勺+圆一720五J,2妒'3axcos('9)一45f,y222五一f1)xsin(2妒3一O)+3840fy2妒3觚COs(吉,9)+960fly222石2口cos(三口)一360rgy2(a一∥)石口cos(2y3+*9)-2880f2y2分≥2*gsin(O)一90rgy222(a上-f1)x2cos(2Iv3+口)一1440rgy2y,(fl-a)n2sin(O)+360f2y2弘,3厅19(口一f1)z2sin(*9)+3840Ay2虬够一口)挣cos暗印+1920石兀J,庇sin(寺Z'9一去回+1920ffzyxZgin(128+妻t9)一720f222ay∥:rsin(口)一3296fy222嘶sin(寺.9)一8640y3∥3卢2疗+57‘60(盯一磊)/曰靠2y3+2160y322卢2V-'3万一1440(a一口)肛2y322—720atrz2y322iZli霎≥髭526帆45茅y83掣V1"261触0x手Y专;芝?瑟527y6、0zj8Y;警?1弩44砸0石;移Z2呸)y32毒2壶j+石2口2—2322口2)“一23Q2+2y32Q2)附录II弯扭耦合模型分析1.考虑Q.:1的系数。A.=(砰N3Fl+N;P4Ml一砰F:t—E∞?∞:+Ⅳ]E珊:+N2N1M109t∞;一N:鼻只E甜;+Ⅳ。B们l国?一N。BN3M。一N4只只鼻+Ⅳ一只E珊;+2P3FloJ?tV;一日2E∞?∞;+只2国22Ⅳ3E-2P3N3F,co;一N4Mlco,co:+Ⅳ4Ⅳ3M103;+N;P4M1国;一N4P,P2F:01092+Ⅳ4Ⅳl鼻肘’l国l国2一N2e2F,eoltO;一N2N3e,MI国;+Ⅳ2鼻M1m12叩22+Ⅳ2P2只E∞l国2一N:N。B^f。Ca).珊:),(P32Nto,;一2P3国?∞;一2Ⅳ;P3∞;+Ⅳ;巧∞?一2碍Ⅳ3∞?+Ⅳ;鼻2国;+Ⅳ?国?国;一2N30)?oJ:+Ⅳ;P4%o;+E2∞?甜;一2N4N,P2P3co?+2NIN,P2r-otco;十P;m:+Nj∞:+P;Nj+2N4N1只置∞1卯2+oJ:oJ;+Ⅳ;马2国;∞;+Ⅳ?鼻2国?∞;一2Ⅳ?只卯?珊;+2N。只B∞?一2N4N3只B+2N4N3只国;L2N4只珊?出;一27,2N3co,co;+4N3只珊?国;+2N2B国?脚;+砰N;+2N2P4鼻∞?吐);一2N4N3B最脚l∞2+2NdB鼻∞2山?+2NlN2只脚l山i一2NlN2只只国1032—2N2BB国2国?+2N3N2只P2∞l卯2-2N1N2最鼻国?甜;一2N3N2只珊l国;一2N3N2只只m;)附录马=(Ⅳl一日O)t+只Ⅳ2Ⅳ3Mt脚:一只与Ⅳ4巧印・一只Ⅳ4ⅣlM,q一只只Ⅳ:鼻棚:Cl=一(Ⅳ;MI∞;+Ml卯?国;一埘M1只一Ⅳ4曩牟一Ⅳ1只只鼻国l彩2+N3丘鼻曩国l国2一2P3NlF,c01∞;一BⅣ2M】∞?国2一Ⅳ3Ⅳ4P,Ml国2一Ⅳ3.v2M1031;+B^彳lⅣ:q+Ⅳ4鼻只F,c02十Ⅳ4EMl珊知2+N4e2F,col珊;+Ⅳ4ⅣlMICOl埘≯N2&FlCO;+Nl曩∞1CO:+Ⅳ2Ml国?∞i—P,e2N2FI国1CO:一N1N2P,Ml(70l国;+Ⅳ;只^彳l∞l∞;+Ⅳ1鼻2F,col珊;)/(坪N÷co?一2B彩?∞;一2—V;Em;+Ⅳ;碍缈?一282N3co?+置2醒《+研国12∞24—2N3∞?耐+孵砰∞;+最2国;研一2马Ⅳ。N。B砰+2N4N1P2∞?功;+砰∞?+Ⅳ;趔+只2N;+2N4NI鼻P4∞l国2+∞?∞;+Ⅳ;碍脚?国;+Ⅳ?鼻2国12Ⅲ22—2N?B珊?∞;+2N4E只国?一2N4N3B只+2N,IN3只印;-2N4只珊?出;-2Ⅳ382珊12∞22+4N3B国?《+2N2B国?曲;+巧Ⅳ:+2N2只只研国;一2Ⅳ4N3鼻只q国2+2N4只B国?国2+2N1只Ⅳ2%3∞l一2N1N2B只珊l印2—2_Ⅳ2BB埘?∞2+2N2N3只县国l珊2—2NiN2B最出?印;一2Ⅳ3BⅣ2∞;国I一2N2N3弓只国;)一只最曩彩?国:一N。BE∞?国;+Ⅳ1尸2只只脚?一N4只^f1茸+.Ⅳ4.Ⅳ1只M。珊?+N掸t&M1∞p2+N2&Mt∞i国l—N拶;&Mt0.1pl一冀Mt出:一2N}ML∞;∞12+只曩国?∞;一Ⅳ3只E国;+N?M,国?脚;一N4E印?口一只与只砰+N,Ⅳ。只^4+N,只只互+2N3BM。∞?一Ⅳ?马M。∞;),(∞—fN?一2B国?∞;一2Ⅳ;只国;+Ⅳ;碍国?一2口Ⅳ3珊h埘鼻2彩;+Ⅳ12珊?国;一2N3奸国:十Ⅳ;日珊;+E2硝嗣一2玑Ni只B国卜2N4N1只印扣;+耳研+孵∞:+口埘+2N4NI鼻只国1092+硝心4+Ⅳ;碍砰国;+Ⅳ÷B2∞?0-)22—2研B砰脚;+2N4只B砰一2Ⅳ3Ⅳ4只0+2Ⅳ3N4只国;一2Ⅳj只国?珊;一2日2N3珊?∞;+4N3B国?珊;+2N2马脚?卯;+砰川+2Ⅳ2只只珊知;一2N3N4e,&col珊2+2Ⅳ4鼻马国?∞2+2&UlⅣ2棚;q一2NlN2只只脚l国2—2Ⅳ2B只国?国2+2Ⅳ3N2只B国l∞2—2N1N2E最国12删22—2Ⅳ3N2只∞;珊I一2N3N2只置印;)附录D。=一(一N3只BE国。一鼻M。∞?国:一只互国?卯;一E只鼻珊?脚:+B只E∞?一Ⅳ;互^,。∞:一Ⅳ2曩只2040,2+N2N3只Ml印2一N2e,Ml脚扣2+NIN2只Ml脚知2一只Ⅳ2砰国2Az=一(4碍E∞?一只Ⅳ:哪一碍EⅣ3一16EⅣ3∞:一4N,N,M:CO;+Ⅳ4Ⅳ3B时。+N1F,e,e3col—N1E只∞;国l+NIEE最研国2一Ⅳ12EM研卯2一NIⅣ4只吖lO)1+2N3MtPKo(c02+N]只B脚;由I+Ⅳ3E鼻只印2+Ⅳ3|v4M1最珊l—N4嵋只研)“只2研砰一2与国拟2-2P3N;呸2+Ⅳ;只2研一2坪N,CO?+墨2孵西+研∞?脚:一2N3co,co:+砰Ⅳ;国;+#2研∞;一2NIN4BB甜?+2Ⅳ1Ⅳ4B∞12Ⅲ22+碍司+埘国;+Ⅳ;碍+2N1N4PlP4(o。曲:+翻;∞;+孵碍卯?国;+Ⅳ?点2国。2脚;一2研B研珊;+2N4B只国?一2Ⅳ4N3只只+2Ⅳ4N3只∞;一2N。只∞。2w22—2N3P12co,co;+4N3P3Co?60;+2N2马珊?棚i+|Ⅳ:碍+2N2只E国?棚;一2N3N4P2Ptcal03.2+2N4P:,ogo,2+2N1N2只国;∞l一2B只%Ⅳi彩】∞2—2N2BB国?国2+2N3N2只只曲l珊2—2NlN2只鼻奸∞;一2Ⅳ3N2最∞23∞I一2Ⅳ3N2鼻只脚;)+4只ENaE国22+4鼻Jv2N3M2国22+16FxP2N2∞lm;+4EN4足E∞l国2—4疋Ⅳ2P2B∞l脚2—16PlN2M2国12口22+64‘彩;脚:一4置Ⅳ4NIM20)1(02+4N2NlM2只彩1092—16N2NlM203l∞i+只B—v4E一4只Ⅳ4B∞;一4E2Ⅳ3E司一47,M2孵国;+8BⅣ3E国;一4P3N4M2研+16E2E研∞22—32BE印?国;+16N4M:国?国;),(4口Ⅳ?∞?一128B∞;珊?一8.v;E彩;+4NjPi∞;-8P;N3∞j+4P?Nj畦+64N;∞12∞14-128N3∞;∞:+4孵碍脚;+64P。2q4脚;-SN4Ⅳ,最B珊?+32虬Ⅳ。P,co(co;+16只2国?:硝+16Ⅳ;∞;+呼Ⅳ;+8N4NtE只国l国2+256w+16Ⅳ;碍国?甜;+16N:P12研∞;一32吖Bq2Ⅲ22+8N4B只国?一2Ⅳ4虬只只+8Ⅳ4N3只∞;一32N4只国12Ⅲ22—32N3P12珊?脚;+64N3P3co'co;+128N2B脚?卯;+.V;砰+32Ⅳ2弓只og呸2—8N4N3E最国lO)2+32N4E只国?国2+32Na—Ⅳ2只∞l∞;-8N2NlB只∞l国2—32N2足Bm?国2+8N2Ⅳ3巴B国l珊2—32PtP2N2Nl山12∞22—32BⅣ2N3ml国;一8N2虬#只国;)71附录B2=-2(4N2N3M2m;一Ⅳ2Ⅳ3只^,2珊2+E.v4Ⅳ3M2国2—46M2_v22Ⅲl“22—4Ⅳ2只E国;一16N2M2研国;+4弓Ⅳ2NiM2国1国;+4眉BⅣ2E珊l珊;+4BⅣ2M2国2国?+Ⅳ2B只E脚2—16N1E山l∞;一4e,2NIECOl∞;+8BⅣlEq国22—4只Ⅳ4E∞l∞22—4N4N1M20./l∞;一鼻只N4E国2—4弓Ⅳ4M2国2曲?一学NiE∞】+只与Ⅳ4ECOI+只Ⅳ4NlM2珊l一只朋’2Ⅳ42印】)/(4警Ⅳ12∞12-128B国?国;一8只孵国:2+4孵芹∞卜8只2Ⅳ3珊卜4Jv;舅2珊;+64N12司珊24—128N309国24+4Ⅳ;&20j;+64P,2脚14Ⅲ22—8N4N1只B曰+32虬N1&ogo);+16P32∞?+16N;co;+碍N;+8N4Ni鼻只脚l∞2+256∞?国;+16只2Ⅳ22∞t2印;+16P,2Ⅳ?印?脚;-32P:I研研曲;+8Ⅳ4只只国?一2Jv4Ⅳ,B只+8Ⅳ4N3只出;-32P4N4脚?脚;-32P】2N3珊?国;+64N3只彩?∞;+】28Ⅳ2只辞国;+砰枳+32N:墨只脚?国;一8N。N,#最q国2+32N。#县奶研+32NIN2只脚1脚;一8N2N1只只ml珊2—32N2BB国2国?+8N2N3只只∞l国2—32P1&N2Ⅳl研03;一32MN3B国l∞;一8N2N3只只珊;)C2=一(一Ⅳ;BM2+4N4M2N1B砰一16B五Ⅳl碍国;+4BENl只og+4最‘量.Ⅳ3国l∞2—4只最鼻Ⅳl甜l珊2+4N2只肘2NI(010)2—16最五E∞?缈2+M忍只五一4Ⅳ3只&国;+弼Ⅳ4只Ⅳ:+8N3£嬲:群+16Nt2M:彳呸2—4N?M2只印?-32N3M2印?脚;f16只E脚?∞;一4B只Eog一4N4P4M2∞?一4碍N4E国?一16只M2国?一Ⅳ4E砰+64M:国知;+4.Ⅳ;M:国:2—4N2t'2M2N3∞l卯2+16N2P2^f2回国2),(4口Ⅳ?国?一128Bc01'co;一8Ⅳ;B国;+4碍哪回一8碍虬∞卜4孵E20,1;+64N,2砰国;一128N30)?叻4+4Ⅳ;砰国;+64P,104∞22一毫Pj囊N4N砝+32P2N4Nt西武+166W+16Nj∞24七霹N;+8N4NI&只ro,c02+256∞?印;+16碍^r22埘12国;+16舅2Ⅳ?曲12c[,22—32与Ⅳ?国?国;+8Ⅳ4只B珊?-2N3N4只只+8Ⅳ3N4只国;一32Ⅳ4只国?∞;一32E2.Ⅳ3印?∞;+64ⅣjB∞?曲;+128N2P2霹∞i+巧N;+32N2只弓∞?国:一8N,JjvjB置q国2+32N4£E讲国2+32N2N1只珊l脚;一8N2N1只只∞1092—32N2尸2曩∞?∞21+8BBN2N3091∞2—32N2N1忍‘印扣;一32只Ⅳ2N30)J彩;一8Ⅳ2Ⅳ3只点西)D2=2(与最EN3卯l一只只五N103l一4BBE珊?+elM2Ⅳ;脚2一只^彳2N2N3珊2、一只只EN3山2一马M2N4N3国l一8EM2N3研国2—4t'2EN3091国22+EN2爿0)2+4只只Eog珊2+只M2NIN4印l+4只肘2N2国?埘2+4只最NIEOl∞22+16最‘国?国;A3=(N4NIEM322∞l出2一N4尸2鼻只z2国1脚2+224只E∞12∞22一N4M;z4砰∞;+16鼻M2珊?脚2—4&M2N2N1钾曲2—4E昱Eiv,珊?∞2+4&M2N4∞?+4只M2Ⅳ?研国2+4EⅣ2碍砰国2)/(4/'32析出卜128B国?脚;一8B埘国22+4明只2∞;一8球N3∞?+4Ⅳ;鼻2甜;+64N,2q2Ⅲ24-128N3砰∞;+4孵砰司+64p12脚?脚;一8NlN4B只甜?+32NlN4B国?缈:2+16P]缈?+16Ⅳ;珊;+Ⅳ;只2+SNlN4只P4国l国2+256ro(co:+16N;日国?国;+16Ⅳ?#2∞?国;一32BⅣ?研国;+8N4B只曲?一2Ⅳ3N4B只+8Ⅳ3N4只∞;一32N4只国?∞;一32N3E2国12口22+64N3与ogo);+128N2B∞?∞;+Ⅳ;巧+32N2墨只∞?印;一8N4N3最弓∞l∞2+32B日N4卯?国2+32N2N1只国1国;一8只BⅣ2Nl脚l棚2—32N2最B国?∞2+8最只Ⅳ2N3脚l脚2—32Ⅳ2Nl县舅研∞;一32BN2N3国103i一8N2N3EN4国;)一只2Ez4珊12∞22+Ⅳ4N3M322脚;+异2N3Ez2珊;一2B—v3巧z2珊;+P4Ⅳ4F322国:+BⅣ4M322讲+只2Ⅳ,E+Ⅳ;P4^厶+N1N2M3240.1'1&E3一N2BEz4甜I∞i十Ⅳ2只A厶z4砰∞;一N3N2E^厶Z2∞;一N2E只只z2印;+Ⅳ2BB巧z2国;一NIN2只M322国1国2+只^如Ⅳ;z203;一N3N4BM3一N4只只E一碍Ez2国?一只z6∞?∞;+Ⅳ,Ez4脚;)/(鼻2Ⅳ;z2fO;一2只Ⅳ;z2晚E2+Ⅳ?z6珊12Ⅲ24++丑226回国;+Ⅳ?日z2∞?+明只222∞?一2Ⅳ3P]Z2山?一2N326研国;+Ⅳ;&222钙2—2Bz6硝∞;+碍Z4珊?+Z8脚14删24+Ⅳ;z4∞24—2N3N2最z4031国;+2NIN4鼻只z2∞l国2—2Ⅳ3N4鼻只z2珊I珊2+E2N?z4珊?∞;一2只N?Z4国?脚;+4N3Bz4国?∞;+2Ⅳ3Ⅳ4只z2m:一2只2N324珊?国;一2NiN4BBz209+2N4B只z2∞?+2N4E昱z4甜?国2+2N1N2只z4卯iq+2N2置只z4国?∞;一2NlN2EBz4国?∞;一2N3N2E只z2国;一2N2只只z4国?∞2+2Ⅳ2N3尸2只z2COl棚2+2ⅣlⅣ4最z4回遥+孵巧z4珊?∞;+2Ⅳ2足z6脚汹;一2N,只z4研国22—2NIN2B只z2珊l珊2+口埘一2N3N4P}PL+PjNj、附录岛=一(Ⅳ2虬M32203;一Ⅳ2Ⅳ3M3B珊2+Ⅳ4Ⅳ3M3E∞2一N2M,Z4脚细;C3=一(一孵马M3+Ⅳ1N2hM3z一Nl玛z4脚l∞;一Ⅳ1w3马z203l(11;一N2‘只z2。i+Ⅳ2N1M3鼻z2∞l珊;+N2EE最z20Jto);+N2M3只z2∞2∞?一NiE置222珊1研一Ⅳ4NjM3220)I国22+2N1E与z2COl09;一N4E最z2国l∞;一N4M3异z2国2珊?+Ⅳ2EB只国2一Nt曩互只出2+Ⅳ4ⅣI^厶Bq+Ⅳ4巧8弓q—M3最Ⅳ;国l—NI墨碍q)z/(毋2Ⅳ;z2脚;一2只Jv;z2脚;+N?z6∞?珊:+只226∞?国;+学Ⅳ?z2∞?+霹Njz’武一2《N322武一2N326峨∞:+《N;z2∞:一2鼍z6∞:∞22+P]Z4珊?+z8珊?国;+Ⅳ;z4曲;一2N2N3最z4EOI∞;+2皇只ⅣlNIZ2∞l脚2—2E最Ⅳ3眠Z2(I)10)2+只2N]22409∞;一2N?只z4奸《+4PjN324珊?蛾2+2P4N3N422∞;一2最2N324∞;∞;一2马马ⅣlN422m?+2只P4Ⅳ4z2∞i+2墨BⅣ424∞2珊?+2只ⅣlN224珊l珊;+2只只Ⅳ224研∞;一2只BⅣlⅣ224∞?m;一2鼻只Ⅳ2Ⅳ322谚一2BBⅣ224∞2国?+2只B|Ⅳ3N2Z2(010)2+2P2N1N,Z4砰∞;+Ⅳ;只224∞?国;+2罡Ⅳ226D?甜;一2只Ⅳ424∞?国;一2N1N2只只z203'10)2、+N;霹一2墨P^N3N4+NjP:、2∞10.12+N2P2M324(-013n,2一N3N2P2M3Z2∞l∞2一N1置只Ez2∞l国2+Ⅳ3只最Ez2卯I珊2—2Ⅳ3鸩z4珊?珊;一N1马Ez4珊?∞;+Ⅳ?M,Z4∞?∞;一Ⅳ3只Ez2国;一EP2Ez4珊?国2一N4只^毛Z2(11÷+只巧z4∞?∞;一B只Ez2国?一N4碍只z2CO?+Ⅳ。Ⅳ4卫^厶Z2山?一Ⅳ?B^厶z2出?+2N3只M322国?+Ⅳ1足BEz2EO?一N。碍E+N3B只E+Ⅳ;^厶Z202,;+M,Z6∞?曲;一只^如Z4037+Ⅳ,Ⅳ。只』l如),(E2Ⅳ;z2国:一2只|Ⅳ;z2印;+N?Z6国,2删:4+只226∞?∞:+碍N?z2回+碍.v;z2砰一2Ⅳ3只222∞?一2Ⅳ326∞?国;+坪孵Z2国;一2Bz6司∞;+口Z4珊?+Z8印14棚24+孵z4国:一2N2N3P224∞;q+2NIⅣ4只只z2珊10.12—2.v3Ⅳ4鼻马z20)10)2+鼻2Ⅳ?z4∞12鲫22—2只Ⅳ?Z4脚?国;+4Ⅳ3与z4090);+2Ⅳ3Ⅳ4只z2国;一2只2N,Z4∞?国;一2NiN4B只z2国?+2玛只Ⅳ422印12+2Ⅳ4B鼻z4国?∞2+2NIN2只z4珊;国l+2N2只只z4∞?nE2—2NIN2E最z4斫彩;-2N2N3只只z2∞;一2N2B只z4国?吐+2Ⅳ2Ⅳ3最Bz2q吐+2NIⅣ4最z4奸埘;+碍孵Z4q2甜22+2Ⅳ2县z6研国;一2N4尸424国?珊;-2N1N2B只z2珊l脚2+口N;一2N4N,只B+碍Ⅳ;)74.:一—堕重——————————————一D,:z(一N2N3P4M3棚2一N2N1只M322co,c02+N2F3巧∞2+Ⅳ2E碍Z‘国扣2+N2只^f322co',c02一NlB只E国l+N4N1e,M3国1+ⅣlP4Ez2国;∞I一2N,EM。Z2国?国2+N(PlM322co,c02一Ⅳ3置只羁吐一N3P2Ez2够;吼+N3P2P3巧∞】一N4NsP2M3∞1-N1EBEz2研珊2+N4P2MsZ2∞?十丑只Ez2研∞2一只只只z2研+最‘z4四∞;+PIM324科0)2+Ⅳ;鼻^彳3∞2)/(最2Ⅳ;z2∞;一2P#;z2脚;十N?Z6国?珊;+P1226卅国;+z:霹N;武+z2琏N2s武一223霹N蠢一2N,z6武吼4+z1N;瑶∞;一2只z6缈14∞22+碍z4国?+Z8co(co;+嵋z4∞;一2N2N3最z4co,co,+2Ⅳ{虬只只z2(-Oi032—2Ⅳ3虬只£z2q必+只2胛Z4国?q2—2e#?z4脚知;+4NsP324脚?国;+2N3N4只z2CO;一2E2N3Z4珊?珊;一2N1N4e:p3z2印i+2N4只只z2茸十2N4B置z4co,c02+2NlN2只z4∞;∞+2N2罡P424群国;一2NlN2P,p224∞;∞;-2PlP4N3N222∞;一2N2BP224cai}c02+2N3N2P3P222(91(-02+2N,N。马z4co,co;+碍Ⅳ;z4∞?∞;+2N2马z6司珊;一2N4P4z4国?脚;一2Ⅳ1N2只P422q国2+Ⅳ;尸32—2N3N4只尸4+Ⅳ:巧)以上各式子中有q=国:=Q,2.考虑Q:=l的系数A。=一(E珊?国;一64碍Ⅳ3一十16口E∞卜4Ⅳ3E印;一64孵e,Ml+16N2NiMlB国l哪2+16N2NsMlPfo;+64NsN4B^毛+64N4只只E一16N4只且f103i一8只E国÷珊;一16N。N3Ml功;一16N3EE2珊;一16N4只E珊;+32N3只E∞;+4只2E研国;+4N4Ml∞?∞;一16Ⅳ;只—瞄l∞;+16N4置P2E国l∞2—16N4只Ⅳ1MI(-OIO)2。+4N2P2Fln)1∞23—409l∞;Ⅳ1N2Ml一4N2只Ml棚;国;+16Nze,e,F1co;一16N:马最E国。珊:),(一128P3Ⅳ;国;+64P、2Ⅳ;∞;一8B国?棚;一128碍N,∞?+4jv?国?国:+4只2国14∞;+64N;碍彩;一8N3cot2co:+64P于.v?∞?+64P字川珊?+16Ⅳ;棚24+SNaP.co,co;一32P,2Ⅳ3珊?国;+12"SN4P4Nsco;+16Ⅳ?最2国?山;一512N4PsP4Ns+128N4P3P40)?一32Ⅳj只脚?脚;+64弓Ⅳ3研国;一32Ⅳ?马国?国;+16P2zNt2∞;∞22十16砖回:+256P42Nj+256P32Nj+国14∞24—32}IN2N3∞l∞;一128N4Nl马只∞?+32N4NlB研国:+32Ⅳ4丑£国2卯?一128Ⅳ一Ⅳ,BE∞l国2+128Ⅳ4Ⅳlp,elOal国2+32Ⅳ2Ⅳ1只q嗣一32Ⅳ1Ⅳ2最舅研喀+32N2墨只q2印22—128N2NsP4Ptm;+128Ⅳ3N2BB甜l∞2—128N1N2B只卯lC02—32N2B最甜2珊?)75附录B1=-8(4N2N3Ml脚;一16N2N3P39ic02+16PlN4N3M1092—4最MlⅣ;f010)jCl=-4(一16N2N3~4N2只曩∞;一N2Ml∞?∞;+4e,U2NlMl盘、∞;+4BBN2E卯l∞;+16N2只P4曩国2+4只Ⅳ2M1国2国:~N1Eq珊;一4只2NlE国1国;一4N4N1M1观∞;一4最N2E甜l曲;+8只NtE∞l扣;一4鼻N4Ml珊2国?一16只只N4E珊2—16P2M1.Ⅳ:国l+16BBⅣ4E国l一16口Ⅳ1Eq+16P3N4NlMl国1)/(一128P3埘珊:+64.Ⅳ;鼻2出;一8只卯?棚;一128只2N3∞?+4Ⅳ?砰甜;+4p,2珊14∞22+64Ⅳ:砰m;一8N3∞12∞24+64P32Ⅳ?∞?+64Ⅳ:只2国?+16N知o;+8N2马ogo,i一32弓2N3∞12∞22+128N4N3只珊;+16E2U?og&呸2—512N4N3B只+128N4B只彩?一32只N4090,;+64Ⅳ3只o,120,;一32只研ogo,;+16碍Ⅳ22q2钙2+16口硝+256巧Ⅳ;+256碍孵+q4印24—32N2虬Bq国;一128N4Nl最与国?+32.N4NlBo,120.7;+32N4鼻最国2国?一128N4N3鼻与国l国2+128N4NIE只珊10)2+32Ⅳ1N2/',col09;一32日BⅣ2N1研国22+32N2鼻只∞?脚;一128N2Ⅳ3弓尸4∞;+128N2N3最只∞1国2—128N2Nt只只脚l珊2—32N2只B国2国?)一16蝴og+32N3只Ml∞?一16Ⅳ12£MI珊?+64N3B只E+64Ⅳ3Ⅳ4只MlP2^彳1∞l国2—64N;MlB一4最置E∞?∞2—64N4E坪一4只M1科-16N4只Mlco?+4只巧o,120,;+4Ⅳ12Ml研国;一16N3P4Flco;-8N3MI珊?∞;一16N4曩国?譬+16N3BEE国10.72+16N4NIB^4ICO?一16N1只BE∞10.72-4NlBE∞-%2+Ml脚14卯22+16Ⅳ;Ml珊;+16N1最BE国?+16ⅣjⅣ1只^毛∞I印2+4Ⅳ2BA,l国?曲2)/(-128Ⅳ;B珊;+64Ⅳ;E2m;一8只国14卯22—128P32^r3砰+4Ⅳ?国?∞;+4e,2∞?国;+64妮砰国;一8N3珊?吨4+64∞?学Ⅳ?+64Ⅳ;只2∞?+16N;o:+8N2B司∞;一32鼻2N3国?国;+128N3N4只∞;+16Ⅳ?墨2彩12∞22-512N3N4只只+128N4只只研一32N4只国?吐E2+64N3Bogo,;一32Ⅳ?马国;∞;+16N;巧国?脚;+16碍∞?+256P42Ⅳ;+256/'32孵+棚?国:一32N3N2马珊;∞l一128N4NIBB∞?+32Ⅳ4Ⅳ1B茸脚;+32N4E马卯?∞2—128Ⅳ3N4E只脚10,/2+128Ⅳ4NIE只埘1092+32P4NlN2脚23国l一32NIN2只只印12鲫22+32N2只只∞?∞;一128N3N2只B国;+128N3N2B最曲l珊2—128NlN2B只q国2—32N2B最∞?∞2)76附录Dt=8(一4Ⅳ1N2置Mtco,c02—16N2N3J_肼l啦+4N2曩埘l奸国2+4Ⅳ2互碍群09'2+16FlN2砰∞2—16N3N4M1最∞l+16N3B只巧∞l+BE脚?∞;一8Ⅳ3Ml鼻国?国2一16Ⅳ3鼻异只国2—4鹄巧最国;q+16N;只Ⅳl国2+片J:l毛co?c02+4墨只互国?甜2—4P2P3曩o?+4.Ⅳ?置MI甜?珊2—4NlE鼻B珊?国2+4NlE只珊;∞1—16Nl曩只马国l+16N1N4只凹I劬十4MMfB研)/(.128P3埘曲22十64暑2孵国;-8Bo,:04—128口Ⅳ3国?+4Ⅳ?∞?脚;+4置2∞?瑚;+64P42Ⅳ;卯;一8N3co?co;+64学Ⅳ?珊?+64明只2脚?+16埘脚;+8N2£研∞;一32Ⅳ3皇2甜?q2+128N4Ⅳ3只曲;+16N12P12国?∞;一512N4N3P3P4+128N4只只∞?一32N4只研甜;+64N3只∞?国;一32ⅣJ2B∞?珊;+16婀只2∞?吨2十16学c04+256N:砰+256埘碍+国?嗣一32N3N2B国;国1—128NIN4最B国?+32NlN4尸2国?∞;+32Ⅳ4B鼻讲∞2—128N3N4EE国1092+128NIN4ele,oJl∞2+32NlN2只回;国1—32NlN2尸2只090.7;+32N2只E国12∞22—128N3N2鼻P,co;+128N3N2最只03l由2—128P3P4N2NI∞lC02-32N2P2马∞?∞2)A:=一(碍E研一碍RⅣ3一只M:N;一F2N3co;+EⅣ:Ⅳ,M:棚;+只日N2E印;一EN2BB珊l国2一ⅣzⅣ1M2COI应E3+Ⅳ2NIM2只国1602一只Ⅳ4NlM2COI卯2+EBⅣ2q嗣+F2N4P2e,coI印2+五脚细;一暑Ⅳ2M2霹印;一只M2孵脚;十只B.v4疋一N4—心M2∞;+Ⅳ4Ⅳ3B^,2+2只Ⅳ,&co;一置2N3‘国;-P,U。E∞;一2B疋珊抽;+丑2如研脚;十Ⅳ。M:国?国;一BⅣ。M:og)/(-2N;与∞;+日2埘∞;一2只缈;珊?一2口Ⅳ,砰+研∞?脚:+E2脚?∞;+孵只2∞;一2Ⅳ3国知;+学研国卜埘毋司+孵缈;+2Ⅳ:B研∞;一2Ⅳ3只2国?司+2N4N3只珊;+Ⅳ?E2研m:2—2N4N3只P4+2N4只B国?一2N4只国?m;+4Ⅳ,只∞?∞;一2Ⅳ?只砰国;+Ⅳ;碍∞÷∞:+只2国?+埘巧+碍N;+co:co;一2BⅣ2Ⅳ3∞l出;一2N4NI最B研+2Ⅳ4Ⅳl尸260260;+2N.墨B硝∞2—2N4.Ⅳ3B县∞10.72+2Ⅳ4Ⅳ1E只讲l珊2+2^7j|Ⅳ2只国1∞;一2EBN2Nl国王回22十2N2E尸4衅珊;一2N2N3E只∞;+2Ⅳ2Ⅳ3县只国l国2—2N2NtB只q国2—2N2只B叫3032)B:=一(BⅣ。E国103;+最BⅣ。E棚l一最M:Ⅳ;甜,∞;一尸2^彳2Ⅳ:卯。+只BⅣ2Eq国;C2=一(埘M2∞22一Ⅳ2最M2Ⅳ3q国2一只疋ENlO)1(*02+N2只M2Nl∞1032一只只Ⅳ4E国2+鼻Ⅳ4N3M2∞2一EⅣ4M2902国?一只2N1E珊l∞22+EⅣ2Ⅳ1M2∞l脚;+B^r4ⅣlM2q+只Ⅳ2M2国2脚;+2只N1F:0109;一N2N3B^彳2彩2+Ⅳ2只只置脚2一ⅣIF:ol∞24一N2M20)】2coi+Ⅳ2N3M290;一Ⅳ4N1M2q09;一碍ⅣlEq—N:只岛国;),(之B埘咙2+埘只2CO:一2B国?∞;一2学N3奸+斫q2“24+E2研四+孵砰国:2—2N3Co,co;+呼析出】2+孵只2国卜孵出:+2N2只国?国;一27,2Ⅳ3∞?国;+2N4Ⅳ3P:o;+鼻2Ⅳ?国?脚;一2N4N3B只+2N4只只∞?一2只Ⅳ4薪co;+4N,B∞知;一2B研∞细;+霹州2w.2脚;+晋(-0卜坪川+口—嵋+国?国;一2N2N3只国1国i一2N4N,P3P2(o;+2N4N1最珊?卯;+2N4墨只∞2国?一2N4Ⅳ3弓只珊l国2+2N4N1暑只珊l∞z+2NIN2只∞1∞;一2墨忍Ⅳ2NI奸印;+2N2只蜀国T2∞22—2Ⅳ2N3置墨国;+2N2Ⅳ3与冀(010)2—2N2N1只只q国2—2Ⅳ2BB脚2国?)+B疋鼻N3国l∞2+BENi只∞j一巴EⅣ1脚?《+Ⅳ4J】lf2ⅣlB讲一BEE∞?国2+Ⅳ2最^,2∞?m2一N4只朋’2∞;+2N3只^f2珊?一2N3M≯武∞12一噬N4F2吱÷Nj^4蕊畦一NjMzP3+N3N4P、M!+时2科国;一BM:国?一N,只‘印;+Ⅳ3只只E一只只疋国?+只E∞?珊;一N:M:P:o?一N4E碍)/(一2—V;B国;+Ⅳ;鼻2国;一2B硝国;一2碍鹄国12+Ⅳ12090,;+E2脚?国;+孵日脚;一2Ⅳ3脚?∞;+只2Ⅳ?砰+口Ⅳ;卯;+Ⅳ;国;+2N2B珊?∞;-2e,2N3co,co;+2Ⅳ3N4只国;+只2Ⅳ?脚?q2—2N3眠只B+2Ⅳ4只B甜12—2N,只薪喀+4Ⅳ3只q2“:2—2只Ⅳ12卯知;+gⅣ;国知;+口脚?+砰聪+口埘+co,co;-2P2N2N30l国;一2P2P3N4N1∞?+2足Ⅳ4N1砰曲?十2N4尸2毋∞?国2—2Ⅳ3N4最B田l噍,2十2N4N1只E∞l国2+2N2N1只ml∞i一2N2NlBE∞?&7主2+2Ⅳ2只鼻研珊;一2N2N3只鼻珊;+2B只N2N3091032—2Ⅳ2ⅣI与只国l国2—2N2县B研口k)D2=(e,e2F2U30J'1一只只‘ⅣlC01一P3e2F2co?+置^如N;o)2一只^彳2Ⅳ2Ⅳ3∞2一日只F2N3脚2—2P,丸t2N3印?珊2一e2M2N4N3国1一P2BN3卯l埘;+F21V2砰(02+e,P4Y20J,2c02十只Ⅳ2NlU,colte,92N2翻?嘞+只五Ntq彰一最^f2Ⅳ2Nl彩?∞2+BA42N40J?一E马EⅣl缈?珊2+,2Ⅳ2只2拜国2+只A彳:Ⅳ;∞?国:+最‘珊?印:+只M:国?国:),(一2BⅣ;珊:+Ⅳ;只2出;一2只科∞;一282N,CO?+N?脚?国;+鼻2国14∞22+Ⅳ;砰甜:一2N3国?嘭4+砰N?og+川巧甜?+Ⅳ;∞;十2N2P2国?珊;一2N,B2∞?珊;+2N3N4只国;+Ⅳ?鼻2∞?出:一2N3N4BP4+2N4Be,ro?一2N4e,0960;+4N3&回;国12—2墨N;∞、2哇+Nj心武哇+嫒∞:+Ni磺+N;霹+国?缈;一2&N2N30J103i一2NN4P3P2co?+2Ⅳ1Ⅳ4B∞?∞;十2BBⅣ4珊?国2—2N4N3B只棚1国2+2NlN4P1只O/l∞2+2N2Nl&c01缈;一2N2N1尸2鼻棚12脚22+2N2E只∞;递一2N2Ⅳ3置.Ⅳ403;+2P2只Ⅳ2N3∞IC02—2P4只N2Nl国IfrO'2—2N2Be,Cc02)A3=一4(4N2马只z4山l∞;一4NlN2M324脚ICO;-4N2只^如z4研曲;+16Ⅳ3N2最^如z2CO;+16N2暑只Ez2讲;+16NlⅣ2只Af322∞I∞2—16N2马只Ez2舡,;一16只M3Ⅳ;z203;+4鼻2Ez4珊;国;一8Z4e,Eog面+Ez6m?国;一4.v3Ez4曲;+4N4M324砰国;一16&U4Ez2国;一16N4N]M322国22+32BⅣ3Ez2脚;一16e,2N3Ez2∞;一16Ⅳ4Jvl只坞Z2由l珊2+16Ⅳ4县鼻EZ2∞l印2—16P3N4M322脚?+64Ⅳ4B只巧+64Ⅳ3Ⅳ4马M3+16口Ez2国12—64孵只坞~64学N3^)/(64Ⅳ;曰z2∞;一8N326脚12删24—8只z6国?D;+4e,226国?∞;+4N?z6珊?珊;一128BⅣ;z2脚;+16N;Z4吨4+16P3224科+128帆N4只z2嗣一128ⅣlN4只最z2印?+32Ⅳ4B最z4国沁2—128Ⅳ3N4鼻Bz2叻珊2+64川只222097+64研只222国?一128N3卑z2q2+128N4P3P4220)?+64N3Bz4珊?珊:一32BN?z4国?国;+64矗2Ⅳ;z2国22+128.ⅣlN4E只z2珊l彩2+8N2曼z6棚?m;一32Ⅳ3N2只z4甜l噍蠢-128N,N2只P4z2国;+32N1N2P424∞23峨+32N2只P424∞?∞;-32N2P2Bz4国?呜-32NlN2丑Bz4国12讲22+16Ⅳ;碍z4国;∞;一32E2N324砰珊;+128N2N3BBz2∞l∞2—128NlⅣ2B只z2国l甜2+32NIN4昱.z4埘;《一32N4P_424脚;脚;+16舅2N?Z4国?脚;+z8国14Ⅲ24—512N3Ⅳ4B只+256Pd2峨+256碍Ⅳ;)B3…8z(N2M324砰∞;-Nl‘z4埘l国;-4N4NlM3z2国l珊;一16M3岛Ⅳ;脚l+4N2N3M322棚;+16N4EBB国l一4Ⅳ;^彳3只z2国1CO;一4N4EBzl出l埘22+4N2EE最z20Jl∞;一4NtEE222COl国;+4N2NIM3e,z2C01国22+16N4Ⅳ3M3只国2+16N4NlM3BmI+4N2M3只z2C020)12+16N2‘只只∞2—16N2N3M3B珊2—16NlE碍彩l+8NlEBz2∞1∞;一4N2E只z2∞;一16N4巧鼻只国2—4N4M3鼻z2∞209)/(64砰Ⅳ;z2脚;一8N326脚12Ⅲ24—8Bz6卅m;+4P,226脚14Ⅲ22+4N2,Z6卯?印;一128BⅣ;z2CO;+t6N;Z40.1:+16只224CO?+128只N3N422珊;一128鲳NlN422∞?+32P,最N424C02脚?一128e点U3N422国lC02+64碍Ⅳ;z2∞12+64碍N?Z2国?一128P32N322CO?+128粥N4z2∞?+64P3N3z4脚?∞;-321v?只z4回12∞22+64只2孵z2CO;+128墨只ⅣIN4220)lO)2+8芝Ⅳ22609嗣-32N2Ⅳ3马z40口l/v;一128犯N2N322∞;+32P4NIN2Z4COlCOi+32e,只N224∞?国;一32e2e,N224C02∞?一32鼻最ⅣlN224脚?∞;+16Ⅳ:碍z4珊;珊;一32鼻2N]Z4珊;脚;+128BBⅣ3N2Z2∞l∞2—128NlN2B只z2∞l甜2+32P2N1N424珊12∞22—32P4N424∞J2Ⅲ22+l§最2孵z409国22+Z8珊J4∞24—512P3P心斟4七256N;Pj七256N;P;、+16NlBBEz2∞?一16N4碍Ez20,2—16N4只M322棚?一16N?马M3220,12—8N3M324∞?∞;一16Nl置只玛z2∞I∞2+16N3日B只z2C01C02+4只J_z4∞?珊;一16N3只Ez2珊;+4N2BM,z4∞?国2—16NjN2BM3Z2珊l脚2+32N3£且厶Z2珊?+16ⅣlⅣ4是M3Z2辞+t6NlNz蜀M322国l吐一4置与Ez4col30)2+64N3%只^毛+64弘只只E+16埘M,z2co:+鸠z6耐∞;一4BM,Z4科)“64巧Ⅳ;z2国;一8N326国?国;一8Bz6埘?国;+4鼻226硝国;+4N?Z6珊?棚;一128P,孵z2国;+16埘Z4∞24+16P3224司+128Ⅳ3Ⅳ4只z2国22-128NiⅣ4最弓z2国?+32N4BEz4090,2—128N3N4最Bz2C01国2+64只2Ⅳ;z2妍+64P32N?Z2国?一128N3譬z2fO?+]28P3&N4Z2彩?+64Ⅳ3只z4碍脚;一32只N?Z409q2+6们2叼z2脚22+128NfN4鼻只z2脚f0)2+8N2最26研脚;一32Ⅳ2J^,3县z4埘23q一128M—Ⅳ3鼻只z2∞;+32NIN2尸424嗣国l+32N2鼻只z4国12∞22—32N2BBz4研国2~32NIN2只马z4甜?卯;+16碍明z4回国:一32E2虬z4国?国;+128Ⅳ2Ⅳ3最Bz2国l∞2—128ⅣlN2只只z2国l珊2+32N】N4Bz4印?印;一32N4只z4国?卯;+z8回国;+16E2Ⅳ?z4珊12珊22—512N_|N3只B+256P42Ⅳ42+25682Ⅳ;)C3=一4(一64Ⅳ;B^如一64N。曰E一4N.只巧z4∞?功;+4Ⅳ?J】l毛Z4砰国;一16B只Ez2国?附录D3=8z(一8N3置Jjl如Z20)120)2+16N3B马E国1—16N3鼻只E国2—4N3BEz2棚22q+4PlP4F322∞?∞2+巴‘z4∞?珊;+E^毛Z4脚?珊2—4P2P3F3Z2珊?+4Ⅳ4最^厶Z2国?+16Ⅳ;舅M3∞2+4N2只M322co,c02—4N2Nl马M322co(e02—16N2N3只肘3吼+4N2F3e:z2co,c02+4Nle,F322∞22国l一4Nl鼻最Ez2∞?∞2+4Ⅳ?鼻^厶Z2c012c02—16Ⅳ1B只Eq+16N4N1只M3091+16N2E碍∞2—16N4N3BM3q),(6422Ⅳ;砰脚22-8虬Z6唧12∞24-8P=;Z6国14∞22+4鼻226群国22+4N(Z6m12Ⅲ24一128只Ⅳ;z2国22+16N;Z4∞24+16e?z4CO?+128N3N4只z2珊22—128N1N4与只z2国?+32N4只鼻z4co,c02—128N3N4鼻最z2COl国2+6422露.v;珊?+6422学Ⅳ?砰一128Z2砰N3CO?+128N49e,z2国?+64N3只z4∞?∞;一32e3N?Z4卯?珊;+64P,2—嵋z2∞:+128NlN4鼻只z2COlC02+8N2P226co,co;一32N2N3Bz4珊;∞l一128PtP4N3N222∞;+32N1N2P424珊23国I+32N2鼻尸424∞12卯22~32N2只最z4脚?∞2—32ⅣlⅣ2只芝z4砰哆2+16只2孵z40,2q2—32暑2码z4四面+128Ⅳ3Ⅳ2B昱z20J'10)2—128NlN2B只z2∞l∞2+32N1N4巴z4砰功;-32N4只z4国?国;+167,2N?z40J?co;+Z8硝国;一512N3N4只只+256N4z砰+256N;只2)以上各式子中有国l=∞:=Q:8】致本论文是在我的导师丁千副教授的悉心指导和热情鼓励下完成的。从论文的选题到论文的完成无不凝聚着导师的心血。导师在工作中循循善诱的教诲、时时的鼓励与亲切的关怀陪伴我度过了两年多紧张而充实的学习生活,导师渊博的学识、严谨的科研作风和平易近人的品格给我留下了极为深刻的印象,为我树立了人生的榜样。我的一点一滴的进步都与他的关怀鼓励和辛勤栽培息息相关。能够聆听丁老师的谆谆教诲,是我人生的一大收获。在此表示衷心的感谢!在论文完成的过程中,很多老师提供了许多帮助。感谢叶敏教授、张琪昌教授、曹树谦副教授在论文研究过程中给予我学术上的帮助和教诲。感谢刘习军教授、贾启芬副教授、吴志强副教授给予我的帮助和教诲。感谢郎作贵高级工程师、张文德高级工程师给我提供的学习环境和1帮助,感谢硕士杜颖、谭海波、陶志军、马红城、王立刚、刘海英、任爱娣,在科研工作中的互相帮助、互相鼓励、取长补短,不断提高科研水平,对此深表谢意。感谢舍友于雯、王彦青、邸玉贤对我的帮助和鼓励。感谢师弟王冬力、张建伟在实验室的真诚相处和热情帮助。我要深深的感谢父母,是他们的理解、支持和鼓励让我有力量一直继续下去,每当我在学习和生活上遇到困难和挫折的时候,是他们的关心和爱护使我重新信心百倍,鼓起奋进的勇气。感谢姐姐一直以来对父母的照顾,感谢懂事的弟弟,让我远在他乡少了几许牵挂。借此机会,还要衷心感谢在百忙中抽出宝贵时间对论文进行评审工作的专家和学者们。

因篇幅问题不能全部显示,请点此查看更多更全内容

Top