灰色系统模型在混凝土坝应力分析中的应用.pdfVIP

灰色系统模型在混凝土坝应力分析中的应用 在混凝土坝应力观测资料分析工作中，常采用统计模型和混合数学模 型来分析大坝应力的发展规律以及影响因素。统计模型和混合模型要 求具有典型分布规律的长序列数据，而在施工期和蓄水初期，由于观 测数据很少，采用上述两种模型很难取得较好的效果。灰色系统模型 的特点是可以从较少的离散数据中寻找内在的规律。 灰色模型简称 GM 模型，GM 模型的建模包括因素分析、等时距处理、 数据标准化、数据生成、建立状态模型、模型效果评价等方法和过程。 本文只探讨灰色理论在混凝土坝应力观测资料分析中的应用。 1、因子选择 在研究大坝的应力规律时，应力数列是参考数列，影响应力的行 为因子则为比较数列，它们均为时间数列。混凝土坝应力主要与水压 力、温度、自重、湿涨以及时效等因素有关系。在蓄水初期的观测资 料分析中，自重应力和湿涨应力以常数表示。 (1)库水位因子。选择水深的三次多项式：H、H2、H3。 (2)时效因子。选择线性函数和对数函数的叠加函数： (3)温度因子。选择当日气温 T、前 30 天平均气温 T30、前 60 天 平均气温 T60、前 1～6 个月的某日气温 Tqi(i=1～6)、前 1～6 个月 的月平均气温Tqi(i=1～6)。 于是，参考数列：应力 x1(0)={x1(k)(0)|k=1,2,3m}； 比较数列：库水位、时效、温度 xi(0)={xi(k)(0)|k=1,2,3,...m} (i=2,3,...n)，m 为样本数，n 为变量数。 2、模型建立 在蓄水初期由于水位较低，水压力对应力的影响并不明显，因此 考虑对水压应力分量采用有限元计算的确定性模型。温度和时效应力 分量采用 GM 模型，将以上各种因素的模型求和就得出表征坝体实测 应力变化规律的灰色混合模型。 2.1 水压应力分量确定性模型 为了寻求水位与应力之间的确定性函数关系，可采用有限元法进 行计算，以便提高数学模型值的精度。其一般方程见下式。 式中：alpha;0、alpha;i 为回归系数。 2.2 温度和时效应力分量 GM 模型 经过等时距处理、数据标准化和数据生成可得到生成数列 xit(1)={xit(k)(1)|k=1,2,...,m}　 (i=1,2,...,n) 通过编制 GM 程序可以很方便的建立动态和静态模型，动态模型 x1(1)的时间（离散）近似关系式为： 静态模型 x1(1)的时间（离散）近似关系式为： 对模型值进行数据还原可用于模型效果的评价，通过后验差比值C 和小误差概率 P 可综合评价模型的精度。在实际建模中，模型的精度 与选入的因子有很大关系，可以通过剔除关联度小的因子来提高模型 精度。 2.3 灰色混合模型的建立 将水压应力分量确定性模型与温度时效应力分量GM模型求和就得 出灰色混合模型，其一般方程式如下： 3、工程实例 桃林口水库位于滦河主要支流青龙河上，工程于 1992 年 10 月动 工，1998 年底竣工。水库大坝为碾压混凝土重力坝，最大坝高 74.5 m，坝顶长度 500 m，坝顶宽 7m，坝底最大宽度 63.76 m。5 号坝段为 应力监测的重点坝段，在距建基面 6 m 的 80.00 m 高程选择应力观测 截面一个，在观测截面的中心线上布设 5 组 5 向应变计组测点，每个 测点布设无应力计一个，其中靠上游的第 1组测点位于常态混凝土中， 其余 4 个测点位于碾压混凝土中。本例以第 1、3、5 测点垂直向应力 为分析对象，选取 1998 年 3 月至 1998 年 11 月的观测时段进行建模。 3.1 水压应力分量 选取 35 级不同上游水深，通过有限元计算相应的应力，然后采用 逐步回归的方法进行拟合，计算结果见表 1。 表 1 中各参数正负号有所不同，说明水压对不同部位测点应力的 影响是不同的。各方程的复相关系数都很接近 1，表明全部因子与应 力的相关程度十分密切，剩余标准差都很接近 0，表明模型回归值与 计算值之间的离散程度微小，因此，采用表 1 的模型回归拟合值作为 灰混模型的水压分量是比较理想的。 3.2 温度时效应力分量 对各因子进行一次累加生成，使用 GM 程序计算动态和静态模型参 数，通过对比较数列的舍取使模型精度达到最高，计算结果见表 2。 表 2 中模型精度均为Ⅰ级，