《基于ANSYS的垂直轴风力机塔架的力学性能分析及优化》.doc

基于ANSYS的垂直轴风力机塔架的力学分析及结构优化 （山东科技大学 机电学院，山东 青岛 266510）摘要: 基于ANSYS对大型垂直轴风力机塔架进行静态学分析、动力学分析及结构优化。首先根据风力机的功率确定塔架的大体尺寸，然后对塔架进行静力分析，并对塔架承受最大载荷时的最大应变进行强度校核；动力学分析包括有预应力的模态分析、地震谱响应分析。ANSYS优化通过改变塔架的厚度使其固有频率增大，从而避免因塔架在发电机的激励下共振产生的破坏。通过对塔架的力学分析确保了塔架性能的可靠性，为塔架的安全性提供了理论依据，并且在满足力学性能的前提下使塔架的质量最轻，降低了塔架的制造成本。关键词: 垂直轴风力机塔架；静力学分析；动力学分析；地震谱响应分析；结构优化；ANSYS中图分类号:TK83 文献标志码: A 文章编号:1672-3767(2011)07-0000-07Mechanical Performance Analysis And Structure Optimization Of Tower Of Vertical Axis Wind TurbineBased On ANSYS （College of Mechanical and Electronic Engineering， Shan Dong University of Scienceand Technology，Qing Dao 266510, China）Abstract: Statics analysis, dynamics analysis and structural optimization of large tower of vertical axis wind turbine were presented based on ANSYS. First of all, according to the power of the wind generators the size of the tower was determined. And then statics analysis for the tower was carried out. And the strength check for the maximum strain of the tower was discussed under the maximum load. Dynamics analysis including prestressed modal analysis, seismic response spectrum analysis was performed. Using ANSYS optimization by changing the thickness of the tower of the inherent frequency increases, so as to avoid the generator for tower under the excitation of resonance produced damage. Through to the mechanical analysis of the tower to ensure the reliability of the tower of the performance for the safety of the tower, to provide the theoretical basis for the mechanical properties, and make the quality of the tower of the lightest in meet the premise of mechanical properties, reduce the tower of manufacturing cost.Keywords: tower of vertical axis wind turbine; statics analysis; dynamics analysis; seismseismic response spectrum analysis; structural optimization; ANSYS风能资源主要集中在50米以上的高空，是低空资源的2倍左右，大容量的风力机机组便于大规模风能的开发利用，因此大功率的风机具有很好的发展前景。在风力发电机中，塔架占成本的25%左右1。塔架是风力机的承载部分，塔架的设计是风力发电机组不可缺少的一部分，尤其对于大功率的风力机组2。结构振动导致的材料疲劳通常是系统损坏的主要原因之一，振动还会产生噪音，影响周围环境，所以塔架的力学分析对于风力机的安全可靠性能具有重要的参考价值。目前，对风力机的研究大都集中在叶片、发动机性能等的研究上，对于塔架的研究很少，尤其是针对大塔架的研究更为少见3，垂直轴风力塔架的分析在国内基本是空白。本研究是在前人分析水平轴风力机的基础上，根据垂直轴风力机的实际受力对塔架进行加力并分析，为大型垂直轴的风力机塔架的制造提供了理论依据。本文首先对垂直轴风力机塔架的模型尺寸的计算与建模，接着对模型加载载荷进行静态的受力分析，在受力的情况下对模型进行有预应力的模态分析，得到塔架的前6阶模态及振型。在模态的基础上，加载模拟的地震频率波动进行谱分析，最后针对塔架模型进行结构的优化，优化的参数是塔架底的壁厚和顶壁厚以及塔底和塔顶的外半径，目标函数是塔架的第一阶频率、质量、最大应力和最大位移。1 垂直轴风力机塔架的计算模型大型风力机组大都采用圆锥形塔柱，本研究分析的是1MW的垂直轴风力机的塔柱，根据高耸机构设计规范GB50135-2006，选择塔柱的尺寸如图1所示：塔高H=70m，底外径D=4.5m，顶外径D=3m，底端壁厚B=30mm，顶端壁厚B=20mm，底座长L=0.5m，材料为Q345钢，材料弹性模量E=2.061011，泊松比=0.3，密度=7.850103kg/m3,许用强度=237MPa，许用抗拉强度=425MPa。此塔架为变截面的锥形圆锥体，由于底端是固定的，可以看成悬臂梁结构，塔架的主要载荷有风的水平推力F，机舱、叶片和轮毂的重力G，塔架的自重G，轮毂的扭矩M等，作用在塔身上的载荷共5种，受力简图如图2所示。1.1 塔架产生的主要载荷1）风轮的水平推力： (1) 式中空气密度（1.293kg/m3）； 图1 塔架结构模型V风速（定为30m/s）； Fig.1 The structure model of towerA风轮的扫风面积（取1700m2）；C推力系数（一般处切入风速时取1左右，切出风速取0.5左右，根据贝茨公式计算，C=8/92）机舱、叶片和轮毂的力： G=mg (2) 式中m为机舱、叶片和轮毂的总重量（质量取39000kg）；g为重力加速度（取9.8m/s2）3）塔架自身的重量：G4）轮毂扭矩：M=P/n (3) 式中P为发电机的输出功率 n为风轮的转速（r/s） 图2 塔架的受力简图 Fig.2 The force diagram of tower5）作用在塔架上的风压： (4)式中为空气的密度；Z为塔架的高度；V是高度为Z处的风速;Z为高度；H为高度；V是高度为10m处的风速；a是地面粗糙度和地面风的切变系数。有公式（1）可得：F=879240(N)；有公式（2）可得：G=382200(N)；把P=1106瓦，n=1/3(r/s)带入公式（4）可得：M=3106 (NM)2 用ANSYS对塔架进行静力分析2.1 建模、加载与分析建模：首先根据上面塔架模型数据在DM中画出塔架的二维模型然后以中心轴为旋转轴进行旋转就得到了塔架的模型，建模完成。加载与分析：进入静力分析模块对模型进行网格划分然后施加边界条件和载荷，步骤如下：1）固定塔架的底端；2）在塔架顶部加载力F和机舱、叶片和轮毂的重力G以及轮毂扭矩M；3）在塔架上加载塔架的重力G和风压P，结果设置为查看总体的应力和总体的应变，最后进行求解。由求解的结果图3可以得到塔架的最大位移为0.69965m，由图4得到塔架的最大应力 图3 塔架的位移/m为155.18MPa Fig. 3 The displacement 0f tower /m2.2 结果分析与讨论根据材料的刚度要求塔架的最大变形在高度H=70m的1%一下即：( m ) （5）计算得出的0.69965(m)0.7(m),不超出最大变形的要求，计算得到得最大应力为155.18MPa=237MPa,强度符合要求，所以塔架结构和静力学性能满足正常工作的要求。3 塔架的模态分析 图4 塔架的应力/Pa塔架的模态分析对塔架的稳定性非常重要， Fig.4 The stress of tower /Pa 模态分析主要是计算塔架的固有频率和阵型，它也是动力学分析的基础，同时也是地震谱响应分析的前序和结构参数优化的参数标准，如果塔架受力，就应该进行有预应力的模态分析，分析的结果才有意义。塔架的前几阶模态对稳定性的影响比较大，因此分析塔架的前六阶模态。3.1 塔架的有限元动力学方程动力学分析同静力学分析都是把塔架离散成有限个单元体，动力学分析要考虑单元特性的同时还考虑单元的惯性力和阻尼力等因素。动力学方程： （6）式中：M为塔架的质量矩阵；C为塔架的阻尼矩阵；K为总体刚度矩阵；X为的节点位移；f为塔架受的外力3.2 ANSYS对塔架模态分析的结果从图5到图10可是得到塔架的一到六阶振型、各阶振型下产生的最大位移和固有频率。如在图5中能得到塔架的振型是一阶的，最大的位移在塔架的顶端为5.676910-3m，一阶固有频率为1.0363Hz。 图5 塔架的第一阶模态振型Fig.5 The first order modal shape of tower图7 塔架的第三阶模态振型Fig .7 the third order modal shape of tower图9 塔架的第五阶模态振型Fig .9 the fifth order modal shape of tower图6 塔架的第二阶模态振型Fig.6 The second order modal shape of tower图8 塔架的第四阶模态振型Fig .8 the fourth order modal shape of tower图10 塔架的第六阶模态振型Fig .10 the sixth order modal shape of tower表1 塔架的前六阶固有频率Tab.1 The the first six order natural frequency of towerMode123456Freq(Hz)1.03631.03974.95384.969912.42712.4633.3 结果分析与讨论根据公式可知，当风力发电机以20r/min的转速运行时，风轮旋转产生的频率为0.3333Hz，风叶的个数是3，f=1Hz，由表1可知，这个频率十分接近塔架的第一、二阶固有频率，频率要满足公式( ff ) /f 10%，才能不引起塔架的共振，而f=1.0363，( ff )/ f=0.03630.1,所以塔架在发电机的激励下是产生共振的，共振对结果的破坏非常大，因此要对塔架的结构进行优化，使它的前两阶频率得到提高，降低共振引起的破坏。4 塔架的地震谱响应分析表2 地震模拟数据Tab.2 The earthquake simulation data序号123456789101112131415频率0.100.110.130.140.170.200.250.330.500.671.001.111.251.431.67加速度0.0020.0030.0030.0050.0060.0060.0100.0210.0320.0470.0700.0880.1050.1100.130序号1617181920212223242526272829频率2.002.503.334.005.006.6710.0011.1112.5014.2916.6720.0025.0050.00加速度0.1500.2000.2550.2650.2550.2000.1650.1530.1400.1310.1210.1110.1000.100地震谱响应就是把模拟地震的波动加载到模型得到的响应。地震对于大型塔架有着毁灭行的破坏，进行地震谱响应分析是十分必要的。表2为中国地震多发地带的地震模拟数据。结果与分析：由图11可得到地震谱响应的最大应力值为7.4838Mpa，远远小于材料的许用应力=237MPa，塔架安全可靠。5 塔架结构的优化分析从结构的模态分析可以看出，风力机结构的安全性存在巨大的隐患，需要对塔架进行结构优化，优化应同时满足结构的安全性和经济性。提高钢塔柱的固有频率可以有效降低结构共振的几率。K=P/ （7） 图11 谱响应的塔架应力 f= （8） Fig.11 The tower stress of Spectrum responseK为刚度矩阵；P是作用在塔架上的恒力；为由于力产生的变形；f为固有频率；M为质量矩阵。由公式（8）可知固有频率取决于结构自身的质量和刚度，材料的质量越轻，刚度越大，固有频率就越高。但是，又由公式（7）可知P是固定的，要想通过提高风力机塔架的刚度就得增大质量，质量一增加M就增大，固有频率f还是没有很大的变化。而减小塔架的质量，K增大M减小，f会增大的很快。因此，减小塔架自身的质量是提高固有频率最有效的方法。表格3 优化数据Tab.2 The optimization data类别单位序号BBRR一阶频率塔架质量最大位移应力mmmmHzKgmMPa10.030.022.251.51.03631.6581050.69965155.1820.020.022.251.51.0331.6521050.70006142.3730.040.022.251.51.02391.6641050.69806157.0740.030.012.251.51.03031.6531050.70043154.3150.030.032.251.51.02871.6631050.70026136.4860.030.021.751.50.770031.4241051.2277240.0570.030.022.751.51.26071.8921050.441398.6180.030.022.2511.0961.4511050.95235163.1990.030.022.2520.991991.8681050.55928143.09100.0229580.0129581.89791.14790.875521.3411051.2678202.64110.0370420.0129581.89791.14790.888241.3491051.2637202.47120.0229580.0270421.89791.14790.879661.3461051.2707197.69130.0370420.0270421.89791.14790.879911.3531051.2507209.96140.0229580.0129582.60211.14791.26451.6701050.61606126.38150.0370420.0129582.60211.14791.26271.6801050.61109124.81160.0229580.0270422.60211.14791.25891.6741050.612123.5170.0370420.0270422.60211.14791.25571.6841050.61238124.5180.0229580.0129581.89791.85210.837081.6321050.86374205.74190.0370420.0129581.89791.85210.83151.6391050.87598206.16200.0229580.0270421.89791.85210.821591.6431050.86307206.89210.0370420.0270421.89791.85210.809031.6511050.87177206.35220.0229580.0129582.60211.85211.16681.9601050.43215104.99230.0370420.0129582.60211.85211.16991.9701050.43248107.36240.0229580.0270422.60211.85211.1511.9721050.43175105.04250.037420.0270422.60211.85211.16451.9811050.43588112.35 本次结构优化的目的是在满足结构安全性的前提下，尽可能使塔架的固有频率得到提高并且用的材料还最少，能降低共振引起的共振和降低经济成本。由于风力机的功率定了，所以塔架的高度不变，可以作为设计参数的是塔的外底半径、外顶半径、塔底壁厚和塔顶壁厚,分别用D、D、R和R表示。优化用的是目标搜索的方法，设计的搜索目标是塔架的第一阶频率、塔架的质 图12 各个设计参数对目标函数的敏感度量、塔架的 Fig.12 Each design parameters on the sensitivity of the objective function最大位移和最大应力。由表3可以看出优化搜索的范围和基本的取值点，优化的方法是把表内数据进行数据拟合并搜索出最佳的方案。由图12可以看出各个参数对目标的敏感度，图分为4栏，第一到第四栏分别是设计参数对第一阶频率、塔架质量、最大应力和最大位移的敏感度，比如第一栏可以看出，对第一阶频率敏感的是塔底壁厚R和塔顶壁厚R并且R 与频率成正比，R与频率成反比图13 最佳优化方案Fig.13 The best optimization scheme5.1结果分析图13给出了三个最佳的优化方案，第一个方案虽然频率最大，但塔架的质量太重不经济；第三个方案塔架的第一阶固有频率太接近1，不满足公式，使塔架在放电机的激励下产生共振，影响塔架的使用寿命；第二个优化方案频率符合要求塔架的质量比第一个方案轻不少，最大应力为130Mpa，小于材料的许用强度，最大位移也是在可以接受的范围，并且B1和B2很接近，都取25.5mm便于制造，所以方案二是最好并且可行的方案，所以塔架的结构尺寸取为：塔高H=70m；底外径D1=5.084m；顶外径D2=2.023m；底端壁厚B1=25.5mm；顶端壁厚B2=25.5mm。6 结论应用有限元软件ANSYS对大型垂直轴风力机塔柱进行静态学分析、动力