0 引言
国内低风速地区市场竞争日趋激烈,通过提升塔架高度获得更大的发电量是必然趋势。随着塔架高度的增加,传统钢塔面临制造、运输、频率穿越等问题[1]。钢-混组合塔架将混凝土塔和钢塔两种结构合理组合,能够满足风电机组对高度的要求,并且具有更好的安全性和经济性。塔架作为风电机组重要承载部件,承受着机组运行过程中产生的各种荷载,塔架不仅要有足够的强度、刚度和稳定性,其固有频率需避开风轮运行频率,避免与叶轮发生共振,其安全可靠性尤为重要。混凝土段与钢塔段的连接段是组合塔架结构的关键部件,连接段对结构的性能有重要影响。塔架制造成本约占整机成本的20%,在满足风电机组安全性的前提下,对塔架结构进行优化,降低成本显得十分重要。
近年来,国内外学者对不同类型塔架结构的优化设计进行了诸多探索性研究,推动了优化设计技术在该领域的发展。UYS等[2]以成本为目标函数,以塔筒直径、塔筒壁厚、加劲肋数量与厚度为设计变量,对一个45 m高的传统钢塔架进行优化设计,获得了成本最低时加劲肋数量与厚度的关系。陈俊岭等[3]结合某120 m高度组合塔架,建立了以横截面几何尺寸、混凝土段高度和钢塔高度为设计变量,以整体稳定性和局部稳定性、强度、刚度和频率为约束条件,以塔架建设总成本为目标函数的优化设计模型并进行了优化,指出组合塔架混凝土段高度为塔架总高度的67%时整体造价最低。陈逸杰等[4]对轮毂高度分别为85.0 m、102.5 m和122.0 m的2.5 MW风力机组合塔架进行优化,结果表明组合塔架混凝土段高度为塔架总高度的74%时整体造价最低。许斌等[5]开展了组合塔架混凝土段高度、钢塔段高度以及截面尺寸的优化,指出当混凝土段的高度接近塔架总高度的70%时造价最低。对组合塔架而言,连接段的设计对结构的性能及建设成本均有重要影响,相关学者对此也开展了进一步的研究。许斌等[6,7]提出了混凝土填充钢箱和穿筋连接的方案,研究表明所设计方案的结构受力性能明显优于传统法兰盘连接。谢咏剑[8]对一种嵌入式开孔板和穿孔钢筋结构连接方案进行了有限元计算,综合分析指出在法兰盘下部设置嵌入式钢筒和穿孔钢筋能有效提高连接段刚度并改善应力分布,可以避免传统法兰盘连接因刚度不足而出现与混凝土分离的现象,该研究表明组合塔架连接段需要精心设计,以确保整体结构的安全。上述研究工作为组合塔架设计质量的提升提供了重要的依据,建立包含连接段厚度等参数的优化设计数学模型并进行优化设计,对提升组合塔架设计质量意义重大。
本文以某3 MW风力机组合塔架为研究对象,构建了以各塔段外径、壁厚、连接段厚度和混凝土段高度等关键几何尺寸为设计变量,以塔架固有频率、关键部位应力、塔顶位移和稳定性等关键性态指标为约束条件,以塔架成本为目标函数的优化设计数学模型,采用有限单元法进行结构计算分析,在Workbench平台上进行塔架优化设计,对结果进行详细分析。
1 组合塔架参数
3 MW组合塔架由4段组成,其中底部第Ⅰ段为混凝土段,高度为55 m,Ⅱ、Ⅲ、Ⅳ段为钢塔段,三段钢塔高度合计82.8 m,塔架总高度为137.8 m。上部三段钢塔筒之间通过法兰盘连接,底部混凝土段通过其顶部的钢垫板与第Ⅱ段钢塔筒下部的法兰盘实现连接,沿混凝土塔壁竖向设置24束预应力钢绞线,每束包含9根低松弛钢绞线（1×9φs15.2）,组合塔架结构及连接段如图1所示,各段的内外径及壁厚参数见表1。
Fig. 1 Diagram of hybrid tower图1 组合塔架示意图 |
Table 1 Design parameters of hybrid tower表1 组合塔架设计参数 |
| 塔段 | 高度 / m | 底部外径 / mm | 顶部外径 / mm | 底部厚度 / mm | 顶部厚度 / mm |
|---|---|---|---|---|---|
| IV段 | 26.88 | 4 450 | 3 280 | 24 | 28 |
| III段 | 27.28 | 4 900 | 4 450 | 32 | 24 |
| II段 | 28.64 | 5 100 | 4 900 | 42 | 32 |
| I段 | 55.00 | 7 800 | 5 100 | 330 | 315 |
Table 2 Material parameters of hybrid tower表2 组合塔架材料参数 |
| 材料 | 弹性模量 / GPa | 泊松比 | 密度 / (kg/m3) |
|---|---|---|---|
| 混凝土 | 36 | 0.2 | 2 500 |
| 钢筒、连接段 | 200 | 0.3 | 7 850 |
| 钢绞线 | 195 | 0.3 | 7 850 |
2 组合塔架优化模型
2.1 设计变量
为反映组合塔架整体特征,综合考虑后选取以下关键几何参数为设计变量：混凝土段底部外径（d1）和壁厚（t1）、混凝土段顶部外径（d2）和壁厚（t2）、混凝土段高度（h）、钢塔段底部壁厚（t3）、钢塔段顶部外径（d3）和壁厚（t4）。为了反映连接段这一关键部位,在上述设计变量的基础上增加了连接段的关键几何参数：钢垫板厚度（t5）、法兰盘厚度（t6）。参考工程经验和相关设计规范,设计变量取值范围如表3。
Table 3 Value scope of the design variables表3 设计变量取值范围 |
| 取值范围 | d1 / m | d2 / m | d3 / m | h / m | t1 / mm | t2 / mm | t3 / mm | t4 / mm | t5 / mm | t6 / mm |
|---|---|---|---|---|---|---|---|---|---|---|
| 下限 | 6 | 3 | 2 | 50 | 260 | 260 | 30 | 10 | 40 | 120 |
| 上限 | 9 | 6 | 4 | 120 | 380 | 380 | 60 | 30 | 80 | 170 |
2.2 约束条件
约束条件是组合塔架优化设计中寻找目标函数最值的多个限制条件,塔架结构设计时必须使各部件具有足够的强度、刚度和稳定性,塔架在各种荷载作用下会产生振动,应避免塔架与叶片的共振。本文优化中的约束条件如下。
（1）频率约束。为了防止塔架与叶片发生共振,一般要求塔架的固有频率与叶片转动频率及叶片通过频率都有10%以上的间隔,本文塔架的固有频率须满足下式：
$\left\{ \begin{matrix} {{f}_{0,1}}0.9{{f}_{\text{R}}} \\ 1.1{{f}_{\text{R}}}{{f}_{0,1}}0.9{{f}_{\text{R,3}}} \\ \end{matrix} \right.$ (1)
式中：f0,1为塔架的第一阶固有频率;fR为正常运行运行时叶片的最大旋转频率;fR,3为3倍的叶片旋转频率。
（2）最大位移约束。参考《高耸结构设计规范》[12]规定的塔架许用挠度应控制在塔架总高度的1/150以内,塔架最大位移约束计算公式为：
${{d}_{\max }}-\frac{H}{150}0$ (2)
式中：dmax为塔架顶部最大水平位移;$\frac{H}{150}$为塔架顶部最大容许水平位移。
（3）钢塔段应力约束。钢塔最大应力不超过设计强度。
${{\sigma }_{\text{e,}\max }}-\frac{{{f}_{\text{y}1}}}{{{\gamma }_{\text{m}}}}0$ (3)
式中：${{\sigma }_{\text{e,}\max }}$为钢塔全段内的最大等效应力;fy1为钢材的屈服极限;${{\gamma }_{\text{m}}}$为材料的安全系数。
（4）连接段应力约束。钢垫板和法兰盘的最大应力不超过其设计强度。
${{\sigma }_{\max }}-\frac{{{f}_{\text{y2}}}}{{{\gamma }_{\text{m}}}}0$ (4)
式中：$\sigma _{max}$为连接段法兰盘和钢垫板的最大等效应力;fy2为连接段钢材的屈服极限。
（5）混凝土段应力约束。最大拉压应力不超过强度设计值。
$\left\{\begin{array}{l}\sigma_{t}-f_{t} \leqslant 0 \\ \sigma_{c}-f_{c} \leqslant 0\end{array}\right.$ (5)
式中：$\sigma _{t}$为混凝土段最大拉应力;$\sigma _{c}$为混凝土段最大压应力;ft为混凝土抗拉强度设计值;fc为混凝土抗压强度设计值。
（6）稳定性约束。保证钢塔局部稳定,考虑钢塔段最大应力不超过塔筒局部稳定的临界应力值;保证塔顶在特定的外荷载作用下塔架的整体稳定性,塔架的第一阶临界屈曲荷载要大于轴向荷载的1.1倍。
$\lambda_{1}-\gamma_{s} \geq 0$ (6)
2.3 目标函数
目标函数是判别一个设计方案优劣标准的数学表达式,反映了结构中某一个最重要的特性或指标。工程结构的优化目标一般从施工技术、经济角度等方面进行考虑,比如结构的总体积最小、结构总造价最低等。本文以组合塔架总成本最低为优化目标。其计算公式如式（7）所示。
$F={{A}_{1}}{{V}_{1}}+{{C}_{1}}{{M}_{1}}+{{C}_{2}}{{M}_{2}}$ (7)
式中：F为组合塔架材料总价格;V1为混凝土段总体积;M1为连接段部件总质量;M2为钢塔段总质量;A1为混凝土单价;C1为连接段材料单价;C2为钢塔筒单价。
本文根据某风电项目材料的平均价格来确定所选材料价格,泵送C60混凝土造价为720元/m3;法兰盘和钢垫板材料价格为13 600元/t;钢塔筒材料价格为9 200元/t。
2.4 优化模型
综合上述分析,3 MW风电机组组合塔架的优化设计数学模型如下。
设计变量：X=[d1,d2,d3,h,t1,t2,t3,t4,t5,t6]T
使目标函数：F = 720V1 + 13 600 M1 + 9 200 M2 → min
满足约束条件：$\left\{ \begin{matrix} {{f}_{0,1}}0.189 \\ 0.231{{f}_{0,1}}0.567 \\ {{d}_{\max }}-0.9180 \\ {{\sigma }_{\text{e,}\max }}-313.60 \\ {{\sigma }_{\max }}-268.20 \\ {{\sigma }_{\text{t}}}-2.040 \\ {{\sigma }_{\text{c}}}-27.50 \\ {{\lambda }_{1}}-1.10 \\ \end{matrix} \right.$
3 结构有限元分析
考虑到结构的复杂性以及各部件的细部特征,本文采用有限单元法进行结构分析。
3.1 有限元模型
组合塔架结构通过Creo参数化建模并导入Ansys Workbench中。塔架简化为底端固定的结构。混凝土采用solid65实体单元,每束预应力钢绞线通过多点约束（multi-point constraints, MPC）算法将预应力钢绞线的两端约束在上下两锚具与塔架的接触面上,预应力钢绞线采用link180杆单元,预应力的加载方式采用等效荷载法,即将钢绞线的作用以荷载的形式施加于结构。塔架承受上部结构重量采用MASS21质量单元来考虑,并通过MPC184多点约束单元将质量单元和塔架连接。钢塔、法兰盘、钢垫板采用solid186实体单元。
为了保证网格的质量,网格划分时,通过切割实体避免出现不规则的几何体以确保单元质量,采用多区域划分技术。整个塔架结构剖分后的单元总数为419 110个,节点总数为353 360个,网格划分结果如图2所示。在数值模拟中,钢塔段与段之间采用绑定接触;法兰盘与钢垫板采用面面接触,法向属性采用摩擦接触,接触面之间可以传递压力而不传递拉力,接触面之间的切向属性采用“有限滑动”模型,接触面摩擦系数取0.3[13];钢垫板与混凝土顶面采用绑定接触;预应力混凝土段与段之间通过C80座浆料与预应力钢绞线连接,使其成为预应力混凝土整体结构,混凝土段与段之间采用绑定接触。通过不同网格的尺寸进行网格无关性的验证,该网格尺寸满足精度要求。
Fig. 2 Mesh of hybrid tower图2 组合塔架整体网格 |
Table 4 Load on top of tower under different conditions表4 不同工况下的塔顶荷载 |
| 工况 | ${{F}_{x}}$/ kN | ${{F}_{y}}$/ kN | ${{F}_{z}}$/ kN | ${{M}_{x}}$/ (kN∙m) | ${{M}_{y}}$/ (kN∙m) | ${{M}_{z}}$/ (kN∙m) |
|---|---|---|---|---|---|---|
| 额定风速 | 534.7 | -410.3 | -2 246 | 1 463.4 | -1 090.7 | 1 538.2 |
| 极限风速 | 1 158.7 | -348.2 | -2 490 | 3 458.6 | -2 551.3 | 1 847.1 |
3.2 初始方案的评估
Fig. 3 (a) First order natural frequency; (b) horizontal displacement of the tower under rated working conditions; (c) horizontal displacement of the tower under extreme working conditions图3 （a）一阶固有频率;（b）额定工况下塔架水平位移;（c）极限工况下塔架水平位移 |
图4~ 图7为额定工况和极限工况下各部件应力云图。钢塔段等效应力最大值出现在背风侧与法兰盘接触处,最大值为213.77 MPa,小于材料的设计强度。连接段法兰盘等效应力最大值为265.23 MPa,钢垫板等效应力最大值为111.99 MPa,连接段满足设计强度要求,连接段应力最大值出现在迎风侧预应力钢绞线孔道附近,由于此处布置了预应力钢绞线,容易产生应力集中现象,因此上述结果是合理的。额定工况下混凝土段竖向全部受压,竖向最大压应力为24.02 MPa,均小于C60混凝土强度设计值。极限工况下混凝土段绝大部分竖向受压,局部出现了较小的拉应力;竖向最大压应力为26.45 MPa,出现在背风侧顶部;最大拉应力为1.72 MPa,出现在迎风侧顶部,满足规范要求。
Fig. 4 Equivalent stress of steel tower under two working conditions图4 两种工况下钢塔等效应力 |
Fig. 5 Vertical normal stress of concrete under two working conditions图5 两种工况下混凝土竖向正应力 |
Fig. 6 Equivalent stress of flange plate under two working conditions图6 两种工况下法兰盘等效应力 |
Fig. 7 Equivalent stress of steel plate under two working conditions图7 两种工况下钢垫板等效应力 |
初始方案的计算分析表明,塔架各项性能指标满足规范要求,是一个可行的设计方案,各种性能指标均有一定的余量,存在一定的优化设计空间。
4 优化设计及结果分析
Table 5 Design variables before and after optimization表5 优化前后的设计变量 |
| 参数 | d1 / m | d2 / m | d3 / m | t1 / mm | t2 / mm | t3 / mm | t4 / mm | t5 / mm | t6 / mm | h / m |
|---|---|---|---|---|---|---|---|---|---|---|
| 初始值 | 7.80 | 5.10 | 3.28 | 330 | 315 | 50 | 140 | 42 | 28 | 55.0 |
| 优化值 | 7.67 | 4.63 | 3.15 | 322 | 303 | 43 | 138 | 45 | 26 | 85.3 |
Table 6 Key performance indicators before and after optimization表6 优化前后关键性能指标 |
| 参数 | 一阶固有 频率 / Hz | 塔顶位移 / mm | 钢塔等效应力 最大值 / MPa | 法兰盘等效应力 最大值 / MPa | 钢垫板等效应力 最大值 / MPa | 混凝土段最大 拉应力 / MPa | 混凝土段最大 压应力 / MPa | 屈曲 因子 |
|---|---|---|---|---|---|---|---|---|
| 初始值 | 0.321 | 853.68 | 213.77 | 265.23 | 111.99 | 1.72 | 26.45 | 1.85 |
| 优化值 | 0.308 | 898.59 | 273.73 | 266.91 | 159.75 | 0 | 23.93 | 1.38 |
5 优化方案疲劳校核
Fig. 7 Fatigue cumulative damage diagram of steel components图7 钢制部件疲劳累积损伤图 |
6 结论
为提高组合塔架结构的安全性能以及降低塔架成本,构建了包含连接段厚度的塔架优化设计数学模型,对某3 MW风电机组组合塔架进行了优化设计,主要结论如下：
（1）对组合塔架整体结构的优化设计切实可行,可为实际工程提供更加经济合理的方案。
（2）对组合塔架而言,混凝土段的高度占总高度的比例对整体的性能与成本有着重要影响,综合分析混凝土段高度为塔架总高度的60% ~ 75%时塔架总成本相对较低。
（3）考虑连接段厚度的优化模型能有效调整连接段的受力性能,有利于提高塔架整体优化的效果。