[发明专利]兆瓦级风力发电机组轮毂有限元建模方法

说明书

技术领域

本发明涉及一种兆瓦级风力发电机组轮毂有限元建模方法，属于风力发电技术领域。

背景技术

目前，在风电技术领域，对于风电机组轮毂，主要采用的结构形式是“球形”，由于其结构的复杂性，不能像叶片、主轴、塔筒等细长零部件一样进行等效分析，它的强度分析（极限强度和疲劳强度）工作主要是采用有限元方法进行，如说明书附图1所示。在模型中，主要包括了轮毂主体、变桨轴承、叶片假体、部分主轴几个零部件，建立叶片假体和部分主轴的的依据是圣维南原理，使轮毂的连接部位的应力也接近真实。在主轴末端（1位置）施加全位移约束，在叶片根部（2位置）施加载荷，载荷通过刚性区域传递到叶片假体顶部，从而对轮毂主体施加作用力。风机在实际运行中，轮毂要承受三个叶片载荷的同时作用，且在同一时刻，三个叶片的载荷是不一样的。背景技术对于轮毂的极限强度分析，采用的方法是：当某一个叶片的某一个分量载荷达到最大值时，连同其它叶片位置的其余载荷分量一同施加到模型中，计算完成后去检查轮毂应力情况，根据相关行业标准去检查其是否满足极限强度要求。但是，这样做的缺点是：它假设当某一个叶片的某一个载荷分量达到最大时，轮毂承受了最大的极限载荷。但是实际情况是：当某一个叶片的某一个分量载荷达到最大值时，轮毂不一定承受最大的载荷。风机运行中经历的载荷工况很复杂，可能在某一个工况下，三个叶片载荷都很大，但是没有一个叶片的某一个载荷分量达到最大值。对于疲劳分析，背景技术一般的做法是：在每一个叶片根部施加单位载荷，从而计算轮毂在单位载荷作用下的单位应力，然后结合载荷的时间序列进行分析，根据Miner线性累计损伤理论进行疲劳分析。这样做，需要在每一个叶片的根部施加六个载荷分量（Mx、My、Mz、Fz、Fy、Fz）进行计算，对三个叶片风机，就要对模型进行18次的加载计算，一般来说，轮毂的有限元模型规模很大，计算每个载荷都需要花费一定的时间，这样做起来很费时间，工作效率低。

发明内容

本发明的目的是提供一种兆瓦级风力发电机组轮毂有限元建模方法，对所有极限载荷工况进行计算，并且减少了有限元模型的计算次数。对于疲劳强度分析，只需要在一个叶片位置施加六个载荷分量，计算六次即可，极大的提高了工作效率，解决背景技术存在的上述问题。

本发明的技术方案是：兆瓦级风力发电机组轮毂有限元建模方法，包含如下步骤：①根据球形轮毂的对称性，取其1/3的模型进行建模，形成1/3的有限元模型，然后旋转形成整体的轮毂有限元模型；②对于原1/3模型上的每一单元，在120度和240度方向上都有单元和其对应，即都能找到由这个单元拷贝形成的在120度和240度位置上的单元；③利用球形轮毂的旋转对称性来进行计算和强度校核，旋转对称性是指：在其中一个叶片加载和在另一个叶片加载形成的轮毂应力云图，两者旋转120度或者240度，是能够完全重合的。

对于轮毂极限强度分析，根据线性叠加原理，在原1/3模型对应的叶片根部分别施加六个单位载荷分量（对于Mx、My、Mz取值为1NM，对于Fx、Fy、Fz取值为1N）进行计算。以一个有限元单元为例，在某一个单位载荷下，可以提取此单元的应力，同时，也可以得到120度方向和240度方向位置两个单元的应力（此两个单元是由原单元拷贝生成），利用球形轮毂的旋转对称性，就可以得到在其他两个叶片根部分别施加同样单位载荷时在此单元的应力贡献。这样，有了每一个叶片根部的单位载荷对此单元的应力贡献，就可以得到在三个叶片根部都存在载荷时，在此单元的应力值（线性叠加原理）。据此，就可以计算每一个极限载荷工况下在此单元的应力。

对于疲劳分析，同样根据旋转对称性，就可以只计算在1/3模型上叶片根部的六个载荷分量下的应力，其余叶片根部施加载荷时在同一个位置贡献的应力可以根据旋转对称性得到。这样，相当于减小了2/3的计算工作量，极大的提高了工作效率。

本发明的积极效果：对于轮毂的极限强度分析，可以考虑到所有的极限载荷工况，而不仅仅是局限于某一个载荷分量达到最大时的情况，提高了计算的可靠性、全面性，同时，减小了计算量。对于疲劳强度分析，极大的减小了计算的工作量，提高了工作效率。

附图说明

图1是用于分析的整体有限元模型；

图2是1/3的轮毂有限元模型；

图3是经过旋转拷贝以后生成的整体轮毂有限元模型；

图4是旋转对称选取单元应力说明；

图中：1.叶片1；2.叶片2；3.叶片3；4.有限元单元一；5.有限元单元二；6.有限元单元三。