发明名称 一种基于高温图像处理的裂纹尖端张开位移测量方法
摘要 一种基于高温图像处理的裂纹尖端张开位移测量方法,属于工程材料、结构形变及力学实验技术领域。本发明的技术特点是通过对带预制裂纹的试件进行高温力学加载实验,通过配备滤波装置的摄像机拍摄试件表面的裂纹的扩展过程,以及通过对获取的图像采用边缘检测的方法对试件缺口裂纹进行定位检测并计算裂纹尖端的张开位移,从而获得材料在高温环境下的裂纹尖端张开位移(crack-tip opening displacement,CTOD)。该方法通过图像处理技术直接计算裂纹尖端张开位移,与常规的通过应力-位移曲线计算得到材料裂纹尖端张开位移的方法相比,该方法简单直接,可通过软件辅助计算直接得到结果,减少人为判断常规的应力-位移曲线中的误差。
申请公布号 CN104266900A 申请公布日期 2015.01.07
申请号 CN201410539452.7 申请日期 2014.10.13
申请人 清华大学 发明人 冯雪;方旭飞;张长兴
分类号 G01N3/00(2006.01)I;G01N3/06(2006.01)I;G01B11/02(2006.01)I 主分类号 G01N3/00(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 一种基于高温图像处理的裂纹尖端张开位移测量方法,其特征在于该方法包括如下步骤:1)将带有预制裂纹的试件放入高温恒温环境下进行力学加载,该试件在高温环境中裂纹尖端发生局部氧化,裂纹尖端图像局部灰度值阶跃变化,使用摄像机连续拍摄加载过程中试件裂纹尖端的变化过程,获得同一温度下试件缺口裂纹在不同时刻的扩展图像,并在摄像机镜头上安装至少一个高温滤波片;2)采用二阶导数来检测、计算裂纹尖端边缘,如下式计算:定义二维图像灰度函数f(x,y)的拉普拉斯值为:<maths num="0001" id="cmaths0001"><math><![CDATA[<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>]]></math><img file="FDA0000585842930000017.GIF" wi="1330" he="153" /></maths>其中<img file="FDA0000585842930000011.GIF" wi="295" he="148" />获取图像的面内坐标为(x,y),x,y对应图像像素点的横坐标和纵坐标,上式中在点(x,y)处的<img file="FDA0000585842930000012.GIF" wi="236" he="136" /><img file="FDA0000585842930000013.GIF" wi="213" he="147" />的差分表示分别为:<maths num="0002" id="cmaths0002"><math><![CDATA[<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>]]></math><img file="FDA0000585842930000014.GIF" wi="1409" he="140" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<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>]]></math><img file="FDA0000585842930000015.GIF" wi="1410" he="142" /></maths>将公式(2)和(3)代入公式(1),得到下式(4):<maths num="0004" id="cmaths0004"><math><![CDATA[<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>]]></math><img file="FDA0000585842930000016.GIF" wi="1907" he="152" /></maths>根据公式(4)得到对应的二维图像灰度函数f(x,y)的拉普拉斯值的局部最大值;3)根据裂纹尖端边缘具体位置的扩展变化,计算材料裂纹尖端张开位移:定义函数g(x<sub>max</sub>,y<sub>max</sub>)表示检测到裂纹尖端边缘对应不同时刻的位置,x<sub>max</sub>,y<sub>max</sub>是对应图像不同时刻通过边缘检测的方法获得二维图像灰度函数f(x,y)的拉普拉斯值的局部最大值的横坐标和纵坐标;时刻t<sub>0</sub>对应的裂纹尖端位置为g<sub>0</sub>(x<sub>max</sub>,y<sub>max</sub>),时刻t<sub>1</sub>对应的裂纹尖端位置为g<sub>1</sub>(x<sub>max</sub>,y<sub>max</sub>);记实验进行时间Δt=t<sub>1</sub>‑t<sub>0</sub>,得到裂纹尖端张开位移δ=g<sub>1</sub>(x<sub>max</sub>,y<sub>max</sub>)‑g<sub>0</sub>(x<sub>max</sub>,y<sub>max</sub>)。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室