风力机轴承的有限元分析Word格式文档下载.docx

风力机轴承的有限元分析Word格式文档下载.docx

《风力机轴承的有限元分析Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《风力机轴承的有限元分析Word格式文档下载.docx（26页珍藏版）》请在冰豆网上搜索。

1引言

1.1课题研究背景

现代社会和经济的发展越来越离不开能源。

从能源发展战略需求来看，因为化石燃料是有限性的、不可再生的，并且使用化石燃料生产会产生环境污染，所以要寻找一条新的、可持续发展的能源道路，以满足生产和生活的需要。

开发利用洁净的、无污染的、可再生的能源越来越受到科研人员的重视。

在新型能源里，利用风能发电是开发方案之一。

目前，尽管世界各国的风力发电量还不到世界总耗电量的2%，但随着风力发电装机容量的快速发展和风力发电机技术的成熟和不断完善，风力发电必将成为世界各国更加重视和重点开发的能源之一。

从全球范围来看，风力发电己经取得了很大、很快的发展。

根据相关部门的统计，全球的风力机容量达到了8000MW。

在西方的一些欧美国家，如美国、德国、丹麦、西班牙等把风力发电的研究创新、生产与应用列为重点研究项目，相继加大了人力及资金的投入。

现在风力发电机的研究，属于交叉学科研究，涉及到空气动力学、电力电子技术、机械设计及电子通信技术等方面，研究目的是引入新技术、新型材料、新型电机以优化结构、提高性能和开发新类型发电机，发展建立一个风能资源测量和评估的计算机仿真系统，并利用自动化控制技术，提高单个和多个风力发电机的效率和可靠性。

目前，德国是世界上最大的风力发电国家，从上个世纪九十年代，其风力发电量一直持续增长，装机容量已达到3000MW。

并且德国的风力机研制能力强，技术水平处于世界前列，其国内有两家世界上著名的风机制造商。

美国在加利福尼亚洲建立了大量的风力发电机，其发电量大约占全美发电量的1%。

棕榈泉风电场大约有4000台风机，是美国最大的风力发电厂，装机容量为50万kW,发电量和一个大型火力发电厂相当。

丹麦是第一个使用风力发电的国家，风力发电的规模在世界排第三，其风力发电量已占全国总发电量的3%，总的装机容量为1450MW。

丹麦风力发电机研究水平和制造能力位于世界前列，丹麦有六家世界著名的风力机制造商。

此外，西班牙、印度、意大利等国家的风力发电规模也在世界前列。

我国幅员辽阔，风能资源丰富，发展风力发电机前景广阔、前途光明。

据不完全统计，可利用的风能达25.3亿MW,特别是以东北、华北、西北构成的大面积陆地型风能资源丰富区，应拓宽风力发电机的应用市场。

我国在20世纪80年代初开始风力发电技术的研究，并且通过国内科研人员的努力，研发了小型风力发电机（额定容量为100W〜10kW），并已经有一定规模的生产，取得了广泛的应用。

在“九五”计划中，国家大力投资进行风电的技术研究和建设、政策的规划

。

在国家的大力支持下，经过科研人员的努力，相继研制出了额定功率为250KW、300KW、600KW的风电机组，并逐步实现商业量产化。

到目前为止，全国共建立了近50个风电场，总装机容量达46万KW。

新疆达坂城风电场是目前国内最大的风电场，总装机容量为575MW;

广东南澳风电场的装机容量为42.8MW。

虽然我国在风力发电的研制上取得了一定的成绩，但是从发展速度和制造生产上看，我国的风电机和欧美国家相比还有很大的差距，所以风力机存在一定的质量问题，再加上风力机的运行环境相对比较恶劣，常常发生零部件失效，如风力机中的叶片折断、齿轮箱的轮齿折断、风力机轴承的损坏等。

1．2风电轴承简介

风力发电机用轴承主要包括：

风机主轴轴承、偏航轴承、变浆轴承、变速箱轴承、发机轴承等，轴承的类型主要有四点接触球轴承、交叉滚子轴承、圆柱滚子轴承、调心滚子轴承、深沟球轴承、双列圆锥滚子轴承等。

风能工业属于高技术行业，它的发展存在一定的技术风险性。

目前，国内很多的先进技术都是由国外引进的。

德国等欧洲国家之所以能位居世界领先地位，在于它拥有世界上领先的风机开发技术和配套产品的能力。

但反观国内，虽然我们的配套企业发展得很快，但还都是刚刚起步。

洛轴和瓦轴作为国内轴承行业的领军企业，致力于轴承的研究和生产制造，对风力发电机的偏航和变浆轴承有了一定的生产能力，但是从总体上讲，国内的风电轴承企业无论是生产能力还是所生产轴承的性能方面都较国外有很大的差距

轴承在风机上很多部位都有应用，以下为风力发电机中用到的大型轴承：

·

主轴轴承：

根据不同的设计类型和使用工况，风机转子主轴可以使用一个轴承或者两个轴承，多采用自动调心滚子轴承或圆锥滚子轴承，其中调心滚子轴承应用较多，主轴的偏心角度一般为

变浆轴承：

变浆轴承用于连接轮毂和叶片，并随着风向的变化来改变叶片的角度，其结构为双列四点接触球轴承，有带内齿和无齿两种结构。

转动范围为

，正常在

范围内，受倾覆力矩大。

偏航轴承：

偏航轴承是风机追踪风向调整迎风面的保证，转动范围为，

在

范围上转动频率最高，上面机舱重量可达近百吨，且要承受倾覆力矩，高刚度的轴承结构和设计良好的密封性能，良好的防腐蚀性能是关键，轴承结构采用四点接触球轴承。

转子主轴轴承（图1．1）支撑叶轮和主轴并将风能转矩传输至齿轮箱，因此主轴轴承承受的载荷非常大，并且因风经常变化，作用在主轴上的负载和转速变动也很大。

在起动风速以下的情况下，处于空转状态下的转子主轴轴承进行低速、轻负载旋转。

由于上述原因及主轴较长、自重大、容易变形等原因，风机主轴轴承多选用调心滚子轴承。

图1.1

由于风机大部分安装在偏远地区，工作环境恶劣，而且维修极不方便，因此，以上所述轴承都必须具有高可靠性，一般风机寿命为20年，因此主轴轴承寿命要大于20年。

风机轴承国产化可提高国内轴承工业的设计应用水平，缩小与国外先进水平的差距，促进国内轴承工业的发展和技术进步，另一方面，可以降低风电成本，加快我国新能源和可再生能源的发展。

1．3自动调心滚子轴承及风力发电机主轴轴承研究

由于风力发电的兴起，有关风力发电方面的研究也日益增多，但是大部分研究主要针对叶片（控制、强度分析等）

，变电系统的控制传输以及发电机相关方面的研究等

，对风力发电机主轴轴承的研究很少。

有关自动调心滚子轴承的研究有不少，风力发电机主轴轴承的研究可以参考这些研究成果。

由于双馈式风力发电机主轴轴承结构大多采用双列自动调心滚子轴承结构，这种结构综合了滚子轴承的高承载能力和调心球轴承的可调心性能，对于该类型的轴承国内外研究较多。

Klecknert

对自动调心滚子轴承进行了分析和软件设计，重点考虑了在给定温度下轴承运行时的力学性能，所编制软件SPHERBEA（SPHERERICALBEARINGANALYSIS）可以进行轴承在轴向、径向和联合载荷下的分析，分析还可以考虑弹流润滑、滚子间隙、滚子速度等。

Shelofast

分析了施加外载荷时的非理想自动调心滚子轴承接触问题，载荷施加在滚动体上，接触位移可由平衡公式得出，滚子用多个薄片代替来分析每一片的接触。

并与APMWinBear专用

软件的模拟做了对比，可确定轴承出现故障的位置。

Timken公司的Zantopulos，Harry,Russell，Joel对双列的圆锥滚子轴承和自动调心滚子轴承做了实验，采用的轴承直径均为110mm，施加22％的额定动态径向载荷，结果显示圆锥滚子轴承的寿命要比自动调心滚子轴承的寿命大的多，但这跟IS0281标准不一致。

日本NSK公司的Ueda，Tohru叙述了自动调心滚子轴承表面失效的机理，指出由于滚动体和滚道之间的摩擦产生的切应力严重影响了轴承的寿命。

日本NTN技术研发中心的Akamatsu,Yoshinobut在前人研究径向载荷和联合载荷下的自动调心滚子轴承数值分析以及实验的基础上，提出了一种类似静态分析方法，对对称型和非对称型自动调心滚子轴承做了研究，并指出在润滑良好的情况下对称型自动调心滚子轴承的扭矩大于非对称型，但是在润滑不良的情况下却比非对称型的低。

由于自动调心滚子轴承可以自动调心并具有高的承载能力因此得到了广泛的应用，和深沟球轴承相比，自动调心滚子轴承的动态特性要复杂的多，接触副上的受力也不相同，这将导致轴向力在滚道上的不平衡和角度的变化。

Cao建立了自动调心滚子轴承的动态模型，考虑了滚道的自由度，这种模型可以定量地分析轴的偏斜量，为要求高的轴承安装提供了指导。

在文献中Cao指出了当载荷较小的时候自动调心滚子轴承中滚子和滚道为点接触，当载荷作用在轴承上接触区域长度超过点接触的极限时为线接触，该极限为接触区域长轴的长度大于滚子的有效长度。

Royston和Basdgan

为评定自动调心滚子轴承刚度系数建立了理论模型，但该模型只能进行定性的分析；

在考虑内部几何形状和预加载荷的情况下应用矢量和矩阵的方法描述了轴承圈之间的全部弹性变形，提供了一个可选择的方法进行内、外圈刚度系数的模拟，Bercea的模型主要是静态的模型，没有考虑滚子和滚道在动态非线性状况下自动调心滚子轴承的性能。

对于风力发电机主轴轴承的研究主要有：

Wildinson设计了检测平台，可以检测风力发电机主轴轴承及齿轮箱等，通过对扭矩、速度、振动等信号的检测可以预测轴承及其他部件的工作状况。

作为风力发电机的重要组成部件，轴承经常受到电击的危险，Paolone，Napolitano等提出了数值和解析的模型来计算风机主轴轴承的电阻抗，并设计了按比例缩小的试验平台进行测试，得到良好的效果。

该方法为大型风力发电机主轴轴承试验提供了参考。

陈长征、孙长城

等人提出了风机轴承可能出现的故障，并引用应力波技术来检测低速轴承的故障。

首先，他们以SKF23068CC／W33为例，进行了三维建模，模拟了滚道内圈和外圈缺陷下应力、应变和接触区载荷分布情况，并和无故障情况下的轴承进行了比较。

通过对应力波的分析得出结论：

小波分析理论可以引用到低速轴承的诊断。

自从自调心滚子轴承出现以来一直没有人提出修形的方法，以上文献中对于自调心滚子轴承的研究都没有涉及到滚子或滚道的修形设计，由滚子类轴承凸度设计知识知

，对轴承滚子或者滚道进行适当的设计应该能效降低最大接触应力，故有必要研究修改滚子、滚道的密合度以及采用变曲率修形的方法来降低自调心滚子轴承最大接触应力，提高轴承寿命的设计方法。

1.3本文主要研究内容

2轴承分析应用概述

2.1有限元分析技术基础概念

有限元分析（FEA）是使用有限元方法来分析静态或动态的物理物体或物理系统。

其基本概念是用较简单的问题代替复杂问题后再求解。

在这种方法中一个物体或系统被分解为由多个相互联结的、简单、独立的点组成的几何模型。

有限元法最初的思想是把一个大的结构划分为有限个称为单元的小区域，在每一个小区域里，假定结构的变形和应力都是简单的，小区域内的变形和应力都容易通过计算机求解出来，进而可以获得整个结构的变形和应力。

它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的（较简单的）近似解，然后推导求解这个域总的满足条件（如结构的平衡条件），从而得到问题的解。

有限元法求解一个问题是要求解联立代数方程组，而不是解微分方程。

求解得到的这些数值解给出的是连续体中多个离散点的未知量的近似值。

用有限元法可以解决结构分析、热传导、流体力学的流场分析、电磁场分析等问题[。

有限元法分析计算的思路和做法可归纳如下。

（1）物体离散化

将某个工程结构离散为由各种单元组成的计算模型，这一步称物体离散化。

作单元剖分离散后，单元于单元之间利用单元的节点相互连接起来；

单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。

所以有限元中分析的结构已不是原有的物体或结构物,而是同材料的由众多单元以一定方式连接成的离散物体。

这样，用有限元分析计算所获得的结果只是近似的。

如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。

（2）单元特性分析

选择位移模式。

在有限单元法中，选择节点位移作为基本未知量时称为位移法；

选择节点力作为基本未知量时称为力法；

取一部分节点力和一部分节点位移作为基本未知量时称为混合法。

位移法便于实现计算程序化，因此，在有限单元法中位移法应用最多。

当采用位移法时，物体或结构物离散化之后，就可把单元中的一些物理量如位移，应变和应力等由节点位移来表示。

这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。

对于一般情况，我们就用坐标变量的简单函数来表示。

分析单元的力学性质。

由已知的单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。

此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。

计算等效节点力。

物体离散化后，假设力是通过节点从一个单元传递到另一个单元。

但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。

因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上的力。

（3）单元组集

利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程，求解未知节点位移有限元方程式。

这里，可以根据方程组的具体特点来选择合适的计算方法。

对于不同物理性质和数学模型的问题，有限元求解法的基本步骤是相同的，只是具体公式推导和运算求解不同。

以下是有限元求解问题的基本步骤。

第一步：

问题及求解域定义。

根据实际问题近似确定求解域的物理性质和几何区域。

第二步：

求解域离散化。

将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域，习惯上称为有限元网络划分。

显然单元越小（网络越细）则离散域的近似程度越好，计算结果也越精确，但计算量及误差都将增大，因此求解域的离散化是有限元法的核心技术之一。

第三步：

确定状态变量及控制方法。

一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示，为适合有限元求解，通常将微分方程化为等价的泛函形式。

第四步：

单元推导。

对单元构造一个适合的近似解，即推导有限单元的列式,其中包括选择合理的单元坐标系，建立单元函数，以某种方法给出单元各状态变量的离散关系，从而形成单元矩阵（结构力学中称刚度阵或柔度阵）。

为保证问题求解的收敛性，单元推导有许多原则要遵循。

对工程应用而言，重要的是应注意每一种单元的解题性能与约束。

例如，单元形状应以规则为好，畸形时不仅精度低，而且有缺秩的危险，将导致无法求解。

第五步：

总装求解。

将单元总装形成离散域的总矩阵方程（联合方程组），反映对近似求解域的离散域的要求，即单元函数的连续性要满足一定的连续条件。

总装是在相邻单元节点进行，状态变量及其导数连续性建立在节点处。

第六步：

联立方程组求解和结果解释。

有限元法最终导致联立方程组。

联立方程组的求解可用直接法、迭代法和随机法。

求解结果是单元节点处状态变量的近似值。

对于计算结果的质量，将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算。

总而言之，有限元分析可分成三个阶段，前处理、处理和后处理。

前处理是建立有限元模型，完成单元网格划分；

后处理则是采集处理分析结果，使用户能简便提取信息，了解计算结果。

2.2ANSYS软件简介

ANSYS有限元软件是一个多用途的有限元法计算机设计程序，可以用来求解结构、流体、电力、电磁场及碰撞等问题。

因此它可应用于航空航天、汽车工业、生物医学、桥梁、建筑、电子产品、重型机械、微机电系统、运动器械等工业领域。

鉴于ANSYS软件在本文的计算分析中的地位十分重要，下面将对其进行简单的介绍。

ANSYS软件主要包括三个部分：

前处理模块，求解计算模块和后处理模块。

前处理模块主要包括参数定义和建立有限元模型。

用户可以方便地构造有限元模型。

该模块用于定义求解所需的数据，用户可选择坐标系统、单元类型、定义实常数和材料特性、建立实体模型并对其进行网格剖分、控制节点单元，以及定义藕合和约束方程等。

通过运行一个系统模块，可以预测求解过程所需的硬件大小及内存。

求解计算模块对所建立的有限元模型进行力学分析和有限元求解。

前处理阶段完成建模后，用户在求解阶段通过求解器获得分析结果。

在该模块中，用户可以定义分析类型和分析选项、施加载荷和载荷步选项，然后开始有限元求解。

主要包括以下子模块：

结构静力学分析、结构动力学、结构非线性、动力学分析、热分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力。

后处理模块可用于查看整个模型或选定的部分模型在某一子步或时间步的求解结果。

利用该模块可以获得各种应力场、应变场及温度场等的等值线图形显示、变形形状显示以及检査和解释分析的结果列表。

可将计算结果以彩色等值线、矢量、粒子流迹、透明及半透明等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。

该模块可以通过便捷的用户界面获得求解过程的计算结果并对这些结果进行运算，这些结果可能包括位移、温度、应力、应变、速度及热流体，输出形式有图形显示和数据列表两种。

在交互式的处理过程中，图形可联机输出到显示设备上或脱机输出到绘图仪上

由于后处理阶段完全同ANSYS前处理和求解阶段集成在一起，故求解结果已保存到数据库中且能立即查看。

ANSYS软件提供了一百种以上的单元类型，可供工程中的各种结构和材料来选择。

ANSYS软件含有多种有限元分析的能力，包括从简单线性静态分析到复杂非线性动态分析。

ANSYS的技术特点如下：

可实现多场及多场耦合功能；

实现前后处理、分析求解及多场分析统一数据库的大型FEA软件；

具有流场优化功能的CFD软件；

融前后处理与分析求解于一体；

强大的非线性分析功能；

快速求解器；

可与大多数的CAD软件集成并有接口；

具有智能网格划分;

具有多层次多框架的产品系列。

ANSYS软件的强大功能，为其在工程领域的广泛应用提供了良好的基础，也为本文中齿轮的动态分析提供了便利。

3接触分析

3.1概述

接触问题是一种高度非线性行为，需要较多的计算机资源。

为了进行切实有效的计算，理解问题的物理特性和建立合理的模型是很重要的。

接触问题存在两个较大的难点：

其一，在用户求解问题之前，用户通常不知道接触区域。

随载荷、材料、边界条件和其它因素的不同，表面之间可以接触或者分开，这往往在很大程度上是难以预料的，并且还可能是突然变化的。

其二，大多数的接触问题需要考虑摩擦作用，有几种摩擦定律和模型可供挑选，它们都是非线性的。

摩擦效应可能是无序的，所以摩擦使问题的收敛性成为一个难点。

注意--如果在模型中，不考虑摩擦，且物体之间的总是保持接触，则可以应用约束方程或自由度藕合来代替接触。

约束方程仅在小应变分析（NLGEOM，off）中可用。

见《ANSYSModelingandMeshingGuide》中的§

12，CouplingandConstraintEquations。

除了上面两个难点外，许多接触问题还必须涉及到多物理场影响，如接触区域的热传导、电流等。

3.1.1显式动态接触分析能力

除了本章讨论的隐式接触分析外，ANSYS还在ANSYS/LS-DYNA中提供了显式接触分析功能。

显式接触分析对于短时间接触-碰撞问题比较理想。

关于ANSYS/LS-DYNA的更多的信息参见《ANSYS/LS-DYNAUser"

sGuide》。

3.2一般接触分类

接触问题分为两