发明名称 基于监测数据与温度应力分析的大跨钢箱梁桥的总应力计算方法及安全预警方法
摘要 本发明公开了一种基于监测数据与温度应力分析的大跨钢箱梁桥的总应力计算方法及安全预警方法,其中总应力计算方法包括步骤一、通过健康监测系统采集到的桥梁的应变与温度;步骤二、从实测应变中分离出温度应变;步骤三、计算钢箱梁上截面上的均匀温度和梯度温度;步骤四、轴向约束应力计算;步骤五、弯曲约束应力计算;步骤六、温度自应力计算;步骤七、从有限元模型中获取自重应力;步骤八、大跨桥梁在运营期间总应力计算。相比传统预警方法而言,本发明预警指标明确、力学模型清晰,并考虑了温度荷载对大跨桥梁结构的影响,适合于工程界推广。
申请公布号 CN106092402B 申请公布日期 2017.04.26
申请号 CN201610373912.2 申请日期 2016.05.31
申请人 东南大学 发明人 张建;夏琪;刘森林
分类号 G01L5/00(2006.01)I 主分类号 G01L5/00(2006.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种基于监测数据与温度应力分析的大跨钢箱梁桥的结构总应力计算方法,其特征在于,包括以下步骤:步骤一、通过健康监测系统采集到的桥梁的应变与温度:顶板应变ε<sub>U1</sub>,ε<sub>U2</sub>,…ε<sub>Ui</sub>…ε<sub>Un</sub>,底板应变ε<sub>L1</sub>,ε<sub>L2</sub>,…ε<sub>Li</sub>…ε<sub>Ln</sub>,顶板温度T<sub>U1</sub>,T<sub>U2</sub>,…T<sub>Ui</sub>…T<sub>Un</sub>,底板温度T<sub>L1</sub>,T<sub>L2</sub>,…T<sub>Li</sub>…T<sub>Ln</sub>;式中,ε<sub>Ui</sub>和ε<sub>Li</sub>分别表示第i截面顶板和底板的监测应变数据;T<sub>Ui</sub>和T<sub>Li</sub>表示第i个截面顶板和底板的监测温度数据;步骤二、从实测应变中分离出温度应变:顶板温度应变ε<sub>TU1</sub>,ε<sub>TU2</sub>,…ε<sub>TUi</sub>…ε<sub>TUn</sub>,底板温度应变ε<sub>TL1</sub>,ε<sub>TL2</sub>,…ε<sub>TLi</sub>…ε<sub>TLn</sub>,以及顶板车载应变ε<sub>VU1</sub>,ε<sub>VU2</sub>,…ε<sub>VUi</sub>…ε<sub>VUn</sub>,底板车载应变ε<sub>VL1</sub>,ε<sub>VL2</sub>,…ε<sub>VLi</sub>…ε<sub>VLn</sub>;步骤三、计算钢箱梁上截面上的均匀温度和梯度温度:截面均匀温度T<sub>Avgi</sub>=T<sub>Li</sub>其中,T<sub>Avgi</sub>表示第i截面的均匀温度;截面梯度温度T<sub>yi</sub>=T<sub>Ui</sub>‑T<sub>Li</sub>其中,T<sub>yi</sub>表示第i截面的梯度温度;步骤四、轴向约束应力计算:均匀温度产生的无约束变形δ<sub>T</sub>=αLT<sub>Avgi</sub>其中,δ<sub>T</sub>是钢箱梁在均匀温度下产生的轴向变形;α表示钢箱梁材料的膨胀系数;L表示桥梁结构有效跨径;由实测顶板和底板温度应变计算有轴向约束情况下的变形<maths num="0001"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>l</mi></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>N</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mfrac><mrow><msub><mi>&epsiv;</mi><mrow><mi>T</mi><mi>U</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>&epsiv;</mi><mrow><mi>T</mi><mi>L</mi><mi>i</mi></mrow></msub></mrow><mn>2</mn></mfrac><mi>&Delta;</mi><mi>l</mi></mrow>]]></math><img file="FDA0001196562870000011.GIF" wi="455" he="83" /></maths>其中,δ<sub>l</sub>是钢箱梁在有轴向约束下产生的实际变形;n表示传感器布置了n截面;Δl表示按传感器数目将钢箱梁跨径划分为长度,Δl=L/n;非线性温度梯度产生的无约束轴向变形σ<sub>RT</sub>(y)=E·α·T<sub>yi</sub><maths num="0002"><math><![CDATA[<mrow><msub><mi>N</mi><mrow><mi>T</mi><mi>y</mi></mrow></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><msub><mi>h</mi><mn>0</mn></msub><mo>/</mo><mn>2</mn></mrow><mrow><msub><mi>h</mi><mn>0</mn></msub><mo>/</mo><mn>2</mn></mrow></msubsup><msub><mi>&sigma;</mi><mrow><mi>R</mi><mi>T</mi></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>b</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mi>d</mi><mi>y</mi></mrow>]]></math><img file="FDA0001196562870000012.GIF" wi="622" he="94" /></maths><maths num="0003"><math><![CDATA[<mrow><msubsup><mi>&epsiv;</mi><mi>T</mi><mi>N</mi></msubsup><mo>=</mo><mfrac><msub><mi>N</mi><mi>T</mi></msub><mrow><mi>E</mi><mi>A</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001196562870000013.GIF" wi="175" he="86" /></maths>其中,σ<sub>RT</sub>(y)表示非线性梯度温度完全转化为的温度应力;N<sub>Ty</sub>表示温度应力等效的轴力;<img file="FDA0001196562870000014.GIF" wi="60" he="63" />表示等效轴力下产生的应变;h<sub>0</sub>表示截面上表面距离截面中心轴的距离;b(y)表示截面宽度,且随截面高度变化;A表示截面面积;E表示结构材料的弹性模量;最后非线性温度梯度产生的无约束轴向变形可以计算为;<maths num="0004"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mrow><mi>T</mi><mi>y</mi></mrow></msub><mo>=</mo><msubsup><mi>L&epsiv;</mi><mi>T</mi><mi>N</mi></msubsup></mrow>]]></math><img file="FDA0001196562870000021.GIF" wi="230" he="71" /></maths>轴向约束应力计算<maths num="0005"><math><![CDATA[<mrow><msub><msup><mi>&sigma;</mi><mo>&prime;</mo></msup><mi>N</mi></msub><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><msub><mi>&delta;</mi><mi>T</mi></msub><mo>+</mo><msub><mi>&delta;</mi><mrow><mi>T</mi><mi>y</mi></mrow></msub><mo>-</mo><msub><mi>&delta;</mi><mi>l</mi></msub><mo>)</mo></mrow></mrow><mi>L</mi></mfrac></mrow>]]></math><img file="FDA0001196562870000022.GIF" wi="385" he="94" /></maths>其中,σ′<sub>N</sub>表示钢箱梁的轴向约束应力;步骤五、弯曲约束应力计算非线性温度梯度产生的无约束弯曲应变及梁端转角<maths num="0006"><math><![CDATA[<mrow><msub><mi>M</mi><mi>T</mi></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><msub><mi>h</mi><mn>0</mn></msub><mo>/</mo><mn>2</mn></mrow><mrow><msub><mi>h</mi><mn>0</mn></msub><mo>/</mo><mn>2</mn></mrow></msubsup><msub><mi>&sigma;</mi><mrow><mi>R</mi><mi>T</mi></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>b</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>y</mi><mi>d</mi><mi>y</mi></mrow>]]></math><img file="FDA0001196562870000023.GIF" wi="668" he="95" /></maths><maths num="0007"><math><![CDATA[<mrow><msubsup><mi>&epsiv;</mi><mi>T</mi><mi>M</mi></msubsup><mo>=</mo><mfrac><mrow><msub><mi>M</mi><mi>T</mi></msub><mo>&CenterDot;</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><mrow><mi>E</mi><mi>I</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001196562870000024.GIF" wi="230" he="87" /></maths>θ<sub>T</sub>=M<sub>T</sub>L/2EI其中,M<sub>T</sub>表示非线性温度梯度产生的温度应力等效为弯矩;<img file="FDA0001196562870000025.GIF" wi="61" he="63" />表示等效弯矩产生的应变;θ<sub>T</sub>表示等效弯矩在无约束状态下产生的梁端转角;EI表示截面抗弯刚度;由实测顶板和底板温度应变计算梁端的实际转角<maths num="0008"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>U</mi></msub><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>&epsiv;</mi><mrow><mi>T</mi><mi>U</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>&epsiv;</mi><mrow><mi>T</mi><mi>L</mi><mi>i</mi></mrow></msub><mo>)</mo><mi>L</mi></mrow><mrow><mn>2</mn><mi>h</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001196562870000026.GIF" wi="349" he="95" /></maths>其中,h表示钢箱梁截面高度;那么,实际被约束的转角以及与边界约束弯矩的关系表示为θ′<sub>R</sub>=θ<sub>T</sub>‑θ<sub>U</sub><maths num="0009"><math><![CDATA[<mrow><msubsup><mi>&theta;</mi><mi>R</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><msub><mi>M</mi><mi>S</mi></msub><mi>L</mi></mrow><mrow><mn>2</mn><mi>E</mi><mi>I</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001196562870000027.GIF" wi="190" he="86" /></maths>其中,θ′<sub>R</sub>表示被约束住的梁端转角;M<sub>s</sub>表示边界约束引起的转动弯矩,那么转动弯矩引起的弯曲约束应力<maths num="0010"><math><![CDATA[<mrow><msub><msup><mi>&sigma;</mi><mo>&prime;</mo></msup><mi>M</mi></msub><mo>=</mo><mfrac><mrow><mi>h</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>T</mi></msub><mo>-</mo><msub><mi>&theta;</mi><mi>U</mi></msub><mo>)</mo></mrow></mrow><mi>L</mi></mfrac></mrow>]]></math><img file="FDA0001196562870000028.GIF" wi="317" he="94" /></maths>其中,σ′<sub>M</sub>表示转动约束应力;步骤六、温度自应力计算<maths num="0011"><math><![CDATA[<mrow><msub><msup><mi>&sigma;</mi><mo>&prime;</mo></msup><mrow><mi>S</mi><mi>E</mi></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>E&epsiv;</mi><mi>T</mi><mi>N</mi></msubsup><mo>+</mo><msubsup><mi>E&epsiv;</mi><mi>T</mi><mi>M</mi></msubsup><mo>-</mo><msub><mi>&sigma;</mi><mrow><mi>R</mi><mi>T</mi></mrow></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001196562870000029.GIF" wi="670" he="63" /></maths>车载应力计算σ′<sub>V</sub>=Eε<sub>VLi</sub>步骤七、从有限元模型中获取自重应力σ′<sub>G</sub>;步骤八、大跨桥梁在运营期间总应力计算:σ′<sub>Total</sub>=σ′<sub>N</sub>+σ′<sub>M</sub>+σ′<sub>V</sub>其中,σ′<sub>Total</sub>为总应力。
地址 210096 江苏省南京市四牌楼2号