专利名称：一种基于光学剪切的二维光学应变花测量方法
技术领域：
本发明涉及一种基于光学剪切的二维光学应变花测量方法。
背景技术：
随着材料科学技术的发展，精确的测量材料的力学性能成为了新型材料应用的必要条件， 如何无损的定量测量材料的变形场成为了研究的重点。由于光学测量方法的非接触无损测量 等特性，各种光学方法已经被广泛应用于材料变形场的测量。
光学方法中主要有干涉方法和非干涉方法。其中干涉方法包括全息千涉法、云纹干涉
法、散斑干涉法等，非干涉方法主要是数字相关方法，这些方法中，干涉方法主要是利用激 光干涉技术，通过分析变形引起的干涉条纹得到试件表面的变形场信息，而数字相关方法主 要是通过对变形前后的图像进行相关处理分析变形场信息。对于干涉方法，由于其采用千涉 原理，因此需要在隔振台上进行实验，受周围环境影响比较大，无法进行复杂工况下的现场 测量。而数字相关方法虽然受环境影响小，但是此方法位移求解精度同图像放大倍数相关， 随着图像放大倍数提高，位移求解精度提高，但是相应的应变计算的区域会减小。
本发明人曾在2007年4月申请了名为一种基于错位相关原理的物体应变测量方法及装置 的发明专利，在该专利中，利用一块楔形镜改变物体成像位置从而改变应变计算标距，提高 了应变的计算精度。但是该方法主要解决的是光学方法中应变求解精度问题，且只能进行一 点的应变计算，而我们对材料进行分析的时候常常需要知道某个区域的力学性能，因此需要 发展一种大区域应变测量方法
上述各种方法，在实际应用中都存在着一定的问题，无法实现操作简单方便，在复杂工 作环境下的现场实时的高精度大区域的无损测量。
发明内容
本发明的目的是提供一种基于光学剪切的二维光学应变花测量方法，可以实现操作简单 方便，在复杂工作环境下进行现场实时的高精度大区域的无损测量。
本发明所采用的技术方案如下 一种基于光学剪切的二维光学应变花测量方法，其特征在 于该方法按如下步骤进行
1) 首先在被测物体表面制作正方形光学应变花标记(a，b,c,d)，边长为2/;然后在成像 光路中放置一块四棱镜，引入光学剪切量A，使物体所成的像变为四个完全相同且中心对称
的像，分别选取该四个像中的正方形光学应变花标记的一个点，组成新的正方形光学应变花
标记，最后使用数字成像设备分别记录变形前后该新的光学应变花标记像；
2) 在所纪录的图像中，沿着新的正方形光学应变花标记点建立直角坐标系，得到变形前 该四个标记点的坐标分别为A点(x^,h)， B点C点(xc,jc)， D点(xD,j;D)，通过数字相关方法可以得到变形后该四个标记点的坐标分别为A点(《，K)， B点(《，乂)， C点 (x:,《)，D点(x:,;;:);
3) 此四个标记点的x，y方向的位移分别为A点(^=《-力)，B点
(w5 =x:_xB,v￡ >"B )， C,点(wc =x;—xc,vc >>c )， D,点(Wd =x:—xD,vD );
其中^,MB,"C,MD为X方向^立移，^,V;j,Vc,Vd为y方向位移；
4) 根据有限单元方法中四节点矩形单元的应变计算公式即可得到正方形光学应变花区域 中任意一点(x，y)处的应变分布
x方向应变为^ =^r[^(/ + >0 —&(/ + ;;)-&(/-力+ (/一;/)] y方向应变为~ = +(/+力+ vB (/— x) - vc (/— x)— (/+ x)] xy夹角方向应变为
、^"^T[ ("x) + ^G +力+ "s(Z-工)—vb(/ + >0 —wc" —x) —VcG —力一"z)G + x) + Vd(7-力]。
上述测量方法中光学剪切量A同四棱镜的棱角"成正比，且同四棱镜距成像透镜的距
离4成正比。
上述测量方法中所用正方形光学应变花标记采用热压印方法制作。
本发明同现有方法和装置相比，具有如下优点及突出性效果采用热压印方法制作正方 形光学应变花标记，制作方便省吋省力；在测量光路中放置四棱镜引入光学剪切量，形成具 有剪切量的图像，增大应变求解区域；根据数字相关方法计算具有剪切量的图像中光学应变 花的四个标记点的位移，根据有限兀方法计算应变花区域的应变，可以实现操作简单方便， 在复杂工作环境下进行现场实时的高精度大区域的无损测量。
图1在被测物体表面制作的光学应变花标记的示意图。 图2无四棱镜成像光路示意图。 图3有四棱镜成像光路示意图。 图4四棱镜立体结构示意图。 图5位移计算示意图。 图6应变花标记点对应关系示意图 l-成像透镜；2-四棱镜；
具体实施例方式
下面结合附图对本发明的原理和具体实施方法作进一步的说明。
图1为本发明使用的正方形光学应变花标记示意图，本发明采用热压印工艺在试件表面 制作标记，该方法首先使用机械加工或者光学刻蚀等方法制作出应变花标记模板，然后将试件清洗干净，烘干后均匀涂抹光刻胶，接着将涂有光刻胶的试件放置在热压印仪上，将加工
好的应变花标记模板放置在试件需要制作标记的位置，最后用热压印仪升温加压即可在试件
表面复制出应变花标记(a》，c,d)。
图2为无四棱镜成像系统的光路示意图，试件表面的正方形光学应变花标记(",6，c,d)
经过成像透镜1，成像在成像面上，原始标记间距为2/，成像面的正方形光学应变花标记的
像(a)'，c',cT )的间距为若假设采用的数字成像设备的靶面大小为f，则成像要求2&<"
成像系统放大倍数为^-〃t。对于有四棱镜的成像系统，其成像光路如图3所示，试件表
面的正方形光学应变花标记(",6,c^)首先经过四棱镜，引入光学剪切，然后再通过成像透
镜成像在成像面匕由于四棱镜引入的光学剪切，使得正方形光学应变花在成像面上的像变
为了具有剪切量的4个像，而本发明所采用的四棱镜如图4所示，为对称结构的四棱镜，所
以成像面上的像为完全对称的剪切量相同的4个像。如图5示，由四棱镜产生具有剪切量的
像后，分别从4个对称的像中选择1个标记点，这4个新的标记点A、 B、 C、 D组成一个新的
正方形应变花标记，该正方形的边长等于剪切量A，其值同四棱镜的棱角"以及四棱镜距成
像透镜的距离/6成正比。因此通过更换不同角度的四棱镜或者调整四棱镜和成像透镜之间的
距离就可以改变剪切量的大小。同样假设用的数字成像设备的耙面大小为"则成像要求A < "
成像系统放大倍数为附2=2〃厶，由于A的值可以通过棱角"和距离4调整，因此可以使
A<2L即对于有四棱镜的成像系统，我们所能采用的成像系统的放大倍数^<^，采用数
字相关方法计算应变花标记点的位移时，其精度也会比无四棱镜的成像系统高；另一方面，
若两套系统的放大倍数相同，由于A〈2"则有四棱镜的成像系统的正方形光学应变花的间
距将会大于无四棱镜的成像系统的正方形光学应变花的间距，即有四棱镜的成像系统所能测
量的应变的区域将会大于无四棱镜的成像系统所能测量的应变的区域。
这个新的应变花标记的4个标记点和原始标记的4个标记点之间存在一个对应关系(图6)，
即A点对应a点，B点对应b点，C点对应c点，D点对应d点。沿着新的正方形应变花标记的边建
立直角坐标系，则可以得到变形前A、 B、 C、 D这4个点的坐标值(A,h)， (xs，_ys)， (xc，;ve)，
采用数字相关方法计算变形后A、 B、 C、 D这4个点的坐标值(x二,K)， (x:,K)，
(xL,)。此四个标记点的x,y方向的位移分别为A点^ =《-=乂-力，B
点^ =《-xs,vfl &， C点Wc二《-xc，vc=;/:_;(;c， D点 =《--。根据
有限单元方法中四节点矩形单元的应变计算公式即可得到光学应变花区域中任意一点(x, y) 处的应变分布
& =h (,") - "s(7 ") _ wc G _力+wd (7 _ >0]
~ =+[、 G + x)十vb (z—x)—vc <7—x)— "+力]
=;[^i(7+x)+^G+_y)+"s(7-^)—vfl(7+>0—"c(/—x)-vc(/—力-wD(/ + x)+vD(/—力]
5有上述分析可以知道，在有四棱镜的成像系统中通过调整棱角"和四棱镜距成像透镜的距离 /6，可以实现高放大倍数大区域的应变测量。
本发明的具体测量过程为，首先在试件表面需要测试区域，用热压印方法制作光学应变
花标记，然后将试件放入带有四棱镜的光学成像光路，由于四棱镜产生光学剪切，形成新的 应变花标记区域，接着使用数字成像设备分别记录变形前后该区域的图像，在变形前的区域 中建立坐标系得到变形前应变花标记点的坐标，利用数字图像相关方法计算变形后光学应变 花各个标记的坐标，根据变形前后的坐标可以得到应变花标记点的位移值，最后根据有限单 元方法中四节点矩形争元的应变计算公式计算出该区域的应变值。
权利要求
1.一种基于光学剪切的二维光学应变花测量方法，其特征在于该测量方法按如下步骤进行1)首先在被测物体表面制作正方形光学应变花标记，边长为2l；然后在成像光路中放置一块四棱镜，引入光学剪切量Δ，使物体所成的像变为四个完全相同且中心对称的像，分别选取四个像中的正方形光学应变花标记中的一个点，组成新的正方形光学应变花标记，最后使用数字成像设备分别记录变形前后所述新的正方形光学应变花标记像；2)在所纪录的图像中，沿着新的正方形光学应变花标记点建立直角坐标系，得到变形前四个标记点的坐标分别为A点(xA，yA)，B点(xB，yB)，C点(xC，yC)，D点(xD，yD)，通过数字相关方法可以得到变形后所述四个标记点的坐标分别为A点(xA*，yA*)，B点(xB*，yB*)，C点(xC*，yC*)，D点(xD*，yD*)；3)新的正方形光学应变花标记的四个标记点的x，y方向的位移分别为A点<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><mrow> <mo>(</mo> <msub><mi>u</mi><mi>A</mi> </msub> <mo>=</mo> <msubsup><mi>x</mi><mi>A</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>x</mi><mi>A</mi> </msub> <mo>,</mo> <msub><mi>v</mi><mi>A</mi> </msub> <mo>=</mo> <msubsup><mi>y</mi><mi>A</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>y</mi><mi>A</mi> </msub> <mo>)</mo></mrow><mo>,</mo> </mrow>]]></math> id="icf0001" file="A2009101319220002C1.tif" wi="46" he="4" top= "115" left = "24" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>B点<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><mrow> <mo>(</mo> <msub><mi>u</mi><mi>B</mi> </msub> <mo>=</mo> <msubsup><mi>x</mi><mi>B</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>x</mi><mi>B</mi> </msub> <mo>,</mo> <msub><mi>v</mi><mi>B</mi> </msub> <mo>=</mo> <msubsup><mi>y</mi><mi>B</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>y</mi><mi>B</mi> </msub> <mo>)</mo></mrow><mo>,</mo> </mrow>]]></math> id="icf0002" file="A2009101319220002C2.tif" wi="46" he="4" top= "115" left = "81" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>C点<maths id="math0003" num="0003" ><math><![CDATA[ <mrow><mrow> <mo>(</mo> <msub><mi>u</mi><mi>C</mi> </msub> <mo>=</mo> <msubsup><mi>x</mi><mi>C</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>x</mi><mi>C</mi> </msub> <mo>,</mo> <msub><mi>v</mi><mi>C</mi> </msub> <mo>=</mo> <msubsup><mi>y</mi><mi>C</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>y</mi><mi>C</mi> </msub> <mo>)</mo></mrow><mo>,</mo> </mrow>]]></math> id="icf0003" file="A2009101319220002C3.tif" wi="47" he="4" top= "115" left = "137" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>D点<maths id="math0004" num="0004" ><math><![CDATA[ <mrow><mrow> <mo>(</mo> <msub><mi>u</mi><mi>D</mi> </msub> <mo>=</mo> <msubsup><mi>x</mi><mi>D</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>x</mi><mi>D</mi> </msub> <mo>,</mo> <msub><mi>v</mi><mi>D</mi> </msub> <mo>=</mo> <msubsup><mi>y</mi><mi>D</mi><mo>*</mo> </msubsup> <mo>-</mo> <msub><mi>y</mi><mi>D</mi> </msub> <mo>)</mo></mrow><mo>;</mo> </mrow>]]></math> id="icf0004" file="A2009101319220002C4.tif" wi="48" he="4" top= "126" left = "31" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>其中uA，uB，uC，uD为x方向位移，vA，vB，vC，vD为y方向位移；4)根据有限单元方法中四节点矩形单元的应变计算公式，即得到原正方形光学应变花区域中任意一点(x，y)处的应变分布x方向应变为<maths id="math0005" num="0005" ><math><![CDATA[ <mrow><msub> <mi>ϵ</mi> <mi>x</mi></msub><mo>=</mo><mfrac> <mn>1</mn> <msup><mrow> <mn>4</mn> <mi>l</mi></mrow><mn>2</mn> </msup></mfrac><mo>[</mo><msub> <mi>u</mi> <mi>A</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>y</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>u</mi> <mi>B</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>y</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>u</mi> <mi>C</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>u</mi> <mi>D</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo></mrow><mo>]</mo> </mrow>]]></math> id="icf0005" file="A2009101319220002C5.tif" wi="101" he="8" top= "145" left = "51" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>y方向应变为<maths id="math0006" num="0006" ><math><![CDATA[ <mrow><msub> <mi>ϵ</mi> <mi>y</mi></msub><mo>=</mo><mfrac> <mn>1</mn> <msup><mrow> <mn>4</mn> <mi>l</mi></mrow><mn>2</mn> </msup></mfrac><mo>[</mo><msub> <mi>v</mi> <mi>A</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>x</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>v</mi> <mi>B</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>v</mi> <mi>C</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>v</mi> <mi>D</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>x</mi> <mo>)</mo></mrow><mo>]</mo> </mrow>]]></math> id="icf0006" file="A2009101319220002C6.tif" wi="81" he="8" top= "155" left = "51" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>xy夹角方向应变为<maths id="math0007" num="0007" ><math><![CDATA[ <mrow><msub> <mi>ϵ</mi> <mi>xy</mi></msub><mo>=</mo><mfrac> <mn>1</mn> <msup><mrow> <mn>4</mn> <mi>l</mi></mrow><mn>2</mn> </msup></mfrac><mo>[</mo><msub> <mi>u</mi> <mi>A</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>x</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>v</mi> <mi>A</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>y</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>u</mi> <mi>B</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>v</mi> <mi>B</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>y</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>u</mi> <mi>C</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>v</mi> <mi>C</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>u</mi> <mi>D</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>+</mo> <mi>x</mi> <mo>)</mo></mrow><mo>+</mo><msub> <mi>v</mi> <mi>D</mi></msub><mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo></mrow><mo>]</mo><mo>.</mo> </mrow>]]></math></maths>
2. 上如权利要求l所述的基于光学剪切的二维光学应变花测量方法，其特征在于光学剪切 量A同四棱镜的棱角《成正比，且同四棱镜距成像透镜的距离/6成正比。
3. 如权利要求l所述的基于光学剪切的二维光学应变花测量方法，其特征在于所用正方形 光学应变花标记采用热压印方法制作。
全文摘要
一种基于光学剪切的二维光学应变花测量方法，该测量方法首先采用热压印方法在被测物体表面制作矩形光学应变花标记，在测量光路中放置四棱镜，引入光学剪切量，形成具有剪切量的图像，根据数字相关方法计算具有剪切量的图像中光学应变花的4个标记点的位移，根据有限元方法计算应变花区域的应变。本发明有效提高了图像质量，改善了数字相关方法的应变求解区域不大的问题，可以实现操作简单方便，在复杂工作环境下进行现场实时的高精度大区域的无损测量。
文档编号G01B11/16GK101514890SQ20091013192
公开日2009年8月26日 申请日期2009年3月27日 优先权日2009年3月27日
发明者戴福隆, 韬 花, 谢惠民 申请人:清华大学