CN103528599B - 基于非均匀光纤布拉格光栅分布式传感的光谱识别方法 - Google Patents
基于非均匀光纤布拉格光栅分布式传感的光谱识别方法 Download PDFInfo
- Publication number
- CN103528599B CN103528599B CN201310455613.XA CN201310455613A CN103528599B CN 103528599 B CN103528599 B CN 103528599B CN 201310455613 A CN201310455613 A CN 201310455613A CN 103528599 B CN103528599 B CN 103528599B
- Authority
- CN
- China
- Prior art keywords
- epsiv
- grating
- strain
- under
- extraneous
- 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
Landscapes
- Optical Transform (AREA)
Abstract
基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,涉及基于非均匀光纤布拉格光栅分布式传感的光谱识别技术领域。解决了现有光谱识别方法在进行光谱识别时,当光纤布拉格光栅长度与被测物体的尺寸相当,同时被测物体体内存在非均匀的传感参数的作用场时,导致光纤布拉格光栅的反射谱形状发生改变,进而影响光谱识别效率的问题。基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,它包括以下步骤:计算非均匀应变场作用下的光栅折射率分布;建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型;传输矩阵模型得出非均匀外界应变场作用下直角三角谱光纤布拉格的光谱变化。本发明适用于非均匀光纤布拉格光栅的光谱识别。
Description
技术领域
本发明涉及光谱识别技术领域,具体涉及基于非均匀光纤布拉格光栅分布式传感的光谱识别技术领域。
背景技术
光纤布拉格光栅(FBG,即,fiberBragggrating)作为传感元件广泛地应用于传感系统中,光纤布拉格光栅传感器是一种基于波长调制的新型光纤传感器。通常的传感应用往往都是准分布式传感。所谓的准分布式传感就是在光栅长度远远小于被测物尺寸的时候,可以把每个光栅看成是一个测量点,为了分析外场对被测物的影响,一般在被测物上设置多个光栅,然后通过同时测量所有测量点的外场来实现对被测物所受到的外场进行分析的一种传感方式。在这种情况下,可以认为沿光栅轴向施加的传感参数是常量,光栅的光谱没有发生形变,而仅仅是光谱出现了整体移动。因此,传感参数的作用就可以通过光谱的峰值反射率对应的布拉格波长的偏移来反映,进而实现对传感参数的传感。
与FBG在传感系统中的准分布式传感不同,当光栅长度与被测物的尺寸相当,同时被测物体内又存在一个非均匀的传感参数的作用场时,那么FBG的反射谱形状会发生改变,原有的光栅布拉格波长同传感参数之间的线性对应关系也就不成立了,所以此时不能再通过测量布拉格波长的偏移量来进行传感参数的测量,而是应该对FBG的反射谱进行更加深入的理解和分析。通常情况下,为了实现各种不均匀传感参量在传感上的应用,采用的基本途径是对非均匀外界作用场下的FBG进行合成,通过光谱合成可以得到光栅上传感参数的分布情况,人们把这个过程称作光栅的分布式传感应用。所有这些基于光栅的传感系统其空间分辨率都受限于光栅的长度,并且光栅上非均匀应力场也影响应力的分辨率,特别是用于材料中有破坏或者应力集中的地方,光栅的光谱形状常常会出现有意义畸变。这种畸变往往是由于光栅的光纤轴上存在应力梯度、非均匀或者不连续的应力场等情况而引起的,因而可以通过实验中测量得到的畸变光谱来判断FBG所受到的外界力场作用,这一点在材料健康监测上是有潜在应用的。
发明内容
本发明为了解决现有光谱识别方法在进行光谱识别时,当光纤布拉格光栅长度与被测物体的尺寸相当,同时被测物体体内存在非均匀的传感参数的作用场时,导致光纤布拉格光栅的反射谱形状发生改变,进而影响光谱识别效率的问题,提出了基于非均匀光纤布拉格光栅分布式传感的光谱识别方法。
基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,它包括以下步骤:
步骤一、计算非均匀应变场作用下的光栅折射率分布;
步骤二、建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型;
步骤三、根据步骤二得到的传输矩阵模型得出非均匀外界应变场作用下直角三角谱光纤布拉格的光谱变化。
建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型,它是由以下步骤实现的:
步骤二一、对没有加入外界应变之前光栅的有效周期进行计算;
步骤二二、通过步骤二一得到的没有加入外界应变之前光栅的有效周期和外界应变引起的光栅周期变化与外界应变的关系,计算出非均匀力场作用下的相位变化;
步骤二三、计算在非均匀外界应变场作用下,光纤中前向传播光场包络的耦合模和后向传播光场包络的耦合模;
步骤二四、计算出外界应变作用下,光传输到光栅位置z处的波数失谐量与外界应变之间的关系;
步骤二五、计算出微应变下光栅耦合系数与非均匀应变场作用下的交流折射率、非均匀应变场作用下的折射率交流分量的变化量和光栅的相位之间的关系;
步骤二六、计算出在外界应变作用下,在坐标位置z处与z+△lε(z)处的第d段光场的传输矩阵,△lε(z)为非均匀应变场作用下光栅第d段光场长度,即△lε(z)=△l(1+ε(z)),△l为没有非均匀应变场作用时光栅第d段光场长度,△l=L/N,N为没有外界应变作用之前光栅长度L平均分成的层数,ε(z)为外界应变;
步骤二七、通过步骤二六得到在外界应变作用下总的传输矩阵;
步骤二八、通过步骤二七得到的在外界应变作用下总的传输矩阵计算出微应变下光栅光谱的反射系数。
有益效果:本发明提出的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法通过传输矩阵法对非均匀光纤布拉格光栅在外界应变场作用下的传输矩阵进行推导,能够使光纤布拉格光栅长度与被测物体的尺寸相当,同时被测物体体内存在非均匀的传感参数的作用场时,光线布拉格光栅的反射谱形状不会发生改变,同时使本发明的光谱识别方法更加简单和便于计算,进而使光谱识别效率提高了5%以上。
附图说明
图1为基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的流程图;
图2为建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型的流程图;
图3为应变场作用下的传输矩阵法示意图;
图4为长度为12mm的直角三角谱光纤布拉格光栅的反射谱图;
图5为长度为12mm的直角三角谱光纤布拉格光栅的交流折射率调制;
图6为长度为12mm的直角三角谱光纤布拉格光栅的相位调制;
图7为均匀应变场作用下直角三角谱光纤布拉格光栅的反射谱;
图8为线性分布的非均匀应变作用下直角三角谱光纤布拉格光栅的反射谱;
图9为二次曲线分布的非均匀应变场作用下直角三角谱光纤布拉格光栅的反射谱;
图10为存在突变的非均匀应变场作用下直角三角谱光纤布拉格光栅的反射谱。
具体实施方式
具体实施方式一、结合图1说明本具体实施方式,基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,它包括以下步骤:
步骤一、计算非均匀应变场作用下的光栅折射率分布;
步骤二、建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型;
步骤三、根据步骤二得到的传输矩阵模型得出非均匀外界应变场作用下直角三角谱光纤布拉格的光谱变化。
具体实施方式二、本具体实施方式与具体实施方式一所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤一中所述的计算非均匀应变场作用下的光栅折射率分布由以下公式实现:
其中,nε(z)为非均匀应变场作用下的光栅折射率分布,neff为没有紫外曝光前光纤纤芯的有效折射率,△ndc(z)为非均匀应变场作用下光栅的直流折射率,△nac(z)为非均匀应变场作用下光栅的交流折射率,△nac,ε(z)为非均匀应变场作用下光栅的折射率交流分量的变化量;△ndc,ε(z)为非均匀应变场作用下光栅的折射率直流分量的变化量;z为非均匀应变场作用前的光栅坐标位置,为非均匀应变场作用下的相位变化,为非均匀应变场作用前光栅的相位,Λ为非均匀应变场作用前光栅周期。
具体实施方式三、结合图2和图3说明本具体实施方式,本具体实施方式与具体实施方式一所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二中所述的建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型,它是由以下步骤实现的:
步骤二一、对没有加入外界应变之前光栅的有效周期进行计算;
步骤二二、通过步骤二一得到的没有加入外界应变之前光栅的有效周期和外界应变引起的光栅周期变化与外界应变的关系,计算出非均匀力场作用下的相位变化;
步骤二三、计算在非均匀外界应变场作用下,光纤中前向传播光场包络的耦合模和后向传播光场包络的耦合模;
步骤二四、计算出外界应变作用下,光传输到光栅位置z处的波数失谐量与外界应变之间的关系;
步骤二五、计算出微应变下光栅耦合系数与非均匀应变场作用下的交流折射率、非均匀应变场作用下的折射率交流分量的变化量和光栅的相位之间的关系;
步骤二六、计算出在外界应变作用下,在坐标位置z处与z+△lε(z)处的第d段光场的传输矩阵,△lε(z)为非均匀应变场作用下光栅第d段光场长度,即△lε(z)=△l(1+ε(z)),△l为没有非均匀应变场作用时光栅第d段光场长度,△l=L/N,N为没有外界应变作用之前光栅长度L平均分成的层数,ε(z)为外界应变;
步骤二七、通过步骤二六得到在外界应变作用下总的传输矩阵;
步骤二八、通过步骤二七得到的在外界应变作用下总的传输矩阵计算出微应变下光栅光谱的反射系数。
具体实施方式四、本具体实施方式与具体实施方式三所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二一中所述的没有加入外界应变之前光栅的有效周期的表达式为Λeff(z)为没有加入外界应变之前光栅的有效周期。
具体实施方式五、本具体实施方式与具体实施方式三所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二二中所述的非均匀力场作用下的相位变化的表达式为 为非均匀力场作用下的相位变化。
具体实施方式六、本具体实施方式与具体实施方式三所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二三中所述的在非均匀外界应变场作用下,光纤中前向传播光场包络的耦合模的表达式为光纤中后向传播光场包络的耦合模的表达式为为外界应变作用下,光传输到光栅位置z处的波数失谐量,qε(z)为外界应变作用下,光传输到光栅位置z处时的耦合系数,为外界应变作用下光栅第d段光场的耦合系数。
具体实施方式七、本具体实施方式与具体实施方式三所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二四中所述的光传输到光栅位置z处的波数失谐量与外界应变之间的关系的表达式为 λ为布拉格波长,为外界应变作用下传输矩阵中第z段的中心波长,且
本实施方式中,△ndc,ε(z)是依据非均匀应变场下光栅每段的折射率变化△nε(z)与外界应变ε(z)的关系,以及△ndc,ε(z)的定义计算出来的,即△nε(z)=-peneffε(z),pe为有效弹光系数,通常的取值为0.22,因此,当ε(z)为一个常数时,△nac,ε(z)=0且△ndc,ε(z)=△nε(z),但是,当ε(z)不是一个常数时,光栅的传输矩阵模型中每个离散段上不同微应变ε(z)引起的折射率变化△nε(z)是不同的,再根据△ndc,ε(z)的定义从而计算出由此得出△nac,ε(z)=△nε(z)-△ndc,ε(z)。
具体实施方式八、本具体实施方式与具体实施方式三所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二五中所述的微应变下光栅耦合系数与非均匀应变场作用下的交流折射率、非均匀应变场作用下的折射率交流分量的变化量和光栅的相位之间的关系的表达式为
具体实施方式九、本具体实施方式与具体实施方式三所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二六中所述的在外界应变作用下,在坐标位置z处与z+△l处的第d段光场的传输矩阵的表达式为:
具体实施方式十、本具体实施方式与具体实施方式三所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法的区别在于,步骤二七中所述的在外界应变作用下总的传输矩阵的表达式为
由于为2×2矩阵,且因此,Tε也为2×2矩阵,即
当Tε已知,则微应变下光栅光谱的反射系数可以通过公式rε(δ)=-T21/T22计算出。
直角三角谱光纤布拉格光栅是典型的非均匀光栅,沿着光轴上的折射率分布既存在交流折射率调制△nac,ε(z)也存在着相位调制图3给出了长度为12mm的直角三角谱光纤布拉格光栅的反射谱图,图4给出了对应图3的直角三角谱光纤布拉格光栅的交流折射率调制,图5给出了对应图3的直角三角谱光纤布拉格光栅的相位调制。
下面给出在四种不同的应变场作用下的光栅光谱的变化情况:
1、均匀应变场
理想的均匀应变场的分布表达式:ε(z)=400(με)0≤z≤12(mm);
2、线性分布的非均匀应变场
理想的线性分布的非均匀应变场的分布表达式:ε(z)=5z+100(με)0≤z≤12(mm);
3、二次曲线分布的非均匀应变场
理想的二次曲线分布的非均匀应变场的分布表达式:ε(z)=2z2+10z+100(με)0≤z≤12(mm);
4、存在突变的非均匀应变场
理想的存在突变的非均匀应变场的分布表达式:
图7-图10分别为均匀应变场作用下直角三角谱光纤布拉格光栅的反射谱、线性分布的非均匀应变作用下直角三角谱光纤布拉格光栅的反射谱、二次曲线分布的非均匀应变场作用下直角三角谱光纤布拉格光栅的反射谱和存在突变的非均匀应变场作用下直角三角谱光纤布拉格光栅的反射谱;从这四幅图中可以看出,在均匀应变场作用下,反射谱的变化与均匀FBG的情况相同,即反射谱的形状不变,仅仅光栅布拉中心波长与反射谱同步移动;在非均匀应变场的作用下,反射谱出现了整体偏移的同时,其形状也发生改变,外界应变越大越复杂,反射谱形变越明显,并且带宽增加也就越大。
Claims (8)
1.基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,它包括以下步骤:
步骤一、计算非均匀应变场作用下的光栅折射率分布;
步骤二、建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型;
所述的建立非均匀外界应变场作用下光栅光谱识别的传输矩阵模型,它是由以下步骤实现的:
步骤二一、对没有加入外界应变之前光栅的有效周期进行计算;
步骤二一、步骤二二、通过步骤二一得到的没有加入外界应变之前光栅的有效周期和外界应变引起的光栅周期变化与外界应变的关系,计算出非均匀力场作用下的相位变化;
步骤二三、计算在非均匀外界应变场作用下,光纤中前向传播光场包络的耦合模和后向传播光场包络的耦合模;
步骤二四、计算出外界应变作用下,光传输到光栅位置z处的波数失谐量与外界应变之间的关系;
步骤二五、计算出微应变下光栅耦合系数与非均匀应变场作用下的交流折射率、非均匀应变场作用下的折射率交流分量的变化量和光栅的相位之间的关系;
步骤二六、计算出在外界应变作用下,在坐标位置z处与z+△lε(z)处的第d段光场的传输矩阵,△lε(z)为非均匀应变场作用下光栅第d段光场长度,即△lε(z)=△l(1+ε(z)),△l为没有非均匀应变场作用时光栅第d段光场长度,△l=L/N,N为没有外界应变作用之前光栅长度L平均分成的层数,ε(z)为外界应变;
步骤二七、通过步骤二六得到在外界应变作用下总的传输矩阵;
步骤二八、通过步骤二七得到的在外界应变作用下总的传输矩阵计算出微应变下光栅光谱的反射系数;
步骤三、根据步骤二得到的传输矩阵模型得出非均匀外界应变场作用下直角三角谱光纤布拉格的光谱变化;
其特征在于,步骤一中所述的计算非均匀应变场作用下的光栅折射率分布由以下公式实现:
其中,nε(z)为非均匀应变场作用下的光栅折射率分布,neff为紫外曝光前光纤纤芯的有效折射率,△ndc(z)为非均匀应变场作用下光栅的直流折射率,△nac(z)为非均匀应变场作用下光栅的交流折射率,△nac,ε(z)为非均匀应变场作用下光栅的折射率交流分量的变化量;△ndc,ε(z)为非均匀应变场作用下光栅的折射率直流分量的变化量;z为非均匀应变场作用前的光栅坐标位置,为非均匀应变场作用下的相位变化,为非均匀应变场作用前光栅的相位,Λ为非均匀应变场作用前光栅周期。
2.根据权利要求1所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,其特征在于,步骤二一中所述的没有加入外界应变之前光栅的有效周期的表达式为Λeff(z)为没有加入外界应变之前光栅的有效周期。
3.根据权利要求1所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,其特征在于,步骤二二中所述的非均匀力场作用下的相位变化的表达式为 为非均匀力场作用下的相位变化。
4.根据权利要求1所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,其特征在于,步骤二三中所述的在非均匀外界应变场作用下,光纤中前向传播光场包络的耦合模的表达式为光纤中后向传播光场包络的耦合模的表达式为δε(z)为外界应变作用下,光传输到光栅位置z处的波数失谐量,qε(z)为外界应变作用下,光传输到光栅位置z处时的耦合系数,为外界应变作用下光栅第d段光场的耦合系数。
5.根据权利要求1所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,其特征在于,步骤二四中所述的光传输到光栅位置z处的波数失谐量与外界应变之间的关系的表达式为 λ为布拉格波长,为外界应变作用下传输矩阵中第z段的中心波长,且
6.根据权利要求1所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,其特征在于,步骤二五中所述的微应变下光栅耦合系数与非均匀应变场作用下的交流折射率、非均匀应变场作用下的折射率交流分量的变化量和光栅的相位之间的关系的表达式为
7.根据权利要求1所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,其特征在于,步骤二六中所述的在外界应变作用下,在坐标位置z处与z+△l处的第d段光场的传输矩阵的表达式为:
8.根据权利要求1所述的基于非均匀光纤布拉格光栅分布式传感的光谱识别方法,其特征在于,步骤二七中所述的在外界应变作用下总的传输矩阵的表达式为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310455613.XA CN103528599B (zh) | 2013-09-29 | 2013-09-29 | 基于非均匀光纤布拉格光栅分布式传感的光谱识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310455613.XA CN103528599B (zh) | 2013-09-29 | 2013-09-29 | 基于非均匀光纤布拉格光栅分布式传感的光谱识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103528599A CN103528599A (zh) | 2014-01-22 |
CN103528599B true CN103528599B (zh) | 2016-06-08 |
Family
ID=49930798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310455613.XA Expired - Fee Related CN103528599B (zh) | 2013-09-29 | 2013-09-29 | 基于非均匀光纤布拉格光栅分布式传感的光谱识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103528599B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109000694B (zh) * | 2018-10-19 | 2020-11-10 | 武汉烽理光电技术有限公司 | 光栅制备在线监测方法及系统 |
CN109521518B (zh) * | 2018-12-14 | 2020-06-02 | 南京师范大学 | 基于直流调制的多通道光纤光栅滤波器及其制作方法 |
CN114923420B (zh) * | 2022-05-18 | 2022-12-27 | 中南大学 | 基于光纤布拉格光栅的裂纹诊断方法、系统及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7664347B2 (en) * | 2006-06-07 | 2010-02-16 | Baker Hughes Incorporated | Multi-core optical fiber sensor |
CN102032954B (zh) * | 2010-12-28 | 2012-07-11 | 上海大学 | 光纤布拉格光栅快速光谱反演解调温度场方法 |
-
2013
- 2013-09-29 CN CN201310455613.XA patent/CN103528599B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103528599A (zh) | 2014-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103712575B (zh) | 光学的弯曲曲率测试方法及传感器 | |
US7520176B1 (en) | Method for real-time structure shape-sensing | |
CN103528599B (zh) | 基于非均匀光纤布拉格光栅分布式传感的光谱识别方法 | |
US9240262B1 (en) | Systems and methods for distributed pressure sensing | |
CN105136359A (zh) | 基于光纤传感器测得的梁表面应变值计算工作载荷的方法 | |
CN105403161B (zh) | 一种利用光纤传感器检测混凝土结构裂缝宽度的方法 | |
CN104132756B (zh) | 一种利用正交偏振模双峰反射谱的光子晶体光纤光栅的压力传感方法 | |
CN105526879A (zh) | 基于光纤光栅的卫星大阵面天线变形在轨测量系统及方法 | |
CN104154874B (zh) | 基于光纤传感的钢筋混凝土锈胀开裂的监测装置及方法 | |
CN102645245B (zh) | 基于光纤布里渊散射的分布式流体压力和温度同时测量方法 | |
Floris et al. | Experimental study of the influence of FBG length on optical shape sensor performance | |
CN106872071A (zh) | 一种基于光频域反射技术的测温方法 | |
CN103913592A (zh) | 一种基于光纤Bragg光栅的应力式风速风向传感器 | |
Li et al. | FBG-based online monitoring for uncertain loading-induced deformation of heavy-duty gantry machine tool base | |
CN103148894A (zh) | 基于光纤布拉格光栅的角度应力传感器 | |
CN101797928A (zh) | 基于半自由度封装fbg的轨道交通计轴装置 | |
CN104279986A (zh) | 一种活塞式水银光纤Bragg光栅倾角传感器及其使用方法 | |
CN106289600A (zh) | 一种光纤应力传感器件 | |
Jiao et al. | Investigation into coaxial cable Fabry–Perot interferometers for strain measurement and crack detection in RC structures | |
CN105890533B (zh) | 一种材料表面应变光纤光栅反向差动检测传感器件 | |
CN102175364A (zh) | 一种相似模拟实验用垂直应力光纤光栅压力传感器 | |
CN104237561A (zh) | 基于光纤光栅曲率感知网络的空间加速度测量方法和装置 | |
CN202748004U (zh) | 光纤光栅位移传感器 | |
CN201945414U (zh) | 一种相似模拟实验用垂直应力光纤光栅压力传感器 | |
CN203337110U (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: 20160608 Termination date: 20170929 |
|
CF01 | Termination of patent right due to non-payment of annual fee |