[发明专利]一种基于高温图像处理的裂纹尖端张开位移测量方法有效

专利信息
申请号: 201410539452.7 申请日: 2014-10-13
公开(公告)号: CN104266900A 公开(公告)日: 2015-01-07
发明(设计)人: 冯雪;方旭飞;张长兴 申请(专利权)人: 清华大学
主分类号: G01N3/00 分类号: G01N3/00;G01N3/06;G01B11/02
代理公司: 北京鸿元知识产权代理有限公司 11327 代理人: 邸更岩
地址: 100084 北京市海淀区1*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要: 一种基于高温图像处理的裂纹尖端张开位移测量方法,属于工程材料、结构形变及力学实验技术领域。本发明的技术特点是通过对带预制裂纹的试件进行高温力学加载实验,通过配备滤波装置的摄像机拍摄试件表面的裂纹的扩展过程,以及通过对获取的图像采用边缘检测的方法对试件缺口裂纹进行定位检测并计算裂纹尖端的张开位移,从而获得材料在高温环境下的裂纹尖端张开位移(crack-tip opening displacement,CTOD)。该方法通过图像处理技术直接计算裂纹尖端张开位移,与常规的通过应力-位移曲线计算得到材料裂纹尖端张开位移的方法相比,该方法简单直接,可通过软件辅助计算直接得到结果,减少人为判断常规的应力-位移曲线中的误差。
搜索关键词: 一种 基于 高温 图像 处理 裂纹 尖端 张开 位移 测量方法
【主权项】:
一种基于高温图像处理的裂纹尖端张开位移测量方法,其特征在于该方法包括如下步骤:1)将带有预制裂纹的试件放入高温恒温环境下进行力学加载,该试件在高温环境中裂纹尖端发生局部氧化,裂纹尖端图像局部灰度值阶跃变化,使用摄像机连续拍摄加载过程中试件裂纹尖端的变化过程,获得同一温度下试件缺口裂纹在不同时刻的扩展图像,并在摄像机镜头上安装至少一个高温滤波片;2)采用二阶导数来检测、计算裂纹尖端边缘,如下式计算:定义二维图像灰度函数f(x,y)的拉普拉斯值为:<mrow><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>&PartialD;</mo><mi>y</mi></mrow><mn>2</mn></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>其中获取图像的面内坐标为(x,y),x,y对应图像像素点的横坐标和纵坐标,上式中在点(x,y)处的的差分表示分别为:<mrow><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>2</mn></msup></mfrac><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>2</mn></msup></mfrac><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>将公式(2)和(3)代入公式(1),得到下式(4):<mrow><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>&PartialD;</mo><mi>y</mi></mrow><mn>2</mn></msup></mfrac><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mn>4</mn><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow>根据公式(4)得到对应的二维图像灰度函数f(x,y)的拉普拉斯值的局部最大值;3)根据裂纹尖端边缘具体位置的扩展变化,计算材料裂纹尖端张开位移:定义函数g(xmax,ymax)表示检测到裂纹尖端边缘对应不同时刻的位置,xmax,ymax是对应图像不同时刻通过边缘检测的方法获得二维图像灰度函数f(x,y)的拉普拉斯值的局部最大值的横坐标和纵坐标;时刻t0对应的裂纹尖端位置为g0(xmax,ymax),时刻t1对应的裂纹尖端位置为g1(xmax,ymax);记实验进行时间Δt=t1‑t0,得到裂纹尖端张开位移δ=g1(xmax,ymax)‑g0(xmax,ymax)。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

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

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

×

专利文献下载

说明:

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

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

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

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

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

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

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

钻瓜专利网在线咨询

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

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