CN105157624B - A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail - Google Patents

A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail Download PDF

Info

Publication number
CN105157624B
CN105157624B CN201510260240.XA CN201510260240A CN105157624B CN 105157624 B CN105157624 B CN 105157624B CN 201510260240 A CN201510260240 A CN 201510260240A CN 105157624 B CN105157624 B CN 105157624B
Authority
CN
China
Prior art keywords
string
points
rail
strings
irregularity
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
CN201510260240.XA
Other languages
Chinese (zh)
Other versions
CN105157624A (en
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.)
Hunan University
Original Assignee
Hunan 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 Hunan University filed Critical Hunan University
Priority to CN201510260240.XA priority Critical patent/CN105157624B/en
Publication of CN105157624A publication Critical patent/CN105157624A/en
Application granted granted Critical
Publication of CN105157624B publication Critical patent/CN105157624B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

The invention discloses a kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail, on the premise of system cost is not increased, by comparing from the specific chord lengths of 1m, with reference to 3:97 chord lengths divide ratio, optimal FIR inverse filter designs method, combine that (3 1 dimension sensors are simultaneously using being 3 points of inclined strings with the flexible of 3 sensors of group, only it is 2 strings with 1# and 3# sensors), short wavelength and long-wavelength detection are combined together, the actual irregularity waveform of track that can be in accurate recovery full recovery band.

Description

A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail
Technical field
The present invention relates to a kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail.
Background technology
Theoretical research and engineering practice are it has been proved that track irregularity is to cause rolling stock to produce vibration, cause wheel track Active force increase, the principal element for causing circuit aspect directly limitation road speed.Therefore, track irregularity is grasped in accurate measurement Actual conditions be to realize the precondition that Scientific Assessment and monitoring management are carried out to rail smooth state.
Longitudinal irregularity waveform present on track be spatial distribution randomness deformation, wave-length coverage from tens millimeters to Over one hundred rice, wavelength is longer, and amplitude is bigger, is highly difficult by being detected to irregularity entirely accurate present on track.By Different in the performance of dynamics of vehicle, vehicle only has response to the track irregularity waveform of a wavelength range.Therefore, detection system System only needs that the irregularity waveform in particular range of wavelengths is restored and measured.Under high-speed condition, 20m~70m wavelength Irregularity, will make the relatively low car body of intrinsic frequency occur exciting, general railway or heavy haul railway (speed 80Km/h~120Km/h) Only medium wavelength (5m~12m) irregularity is any limitation as, high-speed railway is especially paid attention to uneven to short wavelength (30mm~1000mm) Suitable control.Set accordingly and recover band for 30mm~60m.
At present, countries in the world, which are used for measuring the method for longitudinal irregularity, can be summarized as inertial reference method and chord measurement two is big Class.
Inertial reference method is by acceleration transducer to set up an inertial reference benchmark in the car body of motion, utilizes car Displacement transducer on body carrys out relative position of the measure track relative to benchmark, so as to obtain rail surface in inertial coodinate system Relative position.Principle is as shown in Figure 1.
If A (t) represents the vertical vibration acceleration of car body, it is measured by the accelerometer on car body floor , the vibration and impact of car body are relatively small, and maximum is no more than several g (acceleration of gravity).D (L) is represented between axle box and car body Vertical relative displacement, it is measured by the displacement transducer between axle box and accelerometer, to acceleration make Corresponding Integral Processing just can obtain the movement locus of accelerometer mount point, be further that can obtain as algebraic operation with D (L) The track P (L) of rail top face.Its operation expression is:
P (L)=∫ ∫ A (t) dt2-D(L) (1)
Chord measurement drives several displacement transducers along track longitudinal direction according to the spacing word order being pre-designed, and is passed by displacement Sensor measures light source incidence point and rail level distance, is then divided according to chord length than being derived by string measured value, most basic string surveys system System is 2 strings or 3 equal string measuring systems.By taking the string such as 3 points as an example, principle is as shown in Figure 2.
String measured value bo and the actual value bb' for being not equal to track irregularity, but actual value and the product of ssystem transfer function.
Inertial reference method disclosure satisfy that measurement request in theory.But it is due to axle box acceleration caused by track irregularity Dynamic range is very big, and to measure 0.lm to the irregularity of 50m wavelength, resolving accuracy is 1 millimeter, then needs measurement acceleration to move State scope is 0.0013g-311g, and frequency range is 0.56Hz-278Hz.Current sensor and Electrical Measuring Instrument are difficult at this Ensure such resolving accuracy in the big dynamic range of sample.Meanwhile, its accuracy of detection is influenceed by train running speed, and shortwave is examined Survey has certain blind area, is mainly used on large-scale track detection vehicle.
Chord measurement restores accuracy and restricted by transmission function, is broadly divided into 2 strings, the strings such as 3 points, 3 points of inclined strings and at 4 points 4 kinds of measuring methods such as inclined string, (chord length takes 330mm, 3 points of inclined string strings to the amplitude-versus-frequency curve of respective transfer functions as shown in Figure 3 Length is divided into 1:10,4 points of inclined string chord lengths are divided into 3:17:1).
Wherein, string (versed sine method) amplitude gain such as 2 strings and 3 points is very fast with irregularity wavelength change, is mainly used in pair Long wavelength's composition is restored, and now its chord length is generally tens of rice.2 strings can only be to the irregularity waveforms of wavelength > chord lengths Restored, the strings such as 3 points can only restore to the irregularity waveform of the chord lengths of wavelength > 1/2.3 points of inclined strings and 4 points of inclined strings (swearing method partially) transmission function is complex, and amplitude gain is gentle with wavelength change, is mainly used in restoring short wavelength's composition, Its chord length is general less than 1 meter.
Recover band for 30mm~60m of setting, 4 kinds of methods of the above can not effectively meet demand when being used alone.
Inertial reference method is detected and data handling system is complex and expensive, and chord measurement can only be to specific long wave Long or short wavelength's composition is restored, it is impossible to which 30mm~60m full recovery bands of setting are restored.
The content of the invention
The technical problems to be solved by the invention are, in view of the shortcomings of the prior art, providing a kind of for measuring rail longitudinal direction The compound chord measurement of longitudinal irregularity.
In order to solve the above technical problems, the technical solution adopted in the present invention is:One kind is used to measure rail longitudinal direction height The compound chord measurement of irregularity, comprises the following steps:
1) three laser sensors are set respectively on the left rail and right rail of track;
2) i=0 is made;
3) the light source incidence point that three laser sensor synchronous acquisitions on left rail or right rail are arrived to Rail Surface away from From data, the range data of three laser sensor collections is stored in array 1# (i), 2# (i), 3# (i) respectively;
4) i=i+1 is made, is repeated the above steps 3), until i=50000;
5) by step 4) obtained three number groups are assigned to array 11# (j), 12# (j), 13# (j) respectively;
6) exceptional value in array 11# (j), 12# (j), 13# (j) is removed, new array 111# (j), 112# (j) is obtained, 113#(j);
7) 3 points of inclined strings are calculated:Make s1=111# (j), s2=112# (j), s3=113# (j), is utilized3 points of inclined string string measured value y (x) ' are calculated, it is inverse to design 3 points of inclined strings using 3 points of inclined string string measured values Wave filter, recovers string and surveys waveform 30mm~1200mm shortwave long components;Wherein s=1000mm, a=30mm, b=970mm;
8) 2 points of inclined strings are calculated:Make k=0.Whether be zero, if so, 21# (k)=111# (j), 23# (k) if judging j%100 =113# (j), k=k+1;Otherwise, j=j+1 is made, repeats the above steps and remainder judgement is carried out to j values, until j=50000;
9) s is made1'=21# (k), s2'=23# (k), utilizes formula y (x)=s1'-s3' 2 string string measured value y (x) are calculated, And 2 points of inclined string inverse filters are designed using 2 string string measured value y (x), recover string and survey waveform 1200mm~60m long wave long components.
Compared with prior art, the advantageous effect of present invention is that:The present invention is used merely in the traditional chord measurement of improvement On the basis of short wavelength (30mm~1m) or long wavelength (1m~60m) string survey waveform recovery, by comparing from the specific strings of 1m It is long, with reference to 3:97 chord lengths divide ratio, optimal FIR inverse filter designs method, the flexibly means such as combination with 3 sensors of group, Short wavelength and long-wavelength detection are combined together, the actual irregularity waveform of track that can be in accurate recovery full recovery band.
Brief description of the drawings
Fig. 1 is inertial reference method schematic diagram;
Fig. 2 is chord measurement schematic diagram.
Fig. 3 is the transmission function amplitude versus frequency characte figure that four kinds of strings survey method;
Fig. 4 is detecting system hardware block diagram of the present invention;
Fig. 5 (a) is s=330mm, the inclined string amplitude frequency curve figures of 3 points of a=30mm, 40mm, 50mm;Fig. 5 (b) is a=30mm, b:A=1, b:A=5, b:A=10, b:A=15, b:A=20, b:The inclined string amplitude frequency curve figures of 3 points of a=30;Fig. 5 (c) is a= The inclined string amplitude frequency curve figures of 3 points of 30mm, s=330mm, 500mm, 800mm, 1000mm;
Fig. 6 is 2 string Frequency Sampling Method design flow diagrams of the invention;
Fig. 7 is 3 string advantest method design flow diagrams of the invention;
Fig. 8 is that string of the present invention surveys waveform recovery flow chart;
Fig. 9 is three kinds of waveform spatial domain figures of the invention;
Figure 10 is original waveform and the spectral contrast figure for restoring waveform.
Embodiment
Occur to avoid projecting situation beyond rail top by light during bend caused by car body mounting means, rail is indulged It is placed in longitudinal irregularity detecting system on a detection dolly, the dynamic that the progress irregularity that advances can be drawn with Operation Van is surveyed Amount, fully improves operating efficiency.Detecting system is tieed up laser displacement sensor, concentrator, data and adopted by machinery mount, high accuracy 1 A few part compositions such as collecting system, hardware system block diagram is as shown in Figure 4.
Detecting system uses 2 groups of totally 6 dimension laser displacement sensors of high accuracy 1, every group 3, completes side rail longitudinal direction The measurement of longitudinal irregularity.It is placed in 3 sensors of group in a mounting box, light source incidence fixed fire is on one wire, and incident Light is each perpendicular to rail top center line.To ensure 3 sensor synchronous acquisition rail top irregularity data with Spatial zonation from sampling Realization, a rotary photoelectric encoder is installed on wheel shaft, equidistant square-wave pulse signal is exported, it is same that signal is accessed into sensor Step triggering input, sensor can pulse effective edge along synchronism output measurement data, it is ensured that between multiple sensors are synchronously waited Away from sampling.
Implementation process of the present invention is as follows:
The first step:Select suitable chord length and divide ratio
The selected minimum 30mm of laser displacement sensor light source center spacing.Ratio and overall string are divided to compare chord length The long influence to transmission function amplitude versus frequency characte, has carried out following 3 experiments and (has set overall chord length as s, chord length is divided than being a:b):
1. s=330mm, a=30mm, 40mm, 50mm are taken, shown in 3 points of inclined string amplitude frequency curve such as Fig. 5 (a);
2. a=30mm, b are taken:A=1, b:A=5, b:A=10, b:A=15, b:A=20, b:A=30,3 points of inclined string width Shown in frequency curve such as Fig. 5 (b);
3. a=30mm, s=330mm, 500mm, 800mm, 1000mm, 2 strings and 3 points of inclined string amplitude frequency curve such as Fig. 5 are taken (c) shown in.
If transmission function is too small, inverse filter, which excessively amplifies string measured value, can cause violent vibration, therefore need to set width in figure It is worth the blue horizontal warning line of gain 0.1.Drawn a conclusion by experimental result as follows:
1. during wavelength=a, 3 points of inclined string transmission function amplitudes are 0;
2. 3 points of inclined string amplitude gains increase with wavelength and are gradually reduced, b:A ratios are bigger, and gain reduction is slower;
3. a=30 is fixed, during chord length s=1000mm, in cordon bleu line, it is 30mm- that 3 points of inclined strings, which recover wavelength, 2000mm, it is 1.2m-63m that 2 strings, which recover wavelength, and both recover wavelength and effectively overlapped, and meeting actual recovery band 30mm-60m will Ask.
Therefore, chord length s=1000mm, a=30mm, b=970mm are selected.
Second step:Design corresponding inverse filter
2 string transmission functions:
H (ω)=1-e-j*ω*s, s is chord length (3)
3 points of inclined string transmission functions
S is chord length, and chord length is divided than being a:b (4)
The chord length s and its chord length that the first step is selected, which are divided, compares a:B substitutes into formula (3) and formula (4), to design corresponding inverse filter Ripple device.Wherein, 2 string transmission functions are relatively simple, are designed using Frequency Sampling Method, uneven for recovering 1.2m~60m Waveform long wavelength's composition is surveyed along string;3 points of inclined string transmission functions are complex, and phase-frequency characteristic is difficult to calculate, and is carried out using advantest method Design, for recovering 30mm~1200mm short wavelength's compositions.
2 string design flow diagrams are as shown in Figure 6.
3 points of inclined string design flow diagrams are as shown in Figure 7.
3rd step:String surveys the acquisition of data and waveform restores
Detecting system receives the trigger pulse that first encoding device is sent, displacement of the output with 31 dimension sensors of group per 2mm Value.3 points of inclined strings are 2mm per the sampling interval between twice, and 2 strings are 200mm.
2 string string measured values
Y (x)=s1-s3, s1、s3The shift value (5) measured for respective sensor
3 points of inclined string string measured values
s1、s2、s3The shift value (6) measured for respective sensor
The string measured value of 2 strings and 3 points of inclined strings is substituted into the corresponding inverse filter of second step design, you can ripple is surveyed to string Shape is restored accordingly.Provide, rail smooth degree is once evaluated per 100m, overall string according to EN13231-3-2000 It is as shown in Figure 8 that flow is restored in the acquisition and waveform for surveying data.
For ease of observation checking, design one section of long 100m irregularity original waveform, by wavelength be respectively 50mm, 200mm, 1m, 5m, 10m sine-wave superimposed are formed.Remove point of the 2# sensors without string measured value, obtain string and survey data.Respectively through 2 strings With 3 points of inclined string lifterings, obtain short wavelength's composition and long wavelength's composition restores waveform.Three and corresponding frequency spectrum are as shown in Figure 9.
The string that Fig. 9 is actual to be adopted is surveyed waveform and differed greatly with original waveform, and the long wavelength's composition obtained through liftering is just For the low-frequency range of original waveform, its a length of high band of shortwave.Figure 10 restores waveform and essentially coincided with the original waveform of corresponding section, former The spectral components that beginning waveform contains accurately are embodied in waveform is restored, it was demonstrated that the feasibility of this method.

Claims (1)

1. a kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail, it is characterised in that comprise the following steps:
1) three laser sensors are set respectively on the left rail and right rail of track;
2) i=0 is made;
3) distance number of the light source incidence point that three laser sensor synchronous acquisitions on left rail or right rail are arrived to Rail Surface According to the range data of three laser sensor collections is stored in array 1# (i), 2# (i), 3# (i) respectively;
4) i=i+1 is made, is repeated the above steps 3), until i=50000;
5) by step 4) obtained three number groups are assigned to array 11# (j), 12# (j), 13# (j) respectively;
6) exceptional value in array 11# (j), 12# (j), 13# (j) is removed, new array 111# (j), 112# (j), 113# is obtained (j);
7) 3 points of inclined strings are calculated:Make s1=111# (j), s2=112# (j), s3=113# (j), is utilized3 points of inclined string string measured value y (x) ' are calculated, it is inverse to design 3 points of inclined strings using 3 points of inclined string string measured values Wave filter, recovers string and surveys waveform 30mm~1200mm shortwave long components;Wherein s=1000mm, a=30mm, b=970mm;
8) 2 points of inclined strings are calculated:K=0 is made, whether judge j%100 is zero, if so, 21# (k)=111# (j), 23# (k)= 113# (j), k=k+1;Otherwise, j=j+1 is made, repeats the above steps and remainder judgement is carried out to j values, until j=50000;
9) s is made1'=21# (k), s3'=23# (k), utilizes formula y (x)=s1'-s3' calculate 2 string string measured value y (x), and profit 2 points of inclined string inverse filters are designed with 2 string string measured value y (x), recovers string and surveys waveform 1200mm~60m long wave long components.
CN201510260240.XA 2015-05-20 2015-05-20 A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail Expired - Fee Related CN105157624B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510260240.XA CN105157624B (en) 2015-05-20 2015-05-20 A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510260240.XA CN105157624B (en) 2015-05-20 2015-05-20 A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail

Publications (2)

Publication Number Publication Date
CN105157624A CN105157624A (en) 2015-12-16
CN105157624B true CN105157624B (en) 2017-11-03

Family

ID=54798568

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510260240.XA Expired - Fee Related CN105157624B (en) 2015-05-20 2015-05-20 A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail

Country Status (1)

Country Link
CN (1) CN105157624B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109798850B (en) * 2019-02-01 2020-10-20 湖南大学 Rail corrugation measuring method
CN110488042B (en) * 2019-07-22 2021-12-07 中车青岛四方机车车辆股份有限公司 Train acceleration detection method and system, electronic equipment and storage medium
CN111307079A (en) * 2020-03-23 2020-06-19 中国铁建重工集团股份有限公司 Detection device
CN112590559B (en) * 2021-01-04 2022-03-08 西南交通大学 Configuration method for magnetic suspension track smoothness deviation chord survey
CN112902871A (en) * 2021-02-04 2021-06-04 暨南大学 Tramcar track irregularity detection method
CN113071529B (en) * 2021-04-09 2022-09-16 中国铁道科学研究院集团有限公司 Track elasticity detection method and device
CN113111554B (en) * 2021-04-09 2024-03-08 中国铁道科学研究院集团有限公司铁道建筑研究所 Deformation index calculation method for temporary overhead line of high-speed railway passing 80km per hour
CN113319152B (en) * 2021-04-21 2022-10-28 邯郸钢铁集团有限责任公司 Method for realizing intelligent pressure straightening of steel rail
CN113548068B (en) * 2021-07-23 2023-09-12 中车长春轨道客车股份有限公司 Rail surface irregularity detection device and detection method
CN113959354B (en) * 2021-11-26 2022-10-11 中国铁道科学研究院集团有限公司 Device and method for measuring structural deformation based on optical fiber sensing
CN114104034A (en) * 2021-12-24 2022-03-01 湖南凌翔磁浮科技有限责任公司 Permanent magnet track geometric irregularity detection system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101758835A (en) * 2008-12-23 2010-06-30 中国铁道科学研究院基础设施检测研究所 Framework type track gauge and track direction measuring device
CN101929851A (en) * 2009-06-18 2010-12-29 上海友馨数字技术有限公司 Laser high precision and high speed detection apparatus for vertical smoothness of rail
JP5045908B2 (en) * 2007-06-12 2012-10-10 綜合計測株式会社 Railroad track displacement measuring method and railroad track displacement measuring device
CN202644332U (en) * 2012-05-09 2013-01-02 中国铁道科学研究院基础设施检测研究所 Rail top surface irregularity detecting device
CN103174072A (en) * 2011-12-26 2013-06-26 同济大学 Rail surface detection method based on asymmetrical chord offset method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5045908B2 (en) * 2007-06-12 2012-10-10 綜合計測株式会社 Railroad track displacement measuring method and railroad track displacement measuring device
CN101758835A (en) * 2008-12-23 2010-06-30 中国铁道科学研究院基础设施检测研究所 Framework type track gauge and track direction measuring device
CN101929851A (en) * 2009-06-18 2010-12-29 上海友馨数字技术有限公司 Laser high precision and high speed detection apparatus for vertical smoothness of rail
CN103174072A (en) * 2011-12-26 2013-06-26 同济大学 Rail surface detection method based on asymmetrical chord offset method
CN202644332U (en) * 2012-05-09 2013-01-02 中国铁道科学研究院基础设施检测研究所 Rail top surface irregularity detecting device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
三点偏弦法复原轨面不平顺波形的理论及研究;程樱等;《华东交通大学学报》;20110228;第28卷(第1期);第42-46页 *
数字逆滤波技术在轨道不平顺检测中的应用;张国锋等;《信号处理》;20041231;第20卷(第6期);第667-670页第1-3节 *
高速铁路无砟轨道不平顺谱;康熊等;《中国科学:技术科学》;20141231;第44卷(第7期);第688页第2节 *

Also Published As

Publication number Publication date
CN105157624A (en) 2015-12-16

Similar Documents

Publication Publication Date Title
CN105157624B (en) A kind of compound chord measurement for being used to measure the longitudinal longitudinal irregularity of rail
CN102901550B (en) Method for implementing vehicle-mounted dynamic weighing
CN203349914U (en) Sensor module for measuring axle speed and axle weight of double wheel rut vehicle
CN107244335B (en) Wave abrasion detection device and method
CN103335601B (en) Quick and automatic detection device for overall dimension of moving vehicle
CN108515984B (en) Wheel damage detection method and device
CN104032629B (en) A kind of vertical track long wave irregularity on-line monitoring method and system
CN102498000B (en) Method and apparatus for estimating the shape of a tyre contact patch
CN102874278B (en) Vehicle-mounted detection method fusing vehicle speed information and axle box vertical acceleration information for wheel flats
CN104215421A (en) Quick bridge impact coefficient determination method
CN202298394U (en) Detection system for dynamically detecting track gauge of railway track
CN105000033B (en) A kind of check and evaluation system of track geometry irregularities
CN202644332U (en) Rail top surface irregularity detecting device
CN102680064B (en) Vehicle dynamic weighing system
CN104120644A (en) Gravitational acceleration sensor based road surface flatness detection method
CN103322936A (en) Wheel set size on-line monitoring method
CN206781779U (en) Laser scanning rail gauge measuring apparatus under a kind of track checking car
CN105929025B (en) One kind is based on the continuous wheel tread of the time and space and rail deformation detection method
CN103473948B (en) Piezoelectric cable laying structure and the trace lateral attitude recognition methods of overloaded vehicle wheel
CN104598753A (en) Bridge moving vehicle load recognition method based on Brakhage V method
CN104195930A (en) System and method for detecting pavement evenness based on multiple sensors
CN107423692A (en) A kind of rail corrugation fault detection method based on wavelet-packet energy entropy
CN103307982A (en) On-line monitoring device for dimensions of wheel sets
CN104005324B (en) A kind of detection system of pavement structure information
CN104864949A (en) Vehicle dynamic weighing method and device thereof

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
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: 20171103