发明名称 悬索桥施工过程中加劲梁测量坐标修正方法
摘要 本发明公开了一种悬索桥施工过程中考虑风致变形的加劲梁测量坐标修正方法,首先测量出含有风致变形的加劲梁坐标,并同步测试桥位风速、风攻角和风向角。利用指数型风剖面函数求得各点的风速,并将风速分解成横桥向风速和纵桥向风速。将加劲梁、桥塔、主缆和吊杆风荷载施加在有限元模型上,并计算出风致变形值;将其在加劲梁坐标测量值中扣除,即可获得零风速时的加劲梁坐标,即真实的加劲梁坐标。通过该方法可以实现对悬索桥加劲梁坐标测量值进行风致变形的修正,因此测量坐标不用再等待无风或微风时段,可以避免延误工期。
申请公布号 CN104715149A 申请公布日期 2015.06.17
申请号 CN201510113238.X 申请日期 2015.03.16
申请人 东南大学 发明人 张文明
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 张学彪
主权项 一种悬索桥施工过程中加劲梁测量坐标修正方法,所述悬索桥包括桥塔、吊杆和主缆,其特征在于,包括以下步骤:1)、根据桥位地形确定地表粗糙度系数γ;2)、识别出桥梁加劲梁断面的三分力系数以及桥塔、吊杆和主缆的阻力系数;3)、在桥塔顶部安装风速测试仪,测量桥塔顶部的风速、风攻角和风向角;4)、测量加劲梁坐标(X<sub>0</sub>,Y<sub>0</sub>,Z<sub>0</sub>),并同步记录桥塔顶部的风速U、风攻角α<sub>0</sub>和风向角β;5)、利用下式将风速U分解并利用指数型风速剖面求得加劲梁处的纵桥向风速U<sub>Si</sub>和横桥向风速U<sub>Hi</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>U</mi><mi>Si</mi></msub><mo>=</mo><mi>U</mi><msup><mrow><mo>(</mo><mfrac><msub><mi>H</mi><mi>i</mi></msub><msub><mi>H</mi><mi>T</mi></msub></mfrac><mo>)</mo></mrow><mi>&gamma;</mi></msup><mi>cos</mi><mi>&beta;</mi></mtd></mtr><mtr><mtd><msub><mi>U</mi><mi>Hi</mi></msub><mo>=</mo><mi>U</mi><msup><mrow><mo>(</mo><mfrac><msub><mi>H</mi><mi>i</mi></msub><msub><mi>H</mi><mi>T</mi></msub></mfrac><mo>)</mo></mrow><mi>&gamma;</mi></msup><mi>sin</mi><mi>&beta;</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000682437150000011.GIF" wi="475" he="346" /></maths>式中,U<sub>Si</sub>为加劲梁处的纵桥向风速,U<sub>Hi</sub>为加劲梁处的横桥向风速,H<sub>i</sub>为加劲梁的高度,H<sub>T</sub>为桥塔顶部的高度,U为桥塔顶部的风速,γ为地表粗糙度系数,β为桥塔顶部的风向角;6)、计算纵桥向风速引起的作用在加劲梁上的摩擦力:F=ηSLU<sub>Si</sub>式中,η为摩擦系数;S为加劲梁截面周长;L为加劲梁长度;U<sub>Si</sub>为加劲梁处的顺桥向风速;7)、分别计算纵桥向风速和横桥向风速作用在主缆、桥塔和吊杆上的阻力:D<sub>Si</sub>=ρU<sup>2</sup>C<sub>H</sub>(β)w/2D<sub>Hi</sub>=ρU<sup>2</sup>C<sub>V</sub>(β)w/2式中,D<sub>Si</sub>和D<sub>Hi</sub>分别表示桥塔、吊杆或主缆的顺桥向阻力和横桥向阻力;C<sub>H</sub>(β)和C<sub>V</sub>(β)分别表示桥塔、吊杆或主缆的顺桥向阻力系数和横桥向阻力系数;w表示主缆、桥塔或吊杆的宽度;ρ为空气密度,β为桥塔顶部的风向角;8)、建立桥梁的有限元模型,在自重作用下非线性求解,提取各加劲梁单元扭转角,计算该状态下的三分力系数,此时,加劲梁的有效攻角均等于初始攻角α<sub>0</sub>;9)、计算在横桥向风速引起的作用在加劲梁单位长度的横向风荷载P<sub>H</sub>、竖向风荷载P<sub>V</sub>和扭转力矩P<sub>M</sub>:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>H</mi></msub><mo>=</mo><msubsup><mi>&rho;H</mi><mi>Hi</mi><mn>2</mn></msubsup><msub><mi>C</mi><mi>H</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mi>h</mi><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><msub><mi>P</mi><mi>V</mi></msub><mo>=</mo><msubsup><mi>&rho;U</mi><mi>Hi</mi><mn>2</mn></msubsup><msub><mi>C</mi><mi>V</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mi>b</mi><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><msub><mi>P</mi><mi>M</mi></msub><mo>=</mo><msubsup><mi>&rho;U</mi><mi>Hi</mi><mn>2</mn></msubsup><msub><mi>C</mi><mi>M</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><msup><mi>b</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000682437150000021.GIF" wi="508" he="255" /></maths>式中,C<sub>H</sub>(α)、C<sub>V</sub>(α)、C<sub>M</sub>(α)分别表示在有效攻角下加劲梁的阻力、升力、升力矩系数,所谓有效攻角是指静风初始攻角与静风作用引起的加劲梁扭转角之和;h和b分别表示加劲梁的侧向投影高度和宽度;ρ为空气密度;U<sub>Hi</sub>为横桥向风速;10)、在桥梁各单元上施加步骤6)、步骤7)和步骤9)计算的静风荷载,采用牛顿迭代法进行结构几何非线性求解,获得收敛解;11)、提取桥梁加劲梁各单元扭转角,检查扭转角的欧几里得范数是否小于允许误差ε<sub>k</sub>;其中,所述的扭转角的欧几里得范数用下式表达:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mrow><mo>{</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>Na</mi></munderover><msup><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>Na</mi></munderover><msubsup><mi>&theta;</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow><mn>2</mn></msubsup></mrow></mfrac><mo>}</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>&le;</mo><msub><mi>&epsiv;</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0000682437150000022.GIF" wi="509" he="322" /></maths>式中,N<sub>a</sub>为加劲梁单元总数;θ为扭转角;i为当前荷载步编号;j为梁单元序号;ε<sub>k</sub>为允许误差;12)、如果欧几里得范数大于允许误差ε<sub>k</sub>,则重复步骤9)‑11);如果欧几里得范数小于等于允许误差ε<sub>k</sub>,说明本级风速计算结果收敛,输出测量点的位移值(x<sub>wind</sub>,y<sub>wind</sub>,z<sub>wind</sub>);13)、在坐标测量值(X<sub>0</sub>,Y<sub>0</sub>,Z<sub>0</sub>)中将步骤12)计算得到的测量点的位移值扣除:X<sub>real</sub>=X<sub>0</sub>‑x<sub>wind</sub>Y<sub>real</sub>=Y<sub>0</sub>‑y<sub>wind</sub>Z<sub>real</sub>=Z<sub>0</sub>‑z<sub>wind</sub>即得到真实的加劲梁测量坐标(X<sub>real</sub>,Y<sub>real</sub>,Zz<sub>real</sub>)。
地址 210096 江苏省南京市四牌楼2号