垂直轴风力机塔架稳定性研究.docx

垂直轴风力机塔架稳定性研究梅毅1，曲建俊1，许明伟1(1．哈尔滨工业大学机电工程学院，哈尔滨150001) 擅蔓i以一种ikW垂直轴风力机钢制管状塔架为研究对象，采用欧拉法、工程算法和基于有限元的非线性屈曲分析方法计算了塔架屈曲的临界载荷。研究结果表明，是否考虑几何非线性对l临界载荷的计算由重要影响，欧拉法和工程算法计算得到的临界载荷值偏大。在垂直轴风力机管状塔架设计时，为保证结构的稳定性，适合采用有限元方法进行非线性屈曲分析。关■阅：垂直轴风力机：塔架：有限元；屈曲：非线性。引言风能是一种具有较好开发利用前景的可再生资源。水平轴风力机是目前技术比较成熟，市场占有 率较高的风能利用设备，但水平轴风力机并不适应所有的场合。直叶片垂直轴风力机由于具有无需对 风装置，噪声小以及加工制造相对简单等优势，更适合作为小型离网型机组安置于边远山村、边防哨 所等地域，因而引起学者们对垂直轴风力机相关技术问题的广泛研究【1—3】。钢制管状塔架是小型离网型垂直轴风力机发电机组采用的主要支撑形式。管状塔架是一种典型的 薄壳结构，其特点是跨度大，厚度小，自重相对较轻，这种薄壳结构在工程中应用十分广泛，其极限 承载力主要取决于结构的稳定承载能力，因此在工程设计中对此类结构的整体屈曲分析是不可或缺的。风力机塔架失稳是指当处于受载状态的塔架在外界干扰微小时，塔架状态的扰动发生较大变化(即失去原有的平衡形式或几何形状)的现象。研究塔架稳定性的目的在于确定塔架的临界载荷，对保障 直叶片垂直轴风力机塔架及机组的安全具有重要意义。鉴于目前已有文献中有关小型垂直轴风力机塔架稳定性研究较少的情况，本文以一种lkW垂直 轴风力机塔架为研究对象，分别采用欧拉法、风力机塔架稳定性的工程算法和基于有限元的非线性屈曲分析法计算了该塔架的屈曲l临界载荷，并比较和分析了三种方法计算结果的差异，研究结论可为垂直轴风力机的工程设计提供参考。1研究模型普通垂直轴风力机的发电机置于地面，而根据国家标准【4】，离网型风力机塔架的高度一般设计在 5—12m之间，这将导致发电机的主轴过长，给加工和装配带来极大困难。哈尔滨工业大学机电学院风电技术装备研究所研制的lkW离网型垂直轴风力机，采用外转子发电机置于塔架之上，有效缩短了传动链的长度，使机组传动结构更加紧凑，有利于零部件的加工和机组现场安装。风轮结构如图1所 示，叶片为玻璃钢材料，其它结构如支撑臂和主轴等均为钢制结构，总重量约为500kg。机组塔架采 用无缝热轧钢管，如图2所示，材料为O．235，钢管长5m，外径0。168m，壁厚9mm。在下文的屈曲分析中，欧拉法将该管状塔架视作轴心理想压杆，从压杆挠曲的近似线微分方程出 发，导出屈曲临界载荷的计算公式。后两种方法则将塔架作为轴心受压的薄壁圆柱壳来研究其屈曲临界载荷。图11kw离网型垂直轴风力机风轮图2风力机组塔梁所用无缝钢雷2塔架屈曲临界载荷研究2．1欧拉法欧拉法通过研究结构在载荷作用下是否存在两种不同平衡位形的可能性问题，来求解屈曲临界参 数，相对简单。本文研究的塔架可抽象为一端固定，一端自由压杆，如图3所示。为验证垂直轴风力机塔架的稳定性，首先应导出计算压杆临界压力的欧拉公式。Fla}。rⅪl|一月■啊％％图3压杆受力示意图压杆微弯曲时的挠度方程为：日万d2m=M(x)=F(d-y)令昙=∥，(1)式变为：垡+k2m=k2万(2)出2对于二阶常系数微分方程，其通解为m=Asinkx+BCOSkx+万(3)∞’=Akcoskx—Bksini。c(4)代入边界条件：x=0时，to=0，解得：A=0，8=一6，w=6(1-COSkx)将x：I时，y=6代入可得：coskl=0，kl=O．5(2n+1l，r则有露2：委：育(2n+1)2n2E，(5)～r2，)2F：—(2n_+忑1)2f7r2一EI(6)‘(2Zy使杆件保持微弯曲状态的最小压力为临界压力，故取n=0，则压杆的欧拉公式为：吒：_zr2rEl吒2可(7)o7，由于欧拉公式是由弯曲变形的微分方程导出的，而材料服从胡克定律是该微分方程的基础，所以只有当临界应力小于比例极限op时，欧拉公式才是正确的，即压杆的柔度入需大于或等于极限值kl。 入和k1分别由式f8l和式(9}求得。州S，／万扣√百塔架的材料为Q235，E=206X109Pa，op=200×106Pa，1=2．86x10．5m4，A=0．004m2，pt=2，1=5。经计算可得2k=118．3，kl=100．8。kkl，可以使用欧拉公式。根据式(7)可得Fcr：5．8X105N，即塔架屈曲的 临界压力为5．8×105N。2．2上丰茔算法风力机塔架屈曲临界载荷的工程算法由薄壁圆柱壳的线性屈曲理论导出。根据线性屈曲理论有：删国+等窘岬4c》㈣、2NV4(意MV40舻209)_0a(