发明名称 基于Kaczmarz代数迭代重建方法的桥梁移动车辆荷载识别方法
摘要 本发明公开了一种基于Kaczmarz代数迭代重建方法的桥梁移动车辆荷载识别方法,包括以下步骤:(1)、在桥面预设测点位置均布置光电传感器和加速度计,(2)、建立桥梁的简化物理力学模型(3)、通过卷积积分在时域求解                                               <img file="dest_path_image002.GIF" wi="48" he="32" />,得到形如<img file="dest_path_image004.GIF" wi="48" he="22" />的系统方程,则<img file="dest_path_image006.GIF" wi="16" he="16" />为已知的系统矩阵,<img file="dest_path_image008.GIF" wi="14" he="16" />为已知的桥面响应,<img file="dest_path_image010.GIF" wi="17" he="22" />即为所求的移动车辆荷载;(4)利用Kaczmarz代数迭代重建方法通过选取合适的迭代次数<img file="dest_path_image012.GIF" wi="14" he="16" />即可由桥梁加速度响应识别桥面移动车辆荷载<img file="286055dest_path_image010.GIF" wi="17" he="22" />。本发明提出的Kaczmarz代数迭代重建方法只需简单迭代即可实现移动荷载识别,迭代公式中不包含线性方程组以外的其他辅助向量,输入方便,理论清晰,易于推广,非常适用于移动荷载识别。
申请公布号 CN104598931A 申请公布日期 2015.05.06
申请号 CN201510062912.6 申请日期 2015.02.06
申请人 华北水利水电大学 发明人 陈震
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 郑州立格知识产权代理有限公司 41126 代理人 田小伍
主权项 一种基于Kaczmarz代数迭代重建方法的桥梁移动车辆荷载识别方法,其特征在于:包括以下步骤:(1)、在桥面预设测点位置布置光电传感器和加速度计,由光电传感器测量获得移动车辆荷载的行驶速度,加速度计测量获得桥梁在移动车辆荷载作用下的加速度;(2) 建立桥梁的简化物理力学模型,取桥梁长度为<i>L</i>,桥面移动车辆荷载<img file="dest_path_image002.GIF" wi="17" he="22" />以均匀速度<i>c</i>沿着桥面移动,桥梁单位长度质量为<i>ρ</i>,考虑粘性阻尼并取阻尼系数为<i>C</i>,简支梁抗弯刚度为EI;桥面上<i>t</i>时刻移动车辆荷载动力学方程如下形式:<img file="dest_path_image004.GIF" wi="372" he="46" />上式中:<img file="dest_path_image006.GIF" wi="35" he="21" />为桥梁的第<i>n</i>阶模态振型;<img file="dest_path_image008.GIF" wi="45" he="27" />为<img file="590988dest_path_image006.GIF" wi="35" he="21" />的二阶导数,<img file="dest_path_image010.GIF" wi="45" he="27" />为<img file="281995dest_path_image006.GIF" wi="35" he="21" />的一阶导数;<img file="dest_path_image012.GIF" wi="22" he="28" />为粘性阻尼比,<img file="dest_path_image014.GIF" wi="88" he="52" />;C为阻尼系数,不是一个定值,和桥梁自身特征及材料属性有关,可以根据仪器测得;<img file="dest_path_image016.GIF" wi="20" he="23" />为圆频率,<img file="dest_path_image018.GIF" wi="98" he="44" />;<img file="dest_path_image020.GIF" wi="39" he="23" />为桥面移动车辆荷载模态表达式,<img file="dest_path_image022.GIF" wi="140" he="39" />;<img file="dest_path_image024.GIF" wi="33" he="20" />表示车辆荷载<img file="dest_path_image026.GIF" wi="19" he="25" />在t时刻的值;(3)通过卷积积分在时域求解<img file="dest_path_image028.GIF" wi="31" he="21" />:<img file="dest_path_image030.GIF" wi="209" he="52" />这里<img file="dest_path_image032.GIF" wi="169" he="46" />(<img file="dest_path_image034.GIF" wi="36" he="20" />),<img file="dest_path_image036.GIF" wi="109" he="32" />;<img file="dest_path_image038.GIF" wi="47" he="33" />没有具体含义,只是代替其右端项的公式;<img file="dest_path_image040.GIF" wi="24" he="30" />是考虑粘性阻尼后的圆频率;<img file="dest_path_image042.GIF" wi="17" he="18" />是积分符号;由卷积积分得到桥梁上<i>x</i>处,<i>t</i>时刻的位移<i>v</i>为:<img file="dest_path_image044.GIF" wi="440" he="50" />由卷积积分得到桥梁上<i>x</i>处,<i>t</i>时刻的加速度<img file="dest_path_image046.GIF" wi="12" he="17" />为:<img file="dest_path_image048.GIF" wi="330" he="46" />这里<img file="dest_path_image050.GIF" wi="402" he="46" /><img file="dest_path_image052.GIF" wi="39" he="30" />是<img file="556112dest_path_image038.GIF" wi="47" he="33" />的二阶导数;加速度<img file="658060dest_path_image046.GIF" wi="12" he="17" />写成离散形式:<img file="dest_path_image054.GIF" wi="454" he="52" />其中第<i>n</i>阶模态的加速度响应为:<img file="dest_path_image056.GIF" wi="442" he="52" />取<img file="dest_path_image058.GIF" wi="118" he="45" />,<img file="dest_path_image060.GIF" wi="149" he="28" />,<img file="dest_path_image062.GIF" wi="140" he="42" />取<img file="dest_path_image064.GIF" wi="69" he="42" />,<img file="dest_path_image066.GIF" wi="25" he="25" />表示测量响应的采集数据个数,如果有<i>N</i>个测点的加速度响应(<i>N</i>&gt;1) ,加速度<img file="86637dest_path_image046.GIF" wi="12" he="17" />离散形式表示为:<img file="dest_path_image068.GIF" wi="503" he="104" />设在上桥和下桥时刻<img file="dest_path_image070.GIF" wi="61" he="22" />,<img file="dest_path_image072.GIF" wi="74" he="24" />,则有<img file="dest_path_image074.GIF" wi="57" he="22" />;加速度<img file="253438dest_path_image046.GIF" wi="12" he="17" />离散形式改写为:<img file="dest_path_image076.GIF" wi="503" he="82" />由加速度<img file="631330dest_path_image046.GIF" wi="12" he="17" />离散形式得到系统矩阵<i>A<sub>n</sub></i>,可求解如下方程组得到荷载向量<i>f</i> :<img file="dest_path_image078.GIF" wi="155" he="44" />(4)将<img file="232075dest_path_image078.GIF" wi="155" he="44" />看做形如<img file="dest_path_image080.GIF" wi="49" he="22" />的系统方程,则<img file="dest_path_image082.GIF" wi="17" he="18" />为步骤(2)中已知的系统矩阵,<img file="dest_path_image084.GIF" wi="14" he="20" />为已知的桥面响应(加速度响应),<img file="633101dest_path_image002.GIF" wi="17" he="22" />即为所求的移动车辆荷载;(5)取<img file="dest_path_image086.GIF" wi="16" he="25" />为方程<img file="164445dest_path_image080.GIF" wi="49" he="22" />右端项的第<img file="dest_path_image088.GIF" wi="10" he="18" />列向量,<img file="dest_path_image090.GIF" wi="21" he="26" />为系统矩阵<img file="334526dest_path_image082.GIF" wi="17" he="18" />的第<img file="106173dest_path_image088.GIF" wi="10" he="18" />行转置向量;Kaczmarz代数迭代重建方法第<img file="dest_path_image092.GIF" wi="14" he="20" />步迭代解表示为:<img file="dest_path_image094.GIF" wi="237" he="64" />通过选取合适的迭代次数<img file="745227dest_path_image092.GIF" wi="14" he="20" />即可由桥梁加速度响应识别桥面移动车辆荷载<img file="627733dest_path_image002.GIF" wi="17" he="22" />的第<img file="980217dest_path_image092.GIF" wi="14" he="20" />步迭代解<img file="dest_path_image096.GIF" wi="30" he="25" />。
地址 450011 河南省郑州市北环路36号