灰色模型理论在地基沉降预测中的应用.docxVIP

灰色模型理论在地基沉降预测中的应用 1 灰色模型理论在n阶常系数测定中的应用 由于不同理论计算方法本身的局限性和工程条件的复杂性，基站沉降完全取决于理论计算，有时不准确。因此，早期测量沉降是评价后期沉降乃至最终沉降的重要意义。迄今为止,已经有许多预测沉降的方法,Asaoka法是其中应用最广泛者之一,而灰色模型理论近年来也开始得到应用。许多学者在这方面已作了研究。本文作者曾针对灰色理论方法与Asaoka法作了些比较,发现灰色GM(n,1)模型与Asaoka法都是n阶常系数微分方程,并指出当应用灰色理论时,若数据依据以下三个原则选取:①每个数据的时间间隔是相等的;②建模的数据是时间步长中的沉降增量;③数据均选自加荷结束后的观测数据,这时灰色理论方法与Asaoka法是一致的。在满足这三个条件下,本文进一步从数学上严格证明了灰色模型理论和Asaoka法两者的一阶微分方程是相等的,这对灰色理论在沉降预测中的推广应用具有重要意义。 2 证明 2.1 线性直线方程的建立 Asaoka根据Mikasa一维固结方程推导出沉降预测的一阶微分方程: s+cdsdt=d(1)s+cdsdt=d(1) 式中s为沉降,由于方程只取一阶,沉降中不宜包括次固结沉降,当数据取自加荷结束以后,瞬时沉降已完成,所以沉降s中可以包括瞬时沉降;t为时间;c,d为常系数。将式(1)写成差分格式,时间步长取单位1,得 sj=β′0+β′1sj-1(2)sj=β′0+β′1sj?1(2) 或者写成 sj-1=β0+β1sj(3)sj?1=β0+β1sj(3) 式中sj表示s(tj),β′0,β′1,β0,β1为常系数,与式(1)中的系数c,d有关: β′0=d1+c,β′1=c1+c,β0=-dc,β1=1+cc(4) 式(2),式(3)是线性直线方程,可用图解法确定方程系数,也可用最小二乘法来确定。以式(3)为例,设有N+1个等时间间隔的沉降观测数据(s0,s1, \:,sN),可以组成N组数据(sj,sj-1)(j=1,\:,N),为表达方便,记xj=sj,yj=sj-1,N组数据(xj,yj)在直角坐标系中可用直线拟合,该直线的系数根据最小二乘法来确定,则式(3)中的方程系数按下式计算: β0=Σx2Σy-ΣxΣxyΝΣx2-(Σx)2,β1=ΝΣxy-ΣxΣyΝΣx2-(Σx)2(5) 式中“Σ”表示“Ν∑j=1”。 2.2 沉降增量观测 灰色模型理论的一阶微分方程即GM(1,1)模型,表示如下: dsdt+as=b(6) 式中a,b为常系数。设有N+1个等时间间隔的沉降增量观测数据(Δs0,Δs1, \:,ΔsN),作一次累加后得各时刻的总沉降(s0,s1, \:,sN),根据灰色建模理论,方程系数按下式确定: 式中 B=[-s11\~\:\~\:-sΝ1],yΝ=[Δs1,Δs2,\:,ΔsΝ]Τ 2.3 asoaka法、gm1,1模型的初步比较见表1 如果两种方法的一阶微分方程是相等的,则式(1)和式(6)中的系数应存在以下关系: c=1a,d=ba(8) 要证明两方程是等价的,只要证明式(8)成立即可。根据式(7)得 式中δ=-ΣsjΣsj+NΣs2j,“Σ”表示“Ν∑j=1”。同样记xj=sj,yj=sj-1注意到 ΣsjΔsj=Σsj(sj-sj-1)=Σs2j-Σsjsj-1=Σx2-ΣxyΣΔsj=Σ(sj-sj-1)=Σsj-Σsj-1=Σx-Σy 所以 {ab}=1ΝΣx2+(Σx)2[-Ν(Σx2-Σxy)+Σx(Σx-Σy)-Σx(Σx2-Σxy)+Σx2(Σx-Σy)]=1ΝΣx2+(Σx)2[(Σx)2-ΣxΣy-ΝΣx2+ΝΣxyΣxΣxy-Σx2Σy](9) 根据式(4),式(5)和式(9),很快可得到 c=1β1-1=ΝΣx2-(Σx)2(Σx)2-ΣxΣy-ΝΣx2+ΝΣxy=1ad=β01-β1=ΣxΣxy-Σx2Σy(Σx)2-ΣxΣy-ΝΣx2+ΝΣxy=ba 因此式(8)成立,说明两种方法的控制方程完全一致。 Asaoka法所求的最终沉降即系数d,而GM(1,1)模型所求的最终沉降即ba,两者所预测的最终沉降将相等。 两种方法虽在本质上一致,但在应用时各有特色。Asaoka法计算相对要简单些,而灰色GM(1,1)模型计算更复杂些。在应用的灵活性上,灰色GM(1,1)模型却更有优势,可以方便地处理非等时间步长的数据,只需修改式(7)中的B矩阵: B=[Δt10?0ΔtΝ][-s11\~\:\~\:-sΝ1],Δtj=tj-tj-1,j=1,\:,Ν 当预测精度不够时,灰色理论还可建立残差模型进行修正。对于更高阶次的微分方程(n≥2),灰色理论也比Asaoka法更好处理。另外将灰色GM(1,1)模型扩充,也可方便地建立多变量