[发明专利]一种薄壁件靠模装夹贴合状态检测与评估方法有效

专利信息
申请号: 201410476098.8 申请日: 2014-09-16
公开(公告)号: CN104385057A 公开(公告)日: 2015-03-04
发明(设计)人: 刘海波;王永青;郭东明;盛贤君;杨睿;马跃;李亚鹏;王克欣;廉盟 申请(专利权)人: 大连理工大学
主分类号: B23Q17/00 分类号: B23Q17/00
代理公司: 大连理工大学专利中心 21200 代理人: 关慧贞
地址: 116024 辽*** 国省代码: 辽宁;21
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明一种薄壁件靠模装夹贴合状态检测与评估方法属于检测技术领域,特别涉及一种薄壁件靠模装夹贴合状态检测与评估方法。该方法中,先将薄壁件靠模装夹到支撑夹具上,并吊装至数控加工机床上;接着,数控机床驱动喷流水浸超声探头按照预先规划的测量路径,分段线性定向扫描测量被测件,上层测量系统自动采集、存储测点波形数据;最后,通过反射子波辨识与峰值提取、灰度比较、数据滤波、灰度匀化等一系列数据处理操作,完成贴合状态的数字化评估。本发明将喷流水浸超声在机检测与数字化灰度评估相结合的方法,实现了薄壁件靠模装夹贴合状态的准确判断,适于复杂曲面薄壁件靠模装夹贴合状态的在机检测与评估,操作简单、方法可靠。
搜索关键词: 一种 薄壁 件靠模装夹 贴合 状态 检测 评估 方法
【主权项】:
一种薄壁件靠模装夹贴合状态检测与评估方法,其特征在于,检测与评估方法中,首先将薄壁件靠模装夹到支撑夹具上,并吊装至数控加工机床上;接着,数控机床驱动超声探头按照预先规划的测量路径分段线性定向扫描测量被测件,上层测量系统自动采集、存储测点波形数据;最后,通过反射子波峰值提取、灰度比较、数据滤波、灰度匀化等一系列数据处理操作,完成贴合状态评估;薄壁件靠模装夹贴合状态检测与评估方法的具体步骤如下:第一步,测量就位将薄壁件(4)放置在支撑夹具(3)上,带螺栓的第一压板组(5)压紧薄壁件(4);将支撑夹具(3)吊装至工作台(1)上,带螺栓的第二压板组件(2)压紧支撑夹具(3)的法兰,完成薄壁件(4)的测量就位;第二步,超声线性分段定向扫描测量首先,将超声测头(6)安装至机床主轴(7)上;超声测头(6)在机床主轴(7)驱动下运动至第一截面轨迹L1内第一线性运动段τ11的起始控制点C11,并自动调整超声测头(6)的中心线(a)与第一线性运动段τ11的平均法矢N11重合;按如下公式计算第一线性运动段τ11的平均法矢N11<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>N</mi><mn>11</mn></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>m</mi><mn>11</mn></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>m</mi><mn>11</mn></msub></munderover><msubsup><mi>n</mi><mn>11</mn><mi>k</mi></msubsup></mtd></mtr><mtr><mtd><munder><mi>max</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn><mo>~</mo><msub><mi>m</mi><mn>11</mn></msub></mrow></munder><mo>{</mo><msubsup><mi>&theta;</mi><mn>11</mn><mi>k</mi></msubsup><mo>|</mo><msubsup><mi>&theta;</mi><mn>11</mn><mi>k</mi></msubsup><mo>=</mo><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msubsup><mi>n</mi><mn>11</mn><mi>k</mi></msubsup><mo>&CenterDot;</mo><msub><mi>N</mi><mn>11</mn></msub><mo>)</mo></mrow><mo>}</mo><mo>&le;</mo><msub><mi>&delta;</mi><mi>&theta;</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>式中,为第一线性运动段τ11第i点的法矢,与N11的夹角,δθ为法矢阈值,m11为第一线性运动段τ11离散点数;在第一截面轨迹L1对应的数字化测量面内,超声测头(6)连续运动至第二线性运动段τ21的起始控制点C21,并使得超声测头(6)的中心线(a)与第二线性运动段τ21的平均法矢N21重合;超声测头(6)沿第一截面轨迹L1内的各线性运动段依次扫描测量,同时上位机的数据采集系统高频存储各测点坐标及反射回波数据,测量获得第一截面轨迹L1的测量子集Ω1={(pi1,Ui1),i∈[1,r]},pi1为第一截面轨迹L1第i测点坐标,Ui1为第一截面轨迹L1第i测点反射回波,r为截面轨迹内的测点数;接着,超声测头(6)运动至第二截面轨迹L2第一线性运动段τ12的起始控制点C12,并使得超声测头(6)的中心线(a)与第二截面轨迹L2的第一线性运动段τ12的平均法矢N12重合;超声测头(6)沿第二截面轨迹L2内的各线性运动段依次扫描测量,获得第二截面轨迹测L2的量子集Ω2;最后,超声测头(6)沿Z字形扫描路径对薄壁件(4)双向往复扫描测量,测量获得薄壁件(4)的测量总集Ω={Ωj,j∈[1,s]},s为截面轨迹总数,Ωj为第j截面轨迹Lj的测量子集;第三步,灰度计算取第j截面轨迹Lj第i测点的反射回波Uij={uk,k∈[1,tij]},进行第二反射子波区间辨识与峰值提取;uk为第k点反射幅值,tij为第j截面轨迹Lj第i测点反射回波Uij的离散数;首先,遍历反射回波内所有数据,若uk<δu,则uk被判定为波谷点,并记录当前波谷序号Tl=k,完成波谷数值提取,相邻两个波谷点决定着一个反射子波;δu为波谷判定阈值;其次,通过波谷序号计算反射子波区间宽度,Wl=Tl+1‑Tl,若Wl<δW,则Wl=Wl+Wl+1,即将Wl对应的第l反射子波与Wl+1对应的第(l+1)反射子波合并,完成反射子波区间辨识;δW为反射子波辨识阈值;提取第j截面轨迹Lj第i测点的第二反射子波峰值Vij=max{uk,k∈[T2,T3]};对第j截面轨迹Lj第i测点的第二反射子波峰值Vij进行贴合灰度计算,模型如下,<mrow><mi>Gray</mi><msub><mi>U</mi><mi>ij</mi></msub><mo>=</mo><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>ij</mi></msub><msub><mi>V</mi><mn>0</mn></msub></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>&times;</mo><mn>255</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>式中,GrayUij为第j截面轨迹第i测点的贴合灰度值,V0为贴合标定值;按上述操作,对第j测量子集Ωj中的反射回波逐点灰度处理,得到第j灰度子集Gj={(pij,GrayUij),i∈[1,r]},pij为第j截面轨迹第i点坐标;对测量总集Ω,按截面轨迹逐条灰度处理,得到薄壁件(4)的灰度总集G={Gj,j∈[1,s]};第四步,灰度数据模糊加权均值滤波构造第j灰度子集Gj第i测点的模糊滤波窗口Aij={GrayUij,m|m∈[i‑Δ,i+Δ]},Δ为模糊滤波窗口宽度;接着,按如下模型滤波,<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&mu;</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msup><mrow><mo>(</mo><mrow><mo>(</mo><mi>Gray</mi><msub><mi>U</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>-</mo><mi>Gray</mi><msub><mi>U</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><mi>b</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>Gray</mi><msub><mi>U</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&Delta;</mi><mo>+</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mi>i</mi><mo>-</mo><mi>&Delta;</mi></mrow><mrow><mi>i</mi><mo>+</mo><mi>&Delta;</mi></mrow></munderover><mi>Gray</mi><msub><mi>U</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>=</mo><msub><mi>&mu;</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mi>i</mi><mo>-</mo><mi>&Delta;</mi></mrow><mrow><mi>i</mi><mo>+</mo><mi>&Delta;</mi></mrow></munderover><msub><mi>&mu;</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>Gray</mi><msubsup><mi>U</mi><mi>ij</mi><mo>&prime;</mo></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mi>i</mi><mo>-</mo><mi>&Delta;</mi></mrow><mrow><mi>i</mi><mo>+</mo><mi>&Delta;</mi></mrow></munderover><msub><mi>w</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub><mi>Gray</mi><msub><mi>U</mi><mrow><mi>ij</mi><mo>,</mo><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>式中,μij,m为GrayUij,m隶属度值,b为模糊尺度因子,GrayU0为滤波窗口Aij内的灰度均值,wij,m为GrayUij,m权值,GrayU′ij为第(i,j)滤波灰度值;按上述操作,对第j灰度子集Gj中的灰度数据,逐点处理,得到滤波子集Fj;对灰度总集G中的灰度数据,按截面轨迹逐条处理,得到薄壁件(4)的滤波总集F={Fj,j∈[1,s]};第五步,装夹贴合状态数字化评估首先,在薄壁件(4)被测坐标范围内进行等参数坐标密化,得到坐标密化点集Q={qIJ,I∈[1,R],J∈[1,S]};其中,qIJ为第(I,J)密化点,R为截面轨迹内密化点数,S为截面轨迹密化数;接着,对滤波总集F进行四边剖分处理,并进行序列化标记,记由四个测点坐标pij、p(i+1)j、pi(j+1)和p(i+1)(j+1)确定的剖腔为第(i,j)微元剖腔Θij;利用密化点与微元剖腔的坐标,判断二者空间位置关系;若密化点位于两微元剖腔的公共边界和多微元剖腔的公共点上,则将密化点判断在低序列微元剖腔中;利用高维插值技术计算各密化点的灰度值;以第(I,J)密化点qIJ为例,经判断,第(I,J)密化点qIJ位于第(i,j)微元剖腔Θij中,则按如下公式计算第(I,J)密化灰度值<mrow><mo>|</mo><mfenced open='' close='' separators='|'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>p</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>(</mo><msub><mi>X</mi><mi>IJ</mi></msub><mo>-</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>p</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>-</mo><msub><mi>p</mi><mi>ij</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>(</mo><msub><mi>Y</mi><mi>IJ</mi></msub><mo>-</mo><msub><mi>y</mi><mi>ij</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>p</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>-</mo><msub><mi>p</mi><mi>ij</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mtd><mtd><mrow><mo>(</mo><msub><mi>Z</mi><mi>IJ</mi></msub><mo>-</mo><msub><mi>z</mi><mi>ij</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>Gray</mi><msubsup><mi>U</mi><mi>ij</mi><mo>&prime;</mo></msubsup></mtd><mtd><mi>Gray</mi><msubsup><mi>U</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>j</mi></mrow><mo>&prime;</mo></msubsup></mtd><mtd><mi>Gray</mi><msubsup><mi>U</mi><mrow><mi>i</mi><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mo>&prime;</mo></msubsup></mtd><mtd><mi>Gray</mi><msubsup><mi>U</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mo>&prime;</mo></msubsup></mtd><mtd><mi>Gray</mi><msub><mover><mi>U</mi><mo>~</mo></mover><mi>IJ</mi></msub></mtd></mtr></mtable><mrow><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></mfenced></mrow>式中,p(i+1)j为第j截面轨迹第(i+1)点坐标,pi(j+1)为第(j+1)截面轨迹第i点坐标,p(i+1)(j+1)为第(j+1)截面轨迹第(i+1)点坐标;GrayU′(i+1)j为第(i+1,j)滤波灰度值,GrayU′i(j+1)为第(i,j+1)滤波灰度值,GrayU′(i+1)(j+1)为第(i+1,j+1)滤波灰度值;XIJ、YIJ、ZIJ分别为第(I,J)密化点qIJ的X坐标、Y坐标和Z坐标;xij、yij、zij分别为第(i,j)测点pij的x坐标、y坐标和z坐标;‖p(i+1)j‑pi,j2、‖pi(j+1)‑pij2与‖p(i+1)(j+1)‑pij2分别表示p(i+1)j‑pi,j、pi(j+1)‑pij与p(i+1)(j+1)‑pij的二范数,用以计算空间两点的Hausdorff距离;逐点依次计算,得到坐标密化点集Q内各密化点的灰度值;然后,对密化点的贴合状态灰度值进行匀化处理;若当前点的贴合状态灰度值均不大于临近点的贴合状态灰度值,则按如下公式匀化处理,<mrow><mi>Gray</mi><msub><mover><mi>U</mi><mo>~</mo></mover><mi>IJ</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><mrow><mo>(</mo><mi>Gray</mi><msub><mover><mi>U</mi><mo>~</mo></mover><mrow><mrow><mo>(</mo><mi>I</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>J</mi></mrow></msub><mo>+</mo><mi>Gray</mi><msub><mover><mi>U</mi><mo>~</mo></mover><mrow><mi>I</mi><mrow><mo>(</mo><mi>J</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>+</mo><mi>Gray</mi><msub><mover><mi>U</mi><mo>~</mo></mover><mrow><mrow><mo>(</mo><mi>I</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>J</mi></mrow></msub><mo>+</mo><mi>Gray</mi><msub><mover><mi>U</mi><mo>~</mo></mover><mrow><mi>I</mi><mrow><mo>(</mo><mi>J</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>式中,为第(I‑1,J)密化灰度值,为第(I,J‑1)密化灰度值,为第(I+1,J)密化灰度值,为第(I,J+1)密化灰度值;最后,设第(I,J)密化点qIJ四边邻域网格内贴合状态具有一致性,贴合状态灰度值均为逐点填充,输出贴合状态灰度云图,完成薄壁件(4)靠模装夹贴合状态数字化评估。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于大连理工大学,未经大连理工大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201410476098.8/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top