发明名称 一种混凝土坝损伤场演化状态的预警方法
摘要 本发明公开了一种混凝土坝损伤场演化状态的预警方法,按如下步骤进行:(1)构建基于熵的损伤场指示变量,系统的熵定义为:H=F(D,D<sub>i</sub>,L,T);(2)计算损伤场指示变量,坝体损伤场的信息熵为:<img file="DDA00002815004600011.GIF" wi="772" he="139" />(3)依据步(2)计算得到的损伤场指示变量时间序列,基于尖点突变模型,求解熵变指标的突变特征值;(4)依据步骤(3)确定的熵变指标突变特征值作为控制准则来判断损伤场的演化情况,突变特征值Δ=4u<sup>3</sup>+27v<sup>2</sup>,Δ>0时,损伤场总是处于稳定状态;Δ=0时,损伤场处于临界状态,外界微小的干扰有可能造成破坏;Δ<0时,损伤场发生突变,远离平衡态,发出预警。本发明对多测点测值变化分析,实现准确预警,指导工程施工,减少安全事故。
申请公布号 CN103114554A 申请公布日期 2013.05.22
申请号 CN201310044199.3 申请日期 2013.02.04
申请人 河海大学;南京河海科技有限公司 发明人 苏怀智;顾冲时;王校利;杨孟;李金友;胡江;欧斌;秦鹏;傅兆庆
分类号 E02B3/04(2006.01)I 主分类号 E02B3/04(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 缪友菊
主权项 1.一种混凝土坝损伤场演化状态的预警方法,其特征在于,按如下步骤进行:(1)构建基于熵的损伤场指示变量选取混凝土坝的应力应变计测点,在测点上设置应变计,通过应力应变计采集应变测值,同时记录测点的位置坐标和监测的时间,依据测点的位置坐标和应变测值,计算出坝体损伤场的各个测点损伤度和总损伤度,构建基于熵的损伤场指示变量,系统的熵定义为:H=F(D,D<sub>i</sub>,L,T),式中:D为坝体损伤场的总损伤度,D<sub>i</sub>为系统各个测点的损伤度,L为测点的位置坐标,T为监测的时间;(2)计算损伤场指示变量根据应力应变计测点的监测资料采用指数函数模型计算各测点的损伤度D<sub>i</sub>,<img file="FDA00002815004300011.GIF" wi="487" he="118" />D<sub>0</sub>为初始损伤值,D<sub>f</sub>为损伤阈值,即峰值应变ε<sub>f</sub>对应的损伤值,构建出坝体局部区域的损伤场,定义坝体损伤场的信息熵为:<img file="FDA00002815004300012.GIF" wi="774" he="138" />式中,D<sub>i</sub>为各测点的损伤度,i=1,2,……n;D为各测点的损伤度D<sub>i</sub>的加权迭加;λ<sub>i</sub>≥0(i=1,2,...n)表示第i个测点的损伤值在系统总损伤中所占的份额,即λ<sub>i</sub>=D<sub>i</sub>/D(i=1,2,...,n);(3)依据步(2)计算得到的损伤场指示变量时间序列,基于尖点突变模型,求解熵变指标的突变特征值设系统的状态由m个状态变量x<sub>1</sub>,x<sub>2</sub>,...,x<sub>m</sub>来表示,外界控制条件由r个控制变量u<sub>1</sub>,u<sub>2</sub>,...,u<sub>r</sub>表征,则系统的动力学性质由下述势函数来描述:P(x)=f(x<sub>1</sub>,x<sub>2</sub>,...,x<sub>m</sub>,u<sub>1</sub>,u<sub>2</sub>,...,u<sub>r</sub>)=f(x,u),u<sub>1</sub>,u<sub>2</sub>,...,u<sub>r</sub>构成r维控制空间,而x<sub>1</sub>,x<sub>2</sub>,...,x<sub>m</sub>,u<sub>1</sub>,u<sub>2</sub>,...,u<sub>r</sub>构成m+r维相空间;根据截断代数可对P(x)的Taylor级数展开式f<sub>u</sub>(x)进行k次截断,记为<img file="FDA00002815004300013.GIF" wi="205" he="118" />忽略高次项,得到<maths num="0001"><![CDATA[<math><mrow><mi>P</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mover><mo>=</mo><mo>&CenterDot;</mo></mover><msup><mover><mrow><msub><mi>f</mi><mi>u</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mi>k</mi></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>r</mi><mo>=</mo><mn>0</mn></mrow><mi>k</mi></munderover><mfrac><mn>1</mn><mrow><mi>r</mi><mo>!</mo></mrow></mfrac><msup><mi>D</mi><mi>r</mi></msup><msub><mi>f</mi><mi>u</mi></msub><msub><mo>|</mo><mn>0</mn></msub><msup><mi>x</mi><mi>r</mi></msup><mo>,</mo></mrow></math>]]></maths>式中:f<sub>u</sub>(x)=f(x,u);由分类定理知,若P(x)为无限次连续可导的光滑函数,那么,DP(x)=0,R<sup>m+r</sup>中的微分流形M为:<maths num="0002"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>P</mi></mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>1</mn></msub></mfrac><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mfrac><mrow><mo>&PartialD;</mo><mi>P</mi></mrow><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mi>m</mi></msub></mfrac><mo>=</mo><mn>0</mn><mo>;</mo></mrow></math>]]></maths>奇点集S是M的一个子集,对应于将会发生突变体系的平衡状态,S在控制空间R<sup>r</sup>中的投影为分叉集B,在M的任意点上,B关于P(x)所受到的摄动是局部稳定的,在非临界点和非退化临界点不是突变的形式,在退化临界点附近为突变的形式;正则尖点突变的标准势函数为:<img file="FDA00002815004300023.GIF" wi="464" he="108" />对标准势函数求导,并令<img file="FDA00002815004300024.GIF" wi="167" he="111" />即:Z<sup>3</sup>+uZ+v=0;因此,即奇点集S的方程为:<img file="FDA00002815004300025.GIF" wi="383" he="118" />分叉集B的方程为:27v<sup>2</sup>+4u<sup>3</sup>=0;根据突变理论,选取如下尖点突变平衡曲面方程为:Z<sup>3</sup>+uZ+v=0,式中的控制参数u、v随着时间变化,为时间变量t的连续可导函数,即u=u(t),v=v(t);(4)依据步骤(3)确定的熵变指标突变特征值作为控制准则来判断损伤场的演化情况建立一个连续函数H<sub>j</sub>=f(D<sub>j</sub>)来表示系统熵值H<sub>j</sub>(m)随损伤值D<sub>j</sub>的变化规律,将此函数进行泰勒级数展开至5次项:<maths num="0003"><![CDATA[<math><mrow><msub><mi>H</mi><mi>j</mi></msub><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>D</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>4</mn></munderover><msub><mi>a</mi><mi>i</mi></msub><msubsup><mi>D</mi><mi>j</mi><mi>i</mi></msubsup></mrow></math>]]></maths>其中,<maths num="0004"><![CDATA[<math><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mi>i</mi></msup><mi>f</mi></mrow><msubsup><mrow><mo>&PartialD;</mo><mi>D</mi></mrow><mi>j</mi><mi>i</mi></msubsup></mfrac><mo>|</mo><mi>D</mi><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>令<img file="FDA00002815004300028.GIF" wi="358" he="132" />则上式的尖点突变的标准势函数形式:H<sub>j</sub>′=x<sup>4</sup>+ux<sup>2</sup>+vx+ω式中:<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>H</mi><mi>j</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><msub><mi>H</mi><mi>j</mi></msub><msub><mi>a</mi><mn>4</mn></msub></mfrac><mo>;</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>u</mi><mo>=</mo><mfrac><msub><mi>a</mi><mn>2</mn></msub><msub><mi>a</mi><mn>4</mn></msub></mfrac><mo>-</mo><mfrac><msup><msub><mrow><mn>3</mn><mi>a</mi></mrow><mn>3</mn></msub><mn>2</mn></msup><msup><msub><mrow><mn>8</mn><mi>a</mi></mrow><mn>4</mn></msub><mn>2</mn></msup></mfrac><mo>;</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mi>v</mi><mo>=</mo><mfrac><msub><mi>a</mi><mn>1</mn></msub><msub><mi>a</mi><mn>4</mn></msub></mfrac><mo>-</mo><mfrac><mrow><msub><mi>a</mi><mn>2</mn></msub><msub><mi>a</mi><mn>3</mn></msub></mrow><msup><msub><mrow><mn>2</mn><mi>a</mi></mrow><mn>4</mn></msub><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><msup><msub><mi>a</mi><mn>3</mn></msub><mn>3</mn></msup><msup><msub><mrow><mn>8</mn><mi>a</mi></mrow><mn>4</mn></msub><mn>3</mn></msup></mfrac><mo>;</mo></mrow></math>]]></maths>ω为常数;突变特征值Δ=4u<sup>3</sup>+27v<sup>2</sup>,Δ>0时,损伤场总是处于稳定状态;Δ=0时,损伤场处于临界状态,外界微小的干扰有可能造成破坏;Δ<0时,损伤场发生突变,远离平衡态,发出预警。
地址 210009 江苏省南京市鼓楼区西康路1号