CN100590425C - 一种透明平板材料应力检测方法 - Google Patents

一种透明平板材料应力检测方法 Download PDF

Info

Publication number
CN100590425C
CN100590425C CN200610031731A CN200610031731A CN100590425C CN 100590425 C CN100590425 C CN 100590425C CN 200610031731 A CN200610031731 A CN 200610031731A CN 200610031731 A CN200610031731 A CN 200610031731A CN 100590425 C CN100590425 C CN 100590425C
Authority
CN
China
Prior art keywords
centerdot
stress
light
glass
light beam
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
CN200610031731A
Other languages
English (en)
Other versions
CN101082590A (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.)
CHANGSHA CREATOR SYSTEM INTEGRATION Co Ltd
Central South University
Original Assignee
CHANGSHA CREATOR SYSTEM INTEGRATION Co Ltd
Central South 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 CHANGSHA CREATOR SYSTEM INTEGRATION Co Ltd, Central South University filed Critical CHANGSHA CREATOR SYSTEM INTEGRATION Co Ltd
Priority to CN200610031731A priority Critical patent/CN100590425C/zh
Publication of CN101082590A publication Critical patent/CN101082590A/zh
Application granted granted Critical
Publication of CN100590425C publication Critical patent/CN100590425C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

一种透明平板或模型材料应力检测的方法和装置,由光源1、光学系统2、接收器4,5和计算机——图象处理系统6组成。本发明根据两个参照光束与两个不同相位测量光束的光强,利用它们的和差比值求得被测应力角,该方法自动约除了光源波动、光路扰动和光电变换漂移等因素对测量结果的影响。装置的所有部件都在被测玻璃板的一侧,故无须传统方法中的同步驱动和精密补偿机构,运行可靠,实时性强,测量范围大,精度高。

Description

一种透明平板材料应力检测方法
技术领域
本发明涉及一种测量透明平板材料应力的方法和装置。
背景技术
众所周知,透明平板或模型材料如玻璃、有机玻璃,聚碳酸脂、环氧树脂等材料具有暂时双折射性能。国内外的透明平板或模型材料非接触应力测量装置,均采用光弹性应力分析方法。以透明平板玻璃为例,试验测试表明,透明平板玻璃主应力差通常小于I级。故对玻璃表面应力的测量问题,要从光强来确定它的主应力差(σ12)大小,属于非整数条纹级处理范畴,不能应用通常光弹方法中的整数条纹级分析技术。
对于这种非整数条纹级数处理,国内外出现的各种装置,无论在线或离线的玻璃表面应力测量装置,大都采用补偿法。
补偿法就是在线偏振或圆偏振布置的应力光学检测系统中,当光路中未插入被测应力物体时,系统处于平衡(如消光)状态。一旦插入被应力物体,就破坏了原来的平衡,可旋转偏振片或插入一个应力光学补偿器到检测系统中,直至系统恢复到平衡状态(消光状态)。其偏振片旋转的角度或补偿器插入的距离就是对被测对象主应力差的直接描述,如赛纳蒙(Senarmont)法、巴比纳特(Babinet)法等方法。
补偿法的共同特点是必须具有一个旋转偏振片或移动补偿器的精密机械装置,它的机械精密程度直接影响测量精度。
另外,对于在线大视场测量玻璃带的应力,玻璃带宽度通常在3米以上,通常将偏振布置的光学系统分为发射和接受两部分,分别安装在被测玻璃带上、下两侧,由电动机带动它们横向移过玻璃带,得到对玻璃带一次横扫描过程的应力测量。对这种在线应力检测系统,必需有一套在移动过程中保证玻璃带两侧的设备对中并同步横向移动的机械装置。这不仅增加了装置的复杂程度,又降低了它的可靠性。
总之,国内外应用于透明平板如玻璃工业的应力测量装置都需要精密机械补偿机构和伺服电机调零机构;被测玻璃板上下表面都有部件,需要上下部件对中同步移动的同步传动机构,由此带来维护困难,运行可靠性差,测量精度难以保证。
发明内容
为了提高玻璃应力测量装置的可维护性和测量精度,解决大视场应力测量问题,本发明提供了一种透明平板玻璃应力在线自动检测方法和装置,该装置方法将微型计算机-数字图象处理技术应用于工业现场的在线应力检测。
本发明装置可分为四部分组成:光源1、光学系统2、接收器4,5和微型计算机——图象处理系统6。如附图1所示。
光学系统2由起偏器21、1/4波片22、偏振保持光栅23、空间滤波器24、分析镜25、26、聚焦透镜27、28组成;接收器4由CCD感光单元41,接收器5CCD感光单元51组成;计算机图象采集处理系统6由数字图象采集模块、通用或嵌入式计算机及相关显示与通信模块组成。
放置在起偏器21与被测应力的玻璃3之间的1/4波片22,其快轴与起偏器21偏振轴成45°;穿过1/4波片22的光束经过偏振保持光栅23和空间滤波器24后,分裂成与被测玻璃3表面法线呈左右对称的两束光以小角度的方式入射,经被测玻璃3的上、下表面反射后变为四束光,射向光学系统2光路中对应的两个分析镜25和26后,分别投射到CCD感光单元41和51;分析镜25和26的偏振轴分别与起偏器2偏振轴线的夹角为+θ和-θ,并取θ=45°分析镜25的偏振轴分别与起偏器偏振轴的夹角为+θ和-θ,并取θ=45°;接收器4,5使用光学系统2中的聚焦透镜27,28,将近似共路的参照光束与测量光束在空间上分离,并采用同一CCD感光变换单元41、51对同一侧的参照光束与测量光束进行光电变换;接收器4,5同时摄取两幅偏振光图像,经A/D转换后,形成数字图象存入计算机图象处理系统6,计算机对采集的数字图像进行运算处理,得到玻璃带横向各点的应力差;根据投射到接收器4,5、CCD感光变换单元41、51上的光束强度I1u、I1d和I2u、I2d,其中I1u、I1d分别对应于左侧入射经透明平板材料(3)上下表面反射光束,I2u、I2d对应于右侧入射经透明平板3上下表面反射光束,计算从分析镜出来的光强为:
I 1 u : = k 1 · A 1 2 2
I 1 d : = k 1 · τ 1 · A 1 2 2 · [ 1 - sin · ( δ ) ]
I 2 u : = k 2 · A 2 2 2
I 2 d : = k 2 · τ 2 · A 2 2 2 · [ 1 + sin · ( δ ) ] - - - ( 1 )
式中A1、A2为滤波器的出射光幅度,k1为左光束到达CCD感光单元41时的光路衰减系数,经玻璃3上下表面反射光光束可以近似认为共光路,忽略其差别,取相同的光路衰减系数,同理,k2为右光束到达CCD感光单元51时的光路衰减系数。τ1、τ2为玻璃3对左右下表面反射光光束的吸收系数,相对生产线上3米多宽玻璃板来说,反射发生在很小的区域,且左右光束对称入射,因此我们可以假定τ1=τ2。δ为被测玻璃体内主应力差引起的相位差亦称应力角。
由(1)式可得:
sin ( δ ) : = I 1 d · I 2 u - I 2 d · I 1 u I 1 d · I 2 u + I 2 d · I 1 u - - - ( 2 )
由此得到应力角分析式:
δ : = a sin · ( I 1 d · I 2 u - I 2 d · I 1 u I 1 d · I 2 u + I 2 d · I 1 u ) - - - ( 3 )
光束强度I1u、I1d和I1u、I2d经CCD感光单元41和51光电变换成视频信号,视频信号计算机图象采集处理系统处理后,得到相对应的四个灰度值,记为g1u、g1d和g2u、g2d。由于(3)式中的比例关系,光电变换增益系数被自动约除,有:
δ : = a sin · [ g 1 d · g 2 u - g 2 d · g 1 u ( g 1 d · g 2 u + g 2 d · g 1 u ) ] - - - ( 4 )
本发明利用被测平板玻璃3本身的反射分束特性,上表面反射产生参照光束,下表面反射产生测量光束;参照光束与测量光束共光源共光路,因而能自动消除光源的波动和光路的扰动对测量结果的影响。采用分析镜两个不同安置角度(±45°),测取两个不同相位的偏振光的光强,由(4)式运算出被测应力角大小。采用聚焦透镜27和28在空间上分离参照光束与测量光束,采用同一CCD感光单元对一侧的参照光束与测量光束进行感光;保证其光电变换的一致性,同时允许测量过程中光束的空间取向有一定变化。
我们称基于(3)和(4)式的测量原理为相位法。采用(4)式的相位测量方法,它不再含有光幅度A这个量,即其测量精度不受光功率的变化影响,亦不受光电变换增益系数影响。克服了光度法的困难,即消除了光功率随机波动和光电变换增益漂移的问题。
从式(4)可以看出,当g1d*g2u-g2d*g1u>0时,应力角为正值,它表示为张应力;反之则为压应力;两者相等时,则应力为零。在不旋转偏振片的情况下就能判别出拉张或压应力,无需精密补偿机构和同步机构,就可以测得应力值。
相位测量法还有另一优点,就是测量系统的传感器工作区刚好在它的灵敏区,而不是象补偿法那样传感器工作在(消光状态)光强等于零附近的不灵敏区。而相位测量法从式(4)可知,应力角为零时必g1d*g2u=g2d*g1u。又从式(1)可知,该时测量光束光强I1d和I2d为其最大值的一半,参考光束光强I1u和I2u与值应力值无关。通常玻璃表面应力值在±45°以内,故它们必然处在装置的测量灵敏线性区。
图象处理系统6包括图象采集处理器、计算机、打印机等组成。图象采集处理器同时采集来自CCD感光单元41和51的两幅偏振光强图像,这两幅偏振光强图像分别包含四束光光强度I1u、I1d和I2u、I2d信号。计算机图象采集处理系统6从数字图像中提取出相对应的四个灰度值g1u、g1d和g2u、g2d,按(4)式进行运算处理,得出玻璃带横向各点的应力差即应力角δ,在监示器上显示出被测玻璃应力差分布曲线,并可由打印机打印记录。
为了实现对玻璃带宽3~4米以上的大视场应力测量,本发明将所有测量器件封装为应力测量盒,由导轨及驱动系统驱动扫描过被检玻璃板,得到整个玻璃断面应力分布。由于在被检玻璃板另一侧没有测量部件,故不需要同步对中机构。
在预先编制好的程序控制下,整套应力检测装置可连续自动在线测量。
本发明将计算机——数字图象处理技术应用于工业生产现场在线应力自动检测,测量快速,操作方便,实时性强,显示直观。无需精密补偿机构和同步传动机构,就能实时地在线测量大视场的应力大小和方向。并且运行可靠性高,使用维护简单,造价低廉。测量时与光幅度无关,综合误差小,而且系统误差又能很方便进行补偿,装置定标方便,测量精度高。测量光束两次通过被测玻璃板,检测灵敏度增加一倍。本测量方法也同样适应于其他透明平板和模型材料的应力测量。
附图说明
图1为本发明光路原理图;
图2为实施例图。
实施方式
以玻璃工业生产应力测量为例。如附图2所示,为了实现对板宽为3~4米以上玻璃的大视场应力测量,首先在被测玻璃的输送辊道7上架设测量桥9,输送辊道7上有工业生产的被测玻璃3,将附图1中除被测玻璃3以外的所有器件封装为附图2中的应力测量盒10,由导轨及驱动系统9驱动扫描过被检玻璃板3,完成整个玻璃断面应力分布测量。
如附图1所示,光源1采用准直半导体激光器,波长650nm;波片22为对应于650nm波长的1/4波片,偏振保持光栅23采用普通金属光栅,滤波器24采用带两孔的不透明片,只允许光栅的±1级光透过。起偏器21和分析镜25和26采用高比偏振片,聚焦透镜27,28采用大口径低渐晕透镜、CCD感光单元41,51为两台工业摄像机(取γ值为1),选用每一行像素不少500的面阵(或线阵)摄像机,例如TM765型CCD摄像机,矩形像素765×581,低照度灵敏度0.3Lux。它放置光学系统2聚焦透镜27,28的后面。
设起偏器21的偏振轴与玻璃带3的运动方向相同,放置在起偏器21与被测应力的玻璃3之间的1/4波片22,其快轴与起偏器21偏振轴成45°;穿过波片22的光束经过偏振保持光栅23和空间滤波器24后,分裂成与玻璃3表面法线呈左右对称的两束光以小角度的方式入射,经玻璃3的上下表面反射后变为四束光,分别射向光路中对应的分析镜25、26后,分别投射到CCD感光单元41和51,分析镜25、26的偏振轴分别与X轴(即输送辊子轴线,如图2所示)夹角为+θ和-θ,并令θ=45°。
投射到CCD感光单元41和51的光束强度分别记为I1u、I1d和I2u、I2d,I1u、I1d分别对应于左侧入射经玻璃3上下表面反射光束;I1u、I2d对应于右侧入射经玻璃3上下表面反射光束。从检偏镜出来的光强为:
I 1 u : = k 1 · A 1 2 2
I 1 d : = k 1 · τ 1 · A 1 2 2 · [ 1 - sin · ( δ ) ]
I 2 u : = k 2 · A 2 2 2
I 2 d : = k 2 · τ 2 · A 2 2 2 · [ 1 + sin · ( δ ) ] - - - ( 1 )
(1)式中A1、A2为附图1中滤波器24的出射光幅度,k1为左光束到达CCD感光单元41时的光路衰减系数,经玻璃3上下表面反射光光束可以近似认为共光路,忽略其差别,取相同的光路衰减系数,同理,k2为光束到达CCD感光单元51时的光路衰减系数。τ1、τ2为玻璃4对左右光束的吸收系数,相对生产线上3米多宽玻璃板来说,反射发生在很小的区域,且左右光束对称入射,因此可以假定τ1=τ2。δ为被测玻璃体内主应力差引起的相位差即应力角。
由(1)式可得:
sin ( δ ) : = I 1 d · I 2 u - I 2 d · I 1 u I 1 d · I 2 u + I 2 d · I 1 u - - - ( 2 )
δ : = a sin · ( I 1 d · I 2 u - I 2 d · I 1 u I 1 d · I 2 u + I 2 d · I 1 u ) - - - ( 3 )
为了在同一时刻对同一干涉图象摄像,两台摄像机必须在计算机的控制下同步采集。
当摄像机将摄取的光束强度I1u、I1d和I2u、I2d转换成电量后,经视频信号输送给计算机图象处理系统6,计算机从数字图像中提取出相对应的四个灰度值g1u、g1d和g2u、g2d,按(4)式进行运算处理,得出玻璃带横向各点的应力差即应力角δ,同时计算机控制附图2中的驱动系统4完成对玻璃板的扫描,就可以得到玻璃带横向分布的应力值,在监示器上显示出被测玻璃应力差分布曲线,并可由打印机打印记录。

Claims (1)

1.一种透明平板应力的检测方法,其特征在于:采用放置在起偏器(21)与被测平板材料(3)之间的1/4波片(22),其快轴与起偏器偏振轴成45°,穿过1/4波片(22)的光束经过偏振保持光栅(23)和空间滤波器(24)后,分裂成与平板材料(3)表面法线呈左右对称的两束光以小角度的方式入射,经透明平板材料(3)上下表面反射后变为四束光,射向光学系统(2)中对应的两个分析镜(25),(26)后分别投射到不同的接收器(4),(5);分析镜(25)的偏振轴分别与起偏器偏振轴的夹角为+θ和-θ,并取θ=45°;接收器(4),(5)使用光学系统(2)中的聚焦透镜(27),(28),将近似共路的参照光束与测量光束在空间上分离,并采用同一CCD感光变换单元(41)、(51)对同一侧的参照光束与测量光束进行光电变换;接收器(4),(5)同时摄取两幅偏振光图像,经A/D转换后,形成数字图象存入计算机图象处理系统(6),计算机对采集的数字图像进行运算处理,得到玻璃带横向各点的应力差;根据投射到接收器(4),(5)、CCD感光变换单元(41)、(51)上的光束强度I1u、I1d和I2u、I2d,其中I1u、I1d分别对应于左侧入射经透明平板材料(3)上下表面反射光束,I2u、I2d对应于右侧入射经透明平板(3)上下表面反射光束,计算从分析镜出来的光强为:
I 1 u : = k 1 · A 1 2 2 I 1 d : = k 1 · τ 1 · A 1 2 2 · [ 1 - sin · ( δ ) ]
I 2 u : = k 2 · A 2 2 2 I 2 d : = k 2 · τ 2 · A 2 2 2 · [ 1 + sin · ( δ ) ] - - - ( 1 )
式中A1、A2为滤波器的出射光幅度,k1为左光束到达CCD8时的光路衰减系数,k2为光束到达CCD8时的光路衰减系数,τ1、τ2为玻璃对左右光束的吸收系数,δ为被测玻璃体内主应力差引起的相位差,再根据(1)式得到应力公式:
sin ( δ ) : = I 1 d · I 2 u - I 2 d · I 1 u I 1 d · I 2 u + I 2 d · I 1 u
δ : = a sin · [ g 1 d · g 2 u - g 2 d · g 1 u ( g 1 d · g 2 u + g 2 d · g 1 u ) ] .
CN200610031731A 2006-05-30 2006-05-30 一种透明平板材料应力检测方法 Expired - Fee Related CN100590425C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200610031731A CN100590425C (zh) 2006-05-30 2006-05-30 一种透明平板材料应力检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200610031731A CN100590425C (zh) 2006-05-30 2006-05-30 一种透明平板材料应力检测方法

Publications (2)

Publication Number Publication Date
CN101082590A CN101082590A (zh) 2007-12-05
CN100590425C true CN100590425C (zh) 2010-02-17

Family

ID=38912273

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200610031731A Expired - Fee Related CN100590425C (zh) 2006-05-30 2006-05-30 一种透明平板材料应力检测方法

Country Status (1)

Country Link
CN (1) CN100590425C (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806966A (zh) * 2010-03-31 2010-08-18 苏州达信科技电子有限公司 检测装置及检测方法

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101592537B (zh) * 2009-07-10 2012-12-12 成都光明光电股份有限公司 光学玻璃应力测量装置及其测量方法
CN102062656A (zh) * 2010-12-31 2011-05-18 肖天长 玻璃表面应力测试仪
CN103443603B (zh) * 2011-03-18 2016-05-25 旭硝子株式会社 玻璃的表面应力测定装置和玻璃的表面应力测定方法
JP5751036B2 (ja) * 2011-06-09 2015-07-22 旭硝子株式会社 強化ガラス及びその製造方法、該強化ガラスの表面応力測定方法
CN104316233A (zh) * 2014-11-03 2015-01-28 苏州精创光学仪器有限公司 表面应力测量装置
CN104535232B (zh) * 2014-11-18 2017-02-01 湖北新华光信息材料有限公司 一种红外光学材料应力测试装置及其测试方法
CN105509939B (zh) * 2015-12-01 2018-07-27 长春理工大学 一种实时检测激光损伤非透明材料的应力的装置与方法
CN105675186B (zh) * 2016-01-26 2018-09-18 中国科学院声学研究所 基于动态光弹性系统的应力测量方法
CN108106759A (zh) * 2016-11-25 2018-06-01 北京自动化控制设备研究所 一种原子气室表面应力测试装置及测试方法
CN107991007A (zh) * 2018-01-02 2018-05-04 洛阳兰迪玻璃机器股份有限公司 一种玻璃板应力检测装置
CN107677402A (zh) * 2017-08-22 2018-02-09 洛阳兰迪玻璃机器股份有限公司 一种玻璃板应力在线检测方法
CN108051443B (zh) * 2017-11-21 2020-06-30 中国科学院上海光学精密机械研究所 一种玻璃应力及缺陷的测量装置和测量方法
CN108760112B (zh) * 2018-05-25 2020-06-30 中国科学院上海光学精密机械研究所 基于衍射重叠迭代算法的应力测量装置和方法
CN110530800B (zh) * 2018-05-25 2022-03-29 深圳市杰普特光电股份有限公司 玻璃应力缺陷的检测方法和装置
CN112304904B (zh) * 2019-07-15 2023-11-03 松山湖材料实验室 基于滤波阵列的硅片反射率检测方法
CN110376213B (zh) * 2019-08-14 2022-01-25 深圳市麓邦技术有限公司 光学检测系统及方法
CN110793977B (zh) * 2019-10-28 2022-05-31 芜湖东旭光电科技有限公司 玻璃基板对比度检测装置
CN112729143B (zh) * 2020-12-16 2022-03-22 江苏大学 基于比率荧光技术的发光光弹性涂层应变测量系统及方法
CN113514178B (zh) * 2021-04-16 2022-04-15 江南大学 一种基于三波长照明光源的光弹性应力测量系统与方法
CN113447173B (zh) * 2021-05-17 2022-08-30 东北石油大学 一种在光弹性实验中直接测量有弧形缺陷连续油管应力的方法
CN115077755A (zh) * 2022-06-14 2022-09-20 秦皇岛玻璃工业研究设计院有限公司 玻璃生产线冷端应力检测系统及玻璃应力检测方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806966A (zh) * 2010-03-31 2010-08-18 苏州达信科技电子有限公司 检测装置及检测方法

Also Published As

Publication number Publication date
CN101082590A (zh) 2007-12-05

Similar Documents

Publication Publication Date Title
CN100590425C (zh) 一种透明平板材料应力检测方法
CN102967263B (zh) 一种桥梁挠度-转角一体化测量方法
CN103776445B (zh) 分振幅偏振导航角度传感设计方法及装置
CN103822714B (zh) 一种基于偏振分光器的高空间分辨率快照式成像光谱仪与成像方法
CN102426058B (zh) 一种静态干涉成像偏振仪及获得目标的偏振信息的方法
CN101799318A (zh) 一种激光零差测振光学系统及其信号处理方法
CN103822715B (zh) 一种基于偏振分光器的快照式成像光谱仪与成像方法
CN104457600A (zh) 一种光纤准直器阵列的测试装置
CN103017908A (zh) 基于四路分光模块的偏振光特性实时测量装置及方法
CN102288104A (zh) 六轴4细分干涉仪
CN102538971A (zh) 全光场全斯托克斯参量检测装置和检测方法
CN109579744B (zh) 基于光栅的跟随式三维光电自准直方法与装置
CN106247989A (zh) 一种导轨滚转角现场标定及测量装置及方法
CN106323198B (zh) 一种高精度、宽范围和大工作距激光自准直装置与方法
CN202420680U (zh) 全光场全斯托克斯参量检测装置
CN106247992B (zh) 一种高精度、宽范围和大工作距自准直装置与方法
CN106225730B (zh) 便携式组合调零高精度激光大工作距自准直装置与方法
CN107036527A (zh) 同步测量绝对寻址距离与偏摆角度的光学系统与方法
CN103900693B (zh) 一种差分快照式成像光谱仪与成像方法
CN106323197B (zh) 便携式阵列调零高精度激光大工作距自准直装置与方法
CN103176297A (zh) 液晶单元特性测定装置以及液晶单元特性测定方法
CN106017441A (zh) 一种便携式高精度激光大工作距自准直装置与方法
CN105136295A (zh) 一种aotf同一幅图中光谱不均匀的解决方法及装置
CN101561317B (zh) 一种高精度自动测量目标对传输光束偏振态影响的系统
CN109187368A (zh) 一种基于多通道测量的液晶偏振特性实时检测系统及方法

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: 20100217

Termination date: 20100530