塔鬻风载荷的近似计算
郭庆军韩花丽
(中船重工(重庆)海装风电设备有限公司
摘要
重庆400021)
阐述稳定风和阵风时塔架栽荷的影响.在对塔架风载荷计算方面则结合公司的资源,对比GHBladed的仿真结果,
并参考中国建筑规范.总结出一套行之有效的塔架风载荷计算方法。
关键词
塔架风力发电机组风载荷
文献标识码:A
文章编号:1672—9064(2009)04—0005-03
变化.因此需要对塔架不同高度的圆柱截面直径进行定义。海装风电设备有限公司某型号的塔架部分参数在Bladed软件中定义.如表l。
中图分类号:TK83
在已经给定了塔架高度的情况下.塔架顶部和塔架底部载荷在风力发电机机组塔架设计中至关重要。塔架顶部的载荷主要来自于叶轮的空气动力、离心力、机舱和叶轮的重力等因素:塔架底部的载荷主要来自于塔架顶部载荷以及塔架的风载荷.在已知塔架顶部载荷的情况下通过载荷坐标系转化可以求得塔架顶部载荷对塔架底部载荷的影响.那么如果能够求得塔架的风载荷.将2部分的载荷叠加即可求得塔基的载荷.
目前风电行业对风机运行风况的的定义主要依据是欧洲标准IEC61400—1.此标准中将风的种类定义为稳定风、阵风、湍流风3大类,本文主要阐述了稳定风和阵风对塔架载荷的影响。Bladed软件是风电行业公认的风力发电机整机设计及载荷分析的权威软件.其中包括塔架载荷计算的功能.本文以海装风电设备有限公司某种型号的塔架风载荷计算为例.用本文所介绍的塔架风载荷计算方法所计算的结果同Bladed软件计算的结果相比较以验证塔架风载荷计算方法的正确性.
2稳态风、阵风模型
在同一地点.垂直于地面方向的不同高度处风速的大小不同.多数情况下风速随着高度的增大而增大。在载荷计算中风速随高度的分布遵循常规风廓线模型(NWP),公式如下:
V(:)=y,(娩,)。
考高度:y广参考高度处的风速;a一风切变系数。
(1)
式中:V(z)一高度Z处的风速;z-离地面高度;Z一参
不同地面粗糙度的风场风切变系数不同.如果风场测试参数中没有给出风切变系数时.在塔架载荷计算中风切变系数取标准值.通常国内风切变系数的标准值为l/7。
1塔架的结构
塔架的结构有张线支撑式和悬臂梁2种基本形式。塔架所用的材料可以是木杆、铁管或其他圆柱结构.也可以是钢材做成的珩架结构.大型风力机的塔架基本上是锥形圆柱钢塔架。锥形圆柱钢塔架的特点是塔架的直径随高度的变化而
表l
稳定风风况下塔架风载荷计算结果
图1
不同风切变系数的风廓线图
图1为Zr=70m.Vr=
8.5“s.风切变系数分别为
0.1、1,7、0.2的稳定风的风廓线图.风切变系数越大地表粗糙度越低.在较低的高度就可以达到较高的风速。
图2为初始值为25rrds。幅值为lOm/s.阵风持续时间为10.4s的极端工作阵风。阵风的定义可参考欧洲标准
IEC
61400—l一200l。
作者简介:郭庆军(1982一),男,黑龙江大庆人,硕士,-L:丰a-师。从事风力发电方面的技术研究工作。
万方数据
5
风压/(kN/m2)。
在稳态风和极端工作阵风情况下。风振系数废取1,体形系数“取0.6。
{
=
胁・舻rb2[/30(鼢】协
毂高度。
则塔架整体受力为:
(3)
式中:口广轮毂高度处的风速;p一空气密度,没有给出具体
的空气密度时取标准值1.225kg/m3;^一塔架总高度;k一轮
型
形=r扣[1)o(等)。]礼八^)砒
0
(4)
2‘681012lt1618∞
时间^
式中:“^)一高度h处塔架的迎风面尺度即此处的塔架截面直径:
将(4)离散化
图2极端工作阵风
r
形=;扣h(等)4]红八¨‰
式中:^.一不同直径塔架段的长度。
在BIaded中塔架坐标系定义为:
(5)
2
一一
燃2二L乓::踏
~m§]f一加
圈3
】f
/
●'。
XT指向南;zT垂直向上;YT指向东。(风的偏航角度为Oo时.风向为正北)
在BIaded软件中设定塔架参数.定义风速为8.5m/s.入流角度为00.偏航角度为0。的稳定风和初始速度为25m/s.阵风幅度为15m/s.人流角度为00,偏航角度为0。的阵风。在这2种风况下Bladed软件的塔架风载荷计算结果见表1,表2。
么
-_
由于模拟的稳定风和阵风的偏航角度为00.塔架的Y坐标方向没有风载荷.风载荷全部作用在X坐标方向。这个塔架迎风面的风载荷为1m高度处的X坐标方向的风载荷
Bladed中塔架坐标系
3塔架风载荷计算
垂直于建筑物表面上的风载荷标准值应按下述公式计算:
W‘:忍肌他tlJD
(2)
与68m高度处的X坐标方向的风载荷差值.对时间为ls处的稳定风和时间为8s、15s、22s处的阵风进行塔架风载荷计算.Bladed软件计算结果如下:
1s处稳定风:4.446kN(风速8.5m/s)8s处阵风:38.455kN(瞬时风速25m/s)
11s处阵风:36.402kN(瞬时风速24.4m/s)
13s处阵风:33.202kN(瞬时风速23.5)
15s处阵风:70.267kN(瞬时风速32.4m/s)
部分Matlab程序如下:
loadd:\ys.txt;thita=0
X
式中洲k一风载荷标准值/(kN/m2);区一高度Z处的风振
系数;地一风压高度变化系数;风一风载体型系数;埘『.基本
表2
阵风风况下塔架风载荷计算结果
3.1415926/180;apha=O.2:a=O;for
i=l:15;v=8.5x((ys(i+
l,1)+ys(i,1))/2/70)'aphax
COS(thim);a=a+O.5x0.6×
1.225x(v"2x(ys(i,2)+ys(i+l,2))x(ys(i+1,1)-ys(i,1)))/(2x1000);end
Matlab程序计算结果如
(下转第9页)
万方数据
6
是否能形成类似自然生态系统的生态链或食物链,只有这样才能实现物质与能量的闭路循环和废物最少化。园区成员间是否具备市场规范的供需关系以及供需规模、供需的稳定性等都是在成员的选择中要重点考虑的问题。此外,政府也应该加快制订适宜的保障措施和激励机制.建立和完善市场机制下的管理和服务运作模式.引导有利于循环经济的消费和市场行为.积极整合各种资源优势.探索一条适合中国国情的生态工业园区建设之路.参考文献
1何劲.论可持续发展与我国工业生态化建设.湖南商学院学报.1998
6
(5)
2李运帷.生态工业系统构建及模式研究.中南大学硕士研究生毕业
论文.20073
刘茂福.工业生态化战略模型研究.福建经济管理干部学院学报.
2003(2)
4王东.工业生态化与可持续发展.经济研究,2007(1)
5张艳.我国发展生态工业的策略与措施探讨.工业安全与环保.2005
(4)
肖泉.清洁生产与循环经济知识问答.化学工业出版社.2007
7刘少丽.生态工业园——我国工业发展的新选择.江苏商论.2007
(6)
(上接第4页)
有钻孔揭露情况.深部煤层厚度有变厚的规律和趋势.该区的主要可采的有童子岩组l段33、39、4l号煤层.往东部可能有童子岩组3段的19号等煤层。因此.在“盲区”还潜在一定的找煤空间.
开发利用此类的资源.可采用硐探与钻探相结合的办法。①利用现有最低开采水平的生产系统.往目标区域施工井巷工程.探明构造情况和煤层赋存规律;⑨在矿井边缘块段采用钻探来控制深部资源情况和水文情况.为今后的开发利用提供地质依据.预测上述2个块段有可能增加数百万t的资源量,将为老矿井注入新鲜血液.延长矿井寿命。1.4利用水文治理手段.解放“充水区”的压煤
仙亭煤矿+500m后洋运输大巷于于1991年4月掘进施工时遇见突水.最高水量达423m3/11.稳定水量在230m3/h左右.从而影响了矿井的正常生产.当时在巷道施工挡水墙将突水点封在里面.从此+500m后洋运输大巷停止施工.造成了“充水区”压煤。1999年委托地质队进行水文调查.并经专家审查形成意见:以ZK56与ZK84连线为界.以东划为突水危险区。根据2006年末储量动态检测报告提供数据显示.有2040万t的煤炭资源量暂不能利用.将大大缩短矿井服务年限。
为充分利用煤炭资源.解放“充水区”压煤应作为一个重要的课题进行攻关。按矿井防治水的相关规定。委托由资
质单位编制治理方案,组织专家重新论证突水、导水因素,对症下药,堵住导水通道。除了对未封孔老ZK3可能导水外,对突水点附近的ZK7、ZK3封孑L质量也要进行论证.以及F8断层或派生断层等有可能引起突水的因素进行全面分析论证,找出突水的结点,采用先进的堵水办法.堵住与下部灰岩水体联系的导水通道.以解放被滞压的煤炭资源.还原矿井服务年限.
2结束语
(1)老区复采是老矿井找煤最简单的途径,虽然开采难度大、成本高.但老矿井对复采、薄硬煤层开采有较大的依存度,其采出量占总产量的比例逐年提高。缓解了生产接替,实现了资源同收最大化。
(2)通过对地质构造特征的分析,对原认为暂不能利用的深部资源,投入一定的探巷工程.有望圈出数百万t的资源储量.延长服务年限。将是一条投资省、见效快的找煤途径。
(3)对老矿井边缘“盲区”通过补充勘探(硐探、钻探相结合)手段.可能增加数百万吨的煤炭资源.能使一个资源即将枯竭的老矿井重获新生.
(4)通过矿井水文治理。重新论证突水、导水因素,找出突水的结点.堵住导水通道.希望能将2040万t煤炭资源量解放出来,还原矿井服务年限.具有重大的意义。
(上接第6页)
下:
备有限公司某型号塔架进行计算.其结果同Bladed软件的
ls处稳定风:4.43kN(风速8.5m/s);8s处阵风:38.31kN
塔架风载荷仿真计算结果相比较虽然存在着一定的误差.但在工程中是可以接受的.因此本文所提出的塔架计算方法用作风力发电机组塔架风载荷的近似计算是切实可行的。
(瞬时风速25m/s);lls处阵风:36.5kN(瞬时风速24.4m/s);13s处阵风:33.85kN(瞬时风速23.5);15s处阵风:64.34kN(瞬时风速32.4m/s)。
从计算结果可以看出.利用风载荷计算公式得出的塔架风载荷同Bladed软件的输出结果基本相同.在15s时刻处阵风对塔架载荷的Matlab计算结果同BIaded软件的计算结果有5.927kN的误差,这个误差在工程中是可以接受的。
参考文献
l
中华人民共和回国家标准.风力发电机组安全要求(GB18451.1—
2001).2001
2
中华人民共和园国家标准.建筑结构栽荷规范(GB50009—2001).
2002
4结束语
本文所提出的塔架风载荷计算方法通过对海装风电设
3王承熙.张源.风力发电.中国电力出版社,2003
万方数据
9
塔鬻风载荷的近似计算
郭庆军韩花丽
(中船重工(重庆)海装风电设备有限公司
摘要
重庆400021)
阐述稳定风和阵风时塔架栽荷的影响.在对塔架风载荷计算方面则结合公司的资源,对比GHBladed的仿真结果,
并参考中国建筑规范.总结出一套行之有效的塔架风载荷计算方法。
关键词
塔架风力发电机组风载荷
文献标识码:A
文章编号:1672—9064(2009)04—0005-03
变化.因此需要对塔架不同高度的圆柱截面直径进行定义。海装风电设备有限公司某型号的塔架部分参数在Bladed软件中定义.如表l。
中图分类号:TK83
在已经给定了塔架高度的情况下.塔架顶部和塔架底部载荷在风力发电机机组塔架设计中至关重要。塔架顶部的载荷主要来自于叶轮的空气动力、离心力、机舱和叶轮的重力等因素:塔架底部的载荷主要来自于塔架顶部载荷以及塔架的风载荷.在已知塔架顶部载荷的情况下通过载荷坐标系转化可以求得塔架顶部载荷对塔架底部载荷的影响.那么如果能够求得塔架的风载荷.将2部分的载荷叠加即可求得塔基的载荷.
目前风电行业对风机运行风况的的定义主要依据是欧洲标准IEC61400—1.此标准中将风的种类定义为稳定风、阵风、湍流风3大类,本文主要阐述了稳定风和阵风对塔架载荷的影响。Bladed软件是风电行业公认的风力发电机整机设计及载荷分析的权威软件.其中包括塔架载荷计算的功能.本文以海装风电设备有限公司某种型号的塔架风载荷计算为例.用本文所介绍的塔架风载荷计算方法所计算的结果同Bladed软件计算的结果相比较以验证塔架风载荷计算方法的正确性.
2稳态风、阵风模型
在同一地点.垂直于地面方向的不同高度处风速的大小不同.多数情况下风速随着高度的增大而增大。在载荷计算中风速随高度的分布遵循常规风廓线模型(NWP),公式如下:
V(:)=y,(娩,)。
考高度:y广参考高度处的风速;a一风切变系数。
(1)
式中:V(z)一高度Z处的风速;z-离地面高度;Z一参
不同地面粗糙度的风场风切变系数不同.如果风场测试参数中没有给出风切变系数时.在塔架载荷计算中风切变系数取标准值.通常国内风切变系数的标准值为l/7。
1塔架的结构
塔架的结构有张线支撑式和悬臂梁2种基本形式。塔架所用的材料可以是木杆、铁管或其他圆柱结构.也可以是钢材做成的珩架结构.大型风力机的塔架基本上是锥形圆柱钢塔架。锥形圆柱钢塔架的特点是塔架的直径随高度的变化而
表l
稳定风风况下塔架风载荷计算结果
图1
不同风切变系数的风廓线图
图1为Zr=70m.Vr=
8.5“s.风切变系数分别为
0.1、1,7、0.2的稳定风的风廓线图.风切变系数越大地表粗糙度越低.在较低的高度就可以达到较高的风速。
图2为初始值为25rrds。幅值为lOm/s.阵风持续时间为10.4s的极端工作阵风。阵风的定义可参考欧洲标准
IEC
61400—l一200l。
作者简介:郭庆军(1982一),男,黑龙江大庆人,硕士,-L:丰a-师。从事风力发电方面的技术研究工作。
万方数据
5
风压/(kN/m2)。
在稳态风和极端工作阵风情况下。风振系数废取1,体形系数“取0.6。
{
=
胁・舻rb2[/30(鼢】协
毂高度。
则塔架整体受力为:
(3)
式中:口广轮毂高度处的风速;p一空气密度,没有给出具体
的空气密度时取标准值1.225kg/m3;^一塔架总高度;k一轮
型
形=r扣[1)o(等)。]礼八^)砒
0
(4)
2‘681012lt1618∞
时间^
式中:“^)一高度h处塔架的迎风面尺度即此处的塔架截面直径:
将(4)离散化
图2极端工作阵风
r
形=;扣h(等)4]红八¨‰
式中:^.一不同直径塔架段的长度。
在BIaded中塔架坐标系定义为:
(5)
2
一一
燃2二L乓::踏
~m§]f一加
圈3
】f
/
●'。
XT指向南;zT垂直向上;YT指向东。(风的偏航角度为Oo时.风向为正北)
在BIaded软件中设定塔架参数.定义风速为8.5m/s.入流角度为00.偏航角度为0。的稳定风和初始速度为25m/s.阵风幅度为15m/s.人流角度为00,偏航角度为0。的阵风。在这2种风况下Bladed软件的塔架风载荷计算结果见表1,表2。
么
-_
由于模拟的稳定风和阵风的偏航角度为00.塔架的Y坐标方向没有风载荷.风载荷全部作用在X坐标方向。这个塔架迎风面的风载荷为1m高度处的X坐标方向的风载荷
Bladed中塔架坐标系
3塔架风载荷计算
垂直于建筑物表面上的风载荷标准值应按下述公式计算:
W‘:忍肌他tlJD
(2)
与68m高度处的X坐标方向的风载荷差值.对时间为ls处的稳定风和时间为8s、15s、22s处的阵风进行塔架风载荷计算.Bladed软件计算结果如下:
1s处稳定风:4.446kN(风速8.5m/s)8s处阵风:38.455kN(瞬时风速25m/s)
11s处阵风:36.402kN(瞬时风速24.4m/s)
13s处阵风:33.202kN(瞬时风速23.5)
15s处阵风:70.267kN(瞬时风速32.4m/s)
部分Matlab程序如下:
loadd:\ys.txt;thita=0
X
式中洲k一风载荷标准值/(kN/m2);区一高度Z处的风振
系数;地一风压高度变化系数;风一风载体型系数;埘『.基本
表2
阵风风况下塔架风载荷计算结果
3.1415926/180;apha=O.2:a=O;for
i=l:15;v=8.5x((ys(i+
l,1)+ys(i,1))/2/70)'aphax
COS(thim);a=a+O.5x0.6×
1.225x(v"2x(ys(i,2)+ys(i+l,2))x(ys(i+1,1)-ys(i,1)))/(2x1000);end
Matlab程序计算结果如
(下转第9页)
万方数据
6
是否能形成类似自然生态系统的生态链或食物链,只有这样才能实现物质与能量的闭路循环和废物最少化。园区成员间是否具备市场规范的供需关系以及供需规模、供需的稳定性等都是在成员的选择中要重点考虑的问题。此外,政府也应该加快制订适宜的保障措施和激励机制.建立和完善市场机制下的管理和服务运作模式.引导有利于循环经济的消费和市场行为.积极整合各种资源优势.探索一条适合中国国情的生态工业园区建设之路.参考文献
1何劲.论可持续发展与我国工业生态化建设.湖南商学院学报.1998
6
(5)
2李运帷.生态工业系统构建及模式研究.中南大学硕士研究生毕业
论文.20073
刘茂福.工业生态化战略模型研究.福建经济管理干部学院学报.
2003(2)
4王东.工业生态化与可持续发展.经济研究,2007(1)
5张艳.我国发展生态工业的策略与措施探讨.工业安全与环保.2005
(4)
肖泉.清洁生产与循环经济知识问答.化学工业出版社.2007
7刘少丽.生态工业园——我国工业发展的新选择.江苏商论.2007
(6)
(上接第4页)
有钻孔揭露情况.深部煤层厚度有变厚的规律和趋势.该区的主要可采的有童子岩组l段33、39、4l号煤层.往东部可能有童子岩组3段的19号等煤层。因此.在“盲区”还潜在一定的找煤空间.
开发利用此类的资源.可采用硐探与钻探相结合的办法。①利用现有最低开采水平的生产系统.往目标区域施工井巷工程.探明构造情况和煤层赋存规律;⑨在矿井边缘块段采用钻探来控制深部资源情况和水文情况.为今后的开发利用提供地质依据.预测上述2个块段有可能增加数百万t的资源量,将为老矿井注入新鲜血液.延长矿井寿命。1.4利用水文治理手段.解放“充水区”的压煤
仙亭煤矿+500m后洋运输大巷于于1991年4月掘进施工时遇见突水.最高水量达423m3/11.稳定水量在230m3/h左右.从而影响了矿井的正常生产.当时在巷道施工挡水墙将突水点封在里面.从此+500m后洋运输大巷停止施工.造成了“充水区”压煤。1999年委托地质队进行水文调查.并经专家审查形成意见:以ZK56与ZK84连线为界.以东划为突水危险区。根据2006年末储量动态检测报告提供数据显示.有2040万t的煤炭资源量暂不能利用.将大大缩短矿井服务年限。
为充分利用煤炭资源.解放“充水区”压煤应作为一个重要的课题进行攻关。按矿井防治水的相关规定。委托由资
质单位编制治理方案,组织专家重新论证突水、导水因素,对症下药,堵住导水通道。除了对未封孔老ZK3可能导水外,对突水点附近的ZK7、ZK3封孑L质量也要进行论证.以及F8断层或派生断层等有可能引起突水的因素进行全面分析论证,找出突水的结点,采用先进的堵水办法.堵住与下部灰岩水体联系的导水通道.以解放被滞压的煤炭资源.还原矿井服务年限.
2结束语
(1)老区复采是老矿井找煤最简单的途径,虽然开采难度大、成本高.但老矿井对复采、薄硬煤层开采有较大的依存度,其采出量占总产量的比例逐年提高。缓解了生产接替,实现了资源同收最大化。
(2)通过对地质构造特征的分析,对原认为暂不能利用的深部资源,投入一定的探巷工程.有望圈出数百万t的资源储量.延长服务年限。将是一条投资省、见效快的找煤途径。
(3)对老矿井边缘“盲区”通过补充勘探(硐探、钻探相结合)手段.可能增加数百万吨的煤炭资源.能使一个资源即将枯竭的老矿井重获新生.
(4)通过矿井水文治理。重新论证突水、导水因素,找出突水的结点.堵住导水通道.希望能将2040万t煤炭资源量解放出来,还原矿井服务年限.具有重大的意义。
(上接第6页)
下:
备有限公司某型号塔架进行计算.其结果同Bladed软件的
ls处稳定风:4.43kN(风速8.5m/s);8s处阵风:38.31kN
塔架风载荷仿真计算结果相比较虽然存在着一定的误差.但在工程中是可以接受的.因此本文所提出的塔架计算方法用作风力发电机组塔架风载荷的近似计算是切实可行的。
(瞬时风速25m/s);lls处阵风:36.5kN(瞬时风速24.4m/s);13s处阵风:33.85kN(瞬时风速23.5);15s处阵风:64.34kN(瞬时风速32.4m/s)。
从计算结果可以看出.利用风载荷计算公式得出的塔架风载荷同Bladed软件的输出结果基本相同.在15s时刻处阵风对塔架载荷的Matlab计算结果同BIaded软件的计算结果有5.927kN的误差,这个误差在工程中是可以接受的。
参考文献
l
中华人民共和回国家标准.风力发电机组安全要求(GB18451.1—
2001).2001
2
中华人民共和园国家标准.建筑结构栽荷规范(GB50009—2001).
2002
4结束语
本文所提出的塔架风载荷计算方法通过对海装风电设
3王承熙.张源.风力发电.中国电力出版社,2003
万方数据
9