CN103234496B - 一种三坐标测量机二维平台误差的高精度校正方法 - Google Patents
一种三坐标测量机二维平台误差的高精度校正方法 Download PDFInfo
- Publication number
- CN103234496B CN103234496B CN201310106750.2A CN201310106750A CN103234496B CN 103234496 B CN103234496 B CN 103234496B CN 201310106750 A CN201310106750 A CN 201310106750A CN 103234496 B CN103234496 B CN 103234496B
- Authority
- CN
- China
- Prior art keywords
- coordinate
- measuring machine
- error
- coordinate measuring
- turbogrid plates
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000005259 measurement Methods 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000013519 translation Methods 0.000 claims description 16
- 230000036544 posture Effects 0.000 abstract 1
- 238000009826 distribution Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000009966 trimming Methods 0.000 description 1
- 238000007514 turning Methods 0.000 description 1
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- A Measuring Device Byusing Mechanical Method (AREA)
Abstract
本发明公开了一种三坐标测量机二维平台误差的高精度校正方法。利用精度要求低于或等于待测三坐标测量机二维平台的刚性栅格板作为辅助测量装置,并根据测得的六位姿状态下坐标测量机上各个标记点的坐标,运用基于最小二乘法的自校正算法将待测二维平台误差以及所使用的栅格板标尺误差从原始测量数据中分离出来,由此可实现对三坐标测量机二维平台的高精度校正。本发明可有效地获得待测的三坐标测量机二维平台误差,其测量精度可达到亚微米量级;同时该发明无需昂贵的专用高精度辅助装置,具有较高的可靠性,为三坐标测量机二维平台误差提供了一种高精度的校正方法,并具有极高的实际应用价值。
Description
技术领域
本发明属于精密测量领域,尤其涉及一种三坐标测量机二维平台误差的高精度校正方法。
背景技术
超精密加工技术的高速发展对相应的三维检测设备的测量精度提出了更高的要求。三坐标测量机(Coordinate Measuring Machining,CMM)作为传统的通用型高精度测量仪器,在工件形位误差的检测中有着不可替代的作用,并正朝着尺寸小型化以及纳米级精度的方向发展。超精密三坐标测量机研究的关键是对其进行形位误差测量与不确定度评定,进而保证其微纳米级的测量精度。
三坐标测量机是一个具有较多误差源的复杂***,国内外学者对其误差修正方法进行了广泛的研究。为实现测量误差的高精度校正,必须准确得到测量机各个单项误差。根据坐标测量机校准规范JJF1064-2010的要求,校准尺寸测量示值误差时一般使用量块或激光干涉仪(大型坐标测量机的补充测量使用激光干涉仪进行位置示值误差测量)作为标准器。通过测量标准器上标记点的测量值与标准值之间的差异得出被标定的示值误差,并进行拟合和补偿,进而实现对待测坐标测量机的校正。但该方法使用一维标准检具进行单项误差检测,耗时长、专用仪器及装置造价昂贵、数据处理繁琐。尤其是当待测的精密平台为超大规模集成电路制造、高密度存储设备加工以及光纤对接等超精密加工领域使用的纳米级或者亚纳米级超精密工作台时,传统的坐标测量机标定技术因无法找到更高精度的标准计量仪器而无法应用。
发明内容
本发明要解决的技术问题是利用精度等级较低的栅格板有效地分离出三坐标测量机的二维平台误差以及所使用的栅格板标尺误差,进而实现二维平台的高精度校正。
一种三坐标测量机二维平台误差的高精度校正方法如下:
1)选用n×n的方形栅格板,n为自然数,将该栅格板固定在三坐标测量机平台上,并且栅格方向与坐标测量机的x、y轴导轨运动方向对齐;
2)以栅格板为基准建立坐标系,利用该三坐标测量机测得原始位姿状态下对应栅格板上所有标记点中心位置的坐标数据M1;
3)将栅格板相对于原始位姿分别沿x轴正、负方向各平移一个栅格距离,利用三坐标测量机分别测得栅格板上所有标记点中心位置的坐标数据M2和M3;
4)将栅格板平移回步骤2)中的原始位姿,并将其分别相对原始位姿逆时针旋转至90°、180°和270°位置,利用三坐标测量机分别测得对应栅格板上所有标记点中心位置的坐标数据为M4、M5、M6;
5)将步骤2)到步骤4)中所测得的六组坐标数据M1、M2、M3、M4、M5和M6代入相应的六位姿误差分离方程组,即:
其中I为主对角线为1的矩阵,I1、I2为I的x轴方向平移矩阵,R90、R180、R270为I的逆时针90°、180°和270°旋转矩阵,Nx、Ny为栅格板上测量点x和y坐标的标称值,Nxs和Nys以及Nyt和Nxt分别为Nx、Ny由于x轴正负方向平移而形成的矩阵;E为所有元素为1的矩阵。Ax、Ay为测量点的二维平台x和y坐标误差,Gx、Gy为测量点的栅格板x和y坐标标尺误差,Vi、Wi为两理想坐标系原点的x和y坐标偏移(i=1,2,...,6,下同),θi为两理想坐标系的角度偏移,并且有:
6)利用步骤2)到步骤4)所测得的六位姿数据M1、M2、M3、M4、M5、M6以及步骤5)中的方程组,得到二维平台误差Ax、Ay。
7)对三坐标测量机二维平台误差进行校正:
其中[Tx,Ty]T为校正之前三坐标测量机所测得的数据,T0则为校正之后的数据。
本发明的有益效果:本发明利用精度等级较低的栅格板来进行待测三坐标测量机二维平台上多个位姿点坐标的测量,并运用自校正算法有效地分离出三坐标测量机二维平台误差和所使用的栅格板标尺误差。本发明所提出的测量方法无需昂贵的专用高精度校正仪器及装置,具有较高的可靠性,并且测量方法简单,不涉及复杂的实际操作过程,适合于具有精密平台的坐标测量机二维平台误差的高精度校正。
附图说明
图1是三坐标测量机二维平台误差校正中的栅格板辅助装置图;
图2是待测三坐标测量机二维平台和栅格板的理想坐标系关系及***误差示意图;
图3是本发明实施例中在原始位姿状态下利用栅格板测得的三坐标测量机总体误差分布;
图4是本发明实施例中利用六位姿数测得的三坐标测量机二维平台误差分布。
具体实施方式
一种三坐标测量机二维平台误差的高精度校正方法如下:
1)选用n×n的方形栅格板,n为自然数,该栅格板的位置精度可低于待测三坐标测量机二维平台的位置精度,将该栅格板固定在三坐标测量机平台上,并且栅格方向与坐标测量机的x、y轴导轨运动方向对齐,如图1所示;
2)以栅格板为基准建立坐标系,利用该三坐标测量机测得原始位姿状态下对应栅格板上所有标记点中心位置的坐标数据M1;
3)将栅格板相对于原始位姿分别沿x轴正、负方向各平移一个栅格距离,利用三坐标测量机分别测得栅格板上所有标记点中心位置的坐标数据M2和M3;
4)将栅格板平移回步骤2)中的原始位姿,并将其分别相对原始位姿逆时针旋转至90°、180°和270°位置,利用三坐标测量机分别测得对应栅格板上所有标记点中心位置的坐标数据为M4、M5、M6;
5)根据误差源构建两理想坐标系,如图2所示,每一个标记点的误差可以表示如下:
其中Mx、My为测量值的横纵坐标,Nx、Ny为栅格板上测量点x和y坐标的标称值,Ax、Ay为测量点的二维平台x和y坐标误差,Gx、Gy为测量点的栅格板x和y坐标标尺误差,V、W为两理想坐标系原点的x和y坐标偏移,θ为两理想坐标系的角度偏移,白点为实际测量标记点,黑点为理想测量标记点。在上式中取:
根据原始位姿测得的数据,可以得到如下方程组:
式中,Q1x和Q1y为原始位姿测得数据,[·]T为该矩阵的转置。由此可知,若测量原始位姿上的n×n个点,未知数的个数为4n2+3,而方程数为2n2+7,此时方程数小于未知数的个数,因此方程组有无穷解。当位姿数为3时,未知数的个数为4n2+9,而方程数为6n2-2n+7,此时方程组应用最小二乘法有解。当位姿数继续增加时,方程数远大于未知数的个数,故使用最小二乘法计算得到的方程组解的精度越高。因此,实际应用中,位姿数不能少于3个,且位姿数的增加可以有效地提高计算精度。本发明选择位姿数为6进行误差校正。将步骤2)到步骤4)中所测得的六组坐标数据M1、M2、M3、M4、M5和M6代入相应的六位姿误差分离方程组,即:
其中I为主对角线为1的矩阵,I1、I2为I的x轴方向平移矩阵,R90、R180、R270为I的逆时针90°、180°和270°旋转矩阵,Nx、Ny为栅格板上测量点x和y坐标的标称值,Nxs和Nys以及Nyt和Nxt分别为Nx、Ny由于x轴正负方向平移而形成的矩阵;E为所有元素为1的矩阵。Ax、Ay为测量点的二维平台x和y坐标误差,Gx、Gy为测量点的栅格板x和y坐标标尺误差,Vi、Wi为两理想坐标系原点的x和y坐标偏移(i=1,2,...,6,下同),θi为两理想坐标系的角度偏移,并且有:
6)利用步骤2)到步骤4)所测得的六位姿数据M1、M2、M3、M4、M5、M6和步骤5)的方程组,可以得到二维平台误差Ax、Ay。
7)对三坐标测量机二维平台误差进行校正:
其中[Tx,Ty]T为校正之前三坐标测量机所测得的数据,T0则为校正之后的数据。
实施例
实施例中采用的精密二维平台为海克斯康的坐标测量机Global Classical的工作平台,刚性栅格板为海克斯康的配套栅格板,栅格板上标记点误差为±1mm,实验环境温度为20℃。三坐标测量机二维平台误差的校正过程为:
1)将栅格板固定在三坐标测量机平台上,并且栅格方向与坐标测量机的x、y轴导轨运动方向对齐,如图1所示。
2)以栅格板1为基准建立坐标系。对栅格板上120mm×120mm区域范围内的4×4个栅格点按六位姿数进行测量,每两个栅格间隔为L1=40mm,得到原始位姿直接测得值的总体误差如图3所示。为了实现对二维平台误差高精度校正,需要进一步采用本发明提出的高精度校正方法进行处理。利用该三坐标测量机测得栅格板在原始位姿状态下对应所有标记点中心位置的坐标数据M1。
3)将栅格板相对于原始位姿分别沿x轴正、负方向各平移一个栅格距离,利用三坐标测量机分别测得栅格板上所有标记点中心位置的坐标数据M2和M3。
4)将栅格板平移回步骤2)中的原始位姿,并将其分别相对原始位姿逆时针旋转至90°、180°和270°位置,利用三坐标测量机分别测得对应栅格板上所有标记点中心位置的坐标数据为M4、M5、M6。
5)将步骤2)到步骤4)中所测得的六组坐标数据M1、M2、M3、M4、M5和M6代入相应的六位姿误差分离方程组,即:
其中I为主对角线为1的矩阵,I1、I2为I的x轴方向平移矩阵,R90、R180、R270为I的逆时针90°、180°和270°旋转矩阵,Nx、Ny为栅格板上测量点x和y坐标的标称值,Nxs和Nys以及Nyt和Nxt分别为Nx、Ny由于x轴正负方向平移而形成的矩阵;E为所有元素为1的矩阵。Ax、Ay为测量点的二维平台x和y坐标误差,Gx、Gy为测量点的栅格板x和y坐标标尺误差,Vi、Wi为两理想坐标系原点的x和y坐标偏移(i=1,2,...,6,下同),θi为两理想坐标系的角度偏移,并且有:
6)经过处理得到三坐标测量机二维平台误差A在各个位置点的测量值如图4所示,其中1-16点为x轴误差值,17-32点为y轴误差值。通过对图4中各个测点误差的分析可知,三坐标测量机二维平台误差限为±0.002mm,栅格板标尺误差限为±0.8mm,与标称值相符。
当间隔为L1=40mm时,测量误差的概率分布由于近似为正态分布,扩展不确定度为U=k·u。二维平台误差只考虑示值误差引起的不确定度,则u=σ。取k=2,得到二维平台x、y方向误差不确定度分别为1.2892μm和1.4248μm,该值与MCV-500多普勒式激光干涉仪得到待测三坐标测量机二维平台测量值偏差分别为0.07μm(x轴)、0.03μm(y轴)。
7)对三坐标测量机二维平台误差进行校正:
其中 为校正之前三坐标测量机所测得的数据,T0则为校正之后的数据。
Claims (1)
1.一种三坐标测量机二维平台误差的高精度校正方法,其特征在于它的步骤如下:
1)选用n×n的方形栅格板,n为自然数,将该栅格板固定在三坐标测量机平台上,并且栅格方向与坐标测量机的x、y轴导轨运动方向对齐;
2)以栅格板为基准建立坐标系,利用该三坐标测量机测得原始位姿状态下对应栅格板上所有标记点中心位置的坐标数据M1;
3)将栅格板相对原始位姿分别沿x轴正、负方向各平移一个栅格距离,利用三坐标测量机分别测得栅格板上所有标记点中心位置的坐标数据M2和M3;
4)将栅格板平移回步骤2)中的原始位姿,并将其分别相对原始位姿逆时针旋转至90°、180°和270°位置,利用三坐标测量机分别测得对应栅格板上所有标记点中心位置的坐标数据为M4、M5、M6;
5)将步骤2)到步骤4)中所测得的六组坐标数据M1、M2、M3、M4、M5和M6代入相应的六位姿误差分离方程组,即:
其中I为主对角线为1的矩阵,I1、I2为I的x轴方向平移矩阵,R90、R180、R270为I的逆时针90°、180°和270°旋转矩阵,Nx、Ny为栅格板上测量点x和y坐标的标称值,Nxs和Nys以及Nyt和Nxt分别为Nx、Ny由于x轴正负方向平移而形成的矩阵;E为所有元素为1的矩阵;Ax、Ay为测量点的二维平台x和y坐标误差,Gx、Gy为测量点的栅格板x和y坐标标尺误差,Vi、Wi为两理想坐标系原点的x和y坐标偏移,i=1,2,...,6,θi为两理想坐标系的角度偏移,Qix、Qiy分别为第i个位姿中栅格板上测量点的x和y方向的误差,并且有:
6)利用步骤2)到步骤4)所测得的六位姿数据M1、M2、M3、M4、M5、M6以及步骤5)的方程组,得到二维平台误差Ax、Ay;
7)对三坐标测量机二维平台误差进行校正:
其中[Tx,Ty]T为校正之前三坐标测量机所测得的数据,T0则为校正之后的数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106750.2A CN103234496B (zh) | 2013-03-28 | 2013-03-28 | 一种三坐标测量机二维平台误差的高精度校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106750.2A CN103234496B (zh) | 2013-03-28 | 2013-03-28 | 一种三坐标测量机二维平台误差的高精度校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103234496A CN103234496A (zh) | 2013-08-07 |
CN103234496B true CN103234496B (zh) | 2015-09-30 |
Family
ID=48882548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310106750.2A Expired - Fee Related CN103234496B (zh) | 2013-03-28 | 2013-03-28 | 一种三坐标测量机二维平台误差的高精度校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103234496B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103884311B (zh) * | 2014-03-03 | 2017-03-22 | 广东赛因迪科技股份有限公司 | 平面检测机的各个检测器之间的高度差的确定方法 |
CN104374353B (zh) * | 2014-11-20 | 2017-01-25 | 西安电子科技大学 | 一种用于三点法圆度误差分离技术的传感器安装角度选择方法 |
CN105043314B (zh) * | 2015-04-30 | 2017-11-21 | 东莞市神州视觉科技有限公司 | 一种用于锡膏检测精度检验的棱台参数测量方法和*** |
JP2018031754A (ja) * | 2016-08-26 | 2018-03-01 | 株式会社ミツトヨ | 三次元測定装置及び座標補正方法 |
CN109269426B (zh) * | 2018-11-27 | 2024-05-03 | 国家电网有限公司 | 三维位移标定架及三维位移校正系数测量方法 |
CN110057304B (zh) * | 2019-03-01 | 2020-10-20 | 清华大学 | 平面光栅六自由度位移测量***中误差分离与补偿方法 |
CN111982071B (zh) * | 2019-05-24 | 2022-09-27 | Tcl科技集团股份有限公司 | 一种基于tof相机的3d扫描方法及*** |
CN113701798B (zh) * | 2021-08-29 | 2024-01-19 | 孙恺欣 | 一种二维栅格标准板的校准装置和方法 |
CN115876148B (zh) * | 2023-02-07 | 2023-05-19 | 昂坤视觉(北京)科技有限公司 | 一种晶圆平面度测量平台运行误差补偿方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4006835A1 (de) * | 1990-03-05 | 1991-09-12 | Eugen Dr Ing Trapet | Verfahren zur kalibrierung von zweidimensionalen pruefkoerpern mit elimination systematischer fehler des messgeraets |
CN1172246A (zh) * | 1996-07-29 | 1998-02-04 | 清野慧 | 移动物体检测用栅尺及用此栅尺进行移动物体检测的装置 |
CN1991333A (zh) * | 2005-12-30 | 2007-07-04 | 财团法人工业技术研究院 | 零阿贝误差测量***及其方法 |
CN102305608A (zh) * | 2011-05-13 | 2012-01-04 | 哈尔滨工业大学 | 多目标二维交叉运动模拟***误差测量补偿方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3805143B2 (ja) * | 1999-07-14 | 2006-08-02 | 株式会社ミツトヨ | トレンド補正方法 |
-
2013
- 2013-03-28 CN CN201310106750.2A patent/CN103234496B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4006835A1 (de) * | 1990-03-05 | 1991-09-12 | Eugen Dr Ing Trapet | Verfahren zur kalibrierung von zweidimensionalen pruefkoerpern mit elimination systematischer fehler des messgeraets |
CN1172246A (zh) * | 1996-07-29 | 1998-02-04 | 清野慧 | 移动物体检测用栅尺及用此栅尺进行移动物体检测的装置 |
CN1991333A (zh) * | 2005-12-30 | 2007-07-04 | 财团法人工业技术研究院 | 零阿贝误差测量***及其方法 |
CN102305608A (zh) * | 2011-05-13 | 2012-01-04 | 哈尔滨工业大学 | 多目标二维交叉运动模拟***误差测量补偿方法 |
Non-Patent Citations (7)
Title |
---|
Calibration of geometrical systematic error in high-precision spherical surface measurement;daodang wang;《Optics Communications》;20111231;第3878-3884页 * |
Comparative experimental study on absolute measurement of spherical surface with two-sphere method;xi hou;《Optics and Lasers in Engineering》;20111231;第833-840页 * |
IC加工及检测装备超精密工作台自标定技术研究;朱煜 朱立伟;《电子工业专用设备》;20050331;第20-24页 * |
三坐标测量机的误差分析与修正;王琦 陈晓怀;《宇航计测技术》;20060430;第26卷(第2期);第1-3页 * |
三坐标测量机误差测量的新方法;王东升 李书和 张国雄;《天津大学学报》;19990731;第32卷(第4期);第515-517页 * |
新型纳米三坐标测量机误差检定方法的研究;夏瑞雪 陈晓怀;《电子测量与仪器学报》;20100331;第24卷(第3期);第251-256页 * |
超精密级二维工作台的自标定;崔继文 刘雪明 谭久彬;《光学 精密工程》;20120930;第20卷(第9期);第1960-1966页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103234496A (zh) | 2013-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103234496B (zh) | 一种三坐标测量机二维平台误差的高精度校正方法 | |
CN107186548B (zh) | 一种五轴数控机床回转轴几何误差检测方法 | |
CN103389038B (zh) | 激光跟踪仪定目标多站测量的数控机床几何精度检测方法 | |
CN107042528A (zh) | 一种工业机器人的运动学标定***及方法 | |
CN105290968B (zh) | 一种卧式精密磨抛机床工件主轴与工具***的对心方法 | |
CN105486289B (zh) | 一种激光摄影测量***及相机标定方法 | |
Zhang et al. | A general strategy for geometric error identification of multi-axis machine tools based on point measurement | |
CN103499293B (zh) | 一种数控机床的激光跟踪仪虚拟多站式测量方法 | |
CN1667359A (zh) | 超精密工作台自标定方法及装置 | |
CN103026310B (zh) | 实现从机床加工点到安装测量基准点的空间转换方法 | |
CN106363465A (zh) | 多轴数控机床平动轴与旋转轴相互位置关系辨识方法 | |
CN104457688B (zh) | 卫星上批量设备姿态角度矩阵的高精度自动化测量装置 | |
CN102200429A (zh) | 基于激光跟踪复合式测量的数控机床精度检测方法 | |
CN103673976A (zh) | 复合式精度测量坐标系的转换与统一方法及*** | |
CN107101597A (zh) | 一种旋转角测量***的误差标定方法 | |
CN106989670B (zh) | 一种机器人协同的非接触式高精度大型工件跟踪测量方法 | |
CN108801193B (zh) | 一种基于误差与变异规律的三坐标测量机误差测量方法 | |
CN103791878A (zh) | 数控机床几何精度辨识方法 | |
Zha et al. | Four-station laser tracer-based geometric error measurement of rotary table | |
CN106225778A (zh) | 多介质传播中随机波动的视觉定位研究 | |
CN103737433A (zh) | 大尺寸构件复合式精度测量方法 | |
CN105519264B (zh) | 非接触式综合测量大型平板裂缝天线表面三维形貌和缝隙尺寸的方法及其*** | |
CN104977816A (zh) | 基于Compact Particle Swarm Optimization算法的光刻机掩模台微动台的机械参数软测量方法 | |
Hu et al. | A holistic self-calibration algorithm for $ xy $ precision metrology systems | |
CN104634245A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150930 |