发明名称 基于双目立体视觉的结构模态参数识别系统
摘要 本发明公开了一种基于双目立体视觉识别结构模态参数的系统。其包括:搭建双目立体视觉系统,标定获取两个摄像单元的内外参数,进行三维立体校正;对被测结构表面建立易识别的关键点;加载冲击力让被测结构自由运动,获取同步的图像序列文件,依次对每一帧图像识别所有关键点的三维坐标,最终得到整个测试时间内所有关键点的三维运动数据;通过所有关键点的三维运动数据使用特征系统实现算法进行模态识别,获得被测结构的模态参数,输出结果。与传统技术相比,本发明能够非接触地或对被测结构施加较小影响地测量被测结构的模态参数,其成本低、速度快、精度高,尤其适合不方便布置传感器的大型质轻结构的模态参数测量。
申请公布号 CN103323209B 申请公布日期 2016.04.06
申请号 CN201310272079.9 申请日期 2013.07.02
申请人 清华大学 发明人 何家玮;顾永鹏;赵治华;任革学;马彦会;杨程;王宁羽
分类号 G01M7/08(2006.01)I;G06T7/00(2006.01)I 主分类号 G01M7/08(2006.01)I
代理机构 代理人
主权项 一种使用基于双目立体视觉的结构模态参数识别系统进行模态测量的方法,所述基于双目立体视觉的结构模态参数识别系统包括摄像装置、数据采集和处理装置;所述摄像装置包括两个摄像单元和一个触发信号发生装置,所述两个摄像单元基本以平行方式设置,左右相隔一定距离,固定在支撑单元上;两个摄像单元能够在所述触发信号发生装置的同步下对被测结构进行同步拍摄,从而生成同步图像序列文件;所述数据采集和处理装置能够从所述摄像装置获得所述同步图像序列文件,并进行处理,获得所述被测结构的结构模态参数;其特征在于,所述方法包括如下步骤:步骤一:对两个摄像单元进行标定,获得两个摄像单元各自的内参数和它们之间的外参数;通过三维立体校正,获得图像的消除畸变的行对齐映射和重投影矩阵;步骤二:建立被测结构表面的关键点,加载冲击载荷让被测结构自由运动,同时使用标定完的双摄像单元摄像被测结构受冲击载荷后的自由运动,获得同步图像序列文件;步骤三:处理摄像单元的同步图像序列文件,获得测试时间内,关键点的三维运动数据;步骤四:通过所有关键点的三维运动数据,以特征系统实现算法识别被测结构的模态;步骤五:根据识别结果,图形化显示各阶模态的无阻尼频率、阻尼比、阵型、模态置信度(MAC)和模态相位共线性(MPC);所述步骤一中,记Q为某点在摄像单元坐标系中的三维坐标;<img file="FSB0000144901390000011.GIF" wi="59" he="72" />为它在此摄像单元上的以齐次坐标表示的理想成像位置,x<sub>p</sub>和y<sub>p</sub>是图像上的水平方向x坐标和竖直方向y坐标,其中水平方向以向右为正,竖直方向以向上为正;<img file="FSB0000144901390000025.GIF" wi="51" he="65" />为考虑成畸变后真实的位置,x<sub>d</sub>和y<sub>d</sub>是图像上的x坐标和y坐标;M为内参数矩阵,f<sub>x</sub>,f<sub>y</sub>,c<sub>x</sub>,c<sub>y</sub>为摄像单元的x向焦距、y向焦距和光轴中心坐标;F为畸变向量;T,R为右摄像单元相对左摄像单元的平移向量和旋转向量;s是一个尺度比例:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>Q</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>X</mi></mtd></mtr><mtr><mtd><mi>Y</mi></mtd></mtr><mtr><mtd><mi>Z</mi></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>p</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>x</mi><mi>p</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>p</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>d</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>x</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>d</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>M</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>f</mi><mi>x</mi></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>c</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>f</mi><mi>y</mi></msub></mtd><mtd><msub><mi>c</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FSB0000144901390000021.GIF" wi="971" he="215" /></maths>F=[k<sub>1</sub> k<sub>2</sub> p<sub>1</sub> p<sub>2</sub> k<sub>3</sub>]T=[t<sub>x</sub> t<sub>y</sub> t<sub>z</sub>]<sup>T</sup>R=[R<sub>1</sub> R<sub>2</sub> R<sub>3</sub>]这些参数间满足关系式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>~</mo></mover><mi>p</mi></msub><mo>=</mo><mi>s</mi><mi>M</mi><mi>Q</mi></mrow>]]></math><img file="FSB0000144901390000022.GIF" wi="220" he="89" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>x</mi><mi>p</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>p</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><msup><mi>r</mi><mn>2</mn></msup><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><msup><mi>r</mi><mn>4</mn></msup><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><msup><mi>r</mi><mn>6</mn></msup><mo>)</mo></mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>x</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>d</mi></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mn>2</mn><msub><mi>p</mi><mn>1</mn></msub><msub><mi>x</mi><mi>d</mi></msub><msub><mi>y</mi><mi>d</mi></msub><mo>+</mo><msub><mi>p</mi><mn>2</mn></msub><mrow><mo>(</mo><mrow><msup><mi>r</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><msub><mi>x</mi><mi>d</mi></msub><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mn>2</mn><msub><mi>p</mi><mn>2</mn></msub><msub><mi>x</mi><mi>d</mi></msub><msub><mi>y</mi><mi>d</mi></msub><mo>+</mo><msub><mi>p</mi><mn>1</mn></msub><mrow><mo>(</mo><mrow><msup><mi>r</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><msub><mi>y</mi><mi>d</mi></msub><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FSB0000144901390000023.GIF" wi="1193" he="178" /></maths>通过标定获得内外参数:M<sub>left</sub>,M<sub>right</sub>,F<sub>left</sub>,F<sub>right</sub>,R,T。
地址 100084 北京市海淀区清华大学