1 000 kV特高压变电构架风荷载特性

2021-03-29李方慧支旭东

哈尔滨工业大学学报 2021年4期

李方慧, 唐浩, 支旭东

(1.黑龙江大学建筑工程学院，哈尔滨150086；2.结构工程灾变与控制教育部重点实验室(哈尔滨工业大学)，哈尔滨150090；3.土木工程智能防灾减灾工业和信息化部重点实验室(哈尔滨工业大学)，哈尔滨150090)

1 000 kV特高压变电构架是十分重要的生命线工程，输电能力强，效率高，枣庄等一批1 000 kV输变电工程逐步兴建.由于其轻质、高柔、阻尼小的特性，属于典型的风敏感结构。在强风和地震作用下，其结构会产生剧烈的振动，引起杆件变形或断裂，严重的会导致整个结构倒塌[1]。目前，DL/T 5154—2012《架空输电线路杆塔结构设计技术规定》[2]等规范尚未对1 000 kV特高压变电构架风荷载给出明确的计算条款。因此，开展该类结构的风荷载特性研究具有重要意义。

国内外的一些专家学者对变电构架抗风关键问题进行了一系列研究。文献[3-8]利用气弹模型风洞试验、高频天平测力试验和数值模拟等手段对500 kV输电塔架风荷载体型系数、风振响应以及塔线体系耦合振动效应进行了全面而细致的研究，形成了较为成熟的计算理论和实用设计方法。相比而言，1 000 kV输电构架塔身较高、刚度沿高度分布不均匀，结构抗风分析更加复杂。文献[9-11]对1 000 kV特高压变电构架风荷载分布特性及等效静风荷载进行研究。文献[12]对输电线塔架顺风向荷载效应进行研究。文献[13-16]应用高频测力天平技术对800 kV和1 100 kV特高压输电塔架体型系数、风振响应等开展研究。此外，输电构架的塔-线气动阻尼、结构疲劳以及下击暴流作用[17]等方面也是研究热点，本文不展开叙述。

上述研究主要集中在500 kV和800 kV特高压变电构架风工程研究，而对1 000 kV特高压输电塔架风荷载分布特性的研究正在逐步展开，目前尚缺乏系统的研究。本文采用3D打印技术制作1 000 kV特高压变电构架整体和节段的风洞试验刚性缩尺模型，在均匀流、A类和B类三种地貌下完成高频天平测力风洞试验，利用采集的气动力系数时程数据详细考察1 000 kV变电构架整体以及各节段结构在不同风向、地貌类型下气动力系数均值、均方根值、峰度和偏度等统计特性及功率谱变化规律并分析作用机理，从而指导工程实践。

1 风洞试验

1.1 试验设备

高频底部天平测力风洞试验在哈尔滨工业大学风洞与浪槽联合实验室完成，小试验段截面宽4.0 m，高3.0 m，长25 m，配有自动转盘系统，转盘直径2.4 m，可实现0°～360°任意角度调节，风洞气动轮廓如图1所示。高频测力天平试验采用ATI六分力传感器(图2)，型号为ATI delta ip68 si-660-60，采样频率为1 000 Hz，每个样本采样时间60 s。试验中采用固定于自动上下移动支架上的Cobra Probe眼镜蛇三维脉动风速仪测量风速剖面及参考风速。

图1 风洞气动轮廓Fig.1 Aerodynamic profile of wind tunnel

图2 ATI六分力传感器Fig.2 ATI six-axis force sensor

1.2 试验模型

1 000 kV特高压变电构架为空间网格结构，由左右两个独立塔架和中间横梁组成，杆件截面形式为圆钢管，原型结构轮廓尺寸为高70 m，宽49 m(图3)。为了详细考察不同高度塔架的风荷载变化特性分别在均匀流场、A类地貌和B类地貌三种流场条件下对整体模型和部分节段模型进行高频底部测力风洞试验，将左侧塔架沿不同高度分别A、B、C、D四个节段模型(图4)。风洞试验刚性缩尺整体模型和节段模型采用SLA光敏树脂3D打印而成，缩尺比分别为1∶100和1∶50。风洞试验中整体模型及A节段模型见图5。

图3 变电构架结构Fig.3 Structure of substation frame

图4 分段示意(m)Fig.4 Segmentation diagram (m)

图5 风洞试验模型Fig.5 Wind tunnel test models

1.3 风场模拟与试验工况

采用尖塔和粗糙元、地毯对A类和B类地貌进行风场模拟，为确保风洞流场品质，将风洞试验段中风剖面和湍流度曲线与荷载规范分别进行对比分析，结果显示风场模拟效果良好。为详细考察地貌和风向对气动力系数的影响，利用结构的对称性，在0°～90°范围内每间隔15°测量一次。由于篇幅限制，以A节段模型为例详细给出风洞测量气动力系数对比分析结果。

1.4 主要研究内容

利用高频底部测力天平同步采集的力和力矩数据，详细考察x向和y向的气动力系数Cx和Cy以及扭矩系数Cmz，利用如下计算获得。

式中Fx、Fy、Mz分别为x、y方向的力和扭矩。U为参考风速，均匀流场U=14.1 m/s，A类和B类地貌参考风速取模型顶部风速，分别为7.6 m/s和8.7 m/s。ρ为空气密度，S为迎风面积，B为结构特征宽度(表1)。

表1 各风向角下A节段模型及整体模型特征宽度Tab.1 Characteristic width of segment model A and overall model at different wind directions

2 A节段模型气动力系数分析

将A节段模型气动力系数均值、均方根、偏度和峰度以及功率谱等进行对比分析，详细考察风向、地貌等因素对气动力系数的影响规律，通过每个节段模型气动力特性对比分析，可为整体结构风荷载特性的精细化研究提供依据。

2.1 统计特征分析

对比分析不同地貌、风向下的气动力系数的均值、均方根、峰度(kurtosis)和偏度(skewness)等统计特性。三种不同流场下A节段模型Cx和Cy的均值和均方根对比见图6，在0°～90°风向角下，Cx均值逐渐递增，在90°达到最大，Cx均方根值先减后增，在0°达到最大，75°取到最小值。Cy均值逐渐递减，在0°达到最大，Cy均方根逐渐递增，在90°达到最大，主要原因是随着风向角的变化，挡风面积、气体绕流和分离方式不同导致结构的受力状态发生变化。0°风向均匀流场、A类地貌和B类地貌A节段模型Cx均方根之比为1∶1.60∶1.59，Cy均方根之比为1∶1.53∶1.30。均匀流场中各风向下Cx最大均值和Cy最大均值之比为1∶1.10，说明两个垂直方向气动力特性较为接近。

图6 A节段气动力系数均值与均方根值对比Fig.6 Comparison of mean values and RMS values of aerodynamic coefficients of segment A

图7为三种流场条件下A节段模型气动力系数skewness和kurtosis对比，三类流场条件下Cx的skewness波动范围-0.38～-0.01，B类地貌下各个风向下skewness绝对值最大。同时可见，Cy的skewness波动范围0.01～0.55，在30°向角下B类地貌的skewness最大。kurtosis的变化规律见图7(b)，kurtosis值均在3刻度线以上。均匀流、A类地貌、B类地貌三种流场条件下Cx的skewness最大值之比和kurtosis最大值之比分别为1∶1.92∶3.17和1∶1.08∶1.19，同时Cy的skewness最大值之比和kurtosis最大值之比分别为1∶3.57∶7.86和1∶1.11∶1.20。

图7 A节段气动力系数偏度与峰度对比Fig.7 Comparison of skewness and kurtosis of aerodynamic coefficients of segment A

2.2 功率谱分析

从频域角度详细考察地貌、风向对气动力系数的影响规律，图8为0°风向三类流场下A节段模型气动力系数功率谱对比分析。

1) 地貌的影响

图8对比可见，0°风向A节段模型Cx和Cy功率谱在A类地貌和B类地貌高频部分大于均匀流场，这是由于紊流场对流动量高，耗散动量低。

图8 A节段气动力系数功率谱密度地貌对比Fig.8 Comparison of PSD of aerodynamic coefficients of segment A under different terrains

2) 风向角的影响

图9为输电构架A节段模型在均匀流场不同风向角时气动力系数Cx和Cy功率谱对比，在低频部分，对Cx和Cy功率谱密度影响最大的风向角分别是0°和90°。在高频部分，对Cx功率谱密度影响最大的风向角为0°和75°，对Cy功率谱密度影响最大的为0°和90°。主要原因是随着风向角的不同，气体绕流模式发生变化导致流场特性改变，能量传递随之改变。

图9 A节段气动力系数功率谱密度对比Fig.9 Comparison of PSD of aerodynamic coefficients of segment A

3 整体和节段模型气动力系数对比分析

3.1 统计特征对比

3.1.1 均值与均方根

均匀流场不同风向时各节段模型与整体模型气动力系数Cx和Cy均值对比见图10。对比发现，Cx从0°到90°递减，0°时最大，90°时最小，Cy的变化规律相反。90°风向均匀流场下整体模型与四个节段模型Cx的均值之比为1∶0.02∶0.26∶0.22∶0.38，Cy的均值之比为1∶0.46∶0.57∶0.43∶0.69。对比发现各节段模型均值D节段占比最大。

图10 整体模型及节段模型气动力系数均值对比Fig.10 Comparison of mean values of aerodynamic coefficients of overall model and segment models

图11为均匀流场下节段及整体模型气动力系数Cx和Cy均方根对比。整体模型节段模型气动力系数均方根变化趋势基本相同，Cx从0°到90°递减，0°时最大，90°时最小，Cy的变化规律正好相反。90°风向均匀流场下整体模型与A～D节段模型的Cx和Cy均方根之比分别为1∶0.06∶0.37∶0.32∶0.41和1∶0.47∶0.58∶0.44∶0.70。

图11 整体模型及节段模型气动力系数均方根对比Fig.11 Comparison of RMS values of aerodynamic coefficients of overall model and segment models

3.1.2 偏度与峰度

均匀流场输电构架节段模型及整体模型气动力系数Cx和Cy偏度对比见图12。整体模型与节段模型skewness值变化范围为-0.10～0.09，接近于0，表明气动力系数的数据分布接近对称分布。对Cx而言，在0°～90°风向角下，A节段和B节段模型skewness值均在0刻度线以下。对Cy而言，只有D节段的skewness值在0刻度线以下。

图12 整体模型及节段模型气动力系数偏度对比Fig.12 Comparison of skewness of aerodynamic coefficients of overall model and segment models

图13为均匀流场下节段及整体模型气动力系数Cx和Cy峰度对比。整体模型与节段模型气动力系数Cx的kurtosis值变化范围为2.90～3.44，整体模型在75°风向kurtosis达到3.44。对Cy而言，A节段、B节段和D节段的kurtosis值在3刻度线以上，D节段模型45°风向kurtosis达到3.23。