专利名称：：基于ccd图像特征的表面疲劳裂纹检测方法
技术领域：
：本发明涉及一种表面疲劳裂纹的检测方法，尤其是一种综合了数字图像技术和MATLAB分析技术的表面疲劳裂纹的检测方法，具体地说是一种基于CCD图像特征的表面疲劳裂纹检测方法。
背景技术：
：随着现代科学技术的迅猛发展，工业设备要求在高温、高压、高速等条件下使用的场合越来越多，承受的交变应力越来越高，疲劳破坏问题日益突出。机械零件的疲劳破坏是导致其结构强度破坏的主要因素，而疲劳裂纹是疲劳破坏的主要形式。在交变应力作用下构件会萌生微观裂纹，尤其是主要受力构件。若没有得到及时处理，这些微裂纹就可能扩展为宏观裂纹，而宏观裂纹一旦扩展到其临界值会发生骤然失稳，而引起构件的疲劳断裂，造成不可想象的后果。因此，开展疲劳裂纹扩展规律的研究，对我国现代工业的发展有重要的指导意义，并带来极大的社会效益。疲劳裂纹萌生与扩展的测试及分析是金属结构疲劳设计与寿命预测的主要任务之一。曾有学者提出根据材料的性能及其热处理方式等预测疲劳裂纹的da/dN曲线，进而预测其疲劳寿命，但这种估算方法精度并不高。目前，国内外已有的疲劳裂纹检测方法(1)离位方式，主要为目测法，两次测试时间间隔很长，忽略了疲劳裂纹扩展的某些信息，因此，离位方式不能实时连续的跟踪、记录裂纹扩展的全过程；(2)停机在位方式，普遍采用的有表面复型技术(R印lica法)、电位法、磁粉法、超声法、射线法、电磁检测法等，这些方法普遍存在设备复杂、需要接触试样、实施繁琐等问题；(3)在线方式，主要采用声发射法、模态声发射法，这两方法设备昂贵、实施困难、对测量环境要求高。图像识别技术具有精度高、对测量环境要求低、非接触式并可实现实时监测等优点。因此将图像识别技术应用于表面疲劳裂纹的检测，可实现无接触检测，而且设备简单易操作、成本低，节约劳动力，精度可靠，为实现加工自动化、智能化奠定基础。
发明内容本发明的目的是针对现有的疲劳裂纹检测方法设备昂贵、实施繁琐、需要接触试样等问题，提出一种基于CCD图像特征的表面疲劳裂纹检测方法。因为，申请人发现通过对表面疲劳裂纹图像进行预处理，可以获得其长度特征和宽度特征，进而可以获得表面疲劳裂纹的长度曲线(a-N曲线)、宽度曲线(w-N曲线)、长度扩展速率曲线(da/dN-N曲线)和宽度扩展速率曲线(dw/dN-N曲线)。有鉴于此，本发明提出了以下的技术解决方案—种表面疲劳裂纹的图像检测方法，其特征主要包括以下步骤第一步对试样进行疲劳裂纹扩展试验；对试样进行疲劳裂纹扩展试验，利用图像采集设备连续或间隙采集其表面疲劳裂纹的数字图像，将获得的数字图像输入计算机；图像采集的次数可根据采集设备的数据处理能力加以确定，此外，为了验证本发明的精度可同时采用表面复型法对该试样表面疲劳裂纹的扩展进行研究。第二步对表面疲劳裂纹图像进行预处理；将该试样的疲劳裂纹分为微小裂纹和长裂纹，利用MATLAB图像处理技术，分别对微小裂纹和长裂纹图像迸行预处理，通过裁截、縮放操作获得微小裂纹区域图像，通过灰度化、二值化、膨胀腐蚀、开闭操作、细化操作获得长裂纹骨架图像；第三步获取试样表面疲劳裂纹的长度特征和宽度特征参数；(1)利用鼠标点击微小裂纹区域，获得点的坐标值，进而计算出微小裂纹的像素长度(单位pixel);利用树结构存储长裂纹骨架信息，采用递归的方法去除骨架分支，提取出长裂纹主千，进而计算出长裂纹的像素长度(单位pixel)。(2)在所获取的试样的疲劳裂纹扩展过程的全部图像中，以表面疲劳裂纹的扩展初始位置为研究对象，利用鼠标点击获得点的坐标的方法计算表面疲劳裂纹的像素宽度，像素宽度单位为pixel;(3)对图像采集设备进行标定，计算放大因子(单位mm/pbcd),进而计算表面疲劳裂纹的实际长度(单位mm)和实际宽度(单位mm)。实际长度(mm)-像素长度(pixel)X放大因子(mm/pixel),实际宽度(mm)=像素长度(pixel)X放大因子(mm/pixel)。(4)利用指数曲线拟合方法，拟合表面疲劳裂纹扩展的长度曲线(a-N曲线)、宽度曲线(w-N曲线)、长度扩展速率曲线(da/dN-N曲线)和宽度扩展速率曲线(dw/dN-N曲线)。将上述图像检测方法得到的表面疲劳裂纹长度扩展曲线(a-N曲线和da/dN-N曲线)和表面复型法得到的表面疲劳裂纹的长度扩展曲线(a-N曲线和da/dN-N曲线)进行比较，验证本发明的图像捡测方法的精度。通过对误差进行分析得出，误差满足随机分布，分布范围为-4%~4%。比较表征表面疲劳裂纹扩展的长度特征参数和宽度特征参数；将试样的表面疲劳裂纹长度扩展速率曲线(da/dN-N曲线)和宽度扩展速率曲线(dw/dN-N曲线)进行比较，比较结果表明长度扩展特征参数更适合表征表面疲劳裂纹的扩展情况。本发明的有益效果1.图像采集系统拍摄的图像可靠度高，具备成本低、不接触试样和易实现等突出优点。2.充分发挥利用了MATLAB中图像处理工具在提取数字图像特征参数方面的强大优势，对表面疲劳裂纹图像进行有效的预处理，为表面疲劳裂纹长度扩展速率和宽度扩展速率的计算奠定了坚实的基础。3.本发明涉及到的计算方法及图像处理技术易于实现。4.本发明将数据结构、曲线拟合有效应用到表面疲劳裂纹扩展特征参数的计算当中，排除了实验中由于偶然因素而造成的误差，大大提高了数据计算的准确度。5.可实现在线处理，为检测人员提供了有利工具，降低了劳动强度。6.有助于提高对工件疲劳强度、使用寿命进行评价的自动化水平和智能化水平。7.本发明可用于任何材料工件在任何机械加工方式下的表面疲劳裂纹检8.为机械加工参数选择及条件的确定提供了依据。9.为下一步机械加工过程中故障诊断提供了一种可供借鉴且行之有效的方法。10.应用范围广，可用于各种金属、非金属材料加工表面的检测，尤其适用于难加工材料和难以观察到的深孔表面疲劳裂纹的检测。11.有利于降低成本，可减少因微小裂纹而检测不准确出现的将其报废的现象，降低企业材料成本，且检测速度很快，可提高加工速度。图l是利用本发明进行疲劳裂纹扩展试验的试样零件图。图2是本发明中的进行疲劳裂纹扩展试验的试样装夹示意图，图a为夹头的结构示意图，图b为单侧装夹状态示意图。图3是本发明中的图像采集及处理流程示意图。图4是本发明中表面疲劳裂纹长度的表面复型法研究结果。其中图4(a)表示复型法-疲劳裂纹的a-N曲线，图4(b)表示复型法-疲劳裂纹的da/dN-N曲线。图5是本发明中采用的表面疲劳裂纹的图像检测系统流程图。图6是本发明中表面疲劳裂纹长度的图像检测方法研究结果，其中图6(a)为图像识别法-疲劳裂纹的a-N曲线，图6(b)为图像识别法-疲劳裂纹的da/dN-N曲线。图7是本发明中图像检测方法精度验证结果，其中图7(a)a-N曲线对比，图7(b)da/dN-N曲线对比。图8是本发明中表面疲劳裂纹宽度的图像检测方法研究结果，其中图8(a)为图像识别法-疲劳裂纹的w-N曲线，图8(b)为图像识别法-疲劳裂纹的dw/dN-N曲线。图9是本发明中宽度与长度扩展速率曲线比较结果。图10是本发明中长裂纹主干提取相关示意图，其中图10(a)为裂纹骨架示意图，图10(b)为树结构示意图，图IO节点信息结构示意图。具体实施例方式下面结合附图和实施例对本发明作进一步的说明。如图1-4所示。一种表面疲劳裂纹的图像检测方法，它包括以下步骤-第一步对试样进行疲劳裂纹扩展试验并采集数字图像；试样如图l所示，装夹方案如图2所示。利用自行设计的图像采集系统(如图3所示)采集(可在整个试验过程中连续采集或间隔一定周期定期采集)表面疲劳裂纹扩展过程图像，对表面疲劳裂纹采用表面复型法进行分析，获得表面疲劳裂纹扩展各阶段的长度，根据疲劳试验机振动次数和对应的表面疲劳裂纹长度，利用指数曲线拟合得到表面疲劳裂纹的长度扩展曲线(a-N曲线和da/dN-N曲线)，根据表面复型法得到的疲劳裂纹扩展曲线将长度扩展速率小于2X10—4mm/CyCle的归为微小裂纹，长度扩展速大于等于2X10—4mm/cycle的归为长裂纹，如图4所示。第二步，对采集到的数字图像进行处理；如图59所示。表面疲劳裂纹图像的预处理；对于微小裂纹，先裁截包含表面疲劳裂纹的区域，然后进行适当的縮放，获得清晰的表面疲劳裂纹图像，MATLAB图像处理工具箱中利用imcrop()进行裁截，利用imsize()进行縮放；对于长裂纹，(a)采用最大值法进行灰度化，(b)采用区域提取法进行二值化，(c)采用膨胀与腐蚀操作，而后采用开、闭操作，解决表面疲劳裂纹区域不连通的问题，(d)采用细化算法提取表面疲劳裂纹骨架，MATLAB图像处理工具箱中利用imdilate()进行膨胀，imerode()进行腐蚀，bwmo卬h()进行开、闭以及细化操作，。表面疲劳裂纹长度特征研究；(a)对于微小裂纹，利用鼠标点击裂纹上点，同时获得点的坐标，进而计算出微小裂纹的像素长度(pixel);(b)对于长裂纹，利用树结构存储其骨架信息，采用递归的方法去除骨架分支，提取出主干，进而计算出长裂纹的像素长度(pixel);(c)对图像采集系统进行标定，计算出lmm长度的像素长度，获得放大因子(单位mm/pkel),进而计算表面疲劳裂纹的实际长度(单位mm)。实际长度(mm)=像素长度(pixel)X放大因子(mm/pixel)。(d)利用指数曲线拟合方法，拟合表面疲劳裂纹扩展的长度曲线(a-N曲线)和长度扩展速率曲线(da/dN-N曲线)。获取表面疲劳裂纹的宽度特征参数；利用鼠标点击裂纹上点，同时获得点的坐标，计算出表面疲劳裂纹的像素宽度(pixel);计算表面疲劳裂纹的实际宽度(单位mm)，实际宽度(mm)=像素长度(pixel)X放大因子(mm/pkel);利用指数曲线拟合方法，拟合表面疲劳裂纹扩展的宽度曲线(w-N曲线)和宽度扩展速率曲线(dw/dN-N曲线)。对本发明的检测方法的精度进行验证将图像检测方法得到的表面疲劳裂纹长度扩展曲线(a-N曲线和da/dN-N曲线)和表面复型法得到的表面疲劳裂纹的长度扩展曲线(a-N曲线和da/dN-N曲线)进行比较，验证本研究中图像检测方法的精度。通过对误差进行分析得出，误差满足随机分布，分布范围为-4%~4%。比较表征表面疲劳裂纹扩展的长度特征参数和宽度特征参数；将该试样的表面疲劳裂纹长度扩展速率曲线(da/dN-N曲线)和宽度扩展速率曲线(dw/dN-N曲线)进行比较，比较结果表明长度扩展特征参数更适合表征表面疲劳裂纹的扩展情况。图像处理过程中有关参数的说明。(a)最大值法灰度化彩色图像的灰度化就是使彩色图像的R、G、B分量值相等的过程。最大值法使灰度值等于原始图像像素R、G、B三个值中最大的一个，即Gray(灰度值)=max(R，GB)。(b)区域提取法二值化首先，采用鼠标点取的方式获得裂纹区域若干点的灰度值g。，g,，…，g。，取均值记为^^，以^，为阈值对裂纹表面图像进行二值化；而后，提取裂纹区域，MATLAB图像处理工具箱中利用bwselect()。(c)长裂纹主干提取算法对一副疲劳裂纹骨架，如图10所示，对点和边的规定如表13所示。定义如下两个数组I.节点编号数组该数组是一个和裂纹骨架图像大小相等的二维数组，它的每一个元素对应裂纹骨架图像上相应位置的像素点，如果该像素点是节点，则元素记录像素点的节点编号，否则记为o。使用这个数组可以通过像素点的坐标值来查询它的节点编号。II.边线起始点数组用来存放在裂纹骨架树中从一个节点到另一个节点连线上第一个枝干点的坐标。它的行号对应分支节点号，列号对应另一节点号，纵向编号对应像素的x坐标和y坐标。比如数组下标[3][4][2]表示骨架上从3号分支节点到4号节点连线上第一个枝干点的y坐标。使用这个数组可以査找到一条边线的起始点，只要将该起始点设为背景点，就可以将边线和裂纹骨架断开。算法步骤如下(1)对每个未结束分支节点进行处理，若该节点属于I类分支节点，则保留其最长的外边线，将其它短的外边线在裂纹骨架中断开，修改节点信息数组；(2)对每个未结束分支节点进行处理，若该节点属于II类分支节点，将该节点退化为枝干点，与该节点相连的分支节点成为新的末端分支节点，即m类分支节点退化为n类分支节点，新的外边线长度为保留的最长外边线与内边线长度之和，结束该分支节点；G)重复(1)(2)的处理直至剩下一个IV类分支节点；(4)比较该IV类分支节点的外边线长度，最长的两条构成裂纹主干，断开其它外边线。(5)去掉多余前景点，将裂纹主干上被标记为分支节点的前景点置为背景点，若裂纹主干连通性不变，则继续，若改变，则该节点保留。表l裂纹骨架像素点规定<table>tableseeoriginaldocumentpage10</column></row><table>表2裂纹骨架数的边线规定<table>tableseeoriginaldocumentpage10</column></row><table>表3裂纹骨架树的节点规定<table>tableseeoriginaldocumentpage10</column></row><table>本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。权利要求1、一种CCD图像特征的表面疲劳裂纹检测方法，其特征是第一步对试样进行疲劳裂纹扩展试验；对试样进行疲劳裂纹扩展试验，利用图像采集设备连续或间隙采集其表面疲劳裂纹产生过程中的数字图像，将获得的数字图像输入计算机；第二步对表面疲劳裂纹图像进行预处理；将试样的疲劳裂纹分为微小裂纹和长裂纹，利用MATLAB图像处理技术，分别对微小裂纹和长裂纹图像进行预处理，通过裁截、缩放操作获得微小裂纹区域图像，通过灰度化、二值化、膨胀腐蚀、开闭操作、细化操作获得长裂纹骨架图像；第三步获取试样表面疲劳裂纹的长度特征和宽度特征；(1)利用鼠标点击微小裂纹区域，获得点的坐标值，进而计算出微小裂纹的像素长度，像素长度单位为pixel；利用树结构存储长裂纹骨架信息，采用递归的方法去除骨架分支，提取出长裂纹主干，进而计算出长裂纹的像素长度，像素长度单位为pixel；(2)在所获取的试样的疲劳裂纹扩展过程的全部图像中，以表面疲劳裂纹的扩展初始位置为研究对象，利用鼠标点击获得点的坐标的方法计算表面疲劳裂纹的像素宽度，像素宽度单位为pixel；(3)对图像采集设备进行标定，计算放大因子(单位mm/pixel)，进而计算表面疲劳裂纹的实际长度(单位mm)和实际宽度(单位mm)；实际长度(mm)＝像素长度(pixel)×放大因子(mm/pixel)，实际宽度(mm)＝像素长度(pixel)×放大因子(mm/pixel)；(4)利用指数曲线拟合方法，拟合得到表面疲劳裂纹扩展的长度曲线(a-N曲线)、宽度曲线(w-N曲线)、长度扩展速率曲线(da/dN-N曲线)和宽度扩展速率曲线(dw/dN-N曲线)。全文摘要本发明针对现有的疲劳裂纹检测手段精度不高，需要接触试样，设备昂贵，实施繁琐等缺陷，公开了一种CCD图像特征的表面疲劳裂纹检测方法，它利用MATLAB图像处理技术，研究了表征表面疲劳裂纹扩展的长度与宽度特征参数，比较了两种表征方法，并验证了该表面疲劳裂纹图像检测方法的精度。该方法具有高精度、非接触式、设备简单、易实现等优点，通过该方法可以实现故障诊断的智能化，为生产实际提供指导性建议。文档编号G01N21/88GK101413901SQ200810235410公开日2009年4月22日申请日期2008年12月1日优先权日2008年12月1日发明者丁子昀,左敦稳,莫国影,陈兴峰,黎向锋申请人:南京航空航天大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
