CN113532307B - 基于Michelson光纤干涉仪的大范围应变传感器 - Google Patents
基于Michelson光纤干涉仪的大范围应变传感器 Download PDFInfo
- Publication number
- CN113532307B CN113532307B CN202111053113.4A CN202111053113A CN113532307B CN 113532307 B CN113532307 B CN 113532307B CN 202111053113 A CN202111053113 A CN 202111053113A CN 113532307 B CN113532307 B CN 113532307B
- Authority
- CN
- China
- Prior art keywords
- range
- arm
- sensing
- strain
- optical fiber
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/161—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means
Abstract
本发明公开了一种基于Michelson光纤干涉仪的大范围应变传感器,包括宽带光源、2×2光纤耦合器、传感臂、参考臂和光谱分析仪;宽带光源发出的光束经2×2光纤耦合器分为两束相干光分别传输到传感臂和参考臂中,并在两个单模光纤的端面处发生反射,两束反射光发生干涉;光谱分析仪对接收到的干涉光进行分析,采用波长解调或者相位解调的方式对干涉光数据进行处理,根据谱线移动量或空间频率移动量推算得到实际应变量。本发明提出的大范围应变传感器具有高灵敏度,大动态范围,制作过程简单与成本低廉的优点,利于大规模市场应用。
Description
技术领域
本发明涉及应变传感器技术领域,具体而言涉及一种基于Michelson光纤干涉仪的大范围应变传感器。
背景技术
光纤传感器作为一种将被测对象转换为可测的光信号的仪器,具有极高的灵敏度
和精度,并且抗电磁干扰,耐高温高压,耐腐蚀,高绝缘强度,可用于高压、电器噪声、高温、
腐蚀和其他恶劣环境,且几何形状具有多方面的适应性,可制成任意形状的传感各种不同
物理信息(声、磁、温度、旋转、应变等)的光纤传感器。现有常见的应变传感器大多基于法布
里-珀罗干涉仪(FPI),根据已有文献,有人通过二氧化碳激光器将光纤端面形成弧形,然后
通过电弧放电的方法制作了FPI,此方法制成的应变传感器虽然灵敏度能够达到几十皮米
每微应变,但是二氧化碳激光加工手段复杂,设备昂贵,成本太高不利于推广使用。有人通
过放电的手段将两端高敏单模光纤的端面制作成内凹陷型,然后通过熔接机将两根光纤的
内陷端面熔接在一起,中间为空气腔,整个结构可以看成是一个法布里-珀罗腔(F-P),灵敏
度仅可以达到,虽然成本大大降低,但是传感器的灵敏度也随之下降。目前为止光
纤干涉传感器的灵敏度尚未突破;并且高灵敏度传感器的测量范围往往很小,通
常为数十至数百微应变,导致实际应用中存在困难。而大范围的应变传感器灵敏度为了实
现较大的应变长度,其灵敏度通常不高。
例如,专利号为CN106568466A的发明中提出了一种细芯微结构光纤干涉仪传感器及其温度、应变检测方法,包括宽带光源、传感头和光谱仪,传感头为带空气包层的细芯微结构光纤,细芯微结构光纤两端分别通过单模光纤连接宽带光源和光谱仪。宽带光源出射的光经过单模光纤后进入到细芯微结构光纤,由于单模光纤与细芯微结构光纤模场失配,光谱仪输出光谱包括干涉条纹光谱,温度变化或者应变会引起干涉条纹波谷波长漂移,因此在获知干涉条纹波谷波长漂移的情况下,可对应计算出温度变化或应变。该发明中的传感器需要细芯微结构光纤支持,且测量范围较小。
综上所示,目前所报道的光纤应变传感器普遍存在造价昂贵、灵敏度低、测量范围小等缺点。
发明内容
本发明针对现有技术中的不足,提供一种基于Michelson光纤干涉仪的大范围应变传感器,具有高灵敏度,大动态范围,制作过程简单与成本低廉的优点,利于大规模市场应用。
为实现上述目的,本发明采用以下技术方案:
第一方面,本发明实施例提出了一种基于Michelson光纤干涉仪的大范围应变传感器,所述应变传感器包括宽带光源、2×2光纤耦合器、传感臂、参考臂和光谱分析仪;
所述宽带光源、光谱分析仪、传感臂和参考臂分别连接在2×2光纤耦合器的四个接口上;所述传感臂和参考臂均采用单模光纤;
所述宽带光源发出的光束经2×2光纤耦合器分为两束相干光分别传输到传感臂和参考臂中,并在两个单模光纤的端面处发生反射,重新回到光纤耦合器;
所述传感臂的长度大于参考臂的长度,且两者之间的长度差满足以下条件:使经两者反射后重新回到光纤耦合器的两束反射光发生干涉,且干涉光经光纤耦合器传输至光谱分析仪;
所述光谱分析仪对接收到的干涉光进行分析,对传感臂探测到的应变量所属范围进行评估,如果应变量属于第一范围,采用波长解调的方式对干涉光数据进行处理,根据谱线移动量推算得到实际应变量;如果应变量属于第二范围,采用相位解调的方式对干涉光数据进行处理,根据空间频率移动量推算得到实际应变量;所述第一范围和第二范围不重叠且第一范围的取值小于第二范围的取值。
可选地,所述传感臂和参考臂的长度差为50um。
所述对传感臂探测到的应变量所属范围进行评估的过程包括:
根据下述公式计算得到自由光谱范围 FSR:
根据下述公式计算波长的移动范围:
式中,N代表施加于传感光纤上的应变量;
如果波长的移动范围在自由光谱范围内,则应变量属于第一范围,否则,应变量属于第二范围。
可选地,当应变量属于第二范围时,根据下述公式计算得到实际应变量:
式中,f代表反射谱线的空间频率。
第二方面,本发明实施例提出了一种基于Michelson光纤干涉仪的大范围应变传感器的工作方法,所述工作方法包括以下步骤:
S1,取两段长度完全相同的单模光纤放在微位移平台上,在显微镜下观察确定一个起点,从起点向同一个方向截取合适的长度,两段单模光纤截断点相距50um,使两段光纤产生初始的光程差;将较长的那一段单模光纤作为传感臂,短的作为参考臂;
S2,将宽带光源和光谱分析仪连接到2×2光纤耦合器一端的两个接口上,传感臂和参考臂分别连接到光纤耦合器另一端的两个接口上;
S3,通过观察光谱分析仪中所显示谱线的自由光谱范围以确定传感光纤和参考光纤的长度差,判断切割的是否合适,如果不合适,返回步骤S1重新切割,否则,进入步骤S4;
S4, 根据下述公式计算得到自由光谱范围FSR:
S5,根据下述公式计算波长的移动范围:
式中,N代表施加于传感光纤上的应变量;
如果波长的移动范围在自由光谱范围内,结束流程,将计算得到的N作为实际应变量;否则,进入步骤S6;
S6,根据下述公式计算得到实际应变量N:
式中,f代表反射谱线的空间频率。
本发明的有益效果是:
(1)本发明提出的大范围应变传感器利用两根独立的单模光纤分别作为Michelson干涉仪的传感臂与参考臂,通过光纤切割刀与微位移平台给予两根单模光纤不同的初始长度,由此产生恒定的光程差,从而导致干涉现象的发生。在进行应变测量时,传感臂长度的改变导致光程差的改变,造成光的干涉谱线发生移动,从而做到对应变进行高精度测量。
(2)本申请提出的大范围应变传感器在较小应变范围内采用高精度波长解调的方式,拥有极高的灵敏度,根据理论测算,约为同类光纤干涉传感器的数千倍;在大范围测量时内采用相位解调的方法,不仅可以实现零至数万微应变的测量范围,还可以规避大范围应变带来的波长漂移与噪声干扰。
(3)本发明提出的大范围应变传感器只需要利用单模光纤与简单的熔接、拼接技术就可以实现,成本低廉。
附图说明
图1是本发明实施例的基于Michelson光纤干涉仪的大范围应变传感器的结构示意图。
图2(a)是本发明实施例的传感臂的结构示意图;图2(b)是本发明实施例的标注了截断点的参考臂的原始单模光纤的结构示意图。
图3为本发明实施例的光在光纤末端的反射过程示意图。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
需要注意的是,发明中所引用的如“上”、“下”、“左”、“右”、“前”、“后”等的用语,亦仅为便于叙述的明了,而非用以限定本发明可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本发明可实施的范畴。
实施例一
图1是本发明实施例的基于Michelson光纤干涉仪的大范围应变传感器的结构示意图。参见图1,该应变传感器包括宽带光源、2×2光纤耦合器、传感臂(传感光纤)、参考臂(参考光纤)和光谱分析仪。
宽带光源、光谱分析仪、传感臂和参考臂分别连接在2×2光纤耦合器的四个接口上;所述传感臂和参考臂均采用单模光纤。
宽带光源发出的光束经2×2光纤耦合器分为两束相干光分别传输到传感臂和参考臂中,并在两个单模光纤的端面处发生反射,重新回到光纤耦合器。
传感臂的长度大于参考臂的长度,且两者之间的长度差满足以下条件:使经两者反射后重新回到光纤耦合器的两束反射光发生干涉,且干涉光经光纤耦合器传输至光谱分析仪。当应力作用在传感臂上时,这段光纤的腔长发生改变,传感臂与参考臂的光程差随之改变,引起干涉谱线移动,此时,即根据干涉谱线的移动量从而对应变进行高精度测量。另外,该应变传感器仅利用单模光纤,价格低廉,减少成本;并且应变造成的腔长变化,完全作用在传感臂的腔长上,大幅度提高了装置的灵敏度。
本实施例创新地利用了两段单模光纤实现Michelson干涉仪,一根光纤为传感臂,一根光纤为参考臂。一束光经过2×2光纤耦合器之后,分为两束相干光分别传输到传感臂与参考臂中。两束光在光纤的端面处发生反射后,重新回到光纤耦合器并发生干涉,干涉光从耦合器的另一端输出,连接到光谱分析仪。当传感光纤感受到应变时,两束光的光程差随之改变,反映为谱线的移动(小范围)或者自由光谱范围的变化(大范围)。针对两种不同的情况,分别采用波长解调与相位解调对干涉数据进行处理。具体的,光谱分析仪对接收到的干涉光进行分析,对传感臂探测到的应变量所属范围进行评估,如果应变量属于第一范围,采用波长解调的方式对干涉光数据进行处理,根据谱线移动量推算得到实际应变量;如果应变量属于第二范围,采用相位解调的方式对干涉光数据进行处理,根据空间频率移动量推算得到实际应变量;所述第一范围和第二范围不重叠且第一范围的取值小于第二范围的取值。本实施例通过将应力作用产生的微小的形变施加到传感光纤上,使得Michelson干涉仪的光程差发生变化。根据干涉图谱的移动量测出两束反射光光程差的改变进而得出外界施加的应力。波长解调与相位解调同时应用的灵活手段既顾及到了小范围精确测量的需求,同时满足了大范围动态测的场景。
当光被耦合器分束后,在参考臂与传感臂的末端分别发生反射,如图3所示。两束反射光回到耦合器的位置时由于固定的相位差,发生干涉,干涉光强可以被描述为:
这是一个典型的Michelson干涉仪,他的自由光谱范围(FSR)可以被描述为m级次的波谷与其相邻波谷的距离:
当小范围应变作用在传感光纤上时,其波谷的移动代表了传感器的灵敏度,并且可以被描述为:
式中,N代表施加于传感光纤上的应变量,它会导致光纤的拉伸,从而造成波谷的平移。这一公式可以被用于描述小范围应变时,波长解调的方式。当应变范围扩大,波长的移动超出仪器的测量范围,则应该采用相位解调的方式,干涉谱线的空间频率可以被描述为:
式中,f代表反射谱线的空间频率因此,应变导致空间频率的移动可以被描述为:
实施例二
本发明实施例提出了一种基于Michelson光纤干涉仪的大范围应变传感器的工作方法,该工作方法包括以下两个阶段:
(一)制备过程
制备材料需要单模光纤(SMF-28, Corning),光谱分析仪(AQ6370D, Yokogawa,Optical Spectrum Analyzer, OSA)、宽带光源(Benchtop Broadband Source, BBS)、光纤切割刀(CKFC-1, CommKing)微位移平台、光纤夹具,显微镜,拉力传感器,2×2光纤耦合器。
首先如图2(a)和图2(b)所示,取两段长度完全相同的单模光纤放在微位移平台上面,在显微镜下观察确定一个起点,从起点向同一个方向截取合适的长度,两段单模光纤截断点相距50um。两段光纤会产生初始的光程差。同时,将较长的那一段单模光纤作为传感臂,短的作为参考臂。其中,图2(a)是传感臂的结构示意图。图2(b)是标注了截断点的参考臂的原始单模光纤的结构示意图。三角形标记所在位置就是截断点,图2(a)中,截断点维持在原始单模光纤顶端,图2(b)中,截断点距离原始单模光纤顶端50um,从而截取后的参考臂的长度相对于传感臂短了50um。
2×2光纤耦合器两端各有两个接口,连接时,将宽带光源(BBS)和光谱分析仪(OSA)连接到一端的两个接口上。传感光纤和参考光纤分别连接到另一端的两个接口上。光会通过耦合器1:1的分配到传感光纤和参考光纤上。两路反射光返回到耦合器中时,发生干涉,并且传播到光谱分析仪(OSA)中,如图3所示。由于传感光纤比参考光纤长50um,通过观察光谱分析仪中所显示谱线的自由光谱范围(FSR)确定实际上传感光纤和参考光纤的长度差,判断切割的是否合适,不合适的话,重新切割。
该大范围应变传感器中未利用到熔接机等设备,只需要将单模光纤进行割断,加工手段简单,有利于大规模市场应用。
(二)测量过程
经过上述步骤之后,实验装置连接完成,光从宽带光源(BBS)中发出,经过光耦合器,进入两段光纤,将传感臂去除两侧的SMF涂层,将其一侧粘在固定台阶上,另一侧粘贴到距离为20cm的平移阶段,如图1所示。平移阶段移走一段距离,通过光谱分析仪得出在应变范围内的光谱变化,根据干涉图谱的谷值呈现的红移得出应变的大小。
两束反射光回到耦合器的位置时由于固定的相位差,发生干涉,干涉光强可以被描述为:
式中,n为本装置中所使用到的单模光纤的折射率,λ为被追踪波谷的波长,L为参考臂与传感臂之间的长度差。这是一个典型的Michelson干涉仪,他的自由光谱范围(FSR)可以被描述为m级次的波谷与其相邻波谷的距离:
当小范围应变作用在传感光纤上时,其波谷的移动代表了传感器的灵敏度,并且可以被描述为:
式中,N代表施加于传感光纤上的应变量,它会导致光纤的拉伸,从而造成波谷的平移。这一公式可以被用于描述小范围应变时,波长解调的原理。当应变范围扩大,波长的移动超出仪器的测量范围,则应该采用相位解调的方式,干涉谱线的空间频率可以被描述为:
式中,f代表反射谱线的空间频率因此,应变导致空间频率的移动可以被描述为:
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (4)
1.一种基于Michelson光纤干涉仪的大范围应变传感器,其特征在于,所述应变传感器包括宽带光源、2×2光纤耦合器、传感臂、参考臂和光谱分析仪;
所述宽带光源、光谱分析仪、传感臂和参考臂分别连接在2×2光纤耦合器的四个接口上;所述传感臂和参考臂均采用单模光纤;
所述宽带光源发出的光束经2×2光纤耦合器分为两束相干光分别传输到传感臂和参考臂中,并在两个单模光纤的端面处发生反射,重新回到光纤耦合器;
所述传感臂的长度大于参考臂的长度,且两者之间的长度差满足以下条件:使经两者反射后重新回到光纤耦合器的两束反射光发生干涉,且干涉光经光纤耦合器传输至光谱分析仪;
所述光谱分析仪对接收到的干涉光进行分析,对传感臂探测到的应变量所属范围进行评估,如果应变量属于第一范围,采用波长解调的方式对干涉光数据进行处理,根据谱线移动量推算得到实际应变量;如果应变量属于第二范围,采用相位解调的方式对干涉光数据进行处理,根据空间频率移动量推算得到实际应变量;所述第一范围和第二范围不重叠且第一范围的取值小于第二范围的取值;
所述对传感臂探测到的应变量所属范围进行评估的过程包括:
根据下述公式计算得到自由光谱范围FSR:
式中,λ为被追踪波谷的波长,L为参考臂与传感臂之间的长度差,n为传感臂和参考臂所使用到的单模光纤的折射率;λdip(m)为m级次的波谷;λdip(m-1)为m-1次级的波谷;
根据下述公式计算波长的移动范围:
式中,N代表施加于传感光纤上的应变量;
如果波长的移动范围在自由光谱范围内,则应变量属于第一范围,否则,应变量属于第二范围。
2.根据权利要求1所述的基于Michelson光纤干涉仪的大范围应变传感器,其特征在于,所述传感臂和参考臂的长度差为50um。
4.一种基于Michelson光纤干涉仪的大范围应变传感器的工作方法,其特征在于,所述工作方法包括以下步骤:
S1,取两段长度完全相同的单模光纤放在微位移平台上,在显微镜下观察确定一个起点,从起点向同一个方向截取合适的长度,两段单模光纤截断点相距50um,使两段光纤产生初始的光程差;将较长的那一段单模光纤作为传感臂,短的作为参考臂;
S2,将宽带光源和光谱分析仪连接到2×2光纤耦合器一端的两个接口上,传感臂和参考臂分别连接到光纤耦合器另一端的两个接口上;
S3,通过观察光谱分析仪中所显示谱线的自由光谱范围以确定传感光纤和参考光纤的长度差,判断切割的是否合适,如果不合适,返回步骤S1重新切割,否则,进入步骤S4;
S4,根据下述公式计算得到自由光谱范围FSR:
式中,λ为被追踪波谷的波长,L为参考臂与传感臂之间的长度差,n为传感臂和参考臂所使用到的单模光纤的折射率;λdip(m)为m级次的波谷;λdip(m-1)为m-1次级的波谷;
S5,根据下述公式计算波长的移动范围:
式中,N代表施加于传感光纤上的应变量;
如果波长的移动范围在自由光谱范围内,结束流程,将计算得到的N作为实际应变量;否则,进入步骤S6;
S6,根据下述公式计算得到实际应变量N:
式中,f代表反射谱线的空间频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111053113.4A CN113532307B (zh) | 2021-09-09 | 2021-09-09 | 基于Michelson光纤干涉仪的大范围应变传感器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111053113.4A CN113532307B (zh) | 2021-09-09 | 2021-09-09 | 基于Michelson光纤干涉仪的大范围应变传感器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113532307A CN113532307A (zh) | 2021-10-22 |
CN113532307B true CN113532307B (zh) | 2021-11-30 |
Family
ID=78123040
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111053113.4A Active CN113532307B (zh) | 2021-09-09 | 2021-09-09 | 基于Michelson光纤干涉仪的大范围应变传感器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113532307B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114414010B (zh) * | 2022-01-20 | 2023-03-24 | 福州大学 | 双弯曲梁称重传感器 |
CN114858052B (zh) * | 2022-04-26 | 2023-06-06 | 厦门大学 | 基于虚拟参考腔和游标效应的高灵敏大范围干涉测量方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020196993A1 (en) * | 2001-06-26 | 2002-12-26 | Schroeder Robert J. | Fiber optic supported sensor-telemetry system |
CN105466349B (zh) * | 2016-01-18 | 2018-07-10 | 天津大学 | 一种光频域反射中用薄包层光纤提高应变测量灵敏度方法 |
CN106644036A (zh) * | 2016-12-26 | 2017-05-10 | 华中科技大学 | 一种基于聚合物薄膜的声波探测器及双波长解调方法 |
CN108168728B (zh) * | 2017-12-11 | 2019-12-10 | 哈尔滨工程大学 | 非平衡保偏光纤双干涉仪温度应变同时测量装置及方法 |
CN112525373B (zh) * | 2020-11-10 | 2023-05-12 | 广东工业大学 | 一种基于双波长保偏光纤干涉仪的应变温度同时测量装置 |
CN112525372B (zh) * | 2020-11-10 | 2023-06-27 | 广东工业大学 | 基于保偏光纤双臂异轴干涉仪的应变温度同时测量装置及方法 |
-
2021
- 2021-09-09 CN CN202111053113.4A patent/CN113532307B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113532307A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100332833B1 (ko) | 투과형 외인성 패브리-페롯 광섬유 센서와, 이를 이용한 변형률 및 온도 측정방법 | |
CN113532307B (zh) | 基于Michelson光纤干涉仪的大范围应变传感器 | |
US5392117A (en) | Fabry-Perot optical sensing device for measuring a physical parameter | |
CN105865752B (zh) | 采用分布式偏振串扰分析仪全面评判保偏光纤特性的方法和装置 | |
US11346770B2 (en) | Optical fiber sensor for salinity and temperature measurement | |
Blanchard et al. | Two-dimensional bend sensing with a single, multi-core optical fibre | |
US6513390B1 (en) | Temperature insensitive fiber-optic torque and strain sensor | |
CA2695587A1 (en) | Physical quantity measuring apparatus utilizing optical frequency domain reflectometry and method for temperature and strain measurement using the apparatus | |
CN100340839C (zh) | 光纤应变测量仪及其测量方法 | |
CN1746641A (zh) | 光纤温度测量仪及其测量方法 | |
CN110470240A (zh) | 一种光纤曲率测量传感器及其制作方法、测量系统 | |
MacPherson et al. | Tunnel monitoring using multicore fibre displacement sensor | |
CN111220188A (zh) | 基于级联光纤错位移对增强游标效应超灵敏折射率传感器 | |
CN110375780B (zh) | Ofdr断纤续接测量的方法 | |
CN103697921A (zh) | 一种光纤传感头和基于该传感头的测量应变、应力和压力的光纤传感系统及方法 | |
CN106289600A (zh) | 一种光纤应力传感器件 | |
US7280220B2 (en) | Physical quantity measuring method and device therefor | |
Vallan et al. | Static characterization of curvature sensors based on plastic optical fibers | |
CN112710633A (zh) | 一种制备拼接结构马赫-曾德尔折射率传感器方法 | |
Liehr et al. | A novel quasi-distributed fibre optic displacement sensor for dynamic measurement | |
JP2002257520A (ja) | 光ファイバひずみ検出装置 | |
AU2020103491A4 (en) | A twin array Michelson fiber optic white light interferometry strain gauge | |
CN113494890B (zh) | 基于fpi干涉仪的光纤光栅应变传感器精度测量装置及方法 | |
KR100324117B1 (ko) | 전반사형 외부 패브리-페로 광섬유 센서와 이를 이용한변형률 측정방법 | |
Li et al. | A highly sensitive curvature sensor based on Omega shaped long-period fiber grating |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |