CN108928476B - Differential calculation method for pitch angle of helicopter - Google Patents

Differential calculation method for pitch angle of helicopter Download PDF

Info

Publication number
CN108928476B
CN108928476B CN201810678996.XA CN201810678996A CN108928476B CN 108928476 B CN108928476 B CN 108928476B CN 201810678996 A CN201810678996 A CN 201810678996A CN 108928476 B CN108928476 B CN 108928476B
Authority
CN
China
Prior art keywords
axis
angle
angle sensor
straight line
pitch angle
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
CN201810678996.XA
Other languages
Chinese (zh)
Other versions
CN108928476A (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.)
China Helicopter Research and Development Institute
Original Assignee
China Helicopter Research and Development Institute
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 China Helicopter Research and Development Institute filed Critical China Helicopter Research and Development Institute
Priority to CN201810678996.XA priority Critical patent/CN108928476B/en
Publication of CN108928476A publication Critical patent/CN108928476A/en
Application granted granted Critical
Publication of CN108928476B publication Critical patent/CN108928476B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64CAEROPLANES; HELICOPTERS
    • B64C27/00Rotorcraft; Rotors peculiar thereto
    • B64C27/54Mechanisms for controlling blade adjustment or movement relative to rotor head, e.g. lag-lead movement
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60Testing or inspecting aircraft components or systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Toys (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

The invention relates to the field of helicopter pitch angle measurement, in particular to a differential calculation method for a helicopter pitch angle, which comprises the following steps: disposing a first angle sensor and a second angle sensor on the blade mounting surface a1 and the centerpiece plane B1, respectively; establishing a rectangular coordinate system, and establishing and simplifying a model; by controlling the three-axis rotation directions of the blade mounting surface A1 and the center piece plane B1 in the rectangular coordinate system X, Y, Z, the condition that the main rotor wing axis is not vertical is simulated, and the pitch angle of the corresponding condition is calculated. The differential calculation method for the pitch angle of the helicopter can realize the measurement of the pitch angle under the condition that the main rotor shaft is not perpendicular, has stable precision and high efficiency compared with the traditional method, and can also ensure the personal safety of measuring personnel and reduce the labor intensity.

Description

Differential calculation method for pitch angle of helicopter
Technical Field
The invention relates to the field of helicopter pitch angle measurement, in particular to a differential calculation method for a helicopter pitch angle.
Background
After a helicopter flight control system installs or replaces components (such as a booster, a steering engine, and the like) on the helicopter, the system needs to be adjusted and checked to determine whether the position of a specific control device (a steering column, a collective lever, and a pedal) of the helicopter corresponds to a specific blade installation angle. Especially the blade mounting angle, the accuracy of the measurement is directly related to the flight performance and the flight quality of the helicopter.
Conventional blade setting angle measurements are typically made using a level gauge, by which the angle between the blade mounting surface at the root of the blade and the plane of the centerpiece is measured, while the helicopter main rotor shaft is "plumbed". According to the spatial geometrical knowledge, the included angle of the spatial plane is equal to the included angle of a straight line which is respectively positioned in two planes and is perpendicular to the intersection line of the two planes, the angle measured by the level is the included angle between the straight line and the horizontal plane, and on the basis, the measurement principle of the traditional method is that the angle of the straight line which is perpendicular to the intersection line on the blade installation surface, namely the pitch angle, is measured on the basis of ensuring the plane of the central component to be horizontal. The calculation needs to ensure that the plane of the central component is horizontal, and the central component needs to be vertical before measurement, so that the method has the disadvantages of multiple guarantee conditions, low efficiency and poor safety.
Disclosure of Invention
The invention aims to design a differential calculation method for a pitch angle of a helicopter, which is used for solving at least one problem existing in the existing blade installation angle measurement method.
The technical scheme of the invention is as follows:
a differential calculation method for a helicopter pitch angle comprises the following steps:
the method comprises the following steps of firstly, respectively arranging a first angle sensor and a second angle sensor on a blade mounting surface A1 and a centerpiece plane B1, wherein the X-axis direction of the first angle sensor and the X-axis direction of the second angle sensor are parallel to the axis of a blade, and the Y-axis of the first angle sensor and the Y-axis of the second angle sensor are in the same direction;
establishing a rectangular coordinate system, and establishing and simplifying a model according to the structural form of the step one, wherein the axis direction of the blade is the X-axis direction of the rectangular coordinate system, a straight line Y1 represents the Y-axis of the first angle sensor, the included angle between the straight line Y1 and the horizontal plane is alpha, a straight line Y2 represents the Y-axis of the second angle sensor, and the included angle between the straight line Y2 and the horizontal plane is beta;
thirdly, simulating the condition that the main rotor shaft is not vertical by controlling the three-axis rotation directions of the blade mounting surface A1 and the center piece plane B1 in the rectangular coordinate system X, Y, Z, and calculating the pitch angle of the corresponding condition, wherein the conditions comprise the following conditions:
when the blade mounting surface A1 and the center piece plane B1 rotate around the X axis of the rectangular coordinate system, the pitch angle theta is the included angle between the Y axis of the first angle sensor and the Y axis of the second angle sensor;
when the blade mounting surface A1 and the center piece plane B1 rotate around the Z axis of the rectangular coordinate system, the pitch angle theta is the included angle between the Y axis of the first angle sensor and the Y axis of the second angle sensor;
when the blade mounting surface a1 and the center piece plane B1 rotate around the Y axis of the rectangular coordinate system by an angle γ, the line Y1 'is the Y axis of the first rotated angle sensor, the angle between the line Y1' and the horizontal plane is α ', the line Y2' is the Y axis of the second rotated angle sensor, the angle between the line Y2 'and the horizontal plane is β', and the pitch angle θ is calculated according to the following formula (1):
Figure GDA0003107762540000021
optionally, in the third step, the formula (1) is obtained from the following formulas (2) and (3):
Figure GDA0003107762540000031
Figure GDA0003107762540000032
wherein, formula (2) is obtained from formula (4) and formula (5):
Figure GDA0003107762540000033
Figure GDA0003107762540000034
wherein R is the length from a predetermined point on the straight line Y1 to the perpendicular line of the XY surface of the rectangular coordinate system, and L is the length from a point on the straight line Y1' corresponding to the position of the predetermined point on the straight line Y1 to the perpendicular line of the XY surface of the rectangular coordinate system;
from equation (2), equation (3) can be derived in the same way.
The invention has the following effects:
the differential calculation method for the pitch angle of the helicopter can realize the measurement of the pitch angle under the condition that the main rotor shaft is not vertically drilled, has stable precision and high efficiency compared with the traditional method, can also ensure the personal safety of measuring personnel, and reduces the labor intensity.
Drawings
FIG. 1 is a simplified model of pitch angle measurement in the differential calculation method of the pitch angle of the helicopter of the present invention;
FIG. 2 is a main rotor shaft inclination model 1 in the differential calculation method of the pitch angle of the helicopter of the present invention;
FIG. 3 is a main rotor shaft inclination model 2 in the differential calculation method of the pitch angle of the helicopter of the present invention.
Detailed Description
In order to make the implementation objects, technical solutions and advantages of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in more detail below with reference to the accompanying drawings in the embodiments of the present invention. In the drawings, the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The described embodiments are only some, but not all embodiments of the invention. The embodiments described below with reference to the drawings are illustrative and intended to be illustrative of the invention and are not to be construed as limiting the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention. Embodiments of the present invention will be described in detail below with reference to the accompanying drawings.
In the description of the present invention, it is to be understood that the terms "center", "longitudinal", "lateral", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc., indicate orientations or positional relationships based on those shown in the drawings, and are used merely for convenience in describing the present invention and for simplifying the description, and do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore, should not be taken as limiting the scope of the present invention.
The method for calculating the difference between the helicopter pitch angles according to the present invention will be described in further detail with reference to fig. 1 to 3.
The differential calculation method for the helicopter pitch angle respectively calculates the included angles between the plane of the central piece and the horizontal plane and between the blade mounting surface and the horizontal plane, and obtains the pitch angle by calculating the difference value of the two included angles. When the sensor is installed, the X axis is parallel to the axis of the blade, and the intersecting line of the plane of the central part and the blade installation surface is the axis of the blade, so that the pitch angle is converted into the included angle of the Y axis of the sensor. The data measured by the sensor is the included angle between the shaft and the horizontal plane, and when the main rotor shaft of the helicopter is not vertically drilled, the size of the pitch angle is not equal to the difference value of the angle of the Y shaft, and space calculation is needed.
Specifically, the differential calculation method for the pitch angle of the helicopter of the invention comprises the following steps:
the method comprises the following steps of firstly, respectively arranging a first angle sensor and a second angle sensor on a blade mounting surface A1 and a centerpiece plane B1, wherein the X-axis direction of the first angle sensor and the X-axis direction of the second angle sensor are parallel to the axis of a blade, and the Y-axis of the first angle sensor and the Y-axis of the second angle sensor are in the same direction.
And step two, establishing a rectangular coordinate system, and establishing and simplifying a model according to the structural form of the step one, wherein the axis direction of the blade is the X-axis direction of the rectangular coordinate system, a straight line Y1 represents the Y-axis of the first angle sensor, the included angle between the straight line Y1 and the horizontal plane is alpha, a straight line Y2 represents the Y-axis of the second angle sensor, and the included angle between the straight line Y2 and the horizontal plane is beta. As can be seen from the spatial geometry knowledge, the included angle theta between Y1 and Y2 is equal to the included angle between the A1 and B1 planes, namely the pitch angle.
Step three, the situation that the main rotor shaft is not vertically arranged can be represented as A, B surfaces rotating around the X, Y, Z shaft together; by controlling the three-axis rotation directions of the blade mounting surface A1 and the center piece plane B1 in the rectangular coordinate system X, Y, Z, the non-vertical condition of the main rotor shaft is simulated, and the pitch angle of the corresponding condition is calculated, including the following conditions:
when the blade mounting surface A1 and the center piece plane B1 rotate around the X axis of the rectangular coordinate system, the pitch angle theta is the included angle between the Y axis of the first angle sensor and the Y axis of the second angle sensor;
when the blade mounting surface A1 and the center piece plane B1 rotate around the Z axis of the rectangular coordinate system, the pitch angle theta is the included angle between the Y axis of the first angle sensor and the Y axis of the second angle sensor;
when the blade mounting surface a1 and the center piece plane B1 rotate by an angle γ around the Y-axis of the rectangular coordinate system, the straight line Y1 ' is the Y-axis of the first rotated angle sensor, the included angle between the straight line Y1 ' and the horizontal plane (XY plane) is α ', the straight line Y2 ' is the Y-axis of the second rotated angle sensor, the included angle between the straight line Y2 ' and the horizontal plane (XY plane) is β ', the straight lines Y1 and Y1 ' are located on a cone with the point O as a dot and the circle a as a bottom, and the pitch angle θ is calculated according to the following formula (1):
Figure GDA0003107762540000061
further, in the differential calculation method for the pitch angle of the helicopter of the present invention, the step three formula (1) is obtained by the following formula (2) and formula (3):
Figure GDA0003107762540000062
Figure GDA0003107762540000063
wherein, formula (2) is obtained from formula (4) and formula (5):
Figure GDA0003107762540000064
Figure GDA0003107762540000065
where R is the length from a predetermined point (i.e., an arbitrary point) on the straight line Y1 (before rotation) to the XY-plane perpendicular line of the rectangular coordinate system, and L is the length from a point on the straight line Y1' corresponding to the position of the predetermined point on the straight line Y1 to the XY-plane perpendicular line of the rectangular coordinate system.
From equation (2), equation (3) can be derived from the same reasoning, as follows:
when the radius of the circle A is R, the rotation angle of the A, B plane around the Y axis is gamma, and the rotated straight lines Y1 and Y2 are the straight lines Y1 'and Y2', the included angles between the straight lines Y1 'and Y2' and the horizontal plane (XY plane) are respectively alpha 'and beta'.
The straight lines Y1 and Y1' are located on a cone with O as a circle point and the circle a as a bottom, and the geometrical relationship shows that:
Figure GDA0003107762540000071
Figure GDA0003107762540000072
the two formulas are combined to obtain:
Figure GDA0003107762540000073
the radius of the circle B of the bottom surface of the cone rotated by the straight line Y2 around the Y axis is R ', the length from any point on the back line of the Y2 rotation gamma to the vertical line of the horizontal plane is L', and the geometrical relationship shows that:
Figure GDA0003107762540000074
Figure GDA0003107762540000075
the two formulas are combined to obtain:
Figure GDA0003107762540000076
from the expressions (2) and (3), the pitch angle θ is, in any case:
Figure GDA0003107762540000077
in a specific example of the method for calculating the difference of the pitch angle of the helicopter of the present invention:
alpha and beta are respectively 30 degrees and 45 degrees, namely the pitch angle is 15 degrees. When the A1 and B1 planes rotate 20 degrees around the Y axis, namely gamma is 20 degrees, alpha 'and beta' are 28.02 degrees and 41.64 degrees, the difference between alpha 'and beta' is 16.12 degrees, and the value is substituted into the formula to obtain the value of theta of 15.00 degrees.
The differential calculation method for the pitch angle of the helicopter can realize the measurement of the pitch angle under the condition that the main rotor shaft is not vertically drilled, has stable precision and high efficiency compared with the traditional method, can also ensure the personal safety of measuring personnel, and reduces the labor intensity.
The above description is only for the specific embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.

Claims (2)

1. A method for calculating the difference of the pitch angle of a helicopter is characterized by comprising the following steps:
the method comprises the following steps of firstly, respectively arranging a first angle sensor and a second angle sensor on a blade mounting surface (A1) and a centerpiece plane (B1), wherein the X-axis direction of the first angle sensor and the X-axis direction of the second angle sensor are parallel to the axis of a blade, and the Y-axis of the first angle sensor and the Y-axis of the second angle sensor are in the same direction;
establishing a rectangular coordinate system, and establishing and simplifying a model according to the structural form of the step one, wherein the axis direction of the blade is the X-axis direction of the rectangular coordinate system, a straight line Y1 represents the Y-axis of the first angle sensor, the included angle between the straight line Y1 and the horizontal plane is alpha, a straight line Y2 represents the Y-axis of the second angle sensor, and the included angle between the straight line Y2 and the horizontal plane is beta;
thirdly, simulating the condition that the main rotor shaft is not vertical by controlling the three-axis rotation directions of the blade mounting surface (A1) and the center piece plane (B1) in the rectangular coordinate system X, Y, Z, and calculating the pitch angle of the corresponding condition, wherein the conditions comprise the following conditions:
a pitch angle θ is an angle between a Y-axis of the first angle sensor and a Y-axis of the second angle sensor when the blade mounting surface (A1) and a centerpiece plane (B1) are rotated about the X-axis of the rectangular coordinate system;
a pitch angle θ is an angle between a Y-axis of the first angle sensor and a Y-axis of the second angle sensor when the blade mounting surface (A1) and a centerpiece plane (B1) are rotated about a Z-axis of the rectangular coordinate system;
when the blade mounting surface (a1) and the centerpiece plane (B1) rotate around the Y axis of the rectangular coordinate system by an angle γ, the straight line Y1 'is the Y axis of the first rotated angle sensor, the included angle between the straight line Y1' and the horizontal plane is α ', the straight line Y2' is the Y axis of the second rotated angle sensor, the included angle between the straight line Y2 'and the horizontal plane is β', and the pitch angle θ is calculated according to the following formula (1):
Figure FDA0003107762530000021
2. a method of differential calculation of a helicopter pitch angle according to claim 1 wherein in said third step, equation (1) is derived from equations (2), 3) as follows:
Figure FDA0003107762530000022
Figure FDA0003107762530000023
wherein, formula (2) is obtained from formula (4) and formula (5):
Figure FDA0003107762530000024
Figure FDA0003107762530000025
wherein R is the length from a predetermined point on the straight line Y1 to the perpendicular line of the XY surface of the rectangular coordinate system, and L is the length from a point on the straight line Y1' corresponding to the position of the predetermined point on the straight line Y1 to the perpendicular line of the XY surface of the rectangular coordinate system;
from equation (2), equation (3) can be derived from the same reasoning, as follows:
the radius of the circle A is R, the rotation angle of the A, B surface around the Y axis is gamma, the rotated straight lines Y1 and Y2 are straight lines Y1 'and Y2', and the included angles between the straight lines Y1 'and Y2' and a horizontal plane (XY plane) are respectively alpha 'and beta';
the straight lines Y1 and Y1' are located on a cone with O as a circle point and the circle a as a bottom, and the geometrical relationship shows that:
Figure FDA0003107762530000026
Figure FDA0003107762530000031
the two formulas are combined to obtain:
Figure FDA0003107762530000032
the radius of the circle B of the bottom surface of the cone rotated by the straight line Y2 around the Y axis is R ', the length from any point on the back line of the Y2 rotation gamma to the vertical line of the horizontal plane is L', and the geometrical relationship shows that:
Figure FDA0003107762530000033
Figure FDA0003107762530000034
the two formulas are combined to obtain:
Figure FDA0003107762530000035
from the expressions (2) and (3), the pitch angle θ is, in any case:
Figure DEST_PATH_BDA0001710490890000077
CN201810678996.XA 2018-06-27 2018-06-27 Differential calculation method for pitch angle of helicopter Active CN108928476B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810678996.XA CN108928476B (en) 2018-06-27 2018-06-27 Differential calculation method for pitch angle of helicopter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810678996.XA CN108928476B (en) 2018-06-27 2018-06-27 Differential calculation method for pitch angle of helicopter

Publications (2)

Publication Number Publication Date
CN108928476A CN108928476A (en) 2018-12-04
CN108928476B true CN108928476B (en) 2021-09-21

Family

ID=64446379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810678996.XA Active CN108928476B (en) 2018-06-27 2018-06-27 Differential calculation method for pitch angle of helicopter

Country Status (1)

Country Link
CN (1) CN108928476B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113525712A (en) * 2021-06-23 2021-10-22 中国航空工业集团公司上海航空测控技术研究所 Helicopter rotor balance real-time monitoring and adjusting device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101376433A (en) * 2008-10-10 2009-03-04 南京航空航天大学 Helicopter rotor operation method and system
CN102678451A (en) * 2011-06-15 2012-09-19 华锐风电科技(集团)股份有限公司 Variable pitch propeller control method for wind generating set
CN104908976A (en) * 2015-05-19 2015-09-16 北京航空航天大学 Simple rotor mechanism of coaxial dual-rotor helicopter test stand
CN106786784A (en) * 2017-03-22 2017-05-31 广东工业大学 A kind of wind-powered electricity generation group of planes Poewr control method and system
CN106989812A (en) * 2017-05-03 2017-07-28 湖南科技大学 Large fan blade modal method of testing based on photogrammetric technology

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101376433A (en) * 2008-10-10 2009-03-04 南京航空航天大学 Helicopter rotor operation method and system
CN102678451A (en) * 2011-06-15 2012-09-19 华锐风电科技(集团)股份有限公司 Variable pitch propeller control method for wind generating set
CN104908976A (en) * 2015-05-19 2015-09-16 北京航空航天大学 Simple rotor mechanism of coaxial dual-rotor helicopter test stand
CN106786784A (en) * 2017-03-22 2017-05-31 广东工业大学 A kind of wind-powered electricity generation group of planes Poewr control method and system
CN106989812A (en) * 2017-05-03 2017-07-28 湖南科技大学 Large fan blade modal method of testing based on photogrammetric technology

Also Published As

Publication number Publication date
CN108928476A (en) 2018-12-04

Similar Documents

Publication Publication Date Title
CN103471619B (en) A kind of laser strapdown inertial navigation system prism ridge orientation installation error calibration
CN104792323B (en) Dynamic Water level and attitude angle computational methods
CN103712598B (en) Attitude determination method of small unmanned aerial vehicle
CN108204879B (en) A kind of measuring method and system of rotary inertia
CN107389266A (en) A kind of aircraft center of gravity flexible measuring method
US10655983B2 (en) Method for calibrating aircraft tri-axial balance and device thereof
CN105606129A (en) Measurement and calibration method for assisting in mounting of airplane inertial navigation finished product assembly
CN106403993A (en) Measuring method for installation errors of alignment prism
CN112129322B (en) Method for detecting and correcting installation error of strapdown inertial measurement unit and three-axis rotary table
CN109827593A (en) A kind of error self-calibrating method, system and storage medium based on more IMU
CN105758364A (en) Method for establishing collimation axis dynamic model of non-orthogonal axes laser theodolite
CN102607591B (en) Track data generation method for testing strap-down inertial navigation software
CN103925854A (en) Method for horizontally measuring missiles on basis of three-dimensional coordinate transformation principle
CN108928476B (en) Differential calculation method for pitch angle of helicopter
CN106403990A (en) Calibration apparatus for consistency of optical axes
CN106703787B (en) A kind of the well track parameter calculation method and device of radially horizontal well
CN102748010B (en) Attitude measurement system and method and oil well well track measuring system and method
CN111046584A (en) Precise adjustment method of satellite-borne instrument and equipment
CN113184211A (en) Rocket-assisted unmanned aerial vehicle thrust line adjusting method based on centroid test and balancing
CN105627982A (en) Remote vehicle inclined aiming method
WO2018137465A1 (en) Antenna posture data acquisition device and acquisition method, and antenna device
CN111141285A (en) Aviation gravity measuring device
CN110345902B (en) Method for testing control surface corner in non-horizontal state of airplane
CN111966953A (en) Aircraft attitude angle calculation method based on coordinate transformation
CN206930293U (en) Two wheeler

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