CN114397074B - Sinusoidal vibration table control method, system and device - Google Patents

Sinusoidal vibration table control method, system and device Download PDF

Info

Publication number
CN114397074B
CN114397074B CN202210068267.9A CN202210068267A CN114397074B CN 114397074 B CN114397074 B CN 114397074B CN 202210068267 A CN202210068267 A CN 202210068267A CN 114397074 B CN114397074 B CN 114397074B
Authority
CN
China
Prior art keywords
dead zone
module
displacement
signal
vibration table
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
Application number
CN202210068267.9A
Other languages
Chinese (zh)
Other versions
CN114397074A (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.)
Sun Yat Sen University
Sun Yat Sen University Shenzhen Campus
Original Assignee
Sun Yat Sen University
Sun Yat Sen University Shenzhen Campus
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 Sun Yat Sen University, Sun Yat Sen University Shenzhen Campus filed Critical Sun Yat Sen University
Priority to CN202210068267.9A priority Critical patent/CN114397074B/en
Publication of CN114397074A publication Critical patent/CN114397074A/en
Application granted granted Critical
Publication of CN114397074B publication Critical patent/CN114397074B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G01M7/022Vibration control arrangements, e.g. for generating random vibrations

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a sinusoidal vibration table control method, a sinusoidal vibration table control system and a sinusoidal vibration table control device, wherein the sinusoidal vibration table control method comprises the following steps: measuring the displacement of the vibrating table; calculating a vibration table top displacement feedback signal according to the displacement measured value; introducing a dead zone relay module, and inputting a table surface displacement feedback signal for calculating vibration to the dead zone relay module; introducing a low-pass filter, and inputting the output signal of the dead zone relay module to the low-pass filter; measuring the displacement of the actuator and calculating a deviation signal in combination with the filtered feedback signal; and introducing a PID controller, and calculating by combining the deviation signal to obtain an actuator control signal. The system comprises: the device comprises a displacement measurement module, a displacement feedback calculation module, a dead zone relay module, a low-pass filter introduction module, a deviation calculation module and a control calculation module. The apparatus includes a memory and a processor for performing the sinusoidal vibration table control method described above. By using the invention, the distortion degree of the output waveform of the vibrating table is reduced. The invention can be widely applied to the field of vibration control.

Description

Sinusoidal vibration table control method, system and device
Technical Field
The invention relates to the field of vibration control, in particular to a sinusoidal vibration table control method, a sinusoidal vibration table control system and a sinusoidal vibration table control device.
Background
The vibration environment simulation test is to examine the reliability, stability and other indexes of the product in the environment by simulating the vibration environment experienced by the product in the transportation and use processes in the service life period under the laboratory condition, and is widely applied to the industrial fields and departments of aviation, aerospace, weapons, ships, automobiles, buildings and the like at present. Sinusoidal vibration is a typical vibration form, and various periodic vibrations can be studied by using the form of sinusoidal vibration, and by performing vibration test on a product, the reliability thereof can be checked by analyzing test data, and the exposed product defects thereof can be improved.
Most of the existing sinusoidal vibration control algorithms send out main control signals by a computer, and after the feedback signals are correctly measured and compared, deviation signals are obtained and the system is adjusted by a controller. However, there are still some problems that the output signal of the vibrating table has a certain amplitude error and phase lag with the input signal, and the output signal often contains higher harmonics, which results in serious distortion of the output waveform.
Disclosure of Invention
In order to solve the technical problems, the invention aims to provide a sinusoidal vibration table control method, a sinusoidal vibration table control system and a sinusoidal vibration table control device, which can reduce the distortion degree of the output waveform of the vibration table.
The first technical scheme adopted by the invention is as follows: a sinusoidal vibration table control method comprising the steps of:
measuring the displacement of the vibrating table to obtain a displacement measurement value;
calculating a vibration table top displacement feedback signal according to the displacement measured value;
introducing a nonlinear module-dead zone relay module, and inputting a vibration table surface displacement feedback signal to the nonlinear module-dead zone relay module to obtain an output signal of the dead zone relay module;
introducing a low-pass filter, inputting the output signal of the dead zone relay module to the low-pass filter, and outputting a filtering feedback signal;
measuring the displacement of the actuator and calculating a deviation signal in combination with the filtered feedback signal;
and introducing a PID controller, and calculating by combining the deviation signal to obtain an actuator control signal.
Further, the calculation formula of the vibration table displacement feedback signal is as follows:
in the above, y i Representing displacement measurements, n represents the number of sensors that acquire the displacement measurements.
Further, the formula of the nonlinear module-dead zone relay module is as follows:
in the above, a (t) represents a vibration table displacement feedback signal, a s (t) represents a dead zone relay characteristic module output signal, h represents a dead zone width, and M represents a relay output upper limit.
Further, the low-pass filter adopts a first-order low-pass filter, and the formula is as follows:
in the above equation, s represents an operator generated by the laplace transform, and ω represents a filter cut-off frequency.
Further, the calculation formula of the deviation signal is as follows:
e=-x 0 -x
in the above, x 0 Representing the filtered feedback signal, x representing the actuator displacement, e representing the deviation signal.
Further, the calculation formula of the control signal is as follows:
in the above, K i Represent the integral constant, K p Representing the proportionality coefficient, K d Representing the differential constant.
The second technical scheme adopted by the invention is as follows: a sinusoidal vibration table control system, comprising:
the displacement measuring module is used for measuring the displacement of the vibrating table to obtain a displacement measuring value;
the displacement feedback calculation module is used for calculating a vibration table top displacement feedback signal according to the displacement measured value;
the dead zone relay module is used for introducing the nonlinear module-dead zone relay module and inputting the vibration table surface displacement feedback signal to the nonlinear module-dead zone relay module to obtain an output signal of the dead zone relay module;
the low-pass filter introduction module is used for introducing the low-pass filter, inputting the output signal of the dead zone relay module to the low-pass filter and outputting a filtering feedback signal;
the deviation calculation module is used for measuring the displacement of the actuator and calculating a deviation signal by combining the filtered feedback signal;
and the control calculation module is used for introducing a PID controller and calculating to obtain an actuator control signal by combining the deviation signal.
The third technical scheme adopted by the invention is as follows: a sinusoidal vibration table control system, comprising:
at least one processor;
at least one memory for storing at least one program;
the at least one program, when executed by the at least one processor, causes the at least one processor to implement a sinusoidal vibration table control method as described above.
The method, the system and the device have the beneficial effects that: according to the invention, a computer master control signal generating device is not needed, the self-oscillation characteristic of the system is utilized, the sinusoidal vibration table output signal control is performed by generating self-oscillation by the system based on a nonlinear stability analysis method, the amplitude error and the phase lag of the output waveform are avoided, and the distortion degree of the output waveform is reduced.
Drawings
FIG. 1 is a flow chart of the steps of a sinusoidal vibration table control method of the present invention;
FIG. 2 is a schematic diagram of a flow framework in accordance with an embodiment of the present invention;
fig. 3 is a block diagram of a sinusoidal vibration table control system of the present invention.
Detailed Description
The invention will now be described in further detail with reference to the drawings and to specific examples. The step numbers in the following embodiments are set for convenience of illustration only, and the order between the steps is not limited in any way, and the execution order of the steps in the embodiments may be adaptively adjusted according to the understanding of those skilled in the art.
Referring to fig. 1 and 2, the present invention provides a sinusoidal vibration table control method comprising the steps of:
s1, measuring displacement of a vibrating table to obtain a displacement measurement value;
specifically, the table displacement y is measured by displacement sensors distributed over the table 1 ,y 2 ,...,y n The number n of the sensors is determined by the configuration of the vibrating table, the number of the vibrating table is 1 under the conditions that the vibrating table surface and the test piece are smaller and the table surface uniformity is better, and the number of the sensors is two or more under the conditions that the vibrating table surface and the test piece are larger and the table surface uniformity is worse.
S2, calculating a vibration table top displacement feedback signal according to the displacement measured value;
specifically, the calculation formula of the vibration table displacement feedback signal is as follows by using a multipoint average control method:
in the above, y i Representing displacement measurements, n represents the number of sensors that acquire the displacement measurements.
S3, introducing a nonlinear module-dead zone relay module, and inputting a vibration table surface displacement feedback signal to the nonlinear module-dead zone relay module to obtain an output signal of the dead zone relay module;
specifically, the formula of the nonlinear module-dead zone relay module is as follows:
in the above, a (t) represents a vibration table displacement feedback signal, a s (t) represents the dead band relay characteristic module output signal,h represents the dead zone width, and M represents the relay output upper limit.
The dead zone relay module is introduced in the step, and the input signal of the dead zone relay module can enable the dead zone relay module to have output only when the dead zone width is crossed, so that the dead zone width h has a threshold value, and the system can be primarily judged to have output. In addition, the relay output upper limit M jointly determines whether the system has stable output under the condition of different dead zone widths h, and the threshold value of the relay output upper limit M changes along with the different dead zone widths h. Therefore, if the system needs to judge whether stable output exists, the dead zone width h and the threshold value of the relay output upper limit M are required to be solved, and the solving method is as follows:
first, the Nyquist diagram of the open-loop transfer function G(s) is plotted to find the intersection point H (X) with the real axis H ,0);
Secondly, drawing a dead zone relay module negative inverse description function image, wherein the dead zone relay module description function is as followsA is the amplitude of the input signal of the dead zone relay characteristic module;
and if the system does not have stable self-oscillation, the dead zone relay characteristic module input signal a (t) is a sine attenuation signal, and the attenuation process |a (t) | has a maximum value |a (t) | max The specific numerical value is related to the initial state of the system and the environment;
again, the dead zone width h should be greater than |a (t) | max I.e. dead zone width threshold h 0 =|a(t)| max
Again, at h > |a (t) | max Under the condition that the system has stable self-oscillation, the condition is satisfiedI.e. relay output upper threshold +.>
Under the condition that the dead zone width threshold value and the relay output upper limit threshold value are obtained, the dead zone width h and the relay output upper limit M are corrected, so that the system has stable self-oscillation.
When the system has stable self-oscillation motion, on the premise of not changing the value of the dead zone width h, if the vibration amplitude of the actuator is required to be changed, only the upper limit M of the relay is required to be changed.
S4, introducing a low-pass filter, inputting the output signal of the dead zone relay module to the low-pass filter, and outputting a filtering feedback signal;
specifically, the low-pass filter adopts a first-order low-pass filter, and the formula is as follows:s denotes an operator generated by the laplace transform, ω denotes the filter cut-off frequency.
S5, measuring the displacement of the actuator and calculating a deviation signal by combining the filtered feedback signal;
specifically, the calculation formula of the deviation signal: e= -x 0 -x,x 0 Representing the filtered feedback signal, x representing the actuator displacement, e representing the deviation signal.
S6, introducing a PID controller, calculating by combining the deviation signal to obtain an actuator control signal, and enabling the actuator to displace and output force after obtaining the control signal so as to push the vibrating table to displace.
Specifically, the calculation formula of the control signal is as follows:K i represent the integral constant, K p Representing the proportionality coefficient, K d Representing the differential constant.
Aiming at the problems of amplitude error, phase lag, serious distortion degree of output waveform and the like in the application process of the conventional sinusoidal vibration table control algorithm, the invention removes a computer main control signal generating device by designing a novel sinusoidal vibration table control algorithm based on a nonlinear stability analysis method, completes sinusoidal vibration signal output by utilizing the self-oscillation characteristic of a system, reduces state feedback variables and reduces the distortion degree of the output waveform of the vibration table.
As shown in fig. 3, a sinusoidal vibration table control system, comprising:
the displacement measuring module is used for measuring the displacement of the vibrating table to obtain a displacement measuring value;
the displacement feedback calculation module is used for calculating a vibration table top displacement feedback signal according to the displacement measured value;
the dead zone relay module is used for introducing the nonlinear module-dead zone relay module and inputting the vibration table surface displacement feedback signal to the nonlinear module-dead zone relay module to obtain an output signal of the dead zone relay module;
the low-pass filter introduction module is used for introducing the low-pass filter, inputting the output signal of the dead zone relay module to the low-pass filter and outputting a filtering feedback signal;
the deviation calculation module is used for measuring the displacement of the actuator and calculating a deviation signal by combining the filtered feedback signal;
and the control calculation module is used for introducing a PID controller and calculating to obtain an actuator control signal by combining the deviation signal.
A sinusoidal vibration table control device:
at least one processor;
at least one memory for storing at least one program;
the at least one program, when executed by the at least one processor, causes the at least one processor to implement a sinusoidal vibration table control method as described above.
The content in the method embodiment is applicable to the embodiment of the device, and the functions specifically realized by the embodiment of the device are the same as those of the method embodiment, and the obtained beneficial effects are the same as those of the method embodiment.
While the preferred embodiment of the present invention has been described in detail, the invention is not limited to the embodiment, and various equivalent modifications and substitutions can be made by those skilled in the art without departing from the spirit of the invention, and these modifications and substitutions are intended to be included in the scope of the present invention as defined in the appended claims.

Claims (7)

1. A sinusoidal vibration table control method, characterized by comprising the steps of:
measuring the displacement of the vibrating table to obtain a displacement measurement value;
calculating a vibration table top displacement feedback signal according to the displacement measured value;
introducing a nonlinear module-dead zone relay module, and inputting a vibration table surface displacement feedback signal to the nonlinear module-dead zone relay module to obtain an output signal of the dead zone relay module;
introducing a low-pass filter, inputting the output signal of the dead zone relay module to the low-pass filter, and outputting a filtering feedback signal;
measuring the displacement of the actuator and calculating a deviation signal in combination with the filtered feedback signal;
introducing a PID controller, and calculating by combining the deviation signal to obtain an actuator control signal;
the formula of the nonlinear module-dead zone relay module is expressed as follows:
in the above, a (t) represents a vibration table displacement feedback signal, a s (t) represents a dead zone relay characteristic module output signal, h represents a dead zone width, and M represents a relay output upper limit;
the solving method of h and M is as follows:
first, the Nyquist diagram of the open-loop transfer function G(s) is plotted to find the intersection point H (X) with the real axis H ,0);
Secondly, drawing a dead zone relay module negative inverse description function image, wherein the dead zone relay module description function is as followsA is the amplitude of the input signal of the dead zone relay characteristic module;
and if the system does not have stable self-oscillation, the dead zone relay characteristic module input signal a (t) is a sine attenuation signal, and the attenuation process |a (t) | has a maximum value |a (t) | max
Again, the dead zone width h should be greater than |a (t) | max I.e. dead zone width threshold h 0 =|a(t)| max
Again, at h > |a (t) | max Under the condition that the system has stable self-oscillation, the condition is satisfiedI.e. relay output upper threshold +.>
2. The sinusoidal vibration table control method of claim 1, wherein the vibration table displacement feedback signal is calculated as follows:
in the above, y i Representing displacement measurements, n represents the number of sensors that acquire the displacement measurements.
3. The method for controlling a sinusoidal vibration table according to claim 2, wherein the low-pass filter is a first-order low-pass filter, and the formula is as follows:
in the above equation, s represents an operator generated by the laplace transform, and ω represents a filter cut-off frequency.
4. A sinusoidal vibration table control method according to claim 3, characterized in that the calculation formula of the deviation signal is as follows:
e=-x 0 -x
in the above, x 0 Representing the filtered feedback signal, x representing the actuator displacement, e representing the deviation signal.
5. The sinusoidal vibration table control method of claim 4, wherein the control signal is calculated as:
in the above, K i Represent the integral constant, K p Representing the proportionality coefficient, K d Representing the differential constant.
6. A sinusoidal vibration table control system, comprising:
the displacement measuring module is used for measuring the displacement of the vibrating table to obtain a displacement measuring value;
the displacement feedback calculation module is used for calculating a vibration table top displacement feedback signal according to the displacement measured value;
the dead zone relay module is used for introducing the nonlinear module-dead zone relay module and inputting the vibration table surface displacement feedback signal to the nonlinear module-dead zone relay module to obtain an output signal of the dead zone relay module;
the low-pass filter introduction module is used for introducing the low-pass filter, inputting the output signal of the dead zone relay module to the low-pass filter and outputting a filtering feedback signal;
the deviation calculation module is used for measuring the displacement of the actuator and calculating a deviation signal by combining the filtered feedback signal;
the control calculation module is used for introducing a PID controller and calculating to obtain an actuator control signal by combining the deviation signal;
the formula of the nonlinear module-dead zone relay module is expressed as follows:
on the upper partWherein a (t) represents a vibration table displacement feedback signal, a s (t) represents a dead zone relay characteristic module output signal, h represents a dead zone width, and M represents a relay output upper limit;
the solving method of h and M is as follows:
first, the Nyquist diagram of the open-loop transfer function G(s) is plotted to find the intersection point H (X) with the real axis H ,0);
Secondly, drawing a dead zone relay module negative inverse description function image, wherein the dead zone relay module description function is as followsA is the amplitude of the input signal of the dead zone relay characteristic module;
and if the system does not have stable self-oscillation, the dead zone relay characteristic module input signal a (t) is a sine attenuation signal, and the attenuation process |a (t) | has a maximum value |a (t) | max
Again, the dead zone width h should be greater than |a (t) | max I.e. dead zone width threshold h 0 =|a(t)| max
Again, at h > |a (t) | max Under the condition that the system has stable self-oscillation, the condition is satisfiedI.e. relay output upper threshold +.>
7. A sinusoidal vibration table control device, comprising:
at least one processor;
at least one memory for storing at least one program;
when the at least one program is executed by the at least one processor, the at least one processor is caused to implement a sinusoidal vibration table control method as claimed in any one of claims 1-5.
CN202210068267.9A 2022-01-20 2022-01-20 Sinusoidal vibration table control method, system and device Active CN114397074B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210068267.9A CN114397074B (en) 2022-01-20 2022-01-20 Sinusoidal vibration table control method, system and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210068267.9A CN114397074B (en) 2022-01-20 2022-01-20 Sinusoidal vibration table control method, system and device

Publications (2)

Publication Number Publication Date
CN114397074A CN114397074A (en) 2022-04-26
CN114397074B true CN114397074B (en) 2024-02-13

Family

ID=81233286

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210068267.9A Active CN114397074B (en) 2022-01-20 2022-01-20 Sinusoidal vibration table control method, system and device

Country Status (1)

Country Link
CN (1) CN114397074B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE8104989L (en) * 1981-08-24 1983-02-25 Naf Ab PROCEDURE FOR SETTING A PID REGULATOR FOR A PROCESS
CA2077159A1 (en) * 1991-08-30 1993-03-01 Sheldon Gilbert Lloyd Tuning arrangement for process controllers
WO1996034323A1 (en) * 1995-04-28 1996-10-31 Centre National De La Recherche Scientifique Method and device for adjusting a pid controller
RU2475797C1 (en) * 2011-09-20 2013-02-20 Владимир Романович Сабанин Monitor unit for extremal controller
CN106383523A (en) * 2016-09-30 2017-02-08 湖北航天技术研究院总体设计所 Aircraft nonlinear attitude control system stability analysis method
JP2019168777A (en) * 2018-03-22 2019-10-03 三菱重工業株式会社 Control system design support device, control system design support method, and control system design support program
CN111687841A (en) * 2020-06-15 2020-09-22 中山大学 Robot bounce height control method, system, device and storage medium
JP2020190880A (en) * 2019-05-21 2020-11-26 国立研究開発法人防災科学技術研究所 Vibration table control apparatus and vibration table control method
CN112445234A (en) * 2020-11-27 2021-03-05 航天科工火箭技术有限公司 Attitude control method and device for spacecraft

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE8104989L (en) * 1981-08-24 1983-02-25 Naf Ab PROCEDURE FOR SETTING A PID REGULATOR FOR A PROCESS
US4549123A (en) * 1981-08-24 1985-10-22 Naf Controls Ab Method and an apparatus in tuning a PID-regulator
CA2077159A1 (en) * 1991-08-30 1993-03-01 Sheldon Gilbert Lloyd Tuning arrangement for process controllers
WO1996034323A1 (en) * 1995-04-28 1996-10-31 Centre National De La Recherche Scientifique Method and device for adjusting a pid controller
RU2475797C1 (en) * 2011-09-20 2013-02-20 Владимир Романович Сабанин Monitor unit for extremal controller
CN106383523A (en) * 2016-09-30 2017-02-08 湖北航天技术研究院总体设计所 Aircraft nonlinear attitude control system stability analysis method
JP2019168777A (en) * 2018-03-22 2019-10-03 三菱重工業株式会社 Control system design support device, control system design support method, and control system design support program
JP2020190880A (en) * 2019-05-21 2020-11-26 国立研究開発法人防災科学技術研究所 Vibration table control apparatus and vibration table control method
CN111687841A (en) * 2020-06-15 2020-09-22 中山大学 Robot bounce height control method, system, device and storage medium
CN112445234A (en) * 2020-11-27 2021-03-05 航天科工火箭技术有限公司 Attitude control method and device for spacecraft

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于相角裕度的PID 控制器自整定改进算法的研究;郑艳秋等;《核动力工程》;第41卷(第增刊2期);正文第1-3部分 *

Also Published As

Publication number Publication date
CN114397074A (en) 2022-04-26

Similar Documents

Publication Publication Date Title
CN103197542B (en) Time lag system PID controller based on data-driven is calmed method
CN106444372A (en) Sliding mode repetitive controller for motor servo system
Ozana et al. PID controller design based on global optimization technique with additional constraints
CN104199993A (en) Method for designing two-dimension dynamic Kalman filter for intermittent process
CN105629766B (en) Multivariable time delay system discrimination method based on step test
KR20080043941A (en) Integrals of relay feedback responses for extracting process information
CN114397074B (en) Sinusoidal vibration table control method, system and device
Zhou et al. Robust adaptive flexible prescribed performance tracking and vibration control for rigid–flexible coupled robotic systems with input quantization
CN108845594B (en) Vibration harmonic iterative control method based on transfer function
EP4339621A1 (en) Mems accelerometer systems
CN108802153A (en) A kind of calibration method of periodic attenuation, electronic equipment and storage medium
CN115562019B (en) Tower crane periodic sliding mode control method and system with unmatched interference
CN114310874B (en) Cooperative robot joint moment control method, system and computer equipment
CN113742967B (en) Intermittent finite element artificial viscous shock wave processing method based on strong residual error
Gerov et al. Parameter estimation method for the unstable time delay process
Dovetta et al. Uncertainty propagation in model extraction by system identification and its implication for control design
Eek et al. Closed‐loop identification of a continuous crystallization process
CN109977493B (en) Stability determining method of absolute value system based on double sine input description function
Koehler et al. Determining the fidelity of hardware-in-the-loop simulation coupling systems
CN110967045A (en) Temperature effect error compensation system and design method of optical fiber sensor
Balasevicius et al. Design of the deadbeat controller with limited output
Soudani et al. On an internal multimodel control for non-linear systems: a comparative study
JP2013161206A (en) Feedback control device having reset windup countermeasure
Jia et al. A novel polynomial tracking differentiator
Apkarian et al. Robust control of systems with hyperbolic partial differential equations

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