自复位摩擦耗能支撑模型参数的识别

2021-09-26孟丽岩刘吉胜许国山

黑龙江科技大学学报 2021年5期

王涛，李勐，孟丽岩，刘吉胜，许国山

(1.黑龙江科技大学建筑工程学院，哈尔滨 150022， 2.哈尔滨工业大学土木工程学院，哈尔滨 150090)

0 引言

自2003年，Bruneau 等[1]首次提出“可恢复功能抗震结构”的概念以来，具有恢复功能的防震结构已成世界地震工程领域的研究热点。自复位摩擦耗能支撑(Self-centering energy dissipation braces, SCED)是一种典型的具有可恢复功能结构，具有良好的变形能力、能量耗散能力及恢复能力[2]。近年来，研究者们陆续采用了一系列滞回模型描述自复位摩擦耗能支撑的力学性能。2011年，Ma等[3]使用Bouc-Wen模型[4-5]研究自复位阻尼器响应，该模型可以合理地描述支撑力-位移滞回特性。文献[6-7]基于预压弹簧自复位耗能支撑工作原理和原Bouc-Wen 模型，提出了分段简化恢复力模型，实验验证该模型能准确描述各阶段的滞回响应。2016年，Zhou等[8]在流变学模型的基础上，提出了一个预测自复位屈曲约束支撑滞回行为的弹塑性模型。樊晓伟等[9]提出了由双Bouc-Wen 模型和一组线弹性模型组成的支撑非线性原理模型，通过模拟确定了模型参数。2020年，Zhang等[10]提出了经改进的旗形模型(Modified flag-shaped model，MFS)，MFS模型包括线性部分、双线性部分和具有滑移区的弹塑性部分，可以更有效反应自复位摩擦耗能支撑的滞回响应特性。

为了能够准确给出模型参数影响规律，确定模型参数取值，笔者分析MFS模型参数的敏感性，给出模型参数合理取值范围，采用CKF算法识别MFS模型在线参数，验证该算法模型参数识别精度与计算效率。

1 自复位摩擦耗能支撑系统

1.1 SCED支撑滞回模型

SCED支撑作为一种耗能构件，具有良好的耗能效果和结构残余位移控制等优点。其中，MFS模型可以准确描述SCED支撑滞回特性，MFS模型滞回曲线与SCED支撑滞回曲线吻合程度较高，更接近SCED支撑的滞回现象，有效扩展耗能构件在抗震领域的应用。以MFS[10]模型为研究对象，MFS模型由线性部分、双线性弹性部分、具有可滑移弹塑性部分组成，自复位摩擦耗能支撑机理如图1所示。初始阶段刚度K为模型线性、双线性、具有可滑移弹塑性三部分的刚度之和，即K=K1+K2+K3。当摩擦装置启动后，支撑刚度减小至预应力筋刚度K1，卸载初始阶段的卸载刚度等于K1+K3，当力到达激活点b时，刚度变成K1+K2。

图1 自复位摩擦耗能支撑机理Fig. 1 Self-centering energy dissipation brace mechanics

MFS的恢复力可表达为

F=K1x+K2R(x)+K3z，

式中，R(x)——双线性弹性模型恢复力，kN。

R(x)具体表达式为

R(x)=x[1-H(x-b)-H(-x-b)]+

b[H(x-b)-H(x-b)]，

(1)

式中：H(x)——海维赛德阶跃函数；

B——激活位移，mm；

z——模型滞回位移，mm。

z微分表达式为

(2)

K1、K2、K3由式(3)～(5)给出

K1=AK，

(3)

K2=(1-A)(1-Q)K，

(4)

K3=(1-A)QK，

(5)

式中：A——激活刚度比，A=(1-Q)b；

Q——耗能率，是能量耗散的重要指标，Q=K3/(K2+K3)。

1.2 MFS模型有效性验证

为验证本文模拟MFS模型有效性，对MFS模型采用低周期往复位移加载，加载峰值为10 mm，加载步数为3 200步，模型加载制度如图2所示。MFS模型真实参数依次设置为：激活刚度比A=0.015、耗能率Q=0.9、激活位移b=1.4 mm、初始刚度K=300 kN/mm。MFS模型的滞回位移利用四阶Runge-Kutta积分算法得出，计算步长为0.01 s。模拟MFS模型线性部分、双折线部分、具有滑移区弹塑性部分及整体部分模拟结果如图3所示。

图2 加载制度Fig. 2 Loading protocol

由图3可以看出，通过对MFS模型进行低周期往复位移加载，MFS模型线性部分为一条斜直线、模型双线性部分为双折线、模型具有滑移区弹塑性部分及模型整体滞回曲线可以产生一个旗形滞回曲线，得出的计算结果可以反映MFS模型滞回特点，验证了所模拟MFS模型的有效性。

图3 MFS模型滞回曲线Fig. 3 Hysteresis curve of MFS model

2 MFS参数敏感性分析

在验证MFS模型的有效性的基础上，分别选取MFS模型参数激活点b、初始刚度K、模型耗能率Q、模型激活后刚度比A四个参数作为分析对象，以1.2节参数值作为参数基值，通过MFS模型低周期往复加载分别对四个参数分别进行敏感性分析，分析结构如图4～7所示。给出四个参数在定量处理后的变化规律，从而得到参数取值范围。

图4 通过对b定量分析模型各部分及整体滞回曲线Fig. 4 Quantitative analysis of each part and overall hysteretic curve of b

由图4可见，以b=1.4 mm为基值，单一调整参数b值大小，b的取值范围为0～3.5 mm，分别研究MFS模型线性部分、双线性部分、具有滑移区弹塑性部分、整体部分恢复力随b值变化规律。

由图4可以看出，随着激活点b值的增加，MFS模型线性部分没有影响，MFS模型双线性部分、MFS模型具有滑移区弹塑性部分、MFS模型整体的轴向恢复力不断增大；随着支撑轴向力的增大，支撑整体滞回耗能不断增加。

随着激活点b值的增大，MFS模型具有滑移区弹塑性部分，激活位移、激活力逐渐增大，当b=0时，整体滞回曲线近似为线性，耗能几乎为0，当b=3.5 mm时，推迟了模型进入耗能阶段，导致滞回环与真实滞回环卸载路径产生明显差异。因此，建议b的取值范围为0

由图5可以看出，针对K进行敏感性分析。以K=300 kN/mm为基值，单一调整参数K值大小，K的取值范围为0～750 kN/mm，分别研究MFS模型线性部分、双线性部分、具有滑移区弹塑性部分、整体部分恢复力随K值变化规律。

图5 通过对K定量分析模型各部分及整体滞回曲线Fig. 5 Quantitative analysis of each part and overall hysteretic curve of K

随着初始刚度K值的增加，MFS模型线性部分的图像斜率增大，MFS模型双线性部分、MFS模型具有滑移区弹塑性部分、MFS模型整体的轴向恢复力不断增大；初始刚度K与MFS模型轴向力存在正相关关系。因此，建议K的取值范围为K>0。

由图6可以看出，以Q=0.9为基值，单一调整参数Q值大小，Q的取值范围为0～1，分别研究MFS模型线性部分、双线性部分、具有滑移区弹塑性部分、整体部分恢复力随Q值变化规律

图6 通过对Q定量分析模型各部分及整体滞回曲线Fig. 6 Quantitative analysis of each part and overall hysteretic curve of Q

随着耗能率Q值的增加，MFS模型线性部分没有影响，MFS模型双线性部分、具有滑移区弹塑性部分、整体的轴向恢复力不断增大；当Q增加至1时，滞回曲线出现在二、四象限，模型没有意义。因此，建议Q的取值范围为0

由图7可以看出，以A=0.015为基值，单一调整参数A值大小，A的取值范围为0～1，分别研究MFS模型线性部分、双线性部分、具有滑移区弹塑性部分、整体部分恢复力随A值变化规律。

图7 通过对A定量分析模型各部分及整体滞回曲线Fig. 7 Quantitative analysis of each part and the overall hysteretic curve of A

随着激活点A值的增加，MFS模型线性部分的图像斜率增加，MFS模型双线性部分、具有滑移区弹塑性部分、整体的轴向恢复力不断增大；MFS模型双线性部分、具有滑移区弹塑性部分、整体的轴向力会随着耗能率A值的增加而增加，当A增加至1时，滞回曲线变成一条直线，A>1时，模型没有意义。因此，建议A值取值范围为0

3 CKF参数识别

3.1 CKF识别更新MFS模型参数

(6)

式中：x——加载位移；

H(x)——海维赛德阶跃函数。

y=KxA+(1-A)(1-Q)KR(x)+

(1-A)QKx1+WR。

(7)

V=10-25I3，

(8)

式中，I3——3×3的单位矩阵。

W=10-11，

(9)

(10)

(11)

本次模拟加载速度由式(12)近似计算得到

(12)

式中，dt——积分步长，其值为0.01 s。

3.2 仿真结果与分析

通过CKF对MFS模型参数进行识别，通过识别结果与模型参数真实值进行对比，识别仿真结果如图8所示。仿真计算耗时1.058 s，可见该算法具有较高的计算效率。从图8可以看出，CKF可以快速识别模型参数，识别的模型参数滞回位移、加载位移、b、支撑轴向恢复力收敛速度快，且真实值与识别值误差较小，CKF识别模型参数值滞回位移、加载位移、b、支撑轴向恢复力相对误差分别为23.2%、3.6%、2.2%、14.8%，识别参数可以有效反应模型真实性能。