考虑环境参数的风电齿轮箱传动系统疲劳性能优化.pdf
《考虑环境参数的风电齿轮箱传动系统疲劳性能优化.pdf》由会员分享,可在线阅读,更多相关《考虑环境参数的风电齿轮箱传动系统疲劳性能优化.pdf(11页珍藏版)》请在咨信网上搜索。
1、2023年 第47卷 第8期Journal of Mechanical Transmission考虑环境参数的风电齿轮箱传动系统疲劳性能优化谭建军1,2,3 伍 源2 白 杨2 朱才朝2 宋朝省2 董晔弘1(1 中国船舶重工集团 海装风电股份有限公司,重庆 401122)(2 重庆大学 机械传动国家重点实验室,重庆 400044)(3 重庆大学 电气工程学院,重庆 400044)摘要 随机风速会使风电齿轮箱传动系统出现频繁的载荷波动,造成复杂的结构变形,容易产生齿轮偏载,加剧其接触疲劳失效风险。提出了一种考虑环境参数随机不确定性的风电齿轮箱传动系统疲劳性能优化方法,在建立计入全局载荷的风电齿轮
2、箱传动系统动力学模型基础上,利用代理模型方法重构“平均风速、湍流强度-齿轮修形参数-齿轮长期接触疲劳损伤”映射关系,建立考虑风速概率分布的多级齿轮修形参数优化函数,对比了风电齿轮箱传动系统疲劳性能优化效果。结果表明,低速级太阳轮长期接触疲劳损伤值大于0.7,是风电齿轮箱传动系统高可靠设计的薄弱环节之一;优化后的风电齿轮箱齿轮长期接触疲劳损伤值明显降低,其中低速级太阳轮长期接触疲劳损伤值降低了11.37%,疲劳性能提升效果显著。关键词 风电齿轮箱 传动系统 环境参数 疲劳损伤 性能优化Fatigue Performance Optimization of the Wind Turbine Gea
3、rbox Transmission System Considering Environmental ParametersTan Jianjun1,2,3 Wu Yuan2 Bai Yang2 Zhu Caichao2 Song Chaosheng2 Dong Yehong1(1 CSIC(Chongqing)Haizhuang Windpower Equipment Co.,Ltd.,Chongqing 401122,China)(2 State Key Laboratory of Mechanical Transmissions,Chongqing University,Chongqing
4、 400044,China)(3 School of Electrical Engineering,Chongqing University,Chongqing 400044,China)Abstract Random wind speed causes frequent load fluctuations in the wind turbine gearbox transmission system,resulting in complex structural deformations and gear deflections,which aggravate the risk of gea
5、r contact failure.A novel method to improve the fatigue performance of wind turbine gearbox transmission systems considering the stochastic uncertainty of environmental parameters is proposed.The mapping relationship among wind turbulence intensity-gear modification parameter-gear long-term contact
6、fatigue damage is reconstructed by the surrogate model method based on the dynamic model of the wind turbine transmission system considering global loads.And then,a multi-stage gear modification parameter optimization function considering the probability distribution of wind speed is established and
7、 the optimizing fatigue performance is compared.The results show that the long-term gear contact fatigue damage value of the low-speed stage sun gear is more than 0.7,which is one of the weak parts of the high-reliability design of the wind turbine gearbox transmission system.The long-term gear cont
8、act fatigue damage values of optimizing gears have been reduced significantly,in which the low-speed stage sun gear has been reduced by 11.37%,so the fatigue performance has been improved significantly.Key words Wind turbine gearbox Transmission system Environmental parameters Fatigue damage Perform
9、ance optimization文章编号:1004-2539(2023)08-0053-11DOI:10.16578/j.issn.1004.2539.2023.08.00853第47卷0 引言我国在“十四五”期间风电装机容量预期新增250 GW138,发展前景巨大,然而随着国家风电补贴政策的全面取消,对风电机组的度电成本与可靠性提出了严苛要求。风电齿轮箱是风电机组传递力与运动的关键传动装置,成本约占整机的16%,是风电机组提质增效的关键环节138-39。风电齿轮箱设计寿命通常为20年(陆地)或25年(海上),全寿命周期内平均风速与湍流强度的随机变化会使风电齿轮箱面临频繁的输入载荷波动,容易
10、造成多级齿轮早期接触疲劳失效,实际服役寿命远低于设计值2,如图1 所示。因此,开展考虑环境参数的风电齿轮箱传动系统疲劳性能优化,对于提高其对环境工况的适应性、保障风电机组安全运行具有重要意义。近年来,国内外学者针对风电齿轮箱传动系统的疲劳性能分析及优化开展了大量研究。建立高保真的齿轮传动系统动力学模型,是准确计算齿面动载荷、实现疲劳损伤分析及性能优化的重要前提。Helsen等3建立了3种保真度风电齿轮箱传动系统动力学模型,即纯扭转模型、刚性6自由度模型和柔性多体动力学模型,分析了结构柔性化程度对系统振动特性的影响。Guo等4对比不同保真度仿真模型的求解精度与效率,给出了风电齿轮箱传动系统动力学
11、建模时部件刚/柔、齿轮修形以及制造安装误差等建模参考标准。许华超等5-6考虑内齿圈柔性,分析了齿轮传动系统振动模态与动态响应。翟洪飞等7采用弯曲Timoshenko梁单元,建立了箱体子结构模型,分析了其对齿轮传动系统振动特性的影响。Tan等8-9将主轴和齿轮箱传动轴柔性纳入传动系统动力学模型中,分析了其动态特性。孙秋云等10308-31511分析了行星轮销轴误差对风电齿轮箱传动系统动态响应和疲劳寿命的影响。风电齿轮箱传动系统常在随机气动载荷驱动下运行,属于受迫振动,其运行状态不仅与构件柔性变形和内部误差激励等相关,还取决于风速工况。在风电机组全寿命周期内,长期风速概率分布会使风电齿轮箱传动系统
12、出现复杂的载荷波动与结构变形,容易造成齿面载荷分布不均1211,加剧齿轮疲劳失效风险。Xiang等13分析了平均风速和湍流强度对风电齿轮箱各级齿轮短期疲劳损伤的影响,指出平均风速变化会显著改变齿轮疲劳损伤值。Jorgensen等14结合平均风速威布尔分布,分析了风电齿轮箱各级齿轮长期疲劳损伤情况,探讨了湍流强度对风电齿轮长期疲劳损伤的影响。Wang等15通过建立10 MW级风电齿轮箱传动系统动力学模型,结合平均风速威布尔分布,发现行星轮系太阳轮和定轴轮系小齿轮的长期疲劳损伤值较大。为了改善风电齿轮箱传动系统动态性能,Park等16147-151分析了齿轮修形参数对齿面载荷分布的影响,得到了额定
13、工况下的最佳齿轮修形方案。Wei等17314-335分析了齿轮修形参数对啮合刚度与系统振动响应的影响,并通过台架试验对齿轮修形结果进行了验证。Shu等18797-816采用代理模型技术,获取了适用于多种稳定工况的齿轮传动系统最优轴承支承刚度。刘华朝等19158-163以齿轮接触应力最小为优化目标,基于载荷谱对风电齿轮箱各级齿轮修形参数进行了优化。陈岩松等121-14建立了考虑不同风速工况组合的风电齿轮箱行星级均载优化模型,优化后的风电齿轮箱行星级在多种风速工况下均有良好的均载性能。虽然上述研究1213-1416151173341881519163对风电齿轮箱传动系统高可靠设计提供了一定的理论指
14、导,但相关研究多基于单一或多种恒定输入转矩工况开展齿轮传动系统减振优化设计,忽略了因长期风速概率分布造成的风电齿轮箱输入转矩“短期随机波动”和“长期概率分布”的风电特征,难以保证优化后的风电齿轮箱传动系统在全寿命周期内性能达到最佳;同时,全寿命周期内工况数量多,开展全工况仿真计算疲劳损伤耗时长,计算效率低。本文以某型 5 MW 级海上风电齿轮箱为研究对象,考虑随机风速与传动系统拓扑结构,建立了计入全局载荷的风电齿轮箱传动系统动力学模型,基于代理模型重构了“平均风速、湍流强度-齿轮修形参数-齿轮长期接触疲劳损伤”映射关系,建立了考虑风速概率分布的多级齿轮修形参数优化函数,对比了风电齿轮箱传动系统
15、疲劳性能优化效果。1 风电齿轮箱传动系统运行原理如图2所示,风电齿轮箱主要由低速级、中间级和高速级组成。其中,低速级和中间级均为斜齿行图1风电齿轮箱内齿圈齿面剥落Fig.1Early fatigue peeling of ring gears in the wind turbine gearbox54第8期谭建军,等:考虑环境参数的风电齿轮箱传动系统疲劳性能优化星轮系,主要由太阳轮(s1、s2)、行星轮(p1、p2)和行星架(c1、c2)等组成;高速级为斜齿定轴轮系,主要由大齿轮(g1)和小齿轮(g2)组成。叶轮的气动转矩、弯矩、径向载荷和推力等气动载荷通过主轴传递到低速级行星架,进而依次通过
16、低速级行星轮-太阳轮啮合副、低速级太阳轮-中间级行星架内,花键连接将力与运动传递到中间级行星架;同理,驱动中间级太阳轮将功率传递到高速级大齿轮;最后高速级大齿轮-小齿轮啮合副使高速级输出轴驱动发电机进行发电。2 计入全局载荷的风电齿轮箱传动系统动力学模型图3所示为计入全局载荷的风电齿轮箱传动系统动力学模型,主要由两部分构成:OpenFast风电机组整机系统动力学模型和Simpack风电齿轮箱传动系统动力学模型。采用OpenFast建立风电机组整机系统动力学模型,主要由空气动力学模块、结构动力学模块和控制模块组成。空气动力学模块根据平均风速和湍流强度等环境参数,结合叶片气动参数和桨距角等计算作用
17、在叶轮的转矩、弯矩等气动载荷;结构动力学模块包括叶片、塔筒、齿轮箱和发电机等,其中齿轮箱被简化为传动比;控制模块包括偏航控制、桨距角控制和发电机控制,详细模型介绍见文献20。本文通过改变空气动力学模块的平均风速和湍流强度,计算叶轮的转矩、弯矩等气动载荷。采用Simpack建立风电齿轮箱传动系统动力学模型,各齿轮参数如表1所示。其中,传动轴、行星架和箱体采用有限元获取全局质量与刚度矩阵,然后利用模态综合法进行缩聚建模21;齿轮啮合采用可以考虑修形的FE225号力元模拟,轴承支撑采用FE41号力元模拟,其中,轴承支承刚度利用Romax软件计算10308-309,齿轮啮合刚度则根据 ISO 6336
18、 系列标准计算。如图 4 所示,为了将 OpenFast 风电机组整机系统动力学模型和Simpack风电齿轮箱传动系统动力学模型进行联合仿真,本文通过编制工况批处理Matlab程序,动态修改OpenFast空气动力学模块的输入工况.inp文件,计算给定平均风速和湍流强度条件下叶轮的转矩和弯矩等6自由度气动载荷;然后利用Matlab将气动载荷文件格式转为Simpack风电齿轮箱传动系统动力学模型的载荷输入.afs文件,再通过调用宏命令.sjs文件修改风电齿轮箱传动系统动力学模型中齿轮修形参数并控制模型仿真,得到风电齿轮箱传动系统齿轮副动态啮合力;最后,计算了各级齿轮接触疲劳损伤。图2风电齿轮箱传
19、动系统Fig.2Wind power gearbox transmission system表1某型5 MW级风电齿轮箱齿轮基本设计参数Tab.1 Basic design parameters of a 5 MW wind turbine gearbox参数齿数模数/mm螺旋角/()压力角/()低速级内齿圈9324520行星轮2924太阳轮3224中间级内齿圈118178.520行星轮4717太阳轮2317高速级大齿轮12112920小齿轮2412图3计入全局载荷的风电齿轮箱传动系统动力学模型Fig.3Dynamic model of the wind turbine gearbox tra
20、nsmission system considering global loads55第47卷3 考虑风速概率分布的齿轮长期疲劳损伤计算3.1平均风速-湍流强度联合概率分布利用二元Copula函数c(,)构建平均风速-湍流强度联合概率密度函数fvI(u,)22,其表达式为fvI(u,)=c Fv(u),FI()fv(u)fI()(1)式中,fv(u)为10 min时序风速的平均风速概率密度函数,通常服从威布尔分布,fv(u)=(u)-1e-(u);u为平均风速;为形状参数;为尺度参数;fI()为10 min时序风速变化程度的概率密度函数,通常服从伽马分布,fI()=()-1e-;为湍流强度;为
21、尺度参数;为形状参数;()为伽马函数;Fv(u)、FI()分别为fv(u)和fI()的累积分布函数;c(,)为二元Copula函数,本文选用Gumbel Copula函数。3.2齿轮修形方案齿轮修形通常包括齿廓修形和螺旋线修形两种形式。齿廓修形图5(a)主要由修形长度、最大修形量和修形曲线确定239。齿顶最大修形量ek为ek=(4+0.04Ft/bca)4(2)式中,Ft为与齿轮分度圆相切的圆周力;bca为有效齿宽。采用长修形方式,可得修形长度为long=Pb(-1)(3)式中,long为修形长度;、Pb分别为齿轮端面重合度和基圆节距。连接修形起始点到最大修形量ek的曲线为修形曲线,如式(4)
22、所示。e=ek(xl)b(4)式中,l为修形起始点与单对齿啮合上界点之间的距离;x为啮合位置相对坐标值;b为指数,取b=2。螺旋线修形采用修鼓形方式,以齿宽中线为对称中心,采用二次曲线进行修形2311-12,齿宽中间修形量最小,齿宽两端修形量最大,呈对称分布,如图5(b)所示。螺旋线修鼓量表达式为=0.25bca 10-3+0.5fg(5)fg=A(0.1bca+10)(6)式中,fg为齿轮螺旋线误差;A为系数。3.3基于代理模型的齿轮长期接触疲劳损伤计算采用拉丁超立方抽样法与最大差异法24抽取关键环境参数组合和齿形设计参数组合,作为第2节计入全局载荷的风电齿轮箱传动系统动力学模型的参数输入,
23、可以计算出任意环境工况(u,)与齿形设计参数(ek,)组合下的齿轮副时序动态啮合力,并利用疲劳损伤线性累加,计算齿轮接触疲劳损伤25,进而可得在给定u和时具有确定齿形设计参数的齿轮短期接触疲劳损伤,为Dc(u,es1k,s1,s1,eg2k,g2,g2)=ni(u,es1k,s1,s1,eg2k,g2,g2)C-mi=NbC0mif()d(7)图4OpenFast风电机组整机-Simpack风电齿轮箱联合仿真流程Fig.4OpenFast wind turbine model-Simpack wind turbine gearbox model cosimulation process(a)齿
24、廓修形(b)螺旋线修形图5齿轮修形Fig.5Gear modification56第8期谭建军,等:考虑环境参数的风电齿轮箱传动系统疲劳性能优化式中,i为第i个应力块幅值;ni为应力水平i作用在齿轮的实际次数;C、m均为齿轮接触疲劳S-N曲线参数;Nb为应力块数量;f()为应力块分布函数。在风电齿轮箱全寿命周期内,考虑u和出现的概率及持续时间,基于式(1)可得齿轮长期接触疲劳损伤为Dtol(es1k,s1,s1,eg2k,g2,g2)=u1u212TfvI(u,)Dc(u,es1k,s1,s1,eg2k,g2,g2)dud=u1u212TfvI(u,)Dc(u,es1k,s1,s1,eg2k,
25、g2,g2)u(8)式中,T为风电齿轮箱设计寿命。由 于 式(7)中 同 时 存 在 环 境 参 数 随 机 变 量(u,)与齿轮修形区间变量(-pc)(表示随机变量,-表示区间变量),难以直接求解,因此,本文基于代理模型技术分别在“内层”和“外层”构建“齿轮修形变量(-pc)-齿轮短期接触疲劳损伤Dc(u,-pc)”“随机环境变量(u,)-齿轮短期接触疲劳损伤Dc(u,pc)”近似函数,具体步骤(图6)为:Step 1:初始化确定随机环境变量(u,)和齿轮修形区间变量(-pc)范围。Step 2:在“外层”采用拉丁超立方抽样法和最大差异法对u和进行抽样,生成环境参数组合工况样本矩阵X|2 M
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 考虑 环境参数 齿轮箱 传动系统 疲劳 性能 优化
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。