基于计算流体动力学,采用STAR-CCM+软件建立垂直轴风力机的空气动力学模型,计算风机旋转时叶片所受气动荷载.将气动荷载在MATLAB软件中进行拟合,获得时序性风机叶片荷载函数,并将其施加到基于ANSYS软件建立的垂直轴风力机的结构模型上,从而进行动力学分析,得到风机结构的位移响应时程.为抑制风力机的振动,在风机不同位置处施加基于COMBIN40单元而模拟的阻尼器,计算阻尼器耗能减振的能力.
1 风力机力学模型
本文的技术路线如图1所示.H型垂直轴风力机的几何模型如图2所示(其中数字为节点编号),其通常由2~6片叶片组成.目前2叶片或3叶片风机在风电市场上应用较广[1],兹选取典型3叶片H型垂直轴风力发电机为对象,研究风机旋转时风荷载对风机振动稳定性的影响.垂直轴风机多采用对称翼型,故本文中风机的叶片翼型截面选择NACA-0021[12],D型主梁处材料采用碳纤维,叶片蒙皮处材料采用玻璃纤维.主轴与支杆材料采用Q345钢材[13],弹性模量E=200 GPa,泊松比μ=0.17,密度ρ=7.8×103 kg/m3.为减小计算量,选用某大型10 MW风机的100∶1缩尺模型.风机模型的几何尺寸为:弦长C=0.265 m,支撑杆长D=0.8 m,叶片高度H=1.2 m,主轴和撑杆截面皆取圆截面,主轴直径为0.2 m,撑杆直径为0.1 m.
图1
图1
垂直轴风力机风振与减振技术路线图
Fig.1
Flowchart of wind vibration and vibration reduction of VAWT
图2
图2
H型垂直轴风力发电机的几何模型
Fig.2
Geometric model of a H-rotor type vertical axis wind turbine
2 分析方法
2.1 风力机风荷载模拟
式中:ui 、uj分别为i、j方向的速度分量,i、j为坐标轴方向;t为时间;xi、xj为坐标;ρ为流体密度;p为压力;ν为动力黏度系数;上标“’”代表来流的脉动值.
由于本文的研究目的侧重于风机减振,所以对于风机的气动荷载采用一些简化计算,具体如下:① 风力机处于自然大气边界层中,因其叶片高度不大,故可不考虑风速梯度的影响.② 考虑到在实际工程中,近地面的湍流强度较大.且当风力机转速恒定时,来流风速越低,风机叶尖速越大,因风轮转动引起的湍流强度随之增大.在这种高叶尖速比的条件下,湍流强度的大小会对风力机的气动力造成一定影响[16].为简化模型,本文的计算条件参考相关风洞实验条件设置[16],其中来流风速为8 m/s,湍流度为0.005.③ 因本文的研究对象是风力机的缩尺模型(1∶100),故本文的雷诺数是原模型的1/100,根据文献[16],雷诺数越大,风机的空气动力学性能越好,同时,由于叶片升力系数的增大,转子转矩也更大.因此本文所采用的缩尺模型的气动载荷特性与原模型存在一定的偏差,模型的折算气动载荷略小于实际载荷.④ 在CFD数值模拟过程中,为简化计算,未考虑叶片结构变形的影响.
图3
图4
图4
叶片气动荷载时程曲线与函数拟合结果对比
Fig.4
Comparison of blade aerodynamic load time history curve and function fitting results
由于风机的气动荷载呈现周期性,所以可采用周期函数对其进行拟合.在MATLAB中,采用多阶正弦函数对风机的切向力、法向力及转矩时程曲线进行拟合,获得相应的荷载函数,便于在ANSYS中对结构进行瞬态分析.以叶片的法向力为例,任一时间序列的荷载函数可近似模拟为
式中:fn(t)为t时间下的叶片法向力;ai、bi及ci分别为第i个正弦函数的各项系数.
如图4所示,法向力、切向力及转矩的荷载函数曲线的拟合效果较好,几乎与气动荷载时程曲线保持一致.因此,采用多阶正弦函数对气动荷载数据进行拟合的方法是可行的.
2.2 结构风致振动响应分析
图5
风机结构在荷载激励下的运动方程可以表示为
式中:M、c及K分别为结构的质量矩阵、阻尼矩阵及刚度矩阵;
用有限元方法求解上述方程,可以得到各节点的风致响应.
2.3 阻尼器模拟
图6
图6
加阻尼器后风力机模型示意图与弹簧-阻尼单元
Fig.6
VAWT model with dampers and spring-damping unit
阻尼元件的刚度和阻尼系数对于减振效果的影响非常明显.具体体现为:当阻尼元件刚度不变时,峰值位移在整体范围内随着阻尼系数的增大而减小,当刚度足够大时,阻尼系数的改变对峰值位移的影响不大.当阻尼系数不变时,峰值位移在整体范围内随着刚度的增大而减小,当阻尼系数足够大时,刚度的改变对峰值位移的影响不大.基于上述规律并经过反复试算,最终假定阻尼元件刚度为5×107 N/m,阻尼系数为2.5×107 N/(m·s).
3 数值算例验证
3.1 动力作用下悬臂梁变形验证
图7
图7
悬臂梁受动力荷载的时程响应
Fig.7
Time-history response of cantilever beam at dynamic load
3.2 时间步长选取
对于运用ANSYS分析瞬态动力学问题,选取不同的时间步长和时间步会对结果产生影响.时间步长越小,计算精度越高,但相应会造成计算资源的浪费,因此在数值模拟前应当先选取合适的时间步长.兹对计算模型进行在3种时间步长下的模拟,即:0.004、0.002及0.001s,分别得到对应的结构关于x轴的底部弯矩M0(见表1):
式中: hi 为结构物的高度;P(x)为x高度处结构表面压力;A,Ai为受力面积.
表1 不同时间步长下底部弯矩比较
Tab.1
| 时间步长/s | 底部弯矩/(N·m) |
|---|---|
| 0.004 | 117.56 |
| 0.002 | 107.47 |
| 0.001 | 106.68 |
从表1看出,当时间步长为0.002 s时,Mx为107.47 N·m,与时间步长取0.001 s时所得的结果相差仅0.7%,可看作近似相等.故时间步长为0.002 s时,既能保证所得结果的精确性,也能节省计算资源.兹选取时间步长为0.002 s.
4 数值结果与分析
4.1 风荷载模拟
图8
4.2 风机结构的风致响应分析
在运用ANSYS软件计算垂直轴风力机的振动响应时,由于垂直轴风力机旋转过程中叶片的气动荷载随时间呈周期性变化,所以,将周期性变化的气动荷载加载到风机结构上,即是把原来的风载恒定,风机转动的问题转化成了风机不动,而风载随时间变化,从而计算风机的振动响应.在风机结构响应分析中,为简化计算,未考虑结构整体转动的自由度.同时,本文中未考虑风机转动对风振响应的影响,风机转动可能会使风机的风振响应更为复杂,在后续研究中将进一步展开分析.
分析发现t=10 s时,结构的振动幅值较大,故选取t=10 s时刻下的总位移云图分析,如图9所示, 图中Usum为风力机结构的总位移.图9(a)所示为t=10 s时,原结构(模型1)的总位移云图.可见,模型1的最大位移发生在位于叶片顶端的节点16处,达到2.7 mm.即在该时刻下,风机的危险点为节点16,由于风机叶片的对称性,可重点关注叶片顶端节点.图9(b)所示为t=10 s时,在上部支撑杆处加了阻尼器后的模型2总位移图,此时最大位移发生在位于叶片底端的节点9处,数值为1.5 mm,下降达44.4%.图9(c)所示为t=10 s时,在下部支撑杆处加了阻尼器后模型3的总位移图,此时最大位移仍旧发生在位于叶片顶端,数值为1.6 mm,下降达40.7%.
图9
图9
不同垂直轴风力机模型总位移云图对比
Fig.9
Comparison of total displacement contour map of different VAWT models
为更清晰地对比结构的三向位移,用图10表示t=10 s时3个模型在x、y及z向的位移云图,其中,Ux、Uy及Uz分别对应模型在x、y及z向的位移.由图10对比x向位移数值,可见模型1的Ux分布在-2.7~1.0 mm,模型2的Ux分布在-0.5~0.8 mm,模型3的
图10
图10
风力机结构三向位移云图对比
Fig.10
Comparison of three-way displacement contour maps of VAWT
由以上分析可知,原模型危险点出现在叶片顶端,故取节点8为对象,研究3个模型在风机叶片顶端的风致响应.又因x、y向位移数值较大,z向位移较小,故仅对比节点8的x、y向位移响应.图11所示为3个模型节点8的风致响应对比.图中Ux,8表示风力机节点8在x方向的位移,Uy,8表示风力机节点8在y方向的位移.可见,模型2和模型3在叶片顶端节点的位移抑制效果较为显著.表2所示为3个模型的节点8的三向位移改善率统计数据.可见,模型2在x方向的平均位移改善率约为44%,在y方向的平均位移改善率约为54%,在z方向的平均位移改善率约为5%.模型3在x方向的平均位移改善率约为40%, y方向的平均位移改善率约为53%,但在z方向未见明显改善.
图11
图11
3个模型节点8的风致响应对比
Fig.11
Comparison of wind-induced response of three models at node 8
表2 3个模型节点8的三向风致响应统计
Tab.2
| 模型 | i | Ui,8/mm | 改善率/mm | |||
|---|---|---|---|---|---|---|
| 最大 | 最小 | 最大位移 | 最小位移 | |||
| 1 | x | 1.48 | -0.61 | — | — | |
| y | 0.92 | -2.39 | — | — | ||
| z | 0.29 | -0.29 | — | — | ||
| 2 | x | 0.61 | -0.43 | 58.7% | 29.5% | |
| y | 0.48 | -0.93 | 47.8% | 61.1% | ||
| z | 0.26 | -0.29 | 10.3% | 0.0% | ||
| 3 | x | 0.82 | -0.40 | 44.6% | 34.4% | |
| y | 0.45 | -1.08 | 51.1% | 54.8% | ||
| z | 0.31 | -0.31 | -6.9% | -6.9% | ||
5 结论
本文针对H型垂直轴风力机振动问题,基于计算流体动力学方法,运用STAR-CCM+软件数值模拟风机旋转时叶片所受气动荷载,将气动荷载在MATLAB软件中实施函数拟合,获得时序性风机叶片荷载函数,再基于ANSYS软件开展时程分析,得到风机结构的三向位移响应时程.在风机不同位置处施加基于COMBIN40单元而模拟的阻尼器,数值模拟计算其耗能减振能力.主要结论如下:
(1) H型垂直轴风力机的风致位移响应曲线具有一定的周期性.t=10 s时,风机最大位移出现在叶片的顶端,风致x、y向位移处于同一数量级,z向位移较小.
(2) 在H型垂直轴风力机的主轴与支杆连接处布置阻尼器可降低结构位移响应,总位移最大降幅达44%.
(3) 阻尼器布置的位置与结构的位移改变密切相关.t=10 s时,在近风机叶片顶端连杆处布置阻尼器,结构最大位移出现在风机叶片底端,位移下降达44%;在近风机叶片底端连杆处布置阻尼器,最大位移则依旧出现在了风机叶片顶端,位移下降达40.7%.
本文成果可为垂直轴风力机减振研究提供一定的技术参考.
参考文献
全球十大风电国家的装机容量
[J].
Installed capacity of ten wind power countries in the world
[J].
Investigation of wake characteristics for the offshore floating vertical axis wind turbines in pitch and surge motions of platforms
[J].DOI:10.1016/j.energy.2018.10.101 URL [本文引用: 1]
垂直轴风力机在风力发电中的应用现状及展望
[J].
Current situation and prospect on vertical axis wind turbine in wind power generation
[J].
海上风电机组系统动力学建模及仿真分析研究
[D].
Study on systematic dynamic model and simulation for offshore wind turbine
[D].
海上风力机气动特性及新型浮式系统
[D].
Offshore wind turbine aerodynamic performance and novel floating system
[D].
竖轴H型叶轮及导流罩流体动力性能数值模拟
[D].
Numerical simulation on fluid dynamic performance of H-shaped vertical axis turbine and duct
[D].
基于CFD方法的竖轴叶轮性能模拟及翼型优化
[D].
Capability simulation and airfoil optimization of vertical-axis turbine based on CFD method
[D].
垂直轴风机的结构受力分析及其优化
[J].
Mechanical analysis and optimization of structures in vertical wind turbine
[J].
A 3D FEM model for floating wind turbines support structures
[J].DOI:10.1016/j.egypro.2017.10.344 URL [本文引用: 1]
纵向阻尼器对大跨度斜拉桥桥塔的抗震性能影响分析
[J].
Effect of longitudinal damper on seismic performance of long-span cable-stayed bridge tower
[J].
昆明某高层建筑消能减震设计
[J].
Energy dissipation design for a high-rise building in Kunming
[J].
垂直轴风力机锥齿轮式变桨距结构的优化研究
[D].
The cone gear pitch structure optimal design of vertical axis wind turbine
[D].
垂直轴风力机主轴结构优化设计
[J].
Structure optimization design for cylindrical vertical-axis wind turbines
[J].
基于CFD与LS-DYNA耦合技术的风致冷却塔倒塌破坏分析
[C]//陆新征.第28届全国结构工程学术会议论文集(第I册).
Collapse and failure analysis of wind-induced cooling tower based on CFD and LS-DYNA coupling technology
[C]//LU Xinzheng. Proceedings of the 28 th National Conference on Structural Engineering(I).
基于标准k-ε模型的平衡大气边界层模拟
[J].
Application of standard k-ε model to simulate the equilibrium ABL
[J].
Study on power performance for straight-bladed vertical axis wind turbine by field and wind tunnel test
[J].DOI:10.1016/j.renene.2016.01.002 URL [本文引用: 3]
Aerodynamic performance prediction of straight-bladed vertical axis wind turbine based on CFD
[J].
小型风力机叶片建模及其双向耦合仿真分析
[D].
Modeling of small wind turbine blades and its bidirectional coupling simulation analysis
[D].
Three-dimensional Improved Delayed Detached Eddy Simulation of a two-bladed vertical axis wind turbine
[J].DOI:10.1016/j.enconman.2016.11.067 URL [本文引用: 2]
A method of computation for structural dynamics
[J].
应用液体黏滞阻尼器的高层结构抗风工程设计方法与实例
[J].
Design methods and examples of wind resistant engineering for high-rise building with liquid viscous damper
[J].
基于增量动力分析法的高层建筑-阻尼器系统地震易损性分析
[J].
IDA based seismic fragility analysis of high-rise building-damper system
[J].
多维减振阻尼器力学性能研究
[J].
Research on mechanical properties of multi-dimensional vibration control damper
[J].
机翼模拟梁的动力学仿真分析
[D].
Dynamic simulation analysis of a wing analogy beam
[D].
Aerodynamic noise prediction of a Horizontal Axis Wind Turbine using Improved Delayed Detached Eddy Simulation and acoustic analogy
[J].DOI:10.1016/j.enconman.2015.04.011 URL [本文引用: 1]
