CN101870076B - 一种基于性能退化模型的数控机床导轨副寿命预测方法 - Google Patents

一种基于性能退化模型的数控机床导轨副寿命预测方法 Download PDF

Info

Publication number
CN101870076B
CN101870076B CN2010102163927A CN201010216392A CN101870076B CN 101870076 B CN101870076 B CN 101870076B CN 2010102163927 A CN2010102163927 A CN 2010102163927A CN 201010216392 A CN201010216392 A CN 201010216392A CN 101870076 B CN101870076 B CN 101870076B
Authority
CN
China
Prior art keywords
time
life
guideway
neural network
residual life
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
CN2010102163927A
Other languages
English (en)
Other versions
CN101870076A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN2010102163927A priority Critical patent/CN101870076B/zh
Publication of CN101870076A publication Critical patent/CN101870076A/zh
Application granted granted Critical
Publication of CN101870076B publication Critical patent/CN101870076B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

一种基于性能退化模型的数控机床导轨副寿命预测方法,它采集振动信号并进行时频域分析,提取导轨副性能退化敏感特征数据向量并以时间序列方式构成敏感特征矩阵;同时计算导轨副负荷Pi,记录运行时间ti;根据Pi计算额定寿命时间Lhi及导轨副在当前工况下的已经运行的总时间t’和期望剩余寿命LDi,将期望剩余寿命以时间序列方式构成期望剩余寿命向量T;再由双层动态模糊神经网络构成的退化模型对输入的敏感特征矩阵和期望剩余寿命向量T之间的映射关系进行拟合并输出寿命预测结果。该方法考虑数控机床不同工况下导轨副负荷变化对导轨副性能退化的影响,实现导轨副使用过程中剩余寿命的预测,预测精度高,实际使用价值高。

Description

一种基于性能退化模型的数控机床导轨副寿命预测方法
技术领域
本发明属于机电一体化的数控制造技术领域,具体地讲,涉及数控机床导轨副性能退化评估与寿命预测。
背景技术
数控机床的故障严重影响零件的加工质量,给企业带来巨大的经济损失。数控机床的故障包括数控系统故障、电气系统故障、机械系统故障及气动液压系统故障等,其中机械系统的性能退化不可避免,是影响产品质量的关键因素。因此,实时监测机械系统,发现机床运行过程中关键部件如丝杠副、导轨副、主轴组件等的性能退化规律,实时评估部件的剩余寿命及健康状态,有利于建立合理有效的维修计划,减少不必要的停机时间,节省大量的维护费用。
对数控机床而言,由于其加工工况的多样性、负载的动态性、疲劳损伤的不确定性、外界冲击的随机性等,导致机械系统性能退化规律更为复杂,很难建立确切的数学、物理模型来评估机床机械部件的性能退化规律。作为机床关键部件的导轨副,其精度在一定程度上决定了机床的加工精度,其性能随使用时间推移而退化的规律呈现动态变化,目前国内外各机床生产厂家仅限于对导轨副在额定条件下的寿命预测与评估,而导轨副在实际使用过程中大多数情况下均与额定条件不符,因此其寿命预测与评估结果误差大,工程实际应用价值低。
发明内容
本发明的目的就是提供一种基于性能退化模型的数控机床导轨副寿命预测方法,该方法考虑数控机床不同工况下导轨副负荷变化对导轨副性能退化的影响,实现对导轨副使用过程中剩余寿命的评估与预测,预测精度高,工程实际使用价值高。
本发明实现其发明目的所采用的技术方案是,一种基于性能退化模型的数控机床导轨副寿命预测方法,其步骤依次是:
1)振动信号采集:用振动传感器采集导轨副的振动信号,振动信号经信号放大器放大后再通过采集卡送入机床PC端中控服务器进行处理;
2)振动信号时频域分析:机床PC端中控服务器的信号处理模块对采集的振动信号先进行滤波处理,再进行时频域分析得到时频特征数据;
3)时频特征数据提取:对得到的时频特征数据,采用主成分分析法提取出当前工况下对导轨副性能退化敏感的特征数据向量,所述的对导轨副性能退化敏感的特征数据向量是指由主成分分析法分析得到的累计贡献率大于0.9的时频特征分量,将当前工况及以前工况下的敏感特征数据向量以时间序列方式构成敏感特征矩阵,作为退化模型的一个输入量;
4)导轨副负荷计算:由加工参数计算得到当前工况下的切削力Fi
Figure GDA0000080266530000021
式中,下标i表示当前工况,api表示背吃刀量,fi表示进给速度,vi表示切削速度,KF表示修正系数,CF为系数,xF、yF、nF为指数;再将切削力Fi及工作台的重力G分别进行三个方向分解,计算得到二者的合力,作为导轨副负荷Pi,同时记录当前工况下的累计运行时间ti
5)导轨副期望剩余寿命计算:将当前负荷Pi代入导轨副经验寿命计算公式
Figure GDA0000080266530000022
求得当前工况下的导轨副的额定寿命Li,式中,Li表示导轨副额定寿命,C表示额定动负荷,fH表示硬度系数,fT表示温度系数,fC表示接触系数,fw表示负荷系数;利用公式
Figure GDA0000080266530000023
转换为当前工况下的额定寿命时间Lhi,式中,ls表示导轨副行程,nl表示每分钟往返次数;将前一工况下导轨副已运行的时间转换为当前工况下导轨副已运行时间t’i-1=ti-1×Pi-1/Pi,则导轨副在当前工况下的已经运行的总时间t’为:
Figure GDA0000080266530000024
其中k为工况的序号,再将当前工况下的额定寿命时间Lhi减去当前工况下的已运行总时间t’得到当前工况下的期望剩余寿命LDi,即LDi=Lhi-t′;将当前工况及以前工况下的期望剩余寿命以时间序列方式构成期望剩余寿命向量T,作为退化模型的另一输入量;
6)导轨副性能退化评估与剩余寿命预测:退化模型由两个动态模糊神经网络构成,第一个动态模糊神经网络用来在线训练拟合敏感特征数据与期望剩余寿命之间的映射关系并输出寿命预测结果,第二个动态模糊神经网络用于在线学习和权值更新。将步骤3)得到的敏感特征矩阵与步骤5)得到的期望剩余寿命向量同时输入退化模型,退化模型自动选择第一个动态模糊神经网络进行训练以拟合输入的敏感特征值与期望剩余寿命之间的映射关系,第一个动态模糊神经网络训练好后输出预测剩余寿命结果,此时第二个动态模糊神经网络处于备用状态;
若预测剩余寿命结果与期望剩余寿命误差超过设定阈值,则将第一个动态模糊神经网络的权值复制给第二个动态模糊神经网络,第一个动态模糊神经网络继续预测并输出预测剩余寿命结果,而第二个动态模糊神经网络开始学习,快速训练网络得到新的权值,并把新的权值复制给第一个动态模糊神经网络,第一个动态模糊神经网络将得到的新权值覆盖原来的权值再进行寿命预测,直到预测剩余寿命结果与期望剩余寿命误差小于误差阈值,第二个动态模糊神经网络即停止学习进入备用状态,第一个动态模糊神经网络迅速恢复正常预测并输出预测剩余寿命结果。
与现有技术相比,本发明的有益效果是:
本发明将导轨副的振动状态信息特征与时间变量结合起来,考虑不同工况下负荷的变化对导轨副寿命的影响,与机床实际使用情况更符合,使其预测结果更准确;利用动态模糊神经网络构建的退化模型,能够更好的拟合振动信号特征与导轨副寿命之间的复杂映射关系,实现对导轨副实际使用过程中剩余寿命进行更精确的实时评估与预测;双层动态模糊神经网络能够实现边学习边预测功能,使系统抗干扰能力强,鲁棒性好。本发明有利于数控机床使用企业建立合理有效的维修计划,减少不必要的停机时间,节省大量的维护费用,工程实际使用价值高。
具体实施方式
实施例
一种基于性能退化模型的数控机床导轨副寿命预测方法,其步骤依次是:
1)振动信号采集:用振动传感器采集导轨副的振动信号,振动信号经信号放大器放大后再通过采集卡送入机床PC端中控服务器进行处理。
振动传感器和信号放大器可选用各种现有的传感器,如选用丹麦B&K公司的BK4321三向加速度振动传感器、BK2365振动信号放大器。振动传感器具体安装在导轨副的滑块外侧,其数目根据数控机床的导轨副的滑块数而定,采集导轨副三向振动信息,每个传感器输出三通道振动信号。采集卡将所有传感器全部通道振动信号传输给机床PC端中控服务器。
2)振动信号时频域分析:机床PC中控服务器对采集的振动信号先进行滤波处理,再进行时频域分析,得到时频特征数据。
本例的2)步的具体操作是:机床PC中控服务器将所有通道振动信号进行小波包分解,确定最优小波包基,对小波包分解系数的阈值进行量化,再对信号进行小波包重构,得到振动信号的时域图,对重构后的信号进行傅里叶变换,得到信号频域图,根据时频域信息提取包括均值、均方根值、敏感维数、盒维数、倒谱峰值、均方根率、频率标准差、波形参数、峰值因子、脉冲因子、峭度系数等时频特征分量的时频特征数据,对每个传感器的三通道信号的时频特征值进行比较,选取特征值大的通道的时频特征值作为该传感器的时频特征数据;再对上述比较所得的所有传感器的时频特征数据进行平均作为时频域分析得到的时频特征数据。
3)时频特征数据提取:对得到的时频特征数据,采用主成分分析法提取出当前工况下对导轨副性能退化敏感的特征数据向量,将当前工况及以前工况下的敏感特征数据向量以时间序列方式构成敏感特征矩阵,作为退化模型的一个输入量。
本例中,将累计贡献率大于0.9的时频特征分量作为对导轨副性能退化敏感的特征数据向量。
4)导轨副负荷计算:由加工参数计算得到当前工况下的切削力Fi式中,下标i表示当前工况,api表示背吃刀量,fi表示进给速度,vi表示切削速度,KF表示修正系数,CF为系数,xF、yF、nF为指数;再将切削力Fi及工作台的重力G分别进行三个方向分解,计算得到二者的合力,作为导轨副负荷Pi,同时记录当前工况下的累计运行时间ti
本例中,当前工况下切削力Fi及丝杠副负荷Pi的最终值,也即以下第(5)步使用的Fi及Pi值:为在当前工况下各采样时刻计算得出的瞬时值的平均值。
同时,本例中记录当前工况下的累计运行时间ti的具体做法是:
机床运行时,PC中控服务器的计时器开始计时,当PC中控服务器得到的导轨副负荷Pi变化率大于设置的阈值(如5%)时,PC中控服务器认定新的工况出现,计时器归零,按下一个工况重新开始计时(此时PC中控服务器计算出的负荷也为下一个工况的负荷),否则计时器连续计时。
5)导轨副期望剩余寿命计算:将当前负荷Pi代入导轨副经验寿命计算公式求得当前工况下的导轨副的额定寿命Li,式中,Li表示导轨副额定寿命(运行距离),C表示额定动负荷,fH表示硬度系数,fT表示温度系数,fC表示接触系数,fw表示负荷系数;利用公式
Figure GDA0000080266530000053
转换为当前工况下的额定寿命时间Lhi,式中,ls表示导轨副行程,nl表示每分钟往返次数;将前一工况下导轨副已运行的时间转换为当前工况下导轨副已运行时间t’i-1=ti-1×Pi-1/Pi,则导轨副在当前工况下的已经运行的总时间t’为:
Figure GDA0000080266530000061
其中k为工况的序号,再将当前工况下的额定寿命时间Lhi减去当前工况下的已运行总时间t’得到当前工况下的期望剩余寿命LDi,即LDi=Lhi-t′;将当前工况及以前工况下的期望剩余寿命以时间序列方式构成期望剩余寿命向量T,作为退化模型的另一输入量。
6)导轨副性能退化评估与剩余寿命预测:退化模型由两个动态模糊神经网络(Dynamic Fuzzy Neural Network,DFNN)构成,第一个动态模糊神经网络用来在线训练拟合敏感特征数据与期望剩余寿命之间的映射关系并输出寿命预测结果,第二个动态模糊神经网络用于在线学习和权值更新;即,步骤3)得到的敏感特征矩阵与步骤5)得到的期望剩余寿命向量同时输入退化模型,退化模型自动选择第一个动态模糊神经网络进行训练以拟合输入的敏感特征值与期望剩余寿命之间的映射关系,第一个动态模糊神经网络训练好后输出预测剩余寿命结果,此时第二个动态模糊神经网络处于备用状态;
若预测剩余寿命结果与期望剩余寿命误差超过设定阈值,则将第一个动态模糊神经网络的权值复制给第二个动态模糊神经网络,第一个动态模糊神经网络继续预测并输出预测剩余寿命结果,而第二个动态模糊神经网络开始学习,快速训练网络得到新的权值,并把新的权值复制给第一个动态模糊神经网络,第一个动态模糊神经网络将得到的新权值覆盖原来的权值再进行寿命预测,直到预测剩余寿命结果与期望剩余寿命误差小于误差阈值,第二个动态模糊神经网络即停止学习进入备用状态,第一个动态模糊神经网络迅速恢复正常预测并输出预测剩余寿命结果。
本领域的普通技术人员显然清楚并且理解,本发明方法所举的以上实施例仅用于说明本发明方法,而并不用于限制本发明方法。虽然通过实施例有效描述了本发明,本发明存在许多变化而不脱离本发明的精神。在不背离本发明方法的精神及其实质的情况下,本领域技术人员当可根据本发明方法做出各种相应的改变或变形,但这些相应的改变或变形均属于本发明方法要求的保护范围。

Claims (1)

1.一种基于性能退化模型的数控机床导轨副寿命预测方法,其步骤依次是:
1)振动信号采集:用振动传感器采集导轨副的振动信号,振动信号经信号放大器放大后再通过采集卡送入机床PC端中控服务器进行处理;
2)振动信号时频域分析:机床PC端中控服务器的信号处理模块对采集的振动信号先进行滤波处理,再进行时频域分析得到时频特征数据;
3)时频特征数据提取:对得到的时频特征数据,采用主成分分析法提取出当前工况下对导轨副性能退化敏感的特征数据向量,所述的对导轨副性能退化敏感的特征数据向量是指由主成分分析法分析得到的累计贡献率大于0.9的时频特征分量,将当前工况及以前工况下的敏感特征数据向量以时间序列方式构成敏感特征矩阵,作为退化模型的一个输入量;
4)导轨副负荷计算:由加工参数计算得到当前工况下的切削力Fi
Figure FDA0000080266520000011
式中,下标i表示当前工况,api表示背吃刀量,fi表示进给速度,vi表示切削速度,KF表示修正系数,CF为系数,xF、yF、nF为指数;再将切削力Fi及工作台的重力G分别进行三个方向分解,计算得到二者的合力,作为导轨副负荷Pi,同时记录当前工况下的累计运行时间ti
5)导轨副期望剩余寿命计算:将当前负荷Pi代入导轨副经验寿命计算公式求得当前工况下的导轨副的额定寿命Li,式中,Li表示导轨副额定寿命,C表示额定动负荷,fH表示硬度系数,fT表示温度系数,fC表示接触系数,fw表示负荷系数;利用公式
Figure FDA0000080266520000013
转换为当前工况下的额定寿命时间Lhi,式中,ls表示导轨副行程,nl表示每分钟往返次数;将前一工况下导轨副已运行的时间转换为当前工况下导轨副已运行时间t’i-1=ti-1×Pi-1/Pi,则导轨副在当前工况下的已经运行的总时间t’为:
Figure FDA0000080266520000014
其中k为工况的序号,再将当前工况下的额定寿命时间Lhi减去当前工况下的已运行总时间t’得到当前工况下的期望剩余寿命LDi,即LDi=Lhi-t′;将当前工况及以前工况下的期望剩余寿命以时间序列方式构成期望剩余寿命向量T,作为退化模型的另一输入量;
6)导轨副性能退化评估与剩余寿命预测:退化模型由两个动态模糊神经网络构成,第一个动态模糊神经网络用来在线训练拟合敏感特征数据与期望剩余寿命之间的映射关系并输出寿命预测结果,第二个动态模糊神经网络用于在线学习和权值更新;将步骤3)得到的敏感特征矩阵与步骤5)得到的期望剩余寿命向量同时输入退化模型,退化模型自动选择第一个动态模糊神经网络进行训练以拟合输入的敏感特征值与期望剩余寿命之间的映射关系,第一个动态模糊神经网络训练好后输出预测剩余寿命结果,此时第二个动态模糊神经网络处于备用状态;
若预测剩余寿命结果与期望剩余寿命误差超过设定阈值,则将第一个动态模糊神经网络的权值复制给第二个动态模糊神经网络,第一个动态模糊神经网络继续预测并输出预测剩余寿命结果,而第二个动态模糊神经网络开始学习,快速训练网络得到新的权值,并把新的权值复制给第一个动态模糊神经网络,第一个动态模糊神经网络将得到的新权值覆盖原来的权值再进行寿命预测,直到预测剩余寿命结果与期望剩余寿命误差小于误差阈值,第二个动态模糊神经网络即停止学习进入备用状态,第一个动态模糊神经网络迅速恢复正常预测并输出预测剩余寿命结果。
CN2010102163927A 2010-07-02 2010-07-02 一种基于性能退化模型的数控机床导轨副寿命预测方法 Expired - Fee Related CN101870076B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102163927A CN101870076B (zh) 2010-07-02 2010-07-02 一种基于性能退化模型的数控机床导轨副寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102163927A CN101870076B (zh) 2010-07-02 2010-07-02 一种基于性能退化模型的数控机床导轨副寿命预测方法

Publications (2)

Publication Number Publication Date
CN101870076A CN101870076A (zh) 2010-10-27
CN101870076B true CN101870076B (zh) 2012-03-21

Family

ID=42995256

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102163927A Expired - Fee Related CN101870076B (zh) 2010-07-02 2010-07-02 一种基于性能退化模型的数控机床导轨副寿命预测方法

Country Status (1)

Country Link
CN (1) CN101870076B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622496B (zh) * 2011-01-26 2016-07-06 中国科学院大气物理研究所 一种嵌入模糊集状态的自适应多步预报方法和系统
CN102262700B (zh) * 2011-08-01 2013-01-30 北京航空航天大学 基于小波分析的退化数据预处理的产品寿命预测方法
CN102689231B (zh) * 2012-06-12 2013-12-11 西南交通大学 一种丝杠副、导轨副性能退化的图示化分析方法
CN102789545B (zh) * 2012-07-12 2015-08-19 哈尔滨工业大学 基于退化模型匹配的涡轮发动机剩余寿命的预测方法
CN103679280B (zh) * 2012-09-26 2016-12-21 中国人民解放军第二炮兵工程大学 一种性能缓变退化的设备最优维护方法
CN103488881B (zh) * 2013-09-06 2017-05-17 中国人民解放军第二炮兵工程大学 一种不确定退化测量数据下的设备剩余寿命预测方法
CN103793752B (zh) * 2013-09-13 2018-02-06 中国人民解放军第二炮兵工程大学 一种基于退化建模的设备失效次数预测方法
CN106842922B (zh) * 2017-01-14 2020-07-17 合肥工业大学 一种数控加工误差优化方法
CN107703752A (zh) * 2017-10-22 2018-02-16 成都具鑫机械设备有限公司 一种基于模糊神经网络的数控机床检测方法
CN108629073B (zh) * 2018-03-14 2019-05-03 山东科技大学 一种多模式的退化过程建模及剩余寿命预测方法
CN109101753B (zh) * 2018-08-31 2023-06-27 周建全 一种改善机械设备检修系统复杂性的方法
CN109211564B (zh) * 2018-10-24 2020-05-01 哈工大机器人(山东)智能装备研究院 一种用于滚珠丝杠副健康评估的自适应阈值检测方法
CN109255201B (zh) * 2018-10-24 2023-07-14 哈工大机器人(山东)智能装备研究院 一种基于som-mqe的滚珠丝杠副健康评估方法
CN109739183B (zh) * 2019-01-14 2020-05-22 上海赛卡精密机械有限公司 一种基于多传感器的大型数控机床故障监测系统
CN111813046A (zh) * 2020-08-11 2020-10-23 深圳市玄羽科技有限公司 数控机床撞机的监控方法及系统、电子设备及存储介质
CN113870458B (zh) * 2021-10-26 2023-05-30 中冶赛迪工程技术股份有限公司 结晶器振动装置振动性能退化数据构建方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3151655B2 (ja) * 1996-02-19 2001-04-03 オークマ株式会社 工作機械の熱変位推定方法
JP4876528B2 (ja) * 2005-10-25 2012-02-15 株式会社ジェイテクト 工作機械における部品寿命管理システム
DE102006055460A1 (de) * 2006-11-24 2008-05-29 Robert Kuchel Verfahren zur Schmierung von gesteuert verfahrbaren Werkzeugen sowie Schmiermittelspender
CN101373495B (zh) * 2007-08-24 2010-09-29 西门子公司 使用寿命终点判定及当前历史使用寿命估计的方法和系统
CN101960283B (zh) * 2008-03-04 2014-02-19 新日铁住金化学株式会社 层叠体的弯曲寿命预测方法和层叠体的弯曲寿命预测装置
CN101354311A (zh) * 2008-09-05 2009-01-28 重庆大学 汽车后桥寿命预测系统

Also Published As

Publication number Publication date
CN101870076A (zh) 2010-10-27

Similar Documents

Publication Publication Date Title
CN101870076B (zh) 一种基于性能退化模型的数控机床导轨副寿命预测方法
CN101870075B (zh) 一种基于性能退化模型的数控机床丝杠副寿命预测方法
CN109933004B (zh) 基于边缘计算和云协同的机床故障诊断与预测方法及系统
CN109396953B (zh) 基于信号融合的机床工作状态智能辨识系统
CN108038553B (zh) 轧机设备状态在线监测与诊断系统及监测诊断方法
CN101520651B (zh) 一种基于隐马尔科夫链的数控装备可靠性分析方法
CN110355608B (zh) 基于自注意力机制与深度学习的刀具磨损量预测方法
CN104808585B (zh) 一种机床健康状态快速检查方法
CN103034170B (zh) 一种基于区间的数控机床性能预测方法
CN104390657B (zh) 一种发电机组运行参数测量传感器故障诊断方法及系统
CN105109094B (zh) 基于计算机的高精度伺服压装方法
CN103631681A (zh) 一种在线修复风电场异常数据的方法
CN102848266B (zh) 一种机床主轴精度预测方法
CN104808587A (zh) 一种基于机加工设备运行状态的稼动率统计方法
CN109376401A (zh) 一种自适应多源信息优选与融合的机械剩余寿命预测方法
CN108520125A (zh) 一种预测刀具磨损状态的方法及系统
CN112668105B (zh) 一种基于sae与马氏距离的直升机传动轴异常判定方法
CN111458629B (zh) 一种高压开关机械故障的反演方法及装置
CN103203670B (zh) 基于最大信息熵与方向散度的磨削颤振预测方法
CN114135477A (zh) 一种机泵设备状态监测动态阈值预警方法
CN116893643A (zh) 一种基于数据分析的智能机器人行驶轨迹安全管控系统
CN109598309B (zh) 一种金属包装冲压机的检测系统及其监测方法
CN102879192B (zh) 一种直线导轨副精度预测方法
CN101529347A (zh) 识别有延时的调节对象的方法和装置、调节装置和计算机程序产品
CN116703254A (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

Granted publication date: 20120321

Termination date: 20140702

EXPY Termination of patent right or utility model