CN1320340C - 一种提取相干梯度敏感干涉条纹级数的方法及其装置 - Google Patents

一种提取相干梯度敏感干涉条纹级数的方法及其装置 Download PDF

Info

Publication number
CN1320340C
CN1320340C CNB2006100013666A CN200610001366A CN1320340C CN 1320340 C CN1320340 C CN 1320340C CN B2006100013666 A CNB2006100013666 A CN B2006100013666A CN 200610001366 A CN200610001366 A CN 200610001366A CN 1320340 C CN1320340 C CN 1320340C
Authority
CN
China
Prior art keywords
interference fringe
grating
light
coherent gradient
gradient sensing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CNB2006100013666A
Other languages
English (en)
Other versions
CN1800835A (zh
Inventor
姚学锋
许蔚
刘栋梁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CNB2006100013666A priority Critical patent/CN1320340C/zh
Publication of CN1800835A publication Critical patent/CN1800835A/zh
Application granted granted Critical
Publication of CN1320340C publication Critical patent/CN1320340C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明涉及一种提取相干梯度敏感干涉条纹级数的方法及其装置,属于断裂力学参量测量的技术领域。首先将加载试件出射的平行光束分为两路光线,分别产生相干梯度敏感干涉条纹;在两幅相干梯度敏感干涉条纹图中减去相应背景光强,并对其进行直方图变换,使两幅图像的条纹光强的幅度相同,根据光栅间距和光强,计算第一幅干涉条纹图的条纹级数m1。本发明的装置包括激光器、扩束镜、凸透镜、半反半透镜、四个光栅、两组凸透镜和光阑、两个摄像机和计算机。本发明的优点是能够明显提高相干梯度敏感方法的计算精度,不要求预先知道每幅干涉条纹图像对应的载荷,因此便于在动态断裂实验中应用。

Description

一种提取相干梯度敏感干涉条纹级数的方法及其装置
技术领域
本发明涉及一种提取相干梯度敏感干涉条纹级数的方法及其装置,属于断裂力学参量测量的技术领域。
背景技术
相干梯度敏感方法(Coherent Gradient Sensing,以下简称CGS)是一种具有在线空间滤波效果的非接触、全场横向双栅剪切干涉光测方法,它利用两个平行光栅重组由试件变形导致的扭曲光束从而产生干涉条纹,通过光学干涉控制方程建立面内应力梯度或离面位移梯度与条纹级数之间的关系。该方法能够在透射和反射情况下分别用于研究透明和非透明物体的变形信息,其条纹分布直接反映了透射情况下的面内应力梯度场和反射情况下的离面位移梯度场。由于该方法对应力或位移梯度敏感,所以尤其适合于研究具有高应力集中现象的断裂问题。相干梯度敏感方法分为透射式和反射式两种,图1所示为典型相干梯度敏感干涉条纹图,实际拍摄图像中黑条纹代表半级数条纹,最外圈为0.5级条纹。
上述已有技术的不足之处:与其它处理条纹图的光测方法类似,要获得断裂力学参量必须先提取相干梯度敏感干涉条纹级数,而要确定相干梯度敏感干涉条纹级数就必须精确提取条纹中心线。从图1可以看出干涉条纹较粗,并且形状不规则,如果仅仅依靠图像处理的方法提取条纹中心线难以满足精度要求,甚至会带来很大的误差,尤其是对于低级数条纹。
发明内容
本发明的目的是提出一种提取相干梯度敏感干涉条纹级数的方法及其装置,以提高相干梯度敏感方法的计算精度。
本发明提出的提取相干梯度敏感干涉条纹级数的方法,包括以下步骤:
(1)对一束激光进行扩束,得到发散光,再经凸透镜成为平行光,入射到试件表面;
(2)将上述试件出射的平行光束分为两路光线,第一路光线经过间距为Δ1的第一光栅和第二光栅,产生相干梯度敏感干涉条纹,第二路光线经过间距为Δ2的第三光栅和第四光栅,产生相干梯度敏感干涉条纹;
(3)对上述两个干涉条纹场成像,并拍摄干涉条纹场;
(4)在试件未加载前,分别从上述两个光路拍摄试件的数字图像,获得干涉条纹图的背景光强分别为a1、a2
(5)对试件施加载荷,分别从上述两个光路拍摄加载试件的相干梯度敏感干涉条纹图,两个干涉条纹场的光强分布分别为:I1、I2
(6)从上述步骤(5)的加载试件的两幅相干梯度敏感干涉条纹图中减去上述步骤(4)的相应背景光强,并对其进行直方图变换,使两幅图像条纹光强的幅度相同,变换获得的两个光路的干涉条纹图光强分布为I1 *和I2 *
(7)根据上述光栅间距和光强,计算第一幅干涉条纹图的条纹级数m1为:
I 1 * I 2 * = cos ( 2 πm 1 ) cos ( 2 πm 1 Δ 2 / Δ 1 ) .
本发明提出的提取相干梯度敏感干涉条纹级数的装置,包括:
(1)用于产生激光的激光器,用于对激光进行扩束的扩束镜,使发散光成为平行光的凸透镜;
(2)用于将试件出射的平行光束分为两路光线的半反半透镜,用于产生第一路相干梯度敏感干涉条纹的第一光栅和第二光栅,第一光栅和第二光栅的间距为Δ1,用于产生第二路相干梯度敏感干涉条纹的第三光栅和第四光栅,第三光栅和第四光栅的间距为Δ2
(3)用于对上述两个干涉条纹场成像的第一组凸透镜和光阑以及第二组凸透镜和光阑,以及用于拍摄干涉条纹场的第一摄像机、第二摄像机和用于图像处理和干涉条纹级数提取的计算机;
上述激光器、扩束镜、凸透镜、试件、半反半透镜、第一光栅、第二光栅、第一组凸透镜和光阑和第一摄像机依次处于第一光轴上;上述第三光栅、第四光栅、第二组凸透镜和光阑和第二摄像机依次处于第二光轴上,第一光轴与第二光轴于半反半透镜处垂直相交,计算机分别与第一摄像机和第二摄像机相连接。
本发明提出的提取相干梯度敏感干涉条纹级数的方法及其装置,其优点是:可以精确提取相干梯度敏感干涉条纹图任意位置点的条纹级数,实验结果表明该方法能够显著减少由于条纹提取带来的误差,从而提高相干梯度敏感方法的精度。具体优点是:a)能够明显提高相干梯度敏感方法的计算精度;b)不要求预先知道每幅干涉条纹图像对应的载荷,因此便于在动态断裂实验中应用。
附图说明
图1是已有的典型相干梯度敏感干涉条纹图,其中(a)是理论模拟图;(b)是实际拍摄条纹图。
图2是已有技术中产生干涉条纹的光路图。
图3是已有技术中产生干涉条纹的原理图。
图4是本发明提出的提取相干梯度敏感干涉条纹级数的光路示意图。
图5是采用本发明方法,在同一试件上不同光栅间距时产生的裂纹尖端反射相干梯度敏感干涉条纹图,其中(a)是光栅间距为Δ1=50mm,(b)是光栅间距为Δ2=100mm。
图6是本发明方法与已有技术计算得到的应力强度因子结果对比图,纵轴是实验计算值与理论值的相对误差,实线代表已有技术的计算结果,虚线代表本发明方法的计算结果。其中(a)为在干涉条纹图-40°~-30°之间的计算结果对比图,(b)为在干涉条纹图30°~40°之间的计算结果对比图。
图2、图3和图4中,1是激光器,2是扩束镜,3是凸透镜,4是试件,5是第一光栅,6是第二光栅,7是第一滤波透镜,8是第一光阑,9是第一摄像机,10是计算机,11是半反半透镜,12是第三光栅,13是第四光栅,14是第二滤波透镜,15是第二光阑,16是第二摄像机。
具体实施方式
本发明提出的提取相干梯度敏感干涉条纹级数的方法,首先对一束激光进行扩束,得到发散光,再经凸透镜成为平行光,入射到试件表面;将试件出射的平行光束分为两路光线,第一路光线经过间距为Δ1的第一光栅和第二光栅,产生相干梯度敏感干涉条纹,第二路光线经过间距为Δ2的第三光栅和第四光栅,产生相干梯度敏感干涉条纹;对两个干涉条纹场成像,并拍摄干涉条纹场;在试件未加载前,分别从两个光路拍摄试件的数字图像,获得干涉条纹图的背景光强分别为a1、a2;对该试件施加载荷,分别从两个光路拍摄加载试件的相干梯度敏感干涉条纹图,两个干涉条纹场的光强分布分别为:I1、I2;从加载试件的两幅相干梯度敏感干涉条纹图中减去相应背景光强,并对其进行直方图变换,使两幅图像的条纹光强的幅度相同,变换获得的两个光路的干涉条纹图光强分布为I1 *和I2 *;根据光栅间距和光强,计算第一幅干涉条纹图的条纹级数m1为: I 1 * I 2 * = cos ( 2 πm 1 ) cos ( 2 π m 1 Δ 2 / Δ 1 ) .
本发明提出的提取相干梯度敏感干涉条纹级数的装置,如图4所示,包括:
(1)用于产生激光的激光器1,用于对激光进行扩束的扩束镜2,使发散光成为平行光的凸透镜3;
(2)用于将试件4出射的平行光束分为两路光线的半反半透镜11,用于产生第一路相干梯度敏感干涉条纹的第一光栅5和第二光栅6,第一光栅5和第二光栅6的间距为Δ1,用于产生第二路相干梯度敏感干涉条纹的第三光栅12和第四光栅13,第三光棚12和第四光栅13的间距为Δ2,两个干涉条纹场的光强分布分别为:I1、I2
(3)用于对上述两个干涉条纹场成像的第一组凸透镜7和光阑8以及第二组凸透镜14和光阑15,以及用于拍摄干涉条纹场的第一摄像机9、第二摄像机16和用于图像处理和干涉条纹级数提取的计算机10;
上述激光器1、扩束镜2、凸透镜3、试件4、半反半透镜11、第一光栅5、第二光栅6、第一组凸透镜7和光阑8、第一摄像机9依次处于第一光轴上;上述第三光栅12、第四光栅13、第二组凸透镜14和光阑15、第二摄像机16依次处于第二光轴上,第一光轴与第二光轴于半反半透镜11处垂直相交,计算机10分别与第一摄像机9和第二摄像机16相连接。
以下结合附图,详细介绍本发明的内容:
已有的相干梯度敏感方法光路如图2所示:一束相干激光束经过带有滤波孔的扩束镜后通过光学透镜形成平行光,平行光束垂直通过带裂纹的试件,受到试件变形和应力场的影响,平行光束的方向和相位将产生变化。由试件出射的光束虽然不再是平行光束,但是由于方向偏离很小,为简化起见,仍然将带有试件变形信息的光束视为平行光束。由试件上(x,y)点出射的光线由于试件变形具有光程差δS(x,y),类似的由试件上(x,y+ε)点出射的光线由于试件变形具有光程差δS(x,y+ε)。如图3所示,这些平行光束从试件出射后经过两个互相平行、间距为Δ的高密度Ronchi光栅(节距相同为p),产生衍射,如果在光阑后放置摄像机,将可以拍摄到由于试件变形导致光程差引起的干涉条纹,通过计算机对条纹进行处理可以提取相关的力学参数。
如果是透射相干梯度敏感方法,光程差δS(x,y)主要由应力改变导致的试件折射率变化和泊松效应导致的试件厚度变化两个因素产生,所以透射相干梯度敏感方法的光力学控制方程为:
cd ∂ ( σ ~ x + σ ~ y ) ∂ x = mp Δ - - - ( 1 )
其中m是条纹级数,p是光栅节距,Δ是两个光栅之间的间距,c是材料的光学常数,d是试件厚度,
Figure C20061000136600062
代表试件厚度方向上的平均应力分量。
光力学控制方程(1)表明当光栅主方向平行于x轴时,在透射情况下相干梯度敏感条纹代表的物理意义为试件主应力和对x方向梯度的等值线。
对于反射相干梯度敏感方法,光程差δS(x,y)仅由泊松效应导致的试件厚度变化产生,所以反射相干梯度敏感方法的光力学控制方程为:
∂ w ∂ x = mp 2 Δ - - - ( 2 )
其中m是条纹级数,p是光栅节距,Δ是两个光栅之间的间距,w代表试件的离面位移。
光力学控制方程(2)表明当光栅主方向平行于x轴时,在反射情况下相干梯度敏感干涉条纹代表的物理意义为试件离面位移对x方向梯度的等值线。
本发明方法的光路如图4所示,利用半反半透镜将试件4出射的平行光束分为两路光线,每一路光线都经过两个光栅产生相干梯度敏感干涉条纹,这两个干涉条纹场的光强分布分别为:
                   I1=a1+b1cos2πm1                         (3)
                   I2=a2+b2cos2πm2                         (4)
其中a代表背景光强,条纹图按照余弦规律bcos2πm变化。由于是不同光路形成的干涉条纹场,所以它们的背景光强a和条纹光强b互相不同。为了便于后续处理,将光强分布变换为以下形式:
                   I1-a1=b1cos2πm1                         (5)
                   I2-a2=b2cos2πm2                         (6)
其中背景光强a1、a2通过试件未加载时两个光路拍摄的数字图像获得。当相干梯度敏感干涉条纹图减去背景光强的影响以后,通过对数字图像进行直方图变换处理,可以改变两幅图像的条纹光强b1、b2使它们相同为b,但是条纹级数m没有改变。经过数字图像方法处理以后两个光路的干涉条纹图光强分布为I1 *和I2 *
                    I1 *=bcos2πm1                 (7)
                    I2 *=bcos2πm2                 (8)
根据相干梯度敏感方法的光力学控制方程(1)与(2)式,两个平行光栅的间距Δ与干涉条纹级数m成正比:
Δ 1 Δ 2 = m 1 m 2 - - - ( 9 )
如果采用图4光路拍摄同一试件在不同光栅间距下的两幅相干梯度敏感干涉条纹图,结合式(7)、(8)与(9)可以得到两个光栅间距Δ与经过数字图像方法处理后的干涉条纹光强I*和条纹级数m1之间的关系:
I 1 * I 2 * = cos ( 2 π m 1 ) cos ( 2 π m 1 Δ 2 / Δ 1 ) - - - ( 10 )
根据拍摄两幅条纹图时的光栅间距Δ和两幅数字图像方法处理后的干涉条纹图的光强分布I*,利用(10)式就可以获得第一幅条纹图上任一点的精确条纹级数m1
以反射相干梯度敏感方法实验为例对光栅间距法提取条纹级数的效果进行检验。图5所示为同一试件采用图4光路拍摄的两幅相干梯度敏感干涉条纹图,光栅间距分别是Δ1=50mm和Δ2=100mm,此时根据(10)式得:
I 1 * I 2 * = cos ( 2 π m 1 ) cos ( 4 π m 1 ) - - - ( 11 )
对图5中两幅干涉条纹图像进行图像处理,减去背景光强的影响并且使它们条纹光强相同,然后根据(11)式计算获得Δ1=50mm光路拍摄的条纹图中任意一点的条纹级数。为了与直接采用图像处理提取条纹中心线的结果进行对比,图6显示了采用光栅间距法前后的误差情况。由图6可以看到,与直接通过图像处理提取条纹中心线的方法相比光栅间距法明显减小了条纹提取带来的误差,改善了相干梯度敏感方法的计算精度,并且光栅间距法对于低级数条纹效果明显。

Claims (2)

1、一种提取相干梯度敏感干涉条纹级数的方法,其特征在于该方法包括以下步骤:
(1)对一束激光进行扩束,得到发散光,再经凸透镜成为平行光,入射到试件表面;
(2)将上述试件出射的平行光束分为两路光线,第一路光线经过间距为Δ1的第一光栅和第二光栅,产生相干梯度敏感干涉条纹,第二路光线经过间距为Δ2的第三光栅和第四光栅,产生相干梯度敏感干涉条纹;
(3)对上述两个干涉条纹场成像,并拍摄干涉条纹场;
(4)在试件未加载前,分别从上述两个光路拍摄试件的数字图像,获得干涉条纹图的背景光强分别为a1、a2
(5)对试件施加载荷,分别从上述两个光路拍摄加载试件的相干梯度敏感干涉条纹图,两个干涉条纹场的光强分布分别为:I1、I2
(6)从上述步骤(5)的加载试件的两幅相干梯度敏感干涉条纹图中减去上述步骤(4)的相应背景光强,并对其进行直方图变换,使两幅图像条纹光强的幅度相同,变换获得的两个光路的干涉条纹图光强分布为I1 *和I2 *
(7)根据上述光栅间距和光强,计算第一幅干涉条纹图的条纹级数m1为:
I 1 * I 2 * = cos ( 2 πm 1 ) cos ( 2 πm 1 Δ 2 / Δ 1 ) .
2、一种提取相干梯度敏感干涉条纹级数的装置,其特征在于包括:
(1)用于产生激光的激光器,用于对激光进行扩束的扩束镜,使发散光成为平行光的凸透镜;
(2)用于将试件出射的平行光束分为两路光线的半反半透镜,用于产生第一路相干梯度敏感干涉条纹的第一光栅和第二光栅,第一光栅和第二光栅的间距为Δ1,用于产生第二路相干梯度敏感干涉条纹的第三光栅和第四光栅,第三光栅和第四光栅的间距为Δ2
(3)用于对上述两个干涉条纹场成像的第一组凸透镜和光阑以及第二组凸透镜和光阑,以及用于拍摄干涉条纹场的第一摄像机、第二摄像机和用于图像处理和干涉条纹级数提取的计算机;
上述激光器、扩束镜、凸透镜、试件、半反半透镜、第一光栅、第二光栅、第一组凸透镜和光阑和第一摄像机依次处于第一光轴上;上述第三光栅、第四光栅、第二组凸透镜和光阑和第二摄像机依次处于第二光轴上,第一光轴与第二光轴于半反半透镜处垂直相交,计算机分别与第一摄像机和第二摄像机相连接。
CNB2006100013666A 2006-01-20 2006-01-20 一种提取相干梯度敏感干涉条纹级数的方法及其装置 Expired - Fee Related CN1320340C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006100013666A CN1320340C (zh) 2006-01-20 2006-01-20 一种提取相干梯度敏感干涉条纹级数的方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006100013666A CN1320340C (zh) 2006-01-20 2006-01-20 一种提取相干梯度敏感干涉条纹级数的方法及其装置

Publications (2)

Publication Number Publication Date
CN1800835A CN1800835A (zh) 2006-07-12
CN1320340C true CN1320340C (zh) 2007-06-06

Family

ID=36810965

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006100013666A Expired - Fee Related CN1320340C (zh) 2006-01-20 2006-01-20 一种提取相干梯度敏感干涉条纹级数的方法及其装置

Country Status (1)

Country Link
CN (1) CN1320340C (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101683407B1 (ko) * 2012-10-05 2016-12-06 고쿠리츠다이가쿠호우징 카가와다이가쿠 분광 측정 장치
US9482576B2 (en) 2012-10-05 2016-11-01 National University Corporation Kagawa University Spectroscopic measurement device having transmissive optical member with a sloped face
CN105547540B (zh) * 2015-12-09 2019-03-05 清华大学 实时空间相移的相干梯度敏感干涉方法
CN110672027B (zh) * 2019-09-30 2020-07-03 清华大学 一种基于相干梯度敏感干涉的条纹自动处理装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020012122A1 (en) * 2000-03-27 2002-01-31 Boyd David A. Coherent gradient sensing ellipsometer
CN1598541A (zh) * 2004-08-03 2005-03-23 中国海洋大学 一种测量生油岩成熟度的装置与方法
CN1639540A (zh) * 2001-12-31 2005-07-13 沈阳同联集团高新技术有限公司 投影莫尔条纹三维表面形状测量的方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020012122A1 (en) * 2000-03-27 2002-01-31 Boyd David A. Coherent gradient sensing ellipsometer
CN1639540A (zh) * 2001-12-31 2005-07-13 沈阳同联集团高新技术有限公司 投影莫尔条纹三维表面形状测量的方法和装置
CN1598541A (zh) * 2004-08-03 2005-03-23 中国海洋大学 一种测量生油岩成熟度的装置与方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"一种确定压剪全息干涉条纹级数的理论与实验研究" 陈顺云等,固体力学学报,第25卷第3期 2004 *

Also Published As

Publication number Publication date
CN1800835A (zh) 2006-07-12

Similar Documents

Publication Publication Date Title
CN103528524B (zh) 透视测量树脂基复合材料内部离面位移场分布的方法
CN109269438B (zh) 用于多层复杂微纳结构检测的结构光照明显微测量方法
CN104160241A (zh) 利用高维亮度信息的条纹图像的相位分布分析方法、装置及其程序
CN107144525A (zh) 复杂非均质结构体的动态应力场演化规律的测量方法
Larsson et al. Microscopic 3-D displacement field measurements using digital speckle photography
CN1320340C (zh) 一种提取相干梯度敏感干涉条纹级数的方法及其装置
Sciammarella Overview of optical techniques that measure displacements: Murray lecture
CN1975321A (zh) 移相干涉图像的信息处理系统及其处理方法
CN102865810B (zh) 基于正交双光栅的同步相移共光路干涉检测装置及检测方法
CN1841007A (zh) 条纹计测装置的变换系数校正方法和装置及条纹计测装置
Pástor et al. The use of optical methods in the analysis of the areas with stress concentration
CN101320003B (zh) 一种提高相干梯度敏感测量精度的方法
Ramesh Photoelasticity
CN1553139A (zh) 微透镜结构参数及面形畸变的快速检测方法
CN100427930C (zh) 一种提取相干梯度敏感干涉条纹级数的方法
Chaudhary et al. Measurement of the size and spacing of standard wire sieves using an image processing system and wavelet transform
CN2826363Y (zh) 双波长二维空间相移电子散斑干涉仪
CN2651745Y (zh) 具有相移功能的三维电子散斑干涉仪
Kim Identification of the local stiffness reduction of a damaged composite plate using the virtual fields method
Yu et al. Effect of average speckle size on speckle interference fringe pattern
CN1283976C (zh) 双波长二维空间相移电子散斑干涉仪
Rashidifar et al. Strain Measurements by Full-Field Optical Techniques
Agrawal Full-field deformation measurement in wood using digital image processing
CN1609599A (zh) 具有相移功能的三维电子散斑干涉仪
Hao et al. Mechanical Performance Characterization of EMI Shielding Materials Using Optical Experimental Techniques

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20070606

Termination date: 20100220