发明名称 一种桥梁腐蚀监测与寿命预测方法
摘要 本发明公开了一种桥梁腐蚀监测与寿命预测方法,该方法通过三组安装在桥梁不同深度的氯离子浓度传感器以及腐蚀速率传感器所反馈的数据利用菲克第二定律等方法综合判断该桥梁当前状态距离桥梁腐蚀开裂阶段的时间,得出桥梁在没有维护情况下的剩余寿命,实现了对桥梁腐蚀状态的监测和对桥梁寿命的预测。本发明主要解决了对于已建成桥梁腐蚀状态的监测问题,并能通过传感器数据预测桥梁结构的剩余寿命,提醒维护人员及时进行桥梁维护,保障了桥梁结构的安全性和交通的正常通行。
申请公布号 CN103852414A 申请公布日期 2014.06.11
申请号 CN201410098485.2 申请日期 2014.03.17
申请人 东南大学 发明人 叶智锐;许跃如
分类号 G01N17/00(2006.01)I 主分类号 G01N17/00(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 杨晓玲
主权项 1.一种桥梁腐蚀监测与寿命预测方法,其特征在于,该方法包括以下步骤:步骤1)在桥梁混凝土桥墩中设置等距直线排列的至少三个氯离子浓度传感器,记氯离子传感器个数为N个,同时在直线两端的氯离子浓度传感器之间连线的中点处设置腐蚀速率传感器;通过氯离子浓度传感器分别采集时刻t<sub>1</sub>,t<sub>2</sub>,…t<sub>i</sub>…t<sub>S</sub>下桥梁各个点位X<sub>1</sub>、X<sub>2</sub>、…X<sub>n</sub>…X<sub>N</sub>处的氯离子浓度,其中t<sub>i</sub>为第i个时刻,i为时刻序号,i=1,2,…,S,t<sub>1</sub>为第一个数据采集时刻,t<sub>S</sub>为进行寿命预测的时刻,即当前时刻,S为进行寿命预测的时刻序号,X<sub>n</sub>为第n处点位,n=1,2,…,N,在t<sub>i</sub>时刻下所采集的各个点位的氯离子浓度值分别记为C(X<sub>1</sub>,t<sub>i</sub>),C(X<sub>2</sub>,t<sub>i</sub>),…,C(X<sub>n</sub>,t<sub>i</sub>),…C(X<sub>N</sub>,t<sub>i</sub>);同时通过腐蚀速率传感器采集时刻t<sub>1</sub>,t<sub>2</sub>,…t<sub>i</sub>…t<sub>S</sub>下的腐蚀速率,分别记为CR<sub>1</sub>,CR<sub>2</sub>,…CR<sub>i</sub>…CR<sub>S</sub>;步骤2)按照如下方法分别对各时刻的氯离子浓度C(X<sub>1</sub>,t<sub>i</sub>),C(X<sub>2</sub>,t<sub>i</sub>),…,C(X<sub>n</sub>,t<sub>i</sub>),…C(X<sub>N</sub>,t<sub>i</sub>),以及腐蚀速率CR<sub>i</sub>进行初次筛选:a)如果t<sub>i</sub>时刻深度X<sub>1</sub>、X<sub>2</sub>、…X<sub>n</sub>…X<sub>N</sub>处所对应的一组氯离子浓度值C(X<sub>1</sub>,t<sub>i</sub>),C(X<sub>2</sub>,t<sub>i</sub>),…,C(X<sub>n</sub>,t<sub>i</sub>),…C(X<sub>N</sub>,t<sub>i</sub>)有以下情形之一,则剔除该组数据:任意一个氯离子浓度值为负值,任意一个氯离子浓度值&gt;2.0M,该组氯离子浓度值的相关系数r的绝对值|r|&lt;0.75;b)如果t<sub>i</sub>时刻的腐蚀速率CR<sub>i</sub>有以下情形之一,则剔除该数据:腐蚀速率CR<sub>i</sub>为负值,腐蚀速率CR<sub>i</sub>小于上一时刻t<sub>i-1</sub>的测量值CR<sub>i-1</sub>,CR<sub>i-1</sub>&lt;0.1uA/cm<sup>2</sup>且CR<sub>i</sub>&gt;1.0uA/cm<sup>2</sup>;步骤3)首先,根据所述步骤2)初次筛选过后的数据,对于每一时刻t<sub>i</sub>的一组氯离子浓度值C(X<sub>1</sub>,t<sub>i</sub>),C(X<sub>2</sub>,t<sub>i</sub>),…,C(X<sub>n</sub>,t<sub>i</sub>),…C(X<sub>N</sub>,t<sub>i</sub>),利用菲克第二定律计算时刻t<sub>i</sub>对应的表面氯离子浓度C<sub>S(i)</sub>和扩散系数D<sub>i</sub>;然后根据下列方程,求解出每个时刻t<sub>i</sub>对应的腐蚀开始时刻T<sub>th(i)</sub>:<maths num="0001"><![CDATA[<math><mrow><msub><mi>C</mi><mi>th</mi></msub><mo>=</mo><msub><mi>C</mi><mn>0</mn></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>C</mi><mrow><mi>S</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><msub><mi>C</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>[</mo><mn>1</mn><mo>-</mo><mi>erf</mi><mrow><mo>(</mo><mfrac><msub><mi>X</mi><mi>N</mi></msub><mrow><mn>2</mn><msqrt><msub><mi>D</mi><mi>i</mi></msub><mi>Tth</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msqrt></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>其中:C<sub>th</sub>为氯离子浓度临界值,C<sub>0</sub>为初始氯离子浓度值,erf为误差函数,T<sub>th(i)</sub>为腐蚀开始时刻;步骤4)按照以下方法对每个时刻t<sub>i</sub>下的腐蚀开始时刻T<sub>th(i)</sub>和经过步骤2)初次筛选的腐蚀速率CR<sub>i</sub>进行检验:如果满足T<sub>th(i)</sub>&gt;t<sub>i</sub>,且CR<sub>i</sub>&gt;1,则剔除该组数据,否则进一步判断是否满足T<sub>th(i)</sub>≤t<sub>i</sub>,且CR<sub>i</sub>&lt;0.1,如是,则剔除该组数据,否则保留该组数据;步骤5)根据所述步骤4)中检验后的数据,按照以下方法计算每个时刻t<sub>i</sub>距离腐蚀开裂阶段开始时刻的时间,即时刻t<sub>i</sub>的腐蚀开裂预测时间T<sub>i</sub>:如果T<sub>th(i)</sub>&gt;t<sub>i</sub>,则根据公式T<sub>i</sub>=T<sub>th(i)</sub>-t<sub>i</sub>+T<sub>add</sub>计算时刻t<sub>i</sub>的腐蚀开裂预测时间T<sub>i</sub>;否则根据公式T<sub>i</sub>=T<sub>add</sub>计算时刻t<sub>i</sub>的腐蚀开裂预测时间T<sub>i</sub>,其中T<sub>add</sub>为从腐蚀开始时刻到腐蚀开裂阶段开始时刻的时间;步骤6)判断当前时刻t<sub>S</sub>对应的腐蚀开裂预测时间T<sub>S</sub>是否满足T<sub>S</sub>=T<sub>add</sub>,如是,则当前时刻t<sub>S</sub>对应的腐蚀开裂最终预测时间T=T<sub>S</sub>;否则,利用所述步骤5)中得到的T<sub>1</sub>,T<sub>2</sub>,…T<sub>i</sub>,…T<sub>S</sub>,根据下列两式拟合出线性回归方程T<sub>i</sub>=bt<sub>i</sub>+a的常数项a,和时间t<sub>i</sub>的系数b;<img file="FDA0000477785470000021.GIF" wi="688" he="139" />其中<img file="FDA0000477785470000022.GIF" wi="26" he="58" />和<img file="FDA0000477785470000023.GIF" wi="32" he="64" />分别为t<sub>1</sub>,t<sub>2</sub>,…t<sub>i</sub>,…t<sub>S</sub>和T<sub>1</sub>,T,<sub>2</sub>,…,Ti,…Ts的平均值;然后根据线性回归方程T<sub>i</sub>=bt<sub>i</sub>+a,计算出当前时刻t<sub>S</sub>对应的腐蚀开裂最终预测时间T=T<sub>S</sub>=bt<sub>S</sub>+a。
地址 211189 江苏省南京市江宁区东南大学路2号