网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于代理模型的输电塔半刚性节点弯矩-转角曲线预测方法  PDF

  • 刘泉 1
  • 李正良 2
  • 彭思思 2
  • 王涛 2,3
1. 国网经济技术研究院有限公司,北京102209; 2. 重庆大学 土木工程学院,重庆 400045; 3. 哈尔滨工业大学 交通科学与工程学院,哈尔滨150090

中图分类号: TM754

最近更新:2024-05-06

DOI:10.11835/j.issn.1000.582X.2024.04.007

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

为了准确、高效地评估输电塔半刚性节点的力学特性,提出了一种基于代理模型的输电塔半刚性节点弯矩-转角曲线预测方法,通过引入代理模型方法近似半刚性节点几何尺寸与极限抗弯承载力、初始转动刚度之间的函数关系,建立具有较高精度的预测模型,进而结合Kish-Chen幂函数模型拟合输电塔半刚性节点的弯矩-转角曲线。结果表明,提出的基于代理模型的输电塔半刚性节点弯矩-转角曲线预测方法能减少实验和数值模拟的成本,较好地模拟输电塔半刚性节点实际受力-变形情况,为输电塔半刚性节点的工程设计和理论研究提供了参考。

近年来,特高压输电线路在我国得到了迅速发[

1]。作为电网工程中的重要基础设施,输电塔结构中各构件通过节点相互连接,由连接节点实现杆件之间的传力。输电塔节点既可以传递部分弯矩,也会产生一定的相对转角,是介于铰接与刚接之间的半刚性连接,受其构造影响,输电塔半刚性节点的力学性能表现出明显的非线性特征。输电塔节点的破坏往往会导致相连杆件的失效甚至整个结构的破坏,因此,节点的安全可靠对于整个输电塔结构至关重要。近年来,许多学者通过理论分[2]、实验研[3]和数值模[4⁃5]等方法对不同类型的输电塔连接节点展开了研究,其主要内容为:对节点极限承载力的研究以及对节点弯矩-转角(M-θ)曲线关系的研究。在输电塔节点的极限承载力方面,学者们对不同类型节点的实际极限承载力值进行测[4⁃5],并推导承载力理论计算公[2,6] ,研究了影响节点极限承载力大小的因[3]。相较于前者,关于输电塔节点弯矩-转角曲线关系的研究尚在起步阶段,李正良[6]对特高压角钢输电塔K形半刚性连接节点进行了足尺试验,考察节点的弯矩-转角曲线关系;在试验的基础上,焦安亮[7]推导了输电塔半刚性节点弯矩-转角曲线的理论计算公式;王开[8]通过建立精细化的角钢塔节点有限元模型,获取了节点对角钢的弯矩-转角曲线。总体而言,现有的输电塔半刚性节点弯矩-转角曲线的相关研究往往通过昂贵的实验或耗时的数值模拟等方法获[2,6],存在实验成本高、模拟耗时久、计算效率低等不足。

鉴于此,文中提出了一种基于代理模型的输电塔半刚性节点弯矩-转角曲线预测方法,即采用代理模型方法拟合节点几何尺寸与极限抗弯承载力、初始转动刚度之间的函数关系,并结合Kish-Chen幂函数模型拟合输电塔半刚性节点的弯矩-转角曲线。

1 代理模型方法及模型精度评价指标

1.1 代理模型方法

代理模型方法通过获取少量样本点及其响应,构建输入输出系统的近似映射关系,可用于替代复杂工程结构的响应分析计算过程,能够大幅降低计算消耗。目前,常见的代理模型方法有多项式响应面 (Polynomial response surface, PRS)模型、径向基函数 (Radial basis function, RBF)模型、克里金 (Kriging)模型、支持向量回归 (Support vector regression, SVR)模型、响应面(Response surface method, RSM)模型等。

1.1.1 PRS模型

PRS模[

9]应用了数学和统计学相关的理论知识,在已有的输入输出数据基础上,通过回归方法求出基函数中的各个多项式的权重系数,得到输入输出之间的近似数学表达。PRS模型的基本形式为

y˜(x)=β0+i=1dβixi+i1=1di2=1dβi1i2xi1xi2++i1=1di2=1dim=1dβi1i2...imxi1xi2...xim (1)

式中:d为变量的维度;m表示多项式的次数;β0, βi, βij, …, βi1i2im表示回归系数,一般可通过最小二乘法获取;xi, xi1, xi2, …, xim表示输入变量;y˜(x)表示输出响应变量;i, i1, i2, …, im表示变量具体对应的维数。

PRS模型的优势在于构造简单,建模所需样本少,能够达到较好的全局近似效果,还能对系统各个变量进行灵敏度分析,具有较好的透明性。实际工程中应用较为广泛的是线性响应面和二次响应面模型:

y˜(x)=β0+i=1dβixi (2)
y˜(x)=β0+i=1dβixi+i1=1di2=1dβi1i2xi1xi2 (3)

式中,各物理量的含义同式(1),不再赘述。如式(2)所示,线性响应面虽然待定系数和所需样本较少,但不能反映函数的非线性特征;如式(3)所示,二次响应面能够体现函数的非线性特征,计算时通常采用加权最小二乘法确定系数,但对于高维非线性问题,PRS模型的局部拟合效果较[

10]

1.1.2 RBF模型

RBF模[

11]作为一种应用广泛的非线性代理模型方法,具有结构简单、精确插值、鲁棒性好等特点。RBF模型假设输出为基函数的线性组合,采用插值方法进行逼近,通过一元函数对多变量问题进行描[12]。对于一个样本容量为n的训练样本集合X={x1, x2, …, xn}及其响应Y={y1, y2, …, yn},RBF模型的一般表达式为

y˜x=i=1nλiϕx-xi+j=1mβjPjx (4)

式中:λi表示权重系数,可通过求解线性插值方程组获得; ||∙||表示用于计算向量的欧氏距离的二范数;ϕ(·)为基函数,其常用的形式包括线性函数、高斯函数、二次函数、逆二次函数、薄板样条函数等,文中选用的二次函数为基函数为

ϕr=c2+r2 (5)

式中:r表示样本点的欧氏距离;c为实常数,文中取c=1。

式(4)中第二项表示m阶的多项式回归,其中,βj表示回归系数,Pj(x)表示多项式且满足如下正交条件:

i=1nλiPjxi=0 (6)

因此,系数λβ可通过求解下列线性方程组得到:

ΦPTP0λβ=Y0 (7)

式中:Φn阶方阵,Φ(i, j)= ϕ(||xi-xj||); Pn×m阶矩阵,P(i, j)=pj(xi)。

1.1.3 Kriging模型

Kriging模[

13]是一种基于贝叶斯数学模型的线性无偏估计模型,显著特点在于其数学表达式中包括了线性回归参数部分和非参数部分,即:

y˜(x)=G(β,x)+Z(x) (8)

式中:G(β, x)被称为全局近似函数,为Kriging模型的数学期望;β为回归系数;Z(x)为一个均值为零、方差为σ2的静态随机过程,其协方差不等于零,可表示为

Cov(Z(xi),Z(xj))=σ2R(xi,xj) (9)

式中,R(xi, xj)为空间相关函数,其相关性随着距离增大而减小,当距离为0时,等于1;当距离无穷大时,等于0。目前,常用的空间相关函数包括高斯函数、线性函数、指数函数等,文中采用的高斯相关函数为

R(xi,xj)=exp(-k=1dθkxik-xjkδk) (10)

式中:d表示设计变量的维数;xikxjk分别为任意2个样本xixj在第k个维度上的分量;θk为相关性参数;δk为光滑程度参数。

1.2 模型精度评价指标

为评估代理模型的全局及局部近似程度,采用4种通用的精度评价指标来衡量代理模型的精度,包括确定性系数R2、标准均方根误差NRMSE、相对平均绝对误差RAAE、相对最大绝对误差RMAE,计算公式分别为

R2=1-i=1Nyxi-y˜xi2i=1Nyxi-y¯xi2 (11)
NRMSE=1ymax-ymini=1Nyxi-y˜xiN (12)
RAAE=i=1Nyxi-y˜xiN×STD (13)
RMAE=maxyxi-y˜xiSTD (14)

式中:N为测试样本点数量;y(xi)是测试样本点所对应的函数输出结果;ymaxymin分别为函数最大值、最小值;y¯(xi)是函数输出响应的平均值;(xi)为模型近似值;STD表示实际函数值的标准差。

上述模型精度指标中,R2整体反映了代理模型的精度,当值越接近1,模型越精确;NRMSE和RAAE为全局精度指标,其值越小,模型精度越高;RMAE则是1个局部指标,描述的是设计空间内某个局部区域的误差,故其值越小越好。

2 Kish-Chen模型

国内外关于半刚性连接节点弯矩-转角曲线的研究主要集中于框架结构领域,学者们提出了线性模型、多项式模型、B样条模型、幂函数模型、指数模型等拟合半刚性节点的弯矩-转角(M-θ)曲线。其中,Kishi[

14⁃15]在大量实验数据基础上,对Lui[16]、Ang[17]建议的多参数模型进行了改进,提出了Kish-Chen三参数幂函数模型,模型因形式直观、易于实现且具有明确的物理含义而颇受研究者青[18],是应用最为广泛的半刚性节点M-θ曲线预测模型之一。

Kish-Chen幂函数模型使用初始刚度、极限弯矩承载力、形状系数3个参数描述梁柱节点的半刚性,其数学表达式为

M=Kiθ1+(θθ0)ω1ω (15)

式中:Ki为初始转动刚度;ω是曲线拟合得到的形状系数;θ0为相对塑性转角,可通过极限抗弯承载力Mu和初始转动刚度Ki计算得到:

θ0=MuKi (16)

文中采用Kish-Chen模型拟合输电塔半刚性节点的M-θ曲线,由式(16)可知,输电塔半刚性节点的极限抗弯承载力Mu和初始转动刚度Ki是建立M-θ曲线预测模型的关键。

3 输电塔半刚性节点弯矩-转角曲线预测方法

输电塔半刚性节点的极限抗弯承载力Mu和初始转动刚度Ki与节点的几何尺寸密切相关,一般通过实验测试或有限元分析计算得到。由于实验研究成本高且资源相对有限,精细化数值模型计算较为耗时,连接节点的构造又较为复杂,完全依赖实验或者数值模拟对半刚性节点进行参数分析难以实现。因此,提出一种更加简单易行、精度高的方法对于深入研究输电塔半刚性节点几何尺寸与其力学性能之间的关系具有重要现实意义。代理模型方[

19]是解决这类问题的有效途径之一。为进一步降低成本、提高计算效率,文中采用代理模型方法拟合输电塔半刚性节点几何尺寸与极限抗弯承载力Mu、初始转动刚度Ki之间的近似函数关系,并将预测的ui带入Kish-Chen模型中拟合输电塔半刚性节点的M-θ曲线,发展了基于代理模型的输电塔半刚性节点弯矩-转角曲线预测方法。其步骤如下:

1) 选取输电塔半刚性节点几何尺寸作为设计变量,抽取n个训练样本点,通过实验或数值分析模型计算相应的极限抗弯承载力Mu、初始转动刚度Ki响应值;

2) 基于已知训练样本点及其响应,分别建立输电塔半刚性节点MuKi的PRS、Kriging、RBF预测模型;

3) 选取N个测试样本点及其响应值,并根据式(11)~式(14)计算PRS、Kriging、RBF模型的精度;

4) 选择当前精度最高的预测模型,输入某一输电塔半刚性节点的几何尺寸,得到相应的输电塔半刚性节点的极限抗弯承载力u和初始转动刚度i

5) 将代理模型预测结果ui带入式(15)~式(16)所示Kish-Chen幂函数模型,拟合输电塔半刚性节点M-θ关系曲线。

4 算例分析

为验证文中提出的基于代理模型的输电塔半刚性节点弯矩-转角曲线预测方法,以参考文献[

20]中输电塔K形连接节点为例,构造如图1所示,根据文中方法计算其弯矩-转角曲线关系。

图1  输电塔K形连接节[

20]

Fig. 1  K-joint of transmission towers[

20]

文献[

20]中给出了此类节点参数数据,选取其中钢材屈服强度fy=420 MPa、螺栓等级强度ng=10.9级、螺栓个数nl=4、节点板厚tj=10的41组数据,将其分为用于构建代理模型的训练样本以及用于检验模型精度的测试样本,分别如表1表2所示,共涉及3个节点几何尺寸设计变量:肢宽(b)、肢厚(t)、螺栓直径(d)。

表1  训练样本[20]
Table 1  The set of training samples[20]
b/mmt/mmd/mmMu//(kN·m)Ki//(kN·m·rad-1)
140 10 20 72.63 3 250
140 12 20 72.49 3 320
160 10 20 71.91 3 300
160 14 20 71.40 3 410
160 16 20 71.16 3 450
180 14 20 71.94 3 450
180 16 20 71.92 3 480
180 18 20 72.28 3 510
200 14 20 71.25 3 480
200 16 20 71.58 3 520
200 20 20 70.72 3 540
200 24 20 70.82 3 510
160 12 16 59.54 3 200
160 12 22 73.91 3 410
160 12 24 74.32 3 440
160 12 27 74.63 3 460
160 12 30 77.47 3 550
140 10 16 63.87 3 140
140 10 18 64.83 3 180
140 10 22 73.38 3 283
140 10 24 73.59 3 285
140 10 27 73.37 3 309
140 10 30 74.96 3 376
180 14 18 67.00 3 365
180 14 22 74.19 3 511
180 14 24 75.60 3 529
200 16 16 62.38 3 351
200 16 18 66.99 3 412
200 16 22 73.58 3 568
200 16 24 74.14 3 600
200 16 27 74.59 3 647
表2  测试样本[20]
Table 2  The set of test samples[20]
b/mmt/mmd/mmMu/(kN·m)Ki/(kN·m·rad-1)
160 12 18 69.14 3 300
140 13 20 71.15 3 390
180 12 20 72.56 3 410
180 14 30 77.47 3 664
160 12 20 71.88 3 370
180 14 16 62.84 3 294
200 16 30 76.87 3 758
180 14 27 75.68 3 593
140 12 20 71.87 3 360
200 18 20 70.86 3 534
160 12 18 69.14 3 300

根据表1所示的训练样本集合,分别采用PRS、Kriging、RBF 3种代理模型方法建立输电塔半刚性节点极限抗弯承载力Mu、初始转动刚度Ki的预测模型,并依据表2中测试样本集合比较当前代理模型精度。

根据式(11)~式(14)计算3种不同代理模型方法的建模精度,结果如表3所示。由表3可知,采用PRS模型方法建立的输电塔半刚性节点极限抗弯承载力Mu预测模型的整体精度和局部精度优于其他方法;然而,对于输电塔半刚性节点初始转动刚度Ki,使用RBF模型所得结果更佳,精度更高。

表3  代理模型精度
Table 3  The accuracy of surrogate models
代理模型方法PRSKrigingRBF
Mu/(kN·m)Ki/(kN·m·rad-1)Mu/(kN·m)Ki/(kN·m·rad-1)Mu/(kN·m)Ki/(kN·m·rad-1)
R2 0.874 0.961 0.712 0.697 0.773 0.967
NRMSE 0.098 0.065 0.147 0.180 0.131 0.060
RAAE 0.286 0.162 0.398 0.381 0.380 0.139
RMAE 0.649 0.340 1.021 1.345 0.852 0.331

文中分别采用PRS和RBF模型方法建立输电塔K形连接节点的极限抗弯承载力Mu和初始转动刚度Ki的预测模型,进一步将计算结果带入式(15)~式(16)所示的Kish-Chen幂函数模型,最终得到输电塔半刚性节点的近似M-θ关系曲线,将部分节点的拟合结果与文献[

20]中有限元模拟的M-θ关系实际曲线进行比较,如图2所示。

图2  M-θ关系曲线对比

Fig. 2  Comparison of M-θ relationship

图2分别为2种不同规格的输电塔半刚性K形节点弯矩-转角关系曲线,节点的tj =10 mm,d=20 mm,模型中形状系数ω的取值为文献[

20]中建议的平均值。由图2可知,文中基于代理模型以及Kish-Chen模型的预测方法拟合的输电塔半刚性K形节点M-θ关系曲线与实际曲线整体吻合效果良好,表明文中方法具有一定的可靠性,适用于输电塔半刚性节点的弯矩-转角曲线计算。

5 结 论

文中提出了一种普适的输电塔半刚性节点弯矩-转角关系预测方法,该方法基于代理模型建立输电塔半刚性节点极限抗弯承载力和初始转动刚度的预测模型,并结合Kish-Chen模型拟合输电塔半刚性节点的弯矩-转角曲线关系。结果表明,使用PRS和RBF代理模型方法所得输电塔半刚性K形节点几何尺寸与极限抗弯承载力初始转动刚度的预测模型精度较高,最终拟合的输电塔半刚性节点弯矩-转角曲线效果良好。文中方法通过引入代理模型能大幅减少计算消耗,同时能够有效预测输电塔半刚性连接节点的受力-变形性能,广泛适用于各类输电塔半刚性节点的弯矩-转角关系计算,为输电塔节点受力分析及设计提供了参考。

参考文献

1

俞登科, 李正良, 施菁华, . ±800 kV直流双柱悬索拉线塔塔线体系风振响应的风洞试验研究[J]. 中国电机工程学报, 2015, 35(4): 1009-1013. [百度学术] 

Yu D K, Li Z L, Shi J H, et al. Wind tunnel test on wind-induced response of ±800 kV DC cross-rope suspension tower-line[J]. Proceedings of the CSEE, 2015, 35(4): 1009-1013.(in Chinese) [百度学术] 

2

余世策, 孙炳楠, 叶尹, . 高耸钢管塔节点极限承载力的试验研究与理论分析[J]. 工程力学, 2004, 21(3): 155-161. [百度学术] 

Yu S C, Sun B N, Ye Y, et al. Experimental study and theoretical analysis of ultimate strength for steel tubular joint of tall towers[J]. Engineering Mechanics, 2004, 21(3): 155-161.(in Chinese) [百度学术] 

3

李正良, 刘红军. 输电塔钢管-插板连接节点板承载力研究[J]. 土木工程学报, 2011, 44(S1): 52-58. [百度学术] 

Li Z L, Liu H J. Study on bearing capacity of steel tube-plug connection joint plate of transmission tower[J]. China Civil Engineering Journal, 2011, 44(S1): 52-58.(in Chinese) [百度学术] 

4

邓洪洲, 朱雯瑞, 黄斌, . 输电钢管塔典型塔脚节点承载力研究[J]. 工程力学, 2017, 34(9): 175-183. [百度学术] 

Deng H Z, Zhu W R, Huang B, et al. Study on the capacity of typical tower leg joint of transmission steel tubular tower[J]. Engineering Mechanics, 2017, 34(9): 175-183.(in Chinese) [百度学术] 

5

邓洪洲, 姜琦, 黄誉. 输电钢管塔K型管板节点承载力试验及参数[J]. 同济大学学报(自然科学版), 2014, 42(2): 226-231, 314. [百度学术] 

Deng H Z, Jiang Q, Huang Y. Experimental investigation and parametric studies on ultimate strength of K-type tubegusset plate connections in transmission towers[J]. Journal of Tongji University (Natural Science), 2014, 42(2): 226-231, 314.(in Chinese) [百度学术] 

6

李正良, 赵楠, 尤军, . 角钢塔K形节点弯矩-转角曲线研究[J]. 四川大学学报(工程科学版), 2013, 45(4): 46-54. [百度学术] 

Li Z L, Zhao N, You J, et al. Development of moment-rotation curves of angle steel tower K-joints[J]. Journal of Sichuan University (Engineering Science Edition), 2013, 45(4): 46-54.(in Chinese) [百度学术] 

7

焦安亮, 李正良, 刘红军, . 特高压输电塔半刚性连接K形节点受力性能研究[J]. 建筑结构学报, 2014, 35(7): 53-60. [百度学术] 

Jiao A L, Li Z L, Liu H J, et al. Study on ultimate strength of semi-rigid K-type joints in UHV transmission tower[J]. Journal of Building Structures, 2014, 35(7): 53-60.(in Chinese) [百度学术] 

8

王开源. 输电塔节点半刚性及角钢稳定承载力分析[D]. 哈尔滨: 哈尔滨工业大学, 2016. [百度学术] 

Wang K Y. The analysis of semi-rigid joint and single angle stability bearing capacity of transmission tower[D]. Harbin: Harbin Institute of Technology, 2016. (in Chinese). [百度学术] 

9

Booker A J, Dennis J E, Frank P D, et al. A rigorous framework for optimization of expensive functions by surrogates[J]. Structural Optimization, 1999, 17(1): 1-13. [百度学术] 

10

Barton R R. Metamodels for simulation input-output relations[C]//Proceedings of the 24th conference on Winter simulation. December 13 - 16, 1992, Arlington, Virginia, USA. ACM, 1992: 289-299. [百度学术] 

11

Dyn N, Levin D, Rippa S. Numerical procedures for surface fitting of scattered data by radial functions[J]. SIAM Journal on Scientific and Statistical Computing, 1986, 7(2): 639-659. [百度学术] 

12

胡凡, 武泽平, 王东辉, . 序列近似优化方法[J]. 国防科技大学学报, 2017, 39(1): 92-101. [百度学术] 

Hu F, Wu Z P, Wang D H, et al. Sequential approximate optimization method[J]. Journal of National University of Defense Technology, 2017, 39(1): 92-101.(in Chinese) [百度学术] 

13

Kleijnen J P C. Kriging metamodeling in simulation: a review[J]. European Journal of Operational Research, 2009, 192(3): 707-716. [百度学术] 

14

Kishi N , Matsuoka K I , Chen W F ,et al. Moment-rotation relation of top and seat-angle connections[J]. Memoirs of the Muroran Institute of Technology, 1987, 37: 163-172. [百度学术] 

15

Kishi N, Chen W F. Moment-rotation semi-rigid connections[R]. Structural Engineering Report No. CE -STR-87-29. West Lafayette (IN): Purdue University, 1987. [百度学术] 

16

Lui E M, Chen W F. Analysis and behaviour of flexibly-jointed frames[J]. Engineering Structures, 1986, 8(2): 107-118. [百度学术] 

17

Ang K M, Morris G A. Analysis of three-dimensional frames with flexible beam–column connections[J]. Canadian Journal of Civil Engineering, 1984, 11(2): 245-254. [百度学术] 

18

Abolmaali A, Matthys J H, Farooqi M, et al. Development of moment–rotation model equations for flush end-plate connections[J]. Journal of Constructional Steel Research, 2005, 61(12): 1595-1612. [百度学术] 

19

Jin R, Chen W, Simpson T W. Comparative studies of metamodelling techniques under multiple modelling criteria[J]. Structural and Multidisciplinary Optimization, 2001, 23(1): 1-13. [百度学术] 

20

赵楠. 特高压输电塔半刚性连接节点受力性能及结构非线性分析研究[D]. 重庆: 重庆大学, 2014. [百度学术] 

Zhao N. Study on load-bearing performances of semi-rigid joints and structure nonlinear analysis of UHV transmission tower[D].Chongqing: Chongqing University, 2014. (in Chinese) [百度学术]