CN102831260A - 一种风力发电机组轮毂强度计算方法 - Google Patents
一种风力发电机组轮毂强度计算方法 Download PDFInfo
- Publication number
- CN102831260A CN102831260A CN2012102620601A CN201210262060A CN102831260A CN 102831260 A CN102831260 A CN 102831260A CN 2012102620601 A CN2012102620601 A CN 2012102620601A CN 201210262060 A CN201210262060 A CN 201210262060A CN 102831260 A CN102831260 A CN 102831260A
- Authority
- CN
- China
- Prior art keywords
- wheel hub
- finite element
- blade
- fatigue
- calculation method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Wind Motors (AREA)
Abstract
一种大功率水平轴紧凑型风力发电机组轮毂计算方法。在有限元模型中考虑了变桨驱动、变桨轴承的滚子等,并使用了一段叶片假体。让模型的受力传递更加准确,模型计算结果的准确性提高。根据轮毂的受力特点,只对四个极限工况进行计算,提高了计算的效率。对轮毂体的疲劳分析,在轮毂体的表面铺设了一层很薄的壳单元,根据有限元分析中所获得的应力应变结果,结合材料的疲劳寿命性能,应用疲劳理论,计算出零部件的疲劳寿命。
Description
技术领域
本发明涉及一种发电机组关键零件的强度计算方法。具体说是一种风力发电机组轮毂的强度计算方法,主要用于大功率水平轴紧凑型风力发电机组。
背景技术
目前主流的大功率水平轴风力发电机组中,风载荷通过叶片和变桨轴承传递到轮毂,轮毂通过主轴将动力传递到主传动链进行发电。在轮毂上安装了变桨驱动,使风机在风力变化的情况下改变叶片的桨距角度,从而确保风机安全高效的运行。轮毂是风机中的关键零件之一,在设计过程中为确保风机的安全和使用寿命,根据德国劳埃德(GL)认证规范必须对其强度进行计算。目前,对轮毂强度计算方法主要是采用现有的有限元软件进行模拟,现有的有限元分析软件可用于覆盖结构、温度、流体、电磁场和多物理场耦合等研究领域,但是模型的简化方法对计算精度的影响较大,如简化不当甚至不能正确计算出轮毂的强度。
在专利检索中指发现有如下一些专利与本发明有一定的相关:
1、专利申请号为CN201010617957.2,发明名称为“一种风力发电机组主轴强度计算方法”的发明专利,该专利公开了一种风力发电机组主轴强度计算方法,通过现有的有限元软件为平台实施建模,最后定义各个部件材料属性,进行计算,得出变形及应力云图,建模过程中主要包括将主轴、轴承和锁紧盘用实体单元模拟设置合并为一体结构,锁紧盘外部设置四个齿轮箱弹性支承并且每个齿轮箱弹性支撑由弹簧单元组模拟,风轮节点与主轴的各个法兰螺栓孔通过刚性梁单元I连接,锁紧盘与弹簧单元组采用刚性梁单元II连接;于风轮节点施加外部载荷,对弹簧单元组施加全部约束,对轴承施加位置约束。
2、专利申请号为CN201110352047.0,发明名称为“兆瓦级风力发电机组轮毂有限元建模方法”的发明专利,该专利公开了一种兆瓦级风力发电机组轮毂有限元建模方法,根据球形轮毂的对称性,取其1/3的模型进行建模,形成1/3的有限元模型,然后旋转形成整体的轮毂有限元模型;利用球形轮毂的旋转对称性来进行计算和强度校核,旋转对称性是指:在其中一个叶片加载和在另一个叶片加载形成的轮毂应力云图,两者旋转120度或者240度,是能够完全重合的。
但是专利1只是针对主轴强度进行计算的,并不能适用于轮毂强度分析计算;专利2的不足是:1)没有对几何模型的简化进行说明;2)在建立有限元模型时没有对轴承和变桨驱动的模拟进行考虑,导致有限元模型不够准确;3)没有对需要计算的极限载荷工况进行总结筛选,导致计算工作量较大;4)没有对轮毂疲劳寿命的具体计算进行介绍 ,因此仍有必要对如何改善现有轮毂强度分析计算方法中的不足进行分析研究。
发明内容
本发明的目的在于针对现有轮毂强度分析计算方法中的不足,提出一种新的轮毂强度分析计算方法,此方法更为简单,精度更高。
为实现上述目的,本发明通过以下技术方案实现:
一种风力发电机组轮毂强度计算方法,主要用于大功率水平轴紧凑型风力发电机组,通过现有的有限元软件为平台实施建模,定义各个部件材料属性,进行计算,得出应力云图;建模过程中主要包括以下步骤:
(1)建立的有限元模型;包括:叶片根部的一段、变浆轴承、变桨驱动、轮毂、主轴靠近轮毂的一段,分别建立有限元模型;其中,在叶片根部创建节点并将其转换到德国劳埃德GL规范中规定的叶根坐标系下,然后刚性耦合到叶片上;变桨轴承的滚子采用只受压力不受拉力的单元进行模拟;叶片和变桨轴承、变桨轴承和轮毂、轮毂和主轴之间的接触属性均定义为绑定;在轮毂体表面铺设一层很薄的壳单元。
(2)对相关零部件所使用的材料在有限元前处理软件中定义其弹性模量、泊松比等属性。
(3)建模完成后,使用载荷软件(GH bladed)输出GL叶根坐标系下叶根极限载荷,将载荷表中Mx、My、Mz和Myz绝对值最大工况分别施加在叶根中心,并对主轴端部施加全约束。
(4)将有限元模型提交有限元软件进行分析计算,使用有限元软件后处理模块看结果云图并提取相关结果数据。
该发明所述的风力发电机组轮毂强度计算的有益效果为:1、在有限元模型中考虑了变桨驱动、变桨轴承的滚子等,并使用了一段叶片假体。让模型的受力传递更加准确,模型计算结果的准确性提高。2、根据轮毂的受力特点,只对四个极限工况进行计算,提高了计算的效率。3、在轮毂体的表面铺设了一层很薄的壳单元,供对轮毂进行疲劳分析时使用。
附图说明
图1为本发明计算方法所用几何模型结构。
图2为本发明计算方法所用有限元模型结构。
图3为有限元模型中变桨轴承和变桨驱动模拟的细节特征。
图4为轮毂在某极限工况下的应力云图。
图5为轮毂的疲劳损伤结果图。
具体实施方式
下面将结合附图和具体实施例对本发明做进一步的描述。
通过附图可以看出本发明涉及一种风力发电机组轮毂强度计算方法,主要用于大功率水平轴紧凑型风力发电机组,通过现有的有限元软件为平台实施建模,定义各个部件材料属性,进行计算,得出德国劳埃德GL规范中规定的叶根坐标系下叶根极限载荷Mx、My、Mz和Myz绝对值最大四个极限工况下的应力云图;建模过程中主要包括以下步骤:
(1)建立的有限元模型;包括:叶片根部的一段、变浆轴承、变桨驱动、轮毂、主轴靠近轮毂的一段,分别建立有限元模型;其中,在叶片根部创建节点并将其转换到德国劳埃德GL规范中规定的叶根坐标系下,然后刚性耦合到叶片上;变桨轴承的滚子采用只受压力不受拉力的单元进行模拟;叶片和变桨轴承、变桨轴承和轮毂、轮毂和主轴之间的接触属性均定义为绑定;在轮毂体表面铺设一层很薄的壳单元。
(2)对相关零部件所使用的材料在有限元前处理软件中定义其弹性模量、泊松比等属性。
(3)建模完成后,使用载荷软件(GH bladed)输出德国劳埃德GL规范中规定的叶根坐标系下叶根极限载荷,将载荷表中Mx、My、Mz和Myz绝对值最大工况分别施加在叶根中心,并对主轴端部施加全约束。
(4)将有限元模型提交有限元软件进行分析计算,使用有限元软件后处理模块看结果云图并提取相关结果数据。
实施例一
一种风力发电机组轮毂强度计算方法,使用三维软件建立几何模型,在建模前需要对轮毂的外部接口零部件和受力要有一个总体的认识。与轮毂1连接的主要零部件是变浆轴承3、变浆驱动、叶片2、楔形盘和主轴4。在建立几何模型的过程中需对强度影响很小的特征,如小的倒角、凸台等进行简化。建立完毕的几何模型如图1所示。将几何模型转化为中间格式的文件后,导入专业的有限元前处理软件,对模型进行网格划分并对相关零部件定义材料属性,以建立分析所需的有限元模型,所建立的有限元模型如图2和图3所示。为了便于载荷的施加,建立了叶片的假体2,根据圣维南原理,叶片假体2的长度保留大于叶根直径的一倍,并且分别建立了主轴的假体4和变桨轴承3的假体。楔形盘位于叶片假体2和变桨轴承3之间,予以简化,不建立假体。简化掉叶片假体2与变桨轴承假体3、变桨轴承假体3与轮毂1、轮毂1与主轴假体4之间的螺栓连接,在有限元模型中通过设置绑定来耦合各部件。对主轴端部进行全约束,施加载荷和约束后的模型如图2所示。然后将有限元模型使用求解器进行计算,并使用有限元软件的后处理模块查看分析结果,某极限工况下轮毂的应力云图如图4所示。根据德国劳埃德GL规范要求对计算的结果进行评估。
对轮毂体的疲劳分析,需根据德国劳埃德GL规范中规定的叶根坐标系下叶根位置单位载荷有限元分析中所获得的轮毂表面壳单元的应力应变结果,结合材料的疲劳寿命性能,应用疲劳理论,计算出零部件的疲劳寿命。所以对铸件疲劳寿命的计算分为两个部分,即单位载荷下铸件的有限元应力分析计算和根据铸件所承受载荷的疲劳载荷谱和单位载荷下应力值进行疲劳损伤值的计算。图5中的最小寿命为102.899=792.501,满足风机20年的设计寿命要求。
对轮毂体的疲劳分析,需根据有限元分析中所获得的应力应变结果,结合材料的疲劳寿命性能,应用疲劳理论,计算出零部件的疲劳寿命。所以对铸件疲劳寿命的计算分为两个部分,即单位载荷下铸件的有限元应力分析计算和根据铸件所承受载荷的疲劳载荷谱和单位载荷下应力值进行疲劳损伤值的计算。
名词解释:
1、德国劳埃德GL规范——又称德国劳埃德船级社规范,其中GL德国劳埃德船级社风力机设计标准,是由世界权威的叶片认证机构德国劳埃德船级社制定的专门关于风力设备的标准。
2、叶根坐标系——基于叶根处所建立的坐标系,见德国劳埃德GL规范中的规定。
圣维南原理——是弹性力学的基础性原理,是法国力学家A.J.C.B.de圣维南于1855年提出的。
Claims (4)
1.一种风力发电机组轮毂强度计算方法,主要用于大功率水平轴紧凑型风力发电机组,通过现有的有限元软件为平台实施建模,定义各个部件材料属性,进行计算,得出的应力云图,其特征在于,建模过程中主要包括以下步骤:
(1)叶片根部的一段、变浆轴承、变桨驱动、轮毂、主轴靠近轮毂的一段,分别建立有限元模型;
(2)在叶片根部创建节点并将其转换到德国劳埃德GL规范中规定的叶根坐标系下,然后刚性耦合到叶片上;
(3)变桨轴承的滚子采用只受压力不受拉力的单元进行模拟;
(4)叶片和变桨轴承、变桨轴承和轮毂、轮毂和主轴之间均定义为绑定接触。
2.根据权利要求1所述的一种风力发电机组轮毂强度计算方法,其特征在于:建模完成后,在轮毂体表面铺设一层很薄的壳单元。
3.根据权利要求1所述的一种风力发电机组轮毂强度计算方法,其特征在于:根使用载荷软件(GH bladed)输出GL叶根坐标系下叶根极限载荷,将载荷表中Mx、My、Mz和Myz绝对值最大工况分别施加在叶根中心,并对主轴端部施加全约束,得出四个极限工况下的应力云图。
4.根据权利要求1所述的一种风力发电机组轮毂强度计算方法,其特征是:对轮毂体的疲劳分析,需根据有限元分析中所获得的应力应变结果,结合材料的疲劳寿命性能,应用疲劳理论,计算出零部件的疲劳寿命;所以对铸件疲劳寿命的计算分为两个部分,即单位载荷下铸件的有限元应力分析计算和根据铸件所承受载荷的疲劳载荷谱和单位载荷下应力值进行疲劳损伤值的计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210262060.1A CN102831260B (zh) | 2012-07-27 | 2012-07-27 | 一种风力发电机组轮毂强度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210262060.1A CN102831260B (zh) | 2012-07-27 | 2012-07-27 | 一种风力发电机组轮毂强度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102831260A true CN102831260A (zh) | 2012-12-19 |
CN102831260B CN102831260B (zh) | 2015-02-04 |
Family
ID=47334395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210262060.1A Expired - Fee Related CN102831260B (zh) | 2012-07-27 | 2012-07-27 | 一种风力发电机组轮毂强度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102831260B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103366049A (zh) * | 2013-06-25 | 2013-10-23 | 沈阳华创风能有限公司 | 风力发电机组轮毂极限强度的计算方法 |
CN104573172A (zh) * | 2014-11-19 | 2015-04-29 | 新疆金风科技股份有限公司 | 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置 |
CN104699912A (zh) * | 2015-03-25 | 2015-06-10 | 北京博比风电科技有限公司 | 一种风力发电机组轮毂与主轴连接螺栓强度计算方法 |
CN105320785A (zh) * | 2014-07-09 | 2016-02-10 | 南车株洲电力机车研究所有限公司 | 一种用于风电机组载荷计算的工况全自动生成方法及装置 |
CN106014857A (zh) * | 2016-05-16 | 2016-10-12 | 国网冀北电力有限公司秦皇岛供电公司 | 抑制风电机组载荷的独立变桨控制方法及装置 |
CN107577846A (zh) * | 2017-08-07 | 2018-01-12 | 许继集团有限公司 | 一种风力发电机组机舱罩强度确定方法 |
CN109086465A (zh) * | 2017-06-13 | 2018-12-25 | 许继集团有限公司 | 一种风机轮毂与主轴连接螺栓的强度校核方法 |
CN109657320A (zh) * | 2018-12-11 | 2019-04-19 | 石河子大学 | 一种风电机组传动链建模与扭振特性分析方法 |
CN109726414A (zh) * | 2017-10-30 | 2019-05-07 | 北京万源工业有限公司 | 一种风力发电机轮毂结构疲劳计算方法及疲劳分析装置 |
CN109726410A (zh) * | 2017-10-31 | 2019-05-07 | 北京万源工业有限公司 | 风力机偏航连接螺栓强度校核的计算方法 |
CN109766634A (zh) * | 2019-01-11 | 2019-05-17 | 徐州徐工矿业机械有限公司 | 一种矿用大型铸钢件数字化正向研发方法 |
CN110598241A (zh) * | 2019-07-22 | 2019-12-20 | 许昌许继风电科技有限公司 | 一种风机轮毂强度校核方法 |
CN110879925A (zh) * | 2018-09-05 | 2020-03-13 | 江苏金风科技有限公司 | 风力发电机组的顶法兰的疲劳损伤确定方法和设备 |
CN116306179A (zh) * | 2023-05-22 | 2023-06-23 | 湖南云箭科技有限公司 | 机载外挂物起落冲击损伤评估用载荷筛选方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073761A (zh) * | 2010-12-31 | 2011-05-25 | 锋电能源技术有限公司 | 一种风力发电机组主轴强度计算方法 |
CN102508970A (zh) * | 2011-11-09 | 2012-06-20 | 保定天威风电科技有限公司 | 兆瓦级风力发电机组轮毂有限元建模方法 |
-
2012
- 2012-07-27 CN CN201210262060.1A patent/CN102831260B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073761A (zh) * | 2010-12-31 | 2011-05-25 | 锋电能源技术有限公司 | 一种风力发电机组主轴强度计算方法 |
CN102508970A (zh) * | 2011-11-09 | 2012-06-20 | 保定天威风电科技有限公司 | 兆瓦级风力发电机组轮毂有限元建模方法 |
Non-Patent Citations (5)
Title |
---|
吴树梁 等: "水平轴风电机组轮毂疲劳损伤计算方法", 《风机技术》 * |
孙传宗 等: "MW级风力发电机轮毂强度分析", 《沈阳工业大学学报》 * |
李秀珍 等: "风电轮毂变桨加强板的结构分析", 《风能》 * |
邓良 等: "大型风力发电机组轮毂强度数值分析", 《电气技术》 * |
马震岳 等: "水电站厂房机组动荷载施加方式研究", 《水力发电学报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103366049A (zh) * | 2013-06-25 | 2013-10-23 | 沈阳华创风能有限公司 | 风力发电机组轮毂极限强度的计算方法 |
CN105320785A (zh) * | 2014-07-09 | 2016-02-10 | 南车株洲电力机车研究所有限公司 | 一种用于风电机组载荷计算的工况全自动生成方法及装置 |
CN105320785B (zh) * | 2014-07-09 | 2018-05-29 | 南车株洲电力机车研究所有限公司 | 一种用于风电机组载荷计算的工况全自动生成方法及装置 |
CN104573172A (zh) * | 2014-11-19 | 2015-04-29 | 新疆金风科技股份有限公司 | 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置 |
CN104573172B (zh) * | 2014-11-19 | 2017-10-31 | 新疆金风科技股份有限公司 | 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置 |
CN104699912A (zh) * | 2015-03-25 | 2015-06-10 | 北京博比风电科技有限公司 | 一种风力发电机组轮毂与主轴连接螺栓强度计算方法 |
CN106014857B (zh) * | 2016-05-16 | 2019-11-15 | 国网冀北电力有限公司秦皇岛供电公司 | 抑制风电机组载荷的独立变桨控制方法及装置 |
CN106014857A (zh) * | 2016-05-16 | 2016-10-12 | 国网冀北电力有限公司秦皇岛供电公司 | 抑制风电机组载荷的独立变桨控制方法及装置 |
CN109086465A (zh) * | 2017-06-13 | 2018-12-25 | 许继集团有限公司 | 一种风机轮毂与主轴连接螺栓的强度校核方法 |
CN109086465B (zh) * | 2017-06-13 | 2022-12-13 | 许继集团有限公司 | 一种风机轮毂与主轴连接螺栓的强度校核方法 |
CN107577846A (zh) * | 2017-08-07 | 2018-01-12 | 许继集团有限公司 | 一种风力发电机组机舱罩强度确定方法 |
CN109726414A (zh) * | 2017-10-30 | 2019-05-07 | 北京万源工业有限公司 | 一种风力发电机轮毂结构疲劳计算方法及疲劳分析装置 |
CN109726410A (zh) * | 2017-10-31 | 2019-05-07 | 北京万源工业有限公司 | 风力机偏航连接螺栓强度校核的计算方法 |
CN109726410B (zh) * | 2017-10-31 | 2022-12-27 | 广东省风力发电有限公司 | 风力机偏航连接螺栓强度校核的计算方法 |
CN110879925A (zh) * | 2018-09-05 | 2020-03-13 | 江苏金风科技有限公司 | 风力发电机组的顶法兰的疲劳损伤确定方法和设备 |
CN110879925B (zh) * | 2018-09-05 | 2023-07-11 | 江苏金风科技有限公司 | 风力发电机组的顶法兰的疲劳损伤确定方法和设备 |
CN109657320A (zh) * | 2018-12-11 | 2019-04-19 | 石河子大学 | 一种风电机组传动链建模与扭振特性分析方法 |
CN109766634A (zh) * | 2019-01-11 | 2019-05-17 | 徐州徐工矿业机械有限公司 | 一种矿用大型铸钢件数字化正向研发方法 |
CN110598241A (zh) * | 2019-07-22 | 2019-12-20 | 许昌许继风电科技有限公司 | 一种风机轮毂强度校核方法 |
CN110598241B (zh) * | 2019-07-22 | 2022-12-13 | 许昌许继风电科技有限公司 | 一种风机轮毂强度校核方法 |
CN116306179A (zh) * | 2023-05-22 | 2023-06-23 | 湖南云箭科技有限公司 | 机载外挂物起落冲击损伤评估用载荷筛选方法及系统 |
CN116306179B (zh) * | 2023-05-22 | 2023-09-12 | 湖南云箭科技有限公司 | 机载外挂物起落冲击损伤评估用载荷筛选方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102831260B (zh) | 2015-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102831260B (zh) | 一种风力发电机组轮毂强度计算方法 | |
Chehouri et al. | Review of performance optimization techniques applied to wind turbines | |
Grogan et al. | Design of composite tidal turbine blades | |
Gasch et al. | Wind power plants: fundamentals, design, construction and operation | |
CN105160067A (zh) | 风机变桨轴承与轮毂和叶片连接螺栓强度性能的校核方法 | |
Ageze et al. | Wind turbine aeroelastic modeling: basics and cutting edge trends | |
CN103063425A (zh) | 风力机叶片疲劳寿命实时监测方法 | |
CN106286151B (zh) | 一种风电机组低速轴扭转载荷监测方法及载荷分析方法 | |
Keerthana et al. | Aerodynamic analysis of a small horizontal axis wind turbine using CFD | |
Mendoza et al. | Aero-structural rapid screening of new design concepts for offshore wind turbines | |
CN102508970A (zh) | 兆瓦级风力发电机组轮毂有限元建模方法 | |
Sompong et al. | Effects of bend-twist coupling deformation on the aerodynamic performance of a wind turbine blade | |
Maalawi | Special issues on design optimization of wind turbine structures | |
Ferrari | Development of an aeroelastic simulation for the analysis of vertical-axis wind turbines | |
Riziotis et al. | Advanced aeroelastic modelling of complete wind turbine configurations in view of assessing stability characteristics | |
Kesikbaş | Investigations of upscaling effects for aerodynamic design of large wind turbine rotors by using bem theory and optimization | |
Roth-Johnson | Aero-Structural Design Investigations for Biplane Wind Turbine Blades | |
CN105320793A (zh) | 一种从动力学角度评估风电机组变浆控制模型的方法 | |
Ufaira et al. | Numerical Study of Deflection and Stress Distribution on Composite Box Spar Structure–Application In Wind Turbine Blade | |
Goroshko et al. | Simulation of the dynamic response in a storm scenario | |
Choudhury et al. | Effect of orthotropy ratio of the shear web on the aero-elasticity and torque generation of a hybrid wind turbine blade | |
Sargın | Comparative study of transient and quasi-steady aeroelastic analysis of composite wind turbine blade in steady wind conditions | |
Silva Junior | Stress analysis on a thin-walled composite blade of a large wind turbine | |
de Saint-Exupéry | Design Principles, Scaled Design, and Optimization | |
Nigam et al. | Review on design of horizontal axis wind turbine blades |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150204 Termination date: 20190727 |