主权项 |
一种基于高温图像处理的裂纹尖端张开位移测量方法,其特征在于该方法包括如下步骤:1)将带有预制裂纹的试件放入高温恒温环境下进行力学加载,该试件在高温环境中裂纹尖端发生局部氧化,裂纹尖端图像局部灰度值阶跃变化,使用摄像机连续拍摄加载过程中试件裂纹尖端的变化过程,获得同一温度下试件缺口裂纹在不同时刻的扩展图像,并在摄像机镜头上安装至少一个高温滤波片;2)采用二阶导数来检测、计算裂纹尖端边缘,如下式计算:定义二维图像灰度函数f(x,y)的拉普拉斯值为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mo>▿</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>∂</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>∂</mo><mi>x</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</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>∂</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>∂</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>∂</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>∂</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>∂</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>▿</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>∂</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>∂</mo><mi>x</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</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>∂</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>)。 |