CN105320057A - Synchronous movement control method of double vibration tables on the basis of coordinate-transformation matrix - Google Patents

Synchronous movement control method of double vibration tables on the basis of coordinate-transformation matrix Download PDF

Info

Publication number
CN105320057A
CN105320057A CN201510282045.7A CN201510282045A CN105320057A CN 105320057 A CN105320057 A CN 105320057A CN 201510282045 A CN201510282045 A CN 201510282045A CN 105320057 A CN105320057 A CN 105320057A
Authority
CN
China
Prior art keywords
transition matrix
freedom
output
input
coordinate
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.)
Granted
Application number
CN201510282045.7A
Other languages
Chinese (zh)
Other versions
CN105320057B (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.)
Beijing Institute of Spacecraft Environment Engineering
Original Assignee
Beijing Institute of Spacecraft Environment Engineering
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 Beijing Institute of Spacecraft Environment Engineering filed Critical Beijing Institute of Spacecraft Environment Engineering
Priority to CN201510282045.7A priority Critical patent/CN105320057B/en
Publication of CN105320057A publication Critical patent/CN105320057A/en
Application granted granted Critical
Publication of CN105320057B publication Critical patent/CN105320057B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Automation & Control Theory (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The present invention discloses a synchronous movement control method of double vibration tables on the basis of a coordinate-transformation matrix. The synchronous movement control method comprises: converting linear measurement of two control sensors on the double vibration tables to corresponding two degrees of freedom such as translation and rotation, and inputting the coordinate-transformation matrix; converting signals outputted by a multiple-input-multiple-output system to driving signals actually inputted to the power amplifiers of the double vibration tables, and outputting the coordinate-transformation matrix; and synchronously controlling the two-degree of freedom motion, specifically, performing closed-loop control according to the translation degree of freedom reference test magnitude, the rotation degree of freedom reference test magnitude, the transfer function of the system and the measured two-degree of freedom response through the multiple-input-multiple-output system, and realizing the synchronous movement of the double vibration tables. The synchronous movement control method of double vibration tables on the basis of a coordinate-transformation matrix has high control precision, and may preferably solve the technical question of synchronous movement control of the double vibration tables.

Description

Based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix
Technical field
The invention belongs to spacecraft dynamic environmental test field, be specifically related to a kind of Double earthquakes Synchronous motion control method based on coordinate conversion matrix.
Background technology
Double earthquakes is synchronized with the movement and encourages pilot system to mainly contain two kinds of control modes at present.A kind of is power synchronisation control means based on Double earthquakes moving-coil current feedback; A kind of is standard square formation control method based on two reference mark on Double earthquakes.Double earthquakes utilizes phase controller to regulate the drive singal inputing to power amplifier according to the phase differential of two power amplifier feedback signals based on the power synchronisation control means of moving-coil current feedback usually, coordinates two shaking table synchronized movement; And utilize multi-input multi-output control system directly to carry out the control of standard square formation based on the reference test condition that the standard square formation controlling party rule at two reference mark on Double earthquakes is identical according to two reference mark, make two shaking table synchronized movement.The current power synchronisation control means based on Double earthquakes moving-coil current feedback needs single shaft vibrating controller to carry out two platform synchronized movement control of shaking in conjunction with phase controller, system complex; And all there is the not high problem of control accuracy with the standard square formation control method based on two reference mark on Double earthquakes in the method.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of Double earthquakes Synchronous motion control method based on coordinate conversion matrix to make control accuracy high, and can solve Double earthquakes Synchronous motion control technical matters preferably.
In order to solve the problems of the technologies described above, the invention provides a kind of Double earthquakes synchronisation control means based on coordinate conversion matrix, comprising the following steps:
(1) input coordinate transition matrix is built, described structure input coordinate transition matrix step comprises according to the distance determination input coordinate transition matrix on Double earthquakes between two control sensors, and linear measurement Double earthquakes controlling sensor is converted to corresponding translation and rotation two degree of freedom by acting as of input coordinate transition matrix;
(2) output coordinate transition matrix is built, described structure output coordinate transition matrix step comprise according to Double earthquakes excitation between distance determination output coordinate transition matrix, output coordinate transition matrix act as by multi-input multi-output control system export signal be converted to the actual drive singal inputing to Double earthquakes power amplifier;
(3) two degree-of-freedom motion synchro control, described two degree-of-freedom motion synchro control step is specially the white noise random signal that multi-input multi-output control system first sends the low magnitude of two-way, and this white noise random signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; On Double earthquakes, the response of two measurement points enters multi-input multi-output control system after input coordinate transition matrix, and control system calculates the transport function of whole system according to input and output signal; Multi-input multi-output control system calculates initial two-way drive singal according to the translational degree of freedom reference test magnitude arranged, rotational freedom reference test magnitude, system transter subsequently, and this two-way drive singal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; Multi-input multi-output control system is according to two degrees of freedom response level and the translation of measuring and rotates two degrees of freedom reference test magnitude and compare and the error of calculation, revises the two-way drive singal of transmission in real time, thus realizes the control of Double earthquakes synchronized movement.
Double earthquakes Synchronous motion control method based on coordinate conversion matrix provided by the invention is converted to corresponding translation according to the control linear sensor measurement on two shaking tables and rotates two degree of freedom, test magnitude and the standard side of translational degree of freedom Control experiment magnitude of shaking is identical, the test magnitude of rotational freedom is then set to the lower numerical value of a relative translational movement degree of freedom, translational degree of freedom is realized by multi-input multi-output control system, and rotational freedom is suppressed by multi-input multi-output control system, thus realize Double earthquakes synchronized movement.Control method of the present invention controls to compare with direct standard square formation with two platform current feedback synchro control of shaking, and control accuracy is higher, solves Double earthquakes Synchronous motion control technical matters preferably.
Accompanying drawing explanation
Fig. 1 is Double earthquakes synchro control of the present invention response geometric layout figure;
Wherein, 1-first acceleration transducer, 2-second acceleration transducer, the two platform geometric center of shaking of 3-.
Fig. 2 is Double earthquakes synchro control of the present invention excitation geometric layout figure;
Wherein, 4-first shaking table excitation, the excitation of 5-second shaking table.
Fig. 3 is input coordinate transition matrix design of graphics of the present invention.
Fig. 4 is output coordinate transition matrix design of graphics of the present invention.
Fig. 5 is Double earthquakes two degree-of-freedom motion synchro control process flow diagram of the present invention.
Embodiment
Below in conjunction with accompanying drawing, a kind of Double earthquakes Synchronous motion control method based on coordinate conversion matrix of the present invention is further described.
Fig. 1 is Double earthquakes synchro control of the present invention response geometric layout figure.As shown in Figure 1, the distance between the first acceleration transducer 1 shaking table placed and the second acceleration transducer 2 is 2r.First acceleration transducer 1 and the second acceleration transducer 2 are r from the distance of Double earthquakes geometric center 3.How the linear acceleration of two Vertical dimension is measured the translational degree of freedom T and the rotational freedom R that are converted to test specimen by input coordinate transition matrix by describing below.
For translation, the first acceleration transducer 1, second acceleration transducer 2 need move in same direction, mathematic(al) mean method can be adopted to calculate, that is:
T control=0.5A 1+0.5A 2..................................(1)
Wherein A1 represents the response of the first acceleration transducer; A2 represents the response of the second acceleration transducer; T controlrepresent that translational degree of freedom controls magnitude.
For rotation, the first acceleration transducer 1, second acceleration transducer 2 needs direction of motion contrary, also can adopt mathematic(al) mean method to calculate, need consider the impact of acceleration of gravity constant G when linear acceleration being converted to rotary acceleration.That is:
R c o n t r o l = A 1 - A 2 2 r G ... ( 2 )
Convert upper two formula to matrix form, then have:
T c o n t r o l R c o n t r o l = 0.5 0.5 0.5 G r - 0.5 G r A 1 A 2 ... ( 3 )
Input coordinate transition matrix is expressed as follows:
[ T M ] i n p u t = 0.5 0.5 0.5 G r - 0.5 G r ... ( 4 )
Fig. 2 is Double earthquakes synchro control of the present invention excitation geometric layout figure.As shown in Figure 2, for output coordinate transition matrix, by the method for observing, output degree of freedom can be converted to linear incentive signal.Suppose the first shaking table excitation 4, second shaking table excitation 5 apart from the distance of two platform geometric center 3 of shaking be r ~.Utilize input transition matrix derivation method, can show that the relation that excitation degree of freedom and linear incentive export is as follows:
T d r i v e R d r i v e = 0.5 0.5 0.5 G r ~ - 0.5 G r ~ D 1 D 2 ... ( 5 )
Wherein D1 represents the drive singal of the first shaking table excitation 4; D2 represents the drive singal of the second shaking table excitation 5; T driverepresent that translational degree of freedom drives magnitude; R driverepresent that rotational freedom drives magnitude.
Therefore the pumping signal of reality is by obtaining above formula matrix inversion, that is:
D 1 D 2 = 0.5 0.5 0.5 G r ~ - 0.5 G r ~ - 1 T d r i v e R d r i v e ... ( 6 )
Finally show that output linearity pumping signal is as follows with the relation of excitation degree of freedom:
D 1 D 2 = 1 r ~ / G 1 - r ~ / G T d r i v e R d r i v e ... ( 7 )
Output coordinate expression matrix is as follows:
[ T M ] o u t p u t = 1 r ~ / G 1 - r ~ / G ... ( 8 )
Double earthquakes Synchronous motion control method based on coordinate conversion matrix comprises the steps:
(1) input coordinate transition matrix is built.Particularly, build input coordinate transition matrix to comprise according to the distance determination input coordinate transition matrix (see formula 4) (as shown in Figure 3) on Double earthquakes between two control sensors.Linear measurement Double earthquakes controlling sensor is converted to corresponding translation and rotation two degree of freedom by acting as of input coordinate transition matrix.
(2) output coordinate transition matrix is built.Particularly, build output coordinate transition matrix to comprise according to Double earthquakes excitation between distance determination output coordinate transition matrix (see formula 8) (as shown in Figure 4).The signal act as multi-input multi-output control system exports of output coordinate transition matrix is converted to the actual drive singal inputing to Double earthquakes power amplifier.
(3) two degree-of-freedom motion synchro control.Particularly, two degree-of-freedom motion synchro control is the white noise random signal that multi-input multi-output control system first sends the low magnitude of two-way, and this white noise random signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; On Double earthquakes, the response of two measurement points enters multi-input multi-output control system after input coordinate transition matrix, and control system calculates the transport function of whole system according to input and output signal; According to the translational degree of freedom reference test magnitude arranged, rotational freedom reference test magnitude, (rotational freedom reference test magnitude requires for translational degree of freedom test magnitude 10 to multi-input multi-output control system subsequently -3doubly following), system transter calculates initial two-way output signal, this two-way output signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; Multi-input multi-output control system is according to two degrees of freedom response level and the translation of measuring and rotate two degrees of freedom reference test magnitude and compare the also error of calculation, the real-time two-way drive singal revising transmission, thus realize Double earthquakes synchronized movement control (as shown in Figure 5).
Although give detailed description and explanation to the specific embodiment of the present invention above; but what should indicate is; we can carry out various equivalence according to conception of the present invention to above-mentioned embodiment and change and amendment; its function produced do not exceed that instructions and accompanying drawing contain yet spiritual time, all should within protection scope of the present invention.

Claims (5)

1., based on a Double earthquakes Synchronous motion control method for coordinate conversion matrix, comprise following three steps:
(1) input coordinate transition matrix is built, described structure input coordinate transition matrix step comprises according to the distance determination input coordinate transition matrix on Double earthquakes between two control sensors, and linear measurement Double earthquakes controlling sensor is converted to corresponding translation and rotation two degree of freedom by acting as of input coordinate transition matrix;
(2) output coordinate transition matrix is built, described structure output coordinate transition matrix step comprise according to Double earthquakes excitation between distance determination output coordinate transition matrix, output coordinate transition matrix act as by multi-input multi-output control system export signal be converted to the actual drive singal inputing to Double earthquakes power amplifier;
(3) two degree-of-freedom motion synchro control, described two degree-of-freedom motion synchro control step is specially the white noise random signal that multi-input multi-output control system first sends the low magnitude of two-way, and this white noise random signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; On Double earthquakes, the response of two measurement points enters multi-input multi-output control system after input coordinate transition matrix, and control system calculates the transport function of whole system according to input and output signal; Multi-input multi-output control system calculates initial two-way output signal according to the translational degree of freedom reference test magnitude arranged, rotational freedom reference test magnitude, system transter subsequently, and this two-way output signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; Multi-input multi-output control system is according to two degrees of freedom response level and the translation of measuring and rotates two degrees of freedom reference test magnitude and compare and the error of calculation, revises the two-way drive singal of transmission in real time, thus realizes the control of Double earthquakes synchronized movement.
2., as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, described input coordinate transition matrix and described output coordinate transition matrix are the square formation of 2 × 2.
3. as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, rotational freedom test magnitude requires for translational degree of freedom test magnitude 10 -3doubly.
4., as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, described input coordinate transition matrix is:
[ TM ] input = 0.5 G r 0.5 - 0.5 G r 0.5 , Wherein, G is acceleration of gravity constant, and r is the distance of the first acceleration transducer, the second acceleration transducer and Double earthquakes geometric center.
5., as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, described output coordinate transition matrix is:
[ TM ] output = 1 r ~ / G 1 - r ~ / G , Wherein, G is acceleration of gravity constant, be that the first shaking table encourages, the second shaking table encourages the distance with two platform geometric center of shaking.
CN201510282045.7A 2015-05-28 2015-05-28 Double earthquakes Synchronous motion control method based on coordinate conversion matrix Expired - Fee Related CN105320057B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510282045.7A CN105320057B (en) 2015-05-28 2015-05-28 Double earthquakes Synchronous motion control method based on coordinate conversion matrix

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510282045.7A CN105320057B (en) 2015-05-28 2015-05-28 Double earthquakes Synchronous motion control method based on coordinate conversion matrix

Publications (2)

Publication Number Publication Date
CN105320057A true CN105320057A (en) 2016-02-10
CN105320057B CN105320057B (en) 2018-04-06

Family

ID=55247657

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510282045.7A Expired - Fee Related CN105320057B (en) 2015-05-28 2015-05-28 Double earthquakes Synchronous motion control method based on coordinate conversion matrix

Country Status (1)

Country Link
CN (1) CN105320057B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106248327A (en) * 2016-07-11 2016-12-21 北京强度环境研究所 A kind of firing box Transport Vibration assay device and method
CN107543671A (en) * 2016-07-15 2018-01-05 北京卫星环境工程研究所 More and exciting movement controller and digit phase isochronous controller matching process
CN108254145A (en) * 2017-12-29 2018-07-06 苏州东菱智能减振降噪技术有限公司 A kind of control method for realizing multi shaker synchronous vibration
CN112444367A (en) * 2020-12-18 2021-03-05 中国工程物理研究院总体工程研究所 Multi-vibration-table parallel-pushing single-shaft vibration test control method
CN112629833A (en) * 2019-09-25 2021-04-09 上海汽车集团股份有限公司 Load acquisition method and device
CN114646440A (en) * 2022-03-16 2022-06-21 北京卫星环境工程研究所 Three-axis six-degree-of-freedom vibration control input and output coordinate transformation matrix construction method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005121371A (en) * 2003-10-14 2005-05-12 National Institute Of Advanced Industrial & Technology Dynamic matrix sensitivity measuring method for force sensor, and instrument therefor
CN101173876A (en) * 2007-11-30 2008-05-07 哈尔滨工业大学 Control method for three-axis six-freedom hydraulic vibration table based on kinematics
CN101196429A (en) * 2007-12-07 2008-06-11 哈尔滨工业大学 Internal force decoupling control method for driving vibration table with redundant 6 freedom
CN104133479A (en) * 2014-08-08 2014-11-05 上海新跃仪表厂 Test system and method for simulating flexible satellite three-axis attitude coupling movement with single-axis air bearing table

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005121371A (en) * 2003-10-14 2005-05-12 National Institute Of Advanced Industrial & Technology Dynamic matrix sensitivity measuring method for force sensor, and instrument therefor
CN101173876A (en) * 2007-11-30 2008-05-07 哈尔滨工业大学 Control method for three-axis six-freedom hydraulic vibration table based on kinematics
CN101196429A (en) * 2007-12-07 2008-06-11 哈尔滨工业大学 Internal force decoupling control method for driving vibration table with redundant 6 freedom
CN104133479A (en) * 2014-08-08 2014-11-05 上海新跃仪表厂 Test system and method for simulating flexible satellite three-axis attitude coupling movement with single-axis air bearing table

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHEN G,ZHUANG H,DU X等: "Analysis of large-scale shaking table test of dynamic soil-subway station interaction", 《EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION》 *
樊世超等: "多维动力学环境模拟试验技术研究", 《航天器环境工程》 *
邱汉平等: "双振动台正弦振动同步控制新方法研究", 《航天器环境工程》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106248327A (en) * 2016-07-11 2016-12-21 北京强度环境研究所 A kind of firing box Transport Vibration assay device and method
CN107543671A (en) * 2016-07-15 2018-01-05 北京卫星环境工程研究所 More and exciting movement controller and digit phase isochronous controller matching process
CN107543671B (en) * 2016-07-15 2019-02-19 北京卫星环境工程研究所 More and exciting movement controller and digit phase isochronous controller matching process
CN108254145A (en) * 2017-12-29 2018-07-06 苏州东菱智能减振降噪技术有限公司 A kind of control method for realizing multi shaker synchronous vibration
CN108254145B (en) * 2017-12-29 2020-06-26 苏州东菱智能减振降噪技术有限公司 Control method for realizing synchronous vibration of multiple vibration tables
CN112629833A (en) * 2019-09-25 2021-04-09 上海汽车集团股份有限公司 Load acquisition method and device
CN112629833B (en) * 2019-09-25 2024-04-19 上海汽车集团股份有限公司 Load acquisition method and device
CN112444367A (en) * 2020-12-18 2021-03-05 中国工程物理研究院总体工程研究所 Multi-vibration-table parallel-pushing single-shaft vibration test control method
CN114646440A (en) * 2022-03-16 2022-06-21 北京卫星环境工程研究所 Three-axis six-degree-of-freedom vibration control input and output coordinate transformation matrix construction method
CN114646440B (en) * 2022-03-16 2023-06-20 北京卫星环境工程研究所 Three-axis six-degree-of-freedom vibration control input/output coordinate transformation matrix construction method

Also Published As

Publication number Publication date
CN105320057B (en) 2018-04-06

Similar Documents

Publication Publication Date Title
CN105320057A (en) Synchronous movement control method of double vibration tables on the basis of coordinate-transformation matrix
CN105196294B (en) Reconfigurable mechanical arm decentralized control system and control method adopting position measuring
CN103558079B (en) Based on the multiple degrees of freedom loading method of parallel institution driving force closed loop
US20180031457A1 (en) Shear test device and test method of rock mass discontinuities under constant normal stiffness condition
Kang Closed-form force sensing of a 6-axis force transducer based on the Stewart platform
CN104006964B (en) Harmonic drive tooth mesh and deformation of flexible wheel dynamic characterization measurement system
JPS61501413A (en) Load dynamics compensation circuit for servo hydraulic control system
CN101763038A (en) Method for controlling structural modal vibration of dual-frame magnetic levitation control moment gyroscope
CN105015805B (en) A kind of secondary plateform system
CN104407639A (en) Servo control method and servo control device for high-frequency hydraulic shock excitation system
CN104807606B (en) A kind of MATLAB STM32 hybrid powers pilot system and its test method
CN106052986A (en) Three-dimensional six-degree-of-freedom force-displacement hybrid control method
CN108241289B (en) Force-displacement hybrid control method for three-degree-of-freedom loading system
CN206348142U (en) Double-wall pipe elastic support dynamic fatigue test device
CN202520846U (en) Low-frequency vibration isolator aiming at foundation excitation
CN109668704B (en) Separation type hybrid test system and test method
CN106202735A (en) A kind of local nonlinearity foundation soil structural interaction sub-structural test method
CN105508823A (en) Motion control method for flexible six-degree-of-freedom parallel precision adjustment device
CN103344506B (en) A kind of structural face shear strength size effect Combined Trials method
CN109709822B (en) Whole-satellite ground semi-physical test device for large-scale flexible spacecraft
CN102540903A (en) Modelica-language-based simulation modeling method for pump truck boom system
CN209820722U (en) Separated mixing test system
Jiang et al. Analysis of dynamic characteristics of water hydraulic rotating angle self-servo robot joint actuator
Feng et al. A novel adaptive balance-drive mechanism for industrial robots using a series elastic actuator
CN107830980A (en) A kind of unidirectional shear model box of hollow variable cross-section

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

Granted publication date: 20180406

Termination date: 20210528

CF01 Termination of patent right due to non-payment of annual fee