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

基于ANSYS的垂直轴风力机塔架的力学分析及结构优化（山东科技大学机电学院，山东青岛266510）摘要:基于ANSYS对大型垂直轴风力机塔架进行静态学分析、动力学分析及结构优化。首先根据风力机的功率确定塔架的大体尺寸，然后对塔架进行静力分析，并对塔架承受最大载荷时的最大应变进行强度校核；动力学分析包括有预应力的模态分析、地震谱响应分析。ANSYS优化通过改变塔架的厚度使其固有频率增大，从而避免因塔架在发电机的激励下共振产生的破坏。通过对塔架的力学分析确保了塔架性能的可靠性，为塔架的安全性提供了理论依据，并且在满足力学性能的前提下使塔架的质量最轻，降低了塔架的制造成本。关键词:垂直轴风力机塔架；静力学分析；动力学分析；地震谱响应分析；结构优化；ANSYS中图分类号:TK83文献标志码:A文章编号:1672-3767(2011)07-0000-07MechanicalPerformanceAnalysisAndStructureOptimizationOfTowerOfVerticalAxisWindTurbineBasedOnANSYS（CollegeofMechanicalandElectronicEngineering，ShanDongUniversityofScienceandTechnology，QingDao266510,China）Abstract:Staticsanalysis,dynamicsanalysisandstructuraloptimizationoflargetowerofverticalaxiswindturbinewerepresentedbasedonANSYS.Firstofall,accordingtothepowerofthewindgeneratorsthesizeofthetowerwasdetermined.Andthenstaticsanalysisforthetowerwascarriedout.Andthestrengthcheckforthemaximumstrainofthetowerwasdiscussedunderthemaximumload.Dynamicsanalysisincludingprestressedmodalanalysis,seismicresponsespectrumanalysiswasperformed.UsingANSYSoptimizationbychangingthethicknessofthetoweroftheinherentfrequencyincreases,soastoavoidthegeneratorfortowerundertheexcitationofresonanceproduceddamage.Throughtothemechanicalanalysisofthetowertoensurethereliabilityofthetoweroftheperformanceforthesafetyofthetower,toprovidethetheoreticalbasisforthemechanicalproperties,andmakethequalityofthetowerofthelightestinmeetthepremiseofmechanicalproperties,reducethetowerofmanufacturingcost.Keywords:towerofverticalaxiswindturbine;staticsanalysis;dynamicsanalysis;seismseismicresponsespectrumanalysis;structuraloptimization;ANSYS风能资源主要集中在50米以上的高空，是低空资源的2倍左右，大容量的风力机机组便于大规模风能的开发利用，因此大功率的风机具有很好的发展前景。在风力发电机中，塔架占成本的25%左右[1]。塔架是风力机的承载部分，塔架的设计是风力发电机组不可缺少的一部分，尤其对于大功率的风力机组[2]。结构振动导致的材料疲劳通常是系统损坏的主要原因之一，振动还会产生噪音，影响周围环境，所以塔架的力学分析对于风力机的安全可靠性能具有重要的参考价值。目前，对风力机的研究大都集中在叶片、发动机性能等的研究上，对于塔架的研究很少，尤其是针对大塔架的研究更为少见[3]，垂直轴风力塔架的分析在国内基本是空白。本研究是在前人分析水平轴风力机的基础上，根据垂直轴风力机的实际受力对塔架进行加力并分析，为大型垂直轴的风力机塔架的制造提供了理论依据。本文首先对垂直轴风力机塔架的模型尺寸的计算与建模，接着对模型加载载荷进行静态的受力分析，在受力的情况下对模型进行有预应力的模态分析，得到塔架的前6阶模态及振型。在模态的基础上，加载模拟的地震频率波动进行谱分析，最后针对塔架模型进行结构的优化，优化的参数是塔架底的壁厚和顶壁厚以及塔底和塔顶的外半径，目标函数是塔架的第一阶频率、质量、最大应力和最大位移。1垂直轴风力机塔架的计算模型大型风力机组大都采用圆锥形塔柱，本研究分析的是1MW的垂直轴风力机的塔柱，根据《高耸机构设计规范》GB50135-2006，选择塔柱的尺寸如图1所示：塔高H=70m，底外径D1=4.5m，顶外径D2=3m，底端壁厚B1=30mm，顶端壁厚B2=20mm，底座长L=0.5m，材料为Q345钢，材料弹性模量E=2.06×1011，泊松比µ=0.3，密度ρ=7.850×103kg/m3,许用强度σ=237MPa，许用抗拉强度σ1=425MPa。此塔架为变截面的锥形圆锥体，由于底端是固定的，可以看成悬臂梁结构，塔架的主要载荷有风的水平推力F，机舱、叶片和轮毂的重力G1，塔架的自重G2，轮毂的扭矩M等，作用在塔身上的载荷共5种，受力简图如图2所示。1.1塔架产生的主要载荷1）风轮的水平推力：AVCF221(1)式中ρ——空气密度（1.293kg/m3）；图1塔架结构模型V——风速（定为30m/s）；Fig.1ThestructuremodeloftowerA——风轮的扫风面积（取1700m2）；C——推力系数（一般处切入风速时取1左右，切出风速取0.5左右，根据贝茨公式计算，C=8/92）机舱、叶片和轮毂的力：G1=mg(2)式中m为机舱、叶片和轮毂的总重量（质量取39000kg）；g为重力加速度（取9.8m/s2）3）塔架自身的重量：G24）轮毂扭矩：M=P/n(3)式中P为发电机的输出功率n为风轮的转速（r/s）图2塔架的受力简图Fig.2Theforcediagramoftower5）作用在塔架上的风压：20202)(2121VHZVPaZZ]4[(4)式中ρ为空气的密度；Z为塔架的高度；V0是高度为Z处的风速;Z为高度；H0为高度；V0是高度为10m处的风速；a是地面粗糙度和地面风的切变系数。有公式（1）可得：F=879240(N)；有公式（2）可得：G1=382200(N)；把P=1×106瓦，n=1/3(r/s)带入公式（4）可得：M=3×106(NM)2用ANSYS对塔架进行静力分析2.1建模、加载与分析建模：首先根据上面塔架模型数据在DM中画出塔架的二维模型然后以中心轴为旋转轴进行旋转就得到了塔架的模型，建模完成。加载与分析：进入静力分析模块对模型进行网格划分然后施加边界条件和载荷，步骤如下：1）固定塔架的底端；2）在塔架顶部加载力F和机舱、叶片和轮毂的重力G1以及轮毂扭矩M；3）在塔架上加载塔架的重力G2和风压PZ，结果设置为查看总体的应力和总体的应变，最后进行求解。由求解的结果图3可以得到塔架的最大位移为0.69965m，由图4得到塔架的最大应力图3塔架的位移/m为155.18MPaFig.3Thedisplacement0ftower/m2.2结果分析与讨论根据材料的刚度要求塔架的最大变形在高度H=70m的1%一下即：7.0100maxHD(m)]5[（5）计算得出的0.69965(m)0.7(m),不超出最大变形的要求，计算得到得最大应力为155.18MPaσ=237MPa,强度符合要求，所以塔架结构和静力学性能满足正常工作的要求。3塔架的模态分析图4塔架的应力/Pa塔架的模态分析对塔架的稳定性非常重要，Fig.4Thestressoftower/Pa模态分析主要是计算塔架的固有频率和阵型，它也是动力学分析的基础，同时也是地震谱响应分析的前序和结构参数优化的参数标准，如果塔架受力，就应该进行有预应力的模态分析，分析的结果才有意义。塔架的前几阶模态对稳定性的影响比较大，因此分析塔架的前六阶模态。3.1塔架的有限元动力学方程动力学分析同静力学分析都是把塔架离散成有限个单元体，动力学分析要考虑单元特性的同时还考虑单元的惯性力和阻尼力等因素。动力学方程：fKXXCXM,,,（6）式中：M为塔架的质量矩阵；C为塔架的阻尼矩阵；K为总体刚度矩阵；X为的节点位移；f为塔架受的外力]6[3.2ANSYS对塔架模态分析的结果从图5到图10可是得到塔架的一到六阶振型、各阶振型下产生的最大位移和固有频率。如在图5中能得到塔架的振型是一阶的，最大的位移在塔架的顶端为5.6769×10-3m，一阶固有频率为1.0363Hz。图5塔架的第一阶模态振型Fig.5Thefirstordermodalshapeoftower图7塔架的第三阶模态振型Fig.7thethirdordermodalshapeoftower图9塔架的第五阶模态振型Fig.9thefifthordermodalshapeoftower图6塔架的第二阶模态振型Fig.6Thesecondordermodalshapeoftower图8塔架的第四阶模态振型Fig.8thefourthordermodalshapeoftower图10塔架的第六阶模态振型Fig.10thesixthordermodalshapeoftower表1塔架的前六阶固有频率Tab.1ThethefirstsixordernaturalfrequencyoftowerMode123456Freq(Hz)1.03631.03974.95384.969912.42712.4633.3结果分析与讨论根据公式Nf2可知，当风力发电机以20r/min的转速运行时，风轮旋转产生的频率为0.3333Hz，风叶的个数是3，f=1Hz，由表1可知，这个频率十分接近塔架的第一、二阶固有频率，频率要满足公式(f1―f)/f10%，才能不引起塔架的共振，而f1=1.0363，(f1―f)/f=0.03630.1,所以塔架在发电机的激励下是产生共振的，共振对结果的破坏非常大，因此要对塔架的结构进行优化，使它的前两阶频率得到提高，降低共振引起的破坏。4塔架的地震谱响应分析表2地震模拟数据Tab.2Theearthquakesimulationdata序号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可得到地震谱响应