东海大桥基于温度序列挠度研究.doc

东海大桥基于温度序列挠度研究 　　摘 要 东海大桥工程是上海国际航运中心洋山深水港区一期工程的重要配套工程，为洋山深水港区集装箱陆路集疏运和供水、供电、通讯等需求提供服务。 挠度是指建筑体在一定条件下发生弯曲形变时横截面形心沿与轴线垂直方向的线位移。本文主要对挠度与温度之间的关系进行讨论和探究。力求最终所得模型对数据有较好的拟合，并用于对桥梁挠度的预测分析。 关键词 桥梁 挠度 温度 时间序列分析 主成分分析 中图分类号：U448 文献标识码：A 在桥梁建成后，由于受到气候、环境因素的影响，会被逐渐腐蚀和老化，长期的荷载作用，使其强度和刚度降低。不仅会影响行车安全，也会使使用寿命缩短。对结构布局十分复杂的大型桥梁上若仍然沿用传统方法的进行外观检查、养护以及局部检测，已远远难以完备的反映其健康状况，尤其是桥梁的安全问题难以保障。建立一种用以提供全面检测，实时观测大桥结构的承载能力，保证桥体的安全运营。 综合已有研究资料，为建立挠度和温度关系数量模型，研究研究思路如下：首先对温度序列与挠度序列进行获取和整理，对于其中的异常值选定规则予以剔除，其次尝试建立若干可行模型，通过对拟合优度，多重共线性分析和残差序列分析等方法对所建模型进行合理性评价最后，对最终所提出模型进行优缺点评析并提出应用办法。（由于版面限制，输出图表基本略去，仅保留最终拟合结果及少量必要图表） 一、数据集概述及初步分析 在日常观测中，温度随时间的推移而周期性变化，而温度又为挠度的主要影响因素之一，故尝试将数据集DATA视为时间面板数据进行分析并建立适当模型对挠度（GP）进行估计并对数据集DATA中的缺失值进行填补。 （一）数据的性质检验及筛选。 根据时间序列分析的要求，不论自变量或因变量，均不得为白噪声序列。经分析，挠度（GP）与18个温度序列的自相关系数与偏自相关系数显著不为零，因而不是白噪声序列。又通过ADF检验方法，在1%、5%、10%的置信度下，GP（挠度序列，下文中挠度序列可用GP代指）及温度序列均显著拒绝存在单位根的原假设，即GP平稳。最后由6 法删去由仪器失常及其他因素造成的异常值。 （二）缺失数据的填补。 挠度的具体表现为桥梁在温湿度及其他自然环境的共同作用下，由于不同结构体导热性能以及膨胀系数不同而导致的桥体中心点具有连续性的位移变化（隆起或下沉）。因而在此尝试用自回归模型对序列GP进行拟合。 通过自回归图，发现GP在短期内存在的偏自回归并一阶截尾，自回归系数拖尾，即GP有明显的一阶自回归趋势。尝试用模型1对其进行拟合。 其中C1为一维系数向量， t为白噪声项。 其自回归图亦表明残差并未完全呈现白噪声的特征，换用表达式2对GP进行拟合。 其中，GP的一阶差分系数显著，且差分项D（GP）的增加使得R2大幅提高并达标准因而选用模型2对空白数据进行填补。 二、线性时间序列模型的建立 线性参数模型的建立： 协方差分析显示十八个温度序列之间存在着强相关；而任意温度序列与挠度的相关系数均小于0.2，可以认为无线性相关。于是放弃变量线性拟合，寻求温度与挠度之间的关系。尝试下列模型： 其中TEM=（CT1，CT2，…，CT6，ST1，ST2，…ST12）’为18??980的温度矩阵，C1为1??系数向量，C2为1??的系数向量，C0表示残差项，LNGP=ln（-GP）。 拟合后，前三个方程拟合程度极弱，不能有效估计GP的数值和波动情况均予以舍去，其相关输出与分析由于篇幅限制均不列示。仅有LNGP=C?TEM在无自回归项时拟合优度既已达到R2=0.49，推测此模型在经过修正和扩大样本量之后可能会较好的预测LNGP。在进行修正之前舍弃 t 检验参数不显著的变量。 步骤如下，记为操作1： 1）删去模型中最不显著序列（相应 p 值最大的序列）； 2）舍弃此序列对模型重新拟合； 3）重复步骤1）、2），直到所有变量系数都显著； 4）对模型进行最终拟合。 所得残差的自回归系数拖尾，偏自回归系数一阶截尾。表明残差序列存在显著地一阶自回归。对模型加入一阶自回归并重复操作1中1）～4）步， 得R2=0.615148，较上一模型虽有提高但仍不完善。对残差进行同样操作，发现无短期相关。回忆本文先前所提及的挠度受到周期性序列温度的影响。绘制100阶滞后自相关图得残差序列在滞后96阶时对当前阶影响显著，向模型中加入AR（96），并重复操作步骤1）～4）。重新拟合后可以看出，通过自回归项的添加，自变量的数目的显著减少和R2的上升证明以该方向探究温度序列与挠度序列之间的参数关系的方法是可取的。再观察相应的残差序列100阶滞后的自回归图