基于灰色系统的橡胶隔震支座性能老化预测*

罗佳润，马玉宏，崔 杰，沈朝勇

0 概述

橡胶隔震支座是由一层钢板一层橡胶经高温硫化而成的 (日本建筑学会，2006)，其外表面通常设置一定厚度的橡胶保护层来保护内部结构。影响橡胶隔震支座耐久性的因素非常多，主要包括:橡胶材料的老化、橡胶隔震支座的徐变、疲劳 (GB/T 20688.1－2007;全国橡胶与橡胶制品标准化技术委员会橡胶杂品分技术委员会，广州大学工程抗震研究中心，2009)。支座性能的改变是由长期恒定作用、地震荷载作用、行车荷载作用及使用环境中的各种光、热、氧、化学介质等因素引起的 (周福霖，1997;苏经宇等，2012)。橡胶隔震支座的力学性能直接关系到隔震结构在遭遇地震时的响应及抗震性能，因此对橡胶隔震支座性能随时间老化变化规律开展研究是非常必要的。

1 橡胶隔震支座老化性能研究现状

目前，对于橡胶隔震支座老化性能的研究主要是依靠热老化试验 (裴若娟，唐家祥，1993;许斌，唐家祥，1995)和试样实测 (李杨海，2009)两种方法。但对橡胶隔震支座性能随时间老化的实测，非常耗时耗力，试验成本非常昂贵，获得的数据却很有限。另外，对于橡胶隔震支座的老化性能开展的试验研究，通常只是采集初始状态 (老化前)和终止状态(老化后)时的力学性能指标且均采用人工加速老化的试验方法。对试验数据的处理，通常是采用回归分析的方法，在得到回归公式后进行预测，但是对于隔震支座老化性能这种“贫数据”的情况，回归分析的方法在老化性能预测方面具有明显的局限性。

2 灰色系统预测方法简介

邓聚龙 (1993)于20世纪80年代提出了灰色理论，该理论已经广泛应用于农业、经济、医疗、工业、管理、气象、教育、历史等几十个领域，取得较好的经济效益，并取得了许多重大成果 (刘思峰，党耀国，2010)。灰色理论以“部分信息已知、部分信息未知”的“小样本”、“贫信息”不确定性系统为研究对象，主要通过对“部分”已知信息的生成、开发，提取有价值的信息，实现对系统运行行为、演化规律的正确描述和有效监控 (党耀国，刘思峰，2009)。橡胶隔震支座老化性能试验数据符合这种“小样本”、“贫信息”的特点，预测其老化性能可以采用此方法。

3 基于灰色理论的橡胶隔震支座性能老化预测

本文采用的研究数据来源于Morita等 (2009)对直径500 mm的橡胶隔震支座 (表1)跟踪20年研究所得到的性能指标 (表2)。

表1 橡胶隔震支座参数Tab.1 Parameters of laminated rubber bearing

表2 橡胶隔震支座刚度Tab.2 Stiffness of laminated rubber bearing

为对橡胶隔震支座老化性能 (竖向和水平刚度)进行灰色预测，建立传统GM(1，1)模型:采用初始、5年后、10年后、15年后的竖向及水平刚度作为原始数据序列，用20年后的隔震支座竖向及水平刚度值来检验预测模型，预测步骤及结果如下:

初始值:

(1)对 X(0)作一次累加生成 (1－AGO)，得到:

(2)对X(1)作紧邻值生成。令

得到:Z(1){=z(1)(2)，z(1)(3)，z(1)(4}).

(3)对参数a^=(a，b)T进行最小二乘估计:得到:

a^=(BTB)－1BTY.

(5)还原数列序列x^(k+1)=x^(1)(k)－x^(0)(k－1)，依次得到:x^(0)(1)，x^(0)(2)，x^(0)(3)，x^(0)(4)。

(6)对结果进行误差检验:用灰色预测方法得到的老化性能 (水平与竖向刚度)预测值与实测值的对比见表3。

表3 误差检验Tab.3 Deviation checking

以上采用前15年的数据建立GM(1，1)模型，以建立的GM(1，1)模型对隔震支座使用了20年后的竖向刚度及水平刚度进行预测，并将实测结果与预测结果作了比较。由表3可知:对于隔震支座在使用了20年后的竖向刚度预测存在5.731%的误差，而对于水平刚度的预测则相对误差较小，仅相差2.815%。从比较结果来看，采用灰色系统理论对隔震支座的力学性能预测具有较高的精度。

4 灰色系统不同的GM(1，1)模型对预测结果的对比研究

在作灰色系统预测时，建立合理的GM(1，1)灰色模型，确保 GM(1，1)在原点的精度非常重要。为了提高GM(1，1)在原点的精度通常可以采用全数据GM(1，1)模型、新息 GM(1，1)模型及新陈代谢GM(1，1)模型。

新息GM(1，1)模型可将不断进入系统的扰动因素考虑进去，GM(1，1)模型要将每一个新得到的数据送入X(0)中，重新建立GM(1，1)模型，即原始数列为X(0){=x(0)(1)，x(0)(2)，…，x(0)(n}) ，当获得第n+1个数据后，重新构造数列:

随时间的推移，老的信息将越来越不适应新情况，新陈代谢GM(1，1)模型将最新获得的信息考虑进来，同时舍弃最老的信息，即原始数列为X(0){= x(0)(1)，x(0)(2)，…，x(0)(n}) ，当获得第n+1个数据后，舍弃最老信息，加入新获得的信息，重新构造数列:

就本文算例来说，以使用了20年的隔震支座力学性能参数作为新获得的信息，分别建立新息GM(1，1)模型和新陈代谢GM(1，1)模型。

为了对使用了30年后的隔震支座力学性能情况作预测，分别采用新息GM(1，1)模型 (即在上述GM(1，1)模型中加入使用了20年的隔震支座力学性能参数)和新陈代谢模型 (即初始数据舍去，并新加使用了20年的隔震支座力学性能参数)，并将得到的结果作相应的对比，结果见表4，表5。可见，相对来说由新陈代谢模型求得的模拟值误差小于新息模型。

为预测橡胶隔震支座使用25年、30年、35年、40年的力学性能，作发展带预测。将上述建立的3种GM(1，1)模型所对应的估计参数序列a^(1)=［ai，bi］T(1，2，3)，列出如下:

表4 新息GM(1，1)模型及新陈代谢GM(1，1)模型对竖向刚度的模拟值对比Tab.4 Comparison of simulated value of vertical stiffness by using new information GM(1，1)model and metabolism GM(1，1)model

表5 新息GM(1，1)模型及新陈代谢GM(1，1)模型对水平刚度的模拟值对比Tab.5 Comparison of simulated value of horizontal stiffness by using new information GM(1，1)model and metabolism GM(1，1)model

因为

所以发展带上界序列响应式为

因X(0)3的起点为x(0)(2)，故发展带下界序列时间响应式为

由此可得竖向刚度最低预测值为

由竖向刚度最高值和竖向刚度最低值得到基本预测值:

对于水平刚度数据，按照相同的方法，得到水平刚度最高预测值为

水平刚度最低预测值为

由水平刚度最高值和水平刚度最低值得到基本预测值:

为了直观清晰表示，将3种模型的模拟、预测结果绘图如图1、2所示。

图1 竖向刚度模拟及预测值对比Fig.1 Comparison of simulated and predicted value of vertical stiffness

图2 水平刚度模拟及预测值对比Fig.2 Comparison of simulated and predicted value of horizontal stiffness

从图1可以看出，由于在第20年的竖向刚度相比第15年的降低，出现了一个拐点。采用前15年的数据建立的GM(1，1)模型，缺乏对于新数据的利用，造成在第20年的模拟值差别较大，而采用新息GM(1，1)模型及新陈代谢GM(1，1)模型的模拟、预测结果与实测的数据具有较高的拟合度，尤其在原点 (即第20年的竖向刚度)具有很高的模拟精度，新息GM(1，1)模型所得结果仅与实测值相差1.663%，新陈代谢GM(1，1)模型结果则相差1.562%。从图2来看，采用新息GM(1，1)M模型在原点处的模拟值与实测值仅相差0.815 5%，采用新陈代谢GM(1，1)模型在原点处的模拟值与实测值误差更是低至0.384 2%，模拟精度非常高。

与此同时，本文还结合上述几种灰色预测模型，对隔震支座的竖向刚度及水平刚度进行了使用40年后的发展带预测，得出的竖向刚度基本预测值为1 222.96 kN/mm，比初始值增大了15.592%，水平刚度基本预测值为978.21 N/mm，比初始值增大了15.355%。这个预测结果与日本试验得到的结论经过30～60年老化后，水平剪切刚度提高约20%(周福霖，1997)相一致。

5 结论

本文采用3种GM(1，1)模型对橡胶隔震支座的性能随时间的老化作了相应的模拟、预测，由结果可知:灰色系统预测理论用于隔震支座老化性能的模拟预测具有较高的精度。对于一些重要的隔震建筑、隔震桥梁的监测，可以抽样检测隔震支座的力学性能，并在这基础上采用灰色系统预测理论对隔震支座力学性能随时间的老化进行预测，从而为预测隔震结构的抗震性能提供必要的依据。

党耀国，刘思峰.2009.灰色预测与决策模型研究［M］．北京:科学出版社.

邓聚龙.1993.灰色控制系统(第二版)［M］．武昌:华中理工大学出版社.

李杨海.2009.公路桥梁支座实用手册［M］．北京:人民交通出版社，50－52.

刘思峰，党耀国.2010.灰色系统理论及其应用背景［M］．北京:科学出版社.

日本建筑学会.2006.隔震结构设计［M］．刘文光译.冯德民校.北京:地震出版社.

裴若娟，唐家祥.1993.叠层橡胶隔震器的耐久性［J］．工程抗震，(3):36－39.

全国橡胶与橡胶制品标准化技术委员会橡胶杂品分技术委员会，广州大学工程抗震研究中心.2009.橡胶支座国家标准理解与实施［M］．北京:中国标准出版社，71.

苏经宇，曾德民，田杰.2012.隔震建筑概论［M］．北京:冶金工业出版社.

许斌，唐家祥.1995.基础隔震叠层橡胶支座耐久性试验研究［J］．工程抗震，(4):41－44.

周福霖.1997.工程结构减震控制［M］．北京:科学出版社.

GB/T 20688.1－2007，橡胶支座第1部分:隔震橡胶支座试验方法［S］．

Morita K.，Yamagami S.，Takayama M.Long-term performance test of laminated rubber bearing for seismic isolation system［J］//rd International Conference On Advances in Experimental structural Engineering，CD －ROM(2009.10).